Next Article in Journal
Modelling Short- and Long-Term Dependencies of Clustered High-Threshold Exceedances in Significant Wave Heights
Next Article in Special Issue
Mathematical Modeling of the Electrophysical Properties of a Layered Nanocomposite Based on Silicon with an Ordered Structure
Previous Article in Journal
Real-World Data-Driven Machine-Learning-Based Optimal Sensor Selection Approach for Equipment Fault Detection in a Thermal Power Plant
Previous Article in Special Issue
Acoustics of Fractal Porous Material and Fractional Calculus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

General Fractional Vector Calculus

by
Vasily E. Tarasov
Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, 119991 Moscow, Russia
Submission received: 26 October 2021 / Revised: 3 November 2021 / Accepted: 4 November 2021 / Published: 5 November 2021
(This article belongs to the Special Issue Fractals, Fractional Calculus and Applied Statistics)

Abstract

:
A generalization of fractional vector calculus (FVC) as a self-consistent mathematical theory is proposed to take into account a general form of non-locality in kernels of fractional vector differential and integral operators. Self-consistency involves proving generalizations of all fundamental theorems of vector calculus for generalized kernels of operators. In the generalization of FVC from power-law nonlocality to the general form of nonlocality in space, we use the general fractional calculus (GFC) in the Luchko approach, which was published in 2021. This paper proposed the following: (I) Self-consistent definitions of general fractional differential vector operators: the regional and line general fractional gradients, the regional and surface general fractional curl operators, the general fractional divergence are proposed. (II) Self-consistent definitions of general fractional integral vector operators: the general fractional circulation, general fractional flux and general fractional volume integral are proposed. (III) The general fractional gradient, Green’s, Stokes’ and Gauss’s theorems as fundamental theorems of general fractional vector calculus are proved for simple and complex regions. The fundamental theorems (Gradient, Green, Stokes, Gauss theorems) of the proposed general FVC are proved for a wider class of domains, surfaces and curves. All these three parts allow us to state that we proposed a calculus, which is a general fractional vector calculus (General FVC). The difficulties and problems of defining general fractional integral and differential vector operators are discussed to the nonlocal case, caused by the violation of standard product rule (Leibniz rule), chain rule (rule of differentiation of function composition) and semigroup property. General FVC for orthogonal curvilinear coordinates, which includes general fractional vector operators for the spherical and cylindrical coordinates, is also proposed.

1. Introduction

Vector differential and integral operators are important in various fields of mechanics and physics. In standard vector calculus, integrals and derivatives of integer order are used, and therefore this mathematical tool cannot be used to describe systems, media and fields with nonlocality in space. In this regard, it is important to generalize the vector calculus for applications to the description of non-local media and systems. Nonlocality in space means a dependence of the process, states and variable at the current point of space on the changes on other points of the space.
To describe nonlocality in space, we can use derivatives and integrals of non-integer orders. There operators form a calculus if these operators of non-integer orders satisfy nonlocal analogs of fundamental theorems of standard calculus. These theorems connect the integral and differential operators of non-integer orders. Such calculus is called fractional calculus (FC), and operators are called fractional derivatives (FD) and fractional integrals (FI) [1,2,3,4,5,6,7]. The FD and FI have different nonstandard properties. For example, the standard product rule, the standard chain rule, and standard semigroup rule are violated for FD of non-integer order [8,9,10]. Fractional derivatives and integrals are actively used to describe non-standard properties of systems, media and fields with nonlocality in space and time in varios subjects of mechanics and physics [11,12,13,14,15], economics [16,17], and biology [18].
Attempts to generalize some differential operators of vector calculus have been made since the beginning of the 21st century. The history of the development of fractional vector calculus can be conditionally divided into three stages.
(1)
At the first stage, definitions of various fractional generalizations of differential vector operators (gradient, divergence, rotor, Laplace operator) were constructed. This stage began with the work of Riesz published in 1949 [19,20,21] (see also [22,23]). The definition of fractional generalizations of the gradient was proposed in the works of Adda in 1998 [24], which is actually based on the Sonin-Letnikov fractional derivative [25], and the works by Tarasov in 2005 [26,27], which are based on the Caputo fractional derivative. The definition of a fractional curl operators was proposed by Engheta in 1998 [28,29] and Meerschaert, Mortensen and Wheatcraft in 2006 [30]. The definition of a fractional divergence was proposed by Meerschaert, Mortensen and Wheatcraft in 2006 [30].
At the first stage, the proposed definitions of fractional vector differential operators were usually not consistent with each other. Fractional generalizations of integral vector operators (fractional circulation, fractional flux and fractional volume integral) have not been proposed. Fractional generalizations of fundamental theorems of vector calculus (such as the Green’s, Stokes’ and Gauss’s theorems) have not been suggested.
(2)
At the second stage, definitions of fractional generalizations of differential and integral vector operators, which are consistent with each other, were suggested. Fractional generalizations of fundamental theorems of vector calculus were proved. This phase began with the work of Tarasov published in 2008 [31] and book [13], pp. 241–264, where the power-law spatial non-locality is considered by using the Caputo fractional derivatives and Riemann-Liouville fractional integrals. The fractional generalizations of the Green’s, Stokes’ and Gauss’s theorems are formulated and proved in [13,31].
After 2008, other articles began to appear, in which special aspects of self-consistent formulations of the fractional vector calculus are discussed. Let us note these aspects: (a) the product rule for FVC is discussed by Bolster, Meerschaert and Sikorskii in 2012 [32]; (b) applications of fractional gradient to the fractional advection are considered by D’Ovidio, and Garra in 2014 [33]; (c) the discrete fractional vector calculus on lattices is proposed by Tarasov in 2014 [34]; (d) fractional generalizations of the Helmholtz decomposition are proposed by Ortigueira, Rivero and Trujillo in 2015 [35]; (e) the fractional vector operators are considered on convex domain by Agrawal and Xu in 2015 [36]; (f) the fractional vector calculus that is based on the Grunwald-Letnikov derivatives is discussed by Tarasov in 2015 [37], and then by Ortigueira and Machado in 2018 [38]; (g) the fractional Green and Gauss formulas are considered by Cheng and Dai in 2018 [39].
(3)
The third stages of the development of fractional vector calculus actually began in 2021. At this stage, the fractional vector calculus as a self-consistent mathematical theory is generalized for general form of non-locality and general form of kernels of fractional vector differential and integral operators. Self-consistent mathematical theory involves proving generalizations of all fundamental theorems of vector calculus for generalized kernels of operators. This stage began with the work of D’Elia, Gulian, Olson and Karniadakis [40] published in 2021 and based on the generalization of the Meerschaert, Mortensen and Wheatcraft approach to FVC [30]. We can also state that this stage began with our proposed work, which generalizes the approach to formulation of FVC that is proposed in basic paper [31] (see also Chapter 11 in book [13], pp. 241–264) and gave first self-consistent formulation of FVC. In the generalization of FVC from power-law nonlocality to the general form of nonlocality in space, we proposed to use the general fractional calculus (GFC) in the Luchko approach [41,42,43]. In our paper, we proposed the following:
(A)
Self-consistent definitions of general fractional differential vector operators: the regional and line general fractional gradients, the regional and surface general fractional curl operators, the general fractional divergence are proposed.
(B)
Self-consistent definitions of general fractional integral vector operators: the general fractional circulation, general fractional flux and general fractional volume integral are proposed.
(C)
The general fractional gradient, Green’s, Stokes’ and Gauss’s theorems as fundamental theorems of general fractional vector calculus are proved for simple and complex regions. The fundamental theorems (Gradient, Green, Stokes, Gauss theorems) of the proposed general FVC are proved for a wider class of domains, surfaces and curves.
All these three parts allow us to state that in this paper, for the first time, a calculus is proposed, which is a general fractional vector calculus (General FVC).
Below we describe the approach used in this paper.
In fractional calculus, nonlocality is described by the kernel of the operators, which are fractional integrals (FI) and fractional derivatives (FD) of non-integer orders. To take into account various types of nonlocality in space, we can use operators with various types of kernels. It is important to have a general fractional calculus that allows us to describe nonlocality in the most general form. We proposed [44,45] to use general fractional calculus (GFC) to describe systems, media and fields with general form of nonlocality in space.
In this paper, we proposed to use the general fractional calculus (GFC) to formulate a general fractional vector calculus (General FVC). The term “general fractional calculus” (GFC) was proposed by Kochubei in work [46] in 2011 (see also [47,48,49]). In the papers [46,47], the general fractional derivatives (GFDs) and general fractional integral (GFIs) are defined, and the fundamental theorems of the GFC are proved. The GFC is based on the concept of kernel pairs, which was proposed by N. Y. Sonin (1849–1915) in 1884 article [50] (see also [51]). (“Sonin” is more correct name of the Russian sientist [52] instead of “Sonine” that is use in French [50]). The very important form of the GFC was proposed by Luchko in 2021 [41,42,43]. In works [41,42], GFD and GFI of arbitrary order are suggested, and the general fundamental theorems for the GFI and GFDs are proved. Operational calculus for equations with general fractional derivatives is proposed in [43]. The GFC is also developed and applied in physics in works [44,45,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68].
In this paper, we use general fractional derivatives (GFDs), general fractional integrals (GFIs) and fundamental theorems of GFC as mathematical tools to formulate General FVC. The proposed General FVC is built on the basis of the results in the GFC obtained in Luchko’s works in 2021 [41,42,43].
In this paper, the proofs are detailed. This is due to the fact that many standard rules are violated for fractional integrals and derivatives, including Leibniz’s rule, the chain rule and the semigroup rule. These rules are often used in standard vector analysis to prove theorems. For example, the Stokes’ theorem is usually proved by using the chain rule, which cannot be applied to fractional derivatives and general fractional derivatives. Moreover, vector differential operators in fractional calculus become nonlocal, which creates additional difficulties for the accurate formulation and proof of theorems. Nonlocality also leads to the possibility of defining different vector differential operators, instead of one operator in the standard vector calculus. For example, we can define Regional GF Gradient, the Surface GF Gradient and Line GF Gradient. The situation is similar for GF Divergence and GF curl operators. Moreover, fractional vector analogs of fundamental theorems are not fulfilled for all these general fractional vector operations. In the general case, the general fractional (GF) gradient theorem should be considered for the line GF Gradient, the FG Stockes theorem should be considered for the surface GF Curl operator, and the GF Green theorem should be considered for the regional GF Divergence. This is due to violation of the chain rule for general fractional derivatives. In addition, violation of the chain rule leads to the fact that operators defined in different coordinate systems (Cartesian, cylindrical and spherical) are not related to each other by coordinate transformations.
In Section 2, we consider and proposed general fractional integrals and derivatives. In Section 3, we consider and proposed general fractional integral and derivative for [ a , b ] . In Section 4, we consider and proposed line general fractional integral (Line GFI). In Section 5, we give formulation and proof for general fractional gradient theorem. In Section 6, we proof general fractional Green theorem. In Section 7, we proposed double and surface general fractional integral, and flux. In Section 8, we give formulation and proof general fractional Stokes theorem. In Section 9, the general fractional Gauss theorem is formulated and proved. In Section 10, the equalities for general fractional differential vector operators are proposed. Basic interpretations of general fractional differential vector operators are described. In Section 11, we consider and proposed General FVC for orthogonal curvilinear coordinates (OCC), which includes general fractional vector operators for the spherical and cylindrical coordinates.

2. General Fractional Integrals and Derivatives

2.1. Definitions of GFI and GFD

Let us assume that the functions M ( x ) belongs to the space C 1 , 0 ( 0 , ) , and suppose that there exists a function K ( x ) C 1 , 0 ( 0 , ) , such that the Laplace convolution of these functions is equal to one for all x ( 0 , ) . The function f ( x ) belongs to the space C 1 , 0 ( 0 , ) , if this function can be represented in the form f ( x ) = x p g ( x ) , where 1 < p < 0 , and g ( x ) C [ 0 , ) .
Definition 1.
Let the functions ( M ( x ) , K ( x ) ) satisfy the following conditions.
(1) 
The Sonin condition for the kernels M ( x ) and K ( x ) requires that the equation
0 x M ( x x ) K ( x ) d x = 1
holds for all x ( 0 , ) .
(2) 
The functions M ( x ) , K ( x ) belong to the space C 1 , 0 ( 0 , ) ,
M ( x ) , K ( x ) C 1 , 0 ( 0 , ) ,
where
C 1 , 0 ( 0 , ) = { f ( x ) : f ( x ) = x p g ( x ) , x > 0 , 1 < p < 0 , g ( x ) C [ 0 , ) } .
The set of the pairs ( M ( x ) , K ( x ) ) that satisfy condition (1) is called the Sonin set and is denoted by S .
The set of the pairs ( M ( x ) , K ( x ) ) that satisfy conditions (1) and (2) is called the Luchko set and is denoted by L .
The definition of the set L was proposed in [41] to build a self-consistent formulation of a calculus of general fractional integrals and derivatives.
To define GFI and GFD, we use the Luchko’s approach to general fractional calculus, which is proposed in [41,42].
Definition 2.
Let M ( x ) L and f ( x ) C 1 ( 0 , ) = C 1 , ( 0 , ) . The general fractional integral (GFI) with the kernel M ( x ) C 1 , 0 ( 0 , ) is the operator on the space C 1 ( 0 , ) :
I ( M ) x : C 1 ( 0 , ) C 1 ( 0 , ) ,
that is defined by the equation
I ( M ) x [ x ] f ( x ) = ( M X ) ( x ) = 0 x d x M ( x x ) f ( x ) .
If the functions M ( x ) and K ( x ) belong to the Luchko set, then we can define general fractional derivatives D ( K ) x and D ( K ) x , that are associated with GFI I ( M ) x .
Definition 3.
Let M ( x ) , K ( x ) L and f ( x ) C 1 1 ( 0 , ) , i.e., f ( 1 ) ( x ) C 1 ( 0 , ) . The general fractional derivatives (GFD) with kernel K ( x ) C 1 , 0 ( 0 , ) , which is associated with GFI (5), is defined as
D ( K ) x , [ x ] f ( x ) = ( K f ( 1 ) ) ( x ) = 0 x d x K ( x x ) f ( 1 ) ( x )
for x ( 0 , ) . The GFD D ( K ) x is defined by the equation
D ( K ) x [ x ] f ( x ) = d d x ( K f ( 1 ) ) ( x ) = d d x 0 x d x K ( x x ) f ( x )
for x ( 0 , ) .
Operator (6) is called the GFD of the Caputo type, and operator (7) is called the GFD of the Riemann-Liouville type.
Remark 1.
If the kernel pair ( M ( x ) , K ( x ) ) belongs to the Luchko set L [42], the kernel K ( x ) is called associated kernel to M ( x ) . Note that if K ( x ) is associated kernel to M ( x ) , then M ( x ) is associated kernel to K ( x ) . Therefore, if ( M ( x ) , K ( x ) ) belongs to the set L , then both I ( M ) x [ x ] , D ( K ) x [ x ] and I ( K ) x [ x ] , D ( M ) x [ x ] can be used as the general fractional integrals (GFI) and general fractional derivatives (GFD).
Remark 2.
As proved in [41,42], operators (6) and (7) are connected (see Equation (47) in Definition 4 of ([41], p. 8)) by the equation
D ( K ) x , [ x ] f ( x ) = D ( K ) x [ x ] f ( x ) K ( x ) f ( 0 ) .
The GFI and GFD are connected by the fundamental theorems of general fractional calculus (FT of GFC).
Theorem 1.
(Firts fundamental theorem of GFC). If a pair of kernels ( M ( x ) , K ( x ) ) belongs to the Luchko set L , then the equality
I ( M ) x [ x ] D ( K ) x , [ x ] f ( x ) = f ( x ) f ( 0 )
holds for f ( t ) C 1 , ( K ) ( 0 , ) , where
C 1 , ( K ) ( 0 , ) : = { f : f ( x ) = I ( K ) x [ x ] g ( x ) , g ( x ) C 1 ( 0 , ) } .
Proof. 
This theorem is proved as Theorem 3 in [41], p. 9, (see also Theorem 1 in [42], p. 6). □
Theorem 2.
(Second fundamental theorem of GFC). If a pair of kernels ( M ( x ) , K ( x ) ) belongs to the Luchko set L , then the equality
I ( M ) x [ x ] D ( K ) x , [ x ] f ( x ) = f ( x ) f ( 0 )
holds for f ( x ) C 1 1 ( 0 , ) , where
C 1 1 ( 0 , ) : = { f : f ( 1 ) ( x ) C 1 ( 0 , ) } .
Proof. 
This theorem is proved as Theorem 4 in [41], p. 11, (see also Theorem 2 in [42], p. 7). □

2.2. Notations for GFI and GFD Operators

To define fractional vector operations, we introduce the operators that correspond to the fractional differentiation and general fractional integration.
The general fractional integral operator (GFI operator) is
I ( M ) x [ x ] : = 0 x d x M ( x x ) .
The operator (11) acts on functions f ( x ) C 1 ( 0 , ) by
I ( M ) x [ x ] f ( x ) = 0 x d x M ( x x ) f ( x ) .
We define the general fractional differential operator (GFD operator) in the form
D ( K ) x , [ x ] : = 0 x d x K ( x x ) d d x .
The Caputo operator (13) acts on functions on f ( x ) C 1 1 ( 0 , ) by
D ( K ) x , [ x ] f ( x ) = 0 x d x K ( x x ) f ( 1 ) ( x ) .
We note that operator (13) can be represented in the form
D ( K ) x , [ x ] : = I ( K ) x [ x ] D 1 [ x ] .
The proposed notations of GFI and GFD are more convenient, since they allow us to take into account the variables of integration and differentiations, the range of change for these variables.

2.3. Examples of Kernel Pairs from Sonin Set and Luchko Set

Let us give examples of the pair of kernels that belongs to the Sonin set S and the Luchko set S .
Note that if the kernel M ( t ) is associated to kernel K ( t ) , then the kernel K ( t ) is associated to M ( t ) . Therefore, if we have the operators
I ( M ) t [ τ ] X ( τ ) = ( M X ) ( t ) , D ( K ) t , [ τ ] X ( τ ) = ( K X ( 1 ) ) ( t ) ,
where the kernel pair ( M ( t ) , K ( t ) ) belong to the Sonin set S , then we can use the operators
I ( K ) t [ τ ] X ( τ ) = ( K X ) ( t ) , D ( M ) t , [ τ ] X ( τ ) = ( M X ( 1 ) ) ( t ) .
A similar situation with the kernel pairs ( M ( t ) , K ( t ) ) that belong to the Luchko set L .
Example 1.
The pair of the operator kernel
M ( t ) = h α ( t ) = t α 1 Γ ( α ) ,
K ( t ) = h 1 α ( t ) = t α Γ ( 1 α ) ,
where 0 < α < 1 , belongs to the Sonin set S and the Luchko set L . These kernels define the well-known Riemann-Liouville fractional integral and the Caputo fractional derivative [4].
Example 2.
The pair of the kernels [59], p. 3628:
M ( t ) = h α , λ ( t ) = t α 1 Γ ( α ) e λ t
K ( t ) = h 1 α , λ ( t ) + λ α Γ ( 1 α ) γ ( 1 α , λ t ) ,
belongs to the Sonin set S and the Luchko set L , where 0 < α < 1 , and λ 0 , t > 0 , and γ ( β , t ) is the incomplete gamma function
γ ( β , t ) = 0 t τ β 1 e τ d τ .
Example 3.
The pair of kernels
M ( t ) = ( t ) α 1 J α 1 ( 2 t ) ,
K ( t ) = ( t ) α I α ( 2 t ) ,
belongs to the Sonin set S and the Luchko set L (see [50,51], ([59], p. 3627)), if 0 < α < 1 , where
J ν ( t ) = k = 0 ( 1 ) k ( t / 2 ) 2 k + ν k ! Γ ( k + ν + 1 ) , I ν ( t ) = k = 0 ( t / 2 ) 2 k + ν k ! Γ ( k + ν + 1 )
are the Bessel and the modified Bessel functions, respectively.
For α = 0.5 , kernels (22) and (23), takes the form
M ( t ) = cos ( 2 t ) π t ,
K ( t ) = cosh ( 2 t ) π t .
Example 4.
The pair of the kernels (see Equation (7.15) in ([59], p. 3629)):
M ( t ) = t α 1 Φ ( β , α ; λ t ) ,
K ( t ) = sin ( π α ) π t α Φ ( β , 1 α ; λ t ) ,
belongs to the Sonin set S and the Luchko set L , if 0 < α < 1 , and Φ ( β , α ; z ) is the Kummer function
Φ ( β , α ; z ) = k = 0 ( β ) k ( α ) k z k k ! .
Example 5.
The pair of the kernels (see Equations (7.16) and (7.18) in ([59], p. 3629)):
M ( t ) = 1 + λ Γ ( α ) t ,
K ( t ) = 1 π t λ e λ 2 t e r f c ( λ ( t ) ) ,
belongs to the Sonin set S and the Luchko set L , if λ > 0 , and e r f c ( z ) is the complementary error function
e r f c ( z ) = 1 e r f ( z ) = 1 2 π 0 t e z 2 d z .
Example 6.
The pair of the kernels (see Equations (7.17) and (7.19) in ([59], pp. 3629–3630)):
M ( t ) = 1 λ Γ ( α ) t α 1 ,
K ( t ) = λ t α E 1 α , 1 α [ λ t 1 α ] ,
belongs to the Sonin set S and the Luchko set L , where λ > 0 .
Example 7.
An important subset H of the Sonine set was proposed by Andrzej Hanyga in work [63]. In this paper, it was proved that any singular (unbounded in a neighborhood of the point zero) locally integrable completely monotone function is a Sonine kernel. If M ( x ) belongs to the set H then its associate kernel K ( x ) also does belong to H . As an example of a kernel pair from H , we can give the kernel pair [63]
M ( t ) = h 1 β + α ( t ) + h 1 β ( t ) ,
K ( t ) = t β 1 E α , β [ t α ] ,
that belongs to the Sonin set S , if 0 < α < β < 1 , and E α , β [ z ] is the two-parameters Mittag-Leffler function
E α , β [ z ] = k = 0 z k Γ ( α k + β ) ,
where α > 0 , and β , z C . The function (35) is a singular locally integrable completely monotone function.
Example 8.
Some functions with the power-logarithmic singularities at the origin can belongs to the Sonin set S [59], pp. 3627–3630. An example of such pair of the kernels (see Equations (7.22)–(7.24) in ([59], p. 3630)) with the power-logarithmic singularities are the following
M ( t ) = A ln ( t ) Γ ( α ) t α 1 ,
K ( t ) = μ α , h ( t ) = 0 t z α e h z Γ ( z + 1 α ) d z ,
where μ α , h ( t ) is the Volterra function with
h = Γ ( α ) Γ ( α ) A .
Such kernels do not belong to the set C 1 ( 0 , ) .
Example 9.
An important subset K of the Sonine set was proposed by Anatoly N. Kochubei in works [46,47]. The kernels, which belong to the subset K , were defined in terms of their Laplace transforms and using the Sonine condition in the Laplace domain.

3. General Fractional Integral and Derivative for [ a ,   b ]

3.1. Definition of GFI and GFD for [ a ,   b ]

Let us define the general fractional integration and differentiation on [ a ,   b ] with 0 a < b .
Definition 4.
Let function F ( x ) belongs to the function space C 1 ( 0 , ) and ( M ( x ) , K ( x ) ) L .
Then, the general fractional integration on [ a , x ] with 0 a < x is defined by the equation
I [ a , x ] ( M ) [ x ] F ( x ) : = I ( M ) x [ x ] F ( x ) I ( M ) a [ x ] F ( x ) x > a > 0 , I ( M ) x [ x ] F ( x ) x > a = 0 .
For b > a > 0 and functions F ( x ) C 1 ( 0 , ) , we define the general fractional integration on [ a ,   b ] by the equation
I [ a , b ] ( M ) [ x ] F ( x ) : = I ( M ) b [ x ] F ( x ) I ( M ) a [ x ] F ( x ) =
0 b d x M ( b x ) F ( x ) 0 a d x M ( a x ) F ( x ) .
If a = 0 and b > 0 , then
I [ a , b ] ( M ) [ x ] F ( x ) : = I ( M ) b [ x ] F ( x ) .
The FGD for [ a ,   b ] are defined similarly.
Definition 5.
Let function F ( x ) belongs to the funtion space C 1 1 ( 0 , ) and ( M ( x ) , K ( x ) ) L .
Then the general fractional differentiation on [ a , x ] with 0 a < x is defined by the equation
D [ a , x ] ( K ) [ x ] F ( x ) : = D ( K ) x , [ x ] F ( x ) D ( K ) a , [ x ] F ( x ) x > a > 0 , D ( K ) x , [ x ] F ( x ) x > a = 0 .
For b > a > 0 and functions F ( x ) C 1 1 ( 0 , ) , we define the general fractional derivative on [ a ,   b ] by the equation
D [ a , b ] ( K ) [ x ] F ( x ) : = D ( K ) b , [ x ] F ( x ) D ( K ) a , [ x ] F ( x ) =
0 b d x K ( b x ) F ( 1 ) ( x ) 0 a d x M ( a x ) F ( 1 ) ( x ) .
If a = 0 and b > 0 , then
D [ a , b ] ( K ) [ x ] F ( x ) : = D ( K ) b , [ x ] F ( x ) ,
where F ( 1 ) ( x ) = d F ( x ) / d x .
The GFI operators are the map
I [ a , b ] ( M ) [ x ] : C 1 ( 0 , ) C 1 ( 0 , ) .
Remark 3.
For arbitrary a , b 0 , we can define the GFI operator as
I [ a , b ] ( M ) [ x ] F ( x ) : = sgn ( b a ) I ω [ a , b ] ( M ) [ x ] F ( x )
D [ a , b ] ( K ) [ x ] F ( x ) : = sgn ( b a ) D ω [ a , b ] ( K ) [ x ] F ( x )
with the signum function
sgn ( x ) : = 1 ( x < 0 ) , 0 ( x = 0 ) , 1 ( x > 0 ) ,
and
ω [ a , b ] : = [ a , b ] ( b > a ) , { a } ( b = a ) , [ b , a ] ( b < a ) .
We will use the notation
R + n = { ( x 1 , , x n ) : x 1 > 0 , , x n > 0 } ,
R 0 , + n = { ( x 1 , , x n ) : x 1 0 , , x n 0 } ,
where R + 1 = ( 0 , ) .

3.2. Properties of GFI and GFD on [ a   , b ]

Let us prove the additive property of the proposed GFI and GFD on [ a ,   b ] .
Theorem 3.
Let f ( x ) belongs to the space C 1 ( 0 , ) and c > b > a 0 . Then, the property is satisfied
I [ a , b ] ( M ) [ x ] f ( x ) + I [ b , c ] ( M ) [ x ] f ( x ) = I [ a , c ] ( M ) [ x ] f ( x ) .
Let f ( x ) belongs to the space C 1 1 ( 0 , ) and c > b > a 0 . Then, the equations are satisfied
D [ a , b ] ( K ) [ x ] f ( x ) + D [ b , c ] ( K ) [ x ] f ( x ) = D [ a , c ] ( K ) [ x ] f ( x ) .
Proof. 
If f ( x ) C 1 ( 0 , ) , then there are I [ a , b ] ( M ) [ x ] f ( x ) and I [ b , c ] ( M ) [ x ] f ( x ) .
For c > b > a > 0 , using the definition, we get
I [ a , b ] ( M ) [ x ] f ( x ) + I [ b , c ] ( M ) [ x ] f ( x ) =
I ( M ) b [ x ] f ( x ) I ( M ) a [ x ] f ( x ) + I ( M ) c [ x ] f ( x ) I ( M ) b [ x ] f ( x ) =
I ( M ) c [ x ] f ( x ) I ( M ) a [ x ] f ( x ) = I [ a , c ] ( M ) [ x ] f ( x ) .
For c > b > a = 0 , using the definition, we have
I [ 0 , b ] ( M ) [ x ] f ( x ) + I [ b , c ] ( M ) [ x ] f ( x ) =
I ( M ) b [ x ] f ( x ) + I ( M ) c [ x ] f ( x ) I ( M ) b [ x ] f ( x ) =
I ( M ) c [ x ] f ( x ) = I [ 0 , c ] ( M ) [ x ] f ( x ) .
If f ( x ) C 1 1 ( 0 , ) , then there are D [ a , b ] ( K ) [ x ] f ( x ) and D [ b , c ] ( K ) [ x ] f ( x ) . The proof of the additivity of the operators is proved in a similar way. □

3.3. Fundamental Theorems for GFI and GFD on [a, b]

Let us formulate the second FT of GFC for the operator I [ a , x ] ( M ) [ s ] .
Theorem 4.
Let f ( x ) belongs to the space C 1 1 ( 0 , ) , and the pair of kernels ( M ( t ) , K ( t ) ) belongs to the Luchko set L . Then
I [ a , x ] ( M ) [ s ] D ( K ) s , [ x ] f ( x ) = f ( x ) f ( a ) ,
where x > a 0 .
Proof. 
Let us assume that f ( s ) C 1 1 ( 0 , ) and
F ( s ) = D ( K ) s , [ x ] f ( x ) C 1 ( 0 , ) .
Then using the second fundamental theorem GFC in the form
I ( M ) x [ s ] D ( K ) s , [ x ] f ( x ) = f ( x ) f ( 0 ) ,
which holds for all x > 0 , we obtain
I [ a , x ] ( M ) [ s ] D ( K ) s , [ s ] f ( s ) : = I ( M ) x [ s ] D ( K ) s , [ x ] f ( x ) I ( M ) a [ s ] D ( K ) s , [ x ] f ( x ) =
f ( x ) f ( 0 ) f ( a ) f ( 0 ) = f ( x ) f ( a ) .
 □
Corollary 1.
For x = b , Equation (42) takes the form
I [ a , b ] ( M ) [ x ] D ( K ) x , [ x ] f ( x ) = f ( b ) f ( a ) .
Let us formulate the first FT of GFC for the operator I [ a , x ] ( M ) [ s ] .
Theorem 5.
Let us assume that the conditions of the first FT of GFC for the operators I ( M ) x and GFD D ( K ) x , are satisfied. Then
D ( K ) x , [ s ] I [ a , s ] ( M ) [ x ] F ( x ) = F ( x ) ,
where we use
D ( K ) x , [ s ] 1 = 0 .
Proof. 
Then
D ( K ) x , [ s ] I [ a , s ] ( M ) [ x ] F ( x ) : = D ( K ) x , [ s ] I ( M ) s [ x ] F ( x ) D ( K ) x , [ s ] I ( M ) a [ x ] F ( x ) =
F ( x ) I ( M ) a [ x ] F ( x ) D ( K ) x , [ s ] 1 = F ( x ) ,
where we use
D ( K ) x , [ s ] 1 = 0 .
 □
For functions F ( x ) C 1 1 ( 0 , ) , we define the general fractional derivative on [ a ,   b ] by the equation
D [ a , x ] ( K ) [ x ] F ( x ) : = D ( K ) x , [ x ] F ( x ) D ( K ) a , [ x ] F ( x ) .
For the operator D [ a , x ] ( K ) and GFI I ( M ) x , the first fundamental theorem GFC is satisfied in the following form.
Theorem 6.
Let us assume that the conditions of the second FT of GFC for the operators I ( M ) x and GFD D ( K ) x , are satisfied. Let us assume that f ( s ) C 1 ( 0 , ) and
F ( s ) = I ( M ) s [ x ] f ( x ) C 1 1 ( 0 , ) .
Then
D [ a , x ] ( K ) [ s ] I ( M ) s [ x ] f ( x ) = f ( x ) f ( a ) .
Proof. 
Using the first fundamental theorem GFC in the form
D ( K ) x , [ x ] I ( M ) x , [ s ] f ( s ) = f ( x ) ,
which holds for all x > 0 , we obtain
D [ a , x ] ( K ) [ s ] I ( M ) s [ x ] f ( x ) : = D ( K ) x , [ s ] I ( M ) s [ x ] f ( x ) D ( K ) a , [ s ] I ( M ) s [ x ] f ( x ) =
f ( x ) f ( a ) ,
where x > a 0 . □
Remark 4.
The second fundamental theorem GFC does not hold for the operator D [ a , x ] ( K ) with a > 0 , since I ( M ) x [ x ] 1 0 .

4. Line General Fractional Integral (Line GFI)

4.1. Simple Line in R 0 , + 2

Let us define concept of simple line in R 0 , + 2 of the X Y -plane.
Definition 6.
Let a line L R 0 , + 2 be described by the function
y = y ( x ) 0 , x [ a , b ] R 0 , + ,
which is a continuously differentiable function for all x ( a , b ) R 0 , + , i.e., y ( x ) C 1 ( a , b ) . Then the line L will be called Y-simple line on the X Y plane.
Remark 5.
It is possible to weaken the condition of continuous differentiability of the function y ( x ) at all points of the interval ( a , b ) , i.e., y ( 1 ) ( x ) C ( a , b ) . We can consider the Y-simple lines, for which
y ( 1 ) ( x k 0 ) y ( 1 ) ( x k + 0 )
for finite number of points x k [ a , b ] with y ( x k 0 ) = y ( x k + 0 ) . For example, these conditions can be used for broken lines (polygonal chains).
If L is Y-simple line, then every line, which is parallel to the y-axis, intersects the line L at most once for x [ a , b ] .
Similarly, line L R 0 , + 2 is called X-simple line on the X Y plane, if L can be described by the function
x = x ( y ) 0 ,
which is a continuously differentiable function for y ( c , d ) R 0 , + , x ( y ) C 1 ( c , d ) .
Definition 7.
The line L in R 0 , + 2 of the X Y plane is called simple line in R 0 , + 2 of the X Y plane, if L is X-simple and Y-simple line.
Let L R 0 , + 2 is X-simple line on the X Y plane that is described by the single valued function y = y ( x ) for x [ a , b ] . If the functions has inverse function, then L is simple line on the X Y plane. It is known that if the function y = f ( x ) has the derivative f ( 1 ) ( x ) = d f ( x ) / d x for x [ a , b ] such that f ( 1 ) ( x ) > 0 (or f ( 1 ) ( x ) < 0 ) for all x [ a , b ] , then there exists an inverse function x = f 1 ( y ) for y [ c , d ] , where c = y ( a ) and d = y ( b ) .
Therefore we can formulate the following theorem.
Theorem 7.
Let L R 0 , + 2 is Y-simple line on the X Y plane, where y = y ( x ) C 1 ( a , b ) and y ( 1 ) ( x ) > 0 (or y ( 1 ) ( x ) < 0 ) for all x [ a , b ] .
Then there is an inverse function x = x ( y ) C 1 ( c , d ) , the line L is X-simple line with y [ c , d ] , and L is simple line.
Remark 6.
In general FVC, we can consider smooth lines on R 0 , + that consist of lines, which are simple lines with respect to one of the axes and lines parallel to one of the axes (X, Y, Z).

4.2. Simple Line in R 0 , + 3

Let us define concepts of simple line in R 0 , + 3 .
Definition 8.
Let L R 0 , + 3 be a line that is described by the functions
y = y ( x ) 0 , z = z ( x ) 0 , x [ a , b ] R 0 , + ,
which are continuously differentiable functions for x ( a , b ) R 0 , + , i.e., y ( x ) , z ( x ) C 1 ( a , b ) . Then this line will be called Y Z -simple line.
If L is Y Z -simple line, then every plane, which is parallel to the Y Z -plane, intersects the line L at most once for [ a ,   b ] .
Definition 9.
The line L is called simple line in R 0 , + 3 , if L is X Y -, X Z - and Y Z -simple line.
If the derivatives y x ( 1 ) ( x ) and z x ( 1 ) ( x ) are nonzero and do not change sign on the interval x [ a , b ] , then there exist the inverse functions x = x ( y ) C 1 ( c , d ) and x = x ( z ) C 1 ( e , f ) . The derivatives of the functions x ( y ) and x ( z ) are defined by the equations
x y ( 1 ) = 1 y x ( 1 ) : x y ( 1 ) ( y ) = 1 y x ( 1 ) ( x ) x = x ( y ) .
It is obvious that these derivatives of the functions are nonzero and do not change sign also.
We can also state that there exist the functions
y = y ( z ) : = y ( x ( z ) ) C 1 ( e , f ) , z = z ( y ) : = z ( x ( y ) ) C 1 ( c , d ) ,
and the derivatives of these functions that are are nonzero and do not change sign.
As a result, we proved the following theorem.
Theorem 8.
The Y Z -simple line L R 0 , + 3 , for which the derivatives of the functions y = y ( x ) and z = z ( x ) are nonzero and do not change sign on the interval x [ a , b ] , is X Y -simple, X Z -simple line, and therefore simple line in R 0 , + 3 .
Let L be a simple line in R 0 , + 3 and the lines L x y , L x z and L y z are projections of L on the X Y -, X Z -, Y Z planes. Then L x y , L x z and L y z are simple lines in R 0 , + 2 of the X Y -, X Z -, Y Z planes.
Remark 7.
The simple line L R 0 , + 3 , which connects the points A ( a , c , e ) and B ( b , d , f ) , can be defined by three equivalent forms
L = { ( x , y , z ) : x [ a , b ] , y = y ( x ) C 1 [ a , b ] , z = z ( x ) C 1 [ a , b ] } ,
L = { ( x , y , z ) : y [ c , d ] , x = x ( y ) C 1 [ c , d ] , z = z ( y ) C 1 [ c , d ] } ,
L = { ( x , y , z ) : z [ e , f ] , x = x ( z ) C 1 [ e , f ] , y = y ( z ) C 1 [ e , f ] } ,
where y ( x = a ) = c , y ( x = b ) = d , z ( x = a ) = e , z ( x = b ) = f .
Remark 8.
We can also consider the simple lines, for which the function that are not differentiable at a finite number of points. These conditions allows us to use for broken lines (polygonal chains).

