Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

DEMO-FTES: Development, Monitoring, and Control of Fracture Thermal Energy Storage in Crystalline Rock Formations (CRADA Final Report)

The DEMO-FTES project investigated the feasibility of Fracture Thermal Energy Storage (FTES) as a seasonal energy storage solution in crystalline rock formations. FTES leverages hydraulically induced fractures to exchange heat between circulating fluids and the surrounding rock mass, enabling long-term thermal energy retention due to the high specific heat and low thermal conductivity of rock. This approach has the potential to reduce heating and cooling energy demands and enhance building energy resilience. The project combined dimensional analysis, numerical modeling, laboratory experiments, and meso-scale field tests to evaluate FTES performance and advance its technology readiness level from 3 to 5. Scaling analysis identified key dimensionless parameters governing heat transfer and fluid flow, ensuring laboratory and field tests were representative of larger-scale systems. Numerical simulations using TOUGH and iTOUGH2 frameworks supported experiment design and interpretation, modeling fracture geometry, thermal-hydraulic behavior, and thermo-mechanical coupling. Laboratory tests at EPFL involved creating single and multiple fractures in 25 cm cubic samples of Gabbro and Granite under true triaxial stress.

25 ENERGY STORAGE↗

Phase Field Modeling of Near-Wellbore Hydraulic Fracture Nucleation and Propagation

The goal of Task 2 of the project titled "Closing the loop between in situ stress complexity and EGS fracture complexity" is to model fracture nucleation and propagation in the near wellbore region and, in combination with the experimental study conducted as part of Task 4, investigate the relationship between the in situ stress and the near-wellbore fracture patterns. To this end, We have identified phase field as a promising modeling approach to model the nucleation and propagation process in the near wellbore region. In phase field, fractures are represented by a diffuse variable, i.e the damage. Compared to approaches that model fractures as sharp interfaces, phase field does not require any re-meshing or element insertion algorithm, which makes the treatment of complex fracture geometries straightforward, even with simple meshes.

15 GEOTHERMAL ENERGY↗

Utah FORGE 5-2565: Hydrothermal Evolution of Fracture Properties - Workshop Presentation

This is a presentation on the Evolution of Permeability and Strength Recovery of Shear Fractures Under Hydrothermal Conditions project by the U.S. Geological Survey, presented by Dr. David Lockner. The project's objective was to determine how thermal, hydraulic, mechanical, and chemical processes affect the sustainability of fracture networks in geothermal reservoirs and provide strategies for improved EGS techniques that maximize thermal coupling and increase reservoir longevity. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 8, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗

Investigation of fracture in porous materials: a phase-field fracture study informed by ReaxFF

Microscopic features (e.g, pore shapes, sizes, and distribution) in porous material substantially affect the overall mechanical properties such as stiffness and strength. In turn, these material properties determine the macroscopic behaviors of fracture in the porous material. In certain cases, macroscopic properties can be derived from the porous skeleton and void ratio (i.e., porosity), but in many other cases, such derivation is a challenging task. This paper presents a numerical investigation of microporosity and micropore shapes effect on the macrofracture behavior in porous amorphous silica. For this study, we extend the recently-proposed combined molecular dynamic (MD) and phase-field (PF) fracture modeling approach by including different pore shapes in the atomistic domain. In the MD simulations, we adopt ReaxFF to evaluate the material properties, where four different micropore cases are considered. Based on the material properties derived from MD simulations, the macrofracture propagation of porous media is studied using hybrid PF simulation. In the characterization of the pore structure, the concept of pore ligament is proposed to relate the pore shape and the critical energy release rate. Two classical fracture problems were used to evaluate the effect of pore shape on the macrofracture behavior. The results of the case study show that although the micropore shapes change the macrofracture behaviors, these effects vary with the geometry and loading conditions of macroscopic boundary value problems. Furthermore, the case study also shows that the influence of micropore structure can be captured at the macroscopic level through the material properties derived from the MD simulations.

42 ENGINEERING↗

Characterization of an electrically conductive proppant for fracture diagnostics

