Search NASA⌕ Search

SEARCH · Search NASA

Results for “Atmospheric blocking”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

kynema-sgf [SWR-20-85]

Kynema-SGF (formerly AMR-Wind), wherein SGF stands for structured-grid fluid dynamics, is a massively parallel, block-structured adaptive-mesh, incompressible flow solver. The codebase was initiated in 2019 from incflo. The solver is built on top of the AMReX library. AMReX library provides the mesh data structures, mesh adaptivity, as well as the linear solvers used for solving the governing equations. Kynema-SGF is actively developed and maintained by a dedicated multi-institutional team from Lawrence Berkeley National Laboratory, National Laboratory of the Rockies, and Sandia National Laboratories. The primary applications for Kynema-SGF are: performing large-eddy simulations (LES) of atmospheric boundary layer (ABL) flows, simulating wind farm turbine-wake interactions using actuator disk or actuator line models for turbines, and as a background solver when coupled with a near-body solver (e.g., Kynema-UGF) with overset methodology to perform blade-resolved simulations of multiple wind turbines within a wind farm. For offshore applications, the ability to model the air-sea interaction effects and its impact on the ABL characteristics is another focus for the code development effort. As with other codes in the Kynema ecosystem, Kynema-SGF shares the following objectives: *an open, well-documented implementation of the state-of-the-art computational models for modeling wind farm flow physics at various fidelities that are backed by a comprehensive verification and validation (V&V) process; *be capable of performing the highest-fidelity simulations of flow fields within wind farms; and *be able to leverage the high-performance leadership class computing facilities available at DOE national laboratories.

Ananthan, Shreyas↗
DOE OSTI2025

Unsupervised discovery of extreme weather events using universal representations of emergent organization

Spontaneous self-organization is ubiquitous in systems far from thermodynamic equilibrium. While organized structures that emerge dominate transport properties, universal representations that identify and describe these key objects remain elusive. Here, we introduce a theoretically grounded framework for describing emergent organization that, via data-driven algorithms, is constructive in practice. Its building blocks are spacetime lightcones that embody how information propagates across a system through local interactions. We show that predictive equivalence classes of lightcones—local causal states—capture organized behaviors in complex spatiotemporal systems. Employing an unsupervised physics-informed machine learning algorithm and a high-performance computing implementation, we demonstrate automatically discovering organized structures in two real-world domain science problems. We show that local causal states identify vortices and track their power-law decay behavior in two-dimensional fluid turbulence. We then show how to detect and track familiar extreme weather events—hurricanes and atmospheric rivers—and discover other novel structures associated with precipitation extremes in high-resolution climate data at the grid-cell level.

Rupe, Adam [Pacific Northwest National Laboratory ↗
DOE OSTI2021

CACTI ARM Aerial Facility Measurements of Ice Nucleating Particles

The dataset comprises measures, using Colorado State University's Ice Spectrometer (IS, an immersion freezing device with a range from 0°C down to -26 to -29°C) of atmospheric ice nucleating particle (INP) concentrations taken on the Atmospheric Radiation Measurement (ARM) program Aerial Facilty (AAF) G-1 aircraft. INP measurements on the G-1 were collected from varied altitudes on different flights over the region of the Sierras de Córdoba mountain range of north-central Argentina, centred over ARM's Mobile Facility (AMF-1) near Villa Yacanto, where ground-based INP measures were being taken concurrently. Both studies took place as part of the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) Experiment. A filter sampling system was deployed on the G-1 to collect aerosol particles for post-processing of INPs collected on filters (measuring their immersion freezing ability) once returned to Colorado State University (CSU). Filter holders used were pre-cleaned aluminum in-line units loaded with pre-cleaned and pre-sterilized 47 mm dia. Nuclepore polycarbonate filters (0.2 mm pore size). Filters were drawn for varied times, resulting in varied volumes collected (51 to 1667 SL). Mass flow rate was recorded in real-time so that total sampled volume (at standard temperature and pressure) could be determined. A total of 34 sample filters were collected over the IOP, including 5 blanks. Filters were stored at -20°C freezer prior to frozen return to Colorado State University (CSU). . Processing to obtain spectra of INP number concentration active via the immersion freezing mechanism versus temperature was conducted using CSU's IS instrument (McCluskey et al., 2018). For measurment of INPs, collected aerosol particles were re-suspended in 7 mL of 0.02 µm-filtered deionized water. Aliquots of each suspension, and serial dilutions, were dispensed into trays which were fit into aluminum blocks in the IS. Samples are cooled at 0.33°C min-1 and the freezing temperatures of wells recorded automatically. Cumulative INP concentrations were determined by first calculating the INPs per mL of suspension based on Vali (1971) and then converting to concentration per standard liter of air using the proportion of the total liquid sample dispensed and the air sample volumes. Aliquots of suspensions from selected samples were also heat treated (95°C for 20 min) to denature and deactivate biological INPs, and digested in 10% H2O2 at 95°C under UV-B to remove all organic carbon INPs. McCluskey, C. S., J. Ovadnevaite, M. Rinaldi, J. Atkinson, F. Belosi, D. Ceburnis, S. Marullo, T. C. J. Hill, U. Lohmann, Z. A. Kanji, C. O’Dowd, S. M. Kreidenweis, P. J. DeMott, 2018: Marine and Terrestrial Organic Ice Nucleating Particles in Pristine Marine to Continentally-Influenced Northeast Atlantic Air Masses, Journal of Geophysical Research: Atmospheres, 123, 6196–6212, https://doi.org/10.1029/2017JD028033. Vali, G., 1971: Quantitative evaluation of experimental results on the heterogeneous freezing nucleation of supercooled liquids. J. Atmos. Sci., 28, 402–409.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Cost analysis of alternative large-scale high-temperature solid oxide electrolysis hydrogen production facilities

We extend our past cost analysis of gigawatt-scale solid oxide electrolysis (SOE) facilities that produce high purity hydrogen gas from water by estimating construction and operating costs for three new alternative design cases: (1) offsite feed steam generation; (2) near-atmospheric pressure (NAP) stack; and (3) onsite electric boiler feed steam generation. Pressure effects on hydrogen electrode-(cathode-)supported SOE cell (SOEC) stack performance are estimated for the same assumed cell and stack construction and used to determine facility-wide stack capital costs for achieving a fixed H2 production at different pressures. Modular balance of plant (BOP) process equipment capital costs are estimated for each new alternative design case using our past equipment sizing, design, and cost data and scaling relationships. Furthermore, we update BOP equipment sizing and design for the NAP case using Aspen®. Vendor quotes for electric boilers are used to estimate costs for the electric boiler design case. Factory and onsite assembly and installation costs for SOEC stacks and BOP equipment are calculated using our past simplified first-principles approach. First-of-a-kind (FOAK) and N th -of-a-kind (NOAK) production maturity cost estimates are included for all results. The case with NAP stacks offers the lowest facility total capital cost (TCC, ~23% lower than base) while use of small electric boilers requires the highest TCC (~3% higher than base). H 2 production prices decrease from the base of ~$\$2.17$ /kgH 2 to ~$\$1.92$/kgH 2 for 1 GW e DC SIP facilities utilizing NAP stacks supplied by offsites steam situated in large modules and blocks for $\$0.030$/kWh e and $0.009/kWh t prices for electricity and thermal energy, respectively. We report all costs in 2021 US dollars.

Balance of plant (BOP) process equipment↗
DOE OSTI2023

Synthesis of Free-Standing Pd-Ni-P Metallic Glass Nanoparticles with Durable Medium-Range Ordered Structure for Enhanced Electrocatalytic Properties

Topologically disordered metallic glass nanoparticles (MGNPs) with highly active and tailorable surface chemistries have immense potential for functional uses. The synthesis of free-standing MGNPs is crucial and intensively pursued because their activity strongly depends on their exposed surfaces. Herein, a novel laser-evaporated inert-gas condensation method is designed and successfully developed for synthesizing free-standing MGNPs without substrates or capping agents, which is implemented via pulse laser-induced atomic vapor deposition under an inert helium atmosphere. In this way, the metallic atoms vaporized from the targets collide with helium atoms and then condense into short-range-order (SRO) clusters, which mutually assemble to form the MGNPs. Using this method, free-standing Pd 40 Ni 40 P 20 MGNPs with a spherical morphology are synthesized, which demonstrates satisfactory electrocatalytic activity and durability in oxygen reduction reactions. Moreover, local structure investigations using synchrotron pair distribution function techniques reveal the transformation of SRO cluster connection motifs of the MGNPs from face-sharing to edge-sharing modes during cyclic voltammetry cycles, which enhances the electrochemical stability by blocking crystallization. Further, this approach provides a general strategy for preparing free-standing MGNPs with high surface activities, which may have widespread functional applications

36 MATERIALS SCIENCE↗
DOE OSTI2020

Arctic Coastal Erosion: Modeling and Experimentation

Increasing Arctic coastal erosion rates have put critical infrastructure and native communities at risk while also mobilizing ancient organic carbon into modern carbon cycles. Although the Arctic comprises one-third of the global coastline and has some of the fastest eroding coasts, current tools for quantifying permafrost erosion are unable to explain the episodic, storm-driven erosion events. Our approach, mechanistically coupling oceanographic predictions with a terrestrial model to capture the thermo-mechanical dynamics of erosion, enables this much needed treatment of transient erosion events. The Arctic Coastal Erosion Model consists of oceanographic and atmospheric boundary conditions that force a coastal terrestrial permafrost environment in Albany (a multi-physics based finite element model). An oceanographic modeling suite (consisting of WAVEWATCH III, Delft3D-FLOW, and Delft3D-WAVE) produced time-dependent surge and run-up boundary conditions for the terrestrial model. In the terrestrial model, a coupling framework unites the mechanical and thermal aspects of erosion. 3D stress/strain fields develop in response to a plasticity model of the permafrost that is controlled by the frozen water content determined by modeling 3D heat conduction and solid-liquid phase change. This modeling approach enables failure from any allowable deformation (block failure, slumping, etc.). Extensive experimental work has underpinned the ACE Model development including field campaigns to measure in situ ocean and erosion processes, strength properties derived from thermally driven geomechanical experiments, as well as extensive physical composition and geochemical analyses. Combined, this work offers the most comprehensive and physically grounded treatment of Arctic coastal erosion available in the literature. The ACE model and experimental results can be used to inform scientific understanding of coastal erosion processes, contribute to estimates of geochemical and sediment land-to-ocean fluxes, and facilitate infrastructure susceptibility assessments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Low-temperature formation of pyridine and (iso)quinoline via neutral–neutral reactions

Aromatic molecules represent fundamental building blocks in prebiotic chemistry and are contemplated as vital precursors to DNA and RNA nitrogen bases. However, despite the identification of some 300 molecules in extraterrestrial environments, the pathways to pyridine (C 5 H 5 N), pyridinyl (C 5 H 4 N•), and (iso)quinoline (C 9 H 7 N) – the simplest represen-tative of mono and bicyclic aromatic molecule carrying nitrogen – are elusive. Here, we afford compelling evidence on the gas-phase formation of methylene amidogen (H 2 CN•) and cyanomethyl (H 2 CCN•) radicals via molecular beam studies and electronic structure calcula-tions. The modeling of the chemistries of Taurus Molecular Cloud (TMC-1) and Titan’s atmosphere contemplates a complex chain of reactions synthesizing pyridine, pyridinyl, and (iso)quinoline from H 2 CN• and H 2 CCN• at levels of up to 75%. Furthermore, this study affords unique entry points to precursors of DNA and RNA nitrogen bases in hydrocarbon-rich extraterrestrial environments thus changing the way we think about the origin of prebiotic molecules in our Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Where does the dust deposited over the Sierra Nevada snow come from?

Abstract. Mineral dust contributes up to one-half of surface aerosol loading in spring over the southwestern United States, posing an environmental challenge that threatens human health and the ecosystem. Using self-organizing map (SOM) analysis with dust deposition and flux data from WRF-Chem and Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), we identify four typical dust transport patterns across the Sierra Nevada, associated with the mesoscale winds, Sierra barrier jet (SBJ), North Pacific High (NPH), and long-range cross-Pacific westerlies, respectively. We find that dust emitted from the Central Valley is persistently transported eastward, while dust from the Mojave Desert and Great Basin influences the Sierra Nevada during mesoscale transport occurring mostly in winter and early spring. Asian dust reaching the mountain range comes either from the west through straight isobars (cross-Pacific transport) or from the north in the presence of the NPH. Extensive dust depositions are found on the west slope of the mountain, contributed by Central Valley emissions and cross-Pacific remote transport. In particular, the SBJ-related transport produces deposition through landfalling atmospheric rivers, whose frequency might increase in a warming climate.

Huang, Huilin (ORCID:0000000273286738)↗
DOE OSTI2020

A Fluid Solver for Studying Torsional Galloping in Solar-Tracking PV Panel Arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon, and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. These simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2021

Using ground-based measurements to recover the spectra of radiation escaping from distant light-pollution sources

ABSTRACT A professional dark sky assessment normally requires the numerical modelling of light propagation from ground-based sources into the atmosphere. The spectral power distribution (SPD) of light escaping from an artificial source is one of the key parameters needed in predicting models to simulate the optical signatures of the night sky as a function of distance. SPD at zero distance is due to a non-trivial superposition of light beams, but does not necessarily coincide with the weighted SPD of individual lights. This is because most photons can be removed by obstacles or can interact with many surfaces before being directed into the atmosphere. We have developed an experimental and theoretical method for obtaining the initial SPD from ground-based spectroscopy of the night sky. The method is applicable to distant sources of light and is demonstrated in field measurements conducted on the light dome of Chicago with a new sky glow spectrometer. We show that the method requires analysing the relative impact of local light sources scattered around the measuring site on the radiance of the light dome. Our theoretical method along with ground-based spectroscopy represents a new approach to characterizing light-pollution sources which are otherwise difficult to study due to lack of information on either light-source inventory or a heterogeneous 3D structure of light-emitting and blocking environments.

Kocifaj, M.↗
DOE OSTI2022

On the inconsistency of [C/Fe] abundances and the fractions of carbon-enhanced metal-poor stars among various stellar surveys

Carbon-enhanced metal-poor (CEMP) stars are a unique resource for Galactic archaeology because they probe the properties of the First Stars, early chemical evolution, and binary interactions at very low metallicity. Comparing the fractions and properties of CEMP stars in different Galactic environments can provide us with unique insights into the formation and evolution of the Milky Way halo and its building blocks. In this work, we investigate whether directly comparing fractions of CEMP stars from different literature samples of very metal-poor [Fe/H] < -2.0 stars is valid. We compiled published CEMP fractions and samples of Galactic halo stars from the past 25 years, and find that they are not all consistent with each other. Focusing on giant stars, we find significant differences between various surveys when comparing their trends of [Fe/H] versus [C/Fe] and their distributions of CEMP stars. To test the role of the analysis pipelines for low-resolution spectroscopic samples, we re-analysed giant stars from various surveys with the sspp and ferre pipelines. We found systematic differences in [C/Fe] of ~0.1-0.4 dex, partly independent of degeneracies with the stellar atmospheric parameters. These systematics are likely due to the different pipeline approaches, different assumptions in the employed synthetic grids, and/or the comparison of different evolutionary phases. We conclude that current biases in (the analysis of) very metal-poor samples limit the conclusions one can draw from comparing different surveys. We provide some recommendations and suggestions that will hopefully aid the community to unlock the full potential of CEMP stars for Galactic archaeology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Extreme Precipitation Over the Southern Slope of the Tibetan Plateau and the Associated Atmospheric Circulation Anomalies

The southern slope of the Tibetan Plateau (SSTP) is one of the rainiest regions in the world where geological hazards caused by extreme precipitation often occur. This study investigates the characteristics and mechanisms of extreme precipitation over SSTP from June to September during 2001–2020 using Global Precipitation Measurement satellite observation. The extreme precipitation days are defined as the days with top 5% of regional-mean daily precipitation over SSTP in this period, which has an average precipitation of 27.2 mm/d. Averaging over the extreme precipitation days, precipitation peaks at an altitude of about 300 m, coinciding with the climatological maximum precipitation, but with a much larger value of 37.2 mm/d than the climatology of 11.9 mm/d. Composite analysis of circulations on extreme days reveals significant circulation anomalies in both the lower and upper troposphere. Specifically, the lower-tropospheric circulations are characterized by significant westerly anomalies over northern India, and the upper-tropospheric circulations are characterized by northerly anomalies over the central Tibetan Plateau, which are statistically independent. The lower-tropospheric westerly anomalies blowing toward SSTP are blocked by the topography, favoring extreme precipitation over SSTP. The upper-tropospheric northerly anomalies, on the other hand, correspond to anomalous northeasterlies north of SSTP in the middle troposphere and southeasterlies to the south in the lower troposphere, and the convergence of these circulation anomalies favors extreme precipitation over SSTP. Lastly, the lower-tropospheric westerly and upper-tropospheric northerly anomalies, respectively, correspond to less precipitation over the South Asian monsoon region and the Tibetan Plateau.

Extreme precipitation, Tibetan Plateau↗
DOE OSTI2022

Incorporating the effects of objects in an approximate model of light transport in scattering media

A computationally efficient radiative transport model is presented that predicts a camera measurement and accounts for the light reflected and blocked by an object in a scattering medium. The model is in good agreement with experimental data acquired at the Sandia National Laboratory Fog Chamber Facility (SNLFC). The model is applicable in computational imaging to detect, localize, and image objects hidden in scattering media. Here, a statistical approach was implemented to study object detection limits in fog.

42 ENGINEERING↗
DOE OSTI2022

Simulating Atmospheric Boundary Layer Turbulence with Nek5000/RS

We present large-eddy-simulation (LES) modeling approaches for the simulation of atmospheric boundary layer turbulence that are of direct relevance to wind energy production. In this report, we study a GABLS benchmark problem using high-order spectral element code Nek5000/RS, which is supported under the DOE’s Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms [1, 2]. We demonstrate our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. The model fidelity and scaling performance of Nek5000/RS on DOE’s leadership computing platforms in comparison to those of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, are discussed in [3].

17 WIND ENERGY↗
DOE OSTI2020

AEROICESTUDY-Colorado State University Ice Spectrometer

The dataset comprises measures, using Colorado State University's Ice Spectrometer (IS, an immersion freezing device with a range from 0 °C down to -26 to -29 °C) of atmospheric ice nucleating particle (INP) concentrations taken at ARM's Southern Great Plains (SGP) Site Guest Instrument Facility (GIF), in Lamont, Oklahoma, during the Aerosol-Ice Formation Closure Pilot Study (AEROICESTUDY). A filter sampler was mounted approximately 6 m above ground level atop the GIF rooftop platform. Single-use filter units open to the atmosphere fitted with pre-cleaned and pre-sterilized, 47-mm diameter Nuclepore polycarbonate filters were used. Filters were typically drawn for either a full day period (9AM-5PM CDT) or broken up into morning (9AM-12PM CDT) and afternoon (1PM-5PM CDT) periods. Additionally, 4 overnight filters were taken (6PM-8AM CDT). The volume of air filtered was typically between 3000 and 8000 standard liters for daily samples. After collection, filters were stored at -20 °C until processed. Initial processing to obtain spectra of INP number concentration active via the immersion freezing mechanism versus temperature was conducted using CSU's IS instrument (McCluskey et al., 2018). For measurement of INPs, collected aerosol particles were re-suspended in 8 mL of 0.1 µm-filtered deionized (DI) water. Aliquots of each suspension, and serial dilutions, were dispensed into trays which were fit into aluminum blocks in the IS. Samples were cooled at 0.33 °C per min and the freezing temperatures of wells recorded automatically. Cumulative INP concentrations were determined by first calculating the INPs per mL of suspension based on Vali (1971) and then converting to concentration per standard liter of air using the proportion of the total liquid sample dispensed and the air sample volumes. Two-tailed, 95% confidence intervals for binomial sampling were based on Equation 2 from Agresti and Coull (1998). Aliquots of suspensions from selected samples were also heat treated (95 °C for 20 min) to denature and deactiv

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Cloud condensation nuclei concentrations derived from the CAMS reanalysis

