Next Article in Journal / Special Issue
On Weighted (k, s)-Riemann-Liouville Fractional Operators and Solution of Fractional Kinetic Equation
Previous Article in Journal
On Discrete Delta Caputo–Fabrizio Fractional Operators and Monotonicity Analysis
Previous Article in Special Issue
Analytic Fuzzy Formulation of a Time-Fractional Fornberg–Whitham Model with Power and Mittag–Leffler Kernels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the General Solutions of Some Non-Homogeneous Div-Curl Systems with Riemann–Liouville and Caputo Fractional Derivatives

by
Briceyda B. Delgado
1 and
Jorge E. Macías-Díaz
1,2,*
1
Departamento de Matemáticas y Física, Universidad Autónoma de Aguascalientes, Aguascalientes 20100, Mexico
2
Department of Mathematics and Didactics of Mathematics, Tallinn University, 10120 Tallinn, Estonia
*
Author to whom correspondence should be addressed.
Fractal Fract. 2021, 5(3), 117; https://doi.org/10.3390/fractalfract5030117
Submission received: 9 August 2021 / Revised: 3 September 2021 / Accepted: 4 September 2021 / Published: 10 September 2021

Abstract

:
In this work, we investigate analytically the solutions of a nonlinear div-curl system with fractional derivatives of the Riemann–Liouville or Caputo types. To this end, the fractional-order vector operators of divergence, curl and gradient are identified as components of the fractional Dirac operator in quaternionic form. As one of the most important results of this manuscript, we derive general solutions of some non-homogeneous div-curl systems that consider the presence of fractional-order derivatives of the Riemann–Liouville or Caputo types. A fractional analogous to the Teodorescu transform is presented in this work, and we employ some properties of its component operators, developed in this work to establish a generalization of the Helmholtz decomposition theorem in fractional space. Additionally, right inverses of the fractional-order curl, divergence and gradient vector operators are obtained using Riemann–Liouville and Caputo fractional operators. Finally, some consequences of these results are provided as applications at the end of this work.

1. Introduction

From an analytical point of view, the literature offers a wide range of reports that focus on the extension of integer-order methods and results for the fractional case. From a more particular point of view, the fractional generalization of the classical vector calculus operators (that is, the gradient, divergence, curl and Laplacian operators) has been also an active topic of research, which has been developed from different approaches. Some of the first attempts to extend these operators to the fractional scenario are described in [1,2] using the Nishimoto fractional derivative. These operators were used later on in [3] to provide a physical interpretation for the fractional advection-dispersion equation for flow in heterogeneous porous media. In 2008, Vasily E. Tarasov described different approaches to formulate a fractional form of vector calculus with physical applications in [4] (see also references therein). More recently, a new generalization of the Helmholtz decomposition theorem for both fractional time and space was proposed in [5,6] using the discrete Grünwald–Letnikov fractional derivative. Another related work is [7], where the authors investigate the dynamic creation of fractionalized half-quantum vortices in Bose–Einstein condensates of sodium atoms.
In this work, we consider fractional derivatives of the Riemann–Liouville and the Caputo types and provide extensions of the definitions of the main differential operators from vector calculus using these fractional operators. In such a way, we present fractional forms of the divergence, the rotational and the gradient operators. Moreover, we also consider generalized forms of the Dirac and the Laplace operators using fractional derivatives. Our goal in this work is to extend quaternionic analysis to consider fractional forms of the classical differential operators. Analogues of the properties satisfied by the classical operators will be mathematically established in this work. For instance, we will show that
RL D a + α [ w ] = RL div a + α w + RL grad a + α w 0 + RL curl a + α w ,
C D a + α [ w ] = C div a + α w + C grad a + α w 0 + C curl a + α w .
In other words, when we apply the fractional Dirac operator to a quaternion valued function w = w 0 + w , this expression can be decomposed in terms of a fractional divergence, a fractional gradient and a fractional rotational operator. Based on this, we will also provide explicit expressions of general weak solutions for some fractional forms of the div-curl system, considering various analytical hypotheses. More precisely, we will prove that if we restrict ourselves to the class of functions with fractional divergence zero and whose Riemann–Liouville fractional integral has zero normal trace, then the fractional Teodorescu transform represents a solution of the above-mentioned div-curl systems (see Theorems 2 and 3).
This manuscript is organized as follows. In Section 2, we recall some important definitions from the literature, including those of the left and right Riemann–Liouville fractional derivatives, the left and right Caputo fractional derivatives and the two-parameter Mittag–Leffler function. A useful characterization of the functions with summable fractional derivatives is also recalled. In Section 2.2, the Riemann–Liouville and the Caputo Dirac and Laplace operators are introduced. Moreover, some fundamental solutions for the fractional Dirac operators are recalled in Section 2.3. In Section 3, we introduce fractional extensions of the divergence, rotational and gradient differential operators. Some properties among these operators are established, and a useful factorization theorem for the fractional Laplace operators is proven. It is worth mentioning that this factorization is not new; however, we were able to derive it only using the identities preserving the fractional gradient, divergence and rotational operators. Among the most important results, we establish that the fractional Teodorescu transform is a right inverse of the fractional Dirac operator under suitable analytical conditions, and we prove a fractional form of the Divergence Theorem.
Section 4 is devoted to establishing the existence of weak solutions for Riemann–Liouville and Caputo fractional div-curl systems. The explicit form of the operators involved in the solution, as well as some of their properties, allow the solution to be re-expressed as the sum of the fractional gradient of a scalar potential plus a fractional curl of a vector potential; we can say that our solutions preserve a Helmholtz-type decomposition (see Propositions 6 and 7). As a consequence, right inverses of the fractional rotational and divergence operators are provided in a subclass of the fractional divergence-free vector fields. In turn, Section 5 provides some consequences of the factorization results proven in Section 2 to the construction of fractional hyper-conjugate pairs. A theorem providing necessary and sufficient conditions for the existence of Caputo fractional hyper-conjugate pairs is proven in this stage, along with a result of the existence of a right inverse for the fractional gradient. Finally, this manuscript closes with a section of concluding remarks.

2. Background

2.1. Fractional Calculus

