[go: up one dir, main page]
More Web Proxy on the site http://driver.im/ Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter February 8, 2023

Simultaneous Inversion of the Space-Dependent Source Term and the Initial Value in a Time-Fractional Diffusion Equation

  • Shuang Yu , Zewen Wang and Hongqi Yang EMAIL logo

Abstract

The inverse problem for simultaneously identifying the space-dependent source term and the initial value in a time-fractional diffusion equation is studied in this paper. The simultaneous inversion is formulated into a system of two operator equations based on the Fourier method to the time-fractional diffusion equation. Under some suitable assumptions, the conditional stability of simultaneous inversion solutions is established, and the exponential Tikhonov regularization method is proposed to obtain the good approximations of simultaneous inversion solutions. Then the convergence estimations of inversion solutions are presented for a priori and a posteriori selections of regularization parameters. Finally, numerical experiments are conducted to illustrate effectiveness of the proposed method.

MSC 2010: 35R30; 35R11; 47A52

Award Identifier / Grant number: 2020B1212060032

Award Identifier / Grant number: 11571386

Award Identifier / Grant number: 11961002

Award Identifier / Grant number: 12171248

Award Identifier / Grant number: 20212ACB201001

Funding statement: The work of S. Yu is supported in part by the NSF of China under grant 11961002. The work of S. Yu and H. Yang is supported in part by the Key-Area Research and Development Program of Guangdong Province (No. 2021B0101190003), by Guangdong Province Key Laboratory of Computational Science at the Sun Yat-sen University (2020B1212060032), by the NSF of China under grant 11571386. The work of Z. Wang is supported in part by the NSF of China under grants 11961002 and 12171248, and by Jiangxi Provincial Natural Science Foundation (20212ACB201001).

A Proofs of Theorems 1–3

Proof of Theorem 1

From (3.3) and the Hölder inequality, we have

f 2 = n = 1 + ( λ n E α , 1 ( λ n T 2 α ) E α , 1 ( λ n T 1 α ) ( E α , 1 ( λ n T 2 α ) g 1 , n E α , 1 ( λ n T 1 α ) g 2 , n ) ) 2 ( n = 1 + ( λ n E α , 1 ( λ n T 2 α ) E α , 1 ( λ n T 2 α ) E α , 1 ( λ n T 1 α ) ) γ + 2 ( g 1 , n E α , 1 ( λ n T 1 α ) E α , 1 ( λ n T 2 α ) g 2 , n ) 2 ) 2 γ + 2 × ( n = 1 + ( g 1 , n E α , 1 ( λ n T 1 α ) E α , 1 ( λ n T 2 α ) g 2 , n ) 2 ) γ γ + 2 ( n = 1 + 1 ( 1 E α , 1 ( λ n T 1 α ) E α , 1 ( λ n T 2 α ) ) γ λ n γ f n 2 ) 2 γ + 2 ( ( n = 1 + ( g 1 , n ) 2 ) 1 2 + ( n = 1 + ( E α , 1 ( λ n T 1 α ) E α , 1 ( λ n T 2 α ) g 2 , n ) 2 ) 1 2 ) 2 γ γ + 2 .

From Lemmas 13, we get

E α , 1 ( λ 1 T 1 α ) E α , 1 ( λ 1 T 2 α ) E α , 1 ( λ n T 1 α ) E α , 1 ( λ n T 2 α ) and E α , 1 ( λ n T 1 α ) E α , 1 ( λ n T 2 α ) lim n E α , 1 ( λ n T 1 α ) E α , 1 ( λ n T 2 α ) = ( T 2 T 1 ) α .

Obviously, λ n γ exp ( λ n γ ) holds. By the assumption of a priori condition (3.4) and the notation

C 3 := ( 1 E α , 1 ( λ 1 T 1 α ) E α , 1 ( λ 1 T 2 α ) 1 ) 2 γ γ + 2 ,

we have

f 2 C 3 M 4 γ + 2 ( g 1 + ( T 2 T 1 ) α g 2 ) 2 γ γ + 2 .

Similarly, we get

