Next Article in Journal
Two Regularization Methods for the Variational Inequality Problem over the Set of Solutions of the Generalized Mixed Equilibrium Problem
Previous Article in Journal
Convergence of Special Sequences of Semi-Exponential Operators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Solvability of Equations with a Distributed Fractional Derivative Given by the Stieltjes Integral

by
Sergey M. Sitnik
1,*,
Vladimir E. Fedorov
2,3,
Nikolay V. Filin
2,4 and
Viktor A. Polunin
1
1
Applied Mathematics and Computer Modelling, Belgorod State National Research University, Pobedy St. 85, Belgorod 308015, Russia
2
Department of Differential Equations, N.N. Krasovskii Institute of Mathematics and Mechanics of the Ural Branch of the Russian Academy of Sciences, 16, S.Kovalevskaya St., Yekaterinburg 620108, Russia
3
Department of Mathematical Analysis, Faculty of Mathematics, Chelyabinsk State University, 129, Kashirin Brothers St., Chelyabinsk 454001, Russia
4
Department of Differential Equations, Yugra State University, 16, Chekhov St., Khanty-Mansiysk 628012, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(16), 2979; https://doi.org/10.3390/math10162979
Submission received: 26 July 2022 / Revised: 8 August 2022 / Accepted: 16 August 2022 / Published: 18 August 2022
(This article belongs to the Section Computational and Applied Mathematics)

Abstract

:
Linear equations in Banach spaces with a distributed fractional derivative given by the Stieltjes integral and with a closed operator A in the right-hand side are considered. Unlike the previously studied classes of equations with distributed derivatives, such kinds of equations may contain a continuous and a discrete part of the integral, i.e., a standard integral of the fractional derivative with respect to its order and a linear combination of fractional derivatives with different orders. Resolving families of operators for such equations are introduced into consideration, and their properties are studied. In terms of the resolvent of the operator A, necessary and sufficient conditions are obtained for the existence of analytic resolving families of the equation under consideration. A perturbation theorem for such a class of operators is proved, and the Cauchy problem for the inhomogeneous equation with a distributed fractional derivative is studied. Abstract results are applied for the research of the unique solvability of initial boundary value problems for partial differential equations with a distributed derivative with respect to time.

1. Introduction

Equations with fractional derivatives of various forms attract the attention of researchers both from a theoretical point of view and because of their widespread use in applied problems, see, e.g., recent papers [1,2,3] and many other works. The distributed derivatives (other names are continual derivatives [4], mean derivatives [5]) are used for the investigation of some real phenomena and processes when an order of a fractional derivative in a model continuously depends on the process parameters: in the theory of viscoelasticity [5], in modeling dielectric induction and diffusion [6,7], in the kinetic theory [8], and in other scientific fields [4,9,10,11,12]. This fact initiated the interest in equations with distributed derivatives of specialists in computational mathematics [13,14], of researchers in the qualitative theory of differential equations [15,16,17,18,19,20,21].
In the mentioned above works, researchers study specific equations or systems of them with some possible arbitrariness in the choice of parameters. The idea of the present work is to investigate the Cauchy problem for a class of abstract equations with distributed derivatives in order to be able to reduce many initial boundary value problems for partial differential equations or systems of equations of various forms to such a problem and study them through the obtained general results.
Let Z be a Banach space, D β be the fractional Gerasimov–Caputo derivative for β > 0 and the fractional Riemann–Liouville integral for β 0 , and A be a linear closed densely defined in Z operator. For an unknown function z : R + Z consider the distributed order equation
b c D α z ( t ) d μ ( α ) = A z ( t ) , t > 0 ,
with the Cauchy conditions
z ( k ) ( 0 ) = z k , k = 0 , 1 , , m 1 .
Here b , c R , b < c , m 1 < c m N , μ : [ b , c ] C is a function with a bounded variation, c is a variation point of the measure d μ ( t ) , the integral in Equation (1) is of Riemann–Stieltjes.
The Cauchy problem for the distributed order equation
b c ω ( α ) D α z ( t ) d α = A z ( t )
with a given scalar function ω : ( b , c ) C and a bounded operator A, or with an infinitesimal generator A of an analytic semigroup was studied in [22,23,24]. In the case b = 0 , c ( 0 , 1 ] , necessary and sufficient conditions on a linear closed operator A for the existence of an analytic resolving family of operators for Equation (3) were found in [25]. It allowed us to obtain a theorem on the existence of a unique solution for the corresponding inhomogeneous equation. In [26] these results were extended for the case c > 1 and a theorem on perturbations of generators of analytic resolving operators families for Equation (3) was proved. Paper [27] contains analogous results for equation with a discretely distributed Gerasimov–Caputo derivative
k = 1 n ω k D α k z ( t ) = A z ( t ) .
In works [28,29] initial value problems for equations in Banach spaces with distributed Riemann–Liouville derivatives were studied.
Equation (1) with the Riemann–Stieltjes integral in the present work includes the listed above classes of equations as partial cases. Indeed, an arbitrary function μ with a bounded variation has the form μ = μ c + μ d , where μ c is a continuous function with a bounded variation, μ d is a jumps function. Therefore, Equation (1) has the form
b c D α z ( t ) d μ ( α ) = b c D α z ( t ) d μ c ( α ) + b c D α z ( t ) d μ d ( α ) =
= b c μ c ( α ) D α z ( t ) d α + k = 1 n ω k D α k z ( t ) = A z ( t ) ,
if there exists an appropriate derivative μ c , α k are points of jumps of the function μ d , ω k are values of jumps, k = 1 , 2 , , n . Moreover, here we consider b 0 in the lower limit of integration and abandon the additional conditions of the fulfillment of some inequalities related to the integral from the equation (see [25,26]) since we prove that these inequalities follow from our general assumptions on parameters of the problem (see Lemma 1, Lemma 2).
The second section contains the study of the properties of some analytic functions, which are associated with an integral from Equation (1). Then the notion of a k-resolving family, k = 0 , 1 , , m 1 , is introduced and the properties of such families are researched. It is shown that the existence of a 0-resolving family implies the existence of k-resolving families, k = 1 , 2 , , m 1 . In the third section it is proved that inclusion A A W ( θ 0 , a 0 ) for some θ 0 ( π / 2 , π ) , a 0 0 is necessary and sufficient for the existence of analytic k-resolving families for distributed order Equation (1). This result allows us to obtain a unique solvability theorem for Cauchy problem Equations (1) and (2). In the fourth section a theorem on the perturbations for operators from the class A W ( θ 0 , a 0 ) is proved. The unique solvability of the inhomogeneous equation with the distributed derivative and with an operator A A W ( θ 0 , a 0 ) is studied in the fifth section. The last section concerns an application of abstract results to the investigation of an initial boundary value problem for the phase field system of equations with the distributed order time derivative.

2. Resolving Families of Operators and Their Properties

