Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fault friction”

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 73 records · Page 4

Seismological Stress Drops for Confined Ruptures Are Invariant to Normal Stress

Seismic moment and rupture length can be combined to infer stress drop, a key parameter for assessing earthquakes. In natural earthquakes, stress drops are largely depth-independent, which is surprising given the expected dependence of frictional stress on normal stresses and hence overburden. We have developed a transparent experimental fault that allows direct observation of thousands of slip events, with ruptures that are fully contained within the fault. Surprisingly, the observed stress drops are largely independent of both the magnitude of normal stress and its heterogeneity, capturing the independence seen in nature. However, we observe larger, normal stress-dependent stress drops when the fault area is reduced, which allows slip events to frequently reach the edge of the interface. We conclude that confined ruptures have normal stress independent stress drops, and thus the depth-independent stress drops of tectonic earthquakes may be a consequence of their confined nature.

58 GEOSCIENCES↗

A workflow to assess the efficacy of brine extraction for managing injection-induced seismicity potential using data from a CO 2 injection site near Decatur, Illinois

Injection of CO 2 for storage in a geologic formation increases pore pressure and alters in situ stresses. Depending on the orientation of any existing fault and fracture planes, such as critically stressed planes, this stress alteration will modify normal stresses acting on planes and could result in frictional sliding and release stored energy in the form of seismicity. Brine extraction (BE) is a technique that can be applied prior to, or during, CO 2 injection to reduce pore pressure for increasing storage capacity and, potentially, for reducing the likelihood of frictional sliding. Here a workflow is described to assess the efficacy of BE for mitigating frictional sliding (i.e., seismicity) during injection and entails: site characterization, stress calculations and failure assessment, static and dynamic modeling, and BE operational planning. Site characterization describes the stress field used to calculate the Coulomb Failure Function (CFF) that constrains allowable pore pressure changes and injection rates in the numerical simulation of CO 2 injection scenarios. The inclusion of BE in the workflow allows for determination of the potential need for pressure reduction, and evaluation of the effectiveness of this operation. Example application of the workflow using an injection field dataset near Decatur, IL, provides insight on fracture planes and stresses at the site, formation properties and the impact of variable CO 2 injection-rate targets on whether BE plans are required. The study workflow indicates that BE could enhance CO 2 injection rate by 39% and correspondingly reduce the potential for injection-induced seismicity as indicated by a reduction in CFF.

58 GEOSCIENCES↗

Degradation of Phonons in Disordered Moiré Superlattices

The elastic collective modes of a moiré superlattice arise not from vibrations of a rigid crystal but from the relative displacement between the constituent layers. Despite their similarity to acoustic phonons, these modes, called phasons, are not protected by any conservation law. Furthermore, we show that disorder in the relative orientation between the layers and thermal fluctuations associated with their sliding motion degrade the propagation of sound in the moiré superlattice. Specifically, the phason modes become overdamped at low energies and acquire a finite gap, which displays a universal dependence on the twist-angle variance. Thus, twist-angle inhomogeneity is manifested not only in the non-interacting electronic structure of moiré systems, but also in their phonon-like modes. More broadly, our results have important implications for the electronic properties of twisted moiré systems that are sensitive to the electron-phonon coupling.

2-dimensional systems↗

Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties

This project developed machine learning (ML) methods, lab data sets, and field data to advance geothermal exploration and geothermal energy production. The work had three focus areas. One involved the development of ML methods to use microearthquakes (MEQs) for imaging geothermal reservoir properties and improving subsurface characterization – most importantly the evolution of permeability within the evolving reservoir. This part of the work included development of ML approaches for automated MEQ location, focal mechanism determination and identification of earthquake precursors. The second area focused on using MEQ signals generated by geothermal exploration and production to predict the relationship between fluid injection and seismicity. Here, we extended to reservoir scale our success in using ML to predict laboratory earthquakes and fault zone stress state. The third focus area was on lab experiments. Here, we developed new ML models for lab earthquake prediction and identification of precursors to failure to improve earthquake forecasting and early warning in geothermal settings. Major outcomes of our work include ML models that learn from MEQ signals during geothermal exploration and production to predict induced seismicity. MEQs occur naturally in connection with drilling and energy production. We developed ML methods to use the seismic waves from these events to characterize the elastic, hydraulic and poromechanical properties of reservoirs. Our work illuminated fracture geometry and the evolution of fracture permeability by incorporating seismic coda wave analysis and ML methods to relate fluid injection and seismicity. We significantly expanded laboratory earthquake prediction to include methods that use both passive measurements of microearthquakes within the lab fault zones and also active source acoustic measurements of fault zone elastic properties. These methods can now predict fault zone stress state, time to failure and the magnitude of lab earthquakes. Our work showed that repetitive stick- slip failure events during frictional sliding (the lab equivalent of earthquakes) are preceded by a cascade of micro-failure events that radiate energy in a manner that foretells unstable failure – manifest as laboratory MEQs. We documented a mapping between fracture properties and statistical attributes of elastic radiation. We extended existing works to geothermal reservoir scale and developed ML methods to determine reservoir permeability, fracture properties, and their evolution during geothermal energy production. An attractive feature of ML algorithms is their ability to handle big datasets and reveal patterns and correlations that may remain invisible to conventional analyses. Our work connected data from field, laboratory and intermediate scales to study permeability, stress, strength, fracture stiffness and geometry. At the field scale we used data from the Newberry Volcano field site, UtahFORGE, EGS Collab, and also the Bedretto underground research lab in Switzerland. These data sets are bridging the gap between the lab scale, theory, and reservoir scale. Our work produced plain language summaries to improve public understanding of DOE research. We also developed openly distributed ML and seismicity datasets for use by all researchers and we published connections between induced seismicity in geothermal areas and reservoir properties including permeability, fracture properties, and stress state. Our models are designed for the large data sets of induced seismicity typically associated with geothermal sites. We produced labeled event catalogs and used them on geothermal data to assess how ML can facilitate geothermal production and exploration. All datasets are available on the GDR Productivity: The project produced 32 publications in peer reviewed journals (two are in review). It supported the work of 6 PhD students, 40 conference presentations, 6 keynote talks at national meetings, and mentoring and professional development for 4 postdoctoral fellows.

15 GEOTHERMAL ENERGY↗

Scalable preconditioning for the stabilized contact mechanics problem

We present a family of preconditioning strategies for the contact problem in fractured and faulted porous media. We combine low-order continuous finite elements to simulate the bulk deformation with piecewise constant Lagrange multipliers to impose the frictional contact constraints. This formulation is not uniformly inf-sup stable and requires stabilization. We improve previous work by Franceschini et al. (2020) by introducing a novel jump stabilization technique that requires only local geometrical and mechanical properties. We then design scalable preconditioning strategies that take advantage of the block structure of the Jacobian matrix using a physics-based partitioning of the unknowns by field type, namely displacement and Lagrange multipliers. The key to the success of the proposed preconditioners is a pseudo-Schur complement obtained by eliminating the Lagrange multiplier degrees of freedom, which can then be efficiently solved using an optimal multigrid method. Numerical results, including complex real-world problems, are presented to illustrate theoretical properties, scalability and robustness of the preconditioner. A comparison with other approaches available in the literature is also provided.

58 GEOSCIENCES↗

The Stability Transition From Stable to Unstable Frictional Slip With Finite Pore Pressure

Pore fluids are ubiquitous throughout the lithosphere and are commonly invoked as the cause of induced seismicity and slow earthquakes. We perform lab experiments to address these questions for drained fault conditions and low pore pressure. We shear simulated faults at effective normal stress (σ' n ) of 20 MPa and pore pressures Pp from 1 to 4 MPa. We document the full range of lab earthquake behaviors from slow slip to elasto-dynamic rupture and show that slow slip can be explained by the slip rate dependence of the critical rheologic stiffness without dilatancy hardening or other fluid effects. Our fault permeabilities ranges from 10 –18 to 10 –17 m 2 with an initial porosity of 0.1 and estimated fluid diffusion time ≈1 s. Slow slip and quasi-dynamic fault motion may arise from high P p at higher pressures but dilatancy strengthening is not a general requirement.

58 GEOSCIENCES↗

The High-Frequency Signature of Slow and Fast Laboratory Earthquakes

Tectonic faults fail through a spectrum of slip modes, ranging from slow aseismic creep to rapid slip during earthquakes. Understanding the seismic radiation emitted during these slip modes is key for advancing earthquake science and earthquake hazard assessment. In this work, we use laboratory friction experiments instrumented with ultrasonic sensors to document the seismic radiation properties of slow and fast laboratory earthquakes. Stick-slip experiments were conducted at a constant loading rate of 8 μm/s and the normal stress was systematically increased from 7 to 15 MPa. We produced a full spectrum of slip modes by modulating the loading stiffness in tandem with the fault zone normal stress. Acoustic emission data were recorded continuously at 5 MHz. We demonstrate that the full continuum of slip modes radiate measurable high-frequency energy between 100 and 500 kHz, including the slowest events that have peak fault slip rates <100 μm/s. The peak amplitude of the high-frequency time-domain signals scales systematically with fault slip velocity. Stable sliding experiments further support the connection between fault slip rate and high-frequency radiation. Experiments demonstrate that the origin of the high-frequency energy is fundamentally linked to changes in fault slip rate, shear strain, and breaking of contact junctions within the fault gouge. Our results suggest that having measurements close to the fault zone may be key for documenting seismic radiation properties and fully understanding the connection between different slip modes.

58 GEOSCIENCES↗

Rupture Dynamics Simulations of Shallow Crustal Earthquakes on Reverse Slip Faults

The main goal of this study was to perform dynamic rupture simulations in order constrain the characterization of kinematic rupture models implemented in techniques for simulating strong ground motion for crustal earthquakes. First, we investigated the rupture process of the 2008 Iwate-Miyagi Nairiku earthquake by developing a dynamic rupture model on a reverse fault using a trial-and-error technique that produced a slip distribution and near-fault ground motion that matched the recorded ones. The simulations were performed in the frequency range 0-2 Hz, using a 3D staggered grid finite-difference method and a linear slip weakening friction law. Constrained by the observed slip distribution and consistent with dynamic rupture models, the derived kinematic rupture model of the Iwate-Miyagi Nairiku earthquake contains areas with relatively large slip rate, representing strong motion generation areas, set against lower amplitude heterogeneous background slip, and a relatively low slip rate in the weak zone of the top-most crust (upper 3km). Second, we performed rupture dynamics modeling to constrain shallow slip characterization in rupture models for strike-slip crustal earthquakes. The objective was to establish general rules about the characterization of slip rate function and slip at shallow depths and in the strong motion generation areas (SMGAs). The simulations of spontaneous rupture were performed in the frequency range 0-2.0Hz, using a 3D staggered-grid finite-difference method and a layer over half-space 1D crustal velocity model with a minimum shear-wave velocity of 2.8 km/s. In order to account for changes in material ductility and reduction of stress drop, observed in the shallow crust (upper 3-5 km), and the transition from ductile state to brittle state in the upper seismogenic zone, in our stress models we included a shallow weak zone (<4km). In this zone the stress drop was set to zero at the free surface and gradually increased with depth, while the slip weakening distance was set to 75 cm at the free surface and decreased to 50 cm at the base of the weak zone. As in the case of reverse faulting during the Iwate-Miyagi Nairiku earthquake, from these computations we found a systematic change in the shape of the slip-rate function, from Kostrov-type in the deeper part of the fault to a more symmetric cosine-type in the upper few kilometers, near the free-surface. Moreover, the average slip duration in the weak zone, with respect to slip duration in the deeper parts of the fault, increases by at most a factor of 1.5. We found a systematic and gradual change in the shape of the slip-rate function from Kostrov-type in the asperity areas (SMGAs), to more symmetric cosine-type in the upper few km near the free-surface, on the long period motion generation area (LMGA). The effective rise time in the LMGA, with respect to that in the SMGAs, increases by at most a factor of 2. Effective rise time is the time difference between the time at which the slip rate drops to a level that is equal to 25% of its peak and the onset time of the signal. In addition, the slip in the LMGA located above the SMGAs is almost the same as the one in the SMGAs, and about 1.5 times larger than the average slip.

58 GEOSCIENCES↗

Failure prediction and estimation of failure parameters

Machine-learning methods and apparatus are disclosed to determine frictional state or other parameters in an earthquake zone or other failing medium, using acoustic emission, seismic waves, or other detectable indicators of microscopic processes. Predictions of future failures are demonstrated in different regimes. A classifier is trained using time series of acoustic emission data along with historic data of frictional state or failure events. In disclosed examples, random forests and gradient boost trees are used, and grid-search or EGO procedures are used for hyperparameter tuning. Once trained, the classifier can be applied to testing or live data in order to assess a frictional state, assess seismic hazard, or make predictions regarding a future failure event. The technology has been developed in a double direct shear apparatus, but can be widely applied to seismic faults, other terrestrial failures, or failures in man-made structures. Variations are disclosed.

Johnson, Paul Allan↗

Additive friction stir deposition of SS316: Effect of process parameters on microstructure evolution

Solid state nature of additive friction stir deposition (AFSD) additive manufacturing process is very advantageous in terms of defect formation and microstructural refinement in the material. Current study presents the process optimization, microstructural evolution and kinetics of recrystallization for AFSD deposited low stacking fault energy material - SS316. As deposited microstructure shows equiaxed ultra fine grains with an average grain size of ~5.0 ± 0.5 μm. Shear deformation at high temperature during processing leads to the operation of restoration mechanisms. Observation of necklace type microstructure in the as deposited SS316 is attributed to discontinuous dynamic recrystallization during processing. Recrystallization kinetics of the AFSD SS316 is characterized using Johnson-Mehl-Avarami-Kolmogorov (JMAK) model. Deformation – thermal cycling during AFSD process resulted in inconsistent recrystallization kinetics. Variation in strain, strain rate and temperature during processing with processing parameters result in varying microstructure and tool wear. Further, high strength-ductility combination and sustained work hardening in as deposited SS316 appear to arise from transformation and twinning during deformation, leading to the formation of hierarchical twins and martensitic phase after deformation.

36 MATERIALS SCIENCE↗

Frictional Behavior Downdip Along the Subduction Megathrust: Insights From Laboratory Experiments on Exhumed Samples at In Situ Conditions

Subduction megathrusts exhibit a spectrum of slip modes, including catastrophic earthquakes. Although the mechanical and frictional properties of materials sampled from subduction zones have been studied extensively, few datasets have been collected for compositions and at pressure and temperature conditions representative of those in situ. The Nankai subduction zone in southwest Japan is a well-studied margin, and abundant data provide an opportunity to advance our understanding of fault and earthquake physics. Here, we use samples exhumed in the Shimanto and Sanbagawa Belts on Shikoku Island of southwest Japan that represent analogs for materials along the present-day megathrust at depths of ~5–>25 km, and we shear these at their peak in situ pressure-temperature (P-T) conditions. Effective normal stresses range from 28 to 192 MPa, and temperatures from 105°C to 470°C. We used pore fluid pressures of 45–240 MPa, corresponding to fluid overpressure ratios λ of 0.65 and 0.90. Slip velocities of 0.1–100 μm/s were used, in order to focus on the nucleation of instability and earthquakes. We found predominantly velocity-strengthening (inherently stable) behavior under all conditions for λ = 0.65. For λ = 0.90, velocity-weakening behavior was observed at 350°C, with velocity-strengthening behavior at lower and higher temperatures. The rate/state frictional stability parameter (a-b) increases with slip velocity at temperatures up to ~200°C and remains constant or decreases with slip velocity at higher temperatures. Altogether, our results demonstrate the potentially important roles of both temperature and slip velocity in controlling the distribution of stress and frictional rheology along subduction thrusts.

58 GEOSCIENCES↗

Refinements to the Graves–Pitarka Kinematic Rupture Generator, Including a Dynamically Consistent Slip-Rate Function, Applied to the 2019 Mw 7.1 Ridgecrest Earthquake

The main objective of this study is to develop physics–based constraints on the spatiotemporal variation of the slip–rate function using a simplified dynamic rupture model. First, we performed dynamic rupture modeling of the 2019 M w 7.1 Ridgecrest, California, earthquake, to analyze the effects of depth–dependent stress and material friction on slip rate. Then, we used our modeling results to guide refinements to the slip–rate function that were implemented in the Graves–Pitarka kinematic rupture generation technique. The dynamic ruptures were computed on a surface–rupturing, planar strike–slip fault that includes a weak (negative to low–stress–drop) zone in the upper 4 km of the crust. Below the weak zone, we placed high–stress–drop patches designed to mirror the large–slip areas seen in various rupture model inversions of the event. The locations of the high–stress–drop patches and the hypocenter were varied in multiple realizations to investigate how changing the dynamic conditions affected the resulting rupture kinematics, in particular, the slip rate. From these simulations, we observed a systematic change in the shape of the slip–rate function from Kostrov type below the weak zone to a predominantly symmetric shape within the weak zone, along with a depth–dependent reduction of peak slip rate. We generalized these shallow rupture features into a depth–dependent parametric variation of the slip–rate function and implemented it in the Graves–Pitarka kinematic rupture model generator. The performance of the updated kinematic approach was then verified in 0–4 Hz simulations of the M w 7.1 Ridgecrest earthquake, which showed that incorporating the depth–dependent variation in the shape of the slip–rate function improves the fit to the observed near–fault ground motions in the 0.5–3 s period range.

58 GEOSCIENCES↗

Friction stir gradient alloying: A high-throughput method to explore the influence of V in enabling HCP to BCC transformation in a γ-FCC dominated high entropy alloy

The compositional possibilities in high entropy alloys (HEAs) is very vast and effective strategies are needed to establish potential alloy chemistries. In this study, friction stir gradient alloying (FSGA), a recently-introduced high-throughput (HT) technique incorporating compositional and microstructural gradients, was used to explore the possibility of introducing a bcc transformation domain in γ-fcc dominated TRIP HEA Fe 38.5 Mn 20 Cr 15 Co 20 Si 5 Cu 1.5 (at.%) by vanadium addition. FSGA has accelerated the process of attaining the composition-microstructure library and the results suggest that vanadium addition of ~1 at.% results in nucleation of α-bcc. Additionally, this is the first observation of ε-hcp to α-bcc transformation in a TRIP HEA and supports the Olson-Cohen model of martensitic transformation (γ-fcc → ε-hcp → α-bcc). ε-hcp is nucleated with a pyramidal orientation by the formation of planar faults on the {111} close-packed planes of γ-fcc, and α-bcc is nucleated in the vicinity of ε-hcp along the {111} γ-fcc trace. There are indications for the formation of basal orientated ε-hcp from intersecting ε-hcp planes orientated for the pyramidal slip. Thus, the study shows the nucleation and growth of the α-bcc phase attributed to the dual effects of chemistry and strain and potential co-existence of all the three phases.

36 MATERIALS SCIENCE↗

Magnitude and Timing of the Tiltill Rockslide in Yosemite National Park, California

Yosemite National Park, California, is one of the best-documented sites of historical rockfalls and other rock slope failures; however, past work shows that this record does not capture the infrequent largest occurrences, prehistoric events orders of magnitude larger than the largest historic ones. These large prehistoric events are evident as voluminous bouldery landslide deposits, permitting volume and age quantification to better understand local volume–frequency relationships, potential triggering mechanisms, and the hazard such events might pose. The Tiltill rockslide in northern Yosemite is one such example, consisting of 2.1 × 10 6 m 3 ± 1.6 × 10 6 m 3 of talus (1.5 × 10 6 m 3 original volume of rock mass) that slid across the floor of Tiltill Valley, partially damming Tiltill Creek to create a seasonal pond that drains through and around the rockslide mass. This volume and the rockslide's effective coefficient of friction, 0.47, place it near the boundary between long-runout landslides and ordinary Coulomb failure. Although the rockslide superficially appears to consist of two separate lobes, statistically indistinguishable 10 Be exposure dates from eight samples indicate a single event that occurred at 13.0 ± 0.8 ka. The age of the Tiltill rockslide and its relatively low elevation compared to equilibrium line altitudes at this place and time make glacial debutressing a highly unlikely triggering mechanism. Additionally, seismic shaking associated with fault rupture along the eastern Sierra Nevada is shown to be a plausible but unverified trigger.

10Be cosmogenic nuclide exposure dating↗

Spatio-Temporal Clustering of Seismicity Enabled by Off-Fault Plasticity

While significant progress has been made in understanding earthquake source processes in linear elastic domains, the effect of more realistic rheologies including plasticity is poorly understood. Here, in this work, we simulate the sequence of earthquake and aseismic slip of a 2D antiplane rate-and-state fault embedded in a full-space elastic-plastic bulk. We show that off-fault plasticity may lead to partial ruptures as well as temporal clustering of seismic events. Furthermore, the interaction of fault slip and off-fault plasticity results in pockets of slip deficit. While the energy dissipated through plastic deformation remains a small fraction of the total energy budget, its impact on the source characteristics is disproportionally large through the redistribution of stresses and viscous relaxation. Our results suggest a new mechanism of dynamic heterogeneity in earthquake physics that may have important implications on earthquake size distribution and energy budget.

58 GEOSCIENCES↗

Co-introduction of precipitate hardening and TRIP in a TWIP high-entropy alloy using friction stir alloying

Tuning deformation mechanisms is imperative to overcome the well-known strength-ductility paradigm. Twinning-induced plasticity (TWIP), transformation-induced plasticity (TRIP) and precipitate hardening have been investigated separately and have been altered to achieve exceptional strength or ductility in several alloy systems. In this study, we use a novel solid-state alloying method—friction stir alloying (FSA)—to tune the microstructure, and a composition of a TWIP high-entropy alloy by adding Ti, and thus activating site-specific deformation mechanisms that occur concomitantly in a single alloy. During the FSA process, grains of the as-cast face-centered cubic matrix were refined by high-temperature severe plastic deformation and, subsequently, a new alloy composition was obtained by dissolving Ti into the matrix. After annealing the FSA specimen at 900 °C, hard Ni–Ti rich precipitates formed to strengthen the alloy. An additional result was a Ni-depleted region in the vicinity of newly-formed precipitates. The reduction in Ni locally reduced the stacking fault energy, thus inducing TRIP-based deformation while the remaining matrix still deformed as a result of TWIP. Our current approach presents a novel microstructural architecture to design alloys, an approach that combines and optimizes local compositions such that multiple deformation mechanisms can be activated to enhance engineering properties.

36 MATERIALS SCIENCE↗

Downhole Sensing and Event-Driven Sensor Fusion for Depth-of-Cut Based Autonomous Fault Response and Drilling Optimization

Achieving robust and efficient drilling is a critical part of reducing the cost of geothermal energy exploration and extraction. Drilling performance is often evaluated using one or more of three key metrics: depth of cut (DOC), rate of penetration (ROP), and mechanical specific energy (MSE). All three of these quantities are related to each other. DOC refers to the depth a bit penetrates into rock during drilling. This is an important quantity for estimating bit behavior. ROP is the simply the DOC multiplied by the rotational rate, and represents how quickly the drill bit is advancing through the ground. ROP is often the parameter used for drilling control and optimization. Finally, MSE provides insight into drilling efficiency and rock type. MSE calculations rely on ROP, drilling force, and drilling torque. Surface-based sensors at the top of the drill are often used to measure all these quantities. However, top-hole measurements can deviate substantially from the behavior at the bit due to lag, vibrations, and friction. Therefore, relying only on top-hole information can lead to suboptimal drilling control. In this work, we describe recent progress towards estimating ROP, DOC, and MSE using down-hole sensing. We assume down-hole measurements of torque, weight-on-bit (WOB). Our hypothesis is that these measurements can provide more rapid and accurate measures of drilling performance. We show how a multi-layer perceptron (MLP) machine learning algorithm can provide rapid and accurate performance when evaluated on experimental data taken from Sandia’s Hard Rock Drilling Facility. In addition, we implement our algorithms on an embedded system intended to emulate a bottom-hole-assembly for sensing and estimation. Our experimental results show that DOC can be estimated accurately and in real-time. These estimates when combined with measurements for rotary speed, torque, and force can provide improved estimates for ROP and MSE. These results have the potential to enable better drilling assessment, improved control, and extended component lifetimes.

15 GEOTHERMAL ENERGY↗

Utah FORGE 5-2557: Trixial Shear and Wellbore Breakout Tests Report and Data

Included are a report and data from triaxial shear tests on simulated Sierra White fault gouge and a borehole stimulation experiment. This serves as a progress report for the Role of Fluid and Temperature in Fracture Mechanics and Coupled Thermo-Hydro-Mechanical-Chemical (THMC) Processes for Enhanced Geothermal Systems project by Purdue University (Utah FORGE R&D project 5-2557). Tests were done with Sierra White granite (SWG) to investigate slip mechanisms and to evaluate the fracture initiation and propagation from a geothermal borehole. Raw data is provided alongside a report on all experimental testing protocols, experimental design, and results.

15 GEOTHERMAL ENERGY↗