Tikhonov regularization method for identifying the space-dependent source for time-fractional diffusion equation on a columnar symmetric domain

In this paper, we deal with the inverse problem of identifying the unknown source of time-fractional diffusion equation on a columnar symmetric domain. This problem is ill-posed. Firstly, we establish the conditional stability for this inverse problem. Then the regularization solution is obtained by using the Tikhonov regularization method and the error estimates are derived under the a priori and a posteriori choice rules of the regularization parameter. Three numerical examples are presented to illustrate the validity and effectiveness of our method.


Introduction
Diffusion equations with fractional order derivatives have been playing more and more important roles. For instance, they appear in mechanics, chemistry, electrical engineering and medicine [1][2][3][4][5][6][7][8][9][10][11]. Time-fractional diffusion equations can be used to describe some anomalous diffusion phenomena in many fields of science [12][13][14][15]. These fractional order models are more adequate than the integer order models, because the fractional order derivatives enable the description the properties of different substances [16].
In the past years, many regularization methods have been proposed to deal with the inverse problem for time-fractional diffusion equation, such as the backward problems [17][18][19][20][21][22], the inverse source problems [23][24][25][26][27], the Cauchy problem [28][29][30][31], the initial value problem [32,33]. In [34], the authors used the quasi-boundary value method to identify the initial value of heat equation on a columnar symmetric domain. In [35], the authors used a quasi-boundary regularization method to identify the initial value of timefractional diffusion equation on spherically symmetric domain. In this work, we will use the Tikhonov regularization method to identify the space-dependent source for the timefractional diffusion equation on a columnar symmetric domain. We not only give the a priori choice of the regularization parameter, but also give the a posteriori choice of the regularization parameter which only depends on the measurable date. To the best of our knowledge, there are few papers for the time-fractional diffusion equation on a columnar symmetric domain. In this work, we focus on an inverse problem for the following timefractional diffusion equation on a columnar symmetric domain: where r 0 is the radius of the cylinder, g(r) is given, f (r) is unknown source. D α t is the Caputo fractional derivative of order 0 < α < 1 defined by We use the final time to identify the unknown source f (r). In applications, the input function g(r) can be measured, and we assume the function g δ (r) as the measurable data, which satisfies where δ is a noise level of input data. The rest of this paper is organized as follows. In Sect. 2, we give some preliminaries. In Sect. 3, we analyze the ill-posedness of this problem and give the conditional stability. In Sect. 4, the error estimates are obtained under the a priori and a posteriori parameter choice rules. Three numerical examples are presented to demonstrate the effectiveness of our proposed method in Sect. 5.

Preliminaries
In this section, we will give some preliminaries which are very useful for our main conclusion.
If p ≥ 4, we can obtain Lemma 2.5 is proved.
If p ≥ 2, we can obtain Lemma 2.6 is proved.

Ill-posedness and a conditional stability
We can obtain the solution of the problem (1.1) by the method of separation of variables, as follows [39]: where f n = (f (r), ω n (r)), ω n (r) = √ 2 r 0 J 1 (μ n ) J 0 ( μ n r 0 ), n = 1, 2, 3, . . . , and {ω n (r)} ∞ n=1 is an orthonormal basis of L 2 [0, r 0 ], J 0 (·) and J 1 (·) are the zeroth order Bessel function and the first order Bessel function, respectively. μ n is the infinite number real root of the equation and it satisfies We consider the condition u(r, T) = g(r), then we have g n ω n (r), (3.5) where g n = (g(r), ω n (r)). Defining the operator K : f → g, we get The operator K is a linear self-adjoint compact operator, its singular value is {σ n } ∞ n=1 and and we also have Then we can obtain By using Lemma 2.4 and (3.7), we have we see that μ n → ∞ when n → ∞, so 1 σ n → ∞. So the exact data function g(r) must satisfy the property that (g, ω n (r)) decays rapidly. But we cannot ensure the function g(r) decreases, because the function g(r) concerns measurable data, a tiny disturbance of g(r) will cause a great error. So problem (1.1) is ill-posed. Assume for the unknown source f (r) there exists an a priori bound as follows: (3.12) where E > 0 is a constant and · H p denotes the norm in Hilbert space which is defined as follows [40]: The conditional stability of the inverse source problem can be obtained from Theorem 3.1.
Theorem 3.1 Let the a priori bound condition (3.12) hold, then we have Proof Due to (3.10), (3.11), (3.12), and the Hölder inequality, we obtain This completes the proof of Theorem 3.1.

