Dec 12 – 13, 2016
Institut Henri Poincaré
Europe/Paris timezone

Lagrange-Euler Lattice-Boltzmann Method And Its Application to Two-Fluid Flows Dynamics With Possibly High Density Ratio

Dec 13, 2016, 11:45 AM
Amphi Hermite (Institut Henri Poincaré)

Amphi Hermite

Institut Henri Poincaré

11 Rue Pierre et Marie Curie, 75005 Paris, France


Mrs Marie Béchereau (Ecole normale supérieure de Cachan (Paris-Saclay))


Two-fluid extensions of Lattice Boltzmann methods with free boundaries usually consider "microscopic'' pseudopotential interface models. In this paper, we rather propose an interface-capturing Lattice Boltzmann approach where the mass fraction variable is considered as an unknown and is advected. Several works have reported the difficulties of LBM methods to deal with such two-fluid systems especially for high-density ratio configurations. This is due to the mixing nature of LBM, as with Flux vector splitting approaches for Finite Volume methods. We here give another explanation of the lack of numerical diffusion of Lattice Boltzmann approaches to accurately capture contact discontinuities. To fix the problem, we propose an arbitrary Lagrangian-Eulerian (ALE) formulation of Lattice-Boltzmann methods. In the Lagrangian limit, it allows for a proper separated treatment of pressure waves and advection phenomenon. After the ALE solution, a remapping (advection) procedure is necessary to project the variables onto the Eulerian Lattice-Boltzmann grid. We explain how to derive this remapping procedure in order to get second-order accuracy and achieve sharp stable oscillation-free interfaces. It has been shown that mass fractions variables satisfy a local discrete maximum principle and thus stay in the range $[0,1]$. The theory is supported by numerical computations of the free fall of an initial square block of a dense fluid surrounded by a lighter fluid into a box. Figures 1 and 2 are showing the mass fraction field of the light fluid at two successive instants. The density ratio equal to 4 and the computational lattice grid is 400x400. One can appreciate the thickness of the numerical diffuse interface, the capture of complex structures and the capability to compute strong changes of free boundary topology. Even if our methods are currently used for inviscid flows (Euler equations) by projecting the discrete distributions onto equilibrium ones at each time step, we believe that it is possible to extend the framework formulation for multifluid viscous problems. This will be at the aim of a next work.

Presentation materials