Lattice-Boltzmann method for the simulation of multiphase mass transfer and reaction of dilute species

Electrokinetics at liquid-liquid interfaces: Physical models and transport mechanisms

The electrification effects and electrokinetic flow phenomena at immiscible liquid-liquid interfaces have been a subject of scientific inquiry for over a century. Unlike solid-liquid interfaces, liquid-liquid interfaces exhibit not only multiphysical and cross-scale characteristics but also diffuse soft properties, including finite thickness, fluidity, ion adsorbability, and permeability, which introduces diverse interfacial charging mechanisms and conductive dielectric properties, imparting unique characteristics to electrokinetic multiphase flow systems. Electrokinetic multiphase hydrodynamics (EKmHD), grounded in electrochemistry and colloid and interface science, has experienced renewed interest in recent years. This is particularly evident in systems such as the interface between two immiscible electrolyte solutions (ITIES) in electrochemistry, self-propelling droplets in physicochemical hydrodynamics, and digital microfluidics in electromechanics. The multiphase diffuse soft nature of charged liquid-liquid interfaces introduces novel physical scales and theoretical dimensions, positioning EKmHD as a potential foundation for a new interdisciplinary field rather than merely a cross-disciplinary area. This review highlights the need for an integrated research approach that combines interfacial charging mechanisms with electrokinetic flows, alongside a cross-scale modeling framework for interfacial multiphysical transport. It systematically organizes the characteristics of liquid-liquid interfaces from the perspectives of charging mechanisms and electrokinetic behaviors, with particular emphasis on spontaneous partition- and adsorption-induced charging at the interface, and the strong coupling between multiphase diffuse soft interface flow and ion transport. Furthermore, the paper comprehensively summarizes the transport mechanisms of electrokinetic multiphase flows concerning interfacial ion transport and fluid flow, while refining the corresponding dominant dimensionless parameters. Additionally, it systematically consolidates current understanding of typical electrokinetic multiphase flow scenarios, with special focus on potential future research directions. These include the electrokinetic double-sided coupling effects in ITIES systems, solidification and nonlinear effects in droplet/bubble electrophoresis, the validity of the leaky dielectric model, electrokinetic instabilities of jets and ion-selective soft interfaces, and the active and passive control of two-phase electrokinetic wetting dynamics and displacement.

Phase-field-based lattice Boltzmann method for two-phase flows with interfacial mass or heat transfer

In this work, we develop a phase-field-based lattice Boltzmann (LB) method for a two-scalar model of the two-phase flows with interfacial mass or heat transfer. Through the Chapman-Enskog analysis, we show that the present LB method can correctly recover the governing equations for phase field, flow field, and concentration or temperature field. In particular, to derive the two-scalar equations for the mass or heat transfer, we propose a new LB model with an auxiliary source distribution function to describe the extra flux terms, and the discretizations of some derivative terms can be avoided. The accuracy and efficiency of the present LB method are also tested through several benchmark problems, and the influence of mass or heat transfer on the fluid viscosity is further considered by introducing an exponential relation. The numerical results show that the present LB method is suitable for the two-phase flows with interfacial mass or heat transfer.

The Surfactant Role on a Droplet Passing through a Constricted Microchannel in a Pressure-Driven Flow: A Lattice Boltzmann Study

The role of surfactants in the flow of a droplet driven by a pressure gradient through a constricted microchannel is simulated by using our recently developed lattice Boltzmann method. We first study the surfactant role on a droplet flowing through a microchannel with a shrunken square section under different surfactant concentrations and capillary numbers (i.e., imposed pressure gradients). As the surfactant concentration increases, the droplet flow regime first changes from the flow regime I of the droplet getting stuck at the entrance of the constricted channel to the flow regime II of the droplet flowing through the constricted channel with breakup, and then to the flow regime III of the droplet flowing through the constricted channel without breakup. As the capillary number increases, the surfactant role on the number of mother droplets breaking up and the time of mother droplets completely flowing through the constricted section tend to decrease, suggesting that the surfactant effects are gradually weakened. Then, a phase diagram describing how the surfactant concentration and capillary number affect the droplet flow regime is presented. As the surfactant concentration increases, the critical capillary number that distinguishes droplet flow regimes I from II gradually decreases, while the critical capillary number that distinguishes droplet flow regimes II from III first increases and then decreases.

Comparative investigation of a lattice Boltzmann boundary treatment of multiphase mass transport with heterogeneous chemical reactions

Multiphase reactive transport in porous media is an important component of many natural and engineering processes. In the present study, boundary schemes for the continuum species transport-lattice Boltzmann (CST-LB) mass transport model and the multicomponent pseudopotential model are proposed to simulate heterogeneous chemical reactions in a multiphase system. For the CST-LB model, a lattice-interface-tracking scheme for the heterogeneous chemical reaction boundary is provided. Meanwhile, a local-average virtual density boundary scheme for the multicomponent pseudopotential model is formulated based on the work of Li et al. [Li, Yu, and Luo, Phys. Rev. E 100, 053313 (2019)10.1103/PhysRevE.100.053313]. With these boundary treatments, a numerical implementation is put forward that couples the multiphase fluid flow, interfacial species transport, heterogeneous chemical reactions, and porous matrix structural evolution. A series of comparison benchmark cases are investigated to evaluate the numerical performance for different pseudopotential wetting boundary treatments, and an application case of multiphase dissolution in porous media is conducted to validate the present models' ability to solve complex problems. By applying the present LB models with reasonable boundary treatments, multiphase reactive transport in various natural or engineering scenarios can be simulated accurately.

Rayleigh-Plateau Instability of a Particle-Laden Liquid Column: A Lattice Boltzmann Study

Colloidal particles known to be capable of stabilizing fluid-fluid interfaces have been widely applied in emulsion preparation, but their precise role and underlying influencing mechanism remain poorly understood. In this study, a perturbed liquid column with particles evenly distributed on its surface is investigated using a three-dimensional lattice Boltzmann method, which is built upon the color-gradient two-phase flow model but with a new capillary force model and a momentum exchange method for particle dynamics. The developed method is first validated by simulating the wetting behavior of a particle on a fluid interface and the classic Rayleigh-Plateau instability and is then used to explore the effects of particle concentration and contact angle on the capillary instability of the particle-laden liquid column. It is found that increasing the particle concentration can enhance the stability of the liquid column and thus delay the breakup, and the liquid column is most stable under slightly hydrophobic conditions, which corresponds to the lowest initial liquid-gas interfacial free energy. Due to different pressure gradients inside and outside the liquid column and the capillary force being directed away from the neck, hydrophobic particles tend to assemble in a less compact manner near the neck of the deformed liquid column, while hydrophilic particles prefer to gather far away from the neck. For hydrophobic particles, in addition to the influence of the initial liquid-gas interfacial free energy, the self-assembly of particles in a direction opposite to the liquid flow also contributes to opposing the rupture of the liquid column.

Numerical Studies on Cellulose Hydrolysis in Organic-Liquid-Solid Phase Systems with a Liquid Membrane Catalysis Model

The catalytic hydrolysis of cellulose to produce 5-hydroxymethylfurfural (HMF) is a powerful means of biomass resources. The current efficient hydrolysis of cellulose to obtain HMF is dominated by multiphase reaction systems. However, there is still a lack of studies on the synergistic mechanisms and component transport between the various processes of cellulose hydrolysis in a complex multiphase system. In this paper, a liquid membrane catalytic model was developed to simulate the hydrolysis of cellulose and its further reactions, including the adsorption of the liquid membrane on cellulose particles, the consumption of cellulose solid particles, the complex chemical reactions in the liquid membrane, and the transfer of HMF at the phase interface. The simulations show the synergistic effect between cellulose hydrolysis and multiphase mass transfer. We defined an indicator () to characterize the sensitivity of HMF yield to the initial liquid membrane thickness at different reaction stages. decreased gradually when the glucose conversion increased from 0 to 80%, and increased with the thickening of the initial liquid membrane thickness. It was shown that the thickening of the initial liquid membrane thickness promoted the HMF yield under the same glucose conversion. In summary, our results reveal the mechanism of the interaction between multiple physicochemical processes of the cellulose liquid membrane reaction system.

