Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture characterization”

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 1,225 records · Page 68

Nondestructive Evaluation of Stress Corrosion Cracking in a Welded Steel Plate Using Guided Ultrasonic Waves

Stress corrosion cracking (SCC) had occurred in early-generation high-level nuclear waste tanks constructed by welding carbon steel. This paper describes an ultrasonic inspection system and its fundamental ability to detect and quantify the length of SCC on thick welded steel plates. The finite element method (FEM) was applied to simulate the welding process to estimate the welding residual stress field. Growth of stress corrosion cracks is driven by crack stress intensities exceeding the subcritical cracking threshold intensity. The subject plate was experimentally inspected with ultrasonic nondestructive evaluation (NDE) techniques to characterize the extent of SCC. The NDE system uses a piezoelectric transducer to generate guided waves in the thick steel plate, and a scanning laser Doppler vibrometer (SLDV) to measure multidimensional time–space wavefield data over a user-defined scanning area in the plate surface. The measured wavefield data can show wave interactions in a localized area in the plate due to the presence of the discontinuities of the SCC. To generate an inspection image that can precisely show the crack’s location and/or the dimension, the wavefield data are further processed to generate inspection image that maps the entire sample plate so the crack can be clearly identified in the plate while its length can be readily estimated. In conclusion, the ultrasonic test results for crack length agree well with the visually estimated length and are close to that predicted by the FEM for cracks in the weld residual stress field.

36 MATERIALS SCIENCE↗

Fast strength model characterization using Bayesian statistics

A variety of flow stress models exist with new models constantly being developed. These models aim to approximate the strength of materials in a variety of regimes from quasistatic loading through shock scenarios. All models contain an array of parameters which need to be tuned to the material under study. Some models perform well under limited conditions, requiring adjustment of the parameters when venturing outside of those predefined ranges. Other models perform well over a wide range of conditions with a set of parameters, but may be outperformed by other models optimized on a tighter range of conditions. Recent research by Los Alamos demonstrated the ability to optimize the Johnson Cook (JC) model using a set of 3 plate-impact experiments on Aluminum. They utilized Bayesian statistics and emulation to determine optimal parameters for the model with a quantification of parameter uncertainty. We present an extension of this capability to incorporate velocimetry from plate-impact tests, stress-strain data from split Hopkinson pressure bar and quasistatic compression tests, plus profiles from Taylor cylinders in a unified fashion. Statistically robust comparisons of the performance and uncertainty of different realizations of the JC flow stress model were carried out based on calibration to several possible combinations of these three different experiment types.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗

Reactive transport modeling of the Aquifer Thermal Energy Storage (ATES) system at Stockton University, New Jersey during seasonal operations

Hydrogeochemical processes associated with Aquifer Thermal Energy Storage (ATES) operations can often impact the system performance owing to mineral precipitation either at the wellbore or in the aquifer owing to changes in temperature and fluid disequilibria. Although failure of ATES systems due to mineral precipitation ("fouling") is common, predictive reactive-transport models have rarely been applied to plan their design and operation. Here, the objective of this study is to develop a reactive-transport model by coupling thermal, hydrological, and chemical (THC) processes to evaluate effects of introduced atmospheric oxygen on water chemistry, mineral precipitation/dissolution, porosity, and permeability changes associated with an ATES system at Stockton University (New Jersey, USA). The THC model builds on a Thermal-Hydrological-Mechanical (THM) model of the site that evaluated system failure owing to possible fracturing in the caprock or around the wellbore. The causes of the system failure are not known – potential causes include hydraulic fracturing owing to elevated pump pressures that took place, a flow pathway created by one of the boreholes, or a pre-existing natural hydrologic connection between the upper unconfined aquifer and the ATES aquifer, any of which could have led to oxygenated water entering the reservoir and causing the observed Fe-oxide fouling on well screens. The THC model is used to evaluate some of the hypotheses and observations regarding system failure owing to geochemical processes. The reactive-transport code TOUGHREACT V4 was used to model the THC processes during seasonal heating and cooling operations at the Stockton ATES site over 6 years of operation. In the THC simulations, the primary effects on geochemistry were observed when the injection water is saturated with atmospheric oxygen. Simulations show greater precipitation of goethite near the cold wells as compared to the warm wells. Although volume fractions of Fe-hydroxides were relatively small, the model was aimed at processes in the aquifer at the scale of meters and larger rather than at the scale of mm or cm (i.e., a well screen). Kaolinite is the dominant precipitating phase, also around the cold wells. Illite dissolves near the cold wells and precipitates near the warm wells. There is a net decrease in the porosity near the cold wells and increase near the warm wells, although a slight amount of thermal contraction near the cold wells and expansion near the warm wells is responsible for a significant proportion of the porosity change. Owing to the coarse discretization of the numerical grid near the wells (compared to the screen thickness) the magnitude of permeability changes at the wellbore are likely underestimated. The reactive transport model in this study can be used for characterization of aquifers, optimizing the operational parameters (temperature, pressure, pH etc.), and planning of mitigation strategies for ATES systems.

