Subordination Principle for a Class of Fractional Order Differential Equations

The fractional order differential equation u′(t) = Au(t) + γD t Au(t) + f(t), t > 0, u(0) = a ∈ X is studied, whereA is an operator generating a strongly continuous one-parameter semigroup on a Banach space X , D t is the Riemann–Liouville fractional derivative of order α ∈ (0, 1), γ > 0 and f is an X-valued function. Equations of this type appear in the modeling of unidirectional viscoelastic flows. Well-posedness is proven, and a subordination identity is obtained relating the solution operator of the considered problem and the C0-semigroup, generated by the operator A. As an example, the Rayleigh–Stokes problem for a generalized second-grade fluid is considered.


Introduction
Consider the fractional order differential equation: where D α t is the Riemann-Liouville fractional derivative of order α ∈ (0, 1), γ > 0, A is an unbounded closed linear operator defined on a Banach space X and f is an X-valued function.
Our motivation for studying this equation comes from recent works where related problems appear in the modeling of unidirectional viscoelastic flows.For example, if the operator A is some realization of the Laplace operator, then Problem (1) is the Rayleigh-Stokes problem for a generalized second-grade fluid; see e.g., [1][2][3].Exact solutions of this problem in the form of eigenfunction expansion on a bounded space domain are obtained in [4,5].Numerical analysis of Problem (1) with A being the one-or two-dimensional Laplacian with Dirichlet boundary conditions is carried out in [6][7][8][9][10][11].In [12], a compact Duhamel-type representation of the solution is obtained and used for its numerical computation.
In addition, let us note that the governing equation in (1) can be also considered as a distributed order equation in the Riemann-Liouville sense: t Au(t) dβ, t > 0 with weight function µ(β) = δ(β) + γδ(β − α); see [13] and the references cited there.If γ = 0 and f ≡ 0, Problem (1) reduces to the classical abstract Cauchy problem: There is a vast amount of literature devoted to this problem and its equivalent formulation: the theory of strongly continuous (C 0 ) one-parameter semigroups of operators (see, e.g., [14]).In the present paper, it is assumed that the operator A is a generator of a C 0 -semigroup, i.e., that the classical Cauchy Problem (2) is well-posed.Under this assumption, it will be proven that Problem (1) is well-posed, and a relationship between the corresponding solution operators of these two problems will be established in the form: where S(t) and T (t) are the solution operators of the considered Problem (1) and the classical Cauchy Problem (2), respectively.It appears that the function ϕ(t, τ ) is a probability density with respect to both t and τ .Such a relationship between two problems is called the subordination principle; see [15], Ch. 4. Representations of the form of (3) are useful in the study of differential equations of fractional order, in particular for the understanding of the regularity and the asymptotic behavior of the solution of the subordinate problem.For example, the subordination principle for fractional evolution equations with the Caputo derivative (see [16], Ch. 3) has been successfully applied to inverse problems [17], for asymptotic analysis of diffusion wave equations [18], for the study of stochastic solutions [19], semilinear equations of fractional order [20], systems of fractional order equations [21], nonlocal fragmentation models with the Michaud time derivative [22], etc.This gives the author the motivation to present an analogous principle for Problem (1).
The paper is organized as follows.Section 2 contains preliminaries.In Section 3, Problem (1) is reformulated as a Volterra integral equation, and the properties of its kernel are studied, as well as the solution of the scalar version of the problem, where A is a negative constant.Section 4 contains the main result: the subordination principle for Problem (1) and some corollaries.As an example, the Rayleigh-Stokes problem for a generalized second-grade fluid is considered in Section 5.

