Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stress fields”

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 181 records · Page 10

Residual lattice strain in quartzites as a potential palaeo-piezometer

SUMMARY If a crystal lattice is subjected to a stress, it becomes distorted and no longer represents the ideal crystal symmetry, and if the stress introduces defects such as dislocations, some of this distortion is preserved after the applied stress is removed. In this study, we investigate lattice distortion in quartz at the micron scale with synchrotron X-ray Laue diffraction. From Laue images the local deviatoric strain tensor is derived and corresponding stresses are calculated based on elastic properties. The method is applied to metasedimentary quartzites from the Bergell Alps that were deformed at conditions of greenschist facies metamorphism. The residual palaeostrain is represented in maps of the deviatoric strain tensor components and with deviatoric strain axis pole figures. Data suggest overall shortening perpendicular to the schistosity plane but with considerable asymmetry relative to foliation and lineation, probably attributed to simple shear. Crystallographic pole figures from Laue diffraction agree with neutron diffraction and EBSD measurements and display quartz c-axes girdle distributions with maxima also perpendicular to schistosity. The method shows promise to be used as a palaeo-piezometer to unravel the stress field during tectonic deformation.

58 GEOSCIENCES↗

Nanometer-Scale Acoustic Wave Packets Generated by Stochastic Core-Level Photoionization Events

We demonstrate that the absorption of femtosecond hard x-ray pulses excites quasispherical, high-amplitude, and high-wave-vector coherent acoustic phonon wave packets using an all hard-x-ray pump-probe scattering experiment. The time- and momentum-resolved diffuse scattering signal is consistent with an ensemble of 3D strain wave packets induced by the rapid electron cascade dynamics following photoionization at uncorrelated excitation centers. We quantify key parameters of this process, including the localization size of the stress field and the photon energy conversion efficiency into elastic energy. The parameters are determined by the photoelectron and Auger electron cascade dynamics, as well as the electron-phonon interaction. In particular, we obtain the localization size of the observed strain wave packet to be 1.5 and 2.5 nm for bulk SrTiO 3 and KTaO 3 single crystals, respectively. The results provide crucial information on the mechanism of x-ray energy deposition into matter and shed light on the shortest collective length scales accessible to coherent acoustic phonon generation using x-ray excitation. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗

Inelastic relaxation processes in amorphous sodium silicates

During fracture amorphous oxides exhibit irreversible processes, including inelastic and nonrecoverable relaxation effects in the process zone surrounding the crack tip. Here, classical molecular dynamics simulations were used with a reactive forcefield to evaluate inelastic relaxation processes in five amorphous sodium silicate compositions. Overall, the 20% Na 2 O-SiO 2 (NS20) composition exhibited the most inelastic relaxation, followed by the 15% Na 2 O-SiO 2 (NS15) composition, the 25% Na 2 O-SiO 2 (NS25) composition, and finally the 10% (NS10) and 30% (NS30) Na 2 O-SiO 2 compositions. Coordination analysis of the Na + ions identified that during inelastic relaxation the Na+ ions were increasingly coordinated by nonbridging oxygens (NBOs) for the NS10 and NS15 compositions, which was supported by radial analysis of the O-Na-O bond angles surrounding the crack tip. Across the sodium silicate compositional range, two different inelastic relaxation mechanism were identified based on the amount of bridging oxygens (BOs) and NBOs in the Na+ ion coordination shell. At lower (NS10) and higher (NS30) sodium compositions, the entire structured relaxed toward the crack tip. In contrast at intermediate sodium concentrations (NS20) the Na + ion migrates toward the crack tip separately from the network structure. By developing a fundamental understanding of how modified silica systems respond to static stress fields, we will be able to predict how varying amorphous silicate systems exhibit slow crack growth.

36 MATERIALS SCIENCE↗

Deconstructing the Soil Microbiome into Reduced-Complexity Functional Modules

