SUBJECT: Ph.D. Dissertation Defense
   
BY: Irfan Khan
   
TIME: Tuesday, June 22, 2010, 2:30 p.m.
   
PLACE: MRDC Building, 4211
   
TITLE: Direct Numerical Simulation and Analysis of Deformation in Saturated Porous Media
   
COMMITTEE: Dr. Cyrus K. Aidun, Chair (ME)
Dr. Mostafa Ghiaasiaan (ME)
Dr. Suresh Sitaraman (ME)
Dr. George Biros (CC)
Dr. Thorsten Stoesser (CE)
 

SUMMARY

Applications of fluid infiltrated deformable porous media are found in soil consolidation, filtration and absorbency products, fabric and textiles, human tissue and bone. In such media, the geometrical complexity of the solid structure coupled with the need for a unified treatment of all the interacting phases makes a micro-mechanical approach of the problem intractable. The the existing techniques, such as “Biot’s theory” or the “Theory of Porous Media” are based on homogenization techniques, in the sense that they assume the solid and all the fluid phases as a smeared media. These kind of approaches, inspite of their mathematical rigour, are not suitable for micro-mechanical investigations such as the effect of micro-structure on the deformational behaviour or the constitutive relations in the media during deformation. In order to perform such investigations an approach based on direct numerical simulations, capable of leveraging the ever increasing computational power would be ideally suited. In this research work, a numerical scheme based on the hybrid lattice-Boltzmann finite-element method is developed to carry out the direct numerical simulations of deformable porous media. The method has been parallelized to make use of distributed computing. The efficiency and inherent parallel nature of the lattice Boltzmann method coupled with the robustness of finite element method implemented in a parallel framework provided by the highly scalable toolkit of PETSc lend as a powerful tool to not only tackle the problem of deformable porous media but also to undertake any problem involving complex interaction of fluid and solid phases in the linear elastic limit of deformation. The method has been used to understand the deformational characteristics of model porous media made up of spheres and cylinders. The deformational behaviour is seen to match with the existing analytical solution closely. Thus it is found that macroscopic behaviour of a generic porous media can be recovered with the use of model porous media constructed with simplified geometries. This finding motivated research in using model porous geometries to represent more complex real porous geometries in order to perform investigations of deformation on the latter. An attempt has been made to apply this technique to the complex geometries of “felt”, (a fibrous mat used in paper industries). These investigations lead to new understanding on the effect of fiber diameter on the bulk properties of a fibrous media and subsequently on the deformational behaviour of the media. Further the method has been used to investigate the constitutive relationships in deformable porous media. Particularly the relationship between permeability and porosity during the deformation of the media is investigated. Results show the need of geometry specific investigations.