will be discussed.
The fractional Λ-space (X, Y, Z) is defined by,
ith a=b=0 , Eq.(17) yields,
For γ=0.6, the surface Z in the Λ-fractional space is defined by
Z=0.947X1.714Y1.714 (Figure 1) (19)
and it is shown in Figure 2.
Figure 2: The surface Z in the Λ- fractional space.
Further, the tangent space of the surface with γ=0.6, at the point X=Y=0.6 is defined by,
and finally, the equation of the tangent space in the Λ-fractional space,
Z=0.164+0.469(Χ-0.6)+0.469(Υ-0.6) (Figure 3). (21)
Figure 3: The surface with the tangent space in the Λ-fractional space.
The corresponding surface in the initial space to the tangent plane in the Λ-fractional space is defined by,
The surface defined by Eq.(22) is shown in Figure 4
Figure 4: The surface with its tangent surface at the point (x=y=0.8106) at the initial space.
It seems that the initial surface and the tangent surface corresponding to the tangent space at the Λ-space have almost a common tangent plane in the initial space since their mathematical expressions are different but quite close.
The fractional field theorems
The conventional field theorems are expressed by:
- Green’s theorem: Let Qx(x,y), Qy(x,y), be smooth real functions in a domain Ω, with its boundary a smooth closed curve
. Then,
+Corollary: When Qx(x,y), Qy(x,y), are derived by a potential function Φ(x,y) with
the RHS of Eq.(23) becomes zero. That means that the curvilinear integral along a closed smooth boundary is zero.
b. Stoke’s theorem: For a smooth vector field F defined on a simple surface Ω with the boundary , Stoke’s theorem is expressed by,
where,
denotes the scalar product.
c. The Gauss’ (divergence) theorem: For a space region Ω with smooth surface boundary
, the volume integral of the divergence of a vector field F over Ω is equal to the surface integral of F over the boundary
:
Although the field theorems are valid in the fractional Λ-space, they are not necessarily valid in the initial space. Nevertheless, the Λ- space results may be pulled back into the initial space.
The balance principles
Before discussing the general problem of hydrodynamics with its specific principles, it is pointed out that geometry and, consequently, mechanics have meaning in theΛ-fractional space. The results of the various analyses in the Λ-fractional space are transferred as functions in the initial space. No derivatives of the Λ-space have any meaning in the initial space as derivatives. However, they may have a meaning as functions.
Almost all balance principles are based on Reynold’s transport theorem. Hence the modification of that theorem, just to conform to fractional analysis is presented. The conventional Reynold’s transport theorem is expressed by:
For a vector field A applied upon region W with boundary ∂W and Vn is the normal velocity of the boundary ∂W.
a. The balance of mass: The conventional balance of mass, expressing the mass preservation is expressed by:
Recalling the fractional Reynold’s Transport Theorem, we get:
Since Eq. (28) is valid for any volume V, the continuity equation is:
where Div is defined in the Λ-fractional space. That is the continuity equation expressed in fractional form.
b. Balance of linear momentum principle: It is reminded that the conventional balance of linear momentum is expressed in continuum mechanics by:
where V is the velocity, T is the traction on the boundary and B is the body force per unit mass. Likewise, that principle in fractional form is expressed by:
Hence the equation of linear motion, expressing the balance of linear momentum is defined by,
In Eq.(32) the term T is introduced by the pressure P. That is affected by adding the term –PI. Therefore Eq.(32) becomes:
Following similar steps as in the conventional case, the balance of rotational momentum yields the symmetry of the Cauchy stress tensor.
First law of thermodynamics
This law occurs from Eq.(42). Firstly we replace
with,
Then the equation of conservation of energy, taking into consideration the principle of conservation of linear momentum, Eq.(33), yields the rate of change of internal energy in the Λ-fractional space as the sum of stress power plus the heat added by the continuum. The vector C is defined in the Λ-fractional space is defined as the heat flux per unit area per unit time by conduction and Z is per the radiant heat constant unit mass per unit time. Further, The caloric equation of state is expressed by e = e(R,T). Finally, the equation of the conservation of energy is defined by, see Ref. [46 ].
The various results in the Λ-fractional space should be transferred as functions in the initial space.
Navier-stokes equations
The Navier-Stokes Equations consist of the following equations:
- Balance of mass (continuity) Eq.(29).
- Balance of linear momentum, Eq.(33)
- The first law of thermodynamics, Eq (35)
- Constitutive equations
e) The kinetic equation of state:
f) The Fourier law of heat conduction:
g) The caloric equation of state:
The system of these sixteen equations contains 16 unknowns therefore it is determinate. Usually, the studies concerning the Newtonian fluid, are restricted to the equations of conservation of mass (29), linear momentum (33), and constitutive equations (36). The solution in the Λ-fractional space is transferred to the initial space.
The Euler and Bernoulli equations
The Euler equations occur from Eq.(30) for non-viscous fluid, for which T=Ο. This holds because the Lamé constants λ, and μ are considered zero in this case. Therefore, Eq.(33) takes the form:
Let's assume a barotropic condition R=R(P). A pressure function may be defined as:
Moreover, the body force may be described by a potential function:
with these two conditions the conservation of linear momentum, Eq. (33), becomes:
If Eq.(43) is integrated along a streamline, the result is the Bernoulli equation in the fractional form:
Any problem may be solved following the conventional way in the Λ-fractional space, however, the results should be transferred into the initial space.
Darcy’s law
Consider a viscous fluid, flowing in a straight pipe of a constant circular cross-section of inner radius R and cross-section of
with perimeter
If L is the length of the considered segment of the pipe defined by positions x=0 and x=a, then l=b-α. Transferring the length into the Λ-fractional space,
Denoting Q(X, T) the discharge of the fluid in the Λ-space, continuity of flow demands that Q(X, T) be the same at any cross-section of the pipe segment in the Λ-space, depending only upon time
Q(X,T)=Q(T) (46)
The cross-sectional velocity of the flow in the Λ-space may be computed by:
Due to the fluid’s internal friction, shear stresses τ are developed, between the fluid and the inner pipe-wall. Similarly to fractional viscoelasticity, the interface friction shear stress in circular pipes is equal to:
Where μf is fluid viscosity. The friction stresses may be substituted by a fictitious body force per unit length of the pipe in the Λ-fractional space,
The conservation of linear momentum equation in the Λ-fractional space yields:
or
Since, according to Poiseille’s law, F is proportional to the mean flow velocity
Furthermore, the coefficient c of viscous friction is proportional to the fluid viscosity and inversely proportional to the square of the pipe radius:
Thus we obtain the differential equation:
Considering:
The governing Eq.(54) becomes:
For the case that J(t)=J (constant) between the segments x=0 and x=l, the fluid velocity is defined by:
Solution of Fractional D.E:
In case Δp is the pressure drop from one end to the other of the pipe segment with:
Then the total discharge is defined by:
Eq.(58b) is the famous Darcy’s law in the Λ-fractional space. Transferring the flow rate from the Λ-fractional space to the initial space, considering the action of two fractional variables, with fractional order γ1 of time t and γ2 of space x, the flow rate along the pipe is defined by,
Fractional flow in porous media
Considering the difference ΔP to the hydrostatic pressure we get:
Then, we may assume that for the wetted area Av (The area of the voids), Darcy’s law may be applied with:
Assuming that the surface porosity ΦΑ is equal to volume porosity Φ we get:
And thus:
For the specific fluid discharge:
we get:
With
Finally, Darcy’s law may be generalized by:
with P: The pore fluid pressure.
QA: The specific fluid discharge.
K: The permeability of the porous medium.
μf: The viscosity of the fluid.
Transferring the result into the initial space with fractional order γ1 of time t and γ2 of space x, the specific flow is defined by,
Fractional flow in elastic tubes
Consider a thin elastic tube, at its (unstressed) reference placement in the Λ-fractional space, with its inner radius R and its thickness Δ<E. Loading the tube with pressure P, its radius increases by ΔR. Hence its cross-sectional area becomes:
and its current radius
is given by:
Restricting to the linear elasticity where the changes of the radius are infinitesimal, we get:
Therefore the relation between the variable cross-section area and the fluid pressure and cross-section area:
Since mass balance is expressed by:
Recalling Darcy’s law with:
or
Eliminating the discharge Q, we end up to:
Since:
Consequently, Eq. (75) yields:
However, for small changes in pressure, the linear problem is considered, resulting in:
(78)
where:
Solution of parabolic equation
Following just the same procedure, but for the fractional time derivative fields, we get the fractional parabolic equation,
Furthermore, the initial condition expresses the constant pressure value,
Further, at the time = 0+, the pressure at the entry point of the tube is increased by ΔP, keeping constant the pressure value at the exit point of the tube. Consequently, the boundary conditions for T>0 become,
Nondimensionaling the variables we get:
Omitting the upper stars the nondimensional governing equation of the fractional flow becomes,
with the initial condition,
and the boundary conditions,
The initial condition for the pressure,
At T = 0+, the pressure at the entry point is increased by Δp with constant pressure at the exit point. Hence the boundary conditions for T >0 become:
Non-dimensional variables become:
The non-dimensional equation becomes:
With the initial conditions.
And boundary conditions:
Considering first the steady solution:
and introducing a renormalized pressure:
from the above formulas occur:
The initial condition is:
and the boundary conditions:
Applying the separation of variables technique for the fractional diffusion equation we get:
Following the conventional procedure for the diffusion equation, we get:
Due to the appendix, the solution to the Εq. (100) is given by:
Where:
For H=X-1, (Appendix)
Then, the results may be transferred into the initial space, through the transformations,
On Λ-fractional geodesics with corners
The presented analysis on fractional hydrodynamics assumes locally stable fields. Nevertheless, Λ-fractional analysis is inherently global, and consequently, non-smooth fields should be accepted. Continuous fields with non-smooth derivatives may be considered in various fields. yielding smooth geodesics (fields), Abraham & Marsden. Nevertheless, continuous fields with possible corner conditions are acceptable, since only globally stable fields with possible non-smooth geodesics are allowed in the Λ-fractional continuum mechanics. The various variational procedures may be globally considered with the consideration of the Erdmann-Weierstrass conditions, Gelfand & Fomin [42-44].
Proceeding to the analysis presented in the preceding paragraph, the balance laws, are described in the general form,
where Α is the surface inside the material body with the corners of the geodesics and V(n) is the normal velocity of the singular surface. Following Chadwick [45-55], the general jump condition is expressed by,
where
The basic jump conditions, corresponding to the basic equation concerning the mass, linear momentum, and energy conservation are expressed by,
[ρV]=0, (109)
[ρVv+t(n)]=0, (110)
These equations define the shocks in the continuous media. In the fractional analysis, those equations should be satisfied into the Λ-fractional space and the results should be transferred into the initial space.
Fractional shocks in elastic tubes
Let us consider a long cylinder filled with perfect gas and closed at one end by a plane piston. Initially, the gas is considered at rest in the Λ-fractional space with pressure Po and density Ro. If the piston is moving at a constant speed, U determines the speed of propagation of the shock wave in the initial space along with the distributions of the gas density and pressure. Let us point out that a perfect gas is a compressible ideal fluid the pressure being proportional to Rg where g(>1) is a constant.
The shock wave front in the cylinder separates it into two parts, the part with the distributed gas and the part with stationary gas. With the use of Eq.(108), the basic jumping conditions (109-111) are applied to the Λ-fractional space, and if the velocity of the gas behind the piston is U, the mass balance law is expressed by,
R(Vn-U)=RoVn. (112)
Further, the linear momentum balance law is expressed by,
Po-P=-RoVnU , (113)
and the energy balance,
Vn is the local speed of propagation of the shock wave, and P and R are the uniform values of the pressure and density to the rear of the shock. Following Chadwick [45] the shock motion in the Λ-fractional space is defined by the equation
with c0=(gPo/Ro)1/2.
(115)
Finally, the transformation transfers the local shock speed to the initial space.
Conclusion
Since the preliminary elements have been defined (Leibniz Fractional Derivative, Fractional Gradient, Fractional Rotation, Fractional Divergence, Fractional Circulation, and Fractional Gauss’ Theorem), the basic equations of fluid mechanics are reinstated and analyzed. Further, Fractional Darcy’s flow was studied as an application of the fractional flows into porous media. Further, the globally stable flows generate shocks which are described in the context of Λ-fractional analysis. The main issue in our case is the experimental validation of the occurring equations.
Appendix