4.3. Problems with Definition of Line GFI of Vector Field

Let us consider the Y-simple line L 1 R 0 , + 2 in the X Y -plane that is described by the equation
y = y ( x ) C 1 [ a , b ] ,
and the vector field
F ( x , y ) = e x F x ( x , y ) + e y F y ( x , y ) ,
where ( x , y ) R 0 , + 2 . Then the standard line (curve) integral of second kind can be defined by the equation
L 1 F x ( x , y ) d x + F y ( x , y ) d y : = a b F x ( x , y ( x ) ) d x + F y ( x , y ( x ) ) y ( 1 ) ( x ) d x .
If the X-simple line L 2 R 0 , + 2 in the X Y -plane is described by the equation
x = x ( y ) C 1 [ c , d ] ,
then the line integral can be defined by the equation
L 2 F x ( x , y ) d x + F y ( x , y ) d y : = c d F x ( x ( y ) , y ) x ( 1 ) ( y ) d y + F y ( x ( y ) , y ) d y .
If y = y ( x ) C 1 ( a , b ) and y ( 1 ) ( x ) > 0 (or y ( 1 ) ( x ) < 0 ) for all x [ a , b ] . Then there is an inverse function x = x ( y ) C 1 ( c , d ) , and the Y-simple line L 1 is X-simple line. Therefore the line integral (49) can be represented as (50), and we have the equality
a b F x ( x , y ( x ) ) d x + F y ( x , y ( x ) ) y ( 1 ) ( x ) d x = c d F x ( x ( y ) , y ) x ( 1 ) ( y ) d y + F y ( x ( y ) , y ) d y .
This equality is based on the property
c d F y ( x ( y ) , y ) d y = a b F y ( x , y ( x ) ) y ( 1 ) ( x ) d x
a b F x ( x , y ( x ) ) d x = c d F x ( x ( y ) , y ) x ( 1 ) ( y ) d y .
Equations (49) and (50) cannot be used to define the line general fractional integral (line GFI) since property (52) and (53) does not hold for general fractional integrals. In general fractional calculus, we have the inequality
I [ c , d ] ( M 2 ) [ y ] F y ( x ( y ) , y ) I [ a , b ] ( M 1 ) [ x ] F y ( x , y ( x ) ) y ( 1 ) ( x )
that has the form
0 d M 2 ( d y ) F y ( x ( y ) , y ) d y 0 c M 2 ( c y ) F y ( x ( y ) , y ) d y
0 b M 1 ( b x ) F y ( x , y ( x ) ) y ( 1 ) ( x ) d x 0 a M 1 ( a x ) F y ( x , y ( x ) ) y ( 1 ) ( x ) d x .
To solve this problem of definition of line GFI, we can use the fact that the line integral of second kind over simple line can be defined by the equation
L 1 F x ( x , y ) d x + F y ( x , y ) d y : = a b F x ( x , y ( x ) ) d x + c d F y ( x ( y ) , y ) d y .
Therefore the line GFI in R 0 , + 2 can be defined by the following definition.
Definition 10.
Let L be a simple line in R 0 , + 2 of the X Y -plane. Let the functions f x ( x ) : = F x ( x , y ( x ) ) and g y ( y ) : = F y ( x ( y ) , y ) belong to the function space C 1 ( 0 , ) .
Then line GFI for the line L is defined by the equation
I L ( M ) , F = I [ a , b ] ( M 1 ) [ x ] F x ( x , y ( x ) ) + I [ c , d ] ( M 2 ) [ y ] F y ( x ( y ) , y ) =
0 b M 1 ( b x ) F x ( x , y ( x ) ) d x 0 a M 1 ( a x ) F x ( x , y ( x ) ) d x +
0 d M 2 ( d y ) F y ( x ( y ) , y ) d y 0 c M 2 ( c y ) F y ( x ( y ) , y ) d y .
This line GFI exists, if the kernels ( M 1 ( x ) , K 1 ( x ) ) and ( M 2 ( y ) , K 2 ( y ) ) belong to the Luchko set L .
The proposed approach to define line GFI for lines in R 0 , + 2 can be used to define line GFI for lines in R 0 , + 3 . Let L be a simple line in R 0 , + 3 , which is defined in form (46)–(48), and the vector field
F ( x , y , z ) = e x F x ( x , y , z ) + e y F y ( x , y , z ) + e z F z ( x , y , z ) ,
where ( x , y , z ) R 0 , + 3 . The standard line integral of second kind over a simple line L in R 0 , + 3 can be defined by the equation
L F x ( x , y , z ) d x + F y ( x , y , z ) d y + F z ( x , y , z ) d z : = a b F x ( x , y ( x ) , z ( x ) ) d x + c d F y ( x ( y ) , y , z ( y ) ) d y + e f F z ( x ( z ) , y ( z ) , z ) d z .
The line GFI over a simple line L in R 0 , + 3 can be defined by the equation
I L ( M ) , F = I [ a , b ] ( M 1 ) [ x ] F x ( x , y ( x ) , z ( x ) ) + I [ c , d ] ( M 2 ) [ y ] F y ( x ( y ) , y , z ( y ) ) + I [ e , f ] ( M 3 ) [ x ] F z ( x ( z ) , y ( z ) , z ) .
This line GFI exists, if the kernels ( M 1 ( x ) , K 1 ( x ) ) , ( M 2 ( y ) , K 2 ( y ) ) and ( M 3 ( z ) , K 3 ( z ) ) belong to the Luchko set L , and the functions F 1 ( x ) : = F x ( x , y ( x ) , z ( x ) )   F 2 ( y ) : = F y ( x ( y ) , y , z ( y ) ) , and F 3 ( z ) : = F z ( x ( z ) , y ( z ) , z ) belong to the function space C 1 ( 0 , ) .
Remark 9.
The line integral can also be defined for wide class of lines in R 0 , + 3 that are not simple lines, if the lines L can be split into finite number of simple lines in R 0 , + 3 . These lines will be called piecewise simple lines.

4.4. Definition of Line GFI for Vector Field in R 0 , + 3

Let us define some conditions on the vector field. We will assume that the vector field
F = F ( x , y , z ) = e x F x ( x , y , z ) + e y F y ( x , y , z ) + e z F z ( x , y , z )
on the simple line L R 0 , + 3 is described by the following functions that belong to the space C 1 ( 0 , ) :
F 1 ( x ) : = F x ( x , y ( x ) , z ( x ) ) C 1 ( 0 , ) ,
F 2 ( y ) : = F y ( x ( y ) , y , z ( y ) ) C 1 ( 0 , ) ,
F 3 ( z ) : = F z ( x ( z ) , y ( z ) , z ) C 1 ( 0 , ) .
If these conditions are satisfied, then we will write F ( x , y , z ) F 1 ( L ) .
In the case F F 1 ( R + 3 ) , the line general fractional integral for the vector field F and the line L = A B with endpoints A ( a 1 , a 2 , a 3 ) and B ( b 1 , b 1 , b 3 ) with all b k > a k 0 is defined by the equation
I L ( M ) , F = k = 1 3 I [ a k , b k ] ( M k ) [ x k ] F k ( x k ) =
I [ a 1 , b 1 ] ( M 1 ) [ x ] F 1 ( x ) + I [ a 2 , b 2 ] ( M 2 ) [ y ] F 2 ( y ) + I [ a 3 , b 3 ] ( M 3 ) [ z ] F 3 ( z ) =
I ( M 1 ) b 1 [ x ] I ( M 1 ) a 1 [ x ] F x ( x , y ( x ) , z ( x ) ) +
I ( M 2 ) b 2 [ x ] I ( M 2 ) a 2 [ y ] F y ( x ( y ) , y , z ( y ) ) +
I ( M 3 ) b 3 [ z ] I ( M 3 ) a 3 [ z ] F z ( x ( z ) , y ( z ) , z ) =
0 b 1 d x M 1 ( b 1 x ) F x ( x , y ( x ) , z ( x ) ) 0 a 1 d x M 1 ( a 1 x ) F x ( x , y ( x ) , z ( x ) ) +
0 b 2 d y M 2 ( b 2 y ) F y ( x ( y ) , y , z ( y ) ) 0 a 2 d y M 2 ( a 2 y ) F y ( x ( y ) , y , z ( y ) ) +
0 b 3 d z M 3 ( b 3 z ) F z ( x ( z ) , y ( z ) , z ) 0 a 3 d z M 3 ( a 3 z ) F z ( x ( z ) , y ( z ) , z ) .
For the proposed definition of a line general fractional integral, the kernels of the integral operator remain dependent on the difference of variables
M k = M k ( b k x k ) , M k = M k ( a k x k ) ,
and the line GFI itself is expressed through the Laplace convolution as a product in the ring R 1 . Using the Laplace convolution, the line GFI can be written as
I L ( M ) , F = k = 1 3 ( M k F k ) ( b k ) ( M k F k ) ( a k ) ,
where b k > a k 0 , and
( M k F k ) ( c k ) = 0 c k d x k M ( c k x k ) F k ( x k ) ,
where c k = a k or c k = b k with M k ( x k ) C 1 , 0 ( 0 , ) and F k ( x k ) C 1 ( 0 , ) . Using the property of GFI, we have ( M k F k ) ( c k ) C 1 ( 0 , ) .
We can consider the variables x k > a k instead of the numbers b k . Then
I L ( M ) , F = k = 1 3 ( M k F k ) ( x k ) ( M k F k ) ( a k ) ,
where x k > a k 0 .
Remark 10.
In general FVC, we can consider lines that consist of lines, which are simple lines and lines parallel to one of the axes (X, Y, Z). As an example of this type of lines, we can consider polygonal chains (broken lines).

4.5. Line General Fractional Integral for Polygonal Chains

Let us define a line GFI for the polygonal chains. A polygonal chain (broken line) is a geometric figure consisting of line segments connected in series by their endpoints. A polygonal chain L is a sequence of points A 1 , A 2 , …, A n that forms the successively connected line segments [ A 1 , A 2 ] , [ A 2 , A 3 ] , …, [ A n 1 , A n ] . The points A 1 , A 2 , …, A n are called the vertices of polygonal chain.
Let us consider the polygonal chain
L : = k = 1 n 1 [ A k , A k + 1 ] ,
where the vertices A k have coordinates ( x k , y k , z k ) , such that x k , y k , z k 0 for all k = 1 , 2 , , n 1 .
For polygonal chain (57), the line GFI is defined by the equations
I L ( M ) [ x , y , z ] = k = 1 n 1 I [ A k , A k + 1 ] ( M ) [ x , y , z ] ,
where
I [ A k , A k + 1 ] ( M ) [ x , y , z ] : = sgn ( x k + 1 x k ) e x I ω [ x k , x k + 1 ] ( M 1 ) [ x ] +
sgn ( y k + 1 y k ) e y I ω [ y k , y k + 1 ] ( M 2 ) [ y ] + sgn ( z k + 1 z k ) e z I ω [ z k , z k + 1 ] ( M 3 ) [ z ] ,
with the signum function
sgn ( x ) : = 1 ( x < 0 ) , 0 ( x = 0 ) , 1 ( x > 0 ) ,
and
ω [ x k , x k + 1 ] : = [ x k , x k + 1 ] ( x k + 1 > x k ) , { x k } ( x k + 1 = x k ) , [ x k + 1 , x k ] ( x k + 1 < x k ) .
If x k < x k + 1 , y k < y k + 1 , z k < z k + 1 , then
I [ A k , A k + 1 ] ( M ) [ x , y , z ] : = e x I [ x k , x k + 1 ] ( M 1 ) [ x ] + e y I [ y k , y k + 1 ] ( M 2 ) [ y ] + e z I [ z k , z k + 1 ] ( M 3 ) [ z ] ,
If the opposite inequality holds, then a minus sign is put in front of the integral, and in the interval the numbers are set in ascending order. If equality ( x k = x k + 1 or y k = y k + 1 or z k = z k + 1 ) holds, then the integral is considered equal to zero.
Let vector field F be defined by the equation
F = F ( x , y , z ) = e x F x ( x , y , z ) + e y F y ( x , y , z ) + e z F z ( x , y , z ) ,
where
F x ( x , y k ( x ) , z k ( x ) ) , F y ( x k ( y ) , y , z k ( y ) ) , F y ( x k ( y ) , y , z k ( y ) ) C 1 ( 0 , )
for all k = 1 , , n .
Then the line GFI for the polygonal chain (broken line) is defined by Equation (58), where
I [ A k , A k + 1 ] ( M ) [ x , y , z ] F ( x , y , z ) : = sgn ( x k + 1 x k ) I ω [ x k , x k + 1 ] ( M 1 ) [ x ] F x ( x , y k ( x ) , z k ( x ) ) +
sgn ( y k + 1 y k ) I ω [ y k , y k + 1 ] ( M 2 ) [ y ] F y ( x k ( y ) , y , z k ( y ) ) +
sgn ( z k + 1 z k ) I ω [ z k , z k + 1 ] ( M 3 ) [ z ] F z ( x k ( z ) , y k ( z ) , z ) ,
where the functions y = y k ( x ) , x = x k ( y ) , z = z k ( x ) , x = x k ( z ) , z = z ( y ) , y = y k ( x ) are defined by the equations
x x k x k + 1 x k = y y k y k + 1 y k = z z k z k + 1 z k ,
if the denominators are not zero. For example,
y = y k + 1 y k x k + 1 x k ( x x k ) + y k ,
z = z k + 1 z k x k + 1 x k ( x x k ) + z k .
Here
I [ a , b ] ( M 1 ) [ x ] F x ( x , y ( x ) , z ( x ) ) =
0 b d x M 1 ( b x ) F x ( x , y ( x ) , z ( x ) ) 0 a d x M 1 ( b x ) F x ( x , y ( x ) , z ( x ) ) ,
where b > a 0 ,
I [ c , d ] ( M 2 ) [ x ] F y ( x ( y ) , y , z ( y ) ) =
0 d d y M 2 ( d x ) F y ( x ( y ) , y , z ( y ) ) 0 c d y M 2 ( c x ) F y ( x ( y ) , y , z ( y ) ) ,
where d > c 0 ,
I [ e , f ] ( M 3 ) [ x ] F y ( x ( z ) , y ( z ) , z ) =
0 f d z M 3 ( f x ) F z ( x ( z ) , y ( z ) , z ) 0 e d z M 3 ( e x ) F z ( x ( z ) , y ( z ) , z ) ,
where f > e 0 .

4.6. Line GFI for Piecewise Simple Lines

Similarly to the case of a broken line (the polygonal chains), we can define line GFI for line, which consists of simple lines and lines parallel to the axes.
Let us consider a line L R 0 , + 3 , which can be divided into several lines L k = L k [ A k , A k + 1 ] , k = 1 , , n that are simple lines or lines parallel to one of the axes:
L : = k = 1 n 1 L k [ A k , A k + 1 ] ,
where the line L k connects the points A k ( x k , y k , z k ) , and A k + 1 ( x k + 1 , y k + 1 , z k + 1 ) with x k , y k , z k 0 for all k = 1 , 2 , , n 1 . Lines of this kind will be called the piecewise simple lines.
For piecewise simple line (61) in R 0 , + 3 , and the vector field F F 1 ( L ) , the line GFI is defined by the equation
I L ( M ) [ x , y , z ] F ( x , y , z ) : = k = 1 n 1 I [ A k , A k + 1 ] ( M ) [ x , y , z ] F ( x , y , z ) ,
where
I [ A k , A k + 1 ] ( M ) [ x , y , z ] F ( x , y , z ) : = sgn ( x k + 1 x k ) I ω [ x k , x k + 1 ] ( M 1 ) [ x ] F x ( x , y k ( x ) , z k ( x ) ) +
sgn ( y k + 1 y k ) I ω [ y k , y k + 1 ] ( M 2 ) [ y ] F y ( x k ( y ) , y , z k ( y ) ) +
sgn ( z k + 1 z k ) I ω [ z k , z k + 1 ] ( M 3 ) [ z ] F z ( x k ( z ) , y k ( z ) , z ) ,
where the functions y = y k ( x ) , x = x k ( y ) , z = z k ( x ) , x = x k ( z ) , z = z k ( y ) , y = y k ( x ) define the lines L k that are simple or parallel to one of the axes.

4.7. General Fractional Circulation for Rectangle

Let us define concepts of general fractional circulation by using the GFC. Note that these concepts for the kernels
M n ( x n ) = h α n ( x n ) = x n α n 1 Γ ( α n ) , ( for all n = 1 , 2 , 3 )
were proposed in [31] (see also [13], pp. 241–264).
We will consider the piecewise simple line (61), where A n ( a n , c n , e n ) = A 1 ( a 1 , c 1 , e 1 ) , i.e., L is closed line. Let L be the piecewise simple line
L : = k = 1 n 1 L k ,
where the lines L k = ( A k , A k + 1 ) with points A k ( a k , c k , d k ) ( a k + 1 = b k , c k + 1 = d k , e k + 1 = f k ) are simple that are described as
L k : = { ( x , y , z ) : x [ a k , b k ] , y = y k ( x ) C 1 ( 0 , ) , z = z k ( x ) C 1 ( 0 , ) } ,
where the derivatives of the functions y k ( x ) , z k ( x ) are nonzero and do not change sign on the interval [ a k , b k ] foe all k = 1 , , n .
The last requirement for derivatives of the functions y k ( x ) , z k ( x ) allows us to represent the line in the form
L k : = { ( x , y , z ) : y [ c k , d k ] , x = x k ( y ) C 1 ( 0 , ) , z = z k ( y ) C 1 ( 0 , ) } ,
L k : = { ( x , y , z ) : z [ e k , f k ] , x = x k ( z ) C 1 ( 0 , ) , y = y k ( z ) C 1 ( 0 , ) } ,
where y k ( a k ) = c k , y k ( b k ) = d k , z k ( a k ) = e k , z k ( b k ) = f k .
Then, we can define the general fractional circulation in the following form
E L ( M ) ( F ) = I L ( M ) [ x , y , z ] F ( x , y , z ) = k = 1 n 1 I L k ( M ) [ x , y , z ] F ( x , y , z ) = k = 1 n 1 I L k ( M 1 ) [ x ] F x ( x , y k ( x ) , z k ( x ) ) + k = 1 n 1 I L k ( M 2 ) [ y ] F y ( x k ( y ) , y , z k ( y ) ) + k = 1 n 1 I L k ( M 3 ) [ z ] F z ( x k ( z ) , y k ( z ) , z ) = k = 1 n 1 I [ a k , b k ] ( M 1 ) [ x ] F x ( x , y k ( x ) , z k ( x ) ) + k = 1 n 1 I [ c k , d k ] ( M 2 ) [ y ] F y ( x k ( y ) , y , z k ( y ) ) + k = 1 n 1 I [ e k , f k ] ( M 3 ) [ z ] F z ( x k ( z ) , y k ( z ) , z ) .
If you include segments parallel to some axes, then for such segments the integrals are equal to zero.
As a result, we can formulate the following definition.
Definition 11.
A general fractional circulation is a general fractional line integral of the vector field F F 1 ( L ) along a piecewise simple closed line L that is defined by
E L ( M ) ( F ) = I L ( M ) [ x , y , z ] F ( x , y , z ) = k = 1 n 1 I [ a k , b k ] ( M 1 ) [ x ] F x ( x , y k ( x ) , z k ( x ) ) + k = 1 n 1 I [ c k , d k ] ( M 2 ) [ y ] F y ( x k ( y ) , y , z k ( y ) ) + k = 1 n 1 I [ e k , f k ] ( M 3 ) [ z ] F z ( x k ( z ) , y k ( z ) , z ) ,
where a n = a 1 , b n = b 1 , c n = c 1 , d n = d 1 , e n = e 1 , f n = f 1 .
Example 10.
Let us consider the piecewise simple closed line L R 0 , + 2 in X Y -plane, which consists of the simple lines A B , B C , C D , and D A , with the coordinates of the points
A ( a , y 1 ( a ) ) , B ( x 1 ( c ) , c ) , C ( b , y 1 ( b ) ) , D ( x ( d ) , d ) ,
where
y 1 ( a ) = y 2 ( a ) , y 1 ( b ) = y 2 ( b ) , x 1 ( c ) = x 2 ( c ) , x 1 ( d ) = x 2 ( d ) .
The lines A B C and A D C are Y-simple lines
A B C : = { ( x , y ) : 0 a x b , y = y 1 ( x ) 0 } ,
A D C : = { ( x , y ) : 0 a x b , y = y 2 ( x ) 0 } ,
where the functions y 1 ( x ) and y 2 ( x ) belong to the space C 1 ( a , b ) , whose derivatives are nonzero and do not change sign on the interval [ a , b ] .
The lines B A D and B C D are X-simple lines
B A D : = { ( x , y ) : 0 c y d , x = x 1 ( y ) 0 } ,
B C D : = { ( x , y ) : 0 c y c , x = x 2 ( y ) 0 } ,
where the functions x 1 ( x ) and x 2 ( x ) belong to the space C 1 ( c , d ) , whose derivatives do not change sign on the interval [ c , d ] .
Let us assume that the functions F x ( x , y ) and F y ( x , y ) satisfy the condition
F x ( x , y 1 ( x ) ) , F x ( x , y 2 ( x ) ) C 1 ( 0 , ) ,
F y ( x 1 ( y ) , y ) , F y ( x ( y ) , y ) C 1 ( 0 , ) ,
Then general fractional circulation along the piecewise simple closed line
L = A B B C C D D A
is defined by the equation
E L ( M ) ( F ) = I L ( M ) [ x , y , z ] F ( x , y , z ) = I L ( M 1 ) [ x ] F x + I L ( M 2 ) [ y ] F y = I A B C ( M 1 ) [ x ] F x ( x , y 1 ( x ) ) I A D C ( M 1 ) [ x ] F x ( x , y 2 ( x ) ) + I B C D ( M 2 ) [ y ] F y ( x 2 ( y ) , y ) I B A D ( M 2 ) [ y ] F y ( x 1 ( y ) , y ) .
Example 11.
Let us consider the piecewise simple line (61), where A n = A 1 , i.e., L is closed line. For example, we consider the line GFI for the rectangle on R 0 , + 2 with vertices at the points
A ( 0 , 0 ) , B ( a , 0 ) , C ( a , b ) , D ( 0 , b ) .
The sides A B , B C , C D , D A of the rectangle form the line L. For closed line L = A B C D , the line GFI operator I L ( M ) is written as
I L ( M ) = e x I A B ( M 1 ) [ x ] + e y I B C ( M 2 ) [ y ] + e x I C D ( M 1 ) [ x ] + e y I D A ( M 2 ) [ y ] = e x I A B ( M 1 ) [ x ] + e y I B C ( M 2 ) [ y ] e x I D C ( M 1 ) [ x ] e y I A D ( M 2 ) [ y ] ,
where we used that
I C D ( M 1 ) [ x ] = I D C ( M 1 ) [ x ] , I D A ( M 2 ) [ y ] = I A D ( M 2 ) [ y ] .
For the vector field
F ( x , y ) = e x F x ( x , y ) + e y F y ( x , y ) ,
the line GFI has the form
I L ( M ) , F = I A B ( M 1 ) [ x ] F x + I B C ( M 2 ) [ y ] F y I D C ( M 1 ) [ x ] F x I A D ( M 2 ) [ y ] F y , I [ 0 , a ] ( M 1 ) [ x ] F x ( x , 0 ) + I [ 0 , b ] ( M 2 ) [ y ] F y ( a , y ) I [ 0 , a ] ( M 1 ) [ x ] F x ( x , b ) I [ 0 , b ] ( M 2 ) [ y ] F y ( 0 , y ) = I [ 0 , a ] ( M 1 ) [ x ] F x ( x , 0 ) F x ( x , b ) + I [ 0 , b ] ( M 2 ) [ y ] F y ( a , y ) F y ( 0 , y ) .
Example 12.
The general fractional circulation for the line L that is rectangle with sides A B , B C , C D , D A , where the points have coordinates (73), is written as
E L ( M ) ( F ) = I L ( M ) , F = I A B ( M 1 ) [ x ] F x + I B C ( M 2 ) [ y ] F y I C D ( M 1 ) [ x ] F z I D A ( M 2 ) [ y ] F z = I ( M 1 ) a [ x ] F x ( x , 0 ) + I ( M 2 ) b [ y ] F y ( a , y ) I ( M 1 ) a [ x ] F x ( x , b ) I ( M 2 ) b [ y ] F y ( 0 , y ) = I ( M 1 ) a [ x ] F x ( x , 0 ) F x ( x , b ) + I ( M 2 ) b [ y ] ( F y ( a , y ) F y ( 0 , y ) ) .
As a result, we get
E L ( M ) ( F ) = 0 a d x M 1 ( a x ) F x ( x , 0 ) F x ( x , b ) + 0 b d y M 2 ( b y ) F y ( a , y ) F y ( 0 , y ) .
For kernels (62), the general fractional circulation (77) has the form
E L ( M ) ( F ) = 0 a d x ( a x ) α 1 1 Γ ( α 1 ) F x ( x , 0 ) F x ( x , b ) + 0 b d y ( b y ) α 2 1 Γ ( α 2 ) F y ( a , y ) F y ( 0 , y ) ,
which was proposed in [31] (see also [13], pp. 241–264), where the Riemann-Liuoville fractional integrals are used. For α 1 = α 2 = 1 , we get the standard circulation.

5. General Fractional Gradient Theorem

5.1. General Fractional Gradient

Let us give definitions of a set of scalar fields and a general fractional gradient for R 0 , + 3 .
Definition 12.
Let U ( x , y , z ) be a scalar field that satisfies the conditions
D ( K 1 ) x , [ x ] U ( x , y , z ) C 1 ( R + 3 ) ,
D ( K 2 ) y , [ y ] U ( x , y , z ) C 1 ( R + 3 ) ,
D ( K 3 ) z , [ z ] U ( x , y , z ) C 1 ( R + 3 ) .
Then the set of such scalar fields will be denoted as C 1 1 ( R + 3 ) .
Definition 13.
Let U ( x , y , z ) be a scalar field that belongs to the set C 1 1 ( R + 3 ) .
Then the general fractional gradient Grad W ( K ) for the region W = R 0 , + 3 is defined as
Grad W ( K ) U ( x , y , z ) =
e x D ( K 1 ) x , [ x ] U ( x , y , z ) + e y D ( K 2 ) y , [ y ] U ( x , y , z ) + e z D ( K 3 ) z , [ z ] U ( x , y , z ) .
This operator will be called the regional general fractional gradient (regional GF gradient).
Remark 11.
The formula defining the operator can be written in compact form. If the scalar field U ( x , y , z ) belongs to the function space C 1 1 ( R 0 , + 3 ) , then the general fractional gradient for the region W = R + 3 is defined as
Grad W ( K ) U ( x , y , z ) = D W ( K ) U ( x , y , z ) = k = 1 3 e k D ( K k ) x k , [ x k ] F k ,
where
D W ( K ) : = e x D ( K 1 ) x , [ x ] + e y D ( K 2 ) y , [ y ] + e z D ( K 3 ) z , [ z ] .
Remark 12.
The general fractional gradients can be defined not only for W = R 0 , + 3 , but also for regions W R 0 , + 3 , surfaces S R 0 , + 3 and line L R 0 , + 3 .
The gradient theorem is very important for the vector calculus and its generalizations, since this theorem is actually a fundamental theorem for standard vector calculus for line integral and gradient, and their generalizations.
In the following sections, we will analyze the differences between regional and line general fractional gradients to formulate general gradient theorems.
Note that the general fractional gradient for line L R 0 , + 3 allows us to prove the general fractional gradient theorem for a wider class of lines R 0 , + 3 and R 0 , + 2 .

5.2. Difficulties in Generalization of Gradient Theorem

Let L be simple line in R 0 , + 3 that is described by the equations y = y ( x ) and z = z ( x ) for x [ a , b ] . Since the line is simple, the derivatives of the functions y = y ( x ) and z = z ( x ) are not equal to zero and do not change sign on the interval [ a ,   b ] . By definition, the linear integral of a vector field F C 1 ( R + 3 ) for the simple line L can be given by the equation
L ( F , d r ) = a b F x ( x , y ( x ) , z ( x ) ) d x +
c d F y ( x ( y ) , y , z ( y ) ) d y + e f F z ( x ( z ) , y ( z ) , z ) d z .
The standard gradient theorem is proved by using this definition of the line integral and the standard chain rule.
For the vector field F = Grad U , we have
L ( Grad U , d r ) = a b U x ( 1 ) ( x , y ( x ) , z ( x ) ) d x +
c d U y ( 1 ) ( x ( y ) , y , z ( y ) ) d y + e f U z ( 1 ) ( x ( z ) , y ( z ) , z ) d z =
c d U y ( 1 ) ( x ( y ) , y , z ( y ) ) d y + e f U z ( 1 ) ( x ( z ) , y ( z ) , z ) d z =
a b U x ( 1 ) ( x , y ( x ) , z ( x ) ) + U y ( 1 ) ( x , y ( x ) , z ( x ) ) y x ( 1 ) ( x ) + U z ( 1 ) ( x , y ( x ) , z ( x ) ) z x ( 1 ) ( x ) d x .
Using the standard chain rule
d U ( x , y ( x ) , z ( x ) ) d x = U ( x , y , z ) x y = y ( x ) , z = z ( x ) +
U ( x , y , z ) y y = y ( x ) , z = z ( x ) d y ( x ) d x + U ( x , y , z ) z y = y ( x ) , z = z ( x ) d z ( x ) d x ,
we get
L ( Grad U , d r ) = a b d U ( x , y ( x ) , z ( x ) ) d x =
U ( b , y ( b ) , z ( b ) ) U ( a , y ( a ) , z ( a ) ) = U ( b , d , f ) U ( a , c , e ) .
In the fractional calculus and GFC, the standard chain rule is violated.
The line GFI for the simple line of the vector field F is defined by the equation
I L ( M ) , F = I [ a , b ] ( M 1 ) [ x ] F x ( x , y ( x ) , z ( x ) ) +
I [ c , d ] ( M 2 ) [ y ] F y ( x ( y ) , y , z ( y ) ) + I [ e , f ] ( M 3 ) [ z ] F z ( x ( z ) , y ( z ) , z ) .
For the vector field F = Grad W ( K ) U , where
F x ( x , y , z ) = ( Grad W ( K ) U ) x ( x , y , z ) = D ( K 1 ) x , [ x ] U ( x , y , z ) ,
F y ( x , y , z ) = ( Grad W ( K ) U ) y ( x , y , z ) = D ( K 2 ) y , [ y ] U ( x , y , z ) ,
F y ( x , y , z ) = ( Grad W ( K ) U ) y ( x , y , z ) = D ( K 3 ) z , [ z ] U ( x , y , z ) ,
we have
I L ( M ) , F = I [ a , b ] ( M 1 ) [ x ] D ( K 1 ) x , [ x ] U ( x , y , z ) y = y ( x ) , z = z ( x ) +
I [ c , d ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] U ( x , y , z ) x = x ( y ) , z = z ( y ) + I [ e , f ] ( M 3 ) [ z ] D ( K 3 ) z , [ z ] U ( x , y , z ) x = x ( z ) , y = y ( z ) .
In this case, we should be emphasized that the fundamental theorem of GFC cannot be used. This fact is based on the following inequalities
I [ a , b ] ( M 1 ) [ x ] D ( K 1 ) x , [ x ] U ( x , y , z ) y = y ( x ) , z = z ( x ) I [ a , b ] ( M 1 ) [ x ] D ( K 1 ) x , [ x ] U ( x , y , z ) y = y ( x ) , z = z ( x ) ,
I [ a , b ] ( M 1 ) [ x ] D ( K 1 ) x , [ x ] U ( x , y , z ) y = y ( x ) , z = z ( x ) I [ a , b ] ( M 1 ) [ x ] D ( K 1 ) x , [ x ] U ( x , y , z ) y = y ( x ) , z = z ( x ) .
which are satisfied if the functions y = y ( x ) , z = z ( x ) depend on the coordinate x. Similarly, the fundamental theorem does not hold for other variables.
As a result, for the regional GF gradient Grad W ( K ) with W = R 0 , + 3 , the general fractional theorem can be proved only if the line consists of sections parallel to the axes.

5.3. General Fractional Gradient Theorem for Regional GF Gradient

The gradient theorem can be considered as a fundamental theorem of standard calculus for line integrals and gradients.
Let us consider a line L that consists of line segments that are parallel to the axes. To prove the general gradient theorem for such a broken line (polygonal chain), it is convenient to use the concept of an elementary broken line. We can state that any continuous line L that consists only of lines parallel to the axes can be represented as a sequence of elementary lines L k .
An elementary broken line is a line consisting of three (no more than three) segments parallel to different axes. There are 48 such elementary lines, eight of which differ in the directions of the segments along the axis or against the axis and six different different orders.
I L k ( M ) [ x , y , z ] : = sgn ( x k + 1 x k ) e x I ω [ x k , x k + 1 ] ( M 1 ) [ x ] + sgn ( y k + 1 y k ) e y I ω [ y k , y k + 1 ] ( M 2 ) [ y ] + sgn ( z k + 1 z k ) e z I ω [ z k , z k + 1 ] ( M 3 ) [ z ] ,
where ω [ x k , x k + 1 ] is defined by (60) and sgn ( x ) is defined by Equation (59).
Let us prove the following general gradient theorem for the regional GF gradient Grad W ( K ) with W = R 0 , + 3 .
Theorem 9.
(General Gradient Theorem for Regional GF Gradient)
Let L be continuous line in R 0 , + 3 that consists only of lines parallel to the axes can be represented as a sequence of elementary lines, for which A k ( a , c , e ) is the initial point and D ( b , d , f ) is the final point.
Let U ( x , y , z ) be a scalar field that belongs to the set C 1 1 ( R + 3 ) .
Then the line GFI for the vector field Grad W ( K ) U with W = R 0 , + 3 satisfies the equation
I L ( M ) , Grad W ( K ) U = U ( b , d , f ) U ( a , c , e ) .
Proof. 
For simplicity, we will consider only an elementary line L k , when moving along which from the initial point A k ( a k , c k , e k ) to the final point D k ( b k , d k , f k ) coordinate values do not decrease, and with order X Y Z . Let the line L k be represented by the points
A k ( a k , c k , e k ) , B k ( b k , c k , e k ) , C k ( b k , d k , e k ) , D ( b k , d k , f k ) .
For the vector field F F 1 , L 0 ( R + 3 ) , for which
F x ( x , c k , e k ) , F y ( b k , y , e k ) , F z ( b k , d k , z ) C 1 ( 0 , ) ,
the line GFI for the elementary line L k = A k B k C k D k is defined by the equation
I L k ( M ) , F =
I A k B k ( M 1 ) [ x ] F x ( x , c k , e k ) + I B k C k ( M 2 ) [ y ] F y ( b k , y , e k ) + I C k D k ( M 3 ) [ z ] F z ( b k , d k , z ) =
I [ a k , b k ] ( M 1 ) [ x ] F x ( x , c k , e k ) + I [ c k , d k ] ( M 2 ) [ y ] F y ( b k , y , e k ) + I [ e k , f k ] ( M 3 ) [ z ] F z ( b k , d k , z ) .
Let us condider the vector field
F = Grad W ( K ) U ,
where W = R 0 , + 3 , and
F x ( x , y , z ) = D ( K 1 ) x , [ x ] U ( x , y , z ) ,
F y ( x , y , z ) = D ( K 2 ) y , [ y ] U ( x , y , z ) ,
F z ( x , y , z ) = D ( K 3 ) z , [ z ] U ( x , y , z ) .
Then line GFI for the elementary line L k = A k B k C k D k is given as
I L k ( M ) , Grad W ( K ) U = I [ a k , b k ] ( M 1 ) [ x ] D ( K 1 ) x , [ x ] U ( x , c k , e k ) +
I [ c k , d k ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] U ( b k , y , e k ) + I [ e k , f k ] ( M 3 ) [ z ] D ( K 3 ) z , [ z ] U ( b k , d k , z ) .
Using the fundamental theorem fo GFC, we obtain
I L k ( M ) , Grad W ( K ) U = U ( b k , c k , e k ) U ( a k , c k , e k ) +
U ( b k , d k , e k ) U ( b k , c k , e k ) ) + U ( b k , d k , f k ) U ( b k , d k , e k ) =
U ( b k , d k , f k ) U ( a k , c k , e k ) .
Equations for other elementary lines are proved in a similar way.
Using that any continuous line L that consists only of lines parallel to the axes can be represented as a sequence of elementary lines L k , we have that general gradient theorem holds. □
Corollary 2.
The line GFI of the vector field F is independent on the path, which is described by the lines L that consist only of line segments parallel to the axes, if the vector field F can be represented as the regional GF gradient of a function U ( x , y , z ) C 1 1 ( R + 3 ) .

