Axisymmetric multiphase lattice Boltzmann method

Study of phase-field lattice Boltzmann models based on the conservative Allen-Cahn equation

Conservative phase-field (CPF) equations based on the Allen-Cahn model for interface tracking in multiphase flows have become more popular in recent years, especially in the lattice-Boltzmann (LB) community. This is largely due to their simplicity and improved efficiency and accuracy over their Cahn-Hilliard-based counterparts. Additionally, the improved locality of the resulting LB equation (LBE) for the CPF models makes them more ideal candidates for LB simulation of multiphase flows on nonuniform grids, particularly within an adaptive-mesh refinement framework and massively parallel implementation. In this regard, some modifications-intended as improvements-have been made to the original CPF-LBE proposed by Geier et al. [Phys. Rev. E 91, 063309 (2015)PLEEE81539-375510.1103/PhysRevE.91.063309] which require further examination. The goal of the present study is to conduct a comparative investigation into the differences between the original CPF model proposed by Geier et al. [Phys. Rev. E 91, 063309 (2015)PLEEE81539-375510.1103/PhysRevE.91.063309] and the so-called improvements proposed by Ren et al. [Phys. Rev. E 94, 023311 (2016)2470-004510.1103/PhysRevE.94.023311] and Wang et al. [Phys. Rev. E 94, 033304 (2016)2470-004510.1103/PhysRevE.94.033304]. Using the Chapman-Enskog analysis, we provide a detailed derivation of the governing equations in each model and then examine the efficacy of the above-mentioned models for some benchmark problems. Several test cases have been designed to study different configurations ranging from basic yet informative flows to more complex flow fields, and the results are compared with finite-difference simulations. Furthermore, as a development of the previously proposed CPF-LBE model, axisymmetric formulations for the proposed model by Geier et al. [Phys. Rev. E 91, 063309 (2015)PLEEE81539-375510.1103/PhysRevE.91.063309] are derived and presented. Finally, two benchmark problems are designed to compare the proposed axisymmetric model with the analytical solution and previous work. We find that the accuracy of the model for interface tracking is roughly similar for different models at high viscosity ratios, high density ratios, and relatively high Reynolds numbers, while the original CFP-LBE without the additional time-dependent terms outperforms the so-called improved models in terms of efficiency, particularly on distributed parallel machines.

Axisymmetric lattice Boltzmann model for simulating the freezing process of a sessile water droplet with volume change

Droplet freezing not only is of fundamental interest but also plays an important role in numerous natural and industrial processes. However, it is challenging to numerically simulate the droplet freezing process due to its involving a complex three-phase system with dynamic phase change and heat transfer. Here we propose an axisymmetric lattice Boltzmann (LB) model to simulate the freezing process of a sessile water droplet with consideration of droplet volume expansion. Combined with the multiphase flow LB model and the enthalpy thermal LB model, our proposed approach is applied to simulate the sessile water droplet freezing on both hydrophilic and hydrophobic surfaces at a fixed subcooled temperature. Through comparison with the experimental counterpart, the comparison results show that our axisymmetric LB model can satisfactorily describe such sessile droplet freezing processes. Moreover, we use both LB simulations and analytical models to study the effects of contact angle and volume expansion on the freezing time and the cone shape formed on the top of frozen droplets. The analytical models are obtained based on heat transfer and geometric analyses. Additionally, we show analytically and numerically that the freezing front-to-interface angle keeps nearly constant (smaller than 90°).

Axisymmetric compact finite-difference lattice Boltzmann method for blood flow simulations

An axisymmetric compact finite-difference lattice Boltzmann method is proposed to simulate both Newtonian and non-Newtonian flow of blood through a lumen. The curvature of the arteries could be accurately resolved using body-fitted mesh owing to the proposed finite-difference formulation. The axisymmetric nature of the flow, as well as the non-Newtonian nature of blood, are incorporated into the lattice Boltzmann equation using separate source terms. Using Chapman-Enskog expansion it is shown that the resulting lattice Boltzmann equation with these additional source terms recovers the macroscopic axisymmetric hydrodynamic equations. The solver is verified for (1) steady inflow of a Newtonian fluid through a stenosed lumen, (2) temporally developing pulsatile flow (Womersley flow) through a straight lumen with Newtonian fluid, and (3) steady inflow of a non-Newtonian fluid through a straight lumen. The solver is then applied to simulate the steady flow of a non-Newtonian fluid through a stenosed lumen, and it was found that a smaller recirculation zone and lower WSS values are obtained when compared with the flow of a Newtonian fluid. The capability of the solver to simulate spatially developing (velocity-driven) pulsatile flow is then demonstrated by simulating physiological pulsatile flow through an axisymmetric abdominal aortic aneurysm. From this simulation, the cycle-averaged wall shear stress is observed to have a steep gradient going from a minimum (negative) to a maximum (positive) value towards the distal end of the aneurysm, which is prone to the risk of rupture. An iterative procedure to select the geometric and flow parameters for unsteady inflow condition in the lattice Boltzmann method framework is demonstrated that accurately resolves all the timescales to achieve incompressibility. Overall, the present solver seems to be promising to simulate axisymmetric flow of blood with steady and pulsatile inflows while considering the blood rheology.

Hydrodynamic behavior of the pseudopotential lattice Boltzmann method for interfacial flows

The lattice Boltzmann method (LBM) is routinely employed in the simulation of complex multiphase flows comprising bulk phases separated by nonideal interfaces. The LBM is intrinsically mesoscale with a hydrodynamic equivalence popularly set by the Chapman-Enskog analysis, requiring that fields slowly vary in space and time. The latter assumptions become questionable close to interfaces where the method is also known to be affected by spurious nonhydrodynamical contributions. This calls for quantitative hydrodynamical checks. In this paper, we analyze the hydrodynamic behavior of the LBM pseudopotential models for the problem of the breakup of a liquid ligament triggered by the Plateau-Rayleigh instability. Simulations are performed at fixed interface thickness, while increasing the ligament radius, i.e., in the "sharp interface" limit. The influence of different LBM collision operators is also assessed. We find that different distributions of spurious currents along the interface may change the outcome of the pseudopotential model simulations quite sensibly, which suggests that a proper fine-tuning of pseudopotential models in time-dependent problems is needed before the utilization in concrete applications. Taken all together, we argue that the results of the proposed paper provide a valuable insight for engineering pseudopotential model applications involving the hydrodynamics of liquid jets.

Effects of thermal fluctuations in the fragmentation of a nanoligament

We study the effects of thermally induced capillary waves in the fragmentation of a liquid ligament into multiple nanodroplets. Our numerical implementation is based on a fluctuating lattice Boltzmann (LB) model for nonideal multicomponent fluids, including nonequilibrium stochastic fluxes mimicking the effects of molecular forces at the nanoscales. We quantitatively analyze the statistical distribution of the breakup times and the droplet volumes after the fragmentation process at changing the two relevant length scales of the problem, i.e., the thermal length scale and the ligament size. The robustness of the observed findings is also corroborated by quantitative comparisons with the predictions of sharp interface hydrodynamics. Beyond the practical importance of our findings for nanofluidic engineering devices, our study also explores a novel application of LB in the realm of nanofluidic phenomena.

Alternative extrapolation-based symmetry boundary implementations for the axisymmetric lattice Boltzmann method

In this study, alternative symmetry boundary implementations for the axisymmetric lattice Boltzmann (LB) method are proposed based on the nonequilibrium extrapolation and the direct extrapolation schemes. The proposed boundary schemes are directly implemented on the symmetry axis, and the postcollision distribution function and the macroscopic variables at the boundary nodes are extrapolated from the inner fluid nodes; thereby, the singularities arising at the symmetry axis (r=0) during the collision and the macroscopic variable calculations are completely avoided. The accuracy of the present schemes is consistent with the well-established axisymmetric LB model. Moreover, in comparison with previous symmetry boundary schemes, the present implementations are slightly more accurate than the symmetry scheme by Guo et al. [Phys. Rev. E 79, 046708 (2009)10.1103/PhysRevE.79.046708] and numerically more stable than the specular reflection-based schemes.

Forcing scheme analysis for the axisymmetric lattice Boltzmann method under incompressible limit

Because the standard lattice Boltzmann (LB) method is proposed for Cartesian Navier-Stokes (NS) equations, additional source terms are necessary in the axisymmetric LB method for representing the axisymmetric effects. Therefore, the accuracy and applicability of the axisymmetric LB models depend on the forcing schemes adopted for discretization of the source terms. In this study, three forcing schemes, namely, the trapezium rule based scheme, the direct forcing scheme, and the semi-implicit centered scheme, are analyzed theoretically by investigating their derived macroscopic equations in the diffusive scale. Particularly, the finite difference interpretation of the standard LB method is extended to the LB equations with source terms, and then the accuracy of different forcing schemes is evaluated for the axisymmetric LB method. Theoretical analysis indicates that the discrete lattice effects arising from the direct forcing scheme are part of the truncation error terms and thus would not affect the overall accuracy of the standard LB method with general force term (i.e., only the source terms in the momentum equation are considered), but lead to incorrect macroscopic equations for the axisymmetric LB models. On the other hand, the trapezium rule based scheme and the semi-implicit centered scheme both have the advantage of avoiding the discrete lattice effects and recovering the correct macroscopic equations. Numerical tests applied for validating the theoretical analysis show that both the numerical stability and the accuracy of the axisymmetric LB simulations are affected by the direct forcing scheme, which indicate that forcing schemes free of the discrete lattice effects are necessary for the axisymmetric LB method.

Analytical relations for long-droplet breakup in asymmetric T junctions

We develop accurate analytical relations for the droplet volume ratio, droplet length during breakup process, and pressure drop of asymmetric T junctions with a valve in each of the branches for producing unequal-sized droplets. An important advantage of this system is that after manufacturing the system, the size of the generated droplets can be changed simply by adjusting the valves. The results indicate that if the valve ratio is smaller than 0.65, the system enters a nonbreakup regime. Also the pressure drop does not depend on the time and decreases by increasing the valve ratio, namely, opening the degree of valve 1 to valve 2. In addition, the results reveal that by decreasing (increasing) the valve ratio, the droplet length of branch 1 decreases (increases) and the droplet length of branch 2 increases (decreases) linearly while the whole length of the droplet remains unchanged.

Phase-field-based lattice Boltzmann model for axisymmetric multiphase flows

In this paper, a phase-field-based lattice Boltzmann (LB) model is proposed for axisymmetric multiphase flows. Modified equilibrium distribution functions and some source terms are properly added into the evolution equations such that multiphase flows in the axisymmetric coordinate system can be described. Different from previous axisymmetric LB multiphase models, the added source terms that arise from the axisymmetric effect contain no additional gradients, and therefore the present model is much simpler. Furthermore, through the Chapmann-Enskog analysis, the axisymmetric Chan-Hilliard equation and Navier-Stokes equations can be exactly derived from the present model. The model is also able to deal with flows with density contrast. A variety of numerical experiments, including planar and curve interfaces, an elongation field, a static droplet, a droplet oscillation, breakup of a liquid thread, and dripping of a liquid droplet under gravity, have been conducted to test the proposed model. It is found that the present model can capture accurate interface and the numerical results of multiphase flows also agree well with the analytical solutions and/or available experimental data.