Next Article in Journal
Interval Model of the Efficiency of the Functioning of Information Web Resources for Services on Ecological Expertise
Next Article in Special Issue
Unified Approach to Fractional Calculus Images of Special Functions—A Survey
Previous Article in Journal
The Consistency of the CUSUM-Type Estimator of the Change-Point and Its Application
Previous Article in Special Issue
Intrinsic Discontinuities in Solutions of Evolution Equations Involving Fractional Caputo–Fabrizio and Atangana–Baleanu Operators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

The General Fractional Derivative and Related Fractional Differential Equations

1
Department of Mathematics, Physics, and Chemistry, Beuth Technical University of Applied Sciences Berlin, Luxemburger Str. 10, 13353 Berlin, Germany
2
Graduate School of Mathematical Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8914, Japan
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(12), 2115; https://doi.org/10.3390/math8122115
Submission received: 29 September 2020 / Revised: 9 November 2020 / Accepted: 23 November 2020 / Published: 26 November 2020
(This article belongs to the Special Issue Fractional Integrals and Derivatives: “True” versus “False”)

Abstract

:
In this survey paper, we start with a discussion of the general fractional derivative (GFD) introduced by A. Kochubei in his recent publications. In particular, a connection of this derivative to the corresponding fractional integral and the Sonine relation for their kernels are presented. Then we consider some fractional ordinary differential equations (ODEs) with the GFD including the relaxation equation and the growth equation. The main part of the paper is devoted to the fractional partial differential equations (PDEs) with the GFD. We discuss both the Cauchy problems and the initial-boundary-value problems for the time-fractional diffusion equations with the GFD. In the final part of the paper, some results regarding the inverse problems for the differential equations with the GFD are presented.

1. Introduction

In functional analysis, the integral operators with the weakly singular kernels have been an important topic for research for many years. They are defined in the form
( T f ) ( t ) = Ω K ( t , τ ) f ( τ ) d τ ,
where Ω is an open subset of R n and the kernel K = K ( t , τ ) is a real or complex valued continuous function on Ω ¯ × Ω ¯ \ D , D = { ( t , t ) : t Ω ¯ } being the diagonal of Ω ¯ × Ω ¯ that satisfies the condition
| K ( t , τ ) | c | t τ | α with α < n .
In case Ω is a bounded domain, the operator (1) is called the Schur integral operator. It is compact from C ( Ω ¯ ) to C ( Ω ¯ ) [1]. If, additionally, α < n q , 1 p + 1 q = 1 , 1 p < + , then (1) is a Hilbert-Schmidt operator that is compact from L p ( Ω ) to C ( Ω ¯ ) .
In Fractional Calculus (FC), the operators of type (1) with special weakly singular kernels are studied on both bounded and unbounded domains. For instance, the classical left-hand sided Riemann-Liouville fractional integral of order α R + of a function f on a finite or infinite interval ( a , b ) is defined as follows:
I a + α f ( t ) = 1 Γ ( α ) a t ( t τ ) α 1 f ( τ ) d τ , t ( a , b ) .
In the case α = 0 , this integral is interpreted as the identity operator
I a + 0 f ( t ) = f ( t )
because of the relation ([2])
lim α 0 + I a + α f ( t ) = f ( t ) ,
that is valid in particular for f L 1 ( a , b ) in every Lebesgue point of f, i.e., almost everywhere on ( a , b ) .
Evidently, the Riemann-Liouville fractional integral is a generalization of the well-known formula for the n-fold definite integral
I a + n f ( t ) = a t d τ a t d τ a t f ( τ ) d τ = 1 ( n 1 ) ! a t ( t τ ) n 1 f ( τ ) d τ , n N .
In a certain sense, this generalization is unique. In the case of a finite interval (without any restriction of generality, we fix the interval ( 0 , 1 ) ), the following result was proved in [3]:
Let E be the space L p ( 0 , 1 ) , 1 p < + , or C [ 0 , 1 ] . Then there exists precisely one family I α , α > 0 of operators on E satisfying the following conditions:
(CM1)
( I 1 f ) ( t ) = 0 t f ( τ ) d τ , f E (interpolation condition),
(CM2)
( I α I β f ) ( t ) = ( I α + β f ) ( t ) , α , β > 0 , f E (index law),
(CM3)
α I α is a continuous map of R + into the space L ( E ) of the linear bounded operators from E to E for some Hausdorff topology on L ( E ) , weaker than the norm topology (continuity),
(CM4)
f E and f ( t ) 0 (a.e. for E = L p ( 0 , 1 ) ) ⇒ ( I α f ) ( t ) 0 (a.e. for E = L p ( 0 , 1 ) ) for all α > 0 (non-negativity).
That family is given by the Riemann-Liouville Formula (3) with a = 0 and b = 1 . From the present viewpoint ([4]), the conditions (CM1)–(CM4) are very natural for any definition of the fractional integrals defined on a finite interval. As proved in [3], they are also sufficient for uniqueness of the family of the Riemann-Liouville fractional integrals. Thus, in this sense, the Riemann-Liouville fractional integrals are the only “right” one-parameter fractional integrals defined on a finite one-dimensional interval.
The problem regarding the “right” fractional derivatives is more delicate and has no unique solution. Presently, the main approach for introducing the fractional derivatives is to define them as the left-inverse operators to the fractional integrals ([4,5,6]). However, even for the Riemann-Liouville fractional integral, there exist infinitely many different families of operators that fulfill this property ([6]). In particular, for 0 < α 1 , the Riemann-Liouville fractional derivative
( D R L α f ) ( t ) = d d t ( I 0 + 1 α f ) ( t ) ,
the Caputo fractional derivative
( D C α f ) ( t ) = ( I 0 + 1 α d f d τ ) ( t ) ,
and the Hilfer fractional derivative
( D H α f ) ( t ) = ( I 0 + β ( 1 α ) d d τ I 0 + ( 1 α ) ( 1 β ) f ) ( t ) , 0 β 1
are the left-inverse operators to the Riemann-Liouville fractional integral on the suitable nontrivial spaces of functions including the space of the absolutely continuous functions on [ 0 , 1 ] , i.e., the Fundamental Theorem of FC holds true ([6]):
( D X α I 0 + α f ) ( t ) = f ( t ) , t [ 0 , 1 ] , X { R L , C , H } .
Moreover, in [6], infinitely many other families of the fractional derivatives in the sense of Formula (9) called the nth level fractional derivatives were introduced. Let the parameters γ 1 , γ 2 , , γ n R satisfy the conditions
0 γ k and α + s k k with s k : = i = 1 k γ i , k = 1 , 2 , , n .
The nth level fractional derivative of order α , 0 < α 1 and type γ = ( γ 1 , γ 2 , , γ n ) is defined as follows:
( D n L α , ( γ ) f ) ( t ) = k = 1 n ( I 0 + γ k d d τ ) ( I 0 + n α s n f ) ( t ) .
This derivative satisfies the Fundamental Theorem of FC, i.e., the relation
( D n L α , ( γ ) I 0 + α f ) ( t ) = f ( t ) , t [ 0 , 1 ]
holds true on a nontrivial space of functions (see [6] for details).
To keep an overview of these and many other fractional derivatives, it is very natural to consider some general integro-differential operators of convolution type and to clarify the question under what conditions can they be interpreted as a kind of the fractional derivatives. In particular, one expects that for these derivatives and the appropriate defined fractional integrals the Fundamental Theorem of FC holds true. Moreover, for the sake of possible applications, one wants to keep some fundamental properties of solutions to the differential equations with these derivatives. In particular, the property of complete monotonicity of solutions to the appropriate relaxation equation or the positivity of the fundamental solution to the Cauchy problem for the fractional diffusion equation with the time-derivatives of this type.
In the theory of the abstract Volterra integral equations in the Banach spaces, the evolution equations including the integro-differential operators of convolution type
( D k u ) ( t ) = d d t 0 t k ( t τ ) u ( τ ) d τ , t [ 0 , T ] , 0 < T +
have been a subject for research for more than a half century. In particular, in [7] (see also the references therein), the abstract Volterra integral equations including the operators (13) with the completely positive kernels k L 1 ( 0 , T ) have been studied. This class of the kernels can be characterized as follows: A function k L 1 ( 0 , T ) is completely positive on [ 0 , T ] if and only if there exist a 0 and l L 1 ( 0 , T ) , non-negative and non-increasing, satisfying the relation
a k ( t ) + 0 t k ( t τ ) l ( τ ) d τ = 1 , t ( 0 , T ] .
In particular, the completely monotone kernels are completely positive. The notion of completely positive kernels originated from the “positivity preserving property” that is valid for the corresponding Volterra integral equations in the case of the Banach space X = R with the usual norm.
In [8,9], the properties of the appropriate defined weak solutions to the linear and quasi-linear evolutionary partial integro-differential equations of second order with the time-operators of type (13) in the form
( D k u ) ( t ) = d d t 0 t k ( t τ ) ( u ( τ ) u 0 ) d τ , t R +
were addressed in the case of the kernels k L l o c ( R + ) that satisfy the following conditions:
(Z1)
The kernel k is non-negative and non-increasing on R + ,
(Z2)
There exists a kernel l L l o c ( R + ) such that 0 t k ( t τ ) l ( τ ) d τ = 1 , t R + .
An important example of a kernel k that satisfies the conditions (Z1)–(Z2) is the following generalization of the Riemann-Liouville kernel ([8]):
k ( t ) = h 1 α ( t ) exp ( μ t ) , μ 0 , 0 < α < 1 ,
where
h β ( t ) = t β 1 Γ ( β ) , t > 0 , β > 0 .
For this kernel k, the kernel l from the property (Z2) takes the form
l ( t ) = h α ( t ) exp ( μ t ) + μ 0 t h α ( τ ) exp ( μ τ ) d τ , t > 0 .
However, in [8,9] and in earlier publications, the operators of type (13) were not interpreted as a kind of the generalized fractional derivatives. In particular, no construction of the corresponding fractional integral was presented and no conditions that ensure the physically relevant properties of solutions to the time-fractional differential equations including these derivatives were suggested. Both tasks along with a series of other useful properties were addresses in [10,11,12,13], where a very nice theory of the general fractional derivative of type (13) was developed and applied for studying properties of the ordinary and partial differential equations with this derivative. In this survey, we present some selected results obtained in these and other related publications.
The rest of the paper is organized as follows. In Section 2, we introduce the GFD and the related fractional integral and discuss some of their basic properties. Section 3 is devoted to the Cauchy problems for the fractional ODEs with the GFD. In particular, the fractional relaxation equation and properties of its solution and the fractional growth equation and long time asymptotic of its solution are considered [10,11,13]. Moreover, existence and uniqueness of solutions to the Cauchy problem for the nonlinear fractional ODEs with the GFD and their continuous dependence on the problem data are also addressed following [14]. In Section 4, we present some results regarding the fractional PDEs with the GFD. We start with the Cauchy problem for the linear fractional diffusion equation and address its well-posedness with the focus on an interpretation of its fundamental solution as a probability density function [10,11]. Then we proceed with a treatment of the initial-boundary-value problems for the time-fractional diffusion equation including the GFD. Based on a suitable estimate for the GFD of a function at its maximum point, a weak maximum principle for the general time-fractional diffusion equation with the GFD is deduced. Then, following [15], the maximum principle is employed to show uniqueness of the strong and the weak solutions to the initial-boundary-value problems for the general time-fractional diffusion equations. Existence of a weak solution in the sense of Vladimirov [16] is also discussed. Finally, some important results from the recent publications [17,18,19] regarding inverse problems for the fractional differential and integral equations with the GFD are shortly presented.

