Finite Difference-Peridynamic Differential Operator for Solving Transient Heat Conduction Problems

Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.

A Virtual Boundary Element Method for Three-Dimensional Inverse Heat Conduction Problems in Orthotropic Media

This paper aims to apply a virtual boundary element method(VBEM)to solve the inverse problems of three-dimensional heat conduction in orthotropic media.This method avoids the singular integrations in the conventional boundary element method,and can be treated as a potential approach for solving the inverse problems of the heat conduction owing to the boundary-only discretization and semi-analytical algorithm.When the VBEM is applied to the inverse problems,the numerical instability may occur if a virtual boundary is not properly chosen.The method encounters a highly illconditioned matrix for the larger distance between the physical boundary and the virtual boundary,and otherwise is hard to avoid the singularity of the source point.Thus,it must adopt an appropriate regularization method to deal with the ill-posed systems of inverse problems.In this study,the VBEM and different regularization techniques are combined to model the inverse problem of three-dimensional heat conduction in orthotropic media.The proper regularization techniques not only make the virtual boundary to be allocated freer,but also solve the ill-conditioned equation of the inverse problem.Numerical examples demonstrate that the proposed method is efficient,accurate and numerically stable for solving the inverse problems of three-dimensional heat conduction in orthotropic media.

PRECISE INTEGRAL ALGORITHM BASED SOLUTION FORTRANSIENT INVERSE HEAT CONDUCTION PROBLEMSWITH MULTI-VARIABLES

By modeling direct transient heat conduction problems via finite element method (FEM) and precise integral algorithm, a new approach is presented to solve transient inverse heat conduction problems with multi-variables. Firstly, the spatial space and temporal domain are discretized by FEM and precise integral algorithm respectively. Then, the high accuracy semi-analytical solution of direct problem can be got. Finally, based on the solution, the computing model of inverse problem and expression of sensitivity analysis are established. Single variable and variables combined identifications including thermal parameters, boundary conditions and source-related terms etc. are given to validate the approach proposed in 1-D and 2-D cases. The effects of noise data and initial guess on the results are investigated. The numerical examples show the effectiveness of this approach.

Newton-type methods and their modifications for inverse heat conduction problems

This paper studies to numerical solutions of an inverse heat conduction problem.The effect of algorithms based on the Newton-Tikhonov method and the Newton-implicit iterative method is investigated,and then several modifications are presented.Numerical examples show the modified algorithms always work and can greatly reduce the computational costs.

Improved Particle Swarm Optimization for Solving Transient Nonlinear Inverse Heat Conduction Problem in Complex Structure

Accurately solving transient nonlinear inverse heat conduction problems in complex structures is of great importance to provide key parameters for modeling coupled heat transfer process and the structure’s optimization design.The finite element method in ABAQUS is employed to solve the direct transient nonlinear heat conduction problem.Improved particle swarm optimization(PSO)method is developed and used to solve the transient nonlinear inverse problem.To investigate the inverse performances,some numerical tests are provided.Boundary conditions at inaccessible surfaces of a scramjet combustor with the regenerative cooling system are inversely identified.The results show that the new methodology can accurately and efficiently determine the boundary conditions in the scramjet combustor with the regenerative cooling system.By solving the transient nonlinear inverse problem,the improved particle swarm optimization for solving the transient nonlinear inverse heat conduction problem in a complex structure is verified.

NUMERICAL METHOD OF MIXED FINITE VOLUME-MODIFIED UPWIND FRACTIONAL STEP DIFFERENCE FOR THREE-DIMENSIONAL SEMICONDUCTOR DEVICE TRANSIENT BEHAVIOR PROBLEMS

Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.

A new complex variable meshless method for transient heat conduction problems

In this paper, based on the improved complex variable moving least-square （ICVMLS） approximation, a new complex variable meshless method （CVMM） for two-dimensional （2D） transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.

A complex variable meshless local Petrov-Galerkin method for transient heat conduction problems

On the basis of the complex variable moving least-square （CVMLS） approximation, a complex variable meshless local Petrov-Galerkin （CVMLPG） method is presented for transient heat conduction problems. The method is developed based on the CVMLS approximation for constructing shape functions at scattered points, and the Heaviside step function is used as a test function in each sub-domain to avoid the need for a domain integral in symmetric weak form. In the construction of the well-performed shape function, the trial function of a two-dimensional （2D） problem is formed with a one-dimensional （1D） basis function, thus improving computational efficiency. The numerical results are compared with the exact solutions of the problems and the finite element method （FEM）. This comparison illustrates the accuracy as well as the capability of the CVMLPG method.

Resolving double-sided inverse heat conduction problem using calibration integral equation method

