A Precise Integration Boundary-Element Method for Solving Transient Heat Conduction Problems with Variable Thermal Conductivity

In this article, a combined approach of the radial integration boundary element method (RIBEM) and the precise integration method is presented for solving transient heat conduction problems with variable thermal conductivity. First, the system of ordinary differential equations on the boundary integral equation can be obtained by the RIBEM. Then, the precise integration method is adopted to solve the system of ordinary differential equations. Finally, three numerical examples are presented to demonstrate the performance of the present method. The results show that the present approach can obtain satisfactory performance even for very large time-step size.     

A Precise Time-Domain Expanding Boundary-Element Method for Solving Three-Dimensional Transient Heat Conduction Problems with Variable Thermal Conductivity

A combined approach of the radial integration boundary-element method (RIBEM) and the precise algorithm in the time domain is presented for solving three-dimensional transient heat conduction problems with variable thermal conductivity. First, by expanding physical quantities at discrete time intervals, the recursive formulation of the governing equation is derived. Then, the recursive equation is solved by the RIBEM, and a self-adaptive check technique is carried out to estimate how many expansion terms are needed in a time step. Finally, three numerical examples show that the present approach can obtain very stable and accurate results for different time-step size.     

A Combined Approach of RIBEM and Precise Time Integration Algorithm for Solving Transient Heat Conduction Problems

A combined approach of the radial integration boundary-element method (RIBEM) and the precise algorithm in the time domain is presented for solving transient heat conduction problems with heat sources. First, by expanding physical quantities at discrete time intervals, the recursive formulations of the governing equation are derived. Then, the recursive equation is solved by the RIBEM. Finally, a self-adaptive check technique is carried out to estimate how many expansion terms are needed in a time-step size. Numerical results show that the present approach can obtain satisfactory performance.     

Radial integration boundary element method for heat conduction problems with convective heat transfer boundary

A new boundary domain integral equation with convective heat transfer boundary is presented to solve variable coefficient heat conduction problems. Green’s function for the Laplace equation is used to derive the basic integral equation with varying heat conductivities, and as a result, domain integrals are included in the derived integral equations. The existing domain integral is converted into an equivalent boundary integral using the radial integration method by expressing the normalized temperature as a series of radial basis functions. This treatment results in a pure boundary element analysis algorithm and requires no internal cells to evaluate the domain integral. Numerical examples are presented to demonstrate the accuracy and efficiency of the present method.     

The BEM for point heat source estimation: application to multiple static sources and moving sources

This paper deals with an inverse problem, which consists of the identification of point heat sources in a homogeneous solid in transient heat conduction. The location and strength of the line heat sources are both unknown. For a single source we examine the case of a source which moves in the system during the experiment. The two-dimensional and three-dimensional linear heat conduction problems are considered here. The identification procedure is based on a boundary integral formulation using transient fundamental solutions. The discretized problem is non-linear if the location of the line heat sources is unknown. In order to solve the problem we use an iterative procedure to minimize a quadratic norm. The proposed numerical approach is applied to experimental 2D examples using measurements provided by an infrared scanner for surface temperatures and heat fluxes. A numerical example is presented for the 3D application.     

Relativistic moving heat source

The relativistic heat conduction (RHCE) model is particularly important in the analysis of processes involving moving heat sources (MHS) at speeds or frequencies comparable with those of heat propagation in the medium. This paper establishes a unified framework for solving heat conduction problems using the RHCE model. It offers “Fundamental Solutions” in one, two, and three spatial dimensions, for the transient response due to an instantaneous point MHS. Moreover, it presents the transient response due to a continuous point MHS, the quasi-steady response due a periodic point MHS, as well as guidelines for solving the RHCE equation under various loading and boundary conditions.     

Precise Time-Domain Expanding BEM for Solving Non-Fourier Heat Conduction Problems

The precise time-domain expanding boundary-element method (BEM) is presented for solving non-Fourier heat conduction problems. The recursive boundary integral equation is obtained via the precise time-domain expanding method and solved by the BEM, where the radial integral method is used to transform the domain integral into the boundary integral. Also, a self-adaptive judging criterion is used in the solving process. The transformation matrices of domain integrals need to be computed only once, except those related to the heat source. Finally, numerical results show that the present method can obtain stable and accurate results with different time steps.     

A new boundary meshfree method for calculating the multi-domain constant coefficient heat conduction with a heat source problem