2. General Fractional Derivative and Integral

Following [10], in this paper we consider the GFD of the Riemann-Liouville type in the form (compare to (13))
( D k R L f ) ( t ) = d d t 0 t k ( t τ ) f ( τ ) d τ
and of the Caputo type in the form
( D k C f ) ( t ) = 0 t k ( t τ ) f ( τ ) d τ ,
where k is a non-negative locally integrable function. For an absolutely continuous function f satisfying f L l o c 1 ( R + ) , the relation (compare to (15))
( D k C f ) ( t ) = ( D k R L ( f f ( 0 ) ) ) ( t ) = ( D k R L f ) ( t ) k ( t ) f ( 0 ) = d d t 0 t k ( t τ ) f ( τ ) d τ k ( t ) f ( 0 )
between the Caputo and Riemann-Liouville types of GFDs holds true. In [10], the Caputo type GFD was introduced in form (21) that is well defined for a lager class of functions (in particular, for absolutely continuous functions) compared to the definition (20) that requires the inclusion f L l o c 1 ( R + ) . In what follows, we mainly address the GFD of the Caputo type in the sense of the right-hand side of Formula (21).
The Riemann-Liouville and Caputo fractional derivatives defined by (6) and (7), respectively, are particular cases of the GFDs (19) and (20) with the kernel
k ( t ) = h 1 α ( t ) , 0 < α < 1 ,
the power function h β being defined by (17). Other important particular cases of (19) and (20) are the multi-term fractional derivatives and the fractional derivatives of the distributed order. They are generated by (19) and (20) with the kernels
k ( t ) = k = 1 n a k h 1 α k ( t ) , 0 < α 1 < < α n < 1 , a k R , k = 1 , , n
and
k ( t ) = 0 1 h 1 α ( t ) d ρ ( α ) ,
where ρ is a Borel measure on [ 0 , 1 ] .
Even if the operators (19) and (20) have been employed in the theory of the abstract Volterra integral equations for many years, the main advantage of the Kochubei’s approach was to establish a connection of these operators to FC and to introduce a special class of the kernels that ensures both existence of the corresponding fractional integrals and physically relevant properties of the fractional differential equations with these time-fractional derivatives. Moreover, the results presented in [10] and in the subsequent publications were derived using a completely different technique, namely the theory of the complete Bernstein functions ([20]).
The kernels of the GFDs (19) and (21) considered in [10] satisfy the following conditions:
(K1)
The Laplace transform k ˜ of k,
k ˜ ( p ) = ( L k ) ( p ) = 0 k ( t ) e p t d t
exists for all p > 0 ,
(K2)
k ˜ ( p ) is a Stieltjes function,
(K3)
k ˜ ( p ) 0 and p k ˜ ( p ) as p ,
(K4)
k ˜ ( p ) and p k ˜ ( p ) 0 as p 0 .
In what follows, we denote the set of the kernels that satisfy the conditions (K1)–(K4) by K . As we see, the condition of type (Z2) (Sonine condition) does not belong to the set of the conditions (K1)–(K4). However, it is one of the consequences of these conditions and especially of the strong condition (K2). Roughly speaking, a function defined on R + is a Stieltjes function if it can be represented as a restriction of the Laplace transform of a completely monotone function to the real positive semi-axis. Any completely monotone function is non-negative and thus any Stieltjes function is completely monotone as the Laplace transform of a non-negative function. For the strict definition and properties of the Stieltjes functions see e.g., [20,21]. The kernel functions (22) and (23) as well as the function (24) under some suitable conditions on the measure ρ belong to the class K . In [10], another example of a function k K was introduced in terms of its Laplace transform k ˜ ( p ) = p 1 log ( 1 + p β ) , 0 < β < 1 .
As shown in [10], for each k K , there exists a completely monotone function κ such that the Sonine condition holds true:
( k κ ) ( t ) = 0 t k ( t τ ) κ ( τ ) d τ = 1 .
Henceforth by
( I k f ) ( t ) = 0 t κ ( t τ ) f ( τ ) d τ
we denote a general fractional integral (GFI) with the kernel κ associated with the kernel k of the GFD by means of the relation (25).
The notion of the GFI is justified by the following Fundamental Theorem of FC:
Theorem 1
([10]). If f is a locally bounded measurable function on R + , then
( D k C I k f ) ( t ) = f ( t ) .
If f is absolutely continuous on [ 0 , + ) , then
( I k D k C f ) ( t ) = f ( t ) f ( 0 ) .
The Formula (27) and the relation (21) between the GFDs of the Caputo and the Riemann-Liouville types lead to the identity
( D k R L I k f ) ( t ) = f ( t ) ,
i.e., the Riemann-Liouville GFD is also a left inverse operator to the GFI defined by (26).
In the case of the Riemann-Liouville and the Caputo fractional derivatives that are particular cases of the GFDs (19) and (20), respectively, with the power function kernel k defined by (22), the kernel κ in the GFI (25) is also the power function κ ( t ) = h α ( t ) and thus in this case the GFI (26) is nothing else as the conventional Riemann-Liouville fractional integral.
As shown in [22], the functions that satisfy the Sonine condition (25) cannot be continuous at the point t = 0 and thus the “new fractional derivatives” with the continuous kernels introduced recently in the FC literature do not belong to the class of the GFDs that are discussed in this paper.
In the next sections, we consider other physically relevant properties of the GFD including complete monotonicity of the solutions to the fractional relaxation equation with this derivative, positivity of the fundamental solution to the Cauchy problem for the fractional diffusion equation with the time-derivative in form of the GFD, and a maximum principle for the initial-boundary-value problems for the fractional diffusion equation with this derivative.

3. Fractional ODEs with the GFD

3.1. Fractional Relaxation Equation

In this subsection, we consider the fractional relaxation equation
( D k C u ) ( t ) = λ u ( t ) , λ > 0 , t > 0
subject to the initial condition
u ( 0 ) = 1 .
As discussed in [23] (see also references therein), in the framework of the linear viscoelasticity models, the solutions to the relaxation equations are expected to be completely monotone. Only in this case the relaxation processes can be interpreted as superpositions of (infinitely many) elementary, i.e., exponential, relaxation processes. For the fractional relaxation equation with the GFD of the Caputo type, the following result holds true:
Theorem 2
([10]). Let the kernel k of the Caputo type GFD belong to the class K .
Then the Cauchy problem (30), (31) has a unique solution u λ = u λ ( t ) , continuous on [ 0 , + ) , infinitely differentiable and completely monotone on R + .
We remind the readers that a function u : R + R is called completely monotone if the conditions
( 1 ) n u ( n ) ( t ) 0 , t > 0 , n = 0 , 1 , 2
hold true.
In the case of the Cauchy problem (30), (31) with the Caputo fractional derivative (7), the solution can be expressed in terms of the Mittag-Leffler function
u λ ( t ) = E α , 1 ( λ t α ) ,
where E α , β stands for the two-parameters Mittag-Leffler function that is defined by the following convergent series:
E α , β ( z ) = k = 0 z k Γ ( α k + β ) , α > 0 , β , z C .
It is worth mentioning that the fractional relaxation equations with the Riemann-Liouville fractional derivative (6) and the Hilfer derivative (8) have the solutions
u λ ( t ) = t α 1 E α , α ( λ t α ) ,
and
u λ ( t ) = t α + γ 1 1 E α , α + γ 1 ( λ t α ) ,
respectively, provided that we set suitable initial conditions. These solutions are continuous, infinitely differentiable and completely monotone on R + , but have an integrable singularity at the point t = 0 . As to the relaxation equation with the nth level fractional derivative (11), its solution is given by the following theorem:
Theorem 3
([24]). The fractional relaxation equation
( D n L α , ( γ ) u ) ( t ) = λ u ( t ) , λ > 0 , t > 0
subject to the initial conditions
i = k + 1 n I 0 + γ i d d t I 0 + n α s n u ( 0 ) = u k , k = 1 , , n
has a unique solution, continuous and infinitely differentiable on R + , given by the formula
u ( t ) = k = 1 n u k t α + s k k E α , α + s k k + 1 ( λ t α ) , s k = i = 1 k γ i .
If the initial conditions are non-negative ( u k 0 , k = 1 , , n in (35)) and the inequalities
k 1 s k = i = 1 k γ i , k = 1 , , n
hold true, the solution (36) is completely monotone on R + .
In [10], an important probabilistic interpretation of Theorem 2 has been provided. Let D ( t ) be a subordinator of the Lévy process with the Laplace exponent Ψ = Ψ ( s ) ([25]):
E e s D ( t ) = e t Ψ ( s ) ,
where Ψ is a Bernstein function ([20]) with the representation
Ψ ( s ) = b s + 0 + ( 1 e s τ ) Φ ( d τ ) ,
where b 0 is the drift coefficient and Φ is the Lévy measure, such that either b > 0 , or Φ ( R + ) = + , or both.
Because the process D is strictly increasing, it possesses an inverse function
E ( t ) = inf { r > 0 : D ( r ) > t } .
Now we consider a Poisson process N ( t ) with the intensity λ . It is known that N ( E ( t ) ) is a renewal process with the waiting times J n and
P [ J n > t ] = E e λ E ( t ) .
As shown in [10], if the restriction Ψ = Ψ ( p ) , p > 0 of the Laplace exponent Ψ = Ψ ( s ) to the real semi-axes R + is a complete Bernstein function that satisfies the conditions (K3) and (K4), then the right-hand side of Formula (38) can be interpreted as the solution to the Cauchy problem (30), (31). According to Theorem 2, it is continuous on [ 0 , + ) and completely monotone.

3.2. Fractional Growth Equation

In this subsection, some of the results derived in [13] are shortly addressed. We consider the fractional growth equation with the GFD of Caputo type
( D k C u ) ( t ) = λ u ( t ) , λ > 0 , t > 0
subject to the initial condition
u ( 0 ) = 1 .
In the case of the conventional growth equation, the solution is u λ ( t ) = exp ( λ t ) . The problem (39), (40) with the Caputo fractional derivative (7) is solved by the function u λ ( t ) = E α , 1 ( λ t α ) that is known to be of exponential growth as t + ([26]). It turns out that the solution to the problem (39), (40) with the GFD with a kernel k K (k fulfills the conditions (K1)–(K4)) is also of exponential growth as t + .
For formulation of the corresponding result, the notation
Φ ( p ) = p k ˜ ( p )
is introduced, k ˜ being the Laplace transform of k. Because k ˜ is a Stieltjes function, Φ is a Bernstein function and Φ is completely monotone ([20]). The made assumptions ensure that Φ is strictly monotone and thus for every λ > 0 there exists a unique p 0 = p 0 ( λ ) such that Φ ( p 0 ) = λ . Then we have the following result regarding the asymptotic of the unique solution u λ of the Cauchy problem (39), (40) as t + .
Theorem 4
([13]). Let k K and the condition
1 + d p p Φ ( p ) < +
hold true. Then the solution u λ of the Cauchy problem (39), (40) has the following exponential asymptotic behavior:
u λ ( t ) = λ Φ ( p 0 ( λ ) ) p 0 ( λ ) e p 0 ( λ ) t + o e p 0 ( λ ) t , t + .
In the case of the Cauchy problem (39), (40) with the Caputo derivative (7), the function Φ takes the form Φ ( p ) = p α and p 0 ( λ ) = λ 1 / α . The condition (41) is evidently fulfilled and Formula (42) takes the form
u λ ( t ) = 1 α e λ 1 α t + o e λ 1 α t ,
that corresponds to the main term of the asymptotic of the Mittag-Leffler function u λ ( t ) = E α , 1 ( λ t α ) as t + (see e.g., [26]).
Another important particular case is the Cauchy problem (39), (40) with the distributed order derivative (the GFD with the kernel (24)). It turns out that under some standard assumptions Theorem 4 is applicable also in this case (see [13] for details).