15 GEOTHERMAL ENERGY↗

DEEPEN: Final 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

Part of the DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments) project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology, detailed in the journal article below, is based on the method proposed by Poux & O'brien (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., (2022). For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty.

15 GEOTHERMAL ENERGY↗

First-results from the Perseverance SHERLOC Investigation: Aqueous Alteration Processes and Implications for Organic Geochemistry in Jezero Crater, Mars

The Perseverance rover landed in Jezero crater, a site selected to fulfill the Mars-2020 mission goals of characterizing the geology of habitable environments and searching for signs of life while collecting samples for return to Earth [1]. Jezero hosted an open-basin lake during the late Noachian/early Hesperian (~3.7 Ga) [1-2], has units associated with the largest carbonate deposit identified on Mars [3-4], and has a well-preserved delta with clay and carbonate-bearing sediments, well-suited to preservation of organics [1,3-4]. Investigating the nature of organics and aqueous environments within their geologic con-text allows us to understand important aqueous processes and determine habitability within Jezero crater. Previous in situ landed measurements of organics could not resolve their spatial and mineralogical con-text [5-6]. Although Martian meteorites lack geological context, the spatial distribution of organic compounds in Martian meteorites have allowed recognition of an association between aqueous processes and organics [7-8]. Here, we show for the first time in-situ associations between carbonate-forming ultramafic alteration process, later stage aqueous sulfate and perchlorate formation, and organics on the Martian surface. Methodology and geological context: We use the Perseverance rover’s SHERLOC instrument (Scanning Habitable Environments with Raman and Lumines-ence of Organics and Chemicals), a deep-ultraviolet fluorescence and Raman scattering spectrometer capable of mapping the organic and mineral composition with a spatial resolution of 100 μm resolution to report the presence of organics and aqueously formed minerals at Jezero crater [9]. These spectral detections were compared with co-located images obtained with the autofocus context imager (ACI) and the WATSON camera for textural analysis [9]. As of writing, the Per-severance rover has abraded five targets that were measured with the SHERLOC instrument. The five targets are located in two different orbitally-identified geological units within the floor of Jezero crater; the Crater Floor Fractured Rough unit (CF-Fr) and the Séítah region within the Crater Floor Fractured 1 unit (CF-F1) [10]. In orbital infrared spectroscopic data, the CF-Fr unit is associated with pyroxene spectral signatures and minor alteration, while the Séítah region is associated with olivine and minor Mg-rich carbonates and clays [3-4,10]. Carbonation of ultramafic protolith recorded within Jezero crater: All scans of abraded targets within the Séítah region reveal strong peaks at 1080–1090 cm−1 consistent with carbonate and peak singlets or doublets at 820–840 cm−1 attributed to olivine (Fig. 1), consistent with orbital infrared observations. Our detailed micron-scale petrographic and spectroscopic evidence shows that these carbonates formed through carbonation of an ultramafic protolith. The supporting observations include: (1) Carbonate cation compositions match those of olivine, suggesting mixed Fe- and Mg-olivine gave rise to mixed Fe- and Mg-carbonates, similar to observations of ultramafic systems on Earth and within Martian meteorites [3-4,7-8]. (2) The ob-served carbonates co-occur with hydrated materials, gypsum, and potentially aqueously-formed phases, amorphous silicates and phosphate. (3) The spectral and textural variation of olivine and carbonate dominated zones and olivine-carbonate mixtures within both primary grains and interstitial zones are expected for carbonated ultramafic protoliths. (4) These mineral associations and textures closely resemble those observed within the ALH84001 and Nakhlite meteorites attributed to olivine carbonation on Mars [6-7]. Taken together, micron-scale SHERLOC documentation of these phenomena bridge previous orbital and meteorite observations and demonstrate in-situ regionally extensive (~106 km2) ultramafic alteration resulting in geo-logical deposition of carbonates. Furthermore, we observe that olivine carbonation was involved in preserving and possibly synthesizing organics, which makes this environment potentially habitable, as previously suggested in [1,3-4] (Fig. 1). Late-stage aqueous perchlorate and sulfate in Jezero crater: An abrasion target within the CF-Fr unit contains combinations of high intensity 950-955 cm−1 peaks and minor 1090-1095 cm−1 and 1150-1155 cm−1 peaks that are spectral fits to anhydrous perchlorate (Fig. 1). Some spectra show a combination of 950-955 cm−1 peaks with equally strong 1010-1020 cm−1 peaks, low intensity broad features at 1120 cm−1, and occasional broad 3450 cm−1 hydration (-OH) features, indicating a mixture of Ca-sulfate and perchlorate that is minimally hydrated (Fig. 1). The detections of per-chlorates within Jezero crater differ from previous measurements (e.g. Phoenix lander, Curiosity rover, Tissint meteorite [7,11]) because they are observed to be intimately related to aqueous processes including sulfate formation, they present as a secondary white void-fill occurring within the interior of the rock, and they are found to likely be Na-perchlorate. implications for their formation: Three different types of organics embedded within three different lithologies were observed within the abraded targets. Organics associated with low intensity ~340 nm fluorescence were widespread within targets with no apparent association to particular minerals (Fig. 1). Organics associated with ~305 nm and ~275 nm fluorescence correlated with sulfates within the Bellegarde target in the CF-Fr unit, while organics associated with high intensity ~340 nm fluorescence correlated with carbonate, phosphate, and amorphous silicate mixtures within the Garde target in the Séítah region (Fig. 1). Although assignment of fluorescence signatures to specific organic compounds is not conclusive, ~340 nm fluorescence is generally more consistent with 2-ring aromatic organics, ~275 nm fluorescence is more consistent with 1-ring aromatic organics, and ~305 nm fluorescence can be created by either 2- or 1-ring aromatics [12]. These observations indicate that the strongest fluorescence signatures interpreted as organics were found in materials associated with aqueous processes, i.e. sulfate- and carbonate-bearing materials, suggesting both brines and ultramafic carbonation aqueous environments were capable of preserving organics on ancient Mars. In Martian meteorites, simple aromatic organics proposed to have been synthesized through aqueous processes can be found within minerals associated with olivine carbonation and in spatial association with perchlorate and sulfate materials [7-8], similar to SHERLOC observations. Hence, we advance an abiotic aqueous synthesis origin for the organics although we cannot rule out the presence of organics from meteoritic in-fall or putative organic biosignatures. Detailed analyses will be required upon return of these materials to Earth.

