Preprint
Article

Fractional Integro-Differential Diffusion Equations with Nonlocal Initial Conditions via the Resolvent Family

Altmetrics

Downloads

137

Views

31

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

05 September 2023

Posted:

07 September 2023

You are already at the latest version

Alerts
Abstract
The fractional diffusion equations with integrals play a significant role in describing anomalous diffusion phenomena. We first construct an appropriate resolvent family, and use the convolution theorem of Laplace transform, the probability density function, Cauchy integral formula and Fubini theorem to study its related equicontinuity, strongly continuity, compactness, etc. Then, a reasonable mild solution is defined. Finally, by combining some fixed point theorems, we obtain some conclusions on the existence and uniqueness of mild solutions.
Keywords: 
Subject: Computer Science and Mathematics  -   Mathematics

1. Introduction

In this paper, we consider the existence and uniqueness of mild solutions to the following nonlinear fractional integro-differential diffusion equations with nonlocal initial conditions:
t α u ( x , t ) = A u ( x , t ) + I t β f ( x , t , u ( x , t ) ) , ( x , t ) Ω × ( 0 , T ) , u ( x , t ) = 0 , ( x , t ) Ω × ( 0 , T ) , u ( x , 0 ) + g ( u ) = u 0 ( x ) , ( x , t ) Ω × ( 0 , T ] ,
where α ( 0 , 1 ) , β ( 0 , ) , t α and I t β are the α order partial Caputo derivative and β order partial Riemann-Liouville integral with respect to t. u 0 ( x ) H 1 ( Ω ) , which is given in Section 2, Ω R n with smooth boundary Ω , T > 0 , the nonlinear term f and the nonlocal term g are given functions. The coefficient linear operator A is defined by:
A u ( x , t ) = i = 1 n j = 1 n x i a i j ( x ) x j u ( x , t ) p ( x ) u ( x , t ) ,
where a i j ( x ) ( 1 i , j n ) are some real valued functions in Ω satisfying conditions:
a i j ( x ) L ( Ω ) , 1 i , j n ,
i , j = 1 n a i j ( x ) ϑ i ϑ j ς | ϑ | 2 , ϑ R n , a . e . x Ω ,
with some constants ς > 0 , p ( x ) Ω is also a real valued function satisfying
p ( x ) L ( Ω ) , p ( x ) p 0 > 0 , a . e . x Ω ,
where L ( Ω ) denotes the space composed of measurable functions defined on Ω and almost bounded everywhere.
Fractional differential equations originated in 1695, see [1,2,3]. As we have known that they can provide such excellent descriptive models to resolve various problems in reality, to the extent that they are applied in various fields, such as control engineering [4,5], viscoelastic materials [6], fluid mechanics [7], electrochemistry [8], the analysis of epidemic [9] and complex networks [10], and statistical mechanics [11], etc. The relevant problems for the diffusion equations have been studied by many scholars, see [12,13,14] and recent relevant references [15,16,17].
The existence results to nonlocal initial problems in Banach space was initiated by Byszewski and Lakshmikantham [18]. The motivation for these studies is that the nonlocal condition is better to describe the diffusion phenomenon than using the usual local Cauchy problem u ( 0 ) = u 0 . For example, g ( u ) can be given by
g ( u ) = i = 1 l c i u ( τ i ) ,
where c i ( i = 1 , 2 , , l ) are given constants and 0 < τ 1 < τ 2 < < τ l < T . In addition, for some applications of nonlocal conditions, please refer to [19,20,21,22,23].
Fractional equations containing only differential terms have been studied widely. For instance, in [24], Mu et al. considered the initial-boundary value problem of fractional diffusion equations in Caputo sense:
t α u ( x , t ) = A u ( x , t ) + f ( x , t ) in Ω × ( 0 , b ) , u = 0 on Ω × ( 0 , b ) , u ( x , 0 ) = u 0 ( x ) in Ω ,
where f is weighted H o ¨ lder continuous. The existence, uniqueness and regularity of solutions of (3) be established in C ( [ 0 , b ] , H 1 ( Ω ) ) C ( ( 0 , b ] , D ( A ) ) C α ( ( 0 , b ] , H 1 ( Ω ) ) with some assumptions. Nevertheless, a lot of practical phenomena can be depicted via appropriate models, which include differential and integral terms. Hence the appearance of integro-differential equations shows that their excellent applicability in some physical or engineering areas. In [25], Amin et al. obtained the solutions of a integro-differential equation with initial condition
t ξ w ( t ) = r ( t ) + b ( t ) w ( t ) + 0 t W ( t , τ ) w ( τ ) d τ , t [ 0 , 1 ] , w ( 0 ) = w 0 ,
where ξ ( 0 , 1 ] , W is the kernel of integral, r ( t ) and b ( t ) are known. This integral term has certain limitations when describing nonlocal diffusion phenomena, and the applicability of the initial condition is also relatively weak.
According to our knowledge, the mild solutions of fractional differential equations are usually represented by the Mittag-Leffler function or probability density function, see [26,27,28,29,30,31,32]. In [32], Zhou et al. obtained a mild solution, where
λ q 1 ( λ q I A ) 1 = 0 e λ t 0 ψ q ( θ ) T t q θ q d θ d t , q ( 0 , 1 ) ,
ψ q ( θ ) is the probability density function defined on ( 0 , ) and { T ( t ) } t 0 is a C 0 -semigroup generated by A. On the other hand, if f ( x , t , u ( x , t ) ) = g ( u ( x , t ) ) = 0 , then a formal solution of (1) is
u ( x , t ) = i = 1 ( u 0 , X i ) E α ( λ i t α ) X i ( x ) .
Where X i is the eigenfunction related to the eigenvalue λ i of the corresponding eigenvalue problem, that is, A ( X i ) = λ i X i , i = 1 , 2 , , see further details in [29]. Obviously, using these technique directly to solve the problem (1) is quite difficult.
Based on the above, in this paper we apply the ( α , β ) -resolvent family to discuss the mild solutions. Resolvent families are powerful for studying solutions to fractional diffusion equations. Chen et al. [33] established the existence and controllability estimation of mild solutions for a class of evolution equations with nonlocal conditions through a resolvent family. Ponce [23] obtained properties on the behavior of mild solutions for fraction Cauchy problems by a resolvent family. Later, Chang et al. [19] proved that if the source function of a diffusion equation has vector value periodicity or almost periodicity or almost automorphism, then the diffusion equation has a mild solution through a resolvent family. Although there is also an integral of f in [19], some of the proof techniques in which are not applicable to this article due to the derivative order of u being different.
This paper is organized as follows. In Section 2, by selecting the appropriate space, we transform (1) into an abstract Cauchy problem and provide some necessary definitions and preliminary results that will be used in the sequel. Afterwards, we define the mild solutions for (9) by Laplace transform. In Section 3, the existence and uniqueness of the mild solutions be established by several fixed point theorems under some assumptions.

