Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “significant quantity”

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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Supplier's Status for Critical Solid Propellants, Explosive, and Pyrotechnic Ingredients

In the early 1970's a program was initiated at the Naval Surface Warfare Center/Indian Head Division (NSWC/IHDIV) to address the well-known problems associated with availability and suppliers of critical ingredients. These critical ingredients are necessary for preparation of solid propellants and explosives manufactured by the Navy. The objective of the program was to identify primary and secondary (or back-up) vendor information for these critical ingredients, and to develop suitable alternative materials if an ingredient is unavailable. In 1992 NSWC/IHDIV funded Chemical Propulsion Information Agency (CPIA) under a Technical Area Task (TAT) to expedite the task of creating a database listing critical ingredients used to manufacture Navy propellant and explosives based on known formulation quantities. Under this task CPIA provided employees that were 100 percent dedicated to the task of obtaining critical ingredient suppliers information, selecting the software and designing the interface between the computer program and the database users. TAT objectives included creating the Explosive Ingredients Source Database (EISD) for Propellant, Explosive and Pyrotechnic (PEP) critical elements. The goal was to create a readily accessible database, to provide users a quick-view summary of critical ingredient supplier's information and create a centralized archive that CPIA would update and distribute. EISD funding ended in 1996. At that time, the database entries included 53 formulations and 108 critical used to manufacture Navy propellant and explosives. CPIA turned the database tasking back over to NSWC/IHDIV to maintain and distribute at their discretion. Due to significant interest in propellant/explosives critical ingredients suppliers' status, the Propellant Development and Characterization Subcommittee (PDCS) approached the JANNAF Executive committee (EC) for authorization to continue the critical ingredient database work. In 1999, JANNAF EC approved the PDCS panel task. This paper is designed to emphasize the necessity of maintaining a JANNAF community supported database, which monitors PEP critical ingredient suppliers' status. The final product of this task is a user friendly, searchable database that provides a quick-view summary of critical ingredient supplier's information. This database must be designed to serve the needs of JANNAF and the propellant and energetic commercial manufacturing community as well. This paper provides a summary of the type of information to archive each critical ingredient.

Sims, B. L.↗
DOE OSTI2022

A Physics-Based DNI Model for Advancing Solar Resource Assessment and Forecasting: Preprint

Direct Normal Irradiance (DNI) is one of the most used quantities to quantify the magnitude of solar energy resource. The concept of DNI is often interpreted differently for ground measurements and solar forecasting by numerical weather prediction (NWP) models, leading to substantial bias during evaluation of DNI forecasts especially under cloudy-sky conditions. To eliminate the bias, we use the Fast All-sky Radiation Model for Solar applications with DNI (FARMS-DNI) to provide a physics-based solution of solar radiation in the circumsolar region. The FARMS-DNI is implemented in the Weather Research and Forecasting model with solar extensions (WRF-Solar) to forecast day-ahead DNI in the north America. By comparing with conventional predictions from WRF-Solar and satellite observations from the National Solar Radiation Data Base (NSRDB), we found significant improvements in our prediction of DNI.

DNI↗
DOE OSTI2023

A high-resolution large-eddy simulation framework for wildland fire predictions using TensorFlow

Background: Wildfires are becoming more severe, so we need improved tools to predict them over a wide range of conditions and scales. One approach towards this goal entails the use of coupled fire/atmosphere modelling tools. Although significant progress has been made in advancing their physical fidelity, existing tools have not taken full advantage of emerging programming paradigms and computing architectures to enable high-resolution wildfire simulations. Aims: The aim of this study was to present a new framework that enables landscape-scale wildfire simulations with physical representation of combustion at an affordable cost. Methods: We developed a coupled fire/atmosphere simulation framework using TensorFlow, which enables efficient and scalable computations on Tensor Processing Units. Key Results: Simulation results for a prescribed fire were compared with experimental data. Predicted fire behavior and statistical analysis for fire spread rate, scar area, and intermittency showed overall reasonable agreement. Scalability analysis was performed, showing close to linear scaling. Conclusions: While mesh refinement was shown to have less impact on global quantities, such as fire scar area and spread rate, it benefits predictions of intermittent fire behavior, buoyancy-driven dynamics, and small-scale turbulent motion. Implications: This new simulation framework is efficient in capturing both global quantities and unsteady dynamics of wildfires at high spatial resolutions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Reduced-Order Modeling of Hydrogen Releases from Vent Stacks and with Wind Effects