Fracture diagnosis with electromagnetic (EM) and electrical tools requires proppants with high electrical conductivity and mechanical strength. Lab measurements of the electrical and hydraulic conductivity of proppants are critical for selecting the best candidates. Such measurements greatly benefit simulations, field tests, and the ultimate application of such proppants in the field. To that end, a new lab protocol is developed for measuring the electrical and hydraulic conductivity of proppants. The lab setup, which mainly includes a resistivity core holder and a Hassler sleeve core holder, allows for simulation of realistic pressure and temperature conditions when making measurements. Petroleum coke (PC) is proposed as a candidate proppant because of its widespread availability and low cost. Lab measurements show that the effective electrical conductivity of pure PC in a model fracture is approximately 5000 S/m, under a closure stress greater than [Formula: see text] (4000 psi). When PC is mixed with sand, the effective electrical conductivity of the mixture decreases with an increasing weight percentage of sand. Although sand degrades the contact between PC particles, the electrical conductivity stays reasonably high (approximately 1700 S/m) when 50% sand is added. Hydraulic conductivity measurements show that when a fracture is propped with pure PC, the measured fracture conductivity is greater than [Formula: see text] ([Formula: see text]) (dimensionless fracture conductivity greater than 100 for a shale with [Formula: see text] or 100 nD permeability) under a confining pressure of [Formula: see text] (6000 psi). This means that a fracture propped with PC is infinitely conductive in a typical shale formation. When sand is added, the fracture’s hydraulic conductivity becomes even higher, which clearly shows PC’s ability of sustaining high stresses. The proposed protocol provides a robust and effective method that can be generalized for lab testing for other candidate proppants. The data presented clearly show that PC has the potential for field-scale applications in EM hydraulic fracture diagnostics.

Geochemistry & Geophysics↗

The mechanics of brittle granular materials with coevolving grain size and shape

The influence of particle shape on the mechanics of sand is widely recognized, especially in mineral processing and geomechanics. However, most existing continuum theories for engineering applications do not encompass the morphology of the grains and its evolution during comminution. Similarly, the relatively few engineering models accounting for grain-scale processes tend to idealize particles as spheres, with their diameters considered as the primary and sole geometric descriptor. This paper inspires a new generation of constitutive laws for crushable granular continua with arbitrary, yet evolving, particle morphology. We explore the idea of introducing multiple grain shape descriptors into Continuum Breakage Mechanics (CBM), a theory originally designed to track changes in particle size distributions during confined comminution. We incorporate the influence of these descriptors on the elastic strain energy potential and treat them as dissipative state variables. In analogy with the original CBM, and in light of evidence from extreme fragmentation in nature, the evolution of the additional shape descriptors is postulated to converge towards an attractor. Comparisons with laboratory experiments, discrete element analyses and particle-scale fracture models illustrate the encouraging performance of the theory. The theory provides insights into the feedback among particle shape, compressive yielding and inelastic deformation in crushable granular continua. These results inspire new questions that should guide future research into crushable granular systems using particle-scale imaging and computations.

36 MATERIALS SCIENCE↗

Crack tip microplasticity mediated by microstructure gradients

Abstract Traditional fracture theories infer damage at cracks (local field) by surveying loading conditions away from cracks (far field). This approach has been successful in predicting ductile fracture, but it normally assumes isotropic and homogeneous materials. However, myriads of manufacturing procedures induce heterogeneous microstructural gradients that can affect the accuracy of traditional fracture models. This work presents a microstructure‐sensitive finite element approach to explore the shielding effects of grain size and crystallographic orientation gradients on crack tip microplasticity and blunting. A dislocation density‐based crystal plasticity model conveys texture evolution, grain size effects, and directional hardening by computing the constraint from dislocation structures. The results demonstrate that the microstructure can act as a buffer between the local and far fields that affects the crack tip microplasticity variability. For nominal opening loading, grain size and texture affect the local ductility and induce a non‐negligible multiaxial plastic deformation. Furthermore, driving forces based on measuring displacements away from the crack tip are less affected by the microstructure, which suggests that traditional experimental methods smear out important crack tip variability.

Castelluccio, Gustavo M.↗

Meso-scale modelling of compressive fracture in concrete with irregularly shaped aggregates

This paper presents a meso-scale modelling framework to investigate the fracture process in concrete subjected to uniaxial and biaxial compression accounting for its mesostructural characteristics. 3D mesostructure of concrete consisting of coarse aggregates, mortar and interfacial transition zone between them was developed using an in-house code based on the Voronoi tessellation and splining method, which enables to generate the realistic-look aggregates with controllable structural features such as content, location, size and shape. Based on the generated 3D mesostructure, the concrete damage plasticity approach was employed to simulate the compressive fracture behaviour of concrete in terms of crack morphology and stress-strain response against the shape parameters of aggregate. Results indicate that the shape of aggregate has a negligible effect on compressive strength of concrete, which is highly associated with the random location and size distribution of aggregate. The aggregate irregularity has a significant influence on crack initiation and growth of concrete.

36 MATERIALS SCIENCE↗

Implementation of Creep-Damage Model for Concrete Fracture in MOOSE

