3-Omega Method for Thin Films with the Non-Fourier Boundary Condition

The effect of non-Fourier boundary condition on the 3-omega method for measuring the thermal conductivity of microscale thin films using the hyperbolic heat conduction equation and the Fourier equation is examined. Non-Fourier boundary condition with the Fourier equation leads to 80% error in the temperature oscillations and increases the error to 85% in the case of non-Fourier boundary condition with the hyperbolic heat conduction equation. The solution of the non-Fourier boundary condition with the hyperbolic heat conduction equation gives the most accurate thermal conductivity expression. The analysis also provides a method for determining the relaxation time of thin films.     

Temperature Prediction for Tumor Hyperthermia with the Behavior of Thermal Wave

In magnetic nanoparticle hyperthermia for cancer treatment, controlling the heat distribution and temperature elevations is an immense challenge in clinical applications. It is expected for treatment quality to understand the heat transport occurring in biological tissue. The non-Fourier thermal behavior in biological tissue has been experimentally observed. This work uses the thermal wave model to predict the temperature excess occurring in a two-layer concentric spherical tissue with the heat source of Gaussian distribution. The solutions to the hyperbolic bio-heat equation with the space-dependent source term in the spherical coordinate system are presented. The influences of relaxation time, blood perfusion rate, and heating strength on the thermal response in tumor and normal tissue are discussed.     

BEHAVIOR OF THERMAL STRESSES IN A RAPIDLY HEATED THIN PLATE

Using the hyperbolic heat conduction model, thermal stresses generated within a rapidly heated thin metal plate are investigated numerically. The effects of different parameters such as the form, duration, amplitude, and penetration depth of the heating source on the temperature, thermal moment, deflection, and thermal stresses are studied. It is found that under ultra-fast heating of very thin plates, the hyperbolic heat conduction model must be adopted to model the thermal behavior.     

Ultra-fast pulse-laser heating on a two-layered semi-infinite material with interfacial contact conductance

The temperature distribution of a two-layered material in a semi-infinite domain subjected to ultra-fast pulse-laser heating at the front surface is solved in this study. The interfacial contact conductance existing at the interface is included in the analysis. The dual-phase-lag (DPL) model is applied to examine the lagging behavior of the two-layer material. The effects of interfacial contact conductance and the thickness of the surface layer on the temperature solution of the material are also investigated.     

SPECIALLY TAILORED TRANS FINITE-ELEMENT FORMULATIONS FOR HYPERBOLIC HEAT CONDUCTION INVOLVING NON-FOURIER EFFECTS

The phenomenon of hyperbolic heat conduction in contrast to the classical (parabolic) form of Fourier heat conduction involves thermal energy transport that propagates only at finite speeds as opposed to an infinite speed of thermal energy transport. To accommodate the finite speed of thermal wave propagation, a more precise form of heat flux law is involved, thereby modifying the heat flux originally postulated in the classical theory of heat conduction. As a consequence, for hyperbolic heat conduction problems, the thermal energy propagates with very sharp discontinuities at the wave front. The primary purpose of the present paper is to provide accurate solutions to a class of one-dimensional hyperbolic heat conduction problems involving non-Fourier effects that can precisely help understand the true response and furthermore can be used effectively for representative benchmark tests and for validating alternate schemes. As a consequence, the present paper purposely describes modeling/analysis formulations via specially tailored hybrid computations for accurately modeling the sharp discontinuities of the propagating thermal wave front. Comparative numerical test models are presented for various hyperbolic heat conduction models involving non-Fourier effects to demonstrate the present formulations.     

Non-Fourier heat conduction in a hollow sphere with periodic surface heat flux

Presented is the analysis of non-Fourier effect in a hollow sphere exposed to a periodic boundary heat flux. The problem is studied by deriving an analytical solution of the hyperbolic heat conduction equation. Using the obtained analytical expression, the temperature profiles at outer and inner surfaces of the sphere are evaluated for various thermal relaxation times. By comparing the results of non-Fourier model with those obtained from Fourier heat conduction equation, the transition process from parabolic model to hyperbolic one is shown. The phase difference and amplitude ratio of boundary surfaces are calculated as functions of the thermal relaxation time and the results are depicted graphically.     

Transient hygrothermoelastic response in a porous cylinder subjected to ramp-type heat-moisture loading

Abstract

