1D Verification Examples
|
|
- Pamela Gilmore
- 5 years ago
- Views:
Transcription
1 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 analyses to well-established closedform solutions for four different transport scenarios. The transport scenarios include: Case 1: no adsorption and no decay Case 2: adsorption included Case 3: decay included Case 4: adsorption and decay included. Equations for the closed-form solutions are solved using a Microsoft Excel spreadsheet. The solutions require the use of a complimentary error function (erfc). This function is available in MS Excel 2007 as part of the Analysis ToolPak Add-in, which can be installed under the Excel Options Add-Ins dialogue box. 2 Feature Highlights GeoStudio feature highlights include: 1. Comparing CTRAN/W to closed-form solutions for the advection-dispersion equation; 2. Including adsorption and decay processes; and, 3. Comparison of backward difference and central difference iteration schemes. 3 Closed-Form Solutions Closed-form analytical solutions to the advection-dispersion equation are available in the literature for one-dimensional problems involving steady-state seepage flow. The advection-dispersion equation is written as: 2 C C C S ndl nv 2 x = n + ρd nλc nλs x x t t [1] C = concentration of solute in liquid phase, t = time, D L = longitudinal hydrodynamic dispersion, ρ d = dry density, n = porosity or volumetric water content (θ), S = amount of solute sorbed per unit weight of soil, and λ = coefficient of decay. CTRAN Example File: 1D Verification Examples (pdf)(gsz) Page 1 of 12
2 The first and second terms on the left side of Equation 1 represent the dispersive and advective transport of the solute, respectively. All of the terms on the right side account for changes in mass or concentration that may occur with time due to mass flux (term 1), sorption processes (term 2), and radioactive decay of mass in solution (term 3) and mass sorbed onto the solid (term 4). 3.1 Case 1: No Adsorption and No Decay Ogata (1970) provided an analytical solution to the advection-dispersion equation for a homogeneous, isotropic, and saturated porous geological medium. The concentration at any lateral distance for a given elapsed time can be determined: C o x vt vx x + vt C = erfc exp erfc Dt D 2 Dt [2] C = concentration, C o = specified concentration in the source boundary, D = hydrodynamic dispersion coefficient, v = average linear velocity, t = elapsed time, x = distance from the source boundary, and erfc = complementary error function. The solution assumes that the concentration at a distance of x = 0 is maintained at C o for all time (i.e. C(0,t) = C o ), the initial concentration every in the flow domain is zero (i.e. C(x,0) = 0), and the flow domain is infinitely long with a concentration of zero at the far boundary (i.e. C(,t) = 0). It should be noted that the complimentary error function is related to the error function (erf) by the following Erfc(x) = 1 erf(x) and that Erf(0) = 0; erf( ) = 1; erf( x) = erf(x) 3.2 Case 2: Adsorption Included Adsorption is the physical process by which a solute adheres to a solid surface. The relationship between the mass of solute sorbed onto the solid (S) and the concentration of the solute (C) can take a variety of linear and nonlinear forms (Figure 1). The slope of the linear sorption isotherm is often referred to as the distribution coefficient (K d ). The amount of solute sorbed per dry unit weight of solid is given by S = K C [3] d CTRAN Example File: 1D Verification Examples (pdf)(gsz) Page 2 of 12
3 Linear Isotherm Non Linear Isotherm S (g/g) Concentration (g/m3) Figure 1 - Linear and nonlinear sorption isotherms If the equation for the linear sorption isotherm is substituted into the advection-dispersion equation [1] and ignoring the decay terms, the governing partial differential equation is ( KC) 2 C C C d ndl nv 2 x = n + ρ d x x t t Ignoring dispersive transport and dividing by the porosity (n) yields the following [4] C ρd C vx = 1+ Kd x n t [5] ρ + n 1 d K d is referred to as the retardation factor (R). For this specific case, the sorption process reduces the velocity by a factor of R (i.e., v x /R). Accordingly, the arrival time of the C/C o = 0.5 front for a transport problem involving advection-only is reduced by the factor R. The solution to Equation 4 for the same boundary conditions discussed above is given by (Bear, 1972) ( v ) ( v vx ) C x t x t o C erfc R + = exp erfc R Dt D 2 Dt R R [6] CTRAN Example File: 1D Verification Examples (pdf)(gsz) Page 3 of 12
4 It should be noted that the governing partial differential equation is formulated in CTRAN using the more general form shown in Equation 1. Accordingly, the formulation can accommodate unsaturated soil and non-linear adsorption functions. The retardation factor is therefore given by S R = θ + ρ d C [7] θ = volumetric water content, and S C = slope of the sorption isotherm. 3.3 Case 3: Decay Included Radioactive decay is the process in which an unstable atomic nucleus loses energy by emitting radiation. This process will reduce the concentration of radionuclides in both the dissolved and adsorbed phases. The coefficient of decay (λ) in equation [1] is given by: ln 2 λ = [8] t 1/2 t 1/2 = half-life of the radionuclide. Bear (1972, 1979) provided the following analytical solution subject to the same boundary conditions discussed above Co vx C = exp exp x erfc exp x erfc + 2 2D 2 Dt 2 Dt 2 2 x t ( v) + 4λD x+ t ( v) + 4λD ( β ) ( β ) [9] v x λ β = + 2D D 2 [10] 3.4 Case 4: Adsorption and Decay Included If radioactive decay and adsorption are included, the analytical solution becomes (Bear, 1972, 1979) ( ) ( ) 2 2 x t v 4λD x t v 4λD Co vx R R R R C = exp exp( βx) erfc + exp( βx) erfc 2 2D 2 Dt 2 Dt R R CTRAN Example File: 1D Verification Examples (pdf)(gsz) Page 4 of 12
5 v λr β = + 2D D 2 4 Boundary Conditions and Material Properties Figure 2 presents the geometry and mesh of the model domain. The model is comprised of a onedimensional column that is 3 m in length and 0.2 m in height. The mesh consists of 80 elements and 162 nodes. Although the dimensions are arbitrary, the column length was chosen to ensure that the far-field boundary condition has no effect on the results. This condition is in keeping with the analytical solutions, which assume that the C = 0 boundary is located at infinite distance. 1D Advection-Dispersion Steady-State Seepage Figure 2 Model geometry and mesh A screen capture of the KeyIn Analyses dialogue box is presented in Figure 3. A steady-state seepage analysis forms the parent analysis for each transport model. In the seepage analysis, a unit flux of m/sec was applied to the left boundary and the right boundary was assigned a hydraulic head of 1 m. The soil was assigned a hydraulic conductivity and saturated volumetric water content of m/sec and 0.5, respectively. Accordingly, the average linear velocity (v) in the flow domain is m/sec (i.e. v = q/n). Figure 3 Model structure for the 1D verification example CTRAN Example File: 1D Verification Examples (pdf)(gsz) Page 5 of 12
6 For each transport model, the left and right boundaries were set to constant concentrations of 1.0 and 0.0 g/m 3, respectively. Table 1 presents the material properties used in the analyses. The coefficient of diffusion (D) and dispersivity were set to zero and 0.1 m for all cases, yielding a hydrodynamic dispersion (αv) of m 2 /sec. A distribution coefficient K d of m 3 /g was used for the cases with adsorption, producing a retardation factor R equal to 4.0. Note that the activation concentration is set to 0 g/m 3 for each material under the KeyIn Materials dialogue box. Each model was run for an elapsed time of 6000 seconds using 60 time steps. Table 1 Material properties for transport analyses Case Disp. (m) K d (m 3 /g) Dry Density (g/m 3 ) Decay Half-life (sec) Case 1: No Adsorp./No Decay 0.1 Case 2: Adsorp. Included Case 3: Decay Included Case 4: Adsorp. & Decay Results and Discussion 5.1 Case 1: No Adsorption and No Decay Figure 4 presents results for the Case 1 analyses (no adsorption or decay) at elapsed times of 2000, 4000, and 6000 seconds. CTRAN/W was solved using the backward difference time integration scheme. The CTRAN/W results compare very well to the analytical solution; however, the analytical solution is slightly steeper than CTRAN/W. In other words, the CTRAN/W solution is slightly more spread out or smeared compared to the closed-form solution. This phenomenon is due to numerical dispersion, which is inherent in the finite element solution of the transport equation. A better match can be achieved when CTRAN/W is solved using the central difference time integration scheme (Figure 5). In general, the central difference technique provides a better solution than using backward difference for most transport problems. However, the central difference technique is susceptible to numerical oscillation, which can cause the computed concentrations to be larger or smaller than the specified maximum or minimum concentrations. Figure 6 shows a more extreme case of both numerical dispersion and oscillation. Figure 1Numerical dispersion and oscillation can only be minimized, not eliminated. Techniques for minimizing numerical dispersion and oscillation are presented in the CTRAN/W documentation. CTRAN Example File: 1D Verification Examples (pdf)(gsz) Page 6 of 12
7 Case 1: Backward Difference 2000 seconds 4000 seconds 6000 seconds CTRAN 2000 seconds CTRAN 4000 seconds CTRAN 6000 seconds Figure 4 Case 1 (no adsorption or decay) results: CTRAN/W solved using backward difference time integration scheme CTRAN Example File: 1D Verification Examples (pdf)(gsz) Page 7 of 12
8 Case 1: Central Difference 2000 seconds 4000 seconds 6000 seconds CTRAN 2000 seconds CTRAN 4000 seconds CTRAN 6000 seconds Figure 5 Case 1 results: CTRAN/W solved using central difference time integration scheme CTRAN Example File: 1D Verification Examples (pdf)(gsz) Page 8 of 12
9 seconds 0.8 Backward Difference Central Difference Figure 6 Example of numerical dispersion and oscillation 5.2 Case 2: Adsorption Included Figure 7 presents results for the Case 2 analyses, in which adsorption (R = 4.0) is included. CTRAN/W was solved using the central difference time integration scheme. Note that the scale of the x-axis ranges from 0 to 1 m. Again, CTRAN/W is in close agreement with the analytical solution. The adsorption process slows the arrival front as anticipated. For example, the position of C/C o = 0.5 front at a time of 4000 seconds is about 0.27 m, compared to 0.90 m without adsorption (i.e., approximately one-quarter; v*t/r). The small discrepancy is due to the inclusion of dispersivity in this example. CTRAN Example File: 1D Verification Examples (pdf)(gsz) Page 9 of 12
10 Case 2: Central Difference 2000 seconds 4000 seconds 6000 seconds CTRAN 2000 sec CTRAN 4000 sec CTRAN 6000 sec Figure 7 Case 2 results (adsorption included) 5.3 Case 3: Decay Included Figure 8 presents the results for Case 3 (decay included) at an elapsed time of 4000 seconds. The results from Case 1 are also included for comparison. CTRAN/W compares very well with closed-form analytical solution. The largest effect of decay occurs near the source boundary the concentration is the highest. As the concentration approaches zero, the decay component has less effect as there is little mass to decay. The affects of decay are more apparent when a slug of contaminant is introduced into the system. This was modeled by creating a region 0.1 m in length on the left side of the model domain. The region was assigned the same material used for Case 3, but with an activation concentration of 1.0 g/m 3. Accordingly, the initial mass in the system is calculated as 0.1 m 0.2 m 1 m g/m 3 = 0.01 g (i.e. length height unit width porosity Concentration). The column length was extended to 6 m and the model was run for a time of 18,000 seconds. Figure 9 presents the concentration verses distance profiles for three elapsed times. The amount of mass in the model domain (i.e. area under the curve) decreases with time due to decay. This can be checked by hand-calculation via the relationship M=M o e -λt, M o is the initial mass and λ is the coefficient of decay. For example, the calculated mass remaining in the system after an elapsed time of 6900 seconds should be about g using a λ = s -1. CTRAN/W reports the same value for the Total System Mass under View Mass Accumulation. CTRAN Example File: 1D Verification Examples (pdf)(gsz) Page 10 of 12
11 Case 3: Central Difference Case 1: 4000 Seconds seconds CTRAN 4000 sec Figure 8 Case 3 results (decay included) sec 6900 sec sec Figure 9 Effects decay on the transport of a contaminant slug CTRAN Example File: 1D Verification Examples (pdf)(gsz) Page 11 of 12
12 5.4 Case 4: Adsorption and Decay Included A comparison between CTRAN/W and the analytical solution for Case 4 is presented in Figure 10, along with the results from Case 1. There is excellent agreement between the CTRAN/W solution and the analytical solution. In fact, the two solutions are almost identical for all three elapsed times. Furthermore, a comparison between Case 1 and 4 demonstrates that CTRAN/W correctly captures the effect of both adsorption and decay. Mass is lost due to radioactive decay and the contaminant front is slowed due to adsorption Case 4: Central Difference 2000 seconds 4000 seconds 6000 seconds CTRAN 2000 sec CTRAN 4000 sec CTRAN 6000 sec Case 1: 2000 sec Case 2: 4000 sec Case 3: 6000 sec Figure 10 Case 4 results (adsorption and decay) compared to Case 1 (no adsorption/no decay) 6 Concluding Remarks In this example, results from CTRAN/W are compared to the closed-form solution of the advectiondispersion equation for four different transport scenarios. The results demonstrate that CTRAN/W is capable of modeling geochemical processes such as adsorption and decay. In general, the results from CTRAN/W provide a better-match to the analytical solution using the central difference time integration scheme. CTRAN Example File: 1D Verification Examples (pdf)(gsz) Page 12 of 12
Contaminant Modeling
Contaminant Modeling with CTRAN/W An Engineering Methodology February 2012 Edition GEO-SLOPE International Ltd. Copyright 2004-2012 by GEO-SLOPE International, Ltd. All rights reserved. No part of this
More informationGMS 8.0 Tutorial MT3DMS Advanced Transport MT3DMS dispersion, sorption, and dual domain options
v. 8.0 GMS 8.0 Tutorial MT3DMS dispersion, sorption, and dual domain options Objectives Learn about the dispersion, sorption, and dual domain options in MT3DMS Prerequisite Tutorials None Required Components
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 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 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 informationCode-to-Code Benchmarking of the PORFLOW and GoldSim Contaminant Transport Models using a Simple 1-D Domain
Code-to-Code Benchmarking of the PORFLOW and GoldSim Contaminant Transport Models using a Simple 1-D Domain - 11191 Robert A. Hiergesell and Glenn A. Taylor Savannah River National Laboratory SRNS Bldg.
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 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 informationModelling of decay chain transport in groundwater from uranium tailings ponds
Modelling of decay chain transport in groundwater from uranium tailings ponds Nair, R.N., Sunny, F., Manikandan, S.T. Student : 曹立德 Advisor : 陳瑞昇老師 Date : 2014/12/04 Outline Introduction Model Result and
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 informationNATURAL ZEOLITE AS A PERMEABLE REACTIVE BARRIER
NATURAL ZEOLITE AS A PERMEABLE REACTIVE BARRIER PREDICTION OF LEAD CONCENTRATION PROFILE THROUGH ZEOLITE BARRIER N. Vukojević Medvidović, J. Perić, M. Trgo, M. Ugrina, I. Nuić University of Split, Faculty
More informationCubzac-les-Ponts Experimental Embankments on Soft Clay
Cubzac-les-Ponts Experimental Embankments on Soft Clay 1 Introduction In the 197 s, a series of test embankments were constructed on soft clay at Cubzac-les-Ponts in France. These full-scale field tests
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 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 numerical methods for transport and reaction in porous media. Peter Frolkovič University of Heidelberg
Advanced numerical methods for transport and reaction in porous media Peter Frolkovič University of Heidelberg Content R 3 T a software package for numerical simulation of radioactive contaminant transport
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 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 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 informationPermafrost Thawing and Deformations
1 Introduction Permafrost Thawing and Deformations There is a huge interest in the effects of global warming on permafrost regions. One aspect of this is what degree of deformation may be expected if historically
More informationAnalytical solutions for water flow and solute transport in the unsaturated zone
Models for Assessing and Monitoring Groundwater Quality (Procsedines of a Boulder Symposium July 1995). IAHS Publ. no. 227, 1995. 125 Analytical solutions for water flow and solute transport in the unsaturated
More informationVISUAL SOLUTE TRANSPORT: A COMPUTER CODE FOR USE IN HYDROGEOLOGY CLASSES
VISUAL SOLUTE TRANSPORT: A COMPUTER CODE FOR USE IN HYDROGEOLOGY CLASSES Kathryn W. Thorbjarnarson Department of Geological Sciences, San Diego State University, 5500 Campanile Drive, San Diego, California
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 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 informationCharging and Transport Dynamics of a Flow-
Charging and Transport Dynamics of a Flow- Through Electrode Capacitive Deionization System Supporting information Yatian Qu, a,b Patrick G. Campbell, b Ali Hemmatifar, a Jennifer M. Knipe, b Colin K.
More informationv. 8.0 GMS 8.0 Tutorial RT3D Double Monod Model Prerequisite Tutorials None Time minutes Required Components Grid MODFLOW RT3D
v. 8.0 GMS 8.0 Tutorial Objectives Use GMS and RT3D to model the reaction between an electron donor and an electron acceptor, mediated by an actively growing microbial population that exists in both soil
More informationRT3D Double Monod Model
GMS 7.0 TUTORIALS RT3D Double Monod Model 1 Introduction This tutorial illustrates the steps involved in using GMS and RT3D to model the reaction between an electron donor and an electron acceptor, mediated
More information1. Introduction. 2. Model Description and Assumptions
Excerpt from the Proceedings of the COMSOL Conference 2010 Boston The Dissolution and Transport of Radionuclides from Used Nuclear Fuel in an Underground Repository Y. Beauregard *1, M. Gobien 2, F. Garisto
More informationSeepage through a dam embankment
Elevation GEO-SLOPE International Ltd, Calgary, Alberta, Canada www.geo-slope.com Introduction Seepage through a dam embankment The objective of this example is look at a rather simple case of flow through
More informationRT3D Double Monod Model
GMS TUTORIALS RT3D Double Monod Model This tutorial illustrates the steps involved in using GMS and RT3D to model the reaction between an electron donor and an electron acceptor, mediated by an actively
More informationRT3D Rate-Limited Sorption Reaction
GMS TUTORIALS RT3D Rate-Limited Sorption Reaction This tutorial illustrates the steps involved in using GMS and RT3D to model sorption reactions under mass-transfer limited conditions. The flow model used
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 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 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 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 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 information5. Advection and Diffusion of an Instantaneous, Point Source
1 5. Advection and Diffusion of an Instantaneous, Point Source In this chapter consider the combined transport by advection and diffusion for an instantaneous point release. We neglect source and sink
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 informationTemperature Change for Uniform Illumination
Temperature Change for Uniform Illumination Assume that the surface is uniformly illuminated by the laser Energy absorbed at the surface in a ery small depth H=I(1-R) where R = reflectiity I = light intensity
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 informationGEO-SLOPE International Ltd, Calgary, Alberta, Canada Convective air flow in heated column
1 Introduction Convective air flow in heated column The objective of this illustration is show the how coupled thermal / air flow module in AIR/W and TEMP/W can result in convective cells of air movement
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
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 informationRadionuclide Migration: Prediction Experience
Radionuclide Migration: Prediction Experience V.V. Martianov, M.Yu. Sheglov, A.V. Guskov State Unitary Enterprise MosSIA Radon 2/14, 7th Rostovsky pereulok, Moscow 119121 Russia ABSTRACT Many different
More informationTEMPORALLY DEPENDENT DISPERSION THROUGH SEMI-INFINITE HOMOGENEOUS POROUS MEDIA: AN ANALYTICAL SOLUTION
IJRRAS 6 () February www.arpapress.com/volumes/vol6issue/ijrras_6 5.pdf TEMPORALLY EPENENT ISPERSION THROUGH SEMI-INFINITE HOMOGENEOUS POROUS MEIA: AN ANALYTICAL SOLUTION R.R.Yadav, ilip Kumar Jaiswal,
More informationWM 00 Conference, February 27 March 2, 2000, Tucson, AZ DIFFUSION COEFFICIENTS OF CRITICAL RADIONUCLIDES FROM RADIOACTIVE WASTE IN GEOLOGICAL MEDIUM
DIFFUSION COEFFICIENTS OF CRITICAL RADIONUCLIDES FROM RADIOACTIVE WASTE IN GEOLOGICAL MEDIUM ABSTRACT: C. Bucur, A.Popa, C. Arsene and M.Olteanu Institute for Nuclear Research, P.O. Box 78, 0300 Pitesti
More informationSnow Melt with the Land Climate Boundary Condition
Snow Melt with the Land Climate Boundary Condition 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
More informationModelling of Advection-Dominated Transport in a Porous Column with a Time-Decaying Flow Field
Transport in Porous Media (2006) 63: 261 288 Springer 2006 DOI 10.1007/s11242-005-5392-0 Modelling of Advection-Dominated Transport in a Porous Column with a Time-Decaying Flow Field A. P. S. SELVADURAI
More informationCoupling TRIGRS and TOPMODEL in shallow landslide prediction. 1 Presenter: 王俊皓 Advisor: 李錫堤老師 Date: 2016/10/13
Coupling TRIGRS and TOPMODEL in shallow landslide prediction 1 Presenter: 王俊皓 Advisor: 李錫堤老師 Date: 016/10/13 Outline Introduction Literature review Methodology Pre-result Future work Introduction 3 Motivation
More informationJ. Environ. Res. Develop. Journal of Environmental Research And Development Vol. 8 No. 1, July-September 2013
SENSITIVITY ANALYSIS ON INPUT PARAMETERS OF 1-D GROUNDWATER FLOW GOVERNING EQUATION : SOLVED BY FINITE-DIFFERENCE AND THOMAS ALGORITHM IN MICROSOFT EXCEL Goh E.G.* 1 and Noborio K. 2 1. Department of Engineering
More informationv GMS 10.4 Tutorial RT3D Double-Monod Model Prerequisite Tutorials RT3D Instantaneous Aerobic Degradation Time minutes
v. 10.4 GMS 10.4 Tutorial RT3D Double-Monod Model Objectives Use GMS and RT3D to model the reaction between an electron donor and an electron acceptor, mediated by an actively growing microbial population
More informationIOSR Journal of Mathematics (IOSRJM) ISSN: Volume 2, Issue 1 (July-Aug 2012), PP
IOSR Journal of Mathematics (IOSRJM) ISSN: 78-578 Volume, Issue (July-Aug ), PP - Analytical Solutions of One-Dimensional Temporally Dependent Advection-Diffusion Equation along Longitudinal Semi-Infinite
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 informationA new approach to account for fracture aperture variability when modeling solute transport in fracture networks
WATER RESOURCES RESEARCH, VOL. 49, 2241 2252, doi:10.1002/wrcr.20130, 2013 A new approach to account for fracture aperture variability when modeling solute transport in fracture networks Martin Larsson,
More informationSettlement due to Pumping Well
1 Introduction Settlement due to Pumping Well A powerful and flexible feature in SIGMA/W is the ability to compute the volume change arising from an independently computed change in pore-pressure. The
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 informationEvaluation of hydrodynamic dispersion parameters in fractured rocks
Journal of Rock Mechanics and Geotechnical Engineering. 2010, 2 (3): 243 254 Evaluation of hydrodynamic dispersion parameters in fractured rocks Zhihong Zhao 1, anru Jing 1, Ivars Neretnieks 2 1 Department
More informationTHE ESTIMATION OF THE RADIONUCLIDE TRANSPORT PARAMETERS USING THE CDE MODEL
Romanian Reports in Physics, Vol. 59, No. 4, P. 1193 104, 007 Dedicated to Prof. Dumitru Barbu Ion s 70th Anniversary THE ESTIMATION OF THE RADIONUCLIDE TRANSPORT PARAMETERS USING THE CDE MODEL C. BUCUR
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 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 informationGEO-SLOPE International Ltd, Calgary, Alberta, Canada Wick Drain
1 Introduction Wick Drain This example is about modeling the behavior of a wick drain. The primary purpose here is to illustrate how interface elements can conveniently be used to include the effects of
More informationDiffusion and Ion implantation Reference: Chapter 4 Jaeger or Chapter 3 Ruska N & P Dopants determine the resistivity of material Note N lower
Diffusion and Ion implantation Reference: Chapter 4 Jaeger or Chapter 3 Ruska N & P Dopants determine the resistivity of material Note N lower resistavity than p: due to higher carrier mobility Near linear
More informationApplication of the TRT LB model to simulate pesticide transport in cultivated soils
Application of the TRT LB model to simulate pesticide transport in cultivated soils V. POT, N. ELYEZNASNI EGC INRA AgroParisTech,Thiverval-Grignon, France I. GINZBURG HBAN Cemagref, Antony, France H. HAMMOU
More informationNodalization. The student should be able to develop, with justification, a node-link diagram given a thermalhydraulic system.
Nodalization 3-1 Chapter 3 Nodalization 3.1 Introduction 3.1.1 Chapter content This chapter focusses on establishing a rationale for, and the setting up of, the geometric representation of thermalhydraulic
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 informationFilling Pond Head vs Volume Functions
Introduction Filling Pond Head vs Volume Functions The objective of this illustration is show how to model the filling of a pond where the water is seeping into the pond from the soil. The head in the
More informationChapter 2 Theory. 2.1 Continuum Mechanics of Porous Media Porous Medium Model
Chapter 2 Theory In this chapter we briefly glance at basic concepts of porous medium theory (Sect. 2.1.1) and thermal processes of multiphase media (Sect. 2.1.2). We will study the mathematical description
More informationColloid and colloid-facilitated transport modeling using HPx
Colloid and colloid-facilitated transport modeling using HPx D. Jacques 1, D. Zhou 2, J. Makselon 3, I. Engelhardt 2, and S. Thiele-Bruh 4 1 SCK-CEN, 2 TU Berlin, 3 Forschungszentrum Jülich GmbH, 4 University
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 informationChemical Kinetic Retention and Transport in Soils and Geological Media
Software Models for Chemical Kinetic Retention and Transport in Soils and Geological Media htpp:/www.spess.lsu.edu/chem_transport User s Manual Magdi Selim School of Plant, Environmental and Soil Science
More informationA Decomposition Method for Solving Coupled Multi-Species Reactive Transport Problems
Transport in Porous Media 37: 327 346, 1999. c 1999 Kluwer Academic Publishers. Printed in the Netherlands. 327 A Decomposition Method for Solving Coupled Multi-Species Reactive Transport Problems YUNWEI
More informationNUMERICAL SOLUTION OF TWO-REGION ADVECTION-DISPERSION TRANSPORT AND COMPARISON WITH ANALYTICAL SOLUTION ON EXAMPLE PROBLEMS
Proceedings of ALGORITMY 2002 Conference on Scientific Computing, pp. 130 137 NUMERICAL SOLUTION OF TWO-REGION ADVECTION-DISPERSION TRANSPORT AND COMPARISON WITH ANALYTICAL SOLUTION ON EXAMPLE PROBLEMS
More informationImpact of bentonite colloids on radionuclide transport in fractured systems results from field experiments and modelling
Ulrich Noseck, Judith Flügge, Thorsten Schäfer Impact of bentonite colloids on radionuclide transport in fractured systems results from field experiments and modelling Role of colloids in the Safety Case
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 informationEarth dam steady state seepage analysis
Engineering manual No. 32 Updated 3/2018 Earth dam steady state seepage analysis Program: FEM Water Flow File: Demo_manual_32.gmk Introduction This example illustrates an application of the GEO5 FEM module
More information6. Continuous Release - Point Source
1 6. Continuous Release - Point Source A scalar released continuously into a moving fluid and from a discrete point will form a plume that grows in the lateral dimension through diffusion and extends downstream
More informationpifreeze A Freeze / Thaw Plug-in for FEFLOW User Guide
pifreeze A Freeze / Thaw Plug-in for FEFLOW User Guide MIKE 2016 DHI headquarters Agern Allé 5 DK-2970 Hørsholm Denmark +45 4516 9200 Telephone +45 4516 9333 Support +45 4516 9292 Telefax mike@dhigroup.com
More informationAPPLICATION OF ADSORPTION PACKED-BED REACTOR MODEL FOR PHENOL REMOVAL
Tenth International Water Technology Conference, IWTC10 2006, Alexandria, Egypt 131 APPLICATION OF ADSORPTION PACKED-BED REACTOR MODEL FOR PHENOL REMOVAL M. T. Sorour, F. Abdelrasoul and W.A. Ibrahim Sanitary
More informationNUMERICAL STUDY OF NUCLIDE MIGRATION IN A NONUNIFORM HORIZONTAL FLOW FIELD OF A HIGH-LEVEL RADIOACTIVE WASTE REPOSITORY WITH MULTIPLE CANISTERS
NUMERICAL STUDY OF NUCLIDE MIGRATION IN A NONUNIFORM HORIZONTAL FLOW FIELD OF A HIGH-LEVEL RADIOACTIVE WASTE REPOSITORY WITH MULTIPLE CANISTERS RADIOACTIVE WASTE MANAGEMENT AND DISPOSAL KEYWORDS: geological
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat
More informationConsolidation of a 1D column
1 Introduction Consolidation of a 1D column This example simulates a one-dimensional consolidation test. The sides are constrained as in an odometer test. Analytical solutions are available for a case
More informationFinite difference method for solving Advection-Diffusion Problem in 1D
Finite difference method for solving Advection-Diffusion Problem in 1D Author : Osei K. Tweneboah MATH 5370: Final Project Outline 1 Advection-Diffusion Problem Stationary Advection-Diffusion Problem in
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 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 informationDistribution of pore water pressure in an earthen dam considering unsaturated-saturated seepage analysis
E3S Web of Conferences 9, 194 (16) DOI: 1.11/ e3sconf/169194 E-UNSAT 16 Distribution of pore water in an earthen dam considering unsaturated-saturated seepage analysis 1a Kumar Venkatesh, Siva Ram Karumanchi
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 information16 Rainfall on a Slope
Rainfall on a Slope 16-1 16 Rainfall on a Slope 16.1 Problem Statement In this example, the stability of a generic slope is analyzed for two successive rainfall events of increasing intensity and decreasing
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 informationTRACKING DYNAMIC HOLD-UP OF JUICE IN A CANE BED
REFEREED PAPER TRACKING DYNAMIC HOLD-UP OF JUICE IN A CANE BED LOUBSER R C Sugar Milling Research Institute NPC, c/o University of KwaZulu-Natal, Durban, South Africa, rloubser@smri.org Abstract Effective
More informationApplication of the random walk method to simulate the transport of kinetically adsorbing solutes
Groundwater Contamination (Proceedings of the Symposium held during the Third IAHS Scientific Assembly, Baltimore, MD, May 1989), IAHS Publ. no. 185, 1989 Application of the random walk method to simulate
More informationarxiv: v1 [physics.flu-dyn] 18 Mar 2018
APS Persistent incomplete mixing in reactive flows Alexandre M. Tartakovsky 1 and David Barajas-Solano 1 1 Pacific Northwest National Laboratory, Richland, WA 99352, USA arxiv:1803.06693v1 [physics.flu-dyn]
More informationPoint Propagation Model of Pollution Sources of Heavy Metals in the Surface Soil
Int. J. Electrochem. Sci., 8 (013) 5153-5160 International Journal of ELECTROCHEMICAL SCIENCE www.electrochemsci.org Point Propagation Model of Pollution Sources of Heavy Metals in the Surface Soil Hai-ming
More informationLab 14. RADIOACTIVITY
Lab 14. RADIOACTIVITY 14.1. Guiding Question What are the properties of different types of nuclear radiation? How does nucelar decay proceed over time? 14.2. Equipment 1. ST360 Radiation Counter, G-M probe
More informationModified Cam-clay triaxial test simulations
1 Introduction Modified Cam-clay triaxial test simulations This example simulates a series of triaxial tests which can be used to verify that Modified Cam-Clay constitutive model is functioning properly.
More informationFlow and Transport. c(s, t)s ds,
Flow and Transport 1. The Transport Equation We shall describe the transport of a dissolved chemical by water that is traveling with uniform velocity ν through a long thin tube G with uniform cross section
More informationEXTRAPOLATION STUDIES ON ADSORPTION OF THORIUM AND URANIUM AT DIFFERENT SOLUTION COMPOSITIONS ON SOIL SEDIMENTS Syed Hakimi Sakuma
EXTRAPOLATION STUDIES ON ADSORPTION OF THORIUM AND URANIUM AT DIFFERENT SOLUTION COMPOSITIONS ON SOIL SEDIMENTS Syed Hakimi Sakuma Malaysian Institute for Nuclear Technology Research (MINT), Bangi, 43000
More informationTreatment of Colloids in the Safety Case
Treatment of Colloids in the Safety Case Clay Colloids in Aqueous Systems 3 4 February 2016, Berlin Dr Amy Shelton, Radioactive Waste Management (RWM) KBS-3 Concept Based on the multi- barrier principles
More informationNumerical Simulations of Radionuclide Transport through Clay and Confining Units in a Geological Repository using COMSOL
Numerical Simulations of Radionuclide Transport through Clay and Confining Units in a Geological Repository using COMSOL J. Hansmann *1, M. L. Sentis 1, C. Belardinelli 2, B. J. Graupner 1, M. Hugi 1 and
More informationThe HYDRUS-2D Software Package for Simulating the. Two-Dimensional Movement of Water, Heat, and. Multiple Solutes in Variably-Saturated Media
The HYDRUS-2D Software Package for Simulating the Two-Dimensional Movement of Water, Heat, and Multiple Solutes in Variably-Saturated Media Version 2.0 April 1999 U. S. SALINITY LABORATORY AGRICULTURAL
More informationAppendix A. Longwave radiation in a canopy stand
213 Appendix A. ongwave radiation in a canopy stand The longwave radiation transfer equations can be written as (Ross, 1981): di ld (,r) d cosθ = G(,r)i ld (,r)+η,r ( ) (A.1) di lu ( ),r d cosφ = G,r (
More information