Three-dimensional lattice-Boltzmann model of van der Waals fluids

Thermodynamic consistency of a pseudopotential lattice Boltzmann fluid with interface curvature

Thermodynamic consistency of pseudopotential lattice Boltzmann models is a major topic that needs comprehensive evaluations. When interface is flat, pseudopotential models can give density-pressure isotherms in excellent agreement with those from equation of state. When interface is curved, thermodynamic equilibriums are affected by interface curvature, and consistency of pseudopotential models has not been systematically evaluated. In this study, we show that the effect of Laplace pressure on phase equilibrium is quantitatively consistent with Kelvin equation at high reduced temperatures (≥0.7). At low temperatures, inconsistency that can be attributed to the effect of orientation of the interface was noted, and it can be improved by tuning of the pseudopotential. By relating interfacial tension of a simulated fluid to that of a real fluid, the lattice spacing of pseudopotential model is found to be on the order of several molecular diameters, the typical range of intermolecular interactions. Interfacial thickness at different temperatures in pseudopotential model compared well with experiments and molecular dynamics simulations, which confirms that the calculated length scale is reasonable. Evaluation of a free energy lattice Boltzmann model indicate that it is consistent with Kelvin equation at high temperatures. The free energy model, however, is not as accurate as the tested pseudopotential model, and discrepancies may come from the relative inaccuracies in the predictions of vapor densities and the thinner interfaces.

Lattice Boltzmann model with self-tuning equation of state for multiphase flows

A lattice Boltzmann (LB) model for multiphase flows is developed that complies with the thermodynamic foundations of kinetic theory. By directly devising the collision term for the LB equation at the discrete level, a self-tuning equation of state is achieved, which can be interpreted as the incorporation of short-range molecular interaction. A pairwise interaction force is introduced to mimic the long-range molecular interaction, which is responsible for interfacial dynamics. The derived pressure tensor is naturally consistent with thermodynamic theory, and surface tension and interface thickness can be independently prescribed. Numerical tests, including static and dynamic cases, are carried out to validate the present model and good results are obtained. As a further application, head-on collision of equal-sized droplets is simulated and the elusive "bouncing" regime is successfully reproduced.

Wettability effects on contact line dynamics of droplet motion in an inclined channel

The present work discusses the implications of wall wettability and inclination of the surface on droplet dynamics. In this work, mesoscopic illustration of droplet dynamics in a duct with different inclination angles, based on the two-phase lattice Boltzmann model, is reported. Temporal evolution of wetted length, wetted area, and maximum height of the droplet for surfaces with different inclination angles and wettabilities is furnished in detail in order to elucidate the droplet displacement dynamics. It has been observed that the effect of inclination of the surface on droplet dynamics is more pronounced on a hydrophobic surface as compared to a hydrophilic surface. The time evolution of height and contact line motion of the droplet shows that higher angle of inclination of substrate affects the dynamics strongly irrespective of wettability.

Phase-field-based lattice Boltzmann model for incompressible binary fluid systems with density and viscosity contrasts

A lattice Boltzmann model (LBM) is proposed based on the phase-field theory to simulate incompressible binary fluids with density and viscosity contrasts. Unlike many existing diffuse interface models which are limited to density matched binary fluids, the proposed model is capable of dealing with binary fluids with moderate density ratios. A new strategy for projecting the phase field to the viscosity field is proposed on the basis of the continuity of viscosity flux. The new LBM utilizes two lattice Boltzmann equations (LBEs): one for the interface tracking and the other for solving the hydrodynamic properties. The LBE for interface tracking can recover the Chan-Hilliard equation without any additional terms; while the LBE for hydrodynamic properties can recover the exact form of the divergence-free incompressible Navier-Stokes equations avoiding spurious interfacial forces. A series of 2D and 3D benchmark tests have been conducted for validation, which include a rigid-body rotation, stationary and moving droplets, a spinodal decomposition, a buoyancy-driven bubbly flow, a layered Poiseuille flow, and the Rayleigh-Taylor instability. It is shown that the proposed method can track the interface with high accuracy and stability and can significantly and systematically reduce the parasitic current across the interface. Comparisons with momentum-based models indicate that the newly proposed velocity-based model can better satisfy the incompressible condition in the flow fields, and eliminate or reduce the velocity fluctuations in the higher-pressure-gradient region and, therefore, achieve a better numerical stability. In addition, the test of a layered Poiseuille flow demonstrates that the proposed scheme for mixture viscosity performs significantly better than the traditional mixture viscosity methods.