φ 2 = n = 1 + ( ( 1 E α , 1 ( λ n T 1 α ) ) g 2 , n ( 1 E α , 1 ( λ n T 2 α ) ) g 1 , n E α , 1 ( λ n T 2 α ) E α , 1 ( λ n T 1 α ) ) 2 ( n = 1 + ( 1 E α , 1 ( λ n T 1 α ) E α , 1 ( λ n T 2 α ) E α , 1 ( λ n T 1 α ) ) γ + 2 ( g 2 , n 1 E α , 1 ( λ n T 2 α ) 1 E α , 1 ( λ n T 1 α ) g 1 , n ) 2 ) 2 γ + 2 × ( n = 1 + ( g 2 , n 1 E α , 1 ( λ n T 2 α ) 1 E α , 1 ( λ n T 1 α ) g 1 , n ) 2 ) γ γ + 2 ( n = 1 + ( 1 E α , 1 ( λ n T 1 α ) E α , 1 ( λ n T 2 α ) E α , 1 ( λ n T 1 α ) ) γ φ n 2 ) 2 γ + 2 ( ( n = 1 + ( g 2 , n ) 2 ) 1 2 + ( n = 1 + ( 1 E α , 1 ( λ n T 2 α ) 1 E α , 1 ( λ n T 1 α ) g 1 , n ) 2 ) 1 2 ) 2 γ γ + 2 .

In addition, from Lemma 4, there exists a constant C > 0 such that

( n = 1 + ( 1 E α , 1 ( λ n T 1 α ) E α , 1 ( λ n T 2 α ) E α , 1 ( λ n T 1 α ) ) γ φ n 2 ) 2 γ + 2 ( n = 1 + ( Γ ( 1 α ) ( 1 + λ n T 2 α ) C 1 ( 1 E α , 1 ( λ n T 1 α ) E α , 1 ( λ n T 2 α ) ) ) γ φ n 2 ) 2 γ + 2 ( n = 1 + ( Γ ( 1 α ) ( C + T 2 α ) C 1 ( 1 E α , 1 ( λ n T 1 α ) E α , 1 ( λ n T 2 α ) ) ) γ λ n γ φ n 2 ) 2 γ + 2 .

Denoting

C 4 := ( Γ ( 1 α ) ( C + T 2 α ) C 1 ( E α , 1 ( λ 1 T 1 α ) E α , 1 ( λ 1 T 2 α ) 1 ) ) 2 γ γ + 2 ,

by the inequality λ n γ exp ( λ n γ ) , we have the estimate

φ 2 C 4 M 4 γ + 2 ( g 2 + 1 1 E α , 1 ( λ 1 T 1 α ) g 1 ) 2 γ γ + 2 .

The proof is completed. ∎

Proof of Theorem 2

By the triangular inequality, we have

(A.1) f β , γ δ f f β , γ δ f β , γ + f β , γ f ,
(A.2) φ β , γ δ φ φ β , γ δ φ β , γ + φ β , γ φ .

We first estimate f β , γ δ f . For the first term on the right-hand side of inequality (A.1), there exists a constant C > 0 such that

(A.3) f β , γ δ f β , γ 2 n = 1 + ( β exp ( λ n γ ) m 1 , n 1 β exp ( λ n γ ) ( m 1 , n 1 ) 2 + β 2 exp ( 2 λ n γ ) ) ( g 1 , n δ g 1 , n ) 2 + n = 1 + ( β exp ( λ n γ ) m 1 , n 2 β exp ( λ n γ ) ( m 1 , n 2 ) 2 + β 2 exp ( 2 λ n γ ) ) ( g 2 , n δ g 2 , n ) 2 sup n | λ n C 2 + β λ n γ + 2 | 2 ( n = 1 | g 1 , n δ g 1 , n | 2 + n = 1 | g 2 , n δ g 2 , n | 2 ) | ( 1 + γ ) 1 + γ 2 + γ ( 2 + γ ) C 2 ( 1 + γ ) 2 + γ β 1 2 + γ | 2 ( 2 δ 2 ) .

Denote

C 6 := 2 ( 1 + γ ) 1 + γ 2 + γ ( 2 + γ ) C 2 ( 1 + γ ) 2 + γ .

From inequality (A.3), we obtain

(A.4) f β , γ δ f β , γ C 6 δ β 1 2 + γ .

On the other hand, we get

exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 ) + ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 β exp ( λ n γ ) + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 m 1 , n 1 m 2 , n 1 + m 1 , n 2 m 2 , n 2 .

Lemma 3 yields