5.4. Line General Fractional Gradient in R 0 , + 2

The general fractional gradients can be defined not only for the region W = R 0 , + 3 . Using the fact that GFD is integro-differential operator, we can define general fractional gradients for line L R 0 , + 2 and L R 0 , + 3 .
Let us define a general fractional gradient for a simple line on the X Y -plane. In this definition of line general fractional gradient (line GF Gradienr), we can use the fact that GFD can be represented as a sequential action of a first-order derivative and a general fractional integral:
D ( K 1 ) x , [ x ] = I ( K 1 ) x , [ x ] x .
The general fractional vector operators can also be defined as a sequential action of first-order derivatives and general fractional integrals.
Definition 14.
Let L be a simple line in R 0 , + 2 of the X Y -plane, and U ( x , y ) F 1 , L 1 ( R + 2 ) which means
U x ( 1 ) ( x , y ( x ) ) C 1 ( 0 , ) , U y ( 1 ) ( x ( y ) , y ) C 1 ( 0 , ) .
Then the line general fractional gradient for the line L is defined by the equation
Grad L ( K ) U ( x , y ) = D L ( K ) U ( x , y ) : =
e x I ( K 1 ) x , [ x ] U x ( 1 ) ( x , y ( x ) ) + e y I ( K 2 ) y , [ y ] U y ( 1 ) ( x ( y ) , y ) =
e x I ( K 1 ) x , [ x ] U ( x , y ) x y = y ( x ) + e y I ( K 2 ) y , [ y ] U ( x , y ) y x = x ( y ) ,
where the kernels ( M 1 ( x ) , K 1 ( x ) ) and ( M 2 ( y ) , K 2 ( y ) ) belong to the Luchko set L .
We can use the definitions of the line GFI and line GF Gradient to prove the following theorem.
Theorem 10.
(General Fractional Gradient Theorem for line GF Gradient)
Let L be a simple line in R 0 , + 2 of the X Y -plane, which connects the points A ( a , b ) and B ( c , d ) , and the scalar field U ( x , y ) belongs to the set F 1 , L 1 ( R + 2 ) .
Then, the equality
I L ( M ) , Grad L ( K ) U = U ( c , d ) U ( a , b )
holds.
Proof. 
Using that the line GFI of a vector field F for the simple line L R 0 , + 2 of the X Y -plane is defined by the equation
I L ( M ) , F = I [ a , b ] ( M 1 ) [ x ] F x ( x , y ( x ) ) + I [ c , d ] ( M 2 ) [ y ] F y ( x ( y ) , y ) ,
we can consider the line GFI of the vector field that is defined by the line general fractional gradient
F x ( x , y ( x ) ) : = Grad L ( K ) U ( x , y ) x = I ( K 1 ) x , [ x ] U x ( 1 ) ( x , y ( x ) ) ,
F y ( x ( y ) , x ) : = Grad L ( K ) U ( x , y ) y = I ( K 2 ) y , [ y ] U y ( 1 ) ( x ( y ) , y ) .
Therefore, we can get the line GFI of the line FG Gradient
I L ( M ) , Grad L ( K ) U =
I [ a , b ] ( M 1 ) [ x ] I ( K 1 ) x , [ x ] U x ( 1 ) ( x , y ( x ) ) + I [ c , d ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] U y ( 1 ) ( x ( y ) , y ) .
Then, we can use the fact that the kernels ( M 1 ( x ) , K 1 ( x ) ) and ( M 2 ( y ) , K 2 ( y ) ) belong to the Luchko set L . In this case, we have the property
I [ a , b ] ( M 1 ) [ x ] I ( K 1 ) x , [ x ] f ( x ) = ( M 1 ( K 1 f ) ) ( b ) ( M 1 ( K 1 f ) ) ( a ) =
( ( M 1 K 1 ) f ) ( b ) ( ( M 1 K 1 ) f ) ( a ) = ( { 1 } f ) ( b ) ( { 1 } f ) ( a ) =
0 b f ( x ) d x 0 a f ( x ) d x = a b f ( x ) d x .
Similarly, we obtain the equation
I [ c , d ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] g ( y ) = ( { 1 } g ) ( d ) ( { 1 } g ) ( c ) = c d g ( y ) d y .
Therefore, we get
I [ a , b ] ( M 1 ) [ x ] I ( K 1 ) x , [ x ] U x ( 1 ) ( x , y ( x ) ) + I [ c , d ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] U y ( 1 ) ( x ( y ) , y ) =
a b U ( x , y ) x y = y ( x ) d x + c d U ( x , y ) y x = x ( y ) d y =
a b U ( x , y ) x y = y ( x ) d x + a b U ( x , y ) y y = y ( x ) ) y ( x ) x d x =
a b d U ( x , y ( x ) ) d x d x = U ( c , d ) U ( a , b ) ,
where we use the standard gradient theorem.
As a result, we proved the general fractional gradient theorem for the line GF Gradient with the simple line in R 0 , + 2 . □

5.5. Line General Fractional Gradient in R 0 , + 3

The proposed approach to define line GF Gradien for lines in R 0 , + 2 can be used to define line GF Gradient for lines in R 0 , + 3 . Let L be a simple line in R 0 , + 3 , which is defined in form (46)–(48), and the vector field F ( x , y , z ) belongs to the sset F 1 0 ( L ) . Then, the line GFI over a simple line L in R 0 , + 3 can be defined by Equation (56). The standard line integral of second kind over a simple line L in R 0 , + 3 can be defined by the Equation (55).
Definition 15.
Let L be a simple line in R 0 , + 3 , which is defined in form (46)–(48), and a scalar field U ( x , y , z ) satisfies the conditions
U 1 ( x ) : = U x ( 1 ) ( x , y ( x ) , z ( x ) ) = U ( x , y , z ) x y = y ( x ) , z = z ( x ) C 1 n ( 0 , ) ,
U 2 ( y ) : = U y ( 1 ) ( x ( y ) , y , z ( y ) ) = U ( x , y , z ) y x = x ( y ) , z = z ( y ) C 1 n ( 0 , ) ,
U 3 ( z ) : = U z ( 1 ) ( x ( z ) , y ( z ) , z ) = U ( x , y , z ) z x = x ( z ) , y = y ( z ) C 1 n ( 0 , ) .
Then the set of such scalar fields U ( x , y , z ) will be denoted as U ( x , y , z ) F 1 n ( L ) .
Definition 16.
Let L be a simple line in R 0 , + 3 , which is defined in form (46)–(48), and a scalar field U ( x , y , z ) belongs to the set F 1 1 ( L ) .
Then the line general fractional gradient for the line L R 0 , + 3 is defined by the equation
Grad L ( K ) U ( x , y , z ) = D L ( M ) U ( x , y , z ) = I [ a , b ] ( M 1 ) [ x ] U x ( 1 ) ( x , y ( x ) , z ( x ) ) + I [ c , d ] ( M 2 ) [ y ] U y ( 1 ) ( x ( y ) , y , z ( y ) ) + I [ e , f ] ( M 3 ) [ z ] U z ( 1 ) ( x ( z ) , y ( z ) , z ) .
where the pairs of the kernels ( M 1 ( x ) , K 1 ( x ) ) , ( M 2 ( y ) , K 2 ( y ) ) and ( M 3 ( z ) , K 3 ( z ) ) belong to the Luchko set L .
Let us prove the general fractional gradient theorem for line GF Gradient with simple lines.
Theorem 11.
(General Fractional Gradient Theorem for Line GF Gradient)
Let L be a simple line in R 0 , + 3 , which is defined in form (46)–(48), and connects the points A ( a , c , e ) and B ( b , d , f ) .
Let U ( x , y , z ) be a scalar filed that belongs to the set F 1 , L 1 ( R + 3 ) .
Then, the equality
I L ( M ) , Grad L ( K ) U = U ( a , c , e ) U ( b , d , f )
holds.
Proof. 
Using that the line GFI of a vector field F for the simple line L R 0 , + 3 is defined by the equation
I L ( M ) , F = I [ a , b ] ( M 1 ) [ x ] F x ( x , y ( x ) , z ( x ) ) +
I [ c , d ] ( M 2 ) [ y ] F y ( x ( y ) , y , z ( y ) ) + I [ e , f ] ( M 3 ) [ z ] F z ( x ( z ) , y ( z ) , z ) ,
we can consider the line GFI of the vector field that is defined by the line general fractional gradient
F x ( x , y ( x ) , z ( x ) ) : = Grad L ( K ) U ( x , y , z ) x = I ( K 1 ) x , [ x ] U x ( 1 ) ( x , y ( x ) , z ( x ) ) ,
F y ( x ( y ) , y , z ( y ) ) : = Grad L ( K ) U ( x , y , z ) y = I ( K 2 ) y , [ y ] U y ( 1 ) ( x ( y ) , y . z ( y ) ) .
F z ( x ( z ) , y ( z ) , z ) : = Grad L ( K ) U ( x , y , z ) z = I ( K 3 ) z , [ z ] U z ( 1 ) ( x ( z ) , y ( z ) , z ) .
Therefore, we can get the line GFI of the line FG Gradient
I L ( M ) , Grad L ( K ) U = I [ a , b ] ( M 1 ) [ x ] I ( K 1 ) x , [ x ] U x ( 1 ) ( x , y ( x ) , z ( x ) ) +
I [ c , d ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] U y ( 1 ) ( x ( y ) , y . z ( y ) ) +
I [ e , f ] ( M 3 ) [ z ] I ( K 3 ) z , [ z ] U z ( 1 ) ( x ( z ) , y ( z ) , z ) .
Then, we can use the fact that the kernels ( M 1 ( x ) , K 1 ( x ) ) , ( M 2 ( y ) , K 2 ( y ) ) and ( M 3 ( x ) , K 3 ( x ) ) belong to the Luchko set L . In this case, we have the property
I [ a , b ] ( M 1 ) [ x ] I ( K 1 ) x , [ x ] f ( x ) = ( M 1 ( K 1 f ) ) ( b ) ( M 1 ( K 1 f ) ) ( a ) =
( ( M 1 K 1 ) f ) ( b ) ( ( M 1 K 1 ) f ) ( a ) = ( { 1 } f ) ( b ) ( { 1 } f ) ( a ) =
0 b f ( x ) d x 0 a f ( x ) d x = a b f ( x ) d x .
Similarly, we obtain the equations
I [ c , d ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] g ( y ) = ( { 1 } g ) ( d ) ( { 1 } g ) ( c ) = c d g ( y ) d y ,
I [ e , f ] ( M 3 ) [ z ] I ( K 3 ) z , [ z ] h ( z ) = ( { 1 } h ) ( f ) ( { 1 } h ) ( e ) = e f h ( z ) d z .
Therefore, we get
I [ a , b ] ( M 1 ) [ x ] I ( K 1 ) x , [ x ] U x ( 1 ) ( x , y ( x ) . z ( x ) ) +
I [ c , d ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] U y ( 1 ) ( x ( y ) , y , z ( y ) ) +
I [ e , f ] ( M 3 ) [ z ] I ( K 3 ) z , [ z ] U z ( 1 ) ( x ( z ) , y ( z ) , z ) =
a b U ( x , y , z ) x y = y ( x ) , z = z ( x ) d x + c d U ( x , y , z ) y x = x ( y ) , z = z ( y ) d y +
e f U ( x , y , z ) z x = x ( z ) . y = y ( z ) d z =
a b U ( x , y , z ) x y = y ( x ) , z = z ( x ) d x + a b U ( x , y , z ) y y = y ( x ) , z = z ( x ) y ( x ) x d x +
z b U ( x , y , z ) z y = y ( x ) , z = z ( x ) z ( x ) x d x =
a b d U ( x , y ( x ) , z ( x ) ) d x d x = U ( b , d , f ) U ( a , c , e ) ,
where we use the standard gradient theorem.
As a result, we proved the general fractional gradient theorem for the line GF Gradient with the simple line in R 0 , + 3 . □
Corollary 3.
The line GFI of the vector field F is independent on the path, which is described by the lines L that consist of simple lines, if the vector field F can be represented as the line GF Gradient of a function U ( x , y , z ) F 1 1 ( L ) .
Remark 13.
The general fractional gradient theorem for line GF Gradient can be proves for the Y Z -simple line, which can be considered as a union of X Z -simple lines and X Y -simple lines.

6. General Fractional Green Theorem

It is known that the standard Green’s theorem gives the relationship between a line integral around a simple closed curve S and a double integral over the plane region S bounded by S . Let us formulate a generalization of the Green’s theorem for the general fractional vector calculus.

6.1. General Fractional Green Theorem for X and Y Simple Regions

Let us define a simple region on the plane.
Definition 17.
Let S be a region in R 0 , + 2 of the X Y -plane that is bounded above and below by smooth lines L 1 , x , L 2 , x , and lines L 1 , y , L 2 , y , which are parallel to the Y-axis. We will assume that the lines L 1 , x , L 2 , x are Y-simple lines, which are described by the equations
L 1 , x = { ( x , y ) : x [ a , b ] , y = y 1 ( x ) 0 } ,
L 2 , x = { ( x , y ) : x [ a , b ] , y = y 2 ( x ) 0 } ,
where the functions y 1 ( x ) and y 3 ( x ) belong to the function space C 1 ( a , b ) in the closed domain [ a , b ] that is a projection of the region S onto the X-axis, and y 2 ( x ) y 1 ( x ) for all x [ a , b ] .
The lines L 1 , y and L 2 , y parallel to the Y axis are described as
L 1 , y = { ( x , y ) : x = a 0 , 0 y 1 ( a ) y y 2 ( a ) } ,
L 2 , y = { ( x , y ) : x = b 0 , 0 y 1 ( b ) y y 2 ( b ) } .
If the conditions y 1 ( a ) = y 2 ( a ) , y 1 ( b ) = y 2 ( b ) are satisfied, then the lines L 1 , y and L 2 , y are absent.
Then, such regions S on the plane will be called the Y-simple region (simple area along the y-axis).
Definition 18.
Let S be a region in R 0 , + 2 of the X Y -plane. The region S R 0 , + 2 is called simple, if S is simple along X and Y axes. The region S R 0 , + 2 is called piecewise simple, if S can be divided into a finite number of simple regions with respect to these two axes.
Let us consider the vector field
F ( x , y ) = e x F x ( x , y ) + e y F y ( x , y ) ,
for the simple region S, where the functions F x = F x ( x , y ) and F y = F y ( x , y ) satisfy the conditions
F x ( x , y 1 ( x ) ) , F x ( x , y 2 ( x ) ) C 1 ( 0 , ) ,
F y ( x 1 ( y ) , y ) , F y ( x 2 ( y ) , y ) C 1 ( 0 , ) ,
and
F x ( ξ , y ) , F y ( x , ξ ) C 1 1 ( 0 , ) for each ξ [ 0 , ) .
In this case, we will use notation F ( x , y ) F 1 1 ( S ) for S R + 2 .
Theorem 12.
(General Fractional Green’s Theorem for X and Y Simple Region)
Let S be a Y-simple region in R 0 , + 2 of the X Y -plane that can be described as
S : = { ( x , y ) : 0 a x b , 0 y 1 ( x ) y y 2 ( x ) } ,
where y 1 ( x ) and y 2 ( x ) belong to the space C ( a , b ) .
Let also region S be represented by the union of the finite number of X-simple regions
S = k = 1 n S k
of a finite number X-simple regions
S k : = { ( x , y ) : 0 x 1 , k ( y ) x x 2 , k ( y ) , 0 c k y d k ) } ,
where x 1 , k ( y ) and x 2 , k ( y ) belong to the space C ( a , b ) for all k = 1 , , n .
Let the vector field
F ( x , y ) = e x F x ( x , y ) + e y F y ( x , y ) ,
satisfy the conditions
F x ( x , y 1 ( x ) ) , F x ( x , y 2 ( x ) ) C 1 ( 0 , ) ,
F y ( x 1 , k ( y ) , y ) , F y ( x 2 , k ( y ) , y ) C 1 ( 0 , ) ,
and
F x ( ξ , y ) , F y ( x , ξ ) C 1 1 ( 0 , ) f o r e a c h ξ [ 0 , ) .
Then, the equation
I S ( M ) [ x ] F x ( x , y ) + I S ( M ) [ y ] F y ( x , y ) = I S ( M ) [ x , y ] D ( K 1 ) x , [ x ] F y ( x , y ) D ( K 2 ) y , [ y ] F x ( x , y )
holds.
Proof. 
Let us consider the Y-simple region S on the X Y -plane that is described as
S : = { ( x , y ) : 0 a x b , 0 y 1 ( x ) y y 2 ( x ) } ,
where y 1 ( x ) and y 2 ( x ) are continuous functions on [ a , b ] . The boundary L = S of this area can be divided into the lines A B , B C , C D , D A , where the coordinates of the points are the following:
A ( a , y 1 ( a ) ) , B ( b , y 1 ( b ) ) , C ( b , y 2 ( b ) ) , D ( a , y 2 ( a ) ) .
The lines A B and C D are Y-simple lines ( L 1 x = A B , L 2 x = C D ), and B C , D A are lines parallel to the Y axis.
Let us consider the vector field F as a sum of two vectors F x and F y in the form
F ( x , y ) = F x ( x , y ) + F y ( x , y ) : F x = e x F x ( x , y ) , F y = e y F y ( x , y ) ,
where the functions F x ( x , y ) and F r ( x , y ) satisfy conditions (87)–(89).
(1)
The line GFI for the vector field F x = e x F x ( x , y ) and the line L = S is given by the equation
I L ( M ) [ x , y ] , F x ( x , y ) = I A B ( M 1 ) [ x ] F x ( x , y 1 ( x ) ) +
I B C ( M 1 ) [ x ] F x ( b , y ) + I C D ( M 1 ) [ x ] F x ( x , y 2 ( x ) ) + I D A ( M 1 ) [ x ] F x ( a , y ) =
I A B ( M 1 ) [ x ] F x ( x , y 1 ( x ) ) I D C ( M 1 ) [ x ] F x ( x , y 2 ( x ) ) =
I [ a , b ] ( M 1 ) [ x ] F x ( x , y 1 ( x ) ) F x ( x , y 2 ( x ) ) ,
where we use
I B C ( M 1 ) [ x ] F x ( b , y ) = 0 , I D A ( M 1 ) [ x ] F x ( a , y ) = 0 ,
since d x = 0 for x = b and x = a .
Using the second fundamental theorem of GFC in the form
F x ( x , y 1 ( x ) ) F x ( x , y 2 ( x ) ) = I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] F x ( x , y ) ,
where we assume that conditions (87) and (89) are satisfied. Then, we obtain
I L ( M ) [ x , y ] , F x ( x , y ) = I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] F x ( x , y ) .
As a result, we obtain
I L ( M ) [ x , y ] , F x ( x , y ) = I S ( M ) [ x , y ] D ( K 2 ) y , [ y ] F x ( x , y ) ,
where we use the definition of double GFI by the iterated (repeated) GFI
I S ( M ) [ x , y ] = I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] .
(2)
Let us consider the vector field F y = e y F y ( x , y ) and the Y-simple region S on the X Y -plane, which can be described by (85) and (86). Using the similar transformations and equations for the line GFI for L = S , we obtain
I L ( M ) [ x , y ] , F y ( x , y ) = k = 1 n l = 1 m I L k l ( M ) [ x , y ] , F y ( x , y ) =
k = 1 n l = 1 m I [ c k l , d k l ] ( M 2 ) [ y ] I [ x 1 , k l ( y ) , x 2 , k l ( y ) ] ( M 1 ) F y ( x , y ) = I S ( M ) [ x , y ] D ( K 1 ) x , [ x ] F y ( x , y ) ,
where we use the definition of double GFI by the iterated (repeated) GFI
I S ( M ) [ x , y ] = k = 1 n l = 1 m I S k l ( M ) [ x , y ] = k = 1 n l = 1 m I [ c k l , d k l ] ( M 2 ) [ y ] I [ x 1 , k l ( y ) , x 2 , k l ( y ) ] ( M 1 ) [ x ] .
(3)
Therefore, the Green formula for the vector field F = F x + F y takes the form
I L ( M ) [ x , y ] , F ( x , y ) = I S ( M ) [ x , y ] D ( K 1 ) x , [ x ] F y ( x , y ) D ( K 2 ) y , [ y ] F x ( x , y ) ,
if conditions (87)–(89) are satisfied. □
Remark 14.
Using the proof of the theorem, we see that it is not required to use an approach similar to the approach that is proposed for the linearly general fractional gradient. WE should also note that this approach should be used for the general fractional Sokes theorem.
Remark 15.
Using the proof, we can see that the general fractional Green’s theorem can be proved for the region S that is union of finite number of Y-simple regions S j ( j = 1 , , m ), each of which can be represented as a union of a finite number of simple regions S j k , k = 1 , , n .

6.2. General Fractional Green Theorem for Simple Region

Let us consider a simple region S R 0 , + 2 of the X Y -plane that can be described as
S : = { ( x , y ) : 0 a x b , 0 y 1 ( x ) y y 2 ( x ) } ,
S : = { ( x , y ) : 0 x 1 ( y ) x x 2 ( y ) , 0 c y d ) } ,
where y 1 ( x ) and y 2 ( x ) belong to the space C 1 ( a , b ) , the derivatives of these functions are nonzero and do not change sign on the interval [ a ,   b ] . Then the functions x 1 ( y ) and x 2 ( y ) exit, belong to the space C 1 ( c , d ) , where c = x 1 ( a ) = x 2 ( a ) and d = x 1 ( b ) = x 2 ( b ) . In this case, we will use notation S G 1 ( R 0 , + 2 ) .
The boundary of the simple region S G 1 ( R 0 , + 2 ) is a piecewise simple closed line S in the X Y -plane that can be described as the union of two Y-simple lines
S : = { ( x , y ) : x [ a , b ] , y = y 1 ( x ) } { ( x , y ) : x [ a , b ] , y = y 2 ( x ) } ,
or as the union of two X-simple lines:
S : = { ( x , y ) : y [ c , d ] , x = x 1 ( y ) } { ( x , y ) : y [ c , d ] , x = x 2 ( y ) } .
For S G 1 ( R 0 , + 2 ) , we will consider the vector field
F ( x , y ) = e x F x ( x , y ) + e y F y ( x , y ) ,
where the functions F x = F x ( x , y ) and F y = F y ( x , y ) satisfy the conditions
F x ( x , y 1 ( x ) ) , F x ( x , y 2 ( x ) ) C 1 ( 0 , ) ,
F y ( x 1 ( y ) , y ) , F y ( x 2 ( y ) , y ) C 1 ( 0 , ) ,
and
F x ( ξ , y ) , F y ( x , ξ ) C 1 1 ( 0 , ) for each ξ [ 0 , ) .
In this case, we will use notation F ( x , y ) F 1 ( S ) for S R + 2 .
Theorem 13.
(General Fractional Green’s Theorem for Simple Region)
Let S be a simple region in R 0 , + 2 of the X Y -plane such that S G 1 ( R 0 , + 2 ) , and the vector field F [ x , y ] belongs to the set F 1 1 ( R + 2 ) .
Then, the equation
I S ( M ) [ x ] F x ( x , y ) + I S ( M ) [ y ] F y ( x , y ) =
I S ( M ) [ x , y ] D ( K 1 ) x , [ x ] F y ( x , y ) D ( K 2 ) y , [ y ] F x ( x , y )
are satisfied.
Proof. 
Let us consider the piecewise simple closed line S R 0 , + 2 in X Y -plane, which consists of the simple lines A B , B C , C D and D A ,
S = A B B C C D D A
with the coordinates of the points
A ( a , y 1 ( a ) ) , B ( x 1 ( c ) , c ) , C ( b , y 1 ( b ) ) , D ( x 1 ( d ) , d ) ,
where
y 1 ( a ) = y 2 ( a ) , y 1 ( b ) = y 2 ( b ) , x 1 ( c ) = x 2 ( c ) , x 1 ( d ) = x 2 ( d ) .
The lines A B C and A D C are Y-simple lines
A B C : = { ( x , y ) : 0 a x b , y = y 1 ( x ) 0 } ,
A D C : = { ( x , y ) : 0 a x b , y = y 2 ( x ) 0 } ,
where the functions y 1 ( x ) and y 2 ( x ) belong to the space C 1 ( a , b ) , and derivatives of these functions are nonzero and do not change sign on the interval [ a , b ] .
The lines B A D and B C D are X-simple lines
B A D : = { ( x , y ) : 0 c y d , x = x 1 ( y ) 0 } ,
B C D : = { ( x , y ) : 0 c y c , x = x 2 ( y ) 0 } ,
where the functions x 1 ( x ) and x 2 ( x ) belong to the space C 1 ( c , d ) , and derivatives of these functions do not change sign on the interval [ c , d ] .
Using that the functions F x ( x , y ) and F y ( x , y ) satisfy the condition
F x ( x , y 1 ( x ) ) , F x ( x , y 2 ( x ) ) C 1 ( 0 , ) ,
F y ( x 1 ( y ) , y ) , F y ( x ( y ) , y ) C 1 ( 0 , ) ,
the general fractional circulation along the piecewise simple closed line S is described by the equation
I L ( M ) [ x , y ] F ( x , y ) =
I A B C ( M 1 ) [ x ] F x ( x , y 1 ( x ) ) + I C D A ( M 1 ) [ x ] F x ( x , y 2 ( x ) ) +
I B C D ( M 2 ) [ y ] F y ( x 2 ( y ) , y ) + I D A B ( M 2 ) [ y ] F y ( x 1 ( y ) , y ) =
I A B C ( M 1 ) [ x ] F x ( x , y 1 ( x ) ) I A D C ( M 1 ) [ x ] F x ( x , y 2 ( x ) ) +
I B C D ( M 2 ) [ y ] F y ( x 2 ( y ) , y ) I B A D ( M 2 ) [ y ] F y ( x 1 ( y ) , y ) =
I [ a , b ] ( M 1 ) [ x ] F x ( x , y 1 ( x ) ) I [ a , b ] ( M 1 ) [ x ] F x ( x , y 2 ( x ) ) +
I [ c , d ] ( M 2 ) [ y ] F y ( x 2 ( y ) , y ) I [ c , d ] ( M 2 ) [ y ] F y ( x 1 ( y ) , y ) =
I [ a , b ] ( M 1 ) [ x ] F x ( x , y 1 ( x ) ) F x ( x , y 2 ( x ) ) +
I [ c , d ] ( M 2 ) [ y ] F y ( x 2 ( y ) , y ) F y ( x 1 ( y ) , y ) =
I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] F x ( x , y )
I [ c , d ] ( M 2 ) [ y ] I [ x 1 ( y ) , x 2 ( x ) ] ( M 1 ) [ x ] D ( K 1 ) x , [ x ] F y ( x , y ) =
I [ S ] ( M ) [ x , y ] D ( K 1 ) x , [ x ] F y ( x , y ) D ( K 2 ) y , [ y ] F x ( x , y ) ,
where we used the second fundamental theorem of general fractional calculus. □

6.3. General Fractional Green Theorem for Rectangle

Let us formulate the general fractional Green’s theorem for the region in the form of a rectangle.
Theorem 14.
(General Fractional Green’s Theorem for Rectangle)
Let functions F x = F x ( x , y ) and F y = F y ( x , y ) belong to the function space S 1 1 ( R + 2 ) , and S is the rectangle
S : = { ( x , y ) : x [ a , b ] , y [ c , d ] } ,
and the boundary of S be the closed line S . Then
I S ( M ) [ x ] F x ( x , y ) + I S ( M ) [ y ] F y ( x , y ) =
I S ( M ) [ x , y ] D ( K 2 ) y , [ y ] F x ( x , y ) D ( K 1 ) x , [ x ] F y ( x , y ) .
Proof. 
(1) To prove Equation (108), the double fractional integral I S ( M ) [ x , y ] is written as the iterated (repeated) fractional integrals and then the fundamental theorems of general fractional calculus is used.
Let S be the rectangular domain (107) with vertexes in the points
A ( a , c ) , B ( b , c ) , C ( b , d ) , D ( a , d ) .
The sides A B , B C , C D , D A of the rectangular domain (107) form the boundary S of S. For the rectangular region S defined by a x b , and c y d , the iterated (repeated) integral is
I S ( M ) [ x , y ] = I [ a , b ] ( M 1 ) [ x ] I [ c , d ] ( M 2 ) [ y ] ,
where
I [ a , b ] ( M 1 ) [ x ] = I ( M 1 ) b [ x ] I ( M 1 ) a [ x ] , I [ c , d ] ( M 2 ) [ y ] = I ( M 2 ) d [ y ] I ( M 2 ) c [ y ] .
(2) To prove the general fractional Green’s formula, we realize the following transformations
I S ( M ) [ x ] F x + I S ( M ) [ y ] F y = I A B ( M 1 ) [ x ] F x + I C D ( M 1 ) [ x ] F x + I B C ( M 2 ) [ y ] F y + I D A ( M 2 ) [ y ] F y = I [ a , b ] ( M 1 ) [ x ] F x ( x , c ) I [ a , b ] ( M 1 ) [ x ] F x ( x , d ) + I [ c , d ] ( M 2 ) [ y ] F y ( b , y ) I [ c , d ] ( M 2 ) [ y ] F y ( a , y ) = I [ a , b ] ( M 1 ) [ x ] F x ( x , c ) F x ( x , d ) + I [ c , d ] ( M 2 ) [ y ] F y ( b , y ) d y F y ( a , y ) .
Let us use the second fundamental theorem of GFC in the from
F x ( x , c ) F x ( x , d ) = I [ c , d ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] F x ( x , y ) ,
F y ( b , y ) F y ( a , y ) = I [ a , b ] ( M 1 ) [ x ] D ( K 1 ) x , [ x ] F y ( x , y ) .
Then, expression (110) can be written as
I [ a , b ] ( M 1 ) [ x ] I [ c , d ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] F x ( x , y ) + I [ c , d ] ( M 2 ) [ y ] I [ a , b ] ( M 1 ) [ x ] D ( K 1 ) x , [ x ] F y ( x , y ) = I [ a , b ] ( M 1 ) [ x ] I [ c , d ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] F x ( x , y ) D ( K 1 ) x , [ x ] F y ( x , y ) = I S ( M ) [ x , y ] D ( K 2 ) y , [ y ] F x ( x , y ) D ( K 1 ) x , [ x ] F y ( x , y ) ,
where we use (109) for the operator I S ( M ) [ x , y ] . As a result, we obtain the right-hand side of Equation (108). □
Remark 16.
For the kernels
M 1 ( x ) = h α 1 ( x ) , M 2 ( y ) = h α 2 ( y ) , K 1 ( x ) = h 1 α 1 ( x ) , K 2 ( y ) = h 1 α 2 ( y ) ,
a general fractional Green’s theorem was proved in [31], where the Riemann-Liuoville fractional integrals and the Caputo fractional derivatives are used. For α 1 = α 2 = 1 , we get the standard Green’s theorem.

7. Double and Surface General Fractional Integral, and Flux

7.1. Definition of Double GFI by Iterated (Repeated) GFI

Let us consider a region S in the first quarter of the coordinate plane R + 2 , and assume that the area S is bounded such that
S : = { ( x , y ) : 0 a x b , 0 y 1 ( x ) y y 2 ( x ) } ,
where y = y 1 ( x ) and y = y 2 ( x ) are continuous functions, and the following two requirements are satisfied: (1) the projection of S onto the X-axis is bounded by the two values x = a and x = b ; (2) any line parallel to the Y-axis that passes between these two values x = a and x = b intersects the line y = y 1 ( x ) and y = y 2 ( x ) in the interval ( a , b ) at no more than one point. Then the region S will be called the Y-simple domain.
Let us define the double GFI through the iterated (repeated) GFI.
Definition 19.
If the function f = f ( x , y ) C 1 ( R + 1 2 ) satisfies the condition
J ( x ) : = I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] f ( x , y ) C 1 ( 0 , ) ,
where
I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] f ( x , y ) : = I ( M 2 ) y 2 ( x ) [ y ] f ( x , y ) I ( M 2 ) y 1 ( x ) [ y ] f ( x , y ) ,
and S is the Y-simple domain (114), then the iterated (repeated) GFI of the form
I S ( M ) [ x , y ] : = I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] ,
is called the double GFI, where
I [ a , b ] ( M 1 ) [ x ] = I ( M 1 ) b [ x ] I ( M 1 ) a [ x ] .
Using the function
F ( x , y ) = I ( M 2 ) y [ y ] f ( x , y ) = 0 y M 2 ( y y ) f ( x , y ) d y ,
Equation (115) can be written in the form
I S ( M ) [ x , y ] f ( x , y ) : = 0 b M 1 ( b x ) F ( x , y 2 ( x ) ) F ( x , y 1 ( x ) ) d x 0 a M 1 ( a x ) F ( x , y 2 ( x ) ) F ( x , y 1 ( x ) ) d x .
Example 13.
For the rectangular area
S : = { ( x , y ) : 0 a x b , 0 c y d } ,
double GFI (115) has the form
I S ( M ) [ x , y ] : = I [ a , b ] ( M 1 ) [ x ] I [ c , d ] ( M 2 ) [ y ] .
If the function f = f ( x , y ) C 1 ( R + 1 2 ) and a = c = 0 , then
I S ( M ) [ x , y ] f ( x , y ) = I ( M 1 ) b [ x ] I ( M 2 ) d [ y ] f ( x , y ) = 0 b 0 d M 1 ( b x ) M 2 ( d y ) f ( x , y ) d x d y .
Example 14.
For the area
S : = { ( x , y ) : 0 x b , 0 y g ( x ) } ,
and the function f = f ( x , y ) C 1 ( R + 1 2 ) , the double GFI has the form
I S ( M ) [ x , y ] f ( x , y ) : = I ( M 1 ) b [ x ] I ( M 2 ) g ( x ) [ y ] f ( x , y ) = 0 b 0 g ( x ) M 1 ( b x ) M 2 ( g ( x ) y ) f ( x , y ) d x d y .
Remark 17.
Let us consider double GFI (115) for the function
f ( x , y ) = D ( K 2 ) y , [ y ] F x ( x , y ) = 0 y d y K 2 ( y y ) F x ( x , y ) ,
where D ( K 2 ) y , [ y ] is the GFD with respect to the variable y with kernel K 2 ( y ) that is associated to the kernel M 2 ( y ) . Then
I S ( M ) [ x , y ] f ( x , y ) = I ( M 1 ) a [ x ] I ( M 2 ) g ( x ) [ y ] D ( K 2 ) y , [ y ] F x ( x , y ) .
For area (120), we obtain
I S ( M ) [ x , y ] D ( K 2 ) y , [ y ] F x ( x , y ) = I ( M 1 ) a [ x ] I ( M 2 ) g ( x ) [ y ] D ( K 2 ) y , [ y ] F x ( x , y ) =
I ( M 1 ) a [ x ] F 1 ( x , g ( x ) ) F x ( x , 0 ) ,
where the second fundamental theorem of GFC is used. In the general case, we obtain
I S ( M ) [ x , y ] D ( K 2 ) y ) [ y ] F x ( x , y ) = I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] F x ( x , y ) =
I [ a , b ] ( M 1 ) [ x ] F x ( x , y 2 ( x ) ) F x ( x , y 1 ( x ) ) ,
0 b d x M 1 ( b x ) F x ( x , y 2 ( x ) ) F x ( x , y 1 ( x ) )
0 a d x M 1 ( a x ) F x ( x , y 2 ( x ) ) F x ( x , y 1 ( x ) ) ,
where y 2 ( x ) y 1 ( x ) > 0 for all 0 x b (instead of a x b ). Here we assume that
F x ( ξ , y ) C 1 1 ( 0 , ) for all ξ ( 0 , ) ,
F x ( x , y 1 ( x ) ) , F x ( x , y 2 ( x ) ) C 1 ( 0 , ) .