The present section is devoted to recalling some useful definitions from fractional calculus. Throughout, we assume that a , b , α R satisfy α > 0 and a < b . Meanwhile, we suppose that f : R R is a sufficiently smooth function, with the property that f is identically equal to zero outside of the interval [ a , b ] .
Definition 1.
The left and right Riemann–Liouville fractional integrals of f of order α with respect to the interval [ a , b ] (whenever they exist) are the functions I a + α [ f ] and I b α [ f ] , defined, respectively, by (see [8])
I a + α [ f ] ( x ) = 1 Γ ( α ) a x f ( t ) ( x t ) 1 α d t , x > a ,
I b α [ f ] ( x ) = 1 Γ ( α ) x b f ( t ) ( t x ) 1 α d t , x < b .
Let m = [ α ] + 1 Z . The left Riemann–Liouville and the left Caputo fractional derivatives of order α with respect to the interval [ a , b ] are, respectively, defined as follows:
R L D a + α [ f ] ( x ) = D m I a + m α [ f ] ( x ) = 1 Γ ( m α ) d m d x m a x f ( t ) ( x t ) α m + 1 d t , x > a ,
C D a + α [ f ] ( x ) = I a + m α D m [ f ] ( x ) = 1 Γ ( m α ) a x f ( m ) ( t ) ( x t ) α m + 1 d t , x > a .
Finally, we define the right Riemann–Liouville and the right Caputo fractional derivatives of order α with respect to the point a, respectively, as the functions
R L D b α [ f ] ( x ) = D m I b m α [ f ] ( x ) = ( 1 ) m Γ ( m α ) d m d x m x b f ( t ) ( t x ) α m + 1 d t , x < b ,
C D b α [ f ] ( x ) = I b m α D m [ f ] ( x ) = ( 1 ) m Γ ( m α ) x b f ( m ) ( t ) ( t x ) α m + 1 d t , x < b .
For the sake of convenience, we will employ the notation D a ± α when we present properties satisfied by both fractional derivatives RL D a ± α and C D a ± α . Throughout, we let I a + α ( L 1 ) denote the class of all functions f that are represented by the fractional integral (3) of some integrable function, i.e., f = I a + α [ g ] , for some g L 1 ( a , b ) . Using this notation, the following result provides a characterization of these functions.
Theorem 1
(Samko et al. [9]). Let α > 0 and m = [ α ] + 1 . Then, the function f belongs to I a + α ( L 1 ) if and only if I a + m α [ f ] A C m ( [ a , b ] ) , and ( I a + m α [ f ] ) k ( a ) = 0 , for each k { 0 , 1 , , m 1 } .
Definition 2.
If ( I a + m α [ f ] ) ( k ) ( a ) = 0 , for each k { 0 , 1 , m 1 } , then it follows that f ( k ) ( a ) = 0 holds, for each k { 0 , , m 1 } (see [8,9]). In light of this fact, we will say that f has a summable fractional derivative D a + α [ f ] of order α on [ a , b ] if I a + m α [ f ] A C m ( [ a , b ] ) .
In the following discussion, suppose that α > 0 , f admits a summable fractional derivative of order α > 0 on [ a , b ] , and let m = [ α ] + 1 . Then, the following composition rules are satisfied:
I a + α RL D a + α [ f ] ( x ) = f ( x ) k = 0 m 1 ( x a ) α k 1 Γ ( α k ) ( I a + m α [ f ] ) m k 1 ( a ) ,
I a + α C D a + α [ f ] ( x ) = f ( x ) k = 0 m 1 f ( k ) ( a ) k ! ( x a ) k .
On the other hand, we know that both fractional derivatives RL D a + α and C D a + α satisfy the one-sided invertibility property D a + α I a + α [ f ] = f . It is worth noting that this identity is a particular case of the property D a + α I a + β [ f ] = D a + α β [ f ] , which holds for each α , β R satisfying α β > 0 .
It is important to recall also that the following semi-group property for the composition of fractional derivatives is not generally satisfied:
C D a + α C D a + β [ f ] = C D a + β C D a + α [ f ] = C D a + α + β [ f ] .
However, if f ( k ) ( a ) = 0 for k = 0 , 1 , , max { [ α ] + 1 , [ β ] + 1 } 1 , then (11) does hold; see Section 2.2.6 [8]. An analogous condition for the semi-group property in the context of the Riemann–Liouville derivative is found in Section 2.3.6 [8]. Finally, the following relation between the Riemann–Liouville and Caputo fractional derivatives is valid:
C D a + α [ f ] ( x ) = RL D a + α [ f ] ( x ) j = 0 m 1 ( x a ) j α Γ ( j α + 1 ) D x ( j ) [ f ] ( a ) .
Definition 3
(Gorenflo et al. [10]). Let μ , ν R be such that μ > 0 . We define the two-parameter Mittag–Leffler function E μ , ν : C C with parameters μ and ν in terms of the following power series:
E μ , ν ( z ) = n = 0 z n Γ ( μ n + ν ) , z C .

2.2. Fractional Quaternionic Analysis

In this section, we will mention some recent results in fractional Clifford analysis. The Dirac operator in Clifford analysis, also known as the Moisil–Teodorescu differential operator, represents the cornerstone of the analysis in higher dimensions. A remarkable number of systems of differential equations have been analyzed using this operator or a perturbation of it, and the monographs [11,12,13] of the authors Gürlebeck and Sprößig contain many examples of the applications that have been made over the years. See also [14], where the authors introduced the fractional Dirac operator with Caputo derivatives as well as the basic tool of a fractional function theory in more dimensions.
More precisely, this section is devoted to the collection of some recent results of the authors Ferreira et al. [15,16,17], by whom fundamental solutions of the fractional Laplacian were found, where the derivatives are of Riemann–Liouville and Caputo types, as well as of the fractional Dirac operators.
For the remainder, let a i , b i R satisfy a i < b i , for each i = 1 , 2 , 3 . We will suppose that Ω = Π i = 1 3 ( a i , b i ) is a bounded open rectangular domain in R 3 , and let α = ( α 1 , α 2 , α 3 ) , with α i ( 0 , 1 ] , for all i = 1 , 2 , 3 .
Definition 4
(Ferreira et al. [15,16,17]). The fractional Riemann–Liouville and fractional Caputo Dirac operators are represented by R L D a + α and C D a + α , respectively, and they are defined as
R L D a + α = i = 1 3 e i R L x i , a i 1 + α i 2 ,
C D a + α = i = 1 3 e i C x i , a i 1 + α i 2 .
Here, R L x i , a i 1 + α i 2 and C x i , a i 1 + α i 2 are, respectively, the Riemann–Liouville and the Caputo fractional derivative operators of order ( 1 + α i ) / 2 with respect to the variable x i ( a i , b i ) , for each i = 1 , 2 , 3 . We define the fractional Laplace operators R L Δ a + α and R L Δ a + α , respectively, by
R L Δ a + α = i = 1 3 R L x i , a i 1 + α i ,
C Δ a + α = i = 1 3 C x i , a i 1 + α i ,
where R L x i , a i 1 + α i and C x i , a i 1 + α i are, respectively, the Riemann–Liouville and Caputo fractional derivatives of order 1 + α i with respect to the variable x i ( a i , b i ) , for each i = 1 , 2 , 3 .

2.3. Fundamental Solutions