E L Scheller↗

Simulated meteorite impacts and volcanic explosions: Ejecta analyses and planetary implications

Past cratering studies have focused primarily on crater morphology. However, important questions remain about the nature of crater deposits. Phenomena that need to be studied include the distribution of shock effects in crater deposits and crater walls; the origin of mono- and polymict breccia; differences between local and distal ejecta; deformation induced by explosive volcanism; and the production of unshocked, high-speed ejecta that could form the lunar and martian meteorites found on the Earth. To study these phenomena, one must characterize ejecta and crater wall materials from impacts produced under controlled conditions. New efforts at LLNL simulate impacts and volcanism and study resultant deformation. All experiments use the two-stage light-gas gun facility at LLNL to accelerate projectiles to velocities of 0.2 to 4.3 km/s, including shock pressures of 0.9 to 50 GPa. We use granite targets and novel experimental geometries to unravel cratering processes in crystalline rocks. We have thus far conducted three types of simulations: soft recovery of ejecta, 'frozen crater' experiments, and an 'artificial volcano. Our ejecta recovery experiments produced a useful separation of impactites. Material originally below the projectile remained trapped there, embedded in the soft metal of the flyer plate. In contrast, material directly adjacent to the projectile was jetted away from the impact, producing an ejecta cone that was trapped in the foam recovery fixture. We find that a significant component of crater ejecta shows no signs of strong shock; this material comes from the near-surface 'interference zone' surrounding the impact site. This phenomenon explains the existence of unshocked meteorites on the Earth of lunar and martian origin. Impact of a large bolide on neighboring planets will produce high-speed, weakly shocked ejecta, which may be trapped by the Earth's gravitational field. 'Frozen crater' experiments show that the interference zone is highly localized; indeed, disaggregation does not extend beyond approx. 1.5 crater radii. A cone-shaped region extending downward from the impact site is completely disaggregated, including powdered rock that escaped into the projectile tube. Petrographic analysis of crater ejecta and wall material will be presented. Finally, study of ejecta from 0.9- and 1.3-GPa simulations of volcanic explosions reveal a complete lack of shock metamorphism. The ejecta shows no evidence of PDF's, amorphization, high-pressure phases, or mosaicism. Instead, all deformation was brittle, with fractures irregular (not planar) and most intergranular. The extent of fracturing was remarkable, with the entire sample reduced to fragments of gravel size and smaller.

