Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rock”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Using Cosmic Ray Muons to Assess Geological Characteristics in the Subsurface

Cosmic rays are energetic nuclei and elementary particles that originate from stars and intergalactic events. The interaction of these particles with the upper atmosphere produces a wide range of secondary particles that reach the surface of the earth, of which muons are the most prominent. With enough energy, muons can travel up to a few kilometers beneath the surface of the earth before being stopped completely. The terrestrial muon flux profile and associated zenith angle can be utilized to determine geological characteristics of a location (e.g., rock overburden and density) without having to use conventional methods such as boreholes. This work uses a low-power plastic scintillator-based muon detection system as a prototype for this non-destructive geological assay methodology. Four custom designed 102 cm x 51 cm x 5 cm plastic scintillation panels are used to realize two orthogonal detection planes. Optical photons from each scintillation panel are read using OnSemi J-Series 4x4 silicon photomultiplier (SiPM) arrays in conjunction with preamplifiers. Simultaneous triggers between detectors from two planes indicate a coincidence event which is recorded using the QuarkNet data acquisition system (DAQ) from Fermi National Accelerator Laboratory. A custom detector holder was designed to securely mount the detection system and rotate the panels along the zenith to collect data at variable angles. In order to quantify the systematic uncertainties associated with the detector, such as energy depositions and angular resolution of the detector design, a Monte Carlo (MC) simulation using Geant4 is being developed. Cosmic ray flux prediction will be included in the project by adding the CORSIKA MC code to the simulation toolchain. Simulated and experimental data will drive the development and validation of a reconstruction algorithm that, upon completion, is expected to predict average overburden and rock density. Extended detector exposure to muons can be used as a means to understand changes in the surrounding environment like rock porosity. On the experimental front, muons will initially be measured at the surface, establishing the baseline flux. This is followed by recording the muon flux at variable depths and zenith angles, where the data will be used by the reconstruction algorithm to predict the overburden. The result will be benchmarked against geological surveys. The measured flux data will also be used to benchmark independent and established models. Successful proof-of-concept demonstration of this technology can open doors for long term non-invasive geological monitoring. The detector design, experimental methodology, and the benchmarking efforts are detailed in this work.

Gadey, Harish Reddy↗
DOE OSTI2020

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Analcime-wairakite formation during experimental cement-bentonite alteration at 200–300 °C

Engineered barrier system materials in nuclear waste repositories may undergo hydrothermal alteration in response to groundwater saturation and heating events over their long operational timescales. Hydrothermal interactions between engineered materials (e.g., bentonite buffers, cements, waste canister materials) and the host rock environment may drive alteration processes that affect the advantageous properties of some barrier materials. However, such alteration will also promote the formation of zeolites, hydrothermally formed minerals that may themselves act to isolate radionuclides. Understanding the environmental conditions that control radionuclide-sorbing properties of zeolites, such as Si/Al ratio or Na content, will be valuable for assessing the changing properties of engineered barrier materials in the case of an in situ heating event. Here, we present experimental work characterizing the formation of zeolites during hydrothermal interactions between generic clay and cement barrier materials, with a focus on the analcime-wairakite zeolite series, which has known radionuclide sorption and exchange properties. We present the results of hydrothermal experiments combining uncured ordinary Portland cement powder with steel and Wyoming bentonite in the presence of Opalinus clay and a synthesized Opalinus clay groundwater to simulate water-saturated conditions in an argillaceous rock repository. The experiments were conducted isothermally at 200 or 300 °C for 8 to 24 weeks. In contrast to analogous studies that did not include cement reactants, we observed the formation of analcime-group minerals in all experiments. The addition of cement resulted in lower Si/Al ratios in the zeolites compared to similar studies that did not include cement. The presence of cement as a reactant was interpreted to promote analcime formation at 200 °C. At 300 °C, we observed higher calcium and silica concentrations in the aqueous solutions as well as increased wairakite formation and decreased analcime formation compared to the experiments at 200 °C. These results show a fully realized analcime-wairakite solid solution that falls between Si/Al = 2 (ideal) and a trend of analcime minerals that have increasing Si/Al ratios with increasing Na/(Na+Ca). Our results show predictable relationships between the Si/Al ratio and the analcime-wairakite content under hydrothermal conditions and illustrate that repository material interactions may promote the formation of zeolites in the analcime-wairakite solid-solution series during heating events in the subsurface.

58 GEOSCIENCES↗
DOE OSTI2025

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE↗
DOE OSTI2024

Coupled feldspar dissolution and secondary mineral precipitation in batch systems: 6. Labradorite dissolution, calcite growth, and clay precipitation at 60 °C and pH 8.2–8.4

Here, we conducted experiments on concurrent labradorite dissolution, calcite precipitation, and clay precipitation in batch reactor systems and tracked reaction processes using multiple isotope tracers. Labradorite was chosen for its role as a major and reactive component in basalt; the experiments thus directly impact our understanding of CO 2 storage in basalt aquifers and enhanced rock weathering. We doped initial solutions with 29 Si, 43 Ca, and Ca 13 CO 3 (s). Experiments were conducted at 60 °C and pH ~ 8.3 for up to 840 h, with isotope ratios in the experimental aqueous solutions measured using MC-ICP-MS. Unidirectional rates of labradorite dissolution near equilibrium were approximately two orders of magnitude slower than far-from-equilibrium rates reported in the literature. Calcite growth occurred near equilibrium and the rates were limited by the labradorite dissolution rates. In the steady state phase, the interplay of these three heterogeneous reactions—labradorite dissolution, calcite growth, and clay precipitation—results in a coupled system that approaches a near-equilibrium state. The system does not reach true equilibrium because labradorite continues to dissolve, albeit at a much slower rate near equilibrium. The overall reaction can be approximated as, Na 0.4 Ca 0.6 Al 1.6 Si 2.4 O 8 + 0.6HCO 3 - + 1·.7H 2 O + 0.4H + → 0.4Na + + 0.6CaCO 3(s) + 0.5Al 2 Si 2 O 5 (OH) 4(s) + 0.6Al(OH) 4 - + 1.4SiO 2 o (aq). The experimental results show that using short-term far-from-equilibrium rate constants would lead to an overestimation of feldspar weathering rates at the Earth’s surface (e.g., basalt weathering and enhanced rock weathering) and CO 2 mineralization in basalt aquifers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Critical review of the Mogi failure criterion based on true-triaxial laboratory data analysis and theoretical considerations

The influence of the intermediate principal stress, σ 2 , is disregarded when rock strength is analyzed using the Mohr-Coulomb or Hoek-Brown criteria. Several polyaxial strength criteria have been proposed to account for the influence of σ 2 . One of the methods that received considerable attention is the Mogi criterion. The criterion modifies the Nadai or Drucker-Prager approaches by exchanging the normal octahedral stress with what is introduced as the mean effective stress, σ m,2 . This change is based on experimental evidence that the influence of the intermediate principal stress on strength is less pronounced than the influence of the minimum principal stress. However, the Mogi criterion was derived somewhat empirically and consequently it is essential to analyze its potential theoretical deficiencies. In the current note, we analyzed the Mogi failure criterion, describing its characteristics in both τ oct –σ m,2 and σ 1 -σ 2 spaces. We investigated the relation between Mohr-Coulomb and Mogi criteria, which marks the basis for the commonly used Mogi-Coulomb approach. Here, we compared the performance of the Mogi criterion against available polyaxial strength data using three fitting approaches: Mogi-Coulomb based on triaxial compression results only, linear Mogi based on all polyaxial strength results, and power law Mogi also based on all data. The analysis reveals several serious deficiencies of the Mogi approach which all suggest that the Mogi criterion cannot properly represent the σ 2 dependence of rock strength in many cases.

58 GEOSCIENCES↗
DOE OSTI2026

In situ quantification of fracture slip induced by hydraulic injections in a deep borehole: A comparison of two different borehole techniques

In situ measurements of fracture deformation during fluid injection are rare, yet essential for understanding the mechanical response of fractured rock. In this study, we evaluate the reliability of two methods by comparing their slip vector estimates: high-resolution borehole acoustic televiewer images captured before and after injection tests, and displacement data from a three-component borehole deformation probe recorded during the injections. Acoustic televiewer images capture only final in-plane displacement, whereas three-component borehole deformation measurements provide full 3D, transient fracture movement. Four injection tests in a fractured granitic rock mass along an inclined borehole at the Bedretto Underground Laboratory (Switzerland), beneath more than 1100 m of overburden, were analyzed. The two methods yielded consistent kinematics and comparable slip magnitudes, typically in the range of 0.2 – 0.6 mm. Angular differences between estimated slip directions ranged from 0.2° to 28°. For two test intervals with pre-existing fractures optimally oriented for slip, the measured slip vectors agreed with forward-modeled slip vectors based on the regional in situ stress field, supporting the Wallace–Bott hypothesis. Measurable shear displacement was identified in two test intervals containing fractures oriented nearly perpendicular to the least principal stress, including a newly formed hydraulic fracture. This demonstrates that even fractures unfavorably oriented for slip can undergo shear due to local stress perturbations around an inclined borehole.

Acoustic televiewer↗
DOE OSTI2022

Method to account for natural fracture induced elastic anisotropy in geomechanical characterization of shale gas reservoirs

Shale has been usually recognized as a transverse isotropic (TI) medium in conventional geomechanical log interpretation due to its laminated nature. However, when natural fractures exist in the shale rock, additional elastic anisotropy is introduced, converting laminated Shale to an orthorhombic (OB) medium. Previous studies illustrate that neglecting the natural fracture induced anisotropy in shale geomechanical log interpretation could lead to inaccurate evaluations of elastic moduli and in-situ stresses. In this paper, a new method is developed to account for the natural fracture induced anisotropy in geomechanical log interpretation based upon the TI acoustic model developed by the author and a characterization technique of elastic wave anisotropy (Sayers, 1991). The new OB model incorporates the four acoustic log data inputs and five modeling constraints in a nonlinear optimization algorithm to solve for the nine independent stiffness coefficients of an OB rock, and further to solve for the geomechanical properties and in-situ stress profiles in an OB formation. The new method was validated with a Marcellus Gas Shale field case. Both the new OB model and the conventional TI model were applied to interpret the minimum horizontal stress profile for the same formation. By comparing the results, the OB model is more robust than the TI from two aspects. First, the average stress magnitude predicted by the OB model is closer to the one measured by the Diagnostic Fracture Injection Test (DFIT). Second, the OB model predicts a more obvious stress barrier between the lower Marcellus and upper Onondaga Limestone than the TI model does. Finally, the predicted stress barrier is consistent with the observation of the microseismic events of a horizontal well drilled and completed nearby, which reveals that no hydraulic fracture propagates downward through the bottom boundary of Marcellus Shale into the underlying Onondaga Limestone.

03 NATURAL GAS↗
DOE OSTI2024

Enhancement of disposal efficiency for deep geological repositories based on three design factors - Decay heat optimization, increased thermal limit of the buffer and double-layer concept

This study investigates the enhancement of disposal efficiency for deep geological repositories (DGRs) based on three design factors: decay heat optimization, increased thermal limit of the buffer, and double-layer concept using coupled thermo-hydro-mechanical (THM) numerical simulations. Decay heat optimization is achieved by iteratively emplacing spent nuclear fuels having the maximum and minimum decay heat in a canister. Disposal areas can be reduced by 20 % to 40 % compared to the current reference disposal system in Korea (KRS+) in accordance with the combinations of the three design factors, alleviating challenges in site selection for the DGR. This study additionally identifies an optimal layer spacing of 500 m for the double-layer concept in the viewpoint of the buffer temperature, where thermal interaction between the upper and lower layers nearly disappears. However, determining the ultimate disposal and layer spacing requires engineering judgement, considering not only the thermal performance of the DGR but also various factors such as cost and difficulties of the construction and rock mass stability. DGRs designed with an increased thermal limit of the buffer poses a greater probability of rock mass failure around disposal tunnels and deposition holes due to elevated thermal stresses. Densely arranged heat sources for the DGRs with enhanced disposal efficiency lead to larger temperature increase even at the far-field scale, raising a possibility of thermally driven fracture shear activation with associated hydraulic, mechanical, and seismic changes.

58 GEOSCIENCES↗
DOE OSTI2024

Influence of Annealing Temperature on the OER Activity of NiO(111) Nanosheets Prepared via Microwave and Solvothermal Synthesis Approaches