Here, the physical release behavior of hydrogen is important to understand from a safety and design perspective. The consequences of unignited pressurized gaseous hydrogen plumes exiting vent stacks were considered by extending and modifying existing general hydrogen plume models. Entrainment, vent stack backpressure, and the flow regime of hydrogen exiting the vent were found to be significant factors affecting plume shape and size, but further investigation and validation with unchoked, low-Froude-number flows is recommended to improve the model’s robustness. Additionally, models for the effects of wind on unignited plume momentum and entrainment were added to explore this behavior. Wind was assumed to increase mixing of hydrogen with the ambient air, and to affect the momentum of the released jet. Introducing wind into the plume model led to a shorter plume for all wind and jet directions. A high counter-flowing wind led to non-physical results and challenges in interpreting the visualization. The proposed jet plume wind sub-models (specifically entrainment coefficients) were fit and compared to experimental data of different releases of hydrogen into a wind tunnel, but the quantity of data available and experimental conditions were limited. Thus, collection of more empirical data and for a wider range of conditions is recommended for improvement of the proposed computational models. Developing reduced-order models for these physical phenomena can improve accessibility to predicted physical behavior and the rate at which hydrogen systems can be safely designed and deployed.

entrainment↗
DOE OSTI2020

Status, Challenges, and Plans for Consistent Modeling and Evaluation of Fission Data: Nubar, PFNS, and FPY [Slides]

This presentation discusses the importance of consistent evaluations for fission observables. Currently, prompt fission quantities (e.g., neutron multiplicity, neutron energy spectrum, and fission product yields) are evaluated independently of one another, either with separate models or data-only/non-model evaluations (few – if any – shared model inputs) and inconsistencies can arise in evaluated data (e.g. NSE 190, 258 (2018)). They have found that consistent, model-based evaluations lead to more physical constraints on free parameters in the model, more robust predictions for other observables and other isotopes, that there are correlations between the uncertainties on different observables, and that uncertainties can be propagated to unmeasured quantities. They are working on developing models to the point where fission observables are calculated at the quality needed for an evaluation. They are striving to consistently calculate fission observables, such as $\bar{v}$, PFNS, and FPY, to a quality suitable for an evaluation. Current evaluations for fission observables use separate models for each observable (with minimal shared input) or are purely data driven. LANL models, such as the Monte Carlo CGMF and deterministic BeoH, have different strengths; the similar fission fragment initial conditions and decay models leads to consistency between the two codes (FPY can be connected to $\bar{v}$ and PFNS, even when calculated by different codes). Work is underway by Lovell to optimize the models and Neudecker is working to perform detailed uncertainty quantification. The PFNS, in particular, presents a significant challenge.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2006

Experimental Partitioning of Cr(3+) and Sc(3+) into Olivine: Mechanisms and Implications

Olivine (Mg, Fe)Si2O4 does not, by stoichiometry, accept cations such as Sc(3+) or Cr(3+). However, the partition coefficients of Sc and Cr between olivine and liquid are significant 0.2-1.0. We have measured Cr(3+) partition coefficients of near unity and have grown olivines with nearly 3 wt.% Sc2O3. Therefore, there must be a simple means of charge balancing 3+ ions in a crystal structure that was obviously not designed to receive other than 2+ ions on the olivine M sites. The simplest explanation is that two 3+ ions enter the olivine structure by displacing three 2+ ions and creating an M site vacancy. Even this explanation has difficulties. For minor elements in our experiments (~1 wt.%) the odds of a minor element 3+ ion finding a second 3+ for charge balance are of the order of 100:1 against. Because of the reducing conditions of our experiments, Fe(3+) will not suffice; and Al(3+) is not in sufficient quantity in olivine for charge balance. Therefore, Cr or Sc must, in effect, charge balance itself. For true trace elements, the problem is compounded many times. For an ion at the 10 ppm level the chances of finding a second (for example) Sc ion is approx.10(exp 5):1 against. Of course, any other 3+ ion would suffice but comparisons between percent level doping experiments and trace level partitioning indicate that Henry s law is obeyed. This implies that the same substitution mechanism occurs at both the percent and tens of ppm levels. There are two simple solutions to this problem: (i) The electrical conductivity of olivine is such that charge balance need not be local. This requires substantial domains within the olivine crystal in electrical contact by migration of vacancies or electronic defects. (ii) The 3+ cation brings along its own charge-balancing ion because it existed as a dimer in the silicate liquid. Olivine is not a true insulator but is actually a p-type semiconductor. Even so, electrical communication by this means is unlikely over the tens or hundreds of unit cells that would be required for charge balance to be local. Therefore, we cautiously favor the idea that melt speciation is the means by which 3+ ions enter the olivine structure. Possibly this model might be tested by in situ XAFS measurements or by molecular dynamical calculations.

Jones, John↗
DOE OSTI2023

Validation of the single-event method for low-energy electron transport via stopping power calculations with $\mathrm{MCNP}$

Monte Carlo simulations of low-energy ( <50 keV) electron transport in matter are essential for a broad range of application fields. Several Monte Carlo codes have developed specialized treatments for this case, but a comprehensive validation of low-energy electron transport for general-purpose simulations remains lacking in the literature. One approach to accomplish this validation is calculation of stopping power using low-energy electron transport physics, as stopping power is a fundamental radiation transport quantity which must be simulated accurately for nearly any application. Here in this work, we use the Monte Carlo N-Particle (MCNP) radiation transport code with the single-event method for electron transport to calculate stopping powers of low-energy electrons (50 eV to 30 keV) in 41 elemental solids, 14 compound solids, and five rare gas solids, comparing simulation results to published semi-empirical stopping power calculations from optical measurements. In general, the simulations give good agreement (typically within ±10%) with semi-empirical stopping power values at higher energies: 300 eV and above for most elemental solids, 1 keV and above for compound solids, and 400 eV and above for rare gas solids. Agreement between MCNP and semi-empirical values is worse below these energies. The most significant source of error is the EPRDATA14 cross section data, which does not account for changes in electronic structure due to solid-state bonding, particularly in compound materials. The simplistic model of atomic excitation used to generate the EPRDATA14 cross sections is another key source of error. Additionally, the breakdown of the continuous slowing-down approximation introduces significant uncertainty at low energies, although this is a limitation of the calculation method and not of the simulation procedure. Accounting for these and other uncertainty sources, the single event method in MCNP is robust and able to give good accuracy for a variety of low-energy electron transport problems through diverse kinds of materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRSDate not supplied

Exploring the Capabilities of a Machine Learning Algorithm to Detect Space Weather-Significant Emerging Active Regions

Active regions are a source of various phenomena responsible for Space Weather disturbances; therefore, developing a technology for early warning about upcoming magnetic activity is crucial to mitigate its impact. However, observational limitations and the high nonlinearity of processes associated with the accumulation of magnetic flux and its interaction with the surrounding plasma during the emergence through the convection zone make early activity detection a challenging problem. To address these challenges, we developed a physics-driven machine learning model that allows us to detect active regions (ARs) before they become visible on the solar surface by analyzing the power spectra of acoustic oscillations observed by the SDO/HMI instrument. This study is based on a time series of Doppler shift maps of 31x31-degree areas tracked with the Carrington rotation rate for four days before and after the emergence. The Doppler shift time series are processed into the oscillation power maps for four frequency ranges and accompanied by line-of-sight magnetograms and the continuum intensity maps from SDO/HMI. The resulting data are converted into a 1D time series representing the mean temporal variations of these quantities. The redacted time series are used as input to predict AR emergence using the Long Short Term Memory (LSTM) method. The training of the LSTM model is based on 40 ARs, which includes an independent analysis for each sub region that exhibits AR emergence or remains quiet. The emergence of magnetic flux (defined as a decrease of the continuum intensity) was detected with the developed LSTM algorithm from 5 to 48 hours before the reported time by NOAA. The developed model is capable of pointing to the time and location of active region formation. In this presentation, we discuss reasons that impact how early in advance the model can identify the upcoming activity and the possibility of improving the current predictive skills and steps to transition to the operational forecast.

Heliophysics↗
NASA NTRS1975

