Next Article in Journal
Face Recognition via Compact Second-Order Image Gradient Orientations
Next Article in Special Issue
Fractional Integrals and Derivatives: “True” versus “False”
Previous Article in Journal
Deep Adversarial Learning Triplet Similarity Preserving Cross-Modal Retrieval Algorithm
Previous Article in Special Issue
A New Look at the Initial Condition Problem
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Inverse Problem for a Non-Homogeneous Time-Space Fractional Equation

1
Laboratoire LaSIE, La Rochelle Université, 17042 La Rochelle, France
2
Department of Mathematics, Khalifa University of Science and Technology, Abu Dhabi P.O. Box 127788, United Arab Emirates
3
NAAM Research Group, Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(15), 2586; https://doi.org/10.3390/math10152586
Submission received: 24 April 2022 / Revised: 14 July 2022 / Accepted: 18 July 2022 / Published: 25 July 2022

Abstract

:
We consider the inverse problem of finding the solution of a generalized time-space fractional equation and the source term knowing the spatial mean of the solution at any times t ( 0 , T ] , as well as the initial and the boundary conditions. The existence and the continuity with respect to the data of the solution for the direct and the inverse problem are proven by Fourier’s method and the Schauder fixed-point theorem in an adequate convex bounded subset. In the published articles on this topic, the incorrect use of the estimates in the generalized Mittag–Leffler functions is commonly performed. This leads to false proofs of the Fourier series’ convergence to recover the equation satisfied by the solution, the initial data or the boundary conditions. In the present work, the correct framework to recover the decay of fractional Fourier coefficients is established; this allows one to recover correctly the initial data, the boundary conditions and the partial differential equations within the space-time domain.

1. Introduction

Consider the following non-homogeneous fractional evolution boundary value problem:
D t γ u ( x , t ) i = 1 d D x i α i u ( x , t ) = p ( x ) f ( t ) , ( x , t ) Ω × ( 0 , T ) , u ( x , t ) = 0 , ( x , t ) Ω × [ 0 , T ] , u ( x , 0 ) = φ ( x ) , x Ω ¯ ,
where Ω = ( 0 , 1 ) d , d 1 , x = ( x 1 , , x d ) , D t γ u ( x , t ) denotes the (time) fractional derivative in the sense of Caputo, with 0 < γ < 1 , and D x i α i u ( x , t ) ; 1 i d , is the (space) fractional derivative, in the sense of Riemann–Liouville, with 1 < α i < 2 ; and T > 0 denotes the length of the time interval.
Let us fix the framework: When the initial data φ and the source terms ( p , f ) are given, the direct problem consists of looking for the unknown function u in an adequate functional space. On the other hand, when the initial data φ and only the spatial part p are given, the inverse problem consists of looking for the unknown functions ( u , f ) via, for example, the supplementary data:
Ω u ( x , t ) d x = ψ ( t ) , t [ 0 , T ] ,
where ψ is a given function of [ 0 , T ] .
Notice that the fractional operator D x i α i is unbounded, non-self-adjoint with compact resolvent admitting a countable complete basis of eigenfunctions, in L 2 ( Ω , C ) [1]. This will allow us to look for classical solutions via the biorthogonal representation with respect to x. Indeed, the uniform convergence of the resulting series in Ω ¯ × [ 0 , T ] (and of their derivatives) and the uniqueness of the corresponding limits, with respect to the norms · 2 and · , will lead to the existence of classical solutions of the direct problem. For the inverse problems, we will make use of the Schauder fixed point theorem in an adequate convex and bounded subset.
Before we present and prove our results, let us dwell on the existing literature. An analytical and numerical study of inverse problems for fractional diffusion equations has undergone intensive development over the present decade. Kirane and Malik [2] considered at some later date a one-dimensional fractional diffusion equation with a non-local, non-self-adjoint boundary condition. They developed a solution and source function using a bi-orthogonal method and obtained the uniqueness. In addition, they have prove an existence result when both the initial and the final conditions are smooth enough. In [3], the problem in [2] is generalized. Furati et al. used the same bi-orthogonal spatial system to obtain a separable formal solution and source function and to show their uniqueness. They then used the asymptotic expansion of the generalized Mittag–Leffler function to achieve a result of existence under certain smoothness criteria on the initial and final conditions. Aleroev et al. [1] consider a linear heat equation, with a non-local boundary condition, involving a fractional derivative in time. They establish a space independent source term and the distribution of temperature for a problem with integral type over determining condition. They proved the solution’s existence and uniqueness, and its continuous dependence on the data. Among other papers devoted to determining unknown source terms, we may cite [4,5,6].
As was mentioned in the abstract, an incorrect use of the estimates in the generalized Mittag–Leffler functions is commonly performed in the literature on this topic. This leads to false proofs of the Fourier series’ convergence to recover the equation satisfied by the solution, the initial data or the boundary conditions. Remark 1 explains this fairly widespread mistake. In the present work, the correct framework to recover the decay of fractional Fourier coefficients is established; this allows one to correctly recover the initial data, the boundary conditions and the partial differential equation satisfied within the space-time domain.
The manuscript is organized as follows: In Section 2, fundamental concepts of fractional calculus are recalled. In Section 3, we obtain results for the direct problem using Fourier’s method; we determine the existence, uniqueness and continuous dependence on the data. In Section 4, the inverse problem associated with the direct problem is considered, and we show the results of the existence, uniqueness and continuous dependence of the solution on the data.

2. Preliminaries