The classical theory of heat conduction (Fourier theory) predicts an infinite speed for thermal disturbance propagation, which is physically unrealistic. By extending the classical Fourier heat conduction and Fick’s diffusion, this article develops hyperbolic diffusion/heat conduction laws with phase lags of heat/moisture flux to simulate coupled heat-moisture diffusion-propagation behavior with the Defour and Soret effects. A porous cylinder subjected to a ramp-type heating and humidifying at the surface is studied. The Laplace transform is used to obtain a closed-form solution of the temperature, moisture, displacements and stresses in the cylinder. Numerical results are calculated via the inversion of the Laplace transform. Obtained results show that the thermal/moisture relaxation time or phase lag plays a significant role in affecting transient hygrothermoelastic field. For a non-vanishing phase lag, non-Fourier and non-Fickian effects exist and hygrothermal waves have finite propagation speeds. The influences of the phase lag of heat/moisture flux and ramp-type time parameter on the transient response of hygrothermoelastic field are presented graphically. A comparison of the numerical results based on the classical model and the present one is made. The non-Fourier heat conduction and non-Fickian diffusion can effectively avoid the shortcomings induced by the classical Fourier and Fick laws.     

Analysis of transient heat transfer in multilayer thin films with nonlinear thermal boundary resistance

Transient energy transport in thin-layer films with a nonlinear thermal boundary resistance is analyzed theoretically within the framework of the dual-phase-lag heat conduction model. An iterative finite difference numerical method is used and is verified using a derived semi-analytical solution of the problem. Effects of the thermo-physical properties on energy transport when a two-layer film is exposed to a thermal pulse of certain duration and strength are presented. The thermal boundary resistance, the heat flux and temperature gradient phase lags and the thermal conductivities and heat capacities all are important factors that characterize energy transport through the interface and the temperature distribution in the two layers. The maximum interfacial temperature difference that takes place in the transient process of thermal pulse propagation is found to be the proper choice to measure the perfect-ness of the interface with a finite thermal boundary resistance. The results show that even with high values of the thermal boundary resistance the maximum interfacial temperature difference can be very small when the thermal pulse propagates from a high-thermal conductivity and heat capacity layer to a low-thermal conductivity and heat capacity layer. For a certain range of the thermal conductivities and heat capacities, the maximum interfacial temperature difference approaches zero even with high values of the thermal boundary resistance. Thermal conductivities and heat capacities are much more important in characterizing transient heat transfer through the imperfect interface than the phase lags of the heat flux and temperature gradient.     

NON-FOURIER HEAT CONDUCTION PHENOMENA IN POROUS MATERIAL HEATED BY MICROSECOND LASER PULSE

Experiments on porous material heated by a microsecond laser pulse and the corresponding theoretical analysis are carried out. Some non-Fourier heat conduction phenomena are observed in the experimental sample. The experimental results indicate that only if the thermal disturbance is strong enough (i.e., the pulse duration is short enough and the pulse heat flux is great enough) is it possible to observe apparent non-Fourier heat conduction phenomenon in the sample, and evident non-Fourier heat conduction phenomenon can only exist in a very limited region around the thermal disturbance position. The hyperbolic heat conduction (HHC) equation and the dual-phase lag (DPL) model are employed, respectively, to describe the non-Fourier heat condution process happening in the experimental sample, and the finite-difference method (FDM) is used to solve them numerically. The numerical solutions show that both the HHC equation and the DPL model can predict the non-Fourier heat conduction phenomenon emerging in the experimental sample qualitatively. Moreover, if τ_q and τ_T are assumed to have suitable values, the theoretical result of the DPL model is more agreeable to the experimental result.     

Analysis of Non-Fourier Conduction-Radiation Heat Transfer in a Cylindrical Enclosure

This article deals with the analysis of non-Fourier conduction and radiation heat transfer in a participating medium contained between 1-D concentric cylinders. The conducting-radiating medium is radiatively absorbing, emitting, and scattering. The non-Fourier effect is analyzed by suddenly perturbing the temperatures of the concentric cylinders. With radiative information computed using the finite volume method, the finite difference method is used to solve the hyperbolic energy equation. Effects of various parameters such as the extinction coefficient, the scattering albedo, the conduction-radiation parameter, the emissivity and the radius ratio are studied on the temporal evolution of temperature field in the medium. These parameters have been found to significantly influence the temporal temperature field, and non-Fourier effects are captured well. For non-Fourier conduction and Fourier conduction–radiation cases, results have been benchmarked against those available in the literature. A good comparison has been observed. In case of non-Fourier conduction-radiation, for some sample cases, the steady-state temperature distributions have been compared against those available in the literature. Results have been found to agree well.     

