Integral transform solution for hyperbolic heat conduction in a finite slab

An analytical integral transformation of the thermal wave propagation problem in a finite slab is obtained through the generalized integral transform technique (GITT). The use of the GITT approach in the analysis of the hyperbolic heat conduction equation leads to a coupled system of second order ordinary differential equations in the time variable. The resulting transformed ODE system is then numerically solved by Gear's method for stiff initial value problems. Numerical results are presented for the local and average temperatures with different Biot numbers and dimensionless thermal relaxation times, permitting a critical evaluation of the technique performance. A comparison is also performed with previously reported results in the literature for special cases and with those produced through the application of the Laplace transform method (LTM), and the finite volume-Gear method (FVGM).     

An analytical solution for two-dimensional inverse heat conduction problems using Laplace transform

An analytical method has been developed for two-dimensional inverse heat conduction problems by using the Laplace transform technique. The inverse solutions are obtained under two simple boundary conditions in a finite rectangular body, with one and two unknowns, respectively. The method first approximates the temperature changes measured in the body with a half polynomial power series of time and Fourier series of eigenfunction. The expressions for the surface temperature and heat flux are explicitly obtained in a form of power series of time and Fourier series. The verifications for two representative testing cases have shown that the predicted surface temperature distribution is in good agreement with the prescribed surface condition, as well as the surface heat flux.     

Inverse heat conduction problems by using particular solutions

Based on the method of fundamental solutions, we develop in this paper a new computational method to solve two‐dimensional transient heat conduction inverse problems. The main idea is to use particular solutions as radial basis functions (PSRBF) for approximation of the solutions to the inverse heat conduction problems. The heat conduction equations are first analyzed in the Laplace transformed domain and the Durbin inversion method is then used to determine the solutions in the time domain. Least‐square and singular value decomposition (SVD) techniques are adopted to solve the ill‐conditioned linear system of algebraic equations obtained from the proposed PSRBF method. To demonstrate the effectiveness and simplicity of this approach, several numerical examples are given with satisfactory accuracy and stability. © 2011 Wiley Periodicals, Inc. Heat Trans Asian Res; Published online in Wiley Online Library ( wileyonlinelibrary.com ). DOI 10.1002/htj.20335     

Analytical approach with Laplace transform to the inverse problem of one‐dimensional heat conduction transfer: Application to second and third boundary conditions

An analytical method using Laplace transformation has been developed for one‐dimensional heat conduction. This method succeeded in explicitly deriving the analytical solution by which the surface temperature for the first kind of boundary condition can be well predicted. The analytical solutions for the surface temperature and heat flux are applied to the second and third of the boundary conditions. These solutions are also found to estimate the corresponding surface conditions with a high degree of accuracy when the surface conditions smoothly change. On the other hand, when these conditions erratically change such as the first derivative of temperature with time, the accuracy of the estimation becomes slightly less than that for a smooth condition. This trend in the estimation is similar irrespective of any kind of boundary condition. © 2002 Wiley Periodicals, Inc. Heat Trans Asian Res, 32(1): 29–41, 2003; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.10069     

A Novel Variational Formulation of Inverse Problem of Heat Conduction with Free Boundary on an Image

ＡＮｏｖｅｌＶａｒｉａｔｉｏｎａｌＦｏｒｍｕｌａｔｉｏｎｏｆＩｎｖｅｒｓｅＰｒｏｂｌｅｍｏｆＨｅａｔＣｏｎｄｕｃｔｉｏｎｗｉｔｈＦｒｅｅＢｏｕｎｄａｒｙｏｎａｎＩｍａｇｅＰｌａｎｅＧａｏ－ＬｉａｎＬｉｕ（ＳｈａｎｇｈａｉＩｎｓｔｉｔｕｔｅｏｆＭｅｃｈ...     

Improvement of inverse heat conduction solution using Laplace transformation: Method of partial division of time

Any solution for an inverse heat conduction problem makes the estimation of surface temperature and surface heat flux worsen in the case where these values behave like a triangular shape change with time. In order to compensate for this defect, Monde and colleagues, who succeeded in obtaining analytical inverse solutions using the Laplace transform technique, introduce a new idea where these changes over the entire measurement time can be split into several parts depending on the behavior. Therefore, an approximate equation to trace the measured temperature change can be derived, resulting in good estimation of surface temperature and surface heat flux even in the case of the triangular shape change and sharp change. © 2003 Wiley Periodicals, Inc. Heat Trans Asian Res, 32(7): 630–638, 2003; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.10117     

