Analytical Modeling of Nonaqueous Phase Liquid Dissolution with Green's Functions

Equilibrium and bicontinuum nonequilibrium formulations of the advection–dispersion equation (ADE) have been widely used to describe subsurface solute transport. The Green's Function Method (GFM) is particularly attractive to solve the ADE because of its flexibility to deal with arbitrary initial and boundary conditions, and its relative simplicity to formulate solutions for multi‐dimensional problems. The Green's functions that are presented can be used for a wide range of problems involving equilibrium and nonequilibrium transport in semi‐infinite and infinite media. The GFM is applied to analytically model multi‐dimensional transport from persistent solute sources typical of nonaqueous phase liquids (NAPLs). Specific solutions are derived for transport from a rectangular source (parallel to the flow direction) of persistent contamination using first‐, second‐, or third‐type boundary or source input conditions. Away from the source, the first‐ and third‐type condition cannot be expected to represent the exact surface condition. The second‐type condition has the disadvantage that the diffusive flux from the source needs to be specified a priori. Near the source, the third‐type condition appears most suitable to model NAPL dissolution into the medium. The solute flux from the pool, and hence the concentration in the medium, depends strongly on the mass transfer coefficient. For all conditions, the concentration profiles indicate that nonequilibrium conditions tend to reduce the maximum solute concentration and the total amount of solute that enters the porous medium from the source. On the other hand, during nonequilibrium transport the solute may spread over a larger area of the medium compared to equilibrium transport.

[1]  F. Smedt.  Analytical solutions of the one-dimensional convective-dispersive solute transport equation: M.Th. van Genuchten and W.J. Alves. Technical bulletin no. 1661, U.S. Department of Agriculture, Washington, DC, 1982. 151 pp , 1984 .

[2]  J. C. Jaeger,et al.  Conduction of Heat in Solids , 1952 .

[3]  E. Voudrias,et al.  Modeling of contaminant transport resulting from dissolution of nonaqueous phase liquid pools in saturated porous media , 1994 .

[4]  D. Adrian,et al.  Analytical Solution of the Convective-Dispersive Equation for Cation Adsorption in Soils 1 , 1973 .

[5]  G. Roach,et al.  Green's Functions , 1982 .

[6]  M. V. Genuchten,et al.  Mass Transfer Studies in Sorbing Porous Media: II. Experimental Evaluation with Tritium (3H2O)1 , 1977 .

[7]  K. Rathfelder,et al.  Flow and entrapment of dense nonaqueous phase liquids in physically and chemically heterogeneous aquifer formations , 1998 .

[8]  John W. Miles,et al.  Application of Green's functions in science and engineering , 1971 .

[9]  I. Stakgold.  Green's Functions and Boundary Value Problems , 1979 .

[10]  M. V. Genuchten,et al.  Mass transfer studies in sorbing porous media. I. Analytical solutions , 1976 .

[11]  W. J. Alves,et al.  Analytical solutions of the one-dimensional convective-dispersive solute transport equation , 1982 .

[12]  Feike J. Leij,et al.  Analytical solutions for non-equilibrium solute transport in three-dimensional porous media , 1993 .

[13]  Constantinos V. Chrysikopoulos,et al.  Three‐Dimensional Analytical Models of Contaminant Transport From Nonaqueous Phase Liquid Pool Dissolution in Saturated Subsurface Formations , 1995 .

[14]  G. Walker.  Solution to a Class of Coupled Linear Partial Differential Equations , 1987 .

[15]  A. Prakash.  Ground-water contamination due to transient sources of pollution , 1984 .

[16]  G. T. Yeh.  AT123D: analytical transient one-, two-, and three-dimensional simulation of waste transport in the aquifer system , 1981 .

[17]  K. H. Coats,et al.  Dead-End Pore Volume and Dispersion in Porous Media , 1964 .

[18]  I. Javandel,et al.  Analytical solutions for solute transport in a vertical aquifer section , 1997 .

[19]  M. Th. van Genuchten,et al.  Two‐Site/Two‐Region Models for Pesticide Transport and Degradation: Theoretical Development and Analytical Solutions , 1989 .

[20]  A. Haji-sheikh,et al.  Heat Conduction Using Green's Function , 1992 .

[21]  Donald P. Galya,et al.  A Horizontal Plane Source Model for Ground‐Water Transport , 1987 .

[22]  H. Holman,et al.  Evaluation of Transient Dissolution of Slightly Water-Soluble Compounds From a Light Nonaqueous Phase Liquid Pool , 1996 .

[23]  Gour-Tsyh Yeh,et al.  Analytical three‐dimensional transient modeling of effluent discharges , 1976 .

[24]  Murray R. Spiegel,et al.  Schaum's outline of theory and problems of laplace transforms , 1965 .

[25]  F. Molz,et al.  An Analysis of Dispersion in a Stratified Aquifer , 1984 .

[26]  R. Rogers,et al.  An introduction to partial differential equations , 1993 .

[27]  P. J. Wierenga,et al.  A generalized solution for solute flow in soils with mobile and immobile water , 1979 .

[28]  Marcel G. Schaap,et al.  Solute transport modeled with Green's functions with application to persistent solute sources , 2000 .

[29]  T. Ellsworth,et al.  Three‐dimensional analytical solutions to the advection‐dispersion equation in arbitrary Cartesian coordinates , 1993 .

[30]  M. V. Genuchten,et al.  A comprehensive set of analytical solutions for nonequilibrium solute transport with first‐order decay and zero‐order production , 1993 .

[31]  G. Folland.  Introduction to Partial Differential Equations , 1976 .