Analysis of non-Fourier heat conduction in a solid sphere under arbitrary surface temperature change

In this paper, the non-Fourier heat conduction in a solid sphere under arbitrary surface thermal disturbances is solved analytically. Four cases including sudden, simple harmonic periodic, triangular and pulse surface temperature changes are investigated step-by-step. The analytical solutions are obtained using the separation of variables method and Duhamel’s principle along with the Fourier series representation of an arbitrary periodic function and the Fourier integral representation of an arbitrary non-periodic function. Using these analytical solutions, the temperature profiles of the solid sphere are analyzed, and the differences in the temperature response between the “hyperbolic” and “parabolic” are discussed. These solutions can be applicable to all kinds of non-Fourier heat conduction analyses for arbitrary boundary conditions occurred in technology. And as application examples, particular attention is devoted to the cases of triangular surface temperature change and pulse surface temperature change. The examples presented in this paper can be used as benchmark problems for future numerical method validations.     

Analytical solution of non-fourier temperature response in a finite medium under laser-pulse heating

A temperature wave solution predicted by the hyperbolic heat equation is developed for a finite medium exposed to a surface heat flux of laser with an actual temporal profile. By using the analytical solution, the temperature response, the propagation and the reflection of the temperature wave due to such heat pulse are investigated for different pulse duration, thickness of the medium, and energy absorption depth.     

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     

Periodic heat conduction with relaxation time in cylindrical geometry

Steady-periodic heat conduction with relaxation time in an infinitely long hollow cylinder is considered. Four boundary value problems, with boundary conditions of the first and of the second kind, are solved analytically. The solution for a solid cylinder with a sinusoidally varying surface temperature is obtained as a special case of a solution found for the hollow cylinder. The effects of the relaxation time on the steady-periodic temperature field are analysed, in details, for a solid cylinder with a sinusoidally varying surface temperature and for a hollow cylinder with a sinusoidally varying heat flux at the inner surface and with a constant temperature at the outer surface. The results show that thermal resonances may occur and suggest that accurate measurements of the relaxation time could be obtained by means of experiments on steady-periodic heat conduction in cylindrical geometry. Received on 15 April 1997     

Axisymmetric non-Fourier temperature field in a hollow sphere

The non-Fourier axisymmetric (2+1)-dimensional temperature field within a hollow sphere is analytically investigated by the solution of the well-known Cattaneo–Vernotte hyperbolic heat conduction equation. The material is assumed to be homogeneous and isotropic with temperature-independent thermal properties. The method of solution is the standard separation of variables method. General linear time-independent boundary conditions are considered. Ultimately, the presented solution is applied to a (1+1)—as well as a (2+1)—dimensional problem, and their respective non-Fourier thermal behavior is studied. The present solution can be reduced to special cases of interest by choosing appropriate boundary conditions parameters. Dedicated to Prof. Gholamali Atefi, with appreciation and admiration on the occasion of his 65th birthday.     

磁-电-弹性半空间在轴对称热载荷作用下的三维问题研究

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.      

THREE-DIMENSIONAL ANALYSIS OF A MAGNETOELECTROELASTIC HALF-SPACE UNDER AXISYMMETRIC TEMPERATURE LOADING 1)

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.     

Influence of laser pulse shape both on temperature profile and hardened layer depth

This paper presents a mathematical modelling and numerical calculations of heat conduction problems where laser generated heat is assumed as a surface heat source. Also the effect of a laser time structure on a hardened layer depth is examined. Temperature profiles for different laser pulse shapes are determined from the solution of a linear one-dimensional heat conduction equation for semi-infinite medium and discussed in terms of the parameters evolution such as dimensionless: temperature, heat flux, hardening depth, laser impulse duration and increasing time of triangular pulse shape.     

Solution and analysis of hyperbolic heat propagation in hollow spherical objects

