Dynamic Pore Network Modeling of Imbibition in Real Porous Media with Corner Film Flow.

Wetting films can develop in the corners of pore structures during imbibition in a strongly wetting porous medium, which may significantly influence the two-phase flow dynamics. Due to the large difference in scales between main meniscus and corner film, accurate and efficient modeling of the dynamics of corner film remains elusive. In this work, we develop a novel two-pressure dynamic pore network model incorporating the interacting capillary bundle model to analyze the competition between the main meniscus and corner film flow in real porous media. A pore network with four-point star-shaped pore bodies and throat bonds is extracted from the real porous medium based on the pore shape factor and pore cross-sectional area, which is then decomposed into several layers of sub-pore networks, where the first layer of sub-pore network simulates the main meniscus flow while the upper layers characterize the corner film flow. The two-phase flow conductance of throat bonds for different layers of sub-pore networks are determined by high-resolution two-phase lattice Boltzmann modeling, thus inherently considering the viscous coupling effect. In addition, two artificial neural network models are developed to predict the two phases' flow conductance based on the shape of the throat cross section and the fluid properties. The accuracy of the developed model is validated with a lattice Boltzmann simulation of imbibition in a strongly wetting square tube. Then the model is used to simulate imbibition in a strongly wetting sandstone porous medium, and the competition between the main meniscus and the corner film flow is analyzed. The results show that with decreasing capillary number and viscosity ratio between wetting and nonwetting fluids, the development of the wetting corner film becomes more significant.

Pore-Scale Study on Convective Drying of Porous Media.

In this work, a numerical model for isothermal liquid-vapor phase change (evaporation) of the two-component air-water system is proposed based on the pseudopotential lattice Boltzmann method. Through the Chapman-Enskog multiscale analysis, we show that the model can correctly recover the macroscopic governing equations of the multicomponent multiphase system with a built-in binary diffusion mechanism. The model is verified based on the two-component Stefan problem where the measured binary diffusivity is consistent with theoretical analysis. The model is then applied to convective drying of a dual-porosity porous medium at the pore scale. The simulation captures a classical transition in the drying process of porous media, from the constant rate period (CRP, first phase) showing significant capillary pumping from large to small pores, to the falling rate period (FRP, second phase) with the liquid front receding in small pores. It is found that, in the CRP, the evaporation rate increases with the inflow Reynolds number (Re), while in the FRP, the evaporation curves almost collapse at different Res. The underlying mechanism is elucidated by introducing an effective Péclet number (Pe). It is shown that convection is dominant in the CRP and diffusion in the FRP, as evidenced by Pe > 1 and Pe < 1, respectively. We also find a log-law dependence of the average evaporation rate on the inflow Re in the CRP regime. The present work provides new insights into the drying physics of porous media and its direct modeling at the pore scale.

A unified lattice Boltzmann model and application to multiphase flows.

In this work, we develop a unified lattice Boltzmann model (ULBM) framework that can seamlessly integrate the widely used lattice Boltzmann collision operators, including the Bhatnagar-Gross-Krook or single-relation-time, multiple-relaxation-time, central-moment or cascaded lattice Boltzmann method and multiple entropic operators (KBC). Such a framework clarifies the relations among the existing collision operators and greatly facilitates model comparison and development as well as coding. Importantly, any LB model or treatment constructed for a specific collision operator could be easily adopted by other operators. We demonstrate the flexibility and power of the ULBM framework through three multiphase flow problems: the rheology of an emulsion, splashing of a droplet on a liquid film and dynamics of pool boiling. Further exploration of ULBM for a wide variety of phenomena would be both realistic and beneficial, making the LBM more accessible to non-specialists. This article is part of the theme issue 'Progress in mesoscale methods for fluid dynamics simulation'.

Discrete fluidization of dense monodisperse emulsions in neutral wetting microchannels.

The rheology of pressure-driven flows of two-dimensional dense monodisperse emulsions in neutral wetting microchannels is investigated by means of mesoscopic lattice Boltzmann simulations, capable of handling large collections of droplets, in the order of several hundreds. The simulations reveal that the fluidization of the emulsion proceeds through a sequence of discrete steps, characterized by yielding events whereby layers of droplets start rolling over each other, thus leading to sudden drops of the relative effective viscosity. It is shown that such discrete fluidization is robust against loss of confinement, namely it persists also in the regime of small ratios of the droplet diameter over the microchannel width. We also develop a simple phenomenological model which predicts a linear relation between the relative effective viscosity of the emulsion and the product of the confinement parameter (global size of the device over droplet radius) and the viscosity ratio between the disperse and continuous phases. The model shows excellent agreement with the numerical simulations. The present work offers new insights to enable the design of microfluidic scaffolds for tissue engineering applications and paves the way to detailed rheological studies of soft-glassy materials in complex geometries.

