关于热波理论的研究

当热传导服从经典的Fourier定律时,温度场由抛物型方程所控制,热扰动以无限大速度传播。在通常情况下,热波的迅速衰减掩盖了这种佯谬。但对于热爆炸,热核聚变,快速化学反应,强激光与物质相互作用这样一些时间尺度极短(与从非平衡态达到局部平衡态的时间相比)的情形,需对Fourier定律进行修正。本文从Cattaneo方程及其唯象修正,Boltzmann演化方程,分子动力学数值模拟以及连续介质热力学理论四个方面对热波理论近20年的进展进行了评述,并对热波的实验验证及数值分析的某些重要结果进行了介绍,且提出了关于热波非弹性理论的新认识,讨论了可能的应用前景。     

激光脉冲加热的二维广义热弹性问题研究

基于L-S广义热弹性理论,研究了半无限大板局部受到激光脉冲加热时的广义热弹性问题.为避免常规积分变换方法求解带来的精度丢失,采用有限元法直接在时间域进行求解,得到了激光脉冲加热时板中的温度、位移及应力的变化规律.结果表明,直接求解方法可以准确描述热在介质中以有限的速度传播,同时发现,激光脉冲加热过后,结构的最高温度随着时间的推移逐渐降低,且最高温度的位置总在热波波前附近,此处的应力也明显高于其他区域.     

分数阶理论下实心球体受热冲击作用的渐进分析

基于分数阶广义热弹性理论,针对实心球体在外表面受均匀热冲击作用下的一维广义热弹性问题进行研究分析. 利用热冲击的瞬时特征,借助于Laplace 正、反变换技术及柱函数的渐进性质,推导了热冲击作用周期内位移场、温度场和应力场的渐进表达式. 通过计算,得到了不同传热能力下受热冲击作用时热波、热弹性的传播规律以及位移场、温度场及应力场的分布规律. 结果表明:分数阶参数取值的不同,热波、热弹性波的传播以及各物理场的分布均有所不同,分数阶参数可视为延迟时间的影响因子,通过改变延迟效应对热弹性行为的影响来改变热冲击的作用效果.      

分数阶广义热弹性理论下中空柱热弹性分析

基于Laplace变换及特征值法,推导并给出了分数阶广义热弹性理论下中空柱内表面作用有热冲击情况的解析解,通过Laplace数值逆变换法求解得到了位移场、温度场、应力场的分布规律。结果表明:特征值法能准确给出Laplace域内方程组的解;分数阶参数对温度场和应力场有较大影响,对位移场影响较小。作为广义热弹性理论的一种推广,在处理热传导问题时,通过分数阶广义热弹性理论进行研究更科学、全面。     

非稳态分数阶Oldroyd-B流体绕楔形体流动研究

研究了非稳态分数阶Oldroyd-B流体在多孔介质中通过楔形拉伸板的驻点流动问题。基于分数阶Oldroyd-B流体的本构模型建立了动量方程,并在其中引入了浮升力和驻点流动特征。此外,考虑了具有热松弛延迟时间的修正的分数阶Fourier定律,并将其应用于能量方程和对流换热边界条件。接着,采用与L1算法相结合的有限差分法求解控制偏微分方程。最后,分析了相关物理参数对流动的影响。结果表明,随着楔角参数的增加,流体受到的浮升力增大,导致速度加快;达西数越大,介质的孔隙度变大,流体的流动越快;此外,温度分布先略有上升后明显下降,这表明Oldroyd-B流体具有热延迟特性。     

电磁加热条件下皮肤组织的生物热力学行为

研究目的是开发一种数学方法来计算传热过程、热引起的力学响应以及相应的疼痛等级, 从而对临床上应用的各种加热疗法之间的差别进行定量评估. 采用基于有限差分法的数值模拟方法, 基于无限大和均匀化假设, 分析了各种热疗法中皮肤组织的温度、烧伤和热应力分布. 研究发现: 充血对热损伤的影响很小, 但对皮肤的温度分布影响很大, 而这又反过来显著影响由此产生的热应力场; 对于激光加热, 光波越短则峰值温度越高, 但峰值更接近皮肤表面温度; 激光和微波加热所产生的热应力集中于表皮顶层, 因为发热量沿皮肤深度方向呈指数衰减; 薄角质层(常常被忽略)对皮肤组织的热力学响应起主导作用.      

