Entropic lattice boltzmann method for multiphase flows

Two-component lattice Boltzmann model for solute transport in bubbly flows

A free energy lattice Boltzmann model has been developed to describe a binary system consisting of a nonideal solvent and a gaseous ideal solute. A free energy functional to describe this mixture has been derived, and it is used to determine the driving forces for two lattice Boltzmann equations, one for each component. The well-balanced lattice Boltzmann method is used to avoid discretization errors, which allows correct thermodynamic equilibrium to be achieved for both components as well as high density ratios. In this model, the distribution of the solute between the liquid and vapor phases of the solvent follows Henry's law due to the contributions of the two components to the free energy density. The momenta of the two components are coupled through the use of a mixture velocity in the equilibrium distributions of both components. The model also includes surface tension due to gradients of the solvent density and diffusion due to an added mobility term. The effects of model parameters on phase composition, surface tension, and the rate of solute transport are characterized, with examples of static gas bubbles and flat interfaces demonstrated. Mass transfer of the soluble component in the liquid phase of the nonideal component is characterized, and an equation for the solute diffusion coefficient is provided.

A filter inspired by deep-sea glass sponges for oil cleanup under turbulent flow

Oil spill disasters lead to widespread and long-lasting social, economical, environmental and ecological impacts. Technical challenges remain for conventional static adsorption due to hydrodynamic instability under complex water-flow conditions, which results in low oil-capture efficiency, time delay and oil escape. To address this issue, we design a vortex-anchored filter inspired by the anatomy of deep-sea glass sponges (E. aspergillum) by mimicking their exceptional skeletal features and filter-feeding patterns. Results demonstrate that the vortex-anchored filter can retain external turbulent-flow kinetic energy in low-speed vortical flow with small Kolmogorov microscale (85 μm) in the cavity of skeleton, leading to enhanced interfacial mass transfer and residence time by physical field synergy. It improves hydrodynamic stability by reducing Reynolds stresses in nearly quiescent wake flow. The vortex-anchored filter can realize >97% capture of floating, underwater and emulsified oils stably at Reynolds numbers ranging from subcritical to supercritical regimes. This study not only highlights the importance of vortex-anchored mechanism in enhancing interfacial mass transfer and hydrodynamic stability during oil capture beyond previously known benefits of increased residence time, but also represents a paradigm shift to advance biophysically inspired strategies for in-situ, dynamic and robust cleanup of spilled oil, environmental remediation and resource recovery.

Wetting on silicone surfaces

Silicone is frequently used as a model system to investigate and tune wetting on soft materials. Silicone is biocompatible and shows excellent thermal, chemical, and UV stability. Moreover, the mechanical properties of the surface can be easily varied by several orders of magnitude in a controlled manner. Polydimethylsiloxane (PDMS) is a popular choice for coating applications such as lubrication, self-cleaning, and drag reduction, facilitated by low surface energy. Aiming to understand the underlying interactions and forces, motivated numerous and detailed investigations of the static and dynamic wetting behavior of drops on PDMS-based surfaces. Here, we recognize the three most prevalent PDMS surface variants, namely liquid-infused (SLIPS/LIS), elastomeric, and liquid-like (SOCAL) surfaces. To understand, optimize, and tune the wetting properties of these PDMS surfaces, we review and compare their similarities and differences by discussing (i) the chemical and molecular structure, and (ii) the static and dynamic wetting behavior. We also provide (iii) an overview of methods and techniques to characterize PDMS-based surfaces and their wetting behavior. The static and dynamic wetting ridge is given particular attention, as it dominates energy dissipation, adhesion, and friction of sliding drops and influences the durability of the surfaces. We also discuss special features such as cloaking and wetting-induced phase separation. Key challenges and opportunities of these three surface variants are outlined.

Asymptotic freedom in the lattice Boltzmann theory

Asymptotic freedom is a feature of quantum chromodynamics that guarantees its well posedness. We derive an analog of asymptotic freedom enabling unconditional linear stability of lattice Boltzmann simulation of hydrodynamics. We further demonstrate the validity of the derived conditions via the special case of the equilibrium based on entropy maximization, which is shown to be uniquely renormalizable.

Three-dimensional lattice Boltzmann model with self-tuning equation of state for multiphase flows

In this work, the recent lattice Boltzmann (LB) model with self-tuning equation of state (EOS) [Huang et al., Phys. Rev. E 99, 023303 (2019)2470-004510.1103/PhysRevE.99.023303] is extended to three dimensions for the simulation of multiphase flows, which is based on the standard three-dimensional 27-velocity lattice and multiple-relaxation-time collision operator. To achieve the self-tuning EOS, the equilibrium moment is devised by introducing a built-in variable, and the collision matrix is improved by introducing some velocity-dependent nondiagonal elements. Meanwhile, the additional cubic terms of velocity in recovering the Newtonian viscous stress are eliminated to enhance the numerical accuracy. For modeling multiphase flows, an attractive pairwise interaction force is introduced to mimic the long-range molecular interaction, and a consistent scheme is proposed to compensate for the ɛ^{3}-order discrete lattice effect. Thermodynamic consistency in a strict sense is established for the multiphase LB model with self-tuning EOS, and the wetting condition is also treated in a thermodynamically consistent manner. As a result, the contact angle, surface tension, and interface thickness can be independently adjusted in the present theoretical framework. Numerical tests are first performed to validate the multiphase LB model with self-tuning EOS and the theoretical analyses of bulk and surface thermodynamics. The collision of equal-sized droplets is then simulated to demonstrate the applicability and effectiveness of the present LB model for multiphase flows.

Lattice Boltzmann simulation of deformable fluid-filled bodies: progress and perspectives

With the rapid development of studies involving droplet microfluidics, drug delivery, cell detection, and microparticle synthesis, among others, many scientists have invested significant efforts to model the flow of these fluid-filled bodies. Motivated by the intricate coupling between hydrodynamics and the interactions of fluid-filled bodies, several methods have been developed. The objective of this review is to present a compact foundation of the methods used in the literature in the context of lattice Boltzmann methods. For hydrodynamics, we focus on the lattice Boltzmann method due to its specific ability to treat time- and spatial-dependent boundary conditions and to incorporate new physical models in a computationally efficient way. We split the existing methods into two groups with regard to the interfacial boundary: fluid-structure and fluid-fluid methods. The fluid-structure methods are characterised by the coupling between fluid dynamics and mechanics of the flowing body, often used in applications involving membranes and similar flexible solid boundaries. We further divide fluid-structure-based methods into two subcategories, those which treat the fluid-structure boundary as a continuum medium and those that treat it as a discrete collection of individual springs and particles. Next, we discuss the fluid-fluid methods, particularly useful for the simulations of fluid-fluid interfaces. We focus on models for immiscible droplets and their interaction in a suspending fluid and describe benchmark tests to validate the models for fluid-filled bodies.

Particles on demand method: Theoretical analysis, simplification techniques, and model extensions

The particles on demand method [Phys. Rev. Lett. 121, 130602 (2018)0031-900710.1103/PhysRevLett.121.130602] was recently formulated with a conservative finite-volume discretization and validated against challenging benchmarks. In this work, we focus on the properties of the reference frame transformation and its implications on the accuracy of the model. Based on these considerations, we propose strategies that simplify the scheme and generalize it to include a tunable Prandtl number via quasi-equilibrium relaxation. Finally, we adapt concepts from the multiscale semi-Lagrangian lattice Boltzmann formulation to the proposed framework, further improving the potential and the operating range of the kinetic model. Numerical simulations of high Mach compressible flows demonstrate excellent accuracy and stability of the model over a wide range of conditions.

Unrestricted component count in multiphase lattice Boltzmann: A fugacity-based approach

Studies of multiphase fluids utilizing the lattice Boltzmann method (LBM) are typically severely restricted by the number of components or chemical species being modeled. This restriction is particularly pronounced for multiphase systems exhibiting partial miscibility and significant interfacial mass exchange, which is a common occurrence in realistic multiphase systems. Modeling such systems becomes increasingly complex as the number of chemical species increases due to the increased role of molecular interactions and the types of thermodynamic behavior that become possible. The recently introduced fugacity-based LBM [Soomro et al., Phys. Rev. E 107, 015304 (2023)2470-004510.1103/PhysRevE.107.015304] has provided a thermodynamically consistent modeling platform for multicomponent, partially miscible LBM simulations. However, until now, this fugacity-based LB model had lacked a comprehensive demonstration of its ability to accurately reproduce thermodynamic behavior beyond binary mixtures and to remove any restrictions in a number of components for multiphase LBM. In this paper we closely explore these fugacity-based LBM capabilities by showcasing comprehensive, thermodynamically consistent simulations of multiphase mixtures of up to ten chemical components. The paper begins by validating the model against the Young-Laplace equation for a droplet composed of three components. The model is then applied to study mixtures with a range of component numbers from one to six, showing agreement with rigorous thermodynamic predictions and demonstrating linear scaling of computational time with the number of components. We further investigate ternary systems in detail by exploring a wide range of temperature, pressure, and overall composition conditions to produce various characteristic ternary diagrams. In addition, the model is shown to be unrestricted in the number of phases as demonstrated through simulations of a three-component three-phase equilibrium case. The paper concludes by demonstrating simulations of a ten-component, realistic hydrocarbon mixture, achieving excellent agreement with thermodynamics for both flat interface vapor-liquid equilibrium and curved interface spinodal decomposition cases. This study represents a significant expansion of the scope and capabilities of multiphase LBM simulations that encompass multiphase systems of keen interest in engineering.

Entropic equilibrium for the lattice Boltzmann method: Hydrodynamics and numerical properties

The entropic lattice Boltzmann framework proposed the construction of the equilibrium by taking into consideration minimization of a discrete entropy functional. The effect of this entropic equilibrium on properties of the resulting solver has been the topic of discussions in the literature. Here we present a rigorous analysis of the hydrodynamics and numerics of the entropic equilibrium. We demonstrate that the entropic equilibrium features unconditional linear stability, in contrast to the conventional polynomial equilibrium. We reveal the mechanisms through which unconditional linear stability is maintained, most notable of which are adaptive propagation velocity of normal modes and the positive-definite nature of the dissipation rates of hydrodynamic eigenmodes. We further present a simple local correction to considerably reduce the deviations in the effective bulk viscosity.

Equation-of-state-dependent surface free-energy density for wettability in lattice Boltzmann method

In thermodynamic theory, the liquid-vapor fluids can be described by a single multiphase equation of state and the surface wettability is usually characterized by the surface free-energy density. In this work, we propose an equation-of-state-dependent surface free-energy density for the wettability of the liquid-vapor fluids on a solid surface, which can lead to a simple closed-form analytical expression for the contact angle. Meanwhile, the thermodynamically derived equilibrium condition is equivalent to the geometric formulation of the contact angle. To numerically validate the present surface free-energy density, the mesoscopic multiphase lattice Boltzmann model with self-tuning equation of state, which is strictly consistent with thermodynamic theory, is employed, and the two-dimensional wetting condition treatment is extended to the three-dimensional situation with flat and curved surfaces. Two- and three-dimensional lattice Boltzmann simulations of static droplets on flat and curved surfaces are first performed, and the obtained contact angles agree well with the closed-form analytical expression. Then, the three-dimensional lattice Boltzmann simulation of a moving droplet on an inclined wall, which is vertically and sinusoidally oscillated, is carried out. The dynamic contact angles well satisfy the Cox-Voinov law. The droplet movement regimes are consistent with previous experiments and two-dimensional simulations. The dependence of the droplet overall velocity with respect to the dimensionless oscillation strength is also discussed in detail.

Fugacity-based lattice Boltzmann method for multicomponent multiphase systems

The free-energy model can extend the lattice Boltzmann method to multiphase systems. However, there is a lack of models capable of simulating multicomponent multiphase fluids with partial miscibility. In addition, existing models cannot be generalized to honor thermodynamic information provided by any multicomponent equation of state of choice. In this paper, we introduce a free-energy lattice Boltzmann model where the forcing term is determined by the fugacity of the species, the thermodynamic property that connects species partial pressure to chemical potential calculations. By doing so, we are able to carry out multicomponent multiphase simulations of partially miscible fluids and generalize the methodology for use with any multicomponent equation of state of interest. We test this fugacity-based lattice Boltzmann method for the cases of vapor-liquid equilibrium for two- and three-component mixtures in various temperature and pressure conditions. We demonstrate that the model is able to reliably reproduce phase densities and compositions as predicted by multicomponent thermodynamics and can reproduce different characteristic pressure-composition and temperature-composition envelopes with a high degree of accuracy. We also demonstrate that the model can offer accurate predictions under dynamic conditions.

Equations of state in multiphase lattice Boltzmann method revisited

The single-component multiphase fluids can be described by a single equation of state (EOS), and various EOSs have been employed in the multiphase lattice Boltzmann (LB) method. In this work, we revisit five commonly used EOSs, including the van der Waals EOS, the Redlich-Kwong EOS, the Redlich-Kwong-Soave EOS, the Peng-Robinson EOS, and the Carnahan-Starling EOS. The recent multiphase LB model with self-tuning EOS is employed because of its thermodynamic consistency in a strict sense and clear physical picture at the microscopic level. First, the way to incorporate these multiphase EOSs is proposed. Two scaling factors are introduced to independently adjust the surface tension and interface thickness, and the lattice sound speed is EOS-dependent to ensure the numerical stability. Then, numerical tests are conducted to validate the incorporations of these EOSs and compare their numerical performances. The surface tension and interface thickness are set to the same values for different EOSs in the comparisons. The liquid and gas densities, surface tension, and interface thickness by the LB simulation agree well with the thermodynamic results. The maximum density ratios achieved with different EOSs are at the same level and could be very close to each other when the interface thickness is relatively small. The effects of multiphase EOS, density ratio, and dimensionless relaxation time on the spurious current are discussed in detail. It is interesting to find the van der Waals EOS shows the best numerical performance in reducing the spurious current.

Essentially entropic lattice Boltzmann model: Theory and simulations

We present a detailed description of the essentially entropic lattice Boltzmann model. The entropic lattice Boltzmann model guarantees unconditional numerical stability by iteratively solving the nonlinear entropy evolution equation. In this paper we explain the construction of closed-form analytic solutions to this equation. We demonstrate that near equilibrium this analytic solution reduces to the standard lattice Boltzmann model. We consider a few test cases to show that the analytic solution does not exhibit any significant deviation from the iterative solution. We also extend the analytical solution for the Ellipsoidal Statistical (ES)-Bhatnagar-Gross-Krook model to remove the limitation on the Prandtl number for heat transfer problems. The simplicity of the analytic solution removes the computational overhead and algorithmic complexity associated with the entropic lattice Boltzmann models.

Particles on demand for flows with strong discontinuities

Particles-on-demand formulation of kinetic theory [B. Dorschner, F. Bösch and I. V. Karlin, Phys. Rev. Lett. 121, 130602 (2018)0031-900710.1103/PhysRevLett.121.130602] is used to simulate a variety of compressible flows with strong discontinuities in density, pressure, and velocity. Two modifications are applied to the original formulation of the particles-on-demand method. First, a regularization by Grad's projection of particles populations is combined with the reference frame transformations in order to enhance stability and accuracy. Second, a finite-volume scheme is implemented which allows tight control of mass, momentum, and energy conservation. The proposed model is validated with an array of challenging one- and two-dimensional benchmarks of compressible flows, including hypersonic and near-vacuum situations, Richtmyer-Meshkov instability, double Mach reflection, and astrophysical jet. Excellent performance of the modified particles-on-demand method is demonstrated beyond the limitations of other lattice Boltzmann-like approaches to compressible flows.

Multiscale analysis of the particles on demand kinetic model

We present a thorough investigation of the particles on demand kinetic model. After a brief introduction of the method, an appropriate multiscale analysis is carried out to derive the hydrodynamic limit of the model. In this analysis, the effect of the time-space dependent comoving reference frames are taken into account. This could be regarded as a generalization of the conventional Chapman-Enskog analysis applied to the lattice Boltzmann models which feature global constant reference frames. Further simulations of target benchmarks provide numerical evidence confirming the theoretical predictions.

Unified lattice Boltzmann method with improved schemes for multiphase flow simulation: Application to droplet dynamics under realistic conditions

As a powerful mesoscale approach, the lattice Boltzmann method (LBM) has been widely used for the numerical study of complex multiphase flows. Recently, Luo et al. [Philos. Trans. R. Soc. A: Math. Phys. Eng. Sci. 379, 20200397 (2021)10.1098/rsta.2020.0397] proposed a unified lattice Boltzmann method (ULBM) to integrate the widely used lattice Boltzmann collision operators into a unified framework. In this study, we incorporate additional features into this ULBM in order to simulate multiphase flow under realistic conditions. A nonorthogonal moment set [Fei et al., Phys. Rev. E 97, 053309 (2018)10.1103/PhysRevE.97.053309] and the entropic-multi-relaxation-time (KBC) lattice Boltzmann model are used to construct the collision operator. An extended combined pseudopotential model is proposed to realize multiphase flow simulation at high-density ratio with tunable surface tension over a wide range. The numerical results indicate that the improved ULBM can significantly decrease the spurious velocities and adjust the surface tension without appreciably changing the density ratio. The ULBM is validated through reproducing various droplet dynamics experiments, such as binary droplet collision and droplet impingement on superhydrophobic surfaces. Finally, the extended ULBM is applied to complex droplet dynamics, including droplet pancake bouncing and droplet splashing. The maximum Weber number and Reynolds number in the simulation reach 800 and 7200, respectively, at a density ratio of 1000. The study demonstrates the generality and versatility of ULBM for incorporating schemes to tackle challenging multiphase problems.

Lattice Boltzmann model for capillary interactions between particles at a liquid-vapor interface under gravity

A computational technique based on the lattice Boltzmann method (LBM) is developed to simulate the wettable particles adsorbed to a liquid-vapor interface under gravity. The proposed technique combines the improved smoothed-profile LBM for the treatment of moving solid particles in a fluid and the free-energy LBM for the description of a liquid-vapor system. Five benchmark two-dimensional problems are examined: (A) a stationary liquid drop in the vapor phase; a wettable particle adsorbed to a liquid-vapor interface in (B) the absence and (C) the presence of gravity; (D) two freely moving particles at a liquid-vapor interface in the presence of gravity (i.e., capillary flotation forces); and (E) two vertically constrained particles at a liquid-vapor interface (i.e., capillary immersion forces). The simulation results are in good quantitative agreement with theoretical estimations, demonstrating that the proposed technique can reproduce the capillary interactions between wettable particles at a liquid-vapor interface under gravity.

Validation and application of the lattice Boltzmann algorithm for a turbulent immiscible Rayleigh-Taylor system

We develop a multicomponent lattice Boltzmann (LB) model for the two-dimensional Rayleigh-Taylor turbulence with a Shan-Chen pseudopotential implemented on GPUs. In the immiscible case, this method is able to accurately overcome the inherent numerical complexity caused by the complicated structure of the interface that appears in the fully developed turbulent regime. The accuracy of the LB model is tested both for early and late stages of instability. For the developed turbulent motion, we analyse the balance between different terms describing variations of the kinetic and potential energies. Then we analyse the role of the interface in the energy balance and also the effects of the vorticity induced by the interface in the energy dissipation. Statistical properties are compared for miscible and immiscible flows. Our results can also be considered as a first validation step to extend the application of LB model to three-dimensional immiscible Rayleigh-Taylor turbulence. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

Alternative wetting boundary condition for the chemical-potential-based free-energy lattice Boltzmann model

The free-energy lattice Boltzmann (LB) method is a multiphase LB approach based on the thermodynamic theory. Compared with traditional free-energy LB models, which employ a nonideal thermodynamic pressure tensor, the chemical-potential-based free-energy LB model has attracted much attention in recent years as it avoids computing the thermodynamic pressure tensor and its divergence. In this paper, we propose an improved wetting boundary condition for the chemical-potential-based free-energy LB model. Different from the original wetting boundary condition in the literature, the improved wetting boundary condition utilizes a surface chemical potential that is compatible with the chemical potential of the fluid domain. Accordingly, the thermodynamic consistency of the chemical-potential-based free-energy LB model can be retained by the improved wetting boundary condition. Numerical simulations are performed for droplets resting on flat and cylindrical surfaces with different contact angles. The numerical results show that the improved wetting boundary condition yields more reasonable results and the maximum spurious velocities are found to be smaller by 2 ∼ 3 orders of magnitude than those produced by the original wetting boundary condition.

Multiscale semi-Lagrangian lattice Boltzmann method

We present a multi-scale lattice Boltzmann scheme, which adaptively refines particles' velocity space. Different velocity sets of lower and higher order are consistently and efficiently coupled, allowing us to use the higher-order model only when and where needed. This includes regions of high Mach or high Knudsen numbers. The coupling procedure of discrete velocity sets consists of either a projection of the higher-order populations onto the lower-order lattice or lifting of the lower-order populations to the higher-order velocity space. Both lifting and projection are local operations, which enable a flexible adaptive velocity set. The proposed scheme is formulated for both a static and an optimal, co-moving reference frame, in the spirit of the recently introduced Particles on Demand method. The multi-scale scheme is validated with an advection of an athermal vortex and in a jet flow setup. The performance of the proposed scheme is further investigated in the shock structure problem and a high-Knudsen-number Couette flow, typical examples of highly non-equilibrium flows in which the order of the velocity set plays a decisive role. The results demonstrate that the proposed multi-scale scheme can operate accurately, with flexibility in terms of the underlying models and with reduced computational requirements.

Modeling ternary fluids in contact with elastic membranes

We present a thermodynamically consistent model of a ternary fluid interacting with elastic membranes. Following a free-energy modeling approach for the fluid phases, we derive the governing equations for the dynamics of the ternary fluid flow and membranes. We also provide the numerical framework for simulating such fluid-structure interaction problems. It is based on the lattice Boltzmann method for the ternary fluid (Eulerian description) and a finite difference representation of the membrane (Lagrangian description). The ternary fluid and membrane solvers are coupled through the immersed boundary method. For validation purposes, we consider the relaxation dynamics of a two-dimensional elastic capsule placed at a fluid-fluid interface. The capsule shapes, resulting from the balance of surface tension and elastic forces, are compared with equilibrium numerical solutions obtained by surface evolver. Furthermore, the Galilean invariance of the proposed model is proven. The proposed approach is versatile, allowing for the simulation of a wide range of geometries. To demonstrate this, we address the problem of a capillary bridge formed between two deformable capsules.