The purpose of this subsection is to determine fundamental solutions for the fractional Dirac operator and use their properties in the investigation of the solutions of fractional div-curl systems. Beforehand, it is worth recalling that a family of fundamental solutions for the fractional Laplace operators RL Δ a + α and C Δ a + α , and a family of fundamental solutions for the fractional Dirac operators RL D a + α and C D a + α , were obtained in [15,16] for the Riemann–Liouville and Caputo case, respectively. In the case of Riemann–Liouville fractional operators, the authors employed some properties of the Mittag–Leffler function and the Laplace transform in two dimensions.
We will begin this section by recalling some relevant results derived in [15]. To this end, let u be an eigenfunction of the fractional Laplace operator, i.e., suppose that RL Δ a + α u = λ u for some λ C , and assume that u ( x ) admits a summable fractional derivative RL x 1 , a 1 + 1 + α 1 2 in the variable x 1 , and that it belongs to I a i + 1 + α i ( L 1 ) in the variables x 2 and x 3 . In what follows, we will consider the following integral and differential conditions of Cauchy type:
f 0 ( x 2 , x 3 ) = I a 1 + 1 α 1 [ u ] ( a 1 , x 2 , x 3 ) , f 1 ( x 2 , x 3 ) = x 1 , a 1 + α 1 [ u ] ( a 1 , x 2 , x 3 ) .
Lemma 1
(Ferreira et al. [15,17]). A family of eigenfunctions of the fractional Laplace operator R L Δ a + α is given by the function
u λ ( x ) = ( x 1 a 1 ) α 1 1 E 1 + α 1 , α 1 ( x 1 a 1 ) 1 + α 1 R L x 2 , a 2 + 1 + α 2 + R L x 3 , a 3 + 1 + α 3 λ f 0 ( x 2 , x 3 ) + ( x 1 a 1 ) α 1 E 1 + α 1 , 1 + α 1 ( x 1 a 1 ) 1 + α 1 R L x 2 , a 2 + 1 + α 2 + R L x 3 , a 3 + 1 + α 3 λ f 1 ( x 2 , x 3 ) .
Meanwhile, a family of fundamental solutions of the fractional Dirac operator R L D a + α is obtained by considering λ 0 in (19). More precisely, this family of solutions is given by
R L E a + α ( x ) = R L D a + α [ u 0 ] ( x ) ,
where u 0 is a fundamental solution of R L Δ a + α , i.e.,
u 0 ( x ) = ( x 1 a 1 ) α 1 1 E 1 + α 1 , α 1 ( x 1 a 1 ) 1 + α 1 R L x 2 , a 2 + 1 + α 2 + R L x 3 , a 3 + 1 + α 3 f 0 ( x 2 , x 3 ) + ( x 1 a 1 ) α 1 E 1 + α 1 , 1 + α 1 ( x 1 a 1 ) 1 + α 1 R L x 2 , a 2 + 1 + α 2 + R L x 3 , a 3 + 1 + α 3 f 1 ( x 2 , x 3 ) .
Here, f 0 and f 1 satisfy the conditions (18).
Let v be an eigenfunction of the fractional Laplace operator, i.e., suppose that C Δ a + α v = λ v , for some λ C . Assume that v ( x ) admits a summable fractional derivative RL x 1 , a 1 + 1 + α 1 2 in the variable x 1 , and that it belongs to I a i + 1 + α i ( L 1 ) in the variables x 2 and x 3 . By Theorem 1, v ( x 1 , a 2 , x 3 ) = v ( x 1 , x 2 , a 3 ) = 0 . In what follows, we will consider the following Cauchy conditions:
g 0 ( x 2 , x 3 ) = v ( a 1 , x 2 , x 3 ) , g 1 ( x 2 , x 3 ) = v x 1 ( a 1 , x 2 , x 3 ) .
As a consequence, g 0 ( a 2 , x 3 ) = g 0 ( x 2 , a 3 ) = g 1 ( a 2 , x 3 ) = g 1 ( x 2 , a 3 ) = 0 .
Lemma 2
(Ferreira et al. [16,17]). A family of eigenfunctions for the fractional Laplace operator C Δ a + α is given by the function
v λ ( x ) = E 1 + α 1 , 1 ( x 1 a 1 ) 1 + α 1 C x 2 , a 2 + 1 + α 2 + C x 3 , a 3 + 1 + α 3 λ g 0 ( x 2 , x 3 ) + ( x 1 a 1 ) E 1 + α 1 , 2 ( x 1 a 1 ) 1 + α 1 C x 2 , a 2 + 1 + α 2 + C x 3 , a 3 + 1 + α 3 λ g 1 ( x 2 , x 3 ) .
Meanwhile, a family of fundamental solutions of the fractional Dirac operator C D a + α is obtained by considering λ 0 in (23). More precisely, this family of solutions is given by
C E a + α ( x ) = C D a + α [ v 0 ] ( x ) ,
where v 0 is a fundamental solution of C Δ a + α , i.e.,
v 0 ( x ) = E 1 + α 1 , 1 ( x 1 a 1 ) 1 + α 1 C x 2 , a 2 + 1 + α 2 + C x 3 , a 3 + 1 + α 3 g 0 ( x 2 , x 3 ) + ( x 1 a 1 ) E 1 + α 1 , 2 ( x 1 a 1 ) 1 + α 1 C x 2 , a 2 + 1 + α 2 + C x 3 , a 3 + 1 + α 3 g 1 ( x 2 , x 3 ) ,
where g 0 ( x 2 , x 3 ) and g 1 ( x 2 , x 3 ) satisfy (22).

3. Fractional Vector Calculus

