Lattice Boltzmann simulations for invasion patterns during drying of capillary porous media

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.

Lattice Boltzmann Modeling of Drying of Porous Media Considering Contact Angle Hysteresis

Drying of porous media is governed by a combination of evaporation and movement of the liquid phase within the porous structure. Contact angle hysteresis induced by surface roughness is shown to influence multi-phase flows, such as contact line motion of droplet, phase distribution during drainage and coffee ring formed after droplet drying in constant contact radius mode. However, the influence of contact angle hysteresis on liquid drying in porous media is still an unanswered question. Lattice Boltzmann model (LBM) is an advanced numerical approach increasingly used to study phase change problems including drying. In this paper, based on a geometric formulation scheme to prescribe contact angle, we implement a contact angle hysteresis model within the framework of a two-phase pseudopotential LBM. The capability and accuracy of prescribing and automatically measuring contact angles over a large range are tested and validated by simulating droplets sitting on flat and curved surfaces. Afterward, the proposed contact angle hysteresis model is validated by modeling droplet drying on flat and curved surfaces. Then, drying of two connected capillary tubes is studied, considering the influence of different contact angle hysteresis ranges on drying dynamics. Finally, the model is applied to study drying of a dual-porosity porous medium, where phase distribution and drying rate are compared with and without contact angle hysteresis. The proposed model is shown to be capable of dealing with different contact angle hysteresis ranges accurately and of capturing the physical mechanisms during drying in different porous media including flat and curved geometries.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11242-021-01644-9.

Analyzing Impacts of Interfacial Instabilities on the Sweeping Power of Newtonian Fluids to Immiscibly Displace Power-Law Materials

Injection of Newtonian fluids to displace pseudoplastic and dilatant fluids, governed by the power-law viscosity relationship, is common in many industrial processes. In these applications, changing the viscosity of the displaced fluid through velocity alteration can regulate interfacial instabilities, displacement efficiency, the thickness of the static wall layer, and the injected fluid’s tendency to move toward particular parts of the channel. The dynamic behavior of the fluid–fluid interface in the case of immiscibility is highly complicated and complex. In this study, a code was developed that utilizes a multi-component model of the lattice Boltzmann method to decrease the computational cost and accurately model these problems. Accordingly, a 2D inclined channel, filled with a stagnant incompressible Newtonian fluid in the initial section followed by a power-law material, was modeled for numerous scenarios. In conclusion, the results indicate that reducing the power-law index can regulate interfacial instabilities leading to dynamic deformation of static wall layers at the top and the bottom of the channel. However, it does not guarantee a reduction in the thickness of these layers, which is crucial to improve displacement efficiency. The impacts of the compatibility factor and power-law index variations on the filling pattern and finger structure were intensively evaluated.