 Open Access
 Total Downloads : 16
 Authors : Praveen Kumar M, Dattatreya Hegde, Rekha J, S. R. Sudheendra
 Paper ID : IJERTCONV3IS19056
 Volume & Issue : ICESMART – 2015 (Volume 3 – Issue 19)
 Published (First Online): 24042018
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Analytical Solution of the OneDimensional Solute AdvectionDispersion Equation with Spatially Variable Retardation Factor using a Change of Variable and Integral Transform Technique
Praveen Kumar M.1
1Department of Mathematics,
KS School of Engineering & Management, Bangalore, India
Rekha J2
2Department of Mathematics, Cambridge Institute of Technology, Bangalore,
India
Dattatreya Hegde3 3Department of Mathematics, Govt. First Grade College, Sirsi, Karnataka
India.

R. Sudheendra4
4Department of Mathematics,

John Institute of Technology, Bangalore,
India
Abstract: An analytical first order solution to the onedimensional advectiondispersion equation with spatially variable retardation factor is derived using a generalized integral transform method to investigate the transport of sorbing but otherwise nonreacting solutes in hydraulic homogenous but geochemically heterogeneous porous formations. The solution is derived under conditions of steadystate flow and arbitrary initial and inlet boundary conditions. The results obtained by this solution agree well with the results obtained by numerically inverting Laplace transformgenerated solutions previously published in the literature. The solution is developed for a third or flux type inlet boundary condition, which is applicable when considering resident solute concentrations and a semiinfinite porous medium. For mathematical simplicity it is hypothesized that the sorption processes are based on linear equilibrium isotherms and that the local chemical equilibrium assumption is valid. The result from several simulations, compared with predictions based on the classical advectiondispersion equation with constant coefficients, indicate that at early times, spatially variable retardation affects the transport behavior of sorbing solutes. The zeroth moments corresponding to constant and variable retardation are not necessarily equal. The center of mass appears to move more slowly, and solute spreading is enhanced in the variable retardation case. At late time, when the travel distance is much larger than the correlation scale of the retardation factor, the zeroth moment for the variable retardation case is identical to the case of invariant retardation. The analytical solution presented in this paper provides more flexibility with regard to the inlet conditions.
INTRODUCTION
The impact of spatially variable hydraulic parameters on the transport and spreading of conservative, nonreacting solutes in natural subsurface has been the focus of many recent studies. Gelhar et. al. (1979), Sudheendra (2010, 2011). Aral et.al (1996) and others, have provided methodologies for improving the description and prediction
of nonreacting solute transport in complex structured formations, compared with the prediction based on the classical advectiondispersion equation with constant coefficients. On the other hand, the transport of sorbing solutes in geochemically as well as hydraulically heterogeneous porous media has received little attention.
For the importance case of transport of sorbing solutes in geochemically homogeneous porous media, the effects of sorption are commonly accounted for by a dimensionless retardation factor, which may be defined as the ration of the average interstitial fluid velocity to the propagation velocity of the solute. Excluding the possibilities of mass transport limitations and solute transformation or decay, any observed fluctuations on the retardation factor are attributed solely to the variability of the distribution coefficient, which is an experimentally obtained measure of sorption or solute retention by the solid formation. Sorption processes can be complex and depend on many variables, including temperature, pressure, solution pH, and ionic strength, sorbent surface charge, sorbent sorptive capacity, and the presence of species that complete for sorption sites. Spatial or temporal fluctuations in any of these variables accordingly affect the distribution coefficient and, consequently, the movement of sorbing solutes in subsurface porous media. For example, the distribution coefficient of nonpolar organic solutes (synthetic organic chemicals, major constituents of groundwater toxic pollutants) is correlated with the organic carbon content of the sorbent (Karckhoff (1984) and Sudheendra (2014)). Although such a correlation is not fully reliable for every
solutesorbent system (Curtis and Roberts, 1985), it can explain to some extent the variable retardation observed in field experiments (Roberts et.al 1986).
C 2C
t D1 z2