3.3. The Cauchy Problem for a Nonlinear Fractional ODE

Following [14], in this subsection, we address a nonlinear fractional differential equation with the GFD of the Caputo type
( D k C u ) ( t ) = f ( t , u ( t ) ) , t > 0
subject to the initial condition
u ( 0 ) = u 0 R .
In what follows, we again assume that the kernel k of D k C defined by (21) belongs to the class K , i.e., that the conditions (K1)–(K4) are fulfilled.
For derivation of results regarding existence and uniqueness of the solution to the Cauchy problem (43), (44), we first transform it into an integral equation by applying the GFI (26). Let L > 0 and f : [ 0 , L ] × R R be a continuous function. Then the Cauchy problem (43), (44) is equivalent to the integral equation
u ( t ) = u 0 + 0 t κ ( t τ ) f ( τ , u ( τ ) ) d τ , t > 0 ,
where the kernel function κ is determined by the relation (25).
The sufficient conditions for existence of the local and global solutions to the integral Equation (45) and thus to the Cauchy problem (43), (44) are provided in the following two theorems.
Theorem 5
([14]). Let L , Q > 0 , f be a continuous function on the closed domain G = { ( t , τ ) : 0 t L , | τ u 0 | Q } , and l ( 0 , L ] satisfy the inequality
max ( t , τ ) G | f ( t , τ ) | 0 l κ ( τ ) d τ Q .
Then the Cauchy problem (43), (44) has a solution absolutely continuous on the interval [ 0 , l ] .
Theorem 6
([14]). Let L > 0 and f : [ 0 , L ] × R R be a continuous function that satisfies the inequality
| f ( t , τ ) | b 0 + b 1 | τ | p , 0 t L , τ R
with b 0 , b 1 > 0 and 0 < p 1 .
Then the Cauchy problem (43), (44) has a solution absolutely continuous on the interval [ 0 , L ] .
As a corollary from Theorem 6, we get the following result: Let f : [ 0 , + ) × R R be a continuous function that satisfies the inequality
| f ( t , τ ) | b 0 + b 1 | τ | p , 0 t < + , τ R
with b 0 , b 1 > 0 and 0 < p 1 . Then the Cauchy problem (43), (44) has a global solution absolutely continuous on [ 0 , + ) .
As to uniqueness of the solution, it was proved under some stronger conditions compared to the ones required for its existence.
Theorem 7
([14]). Let L > 0 and f : [ 0 , L ] × R R be a continuous function that satisfies the Lipschitz condition with respect to its second variable
| f ( t , s ) f ( t , τ ) | b | s τ | , 0 t L , s , τ R , b > 0 .
Then the Cauchy problem (43), (44) has a unique absolutely continuous solution on the interval [ 0 , L ] .
As before, the result formulated in Theorem 7 can be extended to the case of the infinite interval [ 0 , + ) : Let f : [ 0 , + ) × R R be a continuous function that satisfies the Lipschitz condition
| f ( t , s ) f ( t , τ ) | b | s τ | , 0 t < + , s , τ R , b > 0 .
Then the Cauchy problem (43), (44) has a unique absolutely continuous solution on the interval [ 0 , + ) .
To address a continuous dependence of solutions to the Cauchy problem (43), (44) on the problem data, in [14], a Gronwall-type inequality was derived. It is important by itself and we formulate it below.
Lemma 1
([14]). Let l , v 0 , b > 0 and v C [ 0 , l ] . If the inequality
v ( t ) v 0 + b 0 t κ ( t τ ) v ( τ ) d τ
holds true for t [ 0 , l ] , then
v ( t ) u ( t ) , t [ 0 , l ] ,
where u = u ( t ) is the solution to the fractional relaxation equation
( D k C u ) ( t ) = b u ( t )
subject to the initial condition
u ( 0 ) = v 0 .
Based on Lemma 1, in [14], the continuous dependence of the solution to the Cauchy problem (43), (44) on the problem data was proved.
Theorem 8
([14]). Let the conditions of Theorem 7 be fulfilled, u be the solution to the Cauchy problem (43), (44), and u k be the solution to the Cauchy problem (43), (44) with the initial condition u ( 0 ) = u 0 k . Then
u u k C [ 0 , L ] 0 a s u 0 k u 0 .
Theorem 9
([14]). Let the conditions of Theorem 7 be fulfilled for the functions f and g, u be the solution to the Cauchy problem (43), (44), and v be the solution to the Cauchy problem (43), (44) with the right-hand side g and with the same initial condition u 0 . Then
u v C [ 0 , L ] 0 a s max ( t , τ ) [ 0 , L ] × [ Ω , Ω ] | f ( t , τ ) g ( t , τ ) | 0 ,
where Ω is the upper bound of the supremum norm of the solution u.

4. Time-Fractional PDEs with the GFD

In this section, we present some important results concerning the direct and inverse problems for the time-fractional PDEs with the GFD of the Caputo type.

4.1. Cauchy Problem for the Time-Fractional Diffusion Equation

Following [10], in this subsection we address the properties of solutions to the Cauchy problem for the general time-fractional diffusion equation in the form
( D k C u ( x , · ) ) ( t ) = Δ u ( x , t ) t > 0 , x R n , u ( x , 0 ) = u 0 ( x ) ,
where u 0 is a bounded globally Hölder continuous function on R n .
In [10], solutions to the Cauchy problem (51) were understood in the following sense: Applying formally the Laplace transform in the variable t to the equation in (51), we arrive at the equation
p k ˜ ( p ) u ˜ ( x , p ) k ˜ ( p ) u 0 ( x ) = Δ u ˜ ( x , p ) , p > 0 , x R n
for the Laplace transform u ˜ ( x , p ) of a solution to the Cauchy problem (51). A bounded function u = u ( x , t ) is called an LT-solution of (51), if u is continuous in t on [ 0 , + ) uniformly with respect to x R n , satisfies the initial condition u ( x , 0 ) = u 0 ( x ) , while its Laplace transform u ˜ ( x , p ) is twice continuously differentiable in x, for each p > 0 , and satisfies Equation (52).
Theorem 10
([10]). Let the kernel k of the GFD of the Caputo type belong to the class K .
Then there exist a non-negative function Z = Z ( x , t ) , t > 0 , x R n , x 0 , locally integrable in t and infinitely differentiable in x 0 that satisfies the relation
R n Z ( x , t ) d x = 1 , t > 0 ,
and for any bounded globally Hölder continuous u 0 , the function
u ( x , t ) = R n Z ( x ζ , t ) u 0 ( ζ ) d ζ
is an LT-solution to the Cauchy problem (51).
The function Z = Z ( x , t ) is what is usually called the fundamental solution to the Cauchy problem (51), i.e., the one that formally corresponds to the initial condition u ( x , 0 ) = u 0 ( x ) = δ ( x ) , δ being the Dirac delta function. As stated in Theorem 10, the fundamental solution to the Cauchy problem (51) can be interpreted as a probability density function for each t > 0 and thus the time-fractional diffusion Equation (51) with the GFD could be potentially useful for modeling of the anomalous diffusion processes.
The explicit form of the fundamental solution is as follows ([10]):
Z ( x , t ) = 0 + ( 4 π s ) n / 2 e | x | 2 4 s G ( s , t ) d s , x 0 ,
where
G ( s , t ) = 0 t k ( t τ ) μ s ( d τ )
with the probability measure μ s ( d τ ) that satisfies the relation
e s p k ˜ ( p ) = 0 + e p τ μ s ( d τ ) .
Under the conditions of Theorem 10, the measure μ s ( d τ ) always exists because the function p e s p k ˜ ( p ) is completely monotone on R + for each s 0 .
For the validity of Theorem 10, the condition k K is essential. However, it is worth mentioning that for the well-posedness of the Cauchy problem for the equations with the operators of type D k C much weaker conditions than (K1)–(K4) are sufficient (see e.g., [27]). Here we formulate the uniqueness result for the Cauchy problem (51) in the class of the LT-solutions proved in [10].
Theorem 11
([10]). Let the kernel k of the GFD of the Caputo type be non-negative, locally integrable, nonzero on a set of positive measure, and its Laplace transform k ˜ = k ˜ ( p ) exist for all p > 0 .
If u = u ( x , t ) is a polynomially bounded LT-solution to the Cauchy problem (51) with u 0 ( x ) 0 , then u ( x , t ) 0 .
In the definition of the LT-solutions, their boundedness was required. However, this definition makes sense also for the polynomially bounded solutions, i.e., such solutions u = u ( x , t ) that satisfy the inequality | u ( x , t ) | P u ( | x | ) , where P u are some polynomials independent on t.

