Contaminant Isolation By Cutoff Walls: Reconsideration Of Mass Fluxes

Similar documents
5. TWO-DIMENSIONAL FLOW OF WATER THROUGH SOILS 5.1 INTRODUCTION

Pollution. Elixir Pollution 97 (2016)

The epsilon method: analysis of seepage beneath an impervious dam with sheet pile on a layered soil

F. BEN ABDELGHANI a,b. Tunisie;

CALCULATION OF STEAM AND WATER RELATIVE PERMEABILITIES USING FIELD PRODUCTION DATA, WITH LABORATORY VERIFICATION

Soil Mechanics Permeability of Soils and Seepage page 1 CHAPITRE 9. PERMEABILITY OF SOILS AND SEEPAGE...1

Outils de Recherche Opérationnelle en Génie MTH Astuce de modélisation en Programmation Linéaire

CENG 5210 Advanced Separation Processes. Reverse osmosis

NUMERICAL ANALYSIS OF THE EFFECTS OF STREAMLINING GEOMETRY AND A VECTOR WALL ON THE THERMAL AND FLUID FLOW IN A SRU THERMAL REACTOR.

Vapor Pressure Prediction for Stacked-Chip Packages in Reflow by Convection-Diffusion Model

Computation of turbulent natural convection at vertical walls using new wall functions

Numerical Treatment of Two-Phase Flow in Porous Media Including Specific Interfacial Area

External Flows. Dye streak. turbulent. laminar transition

A New Method for Calculating Oil-Water Relative Permeabilities with Consideration of Capillary Pressure

RADIONUCLIDE DIFFUSION IN GEOLOGICAL MEDIA

18.12 FORCED-DAMPED VIBRATIONS

Heat Transfer Analysis of a Space Radiating Fin with Variable Thermal Conductivity

On Lithium Wall and Performance of Magnetic Fusion Device

T Fluid temperature in the free stream. T m Mean fluid temperature. α Thermal diffusivity. β * Coefficient of concentration expansion

TEMPORALLY DEPENDENT DISPERSION THROUGH SEMI-INFINITE HOMOGENEOUS POROUS MEDIA: AN ANALYTICAL SOLUTION

Towards a Theory of Societal Co-Evolution: Individualism versus Collectivism

PET467E-Analysis of Well Pressure Tests 2008 Spring/İTÜ HW No. 5 Solutions

Effect of Insertion Devices. Effect of IDs on beam dynamics

Evaluation of an Oxygen Injection Technology for In-Situ Hydrocarbon Bioremediation in a Fractured Bedrock Environment

DETERMINING HIGH VOLTAGE CABLE CONDUCTOR TEMPERATURES. Guy Van der Veken. Euromold, Belgium. INVESTIGATIONS. INTRODUCTION.

Effects of Asymmetric Distributions on Roadway Luminance

Exploring the relationship between a fluid container s geometry and when it will balance on edge

BAGNOLD FORMULA REVISITED: INCORPORATING PRESSURE GRADIENT INTO ENERGETICS MODELS

EFFECTS OF SCOUR AND HYDRAULIC GRADIENT ON THE STABILITY OF GRANULAR SOIL SLOPE

Reynolds Averaging. We separate the dynamical fields into slowly varying mean fields and rapidly varying turbulent components.

Estimating vertical velo city on the Equator

Analytical solutions for water flow and solute transport in the unsaturated zone

On the direct kinematics of planar parallel manipulators: special architectures and number of solutions

A turbulence closure based on the maximum entropy method

UNCERTAINTY SCOPE OF THE FORCE CALIBRATION MACHINES. A. Sawla Physikalisch-Technische Bundesanstalt Bundesallee 100, D Braunschweig, Germany

HYDRODYNAMIC STABILITY ANALYSIS OF SHEARED CONVECTIVE BOUNDARY LAYER FLOWS IN STRATIFIED ENVIRONMENTS

MASSACHUSETTS INSTITUTE OF TECHNOLOGY SPRING 2007

O.R. Jimoh, M.Tech. Department of Mathematics/Statistics, Federal University of Technology, PMB 65, Minna, Niger State, Nigeria.

Why do Golf Balls have Dimples on Their Surfaces?

u Velocity in x-direction Velocity in y-direction g- Acceleration due to gravity c - Specific heat at constant presure

1D Verification Examples

Incremental identification of transport phenomena in wavy films

The Comparison between the Effects of Using Two Plane Mirrors Concentrator and that without Mirror on the Flat- Plate Collector

Second-Order Linear ODEs (Textbook, Chap 2)

Department of Mathematic, Ganjdundwara (P.G.) College, Ganjdundwara (Kashiram Nagar) (U.P.)

MEASUREMENTS OF TIME-SPACE DISTRIBUTION OF CONVECTIVE HEAT TRANSFER TO AIR USING A THIN CONDUCTIVE-FILM

Vidyasagar et al., International Journal of Advanced Engineering Technology E-ISSN A.P., India.

Comparison of Heat and Mass Transport at the Micro-Scale

LOW WATER PERMEABILITY MEASUREMENTS OF CLAY SAMPLE. CONTRIBUTION OF STEADY STATE METHOD COMPARED TO TRANSIENT METHODS

Integration of Different Types of Data for Characterization of Reservoir Heterogeneity

INTERACTION OF THE ELECTROMAGNETIC RADIATION WITH THE SURFACE OF PALLADIUM HYDRIDE CATHODES

AN ANALYSIS OF THE LAWS OF DRAINAGE NET COMPOSITION

CHARACTERIZATION OF ULTRASONIC IMMERSION TRANSDUCERS

Data analysis from capillary rheometry can be enhanced by a method that is an alternative to the Rabinowitsch correction

MATRIX-FRACTURE TRANSFER FUNCTIONS FOR PARTIALLY AND COMPLETELY IMMERSED FRACTURES

CHAPTER 4: DYNAMICS: FORCE AND NEWTON S LAWS OF MOTION

Possible Cosmic Influences on the 1966 Tashkent Earthquake and its Largest Aftershocks

Jurassic earthquake sequence recorded by multiple generations of sand blows, Zion National Park, Utah

Transactions on Ecology and the Environment vol 7, 1995 WIT Press, ISSN

