Abstract and Applied Analysis

Volume 2013 (2013), Article ID 206942, 7 pages

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

## Two-Dimension Hydrodynamic Dispersion Equation with Seepage Velocity and Dispersion Coefficient as Function of Space and Time

^{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 26 April 2013; Accepted 5 June 2013

Academic Editor: L. Jódar

Copyright © 2013 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 contamination through the geological formation cannot move and disperse with the same speed and dispersion coefficient, respectively, due to the variability of the geological formation. This paper is therefore first devoted to the description of the hydrodynamic advection dispersion equation with the seepage velocity and dispersion coefficient as function of space and time. Secondly the equation is solved via two analytical techniques: the homotopy decomposition method and the differential transform method. The numerical simulations of the approximated solutions are presented.

#### 1. Introduction

Most of the Earth’s liquid freshwater is found, not in lakes and rivers, but is stored underground in conglomerations of voids, called aquifers, and present in the geological formations constituting the earth’s mantle, where they act as reservoirs of water for rivers and streams, especially during periods of drought or low rainfall—a phenomenon commonly called baseflow. This resource, conventionally referred to as groundwater, therefore forms an essential, in fact major, component of the fresh water resources, as illustrated by the fact that that nearly two billion people depend directly on groundwater or their drinking water, while 40% of the world’s food is produced on farmlands irrigated with groundwater [1].

One reason why groundwater so often constitutes the main source of drinking water in many cities and towns around the world is because it is frequently present in sufficient quantities at the point of demand. However, this seemingly advantage may sometimes be its greatest disadvantage, especially in situations where the groundwater occurs at shallow depths and the area overlying the aquifer is populated densely. This problem is particularly relevant in the present technological age to its vast quantities of waste that is often disposed in an uncontrolled manner. Since space, time, and mass are measurable, and numbers form the basis of all mathematics, the question arises if it is not possible to combine the power of human reason with well-planned observations and describe natural phenomena in abstract mathematical terms? This question leads us to the mathematical formulation of the flow of subsurface water in the geological formations. The development of a theory for natural phenomenon can be briefly summarized in four steps [2, 3].

Apply the rules of mathematical analysis to the observables and try to establish a framework, or hypothesis, able to predict the behaviour of the phenomenon, under different conditions from those used in establishing the hypothesis. The main effort is to establish a suitable mathematical relation between different observables, or what [4] calls a mathematical model for the observables. An ideal mathematical model should not only provide links between different observables, but also lead to a better understanding of the phenomenon. The attractiveness of this approach is that the mathematical model could, in principle, be used to investigate the future behaviour of a given phenomenon and under various conditions. For simplicity, it is always assumed that, the seepage velocity retardation factor and the dispersion coefficient appearing are the mathematical equation describing the movement of contamination in groundwater as constants. However, in the real world observation, these coefficients depend on space and time. One therefore needs to include such behaviour in the mathematical formulation. This paper is therefore devoted to the discursion underpinning the effect of variation in time space of the seepage velocity and dispersion coefficient on one hand, and on the other hand a possible analytical solution of this equation.

#### 2. Mathematical Formulation

The mathematical model, discussed below, has been used extensively in the studies of groundwater pollution problems, with satisfactory results. A relatively complete set of one-dimensional analytical solutions for convective-dispersive solute equations has been recently published by [5]; here we will review a case having a practical application. Let us consider a two-dimensional model consisting of and infinitely ling homogenous isotropic porous media with steady-state uniform flow with seepage velocity . We inject a particular chemical from one end of the model for a period of time such that the input concentration varies as an exponential function of time [5]. The value of that chemical concentration at any time and at a distance from the injection boundary, allowing for the decay and adsorption, may be obtained from the solution of the following set of equations [5]: where is the dispersion coefficient and is the retardation factor and subject to the initial condition: which means that the system is initially free of that chemical, and are constants, and is the radioactive decay constant. And boundaries conditions: In order to include explicitly the variation of the seepage velocity and dispersion coefficient into the mathematical formulation, we replace the constants coefficients , and of (1) by the functions , and to obtain To satisfy the physical behaviour of the aquifer, there exist positives constants , and such that

#### 3. Analytical Solutions

Environmental phenomena, such as groundwater pollution, are highly complex phenomena, which do not lend themselves readily to analysis of analytical models. The discussion presented in this section will therefore be devoted to the derivation of analytical solution to the advection dispersion equation (1). As V. M. Alexandrov wrote in the foreword of a popular science book “Asymtotology: ideas, methods, and applications [sic]”, asymptotic methods belong to the, perhaps, the most romantic area of modern mathematics [6]. Though computer science is growing very fast and numerical simulation is applied everywhere, nonnumerical issues will still play a large role [7–11].

There exist some alternative analytical asymptotic approaches, such as the nonperturbative method [12], modified Lindstedt-Poincare method [13], variational iteration method [14], Adomian decomposition method [14], homotopy perturbation method [15], and book-keeping artificial parameter perturbation method [15]. Recently to solve the groundwater flow equation, Atangana and Botha proposed a modified decomposition method called homotopy decomposition [16]. In this section we make use of two different methods including the HDM and DTM to derive approximate solution of (1). We will start with HDM.

##### 3.1. Useful Tools for the Homotopy Decomposition Method [17–19]

To accommodate readers that are not acquainted with this modified form of decomposition method, we present in this section the basic methodology of the method. To illustrate the fundamental design of this technique, we consider a universal nonlinear nonhomogeneous partial differential equation with initial conditions of the following form: subject to the initial condition: where is a known function, is the general nonlinear differential operator, and represents a linear differential operator. The first step of the method here is to apply the inverse operator of on both sides of (1) to obtain The multi-integral in (2) can be transformed toso that (8) can be reformulated as Using the Homotopy scheme the solution of the above integral equation is given in series form as and the nonlinear term can be represented as where is an embedding parameter. is the polynomials that can be engendered by The homotopy decomposition method is achieved via the polished combination of decomposition method with Abel’s integral and is given by With matching up to the expressions of similar powers of , we can get hold of solutions of assorted orders. The early estimate of the approximation is . It is significant to draw attention to the fact that the early estimate is the Taylor series of order of the exact solution.

##### 3.2. Differential Transform Method (DTM) [17, 20–22]

The necessary definitions of the three-dimensional transform are shown. Given a function which has three components such as , the three-dimensional differential transform function of the function is defined as where is the original function and is the transform function.

The inverse differential transform of is defined as
from (16) and (17) it can be concluded that
Some fundamental mathematical operators performed by three-dimensional differential transform methods are listed below [20]:(1) if , then ,(2) if , then where is a constant,(3) if **,** then ,(4) if **,** then ,(5) if **,** then ,(6) if , then .

#### 4. Applications

In this section, we apply these methods for solving the two-dimension hydrodynamic dispersion equation with seepage velocity and dispersion coefficient as function of space and time.

Since , we first divide (1) by and let to obtain Applying the HDM to the above equation yields Following the DTM method we supposed that the solution of (1) is in the form of with so that

*Example 1. *Let us consider in (1) the following coefficients:
The theoretical parameters used in this simulation are .

The following figure shows the numerical simulations of the approximated solutions obtained via HDM and DTM.

Figures 1 and 2 show the numerical simulation of the approximate solution via HDM and DTM, respectively. Here we fixed and we consider the concentration to be a function of and . It can be observed that both methods give pretty much the same results for the first four components of the series solution.

Figures 3 and 4 show the numerical simulation of the approximate solution via HDM and DTM, respectively. Here we fixed and we consider the concentration to be a function of and .

Figures 5 and 6 are the contour plot obtained via Mathematica. These figures show the possible pathway taken by the pollution through the geological formation called aquifer.

Figure 7 shows the contour plot of the approximate solution as a function of and for a fixed time. In practice this figure shows the behaviour of the contamination through the geological formation for a fixed time . The figure shows that since the seepage velocity and the dispersion coefficient are function of space, the pollution will not have the same intensities everywhere in the aquifer. There will be certain place in the aquifer with more pollution than other, depending on the velocity and the dispersion coefficients associated.

#### 5. Discussion and Conclusion

Natural geological deposits with highly contrasting permeability may form mobile and relatively immobile zones, where the potential mass exchange between mobile and immobile zones results in a wide time distribution for solute trapping. The transport process in groundwater is, by its very nature, always in contact with the matrix of an aquifer. There is thus a possibility that the solutes may interact with the rock matrix and one another. A true mathematical model for groundwater pollution must therefore be able to account for interactions between the dissolved solids and matrix of the aquifer. It will thus be advantageous to look at the nature of the interactions between dissolved solids and a porous medium that may be expected in groundwater pollution. Experimental evidence indicates that when a dissolved solid comes in contact with the matrix of a porous medium, it may (a) pass through the medium with no apparent effect, (b) be absorbed by the porous matrix, and (c) react with the porous matrix and other substances dissolved in the fluid. The dissolved solids encountered in porous flow are, for this reason, often classified as conservative, nonconservative, and reactive tracers [22]. This behaviour implies that the quantity of dissolved solids in a porous medium depends not only on the flow pattern, but also on the nature of the porous matrix and the solution. These situations (a), (b), and (c) can be characterized efficiently by the time-nonlocal model, including the hydrodynamic dispersion equation with seepage velocity and dispersion coefficients as function of space and time as described by Figure 5. If the high-permeable material tends to form preferential flow paths, such as the interconnected paleochannels observed in alluvial depositional systems, then the solute transport may show a heavy leading edge, which can be described by the hydrodynamic dispersion equation with seepage velocity and dispersion coefficients as function of space and time.

#### References

- B. L. Morris, A. R. L. Lawrence, P. J. C. Chilton, B. Adams, R. C. Calow, and B. A. Klinck,
*Groundwater and Its Susceptibility to Degradation: A Global Assessment of the Problem and Options for Management*, Early Warning and Assessment Report Series, RS.03-3, United Nations, 2003. - J. F. Botha,
*Principle of Groundwater Motion*, Institute for Groundwater Study, University of the Free State, Bloemfontein, South Africa, 1996. - J. Davies,
*The Scientific Approach*, Academic Press, London, UK, 2nd edition, 1973. - J. F. Botha, J. P. Verwey, J. Buys, G. Tredoux, J. W. Moodie, and M. Hodgkiss, “Modelling groundwater contamination in the atlantics aquifer,”
*WRC Report*175/1/90, Water Research Commission, Pretoria, South Africa, 1990. View at Google Scholar - M. T. van Genuchten and W. J. Alves,
*Analytical Solutions of One Dimensional Convective-Dispersive Solute Transport Equations*, Technical Bulletin, no. 1661, United State Department of Agriculture, 1982. - I. V. Andrianov and L. I. Manevitch,
*Asymptotology: Ideas, Methods, and Applications*, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2003. - I. Andrianov and J. Awrejcewicz, “Construction of periodic solutions to partial differential equations with non-linear boundary conditions,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 1, no. 4, pp. 327–332, 2000. View at Publisher · View at Google Scholar · View at MathSciNet - C. M. Bender, K. A. Milton, S. S. Pinsky, and L. M. Simmons, Jr., “A new perturbative approach to nonlinear problems,”
*Journal of Mathematical Physics*, vol. 30, no. 7, pp. 1447–1455, 1989. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Delamotte, “Nonperturbative (but approximate) method for solving differential equations and finding limit cycles,”
*Physical Review Letters*, vol. 70, no. 22, pp. 3361–3364, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Atangana, “New class of boundary value problems,”
*Information Sciences Letters*, vol. 1, no. 2, pp. 67–76, 2012. View at Publisher · View at Google Scholar - J. Awrejcewicz and V. A. Krysko,
*Introduction to Asymptotic Methods*, CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - 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. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J.-H. He, “Homotopy perturbation technique,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 178, no. 3-4, pp. 257–262, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Kurulay, A. Secer, and M. A. Akinlar, “A new approximate analytical solution of Kuramoto-Sivashinsky equation using homotopy analysis method,”
*Applied Mathematics & Information Sciences*, vol. 7, no. 1, pp. 267–271, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - J.-H. He, “Bookkeeping parameter in perturbation methods,”
*International Journal of Nonlinear Sciences and Numerical Simulation*, vol. 2, no. 3, pp. 257–264, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Atangana and J. F. Botha, “Analytical solution of the groundwater flow equation obtained via homotopy decomposition method,”
*Journal of Earth Science & Climatic Change*, vol. 3, article 115, 2012. View at Publisher · View at Google Scholar - A. Atangana and A. Secer, “The time-fractional coupled-Korteweg-de-Vries equations,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 947986, 8 pages, 2013. View at Publisher · View at Google Scholar - A. Atangana, O. Aden Ahmed, and N. Bıldık, “A generalized version of a low velocity impact between a rigid sphere and a transversely isotropic strain-hardening plate supported by a rigid substrate using the concept of noninteger derivatives,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 671321, 9 pages, 2013. View at Publisher · 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 - A. Secer, “Solving time-fractional reaction-diffusion equation by reduced differential transform method,”
*International Journal of Science and Emerging Technologies*, vol. 3, no. 1, 2012. View at Google Scholar - A. Secer, “Approximate analytic solution of fractional heat-like and wave-like equations with variable coefficients using the differential transforms method,”
*Advances in Difference Equations*, vol. 2012, article 198, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - A. Atangana,
*A generic assessment of waste disposal at Douala city: principle, practices and uncertainties [Ph.D. thesis]*, Library of the University of the Free State, Bloemfontein, South Africa, 2013.