Mixing‐induced precipitation: Experimental study and multiscale numerical analysis

Water Resources Research - Tập 44 Số 6 - 2008
Alexandre M. Tartakovsky1, G. D. Redden2, Peter C. Lichtner3, Tim Scheibe1, Paul Meakin2,4,5
1Pacific Northwest National Laboratory , Richland, Washington, USA
2Idaho National Laboratory, Idaho Falls, Idaho, USA
3Los Alamos National Laboratory, Los Alamos, New Mexico USA
4Multiphase Flow Assurance Innovation Center, Institute for Energy Technology, Kjeller, Norway
5Physics of Geological Processes, University of Oslo, Oslo, Norway

Tóm tắt

Laboratory experiments, pore‐scale simulations, and continuum (Darcy) ‐scale simulations were used to study mixing‐induced precipitation in porous media. In the experimental investigation, solutions containing Na2CO3 and CaCl2 were each injected in different halves of a quasi two‐dimensional flow cell filled with quartz sand. As a result of the in situ mixing between the two solutions, a narrow calcite precipitate layer (less than 5 mm wide) of more or less uniform width was formed between the individual solutions. Pore‐scale simulations were conducted to help understand the mechanism of precipitation layer formation. The effect of the Peclet number, Pe, and the Damköhler number, Da, on mixing induced precipitation was also investigated. Pore‐scale simulations revealed the presence of large pore‐scale concentration gradients. This, and the presence of features, such as the precipitation layer, with characteristic lengths on the order of the average sand grain diameter, indicate the absence of a clear scale separation required for the strict derivation of Darcy‐scale advection‐dispersion equations. Nevertheless, we found that an adaptive high‐resolution model based on advection‐dispersion equations with grid sizes in the mixing zone smaller than the size of the sand grains can qualitatively reproduce the essential features of the experiment. As an alternative to computationally expensive high‐resolution simulations, we proposed new forms for the homogeneous and heterogeneous reaction terms in Darcy‐scale advection dispersion equations. These terms involve transport and mixing indices that account for highly nonuniform pore‐scale concentration distributions and highly localized reactions. The proposed model accurately estimates the changes in solute concentrations due to homogenous and heterogeneous reactions during precipitation of minerals, observed in the pore‐scale simulations, while conventional low‐resolution advective‐dispersion equations produced erroneous results.

Từ khóa


Tài liệu tham khảo

Allen M. P., 2001, Computer Simulation of Liquids

10.1103/PhysRevLett.72.1384

10.1016/j.advwatres.2004.11.010

Gingold R. A., 1997, Smoothed particle hydrodynamics: Theory and application to non‐spherical stars, Mon. Not. R. Astron. Soc., 181

Krauskopf K. B., 1967, Introduction to Geochemistry

Lichtner P. C., 1996, Reactive Transport in Porous Media, 10.1515/9781501509797

Lichtner P. C.(2001) FLOTRAN user's manual report Los Alamos Natl. Lab. Los Alamos N. M.

10.1086/112164

10.1016/S0016-7037(03)00372-7

10.1088/0034-4885/68/8/R01

10.1006/jcph.1997.5776

Schoenberg I. J., 1946, Contributions to the problem of approximation of equidistant data by analytical functions: Part A, Q. Appl. Math., 4

10.1103/PhysRevE.72.026301

10.2136/vzj2004.0178

10.1016/j.jcp.2005.02.001

10.1016/j.advwatres.2005.11.014

10.1016/j.jcp.2006.08.013

10.1029/2005WR004770

10.1002/aic.690130308

Whitaker S., 1981, Introduction to Fluid Mechanics

10.1007/978-94-017-3389-2

10.1023/A:1010769915901

10.1006/jcph.2002.7189

10.1002/(SICI)1096-9853(19990810)23:9<881::AID-NAG996>3.0.CO;2-K

Thông tin
Thông tin xuất bản

Water Resources Research

Tập 44 Số 6

Thông tin tác giả