Return to search

Pore-Scale Lattice Boltzmann Simulations of Inertial Flows in Realistic Porous Media: A First Principle Analysis of the Forchheimer Relationship

With recent advances in the capabilities of high performance computing (HPC) platforms and the relatively simple representation of complex geometries of porous media, lattice Boltzmann method (LBM) has gained popularity as a means of solving fluid flow and transport problems. In this work, LBM was used to obtain flow parameters of porous media, study the behavior of these parameters at varying flow conditions and quantify the effect of roughness on the parameters by relating the volume averaged flow simulation results to Darcy and Forchheimer equations respectively.
To validate the method, flow was simulated on regular and random sphere arrays in cubic domains, for which a number of analytical solutions are available. Permeability and non-Darcy coefficients obtained from the simulation compared well with Kozeny and Ergun estimates while deviation from the observed constant permeability and tortuosity values occurred aroundRe≈1-10. By defining roughness as hemispherical protrusions on the smooth spheres in the regular array, it was observed from flow streamlines obtained at different roughness heights that the average length of the flow paths increased with increasing roughness height. As such, the medium tortuosity and non-Darcy coefficient increased while the permeability decreased as height of the roughness increased.
Applying the method to a 3D computed tomography image of Castlegate sandstone, the calculated macroscopic permeability and beta factor components were in good agreement with reported experimental values. In addition, LBM beta factors were compared with a number of empirical models for non-Darcy coefficient estimation and were found to be of the same order of magnitude as most of the correlations, although estimates of the models showed wide variation in values. Resolution of the original sample was increased by infilling with more voxels and simulation in the new domain showed better flow field resolution and higher simulated flow regimes compared to those of the original sample, without significant change in the flow parameters obtained. Using the Reynolds number based on the Forchheimer coefficient, the range of transition from Darcy to non-Darcy regime was within the values reported by Ruth and Ma (1993) and Zeng and Grigg (2006).
1

Identiferoai:union.ndltd.org:LSU/oai:etd.lsu.edu:etd-04252011-174119
Date27 April 2011
CreatorsChukwudozie, Chukwudi Paul
ContributorsKarsten, E. Thompson, Mayank, Tyagi, White, D. White
PublisherLSU
Source SetsLouisiana State University
LanguageEnglish
Detected LanguageEnglish
Typetext
Formatapplication/pdf
Sourcehttp://etd.lsu.edu/docs/available/etd-04252011-174119/
Rightsunrestricted, I hereby certify that, if appropriate, I have obtained and attached herein a written permission statement from the owner(s) of each third party copyrighted matter to be included in my thesis, dissertation, or project report, allowing distribution as specified below. I certify that the version I submitted is the same as that approved by my advisory committee. I hereby grant to LSU or its agents the non-exclusive license to archive and make accessible, under the conditions specified below and in appropriate University policies, my thesis, dissertation, or project report in whole or in part in all forms of media, now or hereafter known. I retain all other ownership rights to the copyright of the thesis, dissertation or project report. I also retain the right to use in future works (such as articles or books) all or part of this thesis, dissertation, or project report.

Page generated in 0.0156 seconds