Mesoscopic lattice boltzmann modeling of flowing soft systems

Effect of droplet deformability on shear thinning in a cylindrical channel

Droplets suspended in fluids flowing through microchannels are often encountered in different contexts and scales, from oil extraction down to microfluidics. They are usually flexible and deform as a product of the interplay between flexibility, hydrodynamics, and interaction with confining walls. Deformability adds distinct characteristics to the nature of the flow of these droplets. We simulate deformable droplets suspended in a fluid at a high volume fraction flowing through a cylindrical wetting channel. We find a discontinuous shear thinning transition, which depends on the droplet deformability. The capillary number is the main dimensionless parameter that controls the transition. Previous results have focused on two-dimensional configurations. Here we show that, in three dimensions, even the velocity profile is different. To perform this study, we improve and extend to three dimensions a multicomponent lattice Boltzmann method which prevents the coalescence between the droplets.

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.

Mesoscale perspective on the Tolman length

We demonstrate that the multiphase Shan-Chen lattice Boltzmann method (LBM) yields a curvature dependent surface tension σ as computed from three-dimensional hydrostatic droplets and bubbles simulations. Such curvature dependence is routinely characterized, at first order, by the so-called Tolman length δ. LBM allows one to precisely compute σ at the surface of tension R_{s} and determine the Tolman length from the coefficient of the first order correction. The corresponding values of δ display universality for different equations of state, following a power-law scaling near the critical temperature. The Tolman length has been studied so far mainly via computationally demanding Molecular Dynamics simulations or by means of Density Functional Theory approaches playing a pivotal role in extending Classical Nucleation Theory. The present results open a hydrodynamic-compliant mesoscale arena, in which the fundamental role of the Tolman length, alongside real-world applications to cavitation phenomena, can be effectively tackled. All the results can be independently reproduced through the "idea.deploy" framework.

A unified lattice Boltzmann model and application to multiphase flows

In this work, we develop a unified lattice Boltzmann model (ULBM) framework that can seamlessly integrate the widely used lattice Boltzmann collision operators, including the Bhatnagar-Gross-Krook or single-relation-time, multiple-relaxation-time, central-moment or cascaded lattice Boltzmann method and multiple entropic operators (KBC). Such a framework clarifies the relations among the existing collision operators and greatly facilitates model comparison and development as well as coding. Importantly, any LB model or treatment constructed for a specific collision operator could be easily adopted by other operators. We demonstrate the flexibility and power of the ULBM framework through three multiphase flow problems: the rheology of an emulsion, splashing of a droplet on a liquid film and dynamics of pool boiling. Further exploration of ULBM for a wide variety of phenomena would be both realistic and beneficial, making the LBM more accessible to non-specialists. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

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 simulations capture the multiscale physics of soft flowing crystals

The study of the underlying physics of soft flowing materials depends heavily on numerical simulations, due to the complex structure of the governing equations reflecting the competition of concurrent mechanisms acting at widely disparate scales in space and time. A full-scale computational modelling remains a formidable challenge since it amounts to simultaneously handling six or more spatial decades in space and twice as many in time. Coarse-grained methods often provide a viable strategy to significantly mitigate this issue, through the implementation of mesoscale supramolecular forces designed to capture the essential physics at a fraction of the computational cost of a full-detail description. Here, we review some recent advances in the design of a lattice Boltzmann mesoscale approach for soft flowing materials, inclusive of near-contact interactions (NCIs) between dynamic interfaces, as they occur in high packing-fraction soft flowing crystals. The method proves capable of capturing several aspects of the rheology of soft flowing crystals, namely, (i) a 3/2 power-law dependence of the dispersed phase flow rate on the applied pressure gradient, (ii) the structural transition between an ex-two and ex-one (bamboo) configurations with the associated drop of the flow rate, (iii) the onset of interfacial waves once NCI is sufficiently intense. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.

Discrete fluidization of dense monodisperse emulsions in neutral wetting microchannels