w C
1 z
1C
(1)
Garabedian (1987) & Sudheendra (2012) employed
Thus, the appropriate boundary conditions for the given
model
spectral methods to analyze reactive solute macro dispersion under the assumption that the loghydraulic conductivity is linearly related to both the porosity and the distribution coefficient. His result indicate that solute spreading is enhanced when there is negative correlation
C(z, 0) 0
0
0
C(0, t) C et C(, t) 0
z 0
t 0 . (2)
t 0
between the loghydraulic conductivity and the distribution coefficient. The present work is focused on the transport of pollutants but otherwise nonreacting solutes under local equilibrium conditions in a onedimensional unsaturated porous medium. Analytical solutions are employed to solve
The problem then is to characterize the concentration as a
function of z and t, where the input condition is assumed at the origin and a second type or flux type homogeneous condition is assumed. C0 is initial concentration. To reduce equation (3) to a more familiar form, we take
the onedimensional advectiondispersion equation with uniform, steady fluid flow conditions and spatially variable
Cz, t z, t Exp w1z
2D
w2t
1
1
4D
t (3)
1
1
retardation factor, for a semiinfinite medium and fluxtype inlet boundary condition.
1 1
Substituting equation (3) into equation (1) gives
The main objective of the study is to provide mathematical
D1
2
2
(4)
model for better understanding of transport of pollutant through unsaturated porous media. A mathematical model
t z
The initial and boundary conditions (2) transform to
is an important tool and can play a crucial role in
understanding the mechanism of groundwater pollution
0, t
w2t
C0 Exp 1 1 t t 0
problems. It is a simplified description of physical reality expressed in mathematical terms. Mathematical models that attempt to simulate atmospheric processes involved in groundwater pollution are based, in general, on the
4D
1
1
z, 0 0
, t 0
(5)
z 0
t 0
equation of mass conservation for individual pollutant species. Such models relate in one equation the effects of all the physical aspects and dynamic processes that influence the mass balance on groundwater which include transport, diffusion, removal of pollutants and loss or
Equation (4) may be solved for a time dependent influx of the fluid at z = 0. The solution of equation (4) may be obtained readily by use of Dhamels theorem (Carslaw and Jaeger, 1947).
transformation through chemical reactions.
If C F x, y, z, t
is the solution of the diffusion
MATHEMATICAL MODEL
The AdvectionDispersion equation along with initial condition and boundary conditions can be written as
equation for semiinfinite media in which the initial concentration is zero and its surface is maintained at concentration unity, then the solution of the problem in
which the surface is maintained at temperature t is
t
C (1 n) S
2C
D 2
w C C
C
t
F x, y, z, t d
t n t z z 0
The equilibrium isotherm between solution and adsorbed This theorem is used principally for heat conduction
phase is given by coefficient.
S K
t
C
d t
, Kd is the distribution
problems, but the above has been specialized to fit this specific case of interest. Consider now the problem in which initial concentration is zero and the boundary is maintained at concentration unity. The boundary
C (1 n) K
2
C D C C
C D C C
w
conditions are
t n
d t z2 z C
2
0, t 0
z, 0 1
t 0
z 0
1 (1 n) K
C D C w C C
n
d t
z2 z
, t 0
t 0
C 2C C
The Laplace transform of equation (4) is
R t
D z2

w z

C
2
Let us take D1 = D/R, w1 = w/R, 1=/R. Initially,
L t D1 z2
saturated flow of fluid of concentration, C = 0, takes place in the porous media.
Hence, it is reduced to an ordinary differential equation
2 p 2C w2t
2
(6)
z, t 0 Exp 1 t
1
1
z
The solution of the equation is
D1
Aeqz Beqz
2
4D1
2
2
2
2
2
Exp 2
2 d Exp
d
(9)
where, q p .
0
0
D1
The boundary condition as z requires that B = 0
where,
2
z
D1 t
and
2
2
and boundary condition at z = 0 requires that
A 1
thus
w2
z
p
1 1
.
the particular solution of the Laplace transformed equation is
4D1
D1
1 eq z
p
Evaluation of the integral solution
The integration of the first term of equation (9) gives
2
2
The inversion of the above function is given in any table of
2
Laplace transforms. The result is
. Exp 2 d
e2
(10)
z 2 2
0 2
1 erf (
2
)
D1 t
e d
z
2 Dt
For convenience the second integral may be expressed on terms of error function (Horenstein, 1945), because this function is well tabulated.
Using Duhamels theorem, the solution of the problem
with initial concentration zero and the time dependent surface condition at z = 0 is
Noting that
2
2
2
t
2
2
2
2 .
( ) 2 e 2 d d
0
2
t
z
2 Dt
The second integral of equation (9) may be written as
Since e is a continuous function, it is possible to
2
1
2
differentiate under the integral, which gives
I Exp 2 d
e2 Exp
d
2
2
0
0
2 2 z z
e d Exp
2
2 0
t 2 D t 32 4D t
2
z 1 1
2
(11)
2 D1 t
e Exp
d
The solution to the problem is
0
z t
2 D
z2 d
Exp 4D t t 32
Since the method of reducing integral to a tabulated function is the same for both integrals in the right side of
1 0 1
equation (11), only the first term is considered. Let
(7) a
and the integral may be expressed as
Putting
2
z
D1 t
then the equation (7) can be
2
2
written as
I1 e
Exp
d
2
z2
2
0
2
t 4D 2 e d
(8)
e2
1 Exp a da
z 1
a2
a
2D1 t
w2t
/
2
(12)
Since t C0 Exp 1 1 t the particular
e2
Exp a da
4D1
/
a
solution of the problem be written as
Further, let,
a