Preliminaries
The sets of positive integers, real and complex numbers are denoted as usual by N, R and C, respectively, and N 0 = N ∪ {0}, R + = [0, ∞).Denote by Σ θ the sector: For the sake of brevity, the following notation is used: where Γ(.) is the Gamma function.
Let * denote the Laplace convolution: Denote by J α t the Riemann-Liouville fractional integral of order α: and by D α t the Riemann-Liouville fractional derivative of order α: We use the following standard notations for the Laplace transform: Application of the Laplace transform to the Riemann-Liouville fractional differential operator D α t gives (see, e.g., [23]): For functions continuous for t > 0 and such that f (0) is finite, the second term in (7) vanishes, and the identity reduces to: A C ∞ function f : (0, ∞) → R is said to be completely monotone if: The characterization of completely monotone functions is given by Bernstein's theorem (see, e.g., [24]), which states that a function f : (0, ∞) → R is completely monotone if and only if it can be represented as the Laplace transform of a nonnegative measure.
A C ∞ function f : (0, ∞) → R is called a Bernstein function if it is nonnegative and its first derivative f (t) is a completely monotone function.The classes of completely monotone functions and Bernstein functions will be denoted by CMF and BF.
The asymptotic behavior of the three-parameter Mittag-Leffler function (10) can be obtained from the identity [28]: Other kinds of asymptotic estimates were provided in [29], and on their basis, the convergence of series in Functions (10) in the complex domain, similar to these appearing further in (27) and (34), is studied in the recent papers [25,30].
Recall also the Laplace transform pair [23]: Denote by Φ β (z), β ∈ (0, 1) the following function of the Wright type, also called the Wright M-function or Mainardi function (see, e.g., [31]): The following relationship with the Mittag-Leffler function holds: In particular, this identity implies that Φ β (t) is a probability density function: Let X be a complex Banach space with norm . .Let A be a closed linear unbounded operator in X with dense domain D(A), equipped with the graph norm .A , x A := x + Ax .Denote by (A) the resolvent set of A and by R(s, A) = (sI − A) −1 the resolvent operator of A.
By integrating both sides of the governing equation in (1), we recast Problem (1) with f ≡ 0 into a Volterra integral equation: where the kernel k(t) is specified later; see (20).Here, we recall some definitions and basic theorems, given in [15], concerning abstract Volterra integral equations.
) and ( 15) holds on R + .Problem ( 15) is said to be well-posed if for each a ∈ D(A), there is a unique strong solution u(t; a) of ( 15), and For a well-posed problem, the solution operator S(t) is defined by: Then, the Laplace transform of the solution operator S(t) of Problem ( 15): is given by The generation theorem ( [15], Theorem 1.3) states that Problem ( 15) is well-posed with solution operator S(t) satisfying S(t) ≤ M, t ≥ 0, if and only if: Note that the classical Cauchy Problem (2) can be also rewritten as Volterra Equation ( 15) with kernel k(t) = 1, g(s) = s.In this case, the generation theorem is known as the classical Hille-Yosida theorem.

Integral Reformulation of the Problem
Let first f ≡ 0. Assume u, Au ∈ C(R + , X).Then, integrating both sides of the governing equation in (1), by the use of ( 5) and the identity (J 1−α t Au)(0) = 0, we obtain: that is the Volterra integral Equation ( 15) with kernel k(t) given by: where the function ω α is defined in (4).Conversely, differentiating both sides of (19) and using that: we get back the governing equation in (1).Since (k * Au)(0) = 0, the initial condition is also satisfied.The above observations give rise to following definition.
Definition 2. Problem ( 1) is said to be well-posed if the Volterra integral Equation ( 19) is well-posed.
For well-posed problems, the solution operator of the integral Equation ( 19) is called a solution operator of Problem (1).
Let Problem (1) be well-posed, and let S(t) be the corresponding solution operator.Consider the inhomogeneous Problem (1) with f ∈ L 1 (R + ; X).In the same way as above, it can be rewritten as the integral equation: Then, by the variation of the parameters formula (see e.g., [15]), the solution of the inhomogeneous Problem ( 1) is given by: Therefore, it is essential for the homogeneous, as well as for the inhomogeneous problem, to study the solution operator S(t).
We begin with summarizing some properties of the kernel k(t) defined in (20), relevant for further study.Along with the kernel k(t), the related function: is also of interest, since it appears in (17).Here, k is the Laplace transform of k(t), and the Laplace transform pair ( 6) is used.