The rheology of pressure-driven flows of two-dimensional dense monodisperse emulsions in neutral wetting microchannels is investigated by means of mesoscopic lattice Boltzmann simulations, capable of handling large collections of droplets, in the order of several hundreds. The simulations reveal that the fluidization of the emulsion proceeds through a sequence of discrete steps, characterized by yielding events whereby layers of droplets start rolling over each other, thus leading to sudden drops of the relative effective viscosity. It is shown that such discrete fluidization is robust against loss of confinement, namely it persists also in the regime of small ratios of the droplet diameter over the microchannel width. We also develop a simple phenomenological model which predicts a linear relation between the relative effective viscosity of the emulsion and the product of the confinement parameter (global size of the device over droplet radius) and the viscosity ratio between the disperse and continuous phases. The model shows excellent agreement with the numerical simulations. The present work offers new insights to enable the design of microfluidic scaffolds for tissue engineering applications and paves the way to detailed rheological studies of soft-glassy materials in complex geometries.

Particles on Demand for Kinetic Theory

A novel formulation of fluid dynamics as a kinetic theory with tailored, on-demand constructed particles removes restrictions on flow speed and temperature as compared to its predecessors, the lattice Boltzmann methods and their modifications. In the new kinetic theory, discrete particles are determined by a rigorous limit process which avoids ad hoc assumptions about their velocities. Classical benchmarks for incompressible and compressible flows demonstrate that the proposed discrete-particles kinetic theory opens up an unprecedented wide domain of applications for computational fluid dynamics.

Lattice Kinetic Theory in a Comoving Galilean Reference Frame

We prove that the fully discrete lattice Boltzmann method is invariant with respect to Galilean transformation. Based on this finding, a novel class of shifted lattices is proposed which dramatically increases the operating range of lattice Boltzmann simulations, in particular, for gas dynamics applications. A simulation of vortex-shock interaction is used to demonstrate the accuracy and efficiency of the proposed lattices. With one single algorithm it is now possible to simulate a broad range of applications, from low Mach number flows to transonic and supersonic flow regimes.

Entropic multirelaxation lattice Boltzmann models for turbulent flows

We present three-dimensional realizations of a class of lattice Boltzmann models introduced recently by the authors [I. V. Karlin, F. Bösch, and S. S. Chikatamarla, Phys. Rev. E 90, 031302(R) (2014)] and review the role of the entropic stabilizer. Both coarse- and fine-grid simulations are addressed for the Kida vortex flow benchmark. We show that the outstanding numerical stability and performance is independent of a particular choice of the moment representation for high-Reynolds-number flows. We report accurate results for low-order moments for homogeneous isotropic decaying turbulence and second-order grid convergence for most assessed statistical quantities. It is demonstrated that all the three-dimensional lattice Boltzmann realizations considered herein converge to the familiar lattice Bhatnagar-Gross-Krook model when the resolution is increased. Moreover, thanks to the dynamic nature of the entropic stabilizer, the present model features less compressibility effects and maintains correct energy and enstrophy dissipation. The explicit and efficient nature of the present lattice Boltzmann method renders it a promising candidate for both engineering and scientific purposes for highly turbulent flows.

Modeling the flow in diffuse interface methods of solidification

Fluid dynamical equations in the presence of a diffuse solid-liquid interface are investigated via a volume averaging approach. The resulting equations exhibit the same structure as the standard Navier-Stokes equation for a Newtonian fluid with a constant viscosity, the effect of the solid phase fraction appearing in the drag force only. This considerably simplifies the use of the lattice Boltzmann method as a fluid dynamics solver in solidification simulations. Galilean invariance is also satisfied within this approach. Further, we investigate deviations between the diffuse and sharp interface flow profiles via both quasiexact numerical integration and lattice Boltzmann simulations. It emerges from these studies that the freedom in choosing the solid-liquid coupling parameter h provides a flexible way of optimizing the diffuse interface-flow simulations. Once h is adapted for a given spatial resolution, the simulated flow profiles reach an accuracy comparable to quasiexact numerical simulations.

Internal dynamics and activated processes in soft-glassy materials