The soil microbiome is an invaluable component of the biosphere and critical for ecosystem functions, including biogeochemical cycling, soil-atmosphere gas exchange, degradation of toxic compounds, and promotion of plant growth and stress resistance/resilience. An improved understanding of the soil microbiome will help with predicting how these processes respond to external perturbations, and with harnessing beneficial aspects of the soil microbiome for agronomic applications such as crop amendments. However, the extensive taxonomic and functional diversity inherent within the soil microbiome hinders efficient analysis of this system. Microbial biodiversity in soils is orders of magnitude greater than other commonly studied systems such as the human gut microbiome (Blum, Zechmeister-Boltenstern, and Keiblinger 2019; Berendsen, Pieterse, and Bakker 2012). Thousands of microbial taxa may be found in a single gram of soil (Roesch et al. 2007), and high rates of gene flow and mutation further promote microbial diversification (Sergaki et al. 2018). Concomitantly, functional diversity in soil is similarly extensive. Soil is a heterogeneous mixture of microenvironments with defined physical and chemical attributes (Bach et al. 2018), within which numerous microbial guilds of distinct life-strategies and metabolic capacities can be found (Perez-Garcia, Lear, and Singhal 2016; H.-S. Song et al. 2014). Furthermore, the soil microbiome harbors a significant fraction of rare and/or quiescent members (Blagodatskaya and Kuzyakov 2013) that exist below the threshold of detection of current technologies. Such rare taxa may make significant contributions to process rates (Shade and Gilbert 2015; Dawson et al. 2017), but their scarcity complicates their identification and analysis. Finally, the vast amounts of information generated through holistic analyses of soil communities represent an intense computational burden (Scholz, Lo, and Chain 2012; Prosser 2015) that precludes assessment of the complete functional and taxonomic diversity contained in this ecosystem.

Naylor, Dan T.↗

Stress analysis using direct-S wavelets produced by a vertical vibrator

We have determined how to extract the azimuth of maximum horizontal stress (SHmax) in deep rocks by doing a simple 360° mathematical rotation of a downgoing direct-S wavelet generated at the baseplate of a surface-based vertical vibrator. We worked with direct-S wavelets that travel through stressed rocks to a deep horizontal vertical seismic profile (VSP) geophone. We find that the azimuth where a polarity reversal occurs in mathematical rotations of this downgoing direct-S wavelet defines the azimuth of the SHmax in the overlying rocks. We tested this direct-S wavelet rotation method for determining the SHmax azimuth at a site in the Illinois Basin using legacy VSP data acquired in 2013. SHmax azimuths indicated by this simple wavelet-rotation method were determined when vertical vibrators were stationed at zero offset, at far offset, and at different azimuths around a VSP receiver well. These VSP-based SHmax azimuths agreed with the azimuth of the SHmax found by traditional minifrac tests in the VSP receiver well. This simple VSP data analysis procedure for detecting the azimuth of maximum horizontal stress has never, to our knowledge, been reported or discussed in the geophysical literature. This technical finding should be of interest to the worldwide geophysical community, especially to people who need to monitor how stress fields shift when fluids are injected into, or extracted from, deep porous reservoirs.

Geochemistry & Geophysics↗

Large Plate Experiment of Chloride-induced Stress Corrosion Cracking in Spent Nuclear Fuel Storage Canisters

This report provides a description of the testing (in progress at the Savannah River National Laboratory) to evaluate the corrosion behavior of a large section of welded plate material harvested from the Sandia National Laboratories’ full-size mockup of a spent nuclear fuel (SNF) dry storage canister. The details of the establishment of this plate test, the interim corrosion results (from visual examination), and the general plan for examination of the plate in CY21 to evaluate crack growth behavior (if any) are provided. This work is intended to inform aging management programs for SNF canisters under conditions of extended dry storage. A “large plate” (approximately 51 x 46 cm with thickness 16 mm and weight about 30 kg) and containing a circumferential weldment from the mockup is the test specimen. The knowledge of the salt, temperature, humidity conditions for chloride-induced stress corrosion cracking (CISCC) including the threshold stress intensity factor, K iscc , were used to setup the large plate test in which a set of part-through-wall and through-wall electrical discharge machined (EDM) defects were positioned in the weld residual stress field of the welded plate material from the mockup canister. A marine salt solution was applied to this plate specimen with an air-brush spray, then dried, and then exposed to a sustained 71-73% RH condition at an ambient room temperature of 22 °C. The salt loading and environmental condition is expected to be aggressive to cause CISCC

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Assessing Low-Temperature Geothermal Play Types: Relevant Data and Play Fairway Analysis Methods

