Lattice Boltzmann equation method for multiple immiscible continuum fluids

SPENCER, Timothy, HALLIDAY, Ian and CARE, Chris (2010). Lattice Boltzmann equation method for multiple immiscible continuum fluids. Physical Review E, 82 (066701), 1-22.

Full text not available from this repository.
Official URL:
Link to published version::


This paper generalizes the two-component algorithm of Sec. II, extending it, in Sec. III, to describe N>2 mutually immiscible fluids in the isothermal continuum regime. Each fluid has an independent interfacial tension. While retaining all its computational advantages, we remove entirely the empiricism associated with contact behavior in our previous multiple immiscible fluid models [M. M. Dupin et al., Phys. Rev. E 73, 055701(R) (2006); Med. Eng. Phys. 28, 13 (2006)] while solidifying the physical foundations. Moreover, the model relies upon a fluid-fluid segregation which is simpler, computationally faster, more free of artifacts i.e., the interfacial microcurrent, and upon an interface-inducing force distribution which is analytic. The method is completely symmetric between any numbers of immiscible fluids and stable over a wide range of directly input interfacial tension. We present data on the steady-state properties of multiple interface model, which are in good agreement with theory [R. E. Johnson and S. S. Sadhal, Annu. Rev. Fluid Mech. 17, 289 (1985)], specifically on the shapes of multidrop systems. Section II is an analysis of the kinetic and continuum-scale descriptions of the underlying two-component lattice Boltzmann model for immiscible fluids, extendable to more than two immiscible fluids. This extension requires (i) the use of a more local kinetic equation perturbation which is (ii) free from a reliance on measured interfacial curvature. It should be noted that viewed simply as a two-component method, the continuum algorithm is inferior to our previous methods, reported by Lishchuk et al. [Phys. Rev. E 67, 036701 (2003)] and Halliday et al. [Phys. Rev. E 76, 026708 (2007)]. Greater stability and parameter range is achieved in multiple drop simulations by using the forced multirelaxation-time lattice Boltzmann method developed, along with (for completeness) a forced exactly incompressible Bhatnagar-Gross-Krook lattice Boltzmann model, in the Appendix. These appended schemes closely follow those developed by Guo et al. [Phys. Rev. E 65, 046308 (2002)] for the single-relaxation-time scheme.

Item Type: Article
Research Institute, Centre or Group - Does NOT include content added after October 2018: Materials and Engineering Research Institute > Modelling Research Centre > Materials Modelling group
Identification Number:
Page Range: 1-22
Depositing User: Timothy Spencer
Date Deposited: 11 Apr 2011 13:55
Last Modified: 18 Mar 2021 21:00

Actions (login required)

View Item View Item


Downloads per month over past year

View more statistics