Search NASA⌕ Search

SEARCH · Search NASA

Results for “Contrast Sensitivity”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Pedestal collapse by resonant magnetic perturbations

Pedestal collapse (i.e., the complete loss of the edge transport barrier (ETB)) in DIII-D H-mode plasmas occurs when resonant magnetic perturbations (RMPs) penetrate the steep gradient region at the plasma edge. Normally, RMP driven magnetic islands can occur at the top and bottom of the H-mode pedestal and these islands generate conditions consistent with edge-localized-mode (ELM) suppression and density pump-out, respectively, based on nonlinear two-fluid MHD simulations. In contrast, MHD simulations show that the steep pressure gradient region between the top and bottom of the DIII-D pedestal is generally immune to resonant field penetration due to large local E × B and diamagnetic flows. By this fortuitous circumstance, the edge-transport-barrier and H-mode confinement can be maintained while achieving ELM suppression. However, pedestal collapse can occur in DIII-D when the screening flows are inadequate to prevent field penetration in the steep gradient region of the pedestal. Non-linear two-fluid MHD simulations support the role of resonant field penetration in pedestal collapse for DIII-D H-mode plasmas with weak edge E × B and diamagnetic screening flows. ITER will likely have weaker edge screening flows than present experiments due to its much larger size, making it more susceptible to resonant field penetration in the steep gradient region of the pedestal. Analysis of model ITER equilibria demonstrates that resonant field penetration in the steep pressure gradient region is possible for RMP levels of the order required for ELM suppression. The effect of such penetration on the ITER pedestal will depend sensitively on the resulting degree of island overlap.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Extreme Risk Mitigation in Reinforcement Learning using Extreme Value Theory

Risk-sensitive reinforcement learning (RL) has garnered significant attention in recent years due to the growing interest in deploying RL agents in real-world scenarios. A critical aspect of risk awareness involves modelling highly rare risk events (rewards) that could potentially lead to catastrophic outcomes. These infrequent occurrences present a formidable challenge for data-driven methods aiming to capture such risky events accurately. While risk-aware RL techniques do exist, they suffer from high variance estimation due to the inherent data scarcity. Our work proposes to enhance the resilience of RL agents when faced with very rare and risky events by focusing on refining the predictions of the extreme values predicted by the state-action value distribution. To achieve this, we formulate the extreme values of the state-action value function distribution as parameterized distributions, drawing inspiration from the principles of extreme value theory (EVT). We propose an extreme value theory based actor-critic approach, namely, Extreme Valued Actor-Critic (EVAC) which effectively addresses the issue of infrequent occurrence by leveraging EVT-based parameterization. Importantly, we theoretically demonstrate the advantages of employing these parameterized distributions in contrast to other risk-averse algorithms. Our evaluations show that the proposed method outperforms other risk averse RL algorithms on a diverse range of benchmark tasks, each encompassing distinct risk scenarios.

Wang, Yu↗
DOE OSTI2020

Northstar Mo100 Target Beam Spot Size Sensitivity Study

Northstar Mo100 target for production of the medical isotope Mo99 is particularly sensitive to beam over focusing because of the intended high power density. At 38 MeV and 6.32 mA (240 kW beam power, nominally 160 kW th in the target), small deviation from the nominal 12 mm FWHM beam can significantly increase the target disk temperatures and more importantly the load bearing target window temperature. A rhodotron will be the source of the electron beam to be used in the isotope production. This machine operates in a pulsed beam mode, currently expected to be at 50 Hz, which results in a 100 C temperature oscillation above and below the steady state temperature that is the focus of this study. This must be considered in the determination of minimum beam spot size that can be accepted by the target without risk of failure. In addition to the nominal beam spot size of 12 mm, the analysis of 10, 11 and 14.5 mm beam spot sizes are analyzed. The primary conclusion is that an over focused beam to just 11 mm FWHM is the minimum beam spot allowable, based on window temperature and the Ultimate Tensile Stress (UTS) of the window material Inconel 718 as a function of temperature. In contrast an under focused beam of 14.5 mm FWHM has much reduced target and window temperatures, but there is a resulting 6% drop in isotope production, all other factor staying the same. These results are presented herein.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Development of an Ultrahigh-bandwidth Phase Contrast Imaging System for detection of electron scale turbulence and Gigahertz Radio-Frequency Waves

