1,234 research outputs found
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
Simulation of two-phase flows at large density ratios and high Reynolds numbers using a discrete unified gas kinetic scheme
In order to treat immiscible two-phase flows at large density ratios and high
Reynolds numbers, a three-dimensional code based on the discrete unified gas
kinetic scheme (DUGKS) is developed, incorporating two major improvements.
First, the particle distribution functions at cell interfaces are reconstructed
using a weighted essentially non-oscillatory scheme. Second, the conservative
lower-order Allen-Cahn equation is chosen, instead of the higher-order
Cahn-Hilliard equation, to evolve the free-energy based phase field governing
the dynamics of two-phase interfaces. Five benchmark problems are simulated to
demonstrate the capability of the approach in treating two phase flows at large
density ratios and high Reynolds numbers, including three two dimensional
problems (a stationary droplet, Rayleigh-Taylor instability, and a droplet
splashing on a thin liquid film) and two three-dimensional problems (binary
droplets collision and Rayleigh-Taylor instability). All results agree well
with the previous numerical and the experimental results. In these simulations,
the density ratio and Reynolds number can reach a large value of O(1000). Our
improved approach sets the stage for the DUGKS scheme to handle realistic
two-phase flow problems
Development of lattice boltzmann flux solvers and their applications
Ph.DDOCTOR OF PHILOSOPH
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
- …