4.2. Initial-Boundary-Value Problems for the Time-Fractional Diffusion Equation

In this subsection, we present some results regarding the initial-boundary-value problems for the time-fractional diffusion equation with the GFD of the Caputo type in the form
( D k C u ( x , · ) ) ( t ) = D 2 ( u ) + D 1 ( u ) q ( x ) u ( x , t ) + F ( x , t ) , ( x , t ) Ω × ( 0 , T ] ,
subject to the initial condition
u ( x , t ) | t = 0 = u 0 ( x ) , x Ω ¯
and the boundary condition
u ( x , t ) | ( x , t ) Ω × ( 0 , T ] = v ( x , t ) , ( x , t ) Ω × ( 0 , T ] .
In Equations (55)–(57), Ω is a bounded open domain in R n with a smooth boundary Ω , q C ( Ω ¯ ) , q ( x ) 0 for x Ω ¯ , and
D 1 ( u ) = i = 1 n b i ( x ) u x i , D 2 ( u ) = i , j = 1 n a i j ( x ) 2 u x i x j .
Moreover we assume that D 2 is a uniformly elliptic differential operator.
A function u S ( Ω , T ) is called a strong solution to the initial-boundary-value problem (55)–(57) if it satisfies both Equation (55) and the initial and boundary conditions (56) and (57), respectively. By S ( Ω , T ) , we denoted the space of functions u = u ( x , t ) , ( x , t ) Ω ¯ × [ 0 , T ] that satisfy the inclusions u C ( Ω ¯ × [ 0 , T ] ) , u ( · , t ) C 2 ( Ω ) for any t > 0 , and t u ( x , · ) C ( 0 , T ] L 1 ( 0 , T ) for any x Ω .
First we discuss a maximum principle for the time-fractional diffusion Equation (55). For its validity, we assume that the following conditions for the kernel k of D k C are satisfied:
(LY1)
k C 1 ( R + ) L l o c 1 ( R + ) ,
(LY2)
k ( τ ) > 0 and k ( τ ) < 0 for τ > 0 ,
(LY3)
k ( τ ) = o ( τ 1 ) , τ 0 .
Let us note that the Kochubei’s conditions (K1)–(K4) are not needed for validity of the maximum principle for the general diffusion Equation (55). However, if the condition (K3) holds true, then it follows from the Feller-Karamata Tauberian theorem for the Laplace transform ([21]) that the condition (LY3) is also satisfied.
The maximum principle for the general diffusion Equation (55) is based on an appropriate estimate of the GFD of a function f at its maximum point. It is given in the following theorem.
Theorem 12
([28]). Let the conditions (LY1)–(LY3) be fulfilled, a function f C [ 0 , T ] attain its maximum over the interval [ 0 , T ] at the point t 0 , t 0 ( 0 , T ] , and f C ( 0 , T ] L 1 ( 0 , T ) .
Then the inequality
( D k C f ) ( t 0 ) k ( t 0 ) ( f ( t 0 ) f ( 0 ) ) 0
holds true.
In the case of the Caputo fractional derivative, the inequality (59) takes the known form
( D C α f ) ( t 0 ) t 0 α Γ ( 1 α ) ( f ( t 0 ) f ( 0 ) ) 0 .
In what follows, we use the notation
P k ( u ) : = ( D k C u ) ( t ) D 2 ( u ) D 1 ( u ) + q ( x ) u ( x , t ) .
Theorem 13
([28]). Let the conditions (LY1)–(LY3) be fulfilled and a function u S ( Ω , T ) satisfy the inequality
P k ( u ) 0 , ( x , t ) Ω × ( 0 , T ] .
Then the following maximum principle holds true:
max ( x , t ) Ω ¯ × [ 0 , T ] u ( x , t ) max { max x Ω ¯ u ( x , 0 ) , max ( x , t ) Ω × [ 0 , T ] u ( x , t ) , 0 } .
The maximum principle formulated in Theorem 13 can be applied, among other things, for derivation of some a priori estimates for the strong solutions of the initial-boundary-value problem (55)–(57).
Theorem 14
([28]). Let the conditions (K1)–(K4) and (LY1)–(LY3) be fulfilled and u be a strong solution to the initial-boundary-value problem (55)–(57).
Then
u C ( Ω ¯ × [ 0 , T ] ) max { M 0 , M 1 } + M f ( T ) ,
where
M 0 = u 0 C ( Ω ¯ ) , M 1 = v C ( Ω × [ 0 , T ] ) , M = F C ( Ω × [ 0 , T ] ) ,
and
f ( t ) = 0 t κ ( τ ) d τ ,
where the function κ is the kernel of the GFI defined by (25).
The uniqueness of the strong solution to the initial-boundary-value problem (55)–(57) and its continuous dependence on problem data easily follow from the solution norm estimate (64).
Theorem 15
([28]). The initial-boundary-value problem (55)–(57) possesses at most one strong solution.
This solution—if it exists—continuously depends on the problem data in the sense that if u and u ˜ are strong solutions to the problems with the sources functions F and F ˜ and the initial and boundary conditions u 0 and u ˜ 0 and v and v ˜ , respectively, and
F F ˜ C ( Ω ¯ × [ 0 , T ] ) ϵ ,
u 0 u ˜ 0 C ( Ω ¯ ) ϵ 0 , v v ˜ C ( Ω × [ 0 , T ] ) ϵ 1 ,
then the norm estimate
u u ˜ C ( Ω ¯ × [ 0 , T ] ) max { ϵ 0 , ϵ 1 } + ϵ f ( T )
holds true, where the function f is defined by (66).
In the rest of this subsection, we address uniqueness and existence of a weak solution to the initial-boundary-value problem (55)–(57) defined in the sense of Vladimirov [16]: We call u C ( Ω ¯ × [ 0 , T ] ) a weak solution to the initial-boundary-value problem (55)–(57) in the sense of Vladimirov, if there exist F k C ( Ω ¯ × [ 0 , T ] ) , u 0 k C ( Ω ¯ ) and v k C ( Ω × [ 0 , T ] ) , k = 1 , 2 , satisfying (V1) and (V2) below such that
u k u C ( Ω ¯ × [ 0 , T ] ) 0 as k + .
(V1)
There exist the functions F, u 0 , and v, such that
F k F C ( Ω ¯ × [ 0 , T ] ) 0 as k + ,
u 0 k u 0 C ( Ω ¯ ) 0 as k + ,
v k v C ( Ω × [ 0 , T ] ) 0 as k + .
(V2)
For each k = 1 , 2 , there exists a strong solution u k = u k ( x , t ) to the general time-fractional diffusion equation
( D k C u k ( x , · ) ) ( t ) = D 2 ( u k ) + D 1 ( u k ) q ( x ) u k ( x , t ) + F k ( x , t ) , ( x , t ) Ω × ( 0 , T ] .
subject to the initial condition
u k | t = 0 = u 0 k ( x ) , x Ω ¯
and the boundary condition
u k | Ω × ( 0 , T ] = v k ( x , t ) , ( x , t ) Ω × ( 0 , T ] .
In [28], the correctness of the definition of a weak solution was shown. A weak solution to the problem (55)–(57) in the sense of Vladimirov is a continuous function, not a distribution. However, the weak solutions are not required to be smooth.
Any strong solution to the problem (55)–(57) is evidently also its weak solution. If the problem (55)–(57) possesses a weak solution, then the functions F, u 0 and v from the problem formulation have to belong to the spaces C ( Ω ¯ × [ 0 , T ] ) , C ( Ω ¯ ) and C ( Ω × [ 0 , T ] ) , respectively, as the limits of the sequences of the continuous functions in the uniform norm.
The estimate (64) for the strong solutions holds true also for the weak solutions. To show this, we just let k + in the inequality
u k C ( Ω ¯ T ) max { M 0 k , M 1 k } + M k f ( T ) , k = 1 , 2
with
M 0 k : = u 0 k C ( Ω ¯ ) , M 1 k : = v k C ( Ω × [ 0 , T ] ) , M k : = F k C ( Ω ¯ × [ 0 , T ] ) .
The estimate (64) for the weak solutions is employed to prove the following uniqueness result.
Theorem 16
([28]). The initial-boundary-value problem (55)–(57) possesses at most one weak solution. The weak solution—if it exists—continuously depends on the data given in the problem in the sense of the estimate (67).
In the rest of this subsection, we address the question of existence of a weak solution to the initial-boundary-value problem (55)–(57) in the case of the homogeneous Equation (55) without the first order spatial differential operator D 1 subject to the initial condition (56) and the homogeneous boundary condition (57), i.e., we consider the initial-boundary-value problem
( D k C u ( x , · ) ) ( t ) = D 2 ( u ) q ( x ) u ( x , t ) , ( x , t ) Ω × ( 0 , T ] ,
u ( x , t ) | t = 0 = u 0 ( x ) , x Ω ¯ ,
u ( x , t ) | ( x , t ) Ω × ( 0 , T ] = 0 , ( x , t ) Ω × ( 0 , T ]
under the same conditions on the coefficients of the operator D 2 and the function q that we assumed at the beginning of the subsection. Moreover, we also assume that the kernel k of D k C is from K , i.e., the conditions (K1)–(K4) are satisfied.
First, a formal solution to the initial-boundary-value problem (76)–(78) is constructed in form of the Fourier series
u ( x , t ) = k = 1 ( u 0 , X k ) U k ( t ) X k ( x ) ,
where X k , k = 1 , 2 , are the eigenfunctions corresponding to the eigenvalues λ k of the eigenvalue problem
L ( X ( x ) ) = λ X ( x ) , x Ω ,
X ( x ) | x Ω = 0
for the operator L , L ( x ) = D 2 ( X ) + q ( x ) X ( x ) . Because of the conditions posed on the operator D 2 and the function q, the differential operator L is positive definite and self-adjoint. Thus the eigenvalue problem (80)–(81) has a countable number of the positive eigenvalues 0 < λ 1 λ 2 with the finite multiplicities and—if the boundary Ω of Ω is smooth—any function f M L can be represented through its Fourier series in the form
f ( x ) = k = 1 ( f , X k ) X k ( x ) ,
where X k M L are the eigenfunctions corresponding to the eigenvalues λ k :
L ( X k ) = λ k X k , k = 1 , 2 , .
By M L , the space of the functions f that satisfy the boundary condition (81) and the inclusions f C 1 ( Ω ¯ ) C 2 ( Ω ) , L ( f ) L 2 ( Ω ) is denoted.
As to the functions U k = U k ( t ) , they are solutions to the fractional relaxation equations
( D k C U k ) ( t ) = λ k U k ( t ) , t > 0 , k = 1 , 2 ,
subject to the initial conditions
U k ( 0 ) = 1 , k = 1 , 2 , .
According to Theorem 2 from Section 3, for any λ = λ k > 0 , k = 1 , 2 , this initial-value problem has a unique solution U k = U k ( t ) that belongs to the class C ( R + ) and is a completely monotone function. In particular, any U k is non-negative and non-increasing and thus the inequalities
0 U k ( t ) U k ( 0 ) = 1
hold true. Let us mention that in the case of the single-term time-fractional diffusion equation with the Caputo fractional derivative ( k ( τ ) = τ α Γ ( 1 α ) , 0 < α < 1 ), the solution to the initial-value problem (84), (85) with λ = λ k , k = 1 , 2 , has the form ([29])
U k ( t ) = E α , 1 ( λ k t α ) .
Under some standard assumptions, the formal solution (79) is a weak solution to the initial-boundary-value problem (76)–(78) in the sense of Vladimirov.
Theorem 17
([15]). Let the function u 0 in the initial condition (77) be from the space M L . Then the formal solution (79) of the problem (76)–(78) is its weak solution in the sense of Vladimirov.
For a survey of other results regarding the maximum principles for the time-fractional PDEs of different types see the recent publication [30].