GCr15凹模环形裂纹电磁热局部跨越止裂的模拟研究

裂纹止裂技术是当前工业生产和实际应用中急需解决的关键问题;利用电磁热效应对带有裂纹的金属构件进行止裂,是达到延长工作寿命、提高安全性可靠性的一种行之有效的方法;本文采用有限元方法研究了带有半埋藏环形裂纹的轴对称GCr15凹模电磁热局部跨越止裂;建立了超强脉冲电流放电瞬间,环形裂纹尖端附近电场、温度场和热应力场的耦合关系。数值模拟分成热-电耦合(焦耳热问题)和热-机械耦合两个过程来实现,通过数值解析轴对称模型,给出了脉冲放电瞬间GCr15凹模中的电流-温度-热应力的求解过程,得到了脉冲放电瞬间半埋藏环形裂纹局部跨越止裂时凹模的温度场和热应力场。数值模拟研究结果表明:由于热集中效应,在放电瞬间裂纹尖端金属熔化,并形成局部堆焊,阻止了干线裂纹源的开裂,达到了止裂的目的。     

超短激光脉冲加热薄板的广义热弹扩散问题

由于超短激光脉冲具有功率密度高、持续时间短、加工精度高等优势, 近年来被广泛应用于超精细加工、光学储存和微电子器件制造等领域. 本文基于L-S型广义热弹扩散理论, 建立了考虑材料记忆依赖效应和空间非局部效应的记忆依赖型非局部广义热弹扩散耦合理论, 它能够准确预测几何尺寸与内部特征尺寸相近结构的热弹扩散瞬态响应. 推导了所建理论的控制方程, 并基于拉普拉斯积分变换获得了控制方程的解. 作为算例, 利用所建理论和求解方法研究了半无限大薄板受非高斯激光脉冲加热和化学冲击联合作用下的热弹扩散瞬态响应问题, 得到了薄板的温度、化学势、位移、应力和浓度等随非局部参数、热时间迟滞因子和扩散时间迟滞因子等参数变化的分布规律. 结果表明: 传热对传质影响显著, 传质对传热影响甚微; 非局部参数对位移、应力影响显著, 对温度、化学势和浓度几乎没有影响. 该理论及求解方法的建立, 旨在实现材料在机械、热、化学势等冲击作用下传热传质瞬态响应的准确预测.      

A GENERALIZED THERMOELASTIC DIFFUSION PROBLEM OF THIN PLATE HEATED BY THE ULTRASHORT LASER PULSES 1)

由于超短激光脉冲具有功率密度高、持续时间短、加工精度高等优势, 近年来被广泛应用于超精细加工、光学储存和微电子器件制造等领域. 本文基于L-S型广义热弹扩散理论, 建立了考虑材料记忆依赖效应和空间非局部效应的记忆依赖型非局部广义热弹扩散耦合理论, 它能够准确预测几何尺寸与内部特征尺寸相近结构的热弹扩散瞬态响应. 推导了所建理论的控制方程, 并基于拉普拉斯积分变换获得了控制方程的解. 作为算例, 利用所建理论和求解方法研究了半无限大薄板受非高斯激光脉冲加热和化学冲击联合作用下的热弹扩散瞬态响应问题, 得到了薄板的温度、化学势、位移、应力和浓度等随非局部参数、热时间迟滞因子和扩散时间迟滞因子等参数变化的分布规律. 结果表明: 传热对传质影响显著, 传质对传热影响甚微; 非局部参数对位移、应力影响显著, 对温度、化学势和浓度几乎没有影响. 该理论及求解方法的建立, 旨在实现材料在机械、热、化学势等冲击作用下传热传质瞬态响应的准确预测.     

轴对称横观各向同性热弹性圆柱的分解

