Determination of the time-dependent reaction coefficient and the heat flux in a nonlinear inverse heat conduction problem

ABSTRACT

Diffusion processes with reaction generated by a nonlinear source are commonly encountered in practical applications related to ignition, pyrolysis and polymerization. In such processes, determining the intensity of reaction in time is of crucial importance for control and monitoring purposes. Therefore, this paper is devoted to such an identification problem of determining the time-dependent coefficient of a nonlinear heat source together with the unknown heat flux at an inaccessible boundary of a one-dimensional slab from temperature measurements at two sensor locations in the context of nonlinear transient heat conduction. Local existence and uniqueness results for the inverse coefficient problem are proved when the first three derivatives of the nonlinear source term are Lipschitz continuous functions. Furthermore, the conjugate gradient method (CGM) for separately reconstructing the reaction coefficient and the heat flux is developed. The ill-posedness is overcome by using the discrepancy principle to stop the iteration procedure of CGM when the input data is contaminated with noise. Numerical results show that the inverse solutions are accurate and stable.     

Numerical approximation of solution of nonhomogeneous backward heat conduction problem in bounded region

In this paper we consider a numerical approximation of solution of nonhomogeneous backward heat conduction problem (BHCP) in bounded region based on Tikhonov regularization method. Error estimate at t=0$t=0$ for this method is provided. According to the error estimate, a selection of regularization parameter is given. Meanwhile, a numerical implementation is described and the numerical results show that our algorithm is effective.     

应用QPSO算法求解二维热传导反问题

热传导反问题在国内研究起步较晚,研究方法有很多,但通常方法很难较好地接近全局最优。在经典的微粒群优化算法(PSO)的基础上,通过研究基于量子行为的微粒群优化算法(QPSO)提出了应用基于量子行为的微粒群优化算法进行二维热传导参数优化,具体介绍依据目标函数如何利用上述的算法去寻找最优参数组合。在具体应用中为了提高算法的收敛性和稳定性对算法进行了改进,并进行了大量实验,结果显示在解决热传导反问题优化问题中,基于QPSO算法的性能优越,证明QPSO在热传导领域具有很大的实际应用价值。     

QPSO算法与PSO算法求解二维热传导反问题

热传导反问题在国内研究起步较晚,研究方法有很多,但通常方法很难较好地接近全局最优.在介绍经典的微粒群优化算法(PSO)的基础上,研究基于量子行为的微粒群优化算法(QPSO)的二维热传导参数优化方法,具体介绍依据目标函数如何利用上述的算法去寻找最优参数组合.为了提高算法的收敛性和稳定性,在具体应用中对算法进行了改进,并进行了大量实验,结果显示在解决热传导反问题优化问题中,基于QPSO算法的性能比经典PSO算法更加优越,证明QPSO在热传导领域具有很大的实际应用价值.     

全变差正则化在抛物型方程初始条件反问题的应用

当反问题反演的函数不连续时,一般的正则化算法反演效果不令人满意,用全变差正则化方法对抛物型方程初始条件反问题进行求解,并进行了数值分析和数值模拟,结果显示数值解与真解吻合较好,表明该方法对于不连续函数求解具有高效、稳定等优点.     

Reconstruction of heat transfer coefficients using the boundary element method

This paper describes the reconstruction of the heat transfer coefficient (space, Problem I, or time dependent, Problem II) in one-dimensional transient inverse heat conduction problems from surface temperature or average temperature measurements. Since the inverse problem posed does not involve internal temperature measurements, this means that non-destructive testing of materials can be performed. In the formulation, convective boundary conditions relate the boundary temperature to the heat flux. Numerical results obtained using the boundary element method are presented and discussed.     

A method of fundamental solutions for radially symmetric and axisymmetric backward heat conduction problems

We propose and investigate an application of the method of fundamental solutions (MFS) to the radially symmetric and axisymmetric backward heat conduction problem (BHCP) in a solid or hollow cylinder. In the BHCP, the initial temperature is to be determined from the temperature measurements at a later time. This is an inverse and ill-posed problem, and we employ and generalize the MFS regularization approach [B.T. Johansson and D. Lesnic, A method of fundamental solutions for transient heat conduction, Eng. Anal. Boundary Elements 32 (2008), pp. 697–703] for the time-dependent heat equation to obtain a stable and accurate numerical approximation with small computational cost.     

A method of fundamental solutions for two-dimensional heat conduction

We investigate an application of the method of fundamental solutions (MFS) to heat conduction in two-dimensional bodies, where the thermal diffusivity is piecewise constant. We extend the MFS proposed in Johansson and Lesnic [A method of fundamental solutions for transient heat conduction, Eng. Anal. Bound. Elem. 32 (2008), pp. 697–703] for one-dimensional heat conduction with the sources placed outside the space domain of interest, to the two-dimensional setting. Theoretical properties of the method, as well as numerical investigations, are included, showing that accurate results can be obtained efficiently with small computational cost.     

Numerical solution of a one-dimensional inverse problem by the discontinuous Galerkin method

In this paper, we combine a Tikhonov regularization with a discontinuous Galerkin method to solve an inverse problem in one-dimension. We show that the regularization is simpler than in the case of the inversion using continuous finite elements. We numerically demonstrate that there exist optimal step sizes and polynomial degrees for inversion using the DG method. Numerical results are compared with those obtained by applying the standard finite element method with B-splines as a basis.     

A quasi-reversibility regularization method for the Cauchy problem of the Helmholtz equation

