Axisymmetric multiphase lattice Boltzmann method

A lattice Boltzmann method for axisymmetric multiphase flows is presented and validated. The method is capable of accurately modeling flows with variable density. We develop the classic Shan-Chen multiphase model [ Phys. Rev. E 47 1815 (1993)] for axisymmetric flows. The model can be used to efficiently simulate single and multiphase flows. The convergence to the axisymmetric Navier-Stokes equations is demonstrated analytically by means of a Chapmann-Enskog expansion and numerically through several test cases. In particular, the model is benchmarked for its accuracy in reproducing the dynamics of the oscillations of an axially symmetric droplet and on the capillary breakup of a viscous liquid thread. Very good quantitative agreement between the numerical solutions and the analytical results is observed

Lattice Boltzmann Model for Axisymmetric Multiphase Flows

In this paper, a lattice Boltzmann (LB) model is presented for axisymmetric multiphase flows. Source terms are added to a two-dimensional standard lattice Boltzmann equation (LBE) for multiphase flows such that the emergent dynamics can be transformed into the axisymmetric cylindrical coordinate system. The source terms are temporally and spatially dependent and represent the axisymmetric contribution of the order parameter of fluid phases and inertial, viscous and surface tension forces. A model which is effectively explicit and second order is obtained. This is achieved by taking into account the discrete lattice effects in the Chapman-Enskog multiscale analysis, so that the macroscopic axisymmetric mass and momentum equations for multiphase flows are recovered self-consistently. The model is extended to incorporate reduced compressibility effects. Axisymmetric equilibrium drop formation and oscillations, breakup and formation of satellite droplets from viscous liquid cylindrical jets through Rayleigh capillary instability and drop collisions are presented. Comparisons of the computed results with available data show satisfactory agreement.Comment: 17 pages, 11 figures, to be published in Physical Review

Generalized Lattice Boltzmann Method with multi-range pseudo-potential

The physical behaviour of a class of mesoscopic models for multiphase flows is analyzed in details near interfaces. In particular, an extended pseudo-potential method is developed, which permits to tune the equation of state and surface tension independently of each other. The spurious velocity contributions of this extended model are shown to vanish in the limit of high grid refinement and/or high order isotropy. Higher order schemes to implement self-consistent forcings are rigorously computed for 2d and 3d models. The extended scenario developed in this work clarifies the theoretical foundations of the Shan-Chen methodology for the lattice Boltzmann method and enhances its applicability and flexibility to the simulation of multiphase flows to density ratios up to O(100)

Comparison of multiphase SPH and LBM approaches for the simulation of intermittent flows

Smoothed Particle Hydrodynamics (SPH) and Lattice Boltzmann Method (LBM) are increasingly popular and attractive methods that propose efficient multiphase formulations, each one with its own strengths and weaknesses. In this context, when it comes to study a given multi-fluid problem, it is helpful to rely on a quantitative comparison to decide which approach should be used and in which context. In particular, the simulation of intermittent two-phase flows in pipes such as slug flows is a complex problem involving moving and intersecting interfaces for which both SPH and LBM could be considered. It is a problem of interest in petroleum applications since the formation of slug flows that can occur in submarine pipelines connecting the wells to the production facility can cause undesired behaviors with hazardous consequences. In this work, we compare SPH and LBM multiphase formulations where surface tension effects are modeled respectively using the continuum surface force and the color gradient approaches on a collection of standard test cases, and on the simulation of intermittent flows in 2D. This paper aims to highlight the contributions and limitations of SPH and LBM when applied to these problems. First, we compare our implementations on static bubble problems with different density and viscosity ratios. Then, we focus on gravity driven simulations of slug flows in pipes for several Reynolds numbers. Finally, we conclude with simulations of slug flows with inlet/outlet boundary conditions. According to the results presented in this study, we confirm that the SPH approach is more robust and versatile whereas the LBM formulation is more accurate and faster