4.3. Inverse Problems Involving GFD

The starting point of this subsection is a reconstruction problem for a function based on its values and the values of its GFD in a neighborhood of the final time ([19]):
  • IP1. Let 0 < t 0 < T < + . Given ϕ , g : ( t 0 , T ) R , find a function u : ( 0 , T ) R such that
    u ( t ) | ( t 0 , T ) = ϕ ( t ) , and ( D k C u ) ( t ) | ( t 0 , T ) = g ( t ) .
The inverse problems of type IP1 are potentially useful for applications. For instance, in the framework of the Scott-Blair model of viscoelasticity, the stress is proportional to a time-fractional derivative of the strain ([23]). In this context, the IP1 means a reconstruction of the strain history based on the measurements of strain and stress starting from a certain time t 0 .
In [19], a uniqueness result for the IP1 was proved under the following conditions on the kernel k of D k C :
(KJ1)
μ R : 0 + e μ t | k ( t ) | d t < + ,
(KJ2)
k is real analytic on R + ,
(KJ3)
the Laplace transform k ˜ of k cannot be meromorphically extended to the whole complex plane C .
Theorem 18
([19]). Let the kernel k of D k C fulfill the conditions (KJ1)–(KJ3). Then the following uniqueness results for the IP1 hold true:
(i) 
If u L 1 ( 0 , T ) , k u W 1 , 1 ( 0 , T ) , and u ( t ) | ( t 0 , T ) = ( D k C u ) ( t ) | ( t 0 , T ) = 0 , then u = 0 ,
(ii) 
If u W 1 , 1 ( 0 , T ) and u ( t ) | ( t 0 , T ) = ( D k C u ) ( t ) | ( t 0 , T ) = 0 , then u = 0 .
The results formulated in Theorem 18 were employed in [19] for studying uniqueness of solution to the following source reconstruction problem for the fractional PDEs with the GFD of Caputo type:
  • IP2. Let 0 < t 0 < T < + and Ω R n . Given ϕ , Φ : Ω × ( t 0 , T ) R , find the functions u , F : Ω × ( 0 , T ) R such that they fulfill the equation
    ( D k C B ( u ) ) ( t ) + D l ( u ) A ( u ) = F ( x , t ) , x Ω , t ( 0 , T )
    and the relations
    u ( x , t ) | Ω × ( t 0 , T ) = ϕ , and F | Ω × ( t 0 , T ) = Φ
    hold true.
In the formulation of IP2, D l = j = 1 l q j j t j is a differential operator of order l with respect to the time variable t and with q j R and A and B are some operators that act with respect to the spatial variable x. Moreover, we assume that D ( A ) C ( Ω ) C ( Ω ) , D ( B ) C ( Ω ) C ( Ω ) and the operator B is invertible.
In particular, the time-fractional PDE (89) includes the time-fractional diffusion Equations (51) and (55) that were considered in the previous subsections of this section.
As shown in [19], IP2 can be reduced to IP1. Indeed, let the pair of functions ( u , F ) solve the IP2. The Equation (89) restricted to Ω × ( t 0 , T ) has the form ( D k C B ( u ) ) ( t ) + D l ( ϕ ) A ( ϕ ) = Φ ( x , t ) and thus the function B u is a solution to the following inverse problem of IP1 type:
B u | Ω × ( t 0 , T ) = B ϕ , and D k C B u | Ω × ( t 0 , T ) = g ,
where
g ( x , t ) = Φ ( x , t ) + A ( ϕ ) D l ( ϕ ) , x Ω , t ( t 0 , T ) .
The solution ( u , F ) of IP2 can be explicitly expressed in terms of the solution B u of the IP1 formulated above as follows: u = B 1 B u , F = D k C B u + D l ( u ) A ( u ) . Accordingly, a uniqueness result for the IP2 immediately follows from Theorem 18 (see [19] for details).
It is worth mentioning that the inverse problems IP1 and IP2 are severely ill-posed ([19]) and thus appropriate regularization methods are needed for their numerical treatment.
Next, we consider the following evolutionary integral equation:
u ( x , t ) = 0 t κ ( t τ ) Δ u ( x , τ ) d τ + f ( x , t ) , x R n , t 0 .
Please note that the Cauchy problem (compare to (51))
( D k C u ( x , · ) ) ( t ) = Δ u ( x , t ) + F ( x , t ) t > 0 , x R n , u ( x , 0 ) = u 0 ( x )
can be reduced to an evolutionary integral equation of type (92) by applying the GFI (26) to both sides of this equation and by using Formula (28) from Theorem 1:
u ( x , t ) = 0 t κ ( t τ ) Δ u ( x , τ ) d τ + ( I k C F ( x , · ) ) ( t ) + u 0 ( x ) , x R n , t 0 ,
where the kernel κ is connected with the kernel k of D k C by means of the relation (25).
In [17], an important inverse problem of kernel identification in the boundary value problems for an equation associated with the evolutionary integral Equation (92) was addressed. Let Ω R n be a bounded domain with sufficiently smooth boundary Ω . The direct boundary value problem formulated in [17] is as follows:
u ( x , t ) = 0 t κ ( t τ ) Δ u ( x , τ ) d τ + f ( x , t ) , x Ω , t [ 0 , T ] , B u ( x , t ) = 0 , x Ω , t [ 0 , T ] ,
where B is a boundary operator of the Dirichlet, Neumann, or Robin type, respectively:
B ν ( x ) = ν ( x ) , B ν ( x ) = n ( x ) · ν ( x ) , B ν ( x ) = n ( x ) · ν ( x ) + θ ν ( x ) , θ 0 ,
n ( x ) being the unit outer normal of Ω at the point x Ω .
The inverse problem addressed in [17] is formulated via the so called observation functional Φ that maps the functions defined on Ω ¯ onto R . Usually, the functional Φ is defined in one of the following ways:
Φ [ ν ] = ν ( x 0 ) , Φ [ ν ] = n ( x 0 ) · ν ( x 0 ) , Φ [ ν ] = Ω μ ( x ) ν ( x ) d x ,
where x 0 Ω ¯ and μ : Ω R are given. In the case x 0 Ω , the observation functional has to be different from the boundary operator, i.e., Φ [ ν ] B ν ( x 0 ) .
The inverse problem considered in [17] is as follows:
  • IP3. Given h : ( 0 , T ) R find a kernel κ such that the solution u of the boundary value problem (95) satisfies the condition
    Φ [ u ( t , · ) ] = h ( t ) , t ( 0 , T ) .
In [17], existence, uniqueness, and stability of solutions to the IP3 were studied for a certain class of kernels (see [17] for details).
Finally, we mention that in [18] two other inverse problems for a time-fractional PDE with the GFD of Caputo type were addressed. Let Ω R n be a bounded domain with the boundary Ω . The direct initial-boundary-value problem is formulated as follows:
d d t ( k ( U Φ ) ) ( x , t ) = L x U ( x , t ) + H ( x , t ) , x Ω , t ( 0 , T ) , U ( x , 0 ) = Φ ( x ) , x Ω , B ( U b ) ( x , t ) = 0 , x Ω , t ( 0 , T ) ,
where Φ and b are given functions, the operator L x is a linear second order differential operator with respect to the variable x in the form
L x U ( x , t ) = i , j = 1 n a i j ( x ) 2 U x i x j + j = 1 n a j ( x ) U x j + r ( x ) U ( x , t ) ,
and B is a boundary operator of the Dirichlet or Neumann type, respectively:
B v ( x ) = v ( x ) or B v ( x ) = ω ( x ) · v ( x )
with ω · n > 0 , n being the unite outer normal of Ω at the point x Ω .
The inverse problems considered in [18] are formulated in terms of the given observation function Φ at the final time T in the form
U ( x , T ) = Φ ( x ) , x Ω .
In [18], the following inverse problems were addressed:
  • IP4. (inverse source problem). Let
    H ( x , t ) = g ( x , t ) f ( x ) + h 0 ( x , t ) ,
    where the components g f and h 0 may correspond to different sources or sinks. The factor f is unknown and has to be reconstructed by means of the data given by the observation function Φ from (98). The inverse problem consists in determination of a pair of functions ( f , U ) that satisfies (97), (98), and (99).
Another inverse problem considered in [18] is determination of the coefficient r in the operator L x U :
  • IP5. Determine a pair of functions ( r , U ) that satisfies (97) and (98).
In [18], existence, uniqueness, and stability of solutions to IP4 and to IP5 were shown under some additional conditions posed on the problem data (see [18] for details).
For the surveys of the recent results concerning the inverse problems for the fractional PDEs including different kinds of the conventional fractional derivatives we refer the readers to [31,32,33].
Finally we mention that the theory of the fractional PDEs with the GFDs is still far away from being completed. In particular, the regularity of their solutions is not yet investigated in detail. Another interesting problem for further research would be to address the abstract fractional evolution equation in the form
( D k C u ) ( t ) = A u ( t )
subject to the initial condition u ( 0 ) = x . In (100), A stands for a linear closed unbounded operator densely defined in a Banach space X and the initial condition x belongs to the space X. This problem will be considered elsewhere.

Author Contributions

Y.L. and M.Y. have contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

