The Use of Fractional Order Derivative to Predict the Groundwater Flow

The aim of this work was to convert the Thiem and the Theis groundwater flow equation to the time-fractional groundwater flow model. We first derived the analytical solution of theTheim time-fractional groundwater flow equation in terms of the generalized Wright function. We presented some properties of the Laplace-Carson transform. We derived the analytical solution of the Theistime-fractional groundwater flow equation (TFGFE) via the Laplace-Carson transform method. We introduced the generalized exponential integral, as solution of the TFGFE.This solution is in perfect agreement with the data observed from the pumping test performed by the Institute for Groundwater Study on one of its borehole settled on the test site of the University of the Free State. The test consisted of the pumping of the borehole at the constant discharge rate Q and monitoring the piezometric head for 350 minutes.


Introduction
Groundwater problem is perhaps one of the most difficult real-world problems to be modelled into mathematical formulation.To model this problem accurately, one must know precisely the behavior of the medium through which the water is moving.However, this medium through which the flow occurs can change from one point to another, also from one period to another.For example, the hydraulic conductivity of an aquifer can differ from one direction to another.Several scholars have intensively tried to propose a better model that can be used to predict the movement of water through the aquifer.However, their results still present some lacks.Recently, It was revealed that real problems modelled via fractional order derivative present better results when matching their mathematical representation with experimental data.To test this, Botha and Cloot [1] presented some good results by generalizing the groundwater flow equation to the concept of fractional order derivatives.In the same line of ideas, Atangana [2] examined an approximate solution of the generalized groundwater flow equation via the Frobenius method.The results obtained from his investigation showed better prediction.Recently, Atangana and Botha further extended the fractional groundwater equation to the concept of the fractional-variation order groundwater [3].They presented the stability and the convergence of the numerical scheme via Crank-Nicolson method.Up to now, there is no approximate or exact analytical mathematical expression that can be used to describe the solution of the fractional groundwater flow equation.Therefore, one of the purposes of this work is to present some analytical mathematical expression than can be used as approximate solution of the time-fractional groundwater flow equation.
An aquifer test (or a pumping test) is conducted to evaluate an aquifer by "stimulating" the aquifer through constant pumping and observing the aquifer's "response" (drawdown) in observation wells.Aquifer testing is a common tool that hydrogeologists use to characterize a system of aquifers, aquitards, and flow system boundaries.Aquifer tests are typically interpreted by using an analytical model of aquifer flow (the most fundamental being the Theis solution) to match the data observed in the real-world then assuming that the parameters from the idealized model apply to the realworld aquifer.In more complex cases, a numerical model may be used to analyze the results of an aquifer test, but adding complexity does not ensure better results.For the most part frequently, an aquifer assessment is carried out by propelling water out from one borehole at a fixed speed and for as a minimum of 24 hours at the same time as cautiously evaluating the water levels in the observed borehole.When water is pumped from the pumping well, the pressure in the aquifer that feeds that well declines.This decline in pressure will show up as drawdown (change in hydraulic head) in an observation well.Drawdown decreases with radial distance from the pumping well and drawdown increases with the length of time that the pumping continues.The aquifer characteristics which are evaluated by most aquifer tests are [4,5] as follows.
(i) The hydraulic conductivity is the rate of flow of water through a unit crosses sectional area of an aquifer at a unit hydraulic gradient.In English units the rate of flow is in gallons per day per square foot of cross sectional area.
(ii) Specific storage or storativity being a measure of the amount of water a confined aquifer will give up for a certain change in head.
(iii) The transmissivity is the rate at which water is transmitted through a unit thickness of an aquifer under a unit hydraulic gradient.It is equal to the hydraulic conductivity times the thickness of an aquifer.
The rest of this paper has been presented as follows.In Section 1, we presented the background of the fractional order derivative.We derived the analytical solution of the Theim fractional groundwater flow equation in Section 3.
In Section 4, we presented the derivation of the Theis fractional groundwater flow via the Laplace-Carson transform method.We presented an alternative analytical solution of the Theis fractional groundwater flow equation in terms of the generalized exponential integral in Section 5.The numerical comparisons with experimental data are presented in Section 6 and the conclusion is in Section 7. We will start with the background of the fractional derivative.

Background of the Fractional Order Derivative
There exists the vast literature on different definitions of fractional derivatives.The most popular ones are the Riemann-Liouville and the Caputo derivatives [6][7][8][9][10][11]. Caputo's definition has the form of For the case of the Caputo fractional order derivative, we have the following definition: Each of the pervious fractional order derivatives presents some advantages and disadvantages [6][7][8][9][10].The Riemann-Liouville derivative of a constant is not zero while Caputo derivative of a constant is zero but demands higher conditions of regularity for differentiability [6][7][8][9][10]: to compute the fractional derivative of a function in the Caputo sense, we must first calculate its derivative [11].Caputo derivatives are defined only for differentiable functions while functions that have no first-order derivative might have fractional derivatives of all orders less than one in the Riemann-Liouville sense [12,13].Guy Jumarie has recently modified the Riemann-Liouville derivative (see [14]) The Caputo fractional derivative will be considered in this work due to the applicability of the Caputo derivative in realworld problems [15].

Thiem' Groundwater Flow Equation
The Theim fractional groundwater flow equation is an ordinary differential equation given in the following.The equation describes the change in level of water as function of distance during the pumping test [5]: Subject to the initial condition,  = 2  (Φ(  )).
We will make use of the Laplace transform to derive analytical solution of (4).Thus, multiplying on both sides of (4) by  and secondly applying the Laplace transform, we obtain the following expression: where   =  − 0 + Φ(0 + ) ( = 2, . . ., ).Now, one can derive the solution of the ordinary order differential equation with respect to the Laplace transform of Θ() = (Φ()): with  1 an arbitrary real constant that will be obtained via the initial condition.We next expand the exponential function in the integrand in a series, and using term-by-term integration, we arrive at the following expression: with, of course, Now, applying the inverse Laplace transform on Θ 1 () and using the fact that we obtain with Ψ 1 [⋅] the generalized Wright function [16] for  = 1 and  = 2.We next expand the exponential function exp[− 1− /(1 − )] in power series; multiplying the resulting two series; in addition of this if we consider the number   (, ) defined for  > 0,  = 2, . . ., ( ̸ = (+−1)/,  ∉ N) and  ∈ N 0 , The previous family of number possesses satisfies the following recursive formula: which produces the explicit expression for   (, ) in the form of Now, having the previous expression on hand, we can derive that However, remembering (9) with  = ( − 1) +  + 1 − , we can further derive the following expression for Φ *  () as or in the simplified version we have where It follows that the solution of the fractional-Thiem groundwater flow equation is in the form of In our case  = 2.

Time-Fractional Theis Groundwater Flow Equation
The easiest sweeping statement of subsurface water flow equation, which while we are on the subject is in addition in harmony in the midst of the real physics of the observed fact, is to presume that water level is not in a balanced state but momentary state.Theis (1935) [17] was the first to develop a formula for unsteady-state flow that introduces the time factor and the storativity.He noted that when a well-penetrating extensive confined aquifer is pumped at a constant rate, the influence of the discharge extends outward with time.The rate of decline of head, multiplied by the storativity and summed over the area of influence, equals the discharge.The unsteady-state (or Theis) equation, which was derived from the analogy between the flow of groundwater and the conduction of heat, is perhaps the most widely used partial differential equation in groundwater investigations: The aforementioned equation is classified under parabolic equation.To include explicitly the variability of the medium through which the flow takes place, the standard version of the partial derivative respect to time is replaced here with time-fractional order derivative to obtain with initial condition Φ(, 0) = 0 and boundary condition lim  → ∞ Φ(, ) = 0,  = 2Φ(  , ), here  is the transmissivity of the aquifer,   is the ratio of the borehole, and  is the discharge rate, or the rate at which the water is being taken out of the aquifer.Some few integral transform operators have been intensively used to solve some kind of ordinary and partial differential equations.See, for instance, the Fourier transform, the Laplace transform, the Mellin transform [18], and the Sumudu transform [19][20][21][22].Beside these integral operators, there exists a similar operators called the Laplace-Carson transform [23]; this operator has been neglected.However, this operator has some properties that can be used to solve a kind of ODE, PDE, FODE, and FPDE.The aim of this section is therefore devoted to the discussion underpinning the definition, properties of the Laplace-Carson transform, and its application to the fractional groundwater flow equation.We shall start with the definition and properties.Definition 1.Let () be a continuous function over an open interval (0, ∞) such that its Laplace transform is  time differentiable; then the Laplace-Carson transform  is defined as follows: and the inverse Laplace-Carson transform is defined as where () is the Laplace transform of ().Before we continue, we shall prove that the previous definition is indeed the inverse inverse Laplace-Carson.In fact, from the definition of inverse Laplace-Carson of a function (), we have that thus, It follows that Therefore, the inverse inverse Laplace-Carson is well defined.

Some Properties of Laplace-Carson Transform
In this part of the section, we consider some of the properties of the inverse Laplace-Carson that will enable us to find further transform pairs {(),   ()} without having to compute and consider the following: Let us verify the previous properties.We shall start with (I), by definition, we have the following: and then the first property is verified.
For (II) we have the following by definition: Now, using the property of the Laplace transform [()]() = (1/)(/), we can further obtain and then, the property number (II) is verified.
For number (III), we have the following: Let () = (); then (38) Now, using the property of Laplace transform of the convolution, we obtain the following: and then, using the property of the derivative for the product of two functions, we obtain And then, the property number (VI) is verified.
For number (VII), by definition, we have the following: Now, using the property of the Laplace transform, Now, deriving the previous expression n-time, we obtain the following expression: This completes the proof of number (VI).We shall now use some properties of the Laplace-Carson transform to solve the fractional groundwater flow equation.To achieve this, we shall start by assuming that the solution of the main equation can be separated as follows: then, the separated equations become where  is the separation constant.The first equation (45) can be solve directly by applying on both sides the Laplace transform to obtain Using the inverse formula of Laplace transform of twoparameter Mittag-Leffler function, we get where the Mittag-Leffler function is defined as follows: To solve the second equation, we make use of the Laplace-Carson presented earlier, to obtain Deriving, we obtain the following ordinary differential equation: for which the exact solution is given as Now applying the inverse Laplace transform operator on both sides in the previous equation, we obtain the following in terms of the Bessel function first kind [24]: where the Bessel function first kind is defined as Therefore, the solution of the fractional groundwater flow equation is given as Making use of the initial and boundary conditions, we obtain the constant  to be Then, We shall present an alternative solution in the next section via the Boltzmann variable method.By using the boundary condition, we can determine the Eigen value of (55).

An Alternative Derivation of the Time-Fractional Groundwater Equation
In this section, we present an alternative approximate solution of the groundwater flow equation.A method frequently used to derive some kind of parabolic partial differential equations is the so-called Boltzmann transformation [5], defined for an arbitrary  0 <  by equation.
Let us consider now the following function: with  an arbitrary constant.If we assume that   is the ratio of the borehole from which the water is taken out of the aquifer, then the total volume of the water withdrawn from the aquifer is given by Hence, is the drawdown that will be observed at a distance, , from the pumped borehole after the period Δ 0 .Now suppose that the a previous procedure is repeated  times; that is, water is withdrawn for a short period of time, Δ  , at a consecutive times,  +1 =   + Δ  , ( = 0, 1, . . ., ).Now, since the fractional groundwater flow equation is linear differential equation, the total drawdown at any time  >   will be given by Therefore, if Δ → 0, the definition of the defined integral can be invoked to write ] .(63) However, using the Boltzman variable, we arrive at the following expression: The previous solution will be called the most generalized general solution of fractional groundwater equation.However, this solution can be simplified somewhat under certain conditions.A particularly important solution which arises when  0 is taken at zero and () is a constant independent of time, and then we arrive at Here will be called the generalized exponential integral.It is worth pointing out that if alpha is equal to 1, we recover the exact analytical solution of the groundwater flow equation proposed by Theis.Alternative iterations method [8,11,20] can be used to derive approximate solutions of these problems.

Numerical Simulations
In this section, we investigate the behavior of the analytical solutions of the Theim fractional and the Theis fractional groundwater flow equation.We compare the analytical solution with the experimental data from the pumping test obtained from the experimental site of the Institute for Groundwater Studies, the University of the Free State, Bloemfontein Campus, South Africa.We shall start with the simulation of Theim fractional groundwater flow equation.The analytical solution of the main problem was depicted in Figures 1(a), 1(b), 1(c), and 1(d).It is worth to mention from the figures that the order of the derivative plays an important role in the simulation.We shall present the numerical solution for the Theistime-fractional groundwater flow equation for a fixed distance.The analytical solution of the main problem was depicted in Figures 2(a We shall present in Figure 3  A Chinese proverb says, I quote "An image is equivalent to ten thousand words." With any reservation we can conclude that the fractional order derivative plays an important role while modeling real-world problem into mathematical equation.This solution is in perfect agreement with the data observed from the pumping test performed by the institute for groundwater study on one of their borehole settled on the test site of the University of the Free State.

Conclusion
We have generalized to Theim and Theis groundwater flow equation to the concept of fractional order derivatives.The resulting equations were solved analytically via several techniques, including the Laplace transform method, the Fourier transform method, the Laplace-Carson transform method, and the Boltzmann variable method.The numerical simulations show that the fractional order derivative plays an important role in the simulation process.In addition, we compare the analytical solution with experimental data to access the accuracy of the fractional groundwater model.The analytical solution was in perfect agreement with experimental data.

Figure 2 :
Figure 2: Numerical representation of (66) for different values of alpha.

1 Figure 3 :
Figure 3: Comparison of Theis-time-fractional groundwater flow equation with experimental data from real observation.