Plastic rearrangements play a crucial role in the characterization of soft-glassy materials, such as emulsions and foams. Based on numerical simulations of soft-glassy systems, we study the dynamics of plastic rearrangements at the hydrodynamic scales where thermal fluctuations can be neglected. Plastic rearrangements require an energy input, which can be either provided by external sources, or made available through time evolution in the coarsening dynamics, in which the total interfacial area decreases as a consequence of the slow evolution of the dispersed phase from smaller to large droplets/bubbles. We first demonstrate that our hydrodynamic model can quantitatively reproduce such coarsening dynamics. Then, considering periodically oscillating strains, we characterize the number of plastic rearrangements as a function of the external energy-supply, and show that they can be regarded as activated processes induced by a suitable "noise" effect. Here we use the word noise in a broad sense, referring to the internal non-equilibrium dynamics triggered by spatial random heterogeneities and coarsening. Finally, by exploring the interplay between the internal characteristic time-scale of the coarsening dynamics and the external time-scale associated with the imposed oscillating strain, we show that the system exhibits the phenomenon of stochastic resonance, thereby providing further credit to the mechanical activation scenario.

Gibbs' principle for the lattice-kinetic theory of fluid dynamics

Gibbs' seminal prescription for constructing optimal states by maximizing the entropy under pertinent constraints is used to derive a lattice kinetic theory for the computation of high Reynolds number flows. The notion of modifying the viscosity to stabilize subgrid simulations is challenged in this kinetic framework. A lattice Boltzmann model for direct simulation of turbulent flows is presented without any need for tunable parameters and turbulent viscosity. Simulations at very high Reynolds numbers demonstrate a major extension of the operation range for fluid dynamics.

Direct evidence of plastic events and dynamic heterogeneities in soft-glasses

By using fluid-kinetic simulations of confined and concentrated emulsion droplets, we investigate the nature of space non-homogeneity in soft-glassy dynamics and provide quantitative measurements of the statistical features of plastic events in the proximity of the yield-stress threshold. Above the yield stress, our results show the existence of a finite stress correlation scale, which can be mapped directly onto the cooperativity scale, recently introduced in the literature to capture non-local effects in the soft-glassy dynamics. In this regime, the emergence of a separate boundary (wall) rheology with higher fluidity than the bulk is highlighted in terms of near-wall spontaneous segregation of plastic events. Near the yield stress, where the cooperativity scale cannot be estimated with sufficient accuracy, the system shows a clear increase of the stress correlation scale, whereas plastic events exhibit intermittent clustering in time, with no preferential spatial location. A quantitative measurement of the space-time correlation associated with the motion of the interface of the droplets is key to spot the elastic rigidity of the system.

Achieving tunable surface tension in the pseudopotential lattice Boltzmann modeling of multiphase flows

In this paper, we aim to address an important issue about the pseudopotential lattice Boltzmann (LB) model, which has attracted much attention as a mesoscopic model for simulating interfacial dynamics of complex fluids, but suffers from the problem that the surface tension cannot be tuned independently of the density ratio. In the literature, a multirange potential was devised to adjust the surface tension [Sbragaglia et al., Phys. Rev. E 75, 026702 (2007)]. However, it was recently found that the density ratio of the system will be changed when the multirange potential is employed to adjust the surface tension. An alternative approach is therefore proposed in the present work. The basic strategy is to add a source term to the LB equation so as to tune the surface tension of the pseudopotential LB model. The proposed approach can guarantee that the adjustment of the surface tension does not affect the mechanical stability condition of the pseudopotential LB model, and thus provides a separate control of the surface tension and the density ratio. Meanwhile, it still retains the mesoscopic feature and the computational simplicity of the pseudopotential LB model. Numerical simulations are carried out for stationary droplets, capillary waves, and droplet splashing on a thin liquid film. The numerical results demonstrate that the proposed approach is capable of achieving a tunable surface tension over a very wide range and can keep the density ratio unchanged when adjusting the surface tension.

Hydrodynamically driven colloidal assembly in dip coating

We study the hydrodynamics of dip coating from a suspension and report a mechanism for colloidal assembly and pattern formation on smooth substrates. Below a critical withdrawal speed where the coating film is thinner than the particle diameter, capillary forces induced by deformation of the free surface prevent the convective transport of single particles through the meniscus beneath the film. Capillary-induced forces are balanced by hydrodynamic drag only after a minimum number of particles assemble within the meniscus. The particle assembly can thus enter the thin film where it moves at nearly the withdrawal speed and rapidly separates from the next assembly. The interplay between hydrodynamic and capillary forces produces periodic and regular structures below a critical ratio Ca(2/3)/sqrt[Bo] < 0.7, where Ca and Bo are the capillary and Bond numbers, respectively. An analytical model and numerical simulations are presented for the case of two-dimensional flow with circular particles in suspension. The hydrodynamically driven assembly documented here is consistent with stripe pattern formations observed experimentally in dip coating.

Multiple-component lattice Boltzmann equation for fluid-filled vesicles in flow

We document the derivation and implementation of extensions to a two-dimensional, multicomponent lattice Boltzmann equation model, with Laplace law interfacial tension. The extended model behaves in such a way that the boundary between its immiscible drop and embedding fluid components can be shown to describe a vesicle of constant volume bounded by a membrane with conserved length, specified interface compressibility, bending rigidity, preferred curvature, and interfacial tension. We describe how to apply this result to several, independent vesicles. The extended scheme is completely Eulerian, and it represents a two-way coupled vesicle membrane and flow within a single framework. Unlike previous methods, our approach dispenses entirely with the need explicitly to track the membrane, or boundary, and makes no use whatsoever of computationally expensive and intricate interface tracking and remeshing. Validation data are presented, which demonstrate the utility of the method in the simulation of the flow of high volume fraction suspensions of deformable objects.

Mesoscopic model for microscale hydrodynamics and interfacial phenomena: slip, films, and contact-angle hysteresis

We present a model based on the lattice Boltzmann equation that is suitable for the simulation of dynamic wetting. The model is capable of exhibiting fundamental interfacial phenomena such as weak adsorption of fluid on the solid substrate and the presence of a thin surface film within which a disjoining pressure acts. Dynamics in this surface film, tightly coupled with hydrodynamics in the fluid bulk, determine macroscopic properties of primary interest: the hydrodynamic slip; the equilibrium contact angle; and the static and dynamic hysteresis of the contact angles. The pseudo-potentials employed for fluid-solid interactions are composed of a repulsive core and an attractive tail that can be independently adjusted. This enables effective modification of the functional form of the disjoining pressure so that one can vary the static and dynamic hysteresis on surfaces that exhibit the same equilibrium contact angle. The modeled fluid-solid interface is diffuse, represented by a wall probability function that ultimately controls the momentum exchange between solid and fluid phases. This approach allows us to effectively vary the slip length for a given wettability (i.e., a given static contact angle) of the solid substrate.

Heterogeneous diffuse interfaces: a new mechanism for arrested coarsening in binary mixtures. Heterogeneous diffuse interfaces

We discuss the dynamics of binary fluid mixtures in which surface tension density is allowed to become locally negative within the interface, while still preserving positivity of the overall surface tension (heterogeneous diffuse interface). Numerical simulations of two-dimensional Ginzburg-Landau phase field equations implementing such mechanism and including hydrodynamic motion, show evidence of dynamically arrested domain coarsening. Under specific conditions on the functional form of the surface tension density, dynamical arrest can be interpreted in terms of the collective dynamics of metastable, non-linear excitations of the density field, named compactons, as they are localized to finite-size regions of configuration space and strictly zero elsewhere. Aside from compactons, the heterogeneous diffuse interface scenario appears to provide a robust mechanism for the interpretation of many aspects of soft-glassy behaviour in binary fluid mixtures.

Shear banding from lattice kinetic models with competing interactions

We present numerical simulations based on a Boltzmann kinetic model with competing interactions, aimed at characterizing the rheological properties of soft-glassy materials. The lattice kinetic model is shown to reproduce typical signatures of driven soft-glassy flows in confined geometries, such as Herschel-Bulkley rheology, shear banding and hysteresis. This lends further credit to the present lattice kinetic model as a valuable tool for the theoretical/computational investigation of the rheology of driven soft-glassy materials under confinement.