a
2C
w2 t
in the

a first term of the above equation, then
z, t 0 Exp 1 t
1
1
a
2
2
4D1
I1 e2
e d e2 Exp


a
da .
2
1 2
2 2
2
a
(13)
e
2
2
2 e
e
d e
e
d (15)
However, by definition,
Similar evaluation of the second integral of equation (11)
gives
2
2
e2
e 2
d
e2
2
erfc
I2 e
Exp a a
da
/
2
also,
e2
Exp


a
da
2
/
a
e2
e
d
e2 1 erf
.
Again substituting
a
a
into the first term, the
2
result is
Writing equation (15) in terms of error functions, we get
2
C w2t
I 2 e2
e 2 d e2
Exp

a
da .
z, t 0 Exp 1 1 t .
a
2 4D1
/
e2 erfc e 2 erfc
(16)
Noting that
2
Exp a
a
2 da
Thus, Substitution into equation (3) the solution is
/
2
C 1 Exp w1z
t
Exp
a
2 da
C 2 2D
/
a
0 1
e2 erfc e2 erfc
Substitution into equation (11) gives
1
I e
2
2
e2 d e2
e
2
d . (14)
Resubstituting for and gives
Thus, equation (9) may be expressed as
C 1 Exp w1z t
z w2 4D
C 2 2D
erfc 1 1 1 t
0 1
w2 4D
2 D1t
2
Exp 1 1 1 z
w1
4D1 1
2D1
Exp z
2D
1
z w2 4D
erfc 1 1 1 t (17)
2
D1t
where boundaries are symmetrical the solution of the problem is given by the first term the equation (17). The second term is equation (17) is thus due to the asymmetric boundary imposed in the more general problem. However, it should be noted also that if a point a great distance away from the source is considered, then it is possible to
approximate the boundary condition by C , t C0 ,
which leads to a symmetrical solution.

RESULTS & DISCUSSIONS:
This study presents analytical solutions for one dimensional advectiondispersion equations in unsaturated porous medium in finite domain. The transform method coupled with the generalized integral transform technique is used to obtain the analytical solutions. Solutions are obtained for both first and thirdtype inlet boundary conditions. The developed analytical solutions for finite domain are compared with solutions for the semiinfinite domain to clarify how the exit boundary influences the onedimensional transport in a porous medium system.
Fig. 2: Breakthroughcurve for C/C0 v/s time for z=10m, R=1.0, =0.5 & = 0
Fig. 3: Breakthroughcurve for C/C0 v/s time for z=10m, R=1.0, =0.5 & = 0.25
Fig. 4: Breakthroughcurve for C/C0 v/s time for z=10m, R=1.0, =0.5 & = 0.5
Fig. 5: Breakthroughcurve for C/C0 v/s time for z=10m, R=1.0, =0.5, = 0.75 & 1.0
The main limitations of the analytical methods are that the applicability is for relatively simple problems. The geometry of the problem should be regular. The properties of the soil in the region considered must be homogeneous in the sub region. The analytical method is somewhat more flexible than the standard form of other methods for one dimensional transport model. Figures 1 to 4 represents the concentration profiles verses time in the adsorbing media for depth z = 10m and Retardation factor R=1. It is seen that for a fixed velocity w, dispersion coefficient D and distribution coefficient Kd, C/C0 decreases with depth as porosity n decreases due to the distributive coefficient Kd and if time increases the concentration decreases for different time and decay chain.
Accordingly, the analytical solutions derived for the finite domain will thus be particularly useful for analyzing the onedimensional transport in unsaturated porous medium with a large dispersion coefficient whereas the analytical solution for semiinfinite domain is recommended to be applied for a medium system with a small dispersion coefficient. Moreover, the developed solution is especially useful for validating numerical model simulated solution because realistic problems generally have a finite domain.
From this paper, we conclude that the mathematical solutions have been developed for predicting the possible concentration of a given dissolved substance in steady unidirectional seepage flows through semiinfinite, homogeneous, and isotropic porous media subject to source concentration that vary exponentially with time for spatially variable retardation factor using a change of variable and integral transform technique. The expressions take into account the contaminants as well as mass transfer from the liquid to the solid phase due to adsorption. For simultaneous dispersion and adsorption of a solute, the dispersion system is considered to be adsorbing at a rate proportional to its concentration.