The study of waves and turbulence is vital to the development of future reactor-grade plasma devices developed in the quest for fusion energy. These fluctuations are responsible for moving heat and particles across the magnetic field, and a predictive understanding of them is needed to achieve the density and temperature required to sustain a plasma fusion reaction. While many techniques have been developed for measuring waves and fluctuations, every measurement method has limitations. There are relatively few techniques for measuring very high frequency fluctuations, such as radio frequency waves injected to heat the plasma, unstable waves driven by suprathermal particles, or short wavelength electrostatic waves driven by electron temperature or density gradients. The present project builds upon the proven phase contrast imaging (PCI) technique to extend the response of the diagnostic by orders of magnitude in frequency and almost a factor of ten in spatial resolution. PCI provides a measurement of electron density based on small angle scattering of a CO-2 laser beam by using optical techniques to render a phase shift as an intensity change on a detector. Due to the telecommunications revolution, technological development by manufacturers has focused on components in the near -infrared, so that high- quality lasers and detectors at 1.55 µm are readily available. Shifting PCI design to a new, shorter wavelength has numerous advantages and challenges. Similar detector performance is available with room-temperature arrays with GHz bandwidth, while the detector arrays for 10.6 µm required liquid nitrogen cooling and were therefore limited to a bandwidth of about 1 MHz. Shifting to a shorter wavelength reduces the angle at which the laser beam scatters off of plasma waves, allowing more such scattered components to pass through the aperture of the vacuum vessel port and be collected by the PCI, which increases the spatial resolution. Concomitant with these benefits, various questions of performance arise. At shorter laser wave- length, the sensitivity to mirror and lens quality is increased, the contribution of the laser to the overall system noise is increased, and the sensitivity to vibrations is increased. The custom optical components at the heart of the PCI technique were required to be properly scaled for the shorter wavelength, so fabrication technologies needed to be explored. This project was designed to show that a low noise, high response PCI system at 1.55 µm could be constructed and operated, and then to quantify potential issues to allow extrapolation to a full-size production system providing physics measurements on a large plasma device. The first stage, producing the custom optical component called a Phase Plate, was successfully achieved using two methods. First, an easily reproduced masking and coating technique was able to produce good phase plates with the required parameters. Second, a nanofabrication technique was found to produce extremely high quality phase plates at a competitive price. The PCI constructed with the new phase plates and 1.55 µm laser was found to provide excellent wavelength measurements with the theoretically expected response. The sensitivity to optical surface quality was found to be in line with previous measurements at 10.6 µm. The observed signal-to-noise ratio was similar to the theoretically expected value. The effect of vibrations on PCI was studied with the first measurement of the effect of beam motion on PCI response and comparison to theory, allowing for a quantitative prediction of the effect of vibrations on a production PCI system and the requirements for improved beam stabilization. PCI is an extremely cost-effective method to provide a low noise, absolutely calibrated measurement of plasma fluctuations across a wide spatial scale. This project has shown that a 1.55 µm PCI using modern techniques and components is less expensive than the 10.6 µm PCI of a few years ago, with the largest savings in phase plate fabrication and the infrared detector array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Influence of sea-ice anomalies on Antarctic precipitation using source attribution in the Community Earth System Model

Abstract. We conduct sensitivity experiments using a general circulation model that has an explicit water source tagging capability forced by prescribed composites of pre-industrial sea-ice concentrations (SICs) and corresponding sea surface temperatures (SSTs) to understand the impact of sea-ice anomalies on regional evaporation, moisture transport and source–receptor relationships for Antarctic precipitation in the absence of anthropogenic forcing. Surface sensible heat fluxes, evaporation and column-integrated water vapor are larger over Southern Ocean (SO) areas with lower SICs. Changes in Antarctic precipitation and its source attribution with SICs have a strong spatial variability. Among the tagged source regions, the Southern Ocean (south of 50° S) contributes the most (40 %) to the Antarctic total precipitation, followed by more northerly ocean basins, most notably the South Pacific Ocean (27%), southern Indian Ocean (16 %) and South Atlantic Ocean (11 %). Comparing two experiments prescribed with high and low pre-industrial SICs, respectively, the annual mean Antarctic precipitation is about 150 Gt yr-1 (or 6 %) more in the lower SIC case than in the higher SIC case. This difference is larger than the model-simulated interannual variability in Antarctic precipitation (99 Gt yr-1). The contrast in contribution from the Southern Ocean, 102 Gt yr-1, is even more significant compared to the interannual variability of 35 Gt yr-1 in Antarctic precipitation that originates from the Southern Ocean. The horizontal transport pathways from individual vapor source regions to Antarctica are largely determined by large-scale atmospheric circulation patterns. Vapor from lower-latitude source regions takes elevated pathways to Antarctica. In contrast, vapor from the Southern Ocean moves southward within the lower troposphere to the Antarctic continent along moist isentropes that are largely shaped by local ambient conditions and coastal topography. This study also highlights the importance of atmospheric dynamics in affecting the thermodynamic impact of sea-ice anomalies associated with natural variability on Antarctic precipitation. Our analyses of the seasonal contrast in changes of basin-scale evaporation, moisture flux and precipitation suggest that the impact of SIC anomalies on regional Antarctic precipitation depends on dynamic changes that arise from SIC–SST perturbations along with internal variability. The latter appears to have a more significant effect on the moisture transport in austral winter than in summer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Simulated meteorological impacts of offshore wind turbines and sensitivity to the amount of added turbulence kinetic energy

