Next Article in Journal
Integral Transforms and the Hyers–Ulam Stability of Linear Differential Equations with Constant Coefficients
Previous Article in Journal
Revisiting the Absolute Configuration of Peyssonnoside A Using Vibrational Circular Dichroism Spectroscopy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identifying the Heat Source in Radially Symmetry and Axis-Symmetry Problems

Department of Mathematics, Northwest Normal University, Lanzhou 730070, China
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(2), 134; https://doi.org/10.3390/sym16020134
Submission received: 12 December 2023 / Revised: 14 January 2024 / Accepted: 18 January 2024 / Published: 23 January 2024

Abstract

:
This paper solves the inverse source problem of heat conduction in which the source term only varies with time. The application of the discrete regularization method, a kind of effective radial symmetry and axisymmetric heat conduction problem source identification that does not involve the grid integral numerical method, is put forward. Taking the fundamental solution as the fundamental function, the classical Tikhonov regularization method combined with the L-curve criterion is used to select the appropriate regularization parameters, so the problem is transformed into a class of ill-conditioned linear algebraic equations to solve with an optimal solution. Several numerical examples of inverse source problems are given. Simultaneously, a few numerical examples of inverse source problems are given, and the effectiveness and superiority of the method is shown by the results.

1. Introduction

The following heat equation is a suitable approximation for the transport, diffusion and conduction of natural substances that depend on shapes such as radially symmetric or axisymmetric vessels. We give a spatial one-dimensional version (1) of the model for this type of problem and a two-dimensional variant (2) of the spatial coordinates ρ and h:
u t ( ρ , t ) = 2 u ρ 2 ( ρ , t ) + 1 ρ u ρ ( ρ , t ) + f ( t ) , ρ 0 < ρ < P , 0 < t < T .
u t ( ρ , h , t ) = 2 u ρ 2 ( ρ , h , t ) + 1 ρ u ρ ( ρ , h , t ) + 2 u h 2 ( ρ , h , t ) + f ( t ) , ρ 0 < ρ < P , 0 < h < h 0 , 0 < t < T .
Here u represents the state variable, ρ indicates the radius of a radially symmetric or axisymmetric object under consideration, h is the height of the object, t is for time, and the right side of the equation f denotes physical laws [1,2], which means our case source terms. The problem is that we hope to recover the heat source term f ( t ) from the final time data, but, unfortunately, in many practical problems, the information of the source term f ( t ) is usually unknown, so it is impossible to completely invert the unknown source from the observed data of the real constrained boundary. Since there is no prior information about the function form of the unknown variable, this brings great difficulties to the estimation. It is well known that the thermal inverse source problem is a typical ill-posed problem, and the existence, uniqueness and stability of its solution are often difficult to achieve because the unknown solutions must contain some indirect measurement error that can be determined in the observed data [1]. For instance, the uniqueness and nonuniqueness of the solution to the problem of determining heat sources are discussed in [2], and its conditional stability is described in [3]. The ill-posed problems are divided into continuous ill-posed problems and discrete ill-posed problems; thus, there are generally two ways of thinking in the regularization process of solving ill-posed problems. The first is to effectively optimize the continuous problem via regularization methods to obtain an approximate solution and then discretize it appropriately, so as to obtain an approximate solution of the problem. This idea can be recorded as regularization and then discretization. Using the second method for choosing the appropriate discretization method for pathological algebraic equations and then effectively optimizing the regularization method, the approximate solution of the problem can be obtained. This idea can be remembered as entailing first discretization and then regularization. In the sense of considering practical applications, this paper mainly studies discrete ill-posed problems. The original problem is discretized into a linear system and a least-squares problem via the appropriate discretization methods, and then the problem is solved by combining effective regularization methods. Due to the vigorous development of science and technology, the inverse source problem appears in various fields of life, such as the identification of illegal wells in seawater intrusion, the scalp potential source detection recorded in brain power imaging and the determination of unknown pollution source characteristics in water resources management. Therefore, the research on the inverse source problem has attracted the extensive attention of scholars. Scholars have proposed many techniques to solve inverse source problems, including iterative regularization methods [4,5,6], boundary element methods (BEM) [7] and softening methods [8,9]. In addition, the sequence method [10] and the linear minimum squared error method [11] were also proposed. In these methods, discretizing partial differential equations on a grid is essential. Traditional grid-dependent methods, such as the finite difference method (FDM) and finite element method (FEM), rely on the grid in the region to support the solution to the problem. This method depends heavily on the quality of the grid, resulting in a significant amount of computation time. Although the boundary element method can double the size of the problem and greatly shorten the calculation time, it is very difficult to compute singular integrals on the boundary.
In recent decades, the meshless method has become the fastest growing and most widely used numerical calculation method in the world. It does not require the division of the problem into grids and can avoid the computational difficulties caused by the distortion of the large deformation analysis grid, which gives it special advantages over the grid-based approximation method in dealing with the problems of moving discontinuity, large deformation and high gradients. The solution’s accuracy and numerical stability are also high, so it can be used to solve some complex unsteady flow problems, such as fluid-solid coupling, seepage and so on. The grid-free method has been widely used in practical engineering applications, including the deformation analysis of building structures under earthquake action, the simulation of vehicle load in road traffic and the flow problem of oil-gas-water mixture in oil and gas wells. The meshless method has more than 30, including fundamental solutions (MFS) [12,13,14,15], the radial basis function method (RBFM) [16,17], the reproducing kernel particle method (RKPM) [18], the unit Galerkin method (EFGM) [19], etc. Among them, the MFS does not need grid and numerical integration at all. It is easier to extend from two-dimensional to three-dimensional, and, as a boundary method, it only requires the discretization of the boundary of the solution domain; therefore, it has been widely applied in fluid mechanics, biomechanics and some engineering problems and has achieved certain results. Reference [20] introduces the concept of MFS by solving the boundary value problem of the elliptic homogeneous equation, proving the applicability of MFS to some second- and fourth-order operators, and it conducts a systematic comprehensive evaluation. In order to avoid the singularity of the fundamental solution, the source points need to be distributed on a virtual boundary outside the physical boundary, but this hinders its practical application, so many researchers have proposed improved methods. In the literature [21,22], a non-singular general solution is used to replace the singular fundamental solution, which successfully solves the variable parameter problem and eliminates the well-known defect of offset boundary ambiguity. Another improved method is another boundary type collocation method based on the radial basis function (RBF) proposed by Chen and Tanaka, named BKM, which avoids the artificial boundary in the traditional fundamental solution, can obtain the symmetric matrix structure under certain conditions and speeds up the convergence rate [23]. At the same time, with the introduction of the concept of the origin intensity factor [24], the singular boundary method (SBM) has been proposed to eliminate the singularity of the fundamental solution [25], thus eliminating the fictional boundary of MFS. [26] combines the traditional MFS and localization concepts and proposes the basic solution method of localization (LMFS). LMFS uses the fundamental solution of the governing equation to derive the solution expression so that extremely accurate results can be obtained. In [27], the inverse source problem of radially symmetric and axisymmetric heat transfer equations in space is considered. In [28], the authors systematically introduce solutions for laminar flow and forced convection heat transfer in circular tubes. In the heat conduction problem, since some thermal fluids are incompressible, the method based on artificial compression can achieve a certain effect. A. Rostamzadeh et al. apply new method based on multi-dimensional anthropogenic characteristics (MACB) to a cylinder with heat transfer and forced convection to study steady-state and transient forced convection heat transfer around the cylinder [29]. M. Turkyilmazoglu et al. introduce some new tubular shapes considering a symmetric tubular boundary, show that a conventional pipe with an elliptical cross section can be modified without slip and then analyze the perturbation of the corresponding temperature field with a constant heat flux boundary around the slip parameter [30]. At present, in the field of inverse heat conduction problems [31,32,33,34,35] and Cauchy problems of partial differential equations [36,37,38,39], the combination of MFS and the regularization method has been successfully applied and achieved good research results.
In this paper, we extend the MFS method to solve the time-varying inverse heat source problem under arbitrary geometric conditions, successfully changing the limitation of only one unknown function, and then apply the MFS technique to the obtained equation. On this basis, a time parameter T is introduced to eliminate the singularity. It is proven that the precision of the numerical solution does not depend on the parameters through the analysis of numerical examples. Since the MFS discrete equations have incompatibity, the Tikhonov regularization algorithm is proposed to solve them, and regularization parameters are selected according to L-curve criteria [40,41,42,43,44]. Finally, several numerical examples of inverse source problems are given to verify the effectiveness of the proposed method.

