Permeability estimation on tomographic images using curved boundary schemes in the lattice Boltzmann method

dc.citation.articleNumber103685
dc.citation.journalTitleAdvances in Water Resources
dc.citation.volumeNumber143
dc.contributor.authorRao, Parthib
dc.contributor.authorSchaefer, Laura
dc.date.accessioned2020-08-17T19:26:48Z
dc.date.available2020-08-17T19:26:48Z
dc.date.issued2020
dc.description.abstractThe lattice Boltzmann method (LBM) is a widely-used numerical technique for simulation of single- and two-phase flow in geometries that are obtained using tomographic imaging of natural porous media. Due to ease of implementation and numerical robustness, a vast majority of LBM-based pore-scale simulations employ the so-called bounceback scheme to enforce no-slip velocity boundary condition. Bounceback, however, requires an implicit and tight coupling between the numerical (computational) and image (voxel) grid. This coupling results in large discretization errors, since the pore-matrix interface within the 3D image is rough. This leads to overestimation of the interfacial area, and thereby inaccurate permeability predictions. The use of the bounceback scheme also causes other numerical artifacts, such as viscosity-dependent permeability results. In order to address these deficiencies, in this work, the classical marching cubes algorithm is used to reconstruct a surface mesh from the 3D voxel grid; this mesh approximates the pore-matrix surface with higher accuracy compared to the inherent stair-stepped representation. In addition, (nominally) second-order accurate curved boundary schemes are used to enforce no-slip velocity conditions at the reconstructed pore-matrix interface.The various pre-processing steps, such as surface mesh generation and voxelization, that are necessary to use curved boundary schemes are described in detail. The proposed approach of using curved surfaces and boundary schemes is tested and validated on benchmark pore geometries, including a random packing of monodisperse spheres. We conclude that compared to current methods, curved boundary schemes provide a viable option for obtaining more accurate transport properties for Digital Rock Physics-based applications.
dc.identifier.citationRao, Parthib and Schaefer, Laura. "Permeability estimation on tomographic images using curved boundary schemes in the lattice Boltzmann method." <i>Advances in Water Resources,</i> 143, (2020) Elsevier: https://doi.org/10.1016/j.advwatres.2020.103685.
dc.identifier.doihttps://doi.org/10.1016/j.advwatres.2020.103685
dc.identifier.urihttps://hdl.handle.net/1911/109233
dc.language.isoeng
dc.publisherElsevier
dc.rightsThis is an author's peer-reviewed final manuscript, as accepted by the publisher. The published article is copyrighted by Elsevier.
dc.subject.keywordLattice Boltzmann method
dc.subject.keywordDigital rock physics
dc.subject.keywordPermeability
dc.subject.keywordPore-scale flow simulations
dc.subject.keywordBoundary conditions
dc.titlePermeability estimation on tomographic images using curved boundary schemes in the lattice Boltzmann method
dc.typeJournal article
dc.type.dcmiText
dc.type.publicationpost-print
Files
Original bundle
Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
AWR_revised_manuscript.pdf
Size:
4.34 MB
Format:
Adobe Portable Document Format
Description: