An official website of the United States government
Abstract We present a new computational fluid dynamics approach for simulating two‐phase flow in hybrid systems containing solid‐free regions and deformable porous matrices. Our approach is based on the derivation of a unique set of volume‐averaged partial differential equations that asymptotically approach the Navier‐Stokes Volume‐of‐Fluid equations in solid‐free regions and multiphase Biot Theory in porous regions. The resulting equations extend our recently developed Darcy‐Brinkman‐Biot framework to multiphase flow. Through careful consideration of interfacial dynamics (relative permeability and capillary effects) and extensive benchmarking, we show that the resulting model accurately captures the strong two‐way coupling that is often exhibited between multiple fluids and deformable porous media. Thus, it can be used to represent flow‐induced material deformation (swelling, compression) and failure (cracking, fracturing). The model's open‐source numerical implementation,hybridBiotInterFoam, effectively marks the extension of computational fluid mechanics into modeling multiscale multiphase flow in deformable porous systems. The versatility of the solver is illustrated through applications related to material failure in poroelastic coastal barriers and surface deformation due to fluid injection in poro‐visco‐plastic systems.
more »
« less
This content will become publicly available on November 1, 2026
Interfacial Evolution Explains the Complex Swelling-Shrinkage Responses of Porous Materials from Vacuum-Dry to Full Liquid Saturation
Adsorption-induced swelling occurs in a wide spectrum of natural and engineered porous materials. A key underlying mechanism is the monotonic reduction of solid-fluid surface energy upon fluid adsorption, which lowers the contractive adsorption stress and causes the porous skeleton to swell (Bangham and Fakhoury, 1928). Some mesoporous materials, however, deviate from the monotonic swelling pattern predicted by this mechanism, exhibiting an abrupt shrinkage at intermediate adsorbate partial pressures before swelling resumes and continues to full saturation. This behavior is commonly attributed to capillary condensation of the adsorbate from the vapor to the liquid phase within the pores. Understanding the stresses and the shrinkage induced by capillary condensation is critical in various industrial applications including micro-/nanofabrication, geotechnical engineering in collapsible soils, and sorption-driven actuation technologies. This work aims to develop a unified poromechanics theory that captures the full sequence of adsorption-induced deformation, including initial swelling, contraction during capillary condensation, and resumed expansion near full saturation. The formulation begins with a thermodynamic analysis of an unsaturated deformable porous solid acknowledging the energetics of the solid-fluid (sl), solid-vapor (sv), and liquid-vapor (lv) interfaces. The resulting free energy balance permits the simultaneous derivation of the liquid retention characteristics curve and the coupled mechanical effects driven by adsorption and partial saturation. Within this framework, two strategies for constructing constitutive relations are examined: one explicitly resolves the dynamic evolution of sl-sv-lv interfacial areas to emphasize the underlying physics, while the other partially lumps the surface energies into a macroscopic capillary potential to facilitate model calibration using standard laboratory tests. The models are evaluated using datasets from two markedly different solid-fluid systems: N2 gas adsorption on a hierarchical porous silica at 77 K and water adsorption on a carbon xerogel at 298 K. Both approaches effectively capture the complex, non-monotonic strain isotherms exhibited by the adsorbent. The adsorption-desorption hysteresis is also addressed in a thermodynamically consistent framework. The proposed theory demonstrates both robustness and unifying power in explaining the complex strain isotherms of porous materials along adsorption and desorption paths, covering the entire spectrum from vacuum-dry to fully liquid-saturated states.
more »
« less
- PAR ID:
- 10649088
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- Journal of the Mechanics and Physics of Solids
- ISSN:
- 0022-5096
- Page Range / eLocation ID:
- 106425
- Subject(s) / Keyword(s):
- adsorption capillary condensation poromechanics unsaturated porous media surface energy hysteresis
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Strongly bound surface species like alkylamines adsorbed on the Brønsted acid site of aluminosilicate zeolites exhibit negligible rates of molecular desorption, preventing them from achieving an equilibrated state on experimentally relevant timescales that limit the measurement of their adsorption thermodynamics. Through adsorption-assisted desorption, whereby distinct alkylamines facilitate desorption from Brønsted acid sites, we demonstrate that equilibrated states are achieved. Breakthrough adsorption measurements reveal that while 2-butylammonium on a Brønsted acid site is irreversibly adsorbed, it readily undergoes molecular desorption when exposed to a distinct alkylamine like 2-propanamine. As a result, two-adsorbate equilibrium was achieved when exposing Brønsted acid sites of aluminosilicate zeolites to a binary vapor phase alkylamine mixture. By varying relative vapor phase partial pressures and temperatures, we demonstrate the ability to experimentally measure the adsorption enthalpy and entropy of alkylammonium adsorbates on mostly isolated Brønsted acid sites in H-ZSM-5 (Si/Al = 140). A multi-adsorbate Langmuir isotherm was found to quantitatively describe the co-adsorption of alkylamines varying in size and basicity over a wide range of conditions, through which the relative adsorption enthalpy and entropy of alkylamines were measured. Across a homologous family of sec-alkylamines (C3-C5) adsorbed on isolated Brønsted acid sites, a fixed contribution to the enthalpy (19 ± 4 kJ mol CH2-1) and entropy (25 ± 4 J mol CH2-1 K-1) of adsorption per methylene unit of was found to exist, likely resulting from electrostatic interactions between the alkyl chain and surrounding pore environment.more » « less
-
The aim of this work is to enhance the understanding of the pore structure and adsorption properties of kerogens as applied to organic-rich shales and mudstone rocks. Conventional methods of adsorption characterization from low temperature N2 isotherms rely on the use of the so-called standard isotherms on nonporous substrates (typically silica or amorphous carbons), which may not be accurate for the surfaces of kerogens. In this work, we present a new methodology for pore size characterization of kerogens that relies on a realistic molecular model of kerogen surfaces. Taking advantage of recent advances in modeling the molecular structure of kerogens, we create atomistic three-dimensional (3D) models of amorphous bulk kerogens, rough kerogen surfaces, and mesopores imbedded in the amorphous kerogen matrix. Using grand canonical Monte Carlo (GCMC) simulations, we calculate the reference N2 adsorption isotherms in the micropores of the bulk kerogen matrix, on the kerogen surface, as well as in a series of mesopores confined by rough kerogen walls. Next, we parameterized the quenched solid density functional theory (QSDFT) to reproduce the kerogen surface heterogeneity and GCMC-simulated N2 adsorption isotherms. Furthermore, we approximated the isotherm on the reference kerogen surface by a macroscopic disjoining pressure isotherm, which allows us to use the Derjaguin−Broekhoff−de Boer (DBdB) model to predict adsorption and capillary condensation in meso/macropores. The reference GCMC, QSDFT, and DBdB isotherms are combined into the kernel for calculating the micropore volume, meso- and macropore surfaces, and mesopore size distribution from the experimental adsorption isotherms. The proposed methodology is demonstrated on a typical example of a kerogen II-A sample with a wide mesopore size distribution. The methodology can be extended to other kerogen structures of different maturities to provide a comprehensive characterization of organic porosity in kerogen fractions.more » « less
-
Abstract Adsorption-induced deformation is a change in geometrical dimensions of an adsorbent material caused by gas or liquid adsorption on its surface. This phenomenon is universal and sensitive to adsorbent properties, which makes its prediction a challenging task. However, the pure academic interest is complemented by its importance in a number of engineering applications with porous materials characterization among them. Similar to classical adsorption-based characterization methods, the deformation-based ones rely on the quality of the underlying theoretical framework. This fact stimulates the recent development of qualitative and quantitative models toward the more detailed description of a solid material, e.g. account of non-convex and corrugated pores, calculations of adsorption stress in realistic three-dimension solid structures, the extension of the existing models to new geometries, etc. The present review focuses on the theoretical description of adsorption-induced deformation in micro and mesoporous materials. We are aiming to cover recent theoretical works describing the deformation of both ordered and disordered porous bodies.more » « less
-
Ice nucleation in the atmosphere influences cloud properties, altering precipitation and the radiative balance, ultimately regulating Earth’s climate. An accepted ice nucleation pathway, known as deposition nucleation, assumes a direct transition of water from the vapor to the ice phase, without an intermediate liquid phase. However, studies have shown that nucleation occurs through a liquid phase in porous particles with narrow cracks or surface imperfections where the condensation of liquid below water saturation can occur, questioning the validity of deposition nucleation. We show that deposition nucleation cannot explain the strongly enhanced ice nucleation efficiency of porous compared with nonporous particles at temperatures below −40 °C and the absence of ice nucleation below water saturation at −35 °C. Using classical nucleation theory (CNT) and molecular dynamics simulations (MDS), we show that a network of closely spaced pores is necessary to overcome the barrier for macroscopic ice-crystal growth from narrow cylindrical pores. In the absence of pores, CNT predicts that the nucleation barrier is insurmountable, consistent with the absence of ice formation in MDS. Our results confirm that pore condensation and freezing (PCF), i.e., a mechanism of ice formation that proceeds via liquid water condensation in pores, is a dominant pathway for atmospheric ice nucleation below water saturation. We conclude that the ice nucleation activity of particles in the cirrus regime is determined by the porosity and wettability of pores. PCF represents a mechanism by which porous particles like dust could impact cloud radiative forcing and, thus, the climate via ice cloud formation.more » « less