REFERENCES:

Aral, M.M., Liao, B., 1996. Analytical solutions for two dimensional transport equation with timedependant dispersion coefficients. Journal of Hydrologic Engineering, 1,2032.

Barry, D. A., and Sporito, G. , 1989. Analytical solution of a convectiondispersion model with timedependant transport co efficients. Water Resour.Res.,25,24072416.

Batu, V.,1993. A generalized twodimensional analytical solute transport model in bounded media for fluxtype multiple sources. Water Resour. Res., 29, 28812892.

Bear,J.,and A.Verruijt.,1990.Modelling Groundwater flow and pollution.D Radial Publishing Co.,Tokyo.

Ermak,D.L.,1977.An Analytical Model for Air Pollutant transport and deposition from a point source.Atmos.Environ.,11. [6]J.S.Chen,C.W.Liu,and C.M.Liao., 2003. Two – dimensional Laplace – Transformed Power Series Solution for Solute Transport in a Radially Convergent Flow Field. Adv. Water Res.,26,11131124.

Koch,W.,1989.A Solution of twodimensional atmosphere diffusion equation with heightdependent diffusion coefficient including ground level absorption. Atmos. Environ., 23, 17291732.

Sudheendra S.R., 2010 A solution of the differential equation of longitudinal dispersion with variable coefficients in a finite domain, Int. J. of Applied Mathematics & Physics, Vol.2, No. 2, 193204.

Sudheendra S.R., 2011. A solution of the differential equation of dependent dispersion along uniform and nonuniform flow with variable coefficients in a finite domain, Int. J. of Mathematical Analysis, Vol.3, No. 2, 89105.

Sudheendra S.R. 2012. An analytical solution of onedimensional advectiondiffusion equation in a porous media in presence of radioactive decay, Global Journal of Pure and Applied Mathematics, Vol.8, No. 2, 113124.

Sudheendra S.R., Raji J, & Niranjan CM, 2014. Mathematical Solutions of transport of pollutants through unsaturated porous media with adsorption in a finite domain, Int. J. of Combined Research & Development, Vol. 2, No. 2, 3240.

Sudheendra S.R., Praveen Kumar M. & Ramesh T. 2014. Mathematical Analysis of transport of pollutants through unsaturated porous media with adsorption and radioactive decay, Int. J. of Combined Research & Development, Vol. 2, No. 4, 01 08.

Sudheendra S.R., Raji J, & Niranjan CM, 2014. Mathematical modelling of transport of pollutants in unsaturated porous media with radioactive decay and comparison with soil column experiment, Int. Scientific J. on Engineering & Technology, Vol. 17, No. 5.

Tartakowsky,D.,Di Federico,V.,1997. An analytical solution for contaminant transport in nonuniform flow. Transport in porous media, 27, 8597.

Wexler.E.J.,1992.Analytical solution for one, two and three dimensional solute transport in Ground water systems with uniform flow. U.S. Geological Survey, Techniques of water Resources Investigations, Book 3,Chap. B.7.

Yates, S. R., 1990. An analytical solution for onedimensional transport in heterogeneous porous media. Water Resour. Res., 26, 23312338.

Zoppou,C., and Knight, J.H.., 1997. Analytical solution for advection and advectiondiffusion equation with spatially variable coefficients. Journal of Hydraulic Engineering.,123,144148