Offshore wind energy projects are currently in development off the east coast of the United States and may influence the local meteorology of the region. Wind power production and other commercial uses in this area are related to atmospheric conditions, and so it is important to understand how future wind plants may change the local meteorology. In the absence of measurements of potential wind plant impacts on meteorology, simulations offer the next-best possible insight into wake effects on boundary layer height, temperature, fluxes, and wind speeds. However, simulation tools that capture these effects offer multiple options for representing the amount of turbine-added turbulence that may impact assessments of micrometeorological effects. To explore this sensitivity, we compare 1 year of simulations from the Weather Research and Forecasting (WRF) model with and without wind plants incorporated, focusing on the lease area south of Massachusetts and Rhode Island. The simulations with wind plants are repeated to include both the maximum and minimum amounts of added turbulence to provide bounds on the potential impacts. We assess changes in wind speeds, 2 m temperature, surface heat flux, turbulence kinetic energy (TKE), and boundary layer height during different stability classifications and ambient wind speeds over the entire year and compare results for the degree of added turbulence in the wind plant simulations. Because the wake behavior may be a function of boundary layer stability, in this paper, we also present a machine learning algorithm to quantify the area and distance of the wake generated by the wind plant. This analysis enables us to identify the relationship between wake extent and boundary layer height. We find that hub-height wind speed is reduced within and downwind of the wind plant, with the strongest impacts occurring during stable conditions and faster wind speeds in region 3 of the turbine power curve, although impacts lessen as wind speeds increase past 15 m s−1. In contrast, wind speeds near the surface decrease when no turbine-added turbulence is included but can increase for stably stratified conditions when 100 % of possible TKE is included in the simulations. TKE increases at hub height in the simulations with added TKE for all stability classes, suggesting that atmospheric stability does not immediately modify the TKE generated by turbines. Negligible changes in hub-height TKE manifest in the simulations without the added TKE. At the surface, TKE increases in the simulations with maximum added turbulence only for unstable conditions. In the no-added-turbulence simulations, surface TKE decreases slightly in neutral and unstable simulations. Differences in 2 m temperatures and surface heat fluxes are small but vary considerably with atmospheric stability and the amount of added TKE. Boundary layer heights increase within the wind plant when turbine-added turbulence is included and decrease slightly downwind during stable conditions. In contrast, with no added turbulence, the boundary layer height is in general reduced in stable conditions with wind speeds less than 15 m s −1 and slightly increased in neutral conditions. Finally, shallower upwind boundary layer heights tend to correlate with larger wake areas and distances, though other factors likely also play a role in determining the extent of the wind plant wake. These simulation-based results provide a bound for micrometeorological impacts of wind plant wakes: simulations that couple the atmosphere to the ocean may reduce these impacts, and we await observational verification.

17 WIND ENERGY↗
DOE OSTI2022

Influence of gasoline fuel formulation on lean autoignition in a mixed-mode-combustion (deflagration/autoignition) engine

We report stoichiometric spark-ignition engines suffer efficiency penalties due to throttling losses at low loads, a low specific-heat ratio of the stoichiometric working fluid, and limits on compression ratio due to end-gas autoignition leading to undesirable knocking. Mixed-Mode Combustion (MMC) mitigates these shortcomings by using a lean working fluid where a spark-initiated pilot-stabilized deflagrative flame front is followed by controlled end-gas autoignition. This MMC study investigates the effects of initial conditions (intake air temperature, intake pressure, equivalence ratio, and intake oxygen fraction) on autoignition tendency of four gasoline-range fuels with varying properties and composition. The use of fuels with varying octane sensitivity (S) allowed exploring the importance of low-temperature heat release in triggering autoignition. Fuels with high S were less reactive for conditions that promote low-temperature chemistry (operation at high intake air pressure or without N2 dilution). Conversely, an Alkylate fuel with low S showed a greater autoignition resistance at operating conditions that were unfavorable for low-temperature chemistry. Next, the effect of residual gas composition on autoignition tendency of fuels was examined with a chemical-kinetics model. Among the various molecules in the residual gas, nitric oxide (NO) enhanced the low-temperature chemistry and increased the autoignition tendency most significantly. The fuels’ autoignition response to increasing NO amount corroborates the experimental observations. Next, the sequential autoignition of the end-gas was assessed to be less impacted by thermal stratification because of lean mixtures showing relatively less low-temperature chemistry, when compared to stoichiometric mixtures. Next, the effect of changing equivalence ratio on the autoignition was found to be similar for all fuels, regardless of their S. With changing intake air temperature, the response of fuels’ autoignition tendency depended on the dilution level used. At high dilution (i.e. low intake [O2]), fuels’ reactivity increased with increasing intake air temperature. In contrast, for operation without dilution, the autoignition tendency of the low-S Alkylate fuel decreased with increasing intake air temperature, while that of high-S High Cycloalkane fuel still increased with increasing intake air temperature. In conclusion, conventional octane metrics (RON and MON) have utility in assessing the autoignition tendency under lean MMC operation. Moreover, the fuel requirements for MMC align with that of stoichiometric operation: i.e., high RON and high S fuels are desirable for stable non-knocking operation.

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

Single and dual metal atom catalysts for enhanced singlet oxygen generation and oxygen reduction reaction

We demonstrate rational design of graphene-supported single and dual metal atom catalysts (SACs and DACs) for photocatalytic applications, such as singlet oxygen ( 1 O 2 ) sensitization and H 2 O 2 production. Here we combine density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations with experimental verifications. We found a synergistic effect between triplet sensitization and triplet–triplet (Dexter) energy transfer; both play a role in the photocatalytic activity through the volcano plot of 3d transition metal SACs. More specifically, FeN 4 -SAC exhibits a low ISC energy gap (ΔE ISC ) of 0.039 eV, compared with 0.108 eV for FeNiN8-DACs, both possessing a high Bader charge transfer of 0.366 e- and 0.405 e-, respectively. Guided by these computational results, we synthesized a series of SACs and a DAC and confirmed their structures with scanning transmission electron microscopy (STEM) along with the X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS). We then confirm their band structures with low-energy inverse photoemission spectroscopy (LEIPS) and UV-vis-NIR. Subsequently, we synthesized the catalysts for the photooxygenation of anthracene and two-electron oxygen reduction reaction (ORR) to measure their photocatalytic activity. We found that H 2 O 2 production through the two-electron ORR competes with the 1 O 2 generation through Dexter energy transfer. FeN 4 -SAC demonstrates a high photooxygenation conversion of 86% and a high 1 O 2 quantum yield of 1.04, obtained from electron spin resonance (ESR) spectroscopy, with low H 2 O 2 production. In contrast, NiN 4 -SAC exhibits a low 1 O 2 generation and a high H 2 O 2 production mainly because of the high Gibbs free energy of the OOH* intermediate. In conclusion, this work proposes an effective DFT-guided strategy for designing SACs and DACs for various photocatalytic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Characterization of oceanic post-cold frontal clouds and their model representation. Final report

