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

Visualization of the Sol-Gel Transition in Porous Networks Using Fluorescent Viscosity-Sensitive Probes

The sol-gel transition involves the transformation of a colloidal suspension into a system-spanning, interconnected gel. This process is widely used to reinforce mechanically weakened porous artifacts, such as sculptures but the impact of the restricted geometry (porous network) on the gelation dynamics of the solution remains unclear. Here, using fluorescent viscosity-sensitive molecular rotors, confocal microscopy, and model pores, we visualize the local viscosity changes at the microscale that accompany the sol-gel transition of a methyltriethoxysilane solution into a gel network. We show that, with evaporation of the solvent, a viscosity gradient develops near the free surface, triggering the sol-gel transition inside small pores near the surface. In homogeneous porous media, this leads to skin formation, which reduces the evaporation rate. In heterogeneous porous media, a gradient in gel density develops toward the heart of the porous material, where the gel formation mainly occurs as capillary bridges within smaller pores.

Is Unidirectional Drying in a Round Capillary Always Diffusive?

The unidirectional drying of water in cylindrical capillaries has been described since the discovery of Stefan's solution as a vapor diffusion-controlled process with a square root of time kinetics. Here we show that this well-known process actually depends on the way the capillary is closed. Experiments are performed on the evaporation of water in capillaries closed at one end with a solid material or connected to a fluid reservoir. While we recover Stefan's solution in the first case, we show that in the second situation the water plug evaporates at a constant rate with the water-air meniscus remaining pinned at the exit where evaporation proceeds. The presence of the liquid reservoir closing the capillary combined with a capillary pumping effect induces a flow of the water plug toward the evaporation front leading to a constant-rate drying, substantially faster than the prediction of Stefan's equation. Our results show that a transition from a constant-rate evaporation regime at short times to a diffusion-driven evaporation regime at long times can be observed by increasing the viscosity of the fluid in the reservoir blocking the other end of the capillary. Such transition can also be observed by connecting the capillary end to a solidifying fluid like epoxy glue.

Wicking fingering in electrospun membranes

Pronounced fingering of the waterfront is observed for in-plane wicking in thin, aligned electrospun fibrous membranes. We hypothesize that a perturbation in capillary pressure triggers the onset of fingering, which grows in a non-local manner based on the waterfront gradient. Vertical and horizontal wicking in thin electrospun membranes of poly(ethylene-co-vinyl alcohol) (EVOH) fibers with varying fiber alignment and degree of orientation is studied with backlight photography. A non-local transport model considering the gradient of the waterfront is developed, where fiber orientation is modeled with a correlated random field. The model shows that a transition from straight to highly fingered waterfront occurs during water uptake as observed in the experiment. The size and shape of the fingers depend on fiber orientation. Based on good model agreement, we show that, during wicking in thin electrospun membranes, fingering is initially triggered by a perturbation in capillary pressure caused by the underlying anisotropic and heterogeneous membrane structure which grows in a non-local manner depending on the waterfront gradient.

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.

Implementing Contact Angle Hysteresis in Moving Mesh-Based Two-Phase Flow Numerical Simulations

Contact angle hysteresis is a common phenomenon in nature, which also plays an important role in industrial applications. A numerical model based on the moving mesh two-phase flow method is presented for modeling contact angle hysteresis. The implementation includes a displacement-based penalty method and a state variable method. The pinning, moving, and repinning of the contact lines can be simulated. This method is robust considering both two-dimensional and three-dimensional geometries. To further demonstrate the performance of this method, a fluid-solid interaction model with a cylinder fluctuating on a water surface considering contact angle hysteresis is demonstrated.

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

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. This paper 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 nonlocal 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.