Subordination Principle
Assume the operator A generates a bounded C 0 semigroup T (t).The main goal of this paper is to prove that in this case, Problem (1) is well-posed, and its solution operator S(t) satisfies the relationship: with an appropriate function ϕ(t, τ ).This is the so-called subordination principle: Problem (1) is subordinate to the classical Cauchy Problem (2).Let the function ϕ(t, τ ) be such that its Laplace transform with respect to t satisfies: where g(s) is defined in (22).The reason for this is that then the operator S(t), defined by (28), will satisfy ( 16) and (17).Indeed, by ( 28) and ( 29) and the identity for the Laplace transform of a C 0 -semigroup Then, by the uniqueness of the Laplace transform, S(t) will be the solution operator of Problem (1).Identity (29) implies that the function ϕ(t, τ ) can be found by the inverse Laplace integral: Let us check that the function ϕ(t, τ ) is well defined in this way.According to Theorem 3 (e), {s} > 0 implies {g(s)} > 0.More precisely, if s = re iθ , then: Hence, when r → ∞, |θ| → π/2, the dominant term of {g(s)} is r 1−α sin απ/2 > 0. This together with the estimate (f) of Theorem 3 shows that the integral in ( 31) is absolutely convergent.
We are ready to formulate the main result of this paper.
The definition (28), the estimate (32) and the properties, ( 35) and (36), imply: 33) is established.Next, we deduce the strong continuity of S(t) at the origin from the strong continuity of T (t) at the origin: On the basis of the dominated convergence theorem and by the change of variables σ = t α−1 τ in (28), we obtain: For the function under the integral sign, we get from (34): and thus: where Φ β (z), β ∈ (0, 1), is the Mainardi function (13).Therefore, (40) together with (39) and the integral identity in ( 14) for the Mainardi function imply: In this way, we proved that S(t), defined by (28), is a strongly continuous operator-valued function, satisfying (33).Moreover, in (30), we proved that the Laplace transform of S(t) satisfies: After easily justified differentiation under the integral sign in (41), we obtain from (33) the estimates (18) and, thus, the well-posedness of Problem (1).Then, Identity (16) implies by the uniqueness of the Laplace transform that S(t) is exactly the solution operator of (1).The proof of the theorem is completed.
The positivity of the function ϕ(t, τ ) in (28) has an important implication related to the positivity of the solution operator S(t).
Let X be an ordered Banach space (for a simple introduction, see, e.g., [14], p. 353).For example, such are the spaces of type L p (Ω) or C 0 (Ω) for some Ω ∈ R d , d ∈ N, with the canonical ordering: a function a ∈ X is positive (in symbols: a ≥ 0) if a(x) ≥ 0 for (almost) all x ∈ Ω.
A solution operator S(t) in an ordered Banach space X is called positive if a ≥ 0 implies S(t)a ≥ 0 for any t ≥ 0. In other words, the positivity of a solution operator means that the positivity of the initial condition is preserved in time.
As a direct consequence of the subordination identity (28), we obtain: Corollary 6.If the operator A is a generator of positive C 0 -semigroup, then the solution operator S(t) of Problem (1) is positive.
Another implication of the subordination principle follows from Identities ( 28) and (37): In the case of a contraction solution operator of Problem (1), there holds some inversion of the subordination principle in the following sense.Let us suppose that for some α ∈ (0, 1), Problem ( 1) is well-posed and admits a contraction solution operator S(t), S(t) ≤ 1, t ≥ 0.Then, (18) with n = 0 gives: R(g(s), A) ≤ 1 g(s) , s > 0 Since g(s) is a Bernstein function for s > 0, it is positive and strongly increasing.This, together with g(0) = 0 and g(+∞) = +∞, implies that g(s) : R + → R + is a one-to-one mapping.Therefore, and the Hille-Yosida theorem implies that A generates a contraction semigroup.Then, by the subordination principle in Theorem 5, it follows that for any α ∈ (0, 1), Problem ( 1) is well-posed and admits a contraction solution operator.
Note that, if the C 0 -semigroup generated by the operator A is moreover analytic on some sector Σ θ , then the solution operator S(t) of Problem (1) admits analytic extension to the same sector Σ θ .This result follows from Corollary 2.4 in [15] and is based on Properties (a) and (b) of the kernel k(t), given in Theorem 3.
Next, we compare the subordination principle formulated in Theorem 5 with the analogous subordination principle for fractional evolution equations with the Caputo fractional derivative (see [16]), which can be written in the form: Note that for such equations, the function ϕ(t, τ ) in ( 28) is given by ϕ(t, τ ) = t −α Φ α (τ t −α ), where Φ α is the Mainardi function (13).In the case of Equation (42), the subordinate solution operator S(t) is always analytic in some sector without assuming the analyticity of the C 0 -semigroup T (t).It seems that this is not true for the here considered Equation (1).This difference is due to the following: the kernel corresponding to Equation ( 42) is given by k(t) = ω α (t), and it is sectorial kernel of angle απ/2 < π/2, while the kernel (20) associated with Problem (1) is sectorial of angle π/2.A kernel k(t) is called sectorial of angle θ > 0 if ( [15], Ch. 3): Another difference is that identity (37) does not hold for Equation (42).