为了得到精确的应力场、位移场、温度场,将扭转圆轴的精化理论研究方法推广到轴对称横观各向同性热弹性圆柱。利用Bessel函数以及轴对称横观各向同性热弹性圆柱的通解,给出了轴对称横观各向同性热弹性圆柱的分解定理。根据柱面齐次边界条件获得了精确的精化方程,精化方程可以分解为一阶方程、超越方程、温度方程,从而将横观各向同性热弹性圆柱的轴对称问题分解为轴向拉压问题、超越问题、热-应力耦合问题。超越部分对应端部自平衡情况,可以清晰地了解到端部应力分布对内部应力场的影响,热-应力耦合部分对应无外加应力场时圆柱内部因温度变化引起的热应力。     

Analytical solution of time fractional Cattaneo heat equation for finite slab under pulse heat flux

A fractional Cattaneo model is derived for studying the heat transfer in a finite slab irradiated by a short pulse laser. The analytical solutions for the fractional Cattaneo model, the classical Cattaneo-Vernotte model, and the Fourier model are obtained with finite Fourier and Laplace transforms. The effects of the fractional order parameter and the relaxation time on the temperature fields in the finite slab are investigated. The results show that the larger the fractional order parameter, the slower the thermal wave. Moreover, the higher the relaxation time, the slower the heat flux propagates. By comparing the fractional order Cattaneo model with the classical Cattaneo-Vernotte and Fourier models, it can be found that the heat flux predicted using the fractional Cattaneo model always transports from the high temperature to the low one, which is in accord with the second law of thermodynamics. However, the classical Cattaneo-Vernotte model shows that the unphysical heat flux sometimes transports from the low temperature to the high one.     

结构层热扩散系数对超薄涂层热力性能的影响

将非傅立叶热传导模型(用于超薄热涂层)与傅立叶热传导模型(用于结构层)相结合 求解温度场,运用有限元法求解热涂层热应力和裂纹驱动力,并分析结构层材料热扩散系数 的变化对热涂层的热力学性能(温度场、应力场和断裂性能)的影响. 研究表明,结构层材 料性能变化对温度场的影响主要表现在热冲击后期,对热应力和裂纹尖端驱动力后期的变化 也有一定的影响.     

THERMAL SHOCK OF SEMI-INFINITE BODY WITH MULTI-PULSED INTENSE LASER RADIATION

In order to investigate the thermal shock and the heat conduction property of a target under multi-pulsed laser radiation, analytic expressions of both temperature and thermal stress fields in the target are deduced on the basis of the non-Fourier conduction law and the thermo-elastic theory. Taking a stainless steel target as an example, we can solve the analytic expressions under appropriate boundary conditions by using the finite difference method and MATLAB software, and then reveal the evolution law of both surplus temperaturt, and thermal stress in the target. The results indicate that the temperature curves in the target irradiated by a multi-pulsed laser take on a delayed character in different sections away from the boundary, which is only affected by its relaxation time. The front of the stress wave is very steep in the non-Fourier numerical solutions, which presents an obvious thermal shock, so it is necessary to consider the non-Fourier effect of semi-infinite body under the high energy laser radiation.     

Pulsed laser heating – Fourier and electron kinetic theory approaches

Laser surface pulse heating of engineering metals is in demand in the metal industry and investigation into laser pulse heating becomes fruitful in this regard. Application of Fourier theory to heat conduction due to high power laser irradiation may give closed form solution to the problem. On the other hand, the heat flux through a given plane depends on the electron energy distribution through the material and at the scale of distance required to examine the problem, the material can no longer be considered as being homogeneous continuum, therefore, errors may occur when considering the Fourier theory in laser heating process. The problem requires to be examined in the quantum field. The present study examines the pulse laser heating process when considering both Fourier conduction and electron-kinetic theory approaches. Analytical solution to Fourier conduction equation is obtained for intensity exponential pulses while numerical scheme is introduced to solve the heat transfer equation resulted from kinetic theory approach. It is found that both Fourier and electron kinetic theory approaches result in similar temperature profiles for the pulses having the same energy content. In the case of electron kinetic theory approach the rise time for surface temperature to reach the melting point is shorter than that obtained from the analytical solution. Received on 23 February 1998     

The nonlinear response of Cattaneo-type thermal loading of a laser pulse on a medium using the generalized thermoelastic model

