A boundary identification method for an inverse heat conduction problem with an application in ironmaking

A boundary identification problem in inverse heat conduction is studied, based on data from internal measurement of temperature and heat flux. Formulated as a sideways heat conduction equation, a spatial continuation technique is applied to extend the solution to a known boundary condition at the desired boundary position. Recording the positions traversed in the continuation for each time instant yields the boundary position trajectory and hence the solution of the identification problem. A prospective application of the method can be found in the ironmaking blast furnace, where it is desired to monitor the thickness of the accreted refractory wall based on measurement of its internal state. Simulations featuring noisy measurement data demonstrate the feasibility of the identification method for blast furnace wall thickness estimation.     

共轭梯度法求解非线性多宗量稳态传热反问题

应用共轭梯度法求解非线性多宗量稳态热传导反问题。采用八节点的等参单元在空间上进行离散,建立了便于敏度分析的非线性正演和反演的有限元模型,可直接求导进行敏度分析。给出了相关的数值验证,对测量误差及测点数目的影响作了初步探讨,结果表明,采用的算法能够对非线性稳态热传导中导热系数和边界条件联合反问题进行有效的求解,并具有较高精度。     

Determination of nonuniform heat flux density on boundary in inverse heat conduction problem

This paper develops reducing deviation method to solve inverse heat conduction problem with nonuniform heat flux density on boundary. When the thermocouple is placed a bit far from the boundary, there is no problem on non-convergence in this method. A calculation example is given.     

Whole time domain solution of inverse heat conduction problem in multi-layer media

The unknown surface conditions in composite media is estimated by minimizing the nonlinear least squares error between the computed and measured temperatures over the whole time domain. This approach shown to be stable, efficient and accurate. The unknown surface conditions are assumed to have an abrupt change at unknown time. A sensitivity analysis is conducted to learn more insight into the nature of difficulties that can be encountered in the estimation of the parameters associated with the inverse problem. The stability and accuracy of the method is demonstrated by several numerical examples which provide very strict test conditions. Received on 7 March 1997     

A steady conjugate heat transfer problem with conduction and free convection

Summary A steady conjugate heat transfer problem dealing with conduction in a heat-generating slab and free convection in the surrounding fluid is studied analytically. Free convection is analyzed by a Görtler-type series solution to the boundary-layer equations for non-uniform surface-temperature variations, while conduction is treated by the standard technique of Fourier transforms. Interfacial temperature and heat flux variations from both solutions in series forms are then formally matched to yield algebraic relations for the coefficients in the series. These coefficients can then be simply evaluated in a given problem in terms of three physical parameters. A numerical example is shown.     

An inverse geometry problem in estimating frost growth on an evaporating tube

 When humid air comes into contact with a surface whose temperature is below the dew point of water vapor in air and also below the freezing point, frost deposition takes place over the surface. The phenomena of the frost growth are very complicated and therefore it is very difficult to model mathematically the behavior of frost growth and predict it. In the present study a transient inverse geometry heat conduction problem (shape identification problem) is solved using the conjugate gradient method (CGM) and boundary element method (BEM)-based inverse algorithm to estimate the unknown irregular frost thickness and shape. Results obtained by using the CGM to estimate the frost growth are justified based on the numerical experiments. It is concluded that the accurate frost shape can be estimated by the CGM except for the initial and final time. The reason and improvement of this singularity are addressed. Finally the effects of reducing the number of sensors and increasing the measurement errors on the inverse solutions are discussed. Received on 25 September 2000 / Published online: 29 November 2001     

变几何域传热的表面热流反演方法

变几何域的表面热流反演是一类特殊的热传导逆问题,在再入飞行器烧蚀型防热材料的表面热流反演中具有工程实用价值.本文首先对变几何域传热的正问题计算方法进行了校核验证,然后建立了求解变几何域表面热流反演问题的顺序函数法和共轭梯度法;给出了这两种反演方法的基本思想和算法推导,并针对典型算例进行了仿真.结果表明:两种反演方法都能计算出较好的反演结果,并且算法受测量噪声的影响较小,具有较好的鲁棒性;反演算法能适应不同的几何域变化函数,但几何域变化量的测量误差在表面热流的反演结果中会有较为直接的反映.     

Analytical solution of the overdetermined inverse heat conduction problem with an application to monitoring thermal stresses

The paper presents a technique for determining the transient temperature distribution, when input data as thermocouple responses are known at several interior locations. If the temperature field is known, then the thermal stresses can be calculated. The problem is overdetermined and is solved using a least squares method that minimizes the error between the computed and measured thermocouple temperatures. The present method incorporates the advantages of simplicity and accuracy of analytical solutions. Several numerical examples and measurements are presented as an indication of the accuracy of the presented method. Received on 10 January 1997     

水坝瞬态热传导分析中的拟初始条件边界元法

使用一种时域边界元方法对混凝土水坝进行瞬态热传导分析。在对时间积分进行离散计算时,采用一种拟初始条件法,即在时间步迭代计算的过程中,将之前计算结果对当前时间步的影响都视作当前时间步的初始条件。在所取时间步长较小的情况下,这种处理方法容易导致数值结果不稳定,即每一步的计算误差会累计放大,最终导致计算崩溃。本文提出一种虚拟时刻方法以缓解这类数值不稳定现象,在该方法中,时间步长首先放大至合适尺度,计算某个虚拟时刻(往往在真实计算时刻之后)的温度和流量分布,再通过插值方法换算出真实时刻的温度和流量分布。在虚拟时刻点上的温度和流量计算过程中,边界已知温度或流量由真实时刻的温度或流量进行外插得到。本文简单证明了该方法在温度和流量随时间呈线性变化情况下的正确性,最后给出了两个分析实例,验证了该方法的准确性和稳定性。     