Unified lattice Boltzmann method with improved schemes for multiphase flow simulation: Application to droplet dynamics under realistic conditions.

As a powerful mesoscale approach, the lattice Boltzmann method (LBM) has been widely used for the numerical study of complex multiphase flows. Recently, Luo et al. [Philos. Trans. R. Soc. A: Math. Phys. Eng. Sci. 379, 20200397 (2021)10.1098/rsta.2020.0397] proposed a unified lattice Boltzmann method (ULBM) to integrate the widely used lattice Boltzmann collision operators into a unified framework. In this study, we incorporate additional features into this ULBM in order to simulate multiphase flow under realistic conditions. A nonorthogonal moment set [Fei et al., Phys. Rev. E 97, 053309 (2018)10.1103/PhysRevE.97.053309] and the entropic-multi-relaxation-time (KBC) lattice Boltzmann model are used to construct the collision operator. An extended combined pseudopotential model is proposed to realize multiphase flow simulation at high-density ratio with tunable surface tension over a wide range. The numerical results indicate that the improved ULBM can significantly decrease the spurious velocities and adjust the surface tension without appreciably changing the density ratio. The ULBM is validated through reproducing various droplet dynamics experiments, such as binary droplet collision and droplet impingement on superhydrophobic surfaces. Finally, the extended ULBM is applied to complex droplet dynamics, including droplet pancake bouncing and droplet splashing. The maximum Weber number and Reynolds number in the simulation reach 800 and 7200, respectively, at a density ratio of 1000. The study demonstrates the generality and versatility of ULBM for incorporating schemes to tackle challenging multiphase problems.

Droplet evaporation in finite-size systems: Theoretical analysis and mesoscopic modeling.

The classical D^{2}-Law states that the square of the droplet diameter decreases linearly with time during its evaporation process, i.e., D^{2}(t)=D_{0}^{2}-Kt, where D_{0} is the droplet initial diameter and K is the evaporation constant. Though the law has been widely verified by experiments, considerable deviations are observed in many cases. In this work, a revised theoretical analysis of the single droplet evaporation in finite-size open systems is presented for both two-dimensional (2D) and 3D cases. Our analysis shows that the classical D^{2}-Law is only applicable for 3D large systems (L≫D_{0}, L is the system size), while significant deviations occur for small (L≤5D_{0}) and/or 2D systems. Theoretical solution for the temperature field is also derived. Moreover, we discuss in detail the proper numerical implementation of droplet evaporation in finite-size open systems by the mesoscopic lattice Boltzmann method (LBM). Taking into consideration shrinkage effects and an adaptive pressure boundary condition, droplet evaporation in finite-size 2D/3D systems with density ratio up to 328 within a wide parameter range (K=[0.003,0.18] in lattice units) is simulated, and remarkable agreement with the theoretical solution is achieved, in contrast to previous simulations. The present work provides insights into realistic droplet evaporation phenomena and their numerical modeling using diffuse-interface methods.

Effects of hysteresis window on contact angle hysteresis behaviour at large Bond number.

Contact angle hysteresis, defined as the difference between advancing and receding contact angles, is an important phenomenon in multiphase flow on a wetting surface. In this study, a modified pseudo-potential lattice Boltzmann (LB) multiphase model with tunable surface tension is proposed, which is further coupled with the geometrical formulation contact angle scheme to investigate the motion of droplets invoking the contact angle hysteresis. We focus on the dynamic behaviour of droplets driven by a body force at the Bond number ranging from 1 to 6, which is defined as the ratio of the body force to the capillary force. The droplet morphology change is examined by varying (i) the Bond number and (ii) the hysteresis window. Results show the droplet morphology evolution can be classified into different stages, including stretch, relaxation, and equilibrium. The droplet oscillation phenomenon at large Bond numbers at the equilibrium stage is observed for the first time. In addition, it is found that such oscillation can lead to the breakup and/or coalescence of droplets when the surface waves spread on the top of the droplet. Furthermore, there is slight oscillation of the normalized length, width and height at the equilibrium stage for the neutral hysteresis window while more dramatic oscillation will appear for the hydrophobic hysteresis window.

Three-dimensional cascaded lattice Boltzmann method: Improved implementation and consistent forcing scheme.

The cascaded or central-moment-based lattice Boltzmann method (CLBM) proposed in [Phys. Rev. E 73, 066705 (2006)PLEEE81539-375510.1103/PhysRevE.73.066705] possesses very good numerical stability. However, two constraints exist in three-dimensional (3D) CLBM simulations. First, the conventional implementation for 3D CLBM involves cumbersome operations and requires much higher computational cost compared to the single-relaxation-time (SRT) LBM. Second, it is a challenge to accurately incorporate a general force field into the 3D CLBM. In this paper, we present an improved method to implement CLBM in 3D. The main strategy is to adopt a simplified central moment set and carry out the central-moment-based collision operator based on a general multi-relaxation-time (GMRT) framework. Next, the recently proposed consistent forcing scheme for CLBM [Fei and Luo, Phys. Rev. E 96, 053307 (2017)2470-004510.1103/PhysRevE.96.053307] is extended to incorporate a general force field into 3D CLBM. Compared with the recently developed nonorthogonal CLBM [Rosis, Phys. Rev. E 95, 013310 (2017)2470-004510.1103/PhysRevE.95.013310], our implementation is proved to reduce the computational cost significantly. The inconsistency of adopting the discrete equilibrium distribution functions in the nonorthogonal CLBM is analyzed and validated. The 3D CLBM developed here in conjunction with the consistent forcing scheme is verified through numerical simulations of several canonical force-driven flows, highlighting very good properties in terms of accuracy, convergence, and consistency with the nonslip rule. Finally, the techniques developed here for 3D CLBM can be applied to make the implementation and execution of 3D MRT-LBM more efficient.

Consistent forcing scheme in the cascaded lattice Boltzmann method.

In this paper, we give an alternative derivation for the cascaded lattice Boltzmann method (CLBM) within a general multiple-relaxation-time (MRT) framework by introducing a shift matrix. When the shift matrix is a unit matrix, the CLBM degrades into an MRT LBM. Based on this, a consistent forcing scheme is developed for the CLBM. The consistency of the nonslip rule, the second-order convergence rate in space, and the property of isotropy for the consistent forcing scheme is demonstrated through numerical simulations of several canonical problems. Several existing forcing schemes previously used in the CLBM are also examined. The study clarifies the relation between MRT LBM and CLBM under a general framework.

A multi-component discrete Boltzmann model for nonequilibrium reactive flows.

We propose a multi-component discrete Boltzmann model (DBM) for premixed, nonpremixed, or partially premixed nonequilibrium reactive flows. This model is suitable for both subsonic and supersonic flows with or without chemical reaction and/or external force. A two-dimensional sixteen-velocity model is constructed for the DBM. In the hydrodynamic limit, the DBM recovers the modified Navier-Stokes equations for reacting species in a force field. Compared to standard lattice Boltzmann models, the DBM presents not only more accurate hydrodynamic quantities, but also detailed nonequilibrium effects that are essential yet long-neglected by traditional fluid dynamics. Apart from nonequilibrium terms (viscous stress and heat flux) in conventional models, specific hydrodynamic and thermodynamic nonequilibrium quantities (high order kinetic moments and their departure from equilibrium) are dynamically obtained from the DBM in a straightforward way. Due to its generality, the developed methodology is applicable to a wide range of phenomena across many energy technologies, emissions reduction, environmental protection, mining accident prevention, chemical and process industry.

[Genetic diversity of human Parvovirus B19 VP1 unique region].

OBJECTIVE: Human Parvovirus B19 (HPV B19) is a small (23 nm), non-enveloped DNA virus found in 1974. It has been proved that HPV B19 is associated with a variety of childhood diseases, such as erythema infectious, transient aplastic crisis, aplastic anemia, idiopathic thrombocytopenic purpura and arthropathy, etc. There have been no any effective vaccines to prevent HPV B19 infection so far. The HPV B19 genome is composed of 5.6 kb single strand DNA. This genome encodes a nonstructural protein NS1, two structural proteins VP1 and VP2. Most neutralizing linear epitopes of HPV B19 cluster in the VP1 unique and VP1-VP2 junction regions. Only proteins encoded by genes of the VP1 unique and VP1-VP2 junction regions can stimulate bodies to produce protective antibodies. Aim of the present study was to get the VP1 unique region gene of HPV B19 and to analyze the genetic diversity so as to further study its function and application. METHODS: The VP1 unique region gene of HPV B19 was amplified from the serum of a child with idiopathic thrombocytopenic purpura by PCR. The purified PCR product was cloned into pGEM-T easy vector and transfected into the host strain E. coli (DH5 alpha). Positive clones were chosen and then the target gene was sequenced. RESULTS: The target gene sequence of HPV B19 VP1 unique region was amplified and cloned successfully. It had 705 nucleotides. Compared with the relevant sequences published in Genbank, the sequencing results were revealed with two nucleotides changes in the HPV B19 VP1 unique region, but their coding amino acid were not changed. CONCLUSION: It is suggested that genetic diversity exists in the VP1 unique region of HPV B19. Construction of the recombinant plasmid of HPV B19 VP1 unique region gene might benefit to further study.