Mathematical Description of the Groundwater Flow and that of the Impurity Spread, which Use Temporal Caputo or Riemann–Liouville Fractional Partial Derivatives, is Non ‐ Objective

: In this paper, it is shown that the mathematical description of the bulk fluid flow and that of content impurity spread, which uses temporal Caputo or temporal Riemann–Liouville fractional order partial derivatives, having integral representation on a finite interval, in the case of a horizontal unconfined aquifer is non ‐ objective. The basic idea is that different observers using this type of description obtain different results which cannot be reconciled, in other words, transformed into each other using only formulas that link the numbers representing a moment in time for two different choices from the origin of time measurement. This is not an academic curiosity; it is rather a problem to find which one of the obtained results is correct.


Introduction
The mathematical description of a real-world phenomenon is objective if it is independent of the observer. That is, it is possible to reconcile the observation of a phenomenon into a single coherent description of it. This requirement was pointed out by Galileo Galilee (1564-1642), Isaac Newton (1643-1727), and Albert Einstein  in the context of mathematical description of a mechanical movement: "The mechanical event is independent of the observer". A possible and elementary understanding of the independence of the mechanical event from the observer is the independence of the event from the choice of the frame of reference and from the choice of the moment considered the origin for time measurement, made by an observer. What this means precisely in this paper is presented in the following. To describe mathematically the evolution of a mechanical event, an observer chooses a fixed orthogonal reference frame in the affine Euclidean space, a fixed moment in time (called origin for time measurement), and a unit for time measuring (second). For different observers this choice can be different. In this paper the "objectivity of a mathematical description" means that the description is independent of the choice of the fixed orthogonal reference frame and of the choice of origin for time measuring. This means that the results obtained by two different observers can be reconciled, in other words, transformed into each other, using only formulas that link the coordinates of a point in two fixed orthogonal reference frames and formulas that link the numbers representing a moment in time for two different choices of the origin of time measurement. This concept of "objectivity of a mathematical description" is different from the concept of "objectivity in physics" presented in [1]. The advantage of our kind of understanding of the "objectivity of a mathematical description" used in this paper, is that it is less restrictive than Galilean invariance, Lorentz invariance or Einstein covariance and General covariance, it can be easily applied in a specific case and the reader does not need prior knowledge of either Galilean or Lorentz invariance, or Einstein or General covariance. Mathematical descriptions which depend on the choice of the fixed orthogonal reference frame or on the choice of the origin of time measurement are non-objective in the sense of this paper. In the case of descriptions which are non-objective, two observers who describe the same mechanical event obtain two different results that cannot be reconciled, in other words, cannot be transformed into each other, using only formulas that link the coordinates of a point in two fixed orthogonal reference frames and formulas that link the numbers representing the same moment in time for two different choices of the origin of time measurement. This concept of non-objective description can be easily applied in a specific case and the reader does not need prior knowledge of either Galilean or Lorentz invariance, or Einstein or General covariance. The majority of mathematical descriptions formulated in terms of integer-order derivatives or integer-order partial derivatives, reported in the literature (books of differential equations of mathematical physics), are objective in the sense of this manuscript.
In the following, the objectivity of the description formulated in terms of integer-order derivatives of some phenomena appearing in hydrology is illustrated. The details related to the verification of objectivity in the case of the mathematical description of the groundwater flow to the well and the spreading of impurities were introduced in the manuscript for several reasons: In the literature, accessible to us for free, we have not found the justification of the objectivity of the mathematical description that uses integer-order derivatives.
We wanted to show how the objectivity of a mathematical description in terms of integer-order derivatives can be verified directly by a simple case occurring in nature.
In these details there are formulas, which represent partial results, that are also used in Sections (2)(3)(4)(5) in which the non-objectivity is discussed in the case of the use of temporal Caputo or Riemann-Liouville fractional derivatives having integral representation on a finite interval.
The reasoning that is done to prove the objectivity of the mathematical description in the case of integer-order derivatives, helps in understanding the reasoning that is done to demonstrate the nonobjectivity of the description in the case of using temporal Caputo or Riemann-Liouville fractional derivatives, having integral representation on a finite interval. In the case of using fractional derivatives, the reasoning is by ʺreduction ad absurdumʺ. At the beginning it is assumed that the description is objective. After this, the steps are followed when demonstrating objectivity in the case of using integer-order derivatives. So an equality is obtained, which follows from the assumption of the objectivity. In general, the obtained equality is not true. Hence the conclusion that the hypothesis of the objectivity of the description (in this case of the use of fractional derivatives, having integral representation on a finite interval), is false.
In hydrology [2] and [3] the horizontal unconfined aquifer around the well is represented as a subset  of points of the affine Euclidean space 3 E . A vertical section of  is shown in Figure 1. A particle Q of the bulk fluid at a moment in time M is represented by a point P of . P is the place occupied by the fluid particle Q at the moment in time M . To describe the position of P , observer O chooses a fixed orthogonal reference frame ) , , ; represents the associated velocity of the fluid particle Q , that at the moment in time M t is in the point P of coordinates ) , , ( E . Therefore, for the coordinates the following relations hold: The significance of the quantities appearing in the above mentioned relations are: represent the same piezometric head. Therefore, the following relations hold: Equations (1)-(4) and (6) . This means that the above presented mathematical description of the bulk groundwater flow to the well, in a horizontal unconfined aquifer, is objective.
In classical theory of the 2D flow to the well, in a horizontal unconfined aquifer, [2] and [3], the real valued function that describe the 2D flow in terms of the observer O, verify the equations: where S is the storage coefficient; T is the transmissivity of the aquifer; S Q is the leakage rate; K is the hydraulic conductivity and  is the porosity.
The real valued function , that in terms of the observer * O , describe the bulk fluid 2D flow to the well, verify the equations: where S is the storage coefficient ; T is the transmissivity of the aquifer; S Q is the leakage rate; K is the hydraulic conductivity and  is the porosity.
A short proof of the objectivity of this description is the following: assume that the functions (10) and (11) and consider the functions (14). Note that the following equalities hold: Using Equations (16) and replacing the terms in (12) and (13), follows that the functions (14), verify Equations (12) and (13). The second part of the proof is similar.
Therefore, the description with Equations (10) and (11) is objective. That is, different observers, describing the groundwater flows with these tools, obtain results that can be reconciled, in other words, transformed into each other, using Equations (1)-(4) that link the coordinates of a point in two fixed orthogonal reference frames and formulas that link the numbers representing a moment in time for two different choices of the origin of time measurement.
The spread of an impurity, contained in the bulk fluid flowing in a porous media, is described by the concentration of that impurity. Observer O describes the concentration with the real valued function ) , , , , that verifies the partial differential equation (PDE):  In PDE (17) ) , , ( where ) , , ( *) *, *, ( In case of the 2D flow to a well, in a horizontal unconfined aquifer, for observer O , ) , , and PDE (17) becomes: Under the same hypothesis, for observer and PDE (18) becomes: Objectivity of the impurity spread description means that the solutions of PDEs (19) and (20) verifies Equation (20) verifying Equation (19).
In the following we provide a short proof of the objectivity of this description in the case, when the aquifer is homogeneous, isotropic, and D and S are constants.
In this case Equations (19) and (20) become: , defined by (21). For this function the following equalities hold:   (27) Replacing in (20)  given by (22), then, in a similar way, we obtain that this function is a solution of Equation (19).
So, the description of the impurity spread with (23) is objective. That is, different observers, describing the impurity spread with these tools, obtain results that can be reconciled, in other words, transformed into each other, using Equations (1)-(4).
The objectivity of the above presented descriptions implies that different observers describing the same phenomenon, using integer-order partial derivatives, obtain results that can be reconciled, in other words, transformed into each other, using Equations (1)- (4).
Beside the objective mathematical descriptions of the bulk fluid flow and impurity spread in porous media (see references [1][2][3][4][5][6][7]) formulated in terms of integer-order partial derivatives, there are mathematical descriptions of the bulk fluid flow and impurity spread in porous media which use fractional order temporal or spatial partial derivatives. See, for instance, references [8][9][10][11][12][13][14][15][16][17][18][19][20][21]. In these works, the analysis of the objectivity is missing. At first, we thought that in the case of the description with fractional derivatives, the objectivity is fulfilled and therefore it is ignored. However, curiosity has pushed us to see how the fulfillment of the objectivity condition (in sense of our manuscript) can be proven mathematically. For this special issue, we have chosen fractional behavior in nature for the very simple case of groundwater flow to well and the spread of impurities. Thus were "born" Sections 2 and 3 of the manuscript in which we analyzed the objectivity of the description of the groundwater flow and Sections 4 and 5, in which we analyzed the of objectivity of the description of the spread of impurities, instead of using the integer-order derivatives temporal Caputo or Riemann-Liouville fractional partial derivatives, having integral representation on finite interval. The purpose of the present paper is to show that, in case of an unconfined horizontal aquifer, the mathematical descriptions, that use Caputo or Riemann-Liouville fractional order temporal partial derivatives, having integral representation on a finite interval, generally are non-objective. Two observers, who use fractional derivatives, obtain different results that cannot be reconciled, in other words, transformed into each other, using Equations (1)- (4).
Remember that for a continuously differentiable function , is defined with the following integral representation on a finite interval (see [22]): Note that the derivative defined with (28) was considered by other people before Caputo, like Gherasimov (see [22]). So, the name of Caputo used in this paper may not be appropriate.
For a continuously differentiable function , is defined with the following integral representation on a finite interval (see [22]): In Equations (28)  [ is the integer part of  .

In case of an Unconfined Horizontal Aquifer the Piezometric Head Dynamics Description, Using Temporal Caputo Fractional Order Partial Derivatives, with Integral Representation on a Finite Interval, is Non-Objective
Assume that in the piezometric head dynamics description that the temporal Caputo fractional partial derivative of order  ,   , defined by (14). Note that, in this particular situation, the following equalities hold: Using Equation (32) and replacing the terms in (31), it follows that: if the function (14)) verifies Equation (31), then the following equality holds: Equation (33) is a consequence of the assumption that the mathematical description (30) is objective. However, generally (33) it is not verified. This means that the assumption that the piezometric head dynamics description with Equation (30) (1) and (2), that link the numbers representing the same moment in time for two different choices of the origin of time measurement. The problem is to find which one of the results is correct.

In case of an Unconfined Horizontal Aquifer the Piezometric Head Dynamics Description which Uses Temporal Riemann-Liouville Fractional Order Partial Derivatives, with Integral Representation on a Finite Interval, is Non-Objective
Assume that in the piezometric head dynamics description that the temporal Riemann-Liouville fractional partial derivative of order  , 1 0    , with integral representation on a finite interval, is used. In case of the 2D flow to a well, in a horizontal unconfined aquifer, Equation (10) for observer O and Equation (12) for observer * O become: Objectivity of the piezometric head dynamics description means that the solutions of the fractional partial differential Equations (34)  , defined by (14). Note that, in this particular situation, the following equalities hold: Equation (37) is a consequence of the assumption that the mathematical description (34) is objective. However, generally, (37) is not verified. This means that the assumption that the piezometric head dynamics described by Equation (34) (1) and (2): that link the numbers representing a moment in time for two different choices of the origin of time measurement. The problem is to find which one of the results is correct.

In case of an Unconfined Horizontal Aquifer the Impurity Spread Description which Uses Temporal Caputo fractional Order Partial Derivatives, with Integral Representation on a Finite Interval, is Non-Objective
Assume that in the impurity spread dynamics description the temporal Caputo fractional partial derivative of order  , 1 0    , with integral representation on a finite interval, is used. In case of the 2D flow to a well, in a horizontal unconfined isotropic homogeneous aquifer, Equation (23) Objectivity of the impurity spread description means that the solutions of the fractional partial differential Equations (38)  , defined by (21). Note that, in this particular situation, the following equalities hold: Equation (41) (1) and (2): that link the numbers representing a moment in time for two different choices of the origin of time measurement. The problem is to find which one of the results is correct.

In case of an Unconfined Horizontal Aquifer, the Impurity Spread Description which Uses Temporal Riemann-Liouville Fractional Order Partial Derivatives, with Integral Representation on a Finite Interval, is Non-Objective
Assume that, in the impurity spread dynamics description the temporal Riemann-Liouville fractional partial derivative of order  , fractional order partial derivatives, having integral representation on a finite interval, is necessary. The results obtained by us in Sections 4 and 5 can be instructive for the authors of some of the papers [16][17][18][19][20][21], because they show that the use of temporal Caputo and Riemann-Liouville fractional partial derivatives affect the objectivity of the description of the spread of impurities in porous media. There is an argument for why the analysis of the objectivity of the mathematical description of the spread of impurities in porous media proposed in the papers that use temporal Caputo and Riemann-Liouville fractional order partial derivatives, having integral representation on a finite interval, is necessary. 4. In the early 2000s a discussion started about the initialization problems in [23][24][25][26][27][28]. Some published results in [26] and [27] concluded in the inconsistency of Caputo and Riemann-Liouville's definition to take into account initial conditions if these definitions are used in fractional partial differential equations or in ordinary differential equations. In [23], [26], and [27] a time shift was used to highlight the above mentioned problem. Our approach to the question: why can integer-order derivatives not simply be replaced by fractional-order derivatives to develop the fractional-order theories? is different. What we know from the scientific literature is that the assertion "integer-order derivatives cannot be simply replaced by fractional-order derivatives to develop the fractional-order theories" has not been proven so far. In fact, the general assertion, as formulated, refers to all the equations of mathematical physics and we do not think it will be proven soon. However, what we think is fact is that this statement can be demonstrated in some proper cases. In this paper we actually demonstrate this statement in the case of describing fluid flow in porous media and impurity spread also showing the cause, in other words, that the objectivity of the description is lost. 5. A given mathematical tool is not necessarily appropriate for the mathematical description of a certain real word phenomenon.