( ) 2 ( kg) ( 9.80 m/s 2

Accurate and Estimation Methods for Frequency Response Calculations of Hydroelectric Power Plant

A Boundary Condition for Porous Electrodes

A DIFFERENT APPROACH TO MULTIPLE CORRESPONDENCE ANALYSIS (MCA) THAN THAT OF SPECIFIC MCA. Odysseas E. MOSCHIDIS 1

Multipole expansion of strongly focussed laser beams

Time-Dependent Conduction :

ème Congrès annuel, Section technique, ATPPC th Annual Meeting, PAPTAC

Analytical Study of Airflow Induced by Barometric Pressure and Groundwater Head Fluctuations in a Two-Layered Unsaturated Zone

Sediment yield and availability for two reservoir drainage basins in central Luzon, Philippines

Deep Borehole Disposal Performance Assessment and Criteria for Site Selection

Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.

Poisson s ratio effect of slope stability calculations

Modeling Seepage Control in Hydraulic Structures

The method by which this error has to be combined with other sources of error is I

Safety assessment for disposal of hazardous waste in abandoned underground mines

J. C. Sprott and A. C. England

CHAPTER 3 THE COMMON FACTOR MODEL IN THE POPULATION. From Exploratory Factor Analysis Ledyard R Tucker and Robert C. MacCallum

Capillary rise between closely spaced plates : effect of Van der Waals forces

Pollution des sols et interaction avec la biosphère

Pile Integrity Testing Developments

Fundamentals of DC Testing

Frozen Ground Containment Barrier

Analysis of fluid induced vibration of cryogenic pipes in consideration of the cooling effect

Asymptotics for kink propagation in the discrete Sine-Gordon equation

Solving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions

Erosion and Gas Content in Subsoil

Basic Fluid Mechanics

PLASMA PHYSICS. 1. Charge particle motion in fields of force 2. Plasma particles and their interaction 3. Characteristic propriety of plasma

470 Lecture #7 of 18

Multi-dimensional electro-osmosis consolidation of clays

ROBUST LINEAR DISCRIMINANT ANALYSIS WITH A LAPLACIAN ASSUMPTION ON PROJECTION DISTRIBUTION

DNAPL migration through interbedded clay-sand sequences

ABSTRACT INTRODUCTION

COS 511: Theoretical Machine Learning. Lecturer: Rob Schapire Lecture # 12 Scribe: Indraneel Mukherjee March 12, 2008

Natural Heat Transfer Convection in a Square Cavity Including a Square Heater

Drift & Anisotropy Parameters in Sporadic-E & Normal-E Layers"

REVUE FRANÇAISE D INFORMATIQUE ET DE

A new formulation of immiscible compressible two-phase flow in porous media

6. Continuous Release - Point Source

Thermal Analysis Contents - 1

Plastic Equilibrium in Soil

Transcription:

GeoHalifax29/GéoHalifax29 Contaminant Isolation By Cutoff Walls: Reconsideration Of Mass Fluxes Christopher. Neville S.S. Papadopulos & ssociates, Inc., Waterloo, Ontario BSTRCT Cutoff alls are used frequently to isolate contaminants at both controlled and uncontrolled hazardous aste sites. Neville and ndres (26) presented a containment criterion for contaminant isolation by a cutoff all. Their analysis yields the arcy flux required to achieve containment, based on the condition that long-term advective and diffusive mass fluxes across the all are balanced. In this paper e sho that the condition of zero net mass flux represents only one particular case. Straightforard expressions for the long-term mass fluxes across a cutoff all can also be developed from the same theory for any arcy flux. The expressions for the long-term mass fluxes may be used to estimate the mass flux to the environment in cases here it is difficult to satisfy the criterion of zero net mass flux. RÉSUMÉ es barrières souterraines imperméables sont souvent utilisées pour confiner les sources de contamination présentes sur les terrains contaminés. Neville et ndres (26) ont présenté, pour de telles barrières, un critère de confinement des sources de contamination. Le résultat de leur analyse indique quel est le flux de arcy permettant d atteindre le confinement, en se basant sur l équilibre, à long-terme, entre les flux de masse advectifs et dispersifs à travers la barrière. Cet article démontre que la condition d équilibre, où le flux de masse net de contamination traversant la barrière est nul, représente un cas particulier. es simples expressions mathématiques décrivant les fluxes de masse à travers une barrière peuvent êtres obtenues pour toutes les valeurs de flux de arcy. es expressions quantifiant les débits massiques à long-terme peuvent être utilisée pour estimer la charge de contamination relâchée dans l environnement souterrain, dans les cas où il est difficile d obtenir un parfait confinement. 1 INTROUCTION Cutoff alls are used frequently to isolate contaminants at both controlled and uncontrolled hazardous aste sites. Neville and ndres (26) presented a containment criterion for contaminant isolation by a cutoff all. Their containment criterion builds on the analysis of evlin and Parker (1996) and is developed from the assumption that containment is achieved hen the long-term advective and diffusive mass fluxes across the all are balanced so that the net mass flux is zero. The conceptual model of Neville and ndres (26) is shon schematically in Figure 1. Their analysis may be used to estimate the arcy flux required to achieve containment according to this criterion. In this paper, e sho that the assumption of zero net mass flux represents only one particular case. The same theory may be used to develop a straightforard expression for the long-term net mass flux across a cutoff all for any arcy flux. This expression may be used to estimate the net mass flux to the environment in cases here it may not be possible to satisfy the criterion of zero net mass flux. Figure 1. Conceptual model of source isolation 875

GeoHalifax29/GéoHalifax29 2 NEVILLE N NREWS (26) CONTINMENT CRITERION The analysis of Neville and ndres (26) assumes a constant concentration ithin the source zone and steady one-dimensional flo and dispersive diffusive transport across the cutoff all. In Figure 1, c and c denote the steady concentrations along the inside and outside faces of the all, respectively, and is the thickness of the all. The terms, diff, and mech represent the advective, diffusive, and mechanical-dispersive mass fluxes, respectively. The advective mass flux,, is thus defined: = qc [1] here c is the concentration and q is the arcy flux (positive outards from the source zone). The steadystate dispersive diffusive mass flux,, is given by: dc = θ [2] dx Given a knon concentration of the source and a target concentration of outside of the cutoff all, the Neville- ndres containment criterion yields the arcy flux that is required to achieve a long-term net mass flux of zero: θ c q = ln c 3 GENERL INTERPRETTION OF THE NEVILLE N NREWS (26) NLYSIS The analysis of Neville and ndres (26) accommodates groundater seepage in either direction, and any combination of specified concentrations over both faces of the cutoff all. Hoever, Equation 5 represents only one specific case, in hich the dispersive diffusive and advective mass fluxes are balanced. It is important to note that steadystate conditions may be attained for any arcy flux. To demonstrate this point, e consider the more general case of transient conditions. The governing equation for transient transport is: [5] here θ is the effective porosity and is the dispersion coefficient. The diffusive flux in Equation 2 is a lumped representation of the mechanical-dispersive and diffusive mass fluxes. Both process are assumed to be Fickian processes and additive; therefore, the dispersion coefficient can be interpreted as the sum of the mechanical-dispersive and diffusive fluxes (Bear, 1972): θ = q + θ t x x 2 c c c The initial and boundary conditions are: 2, x [6] c ( x, ) = [7a] α v * = L + [3] c(, t) c = [7b] here α L is the longitudinal dispersivity, v is the groundater velocity, and * is the effective diffusion coefficient. The mechanical-dispersive flux accounts for variations in the groundater velocity across the all that are beneath the scale of resolution of the analysis. For a properly constructed all, this flux should be negligible. Neville and ndres (26) derived a containment criterion for any given all design (all thickness, diffusion coefficient, and arcy flux across the all). Folloing the approach of evlin and Parker (1996), they assumed that containment is achieved hen the net mass flux on the outside face of the all is zero: c(, t) c = [7c] + = [4] 876

GeoHalifax29/GéoHalifax29 The general solution for transient conditions is derived by generalizing the derivation of l-niami and Rushton (1977): qx sinh 2θ (, ) = q sinh c x t c n= 1 2θ q( x ) { } 2 2π q x q t + c 2 2 2θ 4θ 1 n n π t nπ x sin n 2 2 2 2 2θ 2 2 n π q + 2 2 4θ The solution for the steady-state concentration profile can be derived directly from Equation 8 and is given in Neville and ndres (26) as: [8] ifferentiating the general steady-state solution yields: qx dc ( x, ) qv θ = c dx θ q qx q θ c θ q [12] Evaluating the derivative at x = and simplifying yields: q dc(, ) q θ ( c c [13] = ) dx θ q (, ) = ( ) c x c c c qx θ q θ 4 STEY-STTE NET MSS FLUX [9] The expression for the diffusive flux is obtained by substituting Equation 13 into 11: = q c c q θ q [14] The Neville and ndre (26) containment analysis is generalized by deriving expressions for the advective and dispersive diffusive mass fluxes for any magnitude of the arcy flux across the all. The advective mass flux across the outside face of the all is: The net mass flux is defined as the sum of the advective and diffusive mass fluxes: net = + [15] = qc = qc [1] x= The steady-state dispersive diffusive mass flux across the outside face of the all is defined as: dc = θ [11] dx x= Substituting Equations 1 and 14 into Equation [15] yields the expression for the steady-state net mass flux across the outside face of the cutoff all: = qc q c c net q θ q [16] 877

GeoHalifax29/GéoHalifax29 5 EXMPLE CLCULTIONS To illustrate the characteristics of the analysis, transient and steady-state concentration profiles across the all are calculated for to cases. The parameters for the calculations are listed in Table 1. Table 1. Parameters for example calculations Parameter Value * 2.592 1-5 m 2 /d α L. m θ.3 1. m c 1. c.1 Case 1: -6 1-5 m/d (inard flo) q Case 2: +6 1-5 m/d (outard flo) The results for Case 1 are shon in Figure 2. s shon in Figure 2, groundater flo directed inards across the all gives rise to a steady-state concentration profile that is concave-inards. s shon in Figure 3, groundater flo that is directed outards across the all gives rise to a steady-state concentration profile that is concave-outards. The length of time required for steady concentrations to evolve differs for the to cases, but a steady profile is eventually attained in both cases. For Case 1, the inards flo opposes the concentration gradient. The steady-state net mass flux is given by: net ( 6 1 / )(.1) ( 6 1 / )(( 1.) (.1) ) = m d m d ( 6 1 m / d )( 1. m) 2 (.3)( 2.592 1 m / d ) ( 6 1 m / d )( 1. m) 2 (.3)( 2.592 1 m / d ) 7 8 7 = + = [17] 6 1 2.65 1 5.74 1 In this case, the net mass flux is negative, indicating that the inard advective mass flux exceeds the outard diffusive mass flux. The results for Case 2 are shon in Figure 3 In this case, the net mass flux is directed inards across the all. For Case 2, the outard flo is in the same direction as the concentration gradient. The steady-state net mass flux is given by: net ( 6 1 / )(.1) ( 6 1 / )(( 1.) (.1) ) = m d m d ( 6 1 m / d )( 1. m) 2 (.3)( 2.592 1 m / d ) ( 6 1 m / d )( 1. m) 2 (.3)( 2.592 1 m / d ) 7 = + = + [18] 6 1 5.94 1 6. 1 In the second case, diffusion and advection both direct solute toards the outside of the all. The solution for Case 2 does not predict that the entire all attains a concentration of the source c. Rather, solute is dran off continuously at x = to maintain the concentration at the fixed level c. s the outards advective flux increases, the concentration across the all becomes nearly uniform at c and declines abruptly to c at the boundary. The net mass flux for a range of arcy fluxes is plotted in Figure 4. The condition for zero net mass flux is also indicated. These results reinforce the point made previously; the condition of zero net mass flux is only one state along a continuous spectrum of possibilities. If the aquifer outside of the all is extensive, it is reasonable to assume significant dilution of the solute that migrates out of the all. The orst case ith respect to the predicted outard mass flux ill arise hen the outside concentration c is assumed to be zero. In this case, there is no advective flux to counteract the outard diffusion of solute and the expression for the net mass flux reduces to: net q θ = qc q [19] Results ith Equation 2 for a range of values of c are plotted in Figure 5. For reference purposes, results are included from the previous example c /c =.1. To general observations can be made regarding the results. When the arcy flux is directed inards (negative), the net mass flux is only eakly sensitive to the outside concentration. When the arcy flux is directed outards, the net mass flux is essentially independent of the outside concentration. The assumption of an outside concentration c of zero provides an appropriate preliminary estimate of the net mass flux. The solution given by Equation 15 is straightforard in interpretation and implementation, and yields conservative estimates of the net mass flux. 878

GeoHalifax29/GéoHalifax29 1.9.8 Relative concentration, c/c.7.6.5.4.3.2 5, d.1 t = 1, d Steady-state (dashed).1.2.3.4.5.6.7.8.9 1 istance across all, x (m) Figure 2. Concentrations for Case 1, inards flo 1.9 Steady-state.8 5, d Relative concentration, c/c.7.6.5.4.3 t = 1, d 1, d.2.1.1.2.3.4.5.6.7.8.9 1 istance across all, x (m) Figure 3. Concentration for Case 2, outards flo 879

GeoHalifax29/GéoHalifax29.15 Net mass flux net at x =.1.5 Zero net mass flux at q = -3.58 1-5 m/s -.5 -.2 -.1.1.2 arcy flux across all, q (m/d) Figure 4. Net mass flux outside of the all.15 Net mass flux net at x =.1.5 -.5 c /c =.1 c /c =.5 c /c =.1 c /c =. -.2 -.1.1.2 arcy flux across all, q (m/d) Figure 5. Net mass flux outside of the all: Effect of c 88

GeoHalifax29/GéoHalifax29 6 CONCLUSIONS The containment criterion of Neville and ndres (26) is based on the condition that the net mass flux across the cutoff all is zero. steady-state concentration profile ill develop across the all for any magnitude and direction of the arcy flux, and the condition of zero net mass flux may be unnecessarily restrictive to achieve the containment objectives at a particular site. This supplement to the analysis of Neville and ndres (26) presents an expression for the steadystate net mass flux across a cutoff all for any value of the arcy flux. The results of example calculations reveal that hen the arcy flux is directed inards toards the source zone, the net mass flux is eakly dependent upon the concentration on the outside face of the all. When the arcy flux is directed outards from the source zone, the net mass flux is essentially independent of the outside concentration. The assumption of an outside concentration equal to zero yields a simple expression for the net mass flux that is straightforard to interpret and evaluate, yielding a conservative estimate of the net mass flux across the all. 7 REFERENCES l-niami,.n.s. and Rushton, K.R. 1977. nalysis of Flo gainst ispersion in Porous Media, ournal of Hydrology, 33: 87-97. Bear,. 1972. ynamics of Fluids in Porous Media, merican Elsevier Publishing Company, Inc. Ne York, NY, US. evlin,.f., and Parker, B.L. 1996. Optimum Hydraulic Conductivity to Limit Contaminant Flux through Cutoff Walls, Ground Water, 34(4): 719 726. Neville, C. and ndres, C.B. 26. Containment Criterion for Contaminant Isolation by Cutoff Walls, Ground Water, 44(5): 682-686. 881