Kinetic Simulations of Compressible Non-Ideal Fluids: From Supercritical Flows to Phase-Change and Exotic Behavior

We investigate a kinetic model for compressible non-ideal fluids. The model imposes the local thermodynamic pressure through a rescaling of the particle’s velocities, which accounts for both long- and short-range effects and hence full thermodynamic consistency. The model is fully Galilean invariant and treats mass, momentum, and energy as local conservation laws. The analysis and derivation of the hydrodynamic limit is followed by the assessment of accuracy and robustness through benchmark simulations ranging from the Joule–Thompson effect to a phase-change and high-speed flows. In particular, we show the direct simulation of the inversion line of a van der Waals gas followed by simulations of phase-change such as the one-dimensional evaporation of a saturated liquid, nucleate, and film boiling and eventually, we investigate the stability of a perturbed strong shock front in two different fluid mediums. In all of the cases, we find excellent agreement with the corresponding theoretical analysis and experimental correlations. We show that our model can operate in the entire phase diagram, including super- as well as sub-critical regimes and inherently captures phase-change phenomena.

Lattice Boltzmann Solver for Multiphase Flows: Application to High Weber and Reynolds Numbers

The lattice Boltzmann method, now widely used for a variety of applications, has also been extended to model multiphase flows through different formulations. While already applied to many different configurations in low Weber and Reynolds number regimes, applications to higher Weber/Reynolds numbers or larger density/viscosity ratios are still the topic of active research. In this study, through a combination of a decoupled phase-field formulation—the conservative Allen–Cahn equation—and a cumulant-based collision operator for a low-Mach pressure-based flow solver, we present an algorithm that can be used for higher Reynolds/Weber numbers. The algorithm was validated through a variety of test cases, starting with the Rayleigh–Taylor instability in both 2D and 3D, followed by the impact of a droplet on a liquid sheet. In all simulations, the solver correctly captured the flow dynamics andmatched reference results very well. As the final test case, the solver was used to model droplet splashing on a thin liquid sheet in 3D with a density ratio of 1000 and kinematic viscosity ratio of 15, matching the water/air system at We = 8000 and Re = 1000. Results showed that the solver correctly captured the fingering instabilities at the crown rim and their subsequent breakup, in agreement with experimental and numerical observations reported in the literature.

Achieving thermodynamic consistency in a class of free-energy multiphase lattice Boltzmann models

The free-energy lattice Boltzmann (LB) model is one of the major multiphase models in the LB community. The present study is focused on a class of free-energy LB models in which the divergence of thermodynamic pressure tensor or its equivalent form expressed by the chemical potential is incorporated into the LB equation via a forcing term. Although this class of free-energy LB models may be thermodynamically consistent at the continuum level, it suffers from thermodynamic inconsistency at the discrete lattice level owing to numerical errors [Guo et al., Phys. Rev. E 83, 036707 (2010)10.1103/PhysRevE.83.036707]. The numerical error term mainly includes two parts: one comes from the discrete gradient operator and the other can be identified in a high-order Chapman-Enskog analysis. In this paper, we propose an improved scheme to eliminate the thermodynamic inconsistency of the aforementioned class of free-energy LB models. The improved scheme is constructed by modifying the equation of state of the standard LB equation, through which the discretization of ∇(ρc_{s}^{2}) is no longer involved in the force calculation and then the numerical errors can be significantly reduced. Numerical simulations are subsequently performed to validate the proposed scheme. The numerical results show that the improved scheme is capable of eliminating the thermodynamic inconsistency and can significantly reduce the spurious currents in comparison with the standard forcing-based free-energy LB model.

Thermokinetic lattice Boltzmann model of nonideal fluids

We present a kinetic model for nonideal fluids, where the local thermodynamic pressure is imposed through appropriate rescaling of the particle's velocities, accounting for both long- and short-range effects and hence full thermodynamic consistency. The model features full Galilean invariance together with mass, momentum, and energy conservation and enables simulations ranging from subcritical to supercritical flows, which is illustrated on various benchmark flows such as anomalous shock waves or shock droplet interaction.

Compressibility in lattice Boltzmann on standard stencils: effects of deviation from reference temperature

With growing interest in the simulation of compressible flows using the lattice Boltzmann (LB) method, a number of different approaches have been developed. These methods can be classified as pertaining to one of two major categories: (i) solvers relying on high-order stencils recovering the Navier-Stokes-Fourier equations, and (ii) approaches relying on classical first-neighbour stencils for the compressible Navier-Stokes equations coupled to an additional (LB-based or classical) solver for the energy balance equation. In most cases, the latter relies on a thermal Hermite expansion of the continuous equilibrium distribution function (EDF) to allow for compressibility. Even though recovering the correct equation of state at the Euler level, it has been observed that deviations of local flow temperature from the reference can result in instabilities and/or over-dissipation. The aim of the present study is to evaluate the stability domain of different EDFs, different collision models, with and without the correction terms for the third-order moments. The study is first based on a linear von Neumann analysis. The correction term for the space- and time-discretized equations is derived via a Chapman-Enskog analysis and further corroborated through spectral dispersion-dissipation curves. Finally, a number of numerical simulations are performed to illustrate the proposed theoretical study. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.

Chemical-potential multiphase lattice Boltzmann method with superlarge density ratios

The liquid-gas density ratio is a key property of multiphase flow methods to model real fluid systems. Here, a chemical-potential multiphase lattice Boltzmann method is constructed to realize extremely large density ratios. The simulations show that the method reaches very low temperatures, at which the liquid-gas density ratio is more than 10^{14}, while the thermodynamic consistency is still preserved. Decoupling the mesh space from the momentum space through a proportional coefficient, a smaller mesh step provides denser lattice nodes to exactly describe the transition region and the resulting dimensional transformation has no loss of accuracy. A compact finite-difference method is applied to calculate the discrete derivatives in the mesh space with high-order accuracy. These enhance the computational accuracy of the nonideal force and suppress the spurious currents to a very low level, even if the density ratio is up to tens of thousands. The simulation of drop splashing verifies that the present model is Galilean invariant for the dynamic flow field. An upper limit of the chemical potential is used to reduce the influence of nonphysical factors and improve the stability.

Lattice Boltzmann method for simulation of wettable particles at a fluid-fluid interface under gravity

A computational technique was developed to simulate wettable particles trapped at a fluid-fluid interface under gravity. The proposed technique combines the improved smoothed profile-lattice Boltzmann method (iSP-LBM) for the treatment of moving solid-fluid boundaries and the free-energy LBM for the description of isodensity immiscible two-phase flows. We considered five benchmark problems in two-dimensional systems, including a stationary drop, a wettable particle trapped at a fluid-fluid interface in the absence or presence of gravity, two freely moving particles at a fluid-fluid interface in the presence of gravity (i.e., capillary floatation forces), and two vertically constrained particles at a fluid-fluid interface (i.e., capillary immersion forces). The simulation results agreed well with theoretical estimations, demonstrating the efficacy of the proposed technique.

Semi-Lagrangian lattice Boltzmann model for compressible flows on unstructured meshes

Compressible lattice Boltzmann model on standard lattices [M. H. Saadat, F. Bösch, and I. V. Karlin, Phys. Rev. E 99, 013306 (2019).2470-004510.1103/PhysRevE.99.013306] is extended to deal with complex flows on unstructured grid. Semi-Lagrangian propagation [A. Krämer et al., Phys. Rev. E 95, 023305 (2017).2470-004510.1103/PhysRevE.95.023305] is performed on an unstructured second-order accurate finite-element mesh and a consistent wall boundary condition is implemented which makes it possible to simulate compressible flows over complex geometries. The model is validated through simulation of Sod shock tube, subsonic and supersonic flow over NACA0012 airfoil and shock-vortex interaction in Schardin's problem. Numerical results demonstrate that the present model on standard lattices is able to simulate compressible flows involving shock waves on unstructured meshes with good accuracy and without using any artificial dissipation or limiter.

Lattice Boltzmann model for weakly compressible flows

We present an energy conserving lattice Boltzmann model based on a crystallographic lattice for simulation of weakly compressible flows. The theoretical requirements and the methodology to construct such a model are discussed. We demonstrate that the model recovers the isentropic sound speed in addition to the effects of viscous heating and heat flux dynamics. Several test cases for acoustics and thermal and thermoacoustic flows are simulated to show the accuracy of the proposed model.

Pseudoentropic derivation of the regularized lattice Boltzmann method

The lattice Boltzmann method (LBM) facilitates efficient simulations of fluid turbulence based on advection and collision of local particle distribution functions. To ensure stable simulations on underresolved grids, the collision operator must prevent drastic deviations from local equilibrium. This can be achieved by various methods, such as the multirelaxation time, entropic, quasiequilibrium, regularized, and cumulant schemes. Complementing a part of a unified theoretical framework of these schemes, the present work presents a derivation of the regularized lattice Boltzmann method (RLBM), which follows a recently introduced entropic multirelaxation time LBM by Karlin, Bösch, and Chikatamarla (KBC). It is shown that both methods can be derived by locally maximizing a quadratic Taylor expansion of the entropy function. While KBC expands around the local equilibrium distribution, the RLBM is recovered by expanding entropy around a global equilibrium. Numerical tests were performed to elucidate the role of pseudoentropy maximization in these models. Simulations of a two-dimensional shear layer show that the RLBM successfully reproduces the largest eddies even on a 16×16 grid, while the conventional LBM becomes unstable for grid resolutions of 128×128 and lower. The RLBM suppresses spurious vortices more effectively than KBC. In contrast, simulations of the three-dimensional Taylor-Green and Kida vortices show that KBC performs better in resolving small scale vortices, outperforming the RLBM by a factor of 1.8 in terms of the effective Reynolds number.

Wetting boundaries for a ternary high-density-ratio lattice Boltzmann method

We extend a recently proposed ternary free-energy lattice Boltzmann model with high density contrast [Phys. Rev. Lett. 120, 234501 (2018)PRLTAO0031-900710.1103/PhysRevLett.120.234501] by incorporating wetting boundaries at solid walls. The approaches are based on forcing and geometric schemes, with implementations optimized for ternary (and, more generally, higher-order multicomponent) models. Advantages and disadvantages of each method are addressed by performing both static and dynamic tests, including the capillary filling dynamics of a liquid displacing the gas phase and the self-propelled motion of a train of drops. Furthermore, we measure dynamic angles and show that the slip length critically depends on the equilibrium value of the contact angles and whether it belongs to liquid-liquid or liquid-gas interfaces. These results validate the model capabilities of simulating complex ternary fluid dynamic problems near solid boundaries, for example, drop impact solid substrates covered by a lubricant layer.

Multiple-time-scaling lattice Boltzmann method for the convection diffusion equation

