Normalized Dynamic Characterization and Application of Multiple Heat Storage Materials Based on Standard Thermal Resistance北大核心CSCD

基于标准热阻和能量流法,推导出储热材料与换热流体的瞬态换热热阻,通过类比电路分析法,获得了储热-换热过程的瞬态热量流模型及动态响应时间常数。进一步引入节点温度,重新定义换热热阻,获得了储热与换热过程耦合的三阶电路瞬态热量流模型,求解得到了加热、储热和释热三类时间常数,可用于协同表征储热材料中储热与释热的快慢程度,从而实现了多类储热材料的归一化动态表征。通过数值模拟验证与应用对比分析,发现基于多时间常数的归一化动态模型用于表征储热材料的动态特性是可行的,可直接对不同换热、储热材料进行对比分析。案例分析发现与固体储热材料换热时,液态金属的动态换热能力优于熔融盐,而相比于水蒸气和CO₂,空气与陶瓷材料换热能更快达到稳态。     

Unsteady Slip Flow and Heat Transfer Analysis of Oldroyd-B Fluid Over the Stretching Wedge北大核心CSCD

研究了在速度滑移现象存在下,上随体Oldroyd-B流体绕加热的楔形体的非稳态流动。采用松弛-延迟热通量模型,模拟了传热过程和热延迟时间对传热的影响,通过考虑浮升力、热辐射和对流换热边界条件,进一步研究了流动及传热特性。利用同伦分析方法获得常微分方程组的近似解析解,发现滑移参数的增大可以促进流体的流动,以及流体的温度随热辐射参数增大而升高。此外还发现,温度场在热松弛时间和热延迟时间中出现相反的变化趋势。     

Determination of heat transfer properties of media with a single-needle probe

Needle probes with a line heater inside are often used in studying the heat transfer properties of loose rocks. The key problem of contact methods of measuring thermal properties of various media consists in finding thermal contact resistance at the probe/medium interface which must be taken into account in determining the thermal diffusivity of the medium. We describe a mathematical model of heating of a long needle probe in the medium under study, taking into account dimensions and thermal properties of the needle source and assuming that thermal contact between the source and the medium is not ideal. Based on the proposed model, we formulate and solve the inverse problem of finding the thermal diffusivity coefficient of the medium and the heat exchange coefficient at the probe/medium interface. The purpose of the article is to create methodology for determining thermal properties of various media in the field.     

热特性参数可变时抛物线型截面环肋传热的最优化研究

肋片传热的最优化的传统处理方法是当执行一定的传热任务时,具有最少的肋片材料消耗(投资).最少重量的冷却肋片是具有抛物线型的截面.本文应用不变嵌入原理研究了当热特性参数可变时抛物线型截面环肋传热的最优几何尺寸,同时还对两个主要的物理参数,即导热系数变化参数α和放热系数变化指数m对最优几何尺寸的影响进行研究.所得结果对工程设计具有现实指导意义.     

Zonal modeling of radiative heat transfer in industrial furnaces using simplified model for exchange area calculation

The zonal analysis of industrial furnaces is considered with three-dimensional radiative heat transfer, incorporated with the mathematical zone method. In this method exchange areas are determined by simplified numerical integration in three dimensions for surface-surface, surface-gas and gas–gas zones for absorbing and emitting media. By focusing on new strategies to overcome the drawbacks of evaluating direct exchange areas, it is shown that the zone method is an effective numerical method for modeling three-dimensional thermal performance of gas-filled enclosures. Also the developed method for evaluating of exchange area is presented and compared with other methods in both sides of CPU time and accuracy. The method can decrease about 70% in error of calculation of some exchange areas as compared with the other numerical methods.     

Numerical investigation of residual thermal stresses in welded joints of heterogeneous steels with account of technological features of multi-pass welding

The paper describes a two-dimensional mathematical model to evaluate stresses in welded joints formed in multi-pass welding of multi-layered steels. The model is based on a system of equations that includes the Lagrange's variational equation of the incremental theory of plasticity and the Biot's variational principle for heat transfer simulation. In the constitutive equations, the changes in the volume which occur as a result of phase transitions can be taken into account. Therefore, the prehistory and impact of thermal processing of materials on macroscopic properties of the medium can be considered.The variational-difference method is used to solve both the heat transfer equation for calculation of the non stationary temperature field and the quasi-static problem of thermoplasticity at each time-step. The two-dimensional problems were solved to estimate the residual thermal stresses (for the case of plane stress or plane strain) during cooling of welds and assessing their impact on strain localization in the heat-affected zone under tensile and compressive loading considering differences in mechanical properties of welded materials.It is shown that at initial stages of the plastic flow, the residual stresses significantly affect the processes of stress concentration and localization of strains in welded joints. To estimate the model parameters and to verify the modeling results, the available experimental data from scientific literature obtained on the basis of the Satoh test for different welding alloys was used.     

Mixed convection heat transfer by nanofluids in a cavity with two oscillating flexible fins: A fluid–structure interaction approach

