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

Linlin Fei,Jingyu Du,Kai H. Luo,Sauro Succi,Marco Lauricella,Andrea Montessori,Qian Wang

DOI: https://doi.org/10.1063/1.5087266

2023-06-28

Abstract: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.

Fluid Dynamics,Computational Physics

What problem does this paper attempt to address?

This paper attempts to solve some key problems in multiphase flows, especially the simulation of complex interface phenomena under high - density ratio and adjustable surface tension conditions. Specifically, the main objectives of the article include:

1. **Develop an improved computational method**: The author proposes a three - dimensional non - orthogonal multiple - relaxation - time lattice Boltzmann method (MRT - LBM). Compared with the traditional method based on orthogonal basis vectors, this method simplifies the conversion between discrete velocity space and moment space and has better portability on different lattices.

2. **Improve numerical stability and accuracy**: By introducing the non - orthogonal MRT - LBM, the researchers aim to improve numerical stability, especially when dealing with practical multiphase flow problems with extremely high Reynolds numbers and Weber numbers. The stability and accuracy of this method are verified through several benchmark cases in the paper.

3. **Simulate complex phenomena in practical applications**: To demonstrate the capabilities of the new method, the author simulates two challenging practical problems:
- The process of a fuel droplet hitting a dry surface with high Reynolds and Weber numbers, and the spreading film diameter agrees well with the experimental data.
- A droplet hitting a super - hydrophobic wall with a cylindrical obstacle, reproducing the phenomenon of droplet bouncing symmetry breaking on a curved surface by Liu et al. (2015), that is, when the radius of the cylinder is comparable to that of the droplet, the contact time is minimized.

4. **Expand to other fields**: In addition to the above - mentioned specific applications, this method can also be applied to a wider range of multiphase flow problems, such as the complex interface phenomena involved in natural and industrial processes such as raindrop dynamics, plant spraying, water boiling, and gas turbine blade cooling.

In summary, the core issue of this paper is to develop a more efficient and stable computational method to accurately capture and simulate the complex multiphase flow interface behavior, especially under high - density ratio and extreme flow conditions. This not only helps to deeply understand the physical mechanisms behind these phenomena but also provides strong tool support for engineering design in related fields.