Determination of the limit of applicability of the parabolic equation of heat conduction

A criterion is proposed which enables breakdown of the parabolic model of heat conduction for isotropic materials, in which there are no phase transitions, to be revealed, and which defines the conventional boundary between the linear and nonlinear regimes of nonequilibrium thermodynamics in heat problems. __________ Translated from Izmeritel’naya Tekhnika, No. 6, pp. 38–42, June, 2008.     

Effect of radiation heat transfer on thermal diffusivity measurements

Experimental data on thermal conductivity and thermal diffusivity of a semitransparent material generally include an error due to the radiation heat transfer. This error varies in accordance with the experimental conditions such as the temperature level of the sample and the measuring method. In this paper, research on the influence of radiation heat transfer on thermal diffusivity are reviewed, and as an example, the method to correct the radiation component in the apparent thermal diffusivity measured by the stepwise heating technique is presented. The transient heat transfer by simultaneous thermal conduction and radiation in a semitransparent material is analyzed when the front surface is subjected to stepwise heating. The apparent thermal diffusivity, which includes the radiation component, is calculated for various parameters.Paper presented at the Second U.S.-Japan Joint Seminar on Thermophysical Properties, June 23, 1988, Gaithersburg, Maryland, U.S.A.     

A hybrid flux model for heat conduction problems

A hybrid method of solution for the linear problem of heat conduction in a body is presented. The variational support is a two‐field functional whose arguments are heat flux, which meets a priori inner thermal equilibrium, and temperature on the boundary of the body. The stationary conditions of the functional are the Fourier's law and the prescribed boundary conditions. This variational framework allows to develop a finite element model that exhibits good accuracy, especially in the presence of geometry irregularities in a mesh. Copyright © 2000 John Wiley & Sons, Ltd.     

Measurement of the thermal diffusivity by the thermal wave method using low-power heat sources

The errors in measuring the thermal diffusivity by the plane thermal wave method are considered as a function of the thermal flux power density. The minimum values of the thermal flux power density required for measurements with a specified error and the optimum parameters of the samples and of the heat source are determined. __________ Translated from Izmeritel’naya Tekhnika, No. 8, pp. 44–46, August, 2007.     

Direct error in constitutive equation formulation for inverse heat conduction problem

We present a mixed numerical formulation that handles discontinuities well for scalar hyperbolic partial differential equations. The formulation is based on a least‐square error in the constitutive equation. It is motivated by scalar inverse diffusion problems with interior data and applies to convection of a passive scalar in a discontinuous compressible flow field. We motivate the need for a mixed formulation by discretizing using an irreducible finite element method and discuss some of the limitations of that approach. We then develop and prove that the mixed formulation is well posed and verify that it works for problems with continuous and discontinuous thermal conductivity distributions.     

A finite element method for non-Fourier heat conduction in strong thermal shock environments

Non-Fourier effect is important in heat conduction in strong thermal environments. Currently, generally-purposed commercial finite element code for non-Fourier heat conduction is not available. In this paper, we develop a finite element code based on a hyperbolic heat conduction equation, which includes the non-Fourier effect in heat conduction. The finite element space discretization is used to obtain a system of differential equations for the time. The transient responses are obtained by solving the system of differential equations, based on the finite difference, mode superposition, or exact time integral. The code is validated by comparing the numerical results with exact solutions for some special cases. The stability analysis is conducted and it shows that the finite difference scheme is an ideal method for the transient solution of the temperature field. It is found that with mesh refining (decreasing mesh size) and/or high-order elements, the oscillation in the vicinity of sharp change vanishes, and can be essentially suppressed by the finite difference scheme. A relationship between the time step and the space length of the element was identified to ensure that numerical oscillation vanishes.     

Use of the Laplace transformation for data reduction in the flash method of measuring thermal diffusivity