2. Preliminaries

In this section, we provide some definitions and lemmas about fractional calculus and the ( α , β ) -resolvent family that will be used in this paper.
Definition 1 
([2]). The Riemann-Liouville fractional integral of order α ( 0 , 1 ) with respect to t for an integrable function f : [ 0 , ) R is defined as
I t α f ( t ) = 1 Γ ( α ) 0 t ( t s ) α 1 f ( s ) d s , t > 0 ,
where Γ ( · ) is the Gamma function.
Definition 2 
([2]). The Riemann-Liouville fractional derivative of order α ( 0 , 1 ) with respect to t for an absolutely continuous function f : [ 0 , ) R is defined as
R D t α f ( t ) = 1 Γ ( 1 α ) d d t 0 t ( t s ) α f ( s ) d s , t > 0 .
Definition 3 
([2]). The Caputo fractional derivative of order α ( 0 , 1 ) with respect to t for an absolutely continuous function f : [ 0 , ) R can be written as
D t α f ( t ) = R D t α ( f ( t ) f ( 0 ) ) , t > 0 .
If 0 < α < 1 and f : [ 0 , ) R is continuously differentiable, then
D t α f ( t ) = 1 Γ ( 1 α ) 0 t f ( s ) ( t s ) α d s = I t 1 α f ( t ) .
If f is an abstract function with values in a Banach space, then the integrals and derivatives appearing in (5) and (7) are understood in Bochner’s sense.
Definition 4 
([34]). Let P be a metric space, and let S P be a bounded set. The Kuratowski measure of noncompactness is defined by:
ν ( S ) = inf { δ > 0 | S = i = 1 m S i , diam S i δ } .
Lemma 5 
([35]). Let E be a Banach space, and S, S 1 , S 2 be some subsets of E. Then we have the following properties :
(i) ν ( S ) = 0 S is relatively compact;
(ii) ν ( S ) = ν ( S ¯ ) ;
(iii) S 1 S 2 ν ( S 1 ) ν ( S 2 ) ;
(iv) ν ( S 1 + S 2 ) ν ( S 1 ) + ν ( S 2 ) , where S 1 + S 2 = { x | x = y + z , y S 1 , z S 2 } ;
(v) ν ( c S ) = | c | ν ( S ) , c R , where c S = { x | x = c z , z S } ;
(vi) 0 ν ( S ) < + .
Let X = H 1 ( Ω ) with the norm · and J = [ 0 , T ] . Here H 1 ( Ω ) represents the dual space of H 0 1 ( Ω ) , which is the closure of C c ( Ω ) in H 1 ( Ω ) , C c ( Ω ) consists of the functions in C ( Ω ) with compact support, H 1 ( Ω ) is the Sobolev space (see [36]). C ( J , X ) denotes the Banach space of all continuous functions from J into X with the norm
· = sup t J { · ( t ) } .
Similarly, C ( J , R + ) denotes all continuous functions from J into R + with the norm
· = sup t J { · ( t ) } .
We define A : D ( A ) X X with D ( A ) = H 0 1 ( Ω ) and ( A u ) ( t ) x = A u ( x , t ) . Then according to [37], A generates an analytic semigroup { T ( t ) } t 0 on X. Without loss of generality, we can assume that 0 ρ ( A ) , { T ( t ) } t 0 is uniformly bounded, and there exists a constant N > 0 such that
T ( t ) N , for t > 0 .
Denote ρ ( A ) and R ( λ , A ) = ( λ I A ) 1 are the resolvent set and resolvent operator of A, respectively, where I is the identity operator. By [38], there exist δ ( 0 , π 2 ) and C > 0 such that
Σ π 2 + δ : = λ | | arg λ | < π 2 + δ ρ ( A )
and
| | R ( λ , A ) | | C | λ | , λ Σ π 2 + δ { 0 } .
Set u ( t ) ( x ) = u ( x , t ) , f ( t , u ( t ) ) ( x ) = f ( x , t , u ( x , t ) ) , where x X , then (1) can be formulated as an abstract problem with nonlocal initial conditions:
D t α u ( t ) = A u ( t ) + I t β f ( t , u ( t ) ) , t ( 0 , T ] , u ( 0 ) + g ( u ) = u 0 ,
where D t α and I t β denote the α order Caputo derivative and β order Riemann-Liouville integral, respectively. g : C ( J , X ) X , f : J × X X and u 0 X .
Definition 6. 
If A generates a uniformly bounded and analytic semigroup, which satisfies (10) and (11), and for operator-valued function S α , β ( t ) : X X we have
λ β ( λ α I A ) 1 x = 0 e λ t S α , β ( t ) x d t , t 0 , x X ,
then { S α , β ( t ) } t 0 is called the ( α , β ) resolvent family generated by A.
It can be seen from the reference [24,31] that the following results can be obtained.
Remark 1. 
Let Φ ( t ) and Λ ( t ) satisfy
λ α 1 ( λ α I A ) 1 x = 0 e λ t Φ ( t ) x d t , t 0 , x X ,
and
( λ α I A ) 1 x = 0 e λ t Λ ( t ) x d t , t 0 , x X .
Then
Λ ( t ) = t α 1 Ψ ( t ) , Ψ ( t ) = α 0 θ ζ α ( θ ) T ( t α θ ) d θ ,
S α , 1 α ( t ) = Φ ( t ) = 0 ζ α ( θ ) T ( t α θ ) d θ ,
where ζ α ( θ ) is a probability density function defined on ( 0 , ) which satisfies
ζ α ( θ ) 0 , θ ( 0 , ) ,
0 θ v ζ α ( θ ) d θ = Γ ( 1 + v ) Γ ( 1 + α v ) for v ( 1 , ) ,
0 e z θ ζ α ( θ ) d θ = E α ( z ) , α 0 θ e z θ ζ α ( θ ) d θ = E α , α ( z ) , for z C ,
E α , E α , α are the Mittag-Leffler functions [2].
If u : J X is a solution of
D t α u ( t ) = A u ( t ) + f ( t , u ( t ) ) , t ( 0 , T ] , u ( 0 ) = u 0 , u 0 X ,
then
u ( t ) = S α , 1 α ( t ) u 0 + 0 t ( t s ) α 1 Ψ ( t s ) f ( s , u ( s ) ) d s .
Remark 2. 
Assume that { S α , β ( t ) } t 0 is the ( α , β ) resolvent family generated by A. Because the inverse Laplace transform of λ β is t β 1 Γ ( β ) , by (14) and the convolution theorem of Laplace transform, we have
S α , β ( t ) = 1 Γ ( β ) 0 t ( t s ) β 1 s α 1 Ψ ( s ) d s .
Lemma 7. 
If A generates a uniformly bounded and analytic semigroup, which satisfies (10) and (11), then
(i) There exist δ ( 0 , π 2 ) and a constant C > 0 such that
λ α | | arg λ | < π 2 + δ ρ ( A )
and
λ β R ( λ α , A ) C | λ | α + β .
(ii) A generates a ( α , β ) resolvent family { S α , β ( t ) } t 0 , and
S α , β ( t ) = 1 2 π i Γ e λ t λ β R ( λ α , A ) d λ ,
where
Γ = { r e i ( π 2 + δ 1 ) | r [ ρ , ) } { ρ e i δ 2 | | δ 2 | π 2 + δ 1 } { r e i ( π 2 + δ 1 ) | r [ ρ , ) } ,
δ 1 ( 0 , δ ) , δ ( 0 , π 2 ) and ρ > 0 . Moreover,
| | S α , β ( t ) | | M t α + β 1 , t > 0 ,
M > 0 is a constant.
Proof. (i) Let K ( λ ) = λ α , λ = r e i ( π 2 + δ ) with δ ( 0 , π 2 ) and r > 0 , then
| arg ( K ( λ ) ) | = | Im ( ln ( K ( r e i ( π 2 + δ ) ) ) ) | = Im 0 π 2 + δ d d t ln ( K ( r e i t ) ) d t = Im 0 π 2 + δ K ( r e i t ) i r e i t K ( r e i t ) d t α ( π 2 + δ ) < π 2 + δ .
Thus λ α ρ ( A ) , and
λ β R ( λ α , A ) C | λ | α + β ,
where C > 0 is a constant.
(ii) For t > 0 , and δ 1 ( 0 , δ ) , we set
S α , β ( t ) = 1 2 π i Γ e λ t λ β R ( λ α , A ) d λ ,
where
Γ = Γ 1 Γ 2 Γ 3 ,
Γ 1 = { r e i ( π 2 + δ 1 ) | r [ 1 t , ) } ,
Γ 2 = { 1 t e i θ | | θ | π 2 + δ 1 } ,
Γ 3 = { r e i ( π 2 + δ 1 ) | r [ 1 t , ) }
are oriented counterclockwise. From (16) it is easy to see that for t > 0 the integral in (18) converges in the uniform topology. Moreover,
1 2 π i Γ 3 e λ t λ β R ( λ α , A ) d λ ) C 2 π 1 t e r t sin δ 1 r α + β d r C ( sin δ 1 ) α + β 1 2 π sin δ 1 e τ τ α + β d τ t α + β 1 C 1 t α + β 1 .
Similarly, the integral on Γ 1 has the same estimation and on Γ 2 we get
1 2 π i Γ 2 e λ t λ β R ( λ α , A ) d λ ) C 2 π ( π 2 + δ ) π 2 + δ e cos θ d θ t α + β 1 C 2 t α + β 1 .
Hence, (17) holds.
Next we fix λ > 0 , we have
0 T e λ t S α , β ( t ) d t = 1 2 π i Γ μ β R ( μ α , A ) 0 T e ( λ μ ) t d t d μ = λ β R ( λ α , A ) + 1 2 π i Γ e ( λ μ ) T μ β R ( μ α , A ) μ λ d μ ,
where the Cauchy integral formula and Fubini theorem are also used. Due to
Γ e ( λ μ ) T μ β R ( μ α , A ) μ λ d μ M e T λ Γ | d μ | | μ | α + β | λ μ | 0 , T .
Therefore, by taking the limit as T in (19), we have
0 e λ t S α , β ( t ) d t = λ β R ( λ α , A ) .
That is, { S α , β ( t ) } t 0 is generated by A. □
Remark 3. 
Due to different parameters, the method in [19] cannot be directly applied to this paper. By comparing the forms of mild solutions to the studied equations, we find that (17) agrees well with Ψ ( t ) M ω Γ ( α ) in [24].
Lemma 8. 
For t > 0 , R α , β ( t ) = t 1 α β S α , β ( t ) is continuous in the uniform operator topology, where S α , β ( t ) is the ( α , β ) resolvent family generated by A.
Proof. 
Let ϵ > 0 be fixed. Due to the fact that Ψ ( t ) is continuous in the uniform operator topology for t > 0 [39], for arbitrary t 0 > 0 , there exists δ > 0 such that
Ψ ( t 2 ) Ψ ( t 1 ) Γ ( α + β ) Γ ( α ) ϵ ,
for t 2 > t 1 t 0 , and | t 2 t 1 | < δ .
Then, owing to
R α , β ( t 2 ) x R α , β ( t 1 ) x = 1 Γ ( β ) t 2 1 α β 0 t 2 ( t 2 s ) β 1 s α 1 Ψ ( s ) d s t 1 1 α β 0 t 1 ( t 1 s ) β 1 s α 1 Ψ ( s ) d s = 1 Γ ( β ) 0 1 ( 1 τ ) β 1 τ α 1 Ψ ( t 2 τ ) Ψ ( t 1 τ ) d τ ,
and
0 1 ( 1 τ ) β 1 τ α 1 d τ = B ( α , β ) = Γ ( α ) Γ ( β ) Γ ( α + β ) ,
where B ( · , · ) is the Beta function [2], we conclude that
R α , β ( t 2 ) R α , β ( t 1 ) ϵ .
That is, by the arbitrariness of t 0 , R α , β ( t ) is continuous in the uniform operator topology for t > 0 . □
Remark 4. 
{ R α , β ( t ) } t 0 is strongly continuous. That is, for arbitrary x X and 0 t 1 < t 2 T , we have
R α , β ( t 2 ) x R α , β ( t 1 ) x 0 as t 1 t 2 .
Proof. 
Since { Ψ ( t ) } t 0 is strongly continuous [40], there exist δ 1 > 0 such that
Ψ ( t 2 ) x Ψ ( t 1 ) x Γ ( α + β ) Γ ( α ) ϵ ,
for t 2 > t 1 0 , and | t 2 t 1 | < δ 1 .
Due to
R α , β ( t 2 ) x R α , β ( t 1 ) x = 1 Γ ( β ) 0 1 ( 1 τ ) β 1 τ α 1 Ψ ( t 2 τ ) x Ψ ( t 1 τ ) x d τ ϵ ,
we obtain that { R α , β ( t ) } t 0 is strongly continuous. □
It is noteworthy that the strongly continuity of { R α , β ( t ) } t 0 can not be obtained immediately by Lemma 8, in which t > 0 , not t 0 .
Lemma 9. 
If the analytic semigroup { T ( t ) } t 0 generated by A is compact, then R α , β ( t ) is compact for t > 0 .
Proof. 
Set B k = { x X | x k } . In order to show R α , β ( t ) is compact for t > 0 , we need to show that
{ R α , β ( t ) x | x B k }
is relatively compact in X, for any k > 0 and t > 0 .
Let t > 0 be fixed. For any ϵ > 0 , define
R ϵ ( t ) : = α Γ ( β ) t 1 α β 0 t ( t s ) β 1 s α 1 ϵ θ ζ α ( θ ) T ( t α θ ) d θ d s .
Then
R ϵ ( t ) = T ( t α ϵ ) α Γ ( β ) t 1 α β 0 t ( t s ) β 1 s α 1 ϵ θ ζ α ( θ ) T ( t α θ t α ϵ ) d θ d s
and
α Γ ( β ) t 1 α β 0 t ( t s ) β 1 s α 1 ϵ θ ζ α ( θ ) T ( t α θ t α ϵ ) d θ d s N Γ ( β ) 0 t ( 1 τ ) β 1 τ α d τ · α 0 θ ζ α ( θ ) d θ = N Γ ( β ) B ( α , β ) 1 Γ ( α ) = N Γ ( α + β ) ,
where B ( · , · ) is the Beta function. Due to the compactness of T ( t α ϵ ) ( t α ϵ > 0 ) , we obtain that { R ϵ ( t ) x | x B k } is relatively compact in X for arbitrary ϵ > 0 .
In addition, for any x B k , we have
R α , β ( t ) x R ϵ ( t ) x N K Γ ( α + β ) α 0 ϵ θ ζ α ( θ ) d θ .
Therefore, we obtain that { R α , β ( t ) x | x B k } is relatively compact in X. □
Remark 5. 
If B is the infinitesimal generator of a C 0 -semigroup { S ( t ) } t 0 on a Banach space, then S ( t ) is compact if and only if S ( t ) is continuous in the uniform operator topology for t > 0 and R ( λ , A ) is compact for λ ρ ( A ) [38]. If α ( 0 , 1 ) and β = 1 α , S α , β is compact if and only if A generates a compact C 0 semigroup, which is also obtained in [40].
Lemma 10. 
Assume that u C ( J , X ) , D t α u C ( ( 0 , T ] , X ) , u ( t ) D ( A ) for t ( 0 , T ] , and u satisfies (12). Then u satisfies the formal integral equation
u ( t ) = S α , 1 α ( t ) ( u 0 g ( u ) ) + 0 t S α , β ( t s ) f ( s , u ( s ) ) d s .
Proof. 
By the definitions of Caputo derivative and Riemann-Liouville integral [2], one can rewrite (12) as in the equivalent integral equation
u ( t ) = u 0 + 1 Γ ( α ) 0 t ( t s ) α 1 A u ( s ) d s + 1 Γ ( α + β ) 0 t ( t s ) α + β 1 f ( s , u ( s ) ) d s .
For λ Σ π 2 + δ , using the Laplace transform
u ^ ( λ ) = 0 e λ s u ( s ) d s and f ^ ( λ , u ( λ ) ) = 0 e λ s f ( s , u ( s ) ) d s
to (12), we have
λ α u ^ ( λ ) λ α 1 u ( 0 ) = A u ^ ( λ ) + λ β f ^ ( λ , u ( λ ) ) .
Then (23) is equivalent to
u ^ ( λ ) = λ α 1 ( λ α I A ) 1 u ( 0 ) + λ β ( λ α I A ) 1 f ^ ( λ , u ( λ ) ) .
By the inverse Laplace transform and Definition 6, we have
u ( t ) = S α , 1 α ( t ) ( u 0 g ( u ) ) + 0 t S α , β ( t s ) f ( s , u ( s ) ) d s ,
where t J . □
Consequently, we give the definition of a mild solution to (12) as follows.
Definition 11. 
The function u C ( J , X ) is called a mild solution of equation (12) if
u ( t ) = S α , 1 α ( t ) ( u 0 g ( u ) ) + 0 t S α , β ( t s ) f ( s , u ( s ) ) d s .
Lemma 12.([41])The convex closure conv ( Z ) ¯ is compact provided Z is a compact subset of a Banach space.
Lemma 13.([42])Assume that H : Y Y is completely continuous, where Y is a convex subset of a Banach space and 0 Y . Then either there is a fixed point of H or the set y Y | y = μ H ( y ) is unbounded, where μ ( 0 , 1 ) .
Lemma 14.([43]) Suppose that D is a bounded, convex and closed subset of a Banach space, 0 D , N : D D is continuous. If V = conv ¯ N ( V ) or V = N ( V ) { 0 } can obtain ν ( V ) = 0 for every subset V of D, then N has a fixed point.