In this section, we recall basic definitions and notations from fractional calculus [7,8]. Let f : [ a , b ] R be an integrable function. The left-sided Riemann–Liouville fractional integral of order α > 0 is defined by:
I x α f ( x ) = 1 Γ ( α ) a x f ( t ) ( x t ) 1 α d t ,
where Γ ( α ) is the Euler’s Gamma function.
For arbitrary values of α > 0 , the left sided Riemann–Liouville fractional derivative of order α is defined by:
D x α f ( x ) = d n d x n I x n α f ( x ) = 1 Γ ( n α ) d n d t n a x ( x t ) n 1 α f ( t ) d t , n = α + 1 ,
where α denotes the integer part of α . The left-sided Caputo fractional derivative of order α > 0 is defined by:
D x α f ( x ) = I x n α f ( n ) ( x ) = 1 Γ ( n α ) a x ( x t ) n 1 α f ( n ) ( t ) d t , n = α + 1 .
where f ( n ) denotes the n-th classical derivative of f.
The Mittag–Leffler function with two parameters α and β is defined as:
E α , β ( z ) = k = 0 z k Γ ( α k + β ) , for any z C .
In the case β = 1 , the function E α , 1 is the classical Mittag–Leffler function denoted by E α . From [9], we have the following property:
0 z t β 1 E α , β ( λ t α ) d t = z β E α , β + 1 ( λ z α ) , ( α > 0 , β > 0 , λ R ) .
We also need the following result [10].
Theorem 1.
Let θ ( 0 , 2 ) , ν > 0 and μ be such that θ π 2 < μ < min ( π , θ π ) . Then, there is a constant c = c ( θ , ν , μ ) 0 such that:
E θ , ν ( z ) c 1 + | z | , z C , μ | arg ( z ) | π
and
E θ , ν ( z ) c 1 1 + | z | + ( 1 + | z | ) 1 ν θ e R e z 1 θ , z C , | arg ( z ) | μ ,
where R e ( z ) denotes the real part of the complex number z.
Furthermore, from [7], we have:
Lemma 1.
Let ( f i ) i N be a sequence of functions defined on the interval ( a , b ] . Suppose the following conditions are fulfilled:
  • For a given α > 0 , the fractional derivative D z α f i exists for all i N and z ( a , b ] ;
  • Both series i = 1 f i ( z ) and i = 1 D z α f i ( z ) are uniformly convergent on the interval [ a + ε , b ] for any ε > 0 .
Then, the function defined by the series i = 1 f i ( z ) is α-differentiable on ( a , b ] and satisfies:
D z α i = 1 f i ( z ) = i = 1 D z α f i ( z ) .

3. Direct Problem

