Simulating flow and transport with advanced geochemical reactions Recent developments using PHREEQC as reaction engine
|
|
- Buddy Franklin
- 6 years ago
- Views:
Transcription
1 Simulating flow and transport with advanced geochemical reactions Recent developments using PHREEQC as reaction engine Dr. Laurin Wissmeier AF-Consult Switzerland AG, Grundwasserschutz und Entsorgung Täfernstrasse 26 CH-5405 Baden Schweiz ABSTRACT: Reactive transport phenomena are relevant to fields such as groundwater protection, geothermal energy production, heap leaching, mine rehabilitation and nuclear waste disposal. Thus, reactive transport modelling is key in predicting the evolution of the water quality and solute behavior for safety assessments and process optimization. In recent years, couplings of groundwater flow and transport simulators with the geochemical modelling framework PHREEQC have become increasingly popular to describe adequately reaction processes that occur simultaneous to flow and transport. We present the new geochemical reaction engine PhreeqcRM [1], a wrapper class for PHREEQC that has been specifically designed for software couplings with existing flow and transport simulators. We explain the key functionality of PhreeqcRM and demonstrate its potential and performance as a reaction engine for FEFLOW. This implementation is then verified using an analytical solution for a 3D reactive transport example as well as the well-established MoMaS benchmarks for nuclear waste disposal. INTRODUCTION Since the release of IPhreeqc [2], several authors and groups have used this general purpose application programming interface (API) for the geochemical modelling framework PHREEQC [3] to couple it to transport codes in a wide variety of contexts, which indicates the scientific and industrial interest in using PHREEQC as a reaction module. A coupling with COMSOL was developed by Wissmeier and Barry [4], which gave access to the full range of COMSOL s flow and transport capabilities together with the complete set of geochemical reactions in PHREEQC. In the meantime, the outlined general strategy for couplings with different flow and transport simulators in [4] has been followed by Nardi et al. [5], who published another COMSOL-IPhreeqc coupling. Recently, also Korrani et al. [6] followed this strategy in their coupling to UTCHEM. The free environmental flow- and transport-modelling platform OpenGeoSys [7] maintains an interface to IPhreeqc for coupled thermo-hydro-mechanical-chemical (THMC) simulations. Specialized tools using IPhreeqc have been developed by Takahashi and Ishida [8] for cementitious materials and Huber et al. [9] for the paper-making process. Patel et al. [10] implemented IPhreeqc as reaction engine for pore-scale multicomponent reactive transport that used a Lattice-Boltzmann approach. However, IPhreeqc is not specifically tailored to the needs of couplings with transport simulators. Instead, it is designed to maintain the full capabilities of PHREEQC, independent of the context of use. Therefore, the API requires substantial coding in the client software to compose the adequate PHREEQC commands that handle the data exchange and run reactions. In this paper, we present PhreeqcRM, a new reaction module based on IPhreeqc that is specifically designed for couplings with environmental flow and transport simulators. PhreeqcRM provides a high-level interface that allows multicomponent transport codes to implement geochemical reactions with a minimum amount of programming, while maintaining the full functionality of PHREEQC s reaction capabilities. In the first part of the paper, we describe the coupling strategy together with key methods of its API. In the second part, we present the implementation of PhreeqcRM as a reaction engine for the groundwater modelling system FEFLOW. The new coupling is verified through two test cases: a three-dimensional (3D) analytical solution by Sun et al. [11] and the MoMaS (Modeling, Mathematics and numerical Simulations related to nuclear waste management problems) reactive transport benchmark of Groupement de Recherche (GdR) [12].
2 PhreeqcRM is available at Included in the distributions are source files, compilation files (CMake and configure), documentation of all methods of the API, and simple advection examples in C, Fortran, and C++. GENERAL STRATEGY FOR THE IMPLEMENTATION OF PHREEQCRM IN A MULTICOMPONENT TRANSPORT SIMULATOR Figure 1 presents a general strategy for the use of PhreeqcRM as a reaction engine for a transport simulator using common flow chart symbols. The initialization process and the reaction process within the main simulation process are detailed on the right. The flow chart contains the key API tasks that are used in these processes. Initialization process Start simulation Start initialization Load Settings Initialization End of simulation period? False True Create and initialize PhreeqcRM Set initial conditions Set boundary conditions Flow and solute transport Is reaction time step? True Reactions False End initialization Reactions process Start reactions Transfer data to PhreeqcRM Visualization Run PhreeqcRM reactions for t Finalize Transfer data to transport End simulation End reactions Figure 1: Generic program flow for a reactive transport simulator using operator splitting (adapted from [1]).
3 PhreeqcRM API PhreeqcRM is designed as a C++ class with methods for instantiation and destruction, initialization of initial conditions and reactions, determination of boundary conditions and running reaction calculations for a series of model cells and time steps. In addition to C++ methods, wrappers for C and Fortran are provided. Parallelization PhreeqcRM is designed for parallel reaction calculations, using either multi-threading on shared memory systems with OpenMP or multi-processing on distributed memory systems with MPI. The kind of parallelization has to be selected at compile time using the preprocessor definitions USE_OPENMP or USE_MPI. The number of parallel threads for OpenMP has to be provided as a parameter to the constructor at runtime whereas the number of parallel processes for MPI has to be provided as an argument to mpiexec, the command that launches an MPI job. PhreeqcRM has been successfully compiled with OpenMP implementations for Windows (Visual Studio 2010/2012/2013) and Scientific Linux and with MPI implementations for Windows (MSMPI from Microsoft ) and Linux (OPENMPI, version 1.5.4). The efficiency of parallelization is demonstrated by the speedup that increasing numbers of threads (Multithreaded) or processes (Multiprocessing) produce for the easy 1D MoMaS reactive transport benchmark. IMPLEMENTATION OF PHREEQCRM AS THE REACTION ENGINE FOR FEFLOW The implementation of PhreeqcRM as the reaction engine for reactive transport calculations in FEFLOW uses the sequential non-iterative approach for operator splitting, which is ideally suited for parallelization. PhreeqcRM is implemented as a plug-in that calculates reactions at predefined time steps or at every n th transport step, in which case the length of the time step is determined by FEFLOW s automatic time-stepping algorithm for flow and transport. FEFLOW s Interface Manager (IFM) provides the API for the coupling as well as plug-in methods that execute at certain events during program flow. The event-based methods are implemented as callbacks for the FEFLOW main program. For the problem definition, result visualization, and result storage, FEFLOW s existing infrastructure and GUI elements are used as much as possible. New GUI elements for specific coupling settings and for the association of PHREEQC input files with FEFLOW boundary and initial conditions are designed using the platform-independent Qt programing framework. Every node in FEFLOW s finite-element mesh corresponds to a reaction cell in PhreeqcRM. Component concentrations, solution temperature and pressure, and liquid phase saturation are transferred from FEFLOW s transport nodes to the reaction cells in PhreeqcRM before each reaction step. Saturation is transferred from PhreeqcRM to FEFLOW after each reaction step, which allows FEFLOW to simulate saturation-modifying chemical reactions. PLUG-IN OPERATION The starting point of a coupled simulation is a fully functional FEFLOW model for flow and transport with at least one (placeholder) mass transport species. Transport properties of this species (e.g., diffusion coefficient, porosity) are used for all geochemical components. Geochemical components are determined from the associated PHREEQC input files and added automatically at the start of the simulation. For the manual definition of reaction steps a power curve is required that defines reaction steps and changes in boundary conditions. Furthermore, a nodal user data distribution is needed that identifies nodes with the same geochemical initial conditions. Initial liquid phase saturation and temperatures are taken from the FEFLOW model. The FEFLOW plug-in uses separate PHREEQC input files to define the geochemical boundary and initial conditions. A file selection GUI with immediate error checking is used to associate PHREEQC files with FEFLOW nodes that share the same boundary and initial conditions (that is, nodes with the same value for constant concentration boundary or the same nodal user data). For initial conditions, the solution with the highest user number in the file is transferred as the initial solution for all of the specified transport cells; any reactants with the same user number also are transferred. For boundary conditions, the solution with the highest user number in each PHREEQC file is used to define the
4 solution composition at the associated nodes. All coupling-relevant information can be saved in the FEFLOW fem-file through the serialization functionality in the IFM. Output from SELECTED_OUTPUT and USER_PUNCH keywords in the PHREEQC files that are used to define boundary and initial conditions is transferred to FEFLOW as an additional nodal user data definition, which is updated after each reaction step. Nodal user data definitions can be saved together with FEFLOW s result files (dac- and dar-files). This allows for the visualization of geochemical parameters during and after a coupled simulation using FEFLOW s built-in post-processing tools. CODE VERIFICATION A 3D reactive transport analytical solution [11] and a series of MoMaS reactive transport benchmarks [12] were used to test the simulator. Kinetic Decay-Chain Test Case For the verification of the calculation of simple kinetic reactions in PhreeqcRM and its implementation as the reaction engine for FEFLOW the analytical solution of Wexler [13] in combination with the methodology for multi-species transport of Sun et al. [11] is used. The example is adapted from Example 2 of the PHAST manual [14] and simulates a decay chain of four artificial species (C 1 C 4 ) according to the first-order rate expressions dc i (t) dt k i C i (t), i = 1 = { k i 1 C i 1 k i C i (t), 2 i 4 (1) In Eq. (1) C i are species concentrations in mol/l and k 1 = /d, k 2 = /d, k 3 = /d and k 4 = /d are rate constants. The concentration C 1 is introduced over a specified area at the center of the y-z plane with flow in x-direction at 0.2 m/d, and the species form and decay in accordance with the specified reaction rates.
5 Figure 2: Species concentrations (mol/l) for a first-order decay chain in a steady flow field as calculated by the FEFLOW plug-in compared to the analytical solution of Wexler [13] and Sun et al. [11] (adapted from [1]). In addition to the excellent agreement with the analytical solution, as shown in Figure 2, results from FEFLOW are nearly identical to an implementation of PhreeqcRM as reaction engine for PHAST [1]. Reactive Transport Benchmark of MoMaS The MoMaS reactive transport benchmarks [12], commonly known as MoMaS, define test cases for steady-state flow with advection- or dispersion-dominated transient solute transport in 1D and 2D domains together with three versions (easy, medium, and hard) of an artificial reaction network. The general applicability of the SNIA to the MoMaS exercise has been demonstrated for the SPECY code in Carrayrou [15] and Carrayrou et al. [16]. Details of the MoMaS definitions are provided by Carrayrou et al. [12] and Bourgeat et al. [17]. For the verification of the FEFLOW plug-in with the 2D MoMaS cases, an irregular finite-element mesh with 3753 elements and 1970 nodes was used.
6 Figure 3: FEFLOW concentration results for the MoMaS 2D advective cases: (a) X3, easy case at time 1000; (b) X2, medium case at time 1000; (c) CP1, hard case at time 2000; and (d) same as (c) calculated with fine mesh (adapted from [1]). Figure 4: FEFLOW concentration results for the MoMaS 2D dispersive cases; (a) and (b) S, easy case at time 10 with coarse and fine mesh, respectively; (c) X2, medium case at time 10; (d) CP1, hard case at time 1000 (adapted from [1]). Results in Figures 3 and 4 compare well with previously published results in Carrayrou et al. [16], Mayer et al. [18] and Hoffmann et al. [19]. REFERENCES
7 [1] Parkhurst DL, Wissmeier L. PhreeqcRM: A reaction module for transport simulators based on the geochemical model PHREEQC. Advances in Water Resources 2015;83:176\uc0\u8211{}89. doi: /j.advwatres [2] Charlton SR and Parkhurst DL. Modules based on the geochemical model PHREEQC for use in scripting and programming languages. Computers & Geosciences 2011;37(10): [3] Parkhurst DL and Appelo CAJ. Description of Input and Examples for PHREEQC Version 3 A Computer Program for Speciation, Batch-Reaction, One-Dimensional Transport, and Inverse Geochemical Calculations. U.S. Department of the Interior, U.S. Geological Survey Techniques and Methods 6 A43, [4] Wissmeier L and Barry DA. Simulation tool for variably saturated flow with comprehensive geochemical reactions in two- and three-dimensional domains. Environmental Modelling & Software 2011;26(2011): [5] Nardi A, Idiart A, Trinchero P, de Vries LM, and Molinero J. Interface COMSOL-PHREEQC (icp), an efficient numerical framework for the solution of coupled multiphysics and geochemistry. Computers & Geosciences 2014;69: [6] Korrani AKN, Sepehrnoori K and Delshad M. Coupling IPhreeqc with UTCHEM to model reactive flow and transport, Computers & Geosciences : [7] Kolditz O, Bauer S, Bilke L, Böttcher N, Delfs JO, Fischer T, Görke UJ, Kalbacher T, Kosakowski G, McDermott CI, Park CH, Radu F, Rink K, Shao H, Shao HB, Sun F, Sun YY, Singh AK, Taron J, Walther M, Wang W, Watanabe N, Wu Y, Xie M, Xu W, and Zehner B. OpenGeoSys: an open-source initiative for numerical simulation of thermo-hydromechanical/chemical (THM/C) processes in porous media. Environmental Earth Sciences 2012;67(2): [8] Takahashi Y and Ishida T. Modeling of coupled mass transport and chemical equilibrium in cement-solidified soil contaminated with heavy-metal ions. Construction and Building Materials 2014;67, Part A(0): [9] Huber P, Nivelon S, Ottenio P, and Nortier P. Coupling a Chemical Reaction Engine with a Mass Flow Balance Process Simulation for Scaling Management in Papermaking Process Waters. Industrial & Engineering Chemistry Research 2012;52(1): [10] Patel RA, Perko J, Jacques D, De Schutter G, Van Breugel K, and Ye G. A versatile porescale multicomponent reactive transport approach based on lattice Boltzmann method: Application to portlandite dissolution. Physics and Chemistry of the Earth 2014;70 71(0): [11] Sun Y, Petersen JN, and Clement TP. Analytical solutions for multiple species reactive transport in multiple dimensions. Journal of Contaminant Hydrology 1999;35(4): [12] Carrayrou J, Kern M, and Knabner P. Reactive transport benchmark of MoMaS. Computational Geosciences 2010;14(3): [13] Wexler EJ. Analytical solutions for one-, two-, and three-dimensional solute transport in ground-water systems with uniform flow. U.S. Department of the Interior, U.S. Geological Survey, Techniques of Water-Resources Investigations 3 B7, [14] Parkhurst DL, Kipp KL, and Charlton SR. PHAST Version 2--A Program for Simulating Groundwater Flow, Solute Transport, and Multicomponent Geochemical Reactions. U.S. Geological Survey Techniques and Methods 6-A [15] Carrayrou J. Looking for some reference solutions for the reactive transport benchmark of MoMaS with SPECY. Computational Geosciences 2010;14(3): [16] Carrayrou J, Hoffmann J, Knabner P, Kräutle S, de Dieuleveult C, Erhel J, Van der Lee J, Lagneau V, Mayer KU, and MacQuarrie KTB. Comparison of numerical methods for simulating strongly nonlinear and heterogeneous reactive transport problems the MoMaS benchmark case. Computational Geosciences 2010;14(3): [17] Bourgeat A, Bryant S, Carrayrou J, Dimier A, Duijn CJV, Kern M, Knabner P, and Leterrier N. GDR MoMaS Benchmark Reactive Transport. Centre National de la Recherche Scientifique [18] Mayer KU, Frind EO, and Blowes DW. Multicomponent reactive transport modeling in variably saturated porous media using a generalized formulation for kinetically controlled reactions. Water Resources Research 2002;38(9): [19] Hoffmann J. Results of the GdR MoMaS Reactive Transport Benchmark with RICHY2D. Institut für Angewandte Mathematik, 2008.
Numerical modelling of coupled TH processes to facilitate analysis of geothermal reservoirs - Groß Schönebeck as an example-
Numerical modelling of coupled TH processes to facilitate analysis of geothermal reservoirs - Groß Schönebeck as an example- N Watanabe 1, G Blöcher 2, M Cacace 2, H Shao 1, O Kolditz 1 1 Helmholtz Centre
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 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 informationCOUPLING METHODOLOGY WITHIN THE SOFTWARE PLATFORM ALLIANCES
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2005 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) CIMNE, Barcelona, 2005 COUPLING METHODOLOGY WITHIN
More informationA parallelization scheme to simulate reactive transport in the subsurface environment with OGS#IPhreeqc
Geosci. Model Dev., 8, 3333 3348, 2015 doi:10.5194/gmd-8-3333-2015 Author(s) 2015. CC Attribution 3.0 License. A parallelization scheme to simulate reactive transport in the subsurface environment with
More informationDevelopment of Process Coupling System for the Numerical Experiment of High Level Radioactive Waste
Development of Process Coupling System for the Numerical Experiment of High Level Radioactive Waste ATSUSHI NEYAMA Environmental Engineering Group Computer Software Development Co., Ltd. 15-1, Tomihisa-Cho,
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 informationCOUPLED PHYSICAL AND GEOCHEMICAL MODEL FOR EXTERNAL SULFATE ATTACK IN CONCRETE
- Technical Paper - COUPLED PHYSICAL AND GEOCHEMICAL MODEL FOR EXTERNAL SULFATE ATTACK IN CONCRETE Yogarajah ELAKNESWARAN *1, Tetsuya ISHIDA *2 ABSTRACT In this study, a coupled physical and geochemical
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 Suresh C Seetharam 1* and Diederik Jacques 1 1 SCK CEN (Belgian Nuclear Research Centre), Belgium *Corresponding
More informationReactive transport modeling of carbonate diagenesis on unstructured grids
Reactive transport modeling of carbonate diagenesis on unstructured grids Alina Yapparova Chair of Reservoir Engineering Montanuniversität Leoben Dolomitization and geochemical modeling Hollywood may never
More informationLATTICE BOLTZMANN BASED MULTICOMPONENT REACTIVE TRANSPORT MODEL COUPLED WITH GEOCHEMICAL SOLVER FOR PORE SCALE SIMULATIONS
V International Conference on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2013 S. Idelsohn, M. Papadrakakis and B. Schrefler (Eds) LATTICE BOLTZMANN BASED MULTICOMPONENT
More informationLawrence Berkeley National Laboratory Lawrence Berkeley National Laboratory
Lawrence Berkeley National Laboratory Lawrence Berkeley National Laboratory Title TOUGHREACT: a new code of the TOUGH Family for Non-Isothermal multiphase reactive geochemical transport in variably saturated
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 informationKatharina Benisch, Bastian Graupner, Dedong Li, Christof Beyer, Addisalem Mitiku, Sebastian Bauer
Katharina Benisch, Bastian Graupner, Dedong Li, Christof Beyer, Addisalem Mitiku, Sebastian Bauer Institute of Geosciences, University of Kiel, Contact: kb@gpi.uni-kiel.de 6 th Trondheim Carbon Capture
More informationA Newton Krylov method for coupling transport with chemistry in porous media
A Newton Krylov method for coupling transport with chemistry in porous media Laila Amir 1,2 Michel Kern 2 1 Itasca Consultants, Lyon 2 INRIA Rocquencourt Joint work with J. Erhel (INRIA Rennes /IRISA)
More informationCoupling transport with geochemistry Applications to CO2 sequestration
Coupling transport with geochemistry Applications to CO2 sequestration Michel Kern with L. Amir, B. Gueslin, A. Taakili INRIA Paris Rocquencourt 3ème Conférence Internationale de la Société Marocaine de
More informationA COUPLED REACTIVE TRANSPORT MODEL FOR CONTAMINANT LEACHING FROM CEMENTITIOUS WASTE MATRICES ACCOUNTING FOR SOLID PHASE ALTERATIONS
A COUPLED REACTIVE TRANSPORT MODEL FOR CONTAMINANT LEACHING FROM CEMENTITIOUS WASTE MATRICES ACCOUNTING FOR SOLID PHASE ALTERATIONS D. JACQUES *1, J. ŠIMŮNEK **, D. MALLANTS *, M.TH. VAN GENUCHTEN ***,
More informationSearch and Discovery Article #50999 (2014)** Posted August 18, Abstract
Oil Degradation in the Gullfaks Field (Norway): How Hydrogeochemical Modeling can Help to Decipher Organic- Inorganic Interactions Controlling CO 2 Fate and Behavior* Wolfgang van Berk 1, Yunjiao Fu 2,
More informationThe Geochemist s Workbench
The Geochemist s Workbench GWB Essentials The Geochemist's Workbench 9.0 Essentials edition makes quick work of the tasks most commonly encountered in geochemical modeling. With a few mouse clicks, you
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 informationReactive transport in porous media: formulations, non-linear solvers and preconditioners
Reactive transport in porous media: formulations, non-linear solvers and preconditioners Michel Kern with L. Amir, B. Gueslin, A. Taakili Institut National de Recherche en Informatique et Automatique High
More informationPreconditioning the Newton-Krylov method for reactive transport
Preconditioning the Newton-Krylov method for reactive transport Michel Kern, L. Amir, A. Taakili INRIA Paris-Rocquencourt Maison de la Simulation MoMaS workshop Reactive Transport Modeling in the Geological
More informationLabor für Endlagersicherheit
WIR SCHAFFEN WISSEN HEUTE FÜR MORGEN Sergey V. Churakov :: Laboratory for Waste Management :: Paul Scherrer Institut Labor für Endlagersicherheit NES Kompetenzen und Highlights, Oktober 18, 2016, PSI Outline
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 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 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 informationThermo-Hydro-Mechanical-Chemical (THMC) Modelling of the Bentonite Barrier in Final Disposal of High Level Nuclear Waste
Presented at the COMSOL Conference 2008 Hannover Thermo-Hydro-Mechanical-Chemical (THMC) Modelling of the Bentonite Barrier in Final Disposal of High Level Nuclear Waste, Markus Olin, Veli-Matti Pulkkanen,
More informationSolutions of the Elder problem for a wide range of the Rayleigh number
European Water 55: 31-4, 216. 216 E.W. Publications Solutions of the Elder problem for a wide range of the Rayleigh number J.N.E. Papaspyros *, K.N. Moutsopoulos ** and E. Tzavellas Laboratory of Ecological
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 informationSubsurface Environmental Simulation Benchmarking Workshop V
Subsurface Environmental Simulation Benchmarking Workshop V A Coruña (Spain) October 13 th -15 th, 2016 www.turismocoruna.com J. Samper, L. Montenegro & A. Naves Eds. Background The complexity and detail
More informationThermohaline Flow and Reactive Solute Transport in Fractured Rock
Thermohaline Flow and Reactive Solute Transport in Fractured Rock Thomas GRAF Institute of Fluid Mechanics in Civil Engineering Alberta Geological Survey, Edmonton 2008-2009 Basin Hydrogeologist Sasdf
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 information4 Carbonates and CO2
PhreeqcI Introductory Course (Exercises booklet, chapters 4 and 5) by Manuel Prieto Department of Geology, University of Oviedo, Spain mprieto@ @geol.uniovi.es http://wwwbrr.cr.usgs.gov/projects/gwc_coupled/phreeqci/
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 informationUranium in an ISR Wellfield using Kinetic-based
Geochemical Fate and Transport of Dissolved Uranium in an ISR Wellfield using Kinetic-based Modeling Approach and PHREEQC and PHAST Aaron Payne, P.G., M.S. Senior Hydrogeologist PETROTEK ENGINEERING CORPORATION
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 informationHYDROGEOCHEMICAL PROCESSES ASSOCIATED WITH SEAWATER INTRUSION
HYDROGEOCHEMICAL PROCESSES ASSOCIATED WITH SEAWATER INTRUSION 1 1 Departamento de Ingeniería Química PART 1 I. Corrections for solubility calculations II. Species distribution and mineral saturation (PHREEQCI)
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 informationReactive-transport modelling of electrokinetic extraction of heavy metals from marine sediments
EREM63 Reactive-transport modelling of electrokinetic extraction of heavy metals from marine sediments Matteo Masi a, *, Alessio Ceccarini b, Renato Iannelli a a University of Pisa, Department of Energy,
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 information2. Governing Equations. 1. Introduction
Multiphysics Between Deep Geothermal Water Cycle, Surface Heat Exchanger Cycle and Geothermal Power Plant Cycle Li Wah Wong *,1, Guido Blöcher 1, Oliver Kastner 1, Günter Zimmermann 1 1 International Centre
More informationThermal Analysis Contents - 1
Thermal Analysis Contents - 1 TABLE OF CONTENTS 1 THERMAL ANALYSIS 1.1 Introduction... 1-1 1.2 Mathematical Model Description... 1-3 1.2.1 Conventions and Definitions... 1-3 1.2.2 Conduction... 1-4 1.2.2.1
More informationAthena Visual Software, Inc. 1
Athena Visual Studio Visual Kinetics Tutorial VisualKinetics is an integrated tool within the Athena Visual Studio software environment, which allows scientists and engineers to simulate the dynamic behavior
More informationTHERMAL HYDRAULIC REACTOR CORE CALCULATIONS BASED ON COUPLING THE CFD CODE ANSYS CFX WITH THE 3D NEUTRON KINETIC CORE MODEL DYN3D
THERMAL HYDRAULIC REACTOR CORE CALCULATIONS BASED ON COUPLING THE CFD CODE ANSYS CFX WITH THE 3D NEUTRON KINETIC CORE MODEL DYN3D A. Grahn, S. Kliem, U. Rohde Forschungszentrum Dresden-Rossendorf, Institute
More informationMODIFICATION IN ADINA SOFTWARE FOR ZAHORSKI MATERIAL
MODIFICATION IN ADINA SOFTWARE FOR ZAHORSKI MATERIAL Major Maciej Major Izabela Czestochowa University of Technology, Faculty of Civil Engineering, Address ul. Akademicka 3, Częstochowa, Poland e-mail:
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 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 informationDEVELOPMENT OF A CAPE-OPEN 1.0 SOCKET. Introduction
DEVELOPMENT OF A CAPE-OPEN 1. SOCKET Eric Radermecker, Belsim S.A., Belgium Dr. Ulrika Wising, Belsim S.A., Belgium Dr. Marie-Noëlle Dumont, LASSC, Belgium Introduction VALI is an advanced data validation
More informationCALCULATING CATION EXCHANGE WITH PHREEQC (Version 2)
CALCULATING CATION EXCHANGE WITH PHREEQC (Version 2) C.A.J. Appelo Hydrochemical Consultant Valeriusstraat 11 1071 MB Amsterdam, The Netherlands appt@xs4all.nl D.L. Parkhurst U.S. Geological Survey Denver
More informationComparison of numerical schemes for multiphase reactive transport
Énergies renouvelables Production éco-responsable Transports innovants Procédés éco-efficients Ressources durables Comparison of numerical schemes for multiphase reactive transport T. Faney, A. Michel,
More informationComparison in Some Modelling Tasks of Spent Nuclear Fuel Disposal
Presented at the COMSOL Conference 010 Paris COMSOL Multiphysics TOUGHREACT Numerrin Comparison in Some Modelling Tasks of Spent Nuclear Fuel Disposal COMSOL Conference 010, November 17-19, Paris Aku Itälä
More informationREDUCING COST AND ENVIRONMENTAL IMPACT OF GEOTHERMAL POWER THROUGH MODELING OF CHEMICAL PROCESSES IN THE RESERVOIR
PROCEEDINGS, Twenty-Sixth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 29-31, 21 SGP-TR-168 REDUCING COST AND ENVIRONMENTAL IMPACT OF GEOTHERMAL POWER
More informationA Hybrid Deterministic / Stochastic Calculation Model for Transient Analysis
A Hybrid Deterministic / Stochastic Calculation Model for Transient Analysis A. Aures 1,2, A. Pautz 2, K. Velkov 1, W. Zwermann 1 1 Gesellschaft für Anlagen- und Reaktorsicherheit (GRS) ggmbh Boltzmannstraße
More informationNUMERICAL MODELING OF TRANSIENT ACOUSTIC FIELD USING FINITE ELEMENT METHOD
POZNAN UNIVE RSITY OF TE CHNOLOGY ACADE MIC JOURNALS No 73 Electrical Engineering 213 Lukáš KOUDELA* Jindřich JANSA* Pavel KARBAN* NUMERICAL MODELING OF TRANSIENT ACOUSTIC FIELD USING FINITE ELEMENT METHOD
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 informationGEOCHEMISTRY, GROUNDWATER AND POLLUTION,
GEOCHEMISTRY, GROUNDWATER AND POLLUTION, 2 ND EDITION C.A.J. APPELO Hydrochemical Consultant, Amsterdam, the Netherlands D. POSTMA Environment & Resources DTU, Technical University of Denmark, Kgs. Lyngby,
More informationDense Arithmetic over Finite Fields with CUMODP
Dense Arithmetic over Finite Fields with CUMODP Sardar Anisul Haque 1 Xin Li 2 Farnam Mansouri 1 Marc Moreno Maza 1 Wei Pan 3 Ning Xie 1 1 University of Western Ontario, Canada 2 Universidad Carlos III,
More informationINTRODUCTION TO GEOCHEMICAL AND REACTIVE TRANSPORT MODELING. Ondra Sracek
INTRODUCTION TO GEOCHEMICAL AND REACTIVE TRANSPORT MODELING Ondra Sracek Principal types of geochemical modeling: 1. Speciation (for exemple, for lead); 2. Inverse geochemical modeling; 3. Direct geochemical
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 Engineering (CH_ENG)
Chemical Engineering (CH_ENG) 1 Chemical Engineering (CH_ENG) CH_ENG 1000: Introduction to Chemical Engineering Orientation course for freshmen-level students. Introduction to careers and opportunities
More informationA Rational Approach to Pit Lake Chemistry Modeling
Golden CO; USA Reliable Mine Water Technology IMWA 2013 A Rational Approach to Pit Lake Chemistry Modeling Dave Richers¹, C. Doc Richardson¹, Amy L. Hudson² ¹Tetra Tech, 350 Indiana Street, Suite 500,
More informationChapter 2 Thermal Processes
Chapter 2 Thermal Processes Hua Shao, Tianyuan Zheng, Philipp Hein, Son Nguyen, Benoit Garitte, Thomas Nagel and Haibing Shao 2.1 Transient Heat Conduction in a Transversal Isotropic Porous Medium Son
More informationLinking GoldSim with the PHREEQC Geochemical Model with a Dynamic Link Library Element
Linking GoldSim with the PHREEQC Geochemical Model with a Dynamic Link Library Element Ted Eary MWH Americas, Inc. Fort Collins, Colorado (ted.eary@mwhglobal.com or ted_eary@yahoo.com) GoldSim 2007 User
More informationMolecular Dynamics Simulations
Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline LAMMPS AdHiMad Lab 2 What is LAMMPS? LMMPS = Large-scale Atomic/Molecular Massively Parallel Simulator Open-Source (http://lammps.sandia.gov)
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 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 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 informationSIMULATION OF A TIME DEPENDENT 2D GENERATOR MODEL USING COMSOL MULTIPHYSICS
SIMULATION OF A TIME DEPENDENT 2D GENERATOR MODEL USING COMSOL MULTIPHYSICS Kazi Shamsul Arefin,Pankaj Bhowmik, Mohammed Wahiduzzaman Rony and Mohammad Nurul Azam Department of Electrical & Electronic
More informationContents Preface Introduction Model Concepts
Preface xi 1 Introduction 1 1.1 Environmental Problems and Geochemical Modeling 1 1.1.1 High-Level Radioactive Waste Disposal 1 1.1.2 Mining Related Environmental Issues 4 1.1.3 Landfills 8 1.1.4 Deep
More informationPerm State University Research-Education Center Parallel and Distributed Computing
Perm State University Research-Education Center Parallel and Distributed Computing A 25-minute Talk (S4493) at the GPU Technology Conference (GTC) 2014 MARCH 24-27, 2014 SAN JOSE, CA GPU-accelerated modeling
More informationQGIS FLO-2D Integration
EPiC Series in Engineering Volume 3, 2018, Pages 1575 1583 Engineering HIC 2018. 13th International Conference on Hydroinformatics Karen O Brien, BSc. 1, Noemi Gonzalez-Ramirez, Ph. D. 1 and Fernando Nardi,
More informationDiffusion and Reaction in Fe-Based Catalyst for Fischer- Tropsch Synthesis Using Micro Kinetic Rate Expressions
Diffusion and Reaction in Fe-Based Catalyst for Fischer- Tropsch Synthesis Using Micro Kinetic Rate Expressions 3-D CFD Model for Shell & Tube Exchanger with 7 Tubes Ender Ozden and Ilker Tari (2010) Multitubular
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 informationROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT
ROLE OF PORE-SCALE HETEROGENEITY ON REACTIVE FLOWS IN POROUS MATERIALS: VALIDITY OF THE CONTINUUM REPRESENTATION OF REACTIVE TRANSPORT PETER C. LICHTNER 1, QINJUN KANG 1 1 Los Alamos National Laboratory,
More informationVector Lane Threading
Vector Lane Threading S. Rivoire, R. Schultz, T. Okuda, C. Kozyrakis Computer Systems Laboratory Stanford University Motivation Vector processors excel at data-level parallelism (DLP) What happens to program
More informationSYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES
SYMBOLIC AND NUMERICAL COMPUTING FOR CHEMICAL KINETIC REACTION SCHEMES by Mark H. Holmes Yuklun Au J. W. Stayman Department of Mathematical Sciences Rensselaer Polytechnic Institute, Troy, NY, 12180 Abstract
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 informationDynamic Simulation Using COMSOL Multiphysics for Heterogeneous Catalysis at Particle Scale
Dynamic Simulation Using COMSOL Multiphysics for Heterogeneous Catalysis at Particle Scale Ameer Khan Patan *1, Mallaiah Mekala 2, Sunil Kumar Thamida 3 Department of Chemical Engineering, National Institute
More informationMultiphysics Coupling: Hypersonic Flow with Radiation and Ablation
Multiphysics Coupling: Hypersonic Flow with Radiation and Ablation Current Results and Future Strategies Paul T. Bauman, Roy H. Stogner The University of Texas at Austin February 20, 2009 Paul T. Bauman,
More informationDistinct Element Modeling of Coupled Chemo-Mechanical Compaction of Rock Salt
Distinct Element Modeling of Coupled Chemo-Mechanical Compaction of Rock Salt Ki-Bok Min 1, André Niemeijer 1, Derek Elsworth 1, Chris Marone 1 Department of Energy and Mineral Engineering Department of
More informationEnergy Pile Simulation - an Application of THM- Modeling
Energy Pile Simulation - an Application of THM- Modeling E. Holzbecher 1 1 Georg-August University, Göttingen, Germany Abstract Introduction: Energy piles, i.e. heat exchangers located within the foundation
More informationFracture-Matrix Flow Partitioning and Cross Flow: Numerical Modeling of Laboratory Fractured Core Flood
Fracture-Matrix Flow Partitioning and Cross Flow: Numerical Modeling of Laboratory Fractured Core Flood R. Sanaee *, G. F. Oluyemi, M. Hossain, and M. B. Oyeneyin Robert Gordon University *Corresponding
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 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 informationENHANCEMENT OF COMPUTER SYSTEMS FOR CANDU REACTOR PHYSICS SIMULATIONS
ENHANCEMENT OF COMPUTER SYSTEMS FOR CANDU REACTOR PHYSICS SIMULATIONS E. Varin, M. Dahmani, W. Shen, B. Phelps, A. Zkiek, E-L. Pelletier, T. Sissaoui Candu Energy Inc. WORKSHOP ON ADVANCED CODE SUITE FOR
More informationThe Augmented Spherical Wave Method
Introduction Institut für Physik, Universität Augsburg Electronic Structure in a Nutshell Outline 1 Fundamentals Generations 2 Outline 1 Fundamentals Generations 2 Outline Fundamentals Generations 1 Fundamentals
More informationWRF Modeling System Overview
WRF Modeling System Overview Louisa Nance National Center for Atmospheric Research (NCAR) Developmental Testbed Center (DTC) 27 February 2007 1 Outline What is WRF? WRF Modeling System WRF Software Design
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 informationIntroduction to Benchmark Test for Multi-scale Computational Materials Software
Introduction to Benchmark Test for Multi-scale Computational Materials Software Shun Xu*, Jian Zhang, Zhong Jin xushun@sccas.cn Computer Network Information Center Chinese Academy of Sciences (IPCC member)
More informationDesign and implementation of a new meteorology geographic information system
Design and implementation of a new meteorology geographic information system WeiJiang Zheng, Bing. Luo, Zhengguang. Hu, Zhongliang. Lv National Meteorological Center, China Meteorological Administration,
More informationAccess tofaydrological data from GIS applications by graphical software tools an example from the Hydrological Atlas of Germany (HAD)
Remote Sensing and Geographic Information Systems for Design and Operation of Water Resources Systems (Proceedings of Rabat Symposium S3, April 1997). IAHS Publ. no. 242, 1997 255 Access tofaydrological
More informationComputational Study of Chemical Kinetics (GIDES)
Computational Study of Chemical Kinetics (GIDES) Software Introduction Berkeley Madonna (http://www.berkeleymadonna.com) is a dynamic modeling program in which relational diagrams are created using a graphical
More informationCO 2 storage capacity and injectivity analysis through the integrated reservoir modelling
CO 2 storage capacity and injectivity analysis through the integrated reservoir modelling Dr. Liuqi Wang Geoscience Australia CO 2 Geological Storage and Technology Training School of CAGS Beijing, P.
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 informationSTEAMEST: A Software Tool for Estimation of Physical Properties of Water and Steam
226 JOURNAL OF SOFTWARE, VOL. 4, NO. 3, MAY 2009 STEAMEST: A Software Tool for Estimation of Physical Properties of Water and Steam Muhammad Faheem Department of Chemical Engineering, University of Engineering
More informationDistribution of radionuclides in soils modelling the dependence on soil parameters
Landesmessstelle für Radioaktivität Institut für Umweltphysik Fachbereich Physik/Elektrotechnik Distribution of radionuclides in soils modelling the dependence on soil parameters Volker Hormann and Helmut
More informationNumerical Analysis And Experimental Verification of a Fire Resistant Overpack for Nuclear Waste
Numerical Analysis And Experimental Verification of a Fire Resistant Overpack for Nuclear Waste Author P. Geraldini 1, Co-author A. Lorenzo 2 1,2 Sogin S.p.A., Rome, Italy Via Marsala 51/C 00185 Rome 1
More informationComparison of Heat and Mass Transport at the Micro-Scale
Comparison of Heat and Mass Transport at the Micro-Scale Ekkehard Holzbecher, Sandra Oehlmann October 10 th, 2012 Excerpt from the Proceedings of the 2012 COMSOL Conference in Milan Heat & Mass Transfer
More information