In this paper, a novel calibration integral equation is derived for resolving double-sided, two-probe inverse heat conduction problem of surface heat flux estimation. In contrast to the conventional inverse heat conduction techniques, this calibration approach does not require explicit input of the probe locations, thermophysical properties of the host material and temperature sensor parameters related to thermal contact resistance, sensor capacitance and conductive lead losses. All those parameters and properties are inherently contained in the calibration framework in terms of Volterra integral equation of the first kind. The Laplace transform technique is applied and the frequency domain manipulations of the heat equation are performed for deriving the calibration integral equation. Due to the ill-posed nature, regularization is required for the inverse heat conduction problem, a future-time method or singular value decomposition (SVD) can be used for stabilizing the ill-posed Volterra integral equation of the first kind.

Inverse Problem on Heat Conduction in Heterogeneous Medium

Under consideration is a nonclassical stationary problem on heat conduction in a body with the pre-set surface temperature and heat flow. The body contains inclusions at unknown locations and with unknown boundaries. The body and inclusions have different constant thermal conductivities. The author explores the possibility of locating inclusions. The article presents an integral criterion based on which a few statements on identification of inclusions in a body are proved.

Three-dimensional numerical manifold method for heat conduction problems with a simplex integral on the boundary

The three-dimensional numerical manifold method(3D-NMM),which is based on the derivation of Galerkin's variation,is a powerful calculation tool that uses two cover systems.The 3D-NMM can be used to handle continue-discontinue problems and extend to THM coupling.In this study,we extended the 3D-NMM to simulate both steady-state and transient heat conduction problems.The modelling was carried out using the raster methods(RSM).For the system equation,a variational method was employed to drive the discrete equations,and the crucial boundary conditions were solved using the penalty method.To solve the boundary integral problem,the face integral of scalar fields and two-dimensional simplex integration were used to accurately describe the integral on polygonal boundaries.Several numerical examples were used to verify the results of 3D steady-state and transient heat-conduction problems.The numerical results indicated that the 3D-NMM is effective for handling 3D both steadystate and transient heat conduction problems with high solution accuracy.

Statistical Multiscale Analysis of Transient Conduction and Radiation Heat Transfer Problem in Random Inhomogeneous Porous Materials

This paper is devoted to the homogenization and statistical multiscale analysis of a transient heat conduction problem in random porous materials with a nonlinear radiation boundary condition.A novel statistical multiscale analysis method based on the two-scale asymptotic expansion is proposed.In the statistical multiscale formulations,a unified linear homogenization procedure is established and the second-order correctors are introduced for modeling the nonlinear radiative heat transfer in random perforations,which are our main contributions.Besides,a numerical algorithm based on the statistical multiscale method is given in details.Numerical results prove the accuracy and efficiency of our method for multiscale simulation of transient nonlinear conduction and radiation heat transfer problem in random porous materials.

The Solving of the Inverse Thermal Conductivity Problem for Study the Short Linear Heat Pipes

The results of studies by solving the inverse thermal conductivity problem of the heat capacity of evaporator of the short linear heat pipes (HP’s) with a Laval nozzle-liked vapour channel and intended for cooling spacecraft and satellites with strict take-off mass regulation are presented. Mathematical formulation of the inverse problem for the HP’s thermal conductivity in one-dimensional coordinate system is accompanied by the measurement results using the monotonic heating method in a vacuum adiabatic calorimeter the HP’s surface temperatures along the longitudinal axis over the entire temperature load range, thermal resistance, and arrays of thermal power data on the evaporator Qev and vortex flow calorimeter Qcond for the condensation surface allow us to estimate the average value of the evaporator heat capacity Cev by solving the inverse thermal conductivity problem in the HP’s evaporator region. Since at the beginning of working fluid boiling for a certain time interval, the temperature of the capillary-porous evaporator remains close to constant, and with the continuation of heating and by solving the inverse thermal conductivity problem, it becomes possible to calculate the heat capacity of the working evaporator and the evaporation specific heat of the boiling working fluid and compare it with the table values.

A meshless model for transient heat conduction analyses of 3D axisymmetric functionally graded solids

A meshless numerical model is developed for analyzing transient heat conductions in three-dimensional （3D） axisymmetric continuously nonhomogeneous functionally graded materials （FGMs）. Axial symmetry of geometry and boundary conditions reduces the original 3D initial-boundary value problem into a two-dimensional （2D） problem. Local weak forms are derived for small polygonal sub-domains which surround nodal points distributed over the cross section. In order to simplify the treatment of the essential boundary conditions, spatial variations of the temperature and heat flux at discrete time instants are interpolated by the natural neighbor interpolation. Moreover, the using of three-node triangular finite element method （FEM） shape functions as test functions reduces the orders of integrands involved in domain integrals. The semi-discrete heat conduction equation is solved numerically with the traditional two-point difference technique in the time domain. Two numerical examples are investigated and excellent results are obtained, demonstrating the potential application of the proposed approach.

