An adaptive mesh refinement‐multiphase lattice Boltzmann flux solver for three‐dimensional simulation of droplet collision

Yuan, Haizhuan; Liu, Qing; Zeng, Gang

doi:10.1002/fld.5061

Cited by 5 publications

(3 citation statements)

References 52 publications

(96 reference statements)

Supporting

Mentioning

Contrasting

Order By: Relevance

“…One approach was to consolidate adaptive mesh technique which can refine the local mesh sizes in interfacial regions. 32,33 In this way, the interfacial tracking algorithm, or the Cahn-Hilliard equation itself, does not need to be altered. The other way is to introduce terms that can compensate loss of mass during the computation, which requires modification of the interfacial tracking models.…”

Section: Introductionmentioning

confidence: 99%

See 1 more Smart Citation

Multiphase lattice Boltzmann flux solver with conservative Allen‐Cahn model for modeling high‐density‐ratio flows

Chen,

Sun

2024

Numerical Methods in Fluids

View full text Add to dashboard Cite

In this paper, the Allen‐Cahn‐Multiphase lattice Boltzmann flux solver (AC‐MLBFS) is proposed as a new and effective numerical simulation method for multiphase flows with high density ratios. The MLBFS resolves the macroscopic governing equations with the finite volume method and reconstructs numerical fluxes on the cell interface from local solutions to the lattice Boltzmann equation, which combines the advantages of conventional Navier–Stokes solvers and lattice Boltzmann methods for simulating incompressible multiphase flows while alleviating their limitations. Previous MLBFS‐based multiphase solvers performed poorly in mass conservation, which might be caused by the excessive numerical diffusion in the Cahn‐Hilliard (CH) model used as the interface tracking algorithm. To resolve this problem, the present method proposes using the conservative Allen‐Cahn (AC) model as the interfacial tracking algorithm, which can ease the numerical implementation by removing high order derivative terms and alleviate mass leakage by enforcing local mass conservation in the physical model. Numerical validations will be carried out through benchmark tests at high density ratios and in extreme conditions with large Reynolds or Weber numbers. Through these examples, the accuracy and robustness as well as the mass conservation characteristics of the proposed method are demonstrated.

show abstract

Section: Introductionmentioning

confidence: 99%

“…To alleviate this limitation, two types of strategies have been proposed. One approach was to consolidate adaptive mesh technique which can refine the local mesh sizes in interfacial regions 32,33 . In this way, the interfacial tracking algorithm, or the Cahn‐Hilliard equation itself, does not need to be altered.…”

Section: Introductionmentioning

confidence: 99%

Multiphase lattice Boltzmann flux solver with conservative Allen‐Cahn model for modeling high‐density‐ratio flows

Chen,

Sun

2024

Numerical Methods in Fluids

View full text Add to dashboard Cite

show abstract

“…Existing methods to approximate the distribution function for non-equilibrium flows stem from two ideas: the Chapman-Enskog (CE) expansion [6] and the Hermite polynomial expansion [7]. Specifically, the CE expansion performs power series expansion on the equilibrium distribution function [8][9][10]. With more expansion terms, the constructed distribution function is expected to more precisely describe the deviation from equilibrium state.…”

Section: Introductionmentioning

confidence: 99%

Construction of Non-Equilibrium Gas Distribution Functions Through Expansions in Peculiar Velocity Space

Yuan¹,

Chen²,

Shu³

et al. 2023

AAMM

View full text Add to dashboard Cite

Gas distribution function plays a crucial role in the description of gas flows at the mesoscopic scale. In the presence of non-equilibrium flow, the distribution function loses its rotational symmetricity, making the mathematical derivation difficult. From both the Chapman-Enskog expansion and the Hermite polynomial expansion (Grad's method), we observe that the non-equilibrium effect is closely related to the peculiar velocity space (C). Based on this recognition, we propose a new methodology to construct the non-equilibrium distribution function from the perspective of polynomial expansion in the peculiar velocity space of molecules. The coefficients involved in the non-equilibrium distribution function can be exactly determined by the compatibility conditions and the moment relationships. This new framework allows constructing non-equilibrium distribution functions at any order of truncation, and the ones at the third and the fourth order have been presented in this paper for illustration purposes. Numerical validations demonstrate that the new method is more accurate than the Grad's method at the same truncation error for describing non-equilibrium effects. Two-dimensional benchmark tests are performed to shed light on the applicability of the new method to practical engineering problems.

show abstract

Direct numerical simulation of coalescence and separation of binary droplet collision with lamella stabilization

Yu,

Chang

2024

Physics of Fluids

View full text Add to dashboard Cite

Binary droplet collision is a fundamental aspect of various natural phenomena and industrial applications. In this work, direct numerical simulation of coalescence and separation of binary droplet collision is performed over a wide range of Weber numbers and impact factors. The incompressible Navier–Stokes equations are solved by the finite volume approach, coupled with the volume of fluid method. To address the inaccurate prediction of thin lamella in simulation, a lamella stabilization method is introduced to resolve the lamella by adjusting the grid resolution. Compared with experimental data, it is validated that the lamella can be accurately and fully captured with this lamella stabilization method. Moreover, the analysis of shape and energy during the collision is conducted, and the variation of lamella is described in detail, particularly the evolution of the thickness of lamella. The results suggest that for obtaining the full variation of lamella, the maximum refinement size of the grid can reach D/4096. It is also found that without lamella stabilization, excessive dissipation can lead to the failure of predicting coalescence and separation, especially for the cases in the transition between coalescence and separation. Furthermore, even if the same collision outcome can be obtained without lamella stabilization, the number and size of droplets have obvious differences.

show abstract

An adaptive mesh refinement‐multiphase lattice Boltzmann flux solver for three‐dimensional simulation of droplet collision

Cited by 5 publications

References 52 publications

Multiphase lattice Boltzmann flux solver with conservative Allen‐Cahn model for modeling high‐density‐ratio flows

Multiphase lattice Boltzmann flux solver with conservative Allen‐Cahn model for modeling high‐density‐ratio flows

Construction of Non-Equilibrium Gas Distribution Functions Through Expansions in Peculiar Velocity Space

Direct numerical simulation of coalescence and separation of binary droplet collision with lamella stabilization

Contact Info

Product

Resources

About