On simulating reactive chemical transport with dominating precipitated species controlling chemical equilibrium

Computational Geosciences - Tập 2 - Trang 151-170 - 1998
Hwai‐Ping Cheng1, Gour‐Tsyh Yeh
1Department of Civil and Environmental Engineering, The Pennsylvania State University, University Park, USA

Tóm tắt

We present an approach developed to compute chemical equilibrium and its corresponding reactive chemical transport when dominating precipitated species (DPS) exist. In computing chemical equilibrium, most models take the concentrations or activities of component species and precipitated species as the master variables. However, when the amount of a precipitated species is much larger than those of other species, small computational errors on this DPS concentration might introduce large errors on the concentrations of other species and would cause non‐mass‐conserved numerical results. To deal with the existence of DPS, we pick as master variables the concentration change, rather than the concentration, of DPS to compute chemical equilibrium. Since the concentration changes of DPS will no longer be much larger than the concentrations of other species in determining equilibrium, our approach is able to provide correct numerical results. We also employ the modified total analytical concentrations, rather than the total analytical concentrations, of aqueous components as the dependent variables in presenting and solving corresponding transport equations. Several examples are given to reveal the numerical problems associated with DPS and to demonstrate the success of our approach.

Tài liệu tham khảo

C.M. Bethke, Geochemical Reaction Modeling (Oxford University Press, Oxford, 1996). H.P. Cheng, Development and application of a three-dimensional finite element model of subsurface flow, heat transfer, and reactive chemical transport, Ph.D. thesis, Department of Civil and Environmental Engineering, The Pennsylvania State University, University Park, PA (1995). H.P. Cheng and G.T. Yeh, An approach on simulating reactive chemical transport with precipitated species dominating chemical equilibrium, in: Proc. 10th International Conference on Computational Methods in Water Resources, Heidelberg, Germany (July 19–22, 1994). P. Deuflhard, Global inexact Newton methods for very large scale nonlinear problems, Impact of Computing in Science and Engineering 3 (1991) 366–393. A.R. Felmy, D. Girvin and E.A. Jenne, MINTEQ: A computer program for calculating aqueous geochemical equilibria, Report, U.S. Environ. Prot. Agency, Washington, DC (1983). R.Y. James and J. Rubin, Applicability of the local equilibrium assumption to transport through soil of solute affected by ion exchange, in: Chemical Modeling in Aqueous Systems, ed. E.A. Jenne (ACS, Washington, DC, 1971). C.T. Kelley, Iterative Methods for Linear and Nonlinear Equations (SIAM, Philadelphia, PA, 1995). P.C. Lichtner, Continuum formulation of multicomponent-multiphase reactive transport, in: Reactive Transport in Porous Media, Reviews in Mineralogy, Vol. 34, eds. P.C. Lichtner et al. (The Mineralogical Society of America, 1996) Chapter 1, pp. 1–81. D.C. Mangold and C.F. Tsang, A summary of subsurface hydrological and hydrochemical models, Reviews of Geophysics 29(1) (1991) 51–79. J. Šimůnek and D.L. Suarez, Two-dimensional transport model for variably saturated porous media with major ion chemistry, Water Resources Research 30(4) (1994) 1115–1133. W. Stumm and J.J. Morgan, Aqueous Chemistry (Wiley, New York, 1981). T.J. Wolery, Calculation of chemical equilibrium between aqueous solution and minerals: The EQ3/6 software package, Report UCRL-52658, Lawrence Livermore National Laboratory, Livermore, CA (1979). G. Wu and S.T. Chieng, Modeling multicomponent reactive chemical transport in nonisothermal unsaturated/saturated soils, Part 1. Mathematical modeling development, Transactions of the ASAE 38(3) (1995) 817–826. G.T. Yeh, EQMOD: A chemical EQuilibrium MODel of complexation, adsorption, ion-exchange, precipitation/dissolution, oxidation/reduction, and acid-base reactions, Short course notes of modeling reactive multispecies multicomponents transport through subsurface media, The Pennsylvania State University, University Park, PA (November 1992). G.T. Yeh, J.R. Chang and T.E. Short, An exact peak capturing and essentially oscillation-free scheme to solve advection-dispersion-reactive transport equations, Water Resources Research 28(11) (1992) 2937–2951. G.T. Yeh, G. Iskra, J.M. Zachara and J.E. Szecsody, Development and verification of a mixed chemical kinetic and equilibrium model, Advances in Environmental Research 2(1) (1998) 24–56. G.T. Yeh and V.S. Tripathi, HYDROGEOCHEM: A coupled model of HYDROlogical transport and GEOCHEMical equilibria in reactive multicomponent systems, ORNL-6371, Oak Ridge National Laboratory, Oak Ridge, TN (1990).