Search

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.

Nanoporous adsorbents can mechanically swell or shrink once upon the accumulation of guest fluid molecules at their internal surfaces or in their cavities. Existing theories in this field attribute such sorption-induced swelling to a tensile force, while shrinkage is always associated with a contractive force. In this study, however, we propose that the sorption-induced deformation of a porous architecture is not solely dictated by the stress conditions but can also be largely influenced by its mechanical anisotropy. In more detail, the sorption-induced deformation of a polymeric slab is investigated using a hybrid molecular dynamics and Monte Carlo algorithm. When subjected to water loading, the slab is found to swell along its normal direction and display an overall positive volumetric strain. Moreover, the surface roughness is enhanced as a response to the surface energy decrease induced by the water covering the slab external surface. Unexpectedly, the in-plane deformation of the slab material seems to be highly constrained, so that it is far below its normal counterpart. This anisotropy is enhanced when the slab thickness decreases. With a thickness of around 1.35 nm, an in-plane shrinkage is observed throughout the entire hygroscopic range. A theoretical analysis based on a poromechanical model suggests that the anisotropic mechanical properties, which are common for a slab material, are the essence of the constrained in-plane swelling or even shrinkage under the isotropic sorption-induced tensile forces. This study, unveiling overlooked mechanisms of sorption-induced shrinkage in mechanically anisotropic materials, provides new insights into this field.

Prevailing absorbents like wood-derived porous scaffolds or polymeric aerogels are normally featured with hierarchical porous structures. In former molecular simulation studies, sorption, deformation, and coupled sorption-deformation have been studied for single-scale materials, but scarcely for materials where micropores (<2 nm) and mesopores (2-50 nm) coexist. The present work, dealing with a mesoscopic slit pore between two slabs of microporous amorphous cellulose (AC), aims at modeling sorption-deformation interplay in hierarchical porous cellulosic structures inspired by polymeric modern adsorbents. Specifically, the atomic system is modeled by a hybrid workflow combining molecular dynamics (MD) and grand canonical Monte Carlo (GCMC) simulations. The results clarify the multiple sorption/deformation mechanisms in porous materials with different slit-pore sizes, including water filling in micropores, surface covering at the solid-air interface, and subsequent capillary condensation in mesopores. In particular, before the onset of capillary condensation, the sorption behavior of the AC matrix in the hybrid system is almost the same as that of bulk AC, in which sorption and deformation enhance each other through sorption-induced swelling and additional sorption in the newly created voids. Upon capillary condensation, the interaction between the micropores and the mesopore emerges. Water molecules in the mesopore exert a negative hydrostatic pressure perpendicular to the slab surface on the matrices, resulting in an increase in porosity and water content, a decrease in distance between the centers of mass (COMs) of the slabs, and thus a thinning of the slit pore. As described by Bangham's Law, the surface area of the rough slit-pore slab increases proportionally to the surface energy variation during surface covering. For a system composed of a compliant polymer like AC, however, the surface area enlargement does not result in an in-plane swelling as expected but instead in an in-plane shrinkage along with an increase in local roughness or irregularity (an accordion effect).

Diffusion of molecules in porous media is a critical process that is fundamental to numerous chemical, physical, and biological applications. The prevailing theoretical frameworks are challenged when explaining the complex dynamics resulting from the highly tortuous host structure and strong guest-host interactions, especially when the pore size approximates the size of diffusing molecule. This study, using molecular dynamics, formulates a semiempirical model based on theoretical considerations and factorization that offer an alternative view of diffusion and its link with the structure and behavior (sorption and deformation) of material. By analyzing the intermittent dynamics of water, microscopic self-diffusion coefficients are predicted. The apparent tortuosity, defined as the ratio of the bulk to the confined self-diffusion coefficients, is found to depend quantitatively on a limited set of material parameters: heat of adsorption, elastic modulus, and percolation probability, all of which are experimentally accessible. The proposed sorption-deformation-percolation model provides guidance on the understanding and fine-tuning of diffusion.

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.

When dealing with reactions of a liquid reactant and a solid catalyst, macroreactors with vigorous stirring equipment may be dangerous and cause wastage of energy. Reducing the diffusion distance and promoting reactants to reach the catalyst surface for efficient reaction remain the key challenges. Here, inspired by capillary-driven water motion in plants, we propose to implement a self-driven multiplex reaction (SMR) in nanocatalyst-loaded microchannels. Unlike the classical capillary rise, the droplet in SMR has variable pressure difference, leading to tunable flow velocity for controlling the reaction rate without any auxiliary equipment. The SMR in microchannels contributes to an increase in the reaction rate by more than 2 orders of magnitude compared to that in macroreactors. Specifically, this strategy reduces the reaction volume by 170 times, the catalyst usage by about 12 times, and the energy consumption by 50 times. This apparatus with a small volume and less catalyst content promises to provide an efficient strategy for the precise manipulation of chemical reactions.

Sorption hysteresis in nanoporous polymer is an intriguing phenomenon that involves coupling between sorption and deformation. Based on the mechanism revealed at the microscopic level by use of molecular simulation, a poromechanical model is developed capturing all relevant physics and yielding a quantitative description. In this model, the coupling between sorption and deformation is described by a poromechanics framework. More in detail, an upscaling process from the molecular mechanism is implemented to model the hysteresis through the state change of each element upon deformation. We provide two solutions of the model: a numerical one based on the finite element method and an analytical one based on uniform strain assumption. The results from both solutions agree well with the molecular simulation and experimental results, therefore capturing and describing adequately sorption hysteresis. The developed model illustrates that water forms different structural distributions upon adsorption and desorption. A parametric study shows that sorption hysteresis is influenced by material properties. We find that a softer material with stronger adsorbent-adsorbate interaction tends to exhibit more profound sorption hysteresis. The developed model, which relies on the concepts of sorption-deformation coupling and multiscale modeling from atomistic simulations to domain dependent theory, paves the way for a new direction of modeling sorption hysteresis.

An atomistic slit pore model is built to study the sorption-induced deformation of nanoporous materials with the help of molecular simulation. Both sorption and strain isotherms are determined to probe the anisotropic deformation behavior induced upon molecular adsorption. A detailed analysis shows that the driving microscopic mechanisms at different sorption stages are different. At high relative pressure, as expected from the classical macroscopic picture, the pore deformation is governed by the Laplace pressure as the pore gets filled with liquid because of capillary condensation. In such situation, the strain in normal and longitudinal directions can be predicted from the stiffness modulus in the corresponding direction. At low pressure, when liquid films are adsorbed at the pore surfaces and separated by the vapor phase, the strain is driven by the attractive solid-fluid forces and in-plane pressure within the film, and the deformation is confined in the direction parallel to the film-solid interface. Because of the interplay of the two factors, the strain changes from shrinkage to expansion upon increase of pressure. Analysis of isosteric heat of adsorption shows that the contribution arising from the deformation is small compared to the sorption contribution, which indicates that the influence of deformation on the sorption process is limited.

The contact line dynamics of evaporating droplets deposited on a set of parallel microribs is analyzed with the use of a recently developed entropic lattice Boltzmann model for two-phase flow. Upon deposition, part of the droplet penetrates into the space between ribs because of capillary action, whereas the remaining liquid of the droplet remains pinned on top of the microribs. In the first stage, evaporation continues until the droplet undergoes a series of pinning-depinning events, showing alternatively the constant contact radius and constant contact angle modes. While the droplet is pinned, evaporation results in a contact angle reduction, whereas the contact radius remains constant. At a critical contact angle, the contact line depins, the contact radius reduces, and the droplet rearranges to a larger apparent contact angle. This pinning-depinning behavior goes on until the liquid above the microribs is evaporated. By computing the Gibbs free energy taking into account the interfacial energy, pressure terms, and viscous dissipation due to drop internal flow, we found that the mechanism that causes the unpinning of the contact line results from an excess in Gibbs free energy. The spacing distance and the rib height play an important role in controlling the pinning-depinning cycling, the critical contact angle, and the excess Gibbs free energy. However, we found that neither the critical contact angle nor the maximum excess Gibbs free energy depends on the rib width. We show that the different terms, that is, pressure term, viscous dissipation, and interfacial energy, contributing to the excess Gibbs free energy, can be varied differently by varying different geometrical properties of the microribs. It is demonstrated that, by varying the spacing distance between the ribs, the energy barrier is controlled by the interfacial energy while the contribution of the viscous dissipation is dominant if either rib height or width is changed. Main finding of this is study is that, for microrib patterned surfaces, the energy barrier required for the contact line to depin can be enlarged by increasing the spacing or the rib height, which can be important for practical applications.