This paper presents a new way to reduce data in the laser flash method of measuring thermal diffusivity. Experimental temperature vs time data are first transformed by using the Laplace transformation, and then they are fitted with an appropriate theoretical formula. The data reduction procedure is more efficient and enables the use of more realistic models of heat conduction in the sample, because the theoretical formulae for transformed temperatures have a simpler form than those for nontransformed ones. Some examples of the theoretical formulae of transformed temperatures are included here for one- and two-dimensional heat transfer, respectively. The models described take into account a finite pulse time and heat losses from the sample. Two fitting algorithms are proposed. Experimentally, the data reduction procedure has been tested for a correction of the finite pulse time effect in the flash method. The results show that the accuracy of our procedure is comparable with other data reduction methods. Provided that the shape and duration of the pulse are known, this procedure allows elimination of the finite pulse time effect on calculation of the thermal diffusivity for any transformable heat pulse time function, even in cases where the other specialized data reduction procedures have failed.     

Global space–time multiquadric method for inverse heat conduction problem

In this paper, a radial basis collocation method (RBCM) based on the global space–time multiquadric (MQ) is proposed to solve the inverse heat conduction problem (IHCP). The global MQ is simply constructed by incorporating time dimension into the MQ function as a new variable in radial coordinate. The method approximates the IHCP as an over‐determined linear system with the use of two sets of collocation points: one is satisfied with the governing equation and another is for the given conditions. The least‐square technique is introduced to find the solution of the over‐determined linear system. The present work investigates two types of the ill‐posed heat conduction problems: the IHCP to recover the surface temperature and heat flux history on a source point from the measurement data at interior locations, and the backward heat conduction problem (BHCP) to retrieve the initial temperature distribution from the known temperature distribution at a given time. Numerical results of four benchmark examples show that the proposed method can provide accurate and stable numerical solutions for one‐dimensional and two‐dimensional IHCP problems. The sensitivity of the method with respect to the measured data, location of measurement, time step, shape parameter and scaling factor is also investigated. Copyright © 2010 John Wiley & Sons, Ltd.     

Boundary element analysis for axisymmetric heat conduction and thermal stress in steady state

In this paper, axisymmetric heat conduction and thermal stress problems with three types of boundary conditions are analysed by the boundary element method. The temperature and thermal stress fields for the piston of a diesel engine are calculated using triangular finite elements and constant boundary elements, respectively, and the two results agree. However, BEM needs fewer data, less computer time (about one-sixth that of FEM) and storage volume. The advantages of BEM are sufficiently demonstrated.     

Application of photoacoustic and photothermal techniques for heat conduction measurements in a free-standing chemical vapor-deposited diamond film

Heat conduction in a free-standing chemical vapor-deposited polycristalline diamond film has been investigated by means of combined front and rear photoacoustic signal detection techniques and also by means of a “mirage” photothermal beam deflection technique. The results obtained with the different techniques are consistent with a value of α=(5.5±0.4)×10⁻⁴ m² · s⁻¹ for thermal diffusivity, resulting in a value ofκ=(9.8±0.7)×10² W·m⁻¹·K⁻¹ for thermal conductivity when literature values for the density and heat capacity for natural diamond are used.     

Application of cubic Hermitian algorithms to boundary element analysis of heat conduction

This paper presents the application of cubic Hermitian interpolation based finite element schemes for the time integration of the differential-algebraic system arising in the dual reciprocity boundary element formulation of transient diffusion problems. Weighted residual procedure is used to obtain the desired recurrence relations. Numerical results presented for three representative problems involving different types of boundary conditions amply demonstrate the high accuracy of the cubic Hermitian schemes.     

Determination of the thermal diffusivity and specific heat using an exponential heat pulse,including heat-loss effects

