Lattice Boltzmann model for simulating flows with multiple phases and components

Sliding drops across alternating hydrophobic and hydrophilic stripes

We perform a joint numerical and experimental study to systematically characterize the motion of 30 μl drops of pure water and of ethanol in water solutions, sliding over a periodic array of alternating hydrophobic and hydrophilic stripes with a large wettability contrast and a typical width of hundreds of microns. The fraction of the hydrophobic areas has been varied from about 20% to 80%. The effects of the heterogeneous patterning can be described by a renormalized value of the critical Bond number, i.e., the critical dimensionless force needed to depin the drop before it starts to move. Close to the critical Bond number we observe a jerky motion characterized by an evident stick-slip dynamics. As a result, dissipation is strongly localized in time, and the mean velocity of the drops can easily decrease by an order of magnitude compared to the sliding on the homogeneous surface. Lattice Boltzmann numerical simulations are crucial for disclosing to what extent the sliding dynamics can be deduced from the computed balance of capillary, viscous, and body forces by varying the Bond number, the surface composition, and the liquid viscosity. Beyond the critical Bond number, we characterize both experimentally and numerically the dissipation inside the droplet by studying the relation between the average velocity and the applied volume forces.

Numerical simulations of two-fluid turbulent mixing at large density ratios and applications to the Rayleigh-Taylor instability

A tentative review is presented of various approaches for numerical simulations of two-fluid gaseous mixtures at high density ratios, as they have been applied to the Rayleigh-Taylor instability (RTI). Systems exhibiting such RTI behaviour extend from atomistic sizes to scales where the continuum approximation becomes valid. Each level of description can fit into a hierarchy of theoretical models and the governing equations appropriate for each model, with their assumptions, are presented. In particular, because the compressible to incompressible limit of the Navier-Stokes equations is not unique and understanding compressibility effects in the RTI critically depends on having the appropriate basis for comparison, two relevant incompressible limits are presented. One of these limits has not been considered before. Recent results from RTI simulations, spanning the levels of description presented, are reviewed in connection to the material mixing problem. Owing to the computational limitations, most in-depth RTI results have been obtained for the incompressible case. Two such results, concerning the asymmetry of the mixing and small-scale anisotropy anomaly, as well as the possibility of a mixing transition in the RTI, are surveyed. New lines for further investigation are suggested and it is hoped that bringing together such diverse levels of description may provide new ideas and increased motivation for studying such flows.

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

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

Ultrafine particles dispersion modeling in a street canyon: development and evaluation of a composite lattice Boltzmann model

Recently, a branch of particulate matter research concerns on ultrafine particles found in the urban environment, which originate, to a significant extent, from traffic sources. In urban street canyons, dispersion of ultrafine particles affects pedestrian's short term exposure and resident's long term exposure as well. The aim of the present work is the development and the evaluation of a composite lattice Boltzmann model to study the dispersion of ultrafine particles, in urban street canyon microenvironment. The proposed model has the potential to penetrate into the physics of this complex system. In order to evaluate the model performance against suitable experimental data, ultrafine particles levels have been monitored on an hourly basis for a period of 35 days, in a street canyon, in Athens area. The results of the comparative analysis are quite satisfactory. Furthermore, our modeled results are in a good agreement with the results of other computational and experimental studies. This work is a first attempt to study the dispersion of an air pollutant by application of the lattice Boltzmann method.

High-resolution MRI velocimetry compared with numerical simulations

Alterations of the blood flow are associated with various cardiovascular diseases. Precise knowledge of the velocity distribution is therefore important for understanding these diseases and predicting the effect of different medical intervention schemes. The goal of this work is to estimate the precision with which the velocity field can be measured and predicted by studying two simple model geometries with NMR micro imaging and computational fluid dynamics. For these initial experiments, we use water as an ideal test medium. The phantoms consist of tubes simulating a straight blood vessel and a step between two tubes of different diameters, which can be seen as a minimal model of the situation behind a stenosis. For both models, we compare the experimental data with the numerical prediction, using the experimental boundary conditions. For the simpler model, we also compare the data to the analytical solution. As an additional validation, we determine the divergence of the velocity field and verify that it vanishes within the experimental uncertainties. We discuss the resulting precision of the simulation and the outlook for extending this approach to the analysis of specific cases of arteriovascular problems.

