Abstract
A numerical method for the calculation of equilibrium distributions of chemical species in aqueous solutions of electrolytes is presented. This method is constructed by transforming the problem of determining the set of unknown concentrations satisfying the mass balance and mass action equations into the equivalent problem of finding the limits of certain well-behaved mathematical sequences in a multivariable direct iteration scheme. The total (analytical) concentrations are taken as the starting estimates, and the recursive equations are constructed from the starting equations. It is shown that the sequences so constructed are both monotonic and bounded, hence convergent in an orderly fashion by a mathematical axiom. Each unknown free ion concentration (that is, each limit) is approached simultaneously from above and from below, being effectively “sandwiched” in an ever tightening manner. Strict error bounds therefore are easily constructed. The method has been found to be exceedingly efficient in practice, with the error bound (maximum fractional error in taking the upper estimate as the final answer) decreasing in an approximately exponential fashion with respect to iteration number, commonly 0.5–1.5 orders of magnitude per iteration. Correction for nonideality is presented, and the possibility of this giving rise to more than one solution is discussed in connection with three examples of natural waters of low-, medium-, and high-ionic strengths.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
References
Amit, O., and Bentor, Y. K., 1971, pH dilution curves of saline waters: Chem. Geology, v. 7, no. 4, p. 307–313.
Bos, M., and Meershoek, H. Q. J., 1972, A computer program for the calculation of equilibrium in complex systems: Anal. Chim. Acta, v. 61, no. 2, p. 185–199.
Brinkley, S. R., 1947, Calculation of the equilibrium composition of systems of many constituents: Jour. Phys. Chem., v. 51, p. 107–110.
Carpenter, A. B., and Miller, J. C., 1969, Geochemistry of saline subsurface water, Saline County (Missouri): Chem. Geology, v. 4, nos. 1/2, p. 135–167.
Garrels, R. M., and Christ, C. L., 1965, Solutions, minerals, and equilibria: Harper and Row, New York, 450 p.
Garrels, R. M., and Thompson, M. E., 1962, A chemical model for sea water at 25°C and one atmosphere total pressure: Am. Jour. Sci., v. 260, no. 1, p. 57–66.
Hamming, R. W., 1962, Numerical methods for scientists and engineers: McGraw-Hill, New York, 411 p.
Hancock, H. J., and Motzkin, T. S., 1960, Analysis of the mathematical model for chemical equilibrium,in Proc. 1st Conf. High Temp. Systems, p. 82–89.
Henrikson, M.. and Lees, M., 1967, Calculus: with an introduction to numerical methods, differential equations, and linear algebra (prelim. ed., pt. I): Worth Publishers, New York, 343 p.
Ingri, N., Kakolowicz, W., Sillen, L. G., and Warnquist, B., 1967, High-speed computers as a supplement to graphical methods—V. HALTAFALL, a general program for calculating the composition of equilibrium mixtures: Talanta, v. 14, no. 10, p. 1261–1286.
Karpov, I. K., and Kaz'min, L. A., 1972, Calculation of geochemical equilibria in heterogeneous multicomponent systems: Geochemistry Intern., v. 9, no. 2, p. 252–262. (Translated from Geokhimiya, v. 4, p. 402–414.)
Krieger, F. J., and White, W. B., 1948, A simplified method for computing the equilibrium composition of gaseous systems: Jour. Chem. Phys., v. 16, no. 4, p. 358–360.
Morel, F., and Morgan, J., 1972, A numerical method for computing equilibria in aqueous chemical systems: Environ. Sci. Technol., v. 6, p. 58–67.
Perrin, D. D., 1965, Multiple equilibria in assemblages of metal ions and complexing species: A model for biological systems: Nature, v. 206, no. 4980, p. 170–171.
Robinson, R. R., and Stokes, R. H., 1959, Electrolyte solutions (2nd ed. revised): Butterworths Publ. Ltd., London, 550 p.
Shapiro, N. Z., and Shapley, L. S., 1965, Mass action laws and the Gibbs free energy function: Jour. Soc. Ind. Appl. Math., v. 13, p. 353–375.
Ting-Po, I., and Nancollas, G. H., 1972, EQUIL—a general computational method for the calculation of solution equilibria: Anal. Chem., v. 44, p. 1940–1950.
Van Zeggeren, F., and Storey, S. H., 1970, The computation of chemical equilibria: Cambridge Univ. Press, Cambridge, England, 176 p.
Walters, L. J., Jr., and Wolery, T. J., 1975, A monotone-sequences algorithm and FORTRAN IV program for calculation of equilibrium distributions of chemical species: Computers and Geosciences, v. 1, no. 1 (in press).
White, W. B., Johnson, S. M., and Dantzig, G. B., 1958, Chemical equilibria in complex mixtures: Jour. Chem. Phys., v. 28, no. 5, p. 751–755.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Wolery, T.J., Walters, L.J. Calculation of equilibrium distributions of chemical species in aqueous solutions by means of monotone sequences. Mathematical Geology 7, 99–115 (1975). https://doi.org/10.1007/BF02080487
Received:
Issue Date:
DOI: https://doi.org/10.1007/BF02080487