TY - JOUR
T1 - Intercomparison of 3D pore-scale flow and solute transport simulation methods
AU - Yang, Xiaofan
AU - Mehmani, Yashar
AU - Perkins, William A.
AU - Pasquali, Andrea
AU - Schönherr, Martin
AU - Kim, Kyungjoo
AU - Perego, Mauro
AU - Parks, Michael L.
AU - Trask, Nathaniel
AU - Balhoff, Matthew T.
AU - Richmond, Marshall C.
AU - Geier, Martin
AU - Krafczyk, Manfred
AU - Luo, Li Shi
AU - Tartakovsky, Alexandre M.
AU - Scheibe, Timothy D.
N1 - Funding Information:
Sandia National Laboratories is a multiprogram laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-AC04-94AL85000 . This SPH work was supported by the DOE Office of Science Advanced Scientific Computing Research (ASCR) Applied Mathematics program as part of the Collaboratory on Mathematics for Mesoscopic Modeling of Materials (CM4). This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.
Funding Information:
Research at Pacific Northwest National Laboratory (PNNL) was supported by the U. S. Department of Energy (DOE) Office of Biological and Environmental Research (BER) through the PNNL Subsurface Biogeochemical Research Scientific Focus Area project. Computations described here were performed using computational facilities of the Environmental Molecular Sciences Laboratory (EMSL), a national scientific user facility sponsored by DOE-BER and located at PNNL, computational facilities of PNNL’s Institutional Computing program, and computational facilities of the National Energy Research Supercomputing Center, which is supported by the DOE Office of Science under Contract No. DE-AC02-05CH11231. PNNL is operated for the DOE by Battelle Memorial Institute under Contract No. DE-AC06-76RLO 1830 .
Funding Information:
The pore-network modeling was carried out under funding from the Center for Frontiers of Subsurface Energy Security, an Energy Frontier Research Center funded by the U.S. Department of Energy , Office of Science, Office of Basic Energy Sciences, DOE Project No. DE-SC0001114 . We would also like to thank Karsten Thompson and LSU for providing the network extraction algorithm.
Funding Information:
Institute for Computational Modelling in Civil Engineering (iRMB) at Technische Universität Braunschweig gratefully acknowledges financial support by the Deutsche Forschungsgemeinschaft (DFG) for funding the Research Training Group MUSIS (FOR 1083).
Publisher Copyright:
© 2015
PY - 2016/9/1
Y1 - 2016/9/1
N2 - Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based on the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This study provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods.
AB - Multiple numerical approaches have been developed to simulate porous media fluid flow and solute transport at the pore scale. These include 1) methods that explicitly model the three-dimensional geometry of pore spaces and 2) methods that conceptualize the pore space as a topologically consistent set of stylized pore bodies and pore throats. In previous work we validated a model of the first type, using computational fluid dynamics (CFD) codes employing a standard finite volume method (FVM), against magnetic resonance velocimetry (MRV) measurements of pore-scale velocities. Here we expand that validation to include additional models of the first type based on the lattice Boltzmann method (LBM) and smoothed particle hydrodynamics (SPH), as well as a model of the second type, a pore-network model (PNM). The PNM approach used in the current study was recently improved and demonstrated to accurately simulate solute transport in a two-dimensional experiment. While the PNM approach is computationally much less demanding than direct numerical simulation methods, the effect of conceptualizing complex three-dimensional pore geometries on solute transport in the manner of PNMs has not been fully determined. We apply all four approaches (FVM-based CFD, LBM, SPH and PNM) to simulate pore-scale velocity distributions and (for capable codes) nonreactive solute transport, and intercompare the model results. Comparisons are drawn both in terms of macroscopic variables (e.g., permeability, solute breakthrough curves) and microscopic variables (e.g., local velocities and concentrations). Generally good agreement was achieved among the various approaches, but some differences were observed depending on the model context. The intercomparison work was challenging because of variable capabilities of the codes, and inspired some code enhancements to allow consistent comparison of flow and transport simulations across the full suite of methods. This study provides support for confidence in a variety of pore-scale modeling methods and motivates further development and application of pore-scale simulation methods.
UR - http://www.scopus.com/inward/record.url?scp=84951818341&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84951818341&partnerID=8YFLogxK
U2 - 10.1016/j.advwatres.2015.09.015
DO - 10.1016/j.advwatres.2015.09.015
M3 - Article
AN - SCOPUS:84951818341
SN - 0309-1708
VL - 95
SP - 176
EP - 189
JO - Advances in Water Resources
JF - Advances in Water Resources
ER -