The classical Darcy law is generalized by regarding the water flow as a function of a noninteger order derivative of the piezometric head. This generalized law and the law of conservation of mass are then used to derive a new equation for groundwater flow. Two methods including Frobenius and Adomian decomposition method are used to obtain an asymptotic analytical solution to the generalized groundwater flow equation. The solution obtained via Frobenius method is valid in the vicinity of the borehole. 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 boreholes 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. Numerical solutions obtained via Adomian method are compared with the Barker generalized radial flow model for which a fractal dimension for the flow is assumed. Proposition for uncertainties in groundwater studies was given. 1. Introduction The real problem encounter in groundwater studies up to now is the real shape of the geological formation in which water flows in the aquifer under investigation. However, there are many fractured rock aquifers where the flow of groundwater does not fit conventional geometries [1], for example, in South Africa, the Karoo aquifers, characterized by the presence of a very few bedding parallel fractures that serve as the main conduits of water in the aquifers [2]. With a challenge to fit the solution of the groundwater flow equation with experimental data from field observation in particular, the observed drawdown see [3], for all time yields a fit that undervalues the observed drawdown at early times and overvalues it at later times. The variation of observations from theoretically predictable values is usually an indication of uncertainties in the predictable. To investigate the first possibility Botha et al. [2] developed a three-dimensional model for the Karoo aquifer on the campus of the University of the Free State. This model is based on the conventional, saturated groundwater flow equation for density-independent flow: where is the specific storativity, the hydraulic conductivity tensor of the aquifer, the piezometric head, ) the strength of any sources or sinks, with and the usual spatial and time coordinates; the gradient operator, and the time derivative. This model showed that the dominant flow field in these aquifers is vertical and linear and not horizontal and
References
[1]
J. H. Black, J. A. Barber, and D. J. Noy, “Crosshole investigations: the method, theory and analysis of crosshole sinusoidal pressure tests in fissured rocks,” Stripa Projects Internal Reports 86-03, SKB, Stockholm, Sweden.
[2]
J. F. Botha, I. Verwey van Voort, J. J. P. Viviers, W. P. Collinston, and J. C. Loock, “Karoo aquifers. Their geol-ogy, geometry and physical behaviour,” WRC Report 487/1/98, Water Research Commission, Pretoria, 1998.
[3]
G. J. van Tonder, J. F. Botha, W.-H. Chiang, H. Kunstmann, and Y. Xu, “Estimation of the sustainable yields of boreholes in fractured rock formations,” Journal of Hydrology, vol. 241, no. 1-2, pp. 70–90, 2001.
[4]
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.
[5]
J. Bear, Dynamics of Fluids in Porous Media, American Else-vier Environmental Science Series, Elsevier, New York, NY, USA, 1972.
[6]
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.
[7]
R. Courant and F. John, Introduction to Calculus and Analysis, vol. 2, John Wiley & Sons, New York, NY, USA, 1974.
[8]
Y. Cherruault and G. Adomian, “Decomposition methods: a new proof of convergence,” Mathematical and Computer Modelling, vol. 18, no. 12, pp. 103–106, 1993.
[9]
I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, Calif, USA, 1999.
[10]
I. Podlubny, “Geometric and physical interpretation of fractional integration and fractional differentiation,” Fractional Calculus & Applied Analysis, vol. 5, no. 4, pp. 367–386, 2002.
[11]
K. Adolfsson, “Nonlinear fractional order viscoelasticity at large strains,” Nonlinear Dynamics, vol. 38, no. 1–4, pp. 233–246, 2004.
[12]
O. P. Agrawal, “Application of fractional derivatives in thermal analysis of disk brakes,” Nonlinear Dynamics, vol. 38, pp. 191–206, 2004.
[13]
G. Afken, Mathematical Methods for Physicists, Academic Press, London, UK, 1985.
[14]
K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley & Sons, New York, NY, USA, 1993.
[15]
K. B. Oldham and J. Spanier, The Fractional Calculus, Academic Press, New York, NY, USA, 1974.
[16]
G. Adomian, “A review of the decomposition method in applied mathematics,” Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501–544, 1988.
[17]
A. Atangana, “New class of boundary value problems,” Information Science Letters, vol. 1, no. 2, pp. 67–76, 2012.
[18]
G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, Kluwer Academic, Dodrecht, The Netherlands, 1994.
[19]
Y. Cherruault, “Convergence of Adomian's method,” Kybernetes, vol. 18, no. 2, pp. 31–38, 1989.
