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.

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

处于一个多维的热传导问题的热转移系数基于 thermographic 温度测量从反的热传导问题的答案被获得。与有限体积方法一起，修改一个维的修正方法(MODCM ) 为两 twoand 被采用三维的反的问题。一系列数字实验被进行以便验证方法的有效性。另外，温度测量错误的效果，重复的结束标准，反的问题的结果上的等等被调查。方法是能成功地解决反的热传导问题的简单、稳定、精确的，这被证明。

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.

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

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

Fourier正则化同时反演热传导方程的初值和源项

研究边值问题中从已知的给定数据中同时识别初始条件和未知源的问题。用Fourier变换求出问题的精确解,发现该反问题是严重不适定的。为了解决所提出的问题,本文中用Fourier正则化方法构造出正则化解近似解,并在边界条件的假设下,得到先验参数选取下的精确解与其正则化近似解之间的对数型误差估计。

基于局部观测数据的一维热传导方程系数与初值同时反演问题的数值计算

利用局部观测数据研究了一维热传导方程系数与初值同时反演问题。首先利用函数代换及热传导方程解的解析延拓性,证明了基于局部观测数据下系数与初值同时反演问题的唯一性;然后基于深度神经网络,构造了同时反演问题的反演算法;最后通过构造具有解析解和不具有解析解的数值算例,验证了反演算法的有效性。

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.

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.

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.

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

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

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.

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 complementary variational principles with a variable domain in terms of temperature and heat-streamfunction are flrst established.Based on them,two methods of solution-generalized Ritz method and variable-domain FEM-both capable of handiling 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.

基于稳态传热分析的高炉炉壁温度场建模方法

高炉炉壁温度场是掌握炉内热负荷状态和燃料燃烧情况、及时发现高炉异常工况、保障高炉稳定顺行的重要信息。针对现有炉壁温度场获取难及温度场建模过程中未考虑渣皮脱落与生成现象等问题,提出一种基于稳态传热分析的高炉炉壁温度场建模方法。首先,利用COMSOL仿真软件搭建炉壁三维仿真模型;其次,充分考虑渣皮对炉壁温度变化的影响,建立包含渣皮热源项的稳态传热模型,获得初始炉壁内表面温度场;然后,为提高温度场模型精度,基于热传导反问题理论,以冷却壁中热电偶的实测温度为参考,构建炉壁温度场修正模型,实现炉壁温度场数据的迭代修正。实验表明经温度场修正模型修正后,传热模型的计算误差控制在5%以内,满足现场检测精度要求,并能可视化展示炉壁温度场,为高炉调控操作提供了可靠的炉壁温度场数据。

二维时间反向热传导问题的正则化方法及误差估计

讨论二维时间反向热传导问题,从终值时刻t=T(T>0)的温度分布来反演初始时刻的温度分布.该问题在图像处理方面有重要应用.这是一个严重不适定问题,它的解在一定条件下不连续依赖于数据.针对传统正则化方法的缺陷,采用拟逆正则化方法和分数次Tikhonov正则化方法恢复解对数据的依赖性.同时,还给出2种方法相应的先验参数选取规则及其正则解与精确解的误差估计.