Lattice Boltzmann model for simulating flows with multiple phases and components

Multicomponent interparticle-potential lattice Boltzmann model for fluids with large viscosity ratios

This work focuses on an improved multicomponent interparticle-potential lattice Boltzmann model. The model results in viscosity-independent equilibrium densities and is capable of simulating kinematic viscosity ratios greater than 1000. External forces are incorporated into the discrete Boltzmann equation, rather than through an equilibrium velocity shift as in the original Shan and Chen (hereafter, SC) model. The model also requires the derivation of a momentum conserving effective velocity, which is substituted into the equilibrium distribution function and applies to both the single- and multiple-relaxation-time formulations. Additionally, higher-order isotropy is used in the calculation of the fluid-fluid interaction forces to reduce the magnitude of spurious currents (i.e., numerical errors) in the vicinity of interfaces. First, we compare the model to the SC model for static bubble simulations. We demonstrate that the model results in viscosity-independent equilibrium bubble densities for a wide range of kinematic viscosities, which is not the case for the SC model. Furthermore, we show that the model is capable of simulating stable bubbles for kinematic viscosity ratios greater than 1000 (when higher-order isotropy is used), whereas the SC model is known to be limited to kinematic viscosity ratios on the order of 10. Next we verify the model for surface tension via Laplace's law and show that the model results in the same surface tension values for a range of kinematic viscosities and kinematic viscosity ratios of 10, 100, and 1000. The model is also verified for layered cocurrent flow though parallel plates. We show that the simulated velocity profiles preserve continuity at the interface for kinematic viscosity ratios ranging from 0.001 to 1000 and that the model accurately predicts nonwetting and wetting phase relative permeability for kinematic viscosity ratios of 0.01 to 100.

Forcing scheme in pseudopotential lattice Boltzmann model for multiphase flows

The pseudopotential lattice Boltzmann (LB) model is a widely used multiphase model in the LB community. In this model, an interaction force, which is usually implemented via a forcing scheme, is employed to mimic the molecular interactions that cause phase segregation. The forcing scheme is therefore expected to play an important role in the pseudoepotential LB model. In this paper, we aim to address some key issues about forcing schemes in the pseudopotential LB model. First, theoretical and numerical analyses will be made for Shan-Chen's forcing scheme [Shan and Chen, Phys. Rev. E 47, 1815 (1993)] and the exact-difference-method forcing scheme [Kupershtokh et al., Comput. Math. Appl. 58, 965 (2009)]. The nature of these two schemes and their recovered macroscopic equations will be shown. Second, through a theoretical analysis, we will reveal the physics behind the phenomenon that different forcing schemes exhibit different performances in the pseudopotential LB model. Moreover, based on the analysis, we will present an improved forcing scheme and numerically demonstrate that the improved scheme can be treated as an alternative approach to achieving thermodynamic consistency in the pseudopotential LB model.

Effect of the forcing term in the multiple-relaxation-time lattice Boltzmann equation on the shear stress or the strain rate tensor

In this work, the effect of the forcing term (or external force) in the multiple-relaxation-time lattice Boltzmann equation (MRTLBE) on the shear stress or the strain rate tensor is studied theoretically and numerically. Through a Chapman-Enskog analysis and numerical simulations, we show that the shear stress (or the strain rate tensor) derived from the MRTLBE is second-order accurate in space. We then examine the influence of the forcing term on the shear stress or the strain rate tensor, and demonstrate that the forcing term effect must be included when the shear stress or the strain rate tensor is computed with the nonequilibrium part of the distribution function.

High-order thermal lattice Boltzmann models derived by means of Gauss quadrature in the spherical coordinate system

We use the spherical coordinate system in the momentum space and an appropriate discretization procedure to derive a hierarchy of lattice Boltzmann (LB) models with variable temperature. The separation of the integrals in the momentum space into angular and radial parts allows us to compute the moments of the equilibrium distribution function by means of Gauss-Legendre and Gauss-Laguerre quadratures, as well as to find the elements of the discrete momentum set for each LB model in the hierarchy. The capability of the high-order models in this hierarchy to capture specific effects in microfluidics is investigated through a computer simulation of Couette flow by using the Shakhov collision term to get the right value of the Prandtl number.

Three-dimensional lattice Boltzmann model for immiscible two-phase flow simulations

We present an improved three-dimensional 19-velocity lattice Boltzmann model for immisicible binary fluids with variable viscosity and density ratios. This model uses a perturbation step to generate the interfacial tension and a recoloring step to promote phase segregation and maintain surfaces. A generalized perturbation operator is derived using the concept of a continuum surface force together with the constraints of mass and momentum conservation. A theoretical expression for the interfacial tension is determined directly without any additional analysis and assumptions. The recoloring algorithm proposed by Latva-Kokko and Rothman is applied for phase segregation, which minimizes the spurious velocities and removes lattice pinning. This model is first validated against the Laplace law for a stationary bubble. It is found that the interfacial tension is predicted well for density ratios up to 1000. The model is then used to simulate droplet deformation and breakup in simple shear flow. We compute droplet deformation at small capillary numbers in the Stokes regime and find excellent agreement with the theoretical Taylor relation for the segregation parameter β=0.7. In the limit of creeping flow, droplet breakup occurs at a critical capillary number 0.35<Ca(c)<0.4 for the viscosity ratio of unity, consistent with previous numerical simulations and experiments. Droplet breakup can also be promoted by increasing the Reynolds number. Finally, we numerically investigate a single bubble rising under buoyancy force in viscous fluids for a wide range of Eötvös and Morton numbers. Numerical results are compared with theoretical predictions and experimental results, and satisfactory agreement is shown.

Convection in multiphase fluid flows using lattice Boltzmann methods

We present high-resolution numerical simulations of convection in multiphase flows (boiling) using a novel algorithm based on a lattice Boltzmann method. We first study the thermodynamical and kinematic properties of the algorithm. Then, we perform a series of 3D numerical simulations changing the mean properties in the phase diagram and compare convection with and without phase coexistence at Rayleigh number Ra∼10(7). We show that in the presence of nucleating bubbles non-Oberbeck-Boussinesq effects develop, the mean temperature profile becomes asymmetric, and heat-transfer and heat-transfer fluctuations are enhanced, at all Ra studied. We also show that small-scale properties of velocity and temperature fields are strongly affected by the presence of the buoyant bubble leading to high non-gaussian profiles in the bulk.

Towards a mesoscopic model of water-like fluids with hydrodynamic interactions

We present a mesoscopic lattice model for non-ideal fluid flows with directional interactions, mimicking the effects of hydrogen bonds in water. The model supports a rich and complex structural dynamics of the orientational order parameter, and exhibits the formation of disordered domains whose size and shape depend on the relative strength of directional order and thermal diffusivity. By letting the directional forces carry an inverse density dependence, the model is able to display a correlation between ordered domains and low density regions, reflecting the idea of water as a denser liquid in the disordered state than in the ordered one.

Additional interfacial force in lattice Boltzmann models for incompressible multiphase flows

The existing lattice Boltzmann models for incompressible multiphase flows are mostly constructed with two distribution functions: one is the order parameter distribution function, which is used to track the interface between different phases, and the other is the pressure distribution function for solving the velocity field. In this paper, it is shown that in these models the recovered momentum equation is inconsistent with the target one: an additional force is included in the recovered momentum equation. The additional force has the following features. First, it is proportional to the macroscopic velocity. Second, it is zero in every single-phase region but is nonzero in the interface. Therefore it can be interpreted as an interfacial force. To investigate the effects of the additional interfacial force, numerical simulations are carried out for the problem of Rayleigh-Taylor instability, droplet splashing on a thin liquid film, and the evolution of a falling droplet under gravity. Numerical results demonstrate that, with the increase of the velocity or the Reynolds number, the additional interfacial force will gradually have an important influence on the interface and affect the numerical accuracy.

Forcing term in single-phase and Shan-Chen-type multiphase lattice Boltzmann models

Numerous schemes have been proposed to incorporate a bulk forcing term into the lattice Boltzmann equation. In this paper we present a simple and straightforward comparative analysis of five popular schemes [Shan and Chen, Phys. Rev. E 47, 1815 (1993); Phys Rev Lett. 81, 1618 (1998); He et al., Phys. Rev. E 57, R13 (1998); Guo et al., Phys. Rev. E 65, 046308 (2002); Kupershtokh et al., Comput. Math. Appl. 58, 965 (2009)] in which their differences and similarities are identified. From the analysis we classify the schemes into two groups; the behaviors of the schemes in each group are proven to be identical up to second order. Numerical test simulating the two-dimensional unsteady Taylor-Green vortex flow problem demonstrate that all five schemes are of comparable accuracy for single-phase flow. However, for two-phase flow the situation is different, which is demonstrated by incorporating these schemes into different Shan-Chen-type multiphase models. The forcing scheme in the original Shan-Chen (SC) multiphase model turns out to be inaccurate in terms of the resulting surface tension for different density ratios and relaxation times. In the numerical tests, a typical equation of state and interparticle interactions including next-nearest neighbors were incorporated into the SC model. Our results confirm that the surface-tension values obtained from the original SC lattice Boltzmann method (LBM) simulation depend on the value of the relaxation time τ. For τ<0.7Δt, the surface tension agree well with the analytical solutions. However, when τ>0.7Δt, the surface tension turns out to be systematically larger than the analytical one, exceeding it by more than a factor of 2 for τ=2Δt. In contrast, with the application of the scheme proposed by He et al., the SC LBM produces very accurate surface tensions independent of the value of τ. We also found that the densities of the coexisting liquid and gas can be adjusted to match those at thermodynamic equilibrium if the particle interaction term includes next-nearest-neighbor contributions. The obtained results will be useful for further studies of two-phase flow with high density ratios using the SC LBM approach.

Phase separation in thermal systems: a lattice Boltzmann study and morphological characterization

We investigate thermal and isothermal symmetric liquid-vapor separations via a fast Fourier transform thermal lattice Boltzmann (FFT-TLB) model. Structure factor, domain size, and Minkowski functionals are employed to characterize the density and velocity fields, as well as to understand the configurations and the kinetic processes. Compared with the isothermal phase separation, the freedom in temperature prolongs the spinodal decomposition (SD) stage and induces different rheological and morphological behaviors in the thermal system. After the transient procedure, both the thermal and isothermal separations show power-law scalings in domain growth, while the exponent for thermal system is lower than that for isothermal system. With respect to the density field, the isothermal system presents more likely bicontinuous configurations with narrower interfaces, while the thermal system presents more likely configurations with scattered bubbles. Heat creation, conduction, and lower interfacial stresses are the main reasons for the differences in thermal system. Different from the isothermal case, the release of latent heat causes the changing of local temperature, which results in new local mechanical balance. When the Prandtl number becomes smaller, the system approaches thermodynamical equilibrium much more quickly. The increasing of mean temperature makes the interfacial stress lower in the following way: σ=σ(0)[(T(c)-T)/(T(c)-T(0))](3/2), where T(c) is the critical temperature and σ(0) is the interfacial stress at a reference temperature T(0), which is the main reason for the prolonged SD stage and the lower growth exponent in the thermal case. Besides thermodynamics, we probe how the local viscosities influence the morphology of the phase separating system. We find that, for both the isothermal and thermal cases, the growth exponents and local flow velocities are inversely proportional to the corresponding viscosities. Compared with the isothermal case, the local flow velocity depends not only on viscosity but also on temperature.

Consistent pseudopotential interactions in lattice Boltzmann models

We outline a systematic procedure on how to construct and derive interactions for nonideal lattice fluids. Using a mesoscopic approach based on exact lattice theories, we prove the consistency of the resulting diffuse interface theory with continuum thermodynamics. Translated in the framework of the so called "Shan-Chen" model for nonideal lattice fluids, this paper shows how to adjust the associated pseudopotentials to reproduce a free energy model based on a square gradient theory of equilibrium interfaces.

Lattice Boltzmann method for one-dimensional radiation transfer

The macroscopic conservation equations of radiation energy and radiation momentum are derived on the basis of radiation hydrodynamics. Based on the Chapman-Enskog method, the lattice Boltzmann model for one-dimensional radiative transfer is proposed from the Boltzmann equation. The numerical simulation results agree well with the exact solution and show that the lattice Boltzmann method developed in this paper has good accuracy and stability for solving one-dimensional radiative transfer problems.

Lattice Boltzmann simulations of incompressible liquid-gas systems on partial wetting surfaces

A three-dimensional Lattice Boltzmann two-phase model capable of dealing with large liquid and gas density ratios and with a partial wetting surface is introduced. This is based on a high density ratio model combined with a partial wetting boundary method. The predicted three-dimensional droplets at different partial wetting conditions at equilibrium are in good agreement with analytical solutions. Despite the large density ratio, the spurious velocity around the interface is not substantial, and is rather insensitive to the examined liquid and gas density and viscosity ratios. The influence of the gravitational force on the droplet shape is also examined through the variations of the Bond number, where the droplet shape migrates from spherical to flattened interface in tandem with the increase of the Bond number. The predicted interfaces under constant Bond number are also validated against measurements with good agreements.

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.

Thermodynamic consistency in deriving lattice Boltzmann models for describing segregation in non-ideal mixtures

The thermodynamic consistency of kinetic models for non-ideal mixtures in non-isothermal conditions is investigated. A kinetic model is proposed that is suitable for deriving high-order lattice Boltzmann equations by an appropriate discretization of the velocity space, satisfying the Galilean invariance condition and free of spurious terms in the first moment equations.

A local lattice Boltzmann method for multiple immiscible fluids and dense suspensions of drops

The lattice Boltzmann method (LBM) for computational fluid dynamics benefits from a simple, explicit, completely local computational algorithm making it highly efficient. We extend LBM to recover hydrodynamics of multi-component immiscible fluids, while retaining a completely local, explicit and simple algorithm. Hence, no computationally expensive lattice gradients, interaction potentials or curvatures, that use information from neighbouring lattice sites, need to be calculated, which makes the method highly scalable and suitable for high performance parallel computing. The method is analytical and is shown to recover correct continuum hydrodynamic equations of motion and interfacial boundary conditions. This LBM may be further extended to situations containing a high number (O(100)) of individually immiscible drops. We make comparisons of the emergent non-Newtonian behaviour with a power-law fluid model. We anticipate our method will have a range applications in engineering, industrial and biological sciences.

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.

Microfluidic emulsion separation-simultaneous separation and sensing by multilayer nanofilm structures

Emulsion separation is of high relevance for filtration applications, liquid-liquid-partitioning of biomolecules like proteins and recovery of products from droplet microreactors. Selective interaction of various components of an emulsion with substrates is used to design microfluidic flow chambers for efficient separation of emulsions into their individual components. Our lab-on-a-chip device consists of an emulsion separation cell with an integrated silicon sensor chip, the latter allowing the detection of liquid motion via the field-effect signal. Thus, within our lab-on-a-chip device, emulsions can be separated while the separation process is monitored simultaneously. For emulsion separation a surface energy step gradient, namely a sharp interface between the hydrophobic and hydrophilic parts of the separation chamber, is used. The key component of the lab-on-a-chip system is a multilayer and multifunctional nanofilm structure which not only provides the surface energy step gradient for emulsion separation but also constitutes the functional parts of the field-effect transistors. The proof-of-principle was performed using a model emulsion consisting of immiscible aqueous and organic solvent components. Droplet coalescence was identified as a key aspect influencing the separation process, with quite different effects during separation on open surfaces as compared to slit geometry. For a detailed description of this observation, an analytical model was derived and lattice Boltzmann computer simulations were performed. By use of grazing incidence small angle x-ray scattering (GISAXS) interfacial nanostructures during gold nanoparticle deposition in a flow field were probed to demonstrate the potential of GISAXS for in situ investigations during flow.

Simulations of slip flow on nanobubble-laden surfaces

On microstructured hydrophobic surfaces, geometrical patterns may lead to the appearance of a superhydrophobic state, where gas bubbles at the surface can have a strong impact on the fluid flow along such surfaces. In particular, they can strongly influence a detected slip at the surface. We present two-phase lattice Boltzmann simulations of a flow over structured surfaces with attached gas bubbles and demonstrate how the detected slip depends on the pattern geometry, the bulk pressure, or the shear rate. Since a large slip leads to reduced friction, our results give assistance in the optimization of microchannel flows for large throughput.

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.

Advances in modelling of biomimetic fluid flow at different scales

The biomimetic flow at different scales has been discussed at length. The need of looking into the biological surfaces and morphologies and both geometrical and physical similarities to imitate the technological products and processes has been emphasized. The complex fluid flow and heat transfer problems, the fluid-interface and the physics involved at multiscale and macro-, meso-, micro- and nano-scales have been discussed. The flow and heat transfer simulation is done by various CFD solvers including Navier-Stokes and energy equations, lattice Boltzmann method and molecular dynamics method. Combined continuum-molecular dynamics method is also reviewed.

Simulation of evaporation by an extension of the pseudopotential lattice Boltzmann method: a quantitative analysis

An extension of the pseudopotential lattice Boltzmann method is introduced to simulate heat transfer problems involving phase transition. Using this model, evaporation through a plane interface and two-phase Poiseuille flow were simulated and the macroscopic jump conditions were utilized to evaluate the accuracy of the method. We have found that the simulation results are in very good agreement with the analytical solutions as far as we take into account the extent of the interface during the evaluation. Using the same model heterogeneous boiling was simulated taking into account the geometry of a cavity and the important features of the boiling process could be observed in the simulation results.

From bijels to Pickering emulsions: a lattice Boltzmann study

Particle stabilized emulsions are ubiquitous in the food and cosmetics industry, but our understanding of the influence of microscopic fluid-particle and particle-particle interactions on the macroscopic rheology is still limited. In this paper we present a simulation algorithm based on a multicomponent lattice Boltzmann model to describe the solvents combined with a molecular dynamics solver for the description of the solved particles. It is shown that the model allows a wide variation of fluid properties and arbitrary contact angles on the particle surfaces. We demonstrate its applicability by studying the transition from a "bicontinuous interfacially jammed emulsion gel" (bijel) to a "Pickering emulsion" in dependence on the contact angle, the particle concentration, and the ratio of the solvents.

The dynamics of the impact and coalescence of droplets on a solid surface

A simple experimental setup to study the impact and coalescence of deposited droplets is described. Droplet impact and coalescence have been investigated by high-speed particle image velocimetry. Velocity fields near the liquid-substrate interface have been observed for the impact and coalescence of 2.4 mm diameter droplets of glycerol∕water striking a flat transparent substrate in air. The experimental arrangement images the internal flow in the droplets from below the substrate with a high-speed camera and continuous laser illumination. Experimental results are in the form of digital images that are processed by particle image velocimetry and image processing algorithms to obtain velocity fields, droplet geometries, and contact line positions. Experimental results are compared with numerical simulations by the lattice Boltzmann method.

Lattice Boltzmann method for n-dimensional nonlinear hyperbolic conservation laws with the source term

It is important for nonlinear hyperbolic conservation laws (NHCL) to own a simulation scheme with high order accuracy, simple computation, and non-oscillatory character. In this paper, a unified and novel lattice Boltzmann model is presented for solving n-dimensional NHCL with the source term. By introducing the high order source term of explicit lattice Boltzmann method (LBM) and the optimum dimensionless relaxation time varied with the specific issues, the effects of space and time resolutions on the accuracy and stability of the model are investigated for the different problems in one to three dimensions. Both the theoretical analysis and numerical simulation validate that the results by the proposed LBM have second-order accuracy in both space and time, which agree well with the analytical solutions.

Force imbalance in lattice Boltzmann equation for two-phase flows

The capability of modeling and simulating complex interfacial dynamics of multiphase flows has been recognized as one of the main advantages of the lattice Boltzmann equation (LBE). A basic feature of two-phase LBE models, i.e., the force balance condition at the discrete lattice level of LBE, is investigated in this work. An explicit force-balance formulation is derived for a flat interface by analyzing the two-dimensional nine-velocity (D2Q9) two-phase LBE model without invoking the Chapman-Enskog expansion. The result suggests that generally the balance between the interaction force and the pressure does not hold exactly on the discrete lattice due to numerical errors. It is also shown that such force imbalance can lead to some artificial velocities in the vicinity of phase interface.

Lattice Boltzmann method on quadtree grids

A lattice Boltzmann method on nonuniform quadtree grids is proposed. Our method employs the interpolation-supplemented lattice Boltzmann model. The advantages of the quadtree grid are preserved by using linear interpolation instead of quadratic interpolation to complete the streaming step in the lattice Boltzmann method. The back-and-forth error compensation and correction (BFECC) method is used to improve the accuracy, so that the second-order accuracy of the conventional lattice Boltzmann method is maintained. Several numerical cases, including a BFECC streaming test, lid-driven cavity flow, and flow over an asymmetrically placed cylinder in a channel, are carried out to demonstrate the accuracy and efficiency of our method.

Lattice-Boltzmann method combined with smoothed-profile method for particulate suspensions

We developed a simulation scheme based on the coupling of the lattice-Boltzmann method with the smoothed-profile method (SPM) to predict the dynamic behavior of colloidal dispersions. The SPM provides a coupling scheme between continuum fluid dynamics and rigid-body dynamics through a smoothed profile of the fluid-particle interface. In this approach, the flow is computed on fixed Eulerian grids which are also used for the particles. Owing to the use of the same grids for simulation of fluid flow and particles, this method is highly efficient. Furthermore, an external boundary is used to impose the no-slip boundary condition at the fluid-particle interface. In addition, the operations in the present method are local; it can be easily programmed for parallel machines. The methodology is validated by comparing with previously published data.

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.

Thermodynamically Constrained Averaging Theory Approach for Modeling Flow and Transport Phenomena in Porous Medium Systems: 8. Interface and Common Curve Dynamics

This work is the eighth in a series that develops the fundamental aspects of the thermodynamically constrained averaging theory (TCAT) that allows for a systematic increase in the scale at which multiphase transport phenomena is modeled in porous medium systems. In these systems, the explicit locations of interfaces between phases and common curves, where three or more interfaces meet, are not considered at scales above the microscale. Rather, the densities of these quantities arise as areas per volume or length per volume. Modeling of the dynamics of these measures is an important challenge for robust models of flow and transport phenomena in porous medium systems, as the extent of these regions can have important implications for mass, momentum, and energy transport between and among phases, and formulation of a capillary pressure relation with minimal hysteresis. These densities do not exist at the microscale, where the interfaces and common curves correspond to particular locations. Therefore, it is necessary for a well-developed macroscale theory to provide evolution equations that describe the dynamics of interface and common curve densities. Here we point out the challenges and pitfalls in producing such evolution equations, develop a set of such equations based on averaging theorems, and identify the terms that require particular attention in experimental and computational efforts to parameterize the equations. We use the evolution equations developed to specify a closed two-fluid-phase flow model.

Interpolated lattice Boltzmann boundary conditions for surface reaction kinetics

This paper describes a method for implementing surface reaction kinetics in lattice Boltzmann simulations. The interpolated boundary conditions are capable of simulating surface reactions and dissolution at both stationary and moving solid-fluid and fluid-fluid interfaces. Results obtained with the boundary conditions are compared to analytical solutions for first-order and constant-flux kinetic surface reactions in a one-dimensional half space, as well as to the analytical solution for evaporation from the surface of a cylinder. Excellent agreement between analytical and simulated results is obtained for a wide range of diffusivities, lattice velocities, and surface reaction rates. The boundary model's ability to represent dissolution in binary fluid mixtures is demonstrated by modeling diffusion from a rising bubble and dissolution of a droplet near a flat plate.

Lattice Boltzmann method for Lennard-Jones fluids based on the gradient theory of interfaces

In the present study we propose a lattice Boltzmann equation (LBE) model derived from density gradient expansions of the discrete BBGKY evolution equations. The model is based on the mechanical approach of the gradient theory of interfaces. The basic input is the radial distribution function, which is related exclusively to the molecular interaction potential, rather than semiempirical equations of state used in previous LBE models. This function can be provided from independent molecular simulations or from approximate theories. Evidently the accuracy of the interaction potential, and thus the radial distribution function, reflects on the accuracy of the thermodynamic properties and consistency of the derived LBE model. We have applied the proposed model to obtain equilibrium bulk and interfacial properties of a Lennard-Jones fluid at different temperatures, T, close to critical, T(c). The results demonstrate that the LBE model is in excellent agreement with gradient theory as well as with independent literature results based on different molecular simulation approaches. Hence the proposed LBE model can recover accurately bulk and interfacial thermodynamics for a Lennard Jones fluid at T/T(c)>0.9.

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.