The U.S. Department of Energy (DOE) Geothermal Technologies Office (GTO) is supporting the Geothermal Heating and Cooling Geospatial Datasets and Analysis project conducted by the National Renewable Energy Laboratory (NREL) as part of a broader effort to demonstrate the multi-faceted value of integrating geothermal power and geothermal heating and cooling (GHC) technologies into national decarbonization plans and community energy plans. Currently, there is a need to establish baseline low-temperature geothermal resource datasets and evaluate methods for deploying these technologies to provide the basis for supporting private sector investment. This project is focused on collecting baseline datasets, updating conceptual models, and creating Play Fairway Analysis (PFA) workflows for low-temperature (<150 degrees Celsius) geothermal resources of different geothermal play types (i.e., sedimentary basin, orogenic belts, and radiogenic geothermal play types) that could be used for geothermal heating and cooling (GHC), combined heat and power (CHP), and other geothermal direct uses (GDU) applications. Low-temperature geothermal resources are defined as reservoirs - natural or engineered - with temperatures <150 degrees Celsius. While the focus in the NREL effort is on GHC, resources at the upper end of this temperature range can also be used for small-scale power generation. This project does not include Ground Source Heat Pumps (GSHPs) technologies because they can be effectively developed almost anywhere. Low-temperature geothermal resources have not been studied as extensively as higher- to medium-temperature geothermal resources, but there is recent interest in improving understanding of these types of resources with an uptick of interest in geothermal technologies for decarbonizing heating and cooling systems. In addition, Enhanced Geothermal Systems (EGS) and other emerging technologies for exploiting petrothermal resources have opened the possibility of utilizing deep sedimentary basin systems, where porous media provide permeability and high temperatures can be reached at great depths. This project takes the approach of classifying low- temperature geothermal resources by geothermal play type (GPT). We defined and characterized three major classes of low-temperature GPT: sedimentary basins, orogenic systems, and radiogenic systems. We develop methodologies for evaluating and analyzing the potential for these resources building off the PFA approach to de-risking geothermal exploration and characterization. The proposed PFA approach for low-temperature geothermal resources includes: 1) identifying relevant data (e.g., datasets such bottom-hole temperatures from oil and gas wells, heat flow data, Quaternary faults and stress field data, geophysical data, etc.); 2) grouping and weighting of relevant datasets into PFA criteria (e.g., geological, risk, and economic criteria); 3) uncertainty quantification; 4) developing favorability or common risk maps for low-temperature geothermal resources to identify potential locations for more focused data collection; and 5) estimating electric power generation and heating potential at those locations using the GeoRePORT Resource Size Assessment Tool (RSAT). This project should facilitate future deployment of GHC, CHP, and GDU by providing data, tools, and a workflow applicable to low-temperature geothermal resources. Increased deployment of GHC and GDU will help achieve national and local decarbonization goals.

15 GEOTHERMAL ENERGY↗

Quantification of the Crack Evolution Process by Extracting Relevant Signal Components from Wave Propagation and Diffusive Transport Front Measurements

Wave propagation and diffusive transport phenomena in a geological rock sample undergoing crack evolution process are expected to interact with the mechanical discontinuities in the medium. The measurements of the signals associated with these phenomena can be used to assess and monitor the crack-driven micromechanical alterations in the rock. Different wave/diffusion phenomena, such as sonic propagation, pressure diffusion, and acoustic emission (AE), are sensitive to different elements of the mechanical discontinuities generated during the evolution of the crack clusters from initiation to coalescence. Sonic propagation, AE, and pressure diffusion monitoring have the potential to map the crack evolution because the transmitter-receiver arrays can be designed, arranged and tuned to (1) achieve maximum recovery of the scattered waveforms and travel times, (2) capture the later arrivals and multiple reflections, and (3) illuminate large rock volume. However, the structural/topological complexities of the mechanical discontinuities, complex distribution of the stress fields, complex mechanical alterations in media, and fluid redistribution in the crack system pose serious challenges for the detection and modeling of the crack evolution process (from here on, we will use the term ‘crack evolution process’ to mean that the crack evolution occurred under shallow crustal conditions). For purposes of accurately accounting such complexities and heterogeneities in the absence of reliable physical laws, simulation methods, and signal processing techniques, my early-career research proposal will develop and apply novel data-driven machine learning methods to: (1) extract signal components relevant to the various phases of crack evolution and (2) generate a 2D visual map of the crack evolution process.

58 GEOSCIENCES↗

Fracture Network Quantification during CO2 Injection

This is the presentation prepared for the ARMA 2025 (59th US Rock Mechanics/Geomechanics Symposium) Conference held in Santa Fe, New Mexico, June 8-11, 2025. Accurate mapping and quantification of these networks are essential to ensure the integrity of CO2 storage reservoirs, understand and reduce potential leakage, and maintain long-term environmental safety. This study presents a novel machine learning-driven approach, integrated with geomechanical analysis, to quantify fracture networks and assess their spatial distribution during CO2 injection. This paper combines microseismic monitoring data with principles of hydraulic diffusivity and geomechanical analysis to characterize reservoir scale fracture network. The novelty of our approach lies in its capacity to assimilate time-dependent pressure data and microseismicity into a cohesive framework, which not only identifies microseismic triggering fronts but also tracks fracture distribution during active injection. Besides, leveraging image log data and analysis our approach also provides another angle of the insights to solidate the fracture networks understanding and geomechanical impacts. Key results from our study include the detection of over 100 distinct fracture clusters across the injection site, with fracture orientations strongly correlated with the prevailing in-situ stress field.

CO2 storage and sequestration↗

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Geothermal Deep Direct-Use for Low-Carbon Heating: A Case Study at Cornell University

Cornell University is developing a geothermal deep direct-use system to provide reliable, low carbon, and cost-competitive heating for its main campus in Ithaca, New York. This geothermal system, also referred to as Earth Source Heat, would provide the majority of the campus heating load and represents a key pillar for the university to become carbon neutral by 2035. Over the last 10 years, numerous desktop, seismic, and geophysical studies have been undertaken at Cornell to characterize the local subsurface, predict long-term reservoir performance, design the integration into the existing district heating network, and estimate capital and levelized costs. Recently, the university received funding for drilling an observatory well to obtain in-situ subsurface data (e.g., reservoir temperatures, pressures, fluid chemistry, stress field), significantly lowering subsurface uncertainty and de-risking the project. Drilling of the test well is planned for fall 2021. This paper provides an overview and status of the project, including a summary of past studies and a description of the proposed well and planned testing activities.

Cornell University↗

Fracture Network Quantification during CO2 Injection

This is the conference paper accompanying an oral presentation at the ARMA 2025 (59th US Rock Mechanics/Geomechanics Symposium) Conference held in Santa Fe, New Mexico, June 8-11, 2025. Accurate mapping and quantification of these networks are essential to ensure the integrity of CO2 storage reservoirs, understand and reduce potential leakage, and maintain long-term environmental safety. This study presents a novel machine learning-driven approach, integrated with geomechanical analysis, to quantify fracture networks and assess their spatial distribution during CO2 injection. This paper combines microseismic monitoring data with principles of hydraulic diffusivity and geomechanical analysis to characterize reservoir scale fracture network. The novelty of our approach lies in its capacity to assimilate time-dependent pressure data and microseismicity into a cohesive framework, which not only identifies microseismic triggering fronts but also tracks fracture distribution during active injection. Besides, leveraging image log data and analysis our approach also provides another angle of the insights to solidate the fracture networks understanding and geomechanical impacts. Key results from our study include the detection of over 100 distinct fracture clusters across the injection site, with fracture orientations strongly correlated with the prevailing in-situ stress field.

