# An Analytical Solution of Transport of Pollutants in Unsaturated Porous Media with and Without Adsorption

Text Only Version

#### An Analytical Solution of Transport of Pollutants in Unsaturated Porous Media with and Without Adsorption

Niranjan C. M.1

1Dept. of Mathematics, Acharya Institute of Technology,

Bangalore, Karnataka India.

Anitha R2.

2Dept. of Mathematics,

Nitte Meenakshi Institute of Technology, Bangalore, Karnataka

India.

3Dept. of Mathematics Govt. First Grade College Hethur, Sakleshpur(tq) India

1. R. Sudheendra4

4Dept. of Mathematics,

2. John Institute of Technology, Bangalore, Karnataka

India.

Abstract: Most of the investigators use the coordinate transformation (z – ut) in order to solve the equation for dispersion of a moving fluid in porous media. Further, the boundary conditions C = 0 at z = and C = C0 at z = for t

0 are used, which results in a symmetrical concentration distribution. In this paper, the effect of adsorption has been studied for one-dimensional transport of pollutants through the unsaturated porous media. In this study, the advection- dispersion equation has been solved analytically to evaluate the transport of pollutants which takes into account of dissipation coefficient and porosity by considering input concentrations of pollutants that vary with time and depth. The solution has been obtained using Laplace transform, moving coordinates and Duhamels theorem is used to get the solution in terms of complementary error function.

Key words: Advection, dispersion, adsorption, Integral transforms, Ficks law, Moving coordinates, Duhamels theorem

1. INTRODUCTION:

In recent years considerable interest and attention have been directed to dispersion phenomena in flow through porous media. Scheidegger (1954), deJong (1958), and Day (1956) have presented statistical means to establish the concentration distribution and the dispersion coefficient. Advectiondiffusion equation describes the solute transport due to combined effect of diffusion and convection in a medium. It is a partial differential equation of parabolic type, derived on the principle of conservation of mass using Ficks law. Due to the growing surface and subsurface hydro environment degradation and the air pollution, the advectiondiffusion equation has drawn significant attention of hydrologists, civil engineers and mathematical modelers. Its analytical/numerical solutions along with an initial condition and two boundary conditions help to understand the contaminant or pollutant concentration distribution behavior through an open

medium like air, rivers, lakes and porous medium like aquifer, on the basis of which remedial processes to reduce or eliminate the damages may be enforced. It has wide applications in other disciplines too, like soil physics, petroleum engineering, chemical engineering and biosciences.

In the initial works while obtaining the analytical solutions of dispersion problems in ideal conditions, the basic approach was to reduce the advectiondiffusion equation into a diffusion equation by eliminating the convective term(s). It was done either by introducing moving co-ordinates (Ogata and Banks 1961; Harleman and Rumer 1963; Bear 1972; Guvanasen and Volker 1983; Aral and Liao 1996; Marshal et al 1996) or by introducing another dependent variable (Banks and Ali 1964 Ogata 1970; Lai and Jurinak 1971; Marino 1974 and Al-Niami and Rushton 1977). Then Laplace transformation technique has been used to get desired solutions.

Some of the one-dimensional solutions have been given (Tracy 1995) by transforming the non-linear advection diffusion equation into a linear one for specific forms of the moisture content vs. pressure head and relative hydraulic conductivity vs. pressure head curves which allow both two-dimensional and three-dimensional solutions to be derived. A method has been given to solve the transport equations for a kinetically adsorbing solute in a porous medium with spatially varying velocity field and dispersion coefficients (Van Kooten 1996, Sudheendra et.al. 2014).

Later it has been shown that some large subsurface formations exhibit variable dispersivity properties, either as a function of time or as a function of distance (Matheron and deMarsily 1980; Sposito et al

1986; Gelhar et al 1992). Analytical solutions were developed for describing the transport of dissolved substances in heterogeneous semi infinite porous media

C 2C

t D z2

• w C

z

(1 n) n

Kd C

(1)

with a distance dependent dispersion of exponential nature along the uniform flow (Yates 1990, 1992). The temporal moment solution for one dimensional advective-dispersive solute transport with linear equilibrium sorption and first order degradation for time pulse sources has been applied to analyze soil column experimental data (Pang et al 2003, Sudheendra et.al. 2014). An analytical approach was developed for non-equilibrium transport of reactive solutes in the unsaturated zone during an infiltrationredistribution cycle (Severino and Indelman 2004).

The solute is transported by advection and obeys

where C is the constituent concentration in the soil solution, t is the time in minutes, D is the hydrodynamic dispersion coefficient, z is the depth, u is the average pore-

1 n

water velocity and n Kd C is the adsorption term.

Initially saturated flow of fluid of concentration, C = 0, takes place in the medium. At t = 0, the concentration of the plane source is instantaneously changed to C = C0. Then the initial and boundary conditions (Fig. 1) for a semi-infinite column and for a step input are

linear kinetics. Analytical solutions were presented for

solute transport in rivers including the effects of transient

C z, 0 0 ;

z 0