In this paper, a Cauchy problem for the Helmholtz equation is considered. It is known that such a problem is severely ill-posed, i.e. the solution does not depend continuously on the given Cauchy data. We propose a quasi-reversibility regularization method to solve it. Convergence estimates are established under two different a priori assumptions for an exact solution. Numerical results obtained by two different schemes show that our proposed methods work well.     

A numerical study of variable coefficient elliptic Cauchy problem via projection method

In this paper, we investigate a numerical method for the solution of an inverse problem of recovering lacking data on some part of the boundary of a domain from the Cauchy data on other part for a variable coefficient elliptic Cauchy problem. In the process, the Cauchy problem is transformed into the problem of solving a compact linear operator equation. As a remedy to the ill-posedness of the problem, we use a projection method which allows regularization solely by discretization. The discretization level plays the role of regularization parameter in the case of projection method. The balancing principle is used for the choice of an appropriate discretization level. Several numerical examples show that the method produces a stable good approximate solution.     

A spacetime discontinuous Galerkin method for hyperbolic heat conduction

Non-Fourier conduction models remedy the paradox of infinite signal speed in the traditional parabolic heat equation. For applications involving very short length or time scales, hyperbolic conduction models better represent the physical thermal transport processes. This paper reviews the Maxwell-Cattaneo-Vernotte modification of the Fourier conduction law and describes its implementation within a spacetime discontinuous Galerkin (SDG) finite element method that admits jumps in the primary variables across element boundaries with arbitrary orientation in space and time. A causal, advancing-front meshing procedure enables a patch-wise solution procedure with linear complexity in the number of spacetime elements. An h-adaptive scheme and a special SDG shock-capturing operator accurately resolve sharp solution features in both space and time. Numerical results for one spatial dimension demonstrate the convergence properties of the SDG method as well as the effectiveness of the shock-capturing method. Simulations in two spatial dimensions demonstrate the proposed method’s ability to accurately resolve continuous and discontinuous thermal waves in problems where rapid and localized heating of the conducting medium takes place.     

A meshless average source boundary node method for steady-state heat conduction in general anisotropic media

The average source boundary node method (ASBNM) is a recent boundary-type meshless method, which uses only the boundary nodes in the solution procedure without involving any element or integration notion, that is truly meshless and easy to implement. This paper documents the first attempt to extend the ASBNM for solving the steady-state heat conduction problems in general anisotropic media. Noteworthily, for boundary-type meshless/meshfree methods which depend on the boundary integral equations, whatever their forms are, a key but difficult issue is to accurately and efficiently determine the diagonal coefficients of influence matrices. In this study, we develop a new scheme to evaluate the diagonal coefficients via the pure boundary node implementation based on coupling a new regularized boundary integral equation with direct unknowns of considered problems and the average source technique (AST). Seven two- and three-dimensional benchmark examples are tested in comparison with some existing methods. Numerical results demonstrate that the present ASBNM is superior in the light of overall accuracy, efficiency, stability and convergence rates, especially for the solution of the boundary quantities.     

Richardson extrapolation applied to boundary element method results in a Dirichlet problem for the Laplace equation

Richardson extrapolation is used to improve the accuracy of the numerical solutions for the normal boundary flux and for the interior potential resulting from the boundary element method. The boundary integral equations arise from a direct boundary integral formulation for solving a Dirichlet problem for the Laplace equation. The Richardson extrapolation is used in two different applications: (i) to improve the accuracy of the collocation solution for the normal boundary flux and, separately, (ii) to improve the solution for the potential in the domain interior. The main innovative aspects of this work are that the orders of dominant error terms are estimated numerically, and that these estimates are then used to develop an a posteriori technique that predicts if the Richardson extrapolation results for applications (i) and (ii) are reliable. Numerical results from test problems are presented to demonstrate the technique.     

A complex variable boundary element method for a class of boundary value problems in anisotropic thermoelasticity

A boundary element method based on the Cauchy's integral formulae, called the complex variable boundary element method (CVBEM), is proposed for the numerical solution of boundary value problems governing plane thermoelastic deformations of anisotropic elastic bodies. The method is applicable for a wide class of problems which do not involve inertia or coupling effects and can be easily and efficiently implemented on the computer. It is applied to solve specific test problems.     

A truncated Legendre spectral method for solving numerical differentiation

A numerical differentiation problem for a given function with noisy data is discussed in this paper. A truncated spectral method has been introduced to deal with the ill-posedness of the problem. The theoretical analysis shows that the smoother the genuine solution, the higher the convergence rate of the numerical solution by our method. Numerical examples are also given to show the efficiency of the method.     

Worst case bounds on the point-wise discretization error in boundary element method for the elasticity problem

In this work, point-wise discretization error is bounded via interval approach for the elasticity problem using interval boundary element formulation. The formulation allows for computation of the worst case bounds on the boundary values for the elasticity problem. From these bounds the worst case bounds on the true solution at any point in the domain of the system can be computed. Examples are presented to demonstrate the effectiveness of the treatment of local discretization error in elasticity problem via interval methods.     

The backward problem for a time-fractional diffusion-wave equation in a bounded domain

This paper is devoted to solve the backward problem for a time-fractional diffusion-wave equation in a bounded domain. Based on the series expression of the solution for the direct problem, the backward problem for searching the initial data is converted into solving the Fredholm integral equation of the first kind. The existence, uniqueness and conditional stability for the backward problem are investigated. We use the Tikhonov regularization method to deal with the integral equation and obtain the series expression of the regularized solution for the backward problem. Furthermore, the convergence rate for the regularized solution can be proved by using an a priori regularization parameter choice rule and an a posteriori regularization parameter choice rule. Numerical results for five examples in one-dimensional case and two-dimensional case show that the proposed method is efficient and stable.