CO2 storage and sequestration↗

Fuel performance analysis of fully-resolved TRISO compact

The TRi-structural ISOtropic (TRISO) fuel multilayered coating structure offers multiple barriers to fission product release, enhancing safety and performance. The heterogeneous nature of TRISO fuel compacts, comprising thousands of randomly distributed coated fuel particles embedded in a graphite matrix, creates intricate stress fields and thermal gradients that cannot be accurately modeled using simplified one-dimensional or homogenized approaches. Consequently, three-dimensional modeling enables the prediction of fuel compact dimensional changes, internal pressure buildup, and fission product transport pathways under diverse irradiation and thermal conditions. This capability facilitates detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating failures, which directly impact fuel performance and safety margins. This capability is particularly critical for advanced reactors, such as high-temperature gas-cooled reactors and other Generation IV reactor designs where TRISO fuel operates at elevated temperatures and burn-up levels. This work introduces a novel method to generate an optimized packing of TRISO compacts and a complete 3D mesh with random distribution of TRISO particles, which are discretized into each coating component layer.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Unraveling magneto-structural coupling of Ni 2 MnGa alloy under the application of stress and magnetic field using in situ polarized neutron diffraction

The magneto-structural evolution of a Ni 2 MnGa single crystal was investigated by in situ polarized neutron diffraction under both magnetic field and mechanical loading. The ability to separate the nuclear and magnetic scatterings by polarized neutrons enables simultaneous observations of the twin reorientation (through twin boundary motion) and magnetic moment configuration. It is found that under a 1.0 T saturating magnetic field, twin reorientation can be activated by compressive stress and is accompanied by a realignment of the magnetic moment, which follows the easy axis of the twin variant. On the other hand, compressive stress as small as 1.0 MPa can freeze the twin reorientation. The evolution of the magnetic scattering amplitudes suggests that the change in magnetization was mainly achieved through the magnetic domain wall motion, while the twin boundary motion was blocked.

36 MATERIALS SCIENCE↗

Stress-hybrid virtual element method on six-noded triangular meshes for compressible and nearly-incompressible linear elasticity

In this paper, we present a first-order Stress-Hybrid Virtual Element Method (SH-VEM) on six-noded triangular meshes for linear plane elasticity. Here, we adopt the Hellinger–Reissner variational principle to construct a weak equilibrium condition and a stress based projection operator. In each element, the stress projection operator is expressed in terms of the nodal displacements, which leads to a displacement based formulation. This stress-hybrid approach assumes a globally continuous displacement field while the stress field is discontinuous across each element. The stress field is initially represented by divergence-free tensor polynomials based on Airy stress functions, but we also present a formulation that uses a penalty term to enforce the element equilibrium conditions, referred to as the Penalty Stress-Hybrid Virtual Element Method (PSH-VEM). Numerical results are presented for PSH-VEM and SH-VEM, and we compare their convergence to the composite triangle FEM and B-bar VEM on benchmark problems in linear elasticity. The SH-VEM converges optimally in the L 2 norm of the displacement, energy seminorm, and the L 2 norm of hydrostatic stress. Furthermore, the results reveal that PSH-VEM converges in most cases at a faster rate than the expected optimal rate, but it requires the selection of a suitably chosen penalty parameter.

42 ENGINEERING↗

Utah FORGE 2-2439: A Multi-Component Approach to Characterizing In-Situ Stress: Laboratory, Modeling and Field Measurement - Workshop Presentation

This is a presentation on A Multi-Component Approach to Characterizing In-Situ Stress at the U.S DOE FORGE EGS Site: Laboratory, Modeling and Field Measurement project by Battelle [Columbus, OH], presented by Mark Kelley. The project's objective was to characterize stress in the Utah FORGE EGS reservoir using three methods: a laboratory rock-core stress estimation combined with a Machine Learning approach for estimation of in-situ stress from field sonic-log data, a field based in-situ measurement (min-frac) approach, and a modeling approach. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 7, 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↗