A Numerical Method on Inverse Determination of Heat Transfer Coefficient Based on Thermographic Temperature Measurement

The heat transfer coefficient in a multidimensional heat conduction problem is obtained from the solution of the inverse heat conduction problem based on the thermographic temperature measurement. The modified one-dimensional correction method （MODCM）, along with the finite volume method, is employed for both two- and three-dimensional inverse problems. A series of numerical experiments are conducted in order to verify the effectiveness of the method. In addition, the effect of the temperature measurement error, the ending criterion of the iteration, etc. on the result of the inverse problem is investigated. It is proved that the method is a simple, stable and accurate one that can solve successfully the inverse heat conduction problem.

Hlder Stability Estimate for an Inverse Parabolic Problem

This paper deals with a parabolic system in a multi dimentional bounded domain ΩR n with the smooth boundary Ω. We discuss an inverse parabolic problem of determining the indirectly measurable internal heat distribution at any intermediate moment from the heat distribution measurements in arbitrary accessible subdomain ωΩ at some time interval. Our main result is the Hlder stability estimate in the inverse problem and the proof is completed with a Carleman estimate and a eigenfunction expansion for parabolic equations.

A dimension-reduced neural network-assisted approximate Bayesian computation for inverse heat conduction problems

Due to the flexibility and feasibility of addressing ill-posed problems,the Bayesian method has been widely used in inverse heat conduction problems(IHCPs).However,in the real science and engineering IHCPs,the likelihood function of the Bayesian method is commonly computationally expensive or analytically unavailable.In this study,in order to circumvent this intractable likelihood function,the approximate Bayesian computation(ABC)is expanded to the IHCPs.In ABC,the high dimensional observations in the intractable likelihood function are equalized by their low dimensional summary statistics.Thus,the performance of the ABC depends on the selection of summary statistics.In this study,a machine learning-based ABC(ML-ABC)is proposed to address the complicated selections of the summary statistics.The Auto-Encoder(AE)is a powerful Machine Learning(ML)framework which can compress the observations into very low dimensional summary statistics with little information loss.In addition,in order to accelerate the calculation of the proposed framework,another neural network(NN)is utilized to construct the mapping between the unknowns and the summary statistics.With this mapping,given arbitrary unknowns,the summary statistics can be obtained efficiently without solving the time-consuming forward problem with numerical method.Furthermore,an adaptive nested sampling method(ANSM)is developed to further improve the efficiency of sampling.The performance of the proposed method is demonstrated with two IHCP cases.

The dimension splitting element-free Galerkin method for 3D transient heat conduction problems

By transforming a 3D problem into some related 2D problems, the dimension splitting element-free Galerkin(DSEFG) method is proposed to solve 3D transient heat conduction problems. The improved element-free Galerkin(IEFG) method is used for 2D transient heat conduction problems, and the finite difference method is applied in the splitting direction. The discretized system equation is obtained based on the Galerkin weak form of 2D problem; the essential boundary conditions are imposed with the penalty method; and the finite difference method is employed in the time domain. Four exemplary problems are chosen to verify the efficiency of the DSEFG method. The numerical solutions show that the efficiency and precision of the DSEFG method are greater than ones of the IEFG method for 3D problems.

Solution of Inverse Heat Conduction Problem using the Tikhonov Regulariza-tion Method

It is hard to solve ill-posed problems, as calculated temperatures are very sensitive to errors made while calculating "measured" temperatures or performing real-time measurements. The errors can create temperature oscillation, which can be the cause of an unstable solution. In order to overcome such difficulties, a variety of techniques have been proposed in literature, including regularization, future time steps and smoothing digital filters. In this paper, the Tikhonov regularization is applied to stabilize the solution of the inverse heat conduction problem. The impact on the inverse solution stability and accuracy is demonstrated.

A mixed Newton-Tikhonov method for nonlinear ill-posed problems

Newton type methods are one kind of the efficient methods to solve nonlinear ill-posed problems, which have attracted extensive attention. However, computational cost of Newton type methods is high because practical problems are complicated. We propose a mixed Newton-Tikhonov method, i.e., one step Newton-Tikhonov method with several other steps of simplified Newton-Tikhonov method. Convergence and stability of this method are proved under some conditions. Numerical experiments show that the proposed method has obvious advantages over the classical Newton method in terms of computational costs.