Droplet wetting and distortion on flat surfaces with heterogeneous wettability are studied using the 3D Shan-Chen pseudopotential multiphase lattice Boltzmann model (LBM). The contact angles are compared with the Cassie mode, which predicts an apparent contact angle for flat surfaces with different wetting properties, where the droplet size is large compared to the size of the heterogeneity. In this study, the surface studied consists in a regular checkboard pattern with two different Young's contact angles (hydrophilic and hydrophobic) and equal surface fraction. The droplet size and patch size of the checkboard are varied beyond the limit where Cassie's equation is valid. A critical ratio of patch size to droplet radius is found below which the apparent contact angle follows the Cassie mode. Above the critical value, the droplet shape changes from a spherical cap to a more distorted form, and no single contact angle can be determined. The local contact angles are found to vary along the contact line between minimum and maximum values. The droplet is found to wet preferentially the hydrophilic region, and the wetted area fraction of the hydrophilic region increases quasi-linearly with the ratio between patch and droplet sizes. We propose a new equation beyond the critical ratio, defining an equivalent contact angle, where the wetted area fractions are used as weighting coefficients for the maximum and minimum local contact angles. This equivalent contact angle is found to equal Cassie's contact angle.

Separating petroleum hydrocarbons from water is an important problem to address in order to mitigate the disastrous effects of hydrocarbons on aquatic ecosystems. A rational approach to address the problem of marine oil-water separation is to disperse the oil with the aid of surfactants in order to minimize the formation of large slicks at the water surface and to maximize the oil-water interfacial area. Here we investigate the fundamental wetting and transport behavior of such surfactant-stabilized droplets and the flow conditions necessary to perform sieving and separation of these stabilized emulsions. We show that, for water-soluble surfactants, such droplets are completely repelled by a range of materials (intrinsically underwater superoleophobic) due to the detergency effect; therefore, there is no need for surface micro-/nanotexturing or chemical treatment to repel the oil and prevent fouling of the filter. We then simulate and experimentally investigate the effect of emulsion flow rate on the transport and impact behavior of such droplets on rigid meshes to identify the minimum pore opening (w) necessary to filter a droplet with a given diameter (d) in order to minimize the pressure drop across the mesh-and therefore maximize the filtering efficiency, which is strongly dependent on w. We define a range of flow conditions and droplet sizes where minimum droplet deformation is to be expected and therefore find that the condition of w ≈ d is sufficient for efficient separation. With this new understanding, we demonstrate the use of a commercially available filter-without any additional surface engineering or functionalization-to separate oil droplets (d < 100 μm) from a surfactant-stabilized emulsion with a flux of ∼11,000 L m -2 h -1 bar -1 . We believe these findings can inform the design of future oil separation materials.

Droplet impact has been imaged on different rigid, smooth, and rough substrates for three liquids with different viscosity and surface tension, with special attention to the lower impact velocity range. Of all studied parameters, only surface tension and viscosity, thus the liquid properties, clearly play a role in terms of the attained maximum spreading ratio of the impacting droplet. Surface roughness and type of surface (steel, aluminum, and parafilm) slightly affect the dynamic wettability and maximum spreading at low impact velocity. The dynamic contact angle at maximum spreading has been identified to properly characterize this dynamic spreading process, especially at low impact velocity where dynamic wetting plays an important role. The dynamic contact angle is found to be generally higher than the equilibrium contact angle, showing that statically wetting surfaces can become less wetting or even nonwetting under dynamic droplet impact. An improved energy balance model for maximum spreading ratio is proposed based on a correct analytical modeling of the time at maximum spreading, which determines the viscous dissipation. Experiments show that the time at maximum spreading decreases with impact velocity depending on the surface tension of the liquid, and a scaling with maximum spreading diameter and surface tension is proposed. A second improvement is based on the use of the dynamic contact angle at maximum spreading, instead of quasi-static contact angles, to describe the dynamic wetting process at low impact velocity. This improved model showed good agreement compared to experiments for the maximum spreading ratio versus impact velocity for different liquids, and a better prediction compared to other models in literature. In particular, scaling according to We(1/2) is found invalid for low velocities, since the curves bend over to higher maximum spreading ratios due to the dynamic wetting process.

The maximum spreading of an impinging droplet on a rigid surface is studied for low to high impact velocity, until the droplet starts splashing. We investigate experimentally and numerically the role of liquid properties, such as surface tension and viscosity, on drop impact using three liquids. It is found that the use of the experimental dynamic contact angle at maximum spreading in the Kistler model, which is used as a boundary condition for the CFD-VOF calculation, gives good agreement between experimental and numerical results. Analytical models commonly used to predict the boundary layer thickness and time at maximum spreading are found to be less correct, meaning that energy balance models relying on these relations have to be considered with care. The time of maximum spreading is found to depend on both the impact velocity and surface tension, and neither dependency is predicted correctly in common analytical models. The relative proportion of the viscous dissipation in the total energy budget increases with impact velocity with respect to surface energy. At high impact velocity, the contribution of surface energy, even before splashing, is still substantial, meaning that both surface energy and viscous dissipation have to be taken into account, and scaling laws depending only on viscous dissipation do not apply. At low impact velocity, viscous dissipation seems to play an important role in low-surface-tension liquids such as ethanol.

The diffusion of H2O in three amorphous polymer-H2O systems is studied as a function of H2O content using molecular dynamics. A picture of H2O molecule motion comprising alternating steps of being bound at an adsorption site ("stop") and moving ("go") emerges. This picture is made quantitative. The bound time, frequency of stop-go steps, and tortuosity all decrease with H2O content. Fourier analysis of particle motion during bound time segments provides a measure of an attempt frequency that is connected quantitatively to the bound time and an activation energy of a hydrogen bond. For increasing H2O content, the polymer-H2O systems swell, leading to an increase in the diffusion coefficient and porosity and a decrease in activation energy.

A two-phase model of a wood microfibril consisting of crystalline cellulose and amorphous hemicellulose is investigated with molecular dynamics in full range of sorption to understand the molecular origin of swelling and weakening of wood. Water is adsorbed in hemicellulose, and an excess of sorption is found at the interface, while no sorption occurs within cellulose. Water molecules adsorbed on the interface push away polymer chains, forcing the two phases to separate and causing breaking of h-bonds, particularly pronounced on the interface. Existence of two different regions in moisture response is demonstrated. At low moisture content, water is uniformly adsorbed within hemicellulose, breaking a small amount of hydrogen bonds. Microfibril does not swell, and the porosity does not change. As moisture content increases, water is adsorbed preferentially at the interface, which leads to additional swelling and porosity increase at the interface. Young's and shear moduli decrease importantly due to breaking of h-bonds and screening of the long-range interactions.

We investigate the influence of adsorbed water on amorphous cellulose structure and properties, within the full range of moisture content from the dry state to saturation, by molecular dynamics simulation. Increasing water content results in overall swelling, a substantial decrease in stiffness, and higher diffusivity of the water molecules. The obtained sorption curve as well as the range of swelling and weakening are confirmed by experiments. The measured properties undergo a noticeable change at about 10% of moisture content, which suggests that a transition occurs in the porous system, indicating that the sorption process is stepwise. Our analysis of water network formation reveals that the onset of percolation coincides with the moisture content at which a transition in the material properties is observed. An in-depth analysis of the molecular mechanism of hydrogen bonding, van der Waals interactions, and water network in the two regimes enhances the understanding of the adsorption process.