SEARCH · Engineering Papers

Results for “balance function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Adaptive Critical Balance and Firehose Instability in an Expanding, Turbulent, Collisionless Plasma

Abstract Using a hybrid-kinetic particle-in-cell simulation, we study the evolution of an expanding, collisionless, magnetized plasma in which strong Alfvénic turbulence is persistently driven. Temperature anisotropy generated adiabatically by the plasma expansion (and consequent decrease in the mean magnetic-field strength) gradually reduces the effective elasticity of the field lines, causing reductions in the linear frequency and residual energy of the Alfvénic fluctuations. In response, these fluctuations modify their interactions and spatial anisotropy to maintain a scale-by-scale “critical balance” between their characteristic linear and nonlinear frequencies. Eventually the plasma becomes unstable to kinetic firehose instabilities, which excite rapidly growing magnetic fluctuations at ion-Larmor scales. The consequent pitch-angle scattering of particles maintains the temperature anisotropy near marginal stability, even as the turbulent plasma continues to expand. The resulting evolution of parallel and perpendicular temperatures does not satisfy double-adiabatic conservation laws, but is described accurately by a simple model that includes anomalous scattering. Our results have implications for understanding the complex interplay between macro- and microscale physics in various hot, dilute, astrophysical plasmas, and offer predictions concerning power spectra, residual energy, ion-Larmor-scale spectral breaks, and non-Maxwellian features in ion distribution functions that may be tested by measurements taken in high-beta regions of the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Cost- and Risk-Based Seismic Design Optimization of Nuclear Power Plant Safety Systems

Seismic analysis, design, and qualification of systems, structures, and components (SSCs) is a significant contributor to the capital cost of a nuclear power plant. To reduce capital costs of advanced nuclear power plants and make commercial nuclear energy more competitive, innovations are needed in their structural design and construction, and not just in the reactor core and associated systems. Seismic isolation has been identified as an important cost-cutting technology that enables standardization of equipment across various sites. This paper develops and demonstrates a cost- and risk-based seismic design optimization of a representative safety system in a nuclear power plant with the dual goals of minimizing overnight capital cost and meeting safety goals. The design optimization can also include component seismic isolation, in which case, the optimized design includes a set of equipment that needs to be seismically isolated to minimize capital cost. The open-source codes MASTODON and Dakota are used for seismic probabilistic risk assessment and design optimization, respectively. A generic nuclear facility with a safety system comprising SSCs that are common to nuclear power plants is considered for the demonstration of the design optimization and is assumed to be located at the Idaho National Laboratory site. Generic costs and seismic design cost functions are assumed for the SSCs of the safety system. The sum of the costs of the SSCs is minimized in the optimization process, while the risk of failure of the safety system is provided as a constraint. Furthermore, results show that the optimization process reduces capital costs significantly while automatically prioritizing the safety of SSCs that contribute most to the risk of the safety system.

42 ENGINEERING↗

DOE OSTI2023

Yb Substitution and Ultralow Thermal Conductivity of the Ca 3– x Yb x AlSb 3 (0 ≤ x ≤ 0.81(1)) System

Here, a series of Yb-substituted Zintl phases in the Ca 3–x Yb x AlSb 3 (0 ≤ x ≤ 0.81(1)) system has been synthesized by initial arc melting and post-heat treatment, and their isotypic crystal structures were characterized by both powder and single crystal X-ray diffraction analysis. All four title compounds adopted the Ca 3 AlAs 3 -type structure (space group Pnma, Pearson code oP28, Z = 4). The overall structure can be described as a combination of the 1-dimensional (1D) infinite chain of ∞ 1 [Al(Sb 2 Sb 2/2 )] formed by two vertices sharing [AlSb 4 ] tetrahedral moieties and three Ca 2+ /Yb 2+ mixed sites located in between these 1D chains. The charge balance and the resultant independency of the 1D chains in the title system were explained by the Zintl-Klemm formalism [Ca 2+ /Yb 2+ ] 3 [(4b-Al 1– )(1b-Sb 2– ) 2 (2b-Sb 1– ) 2/2 ]. A series of DFT calculations proved that (1) the band overlap between the d-orbital states from two types of cations and the p-orbital states from Sb at the high symmetry Γ point implied a heavily doped degenerate semiconducting behavior of the quaternary Ca 2 YbAlSb 3 model and (2) the site preference of Yb for the M1 site was due to the electronic-factor criterion based on the Q values of each atomic site. The electron localization function calculations also proved that the two different shapes of lone pairs of the Sb atoms—the “umbrella-shape” and the “C-shape”—are determined by local geometry and the coordination environment on the anionic frameworks. Thermoelectric measurements of the quaternary title compound Ca 2.19(1) Yb 0.81 AlSb 3 showed an approximately two times larger ZT value than that of ternary Ca 3 AlSb 3 at 623 K due to increased electrical conductivity and ultralow thermal conductivity originated from Yb substitution for Ca.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Open Architecture for Cost Savings in Advanced Nuclear Reactors

Recently, nuclear power plant build projects in the West have run over budget due to high capital costs and schedule overruns. Compared to other sources of energy, nuclear power plants have higher capital costs. Reactors are often different at every site, resulting in a lack of standardization. Nuclear is expected to compete with other low carbon sources of energy which have lower capital costs making it essential for nuclear to develop ways of reducing costs. Strategies such as standardization, learning rates, modularization, and schedule reduction in advanced reactors can reduce nuclear costs by about 40%. Standardization as a way of cutting capital costs has been explored even in large nuclear power plants. Standardization of certain plant components can result in lower component and installation costs and higher learning from experience. Standardization can be achieved by adopting a criterion of key performance indicators and general design principles for a specific system or component such as the balance of plant. Modularization allows the construction of certain components of SMRs in a factory, which saves time, increases productivity, and encourages higher learning rates. Production learning decreases the time and the cost related to an activity. The potential for modularized components of advanced reactors to be manufactured in factories makes it conducive to achieving higher learning rates. Developing large-capacity nuclear programs through sequential builds cultivates a higher learning rate, which in effect may reduce schedule overruns. Open architecture has been identified as a way to drive standardization among advanced reactor designs and result in cost savings. Open architecture (OA) is defined as a design enabling a diverse supply chain by defining and publishing requirements of systems or equipment in functional and/or interface terms, utilizing technical standards in widespread use. Currently, the nuclear industry’s approach is to use closed architecture, making most designs proprietary. However, collaboration between various advanced reactor vendors and suppliers utilizing the concept of open architecture can result in modular and standardized architecture of subsystems or subcomponents of a nuclear power plant. Completely standardizing nuclear power plants may be impossible, however, certain common subsystems amongst the various reactor designs could be standardized and/or access a wider supply chain and leverage existing learning from other sectors. Open architecture will save time and allocate resources to the parts of the plants that have the most unique features. A key advantage of open architecture is its ability to improve production learning across advanced reactors (AR) types in the industry, by providing and utilizing the same kind of component. Sodium fast reactor (SFR), High Temperature Gas Reactor (HTGR) and Molten Salt Reactor (MSR) are the advanced reactors considered for this project. This paper aims to determine the cost savings in advanced reactor programs due to open architecture learning rate. This work is an extension of work done on light water reactor small modular reactors; the cost methodology was utilized to investigate the impact of open architecture on advanced reactors with a particular focus on sodium fast reactors. The cost data on sodium fast reactors used in the model presented the most adequate information required for the analysis.

Advanced Nuclear Reactors↗

DOE OSTI2025

Engineering PdAu/CeO 2 Alloy/Oxide Interfaces for Selective Methane‐to‐Methanol Conversion with Water

The direct conversion of methane-to-methanol remains a critical challenge in methane valorization. In this study, we unveil the crucial role of PdAu/CeO 2 catalysts in enabling selective methane transformation under mild conditions, using only water as the sole oxidant. Through a combination of experimental techniques, including XPS and catalytic testing, alongside density functional theory (DFT) calculations, we demonstrate that a Pd 0.3 Au 0.7 /CeO 2 catalyst, which predominantly exposes isolated Pd atoms, achieves remarkable methanol selectivity (∼80%) at 500 K with a 1:1 methane-to-water ratio. While Pd/CeO 2 efficiently activates methane, its tendency for overreaction leads to complete methanol decomposition, thereby limiting selectivity. Alloying Pd with Au on ceria mitigates this over-reactivity, preventing methanol degradation while maintaining sufficient catalytic activity. The PdAu/CeO 2 composite exhibits a synergistic effect: Pd in contact with the ceria support facilitates methane activation and water dissociation, while Au fine-tunes reactivity to promote methanol formation. DFT calculations confirm that isolated Pd sites at the PdAu/CeO 2 interface play a key role in balancing activity and selectivity. This work underscores the importance of alloy/oxide interfaces in controlling selective methane conversion with water and offers valuable insights for designing highly efficient catalysts for methanol synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Supramolecular Interactions and Morphology of Self-Assembling Peptide Amphiphile Nanostructures

The morphology of supramolecular peptide nanostructures is difficult to predict given their complex energy landscapes. We investigated peptide amphiphiles containing β-sheet forming domains that form twisted nanoribbons in water. We explained the morphology based on a balance between the energetically favorable packing of molecules in the center of the nanostructures, the unfavorable packing at the edges, and the deformations due to packing of twisted β-sheets. We find that morphological polydispersity of PA nanostructures is determined by peptide sequences, and the twisting of their internal β-sheets. Here, we also observed a change in the supramolecular chirality of the nanostructures as the peptide sequence was modified, although only amino acids with l-configuration were used. Upon increasing charge repulsion between molecules, we observed a change in morphology to long cylinders and then rodlike fragments and spherical micelles. Understanding the self-assembly mechanisms of peptide amphiphiles into nanostructures should be useful to optimize their well-known functions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Accelerated Modeling of Lithium Diffusion in Solid State Electrolytes using Artificial Neural Networks

Abstract Previous efforts to understand structure‐function relationships in high ionic conductivity materials for solid state batteries have predominantly relied on density functional theory (DFT‐) based ab initio molecular dynamics (MD). Such simulations, however, are computationally demanding and cannot be reasonably applied to large systems containing more than a hundred atoms. Here, an artificial neural network (ANN) is trained to accelerate the calculation of high accuracy atomic forces and energies used during such MD simulations. After carefully training a robust ANN for four and five element systems, nearly identical lithium ion diffusivities are obtained for Li 10 GeP 2 S 12 (LGPS) when benchmarking the ANN‐MD results with DFT‐MD. Applying the ANN‐MD approach, the effect of chlorine doping on the lithium diffusivity is calculated in an LGPS‐like structure and it is found that a dopant concentration of 1.3% maximizes ionic conductivity. The optimal concentration balances the competing consequences of effective atomic radii and dielectric constants on lithium diffusion and agrees with the experimental composition. Performing simulations at the resolution necessary to model experimentally relevant and optimal concentrations would be infeasible with traditional DFT‐MD. Systems that require a large number of simulated atoms can be studied more efficiently while maintaining high accuracy with the proposed ANN‐MD framework.

Rao, Karun K.↗

DOE OSTI2023

Machine Learning Prediction of the Experimental Transition Temperature of Fe(II) Spin-Crossover Complexes

Spin-crossover (SCO) complexes are materials that exhibit changes in the spin state in response to external stimuli, with potential applications in molecular electronics. It is challenging to know a priori how to design ligands to achieve the delicate balance of entropic and enthalpic contributions needed to tailor a transition temperature close to room temperature. Here, we leverage the SCO complexes from the previously curated SCO-95 data set [Vennelakanti et al. J. Chem. Phys. 159, 024120 (2023)] to train three machine learning (ML) models for transition temperature (T 1/2 ) prediction using graph-based revised autocorrelations as features. We perform feature selection using random forest-ranked recursive feature addition (RF-RFA) to identify the features essential to model transferability. Of the ML models considered, the full feature set RF and recursive feature addition RF models perform best, achieving moderate correlation to experimental T 1/2 values. We then compare ML T 1/2 predictions to those from three previously identified best-performing density functional approximations (DFAs) which accurately predict SCO behavior across SCO-95, finding that the ML models predict T 1/2 more accurately than the best-performing DFAs. In addition, we study ML model predictions for a set of 18 SCO complexes for which only estimated T 1/2 values are available. Upon excluding outliers from this set, the RF-RFA RF model shows a strong correlation to estimated T 1/2 values with a Pearson’s r of 0.82. In contrast, DFA-predicted T 1/2 values have large errors and show no correlation to estimated T 1/2 values over the same set of complexes. Overall, our study demonstrates slightly superior performance of ML models in comparison with some of the best-performing DFAs, and we expect ML models to improve further as larger data sets of SCO complexes are curated and become available for model training.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Analytical model for quasi-linear flow response to resonant magnetic perturbation in resistive-inertial and viscous-resistive regimes