7.2. Surface General Fractional Integral

The standard surface integral of integer order is a generalization of double integrals to integration over surfaces. This integral can be considered as the double integral analogue of the line integral.
The surface general fractional integral over the surface S in R 0 , + 3 can be defined through the double general fractional integrals over areas S y z , S x z , S x y in the Y Z , X Z , X Y planes, where these areas are projections of the surface S onto these planes.
Let us consider a two-sided smooth surface S, and fix one of its two sides, which is equivalent to choosing a certain orientation on the surface. We also assume that the surface is given by the equation
x = x ( y , z ) 0 ,
where the point ( y , z ) changes in area S y z in the Y Z -plane, bounded by smooth contour S y z , and
S : = { ( x , y , z ) : x = x ( y , z ) 0 , ( y , z ) S y z R 0 , + 2 } .
Let us first consider a vector field of the form
F x : = e x F x ( x , y , z ) .
If the function F x ( x , y , z ) on the surface S belongs to the function space C 1 ( R + 2 ) that is
F x ( x ( y , z ) , y , z ) C 1 ( R + 2 ) ,
then the surface general fractional integral is defined by the equation
I S ( M ) , F x = I S y z ( M ) [ y , z ] F x ( x ( y , z ) , y , z ) .
Remark 18.
If the area S y z in the Y Z -plane is given as
S y z : = { ( y , z ) : 0 c y d , 0 z 1 ( y ) z z 2 ( y ) } ,
where z = z 1 ( y ) and z = z 2 ( y ) are continuous functions, then the double GFI in Equation (124) is defined as
I S y z ( M ) [ y , z ] F x ( x ( y , z ) , y , z ) : = I [ c , d ] ( M 2 ) [ z ] I [ z 1 ( y ) , z 2 ( y ) ] ( M 3 ) [ y ] F x ( x ( y , z ) , y , z ) .
Similar to vector field (122), we can obtain the surface GFIs for the other two projections of the surface S on the X Z and Y Z planes, and functions F y ( x , y , z ) , F z ( x , y , z ) on the surface S. If, instead of the Y Z -plane, we project the surface S onto the X Z -plane or Y Z -plane, then we get two other surface GFI of the second type:
I S ( M ) , F y = I S x z ( M ) [ x , z ] F y ( x , y ( x , z ) , z ) ,
I S ( M ) , F z = I S x y ( M ) [ x , y ] F z ( x , y , z ( x , y ) ) ,
where
F y ( x , y ( x , z ) , z ) C 1 ( R + 2 ) , F z ( x , y , z ( x , y ) ) C 1 ( R + 2 ) .
In this case, the surface GFI is defined by projections onto all three planes ( Y Z , X Z , X Y ) in the form
( I S ( M ) , F ) = I S y z ( M ) [ y , z ] F x ( x ( y , z ) , y , z ) +
I S x z ( M ) [ x , z ] F y ( x , y ( x , z ) , z ) + I S x y ( M ) [ x , y ] F z ( x , y , z ( x , y ) ) ,
where conditions (123) and (125) are satisfied.
Let us define a set of piecewise simple surfaces, a set of vector fields on these surfaces and the surface GFI.
Definition 20.
Let S be an oriented compact smooth surface in the region
R 0 , + 3 = { ( x , y , z ) : x 0 , y 0 , z 0 } .
Let us choose a side of the surface S.
Let the surface S be represented as the union of a finite number of X-simple surfaces S X , i , as well as a finite number of Y-simple surfaces S X , j , and Z-simple surfaces S Z , k , where i = 1 , , n x , j = 1 , , n y , k = 1 , , n z :
S : = k = 1 n x S X , i = j = 1 n y S Y , j = k = 1 n z S Z , k .
Let S i , y z , S j , x z , S k , x y be projections of the surfaces S X , i , S Y , j , S Z , k onto the Y Z , X Z , X Y planes, which can be described by continuous functions
S X , i = { x = x i ( y , z ) 0 i f ( y , z ) S i , y z R 0 , + 2 } ,
S Y , j = { y = y j ( x , z ) 0 i f ( x , z ) S j , x z R 0 , + 2 } ,
S Z , k = { z = z k ( x , y ) 0 i f ( x , y ) S k , x y R 0 , + 2 } .
Such surfaces S will be called piecewise simple surfaces. The set of such surfaces will be denoted by P ( R 0 , + 3 ) .
Definition 21.
Let S be a piecewise simple surface ( S P ( R 0 , + 3 ) ). Let the vector field
F : = e x F x ( x , y , z ) + e y F y ( x , y , z ) + e z F z ( x , y , z )
on the surface S satisfy the conditions
F x ( x i ( y , z ) , y , z ) C 1 ( R + 2 ) ,
F y ( x , y j ( x , z ) , z ) C 1 ( R + 2 ) ,
F z ( x , y , z k ( x , y ) ) C 1 ( R + 2 )
for all i = 1 , , n x , j = 1 , , n y , k = 1 , , n z .
The set of such vector fields F on piecewise simple surface S will be denoted by F S ( R 0 , + 3 ) .
Definition 22.
Let S be a piecewise simple surface ( S P ( R 0 , + 3 ) ) and a vector fields F on this surface S belongs to the set F S ( R 0 , + 3 ) .
Then the surface general fractional vector integral (surface GFI) of the second kind
I S ( M ) = e x I S y z ( M ) [ y , z ] + e y I S x z ( M ) [ z , x ] + e z I S x y ( M ) [ x , y ] = e x i = 1 n x I S i , y z ( M ) [ y , z ] + e y j = 1 n y I S j , x z ( M ) [ z , x ] + e z k = 1 n z I S k , x y ( M ) [ x , y ]
for the vector fiels F F S ( R 0 , + 3 ) is defined by the equation
( I S ( M ) , F ) : = i = 1 n x I S i , y z ( M ) [ y , z ] F x ( x i ( y , z ) , y , z ) +
j = 1 n y I S j , x z ( M ) [ x , z ] F y ( x , y j ( x , z ) , z ) + k = 1 n z I S k , x y ( M ) [ x , y ] F z ( x , y , z k ( x , y ) ) .

7.3. General Fractional Flux

Let us consider the general fractional flux for simple surfaces.
Definition 23.
Let S be an oriented compact smooth simple surface in the region R 0 , + 3 such that S y z , S x z , S x y are the projections of the surface S onto the Y Z , X Z , X Y planes, which can be described by continuous functions
S = { x = x ( y , z ) 0 i f ( y , z ) S y z R 0 , + 2 } ,
S = { y = y ( x , z ) 0 i f ( x , z ) S x z R 0 , + 2 } ,
S = { z = z ( x , y ) 0 i f ( x , y ) S x y R 0 , + 2 } .
Then a general fractional flux of the vector field F F S ( R 0 , + 3 ) across the simple surface S is the surface general fractional integral of the form
Φ S ( M ) ( F ) = I S ( M ) , F = I S y z ( M ) [ y , z ] F x ( x ( y , z ) , y , z ) + I S x z ( M ) [ x , z ] F y ( x , y ( x , z ) , z ) + I S x y ( M ) [ x , y ] F z ( x , y , z ( x , y ) ) ,
where
F x ( x ( y , z ) , y , z ) , F y ( x , y ( x , z ) , z ) , F z ( x , y , z ( x , y ) ) C 1 ( R + 2 ) .
Example 15.
Let us consider a surface S in the form of the rectangle
S = S y z : = { ( y , z ) : x = 0 , 0 c y d , 0 e z f } .
The general fractional flux of the vector field F across this surface is written as
Φ S ( M ) ( F ) = I S y z ( M ) [ y , z ] F x ( 0 , y , z ) = I [ c , d ] ( M 2 ) [ y ] I [ e , f ] ( M 3 ) [ z ] F x ( 0 , y , z ) .
Example 16.
Let us consider a surface S in the form of boundary W of the parallelepiped area
W : = { { ( x , y , z ) : 0 a x b , 0 c y d , 0 e z f } .
The general fractional flux of the vector field
F = e x F x ( x , y , z ) + e y F y ( x , y , z ) + e z F z ( x , y , z )
across the surface S = W is written as
Φ S ( M ) ( F ) = I S ( M ) , F =
I S y z ( M ) [ y , z ] F x ( x , y , z ) + I S x z ( M ) [ x , z ] F y ( x , y , z ) + I S x y ( M ) [ x , y ] F z ( x , y , z ) =
I [ c , d ] ( M 2 ) [ y ] I [ e , f ] ( M 3 ) [ z ] F x ( b , y , z ) I [ c , d ] ( M 2 ) [ y ] I [ e , f ] ( M 3 ) [ z ] F x ( a , y , z ) +
I [ a , b ] ( M 1 ) [ x ] I [ e , f ] ( M 3 ) [ z ] F y ( x , c , z ) I [ a , b ] ( M 1 ) [ x ] I [ e , f ] ( M 3 ) [ z ] F y ( x , d , z ) +
I [ a , b ] ( M 1 ) [ x ] I [ c , d ] ( M 2 ) [ y ] F z ( x , y , f ) I [ a , b ] ( M 1 ) [ x ] I [ c , d ] ( M 2 ) [ y ] F z ( x , y , e ) ,
where
I [ a k , b k ] ( M k ) [ x k ] : = I ( M k ) b k [ x k ] I ( M k ) a k [ x k ] .
Example 17.
For the rectangular area
S : = { 0 x a , 0 y b , z = 0 } ,
the general fractional flux of the vector field F across surface (132) is written as
Φ S ( M ) ( F ) = I S x y ( M ) [ x , y ] F z ( x , y , 0 ) = I S x y ( M 1 ) [ x ] I S x y ( M 2 ) [ y ] F z ( x , y , 0 ) = I ( M 1 ) a [ x ] I ( M 2 ) b [ y ] F z ( x , y , 0 ) = 0 a d x 0 b d y M 1 ( a x ) M 2 ( b y ) F z ( x , y , 0 ) .
For kernels (62), the general fractional flux (133) for surface (132) has the form
Φ S ( M ) ( F ) = 0 a d x 0 b d y ( a x ) α 1 1 Γ ( α 1 ) ( b y ) α 2 1 Γ ( α 2 ) F z ( x , y , 0 ) .
For α 1 = α 2 = 1 , the general fractional flux (134) for surface (132) is the standard flux
Φ S ( F ) = 0 a d x 0 b d y F z ( x , y , 0 ) .

8. General Fractional Stokes Theorem

The Stokes theorem connects the surface integral with the line integral. The Stokes theorem generalizes Green theorem from R 2 to R 3 . If the surface is a flat region lying in the plane, then the Stokes equation gives the Green equation.

8.1. Simple Domain on X Y -Plane

We recall the definitions of simple domains.
Definition 24.
The closed domain D x y R 0 , + 2 on the X Y -plane will be called Y-simple domain, if D x y can be represented in the form
D x y : = { ( x , y ) : 0 a x b , 0 y 1 ( x ) y y 2 ( x ) } ,
where y = y 1 ( x ) , y = y 2 ( x ) are continuous functions for x [ a , b ] .
Similarly, the closed domain D x y on the X Y -plane is called X-simple domain, if D x y can be represented as
D x y : = { ( x , y ) : 0 x 1 ( y ) x x 2 ( y ) , 0 c y d } ,
where x = x 1 ( y ) , x = x 2 ( y ) are continuous functions for y [ c , d ] .
Definition 25.
If D x y is X-simple and Y-simple domain, then D x y is called the simple domain on the X Y -plane.
The simple domain D x y R 0 , + 2 on the X Y -plane can be represented in the form
D x y : = { ( x , y ) : 0 a x b , 0 y 1 ( x ) y y 2 ( x ) } ,
D x y : = { ( x , y ) : 0 x 1 ( y ) x x 2 ( y ) , 0 c y d } ,
where y 1 ( x ) , y 2 ( x ) C [ a , b ] and x 1 ( y ) , x 2 ( y ) C [ c , d ] , and y 1 ( a ) = y 2 ( a ) = c , y 1 ( b ) = y 2 ( b ) = d , x 1 ( c ) = x 2 ( c ) = a , x 1 ( d ) = x 2 ( d ) = b .
The boundary of a Y-simple closed domain D x y in R 0 , + 2 can be represented as a closed line L x y = D x y consisting of the lines
L y 1 : = { ( x , y ) : 0 a x b , y = y 1 ( x ) 0 } ,
L y 2 : = { ( x , y ) : 0 a x b , y = y 2 ( x ) 0 } ,
L x 1 : = { ( x , y ) : x = a 0 0 y 1 ( a ) y y 2 ( a ) } ,
L x 2 : = { ( x , y ) : x = b 0 0 y 1 ( b ) y y 2 ( b ) } .
To simplify the proofs, we can sometimes use the case, when there are no straight lines B a , B b that is, when y 1 ( a ) = y 2 ( a ) and y 1 ( b ) = y 2 ( b ) .
Similarly, we obtain a representation of the boundary for a X-simple domian D x y .

8.2. Simple Surface

Let us give a definition of Z-simple surface.
Definition 26.
Let S be a smooth oriented surface in R 0 , + 3 that the surface S is described by the equation
z = z ( x , y ) 0 ,
where the function z ( x , y ) is continuous in the closed domain D x y ( z ( x , y ) C ( D x y ) ), which is a projection of the surface S onto the X Y -plane.
We will assume that is bounded by an oriented closed smooth line L = S . The boundary D x y of the domain D x y is a projection of the line L = S onto the X Y -plane.
Then the surface S will be called the Z-simple surface, if the domain D x y is simple on the X Y -plane.
Definition 27.
The surface S will be called the simple surface in R 0 , + 3 , if S is simple with respect to the X, Y and Z axes.
The simple surface S can be described by the continuous functions
S : = { ( x , y , z ) : x = x ( y , z ) C ( R + 2 ) , ( y , z ) D y z R 0 , + 2 } ,
S : = { ( x , y , z ) : y = y ( x , z ) C ( R + 2 ) , ( x , z ) D x z R 0 , + 2 } ,
S : = { ( x , y , z ) : z = z ( x , y ) C ( R + 2 ) , ( x , y ) D x , y R 0 , + 2 } ,
where the domains D y z , D x z , D x y are simple domains in R 0 , + 2 , which are projections of the surface S onto the Y Z , X Z , X Y planes.
Remark 19.
We can consider the Z-simple surfaces S in R 0 , + 3 that can be described as a union of X-simple surfaces, as well as a union of Y-simple regions.

8.3. Vector Field on Surface

Let us define the properties of the scalar and vector fields in simple domains in R 0 , + 2 .
Definition 28.
Let D x y R 0 , + 2 be a simple closed domain on the X Y -plane that can be described in the forms
D x y : = { ( x , y ) : 0 a x b , 0 y 1 ( x ) y y 2 ( x ) } ,
D x y : = { ( x , y ) : 0 x 1 ( y ) x x 2 ( y ) , 0 c y d } ,
where y 1 ( x ) , y 2 ( x ) C [ a , b ] and x 1 ( y ) , x 2 ( y ) C [ c , d ] .
Let f ( x , y ) be a scalar field, which is defined in the simple domain D x y R 0 , + 2 on the X Y -plane such that the following conditions are satisfied
f ( x , y 1 ( x ) ) , f ( x , y 2 ( x ) ) C 1 n ( 0 , ) ,
f ( x 1 ( y ) , x ) , f ( x 2 ( y ) , y ) C 1 n ( 0 , ) .
Then this property of the field f ( x , y ) will be denoted as f ( x , y ) F 1 n ( D x y ) .
Let us consider the vector field
F ( x , y ) = e x F x ( x , y ) + e y F y ( x , y )
which is defined in the simple domain D x y R 0 , + 2 . The fact that the field F ( x , y ) satisfies the conditions
F x ( x , y 1 ( x ) ) , F x ( x , y 2 ( x ) ) C 1 n ( 0 , ) ,
F y ( x 1 ( y ) , y ) , F y ( x ( y ) , y ) C 1 n ( 0 , ) ,
will be denoted as F ( x , y ) F 1 n ( D x y ) .
Let us define the properties of the scalar and vector fields on simple surfaces in R 0 , + 3 .
Definition 29.
Let S be a simple surface in R 0 , + 3 that is described by Equations (136)–(138), and D y z , D x z , D x y are simple domains, which are projections of the surface S onto the Y Z , X Z , X Y planes.
Let us consider the field f ( x , y , z ) that is defined on the simple surface S R 0 , + 3 . The fact that the field f ( x , y , z ) satisfies the conditions
f ( x , y , z ( x , y ) ) F 1 n ( D x y ) ,
f ( x ( y , z ) , y , z ) F 1 n ( D y z ) ,
f ( x , y ( x , z ) , z ) F 1 n ( D x z ) ,
will be denoted as f ( x , y , z ) F 1 n ( S ) .
Definition 30.
Let S be a simple surface in R 0 , + 3 that is described by Equations (136)–(138), and D y z , D x z , D x y are projections of the surface S onto the Y Z , X Z , X Y planes.
Let us consider the vector field
F ( x , y , z ) = e x F x ( x , y , z ) + e y F y ( x , y , z ) + e z F z ( x , y , z )
that is defined on the simple surface S R 0 , + 3 .
The fact that the field F ( x , y , z ) satisfies the conditions
F x ( x , y , z ( x , y ) ) , F y ( x , y , z ( x , y ) ) F 1 n ( D x y ) ,
F x ( x , y ( x , z ) , z ) , F z ( x , y ( x , z ) , z ) F 1 n ( D x z ) ,
F y ( x ( y , z ) , y , z ) , F z ( x ( y , z ) , y , z ) F 1 n ( D y z ) ,
will be denoted as F ( x , y , z ) F 1 n ( S ) .

8.4. General Fractional Vector Integrals over Surface

Let us consider a smooth oriented surface S in R 0 , + 3 , which is bounded by an oriented closed smooth line L = S . We will assume that S is the simple surface in R 0 , + 3 , and D y z , D x z , D x y are simple domains, which are projections of the surface S onto the Y Z , X Z , X Y planes.
Definition 31.
The surface GFI opereator over the simple surface S R 0 , + 3 is defined by the equation
I S ( M ) : = I D x y ( M ) + I D x z ( M ) + I D y z ( M ) =
n x I D y z ( M ) + n y I D x z ( M ) + n z I D x y ( M ) ,
where the vectors n x , n y , n z are the normal vector to surfaces D y z , D x z , D x y that are related with orientation of the closed lines L y z = D y z , L x z = D x z , L x y = D x y bounding its.
We can use the vectors e x , e y , e z instead of n x , n y , n z , if all orientation of the closed lines L y z = D y z , L x z = D x z , L x y = D x y is positive.
For the simple domain D x y , the vector operator I D x y ( M ) is defined by the equation
I D x y ( M ) = e z I D x y ( M ) = e z I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) = e z I [ c , d ] ( M 2 ) [ x ] I [ x 1 ( y ) , x 2 ( y ) ] ( M 1 ) ,
where 0 y 1 ( x ) y 2 ( x ) for x [ a , b ] , and 0 x 1 ( y ) x 2 ( y ) for y [ c , d ] . The minus sign in front of the last expression is due to the relationship between the normal vector of the surface D x y and the orientation of the closed contours L x y = D x y .
Similarly, we obtain expressions for the simple domains D x z , D y z .
Definition 32.
The surface GFI of the vector field F ( x , y , z ) F 1 0 ( S ) over the simple surface S R 0 , + 3 is defined by the equation
I S ( M ) , F = I D y z ( M ) [ y , z ] F x ( x ( y , z ) , y , z ) +
I D x z ( M ) [ x , z ] F y ( x , y ( x , z ) , z ) + I D x y ( M ) [ x , y ] F z ( x , y , z ( x , y ) ) ,
if n x = e x , n y = e y , n z = e z .
Definition 33.
Let S R 0 , + 3 be a Z-simple surface that surface can be represented as the union of X-simple surfaces S X , k and Y-simple surfaces S Y , k such that
S : = k = 1 n S X , k = l = 1 m S Y , l ,
and D y z , k , D x z , l are projections of S X , k and S Y , l on the Y Z and X Z planes.
The surface GFI of the vector field F ( x , y , z ) F 1 0 ( S ) over this surface S R 0 , + 3 is defined by the equation
I S ( M ) , F = k = 1 n I D y z , k ( M ) [ y , z ] F x ( x k ( y , z ) , y , z ) +
l = 1 l I D x z , l ( M ) [ x , z ] F y ( x , y l ( x , z ) , z ) + I D x y ( M ) [ x , y ] F z ( x , y , z ( x , y ) ) ,
if n x = e x , n y = e y , n z = e z .

8.5. General Fractional Curl Operators

8.5.1. Regional General Fractional Curl Operator

Let us give the definition of general fractional curl operator for R 0 , + 3 .
Definition 34.
Let F ( x , y , z ) be a vector field that belongs to the set C 1 1 ( R + 3 ) .
Then the regional general fractional curl (Regional GF Curl) for the region W = R 0 , + 3 is defined as
Curl W ( K ) F = e x Curl W ( K ) F x + e y Curl W ( K ) F y + e z Curl W ( K ) F z ,
where
Curl W ( K ) F x ( x , y , z ) = D ( K 2 ) y , [ y ] F z ( x , y , z ) D ( K 3 ) z , [ z ] F y ( x , y , z ) ,
Curl W ( K ) F y ( x , y , z ) = D ( K 3 ) z , [ z ] F x ( x , y , z ) D ( K 1 ) x , [ x ] F z ( x , y , z ) ,
Curl W ( K ) F z ( x , y , z ) = D ( K 1 ) x , [ x ] F y ( x , y , z ) D ( K 2 ) y , [ y ] F x ( x , y , z )
with x , y , z 0 .
Remark 20.
If the vector field
F ( x , y , z ) = k = 1 3 e k F k ( x , y , z )
belongs to the function space C 1 1 ( R 0 , + 3 ) , then the regional general fractional curl operator for the region W = R + 3 can be written in the compact form
( Curl W ( K ) F ) ( x , y , z ) = D W ( K ) , F = e n ε n m k D ( K n ) x n , [ x n ] F k ,
where
D W ( K ) : = e x D ( K 1 ) x , [ x ] + e y D ( K 2 ) y , [ y ] + e z D ( K 3 ) z , [ z ] ,
and ε l m n is Levi-Civita symbol, i.e., it is 1, if ( i , j , k ) is an even permutation of ( 1 , 2 , 3 ) , 0 if any index is repeated, and ( 1 ) if it is an odd permutation, and 0 if any index is repeated.
Remark 21.
The general fractional curl operator can be defined not only for W = R + 3 , but also for regions W R + 3 , surfaces S R + 3 and line L R + 3 . Note that the surface general fractional curl operator is used in the general fractional Stokes theorem to be given and proved in the following sections.
Definition 35.
Let F ( x , y , z ) belong to the space C 1 1 ( R + 3 ) .
Then the regional general fractional curl (the regional GF Curl) for the region
W : = { ( x , y , z ) : x a 0 , y c 0 , z e 0 }
is defined as
Curl W ( K ) F = e x Curl W ( K ) F x + e y Curl W ( K ) F y + e z Curl W ( K ) F z ,
where
Curl W ( K ) F x ( x , y , z ) = D [ c , y ] ( K 2 ) [ y ] F z ( x , y , z ) D [ e , z ] ( K 3 ) [ z ] F y ( x , y , z ) ,
Curl W ( K ) F y ( x , y , z ) = D [ e , z ] ( K 3 ) [ z ] F x ( x , y , z ) D [ a , x ] ( K 1 ) [ x ] F z ( x , y , z ) ,
Curl W ( K ) F y ( x , y , z ) = D [ a , x ] ( K 3 ) [ x ] F y ( x , y , z ) D [ c , y ] ( K 1 ) [ y ] F x ( x , y , z )
with x , y , z 0 .

8.5.2. Surface General Fractional Curl Operator

We have given the definition of the regional GF Curl operator. Let us now define the surface GF operator. In this definition we will use the notations
D z 1 F x ( x , y ( x , z ) , z ) : = F x ( x , y , z ) z y = y ( x , z ) , D y 1 F x ( x , y , z ( x , y ) ) : = F x ( x , y , z ) y z = z ( x , y ) ,
D x 1 F y ( x , y , z ( x , y ) ) : = F y ( x , y , z ) x z = z ( x , y ) , D z 1 F y ( x ( y , z ) , y , z ) : = F y ( x , y , z ) z z = z ( x , y ) ,
D x 1 F z ( x , y ( x , z ) , z ) : = F z ( x , y , z ) x y = y ( x , z ) , D y 1 F z ( x ( y , z ) , y , z ) : = F z ( x , y , z ) y x = x ( y , z ) .
Let us give a definition of the set F S 1 ( R 0 , + 3 ) of vector fields F that are used to define the surface GF Curl operator (compare with Definition 21 of the set F S 1 ( R 0 , + 3 ) that is used in definition of the surface GFI).
Definition 36.
Let S be a piecewise simple surface ( S P ( R 0 , + 3 ) ).
Let the vector field
F : = e x F x ( x , y , z ) + e y F y ( x , y , z ) + e z F z ( x , y , z ) .
on the surface S satisfy the conditions
D z 1 F x ( x , y j ( x , z ) , z ) C 1 ( R + 2 ) , D y 1 F x ( x , y , z k ( x , y ) ) C 1 ( R + 2 ) ,
D x 1 F y ( x , y , z ( x , y ) ) C 1 ( R + 2 ) , D z 1 F y ( x ( y , z ) , y , z ) C 1 ( R + 2 ) ,
D x 1 F z ( x , y ( x , z ) , z ) C 1 ( R + 2 ) , D y 1 F z ( x ( y , z ) , y , z ) C 1 ( R + 2 ) .
for all i = 1 , , n x , j = 1 , , n y , k = 1 , , n z .
The set of such vector fields F on piecewise simple surface S will be denoted by F S 1 ( R 0 , + 3 ) .
Then the surface general fractional curl is defined in the following form. The piecewise simple surface is defined in Definition 20.
Definition 37.
Let S be a piecewise simple surface ( S P ( R 0 , + 3 ) ) and a vector fields F on this surface S belongs to the set F S 1 ( R 0 , + 3 ) .
Then the surface general fractional curl (surface GF Curl) of the vector field F ( x , y , z ) for the piecewise simple surface S is defined as
Curl S ( K ) F = e x Curl S ( K ) F x + e y Curl S ( K ) F y + e z Curl S ( K ) F z ,
where
Curl S ( K ) F x = i = 1 n x I ( K 2 ) y , [ y ] D y 1 F z ( x i ( y , z ) , y , z ) I ( K 3 ) z , [ z ] D z 1 F y ( x i ( y , z ) , y , z ) ,
Curl S ( K ) F y = j = 1 n y I ( K 3 ) z , [ z ] D z 1 F x ( x , y j ( x , z ) , z ) I ( K 1 ) x , [ x ] D x 1 F z ( x , y j ( x , z ) , z ) ,
Curl S ( K ) F z = k = 1 n z I ( K 1 ) x , [ x ] D x 1 F y ( x , y , z k ( x , y ) ) I ( K 2 ) y , [ y ] D y 1 F x ( x , y , z k ( x , y ) )
with x , y , z 0 .

8.6. Difficulties in Generalization of Stokes Theorem

In the standard Stokes theorem, the following equality should be satisfied
S P ( x , y , z ) d x = S P ( x , y , z ) z d x d z P ( x , y , z ) y d x d y .
Equation (140) can be regarded as the Stokes formula for the vector field F = e x P ( x , y . z ) .
If S R 0 , + 3 is Z-simple and Y-simple surface that is described by the equation z = z ( x , y ) for ( x , y ) D x y R 0 , + 2 , and the equation y = y ( x , z ) for ( x , z ) D x z R 0 , + 2 , then
S Q z ( x , y , z ) d x d y = D x y Q ( x , y , z ( x , y ) ) d x d y ,
S Q y ( x , y , z ) d x d z = D x z Q y ( x , y ( x , z ) , z ) d x d z ,
where D x y and D x z are the projections of the surface S on the X Y and X Z -planes.
The values of the function P ( x , y , z ) on the line L R 0 , + 3 are equal to the values of the function P ( x , y , z ( x , y ) ) on the line L x y R 0 , + 2 , which is the projection of the line L onto the X Y -plane
S P ( x , y , z ) d x = D x y P ( x , y , z ( x , y ) ) d x ,
where L = S and L x y = D x y .
Therefore, equality (140) means that
D x y P ( x , y , z ( x , y ) ) d x = D x z P ( x , y , z ( x , y ) ) y d x d z = D x z P ( x , y , z ) y z = z ( x , y ) + P ( x , y , z ) z z = z ( x , y ) z ( x , y ) y d x d z = D x z P ( x , y , z ) z y = y ( x , z ) d x d z D x y P ( x , y , z ) y z = z ( x , y ) d x d y .
Equality (140) is based on the standard chain rule
Q ( x , y , z ( x , y ) ) y = Q ( x , y , z ) y z = z ( x , y ) + Q ( x , y , z ) z z = z ( x , y ) z ( x , y ) y .
For the fractional derivatives and GFD, the standard chain is violated [3], pp. 97–98.
For the general fractional Stokes theorem, we should have the equality
I D x y ( M 1 ) [ x ] P ( x , y , z ) = I D x z ( M ) [ x , z ] D K 3 z , [ z ] P ( x , y , z ) y = y ( x , z ) I D x y ( M ) [ x , y ] D K 2 y , [ y ] P ( x , y , z ) z = z ( x , y ) ,
which is a fractional analog of Equation (141), where S is the Z-simple and Y-simple surface.
Identity (142) does not hold in the general case, due to the inequalities
I [ c , d ] ( M 2 ) [ y ] D y , [ y ] P ( x , y , z ) z = z ( x , y ) I [ c , d ] ( M 2 ) [ y ] D y , [ y ] P ( x , y , z ) z = z ( x , y ) .
and
I [ c , d ] ( M 2 ) [ y ] D y , [ y ] P ( x , y , z ) z = z ( x , y ) I [ c , d ] ( M 2 ) [ y ] D y , [ y ] P ( x , y , z ) z = z ( x , y ) .
We can state that we have the equality only if Z-simple surface S is described by equation z = z ( x ) for all ( x , y ) D x , y . Then the equalities have the form
I [ c , d ] ( M 2 ) [ y ] D y , [ y ] P ( x , y , z ) z = z ( x ) = I [ c , d ] ( M 2 ) [ y ] D y , [ y ] P ( x , y , z ) z = z ( x ) =
P ( x , d , z ( x ) ) P ( x , c , z ( x ) ) .
and
I [ c , d ] ( M 2 ) [ y ] D y , [ y ] P ( x , y , z ) z = z ( x ) =
I [ c , d ] ( M 2 ) [ y ] D y , [ y ] P ( x , y , z ) z = z ( x ) = P ( x , d , z ( x ) ) P ( x , c , z ( x ) ) .
A similar situation for Y-simple surface that is described by the equation y = y ( x , z ) with x , z D x z .
As a result, for the regional GF curl Curl W ( K ) with W = R 0 , + 3 , the general fractional theorem can be proved only if the surface consists of faces parallel to the X Y , X Z , Y Z planes.

8.7. General Fractional Stokes Theorem for Box without Bottom