Slip Boundary Conditions in Ballistic–Diffusive Heat Transport in Nanostructures

Ballistic–diffusive heat conduction, which is predominantly affected by boundaries and interfaces, will occur in nanostructures whose characteristic lengths are comparable to the phonon mean free path (MFP). Here, we demonstrated that interactions between phonons and boundaries (or interfaces) could lead to two kinds of slip boundary conditions in the ballistic–diffusive regime: boundary temperature jump and boundary heat flux slip. The phonon Boltzmann transport equation (BTE) with relaxation time approximation and the phonon tracing Monte Carlo (MC) method were used to investigate these two slip boundary conditions for the ballistic–diffusive heat conduction in nanofilms on a substrate. For cross-plane heat conduction where the boundary temperature jump is the dominant non-Fourier phenomenon, ballistic transport causes the temperature jumps and thus introduces a ballistic thermal resistance. Importantly, when considering the interface effect, the corresponding model was derived based on the phonon BTE and verified by comparing with the MC simulations. In addition, an interface–ballistic coupling effect was identified, which indicates inapplicability of the standard thermal resistance analysis. In contrast, for the in-plane case that is controlled by boundary heat flux slip, both phonon boundary scattering and perturbation of the phonon distribution function induced by the interface can cause heat flux slip, leading to a variation in in-plane thermal resistance. In addition, a model beyond the Fuchs-Sondheimer formula, which can address both the boundary scattering and the interface effects, was derived based on the phonon BTE. The good agreements with the MC simulations indicate its validity.     

Dual-Phase-Lag Heat Conduction in Finite Slabs With Arbitrary Time-Dependent Boundary Heat Flux

Applying a constant or transient heat flux on a plane slab is a common technique in microelectronics technology and material processing, including laser patterning, micromachining, and laser surface treatment processes. Although Fourier's law is typically very precise for evaluating temperatures in solids, a number of experimental observations suggest the existence of non-Fourier transient conduction in these applications. Since the dual-phase-lag (DPL) model of heat conduction can be compatible with the hypothesis of local equilibrium thermodynamics (as shown here), the effects of temperature gradient relaxation time on the non-Fourier hyperbolic conduction in a finite slab subjected to an arbitrary time-dependent surface heat flux is examined by this model. The combination of diffusion- and wave-like features in heat conduction process is properly monitored by the DPL model for two types of heat flow regimes, namely, gradient precedence and flux precedence. The results indicate considerable deviations between the predictions of these regimes.     

RAPID TRANSIENT HEAT CONDUCTION IN MULTILAYER MATERIALS WITH PULSED HEATING BOUNDARY

ABSTRACT

Rapid transient heat conduction in multilayer materials under pulsed heating is solved numerically based on a hyperbolic heat conduction equation and taking into consideration the non-Fourier heat conduction effects. An implicit difference scheme is presented and a stability analysis conducted, which shows that the implicit scheme for the hyperbolic equation is stable. The code is validated by comparing the numerical results with an existing exact solution, and the physically unrealistic conditions placed on the time and space increments are identified. Using the validated model, the numerical solution of thermal wave propagation in multilayer materials is presented. By analyzing the results, the necessary conditions for observing non-Fourier phenomena in the laboratory can be inferred. The results are also compared with the numerical results from the parabolic heat conduction equation. The difference between them is clearly apparent, and this comparison provides new insight for the management of thermal issues in high-energy equipment. The results also illustrate the time scale required for metal films to establish equilibrium in energy transport, which makes it possible to determine a priori the time response and the measurement accuracy of metal film, thermal-resistant thermometers.     

汽轮机转子温度场和热应力的动态分析

将空心的汽轮机转子简化成长圆筒壁模型,在一定的假设条件下得到了描述其导热特性的微分方程。利用复频域分析方法求解了长圆筒壁的一维动态导热特性,并根据热弹性理论得到了转子温度场和热应力的传递函数。进而对内表面绝热、外表面为阶跃温度输入和内、外表面都为阶跃温度输入的两种情况下转子温度场和热应力的变化进行了计算和分析。计算结果表明,采用对转子内、外表面同时加热或冷却的方法,能够在不增加热应力的条件下,加速转子达到温度平衡状态,有利于缩短汽轮机组的开停机时间,并提高跟踪机组负荷变化的灵活性。     

