Multiple-relaxation-time lattice Boltzmann model for incompressible miscible flow with large viscosity ratio and high Péclet number

Phase-field lattice Boltzmann model with singular mobility for quasi-incompressible two-phase flows

In this paper, a lattice Boltzmann for quasi-incompressible two-phase flows is proposed based on the Cahn-Hilliard phase-field theory, which can be viewed as an improved model of a previous one [Yang and Guo, Phys. Rev. E 93, 043303 (2016)2470-004510.1103/PhysRevE.93.043303]. The model is composed of two LBE's, one for the Cahn-Hilliard equation (CHE) with a singular mobility, and the other for the quasi-incompressible Navier-Stokes equations (qINSE). Particularly, the LBE for the CHE uses an equilibrium distribution function containing a free parameter associated with the gradient of chemical potential, such that the variable (and even zero) mobility can be handled. In addition, the LBE for the qINSE uses an equilibrium distribution function containing another free parameter associated with the local shear rate, such that the large viscosity ratio problems can be handled. Several tests are first carried out to test the capability of the proposed LBE for the CHE in capturing phase interface, and the results demonstrate that the proposed model outperforms the original LBE model in terms of accuracy and stability. Furthermore, by coupling the hydrodynamic equations, the tests of double-stationary droplets and droplets falling problems indicate that the proposed model can reduce numerical dissipation and produce physically acceptable results at large time scales. The results of droplets falling and phase separation of binary fluid problems show that the present model can handle two-phase flows with large viscosity ratio up to the magnitude of 10^{4}.

Simplified wetting boundary scheme in phase-field lattice Boltzmann model for wetting phenomena on curved boundaries

In this work, a simplified wetting boundary scheme in the phase-field lattice Boltzmann model is developed for wetting phenomena on curved boundaries. The proposed scheme combines the advantages of the fluid-solid interaction scheme and geometric scheme-easy to implement (no need to interpolate the values of parameters exactly on solid boundaries and find proper characteristic vectors), the value of contact angle can be directly prescribed, and no unphysical spurious mass layer-and avoids mass leakage. Different from previous works, the values of the order parameter gradient on fluid boundary nodes are directly determined according to the geometric formulation rather than indirectly regulated through the order parameters on ghost solid nodes (i.e., ghost contact-line region). For this purpose, two numerical approaches to evaluate the order parameter gradient on fluid boundary nodes are utilized, one with the prevalent isotropic central scheme and the other with a local gradient scheme that utilizes the distribution functions. The simplified wetting boundary schemes with both numerical approaches are validated and compared through several numerical simulations. The results demonstrate that the proposed model has good ability and satisfactory accuracy to simulate wetting phenomena on curved boundaries under large density ratios.

Phase-field lattice Boltzmann model for two-phase flows with large density ratio

In this work, a lattice Boltzmann (LB) model based on the phase-field method is proposed for simulating large density ratio two-phase flows. An improved multiple-relaxation-time (MRT) LB equation is first developed to solve the conserved Allen-Cahn (AC) equation. By utilizing a nondiagonal relaxation matrix and modifying the equilibrium distribution function and discrete source term, the conserved AC equation can be correctly recovered by the proposed MRT LB equation with no deviation term. Therefore, the calculations of the temporal derivative term in the previous LB models are successfully avoided. Numerical tests demonstrate that satisfactory accuracy can be achieved by the present model to solve the conserved AC equation. What is more, the discrete force term of the MRT LB equation for the incompressible Navier-Stokes equations is also simplified and modified in the present work. An alternative scheme to calculate the gradient terms of the order parameter involved in the discrete force term through the nonequilibrium part of the distribution function is also developed. To validate the ability of the present LB model for simulating large density ratio two-phase flows, series of benchmarks, including two-phase Poiseuille flow, droplet impacting on thin liquid film, and planar Taylor bubble are simulated. It is found that the results predicted by the present LB model agree well with the analytical, numerical, and experimental results.

Multiple-relaxation-time lattice Boltzmann method for the Navier-Stokes and nonlinear convection-diffusion equations: Modeling, analysis, and elements

