考虑非傅立叶效应的瞬态导热过程的波动响应分析

瞬态导热分析需要考虑非傅立叶效应.通过对抛物型及双曲型热传导方程,以及耦合热传导方程后的波动方程的数值求解,研究了非傅立叶效应下导热过程中的波动响应.结果表明,双曲型热传导过程具有明显的波动特性,所引起的波动响应前沿值成倍提高,且呈现显著的跃变行为,而波动峰值外的部位围绕着初始值小幅波动.     

Analytical and numerical analysis of time fractional dual-phase-lag heat conduction during short-pulse laser heating

Numerical Algorithms - In this study, we analytically and numerically investigate the non-Fourier heat conduction behavior within a finite medium based on the time fractional dual-phase-lag model....     

A superstatistical model for anomalous heat conduction and diffusion

We propose a superstatistical model for anomalous heat conduction and diffusion, which is formulated by the thermal conductivity distribution, overall temperature and heat flux distributions. Our model obeys Fourier's law and the continuity equation at the individual level. The evolution of the thermal conductivity distribution is described by an advection-diffusion equation. We show that the superstatistical model predict anomalous behaviors including the time-dependent effective thermal conductivity and slow long-time asymptotics. The time-dependence of the effective thermal conductivity is determined by the mean square displacement (MSD), which coincides with existing investigations. The superstatistical structure can also be extended into other non-Fourier models including the Cattaneo and fractional-order heat conduction models.     

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.     

The estimation of the strength of the heat source in the heat conduction problems

A general method is proposed to determine the strength of the heat source in the Fourier and non-Fourier heat conduction problems. A finite difference method, the concept of the future time and a modified Newton–Raphson method are adopted in the problem. The undetermined heat source at each time step is formulated as an unknown variable in a set of equations from the measured temperature and the calculated temperature. Then, an iterative process is used to solve the set of equations. No selected function is needed to represent the undetermined function in advance. Three examples are used to demonstrate the characteristics of the proposed method. The validity of the proposed method is confirmed by the numerical results. The results show that the proposed method is an accurate and stable method to determine the strength of the heat source in the inverse hyperbolic heat conduction problems. Furthermore, the result shows that more future times are needed in the hyperbolic equation than that of parabolic equation. Moreover, the robustness and the accuracy of the estimated results in the non-Fourier problem are not as well as those of the Fourier problem.     

Nonlinear thermoelastic analysis of layered structures

The nonlinear thermoelastic behavior of orthotropic layered slabs and cylinders including radiation boundaries, temperature-dependent material properties, and stress-dependent layer interface conditions is investigated. A one-dimensional finite element formulation employing quadratic layer and linear interface elements is used to perform the analyses. The transient heat conduction portion of the program is temporally discretized via an implicit linear time interpolation algorithm which includes Crank-Nicolson, Galerkin, and Euler backward differencing. The nonlinear heat conduction equations are iteratively evaluated using a modified Newton-Raphson scheme. Direct iteration between heat conduction and stress analysis is employed when stress-dependent interface thermal resistance coefficients are utilized. Verification problems are presented to demonstrate the accuracy of the finite element code.     

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.     

广义热传导方程有限元算法的计算准则

本文作者曾对经典的(抛物型)热传导方程提出了两种单调性的新概念,推导并证明了几组计算准则,可以使其有限元数值解消除很容易出现的振荡和超界现象.本文把上述成果用于广义(双曲型)热传导方程的有限元解中,推导出它的有限元解的计算准则,并获得了一些新结论.     

Inverse determination of thermal conductivity using semi-discretization method

