Singular boundary method for 2D and 3D heat source reconstruction
Introduction
Let be a bounded domain in , where d is the dimensionality of the space and is its boundary. The steady-state heat conduction is governed by Poisson equation with boundary conditions where T is the temperature, x means the spatial coordinates, a bounded domain, q the heat flux, and n the outward normal vector, and the measured data specified on the boundary . In this case, the right hand side denoting a heat source is sometimes unknown.
The identification of the source terms has received considerable attention in science, engineering and bioengineering. Some studies [1] required internal measurements of the temperature as supplement information, while in most of time, these measurements are intrusive. Nevertheless, non-invasive measurement only on the boundary is mathematically insufficient to guarantee the uniqueness of the source terms. Thus, some a priori information on the source should be assumed. El Badia [2] considered an inverse source problem for an anisotropic elliptic equation for which the source is a combination of monopolar and dipolar sources. Collara et al. [3] assumed that there is only one harmonic source, piecewise constant sources in their work, and identified it in non-homogeneous media with boundary measurements. In this study, the heat source is assumed to satisfy a second-order partial differential equation on a physical basis. Thus the source terms can be annihilated and the problem is transformed into a homogeneous fourth-order partial differential equation. Thus we reformulate the problem as
Formulation The source term satisfies a second order operator , namely, , then the problem of Eq. (1) can be reconstructed as The operator can be Laplace, Helmholtz or modified Helmholtz operators.
Then the higher-order direct problem (3), (4) can be solved with numerical methods, such as finite element method, boundary element method, method of fundamental solutions and other meshless methods [4], [5]. In this study, the singular boundary method (SBM) is applied to solve this problem with merely boundary discretization. The SBM is a boundary collocation method, which employs the fundamental solutions as the kernel functions. The key issue in the SBM is to desingularize the source singularities of the fundamental solutions with a finite value called origin intensity factor (OIF). Thus, the SBM avoids the sophisticated mathematics associated with the singular and nearly-singular integration in the BEM and perplexing choice of an optimal auxiliary boundary in the MFS [6], [7]. In the following, the formulation of the SBM will be established to identify the source terms.
Section snippets
Singular boundary method
In the SBM, the solution of partial differential equation (3) is approximated by a linear combination of the fundamental solutions with respect to source points: where and denote the field point and the jth source point, and are the fundamental solutions of the operators and (shown in Table 1). and are the unknown coefficients, which are determined by imposing the boundary
Numerical examples
This section tests two cases to illustrate the applicability and accuracy of the SBM with proposed OIFs to inverse source problem. For the evaluation of numerical accuracy, the root mean square error (RMSE) is employed with definition as where and are the exact and numerical solutions of source terms at the test point , respectively, and denotes the total number of test points in the domain. A relative noise with the intensity is imposed on the
Conclusions
In this study, the SBM in conjunction with the regularization method is employed to identify the source term in the steady-state heat conduction. The inverse problem is transformed to a direct fourth-order problem based on the assumption that the heat source is assumed to satisfy a second-order partial differential equation. Then the SBM approximates the solutions with a linear combination of fundamental solutions with respect to the source points, and the source terms can be derived from
Acknowledgments
This work is supported by the National Natural Science Foundation of China (Grant nos. 11662003, 11602114).
References (12)
- et al.
Application of method of fundamental solutions for elasto-plastic torsion of prismatic rods
Eng. Anal. Bound. Elem.
(2012) - et al.
Identification of piecewise constant sources in non-homogeneous media based on boundary measurements
Appl. Math. Model.
(2015) A high accuracy method for long-time evolution of acoustic wave equation
Appl. Math. Lett.
(2019)- et al.
Optimal sources in the MFS by minimizing a new merit function: Energy gap functional
Appl. Math. Lett.
(2018) - et al.
Singular boundary method for heat conduction problems with certain spatially varying conductivity
Comput. Math. Appl.
(2015) - et al.
A frequency domain formulation of the singular boundary method for dynamic analysis of thin elastic plate
Eng. Anal. Bound. Elem.
(2019)