Let Z be a Banach space, denote for β > 0 , h : R + Z the Riemann–Liouville fractional derivative
J β h ( t ) : = 1 Γ ( β ) 0 t ( t s ) β 1 h ( s ) d s , t > 0 .
Let m 1 < α m N , D m be the derivative of the m-th order, then the Gerasimov–Caputo derivative has the form [30,31,32]
D α h ( t ) : = D m J m α h ( t ) k = 0 m 1 h ( k ) ( 0 ) t k k ! .
We will mean for α < 0 that D α h ( t ) : = J α h ( t ) .
The Laplace transform of a function h : R + Z will be denoted by h ^ or Lap [ h ] , if the expression for h is too large. The Laplace transform of the Gerasimov–Caputo derivative of an order α ( m 1 , m ] satisfies the equality (see, e. g., [33])
D α h ^ ( λ ) = λ α h ^ ( λ ) k = 0 m 1 h ( k ) ( 0 ) λ α 1 k .
Introduce the notations S θ , a : = { μ C : | arg ( μ a ) | < θ , μ a } for θ [ π / 2 , π ] , a R , Σ ψ : = { t C : | arg t | < ψ , t 0 } for ψ ( 0 , π / 2 ] .
Theorem 1
([34], Theorem 0.1, p. 5), ([35], Theorem 2.6.1, p. 84). Let θ 0 ( π / 2 , π ] , a R , X be a Banach space, a map H : ( a , ) X be set. The next assertions are equivalent.
(i)
There exists an analytic function F : Σ θ 0 π / 2 X , for every θ ( π / 2 , θ 0 ) there exists C ( θ ) > 0 , such that for all t Σ θ π / 2 the inequality F ( t ) X C ( θ ) e a Re t is satisfied; F ^ ( λ ) = H ( λ ) at λ > a .
(ii)
The map H is analytically continued on S θ 0 , a , for every θ ( π / 2 , θ 0 ) there exists such a K ( θ ) > 0 , that for all λ S θ , a   H ( λ ) X K ( θ ) | λ a | 1 .
Let by L ( Z ) the Banach space of all linear continuous operators from Z to Z be denoted, denote by C l ( Z ) the set of all linear closed operators, densely defined in Z , acting in the space Z . Endow the domain D A of an operator A C l ( Z ) by the norm of its graph · D A : = · Z + A · Z , then D A is a Banach space.
Consider the Cauchy problem
z ( k ) ( 0 ) = z k , k = 0 , 1 , , m 1 ,
for the distributed order equation
b c D α z ( t ) d μ ( α ) = A z ( t ) , t > 0 ,
where b , c R , b < c , m 1 < c m N , μ : [ b , c ] C is a function with a bounded variation, briefly μ B V ( [ b , c ] ; C ) , c is a variation point of the measure d μ ( t ) . The integral in Equation (6) is understood in the sense of Riemann–Stieltjes.
A solution of the problem in Equations (5) and (6) is a function z C m 1 ( R ¯ + ; Z ) C ( R + ; D A ) , such that b c ω ( α ) D α z ( t ) d α C ( R + ; Z ) and Equalities (5) and (6) are fulfilled. Here R ¯ + : = R + { 0 } .
It is evident, that under the conditions of this section the complex-valued functions
W ( λ ) : = b c λ α d μ ( α ) W k ( λ ) : = k c λ α d μ ( α ) , k = 0 , 1 , , m 1 ,
are analytic on the set S π , 0 . Here the Riemann–Stieltjes integrals are used also.
Lemma 1.
Let b , c R , b < c , m 1 < c m N , μ B V ( [ b , c ] ; C ) , c be a variation point of the measure d μ ( t ) . Then for k = 0 , 1 , , m 1
ε ( 0 , c ) C , ϱ > 0 λ S π , 0 \ { λ C : | λ | < ϱ } | W k ( λ ) | C | λ | c ε ,
ε ( 0 , c ) C , ϱ > 0 λ S π , 0 \ { λ C : | λ | < ϱ } | W ( λ ) | C | λ | c ε .
Proof. 
By the definition of the Riemann–Stieltjes integral for a small δ > 0 there exists a division { α 0 = b , α 1 , , α n 1 , α n = c } of the segment [ b , c ] with a sufficiently small radius r : = max { α k α k 1 : k = 1 , 2 , , n } and with intermediate points ξ k [ α k 1 , α k ] , k = 1 , 2 , , n , we obtain for | λ | > ϱ
b c λ α d μ ( α ) j = 1 n λ ξ j ( μ ( α j ) μ ( α j 1 ) ) δ
| λ | ξ n | μ ( α n ) μ ( α n 1 ) | j = 1 n 1 λ ξ j ξ n ( μ ( α j ) μ ( α j 1 ) ) δ | λ | ξ n
| μ ( α n ) μ ( α n 1 ) | 2 | λ | ξ n .
Here μ ( α n ) μ ( α n 1 ) = μ ( c ) μ ( α n 1 ) 0 for sufficiently small r, since c is a variation point of the measure d μ ( t ) .
For every ε ( 0 , c ) there exists a division of [ b , c ] , such that r ε , then α n α n 1 = c α n 1 ε and ξ n α n 1 c ε .
For W k the proof has the same form. □
Lemma 2.
Let b , c R , b < c , m 1 < c m N , μ B V ( [ b , c ] ; C ) , c be a variation point of the measure d μ ( t ) . Then for all k , l = 0 , 1 , , m 1 , k > l ,
C , ϱ > 0 λ S π , 0 \ { λ C : | λ | < ϱ } | W k ( λ ) W l ( λ ) | C | λ | k ;
C , ϱ > 0 λ S π , 0 \ { λ C : | λ | < ϱ } | W k ( λ ) W ( λ ) | C | λ | k .
Proof. 
Indeed,
| W k ( λ ) W l ( λ ) | = l k λ α d μ ( α ) V b c ( μ ) | λ | k ,
where V b c ( μ ) is the variation of the function μ on the segment [ b , c ] . The second inequality can be proved in the same way. □
Definition 1.
A family of operators { S k ( t ) L ( Z ) : t 0 } , k { 0 , 1 , , m 1 } , is called k-resolving for Equation (6), if:
(i)
S k ( t ) is strongly continuous for 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 of the Cauchy problem
z ( l ) ( 0 ) = 0 , l { 0 , 1 , , m 1 } \ { k } , z ( k ) ( 0 ) = z k
to Equation (6) for any z k D A .
Remark 1.
So, a k-resolving family { S k ( t ) L ( Z ) : t 0 } consists of operators, such that S ( t ) for every fixed t 0 maps any z k D A into the value z ( t ) = S k ( t ) z k at the point t of a solution of Cauchy problem Equations (6) and (7). Thus, the totality of families { S k ( t ) L ( Z ) : t 0 } , k = 0 , 1 , , m 1 , entirely determines the solution of the complete Cauchy problem Equations (5) and (6).
Denote by ρ ( A ) the resolvent set { λ C : ( λ I A ) 1 L ( Z ) } of an operator A C l ( Z ) .
Lemma 3.
Let b , c R , b < c , m 1 < c m N , μ B V ( [ b , c ] ; C ) , c be a variation point of the measure d μ ( t ) , for some k { 0 , 1 , , m 1 } there exist a k-resolving family of operators { S k ( t ) L ( Z ) : t 0 } for Equation (6), such that for all t 0 S k ( t ) L ( Z ) K e a t at some K 1 , a 0 . Then for Re λ > a we have W ( λ ) ρ ( A ) and
S ^ k ( λ ) = W k ( λ ) λ k + 1 ( W ( λ ) I A ) 1 .
Proof. 
For an exponentially bounded solution z of problem Equations (5) and (6), we have
b c D α z ( t ) d μ ( α ) = b 0 D α z ( t ) d μ ( α ) + l = 1 m 1 l 1 l D α z ( t ) d μ ( α ) + m 1 c D α z ( t ) d μ ( α ) ,
hence
Lap b c D α z ( t ) d μ ( α ) = b 0 λ α z ^ ( λ ) d μ ( α ) + l = 1 m 1 l 1 l λ α z ^ ( λ ) k = 0 l 1 z ( k ) ( 0 ) λ α k 1 d μ ( α ) +
+ m 1 c λ α z ^ ( λ ) k = 0 m 1 z ( k ) ( 0 ) λ α k 1 d μ ( α ) = W ( λ ) z ^ ( λ ) k = 0 m 1 z ( k ) ( 0 ) W k ( λ ) λ k + 1 = A z ^ ( λ )
due to the closedness of A. Therefore, points (ii) and (iii) of Definition 1 imply that for every z k D A
W ( λ ) S ^ k ( λ ) z k W k ( λ ) λ k + 1 z k = A S ^ k ( λ ) z k = S ^ k ( λ ) z k A z k .
Consequently, there exists an inverse operator for W ( λ ) I A : D A Z and Equality (8) holds. The right-hand side of this equation is a bounded operator by the assumptions of this lemma for S k , hence W ( λ ) ρ ( A ) . □
Theorem 2.
Let b , c R , b < c , m 1 < c m N , μ B V ( [ b , c ] ; C ) , c be a variation point of the measure d μ ( t ) , for some k { 0 , 1 , , m 1 } there exist a k-resolving family of operators { S k ( t ) L ( Z ) : t 0 } for Equation (6), such that for all t 0 S k ( k ) ( t ) L ( Z ) , S k ( k ) ( t ) L ( Z ) K e a t at some K 1 , a 0 . Then { S k ( k ) ( t ) L ( Z ) : t 0 } is continuous in the point t = 0 in the operator norm in L ( Z ) , if and only if A L ( Z ) .
Proof. 
For Re λ > a due to Lemma 3
0 e λ t ( S k ( k ) ( t ) I ) d t = W k ( λ ) λ ( W ( λ ) I A ) 1 I λ .
Let the function η ( t ) : = S k ( k ) ( t ) I L ( Z ) is continuous on the segment [ 0 , 1 ] and η ( 0 ) = 0 . For ε > 0 take δ > 0 , such that η ( t ) ε for all t [ 0 , δ ] , hence
W k ( λ ) λ ( 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 for t 0 . Therefore, for large enough Re λ > 0 W k ( λ ) ( W ( λ ) I A ) 1 I L ( Z ) < 1 , consequently, the operator ( W ( λ ) I A ) 1 is continuously invertible, [ ( W ( λ ) I A ) 1 ] 1 = W ( λ ) I A L ( Z ) . Thus, A L ( Z ) .
Let A L ( Z ) , R > max { ϱ , ( C 1 A L ( Z ) ) 1 / ( c ε ) } , where constants C , ϱ > 0 , ε ( 0 , c m + 1 ) are taken from Lemma 1. Construct the contour Γ R = Γ 1 , R Γ 2 , R Γ 3 , R , where Γ 1 , R = { R e i φ : φ ( π , π ) } , Γ 2 , R = { r e i π : r [ R , ) } , Γ 3 , R = { r e i π : r [ R , ) } . For t 0 due to Lemma 3 by the inverse Laplace transform we obtain
S k ( k ) ( t ) = 1 2 π i Γ R W k ( λ ) λ ( W ( λ ) I A ) 1 e λ t d λ = 1 2 π i Γ R W k ( λ ) λ l = 0 A l e λ t d λ W ( λ ) l + 1 =
= I + 1 2 π i Γ R 1 λ l = 1 A l e λ t d λ W ( λ ) l + 1 2 π i Γ R W k ( λ ) W ( λ ) λ l = 0 A l e λ t d λ W ( λ ) l + 1 .
The series is convergent, since for λ Γ R by choosing R | W ( λ ) | 1 A L ( Z ) < 1 due to Lemma 1, moreover, by Lemma 1 and Lemma 2 we obtain
A l λ W ( λ ) l L ( Z ) C 1 A L ( Z ) l | λ | ( c ε ) l + 1 , ( W k ( λ ) W ( λ ) ) A l λ W ( λ ) l + 1 L ( Z ) C 2 A L ( Z ) l | λ | ( c ε ) ( l + 1 ) + 1 k .
For small t > 0 take R = 1 / t and obtain
S k ( k ) ( t ) I L ( Z ) C 3 k = 1 3 l = 1 Γ k , R A L ( Z ) l | d λ | | λ | ( c ε ) l + 1 + l = 0 Γ k , R A L ( Z ) l | d λ | | λ | ( c ε ) ( l + 1 ) + 1 k
C 4 l = 1 A L ( Z ) l R ( c ε ) l + l = 0 A L ( Z ) l R ( c ε ) ( l + 1 ) k = C 4 t c ε A L ( Z ) 1 t c ε A L ( Z ) + C 4 t c ε k 1 t c ε A L ( Z ) 0
as t 0 + , since c ε k c ε m + 1 > 0 .  □
Theorem 3.
Let b , c R , b < c , 1 m 1 < c m , μ B V ( [ b , c ] ; C ) , c be a variation point of the measure d μ ( t ) , there exist a 0-resolving family of operators { S 0 ( t ) L ( Z ) : t 0 } for Equation (6), such that for all t 0 S 0 ( 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 , for Equation (6).
Proof. 
Consider for k = 1 , 2 , , m 1 functions
W k ( λ ) λ k W 0 ( λ ) = 1 λ k W 0 ( λ ) W k ( λ ) λ k W 0 ( λ ) .
Due to Lemmas 1 and 2
W 0 ( λ ) W k ( λ ) λ k W 0 ( λ ) C 1 | λ | c ε .
For c > 1 choose ε > 0 , such that c ε > 1 , then there exists the inverse Laplace transform
w k ( t ) : = d i d + i W k ( λ ) e λ t λ k W 0 ( λ ) d λ , d > 0 , k = 1 , 2 , , m 1 ,
moreover, there exists K 1 > 0 , such that for all t 0 | w k ( t ) | K 1 e d t . For k = 1 we take into account that λ 1 = 1 ^ . The functions w k are continuous for t 0 , since integrals converge uniformly with respect to t on every segment. Consequently, functions of the form S k ( t ) : = 0 t w k ( t s ) S 0 ( s ) d s , k = 1 , 2 , , m 1 , are strongly continuous for t 0 ; point (ii) of Definition 1 holds for S k also. Moreover, S k ( t ) L ( Z ) K 2 e ( a + d ) t for all t 0 .
The derivatives
w k ( l ) ( t ) : = d i d + i W k ( λ ) e λ t λ k l W 0 ( λ ) d λ , d > 0 , l = 1 , 2 , , k 1 ; k = 1 , 2 , , m 1 ,
are continuous for t 0 also, since
W k ( λ ) λ k l W 0 ( λ ) = 1 λ k l W 0 ( λ ) W k ( λ ) λ k l W 0 ( λ ) , W 0 ( λ ) W k ( λ ) λ k l W 0 ( λ ) C 1 | λ | c l ε .
From relations Equation (9) it follows 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 k D A
S k ( l ) ( t ) z k = 0 t w k ( l ) ( t s ) S 0 ( s ) z k d s , l = 1 , 2 , , k 1 ,
S k ( l ) ( t ) z k = w k ( k 1 ) ( t ) S 0 ( l k ) ( 0 ) z k + 0 t w k ( k 1 ) ( s ) S 0 ( l k + 1 ) ( t s ) z k d s , l = k , k + 1 , , m 1 ,
therefore, S k ( t ) z k satisfies initial value conditions Equation (7).
From point (iii) of Definition 1, it follows that
W ( λ ) S ^ 0 ( λ ) z k λ 1 W 0 ( λ ) z k = A S ^ 0 ( λ ) z k , z k D A .
Since
S ^ k ( λ ) = w ^ k ( λ ) S ^ 0 ( λ ) = W k ( λ ) λ k W 0 ( λ ) S ^ 0 ( λ ) ,
after the multiplying Equation (10) by w ^ k ( λ ) we obtain
W k ( λ ) λ k W 0 ( λ ) W ( λ ) S ^ 0 ( λ ) z k W k ( λ ) λ k + 1 z k = W k ( λ ) λ k W 0 ( λ ) A S ^ 0 ( λ ) z k ,
W ( λ ) S ^ k ( λ ) z k W k ( λ ) λ k + 1 z k = A S ^ k ( λ ) z k .
Acting by the inverse Laplace transform, we obtain that S k ( t ) z k is a solution of Equation (6). □

3. Analytic Resolving Families

A resolving family of operators is called analytic, if it has an analytic continuation to a sector Σ ψ 0 at some ψ 0 ( 0 , π / 2 ] . An analytic resolving family of operators { S ( t ) L ( Z ) : t 0 } has a type ( ψ 0 , a 0 ) at some ψ 0 ( 0 , π / 2 ] , a 0 R , if for all ψ ( 0 , ψ 0 ) , a > a 0 there exists C ( ψ , a ) > 0 , such that for all t Σ ψ the inequality S ( t ) L ( Z ) C ( ψ , a ) e a Re t is satisfied.
Remark 2.
Analogous notions of analytic resolving families of operators are used in the study of integral evolution equations [34] and fractional differential equations [36]. They generalize the notion of analytic resolving semigroup of operators for the first order equation D t 1 z ( t ) = A z ( t ) (see [37,38,39]).
Following the works [25,26] define a class A W ( θ 0 , a 0 ) as the set of all operators A C l ( Z ) satisfying the following conditions:
(i)
there exist θ 0 ( π / 2 , π ] , a 0 0 , such that W ( λ ) ρ ( A ) for every λ S θ 0 , a 0 ;
(ii)
for every θ ( π / 2 , θ 0 ) , a > a 0 there exists K ( θ , a ) > 0 , such that for all λ S θ , a
( W ( λ ) I A ) 1 L ( Z ) | λ | K ( θ , a ) | W ( λ ) | | λ a | .
Remark 3.
The classes A W ( θ 0 , a 0 ) in works [25,26,27] are partial cases of this class with the same denotation A W ( θ 0 , a 0 ) due to the more general construction of the distributed derivative in the present work. If μ is a constant, excluding a unique jump in the point α = c , class A W ( θ 0 , a 0 ) coincides with the class A α ( θ 0 , a 0 ) , defined in [36]. Operators from the class A 1 ( θ 0 , a 0 ) are generators of an analytic semigroup of operators exactly [37,38,39].
Remark 4.
If A L ( Z ) , then for
| λ | > ( 2 A L ( Z ) ) 1 / ( c ε )
due to Lemma 1
( W ( λ ) I A ) 1 = n = 0 A n ( W ( λ ) ) n + 1 ,
( W ( λ ) I A ) 1 L ( Z ) 1 | W ( λ ) | n = 0 A L ( Z ) n | λ | ( c ε ) n 2 | λ | | W ( λ ) | | λ a | 1 + a | λ | .
We can choose a 0 = 2 ( 2 A L ( Z ) ) 1 / ( c ε ) , θ 0 = π / 4 , then for λ S θ 0 , a 0 we have inequalities Equation (11) and 1 + a 0 | λ | 1 + 2 . Thus, A A W ( θ 0 , a 0 ) . In this reasoning we may take a 0 greater and θ 0 closer to π, if necessary.
For A A W ( θ 0 , a 0 ) the operators
Z k ( t ) : = 1 2 π i Γ W k ( λ ) λ k + 1 W ( λ ) I A 1 e λ t d λ , k = 0 , 1 , , m 1 ,
are defined at t > 0 . Here Γ = Γ + Γ Γ 0 , Γ ± = { μ C : μ = a + r e ± i θ , r ( δ , ) } , Γ 0 = { μ C : μ = a + δ e i φ , φ ( θ , θ ) } for some δ > 0 , a > a 0 , θ ( π / 2 , θ 0 ) .
Theorem 4.
Let b , c R , b < c , m 1 < c m N , μ B V ( [ b , c ] ; C ) , c be a variation point of the measure d μ ( t ) . Then there exists an analytic 0-resolving family of operators of the type ( θ 0 π / 2 , a 0 ) for Equation (6), if and only if A A W ( θ 0 , a 0 ) . In this case, there exists a unique k-resolving family of operators for every k = 0 , 1 , , m 1 , it has the form { Z k ( t ) L ( Z ) : t 0 } .
Proof. 
Let A A W ( θ 0 , a 0 ) , R > δ , Γ 1 , R = Γ 0 , Γ 2 , R = { λ C : λ = a + R e i φ , φ ( θ , θ ) } ,   Γ 3 , R = { λ C : λ = a + r e i θ , r [ δ , R ] } , Γ 4 , R = { λ C : λ = a + r e i θ , r [ δ , R ] } ,   Γ R = k = 1 4 Γ k , R is the positively oriented closed loop, Γ 5 , R = { λ C : λ = a + r e i θ , r [ R , ) } , Γ 6 , R = { λ C : λ = a + r e i θ , r [ R , ) } , then we have Γ = Γ 5 , R Γ 6 , R Γ R \ Γ 2 , R .
For t > 0 , z 0 D A
Z 0 ( t ) z 0 = 1 2 π i Γ W 0 ( λ ) λ W ( λ ) I A 1 z 0 e λ t d λ = 1 2 π i Γ e λ t λ z 0 d λ +
+ 1 2 π i Γ e λ t λ W ( λ ) I A 1 A z 0 d λ + 1 2 π i Γ W 0 ( λ ) W ( λ ) λ W ( λ ) I A 1 e λ t z 0 d λ .
For t [ 0 , 1 ] , λ Γ \ { λ C : | λ | ϱ }
e λ t λ W ( λ ) I A 1 A z 0 Z e a + δ K ( θ , a ) A z 0 Z | W ( λ ) | | λ a | C 1 | λ | c + 1 ε ,
W 0 ( λ ) W ( λ ) λ W ( λ ) I A 1 e λ t z 0 Z C e a + δ K ( θ , a ) z 0 Z | W ( λ ) | | λ a | C 1 | λ | c + 1 ε ,
where we can take any ε ( 0 , c ) , see Lemma 1. Since c + 1 ε > 1 , the integral Z 0 ( t ) converges uniformly with respect to t [ 0 , 1 ] and by the continuity
Z 0 ( 0 ) z 0 = z 0 + 1 2 π i Γ 1 λ W ( λ ) I A 1 A z 0 d λ +
+ 1 2 π i Γ W 0 ( λ ) W ( λ ) λ W ( λ ) I A 1 z 0 d λ =
= z 0 + lim R 1 2 π i Γ R Γ 2 , R + Γ 5 , R + Γ 6 , R 1 λ W ( λ ) I A 1 A z 0 d λ +
+ lim R 1 2 π i Γ R Γ 2 , R + Γ 5 , R + Γ 6 , R W 0 ( λ ) W ( λ ) λ W ( λ ) I A 1 z 0 d λ = z 0
due to Cauchy theorem and inequalities
Γ s , R 1 λ W ( λ ) I A 1 A z 0 d λ Z C 2 R c ε , s = 2 , 5 , 6 ,
Γ s , R W 0 ( λ ) W ( λ ) λ W ( λ ) I A 1 z 0 d λ Z C 2 R c ε , s = 2 , 5 , 6 .
Analogously, for t > 0 , z 0 D A , k { 1 , 2 , , m 1 }
Z 0 ( k ) ( t ) z 0 = 1 2 π i Γ λ k 1 W 0 ( λ ) W ( λ ) I A 1 e λ t z 0 d λ =
= 1 2 π i Γ λ k 1 e λ t d λ z 0 + 1 2 π i Γ λ k 1 W ( λ ) I A 1 A e λ t z 0 d λ +
+ 1 2 π i Γ λ k 1 ( W 0 ( λ ) W ( λ ) ) W ( λ ) I A 1 e λ t z 0 d λ ,
for t [ 0 , 1 ] , λ Γ \ { λ C : | λ | ϱ }
λ k 1 W ( λ ) I A 1 e λ t A z 0 Z | λ | k e a + δ K ( θ , a ) A z 0 Z | W ( λ ) | | λ a | C 1 | λ | c + 1 ε k ,
λ k 1 ( W 0 ( λ ) W ( λ ) ) W ( λ ) I A 1 e λ t z 0 Z C | λ | k e a + δ K ( θ , a ) z 0 Z | W ( λ ) | | λ a | C 1 | λ | c + 1 ε k ,
Z 0 ( k ) ( 0 ) z 0 = 0 + 1 2 π i Γ λ k 1 W ( λ ) I A 1 A z 0 d λ +
+ 1 2 π i Γ λ k 1 ( W 0 ( λ ) W ( λ ) ) W ( λ ) I A 1 z 0 d λ = 0 ,
since c + 1 ε k > 1 for ε ( 0 , c m + 1 ) .
Consider for t > 0 , z k D A , k = 1 , 2 , , m 1 , l = 0 , 1 , , m 1
Z k ( l ) ( t ) z k = 1 2 π i Γ λ l k 1 W k ( λ ) W ( λ ) I A 1 e λ t z k d λ =
= 1 2 π i Γ λ l k 1 e λ t d λ z k + 1 2 π i Γ λ l k 1 W ( λ ) I A 1 e λ t A z k d λ +
+ 1 2 π i Γ λ l k 1 ( W k ( λ ) W ( λ ) ) W ( λ ) I A 1 e λ t z k d λ ,
for t [ 0 , 1 ] , λ Γ \ { λ C : | λ | ϱ }
λ l k 1 W ( λ ) I A 1 e λ t A z k Z | λ | l k e a + δ K ( θ , a ) A z k Z | W ( λ ) | | λ a | C 1 | λ | c + 1 ε l + k ,
λ l k 1 ( W k ( λ ) W ( λ ) ) W ( λ ) I A 1 e λ t z k Z C 2 | λ | l e a + δ K ( θ , a ) | W ( λ ) | | λ a | C 1 | λ | c + 1 ε l .
Then for k l   Z k ( l ) ( 0 ) z k = 0 , Z k ( k ) ( 0 ) z k = z k .
Thus, the functions Z k ( t ) z k C m 1 ( R ¯ + ; Z ) for z k D A satisfy Cauchy conditions Equation (7) with the correponding k { 0 , 1 , , m 1 } . Since the operator A is closed and commutes with ( W ( λ ) I A ) 1 on D A , at z k D A A Z k ( t ) z k C ( R + ; Z ) , so, Z k ( t ) z k C ( R + ; D A ) , k = 0 , 1 , , m 1 .
For Re ν > a we have
Z ^ k ( ν ) = 1 2 π i Γ W k ( λ ) λ k + 1 ( ν λ ) W ( λ ) I A 1 d λ .
For λ S θ , a
W k ( λ ) λ k + 1 ( ν λ ) W ( λ ) I A 1 L ( Z ) K ( θ , a ) | λ | k | λ a | | λ ν | ,
hence,
lim R 1 2 π i Γ s , R W k ( λ ) λ k + 1 ( ν λ ) W ( λ ) I A 1 d λ = 0 , s = 2 , 5 , 6 ,
and by the Cauchy integral formula
Z ^ k ( ν ) = lim R 1 2 π i Γ R W k ( λ ) λ k + 1 ( ν λ ) W ( λ ) I A 1 d λ = W k ( ν ) ν k + 1 W ( ν ) I A 1 .
Take in Theorem 1
H k ( λ ) = W k ( λ ) λ k + 1 ( W ( λ ) I A ) 1 , F k = Z k , k = 0 , 1 , , m 1 ,
then due to the inclusion A A W ( θ 0 , a 0 ) and Lemma 1 for every θ ( π / 2 , θ 0 ) , a > a 0
H k ( λ ) L ( Z ) = 1 | λ | k + 1 ( W k ( λ ) W ( λ ) ) ( W ( λ ) I A ) 1 L ( Z ) + + W ( λ ) ( W ( λ ) I A ) 1 L ( Z ) K ( θ , a ) | λ a | 1 | W ( λ ) | + 1 K 1 ( θ , a ) | λ a |
and by Theorem 1 the mappings Z k : Σ θ 0 π / 2 L ( Z ) are analytic and for every θ ( π / 2 , θ 0 ) , a > a 0 there exists C k ( θ , a ) > 0 such that for all t Σ θ π / 2 Z k ( t ) L ( Z ) C k ( θ , a ) e a Re t , k = 0 , 1 , , m 1 .
For z k D A put x k ( t ) : = Z k ( t ) z k , then A x ^ k ( λ ) = λ k 1 W k ( λ ) W ( λ ) I A 1 A z k . Hence x ^ k ( λ ) D A , A x ^ k ( λ ) = A x ^ k ( λ ) , x ^ k ( λ ) and A x ^ k ( λ ) have analytic extensions on S θ 0 , a 0 , since A A W ( θ 0 , a 0 ) . By Formula (4) of the Laplace transform
Lap b c D α x k ( t ) d μ ( α ) ( λ ) = W ( λ ) W k ( λ ) λ k + 1 W ( λ ) I A 1 z k W k ( λ ) λ k + 1 z k = = W k ( λ ) λ k + 1 W ( λ ) I A 1 A z k = A x ^ k ( λ ) .
Apply the inverse Laplace transform to both sides of the obtained equality and get equality Equation (6) for all points of the function A x k continuity, hence for all t > 0 . Therefore, x k is a solution of problem Equations (6) and (7) and { Z k ( t ) L ( Z ) : t 0 } is an analytic k-resolving family of operators of the type ( θ 0 π / 2 , a 0 ) for Equation (6), k = 0 , 1 , , m 1 .
If there exists an analytic 0-resolving family of operators { S 0 ( t ) L ( Z ) : t 0 } of the type ( θ 0 π / 2 , a 0 ) for Equation (6), by Lemma 3
S ^ 0 ( λ ) = W 0 ( λ ) λ ( W ( λ ) I A ) 1 , Re λ > a 0 .
Theorem 1 implies that A A W ( θ 0 , a 0 ) , S 0 ( t ) Z 0 ( t ) by virtue of the uniqueness of the inverse Laplace transform.
By Theorem 3 there exist k-resolving families of operators { S k ( t ) L ( Z ) : t 0 } for Equation (6), k = 1 , 2 , , m 1 , such that
S ^ k ( λ ) = W k ( λ ) λ k + 1 ( W ( λ ) I A ) 1 , Re λ > a 0 .
This equality implies that S k ( t ) Z k ( t ) , k = 1 , 2 , , m 1 . □
Theorem 5.
Let b , c R , b < c , m 1 < c m N , μ B V ( [ b , c ] ; C ) , c be a variation point of the measure d μ ( t ) , A A W ( θ 0 , a 0 ) . Then for any z k D A , k = 0 , 1 , , m 1 , the function z ( t ) = k = 0 m 1 Z k ( t ) z k is a unique solution of problems (5) and (6). In this case, the solution is analytic in the sector Σ θ 0 π / 2 .
Proof. 
From Theorem 4 and linearity of problems (5) and (6) it follows that this funtion z is an analytic in the sector Σ θ 0 π / 2 solution of the problem.
If there exist two solutions y 1 , y 2 of problems (5) and (6), then their difference y = y 1 y 2 is a solution of Equation (6), which satisfies the initial value conditions y ( k ) ( 0 ) = 0 , k = 0 , 1 , , m 1 . Take T > 0 and redefine the function y on [ T , ) by zero, denote the obtained function by y T . It satisfies Equation (6) on R + , excluding, possibly, the point t = T , where the function y T may be discontinuous. Acting by the Laplace transform on both parts of Equation (6) and due to the initial conditions, we get the equality W ( λ ) y ^ T ( λ ) = A y ^ T ( λ ) . Since A A W ( θ 0 , a 0 ) , for λ S θ 0 , a 0 we get the identity y ^ T ( λ ) = 0 . This means that y T 0 and y 1 ( t ) = y 2 ( t ) for t [ 0 , T ) . We can choose arbitrary large T > 0 , therefore, y 1 y 2 on R ¯ + . □
Theorem 6.
Let b , c R , b < c , 2 < c , μ B V ( [ b , c ] ; C ) , c be a variation point of the measure d μ ( t ) , μ ( t ) R for all t from some left neighbourhood of c, A A W ( θ 0 , a 0 ) for some θ 0 ( π / 2 , π ) , a 0 0 . Then A L ( Z ) .
Proof. 
Due to the definition of the Riemann–Stieltjes integral for a small δ 1 > 0 there exists a division { α 0 = b , α 1 , , α n 1 , α n = c } of the segment [ b , c ] with a sufficiently small radius r : = max { α k α k 1 : k = 1 , 2 , , n } and with any intermediate points ξ k [ α k 1 , α k ] , k = 1 , 2 , , n , such that
W ( λ ) j = 1 n λ ξ j ( μ ( α j ) μ ( α j 1 ) ) < δ 1 ,
for all sufficiently large | λ | due to Lemma 1 | W ( λ ) | > 1 , therefore, for some d > 0
arg W ( λ ) arg j = 1 n λ ξ j ( μ ( α j ) μ ( α j 1 ) ) < d δ 1 .
Since lim | λ | j = 1 n 1 λ ξ j ξ n ( μ ( α j ) μ ( α j 1 ) ) = 0 and μ ( α n ) μ ( α n 1 ) R for a sufficiently small radius r > 0 , for every γ > 0 there exists δ > 0 , such that for all | λ | > δ
ξ n arg λ γ 2 arg j = 1 n λ ξ j ( μ ( α j ) μ ( α j 1 ) ) = = arg λ ξ n + arg j = 1 n 1 λ ξ j ξ n ( μ ( α j ) μ ( α j 1 ) ) + μ ( α n ) μ ( α n 1 ) ξ n arg λ + γ 2 .
Choosing ξ n = c and sufficiently small δ 1 < γ / ( 2 d ) we obtain
| arg W ( λ ) c arg λ | < γ .
Since c > 2 , there exist λ S θ 0 , a 0 , such that | arg W ( λ ) | > π . We have also
C | λ | c ε W ( λ ) C 1 | λ | c .
Further, without limitation of generality, we can assume that C 1 1 .
For arbitrary ν 0 C , such that | ν 0 | > ( C 1 δ ) c + 1 , take λ 0 = ν 0 1 / c , then | λ 0 | > δ , arg λ 0 = arg ν 0 / c ( π / c , π / c ) , λ 0 S θ 0 , a 0 for sufficiently large δ > 0 . The boundary of the region
Ω λ 0 : = λ C : | ν 0 | 1 C 1 / ( c ε ) < | λ | < | ν 0 | + 1 C 1 1 / c , arg λ arg ν 0 c < 2 γ c ,
which belongs to S θ 0 , a 0 for small enough γ > 0 , is mapped by the function ν = W ( λ ) into the contour, for the point of which due to inequalities (12) and (13) | | ν | | ν 0 | | 1 , | arg ν arg ν 0 | γ . Therefore, ν 0 lies inside the contour W [ Ω λ 0 ] and is the image of some point from Ω λ 0 . Thus, { ν C : | ν | > ( C 1 δ ) c + 1 } W [ S θ 0 , a 0 ] ρ ( A ) , since A A W ( θ 0 , a 0 ) . Moreover, for sufficiently large | ν | , where ν = W ( λ ) ,
ν R ν ( A ) L ( Z ) K ( θ , a ) | λ | | λ a | C 2 ,
hence Lemma 5.2 [40] implies the boundedness of the operator A. □

4. Perturbations of Operators of the Class A W ( θ 0 , a 0 )

The result of this section generalizes the perturbation theorem for analytic semigroups of operators [39] and a similar result for generators of resolving families of the distributed fractional derivative of the partial form in [26].
Theorem 7.
Let b , c R , b < c , m 1 < c m N , μ B V ( [ b , c ] ; C ) , c be a variation point of the measure d μ ( t ) , A A W ( θ 0 , a 0 ) for some θ 0 ( π / 2 , π ) , a 0 0 , B C l ( Z ) , for all x D A D B
B x Z β A x Z + γ x Z ,
where β , γ 0 , there exists q ( 0 , 1 ) such that β ( 1 + K ( θ , a ) ) < q for every θ ( π / 2 , θ 0 ) , a > a 0 . Then A + B A W ( θ 0 , a 1 ) for sufficiently large a 1 > a 0 .
Proof. 
Take k > sin 1 θ 0 , λ S θ , k a S θ , a for some θ ( π / 2 , θ 0 ) , a > a 0 , then inequality (14) implies that
B ( W ( λ ) I A ) 1 L ( Z ) β A ( W ( λ ) I A ) 1 L ( Z ) + γ ( W ( λ ) I A ) 1 L ( Z ) β 1 + | λ | K A ( θ , a ) | λ a | + γ | λ | K A ( θ , a ) | λ a | | W ( λ ) | ,
where K A is the constant from the definition of the class A W ( θ 0 , a 0 ) . Note that the value
| λ | | λ a | 1 1 a | λ | 1 1 1 k sin θ 0
is close to 1 and
| λ | | λ a | | W ( λ ) | 1 1 1 k sin θ 0 C ( k a 0 sin θ 0 ) c ε
is close to 0 for a large enough number k. Here C is the constant from Lemma 1. Hence for such a k the inequality
B ( W ( λ ) I A ) 1 L ( Z ) β 1 + K A ( θ , a ) 1 1 k sin θ 0 + γ K A ( θ , a ) 1 1 k sin θ 0 C ( k a 0 sin θ 0 ) c ε q < 1
holds. Further, we have
( W ( λ ) I A B ) 1 ( W ( λ ) I A ) 1 ( I B ( W ( λ ) I A ) 1 ) 1 = = ( W ( λ ) I A ) 1 n = 0 [ B ( W ( λ ) I A ) 1 ] n ,
| λ k a | | λ a | = 1 ( k 1 ) a λ a 1 + ( k 1 ) a | λ a | 1 + 1 sin θ 0 .
Therefore,
( W ( λ ) I A B ) 1 L ( Z ) | λ | K A ( θ , a ) ( 1 q ) | λ a | | W ( λ ) | | λ | K A ( θ , a ) 1 + 1 sin θ 0 ( 1 q ) | λ k a | | W ( λ ) | .
Thus, A + B A W ( θ 0 , a 1 ) , a 1 = k a 0 , we can take for all θ ( π / 2 , θ 0 ) , a > a 1
K A + B ( θ , a ) = K A ( θ , a / k ) 1 q 1 + 1 sin θ 0 .
Remark 5.
Any bounded operator B L ( Z ) satisfies (14) with β = 0 , γ = B L ( Z ) .

5. Inhomogeneous Equation

A solution to the Cauchy problem
z ( k ) ( 0 ) = z k , k = 0 , 1 , , m 1 ,
for the inhomogeneous equation
b c D α z ( t ) d μ ( α ) = A z ( t ) + g ( t ) , t ( 0 , T ] ,
where b < c , m 1 < c m N , μ B V ( [ b , c ] ; C ) , T > 0 , g C ( [ 0 , T ] ; Z ) , is a function z C m 1 ( [ 0 , T ] ; Z ) C ( ( 0 , T ] ; D A ) , such that there exists b c D α z ( t ) d μ ( α ) C ( ( 0 , T ] ; Z ) and equalities (15) and (16) are fulfilled.
Denote
Z ( t ) : = 1 2 π i Γ e λ t W ( λ ) I A 1 d λ .
Lemma 4.
Let b , c R , b < c , m 1 < c m N , μ B V ( [ b , c ] ; C ) , c be a variation point of the measure d μ ( t ) , θ 0 ( π / 2 , π ] , a 0 0 , A A W ( θ 0 , a 0 ) , g C ( [ 0 , T ] ; D A ) . Then the function
z g ( t ) = 0 t Z ( t s ) g ( s ) d s
is a unique solution of problems (15) and (16) with z k = 0 , k = 0 , 1 , , m 1 .
Proof. 
Due to Lemma 1 and the inclusion A A W ( θ 0 , a 0 ) , we have for ε ( 0 , c m + 1 )
Z ( k ) ( t ) L ( Z ) C 1 Γ e t Re λ | d λ | | λ | c ε k , k = 0 , 1 , , m 1 .
Since c ε k c ε m + 2 > 1 for k = 0 , 1 , , m 2 , we have Z ( k ) ( 0 ) = 0 . For k = m 1
Γ 0 e t Re λ | d λ | | λ | c ε m + 1 2 π δ m c + ε e a + δ ,
Γ ± e t Re λ | d λ | | λ | c ε m + 1 δ e r t cos θ d r r c ε m + 1 = ( t cos θ ) c ε m Γ ( m c + ε ) .
Thus, Z ( m 1 ) ( t ) L ( Z ) = O ( t c ε m ) as t 0 + ,
z g ( k ) ( t ) = 0 t Z ( k ) ( t s ) g ( s ) d s ,
z g ( k ) ( 0 ) = 0 , k = 0 , 1 , , m 2 , z g ( m 1 ) ( t ) C 1 t c ε m + 1 0 as t 0 + . Therefore, conditions (15) with z k = 0 , k = 0 , 1 , , m 1 , are fulfilled.
Define g ( t ) 0 at t > T , then z g is defined on R + , z g = Z * g , z ^ g = Z ^ g ^ . Besides, Z ^ ( λ ) = W ( λ ) I A 1 , since due Lemma 1 and the inclusion A A W ( θ 0 , a 0 ) at λ Γ \ { κ C : | κ | < ϱ } for Re ν > a
1 ν λ W ( λ ) I A 1 C 3 | λ | c ε + 1 .
Since c ε + 1 > 1 ,
Lap 0 b D α z g d μ ( α ) ( ν ) = W ( ν ) W ( ν ) I A 1 g ^ ( ν ) = g ^ ( ν ) + A W ( ν ) I A 1 g ^ ( ν ) .
Acting by the inverse Laplace transform, we get
0 b D α z g ( t ) d μ ( α ) = g ( t ) + A ( Z * g ) ( t ) = g ( t ) + A z g ( t ) .
Here due to the inclusion g C ( [ 0 , T ] ; D A ) , the closedness of the operator A and the commutation of operators Z ( t ) , t R + , with A by their construction, the values A ( Z * g ) ( t ) = Z * A g ( t ) , t ( 0 , T ] , are defined.
The proof of the uniqueness of problems (15) and (16) solution is the same as in the proof of Theorem 5. □
Denote by C γ ( [ 0 , T ] ; Z ) for γ ( 0 , 1 ] the class of functions f : [ 0 , T ] Z , such that for all t , s [ 0 , T ] the inequality f ( t ) f ( s ) Z C | t s | γ is satisfied with some C > 0 .
Lemma 5.
Let b , c R , b < c , m 1 < c m N , μ B V ( [ b , c ] ; C ) , c be a variation point of the measure d μ ( t ) , θ 0 ( π / 2 , π ] , a 0 0 , A A W ( θ 0 , a 0 ) , γ ( 0 , 1 ] , g C γ ( [ 0 , T ] ; Z ) . Then function (17) is a unique solution to problems (15) and (16) with z k = 0 , k = 0 , 1 , , m 1 .
Proof. 
Denote for t > 0
Y ( t ) = 1 2 π i Γ ( W ( λ ) W 0 ( λ ) ) ( W ( λ ) I A ) 1 e λ t d λ .
Due to the closedness of A we have for t > 0
A Z ( t ) = 1 2 π i Γ A ( W ( λ ) I A ) 1 e λ t d λ = 1 2 π i Γ W ( λ ) ( W ( λ ) I A ) 1 e λ t d λ = = Y ( t ) + D 1 Z 0 ( t ) , Y ( t ) L ( Z ) C 1 Γ e t Re λ | d λ | | λ | c ε ,
hence im Z ( t ) D A . Arguing as in the previous proof, we obtain that Y ( t ) L ( Z ) = O ( t c ε 1 ) as t 0 + .
Since
Γ W ( λ ) ( W ( λ ) I A ) 1 e λ t d λ L ( Z ) C 2 + K ( θ , a ) Γ ± | λ | e t Re λ | λ a | d s C 3 t cos θ ,
we have A Z ( t ) L ( Z ) = O ( t 1 ) as t 0 + , therefore,
A Z ( t s ) ( g ( s ) g ( t ) ) Z C 3 | t s | γ 1 .
Then the integral
0 t A Z ( t s ) g ( s ) d s = 0 t A Z ( t s ) ( g ( s ) g ( t ) ) d s + 0 t Y ( t s ) g ( t ) d s + ( Z 0 ( t ) I ) g ( t )
converges, since
0 t A Z ( t s ) g ( s ) d s Z C 4 ( t γ + t c ε ) + ( Z 0 ( t ) I ) g ( t ) Z .
Therefore, z g ( t ) D A , z g C ( [ 0 , T ] ; D A ) .
The rest of the proof is the same as for Lemma 4. □
Theorem 5, Lemmas 4 and 5 imply the following assertion.
Theorem 8.
Let b , c R , b < c , m 1 < c m N , μ B V ( [ b , c ] ; C ) , c be a variation point of the measure d μ ( t ) , θ 0 ( π / 2 , π ] , a 0 0 , A A W ( θ 0 , a 0 ) , g C ( [ 0 , T ] ; D A ) C γ ( [ 0 , T ] ; Z ) , γ ( 0 , 1 ] , z k D A , k = 0 , 1 , , m 1 . Then the function
z ( t ) = k = 0 m 1 Z k ( t ) z k + 0 t Z ( t s ) g ( s ) d s
is a unique solution of problems (15) and (16).

6. Application to an Initial-Boundary Value Problem

Let Ω R d be a bounded region with a boundary Ω of the class C , β , γ , δ , ν R , b 0 , c ( 1 , 2 ) , α 1 < α 2 < < α n c , ω j R \ { 0 } , j = 1 , 2 , , n , ω C ( [ b , c ] ; R ) , ω ( c ) 0 in a some left vicinity of c. Consider the initial-boundary value problem
u ( s , 0 ) = u 0 ( s ) , v ( s , 0 ) = v 0 ( s ) , s Ω ,
u t ( s , 0 ) = u 1 ( s ) , v t ( s , 0 ) = v 1 ( s ) , s Ω ,
u ( s , t ) = v ( s , t ) = 0 , ( s , t ) Ω × ( 0 , T ] ,
for the system of equations in Ω × ( 0 , T ]
j = 1 n ω j D t α j u ( s , t ) + b c ω ( α ) D t α u ( s , t ) d α = Δ u ( s , t ) Δ v ( s , t ) + f 1 ( s , t ) ,
j = 1 n ω j D t α j v ( s , t ) + b c ω ( α ) D t α v ( s , t ) d α = ν Δ v ( s , t ) + β u ( s , t ) + γ v ( s , t ) + f 2 ( s , t ) .
The system at ω 2 = ω 3 = = ω n = 0 , α 1 = 1 , ω ( α ) 0 for all α ( b , c ) up to linear replacement of unknown functions u ( s , t ) = u ˜ ( s , t ) + l 2 v ˜ ( s , t ) , v ( s , t ) = l 2 v ˜ ( s , t ) , l R , coincides with the linearization of the phase field system of equations, describing phase transitions of the first kind within the framework of mesoscopic theory [41,42].
Set
Z = ( L 2 ( Ω ) ) 2 , A = Δ Δ β I γ I + ν Δ , D A = ( H 0 2 ( Ω ) ) 2 ,
where H 0 2 ( Ω ) : = z H 2 ( Ω ) : z ( s ) = 0 , s Ω . Hence A C l ( Z ) .
Denote Λ 1 z = Δ z , D Λ 1 = H 0 2 ( Ω ) L 2 ( Ω ) . By { φ k : k N } denote an orthonormal in the sense of the inner product · , · in L 2 ( Ω ) eigenfunctions of the operator Λ 1 , which are enumerated in the non-increasing order of the corresponding eigenvalues { λ k : k N } taking in account their multiplicities.
Theorem 9.
Let c ( 1 , 2 ) , ν > 0 , β , γ , δ R , then there exist θ 0 ( π / 2 , π ) , a 0 0 , such that A A W ( θ 0 , a 0 ) ,
σ ( A ) = λ = 1 2 ( λ k ( 1 + ν ) + γ ± ( γ + λ k ( ν 1 ) ) 2 4 β λ k ) C : k N .
Proof. 
Using decomposition by the basis { φ k : k N } in the space L 2 ( Ω ) and the denotations μ k ± : = 1 2 ( λ k ( 1 + ν ) + γ ± ( γ + λ k ( ν 1 ) ) 2 4 β λ k ) , for
W ( λ ) = j = 1 n ω j λ α j + b c ω ( α ) λ α d α μ k ± , k N ,
obtain the operators
W ( λ ) I A = W ( λ ) I Δ Δ β I ( W ( λ ) γ ) I ν Δ ,
( W ( λ ) I A ) 1 = k = 1 W ( λ ) γ ν λ k λ k β W ( λ ) λ k · , φ k φ k ( W ( λ ) μ k + ) ( W ( λ ) μ k ) .
Since μ k + ν λ k and μ k λ k as k (or inversely, depending on the condition ν 1 or ν ( 0 , 1 ) ), for ν > 0 lim k arg μ k ± = π , therefore, there exists μ 0 = max k N Re μ k ± .
In the proof of Theorem 6 it was shown that for every ε > 0 there exists δ 1 such that for | λ | > δ we have | arg W ( λ ) c arg λ | < ε . Therefore, we can choose sufficiently close to π / 2 θ 0 ( π / 2 , π ) and large enough a 0 > δ sin 1 θ 0 , such that for all k N μ k ± S c θ 0 + 2 ε , a 0 , where c θ 0 + 2 ε < π for small enough ε > 0 . Then due to Lemma 1 | W ( λ ) | > C | λ | c ε > | λ | c ε / 2 , where ε 2 ( c 1 ) , hence W ( λ ) S c θ 0 + ε , a 0 . Then for λ S θ 0 , a 0
W ( λ ) γ ν λ k W ( λ ) μ k + 1 + | γ | | W ( λ ) | sin ε 1 + | γ | ( a 0 sin θ 0 ) c ε / 2 sin ε = C 1 ,
λ k W ( λ ) μ k 1 + 1 sin ε = C 2 ,
β W ( λ ) μ k + | β | ( a 0 sin θ 0 ) c ε / 2 sin ε = C 3 , W ( λ ) λ k W ( λ ) μ k = 1 .
Hence for all λ S θ 0 , a 0 , z ( L 2 ( Ω ) ) 2
R W ( λ ) ( A ) z L ( Z ) 2 C 4 sup k N 1 | W ( λ ) μ k | + 1 | W ( λ ) μ k + | 2 C 4 sin 2 ε | W ( λ ) a 0 | .
Thus, A A W ( θ 0 , a 0 ) . □
By Theorems 8 and 9 we obtain the corollary.
Corollary 1.
Let c ( 1 , 2 ) , ν > 0 , β , γ , δ R , κ ( 0 , 1 ] , f i C ( [ 0 , T ] ; H 2 ( Ω ) ) C κ ( [ 0 , T ] ; L 2 ( Ω ) ) , i = 1 , 2 . Then for all u 0 , u 1 , v 0 , v 1 H 0 2 ( Ω ) there exists a unique solution of problem (18)–(22).
Remark 6.
Analogously, but essentially simpler, we can study the initial boundary value problem
u ( s , 0 ) = u 0 ( s ) , u t ( s , 0 ) = u 1 ( s ) , s Ω ,
u ( s , t ) = 0 , ( s , t ) Ω × ( 0 , T ] ,
for the ultraslow diffusion equation
j = 1 n ω j D t α j u ( s , t ) + b c ω ( α ) D t α u ( s , t ) d α = Δ u ( s , t ) + f ( s , t ) , Ω × ( 0 , T ] .
For this aim, put Z = L 2 ( Ω ) , A = Λ 1 . A similar equation without the sum j = 1 n ω j D t α j u ( s , t ) was investigated in [7].

7. Conclusions

The Cauchy problem for equations in Banach spaces with a distributed fractional derivative and with a linear closed operator A at the unknown function is studied. The derivative is given by the Riemann–Stieltjes integral with respect to the order of the fractional differentiation, therefore, the considered class of equations includes equations with a distributed derivative, defined by a standard integral, or with a discretely distributed derivative, which was researched earlier. The notion of a k-resolving family of operators for the equation is introduced, and properties of such families are studied. It is shown that the existence of a 0-resolving family implies the existence of other k-resolving families, k = 1 , 2 , , m 1 . Necessary and sufficient conditions for the existence of an analytic 0-resolving family of operators in terms of the resolvent of the operator A is the key result of this work. The corresponding class of the operators is denoted by A W ( θ 0 , a 0 ) . The properties of analytic resolving families, generated by operators from this class, are investigated, and a perturbutaion theorem for such class of operators is proved. The unique solvability theorem for the inhomogeneous equation with a distributed fractional derivative and with A A W ( θ 0 , a 0 ) is obtained. Results of the work are applied to the research of an initial boundary value problem for a system of partial differential equations with a distributed time derivative in a general form.
There are a large number of various types of fractional derivatives, and in recent decades, new constructions of them have appeared: Riemann–Liouville derivative, Hadamard derivative, Marchot derivative, Dzhrbashyan–Nersesyan derivative, Prabhakar derivative, Caputo–Fabrizio and Atangana–Baleanu integro-differential operators, etc. Every construction of a fractional derivative corresponds to certain features of considered problems. The Gerasimov–Caputo derivative is the most studied (along with the Riemann–Liouville derivative) and mathematically simplest among such derivatives. The results obtained in this paper allow us to understand the features of equations with distributed derivatives given by the Stieltjes integral. This will allow us to move on to the study of equations with integrals with respect to the differentiation order of fractional derivatives of other types.

Author Contributions

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

Funding

This work is a part of the research carried out at the Ural Mathematical Center and supported by the Ministry of Science and Higher Education of the Russian Federation (agreement no. 075-02-2022-874).

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. The funders had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Rezapour, S.; Addas, M.I.; Etemad, S.; Dien, N.M. On a multi-point p-Laplacian fractional differential equation with generalized fractional derivatives. Math. Methods Appl. Sci. 2022. [Google Scholar] [CrossRef]
  2. Rezapour, S.; Souid, M.S.; Bouazza, Z.; Hussain, A.; Etemad, S. On the fractional variable order thermostat model: Existence theory on cones via piece-wise constant functions. Adv. Nonlinear Anal. Appl. 2022, 2022, 8053620. [Google Scholar] [CrossRef]
  3. Sintunavarat, W.; Turab, A. Mathematical analysis of an extended SEIR model of COVID-19 using the ABC-fractional operator. Math. Comput. Simul. 2022, 198, 65–84. [Google Scholar] [CrossRef]
  4. Nakhushev, A.M. On continual differential equations and their difference analogues. Sov. Math. Dokl. 1988, 37, 729–732. [Google Scholar]
  5. Caputo, M. Mean fractional order derivatives. Differential equations and filters. Ann. Dell’Universita Ferrara Sez. VII Sci. Mat. 1995, 41, 73–84. [Google Scholar] [CrossRef]
  6. Caputo, M. Distributed order differential equations modeling dielectric induction and diffusion. Fract. Calc. Appl. Anal. 2001, 4, 421–442. [Google Scholar]
  7. Kochubei, A.N. Distributed order calculus and equations of ultraslow diffusion. J. Math. Anal. Appl. 2008, 340, 252–280. [Google Scholar] [CrossRef]
  8. Sokolov, I.M.; Chechkin, A.V.; Klafter, J. Distributed-order fractional kinetics. Acta Phys. Pol. B 2004, 35, 1323–1341. [Google Scholar]
  9. 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]
  10. Bagley, R.L.; Torvik, P.J. On the existence of the order domain and the solution of distributed order equations. Part 1. Int. J. Appl. Math. 2000, 2, 865–882. [Google Scholar]
  11. Lorenzo, C.F.; Hartley, T.T. Variable order and distributed order fractional operators. Nonlinear Dyn. 2002, 29, 57–98. [Google Scholar] [CrossRef]
  12. Jiao, Z.; Chen, Y.; Podlubny, I. Distributed-Order Dynamic System. Stability, Simulations, Applications and Perspectives; Springer: London, UK, 2012. [Google Scholar]
  13. 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]
  14. 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]
  15. Pskhu, A.V. On the theory of the continual and integro-differentiation operator. Differ. Equ. 2004, 40, 128–136. [Google Scholar] [CrossRef]
  16. Pskhu, A.V. Partial Differential Equations of Fractional Order; Nauka Publ.: Moscow, Russia, 2005. (In Russian) [Google Scholar]
  17. Pskhu, A.V. Fractional diffusion equation with a discretely distributed differentiation operator. Sib. Elektron. Math. Rep. 2016, 13, 1078–1098. [Google Scholar]
  18. Umarov, S.; Gorenflo, R. Cauchy and nonlocal multi-point problems for distributed order pseudo-differential equations. Z. Für Anal. Ihre Anwend. 2005, 24, 449–466. [Google Scholar]
  19. 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]
  20. Efendiev, B.I. Steklov problem for a second-order ordinary differential equation with a continual derivative. Differ. Equ. 2013, 49, 450–456. [Google Scholar] [CrossRef]
  21. Efendiev, B.I. Lagrange formula for ordinary continual second-order differential equations. Differ. Equ. 2017, 53, 736–744. [Google Scholar] [CrossRef]
  22. 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. (In Russian) [Google Scholar]
  23. Fedorov, V.E.; Streletskaya, E.M. Initial-value problems for linear distributed-order differential equations in Banach spaces. Electron. J. Differ. Equ. 2018, 2018, 176. [Google Scholar]
  24. Fedorov, V.E.; Phuong, T.D.; Kien, B.T.; Boyko, K.V.; Izhberdeeva, E.M. A class of semilinear distributed order equations in Banach spaces. Chelyabinsk Phys. Math. J. 2020, 5, 343–351. [Google Scholar]
  25. Fedorov, V.E. On generation of an analytic in a sector resolving operators family for a distributed order equation. Zap. POMI 2020, 489, 113–129. [Google Scholar] [CrossRef]
  26. Fedorov, V.E. Generators of analytic resolving families for distributed order equations and perturbations. Mathematics 2020, 8, 1306. [Google Scholar] [CrossRef]
  27. Fedorov, V.E.; Filin, N.V. Linear equations with discretely distributed fractional derivative in Banach spaces. Tr. Instituta Mat. I Mekhaniki UrO RAN 2021, 27, 264–280. [Google Scholar]
  28. 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]
  29. Fedorov, V.E.; Du, W.-S.; Kostić, M.; Abdrakhmanova, A.A. Analytic resolving families for equations with distributed Riemann—Liouville derivatives. Mathematics 2022, 10, 681. [Google Scholar] [CrossRef]
  30. Gerasimov, A.N. A generalization of linear laws of deformation and its application to the problems of internal friction. Prikl. Mat. I Mekhanika 1948, 12, 251–260. (In Russian) [Google Scholar]
  31. Caputo, M. Linear model of dissipation whose Q is almost frequancy independent. II. Geophys. J. R. Astron. Soc. 1967, 13, 529–539. [Google Scholar] [CrossRef]
  32. 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]
  33. Podlubny, I. Fractional Differential Equations; Academic Press: Boston, MA, USA, 1999. [Google Scholar]
  34. Prüss, J. Evolutionary Integral Equations and Applications; Springer: Basel, Switzerland, 1993. [Google Scholar]
  35. Arendt, W.; Batty, C.J.K.; Hieber, M.; Neubrander, F. Vector-Valued Laplace Transforms and Cauchy Problems; Springer: Basel, Switzerland, 2011. [Google Scholar]
  36. 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]
  37. Hille, E.; Phillips, R.S. Functional Analysis and Semi-Groups; American Mathematical Society: Providence, RI, USA, 1957. [Google Scholar]
  38. Yosida, K. Functional Analysis; Springer: Berlin/Heidelberg, Germany, 1965. [Google Scholar]
  39. Kato, K. Perturbation Theory for Linear Operators; Springer: Berlin/Heidelberg, Germany, 1966. [Google Scholar]
  40. Goldstein, J.A. Semigroups and second-order differential equations. J. Funct. Anal. 1969, 4, 50–70. [Google Scholar] [CrossRef]
  41. Caginalp, G. An analysis of a phase field model of a free boundary. Arch. Ration. Mech. Anal. 1986, 92, 205–245. [Google Scholar] [CrossRef]
  42. Caginalp, G. Stefan and Hele—Shaw type models as asymptotic limits of the phase-field equations. Phys. Rev. A 1989, 39, 5887–5896. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sitnik, S.M.; Fedorov, V.E.; Filin, N.V.; Polunin, V.A. On the Solvability of Equations with a Distributed Fractional Derivative Given by the Stieltjes Integral. Mathematics 2022, 10, 2979. https://doi.org/10.3390/math10162979

AMA Style

Sitnik SM, Fedorov VE, Filin NV, Polunin VA. On the Solvability of Equations with a Distributed Fractional Derivative Given by the Stieltjes Integral. Mathematics. 2022; 10(16):2979. https://doi.org/10.3390/math10162979

Chicago/Turabian Style

Sitnik, Sergey M., Vladimir E. Fedorov, Nikolay V. Filin, and Viktor A. Polunin. 2022. "On the Solvability of Equations with a Distributed Fractional Derivative Given by the Stieltjes Integral" Mathematics 10, no. 16: 2979. https://doi.org/10.3390/math10162979

APA Style

Sitnik, S. M., Fedorov, V. E., Filin, N. V., & Polunin, V. A. (2022). On the Solvability of Equations with a Distributed Fractional Derivative Given by the Stieltjes Integral. Mathematics, 10(16), 2979. https://doi.org/10.3390/math10162979

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop