Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

An investigation of the effect of pore scale flow on average geochemical reaction rates using direct numerical simulation

Journal Article · · Water Resources Research
DOI:https://doi.org/10.1029/2011WR011404· OSTI ID:1082190
 [1];  [2];  [1];  [2]
  1. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Earth Sciences Division
  2. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Computational Research Division

The scale-dependence of geochemical reaction rates hinders their use in continuum scale models intended for the interpretation and prediction of chemical fate and transport in subsurface environments such as those considered for geologic sequestration of CO2. Processes that take place at the pore scale, especially those involving mass transport limitations to reactive surfaces, may contribute to the discrepancy commonly observed between laboratory-determined and continuum-scale or field rates. In this study we investigate the dependence of mineral dissolution rates on the pore structure of the porous media by means of pore scale modeling of flow and multicomponent reactive transport. The pore scale model is composed of high-performance simulation tools and algorithms for incompressible flow and conservative transport combined with a general-purpose multicomponent geochemical reaction code. The model performs direct numerical simulation of reactive transport based on an operator-splitting approach to coupling transport and reactions. The approach is validated with a Poiseuille flow single-pore experiment and verified with an equivalent 1-D continuum-scale model of a capillary tube packed with calcite spheres. Using the case of calcite dissolution as an example, the high-resolution model is used to demonstrate that nonuniformity in the flow field at the pore scale has the effect of decreasing the overall reactivity of the system, even when systems with identical reactive surface area are considered. In conclusion, the effect becomes more pronounced as the heterogeneity of the reactive grain packing increases, particularly where the flow slows sufficiently such that the solution approaches equilibrium locally and the average rate becomes transport-limited.

Research Organization:
Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
Sponsoring Organization:
USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22); USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR) (SC-21); Energy Frontier Research Center (EFRC)
DOE Contract Number:
AC02-05CH11231
OSTI ID:
1082190
Report Number(s):
LBNL--5385E
Journal Information:
Water Resources Research, Journal Name: Water Resources Research Journal Issue: 3 Vol. 48; ISSN 0043-1397
Publisher:
American Geophysical Union (AGU)
Country of Publication:
United States
Language:
English

Similar Records

Mineral dissolution kinetics at the pore scale
Journal Article · Thu May 24 00:00:00 EDT 2007 · Geochimica et Cosmochimica Acta · OSTI ID:929022

Pore scale modeling of reactive transport involved in geologic CO2 sequestration
Journal Article · Wed Dec 31 23:00:00 EST 2008 · Transport in Porous Media · OSTI ID:956358

Upscaling geochemical reaction rates using pore-scale networkmodeling
Journal Article · Thu May 19 00:00:00 EDT 2005 · Advances in Water Resources · OSTI ID:902804