A CQM-based BEM for transient heat conduction problems in homogeneous materials and FGMs

A method for numerical solution of time-domain boundary integral formulations of transient problems governed by the heat equation is presented. The heat conduction problem is analyzed considering homogeneous and non-homogeneous media. In the case of the non-homogeneous media, the conductor material is assumed to be a functionally graded material, i.e., the material properties vary spatially according to known smooth functions. For some specific spatial variations of the material properties, the fundamental solution and the boundary integral equation of the problem are obtained thanks to a change of variables that transforms the original problem to the standard heat conduction problem for homogeneous materials. For the treatment of time-dependent terms, the convolution quadrature method is adopted to approximate numerically the integral equation of the time-domain boundary element method. In the case that the responses are required at a large number of interior points, the convolution performed to calculate them is very time consuming. It is shown that the discrete convolution of the proposed formulation can be computed by means of the fast Fourier transform technique, which considerably reduces the computational complexity. Results for some transient heat conduction examples are presented to validate the numerical techniques studied.     

Use of 3D-transient analytical solution based on Green’s function to reduce computational time in inverse heat conduction problems

Inverse problems can be found in many areas of science and engineering and can be applied in different ways. Two examples can be cited: thermal properties estimation or heat flux function estimation in some engineering thermal process. Different techniques for the solution of inverse heat conduction problem (IHCP) can be found in literature. However, any inverse or optimization technique has a basic and common characteristic: the need to solve the direct problem solution several times. This characteristic is the cause of the great computational time consumed. In heat conduction problem, the time consumed is, usually, due to the use of numerical solutions of multidimensional models with refined mesh. In this case, if analytical solutions are available the computational time can be reduced drastically. This study presents the development and application of a 3D-transient analytical solution based on Green’s function. The inverse problem is due to the thermal properties estimation of conductors. The method is based on experimental determination of thermal conductivity and diffusivity using partially heated surface method without heat flux transducer. Originally developed to use numerical solution, this technique can, using analytical solution, estimate thermal properties faster and with better accuracy.     

热传导(对流-扩散)方程源项识别的粒子群优化算法

提出了利用粒子群优化(PSO)算法反演热传导方程与对流-扩散方程源项的一种新方法,在已有文献方法的基础上,求解出这两类方程正问题的解析解,再把源项识别问题转化为最优化问题,结合粒子群优化算法寻优求解.通过数值模拟与统计检验,结果表明,此方法可快速有效地实现热传导方程与对流-扩散方程源项的识别,并可推广应用到其它数学物理方程的源项或参数的反演识别.     

A modified quasi-boundary value method for solving the radially symmetric inverse heat conduction problem

In this paper, we consider a radially symmetric inverse heat conduction problem of determining the surface heat flux distribution from a fixed location inside a cylinder. This problem is ill-posed in the Hadamard sense and a conditional stability estimate is given for it. A modified quasi-boundary value regularization method is applied to formulate a regularized solution, and a sharp error estimate between the approximate solution and the exact solution is established by choosing a suitable regularization parameter. A numerical example is presented to verify the efficiency of the regularization method.     

Numerical Solution of the Retrospective Inverse Problem of Heat Conduction with the Help of the Poisson Integral

We consider the retrospective inverse problem that consists in determining the initial solution of the one-dimensional heat conduction equation with a given condition at the final instant of time. The solution of the problem is given in the form of the Poisson integral and is numerically realized by means of a quadrature formula leading to a system of linear algebraic equations with dense matrix. The results of numerical experiments are presented and show the efficiency of the numerical method including the case of the final condition with random errors.     

逆热传导方程的一个正则化方法

考虑了标准的一维逆热传导方程.问题是不适定的,即解不连续地依赖于数据.通过Fourier逼近的方法进行正则化处理,提出了一个新的算法,理论分析和数值实验均表明该算法是稳定的;该算法不仅保留了测量数据的部分高频成份,同时还具有相同的精度和计算复杂性.     

A modified method for a non-standard inverse heat conduction problem

A non-standard inverse heat conduction problem is considered. Data are given along the line x = 1 and the solution at x = 0 is sought. The problem is ill-posed in the sense that the solution (if it exists) does not depend continuously on the data. In order to solve the problem numerically it is necessary to employ some regularization method. In this paper, we study a modification of the equation, where a fourth-order mixed derivative term is added. Error estimates for this equation are given, which show that the solution of the modified equation is an approximation of the heat equation. A numerical implementation is considered and a simple example is given. Some numerical results show the usefulness of the modified method.     

Solution of a three-dimensional boundary-value problem for a fractional differential heat conduction equation

We consider a three-dimensional axisymmetric problem for a differential heat conduction equation with fractional time derivatives. Using the method of homogeneous solutions and integral transformations, we obtain an asymptotic and a numerical solution of the problem. The results of calculations are presented.     

三维热传导型半导体器件问题在局部加密网格上的有限差分格式

局部网格加密技术能很好地解决局部性很强的问题,半导体器件问题的解在半导体的p-n结附近有很强的局部性质.热传导型半导体器件瞬态问题的数学模型由四个方程组成的非线性偏微分方程组的初边值问题决定,电场位势方程是椭圆型的,电子和空穴浓度方程是抛物型的,温度方程是热传导型的.依据实际数值模拟的需要,提出了一类三维热传导型半导体问题在时间上进行局部加密的复合网格上的有限差分格式,并给出了电子、空穴浓度和温度的最大模误差估计以及数值算例.这些研究结果对半导体器件数值模拟的算法理论、实际应用和工程软件系统的研制,均具有重要的价值.     

An exact methodology for solving nonlinear diffusion equations based on integral transforms

This paper develops a new methodology for the solution of nonlinear diffusion equations. The solution technique is based on integral transforms and leads to exact numerical results. We apply the formal methodology to the problem of one-dimensional transient heat conduction. A new form of the heat equation is developed using a generalized expression for temperature-dependent thermal conductivity, based on a power-series expansion, for the three standard orthogonal coordinate systems. The resulting form of the heat equation suggests that the finite integral transform technique may reduce the dimensionality of the heat equation prior to the initiation of any numerical procedure. An example in a slab with linearly varying thermal conductivity is shown to produce exact results for the temperature distribution.     

二维热传导方程源强识别反问题的数值求解

本文研究了一类含有时间变量热源的二维热传导方程.利用有限元方法给出了数值求解过程,并在已知热源位置的前提下,根据某点的温度观测值,利用插值方法,将源强识别问题转化为参数反演问题,通过微分进化方法结合最佳摄动量法对源强识别反问题进行了数值模拟,结果表明所提出的方法是可行有效的.     

A regularization method for solving the radially symmetric backward heat conduction problem

This work is devoted to solving the radially symmetric backward heat conduction problem, starting from the final temperature distribution. The problem is ill-posed: the solution (if it exists) does not depend continuously on the given data. A modified Tikhonov regularization method is proposed for solving this inverse problem. A quite sharp estimate of the error between the approximate solution and the exact solution is obtained with a suitable choice of regularization parameter. A numerical example is presented to verify the efficiency and accuracy of the method.     

A problem on heat conduction in a wire enclosed by a cylindrical material solved by numerical inverse Laplace transform

The present paper describes the solution of a heat conduction problem by use of Laplace transformation and numerical inversion. Numerical results are given and compared with the asymptotic expressions and further the accuracy is discussed by comparing the values obtained from ten-point and eleven-point inversion formulas.     

Tracking a sharp crested wave front in hyperbolic heat transfer

In this article we investigate the numerical oscillations encountered when approximating the solution to the hyperbolic heat conduction equation. We consider a benchmark problem and show that it is not well-posed, unless a jump condition is specified. The alternative is to “smooth” the jump which leads to a sharp crested wave front, but with no discontinuity. To track the wave front we split the problem into auxiliary problems and solve these using different methods. The resulting solution is oscillation-free.     

The quasi‐reversibility method for a final value problem of the time‐fractional diffusion equation with inhomogeneous source

This paper is devoted to discuss a multidimensional backward heat conduction problem for time‐fractional diffusion equation with inhomogeneous source. This problem is ill‐posed. We use quasi‐reversibility regularization method to solve this inverse problem. Moreover, the convergence estimates between regularization solution and the exact solution are obtained under the a priori and the a posteriori choice rules. Finally, the numerical examples for one‐dimensional and two‐dimensional cases are presented to show that our method is feasible and effective.     

Splitting-based conjugate gradient method for a multi-dimensional linear inverse heat conduction problem

In this paper we consider a multi-dimensional inverse heat conduction problem with time-dependent coefficients in a box, which is well-known to be severely ill-posed, by a variational method. The gradient of the functional to be minimized is obtained by the aid of an adjoint problem, and the conjugate gradient method with a stopping rule is then applied to this ill-posed optimization problem. To enhance the stability and the accuracy of the numerical solution to the problem, we apply this scheme to the discretized inverse problem rather than to the continuous one. The difficulties with large dimensions of discretized problems are overcome by a splitting method which only requires the solution of easy-to-solve one-dimensional problems. The numerical results provided by our method are very good and the techniques seem to be very promising.     

一个抛物型方程侧边值问题的正则逼近解在一类Sobolev空间中的最优误差界

逆热传导问题(IHCP)是严重不适定问题,即问题的解(如果存在)不连续依赖于数据．但目前关于逆热传导问题的已有结果主要是针对标准逆热传导问题．文中给出了出现在实际问题中的一个抛物型方程侧边值问题,即一个含有对流项的非标准型逆热传导问题的正则逼近解一类Sobolev空间中的最优误差界．     

A numerical solution of a two-dimensional IHCP

In this paper, we will first study the existence and uniqueness of the solution of a two-dimensional inverse heat conduction problem (IHCP) which is severely ill-posed, i.e., the solution does not depend continuously on the data. We propose a stable numerical approach based on the finite-difference method and the least-squares scheme to solve this problem in the presence of noisy data. We prove the convergence of the numerical solution, then to regularize the resultant ill-conditioned linear system of equations, we apply the Tikhonov regularization 0th, 1st and 2nd method to obtain the stable numerical approximation to the solution. The stability and accuracy of the scheme presented is evaluated by comparison with the Singular Value Decomposition (SVD) method.     

基于时域精细积分算法的瞬态传热多宗量反演

基于有限元法和精细积分算法,提出了一种求解瞬态热传导多宗量反演问题的新方法．采用有限元法和精细积分算法分别对空间、时间变量进行离散,可以得到正演问题高精度的半解析数值模型,由此建立了多宗量反演的计算模式,并给出敏度分析的计算公式．对一维和二维的热物性参数、热源项、边界条件等进行了单宗量和多宗量的反演求解,初步考虑了初值和噪音等对反演结果的影响,数值算例验证了该方法的有效性．     

Parameter selection by discrete mollification and the numerical solution of the inverse heat conduction problem

A procedure for the numerical solution of the one-dimensional inverse heat conduction problem, based on the computaion of the solution associated with a suitable filtered version of the noisy data by discrete mollification is presented and a parameter choice criterion, which automatically determines the radius of mollification as a function of the amount of noise in the data, is introduced. Several numerical examples of interest are also analyzed, showing the accuracy and stability properties of the method.