Lattice-Boltzmann model based on field mediators for immiscible fluids

Pore-scale characteristics of multiphase flow in heterogeneous porous media using the lattice Boltzmann method

This study provides a pore-scale investigation of two-phase flow dynamics during primary drainage in a realistic heterogeneous rock sample. Using the lattice Boltzmann (LB) method, a series of three-dimensional (3D) immiscible displacement simulations are conducted and three typical flow patterns are identified and mapped on the capillary number (Ca)-viscosity ratio(M) phase diagram. We then investigate the effect of the viscosity ratio and capillary number on fluid saturation patterns and displacement stability in Tuscaloosa sandstone, which is taken from the Cranfield site. The dependence of the evolution of saturation, location of the displacement front, 3D displacement patterns and length of the center of mass of the invading fluid on the viscosity ratio and capillary number have been delineated. To gain a quantitative insight into the characteristics of the invasion morphology in 3D porous media, the fractal dimension D_f of the non-wetting phase displacement patterns during drainage has been computed for various viscosity ratios and capillary numbers. The logarithmic dependence of D_f on invading phase saturation appears to be the same for various capillary numbers and viscosity ratios and follows a universal relation.

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.

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.