热泳和Brown运动对太阳能辐射下热分层纳米流体边界层流动的影响

在太阳辐射下的纳米流体中,数值地研究竖向延伸壁面具有可变流条件时的层流运动.使用的纳米流体模型为,在热分层中综合考虑了Brown运动和热泳的影响.应用一个特殊形式的Lie群变换,即缩放群变换,得到相应边值问题的对称群.对平移对称群得到一个精确解,对缩放对称群得到数值解.数值解依赖于Lewis数、Brown运动参数、热分层参数和热泳参数.得到结论:上述参数明显地影响着流场、温度和纳米粒子体积率的分布.显示出纳米流体提高了基流体热传导率和对流的热交换性能,基流体中的纳米粒子还具有改善液体辐射性能的作用,直接提高了太阳能集热器的吸热效率.     

温差影响下的局部滑移接触行为的研究

针对不同温度装配件间接触界面的局部滑移问题,建立了三维稳态热弹性局部滑移接触的半解析求解模型.基于热弹性理论与热传导方程,构建了半空间受热流载荷和力载荷作用下的频响函数并建立了相应的影响系数矩阵.借助离散卷积-快速Fourier变换等数学工具,实现了针对高温压头与热弹性半空间局部滑移接触问题的高效求解.接触界面间的热量传递满足Fourier热传导定律,并且黏/滑状态由Coulomb定律确定.基于该半解析模型分析了不同荷载及温差对表面法向压力分布、摩擦力分布、刚体位移及接触区黏/滑演化行为的影响.研究结果表明,当法向荷载和切向荷载一定时,温差的上升会导致接触区域的减小,引起接触面法向压力及摩擦力的峰值增大,并且会显著影响黏着区与滑移区的分布情况.     

非线性二维导热反问题的混沌-正则化混合解法

考虑热传导系数随温度变化,建立了非线性二维稳态导热反问题数值计算模型。并把混沌优化方法和梯度正则化方法相结合,构成一种混沌-正则化混合算法求该计算模型的全局解。以热传导系数随温度线性变化为例,由布置在结构边界上的观测点温度信息确定了结构材料热传导系数及其随温度变化规律。结果表明混合算法计算结果与初值无关,具有很好的全局寻优性能,而且计算量远比经典遗传算法和单纯采用混沌优化方法小。     

基于分数阶热传导方程激光加热瞬态温度场研究

基于分数阶Taylor(泰勒)级数展开原理,建立单相延迟一阶分数阶近似方程,获得分数阶热传导方程.针对短脉冲激光加热问题建立分数阶热传导方程组,并运用Laplace(拉普拉斯)变换方法进行求解,给出非Gauss(高斯)时间分布的激光内热源温度场解析解.针对具体算例数值研究温度波传播特性.结果表明热传播速度与分数阶阶次有关,分数阶阶次增加,热传播速度减小,温度变化幅度增加.分数阶方程可以用于描述介于扩散方程和热波方程间的热传输过程,且对热传播机制与分数阶热传导方程中分数阶项的关系做了深入剖析.     

考虑热-水-力耦合效应的饱和多孔地基动力响应分析

基于广义热弹性理论,结合Darcy(达西)定律,对Biot波动方程进行了修正,研究了一个受到椭圆余弦波作用的,均质各向同性半无限大饱和多孔地基的热-水-力多场耦合动态响应问题.建立了饱和多孔弹性地基的热-水-力耦合动力响应模型及控制方程,采用正则模态法求解,得到了问题的解析解,分析了地基中渗透系数变化和椭圆余弦波频率变化对饱和多孔地基中各物理量的影响.最终,给出了无量纲的竖向位移、超孔隙水压力、竖向应力和温度等物理量的分布规律.     

基于非傅里叶的有限空圆柱体的温度场解析解及其在谐波均匀的圆柱体上的应用

在有限空圆柱体的热冲击问题中引进非傅里叶分析,考虑了在极端热传递条件下的非稳态热传递过程中热量传播速度的有限性,建立了有限空圆柱体的轴对称非傅里叶温度场的数学模型,利用分离变量法和杜哈尔积分求得有限空圆柱体双曲型热传导问题的精确解析解.并将结果应用于外表面是谐波均匀的热通量的有限空圆柱体,得到其瞬态温度场及其径向、轴向温度分布规律.     

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

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

理想传导弹性体中能量耗散的磁-热-弹性波

广义能量耗散弹性理论(TEWED,G-NⅢ理论)广泛应用于均匀磁场作用下的时谐平面波在无限大的理想导电弹性体中传播的研究.提出了更普遍的有复杂参数的色散方程,通过运用Ieguerre 方法解决复杂条件下耦合磁-热-弹性波的问题,表明耦合磁-热-弹性波问题相当于改进的膨胀波及通过有限热波速度、热弹性耦合、热扩散率及外加磁场修正的、有限速度热波的传播问题.在G-NⅢ模型(TEWED)中,耦合磁-热-弹性波传播时发生衰减和色散,扩散的热量由热传播方程中的阻尼项考虑,而在G-N Ⅱ模型没有发生衰减和耗散.最后给出了类铜材料的数值结果.     

表面效应对热电材料中纳米孔周围热应力的影响

基于完整Gurtin-Murdoch(G-M)低阶表面能模型，进一步探讨了纳米尺度下表面效应的影响.建立了合理考虑构型变化的应力边界条件，实现了研究尺度从宏观到微观的转变.利用复变函数理论和保角映射技术，构建了用于纳米尺度下的热-电-力理论框架模型，得到了热电基体中纳米孔周围热场、温度场以及应力场的半解析解.数值结果表明，相对于完整G-M模型，简化G-M模型(忽略孔洞构型变化的影响)往往会高估表面效应和远场热电载荷对热应力分布的影响.此外，表面效应的存在将在一定程度上缓解纳米孔周围的热应力集中.     

分数阶广义扰动热波方程的与泛函映射解

研究了一类分数阶广义非线性扰动热波方程.首先在典型分数阶热波方程情形下得到解,接着用泛函分析映射方法,求出了分数阶广义非线性扰动热波方程初始边值问题的任意次近似解析解.最后简述了它的物理意义.求得的近似解析解,弥补了单纯用数值方法得到的模拟解的不足.     

Use of 3D-transient analytical solution based on Green’s function to reduce computational time in inverse heat conduction problems

Inverse problems can be found in many areas of science and engineering and can be applied in different ways. Two examples can be cited: thermal properties estimation or heat flux function estimation in some engineering thermal process. Different techniques for the solution of inverse heat conduction problem (IHCP) can be found in literature. However, any inverse or optimization technique has a basic and common characteristic: the need to solve the direct problem solution several times. This characteristic is the cause of the great computational time consumed. In heat conduction problem, the time consumed is, usually, due to the use of numerical solutions of multidimensional models with refined mesh. In this case, if analytical solutions are available the computational time can be reduced drastically. This study presents the development and application of a 3D-transient analytical solution based on Green’s function. The inverse problem is due to the thermal properties estimation of conductors. The method is based on experimental determination of thermal conductivity and diffusivity using partially heated surface method without heat flux transducer. Originally developed to use numerical solution, this technique can, using analytical solution, estimate thermal properties faster and with better accuracy.     

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

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

Simulation of thermal disturbances with finite wave speeds using a high order method

The classical heat diffusion theory based on the Fourier’s model breaks down when considering transient heat flow, for short times, extreme thermal gradients or at low temperatures. The hyperbolic heat conduction equation based on the Cattaneo model for the heat flux incorporates a relaxation mechanism in order to gradually adjust to a change in the temperature gradient. A spectral element method is applied for solving the hyperbolic system treating the heat flux as an independent variable in addition to temperature. The numerical solution is based on the time-space least squares spectral method. Numerical examples are included for discussing the effects of the thermal waves.     

Finite integral transform-based analytical solutions of dual phase lag bio-heat transfer equation