2. Formulation of the Problem

In this section, we will give the specific formulas and simplify the procedure for problems (1) and (2).

2.1. One-Dimensional Heat Conduction Problem

We take into account a bounded domain ρ 0 < ρ < P , 0 < t < T . Let us consider the following two-dimensional heat Equation (1) where the source f ( t ) depends only on time.
u t ( ρ , t ) = 2 u ρ 2 ( ρ , t ) + 1 ρ u ρ ( ρ , t ) + f ( t ) , ρ 0 < ρ < P , 0 < t < T
u ( ρ , 0 ) = u 0 ( ρ ) , ρ 0 ρ P ,
u ( ρ 0 , t ) = d 1 ( t ) , 0 t T ,
u ( P , t ) = d 2 ( t ) , 0 t T ,
with the extra condition
u ( ρ 1 , t ) = s ( t ) ,
where ρ 1 ( ρ 0 , P ) , u 0 ( ρ ) , d 1 ( t ) , d 2 ( t ) and s ( t ) are given functions satisfying the compatibility conditions
u 0 ( 0 ) = 0 , u 0 ( 1 ) = d 2 ( 1 ) , u 0 ( ρ 1 ) = s ( 0 ) .
from (7) to (8), we get this solution. The inverse heat source problem is now formulated as follows: reconstruct the temperature u and the heat source function f ( r ) satisfying (3)–(6).
Remark 1.
Under suitable conditions (see [5,6]), the inverse problem (3)–(6) has a unique solution.
α ( t ) = 0 t f ( η ) d η ,
v ( ρ , t ) = u ( ρ , t ) + α ( t ) ,
The Formulas (3)–(7) are converted to the following formula:
v t ( ρ , t ) = 2 v ρ 2 ( ρ , t ) + 1 ρ v ρ ( ρ , t ) , ρ 0 < ρ < P , 0 < t < T ,
v ( ρ , 0 ) = u 0 ( ρ ) + α ( 0 ) , ρ 0 ρ P ,
v ( ρ 0 , t ) = d 1 ( t ) + α ( t ) , 0 t T
v ( P , t ) = d 2 ( t ) + α ( t ) , 0 t T ,
v ( ρ 1 , t ) = s ( t ) + α ( t ) , ρ 0 ρ 1 P ,
Substituting Equation (15) into Equations (11)–(14) yields the following problem:
v t ( ρ , t ) 2 v ρ 2 ( ρ , t ) 1 ρ v ρ ( ρ , t ) = 0 , ρ 0 < ρ < P , 0 < t < T ,
v ( ρ 0 , t ) v ( ρ 1 , t ) = d 1 ( t ) s ( t ) , 0 t T ,
v ( P , t ) v ( ρ 1 , t ) = d 2 ( t ) s ( t ) , 0 t T ,
v ( ρ , 0 ) = u 0 ( ρ ) , ρ 0 ρ P .
When v is obtained, the unknown heat source f ( t ) can be represented by the following formula:
f ( t ) = α t ( t ) = d 1 t ( t ) v t ( ρ 0 , t )

2.2. Two-Dimensional Heat Conduction Problem

On the basis of the above discussion, we can further extend the source identification method to the axisymmetric transient heat conduction problem. For completeness, we have made detailed instructions. We consider a model that also depends on h; that is, it is axisymmetric and the medium is isotropic. In the cylindrical coordinate system, we have the following formula.
u t ( ρ , h , t ) = 2 u ρ 2 ( ρ , h , t ) + 1 ρ u ρ ( ρ , h , t ) + 2 u h 2 ( ρ , h , t ) + f ( t ) , ρ 0 < ρ < P , 0 < h < h 0 , 0 < t < T .
u ( ρ , h , 0 ) = u 0 ( ρ , h ) , ρ 0 < ρ < P , 0 < h < h 0 ,
u ( ρ , h , t ) = d ( ρ , h , t ) , ρ { ρ 0 , P } , h { 0 , h 0 } , 0 < t < T ,
where d ( ρ , h , 0 ) = u 0 ( ρ , h ) for r { ρ 0 , P } , h { 0 , h 0 } to guarantee that the homogeneous initial condition (22) and the boundary condition (23) are compatible. Using the same technique as the separation of variables, if the functions f ( t ) , d ( ρ , h , t ) are known and fall into appropriate function spaces, we can employ an eigenfunction expansion to generate a unique solution to the direct problem.
Since the source function f ( t ) is unknown, prior information must be provided to recover the unknown source term f ( t ) . The extra data measured at a given point ( ρ 1 , h 1 ) is given by
u ( ρ 1 , h 1 , t ) = g ( t ) , ρ 0 < ρ 1 < P , 0 < h 1 < h 0 , 0 < t < T .
where u ( ρ 1 , h 1 , 0 ) = g ( 0 ) f o r ρ 0 < ρ 1 < P , 0 < h 1 < h 0 for compatibility.
In order to get a computational formula for the unknown source term f ( t ) , we introduce a new function v under the substitutions
v ( ρ , h , t ) = u ( ρ , h , t ) + α ( t ) ,
α t ( t ) = f ( t ) .
Then, the Formulas (21)–(24) can be converted to a homogeneous equation with non-local boundary conditions:
v t ( ρ , t ) = 2 v ρ 2 ( ρ , t ) + 1 ρ v ρ ( ρ , t ) + 2 v h 2 ( ρ , h , t ) , ρ 0 < ρ < P , 0 < h < h 0 , 0 < t < T ,
v ( ρ , h , 0 ) = u 0 ( ρ , h ) , ρ 0 ρ P , 0 h h 0 ,
v ( ρ , h , t ) = d ( ρ , h , t ) + α ( t ) , ρ { ρ 0 , P } , h { 0 , h 0 } , 0 < t < T ,
v ( ρ 1 , h 1 , t ) = g ( t ) + α ( t ) , ρ 0 < ρ 1 < R , 0 < h 1 < h 0 , 0 < t < T ,
By subtracting (30) from (29), it can be observed that the new function v satisfies the following problem:
v t ( ρ , t ) 2 v ρ 2 ( ρ , t ) 1 ρ v ρ ( ρ , t ) 2 v h 2 ( ρ , h , t ) = 0 , ρ 0 < ρ < P , 0 h h 0 , 0 < t < T ,
v ( ρ , h , 0 ) = u 0 ( ρ , h ) , ρ 0 ρ P , 0 h h 0 ,
v ( ρ , h , t ) v ( ρ 1 , h 1 , t ) = d ( t ) g ( t ) , ρ 0 < ρ 1 < P , 0 < h 1 < h 0 , ρ { ρ 0 , P } , h { 0 , h 0 } , 0 < t < T .
When v is obtained, the unknown heat source f ( t ) can be represented by the following formula:
f ( t ) = α t ( t ) = g t ( t ) v t ( ρ 1 , h 1 , t ) .

3. Methodology

The numerical approach for solving the source term identification problem in the radially and axisymmetric heat conduction equation is provided in this section. We use the MFS combined with the Tikhonov regularization. The L-curve rule is utilized to select a suitable regularization parameter. Because of the MFS approximation’s density, the numerical technique is believable. Finding the source locations and collocation boundary points is crucial in the MFS. We will provide the specifics. The literature [45] is where the main notion comes from.

3.1. The Methodology for One-Dimensional Problem

According to literature [35], the fundamental solution of equation v t a 2 v x x = 0 , 0 < x < 1 , 0 t t m a x is expressed as
F ( x , t ) = 1 2 a π t e x 2 4 a 2 t H ( t ) ,
Then, the fundamental solution of Equation (11) is expressed as
F 1 ( ρ , t ; ω , τ ) = H ( t τ ) 4 π ( t τ ) e ( ρ 2 + ω 2 ) / 4 ( t τ ) I 0 ρ ω 2 ( t τ ) ,
where H represents the Heaviside function and I 0 is the modified Bessel function of the first kind with order zero. Note that the special function I 0 ( 0 ) = 1 , lim x I 0 ( x ) = .
We approximate the solution of Equation (11) by using the following linear combination of fundamental solutions.
v ( ρ , t ) v N * ( ρ , t ) = i = 1 2 N 1 v i ( 0 ) F 1 ( ρ , t ; ρ 0 d 1 , τ i ) + i = 1 2 N 1 v i ( 1 ) F 1 ( ρ , t ; P + d 2 , τ i ) , ρ 0 ρ P , 0 t T ,
where N = 4 N 1 , d 1 , d 2 > 0 . The source point locations to be placed in time are chosen as
τ i = T + T ( 2 q i 1 ) 2 M 1 , f o r i = 1 , , 2 N 1 where M 1 = q N 1 with j > 0 fixed .
The 4 N 1 constant coefficients v i ( 0 ) and v i ( 1 ) for i = 1 , , 2 N 1 will be determined by the boundary conditions (17)–(19).
We choose the collocation points:
t n = n T M 1 , for n = 1 , , M 1
and
ρ l = ρ 0 + ( P ρ 0 ) l ( L + 1 ) , for l = 1 , , L .
It can be obtained from the above:
v ( ρ 0 , t n ) v ( ρ 1 , 0 ) = d 1 ( t n ) , n = 1 , , M 1 , v ( P , t n ) v ( ρ 1 , 0 ) = d 2 ( t n ) n = 1 , , M 1 , v ( ρ l , 0 ) = u 0 ( ρ l ) , l = 1 , , L .
Now (16)–(19) form a system,
Kv = b ,
of M = ( 2 M 1 + L ) linear equations with N = 4 N 1 unknowns. We use M N to constrain the unique solution of the system. In addition, considering the large number of conditions for K , we will regularize it using the traditional Tikhonov regularization method, namely using the following equations for regularization:
( K tr K + μ I ) v = K tr b
where K tr is the transpose of K , I is the identity matrix, and μ > 0 is a regularization parameter which will be chosen according to the famous L-curve criterion. Via (37), we have
v ( ρ 0 , t ) v N * ( ρ 0 , t ) = α * ( t ) = i = 1 2 N 1 v i ( 0 ) F 1 ( ρ 0 , t ; ρ 0 d 1 , τ i ) + i = 1 2 N 1 v i ( 1 ) F 1 ( ρ 0 , t ; P + d 2 , τ i ) , 0 t T ,
Therefore, the unknown source term is approximated as
f ( t ) d 1 ( t ) v t * ( ρ 0 , t ) .
According to (37), we need to calculate the first derivative of the fundamental solution (36) with respect to t. Simple calculations show that
F 1 t ( ρ , t ; ω , τ ) = H ( t τ ) 4 π ( t τ ) 2 e ( ρ 2 + ω 2 ) 4 ( t τ ) ( ρ 2 + ω 2 4 t 4 τ ) 4 ( t τ ) I 0 ρ ω 2 ( t τ ) ρ ω 2 ( t τ ) I 1 ρ ω 2 ( t τ )

3.2. The Methodology for Two-Dimensional Problem

The process can be extended to solve the axisymmetric heat conduction problem, which is determined by h and ρ and is therefore described by the axisymmetric two-dimensional heat Equation (27) in a space cylinder.
The following relation between the one-dimensional fundamental solution F 1 ( ρ , t ; ω , τ ) and the two-dimensional F 2 ( ρ , h , t ; ω , ξ , τ ) is given by
F 2 ( ρ , h , t ; ω , ξ , τ ) = 1 ( t τ ) e ( h ξ ) 2 4 ( t τ ) F 1 ( ρ , t ; ω , τ ) .
The solution of the axisymmetric heat Equation (27) can be approximated by a linear combination of the fundamental solutions.
v ( ρ , h , t ) v * ( ρ , h , t ) = i = 1 2 N 1 j = 1 2 N 2 ( ρ ) + 2 N 2 ( h ) v i ( j ) F 2 ( ρ , h , t ; ω j , ξ j , τ i ) , ρ 0 ρ P , 0 h h 0 , 0 t T .
Let Ω = ( ρ 0 , P ) × ( 0 , h 0 ) be the spatial domain of the solution and Γ = Ω . Let D E = ( ω i , ω j ) × ( ξ i , ξ j ) be a fictitious bounded domain in R 2 with D ¯ D E and the fictitious boundary D E = Γ E .
Now, we are going to put the source point on the virtual boundary Γ E and then put the juxtaposed point on the real boundary Γ . The time points used for generating the collocation points are given by
t n = n M 1 T , n = 0 , . . . , M 1 1 .
We all know that in order for MFS to get better results, the collocation point and the source point cannot be placed together in time. With that in mind, we will weave these time points together. Let the source points in time be
τ i = T + i q M 1 T T 2 q M 1 , i = 1 , . . . , 2 N 1 i f N 1 = q M 1 M 1 , T + i q M 1 T T 2 q M 1 i = 1 , . . . , 2 N 1 i f N 1 < M 1 = q N 1 .
There are two cases: N 1 M 1 or N 1 < M 1 . For the sake of simplicity, we assume there exists q > 0 such that N 1 = q M 1 or M 1 = q N 1 . This choice of t n and τ i guarantees that each are pairwise different and equally distributed in [ 0 , T ] and ( T , T ) , respectively. On D × { 0 , T } (see Equations (32) and (33)), let
ρ 0 , T ( l ρ , l h ) = ρ 0 + l ρ K ρ 1 ( P ρ 0 ) , l h K h 1 h 0 , l ρ = 0 , . . . , K ρ 1 , l h = 0 , . . . , K h 1 .
and on the boundary Γ (see Equation (32)), let
ρ m = ρ 0 + m M 2 ( ρ ) ( P ρ 0 ) , 0 , m = 1 , . . . , M 2 ( ρ ) , P , m M 2 ( ρ ) M 2 ( h ) h 0 , m = M 2 ( ρ ) + 1 , . . . , M 2 ( ρ ) + M 2 ( h ) , P m M 2 ( ρ ) M 2 ( h ) M 2 ( ρ ) ( P ρ 0 ) , h 0 , m = M 2 ( ρ ) + M 2 ( h ) + 1 , . . . , 2 M 2 ( ρ ) + M 2 ( h ) , ρ 0 , h 0 m 2 M 2 ( ρ ) M 2 ( h ) M 2 ( h ) h 0 , m = 2 M 2 ( ρ ) + M 2 ( h ) + 1 , . . . , 2 M 2 ( ρ ) + 2 M 2 ( h ) .
On the fictitious boundary Γ E , set
( ω j , ξ j ) = r 0 d ρ + j N 2 ( ρ ) ( 2 d ρ + P ρ 0 ) , d h , j = 1 , . . . , L 1 , P + d ρ , d h + j N 2 ( ρ ) N 2 ( h ) ( 2 d h + h 0 ) , j = L 1 + 1 , . . . , L 2 , P + d ρ j N 2 ( ρ ) N 2 ( h ) N 2 ( ρ ) ( 2 d ρ + P ρ 0 ) , h 0 + d h , j = L 2 + 1 , . . . , L 3 , ρ 0 d ρ , h 0 + d h j 2 N 2 ( ρ ) N 2 ( h ) N 2 ( h ) ( 2 d h + h 0 ) , j = L 3 + 1 , . . . , L 4 .
where d ρ , d h > 0 , L 1 = N 2 ( ρ ) , L 2 = N 2 ( ρ ) + N 2 ( h ) , L 3 = 2 N 2 ( ρ ) + N 2 ( h ) and L 4 = 2 N 2 ( ρ ) + 2 N 2 ( h ) . Boundary conditions are then imposed at the collocation points to determine the constant coefficients v i ( j ) , as given in (43). We have M = M 1 ( 2 M 2 ( ρ ) + M 2 ( h ) ) + K ρ K h collocation points and N = 2 N 1 ( 2 N 2 ( ρ ) + 2 N 2 ( h ) ) source points. To obtain a unique solution, we restrict N M . Now, we can form a linear system.
v ( ρ , h , 0 ) = u 0 ( ρ , h , t n ) , ( ρ , h ) = ρ m ; v ( ρ , h , t n ) v ( ρ 1 , h 1 , t n ) = g ( ρ , h ) , ( ρ , h ) = ρ 0 , T l ρ , l h , n = 1 , , M 1 .
Similarly to (35), we can draw a set of linear Equation (38). Then, we use the Tikhonov regularization method, as in (38).
The solution of problems (21)–(23) is simultaneously obtained by
u * ( ρ , h , t ) = v * ( ρ , h , t ) α * ( t )
and
f * ( t ) = α t * = g t ( t ) v t * ( ρ 1 , h 1 , t ) .
F 2 ( ρ , h , t ; ω , ξ , τ ) = H ( t τ ) 4 π ( t τ ) 3 / 2 e ρ 2 + ω 2 + ( h ξ ) 2 4 ( t τ ) I 0 ρ ω 2 ( t τ ) .
F 2 t ( ρ , t ; ω , τ ) = H ( t τ ) 32 π 3 2 ( t τ ) 7 2 e ρ 2 + ω 2 + ( h ξ ) 2 4 ( t τ ) [ ρ 2 + ω 2 + ( h ξ ) 2 6 ( t τ ) ] I 0 ρ ω 2 ( t τ ) 2 ρ ω I 1 ρ ω 2 ( t τ )

4. Numerical Results

In all of the following examples, we will assume that T = 1 . We use the function rand given in Matlab to generate the noisy data g ˜ i = g i ( 1 + δ · r a n d ( i ) ) , where g i is the exact data, and r a n d ( i ) denotes a random number from the uniform distribution of interval ( 1 , 1 ) . The magnitude δ indicates the noise level of the measurement data.
In order to test the accuracy of the approximate solution, the root mean square error (RMS) and the relative root mean square error (RES) are defined as
R M S ( f ) = 1 N t i = 1 N t ( f ( t i ˜ ) f * ( t i ˜ ) ) 2
R E S ( f ) = i = 1 N t ( f ( t i ˜ ) f * ( t i ˜ ) ) 2 i = 1 N t ( f ( t i ˜ ) ) 2
where N t is the total number of testing points in 0 t 1 , and f ( t i ˜ ) , f * ( t i ˜ ) are, respectively, the exact and approximated value at these points.

4.1. Numerical Examples for MFS: Radially Symmetric Heat Equation

Example 1. 
The exact solution to the problems (3)–(7) is given by:
u ( ρ , t ) = ρ 2 + t 2 ,
f ( t ) = 2 t 4 , 0 t 1 .
As can be seen from Table 1, the farther ρ 1 is away from the regional center within a certain extent, the more accurate the calculation results are, but when ρ 1 increases to a certain value, the accuracy of the calculation results decreases. It can be seen that the calculated results are either too far or too close to the center of the region at ρ 1 . In addition, at a fixed radius, the lower the noise δ value added to the temperature data is, the smaller the value of R E S ( f ) is.
For Example 1, Figure 1 shows the comparison between the exact solution and the approximate solution at no noise and at [1%, 3%, 5%] noise levels, respectively. We can clearly see from the figure that the results are ideal under the conditions of no noise and 1 % noise, while the approximate solution affected by noise at 3% and 5% levels will result in excessive error. Figure 2 shows that under different noise levels δ , the Tikhonov method is used to regularize the L-curve function of the problem obtained by the MFS discretization system to determine the regularized parameters. Figure 3, Figure 4 and Figure 5 specifically analyze the influences of d 1 , d 2 , M N l and δ on the results, and Figure 5 shows that the relative noise level is roughly positively correlated with the results error. Figure 6, Figure 7 and Figure 8 of Example 2 is similar to Example 1, and we will not go into too much detail here.
Example 2. 
Here we have another example:
u ( ρ , t ) = e t J 0 ( ρ ) + e t ,
and
f ( t ) = e t .
Example 3. 
We consider another problem given by
u ( ρ , t ) = ρ 2 + 2 t + s i n ( 2 π t ) ,
and
f ( t ) = 2 + 2 π c o s ( 2 π t ) .
From Table 2, similar conclusions can be observed to those from Table 1.
As can be seen from Figure 9, Figure 10 and Figure 11, the result of Example 3 is obviously closer to ideal than the previous two examples, and the resulting error is less sensitive to the changes in d 1 , d 2 , M N L and δ , which indicates that the stability of the solution can be guaranteed. Numerical experiments show that this method is effective.
Example 4. 
We consider the following direct problem:
u ( ρ , 0 ) = 0 , 0 ρ 1
u ( 0 , t ) = 0 , 0 t 1
u ( 1 , t ) = 0 , 0 t 1
with the heat source function
f ( t ) = 1 , t [ 0 , 0.25 ) , 1 , t [ 0.25 , 0.5 ) , 1 , t [ 0.5 , 0.75 ) , 1 , t [ 0.75 , 1 ] .
We solve the problem with MFS, with M 1 = N 1 = 50 , L = 100 , d 1 = d 2 = 1.3 when the solution is u ( ρ 1 = 0.7 , t ) = s ( t ) . Then, we use the data s ( t ) from the forward problem to solve the inverse problem. Figure 12 shows the comparison between the exact solution and the numerical results. It can be seen from the figure that the numerical approximation is not ideal, but it is reasonably consistent with (60). Compared with the weak relationship between the result and the noise level in Figure 12, it can be seen from Figure 13 that M N L has a greater impact on the result. As the value of M N L increases, the result becomes closer to ideal.

4.2. Numerical Examples for MFS: Axisymmetric Heat Equation

It can be seen from the analysis in Section 2 that, in this case, the nature of the problem becomes worse, which also increases the difficulty of numerical calculation.
Example 5. 
In this example, the axisymmetric equation is formulated in a solid cylinder with 0 < ρ < 1 and 0 < h < 1 . We use 400 c o l l o c a t i o n p o i n t s a n d 300 source points to analyze the factors affecting the error of the results at a noise level of δ = 1 % . The exact solution is given by
u ( ρ , h , t ) = ( ρ 2 + h 2 ) + 3 t 2 , 0 ρ 1 , 0 h 1 , 0 t 1 ,
and
f ( t ) = 6 t 6 , 0 t 1 .
Figure 14 illustrates the relationship between the exact and numerical solution for source f ( t ) on d ρ = 3 , d h = 5 , ρ 1 = 0.1 , h 1 = 0.5 . All the measurement data g i are computed via the exact solution (3.10). The relationship between ρ , h and the relative error of the result is shown in Figure 15. Figure 16 considers the relationship between d ρ , d h and the relative error of the result under the condition of ρ 1 = 0.6 and h 1 = 0.3 .
Example 6. 
The domain is in a solid cylinder with 0 < ρ < 1 and 0 < h < 1 . The exact solution is given by
u ( ρ , h , t ) = 1 + h + ( ρ 2 + h 2 ) 1 / 2 sin π 4 ρ 2 + h 2 e π 2 t / 16 + t e t , 0 ρ 1 , 0 h 1 , 0 t 1 ,
and
f ( t ) = t e t + t , 0 t 1 .
Figure 17 illustrates the relationship between the exact and numerical solution for source f ( t ) on d ρ = 3 , d h = 4 , ρ 1 = 0.6 , h 1 = 0.3 . Figure 18 and Figure 19 respectively show the relationship between ρ 1 , h 1 , d ρ , d h and the relative error of the results. It can be seen that the result is closer to ideal than that of Example 5.

5. Conclusions

In this paper, we extend MFS to solve the inverse heat source problem, which only varies with time in the radially symmetric and axisymmetric cases. We solve the problem with only one unknown function and introduce the time parameter T, which does not affect the precision of the numerical solution to avoid the singularity. We apply a MFS technique to transform the problem into a discrete system by selecting suitable location source points and matching boundary points, and realize the source term identification of the regional radial symmetry and axisymmetric heat conduction problem based on the Tikhonov regularization method and L-curve criterion. We study several numerical examples of one- and two-dimensional radially symmetric and axisymmetric heat source conduction problems. For one-dimensional cases, we also consider the continuous and discontinuous source terms, respectively. The results show that the proposed method provides an effective and accurate solution for such inverse heat conduction problems.
This method can extend MFS to combine with non-classical regularization methods, which will be accomplished in subsequent works. Secondly, there is no error estimation for solving the inverse source problem for the MSF method combined with Tikhonov regularization, which needs to be studied further.

Author Contributions