The hyperbolic heat conduction process in a hollow sphere with its two boundary surfaces subject to sudden temperature changes is solved analytically by means of integration transformation. An algebraic analytical expression of the temperature profile is obtained. Accordingly, the non-Fourier hyperbolic heat propagation in hollow spherical medium is analyzed and possible hyperbolic anomalies are discussed.     

Heating process simulation for a body subjected to an intensive thermal effect on its surface

An approximate method of calculating the nonstationary temperature fields for the case of rapid and intensive heating of bodies is considered. The fast heating of a body is simulated by various boundary conditions at its heated surface. An approximate method is proposed in a one-dimensional formulation to solve the unsteady-state heat conduction equation with these boundary conditions. The applicability of this method in practice is confirmed by comparing the approximate temperature fields with the exact solution. An approximate temperature field is constructed for a particular case with consideration of a temperature dependence of thermal parameters, which allows one to describe the heat localization in a heated zone.     

非傅立叶导热的最新研究进展

对迄今为止有关非傅立叶导热的研究成果进行了全面的综述,其中包括作者在该领域的最新研究进展:空心球体介质双曲线非傅立叶导热模型的分析求解,室温条件下多孔材料内非傅立叶导热的实验结果及数值模拟,非傅立叶导热的“瞬时薄层”模型,非傅立叶导热和非费克质量传递的耦合分析,非傅立叶导热的分子动力学模拟等．文中还对下一步的研究工作进行了展望.      

On the non-classical problem of thermal conduction on anisotropic layers in R3

We propose a solution methodology for boundary problems of parabolic and hyperbolic thermal conduction on anisotropic layers in R³. We study the wave nature of heat transfer with pulse thermal effects in bounded bodies with cavity. We compare the solutions of the parabolic and hyperbolic equations of thermal conduction and we show that the assumption on the wave nature of energy transfer is justified under the conditions for high-speed processes.     

Transient thermal cracking associated with non-classical heat conduction in cylindrical coordinate system

This paper studies the fracture behavior of a thermoelastic cylinder subjected to a sudden temperature change on its outer surface within the framework of non-classical heat conduction.The heat conduction equation is solved by separation of variable technique.Closed form solution for the temperature field and the associated thermal stress are established.The critical parameter governing the level of the transient thermal stress is identified.Exact expression for the transient stress intensity factor is obtained for a crack in the cylinder.The difference between the non-classical solutions and the classical solution are discussed.It is found that the traditional classical heat conduction considerably underestimates the transient thermal stress and thermal stress intensity factor.     

Non-fourier heat conduction in a plane slab with prescribed boundary heat flux

The non-stationary heat conduction in an infinitely wide plane slab with a prescribed boundary heat flux is studied. An arbitrary time dependent boundary heat flux is considered and a non-vanishing thermal relaxation time is assumed. The temperature and the heat flux density distributions are determined analytically by employing Cattaneo-Vernotte's constitutive equation for the heat flux density. It is proved that the temperature and the heat flux density distributions can be incompatible with the hypothesis of local thermodynamic equilibrium. A comparison with the solution which would be obtained by means of Fourier's law is performed by considering the limit of a vanishing thermal relaxation time.     

短脉冲激光加热分数阶导热及其热应力研究

短脉冲激光加热引起材料内部复杂的传热过程及热变形,现有的以Fourier定律或Cattaneo-Vernotte松弛方程结合弹性理论为框架建立起来热应力理论在刻画其热物理过程存在严重缺陷. 本文基于分数阶微积分理论, 以半空间为研究对象, 建立了分数阶Cattaneo热传导方程和相应的热应力方程, 给出了问题的初始条件和边界条件, 采用拉普拉斯变换方法, 给出了非高斯时间分布激光热源辐射下温度场和热应力场的解析解, 研究了短脉冲激光加热的温度场及热应力场的热物理行为. 数值计算中, 首先对理论解进行数值验证, 然后取分数阶变量$p=0.5$研究温度场和热应力场的变化特点及激光参数对温度和热应力的影响,最后数值计算分数阶参数对温度和热应力场的影响. 计算结果表明, 分数阶Cattaneo传热方程和热应力方程描述的温度和热应力任然具有波动特性,与经典的Fourier传热模型和标准的Cattaneo传热模型相比, 分数阶阶次越大, 热波波速越小, 热波波动性越明显; 反之, 则热波波速越大, 热扩散性越强.激光加热和冷却的速度越快, 温度上升和下降的速度越快, 压应力和拉应力交替变化越快, 温度变化幅值越小, 热应力幅值影响不明显.      