This study investigated the effects of two oscillating fins on the heat transfer rate and flow characteristics of a nanofluid inside a square enclosure. Both fins were attached to the hot wall and both fins oscillated at the same frequencies and amplitudes. The finite element method implemented in the arbitrary Lagrangian–Eulerian (ALE) technique was used to solve the equations describing the interactions and movements of the nanofluid and fins. Comparisons of our results and those reported in previous studies demonstrated that the modeling and numerical investigations were valid and reliable. The results showed that the increase in the heat transfer rate was due to the oscillation of the fins. In addition, the increasing trend in the heat transfer rate due to the oscillating fins decreased as the ratio of the thermal conductivity of the fins relative to the nanofluid increased. Increasing the thermal conductivity and viscosity parameters enhanced and weakened the heat transfer rate, respectively.     

主动冷却点阵夹层防热结构温度响应计算模型北大核心CSCD

针对点阵夹层结构主动热防护问题,建立了夹层结构面板和芯体导热与冷却剂对流耦合的非稳态传热理论模型,利用有限体积法离散控制方程并在MATLAB中进行了迭代求解.模型首次考虑了面板与夹芯杆之间的收缩热阻,并利用分离变量法得到了收缩热阻的近似解析解.基于单胞模型和周期性边界条件,模拟得到了模型所需的表面对流传热系数h_(b)和h_(fin).最后,选取多单胞计算工况进行数值模拟和理论模型对比,并讨论了收缩热阻对模型预测精度的影响.结果表明:理论模型能够准确预测夹层结构及内部流体的温度变化,理论与仿真之间的最大误差不超过1%;随着外加热流密度不断增大,忽略收缩热阻使得计算结果造成的误差不断增大;与数值模拟相比,理论模型可显著地减少计算时间并节省计算资源,尤其适用于非均匀、非稳态复杂热载荷下点阵夹层结构的温度响应计算.     

From SOI to Abstract Electric-Thermal-1D Multiscale Modeling for First Order Thermal Effects

Self-heating occurs in integrated circuits, specially for SOI-based devices. Naturally, the heat distribution affects the circuit’s functionality. For reliable designs in SOI-chip technology, and other applications, the thermal aspects have to be addressed. Therefore we develop a model, which is based on distributed 1D and lumped 0D elements, and takes into account that heat is stored and slowly conducted between elements. The emerging coupled multiscale system of heat evolution and electric network consists of parabolic partial-differential (thermal part) and differential-algebraic equations (electric network part). For the thermal model, we verify properties as positivity and strict passivity. Since time scales differ largely, the coupled problem exhibits multirate potential.     

Multi Input Dynamical Modeling of Heat Flow With Uncertain Diffusivity Parameter

This paper focuses on the multi-input dynamical modeling of one-dimensional heat conduction process with uncertainty on thermal diffusivity parameter. Singular value decomposition is used to extract the most significant modes. The results of the spatiotemporal decomposition have been used in cooperation with Galerkin projection to obtain the set of ordinary differential equations, the solution of which synthesizes the temporal variables. The spatial properties have been generalized through a series of test cases and a low order model has been obtained. Since the value of the thermal diffusivity parameter is not known perfectly, the obtained model contains uncertainty. The paper describes how the uncertainty is modeled and how the boundary conditions are separated from the remaining terms of the dynamical equations. The results have been compared with those obtained through analytic solution.     

Response of heat transfer from a moving flat plate in a parabolic flow

This paper deals with the solutions of steady as well as unsteady three-dimensional incompressible thermal boundary layer equations and the study of the response of heat transfer when there is a parabolic flow over a moving flat plate. The components of velocity in boundary layer are discussed by Sarma and Gupta and those results are used to analyse thermal boundary layer equations. A general analysis is made from which we deduce (i) Solutions of two-dimensional thermal boundary layer on a moving flat plate, (ii) Solutions of thermal boundary layer on a yawed flat plate, (iii) Solutions of thermal boundary layer when there is a parabolic flow over a moving flat plate by giving different values to β and Cx. Solutions are developed for large and small times and curves are drawn representing the variations of heat transfer from the plate with time for all the cases. The limiting time is also calculated.     

Modeling of thermal phenomena in single laser beam and laser-arc hybrid welding processes using projection method

This paper concerns mathematical and numerical modeling of thermal phenomena accompanying single laser and laser-arc hybrid butt welding of steel sheets. Coupled heat transfer and fluid flow in the fusion zone were described respectively by transient heat transfer equation and Navier–Stokes equation. Laser beam and electric arc heat sources were modeled using different heat source power distributions. Latent heat associated with the material’s state changes, buoyancy forces and liquid material flow through a porous medium were taken into account in considerations. Differential governing equations were numerically solved using projection method combined with finite volume method. Elaborated solution algorithm was implemented into computer solver used for simulation of heat transfer and fluid flow during welding. The geometry of the weld and heat affected zone as well as cooling rates were estimated on the basis of numerically obtained temperature field.     