Earth-abundant transition metal oxides are promising alternatives to precious metal oxides as electrocatalysts for the oxygen evolution reaction (OER) and are intensively investigated for alkaline water electrolysis. OER electrocatalysis, like most other catalytic reactions, is surface-initiated, and the catalyst performance is fundamentally determined by the surface properties. Most transition metal oxide catalysts show OER activities that depend on the predominantly exposed crystal facets/surface structure. Therefore, the design of synthetic strategies to obtain the most active crystal facets is of significant research interest. In this work, rock salt NiO OER catalysts with (111) predominantly exposed facets were synthesized by a solvothermal (ST) method either heated under supercritical or microwave-assisted (MW) conditions. Particular emphasis was placed on the influence of the post annealing temperature on the structural configuration and OER activity to compare their catalytic performances. The as-prepared electrocatalysts are pure α-Ni hydroxides which were converted to rock salt NiO (111) nanosheets with hexagonal pores after heat treatment at different temperatures. The OER activity of the electrodes has been evaluated in 0.1 M KOH using geometric and intrinsic current densities via normalization by the disk area and BET area, respectively. The lowest overpotential at a geometric current density of 10 mA/cm 2 is found for samples pretreated by heating between 400 and 500 °C with a catalyst loading of 115 μg/cm 2 . Despite the very similar nature of the catalysts obtained from the two methods, the ST electrodes show a higher geometric and intrinsic current density for 500 °C pretreatment. The MW electrodes, however, achieve an optimal geometric current density for 400 °C pretreatment, while their intrinsic current density requires pretreatment over 600 °C. Interestingly, pretreated electrodes show consistently higher OER activity as compared to the poorly crystalline/less ordered hydroxide as-prepared electrocatalysts. Thus, our study highlights the importance of the synthesis method and pretreatment at an optimal temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Significance of Low‐Velocity Zones on Solute Retention in Rough Fractures

Natural fractures are characterized by high internal heterogeneity. This internal variability is the cause of flow channeling, which in turn leads to contaminant transport taking place primarily along the high-velocity channels. Mass exchange between the high-velocity channels and the low-velocity zones has the potential to enhance contaminant retention, due to solute diffusion into the low-velocity zones and subsequent exposure to additional surface area for diffusion into the bordering rock matrix. Here, we derive a random walk particle tracking method for heterogeneous fractures, which includes an additional term to account for the aperture gradient. The method takes into account advection, diffusion in the fracture and matrix diffusion. The developed numerical framework is applied to assess the effect of low-velocity zones in rough self-affine fractures. The results show that diffusion into low-velocity zones has a visible but modest impact on contaminant retention. The magnitude of this impact does not change considerably, regardless of whether diffusion into the rock matrix is considered in the model, and increases for a decreasing average Péclet number of the fracture.

58 GEOSCIENCES↗
DOE OSTI2024

4D Electrical Resistivity Imaging of Stress Perturbations Induced During High–Pressure Shear Stimulation Tests

Fluid flow through fractured media is typically governed by the distribution of fracture apertures, which are in turn governed by stress. Consequently, understanding subsurface stress is critical for understanding and predicting subsurface fluid flow. Although laboratory–scale studies have established a sensitive relationship between effective stress and bulk electrical conductivity in crystalline rock, that relationship has not been extensively leveraged to monitor stress evolution at the field scale using electrical or electromagnetic geophysical monitoring approaches. In this paper we demonstrate the use time–lapse 3–dimensional (4D) electrical resistivity tomography to image perturbations in the stress field generated by pressurized borehole packers deployed during shear–stimulation attempts in a 1.25 km deep metamorphic crystalline rock formation.

58 GEOSCIENCES↗
DOE OSTI2024

Electrical Resistivity Changes During Heating Experiments Unravel Heterogeneous Thermal‐Hydrological‐Mechanical Processes in Salt Formations

Abstract Rock salt is considered a suitable medium for the permanent disposal of heat‐generating radioactive waste due to its isolation properties. However, excavation damage and heating induce complex and heterogeneous thermal‐hydrological‐mechanical (THM) processes across different zones. Quantifying this heterogeneity is crucial for accurate long‐term performance assessment models, but traditional methods lack the necessary resolution. This study employs 4D electrical resistivity tomography (ERT) monitoring during controlled heating experiments in a salt formation to unravel the spatiotemporal dynamics of THM processes. Advanced time‐lapse inversion and clustering analysis quantify subsurface properties and map the heterogeneity of THM dynamics. The ERT results can estimate subsurface properties and delineate the damaged and intact zones, enabling appropriate parameterization and representation of processes for long‐term modeling. This approach may be used in further improving the predictive models and ensuring the safe long‐term disposal of radioactive waste in rock salt.

58 GEOSCIENCES↗
DOE OSTI2025

Clear as mud redefined: Tunable transparent mineral scaffolds for visualizing microbial processes below ground