Numerical implementation of an isotropic creep-damage model for concrete in multiphysics object-oriented simulation environment (MOOSE) finite element framework is presented in this paper. The constitutive model considers the combined effect of instantaneous and delayed strains on damage propagation. The implementation considers creep using generalized Maxwell or Kelvin-Voigt models. Using strain splitting assumptions, the total mechanical strains are split into elastic and creep components. Damage is considered to evolve as a function of the elastic and creep strains. This work considers damage as a function of fracture energy using the characteristic length of each finite element. This approach preserves the energy release rate of each element and avoids vanishing energy dissipation as the mesh is refined. Further, a creep-damage parameter is used to quantify the effect of creep strain on damage. The model is tested against published results on notched three-point bending specimens involving non-linear creep and predicts that about a third of the creep strain contributes towards damage for the experiments simulated. Results show that the proposed framework has predictive capabilities, and the model can be extended for more complex systems.

36 MATERIALS SCIENCE↗

Using Homogenized Models to Explore the Effect of Fracture Densities on Weathering

Despite its importance, only a few researchers have incorporated the effects of fracturing into models of reactive transport for rock weathering. Here we explore 2D simulations that describe weathering under conditions of diffusive and advective transport within heterogeneous media consisting of rocky blocks and fractures. In our simulations, the Darcy velocities vary in space and time and depend on weathering processes within the rock matrix. We explore simulations with saturated and unsaturated flow for weathering bedrock that consists of blocks separated by inert or weathered material. The simulations show that a simplified homogenized model can approximate exact solutions for some of the simulated columns and hills and can allow exploration of coupling between flow and reaction in fractured rock. These hillslope simulations document that, even in the presence of 2D water flow, i) an increase in fracture density results in faster weathering advance rates; and ii) the water table locates deeper for a rock system that is weathered and fractured rather than weathered and unfractured. Some of these patterns have also been observed for natural systems. But these simulations also highlight how simplified models that do not use appropriate averaging of heterogeneities can be inaccurate in predicting weathering rate for natural systems. For example, if water flows both vertically and laterally through the vadose zone of a hill, then a prediction of the depth of regolith that is based on modeling strictly unidirectional downward infiltration will be unrealistically large. Likewise, if the fracture density observed near the land surface is used in a model to predict depth of weathering for a system where the fracture density decreases downward, the model will overestimate regolith depth. Learning how to develop accurately homogenized models could thus enable better conceptual models and predictions of weathering advance in natural systems.

58 GEOSCIENCES↗

Study of Foam Fracturing for Enhanced Geothermal Systems Using Model Material

Foam fracturing is considered a potential approach to address water concerns with hydraulic fracturing in the development of enhanced geothermal systems (EGS). In many EGS sites, water required for hydraulic fracturing is either unavailable, extremely costly or environmentally unsustainable. This paper presents work performed in the first year of a project sponsored by the U.S. DOE GTO Waterless Stimulation Initiative that investigates foam fracturing as an alternative approach to hydraulic fracturing. One element of this work explores foam formulations that are suitable for EGS conditions. A separate paper submitted to this workshop covers the relevant work on foam characterization (Thakore et al., 2020).This paper describes the development of an experimental setup for studying the applicability of foam fluids for hydraulic fracturing along with enhanced fracture efficiency through cyclic pressurization. The system is equipped with pulse rate automation, dual pulse valves for enhanced pulsation amplitudes, and other functions. It can perform static injection at a given pressurization rate, and pulsed injection at specified pulse size and rate. Both single-phase and foamed fluids can be used. Initial experimental results of foam fracturing using cement as a model material are reported. The observations from the experimental work will be presented and discussed.

Wang, Hong↗

Finite Element Modeling of Diffusion in Fractured Porous Media by Using Hierarchical Material Properties

Abstract Fractured porous media challenge modeling approaches due to high computational costs and excessive mesh refinement imposed by the extreme scale variability of fractures and the heterogeneity of the surrounding porous rock. To overcome such difficulties, we utilize the hierarchical finite element method ( Hi ‐FEM) that has been developed previously to simulate the electrical potential distribution in complex geologic environments. The method employs the hierarchical basis functions in classical finite element analysis to enable representation of material properties on each dimensional component of a given 3D unstructured finite element, thereby inherently allowing for interactions at the boundary between fracture and a host rock. In this study, we extend its application to transient fluid flow and heat conduction in the Laplace domain. Time‐domain flow solutions are obtained by numerical inverse Laplace transform. We evaluate the accuracy of the method using different flow models and demonstrate its robustness for large‐scale, rock mass models featuring complex fracture networks. Moreover, for the computation of nodal Darcian velocity fields in fractured porous media where the fractures are represented as 2D features, a new approach that employs the Yeh's Galerkin model for both volume and facet elements is proposed. Results show that Hi ‐FEM can produce accurate flow solutions for fractured porous media without any need of coupling or transfer mechanism while still being computationally economical and numerically robust, even for large‐scale simulations.