Let us consider a box-shaped surface S without a bottom, i.e., a parallelepiped surface without a bottom face (base). The parallelepiped can be described as
W : = { ( x , y , z ) : 0 a x b , 0 c y d , 0 e z f } .
The vertices of the parallelepiped have the coordinates
A ( a , c , e ) , A ( a , c , f ) ,
B ( b , c , e ) , B ( b , c , f ) ,
C ( b , d , e ) , C ( b , d , f ) ,
D ( a , d , e ) , D ( a , d , f ) .
The boundary of this parallelepiped is
W = S A B C D ,
where the surface A B C D is a bottom face of the parallelepiped region, the surface S consists of the following faces
S = D x y D x z D y z ,
where
D x y = A B C D ,
D x z = A B B A C D D C ,
D y z = A D D A B C C B
without the face A B C D .
The closed line L = S , which is the boundary of the surface S in the form of rectangle A B C D , consists of the segments
L = A B C D = A B B C C D D A .
The surface GFI operator is desribed as
I S ( M ) : = e x I D y z ( M ) [ y , z ] + e y I D x z ( M ) [ x , z ] + e z I D x y ( M ) [ x , y ] .
Then
I S ( M ) , F = I D y z ( M ) [ y , z ] F x ( x , y , z ) + I D x z ( M ) [ x , z ] F y ( x , y , z ) + I D x y ( M ) [ x , y ] F z ( x , y , z ) .
To calculate the expression I S ( M ) , Curl S ( K ) F , we should consider the following surface GFIs
I S ( M ) , Curl S ( K ) F =
I D y z ( M ) [ y , z ] D ( K 2 ) y , [ y ] F z D ( K 3 ) z , [ z ] F y +
I D x z ( M ) [ x , z ] D ( K 3 ) z , [ z ] F x D ( K 1 ) x , [ x ] F z +
I D x y ( M ) [ x , y ] D ( K 1 ) x , [ x ] F y D ( K 2 ) y , [ y ] F x .
Theorem 15.
(General fractional Stokes theorem for parallelepiped surface without bottom) Let S is a smooth oriented surface (144) in R 0 , + 3 , which is bounded by an oriented closed smooth line L = S that is given by (148), and D y z , D x z , D x y are domains (145)–(147), which are projections of the surface S onto the Y Z , X Z , X Y planes.
Then, for the vector field F ( x , y , z ) F 1 1 ( S ) , we have the equation
I S ( M ) , F = I S ( M ) , Curl S ( K ) F ,
which is the general fractional Stokes equation.
For the surface S in the form of the parallelepiped without a bottom face, the general fractional Stokes equation has the form
I D y z ( M ) [ y , z ] D ( K 2 ) y , [ y ] F z D ( K 3 ) z , [ z ] F y +
I D x z ( M ) [ x , z ] D ( K 3 ) z , [ z ] F x D ( K 1 ) x , [ x ] F z +
I D x y ( M ) [ x , y ] D ( K 1 ) x , [ x ] F y D ( K 2 ) y , [ y ] F x =
I A B ( M 1 ) [ x ] F x I D C ( M 1 ) [ x ] F x + I B C ( M 2 ) [ y ] F y I A D ( M 2 ) [ y ] F y .
Proof. 
(YZ) Let us consider the GFI operator I D y z ( M ) . The surface D y z consists of two X-simple surfaces B C C B and A D D A . Therefore, we get
I D y z ( M ) [ y , z ] = I B C C B ( M ) [ y , z ] I A D D A ( M ) [ y , z ] ,
where we took into account the direction of the normals to the outer surface with respect to the direction of the basis vector e x .
(YZ1) For the first X-simple surface B C C B .
I B C C B ( M ) [ y , z ] D ( K 2 ) y , [ y ] F z D ( K 3 ) z , [ z ] F y =
I B C C B ( M ) [ y , z ] D ( K 2 ) y , [ y ] F z ( b , y , z ) I B C C B ( M ) [ y , z ] D ( K 3 ) z , [ z ] F y ( b , y , z ) =
I [ e , f ] ( M 3 ) [ z ] I [ c , d ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] F z ( b , y , z ) I [ c , d ] ( M 2 ) [ y ] I [ e , f ] ( M 3 ) [ z ] D ( K 3 ) z , [ z ] F y ( b , y , z ) =
I [ e , f ] ( M 3 ) [ z ] F z ( b , d , z ) F z ( b , c , z ) I [ c , d ] ( M 2 ) [ y ] F y ( b , y , f ) F y ( b , y , e ) =
I [ e , f ] ( M 3 ) [ z ] F z ( b , d , z ) I [ e , f ] ( M 3 ) [ z ] F z ( b , c , z ) I [ c , d ] ( M 2 ) [ y ] F y ( b , y , f ) + I [ c , d ] ( M 2 ) [ y ] F y ( b , y , e ) =
I C C ( M 3 ) [ z ] F z ( x , y , z ) I B B ( M 3 ) [ z ] F z ( x , y , z ) I B C ( M 2 ) [ y ] F y ( x , y , z ) + I B C ( M 2 ) [ y ] F y ( x , y , z ) .
(YZ2) For the second X-simple surface A D D A .
I A D D A ( M ) [ y , z ] D ( K 2 ) y , [ y ] F z D ( K 3 ) z , [ z ] F y =
I A D D A ( M ) [ y , z ] D ( K 2 ) y , [ y ] F z ( a , y , z ) + I A D D A ( M ) [ y , z ] D ( K 3 ) z , [ z ] F y ( a , y , z ) =
I [ e , f ] ( M 3 ) [ z ] I [ c , d ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] F z ( a , y , z ) + I [ c , d ] ( M 2 ) [ y ] I [ e , f ] ( M 3 ) [ z ] D ( K 3 ) z , [ z ] F y ( a , y , z ) =
I [ e , f ] ( M 3 ) [ z ] F z ( a , d , z ) F z ( a , c , z ) + I [ c , d ] ( M 2 ) [ y ] F y ( a , y , f ) F y ( a , y , e ) =
I [ e , f ] ( M 3 ) [ z ] F z ( a , d , z ) + I [ e , f ] ( M 3 ) [ z ] F z ( a , c , z ) + I [ c , d ] ( M 2 ) [ y ] F y ( a , y , f ) I [ c , d ] ( M 2 ) [ y ] F y ( a , y , e ) =
I D D ( M 3 ) [ z ] F z ( x , y , z ) + I A A ( M 3 ) [ z ] F z ( x , y , z ) + I A D ( M 2 ) [ y ] F y ( x , y , z ) I A D ( M 2 ) [ y ] F y ( x , y , z ) .
(XZ) Let us consider the GFI operator I D x z ( M ) [ x , z ] . The surface D x z consists of two Y-simple surfaces A B B A and C D D C . Therefore, we get
I D x z ( M ) [ x , z ] = I A B B A ( M ) [ x , z ] + I C D D C ( M ) [ x , z ] .
Here we took into account the direction of the normals to the outer surface with respect to the direction of the basis vector e y . The minus in front of the operator I A B B A ( M ) is due to the fact that the normal to the surface is directed in the opposite direction with respect to the vector e y .
(XZ1) For the first Y-simple surface A B B A .
I A B B A ( M ) [ x , z ] D ( K 3 ) z , [ z ] F x D ( K 1 ) x , [ x ] F z =
I A B B A ( M ) [ x , z ] D ( K 3 ) z , [ z ] F x + I A B B A ( M ) [ x , z ] D ( K 1 ) x , [ x ] F z ) =
I [ a , b ] ( M 1 ) [ x ] I [ e , f ] ( M 3 ) D ( K 3 ) z , [ z ] F x ( x , c , z ) + I [ e , f ] ( M 3 ) I [ a , b ] ( M 1 ) [ x ] D ( K 1 ) x , [ x ] F z ( x , c , z ) =
I [ a , b ] ( M 1 ) [ x ] F x ( x , c , f ) F x ( x , c , e ) + I [ e , f ] ( M 3 ) ( F z ( b , c , z ) F z ( a , c , z ) ) =
I A B ( M 1 ) [ x ] F x ( x , y , z ) + I A B ( M 1 ) [ x ] F x ( x , y , z ) + I B B ( M 3 ) [ z ] F z ( x , y , z ) I A A ( M 3 ) [ z ] F z ( x , y , z ) .
(XZ2) For the second Y-simple surface C D D C .
I C D D C ( M ) [ x , z ] D ( K 3 ) z , [ z ] F x D ( K 1 ) x , [ x ] F z =
I C D D C ( M ) [ x , z ] D ( K 3 ) z , [ z ] F x I C D D C ( M ) [ x , z ] D ( K 1 ) x , [ x ] F z ) =
I [ a , b ] ( M 1 ) [ x ] I [ e , f ] ( M 3 ) [ z ] D ( K 3 ) z , [ z ] F x ( x , d , z ) I [ e , f ] ( M 3 ) [ z ] I [ a , b ] ( M 1 ) [ x ] D ( K 1 ) x , [ x ] F z ( x , d , z ) ) =
I [ a , b ] ( M 1 ) [ x ] F x ( x , d , f ) F x ( x , d , e ) I [ e , f ] ( M 3 ) [ z ] F z ( b , d , z ) F z ( a , d , z ) =
I D C ( M 1 ) [ x ] F x ( x , y , z ) I D C ( M 1 ) [ x ] F x ( x , y , z ) I C C ( M 3 ) [ z ] F z ( x , y , z ) + I D D ( M 3 ) [ z ] F z ( x , y , z ) .
(XY) Let us consider the GFI operator I D x y ( M ) [ x , y ] . The surface D x y consists of one Z-simple surface A B C D .
I D x y ( M ) [ x , y ] = I A B C D ( M ) [ x , y ] .
Therefore, we get
I A B C D ( M ) [ x , y ] D ( K 1 ) x , [ x ] F y D ( K 2 ) y , [ y ] F x =
I [ a , b ] ( M 1 ) [ x ] I [ c , d ] ( M 2 ) [ y ] D ( K 1 ) x , [ x ] F y ( x , y , f ) D ( K 2 ) y , [ y ] F x ( x , y , f ) =
I [ c , d ] ( M 2 ) [ y ] I [ a , b ] ( M 1 ) [ x ] D ( K 1 ) x , [ x ] F y ( x , y , f ) I [ a , b ] ( M 1 ) [ x ] I [ c , d ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] F x ( x , y , f ) ) =
I [ c , d ] ( M 2 ) [ y ] F y ( b , y , f ) F y ( a , y , f ) I [ a , b ] ( M 1 ) [ x ] F x ( x , d , f ) F x ( x , c , f ) =
I B C ( M 2 ) [ y ] F y ( x , y , z ) I A D ( M 2 ) [ y ] F y ( x , y , z ) I D C ( M 1 ) [ x ] F x ( x , y , z ) + I A B ( M 1 ) [ x ] F x ( x , y , z ) .
(XYZ) As a result, we get
I S ( M ) , Curl S ( K ) F =
I C C ( M 3 ) [ z ] F z I B B ( M 3 ) [ z ] F z I B C ( M 2 ) [ y ] F y + I B C ( M 2 ) [ y ] F y +
I D D ( M 3 ) [ z ] F z + I A A ( M 3 ) [ z ] F z + I A D ( M 2 ) [ y ] F y I A D ( M 2 ) [ y ] F y +
I A B ( M 1 ) [ x ] F x + I A B ( M 1 ) [ x ] F x + I B B ( M 3 ) [ z ] F z I A A ( M 3 ) [ z ] F z +
I D C ( M 1 ) [ x ] F x I D C ( M 1 ) [ x ] F x I C C ( M 3 ) [ z ] F z + I D D ( M 3 ) [ z ] F z +
I B C ( M 2 ) [ y ] F y I A D ( M 2 ) [ y ] F y I D C ( M 1 ) [ x ] F x + I A B ( M 1 ) [ x ] F x .
Bringing down similar terms, we obtain
I S ( M ) , Curl S ( K ) F =
I A B ( M 1 ) [ x ] F x I D C ( M 1 ) [ x ] F x + I B C ( M 2 ) [ y ] F y I A D ( M 2 ) [ y ] F y =
I L x ( M 1 ) [ x ] F x + I L y ( M 2 ) [ y ] F y = I L ( M ) , F ,
where L = A B C D is the closed line in the X Y -plane, which is the boundary L = S of the surface S. □

8.8. General Fractional Stokes Theorem for Surface GF Curl

Let us prove the general fractional Stokes theorem for surface general fractional gradient, where surface consist of simple surfaces or surfaces parallel to the coordinate planes.
Theorem 16.
(General Fractional Stokes Theorem for Surface GF Curl)
Let S R 0 , + 3 be a simple surface (or a piecewise simple surface S P ( R 0 , + 3 ) ), and the vector field F belongs to the set F S 1 ( R 0 , + 3 ) .
Then, the equation
I S ( M 1 ) , F = I S ( M ) , Curl S ( K ) F
holds, where Curl S ( K ) is the surface GF Curl.
Proof. 
Let us prove the general Stokes theorem for the vector field F x = e x P ( x , y , z ) . The general Stokes equations for the vector fields F y = e y F y ( x , y , z ) and F z = e z F z ( x , y , z ) are proved similarly.
Let S R 0 , + 3 be a Z-simple and Y-simple surface that is described by the equation z = z ( x , y ) for ( x , y ) D x y R 0 , + 2 , and the equation y = y ( x , z ) for ( x , z ) D x z R 0 , + 2 .
The values of the function P ( x , y , z ) on the line L R 0 , + 3 are equal to the values of the function P ( x , y , z ( x , y ) ) on the line L x y R 0 , + 2 , which is a projection of the line L onto the X Y -plane,
I S ( M 1 ) [ x ] P ( x , y , z ) = I D x y ( M 1 ) [ x ] P ( x , y , z ( x , y ) ) ,
where L = S and L x y = D x y .
Let us assume that L x y consists of two Y-simple lines L x y 1 and L x y 2 that are described by equations y = y 1 ( x ) 0 and y = y 2 ( x ) , where y 2 ( x ) y 1 ( x ) 0 . Then
I L x y ( M 1 ) [ x ] P ( x , y , z ( x , y ) ) = I L x y 1 ( M 1 ) [ x ] P ( x , y , z ( x , y ) ) + I L x y 2 ( M 1 ) [ x ] P ( x , y , z ( x , y ) ) = I [ a , b ] ( M 1 ) [ x ] P ( x , y 1 ( x ) , z ( x , y 1 ( x ) ) ) I [ a , b ] ( M 1 ) [ x ] P ( x , y 2 ( x ) , z ( x , y 2 ( x ) ) ) = I [ a , b ] ( M 1 ) [ x ] P ( x , y 1 ( x ) , z ( x , y 1 ( x ) ) ) P ( x , y 2 ( x ) , z ( x , y 2 ( x ) ) ) .
Using the fundamental theorem of general fractional calculus, expression (149) can be written as
I [ a , b ] ( M 1 ) [ x ] P ( x , y 1 ( x ) , z ( x , y 1 ( x ) ) ) P ( x , y 2 ( x ) , z ( x , y 2 ( x ) ) ) = I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] P ( x , y , z ( x , y ) ) .
Then using the definition of the general fractional derivative and the standard chain rule for the first-order derivative, we get
I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] D ( K 2 ) y , [ y ] P ( x , y , z ( x , y ) ) =
I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] P ( x , y , z ( x , y ) ) y =
I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] P ( x , y , z ) y z = z ( x , y ) + P ( x , y , z ) z z = z ( x , y ) z ( x , y ) y =
I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] P ( x , y , z ) y z = z ( x , y )
I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] P ( x , y , z ) z z = z ( x , y ) z ( x , y ) y .
Let us use the property
I [ y 1 , y 2 ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] U ( y ) = ( M 2 ( K 2 U ) ) ( y 2 ) ( M 2 ( K 2 U ) ) ( y 1 ) =
( ( M 2 K 2 ) U ) ( y 2 ) ( ( M 2 K 2 ) U ) ( y 1 ) =
( { 1 } U ) ( y 2 ) ( { 1 } U ) ( y 1 ) = y 1 y 2 U ( y ) d y
for the expression of term (152) in the form
I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] P ( x , y , z ) z z = z ( x , y ) z ( x , y ) y = y 1 ( x ) y 2 ( x ) P ( x , y , z ) z z = z ( x , y ) z ( x , y ) y d y .
Then we can use the equation, which is used in the standard Stokes theorem, ( { 1 } U ) ( z ) = ( M 3 K 3 U ) ( z ) , and the fundamental theorem of GFC in the form
y 1 ( x ) y 2 ( x ) P ( x , y , z ) z z = z ( x , y ) z ( x , y ) y d y = z 1 ( x ) z 2 ( x ) P ( x , y , z ) z y = y ( x , z ) d z = I [ z 1 ( x ) , z 2 ( x ) ] ( M 3 ) [ z ] I ( K 3 ) z , [ z ] P ( x , y , z ) z y = y ( x , z ) .
Using (154), expression (152) can be written in the form
I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] P ( x , y , z ) z z = z ( x , y ) z ( x , y ) y = I [ a , b ] ( M 1 ) [ x ] I [ z 1 ( x ) , z 2 ( x ) ] ( M 3 ) [ z ] I ( K 3 ) z , [ z ] P ( x , y , z ) z y = y ( x , z ) .
Therefore, for two terms (151) and (152), we obtain the equations
I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] I ( K 2 ) y , [ y ] P ( x , y , z ) y z = z ( x , y ) = I D x y ( M ) [ x , y ] Curl S ( K ) F x z ,
I [ a , b ] ( M 1 ) [ x ] I [ z 1 ( x ) , z 2 ( x ) ] ( M 3 ) [ z ] I ( K 3 ) z , [ z ] P ( x , y , z ) z y = y ( x , z ) = I D x z ( M ) [ x , z ] Curl S ( K ) F x y ,
for the vector field F x = e x P ( x , y , z ) .
As a result, we proved the general Stokes theorem for the surface FG Curl and the vector field F x = e x P ( x , y , z ) .
I S ( M 1 ) , F x = I S ( M ) , Curl S ( K ) F x .
Equations for the remaining components of the vector field F y = e y F y ( x , y , z ) and F z = e z F z ( x , y , z ) are proved similarly. □

9. General Fractional Gauss Theorem

9.1. Definition of Triple GFI by Iterated GFI

Let us define concept of the Z-simple region W in R 0 , + 3 .
Definition 38.
Let W be region in R 0 , + 3 that is bounded above and below by smooth surfaces S 2 , x y , S 1 , x y and a lateral surface S z , whose generatrices are parallel to the Z-axis. Let surfaces S 1 , x y , S 2 , x y be described by the equations
z = z 1 ( x , y ) 0 , z = z 2 ( x , y ) 0 ,
where the functions are continuous in the closed domain W x y that is a projection of the region W onto the X Y -plane, and z 2 ( x , y ) z 1 ( x , y ) for all ( x , y ) W x y .
Then, the region W will be called the Z-simple region (simple area along the Z-axis). The region W is called simple, if W is simple along three axes (X, Y, Z). If W can be divided into a finite number of such regions with respect to all three axes, then W will be called piecewise simple region in R 0 , + 3 .
Definition 39.
Let W R 0 , + 3 be Z-simple domain that is that is bounded above and below by smooth surfaces S 2 , x y , S 1 , x y described by Equation (158).
Let scalar field f ( x , y , z ) be satisfy the condition
J 2 ( x , y ) = I [ z 1 ( x , y ) , z 2 ( x , y ) ] ( M 3 ) [ z ] f ( x , y , z ) C 1 ( R + 2 ) .
Then the triple general fractional integral (triple GFI) is defined by the equation
I W ( M ) [ x , y , z ] f ( x , y , z ) : = I W x y ( M ) [ x , y ] I [ z 1 ( x , y ) , z 2 ( x , y ) ] ( M 3 ) [ z ] f ( x , y , z ) ,
where I W x y ( M ) [ x , y ] is double GFI.
Definition 40.
Let W R 0 , + 3 be Z-simple domain that is that is bounded above and below by smooth surfaces S 2 , x y , S 1 , x y described by Equation (158). Let W x y R 0 , + 2 be projection of W on the X Y -plane such that W x y is Y-simple region in X Y -plane that is bounded by the lines y = y 2 ( x ) and y = y 1 ( x ) , where y = y 2 ( x ) and y = y 1 ( x ) are continuous functions on the interval [ a , b ] , b > a 0 and y 2 ( x ) y 1 ( x ) > 0 for all x [ a , b ] .
Let scalar fiels f ( x , y , z ) satisfy condition (159) and
J 1 ( x ) = I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] I [ z 1 ( x , y ) , z 2 ( x , y ) ] ( M 3 ) [ z ] f ( x , y , z ) C 1 ( 0 , ) .
Then the triple general fractional integral (triple GFI) is defined in the form
I W ( M ) [ x , y , z ] f ( x , y , z ) : = I [ a , b ] ( M 1 ) [ x ] I [ y 1 ( x ) , y 2 ( x ) ] ( M 2 ) [ y ] I [ z 1 ( x , y ) , z 2 ( x , y ) ] ( M 3 ) [ z ] f ( x , y , z ) .
A volume general fractional integral (volume GFI) of a scalar field is a triple general fractional integral for the Z-simple region W R 0 , + 3 .
Example 18.
Using the parallelepiped
W : = { 0 x a , 0 y b , 0 z c } ,
the volume general fractional integral can be written as
I W ( M ) [ x , y , z ] f ( x , y , z ) = 0 a d x 0 b d y 0 c d z M 1 ( a x ) M 2 ( b y ) M 3 ( c z ) f ( x , y , z ) .
For the kernels (62), the general fractional flux (161) has the form
I W ( M ) [ x , y , z ] f ( x , y , z ) = 0 a d x 0 b d y 0 c d z ( a x ) α 1 1 Γ ( α 1 ) ( b y ) α 2 1 Γ ( α 2 ) ( c z ) α 3 1 Γ ( α 3 ) f ( x , y , z ) .
For α 1 = α 2 = α 3 = 1 , Equation (162) is given as
I W ( M ) [ x , y , z ] f ( x , y , z ) = W d V f ( x , y , z ) = 0 a d x 0 b d y 0 c d z f ( x , y , z ) ,
which is the standard volume integral for the function f ( x , y , z ) .

9.2. General Fractional Divergence

In this section, we give the definition of general fractional divergence for R 0 , + 3 .
Let us define sets of vector fields that will be used in the definition of the regional general fractional divergence.
Definition 41.
Let F ( x , y , z ) be a vector field that satisfies the conditions
D ( K 1 ) x , [ x ] F x ( x , y , z ) C 1 ( R + 3 ) , D ( K 2 ) y , [ y ] F y ( x , y , z ) C 1 ( R + 3 ) , D ( K 3 ) z , [ z ] F z ( x , y , z ) C 1 ( R + 3 ) .
Then the set of such vector fields will be denoted as F 1 , D i v 1 ( R + 3 ) .
We can also consider vector fields F ( x , y , z ) that belong to the function space C 1 1 ( R + 3 ) .
Definition 42.
Let F ( x , y , z ) be a vector field that satisfies the conditions
F x ( x , y , z ) , F y ( x , y , z ) , F z ( x , y , z ) C 1 1 ( R + 3 ) .
Then the set of such vector fields will be denoted as C 1 1 ( R + 3 ) .
In other words, the condition F ( x , y , z ) C 1 1 ( R + 3 ) means that all general fractional derivatives of all components of the vector field F ( x , y , z ) with respect to all coordinated belong to space C 1 ( R + 3 ) .
Let us define the general fractional divergence.
Definition 43.
Let F ( x , y , z ) be a vector field that belongs to the set F 1 , D i v 1 ( R + 3 ) or F ( x , y , z ) C 1 1 ( R + 3 ) .
Then the general fractional divergence Div W ( K ) for the region W = R 0 , + 3 is defined as
Div W ( K ) F = D W ( K ) , F =
D ( K 1 ) x , [ x ] F x ( x , y , z ) + D ( K 2 ) y , [ y ] F y ( x , y , z ) + D ( K 3 ) z , [ z ] F z ( x , y , z ) .
Remark 22.
The formula defining the operator can be written in compact form. If the vector field
F ( x , y , z ) = k = 1 3 e k F k ( x , y , z )
belongs to the function space C 1 1 ( R + 3 ) , then the general fractional divergence for the region W = R + 3 is defined as
( Div W ( K ) F ) ( x , y , z ) = D W ( K ) , F = k = 1 3 D ( K k ) x k , [ x k ] F k ,
where
D W ( K ) : = e x D ( K 1 ) x , [ x ] + e y D ( K 2 ) y , [ y ] + e z D ( K 3 ) z , [ z ] .
Remark 23.
The general fractional divergence can be defined not only for W = R 0 , + 3 , but also for regions W R 0 , + 3 , surfaces S R 0 , + 3 and line L R 0 , + 3 . Note that the general fractional divergence for regions W R 0 , + 3 is used in the general fractional Gausss theorem to be given and proved in the following sections.
Definition 44.
Let F ( x , y , z ) be a vector field that belongs to the set C 1 1 ( R + 3 ) , and the region W be defined in the form
W : = { ( x , y , z ) : x a 0 , y c 0 , z e 0 } .
Then the general fractional divergence Div W ( K ) for the region W is defined as
Div W ( K ) F = D W ( K ) , F =
D [ a , x ] ( K 1 ) [ x ] F x ( x , y , z ) + D [ c , y ] ( K 2 ) [ y ] F y ( x , y , z ) + D [ e , z ] ( K 3 ) [ z ] F z ( x , y , z ) .

9.3. General Fractional Gauss Theorem for Z-Simple Region

The standard Gauss theorem (the Gauss-Ostrogradsky theorem) relates the flux of a vector field through a closed surface to the divergence of the field in the region enclosed. The Gauss theorem states that the surface integral of a vector field over a closed surface, which is the flux through the surface, is equal to the volume integral of the divergence over the region inside the surface.
Let us define a set of vector field, for which general fractional Gauss theorem will be formulated.
Definition 45.
Let W in R 0 , + 3 be Z-simple region such that W is a piecewise Y-simple and X-simple region
W = k = 1 n W X , k = j = 1 m W Y , j ,
where W X , k are the X-simple regions that is described by x = x k , 1 ( y , z ) , x = x k , 2 ( y , z ) for y , z D y z , and W Y , j are the Y-simple regions that are described by the functions y = y j , 1 ( x , z ) , y = y j , 2 ( x , z ) for ( x , z ) D x z .
Let vector field F ( x , y , z ) satisfy the conditions
F x ( x k , 2 ( y , z ) , y , z ) , F x ( x k , 1 ( y , z ) , y , z ) C 1 ( R + 2 ) ,
F y ( x , y j , 2 ( x , z ) , z ) , F y ( x , y j , 1 ( x , z ) , z ) C 1 ( R + 2 ) ,
F z ( x , y , z 1 ( x , y ) ) , F z ( x , y , z 2 ( x , y ) ) C 1 ( R + 2 ) .
Then the set of such vector fields F will be denoted as F ( x , y , z ) F 1 ( W ) .
f x ( x , y , z ) : = D ( K 1 ) x , [ x ] F x ( x , y , z ) C 1 ( R + 3 ) ,
f y ( x , y , z ) : = D ( K 2 ) y , [ y ] F y ( x , y , z ) C 1 ( R + 3 ) ,
f z ( x , y , z ) : = D ( K 3 ) z , [ z ] F z ( x , y , z ) C 1 ( R + 3 ) .
Then the set of such vector fields F will be denoted as F ( x , y , z ) F 1 1 ( R + 3 ) .
Theorem 17.
(General fractional Gauss theorem for Z-simple region)
Let W in R 0 , + 3 be Z-simple region such that W is a piecewise Y-simple and X-simple region. Let W is bounded above and below by smooth surfaces S 2 , x y , S 1 , x y , which is described by Equation (158), and a lateral surface S z , whose generatrices are parallel to the Z-axis.
Let the vector field F ( x , y , z ) belongs to the sets F 1 ( W ) and F 1 1 ( R + 3 ) .
Then the general fractional Gauss equation has the from
I W ( M ) [ x , y , z ] Div W ( K ) F = I W , F
that can be written as
I W ( M ) [ x , y , z ] D ( K 1 ) x , [ x ] F x ( x , y , z ) + D ( K 2 ) y , [ y ] F y ( x , y , z ) + D ( K 3 ) z , [ z ] F z ( x , y , z ) =
I S ( M ) [ y , z ] F x ( x , y , z ) + I S ( M ) [ x , y ] F y ( x , y , z ) + I S ( M ) [ x , y ] F z ( x , y , z ) .
Proof. 
Let us consider the triple GFI
I W ( M ) [ x , y , z ] = I W x y ( M ) [ x , y ] I [ z 1 ( x , y ) , z 2 ( x , y ) ] ( M 3 ) [ z ]
for the function
f z ( x , y , z ) = D ( K 3 ) z , [ z ] F z ( x , y , z ) C 1 ( 0 , ) ,
where F z ( x , y , z ) C 1 1 ( 0 , ) for each point ( x , y ) D x , y .
Using the second fundamental theorem of GFC, we get
I [ z 1 ( x , y ) , z 2 ( x , y ) ] ( M 3 ) [ z ] D ( K 3 ) z , [ z ] F z ( x , y , z ) =
I ( M 3 ) z 2 ( x , y ) [ z ] D ( K 3 ) z , [ z ] F z ( x , y , z ) I ( M 3 ) z 1 ( x , y ) [ z ] D ( K 3 ) z , [ z ] F z ( x , y , z ) =
F z ( x , y , z 2 ( x , y ) ) F z ( x , y , 0 ) F z ( x , y , z 2 ( x , y ) ) F z ( x , y , 0 ) ) =
F z ( x , y , z 2 ( x , y ) ) F z ( x , y , z 1 ( x , y ) ) .
If we assume that
F z ( x , y , z 1 ( x , y ) ) , F z ( x , y , z 2 ( x , y ) ) C 1 ( R + 2 ) ,
then we get
I W ( M ) [ x , y , z ] D ( K 3 ) z , [ z ] F z ( x , y , z ) = I W x y ( M ) [ x , y ] F z ( x , y , z 2 ( x , y ) ) F z ( x , y , z 1 ( x , y ) ) .
Therefore, the triple GFI can be represented through the surface GFIs in the form
I W ( M ) [ x , y , z ] D ( K 3 ) z , [ z ] F z ( x , y , z ) = I S 2 , x y ( M ) [ x , y ] F z ( x , y , z ) + I S 1 , x y ( M ) [ x , y ] F z ( x , y , z ) ,
where the surface GFIs are represented by the definition in the form
I S 2 , x y ( M ) [ x , y ] F z ( x , y , z ) = I W x y ( M ) [ x , y ] F z ( x , y , z 2 ( x , y ) ) ,
I S 1 , x y ( M ) [ x , y ] F z ( x , y , z ) = I W x y ( M ) [ x , y ] F z ( x , y , z 1 ( x , y ) ) .
Then, taking into account that the surface GFI oves the surface S z is equati to zero
I S z ( M ) [ x , y ] F z ( x , y , z ) = 0 ,
we obtain
I W ( M ) [ x , y , z ] D ( K 3 ) z , [ z ] F z ( x , y , z ) = I S ( M ) [ x , y ] F z ( x , y , z ) ,
where W = S = S 1 , x y S 2 , x y S z is closed surface that contains the region W inside the surface W .
If W is piecewise Y-simple and X-simple region, such that
W = k = 1 n W X , k = j = 1 m W Y , j ,
where W X , k are the X-simple regions, and W Y , j are the Y-simple regions. Then, we have
W y z = k = 1 n W y z , k , W x z = j = 1 m W x z , j ,
where W y z and W x z are projections of the region W into Y Z and X Z planes, W y z , k are projections of the X-simple region W X , k , and W x z , j are projections of the Y-simple region W Y , j .
Then the following equations are proved similarly
I W ( M ) [ x , y , z ] D ( K 1 ) x , [ x ] F x ( x , y , z ) = k = 1 n I W y z , k ( M ) [ y , z ] F x ( x k , 2 ( y , z ) , y , z ) F x ( x k , 1 ( y , z ) , y , z ) ,
I W ( M ) [ x , y , z ] D ( K 2 ) y , [ y ] F y ( x , y , z ) = j = 1 m I W x z , j ( M ) [ x , z ] F y ( x , y j , 2 ( x , z ) , z ) F y ( x , y j , 1 ( x , z ) , z ) .
Therefore, we get
I W ( M ) [ x , y , z ] D ( K 1 ) x , [ x ] F x ( x , y , z ) = I S ( M ) [ y , z ] F x ( x , y , z ) ,
I W ( M ) [ x , y , z ] D ( K 2 ) y , [ y ] F y ( x , y , z ) = I S ( M ) [ x , y ] F y ( x , y , z ) .
Using that general fractional divercence is defined as
Div W ( K ) F = D ( K 1 ) x , [ x ] F x ( x , y , z ) + D ( K 2 ) y , [ y ] F y ( x , y , z ) + D ( K 3 ) z , [ z ] F z ( x , y , z ) ,
we can get the general fractional Gauss equation.
As a result, we derive the equation
I W ( M ) [ x , y , z ] Div W ( K ) F = I W , F .
that can be written as
I W ( M ) [ x , y , z ] D ( K 1 ) x , [ x ] F x ( x , y , z ) + D ( K 2 ) y , [ y ] F y ( x , y , z ) + D ( K 3 ) z , [ z ] F z ( x , y , z ) =
I S ( M ) [ y , z ] F x ( x , y , z ) + I S ( M ) [ x , y ] F y ( x , y , z ) + I S ( M ) [ x , y ] F z ( x , y , z ) .
 □
Remark 24.
The general fractional Gauss theorem is proved similarly for region W R 0 , + 3 that can be represented as unions of the Z-simple regions W k , which are piecewise Y-simple and X-simple regions in R 0 , + 3 .

9.4. General Fractional Gauss Theorem for Parallelepiped

The standard Gauss theorem (the Gauss-Ostrogradsky theorem) states the folllwing. Let W be a region in R 0 , + 3 with boundary W . Then the volume integral of the divergence of vector field F over W and the surface integral of F over the boundary W are related by
I W 1 , F = I W 1 Div F .
For the parallelepiped region, a general fractional Gauss theorem can be formulated in the following form.
Theorem 18.
(Fractional Gauss’s Theorem for a Parallelepiped)
Let F x ( x , y , z ) , F y ( x , y ) , F z ( x , y , z ) belong to the function space F 1 1 ( R + 3 ) , and the region W R 0 , + 3 has the form of the parallelepiped
W : = { ( x , y , z ) : 0 a x b , 0 c y d , 0 e z f } .
If the boundary of W be a closed surface S = W , then
I W ( M ) , F = I W ( M ) Div W ( K ) F .
Proof. 
For Cartesian coordinates, we have the vector field F = F x e x + F y e y + F z e z , and the GFI operators
I W ( M ) = I W ( M ) [ x , y , z ] , I W ( M ) = e x I W ( M ) [ y , z ] + e y I W ( M ) [ x , z ] + e z I W ( M ) [ x , y ] .
Then
I W ( M ) , F = I S y z ( M ) [ y , z ] F x + I S x z ( M ) [ x , z ] F y + I S x y ( M ) [ x , y ] F z ,
and
I W ( M ) Div W ( K ) F = I W ( M ) [ x , y , z ] D S ( K ) [ x ] F x + D S ( K ) [ y ] F y + D S ( K ) [ z ] F z ,
where S x y , S x z , S y z are projections of S = W into X Y , X Z , Y Z planes.
If W is parallelepiped (165), then GFI operators (167) are
I W ( M ) [ x , y , z ] = I [ a , b ] ( M 1 ) [ x ] I [ c , d ] ( M 2 ) [ y ] I [ e , f ] ( M 3 ) [ z ] ,
and
I W ( M ) [ y , z ] = I [ c , d ] ( M 2 ) [ y ] I [ e , f ] ( M 3 ) [ z ] ,
I W ( M ) [ x , z ] = I [ a , b ] ( M 1 ) [ x ] I [ e , f ] ( M 3 ) [ z ] ,
I W ( M ) [ x , y ] = I [ a , b ] ( M 1 ) [ x ] I [ c , d ] ( M 2 ) [ y ] .
Using the fundanental theorem of GFC, we can realize the following transformations
I W ( M ) , F = I W ( M ) [ y , z ] F x + I W ( M ) [ z , x ] F y + I W ( M ) [ x , y ] F z =
I [ c , d ] ( M 2 ) [ y ] I [ e , f ] ( M 3 ) [ z ] F x ( b , y , z ) F x ( a , y , z ) +
I [ a , b ] ( M 1 ) [ x ] I [ e , f ] ( M 3 ) [ z ] F y ( x , d , z ) F y ( x , c , z ) +
I [ a , b ] ( M 1 ) [ x ] I [ c , d ] ( M 2 ) [ y ] F z ( x , y , f ) F z ( x , y , e ) =
I [ a , b ] ( M 1 ) [ x ] I [ c , d ] ( M 2 ) [ y ] I [ e , f ] ( M 3 ) [ z ] { D ( K 1 ) x , [ x ] F x ( x , y , z ) +
D ( K 2 ) y , [ y ] F y ( x , y , z ) + D ( K 3 ) z , [ z ] F z ( x , y , z ) } =
I W ( M ) ( D W ( M ) , F ) = I W ( M ) Div W ( K ) F .
This ends the proof of the fractional Gauss’s formula for parallelepiped region. □

10. Equalities for General Fractional Vector Operators and Interpretations

10.1. Equalities for General Fractional Differential Vector Operators

Let us give the basic relations for the differential vector operators of general fractional calculus.
  • For the scalar field U = U ( x , y , x ) , we have
    Div W ( K ) Grad W α U = j = 1 3 D ( K j ) , x j [ x j ] D ( K j ) , x j [ x j ] U .
    In the general case,
    D ( K j ) , x j [ x j ] D ( K j ) , x j [ x j ] U D ( K j K j ) , x j [ x j ] U .
    This is due to the fact that the general fractional integral and the first-order derivative (in definitions of GFDs) do not commute
    D ( K ) , x [ x ] = I ( K ) x [ x ] f ( 1 ) ( x ) d d x I ( K ) x [ x ] f ( x ) = D ( K ) x [ x ] .
    For example, in the case of power-law kernels K ( x ) = h 1 α ( x ) , we have [31] the inequality
    D C , 0 + α 2 D C , 0 + 2 α .
    As a result, we have the definition of the “scalar” general fractional Laplacian
    Δ W ( K ) U = , Div W ( K ) Grad W α U .
  • For the vector field F = F ( x , y , x ) , we have
    Grad W α Div W ( K ) F = m = 1 3 j = 1 3 e m D ( K m ) , x m [ x m ] D ( K j ) , x j [ x j ] F j .
  • The second relation for the scalar field U = U ( x , y , z ) is
    Curl W ( K ) Grad W ( K ) U = e l ε l m n D W , x m [ x m ] D W , x n [ x n ] U = 0 ,
    where ε l m n is Levi-Civita symbol, i.e., it is 1 if ( i , j , k ) is an even permutation of ( 1 , 2 , 3 ) , ( 1 ) if it is an odd permutation, and 0 if any index is repeated. The fulfillment of equality (180) is due to the properties of the Levi-Civita symbol.
  • For the vector field F = e m F m , it is easy to prove the relation
    Div W ( K ) Curl W ( K ) F ( x , y , z ) = k = 1 3 m = 1 3 D ( K k ) , x k [ x k ] ε k l m D ( K m ) , x m [ x m ] F m ( x , y , z ) = k = 1 3 m = 1 3 ε k l m D ( K k ) , x k [ x k ] D ( K m ) , x m [ x m ] F m ( x , y , z ) = 0 ,
    where we use antisymmetry of ε k l m with respect to l and m. Equality (181) is also satisfied by the properties of the Levi-Civita symbol.
  • There exists a relation for the double curl operator in the form
    Curl W ( K ) Curl W ( K ) F ( x , y , z ) = e l ε l m n D , x m ( K m ) [ x m ] ε n p q D ( K p ) , x p [ x p ] F q ( x , y , z ) = e l ε l m n ε n p q D , x m ( K m ) [ x m ] D ( K p ) , x p [ x p ] F q ( x , y , z ) .
    Using the relation
    ε l m n ε l p q = δ m p δ n q δ m q δ n p ,
    we obtain
    Curl W ( K ) Curl W ( K ) F = Grad W ( K ) Div W ( K ) F Δ W ( K ) F .
    As a result, we have the definition of the “vector” general fractional Laplacian
    Δ W ( K ) F = , Grad W ( K ) Div W ( K ) F Curl W ( K ) Curl W ( K ) F
    that is a generalization of the standard vector Laplacian [69,70].
  • In the general case,
    D ( K ) , x [ x ] ( f ( x ) g ( x ) ) ( D ( K ) , x [ x ] f ( x ) ) g ( x ) + f ( x ) ( D ( K ) , x [ x ] g ( x ) ) .
    For example (see Theorem 15.1 of [1]), if f ( x ) and g ( x ) are analytic functions on [ a , b ] and the kernel K ( x ) = h 1 α ( x ) , then product rule for the Riemann-Liouville fractional deirvative has the form
    D R L , 0 + x [ x ] ( f ( x ) g ( x ) ) = j = 0 a ( α , j ) ( D R L , 0 + α j [ x ] f ( x ) ) ( g ( j ) ( x ) ) ,
    where
    a ( α , j ) = Γ ( α + 1 ) Γ ( j + 1 ) Γ ( α j + 1 ) .
    As a result, we have
    Grad W ( K ) f g Grad W ( K ) f g + f Grad W ( K ) g ,
    Div W ( K ) U F Grad W ( K ) U F + U Div W ( K ) F .
    Curl W ( K ) U F Grad W ( K ) U F + U Curl W ( K ) F .
    These relations state that we cannot use the standard Leibniz rule (the product rule) in the general fractional vector calculus.

10.2. Physical Interpretations of General Fractional Differential Vector Operators

