An official website of the United States government
This paper is concerned with the numerical solution of compressible fluid flow in a fractured porous medium. The fracture represents a fast pathway (i.e., with high permeability) and is modeled as a hypersurface embedded in the porous medium. We aim to develop fast-convergent and accurate global-in-time domain decomposition (DD) methods for such a reduced fracture model, in which smaller time step sizes in the fracture can be coupled with larger time step sizes in the subdomains. Using the pressure continuity equation and the tangential PDEs in the fracture-interface as transmission conditions, three different DD formulations are derived; each method leads to a space-time interface problem which is solved iteratively and globally in time. Efficient preconditioners are designed to accelerate the convergence of the iterative methods while preserving the accuracy in time with nonconforming grids. Numerical results for two-dimensional problems with non-immersed and partially immersed fractures are presented to show the improved performance of the proposed methods.
more »
« less
This content will become publicly available on July 15, 2026
Nonlinear model of a porous body with fluid-driven fracture under cohesion contact conditions and fluid volume control
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
- Award ID(s):
- 2307563
- PAR ID:
- 10644153
- Publisher / Repository:
- World Scientific
- Date Published:
- Journal Name:
- Mathematical Models and Methods in Applied Sciences
- Volume:
- 35
- Issue:
- 11
- ISSN:
- 0218-2025
- Page Range / eLocation ID:
- 2311 to 2328
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT This work proposes a novel approach for coupling non‐isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid‐filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface‐capturing phase‐field fracture method and interface‐tracking thermo‐fluid‐structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase‐field approach to represent fractures and reconstruct the geometry to frame a thermo‐fluid‐structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase‐field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid‐filled fracture with the phase‐field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase‐field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.more » « less
-
Traumatic brain injury remains a significant global health concern, requiring advanced understanding and mitigation strategies. In current brain concussion research, there is a significant knowledge gap: the critical role of transient cerebrospinal fluid (CSF) flow in the porous subarachnoid space (SAS) has long been overlooked. To address this limitation, we are developing a simplified mathematical model to investigate the CSF pressurization in the porous arachnoid trabeculae and the resulting motion of brain matter when the head is exposed to a translational impact. The model simplifies the head into an inner solid object (brain) and an outer rigid shell (skull) with a thin, porous fluid gap (SAS). The CSF flow in the impact side (coup region) and the opposite side (contrecoup region) is modeled as porous squeezing and expanding flows, respectively. The flow through the side regions, which connect these regions, is governed by Darcy's law. We found that the porous arachnoid trabeculae network significantly dampens brain motion and reduces pressure variations within the SAS compared to a SAS without the porous arachnoid trabeculae (AT). This effect is particularly pronounced under high-frequency, periodic acceleration impacts, thereby lowering the risk of injury. The dampening effect can be attributed to the low permeability of the AT, which increases resistance to fluid movement and stabilizes the fluid and pressure responses within the SAS, thereby reducing extreme pressure fluctuations and brain displacement under impact. This work provides a foundational understanding of CSF flow dynamics, paving the way for innovative approaches to brain injury prevention and management.more » « less
-
Summary We present an algorithm to simulate curvilinear hydraulic fractures in plane strain and axisymmetry. We restrict our attention to sharp fractures propagating in an isotropic, linear elastic medium and driven by the injection of a laminar, Newtonian fluid governed by lubrication theory, and we require the existence of a finite lag region between the fluid front and the crack tip. The key novelty of our approach is in how we discretize the evolving crack and fluid domains: we utilize universal meshes (UMs), a technique to create conforming triangulations of a problem domain by only perturbing nodes of a universal background mesh in the vicinity of the boundary. In this way, we construct meshes, which conform to the crack and to the fluid front. This allows us to build standard piecewise linear finite element spaces and to monolithically solve the quasistatic hydraulic fracture problem for the displacement field in the rock and the pressure in the fluid. We demonstrate the performance of our algorithms through three examples: a convergence study in plane strain, a comparison with experiments in axisymmetry, and a novel case of a fracture in a narrow pay zone.more » « less
-
A ubiquitous arrangement in nature is a free-flowing fluid coupled to a porous medium, for example a river or lake lying above a porous bed. Depending on the environmental conditions, thermal convection can occur and may be confined to the clear fluid region, forming shallow convection cells, or it can penetrate into the porous medium, forming deep cells. Here, we combine three complementary approaches – linear stability analysis, fully nonlinear numerical simulations and a coarse-grained model – to determine the circumstances that lead to each configuration. The coarse-grained model yields an explicit formula for the transition between deep and shallow convection in the physically relevant limit of small Darcy number. Near the onset of convection, all three of the approaches agree, validating the predictive capability of the explicit formula. The numerical simulations extend these results into the strongly nonlinear regime, revealing novel hybrid configurations in which the flow exhibits a dynamic shift from shallow to deep convection. This hybrid shallow-to-deep convection begins with small, random initial data, progresses through a metastable shallow state and arrives at the preferred steady state of deep convection. We construct a phase diagram that incorporates information from all three approaches and depicts the regions in parameter space that give rise to each convective state.more » « less