Color-gradient lattice Boltzmann model for simulating droplet motion with contact-angle hysteresis

Lattice Boltzmann method (LBM) is an effective tool for simulating the contact-line motion due to the nature of its microscopic dynamics. In contact-line motion, contact-angle hysteresis is an inherent phenomenon, but it is neglected in most existing color-gradient based LBMs. In this paper, a color-gradient based multiphase LBM is developed to simulate the contact-line motion, particularly with the hysteresis of contact angle involved. In this model, the perturbation operator based on the continuum surface force concept is introduced to model the interfacial tension, and the recoloring operator proposed by Latva-Kokko and Rothman is used to produce phase segregation and resolve the lattice pinning problem. At the solid surface, the color-conserving wetting boundary condition [Hollis et al., IMA J. Appl. Math. 76, 726 (2011)] is applied to improve the accuracy of simulations and suppress spurious currents at the contact line. In particular, we present a numerical algorithm to allow for the effect of the contact-angle hysteresis, in which an iterative procedure is used to determine the dynamic contact angle. Numerical simulations are conducted to verify the developed model, including the droplet partial wetting process and droplet dynamical behavior in a simple shear flow. The obtained results are compared with theoretical solutions and experimental data, indicating that the model is able to predict the equilibrium droplet shape as well as the dynamic process of partial wetting and thus permits accurate prediction of contact-line motion with the consideration of contact-angle hysteresis.

Lattice-Boltzmann-based two-phase thermal model for simulating phase change

A lattice Boltzmann (LB) method is presented for solving the energy conservation equation in two phases when the phase change effects are included in the model. This approach employs multiple distribution functions, one for a pseudotemperature scalar variable and the rest for the various species. A nonideal equation of state (EOS) is introduced by using a pseudopotential LB model. The evolution equation for the pseudotemperature variable is constructed in such a manner that in the continuum limit one recovers the well known macroscopic energy conservation equation for the mixtures. Heats of reaction, the enthalpy change associated with the phase change, and the diffusive transport of enthalpy are all taken into account; but the dependence of enthalpy on pressure, which is usually a small effect in most nonisothermal flows encountered in chemical reaction systems, is ignored. The energy equation is coupled to the LB equations for species transport and pseudopotential interaction forces through the EOS by using the filtered local pseudotemperature field. The proposed scheme is validated against simple test problems for which analytical solutions can readily be obtained.

Stick-slip sliding of water drops on chemically heterogeneous surfaces

We present a comprehensive study of water drops sliding down chemically heterogeneous surfaces formed by a periodic pattern of alternating hydrophobic and hydrophilic stripes. Drops are found to undergo a stick-slip motion whose average speed is an order of magnitude smaller than that measured on a homogeneous surface having the same static contact angle. This motion is the result of the periodic deformations of the drop interface when crossing the stripes. Numerical simulations confirm this view and are used to elucidate the principles underlying the experimental observations.

Mixing and internal dynamics of droplets impacting and coalescing on a solid surface

The coalescence and mixing of a sessile and an impacting liquid droplet on a solid surface are studied experimentally and numerically in terms of lateral separation and droplet speed. Two droplet generators are used to produce differently colored droplets. Two high-speed imaging systems are used to investigate the impact and coalescence of the droplets in color from a side view with a simultaneous gray-scale view from below. Millimeter-sized droplets were used with dynamical conditions, based on the Reynolds and Weber numbers, relevant to microfluidics and commercial inkjet printing. Experimental measurements of advancing and receding static contact angles are used to calibrate a contact angle hysteresis model within a lattice Boltzmann framework, which is shown to capture the observed dynamics qualitatively and the final droplet configuration quantitatively. Our results show that no detectable mixing occurs during impact and coalescence of similar-sized droplets, but when the sessile droplet is sufficiently larger than the impacting droplet vortex ring generation can be observed. Finally we show how a gradient of wettability on the substrate can potentially enhance mixing.

