Next Article in Journal
Numerical Solution of Nonlinear Fractional Diffusion Equation in Framework of the Yang–Abdel–Cattani Derivative Operator
Next Article in Special Issue
Quantum Walks in Hilbert Space of Lévy Matrices: Recurrences and Revivals
Previous Article in Journal
Alternating Inertial and Overrelaxed Algorithms for Distributed Generalized Nash Equilibrium Seeking in Multi-Player Games
Previous Article in Special Issue
Subdiffusive Reaction Model of Molecular Species in Liquid Layers: Fractional Reaction-Telegraph Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Inverse Source Problem for the Generalized Subdiffusion Equation with Nonclassical Boundary Conditions

by
Emilia Bazhlekova
Institute of Mathematics and Informatics, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Bl. 8, 1113 Sofia, Bulgaria
Submission received: 30 April 2021 / Revised: 24 June 2021 / Accepted: 26 June 2021 / Published: 30 June 2021
(This article belongs to the Special Issue Fractional Dynamics: Theory and Applications)

Abstract

:
An initial-boundary-value problem is considered for the one-dimensional diffusion equation with a general convolutional derivative in time and nonclassical boundary conditions. We are concerned with the inverse source problem of recovery of a space-dependent source term from given final time data. Generalized eigenfunction expansions are used with respect to a biorthogonal pair of bases. Existence, uniqueness and stability estimates in Sobolev spaces are established.

1. Introduction

Anomalous (non-Fickian) diffusion processes are modeled by employing different types of fractional partial differential and integro-differential equations [1]. Analytical and numerical solutions to fractional differential equations have been investigated by several authors, see e.g., [2,3]. Recently, generalized subdiffusion equations with different memory kernels are extensively studied as models that unify a wide range of anomalous diffusion patterns [4,5,6].
In this work we consider the generalized subdiffusion equation in the form
D t ( k ) u ( x , t ) = u x x ( x , t ) + F ( x , t ) ,
where D t ( k ) is an integro-differential operator defined by
( D t ( k ) f ) ( t ) = d d t 0 t k ( t τ ) ( f ( τ ) f ( 0 ) ) d τ , t > 0 ,
with a locally integrable memory kernel k ( t ) . For the kernel k ( t ) we also assume that its Laplace transform k ^ ( s ) exists for all s > 0 and
k ^ ( s ) SF , lim s + s k ^ ( s ) = + ,
where SF denotes the class of Stieltjes functions (the definition of this class is given in the next section). The most prominent particular examples of operators D t ( k ) are the first-order derivative d d t , corresponding to k ^ ( s ) = 1 , the Caputo time-fractional derivative of order α ( 0 , 1 ) , where k ^ ( s ) = s α 1 , as well as linear combinations with positive coefficients of such derivatives.
Let us note that assumptions (3) are weaker than those required in the definition of the so-called general fractional derivative, introduced in [7] and studied in detail in [8,9]. The operator D t ( k ) is a general fractional derivative, if, along with (3), the following additional limiting behavior conditions are imposed: k ^ ( s ) 0 as s ; k ^ ( s ) and s k ^ ( s ) 0 as s 0 . To cover some examples of physically meaningful models with corresponding memory kernels, which do not satisfy some of the additional conditions, they are not required in this work. Such are, for instance, the subdiffusion equation with the truncated power-law memory kernel k ( t ) = e γ t t α Γ ( 1 α ) , γ > 0 , α ( 0 , 1 ) , considered in [4,6], the fractal mobile/immobile solute transport equation introduced in [10], and the Jeffreys’ type heat conduction model in the diffusion regime [11,12]. On the other hand, the assumption k ^ ( s ) SF is typical for a subdiffusion model (see e.g., [4,5,6]) and allows the use of the convenient Bernstein functions technique [13]. It implies that the kernel k ( t ) admits the representation
k ( t ) = k 0 δ ( t ) + k 1 ( t ) ,
where k 0 0 , δ ( · ) denotes the Dirac delta function, and k 1 ( t ) L l o c 1 ( R + ) is a completely monotone function, i.e. it is of class C ( R + ) and
( 1 ) n k 1 ( n ) ( t ) 0 , t > 0 , n = 0 , 1 , 2 ,
For example, in the case of the first-order derivative k 0 = 1 and k 1 0 , while k 0 = 0 and k 1 = t α Γ ( 1 α ) for the Caputo time-fractional derivative of order α ( 0 , 1 ) .
A more general setting for the subdiffusion equation (1) has been introduced in [14] and further developed and applied in e.g., [15,16]. In this setting, it is assumed that the function k 1 ( t ) in representation (4) is non-negative and non-increasing (satisfies (5) only for n = 0 , 1 ) and there exists a locally integrable kernel l ( t ) , such that ( k l ) ( t ) = 1 , t > 0 . The kernel l ( t ) with these properties is referred to as completely positive kernel. Let us note that assumptions (3) also ensure the existence of such kernel l ( t ) , which in this case is completely monotone.
In the present work we are concerned with an inverse source problem for the subdiffusion Equation (1). Different kinds of inverse problems for the diffusion equation with the Caputo time-derivative of order α ( 0 , 1 ) are extensively studied recently, see e.g., [17,18,19,20,21]. For a comprehensive tutorial on inverse problems for anomalous diffusion processes we refer to [22]. Identification of a space-dependent source factor h ( x ) in a source function of the form F ( x , t ) = h ( x ) q ( x , t ) from final overdetermination are studied in [17,19,23,24,25], where different assumptions on the known source factor q ( x , t ) are discussed. Concerning the generalized subdiffusion equation, various types of inverse problems for such equations are studied in [26,27,28].
In this work, we consider the problem of identifying a space-dependent source factor h ( x ) and the solution u ( x , t ) to the following nonlocal boundary-value problem with final overdetermination
D t ( k ) u ( x , t ) = u x x ( x , t ) + q ( t ) h ( x ) , x ( 0 , 1 ) , t ( 0 , T ) ,
u ( 1 , t ) = 0 , u x ( 0 , t ) = u x ( 1 , t ) , t ( 0 , T ] ,
u ( x , 0 ) = 0 , u ( x , T ) = g ( x ) , x [ 0 , 1 ] ,
where the operator D t ( k ) acting with respect to the time variable is defined in (2), q ( t ) is a prescribed continuous function, g ( x ) is a known square integrable function, and T > 0 is the final time.
In practical applications, the input data g ( x ) is given by measurement and actually the measured data g ν ( x ) is available, which is merely in L 2 ( 0 , 1 ) and satisfies
g ν g L 2 ( 0 , 1 ) ν ,
where the constant ν represents the noise level.
The study of nonlocal boundary-value problems is motivated by the fact that in many cases a nonlocal condition is more realistic in treating physical problems than the classical local conditions. Inverse source problems with nonlocal boundary conditions are studied e.g., in [29,30,31,32,33,34]. The papers [29,30] are concerned with two particular cases of the inverse source problem (6)–(8). In these works, existence of a unique solution in the classical sense is established for D t ( k ) = d d t and D t ( k ) —the Caputo time-fractional derivative of order α ( 0 , 1 ) , respectively, and q ( t ) = 1 .
In the case of inverse source problem with final overdetermination for the time-fractional subdiffusion equation on Ω × ( 0 , T ) with the classical Dirichlet boundary conditions the following estimates are satisfied (see e.g., [19,22,35])
C 1 g H 2 ( Ω ) h L 2 ( Ω ) C 2 g H 2 ( Ω ) .
The same behavior is observed as well with the classical diffusion equation. Therefore, such inverse problems are moderately ill-posed: the overdetermination function g ( x ) has a better regularity than the source term h ( x ) , as the regularity loss is two spatial derivatives.
The main goal of the present work is to study stability of the inverse source problem (6)–(8) in Sobolev spaces, and to prove analogous estimates for the overdetermination function g ( x ) and the source factor h ( x ) . To this end generalized eigenfunction expansions are used with respect to a biorthogonal pair of bases. The present paper is a continuation of the recent work [36], in which inverse source problem (6)–(8) is solved in the classical sense, together with detailed study of particular cases for the operator D t ( k ) and some numerical examples.
The rest of the paper is organized as follows. In Section 2 the assumptions on the memory kernel are discussed by the use of the Bernstein functions technique. In Section 3 two examples of generalized subdiffusion equations are considered. Properties of the solution of the generalized relaxation equation are summarized in Section 4. Section 5 is concerned a biorthonormal pair of Riesz bases for the considered problem. In Section 6 formal spectral expansions for the solution are derived. In Section 7 we establish uniqueness and stability estimates for the inverse problem. Concluding remarks are given in Section 8.

2. Assumptions on the Memory Kernel

We start with some preliminaries on completely monotone functions, Bernstein, complete Bernstein and Stieltjes functions. Let us denote the Laplace transform of a function by
L { f ( t ) } ( s ) = f ^ ( s ) = 0 e s t f ( t ) d t .
The class of completely monotone functions, defined in (5), is denoted by CMF . The characterization of the class CMF is given by the Bernstein’s theorem stating that a function is completely monotone if and only if it can be represented as the Laplace transform of a non-negative measure (non-negative function or generalized function).
The class of Stieltjes functions ( SF ) consists of all functions defined on R + which can be written as a restriction of the Laplace transform of a completely monotone function to the real positive semi-axis. More precisely, φ SF if and only if it admits the representation [7]
φ ( s ) = φ 0 + 0 e s t ψ ( t ) d t , s > 0 ,
where φ 0 0 , ψ CMF and the Laplace transform of ψ exists for any s > 0 . Moreover, φ 0 = lim s + φ ( s ) , see e.g., [37], Theorem 2.6.
A non-negative function φ on R + is said to be a Bernstein function ( φ BF ) if φ ( s ) CMF ; φ ( s ) is said to be a complete Bernstein function ( CBF ) if and only if
φ ( s ) / s SF , s > 0 .
The inclusions SF CMF and CBF BF are valid. An example of a completely monotone function is the Mittag–Leffler function E α , β ( s ) (for definition see (16)) provided 0 < α 1 and α β . If α [ 0 , 1 ] then s α CBF and s α 1 SF .
A selection of properties of the classes of functions defined above is given next. For proofs and more details on these special classes of functions we refer to [13,15], Chapter 4.
(P1)
The class CMF is closed under point-wise multiplication.
(P2)
If φ CMF and ψ BF then the composite function φ ( ψ ) CMF .
(P3)
φ CBF if and only if 1 / φ SF .
(P4)
If φ BF then it admits a continuous extension to C ¯ + , which is holomorphic in C + and satisfies φ ( s ) > 0 for all s > 0 .
We proceed with a short discussion on the assumptions (3) for the kernel k ( t ) and their implications.
Let us first note that the assumption k ^ ( s ) SF implies the non-negativity of the Green function G ( x , t ) to Equation (1), which is a necessary condition for a diffusion model. Indeed, in Laplace domain the Green function obeys the identity
G ^ ( x , s ) = k ^ ( s ) 2 exp | x | s k ^ ( s ) , x R , s > 0 .
By the Bernstein’s theorem, it is sufficient to prove that G ^ ( x , s ) CMF for s > 0 and | x | considered to be a parameter. According to (11), k ^ ( s ) SF is equivalent to s k ^ ( s ) CBF . Then (P2) implies that exp | x | s k ^ ( s ) CMF as a composition of the completely monotone in t function e | x | t and the Bernstein function s k ^ ( s ) . Moreover, k ^ ( s ) SF CMF . Therefore, G ^ ( x , s ) CMF as a product of two completely monotone functions, see (P1).
Concerning the limiting behavior of the kernel k ( t ) the initial value theorem for the Laplace transform implies lim t 0 + k ( t ) = lim s s k ^ ( s ) = + . Thus, we restrict our attention to kernels, singular at the origin.
The representation (4) for the kernel k ( t ) follows from the assumption k ^ ( s ) SF and the characterization (10) of Stieltjes functions.
Furthermore, we are looking for a corresponding Sonine kernel l ( t ) to the kernel k ( t ) . This is a function satisfying l ( t ) L l o c 1 ( R + ) , such that
( k l ) ( t ) 1 .
Here ∗ denotes the convolution ( k l ) ( t ) = 0 t k ( t τ ) l ( τ ) d τ . Relation (12) is equivalent to k ^ ( s ) l ^ ( s ) = 1 / s and by the use of (P3) assumptions (3) imply l ^ ( s ) SF and lim s + l ^ ( s ) = 0 . Therefore, according to (10), under the assumptions (3) a resolvent kernel l ( t ) exists and l ( t ) CMF .
Let us note that if f is integrable then ( k f ) ( t ) = ( k f ) ( t ) + k ( t ) f ( 0 ) . Therefore, in this case ( D t ( k ) f ) ( t ) = ( k f ) ( t ) and (12) implies that ( D t ( k ) f ) ( t ) = g ( t ) is equivalent to f ( t ) = ( l g ) ( t ) . This enables us to rewrite Equation (1) as a generalized diffusion equation in the so-called modified form, see e.g., [4]. In [38] we presented arguments why this form is more natural when considering multi-component systems, in particular, systems with a source. Assume that there is a source term h ( x ) in the conservation law, which depends only on the space variable. Then the subdiffusion equation in modified form reads
u t ( x , t ) = t l u x x ( x , t ) + h ( x ) .
Applying the convolution operator ( k ) to both sides of (13) we obtain
k u t = k t l u x x + ( k 1 ) ( t ) h ( x ) = t k l u x x k ( t ) l u x x ( x , 0 ) + ( k 1 ) ( t ) h ( x ) ,
which implies by the use of (12)
D t ( k ) u = u x x ( x , t ) + q ( t ) h ( x ) ,
where q ( t ) = 0 t k ( τ ) d τ . Therefore, q ( t ) is a Bernstein function with q ( 0 ) = k 0 0 . In particular, it is continuous, non-negative and non-decreasing function.

3. Two Examples of Generalized Subdiffusion Equations

In [36] we considered basic examples of kernels k ( t ) , satisfying (3), which do not contain a Dirac delta function in their definition, i.e., k 0 = 0 in (4). In this section, two additional physically meaningful examples are given, such that k 0 > 0 .
The two models are defined in terms of derivatives of fractional order and kernels of Mittag–Leffler type.
The Riemann-Liouville and the Caputo fractional derivatives of order α ( 0 , 1 ) , denoted by D t α and C D t α , respectively, are defined by the identities [39]
D t α f ( t ) = d d t J t 1 α f ( t ) , C D t α f ( t ) = d d t J t 1 α ( f ( · ) f ( 0 ) ) ( t ) , 0 < α < 1 ,
where J t β denotes the Riemann-Liouville fractional integral:
J t β f ( t ) = 1 Γ ( β ) 0 t ( t τ ) β 1 f ( τ ) d τ , β > 0 ; J t 0 = I .
The Mittag–Leffler function is defined by the series [39]
E α , β ( z ) = k = 0 z k Γ ( α k + β ) α , β , z C , α > 0 ; E α ( z ) = E α , 1 ( z ) ;
and satisfies the Laplace transform pair
L t β 1 E α , β ( λ t α ) ( s ) = s α β s α + λ .

3.1. Two Time-Scale Diffusion Model

The time-fractional diffusion equation with the fractional Caputo time-derivative of order α ( 0 , 1 ) accurately describes the power-law decaying behavior of the subdiffusive transport of solutes in heterogeneous media. However, the first-order time derivative of the solution to this equation exhibits a nonphysical singularity of order O ( t α 1 ) near the initial time t = 0 . This disadvantage can be avoided if the following two time-scale model is considered
u t ( x , t ) + a C D t α u ( x , t ) = u x x ( x , t ) + F ( x , t ) , α ( 0 , 1 ) .
Equation (18) is introduced in [10] as a model for anomalously diffusive transport of solute in heterogeneous porous media, where it is referred to as Fractal mobile/immobile solute transport equation.
In the case of Equation (18) the corresponding Sonine kernels k ( t ) and l ( t ) satisfy the identities
k ( t ) = δ ( t ) + a t α Γ ( 1 α ) , k ^ ( s ) = 1 + a s α 1 , l ( t ) = E 1 α ( a t 1 α ) , l ^ ( s ) = s α s 1 α + a .

3.2. Fractional Jeffreys-Type Heat Conduction Equation

According to the fractional Jeffreys-type constitutive law (see [11] Chapter 7) the heat flux q ( x , t ) and the temperature u ( x , t ) are related via the equation
1 + a D t α q ( x , t ) = k 1 + b D t α u x ( x , t ) ,
where D t α denotes the Riemann-Liouville fractional time-derivative of order α ( 0 , 1 ] , a > 0 and b > 0 are generalized relaxation times, k is the thermal conductivity. Combining Equation (19) with the energy balance equation
x q ( x , t ) + J ( x , t ) = C u t ( x , t ) ,
where C denotes the heat capacity and J ( x , t ) is a heat source, one finds the following fractional Jeffreys-type heat conduction equation
1 + a D t α u t ( x , t ) = D 1 + b D t α u x x ( x , t ) + C 1 1 + a D t α J ( x , t ) ,
where D = k / C is the thermal diffusivity. For notational simplicity we set D = 1 . Suppose a heat source, which does not change with time and set h ( x ) = C 1 J ( x ) .
Let us assume in addition 0 < a < b . It is proven in [12] that in this case Equation (21) is a generalized subdiffusion equation and admits the form (6) with kernel (see [12])
k ( t ) = a b δ ( t ) + 1 a b 1 b t α 1 E α , α 1 b t α ; k ^ ( s ) = 1 + a s α 1 + b s α .
By the use of (12) and (17) one finds that the corresponding Sonine kernel l ( t ) obeys
l ( t ) = 1 + b a 1 E α 1 a t α ; l ^ ( s ) = 1 + b s α s ( 1 + a s α ) .
For the source factor q ( t ) we derive from (22) q ( t ) = 1 + a b 1 E α 1 b t α . Therefore, in this case q ( t ) is a continuous, positive and increasing function with q ( 0 ) = a / b > 0 .

4. Generalized Relaxation Equation