lim n + λ n 2 ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 m 1 , n 1 m 2 , n 1 + m 1 , n 2 m 2 , n 2 = ( T 1 α T 2 α ) 2 Γ ( 1 α ) T 1 α T 2 α ( T 1 α + T 2 α ) .

Therefore, there exist a constant

C 5 ( 0 , ( ( T 1 α T 2 α ) 2 Γ ( 1 α ) T 1 α T 2 α ( T 1 α + T 2 α ) ) 1 2 )

satisfying

(A.5) ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 m 1 , n 1 m 2 , n 1 + m 1 , n 2 m 2 , n 2 ( C 5 λ n ) 2 ,

with constant C 5 depending on parameters T 1 , T 2 and 𝛼. We now estimate the term f β , γ f . We have

f β , γ f 2 n = 1 + β exp ( λ n γ ) f n X n ( x ) β exp ( λ n γ ) + C 5 2 λ n 2 2 + n = 1 + β exp ( λ n γ ) φ n X n ( x ) β exp ( λ n γ ) + C 5 2 λ n 2 2 n = 1 + β 2 exp ( λ n γ ) ( β exp ( λ n γ ) + C 5 2 λ n 2 ) 2 exp ( λ n γ ) | f n | 2 + n = 1 + β 2 exp ( λ n γ ) ( β exp ( λ n γ ) + C 5 2 λ n 2 ) 2 exp ( λ n γ ) | φ n | 2 sup n ( β exp ( λ n γ 2 ) β exp ( λ n γ ) + C 5 2 λ n 2 ) 2 n = 1 + exp ( λ n γ ) | f n | 2 + sup n ( β exp ( λ n γ 2 ) β exp ( λ n γ ) + C 5 2 λ n 2 ) 2 n = 1 + exp ( λ n γ ) | φ n | 2 .

From Lemma 5 and 0 < β < 1 , we get

f β , γ f C 7 M β 1 4 , where C 7 = max { 3 3 4 4 ( 2 C 5 ) , λ N 0 2 exp ( λ N 0 2 2 ) ( 2 C 5 2 ) } .

This inequality, combined with estimate (A.4), yields the estimate

f β , γ δ f C 6 δ β 1 2 + γ + C 7 M β 1 4 .

By choosing the regularization parameter β = δ 8 + 4 γ 6 + γ , we get

f β , γ δ f ( C 6 + C 7 M ) δ 2 + γ 6 + γ .

Similarly, when 0 < β < 1 , γ > 0 , there exists C 8 , and selecting the regularization parameter β = δ 8 + 4 γ 6 + γ , inequality (A.2) yields

φ β , γ δ φ ( C 6 + C 8 M ) δ 2 + γ 6 + γ .

The proof is completed. ∎

Proof of Theorem 3

For τ δ defined in (4.7), there is the following inequality estimation:

τ δ = K 1 , 1 ( f β δ ) + K 2 , 1 ( φ β δ ) g 1 δ + K 1 , 2 ( f β δ ) + K 2 , 2 ( φ β δ ) g 2 δ
= n = 1 + β exp ( λ n γ ) ( ( m 1 , n 2 ) 2 + ( m 2 , n 2 ) 2 + β exp ( λ n γ ) ) g 1 , n δ + β exp ( λ n γ ) ( m 1 , n 2 m 1 , n 1 + m 2 , n 1 m 2 , n 2 ) g 2 , n δ ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) + β 2 exp ( 2 λ n γ ) X n ( x )
+ n = 1 + β exp ( λ n γ ) ( m 1 , n 1 m 1 , n 2 + m 2 , n 1 m 2 , n 2 ) g 1 , n δ β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 2 , n 1 ) 2 + β exp ( λ n γ ) ) g 2 , n δ ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) + β 2 exp ( 2 λ n γ ) X n ( x )
4 δ + i = 1 2 ( n = 1 + β exp ( λ n γ ) g i , n X n ( x ) β exp ( λ n γ ) + ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 + n = 1 + β exp ( λ n γ ) g i , n X n ( x ) β exp ( λ n γ ) + ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 m 1 , n 1 m 2 , n 1 + m 1 , n 2 m 2 , n 2 ) .
For C 5 defined in (A.5), we obtain