Data structure and movement for lattice-based simulations

We show that for the lattice Boltzmann model, the existing paradigm in computer science for the choice of the data structure is suboptimal. In this paper we use the requirements of physical symmetry necessary for recovering hydrodynamics in the lattice Boltzmann description to propose a hybrid data layout for the method. This hybrid data structure, which we call a structure of an array of structures, is shown to be optimal for the lattice Boltzmann model. Finally, the possible advantages of establishing a connection between group theoretic symmetry requirements and the construction of the data structure is discussed in the broader context of grid-based methods.

Interaction pressure tensor for a class of multicomponent lattice Boltzmann models

We present a theory to obtain the pressure tensor for a class of nonideal multicomponent lattice Boltzmann models, thus extending the theory presented by X. Shan [Phys. Rev. E 77, 066702 (2008)] for single-component fluids. We obtain the correct form of the pressure tensor directly on the lattice and the resulting equilibrium properties are shown to agree very well with those measured from numerical simulations. Results are compared with those of alternative theories.

Effect of hydrodynamic and fluid-solid interaction forces on the shape and stability of a droplet sedimenting on a horizontal wall

In this paper, we study the sedimentation of droplet onto horizontal surfaces using a Shan and Chen multicomponent lattice Boltzmann model [Shan and Chen, Phys. Rev. E 47, 1815 (1993)]. Numerical simulations are performed for static nonwetting droplets and are compared with a theoretical model. Further simulations are performed using wetting droplets and changing the wetting characteristics of the surface. Novel study of droplet sedimentation with different wetting characteristics of wall surface shows three different regimes of droplet shapes. The process of different regimes can be controlled by changing the wetting characteristics of surface. It is observed that the droplet attains nonequilibrium region as the attraction force of surrounding droplet with the wall increases. This can lead to eventual breakup of the droplet.

Lattice Boltzmann modeling of directional wetting: comparing simulations to experiments

Lattice Boltzmann Modeling (LBM) simulations were performed on the dynamic behavior of liquid droplets on chemically striped patterned surfaces, ultimately with the aim to develop a predictive tool enabling reliable design of future experiments. The simulations accurately mimic experimental results, which have shown that water droplets on such surfaces adopt an elongated shape due to anisotropic preferential spreading. Details of the contact line motion such as advancing of the contact line in the direction perpendicular to the stripes exhibit pronounced similarities in experiments and simulations. The opposite of spreading, i.e., evaporation of water droplets, leads to a characteristic receding motion first in the direction parallel to the stripes, while the contact line remains pinned perpendicular to the stripes. Only when the aspect ratio is close to unity, the contact line also starts to recede in the perpendicular direction. Very similar behavior was observed in the LBM simulations. Finally, droplet movement can be induced by a gradient in surface wettability. LBM simulations show good semiquantitative agreement with experimental results of decanol droplets on a well-defined striped gradient, which move from high- to low-contact angle surfaces. Similarities and differences for all systems are described and discussed in terms of the predictive capabilities of LBM simulations to model direction wetting.

Axisymmetric multiphase lattice Boltzmann method

A lattice Boltzmann method for axisymmetric multiphase flows is presented and validated. The method is capable of accurately modeling flows with variable density. We develop the classic Shan-Chen multiphase model [Phys. Rev. E 47, 1815 (1993)] for axisymmetric flows. The model can be used to efficiently simulate single and multiphase flows. The convergence to the axisymmetric Navier-Stokes equations is demonstrated analytically by means of a Chapmann-Enskog expansion and numerically through several test cases. In particular, the model is benchmarked for its accuracy in reproducing the dynamics of the oscillations of an axially symmetric droplet and on the capillary breakup of a viscous liquid thread. Very good quantitative agreement between the numerical solutions and the analytical results is observed.

