Next Article in Journal
Fuel Cell Fractional-Order Model via Electrochemical Impedance Spectroscopy
Next Article in Special Issue
Uncertainty Quantification of Random Microbial Growth in a Competitive Environment via Probability Density Functions
Previous Article in Journal
Two Stage Implicit Method on Hexagonal Grids for Approximating the First Derivatives of the Solution to the Heat Equation
Previous Article in Special Issue
Novel Techniques for a Verified Simulation of Fractional-Order Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Strongly Continuous Resolving Families of Operators for Fractional Distributed Order Equations

1
Department of Mathematical Analysis, Chelyabinsk State University, 454001 Chelyabinsk, Russia
2
Laboratory of Functional Materials, South Ural State University, 454080 Chelyabinsk, Russia
3
Department of Differential Equations, N.N. Krasovskii Institute of Mathematics and Mechanics of the Ural Branch of the Russian Academy of Sciences, 620108 Yekaterinburg, Russia
*
Author to whom correspondence should be addressed.
Submission received: 9 January 2021 / Revised: 24 February 2021 / Accepted: 26 February 2021 / Published: 2 March 2021
(This article belongs to the Special Issue Fractional Order Systems: Deterministic and Stochastic Analysis)

Abstract

:
The aim of this work is to find by the methods of the Laplace transform the conditions for the existence of a strongly continuous resolving family of operators for a linear homogeneous equation in a Banach space with the distributed Gerasimov–Caputo fractional derivative and with a closed densely defined operator A in the right-hand side. It is proved that the existence of a resolving family of operators for such equation implies the belonging of the operator A to the class C W ( K , a ) , which is defined here. It is also shown that from the continuity of a resolving family of operators at t = 0 the boundedness of A follows. The existence of a resolving family is shown for A C W ( K , a ) and for the upper limit of the integration in the distributed derivative not greater than 2. As corollary, we obtain the existence of a unique solution for the Cauchy problem to the equation of such class. These results are used for the investigation of the initial boundary value problems unique solvability for a class of partial differential equations of the distributed order with respect to time.

1. Introduction

In the last couple of decades, a new branch of the theory of differential equations has emerged as equations with distributed fractional derivatives. A distributed derivative is an integral of a fractional derivative with respect to its order. Apparently, the first works on distributed order derivatives were those by Nakhushev [1,2] (they are called continual derivatives here) and Caputo [3] (mean derivatives). With the development of the field of applications of fractional integro-differential calculus, such equations began to arise in various problems in viscoelasticity theory [4], kinetic theory [5], modeling of ultraslow diffusion [6], and other scientific areas (see [7] and references therein). Numerical methods for solving equations with distributed fractional derivatives have been very actively investigated the last twenty years (see, e.g., [8,9]). Among the works in which the qualitative properties of distributed order differential operators and corresponding equations are investigated, in addition to the mentioned works of Nakhushev [1,2], we note the monograph and other works of Pskhu [10,11], series of articles by Atanacković, Oparnica and Pilipović ([12] and others), and works by Kochubei (e.g., [6]).
Let Z be a Banach space and A : D A Z a linear closed densely defined operator, m 1 < b m N , ω : [ 0 , b ] C . Consider the Cauchy problem
z ( k ) ( 0 ) = z k , k = 0 , 1 , , m 1 ,
for the distributed order equation
0 b ω ( α ) D t α z ( t ) d α = A z ( t ) , t > 0 ,
with the Gerasimov–Caputo derivative D t α . In [13,14], this problem is researched for the case of a bounded operator A. The Cauchy problem and other initial value problems for Equation (2) and for an analogous equation with a linear degenerate operator at the distributed order derivative are studied in [15,16] under the condition of the sectoriality of operators, with the Riemann–Liouville or the Gerasimov–Caputo distributed fractional derivatives in the equation.
Sufficient and necessary conditions of an unbounded closed densely defined operator A for the existence of an analytic resolving family of operators for Equation (2) are found in [17,18,19]. Thus, an extension of the theorem on generators of analytic semigroups of operators to the case of distributed order equations is obtained. This allows investigating a unique solvability of problem (1) and (2). A theorem on perturbations of generators for analytic resolving family of operators for distributed order Equation (2) is proved. These results were applied to study of initial boundary value problems for some partial differential equations of a distributed order with respect to time.
In this work, we study by the methods of the Laplace transform the conditions of the existence of a strongly continuous resolving family of operators for distributed order Equation (2) in terms of the operator A. In the first section, properties of such families of operators are studied. Firstly, we consider some properties of certain functions arising after the acting of the Laplace transform on a distributed order derivative. Then, the definitions of a resolving family of operators for Equation (2) and a class C W ( K , a ) of linear closed densely defined operators are introduced. It is proved that under the condition of the existence of a resolving family of operators for the distributed order equation the operator A belongs to C W ( K , a ) , and the form of the Laplace transform of a resolving family is found. It is shown that in contrast to the fixed-order differential equations the primitive of a solution of (2) is not a solution of this equation. Therefore, we need to introduce the so-called k-resolving families of operators, which give solutions of the Cauchy problem z ( l ) ( 0 ) = z l , l = 0 , 1 , , m 1 , l k , z ( k ) ( 0 ) = z k to Equation (2). It is proved that, if there exists a resolving family of operators { S ( t ) L ( Z ) : t 0 } of (2), then there exist k-resolving families of operators for every k = 1 , 2 , , m 1 , which are explicitly expressed in terms of { S ( t ) L ( Z ) : t 0 } . It is also shown that continuity of resolving family at t = 0 implies the boundedness of the operator A.
In the second section, the proposed approximations for operators S ( t ) allow proving the existence of the resolving family for Equation (2) with A C W ( K , a ) and b ( 0 , 2 ] . As a corollary, we obtain the unique solvability theorem for problem (1) and (2) at b ( 0 , 2 ] . In the third section, the abstract results are applied to the study of initial boundary value problems for a class of partial differential equations of the distributed order with respect to time.

2. Properties of Resolving Families of Operators