storage and first order decay (Smedt 2006, Sudheendra

C 0, t C0 ;

t 0

(2)

2011, 2012). Pore flow velocity was assumed to be a non- divergence, free, unsteady and non-stationary random

C , t 0 ;

t 0

function of space and time for ground water contaminant transport in a heterogeneous media (Sirin 2006). A two- dimensional semi-analytical solution was presented to analyze streamaquifer interactions in a coastal aquifer where groundwater level responds to tidal effects (Kim et al 2007).

The problem then is to characterize the concentration as a function of x and t.

To reduce equation (1) to a more familiar form, let

wz w2t (1 n)

A more direct method is presented here for

Cz, t z, t exp 2D 4D n

Kdt (3)

solving the differential equation governing the process of dispersion. It is assumed that the porous medium is homogeneous and isotropic and that no mass transfer

Substitution of equation (3) reduces equation (1) to Ficks law of diffusion equation

occurs between the solid and liquid phases. It is assumed also that the solute transport, across any fixed plane, due to

D

2

(4)

microscopic velocity variations in the flow tubes, may be quantitatively expressed as the product of a dispersion coefficient and the concentration gradient. The flow in the medium is assumed to be unidirectional and the average velocity is taken to be constant throughout the length of the flow field. In this paper, the solutions are obtained for two solute dispersion problems in a longitudinal finite length, respectively. In the first problem time dependent solute dispersion of increasing or decreasing nature along a uniform flow trough a homogeneous domain is studied. The input condition is of uniform and varying nature, respectively.

2. MATHEMATICAL FORMULATION AND MODEL We consider one-dimensional unsteady flow

through the semi-infinite unsaturated porous media in the

x-z plane in the presence of a toxic material. The uniform flow is in the z-direction. The medium is assumed to be isotropic and homogeneous so that all physical quantities are assumed to be constant. Initially the concentration of the contaminant in the media is assumed to be zero and a constant source of concentration of strength C0 exists at the surface. The velocity of the groundwater is assumed to be constant. With these assumptions the basic equation governing the flow is

t z2

Figure 1 : Physical Layout of the Model

The above initial and boundary conditions (2) transform to

w2t

(1 n)

The boundary condition as z requires that B = 0 and

0, t C0

exp

4D

n Kd t ;

t 0

boundary condition at z = 0 requires that

A 1 , thus

z, 0 0 ;

z 0

p

p

the particular solution of the Laplace transform equation is

, t 0 ;

t 0

1

p

e qz

(5)

It is thus required that equation (4) be solved for a time dependent influx of fluid at z = 0. The solution of equation

The inversion of the above function is given in a table of Laplace transforms (Carslaw and Jaeger, 1947). The result is

(4) can be obtained by using Duhamels theorem [Carslaw

x

2 2

and Jeager, 1947].

1 erf

2

e

Dt z

d.

(7)

If C Fx, y, z, t

is the solution of the diffusion

2 Dt

equation for semi-infinite 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

Utilizing Duhamels theorem, the solution of the problem with initial concentration zero and the time dependent surface condition at z = 0 is

t

t

t

C 0 ( ) t F (x, y, z, t ) d .

2

2

This theorem is used principally for heat conduction problems, but above has been specialized to fit this specific

0

t

e

z

d d

case of interest.

2 Dt

since

e 2

is a continuous function, it is possible

Consider now the problem in which initial concentration is

zero and the boundary is maintained at concentration unity.

differentiate under the integral, which gives

z 2

The boundary conditions are

2

t

e 2 d

z 4 Dt

3

2

e

e

0, t 1 ; t 0

z 2

D t

2 Dt

z, 0 0 ; z 0 .

The solution of the problems is

• z 2

• z 2

, t 0 ;

t 0

z t

4 Dt d

This problem can be solved by the application of the Laplace transform. The concentration which is function of t and whatever space coordinates, say z, t, occur in the problem. We write

Letting

t e

D

D

0

z

t 32

z, p

e pt

0

z, t dt

2 Dt

the solution can be written as

Hence, if equation (4) is multiplied by e pt and integrated

2

z 2 2

term by term it is reduced to an ordinary differential equation

z

t

4D

e d. . (8)

d 2 p

2

2

d

d

2 Dt

dz2 D

(6)

Since t C0 exp 4D

(1 n) n

K t the

w2t

w2t

The solution of the equation (6) can be written as

A eqz B eqz

particular solution of the problem can be written as

where q p .

2

2

D

w2

z, t 2C 0 e 4 D

(1n)

n

Kd t

2

2

2

2

exp 2

2 d exp

d

0

0

(9)

2

e2 1

2 exp

• a

da

where

w2

4D

4D

(1 n)

n

K z

d 2 D

and

a

a

2

z e2 exp

• a

da .

.

a

2 Dt

Further, let

3. Evaluation of the integral solution

a

a

a

The integration of the first term of equation (9) gives

(Pierce, 1956)

in the first term of the above equation, then

2 2

2

(1 n)

I e2

e 2 d e2

exp a

da .

e 2 d

0

n Kd C

e2

2

1

a

