Next Article in Journal
Design of the Novel Fractional Order Hybrid Whale Optimizer for Thermal Wind Power Generation Systems with Integration of Chaos Infused Wind Power
Previous Article in Journal
Positive Solutions and Their Existence of a Nonlinear Hadamard Fractional-Order Differential Equation with a Singular Source Item Using Spectral Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Inverse Problem for the Subdiffusion Equation with a Non-Local in Time Condition

by
Ravshan Ashurov
1,2,†,
Marjona Shakarova
2,† and
Sabir Umarov
3,*,†
1
Department of Mathematics, New Uzbekistan University, Movarounnahr Street 1, Tashkent 100000, Uzbekistan
2
Institute of Mathematics, Uzbekistan Academy of Science, Student Town Street 9, Tashkent 100174, Uzbekistan
3
Department of Mathematics, University of New Haven, 300 Boston Post Road, West Haven, CT 06516, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Fractal Fract. 2024, 8(7), 378; https://doi.org/10.3390/fractalfract8070378
Submission received: 10 May 2024 / Revised: 17 June 2024 / Accepted: 25 June 2024 / Published: 27 June 2024

Abstract

:
An inverse problem of determining the right-hand side of the abstract subdiffusion equation with a fractional Caputo derivative is considered in a Hilbert space H. For the forward problem, instead of the Cauchy condition, the non-local in time condition u ( 0 ) = u ( T ) is taken. The right-hand side of the equation has the form g ( t ) f with a given function g ( t ) and an unknown element f H . If the function g ( t ) preserves its sign, then under a over-determined condition u ( t 0 ) = ψ , t 0 ( 0 , T ) , it is proved that the solution of the inverse problem exists and is unique. An example is given showing the violation of the uniqueness of the solution for some sign-changing functions g ( t ) . For such functions g ( t ) , under certain conditions on this function, one can achieve the well-posedness of the problem by choosing properly t 0 . Moreover, we show that for some g ( t ) , for the existence of a solution to the inverse problem, certain orthogonality conditions must be satisfied, but in this case there is no uniqueness. To the best knowledge of authors, the inverse problem with the non-local condition u ( 0 ) = u ( T ) has been considered for the first time. Moreover, all the results obtained are new not only for the subdiffusion equation, but also for the classical diffusion equation.

1. Introduction

In recent decades, it has been discovered that fractional differential equations are an effective mathematical tool for modeling various anomalous phenomena (see, e.g., [1,2,3,4,5]). In this regard, researchers are actively finding new methods for solving fractional-order differential equations and discovering new properties of solutions (see, e.g., [1,6,7,8,9]).
This work is devoted to the study of the inverse problem of determining the source function, the right-hand side of the subdiffusion equation. The appearance of inverse problems is due to the fact that when mathematically modeling various processes, the properties of the medium, which in direct problems are considered given quantities, often turn out to be unknown in practice. Mathematically, this means that, along with solving the direct problem, it is also necessary to determine either the coefficients of the equations, or the initial or boundary conditions, or the properties of the domain in which the process under study occurs, or the properties of the object being studied.
Let H be a separable Hilbert space with the scalar product ( · , · ) , the norm | | · | | , and A be an self-adjoint positive operator acting in H with a domain of definition D ( A ) . Assume that the inverse of A is a compact operator. Then A has a complete system of orthonormal eigenfunctions { v k } and a countable set of positive eigenvalues λ k : 0 < λ 1 λ 2 + .
Let C ( ( a , b ) ; H ) stand for a set of continuous vector-valued functions h ( t ) , defined on t ( a , b ) with values in H. For h : R + H fractional analogs of integrals and derivatives are defined similar to scalar functions (see, e.g., [10]) and understood in the sense of Bochner. Recall that the Caputo fractional derivative of order ρ > 0 of function h ( t ) has the form (see, e.g., [9], p. 14)
D t ρ h ( t ) = 1 Γ ( 1 ρ ) 0 t h ( ξ ) ( t ξ ) ρ d ξ , t > 0 ,
provided the right-hand side exists, where Γ ( · ) is Euler’s gamma function. If ρ = 1 , then we assume that the fractional derivative coincides with the ordinary classical derivative: D t h ( t ) = d d t h ( t ) .
Let ρ ( 0 , 1 ] be a fixed number. Consider the following non-local boundary value problem
D t ρ u ( t ) + A u ( t ) = g ( t ) f , t ( 0 , T ] , u ( T ) = u ( 0 ) ,
where g ( t ) C [ 0 , T ] is a given scalar function, and f H is a known element of H. This case problem (1) of finding an unknown u ( t ) is called a forward problem. The forward problem is studied in [11].
The local boundary condition is assumed to be given at one specific point. The famous example is the Cauchy condition. In problem (1) the relation between the value of the sought function at the initial and final moments of time is specified. Therefore, such problems are usually called non-local problems (see, e.g., [12]). The physical meaning of the problem under consideration can, for example, be interpreted as the following control problem: find that part of the heat source g ( t ) f ( x ) that depends on the physical location of the source. That is f ( x ) , which would provide the same temperature at the final moment of time as the temperature at the initial moment.
The main purpose of our study is the inverse problem. Namely, to determine the right-hand side of Equation (1), namely element f H , when it is unknown. To solve the inverse problem, one needs an extra condition. Following A.I. Prilepko and A.B. Kostin [13] (see also K.B. Sabitov [14]) we consider the additional condition in the form:
u ( t 0 ) = ψ ,
where t 0 is a given fixed point of the interval ( 0 , T ) .
Usually, authors specify an additional condition (2) at the final time t 0 = T (see, e.g., [15,16,17] for classical diffusion equations and [18,19,20,21] for subdiffusion equations). This choice is convenient for real models, since when the process is over it is easy to measure u ( T ) . However, in some cases, choosing t 0 = T may violate the uniqueness of the solution of the inverse problem, and if t 0 < T is chosen, then it is possible to achieve uniqueness in these cases too (see, e.g., [22]).
We call the non-local boundary value problem (1) with unknown f on the right and equipped with the additional condition (2) an inverse problem. The solution to the inverse problem is understood in the following sense:
Definition 1.
A pair of functions { u ( t ) , f } such that D t ρ u ( t ) , A u ( t ) C ( ( 0 , T ] ; H ) , u ( t ) C ( [ 0 , T ] ; H ) , f H , and satisfying conditions (1)–(2), is called a solution of the inverse problem (1), (2).
Inverse problems of determining the right-hand side of various subdiffusion equations have been studied by a number of authors due to the importance of such problems for applications. Note that the methods for studying inverse problems depend on whether f or g ( t ) is unknown. In the case when the function g ( t ) is unknown, but f is known, the inverse problem with various additional conditions is equivalently reduced to solving an integral equation (see, e.g., [23,24,25,26]). For the subdiffusion and diffusion equations, the case g ( t ) 1 and the unknown f has been studied by many authors; in this case, the Fourier method is used to solve the inverse problem (see, e.g., [27,28,29,30,31,32] and references therein). In [33], the authors considered the inverse problem of the simultaneous determination of the order of the Riemann–Liouville fractional derivative and the source function in the subdiffusion equation. Using the Fourier method, the authors proved that the solution to this inverse problem exists and is unique. Note that in works [27,28,33], the additional condition (2) is specified at t 0 = T .
Let us focus on the case where g ( t ) 1 and f is unknown, since the present paper is devoted to this case. We note that in all the papers cited below, the Cauchy problem with the initial condition u ( 0 ) = ϕ is considered instead of a non-local condition u ( 0 ) = u ( T ) in problem (1). For the classical diffusion equation (i.e., ρ = 1 ), such an inverse problem with some additional conditions has been studied in sufficient detail in [34] (Chapter 8) and in papers [13,14,15,16,35]. Since the equation under consideration also includes the diffusion equation, we will dwell on these works in more detail. In works [15,16], abstract diffusion equations in Banach and Hilbert spaces with additional condition (2) at the final time t 0 = T , are considered. In [15], in the case of a Hilbert space, the elliptic part of the equation is self-adjoint, while in [16], non-self-adjoint elliptic operator is also included. The uniqueness criterion is found and the existence of a generalized solution is proved. The elliptic part of the diffusion equation in the work [13] is a second-order differential expression. Both non-self-adjoint and self-adjoint elliptic parts are considered. In this paper, g ( t ) also depends on a space variable: g ( t ) : = g ( x , t ) , and as an additional condition the authors take (2) and in some cases the integral condition. In the case of a self-adjoint elliptic part, the authors managed to find a criterion for the uniqueness of the generalized solution of the inverse problem. In works [14,35], the elliptic part of the equation is u x x , defined on the interval and the Laplace operator on the rectangle, respectively. Considering the overdetermined condition in the form (2), a criterion for the uniqueness of the classical solution is found and the Fourier method is used to construct a classic solution.
A similar inverse problem ( g ( t ) 1 and f is unknown) for the abstract subdiffusion equation (i.e., 0 < ρ < 1 ) with the Caputo derivative was studied in [36]. To define the function f, the authors used the following additional condition 0 T u ( t ) d μ ( t ) = u T . Authors of [18] (see also [19,20]) studied the uniqueness of the solution of the inverse problem for the subdiffusion equation, the elliptic part of which depends on time, with the additional condition u ( T ) = φ . It is proved that if the function g ( t ) does not change sign, then the solution of the inverse problem is unique. It should be especially noted that in [19] the authors constructed an example of a function g ( t ) that changes sign in the domain under consideration, as a result of which the uniqueness of the solution of the inverse problem is lost.
Finally, let us note paper [22], where in their Equation (1) the Laplace operator is taken as the operator A and, as noted above, the Cauchy condition is taken instead of the non-local condition. The authors have taken the additional condition (2) to find the unknowns. In terms of the function g ( t ) , a uniqueness criterion is found, and under the condition that the criterion is satisfied, the existence of a unique solution to the inverse problem is proved.
It is well known that the inverse problem is ill-posed, i.e., the solution does not depend continuously on the given data. Therefore, in papers [37,38] (see also references therein) various regularization methods are proposed for constructing an approximate solution of the inverse problem under the condition u ( T ) = φ .
In the current paper we consider a new inverse problem (1), (2). To the authors’ best knowledge, the above works contain no results concerning the correct solvability of such inverse problems as we study in the present paper.
The present paper is organized as follows. In Section 2, we recall some preliminary facts related to Mittag–Leffler functions used in this paper. Section 3 is devoted to the study of the inverse problem (1), (2) of determining the right-hand side of the equation. This section shows that if the function g ( t ) preserves its sign, then the solution to the inverse problem (1), (2) exists and is unique. An example is given showing the violation of the uniqueness of the solution in the case of some not sign-preserving functions g ( t ) . It is also proven that under certain conditions on the function g ( t ) it is possible to achieve well-posedness of the problem by choosing t 0 properly. The paper ends with the conclusion section.

2. Preliminaries

In this section, we recall some facts about Mittag–Leffler functions and the solution of the forward problem (1), which we will use in the following section.
The action of the abstract operator A under consideration on the element h H can be written as
A h = k = 1 λ k h k v k ,
where h k is the Fourier coefficient of the element h: h k = ( h , v k ) . Obviously, the domain of definition of this operator has the form
D ( A ) = { h H : k = 1 λ k 2 | h k | 2 < } .
For elements of D ( A ) , we introduce the norm
| | h | | 1 2 = k = 1 λ k 2 | h k | 2 = | | A h | | 2 ,
and together with this norm D ( A ) turns into a Hilbert space.
For 0 < ρ < 1 and an arbitrary complex number μ , let E ρ , μ ( z ) denote the two-parameter Mittag–Leffler function of the complex argument z (see, e.g., [5], p. 15):
E ρ , μ ( z ) = k = 0 z k Γ ( ρ k + μ ) .
If the parameter μ = 1 , then we have the classical Mittag–Leffler function: E ρ ( z ) = E ρ , 1 ( z ) .
Recall some properties of the Mittag–Leffler functions.
Lemma 1
(see [39], p. 136). For any t 0 one has
| E ρ , μ ( t ) | C 1 + t ,
where constant C does not depend on t and μ.
Lemma 2
(see [5], p. 47). The classical Mittag–Leffler function of the negative argument E ρ ( t ) is monotonically decreasing function for all 0 < ρ < 1 and
0 < E ρ ( t ) < 1 , E ρ ( 0 ) = 1 .
Lemma 3
(see [5], formula (4.4.17), p. 64). Let μ be an arbitrary complex number. Then, the following asymptotic formula
E ρ , μ ( t ) = k = 1 2 ( t ) k Γ ( μ k ρ ) + O ( t 3 ) , t ,
is valid.
Lemma 4
(see [5], formula (4.2.3), p. 57). For all α > 0 , μ C the following recurrence relation holds:
E ρ , μ ( t ) = 1 Γ ( μ ) t E ρ , μ + ρ ( t ) .
Lemma 5
(see [39], formula (1.16), p. 120). Let ρ > 0 , μ > 0 and λ C . Then, for all positive t one has
1 Γ ( μ ) 0 t ( t η ) μ 1 η ρ 1 E ρ , ρ ( λ η ρ ) d η = t μ + ρ 1 E ρ , ρ + μ ( λ t ρ ) .
To solve the inverse problem we use the following result on the forward problem: If f H and g C [ 0 , T ] then the unique solution of the forward problem has the form (see [11], Theorem 3 and Corollary 1)
u ( t ) = k = 1 ω k ( t ) + ω k ( T ) 1 E ρ ( λ k T ρ ) E ρ ( λ k t ρ ) v k ,
where
ω k ( t ) = f k 0 t ( t η ) ρ 1 E ρ , ρ ( λ k ( t η ) ρ ) g ( η ) d η ,
and f k , k 1 , are Fourier coefficients of the element f.

3. Existence and Uniqueness of the Inverse Problem

We apply the additional condition (2) to Equation (4) and denote by ψ k the Fourier coefficients of the element ψ : ψ k = ( ψ , v k ) . Then
k = 1 f k [ ( 1 E ρ ( λ k T ρ ) ) b k , ρ ( t 0 ) + E ρ ( λ k t 0 ρ ) b k , ρ ( T ) ] v k = k = 1 ψ k ( 1 E ρ ( λ k T ρ ) ) v k ,
where
b k , ρ ( t ) = 0 t ( t s ) ρ 1 E ρ , ρ ( λ k ( t s ) ρ ) g ( s ) d s .
Hence, to find f k , we obtain the following equation
f k Δ ρ ( k , t 0 , T ) = ψ k ( 1 E ρ ( λ k T ρ ) ) ,
where
Δ ρ ( k , t 0 , T ) = ( 1 E ρ ( λ k T ρ ) ) b k , ρ ( t 0 ) + E ρ ( λ k t 0 ρ ) b k , ρ ( T ) .
Of course the case Δ ρ ( k , t 0 , T ) = 0 is critical. Since, according to Lemma 2
0 < 1 E ρ ( λ k T ρ ) < 1 and 0 < E ρ ( λ k t 0 ρ ) < 1 ,
then this can happen only when g ( t ) changes its sign. Note that in this case, the functions b k , ρ ( t 0 ) and b k , ρ ( T ) can also change their signs.
Let us divide the set of natural numbers N into two groups K 0 , ρ and K ρ : N = K ρ K 0 , ρ . Namely, the number k belongs to K 0 , ρ , if Δ ρ ( k , t 0 , T ) = 0 , and k K ρ if Δ ρ ( k , t 0 , T ) 0 . Note that for some t 0 and T the set K 0 , ρ may be empty, then K ρ = N . For example, if g ( t ) is sign-preserving, then K ρ = N for all k regardless of the values t 0 and T.
Now we establish lower bounds for Δ ρ ( k , t 0 , T ) . The lemma below states that if g ( t ) does not change its sign, then K 0 , ρ is empty.
Lemma 6.
Let ρ ( 0 , 1 ] , g ( t ) C [ 0 , T ] and g ( t ) 0 , t [ 0 , T ] . Then, there is a constant C > 0 , depending on t 0 and T, such that for all k:
| Δ ρ ( k , t 0 , T ) | C λ k .
Proof. 
Since g ( t ) 0 on the closed interval [ 0 , T ] , we have either g ( t ) g + = c o n s t > 0 or g ( t ) g = c o n s t < 0 for all t [ 0 , T ] . It is sufficient to consider the first case. The second case can be treated in exactly the same way. For τ ( 0 , T ] , we have (see Lemma 5)
b k , ρ ( τ ) g + 0 τ η ρ 1 E ρ , ρ ( λ k η ρ ) d η = g + τ ρ E ρ , ρ + 1 ( λ k τ ρ ) .
Taking into account (see Lemma 4)
E ρ , ρ + 1 ( t ) = t 1 ( 1 E ρ ( t ) )
we obtain
b k , ρ ( τ ) 1 λ k ( 1 E ρ ( λ k τ ρ ) ) g + 1 λ k ( 1 E ρ ( λ 1 τ ρ ) ) g + C τ λ k , C τ > 0 .
Therefore
Δ ρ ( k , t 0 , T ) ( 1 E ρ ( λ k T ρ ) ) C t 0 λ k + E ρ ( λ k t 0 ρ ) C T λ k ( 1 E ρ ( λ k T ρ ) ) C t 0 λ k ,
which implies the desired assertion due to Lemma 2. □
Theorem 1.
Let ρ ( 0 , 1 ] , g ( t ) C [ 0 , T ] and g ( t ) 0 , t [ 0 , T ] . Assume ψ D ( A ) . Then there exists a unique solution of the inverse problem (1) and (2):
u ( t ) = k = 1 f k b k , ρ ( t ) + b k , ρ ( T ) 1 E ρ ( λ k T ρ ) E ρ ( λ k t ρ ) v k , f = k = 1 f k v k ,
where
f k = ψ k ( 1 E ρ ( λ k T ρ ) ) ( 1 E ρ ( λ k T ρ ) ) b k , ρ ( t 0 ) + E ρ ( λ k t 0 ρ ) b k , ρ ( T ) , k 1 .
Moreover, there is a positive constant C such that the following estimate holds:
| | D t ρ u ( t ) | | + | | A u ( t ) | | + | | f | | C | | ψ | | 1 , t > 0 .
Proof. 
The assumption of the theorem, due to Lemma 6, implies Δ ρ ( k , t 0 , T ) 0 for all k N . Hence, it follows from (5) that
f k = ψ k ( 1 E ρ ( λ k T ρ ) ) Δ ρ ( k , t 0 , T ) ,
and
u k ( t ) = f k b k , ρ ( t ) + b k , ρ ( T ) 1 E ρ ( λ k T ρ ) E ρ ( λ k t ρ ) .
With these Fourier coefficients, for f and u ( t ) we have the following infinite series representations:
f = k = 1 ψ k ( 1 E ρ ( λ k T ρ ) ) Δ ρ ( k , t 0 , T ) v k ,
and
u ( t ) = k = 1 f k b k , ρ ( t ) + b k , ρ ( T ) 1 E ρ ( λ k T ρ ) E ρ ( λ k t ρ ) v k .
Let F j be the partial sum of series (9):
F j = k = 1 j ψ k ( 1 E ρ ( λ k T ρ ) ) Δ ρ ( k , t 0 , T ) v k .
We show that the series F j is convergent in H .
For this, applying Parseval’s identity, we obtain
| | F j | | 2 = k = 1 j | ψ k ( 1 E ρ ( λ k T ρ ) ) Δ ρ ( k , t 0 , T ) | 2 .
The latter, by virtue of Lemmas 2 and 6, implies the following estimate:
| | F j | | 2 C k = 1 j λ k 2 | ψ k | 2 C | | ψ | | 1 2 .
Now, letting j and taking into account that ψ D ( A ) , we obtain f H .
Now, consider series (10) for u ( t ) . Using the representation (7) for f k , this series can be written in the form
u ( t ) = k = 1 ψ k ( 1 E ρ ( λ k T ρ ) ) Δ ρ ( k , t 0 , T ) b k , ρ ( t ) + b k , ρ ( T ) 1 E ρ ( λ k T ρ ) E ρ ( λ k t ρ ) v k .
Let S j ( t ) be the j-th partial sum of series (12). Then, by virtue of Parseval’s identity, we have
| | A S j ( t ) | | 2 = k = 1 j | λ k ψ k ( 1 E ρ ( λ k T ρ ) ) Δ ρ ( k , t 0 , T ) b k , ρ ( t ) + b k , ρ ( T ) 1 E ρ ( λ k T ρ ) E ρ ( λ k t ρ ) | 2 .
Applying Lemmas 2 and 6, we obtain
| | A S j ( t ) | | 2 C k = 1 j λ k 4 | ψ k | 2 | b k , ρ ( t ) | + | b k , ρ ( T ) | 2 ,
with a constant C > 0 . Due to Lemmas 1 and 5, we obtain
| | A S j ( t ) | | 2 C k = 1 j λ k 4 | ψ k | 2 max 0 t T | g ( t ) | 2 [ t ρ E ρ , ρ + 1 ( λ k t ρ ) + T ρ E ρ , ρ + 1 ( λ k T ρ ) ] 2 C k = 1 j λ k 4 | ψ k | 2 max 0 t T | g ( t ) | 2 t ρ 1 + λ k t ρ + T ρ 1 + λ k T ρ 2 C k = 1 j λ k 2 | ψ k | 2 max 0 t T | g ( t ) | 2           C | | ψ | | 1 2 max 0 t T | g ( t ) | 2 C 1 ψ 1 2 , t 0 .
Now, letting j , we obtain A u ( t ) C ( [ 0 , T ] ; H ) and in particular u ( t ) C ( [ 0 , T ] ; H ) . Further, it follows from Equation (1) that
D t ρ S j ( t ) = A S j ( t ) + k = 1 j f k g ( t ) v k , t > 0 .
Therefore, it follows from estimates (11) and (13), that
| | D t ρ S j ( t ) | | 2 C ψ 1 2 , t > 0 ,
and
D t ρ u ( t ) C ( ( 0 , T ] ; H ) .
To prove the uniqueness of the solution, assume that there exist two distinct pairs { u 1 , f 1 } and { u 2 , f 2 } , satisfying the inverse problem (1) and (2). We need to show that u ( t ) u 1 ( t ) u 2 ( t ) 0 , and f f 1 f 2 = 0 . For the pair { u ( t ) , f } , we have the infinitely many non-local boundary value problems:
D t ρ u ( t ) + A u ( t ) = f g ( t ) , t ( 0 , T ] , u ( T ) = u ( 0 ) , u ( t 0 ) = 0 , t 0 ( 0 , T ) .
Let u k ( t ) = ( u ( t ) , v k ) and f k = ( f , v k ) be Fourier coefficients of u ( t ) and f, respectively. Then, due to the self-adjointness of operator A, we obtain
D t ρ u k ( t ) = ( D t ρ u , v k ) = ( A u , v k ) + f k g ( t ) = ( u , A v k ) + f k g ( t ) = λ k u k ( t ) + f k g ( t ) .
Therefore, for u k ( t ) , t > 0 , we have the non-local boundary value problem
D t ρ u k ( t ) + λ k u k ( t ) = f k g ( t ) , t > 0 , u k ( T ) = u k ( 0 ) , k 1 ,
and the additional condition
u k ( t 0 ) = 0 .
For fixed f k , the unique solution of non-local boundary value problem (15) has the form (see Formula (4))
u k ( t ) = f k b k , ρ ( t ) + b k , ρ ( T ) 1 E ρ ( λ k T ρ ) E ρ ( λ k t ρ ) , k 1 .
Now, applying the additional condition in (16), we have
u k ( t 0 ) = f k ( 1 E ρ ( λ k T ρ ) ) b k , ρ ( t 0 ) + b k , ρ ( T ) E ρ ( λ k t 0 ρ ) 1 E ρ ( λ k T ρ ) = f k Δ ρ ( k , t 0 , T ) 1 E ρ ( λ k T ρ ) = 0 .
In accordance with Lemma 6, we have Δ ρ ( k , t 0 , T ) 0 for all k N . Moreover, Lemma 2 implies 1 E ρ ( λ k T ρ ) 0 for all k N . Therefore, Equation (18) is valid only if f k = 0 for all k N and due to the completeness of the set of eigenvectors { v k } in H, we obtain f = 0 . In turn, Equation (17) implies u k ( t ) 0 for all k N , and consequently, u ( t ) 0 . □
Now, let us consider the case when g ( t ) changes its sign. In this case, the function Δ ρ ( k , t 0 , T ) may become zero, and as a result, the set K 0 , ρ may be non-empty. As a consequence, in this case the solution of the inverse problem may not be unique. Consider the following example.
Example 1.
Let the inverse problem have the form
D t ρ u ( x , t ) u x x ( x , t ) = f ( x ) g ( t ) , ( x , t ) ( 0 , π ) × ( 0 , T ] , u ( 0 , t ) = u ( π , t ) = 0 , t [ 0 , T ] , u ( x , 0 ) = u ( x , T ) , x [ 0 , π ] , u ( x , T / 2 ) = 0 , x [ 0 , π ] ,
where 0 < ρ 1 , f L 2 ( 0 , π ) and g C [ 0 , T ] . First, note that this problem has a trivial solution ( u , f ) = ( 0 , 0 ) .
Assume T = 1 . One can easily verify that the pair of functions
u ( x , t ) = ω ( t ) sin ( x ) , f ( x ) = sin ( x ) ,
also satisfy problem (19) with
g ( t ) = D t ρ ω ( t ) + ω ( t ) ,
where ω ( t ) = ( t 1 2 ) 2 . The function g ( t ) is continuous on the interval [ 0 , 1 ] and can be written in the form
g ( t ) = t 2 ρ Γ ( 3 ρ ) t 1 ρ Γ ( 2 ρ ) + t 1 2 2 = t 1 2 2 2 ( t + ρ ) Γ ( 3 ρ ) t 1 ρ .
The reason for existing of a non-unique solution in this example is the function g ( t ) is not sign-preserving. Indeed, one has g ( 0 ) = 1 4 > 0 , and
g 1 2 = 3 2 ρ 2 1 ρ Γ ( 3 ρ ) < 0 ,
for all 0 < ρ 1 .
An important feature of the non-local problem being considered, which comes out from Theorem 1, is the relation of the eigenvalues of the operator A to the set of zeros of the function Δ ρ ( k , t 0 , T ) . In order to further analyze such a connection, we separately consider the cases of diffusion ( ρ = 1 ) and subdiffusion ( 0 < ρ < 1 ) equations.
Lemma 7.
Let ρ = 1 , g ( t ) C 1 [ 0 , T ] and g ( t 0 ) 0 . Then there exists a number k 0 such that for all k k 0 , the following estimate holds:
| Δ 1 ( k , t 0 , T ) |   C λ k ,
where constant C depends on k 0 , t 0 , and T.
Proof. 
By integrating by parts and utilizing the mean value theorem, we have
b k , 1 ( t 0 ) = 0 t 0 e λ k s g ( t 0 s ) d s = 1 λ k g ( t 0 s ) e λ k s | 0 t 0 1 λ k 0 t 0 e λ k s g ( t 0 s ) d s = 1 λ k g ( t 0 ) g ( 0 ) e λ k t 0 + g ( ξ k ) λ k 2 e λ k t 0 1 , ξ k [ 0 , t 0 ] .
Obviously, there exists a number k 0 such that for all k k 0 one has
| 1 2 g ( t 0 ) | | e λ k t 0 g ( 0 ) | + | g ( ξ k ) λ k e λ k t 0 1 | .
Thus
| b k , 1 ( t 0 ) | | g ( t 0 ) | 2 λ k .
Since function g ( t ) is bounded, then for any τ [ 0 , T ] we have
| b k , 1 ( τ ) | C 0 λ k .
Therefore
| Δ 1 ( k , t 0 , T ) |   =   | ( 1 e λ k T ) b k , 1 ( t 0 ) + e λ k t 0 b k , 1 ( T ) |
| b k , 1 ( t 0 ) |     | e λ k T b k , 1 ( t 0 ) |     | e λ k t 0 b k , 1 ( T ) |
| g ( t 0 ) | 2 λ k C 0 e λ k T λ k C 0 e λ k t 0 λ k .
Thus, for sufficiently large k, we obtain the statement of the lemma. □
Corollary 1.
Under conditions of Lemma 7, estimate (20) holds for all k K 1 .
Corollary 2.
Under conditions of Lemma 7, the set K 0 , 1 contains a finite number of elements.
In the case of subdiffusion equation ( ρ ( 0 , 1 ) ) the following assertion holds:
Lemma 8.
Let ρ ( 0 , 1 ) , g ( t ) C 1 [ 0 , T ] and g ( 0 ) 0 . Then, there exist numbers T 0 > 0 and k 0 such that, for all T T 0 and k k 0 , the following estimates holds:
| Δ ρ ( k , t 0 , T ) |   C λ k ,
where constant C depends on k 0 , t 0 , and T.
Proof. 
Let ρ ( 0 , 1 ) and τ ( 0 , T ] . Using equality (3), we integrate by parts and apply the mean value theorem. Then, we have
b k , ρ ( τ ) = 0 τ g ( τ s ) s ρ 1 E ρ , ρ ( λ k s ρ ) d s = 0 τ g ( τ s ) d s ρ E ρ , ρ + 1 ( λ k s ρ ) = g ( τ s ) s ρ E ρ , ρ + 1 ( λ k s ρ ) | 0 τ + 0 τ g ( τ s ) s ρ E ρ , ρ + 1 ( λ k s ρ ) d s = g ( 0 ) τ ρ E ρ , ρ + 1 ( λ k τ ρ ) + g ( ξ k ) 0 τ s ρ E ρ , ρ + 1 ( λ k s ρ ) d s , ξ k [ 0 , τ ] .
Using Formula (3), the integral in (24) can be represented as
0 τ s ρ E ρ , ρ + 1 ( λ k s ρ ) d s = τ ρ + 1 E ρ , ρ + 2 ( λ k τ ρ ) .
The latter due to the asymptotic estimate of the Mittag–Leffler functions (Lemma 3) implies
b k , ρ ( τ ) = g ( 0 ) λ k + g ( ξ k ) λ k τ + O 1 ( λ k t 0 ρ ) 2 .
Now, without loss of generality, we assume that g ( 0 ) > 0 . Then, for sufficiently small τ and sufficiently large k we obtain the lower estimate
b k , ρ ( τ ) g ( 0 ) 2 λ k .
Therefore
Δ ρ ( k , t 0 , T ) ( 1 E ρ ( λ k T ρ ) ) g ( 0 ) 2 λ k + E ρ ( λ k t 0 ρ ) g ( 0 ) 2 λ k ,
and since the classical Mittag–Leffler function of the negative argument is a monotonically decreasing function (see Lemma 2), then the desired result follows. □
Corollary 3.
Under conditions of Lemma 8, estimate (23) holds for all T T 0 and k K ρ .
Corollary 4.
Under conditions of Lemma 8, if T is sufficiently small, then the set K 0 , ρ has a finite number elements.
Theorem 2.
Let g ( t ) C 1 [ 0 , T ] , ψ D ( A ) . Further, we assume that for ρ = 1 the conditions of Lemma 7 are satisfied, t 0 ( 0 , T ) is chosen correspondingly and for ρ ( 0 , 1 ) , the conditions of Lemma 8 are satisfied and T is sufficiently small.
(1) If the set K 0 , ρ is empty, i.e., Δ ρ ( k , t 0 , T ) 0 , for all k, then there exists a unique solution of the inverse problem (1) and (2):
f = k = 1 f k v k ,
u ( t ) = k = 1 f k b k , ρ ( t ) + b k , ρ ( T ) 1 E ρ ( λ k T ρ ) E ρ ( λ k t ρ ) v k ,
where
f k = ψ k ( 1 E ρ ( λ k T ρ ) ) Δ ρ ( k , t 0 , T ) .
(2) If the set K 0 , ρ is not empty, then for the existence of a solution to the inverse problem, it is necessary and sufficient that the following conditions
ψ k = ( ψ , v k ) = 0 , k K 0 , ρ ,
be satisfied. In this case, the solution to the problem (1) and (2) exists, but is not unique. In this case the solution has the representation
f = k K ρ ψ k ( 1 E ρ ( λ k T ρ ) ) Δ ρ ( k , t 0 , T ) v k + k K 0 , ρ f k v k ,
u ( t ) = k = 1 f k b k , ρ ( t ) + b k , ρ ( T ) 1 E ρ ( λ k T ρ ) E ρ ( λ k t ρ ) v k ,
where f k , k K 0 , ρ , are arbitrary real numbers.
Proof. 
The proof of the first part of the theorem is completely analogous to the proof of Theorem 1.
In what concerns the proof of the second part of the theorem, we note the following. If k K ρ , then again it follows from (5) that the relations in (7) and (8) hold. Then the convergence of the corresponding series can be proved in the same way as in Theorem 1.
If k K 0 , ρ , i.e., Δ ρ ( k , t 0 , T ) = 0 , then the solution of Equation (5) with respect to f k exists if and only if condition (25) is satisfied. In this case, the solution of the equation contains arbitrary numbers f k . As shown above (see Corollaries 2 and 4), under the conditions of the theorem, the set K 0 , ρ , ρ ( 0 , 1 ] , contains only a finite number of elements. □
Note that we considered an abstract subdiffusion equation with a self-adjoint operator A in a Hilbert space H. The choice of an abstract operator allows us to generalize various known models. Since the operator A is only required to have a complete orthonormal system of eigenfunctions, then as A one can consider any of the elliptic operators given in the work by Ruzhansky et al. [40]. For example, let us take L 2 ( Ω ) , Ω R N , as the Hilbert space H and let A be the Laplace operator ( Δ ) with the Dirichlet condition. Assume that the boundary Ω of Ω is sufficiently smooth. Then, problem (1) takes the form
D t ρ u ( x , t ) Δ u ( x , t ) = f ( x ) g ( t ) , ( x , t ) Ω × ( 0 , T ] , u ( x , t ) = 0 , x Ω , t [ 0 , T ] , u ( x , 0 ) = u ( x , T ) , x Ω , u ( x , t 0 ) = ψ ( x ) , x Ω , t 0 ( 0 , T ) ,
where g C [ 0 , T ] and ψ L 2 ( Ω ) are given functions and f is an unknown function in L 2 ( Ω ) .
Let v k ( x ) and λ k be eigenfunctions and eigenvalues of the spectral problem:
Δ v ( x ) = λ v ( x ) , v ( x ) = 0 , x Ω .
Assume that all the conditions of Theorem 2 (or Theorem 1) are satisfied and ψ belongs to the Sobolev space W 2 2 ( Ω ) and ψ ( x ) = 0 , x Ω . Then, ψ D ( A ) and the statements of Theorem 2 (or, respectively, Theorem 1) are valid.

4. Conclusions

The work is devoted to one of the important applications of the inverse problem—the problem of determining the source function in subdiffusion equations. The unknown right-hand side is given in the form g ( t ) f , where f H is unknown. In contrast to the well-known works, the non-local time condition is taken instead of the Cauchy condition. If g ( t ) is sign-preserving, then the existence and uniqueness of the solution of the inverse problem with the over-determination condition u ( t 0 ) = ψ is proved with t 0 ( 0 , T ) . An example is constructed showing the lack of uniqueness of the solution if g ( t ) is not sign-preserving. The meaning of taking t 0 ( 0 , T ) instead of T in the over-determined condition is that for some sign-changing functions g ( t ) , the choice of t 0 (as it is shown in Theorem 2) can ensure the uniqueness of the solution. Note that if t 0 = T , then in Theorem 2 there will be no unique solution to the inverse problem. In the absence of uniqueness of the solution, the conditions for the orthogonality of the element ψ to some eigenfunctions of the operator A are found, the fulfillment of which ensures the existence of a solution to the inverse problem. It should be noted that in this case, the solution is not unique.

Author Contributions

Conceptualization, R.A., M.S. and S.U. 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

The data presented in this study are available upon request from the corresponding author.

Acknowledgments

We are grateful to Sh. Alimov for discussions of the obtained results.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Machado, J.A.T. (Ed.) Handbook of Fractional Calculus with Applications; DeGruyter: Berlin, Germany, 2019; pp. 1–8. [Google Scholar]
  2. Kilbas, A.; Srivastava, H.; Trujillo, J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  3. Uchaikin, V.V. Fractional Derivatives for Physicists and Engineers, 1, Background and Theory, 2, Application; Springer: Berlin, Germany, 2013. [Google Scholar]
  4. Umarov, S.R.; Hahn, M.; Kobayashi, K. Beyond the Triangle: Brownian Motion, Itô Calculas, and Fokker-Plank Equation—Fractional Generalizations; World Scientific: Singapore, 2018. [Google Scholar]
  5. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogozin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
  6. Umarov, S.R. Introduction to Fractional and Pseudo-Differential Equations with Singular Symbols; Springer: Berlin, Germany, 2015. [Google Scholar]
  7. Hilfer, R.; Luchko, Y.; Tomovski, Z. Operational method for solution of the fractional differential equations with the generalized Riemann-Liouville fractional derivatives. Frac. Calculus Appl. Anal. 2009, 12, 299–318. [Google Scholar]
  8. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  9. Pskhu, A.V. Fractional Differential Equations; Nauka: Moscow, Russia, 2005. (In Russian) [Google Scholar]
  10. Lizama, C. Abstract linear fractional evolution equations. In Handbook of Fractional Calculus with Applications; Marchado, J.A.T., Ed.; De Gruyter: Berlin, Germany, 2019; Volume 2, pp. 465–497. [Google Scholar]
  11. Ashurov, R.; Fayziev, Y. On the non-local problems in time for time-fractional subdiffusion equations. Fractal Fract. 2022, 6, 41. [Google Scholar] [CrossRef]
  12. Ruzhansky, M.; Tokmagambetov, N.; Torabek, B.T. On a non-local problem for a multi-term fractional diffusion-wave equation. Fract. Calculus Appl. Anal. 2020, 23, 324–355. [Google Scholar] [CrossRef]
  13. Prilepko, A.I.; Kostin, A.B. On certain inverse problems for parabolic equations with final and integral observation. Mat. Sb. 1992, 183, 49–68. [Google Scholar] [CrossRef]
  14. Sabitov, K.B.; Zaynullov, A.R. On the theory of the known inverse problems for the heat transfer equation. Ser. Phys. Math. Sci. 2019, 161, 274–291. [Google Scholar] [CrossRef]
  15. Orlovskii, D.G. On a problem of determining the parameter of an evolution equation. Differ. Uravn. 1990, 26, 1614–1621. [Google Scholar]
  16. Tikhonov, I.V.; Éidelman, Y.S. Uniqueness criterion in an inverse problem for an abstract differential equation with nonstationary inhomogeneous term. Mat. Notes 2005, 77, 273–290. [Google Scholar] [CrossRef]
  17. Bockstal, K.V.; Slodicka, M. Recovery of a space-dependent vector source in thermoelastic systems. Inverse Probl. Sci. Eng. 2015, 23, 956–968. [Google Scholar] [CrossRef]
  18. Slodicka, M. Uniqueness for an inverse source problem of determining a space-dependent source in a non-autonomous time-fractional diffusion equation. Fract. Calc. Appl. Anal. 2020, 23, 1703–1711. [Google Scholar] [CrossRef]
  19. Slodicka, M.; Siskova, K.; Bockstal, V. Uniqueness for an inverse source problem of determining a space dependent source in a time-fractional diffusion equation. Appl. Math. Lett. 2019, 91, 15–21. [Google Scholar] [CrossRef]
  20. Slodicka, M. Some direct and inverse source problems in nonlinear evolutionary PDEs with Volterra operators. Inverse Probl. 2022, 38, 12. [Google Scholar] [CrossRef]
  21. Bockstal, K.V. Uniqueness for inverse source problems of determining a space-dependent source in time-fractional equations with non-smooth solutions. Fractal Fract. 2021, 5, 169. [Google Scholar] [CrossRef]
  22. Ashurov, R.; Shakarova, M. Inverse problem for the subdiffusion equation with fractional Caputo derivative. Ufimskiy Mat. J. 2024, 16, 111–125. [Google Scholar]
  23. Liu, Y.; Li, Z.; Yamamoto, M. Inverse problems of determining sources of the fractional partial differential equations. In Handbook of Fractional Calculus with Applications; Marchado, J.A.T., Ed.; De Gruyter: Berlin, Germany, 2019; Volume 2, pp. 411–430. [Google Scholar]
  24. Ashurov, R.; Shakarova, M. Time-dependent source identification problem for fractional Schrödinger type equations. Lobachevskii J. Math. 2022, 42, 517–525. [Google Scholar] [CrossRef]
  25. Ashurov, R.; Shakarova, M. Time-dependent source identification problem for a fractional Schrödinger equation with the Riemann-Liouville derivative. Ukr. Math. J. 2023, 75, 871–887. [Google Scholar] [CrossRef]
  26. Sakamoto, K.; Yamamoto, M. Initial value boundary value problems for fractional diffusion-wave equations and applications to some inverse problems. J. Math. Anal. Appl. 2011, 382, 426–447. [Google Scholar] [CrossRef]
  27. Furati, K.; Iyiola, O.; Kirane, M. An inverse problem for a generalized fractional diffusion. Appl. Math. Comput. 2014, 249, 24–31. [Google Scholar] [CrossRef]
  28. Ashurov, R.R.; Mukhiddinova, A.T. Inverse Problem of Determining the Heat Source Density for the Subdiffusion Equation. Differ. Equ. 2020, 56, 1550–1563. [Google Scholar] [CrossRef]
  29. Kirane, M.; Malik, 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]
  30. Kirane, M.; Samet, B.; Torebek, B. Determination of an unknown source term and the temperature distribution for the subdiffusion equation at the initial and final data. Electron. J. Differ. Equ. 2017, 217, 1–13. [Google Scholar]
  31. Li, Z.; Liu, Y.; Yamamoto, M. Initial-boundary value problem for multi-term time-fractional diffusion equation with positive constant coefficients. Appl. Math. Comput. 2015, 257, 381–397. [Google Scholar] [CrossRef]
  32. Malik, S.; Aziz, S. An inverse source problem for a two parameter anomalous diffusion equation with non-local boundary conditions. Comput. Math. Appl. 2017, 3, 7–19. [Google Scholar]
  33. Ashurov, R.R.; Fayziev, Y.E. Determination of fractional order and source term in a fractional subdiffusion equation. Eurasian Math. J. 2022, 13, 19–31. [Google Scholar] [CrossRef]
  34. Kabanikhin, S.I. Inverse and Ill-Posed Problems. In Theory and Applications; De Gruyter: Berlin, Germany, 2011. [Google Scholar]
  35. Sabitov, K.B.; Zaynullov, A.R. Inverse problems for a two-dimensional heat equation with unknown right-hand side. Russ. Math. 2021, 3, 75–88. [Google Scholar] [CrossRef]
  36. Fedorov, V.E.; Nagumanova, A.V. Inverse problem for evolutionary equation with the Gerasimov-Caputo fractional derivative in the sectorial case. Bull. Irkutsk. State Univ. Ser. Math. 2019, 28, 123–137. [Google Scholar] [CrossRef]
  37. Liu, S.; Sun, F.; Feng, L. Regularization of inverse source problem for fractional diffusion equation with Riemann-Liouville derivative. Comput. Appl. Math. 2021, 40, 112. [Google Scholar] [CrossRef]
  38. Niu, P.; Helin, T.; Zhang, Z. An inverse random source problem in a stochastic fractional diffusion equation. Inverse Probl. 2020, 36, 045002. [Google Scholar] [CrossRef]
  39. Dzherbashian, M.M. Integral Transforms and Representation of Functions in the Complex Domain; Nauka: Moscow, Russia, 1966. (In Russian) [Google Scholar]
  40. Ruzhansky, M.; Tokmagambetov, N.; Torebek, B.T. Inverse source problems for positive operators. I: Hypoelliptic Diffus. Subdiffusion Equations. J. Inverse Ill-Possed Probl. 2019, 27, 891–911. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ashurov, R.; Shakarova, M.; Umarov, S. An Inverse Problem for the Subdiffusion Equation with a Non-Local in Time Condition. Fractal Fract. 2024, 8, 378. https://doi.org/10.3390/fractalfract8070378

AMA Style

Ashurov R, Shakarova M, Umarov S. An Inverse Problem for the Subdiffusion Equation with a Non-Local in Time Condition. Fractal and Fractional. 2024; 8(7):378. https://doi.org/10.3390/fractalfract8070378

Chicago/Turabian Style

Ashurov, Ravshan, Marjona Shakarova, and Sabir Umarov. 2024. "An Inverse Problem for the Subdiffusion Equation with a Non-Local in Time Condition" Fractal and Fractional 8, no. 7: 378. https://doi.org/10.3390/fractalfract8070378

Article Metrics

Back to TopTop