Abstract and Applied Analysis

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

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

Research Article

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

Abdon Atangana¹ and S. C. Oukouomi Noutchie²

¹Institute for Groundwater Studies, Faculty of Natural and Agricultural Sciences, University of the Free State, Bloemfontein 9300, South Africa

²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 to

so 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.

Figure 1: Simulation of the approximate solution obtained via HDM for a fixed .

Figure 2: Simulation of the approximate solution obtained via DTM for a fixed .

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 .

Figure 3: Simulation of the approximate solution obtained via HDM for a fixed .

Figure 4: Simulation of the approximate solution obtained via DTM for a fixed .

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 5: Contour plot of the approximated solution in x-direction for fixed .

Figure 6: Contour plot of the approximate solution in y-direction for fixed .

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.

Figure 7: Contour plot of the approximated solution for a fixed time .

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.

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.

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.

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.

A. Atangana, “New class of boundary value problems,” Information Sciences Letters, vol. 1, no. 2, pp. 67–76, 2012.

J. Awrejcewicz and V. A. Krysko, Introduction to Asymptotic Methods, CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC, Boca Raton, Fla, USA, 2006.

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.

J.-H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 3-4, pp. 257–262, 1999.

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.

J.-H. He, “Bookkeeping parameter in perturbation methods,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 2, no. 3, pp. 257–264, 2001.

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.

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.

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.

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.

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.

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.

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.