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 415 records · Page 23

Hydraulic Modeling of Induced and Propagated Fractures: Analysis of Flow and Pressure Data From Hydromechanical Experiments in the COSC-1 Deep Borehole in Crystalline Rock Near Åre, Sweden

To characterize the coupled hydromechanical behavior of rock fractures, the step-rate injection method for fracture in-situ properties (SIMFIP) was conducted with a specialized downhole probe developed by Guglielmi et.al. In June 2019, a field campaign was carried out near Åre, Sweden, where the SIMFIP probe was applied in the Collisional Orogeny in the Scandinavian Caledonides-1 scientific borehole to understand the dynamics of injection-induced fracture initiation, fracture opening, and shearing due to water injection-withdrawal in a borehole interval isolated by two packers. Three intervals were investigated at ~500 m depth: (a) an unfractured section (intact rock), (b) a section with non-conductive fractures, and (c) a section with hydraulically conductive fractures. Pressure, injection flow rate, and borehole wall displacement were simultaneously measured during the tests. In the present study, the geometry of the induced fracture and deformation of existing fractures at different time stages of the tests are determined based on a hydrologic model by using the measured pressure and flow data during each time stage of the experiment. A numerical model for the fluid flow within the fracture and the packed-off borehole interval is implemented within COMSOL Multiphysics. By matching model simulations with observed data for all three sections, estimates of the induced and propagated fractures' radius and aperture at successive time stages have been obtained in each case. We could also determine the non-linear relationship between fracture aperture and pressure for values above fracture opening pressures. The model results provide insights for the understanding of pressure-induced fracture initiation and propagation in crystalline rock.

58 GEOSCIENCES↗

Viscous Behavior of Clay‐Rich Rocks and Its Role in Focused Fluid Flow

Abstract Focused fluid flow is common in sedimentary basins worldwide, where flow structures often penetrate through sandy reservoir rocks, and clay‐rich caprocks. To better understand the mechanisms forming such structures, the impacts of the viscoelastic deformation and strongly nonlinear porosity‐dependent permeability of clay‐rich materials are assessed from an experimental and numerical modeling perspective. The experimental methods to measure the poroviscoelastic and transport properties of intact and remolded shale have been developed, and the experimental data is used to constrain the numerical simulations. It is demonstrated that viscoelastic deformation combined with nonlinear porosity‐dependent permeability triggers the development of localized flow channels, often imaged as seismic chimneys. The permeability inside a channel increases by several orders of magnitude compared to the background values. In addition, the propagation time scale and the channel size strongly depend on the material properties of the fluid and the rock. The time‐dependent behavior of the clay‐rich rock may play a key role in the long‐term integrity of the subsurface formations.

58 GEOSCIENCES↗

Measuring Fault Zone and Host Rock Hydraulic Properties Using Tidal Responses

Pore pressure diffusion along faults influences induced seismicity and rupture mechanics. In situ hydraulic diffusivity measurements along faults are rare and generally lower than inferred from seismicity migration. Here we use the tidal response of deep geothermal boreholes to measure fault diffusivity and permeability. Initial interpretations of the observation with a homogeneous confined aquifer model result in diffusivities of 10 -3 –10 -1 m 2 /s. However, this model mixes signals from both the conduit and the host rock. We develop a model for tidal response with a fault passing through the aquifer based on the fault-guided fracture network and solve for hydraulic properties in both the fault and the host rock. The resulting fault permeability is 2 × 10 -14 –7 × 10 -14 m 2 (90% confidence interval[CI]) and fault diffusivity is 0.08–0.33 m 2 /s (90% CI), which is 2 orders of magnitude higher than the host rock diffusivity in some wells, thus highlighting the role of faults as fluid conduits.

58 GEOSCIENCES↗

Nucleation-promoting and growth-limiting synthesis of disordered rock-salt Li-ion cathode materials

Disordered rock-salt oxides and oxyfluorides are promising positive electrode materials for high-performance lithium-ion batteries free of nickel and cobalt. However, conventional synthesis methods rely on post-synthesis pulverization to achieve cycling-appropriate particle sizes, offering limited control over particle microstructure and crystallinity. This accelerates degradation and complicates secondary particle processing. Here we present a synthesis strategy that enhances nucleation while suppressing particle growth and agglomeration across various disordered rock-salt compositions, including lithium–manganese–titanium oxide, lithium–manganese–niobium oxide, and lithium–nickel–titanium oxide systems. Applied to Li 1.2 Mn 0.4 Ti 0.4 O 2 , this method yields highly crystalline, well-dispersed sub-200 nm particles that form homogeneous electrode films with stable cycling behavior. Tested in cells with lithium metal as the counter electrode, these electrodes deliver ~200 mAh/g with 85% capacity retention relative to the first cycle after 100 cycles (20 mA/g, 1.5–4.8 V), and an average discharge voltage loss of 4.8 mV per cycle, compared to 38.6% retention and 7.5 mV loss per cycle for electrodes derived from pulverized solid-state particles. This approach suggests a route to enhance the performance and durability of disordered rock-salt electrodes for sustainable lithium-ion batteries.