This section contains preliminary results on the inhomogeneous generalized relaxation equation, for more details we refer to [35,36].
The solution of the inhomogeneous generalized relaxation equation
( D t ( k ) v ) ( t ) + λ v ( t ) = r ( t ) , λ 0 , t > 0 , v ( 0 ) = a ,
admits the representation
v ( t ) = a S ( t ; λ ) + 0 t G ( t τ ; λ ) r ( τ ) d τ ,
where the functions S ( t ; λ ) and G ( t ; λ ) are defined in Laplace domain as follows
S ^ ( s ; λ ) = k ^ ( s ) s k ^ ( s ) + λ G ^ ( s ; λ ) = 1 s k ^ ( s ) + λ λ 0 .
In the particular case λ = 0
S ( t ; 0 ) 1 , G ( t ; 0 ) = l ( t ) , t 0 ,
where l ( t ) is the corresponding Sonine kernel to the kernel k ( t ) , see (12).
Next, we summarize some properties of the functions S ( t ; λ ) and G ( t ; λ ) , see e.g., [35], Theorem 2.2 and [36], Lemma 4.1, and [40], Theorem 5.2.
Proposition 1.
Assume the kernel k satisfies conditions (3) and λ 0 . Then Equation (24) has a unique solution given by (25), where the functions S ( t ; λ ) and G ( t ; λ ) are locally integrable and completely monotone with respect to t R + , and admit holomorphic extensions to the half-plane C + . Moreover,
d d t S ( t ; λ ) = λ G ( t ; λ ) , t 0 ;
and for any λ λ * > 0
S ( t ; λ ) S ( t ; λ * ) , t 0 ;
and
c λ 0 T G ( t ; λ ) d t < 1 , T > 0 ,
where the constant c = 1 S ( T ; λ * ) > 0 is independent of λ.
Let us note that for any λ 0 and t > 0
0 t G ( τ ; λ ) d τ > 0 .
For λ > 0 this follows from (30). If λ = 0 then (31) reduces to ( 1 l ) ( t ) > 0 , see (27). Since l ( t ) CMF then ( 1 l ) ( t ) BF . This implies that it is positive for all t > 0 , as a continuous, non-negative and non-decreasing function for t 0 , which is analytic for t > 0 , see property (P4) in Section 2.
In the particular case of the relaxation equation with the Caputo fractional derivative of order α ( 0 , 1 ) the representation are valid S ( t ; λ ) = E α ( λ t α ) and G ( t ; λ ) = t α 1 E α , α ( λ t α ) , where E α , β ( · ) is the Mittag–Leffler function. In this case, (29) is the well-known relation
d d t E α λ t α = λ t α 1 E α , α λ t α
and estimates (30) are used in [17,19,20,21] for the proof of uniqueness or/and (conditional) stability results of inverse problems for time-fractional diffusion equations.
Inequalities (30) and (31) are used in [35,36] for the study of inverse source problems for generalized subdiffusion equations, based on eigenfunction expansion.

5. Biorthonormal Pair of Riesz Bases

Let us denote by . , . the inner product in L 2 ( 0 , 1 ) , i.e., f , g = 0 1 f ( x ) g ( x ) d x . The norm . L 2 ( 0 , 1 ) in L 2 ( 0 , 1 ) is f L 2 ( 0 , 1 ) = f , f 1 / 2 . The Sobolev space H 2 ( 0 , 1 ) is defined as the subset of functions f L 2 ( 0 , 1 ) , such that f and its weak derivatives f , f have a finite L 2 norm. Consider the space H 2 ( 0 , 1 ) equipped with the norm f H 2 ( 0 , 1 ) = f L 2 ( 0 , 1 ) + f L 2 ( 0 , 1 ) .
We apply the technique of spectral decomposition with respect to generalized eigenfunctions for the non-selfadjoint operator defined by the boundary conditions (7), see e.g., [41,42,43]. Since the eigenvalues of the spectral problem for the second order differential operator with the boundary conditions (7) are 4 π 2 n 2 , n N 0 , and for n 0 each eigenvalue has multiplicity 2, the system of eigenfunctions is not complete and must be supplemented with associated functions. In this way, the following system of generalized eigenfunctions (eigenfunctions and associated eigenfunctions) { X 0 , n , X 1 , n , X 2 , n } n = 1 is obtained:
X 1 , 0 ( x ) = 2 ( 1 x ) , X 1 , n ( x ) = 4 ( 1 x ) cos λ n x , X 2 , n ( x ) = 4 sin λ n x , λ n = 2 π n , n N .
This system of functions is a basis in L 2 ( 0 , 1 ) . Any function f L 2 ( 0 , 1 ) admits the following unique formal spectral expansion
f ( x ) = f 1 , 0 X 1 , 0 ( x ) + n = 1 f 1 , n X 1 , n ( x ) + f 2 , n X 2 , n ( x )
with coefficients defined by the identities
f j , n = f , Y j , n ,
where the system of functions { Y 0 , n , Y 1 , n , Y 2 , n } n = 1 is
Y 1 , 0 ( x ) = 1 , Y 1 , n ( x ) = cos λ n x , Y 2 , n ( x ) = x sin λ n x , λ n = 2 π n , n N .
It is a basis in L 2 ( 0 , 1 ) , orthonormal to { X 0 , n , X 1 , n , X 2 , n } n = 1 , that is X i , n , Y j , m = δ i j δ n m .
The system of functions { X 0 , n , X 1 , n , X 2 , n } n = 1 defined in (32) is a Riesz basis in L 2 ( 0 , 1 ) [41,42,43]. This means that for any f L 2 ( 0 , 1 ) the following estimates are satisfied
C f L 2 ( 0 , 1 ) 2 f 1 , 0 2 + n = 1 f 1 , n 2 + f 2 , n 2 C + f L 2 ( 0 , 1 ) 2
for some constants C ± > 0 .
Proposition 2.
Assume f H 2 ( 0 , 1 ) and f ( 1 ) = 0 , f ( 0 ) = f ( 1 ) . Then there exist constants C 1 ± , such that the estimates hold true
C 1 f H 2 ( 0 , 1 ) 2 f 1 , 0 2 + n = 1 ( 1 + λ n 4 ) f 1 , n 2 + f 2 , n 2 C 1 + f H 2 ( 0 , 1 ) 2 .
Proof. 
Let f satisfies the assumptions of the proposition. Taking into account the elementary inequalities a 2 + b 2 ( a + b ) 2 2 ( a 2 + b 2 ) for a , b 0 , it is sufficient to prove the estimates
C 2 f L 2 ( 0 , 1 ) 2 n = 1 λ n 4 f 1 , n 2 + f 2 , n 2 C 2 + f L 2 ( 0 , 1 ) 2 .
Applying twice integration by parts and taking into account that f ( 1 ) = 0 , f ( 0 ) = f ( 1 ) , we find by the use of (34) the following expressions for the coefficients of f ( x )
f 1 , 0 = f , Y 1 , 0 = 0 , f 1 , n = f , Y 1 , n = λ n 2 f 1 , n , f 2 , n = f , Y 2 , n = 2 λ n f 1 , n λ n 2 f 2 , n , n N .
In this way, we obtain the expansion
f ( x ) = n = 1 λ n 2 f 1 , n X 1 , n ( x ) + ( 2 λ n f 1 , n λ n 2 f 2 , n ) X 2 , n ( x ) ,
which gives by applying the Riesz basis property (36) the estimates
C f L 2 ( 0 , 1 ) 2 n = 1 λ n 4 f 1 , n 2 + ( 2 λ n f 1 , n λ n 2 f 2 , n ) 2 C + f L 2 ( 0 , 1 ) 2 .
The lower bound in (38) follows easily from the lower bound in (40) using the inequalities ( 2 λ n f 1 , n λ n 2 f 2 , n ) 2 8 λ n 2 f 1 , n 2 + 2 λ n 4 f 2 , n 2 and λ n 2 < λ n 4 , n N .
To deduce the upper bound in (38) from the upper bound in (40), we use the following identities obtained from (34) by integration by parts
f 1 , n = f , cos λ n x = 1 λ n 2 f , cos λ n x ;
f 2 , n = x f ( x ) , sin λ n x = 1 λ n 2 x f ( x ) , sin λ n x 2 λ n 3 f ( x ) , cos λ n x .
By applying the Bessel inequality for the trigonometric series identity (41) yields
n = 1 λ n 4 f 1 , n 2 = n = 1 f , cos λ n x 2 f L 2 ( 0 , 1 ) 2
and, in a similar way, (42) implies
n = 1 λ n 4 f 2 , n 2 C 1 n = 1 x f ( x ) , sin λ n x 2 + C 2 n = 1 f ( x ) , cos λ n x 2 C 1 x f ( x ) L 2 ( 0 , 1 ) 2 + C 2 f ( x ) L 2 ( 0 , 1 ) 2 C 3 f L 2 ( 0 , 1 ) 2 .
In this way, combining (43) and (44), the upper bound in (38) is also established. □