A new high-precision boundary meshfree method, namely virtual boundary meshfree Galerkin method (VBMGM), for calculating the multi-domain constant coefficient heat conduction with a heat source problem is given. In the paper, the radial basis function interpolation is used to solve the virtual source function of virtual boundary and the heat source within each subdomain. Simultaneously, the equation of VBMGM for multi-domain constant coefficient heat conduction with a heat source problem is obtained by the Galerkin method. Therefore, the proposed method has common advantages of the boundary element method, meshfree method, and Galerkin method. Coefficient matrix of this specific expression is symmetrical and the specific expression of VBMGM for the multi-domain constant coefficient heat conduction with a heat source problem is given. Two numerical examples are given. The numerical results are also compared with other numerical methods. The accuracy and feasibility of the method for the multi-domain constant coefficient heat conduction with a heat source problem are proved.     

Non-Fourier effects on transient temperature resulting from periodic on-off heat flux

The transient temperatures resulting from a periodic on-off heat flux boundary condition have many applications, including, among others, the sintering of catalysts frequently found during coke burnoff, and the use of laser pulses for annealing of semiconductors. In such situations, the duration of the pulses is so small (i.e. picosecond-nanosecond) that the classical heat diffusion phenomenon breaks down and the wave nature of energy propagation characterized by the hyperbolic heat conduction equation governs the temperature distribution in the medium. In this work, an explicit analytic solution is presented for a linear transient heat conduction problem in a semi-infinite medium subjected to a periodic on-off type heat flux at the boundary x = 0 by solving the hyperbolic heat conduction equation. The non-linear case allowing for the added effect of surface radiation into an external ambient is studied numerically.     

A regularized Fourier sine series solution of a 3D backward heat conduction problem with extremal long time span

Abstract

In the article, we solve an extremal long time span backward heat conduction problem (BHCP) of a 3D nonhomogeneous heat conduction equation with nonhomogeneous boundary conditions in a cuboid. We first derive a time-dependent 3D homogenization function, such that in terms of the new variable by a variable transformation we can find the expansion coefficients in a closed form by using the Fourier sine series method. After a simple regularization technique, a stable analytic series solution of the 3D BHCP is available. We also develop a regularized Fourier sine series solution of temperature in the whole space-time domain. Numerical tests for the BHCPs in a large space-time domain reveal that the present method is very accurate to recover the initial temperature and the whole solution.     

A heat-flux based “building block” approach for solving heat conduction problems

A numerical approximation of the Green’s function equation based on a heat-flux formulation is given. It is derived by assuming as a functional form of the surface heat flux a stepwise variation with space and time. The obtained approximation is very important in investigation of the inverse heat conduction problems (IHCPs) because it gives a convenient expression for the temperature in terms of the heat flux components. Additionally, it is very important for the unsteady surface element (USE) method which is a modern boundary discretization method. Green’s function approximate solution equation (GFASE) also creates ‘naturally’ fixed groups or modules of work elements called “building blocks” that may be added together to obtain space and time values of temperature. In the current case, they are subject to a partial heating by an applied surface heat flux. The “building block” solution can be derived by using the various analytical and numerical approaches available in heat conduction literature though the exact analysis is preferable, as discussed in the text. Poorly-convergent series deriving from Green’s functions approach are replaced by closed-form algebraic solutions.     

An efficient BEM formulation for three-dimensional steady-state heat conduction analysis of composites

In this work, an efficient boundary element formulation has been presented for three-dimensional steady-state heat conduction analysis of fiber reinforced composites. The cylindrical shaped fibers in the three-dimensional composite matrix are represented by a system of curvilinear line elements with a prescribed diameter which facilitates efficient analysis and modeling together with the reduction in dimensionality of the problem. The variations in the temperature and flux fields in the circumferential direction of the fiber are represented in terms of a trigonometric shape function together with a linear or quadratic variation in the longitudinal direction. The resulting integrals are then treated semi-analytically which reduces the computational task significantly. The computational effort is further minimized by analytically substituting the fiber equations into the boundary integral equation of the material matrix with hole, resulting in a modified boundary integral equation of the composite matrix. An efficient assembly process of the resulting system equations is demonstrated together with several numerical examples to validate the proposed formulation. An example of application is also included.     

On the minimum entropy production in steady state heat conduction processes

On the basis of minimum entropy generation principle, a new formulation of the boundary value problems is proposed. Applying Euler–Lagrange variational formalism, a new mathematical form of heat conduction equation with additional heat source terms has been derived. To obtain a unique solution a special mathematical form of boundary conditions for 2D and 3D problems is required. As a result, entropy generation rate of the process can significantly be reduced, which leads to the decrease of the irreversibility ratio according to the Gouy–Stodola theorem. Minimization of entropy generation in heat conduction process is always possible by introducing additional heat sources.     

Multiple transient point heat sources identification in heat diffusion: application to experimental 2D problems

This paper deals with an inverse problem, which consists of the experimental identification of line heat sources in a homogeneous solid in transient heat conduction. The location and strength of the line heat sources are both unknown. For a single source we examine the case of a source which moves in the system during the experiment. The identification procedure is based on a boundary integral formulation using transient fundamental solutions. The discretized problem is non-linear if the location of the line heat sources is unknown. In order to solve the problem we use an iterative procedure to minimize a cost function comparing the modelled heat source term and the measurements. The proposed numerical approach is applied to experimental 2D examples using measurements provided by an infrared scanner for surface temperatures and heat fluxes. In some particular examples, internal thermocouples can be used. A time regularization procedure associated to future time-steps is used to correctly solve the ill-posed problem.     

Application of the method of fundamental solutions and radial basis functions for inverse heat source problem in case of steady-state

The paper deals with the inverse determination of heat sources in steady 2-D heat conduction problem. The problem is described by Poisson equation in which the function of the right hand side is unknown. The identification of the strength of a heat source is given by using the boundary condition and a known value of temperature in chosen points placed inside the domain. For the solution of the inverse problem of identification of the heat source the method of fundamental solution with radial basis functions is proposed. The accurate results have been obtained for five test problems where the analytical solutions were available.     

A semi-analytical boundary collocation solver for the inverse Cauchy problems in heat conduction under 3D FGMs with heat source

Abstract

In this article, the inverse Cauchy problems in heat conduction under 3D functionally graded materials (FGMs) with heat source are solved by using a semi-analytical boundary collocation solver. In the present semi-analytical solver, the combined boundary particle method and regularization technique is employed to deal with ill-pose inverse Cauchy problems. The domain mapping method and variable transformation are introduced to derive the high-order general solutions satisfying the heat conduction equation of 3D FGMs. Thanks to these derived high-order general solutions, the proposed scheme can only require the boundary discretization to recover the solutions of the heat conduction equations with a heat source. The regularization technique is used to eliminate the effect of the noisy measurement data on the accessible boundary surface of 3D FGMs. The efficiency of the proposed solver for inverse Cauchy problems is verified under several typical benchmark examples related to 3D FGM with specific spatial variations (quadratic, exponential and trigonometric functions).     

A hybrid Green’s function method for the hyperbolic heat conduction problems

The present study is devoted to propose a hybrid Green’s function method to investigate the hyperbolic heat conduction problems. The difficulty of the numerical solutions of hyperbolic heat conduction problems is the numerical oscillation in the vicinity of sharp discontinuities. In the present study, we have developed a hybrid method combined the Laplace transform, Green’s function and ε-algorithm acceleration method for solving time dependent hyperbolic heat conduction equation. From one- to three-dimensional problems, six different examples have been analyzed by the present method. It is found from these examples that the present method is in agreement with the Tsai-tse Kao’s solutions [Tsai-tse Kao, Non-Fourier heat conduction in thin surface layers, J. Heat Transfer 99 (1977) 343–345] and does not exhibit numerical oscillations at the wave front. The propagation of the two- and three-dimensional thermal wave becomes so complicated because it occur jump discontinuities, reflections and interactions in these numerical results of the problem and it is difficult to find the analytical solutions or the result of other study to compare with the solutions of the present method.     

Fractional heat conduction in a space with a source varying harmonically in time and associated thermal stresses

Time-nonlocal generalization of the classical Fourier law with the “long-tail” power kernel can be interpreted in terms of fractional calculus (theory of integrals and derivatives of noninteger order) and leads to the time-fractional heat conduction equation with the Caputo derivative. Fractional heat conduction equation with the harmonic source term under zero initial conditions is studied. Different formulations of the problem for the standard parabolic heat conduction equation and for the hyperbolic wave equation appearing in thermoelasticity without energy dissipation are discussed. The integral transform technique is used. The corresponding thermal stresses are found using the displacement potential.     

Transient analytical solution to heat conduction in composite circular cylinder

An analytical method leading to the solution of transient temperature filed in multi-dimensional composite circular cylinder is presented. The boundary condition is described as time-dependent temperature change. For such heat conduction problem, nearly all the published works need numerical schemes in computing eigenvalues or residues. In this paper, the proposed method involves no such numerical work. Application of ‘separation of variables’ is novel. The developed method represents an extension of the analytical approach derived for solving heat conduction in composite slab in Cartesian coordinates. Close-formed solution is provided and its agreement with numerical result is good which demonstrates a good accuracy of the developed solution form.