Low-level clouds are ubiquitous over the oceans, and their impact on the Earth’ radiative budget is important. However, general circulation models (GCMs) that are used to represent our climate and its evolution still experience problems to represent their correct amount and radiative impact. In the extratropical latitudes (~30-60 N/S), cloud properties are strongly modulated by the occurrence of anti-cyclones (high pressure systems) and cyclones (low pressure systems). While deep clouds dominate in cyclones near the center and along the cold and warm fronts, low-level clouds populate the colder regions of the cyclones, and more specifically the region in the wake of the cold fronts, the post-cold frontal region (PCF). The GCMs mentioned above were found to underestimate the amount of these PCF clouds, and this caused errors in their representation of the amount of solar radiation at the surface of the southern hemisphere oceans, causing an excess in absorption and errors in long-term predictions. To better understand which processes are ill-represented in the GCMs to cause such issues, we used ground-based observations from the Eastern North Atlantic (ENA) ARM site and a high resolution regional model (called WRF) to examine the properties of PCF clouds. With the WRF model, we simulated the passage of a cold front at the ENA site using theoretically distinct representations of convection and the physics of the boundary layer. While these physical representations had little impact on the timing, structure and circulation of the cold front, we found a high sensitivity to how the large scale information was fed into the model. Errors rapidly developed if the incoming flow information was injected too far from the site, but 1000 km was found to be an optimal distance between the outer boundary and the ENA site for a realistic cold front passage. With this model, we then conducted a series of experiment where convection and boundary layer schemes were changed, and revealed that 1) the clouds were more sensitive to the convection than boundary layer schemes and 2) the interaction between the two has a large impact. Furthermore, changing the representation of the physics in the model had greater impacts on PCF clouds than changing the initial conditions. As large-scale models increase their spatial resolution, convection will eventually be resolved but it remains that their representation of the physics of the boundary layer will matter for the correct representation of shallow (low-level) clouds. Observations were used to find which large-scale information had a clear relationship with PCF cloud and precipitation properties, in view of helping with GCM parameterizations and evaluation. Both lower troposphere stability and surface forcing were found to be important for PCF clouds that experience more vigorous dynamics than their tropical or subtropical counterparts. For cloud base and top heights, a simple measure of the change in potential temperature between the surface and the 800 hPa was found to be a good predictor. This relationship was found to hold as well when using data from campaigns recently conducted in the southern ocean south of Tasmania. A recent GCM was found to succeed in representing similar relations. For cloud depth, water content, and different characteristics of precipitation both the temperature contrast between the surface and near-air and the near-surface wind speed had a strong relationship. To summarize, through observational constraints and regional model experiments, the project has provided guidance to evaluate GCMs and tools to improve their representation of low-level clouds in the extratropics. It has also highlighted differences in behavior between these clouds and those typically found in the subtropics, thus opening new avenues of research for shallow clouds in general.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Commercial building HVAC demand flexibility with model predictive control: Field demonstration and literature insights

Model Predictive Control (MPC) for building Heating Ventilation and Air Conditioning (HVAC) systems is beginning to gain traction in the market, with a few controls companies incorporating it into their product offerings. However, it remains difficult to assess whether the energy cost savings are enough to justify the cost of MPC implementation for a particular building, given the limited number of reported demonstrations. For small commercial and residential buildings with relatively uniform systems, standardized approaches can help lower implementation costs. In contrast, for large buildings or district systems, the potential magnitude of cost savings could justify more customized solutions. Estimating the cost-effectiveness of MPC becomes more challenging for medium and large commercial buildings, where a one-size-fits-all solution may not be suitable, and the potential energy cost savings may be insufficient to justify a customized solution. To make MPC technology more appealing, incorporating additional value streams beyond energy efficiency alone can significantly increase its attractiveness. One such revenue stream is demand flexibility, in response to dynamic electricity prices, where MPC can leverage the thermal mass of the building to shift the load and support the grid. Building on an extensive literature review of MPC field studies focused on cost savings and demand flexibility, this paper presents the results of implementing MPC control in a large office building HVAC system in Berkeley, CA. Four different dynamic electricity price profiles were integrated into the MPC objective function to shift building demand while maintaining comfort, and field testing was performed with each price profile across four seasons. The results show potential for 40–65 % demand decrease percentage and up to 61 % annual cost savings compared to the existing rule-based control strategy, under the tested dynamic price scenarios. This paper also presents a sensitivity analysis on the cost savings with respect to the price profile variability, discusses the implementation effort for the price-responsive MPC, and compares the cost savings found in this study to those found in literature on the basis of dynamic price variability, or so-called Electricity Price Relative Standard Deviation.

Zanetti, Ettore↗
DOE OSTI2022

Study of a Proposed fissionTPC Measurement of 233 U(n,f)/ 235 U(n,f) Cross Section Ratio

A review of the available 233 U( n, f )/ 235 U( n, f ) cross section ratio data compared to the ENDF evaluation in the fast neutron energy range, 0.1–20 MeV, shows that discrepancies exist between some of the various measurements and the evaluation. The discrepancies vary and a function of incident neutron energy and range between 1–5%. In this document we review two potential options for performing new measurements of the cross section ratio to potentially resolve these discrepancies. The two options are compared and contrasted in terms of measurement accuracy and cost. The first option would utilize the fissionTPC, which was designed specifically to achieve measurement uncertainties of better than 1% to resolve a similar data–evaluation discrepancy in the 239 Pu( n, f ) cross section. While this would seem to be an ideal option in this case also, there are challenges. The fissionTPC is a complex instrument and the cost of maintenance and operations are significant. Furthermore, some of the fissionTPC equipment is nearing the end of its useful life which presents a risk to completion of a measurement campaign. The second option is to use a dual-fission chamber, which is the more traditional approach to these types of measurements. The cost is estimated to be about half that of a fissionTPC measurement, while the much simpler design of the detector reduces overall risk. The trade-off is the measurement uncertainty, which would be on order of 1–3% with a careful measurement and analysis. This uncertainty however may not be good enough to fully resolve the discrepancies. The objective of a new measurement would be the improvement in the ENDF evaluation with the end goal of improving the accuracy of models and simulations at the end of the nuclear data pipeline. An understanding of the sensitivity of these models to the data uncertainties would help guide the choice of measurement techniques.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Biophysical model of eelgrass and water quality in Coos Bay, OR shows greater mitigation potential for ocean acidification than hypoxia

Seagrass beds provide important ecosystem services and are valued, in part, for their potential to mediate stressors such as ocean acidification and hypoxia (OAH) for sensitive species. However, the susceptibility of seagrasses to anthropogenic impacts and recent declines motivate the need to better understand the drivers of seagrass and the water quality consequences that occur with variation in seagrass abundance. To meet this need, we leveraged existing monitoring data (water quality and seagrass), hydrodynamic circulation model, and biogeochemical model framework with seagrass submodel, to produce a biophysical model of Coos Bay estuary, Oregon, U.S. The model includes biogeochemical processes involving water quality, plankton, seagrass, and sediment-water interactions. Ecosystem models like this are useful for evaluating complex estuarine systems because they allow us to extend our understanding of system dynamics beyond existing observations and perform experiments to identify the processes driving observed patterns. We used the biophysical model of Coos Bay to evaluate the dynamics of water quality and native eelgrass (Zostera marina) under three eelgrass abundance scenarios (zero eelgrass, current extent, and maximum observed extent) to elucidate the relationship between eelgrass and OAH. Including eelgrass in the Coos Bay model produced results that more closely resembled water quality observations - dissolved oxygen (DO) and pH were more dynamic in simulations with eelgrass, often having both higher highs and lower lows. While there were some areas of the estuary where DO improved with the addition of eelgrass to the model there was overall a small net increase in harmful DO conditions (based on a salmon physiological threshold). In contrast, ocean acidification conditions, pH and calcium carbonate saturation state for aragonite (Ω), were improved (based on oyster requirements) with the addition of eelgrass - although the magnitude of improvement differed seasonally and spatially. Our new model represents a useful tool - one which accounts for and controls the relevant physical and biogeochemical processes - to evaluate conditions that confer resilience or enhance vulnerability to OAH in an important Pacific Northwest coastal estuary and results can inform the OAH-related dynamics occurring in other eastern boundary current estuaries.

FVCOM-ICM↗
DOE OSTI2022

Trace Impurities Identified as Forensic Signatures in CMX-5 Fuel Pellets Using X-ray Spectroscopic Techniques

We report that small-particle analysis is a highly promising emerging forensic tool for analysis of interdicted special nuclear materials. Integration of microstructural, morphological, compositional, and molecular impurity signatures could provide significant advancements in forensic capabilities. We have applied rapid, high-sensitivity, hard X-ray synchrotron chemical imaging to analyze impurity signatures in two differently fabricated fuel pellets from the 5th Collaborative Materials Exercise (CMX5) of the IAEA Nuclear Forensics International Working Group. The spatial distributions, chemical compositions, and morphological and molecular characteristics of impurities were evaluated using X-ray absorption near-edge structure (XANES) and X-ray fluorescence chemical imaging to discover principal impurities, their granularity, particle sizes, modes of occurrence (distinct grains vs incorporation in the UO 2 lattice), and sources and mechanisms of incorporation. Differences in UO 2+x stoichiometry were detected at the microscale in nominally identical UO 2 ceramics (CMX5-A and CMX5-B), implying the presence of multiple UO 2 host phases with characteristic microstructures and feedstock compositions. Al, Fe, Ni, W, and Zr impurities and integrated impurity signature analysis identified distinctly different pellet synthesis and processing methods. For example, two different Al, W, and Zr populations in the CMX5-B sample indicated a more complex processing history than the CMX5-A sample. K-edge XANES measurements reveal both metallic and oxide forms of Fe and Ni but with different proportions between each sample. Altogether, these observations suggest multiple sources of impurities, including fabrication (e.g., force-sieving) and feedstock (mineral oxides). This study demonstrates the potential of synchrotron techniques to integrate different signatures across length scales (angstrom to micrometer) to detect and differentiate between contrasting UO 2 fuel fabrication techniques.

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

Colloidal Behavior of Plutonium Oxide in Concentrated Electrolyte Solutions

The Hanford Site in Washington State manages legacy high-level radioactive waste streams that display major chemistry and engineering challenges, including the high salt levels and pH values that correspond to conditions under which many classical concepts describing chemical reactivity cannot be applied. One particular challenge that needs to be tackled at Hanford is that Pu concentrations, [Pu], in the soluble phases of the tank wastes are higher than expected based on the solubility of crystalline PuO2, which is widely accepted to be caused by the formation of PuO2 colloid, consisting of nano- to submicron-sized particles (PuO2 NPs). Fundamental research underpinning the behavior of PuO2 NPs under conditions not only relevant to the Hanford tank waste but at high ionic strength in general is needed to reliably predict the chemical reactivity of PuO2 NPs and develop engineering solutions to safely and efficiently process high-level radioactive waste into forms suitable for long-term storage. In this work, we study the behavior of PuO2 NPs (particle size ~100 nm) under high ionic strength conditions by reacting it with highly concentrated (up to 5 M) salt solutions. We explore different electrolyte compositions to elucidate the impact of different anions (NO3-, Cl-, ClO4-, SO42-, C2O42-, CO32-) on the stability of PuO2 NP in the acidic and alkaline pH regime. PuO2 NP aggregation and precipitation as function electrolyte concentration is tracked by a combination of liquid scintillation counting, dynamic light scattering for determination of particle size distributions, and zeta potentials as a proxy for particle charge. At acidic pH, electrolytes containing non-coordinating anions, such as NaNO3, NaCl, and NaClO4 mostly stabilize PuO2 NPs over a large electrolyte concentration range, showing only subtle differences in their reactivity. Other electrolyte anions show a more pronounced effect on the PuO2 NP stability: SO42-, binds directly to the particles’ surface, reverses the particle charge, and precipitates the PuO2 NPs efficiently even at intermediate sulfate concentrations (>0.1 M). In contrast, C2O42- is found to lead to high [Pu] in solution, in the milli-molar range, even at mildly acidic pH (~4). Thermodynamic modeling of the dissolved Pu concentrations using PHREEQC is unable to predict the observed [Pu] in the acidic pH regime, supporting the influence of colloids in maintaining elevated [Pu]. It is noteworthy that the current thermodynamic databases do not include constants for colloidal Pu phases and cannot accurately predict many of the high ionic strength solutions relevant to this work. The mechanisms and models responsible for these observations will need further investigation in the future. At high pH values (~12), PuO2 NPs exhibits classical sol-gel chemistry, meaning that upon destabilization of the colloidal sol, for example by addition of concentrated NaOH, highly porous and viscid PuO2 coagulates are formed that consist of a three-dimensional network likely held together by physical interactions. The PuO2 NP coagulate shows no significant reversibility of the aggregation when contacted with concentrated brines; however, PuO2 NPs can be efficiently resuspended in solution by addition of diluted electrolytes, alkaline solutions containing high amounts of carbonate, or simple addition of water. Especially carbonate is shown to stabilize PuO2 NPs in solution at high pH, characterized by stable colloidal suspensions that are resistant against sedimentation during centrifugation. Thermodynamic modeling of the carbonate system was able to predict an increasing dissolved Pu concentration with increasing carbonate concentration. However, the model was profoundly sensitive to the fixed redox potential and does not include any thermodynamic constants for colloidal Pu species.

Neumann, Julia↗
DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗
DOE OSTI2023

Dark field X-ray microscopy below liquid-helium temperature: The case of NaMnO 2

Dark field X-ray microscopy (DFXM) is an experimental technique employed to investigate material properties by probing their ‘mesoscale,’ or microscale structures, in a bulk-sensitive manner using hard X-rays at synchrotron radiation sources. However, challenges remain when it comes to applications of this technique to examine low-temperature phenomena in quantum materials, which exhibit complex phase transitions at cryogenic temperatures. One such material is NaMnO 2 , which hosts an antiferromagnetic transition at 45 K that is suspected to coincide with local structural transitions from its majority monoclinic phase to nanoscale triclinic domains. Direct observation of local heterogeneities and this effect at low temperatures in NaMnO 2 is an important step in understanding this material and serves as an ideal candidate study for expanding the DFXM experimental design space. This paper details a foundational high-resolution DFXM study, down to liquid-helium temperature and below, conducted to explore phase transitions in NaMnO 2 . In conclusion, the outlined experiment ushers in the evaluation of other functional materials at low temperatures using this technique.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Influence of physical and chemical hydrology on bioremediation of a U-contaminated aquifer informed by reactive transport modeling incorporating 238 U/ 235 U ratios