Stability and dynamics of droplets on patterned substrates: insights from experiments and lattice Boltzmann simulations

The stability and dynamics of droplets on solid substrates are studied both theoretically and via experiments. Focusing on our recent achievements within the DFG-priority program 1164 (Nano- and Microfluidics), we first consider the case of (large) droplets on the so-called gradient substrates. Here the term gradient refers to both a change of wettability (chemical gradient) or topography (roughness gradient). While the motion of a droplet on a perfectly flat substrate upon the action of a chemical gradient appears to be a natural consequence of the considered situation, we show that the behavior of a droplet on a gradient of topography is less obvious. Nevertheless, if care is taken in the choice of the topographic patterns (in order to reduce hysteresis effects), a motion may be observed. Interestingly, in this case, simple scaling arguments adequately account for the dependence of the droplet velocity on the roughness gradient (Moradi et al 2010 Europhys. Lett. 89 26006). Another issue addressed in this paper is the behavior of droplets on hydrophobic substrates with a periodic arrangement of square shaped pillars. Here, it is possible to propose an analytically solvable model for the case where the droplet size becomes comparable to the roughness scale (Gross et al 2009 Europhys. Lett. 88 26002). Two important predictions of the model are highlighted here. (i) There exists a state with a finite penetration depth, distinct from the full wetting (Wenzel) and suspended (Cassie-Baxter, CB) states. (ii) Upon quasi-static evaporation, a droplet initially on the top of the pillars (CB state) undergoes a transition to this new state with a finite penetration depth but then (upon further evaporation) climbs up the pillars and goes back to the CB state again. These predictions are confirmed via independent numerical simulations. Moreover, we also address the fundamental issue of the internal droplet dynamics and the terminal center of mass velocity on a flat substrate.

Phase-field model of long-time glasslike relaxation in binary fluid mixtures

We present a new phase-field model for binary fluids, exhibiting typical signatures of soft-glassy behavior, such as long-time relaxation, aging, and long-term dynamical arrest. The present model allows the cost of building an interface to vanish locally within the interface, while preserving positivity of the overall surface tension. A crucial consequence of this property, which we prove analytically, is the emergence of free-energy minimizing density configurations, hereafter named "compactons," to denote their property of being localized to a finite-size region of space and strictly zero elsewhere (no tails). Thanks to compactness, any arbitrary superposition of compactons still is a free-energy minimizer, which provides a direct link between the complexity of the free-energy landscape and the morphological complexity of configurational space.

Shear stress in nonideal fluid lattice Boltzmann simulations

The viscous stress in nonideal fluid lattice Boltzmann methods is investigated theoretically and by simulations. Three representative liquid-gas models are compared in a steady-state situation, where an analytical expression for the viscous stress is available. It is shown that, in the presence of nonideal fluid interactions or strong body forces, the accuracy of the computed viscous stress depends strongly on the underlying model.

Thermal fluctuations in the lattice Boltzmann method for nonideal fluids

We introduce thermal fluctuations in the lattice Boltzmann method for nonideal fluids. A fluctuation-dissipation theorem is derived within the Langevin framework and applied to a specific lattice Boltzmann model that approximates the linearized fluctuating Navier-Stokes equations for fluids based on square-gradient free-energy functionals. The obtained thermal noise is shown to ensure equilibration of all degrees of freedom in a simulation to high accuracy. Furthermore, we demonstrate that satisfactory results for most practical applications of fluctuating hydrodynamics can already be achieved using thermal noise derived in the long-wavelength limit.

Graphics processing unit implementation of lattice Boltzmann models for flowing soft systems

A graphic processing unit (GPU) implementation of the multicomponent lattice Boltzmann equation with multirange interactions for soft-glassy materials ["glassy" lattice Boltzmann (LB)] is presented. Performance measurements for flows under shear indicate a GPU/CPU speed up in excess of 10 for 1024(2) grids. Such significant speed up permits to carry out multimillion time-steps simulations of 1024(2) grids within tens of hours of GPU time, thereby considerably expanding the scope of the glassy LB toward the investigation of long-time relaxation properties of soft-flowing glassy materials.