SEARCH · Search NASA

Results for “stress plot”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRS2000

Design and Evaluation of a New Boundary-Layer Rake for Flight Testing

A new boundary-layer rake has been designed and built for flight testing on the NASA Dryden Flight Research Center F-15B/Flight Test Fixture. A feature unique to this rake is its curved body, which allows pitot tubes to be more densely clustered in the near-wall region than conventional rakes allow. This curved rake design has a complex three-dimensional shape that requires innovative solid-modeling and machining techniques. Finite-element stress analysis of the new design shows high factors of safety. The rake has passed a ground test in which random vibration measuring 12 g rms was applied for 20 min in each of the three normal directions. Aerodynamic evaluation of the rake has been conducted in the NASA Glenn Research Center 8 x 6 Supersonic Wind Tunnel at Mach 0-2. The pitot pressures from the new rake agree with conventional rake data over the range of Mach numbers tested. The boundary-layer profiles computed from the rake data have been shown to have the standard logarithmic-law profile. Skin friction values computed from the rake data using the Clauser plot method agree with the Preston tube results and the van Driest II compressible skin friction correlation to approximately +/-5 percent.

Bui, Trong T.↗

DOE OSTI2019

Analysis of strengthening due to grain boundaries and annealing twin boundaries in the CrCoNi medium-entropy alloy

CrCoNi exhibits the best combination of strength and ductility among all the equiatomic single-phase FCC subsets of the CrMnFeCoNi high-entropy alloy. Here, its yield strength was determined in compression as a function of grain size and temperature. Yield strength was also plotted as a function of “crystallite” size, which takes into account both annealing twin boundaries and grain boundaries. The resulting Hall-Petch slopes were straight lines but with different slopes that depend on the number of twin boundaries per grain. Scanning transmission electron microscopy of deformed specimens revealed the formation of dislocation pile-ups at grain and annealing twin boundaries indicating that the latter also act as obstacles to slip and contribute to strength. Using a simple pile-up model, the strengths of the grain and twin boundaries were estimated to lie in the range 900–1250 MPa. Assuming that they have the same strength, in the case of twin boundaries this strength corresponds roughly to the stress required to constrict Shockley partials, which suggests that dissociated dislocations have to become compact before they can cross the annealing twin boundaries.

36 MATERIALS SCIENCE↗

NASA NTRS2015

Insulation Resistance Degradation in Ni-BaTiO3 Multilayer Ceramic Capacitors

Insulation resistance (IR) degradation in NiBaTiO3 multilayer ceramic capacitors has been characterized by the measurement of both time to failure (TTF) and direct current leakage current as a function of stress time under highly accelerated life test conditions. The measured leakage current time dependence data fit well to an exponential form, and a characteristic growth time tau (sub SD) can be determined. A greater value of tau (sub SD) represents a slower IR degradation process. Oxygen vacancy migration and localization at the grain boundary region results in the reduction of the Schottky barrier height and has been found to be the main reason for IR degradation in NiBaTiO3 capacitors. The reduction of barrier height as a function oftime follows an exponential relation of phi (t ) = phi (0) e (exp -2Kt), where 13 the degradation rate constant K Koe (Ek/kT) is inversely proportional to the mean TTF (MTTF) and can be determined using an Arrhenius plot. For oxygen vacancy electromigration, a lower barrier height phi (0) will favor a slow IR degradation process, but a lower phi (0) will also promote electronic carrier conduction across the barrier and decrease the IR. As a result, a moderate barrier height phi (0) (and therefore a moderate IR value) with a longer MTTF (smaller degradation rate constant K) will result in a minimized IR degradation process and the most improved reliability in NiBaTiO3 multilayer ceramic capacitors.

reliability↗

DOE OSTI2026

Multi‐season analysis reveals hundreds of drought‐responsive genes in sorghum