3. Main results

In order to obtain the main results, we should introduce some hypotheses.
( H 1 ) f ( t , z ) is continuous with respect to z X for almost all t J , and strongly measurable with respect to t J for any z X ;
( H 1 ) f ( t , z ) is strongly measurable with respect to any z X and almost all t J ;
( H 2 ) g : C ( J , X ) X is completely continuous, and there exists a constant L > 0 such that
g ( u ) L u
for any u C ( J , X ) ;
( H 2 ) For g : C ( J , X ) X , there exists a constant L > 0 such that
g ( u ) L u
for any u C ( J , X ) ;
( H 3 ) There exists a continuous function φ : J R + such that
f ( t , z ) φ ( t ) z ,
for a.e. t J and each z X ;
( H 4 ) T ( t ) is compact for t > 0 ;
( H 5 ) For any bounded subset X 1 of X and each t J , there exists a 1 > 0 such that
ν ( g ( X 1 ) ) a 1 ν ( X 1 ) ,
and ν ( f ( t , X 1 ) ) φ ( t ) ν ( X 1 ) , where φ ( t ) is defined as in ( H 3 );
( H 6 ) There exist L g > 0 such that
g ( y 1 ) g ( y 2 ) L g y 1 y 2 ,
and f ( t , y 1 ( t ) ) f ( t , y 2 ( t ) ) φ ( t ) ( t ) y 1 y 2 , where φ ( t ) is defined as in ( H 3 ), t J , y 1 , y 2 B k 2 , and
B k 2 = u C ( J , X ) | u k 2 , k 2 = M ( u 0 + g ( 0 ) ) ( α + β ) ( α + β ) ( 1 M L g ) M T α + β φ .
Theorem 15. 
If ( H 1 )-( H 4 ) are satisfied, then (12) has a mild solution provided that
( α + β ) ( 1 M L ) > M T α + β φ .
Proof. 
Let
k 1 = M u 0 ( α + β ) ( α + β ) ( 1 M L ) M T α + β φ .
For any u B k 1 , by Lemma 7 and ( H 2 ), we have
S α , 1 α ( t ) ( u 0 g ( u ) ) M ( u 0 + k 1 L ) .
Furthermore, due to ( H 1 ), f ( t , u ( t ) ) is a measurable function on J. By Lemma 7 and ( H 3 ), we get
0 t S α , β ( t s ) f ( s , u ( s ) ) d s M 0 t ( t s ) α + β 1 φ ( s ) u ( s ) d s M k 1 0 t ( t s ) α + β 1 φ ( s ) d s M T α + β k 1 α + β φ ,
Then, S α , β ( t s ) f ( s , u ( s ) ) is Lebesgue integrable with respect to s J and t J , which implies that S α , β ( t s ) f ( s , u ( s ) ) is Bochner integrable with respect to s J and t J because of Bochner’s theorem.
Now we can define an operator Q on B k 1 as follows:
( Q u ) ( t ) = S α , 1 α ( t ) ( u 0 g ( u ) ) + 0 t S α , β ( t s ) f ( s , u ( s ) ) d s .
We firstly prove that Q is a completely continuous operator. Suppose that
u n , u B k 1 such that u n u as n ,
then
Q u n ( t ) Q u ( t ) M g ( u ) g ( u n ) + M 0 t ( t s ) α + β 1 f ( s , u n ( s ) ) f ( s , u ( s ) ) d s = : I 1 + I 2 .
Obviously, I 1 0 as n by (27) and ( H 2 ). Since ( H 3 ) we have
I 2 M 0 t ( t s ) α + β 1 ( f ( s , u n ( s ) ) + f ( s , u ( s ) ) ) d s 2 M k 1 0 t ( t s ) α + β 1 φ ( s ) d s 2 M k 1 T α + β α + β φ .
Then by Lebesgue dominated convergence theorem and ( H 1 ), we get I 2 0 as n . Therefore,
Q u n ( t ) Q u ( t ) 0 as n .
That is, Q is continuous on B k 1 .
Next, we prove that { Q u | u B k 1 } is relatively compact. It suffices to show that { Q u | u B k 1 } is uniformly bounded and equicontinuous, and { Q u ( t ) | u B k 1 } is relatively compact in X for any t J .
(25) and (26) imply
( Q u ) ( t ) k 1 ,
which means { Q u | u B k 1 } is uniformly bounded. Take u B k 1 and 0 t 1 < t 2 T , then Q u ( t 2 ) Q u ( t 1 ) I 3 + I 4 + I 5 , where
I 3 = [ S α , 1 α ( t 2 ) S α , 1 α ( t 1 ) ] [ u 0 g ( u ) ] , I 4 = t 1 t 2 S α , β ( t 2 s ) f ( s , u ( s ) ) d s , I 5 = 0 t 1 [ S α , β ( t 2 s ) S α , β ( t 1 s ) ] f ( s , u ( s ) ) d s .
It is easy to see that I 3 0 as t 1 t 2 by Remark 4. By ( H 3 ) and (17), we have
I 4 M t 1 t 2 ( t 2 s ) α + β 1 φ ( s ) u ( s ) d s M k 1 φ α + β ( t 2 t 1 ) α + β ,
which implies that I 4 0 as t 1 t 2 .
I 5 k 1 φ 0 t 1 ( t 2 s ) α + β 1 R α , β ( t 2 s ) R α , β ( t 1 s ) d s + M k 1 φ 0 t 1 ( t 1 s ) α + β 1 ( t 2 s ) α + β 1 d s k 1 φ 0 t 1 ϵ ( t 2 s ) α + β 1 d s sup s [ 0 , t 1 ϵ ] R α , β ( t 2 s ) R α , β ( t 1 s ) + 2 M k 1 φ t 1 ϵ t 1 ( t 2 s ) α + β 1 d s + M k 1 φ α + β t 1 α + β t 2 α + β + ( t 2 t 1 ) α + β k 1 φ α + β t 2 α + β ( t 2 t 1 + ϵ ) α + β sup s [ 0 , t 1 ϵ ] R α , β ( t 2 s ) R α , β ( t 1 s ) + 2 M k 1 φ α + β ( t 2 t 1 + ϵ ) α + β ( t 2 t 1 ) α + β + M k 1 φ α + β ( t 2 t 1 ) α + β
Then I 5 0 as t 1 t 2 and ϵ 0 by Lemma 8 and (17). Now, we can concluded that { Q u | u B k 1 } is equicontinuous.
Obviously, due to ( H 2 ), { Q u ( 0 ) | u B k 1 } is relatively compact. By ( H 2 ), ( H 4 ) and Remark 5, we can similarly prove the compactness of { S α , 1 α ( t ) ( u 0 g ( u ) ) | u B k 1 } for t ( 0 , T ] . Due to ( H 1 ), ( H 4 ) and Lemma 9, U is compact for t ( 0 , T ] , then conv ( U ) ¯ is compact for t ( 0 , T ] by Lemma 12, where
U = { R α , β ( t s ) f ( s , u ( s ) ) | s [ 0 , t ) , u B k 1 } .
By Mean-Value Theorem for the Bochner integral [44],
0 t S α , β ( t s ) f ( s , u ( s ) ) d s t α + β α + β conv ( U ) ¯ , t ( 0 , T ] ,
As a consequence, { Q u ( t ) | u B k 1 } ¯ is compact in X for all t ( 0 , T ] . Then { Q u ( t ) | u B k 1 } is relatively compact in X for any t J .
By the Arzela-Ascoli theorem, { Q u | u B k 1 } is relatively compact. By combining with the continuity of Q, we conclude that Q : B k 1 B k 1 is completely continuous.
Set
M 1 = { u B k 1 , u = η Q u , η ( 0 , 1 ) } ,
let us prove the boundedness of M 1 . Apparently, 0 M 1 . For u M 1 , we have
u ( t ) η [ M ( u 0 + g ( u ) ) + M 0 t ( t s ) α + β 1 φ ( s ) u ( s ) d s η M ( u 0 + k 1 L ) + M T α + β k 1 φ α + β < M ( u 0 + k 1 L ) + M T α + β k 1 φ α + β
for t J . Therefore, Q has a fixed point by Lemma 13. That is, (12) has a mild solution. □
Theorem 16. 
Under the assumptions ( H 1 ), ( H 2 ), ( H 3 ) and ( H 5 ), (12) has a mild solution provided M a 1 + T α + β α + β φ < 1 and ( α + β ) ( 1 M L ) > M T α + β φ .
Proof. 
Similar to Theorem 15, Q : B k 1 B k 1 is continuous, { Q u | u B k 1 } is uniformly bounded and equicontinuous, where k 1 is defined as in Theorem 15. Let V B k 1 such that V Q ( V ) ¯ { 0 } , then v ( t ) = ν ( V ( t ) ) is continuous for any t J because of the boundness and equicontinuity of V. Then by Lemma 5, ( H 5 ) and (17), we have
v sup t J ν ( ( Q V ) ( t ) { 0 } ) sup t J ν ( ( Q V ) ( t ) ) ν ( g ( V ) ) sup t J S α , 1 α ( t ) + sup t J 0 t S α , β ( t s ) φ ( s ) ν ( V ( s ) ) d s M a 1 v + M sup t J 0 t ( t s ) α + β 1 φ ( s ) v ( s ) d s M ( a 1 + T α + β α + β φ ) v .
Since M ( a 1 + T α + β α + β φ ) < 1 , we have v = 0 , namely, v ( t ) = ν ( V ( t ) ) = 0 . Consequently, V ( t ) is relatively compact in X. By the Arzela-Ascoli theorem, { Q u | u B k 1 } is relatively compact, which means ν ( V ) = 0 . Thus the proof is completed by Lemma 14. □
Theorem 17. 
If ( H 1 ), ( H 3 ) and ( H 6 ) are satisfied, then (12) has a unique mild solution provided ( α + β ) ( 1 M L g ) > M T α + β φ .
Proof. 
Let
k 2 = M ( u 0 + g ( 0 ) ) ( α + β ) ( α + β ) ( 1 M L g ) M T α + β φ .
By ( H 6 ), we have
S α , 1 α ( t ) ( u 0 g ( u ) ) M ( u 0 + g ( 0 ) + k 2 L g ) ,
and g is continuous. Similar to Theorem 15, for u B k 2 , S α , 1 α ( t ) ( u 0 g ( u ) ) exists, and S α , β ( t s ) f ( s , u ( s ) ) is Bochner integrable with respect to s J and t J . We also have that Q maps B k 2 into itself. Apparently, we only need prove that Q has a unique fixed point on B k 2 . For any u , v B k 2 and t J , according to ( H 6 ) we have
( Q u ) ( t ) ( Q v ) ( t ) S α , 1 α ( t ) [ g ( v ) g ( u ) ] + 0 t S α , β ( t s ) f ( s , u ( s ) ) f ( s , v ( s ) ) d s M L g u v + M u v 0 t ( t s ) α + β 1 φ ( s ) d s ( M L g + M T α + β φ α + β ) u v .
Since M L g + M T α + β φ α + β < 1 , by Banach contraction principle, we conclude that Q has a unique fixed point on B k 2 , the proof is complete. □

4. Conclusions

This article investigates the nonlocal initial value problem for a class of fractional order integro-differential diffusion equations with Dirichlet boundary conditions. By selecting appropriate Banach spaces, we transform the original problem into an abstract Cauchy problem. At the same time, we prove that operator A can generate a resolvent family { S α , β ( t ) } t 0 , and provide proof that it is continuous in the uniform topology and strongly continuous. If the analytical semigroup { T ( t ) } t 0 generated by A is compact, then { S α , β ( t ) } t > 0 is also compact. Through the relationship between Laplace transform and convolution, we also establish the relationship between { S α , β ( t ) } t 0 and the probability density function ζ α ( θ ) . It is worth noting that since A in [23] is ω -sector ( ω <0) and the order of the equation is different from that in this paper, the proof of S α , β ( t ) is bounded cannot be directly applied in this paper. By the Laplace transform, the definition of a mild solution for (12) be given, and finally prove the existence and uniqueness of the mild solution through several fixed point theorems. In the future, we may focus on the following meaningful topics:
  • The condition ( H 4 ) is strict, we hope to obtain the existence of the solution in this paper without requiring T ( t ) to be compact;
  • For α ( 0 , 1 ) and β ( 0 , ) in this paper, the stability and periodicity of S α , β ( t ) ;
  • By Remark 2, we have known S α , β ( t ) can establisn a relationship with ζ α ( θ ) and T ( t ) , then we expect to investigate the regularity of solutions for (12) by { S α , β ( t ) } t 0 if f satisfies a certain degree of continuity.

Author Contributions

J.M.: Developed the concept, methodology, and revised the final paper. Z.Y.: Prepared the draft. Y.Z.: provided valuable feedback on this article. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Fundamental Research Funds for the Central Universities (31920230052), Natural Science Foundation of Gansu Province (21JR1RA201) and Innovation Team of Intelligent Computing and Dynamical System Analysis and Application of Northwest Minzu University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dugowson, S. Les différentielles métaphysiques: histoire et philosophie de la généralisation de l’ordre de la dérivation. Paris 13, 1994.
  2. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and applications of fractional differential equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  3. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives-Theory and Applications; Gordon and Breach: Amsterdam, 1993. [Google Scholar]
  4. Agrawal, O.P. A general formulation and solution scheme for fractional optimal control problems. Nonlinear Dynamics. 2004, 38, 323–337. [Google Scholar] [CrossRef]
  5. Duarte, F.B.M.; Machado, J.A.T. Chaotic phenomena and fractional-order dynamics in the trajectory control of redundant manipulators. Nonlinear Dynamics. 2002, 29, 315–342. [Google Scholar] [CrossRef]
  6. Lederman, C.; Roquejoffre, J.M.; Wolanski, N. Mathematical justification of a nonlinear integro-differential equation for the propagation of spherical flames. Annali di Matematica Pura ed Applicata. 2004, 183, 173–239. [Google Scholar] [CrossRef]
  7. Kulish, V.V.; Lage, J.L. Application of fractional calculus to fluid mechanics. J. Fluids Engineering. 2002, 124, 803–806. [Google Scholar] [CrossRef]
  8. Oldham, K.B. Fractional differential equations in electrochemistry. Advances in Engineering software. 2010, 41, 9–12. [Google Scholar] [CrossRef]
  9. Ma, W.; Zhao, Y.; Guo, L.; Chen, Y. Qualitative and quantitative analysis of the COVID-19 pandemic by a two-side fractional-order compartmental model. ISA transactions. 2022, 124, 144–156. [Google Scholar] [CrossRef]
  10. Ma, W.; Li, Z.; Ma, N. Synchronization of discrete fractional-order complex networks with and without unknown topology. Chaos: An Interdisciplinary Journal of Nonlinear Science. 2022, 32, 013112. [Google Scholar] [CrossRef]
  11. Gorenflo, R.; Mainardi, F. Fractals and fractional calculus in continuum mechanics; Springer: Vienna, Fractional calculus, 1997; Volume 378, pp. 291–348. [Google Scholar]
  12. Eidelman, S.D.; Kochubei, A.N. Cauchy problem for fractional diffusion equations. Journal of differential equations. 2004, 199, 211–255. [Google Scholar] [CrossRef]
  13. Metzler, R.; Klafter, J. Boundary value problems for fractional diffusion equations. Physica A: Statistical Mechanics and its Applications. 2000, 278, 107–125. [Google Scholar] [CrossRef]
  14. Wyss, W. The fractional diffusion equation. Journal of Mathematical Physics. 1986, 27, 2782–2785. [Google Scholar] [CrossRef]
  15. Mohammed, W.W.; Alshammari, M.; Cesarano, C.; Albadrani, S.; Morshedy, M. Brownian motion effects on the stabilization of stochastic solutions to fractional diffusion equations with polynomials. Mathematics. 2022, 10, 1458. [Google Scholar] [CrossRef]
  16. Sin, C.S. Cauchy problem for general time fractional diffusion equation. Fractional Calculus and Applied Analysis. 2020, 23, 1545–1559. [Google Scholar] [CrossRef]
  17. Tuan, N.H.; Ngoc, T.B.; Zhou, Y.; O’Regan, D. On existence and regularity of a terminal value problem for the time fractional diffusion equation. Inverse Problems. 2020, 36, 055011. [Google Scholar] [CrossRef]
  18. Byszewski, L.; Lakshmikantham, V. Theorem about the existence and uniqueness of a solution of a nonlocal abstract Cauchy problem in a Banach space. Applicable analysis. 1991, 40, 11–19. [Google Scholar] [CrossRef]
  19. Chang, Y.K.; Ponce, R. Mild solutions for a multi-term fractional differential equation via resolvent operators. AIMS Mathematics. 2021, 6, 2398–2417. [Google Scholar] [CrossRef]
  20. Dien, N.M. On mild solutions of the generalized nonlinear fractional pseudo-parabolic equation with a nonlocal condition. Fractional Calculus and Applied Analysis. 2022, 25, 559–583. [Google Scholar] [CrossRef]
  21. Dubey, S.; Sharma, M. Solutions to fractional functional differential equations with nonlocal conditions. Fractional Calculus and Applied Analysis. 2014, 17, 654–673. [Google Scholar] [CrossRef]
  22. Gou, H.; Li, Y. Extremal mild solutions to Hilfer evolution equations with non-instantaneous impulses and nonlocal conditions. Fractional Calculus and Applied Analysis. 2023, 26, 1145–1185. [Google Scholar] [CrossRef]
  23. Ponce, R. Asymptotic behavior of mild solutions to fractional Cauchy problems in Banach spaces. Applied Mathematics Letters. 2020, 105, 106322. [Google Scholar] [CrossRef]
  24. Mu, J.; Ahmad, B.; Huang, S. Existence and regularity of solutions to time-fractional diffusion equations. Computers & Mathematics with Applications. 2017, 73, 985–996. [Google Scholar]
  25. Amin, R.; Shah, K.; Asif, M.; Khan, I.; Ullah, F. An efficient algorithm for numerical solution of fractional integro-differential equations via Haar wavelet. Journal of Computational and Applied Mathematics. 2021, 381, 113028. [Google Scholar] [CrossRef]
  26. Ahmadova, A.; Huseynov, I.T.; Fernandez, A.; Mahmudov, N.I. Trivariate Mittag-Leffler functions used to solve multi-order systems of fractional differential equations. Communications in Nonlinear Science and Numerical Simulation. 2021, 97, 105735. [Google Scholar] [CrossRef]
  27. Bazhlekova, E. Completely monotone multinomial Mittag-Leffler type functions and diffusion equations with multiple time-derivatives. Fractional Calculus and Applied Analysis. 2021, 24, 88–111. [Google Scholar] [CrossRef]
  28. Can, N.H.; Luc, N.H.; Baleanu, D.; Zhou, Y.; Long, L.D. Inverse source problem for time fractional diffusion equation with Mittag-Leffler kernel. Advances in Difference Equations. 2020, 2020, 1–18. [Google Scholar] [CrossRef]
  29. Luchko, Y. Some uniqueness and existence results for the initial-boundary-value problems for the generalized time-fractional diffusion equation. Computers & Mathematics with Applications. 2010, 59, 1766–1772. [Google Scholar]
  30. Zhou, Y.; He, J.W. A Cauchy problem for fractional evolution equations with Hilfer’s fractional derivative on semi-infinite interval. Fractional Calculus and Applied Analysis. 2022, 25, 924–961. [Google Scholar] [CrossRef]
  31. Zhou, Y.; He, J.W. New results on controllability of fractional evolution systems with order α∈(1,2). Evolution Equations and Control Theory. 2021, 10, 491–509. [Google Scholar] [CrossRef]
  32. Zhou, Y.; Jiao, F. Nonlocal Cauchy problem for fractional evolution equations. Nonlinear analysis: real world applications. 2010, 11, 4465–4475. [Google Scholar] [CrossRef]
  33. Chen, P.; Zhang, X.; Li, Y. Existence and approximate controllability of fractional evolution equations with nonlocal conditions via resolvent operators. Fractional Calculus and Applied Analysis. 2020, 23, 268–291. [Google Scholar] [CrossRef]
  34. Kuratowski, K. Sur les espaces complets. Fundamenta Mathematicae. 1930, 15, 301–309. [Google Scholar] [CrossRef]
  35. Banaś, J.; Goebel, K. Measures of Noncompactness in Banach Spaces. Marcel Dekker: New York, 1980.
  36. Evans, L.C. Partial differential equations. American Mathematical Society, 2022.
  37. Yagi, A. Abstract parabolic evolution equations and their applications; Springer Science & Business Media, 2009.
  38. Pazy, A. Semigroups of linear operators and applications to partial differential equations; Springer Science & Business Media: Berlin, 2012. [Google Scholar]
  39. Wang, R.N.; Chen, D.H.; Xiao, T.J. Abstract fractional Cauchy problems with almost sectorial operators. Journal of Differential Equations. 2012, 252, 202–235. [Google Scholar] [CrossRef]
  40. Zhou, Y.; Jiao, F. Existence of mild solutions for fractional neutral evolution equations. Computers & Mathematics with Applications. 2010, 59, 1063–1077. [Google Scholar]
  41. Kaloshin, V.Y. An extension of the Artin-Mazur theorem. Annals of Mathematics. 1999, 150, 729–741. [Google Scholar] [CrossRef]
  42. Shanmugam, T.; Muthiah, M.; Radenović, S. Existence of positive solution for the eighth-order boundary value problem using classical version of Leray-Schauder alternative fixed point theorem. Axioms. 2019, 8, 1–13. [Google Scholar] [CrossRef]
  43. Mo¨nch, H. Boundary value problems for nonlinear ordinary differential equations of second order in Banach spaces. Nonlinear Analysis: Theory, Methods & Applications. 1980, 4, 985–999. [Google Scholar]
  44. Diestel, J.; Uhl, J.J. Vector measures. American Mathematical Society: Providence, RI, 1977.
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.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated