Analytical solutions for water flow and solute transport in the unsaturated zone
|
|
- Willa Cross
- 5 years ago
- Views:
Transcription
1 Models for Assessing and Monitoring Groundwater Quality (Procsedines of a Boulder Symposium July 1995). IAHS Publ. no. 227, Analytical solutions for water flow and solute transport in the unsaturated zone M. H. NACHABE, A. L. ISLAS & T. H. ILLANGASEKARE University of Colorado at Boulder, Campus Box 428, Boulder, Colorado, 80309, USA Abstract Analytical solutions for the water flow and solute transport equations in the unsaturated zone are presented. We use the Broadbridge & White non-linear model to solve the Richards' equation for vertical flow under a constant infiltration rate. Then we extend the water flow solution and develop an exact parametric solution for the advectiondispersion equation. The location of a solute front is determined with the method of characteristics. The dispersion component is incorporated into the final solution using a singular perturbation method. The analytical solutions are simple and can be generated without resorting to computationally demanding numerical schemes. Indeed, the analytical solutions provide tools to verify the accuracy of numerical models. Comparison with a finite element numerical solution indicates that a good match for the predicted water content is achieved when the mesh grid is one fourth the capillary length scale of the porous medium. However, numerical dispersion and spatial oscillations were significant in solving the solute transport equation at this level of discretization. INTRODUCTION The movement of water and solutes through the unsaturated zone has been of importance in traditional applications of ground water hydrology, soil physics, and agronomy. In recent years, the need to understand the behavior of hazardous waste and toxic chemicals in soils has resulted in a renewed interest in this subject. One of the primary concerns is that dissolved contaminants may migrate through the unsaturated zone, reach the saturated zone, and contaminate the ground water. The search for analytical solutions to model flow and solute transport continues to be of scientific interest. Typically, water flow and solute transport in unsaturated soils result in transient phenomena, making it a challenging problem. The nature of soil hydraulic properties like diffusivity and hydraulic conductivity renders the governing flow equation non-linear. In this paper, the Broadbridge & White (1988) (B&W, hereafter) non-linear flow model for constant infiltration rate is adopted. The ability of this model to conform to a variety of soil types, as demonstrated by White & Broadbridge (1988) and White & Perroux (1989), is very encouraging in field applications The objective of this paper is to use B&W analytical flow model to simulate solute movement in soils so that a closed form analytical solution for solute
2 126 M. H. Nachabe et al. transport in the unsaturated zone is available. The analytical solutions are compared with a numerical model to assess the accuracy of the numerical approximations. THEORY Two fundamental equations govern water flow and solute transport in unsaturated soils. Richards' equation is obtained by combining Darcy's law for unsaturated flow with the continuity equation. For one dimensional vertical flow, this equation is dt dz Z>(0), 30" -#-Jf(6) (1) dzj dz The transport of a conservative solute (i.e. non-volatile, non-sorbing, non-degrading) in the unsaturated zone is described by the advection-dispersion equation written in the form -Be -(0C) = ^{qc)+^\ EÏf (2) d dt 3z dzl dz In both equations, t is time, z is depth taken zero at the soil surface and positive downward, Q(z,t) is volumetric soil water content, Z>(0) is soil water diffusivity, K(Q) is unsaturated soil hydraulic conductivity, c(z,t) is solute concentration (mass of solute over volume of solution), q(z,t) is Darcy's water flux q(z,t) =K(B)-D(Q) jr (3) dz and is a hydrodynamic dispersion coefficient, a parameter that lumps the impacts of molecular diffusion and mechanical mixing. Equations (1), (2) and (3) aie coupled and they should be solved sequentially. First, the flow equation is solved for 0(z,r) and the subsequent substitution of Q(z,t) into (3) provides the unsaturated Darcy's flux q(z,t). The solutions of Q(z,t) and q(z,t) become input for (2). Equation (2) is reduced to a more practical mathematical form using the continuity equation for water flow 38 dq n Substitution of (4) into (2) yields 9r q dz + dz^t'dz ) {5) Equation (5) applies for transient transport of conservative solutes in unsaturated soils. ANALYTICAL SOLUTIONS The water content Q(z,t) solution of (1) was determined by B&W (1988) for a constant infiltration rate. This solution is presented here for the sake of completeness. The dimensionless form of Richard's equation is
3 Analytical solutions for water flow and solute transport in the unsaturated zone 127 =À- 9r* dz D*(ej^oz* 3 < 6 > oz* where z* = zfk s and f* = t/t s. X s and f* are the capillary length and time scales of the porous medium (Broadbridge & White, 1988; Nachabe et al., 1994; Nachabe & Illangasekare, 1994). Similarly the initial and boundary conditions are expressed in dimensionless form as 0* = 0 r* = 0 z* > 0 (7) K* (0 ) - >*(e*);r- = <7*n z*=0,r*>0 (8) u oz* where q* 0 = (q 0 - K 0 )/AK. The solution for (6) subject to the conditions in (7) and (8) is (B&W, 1988) 6. =CI1 = I (9) 2p+ 1 - u «K z* = C _1 p 2 ( 1 + P _1 )T:+p(2 + p"' 1 )Ç-ln(w) (10) The analytical expressions in (9) and (10) are documented in B&W (1988). Equations (9) and (10) constitute an exact parametric solution with parameter Ç. PARAMETRIC SOLUTION FOR THE ADVECTION-DISPERSION EQUATION We capitalize on the analytical form for the water flow solution and pursue an extension of this result to the movement of a conservative solute in the unsaturated zone. This is achieved by coupling the advection-dispersion equation with the water content solution in (9) and (10). Wilson & Gelhar (1981) developed analytical solutions for the advection-dispersion equation, but as opposed to the result presented here, their solution for water flow in the soil was an adaptation of the Philip & Knight (1974) iterative procedure, which does not provide exact parametric solution. We assume that initially there is zero concentration of solute in the soil profile. A step of solute concentration c 0 is applied at the surface. By introducing the dimensionless variables c* = clc 0 and E* = E[t s /k s 2 ] into (5), the advection dispersion equation is expressed in the mathematical form 9 n"i 9c * 7) (( 0 n"l dc *\ K n dc e * - +^k = ^ll e - + Â5r5^J- (^+^)3i: (11) The dimensionless initial and boundary conditions are written as c* = l z = 0 f>0 (12) c*=0 r* = 0 z*>l (13)
4 128 M. H. Nachabe et al. The solution procedure to derive a parametric solution for (1) is similar to the one adopted by Wilson & Gelhar (1981). This solution procedure solves sequentially for the two mechanisms of solute transport, i.e. convection with the mean water velocity and dispersive mixing. The convection component is treated by the method of characteristics to determine the location of the solute front. The dispersive component is then superimposed on the resulting front using the method of singular perturbation. The convective component is obtained by setting the dispersive in (11) equal to zero with the result 0 \3c* K dc 9 \ n, ^Âëfe+^ + Âf^0 (14) Thus we find that the characteristic r\ is given by the expression (see Nachabe et al., 199 A) ^(à + 1 ) Cz *" c-( â +^o)^ (is) which using the parametric solution in (10) for z* can be written in terms of Ç and the time x as r) = A 1 t: + S 1 Ç-C 1 lnm(ç,x) (16) where A l5 B x and C\ are constants defined symbolically as ^1 = C I(P 2 + P ) - ( ^ 1 2 ) (17) B x = Cj(2p+1) -1 (18) c i = 1 + cfe (19) This solution for T is a parametric expression in Ç. The mathematical procedures to locate the position zj(u) of the sharp front include two steps: (1) set TI equal to zero and determine tf, and (2) set (15) equal to zero and solve for zy in terms of dimensionless time t*. PERTURBATION APPROXIMATION OF DISPERSION Zj(t*) is the depth of the solute sharp front. The solute concentration c*(z*,f*) is one beyond this front and zero behind it. This contrast at the front boundary produces a singular problem that is tackled mathematically using the method of asymptotic expansion. For E* «1, an inner expansion valid in the vicinity of the front is introduced by stretching the coordinate system and rewriting the dispersion equation in the new coordinate system. The solution is (Nachabe et al., 1994)
5 Analytical solutions for water flow and solute transport in the unsaturated zone 129 c* C z *> t*) 2 erfc 2 * 1/2 (z* 9 (20) ILLUSTRATION The solutions for water flow and solute transport in (9), (10), (16) and (20) area analytic. So, the water content profile 8*(z*,r*) and the concentration of solute in the soil c*(z*,t*) are determined without iterative steps commonly used in numerical schemes. The input requirement for analytical simulation includes: (1) soil hydraulic properties: the saturated conductivity K s, the saturated water content Q s, the parameter C, and the capillary length scale X s ; (2) the soil initial moisture content 0 and the water flux boundary condition q* 0 and (3)the dispersion coefficient E*. For the sake of illustration, we choose C=1.2, X S =3M cm, K s =4.5 cm h" 1, q* 0 =0.2, 6^=0.4, 6 =0.08, *=10-3. The values of these parameters reflected the conditions at a site in Golden, Colorado. t s is calculated to be 0.27 h. Numerical simulations with the PRINCETON finite element model (Celia, 1991; Celia et al, 1990) were conducted for the same input conditions. The time step in the PRINCETON code is adjusted internally depending on the number of nonlinear iterations (Celia, 1991). For the finite element numerical scheme, we chose an element size of one fourth the capillary length scale (i.e. V 4 ) or 0.96 cm for this particular case. Figure 1 suggests a good match between the analytical and numerical solutions for the normalized water content 0* at t* equal 2, 10 and 24 at this level of discretization. However, spatial oscillations (overshot and undershoot) and numerical dispersion were significantly large for the advection dispersion equation. Figure 2 u - t* =2 ^ ~ - JS If -2 - * -4 - /7-6 - ' t* = I0 ^- t' = 3^ A/ AS If / / / / h i e. Fig. 1 Analytical (solid line) and numerical (dashed line) solutions of 9*.
6 130 M. H. Nachabe et al. -Z» f = 10-3 t* = 2i Fig. 2 Analytical (solid line) and numerical (dashed line) solutions of c*. C. compares the analytical and numerical solutions when we reduced the mesh grid in the finite element code to one eighth the capillary length scale. While the match between the concentrations is not very satisfactory, especially at t* = 24, the trade-off between numerical dispersion and spatial oscillations are typical for advection dominated solute transport like the case considered here. Indeed, Pinder & Huyakorn (1983, p. 207) cautioned "that because advection-dominated transport problems are difficult to solve numerically, the user of computer models needs to test the adequacy of selected mesh and time step discretizations to avoid obtaining misleading or meaningless solutions. For this reason we recommend the use of available analytical solutions for calibration..." CONCLUSION We extended the Broadbridge and White infiltration model to derive a closed form analytical solution for the advection-dispersion of a conservative solute in the unsaturated zone. The presented analytical solutions are exact and parametric, and the soil water content and the solute concentration solutions are easily generated and programmed. The computation time required to achieve a complete solution is insignificant when compared with iterative numerical schemes. The numerical problems involved in solving advection dominated transport are well recognized among practitioners. The developed analytical solutions should provide a mean for evaluating the accuracy of existing numerical codes.
7 REFERENCES Analytical solutions for water flow and solute transport in the unsaturated zone 131 Broadbridge, P. & White, I. (1988) Constant rate rainfall infiltration: a versatile nonlinear model, 1. Analytical Solution. Wat. Resour. Res. 24, Celia, M. A. (1991) The one-dimensional Princeton unsaturated code. Distributed by Water Resources Program, Princeton, University, Princeton, New Jersey. Celia, M. A., Bouloutas, E. T. & Zabra, R. L. (1990) A general mass-conservative numerical solution for the unsaturated flow equation. Wat. Resour. Res. 26, Huyakom, P. S. & Pinder, G. F. (1983) Computational Methods in Subsurface Flow. Academic Press, Inc. Nachabe, M. H. & Illangasekare, T.H. (1994) Use of tension infiltrometer data with unsaturated hydraulic conductivity models. Ground Water 32, Nachabe, M. H., Islas, A. L. & Illangasekare, T. H. (1995) Analytical solutions for water flow and solute transport in the unsaturated zone. Ground Water 33 (in press). Philip, J. R. & Knight, J. G. (1974) On solving the unsaturated flow equation, 3. New quasi-analytical technique. Soil Sci. 117,1-13. Wilson, J. L. & Gelhar, L. W. (1981) Analysis of longitudinal dispersion in unsaturated flow, 1. The analytical method. Wat. Resour. Res. 17, White, I. & Broadbridge, P. (1988) Constant rate rainfall infiltration: a versatile nonlinear model, 2. Applications of solutions. Wat. Resour. Res. 24, White, I.. & Perroux, K. M. (1989) Estimation of unsaturated hydraulic conductivity from field sorptivity measurements. Soil Sci. Soc. Am. J. 53.
P. Broadbridge. Snippets from Infiltration: where Approximate Integral Analysis is Exact.
P. Broadbridge Snippets from Infiltration: where Approximate Integral Analysis is Exact. Hydrology of 1D Unsaturated Flow in Darcy-Buckingham-Richards approach. Nonlinear diffusion-convection equations
More informationUpscaling of Richards equation for soil moisture dynamics to be utilized in mesoscale atmospheric models
Exchange Processes at the imnd Surface for a Range of Space and Time Scales (Proceedings of the Yokohama Symposium, July 1993). [AHS Publ. no. 212, 1993. 125 Upscaling of Richards equation for soil moisture
More informationψ ae is equal to the height of the capillary rise in the soil. Ranges from about 10mm for gravel to 1.5m for silt to several meters for clay.
Contents 1 Infiltration 1 1a Hydrologic soil horizons...................... 1 1b Infiltration Process......................... 2 1c Measurement............................ 2 1d Richard s Equation.........................
More informationDarcy s Law, Richards Equation, and Green-Ampt Equation
Darcy s Law, Richards Equation, and Green-Ampt Equation 1. Darcy s Law Fluid potential: in classic hydraulics, the fluid potential M is stated in terms of Bernoulli Equation (1.1) P, pressure, [F L!2 ]
More informationAdvanced Hydrology Prof. Dr. Ashu Jain Department of Civil Engineering Indian Institute of Technology, Kanpur. Lecture 6
Advanced Hydrology Prof. Dr. Ashu Jain Department of Civil Engineering Indian Institute of Technology, Kanpur Lecture 6 Good morning and welcome to the next lecture of this video course on Advanced Hydrology.
More information1D Verification Examples
1 Introduction 1D Verification Examples Software verification involves comparing the numerical solution with an analytical solution. The objective of this example is to compare the results from CTRAN/W
More information&
Inverse Problem for Transient Unsaturated Flow Identifiability and Non Uniqueness K. S. Hari Prasad & M. S. Ghidaoui Department of Civil and Structural Engineering The Hong Kong University of Science and
More informationAnalytical approach predicting water bidirectional transfers: application to micro and furrow irrigation
Advances in Fluid Mechanics VI 633 Analytical approach predicting water bidirectional transfers: application to micro and furrow irrigation D. Crevoisier Irrigation Research Unit, Cemagref Montpellier,
More informationABSTRACT INTRODUCTION
Transport of contaminants from non-aqueous phase liquid pool dissolution in subsurface formations C.V. Chrysikopoulos Department of Civil Engineering, University of California, ABSTRACT The transient contaminant
More informationTemperature dependent multiphase flow and transport
Temperature dependent multiphase flow and transport J.F. Sykes, A.G. Merry and J. Zhu Department of Civil Engineering, University of Waterloo, Waterloo, Ontario, Canada N2L 3G1 E-mail: sykesj@uwaterloo.ca
More informationWorkshop Opening Presentation of the DYNAS Project
Workshop Opening Presentation of the DYNAS Project Alexandre Ern, Eric Gaume, Cyril Kao, Jérôme Jaffré CERMICS École Nationale des Ponts et Chausées Marne-la-Vallée, France DYNAS Workshop 06-08/12/04 p.1
More informationEXAMPLE PROBLEMS. 1. Example 1 - Column Infiltration
EXAMPLE PROBLEMS The module UNSATCHEM is developed from the variably saturated solute transport model HYDRUS-1D [Šimůnek et al., 1997], and thus the water flow and solute transport parts of the model have
More informationEVALUATION OF CRITICAL FRACTURE SKIN POROSITY FOR CONTAMINANT MIGRATION IN FRACTURED FORMATIONS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More information1. Water in Soils: Infiltration and Redistribution
Contents 1 Water in Soils: Infiltration and Redistribution 1 1a Material Properties of Soil..................... 2 1b Soil Water Flow........................... 4 i Incorporating K - θ and ψ - θ Relations
More informationComparison of Averaging Methods for Interface Conductivities in One-dimensional Unsaturated Flow in Layered Soils
Comparison of Averaging Methods for Interface Conductivities in One-dimensional Unsaturated Flow in Layered Soils Ruowen Liu, Bruno Welfert and Sandra Houston School of Mathematical & Statistical Sciences,
More informationUnsaturated Flow (brief lecture)
Physical Hydrogeology Unsaturated Flow (brief lecture) Why study the unsaturated zone? Evapotranspiration Infiltration Toxic Waste Leak Irrigation UNSATURATAED ZONE Aquifer Important to: Agriculture (most
More informationResearch Scholar, PACIFIC University, Udaipur, Rajasthan, India. Department of Mathematics, UGC Visiting Professor, HNGU Patan, Gujarat, India
2018 IJSRSET Volume 4 Issue 4 Print ISSN: 2395-1990 Online ISSN : 2394-4099 Themed Section : Engineering and Technology Adomian Decomposition Study of Unidimensional Flow in Unsaturated Porous Media Khatri
More informationBuilding a Robust Numerical Model for Mass Transport Through Complex Porous Media
Presented at the COMSOL Conference 2008 Hannover Building a Robust Numerical Model for Mass Transport Through Complex Porous Media Janez Perko, Dirk Mallants Belgian Nuclear Research Centre SCK CEN Elise
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationdynamics of f luids in porous media
dynamics of f luids in porous media Jacob Bear Department of Civil Engineering Technion Israel Institute of Technology, Haifa DOVER PUBLICATIONS, INC. New York Contents Preface xvii CHAPTER 1 Introduction
More informationComputational modelling of reactive transport in hydrogeological systems
Water Resources Management III 239 Computational modelling of reactive transport in hydrogeological systems N. J. Kiani, M. K. Patel & C.-H. Lai School of Computing and Mathematical Sciences, University
More informationNumerical Solution of the Two-Dimensional Time-Dependent Transport Equation. Khaled Ismail Hamza 1 EXTENDED ABSTRACT
Second International Conference on Saltwater Intrusion and Coastal Aquifers Monitoring, Modeling, and Management. Mérida, México, March 3-April 2 Numerical Solution of the Two-Dimensional Time-Dependent
More informationA THIRD-ORDER NUMERICAL SCHEME WITH UPWIND WEIGHTING FOR SOLVING THE SOLUTE TRANSPORT EQUATION
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, VOL. 40, 1623 1637 (1997) A THIRD-ORDER NUMERICAL SCHEME WITH UPWIND WEIGHTING FOR SOLVING THE SOLUTE TRANSPORT EQUATION KANGLE HUANG, JIR[ I
More informationSimulation of Unsaturated Flow Using Richards Equation
Simulation of Unsaturated Flow Using Richards Equation Rowan Cockett Department of Earth and Ocean Science University of British Columbia rcockett@eos.ubc.ca Abstract Groundwater flow in the unsaturated
More informationA note on benchmarking of numerical models for density dependent flow in porous media
Advances in Water Resources 29 (2006) 1918 1923 www.elsevier.com/locate/advwatres A note on benchmarking of numerical models for density dependent flow in porous media B. Ataie-Ashtiani *, M.M. Aghayi
More informationAn Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation. Ahsan Munir, PhD Tom Spirka, PhD
An Introduction to COMSOL Multiphysics v4.3b & Subsurface Flow Simulation Ahsan Munir, PhD Tom Spirka, PhD Agenda Provide an overview of COMSOL 4.3b Our products, solutions and applications Subsurface
More informationMonte Carlo analysis of macro dispersion in 3D heterogeneous porous media
Monte Carlo analysis of macro dispersion in 3D heterogeneous porous media Arthur Dartois and Anthony Beaudoin Institute P, University of Poitiers, France NM2PourousMedia, Dubrovnik, Croatia 29 Sep - 3
More information1906 Dilip kumar Jaiswal et al./ Elixir Pollution 31 (2011) Available online at (Elixir International Journal)
196 Dilip kumar Jaiswal et al./ Eliir Pollution 31 (211) 196-191 ARTICLE INF O Article history: Received: 21 December 21; Received in revised form: 16 January 211; Accepted: 1 February 211; Keywords Advection,
More informationFinite Element Solution to Groundwater Transport of Solutes Undergoing Decay and Non-Linear Sorption
Hydrology Days 2010 1 Finite Element Solution to Groundwater Transport of Solutes Undergoing Decay and Non-Linear Sorption Domenico Baú Department of Civil & Environmental Engineering Colorado State University,
More informationComparison of Heat and Mass Transport at the Micro-Scale
Comparison of Heat and Mass Transport at the Micro-Scale E. Holzbecher, S. Oehlmann Georg-August Univ. Göttingen *Goldschmidtstr. 3, 37077 Göttingen, GERMANY, eholzbe@gwdg.de Abstract: Phenomena of heat
More informationThe matric flux potential: history and root water uptake
The matric flux potential: history and root water uptake Marius Heinen, Peter de Willigen (Alterra, Wageningen-UR) Jos van Dam, Klaas Metselaar (Soil Physics, Ecohydrology and Groundwater Management Group,
More informationThe generation and mobility of colloids in soils Production et mobilité des colloï des dans les sols
Scientific registration n o : 1112 Symposium n o : 4 Presentation : poster The generation and mobility of colloids in soils Production et mobilité des colloï des dans les sols WATERS Angela, CHITTLEBOROUGH
More informationOn the solution of incompressible two-phase ow by a p-version discontinuous Galerkin method
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2006; 22:741 751 Published online 13 December 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/cnm.846
More informationRADIONUCLIDE DIFFUSION IN GEOLOGICAL MEDIA
GEOPHYSICS RADIONUCLIDE DIFFUSION IN GEOLOGICAL MEDIA C. BUCUR 1, M. OLTEANU 1, M. PAVELESCU 2 1 Institute for Nuclear Research, Pitesti, Romania, crina.bucur@scn.ro 2 Academy of Scientists Bucharest,
More information3D NUMERICAL SIMULATION OF THE TRANSPORT OF EXPLOSIVES FROM UXO s AND LANDMINES
3D NUMERICAL SIMULATION OF THE TRANSPORT OF EXPLOSIVES FROM UXO s AND LANDMINES Maik Irrazábal, Julio G. Briano, Miguel Castro, and Samuel P. Hernández Center for Chemical Sensors Development, University
More informationAbstract. 1 Introduction
Correlations for mass transfer coefficients applicable to NAPL pool dissolution in subsurface formations C.V. Chrysikopoulos & T.-J. Kim Department of Civil and Environmental Engineering, University of
More informationPollution. Elixir Pollution 97 (2016)
42253 Available online at www.elixirpublishers.com (Elixir International Journal) Pollution Elixir Pollution 97 (2016) 42253-42257 Analytical Solution of Temporally Dispersion of Solute through Semi- Infinite
More informationSteady flows in unsaturated soils are stable
Steady flows in unsaturated soils are stable C. J. van Duijn and G. J. M. Pieters (c.j.v.duijn@tue.nl,g.j.m.pieters@tue.nl) Department of Mathematics and Computer Science, Eindhoven University of Technology,
More informationCHARACTERIZATION OF HETEROGENEITIES AT THE CORE-SCALE USING THE EQUIVALENT STRATIFIED POROUS MEDIUM APPROACH
SCA006-49 /6 CHARACTERIZATION OF HETEROGENEITIES AT THE CORE-SCALE USING THE EQUIVALENT STRATIFIED POROUS MEDIUM APPROACH Mostafa FOURAR LEMTA Ecole des Mines de Nancy, Parc de Saurupt, 54 04 Nancy, France
More informationAbstract. This research is supported by US Army Corps of Engineers under Grant No. DACA K-0055 with Penn State.
Modeling Water Flow and Chemical and Sediment Transport in Watershed Systems Hwai-Ping Cheng and Gour-Tsyh Yeh Department of Civil and Environmental Engineering, The Pennsylvania State University, University
More informationModeling moisture transport by periodic homogenization in unsaturated porous media
Vol. 2, 4. 377-384 (2011) Revue de Mécanique Appliquée et Théorique Modeling moisture transport by periodic homogenization in unsaturated porous media W. Mchirgui LEPTIAB, Université de La Rochelle, France,
More informationKEY WORDS heavy metals, geochemistry, solute transport, mathematical model, unsaturated flow, remediation
MODELING OF HEAVY METAL MOVEMENT IN VEGETATED, UNSATURATED SOILS WITH EMPHASIS ON GEOCHEMISTRY K.V. Nedunuri 1, R.S. Govindaraju 1, L.E. Erickson 2 and A.P. Schwab 3, 1 Department of Civil Engineering,
More informationDevelopments in a Coupled Thermal- Hydraulic-Chemical-Geomechanical Model for Soil and Concrete
Developments in a Coupled Thermal- Hydraulic-Chemical-Geomechanical Model for Soil and Concrete S.C. Seetharam * and D. Jacques 24 th October 2013 * suresh.seetharam@sckcen.be (Belgian Nuclear Research
More informationF. BEN ABDELGHANI a,b. Tunisie;
Use of The HydroGeosphere code to simulate Water flow and contaminants transport around mining waste disposed into an open pit within fractured rock mass. F. BEN ABDELGHANI a,b a Departement of Civil,
More informationAnalysis of Multiphase Flow under the Ground Water
Analysis of Multiphase Flow under the Ground Water Pramod Kumar Pant Department of Mathematics, Bhagwant University, Ajmer, Rajasthan, India Abstract The single-phase fluid flow through a porous medium
More informationChemical Hydrogeology
Physical hydrogeology: study of movement and occurrence of groundwater Chemical hydrogeology: study of chemical constituents in groundwater Chemical Hydrogeology Relevant courses General geochemistry [Donahoe]
More informationSafety assessment for disposal of hazardous waste in abandoned underground mines
Safety assessment for disposal of hazardous waste in abandoned underground mines A. Peratta & V. Popov Wessex Institute of Technology, Southampton, UK Abstract Disposal of hazardous chemical waste in abandoned
More information(9C/(9t)t = a(x,t) (92C/3x2)t + b(x,t) (9C/9 x)t + c(x,t)ct + d(x,t)
ABSTRACT Safe management Including disposal of radioactive wastes from the various parts of the nuclear fuel cycle is an important aspect of nuclear technology development. The problem of managing radioactive
More informationPrinciples of soil water and heat transfer in JULES
Principles of soil water and heat transfer in JULES Anne Verhoef 1, Pier Luigi Vidale 2, Raquel Garcia- Gonzalez 1,2, and Marie-Estelle Demory 2 1. Soil Research Centre, Reading (UK); 2. NCAS-Climate,
More informationO.R. Jimoh, M.Tech. Department of Mathematics/Statistics, Federal University of Technology, PMB 65, Minna, Niger State, Nigeria.
Comparative Analysis of a Non-Reactive Contaminant Flow Problem for Constant Initial Concentration in Two Dimensions by Homotopy-Perturbation and Variational Iteration Methods OR Jimoh, MTech Department
More informationFreezing Around a Pipe with Flowing Water
1 Introduction Freezing Around a Pipe with Flowing Water Groundwater flow can have a significant effect on ground freezing because heat flow via convection is often more effective at moving heat than conduction
More informationTowards Seamless Interactions Between Geologic Models and Hydrogeologic Applications
Towards Seamless Interactions Between Geologic Models and Hydrogeologic Applications Ross, M. 1, M. Parent 2, R. Martel 1, and R. Lefebvre 1 1 Institut National de la Recherche Scientifique (INRS-ETE),
More informationCapillary effect on watertable fluctuations in unconfined aquifers
Capillary effect on watertable fluctuations in unconfined aquifers Jun Kong 1, Cheng-Ji Shen, Pei Xin, Zhiyao Song 3, Ling Li,1,#, D.A. Barry 4, D.-S. Jeng 5, F. Stagnitti 6, D.A. Lockington and J.-Y.
More informationEXPERIENCES FROM THE SOURCE-TERM ANALYSIS OF A LOW AND INTERMEDIATE LEVEL RADWASTE DISPOSAL FACILITY
EXPERIENCES FROM THE SOURCE-TERM ANALYSIS OF A LOW AND INTERMEDIATE LEVEL RADWASTE DISPOSAL FACILITY Jin Beak Park, Joo-Wan Park, Eun-Young Lee and Chang-Lak Kim Korea Hydro & Nuclear Power Co., Ltd. (KHNP)
More informationNumerical investigations of hillslopes with variably saturated subsurface and overland flows
Numerical investigations of hillslopes with variably saturated subsurface and overland flows ARC DYNAS H. Beaugendre, T. Esclaffer, A. Ern and E. Gaume DYNAS Workshop 06-08/12/04 DYNAS Workshop 06-08/12/04
More informationSolving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions
International Workshop on MeshFree Methods 3 1 Solving Pure Torsion Problem and Modelling Radionuclide Migration Using Radial Basis Functions Leopold Vrankar (1), Goran Turk () and Franc Runovc (3) Abstract:
More informationA MATLAB Implementation of the. Minimum Relative Entropy Method. for Linear Inverse Problems
A MATLAB Implementation of the Minimum Relative Entropy Method for Linear Inverse Problems Roseanna M. Neupauer 1 and Brian Borchers 2 1 Department of Earth and Environmental Science 2 Department of Mathematics
More informationModelling Non-isothermal Flows in Porous Media: A Case Study Using an Example of the Great Artesian Basin, Australia
New Methods in Applied and Computational Mathematics (NEMACOM'98) Proceedings of the Centre for Mathematics and its Applications Modelling Non-isothermal Flows in Porous Media: A Case Study Using an Example
More information1.061 / 1.61 Transport Processes in the Environment
MIT OpenCourseWare http://ocw.mit.edu 1.061 / 1.61 Transport Processes in the Environment Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 6.
More informationA Generalized Numerical Approach for Modeling Multiphase Flow and Transport in Fractured Porous Media
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 6, No. 1, pp. 85-108 Commun. Comput. Phys. July 2009 A Generalized Numerical Approach for Modeling Multiphase Flow and Transport in Fractured Porous Media Yu-Shu
More informationRhodamine WT as a reactive tracer: laboratory study and field consequences
Tracers and Modelling in Hydrogeology (Proceedings of the TraM'2000 Conference held at Liege, Belgium, May 2000). IAHS Publ. no. 262, 2000. 201 Rhodamine WT as a reactive tracer: laboratory study and field
More informationTHREE DIMENSIONAL SIMILARITY SOLUTIONS OF THE NONLINEAR DIFFUSION EQUATION FROM OPTIMIZATION AND FIRST INTEGRALS
J. Austral. Math. Soc. {Series B) 23 (1982), 297-309 THREE DIMENSIONAL SIMILARITY SOLUTIONS OF THE NONLINEAR DIFFUSION EQUATION FROM OPTIMIZATION AND FIRST INTEGRALS J.-Y. PARLANGE, R. D. BRADDOCK, G.
More informationDEVELOPMENT OF A NUMERICAL APPROACH FOR SIMULATION OF SAND BLOWING AND CORE FORMATION
TMS (The Minerals, Metals & Materials Society), DEVELOPMENT OF A NUMERICAL APPROACH FOR SIMULATION OF SAND BLOWING AND CORE FORMATION G.F. Yao, C. W. Hirt, and
More informationRupelton layer. 1 An introduction to the problem and some basic definitions and an introduction to the problem
Abstract In this work, we use the finite volume code pdelib/thermoha2 of Enchery and Fuhrmann [2] to understand the influence of layers of small thichness on the behaviour of the fluid. To this end consider
More informationPore-scale modeling extension of constitutive relationships in the range of residual saturations
WATER RESOURCES RESEARCH, VOL. 37, NO. 1, PAGES 165 170, JANUARY 2001 Pore-scale modeling extension of constitutive relationships in the range of residual saturations Rudolf J. Held and Michael A. Celia
More informationslow velocity fast velocity
Hydrodynamic disersion Disersive transort of contaminants results from the heterogeneous distribution of water flow velocities within and between different soil ores (Figure 1, left). Disersion can be
More informationThree-dimensional Modelling of Reactive Solutes Transport in Porous Media
151 A publication of CHEMICAL ENGINEERING TRANSACTIONS VOL. 41, 214 Guest Editors: Simonetta Palmas, Michele Mascia, Annalisa Vacca Copyright 214, AIDIC Servizi S.r.l., ISBN 978-88-9568-32-7; ISSN 2283-9216
More informationJoint inversion of geophysical and hydrological data for improved subsurface characterization
Joint inversion of geophysical and hydrological data for improved subsurface characterization Michael B. Kowalsky, Jinsong Chen and Susan S. Hubbard, Lawrence Berkeley National Lab., Berkeley, California,
More informationFrozen Ground Containment Barrier
Frozen Ground Containment Barrier GEO-SLOPE International Ltd. www.geo-slope.com 1200, 700-6th Ave SW, Calgary, AB, Canada T2P 0T8 Main: +1 403 269 2002 Fax: +1 888 463 2239 Introduction Frozen soil barriers
More informationDetermination of the geomorphological instantaneous unit hydrograph using tracer experiments in a headwater basin
Hydrology, Water Resources and Ecology in Headwaters (Proceedings of the HeadWater'98 Conference held at Meran/Merano, Italy, April 1998). 1AHS Publ. no. 248, 1998. 327 Determination of the geomorphological
More informationPoiseuille Advection of Chemical Reaction Fronts
Utah State University DigitalCommons@USU All Physics Faculty Publications Physics 8- Poiseuille Advection of Chemical Reaction Fronts Boyd F. Edwards Utah State University Follow this and additional works
More informationSAFETY ASSESSMENT CODES FOR THE NEAR-SURFACE DISPOSAL OF LOW AND INTERMEDIATE-LEVEL RADIOACTIVE WASTE WITH THE COMPARTMENT MODEL: SAGE AND VR-KHNP
SAFETY ASSESSMENT CODES FOR THE NEAR-SURFACE DISPOSAL OF LOW AND INTERMEDIATE-LEVEL RADIOACTIVE WASTE WITH THE COMPARTMENT MODEL: SAGE AND VR-KHNP J. B. Park, J. W. Park, C. L. Kim, M. J. Song Korea Hydro
More informationJournal of Hydrology, 45 (1980) Elsevier Scientific Publishing Company, Amsterdam -- Printed in The Netherlands
Journal of Hydrology, 45 (1980) 289--303 289 Elsevier Scientific Publishing Company, Amsterdam -- Printed in The Netherlands [4] CONSTANT-RAINFALL INFILTRATION M. KUTILEK Laboratory of Soil Science, Department
More informationCapillary effect on water table fluctuations in unconfined aquifers
WATER RESOURCES RESEARCH, VOL. 49, 3064 3069, doi:10.1002/wrcr.20237, 2013 Capillary effect on water table fluctuations in unconfined aquifers Jun Kong, 1 Cheng-Ji Shen, 2 Pei Xin, 2 Zhiyao Song, 3 Ling
More informationTeaching Unsaturated Soil Mechanics as Part of the Undergraduate Civil Engineering Curriculum
Teaching Unsaturated Soil Mechanics as Part of the Undergraduate Civil Engineering Curriculum Delwyn G. Fredlund, Visiting Professor Kobe University, Kobe, Japan Sapporo, Hokkaido, Japan February 15, 2005
More informationNumerical Study of Heat Transfer and Contaminant Transport in an Unsaturated Porous Soil
Journal of Water Resource and Protection, 2014, 6, 1238-1247 Published Online September 2014 in SciRes. http://www.scirp.org/journal/jwarp http://dx.doi.org/10.4236/jwarp.2014.613113 Numerical Study of
More informationStability Analysis of Landslide Dam under Rainfall
Stability Analysis of Landslide Dam under Rainfall Pei-Hsun Tsai, Zheng-Yi Feng 2, Fan-Chieh Yu 3 and Jian-Han Lin 4 Associate Professor, Department of Construction Engineering, Chaoyang University of
More informationThis section develops numerically and analytically the geometric optimisation of
7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically
More information6.6 Solute Transport During Variably Saturated Flow Inverse Methods
6.6 Solute Transport During Variably Saturated Flow Inverse Methods JIÌÍ ŠIMæNEK AND MARTINUS TH. VAN GENUCHTEN, USDA-ARS, George E. Brown, Jr. Salinity Laboratory, Riverside, California DIEDERIK JACQUES,
More informationChapter 2 Water Movement and Solute Transport in Unsaturated Porous Media
Chapter 2 Water Movement and Solute Transport in Unsaturated Porous Media The unsaturated zone, also termed the vadose zone, is the portion of the subsurface above the groundwater table. It contains air
More informationDNAPL migration through interbedded clay-sand sequences
Groundwater Quality: Natural and Enhanced Restoration of Groundwater Pollution (Proceedings ofthe Groundwater Quality 2001 Conference held al Sheffield. UK. June 2001). IAHS Publ. no. 275. 2002. 455 DNAPL
More informationMHD and Thermal Dispersion-Radiation Effects on Non-Newtonian Fluid Saturated Non-Darcy Mixed Convective Flow with Melting Effect
ISSN 975-333 Mapana J Sci,, 3(22), 25-232 https://doi.org/.2725/mjs.22.4 MHD and Thermal Dispersion-Radiation Effects on Non-Newtonian Fluid Saturated Non-Darcy Mixed Convective Flow with Melting Effect
More informationOn the advective-diffusive transport in porous media in the presence of time-dependent velocities
GEOPHYSICAL RESEARCH LETTERS, VOL. 3, L355, doi:.29/24gl9646, 24 On the advective-diffusive transport in porous media in the presence of time-dependent velocities A. P. S. Selvadurai Department of Civil
More informationSolute dispersion in a variably saturated sand
WATER RESOURCES RESEARCH, VOL. 39, NO. 6, 1155, doi:10.1029/2002wr001649, 2003 Solute dispersion in a variably saturated sand Takeshi Sato Department of Civil Engineering, Gifu University, Yanagido, Gifu,
More informationSecond-Order Linear ODEs (Textbook, Chap 2)
Second-Order Linear ODEs (Textbook, Chap ) Motivation Recall from notes, pp. 58-59, the second example of a DE that we introduced there. d φ 1 1 φ = φ 0 dx λ λ Q w ' (a1) This equation represents conservation
More informationThe Use of COMSOL for Integrated Hydrological Modeling
The Use of for Integrated Hydrological Modeling Ting Fong May Chui *, David L. Freyerg Department of Civil and Environmental Engineering, Stanford University *Corresponding author: 380 Panama Mall, Terman
More informationQuantifying shallow subsurface flow and salt transport in the Canadian Prairies
Quantifying shallow subsurface flow and salt transport in the Canadian Prairies Andrew Ireson GIWS, University of Saskatchewan www.usask.ca/water Uri Nachshon Garth van der Kamp GIWS, University of Saskatchewan
More informationNumerical fractional flow modelling of inhomogeneous
Numerical fractional flow modelling of inhomogeneous air sparging E.F. Kaasschieter\ G.J. Mulder^ & J.D. van der Werff ten Bosch^ Department of Mathematics and Computing Science, hoven University of Technology,
More informationPredicting the soil-water characteristics of mine soils
Predicting the soil-water characteristics of mine soils D.A. Swanson, G. Savci & G. Danziger Savci Environmental Technologies, Golden, Colorado, USA R.N. Mohr & T. Weiskopf Phelps Dodge Mining Company,
More informationMODELING THE SALTWATER INTRUSION PHENOMENON IN COASTAL AQUIFERS - A CASE STUDY IN THE INDUSTRIAL ZONE OF HERAKLEIO IN CRETE
Global NEST Journal, Vol. 7, No 2, pp 197-203, 2005 Copyright 2005 Global NEST Printed in Greece. All rights reserved MODELING THE SALTWATER INTRUSION PHENOMENON IN COASTAL AQUIFERS - A CASE STUDY IN THE
More informationInfiltration from irrigation channels into soil with impermeable inclusions
ANZIAM J. 46 (E) pp.c1055 C1068, 2005 C1055 Infiltration from irrigation channels into soil with impermeable inclusions Maria Lobo David L. Clements Nyoman Widana (Received 15 October 2004; revised 16
More informationInternational Research Journal of Engineering and Technology (IRJET) e-issn: Volume: 04 Issue: 09 Sep p-issn:
International Research Journal of Engineering and Technology (IRJET) e-issn: 395-56 Volume: 4 Issue: 9 Sep -7 www.irjet.net p-issn: 395-7 An analytical approach for one-dimensional advection-diffusion
More informationToward the Mitigation of Spurious Cloud-Edge Supersaturation in Cloud Models
1224 M O N T H L Y W E A T H E R R E V I E W VOLUME 136 Toward the Mitigation of Spurious Cloud-Edge Supersaturation in Cloud Models WOJCIECH W. GRABOWSKI AND HUGH MORRISON National Center for Atmospheric
More informationApplication of Fuzzy Logic and Uncertainties Measurement in Environmental Information Systems
Fakultät Forst-, Geo- und Hydrowissenschaften, Fachrichtung Wasserwesen, Institut für Abfallwirtschaft und Altlasten, Professur Systemanalyse Application of Fuzzy Logic and Uncertainties Measurement in
More informationAbstract. 1 Introduction
An experiment in modelling vertical P-transport in the unsaturated zone of a former sewage farm with a simple compartment model S. Pudenz, G. Nutzmann, J. Gelbrecht Institute of Freshwater and Fish Ecology,
More informationThermal Dispersion and Convection Heat Transfer during Laminar Transient Flow in Porous Media
Thermal Dispersion and Convection Heat Transfer during Laminar Transient Flow in Porous Media M.G. Pathak and S.M. Ghiaasiaan GW Woodruff School of Mechanical Engineering Georgia Institute of Technology,
More informationDepartment of Ocean Engineering, Indian Institute of Technology-Madras, Chennai , India. *Corresponding author.
J. Earth Syst. Sci. (2018) 127:53 c Indian Academy of Sciences https://doi.org/10.1007/s12040-018-0950-3 Interaction of dissolution, sorption and biodegradation on transport of BTEX in a saturated groundwater
More informationIn all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head.
Groundwater Seepage 1 Groundwater Seepage Simplified Steady State Fluid Flow The finite element method can be used to model both steady state and transient groundwater flow, and it has been used to incorporate
More informationNumerical evaluation of a second-order water transfer term for variably saturated dual-permeability models
WATER RESOURCES RESEARCH, VOL. 40, W07409, doi:10.1029/2004wr003285, 2004 Numerical evaluation of a second-order water transfer term for variably saturated dual-permeability models J. Maximilian Köhne
More informationNotes on Spatial and Temporal Discretization (when working with HYDRUS) by Jirka Simunek
Notes on Spatial and Temporal Discretization (when working with HYDRUS) by Jirka Simunek 1. Temporal Discretization Four different time discretizations are used in HYDRUS: (1) time discretizations associated
More information