The general fractional Grad, Curl, and Div operators (in the Regional, Surface and Line form) are nonlocal characteristics of scalar and vector fields in non-local media and continua. Therefore, these characteristics depend on the region, surface and line that are used in considerations in models of non-local media and systems.
Let us give some basic interpretations of general fractional divergence, gradient and curl operators.
(1)
The general fractional gradient specifies the direction of the fastest increase in the scalar field U ( x , y , z ) in a nonlocal medium. The length of the vector of the line general fractional gradient is equal to the rate of increase of this field in this direction, where this increase takes into account changes in this field at the previous points of the line. In other words, the general fractional gradient describes the increase of scalar field, that takes into account all values of all velocities of the highest integer orders (jet) with some weights [71].
(2)
The general fractional curl operators can be interpreted in the following representation. If the velocity field of a non-local medium is used as a vector field in a region, on a surface or a line, the Curl operators characterize the rotational component of this vector field. The regional, surface and line operators Curl characterize the degree (measure) of non-potentiality of a vector field in nonlocal media and systems. The equality to zero of the general fractional curl of a vector field in a certain region, surface or line means the potentiality of the vector field in a nonlocal medium.
(3)
The general fractional divergence characterizes the degree (measure) of non-solenoidality of a vector field in a nonlocal media and systems. The general fractional divergence in a certain area of space characterizes the amount of a nonlocal medium that arises or disappears within the considered area per unit time. We can say that the general fractional divergence of the vector field describes the power of sources and sinks in a nonlocal medium. The equality to zero of the general fractional divergence in a certain region of space means the absence of sources and sinks of the vector field in a nonlocal medium. In other words, the media or fields in the considered region of space do not disappear and does not arise.

11. General FVC for Orthogonal Curvilinear Coordinates

11.1. Orthogonal Curvilinear Coordinates (OCC)

Curvilinear system of coordinates, or curvilinear coordinates, is a coordinate system in the Euclidean space. In standard vector calculus. the curvilinear coordinates are usually used on a plane ( n = 2 ) and in space ( n = 3 ). For such systems, the coordinate lines may be curved.
In Euclidean space, the use of orthogonal curvilinear coordinates (OCC) is of particular importance, since the formulas related look simpler in orthogonal coordinates than in the general case. The orthogonality can simplify the calculations. The well-known examples of such curvilinear coordinate systems in three-dimensional Euclidean space are cylindrical and spherical coordinates.
The curvilinear coordinates ( q 1 , q 2 , q 3 ) may be derived from a set of Cartesian coordinates ( x , y , z ) by using nonlinear coordinate transformations. It should be emphasized that the violation of the standard chain rule leads to the fact that general fractional vector operators defined in different coordinate systems (Cartesian, cylindrical and spherical) are not related to each other by coordinate transformations. Due to this, it is impossible to obtain general fractional vector integral and differential operators in spherical and cylindrical coordinates by using the coordinate transformation. Therefore, the definitions of general fractional vector operators should be formulated separately.
The specific form of standard vector differential operators may differ, but these forms will be equivalent due to the Leibniz (product) rule. For fractional and generalized fractional calculus, the standard Leibniz rule does not hold. Due to this, such forms of notation cannot be equivalent. In this case, the specific form of the cylindrical and spherical general fractional operators must be such that the theorems of Green, Stock and Gauss hold.
It is known that the expressions for the gradient, divergence, curl and line, surface and volume integrals can be directly expressed. For orthogonal curvilinear coordinates (OCC), these integral and differential operators of the vector calculus can be expressed through the functions:
H i ( q ) = x q 1 2 + y q 2 2 + z q 3 2 .
The positive values H i = H i ( q ) = H i ( q 1 , q 2 , q 3 ) , which depend on a point in R 0 , + 3 , are called the Lame coefficients or scale factors.
For the cylindrical coordinates q 1 = r , q 2 = ϕ , q 3 = z , the Lame coefficients are
H 1 = H r = 1 , H 2 = H ϕ = r , H 3 = H z = 1 ,
where r [ 0 , ) , ϕ [ 0 , 2 π ] , z R .
For the spherical coordinates q 1 = r , q 2 = θ , q 3 = ϕ , the Lame coefficients are
H 1 = H r = 1 , H 2 = H θ = r , H 3 = H ϕ = r sin θ ,
where r [ 0 , ) , θ [ 0 , π ] , ϕ [ 0 , 2 π ] .

11.2. General Fractional Vector Differential Operators in OCC

In this section, we give only equations that will be used in definitions of the general fractional gradient, divergence and curl operators in OCC, and their examples for spherical and cylindrical coordinates. Complete definitions of these operators with function sets for which these operators are defined will be given in the following sections. In this section, formulas will be given only for regional GF vector differential operators. Definitions of line and surface general fractional gradient, divergence and curl operators in OCC will be given in the following sections.

11.2.1. General Fractional Grad, Div, Curl in OCC

The regional general fractional gradient for orthogonal curvilinear coordinates (the GF Gradient in OCC) is expressed in the form
Grad W ( K ) U = e 1 Grad W ( K ) U 1 + e 2 Grad W ( K ) U 2 + e 3 Grad W ( K ) U 3
where
Grad W ( K ) U 1 : = 1 H 1 ( q ) D ( K 1 ) q 1 , [ q 1 ] U ( q 1 , q 2 , q 3 ) ,
Grad W ( K ) U 2 : = 1 H 2 ( q ) D ( K 2 ) q 2 , [ q 2 ] U ( q 1 , q 2 , q 3 ) ,
Grad W ( K ) U 3 : = 1 H 3 ( q ) D ( K 3 ) q 3 , [ q 3 ] U ( q 1 , q 2 , q 3 ) ,
where H k = H k ( q ) = H k ( q 1 , q 2 , q 3 ) , and e 1 = e q 1 , e 2 = e q 2 , e 3 = e q 3 , and the function U ( q ) belongs to the set C 1 1 ( R + 3 ) .
The regional general fractional divergence in orthogonal curvilinear coordinates (the GF Divergence in OCC) is expressed in the form
Div W ( K ) F = Div W ( K ) F 1 + Div W ( K ) F 2 + Div W ( K ) F 3 ,
where
Div W ( K ) F 1 : = 1 H 1 H 2 H 3 D ( K 1 ) q 1 , [ q 1 ] F 1 H 2 H 3 ( q 1 , q 2 , q 3 ) ,
Div W ( K ) F 2 : = 1 H 1 H 2 H 3 D ( K 1 ) q 1 , [ q 1 ] F 2 H 1 H 2 ( q 1 , q 2 , q 3 ) ,
Div W ( K ) F 3 : = 1 H 1 H 2 H 3 D ( K 1 ) q 1 , [ q 1 ] F 3 H 1 H 2 ( q 1 , q 2 , q 3 ) ,
if F ( q ) = F 1 ( q ) e 1 + F 2 ( q ) e 2 + F 3 ( q ) e 3 belong to the space C 1 1 ( R + 3 ) .
The regional general fractional curl in orthogonal curvilinear coordinates (the GF Curl in OCC) is expressed by the eqations
Curl W ( K ) F = e 1 Curl W ( K ) F 1 + e 2 Curl W ( K ) F 2 + e 3 Curl W ( K ) F 3 ,
where
Curl W ( K ) F 1 = 1 H 2 H 3 D ( K 2 ) q 2 , [ q 2 ] F 3 H 3 D ( K 3 ) q 3 , [ q 3 ] F 2 H 2 ,
Curl W ( K ) F 2 = 1 H 1 H 3 D ( K 3 ) q 3 , [ q 3 ] F 1 H 1 D ( K 1 ) q 1 , [ q 1 ] F 3 H 3 ,
Curl W ( K ) F 3 = 1 H 1 H 2 D ( K 1 ) q 1 , [ q 1 ] F 2 H 2 D ( K 2 ) q 2 , [ q 2 ] F 1 H 1 .

11.2.2. General Fractional Grad, Div, Curl in Spherical Coorditates

The Cartesian coordinates ( x , y , z ) and spherical coordinates ( r , θ , φ ) are connected by the equations
x = r sin θ cos φ , y = r sin θ sin φ , z = r cos θ ,
where r [ 0 , ) is the length of the radial vector connecting the origin to the point P ( x , y , z ) , θ [ 0 , π ] is the polar angle, φ [ 0 , 2 π ] is the azimuthal angle. The basic vectors of these coordinate systems
e x = e r sin θ cos φ + e θ cos θ cos φ e φ sin φ e y = e r sin θ sin φ + e θ cos θ sin φ + e φ cos φ e z = e r cos θ e θ sin θ .
Let us emphasize that the violation of the standard chain rule leads to the fact that general fractional vector operators defined in Cartesian, and spherical coordinate are not related to each other by coordinate transformations. Because of this, it is impossible to obtain the fractional vector integral and differential operators in spherical coordinated by using coordinate transformation (191). Therefore, the definitions of general fractional vector operators should be formulated separately.
Let us define the regions
Ball + 3 = { ( r , θ , φ ) : r ( 0 , ) , θ ( 0 , π ] , φ ( 0 , 2 π ] } ,
Ball 0 , + 3 = { ( r , θ , φ ) : r [ 0 , ) , θ [ 0 , π ] , φ [ 0 , 2 π ] } .
We can see that
Ball + 3 R 3 + , Ball 0 , + 3 R 0 , + 3 .
We will consider the vector field
F = e r F r ( r , θ , φ ) + e θ F θ ( r , θ , φ ) + e φ F φ ( r , θ , φ )
that belongs to F 1 1 ( R + 3 ) .
The general fractional gradient, divergence and curl operators in spherical coordinates are defined in the following form. For simplicity, we only present the definitions of regional operators. Linear and surface general fractional operators in spherical coordinates are defined similarly.
Definition 46.
If F ( r , θ , φ ) C 1 1 ( R + 3 ) , then the regional general fractional gradient in spherical coordinates for the region W = Ball 0 , + 3 is defined as
Grad W ( K ) F ( r , θ , φ ) = e r Grad W ( K ) F r + e θ Grad W ( K ) F θ + e φ Grad W ( K ) F φ ,
where
Grad W ( K ) F r ( r , θ , φ ) = D ( K r ) r , [ r ] F ( r , θ , φ )
Grad W ( K ) F θ ( r , θ , φ ) = 1 r D ( K θ ) θ , [ θ ] F ( r , θ , φ ) ,
Grad W ( K ) F φ ( r , θ , φ ) = 1 r sin θ D ( K φ ) φ , [ φ ] F ( r , θ , φ ) .
Definition 47.
If F ( r , θ , φ ) C 1 1 ( R + 3 ) , then the regional general fractional divergence in spherical coordinates for the region W = Ball 0 , + 3 is defined as
Div W ( K ) F = Div W ( K ) F r + Div W ( K ) F θ + Div W ( K ) F φ ,
where
Div W ( K ) F r ( r , θ , φ ) = 1 r 2 D ( K r ) r , [ r ] ( r ) 2 F r ( r , θ , φ )
Div W ( K ) F θ ( r , θ , φ ) = 1 r sin θ D ( K θ ) θ , [ θ ] F θ ( r , θ , φ ) sin θ ,
Div W ( K ) F φ ( r , θ , φ ) = 1 r sin θ D ( K φ ) φ , [ φ ] F φ ( r , θ , φ ) .
Let us note the violations of the standard product (Leibniz) rule. For example, we have the inequality
D ( K θ ) θ , [ θ ] F θ ( r , θ , φ ) sin θ
D ( K θ ) θ , [ θ ] F θ ( r , θ , φ ) sin θ + F θ ( r , θ , φ ) D ( K θ ) θ , [ θ ] sin θ .
As a second example
1 r 2 D ( K r ) r , [ r ] ( r ) 2 F r ( r , θ , φ )
D ( K r ) r , [ r ] F r ( r , θ , φ ) + 2 r F θ ( r , θ , φ ) .
Definition 48.
If F ( r , θ , φ ) C 1 1 ( R + 3 ) , then the regional general fractional curl operator in spherical coordinates for the region W = Ball 0 , + 3 is defined as
Curl W ( K ) F = e r Curl W ( K ) F r + e θ Curl W ( K ) F θ + e φ Curl W ( K ) F φ ,
where
Curl W ( K ) F r ( r , θ , φ ) = 1 r sin θ D ( K θ ) θ , [ θ ] ( F φ ( r , θ , φ ) sin θ ) 1 r sin θ D ( K φ ) φ , [ φ ] F φ ( r , θ , φ ) ,
Curl W ( K ) F θ ( r , θ , φ ) = 1 r sin θ D ( K φ ) φ , [ φ ] F r ( r , θ , φ ) 1 r D ( K r ) r , [ r ] r F φ ( r , θ , φ ) ,
Curl W ( K ) F φ ( r , θ , φ ) = 1 r D ( K r ) r , [ r ] r F θ ( r , θ , φ ) 1 r D ( K θ ) θ , [ θ ] F r ( r , θ , φ ) .

11.2.3. General Fractional Grad, Div, Curl in Cylindrical Coorditates

The Cartesian and cylindrical coordinates
x = r cos φ y = r sin φ z = z ,
where r [ 0 , ) , φ [ 0 , 2 π ] , and z [ 0 , ) . The basic vectors of these coordinate systems
e x = e r cos φ e φ sin φ e y = e r sin φ + e φ cos φ e z = e z .
Note that it is impossible to obtain the fractional vector integral and differential operators in cylindrical coordinated by using coordinate transformation (192). Therefore, the definitions of general fractional vector operators should be formulated separately.
Let us define the regions
Cylinder + 3 = { ( r , φ , z ) : r ( 0 , ) , φ ( 0 , 2 π ] , z ( 0 , i n f t y ) } ,
Cylinder 0 , + 3 = { ( r , φ , z ) : r [ 0 , ) , φ [ 0 , 2 π ] , z [ 0 , i n f t y ) } .
We can see that
Cylinder + 3 R 3 + , Cylinder 0 , + 3 R 0 , + 3 .
Let us consider the vector field
F = e r F r ( r , φ , z ) + e φ F φ ( r , φ , z ) + e z F y ( r , φ , z )
that belongs to F 1 1 ( R + 3 )
The general fractional gradient, divergence and curl operators in cylindrical coordinates are defined in the following form. For simplicity, we give definitions of regional operators. Linear and surface general fractional operators in cylindrical coordinates are defined similarly.
Definition 49.
If F ( r , φ , z ) C 1 1 ( R + 3 ) , then the regional general fractional gradient in cylindrical coordinates for the region W = Cylinder 0 , + 3 is defined as
Grad W ( K ) F ( r , φ , z ) = e r Grad W ( K ) F r + e φ Grad W ( K ) F φ + e z Grad W ( K ) F z ,
where
Grad W ( K ) F r ( r , φ , z ) = D ( K r ) r , [ r ] F ( r , φ , z )
Grad W ( K ) F φ ( r , φ , z ) = 1 r D ( K φ ) φ , [ φ ] F ( r , φ , z ) ,
Grad W ( K ) F z ( r , φ , z ) = D ( K z ) z , [ z ] F ( r , φ , z ) .
Definition 50.
If F ( r , φ , z ) C 1 1 ( R + 3 ) , then the regional general fractional divergence in cylindrical coordinates for the region W = Cylinder 0 , + 3 is defined as
Div W ( K ) F = Div W ( K ) F r + Div W ( K ) F φ + Div W ( K ) F z ,
where
Div W ( K ) F r ( r , φ , z ) = 1 r D ( K r ) r , [ r ] r F r ( r , φ , z )
Div W ( K ) F φ ( r , φ , z ) = 1 r D ( K φ ) φ , [ φ ] F φ ( r , φ , z ) ,
Div W ( K ) F z ( r , φ , z ) = D ( K z ) z , [ z ] F z ( r , φ , z ) .
We should note that
1 r D ( K r ) r , [ r ] r F r ( r , φ , z ) D ( K r ) r , [ r ] F r ( r , φ , z ) + 1 r F r ( r , φ , z )
since the standard Leibniz rule does not hold for fractional derivatives of non-integer order and for general fractional derivatives.
Definition 51.
If F ( r , φ , z ) C 1 1 ( R + 3 ) , then the regional general fractional curl operator in cylindrical coordinates for the region W = Cylinder 0 , + 3 is defined as
Curl W ( K ) F = e r Curl W ( K ) F r + e φ Curl W ( K ) F φ + e z Curl W ( K ) F z ,
where
Curl W ( K ) F r ( r , φ , z ) = 1 r D ( K φ ) φ , [ φ ] F z ( r , φ , z ) D ( K z ) z , [ z ] F z ( r , φ , z ) ,
Curl W ( K ) F φ ( r , φ , z ) = D ( K z ) z , [ z ] F r ( r , φ , z ) D ( K r ) r , [ r ] F z ( r , φ , z ) ,
Curl W ( K ) F z ( r , φ , z ) = 1 r D ( K r ) r , [ r ] r F φ ( r , φ , z ) 1 r D ( K φ ) φ , [ φ ] F r ( r , φ , z ) .

11.3. General Fractional Integral Operators in OCC

11.3.1. GFI in OCC

The following GGI operators can be used to define the line, surface and volume general fractional integrals in orthogonal curvilinear coordinates (the GF integrals in OCC).
Definition 52.
Let f ( q ) = f ( q 1 , q 2 , q 3 ) be function in R 0 , + 3 that satisfies the conditions
H 1 ( q 1 , a , b ) f ( q 1 , a , b ) C 1 ( 0 , ) ,
H 2 ( a , q 2 , b ) f ( a , q 2 , b ) C 1 ( 0 , ) ,
H 3 ( a , b , q 3 ) f ( a , b , q 3 ) C 1 ( 0 , ) ,
for all a , b 0 .
Then the set of such functions is denoted by C H , 1 ( R + 3 ) .
Let function f ( q ) = f ( q 1 , q 2 , q 3 ) belong to the set C H , 1 ( R + 3 ) . For the orthogonal curvilinear coordinate q k ( k = 1 , 2 , 3 ), the GFI operators are defined as
I ( M k ) q k ^ [ q k ] f ( q ) : = I ( M k ) q k [ q k ] H k ( q ) f ( q ) = 0 q k M k ( q k q k ) f ( q ) H k ( q ) d q k ,
and
I ( M k ) ^ [ a k , b k ] [ q k ] f ( q ) : = I ( M k ) b k ^ [ q k ] f ( q ) I ( M k ) a k ^ [ q k ] f ( q ) ,
where notations (193) means the following
I ( M 1 ) q 1 ^ [ q 1 ] f ( q 1 , q 2 , q 3 ) = I ( M 1 ) q 1 [ q 1 ] H 1 ( q 1 , q 2 , q 3 ) f ( q 1 , q 2 , q 3 ) =
0 q 1 M 1 ( q 1 q 1 ) f ( q 1 , q 2 , q 3 ) H 1 ( q 1 , q 2 , q 3 ) d q 1 ,
I ( M 2 ) q 2 ^ [ q 2 ] f ( q 1 , q 2 , q 3 ) : = I ( M 2 ) q 2 [ q 2 ] H 2 ( q 1 , q 2 , q 3 ) f ( q 1 , q 2 , q 3 ) =
0 q 2 M 2 ( q 2 q 2 ) f ( q 1 , q 2 , q 3 ) H 2 ( q 1 , q 2 , q 3 ) d q 2 ,
I ( M 3 ) q 3 ^ [ q 3 ] f ( q 1 , q 2 , q 3 ) : = I ( M 3 ) q 3 [ q 3 ] H 3 ( q 1 , q 2 , q 3 ) f ( q 1 , q 2 , q 3 ) =
0 q 3 M 3 ( q 3 q 3 ) f ( q 1 , q 2 , q 3 ) H 3 ( q 1 , q 2 , q 3 ) d q 3 .

11.3.2. GFI in Spherical Coordinates

To define the line spherical GFI, we can use the following GFI operators.
(1)
The radial GFI operator
I ( M r ) r ^ [ r ] f ( r , θ , φ ) : = I ( M r ) r [ r ] r f ( r , θ , φ ) = 0 r M r ( r r ) f ( r , θ , φ ) r d r .
(2)
The polar GFI operator
I M θ θ ^ [ θ ] f ( , θ , φ ) = I M θ θ [ θ ] r f ( r , θ , φ ) =
0 θ M θ ( θ θ ) f ( r , θ , φ ) r d θ = r 0 θ M θ ( θ θ ) f ( r , θ , φ ) d θ .
(3)
The azimuthal GFI operator
I ( M φ ) φ ^ [ φ ] f ( r , θ , φ ) = I ( M φ ) φ [ φ ] r sin θ f ( r , θ , φ ) =
0 φ M φ ( φ φ ) f ( r , θ , φ ) r sin θ d φ = r sin θ 0 φ M φ ( φ φ ) f ( r , θ , φ ) d φ .
We can define the operators on the positive intervals
I [ r 0 , r ] ( M r ) ^ [ r ] f ( r , θ , φ ) = I ( M r ) r ^ [ r ] f ( r , θ , φ ) I ( M r ) r 0 ^ [ r ] f ( r , θ , φ )
for r 0 > 0 . For r 0 = r 1 > 0 , r = r 1 , we have
I [ r 1 , r 2 ] ( M r ) ^ [ r ] f ( r , θ , φ ) = I ( M r ) r 2 [ r ] ^ f ( r , θ , φ ) I ( M r ) r 1 [ r ] ^ f ( r , θ , φ ) .
Integrals with respect to other variables are defined in a similar way.
I [ θ 0 , θ ] M θ ^ [ θ ] f ( r , θ , φ ) = I M θ θ [ θ ] ^ f ( r , θ , φ ) I M θ θ 0 ^ [ θ ] f ( r , θ , φ ) ,
I [ θ 1 , θ 2 ] M θ ^ [ θ ] f ( r , θ , φ ) = I M θ θ 2 [ θ ] ^ f ( r , θ , φ ) I M θ θ 1 ^ [ θ ] f ( r , θ , φ ) ,
and
I [ φ 0 , φ ] ( M φ ) ^ [ φ ] f ( r , θ , φ ) = I ( M φ ) φ ^ [ φ ] f ( r , θ , φ ) I ( M φ ) φ 0 ^ [ φ ] f ( r , θ , φ ) ,
I [ φ 1 , φ 1 ] ( M φ ) ^ [ φ ] f ( r , θ , φ ) = I ( M φ ) φ 2 ^ [ φ ] f ( r , θ , φ ) I ( M φ ) φ 1 ^ [ φ ] f ( r , θ , φ ) .
For intervals r 0 , θ 0 = 0 , and φ = 0 , we define
I [ 0 , r ] ( M r ) ^ [ r ] : = I ( M r ) r ^ [ r ] , I [ 0 , θ ] ( M θ ) ^ [ θ ] : = I ( M θ ) θ ^ [ θ ] , I [ 0 , φ ] ( M φ ) ^ [ φ ] : = I ( M φ ) φ ^ [ φ ] .

11.3.3. GFI in Cylindrical Coordinates

Let us define the cylindrical GFI, we can use the following GFI operators.
(1)
If f ( r ) C 1 ( 0 , ) , then
I ( M r ) r ^ [ r ] f ( r ) : = I ( M r ) r [ r ] f ( r ) : = 0 r M r ( r r ) f ( r ) d r
and
I [ r 1 , r 2 ] ( M r ) ^ [ r ] f ( r ) : = I ( M r ) r 2 ^ [ r ] f ( r ) I ( M r ) r 1 ^ [ r ] f ( r ) =
0 r 2 M r ( r 2 r ) f ( r ) d r 0 r 1 M r ( r 1 r ) f ( r ) d r ,
where r 2 > r 1 > 0 , and I [ 0 , r 2 ] ( M r ) ^ [ r ] : = I ( M r ) r 2 ^ [ r ] f ( r ) .
(2)
If f ( φ ) C 1 ( 0 , ) , then
I ( M φ ) φ ^ [ φ ] f ( φ ) : = I ( M φ ) φ [ φ ] r f ( φ ) : =
0 φ M φ ( φ φ ) f ( φ ) r d φ = r 0 φ M φ ( φ φ ) f ( φ ) d φ ,
and
I [ φ 1 , φ 2 ] ( M φ ) ^ [ φ ] f ( φ ) : = I ( M φ ) φ 2 ^ [ φ ] f ( φ ) I ( M φ ) φ 1 ^ [ φ ] f ( φ ) ,
if 0 < φ 1 < φ 2 < 2 π and
I [ 0 , φ 2 ] ( M φ ) ^ [ φ ] f ( φ ) : = I ( M φ ) φ 2 ^ [ φ ] f ( φ ) .
(3)
If f ( z ) C 1 ( 0 , ) , then
I ( M z ) z ^ [ z ] f ( z ) = 0 z M z ( z z ) f ( z ) d z .
and
I [ z 1 , z 2 ] M z ^ [ z ] f ( z ) = I ( M z ) z 2 ^ [ z ] f ( z ) I ( M z ) z 1 ^ [ z ] f ( z ) ,
if z 1 > 0 and
I [ 0 , z 2 ] ( M z ) ^ [ z ] f ( z ) = I ( M z ) z 2 ^ [ z ] f ( z ) .

11.4. General Fractional Operators in Curvilinear Coordinates

11.4.1. Definition of Line GFI for Vector Field in OCC

Let us define the line GFI in orthogonal curvilinear coordinates (OCC) in R 0 , + 2 of the Q 1 Q 2 -plane.
Definition 53.
Let L be a simple line in R 0 , + 2 of the Q 1 Q 2 -plane. Let the functions
f 1 ( q 1 ) : = H 1 ( q 1 , y ( q 1 ) ) F 1 ( q 1 , y ( q 1 ) ) , f 2 ( q 2 ) : = H 2 ( x ( q 2 ) , q 2 ) F 2 ( x ( q 2 ) , q 2 )
belong to the set C 1 ( 0 , ) .
Then line GFI for the line L is defined by the equation
I ^ L ( M ) , F = I [ a , b ] ( M 1 ) [ q 1 ] f 1 ( q 1 ) + I [ c , d ] ( M 2 ) [ q 2 ] f 2 ( q 2 ) = I [ a , b ] ( M 1 ) [ q 1 ] H 1 ( q 1 , y ( q 1 ) ) F 1 ( q 1 , y ( q 1 ) ) + I [ c , d ] ( M 2 ) [ q 2 ] H 2 ( x ( q 2 ) , q 2 ) F y ( x ( q 2 ) , q 2 ) .
Line GFI (194) exists, if the kernel pairs ( M 1 ( q 1 ) , K 1 ( q 1 ) ) and ( M 2 ( q 2 ) , K 2 ( q 2 ) ) belong to the Luchko set L .
Let us define the line GFI in orthogonal curvilinear coordinates (OCC) in R 0 , + 3 of the Q 1 Q 2 Q 3 -space. To give this definition, we will describe conditions on the vector field F ( q 1 , q 2 , q 3 ) . We will assume that the vector field
F = F ( q ) = e 1 F 1 ( q 1 , q 2 , q 3 ) + e 2 F 2 ( q 1 , q 2 , q 3 ) + e 3 F 3 ( q 1 , q 2 , q 3 )
on the simple line L R 0 , + 3 is described by the following functions that belong to the space C 1 ( 0 , ) :
f 1 ( q 1 ) : = H 1 ( q 1 , y ( q 1 ) , z ( q 1 ) ) F 1 ( q 1 , y ( q 1 ) , z ( q 1 ) ) C 1 ( 0 , ) ,
f 2 ( q 2 ) : = H 2 ( x ( q 2 ) , q 2 , z ( q 2 ) ) F 2 ( x ( q 2 ) , q 2 , z ( q 2 ) ) C 1 ( 0 , ) ,
f 3 ( q 3 ) : = H 3 ( x ( q 3 ) , y ( q 3 ) , q 3 ) F 3 ( x ( q 3 ) , y ( q 3 ) , q 3 ) C 1 ( 0 , ) .
If these conditions are satisfied, then we will write F ( q ) F H , 1 ( L ) .
In the case F F H , 1 ( R + 3 ) , the line general fractional integral for the vector field F and the line L = A B with endpoints A ( a 1 , a 2 , a 3 ) and B ( b 1 , b 1 , b 3 ) with all b k > a k 0 is defined by the equation
I ^ L ( M ) , F = k = 1 3 I [ a k , b k ] ( M k ) [ q k ] f k ( q k ) = I [ a 1 , b 1 ] ( M 1 ) [ q 1 ] f 1 ( q 1 ) + I [ a 2 , b 2 ] ( M 2 ) [ q 2 ] f 2 ( q 2 ) + I [ a 3 , b 3 ] ( M 3 ) [ q 3 ] f 3 ( q 3 ) = I [ a 1 , b 1 ] ( M 1 ) [ q 1 ] H 1 ( q 1 , y ( q 1 ) , z ( q 1 ) ) F 1 ( q 1 , y ( q 1 ) , z ( q 1 ) ) + I [ a 2 , b 2 ] ( M 2 ) [ q 2 ] H 2 ( x ( q 3 ) , q 2 , z ( q 2 ) ) F 2 ( x ( q 3 ) , q 2 , z ( q 2 ) ) + I [ a 3 , b 3 ] ( M 3 ) [ q 3 ] H 3 ( x ( q 3 ) , y ( q 3 ) , q 3 ) F 3 ( x ( q 3 ) , y ( q 3 ) , q 3 ) .
Using the Laplace convolution, the line GFI can be written as
I ^ L ( M ) , F = k = 1 3 ( M k f k ) ( b k ) ( M k f k ) ( a k ) ,
where b k > a k 0 , and
( M k f k ) ( c k ) = 0 c k d q k M ( c k q k ) f k ( q k ) ,
where c k = a k or c k = b k with M k ( q k ) C 1 , 0 ( 0 , ) and f k ( q k ) C 1 ( 0 , ) .

11.4.2. Line GFI for Piecewise Simple Lines in OCC

Let us consider a line L R 0 , + 3 , which can be divided into several lines L k = L k [ A k , A k + 1 ] , k = 1 , , n that are simple lines or lines parallel to one of the axes:
L : = k = 1 n 1 L k [ A k , A k + 1 ] ,
where the line L k connects the points A k ( x k , y k , z k ) , and A k + 1 ( x k + 1 , y k + 1 , z k + 1 ) with x k , y k , z k 0 for all k = 1 , 2 , , n 1 . Lines of this kind will be called the piecewise simple lines.
For piecewise simple line (196) in OCC of R 0 , + 3 , and the vector field F F 1 ( L ) , the line GFI is defined by the equation
I ^ L ( M ) [ q 1 , q 2 , q 3 ] , F ( q 1 , q 2 , q 3 ) : = k = 1 n 1 I [ A k , A k + 1 ] ( M ) [ q 1 , q 2 , q 3 ] , F ( q 1 , q 2 , q 3 ) ,
where
I ^ [ A k , A k + 1 ] ( M ) [ q 1 , q 2 , q 3 ] , F ( q 1 , q 2 , q 3 ) : =
sgn ( q 1 , k + 1 q 1 , k ) I ω [ q 1 , k , q 1 , k + 1 ] ( M 1 ) [ q 1 ] H 1 ( q 1 , y k ( q 1 ) , z k ( q 1 ) ) F 1 ( q 1 , y k ( q 1 ) , z k ( q 1 ) ) +
sgn ( q 2 , k + 1 q 2 , k ) I ω [ q 2 , k , q 2 , k + 1 ] ( M 2 ) [ q 2 ] H 2 ( x k ( q 2 ) , q 2 , z k ( q 2 ) ) F 2 ( x k ( q 2 ) , q 2 , z k ( q 2 ) ) +
sgn ( q 3 , k + 1 q 3 , k ) I ω [ q 3 , k , q 3 , k + 1 ] ( M 3 ) [ q 3 ] H 3 ( x k ( q 3 ) , y k ( q 3 ) , q 3 ) F 3 ( x k ( q 3 ) , y k ( q 3 ) , q 3 ) ,
where the functions y = y k ( q 1 ) , x = x k ( q 2 ) , z = z k ( q 1 ) , x = x k ( q 3 ) , z = z k ( q 2 ) , y = y k ( q 1 ) define the lines L k that are simple or parallel to one of the axes of OCC.

11.4.3. Regional GF Gradient in OCC

Let us give definitions of a set of scalar fields and a general fractional gradient in OCC for R 0 , + 3 .
Definition 54.
Let U ( q ) = U ( q 1 , q 2 , q 3 ) be a scalar field that satisfies the conditions
D ( K 1 ) q 1 , [ q 1 ] U ( q 1 , q 2 , q 3 ) C 1 ( R + 3 ) ,
D ( K 2 ) q 2 , [ q 2 ] U ( q 1 , q 2 , q 3 ) C 1 ( R + 3 ) ,
D ( K 3 ) q 3 , [ q 3 ] U ( q 1 , q 2 , q 3 ) C 1 ( R + 3 ) .
Then the set of such scalar fields U ( q ) = U ( q 1 , q 2 , q 3 ) will be denoted as C 1 1 ( R + 3 ) .
Definition 55.
Let U ( q ) = U ( q 1 , q 2 , q 3 ) be a scalar field that belongs to the set C 1 1 ( R + 3 ) .
Then the general fractional gradient Grad W ( K ) in OCC for the region W = R 0 , + 3 is defined as
Grad W ( K ) U ( x , y , z ) = e 1 1 H 1 ( q ) D ( K 1 ) q 1 , [ x ] U ( q 1 , q 2 , q 3 ) +
e 2 1 H 2 ( q ) D ( K 2 ) q 2 , [ y ] U ( q 1 , q 2 , q 3 ) + e 3 1 H 3 ( q ) D ( K 3 ) q 3 , [ z ] U ( q 1 , q 2 , q 3 ) .
This operator will be called the regional general fractional gradient (regional GF gradient) in OCC.

11.4.4. Line GF Gradient in OCC for R 0 , + 2

The general fractional gradients (Line GF Gradient) in OCC can be defined not only for the region W = R 0 , + 3 . Using the fact that GFD is integro-differential operator, we can define GF Gradients for line L R 0 , + 2 and L R 0 , + 3 . The general fractional vector operators can be defined as a sequential action of first-order derivatives and general fractional integrals.
Definition 56.
Let L be a simple line in R 0 , + 2 of the Q 1 Q 2 -plane, and U ( q 1 , q 2 ) F 1 , L 1 ( R + 2 ) that means
U q 1 ( 1 ) ( q 1 , y ( q 1 ) ) C 1 ( 0 , ) , U q 2 ( 1 ) ( x ( q 2 ) , q 1 ) C 1 ( 0 , ) .
Then, the line general fractional gradient in OCC (Line GF Gradient in OCC) for the line L is defined by the equation
Grad L ( K ) U ( q 1 , q 2 ) : =
e 1 1 H 1 ( q ) I ( K 1 ) q 1 , [ q 1 ] U q 1 ( 1 ) ( q 1 , y ( q 1 ) ) + e 2 1 H 2 ( q ) I ( K 2 ) q 2 , [ q 2 ] U q 2 ( 1 ) ( x ( q 2 ) , q 1 ) =
e 1 1 H 1 ( q 1 , q 2 ) I ( K 1 ) q 1 , [ q 1 ] U ( q 1 , q 2 ) q 1 q 2 = y ( q 1 ) + e 2 1 H 2 ( q 1 , q 2 ) I ( K 2 ) q 1 , [ q 2 ] U ( q 1 , q 2 ) q 2 q 1 = x ( q 2 ) ,
where the kernels ( M 1 ( q 1 ) , K 1 ( q 1 ) ) and ( M 2 ( q 2 ) , K 2 ( q 1 ) ) belong to the Luchko set L .
We can use the definitions of the line GFI and line GF Gradient in OCC to prove the following theorem.
Theorem 19.
(General Fractional Gradient Theorem for line GF Gradient for R 0 , + 2 )
Let L be a simple line in R 0 , + 2 of the Q 1 Q 2 -plane, which connects the points A ( a , b ) and B ( c , d ) , and the scalar field U ( q 1 , q 2 ) belongs to the set F 1 , L 1 ( R + 2 ) .
Then, the equality
I L ( M ) , Grad L ( K ) U = U ( c , d ) U ( a , b )
holds, where b a 0 and d c 0 .

11.4.5. Theorem for Line GF Gradient in OCC of R 0 , + 3