Analytical method of two‐dimensional inverse heat conduction problem using Laplace transformation: Effect of number of measurement points

An analytical method has been developed for the inverse problem of two‐dimensional heat conduction using the Laplace transform technique. The inverse problem is solved for only two unknown surface conditions and the other surfaces are insulated in a finite rectangular body. In actual temperature measurement, the number of points in a solid is usually limited so that the number of temperature measurements required to approximate the temperature change in the solid becomes too small to obtain an approximate function using a half polynomial power series of time and the Fourier series of the eigenfunction. In order to compensate for this lack of measurement points, the third‐order Spline method is recommended for interpolating unknown temperatures at locations between measurement points. Eight points are recommended as the minimum number of temperature measurement points. The calculated results for a number of representative cases indicate that the surface temperature and the surface heat flux can be predicted well, as revealed by comparison with the given surface condition. © 2003 Wiley Periodicals, Inc. Heat Trans Asian Res, 32(7): 618–629, 2003; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.10116     

Application of inverse solution in two‐dimensional heat conduction problem using Laplace transformation

An inverse solution has been explicitly derived for two‐dimensional heat conduction in cylindrical coordinates using the Laplace transformation. The applicability of the inverse solution is checked using the numerical temperatures with a normal random error calculated from the corresponding direct solution. For a gradual temperature change in a solid, the surface heat flux and temperature can be satisfactorily predicted, while for a rapid change in the temperature this method needs the help of a time partition method, in which the entire measurement time is split into several partitions. The solution with the time partitions is found to make an improvement in the prediction of the surface heat flux and temperature. It is found that the solution can be applied to experimental data, leading to good prediction. © 2003 Wiley Periodicals, Inc. Heat Trans Asian Res, 32(7): 602–617, 2003; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/htj.10115     

Finite element method formulation in polar coordinates for transient heat conduction problems

The aim of this paper is the formulation of the finite element method in polar coordinates to solve transient heat conduction problems. It is hard to find in the literature a formulation of the finite element method (FEM) in polar or cylindrical coordinates for the solution of heat transfer problems. This document shows how to apply the most often used boundary conditions. The global equation system is solved by the Crank-Nicolson method. The proposed algorithm is verified in three numerical tests. In the first example, the obtained transient temperature distribution is compared with the temperature obtained from the presented analytical solution. In the second numerical example, the variable boundary condition is assumed. In the last numerical example the component with the shape different than cylindrical is used. All examples show that the introduction of the polar coordinate system gives better results than in the Cartesian coordinate system. The finite element method formulation in polar coordinates is valuable since it provides a higher accuracy of the calculations without compacting the mesh in cylindrical or similar to tubular components. The proposed method can be applied for circular elements such as boiler drums, outlet headers, flux tubes. This algorithm can be useful during the solution of inverse problems, which do not allow for high density grid. This method can calculate the temperature distribution in the bodies of different properties in the circumferential and the radial direction. The presented algorithm can be developed for other coordinate systems. The examples demonstrate a good accuracy and stability of the proposed method.     

Analytical solution of the parabolic and hyperbolic heat transfer equations with constant and transient heat flux conditions on skin tissue

In this article, the parabolic (Pennes bioheat equation) and hyperbolic (thermal wave) bioheat transfer models for constant, periodic and pulse train heat flux boundary conditions are solved analytically by applying the Laplace transform method for skin as a semi-infinite and finite domain. The bioheat transfer analysis with transient heat flux on skin tissue has only been studied by Pennes equation for a semi-infinite domain. For modeling heat transfer in short duration of an initial transient, or when the propagation speed of the thermal wave is finite, there are major differences between the results of parabolic and hyperbolic heat transfer equations. The non-Fourier bioheat transfer equation describes the thermal behavior in the biological tissues better than Fourier equation. The outcome of transient heat flux condition shows that by penetrating into the depths beneath the skin subjected to heat, the amplitude of temperature response decreases significantly. The blood perfusion rate can be predicted using the phase shift between the surface temperature and transient surface heat flux. The thermal damage of the skin is studied by applying both the parabolic and hyperbolic bioheat transfer equations.     

A new numerical method to simulate the non-Fourier heat conduction in a single-phase medium

Many non-equilibrium heat conduction processes can be described by the macroscopic dual-phase lag model (DPL model). In this paper, a numerical method, which combines the dual reciprocity boundary element method (DRBEM) with Laplace transforms, is constructed to solve such mathematical equation. It is used to simulate the non-Fourier phenomenon of heat conduction in a single-phase medium, then numerically predict the differences between the thermal diffusion, the thermal wave and the non-Fourier heat conduction under different boundary conditions including pulse for one- and two-dimensional problems. In order to check this numerical method's reliability, the numerical solutions are still compared with two known analytical solutions.     

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

