Lattice Boltzmann model for simulating flows with multiple phases and components

Recursive regularization step for high-order lattice Boltzmann methods

A lattice Boltzmann method (LBM) with enhanced stability and accuracy is presented for various Hermite tensor-based lattice structures. The collision operator relies on a regularization step, which is here improved through a recursive computation of nonequilibrium Hermite polynomial coefficients. In addition to the reduced computational cost of this procedure with respect to the standard one, the recursive step allows to considerably enhance the stability and accuracy of the numerical scheme by properly filtering out second- (and higher-) order nonhydrodynamic contributions in under-resolved conditions. This is first shown in the isothermal case where the simulation of the doubly periodic shear layer is performed with a Reynolds number ranging from 10^{4} to 10^{6}, and where a thorough analysis of the case at Re=3×10^{4} is conducted. In the latter, results obtained using both regularization steps are compared against the Bhatnagar-Gross-Krook LBM for standard (D2Q9) and high-order (D2V17 and D2V37) lattice structures, confirming the tremendous increase of stability range of the proposed approach. Further comparisons on thermal and fully compressible flows, using the general extension of this procedure, are then conducted through the numerical simulation of Sod shock tubes with the D2V37 lattice. They confirm the stability increase induced by the recursive approach as compared with the standard one.

Insights into the Impact of Surface Hydrophobicity on Droplet Coalescence and Jumping Dynamics

Droplet coalescence jumping on superhydrophobic surfaces attracts much research attention owing to its capability in enhancing condensation for energy and water applications. In this work, we reveal the impact of the finite surface adhesion to explain velocity discrepancies observed in recent droplet jumping studies, particularly when droplet sizes are a few micrometers (1-10 μm). Surface adhesion, which is usually neglected, can significantly affect both droplet coalescence and departure dynamics. It causes oscillations on velocity and contact area in the droplet coalescence process, as observed numerically and experimentally. Comparing the increasing rate of jumping velocity with contact angle for three different droplet sizes, we show that smaller droplets exhibit higher sensitivity to the change of surface hydrophobicity. We also specify the range of surface superhydrophobicity where the jumping velocity monotonically decreases (θ ≳ 170°), increases (θ ≲ 160°), or changes non-monotonically in transition (160° ≲ θ ≲170°) with droplet size. As a result, there exists a broad jumping velocity range for micrometer-sized droplets on a superhydrophobic surface with a slight contact angle variation. This work offers an extended understanding of the droplet coalescence and jumping dynamics to resolve the discrepancies in recent experimental observations.

Lattice Boltzmann solution of the transient Boltzmann transport equation in radiative and neutron transport

Applications of the transient Boltzmann transport equation (BTE) have undergone much investigation, such as radiative heat transfer and neutron transport. This paper provides a lattice Boltzmann model to efficiently resolve the multidimensional transient BTE. For a higher angular resolution, enough transport directions are considered while the transient BTE in each direction is treated as a conservation law equation and solved independently. Both macroscopic equations recovered from a Chapman-Enskog expansion and simulated results of typical benchmark problems show not only the second-order accuracy but also the flexibility and applicability of the proposed lattice Boltzmann model. This approach may contribute a powerful technique for the parallel simulation of large-scale engineering and some alternative perspectives for solving the nonlinear transport problem further.

Lattice Boltzmann modeling and simulation of liquid jet breakup

A three-dimensional color-fluid lattice Boltzmann model for immiscible two-phase flows is developed in the framework of a three-dimensional 27-velocity (D3Q27) lattice. The collision operator comprises the D3Q27 versions of three suboperators: a multiple-relaxation-time (MRT) collision operator, a generalized Liu-Valocchi-Kang perturbation operator, and a Latva-Kokko-Rothman recoloring operator. A D3Q27 version of an enhanced equilibrium distribution function is also incorporated into this model to improve the Galilean invariance. Three types of numerical tests, namely, a static droplet, an oscillating droplet, and the Rayleigh-Taylor instability, show a good agreement with analytical solutions and numerical simulations. Following these numerical tests, this model is applied to liquid-jet-breakup simulations. The simulation conditions are matched to the conditions of the previous experiments. In this case, numerical stability is maintained throughout the simulation, although the kinematic viscosity for the continuous phase is set as low as 1.8×10^{-4}, in which case the corresponding Reynolds number is 3.4×10^{3}; the developed lattice Boltzmann model based on the D3Q27 lattice enables us to perform the simulation with parameters directly matched to the experiments. The jet's liquid column transitions from an asymmetrical to an axisymmetrical shape, and entrainment occurs from the side of the jet. The measured time history of the jet's leading-edge position shows a good agreement with the experiments. Finally, the reproducibility of the regime map for liquid-liquid systems is assessed. The present lattice Boltzmann simulations well reproduce the characteristics of predicted regimes, including varicose breakup, sinuous breakup, and atomization.

Chemical-potential-based lattice Boltzmann method for nonideal fluids

Chemical potential, as an important thermodynamic quantity, has been popularly used in thermodynamic modeling for complex systems, especially for those involving the phase transitions and chemical reactions. Here we present a chemical-potential-based multiphase lattice Boltzmann model, in which the nonideal force is directly evaluated by a chemical potential. The numerical computation is more efficient than the pressure-tensor-based model [Wen et al. Europhys. Lett. 112, 44002 (2015)10.1209/0295-5075/112/44002] because the calculations of the pressure tensor and its divergence are avoided. We have derived several chemical potentials of the popular equations of state from the free-energy density function. The theoretical analyses and numerical results support that the present model satisfies thermodynamics and Galilean invariance. An effective chemical-potential boundary condition is also implemented to investigate the wettability of a solid surface, and the contact angle can be linearly tuned by the surface chemical potential.

Thermohydrodynamics of an evaporating droplet studied using a multiphase lattice Boltzmann method

In this paper, the thermohydrodynamics of an evaporating droplet is investigated by using a single-component pseudopotential lattice Boltzmann model. The phase change is applied to the model by adding source terms to the thermal lattice Boltzmann equation in such a way that the macroscopic energy equation of multiphase flows is recovered. In order to gain an exhaustive understanding of the complex hydrodynamics during evaporation, a single droplet is selected as a case study. At first, some tests for a stationary (non-)evaporating droplet are carried out to validate the method. Then the model is used to study the thermohydrodynamics of a falling evaporating droplet. The results show that the model is capable of reproducing the flow dynamics and transport phenomena of a stationary evaporating droplet quite well. Of course, a moving droplet evaporates faster than a stationary one due to the convective transport. Our study shows that our single-component model for simulating a moving evaporating droplet is limited to low Reynolds numbers.

Internal dynamics of microgels: A mesoscale hydrodynamic simulation study

We analyze the dynamics of polymers in a microgel system under different swelling conditions. A microgel particle consists of coarse-grained linear polymers which are tetra-functionally crosslinked and undergoes conformational changes in response to the external stimuli. Here, a broad range of microgel sizes, extending from tightly collapsed to strongly swollen particles, is considered. In order to account for hydrodynamic interactions, the microgel is embedded in a multiparticle collision dynamics fluid while hydrophobic attraction is modelled by an attractive Lennard-Jones potential and swelling of ionic microgels is described through the Debye-Hückel potential. The polymer dynamics is analyzed in terms of the monomer mean square displacement and the intermediate scattering function S(q, t). The scattering function decays in a stretched-exponential manner, with a decay rate exhibiting a crossover from a collective diffusive dynamics at low magnitudes of the wavevector q to a hydrodynamic-dominated dynamics at larger q. There is little difference between the intermediate scattering functions of microgels under good solvent conditions and strongly swollen gels, but strongly collapsed gels exhibit a faster decay at short times and hydrodynamic interactions become screened. In addition, we present results for the dynamics of the crosslinks, which exhibit an unexpected, semiflexible polymer-like dynamics.

Modeling mass transfer and reaction of dilute solutes in a ternary phase system by the lattice Boltzmann method

In this work, we propose a general approach for modeling mass transfer and reaction of dilute solute(s) in incompressible three-phase flows by introducing a collision operator in lattice Boltzmann (LB) method. An LB equation was used to simulate the solute dynamics among three different fluids, in which the newly expanded collision operator was used to depict the interface behavior of dilute solute(s). The multiscale analysis showed that the presented model can recover the macroscopic transport equations derived from the Maxwell-Stefan equation for dilute solutes in three-phase systems. Compared with the analytical equation of state of solute and dynamic behavior, these results are proven to constitute a generalized framework to simulate solute distributions in three-phase flows, including compound soluble in one phase, compound adsorbed on single-interface, compound in two phases, and solute soluble in three phases. Moreover, numerical simulations of benchmark cases, such as phase decomposition, multilayered planar interfaces, and liquid lens, were performed to test the stability and efficiency of the model. Finally, the multiphase mass transfer and reaction in Janus droplet transport in a straight microchannel were well reproduced.

Analysis of force treatment in the pseudopotential lattice Boltzmann equation method

In this paper, different force treatments are analyzed in detail for a pseudopotential lattice Boltzmann equation (LBE), and the contribution of third-order error terms to pressure tensor with a force scheme is analyzed by a higher-order Chapman-Enskog expansion technique. From the theoretical analysis, the performance of the original force treatment of Shan-Chen (SC), Ladd, Guo et al., and the exact difference method (EDM) are ɛ_{Ladd}<ɛ_{Guo}<ɛ_{EDM}≤ɛ_{SC} with the relaxation time τ≥1, while ɛ_{Ladd}<ɛ_{Guo}<ɛ_{SC}<ɛ_{EDM} with τ<1; here ɛ is a parameter related to the mechanical stability and the subscripts are the corresponding force scheme. To be consistent with the thermodynamic theory, a force term is introduced to modify the coefficients in the pressure tensor. Some numerical simulations are conducted to show that the predictions of modified force treatment of the pseudopotential LBE are all in good agreement with the analytical solution and other predictions.

Direction Dependence of Adhesion Force for Droplets on Rough Substrates

Determination of solid surface free energy is still an open problem. At present, there are two leading theories on how to determine the adhesion of droplets on rough substrates: one theory stresses that the droplet adhesion force lies in the areas of contact and interaction energy between liquid and solid molecules, whereas the other holds that the length of the edge of drops is essential. In this work, we unify the two theories through lattice Boltzmann simulations and demonstrate that the adhesion force could depend on either the contact area or the contact line, depending on the direction of the adhesion force measured, that is, by vertically separating the two materials or laterally sliding the droplet on the substrate. We reveal that for separating droplets away from rough substrates, the vertical adhesion (pull-off) force depends more significantly on the contact area rather than on the contact line. However, for sliding a droplet on substrates, the lateral adhesion force depends on the contact line while being independent of the contact area.

Generalized three-dimensional lattice Boltzmann color-gradient method for immiscible two-phase pore-scale imbibition and drainage in porous media

This article presents a three-dimensional numerical framework for the simulation of fluid-fluid immiscible compounds in complex geometries, based on the multiple-relaxation-time lattice Boltzmann method to model the fluid dynamics and the color-gradient approach to model multicomponent flow interaction. New lattice weights for the lattices D3Q15, D3Q19, and D3Q27 that improve the Galilean invariance of the color-gradient model as well as for modeling the interfacial tension are derived and provided in the Appendix. The presented method proposes in particular an approach to model the interaction between the fluid compound and the solid, and to maintain a precise contact angle between the two-component interface and the wall. Contrarily to previous approaches proposed in the literature, this method yields accurate solutions even in complex geometries and does not suffer from numerical artifacts like nonphysical mass transfer along the solid wall, which is crucial for modeling imbibition-type problems. The article also proposes an approach to model inflow and outflow boundaries with the color-gradient method by generalizing the regularized boundary conditions. The numerical framework is first validated for three-dimensional (3D) stationary state (Jurin's law) and time-dependent (Washburn's law and capillary waves) problems. Then, the usefulness of the method for practical problems of pore-scale flow imbibition and drainage in porous media is demonstrated. Through the simulation of nonwetting displacement in two-dimensional random porous media networks, we show that the model properly reproduces three main invasion regimes (stable displacement, capillary fingering, and viscous fingering) as well as the saturating zone transition between these regimes. Finally, the ability to simulate immiscible two-component flow imbibition and drainage is validated, with excellent results, by numerical simulations in a Berea sandstone, a frequently used benchmark case used in this field, using a complex geometry that originates from a 3D scan of a porous sandstone. The methods presented in this article were implemented in the open-source PALABOS library, a general C++ matrix-based library well adapted for massive fluid flow parallel computation.

Droplet Dynamics of Newtonian and Inelastic Non-Newtonian Fluids in Confinement

Microfluidic droplet technology has been developing rapidly. However, precise control of dynamical behaviour of droplets remains a major hurdle for new designs. This study is to understand droplet deformation and breakup under simple shear flow in confined environment as typically found in microfluidic applications. In addition to the Newtonian–Newtonian system, we consider also both a Newtonian droplet in a non-Newtonian matrix fluid and a non-Newtonian droplet in a Newtonian matrix. The lattice Boltzmann method is adopted to systematically investigate droplet deformation and breakup under a broad range of capillary numbers, viscosity ratios of the fluids, and confinement ratios considering shear-thinning and shear-thickening fluids. Confinement is found to enhance deformation, and the maximum deformation occurs at the viscosity ratio of unity. The droplet orients more towards the flow direction with increasing viscosity ratio or confinement ratio. In addition, it is noticed that the wall effect becomes more significant for confinement ratios larger than 0.4. Finally, for the whole range of Newtonian carrier fluids tested, the critical capillary number above which droplet breakup occurs is only slightly affected by the confinement ratio for a viscosity ratio of unity. Upon increasing the confinement ratio, the critical capillary number increases for the viscosity ratios less than unity, but decreases for the viscosity ratios more than unity.

Pseudopotential lattice Boltzmann equation method for two-phase flow: A higher-order Chapmann-Enskog expansion

In this paper, a higher order Chapmann-Enskog expansion technique is applied to pseudopotential lattice Boltzmann equation (LBE), and the contribution of third order error terms to pressure tensor is analyzed in detail. To be consistent with the thermodynamic theory, a force term is introduced to modify the coefficients in the pressure tensor. Some numerical simulations are conducted to validate the LBE, and the results show that the predictions of the present LBE agree well with the analytical solution and other predictions.

Lattice Boltzmann Simulation of Droplets Impacting on Superhydrophobic Surfaces with Randomly Distributed Rough Structures

Superhydrophobic surfaces have attracted much attention in environmental control because of their excellent water-repellent properties. A successful design of superhydrophobic surfaces requires a correct understanding of the influences of surface roughness on water-repellent behaviors. Here, a new approach, a mesoscale lattice Boltzmann simulation approach, is proposed and used to model the dynamic behavior of droplets impacting on surfaces with randomly distributed rough microstructures. The fast Fourier transformation method is used to generate non-Gaussian randomly distributed rough surfaces, with the skewness and kurtosis obtained from real surfaces. Then, droplets impacting on the rough surfaces are modeled. It is found that the shape of droplet spreading is obviously affected by the distributions of surface asperity. Decreasing the skewness and keeping the kurtosis around 3 is an effective method to enhance the ability of droplet rebound. The new approach gives more detailed insights into the design of superhydrophobic surfaces.

Pore-scale lattice Boltzmann simulation of flow and mass transfer in bioreactor with an immobilized granule for biohydrogen production

The photo bioreaction combined with flow and mass transfer is simulated with pore-scale lattice Boltzmann (LB) method, which is the scenario of a bioreactor filled with a porous granule immobilized photosynthetic bacteria cells for hydrogen production. The quartet structure generation set (QSGS) is used to generate porous structure of the immobilized granule. The effects of porosity of the immobilized granule on flow and concentration fields as well as the hydrogen production performance are investigated. Higher porosity facilitates the substrate solution smoothly flowing through the porous granule with increasing velocity, and thus results in higher product concentration inside the immobilized granule. Additionally, the substrate consumption efficiency increases, while hydrogen yield slightly decreases with increasing porosity, and they tend to stable for the porosity larger than 0.5. Furthermore, the LB numerical results have a good agreement with the experimental results. It is demonstrated that the pore-scale LB simulation method coupling with QSGS is available to simulate the photo hydrogen production in the bioreactor with porous immobilized granules.

Droplet migration on hydrophobic–hydrophilic hybrid surfaces: a lattice Boltzmann study

Droplet migration on textured hydrophobic–hydrophilic surfaces.

Lattice Boltzmann simulations of droplet formation in confined channels with thermocapillary flows

Based on mesoscale lattice Boltzmann simulations with the "Shan-Chen" model, we explore the influence of thermocapillarity on the breakup properties of fluid threads in a microfluidic T-junction, where a dispersed phase is injected perpendicularly into a main channel containing a continuous phase, and the latter induces periodic breakup of droplets due to the cross-flowing. Temperature effects are investigated by switching on-off both positive-negative temperature gradients along the main channel direction, thus promoting a different thread dynamics with anticipated-delayed breakup. Numerical simulations are performed at changing the flow rates of both the continuous and dispersed phases, as well as the relative importance of viscous forces, surface tension forces, and thermocapillary stresses. The range of parameters is broad enough to characterize the effects of thermocapillarity on different mechanisms of breakup in the confined T-junction, including the so-called "squeezing" and "dripping" regimes, previously identified in the literature. Some simple scaling arguments are proposed to rationalize the observed behavior, and to provide quantitative guidelines on how to predict the droplet size after breakup.

Chimaera simulation of complex states of flowing matter

We discuss a unified mesoscale framework (chimaera) for the simulation of complex states of flowing matter across scales of motion. The chimaera framework can deal with each of the three macro-meso-micro levels through suitable 'mutations' of the basic mesoscale formulation. The idea is illustrated through selected simulations of complex micro- and nanoscale flows.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.

Cascaded lattice Boltzmann method with improved forcing scheme for large-density-ratio multiphase flow at high Reynolds and Weber numbers

A recently developed forcing scheme has allowed the pseudopotential multiphase lattice Boltzmann method to correctly reproduce coexistence curves, while expanding its range to lower surface tensions and arbitrarily high density ratios [Lycett-Brown and Luo, Phys. Rev. E 91, 023305 (2015)PLEEE81539-375510.1103/PhysRevE.91.023305]. Here, a third-order Chapman-Enskog analysis is used to extend this result from the single-relaxation-time collision operator, to a multiple-relaxation-time cascaded collision operator, whose additional relaxation rates allow a significant increase in stability. Numerical results confirm that the proposed scheme enables almost independent control of density ratio, surface tension, interface width, viscosity, and the additional relaxation rates of the cascaded collision operator. This allows simulation of large density ratio flows at simultaneously high Reynolds and Weber numbers, which is demonstrated through binary collisions of water droplets in air (with density ratio up to 1000, Reynolds number 6200 and Weber number 440). This model represents a significant improvement in multiphase flow simulation by the pseudopotential lattice Boltzmann method in which real-world parameters are finally achievable.

Gray free-energy multiphase lattice Boltzmann model with effective transport and wetting properties

The paper shows that it is possible to combine the free-energy lattice Boltzmann approach to multiphase modeling of fluids involving both liquid and vapor with the partial bounce back lattice Boltzmann approach to modeling effective media. Effective media models are designed to mimic the properties of porous materials with porosity much finer than the scale of the simulation lattice. In the partial bounce-back approach, an effective media parameter or bounce-back fraction controls fluid transport. In the combined model, a wetting potential is additionally introduced that controls the wetting properties of the fluid with respect to interfaces between free space (white nodes), effective media (gray nodes), and solids (black nodes). The use of the wetting potential combined with the bounce-back parameter gives the model the ability to simulate transport and sorption of a wide range of fluid in material systems. Results for phase separation, permeability, contact angle, and wicking in gray media are shown. Sorption is explored in small sections of model multiscale porous systems to demonstrate two-step desorption, sorption hysteresis, and the ink-bottle effect.

Revised Chapman-Enskog analysis for a class of forcing schemes in the lattice Boltzmann method

In the lattice Boltzmann (LB) method, the forcing scheme, which is used to incorporate an external or internal force into the LB equation, plays an important role. It determines whether the force of the system is correctly implemented in an LB model and affects the numerical accuracy. In this paper we aim to clarify a critical issue about the Chapman-Enskog analysis for a class of forcing schemes in the LB method in which the velocity in the equilibrium density distribution function is given by u=∑_{α}e_{α}f_{α}/ρ, while the actual fluid velocity is defined as u[over ̂]=u+δ_{t}F/(2ρ). It is shown that the usual Chapman-Enskog analysis for this class of forcing schemes should be revised so as to derive the actual macroscopic equations recovered from these forcing schemes. Three forcing schemes belonging to the above class are analyzed, among which Wagner's forcing scheme [A. J. Wagner, Phys. Rev. E 74, 056703 (2006)10.1103/PhysRevE.74.056703] is shown to be capable of reproducing the correct macroscopic equations. The theoretical analyses are examined and demonstrated with two numerical tests, including the simulation of Womersley flow and the modeling of flat and circular interfaces by the pseudopotential multiphase LB model.

Droplet dynamics in rotating flows

This paper deals with investigations of droplet dynamics in rotating flows. In many previous studies droplet dynamics was analyzed in simple unidirectional flows. To fill this gap, the focus of this study is an overview on investigations of droplet dynamics in a complex rotating flow. A Lattice Boltzmann Method with high potential in simulation of two-phase unsteady flows is applied to simulate the physics of the problem in a lid-driven cavity. In spite of its simple geometry, there is a complex rotating flow field containing different vortices and shear regions. The Reynolds number based on the cavity length scale and the upper wall velocity, ReL, is considered to be 1000. We discuss here effects of different parameters such as: density ratios (1, 5, 10, 100, and 1000), droplet sizes (D/L=0.097, 0.114, 0.131 and 0.2), and droplet initial positions (1/8, 2/8, and 3/8 of the cavity length, L, out of center). The results are discussed in terms of global flow physics and its interaction with the droplet, drop deformation during its motion along with the main flow, and droplet trajectories. It is shown that there are strong interactions between the droplet and the main carrying flow. During motion, the droplets pass through different flow regions containing acceleration/deceleration zones. Consequently, the droplets experience different shear forces resulting in stretching, shrinking, rotating and dilatation which all contribute to the dynamics of the droplet.

Pinning-Depinning Mechanism of the Contact Line during Evaporation on Chemically Patterned Surfaces: A Lattice Boltzmann Study

In this paper, the pinning and depinning mechanism of the contact line during droplet evaporation on chemically stripe-patterned surfaces is numerically investigated using a thermal multiphase lattice Boltzmann (LB) model with liquid-vapor phase change. A local force balance in the context of diffuse interfaces is introduced to explain the equilibrium states of droplets on chemically patterned surfaces. It is shown that when the contact line is pinned on a hydrophobic-hydrophilic boundary, different contact angles can be interpreted as the variation of the length of the contact line occupied by each component. The stick-slip-jump behavior of evaporating droplets on chemically patterned surfaces is well captured by the LB simulations. Particularly, a slow movement of the contact line is clearly observed during the stick (pinning) mode, which shows that the pinning of the contact line during droplet evaporation on chemically stripe-patterned surfaces is actually a dynamic pinning process and the dynamic equilibrium is achieved by the self-adjustment of the contact lines occupied by each component. Moreover, it is shown that when the surface tension varies with the temperature, the Marangoni effect has an important influence on the depinning of the contact line, which occurs when the horizontal component (toward the center of the droplet) of the force caused by the Marangoni stress overcomes the unbalanced Young's force toward the outside.

Unidirectional Fast Growth and Forced Jumping of Stretched Droplets on Nanostructured Microporous Surfaces

Superhydrophobic nanostructured surfaces have demonstrated outstanding capability in energy and water applications by promoting dropwise condensation, where fast droplet growth and efficient condensate removal are two key parameters. However, these parameters remain contradictory. Although efficient droplet removal is easily obtained through coalescence jumping on uniform superhydrophobic surfaces, simultaneously achieving fast droplet growth is still challenging. Also, on such surfaces droplets can grow to larger sizes without restriction if there is no coalescence. In this work, we show that superhydrophobic nanostructured microporous surfaces can manipulate the droplet growth and jumping. Microporous surface morphology effectively enhances the growth of droplets in pores owing to large solid-liquid contact area. At low supersaturations, the upward growth rate (1-1.5 μm/s) of these droplets in pores is observed to be around 15-25 times that of the droplets outside the pores. Meanwhile, their top curvature radius increases relatively slowly (∼0.25 μm/s) due to pore confinement, which results in a highly stretched droplet surface. We also observed forced jumping of stretched droplets in pores either through coalescence with spherical droplets outside pores or through self-pulling without coalescence. Both experimental observation and theoretical modeling reveal that excess surface free energy stored in the stretched droplet surface and micropore confinement are responsible for this pore-scale-forced jumping. These findings reveal the insightful physics of stretched droplet dynamics and offer guidelines for the design and fabrication of novel super-repellent surfaces with microporous morphology.

Controlled capillary assembly of magnetic Janus particles at fluid-fluid interfaces

Capillary interactions can be used to direct assembly of particles adsorbed at fluid-fluid interfaces. Precisely controlling the magnitude and direction of capillary interactions to assemble particles into favoured structures for materials science purposes is desirable but challenging. In this paper, we investigate capillary interactions between magnetic Janus particles adsorbed at fluid-fluid interfaces. We develop a pair-interaction model that predicts that these particles should arrange into a side-side configuration, and carry out simulations that confirm the predictions of our model. Finally, we investigate the monolayer structures that form when many magnetic Janus particles adsorb at the interface. We find that the particles arrange into long, straight chains exhibiting little curvature, in contrast with capillary interactions between ellipsoidal particles. We further find a regime in which highly ordered, lattice-like monolayer structures form, which can be tuned dynamically using an external magnetic field.

Rebound suppression of a droplet impacting on an oscillating horizontal surface

The behavior of a droplet impinging onto a solid substrate can be influenced significantly by the horizontal motion of the substrate. The coupled interactions between the moving wall and the impacting droplet may result in various outcomes, which may be different from the usual normal droplet impact on a stationary wall. In this paper, we present a method to suppress drop rebound on hydrophobic surfaces via transverse wall oscillations, normal to the impact direction. The numerical investigation shows that the suppression of droplet rebound has a direct relationship with the oscillation phase, amplitude, and frequency. For a particular range of oscillation frequencies and amplitudes, a lateral shifting of the droplet position is observed along the oscillating direction. While large oscillation amplitude favors the process of droplet deposition, a high frequency promotes droplet rebound from the oscillating wall. A linear trend in the transition region between deposition and rebound is observed from a scaled phase diagram of the oscillation amplitude versus frequency. We provide a systematic investigation of drop deposition and elucidate the mechanism of rebound suppression through the temporal evolution of the nonaxial kinetic energy and the velocity flow field.

Multiple-relaxation-time color-gradient lattice Boltzmann model for simulating two-phase flows with high density ratio

In this paper we propose a color-gradient lattice Boltzmann (LB) model for simulating two-phase flows with high density ratio and high Reynolds number. The model applies a multirelaxation-time (MRT) collision operator to enhance the stability of the simulation. A source term, which is derived by the Chapman-Enskog analysis, is added into the MRT LB equation so that the Navier-Stokes equations can be exactly recovered. Also, a form of the equilibrium density distribution function is used to simplify the source term. To validate the proposed model, steady flows of a static droplet and the layered channel flow are first simulated with density ratios up to 1000. Small values of spurious velocities and interfacial tension errors are found in the static droplet test, and improved profiles of velocity are obtained by the present model in simulating channel flows. Then, two cases of unsteady flows, Rayleigh-Taylor instability and droplet splashing on a thin film, are simulated. In the former case, the density ratio of 3 and Reynolds numbers of 256 and 2048 are considered. The interface shapes and spike and bubble positions are in good agreement with the results of previous studies. In the latter case, the droplet spreading radius is found to obey the power law proposed in previous studies for the density ratio of 100 and Reynolds number up to 500.

Contact angle adjustment in equation-of-state-based pseudopotential model

The single component pseudopotential lattice Boltzmann model has been widely applied in multiphase simulation due to its simplicity and stability. In many studies, it has been claimed that this model can be stable for density ratios larger than 1000. However, the application of the model is still limited to small density ratios when the contact angle is considered. The reason is that the original contact angle adjustment method influences the stability of the model. Moreover, simulation results in the present work show that, by applying the original contact angle adjustment method, the density distribution near the wall is artificially changed, and the contact angle is dependent on the surface tension. Hence, it is very inconvenient to apply this method with a fixed contact angle, and the accuracy of the model cannot be guaranteed. To solve these problems, a contact angle adjustment method based on the geometry analysis is proposed and numerically compared with the original method. Simulation results show that, with our contact angle adjustment method, the stability of the model is highly improved when the density ratio is relatively large, and it is independent of the surface tension.

Grayscale lattice Boltzmann model for multiphase heterogeneous flow through porous media

The grayscale lattice Boltzmann (LB) model has been recently developed to model single-phase fluid flow through heterogeneous porous media. Flow is allowed in each voxel but the degree of flow depends on that voxel's resistivity to fluid motion. Here we extend the grayscale LB model to multiphase, immiscible flow. The new model is outlined and then applied to a number of test cases, which show good agreement with theory. This method is subsequently used to model the important case where each voxel may have a different resistance to each particular fluid that is passing through it. Finally, the method is applied to model fluid flow through real porous media to demonstrate its capability. Both the capillary and viscous flow regimes are recovered in these simulations.

Modeling adsorption with lattice Boltzmann equation

The research of adsorption theory has recently gained renewed attention due to its critical relevance to a number of trending industrial applications, hydrogen storage and shale gas exploration for instance. The existing theoretical foundation, laid mostly in the early twentieth century, was largely based on simple heuristic molecular interaction models and static interaction potential which, although being insightful in illuminating the fundamental mechanisms, are insufficient for computations with realistic adsorbent structure and adsorbate hydrodynamics, both critical for real-life applications. Here we present and validate a novel lattice Boltzmann model incorporating both adsorbate-adsorbate and adsorbate-adsorbent interactions with hydrodynamics which, for the first time, allows adsorption to be computed with real-life details. Connection with the classic Ono-Kondo lattice theory is established and various adsorption isotherms, both within and beyond the IUPAC classification are observed as a pseudo-potential is varied. This new approach not only enables an important physical to be simulated for real-life applications, but also provides an enabling theoretical framework within which the fundamentals of adsorption can be studied.

Lattice Boltzmann method for convection-diffusion equations with general interfacial conditions

In this work, we propose an interfacial scheme accompanying the lattice Boltzmann method for convection-diffusion equations with general interfacial conditions, including conjugate conditions with or without jumps in heat and mass transfer, continuity of macroscopic variables and normal fluxes in ion diffusion in porous media with different porosity, and the Kapitza resistance in heat transfer. The construction of this scheme is based on our boundary schemes [Huang and Yong, J. Comput. Phys. 300, 70 (2015)JCTPAH0021-999110.1016/j.jcp.2015.07.045] for Robin boundary conditions on straight or curved boundaries. It gives second-order accuracy for straight interfaces and first-order accuracy for curved ones. In addition, the new scheme inherits the advantage of the boundary schemes in which only the current lattice nodes are involved. Such an interfacial scheme is highly desirable for problems with complex geometries or in porous media. The interfacial scheme is numerically validated with several examples. The results show the utility of the constructed scheme and very well support our theoretical predications.

Lattice Boltzmann Model of 3D Multiphase Flow in Artery Bifurcation Aneurysm Problem

This paper simulates and predicts the laminar flow inside the 3D aneurysm geometry, since the hemodynamic situation in the blood vessels is difficult to determine and visualize using standard imaging techniques, for example, magnetic resonance imaging (MRI). Three different types of Lattice Boltzmann (LB) models are computed, namely, single relaxation time (SRT), multiple relaxation time (MRT), and regularized BGK models. The results obtained using these different versions of the LB-based code will then be validated with ANSYS FLUENT, a commercially available finite volume- (FV-) based CFD solver. The simulated flow profiles that include velocity, pressure, and wall shear stress (WSS) are then compared between the two solvers. The predicted outcomes show that all the LB models are comparable and in good agreement with the FVM solver for complex blood flow simulation. The findings also show minor differences in their WSS profiles. The performance of the parallel implementation for each solver is also included and discussed in this paper. In terms of parallelization, it was shown that LBM-based code performed better in terms of the computation time required.

Bubble accumulation and its role in the evolution of magma reservoirs in the upper crust

Volcanic eruptions transfer huge amounts of gas to the atmosphere. In particular, the sulfur released during large silicic explosive eruptions can induce global cooling. A fundamental goal in volcanology, therefore, is to assess the potential for eruption of the large volumes of crystal-poor, silicic magma that are stored at shallow depths in the crust, and to obtain theoretical bounds for the amount of volatiles that can be released during these eruptions. It is puzzling that highly evolved, crystal-poor silicic magmas are more likely to generate volcanic rocks than plutonic rocks. This observation suggests that such magmas are more prone to erupting than are their crystal-rich counterparts. Moreover, well studied examples of largely crystal-poor eruptions (for example, Katmai, Taupo and Minoan) often exhibit a release of sulfur that is 10 to 20 times higher than the amount of sulfur estimated to be stored in the melt. Here we argue that these two observations rest on how the magmatic volatile phase (MVP) behaves as it rises buoyantly in zoned magma reservoirs. By investigating the fluid dynamics that controls the transport of the MVP in crystal-rich and crystal-poor magmas, we show how the interplay between capillary stresses and the viscosity contrast between the MVP and the host melt results in a counterintuitive dynamics, whereby the MVP tends to migrate efficiently in crystal-rich parts of a magma reservoir and accumulate in crystal-poor regions. The accumulation of low-density bubbles of MVP in crystal-poor magmas has implications for the eruptive potential of such magmas, and is the likely source of the excess sulfur released during explosive eruptions.