Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 543026, 9 pages

http://dx.doi.org/10.1155/2013/543026

## The Use of Fractional Order Derivative to Predict the Groundwater Flow

^{1}Institute for Groundwater Studies, Faculty of Natural and Agricultural Sciences, University of the Free State, P.O. Box 9300, Bloemfontein, South Africa^{2}Department of Mathematics, Faculty of Art & Sciences, Celal Bayar University, Muradiye Campus, 45047 Manisa, Turkey

Received 2 July 2013; Revised 27 August 2013; Accepted 3 September 2013

Academic Editor: Tirivanhu Chinyoka

Copyright © 2013 Abdon Atangana and Necdet Bildik. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

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 the Theim 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 Theis-time-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.

#### 1. 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 real-world 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.

#### 2. 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–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–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–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 real-world problems [15].

#### 3. 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, .

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 . Now, one can derive the solution of the ordinary order differential equation with respect to the Laplace transform of : with 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 and using the fact that we obtain with the generalized Wright function [16] for and . We next expand the exponential function in power series; multiplying the resulting two series; in addition of this if we consider the number defined for , , and , 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 , 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 .

#### 4. 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 and boundary condition , , 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–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 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.

#### 5. 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 , 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 , we can further obtain and then, the property number (II) is verified.

For number (III), we have the following: Let ; then By the theorem of inverse Laplace transform, we obtain Number (IV) and (V) are obvious to be verified. For number (VI), we have the following by definition: 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 two-parameter 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).

#### 6. 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 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 .

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, , . Now, since the fractional groundwater flow equation is linear differential equation, the total drawdown at any time will be given by Therefore, if , the definition of the defined integral can be invoked to write 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 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.

#### 7. 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 Theis-time-fractional groundwater flow equation for a fixed distance. The analytical solution of the main problem was depicted in Figures 2(a), 2(b), 2(c), and 2(d).

We shall present in Figure 3 the comparison of the analytical solution of the Theis-time-fractional ground water flow equation, with the experimental data 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 and monitoring the piezometric head for 350 minutes, at the distance of m.

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.

#### 8. 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.

#### References

- J. F. Botha and A. H. Cloot, “A generalised groundwater flow equation using the concept of non-integer order derivatives,”
*Water SA*, vol. 32, no. 1, pp. 1–7, 2006. View at Google Scholar · View at Scopus - A. Atangana, “Numerical solution of space-time fractional order derivative of groundwater flow equation,” in
*Proceedings of the International Conference of Algebra and Applied Analysis*, p. 20, Istanbul, Turkey, June 2012. - A. Atangana and J. F. Botha, “Generalized groundwater flow equation using the concept of variable order derivative,”
*Boundary Value Problems*, vol. 2013, article 53, 2013. View at Publisher · View at Google Scholar - J. Boonstra and R. A. L. Kselik,
*SATEM 2002: Software for Aquifer Test Evaluation*, International Institute for Land Reclamation and Improvement, Wageningen, The Netherlands, 2002. - G. P. Kruseman and N. A. de Ridder,
*Analysis and Evaluation of Pumping Test Data*, International Institute for Land Reclamation and Improvement, Wageningen, The Netherlands, 2nd edition, 1990. - H. Jafari and C. M. Khalique, “Homotopy perturbation and variational iteration methods for solving fuzzy differential equations,”
*Communications in Fractional Calculus*, vol. 3, no. 1, pp. 38–48, 2012. View at Google Scholar - S. Duan, R. Rach, D. Buleanu, and A. M. Wazwaz, “A review of the Adomian decomposition method and its applications to fractional differential equations,”
*Communications in Fractional Calculus*, vol. 3, no. 2, pp. 73–99, 2012. View at Google Scholar - A. Atangana and E. Alabaraoye, “Solving system of fractional partial differential equations arisen in the model of HIV infection of CD4+ cells and attractor one-dimensional Keller-Segel equation,”
*Advances in Difference Equations*, vol. 2013, article 94, 2013. View at Publisher · View at Google Scholar - D. Baleanu, K. Diethelm, E. Scalas, and J. J. Trujillo,
*Fractional Calculus Models and Numerical Methods*, Series on Complexity, Nonlinearity and Chaos, World Scientific, Singapore, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - A. Atangana and A. Secer, “A note on fractional order derivatives and table of fractional derivatives of some special functions,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 279681, 8 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - G. C. Wu, “New trends in the variational iteration method,”
*Communications in Fractional Calculus*, vol. 2, pp. 59–75, 2011. View at Google Scholar - G. Jumarie, “On the representation of fractional Brownian motion as an integral with respect to (dt)a,”
*Applied Mathematics Letters*, vol. 18, no. 7, pp. 739–748, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - G. Jumarie, “Modified Riemann-Liouville derivative and fractional Taylor series of nondifferentiable functions further results,”
*Computers and Mathematics with Applications*, vol. 51, no. 9-10, pp. 1367–1376, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - G. Jumarie, “Table of some basic fractional calculus formulae derived from a modified Riemann-Liouville derivative for non-differentiable functions,”
*Applied Mathematics Letters*, vol. 22, no. 3, pp. 378–385, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - A. Atangana, “New class of boundary value problems,”
*Information Sciences Letters*, vol. 1, no. 2, pp. 67–76, 2012. View at Google Scholar - S. G. Samko, A. A. Kilbas, and O. I. Marichev,
*Fractional Integrals and Derivatives, Translated from the 1987 Russian Original*, Gordon and Breach, Yverdon, Switzerland, 1993. View at MathSciNet - C. V. Theis, “The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using ground-water storage,”
*Transactions, American Geophysical Union*, vol. 16, no. 2, pp. 519–524, 1935. View at Google Scholar - P. Flajolet, X. Gourdon, and P. Dumas, “Meilin transforms and asymptotics: harmonic sums,”
*Theoretical Computer Science*, vol. 144, no. 1-2, pp. 3–58, 1995. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - G. K. Watugala, “Sumudu transform: a new integral transform to solve differential equations and control engineering problems,”
*International Journal of Mathematical Education in Science and Technology*, vol. 24, pp. 35–43, 1993. View at Google Scholar - A. Atangana and A. Kılıçman, “The use of Sumudu transform for solving certain nonlinear fractional heat-like equations,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 737481, 12 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - M. G. M. Hussain and F. B. M. Belgacem, “Transient solutions of Maxwell's equations based on sumudu transform,”
*Progress in Electromagnetics Research*, vol. 74, pp. 273–289, 2007. View at Google Scholar · View at Scopus - S. Weerakoon, “The “Sumudu transform” and the Laplace transform—reply,”
*International Journal of Mathematical Education in Science and Technology*, vol. 28, no. 1, pp. 159–160, 1997. View at Google Scholar - F. Oberhettinger and L. Badii,
*Tables of Laplace Transforms*, Springer, Berlin, Germany, 1973. View at MathSciNet - I. Podlubny,
*Fractional Differential Equations*, Academic Press, New York, NY, USA, 1999. View at MathSciNet