Regularization method and convergence estimate
In this section, we will use the Tikhonov regularization method to obtain the regularization solution for problem (1.1). From Sect. 3, we know that {ω n (r)} ∞ n=1 is an orthogonal basis of L 2 [0, r 0 ; r], {σ n } ∞ n=1 is singular value of the linear self-adjoint compact operator K , and We adopt the Tikhonov regularization method to solve the ill-posed problem, which minimizes the following functional: where ν > 0 is the regularization parameter. From Theorem 2.12 of [41], we see that the minimum f ν (r) satisfies By a singular value decomposition of the compact self-adjoint operator, we have [42] f ν (r) = ∞ n=1 σ n σ 2 n + ν g n ω n (r).
Then we give the regularization solution with measurable data as follows:

An a priori parameter choice
In this subsection, we will give error estimates for under the suitable choice for the regularization parameter.
then the following error estimates hold: where C 1 , C 2 are positive constants depending on p, C, μ 1 .
Proof By utilizing triangle inequality, we have Using (4.3) and (1.3), we obtain Considering (3.10) and condition (3.12), we can obtain Using (3.11) and Lemma 2.5, we get We choose the regularization parameter as follows: then we have

An a posteriori selection rule
In this subsection, we will utilize Morozov's discrepancy principle to give an a posteriori regularization parameter choice. That is, we will choose the solution ν of the following equation as an a posterior regularization parameter: where τ > 1 is constant. We need Lemma 4.1 to obtain the existence and uniqueness of (4.5).
The proof of Lemma 4.1 is skipped.

Lemma 4.2
Let ν be the solution of (4.5), we have the following inequality: where C 3 , C 4 are positive constants depending on p, C, μ 1 .
Proof Using (4.5), we have Utilizing (3.11) and Lemma 2.6, we get So we can obtain Next, we will give the error estimate for under the a posteriori choice rule.
( 1 2 ( Proof Utilizing the triangle inequality, we have Using Lemma 4.1 and (4.4), we get For the second part of the right side of (4.6), using (1.3) and (4.5), for 0 < p < 1, we obtain We have Using Theorem 3.1, we get For p ≥ 2, because H p compacts into H 2 , then there exists a m ∈ N such that f (r) It is clear that ( 1 2 ( The proof of Theorem 4.2 is completed.
Then we take g δ as noise by adding a random perturbation, i.e., where ε reflects the relative error level.
Example 1 Take the source function f (r) = r sin(r).
Example 2 Consider the piecewise smooth function Example 3 Consider the discontinuous function (5.4) Figure 1 shows the comparisons between the exact solution and the regularization under the a priori and a posteriori regularization parameters choice when α = 0.2 for ε = 0.001, and ε = 0.0001 with Example 1. Figure 2 shows the comparisons between the exact solution and the regularization under the a priori and a posteriori regularization parameters  Figure 3 shows the comparisons between the exact solution and the regularization under the priori and posteriori regularization parameters choice when α = 0.2 for ε = 0.001, and ε = 0.0001 with Example 2. Figure 4 shows the comparisons between the exact solution and the regularization under the priori and posteriori regularization parameters choice when α = 0.6 for ε = 0.001, and ε = 0.0001 with Example 2. Figure 5 shows the comparisons between the exact solution and the regularization under the a priori and a posteriori regularization pa-  Figure 6 shows the comparisons between the exact solution and the regularization under the a priori and a posteriori regularization parameters choice when α = 0.6 for ε = 0.001, and ε = 0.0001 with Example 3. From Figs. 1-6, we can find the smaller ε, the better the computed approximation is. Moreover, we can also find the smaller α, the results are also better. Finally, we find that the results of Example 1 are better than that of Examples 2, 3, because in Examples 2, 3, the exact solutions are non-smooth and discontinuous functions, the recover data near the non-smooth and discontinuity points are not accurate. This is the well-known Gibbs phenomenon. But for the ill-posed problem, the results presented in Figs. 3-6 are reasonable.

Conclusion
In this paper, we use the Tikhonov regularization method to identify the source of the time-fractional diffusion equation on a columnar symmetric domain. Based on a conditional stability result, the error estimates are obtained under the a priori and a posteriori choice rules of regularization parameter. Meanwhile, the numerical examples verify the efficiency and accuracy of this method. Our original contributions are that we first identify the source of the time-fractional diffusion equation on a columnar symmetric domain. Moreover, we give the a posteriori regularization choice rule which only depends on the measurable data. In future, we will consider the inverse problem of identifying the initial value of the time-fractional diffusion equation on a columnar symmetric domain and give the optimal error estimate analyze. In addition, in this paper, we consider the timefractional derivative is the Caputo fractional derivative of order 0 < α < 1, but in Ref. [45][46][47][48][49][50][51][52], one mentioned the Caputo-Fabrizio fractional integro-differential equation, which is very useful in practice. We will consider the inverse problem of Caputo-Fabrizio fractional integro-differential equations, and use the Tikhonov regularization method to solve this inverse problem.