Batteries↗

Miscibility of rock and ice in the interiors of water worlds

Abstract Super-Earths and sub-Neptunes are the most common planet types in our galaxy. A subset of these planets is predicted to be water worlds, bodies that are rich in water and poor in hydrogen gas. The interior structures of water worlds have been assumed to consist of water surrounding a rocky mantle and iron core. In small planets, water and rock form distinct layers with limited incorporation of water into silicate phases, but these materials may interact differently during the growth and evolution of water worlds due to greater interior pressures and temperatures. Here, we use density functional molecular dynamics (DFT-MD) simulations to study the miscibility and interactions of enstatite (MgSiO 3 ), a major end-member silicate phase, and water (H 2 O) at extreme conditions in water world interiors. We explore pressures ranging from 30 to 120 GPa and temperatures from 500 to 8000 K. Our results demonstrate that enstatite and water are miscible in all proportions if the temperature exceeds the melting point of MgSiO 3 . Furthermore, we performed smoothed particle hydrodynamics simulations to demonstrate that the conditions necessary for rock-water miscibility are reached during giant impacts between water-rich bodies of 0.7–4.7 Earth masses. Our simulations lead to water worlds that include a mixed layer of rock and water.

79 ASTRONOMY AND ASTROPHYSICS↗

Emerging investigator series: kinetics of diopside reactivity for carbon mineralization in mafic–ultramafic rocks

The ongoing use of fossil fuels to supply modern energy demands has necessitated research on combating carbon dioxide (CO 2 ) emissions and climate change. Carbon storage via mineral trapping in basalt and related rocks is a promising strategy. However, mineralization rates depend on the variable minerology that makes up these rock formations. Diopside (CaMgSi 2 O 6 ) is a common pyroxene mineral in ultramafic and mafic rocks including basalt, but relatively little work has been done to understand its carbon mineralization kinetics using hydrated supercritical CO 2 , which induces the formation of reactive nanoscale interfacial water films. Here, in situ XRD experiments at 50–110 °C and 90 bar indicate that diopside transforms into a myriad of Mg/Ca carbonates, including huntite [Mg 3 Ca(CO 3 ) 4 ] and very high magnesium calcite (VHMC, i.e., protodolomite). Through ex situ characterization, we were able to constrain reaction pathways for the dissolution–precipitation transformation process including metastable intermediate precipitates. Experiments performed at variable temperatures enabled Avrami-derived rate constants and an apparent activation energy of 97 ± 16 kJ mol –1 , implying the dissolution of diopside is the rate-limiting step. Density functional theory (DFT) calculations, used to gain molecular insight into the surface stability of the diopside during dissolution, suggest that exposed calcium cations are susceptible to dissolution when put in contact with water given their coordination environment. The collective results point to the high CO 2 mineralization potential of diopside in basalts, which could help guide parameterization of reactive transport models needed to design and permit commercial-scale subsurface carbon storage operations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Epitaxial growth of rock salt MgZrN 2 semiconductors on MgO and GaN

Ternary nitride compound semiconductors have attracted recent attention as electronic materials since their properties can be tuned by cation stoichiometry and ordering. A recently discovered example is MgZrN 2 , a ternary analog to the rock salt semiconductor ScN. MgZrN 2 has a larger bandgap and stronger dielectric response than the binary compound. Polycrystalline thin films of MgZrN 2 have been studied, but demonstration of high-quality growth is still required to establish its suitability for technological applications. Here, we report on epitaxial growth of MgZrN 2 thin films on (100) and (111) MgO substrates and (001) GaN templates. The MgZrN 2 composition is confirmed by Rutherford backscattering spectrometry, showing no oxygen in the film except for a thin surface oxide layer. Epitaxial growth results in MgZrN 2 with x-ray diffraction rocking curves with a full-width at half-maximum in the range of 0.3–3.0°, depending on the substrate. Transmission electron microscopy analysis of the MgZrN 2 film grown on a (111) MgO substrate confirms epitaxial growth and shows a sharp film/substrate interface. In-plane temperature-dependent Hall effect measurements show that the material is an n-type semiconductor with a relatively high concentration ( n 300K ≈ 10 19 –10 20 cm -3 ) of thermally activated electrons. Room-temperature transport measurements show a conductivity of 25 S cm -1 and a Seebeck coefficient of -80 μ V K -1 . Overall, these results provide an important step toward integration of rock salt MgZrN 2 with other technological nitrides for device applications.