Correlation of boundary layer quantities for hypersonic laminar flows with zero pressure gradient for several gases

The laminar boundary layer has been theoretically studied for six gases for flows over cold walls with zero pressure gradient at Mach numbers between 5.5 and 12.5 to correlate boundary layer quantities for the various gases. The flow conditions considered correspond to those that can be generated in test facilities such as the shock tunnel and the expansion tube. Computed results obtained using real gas properties indicate that the Eckert number based on edge conditions serves to correlate the results in terms of the wall shear stress and enthalpy gradient, the Stanton number, and the momentum thickness for the various gases within plus or minus 10 per cent for Te = Tw and Te approximately 3Tw. Computed Reynolds analogy factors exhibit very good agreement with those predicted by the Colburn analogy. Velocity and displacement thicknesses correlate well with Eckert number for Te = Tw, but fail to correlate for Te approximately 3Tw. Differences in results are traced to property variations. Results show that the Eckert number is a significant correlating variable for the flows considered.

Cook, W. J.↗
NASA NTRS1994

Gravity model improvement using the DORIS tracking system on the SPOT 2 satellite

A high-precision radiometric satellite tracking system, Doppler Orbitography and Radio-positioning Integrated by Satellite system (DORIS), has recently been developed by the French space agency, Centre National d'Etudes Spatiales (CNES). DORIS was designed to provide tracking support for missions such as the joint United States/French TOPEX/Poseidon. As part of the flight testing process, a DORIS package was flown on the French SPOT 2 satellite. A substantial quantity of geodetic quality tracking data was obtained on SPOT 2 from an extensive international DORIS tracking network. These data were analyzed to assess their accuracy and to evaluate the gravitational modeling enhancements provided by these data in combination with the Goddard Earth Model-T3 (GEM-T3) gravitational model. These observations have noise levels of 0.4 to 0.5 mm/s, with few residual systematic effects. Although the SPOT 2 satellite experiences high atmospheric drag forces, the precision and global coverage of the DORIS tracking data have enabled more extensive orbit parameterization to mitigate these effects. As a result, the SPOT 2 orbital errors have been reduced to an estimated radial accuracy in the 10-20 cm RMS range. The addition of these data, which encompass many regions heretofore lacking in precision satellite tracking, has significantly improved GEM-T3 and allowed greatly improved orbit accuracies for Sun-synchronous satellites like SPOT 2 (such as ERS 1 and EOS). Comparison of the ensuing gravity model with other contemporary fields (GRIM-4C2, TEG2B, and OSU91A) provides a means to assess the current state of knowledge of the Earth's gravity field. Thus, the DORIS experiment on SPOT 2 has provided a strong basis for evaluating this new orbit tracking technology and has demonstrated the important contribution of the DORIS network to the success of the TOPEX/Poseidon mission.

Nerem, R. S.↗
NASA NTRS2014

How Accurate is Land/Ocean Moisture Transport Variability in Reanalyses?

Quantifying the global hydrological cycle and its variability across various time scales remains a challenge to the climate community. Direct measurements of evaporation (E), evapotranspiration (ET), and precipitation (P) are not feasible on a global scale, nor is the transport of water vapor over the global oceans and sparsely populated land areas. Expanding satellite data streams have enabled development of various water (and energy) flux products, complementing reanalyses and facilitating observationally constrained modeling. But the evolution of the global observing system has produced additional complications--improvements in satellite sensor resolution and accuracy have resulted in "epochs" of observational quasi-uniformity that can adversely affect reanalysis trends. In this work we focus on vertically integrated moisture flux convergence (VMFC) variations within the period 1979 - present integrated over global land. We show that VMFC in recent reanalyses (e.g. ERA-I, NASA MERRA, NOAA CFSR and JRA55) suffers from observing system changes, though differently in each product. Land Surface Models (LSMs) forced with observations-based precipitation, radiation and near-surface meteorology share closely the interannual P-ET variations of the reanalyses associated with ENSO events. (VMFC over land and P-ET estimates are equivalent quantities since atmospheric storage changes are small on these scales.) But the long-term LSM trend over the period since 1979 is approximately one-fourth that of the reanalyses. Additional reduced observation reanalyses assimilating only surface pressure and /or specifying seasurface temperature also have a much smaller trend in P-ET like the LSMs. We explore the regional manifestation of the reanalysis P-ET / VMFC problems, particularly over land. Both principal component analysis and a simple time series changepoint analysis highlight problems associated with data poor regions such as Equatorial Africa and, for one reanalysis, the Equatorial Andes region. Onset of the availability of passive microwave Special Sensor Microwave Imager (SSMI) moisture data in July 1987 and the transition from the Microwave Sounder Unit (MSU) to an advanced version (AMSU) have significant impacts on VMFC variability. Simple accounting for these errors of leading importance results in modified reanalysis VMFC estimates that agree much better with the LSM results. Regional details of the modified reanalysis VMFC and LSM P-ET are related to changes in Pacific Decadal Variability as manifest in SST changes after the late 1990s.

Robertson, F. R.↗
DOE OSTI2021

Crossover Effects in Batteries with High–Nickel Cathodes and Lithium–Metal Anodes

It is well understood that cathode-to-anode crossover, especially of transition-metal ions, can significantly impact the long-term cycling of lithium-ion batteries. The dissolved transition-metal ions in lithium-ion cells deposit on the graphite anode, disrupt the solid-electrolyte interphase (SEI), and catalyze further side reactions. Meanwhile, crossover effects in lithium-metal batteries have rarely been studied. This study is the first to investigate crossover effects in lithium-metal batteries with high-nickel layered-oxide cathodes. It is shown that the crossover of transition-metal ions from LiNi 0.9 Mn 0.05 Co 0.05 O 2 has minimal effect on the lithium-metal anode (LMA) due to the following reasons. The catalytic transition metals 1) have less effect on an inherently reactive LMA, 2) are diluted in a thicker SEI, and 3) are produced in overall lower quantity due to the limited cycle life of the LMA. Conversely, the LMA generates soluble decomposition products that cross over to the cathode even during early cycling. This crossover accelerates impedance growth and capacity fade at the cathode and is partially responsible for the mismatch between the performance of half and full-cells with layered-oxide cathodes. Finally, this study highlights the need for better battery design with LMA, potentially including electrolyte or cell modifications.

25 ENERGY STORAGE↗
DOE OSTI2024

Probing thermodynamics of radiogenic helium and defects in $δ$-plutonium alloys and interactions with adsorbed environmental gases

Differential scanning calorimetry coupled with simultaneous evolved gas analysis (DSC-EGA) on aged δ-Pu samples shows that most radiogenic helium remains trapped within the Pu matrix at temperatures very close to, or slightly above, the melting temperature. Our results indicate that helium release from 50-year-old δ-Pu occurs as a burst just below the melting temperature (>0.994 T m ), with subsequent pressure oscillations as temperature increases. Subordinate quantities of H 2 were also released along with helium. Here, the helium emission tails off and ceases above ~720 – 750°C. In a δ-Pu alloy aged 6 years, the initial helium burst occurs slightly above melting (~1.015 – 1.042 T m ), with a discrete, larger helium spike occurring between 670 and 686°C. The proximity of helium release to the liquidus transition presented challenges in the deconvolution of overlapping process enthalpies, the liquidus endotherm, and the exotherm resulting from bubble collapse, annealing and gas expulsion. Helium’s strong affinity for vacancy binding in a 2He-vac configuration is predicted by Density Functional Theory (DFT) modeling. The measured stored energy associated with the He release events in a 50-year-old δ alloy is on the order of ~ 10–11 J/g, which is significantly higher than stored energies measured in the sub-solidus regimes (~2 J/g) that are related to the solid-state annealing of processing- and radiation-induced defects. This implies that aged δ Pu alloys have a remarkable resilience to accommodate the lattice strain produced by the internal pressure of the helium bubbles and provides further insight into the thermodynamic behavior of aged δ Pu.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Improved Boundary Conditions for Coupled Geospace Models: An Application in Modeling Spacecraft Surface Charging Environment

