Over the last two decades, lattice Boltzmann methods have become an increasingly popular tool to compute the flow in complex geometries such as porous media. In addition to single phase simulations allowing, for example, a precise quantification of the permeability of a porous sample, a number of extensions to the lattice Boltzmann method are available which allow to study multiphase and multicomponent flows on a pore scale level. In this article, we give an extensive overview on a number of these diffuse interface models and discuss their advantages and disadvantages. Furthermore, we shortly report on multiphase flows containing solid particles, as well as implementation details and optimization issues.

During the last decade, lattice-Boltzmann (LB) simulations have been improved to become an efficient tool for determining the permeability of porous media samples. However, well known improvements of the original algorithm are often not implemented. These include for example multirelaxation time schemes or improved boundary conditions, as well as different possibilities to impose a pressure gradient. This paper shows that a significant difference of the calculated permeabilities can be found unless one uses a carefully selected setup. We present a detailed discussion of possible simulation setups and quantitative studies of the influence of simulation parameters. We illustrate our results by applying the algorithm to a Fontainebleau sandstone and by comparing our benchmark studies to other numerical permeability measurements in the literature.

The lattice Boltzmann method has been successfully applied for the simulation of flow through porous media in the creeping regime. Its technical properties, namely discretization, straightforward implementation and parallelization, are responsible for its popularity. However, flow through porous media is not restricted to near zero Reynolds numbers since inertial effects play a role in numerous natural and industrial processes. In this paper we investigate the capability of the lattice Boltzmann method to correctly describe flow in porous media at moderate Reynolds numbers. The selection of the lattice resolution, the collision kernel and the boundary conditions becomes increasingly important and the challenge is to keep artifacts due to compressibility effects at a minimum. The lattice Boltzmann results show an accurate quantitative agreement with Finite Element Method results and evidence the capability of the method to reproduce Darcy's law at low Reynolds numbers and Forchheimer's law at high Reynolds numbers.

Chaotic micromixers such as the staggered herringbone mixer developed by Stroock et al. allow efficient mixing of fluids even at low Reynolds number by repeated stretching and folding of the fluid interfaces. The ability of the fluid to mix well depends on the rate at which ''chaotic advection'' occurs in the mixer. An optimization of mixer geometries is a non-trivial task which is often performed by time consuming and expensive trial and error experiments. In this paper an algorithm is presented that applies the concept of finite-time Lyapunov exponents to obtain a quantitative measure of the chaotic advection of the flow and hence the performance of micromixers. By performing lattice Boltzmann simulations of the flow inside a mixer geometry, introducing massless and noninteracting tracer particles and following their trajectories the finite time Lyapunov exponents can be calculated. The applicability of the method is demonstrated by a comparison of the improved geometrical structure of the staggered herringbone mixer with available literature data.

No abstract