Example
The following problem, which is a particular case of Problem (1), is considered in [11].Let Ω ⊂ R d be a bounded domain with sufficiently smooth boundary ∂Ω, and T > 0. Consider the initial-boundary value problem: ∂ ∂t u(x, t) = ∆u(x, t) + γD α t ∆u(x, t) + f (x, t), x ∈ Ω, t ∈ (0, T ); u(x, t) = 0, x ∈ ∂Ω, t ∈ (0, T ); (43) where ∆ is the Laplace operator acting on space variables.Let X = L 2 (Ω).Define the operator n=1 is the eigensystem of the operator A, then 0 < λ 1 ≤ λ 2 ≤ ..., λ n → ∞ as n → ∞, and {φ n } ∞ n=1 form an orthonormal basis of X = L 2 (Ω).Applying eigenfunction decomposition and Laplace transform, the solution of Problem (43) is obtained in the form: where a n = (a, φ n ), f n (t) = (f (., t), φ n ), and u n (t) = u(t, λ n ): the solution of the scalar Equation ( 23) with λ = λ n , n ∈ N. Therefore, the solution operator for this problem admits the representation: The C 0 -semigroup T (t) generated by the operator A (corresponding to the solution of (43) with γ = 0) is given by: Theorem 5 implies the well-posedness of Problem (43) from the well-posedness of the corresponding initial boundary value problem for the diffusion Equation ((43) with γ = 0) and gives the relationship between the solution operators S(t) and T (t).Note that the subordination Identity (28) in this case can be obtained from Identity (38).Since the C 0 -semigroup T (t) gives the solution of a diffusion equation, it is positive, i.e., it preserves the positivity of the initial function a. Applying Corollary 6, we deduce that the solution operator of Problem (43) is also positive: i.e., if a(x) ≥ 0 for a.a.x ∈ Ω, then the solution of Problem (43) with f ≡ 0 given by u(x, t) = S(t)a is positive for a.a.x ∈ Ω, and all t ≥ 0.

Conclusions
An abstract fractional order differential equation is studied, which contains as particular cases the Rayleigh-Stokes problem for a generalized second-grade fluid with a fractional derivative model and a two-term fractional diffusion equation in the Riemann-Liouville sense.Well-posedness is proven, and a subordination identity is presented relating the solution operator of the considered equation and the solution operator of the classical Cauchy problem.The subordination identity contains a special function represented as a series of three-parameter Mittag-Leffler functions.As a limiting case, the Mainardi function is recovered.
A feasible generalization of the presented results is an extension of the subordination principle to the general case when T (t) is an exponentially-bounded C 0 -semigroup: T (t) ≤ Ce ωt , ω > 0.Then, estimates for the solution of the scalar equation with A = ω > 0 will be necessary.

T0
u n (t) dt < 1/λ n .Applying the Young inequality for the convolution k * f L 2 ≤ k L 1 f L 2 , it follows: