Abstract and Applied Analysis

Volume 2014, Article ID 498381, 7 pages

http://dx.doi.org/10.1155/2014/498381

## A Modified Groundwater Flow Model Using the Space Time Riemann-Liouville Fractional Derivatives Approximation

^{1}Institute for Groundwater Studies, Faculty of Natural and Agricultural Sciences, University of the Free State, Bloemfontein 9300, South Africa^{2}Department of Mathematical Sciences, North-West University, Mafikeng Campus, Mmabatho 2735, South Africa

Received 16 December 2013; Accepted 10 February 2014; Published 8 May 2014

Academic Editor: Dumitru Baleanu

Copyright © 2014 Abdon Atangana and S. C. Oukouomi Noutchie. 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 notion of uncertainty in groundwater hydrology is of great importance as it is known to result in misleading output when neglected or not properly accounted for. In this paper we examine this effect in groundwater flow models. To achieve this, we first introduce the uncertainties functions as function of time and space. The function accounts for the lack of knowledge or variability of the geological formations in which flow occur (aquifer) in time and space. We next make use of Riemann-Liouville fractional derivatives that were introduced by Kobelev and Romano in 2000 and its approximation to modify the standard version of groundwater flow equation. Some properties of the modified Riemann-Liouville fractional derivative approximation are presented. The classical model for groundwater flow, in the case of density-independent flow in a uniform homogeneous aquifer is reformulated by replacing the classical derivative by the Riemann-Liouville fractional derivatives approximations. The modified equation is solved via the technique of green function and the variational iteration method.

#### 1. Introduction

This paper investigates the effects of uncertainty on the predictive accuracy of flow through porous media; it is commonly believed that the problem that occurs in groundwater models is the suitable geometry in which flow occurs on one hand and the deviation of theoretical expected values from observations on the other hand. Therefore, it is important to notice that miniscule effects observed always require the most new modifications of ideas. Scientists in the field of hydrogeology in particular are used to deal with doubt and uncertainty, because it is impossible to understand or to model the phenomena that occur in aquifers exactly. All historical and current theoretical knowledge in groundwater investigations are uncertain and doubtful. This experience with doubt and uncertainty is important. We believe that it is of great importance, and one that extends beyond the theories which are used to interpret the phenomena that take place in aquifers. Doubt is clearly a value that must be analytically included in groundwater flow models. Uncertainty in groundwater hydrology originates from different sources. Neglecting uncertainty in groundwater assessments can lead to incorrect results and misleading output. Generally there are various sources of uncertainty in model outputs, for example, uncertainty associated with lack of knowledge or accuracy of the model inputs as well as the structural uncertainty related to the mathematical interpretation of the model. The assessment and presentation of the effects of uncertainty are now widely recognized as important parts of analyses for complex systems [1–6]. At the simplest level such analyses can be viewed as the study of functions. In order to include explicitly the possible effect of the uncertainties into mathematical models, we introduced in this paper the uncertainties in groundwater models as a function of time and space.

Consider

#### 2. Modification of Groundwater Flow Equation

To be clear, the modification of the classical model for groundwater flow in the case of density independent flow in the uniform and homogeneous aquifer is considered in this paper.

To modify this, we make use of Riemann-Liouville fractional derivatives that were introduced in [7] and attempted by many others, see for example, [8]. These derivatives are defined as Here, is the Euler gamma function; +1, where is the integer part of for that is and for . Following equation (2) we have that and . The integral operator defined above for fractional exponents and depending on coordinates and time can be expressed in terms of ordinary derivative and integral [7] for . For this matter, generalized Riemann-Liouville fractional derivatives satisfy the approximate relations.

Consider The above relations make it possible to describe the flow system, including the effect of uncertainties on the behaviour of physical systems, by means of partial differential and integral equations.