Conceptualization, X.X.; methodology, X.X.; software, Y.S.; validation, Y.S.; formal analysis, Y.S.; writing—original draft preparation, Y.S.; writing—review and editing, Y.S.; visualization, Y.S.; supervision, X.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kirsch, A. An Introduction to the Mathematical Theory of Inverse Problems; Springer: New York, NY, USA, 2011. [Google Scholar]
  2. Denisov, A.M. Uniqueness and nonuniqueness of the solution to the problem of determining the source in the heat equation. Comput. Math. Math. Phys. 2016, 56, 1737–1742. [Google Scholar] [CrossRef]
  3. Yamamoto, M. Conditional stability in determination of force terms of heat equations in a rectangle. Math. Comput. Model. 1993, 18, 79–88. [Google Scholar] [CrossRef]
  4. Johansson, B.T.; Lesnic, D. Determination of a spacewise dependent heat source. J. Comput. Appl. Math. 2007, 209, 66–80. [Google Scholar] [CrossRef]
  5. Johansson, B.T.; Lesnic, D. A variational method for identifying a spacewise-dependent heat source. IMA J. Appl. Math. 2007, 72, 748–760. [Google Scholar] [CrossRef]
  6. Johansson, B.T.; Lesnic, D. A procedure for determining a spacewise dependent heat source and the initial temperature. Appl. Anal. 2008, 87, 265–276. [Google Scholar] [CrossRef]
  7. Farcas, A.; Lesnic, D. The boundary-element method for the determination of a heat source dependent on one variable. J. Eng. Math. 2006, 54, 375–388. [Google Scholar] [CrossRef]
  8. Yi, Z.; Murio, D.A. Source term identification in 1-D IHCP. Comput. Math. Appl. 2004, 47, 1921–1933. [Google Scholar] [CrossRef]
  9. Yi, Z.; Murio, D.A. Source term identification in 2-D IHCP. Comput. Math. Appl. 2004, 47, 1517–1533. [Google Scholar] [CrossRef]
  10. Yang, C.Y. A sequential method to estimate the strength of the heat source based on symbolic computation. Int. J. Heat Mass. Transfer 1998, 41, 2245–2252. [Google Scholar] [CrossRef]
  11. Yang, C.Y. Solving the two-dimensional inverse heat source problem through the linear least-squares error method. Int. J. Heat Mass Transfer 1998, 41, 393–398. [Google Scholar] [CrossRef]
  12. Golberg, M.A.; Chen, C.S. The method of fundamental solutions for potential hemholtz and diffusion problems. In Boundary Integral Methods: Numerical and Mathematical Aspects; Computational Mechanics Publication: Southampton, UK, 1999; pp. 103–176. [Google Scholar]
  13. Fairweather, G.; Karageorghis, A. The method of fundamental solutions for elliptic boundary value problems. Adv. Comput. Math. 1998, 9, 69–95. [Google Scholar] [CrossRef]
  14. Powell, M.D.J. Radial basis functions for multivariable interpolation: A review. In Numerical Analysis, Longman Scientific and Technical; Griffiths, D.F., Watson, G.A., Eds.; Longman Scientific and Technical: Harlow, UK, 1987; pp. 223–241. [Google Scholar]
  15. Kupradze, V.D. A method for the approximate solution of limiting problems in mathematical physics. USSR Comput. Maths. Math. Phys. 1964, 4, 199–205. [Google Scholar] [CrossRef]
  16. Young, D.L.; Tsai, C.C.; Murugesan, K.; Fan, C.M.; Chen, C.W. Time-dependent fundamental solutions for homogeneous diffusion problems. Eng. Anal. Bound. Elem. 2004, 28, 1463–1473. [Google Scholar] [CrossRef]
  17. Wu, Z.; Hermite, C. Birkhoff interpolation for scattered data by radial basis function, Approxim. Theory Appl. 1992, 8, 1–10. [Google Scholar]
  18. Liu, W.K.; Jun, S.; Li, S.; Adee, J.; Belytschko, T. Reproducing kernel particle methods for structural dynamics. Internat. J. Numer. Methods Engrg. 1995, 38, 1655–1679. [Google Scholar] [CrossRef]
  19. Li, X.; Li, S. Effect of an efficient numerical integration technique on the element-free Galerkin method. Appl. Numer. Math. 2023, 192, 204–225. [Google Scholar] [CrossRef]
  20. Bogomolny, A. Fundamental solutions method for elliptic boundary value problems. SIAM J. Numer. Anal. 1985, 22, 644–669. [Google Scholar] [CrossRef]
  21. Chen, J.T.; Chen, I.L.; Chen, K.H.; Lee, Y.T.; Yeh, Y.T. A meshless method for free vibration analysis of circular and rectangular clamped plates using radial basis function. Eng. Anal. Bound. Elem. 2004, 28, 535–545. [Google Scholar] [CrossRef]
  22. Kang, S.W.; Lee, J.M.; Kang, Y.J. Vibration analysis of arbitrary shape membranes using non-dimensional dynamic influence function. J. Sound Vib. 1999, 221, 117–132. [Google Scholar] [CrossRef]
  23. Chen, W.; Tanaka, M. A meshless, integration-free, and boundary-only RBF technique. Comput. Math. Appl. 2002, 43, 379–391. [Google Scholar] [CrossRef]
  24. Chen, W.; Tanaka, M. New insights in boundary-only and domain-type RBF methods. Int. J. Nonlinear Sci. Numer. Simul. 2000, 1, 145–151. [Google Scholar] [CrossRef]
  25. Chen, W. Singular boundary method: A novel, simple, meshfree, boundary collocation numerical method. Chin. J. Solid Mech. 2009, 30, 592–599. [Google Scholar]
  26. Fan, C.M.; Huang, Y.K.; Chen, C.S.; Kuo, S.R. Localized method of fundamental solutions for solving two-dimensional Laplace and biharmonic equations. Eng. Anal. Bound. Elem. 2019, 101, 188–197. [Google Scholar] [CrossRef]
  27. Amin, M.E.; Xiong, X. Source identification problems for radially symmetric and axis-symmetric heat conduction equations. Appl. Numer. Math. 2019, 138, 1–18. [Google Scholar] [CrossRef]
  28. Shah, R.K.; London, A.L. Laminar Flow Forced Convection in Ducts: A Source Book for Compact Heat Exchanger Analytical Data; Academic Press: Cambridge, MA, USA, 1978. [Google Scholar]
  29. Rostamzadeh, A. Towards Multidimensional Artificially Characteristic-Based Scheme for Incompressible Thermo-Fluid Problems. Mechanika 2018, 23, 826–834. [Google Scholar] [CrossRef]
  30. Turkyilmazoglu, M.; Duraihem, F.Z. Full Solutions to Flow and Heat Transfer from Slip-Induced Microtube Shapes. Micromachines 2023, 14, 894. [Google Scholar] [CrossRef]
  31. Hon, Y.C.; Wei, T. A fundamental solution method for inverse heat conduction problem. Eng. Anal. Bound. Elem. 2004, 28, 489–495. [Google Scholar] [CrossRef]
  32. Hon, Y.C.; Wei, T. Numerical computation for multidimensional inverse heat conduction problem. CMES Comput. Model. Eng. Sci. 2005, 7, 119–132. [Google Scholar]
  33. Johansson, B.T.; Lesnic, D. A method of fundamental solutions for transient heat conduction. Eng. Anal. Bound. Elem. 2008, 32, 697–703. [Google Scholar] [CrossRef]
  34. Chantasiriwan, S. Methods of fundamental solutions for time-dependent heat conduction problems. Int. J. Numer. Meth. Eng. 2006, 66, 147–165. [Google Scholar] [CrossRef]
  35. Yan, L.; Fu, C.L.; Yang, F.L. The method of fundamental solutions for the inverse heat source problem. Eng. Anal. Bound. Elem. 2008, 32, 216–222. [Google Scholar] [CrossRef]
  36. Wei, T.; Hon, Y.C.; Ling, L. Method of fundamental solutions with regularization techniques for Cauchy problems of elliptic operators. Eng. Anal.Boundary Elem. 2007, 31, 373–385. [Google Scholar] [CrossRef]
  37. Marin, L.; Lesnic, D. The method of fundamental solutions for the Cauchy problem associated with two-dimensional Helmholtz-type equations. Comput. Struct. 2005, 83, 267–278. [Google Scholar] [CrossRef]
  38. Marin, L.; Lesnic, D. The method of fundamental solutions for inverse boundary value problems associated with the two-dimensional biharmonic equation. Math. Comput. Model. 2005, 42, 261–278. [Google Scholar] [CrossRef]
  39. Marin, L.; Lesnic, D. The method of fundamental solutions for nonlinear functionally graded materials. Int. J. Solids Struct. 2007, 44, 6878–6890. [Google Scholar] [CrossRef]
  40. Hansen, P.C. Rank-Deficient and Discrete Ill-Posed Problems; SIAM: Philadelphia, PA, USA, 1998. [Google Scholar]
  41. Morozov, V.A. Methods for solving Incorrectly Posed Problems; Springer: New York, NY, USA, 1984. [Google Scholar]
  42. Engl, H.W.; Hanke, M.; Neubauer, A. Regularization of Inverse Problems, Mathematics, and Its Applications; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1996; Volume 357. [Google Scholar]
  43. Tautenhahn, U.; Hämarik, U. The use of monotonicity for choosing the regularization parameter in ill-posed problems. Inverse Probl. 1999, 15, 1487–1505. [Google Scholar] [CrossRef]
  44. Golub, G.; Heath, M.; Wahba, G. Generalized cross-validation as a method for choosing a good ridge parameter. Technometrics 1979, 21, 215–223. [Google Scholar] [CrossRef]
  45. Wei, T.; Hon, Y.C. Numerical differentiation by radial basis functions approximation. Adv. Comput. Math. 2007, 27, 247–272. [Google Scholar] [CrossRef]
Figure 1. Example 1. In the exact solution ( ) and numerical solution ( ) for f ( t ) with three different levels of noise [1%,3%,5%] with μ = [ 6.89 × 10 5 , 3.66 × 10 4 , 8.45 × 10 4 ] in (bd), respectively, and with noiseless data (a) with μ = 1.00 × 10 9 , the parameters are M 1 = 10 , N 1 = 10 , L = 20 and d 1 = d 2 = 2.5 , and the RMSEs of (ad) are [ 2.8304 × 10 5 , 0.0199 , 0.0673 , 0.1117 ] , respectively.
Figure 1. Example 1. In the exact solution ( ) and numerical solution ( ) for f ( t ) with three different levels of noise [1%,3%,5%] with μ = [ 6.89 × 10 5 , 3.66 × 10 4 , 8.45 × 10 4 ] in (bd), respectively, and with noiseless data (a) with μ = 1.00 × 10 9 , the parameters are M 1 = 10 , N 1 = 10 , L = 20 and d 1 = d 2 = 2.5 , and the RMSEs of (ad) are [ 2.8304 × 10 5 , 0.0199 , 0.0673 , 0.1117 ] , respectively.
Symmetry 16 00134 g001
Figure 2. Example 1. The L-curve plot for δ = 1 % ( ) , 6.89 × 10 5 , δ = 3 % ( ) , 3.66 × 10 4 and δ = 5 % ( * ) , 8.45 × 10 4 .
Figure 2. Example 1. The L-curve plot for δ = 1 % ( ) , 6.89 × 10 5 , δ = 3 % ( ) , 3.66 × 10 4 and δ = 5 % ( * ) , 8.45 × 10 4 .
Symmetry 16 00134 g002
Figure 3. Example 1. The RES for ( d 1 , d 2 ) [ 1 , 5 ] × [ 1 , 5 ] , M 1 = 10, N 1 = 10, L = 20 and δ = 1 % .
Figure 3. Example 1. The RES for ( d 1 , d 2 ) [ 1 , 5 ] × [ 1 , 5 ] , M 1 = 10, N 1 = 10, L = 20 and δ = 1 % .
Symmetry 16 00134 g003
Figure 4. Example 1. The RES for d 1 = d 2 = 2.5 , L = 2 M 1 , δ = 1 % and ( M 1 , N 1 ) [ 10 , 50 ] × [ 10 , 50 ] , M = 4 M 1 collocation points and N = 4 N 1 source points.
Figure 4. Example 1. The RES for d 1 = d 2 = 2.5 , L = 2 M 1 , δ = 1 % and ( M 1 , N 1 ) [ 10 , 50 ] × [ 10 , 50 ] , M = 4 M 1 collocation points and N = 4 N 1 source points.
Symmetry 16 00134 g004
Figure 5. Example 1. RES versus relative noise levels.
Figure 5. Example 1. RES versus relative noise levels.
Symmetry 16 00134 g005
Figure 6. Example 2. In the exact solution ( ) and numerical solution ( ) for f ( t ) with three different levels of noise [1%,3%,5%] with μ = [ 7.04 × 10 4 , 3.71 × 10 3 , 5.29 × 10 3 ] in (bd), respectively, and with noiseless data (a) with μ = 1.00 × 10 9 , the parameters are M 1 = 10 , N 1 = 10 , L = 20 and d 1 = d 2 = 2.5 , and the RMSEs of (ad) are [ 1.1728 × 10 4 , 0.0160 , 0.0852 , 0.1031 ] , respectively.
Figure 6. Example 2. In the exact solution ( ) and numerical solution ( ) for f ( t ) with three different levels of noise [1%,3%,5%] with μ = [ 7.04 × 10 4 , 3.71 × 10 3 , 5.29 × 10 3 ] in (bd), respectively, and with noiseless data (a) with μ = 1.00 × 10 9 , the parameters are M 1 = 10 , N 1 = 10 , L = 20 and d 1 = d 2 = 2.5 , and the RMSEs of (ad) are [ 1.1728 × 10 4 , 0.0160 , 0.0852 , 0.1031 ] , respectively.
Symmetry 16 00134 g006
Figure 7. Example 2. The RES for ( d 1 , d 2 ) [ 1 , 5 ] × [ 1 , 5 ] M 1 = 10 , N 1 = 10 , L = 20 and δ = 1 % .
Figure 7. Example 2. The RES for ( d 1 , d 2 ) [ 1 , 5 ] × [ 1 , 5 ] M 1 = 10 , N 1 = 10 , L = 20 and δ = 1 % .
Symmetry 16 00134 g007
Figure 8. Example 2. The RES for d 1 = d 2 = 2 , L = 2 M 1 , δ = 1 % and ( M 1 , N 1 ) [ 10 , 50 ] × [ 10 , 50 ] , M = 4 M 1 collocation points and N = 4 N 1 source points.
Figure 8. Example 2. The RES for d 1 = d 2 = 2 , L = 2 M 1 , δ = 1 % and ( M 1 , N 1 ) [ 10 , 50 ] × [ 10 , 50 ] , M = 4 M 1 collocation points and N = 4 N 1 source points.
Symmetry 16 00134 g008
Figure 9. Example 3. In the exact solution ( ) and numerical solution ( ) for f ( t ) with three different levels of noise [1%,3%,5%] with μ = [ 6.89 × 10 5 , 3.66 × 10 4 , 8.45 × 10 4 ] in (bd), respectively, and with noiseless data (a) with μ = 1.00 × 10 9 , the parameters are M 1 = 10 , N 1 = 10 , L = 20 and d 1 = d 2 = 2.5 , and the RMSEs of (ad) are [ 2.8304 × 10 5 , 0.0199 , 0.0673 , 0.1117 ] , respectively. Of note are the (a,b) in Figure 9: although there looks to be no difference, they actually are not equal at some point. Because Example 2 is not sensitive to the change in noise data and the change in the noise is small, the change between them is not too large.
Figure 9. Example 3. In the exact solution ( ) and numerical solution ( ) for f ( t ) with three different levels of noise [1%,3%,5%] with μ = [ 6.89 × 10 5 , 3.66 × 10 4 , 8.45 × 10 4 ] in (bd), respectively, and with noiseless data (a) with μ = 1.00 × 10 9 , the parameters are M 1 = 10 , N 1 = 10 , L = 20 and d 1 = d 2 = 2.5 , and the RMSEs of (ad) are [ 2.8304 × 10 5 , 0.0199 , 0.0673 , 0.1117 ] , respectively. Of note are the (a,b) in Figure 9: although there looks to be no difference, they actually are not equal at some point. Because Example 2 is not sensitive to the change in noise data and the change in the noise is small, the change between them is not too large.
Symmetry 16 00134 g009
Figure 10. Example 3. The RES for ( d 1 , d 2 ) [ 1 , 5 ] × [ 1 , 5 ] M 1 = 10 , N 1 = 10 , L = 20 and δ = 1 % .
Figure 10. Example 3. The RES for ( d 1 , d 2 ) [ 1 , 5 ] × [ 1 , 5 ] M 1 = 10 , N 1 = 10 , L = 20 and δ = 1 % .
Symmetry 16 00134 g010
Figure 11. Example 3. The RES for d 1 = d 2 = 2.5 , L = 2 M 1 δ = 1 % and ( M 1 , N 1 ) [ 10 , 50 ] × [ 10 , 50 ] , M = 4 M 1 collocation points and N = 4 N 1 source points.
Figure 11. Example 3. The RES for d 1 = d 2 = 2.5 , L = 2 M 1 δ = 1 % and ( M 1 , N 1 ) [ 10 , 50 ] × [ 10 , 50 ] , M = 4 M 1 collocation points and N = 4 N 1 source points.
Symmetry 16 00134 g011
Figure 12. Example 4. The analytical f ( t ) ( ) and its approximation f * ( t ) with M 1 = N 1 = 26 , L = 52 , d 1 = d 2 = 0.8 and various noise added into the measured data, namely δ = 1 % ( ) ,   δ = 3 % ( ) and δ = 5 % ( ) .
Figure 12. Example 4. The analytical f ( t ) ( ) and its approximation f * ( t ) with M 1 = N 1 = 26 , L = 52 , d 1 = d 2 = 0.8 and various noise added into the measured data, namely δ = 1 % ( ) ,   δ = 3 % ( ) and δ = 5 % ( ) .
Symmetry 16 00134 g012
Figure 13. Example 4. The analytical f ( t ) ( ) and its approximation f * ( t ) with δ = 1 % , d 1 = d 2 = 0.8 and variable M N L , namely M = N = 10 , L = 20 ( ) , M = N = 20 , L = 40 ( ) , M = N = 30 , L = 60 ( ) , M = N = 40 , L = 80 ( ) , M = N = 50 , L = 100 ( ) .
Figure 13. Example 4. The analytical f ( t ) ( ) and its approximation f * ( t ) with δ = 1 % , d 1 = d 2 = 0.8 and variable M N L , namely M = N = 10 , L = 20 ( ) , M = N = 20 , L = 40 ( ) , M = N = 30 , L = 60 ( ) , M = N = 40 , L = 80 ( ) , M = N = 50 , L = 100 ( ) .
Symmetry 16 00134 g013
Figure 14. Example 5. In the exact solution ( ) and numerical solution ( ) for source f ( t ) ; we use 400 collocation points and 300 source points, and d ρ = 3 , d h = 5 , ρ 1 = 0.1 , h 1 = 0.5 , and noise level δ = 1 % with μ = 2.18 × 10 5 , and RMSE is 0.0745.
Figure 14. Example 5. In the exact solution ( ) and numerical solution ( ) for source f ( t ) ; we use 400 collocation points and 300 source points, and d ρ = 3 , d h = 5 , ρ 1 = 0.1 , h 1 = 0.5 , and noise level δ = 1 % with μ = 2.18 × 10 5 , and RMSE is 0.0745.
Symmetry 16 00134 g014
Figure 15. Example 5. The RES for ( ρ 1 , h 1 ) ( 0.1 , 0.9 ) × ( 0.1 , 0.9 ) , d ρ = 3 , d h = 5 , M = 400 , N = 300 and δ = 1 % .
Figure 15. Example 5. The RES for ( ρ 1 , h 1 ) ( 0.1 , 0.9 ) × ( 0.1 , 0.9 ) , d ρ = 3 , d h = 5 , M = 400 , N = 300 and δ = 1 % .
Symmetry 16 00134 g015
Figure 16. Example 5. The RES for ( d ρ , d h ) [ 1 , 5 ] × [ 1 , 5 ] , ρ 1 = 0.1 , h 1 = 0.5 , M = 400 ,   N = 300 and δ = 1 % .
Figure 16. Example 5. The RES for ( d ρ , d h ) [ 1 , 5 ] × [ 1 , 5 ] , ρ 1 = 0.1 , h 1 = 0.5 , M = 400 ,   N = 300 and δ = 1 % .
Symmetry 16 00134 g016
Figure 17. Example 6. In the exact solution ( ) and numerical solution ( ) for source f ( t ) ; we use 400 collocation points and 300 source points, and d ρ = 3 , d h = 4 , ρ 1 = 0.6 , h 1 = 0.3 , and noise level δ = 1 % with μ = 6.49 × 10 5 , and RMSE is 0.0109.
Figure 17. Example 6. In the exact solution ( ) and numerical solution ( ) for source f ( t ) ; we use 400 collocation points and 300 source points, and d ρ = 3 , d h = 4 , ρ 1 = 0.6 , h 1 = 0.3 , and noise level δ = 1 % with μ = 6.49 × 10 5 , and RMSE is 0.0109.
Symmetry 16 00134 g017
Figure 18. Example 6. The RES for ( ρ 1 , h 1 ) ( 0.1 , 0.9 ) × ( 0.1 , 0.9 ) , d ρ = 3 , d h = 4 , M = 400 , N = 300 and δ = 1 % .
Figure 18. Example 6. The RES for ( ρ 1 , h 1 ) ( 0.1 , 0.9 ) × ( 0.1 , 0.9 ) , d ρ = 3 , d h = 4 , M = 400 , N = 300 and δ = 1 % .
Symmetry 16 00134 g018
Figure 19. Example 6. The RES for ( d ρ , d h ) [ 1.5 , 4.5 ] × [ 1.5 , 4.5 ] , ρ 1 = 0.6 , h 1 = 0.3 ,   M = 400 , N = 300 and δ = 1 % .
Figure 19. Example 6. The RES for ( d ρ , d h ) [ 1.5 , 4.5 ] × [ 1.5 , 4.5 ] , ρ 1 = 0.6 , h 1 = 0.3 ,   M = 400 , N = 300 and δ = 1 % .
Symmetry 16 00134 g019
Table 1. The values of R E S ( f ) for various δ and ρ 1 , d 1 = d 2 = 2.5 .
Table 1. The values of R E S ( f ) for various δ and ρ 1 , d 1 = d 2 = 2.5 .
ρ 1 δ = 10 % δ = 5 % δ = 1 % δ = 0.1 %
0.10.0560 7.0959 × 10 3 3.2408 × 10 3 1.2344 × 10 3
0.20.0559 7.0575 × 10 3 2.4441 × 10 3 4.9399 × 10 4
0.30.0556 6.9314 × 10 3 1.3613 × 10 3 1.4585 × 10 4
0.40.0548 6.5549 × 10 3 8.9771 × 10 4 1.0601 × 10 4
0.50.0534 6.5024 × 10 3 8.6667 × 10 4 1.1133 × 10 4
0.60.0516 6.6017 × 10 3 9.8294 × 10 4 1.2614 × 10 4
0.70.0500 6.8182 × 10 3 1.0828 × 10 3 1.5702 × 10 4
0.80.0497 7.4956 × 10 3 1.2469 × 10 3 1.9143 × 10 4
0.90.0513 9.0394 × 10 3 1.4884 × 10 3 2.5788 × 10 4
Table 2. The values of R E S ( f ) for various of δ and ρ 1 , d 1 = d 2 = 2.5 .
Table 2. The values of R E S ( f ) for various of δ and ρ 1 , d 1 = d 2 = 2.5 .
ρ 1 δ = 10 % δ = 5 % δ = 1 % δ = 0.1 %
0.1 3.5361 × 10 2 4.4787 × 10 3 2.0455 × 10 3 7.7911 × 10 4
0.2 3.5308 × 10 2 4.4545 × 10 3 1.5426 × 10 3 3.1179 × 10 4
0.3 3.5094 × 10 2 4.3749 × 10 3 8.5922 × 10 4 9.2055 × 10 5
0.4 3.4602 × 10 2 4.1372 × 10 3 5.6661 × 10 4 6.6909 × 10 5
0.5 3.3723 × 10 2 4.1041 × 10 3 5.4701 × 10 4 7.0268 × 10 5
0.6 3.2567 × 10 2 4.1668 × 10 3 6.2040 × 10 4 7.9617 × 10 5
0.7 3.1582 × 10 2 4.3035 × 10 3 6.8345 × 10 4 9.9104 × 10 5
0.8 3.1394 × 10 2 4.7310 × 10 3 7.8702 × 10 4 1.2082 × 10 4
0.9 3.2409 × 10 2 5.7054 × 10 3 9.3945 × 10 4 1.6277 × 10 4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shen, Y.; Xiong, X. Identifying the Heat Source in Radially Symmetry and Axis-Symmetry Problems. Symmetry 2024, 16, 134. https://doi.org/10.3390/sym16020134

AMA Style

Shen Y, Xiong X. Identifying the Heat Source in Radially Symmetry and Axis-Symmetry Problems. Symmetry. 2024; 16(2):134. https://doi.org/10.3390/sym16020134

Chicago/Turabian Style

Shen, Yu, and Xiangtuan Xiong. 2024. "Identifying the Heat Source in Radially Symmetry and Axis-Symmetry Problems" Symmetry 16, no. 2: 134. https://doi.org/10.3390/sym16020134

APA Style

Shen, Y., & Xiong, X. (2024). Identifying the Heat Source in Radially Symmetry and Axis-Symmetry Problems. Symmetry, 16(2), 134. https://doi.org/10.3390/sym16020134

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop