Thermal fluctuations in the lattice Boltzmann method for nonideal fluids

Metastable and unstable hydrodynamics in multiphase lattice Boltzmann

Metastability in liquids is at the foundation of complex phase transformation dynamics such as nucleation and cavitation. Intermolecular interaction details, beyond the equation of state, and thermal hydrodynamic fluctuations play a crucial role. However, most numerical approaches suffer from a slow time and space convergence, thus hindering the convergence to the hydrodynamic limit. This work shows that the Shan-Chen lattice Boltzmann model has the unique capability of simulating the hydrodynamics of the metastable state. The structure factor of density fluctuations is theoretically obtained and numerically verified to a high precision, for all simulated wave vectors, reduced temperatures, and pressures, deep into the metastable region. Such remarkable agreement between the theory and simulations leverages the exact implementation at the lattice level of the mechanical equilibrium condition. The static structure factor is found to consistently diverge as the temperature approaches the critical point or the density approaches the spinodal line at a subcritical temperature. Theoretically predicted critical exponents are observed in both cases. Finally, the phase separation in the unstable branch follows the same pattern, i.e., the generation of interfaces with different topology, as observed in molecular dynamics simulations.

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 lattice Boltzmann study of particle settling in a fluctuating multicomponent fluid under confinement

We present mesoscale numerical simulations based on the coupling of the fluctuating lattice Boltzmann method for multicomponent systems with a wetted finite-size particle model. This newly coupled methodologies are used to study the motion of a spherical particle driven by a constant body force in a confined channel with a fixed square cross section. The channel is filled with a mixture of two liquids under the effect of thermal fluctuations. After some validations steps in the absence of fluctuations, we study the fluctuations in the particle's velocity at changing thermal energy, applied force, particle size, and particle wettability. The importance of fluctuations with respect to the mean settling velocity is quantitatively assessed, especially in comparison with unconfined situations. Results show that the expected effects of confinement are very well captured by the numerical simulations, wherein the confinement strongly enhances the importance of velocity fluctuations, which can be one order of magnitude larger than what expected in unconfined domains. The observed findings underscore the versatility of the proposed methodology in highlighting the effects of confinement on the motion of particles in the presence of thermal fluctuations.

Impact of the prequench state of binary fluid mixtures on surface-directed spinodal decomposition

Using lattice Boltzmann simulations we investigate the impact of the amplitude of concentration fluctuations in binary fluid mixtures prior to demixing when in contact with a surface that is preferentially wet by one of the components. We find a bicontinuous structure near the surface for an initial, prequench state of the mixture close to the critical point where the amplitude of concentration fluctuations is large. In contrast, if the initial state of the mixture is not near the critical point and concentration fluctuations are relatively weak, then the morphology is not bicontinuous but remains layered until the very late stages of coarsening. In both cases, it is the morphology of a depletion layer rich in the nonpreferred component that dictates the growth exponent of the thickness of the fluid layer that is in direct contact with the substrate. In the early stages of demixing, we find a growth exponent consistent with a value of 1/4 for a prequench state away from the critical point, which is different from the usual diffusive scaling exponent of 1/3 that we recover for a prequench state close to the critical point. We attribute this to the structure of a depletion layer that is penetrated by tubes of the preferred fluid, connecting the wetting layer to the bulk fluid even in the early stages if the initial state is characterized by concentration fluctuations that are large in amplitude. Furthermore, we find that in the late stages of demixing the flow through these tubes results in significant in-plane concentration variations near the substrate, leading to dropletlike structures with a concentration lower than the average concentration in the wetting layer. This causes a deceleration in the growth of the wetting layer in the very late stages of the demixing. Irrespective of the prequench state of the mixture, the late stages of the demixing process produce the same growth law for the layer thickness, with a scaling exponent of unity usually associated with the impact of hydrodynamic flow fields.

Large Fluctuations in Nonideal Coarse-Grained Systems

Using the recently introduced molecular dynamics lattice gas approach, we test fluctuations of coarse-grained quantities. We show that as soon as the system can no longer be considered an ideal gas fluctuations fail to diminish upon coarse graining as is usually expected. These results suggest that current approaches to simulating fluctuating hydrodynamics may have to be augmented to achieve quantitative results for systems with a nonideal equation of state. The molecular dynamics lattice gas method gives a guidance to the exact nature of the fluctuation in such systems.

Lattice Boltzmann simulations of nonequilibrium fluctuations in a nonideal binary mixture

In recent years the lattice Boltzmann (LB) methodology has been fruitfully extended to include the effects of thermal fluctuations. So far, all studied cases pertain to equilibrium fluctuations, i.e., fluctuations with respect to an equilibrium background state. In this paper we take a step further and present results of fluctuating LB simulations of a binary mixture confined between two parallel walls in the presence of a constant concentration gradient in the wall-to-wall direction. This is a paradigmatic setup for the study of nonequilibrium (NE) fluctuations, i.e., fluctuations with respect to a nonequilibrium state. We analyze the dependence of the structure factors for the hydrodynamical fields on the wave vector q in both the directions parallel and perpendicular to the walls, highlighting the long-range (∼|q|^{-4}) nature of correlations in the NE framework. Results at the small scales (high wave numbers) quantitatively agree with the predictions of fluctuating hydrodynamics without fitting parameters. At larger scales (low wave numbers), however, results show finite-size effects induced by confinement and call for further studies aimed at controlling boundary conditions in the fluctuating LB framework as well as compressibility effects. Moreover, in the presence of a nonideal equation of state of the mixture, we also observe that the (spatially homogeneous) average pressure changes, due to a genuinely new contribution triggered by NE fluctuations. These NE pressure effects are studied at changing the system size and the concentration gradient. Taken all together, we argue that the results of this article are useful and instrumental to boost the applicability of the fluctuating LB methodology in the framework of NE fluctuations, possibly in conjunction with experiments.

Effects of thermal fluctuations in the fragmentation of a nanoligament

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

Integer lattice gas with Monte Carlo collision operator recovers the lattice Boltzmann method with Poisson-distributed fluctuations

We examine a new kind of lattice gas that closely resembles modern lattice Boltzmann methods. This new kind of lattice gas, which we call a Monte Carlo lattice gas, has interesting properties that shed light on the origin of the multirelaxation time collision operator, and it derives the equilibrium distribution for an entropic lattice Boltzmann. Furthermore these lattice gas methods have Galilean invariant fluctuations given by a Poisson statistics, giving further insight into the properties that we should expect for fluctuating lattice Boltzmann methods.

Mesoscopic modelling and simulation of soft matter

The deformability of soft condensed matter often requires modelling of hydrodynamical aspects to gain quantitative understanding. This, however, requires specialised methods that can resolve the multiscale nature of soft matter systems. We review a number of the most popular simulation methods that have emerged, such as Langevin dynamics, dissipative particle dynamics, multi-particle collision dynamics, sometimes also referred to as stochastic rotation dynamics, and the lattice-Boltzmann method. We conclude this review with a short glance at current compute architectures for high-performance computing and community codes for soft matter simulation.

Fluctuating lattice Boltzmann method for the diffusion equation

We derive a fluctuating lattice Boltzmann method for the diffusion equation. The derivation removes several shortcomings of previous derivations for fluctuating lattice Boltzmann methods for hydrodynamic systems. The comparative simplicity of this diffusive system highlights the basic features of this first exact derivation of a fluctuating lattice Boltzmann method.

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.

Fluctuating multicomponent lattice Boltzmann model

Current implementations of fluctuating lattice Boltzmann equations (FLBEs) describe single component fluids. In this paper, a model based on the continuum kinetic Boltzmann equation for describing multicomponent fluids is extended to incorporate the effects of thermal fluctuations. The thus obtained fluctuating Boltzmann equation is first linearized to apply the theory of linear fluctuations, and expressions for the noise covariances are determined by invoking the fluctuation-dissipation theorem directly at the kinetic level. Crucial for our analysis is the projection of the Boltzmann equation onto the orthonormal Hermite basis. By integrating in space and time the fluctuating Boltzmann equation with a discrete number of velocities, the FLBE is obtained for both ideal and nonideal multicomponent fluids. Numerical simulations are specialized to the case where mean-field interactions are introduced on the lattice, indicating a proper thermalization of the system.

Modeling multiphase flow using fluctuating hydrodynamics

Fluctuating hydrodynamics provides a model for fluids at mesoscopic scales where thermal fluctuations can have a significant impact on the behavior of the system. Here we investigate a model for fluctuating hydrodynamics of a single-component, multiphase flow in the neighborhood of the critical point. The system is modeled using a compressible flow formulation with a van der Waals equation of state, incorporating a Korteweg stress term to treat interfacial tension. We present a numerical algorithm for modeling this system based on an extension of algorithms developed for fluctuating hydrodynamics for ideal fluids. The scheme is validated by comparison of measured structure factors and capillary wave spectra with equilibrium theory. We also present several nonequilibrium examples to illustrate the capability of the algorithm to model multiphase fluid phenomena in a neighborhood of the critical point. These examples include a study of the impact of fluctuations on the spinodal decomposition following a rapid quench, as well as the piston effect in a cavity with supercooled walls. The conclusion in both cases is that thermal fluctuations affect the size and growth of the domains in off-critical quenches.