In the theorem we will use the following definitions.
Definition 57.
Let L R 0 , + 3 be a line in OCC that is described by the functions
q 2 = y ( q 1 ) 0 , q 3 = z ( q 1 ) 0 , q 1 [ a , b ] R 0 , + ,
which are continuously differentiable functions for q 1 ( a , b ) R 0 , + , i.e., y ( q 1 ) , z ( q 1 ) C 1 ( a , b ) . Then this line will be called Q 2 Q 3 -simple line. The line L is called simple line in OCC of R 0 , + 3 , if L is Q 1 Q 2 -, Q 1 Q 3 - and Q 2 Q 3 -simple line.
Let us consider a simple line L R 0 , + 3 in OCC, which connects the points A ( a , c , e ) and B ( b , d , f ) , and can be described by following equivalent forms
L = { ( q 1 , q 2 , q 3 ) : q 1 [ a , b ] , q 2 = y ( q 1 ) C 1 [ a , b ] , q 3 = z ( q 1 ) C 1 [ a , b ] } ,
L = { ( q 1 , q 2 , q 3 ) : q 2 [ c , d ] , q 1 = x ( q 2 ) C 1 [ c , d ] , q 3 = z ( q 2 ) C 1 [ c , d ] } ,
L = { ( q 1 , q 2 , q 3 ) : q 3 [ e , f ] , q 1 = x ( Q 3 ) C 1 [ e , f ] , q 2 = y ( q 3 ) C 1 [ e , f ] } ,
where y ( q 1 = a ) = c , y ( 1 1 = b ) = d , z ( q 1 = a ) = e , z ( q 1 = b ) = f .
Definition 58.
Let L be a simple line in R 0 , + 3 , which is defined in form (198)–(200), and a scalar field U ( q 1 , q 2 , q 3 ) satisfies the conditions
U 1 ( x ) : = U x ( 1 ) ( x , y ( x ) , z ( x ) ) = U ( q 1 , q 2 , q 3 ) x y = y ( x ) , z = z ( x ) C 1 n ( 0 , ) ,
U 2 ( y ) : = U y ( 1 ) ( x ( y ) , y , z ( y ) ) = U ( q 1 , q 2 , q 3 ) y x = x ( y ) , z = z ( y ) C 1 n ( 0 , ) ,
U 3 ( z ) : = U z ( 1 ) ( x ( z ) , y ( z ) , z ) = U ( q 1 , q 2 , q 3 ) z x = x ( z ) , y = y ( z ) C 1 n ( 0 , ) .
Then the set of such scalar fields U ( q 1 , q 2 , q 3 ) will be denoted as U ( q 1 , q 2 , q 3 ) F 1 n ( L ) .
Definition 59.
Let L be a simple line in R 0 , + 3 , which is defined in form (198)–(200), and a scalar field U ( q 1 , q 2 , q 3 ) belongs to the set F 1 1 ( L ) .
Then the line general fractional gradient for the line L R 0 , + 3 is defined by the equation
Grad L ( K ) U ( q 1 , q 2 , q 3 ) = D L ( M ) U ( q 1 , q 2 , q 3 ) = 1 H 1 ( q ) I [ a , b ] ( M 1 ) [ q 1 ] U 1 ( 1 ) ( q 1 , y ( q 1 ) , z ( q 1 ) ) + 1 H 2 ( q ) I [ c , d ] ( M 2 ) [ q 2 ] U 2 ( 1 ) ( x ( q 2 ) , q 2 , z ( q 2 ) ) + 1 H 3 ( q ) I [ e , f ] ( M 3 ) [ q 3 ] U 3 ( 1 ) ( x ( q 3 ) , y ( q 3 ) , q 3 ) .
where the pairs of the kernels ( M k ( q k ) , K ( q k ) ) , ( k = 1 , 2 , 3 ) belong to the Luchko set L .
For OCC, the general fractional gradient theorem for the line GF Gradient is fomulated for simple lines in the form.
Theorem 20.
(General Fractional Gradient Theorem for OCC and Line GF Gradient)
Let L be a simple line in OCC for R 0 , + 3 , which is described in form (198)–(200), and connects the points A ( a , c , e ) and B ( b , d , f ) .
Let U ( q 1 , q 2 , q 3 ) be a scalar filed that belongs to the set F 1 , L 1 ( R + 3 ) .
Then, the equality
I ^ L ( M ) , Grad L ( K ) U = U ( a , c , e ) U ( b , d , f )
holds.
Proof. 
Using that the line GFI of a vector field F for the simple line L R 0 , + 3 is defined by the equation
I ^ L ( M ) , F = I [ a , b ] ( M 1 ) [ q 1 ] H 1 ( q 1 , y ( q 1 ) , z ( q 1 ) ) F 1 ( q 1 , y ( q 1 ) , z ( q 1 ) ) +
I [ c , d ] ( M 2 ) [ q 2 ] H 2 ( x ( q 2 ) , q 2 , z ( q 2 ) ) F 2 ( x ( q 2 ) , q 2 , z ( q 2 ) ) +
I [ e , f ] ( M 3 ) [ q 3 ] H 3 ( x ( q 3 ) , y ( q 3 ) , q 3 ) F 3 ( x ( q 3 ) , y ( q 3 ) , q 3 ) ,
we can consider the line GFI of the vector field that is defined by the line general fractional gradient
F 1 ( q 1 , y ( q 1 ) , z ( q 1 ) ) : = Grad L ( K ) U ( q 1 , q 2 , q 3 ) 1 =
1 H 1 ( q 1 , y ( q 1 ) , z ( q 1 ) ) I ( K 1 ) q 1 , [ q 1 ] U q 1 ( 1 ) ( q 1 , y ( q 1 ) , z ( q 1 ) ) ,
F 2 ( x ( q 2 ) , q 2 , z ( q 2 ) ) : = Grad L ( K ) U ( q 1 , q 2 , q 3 ) 2 =
1 H 2 ( x ( q 2 ) , q 2 , z ( q 2 ) ) I ( K 2 ) q 2 , [ q 2 ] U q 2 ( 1 ) ( x ( q 2 ) , q 2 , z ( q 2 ) ) ,
F 3 ( x ( q 3 ) , y ( q 3 ) , q 3 ) : = Grad L ( K ) U ( q 1 , q 2 , q 3 ) 3 =
1 H 3 ( x ( q 3 ) , y ( q 3 ) , q 3 ) I ( K 3 ) q 3 , [ q 3 ] U q 3 ( 1 ) ( x ( q 3 ) , y ( q 3 ) , q 3 ) .
Therefore, we can get the line GFI of the line FG Gradient
I ^ L ( M ) , Grad L ( K ) U =
I [ a , b ] ( M 1 ) ^ [ q 1 ] 1 H 1 ( q 1 , y ( q 1 ) , z ( q 1 ) ) I ( K 1 ) q 1 , [ q 1 ] U q 1 ( 1 ) ( q 1 , y ( q 1 ) , z ( q 1 ) ) +
I [ c , d ] ( M 2 ) ^ [ q 2 ] 1 H 2 ( x ( q 2 ) , q 2 , z ( q 2 ) ) I ( K 2 ) q 2 , [ q 2 ] U q 2 ( 1 ) ( x ( q 2 ) , q 2 , z ( q 2 ) ) +
I [ e , f ] ( M 3 ) ^ [ q 3 ] 1 H 3 ( x ( q 3 ) , y ( q 3 ) , q 3 ) I ( K 3 ) q 3 , [ q 3 ] U q 3 ( 1 ) ( x ( q 3 ) , y ( q 3 ) , q 3 ) =
I [ a , b ] ( M 1 ) [ q 1 ] I ( K 1 ) q 1 , [ q 1 ] U q 1 ( 1 ) ( q 1 , y ( q 1 ) , z ( q 1 ) ) +
I [ c , d ] ( M 2 ) [ q 2 ] I ( K 2 ) q 2 , [ q 2 ] U q 2 ( 1 ) ( x ( q 2 ) , q 2 , z ( q 2 ) ) +
I [ e , f ] ( M 3 ) [ q 3 ] I ( K 3 ) q 3 , [ q 3 ] U q 3 ( 1 ) ( x ( q 3 ) , y ( q 3 ) , q 3 ) ,
where we use
I [ a , b ] ( M 1 ) ^ [ q 1 ] = I [ a , b ] ( M 1 ) [ q 1 ] H 1 ( q 1 , y ( q 1 ) , z ( q 1 ) ) ,
I [ c , d ] ( M 2 ) ^ [ q 2 ] = I [ c , d ] ( M 2 ) [ q 2 ] H 2 ( x ( q 2 ) , q 2 , z ( q 2 ) ) ,
I [ e , f ] ( M 3 ) ^ [ q 3 ] = I [ e , f ] ( M 3 ) [ q 3 ] H 3 ( x ( q 3 ) , y ( q 3 ) , q 3 ) .
Further, transformations are made similar to the transformations performed in the proof of Theorem 20. Using the fact that the kernels ( M k ( q k ) , K k ( q k ) ) , k = 1 , 2 , 3 , belong to the Luchko set L , we obtain
I [ a , b ] ( M 1 ) [ q 1 ] I ( K 1 ) q 1 , [ q 1 ] f 1 ( q 1 ) = ( M 1 ( K 1 f 1 ) ) ( b ) ( M 1 ( K 1 f 1 ) ) ( a ) =
( ( M 1 K 1 ) f 1 ) ( b ) ( ( M 1 K 1 ) f 1 ) ( a ) = ( { 1 } f 1 ) ( b ) ( { 1 } f 1 ) ( a ) =
0 b f 1 ( q 1 ) d q 1 0 a f 1 ( q 1 ) d q 1 = a b f 1 ( q 1 ) d q 1 .
Similarly, we can obtain the equations for q 2 and q 3 .
Therefore, we get
I ^ L ( M ) , Grad L ( K ) U =
a b U ( q 1 , q 2 , q 3 ) q 1 q 2 = y ( q 1 ) , q 3 = z ( q 1 ) d q 1 + c d U ( q 1 , q 2 , q 3 ) q 2 q 1 = x ( q 2 ) , q 3 = z ( q 2 ) d q 2 +
e f U ( q 1 , q 2 , q 3 ) q 3 q 1 = x ( q 3 ) . q 2 = y ( q 3 ) d q 3 =
a b U ( q 1 , q 2 , q 3 ) q 1 q 2 = y ( q 1 ) , q 3 = z ( q 1 ) d q 1 + a b U ( q 1 , q 2 , q 3 ) q 2 q 2 = y ( q 1 ) , q 3 = z ( q 1 ) y ( q 1 ) q 1 d q 1 +
z b U ( q 1 , q 2 , q 3 ) q 3 q 2 = y ( q 1 ) , q 3 = z ( q 1 ) z ( q 1 ) q 1 d q 1 =
a b d U ( q 1 , y ( q 1 ) , z ( q 1 ) ) d q 1 d q 1 = U ( b , d , f ) U ( a , c , e ) ,
where the standard gradient theorem is used.
As a result, we proved the general fractional gradient theorem for OCC with the line GF Gradient and the simple line in R 0 , + 3 . □

11.5. Regional and Surface GF Curl in OCC

In this section, we proposed definitions of two type GF Curl operators in orthogonal curvilinear coordinates (OCC) for R 0 , + 3 .

11.5.1. Regional GF Curl in OCC

Let us give the definition of Regional GF Curl operator in OCC for R 0 , + 3 .
Definition 60.
Let F ( q 1 , q 2 , q 3 ) be a vector field which satisfies conditions
F ^ k ( q 1 , q 2 , q 3 ) : = H k ( q 1 , q 2 , q 3 ) F k ( q 1 , q 2 , q 3 ) C 1 1 ( R + 3 )
for all k = 1 , 2 , 3 in OCC.
Then the regional general fractional curl in OCC for the region W = R 0 , + 3 is defined as
Curl W ( K ) F = e 1 Curl W ( K ) F 1 + e 2 Curl W ( K ) F 2 + e 3 Curl W ( K ) F 3 ,
where
Curl W ( K ) F 1 ( q 1 , q 2 , q 3 ) = 1 H 2 H 3 D ( K 2 ) q 2 , [ q 2 ] F ^ 3 ( q 1 , q 2 , q 3 ) D ( K 3 ) q 3 , [ q 3 ] F ^ 2 ( q 1 , q 2 , q 3 ) ,
Curl W ( K ) F 2 ( q 1 , q 2 , q 3 ) = 1 H 1 H 3 D ( K 3 ) q 3 , [ q 3 ] F ^ 1 ( q 1 , q 2 , q 3 ) D ( K 1 ) q 1 , [ q 1 ] F ^ 3 ( q 1 , q 2 , q 3 ) ,
Curl W ( K ) F 3 ( q 1 , q 2 , q 3 ) = 1 H 1 H 2 D ( K 1 ) q 1 , [ q 1 ] F ^ 2 ( q 1 , q 2 , q 3 ) D ( K 2 ) q 2 , [ q 2 ] F ^ 1 ( q 1 , q 2 , q 3 )
with q 1 , q 2 , q 3 0 .

11.5.2. Surface GF Integral in OCC

Using Definition 20 piecewise simple surface ( S P ( R 0 , + 3 ) ) in OCC, we propose definition of surface GF integral in OCC.
Let us define a set F S , H ( R 0 , + 3 ) . of vector fields that is used in definition of surface GF integral in OCC.
Definition 61.
Let S be a piecewise simple surface ( S P ( R 0 , + 3 ) ). Let the vector field
F : = e 1 F 1 ( q 1 , q 2 , q 3 ) + e 2 F 2 ( q 1 , q 2 , q 3 ) + e 3 F 3 ( q 1 , q 2 , q 3 )
on the surface S satisfy the conditions
F ˜ 1 ( x i ( q 2 , q 3 ) , q 2 , q 3 ) C 1 ( R + 2 ) ,
F ˜ 2 ( q 1 , y j ( q 1 , q 2 ) , q 3 ) C 1 ( R + 2 ) ,
F ˜ 3 ( q 1 , q 2 , z k ( q 2 , q 3 ) ) C 1 ( R + 2 ) ,
for all i = 1 , , n x , j = 1 , , n y , and k = 1 , , n z , where we use the notations
F ˜ 1 ( q 1 , q 2 , q 3 ) : = H 2 ( q 1 , q 2 , q 3 ) H 3 ( q 1 , q 2 , q 3 ) F 1 ( q 1 , q 2 , q 3 ) ,
F ˜ 2 ( q 1 , q 2 , q 3 ) : = H 1 ( q 1 , q 2 , q 3 ) H 3 ( q 1 , q 2 , q 3 ) F 2 ( q 1 , q 2 , q 3 ) ,
F ˜ 3 ( q 1 , q 2 , q 3 ) : = H 1 ( q 1 , q 2 , q 3 ) H 2 ( q 1 , q 2 , q 3 ) F 3 ( q 1 , q 2 , q 3 ) .
The set of such vector fields F on piecewise simple surface S will be denoted by F S , H ( R 0 , + 3 ) .
We can use the representation of the tilde F k in the form
F ˜ k ( q 1 , q 2 , q 3 ) : = H 1 ( q 1 , q 2 , q 3 ) H 2 ( q 1 , q 2 , q 3 ) H 3 ( q 1 , q 2 , q 3 ) F k ( q 1 , q 2 , q 3 ) H k ( q 1 , q 2 , q 3 ) .
Let us give a definition of surface GF integral in OCC.
Definition 62.
Let S be a piecewise simple surface ( S P ( R 0 , + 3 ) ) and a vector fields F on this surface S belongs to the set F S , H ( R 0 , + 3 ) .
Then, the surface general fractional vector integral in OCC (Surface GFI in OCC) of the second kind
I ^ S ( M ) = e 1 I ^ S 23 ( M ) [ q 2 , q 3 ] + e 2 I ^ S 13 ( M ) [ q 1 , q 3 ] + e 3 I ^ S 12 ( M ) [ q 1 , q 2 ] = e 1 i = 1 n 1 I ^ S i , 23 ( M ) [ q 2 , q 3 ] + e 1 j = 1 n 2 I ^ S j , 13 ( M ) [ q 1 , q 3 ] + e 2 k = 1 n 3 I ^ S k , 12 ( M ) [ q 1 , q 2 ]
for the vector fiels F F S , H ( R 0 , + 3 ) is defined by the equation
( I ^ S ( M ) , F ) : = i = 1 n 1 I ^ S i , 23 ( M ) [ q 2 , q 3 ] F 1 ( x i ( y , z ) , y , z ) +
j = 1 n 2 I ^ S j , 13 ( M ) [ q 1 , q 3 ] F 2 ( q 1 , y j ( q 1 , q 3 ) , q 3 ) + k = 1 n 3 I ^ S k , 12 ( M ) [ q 1 , q 2 ] F 3 ( q 1 , q 2 , z k ( q 1 , q 2 ) ) =
i = 1 n 1 I S i , 23 ( M ) [ q 2 , q 3 ] F ˜ 1 ( x i ( q 2 , q 3 ) , q 2 , q 3 ) + j = 1 n 2 I S j , 13 ( M ) [ q 1 , q 3 ] F ˜ 2 ( q 1 , y j ( q 1 , q 3 ) , q 3 ) +
k = 1 n 3 I S k , 12 ( M ) [ q 1 , q 2 ] F ˜ 3 ( q 1 , q 2 , z k ( q 1 , q 2 ) ) .
Here the areas
S 23 = i = 1 n 1 S i , 23 , S 13 = j = 1 n 2 S j , 13 , S 12 = k = 1 n 3 S k , 12
are the projections of the surface S onto the Q 2 Q 3 , Q 1 Q 3 , Q 1 Q 2 planes in OCC, where S i , 23 , S j , 13 and S k , 12 are simple areas in these planes (or simple along some axes).

11.5.3. Surface GF Curl in OCC

Let us now give a definition of the Surface GF Curl in OCC for R 0 , + 3 for piecewise simple surface. The piecewise simple surface is given in Definition 20.
Let us give a definition of the set F S , H 1 ( R 0 , + 3 ) of vector fields F that are used to define the surface GF Curl operator in OCC (compare with Definition 21 of the set F S 1 ( R 0 , + 3 ) that is used in definition of the surface GFI).
Definition 63.
Let S be a piecewise simple surface ( S P ( R 0 , + 3 ) ).
Let the vector field
F : = e 1 F 1 ( q 1 , q 2 , q 3 ) + e 2 F 2 ( q 1 , q 2 , q 3 ) + e 3 F 3 ( q 1 , q 2 , q 3 )
on the surface S satisfy the conditions
D q 3 1 F ^ 1 ( q 1 , y j ( q 1 , q 3 ) , q 3 ) : = F ^ 1 ( q 1 , q 2 , q 3 ) q 3 q 2 = y j ( q 1 , q 3 ) C 1 ( R + 2 ) ,
D q 2 1 F ^ 1 ( q 1 , q 2 , z k ( q 1 , q 2 ) ) : = F ^ 1 ( q 1 , q 2 , q 3 ) q 2 q 3 = z k ( q 1 , q 2 ) C 1 ( R + 2 ) ,
D q 1 1 F ^ 2 ( q 1 , q 2 , z k ( q 1 , q 2 ) ) : = F ^ 2 ( q 1 , q 2 , q 3 ) q 1 q 3 = z k ( q 1 , q 2 ) C 1 ( R + 2 ) ,
D q 3 1 F ^ 2 ( x i ( q 2 , q 3 ) , q 2 , q 3 ) : = F ^ 2 ( q 1 , q 2 , q 3 ) q 3 q 1 = x i ( q 2 , q 3 ) C 1 ( R + 2 ) ,
D q 1 1 F ^ 3 ( q 1 , y j ( q 1 , q 3 ) , q 3 ) : = F ^ 3 ( q 1 , q 2 , q 3 ) q 1 q 2 = y j ( q 1 , q 3 ) C 1 ( R + 2 ) ,
D q 2 1 F ^ 3 ( x i ( q 2 , q 3 ) , q 2 , q 3 ) : = F ^ 3 ( q 1 , q 2 , q 3 ) q 2 q 1 = x i ( q 2 , q 3 ) C 1 ( R + 2 )
for all i = 1 , , n 1 , j = 1 , , n 2 , and k = 1 , , n 3 .
The set of such vector fields F on piecewise simple surface S will be denoted by F S , H 1 ( R 0 , + 3 ) .
Then the surface general fractional curl is defined in the following form.
Definition 64.
Let S be a piecewise simple surface ( S P ( R 0 , + 3 ) ) and a vector fields F on this surface S belongs to the set F S 1 ( R 0 , + 3 ) .
Then the surface general fractional curl in OCC (Surface GF Curl in OCC) of the vector field F ( q 1 , q 2 , q 3 ) for the piecewise simple surface S is defined as
Curl S ( K ) F = e 1 Curl S ( K ) F 1 + e 2 Curl S ( K ) F 2 + e 3 Curl S ( K ) F 3 ,
where
Curl S ( K ) F 1 = i = 1 n 1 1 H 2 ( x i ( q 2 , q 3 ) , q 2 , q 3 ) H 3 ( x i ( q 2 , q 3 ) , q 2 , q 3 )
I ( K 2 ) q 2 , [ q 2 ] D q 2 1 F ^ 3 ( x i ( q 2 , q 3 ) , q 2 , q 3 ) I ( K 3 ) q 3 , [ q 3 ] D q 3 1 F ^ 2 ( x i ( q 2 , q 3 ) , q 2 , q 3 ) ,
Curl S ( K ) F 2 = j = 1 n 2 1 H 1 ( q 1 , y j ( q 1 , q 3 ) , q 3 ) H 3 ( q 1 , y j ( q 1 , q 3 ) , q 3 )
I ( K 3 ) q 3 , [ q 3 ] D q 3 1 F ^ 1 ( q 1 , y j ( q 1 , q 3 ) , q 3 ) I ( K 1 ) q 1 , [ q 1 ] D q 1 1 F ^ 3 ( q 1 , y j ( q 1 , q 3 ) , q 3 ) ,
Curl S ( K ) F 3 = k = 1 n 2 1 H 1 ( q 1 , q 2 , z k ( q 1 , q 2 ) ) H 2 ( q 1 , q 2 , z k ( q 1 , q 2 ) )
I ( K 1 ) q 1 , [ q 1 ] D q 1 1 F ^ 2 ( q 1 , q 2 , z k ( q 1 , q 2 ) ) I ( K 2 ) q 2 , [ q 2 ] D q 2 1 F ^ 1 ( q 1 , q 2 , z k ( q 1 , q 2 ) )
with q 1 , q 2 , q 3 0 .
Let us give the expression of the surface GFI of the vector field in the form of the Surface GF Curl in OCC.
Remark 25.
The surface GFI of the surface GF Curl in OCC has the form
( I ^ S ( M ) , Curl S ( K ) F ) : =
i = 1 n 1 I ^ S i , 23 ( M ) [ q 2 , q 3 ] Curl S ( K ) F 1 ( x i ( y , z ) , y , z ) +
j = 1 n 2 I ^ S j , 13 ( M ) [ q 1 , q 3 ] Curl S ( K ) F 2 ( q 1 , y j ( q 1 , q 3 ) , q 3 ) +
k = 1 n 3 I ^ S k , 12 ( M ) [ q 1 , q 2 ] Curl S ( K ) F 3 ( q 1 , q 2 , z k ( q 1 , q 2 ) ) =
i = 1 n 1 I S i , 23 ( M ) [ q 2 , q 3 ] H 2 ( x i ( q 2 , q 3 ) , q 2 , q 3 ) H 3 ( x i ( q 2 , q 3 ) , q 2 , q 3 ) Curl S ( K ) F 1 ( x i ( y , z ) , y , z ) +
j = 1 n 2 I S j , 13 ( M ) [ q 1 , q 3 ] H 1 ( q 1 , y j ( q 1 , q 3 ) , q 3 ) H 3 ( q 1 , y j ( q 1 , q 3 ) , q 3 ) Curl S ( K ) F 2 ( q 1 , y j ( q 1 , q 3 ) , q 3 ) +
k = 1 n 3 I S k , 12 ( M ) [ q 1 , q 2 ] H 1 ( q 1 , q 2 , z k ( q 1 , q 2 ) ) H 2 ( q 1 , q 2 , z k ( q 1 , q 2 ) ) Curl S ( K ) F 3 ( q 1 , q 2 , z k ( q 1 , q 2 ) ) =
i = 1 n 1 I S i , 23 ( M ) [ q 2 , q 3 ] ( I ( K 2 ) q 2 , [ q 2 ] D q 2 1 F ^ 3 ( x i ( q 2 , q 3 ) , q 2 , q 3 )
I ( K 3 ) q 3 , [ q 3 ] D q 3 1 F ^ 2 ( x i ( q 2 , q 3 ) , q 2 , q 3 ) ) +
j = 1 n 2 I S j , 13 ( M ) [ q 1 , q 3 ] ( I ( K 3 ) q 3 , [ q 3 ] D q 3 1 F ^ 1 ( q 1 , y j ( q 1 , q 3 ) , q 3 )
I ( K 1 ) q 1 , [ q 1 ] D q 1 1 F ^ 3 ( q 1 , y j ( q 1 , q 3 ) , q 3 ) ) +
k = 1 n 3 I S k , 12 ( M ) [ q 1 , q 2 ] ( I ( K 1 ) q 1 , [ q 1 ] D q 1 1 F ^ 2 ( q 1 , q 2 , z k ( q 1 , q 2 ) )
I ( K 2 ) q 2 , [ q 2 ] D q 2 1 F ^ 1 ( q 1 , q 2 , z k ( q 1 , q 2 ) ) ) .
For the case, F 2 = F 3 = 0 , we have
F = F 1 = e 1 P ,
and
( I ^ S ( M ) , Curl S ( K ) F ) : =
j = 1 n 2 I S j , 13 ( M ) [ q 1 , q 3 ] I ( K 3 ) q 3 , [ q 3 ] D q 3 1 P ^ ( q 1 , y j ( q 1 , q 3 ) , q 3 )
k = 1 n 3 I S k , 12 ( M ) [ q 1 , q 2 ] I ( K 2 ) q 2 , [ q 2 ] D q 2 1 P ^ ( q 1 , q 2 , z k ( q 1 , q 2 ) ) .
This case will be considered in the proof of the general fractional Stokes theorem for Surface GF Curl with simple surface, which is given in the next section.

11.6. General Fractional Stokes Theorem for Surface GF Curl in OCC

Let us prove the general fractional Stokes theorem for surface general fractional curl operator, where surface consists of simple surfaces or surfaces parallel to the coordinate planes.
Theorem 21.
(General Fractional Stokes Theorem for Surface GF Curl in OCC)
Let S R 0 , + 3 be a simple surface (or a piecewise simple surface S P ( R 0 , + 3 ) ), and the vector field F belongs to the set F S , H 1 ( R 0 , + 3 ) .
Then, the equation
I ^ S ( M ) , F = I ^ S ( M ) , Curl S F
holds, where Curl S ( K ) is the surface GF Curl that is defined by Definition 64.
Proof. 
Let us prove the theorem for the vector field F 1 = e 1 P ( q 1 , q 2 , q 3 ) . The general Stokes equations for the vector fields F 2 = e 2 F 2 ( q 1 , q 2 , q 3 ) and F 3 = e 3 F 3 ( q 1 , q 2 , q 3 ) are proved similarly. The proof for the vector field F = k = 0 3 e k F k is realized by the sum of the vector fields F k .
Let S R 0 , + 3 be a Q 3 -simple and Q 2 -simple surface that is described by the equation q 3 = z ( q 1 , q 2 ) for ( q 1 , q 2 ) S 12 R 0 , + 2 , and the equation q 2 = y ( q 1 , q 3 ) for ( q 1 , q 3 ) S 13 R 0 , + 2 .
The values of the function P ( q 1 , q 2 , q 3 ) on the line L R 0 , + 3 are equal to the values of the function P ( q 1 , q 2 , z ( q 1 , q 2 ) ) on the line L 12 R 0 , + 2 , which is a projection of the line L onto the Q 1 Q 2 -plane,
I S ( M ) , F = I ^ S ( M 1 ) [ q 1 ] P ( q 1 , q 2 , q 3 ) = I ^ S 12 ( M 1 ) [ q 1 ] P ( q 1 , q 2 , z ( q 1 , q 2 ) ) ,
where L = S and L 12 = S 12 .
Let us assume that L 12 consists of two Q 2 -simple lines L 12 a and L 12 b that are described by equations q 2 = y a ( q 1 ) 0 and q 2 = y b ( q 1 ) , where y b ( q 1 ) y a ( q 1 ) 0 . Then
I ^ L 12 ( M 1 ) [ q 1 ] P ( q 1 , q 2 , z ( q 1 , q 2 ) ) = I ^ L 12 a ( M 1 ) [ q 1 ] P ( q 1 , q 2 , z ( q 1 , q 2 ) ) + I ^ L 12 b ( M 1 ) [ q 1 ] P ( q 1 , q 2 , z ( q 1 , q 2 ) ) = I [ a , b ] ( M 1 ) [ q 1 ] P ^ ( q 1 , y 1 ( q 1 ) , z ( q 1 , y 1 ( q 1 ) ) ) I [ a , b ] ( M 1 ) [ q 1 ] P ^ ( q 1 , y 2 ( q 1 ) , z ( q 1 , y 2 ( q 1 ) ) ) = I [ a , b ] ( M 1 ) [ q 1 ] P ^ ( q 1 , y 1 ( q 1 ) , z ( q 1 , y 1 ( q 1 ) ) ) P ^ ( q 1 , y 2 ( q 1 ) , z ( q 1 , y 2 ( q 1 ) ) ) ,
where
P ^ ( q 1 , q 2 , q 3 ) : = H 1 ( q 1 , q 2 , q 3 ) P ( q 1 , q 2 , q 3 ) .
Using the fundamental theorem of general fractional calculus, expression (210) can be written as
I [ a , b ] ( M 1 ) [ q 1 ] P ^ ( q 1 , y 1 ( q 1 ) , z ( q 1 , y 1 ( q 1 ) ) ) P ^ ( q 1 , y 2 ( q 1 ) , z ( q 1 , y 2 ( q 1 ) ) ) = I [ a , b ] ( M 1 ) [ q 1 ] I [ y 1 ( q 1 ) , y 2 ( q 1 ) ] ( M 2 ) [ q 2 ] D ( K 2 ) q 2 , [ q 2 ] P ^ ( q 1 , q 2 , z ( q 1 , q 2 ) ) .
Then using the definition of the general fractional derivative and the standard chain rule for the first-order derivative, we get
I [ a , b ] ( M 1 ) [ q 1 ] I [ y 1 ( q 1 ) , y 2 ( q 1 ) ] ( M 2 ) [ q 2 ] D ( K 2 ) q 2 , [ q 2 ] P ^ ( q 1 , q 2 , z ( q 1 , q 2 ) ) =
I [ a , b ] ( M 1 ) [ q 1 ] I [ y 1 ( q 1 ) , y 2 ( q 1 ) ] ( M 2 ) [ q 2 ] I ( K 2 ) q 2 , [ q 2 ] P ^ ( q 1 , q 2 , z ( q 1 , q 2 ) ) q 2 =
I [ a , b ] ( M 1 ) [ q 1 ] I [ y 1 ( q 1 ) , y 2 ( q 1 ) ] ( M 2 ) [ q 2 ] I ( K 2 ) q 2 , [ q 2 ] ( P ^ ( q 1 , q 2 , q 3 ) q 2 q 3 = z ( q 1 , q 2 ) +
P ^ ( q 1 , q 2 , q 3 ) q 3 q 3 = z ( q 1 , q 2 ) z ( q 1 , q 2 ) q 2 ) =
I [ a , b ] ( M 1 ) [ q 1 ] I [ y 1 ( q 1 ) , y 2 ( q 1 ) ] ( M 2 ) [ q 2 ] I ( K 2 ) q 2 , [ q 2 ] P ^ ( q 1 , q 2 , q 3 ) q 2 q 3 = z ( q 1 , q 2 )
I [ a , b ] ( M 1 ) [ q 1 ] I [ y 1 ( q 1 ) , y 2 ( q 1 ) ] ( M 2 ) [ q 2 ] I ( K 2 ) q 2 , [ q 2 ] P ^ ( q 1 , q 2 , q 3 ) q 3 q 3 = z ( q 1 , q 2 ) z ( q 1 , q 2 ) q 2 .
Let us use the property
I [ y 1 , y 2 ] ( M 2 ) [ q 2 ] I ( K 2 ) q 2 , [ q 2 ] U ( q 2 ) = ( M 2 ( K 2 U ) ) ( y 2 ) ( M 2 ( K 2 U ) ) ( y 1 ) =
( ( M 2 K 2 ) U ) ( y 2 ) ( ( M 2 K 2 ) U ) ( y 1 ) =
( { 1 } U ) ( y 2 ) ( { 1 } U ) ( y 1 ) = y 1 y 2 U ( q 2 ) d q 2
for the expression of term (213) in the form
I [ y 1 ( q 1 ) , y 2 ( q 1 ) ] ( M 2 ) [ q 2 ] I ( K 2 ) q 2 , [ q 2 ] P ^ ( q 1 , q 2 , q 3 ) q 3 q 3 = z ( q 1 , q 2 ) z ( q 1 , q 2 ) q 2 = y 1 ( q 1 ) y 2 ( q 1 ) P ^ ( q 1 , q 2 , q 3 ) q 3 q 3 = z ( q 1 , q 2 ) z ( q 1 , q 2 ) q 2 d q 2 .
Then we can use the equation, which is used in the standard Stokes theorem, ( { 1 } U ) ( q 3 ) = ( M 3 K 3 U ) ( q 3 ) , and the fundamental theorem of GFC in the form
y 1 ( q 1 ) y 2 ( q 1 ) P ^ ( q 1 , q 2 , q 3 ) q 3 q 3 = z ( q 1 , q 2 ) z ( q 1 , q 2 ) q 2 d q 2 = z 1 ( q 1 ) z 2 ( q 1 ) P ^ ( q 1 , q 2 , q 3 ) q 3 q 2 = y ( q 1 , q 3 ) d q 3 = I [ z 1 ( q 1 ) , z 2 ( q 1 ) ] ( M 3 ) [ q 3 ] I ( K 3 ) q 3 , [ q 3 ] P ^ ( q 1 , q 2 , q 3 ) q 3 q 2 = y ( q 1 , q 3 ) .
Using (215), expression (213) can be written in the form
I [ a , b ] ( M 1 ) [ q 1 ] I [ y 1 ( q 1 ) , y 2 ( q 1 ) ] ( M 2 ) [ q 2 ] I ( K 2 ) q 2 , [ q 2 ] P ^ ( q 1 , q 2 , q 3 ) q 3 q 3 = z ( q 1 , q 2 ) z ( q 1 , q 2 ) q 2 = I [ a , b ] ( M 1 ) [ q 1 ] I [ z 1 ( q 1 ) , z 2 ( q 1 ) ] ( M 3 ) [ q 3 ] I ( K 3 ) q 3 , [ q 3 ] P ^ ( q 1 , q 2 , q 3 ) q 3 q 2 = y ( q 1 , q 3 ) .
Therefore, for two terms (212) and (213), we obtain the equations
I [ a , b ] ( M 1 ) [ q 1 ] I [ y 1 ( q 1 ) , y 2 ( q 1 ) ] ( M 2 ) [ q 2 ] I ( K 2 ) q 2 , [ q 2 ] P ^ ( q 1 , q 2 , q 3 ) q 2 q 3 = z ( q 1 , q 2 ) = I [ a , b ] ( M 1 ) [ q 1 ] I [ y 1 ( q 1 ) , y 2 ( q 1 ) ] ( M 2 ) [ q 2 ] 1 H 1 ( q 1 , q 2 , z ( q 1 , q 2 ) ) H 2 ( q 1 , q 2 , z ( q 1 , q 2 ) ) · I ( K 2 ) q 2 , [ q 2 ] P ^ ( q 1 , q 2 , q 3 ) q 2 q 3 = z ( q 1 , q 2 ) = I ^ S 12 ( M ) [ q 1 , q 2 ] Curl S ( K ) F 1 3 ,
and
I [ a , b ] ( M 1 ) [ q 1 ] I [ z 1 ( q 1 ) , z 2 ( q 1 ) ] ( M 3 ) [ q 3 ] I ( K 3 ) q 3 , [ q 3 ] P ^ ( q 1 , q 2 , q 3 ) q 3 q 2 = y ( q 1 , q 3 ) = I [ a , b ] ( M 1 ) [ q 1 ] I [ z 1 ( q 1 ) , z 2 ( q 1 ) ] ( M 3 ) [ q 3 ] 1 H 1 ( q 1 , q 2 , z ( q 1 , q 2 ) ) H 2 ( q 1 , q 2 , z ( q 1 , q 2 ) ) · I ( K 3 ) q 3 , [ q 3 ] P ^ ( q 1 , q 2 , q 3 ) q 3 q 2 = y ( q 1 , q 3 ) = I ^ S 13 ( M ) [ q 1 , q 3 ] Curl S ( K ) F 1 2 ,
for the vector field F 1 = e 1 P ( q 1 , q 2 , q 3 ) .
As a result, we proved the general Stokes theorem for the surface FG Curl and the vector field F 1 = e 1 P ( q 1 , q 2 , q 3 ) .
I ^ S ( M 1 ) , F 1 = I ^ S ( M ) , Curl S ( K ) F 1 .
Equations for the remaining components of the vector field F 2 = e 2 F 2 ( q 1 , q 2 , q 3 ) and F 3 = e 3 F 3 ( q 1 , q 2 , q 3 ) are proved similarly. □

