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

Qin, Feifei; Zhao, Jianlin; Kang, Qinjun; Derome, Dominique; Carmeliet, Jan

doi:10.1007/s11242-021-01644-9

Cited by 29 publications

(39 citation statements)

References 84 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Considering the important role of wettability on drying patterns in porous media, , we analyze the implementation of an accurate contact angle. In 2007, Ding et al showed that the contact angle in complex geometries in the phase-field method can be accurately implemented based on a geometric-function scheme .…”

Section: Model Validationmentioning

confidence: 99%

Pore-Scale Study on Convective Drying of Porous Media

et al. 2022

Self Cite

View full text Add to dashboard Cite

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 Pećlet 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.

show abstract

Section: Model Validationmentioning

confidence: 99%

Pore-Scale Study on Convective Drying of Porous Media

et al. 2022

Self Cite

View full text Add to dashboard Cite

show abstract

“…This should impact the fluidfluid surface area, which in turn will have impact on the mass transfer process in the interfacial region. A possibility to explore this issue could be to rely on some direct numerical simulations [44,45] on microstructure digital images using the results of the present article as guidelines.…”

Section: Discussionmentioning

confidence: 99%

Coupling between internal and external mass transfer during stage-1 evaporation in capillary porous media: Interfacial resistance approach

Talbi

Prat

2021

Phys. Rev. E

View full text Add to dashboard Cite

The coupling boundary condition to be imposed at the evaporative surface of a porous medium is studied from pore network simulations considering the capillary regime. The study highlights the formation of a thin edge effect region of smaller saturation along the evaporative surface. It is shown that this thin region forms in the breakthrough period at the very beginning of the drying process. The size of this region is studied and shown to be not network size dependent. This region is shown to be the locus of a non-local equilibrium effect. The features lead to the consideration of a coupling boundary condition involving an interfacial mass transfer resistance and an external mass transfer resistance. Contrary to previous considerations, it is shown that both resistances depend on the variation of the saturation, i.e. the fluid topology, and the size of the external mass transfer layer, i.e. the mass transfer rate. This is explained by the evolution of the vapor partial pressure distribution at the surface which becomes increasingly heterogeneous during evaporation and depends on both the evolving fluid distribution in the interfacial region and the mass transfer rate. However, the geometric effects due to the configuration of the fluids can be separated from rate effects that arise due to the nonequilibrium mass transport.

show abstract

“…The implementation of sub-resolution processes, e.g. contact angle hysteresis on rough surfaces, 23,24 is an ongoing challenge for coarse-grain-like models. Theoretical and experimental works on hysteresis and wetting of rough surfaces often discuss the superposition of different effects, e.g.…”

Section: Introductionmentioning

confidence: 99%