36 MATERIALS SCIENCE↗

On various mechanisms of shear wave generation from underground chemical explosions in hard rocks

Significant shear waves observed during underground explosions can make them look similar to natural earthquakes. In the near field, they can be caused by the source geometry and emplacement conditions, heterogeneous and anisotropic properties of the rock mass, relaxation of a non-uniform in situ stress around the source after the departure of the shock wave as well as by activation of nearby faults. Here, the stress relaxation mechanism is important at depth when the difference in main principal stresses becomes significant. Another mechanism of shear motion generation is due to explosive gas flow from the source cavity into the cracks initiated by the shock wave and joints connected to the cavity. Also, the shock wave from the explosion can cause weakening of the rock either due to permanent damage or shock-induced liquefaction at saturated joints. Numerical modelling of underground chemical explosions in hard rock at various confinement stresses has been conducted to compare various mechanisms of shear wave generation and their seismic characteristics with that of a simple fault rupture process observed during an earthquake.

58 GEOSCIENCES↗

An integrated experimental–modeling approach to identify key processes for carbon mineralization in fractured mafic and ultramafic rocks

Abstract Controlling atmospheric warming requires immediate reduction of carbon dioxide (CO2) emissions, as well as the active removal and sequestration of CO2 from current point sources. One promising proposed strategy to reduce atmospheric CO2 levels is geologic carbon sequestration (GCS), where CO2 is injected into the subsurface and reacts with the formation to precipitate carbonate minerals. Rapid mineralization has recently been reported for field tests in mafic and ultramafic rocks. However, unlike saline aquifers and depleted oil and gas reservoirs historically considered for GCS, these formations can have extremely low porosities and permeabilities, limiting storage volumes and reactive mineral surfaces to the preexisting fracture network. As a result, coupling between geochemical interactions and the fracture network evolution is a critical component of long-term, sustainable carbon storage. In this paper, we summarize recent advances in integrating experimental and modeling approaches to determine the first-order processes for carbon mineralization in a fractured mafic/ultramafic rock system. We observe the critical role of fracture aperture, flow, and surface characteristics in controlling the quantity, identity, and morphology of secondary precipitates and present where the influence of these factors can be reflected in newly developed thermo-hydro-mechanical–chemical models. Our findings provide a roadmap for future work on carbon mineralization, as we present the most important system components and key challenges that we are overcoming to enable GCS in mafic and ultramafic rocks.

58 GEOSCIENCES↗

Hydrologic Impacts of a Strike-Slip Fault Zone: Insights from Joint 3D Body-Wave Tomography of Rock Valley

The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. Here, we manually picked ~17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (V P ), S-wave velocity (V S ), and V P /V S ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ~700 m in areas and generally align with geologic and structural expectations. V P and V S are relatively reduced near mapped and inferred faults, with the most prominent lower V P and V S zone around the densest collection of faults. We image V P /V S ratios ranging from ~1.5 to ~2.4, the extremes of which occur at a depth of ~100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ~100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher V P /V S areas and drier rock in lower V P /V S areas. A relatively higher V P /V S area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.

58 GEOSCIENCES↗

Determining Stress Orientation in Rock Valley, Nevada, Using Ambient Seismic Noise

The stress field and the mechanical properties of rocks are important to consider for nuclear explosion monitoring due to their effect on seismic wave radiation from earthquakes and explosions. At the Rock Valley Direct Comparison site, the regional orientation of the maximum horizontal compressive stress (SH max ) is well constrained, but it is unknown whether there are local heterogeneities. Here, I show that stress-induced anisotropy in nonlinear elasticity can be used to estimate the orientation of SH max . Rocks have compliant internal contacts, such as fractures and mineral grain boundaries, that respond to applied strains more strongly than individual mineral crystals. This strain response is asymmetric between compressive and dilatational strains and is affected by anisotropy in the ambient stress field. Traditional seismic velocity measurements are of strain-averaged velocity, which is less sensitive than nonlinear elasticity to the behavior of compliant internal contacts and therefore the stress field and fracture behavior. My results show that the orientation of SH max measured using ambient seismic noise aligns with regional estimates, with some heterogeneity that can also be explained by the limitations of the seismic array. I demonstrate the potential of a passive technique for monitoring the stress field in places that currently lack local measurements.