Let Z be a Banach space. Denote R ¯ + : = R + { 0 } , R ¯ : = R { 0 } , for β > 0 , t > 0 , h : R + Z
J t β h ( t ) : = 1 Γ ( β ) 0 t ( t s ) β 1 h ( s ) d s
is the Riemann–Liouville integral. Let m 1 < α m N , D t m be the usual derivative of the mth order, D t α be the Gerasimov–Caputo derivative of h : R ¯ + Z [20,21,22]:
D t α h ( t ) : = D t m J t m α h ( t ) k = 0 m 1 h ( k ) ( 0 ) t k k ! .
By Lap ( Z ) , denote the set of functions h : R + Z , such that the Laplace transform, which will be denoted by h ^ or Lap [ h ] , is defined. The Laplace transform of the Gerasimov–Caputo derivative of the order α > 0 satisfies the equality (see, e.g., [23])
D t α h ^ ( λ ) = λ α h ^ ( λ ) k = 0 m 1 h ( k ) ( 0 ) λ α 1 k .
Denote by L ( Z ) the Banach space of all linear continuous operators from Z to Z , and by C l ( Z ) denote the set of all linear closed operators, densely defined in Z , acting in the space Z . We supply the domain D A of an operator A C l ( Z ) by the norm of its graph and, thus, we get the Banach space D A .
Consider the Cauchy problem
z ( 0 ) = z 0 , z ( k ) ( 0 ) = 0 , k = 1 , 2 , , m 1 ,
for a distributed order equation
0 b ω ( α ) D t α z ( t ) d α = A z ( t ) , t > 0 ,
where 0 m 1 < b m N , ω : [ 0 , b ] C , A C l ( Z ) . By a solution of problem (4), (5) we call such function z C m 1 ( R ¯ + ; Z ) C ( R + ; D A ) that there exists 0 b ω ( α ) D t α z ( t ) d α C ( R + ; Z ) and equalities (4) and (5) are fulfilled.
Denote W k ( λ ) : = k b ω ( α ) λ α d α , k = 0 , 1 , , m 1 , W ( λ ) : = W 0 ( λ ) .
Lemma 1.
Let b > 0 , ω L 1 ( 0 , b ; C ) . Then, W k , k = 0 , 1 , , m 1 , are analytic on the set C \ { λ C : Re λ 0 } .
Proof. 
This assertion is proved for W 0 in [17]; for other W k , it can be proved similarly. □
Lemma 2.
Let b > 0 , ω : ( 0 , b ) C be a piecewise continuous function, which be continuous from the left at the point b, ω ( b ) 0 . Then, for every k = 0 , 1 , , m 1 ,
ε ( 0 , b ) C , ϱ > 0 λ { λ C : | λ | ϱ } \ R ¯ | W k ( λ ) | | λ | b ε / C .
Proof. 
Since ω ( b ) 0 , we have Re ω ( α ) e i α arg λ 0 for all α ( c , b ) , where c ( 0 , b ) is close enough to b. Hence, we have for some C 1 = C 1 ( ε 1 ) > 0 and for all large enough | λ |
c b ω ( α ) λ α d α Re c b ω ( α ) λ α d α = c b Re ω ( α ) e i α arg λ | λ | α d α
= Re ω ( ξ ) e i ξ arg λ c b | λ | α d α = Re ω ( ξ ) e i ξ arg λ | λ | b | λ | c ln | λ | C 1 | λ | ε 1
for some ξ ( c , b ) and every ε 1 ( c , b ) . Thus, at ε 2 ( c , ε 1 ) for some C 2 > 0 and large enough | λ |
k b ω ( α ) λ α d α C 1 | λ | ε 1 | λ | c k c | ω ( α ) | d α C 2 | λ | ε 2 .
Take ε = b ε 2 , C = C 2 1 . □
Lemma 3.
Let b > 0 , ω : ( 0 , b ) C be a bounded a.e. on ( 0 , b ) function. Then, for all k , l = 0 , 1 , , m 1 , k > l ,
C , ϱ > 0 λ { λ C : | λ | ϱ } \ R ¯ | W k ( λ ) W l ( λ ) | C | λ | k .
Proof. 
We have for λ { λ C : | λ | ϱ } \ R ¯ at some ϱ > 0
| W k ( λ ) W l ( λ ) | = l k ω ( α ) λ α d α C 1 | λ | k | λ | l ln | λ | C | λ | k .
A family of operators { S ( t ) L ( Z ) : t 0 } is called resolving for Equation (5), if the next conditions are satisfied:
(i)
S ( t ) is strongly continuous at t 0 , S ( 0 ) = I .
(ii)
S ( t ) [ D A ] D A , S ( t ) A z = A S ( t ) z for all z D A , t 0 .
(iii)
S ( t ) z 0 is a solution of Cauchy problem (4) and (5) for all z 0 D A .
A resolving family of operators { S ( t ) L ( Z ) : t 0 } is called exponentially bounded, if there exist K 1 , a 0 , such that S ( t ) L ( Z ) K e a t for all t 0 .
Remark 1.
Concepts of the resolving family of operators are used in the study of first order equations [24,25,26] ( C 0 -continuous semigroup of operators), integro-differential equations [27], integral evolution equations [28], and fractional differential equations [23].
Denote N 0 : = N { 0 } . An operator A C l ( Z ) is called the operator of the class C W ( K , a ) for some K 1 , a 0 , if the following two conditions are fulfilled:
(i)
If Re λ > a , then W ( λ ) ρ ( A ) .
(ii)
If Re λ > a , then for all n N 0 d n d λ n W ( λ ) λ ( W ( λ ) I A ) 1 L ( Z ) K n ! ( Re λ a ) n + 1 .
Here, we suppose that ω L 1 ( 0 , b ; C ) ;, hence, due to Lemma 1 and the properties of an operator resolvent, there exist the derivatives in Condition (ii) at Re λ > a .
Theorem 1.
Let b > 0 , ω L 1 ( 0 , b ; C ) ; there exists a resolving family of operators { S ( t ) L ( Z ) : t 0 } for Equation (5), such that for all t 0 S ( t ) L ( Z ) K e a t at some K 1 , a 0 . Then, A C W ( K , a ) and there exists the Laplace transform S ^ ( λ ) = W ( λ ) λ ( W ( λ ) I A ) 1 at Re λ > a .
Proof. 
Let there exist a resolving family of operators { S ( t ) L ( Z ) : t 0 } for Equation (5), for all t 0 S ( t ) L ( Z ) K e a t at some K 1 , a 0 . Then, at Re λ > a , there exists S ^ ( λ ) . From Equation (5), due to Paragraphs (ii) and (iii) of the resolving family definition, we obtain at z 0 D A equalities
0 b ω ( α ) D t α S ( t ) z 0 d α = A S ( t ) z 0 = S ( t ) A z 0 ,
hence, due to the closedness of the operator A at Re λ > a S ^ ( λ ) [ D A ] D A , and by virtue of (3)
Lap 0 b ω ( α ) D t α S ( t ) z 0 d α ( λ ) = W ( λ ) S ^ ( λ ) z 0 W ( λ ) λ z 0 = S ^ ( λ ) A z 0 = A S ^ ( λ ) z 0 .
Therefore, the operator W ( λ ) I A : D A Z is bijective and S ^ ( λ ) = W ( λ ) λ ( W ( λ ) I A ) 1 for all λ C , such that Re λ > a , and at n N
d n d λ n W ( λ ) λ ( W ( λ ) I A ) 1 = ( 1 ) n 0 t n e λ t S ( t ) d t .
Therefore, Conditions (i) and (ii) of the definition of the class C W ( K , a ) are fulfilled. □
Remark 2.
Let b > 1 , z 0 ( t ) be a solution of Cauchy problem (4) and (5) and z 0 ( t ) L ( Z ) K e a t at some K 1 , a 0 . Then, it can be shown that
z 0 ( t ) = 1 2 π i c i c + i W ( λ ) λ ( W ( λ ) I A ) 1 e λ t z 0 d λ : = S 0 ( t ) z 0 , c > a ,
if there exist the inverse operators ( W ( λ ) I A ) 1 and integral (8). If z 1 ( t ) is such solution of the Cauchy problem z ( 0 ) = 0 , z ( 1 ) ( 0 ) = z 1 , z ( k ) ( 0 ) = 0 , k = 2 , 3 , , m 1 , for Equation (5), then
z 1 ( t ) = 1 2 π i c i c + i W 1 ( λ ) λ 2 ( W ( λ ) I A ) 1 e λ t z 1 d λ : = S 1 ( t ) z 1
J t 1 S 0 ( t ) z 1 = 1 2 π i c i c + i W ( λ ) λ 2 ( W ( λ ) I A ) 1 e λ t z 1 d λ , c > a .
Thus, we see that, in contrast to linear homogeneous equations of a fixed order, S 1 ( t ) J t 1 S 0 ( t ) . In other words, for a solution z ( t ) of Equation (5), J t 1 z ( t ) is not a solution, generally speaking. Therefore, we need to consider resolving operators families { S k ( t ) : t 0 } , k = 0 , 1 , , m 1 , corresponding to every initial problem
z ( l ) ( 0 ) = 0 , l = 0 , 1 , , m 1 , l k , z ( k ) ( 0 ) = z k .
A family of operators { S k ( t ) L ( Z ) : t 0 } , k { 1 , 2 , , m 1 } , is called k-resolving for Equation (5), if the next conditions are satisfied:
(i)
S k ( t ) is strongly continuous at t 0 .
(ii)
S k ( t ) [ D A ] D A , S k ( t ) A z = A S k ( t ) z for all z D A , t 0 .
(iii)
S k ( t ) z k is a solution problem (5) and (9) at every z k D A .
Theorem 2.
Let b > 1 , ω L 1 ( 0 , b ; C ) , W satisfy conditions (6) and (7); there exists a resolving family of operators { S ( t ) L ( Z ) : t 0 } for Equation (5), such that for all t 0 S ( t ) L ( Z ) K e a t at some K 1 , a 0 . Then, there exist k-resolving families { S k ( t ) L ( Z ) : t 0 } , k = 1 , 2 , , m 1 , of (5).
Proof. 
Theorem 1 implies that A C W ( k , a ) . Consider at k = 1 , 2 , , m 1 the functions
W k ( λ ) λ k W ( λ ) = 1 λ k W ( λ ) W k ( λ ) λ k W ( λ ) .
Due to (6) and (7)
W ( λ ) W k ( λ ) λ k W ( λ ) C 1 | λ | b ε .
Since we can choose such ε > 0 , that b ε > 1 , there exists the inverse Laplace transform
w k ( t ) : = c i c + i e λ t W k ( λ ) λ k W ( λ ) d λ , c > 0 , k = 1 , 2 , , m 1 .
Here, at k = 1 , we take into account that λ 1 = 1 ^ . These scalar functions w k are continuous at t 0 , since the integrals converge uniformly with respect to t on every segment. Consequently, the operator-valued functions
S k ( t ) = 0 t w k ( t s ) S ( s ) d s , k = 1 , 2 , , m 1 ,
are strongly continuous at t 0 and Paragraph (ii) of the definition of k-resolving family of operators is valid also.
Note that the derivatives
w k ( l ) ( t ) : = c i c + i e λ t W k ( λ ) λ k l W ( λ ) d λ , c > 0 , l = 1 , 2 , , k 1 , k = 1 , 2 , , m 1 .
are continuous at t 0 as well, since
W k ( λ ) λ k l W ( λ ) = 1 λ k l W ( λ ) W k ( λ ) λ k l W ( λ ) , W ( λ ) W k ( λ ) λ k l W ( λ ) C 1 | λ | b l ε .
Moreover, relations (10) imply that w k ( l ) ( 0 ) = 0 , l = 0 , 1 , , k 2 , w k ( k 1 ) ( 0 ) = 1 . Thus, for k = 1 , 2 , , m 1 , z D A ,
S k ( l ) ( t ) z = 0 t w k ( l ) ( t s ) S ( s ) z d s , l = 1 , 2 , , k 1 ,
S k ( l ) ( t ) z = S ( l k ) ( t ) z + 0 t w k ( k 1 ) S ( l k + 1 ) ( t s ) z d s , l = k , k + 1 , , m 1 ,
hence, S k ( t ) z k satisfies initial conditions (9).
We have at Re λ > a
S ^ k ( λ ) = w ^ k ( λ ) S ^ ( λ ) = W k ( λ ) λ k W ( λ ) W ( λ ) λ ( W ( λ ) I A ) 1 = W k ( λ ) λ k + 1 ( W ( λ ) I A ) 1 .
Therefore, at z k D A ,
Lap 0 b ω ( α ) D t α S k ( t ) z k d α ( λ ) = W ( λ ) W k ( λ ) λ k + 1 ( W ( λ ) I A ) 1 z k W k ( λ ) λ k + 1 z k
= A W k ( λ ) λ k + 1 ( W ( λ ) I A ) 1 z k = A S ^ k ( λ ) z k = S ^ k ( λ ) A z k .
Acting by the inverse Laplace transform, we obtain that S k ( t ) z k is a solution of Equation (5). □
Theorem 3.
Let b > 0 , ω L 1 ( 0 , b ; C ) , W satisfy condition (6); there exist an exponentially bounded resolving family of operators of Equation (5). This family is continuous in the norm of L ( Z ) , if and only if A L ( Z ) .
Proof. 
Take Re λ > a and obtain as in the proof of Lemma 4
W ( λ ) λ ( W ( λ ) I A ) 1 I λ = 0 e λ t ( S ( t ) I ) d t .
Let the function η ( t ) = S ( t ) I L ( Z ) be continuous on the segment [ 0 , 1 ] and η ( 0 ) = 0 . For ε > 0 , take δ > 0 such that η ( t ) ε for all t [ 0 , δ ] , then
W ( λ ) λ ( W ( λ ) I A ) 1 I λ L ( Z ) 0 δ e λ t η ( t ) d t + δ e λ t η ( t ) d t ε λ + o 1 λ
as Re λ + , since η ( t ) K e a t + 1 at t 0 . Consequently, for large enough λ
W ( λ ) ( W ( λ ) I A ) 1 I L ( Z < 1 ,
hence, the operator W ( λ ) ( W ( λ ) I A ) 1 is continuously invertible, [ W ( λ ) ( W ( λ ) I A ) 1 ] 1 = W ( λ ) 1 ( W ( λ ) I A ) L ( Z ) . Thus, A L ( Z ) .
Now, let A L ( Z ) . In this case, we perform the resolving family of operators { S ( t ) L ( Z ) : t 0 } as the inverse Laplace transform of W ( λ ) λ ( W ( λ ) I A ) 1 . Indeed, set Γ R = Γ 1 , R Γ 2 , R is a positively oriented closed loop, where
Γ 1 , R = { c + i y C : y [ R , R ] } , Γ 2 , R = { c + R e i φ C : φ ( π / 2 , π / 2 ) } ,
define also the contours Γ 3 , R = { c + i y C : y ( R , ) } , Γ 4 , R = { c + i y C : y ( , R ) } , then { c + i y C : y R } = Γ 3 , R Γ 4 , R Γ R \ Γ 2 , R . Consider for large enough c > 0 the operators at t 0
S ( t ) = 1 2 π i c i c + i e λ t W ( λ ) λ ( W ( λ ) A ) 1 d λ = I + 1 2 π i c i c + i n = 1 A n e λ t d λ λ W ( λ ) n
= I + 1 2 π i Γ R Γ 2 , R + Γ 3 , R + Γ 4 , R n = 1 A n e λ t d λ λ W ( λ ) n .
At ϱ > 0 from (6) for λ { c + i k C : c > a , c 2 + k 2 ϱ 2 }
A n λ W ( λ ) n L ( Z ) C A L ( Z ) n | λ | ( b ε ) n + 1 .
Take small t > 0 and R = 1 / t , then
Z 0 ( t ) I L ( Z ) C e 2 π k = 1 3 n = 1 Γ k , R A L ( Z ) n | λ | ( b ε ) n + 1 | d λ | C 1 n = 1 A L ( Z ) n R ( b ε ) n
= C 1 n = 1 t ( b ε ) n A L ( Z ) n = C 1 t b ε A L ( Z ) 1 t b ε A L ( Z ) 0
as t 0 + . □
Remark 3.
As shown in the proof of Theorem 3, if A L ( Z ) , then resolving operators of Equation (5) has the form
S ( t ) = n = 0 a n ( t ) A n ,
where
a 0 ( t ) 1 , a n ( t ) = 1 2 π i γ e λ t d λ λ W ( λ ) n , n N ,
γ = { R e i φ : φ ( π , π ) } { x e i π , x ( R , ) } { x e i π , x ( , R ) } at R > 0 large enough. For equation D t α z ( t ) = A z ( t ) , we have W ( λ ) = λ α , and we obtain using the Hankel representation for the Euler gamma function that
a n ( t ) = 1 2 π i γ e λ t d λ λ α n + 1 = t α n 2 π i t γ e μ d μ μ α n + 1 = t α n Γ ( α n + 1 ) , n N .
Therefore, S ( t ) = n = 0 t α n A n Γ ( α n + 1 ) = E α , 1 ( t α A ) is the Mittag–Leffler function.

3. Existence of Resolving Families of Operators and the Cauchy Problem Unique Solvability

Let A C W ( K , a ) , for brevity, denote
H ( λ ) : = W ( λ ) λ ( W ( λ ) I A ) 1 , Re λ > a ,
and define at n N , t 0
S n ( t ) : = e n t I + k = 0 ( 1 ) k ( n ( n + a ) t ) k + 1 k ! ( k + 1 ) ! H ( k ) ( n + a ) .
The series converges uniformly on every segment [ 0 , T ] and S n ( t ) L ( Z ) K e a t for all t 0 . It can be shown that S n ( t ) is infinitely differentiable at t > 0 . Put for Re λ > a
H n ( λ ) : = 0 e λ t S n ( t ) d t ,
then for s > a
S n ( t ) = s i s + i e λ t H n ( λ ) d λ .
Here and further, i is the imaginary unit.
Lemma 4.
Let b > 0 , ω L 1 ( 0 , b ; C ) , W satisfy condition (6); A C W ( K , a ) , r = [ 1 / b ] + 1 at b ( 0 , 1 ] , D A r is dense in Z . Then, for any z Z , there exists the limit lim n S n ( t ) z in Z .
Proof. 
It is sufficient to prove the thesis on a dense set in Z . For z D A
H ( λ ) z = z λ + H ( λ ) A z W ( λ ) ,
so that the integral
s i s + i e λ t H ( λ ) z d λ , s > a ,
converges. Indeed, the first term in the right-hand side of (12) is the Laplace transform of the constant z and the second term is absolutely integrable at b > 1 due to conditions (6) and A C W ( K , a ) . If b ( 0 , 1 ] , then r b > 1 and for z D A r
H ( λ ) z = z λ + A z λ W ( λ ) + + A r 1 z λ W ( λ ) r 1 + H ( λ ) A r z W ( λ ) r .
Thus, integral (13) converges, if we choose a small enough number ε > 0 in condition (6).
We have for Re λ > a
H n ( λ ) = 0 e ( n + λ ) t I + k = 0 ( 1 ) k ( n ( n + a ) t ) k + 1 k ! ( k + 1 ) ! H ( k ) ( n + a ) d t
= I n + λ + k = 0 ( 1 ) k ( n ( n + a ) ) k + 1 H ( k ) ( n + a ) k ! ( n + λ ) k + 2
= I n + λ + n ( n + a ) ( n + λ ) 2 k = 0 λ ( n + a ) n + λ ( n + a ) k H ( k ) ( n + a ) k ! = I n + λ + n ( n + a ) ( n + λ ) 2 H λ ( n + a ) n + λ ,
so that lim n H n ( λ ) = H ( λ ) since H ( λ ) is analytic on { λ C : Re λ > a } due to Lemma 1. Moreover, for Re λ > a , z D A due to (12)
H n ( λ ) z z λ = z λ + z n + λ + n n + λ z λ + ( n + a ) H ( λ ( n + a ) n + λ ) A z ( n + λ ) W ( λ ( n + a ) n + λ ) = n ( n + a ) H ( λ ( n + a ) n + λ ) A z ( n + λ ) 2 W ( λ ( n + a ) n + λ ) .
Therefore, for all Re λ > a , n N ,
H n ( λ ) z z λ L ( Z ) C 1 n | n + a | A z Z | n + λ | 2 | λ | n + a | n + λ | b ε Re ( λ ( n + a ) n + λ a ) C 2 | λ | b ε n + λ n + a b ε 2 C 3 | λ | 2 .
We can pass to the limit as n in (11) and the thesis is proved, since D A (or D A r at b ( 0 , 1 ] ) is dense in Z . Note that the strong convergence obtained is uniform with respect to t from every segment [ 0 , T ] . □
Put Z 0 ( t ) : = s lim n S n ( t ) , where s lim n denotes the limit in the strong topology.
Remark 4.
It is proven that, for z D A or z D A r at b ( 0 , 1 ] ,
Z 0 ( t ) z = s i s + i e λ t H ( λ ) z d λ , s > a .
Lemma 5.
Let b > 0 , ω L 1 ( 0 , b ; C ) , W satisfy condition (6), A C W ( K , a ) , r = [ 1 / b ] + 1 at b ( 0 , 1 ] , D A r is dense in Z . Then, for all t 0 Z 0 ( t ) L ( Z ) K e a t , Z ^ ( λ ) = H ( λ ) at Re λ > a .
Proof. 
Due to the construction Z 0 ( t ) we have the inequality Z 0 ( t ) L ( Z ) K e a t by virtue of such inequality for S n ( t ) at every n N . Equality (14) implies that for every z D A or z D A r at b ( 0 , 1 ] Z ^ 0 ( λ ) z = H ( λ ) z . It remains to continue this equality on Z due to uniform boundedness of { H ( λ ) L ( Z ) : Re λ > a } , which follows from the condition A C W ( K , a ) . □
Theorem 4.
Let b ( 0 , 2 ] , ω L 1 ( 0 , b ; C ) , W satisfy condition (6); r = [ 1 / b ] + 1 at b ( 0 , 1 ] , D A r is dense in Z . Then, there exists a resolving family of operators { S ( t ) : t 0 } for Equation (5), such that S ( t ) L ( Z ) K e a t for some K 1 , a 0 , if and only if A C W ( K , a ) .
Proof. 
The direct assertion of this theorem is proved in Theorem 1. Let an operator A C W ( K , a ) . From the proof of Lemma 4, it follows that the family { Z 0 ( t ) L ( Z ) : t 0 } is strongly continuous, and, due to Lemma 5, it is exponentially bounded.
The operators d n d λ n W ( λ ) λ ( W ( λ ) I A ) 1 can be expressed in terms of sums of ( W ( λ ) I A ) 1 powers multiplied by scalar functions, therefore, for every z D A
d n d λ n W ( λ ) λ ( W ( λ ) I A ) 1 z D A ,
A d n d λ n W ( λ ) λ ( W ( λ ) I A ) 1 z = d n d λ n W ( λ ) λ ( W ( λ ) I A ) 1 A z ,
consequently, S n ( t ) z D A , A S n ( t ) z = S n ( t ) A z . Hence, for z D A Z 0 ( t ) z D A , A Z 0 ( t ) z = Z 0 ( t ) A z .
We have at z Z Z 0 ( 0 ) z = lim n S n ( 0 ) z = z , and for b > 1 , z D A
Z 0 ( 0 ) z = lim n S n ( 0 ) z = lim n n [ ( n + a ) H ( n + a ) I ] z = lim n n ( n + a ) W ( n + a ) H ( n + a ) A z = 0
due to (6).
Take z ( t ) = Z 0 ( t ) z 0 at some z 0 D A ; then, z ^ ( λ ) = H ( λ ) z 0 due to Lemma 5. Substitute z ( t ) in in the left-hand side of Equation (5) and act by the Laplace transform on the obtained expression taking into account the equalities obtained above z ( 0 ) = z 0 and at b ( 1 , 2 ] z ( 0 ) = 0 . Therefore, we obtain W ( λ ) z ^ ( λ ) W ( λ ) λ z 0 = W ( λ ) 2 λ ( W ( λ ) I A ) 1 z 0 W ( λ ) λ z 0 = A W ( λ ) λ ( W ( λ ) I A ) 1 z 0 = A z ^ ( λ ) , hence, z ( t ) is a solution of problem (4) ad (5). Thus, { Z 0 ( t ) L ( Z ) : t 0 } is a resolving family of operators for Equation (5). □
Corollary 1.
Let b ( 0 , 1 ] , ω L 1 ( 0 , b ) , W satisfy condition (6); K 1 , a 0 , A C W ( K , a ) , r = [ 1 / b ] + 1 , D A r is dense in Z . Then, for every z 0 D A , there exists a unique in the space Lap ( Z ) solution of the problem z ( 0 ) = z 0 for Equation (5). The solution has the form z ( t ) = Z 0 ( t ) z 0 .
Proof. 
Due to the previous theorem, it remains to prove the uniqueness of the solution. If there exist two solutions z 1 , z 2 of the problem z ( 0 ) = z 0 for Equation (5) from the space Lap ( Z ) , then their difference y = z 1 z 2 Lap ( Z ) is a solution of Equation (5) and satisfies the initial condition y ( 0 ) = 0 . Acting by the Laplace transform on both parts of equation (5) and due to the initial condition we get the equality W ( λ ) y ^ ( λ ) = A y ^ ( λ ) . Since A C W ( K , a ) , at Re λ > a we get the identity y ^ ( λ ) 0 . It means that y 0 . Therefore, z ( t ) = Z 0 ( t ) z 0 is a unique solution at z 0 D A in Lap ( Z ) . □
Corollary 2.
Let b ( 1 , 2 ] , ω L 1 ( 0 , b ) , W, W 1 satisfy conditions (6) and (7), K 1 , a 0 , A C W ( K , a ) , z 0 , z 1 D A . Then, there exists a unique in the space Lap ( Z ) solution of the problem
z ( 0 ) = z 0 , z ( 1 ) ( 0 ) = z 1
for Equation (5). The solution has the form z ( t ) = Z 0 ( t ) z 0 + Z 1 ( t ) z 1 , where
Z 1 ( t ) : = 0 t w 1 ( t s ) Z 0 ( s ) d s .
Proof. 
Recall that the function
w 1 ( t ) : = c i c + i e λ t W 1 ( λ ) λ W ( λ ) d λ , c > 0 .
is defined in the proof of Theorem 2, where it is shown that the function Z 1 ( t ) z 1 is a solution of the Cauchy problem z ( 0 ) = 0 , z ( 1 ) ( 0 ) = z 1 for Equation (5).
Reasoning as in the proof of Corollary 1, it is easy to show the uniqueness of the solution in Lap ( Z ) . □
Remark 5.
If we consider problem (4) and (5) on a segment [ 0 , T ] , then we can continue the function y on ( T , ) by a continuous bounded way and, reasoning analogously, we get the uniqueness of the solution on the segment.

4. Application to a Class of Initial-Boundary Value Problems

Let polynomials P n ( λ ) = j = 0 n c j λ j 0 , Q n ( λ ) = j = 0 n d j λ j be given, c j , d j R , j = 0 , 1 , , n , d n 0 . Suppose that Ω R d is a bounded region with a smooth boundary Ω , operator pencil Λ , B 1 , B 2 , , B p is regularly elliptic [29], where
( Λ u ) ( s ) = | q | 2 p a q ( s ) | q | u ( s ) s 1 q 1 s 2 q 2 s d q d , a q C ( Ω ¯ ) ,
( B l u ) ( s ) = | q | p l b l q ( s ) | q | u ( s ) s 1 q 1 s 2 q 2 s d q d , b l q C ( Ω ) , l = 1 , 2 , , p ,
q = ( q 1 , q 2 , , q d ) N 0 d , | q | = q 1 + + q d . Define the operator Λ 1 C l ( L 2 ( Ω ) ) with the domain D Λ 1 = H { B l } 2 p ( Ω ) [29] by the equality Λ 1 u = Λ u . Suppose that the operator Λ 1 is self-adjoint; then, the spectrum σ ( Λ 1 ) of the operator Λ 1 is real and discrete [29]. Let, moreover, σ ( Λ 1 ) be bounded from the right and do not contain the origin and { φ k : k N } be an orthonormal in L 2 ( Ω ) system of eigenfunctions of the operator Λ 1 , numbered by non-increasing corresponding eigenvalues { λ k : k N } , taking into account their multiplicities.
Consider the initial-boundary value problem
u ( s , 0 ) = u 0 ( s ) , s Ω ,
B l Λ k u ( s , t ) = 0 , k = 0 , 1 , , n 1 , l = 1 , 2 , , p , ( s , t ) Ω × ( 0 , T ) ,
0 b ω ( α ) D t α P n ( Λ ) u ( s , t ) d α = i Q n ( Λ ) u ( s , t ) , ( s , t ) Ω × ( 0 , T ) ,
where 0 < b 1 , ω : ( 0 , b ) R , i is the imaginary unit. Set n 0 = max { j { 0 , 1 , , n } : c j 0 } ,
X = { v H 2 r n 0 ( Ω ) : B l Λ k v ( s ) = 0 , k = 0 , 1 , , n 0 1 , l = 1 , 2 , , p , x Ω } ,
Y = L 2 ( Ω ) , L = P n ( Λ ) , M = i Q n ( Λ ) ,
D M = { v H 2 r n ( Ω ) : B l Λ k v ( s ) = 0 , k = 0 , 1 , , n 1 , l = 1 , 2 , , p , x Ω } .
Then, L L ( X ; Y ) , M C l ( X ; Y ) (moreover, M L ( X ; Y ) , if n 0 = n , i. e. c n 0 ). Let P n ( λ k ) 0 for all k N ; then, there exists an inverse operator L 1 L ( Y ; X ) and problem (16)–(18) is represented as the Cauchy problem z ( 0 ) = z 0 for Equation (5), where Z = X , A = L 1 M C l ( Z ) , D A = D M , z 0 = u 0 ( · ) . We see that D A r is dense in Z , where r = [ 1 / b ] + 1 .
Lemma 6.
Let b ( 0 , 1 ] , ω C ( [ 0 , b ] ; R ) , ω ( α ) 0 at α [ 0 , b ) , ω ( b ) > 0 , the spectrum σ ( Λ 1 ) does not contain the origin and zeros of the polynomial P n ( λ ) , and designations (19)–(21) are valid. Then, A C W ( K , a ) .
Proof. 
Due to Lemma 1, the function W ( λ ) : = 0 b ω ( α ) λ α d α is analytic on C \ { λ C : Re λ 0 } . If Re λ > 1 , then W ( λ ) { Re λ 1 } , since ω ( α ) 0 a.e. on ( 0 , b ) and α 1 .
For v X , Re λ > 0 ,
( W ( λ ) I A ) 1 v X 2 = ( W ( λ ) I i P n ( Λ ) 1 Q n ( Λ ) ) 1 X 2 = k = 1 ( 1 + λ k 2 n 0 ) | v , φ k | 2 | W ( λ ) i Q n ( λ k ) P n ( λ k ) | 2
k = 1 ( 1 + λ k 2 n 0 ) | v , φ k | 2 | W ( λ ) | 2 | 1 i W ( λ ) 1 Q n ( λ k ) P n ( λ k ) | 2 k = 1 ( 1 + λ k 2 n 0 ) | v , φ k | 2 | W ( λ ) | 2 min k N , Re λ > 1 | 1 i W ( λ ) 1 Q n ( λ k ) P n ( λ k ) | 2 v X 2 | W ( λ ) | 2 .
Here, the evident inclusion i W ( λ ) 1 Q n ( λ k ) P n ( λ k ) { Im λ 1 } { Im λ 1 } is used. Thus, we have the inequality λ H ( λ ) L ( Z ) 1 .
Furthermore,
H ( λ ) = W ( λ ) λ 2 ( W ( λ ) I A ) 1 + W ( λ ) λ ( W ( λ ) I A ) 1 W ( λ ) W ( λ ) λ ( W ( λ ) I A ) 2 ,
λ 2 H ( λ ) L ( Z ) W ( λ ) ( W ( λ ) I A ) 1 L ( Z ) + λ W ( λ ) ( W ( λ ) I A ) 1 L ( Z ) +
+ λ W ( λ ) W ( λ ) ( W ( λ ) I A ) 2 L ( Z ) 3 ,
since we can repeat the first part of the proof with λ W ( λ ) instead of W ( λ ) . Consequently, by the result of Prüss (Theorem 0.4 at k = 0 [28] and Corollary 2.5.2 [30]) we have A C W ( K , 1 ) . □
Lemma 6 and Corollary 1 imply the next unique solvability theorem.
Theorem 5.
Let b ( 0 , 1 ] , ω C ( [ 0 , b ] ; R ) , ω ( α ) 0 at α [ 0 , b ) , ω ( b ) > 0 , the spectrum σ ( Λ 1 ) do not contain the origin and zeros of the polynomial P n ( λ ) , u 0 D M . Then, there exists a unique solution of problem (16)–(18).
Remark 6.
From the proof of Lemma 6, it follows that we can obtain the analogous result for a case b ( 1 , 2 ] , if we provide the positivity of the value min k N , Re λ > a | 1 i W ( λ ) 1 Q n ( λ k ) P n ( λ k ) | at some a > 0 .
Example 1.
Take P 1 ( λ ) 1 , Q 1 ( λ ) = λ , A u = Δ u , p = 1 , B 1 = I . Then, (16)–(18) is the initial-boundary value problem for the modified equation of the ultraslow diffusion [6], which here can be interpreted as equation of ultraslow dynamics of non-relativistic systems:
0 b ω ( α ) D t α u ( s , t ) d α = i Δ u ( s , t ) , ( s , t ) Ω × ( 0 , T ) ,
u ( s , t ) = 0 , ( s , t ) Ω × ( 0 , T ) ,
u ( s , 0 ) = u 0 ( s ) , s Ω .

5. Conclusions

In the present work, it is proved that a sufficient (for b ( 0 , 2 ] ) and necessary (for a general case) condition of the existence of a strongly continuous resolving family of operators for differential equation of a distributed order (2) is A C W ( K , a ) . The approximation of such family of operators is proposed. It is proved that, if there exists a resolving family of operators { S ( t ) L ( Z ) : t 0 } of (2), then there exist the so-called k-resolving families of operators for every k = 1 , 2 , , m 1 , which give solutions of Cauchy problem z ( l ) ( 0 ) = z l , l = 0 , 1 , , m 1 , l k , z ( k ) ( 0 ) = z k to Equation (2). They are explicitly expressed in terms of { S ( t ) L ( Z ) : t 0 } by a non-trivial way, not in the form of the primitive as in the case of fixed-order differential equations. Further these results can be used for consideration of the case b > 2 and the study of the unique solvability of the Cauchy problem for the corresponding inhomogeneous equation.

Author Contributions

Conceptualization, V.E.F.; methodology, V.E.F.; software, N.V.F.; validation, V.E.F.; formal analysis, N.V.F.; investigation, V.E.F. and N.V.F.; resources, N.V.F.; data curation, N.V.F.; writing–original draft preparation, V.E.F.; writing—review and editing, V.E.F.; visualization, N.V.F.; supervision, V.E.F.; project administration, V.E.F.; and funding acquisition, V.E.F. All authors have read and agreed to the published version of the manuscript.

Funding

The reported study was funded by the Russian Foundation for Basic Research, project number 21-51-54003; and by Act 211 of Government of the Russian Federation, contract 02.A03.21.0011. It was also performed in the framework of Ural Mathematical Center research.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Nakhushev, A.M. On continual differential equations and their difference analogues. Sov. Math. Dokl. 1988, 37, 729–732. [Google Scholar]
  2. Nakhushev, A.M. Positiveness of the operators of continual and discrete differentiation and integration, which are quite important in the fractional calculus and in the theory of mixed-type equations. Differ. Equ. 1998, 34, 103–112. [Google Scholar]
  3. Caputo, M. Mean fractional order derivatives. Differential equations and filters. Annali dell’Univ. Ferrara Sez. VII. Sci. Mat. 1995, 41, 73–84. [Google Scholar]
  4. Lorenzo, C.F.; Hartley, T.T. Variable order and distributed order fractional operators. Nonlinear Dyn. 2002, 29, 57–98. [Google Scholar] [CrossRef]
  5. Sokolov, I.M.; Chechkin, A.V.; Klafter, J. Distributed-order fractional kinetics. Acta Phys. Pol. B 2004, 35, 1323–1341. [Google Scholar]
  6. Kochubei, A.N. Distributed order calculus and equations of ultraslow diffusion. J. Math. Anal. Appl. 2008, 340, 252–280. [Google Scholar] [CrossRef] [Green Version]
  7. Caputo, M.; Fabrizio, M. The kernel of the distributed order fractional derivatives with an application to complex materials. Fractal Fracttionals 2017, 1, 13. [Google Scholar] [CrossRef] [Green Version]
  8. Diethelm, K.; Ford, N.J. Numerical solution methods for distributed order time fractional diffusion equation. Fract. Calc. Appl. Anal. 2001, 4, 531–542. [Google Scholar]
  9. Diethelm, K.; Ford, N.; Freed, A.D.; Luchko, Y. Algorithms for the fractional calculus: A selection of numerical methods. Comput. Methods Appl. Mech. Eng. 2003, 194, 743–773. [Google Scholar] [CrossRef] [Green Version]
  10. Pskhu, A.V. On the theory of the continual and integro-differentiation operator. Differ. Equ. 2004, 40, 128–136. [Google Scholar] [CrossRef]
  11. Pskhu, A.V. Partial Differential Equations of Fractional Order; Nauka Publ.: Moscow, Russia, 2005. (In Russian) [Google Scholar]
  12. Atanacković, T.M.; Oparnica, L.; Pilipović, S. On a nonlinear distributed order fractional differential equation. J. Math. Anal. Appl. 2007, 328, 590–608. [Google Scholar] [CrossRef]
  13. Streletskaya, E.M.; Fedorov, V.E.; Debbouche, A. The Cauchy problem for distributed order equations in Banach spaces. Math. Notes NEFU 2018, 25, 63–72. [Google Scholar]
  14. Fedorov, V.E.; Streletskaya, E.M. Initial-value problems for linear distributed-order differential equations in Banach spaces. Electron. J. Differ. Equ. 2018, 2018, 1–17. [Google Scholar]
  15. Fedorov, V.E.; Abdrakhmanova, A.A. A class of initial value problems for distributed order equations with a bounded operator. In Stability, Control and Differential Games; Tarasyev, A., Maksimov, V., Filippova, T., Eds.; Springer Nature: Cham, Switzerland, 2020; pp. 251–262. [Google Scholar]
  16. Fedorov, V.E.; Abdrakhmanova, A.A. Distributed order equations in Banach spaces with sectorial operators. In Transmutation Operators and Applications; Kravchenko, V.F., Sitnik, S.M., Eds.; Springer Nature: Cham, Switzerland, 2020; pp. 509–538. [Google Scholar]
  17. Fedorov, V.E. On generation of an analytic in a sector resolving operators family for a distributed order equation. Zapiski POMI 2020, 489, 113–129. [Google Scholar]
  18. Fedorov, V.E. Generators of analytic resolving families for distributed order equations and perturbations. Mathematics 2020, 8, 1306. [Google Scholar] [CrossRef]
  19. Debbouche, A.; Fedorov, V.E. A class of fractional degenerate evolution equations with delay. Mathematics 2020, 8, 1700. [Google Scholar] [CrossRef]
  20. Gerasimov, A.N. A generalization of linear laws of deformation and its application to the problems of internal friction. Prikladnaya Matematika i Mekhanika 1948, 12, 251–260. (In Russian) [Google Scholar]
  21. Caputo, M.; Mainardi, F. A new dissipation model based on memory mechanism. Pure Appl. Geophys. 1971, 91, 134–147. [Google Scholar] [CrossRef]
  22. Novozhenova, O.G. Life and science of Alexey Gerasimov, one of the pioneers of fractional calculus in Soviet Union. Fract. Calc. Appl. Anal. 2017, 20, 790–809. [Google Scholar] [CrossRef]
  23. Bajlekova, E.G. Fractional Evolution Equations in Banach Spaces. Ph.D. Thesis, University Press Facilities, Eindhoven University of Technology, Eindhoven, The Netherlands, 2001. [Google Scholar]
  24. Hille, E.; Phillips, R.S. Functional Analysis and Semi-Groups; American Mathematical Society: Providence, RI, USA, 1957. [Google Scholar]
  25. Yosida, K. Functional Analysis; Springer: Berlin/Heidelberg, Germany, 1965. [Google Scholar]
  26. Kato, K. Perturbation Theory for Linear Operators; Springer: Berlin/Heidelberg, Germany, 1966. [Google Scholar]
  27. Da Prato, G.; Iannelli, M. Linear integro-differential equations in Banach spaces. Rend. Semin. Mat. della Univ. Padova 1980, 62, 207–219. [Google Scholar]
  28. Prüss, J. Evolutionary Integral Equations and Applications; Springer: Basel, Switzerland, 1993. [Google Scholar]
  29. Triebel, H. Interpolation Theory, Functional Spaces, Differential Operators; North Holland Publ.: Amsterdam, The Netherlands, 1978. [Google Scholar]
  30. Arendt, W.; Batty, C.J.K.; Hieber, M.; Neubrander, F. Vector-Valued Laplace Transforms and Cauchy Problems; Springer: Basel, Switzerland, 2011. [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

Fedorov, V.E.; Filin, N.V. On Strongly Continuous Resolving Families of Operators for Fractional Distributed Order Equations. Fractal Fract. 2021, 5, 20. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5010020

AMA Style

Fedorov VE, Filin NV. On Strongly Continuous Resolving Families of Operators for Fractional Distributed Order Equations. Fractal and Fractional. 2021; 5(1):20. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5010020

Chicago/Turabian Style

Fedorov, Vladimir E., and Nikolay V. Filin. 2021. "On Strongly Continuous Resolving Families of Operators for Fractional Distributed Order Equations" Fractal and Fractional 5, no. 1: 20. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5010020

Article Metrics

Back to TopTop