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

and Applied Analysis 3 a universal nonlinear nonhomogeneous partial differential equation with initial conditions of the following form:


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

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 onedimensional 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 V. We inject a particular chemical from one end of the model for a period of time  0 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  0 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 To satisfy the physical behaviour of the aquifer, there exist positives constants  1 ,  2 , and  3 such that (5)

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][8][9][10][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 equation ( 1).We will start with HDM.[17][18][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:

Useful Tools for the Homotopy Decomposition Method
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 +  ( (, )) + ℎ (, ) .(10) 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  ∈ (0, 1] is an embedding parameter.P  () is the polynomials that can be engendered by ,  = 0, 1, 2 . . . .
The homotopy decomposition method is achieved via the polished combination of decomposition method with Abel's integral and is given by 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.[17,[20][21][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

Differential Transform Method (DTM)
where (, , ) is the original function and (, ℎ, ) is the transform function.
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  = 2 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  = 2 and we consider the concentration to be a function of  and .  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  = 10.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.

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.

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

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

Figures 5 and 6
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.Figure7shows 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  = 10.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 4 :
Figure 4: Simulation of the approximate solution obtained via DTM for a fixed  = 2.
of c(x, y = 2, t) in the aquifer

Figure 5 :
Figure 5: Contour plot of the approximated solution in x-direction for fixed .
c(x = 2, y, t) in the aquifer y distance

Figure 6 :
Figure 6: Contour plot of the approximate solution in y-direction for fixed .
of c(x, y, t = 10) in the aquifer y distance

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