Beskardes, G. D.↗

Application of an Embedded Fracture and Borehole Modeling Approach to the Understanding of EGS Collab Experiment 1

EGS Collab is a series of meso-scale experiments and associated numerical simulation activities being funded by the United States Department of Energy, Geothermal Technologies Office (GTO) to investigate enhanced geothermal system processes under in-situ stress and slightly elevated temperature conditions. This project is designed to provide scientists and engineers with immediate access to impermeable rock at scales larger than possible in the laboratory, but generally smaller than those for commercial production. Immediate access to rock is provided via the existing drifts of the former Homestake Gold Mine, now operated as the Sanford Underground Research Facility in Lead, South Dakota. The objectives of the EGS Collab project are to develop well controlled fracture networks between injector and producer boreholes using normal and shear stimulation for permeability enhancement. The first experimental site was located off the West Access drift on the 4850 Level (4850 feet below ground surface) in phyllite of the Precambrian Poorman formation, and involved the creation of a fracture network comprising a combination of hydraulic and natural fractures. A second experimental site is now being considered near the battery alcove on the 4100 Level in amphibolite and rhyolite of the Yates Unit. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. Experiment 1 has comprised a series of successful tests, including a long-term chilled-water circulation test, but the testbed has two atypical EGS elements. Active ventilation in the adjacent drift over a 50-year period cooled the testbed rock mass and hydraulic fracturing intersected monitoring boreholes making them conduits for fluid flow. Numerical simulations executed in support of the design of the EGS Collab Experiment 1, computed a radial temperature distribution and stress gradient orthogonal to the drift axis, resulting in the forecast of an oblong hydraulic fracture geometry, extended in the direction of the drift from the stimulation borehole. This paper describes the numerical simulation of long-term chilled-water circulation test with an embedded fracture and borehole modeling approach. The principal objective of the simulation work is show agreement between experimental observations in terms of production fluid temperatures, monitoring borehole temperatures, injection and production pressures, and tracer recoveries using a single conceptual model for the fracture network, using characterization and monitoring data generally available for EGS. A secondary objective is to improve the numerical simulation result comparisons with additional data available from the broader monitoring equipment within the testbed versus deeper, hotter, and more remote EGS. One complicating factor for comparisons between the numerical simulation results and experimental observations, is the Joule-Thomson heating associated with large pressure drops across the fracture network.

numerical simulation, embedded fractures and boreh↗

Hydro-mechanical coupled modeling of hydraulic fracturing using the mesh fragmentation technique

The inclusion of special interface elements between standard finite elements is becoming a common technique to model discontinuities in porous media using the finite element method (FEM). Recently, it was proved that high aspect ratio (HAR) elements can reproduce very satisfactorily the main features of behavior observed in a single hydraulically-induced fracture. In this paper, the mesh fragmentation technique (MFT) is proposed to extend the application of HAR elements to simulate more general problems involving the formation of hydraulic fractures in rocks by introducing this type of element in-between the standard FE of a typical mesh. The proposed framework makes use of standard FE techniques. Therefore, it is relatively easy to upgrade an existing FE program for continuous porous media to deal with evolving discontinuities. The proposed approach is verified and validated against available analytical and numerical solutions. The crack patterns predicted by the MFT are also compared with the ones observed in the laboratory. The formation of multiple fractures is also analyzed. In all cases, the proposed technique was able to properly capture the main features of rock behavior subjected to hydraulic fracturing.

58 GEOSCIENCES↗

Atomistic measurement and modeling of intrinsic fracture toughness of two-dimensional materials

Quantifying the intrinsic mechanical properties of two-dimensional (2D) materials is essential to predict the long-term reliability of materials and systems in emerging applications ranging from energy to health to next-generation sensors and electronics. Currently, measurements of fracture toughness and identification of associated atomistic mechanisms remain challenging. Herein, we report an integrated experimental–computational framework in which in-situ high-resolution transmission electron microscopy (HRTEM) measurements of the intrinsic fracture energy of monolayer MoS 2 and MoSe 2 are in good agreement with atomistic model predictions based on an accurately parameterized interatomic potential. Changes in crystalline structures at the crack tip and crack edges, as observed in in-situ HRTEM crack extension tests, are properly predicted. Such a good agreement is the result of including large deformation pathways and phase transitions in the parameterization of the inter-atomic potential. The established framework emerges as a robust approach to determine the predictive capabilities of molecular dynamics models employed in the screening of 2D materials, in the spirit of the materials genome initiative. Moreover, it enables device-level predictions with superior accuracy (e.g., fatigue lifetime predictions of electro- and opto-electronic nanodevices).

Zhang, Xu↗