The one-dimensional heat diffusion equation has been solved analytically for the case of a heat pulse of the form F(t) = exp(–t/)/ applied to the front face of a homogeneous body including the effects of heat loss from the front and back faces. Approximate expressions are presented which yield a simple, accurate technique for the determination of the thermal diffusivity and specific heat, suitable to a wide range of heat-pulse time constant and heat-loss parameters, without recourse to graphical techniques or requiring further computer analysis. A procedure is described for the determination of an effective time constant to allow application of the present results to the case of a nonexponential heat pulse. Experimental results supporting the theoretical analysis are presented for five samples of silicon germanium alloys of various thicknesses, determined using a xenon flash tube heat-pulse exhibiting an exponential dependence. Proper consideration of the experimental heat pulse shape is shown to lead to reliable corrections to the apparent thermal diffusivity, even for relatively long heat-pulse times.Paper presented at the Tenth Symposium on Thermophysical Properties, June 20–23, 1988, Gaithersburg, Maryland, U.S.A.     

聚合物基复合材料的导热性能研究

讨论了聚合物基高导热高绝缘纳米复合材料的导热机理与常用的导热理论模型。考虑到填充率、温度等的影响，用不同的理论模型计算了氧化铝纳米颗粒填充环氧树脂的热导率，并结合相关研究实验对不同的导热理论模型进行分析比较。     

A non‐iterative finite element method for inverse heat conduction problems

A non‐iterative, finite element‐based inverse method for estimating surface heat flux histories on thermally conducting bodies is developed. The technique, which accommodates both linear and non‐linear problems, and which sequentially minimizes the least squares error norm between corresponding sets of measured and computed temperatures, takes advantage of the linearity between computed temperatures and the instantaneous surface heat flux distribution. Explicit minimization of the instantaneous error norm thus leads to a linear system, i.e. a matrix normal equation, in the current set of nodal surface fluxes. The technique is first validated against a simple analytical quenching model. Simulated low‐noise measurements, generated using the analytical model, lead to heat transfer coefficient estimates that are within 1% of actual values. Simulated high‐noise measurements lead to h estimates that oscillate about the low‐noise solution. Extensions of the present method, designed to smooth oscillatory solutions, and based on future time steps or regularization, are briefly described. The method's ability to resolve highly transient, early‐time heat transfer is also examined; it is found that time resolution decreases linearly with distance to the nearest subsurface measurement site. Once validated, the technique is used to investigate surface heat transfer during experimental quenching of cylinders. Comparison with an earlier inverse analysis of a similar experiment shows that the present method provides solutions that are fully consistent with the earlier results. Although the technique is illustrated using a simple one‐dimensional example, the method can be readily extended to multidimensional problems. Copyright © 2003 John Wiley & Sons, Ltd.     

Thermal Transport Properties of Layered Materials: Identification by a New Numerical Algorithm for Transient Measurements

Transient methods are widely used to determine thermal transport properties. In some situations they can be used for homogeneous media to measure several properties either simultaneously or separately. In this context an analytic model is available and a well-posed inverse problem of parameter identification has to be solved. The examination of composite media is more complicated. The algorithm proposed here allows simultaneous determination of the thermal conductivity and thermal diffusivity of layered dielectrics by transient measurements. It is based on a plane source that acts both as a resistive heater and temperature sensor. For the technique to be successful two essential aspects have to be considered: firstly, the mathematical modeling of the measured data (the forward problem) and secondly, the problem of ill-posedness of the inverse problem. For the proposed measurement configuration, a new fast data analysis algorithm based on an analytic solution for the forward problem is presented. In principle, a numerical solution such as an FEM solution of the heat conduction equation can be used instead of the analytical one, but the computational effort is much greater. The inverse problem is formulated as an output-least-squares problem, which leads to a transcendent algebraic system of equations. The method was successfully tested for different situations.Paper presented at the Fifteenth Symposium on Thermophysical Properties, June 22--27, 2003, Boulder, Colorado, U.S.A.     

The method of fundamental solutions for inverse source problems associated with the steady‐state heat conduction

