Exact and numerical solutions to the moving boundary problem resulting from reversible heterogeneous reactions and aqueous diffusion in a porous medium

Numerical finite difference calculations for diffusional mass transport incorporating reactions between minerals and a fluid phase, based on a continuum representation of porous media, are compared with the exact solution. The finite difference algorithm is based on the weak formulation of the moving boundary problem in which a fixed grid of node points is used. Mineral reactions are considered to be in local equilibrium with a fluid phase and may take place either at sharp reaction fronts or distributed homogeneously throughout a control volume. The theory is applied to one- and two-component systems. Analytical solutions to the finite difference equations for the first and second node points provide for a detailed comparison with the exact solution. It was found that on a time scale that is small compared with the time required for the solid to completely dissolve at a single node point, the finite difference approximation yields a spurious behavior for the concentration and solid phase volume fraction. However, the finite difference algorithm reproduces the average behavior of the motion of the reaction front and concentration of the reacting species provided the advance of the front is sufficiently slow resulting in a quasi-steady state solution. A numerical example ismore » presented for the dissolution of quartz at 550 /sup 0/C and 1000 bars to illustrate the general theory.« less

[1]  Mark Cross,et al.  An explicit numerical method to , 1983 .

[2]  F. Frank.  Radially symmetric phase growth controlled by diffusion , 1950, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[3]  Peter C. Lichtner,et al.  Continuum model for simultaneous chemical reactions and mass transport in hydrothermal systems , 1985 .

[4]  J. Weare,et al.  Diffusion metasomatism and mineral reaction zones; general principles and application to feldspar alteration , 1976 .

[5]  G. W. Fisher.  Nonequilibrium thermodynamics as a model for diffusion controlled metamorphic processes; errata , 1973 .

[6]  H. Mao,et al.  Bimetasomatism resulting from intergranular diffusion; I, A theoretical model for monomineralic reaction zone sequences , 1976 .

[7]  P. Lichtner,et al.  Interdiffusion with multiple precipitation/dissolution reactions: Transient model and the steady-state limit , 1986 .

[8]  H. Mao,et al.  Bimetasomatism resulting from intergranular diffusion; II, Prediction of multimineralic zone sequences , 1979 .

[9]  J. Ildefonse,et al.  Experimental study of silica diffusion during metasomatic reactions in the presence of water at 550°C and 1000 bars , 1976 .

[10]  R. Joesten.  Evolution of mineral assemblage zoning in diffusion metasomatism , 1977 .

[11]  M. Cross,et al.  Accurate solutions of moving boundary problems using the enthalpy method , 1981 .

[12]  A. Katchalsky,et al.  A simple model of interdiffusion with precipitation , 1970 .

[13]  H. Helgeson,et al.  Calculation of the thermodynamic properties of aqueous silica and the solubility of quartz and its polymorphs at high pressures and temperatures , 1977 .

[14]  James B. Thompson.  Geochemical reaction and open systems , 1970 .