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
-
Abstract Krauklis waves are generated by pressure disturbances in fluid‐filled cavities and travel along the solid‐fluid interface. Their far‐field radiation, observed in seismic data from volcanoes or hydraulic fracturing, is known as long‐period events. Characterized by low velocity and resonance, Krauklis waves help estimate fracture size and discern fluids in saturated fractures. Despite numerous theoretical models analyzing Krauklis waves, the existing paradigms are founded on static flow conditions. However, in geological contexts, the assumption of static flow may not be valid. We developed an experimental apparatus using a tri‐layer model consisting of a pair of aluminum plates to examine the effect of fluid flow on Krauklis waves. We employed an infusion syringe pump to inject fluids into the fracture under different flow rates. We used water, oil, and an aqueous solution of Polyethylene glycol as fracture fluids. We calculated resonant frequency, phase velocity, and quality factor to characterize the Krauklis waves. Our findings reveal that an increase in flow rate leads to a higher phase velocity, higher quality factor, and a shift to higher resonant frequency when the flow is in the direction of initial wave propagation while decreasing amplitude. Additionally, when the flow is in the opposite direction of initial wave propagation, we note higher wave absorption and distortion of the Krauklis waves. Our observations unequivocally affirm that fluid flow leaves strong signatures on the Krauklis waves, providing a robust basis for characterizing fluid dynamics within geological settings through the analysis of Krauklis wave.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