Thermodynamic consistency of liquid-gas lattice Boltzmann methods: interfacial property issues

In the present study we examine the thermodynamic consistency of lattice Boltzmann equation (LBE) models that are based on the forcing method by comparing different numerical treatments of the LBE for van der Waals fluids. The different models are applied for the calculation of bulk and interfacial thermodynamic properties at various temperatures. The effect of the interface density gradient parameter, kappa , that controls surface tension, is related explicitly with the fluid characteristics, including temperature, molecular diameter, and lattice spacing, through the employment of a proper intermolecular interaction potential. A comprehensive analysis of the interfacial properties reveals some important shortcomings of the LBE methods when central finite difference schemes are employed in the directional derivative calculations and proposes a proper treatment that ensures thermodynamically consistent interfacial properties in accord with the van der Waals theory. The results are found to be in excellent quantitative agreement with exact results of the van der Waals theory preserving all the major features of the interfacial characteristics of vapor-liquid systems of different shapes and sizes.

Interface width and bulk stability: requirements for the simulation of deeply quenched liquid-gas systems

Simulations of liquid-gas systems with interface terms evaluated by central difference discretizations are observed to fail to give accurate results for two reasons: the interface can get "stuck" on the lattice or a density overshoot develops around the interface. In the first case, the bulk densities can take a range of values, dependent on the initial conditions. In the second case, inaccurate bulk densities are found. We derived the minimum interface width required for the accurate simulation of liquid-gas systems with a diffuse interface. This criterion is demonstrated for lattice Boltzmann simulations of a van der Waals gas. Combining this criterion with predictions for the bulk stability defines the parameter range for stable and accurate simulation results even for high density ratios of over 1000.

Simulation of droplet formation and coalescence using lattice Boltzmann-based single-phase model

A lattice Boltzmann method-based single-phase free surface model is developed to study the interfacial dynamics of coalescence, droplet formation and detachment phenomena related to surface tension and wetting effects. Compared with the conventional multiphase models, the lattice Boltzmann-based single-phase model has a higher computational efficiency since it is not necessary to simulate the motion of the gas phase. A perturbation, which is given in the same fashion as the perturbation step in Gunstensen's color model, is added to the distribution functions of the interface cells for incorporating the surface tension into the single-phase model. The assignment of different mass gradients along the fluid-wall interface is used to model the wetting properties of the solid surface. Implementations of the model are demonstrated for simulating the processes of the droplet coalescence, the droplet formation and detachment from ceiling and from nozzles with different shapes and different wall wetting properties.

Lattice-Boltzmann model based on field mediators for immiscible fluids

In this paper, a lattice BGK (Bhatnagar-Gross-Krook) model is proposed for immiscible fluids. Collision operator is decoupled considering mutual and cross collisions between lattice particles, with three independent parameters related to the species diffusivity and to the viscosity of each fluid. Field mediator's concept, described by Santos and Philippi [Phys. Rev. E 65, 046305 (2002)], is extended to the framework of the lattice-Boltzmann equation and interference between mediators and particles is modeled by considering that there is a deviation in particles velocity, proportional to the mediators' distribution at the site. A Chapman-Enskog analysis is performed leading to theoretical predictions of the macroscopic equations inside the transition layer and to the transition-layer thickness. Chapman-Enskog analysis is restricted to near-equilibrium states and was unable to predict the correct second-order interfacial tension dependence on the modeled long-range fields intensity. Interfacial tension was, only, correctly retrieved using a nonequilibrium solution. Theoretical predictions are compared with simulation results and the model is tested considering its ability in describing the dynamical behavior of the interface and Galilean invariance.