Microbially-catalyzed reductive immobilization of aqueous uranium (U) as a solid phase has been proposed as a U remediation technique. Both laboratory and field experiments have demonstrated that this reduction reaction alters the 238 U/ 235 U ratio, producing a 238 U-enriched U(IV) solid. In contrast, other major U reactive transport processes fractionate these isotopes much less. This suggests the potential to quantify the extent of bioreduction occurring in groundwater containing U using the 238 U/ 235 U ratio, which would substantially improve upon current practices largely relying on U concentration measurements alone. Current reactive transport models for uranium dynamics include only concentration measurements, which are strongly influenced by highly coupled reactive transformation and aqueous transport processes. The complex physical and chemical behavior of U potentially compromises such quantitative analysis of U storage and release. Here we report the first numerical reactive transport model which explicitly incorporates variations in the 238 U/ 235 U ratio of U and demonstrates improved interpretation of the principal chemical reactions and groundwater transport processes affecting the subsurface mobility and distribution of this widespread contaminant. Recent U bioreduction studies performed in a contaminated aquifer in Rifle, Colorado, USA applied Rayleigh distillation models to interpret U stable isotope fractionation observed as a result of acetate amendment. These simplified models were unable to resolve the spatiotemporal pattern of U isotope fractionation recorded in the aqueous solutes. Here, we employ the multi-component, isotope-enabled CrunchTope reactive transport software to interpret these measured U isotope ratios, and demonstrate accurate reproduction of observed trends in both geochemistry and 238 U/ 235 U ratios for two consecutive years of field experiments. Overall, our results indicate that accurately modeling both the U concentration and isotope ratio distributions greatly constrains the parameter space of the model. We find that the transport properties of U in the Rifle aquifer are governed by the presence of low-permeability regions, which the isotopes are uniquely sensitive to. When U reduction is spatially constrained by these low-permeability regions, the shift in the 238 U/ 235 U ratio becomes more muted. Accurate modeling of observed U isotope ratios thus provides a powerful means to better understand bioremediation, and the current study serves to advance the application of this novel method.

238U/235U↗
DOE OSTI2022

Estimating SHmax azimuth with P sources and vertical geophones: Use P-P reflection amplitudes or use SV-P reflection times?

We compared two methods for extracting the azimuth of maximum horizontal stress (SHmax) from 3D land-based seismic data generated by a P source and recorded with vertical geophones. In the first method, we used the direct-SV mode that is produced by all land-based P sources. P sources generate SV illumination that radiates in all azimuth directions from a source station and creates SV-P reflections that are recorded by vertical geophones. Unless stratigraphy has steep dip, SV-P raypaths recorded by vertical geophones are the reverse of P-SV raypaths recorded by horizontal geophones. Thus, SV-P data provide the same S-wave sensitivity to stress fields as popular P-SV data do. In the second method, we retrieved P-P reflections and then performed an amplitude-variation-with-azimuth (AVA) analysis of the amplitude-gradient behavior of P-P reflection wavelets. We did this analysis in narrow azimuth corridors to determine the gradient of reflection-wavelet amplitudes as a function of azimuth. This P-P AVA amplitude-gradient method has been of great interest in the reflection seismology community since it was introduced in the late 1990s. Each of these methods, AVA analysis of the gradient of P-P reflection amplitudes and azimuth-dependent arrival times of SV-P reflections, can be used to determine the azimuth of SHmax stress. We compare the results of the two methods with ground truth measurements of SHmax azimuth at a CO 2 sequestration site in the Michigan Basin. SHmax azimuths were determined from P-P and SV-P data at three major boundaries at depths of approximately 3500 ft (1067 m), 5500 ft (1676 m), and 7500 ft (2286 m). Two estimates of SHmax azimuth (one using SV-P data and one using P-P data) were made at each stacking bin inside a 24 mi 2 (62 km 2 ) image space. The result was approximately 98,000 estimates of SHmax azimuth across each of these three boundaries for each of these two prediction strategies. Histogram displays of PP AVA gradient estimates had peaks at correct azimuths of SHmax at all three depths, but the spread of the distributions widened with depth and split into two peaks at the deepest boundary. In contrast, each histogram of SHmax azimuth predicted by azimuth-dependent SV-P traveltimes had a single, definitive peak that was positioned at the correct SHmax azimuth at all three boundary depths.

Geochemistry & Geophysics↗