Engineering Papers⌕ Search

Engineering topics

Settgast, Randolph

Publications and source records attributed to Settgast, Randolph.

A multi-resolution approach to hydraulic fracture simulation

Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.

58 GEOSCIENCES↗

Interface and embedded joint methods for modeling dynamic fracture opening by explosive products [Two methods for modeling dynamic fracture opening by explosive products]

Two computational approaches are proposed in the paper to model dynamic fracture opening by explosive products. The first method assumes that the fractures may be modeled using flow elements embedded along the mesh lines. This method models crack opening in a straightforward way by splitting the nodes of the computational grid. It can account for crack branching; however, the crack directions are constrained by existing mesh faces, which may lead to mesh dependence. Also, the stress in flow elements is calculated explicitly separate from the surrounding solid elements that can impose additional limits on the time step stability condition for explicit integration. The second approach uses embedded flow elements to model the cracks. Typical thickness of the cracks is much smaller than the element size. Therefore, gas pressure in the cracks is assumed to be in stress equilibrium with the element stress. To achieve this, the crack thickness and the state of the gas is updated simultaneously with the state of the solid element which contains the crack. Therefore, the time step is controlled by the explicit solver applied for the solid and does not depend on the thickness of the crack. In conclusion, the main disadvantage of the second approach is due to the complexity of modeling multiple intersecting cracks, which go through the same element. We discuss the areas of possible applications of these 2 methods and the ways to improve and enhance them for future practical applications.

58 GEOSCIENCES↗

GEOS Code Development Road Map

GEOS is a massively parallel computational framework designed to enable HPC-based simulations of subsurface reservoir stimulation activities with the goal of optimizing current operations and evaluating innovative stimulation methods. GEOS will enable coupling of different solvers associated with the various physical processes occurring during reservoir stimulation in unique and sophisticated ways, adapted to various geologic settings, materials and stimulation methods. The overall architecture of the framework includes consistent data structures and will allow incorporation of additional physical and materials models as demanded by future applications. Along with predicting the initiation, propagation and reactivation of fractures, GEOS will also generate a seismic source term that can be linked with seismic wave propagation codes to generate synthetic microseismicity at surface and downhole arrays. Similarly, the output from GEOS can be linked with existing fluid/thermal transport codes. GEOS can also be linked with existing, non-intrusive uncertainty quantification schemes to constrain uncertainty in its predictions and sensitivity to the various parameters describing the reservoir and stimulation operations. We anticipate that an implicit-explicit 3D version of GEOS, including a preliminary seismic source model, will be available for parametric testing and validation against experimental and field data by Oct. 1, 2013.

02 PETROLEUM↗

PowderSim: Lagrangian Discrete and Mesh-Free Continuum Simulation Code for Cohesive Soils

PowderSim is a calculation tool that combines a discrete-element method (DEM) module, including calibrated interparticle-interaction relationships, with a mesh-free, continuum, SPH (smoothed-particle hydrodynamics) based module that utilizes enhanced, calibrated, constitutive models capable of mimicking both large deformations and the flow behavior of regolith simulants and lunar regolith under conditions anticipated during in situ resource utilization (ISRU) operations. The major innovation introduced in PowderSim is to use a mesh-free method (SPH-based) with a calibrated and slightly modified critical-state soil mechanics constitutive model to extend the ability of the simulation tool to also address full-scale engineering systems in the continuum sense. The PowderSim software maintains the ability to address particle-scale problems, like size segregation, in selected regions with a traditional DEM module, which has improved contact physics and electrostatic interaction models.

Johnson, Scott↗