Hybrid mimetic finite-difference and virtual element formulation for coupled poromechanics

We present a hybrid mimetic finite-difference and virtual element formulation for coupled single-phase poromechanics on unstructured meshes. The key advantage of the scheme is that it is convergent on complex meshes containing highly distorted cells with arbitrary shapes. We use a local pressure-jump stabilization method based on unstructured macro-elements to prevent the development of spurious pressure modes in incompressible problems approaching undrained conditions. A scalable linear solution strategy is obtained using a block-triangular preconditioner designed specifically for the saddle-point systems arising from the proposed discretization. The accuracy and efficiency of our approach are demonstrated numerically on two-dimensional benchmark problems.Comment: 25 pages, 17 figure

A phase-field model for hydraulic fracture nucleation and propagation in porous media

Many geo-engineering applications, e.g., 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, e.g., 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 (2D) 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 (3D) domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier

Dynamic Soil-Foundation-Structure Interaction Analyses of Large Caissons

Large cellular reinforced concrete caissons exist as foundations of major long-span bridges across waterways in many parts of the country. This study was conducted to evaluate the important factors affecting the seismic response of large caissons. The paper presents the results of equivalent linear and non-linear analyses performed for a typical caisson idealized based on the cellular caisson at Pier W3 of the West San Francisco Bay Bridge subject to ground motion with a peak rock acceleration of 0.6 g. This caisson is 38.7 m (127 fi) long by 22.9 m (75 ft) wide submerged in about 32.6 m (107 ft) of water. It is embedded in 33.5 m (110 fi) of soil deposits and is founded on rock. Equivalent linear 3-D and 2-D analyses conducted in the direction of the short axis (longitudinal) were performed using a modified version of computer program SASSI. The results of these 3-D and 2-D analyses are similar. Non-linear analyses were performed for 2-D models using computer program FLAC. The results indicate that side gapping, base lifting, interface sliding, and soil yielding reduce the earth pressure, base bearing stress, caisson shear and bending moment, and caisson motions. However, the frequency characteristics of the responses appear to be relatively unaffected