The non-Fourier heat conduction equation is solved for a finite medium under arbitrarily chosen surface disturbances by utilizing the principle of superposition along with the Fourier integral representation of arbitrary non-periodic functions. The obtained solution is such that for a given non-periodic disturbance, we just have to introduce the Fourier integral coefficients of that function to the presented general solution and perform some definite integrations, analytically if possible, but in general numerically which is a straightforward computational task.     

A technique for uncertainty analysis for inverse heat conduction problems

A technique is presented for the uncertainty analysis of the linear Inverse Heat Conduction Problem (IHCP) of estimating heat flux from interior temperature measurements. The selected IHCP algorithm is described. The uncertainty in thermal properties and temperature measurements is considered. A propagation of variance equation is used for the uncertainty analysis. An example calculation is presented. Parameter importance factors are defined and computed for the example problem; the volumetric heat capacity is the dominant parameter and an explanation is offered. Thoughts are presented on extending the analysis to include the non-linear problem of temperature dependent properties.     

Lack of oscillations in Dual-Phase-Lagging heat conduction for a porous slab subject to imposed heat flux and temperature

This study shows that the physical conditions necessary for thermal waves to materialize in Dual-Phase-Lagging porous media conduction are not attainable in a porous slab subject to a combination of constant heat flux and temperature (Neumann and Dirichlet) boundary conditions. It is demonstrated that the approximate equivalence between Dual-Phase-Lagging (DuPhlag) heat conduction model and the Fourier heat conduction in porous media subject to Lack of Local Thermal Equilibrium (La Lotheq) that suggested the possibility of thermal oscillations and resonance reveals a condition that cannot be fulfilled because of physical constraints.     

虚拟激励法求解激励随机热传导问题

为更有效地求解激励随机热传导问题,利用有限元技术和虚拟激励法,建立了激励随机热传导问题的数值求解模型,可考虑平稳与非平稳的随机激励。数值算例中的对比分析表明:通过虚拟激励技术,可有效计算真实响应的自/互功率谱密度和响应偏差。此外,对激励随机热传导逆问题进行了初步探讨。     

A finite element analysis for inverse heat conduction problems

This paper presents an efficient inverse analysis technique based on a sensitivity coefficient algorithm to estimate the unknown boundary conditions of multidimensional steady and transient heat conduction problems. Sensitivity coefficients were used to represent the temperature response of a system under unit loading conditions. The proposed method, coupled with the sensitivity analysis in the finite element formulation, is capable of estimating both the unknown temperature and heat flux on the surface provided that temperature data are given at discrete points in the interior of a solid body. Inverse heat conduction problems are referred to as ill-posed because minor inaccuracy or error in temperature measurements cause a drastic effect on the predicted surface temperature and heat flux. To verify the accuracy and validity of the new method, two-dimensional steady and transient problems are considered. Their surface temperature and heat flux are evaluated. From a comparison with the exact solution, the effects of measurement accuracy, number and location of measuring points, a time step, and regularization terms are discussed. © 1998 Scripta Technica. Heat Trans Jpn Res, 26(6): 345–359, 1997     

Peridynamic simulation of transient heat conduction problems in functionally gradient materials with cracks

In this article, a modified state-based peridynamic (PD) method is proposed to solve transient heat conduction problems in functionally gradient materials (FGMs) with extending insulated cracks. The PD formulation of transient heat conduction has been derived by using the time integration through the forward difference technique. Numerical simulations have been performed to verify the accuracy and effectiveness of the proposed method. The analytical solution and the finite element method results are used for comparison. In this work, the material properties of functionally graded materials are assumed to vary exponentially in z-direction. Our PD results show good agreement with analytical solutions and results from the finite element method. Hence the proposed PD method is suitable to deal with the transient heat conduction problem in FGMs with extending insulated cracks.     

Integral transform methods for inverse problem of heat conduction with known boundary of a thin rectangular object and its stresses

The three-dimensional inverse transient thermoelastic problem for a thin rectangular object is considered within the context of the theory of generalized thermoelasticity. The upper surface of the rectangular object occupying the space D: a≤x≤a; b≤y≤b; 0≤z≤h; with the known boundary conditions. Laplace and Finite Marchi-Fasulo transform techniques are used to determine the unknown temperature, temperature distribution, displacement and thermal stresses on upper plane surface of a thin rectangular object. The distributions of the considered physical variables are obtained and represented graphically.