Persistent drought affects global crop production and is becoming more severe in many parts of the world in recent decades. Deciphering how plants respond to drought will facilitate the development of flexible mitigation strategies. Sorghum bicolor L. Moench (sorghum), a major cereal crop and an emerging bioenergy crop, exhibits remarkable resilience to drought. To better understand the molecular traits that underlie sorghum's remarkable drought tolerance, we undertook a large-scale sorghum gene expression profiling effort, totaling nearly 1500 transcriptome profiles, across a 3-year field study with replicated plots in California's Central Valley. This study included time-resolved gene expression data from roots and leaves of two sorghum genotypes, BTx642 and RTx430, with different pre-flowering and post-flowering drought-tolerance adaptations under control and drought conditions. Quantification of genotype-specific drought tolerance effects was enabled by de novo sequencing, assembly, and annotation of both BTx642 and RTx430 genomes. These reference-quality genomes were used to construct a pangene set for characterizing conserved and genotype-specific expression. By integrating time-resolved transcriptomic responses to drought in the field across three consecutive years, we identified a set of 726 drought-responsive genes that responded similarly in all 3 years of our field study. Functional enrichment analysis identified abiotic stress, secondary cell wall-related processes and metabolism as particularly affected under both types of drought stress. We also found that some glyoxylate cycle pathway genes, including malate synthase and isocitrate lyase, are differentially regulated particularly during post-flowering drought stress, implicating this pathway as potentially important for drought responsiveness. This expansive dataset represents a unique resource for sorghum and drought research communities and provides a methodological framework for the integration of multi-faceted time-resolved transcriptomic datasets.