In this section, we will study the forward problem:
D t γ u ( x , t ) i = 1 d D x i α i u ( x , t ) = p ( x ) f ( t ) , ( x , t ) Ω × ( 0 , T ) , u ( x , t ) = 0 , ( x , t ) Ω × [ 0 , T ] , u ( x , 0 ) = φ ( x ) , x Ω ¯ ,
where p, f and φ are given functions such that p , φ C Ω ¯ , R and f C [ 0 , T ] , R .
We are interested in classical solutions to Problem (10), that is, solutions u ( x , t ) which are continuous on Ω ¯ × [ 0 , T ] and satisfy D t γ u ( x , t ) L 2 ( Ω × ( 0 , T ) ) and D x i α i u ( x , t ) L 2 ( Ω × ( 0 , T ) ) , for any ( x , t ) Ω × ( 0 , T ) .
Example 1.
Consider the case d = 2 and Ω = ( 0 , π ) × ( 0 , π ) , with 0 < γ < 1 and 1 < α 1 , α 2 < 2 . We fix the initial condition φ ( x 1 , x 2 ) , the spatial and temporal source terms p ( x 1 , x 2 ) and f ( t ) for any ( x 1 , x 2 ) Ω and t > 0 as the following:
φ ( x 1 , x 2 ) = 0 , p ( x 1 , x 2 ) = sin ( x 1 ) sin ( x 2 ) + x 1 1 α 1 E 1 2 , 2 α 1 ( x 1 2 ) sin ( y ) + sin ( x ) x 2 1 α 2 E 1 2 , 2 α 2 ( x 2 2 ) , f ( t ) = E γ , 1 t γ .
The exact solution is then u ( x 1 , x 2 , t ) = sin ( x 1 ) sin ( x 2 ) E γ , 1 t γ , since it is known that
D t γ E γ , 1 t γ = E γ , 1 t γ , t > 0 , D x i α i sin ( x i ) = x i 1 α i E 1 2 , 2 α i ( x i 2 ) , i = 1 , 2 .
Since the problem (10) is linear, we can split its solution as follows:
u ( x , t ) = v ( x , t ) + w ( x , t ) ,
where v ( x , t ) is the solution of the homogeneous problem:
D t γ v ( x , t ) i = 1 d D x i α i v ( x , t ) = 0 , ( x , t ) Ω × ( 0 , T ) , v ( x , t ) = 0 , ( x , t ) Ω × [ 0 , T ] , v ( x , 0 ) = φ ( x ) , x Ω ¯ ,
and w ( x , t ) is the solution of the problem:
D t γ w ( x , t ) i = 1 d D x i α i w ( x , t ) = p ( x ) f ( t ) , ( x , t ) Ω × ( 0 , T ) , w ( x , t ) = 0 , ( x , t ) Ω × [ 0 , T ] , w ( x , 0 ) = 0 , x Ω ¯ .
For every direction x i , 1 i d , the spectral problem associated with Problem (12) is:
D x i α i X i ( x i ) = λ X i ( x i ) , 0 x i 1 , with X i ( 0 ) = X i ( 1 ) = 0 .
The eigenfunctions of (14) are given by [9]:
X n i ( x ) = x i α i 1 E α i , α i ( λ n i x α i ) , n 1 ,
where λ n i are the associated eigenvalues. Notice that these eigenvalues form a countable set that is denoted by λ n i n 1 . Moreover, they satisfy the following:
(P1)
There are only finitely many real eigenvalues and the rest appears as complex conjugate pairs;
(P2)
| λ n i | ( 2 π n ) α i as n + ;
(P3)
α i π 2 < arg λ n i π for n sufficiently large and arg λ n i α i π 2 as n + .
Now, consider a real number μ such that:
μ γ π 2 , min α i π 2 , γ π ,
and this is possible since α i > γ and consequently the interval γ π 2 , min α i π 2 , γ π is not empty. Therefore, the property (P3) implies that
μ < α i π 2 < arg λ n i π for n large enough .
Since γ π 2 , min α i π 2 , γ π γ π 2 , min γ π , π , then applying Theorem 1, we obtain the estimate:
c 0 , t [ 0 , T ] , E γ , 1 ( λ n i t γ ) c 1 + | λ n i | t γ .
Remark 1.
In the published papers on this topic, the mistake commonly made is due to a misapplication of Theorem 1. Indeed, the following estimate
c 0 , n 1 , z [ 0 , 1 ] , E α i , β ( λ n i z α i ) c 1 + | λ n i | z α i
cannot hold true for all the eigenvalues λ n i , n 1 , since lim n arg λ n i = α i π 2 . Consequently, there is no real number μ α i π 2 , min ( π , α i π ) such that
n 1 , μ arg λ n i π .
Therefore, the conditions for applying Theorem 1 are not satisfied to obtain (18).
Based on [11], for every 1 i d , the family ( X n i ) n 1 is a basis of L 2 ( ( 0 , 1 ) , C ) , which is not orthogonal, because the operator D x i α i is not self-adjoint. Consequently, for each n 1 , we introduce the family:
X ^ n i ( x i ) = ( 1 x i ) α i 1 E α i , α i λ n i ( 1 x i ) α i ,
corresponding to the eigenfunctions of the adjoint operator of D x i α i , with respect to the inner product in L 2 ( ( 0 , 1 ) , C ) :
f , g L 2 ( ( 0 , 1 ) , C ) = 0 1 f ( x i ) g ¯ ( x i ) d x i ,
where g ¯ ( x i ) denotes the conjugate of the complex number g ( x i ) . In the sequel, the inner products in L 2 ( ( 0 , 1 ) , C ) and L 2 ( [ 0 , T ] , C ) will be simply denoted by · , · L 2 ( 0 , 1 ) and · , · L 2 ( 0 , T ) , respectively.
Notice that the adjoint operator of D x i α i is the right-sided Riemann–Liouville fractional derivative of order α i > 0 defined by:
x i D α i f ( x ) = 1 Γ ( n α i ) d n d t n x i 1 ( t x i ) n 1 α i f ( t ) d t , n = α i + 1 .
Moreover, the family ( X ^ n i ) n 1 is biorthogonal to ( X n i ) n 1 in L 2 ( ( 0 , 1 ) , C ) and satisfies
X n i , X ^ m i L 2 ( 0 , 1 ) = 0 for n m and X n i , X ^ n i L 2 ( 0 , 1 ) = C n i > 0 .
On the other hand, for any integers n 1 1 , , n d 1 , the resulting problem in the variable t is:
D t γ T ( t ) = i = 1 d λ n i i T ( t ) , 0 < t < T , with T ( 0 ) = φ n 1 n d ,
whose solution is given by:
T n 1 n d ( t ) = φ n 1 n d E γ , 1 i = 1 d λ n i i t γ ,
where:
φ n 1 n d = φ , i = 1 d X ^ n i i L 2 ( Ω ) i = 1 d X n i i , i = 1 d X ^ n i i L 2 ( Ω ) ,
since φ L 2 ( Ω ) . The solution of Problem (12) is then formally given by:
v ( x , t ) = n 1 1 n d 1 φ n 1 n d E γ , 1 i = 1 d λ n i i t γ i = 1 d X n i i .
In the same way, we can look for the classical solution of the problem (13) as:
w ( x , t ) = n 1 1 n d 1 w n 1 n d ( t ) i = 1 d X n i i ,
where the coefficients w n 1 n d ( t ) are to be found. To this end, if p L 2 ( Ω ) , then:
p ( x ) = n 1 1 n d 1 p n 1 n d i = 1 d X n i i ,
with
p n 1 n d = p , i = 1 d X ^ n i i L 2 ( Ω ) i = 1 d X n i i , i = 1 d X ^ n i i L 2 ( Ω ) ·
Substituting the solution w given by (25) in Equation (13), we find formally:
n 1 1 n d 1 D t γ i = 1 d λ n i i w n 1 n d ( t ) i = 1 d X n i i ( x i ) = f ( t ) n 1 1 n d 1 p n 1 n d i = 1 d X n i i ( x i ) .
Therefore, for each n 1 , , n d , the coefficient w n 1 n d ( t ) satisfies the equation:
D t γ w n 1 n d ( t ) i = 1 d λ n i i w n 1 n d ( t ) = f ( t ) p n 1 n d .
Equation (28) is supplemented with the initial condition:
w n 1 n d ( 0 ) = w ( · , 0 ) , i = 1 d X ^ n i i L 2 ( Ω ) i = 1 d X n i i , i = 1 d X ^ n i i L 2 ( Ω ) = 0 .
Following [12], the solution of (28) is:
w n 1 n d ( t ) = p n 1 n d 0 t s γ 1 E γ , γ i = 1 d λ n i i s γ f ( t s ) d s .
Therefore, the solution of problem (13) is written as:
w ( x , t ) = n 1 1 n d 1 p n 1 n d i = 1 d X n i i ( x i ) 0 t s γ 1 E γ , γ i = 1 d λ n i i s γ f ( t s ) d s .
In what follows, c is a constant that might change from line to line.
Theorem 2.
Let p C Ω ¯ , R , f C ( [ 0 , T ] , R ) and φ C Ω ¯ , R such that D x 1 α 1 D x d α d p C ( Ω ¯ , R ) , D x 1 α 1 D x d α d φ C ( Ω ¯ , R ) and φ | Ω = 0 . Then, there exists a unique classical solution of Problem (10).
Proof. 
Existence of a solution:
To lighten the proof without loss of generality, we present the case d = 1 and explain below the general case d 1 . In this situation, we consider the problem:
D t γ u ( x , t ) D x α u ( x , t ) = p ( x ) f ( t ) , ( x , t ) ( 0 , 1 ) × ( 0 , T ) , u ( x , t ) = 0 , ( x , t ) { 0 , 1 } × [ 0 , T ] , u ( x , 0 ) = φ ( x ) , x [ 0 , 1 ] ,
with 0 < γ < 1 and 1 < α < 2 . We need to prove the convergence of series (24), (31) and the series associated with their fractional derivatives in space and time. We set Φ : = D x α φ L 2 ( 0 , 1 ) . Then:
Φ = n 1 Φ n X n
where
Φ n = Φ , X ^ n L 2 ( 0 , 1 ) X n , X ^ n L 2 ( 0 , 1 ) = D x α φ , X ^ n L 2 ( 0 , 1 ) X n , X ^ n L 2 ( 0 , 1 ) = φ , x D α X ^ n L 2 ( 0 , 1 ) X n , X ^ n L 2 ( 0 , 1 ) = λ n ¯ φ , X ^ n L 2 ( 0 , 1 ) X n , X ^ n L 2 ( 0 , 1 ) = λ n ¯ φ n
Then, there is a constant c 0 such that, for every x [ 0 , 1 ] , t [ 0 , T ] and n are sufficiently large such that:
φ n X n ( x ) E γ , 1 ( λ n t γ ) = Φ n X n ( x ) λ n ¯ E γ , 1 ( λ n t γ ) c λ n c n α .
The series of general terms a n = c n α is convergent since α > 1 . Using the M-test of Weierstrass, the series (24) is normally convergent on [ 0 , 1 ] × [ 0 , T ] . Consequently, the function v ( x , t ) defined by (24) is continuous on [ 0 , 1 ] × [ 0 , T ] . The same idea can be used for the series (31) by setting P = D x α p L 2 ( 0 , 1 ) . Hence, the functions v and w verify the initial and the boundary conditions of (12) and (13), respectively.
Finally, to prove the convergence of the series associated to the fractional derivatives of v and w, it suffices to prove the convergence of series associated with: D t γ v , D t γ w , D x α v and D x α w . To this end, let
v n ( x , t ) = φ n X n ( x ) E γ , 1 ( λ n t γ ) ,
and
w ˜ n ( x , t ) = p n X n ( x ) 0 t s γ 1 E γ , γ ( λ n s γ ) f ( t s ) d t ;
so v ( x , t ) = n = 1 v n ( x , y , t ) and w ( x , t ) = n = 1 w ˜ n ( x , t ) .
From above, for any n 1 , x ( 0 , 1 ) and t ( 0 , T ) , one has:
D t γ v n ( x , t ) = φ n X n ( x ) λ n E γ , 1 ( λ n t γ ) = Φ n λ n ¯ X n ( x ) × λ n E γ , 1 ( λ n t γ ) c n α × | λ n | 1 + | λ n | t γ ( c is independent of x , t , and n ) c n α t γ ·
and
D x α v n ( x , t ) = λ n φ n X n ( x ) E γ , 1 ( λ n t γ ) = D t γ v n ( x , t ) c n α t γ ·
Therefore, the series n 1 D t γ v n ( x , t ) and n 1 D x α v n ( x , t ) converges normally on [ 0 , 1 ] × [ ε , T ] , for every ε ( 0 , T ) .
Applying the same ideas, we obtain the normal convergence of the series n 1 D t γ w ˜ n ( x , t ) and n 1 D x α w ˜ n ( x , t ) on [ 0 , 1 ] × [ ε , T ] , for every ε ( 0 , T ) . Whence, we obtain the desired convergence on the whole set [ 0 , 1 ] × ( 0 , T ] by Lemma 1, which achieves the claim.
Uniqueness of the solution:
Let u 1 , u 2 be two classical solutions of Problem (10) and set u * u 1 u 2 . Then, u * satisfies the following problem:
D t γ u * ( x , t ) D x α u * ( x , t ) = 0 , ( x , t ) ( 0 , 1 ) × ( 0 , T ) u * ( 0 , t ) = u * ( 1 , t ) = 0 , t [ 0 , T ] , u * ( x , 0 ) = 0 , x [ 0 , 1 ] .
Proceeding as above with φ 0 , we obtain that φ n = 0 , for every n 1 , which implies that u * 0 . □
Remark 2.
The generalization of the previous proof to the case d 1 proceeds as follows: We set Φ : = D x 1 α 1 D x d α d φ L 2 ( Ω ) . Then,
Φ ( x ) = n 1 1 n d 1 Φ n 1 n d i = 1 d X n i i ( x i ) .
Direct computations give:
Φ n 1 n d = i = 1 d λ n i i ¯ × φ n 1 n d
  • The analogous of the estimate (33) becomes, for any x Ω ¯ , t [ 0 , T ] and n 1 , , n d sufficiently large:
    φ n 1 n d i = 1 d X n i i ( x i ) E γ , 1 i = 1 d λ n i i t γ c i = 1 d λ n i i c i = 1 d n i α i ·
  • The analogous of the time derivative estimate (34) becomes, for any x Ω , t ( 0 , T ) and n 1 , , n d sufficiently large:
    D t γ v n 1 n d ( x , t ) = φ n 1 n d k = 1 d X n k k ( x k ) k = 1 d λ n k k E γ , 1 k = 1 d λ n k k t γ Φ n 1 n d k = 1 d λ n k k ¯ k = 1 d X n k k ( x k ) × k = 1 d λ n k k E γ , 1 k = 1 d λ n k k t γ c k = 1 d n k α k × k = 1 d λ n k k 1 + t γ k = 1 d λ n k k c t γ k = 1 d n k α k ·
  • The analogues of the space-derivative estimate (34) becomes, for any x Ω , t ( 0 , T ) and n 1 , , n d sufficiently large:
    k = 1 d D x k α k v n 1 n d ( x , t ) = k = 1 d λ n k k φ n 1 n d k = 1 d X n k k ( x k ) E γ , 1 k = 1 d λ n k k t γ = D t γ v n 1 n d ( x , t ) c t γ k = 1 d n k α k ·
    These new estimates lead to the same conclusions about the Fourier series in the case d 1 as in the case d = 1 .
Theorem 3.
Under the same conditions as in Theorem 2, the solution of the direct problem (10) depends continuously on the given data.
Proof. 
The arguments developed in Remark 2 allow, without loss of generality, to restrict ourselves to the case d = 1 .
Let u and u ^ be the solutions to the direct problem, corresponding to the data { p , f , φ } and { p ^ , f ^ , φ ^ } , respectively. Using (11), we find
| u ( x , t ) u ^ ( x , t ) | | v ( x , t ) v ^ ( x , t ) | + | w ( x , t ) w ^ ( x , t ) | ,
where { v , w } and { v ^ , w ^ } correspond to the data { p , f , φ } and { p ^ , f ^ , φ ^ } , respectively. The same arguments as applied to (24) before lead to
| v ( x , t ) v ^ ( x , t ) | c φ φ ^ L ( 0 , 1 ) ,
where c is a positive constant. On the other hand, a direct computation gives
w n w ^ n 0 t | ( t s ) γ E γ , γ ( λ n ( t s ) γ ) ( f ( s ) p n f ^ ( s ) p ^ n ) | d s .
It holds that
| f ( s ) p n f ^ ( s ) p ^ n | = | f ( s ) p n f ( s ) p ^ n + f ( s ) p ^ n f ^ ( s ) p ^ n | f L ( ( 0 , T ] ) | p n p ^ n | + | p ^ n | f f ^ L ( 0 , T ) .
Therefore, there is a constant c > 0 such that
| w ( x , t ) w ^ ( x , t ) | c ( f L ( 0 , 1 ) p p ^ L ( 0 , 1 ) + p ^ L ( 0 , 1 ) f f ^ L ( 0 , T ) ) ,
which achieves the proof. □

4. The Inverse Problem

Let us consider the problem:
D t γ u ( x , t ) i = 1 d D x i α i u ( x , t ) = p ( x ) f ( t ) , ( x , t ) Ω × ( 0 , T ) , u ( x , t ) = 0 , ( x , t ) Ω × [ 0 , T ] , u ( x , 0 ) = φ ( x ) , x Ω ¯ , Ω u ( x , t ) d x = ψ ( t ) , t [ 0 , T ] .
where the functions p, φ and ψ are given, and the functions u and f are unknown.
Example 2.
As above, consider the case d = 2 and Ω = ( 0 , π ) × ( 0 , π ) , with 0 < γ < 1 and 1 < α 1 , α 2 < 2 . We fix the initial condition φ ( x 1 , x 2 ) , the spatial source term p ( x 1 , x 2 ) for any ( x 1 , x 2 ) Ω and the “mean value” ψ ( t ) , t [ 0 , T ] of the solution u ( x 1 , x 2 , t ) on Ω:
φ ( x 1 , x 2 ) = 0 , p ( x 1 , x 2 ) = sin ( x 1 ) sin ( x 2 ) + x 1 1 α 1 E 1 2 , 2 α 1 ( x 1 2 ) sin ( y ) + sin ( x ) x 2 1 α 2 E 1 2 , 2 α 2 ( x 2 2 ) , ψ ( t ) = 4 E γ , 1 t γ , t [ 0 , T ] .
The unknown temporal source term is given by f ( t ) = E γ , 1 t γ and the exact solution u ( x 1 , x 2 , t ) = sin ( x 1 ) sin ( x 2 ) E γ , 1 t γ .
To solve Problem (35), we start by expanding u ( x , t ) and p ( x ) , using systems (15) and obtain:
u ( x , t ) = n 1 1 n d 1 c n 1 n d ( t ) i = 1 d X n i i ( x i ) ,
with
c n 1 n d ( t ) = u ( · , t ) , i = 1 d X ^ n i i L 2 ( Ω ) i = 1 d X n i i , i = 1 d X ^ n i i L 2 ( Ω ) ,
p ( x ) = n 1 1 n d 1 p n 1 n d i = 1 d X n i i ( x i )
and
p n 1 n d = p , i = 1 d X ^ n i i L 2 ( Ω ) i = 1 d X n i i , i = 1 d X ^ n i i L 2 ( Ω ) ·
By substitution in (35), we formally obtain:
n 1 1 n d 1 D t γ i = 1 d λ n i i c n 1 n d ( t ) i = 1 d X n i i ( x i ) = f ( t ) n 1 1 n d 1 p n 1 n d i = 1 d X n i i ( x i ) .
It follows that
D t γ c n 1 n d ( t ) i = 1 d λ n i i c n 1 n d ( t ) = f ( t ) p n 1 n d , 0 < t < T ,
with the initial condition c n 1 n d ( 0 ) = φ n 1 n d and φ n 1 n d is defined by (23). Using (36) and (6) we obtain:
ψ ( t ) = Ω u ( x , t ) d x = n 1 1 n d 1 c n 1 n d ( t ) i = 1 d 0 1 x i α i 1 E α i , α i k = 1 d λ n k k x α k d x i = n 1 1 n d 1 i = 1 d E α i , α i + 1 k = 1 d λ n k k × c n 1 n d ( t ) .
Let us set
Γ n 1 n d = i = 1 d E α i , α i + 1 k = 1 d λ n k k .
Using (41), we formally find:
D t γ ψ ( t ) n 1 1 n d 1 i = 1 d λ n i i Γ n 1 n d × c n 1 n d ( t ) = f ( t ) n 1 1 n d 1 Γ n 1 n d p n 1 n d ·
Therefore, we obtain:
f ( t ) = D t γ ψ ( t ) n 1 1 n d 1 i = 1 d λ n i i Γ n 1 n d × c n 1 n d ( t ) n 1 1 n d 1 Γ n 1 n d p n 1 n d ·
According to [12], the solution of (41) is given by:
c n 1 n d ( t ) = φ n 1 n d ( t ) E γ , 1 i = 1 d λ n i i t γ + p n 1 n d 0 t s γ 1 E γ , γ i = 1 d λ n i i s γ f ( t s ) d s .
Replacing c n 1 n d ( t ) given by (44) in (43), we then find the integral equation satisfied by the unknown function f ( t ) :
f ( t ) = 1 a D t γ ψ ( t ) ϕ ( t ) 0 t s γ 1 f ( t s ) ρ ( s ) d s = 1 a D t γ ψ ( t ) ϕ ( t ) 0 t ( t s ) γ 1 f ( s ) ρ ( t s ) d s ,
where
a = n 1 1 n d 1 Γ n 1 n d p n 1 n d
ϕ ( t ) = n 1 1 n d 1 i = 1 d λ n i i φ n 1 n d E γ , 1 i = 1 d λ n i i t γ Γ n 1 n d
and
ρ ( t ) = n 1 1 n d 1 i = 1 d λ n i i p n 1 n d E γ , γ i = 1 d λ n i i t γ Γ n 1 n d ·
Theorem 4.
Let
  • p C Ω ¯ , R such that D x 1 α 1 D x d α d p C Ω ¯ , R and Ω p ( x ) d x 0 .
  • Let ψ C ( [ 0 , T ] , R ) such that D t γ ψ C ( [ 0 , T ] , R ) .
  • φ C Ω ¯ , R such that D x 1 α 1 D x d α d φ C Ω ¯ , R and φ | Ω = 0 .
Then, there exists at least one solution to Problem (35).
Before proving Theorem 4, we introduce the functional framework for the fixed point integral Equation (45). First, applying the same arguments as before, we can show that the functions ρ and ϕ defined by (48) and (47) are γ -differentiable on ] 0 , T [ . Moreover, the condition Ω p ( x ) d x 0 implies that the real number a 0 .
At this stage, we will establish the existence of at least a solution to the inverse problem (35) and that this solution depends continuously on the data. First of all, we will use the Schauder fixed point theorem in Banach spaces with the Arzela–Ascoli compactness result. To this end, we start by defining the following operator:
B : C ( [ 0 , T ] , R ) C ( [ 0 , T ] , R ) f B ( f ) : t 1 a D t γ ψ ( t ) ϕ ( t ) 0 t ( t s ) γ 1 f ( s ) ρ ( t s ) d s .
To prove that the operator B admits a fixed point, start by showing that B maps a certain closed convex set into itself, in the space C ( [ 0 , T ] , R ) equipped with the Bielecki norm. For every δ > 0 , we introduce the Bielecki norm:
u δ = sup t [ 0 , T ] e δ t | u ( t ) | .
The space C ( [ 0 , T ] , R ) , · δ is a Banach space, and the two norms · δ and · are equivalent.
Lemma 2.
Under the above notations, there exists a positive constant δ * > 0 such that for any δ > δ * , there is a radius R > 0 such that the closed convex ball
K = { f C ( [ 0 , T ] , R ) : f δ R }
is stable by the operator B; that is, B ( K ) K .
Proof. 
For any f C ( [ 0 , T ] , R ) , t [ 0 , T ] and δ > 0 , we have:
e δ t | B ( f ) ( t ) | = e δ t | a | D t γ ψ ( t ) ϕ ( t ) 0 t ( t s ) γ 1 ρ ( t s ) f ( s ) d s , = 1 | a | e δ t D t γ ψ ( t ) ϕ ( t ) 0 t e δ ( t s ) ( t s ) γ 1 ρ ( t s ) e δ s f ( s ) d s , 1 | a | D t γ ψ δ + ϕ δ + f δ 0 T e δ s s γ 1 ρ ( s ) d s .
At this stage, recall that for any γ > 0 , the function E δ : s 1 | a | e δ s s γ 1 ρ ( s ) is integrable on [ 0 , T ] . Moreover, it holds that lim δ + E δ L 1 ( [ 0 , T ] ) = 0 ; consequently:
δ * > 0 , δ δ * , 1 | a | 0 T e δ s s γ 1 ρ ( s ) d s 1 | a | 0 T e δ * s s γ 1 ρ ( s ) d s < 1 .
That is, there are two real constants c 1 ( δ ) : = 1 | a | D t γ ψ δ + ϕ δ and c 2 : = 1 | a | 0 T e δ * s s γ 1 ρ ( s ) d s satisfying c 1 ( δ ) 0 and 0 c 2 < 1 such that:
δ δ * , B ( f ) δ c 1 ( δ ) + c 2 f δ .
Whence, for any δ δ * , there exists R > 0 such that the closed ball K of radius R in C ( [ 0 , T ] , R ) , · δ is stable by the operator B; this achieves the proof. □
Remark 3.
Using the Bielecki norm allows one to have no constraint on the maximum value that T can take. On the other hand, if we use the classical infinite norm, then T must be less than a finite quantity depending on the data of the problem.
Lemma 3.
The family B ( f ) f K is equicontinuous, that is:
t [ 0 , T ] , ε > 0 , r > 0 : t [ 0 , T ] ] t r , t + r [ , f K , | B ( f ) ( t ) B ( f ) ( t ) | ε .
Proof. 
It suffices to prove that the family B 1 ( f ) f K is equicontinuous, where
B 1 : C ( [ 0 , T ] , R ) C ( [ 0 , T ] , R ) f B 1 ( f ) : t 0 t ( t s ) γ 1 ρ ( t s ) f ( s ) d s .
We introduce the functions W 1 : s s γ 1 ρ ( s ) and W 2 : s e δ s W 1 ( s ) , which are clearly integrable on [ 0 , T ] . Let t and t in ] 0 , T ] . We can assume, without loss of generality, that 0 < t t T .
B 1 ( f ) ( t ) B 1 ( f ) ( t ) = 0 t W 1 ( t s ) W 1 ( t s ) f ( s ) d s + e δ t t t W 2 ( t s ) e δ s f ( s ) d s e δ T f δ 0 t | W 1 ( t s ) W 1 ( t s ) | d s + 0 t t | W 2 ( s ) | d s .
First, for every s ( 0 , t ) , we have
| W 1 ( t s ) W 1 ( t s ) | ( t s ) γ 1 | ρ ( t s ) | + ( t s ) γ 1 | ρ ( t s ) | , 2 ρ ( t s ) γ 1 .
Using the fact that the map s ( t s ) γ 1 is in L 1 ( 0 , t ) and applying the dominated convergence theorem of Lebesgue, we deduce that
lim t t 0 t | W 1 ( t s ) W 1 ( t s ) | d s = 0 .
Fixing now an ε > 0 and t ( 0 , T ) , it follows that
r 1 > 0 , t [ 0 , T ] ] t r 1 , t + r 1 [ , 0 t | W 1 ( t s ) W 1 ( t s ) | d s e δ T 4 R ρ ε .
On the other hand, using the fact that the map W 2 is in L 1 ( 0 , T ) and applying the dominated convergence theorem of Lebesgue, we deduce that
lim t t 0 t t | W 2 ( s ) | d s = 0 .
That is:
r 2 > 0 , t [ 0 , T ] ] t r 2 , t + r 2 [ , 0 t t | W 2 ( s ) | d s e δ T 2 R ε .
Let us point out that radii r 1 and r 2 depend only on t and ε , but are independent of f. Whence, for any t [ 0 , T ] and ε > 0 , there exists r : = min { r 1 , r 2 } > 0 such that
t [ 0 , T ] ] t r , t + r [ , f K , | B 1 ( f ) ( t ) B 1 ( f ) ( t ) | ε ,
which achieves the proof. □
By Arzela–Ascoli theorem, B ( K ) is relatively compact. Now, we are able to use the Schauder fixed-point theorem [13].
Theorem 5
(Schauder fixed-point theorem). Let ( E , . ) be a Banach space, and let K E be convex and closed. Let L : K K be a continuous operator such that L ( K ) is relatively compact. Then, L has a fixed point in K.
Proof of Theorem 4.
We know from Lemmas 2 and 3 that B is a continuous operator and B ( K ) is relatively compact, with K convex and closed. By Theorem 5, the operator B admits a fixed point in K, so the integral Equation (45) has a solution in [ 0 , T ] .
Once the existence of the source term f is established, the existence of u follows in a similar way with the arguments developed above for the direct problem. □
Theorem 6.
Under the conditions of Theorem 4, every solution of the inverse problem (35) depends continuously on the given data φ, p and ψ.
Proof. 
Let { u ( x , t ) , f ( t ) } , { u ˜ ( x , t ) , f ˜ ( t ) } be two solution sets of the inverse problem, corresponding to the data { p , φ , ψ } , { p ˜ , φ ˜ , ψ ˜ } , respectively. Using (45), we have:
f ( t ) = D t γ ψ ( t ) a ϕ ( t ) a 0 t ( t s ) γ 1 f ( s ) ρ ( t s ) a d s ,
and
f ˜ ( t ) = D t γ ψ ˜ ( t ) a ˜ ϕ ˜ ( t ) a ˜ 0 t ( t s ) γ 1 f ˜ ( s ) ρ ˜ ( t s ) a ˜ d s ,
where { a , ϕ , ρ } , { a ˜ , ϕ ˜ , ρ ˜ } correspond to the data { p , φ , ψ } , { p ˜ , φ ˜ , ψ ˜ } , respectively.
Remark first that ρ , a and ϕ depend continuously on p and φ . Indeed, for the function ρ , using (39) and (48), we have the estimate:
ρ ρ ˜ L ( 0 , T ) c p p ˜ L ( Ω ) .
For the real constant a, using (46) and (39), we have the estimate:
| a a ˜ | c p p ˜ L ( Ω ) .
For the function ϕ , using (47) and (23), we have the estimate:
| ϕ ( t ) ϕ ˜ ( t ) | c φ φ ˜ L ( Ω ) .
Let g ( t ) and g ˜ ( t ) be such that
g ( t ) = 0 t ( t s ) γ 1 f ( s ) ρ ( t s ) a d s , g ˜ ( t ) = 0 t ( t s ) γ 1 f ˜ ( s ) ρ ˜ ( t s ) a ˜ d s .
Then, for any δ > 0 :
e δ t | f ( t ) f ˜ ( t ) | e δ t D t γ ψ ( t ) a D t γ ψ ˜ ( t ) a ˜ + ϕ ( t ) a ϕ ˜ ( t ) a ˜ + | g ( t ) g ˜ ( t ) | .
Now,
e δ t D t γ ψ ( t ) a D t γ ψ ˜ ( t ) a ˜ D t γ ψ ( t ) a D t γ ψ ( t ) a ˜ + D t γ ψ ( t ) a ˜ D t γ ψ ˜ ( t ) a ˜ D t γ ψ L ( 0 , T ) 1 a 1 a ˜ + 1 | a ˜ | D t γ ψ D t γ ψ ˜ L ( 0 , T ) ,
and consequently,
e δ t D t γ ψ ( t ) a D t γ ψ ˜ ( t ) a ˜ c M 1 | a a ˜ | p p ˜ L ( Ω ) + 1 | a ˜ | D t γ ψ D t γ ψ ˜ δ ,
where M 1 is a positive constant such that D t γ ψ L ( 0 , T ) M 1 . Similarly, we obtain:
e δ t ϕ ( t ) a ϕ ˜ ( t ) a ˜ c ϕ L ( 0 , T ) | a a ˜ | p p ˜ L ( Ω ) + 1 | a ˜ | ϕ ϕ ˜ L ( Ω ) .
It follows that
e δ t ϕ ( t ) a ϕ ˜ ( t ) a ˜ c M 2 | a a ˜ | p p ˜ L ( Ω ) + c | a ˜ | φ φ ˜ L ( Ω ) ,
where M 2 is a positive constant such that ϕ L ( Ω ) M 2 . Finally, we obtain:
e δ t | g ( t ) g ˜ ( t ) | = e δ t 0 t ( t s ) γ 1 f ( t ) ρ ( t s ) a f ˜ ( t ) ρ ˜ ( t s ) a ˜ d s e δ t 0 t ( t s ) γ 1 f ( t ) ρ ( t s ) a f ˜ ( t ) ρ ( t s ) a d s + e δ t 0 t ( t s ) γ 1 f ˜ ( t ) ρ ( t s ) a f ˜ ( t ) ρ ˜ ( t s ) a ˜ d s f ˜ L ( 0 , T ) 0 t ( t s ) γ 1 ρ ( t s ) a ρ ˜ ( t s ) a ˜ d s + M 3 | a | f f ˜ δ I ( δ ) ,
where M 3 is a positive constant such that ρ L ( 0 , T ) M 3 and
I ( δ ) = 0 T ( t s ) γ 1 e δ ( t s ) d s .
Similarly, we obtain:
ρ ( t s ) a ρ ˜ ( t s ) a ˜ c M 3 | a a ˜ | p p ˜ L ( Ω ) + 1 | a ˜ | ρ ρ ˜ L ( 0 , T ) , c ( M 3 + 1 ) | a a ˜ | p p ˜ L ( Ω ) .
Then,
e δ t 0 t ( t s ) γ 1 ρ ( t s ) a ρ ˜ ( t s ) a ˜ d s c ( M 3 + 1 ) | a a ˜ | p p ˜ L ( Ω ) T γ γ .
Therefore, we obtain the estimate:
| g ( t ) g ˜ ( t ) | T γ γ c ( M 3 + 1 ) | a a ˜ | p p ˜ L ( Ω ) + M 3 I ( δ ) | a | f f ˜ L ( 0 , T ) .
Combining the last estimate and (55), we deduce the estimate on the source term:
1 M 3 I ( δ ) | a | f f ˜ δ c ˜ D t γ ψ D t γ ψ ˜ L ( 0 , T ) + φ φ ˜ L ( Ω ) + p p ˜ L ( Ω ) ,
where c ˜ is a positive constant depending on a, a ˜ , T, M 1 , M 2 and M 3 .
Since lim δ + I ( δ ) = 0 , there is δ * > 0 such that
1 M 3 I ( δ ) | a | > 0 , for any δ δ *
holds true. Whence, there is a positive constant C > 0 , independent of f and f ˜ such that
f f ˜ δ C D t γ ψ D t γ ψ ˜ L ( 0 , T ) + φ φ ˜ L ( Ω ) + p p ˜ L ( Ω ) .
That is, the solution f depends continuously on the data.
The continuity of the part u of the solution follows with classical arguments for linear problems with the norm: u , δ = sup x Ω ¯ sup t [ 0 , T ] e δ t u ( x , t ) . □

5. Conclusions

We studied the inverse problem of finding the solution of a generalized time-space fractional equation and the source term, knowing the spatial mean of the solution at any time. The existence and the continuity with respect to the data of the solution for the direct and the inverse problem are proved by Fourier’s method and the Schauder fixed-point theorem, in an adequate convex bounded subset.
In the literature concerning this topic, an incorrect use of the estimates in the generalized Mittag–Leffler functions is commonly performed. This leads to false proofs of the Fourier series’ convergence to recover the partial differential equation satisfied by the solution, the initial data or the boundary conditions. In the present manuscript, the correct framework to recover the decay of fractional Fourier coefficients is established; this allows one to correctly recover the initial data, the boundary conditions and the partial differential equation satisfied within the space-time domain.

Author Contributions

Formal analysis, A.E.H. and M.K.; Investigation, A.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are very grateful to the anonymous referees for their valuable comments and suggestions and to Rafik IMEKRAZ for his careful proofreading and interesting remarks.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Aleroev, T.S.; Kirane, M.; Malik, S. Determination of a source term for a time fractional diffusion equation with an integral type over-determining condition. Electron. J. Differ. Equ. (EJDE) 2013, 2013, 1–16. [Google Scholar]
  2. Kirane, M.; Malik, S.A. Determination of an unknown source term and the temperature distribution for the linear heat equation involving fractional derivative in time. Appl. Math. Comput. 2011, 218, 163–170. [Google Scholar] [CrossRef] [Green Version]
  3. Furati, K.M.; Iyiola, O.S.; Kirane, M. An inverse problem for a generalized fractional diffusion. Appl. Math. Comput. 2014, 249, 24–31. [Google Scholar] [CrossRef]
  4. Tatar, S.; Ulusoy, S. An inverse source problem for a one-dimensional space–time fractional diffusion equation. Appl. Anal. 2014, 94, 2233–2244. [Google Scholar] [CrossRef]
  5. Tatar, S.; Tinaztepe, R.; Ulusoy, S. Simultaneous inversion for the exponents of the fractional time and space derivatives in the space-time fractional diffusion equation. Appl. Anal. 2014, 95, 1–23. [Google Scholar] [CrossRef]
  6. Kirane, M.; Malik, S.; Al-Gwaiz, M. An inverse source problem for a two dimensional time fractional diffusion equation with nonlocal boundary conditions. Math. Methods Appl. Sci. 2013, 36, 1056–1069. [Google Scholar] [CrossRef]
  7. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives; Gordon and Breach Science Publishers: Yverdon-les-Bains, Switzerland, 1993; Volume 1993. [Google Scholar]
  8. Diethelm, K. The Analysis of Fractional Differential Equations: An Application-Oriented Exposition Using Differential Operators of Caputo Type; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  9. Rogosin, S.; Gorenflo, R.; Kilbas, A.; Mainardi, F. Mittag-Leffler Function, Related Topics and Applications; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar] [CrossRef]
  10. Podlubny, I. Fractional Differential Equations; Number 198 in Mathematics in science and engineering; Academic Press: San Diego, CA, USA, 1999; OCLC: 245654622. [Google Scholar]
  11. Aleroev, T.S.; Khasambiev, M.V. Kraevaya zadacha dlya odnomernogo drobnogo differentsial’nogo uravneniya advektsii-diffuzii (in Russian) [Boundary Value Problem for one-dimensional differential advection-dispersion equation]. Vestnik MGSU [Proc. Mosc. State Univ. Civ. Eng.] 2014, 6, 71–76. [Google Scholar] [CrossRef]
  12. Luchko, Y.; Gorenflo, R. An operational method for solving fractional differential equations with the Caputo derivatives. Acta Math. Vietnam 1999, 24, 207–233. [Google Scholar]
  13. Zeidler, E. Fixed-Point Theorems, 2nd ed.; Number Eberhard Zeidler; 1 in Nonlinear Functional Analysis and Its Applications; Springer: New York, NY, USA, 1993; OCLC: 180672677. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

El Hamidi, A.; Kirane, M.; Tfayli, A. An Inverse Problem for a Non-Homogeneous Time-Space Fractional Equation. Mathematics 2022, 10, 2586. https://doi.org/10.3390/math10152586

AMA Style

El Hamidi A, Kirane M, Tfayli A. An Inverse Problem for a Non-Homogeneous Time-Space Fractional Equation. Mathematics. 2022; 10(15):2586. https://doi.org/10.3390/math10152586

Chicago/Turabian Style

El Hamidi, Abdallah, Mokhtar Kirane, and Ali Tfayli. 2022. "An Inverse Problem for a Non-Homogeneous Time-Space Fractional Equation" Mathematics 10, no. 15: 2586. https://doi.org/10.3390/math10152586

APA Style

El Hamidi, A., Kirane, M., & Tfayli, A. (2022). An Inverse Problem for a Non-Homogeneous Time-Space Fractional Equation. Mathematics, 10(15), 2586. https://doi.org/10.3390/math10152586

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