58 GEOSCIENCES↗

Rock Fracturing Using High-Pressure Ethylene/Nitrous Oxide Detonations

The present work investigates high initial pressure detonations of a stoichiometric mixture of ethylene and nitrous oxide (C2H4 + 6N2O) as a method of fracturing rock beneath the ground surface. These tests were conducted at a test site operated by the Energetic Materials Research and Testing Center (EMRTC), Socorro, New Mexico. The volume under the surface used for testing (called the Down Hole Assembly) consists of a 0.438 in. ID x 50 ft. long stainless-steel tube running down from the test site to a well bore which is 3 in. ID x 10 ft. long and the rock in the well bore is exposed to the propagating combustion wave. The testing carried out at Zucrow Laboratories in the smaller, alloy steel combustion vessel provided a scaling of pressures expected in the well bore. The combustion is initiated by energizing an EBW (Exploding Bridge Wire) above the ground surface. The experimental setup accommodates one high pressure (100,000 psia) transducer to measure the pressure peak and is placed approximately 5 ft. above the ground surface and 5 ft. downstream of the EBW. The focus of this series of experiments is to investigate the dependence of fracture to the rock beneath the surface on initial pressures of the mixture of ethylene and nitrous oxide. Experiments were carried out at initial pressures varying between 125 psia and 300 psi. The transducer recorded elevated pressures, which were 2.3 to 2.6 times in excess of the CJ values. The experimental results are discussed and explained in this report.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Legacy surface change analyses of the 1993 Rock Valley earthquake sequence for direct comparison to planned NA-22 underground conventional high-explosive experiments (Source Physics Experiment 3 (SPE3) - RV/DC Task 1.3 FY22 Final Report)

Recent work under two previous phases of the NNSA NA-22 Source Physics Experiment (SPE) have shown that underground chemical high-explosive experiments can produce detectable surface changes that differ in spatial extent and vertical magnitude depending on the geologic media at the site (Schultz-Fellenz et al., 2018; 2020; Crawford et al., 2021). Neither of these two prior phases of SPE identified natural earthquake-related surface effects in the same region that occurred at a similar depth as the explosive experiments for direct comparison. While earthquakes can also produce surface changes that are detectable using remote sensing data analyses, it is expected that the pattern and spatial extent of surface changes would vary between earthquakes and explosions. However, no direct-observed surface-change signature comparison between earthquakes and explosions has ever been performed. The SPE Phase 3 Rock Valley/Direct Comparison (RV/DC) program presents a unique opportunity to investigate and characterize co-occurrence of both earthquakes and explosions. In this project, we worked to address five tasks in a workflow, as follows: 1. Identify and obtain existing high-resolution legacy satellite and aerial imagery as close in time before and after the 1993 Rock Valley earthquake sequence to temporally constrain the analyses. 2. Transform these pre-earthquake and post-earthquake datasets into digital elevation models (DEMs) using geospatial analysis software packages (e.g., ArcGIS, Agisoft Metashape, and Google Earth Engine). 3. Perform DEM differencing analyses to assess and quantify earthquake-related changes from the 1993 sequence, and develop map products that visualize these analyses. 4. Use the analyses from (3) to: (a) assess spatial distribution and magnitude of surface changes due to the 1993 Rock Valley earthquake sequence, and (b) determine parameters of forthcoming, planned explosion-related surface change data collection from sensors mounted on unmanned aerial vehicles (UAVs) (e.g., spatial extent of collection, design and density of survey control, sensors to deploy, forward speed and line spacing of UAV flight lines, flight altitude). 5. Develop a summary report on the analyses, including how the analyses define parameters and identify focus areas for any future surface change analytical field campaigns related to the explosive experiment. Analyzing these legacy data and identifying whether they can detect any surface changes related to the earthquake sequence facilitates opportunities for direct signature comparison of surface change from explosions at one location, which has never previously been performed. Comparing the surface change signatures from a co-located and depth-equivalent earthquake and an explosion could help to advance remote sensing event discrimination techniques. This report summarizes the work completed toward this ambitious goal.

42 ENGINEERING↗

Microstructural Simulation of Ion-Irradiated Natural Rocks and Minerals

