Title page for ETD etd-04252011-174119

Type of Document Master's Thesis
Author Chukwudozie, Chukwudi Paul
Author's Email Address chuckdii2002@yahoo.com,cchukw1@tigers.lsu.edu
URN etd-04252011-174119
Title Pore-Scale Lattice Boltzmann Simulations of Inertial Flows in Realistic Porous Media: A First Principle Analysis of the Forchheimer Relationship
Degree Master of Science in Petroleum Engineering (M.S.P.E.)
Department Petroleum Engineering
Advisory Committee
Advisor Name Title
Mayank, Tyagi Committee Chair
Karsten, E. Thompson Committee Member
White, D. White Committee Member
  • flow regimes in porous media
  • flow streamlines
  • flow visualization
  • beta factor
  • tortuosity
  • permeability
  • non-darcy coefficient
  • LBM
Date of Defense 2010-11-23
Availability unrestricted
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).


  Filename       Size       Approximate Download Time (Hours:Minutes:Seconds) 
 28.8 Modem   56K Modem   ISDN (64 Kb)   ISDN (128 Kb)   Higher-speed Access 
  finalthesis.pdf 12.51 Mb 00:57:55 00:29:47 00:26:03 00:13:01 00:01:06

Browse All Available ETDs by ( Author | Department )

If you have questions or technical problems, please Contact LSU-ETD Support.