Lattice Boltzmann modeling of interfacial mass transfer in a multiphase system

In the present study, a numerical model based on the lattice Boltzmann method (LBM) is proposed to simulate multiphase mass transfer, referred to as the CST-LB model. This model introduced continuum species transfer (CST) formulation by an additional collision term to model the mass transfer across the multiphase interface. The boundary condition treatment of this model is also discussed. In order to verify the applicability, the CST-LB model is combined with the pseudopotential multiphase model to simulate a series of benchmark cases, including concentration jump near the interface, gas dissolution in a closed system, species transport during drainage in a capillary tube, and multiphase species transport in the porous media. This CST-LB model can also be coupled with other multiphase LBMs since the model depends on the phase fraction field, which is not explicitly limited to specified multiphase models.

Pore-scale numerical simulation of low salinity water flooding using the lattice Boltzmann method

HYPOTHESIS

The change of wettability toward more water-wet by the injection of low salinity water can improve oil recovery from porous rocks, which is known as low salinity water flooding. To simulate this process at the pore-scale, we propose that the alteration in surface wettability mediated by thin water films which are below the resolution of simulation grid blocks has to be considered, as observed in experiments. This is modeled by a wettability alteration model based on rate-limited adsorption of ions onto the rock surface.

SIMULATIONS

The wettability alteration model is developed and incorporated into a lattice Boltzmann simulator which solves both the Navier-Stokes equation for oil/water two-phase flow and the advection-diffusion equation for ion transport. The model is validated against two experiments in the literature, then applied to 3D micro-CT images of a rock.

FINDINGS

Our model correctly simulated the experimental observations caused by the slow wettability alteration driven by the development of water films. In the simulations on the 3D rock pore structure, a distinct difference in the mixing of high and low salinity water is observed between secondary and tertiary low salinity flooding, resulting in different oil recoveries.

Liquid membrane catalytic model of hydrolyzing cellulose into 5-hydroxymethylfurfural based on the lattice Boltzmann method

Conversion of cellulose to 5-hydroxymethylfurfural (HMF) is an important means of biomass utilization.

Numerical Study of Surfactant Dynamics during Emulsification in a T-Junction Microchannel

Microchannel emulsification requires large amounts of surfactant to prevent coalescence and improve emulsions lifetime. However, most numerical studies have considered surfactant-free mixtures as models for droplet formation in microchannels, without taking into account the distribution of surfactant on the droplet surface. In this paper, we investigate the effects of nonuniform surfactant coverage on the microfluidic flow pattern using an extended lattice-Boltzmann model. This numerical study, supported by micro-particle image velocimetry experiments, reveals the likelihood of uneven distribution of surfactant during the droplet formation and the appearance of a stagnant cap. The Marangoni effect affects the droplet breakup by increasing the shear rate. According to our results, surfactant-free and surfactant-rich droplet formation processes are qualitatively different, such that both the capillary number and the Damköhler number should be considered when modeling the droplet generation in microfluidic devices. The limitations of traditional volume and pressure estimation methods for determining the dynamic interfacial tension are also discussed on the basis of the simulation results.

Passive Mixing inside Microdroplets

Droplet-based micromixers are essential units in many microfluidic devices for widespread applications, such as diagnostics and synthesis. The mixers can be either passive or active. When compared to active methods, the passive mixer is widely used because it does not require extra energy input apart from the pump drive. In recent years, several passive droplet-based mixers were developed, where mixing was characterized by both experiments and simulation. A unified physical understanding of both experimental processes and simulation models is beneficial for effectively developing new and efficient mixing techniques. This review covers the state-of-the-art passive droplet-based micromixers in microfluidics, which mainly focuses on three aspects: (1) Mixing parameters and analysis method; (2) Typical mixing element designs and the mixing characters in experiments; and, (3) Comprehensive introduction of numerical models used in microfluidic flow and diffusion.

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.