| ex9 {phreeqc} | R Documentation |
Example 9–Kinetic Oxidation of Dissolved Ferrous Iron With Oxygen
Description
Kinetic rate expressions can be defined in a completely general
way in PHREEQC by using Basic statements in the RATES data block. The rate
expressions can be used in batch-reaction and transport calculations with
the KINETICS data block. For transport calculations (ADVECTION or TRANSPORT),
kinetic reactions can be defined cell by cell by the number range following
the KINETICS keyword (KINETICS m-n). The rate expressions are integrated with
an embedded (up to) 5th-order Runge-Kutta-Fehlberg algorithm, or with a
stiff, variable-order, variable-step multistep solver (Cohen and Hindmarsh,
1996). Equilibrium is calculated before a kinetic calculation is initiated
and again when a kinetic reaction increment is added. Equilibrium includes
solution species equilibrium; exchange-, equilibrium-phase-, solid-solution-,
and surface-assemblage equilibrium; and gas-phase equilibrium. A check is
performed to ensure that the difference between estimates of the integrated
rate over a time interval is smaller than a user-defined tolerance. If the
tolerance is not satisfied, then the integration over the time interval is
automatically restarted with a smaller time interval. The example can be run
using the phrRunString routine.
Source
http://wwwbrr.cr.usgs.gov/projects/GWC_coupled/phreeqc
References
https://pubs.usgs.gov/tm/06/a43/pdf/tm6-A43.pdf
See Also
Other Examples:
ex10,
ex11,
ex12,
ex13a,
ex14,
ex15,
ex16,
ex17,
ex18,
ex19,
ex1,
ex20a,
ex21,
ex22,
ex2,
ex3,
ex4,
ex5,
ex6,
ex7,
ex8
Examples
phrLoadDatabaseString(phreeqc.dat)
phrSetOutputStringsOn(TRUE)
phrRunString(ex9)
phrGetOutputStrings()