General Non-Local Continuum Mechanics: Derivation of Balance Equations

In this paper, mechanics of continuum with general form of nonlocality in space and time is considered. Some basic concepts of nonlocal continuum mechanics are discussed. General fractional calculus (GFC) and general fractional vector calculus (GFVC) are used as mathematical tools for constructing mechanics of media with general form of nonlocality in space and time. Balance equations for mass, momentum, and energy, which describe conservation laws for nonlocal continuum, are derived by using the fundamental theorems of the GFC. The general balance equation in the integral form are derived by using the second fundamental theorems of the GFC. The first fundamental theorems of GFC and the proposed fractional analogue of the Titchmarsh theorem are used to derive the differential form of general balance equations from the integral form of balance equations. Using the general fractional vector calculus, the equations of conservation of mass, momentum, and energy are also suggested for a wide class of regions and surfaces.


Introduction
A general theory of nonlocal continuum mechanics was formally initiated by the papers of Kröner [1] in 1967, and then, by the works of Eringen [2][3][4] in 1972. Nonlocal continuum theory [5,6] is based on the assumption that the forces between material points are a long-range type, thus reflecting the long-range character of inter-atomic interactions. The nonlocal properties of media can manifest themselves in the form of spatial non-locality, fading memory and distributed lag (non-locality in time), fractional long-range interactions and fractional dispersions.
At the present time, three approaches to formulating the mechanics of nonlocal media can be distinguished from a mathematical point of view.

1.
The first approach is based on the use of various integral and integro-differential operators of general form, which are not self-consistent with each other and do not form a calculus. Eringen's book [5] and Rogula's work [6] can be attributed to this approach. As an example, it should also be noted the works [7][8][9] that described continuous media with non-locality in time. Note some modern reviews of various aspects of nonlocal mechanics in [10][11][12][13][14]. The models of non-local media described in these reviews do not assume that integral and integro-differential operators form a calculus. It should be emphasized that the so-called differential (gradient) models, which are based on the differential equations of integer orders, cannot be considered as a tool for describing non-local media from a mathematical point of view. This is due to the wellknown fact that integer-order differential equations are defined in an infinitesimally small neighborhood of the point under consideration and, therefore, are a tool for describing only local media. However, differential (gradient) models can be used to In the framework of the second approach, there are many fractional differential and integral equations that have exact and approximate solutions that describe various problems of non-local mechanics. This fact is an important advantage of the second approach in comparison with the first approach to nonlocal mechanics. Despite this advantage, the second approach has a serious disadvantage. In the models used in the framework of the second approach, only media and materials with power-law nonlocality in space and time were considered. 3.
The third approach, which is proposed in this article, is designed to combine the advantages of the first and second approaches. To describe a wide class of nonlocal continua and media, it is important to use a wide class of operator kernels for which integral and integro-differential operators would form a calculus. Self-consistency of integral and differential operators is provided by non-local analogues of the first and second fundamental theorems of some general calculus. It is also important to have a non-local analogue of the vector calculus, in which non-local analogues of the theorems of standard vector calculus would hold. The possibility of formulating the third approach to nonlocal continuum mechanics is based on the general fractional calculus. The term general fractional calculus has been suggested in article [74] by Kochubei in 2011 (see also [75][76][77]). In works [74,75], the concepts of general fractional derivative (GFD) and general fractional integral (GFI) are proposed and the fundamental theorems of the GFC were proved. The GFC is based on the concept of kernel pairs, which was proposed by Sonin in 1884 work [78,79] (note that sometimes a French transliteration of his surname as "Sonine" is used instead of English transliteration [80]). In this paper, it is proposed to use the Luchko form of the general fractional calculus (GFC) to formulate a continuum mechanics of general nonlocality in space and time. This very important form of the GFC was proposed by Luchko in 2021 [81][82][83] (see also [84,85]). In works [81,82], the GFD and GFI are suggested and the general fundamental theorems of the GFC are proved. Then, in work [86], a general fractional vector calculus (GFVC) was proposed based on the general fractional calculus in the form of Luchko as a generalization of the approach that is suggested in paper [65] (see also Chapter 11 in book [27] (pp. 241-264)).
The GFVC allows us to formulate continuum mechanics with general form of spatial nonlocality. This approach is proposed to be used to derive the balance equations for general non-local media in this article. It should be noted that in 2021 a non-local vector calculus was also proposed by D'Elia, Gulian, Olson and Karniadak [87] as a generalization of the Meerschaert, Mortensen and Wheatcraft approach to FVC. However, this calculus is not related to the general fractional calculus.
In the paper, the Euler approach for description of nonlocal continuum is used. Motion of a continuum obeys the conservation laws for mass, momentum, and total energy. Balance equations for mass, momentum, and total energy of general nonlocal continuum are mathematical formulations of the conservation laws applied to the fixed region of the nonlocal continuum. In this paper, mathematical equations that describe these conservation laws for nonlocal continuum are derived. For this purpose, the fundamental theorems of the GFC are used: (a) To derive balance equations in general integral form, the second fundamental theorem of the GFC is used; (b) The fractional analogue of the Titchmarsh theorem and the first fundamental theorem of the GFC are used to derive balance equation of general nonlocal continuum in the general fractional differential form from the GF integral form; (c) Using the General FVC, the general fractional differential equations for conservation of mass, momentum, and energy are proposed for a wide class of regions and surfaces in general nonlocal media.
In Section 2, properties of the general fractional derivatives and integrals are considered. This section also briefly describes some elements of the general fractional vector calculus. In Section 3, some basic concepts of general nonlocal continuum are discussed. In Section 4, the general fractional continuity equation, which describes a conservation law for mass of general nonlocal continuum, is derived. In Section 5, the general transport equation for momentum that describes conservation law for momentum of general nonlocal continuum is obtained. A derivation of the fractional equilibrium equations for stresses is also suggested. In Section 6, the general fractional transport equation for total energy of continua with general nonlocality is derived. In Section 7, a brief conclusion and a list of main results are given.

Definitions of GFI and GFD
Let us assume that the functions M(t) belongs to the space C −1,0 (0, ∞) and suppose that there exists a function K(t) ∈ C −1,0 (0, ∞), such that the Laplace convolution of these functions is equal to one for all t ∈ (0, ∞). The function ρ(t) belongs to the space C −1,0 (0, ∞), if this function can be represented in the form ρ(t) = t p g(t), where −1 < p < 0 and g(t) ∈ C[0, ∞).

Definition 1.
Let the functions M(t) and K(t)) satisfy the following conditions.
The set of the pairs (M(t), K(t)) that satisfy conditions (1) and (2) is called the Luchko set and is denoted by L.
To define the GFI and GFD, the Luchko approach to general fractional calculus, which is proposed in [81,82], is used.

Definition 2.
Let the pair of the kernels (M(t), K(t)) belongs to the Luchko set L.
The GFI and GFD are connected by the fundamental theorems of general fractional calculus (FT of GFC).

Examples of Kernel Pairs
Let us give examples of the pair of kernels that belongs to the Sonin set S and the Luchko set L. 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 the operators are considered for the kernel pair (M(t), K(t)) belongs to the Sonin set S, then the following operators can also be used A similar situation with the kernel pairs (M(t), K(t)) that belong to the Luchko set L. Let us describe the notation that are used in the Table 1. The function γ(β, t) with β > 0 is the incomplete gamma function The functions J ν (t) and, I ν (t) with ν > 0 are the Bessel and the modified Bessel functions, respectively, Mathematics 2022, 10, 1427 6 of 43 The function Φ(β, α; z) is the Kummer function where (β) k and (α) k are Pochhammer symbols. The function erfc(z) with z ∈ C is the complementary error function The function E α,β [z] is the two-parameters Mittag-Leffler function where α > 0 and β, z ∈ C.  Let us define the general fractional integration and differentiation on [a, b] with 0 ≤ a < b. Definition 3. Let ρ(x ) ∈ 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 is given by the equation The general fractional derivative on [a, b] is given as Let us formulate the fundamental theorems of the GFC for the operator I Theorem 3 (The second fundamental theorem of GFC for I x (M) and GFD D x, * (K) ). Let ρ(x) belongs to the space C 1 −1 (0, ∞), and the pair of kernels (M(t), K(t)) the Luchko set L. Then where x > a ≥ 0.
Theorem 4 (The first fundamental theorem of GFC for I x (M) and GFD D x, * (K) ). Let us assume that the conditions of the first FT of GFC for the operators I x (M) and GFD D x, * (K) are satisfied. Then where it is used vanishes identically, i.e., condition (25) holds for all x ∈ [0, ∞). Then, either one of M(x) or f (x) must vanish identically, i.e., f (x) = 0 for all x ∈ [0, ∞).
Proof. The proof of this theorem without using the theory of functions of a complex variable is proposed in Part VI, Chapter VI, Section 8 in. [88] (p. 180) and Chapter VI, Section 5 in [89] (p. 166). This ends the proof.
The first fundamental theorems of the GFC can be used to derive the GF differential form of general balance equations from the GF integral form of general balance equations. Let us formulate and prove a basic theorem that allows us to implement this derivation of GF differential form of nonlocal balance equations. This theorem can be regarded as an analogue of Titchmarsh's theorem for general fractional calculus.
Theorem 6 (General fractional Titchmarsh's theorem). Let the pair of kernels (M(x), K(x)) belongs to the Luchko set. Let f (x) be real valued function such that f (x) Proof. The action of the generalized fractional derivative on equality (26) leads to the equation since the action of the generalized fractional derivative on zero is equal to zero. Let us use the first fundamental theorem of GFC in the form which holds for all x ∈ [0, ∞), if f (x) ∈ C −1,(K) (0, ∞). Using (28), Equation (26) gives the equality that holds for x ∈ [0, ∞). This ends the proof.

Triple GFI
The following notation will be used where R 1 + = (0, ∞). Let us define concept of the Z-simple region W in R 3 0,+ . Definition 4. Let W be region in R 3 0,+ that is bounded above and below by smooth surfaces S 2,xy , S 1,xy and a lateral surface S z , whose generatrices are parallel to the Z-axis. Let surfaces S 1,xy , S 2,xy be described by the equations where the functions are continuous in the closed domain W xy that is a projection of the region W onto the XY-plane and z 2 (x, y) ≥ z 1 (x, y) for all (x, y) ∈ W xy . 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 3 0,+ .
Definition 5. Let W ⊂ R 3 0,+ be Z-simple domain that is bounded above and below by smooth surfaces S 2,xy , S 1,xy described by Equation (32). Let W xy ⊂ R 2 0,+ be projection of W on the XY-plane such that W xy is Y-simple region in XY-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 ( Let the scalar field ρ(x, y, z) be satisfy the conditions and I Then, the triple general fractional integral (triple GFI) is defined in the form

Surface GFI
The surface GFI over the surface S in R 3 0,+ is defined through the double GFI over areas S xy , S xz , S yz in the XY, XZ, YZ planes, where these areas are projections of the surface S onto these planes.
Let us give definitions of piecewise simple surfaces, vector fields on these surfaces and the surface GFI. Definition 6. Let S be an oriented compact smooth surface in the region 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 and k = 1, . . . , n z such that Let S i,yz , S j,xz , S k,xy be projections of the surfaces S X,i , S Y,j , S Z,k onto the YZ, XZ, XY planes, which can be described by continuous functions Such surfaces S is called piecewise simple surfaces. The set of such surfaces is denoted by P(R 3 0,+ ).
The set of such vector fields F on piecewise simple surface S is denoted by F S (R 3 0,+ ).
Definition 8. Let S be a piecewise simple surface (S ∈ P(R 3 0,+ )) and a vector fields F on this surface S belongs to the set F S (R 3 0,+ ). Then, the surface general fractional vector integral (surface GFI) of the second kind for the vector fiels F ∈ F S (R 3 0,+ ) is defined by the equation Example 1. 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. Let us also assume that the surface S is given by the equation where the point (x, y) ∈ R 2 +,0 changes in area S xy in the XY-plane, bounded by smooth contour ∂S xy and S := {(x, y, z) : z = z(x, y) ≥ 0, (x, y) ∈ S xy ⊂ R 2 0,+ }.
Let F be a vector field F := e z F z (x, y, z), where F z (x, y, z(x, y)) ∈ C −1 (R 2 + ). Then, the surface GFI is defined as If S xy is region of the XY-plane such that where y = y 1 (x) and y = y 2 (x) are continuous functions, Then, the double GFI in Equation (49) is represented as

General Fractional Divergence
In this subsection, the definition of general fractional divergence for R 3 0,+ is given. Let us define sets of vector fields that is used in the definition of the regional general fractional divergence. Definition 9. Let J(x, y, z) be a vector field that satisfies the conditions Then, the set of such vector fields is denoted as F 1 −1,Div (R 3 + ). Let J(x, y, z) be a vector field that satisfies the conditions Then, the set of such vector fields is denoted as In other words, the condition J(x, y, z) ∈ C 1 −1 (R 3 + ) means that all general fractional derivatives of all components of the vector field J(x, y, z) with respect to all coordinates belong to space C −1 (R 3 + ). Let us define the regional general fractional divergence.

Definition 10.
Let J(x, y, z) be a vector field that belongs to the set F 1 −1,Div (R 3 + ) or C 1 −1 (R 3 + ). Then, the general fractional divergence Div (K) W for the region W = R 3 0,+ is defined as One can consider the general fractional nabla operator (del operator) that is defined by equation Using this operator, the GF divergence can be written as where F ∈ F 1 −1,Div (R 3 + ) or F ∈ C 1 −1 (R 3 + ).

Remark 2.
It should be noted that the general fractional divergence can be defined not only as a regional operator for regions W ⊂ R 3 0,+ , but also can be defined for surfaces S ⊂ R 3 0,+ and lines L ⊂ R 3 0,+ . This possibility is due to the fact that this operator is non-local and, in fact, it is an integro-differential operator. The general fractional nabla operator can also be defined for the L, S, W ⊂ R 3 0,+ , that corresponds to the lines, surfaces and regions in R 3 0,+ .

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 fields, for which general fractional Gauss theorem is formulated.
Definition 11. Let W in R 3 0,+ be Z-simple region such that W is a piecewise Y-simple and Xsimple region 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 yz and W Y,j are the Y-simple regions that are described by the functions y = y j,1 (x, z), Then, the set of such vector fields J is denoted as J(x, y, z) ∈ F −1 (∂W). Let vector field J(x, y, z) satisfy the conditions Then, the set of such vector field J is denoted as J(x, y, z) ∈ F 1 −1 (R 3 + ).
Theorem 7 (General fractional Gauss theorem for Z-simple region). Let W in R 3 0,+ be Zsimple 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,xy , S 1,xy , which is described by Equation (32) and a lateral surface S z , whose generatrices are parallel to the Z-axis.
Let the vector field J(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 that can be written as Proof. This theorem was proved in [86]. The general fractional Gauss theorem is also satisfied for regions W ⊂ R 3 0,+ that can be represented as unions of the Z-simple regions W k , which are piecewise Y-simple and X-simple regions in R 3 0,+ .
For parallelepiped regions, a general fractional Gauss theorem can be formulated in the following form.
Theorem 8 (General fractional Gauss theorem for Gauus's theorem for parallelepiped region). Let J x (x, y, z), J y (x, y, z), J z (x, y, z) belong to the function space F 1 −1 (R 3 + ) and the region W ⊂ R 3 0,+ has the form of the parallelepiped If the boundary of W be a closed surface S = ∂W, then For the Cartesian coordinates, the vector field J = J x e x + J y e y + J z e z and the GFI operators expressions of Equation (70) have the form and I where S xy , S xz , S yz are projections of S = ∂W into XY, XZ, YZ planes.
Proof. This theorem was proved in [86].

General Density Function
At the beginning, for simplicity, the one-dimensional case of mass distribution is considered. Let the medium be distributed along the line (beam) R 0,+ = [0, ∞). One can consider the function m(x) ∈ C 1 [0, ∞) that describes the mass of the region describes the standard linear density of the medium. The function ρ s (x) ∈ C[0, ∞) allows us to describe the mass of the region W = [a, b] ⊂ R 0,+ and the mass m(x) by the equations and Substitution of Equation (76) into Equation (74) gives identity by virtue of the first fundamental theorem of calculus, that states the equality Mathematics 2022, 10, 1427 Let us define the real-valued function if the function m(x) belongs to the space Then, Equation (78) takes the standard form (74). If the kernel K(x) belongs the function space C −1,0 (0, ∞) and together with the kernel M x (x) form a kernel pair from the Luchko set. Then, Equation (78) can be written in the form where is the general fractional derivative. The function ρ(x), which is defined by Equation (80), can be interpreted as a general density function. In this case, the function ρ(x) belongs to the space C −1 (0, ∞).
The function (78) can be used to describe mass of the region W = [a, b] ⊂ R 0,+ and the mass m(x) by the equations and If the kernel M x (x) is equal to 1. Then, Equations (81) and (82) take the standard form (75) and (76), respectively.
It should be note that substitution of Equation (82) into Equation (78) can give an identity only under certain conditions on the kernels M x (x) and K x (x) of these operators. This condition is that these kernels form a pair belonging to the Luchko set. In this case, the integral operators are general fractional integrals. Therefore, it is necessary to use function (80) that can be used to describe mass of the region W = [a, b] ⊂ R 0,+ and the mass m(x) by equations and Substitution of Equation (84) into Equation (80) gives identity by virtue of the first fundamental theorem of the general fractional calculus (see Theorems 1 and 4).

General Nonlocal Continuum
The concept of the general density function (80) allows us to consider a general nonlocal continuum (GNC). From the point of view of the mathematical approach, the onedimensional linear GNC is a continuum, nonlocality of which can be described by the pair of kernels from the Luchko sets, but cannot be described by the kernel pairs The general density and other fields of continuum are functions of the time variable t and the space coordinates (x, y, z). In this case, the general density function is ρ = ρ(t, x, y, z).

For the region
the mass of the continuum in the region W 0 is described by equation x , y , z ). (86) Then, the density function is described by the expression Proposed Equations (86) and (87) can be written for the wide class of regions W ⊂ R 3 0,+ . For example, one can use the regions W that can be represented as a union of all the Zsimple region.
As a result, it is possible to formulate the concept of a general nonlocal continuum (GNC) with nonlocality in space and time. Note that the nonlocality in time can be interpreted as a memory in many cases.

Definition 12.
From the point of view of mathematics, the GNC is a continuum, nonlocality of which can be described by the pair of kernels from the Luchko sets, but cannot be described by the kernel pairs From physical point of view, nonlocal continuum is a medium whose behavior at any interior point in space and time depends on the state of all other points in the medium in addition to its own state and the state of external fields.
Using the concept of GNC and the concepts of the generalized density of physical quantities (fields of a non-local continuum medium), balance equations for mass, momentum, and energy can be derived for such media.
Note that a microstructural basis of nonlocal properties of continuum is the assumption that the forces between material points are a long-range type, thus reflecting the long-range character of inter-atomic forces. Microstructural models of general nonlocal continuum can be formulated as models of discrete systems with long-range interactions, frequency and spatial dispersion [5,27].

Nonlocality of Mass in Continuum
For the first time, the continuity equation, which expresses the law of conservation of mass for nonlocal media, was obtained by Wheatcraft and Meerschaert in [90] for the media with power-law nonlocality (see also [91,92]). The proposed equation of mass conservation is a fractional differential equation with the Caputo fractional derivative with respect to space coordinates. The nonlocality in time variable, which can be interpreted as a memory, is not considered in this work and the derivative with respect to time has the first order. In paper [90], authors used a generalization of the standard derivation method for the continuity equation, that is based on the Taylor series with integer derivatives (for example, see Chapters 2 and 3 in [93,94]). In work [90], the fractional generalization of the Taylor formula in the Odibat-Shawagfeh form [95] is used.
The main restrictions of the standard mass the continuity equation are the following: (A) First, the standard derivation method that applies to local media is valid for flux fields, in which flux changes are small and linear or piecewise linear within the fixed region W. (B) The size of the fixed area and the scale of the measurements must be large compared to the scale of the heterogeneity in the medium. These restrictions are necessary due to the fact that changes in the flow in the fixed region W can be approximated by the Taylor series with derivatives of first orders.
It was shown in [90] that the fractional Taylor series can be an exact representation of a nonlinear flow with a power-law nonlinearity, using only the first two terms of this series. Using the fractional Taylor series to describe the change in flux through a fixed region, fractional differential equations describing the conservation of mass with power nonlocality were obtained in [90]. The proposed equation does not contain the described restrictions of the standard mass conservation equation.
However, it should be noted that the use of fractional Taylor series to obtain balance equations for non-local media has significant drawbacks. To do this, consider the generalization of the Taylor series, which was proposed in [95].
The function f (x) with x ∈ [a, b] can be expanded by using the generalized Taylor series with the left-sided Caputo derivatives of order 0 < α ≤ 1 in the form where (89) and R Nα (x, a+) is the remainder term, which can be represented in the form where a ≤ ξ + ≤ x. The generalized Taylor series (88) is applicable to functions f (x) that satisfy the condition for k = 0, 1, 2.

Remark 3.
Note that there is a problem in application of Equation (88), which is caused by the coincidence of the initial and final values in the derivatives D α C;a+ f (a) (see Equation (89)). To derive correct balance equations by using the Taylor series, it is necessary to use a fractional Taylor series at an arbitrary point x 0 , which does not coincide with the initial point of fractional derivative D α C;a+ , i.e., x 0 = a. In this case, it is possible to avoid some restrictions on a possible application of Equation (88). To solve this problem, a generalization of the Taylor formulas (88) is proposed to the case when the Caputo derivative is considered at an arbitrary point x 0 ≥ a (see Sections 3.3 and 4.2 in book [96]). However, in this case, the remainder term will be represented as the sum of the remainder terms of form (90) for different points and therefore such a term cannot be considered as a small value and it cannot be neglected in the general case.
As a result, the proposed approach, which is based on the fractional Taylor series, has some drawbacks due to the properties of the coefficients of the fractional Taylor series and the properties of the media.
(1) First, for local media, it is quite logical that the size of a fixed region tends to zero, and thus an infinitely small region of medium can be used. However, for nonlocal media, this procedure is not entirely logical. (2) Second, the coefficients of the Taylor series are described as fractional derivatives from a to some upper limit, which tends to a even for finite fixed regions. This leads to the fact that the mass balance equation for a nonlocal medium should be described only by fractional derivatives on the infinitesimal interval [a, a + ) with → 0. A possible approach to solve this problem is proposed in [92] for power-law nonlocality. If using the Taylor series to express in terms of fractional derivatives on a finite interval. Then, the remainder (more precisely, the difference between the two remainder terms) will not tend to zero and these terms cannot be neglected (see Remark 3).
(3) Third, the equations were obtained only for the power-law type of nonlocality of the medium. In addition to this, it should also be emphasized that in article [90] other laws of conservation of nonlocal media were not derived.

Remark 4.
It should be noted that the use of the fractional derivative of the non-integer order is actually equivalent to using an infinite number of derivatives of integer orders, which assume arbitrarily large values. For example, the Riemann-Liouville derivative [15,18] can be represented in the form of the infinite series for analytic functions on (a, b), (see Lemma 15.3 in [15]). Using Equation 2.4.6 in [18], of the form and Equation (92), the Caputo derivative of the non-integer order can also be represented by the infinite series of derivatives of integer orders. Therefore, the use of fractional derivatives means that a contribution of derivatives of all integer orders is taken into account with weights of the power-law type. This fact creates additional difficulties for using the fractional Taylor series to obtain balance equations for non-local media. Moreover, this approach does not allow one to obtain balance equations in integral form and for a finite region of a non-local medium.
By virtue of the above arguments, it is important to mathematically and consistently derive all the conservation laws and obtain balance equations for nonlocal media with a general form of nonlocality. This conclusion became possible due to the creation of a general fractional calculus in the form proposed by Luchko in [81][82][83][84][85] and a general fractional vector calculus in [86].

General Nonlocal Continuity Equation
In this section, the motion of a nonlocal medium in a stationary reference frame with a fixed Cartesian coordinate system is considered. The Euler approach for description of continuum is used. Balance equation for mass of general nonlocal continuum is a mathematical formulation of the conservation law applied to the fixed region of the nonlocal continuum.
Let us select an region W ⊂ R 3 0,+ of continuum in the form of a parallelepiped AB . . .
Let point S ∈ W have coordinates (x, y, z) such that x ∈ [a, b], y ∈ [c, d] and z ∈ [e, f us consider the density ρ = ρ(t, x, y, z), the velocity V = V(t, x, y, z) and the flow J(t, x, y ρ(t, x, y, z) V(t, x, y, z) of mass at the point S(x, y, z).
Let us consider the density ρ = ρ(t, x, y, z), the velocity V = V(t, x, y, z) and the flow J(t, x, y, z) = ρ(t, x, y, z) V(t, x, y, z) of mass at the point S(x, y, z).
In this section, a connection between the density and the velocity of a continuous medium in the nonlocal case is derived by assuming the velocity V(t, x, y, z) and density ρ(t, x, y, z) to be continuous functions of time and coordinates that belong to the function space C 1 (R 4 0,+ ). In general, the following conditions are assumed and where J k (t, x, y, z) = ρ(t, x, y, z) V k (t, x, y, z).
The function spaces C {n,m} −1 (R 4 + ) are defined in the next subsection. Note that C(0, ∞) ⊂ C −1 (0, ∞). Therefore the continuous functions of time and coordinates can be considered as special cases of the general approach.

Function Spaces for Derivation of Balance Equations
Let us define two function spaces, which will be used in the theorems about the balance equations for general nonlocal media. The first type of spaces (C {n,m} −1 (R 4 + )) will be used in derivation of the balance equations in the general integral form. The second function space (C −1,(K) (R 4 + ) will be used in derivation equations in the general differential forms from integral form. (1) Property with respect to time variable t: This condition means that f (t, x, y, z) ∈ C n −1 (0, ∞) for all (x, y, z) ∈ R 3 0,+ . (2) Property with respect to the space variables (x, y, z): This condition means that f (t, x, y, z) ∈ C m −1 (R 3 + ) for all t ∈ R 0,+ .

Mass of Nonlocal Continuum
Let the density ρ = ρ(t, x, y, z) as a function of time t ≥ 0 belong to the space C 1 −1 (0, ∞) for all (x, y, z) ∈ R 3 0,+ . The mass of nonlocal medium in the parallelepiped region W ⊂ R 3 0,+ is given by the equation A discussion of the justification for equations of mass nonlocality from the physical point of view is considered in Section 3.3.
If the kernels M k (x k ) with k = 1, 2, 3 have the form where 0 < α k ≤ 1, Equation (109) The mass of the medium in the region W at the time t = t 1 ≥ 0 is determined by equation At the next moment in time t = t 2 = t 1 + ∆t the mass of the selected region is equal to Thus, due to a change in the density of the medium in a given fixed region W, the mass during the time ∆t = t 2 − t 1 is changed by the value The fundamental theorem of the GFC [81] in the form can be used in Equation (114) to obtain the expression Equation (116) describes the changes of mass due to change of medium density in the fixed region W ⊂ R 3 +,0 during the time ∆t.

Mass Flow of Nonlocal Continuum
Let us consider the density ρ(t, x, y, z), the velocity V(t, x, y, z) and the flow J(t, x, y, z) = ρ(t, x, y, z) V(t, x, y, z) of mass at the point S(x, y, z) ∈ W.
The mass change in the fixed region W will occur due to its mass transfer across the boundaries of the region. Let us now calculate the change in mass in the area W during the time ∆t = t 2 − t 1 , when it is transferred across the boundaries of the region W.
The medium mass, which is flow through the face W a = ADD 1 A 1 for the time ∆t in the direction of the OX-axis, is given by equation where J x (t, a) is the flow of mass through the face W a of the parallelepiped W such that Here J x (t, x, y, z) = ρ V x , where ρ = ρ(t, x, y, z) and V x = V x (t, x, y, z). The mass, which is transported through the opposite face W b = BCC 1 B 1 in time ∆t in the direction of the OX-axis, is given by the equation Let us assume that the mass flowing into the region is positive and the mass flowing out negative. Then, the change in mass in the region W, when the mass is transferred through the W a and W b faces (see Figure 1) perpendicular to the OX-axis, is equal to where The fundamental theorem of the GFC in the form Considering similarly, it is possible to obtain the mass change in the region W through the faces that are perpendicular to OY and OZ axes. Reasoning in a similar way, the mass that is transferred through the faces perpendicular to the axes OY and OZ, has the form, respectively: Equations (123)-(125) describe the change of mass in the region W ⊂ R 3 0,+ .

Mass Balance Equation
The sum of values ∆m t , ∆m x , ∆m y and ∆m z gives the equation that describes the mass conservation. Using Equations (116), (123)-(125), the balance Equation (126) can be represented in the form As a result, the following theorem is proved.
Equations (127) and (128) must be satisfied for any finite time intervals ∆t = t 2 − t 1 , where t 2 > t 1 ≥ 0 and any finite sizes of the parallelepiped, This fact can be used to obtain the balance equation for general nonlocal continuum in a differential form.
To derive balance equation of general nonlocal continuum in the fractional differential form, it is possible to use the fractional analogue of the Titchmarsh Theorem (see Theorem 6). This theorem and the first fundamental theorems of GFC are used to derive GF differential form of general balance equations from the GF integral form of general balance equations. As a result, the balance equation is obtained in the form which can be called the general nonlocal continuity equation. Equation (129) is the desired equation that describes the balance of mass in general non-local media. As a result, the following theorem is proved for a differential form of the mass balance in general nonlocal continuum.
Theorem 10 (Mass balance in GF differential form). Let the function ρ = ρ(t, x 1 , x 2 , x 3 ) and J n = J n (t, x 1 , x 2 , x 3 ) (n = 1, 2, 3) satisfy the conditions for all n = 1, 2, 3. Then, the GF differential balance equation of mass of general nonlocal continuum in the parallelepiped region W ⊂ R 3 0,+ has the form Remark 5. It should be emphasized that for the non-local media the differential balance equation depends on the region W and in fact is an integro-differential equation.
Equation (132) is a mathematical expression of the existence of continuum flow without the formation of discontinuity surfaces. The proposed continuity equation is valid for any nonlocal continuous medium, viscous and compressible, with steady and unsteady motions.
Equation (132) is the basic equation of the mechanics of nonlocal medium, the nonlocality of which can be described by kernels belonging to the Luchko set.

General FVC Form of Mass Balance Equation
For general nonlocal media, the continuity equation, which is derived for parallelepiped region W, can be generalized to a wide class of domains and surfaces. This generalization can be realized by using the theorems of the general fractional vector calculus (General FVC).
Using the general FVC, Equation (129) can be written in the vector form where Div

(K)
W is the general fractional divergence, which is defined by Equation (57) and J(t, x, y, z) = e x J x (t, x, y, z) + e y J y (t, x, y, z) + e z J z (t, x, y, z), J(t, x, y, z) = ρ(t, x, y, z) V(t, x, y, z). Equations (133) and (136) is the fractional differential equation of continuity that represents conservation law of mass for the general nonlocal continuum. Equation (133) is valid for fractional nonlocal continuum in steady and unsteady motion, and it is the fundamental equation of general nonlocal continuum mechanics for the wide class of the Sonin kernel pairs from the Luchko set.

Remark 6.
For the general fractional divergence the standard properties of divergence is violated. For example, the well-known equality is violated for the operator Div W in the general case. For the general fractional vector calculus [86], there is the inequality where r = x e x + y e y + z e z . Inequality (138) is a consequence of the violation of the standard product rule (the Leibniz rule) for fractional derivatives of non-integer orders and all the types of the general fractional derivatives. As a result, some equivalent forms of representations of the standard balance equations are not equivalent in the nonlocal case.

General Fractional Equation for Momentum
Let us consider motions of general nonlocal medium, assuming that the velocity, density, pressure and, mass forces are continuous functions of time and position. The Euler approach for description of continuum is used. Balance equation for momentum of general nonlocal continuum is a mathematical formulation of the conservation law applied to the fixed region of the nonlocal continuum. In the Cartesian coordinate system OXYZ, an element of continuum is selected in the form of the parallelepiped with sizes ∆x, ∆y, ∆z (see Figure 1). Let ρ = ρ(t, x, y, z) be the density and V = V(t, x, y, z) be the velocity at the point S(x, y, z) ∈ R 3 +,0 . Then, the momentum of continuum particles in the given fixed parallelepiped area W, at a time t, is desccribed by the equation The changes of the momentum P(t) are caused by changes of the density and the velocity of the particles that came to the region W, at the time t 1 and t 2 , such that In the given fixed parallelepiped area W, the change of momentum during the time ∆t = t 2 − t 1 is written as The fundamental theorem of the GFC in the form x, y, z) = J(t 2 , x, y, z) − J(t 1 , x, y, z) (144) is used to obtain the equation The momentum in the parallelepiped region W is varied by the following: (a) the transfer of momentum across the boundary; (b) the momentum of surface forces; (c) the momentum of mass forces.
These changes are considered in the following subsections.

Momentum Transfer across Boundaries
Let us consider the change of momentum due to its transport across the boundary. The momentum transferred through the face W a during the time ∆t in the positive direction of the axis OX is The momentum transferred through the face W b during the time ∆t can be represented by The change of momentum by moving it across the faces, which is perpendicular to the OX-axis, is found if Equation (147) Then, the fundamental theorem of the GFC in the form is used to obtain the equation Similar arguments is found for the change of momentum in the region W, when moving nonlocal continuum through a pair of faces perpendicular to the axes of OY and OZ, in the form Summing relations (150)-(152), one can obtain the change of momentum in the parallelepiped region W due to its transport through the boundary where

The Momentum of the Mass Force
In classical mechanics of point particles without nonlocality in time (memoryless), the momentum of the force is described as For the case of nonlocality in time, the momentum of the force in given by the equation Let us consider the momentum of the mass force for non-local continuum. The density of the mass (volume) force is denoted as F = F(t, x, y, z). Then, the momentum of the mass forces is the following where it is assumed that the condition ρ(t, x, y, z) F(t, x, y, z) ∈ C −1 (R 4 0,+ ) holds for all (t, x, y, z) ∈ R 4 0,+ .

The Momentum of the Surface Force
In this subsection, the momenta of surface forces are described. Let us consider the stress p x for the face with the normal that is parallel to the OX-axis and positive direction of the face. On the opposite face the stress is p −x . Additionally, p x = −p −x . The momentum of the force on the face W a is On the opposite face W b , the momentum of the force is The total momentum for the faces perpendicular to the OX-axis, is the sum of Equations (158) and (159) in the form: where one can take into account the property p x = −p −x and W x = W a W b . The fundamental theorem of the GFC in the form is used to obtain the equation Similarly, one can obtain the momenta of the surface forces acting on the faces that are perpendicular to the OX, OY and OZ axes in the form for k = 1, 2, 3 and Summing Equations (163), the momentum of the surface forces is described by the equation where p k can be represented through the stress tensor σ kl .

General Fractional Momentum Equation
The balance equation of momentum has the form Summing expressions (153), (156) and (164), equating the resulting sum to relation (145), one can obtain Equation (166) describes transfer of momentum in terms of stresses for fractional nonlocal continuum.
Equation (166) can be written in the compact form by using x 1 = x, x 2 = y, x 3 = z. As a result, the following theorem can be formulated.

Theorem 11 (Momentum balance in GI integral form). Let the conditions
are satisfied for all n = 1, 2, 3. Then, the GF integral balance equation of momentum of general nonlocal continuum in the parallelepiped region W ⊂ R 3 0,+ has the form Equation (166) should be satisfied for all the parallelepiped regions W ⊂ R 3 0,+ and all time intervals [t 1 , t 2 ] ⊂ R 0,+ . To derive balance equation for momentum of general nonlocal continuum in the GF differential form, the fractional analogue of the Titchmarsh Theorem (see Theorem 6) and the first fundamental theorems of GFC should be used. As a result, the general balance equation for momentum in integral form (166) gives the differential form where σ kl is the stress, Equation (171) can be written in the component form for all t > 0 and x k > 0. Equation (174) can be written in the compact form. As a result, the following theorem for momentum balance in general nonlocal continuum was proved.
Theorem 12 (Momentum balance in GF differential form). Let the following conditions be satisfied for all n, l = 1, 2, 3. Then, the GF differential balance equation of momentum of general nonlocal continuum in the parallelepiped region W ⊂ R 3 0,+ has the form For the non-local media the differential balance Equation (178) depend on the region W, since these equations are integro-differential equations. Equations (174) and (178) is the desired equation describing the balance of momentum in general non-local media.

GF Divergence of a Dyadic Product
At the beginning, let us consider the standard local case. The nabla operator and the dyadic product of the vectors can be used to get ∇, where the brackets , denotes the scalar product and the basis vectors e k are fixed.
Then, using the equality and the chain (Leibniz) rule Therefore, The vector notation for local case can be used to get the equality that can be interpreted as a chain rule for the standard nable operator. For general fractional derivatives, the standard product (Leibniz) rule is violated. Therefore, for nonlocal case the following inequality has the form Equations (181) and (183) gives The general fractional analog of Equation (190) has the form Let us consider the nonlocal case and general fractional differential vector operators.
The general fractional divergence of dyad product can be defined by equation where a fixed frame of reference and a fixed Cartesian coordinate system are used.
It should be noted that the general fractional divergence does not satisfy the standard properties of the divergence, such as In the general case, the following inequality shoud be taken into account for general fractional divergence. Inequality (194) is a consequence of the violation of the standard product rule (the Leibniz rule) for fractional derivatives of non-integer orders and general fractional derivatives of all the types.

General FVC Form of Momentum Balance Equation
The continuity equation of thegeneral nonlocal media, which is derived by using parallelepiped region W, can be generalized to a wide class of domains and surfaces by using the theorems of the general fractional vector calculus. Equation (174) contains general fractional derivatives D t, * and and x k ∈ [0, ∞) for k = 1, 2, 3.
Using the general fractional vector calculus, Equation (174) can be written in the vector form where P σ is the stress tensor In a more compact form, Equation (195) can be written as For the local case, Equations (197) give the standard equations with the partial derivative ∂ ∂t and divergence (div) instead of the general fractional operators D t, * (K t ) and Div

General Fractional Equilibrium Equations for Stresses
In this subsection, the general fractional equilibrium equations for the stress tensor of general nonlocal continuum is derived. For simplification, one can consider a twodimensional elastic isotropic non-local continuum. Let us consider a rectangular region W 2 ⊂ R ,+ (see Figure 2) that is given as In this subsection, the general fractional equilibrium equations for the stress tensor of general nonlocal continuum is derived. For simplification, one can consider a two-dimensional elastic isotropic non-local continuum. Let us consider a rectangular region W 2 ⊂ R ,+ (see Figure 2) that is given as In the Figure 2, the following notations are used and σ x+ xy = σ xy (b, c), σ xy = σ xy (x, y), σ y+ xy = σ xy (a, d), σ xy = σ xy (x, y).
The force equilibrium condition gives the equation In the Figure 2, the following notations are used and The force equilibrium condition gives the equation For direction, that is parallel to the OX-axis, the fundamental theorem of GFC for the stress can be used in the form Then, the equation can be written as The fact that Equation (208) must be satisfied for any rectangular regions W 2 ⊂ R ,+ and the general fractional analogue of Titchmarsh's Theorem 6 leads to the general fractional differential equation For the power-law non-locality K j (x j ) = h 1−α j (x j ), Equation (209) takes the form For α x = α y = 1, Equation (210) takes the standard form Similarly for the direction that is parallel to the OY-axis, the force equilibrium condition gives the equation The two-dimensional case can be generalized to the three-dimensional case. As a result, the general fractional equilibrium equations take the form where

General Fractional Equation for Total Energy
The Euler approach for description of continuum is used to derive the energy balance equation. Equation for total energy of the general nonlocal continuum is a mathematical formulation of the conservation law for energy applied to the fixed region of the nonlocal continuum. This law states that the change of the total (kinetic and internal) energy ∆E t in the fixed region of nonlocal continuum is described by equation where ∆E S is the energy transport across the boundaries S = ∂W of the region W ⊂ R 3 0,+ ; A M is a work of mass forces (body forces) and A S is a work of surface forces and Q S is a heat supplied through the surface S. In this section, it will be assumed that the volumetric input amount of energy is equal to zero Q M = 0.
Let us consider a region W ⊂ R 3 0,+ in the form of the elementary parallelepiped that is presented in Figure 1. It will be supposed that the velocity V(t, x, y, z), density ρ(t, x, y, z), stress and mass forces are known functions at the point S(x, y, z) ∈ W, which are continuous functions of coordinates (x, y, z) and time t. In general, it is possible for a wider class of functions or function space.
The value V 2 /2 is the kinetic energy per unit mass and e i is the internal energy per unit mass of the medium in a given time t ≥ 0. The function E(t, x, y, z) = e i (t, x, y, z) describes the total energy per unit mass of the medium at t ≥ 0. The value (ρ V 2 /2)(t, x, y, z) is the density of the kinetic energy of the nonlocal continuum in a given time; (ρ e i )(t, x, y, z) is the density of the internal energy of nonlocal continuum at t ≥ 0. The sum of the densities of internal energy and the kinetic energy is called the density of total energy where ρ, e i and V are functions of the time variable t ≥ 0 and the space coordinates (x, y, z) ∈ R 3 0,+ . The kinetic energy E kin (t) at time t ≥ 0 is defined by equation The internal energy E int (t) at time t ≥ 0 is defined by equation The total energy E(t) of the fixed region W of the medium at time t ≥ 0 is given by the equation W [x, y, z](ρE)(t, x, y, z).
In the general case, the following condition can be considered for Equation (219). For t = t 1 and t = t 2 , the total energy is described as Therefore, the change of the total energy in a time ∆t = t 2 − t 1 is written as Using the fundamental theorem of the GFC in the form Equation (224) takes the form Change of the total energy ∆E t in the region W can be caused by the following processes: (a) the energy transport across the boundaries of the region ∆E S ; (b) the work of the mass forces ∆A F ; (c) the work of the surface forces A S ; (d) the heat flux through the boundary Q S (internal heat sources are not considered).
In the following sub-sections, expressions describing the contribution of these processes to the energy balance in the nonlocal medium will be obtained.

Change of Total Energy Due to the Transfer across Boundaries
Change of the total energy due to the transfer across boundaries of the selected fixed region W can be obtained analogously as was done above for calculating the change in momentum. The energy transfer through the boundary perpendicular to the OX-axis has the form Then, the sum of the terms, that take into account the mass transfer across all the boundaries of the fixed region W, has the form where (ρ E)(t, x 1 , x 2 , x 3 ) is the density of total energy that is given in (217).

Work of Mass Forces
In classical mechanics of point particle with nonlocality in time, which is characterized by the kernel M(t), the work of the force f(t) is described by equation where v(t) is the velocity of the point particle.
In the continuum mechanics, the work ∆A F of the force for ∆t = t 2 − t 1 is described by the equation where F is the force per unit mass fo the medium, ρ F is the density of force and is the scalar product of the vector fields that can be interpreted as a power of the mass force per unit mass (the density of power).

Work of Surface Forces
Let us obtain an expression for the work of the surface forces. The vector p x denotes the stress vector at the face AA 1 D 1 D, Let us denote the stress vector on the face W a with the direction of the normal coinciding with the positive direction of the OX-axis as p x . For the opposite direction of the normal, it should be used p −x . Furthermore, p x = −p −x .
Let us find the total work ∆A S of the surface forces acting on the faces of the parallelepiped region W ⊂ R 3 0,+ . The work of the surface force acting on the area W a with stress p −x is For the face W b , the work of the surface force with stress p x is The sum of the works of the surface forces acting on the faces W a and W b that are perpendicular to OX, can be written as where p −x = −p x and W x = W a W b are taken into account. The fundamental theorem of the GFC in the form can be used to get In Equation (237), the stress vector p x can be represented in the form Therefore, the scalar product can be written as and the work (237) is represented as the sum This allows us to assert that the total work ∆A S , which is caused by the surface forces, is defined by the normal stresses σ xx , σ yy , σ zz and the tangential stresses σ xy , σ yz , σ xz .
Considering the work of the forces acting on the faces W a and W b that are perpendicular to the OX-axis, one can obtain the work of the normal forces for the time interval [t 1 , Therefore, The funfamental theorem of the GFC is used to get where it is assumed that It should be noted that the standard Leibniz rule for GF derivatives is not satisfied. Therefore, the inequality is satisfied in the general case.
The work of tangential forces on the faces W a and W b , which are perpendicular to the OX-axis, is equal to Thus, the total work ∆A x of the surface forces acting on the faces that are perpendicular to OX, for the time interval [t 1 , t 2 ], is equal where Similarly, one can obtain the works for the faces that are perpendicular to the OY-axis and OZ-axis. For all pairs of faces, the total work ∆A j is described by the equation As a result, the total work, which is caused by the surface forces, is described in the form where

Change of Total Energy by the Heat Flow
Let us find a change of the total energy due to the heat flow. We denote the specific vector of the density of heat flow by q(t, x, y, z). Then, the heat flow through the face W a , which is perpendicular to the OX-axis, is For the face W b , the heat flow through the face is The change of the total energy by the heat flux transported through the faces W x = W a W b that are perpendicular to the OX-axis, is given by Using the fundamental theorem of the GFC in the form and Equation (255) can be written as For all pairs of faces, which are perpendicular to all the axes, the change of the total energy is described by the equations Summing Equation (259), the change of the total energy by the heat flow is obtained in the form where q j = q j (t, x 1 , x 2 , x 3 ) (j = 1, 2, 3) are components of the density of heat flow.

Change of Total Energy of All Sources
The balance equation for the total energy has the form Summing relations (228), (230), (251) and (260) that have the form and equating the resulting sum to relation (226) in the form the balance equation for general nonlocal continuum takes the integral form As a result, the following theorem was proved.
Theorem 13 (Energy balance in GF integral form). Let the following conditions be satisfied Then, the GF integral balance equation of total energy of general nonlocal continuum in the parallelepiped region W ⊂ R 3 0,+ has the form The fact that Equation (274) must be satisfied for any regions in the form of a parallelepiped and the general fractional analogue of Titchmarsh's Theorem 6, can be used to obtain the balance equation for the total energy in the general fractional differential form Equation (275) can be rewritten as Equation (276) can be represented in the compact form. As a result, the theorem for the total energy balance in general nonlocal continuum can be formulated in the following form.
Then, the GF differential balance equation of total energy of general nonlocal continuum in the parallelepiped region W ⊂ R 3 0,+ has the form D t, * where ρ E = ρ e i + 1 2 ρV 2 .
Note that for the non-local media the differential balance equation for total energy (279) depends on the region W and in fact is an integro-differential equation.

Remark 8.
Let us consider the case of a memoryless continuum. In this case, the general fractional derivative with respect to time variable is a standard derivative of the first order. Then, Equation (279) takes the form If all general fractional derivatives with respect to coordinates are standard derivatives of the first order, then Equation (281) takes the form of the standard balance equation of total energy for local media.

General FVC Form of Energy Balance Equation
The balance equations of general nonlocal media, which are derived by using parallelepiped region W, can be generalized to a wide class of domains and surfaces by using the theorems of the general fractional vector calculus [86].
Using the general fractional vector calculus (General FVC), Equation (279) can be written in the form where Div (K) The product (P σ , V) of the tensor P σ and the vector V is the vector that has the components The components of the tensor P have the form The tensor P σ consists of nine components σ jk (t, x 1 , x 2 , x 3 ) that define the stress state at the point (x 1 , x 2 , x 3 ) ∈ W ⊂ R 3 0,+ .

Spatial Power-Law Nonlocality without Memory
Let us consider the case of memoryless continuum with the spatial power-law nonlocality. In this case, the general fractional derivatives are the Caputo fractional derivatives and Equation (276) Let us substitute expression (280) into Equation (288). The product rule (the Leibniz rule) for the time-derivative of the first order and the continuity equation, can be used to get Equation (289) describes the change of total energy in continuum with the power-law spatial nonlocality. It is valid for any of the stress-strain rates for the power-law non-local continua without memory.

Conclusions
In the suggested paper, the theory of nonlocal continua with general type of nonlocality is developed. The non-local continuum mechanics is formulated by using the general fractional calculus (GFC) in the Luchko form [81][82][83] and general fractional vector calculus (GFVC) in the form that is proposed in [86].
The integral and differential equations, which describe the conservation laws of mass, momentum, and energy for the general nonlocal continuum, are derived. The integral forms of the balance equations of general nonlocal continuum are derived by using the second fundamental theorems of general fractional calculus. The differential forms of the balance equations of general nonlocal continuum are derived by using the proposed fractional analogue of the Titchmarsh theorem and the first fundamental theorems of the GFC. Using the GFVC, the general equations for conservation of mass, momentum, and energy are obtained for continuum with general form of nonlocality in space and time.
In this paper, some basic concepts of general nonlocal continuum are discussed in Section 3. The main results of this article are derivations of the balance equations for mass, momentum, and energy, which describe conservation laws for general nonlocal continuum. These equations are derived in the two forms: general fractional integral and differential forms. In this article, it is proved that the fundamental theorems of calculus, which describe the relationship between differential and integral operators, make it possible to derive balance equations.
Let us briefly point out the main ideas and methods for deriving balance equations for nonlocal media: • To derive general balance equation in the GF integral form, the second fundamental theorems (Theorems 2 and 3) of the GFC is used; • The first fundamental theorems (Theorems 1 and 4) of the GFC and the proposed fractional analogue of the Titchmarsh theorem (Theorem 6) are used to derive differential form of general balance equations from the integral form of balance equations; • Using the general fractional vector calculus, the balance equations are suggested for a wide class of regions and surfaces of the general nonlocal continuum.
A more detailed list of the main results in the form of theorems and equations obtained in this paper would include:

•
The mass balance equation for continuum with general space and time nonlocality:

-
The mass balance equation in the GF integral form for general nonlocal continuum is given by Theorem 9 for regions in the form of a parallelepiped.

-
The mass balance equation in the GF differential form for general nonlocal continuum is given by Theorem 10 for regions in the form of a parallelepiped.

-
Using the general fractional vector calculus, the mass balance equation in GF differential form is described by Equation (133) for a wide class of regions and surfaces.
• The momentum balance equation for continuum with general space and time nonlocality:

-
The momentum balance equation in the GF integral form for general nonlocal continuum is given by Theorem 11 for regions in the form of a parallelepiped.

-
The momentum balance equation in the GF differential form for general nonlocal continuum is given by Theorem 12 for regions in the form of a parallelepiped.

-
Using the general fractional vector calculus, the momentum balance equation in GF differential form is described by Equation (195) for a wide class of regions and surfaces.
• The energy balance equation for continuum with general space and time nonlocality:

-
The energy balance equation in the GF integral form for general nonlocal continuum is given by Theorem 13 for regions in the form of a parallelepiped.

-
The energy balance equation in the GF differential form for general nonlocal continuum is given by Theorem 14 for regions in the form of a parallelepiped.

-
Using the general fractional vector calculus, the energy balance equation in GF differential form is described by Equation (282) for a wide class of regions and surfaces.
It should be noted that in order to solve the equations of general nonlocal continuum mechanics, one can use the methods of general operational calculus that was proposed by Luchko in [83,85].