Abstract Spacecraft surface charging in the inner magnetosphere often occurs in the pre‐midnight to the dawn sector when electron fluxes of tens of keV increase. Inner magnetosphere ring current models can be used to simulate/predict the spacecraft surface charging environment, with their outer boundary conditions specified either based on observations or provided by other models, such as MHD models. In the latter approach, using MHD quantities, the flux spectrum at the outer boundary is commonly assumed to follow a Kappa or Maxwellian distribution function. Such a method however often departs greatly from the realistic spectrum at E < tens of keV, a crucial energy range in the surface charging anomaly. In order to achieve a better representation of the surface charging environment, we propose to combine the MHD‐parameterized flux spectrum with an empirical electron flux model of E < 40 keV to set the electron flux boundary condition. Results indicate that as opposed to the case where the MHD‐parameterized flux distribution is solely used at the model boundary, simulations with the new boundary condition yields a more intense surface charging environment. The integrated electron flux between 10 < E < 50 keV, a measure of the severity of the surface charging environment, is significantly enhanced by 1‐2 orders of magnitude, leading to a much better agreement with Van Allen Probes measurements. This study hence demonstrates a reasonable solution to the setting of outer boundary conditions for inner magnetosphere models and is recommended for coupled geospace circulation models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Active Learning A Neural Network Model For Gold Clusters & Bulk From Sparse First Principles Training Data

Small metal clusters are of fundamental scientific interest and of tremendous significance in catalysis. These nanoscale clusters display diverse geometries and structural motifs depending on the cluster size; a knowledge of this size-dependent structural motifs and their dynamical evolution has been of longstanding interest. Given the high computational cost of first-principles calculations, molecular modeling and atomistic simulations such as molecular dynamics (MD) has proven to be an important complementary tool to aid this understanding. Classical MD typically employ predefined functional forms which limits their ability to capture such complex size-dependent structural and dynamical transformation. Neural Network (NN) based potentials represent flexible alternatives and in principle, well-trained NN potentials can provide high level of flexibility, transferability and accuracy on-par with the reference model used for training. A major challenge, however, is that NN models are interpolative and requires large quantities (similar to 10 4 or greater) of training data to ensure that the model adequately samples the energy landscape both near and far-from-equilibrium. A highly desirable goal is minimize the number of training data, especially if the underlying reference model is first-principles based and hence expensive. In this work, we introduce an active learning (AL) scheme that trains a NN model on-the-fly with minimal amount of first-principles based training data. Our AL workflow is initiated with a sparse training dataset (similar to 1 to 5 data points) and is updated on-the-fly via a Nested Ensemble Monte Carlo scheme that iteratively queries the energy landscape in regions of failure and updates the training pool to improve the network performance. Using a representative system of gold clusters, we demonstrate that our AL workflow can train a NN with similar to 500 total reference calculations. Using an extensive DFT test set of similar to 1100 configurations, we show that our AL-NN is able to accurately predict both the DFT energies and the forces for clusters of a myriad of different sizes. Our NN predictions are within 30 meV/atom and 40 meV/angstrom of the reference DFT calculations. Moreover, our AL-NN model also adequately captures the various size-dependent structural and dynamical properties of gold clusters in excellent agreement with DFT calculations and available experiments. We finally show that our AL-NN model also captures bulk properties reasonably well, even though they were not included in the training data.

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

Manganese: The overlooked contaminant in the world largest mine tailings dam collapse

Manganese (Mn) is an abundant element in terrestrial and coastal ecosystems and an essential micronutrient in the metabolic processes of plants and animals. Mn is generally not considered a potentially toxic element due to its low content in both soil and water. However, in coastal ecosystems, the Mn dynamic (commonly associated with the Fe cycle) is mostly controlled by redox processes. Here, we assessed the potential contamination of the Rio Doce estuary (SE Brazil) by Mn after the world’s largest mine tailings dam collapse, potentially resulting in chronic exposure to local wildlife and humans. Estuarine soils, water, and fish were collected and analyzed seven days after the arrival of the tailings in 2015 and again two years after the dam collapse in 2017. Using a suite of solid-phase analyses including X-ray absorption spectroscopy and sequential extractions, our results indicated that a large quantity of Mn II arrived in the estuary in 2015 bound to Fe oxyhydroxides. Over time, dissolved Mn and Fe were released from soils when Fe III oxyhydroxides underwent reductive dissolution. Due to seasonal redox oscillations, both Fe and Mn were then re-oxidized to Fe III , Mn III , and Mn IV and re-precipitated as poorly crystalline Fe oxyhydroxides and poorly crystalline Mn oxides. In 2017, redox conditions (Eh: -47 ± 83 mV; pH: 6.7 ± 0.5) favorable to both Fe and Mn reduction led to an increase (~880%) of dissolved Mn (average for 2015: 66 ± 130 μg L -1 ; 2017: 582 ± 626 μg L -1 ) in water and a decrease (~75%, 2015: 547 ± 498 mg kg -1 ; 2017: 135 ± 80 mg kg -1 ) in the total Mn content in soils. The crystalline Fe oxyhydroxides content significantly decreased while the fraction of poorly ordered Fe oxides increased in the soils limiting the role of Fe in Mn retention. The high concentration of dissolved Mn found within the estuary two years after the arrival of mine tailings indicates a possible chronic contamination scenario, which is supported by the high levels of Mn in two species of fish living in the estuary. Our work suggests a high risk to estuarine biota and human health due to the rapid Fe and Mn biogeochemical dynamic within the impacted estuary.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Improving CLM5.0 Biomass and Carbon Exchange Across the Western United States Using a Data Assimilation System

