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

: 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.


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: ( 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 gridbased 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 oilgas-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.

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

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.
with the extra condition where ρ 1 ∈ (ρ 0 , P),u 0 (ρ),d 1 (t),d 2 (t) and s(t) are given functions satisfying the compatibility conditions 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).
The Formulas ( 3)-( 7) are converted to the following formula: Substituting Equation ( 15) into Equations ( 11)-( 14) yields the following problem: When v is obtained, the unknown heat source f (t) can be represented by the following formula:

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.
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 where u(ρ 1 , h 1 , 0) = g(0) for ρ 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 Then, the Formulas ( 21)-( 24) can be converted to a homogeneous equation with non-local boundary conditions: By subtracting ( 30) from ( 29), it can be observed that the new function v satisfies the following problem: When v is obtained, the unknown heat source f (t) can be represented by the following formula:

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.

The Methodology for One-Dimensional Problem
According to literature [35], the fundamental solution of equation Then, the fundamental solution of Equation ( 11) is expressed as 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 We approximate the solution of Equation ( 11) by using the following linear combination of fundamental solutions. where The source point locations to be placed in time are chosen as i and v i for i = 1, • • • , 2N 1 will be determined by the boundary conditions ( 17)- (19).
We choose the collocation points: and It can be obtained from the above: of M = (2M 1 + L) linear equations with N = 4N 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: 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 i F 1 (ρ 0 , t; P Therefore, the unknown source term is approximated as According to (37), we need to calculate the first derivative of the fundamental solution (36) with respect to t. Simple calculations show that

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 twodimensional 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 The solution of the axisymmetric heat Equation ( 27) can be approximated by a linear combination of the fundamental solutions.
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 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 There are two cases: For the sake of simplicity, we assume there exists q > 0 such that N 1 = qM 1 or M 1 = qN 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 and on the boundary Γ (see Equation ( 32)), let On the fictitious boundary Γ E , set where 2 .Boundary conditions are then imposed at the collocation points to determine the constant coefficients v (j) i , as given in (43).We have 2 ) source points.To obtain a unique solution, we restrict N ≤ M. Now, we can form a linear system.
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 and

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 gi = g i (1 + δ • rand(i)), where g i is the exact data, and rand(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 49) 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.

Numerical Examples for MFS: Radially Symmetric Heat Equation
Example 1.The exact solution to the problems ( 3)-( 7) is given by: 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 RES( 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.Figures 3-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.Figures 6-8 of Example 2 is similar to Example 1, and we will not go into too much detail here.Example 3. We consider another problem given by u(ρ, t) From Table 2, similar conclusions can be observed to those from Table 1.As can be seen from Figures 9-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: with the heat source function   We solve the problem with MFS, with 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.

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 collocation points and 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 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 + (ρ

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.

Figure 17 illustrates
Figure17illustrates 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.Figures18 and 19respectively 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.