11.7. General Fractional Gauss Theorem in OCC

11.7.1. Definition of Triple GFI by Iterated GFI in OCC

Let us use the concept Definition 38 of Z-simple region for OCC in R 0 , + 3 .
Definition 65.
Let W R 0 , + 3 be Q 3 -simple domain in OCC that is bounded above and below by smooth surfaces S 12 b , S 12 a described by the equations
q 3 = z 1 ( q 1 , q 2 ) 0 , q 3 = z 2 ( q 1 , q 2 ) 0 ,
where the functions are continuous in the closed domain W 12 that is a projection of the region W onto the Q 1 Q 2 -plane, and z 3 ( q 1 , q 2 ) z 1 ( q 1 , q 2 ) for all ( q 1 , q 2 ) W 12 .
Let scalar field U ( q 1 , q 2 , q 3 ) be satisfy the condition
I [ z 1 ( q 1 , q 2 ) , z 2 ( q 1 , q 2 ) ] ( M 3 ) [ q 3 ] U ˜ ( q 1 , q 2 , q 3 ) C 1 ( R + 2 ) ,
where
U ˜ ( q 1 , q 2 , q 3 ) : = H 1 ( q 1 , q 2 , q 3 ) H 2 ( q 1 , q 2 , q 3 ) H 3 ( q 1 , q 2 , q 3 ) U ( q 1 , q 2 , q 3 ) .
Then the triple general fractional integral in OCC is defined by the equation
I ^ W ( M ) [ q 1 , q 2 , q 3 ] U ( q 1 , q 2 , q 3 ) : = I W 12 ( M ) [ q 1 , q 2 ] I [ z 1 ( q 1 , q 2 ) , z 2 ( q 1 , q 2 ) ] ( M 3 ) [ q 3 ] U ˜ ( q 1 , q 2 , q 3 ) ,
where I W 12 ( M ) [ q 1 , q 2 ] is double GFI.
Definition 66.
Let W R 0 , + 3 be Q 3 -simple domain that is that is bounded above and below by smooth surfaces S 12 b , S 12 a described by Equation (219). Let W 12 R 0 , + 2 be projection of W on the Q 1 Q 2 -plane such that W 12 is Q 2 -simple region in Q 1 Q 2 -plane that is bounded by the lines q 2 = y 1 ( q 1 ) and q 2 = y 2 ( q 1 ) , where q 2 = y 1 ( q 1 ) and q 2 = y 2 ( q 1 ) are continuous functions on the interval [ a , b ] , b > a 0 and y 2 ( q 1 ) y 1 ( q 1 ) > 0 for all q 1 [ a , b ] .
Let scalar field U ( q 1 , q 2 , q 3 ) satisfy condition (220) and
I [ y 1 ( q 1 ) , y 2 ( q 1 ) ] ( M 2 ) [ q 2 ] I [ z 1 ( q 1 , q 2 ) , z 2 ( q 1 , q 2 ) ] ( M 3 ) [ q 3 ] U ˜ ( q 1 , q 2 , q 3 ) C 1 ( 0 , ) .
Then the triple general fractional integral (triple GFI) is defined in the form
I ^ W ( M ) [ q 1 , q 2 , q 3 ] U ( q 1 , q 2 , q 3 ) : = I [ a , b ] ( M 1 ) [ q 1 ] I [ y 1 ( q 1 ) , y 2 ( q 1 ) ] ( M 2 ) [ q 2 ] I [ z 1 ( q 1 , q 2 ) , z 2 ( q 1 , q 2 ) ] ( M 3 ) [ q 3 ] U ˜ ( q 1 , q 2 , q 3 ) .
A volume general fractional integral (volume GFI) of a scalar field is a triple general fractional integral for the Z-simple region W R 0 , + 3 .
Example 19.
Using the parallelepiped region
W : = { 0 q 1 a 1 , 0 q 2 a 2 , 0 q 3 a 3 } ,
the volume general fractional integral can be written as
I W ( M ) [ q 1 , q 2 , q 3 ] f ( q 1 , q 2 , q 3 ) = 0 a 1 d q 1 0 a 2 d q 2 0 a 3 d q 3 M 1 ( a 1 q 1 ) M 2 ( a 2 q 2 ) M 3 ( a 3 q 3 ) f ˜ ( q 1 , q 2 , q 3 ) .
For the power-law kernels, the general fractional flux (223) has the form
I ^ W ( M ) [ q 1 , q 2 , q 3 ] f ( q 1 , q 2 , q 3 ) = 0 a 1 d q 1 0 a 2 d q 2 0 a 3 d q 3 ( a 1 q 1 ) α 1 1 Γ ( α 1 ) ( a 2 q 2 ) α 2 1 Γ ( α 2 ) ( a 3 q 3 ) α 3 1 Γ ( α 3 ) f ˜ ( q 1 , q 2 , q 3 ) .
For α 1 = α 2 = α 3 = 1 , Equation (224) is given as
I ^ W ( M ) [ q 1 , q 2 , q 3 ] f ( q 1 , q 2 , q 3 ) = 0 a 1 d q 1 0 a 2 d q 2 0 a 3 d q 3 f ˜ ( q 1 , q 2 , q 3 ) ,
which is the standard volume integral for the function f ( q 1 , q 2 , q 3 ) .

11.7.2. General Fractional Divergence in OCC

In this section, we give the definition of general fractional divergence in OCC for R 0 , + 3 .
Let us define sets of vector fields that will be used in the definition of the regional general fractional divergence.
Definition 67.
Let F ( q 1 , q 2 , q 3 ) be a vector field that satisfies the conditions
D ( K 1 ) q 1 , [ q 1 ] F ˜ 1 ( q 1 , q 2 , q 3 ) C 1 ( R + 3 ) ,
D ( K 2 ) q 2 , [ q 2 ] F ˜ 2 ( q 1 , q 2 , q 3 ) C 1 ( R + 3 ) ,
D ( K 3 ) q 3 , [ q 3 ] F ˜ 3 ( q 1 , q 2 , q 3 ) C 1 ( R + 3 ) ,
where
F ˜ k ( q 1 , q 2 , q 3 ) = H 1 ( q 1 , q 2 , q 3 ) H 2 ( q 1 , q 2 , q 3 ) H 3 ( q 1 , q 2 , q 3 ) F k ( q 1 , q 2 , q 3 ) H k ( q 1 , q 2 , q 3 ) .
Then the set of such vector fields will be denoted as F 1 , D i v , H 1 ( R + 3 ) .
We can also consider vector fields F ( q 1 , q 2 , q 3 ) that satisfy the condition F ^ k C 1 1 ( R + 3 ) .
Definition 68.
Let F ( q 1 , q 2 , q 3 ) be a vector field that satisfies the conditions
F ˜ k ( q 1 , q 2 , q 3 ) C 1 , H 1 ( R + 3 )
for all k = 1 , 2 , 3 .
Then the set of such vector fields will be denoted as C 1 , H 1 ( R + 3 ) .
Let us define the general fractional divergence in OCC.
Definition 69.
Let F ( q 1 , q 2 , q 3 ) be a vector field that belongs to the set F 1 , D i v , H 1 ( R + 3 ) or C 1 , H 1 ( R + 3 ) .
Then the regional general fractional divergence Div W ( K ) in OCC for the region W = R 0 , + 3 is defined as
Div W ( K ) F = 1 H 1 H 2 H 3 D ( K 1 ) q 1 , [ q 1 ] F ˜ 1 ( q 1 , q 2 , q 3 ) + D ( K 2 ) q 2 , [ q 2 ] F ˜ 2 ( q 1 , q 2 , q 3 ) + D ( K 3 ) q 3 , [ q 3 ] F ˜ 3 ( q 1 , q 2 , q 3 ) ,
where H k = H k ( q 1 , q 2 , q 3 ) .
Divergence (226) is the regional general fractional differential operator in OCC. We can also define the line and surface general fractional divergence in OCC.
Remark 26.
The equation, which defines the regional general fractional divergence in OCC, can be written in compact form. If the vector field
F ( q 1 , q 2 , q 3 ) = k = 1 3 e k F k ( q 1 , q 2 , q 3 )
belongs to the function space C 1 1 ( R + 3 ) , then the regional general fractional divergence in OCC for the region W = R + 3 can be defined as
( Div W ( K ) F ˜ ) ( q 1 , q 2 , q 3 ) = D ^ W ( K ) , F = 1 H 1 H 2 H 3 k = 1 3 D ( K k ) q k , [ q k ] F ˜ k ,
where
D ^ W ( K ) : = e 1 D ( K 1 ) q 1 , [ q 1 ] + e 2 D ( K 2 ) q 2 , [ q 2 ] + e 3 D ( K 3 ) q 3 , [ q 3 ] .
Remark 27.
The general fractional divergence in OCC can be defined not only for W = R 0 , + 3 , but also for regions W R 0 , + 3 , surfaces S R 0 , + 3 and line L R 0 , + 3 . Note that the general fractional divergence in OCC for regions W R 0 , + 3 is used in the general fractional Gausss theorem for OCC to be given and proved in the following sections.

11.8. General Fractional Gauss Theorem for OCC

Let us define a set of vector field, for which general fractional Gauss theorem will be formulated.
Definition 70.
Let W in R 0 , + 3 be Q 3 -simple region such that W is a piecewise Q 2 -simple and Q 1 -simple region
W = k = 1 n W Q 1 , k = j = 1 m W Q 2 , j ,
where W Q 1 , k are the Q ! X -simple regions that is described by q 1 = x k , 1 ( q 2 , q 3 ) , q 1 = x k , 2 ( q 2 , q 3 ) for ( q 2 , q 3 ) D 12 , and W Q 2 , j are the Q 2 -simple regions that are described by the functions q 2 = y j , 1 ( q 1 , q 3 ) , q 2 = y j , 2 ( q 1 , q 3 ) for ( q 1 , q 3 ) D 13 .
Let vector field F ( q 1 , q 2 , q 3 ) satisfy the conditions
F ˜ 1 ( x k , 2 ( q 2 , q 3 ) , q 2 , q 3 ) , F ˜ 1 ( x k , 1 ( q 2 , q 3 ) , q 2 , q 3 ) C 1 ( R + 2 ) ,
F ˜ 2 ( q 1 , y j , 2 ( q 1 , q 3 ) , q 3 ) , F ˜ 2 ( q 1 , y j , 1 ( q 1 , q 3 ) , q 3 ) C 1 ( R + 2 ) ,
F ˜ 3 ( q 1 , q 2 , z 2 ( q 1 , q 2 ) ) , F ˜ 3 ( q 1 , q 2 , z 1 ( q 1 , q 2 ) ) C 1 ( R + 2 ) ,
where
F ˜ k ( q 1 , q 2 , q 3 ) : = H 1 ( q 1 , q 2 , q 3 ) H 2 ( q 1 , q 2 , q 3 ) H 3 ( q 1 , q 2 , q 3 ) F k ( q 1 , q 2 , q 3 ) H k ( q 1 , q 2 , q 3 ) .
Then the set of such vector fields F will be denoted as F ˜ 1 , W ( R + 2 ) .
Definition 71.
Let vector field F ( q 1 , q 2 , q 3 ) satisfy the conditions
f 1 ( q 1 , q 2 , q 3 ) : = D ( K 1 ) q 1 , [ q 1 ] F ˜ 1 ( q 1 , q 2 , q 3 ) C 1 ( R + 3 ) ,
f 2 ( q 1 , q 2 , q 3 ) : = D ( K 2 ) q 2 , [ q 2 ] F ˜ 2 ( q 1 , q 2 , q 3 ) C 1 ( R + 3 ) ,
f 3 ( q 1 , q 2 , q 3 ) : = D ( K 3 ) q 3 , [ q 3 ] F ˜ 3 ( q 1 , q 2 , q 3 ) C 1 ( R + 3 ) .
Then the set of such vector fields F will be denoted as F ˜ 1 1 ( R + 3 ) .
Theorem 22.
(General fractional Gauss theorem for OCC)
Let W in R 0 , + 3 be Q 3 -simple region such that W is a piecewise Q 1 -simple and Q 2 -simple region. Let W is bounded above and below by smooth surfaces S 12 b , S 12 a , which is described by Equation (219), and a lateral surface S 3 , whose generatrices are parallel to the Q 3 -axis.
Let the vector field F ( q 1 , q 2 , q 3 ) belongs to the sets F ˜ 1 ( W ) and F ˜ 1 1 ( R + 3 ) .
Then the general fractional Gauss equation has the from
I ^ W ( M ) [ q 1 , q 2 , q 3 ] Div W ( K ) F = I ^ W , F
that can be written as
I W ( M ) [ q 1 , q 2 , q 3 ] ( D ( K 1 ) q 1 , [ q 1 ] F ˜ 1 ( q 1 , q 2 , q 3 ) +
D ( K 2 ) q 2 , [ q 2 ] F ˜ 2 ( q 1 , q 2 , q 3 ) + D ( K 3 ) q 3 , [ q 3 ] F ˜ 3 ( q 1 , q 2 , q 3 ) ) =
I ^ S ( M ) [ q 2 , q 3 ] F 1 ( q 1 , q 2 , q 3 ) + I ^ S ( M ) [ q 1 , q 2 ] F 2 ( q 1 , q 2 , q 3 ) + I ^ S ( M ) [ q 1 , q 2 ] F 3 ( q 1 , q 2 , q 3 ) .
Proof. 
Let us consider the triple GFI
I ^ W ( M ) [ q 1 , q 2 , q 3 ] U = I W ( M ) [ q 1 , q 2 , q 3 ] U ˜ = I W 12 ( M ) [ q 1 , q 2 ] I [ z 1 ( q 1 , q 2 ) , z 2 ( q 1 , q 2 ) ] ( M 3 ) [ q 3 ] U ˜ ( q 1 , q 2 , q 3 )
for the scalar field
U ( q 1 , q 2 , q 3 ) = 1 H 1 ( q 1 , q 2 , q 3 ) H 2 ( q 1 , q 2 , q 3 ) H 3 ( q 1 , q 2 , q 3 ) D ( K 3 ) q 3 , [ q 3 ] F ˜ 3 ( q 1 , q 2 , q 3 ) C 1 ( 0 , ) ,
where F 3 ( q 1 , q 2 , q 3 ) C 1 1 ( 0 , ) for each point ( q 1 , q 2 ) D 12 .
Using Equations (221), (227) and (228), we get
I W ( M ) [ q 1 , q 2 , q 3 ] U ˜ = I W 12 ( M ) [ q 1 , q 2 ] D ( K 3 ) q 3 , [ q 3 ] F ˜ 3 ( q 1 , q 2 , q 3 ) ,
since H 1 H 2 H 3 ( H 1 H 2 H 3 ) 1 = 1 .
Using the second fundamental theorem of GFC, we get
I [ z 1 ( q 1 , q 2 ) , z 2 ( q 1 , q 2 ) ] ( M 3 ) [ q 3 ] D ( K 3 ) q 3 , [ q 3 ] F ˜ 3 ( q 1 , q 2 , q 3 ) = I ( M 3 ) z 2 ( q 1 , q 2 ) [ q 3 ] D ( K 3 ) q 3 , [ q 3 ] F 3 ( q 1 , q 2 , q 3 ) I ( M 3 ) z 1 ( q 1 , q 2 ) [ q 3 ] D ( K 3 ) q 3 , [ q 3 ] F ˜ 3 ( q 1 , q 2 , q 3 ) = F ˜ 3 ( q 1 , q 2 , z 2 ( q 1 , q 2 ) ) F ˜ 3 ( q 1 , q 2 , 0 ) F ˜ 3 ( q 1 , q 2 , z 2 ( q 1 , q 2 ) ) F ˜ 3 ( q 1 , q 2 , 0 ) ) = F ˜ 3 ( q 1 , q 2 , z 2 ( q 1 , q 2 ) ) F ˜ 3 ( q 1 , q 2 , z 1 ( q 1 , q 2 ) ) .
Substituting expression (230) into Equation (229) and assuming that the condition
F ˜ 3 ( q 1 , q 2 , z 1 ( q 1 , q 2 ) ) , F ˜ 3 ( q 1 , q 2 , z 2 ( q 1 , q 2 ) ) C 1 ( R + 2 ) ,
is satisfied, we obtain the equation
I W ( M ) [ q 1 , q 2 , q 3 ] D ( K 3 ) q 3 , [ q 3 ] F ˜ 3 ( q 1 , q 2 , q 3 ) =
I W 12 ( M ) [ q 1 , q 2 ] F ˜ 3 ( q 1 , q 2 , z 2 ( q 1 , q 2 ) ) F ˜ 3 ( q 1 , q 2 , z 1 ( q 1 , q 2 ) ) .
Therefore, the triple GFI can be represented through the surface GFIs in the form
I W ( M ) [ q 1 , q 2 , q 3 ] D ( K 3 ) q 3 , [ q 3 ] F ˜ 3 ( q 1 , q 2 , q 3 ) = I S 12 b ( M ) [ q 1 , q 2 ] F ˜ 3 ( q 1 , q 2 , q 3 ) + I S 12 a ( M ) [ q 1 , q 2 ] F ˜ 3 ( q 1 , q 2 , q 3 ) ,
where the surface GFIs are represented by the definition in the form
I S 12 b ( M ) [ q 1 , q 2 ] F ˜ 3 ( q 1 , q 2 , q 3 ) = I W 12 ( M ) [ q 1 , q 2 ] F ˜ 3 ( q 1 , q 2 , z 2 ( q 1 , q 2 ) ) ,
I S 12 a ( M ) [ q 1 , q 2 ] F ˜ 3 ( q 1 , q 2 , q 3 ) = I W 12 ( M ) [ q 1 , q 2 ] F ˜ 3 ( q 1 , q 2 , z 1 ( q 1 , q 2 ) ) .
Then, taking into account that the surface GFI over the surface S z is equal to zero
I S 3 ( M ) [ q 1 , q 2 ] F ˜ 3 ( q 1 , q 2 , q 3 ) = 0 ,
we obtain
I W ( M ) [ q 1 , q 2 , q 3 ] D ( K 3 ) q 3 , [ z ] F ˜ 3 ( q 1 , q 2 , q 3 ) = I S ( M ) [ q 1 , q 2 ] F ˜ 3 ( q 1 , q 2 , q 3 ) = I ^ S ( M ) [ q 1 , q 2 ] F 3 ( q 1 , q 2 , q 3 ) ,
where W = S = S 12 a S 12 b S 3 is closed surface that contains the region W inside the surface W .
Further, similarly with the proof of general fractional Gauss theorem for the Cartesian coordinate system and and by analogy with the transformation described above, we obtain the equation of the theorem to be proved. □
Example 20.
As an example, let us give the general fractional Gauss theorem in cylindrical coordinates.
Let us consider the simple domain in the form
W : = { ( r , φ , z ) : 0 r 1 r r 2 , 0 φ 1 φ φ 2 < 2 π , 0 z 1 z z 2 } .
The vertices of the domain are the following points
A 1 ( r 1 , φ 1 , z 1 ) , A 2 ( r 2 , φ 1 , z 1 ) , B 1 ( r 1 , φ 2 , z 1 ) , B 2 ( r 2 , φ 2 , z 1 ) ,
A 1 ( r 1 , φ 1 , z 2 ) , A 2 ( r 2 , φ 1 , z 2 ) , B 1 ( r 1 , φ 2 , z 2 ) , B 2 ( r 2 , φ 2 , z 2 ) .
Let F F 1 1 ( W ) , then
I W ( M ) ^ [ r , φ , z ] Div W ( K ) F = I W ( M ) ^ [ r , φ , z ] , F ( r , φ , z ) .
For simple domain (231) the cylindrical general fractional Gauss equation has the form
I [ r 1 , r 2 ] ( M r ) ^ [ r ] I [ φ 1 , φ 2 ] ( M φ ) ^ [ φ ] I [ z 1 , z 2 ] ( M z ) ^ [ z ] ( 1 r D ( K r ) r , [ r ] r F r ( r , φ , z ) +
1 r D ( K φ ) φ , [ φ ] F φ ( r , φ , z ) + D ( K z ) z , [ z ] F z ( r , φ , z ) ) =
I W φ , z ( M ) ^ [ φ , z ] F r ( r , φ , z ) + I W r , z ( M ) ^ [ r , z ] F φ ( r , φ , z ) + I W r , φ ( M ) ^ [ r , φ ] F z ( r , φ , z ) .

12. Conclusions

The general fractional vector calculus (General FVC) is proposed as a generalization of the fractional vector calculus suggested in [13,31]. The formulation of the General FVC is based on the results in the calculus of general fractional integrals and derivatives that is proposed in Luchko’s work [41].
The formulation of General FVC is self-consistent form, i.e., definitions of fractional generalizations of differential and integral vector operators are consistent with each other, and generalizations of fundamental theorems of vector calculus were proved. In this paper, the definitions of general fractional integral vector operators: the general fractional circulation, general fractional flux and general fractional volume integral are proposed. Definitions of general fractional differential vector operators, including the regional and line general fractional gradients, the regional and surface general fractional curl operators, the general fractional divergence, are suggested. Fundamental theorems of General FVC, which are general analogs of the standard gradient, Green’s, Stokes’ and Gauss’s theorems, are proved for simple and complex regions. Let us emphasize that the fractional vector analogs of fundamental theorems (such as the gradient, Stock’s and Gauss theorems) are not fulfilled for all type (regional, surface and line) of the general fractional vector operators (the gradient, curl and divergence). In the general case, the general fractional (GF) gradient theorem should be considered for the line GF Gradient, the FG Stock’s theorem should be considered for the surface GF Curl operator, and the GF Green theorem should be considered for the regional GF Divergence. This is due to violation of the chain rule for general fractional derivatives. The General FVC for orthogonal curvilinear coordinates, which includes general fractional vector operators for the spherical and cylindrical coordinates, is described.
The proposed General FVC can be used as a mathematical tool in general fractional dynamics (GFDynamics) [44,45], in which non-locality in space is taken into account in general form.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declare no conflict of interest.

References

  1. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives Theory and Applications; Gordon and Breach: New York, NY, USA, 1993; 1006p, ISBN 9782881248641. [Google Scholar]
  2. Kiryakova, V. Generalized Fractional Calculus and Applications; Longman and J. Wiley: New York, NY, USA, 1994; 360p, ISBN 9780582219779. [Google Scholar]
  3. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1998; 340p. [Google Scholar]
  4. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006; 540p. [Google Scholar]
  5. Diethelm, K. The Analysis of Fractional Differential Equations: An Application-Oriented Exposition Using Differential Operators of Caputo Type; Springer: Berlin, Germany, 2010; 247p. [Google Scholar] [CrossRef] [Green Version]
  6. Kochubei, A.; Luchko, Y. (Eds.) Handbook of Fractional Calculus with Applications. Volume 1. Basic Theory; Walter de Gruyter GmbH: Berlin, Germany; Boston, MA, USA, 2019; 481p, ISBN 978-3-11-057081-6. [Google Scholar] [CrossRef]
  7. Kochubei, A.; Luchko, Y. (Eds.) Handbook of Fractional Calculus with Applications. Volume 2. Fractional Differential Equations; Walter de Gruyter GmbH: Berlin, Germany; Boston, MA, USA, 2019; 519p, ISBN 978-3-11-057082-3. [Google Scholar] [CrossRef] [Green Version]
  8. Tarasov, V.E. No nonlocality. No fractional derivative. Commun. Nonlinear Sci. Numer. Simul. 2018, 62, 157–163. [Google Scholar] [CrossRef] [Green Version]
  9. Tarasov, V.E. On chain rule for fractional derivatives. Commun. Nonlinear Sci. Numer. Simul. 2016, 30, 1–4. [Google Scholar] [CrossRef]
  10. 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]
  11. Tarasov, V.E. (Ed.) Handbook of Fractional Calculus with Applications. Volume 4. Application in Physics. Part A; Walter de Gruyter GmbH: Berlin, Germany; Boston, MA, USA, 2019; 306p, ISBN 978-3-11-057088-5. [Google Scholar] [CrossRef]
  12. Tarasov, V.E. (Ed.) Handbook of Fractional Calculus with Applications. Volume 5. Application in Physics. Part B; Walter de Gruyter GmbH: Berlin, Germany; Boston, MA, USA, 2019; 319p, ISBN 978-3-11-057089-2. [Google Scholar] [CrossRef] [Green Version]
  13. Tarasov, V.E. Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media; Springer: Berlin/Heidelberg, Germany, 2010; 516p, ISBN 978-3-642-14003-7. [Google Scholar] [CrossRef] [Green Version]
  14. Mainardi, F. Fractional Calculus and Waves Linear Viscoelasticity; Imperial College Press: Lodon, UK, 2010. [Google Scholar] [CrossRef]
  15. Klafter, J.; Lim, S.C.; Metzler, R. (Eds.) Fractional Dynamics—Recent Advances; World Scientific: Singapore, 2012. [Google Scholar] [CrossRef]
  16. Tarasov, V.E. (Ed.) Mathematical Economics: Application of Fractional Calculus; MDPI: Basel, Switzerland; Beijing, China, 2020; 278p, ISBN 978-3-03936-119-9. [Google Scholar] [CrossRef]
  17. Tarasov, V.E.; Tarasova, V.V. Economic Dynamics with Memory: Fractional Calculus Approach; De Gruyter: Berlin, Germany; Boston, MA, USA, 2021; 576p, ISBN 978-3-11-062460-1. [Google Scholar] [CrossRef]
  18. Ionescu, C.; Lopes, A.; Copot, D.; Tenreiro Machado, J.; Bates, J. The role of fractional calculus in modeling biological phenomena: A review. Commun. Nonlinear Sci. Numer. Simul. 2017, 51, 141–159. [Google Scholar] [CrossRef]
  19. Riesz, M. L’intégrale de Riemann-Liouville et le probléme de Cauchy [The Riemann-Liouville integral and the Cauchy problem]. Acta Math. 1949, 81, 1–222. (In French) [Google Scholar] [CrossRef]
  20. Riesz, M. Marcel Riesz Collected Papers; Springer Collected Works in Mathematics; Lars, G., Hormander, L., Eds.; Springer: Berlin/Heidelberg, Germany, 1988; 897p, ISBN 978-3-642-34603-3. Available online: https://0-www-springer-com.brum.beds.ac.uk/gp/book/9783642346033 (accessed on 12 October 2021).
  21. Prado, H.; Rivero, M.; Trujillo, J.J.; Velasco, M.P. New results from old investigation: A note on fractional m-dimensional differential operators. The fractional Laplacian. Fract. Calc. Appl. Anal. 2015, 18, 290–306. [Google Scholar] [CrossRef]
  22. Kwasnicki, M. Ten equivalent definitions of the fractional Laplace operator. Fract. Calc. Appl. Anal. 2017, 20, 7–51. [Google Scholar] [CrossRef] [Green Version]
  23. Tarasov, V.E. Exact discretization of fractional Laplacian. Comput. Math. Appl. 2017, 73, 855–863. [Google Scholar] [CrossRef]
  24. Ben Adda, F. Geometric interpretation of the differentiability and gradient of real order. Comptes Rendus Acad. Sci. Ser. Math. 1998, 326, 931–934. (In French) [Google Scholar] [CrossRef]
  25. Tarasov, V.E. Remark to history of fractional derivatives on complex plane: Sonine-Letnikov and Nishimoto derivatives. Fract. Differ. Calc. 2016, 6, 147–149. [Google Scholar] [CrossRef] [Green Version]
  26. Tarasov, V.E. Fractional generalization of gradient systems. Lett. Math. Phys. 2005, 73, 49–58. [Google Scholar] [CrossRef] [Green Version]
  27. Tarasov, V.E. Fractional generalization of gradient and Hamiltonian systems. J. Phys. A Math. Gen. 2005, 38, 5929–5943. [Google Scholar] [CrossRef]
  28. Engheta, N. Fractional curl operator in electromagnetics. Microw. Opt. Technol. Lett. 1998, 17, 86–91. [Google Scholar] [CrossRef]
  29. Naqvi, Q.A.; Abbas, M. Complex and higher order fractional curl operator in electromagnetics. Opt. Commun. 2004, 241, 349–355. [Google Scholar] [CrossRef]
  30. 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]
  31. Tarasov, V.E. Fractional vector calculus and fractional Maxwell’s equations. Ann. Phys. 2008, 323, 2756–2778. [Google Scholar] [CrossRef] [Green Version]
  32. Bolster, D.; Meerschaert, M.; Sikorskii, A. Product rule for vector fractional derivatives. Fract. Calc. Appl. Anal. 2012, 15, 463–478. [Google Scholar] [CrossRef] [Green Version]
  33. D’Ovidio, M.; Garra, R. Fractional gradient and its application to the fractional advection equation. Electron. J. Probab. 2014, 19, 1–31. [Google Scholar] [CrossRef]
  34. Tarasov, V.E. Toward lattice fractional vector calculus. J. Phys. A 2014, 47, 355204. [Google Scholar] [CrossRef]
  35. 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]
  36. Agrawal, O.M.; Xu, Y. Generalized vector calculus on convex domain. Commun. Nonlinear Sci. Numer. Simul. 2015, 23, 129–140. [Google Scholar] [CrossRef]
  37. Tarasov, V.E. Fractional-order difference equations for physical lattices and some applications. J. Math. Phys. 2015, 56, 103506. [Google Scholar] [CrossRef] [Green Version]
  38. Ortigueira, M.D.; Machado, J.T. On fractional vectorial calculus. Bull. Pol. Acad. Sci. Tech. Sci. 2018, 66, 389–402. [Google Scholar] [CrossRef]
  39. Cheng, J.; Dai, E. Higher-order fractional Green and Gauss formulas. J. Math. Anal. Appl. 2018, 462, 157–171. [Google Scholar] [CrossRef]
  40. D’Elia, M.; Gulian, M.; Olson, H.; Karniadakis, G.E. Towards a unified theory of fractional and nonlocal vector calculus. Fract. Calc. Appl. Anal. 2021, 24, 1301–1355. [Google Scholar]
  41. Luchko, Y. General fractional integrals and derivatives with the Sonine kernels. Mathematics 2021, 9, 594. [Google Scholar] [CrossRef]
  42. Luchko, Y. General fractional integrals and derivatives of arbitrary order. Symmetry 2021, 13, 755. [Google Scholar] [CrossRef]
  43. Luchko, Y. Operational calculus for the general fractional derivatives with the Sonine kernels. Fract. Calc. Appl. Anal. 2021, 24, 338–375. [Google Scholar] [CrossRef]
  44. Tarasov, V.E. General fractional dynamics. Mathematics 2021, 9, 1464. [Google Scholar] [CrossRef]
  45. Tarasov, V.E. General non-Markovian quantum dynamics. Entropy 2021, 23, 1006. [Google Scholar] [CrossRef]
  46. Kochubei, A.N. General fractional calculus, evolution equations, and renewal processes. Integral Equ. Oper. Theory 2011, 71, 583–600. [Google Scholar] [CrossRef] [Green Version]
  47. Kochubei, A.N. General fractional calculus. In Handbook of Fractional Calculus with Applications. Volume 1. Basic Theory; Kochubei, A., Luchko, Y., Tenreiro Machado, J.A., Eds.; De Gruyter: Berlin, Germany; Boston, MA, USA, 2019; Chapter 5; pp. 111–126. [Google Scholar] [CrossRef]
  48. 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., Tenreiro Machado, J.A., Eds.; De Gruyter: Berlin, Germany; Boston, MA, USA, 2019; Chapter 11; pp. 223–234. [Google Scholar] [CrossRef]
  49. Kochubei, A.N.; Kondratiev, Y. Growth equation of the general fractional calculus. arXiv 2019, arXiv:1907.05290. Available online: https://arxiv.org/abs/1907.05290 (accessed on 12 October 2021).
  50. Sonine, N. Sur la generalisation d’une formule d’Abel [On the generalization of an Abel formula]. Acta Math. 1884, 4, 171–176. [Google Scholar] [CrossRef]
  51. Sonin, N.Y. (Ed.) On the generalization of an Abel formula. In Investigations of Cylinder Functions and Special Polynomials; GTTI: Moscow, Russia, 1954; pp. 148–154. Available online: https://ur.zlibcdn2.com/book/2157173/2a8410 (accessed on 12 October 2021). (In Russian)
  52. Sonin, Nikolai Yakovlevich Math-Net.Ru. 2021. Available online: http://www.mathnet.ru/eng/person40324 (accessed on 12 October 2021).
  53. Kochubei, A.N.; Kondratiev, Y. Fractional kinetic hierarchies and intermittency. Kinetic and related models. Am. Inst. Math. Sci. 2017, 10, 725–740. [Google Scholar] [CrossRef]
  54. Sin, C.-S. Well-posedness of general Caputo-type fractional differential equations. Fract. Calc. Appl. Anal. 2018, 21, 819–832. [Google Scholar] [CrossRef]
  55. 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]
  56. Luchko, Y.; Yamamoto, M. The General Fractional Derivative and Related Fractional Differential Equations. Available online: https://www.researchgate.net/publication/346447307 (accessed on 12 October 2021).
  57. Luchko, Y. Special functions of fractional calculus in form of convolution series and their applications. arXiv 2021, arXiv:2108.05668. Available online: https://arxiv.org/abs/2108.05668 (accessed on 12 October 2021).
  58. Luchko, Y. Convolution series and the generalized convolution Taylor formula. arXiv 2021, arXiv:2107.10198. Available online: https://arxiv.org/abs/2107.10198 (accessed on 12 October 2021).
  59. Samko, S.G.; Cardoso, R.P. Integral equations of the first kind of Sonine type. Int. J. Math. Math. Sci. 2003, 2003, 238394. Available online: https://www.hindawi.com/journals/ijmms/2003/238394/ (accessed on 12 October 2021). [CrossRef] [Green Version]
  60. Samko, S.G.; Cardoso, R.P. Sonine integral equations of the first kind in Lp(0; b). Fract. Calc. Appl. Anal. 2003, 6, 235–258. [Google Scholar]
  61. Kinash, N.; Janno, J. Inverse problems for a generalized subdiffusion equation with final overdetermination. Math. Model. Anal. 2019, 24, 236–262. [Google Scholar] [CrossRef]
  62. 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. Available online: https://www.researchgate.net/publication/337449358 (accessed on 12 October 2021).
  63. Hanyga, A. A comment on a controversial issue: A generalized fractional derivative cannot have a regular kernel. Fract. Calc. Appl. Anal. 2020, 23, 211–223. [Google Scholar] [CrossRef] [Green Version]
  64. Giusti, A. General fractional calculus and Prabhakar’s theory. Commun. Nonlinear Sci. Numer. Simul. 2020, 83, 105114. [Google Scholar] [CrossRef] [Green Version]
  65. Tarasov, V.E. General Fractional Calculus: Multi-Kernel Approach. 2021. Available online: https://0-www-proquest-com.brum.beds.ac.uk/docview/2549481424 (accessed on 12 October 2021).
  66. Luchko, Y.; Gorenflo, R. An operational method for solving fractional differential equations. Acta Math. Vietnam. 1999, 24, 207–234. [Google Scholar]
  67. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin, Germany, 2020; 540p. [Google Scholar] [CrossRef]
  68. Hadid, S.B.; Luchko, Y. An operational method for solving fractional differential equations of an arbitrary real order. Panam. Math. J. 1996, 6, 57–73. [Google Scholar]
  69. Moon, P.; Spencer, D.E. The meaning of the vector Laplacian. J. Frankl. Inst. 1953, 256, 551–558. [Google Scholar] [CrossRef]
  70. Tarasov, V.E. Vector calculus in non-integer dimensional space and its applications to fractal media. Commun. Nonlinear Sci. Numer. Simul. 2015, 20, 360–374. [Google Scholar] [CrossRef] [Green Version]
  71. Tarasov, V.E. Geometric interpretation of fractional-order derivative. Fract. Calc. Appl. Anal. 2016, 19, 1200–1221. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tarasov, V.E. General Fractional Vector Calculus. Mathematics 2021, 9, 2816. https://0-doi-org.brum.beds.ac.uk/10.3390/math9212816

AMA Style

Tarasov VE. General Fractional Vector Calculus. Mathematics. 2021; 9(21):2816. https://0-doi-org.brum.beds.ac.uk/10.3390/math9212816

Chicago/Turabian Style

Tarasov, Vasily E. 2021. "General Fractional Vector Calculus" Mathematics 9, no. 21: 2816. https://0-doi-org.brum.beds.ac.uk/10.3390/math9212816

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