( τ 4 ) 2 δ 2 2 ( n = 1 + β exp ( λ n γ ) m 1 , n 1 f n X n ( x ) β exp ( λ n γ ) + ( C 5 λ n ) 2 2 + n = 1 + β exp ( λ n γ ) m 2 , n 1 φ n X n ( x ) β exp ( λ n γ ) + ( C 5 λ n ) 2 2 ) + 2 ( n = 1 + β exp ( λ n γ ) m 1 , n 2 f n X n ( x ) β exp ( λ n γ ) + ( C 5 λ n ) 2 2 + n = 1 + β exp ( λ n γ ) m 2 , n 2 φ n X n ( x ) β exp ( λ n γ ) + ( C 5 λ n ) 2 2 ) 4 ( ( sup n β exp ( 1 ) λ n 1 γ 2 β exp ( 0 ) λ n 2 + ( C 5 ) 2 ) 2 n = 1 + λ n γ | f n | 2 + ( sup n β exp ( 1 ) λ n 2 γ 2 β exp ( 0 ) λ n 2 + ( C 5 ) 2 ) 2 n = 1 + λ n γ | φ n | 2 ) .

Considering different ranges of 𝛾, we have the following estimates:

sup n β exp ( 1 ) λ n 1 γ 2 β λ n 2 + ( C 5 ) 2 { exp ( 1 ) ( 2 γ ) 2 γ 4 ( 2 + γ ) 2 + γ 4 4 ( C 5 ) 2 + γ 4 β 2 + γ 4 , 0 < γ < 2 , exp ( 1 ) λ 1 1 γ 2 ( C 5 ) 2 β , γ 2 ,
sup n β exp ( 1 ) λ n 2 γ 2 β λ n 2 + ( C 5 ) 2 { exp ( 1 ) γ γ 4 ( 4 γ ) 4 γ 4 4 ( C 5 ) γ 2 β γ 4 , 0 < γ < 4 , exp ( 1 ) λ 1 2 γ 2 ( C 5 ) 2 β , γ 4 .
From λ n γ exp ( λ n γ ) and denoting

E 1 = ( ( exp ( 1 ) ( 2 γ ) 2 γ 4 ( 2 + γ ) 2 + γ 4 2 ( C 5 ) 3 γ + 6 4 ) 2 + ( exp ( 1 ) γ γ 4 ( 4 γ ) 4 γ 4 2 ( C 5 ) γ ) 2 ) 1 2 , E 2 = ( ( 2 exp ( 1 ) λ 1 1 γ 2 ( C 5 ) 4 ) 2 + ( exp ( 1 ) γ γ 4 ( 4 γ ) 4 γ 4 2 ( C 5 ) γ ) 2 ) 1 2 , E 3 = ( ( 2 exp ( 1 ) λ 1 1 γ 2 ( C 5 ) 4 ) 2 + ( 2 exp ( 1 ) λ 1 2 γ 2 ( C 5 ) 4 ) 2 ) 1 2 ,

we get

( τ 4 ) δ { E 1 β 2 + γ 4 M , 0 < γ < 2 , E 2 β M , 2 γ < 4 , E 3 β M , γ 4 .

Therefore, denoting C 9 = max { E 1 , E 2 , E 3 } , we have the estimate

1 β 1 2 + γ { ( C 9 M τ 4 ) 4 ( 2 + γ ) 2 δ 4 ( 2 + γ ) 2 , 0 < γ < 2 , ( C 9 M τ 4 ) 1 2 + γ δ 1 2 + γ , γ 2 .

Then

f β , γ δ f β , γ { C 6 ( C 9 M τ 4 ) 4 ( 2 + γ ) 2 δ γ 2 + γ , 0 < γ < 2 , C 6 ( C 9 M τ 4 ) 1 2 + γ δ 1 + γ 2 + γ , γ 2 .

Similarly, the estimate for initial values is

φ β , γ δ φ β , γ { C 6 ( C 9 M τ 4 ) 4 ( 2 + γ ) 2 δ γ 2 + γ , 0 < γ < 2 , C 6 ( C 9 M τ 4 ) 1 2 + γ δ 1 + γ 2 + γ , γ 2 .

Furthermore,

f β , γ f Exp 2 n = 1 + ( β 2 exp ( 2 λ n γ ) + β exp ( λ n γ ) ( ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) + β 2 exp ( 2 λ n γ ) ) 2 exp ( λ n γ ) | f n | 2 + n = 1 + ( β exp ( λ n γ ) ( m 1 , n 1 m 2 , n 1 + m 1 , n 2 m 2 , n 2 ) ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) + β 2 exp ( 2 λ n γ ) ) 2 exp ( λ n γ ) | φ n | 2 n = 1 + exp ( λ n γ ) | f n | 2 + n = 1 + exp ( λ n γ ) | φ n | 2 2 ( M ) 2 .

Similarly, we have the same result on φ β , γ φ Exp 2 . On the other hand, let

F i ( n ) = m 1 , n i m 2 , n i m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 , i = 1 , 2 .

From noise assumptions (1.2), there exists C 10 such that

K 1 ( f β , γ f , φ β , γ φ ) 4 C 10 δ + n = 1 + β exp ( λ n γ ) ( ( m 1 , n 2 ) 2 + ( m 2 , n 2 ) 2 + β exp ( λ n γ ) ) g 1 , n δ + β exp ( λ n γ ) ( m 1 , n 2 m 1 , n 1 + m 2 , n 1 m 2 , n 2 ) g 2 , n δ X n ( x ) ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) + β 2 exp ( 2 λ n γ ) lim n | F 1 ( n ) m 2 , n 2 m 2 , n 1 | + n = 1 + β exp ( λ n γ ) ( m 1 , n 1 m 1 , n 2 + m 2 , n 1 m 2 , n 2 ) g 1 , n δ β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 2 , n 1 ) 2 + β exp ( λ n γ ) ) g 2 , n δ X n ( x ) ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) + β 2 exp ( 2 λ n γ ) lim n | F 1 ( n ) | + n = 1 + β exp ( λ n γ ) ( ( m 1 , n 2 ) 2 + ( m 2 , n 2 ) 2 + β exp ( λ n γ ) ) g 1 , n δ β exp ( λ n γ ) ( m 1 , n 2 m 1 , n 1 + m 2 , n 1 m 2 , n 2 ) g 2 , n δ X n ( x ) ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) + β 2 exp ( 2 λ n γ ) lim n | F 2 ( n ) m 2 , n 1 m 2 , n 2 | + n = 1 + β exp ( λ n γ ) ( m 1 , n 1 m 1 , n 2 + m 2 , n 1 m 2 , n 2 ) g 1 , n δ + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 2 , n 1 ) 2 + β exp ( λ n γ ) ) g 2 , n δ X n ( x ) ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) + β 2 exp ( 2 λ n γ ) lim n | F 1 ( n ) | .