M.Y. was supported by the Grant-in-Aid for Scientific Research (S) 15H05740 of the Japan Society for the Promotion of Science and by the National Natural Science Foundation of China (No. 11771270, 91730303). This work was prepared under the support of the “RUDN University Program 5-100”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alt, H.W. Lineare Funktionalanalysis; Springer: Berlin, Germany, 2006. [Google Scholar]
  2. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach: New York, NY, USA, 1993. [Google Scholar]
  3. Cartwright, D.I.; McMullen, J.R. A note on the fractional calculus. Proc. Edinb. Math. Soc. 1978, 21, 79–80. [Google Scholar] [CrossRef] [Green Version]
  4. Hilfer, R.; Luchko, Y. Desiderata for Fractional Derivatives and Integrals. Mathematics 2019, 7, 149. [Google Scholar] [CrossRef] [Green Version]
  5. Diethelm, K.; Garrappa, R.; Giusti, A.; Stynes, M. Why fractional derivatives with nonsingular kernels should not be used. Fract. Calc. Appl. Anal. 2020, 23, 610–634. [Google Scholar] [CrossRef]
  6. Luchko, Y. Fractional derivatives and the fundamental theorem of fractional calculus. Fract. Calc. Appl. Anal. 2020, 23, 939–966. [Google Scholar] [CrossRef]
  7. Clément, P. On abstract Volterra equations in Banach spaces with completely positive kernels. In Lecture Notes in Math; Kappel, F., Schappacher, W., Eds.; Springer: Berlin, Germany, 1984; Volume 1076, pp. 32–40. [Google Scholar]
  8. Zacher, R. Boundedness of weak solutions to evolutionary partial integro-differential equations with discontinuous coefficients. J. Math. Anal. Appl. 2008, 348, 137–149. [Google Scholar] [CrossRef] [Green Version]
  9. Zacher, R. Weak solutions of abstract evolutionary integro-differential equations in Hilbert spaces. Funkcial. Ekvac. 2009, 52, 1–18. [Google Scholar] [CrossRef] [Green Version]
  10. Kochubei, A.N. General fractional calculus, evolution equations, and renewal processes. Integr. Equa. Oper. Theory 2011, 71, 583–600. [Google Scholar] [CrossRef] [Green Version]
  11. Kochubei, A.N. General fractional calculus. In Handbook of Fractional Calculus with Applications; Volume 1: Basic Theory; Kochubei, A., Luchko, Y., Eds.; Walter de Gruyter: Berlin, Germany; Boston, MA, USA, 2019; pp. 111–126. [Google Scholar]
  12. Kochubei, A.N. Equations with general fractional time derivatives. Cauchy problem. In Handbook of Fractional Calculus with Applications; Volume 2: Fractional Differential Equations; Kochubei, A., Luchko, Y., Eds.; Walter de Gruyter: Berlin, Germany; Boston, MA, USA, 2019; pp. 223–234. [Google Scholar]
  13. Kochubei, A.N.; Kondratiev, Y. Growth Equation of the General Fractional Calculus. Mathematics 2019, 7, 615. [Google Scholar] [CrossRef] [Green Version]
  14. Sin, C.-S. Well-posedness of general Caputo-type fractional differential equations. Fract. Calc. Appl. Anal. 2018, 21, 819–832. [Google Scholar] [CrossRef]
  15. Luchko, Y.; Yamamoto, M. General time-fractional diffusion equation: Some uniqueness and existence results for the initial-boundary-value problems. Fract. Calc. Appl. Anal. 2016, 19, 675–695. [Google Scholar] [CrossRef]
  16. Vladimirov, V.S. Equations of Mathematical Physics; Nauka: Moscow, Russia, 1971. [Google Scholar]
  17. Janno, J.; Kasemets, K. Identification of a kernel in an evolutionary integral equation occurring in subdiffusion. J. Inverse Ill-Posed Probl. 2017, 25, 777–798. [Google Scholar] [CrossRef]
  18. Kinash, N.; Janno, J. Inverse problems for a generalized subdiffusion equation with final overdetermination. Math. Model. Anal. 2019, 24, 236–262. [Google Scholar]
  19. Kinash, N.; Janno, J. An Inverse Problem for a Generalized Fractional Derivative with an Application in Reconstruction of Time- and Space-Dependent Sources in Fractional Diffusion and Wave Equations. Mathematics 2019, 7, 1138. [Google Scholar] [CrossRef] [Green Version]
  20. Schilling, R.L.; Song, R.; Vondracek, Z. Bernstein Functions. Theory and Application; Walter de Gruyter: Berlin, Germany, 2010. [Google Scholar]
  21. Feller, W. An Introduction to Probability Theory and Its Applications; Wiley: New York, NY, USA, 1966; Volume 2. [Google Scholar]
  22. Hanyga, A. A comment on a controversial issue: A generalized fractional derivative cannot have a regular kernel. Fract. Calc. Anal. Appl. 2020, 23, 211–223. [Google Scholar] [CrossRef] [Green Version]
  23. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity; Imperial College Press: London, UK, 2010. [Google Scholar]
  24. Luchko, Y. On Complete Monotonicity of Solution to the Fractional Relaxation Equation with the nth Level Fractional Derivative. Mathematics 2020, 8, 1561. [Google Scholar] [CrossRef]
  25. Bertoin, J. Lévy Processes; Cambridge University Press: Cambridge, UK, 1996. [Google Scholar]
  26. Gorenflo, R.; Mainardi, F.; Rogosin, S. Mittag-Leffler function: Properties and applications. In Handbook of Fractional Calculus with Applications; Volume 1: Basic Theory; Kochubei, A., Luchko, Y., Eds.; Walter de Gruyter: Berlin, Germany; Boston, MA, USA, 2019; pp. 269–298. [Google Scholar]
  27. Gripenberg, G. Volterra integro-differential equations with accretive nonlinearity. J. Differ. Equ. 1985, 60, 57–79. [Google Scholar] [CrossRef] [Green Version]
  28. Luchko, Y.; Yamamoto, M. On the maximum principle for a time-fractional diffusion equation. Fract. Calc. Appl. Anal. 2017, 20, 1131–1145. [Google Scholar] [CrossRef] [Green Version]
  29. Luchko, Y.; Gorenflo, R. An operational method for solving fractional differential equations with the Caputo derivatives. Acta Math. Vietnam. 1999, 24, 207–233. [Google Scholar]
  30. Luchko, Y.; Yamamoto, M. Maximum principle for the time-fractional PDEs. In Handbook of Fractional Calculus with Applications; Volume 2: Fractional Differential Equations; Kochubei, A., Luchko, Y., Eds.; Walter de Gruyter: Berlin, Germany; Boston, MA, USA, 2019; pp. 299–326. [Google Scholar]
  31. Liu, Y.; Li, Z.; Yamamoto, M. Inverse problems of determining sources of the fractional partial differential equations. In Handbook of Fractional Calculus with Applications; Volume 2: Fractional Differential Equations; Kochubei, A., Luchko, Y., Eds.; Walter de Gruyter: Berlin, Germany; Boston, MA, USA, 2019; pp. 411–430. [Google Scholar]
  32. Li, Z.; Liu, Y.; Yamamoto, M. Inverse problems of determining parameters of the fractional partial differential equations. In Handbook of Fractional Calculus with Applications; Volume 2: Fractional Differential Equations; Kochubei, A., Luchko, Y., Eds.; Walter de Gruyter: Berlin, Germany; Boston, MA, USA, 2019; pp. 431–442. [Google Scholar]
  33. Li, Z.; Yamamoto, M. Inverse problems of determining coefficients of the fractional partial differential equations. In Handbook of Fractional Calculus with Applications; Volume 2: Fractional Differential Equations; Kochubei, A., Luchko, Y., Eds.; Walter de Gruyter: Berlin, Germany; Boston, MA, USA, 2019; pp. 443–464. [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

Luchko, Y.; Yamamoto, M. The General Fractional Derivative and Related Fractional Differential Equations. Mathematics 2020, 8, 2115. https://doi.org/10.3390/math8122115

AMA Style

Luchko Y, Yamamoto M. The General Fractional Derivative and Related Fractional Differential Equations. Mathematics. 2020; 8(12):2115. https://doi.org/10.3390/math8122115

Chicago/Turabian Style

Luchko, Yuri, and Masahiro Yamamoto. 2020. "The General Fractional Derivative and Related Fractional Differential Equations" Mathematics 8, no. 12: 2115. https://doi.org/10.3390/math8122115

APA Style

Luchko, Y., & Yamamoto, M. (2020). The General Fractional Derivative and Related Fractional Differential Equations. Mathematics, 8(12), 2115. https://doi.org/10.3390/math8122115

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