6. Formal Spectral Expansions for the Solution

In this section, we find formal spectral expansions for the unknown functions h ( x ) and u ( x , t ) in the inverse source problem (6)–(8). Suppose
h ( x ) = h 1 , 0 X 1 , 0 ( x ) + n = 1 h 1 , n X 1 , n ( x ) + h 2 , n X 2 , n ( x ) ,
u ( x , t ) = u 1 , 0 ( t ) X 1 , 0 ( x ) + n = 1 u 1 , n ( t ) X 1 , n ( x ) + u 2 , n ( t ) X 2 , n ( x ) .
To find the unknown coefficients in (45) and (46), we insert these expansions in Equation (6) and, by taking into account (39) and the initial condition u ( x , 0 ) = 0 , we obtain by the uniqueness property of the spectral expansion the following system of equations
D t ( k ) u 1 , n ( t ) + λ n 2 u 1 , n ( t ) = h 1 , n q ( t ) , u 1 , n ( 0 ) = 0 , n N 0 ,
D t ( k ) u 2 , n ( t ) + λ n 2 u 2 , n ( t ) = 2 λ n u 1 , n ( t ) + h 2 , n q ( t ) , u 2 , n ( 0 ) = 0 , n N .
We solve first relaxation equation (47) and after that (48). In this way, by applying formula (25) we obtain
u 1 , n ( t ) = h 1 , n A n ( t ) , n N 0 ,
u 2 , n ( t ) = h 2 , n A n ( t ) + 2 λ n h 1 , n B n ( t ) , n N .
Here the functions A n ( t ) and B n ( t ) are defined as follows
A n ( t ) = ( q G n ) ( t ) , B n ( t ) = ( q G n G n ) ( t ) ; G n ( t ) = G ( t ; λ n 2 ) ,
where G ( t ; λ ) is defined in (26). In this way, we obtain
u ( x , t ) = h 1 , 0 A 0 ( t ) X 1 , 0 ( x ) + n = 1 h 1 , n A n ( t ) X 1 , n ( x ) + ( h 2 , n A n ( t ) + 2 λ n h 1 , n B n ( t ) ) X 2 , n ( x ) .
Therefore, the function g ( x ) = u ( x , T ) admits the spectral expansion
g ( x ) = h 1 , 0 A 0 ( T ) X 1 , 0 ( x ) + n = 1 h 1 , n A n ( T ) X 1 , n ( x ) + ( h 2 , n A n ( T ) + 2 λ n h 1 , n B n ( T ) ) X 2 , n ( x ) .
By the uniqueness property of the spectral expansion, from (53) we deduce
g 1 , n = h 1 , n A n ( T ) , n N 0 ; g 2 , n = h 2 , n A n ( T ) + 2 λ n h 1 , n B n ( T ) , n N .
Since A n ( T ) > 0 for all n N 0 (see Proposition 3), it follows from (54)
h 1 , n = 1 A n ( T ) g 1 , n , n N 0 ; h 2 , n = 1 A n ( T ) g 2 , n 2 λ n g 1 , n B n ( T ) A n ( T ) , n N .
Relations (54) and (55) define a mapping between the overdetermination function g ( x ) and the source factor h ( x ) and our main goal is to determine the properties of this mapping.
Plugging the coefficients h j , n in (52), we derive the coefficients u j , n ( t ) in the spectral expansion (46) of the solution u ( x , t )
u 1 , n ( t ) = g 1 , n A n ( t ) A n ( T ) , n N 0 ;
u 2 , n ( t ) = g 2 , n A n ( t ) A n ( T ) + 2 λ n g 1 , n B n ( t ) A n ( T ) A n ( t ) B n ( T ) A n 2 ( T ) , n N .
In this way, inserting (55) in (45), and (56), (57) in (46), we obtain the formal expansions for h ( x ) and u ( x , t ) . The functions A n ( t ) and B n ( t ) in these expansions depend on the specific memory kernel k ( t ) via G n ( t ) and the time-dependent source factor q ( t ) .

7. Uniqueness of Solution and Stability Estimates in Sobolev Spaces

In this section, we prove that under some assumptions on the overdetermination function g ( x ) and the time-dependent source factor q ( t ) , the formal expansions (45) with coefficients (55) and (46) with coefficients (56) and (57) define a unique solution ( h , u ) to the considered inverse source problem.
For the given time-dependent source factor q ( t ) we assume
q C [ 0 , T ] ; q ( t ) 0 , t [ 0 , T ] ; q ( T ) > 0 .
Let us note that since q ( t ) is a continuous function for t ( 0 , T ] , the assumption q ( T ) > 0 implies that there exist q 0 > 0 and T 1 ( 0 , T ) , such that
q ( t ) q 0 > 0 f o r a l l t [ T 1 , T ] .
First, we prove some estimates for the functions A n ( t ) and B n ( t ) .
Proposition 3.
Let T > 0 be arbitrarily fixed. Assume that conditions (3) and (58) for the functions k ( t ) and q ( t ) , respectively, are satisfied. Then the functions A 0 ( t ) , A n ( t ) , B n ( t ) , n N , are continuous and non-negative on [ 0 , T ] , vanish at t = 0 , and A n ( T ) > 0 for all n N 0 . Moreover, the following estimates for t [ 0 , T ] and n N are satisfied:
A n ( T ) C 1 / λ n 2 ,
A n ( t ) C 2 / λ n 2 ,
B n ( t ) C 2 / λ n 4 ,
where the constants C j > 0 , j = 1 , 2 , are independent of n.
Proof. 
The functions A n ( t ) and B n ( t ) are continuous and non-negative as convolutions of functions with these properties. Fix T 1 ( 0 , T ) , such that (59) is satisfied. For n 0 we deduce
A n ( T ) = 0 T q ( t ) G n ( T t ) d t T 1 T q ( t ) G n ( T t ) d t q 0 0 T T 1 G n ( t ) d t ,
due to the non-negativity of the functions under the integral sign and inequality (59). Then (63) together with (31) implies A n ( T ) > 0 , n N 0 .
Estimates (63) and (30) in Proposition 1 imply for n N
λ n 2 A n ( T ) q 0 1 S T T 1 ; λ 1 2 > 0 ,
where λ 1 2 = 4 π 2 is the smallest positive eigenvalue and the function S ( t ; λ ) is defined in (26). In this way, (60) is established with C 1 = q 0 ( 1 S ( T T 1 ; λ 1 2 ) ) > 0 .
Next, we prove estimates (61) and (62). The representation in (51) and the estimate from above in (30) yield
A n ( t ) = q G n ( t ) q C [ 0 , T ] 1 G n ( t ) C 2 λ n 2
which implies (61) with C 2 = q C [ 0 , T ] > 0 . By applying (61) and (30) we derive (62):
B n ( t ) = q G n G n ( t ) = A n G n ( t ) C 2 λ n 2 1 G n ( t ) C 2 λ n 4 .
Since ( 1 G n ) ( 0 ) = 0 , the above two estimates also imply that A n ( 0 ) = 0 and B n ( 0 ) = 0 . □
We are ready to formulate and prove the main result of this work.
Theorem 1.
Let T > 0 be arbitrarily fixed. For any given g H 2 ( 0 , 1 ) , such that g ( 1 ) = 0 , g ( 0 ) = g ( 1 ) , there exists a unique solution ( h , u ) to problem (6)–(8), satisfying h L 2 ( 0 , 1 ) and u C ( [ 0 , T ] ; H 2 ( 0 , 1 ) ) . The functions h ( x ) and u ( x , t ) are defined by the spectral expansions (45) and (46) with coefficients given in (55), (56) and (57). Moreover, there exist constants c ̲ > 0 and c ¯ > 0 , such that
c ̲ g H 2 ( 0 , 1 ) h L 2 ( 0 , 1 ) c ¯ g H 2 ( 0 , 1 ) .
Proof. 
Uniqueness of the solution follows from the uniqueness property of the spectral expansions and the fact that A n ( T ) 0 , for all n 0 (see Proposition 3). Indeed, if g 0 , then g j , n = 0 , and therefore, (55)–(57) imply h j , n = 0 , u j , n = 0 , i.e., all coefficients in the expansions (45) and (46) vanish and, therefore, h 0 and u 0 .
The initial condition u ( x , 0 ) = 0 is satisfied since A n ( 0 ) = 0 and B n ( 0 ) = 0 . The two boundary conditions (7) hold by construction, since they are satisfied by the basis functions X 1 , n ( x ) and X 2 , n ( x ) defined in (32).
Applying estimates (60), (62), and taking into account the inequalities ( a ± b ) 2 2 ( a 2 + b 2 ) and λ n > 1 , n N , identities (55) imply
h 1 , n 2 C 1 λ n 4 g 1 , n 2 , h 2 , n 2 C 2 λ n 4 g 1 , n 2 + g 2 , n 2 , n N .
Estimates (65) together the lower bound in (36) with f = h and the upper bound in (37) with f = g imply the upper bound in (64). By applying similar argument, we obtain from (54) by the use of Proposition 3 the estimates
λ n 4 g 1 , n 2 C 3 h 1 , n 2 , λ n 4 g 2 , n 2 C 4 h 1 , n 2 + h 2 , n 2 , n N .
The last estimates together with the lower bound in (37) with f = g and the upper bound in (36) with f = h yield the lower bound in (64).
Furthermore, Proposition 3 implies for t [ 0 , T ]
A n ( t ) A n ( T ) 2 C 5 , B n ( t ) A n ( T ) A n ( t ) B n ( T ) A n 2 ( T ) 2 2 B n ( t ) A n ( T ) 2 + 2 C 5 B n ( T ) A n ( T ) 2 C 6 λ n 2 .
Representations (56) and (57) for the coefficients of the solution u ( x , t ) yield by the use of estimates (67)
u 1 , n 2 ( t ) C 7 g 1 , n 2 , u 2 , n 2 ( t ) C 8 g 1 , n 2 + g 2 , n 2 , t [ 0 , T ] , n N .
Since g H 2 ( 0 , 1 ) , it follows from (68) that
u C ( [ 0 , T ] ; H 2 ( 0 , 1 ) ) C 9 g H 2 ( 0 , 1 ) < .
Therefore u C ( [ 0 , T ] ; H 2 ( 0 , 1 ) ) . The proof of the theorem is complete. □
It is worth noting that if q ( t ) C 1 [ 0 , T ] then the inverse source problem under consideration can be reduced to a Fredholm equation of second kind, by applying the argument proposed in [17]. In general, this fact is well known for parabolic equations (e.g., [44]). In this case, it is sufficient to prove the uniqueness for the inverse problem, and then the stability follows by the Fredholm alternative. However, in the present work, continuous differentiability of q ( t ) is not assumed, see (58), and a specialized argument is needed.
Adopting the method of [17], the more general case q = q ( x , t ) can be also treated by the use of the generalized eigenfunction expansion, provided the function q ( x , t ) obeys sufficient regularity. Although no explicit series representation of the source term h ( x ) can be derived in this case, the estimates (64) can be established by applying the Fredholm alternative in L 2 ( 0 , 1 ) .
We close this section with a conditional stability result. It concerns the case when the overdetermination function g is only square integrable, g L 2 ( 0 , 1 ) , which is relevant in the study of stability with respect to the noise level applied to the input, see (9). In this case, an a priori bound assumption for the unknown function h ( x ) is needed, see e.g., [21]. Let us assume the following a priori bound: h H 2 ( 0 , 1 ) E . Then
h L 2 ( 0 , 1 ) C E g L 2 ( 0 , 1 ) 1 / 2
where the constant C E depends on E and the constants appearing in the inequalities (36), (37), (65) and (66). Estimate (69) can be derived adapting the technique in [35] to the case of generalized eigenfunction expansions by the use of the estimates (36), (37), (65) and (66) and the Cauchy-Schwarz inequality. Details of the proof are omitted here.

8. Concluding Remarks

This work is concerned with an inverse source problem for the one-dimensional generalized subdiffusion equation on a bounded interval. The differences to the classical case are as follows. First, temporal derivative in the diffusion equation is replaced by a general convolutional derivative in time. This allows the modeling of memory effects of the described physical process. Second, traditional Dirichlet boundary condition is applied at the one side of the interval only and the other boundary condition is nonlocal.
The main result of the paper is Theorem 1 establishing a two-sided estimate of the L 2 norm of the spatial part h ( x ) of the forcing term via the Sobolev H 2 norm of the final overdetermination function g ( x ) . The derived estimates extend earlier obtained results on the classical diffusion problem and the anomalous diffusion problem with the Caputo time-fractional derivative subject to Dirichlet boundary conditions. The proofs are based on generalized eigenfunction expansions with respect to a biorthogonal pair of bases, together with estimates of the special functions in time appearing in the expansions, which inherit relevant properties of Mittag–Leffler functions.
It is feasible to extend the results to multidimensional cases, such as the two-dimensional nonlocal problem in [34].
Another approach to further study of such types of inverse problems is to employ the general theory for parabolic equations [44]. Under appropriate assumptions on the given functions the problem can be reduced to a Fredholm equation of the second kind and then the Fredholm alternative can be applied.

Funding

