Lattice-Boltzmann simulations of the capillary pressure-saturation-interfacial area relationship for porous media

Mark L. Porter, Marcel Schaap, Dorthe Wildenschild

Research output: Contribution to journalArticle

120 Citations (Scopus)

Abstract

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.

Original languageEnglish (US)
Pages (from-to)1632-1640
Number of pages9
JournalAdvances in Water Resources
Volume32
Issue number11
DOIs
StatePublished - Nov 2009

Fingerprint

capillary pressure
porous medium
imbibition
saturation
drainage
simulation
hysteresis
wetting
unsaturated flow
glass
river

Keywords

  • Capillary pressure
  • Computed microtomography
  • Interfacial area
  • Lattice-Boltzmann
  • Multiphase flow
  • Porous media

ASJC Scopus subject areas

  • Water Science and Technology

Cite this

Lattice-Boltzmann simulations of the capillary pressure-saturation-interfacial area relationship for porous media. / Porter, Mark L.; Schaap, Marcel; Wildenschild, Dorthe.

In: Advances in Water Resources, Vol. 32, No. 11, 11.2009, p. 1632-1640.

Research output: Contribution to journalArticle

@article{beb78993fecf4093aae129e32e1ed341,
title = "Lattice-Boltzmann simulations of the capillary pressure-saturation-interfacial area relationship for porous media",
abstract = "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.",
keywords = "Capillary pressure, Computed microtomography, Interfacial area, Lattice-Boltzmann, Multiphase flow, Porous media",
author = "Porter, {Mark L.} and Marcel Schaap and Dorthe Wildenschild",
year = "2009",
month = "11",
doi = "10.1016/j.advwatres.2009.08.009",
language = "English (US)",
volume = "32",
pages = "1632--1640",
journal = "Advances in Water Resources",
issn = "0309-1708",
publisher = "Elsevier Limited",
number = "11",

}

TY - JOUR

T1 - Lattice-Boltzmann simulations of the capillary pressure-saturation-interfacial area relationship for porous media

AU - Porter, Mark L.

AU - Schaap, Marcel

AU - Wildenschild, Dorthe

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

VL - 32

SP - 1632

EP - 1640

JO - Advances in Water Resources

JF - Advances in Water Resources

SN - 0309-1708

IS - 11

ER -