Pattern formation in liquid-vapor systems under periodic potential and shear

Growth regimes in three-dimensional phase separation of liquid-vapor systems

The liquid-vapor phase separation is investigated via lattice Boltzmann simulations in three dimensions. After expressing length and time scales in reduced physical units, we combined data from several large simulations (on 512^{3} nodes) with different values of viscosity, surface tension, and temperature, to obtain a single curve of rescaled length l[over ̂] as a function of rescaled time t[over ̂]. We find evidence of the existence of kinetic and inertial regimes with growth exponents α_{d}=1/2 and α_{i}=2/3 over several time decades, with a crossover from α_{d} to α_{i} at t[over ̂]≃1. This allows us to rule out the existence of a viscous regime with α_{v}=1 in three-dimensional liquid-vapor isothermal phase separation, differently from what happens in binary fluid mixtures. An in-depth analysis of the kinetics of the phase separation process, as well as a characterization of the morphology and the flow properties, are further presented in order to provide clues into the dynamics of the phase-separation process.

Simulation of collapsing cavitation bubbles in various liquids by lattice Boltzmann model coupled with the Redlich-Kwong-Soave equation of state

A computational technique based on the pseudo-potential multiphase lattice Boltzmann method (LBM) is employed to investigate the collapse dynamics of cavitation bubbles of various liquids in the vicinity of the solid surface with different wettability conditions. The Redlich-Kwong-Soave equation of state (EoS) that includes an acentric factor is incorporated to consider the physical properties of water (H_{2}O), liquid nitrogen (LN_{2}), and liquid hydrogen (LH_{2}) in the present simulations. Accuracy and performance of the present multiphase LBM are examined by simulation of the homogenous and heterogeneous cavitation phenomena. The good agreement of the results obtained based on the present solution algorithm in comparison with the available data confirms the validity and capability of the multiphase LBM employed. Then, the cavitation bubble collapse near the solid wall is studied by considering the H_{2}O, LN_{2}, and LH_{2} fluids, and the wettability effect of the surface on the collapse dynamics is investigated. The obtained results demonstrate that the collapse phenomenon for the H_{2}O is more aggressive than that of the LH_{2} and LN_{2}. The cavitation bubble of the water has a shorter collapse time with an intense liquid jet, while the collapse process in the LN_{2} takes a longer time due to the larger radius of its bubble at the rebound. Also, this study demonstrates that the increment of the hydrophobicity of the wall causes less energy absorption by the solid surface from the liquid phase around the bubble that leads to form a liquid jet with higher kinetic energy. Therefore, the bubble collapse process occurs more quickly for hydrophobic surfaces, regardless of the fluids considered. The present study shows that the pseudopotential LBM with incorporating an appropriate EoS and a robust forcing scheme is an efficient numerical technique for simulation of the dynamics of the cavitation bubble collapse in different fluids.

Lattice Boltzmann methods and active fluids

We review the state of the art of active fluids with particular attention to hydrodynamic continuous models and to the use of Lattice Boltzmann Methods (LBM) in this field. We present the thermodynamics of active fluids, in terms of liquid crystals modelling adapted to describe large-scale organization of active systems, as well as other effective phenomenological models. We discuss how LBM can be implemented to solve the hydrodynamics of active matter, starting from the case of a simple fluid, for which we explicitly recover the continuous equations by means of Chapman-Enskog expansion. Going beyond this simple case, we summarize how LBM can be used to treat complex and active fluids. We then review recent developments concerning some relevant topics in active matter that have been studied by means of LBM: spontaneous flow, self-propelled droplets, active emulsions, rheology, active turbulence, and active colloids.

Corner-transport-upwind lattice Boltzmann model for bubble cavitation

Aiming to study the bubble cavitation problem in quiescent and sheared liquids, a third-order isothermal lattice Boltzmann model that describes a two-dimensional (2D) fluid obeying the van der Waals equation of state, is introduced. The evolution equations for the distribution functions in this off-lattice model with 16 velocities are solved using the corner-transport-upwind (CTU) numerical scheme on large square lattices (up to 6144×6144 nodes). The numerical viscosity and the regularization of the model are discussed for first- and second-order CTU schemes finding that the latter choice allows to obtain a very accurate phase diagram of a nonideal fluid. In a quiescent liquid, the present model allows us to recover the solution of the 2D Rayleigh-Plesset equation for a growing vapor bubble. In a sheared liquid, we investigated the evolution of the total bubble area, the bubble deformation, and the bubble tilt angle, for various values of the shear rate. A linear relation between the dimensionless deformation coefficient D and the capillary number Ca is found at small Ca but with a different factor than in equilibrium liquids. A nonlinear regime is observed for Ca≳0.2.

Lattice Boltzmann study of chemically-driven self-propelled droplets

We numerically study the behavior of self-propelled liquid droplets whose motion is triggered by a Marangoni-like flow. This latter is generated by variations of surfactant concentration which affect the droplet surface tension promoting its motion. In the present paper a model for droplets with a third amphiphilic component is adopted. The dynamics is described by Navier-Stokes and convection-diffusion equations, solved by the lattice Boltzmann method coupled with finite-difference schemes. We focus on two cases. First, the study of self-propulsion of an isolated droplet is carried on and, then, the interaction of two self-propelled droplets is investigated. In both cases, when the surfactant migrates towards the interface, a quadrupolar vortex of the velocity field forms inside the droplet and causes the motion. A weaker dipolar field emerges instead when the surfactant is mainly diluted in the bulk. The dynamics of two interacting droplets is more complex and strongly depends on their reciprocal distance. If, in a head-on collision, droplets are close enough, the velocity field initially attracts them until a motionless steady state is achieved. If the droplets are vertically shifted, the hydrodynamic field leads to an initial reciprocal attraction followed by a scattering along opposite directions. This hydrodynamic interaction acts on a separation of some droplet radii otherwise it becomes negligible and droplets motion is only driven by the Marangoni effect. Finally, if one of the droplets is passive, this latter is generally advected by the fluid flow generated by the active one.

Discrete Boltzmann modeling of multiphase flows: hydrodynamic and thermodynamic non-equilibrium effects

A discrete Boltzmann model (DBM) is developed to investigate the hydrodynamic and thermodynamic non-equilibrium (TNE) effects in phase separation processes. The interparticle force drives changes and the gradient force, induced by gradients of macroscopic quantities, opposes them. In this paper, we investigate the interplay between them by providing a detailed inspection of various non-equilibrium observables. Based on the TNE features, we define TNE strength which roughly estimates the deviation amplitude from the thermodynamic equilibrium. The time evolution of the TNE intensity provides a convenient and efficient physical criterion to discriminate the stages of the spinodal decomposition and domain growth. Via the DBM simulation and this criterion, we quantitatively study the effects of latent heat and surface tension on phase separation. It is found that the TNE strength attains its maximum at the end of the spinodal decomposition stage, and it decreases when the latent heat increases from zero. The surface tension effects are threefold, prolong the duration of the spinodal decomposition stage, decrease the maximum TNE intensity, and accelerate the speed of the domain growth stage.