A multiple-time-scaling (MTS) strategy that decouples the time discretization in different domains and enables flexible time-step coarsening, refinement, and stretching in the lattice Boltzmann method (LBM) for the convection diffusion equation is developed. The key in the multiple-time scaling is the satisfaction of physical interface conditions without nested iterations in each time step by implementing appropriate interface schemes for the distribution functions in the LBM. The applicability and second-order accuracy of the MTS-LBM approach is demonstrated with two numerical tests. Our approach greatly expands the avenue and expedites the progress of applying the LBM for modeling complex flows and transport phenomena involving multiphases and multicomponents with large property ratios.

Interaction pressure tensor on high-order lattice Boltzmann models for nonideal fluids

In this work we address the application of pseudopotentials directly on high-order lattice Boltzmann models. We derive a general expression for the pressure tensor on high-order lattices considering all nonideal interactions, including intra- and intermolecular interactions, following the discrete lattice theory introduced by X. Shan [Phys. Rev. E 77, 066702 (2008)PLEEE81539-375510.1103/PhysRevE.77.066702]. From the derived expression, a generalized continuum approximation, truncated at fourth-order isotropy, is obtained that is readily applicable to high-order lattices. With this, we demonstrate that high-order lattice models with pseudopotentials can satisfy thermodynamic consistency. The derived generalized expression and continuum approximation are validated for the case of a flat interface and compared against the standard definition available from the literature. The generalized expression is also shown to accurately reproduce the Laplace experiment for a variety of high-order lattice structures. This work sets the preliminary steps towards the application of high-order lattice models for simulating nonideal fluid mixtures.

Tricoupled hybrid lattice Boltzmann model for nonisothermal drying of colloidal suspensions in micropore structures

A tricoupled hybrid lattice Boltzmann model (LBM) is developed to simulate colloidal liquid evaporation and colloidal particle deposition during the nonisothermal drying of colloidal suspensions in micropore structures. An entropic multiple-relaxation-time multirange pseudopotential two-phase LBM for isothermal interfacial flow is first coupled to an extended temperature equation for simulating nonisothermal liquid drying. Then the coupled model is further coupled with a modified convection diffusion equation to consider the nonisothermal drying of colloidal suspensions. Two drying examples are considered. First, drying of colloidal suspensions in a two-pillar micropore structure is simulated in two dimensions (2D), and the final configuration of colloidal particles is compared with the experimental one. Good agreement is observed. Second, at the temperature of 343.15 K (70^{∘}C), drying of colloidal suspensions in a complex spiral-shaped micropore structure containing 220 pillars is simulated (also in 2D). The drying pattern follows the designed spiral shape due to capillary pumping, i.e., transport of the liquid from larger pores to smaller ones by capillary pressure difference. Since the colloidal particles are passively carried with liquid, they accumulate at the small menisci as drying proceeds. As liquid evaporates at the small menisci, colloidal particles are deposited, eventually forming solid structures between the pillars (primarily), and at the base of the pillars (secondarily). As a result, the particle deposition conforms to the spiral route. Qualitatively, the simulated liquid and particle configurations agree well with the experimental ones during the entire drying process. Quantitatively, the model demonstrates that the evaporation rate and the particle accumulation rate slowly decrease during drying, similar to what is seen in the experimental results, which is due to the reduction of the liquid-vapor interfacial area. In conclusion, the hybrid model shows the capability and accuracy for simulating nonisothermal drying of colloidal suspensions in a complex micropore structure both qualitatively and quantitatively, as it includes all the required physics and captures all the complex features observed experimentally. Such a tricoupled LBM has a high potential to become an efficient numerical tool for further investigation of real and complex engineering problems incorporating drying of colloidal suspensions in porous media.

Entropy production in thermal phase separation: a kinetic-theory approach

Entropy production during the process of thermal phase-separation of multiphase flows is investigated by means of a discrete Boltzmann kinetic model. The entropy production rate is found to increase during the spinodal decomposition stage and to decrease during the domain growth stage, attaining its maximum at the crossover between the two. Such behaviour provides a natural criterion to identify and discriminate between the two regimes. Furthermore, the effects of heat conductivity, viscosity and surface tension on the entropy production rate are investigated by systematically probing the interplay between non-equilibrium energy and momentum fluxes. It is found that the entropy production rate due to energy fluxes is an increasing function of the Prandtl number, while the momentum fluxes exhibit an opposite trend. On the other hand, both contributions show an increasing trend with surface tension. The present analysis inscribes within the general framework of non-equilibrium thermodynamics and consequently it is expected to be relevant to a broad class of soft-flowing systems far from mechanical and thermal equilibrium.

Mass conservative lattice Boltzmann scheme for a three-dimensional diffuse interface model with Peng-Robinson equation of state

Peng-Robinson (P-R) equation of state (EOS) has been widely used in the petroleum industry for hydrocarbon fluids. In this work, a three-dimensional diffuse interface model with P-R EOS for two-phase fluid system is solved by the lattice Boltzmann (LB) method. In this diffuse interface model, an Allen-Cahn (A-C) type phase equation with strong nonlinear source term is derived. Using the multiscale Chapman-Enskog analysis, the A-C type phase equation can be recovered from the proposed LB method. Besides, a Lagrange multiplier is introduced based on the mesoscopic character of the LB scheme so that total mass of the hydrocarbon system is preserved. Three-dimensional numerical simulations of realistic hydrocarbon components, such as isobutane and propane, are implemented to illustrate the effectiveness of the proposed mass conservative LB scheme. Numerical results reach a better agreement with laboratory data compared to previous results of two-dimensional numerical simulations.

Ternary Free-Energy Entropic Lattice Boltzmann Model with a High Density Ratio

A thermodynamically consistent free energy model for fluid flows comprised of one gas and two liquid components is presented and implemented using the entropic lattice Boltzmann scheme. The model allows a high density ratio, up to the order of O(10^{3}), between the liquid and gas phases, and a broad range of surface tension ratios, covering partial wetting states where Neumann triangles are formed, and full wetting states where complete encapsulation of one of the fluid components is observed. We further demonstrate that we can capture the bouncing, adhesive, and insertive regimes for the binary collisions between immiscible droplets suspended in air. Our approach opens up a vast range of multiphase flow applications involving one gas and several liquid components.

Dynamics of Contact Line Pinning and Depinning of Droplets Evaporating on Microribs

The contact line dynamics of evaporating droplets deposited on a set of parallel microribs is analyzed with the use of a recently developed entropic lattice Boltzmann model for two-phase flow. Upon deposition, part of the droplet penetrates into the space between ribs because of capillary action, whereas the remaining liquid of the droplet remains pinned on top of the microribs. In the first stage, evaporation continues until the droplet undergoes a series of pinning-depinning events, showing alternatively the constant contact radius and constant contact angle modes. While the droplet is pinned, evaporation results in a contact angle reduction, whereas the contact radius remains constant. At a critical contact angle, the contact line depins, the contact radius reduces, and the droplet rearranges to a larger apparent contact angle. This pinning-depinning behavior goes on until the liquid above the microribs is evaporated. By computing the Gibbs free energy taking into account the interfacial energy, pressure terms, and viscous dissipation due to drop internal flow, we found that the mechanism that causes the unpinning of the contact line results from an excess in Gibbs free energy. The spacing distance and the rib height play an important role in controlling the pinning-depinning cycling, the critical contact angle, and the excess Gibbs free energy. However, we found that neither the critical contact angle nor the maximum excess Gibbs free energy depends on the rib width. We show that the different terms, that is, pressure term, viscous dissipation, and interfacial energy, contributing to the excess Gibbs free energy, can be varied differently by varying different geometrical properties of the microribs. It is demonstrated that, by varying the spacing distance between the ribs, the energy barrier is controlled by the interfacial energy while the contribution of the viscous dissipation is dominant if either rib height or width is changed. Main finding of this is study is that, for microrib patterned surfaces, the energy barrier required for the contact line to depin can be enlarged by increasing the spacing or the rib height, which can be important for practical applications.

Water ring-bouncing on repellent singularities

Texturing a flat superhydrophobic substrate with point-like superhydrophobic macrotextures of the same repellency makes impacting water droplets take off as rings, which leads to shorter bouncing times than on a flat substrate. We investigate the contact time reduction on such elementary macrotextures through experiment and simulations. We understand the observations by decomposing the impacting drop reshaped by the defect into sub-units (or blobs) whose size is fixed by the liquid ring width. We test the blob picture by looking at the reduction of contact time for off-centered impacts and for impacts in grooves that produce liquid ribbons where the blob size is fixed by the width of the channel.

Numerical simulation of three-component multiphase flows at high density and viscosity ratios using lattice Boltzmann methods

In this paper, we propose a multiphase lattice Boltzmann model for numerical simulation of ternary flows at high density and viscosity ratios free from spurious velocities. The proposed scheme, which is based on the phase-field modeling, employs the Cahn-Hilliard theory to track the interfaces among three different fluid components. Several benchmarks, such as the spreading of a liquid lens, binary droplets, and head-on collision of two droplets in binary- and ternary-fluid systems, are conducted to assess the reliability and accuracy of the model. The proposed model can successfully simulate both partial and total spreadings while reducing the parasitic currents to the machine precision.

Entropic lattice Boltzmann model for charged leaky dielectric multiphase fluids in electrified jets

We present a lattice Boltzmann model for charged leaky dielectric multiphase fluids in the context of electrified jet simulations, which are of interest for a number of production technologies including electrospinning. The role of nonlinear rheology on the dynamics of electrified jets is considered by exploiting the Carreau model for pseudoplastic fluids. We report exploratory simulations of charged droplets at rest and under a constant electric field, and we provide results for charged jet formation under electrospinning conditions.

Fluid-structure interaction with the entropic lattice Boltzmann method

We propose a fluid-structure interaction (FSI) scheme using the entropic multi-relaxation time lattice Boltzmann (KBC) model for the fluid domain in combination with a nonlinear finite element solver for the structural part. We show the validity of the proposed scheme for various challenging setups by comparison to literature data. Beyond validation, we extend the KBC model to multiphase flows and couple it with a finite element method (FEM) solver. Robustness and viability of the entropic multi-relaxation time model for complex FSI applications is shown by simulations of droplet impact on elastic superhydrophobic surfaces.

Essentially Entropic Lattice Boltzmann Model

The entropic lattice Boltzmann model (ELBM), a discrete space-time kinetic theory for hydrodynamics, ensures nonlinear stability via the discrete time version of the second law of thermodynamics (the H theorem). Compliance with the H theorem is numerically enforced in this methodology and involves a search for the maximal discrete path length corresponding to the zero dissipation state by iteratively solving a nonlinear equation. We demonstrate that an exact solution for the path length can be obtained by assuming a natural criterion of negative entropy change, thereby reducing the problem to solving an inequality. This inequality is solved by creating a new framework for construction of Padé approximants via quadrature on appropriate convex function. This exact solution also resolves the issue of indeterminacy in case of nonexistence of the entropic involution step. Since our formulation is devoid of complex mathematical library functions, the computational cost is drastically reduced. To illustrate this, we have simulated a model setup of flow over the NACA-0012 airfoil at a Reynolds number of 2.88×10^{6}.