Date:May 14, 2012

Lattice Boltzmann methods

(snapshots of bubble dynamics during the decompression of a suspension)

The study of complex physical processes involving multi-phase flows, pore-scale processes in porous media, reactive flow problems and diffusion in complex geometries requires the development of novel numerical models. Although we also use finite difference, finite volume or boundary integral methods, our main focus in term of numerical method is the lattice Boltzmann method (LB).

The lattice Boltzmann method is a kinetic theory based approach to solve partial differential equations. It is very successful for fluid dynamics in complex geometries, multiphase flows (especially for immiscible fluids) and flows with evolving solid-fluid boundaries. We have been developing new algorithms for

  • imposing Dirichlet boundary conditions for the diffusion equation in complex geometries (Huber et al., IJMPC, 2011).
  • multiphase flows, free surface flows, bubble dynamics (in prep).
  • melting, dissolution processes coupled with fluid flow (Huber et al., 2008, Parmigiani et al., 2011).
  • pore-scale reactive transport models (Huber and Shafei, submitted).
  • multiphysics (Parmigiani et al., 2008).
  • noble gas diffusion for themochronology (Huber et al., GCA, 2011).
  • multiphase reactive flows in porous media (Parmigiani et al., JFM, 2011).
  • multicomponent coupled diffusion problems (Huber et al., JCP, 2010).


We are proposing a class for graduate student onĀ  lattice Boltzmann methods for geosciences (EAS 8013).