An official website of the United States government
Abstract Dynamic shearing banding and fracturing in unsaturated porous media are significant problems in engineering and science. This article proposes a multiphase micro‐periporomechanics (PPM) paradigm for modeling dynamic shear banding and fracturing in unsaturated porous media. Periporomechanics (PPM) is a nonlocal reformulation of classical poromechanics to model continuous and discontinuous deformation/fracture and fluid flow in porous media through a single framework. In PPM, a multiphase porous material is postulated as a collection of a finite number of mixed material points. The length scale in PPM that dictates the nonlocal interaction between material points is a mathematical object that lacks a direct physical meaning. As a novelty, in the coupled PPM, a microstructure‐based material length scale is incorporated by considering micro‐rotations of the solid skeleton following the Cosserat continuum theory for solids. As a new contribution, we reformulate the second‐order work for detecting material instability and the energy‐based crack criterion and J‐integral for modeling fracturing in the PPM paradigm. The stabilized Cosserat PPM correspondence principle that mitigates the multiphase zero‐energy mode instability is augmented to include unsaturated fluid flow. We have numerically implemented the novel PPM paradigm through a dual‐way fractional‐step algorithm in time and a hybrid Lagrangian–Eulerian meshfree method in space. Numerical examples are presented to demonstrate the robustness and efficacy of the proposed PPM paradigm for modeling shear banding and fracturing in unsaturated porous media.
more »
« less
Energetic formulation of large‐deformation poroelasticity
Abstract The modeling of coupled fluid transport and deformation in a porous medium is essential to predict the various geomechanical process such as CO2 sequestration, hydraulic fracturing, and so on. Current applications of interest, for instance, that include fracturing or damage of the solid phase, require a nonlinear description of the large deformations that can occur. This paper presents a variational energy‐based continuum mechanics framework to model large‐deformation poroelasticity. The approach begins from the total free energy density that is additively composed of the free energy of the components. A variational procedure then provides the balance of momentum, fluid transport balance, and pressure relations. A numerical approach based on finite elements is applied to analyze the behavior of saturated and unsaturated porous media using a nonlinear constitutive model for the solid skeleton. Examples studied include the Terzaghi and Mandel problems; a gas–liquid phase‐changing fluid; multiple immiscible gases; and unsaturated systems where we model injection of fluid into soil. The proposed variational approach can potentially have advantages for numerical methods as well as for combining with data‐driven models in a Bayesian framework.
more »
« less
- Award ID(s):
- 2012259
- PAR ID:
- 10447139
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal for Numerical and Analytical Methods in Geomechanics
- Volume:
- 46
- Issue:
- 5
- ISSN:
- 0363-9061
- Page Range / eLocation ID:
- p. 910-932
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
The quasi-static problem describes a nonlinear porous body with a non-penetrating Barenblatt’s crack driven by the fracturing fluid, and its propagation is under investigation. By this, a bulk modulus of the porous body depends linearly on the density, the fracture faces allow contact with cohesion, and leak-off of the fluid into reservoir is accounted by the model. The mathematical problem consists in finding time-continuous functions of a displacement and a mean fluid pressure in the fracture, which satisfy the coupled system of the variational inequality and the fluid mass balance, which is controlled by the volume of fracking fluid pumped into the fracture. Well-posedness of the governing relations is proved rigorously by applying the method of Lagrange multipliers and using optimality conditions for the constrained minimization problem. As anillustrative example, a numerical benchmark problem of the fluid-driven fracture is presented in one dimension and computed by a Newton-type algorithm.more » « less
-
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
-
Abstract Many geo‐engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase‐field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase‐field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy‐based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase‐field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase‐field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. The comparison of numerical results for two‐dimensional test cases with existing analytical solutions demonstrates that the proposed phase‐field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three‐dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.more » « less
