Generalized equilibria for color-gradient lattice Boltzmann model based on higher-order Hermite polynomials: A simplified implementation with central moments

We propose generalized equilibria of a three-dimensional color-gradient lattice Boltzmann model for two-component two-phase flows using higher-order Hermite polynomials. Although the resulting equilibrium distribution function, which includes a sixth-order term on the velocity, is computationally cumbersome, its equilibrium central moments (CMs) are velocity-independent and have a simplified form. Numerical experiments show that our approach, as in Wen et al. [{Phys. Rev. E \textbf{100}, 023301 (2019)}] who consider terms up to third order, improves the Galilean invariance compared to that of the conventional approach. Dynamic problems can be solved with high accuracy at a density ratio of 10; however, the accuracy is still limited to a density ratio of 1000. For lower density ratios, the generalized equilibria benefit from the CM-based multiple-relaxation-time model, especially at very high Reynolds numbers, significantly improving the numerical stability.Comment: 22 pages, 8 figure

A three-dimensional non-orthogonal multiple-relaxation-time phase-field lattice Boltzmann model for multiphase flows at large density ratios and high Reynolds numbers

This study proposes a three-dimensional non-orthogonal multiple-relaxation-time (NMRT) phase-field multiphase lattice Boltzmann (PFLB) model within a recently established unified lattice Boltzmann model (ULBM) framework [Luo et al., Phil. Trans. R. Soc. A 379, 20200397, 2021]. The conservative Allen-Cahn equation and the incompressible Navier-Stokes (NS) equations are solved. In addition, a local gradient calculation scheme for the order parameter of the Allen-Cahn equation is constructed with the non-equilibrium part of the distribution function. A series of benchmark cases are conducted to validate the proposed model, including the two-phase Poiseuille flow, Rayleigh-Taylor instability, binary liquid/metal droplet collision, and a bubble rise in water. The present simulation results are in good agreement with existing simulation and experimental data. In the simulation of the co-current two-phase Poiseuille flow, the present model is proven to resolve the discontinuity at the phase interface and provide accurate results at extremely high density ratios (i.e., up to ). Finally, the proposed model is adopted to simulate two challenging cases: (1) water droplet splashing during its impacting on a thin liquid film and (2) liquid jet breakup. The simulation results demonstrate an excellent agreement with previous experimental results, both qualitatively and quantitatively. In these simulations, the Weber number and Reynolds number reach 105 and 6000, respectively, and the viscosity can be as low as , in the lattice unit

Lattice Boltzmann methods for multiphase flow and phase-change heat transfer

Over the past few decades, tremendous progress has been made in the development of particle-based discrete simulation methods versus the conventional continuum-based methods. In particular, the lattice Boltzmann (LB) method has evolved from a theoretical novelty to a ubiquitous, versatile and powerful computational methodology for both fundamental research and engineering applications. It is a kinetic-based mesoscopic approach that bridges the microscales and macroscales, which offers distinctive advantages in simulation fidelity and computational efficiency. Applications of the LB method are now found in a wide range of disciplines including physics, chemistry, materials, biomedicine and various branches of engineering. The present work provides a comprehensive review of the LB method for thermofluids and energy applications, focusing on multiphase flows, thermal flows and thermal multiphase flows with phase change. The review first covers the theoretical framework of the LB method, revealing certain inconsistencies and defects as well as common features of multiphase and thermal LB models. Recent developments in improving the thermodynamic and hydrodynamic consistency, reducing spurious currents, enhancing the numerical stability, etc., are highlighted. These efforts have put the LB method on a firmer theoretical foundation with enhanced LB models that can achieve larger liquid-gas density ratio, higher Reynolds number and flexible surface tension. Examples of applications are provided in fuel cells and batteries, droplet collision, boiling heat transfer and evaporation, and energy storage. Finally, further developments and future prospect of the LB method are outlined for thermofluids and energy applications

A unified lattice Boltzmann model and application to multiphase flows

In this work, we develop a unified lattice Boltzmann model (ULBM) framework that can seamlessly integrate the widely used lattice Boltzmann collision operators, including the Bhatnagar–Gross–Krook or single-relation-time, multiple-relaxation-time, central-moment or cascaded lattice Boltzmann method and multiple entropic operators (KBC). Such a framework clarifies the relations among the existing collision operators and greatly facilitates model comparison and development as well as coding. Importantly, any LB model or treatment constructed for a specific collision operator could be easily adopted by other operators. We demonstrate the flexibility and power of the ULBM framework through three multiphase flow problems: the rheology of an emulsion, splashing of a droplet on a liquid film and dynamics of pool boiling. Further exploration of ULBM for a wide variety of phenomena would be both realistic and beneficial, making the LBM more accessible to non-specialists. This article is part of the theme issue ‘Progress in mesoscale methods for fluid dynamics simulation’

A Volume of Fluid Method for Three Dimensional Direct Numerical Simulations of Immiscible Droplet Collisions

An advanced Volume of Fluid (VOF) method is presented that enables performant three-dimensional Direct Numerical Simulations (DNS) of the interaction of two immiscible fluids in a gaseous environment with large topology changes, e.g., binary droplet collisions. One of the challenges associated with the introduction of a third immiscible phase into the VOF method is the reconstruction of the phase boundaries near the triple line in arbitrary arrangements. For this purpose, an efficient method based on a Piecewise Linear Interface Calculation (PLIC) is shown. Moreover, the surface force modeling with the robust Continuous Surface Stress (CSS) model was enhanced to treat such three-phase situations with large topology changes and thin films. A consistent scaling of the fluid properties at the interfaces ensures energy conservation. The implementation of these methods in the multi-phase flow solver Free Surface 3D (FS3D) allowed a successful validation. A qualitative comparison of the morphology in binary collisions of immiscible droplets as well as a quantitative comparison regarding the threshold velocities that distinguish different collision regimes shows excellent agreement with experimental results. These simulations enable the evaluation of experimentally inaccessible data like the contributions of the kinetic, surface and dissipative energy of both immiscible liquids during the collision process. Furthermore, the comparison with binary collisions of the same liquids highlights similarities and differences between the collisions. Both can support the modeling of the immiscible liquid interaction in the future.Comment: dataset on https://doi.org/10.18419/darus-3557; Changes in v2: VOF in title written out as Volume of Fluid, a handful of typos removed and some sentences language polished, keywords added, new citation styl

Modeling realistic multiphase flows using a non-orthogonal multiple-relaxation-time lattice Boltzmann method

In this paper, we develop a three-dimensional multiple-relaxation-time lattice Boltzmann method (MRT-LBM) based on a set of non-orthogonal basis vectors. Compared with the classical MRT-LBM based on a set of orthogonal basis vectors, the present non-orthogonal MRT-LBM simplifies the transformation between the discrete velocity space and the moment space, and exhibits better portability across different lattices. The proposed method is then extended to multiphase flows at large density ratio with tunable surface tension, and its numerical stability and accuracy are well demonstrated by some benchmark cases. Using the proposed method, a practical case of a fuel droplet impacting on a dry surface at high Reynolds and Weber numbers is simulated and the evolution of the spreading film diameter agrees well with the experimental data. Furthermore, another realistic case of a droplet impacting on a super-hydrophobic wall with a cylindrical obstacle is reproduced, which confirms the experimental finding of Liu \textit{et al.} [``Symmetry breaking in drop bouncing on curved surfaces," Nature communications 6, 10034 (2015)] that the contact time is minimized when the cylinder radius is comparable with the droplet cylinder.Comment: 19 pages, 11 figure