The Neumann boundary problem for axisymmetric fractional heat conduction equation in a solid with?cylindrical hole and associated thermal stress

The theory of thermoelasticity based on the heat conduction equation with the Caputo time-fractional derivative of order α is used to study thermal stress in an infinite medium with a cylindrical hole. Two types of Neumann boundary conditions are considered: the constant value of the normal derivative of the temperature and constant heat flux at the surface of a cavity. The solution is obtained applying Laplace and Weber integral transforms. Numerical results are illustrated graphically.     

AN ANALYSIS OF TRANSIENT HEAT CONDUCTION FOR FUNCTIONALLY GRADED SPHERES BASED ON STATE SPACE THEORY1)

基于状态空间理论研究功能梯度圆球的球对称瞬态热传导问题。根据热传导方程和热流密度的定义,取温度场和热流密度为系统的状态向量,通过将圆球分层和在时域内应用差分格式对控制方程进行离散,建立了系统的状态方程,给出了功能梯度圆球瞬态热传导问题的半解析解。算例分析表明：本文解不但结果正确、计算效率高,而且适用于材料参数沿径向任意梯度变化的圆球瞬态热传导分析。     

轴对称自由振动功能梯度材料微圆板中的热弹性阻尼

基于经典弹性薄板理论和单向耦合热传导理论,研究了材料性质沿厚度连续变化的功能梯度微圆板的热弹性阻尼特性．首先,考虑热力耦合效应,建立了功能梯度微圆板轴对称横向自由振动微分方程．然后,忽略温度梯度在面内的变化,建立了单向耦合变系数一维热传导方程．采用分层均匀化近似方法,将变系数热传导方程转化为一系列常系数的微分方程,利用上下表面的热边界条件和层间连续性条件获得了微圆板温度场解析解．将所得温度场代入微圆板的自由振动微分方程,得到了包含热弹性阻尼的复频率,从而获得了反映热弹性阻尼水平的逆品质因子．最后,针对材料性质沿板厚按幂函数变化的陶瓷-金属功能梯度微圆板,定量地分析材料梯度指数、几何尺寸、边界条件、温度环境等对微圆板热弹性阻尼的影响．     

A boundary identification method for an inverse heat conduction problem with an application in ironmaking

A boundary identification problem in inverse heat conduction is studied, based on data from internal measurement of temperature and heat flux. Formulated as a sideways heat conduction equation, a spatial continuation technique is applied to extend the solution to a known boundary condition at the desired boundary position. Recording the positions traversed in the continuation for each time instant yields the boundary position trajectory and hence the solution of the identification problem. A prospective application of the method can be found in the ironmaking blast furnace, where it is desired to monitor the thickness of the accreted refractory wall based on measurement of its internal state. Simulations featuring noisy measurement data demonstrate the feasibility of the identification method for blast furnace wall thickness estimation.     

Transient heat flow in a composite slab-constant flux, zero flux boundary conditions

Summary The time-dependent, one-dimensional equation of heat conduction is solved for a slab of two layers in perfect thermal contact. At one boundary there is a constant heat flux into the slab, and at the other boundary there is a zero flux. The solution for the temperature distribution is obtained with the aid of the Laplace transformation.This work was supported by the U.S. Naval Weapons Evaluation Facility, Albuquerque, N. M., U.S.A.