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.

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.

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.

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.

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.

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 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.

Regularization Method to the Parameter Identification of Interfacial Heat Transfer Coefficient and Properties during Casting Solidification

The accurate material physical properties, initial and boundary conditions are indispensable to the numerical simulation in the casting process, and they are related to the simulation accuracy directly. The inverse heat conduction method can be used to identify the mentioned above parameters based on the temperature measurement data. This paper presented a new inverse method according to Tikhonov regularization theory. A regularization functional was established and the regularization parameter was deduced, the Newton-Raphson iteration method was used to solve the equations. One detailed case was solved to identify the thermal conductivity and specific heat of sand mold and interfacial heat transfer coefficient (IHTC) at the meantime. This indicates that the regularization method is very efficient in decreasing the sensitivity to the temperature measurement data, overcoming the ill-posedness of the inverse heat conduction problem (IHCP) and improving the stability and accuracy of the results. As a general inverse method, it can be used to identify not only the material physical properties but also the initial and boundary conditions' parameters.

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.

Spatial interfacial heat transfer and surface characteristics during gravity casting of A356 alloy

As one of the key boundary conditions during casting solidification process, the interfacial heat transfer coefficient (IHTC) affects the temperature variation and distribution. Based on the improved nonlinear estimation method (NEM), thermal measurements near both bottom and lateral metal-mold interfaces throughout A356 gravity casting process were carried out and applied to solving the inverse heat conduction problem (IHCP). Finite element method (FEM) is employed for modeling transient thermal fields implementing a developed NEM interface program to quantify transient IHTCs. It is found that IHTCs at the lateral interface become stable after the volumetric shrinkage of casting while those of the bottom interface reach the steady period once a surface layer has solidified. The stable value of bottom IHTCs is 750 W/(m^2·℃), which is approximately 3 times that at the lateral interface. Further analysis of the interplay between spatial IHTCs and observed surface morphology reveals that spatial heat transfer across casting-mold interfaces is the direct result of different interface evolution during solidification process.

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.

Numerical Methods for Inverse Problem of Heat Conduction with Unknown Boundary Based on Variational Principles with Variable Domain

In this article a variable-domain variational approach to the entitled problem is presented.A pair of comple- mentary variational principles with a variable domain in terms of temperature and heat-streamfunction are first established.Based on them,two methods of solution—generalized Ritz method and variable-domain FEM— both capable of handling problems with unknown boundaries,are suggested.Then,three sample numerical examples have been tested.The computational process is quite stable,and the results are encouraging.This variational approach can be extended straightforwardly to 3-D inverse problems as well as to other problems in mathematical physics.

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

By introducing an image plane, the inverse heat conduction problem with free boundary is transformed into one with completely known boundaryt which is much simpler to handle.As a by-product, the classical Kirchhoff’s transformation for accounting for variable conductivity is rederived and an invariance property of the inverse problem solution with respect to variable conductivity is indicated. Then a pair of complementary extremum principles are established on the image plane, providing a sound theoretical foundation for the Ritz’s method and finite element method (FEM).An example solved by FEM is also given.

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.

FOURIER REGULARIZATION FOR DETERMINING SURFACE HEAT FLUX FROM INTERIOR OBSERVATION BASED ON A SIDEWAYS PARABOLIC EQUATION

In this paper we consider a non-standard inverse heat conduction problem for determining surface heat flux from an interior observation which appears in some applied subjects. This problem is ill-posed in the sense that the solution (if it exists)does not depend continuously on the data. A Fourier method is applied to formulate a regularized approximation solution, and some sharp error estimates are also given.

热传导方程Robin系数反问题解的唯一性及正则化解的存在性

Robin系数在热传导模型中刻画了热传导区域边界上的热交换,是一类非常重要的参数,本文基于某小时段温度测量值反演热传导模型中的Robin系数.首先,在边界值以及测量值满足一定的光滑性条件时,给出了反问题解的唯一性;其次,基于Tikhonov正则化思想,通过构造目标泛函将反问题转化为求目标泛函的极小值,并证明了泛函极小元的存在性.

反求铝合金淬火边界条件的序列LM梯度算法

准确获取边界条件在淬火仿真中至关重要。为此,根据瞬态热传导模型,利用正则化方法建立反求热通量的目标函数,结合序列函数指定法与LM梯度法,通过最小化目标函数推导出序列LM梯度算法,实现末端淬火实验中热通量的迭代计算。其次,通过训练观测点的平均绝对误差的“3-7-3”神经网络预测模型,建立了序列LM梯度算法关键参数的遗传算法优化方法。最后,以换热系数为边界条件建立淬火仿真模型。结果表明:在相同淬火时间内,一维序列LM梯度算法、一点法和两点法,在观测点的仿真温度平均相对误差分别为5.6%、5.9%、12.6%,17.1%、15.4%、15.3%和17.7%、16.0%、15.6%。因此,与一点法、两点法相比,序列LM梯度算法具有更高的精度。