A finite integral transform (FIT)-based analytical solution to the dual phase lag (DPL) bio-heat transfer equation has been developed. One of the potential applications of this analytical approach is in the field of photo-thermal therapy, wherein the interest lies in determining the thermal response of laser-irradiated biological samples. In order to demonstrate the applicability of the generalized analytical solutions, three problems have been formulated: (1) time independent boundary conditions (constant surface temperature heating), (2) time dependent boundary conditions (medium subjected to sinusoidal surface heating), and (3) biological tissue phantoms subjected to short-pulse laser irradiation. In the context of the case study involving biological tissue phantoms, the FIT-based analytical solutions of Fourier, as well as non-Fourier, heat conduction equations have been coupled with a numerical solution of the transient form of the radiative transfer equation (RTE) to determine the resultant temperature distribution. Performance of the FIT-based approach has been assessed by comparing the results of the present study with those reported in the literature. A comparison of DPL-based analytical solutions with those obtained using the conventional Fourier and hyperbolic heat conduction models has been presented. The relative influence of relaxation times associated with the temperature gradients (τ_T) and heat flux (τ_q) on the resultant thermal profiles has also been discussed. To the best of the knowledge of the authors, the present study is the first successful attempt at developing complete FIT-based analytical solution(s) of non-Fourier heat conduction equation(s), which have subsequently been coupled with numerical solutions of the transient form of the RTE. The work finds its importance in a range of areas such as material processing, photo-thermal therapy, etc.     

Modelling non-Fourier heat conduction with periodic thermal oscillation using the finite integral transform

The non-Fourier heat conduction in a finite medium subjected to a periodic heat flux is modelled using the finite integral transform technique and an analytic solution is obtained. An analogy between thermal oscillation and oscillation of mechanical and electrical systems is drawn. A transition criterion from the non-Fourier heat conduction formulation to the Fourier formulation is obtained and a simple analytical expression of the phase and amplitude of thermal oscillation is derived.     

Transient thermoelastic response in a cracked strip of functionally graded materials via generalized fractional heat conduction

This work is devoted to analyzing a thermal shock problem of an elastic strip made of functionally graded materials containing a crack parallel to the free surface based on a generalized fractional heat conduction theory. The embedded crack is assumed to be insulated. The Fourier transform and the Laplace transform are employed to solve a mixed initial-boundary value problem associated with a time-fractional partial differential equation. Temperature and thermal stresses in the Laplace transform domain are evaluated by solving a system of singular integral equations. Numerical results of the thermoelastic fields in the time domain are given by applying a numerical inversion of the Laplace transform. The temperature jump between the upper and lower crack faces and the thermal stress intensity factors at the crack tips are illustrated graphically, and phase lags of heat flux, fractional orders, and gradient index play different roles in controlling heat transfer process. A comparison of the temperature jump and thermal stress intensity factors between the non-Fourier model and the classical Fourier model is made. Numerical results show that wave-like behavior and memory effects are two significant features of the fractional Cattaneo heat conduction, which does not occur for the classical Fourier heat conduction.     

On dynamic effects in an elastic half-space under “thermal impact” : PMM vol. 38, n 6, 1974, pp. 1105–1113

The general uncoupled dynamical problem of thermoelasticity for a half-space under the condition of a thermal impact with a finite rate of change in temperature on its boundary is solved by the method of principal (fundamental) functions within the framework of a generalized theory of heat conduction.An elastic steel half-space is analyzed as an illustration. The problem on thermal stresses originating in an elastic half-space due to thermal impact produced by a jump change in temperature on the boundary was first analyzed in [1]. Since the temperature change on the boundary occurs at a finite rate, it is generally impossible to realize the thermal impact considered in [1] physically. The dynamic effects in an elastic half-space under a thermal impact with finite rate of change in the temperature on the boundary have been studied in [2]. For high rates of change of the heat flux we obtain a generalized wave equation of heat conduction [3] taking into account the finite velocity of heat propagation. Hence, the solution of the ordinary parabolic heat conduction equation used in [1, 2] does not correspond to the true temperature field. The problems of [1, 2] have been examined in [4, 5], respectively, within the framework of a generalized theory of heat conduction.     

应力降低因子和陶瓷材料抗热冲击阻力参数的性质及适用条件

将对流条件下薄板的瞬态导热问题的解析解引入自由弹性薄板的热应力场模型中,给出了相应应力降低因子的具体表达形式．为了便于比较,进一步定义了一个新的应力降低因子．详细讨论了应力降低因子及分别对应于高Biot模数和低Biot模数的第1个和第2个抗热冲击阻力参数及与中间量级的Biot模数相对应的近似表达式的性质及适用条件．将传热学与热弹性力学或断裂力学相结合的方法及有限元方法是该文所推荐的抗热震性能计算方法．研究表明,采用断裂临界温差和断裂临界无量纲时间相结合的方式能够直观简洁地表征陶瓷材料的抗热震性能．