Cole, Benjamin [USDOE Joint Genome Institute (JGI)↗

NASA NTRS2015

Insulation Resistance Degradation in Ni-BaTiO3 Multilayer Ceramic Capacitors

Insulation resistance (IR) degradation in Ni-BaTiO3 multilayer ceramic capacitors has been characterized by the measurement of both time to failure and direct-current (DC) leakage current as a function of stress time under highly accelerated life test conditions. The measured leakage current-time dependence data fit well to an exponential form, and a characteristic growth time SD can be determined. A greater value of tau(sub SD) represents a slower IR degradation process. Oxygen vacancy migration and localization at the grain boundary region results in the reduction of the Schottky barrier height and has been found to be the main reason for IR degradation in Ni-BaTiO3 capacitors. The reduction of barrier height as a function of time follows an exponential relation of phi (𝑡)=phi (0)e(exp -2Κt), where the degradation rate constant 𝐾=𝐾o𝑒(𝐸𝑘/𝑘𝑇) is inversely proportional to the mean time to failure (MTTF) and can be determined using an Arrhenius plot. For oxygen vacancy electromigration, a lower barrier height phi(0) will favor a slow IR degradation process, but a lower phi(0) will also promote electronic carrier conduction across the barrier and decrease the insulation resistance. As a result, a moderate barrier height phi(0) (and therefore a moderate IR value) with a longer MTTF (smaller degradation rate constant 𝐾) will result in a minimized IR degradation process and the most improved reliability in Ni-BaTiO3 multilayer ceramic capacitors.

dielectric degradation↗

DOE OSTI2023

Characterizing the flux effect on the irradiation embrittlement of reactor pressure vessel steels using machine learning

In-service exposure to high-energy neutrons embrittles reactor pressure vessel (RPV) steels. An increase in the yield stress (Δσ y ) results in a corresponding increase in the brittle to ductile transition temperature (ΔT c ). Most existing models underpredict ΔT c at higher fluence following accelerated irradiations in test reactors. High fluence, up to 10 20 n/cm 2 in some cases, will be reached over extended RPV vessel operation of 80 years, or more, at low service flux. Embrittlement has been extensively studied in accelerated, higher flux test reaction irradiations. However, the use of test reactor data naturally raises the question of flux effects. Here, this study used a machine learning approach trained on a set of hardening data, covering a wide range of flux, fluence, and steel compositions to determine the interactive effects of both irradiation and material variables on Δσ y . The analysis included machine learning-based cross-plots of the variable dependence of Δσ y for six core steels (i.e., CM6, LC, LD, LG, LH, and LI), with controlled differences in their Cu and Ni contents. A primary objective is to evaluate an effective to actual fluence (Φt e /Φt) ratio, as a function of flux, fluence, and steel composition. This is information critical to properly use intermediate flux-high fluence data in calibrating a low flux-high fluence embrittlement model. The predicted Φt e /Φt is reasonably consistent with estimates previously derived from a physics-based solute recombination trap model.

36 MATERIALS SCIENCE↗

DOE OSTI2025

A two-loop four-point form factor at function level

Recently, the maximally-helicity-violating four-point form factor for the chiral stress-energy tensor in planar $\mathcal{N}$ = 4 super Yang-Mills was computed to three loops at the level of the symbol associated with multiple polylogarithms. It exhibits antipodal self-duality, or invariance under the combined action of a kinematic map and reversing the ordering of letters in the symbol. Here we lift the two-loop form factor from symbol level to function level. We provide an iterated representation of the function’s derivatives (coproducts). In order to do so, we find a three-parameter limit of the five-parameter phase space where the symbol’s letters are all rational. We also use function-level information about dihedral symmetries and the soft, collinear, and factorization limits, as well as limits governed by the form-factor operator product expansion (FFOPE). We provide plots of the remainder function on several kinematic slices, and show that the result is compatible with the FFOPE data. We further verify that antipodal self-duality is valid at two loops beyond the level of the symbol.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Enhancing the Quality and Reliability of Machine Learning Interatomic Potentials through Better Reporting Practices

Recent developments in machine learning interatomic potentials (MLIPs) have empowered even nonexperts in machine learning to train MLIPs for accelerating materials simulations. However, reproducibility and independent evaluation of presented MLIP results is hindered by a lack of clear standards in current literature. In this Perspective, we aim to provide guidance on best practices for documenting MLIP use while walking the reader through the development and deployment of MLIPs including hardware and software requirements, generating training data, training models, validating predictions, and MLIP inference. We also suggest useful plotting practices and analyses to validate and boost confidence in the deployed models. Finally, we provide a step-by-step checklist for practitioners to use directly before publication to standardize the information to be reported. Altogether, we hope that our work will encourage the reliable and reproducible use of these MLIPs, which will accelerate their ability to make a positive impact in various disciplines including materials science, chemistry, and biology, among others.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2000

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗

NASA NTRS2004

Diffraction Studies of the Atomic Vibrations of Bulk and Surface Atoms in the Reciprocal and Real Spaces of Nanocrystalline SiC

To describe and evaluate the vibrational properties of nanoparticles it is necessary to distinguish between the surface and the core of the particles. Theoretical calculations show that vibrational density of states of the inner atoms of nanograins is similar to bulk material but shifted to higher energies which can be explained by the fact that the gain core is stressed (hardened) due to the presence of internal pressure. Theoretical calculations also show that there is a difference between vibrational properties of a crystal lattice of the grain interior in isolated particles and in a dense (sintered) nanocrystalline material. This is probably due to a coupling of the modes inside the grains via the grain boundaries in dense nanocrystalline bodies. We examined strains present in the surface shell based on examination of diamond and Sic nanocrystals in reciprocal (Bragg-type scattering) and real (PDF analysis) space analysis of neutron diffraction data. Recently we examined the atomic thermal motions in nanocrystalline Sic based on the assumption of a simple Einstein model for uncorrelated atomic notions. According to this model, the Bragg intensity is attenuated as a function of scattering angle by the Debye-Waller factor. Based on this assumption overall temperature factors were determined from the Wilson plots.

Stelmakh, S.↗

DOE OSTI2021

Root lateral interactions drive water uptake patterns under water limitation

The plasticity of root water uptake determines the maintenance of transpiration during periods of water limitation and drought. However, the mechanistic basis of plant water uptake, as well as the implications of water uptake strategies at the individual and ecosystem scale remain elusive. We model three-dimensional root water uptake under variably saturated conditions for a one-hectare temperate forest plot for a growing season with a pronounced mid-season dry period. Variations in root architecture, hydraulic properties, and degree of lateral interaction between root systems produce divergent local responses to water limitation and provide insights on individual and community response to meteorological conditions. Results demonstrate the plasticity of ecosystem-scale responses to surface drying, where interacting roots shift regions of active uptake to deeper soil layers with less abundant root biomass. These shifts, a product of both root system and soil hydraulic properties, illustrate intimate links between root and soil hydraulics in determining plant water sourcing. Furthermore, we demonstrate that root lateral interactions are beneficial at the ecosystem-scale, even when trees compete for water. Specifically, a more spatially extensive root system facilitates access to a larger soil water reservoir, often ameliorating water limitation and reducing sharp water potential gradients. While the reduction of water stress is a benefit, it can be offset by increased root construction and maintenance costs associated with the larger rooting system. A plausible “viability” region of root communal co-existence is therefore implied where competitive pressures and root production costs are balanced by a potential water benefit.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Analytical Model to Detect Fault Permeability Alteration Induced by Fault Reactivation in Compartmentalized Reservoirs

Underground fluid injection induces changes in in situ stress condition of the target formation and local faults that can potentially lead to fault reactivation, which may result in the leakage of injected and/or native fluids into neighboring formations. In this paper, we introduce an analytical method to detect fault reactivation caused by fluid injection into deep faulted aquifers considering across-fault leakage. The fundamental assumption made in our model is that fault permeability will be altered upon fault slip. Therefore, we model fault reactivation as a sudden change in fault permeability at the onset of fault slip. The fault is modeled as a linear interface between two permeable formations with equal rock and fluid properties. The governing equations are coupled through the fault interface and are solved using the Laplace-Fourier integral transform technique. Based on the analytical solution, we find the characteristic bottomhole pressure and pressure derivative responses that enable detecting fault reactivation using diagnostic plots. We observe that pressure derivative undergoes a rapid change at the onset of fault slip followed by a late-time trend to attain a new equilibrium governed by the altered fault permeability. Furthermore, we discuss the evolution of the across-fault leakage rate upon and after fault slippage. Here, the results from this study are presented in the form diagnostic plots and type curves that may be used for reservoir and fault characterization purposes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

AGR-3/4 TRISO Fuel Compact Ceramography

The combined third and fourth irradiation in the Advanced Gas Reactor (AGR) program (AGR-3/4) contained tristructural isotropic (TRISO)-coated particle fuel and designed-to-fail (DTF) fuel particles. The DTF particles were only coated with highly anisotropic pyrocarbon (PyC) so they would purposely fail during the AGR-3/4 irradiation and provide a source of fission products for measurement. To observe the post-irradiation morphology of these DTF particles and the TRISO-coated particles, three AGR-3/4 compacts were mounted in epoxy, sectioned above their centerlines, and ground/polished. Three rounds of grinding/polishing and optical microscopy were performed so that the particles could be observed at multiple planes. Each compact contained approximately 1,872 TRISO particles and exactly 20 DTF particles. The three AGR-3/4 compacts examined covered a wide temperature range and featured both the hottest average irradiation temperature (1375°C) and the coldest/lowest burnup (872°C and 5.5% fissions per initial metal atom) of any compacts to undergo post-irradiation examination (PIE) in the AGR program to date. A total of 29 DTF particles were located and observed via microscopy across the three compacts. All the DTF particles observed via microscopy had completely failed PyC coatings. Some different DTF kernel and PyC morphologies were observed that have been attributed to the differences in irradiation temperature and/or burnup. At low and medium irradiation temperatures (roughly 850 to 1050°C), it appears irradiation-induced dimensional changes in the DTF PyC caused it to fracture and fold in on itself, and the kernel deformed to accommodate this PyC deformation. The extent of DTF kernel deformation to accommodate DTF PyC buckling tended to be greater for the medium-temperature fuel compared to the cold fuel. In the high-irradiation-temperature compact (1375°C), the DTF PyC layer appears to have completely reacted chemically with the kernel material. The only material surrounding the DTF particles in this hot compact resembles the compact matrix material. Observations of the TRISO-coated driver fuel particles and fuel compact graphitic matrix were also made. The TRISO particle kernels in the medium-temperature compact (1047°C), and many in the high-temperature compact (1375°C), showed morphologies consistent with what was seen in AGR-1 and AGR-2. In other TRISO particles from the high-temperature compact, a spatial gradient in the kernel appearance was evident. Here the cool side of the kernel (away from the center of the compact) had a much darker appearance (like that of the buffer and pyrocarbon). This is the first observation of this kind of kernel spatial gradient in AGR fuel. It is believed that the high irradiation temperature of that compact (an average of 1375°C) activates or accelerates an unidentified chemical reaction, and the temperature gradient within the fuel causes a sharp spatial gradient. The low-temperature/low-burnup compact also showed unique kernel morphologies where the TRISO-coated particles had clearly distinguishable oxide and carbide phases at the center of the kernel, an oxide rind surrounding this, and remnants of a carbide skin surrounding the oxide rind. These are features observed in the as-fabricated fuel that are no longer present in higher burnup fuel. Occasional gaps between the outer pyrolytic carbon (OPyC) and graphitic matrix material were observed in all compacts. Gaps were most often found in the small, matrix-filled spaces between adjacent particles. Finally, buffer morphologies in AGR-3/4 TRISO particles were like those observed in AGR-1 and AGR-2, and the buffer fracture frequencies in AGR-3/4 and AGR-2 TRISO particles were plotted versus fast neutron fluence and irradiation temperature. In AGR-3/4 the buffer fracture rate was highest (23%) for an irradiation temperature of 1047°C and a fast fluence of 5.18E25 n/m 2 . Increasing the irradiation temperature to 1375°C reduced the buffer fracture rate to 14%. Reducing the irradiation temperature to 872°C and the fast fluence by a factor of 3 reduced the buffer fracture rate to 7%. The temperature and fluence dependencies observed for AGR-3/4 buffer fracture were consistent with those observed in AGR-2. Lower fluences and/or higher temperatures significantly reduced the buffer fracture frequencies. Lower buffer fracture frequencies were also observed at lower temperatures as long as the fluence was also significantly reduced. This is believed to be due to lower fluence resulting in less irradiation-induced dimensional change (shrinkage) of the buffer, and a higher temperature promoting enhanced creep relaxation of stresses within the buffer, leading to less buffer fracture.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2024

Usage of ChatGPT for Engineering Design and Analysis Tool Development

ChatGPT, a generative AI large language model, has recently captured significant attention in both the computer science community and the broader public domain. It has demonstrated a wide range of capabilities, from answering simple questions to writing fully functional computer code. This study spotlights both the capabilities and limitations of ChatGPT when addressing engineering problems. The model's capacity to generate practical engineering tools is highlighted through an example of a prompt that leads to an interactive plotting tool, enabling the examination of the fluid boundary layer around a fan blade. Subsequently, the paper also uncovers potential pitfalls in ChatGPT’s application, shown through an unsuccessful attempt to use ChatGPT to automate a process in Ansys Workbench through scripting. The research further investigates ChatGPT's proficiency in addressing inquiries and providing explanations about the functionalities of OpenMDAO, an open-source, multidisciplinary design, analysis, and optimization tool developed at NASA Glenn Research Center. Finally, an optimization methodology, developed with ChatGPT’s help, is applied to the structural optimization of a fan blade. The developed optimization method utilizes T-Blade3 for geometry generation, Ansys Mechanical for meshing and finite element analysis, and sci-kit learn’s MLPRegressor method to generate a trained neural network model of the design space. OpenMDAO is then used to find the optimal point within the design space. The outcome is a significant reduction in stress in the optimized model—less than one-fifth of the stress value in the baseline model.

Design↗

NASA NTRSDate not supplied

Usage of ChatGPT for Engineering Design and Analysis Tool Development

Design↗

DOE OSTI2022

Disturbance legacies regulate coastal forest soil stability to changing salinity and inundation: A soil transplant experiment

Coastal forests worldwide are vulnerable to a dramatic transition from upland to wetland as sea-level rise accelerates and regimes of precipitation and storms change. However, the biogeochemical impacts of shifting salinity and inundation disturbance that foreshadow forest to wetland state transitions are largely unknown. This experiment used a natural salinity gradient in a tidal creek in eastern Maryland, U.S.A., to examine how soil respiration and chemistry may change under novel salinity and inundation disturbance regimes. Soil monoliths were transplanted in a reciprocal design among plots varying in seawater exposure and elevation above the creek. We monitored the monoliths’ carbon dioxide (CO 2 ) flux for two years and performed soil chemical analyses at the end of the experiment. Soil CO 2 flux was affected by changing disturbance regimes and responses were dependent upon the salinity and inundation legacies associated with each study location. Lowland soil CO 2 flux was resistant to changing salinity and inundation disturbance, with transplanted soil monolith chemistry composition in between that of its origin and destination. Conversely, upland soil CO 2 flux was sensitive to changing salinity and inundation disturbance and remained suppressed throughout the 2-year study when exposed to wetter, saline conditions. Additionally, transplanted upland soil chemistry composition was like that of its destination, rather than origin, with upland soils displaying higher pH, base saturation, and nutrient availability relative to lowland soils. Here, we hypothesize that reductions in soil respiration rates were driven by loss of soil nutrients and osmotic and redox stress on microbial communities following exposure to seawater. Together, our results suggest that disturbance legacies shape coastal forest soil responses to changing salinity and inundation disturbance regimes. However, fully understanding the dependence of system responses on disturbance legacies requires future study across a variety of systems and spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

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↗

DOE OSTI2022

Modeling the topographic influence on aboveground biomass using a coupled model of hillslope hydrology and ecosystem dynamics

Abstract. Topographic heterogeneity and lateral subsurface flow at the hillslope scale of ≤1 km may have outsized impacts on tropical forest through their impacts on water available to plants under water-stressed conditions. However, vegetation dynamics and finer-scale hydrologic processes are not concurrently represented in Earth system models. In this study, we integrate the Energy Exascale Earth System Model (E3SM) land model (ELM) that includes the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), with a three-dimensional hydrology model (ParFlow) to explicitly resolve hillslope topography and subsurface flow and perform numerical experiments to understand how hillslope-scale hydrologic processes modulate vegetation along water availability gradients at Barro Colorado Island (BCI), Panama. Our simulations show that groundwater table depth (WTD) can play a large role in governing aboveground biomass (AGB) when drought-induced tree mortality is triggered by hydraulic failure. Analyzing the simulations using random forest (RF) models, we find that the domain-wide simulated AGB and WTD can be well predicted by static topographic attributes, including surface elevation, slope, and convexity, and adding soil moisture or groundwater table depth as predictors further improves the RF models. Different model representations of mortality due to hydraulic failure can change the dominant topographic driver for the simulated AGB. Contrary to the simulations, the observed AGB in the well-drained 50 ha forest census plot within BCI cannot be well predicted by the RF models using topographic attributes and observed soil moisture as predictors, suggesting other factors such as nutrient status may have a larger influence on the observed AGB. The new coupled model may be useful for understanding the diverse impact of local heterogeneity by isolating the water availability and nutrient availability from the other external and internal factors in ecosystem modeling.

54 ENVIRONMENTAL SCIENCES↗