Gratz, A. J.↗

Ultrasonic/Sonic Rotary-Hammer Drills

Ultrasonic/sonic rotary-hammer drill (USRoHD) is a recent addition to the collection of apparatuses based on ultrasonic/sonic drill corer (USDC). As described below, the USRoHD has several features, not present in a basic USDC, that increase efficiency and provide some redundancy against partial failure. USDCs and related apparatuses were conceived for boring into, and/or acquiring samples of, rock or other hard, brittle materials of geological interest. They have been described in numerous previous NASA Tech Briefs articles. To recapitulate: A USDC can be characterized as a lightweight, lowpower, piezoelectrically driven jackhammer in which ultrasonic and sonic vibrations are generated and coupled to a tool bit. A basic USDC includes a piezoelectric stack, an ultrasonic transducer horn connected to the stack, a free mass ( free in the sense that it can bounce axially a short distance between hard stops on the horn and the bit), and a tool bit. The piezoelectric stack creates ultrasonic vibrations that are mechanically amplified by the horn. The bouncing of the free mass between the hard stops generates the sonic vibrations. The combination of ultrasonic and sonic vibrations gives rise to a hammering action (and a resulting chiseling action at the tip of the tool bit) that is more effective for drilling than is the microhammering action of ultrasonic vibrations alone. The hammering and chiseling actions are so effective that unlike in conventional twist drilling, little applied axial force is needed to make the apparatus advance into the material of interest. There are numerous potential applications for USDCs and related apparatuses in geological exploration on Earth and on remote planets. In early USDC experiments, it was observed that accumulation of cuttings in a drilled hole causes the rate of penetration of the USDC to decrease steeply with depth, and that the rate of penetration can be increased by removing the cuttings. The USRoHD concept provides for removal of cuttings in the same manner as that of a twist drill: An USRoHD includes a USDC and a motor with gearhead (see figure). The USDC provides the bit hammering and the motor provides the bit rotation. Like a twist drill bit, the shank of the tool bit of the USRoHD is fluted. As in the operation of a twist drill, the rotation of the fluted drill bit removes cuttings from the drilled hole. The USRoHD tool bit is tipped with a replaceable crown having cutting teeth on its front surface. The teeth are shaped to promote fracturing of the rock face through a combination of hammering and rotation of the tool bit. Helical channels on the outer cylindrical surface of the crown serve as a continuation of the fluted surface of the shank, helping to remove cuttings. In the event of a failure of the USDC, the USRoHD can continue to operate with reduced efficiency as a twist drill. Similarly, in the event of a failure of the gearmotor, the USRoHD can continue to operate with reduced efficiency as a USDC.

Badescu, Mircea↗

Microscopic and Macroscopic Characterization of Grain Boundary Energy and Strength in Silicon Carbide via Machine-Learning Techniques

Predicting the properties of grain boundaries poses a challenge because of the complex relationships between structural and chemical attributes both at the atomic and continuum scales. Grain boundary systems are typically characterized by parameters used to classify local atomic arrangements in order to extract features such as grain boundary energy or grain boundary strength. The present work utilizes a combination of high-throughput atomistic simulations, macroscopic and microscopic descriptors, and machine-learning techniques to characterize the energy and strength of silicon carbide grain boundaries. Additionally, a diverse data set of symmetric tilt and twist grain boundaries are described using macroscopic metrics such as misorientation, the alignment of critical low-index planes, and the Schmid factor, but also in terms of microscopic metrics, by quantifying the local atomic structure and chemistry at the interface. These descriptors are used to create random-forest regression models, allowing for their relative importance to the grain boundary energy and decohesion stress to be better understood. Results show that while the energetics of the grain boundary were best described using the microscopic descriptors, the ability of the macroscopic descriptors to reasonably predict grain boundaries with low energy suggests a link between the crystallographic orientation and the resultant atomic structure that forms at the grain boundary within this regime. For grain boundary strength, neither microscopic nor macroscopic descriptors were able to fully capture the response individually. However, when both descriptor sets were utilized, the decohesion stress of the grain boundary could be accurately predicted. These results highlight the importance of considering both macroscopic and microscopic factors when constructing constitutive models for grain boundary systems, which has significant implications for both understanding the fundamental mechanisms at work and the ability to bridge length scales.