In this work, a quasi-linear model for plasma flow response to the resonant magnetic perturbation (RMP) in a tokamak has been rigorously developed in the resistive-inertial (RI) and viscousresistive (VR) regimes purely from the two-field reduced MHD model. Models for plasma response to RMP are commonly composed of equations for the resonant magnetic field response (i.e. the magnetic island) and the torque balance of plasma flow. However, in previous plasma response models, the magnetic island and the torque balance equations are often derived separately from reduced MHD and full MHD equations, respectively. By contrast, in this work we derive both the magnetic island response and the torque balance equations in a quasi-linear model for plasma flow response entirely from a set of two-field reduced MHD equations. Such a quasi-linear model can recover previous plasma flow response models within certain limits and approximations. Furthermore, the physical origins of quasi-linear forces and moments in the flow response equation are also accurately calculated and clarified self-consistently.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Evolutionary constraints and climate variability jointly shape starch–sugar balance in woody plants

Nonstructural carbohydrates (NSC) buffer plants against carbon imbalances, yet their partitioning between storage and soluble pools remains elusive at global scales. Here, we compiled a dataset of starch to soluble sugar ratio (St : Su) for 308 woody species across 220 sites world-wide and introduce a dimensionless index that integrates storage and demand while minimizing methodological artifacts. St : Su was strongly associated with growth, identifying it as a key axis of carbon allocation. Foliage consistently exhibited lower St : Su than lignified organs, reflecting a division between transient and conservative pools. Conifers accumulated more starch in foliage but less in stems relative to angiosperms, while leaf habits and mycorrhizal associations further modulated organ-specific strategies. Contrary to expectation, foliar and root St : Su varied little among biomes, but stems exhibited higher ratios in tropical rainforests than in boreal or arid regions, reflecting differences in species composition and adaptive storage under disturbance. Phylogeny constrained stem storage, whereas climatic variability, rather than mean conditions, dominated allocation in leaves and roots. These findings establish St : Su as a robust functional trait linking allocation strategies, growth, and resilience, which can be used to improve vegetation model prediction of forest productivity and mortality under climate variability.

Li, Weibin [Lanzhou Univ. (China)] (ORCID:00000001↗

DOE OSTI2025

Shrub Expansion Can Counteract Carbon Losses From Warming Tundra

Arctic warming is causing substantial compositional, structural, and functional changes in tundra vegetation including shrub and tree-line expansion and densification. However, predicting the carbon trajectories of the changing Arctic is challenging due to interacting feedbacks between vegetation composition and structure, and surface characteristics. We conduct a sensitivity analysis of the current-date to 2100 projected surface energy fluxes, soil carbon pools, and CO 2 fluxes to different shrub expansion rates under future emission scenarios (intermediate—RCP4.5, and high—RCP8.5) using the Arctic-focused configuration of E3SM Land Model (ELM). We focus on Trail Valley Creek (TVC), an upland tundra site in the western Canadian Arctic, which is experiencing shrub densification and expansion. We find that shrub expansion did not significantly alter the modeled surface energy and water budgets. However, the carbon balance was sensitive to shrub expansion, which drove higher rates of carbon sequestration as a consequence of higher shrubification rates. Thus, at low shrub expansion rates, the site would become a carbon source, especially under RCP8.5, due to higher temperatures, which deepen the active layer and enhance soil respiration. At higher shrub expansion rates, TVC would become a net CO 2 sink under both Representative Concentration Pathway scenarios due to higher shrub productivity outweighing temperature-driven respiration increase. Our simulations highlight the effect of shrub expansion on Arctic ecosystem carbon fluxes and stocks. We predict that at TVC, shrubification rate would interact with climate change intensity to determine whether the site would become a carbon sink or source under projected future climate.

Yazbeck, Theresia [The Ohio State Univ., Columbus,↗

DOE OSTI2025

Fouling behavior of zwitterionic membranes compared to polyamide membranes

Membrane fouling remains a critical bottleneck for reverse osmosis (RO) desalination, driving energy consumption and reducing membrane lifetime. Here, we employ all-atom molecular dynamics simulations to investigate the antifouling behavior of random zwitterionic amphiphilic copolymer (r-ZAC) membranes composed of sulfobetaine methacrylate (SBMA) and allyl methacrylate (AMA), benchmarked against conventional polyamide (PA) RO membranes. Structural and dynamical analyses—including radial distribution functions, coordination numbers, tetrahedral order parameters, vector orientation, and residence-time correlation functions—reveal that r-ZAC surfaces sustain tightly bound, long-lived hydration layers with preserved tetrahedrality and anisotropic water orientation, in sharp contrast to the weak and disordered hydration of PA. Steered molecular dynamics simulations demonstrate that r-ZAC membranes impose substantial free-energy barriers to foulant approach (alginate ≈ 90 kcal/mol, sucrose ≈ 35 kcal/mol, humic acid ≈ 15 kcal/mol), whereas PA membranes exhibit negligible barriers (< 1 kcal/mol) and thermodynamically favorable adsorption. Detailed foulant–surface interaction analyses show that zwitterionic hydration and electrostatic heterogeneity in r-ZAC suppress adhesion, except in the case of amphiphilic humic acid, which exploits multiple binding modes. Together, these results establish molecular-level design principles for antifouling membranes: the combination of zwitterionic hydration, structured interfacial water, and controlled amphiphilic balance in r-ZAC membranes provides superior resistance to organic fouling relative to PA.

Cross-linked polyamide↗

DOE OSTI2024

Predicting Thermochemical Equilibria with Interacting Defects: Sr 1 − x Ce x Mn O 3 − δ Alloys for Water Splitting

Solar thermochemical hydrogen is one of the few potential routes towards direct fuel production from renewable energy sources, but the thermodynamic boundary conditions for efficient and economic energy conversion are challenging. Success or failure of a given oxide working material depends on the subtle balance between enthalpy and entropy contributions in the redox processes. Developing a mechanistic understanding of the behavior of materials on the basis of atomistic models and first-principles calculations is an important part of advancing the technology. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox-active defects start to interact with each other, thereby impeding the formation of additional defects. This problem is of more general importance to applications that rely on high levels of off-stoichiometry or doping, including, for example, batteries, thermoelectrics, and ceramic fuel cells. To account for such repulsive defect interactions, we introduce a statistical mechanics approach, defining an expression for the free energy of defect interaction based on limited sampling of defect configurations in density functional theory supercell calculations. The parameterization of this energy contribution as a function of defect concentration and temperature allows on-the-fly simulation of thermochemical equilibria. The approach consistently incorporates finite temperature effects by including the leading contributions to the temperature-dependent free energy for the case at hand, i.e., the ideal gas and configurational enthalpies and entropies. We demonstrate the capability and utility of the approach by simulating the water splitting redox processes for Sr 1 − x Ce x Mn O 3 − δ alloys. Published by the American Physical Society 2024

08 HYDROGEN↗

DOE OSTI2022

Combining GWAS and population genomic analyses to characterize coevolution in a legume‐rhizobia symbiosis

Abstract The mutualism between legumes and rhizobia is clearly the product of past coevolution. However, the nature of ongoing evolution between these partners is less clear. To characterize the nature of recent coevolution between legumes and rhizobia, we used population genomic analysis to characterize selection on functionally annotated symbiosis genes as well as on symbiosis gene candidates identified through a two‐species association analysis. For the association analysis, we inoculated each of 202 accessions of the legume host Medicago truncatula with a community of 88 Sinorhizobia (Ensifer) meliloti strains. Multistrain inoculation, which better reflects the ecological reality of rhizobial selection in nature than single‐strain inoculation, allows strains to compete for nodulation opportunities and host resources and for hosts to preferentially form nodules and provide resources to some strains. We found extensive host by symbiont, that is, genotype‐by‐genotype, effects on rhizobial fitness and some annotated rhizobial genes bear signatures of recent positive selection. However, neither genes responsible for this variation nor annotated host symbiosis genes are enriched for signatures of either positive or balancing selection. This result suggests that stabilizing selection dominates selection acting on symbiotic traits and that variation in these traits is under mutation‐selection balance. Consistent with the lack of positive selection acting on host genes, we found that among‐host variation in growth was similar whether plants were grown with rhizobia or N‐fertilizer, suggesting that the symbiosis may not be a major driver of variation in plant growth in multistrain contexts.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Spread spectrum time domain reflectometry (SSTDR) and frequency domain reflectometry (FDR) cable inspection using machine learning

Cables are initially qualified for nuclear power plant use for 40 years. As plants extend their operating license to 60 and 80 years, justification for continued cable use must shift to a condition-based approach since it is cost prohibitive to completely replace cables that are likely still capable of performing their design function. The Pacific Northwest National Laboratory (PNNL) Accelerated and Real Time Experimental Nodal Analysis (ARENA) cable motor test bed was used to test the response of a commercial spread spectrum time domain reflectometry (SSTDR) system, a laboratory instrument software-controlled SSTDR, and a vector network analyzer-based frequency domain reflectometry (FDR) system to various cable anomalies. The three instrument systems were able to interrogate cables over a range of frequency bandwidths that can be helpful for human data analysis. Data were subjected to supervised and unsupervised machine learning (ML) analyses to distinguish normal undamaged cable responses from anomalous cable responses. Both supervised and unsupervised ML approaches produced encouraging results with an undamaged/anomalous prediction accuracy from 0.69% to 0.87%. Recommendations for further development and field implementation include increased and more balanced sample sets particularly including more training data.

SSTDR, FDR, Reflectometry, Machine Learning, ARENA↗

DOE OSTI2022

A Mountaintop View Requires Minimal Sorting: A Faster Contour Tree Algorithm

Consider a scalar field f : M → R, where M is a triangulated simplicial mesh in R d . A level set, or contour, at value v is a connected component of f –1 (v). As v is changed, these contours change topology, merge into each other, or split. Contour trees are concise representations of f that track this contour behavior. The vertices of these trees are the critical points of f, where the gradient is zero. The edges represent changes in the topology of contours. It is a fundamental data structure in data analysis and visualization, and there is significant previous work (both theoretical and practical) on algorithms for constructing contour trees. Suppose M has n vertices, N facets, and t critical points. A classic result of Carr, Snoeyink, and Axen (2000) gives an algorithm that takes O(n log n+Nα(N)) time (where α(·) is the inverse Ackermann function). A further improvement to O(t log t + N) time was given by Chiang et al. All these algorithms involve a global sort of the critical points, a significant computational bottleneck. Unfortunately, lower bounds of Ω(t log t) also exist. We present the first algorithm that can avoid the global sort and has a refined time complexity that depends on the contour tree structure. Intuitively, if the tree is short and fat, we get significant improvements in running time. For a partition of the contour tree into a set of descending paths, P, our algorithm runs in O($\Sigma$ pϵP |p| log |p| + tα(t) + N). This is at most O(t log D + N), where D is the diameter of the contour tree. Moreover, it is O(tα(t) + N) for balanced trees, a significant improvement over the previous complexity. Our algorithm requires numerous ideas: partitioning the contour tree into join and split trees, a local growing procedure to iteratively build contour trees, and the use of heavy path decompositions for the time complexity analysis. There is a crucial use of a family of binomial heaps to maintain priorities, ensuring that any comparison made is between comparable nodes of the contour tree. We also prove lower bounds showing that the $\Sigma$ pϵP |p| log |p| complexity is inherent to computing contour trees.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Climate adaptation and sustainability in switchgrass: exploring plant-microbe-soil interactions across continental scale environmental gradients

Less carbon-intensive energy sources are needed to reduce greenhouse gas emissions and their predicted role in climate change. There is growing interest in the potential of biofuels for meeting this need. A critical question is whether large-scale biofuel production can be sustainable over the time scales needed to mitigate our carbon debt from fossil fuel consumption. The carbon balance and ultimately the sustainability of biofuel feedstock production is the result of complex climate-coupled interactions between carbon fixation, sequestration, and release through combustion. Similarly, the long-term productivity of biofuels depends on the environmental factors limiting plant growth. These factors are often related to soil resources which involve complex interactions at the plant-microbe-soil interface impacting their availability and cycling. Our collaborative project addressed sustainable switchgrass (Panicum virgatum) production by exploring Plant Systems, Plant-Microbiome Interactions, and Ecosystem Processes through the integrating lens of Multi-Scale Modeling. Our research was based on detailed characterization of genetically diverse switchgrass genotypes planted in common gardens across a continental latitudinal gradient. The underlying theme of our Plant Systems research was the use of locally adapted plant material to explore plant function, to understand the mechanistic basis of environmental interactions, and to discover the plant genes important for adaptation and sustainability in the face of climate change. Our Plant-Microbiome Interaction project characterized the microbial communities associated with switchgrass using genomic tools. Our Ecosystem Processes research focused on carbon cycle responses at the ecosystem level using stand level plantings. Finally, our Multi-Scale Modeling helped to define conditions of a sustainable biofuel system and identify key tradeoffs between genetic diversity, productivity, and ecosystem services. Genome-wide association analyses were used to identify alleles that contribute to successful establishment and biomass production across North America. Together, our work provided a baseline analyses of the potential of switchgrass as a biofuel feedstock. Our project resulted in a number of successful outcomes. First, we were successful in collecting switchgrass germplasm across the species range, propagating the material, and establishing common garden experiments across the species range. In collaboration with DOE JGI, we successfully assembled the first tetraploid switchgrass genome and published this resource with an analyses of the genetic basis local adaptation from our gardens (Lowry et al. 2019, Lovell et al. 2021). The gardens were used to characterize the genetic architecture for a number of important plant phenotypes. Our project also conducted extensive sampling and sequencing to characterize the bacterial and fungal associates of switchgrass roots and leaves. We showed that host genotype, location, and harvesting practices can play a role in microbiome assembly (Singer et al. 2019 & 2022, Van Wallendael et al. 2020 & 2022, Edwards et al. 2023). Our ecosystem processes work created baseline dataset of carbon and nutrient cycling in realistic stand plantings of switchgrass. Data from this experiment provided new insight into the role of plant traits, phenology, and local environments in ecosystem processes like soil respiration, net-ecosystem exchange, and dynamics of soil and plant nutrients (Ricketts et al. 2023). Finally, our crop modelling experiments help to characterize the sensitivity of common modeling frameworks to parameters, identify key limiters of productivity across large geographic scales, and leverage patterns of local adaptation in prediction. Ultimately, these studies help to identify critical plant-microbe-soil traits that may be manipulated, through breeding or agronomic management, to improve the sustainability of biofuel feedstocks.

09 BIOMASS FUELS↗

DOE OSTI2022

Evaluated 238 U(n,f) Average Prompt Fission Neutron Multiplicities Including the CGMF Model

This report documents an evaluation of the average prompt fission neutron multiplicity, $\overline{v}_p$, of 238 U from 800 keV to MeV. This evaluation had to be re-done from “scratch” as the input to previous $\overline{v}_p$ evaluations, specifically ENDF/B-VIII.0, was not found. That means that all available experimental data were re-analyzed and uncertainties were re-estimated. The new evaluated 238 U $\overline{v}_p$ based on only experimental data differs distinctly from ENDF/B-VIII.0 $\overline{v}_p$ from 2 to 4.5 MeV, and from 6 to 7 MeV, and is otherwise similar. The difference from 2 to 4.5 MeV stems from the fact that ENDF/B-VIII.0 was tweaked in this energy range to data of Frehaut, while two other, equally trustworthy, data sets would indicate an evaluated 238 U $\overline{v}_p$ that is up to 2% higher. Also, second chance fission in ENDF/B-VIII.0 was smoothed over from 6–7 MeV. Another major difference to ENDF/B-VIII.0 is that one of the evaluations presented here includes model information from the Hauser-Feshbach fission fragment decay code CGMF, while ENDF/B-VIII.0 is based purely on experimental data. CGMF links several fission quantities with each other; $\overline{v}_p$ is predicted by assumptions made on, e.g., pre-neutron emission yields as a function of mass, the total kinetic energy, or spin and parity of fission fragments. This allows to validate the new 238 U $\overline{v}_p$ by using CGMF parameters obtained from fitting to experimental 238 U $\overline{v}_p$ to predict yields as a function of mass, the average total kinetic energy, or the mean energy of the prompt fission neutron spectrum. These model-predicted values can then be compared to experimental and evaluated data. The model-predicted fission-observable values using evaluated parameters obtained here are reasonably close to experimental data for some observables, but are farther away from experimental data related to TKE observables. In addition to that, the evaluated 238 U(n,f) $\overline{v}_p$ shows similar deviations from ENDF/B-VIII.0 as for the evaluation with only experimental data. This difference is expected to lead to changes in simulated effective neutron multiplication factor, $k_{eff}$ of ICSBEP critical assemblies that are sensitive to 238 U in the fast range (BigTen, Flattop, Flattop-Pu). These changes in $k_{eff}$ need to be counter-balanced. Chi-Nu PFNS experimental data are expected to be released in the next few months that might lead to the needed changes in the PFNS. Until then, we hold off in benchmarking the new 238 U(n,f) $\overline{v}_p$ as well as submitting it to ENDF/B-VIII.1. Also, new high-precision 238 U $\overline{v}_p$ are expected to be measured by the CEA in the next two years that will shed further light on question on 238 U $\overline{v}_p$ from 2–4.5 and 6–7 MeV.

238U↗