Let us examine some properties of the above derivative operator [9].(i)Addition. If , and are differentiable in the opened interval** I,** then
(ii)Division. If and are differentiable on the opened interval** I,** then
(iii)Multiplication. If , and are differentiable in the opened interval** I,** then
(iv)Power. If , and are differentiable in the opened interval** I,** then
(v)If , and are two times differentiable in the opened interval** I,** then
It is important to observe that if , we recover the properties of normal derivatives. Recent investigations suggest that the flow is influenced by the geometry of the bedding parallel factures. An attempt to circumvent this problem, Barker introduced a model in which the geometry of the aquifer is regarded as a fractal [10]. In the same direction, the authors in [11] introduced the concept of noninteger fractional derivative to investigate a radially symmetric form of (1); by replacing the classical first order derivative of the piezometric head by a complementary fractional derivative as results of their investigation, they found that there is a close relationship between the fractal and the fractional. Therefore, to include the fractal dimension into the mathematical formulation of the modified groundwater flow equation we next introduce the constant fractal dimension . The classical model for groundwater flow, in the case of density-independent flow in a uniform homogeneous aquifer, can then be reformulated as follows:
where is the hydraulic conductivity of the aquifer, is the specific storativity of the aquifer, is the strength of any sources or sink, here it will be neglected, and finally is the piezometric head.

In order to meet the physical and mathematical requirements we impose the uncertainties function to be a positive function such that Equation (9) makes it possible to describe the flow through the geological formation, and the effect of uncertainties on the behaviour of physical systems, by means of partial differential and integral equations. However there is no analytical solution for this equation, in fact the analytical solution is very difficult to determine. Therefore we need the following approximation to simplify (9): Making use of (9), (3), and (11) we obtain the following equation: Since uncertainties additions to unit are small, the right- and left-hand sides of (12) can be divided by () to obtain the following approximate equation: where term is omitted because it is significantly very small and since we are dealing with approximation here, it needs not to be considered in this case. For simplicity (13a) can be reformulated as

Here the additional term can be roughly approximate to and satisfies the equation of classical model for groundwater flow in the case of density independent flow in the uniform and homogeneous aquifer. It is important to observe that the modified equations (13a) and (13b) differ from the standard form of groundwater flow equation in three properties:

There is a new operator that takes into account the variation in piezometric head and uncertainties function given below as Second, the “force” appears due to the coordinate dependence of uncertainty function. And finally, there is a derivative-free term that depends only on the uncertainties time function and is proportional to the piezometric head and characterises, depending on the coefficient sign, the retardation or enhancement of the flow through the porous media.

It is important to point out that those terms in (13a) that involve fractional additions, and , to the time and space dimensions are small. It follows that this equation can be solved approximately by changing the function by , which satisfies the standard version of the groundwater flow equation which is the left-side of (13a) and (13b), in terms concerning . Now let us suppose that such change is made in the expression (14). Equations (13a) and (13b) become Before solving the above equation, one needs to relate the additional function in the modified equation to physical situation that takes place in the aquifers. Some deterministic models treat the properties of porous media as lumped parameters (essentially, as a black box), but this prevents the representation of heterogeneous hydraulic properties in the model. Heterogeneity or variability in aquifer properties is characteristic of all geologic systems and is now recognised as playing a key role in influencing groundwater flow and solute transport. Thus, it is often preferable to apply distributed-parameter models, which allow the representation of more realistic distributions of system properties.

The lithology of most geological formations tends to vary significantly, both horizontally and vertically. Consequently, geological formations are seldom homogeneous. Figure 1 is an example of layered heterogeneity.

Heterogeneity occurs not only in the way shown in the Figure 1, however, individual layers may pinch out; their grain size may vary in horizontal direction, they may contain lenses of other grain sizes, or they may be discontinuous by faulting or scour-and-fill structures.

The distribution of sedimentary facies controls the heterogeneity of hydrogeological properties of porous sedimentary aquifers at different scales. The arrangement of individual facies and their porosity and permeability determine the path of groundwater flow across sedimentary bodies. Therefore the capability to forecast hydrogeological heterogeneity due to facies changes helps to improve solutions of flow and diffusion problems in this kind of aquifer. When real aquifers are studied, it is impossible to model groundwater flow at a scale such that we can take into account the effects of fine-scale sedimentary heterogeneity; in fact this would require a precise knowledge of the sedimentary bodies that cannot be obtained from sparse data at some wells and this would be prohibitive for the required computing power. Therefore the fine scale heterogeneity is usually “up-scaled” and the heterogeneous real medium is substituted at a larger scale with an equivalent often anisotropic medium, whose parameters allow the reproduction of the average flow of the real heterogeneous sedimentary structure. In this paper the function will be considered to account for the effect of heterogeneity and variability of the geological formation system in which the groundwater flows.

#### 3. Solutions of the Modified Groundwater Flow Equation