Determining number concentrations of cloud condensation nuclei (CCN) is one of the first steps in the chain in analysis of cloud droplet formation, the direct microphysical link between aerosols and cloud droplets, and a process key for aerosol–cloud interactions (ACI). However, due to sparse coverage of in situ measurements and difficulties associated with retrievals from satellites, a global exploration of their magnitude, source as well as temporal and spatial distribution cannot be easily obtained. Thus, a better representation of CCN numbers is one of the goals for quantifying ACI processes and achieving uncertainty-reduced estimates of their associated radiative forcing. Here, we introduce a new CCN dataset which is derived based on aerosol mass mixing ratios from the latest Copernicus Atmosphere Monitoring Service reanalysis (CAMSRA) in a diagnostic model that uses CAMSRA aerosol properties and a simplified kappa-Köhler framework suitable for global models. The emitted aerosols in CAMSRA are not only based on input from emission inventories using aerosol observations, they also have a strong tie to satellite-retrieved aerosol optical depth (AOD) as this is assimilated as a constraining factor in the reanalysis. Furthermore, the reanalysis interpolates for cases of poor or missing retrievals and thus allows for a full spatiotemporal quantification of CCN numbers. The derived CCN dataset captures the general trend and spatial and temporal distribution of total CCN number concentrations and CCN from different aerosol species. A brief evaluation with ground-based in situ measurements demonstrates the improvement of the modelled CCN over the sole use of AOD as a proxy for CCN as the overall correlation coefficient improved from 0.37 to 0.71. However, we find the modelled CCN from CAMSRA to be generally high biased and find a particular erroneous overestimation at one heavily polluted site which emphasises the need for further validation. The CCN dataset (https://doi.org/10.26050/WDCC/QUAERERE_CCNCAMS_v1, Block, 2023), which is now freely available to users, features 3-D CCN number concentrations of global coverage for various supersaturations and aerosol species covering the years 2003–2021 with daily frequency. This dataset is one of its kind as it offers lots of opportunities to be used for evaluation in models and in ACI studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Wood Decomposition: Understanding Processes Regulating Carbon Transfer to Soil Carbon Pools Using FACE Wood at Multiple Scales

Dead wood comprises approximately 20% of forest biomass. While it is widely recognized for its role in ecosystem functions relative to habitat and nutrient cycling, there are considerable uncertainties associated with its role in the forest carbon (C) cycle. The principal uncertainty is whether dead wood contributes to the soil C pool. With little supporting evidence, the presumption is that C in wood is respired back to the atmosphere through decomposition processes. The inability to attribute forest soil C to its source components has been a barrier to addressing this important forest carbon cycle question. The U.S. Dept. of Energy’s Free Air CO 2 Enrichment (FACE) experiment produced wood that has a distinct δ 13 C signature that can be used to trace its fate through organic matter decomposition processes. Using two hardwood species (aspen, birch) from the Rhinelander FACE site and pine from the Duke FACE site, the FACE Wood Decomposition Experiment (FWDE) was established in 2011 to provide the foundation for addressing wood decomposition and its associated processes at the continental scale. Plots were installed on nine experimental forests representing a range in bioclimatic forest conditions within the continental U.S. Each identical plot has logs of the three species incubating on the soil surface and pine logs suspended in the air, representing standing dead wood. Abiotic conditions are monitored continuously, and soil and forest conditions characterized. Logs on the sites have been measured every two years to assess decay. This project was initiated in 2016 with specific objectives to: (a) determine whether dead wood C contributes to the soil C pool, (b) assess the fungal community composition mediating the decay, (c) assess the interactions of termites with microbial decay, and (d) develop a mechanistic model to simulate wood decomposition. The principal field sampling was conducted in 2017; correspondingly, the results reported here represent six years of decomposition. Wood decomposition varied from 13 to 75 % loss across the bioclimatic gradient. Sites in the southeastern U.S. coastal plain and piedmont had fastest decomposition rates for each of the species. In contrast sites in the Rocky Mountains had the slowest decay rates, which were associated with cold temperatures. Standing deadwood decomposes at a much slower rate than logs on the soil surface; a response due to differences in wood moisture content. Wood carbon was measured in the mineral soil across each of the sites; accordingly, this is the first study to affirm that a portion of the wood C does become part of the soil C pool. Averaged across the FWDE sites, approximately 1.4% of the soil C in the upper 10 cm of the mineral soil directly beneath the log was derived from the deadwood; within 50 cm of the log, wood C accounted for an average of 2.1%. These results demonstrate that the contribution of deadwood to the soil C pool is a soil volume surrounding the log; and it’s likely a larger volume than what we sampled. We also hypothesized that there would be more wood C in the mineral soil when termites were present. However, that hypothesis wasn’t supported by measurements on the Santee Experimental Forest termite exclusion (SEF-TX) plots under FACE logs that incubated for 6 years. Across all the sites, the majority of the wood C was transferred to the atmosphere. Approximately 7,000 fungal species were sequenced. The fungal communities were different among each of the sites suggesting local colonization; interestingly, within sites the three species also exhibited differences in fungal community composition. White rot was the dominant fungal form in pine and hardwood, which yields less lignin in the wood residue. Brown rot occurred in approximately 20% of the pine and aspen logs, and it did not occur in birch logs. The occurrence of brown rot was only weakly correlated to site factors. We conducted two additional experiments to assess the effects of termites on wood decomposition and the associated interactions with the fungal community. Pine and birch logs were incubated on the Santee Experimental Forest termite exclusion plots in upland and riparian settings. There was no difference in decay between areas where termites were able colonize and plots where they had been excluded. Highly variable foraging by termites was a primary contributor the result. Correspondingly, fungal communities did not vary on logs incubated among the two treatments, but they did reflect site and species differences. The second experiment was conducted on each of the nine FWDE sites where pine and birch wood blocks were incubated for two years with and without a wrapping designed to exclude termites and other arthropods. The effects of termites on wood decomposition could be confirmed on two of the four sites known to have subterranean termites. On the two southeastern sites, termites contributed approximately 8% to mass loss in pine blocks over two years. A major factor clouding the role of termites in wood decomposition is the inability to assess the time, duration, and intensity of the foraging with the logs. Another uncertainty is related to whether a piece of wood will be foraged. A process-based wood decomposition model was developed representing biological processes mediated by fungi, bacteria, termites, beetles, as well as physical and chemical processes. The Coarse Wood Decomposition Model (CWDM) was parameterized based on a thorough review of the literature and it was validated using data from the FWDE. In the initial setup, the model considers softwoods and hardwoods groups, as opposed to individual species. Sensitivity analyses incorporating climate data from 89 locations in North America demonstrated that CWDM was stable from the tropics to the boreal zone. Model calibration and validation showed it to perform very well. The CWDM is poised for additional testing with subsequent FWDE data as well as other data sets, and it is available for incorporation into forest biogeochemical models or DOEs Earth System Modeling framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Computational Modeling of Graphite Degradation due to Molten Salt Infiltration and Wear

Molten-salt reactors (MSRs) represent a promising next-generation reactor design, with graphite serving as a moderator and/or reflector in several designs. However, due to limited experimental data and operational experience, a technical understanding of the structural integrity of graphite in molten salt environments remains incomplete. This report presents a modeling-based evaluation of graphite degradation in MSR environments, focusing on the effects of salt infiltration in fuel salt-based designs and surface wear in pebble bed reactor designs. The objective of this study is to enhance understanding of the structural integrity challenges posed by these degradation mechanisms and to provide a framework for assessing graphite behavior in MSRs. The first part of the report investigates the phenomenon of molten salt infiltration into graphite. This infiltration occurs when molten salt permeates the interconnected pore structure of the graphite moderator, driven by factors such as pressure differentials and the physical properties of both the salt and graphite. The infiltration process is influenced by characteristics of the pore structure, viscosity of the molten salt, and the interfacial energies between the graphite, salt, and the atmosphere within the graphite pore. Utilizing a coupled multiphysics modeling approach with Grizzly software, the study evaluates the stress induced by internal heat sources due to infiltration, which can lead to structural concerns. This evaluation is crucial for understanding how infiltration affects the mechanical integrity of graphite components in MSRs. The study considers the Molten-Salt Reactor Experiment (MSRE) graphite stringer geometry due to the availability of relevant data. Through detailed finite element analysis, the study examines stress distributions at varying infiltration percentages, revealing that stress levels increase with higher amounts of infiltration. Rare-event simulations, using the parallel subset simulation (PSS) framework, further quantify the failure probabilities under input uncertainties, with a user-specified failure metric. The PSS framework also identifies critical input parameters that significantly affect the stress values, including infiltration amount, thermal conductivity, and power density. Additionally, considering realistic reactor scenarios, the analysis was performed to account for the combined effects of radiation and infiltration, and modeling strategies on how to analyze new reactor designs or new graphite grades are discussed. The second part of the report focuses on wear mechanisms in pebble bed-based MSRs. As graphite fuel pebbles interact with the graphite reflector block, wear can result in material loss and the formation of surface defects, which may act as stress concentrators. A similar multiphysics modeling framework is employed to assess the impact of wear on the structural integrity of graphite components. This study considers a generic fluoride-cooled high-temperature reactor (gFHR) design due to the availability of comprehensive data. Worst-case scenario dimensions of the reflector blocks were analyzed under thermal and radiation conditions. Subsequently, wear in the form of idealized pits and grooves is modeled on the inner surface of the graphite block, with the maximum stress from previous simulations. The simulations show that groove-type defects are more detrimental than pits, leading to higher stress concentrations. Considering worst-case simulation scenarios and experimental wear rates, it was determined that the formation of a surface defect critical enough to affect the stress may not be possible in a gFHR design. Overall, the findings of this research contribute to the development of robust modeling tools for predicting graphite behavior under various operational conditions in MSRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