The nonlinear thermoelastic responses of an elastic medium exposed to laser generated shortpulse heating are investigated in this article. The thermal wave propagation of generalized thermoelastic medium under the impact of thermal loading with energy dissipation is the focus of this research. To model the thermal boundary condition(in the form of thermal conduction),generalized Cattaneo model(GCM) is employed. In the reference configuration, a nonlinear coupled Lord-Shulman-type generalized thermoelasticity formulation using finite strain theory(FST) is developed and the temperature dependency of the thermal conductivity is considered to derive the equations. In order to solve the time-dependent and nonlinear equations, Newmark's numerical time integration technique and an updated finite element algorithm is applied and to ensure achieving accurate continuity of the results, the Hermitian elements are used instead of Lagrangian's. The numerical responses for different factors such as input heat flux and nonlinear terms are expressed graphically and their impacts on the system's reaction are discussed in detail.The results of the study are presented for Green–Lindsay model and the findings are compared with Lord-Shulman model especially with regards to heat wave propagation. It is shown that the nature of the laser's thermal shock and its geometry are particularly determinative in the final stage of deformation. The research also concluded that employing FST leads to achieving more accuracy in terms of elastic deformations; however, the thermally nonlinear analysis does not change the results markedly. For this reason, the nonlinear theory of deformation is required in laser related reviews, while it is reasonable to ignore the temperature changes compared to the reference temperature in deriving governing equations.     

热环境对超高温陶瓷材料抗热震性能的影响

基于超高温陶瓷材料相关材料参数对温度的敏感性,以现有的抗热震断裂评价理论为基础,通过考虑温度对超高温陶瓷材料物理性能的影响,探讨了热震温度范围内热应力衰减系数、表面热传递系数、热冲击初始环境温度等热环境参数对抗热冲击阻力参数、断裂临界温差和断裂临界温度的影响.      

两微极磁热弹性固体分界面上波的反射与折射

给出了磁场、热场和弹性场多场耦合作用下微极广义热弹性固体的一般控制方程.该方 程既包含了磁场、热场和弹性场的耦合作用,又在其广义热传导方程中涵盖了耦合热弹理论 (C-D)及其5类推广(L-S理论,G-L理论,G-N(II,III)理论和C-T理论).运用该微极广义磁热 弹性控制方程,研究了在定常磁场作用下, 具有均匀初始温度的两理想接触微极弹性介质平面分界面上磁热弹性波的反射和折射现象.给出了分别在缺少磁场、热场作用或不同广义热传 导理论下反射或折射热波、纵向位移波、耦合横向和微旋转波与入射纵向位移波的振幅比随 入射角变化的关系曲线.对缺少磁、热和微极性以及热松弛时间时对应的反射、折射系数进 行了对比.结果表明磁、热和微极性以及热松弛时间对振幅比均有不同程度的影 响,与磁、热和微极性一样,热松弛时间对不同类型波的影响能力差别明显,但对同 一类型的反射波和折射波的影响相似.     

Linear and non-linear double diffusive convection in a rotating porous layer using a thermal non-equilibrium model

Double diffusive convection in a fluid-saturated rotating porous layer heated from below and cooled from above is studied when the fluid and solid phases are not in local thermal equilibrium, using both linear and non-linear stability analyses. The Darcy model that includes the time derivative and Coriolis terms is employed as momentum equation. A two-field model that represents the fluid and solid phase temperature fields separately is used for energy equation. The onset criterion for stationary, oscillatory and finite amplitude convection is derived analytically. It is found that small inter-phase heat transfer coefficient has significant effect on the stability of the system. There is a competition between the processes of thermal and solute diffusions that causes the convection to set in through either oscillatory or finite amplitude mode rather than stationary. The effect of solute Rayleigh number, porosity modified conductivity ratio, Lewis number, diffusivity ratio, Vadasz number and Taylor number on the stability of the system is investigated. The non-linear theory based on the truncated representation of Fourier series method predicts the occurrence of subcritical instability in the form of finite amplitude motions. The effect of thermal non-equilibrium on heat and mass transfer is also brought out.