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
A Darcy‐Brinkman‐Biot Approach to Modeling the Hydrology and Mechanics of Porous Media Containing Macropores and Deformable Microporous Regions
Abstract The coupled hydrology and mechanics of soft porous materials (such as clays, hydrogels, membranes, and biofilms) is an important research area in several fields, including water and energy technologies as well as biomedical engineering. Well‐established models based on poromechanics theory exist for describing these coupled properties, but these models are not adapted to describe systems with more than one characteristic length scale, that is, systems that contain both macropores and micropores. In this paper, we expand upon the well‐known Darcy‐Brinkman formulation of fluid flow in two‐scale porous media to develop a “Darcy‐Brinkman‐Biot” formulation: a general coupled system of equations that approximates the Navier‐Stokes equations in fluid‐filled macropores and resembles the equations for poroelasticity in microporous regions. We parameterized and validated our model for systems that contain either plastic (swelling clay) or elastic microporous regions. In particular, we used our model to predict the permeability of an idealized siliciclastic sedimentary rock as a function of pore water salinity and clay content. Predicted permeability values are well described by a single parametric relation between permeability and clay volume fraction that agrees with existing experimental data sets. Our novel formulation captures the coupled hydro‐chemo‐mechanical properties of sedimentary rocks and other deformable porous media in a manner that can be readily implemented within the framework of Digital Rock Physics.
more »
« less
- Award ID(s):
- 1752982
- PAR ID:
- 10443955
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 55
- Issue:
- 10
- ISSN:
- 0043-1397
- Page Range / eLocation ID:
- p. 8096-8121
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In this article, we consider a phase field model with different densities and viscosities for the coupled two-phase porous media flow and two-phase free flow, as well as the corresponding numerical simulation. This model consists of three parts: a Cahn–Hilliard–Darcy system with different densities/viscosities describing the porous media flow in matrix, a Cahn–Hilliard–Navier–Stokes system with different densities/viscosities describing the free fluid in conduit, and seven interface conditions coupling the flows in the matrix and the conduit. Based on the separate Cahn–Hilliard equations in the porous media region and the free flow region, a weak formulation is proposed to incorporate the two-phase systems of the two regions and the seven interface conditions between them, and the corresponding energy law is proved for the model. A fully decoupled numerical scheme, including the novel decoupling of the Cahn–Hilliard equations through the four phase interface conditions, is developed to solve this coupled nonlinear phase field model. An energy-law preservation is analyzed for the temporal semi-discretization scheme. Furthermore, a fully discretized Galerkin finite element method is proposed. Six numerical examples are provided to demonstrate the accuracy, discrete energy law, and applicability of the proposed fully decoupled scheme.more » « less
-
We consider the interaction between a free flowing fluid and a porous medium flow, where the free flowing fluid is described using the time dependent Stokes equations, and the porous medium flow is described using Darcy’s law in the primal formulation. To solve this problem numerically, we use a diffuse interface approach, where the weak form of the coupled problem is written on an extended domain which contains both Stokes and Darcy regions. This is achieved using a phase-field function which equals one in the Stokes region and zero in the Darcy region, and smoothly transitions between these two values on a diffuse region of width (ϵ) around the Stokes-Darcy interface. We prove convergence of the diffuse interface formulation to the standard, sharp interface formulation, and derive rates of convergence. This is performed by deriving a priori error estimates for discretizations of the diffuse interface method, and by analyzing the modeling error of the diffuse interface approach at the continuous level. The convergence rates are also shown computationally in a numerical example.more » « less
-
抄録 It is important to understand the long-term migration of radionuclides when considering long-lasting rock engineering projects such as the geological disposal of radioactive waste. The network of fractures and pores in a rock mass plays a major role in fluid migration as it provides pathways for fluid flow. The geometry of such a network can change due to fracture sealing by fine-grained material over extended periods of time. Groundwater commonly contains fine-grained material such as clay minerals, and it is probable that such minerals accumulate within rock fractures during groundwater flow, thereby decreasing fracture apertures and bulk permeability. It is therefore essential to conduct permeability measurements using water that includes fine-grained minerals in order to understand the evolving permeability characteristics of rock. However, this has not been studied to date in in-situ rock mass. Therefore, in the present study, we perform permeability measurements in a granite rock mass to investigate the change of permeability that occurs under the flow of water that includes clays. Our findings show that clay particles accumulate in fractures and that the permeability (hydraulic conductivity) of the granite rock mass decreases over time. The decrease was more significant in the earlier time. We conclude that the accumulation of clay minerals in the fracture decreases the permeability of a rock mass. Furthermore, we consider that the filling and closure of fractures in rock is possible under the flow of groundwater that contains clay minerals.more » « less
-
We develop a mixed finite element method for the coupled problem arising in the interaction between a free fluid governed by the Stokes equations and flow in deformable porous medium modeled by the Biot system of poroelasticity. Mass conservation, balance of stress, and the Beavers–Joseph–Saffman condition are imposed on the interface. We consider a fully mixed Biot formulation based on a weakly symmetric stress-displacement-rotation elasticity system and Darcy velocity-pressure flow formulation. A velocity-pressure formulation is used for the Stokes equations. The interface conditions are incorporated through the introduction of the traces of the structure velocity and the Darcy pressure as Lagrange multipliers. Existence and uniqueness of a solution are established for the continuous weak formulation. Stability and error estimates are derived for the semi-discrete continuous-in-time mixed finite element approximation. Numerical experiments are presented to verify the theoretical results and illustrate the robustness of the method with respect to the physical parameters.more » « less