Effects of the hierarchical structure of rough solid surfaces on the wetting of microdroplets

We used the lattice Boltzmann method to investigate how the hierarchical structure of a rough solid surface, which in this work is modeled as the microstructure (micropillars) covered with nanostructures (nanopillars), affects the contact angle of microdroplets atop of the solid surface and the wetting transition between the Wenzel and Cassie states. Our simulation results show that the Wenzel-to-Cassie state transition can be achieved by decreasing the fluid-solid attraction, increasing the micropillar spacing, or coating the microstructures with nanostructures. For the effect of the hierarchical structure on the contact angle, we find that the micropillars show a negligible effect on the contact angle, but they may affect the sliding angle. In contrast, it is the nanostructure that determines the contact angle. The contact angle increases with the nanopillar length until reaching a maximal value, but its dependence on the nanopillar spacing becomes more complicated. The contact angle may first increase with the nanopillar spacing and then decreases, or decreases monotonously, depending on whether the liquid enters the nanostructure or not. In this work, we also demonstrate in the presence of contact line pinning, that the pinning effect affects the apparent contact angle.

Gravity-driven thin liquid films over topographical substrates

We investigate the time-dependent evolution of thin liquid films over inclined substrates using a multi-component lattice Boltzmann algorithm. Substrates with and without grooves are considered and the effects of the inclination angle on the dynamics and the coating of the substrates are studied. Our results indicate that the dynamics is enhanced and the ridge height and its displacement are increased by increasing the inclination angle. However, by increasing the inclination angle the maximum depth that can be successfully coated is reduced. Also, although for any given groove depth the width should be larger than a critical value for successful coating, the critical width decreases for smaller inclination angles. For different inclination angles we derive and report the critical sizes of the grooves for successful coating of the substrates.

Evaluation of outflow boundary conditions for two-phase lattice Boltzmann equation

Outflow boundary condition (OBC) is a critical issue in computational fluid dynamics. As a type of numerical method for fluid flows, the lattice Boltzmann equation (LBE) method has gained much success in a variety of complex flows, and certain OBCs have been suggested for the LBE in simulating simple single-phase flows. However, very few discussions on the OBCs have been made for the two-phase LBE method. In this work, three types of OBCs that are widely used in the LBE for single-phase flows, i.e., the Neumann boundary condition, the convective boundary condition, and the extrapolation boundary condition, are extended to a two-phase LBE method and their performances are investigated. The comprehensive results of several two-phase flows show that these boundary conditions behave quite differently in the simulations of two-phase flows. Specifically, it is found that the Neumann boundary condition and the extrapolation boundary condition give rather poor predictions, while the type of convective boundary conditions work well, although the choice of the convection velocity has some slight influences on the results. We also apply these OBC schemes to some other two-phase models, and similar observations are found.

Simulations of droplet coalescence in simple shear flow

Simulating droplet coalescence is challenging because small-scale (tens of nanometers) phenomena determine the behavior of much larger (micrometer- to millimeter-scale) droplets. In general, liquid droplets colliding in a liquid medium coalesce when the capillary number is less than a critical value. We present simulations of droplet collisions and coalescence in simple shear flow using the free-energy binary-liquid lattice Boltzmann method. In previous simulations of low-speed collisions, droplets coalesced at unrealistically high capillary numbers. Simulations of noncoalescing droplets have not been reported, and therefore, the critical capillary number for simulated collisions was unknown. By simulating droplets with radii up to 100 lattice nodes, we determine the critical capillary number for coalescence and quantify the effects of several numerical and geometric parameters. The simulations were performed with a well-resolved interface, a Reynolds number of one, and capillary numbers from 0.01 to 0.2. The ratio of the droplet radius and interface thickness has the greatest effect on the critical capillary number. As in experiments, the critical capillary number decreases with increasing droplet size. A second numerical parameter, the interface diffusivity (Péclet number) also influences the conditions for coalescence: coalescence occurs at higher capillary numbers with lower Péclet numbers (higher diffusivity). The effects of the vertical offset between the droplets and the confinement of the droplets were also studied. Physically reasonable results were obtained and provide insight into the conditions for coalescence. Simulations that match the conditions of experiments reported in the literature remain computationally impractical. However, the scale of the simulations is now sufficiently large that a comparison with experiments involving smaller droplets (≈10 μm) and lower viscosities (≈10(-6) m(2)/s, the viscosity of water) may be possible. Experiments at these conditions are therefore needed to determine the interface thickness and Péclet number that should be used for predictive simulations of coalescence phenomena.

Simulation of intra-aneurysmal blood flow by different numerical methods

The occlusional performance of sole endoluminal stenting of intracranial aneurysms is controversially discussed in the literature. Simulation of blood flow has been studied to shed light on possible causal attributions. The outcome, however, largely depends on the numerical method and various free parameters. The present study is therefore conducted to find ways to define parameters and efficiently explore the huge parameter space with finite element methods (FEMs) and lattice Boltzmann methods (LBMs). The goal is to identify both the impact of different parameters on the results of computational fluid dynamics (CFD) and their advantages and disadvantages. CFD is applied to assess flow and aneurysmal vorticity in 2D and 3D models. To assess and compare initial simulation results, simplified 2D and 3D models based on key features of real geometries and medical expert knowledge were used. A result obtained from this analysis indicates that a combined use of the different numerical methods, LBM for fast exploration and FEM for a more in-depth look, may result in a better understanding of blood flow and may also lead to more accurate information about factors that influence conditions for stenting of intracranial aneurysms.

Hydrodynamically driven colloidal assembly in dip coating

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

Lattice Boltzmann model for thermal binary-mixture gas flows

A lattice Boltzmann model for thermal gas mixtures is derived. The kinetic model is designed in a way that combines properties of two previous literature models, namely, (a) a single-component thermal model and (b) a multicomponent isothermal model. A comprehensive platform for the study of various practical systems involving multicomponent mixture flows with large temperature differences is constructed. The governing thermohydrodynamic equations include the mass, momentum, energy conservation equations, and the multicomponent diffusion equation. The present model is able to simulate mixtures with adjustable Prandtl and Schmidt numbers. Validation in several flow configurations with temperature and species concentration ratios up to nine is presented.

Lattice Boltzmann modeling of multiphase flows at large density ratio with an improved pseudopotential model

Owing to its conceptual simplicity and computational efficiency, the pseudopotential multiphase lattice Boltzmann (LB) model has attracted significant attention since its emergence. In this work, we aim to extend the pseudopotential LB model to simulate multiphase flows at large density ratio and relatively high Reynolds number. First, based on our recent work [Q. Li, K. H. Luo, and X. J. Li, Phys. Rev. E 86, 016709 (2012)], an improved forcing scheme is proposed for the multiple-relaxation-time pseudopotential LB model in order to achieve thermodynamic consistency and large density ratio in the model. Next, through investigating the effects of the parameter a in the Carnahan-Starling equation of state, we find that the interface thickness is approximately proportional to 1/√a. Using a smaller a will lead to a wider interface thickness, which can reduce the spurious currents and enhance the numerical stability of the pseudopotential model at large density ratio. Furthermore, it is found that a lower liquid viscosity can be gained in the pseudopotential model by increasing the kinematic viscosity ratio between the vapor and liquid phases. The improved pseudopotential LB model is numerically validated via the simulations of stationary droplet and droplet oscillation. Using the improved model as well as the above treatments, numerical simulations of droplet splashing on a thin liquid film are conducted at a density ratio in excess of 500 with Reynolds numbers ranging from 40 to 1000. The dynamics of droplet splashing is correctly reproduced and the predicted spread radius is found to obey the power law reported in the literature.

Pore-scale modeling of multiphase reactive transport with phase transitions and dissolution-precipitation processes in closed systems