From Lemma 3, we get

lim n | F 1 ( n ) | = lim n | F 1 ( n ) m 2 , n 2 m 2 , n 1 | = lim n | F 1 ( n ) m 2 , n 1 m 2 , n 2 | T 1 α T 2 α Γ ( 1 α ) T 1 α T 2 α .

Therefore, by the Morozov discrepancy principle (4.7), we get

K 1 ( f β , γ f , φ β , γ φ ) ( 4 C 10 + 2 τ T 1 α T 2 α Γ ( 1 α ) T 1 α T 2 α ) δ .

Similarly, for the operator K 2 , there exists C 11 such that

K 2 ( f β , γ f , φ β , γ φ ) 4 C 11 δ + n = 1 + β exp ( λ n γ ) ( ( m 1 , n 2 ) 2 + ( m 2 , n 2 ) 2 + β exp ( λ n γ ) ) g 1 , n δ + β exp ( λ n γ ) ( m 1 , n 2 m 1 , n 1 + m 2 , n 1 m 2 , n 2 ) g 2 , n δ X n ( x ) ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) + β 2 exp ( 2 λ n γ ) lim n | F 2 ( n ) | + n = 1 + β exp ( λ n γ ) ( m 1 , n 1 m 1 , n 2 + m 2 , n 1 m 2 , n 2 ) g 1 , n δ β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 2 , n 1 ) 2 + β exp ( λ n γ ) ) g 2 , n δ X n ( x ) ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) + β 2 exp ( 2 λ n γ ) lim n | F 2 ( n ) m 1 , n 1 m 1 , n 2 | + n = 1 + β exp ( λ n γ ) ( ( m 1 , n 2 ) 2 + ( m 2 , n 2 ) 2 + β exp ( λ n γ ) ) g 1 , n δ β exp ( λ n γ ) ( m 1 , n 2 m 1 , n 1 + m 2 , n 1 m 2 , n 2 ) g 2 , n δ X n ( x ) ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) + β 2 exp ( 2 λ n γ ) lim n | F 2 ( n ) | + n = 1 + β exp ( λ n γ ) ( m 1 , n 1 m 1 , n 2 + m 2 , n 1 m 2 , n 2 ) g 1 , n δ + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 2 , n 1 ) 2 + β exp ( λ n γ ) ) g 2 , n δ X n ( x ) ( m 1 , n 1 m 2 , n 2 m 1 , n 2 m 2 , n 1 ) 2 + β exp ( λ n γ ) ( ( m 1 , n 1 ) 2 + ( m 1 , n 2 ) 2 + ( m 2 , n 1 ) 2 + ( m 2 , n 2 ) 2 ) + β 2 exp ( 2 λ n γ ) lim n | F 2 ( n ) m 2 , n 1 m 2 , n 2 | .