Flow and heat transfer for a power-law electrically conducting fluid flowing between parallel plates under transverse magnetic field with viscous dissipation

The flow and heat transfer problem with viscous dissipation for electrically conducting non-Newtonian fluids with power-law model in the thermal entrance region of two parallel plates with magnetic field under constant heat flux and constant wall temperature conditions has been studied. The governing equations have been solved numerically using quasilinearization technique and implicit finite-difference scheme. It has been found that the effect of viscous dissipation on heat transfer is quite significant for heating and cooling conditions at the wall.     

Multi-Fracture Reservoir Simulations for Long Time Physical Processes

Reservoirs with multi-fracture techniques are developed and frequently used for oil and gas industry. Recently, they are also used for deep geothermal reservoirs especially for Hot Dry Rock (HDR). The analysis of the reservoir is generally interested in long time physical properties (10–100 years), e.g. fluid flow, heat transport etc. Typical CFD simulations are limited in this context. Here we developed a fluid flow and heat transport modeling in a multi-fracture reservoir based on the so-called Mixed Dimensional Model (MDM), which describes the different characteristic flows and the heat transport in different dimensions. In the mathematical point of view, these models are discretized based on the Cellular Automaton (CA) method combined with other necessary numerical techniques. The different cases of fluid flow and heat transport in multi-fracture reservoirs have been simulated and shown physical results very reasonably with less computational time. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

一种新的变域传热问题的边界条件

回顾古典热传导方程建立的假设条件的基础上,分析了热层材料由于表面烧蚀而引起的传热区域内部的热漏现象.基于能量守恒原理,利用有限元分析法推导出变域热传导方程,并得到了热漏函数的表达式,提出了变域传热问题边界条件的改进形式.为了检验这种边界条件的合理性,利用Crank-Nicholson法对此数学模型进行空间和时间离散化,并进行了数值仿真求解.仿真结果证实,基于边界条件改进形式的数学模型使计算更方便,结果更符合实际,从而为工程应用提供理论分析的依据.同时,该数学模型也为研究动边界发汗冷却控制问题奠定了理论研究的基础.     

热防护服-空气-皮肤热传导模型及其解析解

建立了高温环境下热防护服-空气-皮肤的热传导模型.利用热传导时,层合界面间温度相等和热流量连续的条件,结合微分思想,用分离变量法推导了微小时间段内模型热传导的解析解,然后通过循环得到整个时域内的解析解.利用求得的解析解分析了在80 ℃的环境温度下模型各位置温度和热流密度的变化情况,以及在不同环境温度下皮肤表面温度变化和热损伤情况.该求解方法可用来分析一般层合结构传热问题,计算结果对热防护服的设计和效果评价具有一定的参考意义.     

Chattering in dynamic mathematical two-phase flow models

This paper presents a chattering problem which arises in a dynamic mathematical two-phase flow model. The real system under study is also introduced, the DISS test facility, a parabolic-trough solar thermal power plant. The heat transfer fluid in the DISS facility is the steam-water mixture. A dynamic model of this plant, using Modelica as the modeling language, was previously developed in order to study its behavior. Chattering arises in the pipe model reducing the computational performance and hence limiting the applicability of the model. The problem source is studied and analysed together with an approach to the problem which is based on the smooth interpolation of some thermodynamic properties.     

非常量导热系数几例微分方程的推立

在本文里,曾先后假设物体的导热系数是依直线和指数函数空间地起改变,就这样来建立了六个二阶热传导微分方程;又对于变密度、变比热、变导热系数这样的更一般的情况也推立了六个二阶热传导的微分方程.     

Mathematical modelling of spark-ignition engines

The flow field, scavenging efficiency, power output, heat transfer losses, and unburned hydrocarbon emissions have been numerically studied by means of a two-equation model of turbulence in a four-stroke, homogeneous-charge, spark-ignition engine. The engine is equipped with an intake valve, an exhaust valve, and a constant rate heat source which simulates the spark plug. Combustion has been modelled by means of a one-step irreversible chemical reaction whose rate is controlled by an Arrhenius-type expression. The numerical results indicate that the intake stroke is characterized by the formation of two eddies which persist in the compression stroke. Turbulence is generated at the shear layers of the air jet drawn into the cylinder, but its level decreases in the compression stroke. Due to the heat released by the spark plug and the chemical reaction, a spherical flame kernel is formed. This kernel evolves into a cylindrical flame when the flame front reaches the piston. Fuel remains unburnt at the corner between the cylinder head and the cylinder wall due to heat transfer losses. The numerical results also indicate that despite uncertainties about the turbulence and heat transfer models, an engine model such as the one studied here can be used to understand the flow field, heat transfer losses, scavenging efficiency, and power output in conventional spark-ignition engines. Such capabilities are very helpful in the development and optimization stages of engines. For example, here the engine model thermal and scavenging efficiencies are 15.69% and 94%, respectively. The peak pressure is 33 atm and occurs at 6° ATDC. The unburnt hydrocarbon emissions are 7.41% of the total fuel admitted into the cylinder.