TY - JOUR
T1 - Lattice-Boltzmann simulations of the capillary pressure-saturation-interfacial area relationship for porous media
AU - Porter, Mark L.
AU - Schaap, Marcel G.
AU - Wildenschild, Dorthe
N1 - Funding Information:
Porter and Wildenschild were supported, in part, by NSF-EAR-06101108 and NSF-EAR-0337711. Schaap was supported, in part, by NSF-EAR-0337378 and NSF-EAR-061003. A portion of this work was conducted at GeoSoilEnviroCARS (Sector 13), Advanced Photon Source (APS), Argonne National Laboratory. GeoSoilEnviroCARS is supported by the National Science Foundation – Earth Sciences ( EAR-0622171 ) and Department of Energy – Geosciences ( DE-FG02-94ER14466 ). Use of the Advanced Photon Source was supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences , under Contract No. DE-AC02-06CH11357 .
PY - 2009/11
Y1 - 2009/11
N2 - Hysteresis in the relationship between capillary pressure (Pc), wetting phase saturation (Sw) and nonwetting-wetting interfacial area per volume (anw) is investigated using multiphase lattice-Boltzmann simulations of drainage and imbibition in a glass bead porous system. In order to validate the simulations, the Pc s(-) Sw and anw s(-) Sw main hysteresis loops were compared to experimental data reported by Culligan et al. [Culligan KA, Wildenschild D, Christensen BS, Gray WG, Rivers ML, Tompson AB. Interfacial area measurements for unsaturated flow through porous media. Water Resour Res 2004;40:W12413]. In general, the comparison shows that the simulations are reliable and capture the important physical processes in the experimental system. Pc s(-) Sw curves, anw s(-) Sw curves and phase distributions (within the pores) show good agreement during drainage, but less satisfactory agreement during imbibition. Drainage and imbibition scanning curves were simulated in order to construct Pc s(-) Sw s(-) anw surfaces. The root mean squared error (RMSE) and mean absolute error (MAE) between drainage and imbibition surfaces was 0.10 mm-1 and 0.03 mm-1, respectively. This small difference indicates that hysteresis is virtually nonexistent in the Pc s(-) Sw s(-) anw relationship for the multiphase system studied here. Additionally, a surface was fit to the main loop (excluding scanning curves) of the drainage and imbibition Pc s(-) Sw s(-) anw data and compared to the surface fit to all of the data. The differences between these two surfaces were small (RMSE = 0.05 mm-1 and MAE = 0.01 mm-1) indicating that the Pc s(-) Sw s(-) anw surface is adequately represented without the need for the scanning curve data, which greatly reduces the amount of data required to construct the non-hysteretic Pc s(-) Sw s(-) anw surface for this data.
AB - Hysteresis in the relationship between capillary pressure (Pc), wetting phase saturation (Sw) and nonwetting-wetting interfacial area per volume (anw) is investigated using multiphase lattice-Boltzmann simulations of drainage and imbibition in a glass bead porous system. In order to validate the simulations, the Pc s(-) Sw and anw s(-) Sw main hysteresis loops were compared to experimental data reported by Culligan et al. [Culligan KA, Wildenschild D, Christensen BS, Gray WG, Rivers ML, Tompson AB. Interfacial area measurements for unsaturated flow through porous media. Water Resour Res 2004;40:W12413]. In general, the comparison shows that the simulations are reliable and capture the important physical processes in the experimental system. Pc s(-) Sw curves, anw s(-) Sw curves and phase distributions (within the pores) show good agreement during drainage, but less satisfactory agreement during imbibition. Drainage and imbibition scanning curves were simulated in order to construct Pc s(-) Sw s(-) anw surfaces. The root mean squared error (RMSE) and mean absolute error (MAE) between drainage and imbibition surfaces was 0.10 mm-1 and 0.03 mm-1, respectively. This small difference indicates that hysteresis is virtually nonexistent in the Pc s(-) Sw s(-) anw relationship for the multiphase system studied here. Additionally, a surface was fit to the main loop (excluding scanning curves) of the drainage and imbibition Pc s(-) Sw s(-) anw data and compared to the surface fit to all of the data. The differences between these two surfaces were small (RMSE = 0.05 mm-1 and MAE = 0.01 mm-1) indicating that the Pc s(-) Sw s(-) anw surface is adequately represented without the need for the scanning curve data, which greatly reduces the amount of data required to construct the non-hysteretic Pc s(-) Sw s(-) anw surface for this data.
KW - Capillary pressure
KW - Computed microtomography
KW - Interfacial area
KW - Lattice-Boltzmann
KW - Multiphase flow
KW - Porous media
UR - http://www.scopus.com/inward/record.url?scp=70349948852&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=70349948852&partnerID=8YFLogxK
U2 - 10.1016/j.advwatres.2009.08.009
DO - 10.1016/j.advwatres.2009.08.009
M3 - Article
AN - SCOPUS:70349948852
SN - 0309-1708
VL - 32
SP - 1632
EP - 1640
JO - Advances in Water Resources
JF - Advances in Water Resources
IS - 11
ER -