For the remainder of this section, we will study the fractional divergence, gradient and rotational operators as parts of a decomposition of the fractional Dirac operator in three dimensions. More precisely, recall that if w = w 0 + w is a quaternionic-valued function, then the following decomposition in quaternionic form is satisfied:
D w = div w + grad w 0 + curl w .
Here, D = i = 1 3 e i i is the classical Dirac operator, which is also called the Moisil–Teodorescu differential operator. For more details about quaternionic analysis, see [11,18,19]. Our goal in this section is to provide an extension of this decomposition (27) using fractional operators of the Riemann–Liouville and Caputo types.
Let w = w 0 + i = 1 3 e i w i be a quaternionic-valued function in A C ( Ω ) , whose scalar part is denoted by Sc [ w ] = w 0 and its vector part by Vec [ w ] = w = i = 1 3 e i w i . Then, the action of the operator RL D a + α on w reduces to
RL D a + α [ w ] = i = 1 3 e i RL x i , a i 1 + α i 2 [ w ] = i = 1 3 RL x i , a i 1 + α i 2 [ w i ] + i = 1 3 e i RL x i , a i 1 + α i 2 [ w 0 ] + e 1 RL x 2 , a 2 1 + α 2 2 [ w 3 ] RL x 3 , a 3 1 + α 3 2 [ w 2 ] + e 2 RL x 3 , a 3 1 + α 3 2 [ w 1 ] RL x 1 , a 1 1 + α 1 2 [ w 3 ] + e 3 RL x 1 , a 1 1 + α 2 2 [ w 2 ] RL x 2 , a 2 1 + α 2 2 [ w 1 ] .
The above decomposition of Equation (27) originates a fractional version of the classical divergence, rotational and gradient differential operators from vector calculus. These operators are, respectively, the scalar component of (27), the vector term acting over w and the vector term of the equation acting over w 0 . These facts are the motivation to analyze the following fractional differential operators.
Definition 5.
Define the fractional divergence, curl and gradient operators in the Riemann–Liouville sense by
R L div a + α w = i = 1 3 R L x i , a i 1 + α i 2 [ w i ] ,
R L curl a + α w = e 1 R L x 2 , a 2 1 + α 2 2 [ w 3 ] R L x 3 , a 3 1 + α 3 2 [ w 2 ] + e 2 R L x 3 , a 3 1 + α 3 2 [ w 1 ] R L x 1 , a 1 1 + α 1 2 [ w 3 ] + e 3 R L x 1 , a 1 1 + α 2 2 [ w 2 ] R L x 2 , a 2 1 + α 2 2 [ w 1 ] ,
R L grad a + α [ w 0 ] = i = 1 3 e i R L x i , a i 1 + α i 2 [ w 0 ] .
It is important to point out that the fractional operators (28)–(30) reduce, respectively, to the classical div, curl , and grad operators from vector calculus when α i = 1 , for each i = 1 , 2 , 3 . Moreover, if α > 0 and α i = α , for each i = 1 , 2 , 3 , then the above fractional operators coincide with the divergence, curl and gradient operators defined in [1,2] up to a constant factor. See also [4] and references therein for a historic account of the efforts to formulate a fractional form of vector calculus. Unlike the classical vector calculus operators, these fractional operators are non-local. Consequently, the fractional divergence, curl and gradient depend on the domain Ω .
Notice now that (27) can be rewritten as the following decomposition
RL D a + α [ w ] = RL div a + α w + RL grad a + α w 0 + RL curl a + α w .
Since the specific form of the Riemann–Liouville fractional derivative does not affect the above decomposition, we analogously obtain the following decomposition in terms of Caputo fractional derivatives:
C D a + α [ w ] = C div a + α w + C grad a + α w 0 + C curl a + α w .
Here, the operators C div a + α , C curl a + α and C grad a + α are defined as in (28)–(30), respectively, using Caputo fractional derivatives instead of Riemann–Liouville operators.
We define now a class of functions in AC 1 ( Ω ) where we can apply the semi-group property (11).
Definition 6.
We set Z a + α ( Ω ) = { f A C 1 ( Ω ) : f ( a 1 , x 2 , x 3 ) = f ( x 1 , a 2 , x 3 ) = f ( x 1 , x 2 , a 3 ) = 0 } .
Proposition 1.
If f = f 0 + f Z a + α ( Ω ) , then
(i)
R L div a + α R L curl a + α [ f ] = 0 ,
(ii)
R L curl a + α R L grad a + α [ f 0 ] = 0 ,
(iii)
R L div a + α R L grad a + α [ f 0 ] = R L Δ a + α [ f 0 ] ,
(iv)
R L grad a + α R L div a + α [ f ] R L curl a + α R L curl a + α [ f ] = R L Δ a + α [ f ] .
Moreover, the identities(i)–(iv)also hold for the Caputo fractional operators.
Proof. 
The results readily follow from the identities
RL x i , a i 1 + α i 2 RL x i , a i 1 + α i 2 = RL x i , a i 1 + α i , i = 1 , 2 , 3 .
RL x i , a i 1 + α i 2 RL x j , a j 1 + α j 2 = RL x j , a j 1 + α j 2 RL x i , a i 1 + α i 2 , i , j = 1 , 2 , 3 ,
which are trivially satisfied in Z a + α ( Ω ) . The identities with Caputo fractional operators are established in a similar fashion. □
Similar identities using Caputo fractional derivatives were proven in [4] when all the orders of the fractional derivatives are equal, i.e., when there exists α > 0 such that α i = α , for each i = 1 , 2 , 3 . On the other hand, a direct consequence of Proposition 1 is that the fractional Dirac operator RL D a + α factorizes the fractional Laplace operator RL Δ a + α . A more general factorization for functions taking values in C l 0 , 3 was proven in Section 4 [16]. In light of these remarks, the following result is a direct consequence of Proposition 1.
Corollary 1.
If f = f 0 + f Z a + α ( Ω ) , then the following factorizations of the fractional Laplace operators are satisfied:
R L Δ a + α [ f ] = R L D a + α R L D a + α [ f ] ,
C Δ a + α [ f ] = C D a + α C D a + α [ f ] .
We turn our attention now to the right Caputo fractional Dirac operator, which is given by the expression
C D b α = i = 1 3 e i C x i , b i 1 + α i 2 .
The following fractional Stokes formula was proven in Theorem 10 [17]:
Ω ( [ h ] C D b α ) ( y ) f ( y ) + h ( y ) ( RL D a + α [ f ] ) ( y ) d y = Ω h ( y ) η ( y ) d s y I a + α [ f ] ( y ) .
Here, we require that f , h A C 1 ( Ω ) A C ( Ω ¯ ) and I a + α [ f ] = i = 1 3 I a i + 1 α i 2 [ f ] . It is worth noting here that the operator C D b α acts on the right, while RL D a + α acts on the left. Intuitively, the last formula shows that the left Riemann–Liouville and right Caputo fractional Dirac operators act by ‘intertwining’ to obtain the fractional analogue of the Stokes formula.
The following result is a fractional form of the well-known Divergence Theorem.
Proposition 2
(Fractional Divergence Theorem). If f A C 1 ( Ω ) A C ( Ω ¯ ) , then
Ω R L div a + α [ f ] ( y ) d y = Ω η ( y ) · I a + α [ f ] ( y ) d s y ,
Ω R L curl a + α [ f ] ( y ) d y = Ω η ( y ) × I a + α [ f ] ( y ) d s y .
Proof. 
Taking h 1 in (38) and using that [ 1 ] C D b α = 0 yields that
Ω RL D a + α [ f ] ( y ) d y = Ω η ( x ) I a + α [ f ] ( y ) d s y .
Due to the decomposition (31) and because I a + α [ f ] is purely vectorial, we can readily calculate their scalar and vector parts, respectively. As a consequence, we readily achieve formulas (39) and (40), respectively. □
Proposition 3.
Let f Z a + α ( Ω ) , c 0 A C ( Ω ) and 0 < α i < 1 , for all i = 1 , 2 , 3 . Then, the following identities hold:
R L div a + α [ c 0 f ] = f · R L grad a + α [ c 0 ] ,
R L curl a + α [ c 0 f ] = f × R L grad a + α [ c 0 ] .
Proof. 
We will only calculate RL x 1 , a 1 + 1 + α 1 2 [ f i c 0 ] using integration by parts and Leibniz’ rule, the determination of RL x 2 , a 2 + 1 + α 2 2 [ f i c 0 ] and RL x 3 , a 3 + 1 + α 3 2 [ f i c 0 ] being similar. Beforehand, note that m i = 1 , for all i = 1 , 2 , 3 . Recall now that D a + α is a left inverse of I a + α ; and by hypothesis f i ( a 1 , x 2 , x 3 ) = 0 , for all i = 1 , 2 , 3 . It follows then that
RL x 1 , a 1 + 1 + α 1 2 [ f i c 0 ] ( x ) = x 1 1 Γ ( 1 α 1 2 ) a 1 x 1 f i ( t , x 2 , x 3 ) c 0 ( t , x 2 , x 3 ) ( x 1 t ) 1 + α 1 2 d t = x 1 f i ( t , x 2 , x 3 ) I a 1 + 1 α 1 2 [ c 0 ] ( t , x 2 , x 3 ) | a 1 x 1 I a 1 + 1 I a 1 + 1 α 1 2 [ c 0 ] x 1 f i = x 1 f i ( x 1 , x 2 , x 3 ) I a 1 + 1 α 1 2 [ c 0 ] ( x 1 , x 2 , x 3 ) I a 1 + 1 α 1 2 [ c 0 ] x 1 f i = f i ( x ) R L x 1 , a 1 + 1 + α 1 2 [ c 0 ] ( x ) .
Analogously, RL x j , a j + 1 + α j 2 [ f i c 0 ] = f i RL x j , a j + 1 + α j 2 [ c 0 ] , for each j = 2 , 3 and i = 1 , 2 , 3 . As a consequence,
RL div a + α [ c 0 f ] = i = 1 3 RL x i , a i + 1 + α i 2 [ c 0 f i ] = i = 1 3 f i RL x i , a i + 1 + α i 2 [ c 0 ] = f · RL grad a + α [ c 0 ] .
This establishes the first identity of the conclusion. The proof of the second equation is analogous. □
Before closing this section, it is natural to compare qualitatively the results obtained in traditional vector calculus against those in the fractional case. In classical vector calculus, the following product rules are satisfied:
div ( c 0 f ) = f · grad c 0 + c 0 div f ,
curl ( v 0 f ) = f × grad c 0 + c 0 curl f .
On the other hand, in the fully fractional case considered in Proposition 3, when we restrict f to the class of functions Z a + α ( Ω ) , the first part of these identities is also satisfied, except that the second terms on the right-hand sides of (46) and (47) are not present anymore. Notice that it is not difficult to construct a family of functions belonging to Z a + α ( Ω ) , for instance f ( x ) = ( x 1 a 1 ) γ 1 ( x 2 a 2 ) γ 2 ( x 3 a 3 ) γ 3 g ( x ) for all g AC 1 ( Ω ) and γ i 0 for all i = 1 , 2 , 3 .

4. Fractional Div-Curl Systems

4.1. Properties of the Fractional Teodorescu Transform

As a derivation of the fractional Borel–Pompeiu formula [17], the authors defined the Caputo-type Teodorescu transform in a very similar way to the following definition, the difference being that the kernel is now a fundamental solution of the fractional Dirac operator defined in (20).
Definition 7.
Let x = ( x 1 , x 2 , x 3 ) with x i > a i , for all i = 1 , 2 , 3 . Define the Riemann–Liouville and Caputo fractional Teodorescu transform by
R L T a + α [ w ] ( x ) = Ω R L E a + α ( x + a y ) w ( y ) d y ,
C T a + α [ w ] ( x ) = Ω C E a + α ( x + a y ) w ( y ) d y .
Here, we follow the nomenclature and conventions of Lemma 1. Moreover, the derivatives R L E a + α = R L D a + α [ u 0 ] and C E a + α = C D a + α [ v 0 ] that appear in the kernel of (48) and (49), respecrively, are with respect to the variable x .
In the following, it will be convenient to employ the translation operator, which is defined by T z f ( y ) = f ( y + z ) . Similarly, we will use the reflection operator given by R y f ( y ) = f ( y ) . An important relation between D a + α and T θ is
RL D ( θ + a ) + α [ T θ f ] ( x ) = RL D a + α [ f ] ( x + θ ) ,
where the derivative is taken with respect to the variable x . The proof is analogous to Theorem 11 [17], which is for the Caputo case, but we will give it here for completeness.
Proposition 4.
Let α i ( 0 , 1 ) , for each i = 1 , 2 , 3 , and let α = min 1 i 3 { α i } . The fractional Teodorescu transform R L T a + α is a right inverse of the fractional Dirac operator R L D a + α in L p ( Ω ) , for all p R satisfying 1 < p < 2 1 α .
Proof. 
Observe firstly that the fundamental solution RL E a + α ( x ) of the fractional Dirac operator RL D a + α is defined only for x i > a i , for each i = 1 , 2 , 3 . Moreover, it satisfies the identity RL D a + α RL E a + α ( x ) = δ ( x a ) or, equivalently, RL Δ a + α [ u 0 ] ( x ) = RL D a + α RL D a + α [ u 0 ] ( x ) = δ ( x a ) . In the following, the derivatives RL D y + α and RL D a + α are with respect to the variable x . Using (50) with θ = a y yields
RL D a + α RL T a + α [ w ] ( x ) = Ω RL D y + α RL E a + α ( x + a y ) w ( y ) d y = Ω RL D y + α T a y R L D a + α [ u 0 ] ( x ) w ( y ) d y = Ω T a y R L D a + α R L D a + α [ u 0 ] ( x ) w ( y ) d y = Ω T a y δ ( x a ) w ( y ) d y = Ω δ ( x y ) w ( y ) d y = w ( x ) ,
which we wished to prove. □
In the following, we will employ a key decomposition of the classical Teodorescu operator used in [20,21,22,23] for different kinds of bounded or unbounded domains in R 3 . For the fractional version, we denote the component operators of the fractional Teodorescu transform as follows:
RL T a + α [ w 0 + w ] : = RL T 0 , a + α [ w ] + RL T 1 , a + α [ w 0 ] + RL T 2 , a + α [ w ] .
The first term on the right-hand side of Equation (52) is the scalar part, while the last two summands represent the vector part, and it has been split into two components for the sake of convenience. These three terms are given by
RL T 0 , a + α [ w ] ( x ) = Ω RL E a + α ( x + a y ) · w ( y ) d y ,
RL T 1 , a + α [ w 0 ] ( x ) = Ω RL E a + α ( x + a y ) w 0 ( y ) d y ,
RL T 2 , a + α [ w ] ( x ) = Ω RL E a + α ( x + a y ) × w ( y ) d y .
We will see later in Corollary 2 how some of these component operators themselves represent right inverse operators of the fractional divergence and rotational operators, under certain conditions. Moreover, RL T a + α , as a good generalization of the classical Teodorescu operator in quaternionic analysis, preserves many of its properties. To this end, we use the above decomposition (52), RL T a + α [ g ] = RL T 0 , a + α [ g ] + RL T 2 , a + α [ g ] , to see necessary and sufficient conditions to guarantee that both its scalar part and its vector part belong to the kernel of the fractional Laplacian RL Δ a + α . In order to apply the factorization of Corollary 1 to the fractional Teodorescu transform, we will need to guarantee that the condition RL T a + α [ g ] Z a + α ( Ω ) can be satisfied.
Proposition 5.
If R L T a + α [ g ] Z a + α ( Ω ) , then the following hold for the fractional Teodorescu transform R L T a + α [ g ] :
(i)
The scalar part of R L T a + α [ g ] , R L T 0 , a + α [ g ] , belongs to the kernel of R L Δ a + α if and only if R L div a + α g = 0 .
(ii)
The vector part of R L T a + α [ g ] , R L T 2 , a + α [ g ] , belongs to the kernel of R L Δ a + α if and only if R L curl a + α g = 0 .
Moreover, the statements(i) and (ii) also hold for the Caputo fractional Teodorescu transform.
Proof. 
Using the factorization in Corollary 1, Proposition 4 and the decomposition (31), it readily follows that
RL Δ a + α RL T a + α [ g ] = RL D a + α RL D a + α RL T a + α [ g ] = RL div a + α g RL curl a + α g .
Taking its scalar part or vector part, respectively, we obtain the desired result. □

4.2. Riemann–Liouville System

In the following, we will study a fractional form of the classical div-curl system and construct its solution. More precisely, we fix the domain Ω , and consider the fractional system
RL div a + α w = g 0 , RL curl a + α w = g ,
where g 0 L p ( Ω , R ) and g L p ( Ω , R 3 ) , for some 1 < p < 2 / ( 1 α ) . Notice that if the solution is such that w Z a + α ( Ω ) , then g satisfies RL div a + α g = 0 , i.e., g is a ‘fractional divergence-free’ vector field. Let f 0 A C ( Ω ) . The following relations will be fundamental in the sequel:
RL grad a + α [ f 0 ] ( θ y ) = RL grad ( θ a ) α [ T θ R y [ f 0 ] ] ( y ) , C grad a + α [ f 0 ] ( θ y ) = C grad ( θ a ) α [ T θ R y [ f 0 ] ] ( y ) .
Here, the derivatives are taken with respect to the variable y . In the sequel and for the sake of convenience, we will employ Ker to denote the kernel of operators. As in the previous section, we will let Ω = Π i = 1 3 ( a i , b i ) be a bounded open rectangular domain in R 3 , and assume that α = ( α 1 , α 2 , α 3 ) , with α i ( 0 , 1 ) , for each i = 1 , 2 , 3 . Using this notation, we have the following result.
Theorem 2.
Let g = g 0 + g L p ( Ω ) with 1 < p < 2 / ( 1 α ) . If R L div a + α [ g ] = 0 and the normal trace of I a + α [ g ] vanishes, then a general weak solution of the fractional Riemann–Liouville div-curl system (57) is given by
w = R L T 1 , a + α [ g 0 ] + R L T 2 , a + α [ g ] + R L grad a + α [ h ] ,
where h Ker ( R L Δ a + α ) Z a + α ( Ω ) is an arbitrary scalar function.
Proof. 
By Proposition 4 and the decomposition (31), the fractional Teodorescu transform RL T a + α [ g 0 + g ] is a quaternionic solution of the system (57). To obtain a pure-vector solution, note that the decomposition (52) yields
RL T a + α [ g 0 + g ] = RL T 0 , a + α [ g ] RL T 1 , a + α [ g 0 ] + RL T 2 , a + α [ g ] .
Taking h ( y ) = T x + a R y u 0 ( y ) and f ( y ) = g ( y ) in the Stokes formula (38) (all derivatives with respect to the variable y ) and letting u 0 be a fundamental solution of RL Δ a + α given in Lemma 1, we obtain
Ω T x + a R y [ u 0 ] ( y ) η ( y ) I a + α [ g ] d s y = Ω [ T x + a R y [ u 0 ] ] C D b α ( y ) g ( y ) + T x + a R y [ u 0 ] ( y ) R L D a + α [ g ] ( y ) d y .
Calculate the scalar part of (61) and use the hypotheses Sc D a + α [ g ] = RL div a + α [ g ] = 0 in Ω and I a + α [ g ] | · η = 0 on Ω to reach
Ω C grad b α T x + a R y [ u 0 ] ( y ) · g ( y ) d y = 0 .
By differentiating now under the integral sign and using the traditional Leibniz’ rule, we readily obtain the following fundamental relation:
RL grad a + , x α [ u 0 ] ( x + a y ) = RL grad a + , y α [ u 0 ] ( x + a y ) .
Here, the second sub-index indicates whether we are taking derivatives with respect to the variable x or y . Recall that RL E a + α = RL grad a + , x α [ u 0 ] with respect to the variable x . On the other hand, due to u 0 Z a + α ( Ω ) and relations (12), (63) and (58), we obtain
RL E a + α ( x + a y ) = C grad a + , y α [ u 0 ] ( x + a y ) = C grad x , y α [ T x + a R y [ u 0 ] ] ( y ) .
However, we know that RL E a + α ( x + a y ) is only defined when x i > y i , for all i = 1 , 2 , 3 . As a consequence, we can readily replace the operator C grad b α by C grad x α in (62). Finally, by employing (64) and (62), we readily reach
RL T 0 , a + α [ g ] ( x ) = Ω RL E a + α ( x + a y ) · g ( y ) d y = 0 , in Ω .
This means that RL T a + α [ g ] = RL T 1 , a + α [ g 0 ] + RL T 2 , a + α [ g ] is purely vectorial and, moreover,
RL div a + α + RL curl a + α RL T a + α [ g 0 + g ] = RL D a + α RL T a + α [ g 0 + g ] = g 0 + g .
Setting the scalar and vector parts equal to each other, we obtain that RL T a + α [ g ] = RL T 1 , a + α [ g 0 ] + RL T 2 , a + α [ g ] is a solution of the fractional div-curl system (57). Finally, the fact that the solution is not unique is a consequence of the identities (ii) and (iii) of Proposition 1. □
In the limit α i 1 , the fractional div-curl system (57) reduces to the well-known integer-order system from vector calculus, and the hypothesis RL div a + α [ g ] = 0 reduces to the evident requirement that g be a divergence-free vector field. Moreover, I a i + 1 α i 2 [ u 0 g ] u 0 g . This means that it is sufficient to require that g has zero normal trace. Taking g 0 0 or g 0 in (57), we readily obtain
Corollary 2.
Under the same assumptions of Theorem 2, R L T 2 , a + α is a right inverse operator of R L curl a + α in the class of functions considered in Theorem 2. Meanwhile, R L T 1 , a + α is always a right inverse operator of R L div a + α in L p ( Ω ) .
Another important analogy with the classical vector calculus is that the solution of the non-linear fractional system (57) also admits a Helmholtz-type fractional decomposition as follows.
Proposition 6.
Under the same assumptions of Theorem 2, the general weak solution of the fractional Riemann–Liouville div-curl system (57) admits a fractional Helmholtz decomposition as follows
w = R L grad x α φ 0 R L curl x α φ , in Ω ,
where the scalar potential φ 0 and the vector potential φ are given by
φ 0 ( x ) = Ω T x + a R y [ u 0 ] ( y ) g 0 ( y ) d y , φ ( x ) = Ω T x + a R y [ u 0 ] ( y ) × g ( y ) d y .
Proof. 
Due to R L E a + α = R L grad a + α [ u 0 ] and by relation (58), we can write
R L T 1 , a + α [ g 0 ] ( x ) = Ω R L grad x α T x + a R y [ u 0 ] ( y ) g 0 ( y ) d y , R L T 2 , a + α [ g ] ( x ) = Ω RL grad x α T x + a R y [ u 0 ] ( y ) × g ( y ) d y .
Finally, since the fractional gradient involved in the above re-expressions of the component operators R L T 1 , a + α and R L T 2 , a + α is taken with respect to the variable x , we readily obtain (67). □

4.3. Caputo System

For the corresponding fractional div-curl system in the sense of Caputo derivatives, we will follow the same approach as that used with the Riemann–Liouville div-curl system (57). Let us consider the system
C div a + α w = g 0 , C curl a + α w = g ,
where g 0 L p ( Ω , R ) and g L p ( Ω , R 3 ) , for some 1 < p < 2 / ( 1 p ) . The cornerstone in our analysis will be again the fractional Teodorescu transform associated with Caputo derivatives, which means that its kernel is a fundamental solution of the fractional Dirac operator C D a + α . As in the case of the Riemann–Liouville fractional Teodorescu operator, we define the decomposition
C T a + α [ w 0 + w ] : = C T 0 , a + α [ w ] + C T 1 , a + α [ w 0 ] + C T 2 , a + α [ w ] ,
where C T 0 , a + α , C T 1 , a + α and C T 2 , a + α are given by (53), (54) and (55), respectively, but using now the Caputo kernel C E a + α in the integrand, instead of the Riemann–Liouville kernel R L E a + α .
Analogously to the proof of Theorem 2, it follows that C T a + α [ g 0 + g ] is a quaternionic solution of the fractional div-curl system (68). This is a direct consequence of the fact that C T a + α is a right inverse of C D a + α in L p (see [17], Theorem 11). Apply the decomposition (69) to the quaternion-valued function g = g 0 + g , where g 0 and g are the known data provided by the fractional div-curl system (68). In this way, we notice that
C T a + α [ g 0 + g ] = C T 0 , a + α [ g ] C T 1 , a + α [ g 0 ] + C T 2 , a + α [ g ] .
To obtain a purely vectorial solution, we will impose suitable conditions over g in order to guarantee that the scalar part of C T a + α [ g 0 + g ] vanishes in Ω , i.e., that C T 0 , a + α [ g ] 0 is satisfied in Ω . We will see that the fractional divergence-free functions of the Riemann–Liouville type, whose Riemann–Liouville fractional integral has zero normal trace, belong to the kernel of the operator C T 0 , a + α , in the same way as seen in Theorem 2 for the operator R L T 0 , a + α . Nevertheless, the upcoming proof is relatively more straightforward.
Theorem 3.
Let g = g 0 + g L p ( Ω ) with 1 < p < 2 / ( 1 α ) . If R L div a + α [ g ] = 0 and the normal trace of I a + α [ g ] vanishes, then a general weak solution of the fractional Caputo div-curl system (68) is given by
w = C T 1 , a + α [ g 0 ] + C T 2 , a + α [ g ] + C grad a + α [ h ] ,
where h Ker ( C Δ a + α ) Z a + α ( Ω ) is an arbitrary scalar function.
Proof. 
It suffices to prove that C T 0 , a + α [ g ] = 0 in Ω , in light of the previous discussion. Take h ( y ) = T x + a R y v 0 ( y ) and f ( y ) = g ( y ) in the Stokes formula (38), and let v 0 be the fundamental solution of C Δ a + α (25). It follows that
Ω T x + a R y [ v 0 ] ( y ) η ( y ) I a + α [ g ] d s y = Ω [ T x + a R y [ v 0 ] ] C D b α ( y ) g ( y ) + T x + a R y [ v 0 ] ( y ) R L D a + α [ g ] ( y ) d y .
Taking now the scalar part of (71) and using the hypotheses, we readily obtain
Ω C grad b α T x + a R y [ v 0 ] ( y ) · g ( y ) d y = 0 .
As a consequence of (58), we reach C E a + α ( x + a y ) = C grad x α [ T x + a R y [ v 0 ] ] ( y ) . However, we know that C E a + α ( x + a y ) is only defined for x i > y i , for all i = 1 , 2 , 3 . This implies that we can replace the operator C grad b α with C grad x α in (72). It is easy to see then that
C T 0 , a + α [ g ] ( x ) = Ω C E a + α ( x + a y ) · g ( y ) d y = 0 , in Ω ,
whence the conclusion readily follows. □
Corollary 3.
Under the same assumptions of Theorem 3, C T 2 , a + α is a right inverse operator of C curl a + α in the class of functions considered in Theorem 3. Moreover, C T 1 , a + α is a right inverse operator of C div a + α in L p ( Ω ) .
Analogously to the Riemann–Liouville case, the fractional Caputo div-curl system (68) also admits a fractional Helmholtz decomposition, but now the potential is in terms of v 0 defined in (25), which is a fundamental solution of the fractional Laplace operator of Caputo type.
Proposition 7.
Under the same assumptions of Theorem 3, the general weak solution of the fractional Caputo div-curl system (68) admits a fractional Helmholtz decomposition as follows
w = C grad x α ψ 0 C curl x α ψ , in Ω ,
where the scalar potential ψ 0 and the vector potential ψ are given by
ψ 0 ( x ) = Ω T x + a R y [ v 0 ] ( y ) g 0 ( y ) d y , ψ ( x ) = Ω T x + a R y [ v 0 ] ( y ) × g ( y ) d y .
Proof. 
The proof is analogous to that of Proposition 6 for the Riemann–Liouville case. □

5. Application

Let Ω = Π i = 1 3 ( a i , b i ) be as in the previous sections, and assume 0 < α i < 1 , for all i = 1 , 2 , 3 . The present section provides some consequences of the factorization provided by Corollary 1 to the construction of fractional hyper-conjugate pairs. In addition, we will give an explicit expression of a right inverse of the fractional gradient of Caputo type considering different derivative orders.
Definition 8.
Let w = w 0 + w A C ( Ω ) . We say that ( w 0 , w ) is a Riemann–Liouville fractional hyper-conjugate pair if w Ker ( R L D a + α ) . Analogously, ( w 0 , w ) is a Caputo fractional hyper-conjugate pair if w Ker ( C D a + α ) .
The following result is a straightforward consequence of the factorization of the fractional Laplace operator in the class Z a + α ( Ω ) provided by Corollary 1. For this reason, we omit the proof.
Corollary 4.
Let w = i = 0 3 w i , and suppose that w Ker ( R L D a + α ) Z a + α ( Ω ) (respectively, w Ker ( C D a + α ) Z a + α ( Ω ) ). Then, w i Ker ( R L Δ a + α ) (respectively, w i Ker ( C Δ a + α ) ), for all i = 0 , 1 , 2 , 3 .
By Definition 8, it is obvious that ( w 0 , w ) forms a Riemann–Liouville fractional hyper-conjugate pair if and only if the following fractional div-curl system is satisfied:
R L div a + α w = 0 , R L curl a + α w = R L grad a + α w 0 .
Similarly, ( w 0 , w ) is a Caputo fractional hyper-conjugate pair if and only if
C div a + α w = 0 , C curl a + α w = C grad a + α w 0 .
The above systems (75) and (76) can be considered fractional generalizations of the Moisil–Teodorescu system studied for the first time in [24].
Let us define the following integral operator in terms of the Riemann–Liouville fractional integrals I a i + 1 + α i 2 as
A a + α [ f ] ( x ) = I a 1 + 1 + α 1 2 [ f 1 ] ( x 1 , a 2 , a 3 ) + I a 2 + 1 + α 2 2 [ f 2 ] ( x 1 , x 2 , a 3 ) + I a 3 + 1 + α 3 2 [ f 3 ] ( x 1 , x 2 , x 3 ) .
As the following result shows, it turns out that A a + α behaves as a right-inverse operator of C grad a + α in the class of functions satisfying C curl a + α f = 0 . For this reason, A a + α is called the fractional anti-gradient operator.
Proposition 8.
If C curl a + α f = 0 , then C grad a + α A a + α [ f ] = f .
Proof. 
Using the characterization of Caputo fractional hyper-conjugate pairs given under Corollary 4 and differentiating under the integral sign, we readily obtain
C x 1 , a 1 + 1 + α 1 2 A a + α [ f ] ( x ) = C x 1 , a 1 + 1 + α 1 2 I a 1 + 1 + α 1 2 [ f 1 ] ( x 1 , a 2 , a 3 ) + I a 2 + 1 + α 2 2 [ f 2 ] ( x 1 , x 2 , a 3 ) + I a 3 + 1 + α 3 2 [ f 3 ] ( x 1 , x 2 , x 3 ) = f 1 ( x 1 , a 2 , a 3 ) + I a 2 + 1 + α 2 2 C x 1 , a 1 + 1 + α 1 2 [ f 2 ] ( x 1 , x 2 , a 3 ) + I a 3 + 1 + α 3 2 C x 1 , a 1 + 1 + α 1 2 [ f 3 ] ( x 1 , x 2 , x 3 ) .
Now, by hypothesis C curl a + α f = 0 or, equivalently, the following identities are satisfied:
C x 2 , a 2 1 + α 2 2 [ f 3 ] C x 3 , a 3 1 + α 3 2 [ f 2 ] = C x 3 , a 3 1 + α 3 2 [ f 1 ] C x 1 , a 1 1 + α 1 2 [ f 3 ] = C x 1 , a 1 1 + α 1 2 [ f 2 ] C x 2 , a 2 1 + α 2 2 [ f 1 ] = 0 .
Substituting (79) into (78) and using the composition rule (10), it follows that
C x 1 , a 1 + 1 + α 1 2 A a + α [ f ] ( x ) = f 1 ( x 1 , a 2 , a 3 ) + f 1 ( x 1 , x 2 , a 3 ) f 1 ( x 1 , a 2 , a 3 ) + f 1 ( x 1 , x 2 , x 3 ) f 1 ( x 1 , x 2 , a 3 ) = f 1 ( x ) .
Analogously, one can establish that C x i , a i + 1 + α i 2 A a + α [ f ] ( x ) = f i ( x ) , for i = 2 , 3 . The conclusion readily follows now. □
The next proposition shows that the fractional anti-gradient operator A a + α allows us to construct a Caputo fractional hyper-conjugate pair w 0 when w Ker ( C Δ a + α ) is known beforehand. The proposition is clearly a generalization of [20], Proposition 2.1.
Proposition 9.
Let w Ker ( C Δ a + α ) Z a + α ( Ω ) . A necessary and sufficient condition for the existence of a Caputo fractional hyper-conjugate pair of w is that C div a + α w = 0 . In that case, there is w 0 such that w = w 0 + w Ker ( C D a + α ) .
Proof. 
The necessity is clear due to the characterization of Caputo fractional hyper-conjugate pairs provided by (76). Suppose now that C div a + α w = 0 . By Proposition 1(iv), it follows that C curl a + α C curl a + α w = 0 . On the other hand, Proposition 8 ensures that
C grad a + α A a + α [ C curl a + α w ] = C curl a + α w .
The conclusion of this result follows from (76) if we let w 0 = A a + α [ C curl a + α w ] . □

6. Conclusions