This paper presents the use of the method of fundamental solutions (MFS) for recovering the heat source in steady‐state heat conduction problems from boundary temperature and heat flux measurements. It is well known that boundary data alone do not determine uniquely a general heat source and hence some a priori knowledge is assumed in order to guarantee the uniqueness of the solution. In the present study, the heat source is assumed to satisfy a second‐order partial differential equation on a physical basis, thereby transforming the problem into a fourth‐order partial differential equation, which can be conveniently solved using the MFS. Since the matrix arising from the MFS discretization is severely ill‐conditioned, a regularized solution is obtained by employing the truncated singular value decomposition, whilst the optimal regularization parameter is determined by the L‐curve criterion. Numerical results are presented for several two‐dimensional problems with both exact and noisy data. The sensitivity analysis with respect to two solution parameters, i.e. the number of source points and the distance between the fictitious and physical boundaries, and one problem parameter, i.e. the measure of the accessible part of the boundary, is also performed. The stability of the scheme with respect to the amount of noise added into the data is analysed. The numerical results obtained show that the proposed numerical algorithm is accurate, convergent, stable and computationally efficient for solving inverse source problems in steady‐state heat conduction. Copyright © 2006 John Wiley & Sons, Ltd.     

A meshless model for transient heat conduction in functionally graded materials

A meshless numerical model is developed for analyzing transient heat conduction in non-homogeneous functionally graded materials (FGM), which has a continuously functionally graded thermal conductivity parameter. First, the analog equation method is used to transform the original non-homogeneous problem into an equivalent homogeneous one at any given time so that a simpler fundamental solution can be employed to take the place of the one related to the original problem. Next, the approximate particular and homogeneous solutions are constructed using radial basis functions and virtual boundary collocation method, respectively. Finally, by enforcing satisfaction of the governing equation and boundary conditions at collocation points of the original problem, in which the time domain is discretized using the finite difference method, a linear algebraic system is obtained from which the unknown fictitious sources and interpolation coefficients can be determined. Further, the temperature at any point can be easily computed using the results of fictitious sources and interpolation coefficients. The accuracy of the proposed method is assessed through two numerical examples.     

梯度复合材料热传导分析的梯度单元法

给出了一种适用于梯度复合材料热传导分析的梯度单元, 采用细观力学方法描述材料变化的热物理属性, 通过线性插值和高阶插值温度场分别给出了4节点和8节点梯度单元随空间位置变化的热传导刚度矩阵。推导了在温度梯度载荷和热流密度载荷作用下, 矩形梯度板的稳态温度场和热通量场精确解。基于该精确解对比了连续梯度模型和传统的离散梯度模型的热传导有限元计算结果, 验证了梯度单元的有效性, 并讨论了相关参数对梯度单元的影响。结果表明, 梯度单元和均匀单元得到的温度场基本一致; 当热载荷垂直于材料梯度方向时, 梯度单元能够给出更加精确的局部热通量场; 当热载荷平行于材料梯度方向时, 4节点梯度单元性能恶化, 8节点梯度单元和均匀单元的计算结果与精确解吻合很好。     

Topology optimization of multi-material for the heat conduction problem based on the level set method

This article presents a numerical approach of topology optimization with multiple materials for the heat conduction problem. The multiphase level set model is used to implicitly describe the geometric boundaries of material regions with different conductivities. The model of multi-material representation has no emergence of the intermediate density. The optimization objective is to construct the optimal heat conductive paths which improve the efficiency of heat transfer. The dissipation of thermal transport potential capacity is taken as the objective function. The sensitivity analysis is implemented by the adjoint variable method, which is the foundation of constructing the velocity field of the level set equation. The optimal result is gradually realized by the evolution of multi-material boundaries, and the topological changes are naturally handled during the optimization process. Finally, the numerical examples are presented to demonstrate the feasibility and validity of the proposed method for topology optimization of the heat conduction problem.