Non-Fourier heat conduction in a finite medium subjected to arbitrary periodic surface disturbance

The non-Fourier transient heat conduction in a finite medium under arbitrary periodic surface thermal disturbance is investigated analytically. In order to obtain the desired temperature field from the known solution for non-Fourier heat conduction under a harmonic disturbance, the principle of superposition along with the Fourier series representation of an arbitrary periodic function is employed. The developed method can be applied for more realistic periodic boundary conditions occurred in nature and technology.     

Finite-Volume Solution of a Hyperbolic Two-Step Conduction Equation

The finite-volume method is developed for the solution of two-step hyperbolic conduction equations and applied to fast transient non-Fourier conduction problems in which two energy carriers are not in thermal equilibrium. Numerical methodology has been formulated to be applied to heterogeneous materials with temperature-dependent properties. Application to ultrafast laser heating of thin composite metals shows favorable results in comparison with the existing data.     

Semi-analytical solutions for the dual phase lag heat conduction in multilayered media

A semi-analytical solution procedure for transient heat transfer in composite mediums consisting of multi-layers within the framework of the dual phase lag model is presented. The procedure is then used to derive solutions for the temperature-, temperature gradient-, and heat flux distributions in a two-layer composite planar slab, a bi-layered solid-cylinder and sphere. The solutions obtained are applicable to the classical Fourier heat diffusion, hyperbolic heat conduction, phonon–electron interaction, and phonon scattering models with perfect or imperfect contact and with layers of different materials. The interfacial contact resistance, the heat flux and temperature gradient phase lags, thermal diffusivities and conductivities, initial temperatures of the composite medium and a general time-dependent boundary heat flux enter the solutions as parameters, allowing the solutions obtained to be applicable to a wide range of arrangements including perfect and imperfect contact. Analysis of thermal wave propagation, transmission and reflection in planar, cylindrical and spherical geometries with imperfect interfaces are presented, and geometrical—as well as the temperature gradient phase lag—effects on the thermal lagging behavior in different layered media are discussed.     

Notable physical anomalies manifested in non-Fourier heat conduction under the dual-phase-lag model

In this study, a number of notable physical anomalies concerning non-Fourier heat conduction under the dual-phase-lag (DPL) model are observed and investigated. It is found that, during the transient heat transfer process, the over-diffusion mode predicts a “hyper-active” to “under-active” transition in thermal behavior. The main cause behind it lies in the time-varying effect of τ_T (the phase lag of the temperature gradient) on the thermal response. Also, change of polarity in reflected thermal waves can be observed when a constant-temperature boundary is involved, which hints that a heating process may be followed by a spontaneous cooling effect. A fairly strong connection is present between the τ_T-induced dispersive effect and an unusual thermal accumulation phenomenon in an on–off periodic heating process. Furthermore, a paradox involving a moving medium is detected in the DPL model, which can be solved by replacing the temporal partial derivatives in the DPL equation with the material derivatives. During the process of analysis, a high-order characteristics-based TVD scheme is relied on to provide accurate and reliable numerical simulations to the DPL heat conduction equation under various initial-boundary conditions.     

Thermodynamic stability analysis of non-Fourier model-based bio-heat transfer in laser-irradiated cylindrical-shaped multilayered biological tissue

The present research focuses on examining the thermic response of living tissue in the form of a triple-layered cylindrical structure when subjected to laser light and the compatibility analysis of non-Fourier heat transfer with thermodynamics second law. The temperature field in the triple-layered cylindrical living tissue subjected to laser light is determined by numerically solving the transient radiative transfer equation in conjunction with the dual phase lag (DPL) based bioheat equation. Once the temperature field is known, the equilibrium and nonequilibrium entropy production rate (EPR) is calculated based on the hypothesis of classical irreversible thermodynamics and extended irreversible thermodynamics, respectively. The present results are verified against the data from the literature and found a good match between them. A comparative analysis of the Fourier and non-Fourier models is accomplished. The equilibrium and nonequilibrium EPR values for the Fourier model are found to be positive. While the equilibrium EPR is negative for non-Fourier heat conduction and does not satisfy the thermodynamics second law, nonequilibrium EPR is always a positive value for Fourier, DPL, and hyperbolic models and satisfies thermodynamics second law. It has been investigated how thermal relaxation times affect the temperature field and EPRs in tissue are subjected to laser light.