In this work, we extend some results from vector calculus to the fractional case—for instance, the space fractional Helmholtz Decomposition Theorem provided by Propositions 6 and 7. The key tools used are the decompositions of the fractional Teodorescu transform in the Riemann–Liouville case (52) and in the Caputo case (69) as well as various properties associated with these fractional operators, which are thoroughly established in this manuscript. To this end, we consider fractional derivatives in the senses of Riemann–Liouville and Caputo, and we analyze fractional forms of various integer-order differential operators, including the divergence, the rotational, the gradient, the Dirac and the Laplace operators. As the most important result, we prove an existence theorem for the solutions of a div-curl system, considering fractional differential operators of the Riemann–Liouville and Caputo types (see Theorems 2 and 3). Other important generalizations of well-known theorems from vector calculus are proven in this way. More precisely, we present fractional versions of the classical Divergence and Stokes Theorems for vector fields (see Proposition 2). Furthermore, we focus on the construction of fractional hyper-conjugate pairs, which represent a fractional generalization of the well-known Moisil–Teodorescu system in quaternionic analysis. Finally, we note that we are also able to provide an explicit expression for an inverse of the fractional gradient operator when we restrict ourselves to vector fields whose fractional rotational is zero, when we consider fractional derivatives of the Caputo type (see Proposition 8).

Author Contributions

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

Funding

The corresponding author (J.E.M.-D.) wishes to acknowledge the financial support of the National Council for Science and Technology of Mexico (CONACYT) through grant A1-S-45928.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors wish to thank the anonymous reviewers for their comments and criticisms. All of their comments were taken into account in the revised version of the paper, resulting in a substantial improvement with respect to the original submission.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Adda, F.B. Geometric interpretation of the differentiability and gradient of real order. C. R. L’Acad. Des Sci. Ser. Math. 1998, 8, 931–934. [Google Scholar]
  2. Adda, F.B. The differentiability in the fractional calculus. Nonlinear Anal. 2001, 47, 5423–5428. [Google Scholar] [CrossRef]
  3. Meerschaert, M.M.; Mortensen, J.; Wheatcraft, S.W. Fractional vector calculus for fractional advection–dispersion. Phys. A Stat. Mech. Its Appl. 2006, 367, 181–190. [Google Scholar] [CrossRef]
  4. Tarasov, V.E. Fractional vector calculus and fractional Maxwell’s equations. Ann. Phys. 2008, 323, 2756–2778. [Google Scholar] [CrossRef] [Green Version]
  5. Ortigueira, M.D.; Rivero, M.; Trujillo, J.J. From a generalised Helmholtz decomposition theorem to fractional Maxwell equations. Commun. Nonlinear Sci. Numer. Simul. 2015, 22, 1036–1049. [Google Scholar] [CrossRef]
  6. Ortigueira, M.; Machado, J. On fractional vectorial calculus. Bull. Pol. Acad. Sci. Tech. Sci. 2018, 66, 389–402. [Google Scholar]
  7. Ji, A.C.; Liu, W.M.; Song, J.L.; Zhou, F. Dynamical creation of fractionalized vortices and vortex lattices. Phys. Rev. Lett. 2008, 101, 4. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: Amsterdam, The Netherlands, 1998; Volume 198. [Google Scholar]
  9. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives; Gordon and Breach Science Publishers: Yverdon Yverdon-les-Bains, Switzerland, 1993; Volume 1. [Google Scholar]
  10. Gorenflo, R.; Mainardi, F.; Rogosin, S. Mittag–Leffler function: Properties and applications. Handb. Fract. Calc. Appl. 2019, 1, 269–296. [Google Scholar]
  11. Gürlebeck, K.; Sprössig, W. Quaternionic Analysis and Elliptic Boundary Value Problems; Birkhäuser: Basel, Switzerland, 1990. [Google Scholar]
  12. Gürlebeck, K.; Sprößig, W. Quaternionic and Clifford Calculus for Physicists and Engineers; John Wiley & Sons: Chicheste, UK, 1997. [Google Scholar]
  13. Gürlebeck, K.; Habetha, K.; Sprößig, W. Holomorphic Functions in the Plane and n-Dimensional Space; Birkhäuser: Basel, Switzerland, 2008. [Google Scholar]
  14. Kähler, U.; Vieira, N. Fractional Clifford analysis. In Hypercomplex Analysis: New Perspectives and Applications; Bernstein, S., Kähler, U., Sabadini, I., Sommen, F., Eds.; Trends in Mathematics; Springer: Birkhäuser, Germany, 2014; pp. 191–201. [Google Scholar]
  15. Ferreira, M.; Vieira, N. Eigenfunctions and fundamental solutions of the fractional Laplace and Dirac operators: The Riemann–Liouville case. Complex Anal. Oper. Theory 2016, 10, 1081–1100. [Google Scholar] [CrossRef] [Green Version]
  16. Ferreira, M.; Vieira, N. Eigenfunctions and fundamental solutions of the fractional Laplace and Dirac operators using Caputo derivatives. Complex Var. Elliptic Equ. 2017, 62, 1237–1253. [Google Scholar] [CrossRef]
  17. Ferreira, M.; Kraußhar, R.S.; Rodrigues, M.M.; Vieira, N. A higher dimensional fractional Borel–Pompeiu formula and a related hypercomplex fractional operator calculus. Math. Methods Appl. Sci. 2019, 42, 3633–3653. [Google Scholar] [CrossRef]
  18. Sudbery, A. Quaternionic analysis. Math. Proc. Camb. Phil. Soc. 1979, 85, 99–225. [Google Scholar] [CrossRef] [Green Version]
  19. Kravchenko, V.V. Applied Quaternionic Analysis; Heldermann Verlag: Lemgo, Germany, 2003. [Google Scholar]
  20. Delgado, B.B.; Porter, R.M. General solution of the inhomogeneous div-curl system and consequences. Adv. Appl. Clifford Algebr. 2017, 27, 3015–3037. [Google Scholar] [CrossRef]
  21. Delgado, B.B.; Porter, R.M. Hilbert transform for the three-dimensional Vekua equation. Complex Var. Elliptic Equ. 2019, 64, 1797–1824. [Google Scholar] [CrossRef] [Green Version]
  22. Delgado, B.B.; Kravchneko, V.V. A Right Inverse Operator for curl+λ and Applications. Adv. Appl. Clifford Algebr. 2019, 29, 40. [Google Scholar] [CrossRef]
  23. Delgado, B.B.; Macías-Díaz, J. An Exterior Neumann Boundary-Value Problem for the Div-Curl System and Applications. Mathematics 2021, 9, 1609. [Google Scholar] [CrossRef]
  24. Moisil, G.; Theodorescu, N. Functions holomorphes dans l’espace. Mathematica 1931, 5, 142–159. [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

Delgado, B.B.; Macías-Díaz, J.E. On the General Solutions of Some Non-Homogeneous Div-Curl Systems with Riemann–Liouville and Caputo Fractional Derivatives. Fractal Fract. 2021, 5, 117. https://doi.org/10.3390/fractalfract5030117

AMA Style

Delgado BB, Macías-Díaz JE. On the General Solutions of Some Non-Homogeneous Div-Curl Systems with Riemann–Liouville and Caputo Fractional Derivatives. Fractal and Fractional. 2021; 5(3):117. https://doi.org/10.3390/fractalfract5030117

Chicago/Turabian Style

Delgado, Briceyda B., and Jorge E. Macías-Díaz. 2021. "On the General Solutions of Some Non-Homogeneous Div-Curl Systems with Riemann–Liouville and Caputo Fractional Derivatives" Fractal and Fractional 5, no. 3: 117. https://doi.org/10.3390/fractalfract5030117

Article Metrics

Back to TopTop