Influence of multi-species solute transport on modeling of hydrated Portland cement leaching in strong nitrate solutions

Abstract The differences between the predictions of Fickian diffusion and Nernst-Planck electro-diffusion based reactive transport models are investigated in the context of an important practical problem – leaching of Portland cement in strong nitrate porewater solutions under saturated conditions. Numerical results are presented for three distinct porewater solutions: deionized water, ammonium nitrate solution, and diluted porewater from a cementitious low-activity nuclear waste form. Solute concentrations predicted by the two transport models show significant differences: as much as a factor of two for ammonium nitrate leaching and as much as a factor of four for waste form leaching; whereas, negligible differences are observed for deionized water leaching. The major conclusion is that Fickian diffusion may be accurate enough for deionized water leaching, but electro-diffusion should be considered in the reactive solute transport modeling of ammonium nitrate decalcification and in the degradation modeling of cementitious materials exposed to nuclear waste form solutions.

[1]  Patrick Goblet,et al.  Module-oriented modeling of reactive transport with HYTEC , 2003 .

[2]  Philippe Blanc,et al.  Modeling diffusion of an alkaline plume in a clay barrier , 2004 .

[3]  Alexandra B. Ribeiro,et al.  Computing multi-species chemical equilibrium with an algorithm based on the reaction extents , 2013, Comput. Chem. Eng..

[4]  Konrad B. Krauskopf,et al.  Principles and applications of inorganic geochemistry: Gunter Faure. MacMillan, 1991, xiii + 626p., US $58.00 (ISBN 0-02-336441-6) , 1991 .

[5]  Charles Nicholson,et al.  Diffusion and related transport mechanisms in brain tissue , 2001 .

[6]  E. Samson,et al.  Modeling the transport of ions in unsaturated cement-based materials , 2007 .

[7]  J. Meeussen,et al.  ORCHESTRA: an object-oriented framework for implementing chemical equilibrium models. , 2003, Environmental science & technology.

[8]  Ung Jin Na,et al.  Service life prediction of concrete wharves with early-aged crack: Probabilistic approach for chloride diffusion , 2009 .

[9]  Lars-Olof Nilsson,et al.  Numerical simulation of multi-species transport through saturated concrete during a migration test - MsDiff code , 2000 .

[10]  G. Escadeillas,et al.  Use of ammonium nitrate solution to simulate and accelerate the leaching of cement pastes due to deionized water , 1997 .

[11]  J. Vinograd,et al.  Diffusion of Electrolytes and of the Ions in their Mixtures , 1941 .

[12]  Ray Luo,et al.  On removal of charge singularity in Poisson-Boltzmann equation. , 2009, The Journal of chemical physics.

[13]  Martin Z. Bazant,et al.  Homogenization of the Poisson-Nernst-Planck equations for Ion Transport in Charged Porous Media , 2012, SIAM J. Appl. Math..

[14]  C. Bethke.  Geochemical and Biogeochemical Reaction Modeling , 2007 .

[15]  J. Rasaiah,et al.  Solvent Structure, Dynamics, and Ion Mobility in Aqueous Solutions at 25 °C , 1998 .

[16]  A. Walter,et al.  Modeling of multicomponent reactive transport in groundwater: 1. Model development and evaluation , 1994 .

[17]  J. Strikwerda.  Finite Difference Schemes and Partial Differential Equations, Second Edition , 2004 .

[18]  S. Ben-Yaakov.  Diffusion of sea water ions—I. Diffusion of sea water into a dilute solution , 1972 .

[19]  A. Lasaga.  Multicomponent exchange and diffusion in silicates , 1979 .

[20]  J. Eijkel,et al.  Principles and applications of nanofluidic transport. , 2009, Nature nanotechnology.

[21]  Duncan Herfort,et al.  Thermodynamics and cement science , 2011 .

[22]  E. Oelkers,et al.  Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Aqueous tracer diffusion coefficients of ions to 1000°C and 5 kb , 1988 .

[23]  John Crank,et al.  The Mathematics Of Diffusion , 1956 .

[24]  Samson,et al.  Numerical Solution of the Extended Nernst-Planck Model. , 1999, Journal of colloid and interface science.