36 MATERIALS SCIENCE↗

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗

Compositional and Microstructural Evolution of Olivine During Pulsed Laser Irradiation: Insights Based on a FIB/Field-Emission TEM Study

Introduction: The use of pulsed laser irradiation to simulate the short duration, high-energy conditions characteristic of micrometeorite impacts is now an established approach in experimental space weathering studies. The laser generates both melt and vapor deposits that contain nanophase metallic Fe (npFe(sup 0)) grains with size distributions and optical properties similar to those in natural impact-generated melt and vapor deposits. There remains uncertainty, however, about how well lasers simulate the mechanical work and internal (thermal) energy partitioning that occurs in actual impacts. We are currently engaged in making a direct comparison between the products of laser irradiation and experimental/natural hypervelocity impacts. An initial step reported here is to use analytical TEM is to attain a better understanding of how the microstructure and composition of laser deposits evolve over multiple cycles of pulsed laser irradiation. Experimental Methods: We irradiated pressed-powder pellets of San Carlos olivine (Fo(sub 90)) with up to 99 rastered pulses of a GAM ArF excimer laser. The irradiated surface of the sample were characterized by SEM imaging and areas were selected for FIB cross sectioning for TEM study using an FEI Quanta dual-beam electron/focused ion beam instrument. FIB sections were characterized using a JEOL2500SE analytical field-emission scanning transmission electron microscope (FE-STEM) optimized for quantitative element mapping at less than 10 nm spatial resolutions. Results: In the SEM the 99 pulse pressed pellet sample shows a complex, inhomogeneous, distribution of laser-generated material, largely concentrated in narrow gaps and larger depressions between grains. Local concentrations of npFe0 spherules 0.1 to 1 micrometers in size are visible within these deposits in SEM back-scatter images. Fig. 1 shows bright-field STEM images of a FIB cross-section of a one of these deposits that continuously covers the top and sloping side of an olivine grain. The deposit has 3 microstructurally distinct sub-layers composed of silicate glass with varying modal fractions and size distributions of npFe( sup 0) spherules, along with nanocrystalline silicate material. A relatively thin (50-300 nm) topmost surface layer has a high-concentration of npFe0 spherules 5-20 nm in size. Element mapping shows the layer to be enriched in Fe by a factor of 2.5 relative to the olivine substrate, with Mg and Si depleted by 20% and 10% respectively. This is compositionally complementary to the underlying, middle layer of the deposit that is depleted in Fe, enriched in Mg and has a much lower npFe0 concentration. A third layer of nanocrystalline olivine occurs at the substrate interface. Discussion: The FE-STEM results suggest the topmost layer is a vapor deposit, underlain by a thicker microstructurally complex melt-generated layer. The compositional relations suggest the melt layer was partially vaporized, preferentially losing more volatile elements (e.g., Fe). The vaporized material re-condensed to form the thin, npFe(sup 0)-rich surface deposit during or immediately after the scan cycle. Nanocrystalline olivine that grew within the melt layer as it formed and cooled is similar in volume and microstructure to what we have observed in the impact melt lining of a micrometeorite impact crater in olivine. This suggest the time-temperature relations attained in the laser sample may not be too different from a micrometeorite impact. Our TEM observations, however, do not show evidence for the same level of mechanical dam-age (e.g., fracturing) seen around the natural micrometeorite crater.

Christoffersen, R.↗

Deterministic and Statistical Modeling of a New Thermal Breakout Technology for Measuring the Maximum Horizontal In-Situ Stress