Therefore, we have

K 2 ( f β , γ f , φ β , γ φ ) ( 4 C 11 + 2 τ T 1 α T 2 α Γ ( 1 α ) T 1 α T 2 α ) δ .

From inequality (3.5) in Theorem 1, we get

f β , γ f ( C 1 ) γ γ + 2 f β , γ f Exp 2 γ + 2 ( K 1 ( f β , γ f , φ β , γ φ ) + ( T 2 T 1 ) α K 2 ( f β , γ f , φ β , γ φ ) ) γ γ + 2 C 12 M 1 γ + 2 δ γ γ + 2 ,

where

C 12 = ( 2 1 γ C 1 ( 4 ( C 10 + T 2 α T 1 α C 11 ) + 2 τ Γ ( 1 α ) ( T 1 α + T 2 α ) T 2 α T 1 α T 2 α ) ) γ γ + 2 .

And (3.6) in Theorem 1 yields

φ β , γ φ ( C 2 ) γ γ + 2 φ β , γ φ Exp 2 γ + 2 ( K 2 ( f β , γ f , φ β , γ φ ) + 1 1 E α , 1 ( λ 1 T 1 α ) K 1 ( f β , γ f , φ β , γ φ ) ) γ γ + 2 C 13 M 1 γ + 2 δ γ γ + 2 ,

where

C 13 = ( 2 1 γ C 2 ( 4 ( C 11 + 1 1 E α , 1 ( λ 1 T 1 α ) C 10 ) + 2 E α , 1 ( λ 1 T 1 α ) 1 E α , 1 ( λ 1 T 1 α ) 2 τ T 1 α T 2 α Γ ( 1 α ) T 1 α T 2 α ) ) γ γ + 2 .

Therefore, denoting

C 14 = C 6 ( C 9 M τ 4 ) 4 ( 2 + γ ) 2 and C 15 = C 6 ( C 9 M τ 4 ) 1 2 + γ ,

we have the conclusions

