Mixed transport reaction control of gypsum dissolution kinetics in aqueous solutions and initiation of gypsum karst, 1997, Raines M. A. , Dewers T. A. ,

Experiments with gypsum in aqueous solutions at 25 degrees C, low ionic strengths, and a range of saturation states indicate a mixed surface reaction and diffusional transport control of gypsum dissolution kinetics. Dissolution rates were determined in a mixed flow/rotating disc reactor operating under steady-state conditions, in which polished gypsum discs were rotated at constant speed and reactant solutions were continuously fed into the reactor. Rates increase with velocity of spin under laminar conditions (low rates of spin), but increase asymptotically to a constant rate as turbulent conditions develop with increasing spin velocity, experiencing a small jump in magnitude across the laminar-turbulent transition. A Linear dependence of rates on the square root of spin velocity in the laminar regime is consistent with rates being limited by transport through a hydrodynamic boundary layer. The increase in rate with onset of turbulence accompanies a near discontinuous drop in hydrodynamic boundary layer thickness across the transition. A relative independence of rates on spinning velocity in the turbulent regime plus a nonlinear dependence of rates on saturation state are factors consistent with surface reaction control. Together these behaviors implicate a 'mixed' transport and reaction control of gypsum dissolution kinetics. A rate law which combines both kinetic mechanisms and can reproduce experimental results under laminar flow conditions is proposed as follows: R = k(t) {1 - Omega(b)() zeta [1 - (1 2(1 - Omega(b)())(1/2)]} where k(t) is the rate coefficient for transport control, and Omega(b)() is the mean ionic saturation state of the bulk fluid. The dimensionless parameter zeta(=Dm(eq)()/2 delta k() where m(eq)() = mean ionic molal equilibrium concentration, D is the diffusion coefficient through the hydrodynamic boundary layer, delta equals the boundary layer thickness and k() is the rate constant for surface reaction control) indicates which process, transport or surface reaction, dominates, and is sensitive to the hydrodynamic conditions in the reactor. For the range of conditions used in our experiments, zeta varies from about 1.4 to 4.5. Rates of gypsum dissolution were also determined in situ in a cavern system in the Permian Blaine Formation, southwestern Oklahoma. Although the flow conditions in the caverns were not determinable, there is good agreement between lab- and field-determined rates in that field rate magnitudes lie within a range of rates determined experimentally under zero to low spin velocities A numerical model coupling fluid flow and gypsum reaction in an idealized circular conduit is used to estimate the distance which undersaturated solutions will travel into small incipient conduits before saturation is achieved. Simulations of conduit wall dissolution showed-member behaviors of conduit formation and surface denudation that depend on flow boundary conditions (constant discharge or constant hydraulic gradient and initial conduit radius. Surface-control of dissolution rates. which becomes more influential with higher fluid flow velocity, has the effect that rate decrease more slowly as saturation is approached than otherwise would occur if rates were controlled by transport alone. This has the effect that reactive solutions can penetrate much farther into gypsum-bearing karst conduits than heretofore thought possible, influencing timing and mechanism of karst development as well as stability of engineered structures above karst terrain

An experimental study of calcite and limestone dissolution rates as a function of pH from -1 to 3 and temperature from 25 to 80 degrees C, 1998, Alkattan M, Oelkers Eh, Dandurand Jl, Schott J,

Dissolution rates of single calcite crystals, limestones, and compressed calcite powders were determined from sample weight loss using free-drift rotating disk techniques. Experiments were performed in aqueous HCl solutions over the bulk solution pH range -1 to 3, and at temperatures of 25 degrees, 50 degrees, and 80 degrees C. Corresponding rates of the three different sample types are identical within experimental uncertainty. Interpretation of these data using equations reported by Gregory and Riddiford [Gregory, D.P., Riddiford, A.C., 1956. Transport to the surface of a rotating disc. J. Chem. Sec. London 3, 3756-3764] yields apparent rate constants and H diffusion coefficients. The logarithms of overall calcite dissolution rates (r) obtained at constant disk rotation speed are inversely proportional to the bulk solution pH, consistent with r = k(2') a(H,b), where k(2)' stands for an apparent rate constant and a(H,b) designates the hydrogen ion activity in the bulk solution, This variation of dissolution rates with pH is consistent with corresponding rates reported in the literature and the calcite dissolution mechanism reported by Wollast [Wollast, R., 1990. Rate and mechanism of dissolution of carbonates in the system CaCO3-MgCO3. In: Stumm, W. (Ed.), Aquatic Chemical Kinetics. Wiley, pp. 431-445]. Apparent rate constants for a disk rotation speed of 340 rpm increase from 0.07 0.02 to 0.25 0.02 mol m(-2) s(-1) in response to increasing temperature from 25 degrees to 80 degrees C. H diffusion coefficients increase from (2.9 to 9.2) x 10(-9) m(2) s(-1) over this temperature range with an apparent activation energy of 19 kJ mol(-1). (C) 1998 Elsevier Science B.V. All rights reserved