The current state-of-the-art technology for in-situ stress measurements involves an integrated approach that combines borehole breakout observations, drilling-induced tensile fractures, and hydraulic fracturing tests (i.e., “mini-fracs”). This methodology has achieved wide application in the oil and gas industry but has several limitations that often prevent successful in-situ stress measurements. One major limitation is that breakouts do not appear in all boreholes and are generally only a natural occurrence that cannot easily be controlled. Because borehole breakouts are used to directly measure the maximum horizontal in-situ stress magnitude, the absence of borehole breakouts presents a major data gap for in-situ stress measurements. In response to this data gap, a new US Department of Energy (US DOE)-sponsored thermal breakout technology that will provide a method for thermally inducing borehole breakouts and allow the consistent measurement of the maximum horizontal stress magnitude is currently in development. This thermal breakout technology involves heating the borehole and increasing the thermoelastic compressive stress in the rock until a breakout develops, which can be directly correlated to the maximum horizontal stress magnitude. The first step in this project was an analytical modeling study of the thermal breakout process. Based on the Kirsch solution (Kirsch 1898), a deterministic and statistical analysis was performed on the pertinent parameters that influence the maximum horizontal stress calculation. As a result of the analysis, the findings indicate that the thermal breakout technology is feasible and provides improved accuracy and/or an enhanced ability to measure the maximum horizontal stress. Future work as part of this US DOE-sponsored project includes additional validation through more detailed numerical modeling, laboratory testing, and field testing of the thermal breakout technology.

58 GEOSCIENCES↗

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗

AERoBOND Project Summary

Under NASA’s Convergent Aeronautics Solutions (CAS) project, the Adhesive-Free Bonding of Complex Composites (AERoBOND) project investigated off-stoichiometric epoxy polymers for fast, reliable assembly of epoxy matrix composite structures. The project goal was to demonstrate feasibility of the AERoBOND joining method by demonstrating mechanical properties greater than 80% of conventional co-cured materials while reducing structure weight by 1%. The project consisted of three convergent research areas: material and process development, systems analysis, and material and process modeling. Material and process development was the largest component of AERoBOND with approximately 6 FTE and 1WYE of support to formulate and characterize new resins, prepare carbon fiber prepregs, fabricate laminates, measure mechanical properties, analyze failure results, and select material and process improvements. The systems analysis activity estimated the potential reduction in part count and aircraft weight by comparing models of composite wing boxes with no fasteners (co-cured structure), fasteners in major joints (co-cured stringers), and fasteners in all joints. The materials and process modeling activity included a molecular model of the AERoBOND materials system to predict mechanical properties of resins with offset stoichiometry and a process model to predict the effect of resin formulation and processing conditions on the extent of mixing and degree of cure in a finished joint. As the number of airline passenger trips doubles in the next 20 years (IATA/Tourism Economics Air Passenger Forecasts, April 2019), the increased demand for new commercial aircraft is now the single greatest technical challenge to the airframe manufacturing industry. To meet efficiency requirements, new aircraft must be fabricated primarily from high performance structural composites, but manufacturing processes are inherently slow with the largest bottleneck attributed to assembly and installation of fasteners (NASA/TM–2019-220428). Manufactures of commercial transport aircraft are compelled to install more than 100,000 redundant fasteners into bonded joints to prevent failures due to unpredictable weak bonds. In structural adhesive bonds, the interface between adherend and adhesive is nearly two-dimensional making it susceptible to minute quantities of contamination, which can cause weak bonds. Currently, bond strength assessment is only possible through destructive testing (i.e., breaking the joint). For these reasons, regulatory organizations such as the Federal Aviation Administration (FAA) often require redundant load paths in secondary-bonded, primary-structures to alleviate concerns with bond performance. The AERoBOND process enables reflow of matrix resin during assembly to eliminate the material discontinuity at the interface, thereby eliminating the dependence of mechanical performance on interfacial adhesion. The AERoBOND joint is equivalent to the interlaminar region obtained during a co-cure process, so joint performance depends on the cohesive properties of the matrix resin. Conventional co-cured structures, although too costly and complex for large-scale manufacturing, are trusted by manufacturers and regulators, and are certified for flight with few or no redundant fasteners.Systems analysis performed on a composite wing model at the scale of a single-aisle commercial transport aircraft indicated that >20,000 redundant fasteners per wing could be eliminated by implementing the AERoBOND joining method. A total weight reduction of 15% was predicted in a wing box by eliminating fasteners and thinning components that must no longer support localized fastener loads and accommodate fastener dimensions. Interlaminar shear fracture toughness measured by the end-notched flexure test was greater than 1 kJ/m2 (nearly 140% of the co-cured benchmark property), which is greatly in excess of the project goals for mechanical properties. Testing was planned to measure interlaminar tensile fracture toughness as well as interlaminar tensile and shear strengths using the same AERoBOND configuration, but was delayed due to closure of LaRC facilities during the COVID-19 pandemic. The AERoBOND process model is partially validated and available for experimental use. It allows the user to input AERoBOND process parameters such as material composition, laminate configuration, and cure cycle to predict the final cure state of the AERoBOND joint. A preliminary, multi-scale material model was developed to predict AERoBOND joint mechanical properties (stiffness and strength) based on the cure state of the joint provided by the process model. The timing for transition of this technology within NASA is excellent as NASA initiates new enduring projects to address composites manufacturing rate challenges. AERoBOND technology is well suited to AAVP/AATT objectives for rapid manufacturing of a composite wing. A minimal effort (1 FTE/$15k procurement/0 WYE) is proposed in FY21 to continue a minor mechanical testing effort and maintain a SAA with ASX composites to develop commercial quality prepreg material. An RFI with the composites industry is suggested to quantify the technology gap between the current TRL and the TRL needed for transition to industry. A moderate effort [3-4 FTE/$150k/1 WYE (~$115k)] is proposed in FY22 for the “high rate composites manufacturing” project currently in planning. The partnership with ASX Composites will be expanded to produce material for sub-element/element-scale “panel-off” activities. Industry partnerships with airframe manufacturers is an expected component to explore damage tolerance and environmental stability. Further development of multi-scale modeling tools (process model, meso-scale model, and molecular model) is planned to enhance and deliver tools for rapid manufacturing infusion.

Frank Louis Palmieri↗

Mineralogy and Geochemistry of Sedimentary Rocks and Eolian Sediments in Gale Crater, Mars: A Review After Six Earth Years of Exploration with Curiosity

The Mars Science Laboratory Curiosity rover arrived at Mars in August 2012 with a primary goal of characterizing the habitability of ancient and modern environments. Curiosity was sent to Gale crater to study a sequence of ∼3.5 Ga old sedimentary rocks that, based on orbital visible and near- to short-wave infrared reflectance spectra, contain secondary minerals that suggest deposition and/or alteration in liquid water. The sedimentary sequence in the lower slopes of Mount Sharp in Gale crater preserves a dramatic shift on early Mars from a relatively warm and wet climate to a cold and dry climate, based on a transition from smectite-bearing strata to sulfate-bearing strata. The rover is equipped with instruments to examine the sedimentology and identify compositional changes in the stratigraphy. The Chemistry and Mineralogy (CheMin) instrument is one of two internal laboratories on Curiosity and includes a transmission X-ray diffractometer (XRD) and X-ray fluorescence (XRF) spectrometer. CheMin measures loose sediment samples scooped from the surface and drilled rock powders, and the XRD provides quantitative mineralogy to a detection limit of ∼1 wt.% for crystalline phases. Curiosity has traversed>20 km since landing and has primarily been exploring an ancient lake environment fed by streams and groundwater. Of the 19 drilled rock samples analyzed by CheMin as of sol 2300 (January 2019), 15 are from fluvio-lacustrine deposits that comprise the Bradbury and Murray formations. Most of these samples were drilled from units that did not have a clear mineralogical signature from orbit. Results from CheMin demonstrate an astounding diversity in the mineralogy of these rocks that signifies geochemical variations in source rocks, transportation mechanisms, and depositional and diagenetic fluids. Most detrital igneous minerals are basaltic, but the discovery in a few samples of abundant silicate minerals that usually crystallize from evolved magmas on Earth remains enigmatic. Trioctahedral smectite and magnetite at the base of the section may have formed from low-salinity pore waters with a circumneutral pH in lake sediments. A transition to dioctahedral smectite, hematite, and Ca-sulfate going up section suggests a change to more saline and oxidative aqueous conditions in the lake waters themselves and/or in diagenetic fluids. Perhaps one of the biggest mysteries revealed by CheMin is the high abundance of X-ray amorphous materials (15–73 wt.%) in all samples drilled or scooped to date. CheMin has analyzed three modern eolian sands, which have helped constrain sediment transport and mineral segregation across the active Bagnold Dune Field. Ancient eolian sandstones drilled from the Stimson formation differ from modern eolian sands in that they contain abundant magnetite but no olivine, suggesting that diagenetic processes led to the alteration of olivine to release Fe(II) and precipitate magnetite. Fracture-associated halos in the Stimson and the Murray formations are evidence for complex aqueous processes long after the streams and lakes vanished from Gale crater. The sedimentology and composition of the rocks analyzed by Curiosity demonstrate that habitable environments persisted intermittently on the surface or in the subsurface of Gale crater for perhaps more than a billion years.