A pore-scale model based on the lattice Boltzmann (LB) method is developed for multiphase reactive transport with phase transitions and dissolution-precipitation processes. The model combines the single-component multiphase Shan-Chen LB model [X. Shan and H. Chen, Phys. Rev. E 47, 1815 (1993)], the mass transport LB model [S. P. Sullivan et al., Chem. Eng. Sci. 60, 3405 (2005)], and the dissolution-precipitation model [Q. Kang et al., J. Geophys. Res. 111, B05203 (2006)]. Care is taken to handle information on computational nodes undergoing solid-liquid or liquid-vapor phase changes to guarantee mass and momentum conservation. A general LB concentration boundary condition is proposed that can handle various concentration boundaries including reactive and moving boundaries with complex geometries. The pore-scale model can capture coupled nonlinear multiple physicochemical processes including multiphase flow with phase separations, mass transport, chemical reactions, dissolution-precipitation processes, and dynamic evolution of the pore geometries. The model is validated using several multiphase flow and reactive transport problems and then used to study the thermal migration of a brine inclusion in a salt crystal. Multiphase reactive transport phenomena with phase transitions between liquid-vapor phases and dissolution-precipitation processes of the salt in the closed inclusion are simulated and the effects of the initial inclusion size and temperature gradient on the thermal migration are investigated.

Hydrodynamics of Leidenfrost droplets in one-component fluids

Using the dynamic van der Waals theory [Phys. Rev. E 75, 036304 (2007)], we numerically investigate the hydrodynamics of Leidenfrost droplets under gravity in two dimensions. Some recent theoretical predictions and experimental observations are confirmed in our simulations. A Leidenfrost droplet larger than a critical size is shown to be unstable and break up into smaller droplets due to the Rayleigh-Taylor instability of the bottom surface of the droplet. Our simulations demonstrate that an evaporating Leidenfrost droplet changes continuously from a puddle to a circular droplet, with the droplet shape controlled by its size in comparison with a few characteristic length scales. The geometry of the vapor layer under the droplet is found to mainly depend on the droplet size and is nearly independent of the substrate temperature, as reported in a recent experimental study [Phys. Rev. Lett. 109, 074301 (2012)]. Finally, our simulations demonstrate that a Leidenfrost droplet smaller than a characteristic size takes off from the hot substrate because the levitating force due to evaporation can no longer be balanced by the weight of the droplet, as observed in a recent experimental study [Phys. Rev. Lett. 109, 034501 (2012)].

Phase-field-based lattice Boltzmann model for incompressible binary fluid systems with density and viscosity contrasts

A lattice Boltzmann model (LBM) is proposed based on the phase-field theory to simulate incompressible binary fluids with density and viscosity contrasts. Unlike many existing diffuse interface models which are limited to density matched binary fluids, the proposed model is capable of dealing with binary fluids with moderate density ratios. A new strategy for projecting the phase field to the viscosity field is proposed on the basis of the continuity of viscosity flux. The new LBM utilizes two lattice Boltzmann equations (LBEs): one for the interface tracking and the other for solving the hydrodynamic properties. The LBE for interface tracking can recover the Chan-Hilliard equation without any additional terms; while the LBE for hydrodynamic properties can recover the exact form of the divergence-free incompressible Navier-Stokes equations avoiding spurious interfacial forces. A series of 2D and 3D benchmark tests have been conducted for validation, which include a rigid-body rotation, stationary and moving droplets, a spinodal decomposition, a buoyancy-driven bubbly flow, a layered Poiseuille flow, and the Rayleigh-Taylor instability. It is shown that the proposed method can track the interface with high accuracy and stability and can significantly and systematically reduce the parasitic current across the interface. Comparisons with momentum-based models indicate that the newly proposed velocity-based model can better satisfy the incompressible condition in the flow fields, and eliminate or reduce the velocity fluctuations in the higher-pressure-gradient region and, therefore, achieve a better numerical stability. In addition, the test of a layered Poiseuille flow demonstrates that the proposed scheme for mixture viscosity performs significantly better than the traditional mixture viscosity methods.