Mixed boundary value problem of heat conduction for infinite slab

Summary An axisymmetric steady state heat conduction boundary value problem having mixed boundary conditions on both faces of an infinite slab, is reduced to a pair of Fredholm integral equations of the second kind. For large values of h, the slab thickness, a solution correct to O(h ^–6) is obtained by expanding the kernels in power series.Presently at Imperial College, London.     

A three-dimensional optimal control problem in determining the boundary control heat fluxes

A three-dimensional optimal control algorithm in determining the strength of the unknown optimal boundary control heat fluxes utilizing the Conjugate Gradient Method (CGM) and a general purpose commercial code CFX4.2 is applied successfully in the present study based on the desired domain temperature distributions at the final time of heating. Results obtained by using the conjugate gradient method to solve this three-dimensional optimal control problems are justified based on the numerical experiments. Two different computational domains and two different desired temperature distributions are given and the corresponding optimal control heat fluxes are to be determined. Results show that the optimal control heat fluxes can always be obtained with any arbitrary initial guesses of the boundary fluxes.     

A boundary value problem in extended thermodynamics - one-dimensional steady flows with heat conduction

Abstract. One-dimensional steady flows with heat conduction, treated by the 13-moment theory of extended thermodynamics are considered. The usual well-posed boundary conditions for the corresponding problems in the Navier-Stokes-Fourier (NSF) theory are insufficient to give unique solutions. In order to have unique physically sensible solutions, minimization of the deviation of iterative approximations from the exact solutions, proposed in [1], is used as a criterion. Moreover, the solutions are shown to be invariant with respect to change of consistent boundary conditions - a requirement abide by our physical intuition. In the problems of plane shearing flow and Couette flow, the minimization with respect to two uncontrollable parameters is involved. The examples are carried out numerically, and the results are compared with the classical results of NSF theory.Received: 16 August 2002, Accepted: 20 April 2003, Published online: 5 December 2003PACS: 83.20.Lr, 83.50.Ax Correspondence to: I-S. Liu Dedicated to Professor Ingo Müller on the occasion of his 65th birthday     

Solution of non-linear inverse heat conduction problems using the method of lines

 Two space marching methods for solving the one-dimensional nonlinear inverse heat conduction problems are presented. The temperature-dependent thermal properties and the boundary condition on the accessible part of the boundary of the body are known. Additional temperature measurements in time are taken with a sensor located in an arbitrary position within the solid, and the objective is to determine the surface temperature and heat flux on the remaining part of the unspecified boundary. The methods have the advantage that time derivatives are not replaced by finite differences and the good accuracy of the method results from an appropriate approximation of the first time derivative using smoothing polynomials. The extension of the first method presented in this study to higher dimensions inverse heat conduction problems is straightforward. Received on 3 May 1999     

An inverse problem formulation of the immersed-boundary method

We formulate the immersed-boundary method (IBM) as an inverse problem. A control variable is introduced on the boundary of a larger domain that encompasses the target domain. The optimal control is the one that minimizes the mismatch between the state and the desired boundary value along the immersed target-domain boundary. We begin by investigating a naïve problem formulation that we show is ill-posed: in the case of the Laplace equation, we prove that the solution is unique, but it fails to depend continuously on the data; for the linear advection equation, even solution uniqueness fails to hold. These issues are addressed by two complimentary strategies. The first strategy is to ensure that the enclosing domain tends to the true domain, as the mesh is refined. The second strategy is to include a specialized parameter-free regularization that is based on penalizing the difference between the control and the state on the boundary. The proposed inverse IBM is applied to the diffusion, advection, and advection-diffusion equations using a high-order discontinuous Galerkin discretization. The numerical experiments demonstrate that the regularized scheme achieves optimal rates of convergence and that the reduced Hessian of the optimization problem has a bounded condition number, as the mesh is refined.     

Identification of a moving boundary for a heat conduction problem in a multilayer medium

In this paper, we give a uniqueness theorem for the moving boundary of a heat problem in a composite medium. Through solving the Cauchy problem of heat equation in each subdomain, we finally find an approximation to the moving boundary for one-dimensional heat conduction problem in a multilayer medium. The numerical scheme is based on the use of the method of fundamental solutions and a discrete Tikhonov regularization technique with the generalized cross-validation choice rule for a regularization parameter. Numerical experiments for five examples show that our proposed method is effective and stable.     

Non-fourier heat conduction in a plane slab with prescribed boundary heat flux

The non-stationary heat conduction in an infinitely wide plane slab with a prescribed boundary heat flux is studied. An arbitrary time dependent boundary heat flux is considered and a non-vanishing thermal relaxation time is assumed. The temperature and the heat flux density distributions are determined analytically by employing Cattaneo-Vernotte's constitutive equation for the heat flux density. It is proved that the temperature and the heat flux density distributions can be incompatible with the hypothesis of local thermodynamic equilibrium. A comparison with the solution which would be obtained by means of Fourier's law is performed by considering the limit of a vanishing thermal relaxation time.     

An iterative method to solve the heat transfer problem under the non-linear boundary conditions

The aim of the paper is to determine the approximation of the tangential matrix for solving the non-linear heat transfer problem. Numerical model of the strongly non-linear heat transfer problem based on the theory of the finite element method is presented. The tangential matrix of the Newton method is formulated. A method to solve the heat transfer with the non-linear boundary conditions, based on the secant slope of a reference function, is developed. The contraction mapping principle is introduced to verify the convergence of this method. The application of the method is shown by two examples. Numerical results of these examples are comparable to the ones solved with the Newton method and the commercial software COMSOL for the heat transfer problem under the radiative boundary conditions.