E.B. Rampe↗

Deep Learning–Assisted Multiobjective Optimization of Geological CO 2 Storage Performance under Geomechanical Risks

In geological CO 2 storage, designing the optimal well control strategy for CO 2 injection to maximize CO 2 storage while minimizing the associated geomechanical risks is not trivial. This challenge arises due to pressure buildup, CO 2 plume migration, the highly nonlinear nature of geomechanical responses to rock-fluid interaction, and the high computational cost associated with coupled flow and geomechanics simulations. In this paper, we introduce a novel optimization framework to address these challenges. The optimization problem is formulated as follows: maximize total CO 2 storage while minimizing geomechanical risks by adjusting the injection schedules within bounded constraints. The geomechanical risks are primarily driven by injection-induced pressure build-up, which is characterized by ground displacement and the induced microseismicity. We used the Fourier neural operator (FNO)-based deep learning model to construct surrogate models, replacing the time-consuming coupled flow and geomechanics simulations for evaluating the aforementioned objective functions. The developed surrogate models have been incorporated into a multiobjective optimization framework through a genetic algorithm to reduce the computational burden. The proposed optimization framework reduces the computational cost from approximately 2,400 hours, when using objective function evaluations based on physics-based simulations, to around 20 minutes. A set of Pareto-optimal solutions of the proposed workflow yields nontrivial optimal decisions, reducing the microseismicity potential and the vertical displacement. This Pareto front highlights the optimal trade-offs between CO 2 storage amount, safety, and ground displacement, emphasizing the need for careful optimization and management of injection strategies to achieve a balanced outcome. The novelty of this work is twofold. First, we demonstrate the importance of incorporating the minimization of the geomechanical risks as objective functions into the CO 2 storage optimization workflow to mitigate the potential risk of induced microseismicity and ground displacement. Second, we leverage the FNO-based surrogate models to optimize a real-field CO 2 storage operation.

42 ENGINEERING↗

Biogeochemistry of the Antrim Shale Natural Gas Reservoir

The Antrim Shale, located in the Michigan Basin, United States (U.S.), is a major U.S. shale play having produced over 2.5 Trillion Cubic Feet (Tcf) of unconventional shale natural gas as of 2010. The shallow nature of this formation sets it apart from other, more characterized unconventional shale gas plays. The depth of gas production of the Antrim ranges from approximately 150 to 600 m and it is typically vertically drilled, contrary to deeper, horizontally drilled shales. A thorough understanding of the biogeochemistry and microbiology of this complex system will be advantageous for improving well performance, produced water management, and potential biocidal treatment as microbial community composition can vary substantially even among closely spaced wells. In this study, we analyzed produced water collected from nine different wells in the Antrim Shale by investigating the geochemical and microbial community composition of the produced water to gain greater insight into the overall biogeochemistry of this unique shale system. The majority of the wells from this study had high total dissolved solids (TDS) primarily composed of chloride and sodium, averaging 86,804 mg/L with a maximum 116,223 mg/L; however, three of the wells sampled along the northern margin of the basin exhibited significantly lower TDS ranging from 4932 to 6496 mg/L. Overall, our microbial community analysis revealed relatively low abundance within our samples and high variability of the microbial community among the sampled wells. The majority of bacterial sequences were identified within Proteobacteria, Firmicutes, and Actinobacteria phyla and metagenomic sequencing revealed the low presence of Methanobacteriaceae within each sample. We also investigated potential microbial community drivers and found that TDS, sodium, chloride, iodide, bromide, ammonium, potassium, and strontium were significantly correlated with the observed microbial community. The varying geochemical conditions between wells demonstrate different subsurface environmental niches, potentially driving the heterogeneous microbial communities we observed from well to well. This analysis suggests an important relationship between both well location and geochemistry and the observed microbial community that can persist in the reservoir. Continued studies of the Antrim Shale will improve our understanding of the complex interdependencies of this ecosystem.

03 NATURAL GAS↗