In this work a semi-discretization method is presented for the inverse determination of spatially- and temperature-dependent thermal conductivity in a one-dimensional heat conduction domain without internal temperature measurements. The temperature distribution is approximated as a polynomial function of position using boundary data. The derivatives of temperature in the differential heat conduction equation are taken derivative of the approximated temperature function, and the derivative of thermal conductivity is obtained by finite difference technique. The heat conduction equation is then converted into a system of discretized linear equations. The unknown thermal conductivity is estimated by directly solving the linear equations. The numerical procedures do not require prior information of functional form of thermal conductivity. The close agreement between estimated results and exact solutions of the illustrated examples shows the applicability of the proposed method in estimating spatially- and temperature-dependent thermal conductivity in inverse heat conduction problem.     

四段式有限元法分析热耦合的流体-固体相互作用问题

给出了一种流(体)-热-结构综合的分析方法,固体中的热传导耦合了粘性流体中的热对流,因而在固体中产生热应力.应用四段式有限元法和流线逆风Petrov-Galerkin法分析热粘性流动,应用Galerkin法分析固体中的热传导和热应力.应用二阶半隐式Crank-Nicolson格式对时间积分,提高了非线性方程线性化后的计算效率.为了简化所有有限元公式,采用3节点的三角形单元,对所有的变量:流体的速度分量、压力、温度和固体的位移,使用同阶次的插值函数.这样做的主要优点是,使流体-固体介面处的热传导连接成一体.数个测试问题的结果表明,这种有限元法是有效的,且能加深对流(体)-热-结构相互作用现象的理解.     

An inverse non-Fourier heat conduction problem approach for estimating the boundary condition in electronic device

This study is intended to provide an inverse method for estimating the unknown boundary condition T(0,y,t) in a non-Fourier heat conduction electronic device. In this study, finite-difference methods are employed to discretize the problem domain, and then a linear inverse model is constructed to identify the unknown boundary condition. The present approach is to rearrange the matrix forms of the differential governing equations and to estimate the unknown conditions. Then, the linear least-squares method is adopted to obtain the solution.The results show that one measuring point is sufficient to estimate the unknown boundary condition T(0,y,t) without measurement errors. When considering the measurement errors, the magnitudes of the discrepancies in the boundary condition T(0,y,t) are directly proportional to the size of measurement errors. Due to the complicated reflection and interaction of the thermal waves, this phenomenon reflects the fact that the inverse non-Fourier heat conduction problem is different from the inverse Fourier heat conduction problem.In contrast to the traditional approach, the advantage of applying this method in inverse analysis is that no prior information is needed on the functional form of the unknown quantities. In addition, no initial guess is required and the calculation can be done in only one iteration.     

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.     

非Fourier温度场分布的奇摄动解

应用非Fourier热传导定律构建了单层材料中温度场模型，即一类在无界域上带小参数的奇摄动双曲方程，通过奇摄动展开方法，得到了该问题的渐近解.首先应用奇摄动方法得到了该问题的外解和边界层矫正项，通过对内解和外解的最大模估计和关于时间导数的最大模估计以及线性抛物方程理论，得到了内外解的存在唯一性，从而得到了解的形式渐近展开式.通过余项估计，给出了渐近解的L2估计，得到了渐近解的一致有效性，从而得到了无界域上温度场的分布.通过奇摄动分析，给出了非Fourier 温度场与Fourier 温度场的关系，描述了非Fourier温度场的具体形态.     

Use of geometry in finite element thermal radiation combined with ray tracing

In heat transfer for space applications, the exchanges of energy by radiation play a significant role. In this paper, we present a method which combines the geometrical definition of the model with a finite element mesh. The geometrical representation is advantageous for the radiative component of the thermal problem while the finite element mesh is more adapted to the conductive part. Our method naturally combines these two representations of the model. The geometrical primitives are decomposed into cells. The finite element mesh is then projected onto these cells. This results in a ray tracing acceleration technique. Moreover, the ray tracing can be performed on the exact geometry, which is necessary if specular reflectors are present in the model. We explain how the geometrical method can be used with a finite element formulation in order to solve thermal situation including conduction and radiation. We illustrate the method with the model of a satellite.     

Dual-phase-lag analysis of CNT–MoS2–ZrO2–SiO2–Si nano-transistor and arteriole in multi-layered skin

The numerical simulation of non-Fourier dual-phase-lag heat conduction in three-dimensional bio and nano media has been investigated with two developed solvers using the OpenFOAM C++ finite-volume libraries. Despite the rapid growth in the use of non-Fourier heat transfer model, most researchers conduct simulations over simple geometries due to numerical restrictions. The developed solvers have been verified for four benchmark cases. After that, two complex geometries have been selected to show the accuracy of the solvers for complicated geometries. The first case is the temperature distribution in the human hand skin in the presence of a circular arteriole, while the second case is the thermal analysis of a newly constructed CNT–MoS₂ transistor with 1 nm gate length. Both cases are multilayered structures with curved boundaries. The computed temperature rises in the blood vessel and in the transistor are 63°C after 15 s and 16 °C after 10 ps, respectively. It is found that considering the temperature-dependent electrical properties of ZrO₂ and CNT is essential to present an accurate mathematical model.     

考虑记忆效应及尺寸效应窄长薄板的磁-热弹性耦合动态响应

引入记忆依赖微分的双相滞后热弹性理论能较完善地描述非Fourier导热现象,然而迄今尚未发现该理论综合考虑微尺度效应和磁、热、弹等多场耦合效应对材料力学行为的影响.通过考虑记忆依赖效应和非局部效应修正了双相滞后广义热弹性理论,基于改进后的理论研究了受周期性变化热源作用时窄长薄板的磁-热弹性耦合问题.首先建立问题的控制方程;然后结合边界条件与初值条件,利用Laplace变换和反变换技术对该问题进行求解;最后分别考察了磁场、相位滞后、时间延迟因子、核函数、非局部效应、时间对各无量纲量的影响,为微尺度材料的动态响应提供了有力参考依据.     

The thermal problem for a submerged stream

The general solution of the thermal problem of convective heat conduction with volume heat dissipation caused by viscous dissipation of kinetic energy of the fluid, whose velocity field is determined by the exact solution /1/ of the Navier-Stokes equations, is considered for a submerged stream. The possible formulation of the heat problem and the characteristic behaviour of the solutions are investigated. The solutions obtained have a special feature, namely the existence, under specified conditions, of two regimes of convective heat exchange.     

Numerical analysis of quenching – Heat conduction in metallic materials

Metallic materials present a complex behavior during heat treatment processes. In a certain temperature range, change of temperature induces a phase transformation of metallic structure, which alters physical properties of the material. Indeed, measurements of specific heat and conductivity show strong temperature-dependence during processes such as quenching of steel. Several mathematical models, as solid mixtures and thermal–mechanical coupling, for problems of heat conduction in metallic materials, have been proposed. In this work, we take a simpler approach without thermal–mechanical coupling of deformation, by considering the nonlinear temperature-dependence of thermal parameters as the sole effect due to those complex behaviors. The above discussion of phase transformation of metallic materials serves only as a motivation for the strong temperature-dependence as material properties. In general, thermal properties of materials do depend on the temperature, and the present formulation of heat conduction problem may be served as a mathematical model when the temperature-dependence of material parameters becomes important. For this mathematical model we present the error estimate using the finite element method for the continuous-time case.     

Numerical analysis of a dual-phase-lag model involving two temperatures

In this paper, we numerically analyse a phase-lag model with two temperatures which arises in the heat conduction theory. The model is written as a linear partial differential equation of third order in time. The variational formulation, written in terms of the thermal acceleration, leads to a linear variational equation, for which we recall an existence and uniqueness result and an energy decay property. Then, using the finite element method to approximate the spatial variable and the implicit Euler scheme to discretize the time derivatives, fully discrete approximations are introduced. A discrete stability property is proved, and a priori error estimates are obtained, from which the linear convergence of the approximation is derived. Finally, some one-dimensional numerical simulations are described to demonstrate the accuracy of the approximation and the behaviour of the solution.