f β , γ δ f { ( C 14 + C 12 M 1 γ + 2 ) δ γ 2 + γ , 0 < γ < 2 , ( C 15 + C 12 M 1 γ + 2 ) δ 1 + γ 2 + γ , γ 2 , φ β , γ δ φ { ( C 14 + C 13 M 1 γ + 2 ) δ γ 2 + γ , 0 < γ < 2 , ( C 15 + C 13 M 1 γ + 2 ) δ 1 + γ 2 + γ , γ 2 .

The proof is completed. ∎

References

[1] W. Audeh and F. Kittaneh, Singular value inequalities for compact operators, Linear Algebra Appl. 437 (2012), no. 10, 2516–2522. 10.1016/j.laa.2012.06.032Search in Google Scholar

[2] B. Berkowitz, H. Scher and S. E. Silliman, Anomalous transport in laboratory-scale, heterogeneous porous media, Water Resources Res. 36 (2000), no. 1, 149–158. 10.1029/1999WR900295Search in Google Scholar

[3] J. Cheng, J. Nakagawa, M. Yamamoto and T. Yamazaki, Uniqueness in an inverse problem for a one-dimensional fractional diffusion equation, Inverse Problems 25 (2009), no. 11, Article ID 115002. 10.1088/0266-5611/25/11/115002Search in Google Scholar

[4] M. Giona, S. Cerbelli and H. E. Roman, Fractional diffusion equation and relaxation in complex viscoelastic materials, Phys. A 191 (1992), no. 1–4, 449–453. 10.1016/0378-4371(92)90566-9Search in Google Scholar

[5] G. Floridia, Z. Li and M. Yamamoto, Well-posedness for the backward problems in time for general time-fractional diffusion equation, Atti Accad. Naz. Lincei Rend. Lincei Mat. Appl. 31 (2020), no. 3, 593–610. 10.4171/RLM/906Search in Google Scholar

[6] D. Hou, M. T. Hasan and C. Xu, Müntz spectral methods for the time-fractional diffusion equation, Comput. Methods Appl. Math. 18 (2018), no. 1, 43–62. 10.1515/cmam-2017-0027Search in Google Scholar

[7] D. Jiang, Z. Li, Y. Liu and M. Yamamoto, Weak unique continuation property and a related inverse source problem for time-fractional diffusion-advection equations, Inverse Problems 33 (2017), no. 5, Article ID 055013. 10.1088/1361-6420/aa58d1Search in Google Scholar

[8] S. Jiang, K. Liao and T. Wei, Inversion of the initial value for a time-fractional diffusion-wave equation by boundary data, Comput. Methods Appl. Math. 20 (2020), no. 1, 109–120. 10.1515/cmam-2018-0194Search in Google Scholar

[9] G. Li, D. Zhang, X. Jia and M. Yamamoto, Simultaneous inversion for the space-dependent diffusion coefficient and the fractional order in the time-fractional diffusion equation, Inverse Problems 29 (2013), no. 6, Article ID 065014. 10.1088/0266-5611/29/6/065014Search in Google Scholar

[10] Z. Li, Y. Liu and M. Yamamoto, Initial-boundary value problems for multi-term time-fractional diffusion equations with positive constant coefficients, Appl. Math. Comput. 257 (2015), 381–397. 10.1016/j.amc.2014.11.073Search in Google Scholar

[11] Z. Li and M. Yamamoto, Initial-boundary value problems for linear diffusion equation with multiple time-fractional derivatives, preprint (2013), https://arxiv.org/abs/1306.2778. Search in Google Scholar

[12] K. Liao and T. Wei, Identifying a fractional order and a space source term in a time-fractional diffusion-wave equation simultaneously, Inverse Problems 35 (2019), no. 11, Article ID 115002. 10.1088/1361-6420/ab383fSearch in Google Scholar

[13] J. J. Liu and M. Yamamoto, A backward problem for the time-fractional diffusion equation, Appl. Anal. 89 (2010), no. 11, 1769–1788. 10.1080/00036810903479731Search in Google Scholar

[14] Y. Liu, Z. Li and M. Yamamoto, Inverse problems of determining sources of the fractional partial differential equations, Handbook of Fractional Calculus with Applications. Vol. 2, De Gruyter, Berlin (2019), 411–429. 10.1515/9783110571660-018Search in Google Scholar

[15] Y. Luchko, Some uniqueness and existence results for the initial-boundary-value problems for the generalized time-fractional diffusion equation, Comput. Math. Appl. 59 (2010), no. 5, 1766–1772. 10.1016/j.camwa.2009.08.015Search in Google Scholar

[16] Y. Luchko, Initial-boundary problems for the generalized multi-term time-fractional diffusion equation, J. Math. Anal. Appl. 374 (2011), no. 2, 538–548. 10.1016/j.jmaa.2010.08.048Search in Google Scholar

[17] R. Metzler and J. Klafter, Subdiffusive transport close to thermal equilibrium: From the Langevin equation to fractional diffusion, Phys. Rev. E 61 (2000), no. 6, 6308–11. 10.1103/PhysRevE.61.6308Search in Google Scholar

[18] R. Metzler, J. Klafter and I. M. Sokolov, Anomalous transport in external fields: Continuous time random walks and fractional diffusion equations extended, Phys. Rev. E 58 (1998), no. 2, 1621–33. 10.1103/PhysRevE.58.1621Search in Google Scholar

[19] V. A. Morozov, Methods for Solving Incorrectly Posed Problems, Springer, New York, 2012. Search in Google Scholar

[20] R. R. Nigmatullin, The realization of the generalized transfer equation in a medium with fractal geometry, Phys. Status Solidi (b) 133 (1986), no. 1, 425–430. 10.1002/pssb.2221330150Search in Google Scholar

[21] I. Podlubny, Fractional Differential Equations, Math. Sci. Eng. 198, Academic Press, San Diego, 1999. Search in Google Scholar

[22] H. Pollard, The completely monotonic character of the Mittag-Leffler function E a ( x ) , Bull. Amer. Math. Soc. 54 (1948), 1115–1116. 10.1090/S0002-9904-1948-09132-7Search in Google Scholar

[23] Z. Ruan, J. Z. Yang and X. Lu, Tikhonov regularisation method for simultaneous inversion of the source term and initial data in a time-fractional diffusion equation, East Asian J. Appl. Math. 5 (2015), no. 3, 273–300. 10.4208/eajam.310315.030715aSearch in Google Scholar

[24] Z. Ruan and S. Zhang, Simultaneous inversion of time-dependent source term and fractional order for a time-fractional diffusion equation, J. Comput. Appl. Math. 368 (2020), Article ID 112566. 10.1016/j.cam.2019.112566Search in Google Scholar

[25] Z. Ruan, W. Zhang and Z. Wang, Simultaneous inversion of the fractional order and the space-dependent source term for the time-fractional diffusion equation, Appl. Math. Comput. 328 (2018), 365–379. 10.1016/j.amc.2018.01.025Search in Google Scholar

[26] K. Sakamoto and M. Yamamoto, Initial value/boundary value problems for fractional diffusion-wave equations and applications to some inverse problems, J. Math. Anal. Appl. 382 (2011), no. 1, 426–447. 10.1016/j.jmaa.2011.04.058Search in Google Scholar

[27] R. Schumer, D. A. Benson and M. M. Meerschaert, Fractal mobile/immobile solute transport, Water Resources Res. 39 (2003), 10.1029/2003WR002141. 10.1029/2003WR002141Search in Google Scholar

[28] I. M. Sokolov and J. Klafter, From diffusion to anomalous diffusion: A century after Einstein’s Brownian motion, Chaos 15 (2005), no. 2, Article ID 026103. 10.1063/1.1860472Search in Google Scholar PubMed

[29] Z. Wang, S. Qiu, S. Yu, B. Wu and W. Zhang, Exponential Tikhonov regularization method for solving an inverse source problem of time fractional diffusion equation, J. Comput. Math. 41 (2023), no. 2, 173–190. 10.4208/jcm.2107-m2020-0133Search in Google Scholar

[30] T. Wei, X. L. Li and Y. S. Li, An inverse time-dependent source problem for a time-fractional diffusion equation, Inverse Problems 32 (2016), no. 8, Article ID 085003. 10.1088/0266-5611/32/8/085003Search in Google Scholar

[31] T. Wei and Y. Luo, A generalized quasi-boundary value method for recovering a source in a fractional diffusion-wave equation, Inverse Problems 38 (2022), no. 4, Paper No. 045001. 10.1088/1361-6420/ac50b9Search in Google Scholar

[32] J. Wen, Z.-X. Liu, C.-W. Yue and S.-J. Wang, Landweber iteration method for simultaneous inversion of the source term and initial data in a time-fractional diffusion equation, J. Appl. Math. Comput. 68 (2022), no. 5, 3219–3250. 10.1007/s12190-021-01656-0Search in Google Scholar

[33] W. Wyss, The fractional diffusion equation, J. Math. Phys. 27 (1986), no. 11, 2782–2785. 10.1063/1.527251Search in Google Scholar

[34] X. Xiong and X. Xue, A fractional Tikhonov regularization method for identifying a space-dependent source in the time-fractional diffusion equation, Appl. Math. Comput. 349 (2019), 292–303. 10.1016/j.amc.2018.12.063Search in Google Scholar

[35] X.-B. Yan, Z.-Q. Zhang and T. Wei, Simultaneous inversion of a time-dependent potential coefficient and a time source term in a time fractional diffusion-wave equation, Chaos Solitons Fractals 157 (2022), Paper No. 111901. 10.1016/j.chaos.2022.111901Search in Google Scholar

[36] Y. Zhang, T. Wei and Y.-X. Zhang, Simultaneous inversion of two initial values for a time-fractional diffusion-wave equation, Numer. Methods Partial Differential Equations 37 (2021), no. 1, 24–43. 10.1002/num.22517Search in Google Scholar

Received: 2022-03-10
Revised: 2022-12-05
Accepted: 2022-12-28
Published Online: 2023-02-08
Published in Print: 2023-07-01

© 2023 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 21.1.2025 from https://www.degruyter.com/document/doi/10.1515/cmam-2022-0058/html
Scroll to top button