[25]  Thomas Wagner,et al.  GEM-Selektor geochemical modeling package: revised algorithm and GEMS3K numerical kernel for coupled simulation codes , 2012, Computational Geosciences.

[26]  A. Lasaga,et al.  The treatment of multi-component diffusion and ion pairs in diagenetic fluxes , 1979 .

[27]  Christophe Tournassat,et al.  Influence of reaction kinetics and mesh refinement on the numerical modelling of concrete/clay interactions , 2009 .

[28]  Carl I. Steefel,et al.  Multicomponent reactive transport in discrete fractures: I. Controls on reaction front geometry , 1998 .

[29]  S. Finsterle,et al.  Electrokinetic coupling in unsaturated porous media. , 2007, Journal of colloid and interface science.

[30]  Edward L Cussler,et al.  Diffusion: Mass Transfer in Fluid Systems , 1984 .

[31]  C. Appelo,et al.  Multicomponent diffusion of a suite of tracers (HTO, Cl, Br, I, Na, Sr, Cs) in a single sample of Opalinus Clay , 2010 .

[32]  S. Subramaniam,et al.  Treatment of electrostatic effects in proteins: Multigrid‐based newton iterative method for solution of the full nonlinear poisson–boltzmann equation , 1994, Proteins.

[33]  A. Revil,et al.  Diffusion of ions in unsaturated porous materials. , 2008, Journal of colloid and interface science.

[34]  B. Johannesson.  Comparison between the Gauss' law method and the zero current method to calculate multi-species ionic diffusion in saturated uncharged porous materials , 2010 .

[35]  C. Dekker,et al.  Surface-charge-governed ion transport in nanofluidic channels. , 2004, Physical review letters.

[36]  Sankaran Mahadevan,et al.  Numerical simulation of cementitious materials degradation under external sulfate attack , 2010 .

[37]  D. L. Parkhurst,et al.  User's guide to PHREEQC (Version 2)-a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations , 1999 .

[38]  B. Johannesson,et al.  Development of a multi-species mass transport model for concrete with account to thermodynamic phase equilibriums , 2011 .

[39]  J. Marchand,et al.  Calculation of ionic diffusion coefficients on the basis of migration test results , 2003 .

[40]  Jin Xia,et al.  Numerical simulation of ionic transport in cement paste under the action of externally applied electric field , 2013 .

[41]  M. Bazant,et al.  Steric effects in the dynamics of electrolytes at large applied voltages. II. Modified Poisson-Nernst-Planck equations. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[42]  V. S. Tripathi,et al.  A critical evaluation of recent developments in hydrogeochemical transport models of reactive multichemical components , 1989 .

[43]  Björn Johannesson,et al.  Multi-species ionic diffusion in concrete with account to interaction between ions in the pore solution and the cement hydrates , 2007 .

[44]  J. Bear,et al.  Introduction to Modeling of Transport Phenomena in Porous Media , 1990 .

[45]  J. H. Cushman,et al.  Multicomponent, Multiphase Thermodynamics of Swelling Porous Media with Electroquasistatics: I. Macroscale Field Equations , 1999 .

[46]  Henry S. Frank,et al.  Ion-solvent interaction. Structural aspects of ion-solvent interaction in aqueous solutions: a suggested picture of water structure , 1957 .

[47]  F. Glasser,et al.  Alkali binding in cement pastes: Part I. The C-S-H phase , 1999 .

[48]  P. Renaud,et al.  Transport phenomena in nanofluidics , 2008 .

[49]  Everett L. Shock,et al.  Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Correlation algorithms for ionic species and equation of state predictions to 5 kb and 1000°C , 1988 .

[50]  M. Kurnikova,et al.  Poisson-Nernst-Planck theory approach to the calculation of current through biological ion channels , 2005, IEEE Transactions on NanoBioscience.

[51]  P. Renaud,et al.  Ionic transport phenomena in nanofluidics: experimental and theoretical study of the exclusion-enrichment effect on a chip. , 2005, Nano letters.

[52]  J. H. Cushman,et al.  Multicomponent, Multiphase Thermodynamics of Swelling Porous Media with Electroquasistatics: II. Constitutive Theory , 2000 .

[53]  Karsten Pruess,et al.  TOUGHREACT - A simulation program for non-isothermal multiphase reactive geochemical transport in variably saturated geologic media: Applications to geothermal injectivity and CO2 geological sequestration , 2006, Comput. Geosci..

[54]  M. Bazant,et al.  Steric effects in the dynamics of electrolytes at large applied voltages. I. Double-layer charging. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[55]  Katrien Audenaert,et al.  On the time dependency of the chloride migration coefficient in concrete , 2010 .

[56]  Serdar Kuyucak,et al.  Models of permeation in ion channels , 2001 .

[57]  James J. Beaudoin,et al.  Modeling chemical activity effects in strong ionic solutions , 1999 .

[58]  E. Shock,et al.  Inorganic species in geologic fluids: correlations among standard molal thermodynamic properties of aqueous ions and hydroxide complexes. , 1997, Geochimica et cosmochimica acta.

[59]  Dave Easterbrook,et al.  Multi-phase modelling of ionic transport in concrete when subjected to an externally applied electric field , 2012 .

[60]  Yizhak Marcus,et al.  Effect of ions on the structure of water: structure making and breaking. , 2009, Chemical reviews.

[61]  B. Lothenbach.  Thermodynamic equilibrium calculations in cementitious systems , 2010 .

[62]  J. Molinero,et al.  On the relevance of electrochemical diffusion for the modeling of degradation of cementitious materials , 2010 .

[63]  Kristian Krabbenhoft,et al.  Application of the Poisson–Nernst–Planck equations to the migration test , 2008 .

[64]  Bruno Gérard,et al.  Simplified modelling of calcium leaching of concrete in various environments , 2002 .

[65]  E. Samson,et al.  Modeling Degradation of Cementitious Materials in Aggressive Aqueous Environments , 2013 .

[66]  C. Appelo,et al.  Multicomponent diffusion modeling in clay systems with application to the diffusion of tritium, iodide, and sodium in Opalinus Clay. , 2007, Environmental science & technology.

[67]  Ionel Michael Navon,et al.  VARIATM—A FORTRAN program for objective analysis of pseudostress wind fields using large-scale conjugate-gradient minimization , 1991 .

[68]  T. Luping,et al.  On the mathematics of time-dependent apparent chloride diffusion coefficient in concrete , 2007 .

[69]  B. Lothenbach,et al.  Thermodynamic modelling of the hydration of Portland cement , 2006 .

[70]  Manfred Rudolph,et al.  Digital simulations with the fast implicit finite difference (FIFD) algorithm: Part II. An improved treatment of electrochemical mechanisms with second-order reactions , 1992 .

[71]  Z. Bažant,et al.  Modeling Chloride Penetration in Saturated Concrete , 1999 .

[72]  N. Linde,et al.  Chemico-electromechanical coupling in microporous media. , 2006, Journal of colloid and interface science.

[73]  Kenneth A. Snyder,et al.  Effect of speciation on the apparent diffusion coefficient in nonreactive porous systems , 2001 .

[74]  Isaak Rubinstein.  Electro-diffusion of ions , 1987 .

[75]  B. Johannesson,et al.  Framework for reactive mass transport: Phase change modeling of concrete by a coupled mass transport and chemical equilibrium model , 2014 .

[76]  D. Kulik.  Improving the structural consistency of C-S-H solid solution thermodynamic models , 2011 .

[77]  B. Lothenbach,et al.  Thermodynamic Modelling of the Effect of Temperature on the Hydration and Porosity of Portland Cement , 2008 .

[78]  Christophe Tournassat,et al.  Modelling approaches for anion-exclusion in compacted Na-bentonite , 2011 .

[79]  J. Marchand,et al.  Describing ion diffusion mechanisms in cement-based materials using the homogenization technique , 1999 .

[80]  Barbara Lothenbach,et al.  Thermodynamic properties of Portland cement hydrates in the system CaO–Al2O3–SiO2–CaSO4–CaCO3–H2O , 2007 .

[81]  F. J. Pearson,et al.  Nagra/PSI Chemical Thermodynamic Data Base 01/01 , 2002 .