The Western United States is dominated by natural lands that play a critical role for carbon balance, water quality, and timber reserves. This region is also particularly vulnerable to forest mortality from drought, insect attack, and wildfires, thus requiring constant monitoring to assess ecosystem health. Carbon monitoring techniques are challenged by the complex mountainous terrain, thus there is an opportunity for data assimilation systems that combine land surface models and satellite-derived observations to provide improved carbon monitoring. Here, we use the Data Assimilation Research Testbed to adjust the Community Land Model (CLM5.0) with remotely sensed observations of leaf area and above-ground biomass. The adjusted simulation significantly reduced the above-ground biomass and leaf area, leading to a reduction in both photosynthesis and respiration fluxes. The reduction in the carbon fluxes mostly offset, thus both the adjusted and free simulation projected a weak carbon sink to the land. This result differed from a separate observation-constrained model (FLUXCOM) that projected strong carbon uptake to the land. Simulation diagnostics suggested water limitation had an important influence upon the magnitude and spatial pattern of carbon uptake through photosynthesis. We recommend that additional observations important for water cycling (e.g., snow water equivalent, land surface temperature) be included to improve the veracity of the spatial pattern in carbon uptake. Furthermore, the assimilation system should be enhanced to maximize the number of the simulated state variables that are adjusted, especially those related to the recommended observed quantities including water cycling and soil carbon.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Dependence of the boundary heat flux width on core and edge profiles in Alcator C-Mod

This work presents new evidence that the heat flux width, λ q , in the Alcator C-Mod tokamak scales with the edge electron pressure, as observed in the ASDEX Upgrade (AUG) tokamak, but the scaling with volume-averaged pressure, $\bar{p}$, from the plasma stored energy, is a better predictor of λ q in Alcator C-Mod than the edge electron pressure. These previous studies, which find that λ q decreases with increasing plasma pressure, imply that a high performance core at high pressure will lead to challenging heat and particle exhaust due to very small λ q . This concern has led to our significant enlargement of the C-Mod database with the electron density, temperature, and pressure profile data from the Thomson scattering and electron cyclotron emission diagnostics. Using the C-Mod database augmented with new profile data, we find that λ q decreases with increasing edge electron pressure as ${\lambda }_{q}\propto {p}_{\mathrm{e},95}^{-0.26}$, similar to results from AUG, and showing the strength of cross-machine comparisons. We also find that ${\lambda }_{q}\propto {p}_{\mathrm{e},\mathrm{c}\mathrm{o}\mathrm{r}\mathrm{e}}^{-0.56}$, consistent with the original finding from C-Mod that the heat flux width scales as ${\bar{p}}^{-0.48}$. The scalings of λ q with separatrix pressure and gradient scale length are found to match the AUG results qualitatively. The C-Mod scalings with edge plasma quantities have more scatter than the $\bar{p}$ scaling, and, importantly, show different trends for H-modes relative to L- and I-mode. Furthermore, investigating the source of this discrepancy presents an opportunity for further study that may improve our ability to predict the heat flux width in different confinement scenarios in the pursuit of optimizing core-edge performance in future reactors.

Alcator C-Mod↗