Search

In this work, we present a general second-order phase-field model for the transport of insoluble surfactant in incompressible two-phase flows. In this model, the second-order local Allen-Cahn equation is applied for interface capturing, a general form of the simple scalar transport equation [S. S. Jain, J. Comput. Phys. 515, 113277 (2024)] is adopted for interface-confined surfactant, and the consistent and conservative Navier-Stokes equations with the Marangoni force is used for fluid flows. To solve this model, we further developed a mesoscopic lattice Boltzmann (LB) method, in which the LB model for surfactant transport equation is proposed under the general LB framework for the convection-diffusion type equation, and it can correctly recover the governing equation for surfactant transport. The accuracy of the present LB method is tested by several benchmark problems, and the numerical results show it has a good performance for the transport of the insoluble surfactant in two-phase flows.

In this work, a phase-field model is developed for the dendritic growth with gas bubbles in the solidification of binary alloys. In this model, a total free energy for the complex gas-liquid-dendrite system is proposed through considering the interactions of gas bubbles, liquid melt and solid dendrites, and it can reduce to the energy for gas-liquid flows in the region far from the solid phase, while degenerate to the energy for thermosolutal dendritic growth when the gas bubble disappears. The governing equations are usually obtained by minimizing the total free energy, but here some modifications are made to improve the capacity of the conservative phase-field equation for gas bubbles and convection-diffusion equation for solute transfer. Additionally, through the asymptotic analysis of the thin-interface limit, the present general phase-field model for alloy solidification can match the corresponding free boundary problem, and it is identical to the commonly used models under a specific choice of model parameters. Furthermore, to describe the fluid flow, the incompressible Navier-Stokes equations are adopted in the entire domain including gas, liquid, and solid regions, where the fluid-structure interaction is considered by a simple diffuse-interface method. To test the present phase-field model, the lattice Boltzmann method is used to study several problems of gas-liquid flows, dendritic growth as well as the solidification in presence of gas bubbles, and a good performance of the present model for such complex problems is observed., Comment: 29 pages, 23 figures

In this work, we develop a phase-field-based lattice Boltzmann (LB) method for a two-scalar model of the two-phase flows with interfacial mass/heat transfer. Through the Chapman-Enskog analysis, we show that the present LB method can correctly recover the governing equations for phase field, flow field and concentration/temperature field. In particular, to derive the two-scalar equations for the mass/heat transfer, we propose a new LB model with an auxiliary source distribution function to describe the extra flux terms, and the discretizations of some derivative terms can be avoided. The accuracy and efficiency of the present method is also tested through several benchmark problems, and the influence of mass/heat transfer on the fluid viscosity is further considered by introducing an exponential relation. The numerical results show that the present LB method is suitable for the two-phase flows with interfacial mass/heat transfer., Comment: 30 pages, 12 figures

Modeling and simulation of multiphase flows in complex geomerties are challenging due to the complexity in describing the interface topology changes among different phases and the difficulty in implementing the boundary conditions on the irregular solid surface. In this work, we first developed a diffuse-domain (DD) based phase-field model for multiphase flows in complex geometries. In this model, the irregular fluid region is embedded into a larger and regular domain by introducing a smooth characteristic function. Then, the reduction-consistent and conservative phase-field equation for the multiphase field and the consistent and conservative Navier-Stokes equations for the flow field are reformulated as the diffuse-domain based consistent and conservative (DD-CC) equations where some additional source terms are added to reflect the effects of boundary conditions. In this case, there is no need to directly treat the complex boundary conditions on the irregular solid surface, and additionally, based on a matched asymptotic analysis, it is also shown that the DD-CC equations can converge to the original governing equations as the interface width parameter tends to zero. Furthermore, to solve the DD-CC equations, we proposed a novel and simple lattice Boltzmann (LB) method with a Hermite-moment-based collision matrix which can not only keep consistent and conservation properties, but also improve the numerical stability with a flexible parameter. With the help of the direct Taylor expansion, the macroscopic DD-CC equations can be recovered correctly from the present LB method. Finally, to test the capacity of LB method, several benchmarks and complex problems are considered, and the numerical results show that the present LB method is accurate and efficient for the multiphase flows in complex geomerties., Comment: 22 pages, 9 figures

In this work, a thermodynamically consistent and conservative diffuse-interface model for gas-liquid-solid multiphase flows is proposed. In this model, a novel free energy for the gas-liquid-solid multiphase flows is established according to a ternary phase-field model, and it not only contains the standard bulk and interface free energies for two-phase flows, but also includes some additional terms to reflect the penalty in the solid phase and the wettability on the solid surface. Furthermore, a smooth indicator function of the solid phase is also introduced in the consistent Navier-Stokes equations to achieve a high viscosity in the solid phase and preserve the velocity boundary conditions on the solid surface. Based on the proposed diffuse-interface model, the fluid interface dynamics, the fluid-structure interaction, and the wetting property of the solid surface can be described simply and efficiently. Additionally, the total energy is also proved to be dissipative for the two-phase flows in the stationary geometries. To test the present diffuse-interface model, we develop a consistent and conservative lattice Boltzmann method and conduct some simulations. The numerical results also confirm the energy dissipation and good capability of the proposed diffuse-interface model in the study of two-phase flows in complex geometries and gas-liquid-particle multiphase flows., Comment: 27 pages, 23 figures

In this work, a ternary phase-field model for two-phase flows in complex geometries is proposed. In this model, one of the three components in the classical ternary Cahn-Hilliard model is considered as the solid phase, and only one Cahn-Hilliard equation with degenerate mobility needs to be solved due to the condition of volume conservation, which is consistent with the standard phase-field model with a single-scalar variable for two-phase flows. To depict different wetting properties at the complex fluid-solid boundaries, the spreading parameters in ternary phase-field model are determined based on the Young's law, in which the liquid-solid surface tension coefficient is assumed to be a linear function of gas-liquid surface tension coefficient and related to the contact angle and the minimum curvature of the solid surface. In addition, to achieve a high viscosity in the solid phase and preserve the velocity boundary conditions on the solid surface, the phase-field variable of the solid phase is also used to derive the modified Navier-Stokes equations. To test the present model, we further develop a consistent and conservative Hermite-moment based lattice Boltzmann method where an adjustable scale factor is introduced to improve the numerical stability, and conduct the numerical simulations of several benchmark problems. The results illustrate that present model has the good capability in the study of the two-phase flows in complex geometries., Comment: 23 pages, 25 figures

In this work, we conduct a detailed comparison between two second-order conservative Allen-Cahn (AC) models [\emph{Model A}: Zheng \emph{et al.}, Phys. Rev. E 101, 0433202 (2020) and \emph{Model B}: Mirjalili and Mani, (2023)] for the immiscible $N$-phase flows. Mathematically, these two AC equations can be proved to be equivalent under some approximate conditions. However, the effects of these approximations are unclear from the theoretical point of view, and would be considered numerically. To this end, we propose a consistent and conservative lattice Boltzmann method for the AC models for $N$-phase flows, and present some numerical comparisons of accuracy and stability between these two AC models. The results show that both two AC models have good performances in accuracy, but the \emph{Model B} is more stable for the realistic complex $N$-phase flows, although there is an adjustable parameter in the \emph{Model A}., Comment: 21 pages, 18 figures

In this work, we proposed a diffuse interface model for the dendritic growth with thermosolutal convection. In this model, the sharp boundary between the fluid and solid dendrite is replaced by a thin but nonzero thickness diffuse interface, which is described by the order parameter governed by the phase-field equation for the dendritic growth. The governing equations for solute and heat transfer are modified such that the previous special treatments for source term can be avoided. To solve the model for the dendritic growth with thermosolutal convection, we also developed a diffuse-interface multi-relaxation-time lattice Boltzmann (LB) method. In this method, the order parameter in the phase-field equation is combined into the force caused by the fluid-solid interaction, and the treatment on the complex fluid-solid interface can be avoided. In addition, four LB models are developed for the phase-field equation, concentration equation, temperature equation and the Navier-Stokes equations in a unified framework. Finally, to test the present diffuse-interface LB method, we performed some simulations of the dendritic growth, and found that the numerical results are in good agreements with some previous works., Comment: 20 pages, 14 figures

In this work, we consider a general consistent and conservative phase-field model for the incompressible two-phase flows. In this model, not only the Cahn-Hilliard or Allen-Cahn equation can be adopted, but also the mass and the momentum fluxes in the Navier-Stokes equations are reformulated such that the consistency of reduction, consistency of mass and momentum transport, and the consistency of mass conservation are satisfied. We further develop a lattice Boltzmann (LB) method and show that through the direct Taylor expansion, the present LB method can correctly recover the consistent and conservative phase-field model. Additionally, if the divergence of the extra momentum flux is seen as a force term, the extra force in the present LB method would include another term which has not been considered in the previous LB models. To quantitatively evaluate the incompressibility and the consistency of the mass conservation, two statistical variables are introduced in the study of the deformation of a square droplet, and the results show that the present LB method is more accurate. The layered Poiseuille flow and a droplet spreading on an ideal wall are further investigated, and the numerical results are in good agreement with the analytical solutions. Finally, the problems of the Rayleigh-Taylor instability and dam break with the high Reynolds numbers and/or large density ratios are studied, and it is found that the present consistent and conservative LB method is robust for such complex two-phase flows.

In this paper, a multiple-distribution-function lattice Boltzmann method (MDF-LBM) with multiple-relaxation-time model is proposed for incompressible Navier-Stokes equations (NSEs) which are considered as the coupled convection-diffusion equations (CDEs). Through direct Taylor expansion analysis, we show that the Navier-Stokes equations can be recovered correctly from the present MDF-LBM, and additionally, it is also found that the velocity and pressure can be directly computed through the zero and first-order moments of distribution function. Then in the framework of present MDF-LBM, we develop a locally computational scheme for the velocity gradient where the first-order moment of the non-equilibrium distribution is used, this scheme is also extended to calculate the velocity divergence, strain rate tensor, shear stress and vorticity. Finally, we also conduct some simulations to test the MDF-LBM, and find that the numerical results not only agree with some available analytical and numerical solutions, but also have a second-order convergence rate in space., Comment: 32 pages, 18 figures

In this paper, we propose a lattice Boltzmann (LB) model for the generalized coupled cross-diffusion-fluid system. Through the direct Taylor expansion method, the proposed LB model can correctly recover the macroscopic equations. The cross diffusion terms in the coupled system are modeled by introducing additional collision operators, which can be used to avoid special treatments for the gradient terms. In addition, the auxiliary source terms are constructed properly such that the numerical diffusion caused by the convection can be eliminated. We adopt the developed LB model to study two important systems, i.e., the coupled chemotaxis-fluid system and the double-diffusive convection system with Soret and Dufour effects. We first test the present LB model through considering a steady-state case of coupled chemotaxis-fluid system, then we analyze the influences of some physical parameters on the formation of sinking plumes. Finally, the double-diffusive natural convection system with Soret and Dufour effects is also studied, and the numerical results agree well with some previous works.