This work is supported by the National Scientific Program “Information and Communication Technologies for a Single Digital Market in Science, Education and Security (ICTinSES)”, contract No DO1–205/23.11.2018, financed by the Ministry of Education and Science in Bulgaria.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author is grateful to the anonymous reviewers for the constructive comments. The work is performed in the frames of the bilateral research project “Operators, differential equations and special functions of Fractional Calculus---numerics and applications” between Bulgarian Academy of Sciences and Serbian Academy of Sciences and Arts.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Metzler, R.; Jeon, J.H.; Cherstvy, A.G.; Barkai, E. Anomalous diffusion models and their properties: Non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking. Phys. Chem. Chem. Phys. 2014, 16, 24128. [Google Scholar] [CrossRef] [Green Version]
  2. Kochubei, A.; Luchko, Y. (Eds.) Fractional Differential Equations; De Gruyter: Berlin, Germany; Boston, MA, USA, 2019. [Google Scholar]
  3. Guo, B.; Pu, X.; Huang, F. Fractional Partial Differential Equations and Their Numerical Solutions; World Scientific: Singapore, 2015. [Google Scholar]
  4. Sandev, T.; Sokolov, I.M.; Metzler, R.; Chechkin, A. Beyond monofractional kinetics. Chaos Solitons Fractals 2017, 102, 210–217. [Google Scholar] [CrossRef]
  5. Sandev, T.; Metzler, R.; Chechkin, A. From continuous time random walks to the generalized diffusion equation. Fract. Calc. Appl. Anal. 2018, 21, 10–28. [Google Scholar] [CrossRef] [Green Version]
  6. Sandev, T.; Tomovski, Z.; Dubbeldam, J.L.A.; Chechkin, A. Generalized diffusion-wave equation with memory kernel. J. Phys. Math. Theor. 2019, 52, 015201. [Google Scholar] [CrossRef] [Green Version]
  7. Kochubei, A. General fractional calculus, evolution equations, and renewal processes. Integr. Equ. Oper. Theory 2011, 71, 583–600. [Google Scholar] [CrossRef] [Green Version]
  8. Luchko, Y.; Yamamoto, M. The general fractional derivative and related fractional differential equations. Mathematics 2020, 8, 2115. [Google Scholar] [CrossRef]
  9. Luchko, Y. Operational calculus for the general fractional derivative and its applications. Fract. Calc. Appl. Anal. 2021, 24, 338–375. [Google Scholar] [CrossRef]
  10. Schumer, R.; Benson, D.; Meerschaert, M.; Baeumer, B. Fractal mobile/immobile solute transport. Water Resour. Res. 2003, 39, SBH131–SBH1312. [Google Scholar] [CrossRef]
  11. Atanacković, T.M.; Pilipović, S.; Stanković, B.; Zorica, D. Fractional Calculus with Applications in Mechanics: Vibrations and Diffusion Processes; John Wiley & Sons: London, UK, 2014. [Google Scholar]
  12. Bazhlekova, E.; Bazhlekov, I. Transition from diffusion to wave propagation in fractional Jeffreys-type heat conduction equation. Fractal Fract. 2020, 4, 32. [Google Scholar] [CrossRef]
  13. Schilling, R.; Song, R.; Vondraček, Z. Bernstein Functions: Theory and Applications; De Gruyter: Berlin, Germany, 2010. [Google Scholar]
  14. Clement, P.; Prüss, J. Completely positive measures and Feller semigroups. Math. Ann. 1990, 287, 73–105. [Google Scholar] [CrossRef]
  15. Prüss, J. Evolutionary Integral Equations and Applications; Birkhäuser: Basel, Switzerland, 1993. [Google Scholar]
  16. Vergara, V.; Zacher, R. Optimal Decay Estimates for Time-Fractional and Other NonLocal Subdiffusion Equations via Energy Methods. SIAM J. Math. Anal. 2015, 47, 210–239. [Google Scholar] [CrossRef]
  17. Sakamoto, K.; Yamamoto, M. Inverse source problem with a final overdetermination for a fractional diffusion equation. Math. Control Relat. Fields 2011, 1, 509–518. [Google Scholar] [CrossRef]
  18. Tuan, N.H.; Debbouche, A.; Ngoc, T.B. Existence and regularity of final value problems for time fractional wave equations. Comput. Math. Appl. 2019, 78, 1396–1414. [Google Scholar] [CrossRef]
  19. Wei, T.; Wang, J. A modified quasi-boundary value method for an inverse source problem of the time-fractional diffusion equation. Appl. Num. Math. 2014, 78, 95–111. [Google Scholar] [CrossRef]
  20. Zhang, Y.; Xu, X. Inverse source problem for a fractional diffusion equation. Inverse Probl. 2011, 27, 035010. [Google Scholar] [CrossRef]
  21. Zhang, Z.; Wei, T. Identifying an unknown source in time-fractional diffusion equation by a truncation method. Appl. Math. Comp. 2013, 219, 5972–5983. [Google Scholar] [CrossRef]
  22. Jin, B.; Rundell, W. A tutorial on inverse problems for anomalous diffusion processes. Inverse Probl. 2015, 31, 035003. [Google Scholar] [CrossRef] [Green Version]
  23. Tuan, N.; Kirane, M.; Hoan, L.; Long, L. Identification and regularization for unknown source for a time-fractional diffusion equation. Comput. Math. Appl. 2017, 73, 931–950. [Google Scholar] [CrossRef]
  24. Slodička, M.; Šišková, K.; Bockstal, K.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]
  25. Slodička, 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, 1702–1711. [Google Scholar] [CrossRef]
  26. Kinash, N.; Janno, J. An inverse problem for a generalized fractional derivative with an application in reconstruction of time- and space-dependent sources in fractional diffusion and wave equations. Mathematics 2019, 7, 1138. [Google Scholar] [CrossRef] [Green Version]
  27. Kinash, N.; Janno, J. Inverse problems for a generalized subdiffusion equation with final overdetermination. Math. Model. Anal. 2019, 24, 236–262. [Google Scholar] [CrossRef]
  28. Kinash, N.; Janno, J. Inverse problems for a perturbed time fractional diffusion equation with final overdetermination. Math. Meth. Appl. Sci. 2018, 41, 1925–1943. [Google Scholar] [CrossRef]
  29. Kaliev, I.; Sabitova, M. Problems of determining the temperature and density of heat sources from the initial and final temperatures. J. Appl. Indust. Math. 2010, 4, 332–339. [Google Scholar] [CrossRef]
  30. Kirane, M.; Malik, S. Determination of an unknown source term and the temperature distribution for the linear heat equation involving fractional derivative in time. Appl. Math. Comp. 2011, 218, 163–170. [Google Scholar] [CrossRef] [Green Version]
  31. Furati, K.; Iyiola, O.; Kirane, M. An inverse problem for a generalized fractional diffusion. Appl. Math. Comp. 2014, 249, 24–31. [Google Scholar] [CrossRef]
  32. Furati, K.; Iyiola, O.; Mustapha, K. An inverse source problem for a two-parameter anomalous diffusion with local time datum. Comput. Math. Appl. 2017, 73, 1008–1015. [Google Scholar] [CrossRef] [Green Version]
  33. Ali, M.; Aziz, S.; Malik, S.A. Inverse problem for a multi-term fractional differential equation. Fract. Calc. Appl. Anal. 2020, 23, 799–821. [Google Scholar] [CrossRef]
  34. Kirane, M.; Malik, S.; Al-Gwaiz, M. An inverse source problem for a two dimensional time fractional diffusion equation with nonlocal boundary conditions. Math. Meth. Appl. Sci. 2013, 36, 1056–1069. [Google Scholar] [CrossRef]
  35. Bazhlekova, E. Estimates for a general fractional relaxation equation and application to an inverse source problem. Math. Met. Appl. Sci. 2018, 41, 9018–9026. [Google Scholar] [CrossRef] [Green Version]
  36. Bazhlekova, E.; Bazhlekov, I. Identification of a space-dependent source term in a nonlocal problem for the general time-fractional diffusion equation. J. Comput. Appl. Math. 2021, 386, 113213. [Google Scholar] [CrossRef]
  37. Gripenberg, G.; Londen, S.O.; Staffans, O. Volterra Integral and Functional Equations; Cambridge University Press: Cambridge, UK, 1990. [Google Scholar]
  38. Bazhlekov, I.; Bazhlekova, E. Fractional derivative modeling of bioreaction-diffusion processes. AIP Conf. Proc. 2021, 2333, 060006. [Google Scholar] [CrossRef]
  39. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  40. Bazhlekova, E. Completely monotone multinomial Mittag-Leffler type functions and diffusion equations with multiple time-derivatives. Fract. Calc. Appl. Anal. 2021, 24, 88–111. [Google Scholar] [CrossRef]
  41. Ionkin, N. Solution of a boundary problem for a heat conduction equation with nonclassical boundary condition. Differ. Equ. 1977, 13, 294–304. [Google Scholar]
  42. Ionkin, N.; Moiseev, E. A two-point boundary value problem for a heat conduction equation. Differ. Equ. 1979, 15, 1284–1295. [Google Scholar]
  43. Moiseev, E. The solution of a nonlocal boundary value problem by the spectral method. Differ. Equ. 1999, 35, 1105–1112. [Google Scholar]
  44. Prilepko, A.I.; Orlovsky, D.G.; Vasin, I.A. Methods for Solving Inverse Problems in Mathematical Physics; Marcel-Dekker: New York, NY, USA, 2000. [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

Bazhlekova, E. An Inverse Source Problem for the Generalized Subdiffusion Equation with Nonclassical Boundary Conditions. Fractal Fract. 2021, 5, 63. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5030063

AMA Style

Bazhlekova E. An Inverse Source Problem for the Generalized Subdiffusion Equation with Nonclassical Boundary Conditions. Fractal and Fractional. 2021; 5(3):63. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5030063

Chicago/Turabian Style

Bazhlekova, Emilia. 2021. "An Inverse Source Problem for the Generalized Subdiffusion Equation with Nonclassical Boundary Conditions" Fractal and Fractional 5, no. 3: 63. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5030063

Article Metrics

Back to TopTop