Microbes inhabiting complex porous microenvironments in sediments and aquifers catalyze reactions that are critical to global biogeochemical cycles and ecosystem health. However, the opacity and complexity of porous sediment and rock matrices have considerably hindered the study of microbial processes occurring within these habitats. Here, we generated microbially compatible, optically transparent mineral scaffolds to visualize and investigate microbial colonization and activities occurring in these environments, in laboratory settings and in situ. Using inexpensive synthetic cryolite mineral, we produced optically transparent scaffolds mimicking the complex 3D structure of sediments and rocks by adapting a suspension-based, freeze-casting technique commonly used in materials science. Fine-tuning of parameters, such as freezing rate and choice of solvent, provided full control of pore size and architecture. The combined effects of scaffold porosity and structure on the movement of microbe-sized particles, tested using velocity tracking of fluorescent beads, showed diverse yet reproducible behaviors. The scaffolds we produced are compatible with epifluorescence microscopy, allowing the fluorescence-based identification of colonizing microbes by DNA-based staining and fluorescence in situ hybridization (FISH) to depths of 100 µm. Additionally, Raman spectroscopy analysis indicates minimal background signal in regions used for measuring deuterium and 13 C enrichment in microorganisms, highlighting the potential to directly couple D 2 O or 13 C stable isotope probing and Raman-FISH for quantifying microbial activity at the single-cell level. To demonstrate the relevance of cryolite scaffolds for environmental field studies, we visualized their colonization by diverse microorganisms within rhizosphere sediments of a coastal seagrass plant using epifluorescence microscopy. The tool presented here enables highly resolved, spatially explicit, and multimodal investigations into the distribution, activities, and interactions of underground microbes typically obscured within opaque geological materials until now.

36 MATERIALS SCIENCE↗
DOE OSTI2023

SSGF Mineral Major Elements and Lithium Concentration

Presented are major element and lithium concentrations of minerals from the Salton Sea Geothermal Field (SSGF), located in the Imperial Valley, California. With a recent increase in demand for lithium, the area is now being studied as a source of geothermal brine for lithium extraction. This study was performed to help quantify lithium storage in the SSGF brine and surrounding minerals. Rocks and brines sampled in this study are from surface rhyolitic domes on the South Eastern shore of the Salton Sea, the sedimentary and evaporitic rocks in the Durmid Hills, rhyolitic drill clippings previously studied by Schmitt & Hulen, commercial drill wells, and CA State 2-14 well.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Utah FORGE: Well 16B(78)-32 Drill Core Fracture Analysis Images and Data

This dataset contains drilling core data from well 16B(78)-32, including PDF documents with flattened core images annotated by feature type and core interval, as well as spreadsheets detailing feature morphologies by depth, planar feature measurements, and planar feature orientations rotated to in situ conditions. Core was recovered from three intervals, one per stimulation stage, in the crystalline rocks affected by the stimulation of well 16A(78)-32. Seven core runs were conducted, yielding 135.8 feet of recovered core. Features in the core were categorized into planar fractures, semi-planar fractures, unbroken mineralized fractures, rough fractures, curviplanar fractures, concave-convex surfaces, and planar compositional features such as mylonite or dike-like structures. Planar features were measured while the core was positioned horizontally, with the core axis aligned to a downhole azimuth of 42 degrees. Planar core measurements from stimulations 2 and 3 that could be confidently correlated with FMI data were rotated to in situ orientations. This was done by rotating the planes along vertical and horizontal axes to match the azimuth and inclination data recorded in the directional survey of well 16B(78)-32, as well as applying an axial rotation to resemble the fracture orientations observed in the FMI log at corresponding depths. Coherent sets of planar fracture measurements were made by aligning the core within each 3-foot section of the dissected core barrel, and between adjacent 3-foot sections within a core run by matching rock fabrics, saw cuts and/or tool marks. Where coherent fracture measurements could not be made within a core run, data sets are denoted by a subscript (i.e. 2-Ta and 2-Tb both come from tangent core run number 2).

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Simulation of emergent compaction banding fronts caused by frictional boundaries

This study examines the role of boundary friction in promoting heterogeneous compaction in soft rock specimens loaded at high confining pressure outside the domain of compaction localisation. An elastoplastic constitutive model characterised by tunable hardening/softening behaviour is used to conduct the analyses. Finite-element simulations suggest that material instability is a non-necessary condition for the emergence of compaction fronts. Such fronts propagated as a result of a severe deviation in the local responses induced by frictional constraints. Overall, these findings suggest that boundary effects can bias the assessment of the extent of the compaction localisation domain. Experimental countermeasures and informed model calibration procedures are therefore necessary to minimise such bias and enable more accurate predictions of soft rock compaction.

42 ENGINEERING↗