This report details the contribution of Oak Ridge National Laboratory to the Nuclear Energy University Program (NEUP) project, "Rapid Characterization of Concrete Mineralogy using Multi-Scale Tools" led by the University of Illinois at Urbana-Champaign. To support the subsequent license renewal (SLR) of the US nuclear power plants (NPPs) fleet, the performance of concrete-forming aggregates against neutron irradiation needs to be assessed. Ion irradiation is proposed as a rapid, cost-effective surrogate method to full neutron irradiation testing. To complement the characterization of ion-irradiated rock specimens, ORNL has run finite-elements and fast-Fourier transform (FFT-based) simulations using the codes MARS and Microstructure-Oriented Scientific Analysis of Irradiated Concrete (MOSAIC). The main conclusion of this analysis is that the apparent post–ion-irradiation step-height underestimates the accumulated free radiation-induced volumetric expansion (RIVE) in the ion-implanted depth by about 15% at full amorphization and about 25% toward the beginning of the ion irradiation experiment. The discrepancy is explained by the fact that the step height is proportional to the sum of the RIVE (i.e., one third of the RIVE) and the irradiation-assisted plastic strains in the vertical direction (lower than two thirds of the RIVE). Because of the large lateral compressive stresses caused by the restraining effect of the unirradiated substrate, the stress field in the mineral grains and at the grain boundary (GB) in the ion-implanted layer is much different than the comparable stress field occurring during neutron irradiation. Hence, the mismatch RIVE causing cracks in the rock-forming minerals leads to different cracking patterns. Ion irradiation continues to be used as a rapid technique to assess the RIVE potential of rocks.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Rock Valley Dense Gravity Acquisition

Characterizing the shallow structure of the Rock Valley region of the Nevada National Security Site is a critical component of the Rock Valley Direct Comparison project. Geophysical data of the region is needed for operational decisions, to constrain geologic models used for simulation, and to facilitate the analysis of future explosive source data. Local measurements of gravity are a key piece of geophysical information that helps to resolve the underlying geologic composition, fault structure, and density characteristics, yet, in the Rock Valley region these measurements are sparse on the scale of the testbed. In this report, we present the details of a recent gravity data acquisition survey designed to collect a dense dataset in the region of interest that complements the existing gravity work but greatly enhances our resolution. This dataset will be integrated with a complementary Los Alamos National Laboratory gravity collection and combined with the existing seismic data in a joint inversion. These measurements were conducted over two weeks with a portable gravimeter and high-resolution GPS and include repeat measurements at a USGS base station as well as reoccupation of gravity sites in the regional dataset. This collection of over 100 new dense gravity measurements will facilitate refinement of the existing Geologic Framework Model and directly complement newly acquired dense seismic data, ultimately improving the project’s ability to investigate the direct comparison of shallow earthquake and explosive sources.

58 GEOSCIENCES↗

Life Cycle Analysis of Enhanced Rock Weathering

This study aims to evaluate the life cycle impacts of weathering agricultural fields in the U.S. Midwest with igneous rock (dunite/olivine) and industrial waste minerals using literature and estimated data from a techno-economic analysis (TEA). The GWP impact with igneous rock and industrial waste materials ranges 34–263 and 23–180 kg CO2-eq/tonne CO2 sequestered, respectively, with comminution and transportation as the main environmental hotspots for igneous rock, and transportation for industrial waste.

Izar-Tenorio, Jorge↗

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

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

25 ENERGY STORAGE↗

RockFlow: Fast Generation of Synthetic Source Rock Images Using Generative Flow Models

Image-based evaluation methods are a valuable tool for source rock characterization. The time and resources needed to obtain images has spurred development of machine-learning generative models to create synthetic images of pore structure and rock fabric from limited image data. While generative models have shown success, existing methods for generating 3D volumes from 2D training images are restricted to binary images and grayscale volume generation requires 3D training data. Shale characterization relies on 2D imaging techniques such as scanning electron microscopy (SEM), and grayscale values carry important information about porosity, kerogen content, and mineral composition of the shale. Here, we introduce RockFlow, a method based on generative flow models that creates grayscale volumes from 2D training data. We apply RockFlow to baseline binary micro-CT image volumes and compare performance to a previously proposed model. We also show the extension of our model to 2D grayscale data by generating grayscale image volumes from 2D SEM and dual modality nanoscale shale images. The results show that our method underestimates the porosity and surface area on the binary baseline datasets but is able to generate realistic grayscale image volumes for shales. With improved binary data preprocessing, we believe that our model is capable of generating synthetic porous media volumes for a very broad class of rocks from shale to carbonates to sandstone.

58 GEOSCIENCES↗