In this paper, we first present a unified framework of multiple-relaxation-time lattice Boltzmann (MRT-LB) method for the Navier-Stokes and nonlinear convection-diffusion equations where a block-lower-triangular-relaxation matrix and an auxiliary source distribution function are introduced. We then conduct a comparison of the four popular analysis methods (Chapman-Enskog analysis, Maxwell iteration, direct Taylor expansion, and recurrence equations approaches) that have been used to obtain the macroscopic Navier-Stokes and nonlinear convection-diffusion equations from the MRT-LB method and show that from mathematical point of view, these four analysis methods can give the same equations at the second-order of expansion parameters. Finally, we give some elements that are needed in the implementation of the MRT-LB method and also find that some available LB models can be obtained from this MRT-LB method.

Discrete effects on boundary conditions of the lattice Boltzmann method for fluid flows with curved no-slip walls

The lattice Boltzmann method (LBM) has been formulated as a powerful numerical tool to simulate incompressible fluid flows. However, it is still a critical issue for the LBM to overcome the discrete effects on boundary conditions successfully for curved no-slip walls. In this paper, we focus on the discrete effects of curved boundary conditions within the framework of the multiple-relaxation-time (MRT) model. We analyze in detail a single-node curved boundary condition [Zhao et al., Multiscale Model. Simul. 17, 854 (2019)10.1137/18M1201986] for predicting the Poiseuille flow and derive the numerical slip at the boundary dependent on a free parameter as well as the distance ratio and the relaxation times. An approach by virtue of the free parameter is then proposed to eliminate the slip velocity while with uniform relaxation parameters. The theoretical analysis also indicates that for previous curved boundary schemes only with the distance ratio and the halfway bounce-back (HBB) boundary scheme, the numerical slip cannot be removed with uniform relaxation times virtually. We further carried out some simulations to validate our theoretical derivations, and the numerical results for the case of straight and curved boundaries confirm our theoretical analysis. Finally, for fluid flows with curved boundary geometries, resorting to more degrees of freedom from the boundary scheme may have more potential to eliminate the discrete effect at the boundary with uniform relaxation times.

Stability of the lattice kinetic scheme and choice of the free relaxation parameter

The lattice kinetic scheme (LKS), a modified version of the classical single relaxation time (SRT) lattice Boltzmann method, was initially developed as a suitable numerical approach for non-Newtonian flow simulations and a way to reduce memory consumption of the original SRT approach. The better performances observed for non-Newtonian flows are mainly due to the additional degree of freedom allowing an independent control over the relaxation of higher-order moments, independently from the fluid viscosity. Although widely applied to fluid flow simulations, no theoretical analysis of LKS has been performed. The present work focuses on a systematic von Neumann analysis of the linearized collision operator. Thanks to this analysis, the effects of the modified collision operator on the stability domain and spectral behavior of the scheme are clarified. Results obtained in this study show that correct choices of the "second relaxation coefficient" lead, to a certain extent, to a more consistent dispersion and dissipation for large values of the first relaxation coefficient. Furthermore, appropriate values of this parameter can lead to a larger linear stability domain. At moderate and low values of viscosity, larger values of the free parameter are observed to increase dissipation of kinetic modes, while leaving the acoustic modes untouched and having a less pronounced effect on the convective mode. This increased dissipation leads in general to less pronounced sources of nonlinear instability, thus improving the stability of the LKS.

Theoretical and numerical analysis of the lattice kinetic scheme for complex-flow simulations

The lattice kinetic scheme (LKS) is a modified version of the classical single relaxation time lattice Boltzmann method. Although used for many applications, especially when large variations in viscosity are involved, a thorough analysis of the scheme has not been provided yet. In the context of this work, the macroscopic behavior of this scheme is evaluated through the Chapman-Enskog analysis. It is shown that the additional degree of freedom provided in the scheme allows for an independent control of higher-order moments. These results are further corroborated by numerical simulations. The behavior of this numerical scheme is studied for selected external and internal flows to clarify the effect of the free parameter on the different moments of the distribution function. It is shown that it is more stable than SRT (single relaxation time) when confronted to fully periodic under-resolved simulations (especially for λ≈1). It can also help minimize the error coming from the viscosity-dependence of the wall position when combined with the bounce-back approach; although still present, viscosity-dependence of the wall position is reduced. Furthermore, as shown through the multiscale analysis, specific choices of the free parameter can cancel out the leading-order error. Overall, the LKS is shown to be a useful and efficient alternative to the SRT method for simulating numerically complex flows.

Prediction of the moments in advection-diffusion lattice Boltzmann method. II. Attenuation of the boundary layers via double-Λ bounce-back flux scheme

Impact of the unphysical tangential advective-diffusion constraint of the bounce-back (BB) reflection on the impermeable solid surface is examined for the first four moments of concentration. Despite the number of recent improvements for the Neumann condition in the lattice Boltzmann method-advection-diffusion equation, the BB rule remains the only known local mass-conserving no-flux condition suitable for staircase porous geometry. We examine the closure relation of the BB rule in straight channel and cylindrical capillary analytically, and show that it excites the Knudsen-type boundary layers in the nonequilibrium solution for full-weight equilibrium stencil. Although the d2Q5 and d3Q7 coordinate schemes are sufficient for the modeling of isotropic diffusion, the full-weight stencils are appealing for their advanced stability, isotropy, anisotropy and anti-numerical-diffusion ability. The boundary layers are not covered by the Chapman-Enskog expansion around the expected equilibrium, but they accommodate the Chapman-Enskog expansion in the bulk with the closure relation of the bounce-back rule. We show that the induced boundary layers introduce first-order errors in two primary transport properties, namely, mean velocity (first moment) and molecular diffusion coefficient (second moment). As a side effect, the Taylor-dispersion coefficient (second moment), skewness (third moment), and kurtosis (fourth moment) deviate from their physical values and predictions of the fourth-order Chapman-Enskog analysis, even though the kurtosis error in pure diffusion does not depend on grid resolution. In two- and three-dimensional grid-aligned channels and open-tubular conduits, the errors of velocity and diffusion are proportional to the diagonal weight values of the corresponding equilibrium terms. The d2Q5 and d3Q7 schemes do not suffer from this deficiency in grid-aligned geometries but they cannot avoid it if the boundaries are not parallel to the coordinate lines. In order to vanish or attenuate the disparity of the modeled transport coefficients with the equilibrium weights without any modification of the BB rule, we propose to use the two-relaxation-times collision operator with free-tunable product of two eigenfunctions Λ. Two different values Λ_{v} and Λ_{b} are assigned for bulk and boundary nodes, respectively. The rationale behind this is that Λ_{v} is adjustable for stability, accuracy, or other purposes, while the corresponding Λ_{b}(Λ_{v}) controls the primary accommodation effects. Two distinguished but similar functional relations Λ_{b}(Λ_{v}) are constructed analytically: they preserve advection velocity in parabolic profile, exactly in the two-dimensional channel and very accurately in a three-dimensional cylindrical capillary. For any velocity-weight stencil, the (local) double-Λ BB scheme produces quasi-identical solutions with the (nonlocal) specular-forward reflection for first four moments in a channel. In a capillary, this strategy allows for the accurate modeling of the Taylor-dispersion and non-Gaussian effects. As illustrative example, it is shown that in the flow around a circular obstacle, the double-Λ scheme may also vanish the dependency of mean velocity on the velocity weight; the required value for Λ_{b}(Λ_{v}) can be identified in a few bisection iterations in given geometry. A positive solution for Λ_{b}(Λ_{v}) may not exist in pure diffusion, but a sufficiently small value of Λ_{b} significantly reduces the disparity in diffusion coefficient with the mass weight in ducts and in the presence of rectangular obstacles. Although Λ_{b} also controls the effective position of straight or curved boundaries, the double-Λ scheme deals with the lower-order effects. Its idea and construction may help understanding and amelioration of the anomalous, zero- and first-order behavior of the macroscopic solution in the presence of the bulk and boundary or interface discontinuities, commonly found in multiphase flow and heterogeneous transport.

Boundary scheme for linear heterogeneous surface reactions in the lattice Boltzmann method

A boundary scheme is proposed to treat the linear heterogeneous surface reaction (i.e., general Robin boundary condition) for the lattice Boltzmann method (LBM) in this study. The basic idea of the present scheme is to compute the scalar variable gradient in the boundary condition with the moment of the nonequilibrium distribution functions. The unknown distribution functions on the boundary can then be obtained on the basis of the given boundary condition. For a straight wall, where the lattice nodes are located on the boundaries, both the scalar variable and its gradient can be expressed in terms of the distribution functions at the boundary node, and the scheme is purely local. The scheme is also extended to problems with a curved wall, in which a linear extrapolation is employed to realize the exact boundary position. A common feature of the two schemes lies in the easy treatment of the heterogenous reactions in comparison with existing methods. A number of simulations are performed to test the accuracy of the schemes. The results show that for flat walls the scheme can achieve second-order accuracy, while for curved walls the order of the accuracy is between 1.0 and 1.5.