Numerical methods yield approximate solutions to the governing equation through the discretisation of space and time. Within the discretised problem domain, the variable internal properties, boundaries, and stresses of the system are approximated. Deterministic, distributed-parameter, and numerical models can relax the rigid idealised conditions of analytical models or lumped-parameter models, and they can therefore be more realistic and flexible for simulating fields conditions. Our next concern in this paper is to provide solution of the above equation. To achieve this we will make use of two techniques including: the green function and the variational iteration method. We will start with the variational iteration method.

##### 3.1. Variational Iteration Method

The values of the variational iteration method and its applications for a range of categories of differentials equations can be viewed in [13–15].

Following the work recently done by Theis in 1935 [16], in which they proposed an analytical solution to the standard version of the groundwater flow equation, this solution can be approximated as where is the constant discharge rate, is the transmissivity of the aquifer and the piezometric head. It follows that the right side of (13a) and (13b) is known. On the basis of the above equation and knowing the function , one can derive a solution of (13a) and (13b) where the unknown is the function. To make things simple, we put . And (13a) and (13b) become To solve (20) by means of variational iteration method, we put (20) in the form The correction functional for (21) can be approximately expressed as follows for this matter as where is a general Lagrange multiplier [17], which can be recognized optimally by means of variation assumption [17–19], here , , and are considered as constrained variations. Making the above functional stationary Capitulates the next Lagrange multipliers, giving up to the following Lagrange multipliers for the case where and for . For these matter if , we obtained the following iteration formula: Hence we commerce with Means that before the water is pumped out from the borehole, the water level in the aquifer is the same and is considered here to be zero level.

It is worth noting that if the zeroth component is defined, then the remaining components ≥ 1 can be completely determined such that each term is determined by using the previous terms, and the series solutions are thus entirely determined. Finally, the solution is approximated by the truncated series We follow next with the second component To calculate we first need to define explicitly the function . The following function we define here does not actually have a physical meaning, but we use it as example.

To make thing simple, we suppose that and and the function becomes In this matter two components of the decomposition series were obtained of which was evaluated to have the following expansion:

##### 3.2. Green Function Methods

To solve (13a) and (13b), we go on to construct a suitable green’s function for this case in point. Let ) be the green’s function to be constructed, where and . is chosen so as to satisfy homogeneous boundary conditions corresponding to the boundary conditions. It is important to notice that the homogeneous solution of (13a) and (13b) is similar to the diffusion equation if one replaces by ; therefore, the green function involved here is the green’s function for the diffusion equation. Since the aquifer is said to be infinite, the green function for flow equation for infinite aquifer is given by [20] Here the function is to be determined by using the boundary condition. The above equation satisfies an important integral property which is valid for .

Consider This equation is an expression of groundwater flow. At a time and at a position, the piezometer is introduced in the borehole that taps the aquifer. The water that is pumped out from the aquifer through the borehole is migrating through the porous media, but in such a way that the total amount of water in the aquifer is reduced as time goes on if there is no recharge. Since (14) still holds, we can observe that In addition, the green’s function used for this purpose is a solution to the following equation: The general solution of (13a) and (13b) can then be given as function of the green function as Here Due to lake of experimental data for this situation, no graphical representation will be presented in this paper. One needs to model the function of uncertainties introduced in this paper and use it for computational simulation and the analytical solution of the modified groundwater flow equation can then be compared with experimental data. Or one can from the standard solution measure the aquifer’s parameters and use it to determine the values of the function and this is not done in this paper.

#### 4. Conclusion

In this paper we modify the standard version of groundwater flow by replacing the standard derivative with Riemann-Liouville fractional derivatives approximations. The modified equations (13a) and (13b) differ from the standard form of groundwater flow equation in three properties. There is a new operator that takes into account the variation in piezometric head and uncertainties function; second, the “force” appears due to the coordinate dependence of uncertainty function; and finally, there is a derivative-free term that depends only on the uncertainties time function. The modified equation takes into account both the flow through the porous media and the effect of variability of the aquifer or the concept of heterogeneity of the aquifers [21]. The modified equation was solved via the green function technique and variational iteration method.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Authors’ Contribution

Abdon Atangana wrote the first draft and both authors revised and submitted the final version.

#### Acknowledgment

This investigation was sponsored by the Lean Claude postdoctoral Claude Leon Foundation Postdoctoral Fellowships 2014.

#### References

- Y. L. Klimontovich,
*Statistical Physics of Opezi Sys Levis*, Kluwer Academic Publishers, Dodrecht, The Netherlands, 1995. - B. B. Mandelbrot,
*The Fractal Geometry of Nature*, W. H. Freeman, New York, NY, USA, 1982. View at MathSciNet - R. Metzler and J. Klafter, “The random walk's guide to anomalous diffusion: a fractional dynamics approach,”
*Physics Report*, vol. 339, no. 1, pp. 1–77, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - W. G. Gloecke and T. F. Nonnemacher, “Fox function representation of non-debye relaxation processes,”
*Journal of Statistical Physics*, vol. 71, no. 3-4, pp. 741–757, 1993. View at Publisher · View at Google Scholar - V. V. Yanovsky, A. V. Chechkin, D. Schertzer, and A. V. Tur, “Levy anomalous diffusion and fractional Fokker-Planck equation,”
*Physica A: Statistical Mechanics and Its Applications*, vol. 282, no. 1, pp. 13–34, 2000. View at Publisher · View at Google Scholar · View at Scopus - A. V. Chechkin, R. Gorenflo, and I. M. Sokolov, “Retarding subdiffusion and accelerating superdiffusion governed by distributed-order fractional diffusion equations,”
*Physical Review E—Statistical, Nonlinear, and Soft Matter Physics*, vol. 66, no. 4, Article ID 046129, 7 pages, 2002. View at Publisher · View at Google Scholar · View at Scopus - Y. L. Kobelev and E. P. Romanov, “The effect of surface fractal characteristics of solid electrolytes on temperature dependence for constant-phase-angle elements,”
*Doklady Physics*, vol. 45, no. 9, pp. 439–442, 2000. View at Publisher · View at Google Scholar - L. V. Kobelev, “On high energy physics and field theory,” in
*Proceedings of the 24th International Workshop on High Energy Physics and Field Theory*, p. 126, Protvino, Russia, June 2001. - A. Atangana and A. Kılıçman, “A possible generalization of acoustic wave equation using the concept of perturbed derivative order,”
*Mathematical Problems in Engineering*, vol. 2013, Article ID 696597, 6 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - J. A. Barker, “A generalized radial flow model for hydraulic tests in fractured rock,”
*Water Resources Research*, vol. 24, no. 10, pp. 1796–1804, 1988. View at Google Scholar · View at Scopus - A. Atangana and P. D. Vermeulen, “Analytical solutions of a space-time fractional derivative of groundwater flow equation,”
*Abstract and Applied Analysis*, vol. 2014, Article ID 381753, 11 pages, 2014. View at Publisher · View at Google Scholar · View at MathSciNet - N. Goldscheider and D. Drew, Eds.,
*Methods in Karst Hydrogeology*, Taylor & Francis, London, UK, 2007. - A. Cloot and J. F. Botha, “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 Publisher · View at Google Scholar · View at Scopus - J.-H. He, “Variational theory for linear Magneto-Electro-Elasticity,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 2, no. 4, pp. 309–316, 2001. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - J.-H. He, “Variational principles for some nonlinear partial differential equations with variable coefficients,”
*Chaos, Solitons and Fractals*, vol. 19, no. 4, pp. 847–851, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet · View at Scopus - 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 of the American Geophysical Union*, vol. 16, pp. 519–524, 1935. View at Google Scholar - Z. M. Odibat and S. Momani, “Application of variational iteration method to nonlinear differential equations of fractional order,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 7, no. 1, pp. 27–34, 2006. View at Google Scholar · View at Scopus - H. Jafari and H. Tajadodi, “He's variational iteration method for solving fractional Riccati differential equation,”
*International Journal of Differential Equations*, vol. 2010, Article ID 764738, 8 pages, 2010. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W.-H. Su, D. Baleanu, X.-J. Yang, and H. Jafari, “Damped wave equation and dissipative wave equation in fractal strings within the local fractional variational iteration method,”
*Fixed Point Theory and Applications*, vol. 2013, article 89, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - P. M. Morse and H. Feshbash,
*Methods of Theoretical Physics*, McGraw-Hill, New York, NY, USA, 1953. - A. Atangana, “A generalized advection dispersion equation,”
*Journal of Earth System Science*, vol. 123, no. 1, pp. 101–108, 2014. View at Google Scholar