Thin liquid film flow over substrates with two topographical features

A multicomponent lattice Boltzmann scheme is used to investigate the surface coating of substrates with two topographical features by a gravity-driven thin liquid film. The considered topographies are U- and V-shaped grooves and mounds. For the case of substrates with two grooves, our results indicate that for each of the grooves there is a critical width such that if the groove width is larger than the critical width, the groove can be coated successfully. The critical width of each groove depends on the capillary number, the contact angle, the geometry, and the depth of that groove. The second groove critical width depends on, in addition, the geometry and the depth of the first groove; for two grooves with the same geometries and depths, it is at least equal to that of the first groove. If the second groove width lies between the critical widths, the second groove still can be coated successfully on the condition that the distance between the grooves is considered larger than a critical distance. For considered contact angles and capillary numbers our results indicate that the critical distance is a convex function of the capillary number and the contact angle. Our study also reveals similar results for the case of substrates with a mound and a groove.

Fracture-scale model of immiscible fluid flow

We present a method for modeling immiscible fluid flow in narrow fractures. The method combines a parallel-plate flow model with a recoloration approach adapted from multiple-component lattice-Boltzmann methods. The resulting numerical method is straightforward to implement and accurately reproduces the relevant fluid behavior. To demonstrate the method, single-droplet simulations are compared to analytical solutions that isolate the contributions from the in-plane and normal curvature. Simulations reproducing capillary forces inside a widening aperture are also presented. Excellent agreement is found in all cases considered. Finally, the method's ability to model fracture flow is demonstrated by deriving relative permeability curves for flow through a heterogeneous aperture created between two fractal surfaces.

Scheme for contact angle and its hysteresis in a multiphase lattice Boltzmann method

In this paper, a scheme for specifying contact angle and its hysteresis is incorporated into a multiphase lattice Boltzmann method. The scheme is validated through investigations of the dynamic behaviors of a droplet sliding along two kinds of walls: a smooth (ideal) wall and a rough or chemically inhomogeneous (nonideal) wall. For an ideal wall, the wettability of solid substrates is able to be prescribed. For a nonideal wall, arbitrary contact angle hysteresis can be obtained through adjusting advancing and receding angles. Significantly different phenomena can be recovered for the two kinds of walls. For instance, a droplet on an inclined ideal wall under gravity is impossible to stay stationary. However, the droplet on a nonideal wall may be pinned due to contact angle hysteresis. The steady interface shapes of the droplet on an inclined nonideal wall under gravity or in a shear flow quantitatively agree well with the previous numerical studies. Besides, the complex motion of a droplet creeping like an inchworm could be simulated. The scheme is found suitable for the study of contact line problems with and without contact angle hysteresis.

Phase-field-based lattice Boltzmann finite-difference model for simulating thermocapillary flows

A phase-field-based hybrid model that combines the lattice Boltzmann method with the finite difference method is proposed for simulating immiscible thermocapillary flows with variable fluid-property ratios. Using a phase field methodology, an interfacial force formula is analytically derived to model the interfacial tension force and the Marangoni stress. We present an improved lattice Boltzmann equation (LBE) method to capture the interface between different phases and solve the pressure and velocity fields, which can recover the correct Cahn-Hilliard equation (CHE) and Navier-Stokes equations. The LBE method allows not only use of variable mobility in the CHE, but also simulation of multiphase flows with high density ratio because a stable discretization scheme is used for calculating the derivative terms in forcing terms. An additional convection-diffusion equation is solved by the finite difference method for spatial discretization and the Runge-Kutta method for time marching to obtain the temperature field, which is coupled to the interfacial tension through an equation of state. The model is first validated against analytical solutions for the thermocapillary driven convection in two superimposed fluids at negligibly small Reynolds and Marangoni numbers. It is then used to simulate thermocapillary migration of a three-dimensional deformable droplet and bubble at various Marangoni numbers and density ratios, and satisfactory agreement is obtained between numerical results and theoretical predictions.