a

For convenience the second integral can be expressed in terms of error function (Horenstein, 1945), because this function is well tabulated. Noting that

Similar evaluation of the second integral of equation (10) gives

2 2

2

2

2

2

2

2

2

I2 e exp

a da e2 exp

• a da

2

a a

2

Again substituting

• a into the first term, the

a

The second integral of equation (9) can be written as result is

2

2

I exp 2 d

2

I e2

e 2 d e2 exp a da .

0

1

2

2

2

a

a

e2 exp

d e2 exp

d

Noting that

2 0

(10)

0

2

2

exp a a

2 da exp a a 2 da

Since the method of reducing integral to a tabulated function is the same for both integrals in the right side of

equation (10), only the first term is considered. Let a

and adding and subtracting, we get

Substitute this into equation (10) gives

2

e2 exp a

da .

I e2

e 2 d e2

e 2 d

a

.

The integral can be expressed as

2

Thus, equation (9) can be expressed as

I e2 exp

d

2 (1n)

Kd t

0

2C0

4 D n

w

w

z, t e .

e2 1 e2

e 2 d e2

e 2 d (11)

However, by definition

e2 e 2 d e2 erfc

flexible than the standard form of other methods for one- dimensional transport model. Figures 1 to 4 represents the concentration profiles verses distance along the media for

also,

2

different values of porosity n. 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, whereas

e2

e 2 d e2

erf

concentration profile versus time for different values of

2

2

1

depth z. For a fixed z it is seen that concentration increases in the beginning due to lesser effect of dispersion

Writing equation (11) in terms of the error functions, we get

coefficient D and reaches a steady-state value for larger time.

w2 (1n)

Kd t

z, t C0 e 4D n 2 erfc e2 erfc

e

2

Substitute the value of z, t in equation (3) the solution reduces to

e erfc e erfc .

e erfc e erfc .

C 1 exp wz 2 2

C0 2

2D

(12)

Fig. 1: Break-through-curve for C/C0 v/s depth

Resubstituting the value of and gives

for n=0.5 and Kd=0.4

C 1 wz

w2

(1 n)

z

exp exp

Kd

4D

4D

C0 2 2D n

D

z w2t (1 n)

erfc K t

2

2

n

n

Dt 4D d

w2

(1 n)

z

d

d

exp K .

z

4D

w2t

n

(1 n)

D

Fig. 2: Break-through-curve for C/C0 v/s depth for n=0.5 and Kd=1.0

erfc

K t

(13)

2

2

n

n

Dt 4D d

where boundaries are symmetrical the solution of the problem is given by the first term of equation (13). The second term in equation (13) is thus due to the asymmetric boundary imposed in a general problem. However, it should be noted that if a point a great distance away from the source is considered, then it is possible to appropriate

the boundary conditions by C , t C0 , which leads

to a symmetrical solution.

4. Results & Discussions:

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

Fig. 3: Break-through-curve for C/C0 v/s depth for n=1.0 and Kd=0.4

Fig. 4: Break-through-curve for C/C0 v/s depth for n=1.0 and Kd=1.0

The figures represent C/C0 verses time for different values of distribution coefficient Kd. It is seen that for a fixed Kd, concentration increases slowly up to t=10 days because of the less adsorption of pollutants on the solid surface and then reaches a constant value for larger time where the effect of distribution coefficient Kd is small. We conclude that the integral transform method is a powerful method to derive analytical solutions for solute transport of a adsorption in homogeneous porous media and under different flow conditions. Steady-state concentration distributions and temporal moments can be directly derived from these solutions and transient concentration distribution is accessible via numerical inversion. The derived solutions are of great value for bench-marking numerical reactive transport codes.

REFERENCES:

1. Aral,M.M.,Liao,B.,1996.Analytical solutions for two-dimensional transport equation with time-dependant dispersion co-efficients. Journal of Hydrologic Engineering, 1,20-32.

2. Barry,D.A.,and Sporito,G.,1989. Analytical solution of a convection- dispersion model with time-dependant transport co-efficients. Water Resour.Res.,25,2407-2416.

3. Batu,V.,1993. A genaralized two-dimensional analytical solute transport model in bounded media for flux-type multiple sources. Water Resour.Res.,29,2881-2892.

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

5. Ermak,D.L.,1977.An Analytical Model for Air Pollutant transport and deposition from a point source.Atmos.Environ.,11,231-237.

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,1113-1124.

7. Koch,W.,1989.A Solution of two-dimensional atmosphere diffusion equation with height-dependent diffusion coefficient including ground level absorption. Atmos. Environ., 23, 1729-1732.

8. 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, 193-204.

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

10. Sudheendra S.R. 2012. An analytical solution of one-dimensional advection-diffusion equation in a porous media in presence of radioactive decay, Global Journal of Pure and Applied Mathematics, Vol.8, No. 2, 113-124.

11. 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, 32-40.

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

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

14. Tartakowsky,D.,Di Federico,V.,1997. An analytical solution for contaminant transport in non-uniform flow. Transport in porous media, 27, 85-97.

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