Direct numerical simulation of saturated deformable porous media using parallel hybrid Lattice-Boltzmann and finite element method

Numerical techniques for modeling saturated deformable porous media have mainly been based on mixture theory or homogenization techniques. However, these techniques rely on phenomenological relationships for the constitutive equations along with assumptions of homogeneous and isotopic material properties to obtain closure. Direct numerical simulations of the multiphasic problem for flow in deformable porous media avoid such assumptions and thus can provide significantly accurate understanding of the physics involved. They serve as a tool to investigate the constitutive relationships in complex geometries. They also allow the validation of the existing mixture theory models and determine their limitations. In this work, a parallel hybrid method using Lattice Boltzmann Method (LBM) for fluid phase and Finite Element Method (FEM) for solid phase is used for direct numerical simulation of saturated deformable porous media. The method provides a number of unique features including scalability on distributed computing necessary for such a problem. The method has been validated for modeling fluid-structure interactions in complex geometries against a number of experimental and analytical solutions. Further some challenging problems has been chosen to show the capability of the method. Copyright (C) 2010 John Wiley & Sons, Ltd.