Fluctuating ideal-gas lattice Boltzmann method with fluctuation dissipation theorem for nonvanishing velocities

Current implementations of fluctuating ideal-gas descriptions with the lattice Boltzmann methods are based on a fluctuation dissipation theorem, which, while greatly simplifying the implementation, strictly holds only for zero mean velocity and small fluctuations. We show how to derive the fluctuation dissipation theorem for all k, which was done only for k=0 in previous derivations. The consistent derivation requires, in principle, locally velocity-dependent multirelaxation time transforms. Such an implementation is computationally prohibitively expensive but, with a small computational trick, it is feasible to reproduce the correct FDT without overhead in computation time. It is then shown that the previous standard implementations perform poorly for non vanishing mean velocity as indicated by violations of Galilean invariance of measured structure factors. Results obtained with the method introduced here show a significant reduction of the Galilean invariance violations.

Interfacial roughening in nonideal fluids: dynamic scaling in the weak- and strong-damping regime

Interfacial roughening denotes the nonequilibrium process by which an initially flat interface reaches its equilibrium state, characterized by the presence of thermally excited capillary waves. Roughening of fluid interfaces has been first analyzed by Flekkoy and Rothman [Phys. Rev. Lett. 75, 260 (1995)], where the dynamic scaling exponents in the weakly damped case in two dimensions were found to agree with the Kardar-Parisi-Zhang universality class. We extend this work by taking into account also the strong-damping regime and perform extensive fluctuating hydrodynamics simulations in two dimensions using the Lattice Boltzmann method. We show that the dynamic scaling behavior is different in the weakly and strongly damped case.

Critical dynamics of an isothermal compressible nonideal fluid

A pure fluid at its critical point shows a dramatic slow-down in its dynamics, due to a divergence of the order-parameter susceptibility and the coefficient of heat transport. Under isothermal conditions, however, sound waves provide the only possible relaxation mechanism for order-parameter fluctuations. Here we study the critical dynamics of an isothermal, compressible nonideal fluid via scaling arguments and computer simulations of the corresponding fluctuating hydrodynamics equations. We show that, below a critical dimension of 4, the order-parameter dynamics of an isothermal fluid effectively reduces to "model A," characterized by overdamped sound waves and a divergent bulk viscosity. In contrast, the shear viscosity remains finite above two dimensions. Possible applications of the model are discussed.

Simulation of static critical phenomena in nonideal fluids with the lattice Boltzmann method

A fluctuating nonideal fluid at its critical point is simulated with the lattice Boltzmann method. It is demonstrated that the method, employing a Ginzburg-Landau free energy functional, correctly reproduces the static critical behavior associated with the Ising universality class. A finite-size scaling analysis is applied to determine the critical exponents related to the order parameter, compressibility and specific heat. A particular focus is put on finite-size effects and issues related to the global conservation of the order parameter.

Lattice-Boltzmann-Langevin simulations of binary mixtures

We report a hybrid numerical method for the solution of the Model H fluctuating hydrodynamic equations for binary mixtures. The momentum conservation equations with Landau-Lifshitz stresses are solved using the fluctuating lattice Boltzmann equation while the order parameter conservation equation with Langevin fluxes is solved using stochastic method of lines. Two methods, based on finite difference and finite volume, are proposed for spatial discretization of the order parameter equation. Special care is taken to ensure that the fluctuation-dissipation theorem is maintained at the lattice level in both cases. The methods are benchmarked by comparing static and dynamic correlations and excellent agreement is found between analytical and numerical results. The Galilean invariance of the model is tested and found to be satisfactory. Thermally induced capillary fluctuations of the interface are captured accurately, indicating that the model can be used to study nonlinear fluctuations.

Modelling thermal fluctuations in non-ideal fluids with the lattice Boltzmann method

Recently, we proposed a theoretical framework to include thermal fluctuations into the Lattice Boltzmann (LB) method for non-ideal fluids. Here, we apply a variant thereof to a certain class of force-based non-ideal fluid LB models. We find that ideal-gas-like noise is an exact result of the fluctuation-dissipation theorem in the hydrodynamic regime. It is shown that satisfactory equilibration of the density and fluid momentum can be obtained in a simulation over a wide range of length scales.