Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “robust parameter design”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Advancing Concentrating Solar Thermal Modeling Using System Advisor Model (SAM)

Concentrating solar thermal (CST) technologies play a critical role in enabling dispatchable power and high-temperature industrial heat applications. Accurate and flexible modeling tools are essential for evaluating system performance, guiding technology research and development, and informing investment decisions. The National Laboratory of the Rockies's System Advisor Model (SAM) is a widely used techno-economic simulation platform for CST systems, providing detailed performance and financial modeling capabilities for multiple CST system configurations. SAM integrates physics-based performance models with financial analysis to simulate the behavior of complex energy systems under realistic operating conditions. For CST technologies (including tower, parabolic trough, and linear Fresnel), SAM enables hourly simulations using site-specific weather data that ensure feasible operating conditions and convergence of mass and energy between core system components (i.e., solar field, receiver, thermal energy storage, and power cycle). These capabilities allow researchers and developers to evaluate annual energy production, capacity factors, levelized cost of energy (LCOE), and system dispatch strategies. A key advantage of SAM lies in its flexibility for parametric analysis and large-scale computational studies. Users can vary system design parameters such as heliostat field layout, receiver dimensions, thermal energy storage capacity, power block sizing, and installation cost assumptions to investigate their impact on system performance and financial metrics. When combined with automated scripting through LK, SDKTool, or Python interfaces, SAM enables high-throughput simulation workflows that support sensitivity analysis, technology benchmarking, and optimization studies. These approaches are particularly valuable for next-generation CST concepts, where design spaces are large and system interactions are complex. Another important capability of SAM is its support for dispatch optimization and thermal energy storage modeling, which are central to the value proposition of CST technologies. The ability to simulate integrated storage and flexible power generation allows researchers to explore strategies that maximize grid value, improve capacity utilization, and enhance integration with variable resources such as photovoltaic and wind generation. This poster will present an overview of SAM's thermal system modeling capabilities including concentrating solar. Additionally, we will highlight new feature developments including: 1) implementing Google's OR-Tools optimization platform for faster and more robust dispatch optimization, 2) developing a new power load following controller for modeling behind-the-meter applications, 3) enabling direct modeling of CSP-PV hybrid systems with the inclusion of battery storage, and 4) developing a multi-receiver falling particle Gen3 system model.

14 SOLAR ENERGY↗
DOE OSTI2021

BETO 2021 Peer Review - FCIC Task 6: High Temperature Conversion

The impacts of feedstock variability on pyrolysis processes are significant but poorly defined. Current engineering designs are based on empirical guidelines, useful only over a narrow range of feedstock properties. The objectives of this project are to (1) Develop science-based knowledge of how feedstock attributes and operational parameters impact pyrolysis process reliability and product quality; and (2) Build an experimental and computational toolset that predicts these outcomes, enabling processes to optimize reliability and product quality. Biomass is a complex feedstock. Controlling for and testing the effects of individual attributes is very challenging. This project couples multiscale experimentation and modeling to accurately capture the fundamental physics and chemistry of biomass flow and conversion behavior in feeding and pyrolysis reactor operations. Our focus is on pine residue attributes – anatomical fraction (bark, needles, wood), particle morphology (size/shape distribution, density, porosity), and chemical composition (extractives, biopolymers, alkali metals) – that impact product quality for downstream catalytic upgrading. Because detailed pyrolysis product characterization is limited, cutting-edge analytical techniques are being developed to reveal impactful product attributes. The tools and knowledge developed here will enable integrated pyrolysis-based processes that are more robust, flexible, and market-responsive with respect to feedstock variability.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

On the use of vibronic coherence to identify reaction coordinates for ultrafast excited-state dynamics of transition metal-based chromophores

In this study, the question of whether one can use information from quantum coherence as a means of identifying vibrational degrees of freedom that are active along an excited-state reaction coordinate is discussed. Specifically, we are exploring the notion of whether quantum oscillations observed in single-wavelength kinetics data exhibiting coherence dephasing times that are intermediate between that expected for either pure electronic or pure vibrational dephasing are vibronic in nature and therefore may be coupled to electronic state-to-state evolution. In the case of a previously published Fe(ιι) polypyridyl complex, coherences observed subsequent to 1 A 1 → 1 MLCT excitation were linked to large-amplitude motion of a portion of the ligand framework; dephasing times on the order of 200–300 fs suggested that these degrees of freedom could be associated with ultrafast (~100 fs) conversion from the initially formed MLCT excited state to lower-energy, metal-centered ligand-field excited state(s) of the compound. Incorporation of an electronically benign but sterically restrictive Cu(ι) ion into the superstructure designed to interfere with this motion yielded a compound exhibiting a ~25-fold increase in the compound’s MLCT lifetime, a result that was interpreted as confirmation of the initial hypothesis. However, new data acquired on a different chemical system – Cr(acac') 3 (where acac' represents various derivatives of acetylacetonate) – yielded results that call into question this same hypothesis. Coherences observed subsequent to 4 A 2 → 4 T 2 ligand-field excitation on a series of molecules implicated similar vibrational degrees of freedom across the series, but exhibited dephasing times ranging from 340 fs to 2.5 ps without any clear correlation to the dynamics of excited-state evolution in the system. Taken together, the results obtained on both of these chemical platforms suggest that while identification of coherences can indeed point to degrees of freedom that should be considered as candidate modes for defining reaction trajectories, our understanding of the factors that determine the interplay across coherences, dephasing times, and electronic and geometric structure is insufficient at the present time to view this parameter as a robust metric for differentiating active versus spectator modes for ultrafast dynamics.

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

BSMV-mediated genome editing exhibits host-specific heritability: germline transmission in barley and somatic edits in Nicotiana benthamiana

Plant RNA virus–mediated guide RNA (gRNA) delivery represents a transformative advance in genome editing technologies. Unlike conventional transformation methods that rely on labor-intensive tissue culture and regeneration for each individual gRNA delivery, viral vectors can rapidly and systemically transmit gRNAs into pre-established Cas-expressing plants, providing an accelerated route for functional genomics and trait discovery directly in planta . However, key design parameters, including subgenomic promoter choice, transcript architecture, and their effects on viral fitness and editing outcomes, remain to be elucidated for most viral platforms. We developed five Barley stripe mosaic virus (BSMV) vectors, each with distinct subgenomic promoter elements to drive single gRNA expression. These were initially evaluated in Cas9-expressing transgenic Nicotiana benthamiana plants targeting the Phytoene desaturase ( PDS ) gene to compare their editing efficiencies. Single gRNAs expressed under the duplicated γb subgenomic promoter or when fused directly to the γb genome achieved the highest mutation frequencies (up to 90% at 60 days post-inoculation), whereas β1- and β2-driven sgRNAs produced delayed and reduced editing. Thus, promoter selection critically determines gRNA accumulation and the efficacy of BSMV-mediated genome editing. The top-performing design was then applied to Cas9-expressing barley ( Hordeum vulgare ) targeting HvCMF7 (conferring green-white variegation) and HvGW2.1 (impacts grain width and weight). BSMV spread systemically throughout barley, inducing somatic and heritable mutations at frequencies up to 100%, with virus-free edited progeny. In contrast, despite robust somatic editing in N. benthamiana, no heritable mutations were detected indicating species-dependent limitations in germline transmission. Our systematic comparison of subgenomic promoter architectures establishes clear design principles for optimizing viral vector–mediated delivery. Promoter choice and transcript structure critically shape editing efficiency and viral stability. The host-specific boundary for germline editing, defined by efficient heritable editing in barley but not N. benthamiana , highlights where BSMV offers advantages and where alternative vectors or hybrid strategies are required, guiding rational platform selection for diverse crop species and applications. Collectively, these findings establish BSMV as a promising next-generation vector for rapid, tissue culture–free, and transformation-independent genome editing in cereals and other recalcitrant monocots.

barley↗
DOE OSTI2025

Investigating performance and variability of NIF ICF experiments with deep learning

The parameter space involved in designing an inertial confinement fusion shot at the National Ignition Facility (NIF) is massively multi-dimensional and the cost of a single shot makes a comprehensive set of sensitivity studies in the laboratory impractical. The use of machine learning to overcome these challenges has gained popularity and has had several successful applications by the scientific community. We extend on these efforts by training a neural network (NN) on information about the experimental design, engineering elements, and drive asymmetry to predict with uncertainty the neutron yield of an experiment. We find the measured and model predicted values are in good agreement, with an R 2 value of 0.91 for a randomly selected test dataset. Almost all the predicted 95% credible intervals contain the corresponding measured value for both training and test datasets. We identify correlations picked up by the NN between the shot design, yield, and variability and use them to motivate shot sensitivity studies. The first shot to exceed the Lawson-like ignition criteria (N210808) was conducted at the NIF and subsequent shots studied the design’s robustness. In a follow-up shot to N210808, our model predicts capsule quality to be the main performance degradation mechanism that prevented the shot from repeating previous performance levels. Shot N221204 was the first shot to exceed a target energy gain of 1. Our model predicts increased yield with reduced coast time for a N221204 study and greater variability for designs with lower peak powers at constant yield. The model’s fast prediction speed and uncertainty prediction are useful for identifying interesting design paths that could warrant further investigation with conventional simulations to search for robust high yield designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Sensitivity analysis, surrogate modeling, and optimization of pebble-bed reactors considering normal and accident conditions

This research provides a valuable tool that streamlines the optimization process while significantly increasing its accuracy. This study creates a robust framework for reactor design optimization by incorporating comprehensive modeling using the Comprehensive Reactor Analysis Bundle, or BlueCRAB, within the Multiphysics Object-Oriented Simulation Environment (MOOSE). BlueCRAB is the United States Nuclear Regulatory Commission's code suite for non-light water reactor analysis and includes the Griffin, Pronghorn, and Bison applications. This not only improves the efficiency of the optimization process but also enhances the reliability of the results. Such a tool is essential for advancing the state-of-the-art in pebble-bed reactor technology and is critical for achieving the goals of Generation IV reactors, which aim for safe, sustainable, and economically viable nuclear energy solutions. This work presents and applies this workflow on pebble-bed reactors while considering both normal and off-normal conditions. A representative gas-cooled pebble-bed reactor at equilibrium core conditions serves as the nominal design specification for normal operation and is based on previous research. The depressurized loss-of-forced-cooling accident is deployed for off-normal conditions in this work. After defining design-related parameters and quantities of interest regarding reactor safety and performance, this multiphysics model is sampled using the MOOSE stochastic tools module. The result is a comprehensive dataset of configurations, enabling sensitivity analysis and the generation of surrogate models. Subsequently, the dataset and surrogate models are employed in two optimization studies aimed at maximizing fuel utilization and economic profit while adhering to safety and operational constraints. Performing the optimization process with fuel utilization as the metric leads to an improvement of approximately 10%, compared to engineering-judgment-based nominal conditions. The optimization on economic profit leads to an estimated increase of ~300 million USD over the lifetime of the reactor.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Constructing Data-Driven Predictions at the Far Detector for NOvA's Neutrino Oscillation Analysis.

NOvA, is a two-detector, long-baseline neutrino oscillation experiment located at Fermilab, Batavia, IL, USA. It is designed primarily to constrain neutrino oscillation parameters using $\nu_\mu \ (\bar{\nu}_\mu)$ disappearance and $\nu_e \ (\bar{\nu}_e)$ appearance data. The Neutrinos at Main Injector (NuMI) beamline at Fermilab provides a high purity 900 KW intense beam of neutrinos and anti-neutrinos to NOvA. The NOvA Near Detector, located 100m underground and 1km away from the beam source, observes the un-oscillated $\nu_\mu \ (\bar{\nu}_\mu)$ and beam $\nu_e \ (\bar{\nu}_e)$ event spectrum. The Far Detector, located in Ash River, MN, USA, is 809 km from the ND and records the oscillated $\nu_e \ (\bar{\nu}_e)$ and the un-oscillated $\nu_\mu \ (\bar{\nu}_\mu)$ event spectrum. NOvA uses a data-driven technique called extrapolation to predict the expected number of $\nu_\mu \ (\bar{\nu}_\mu)$ and $\nu_e \ (\bar{\nu}_e)$ events at the Far Detector using the Near Detector data. The use of data from a functionally equivalent Near Detector provides a powerful constraint on the systematic uncertainties in NOvA neutrino oscillation analyses. As NOvA continues to add data statistics, a robust constraint on systematics becomes more crucial for neutrino oscillation analysis. The details of the NOvA neutrino oscillation analysis framework and how it constrains dominant systematic uncertainties using the Near Detector data will be discussed in this poster.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Risk assessment of wellbore leakage during underground hydrogen storage

The expansion of renewable energy sources would require large-scale energy storage options to overcome the intermittent nature of these sources. Underground hydrogen storage (UHS) in depleted hydrocarbon reservoirs offers a scalable and practical energy storage solution. These reservoirs are chosen for their availability and large capacity, but the unique properties of hydrogen raise concerns about potential leakage pathways, particularly through wellbores. In this study, we develop and apply, for the first time, reduced-order models (ROMs) specifically designed for efficient leakage risk prediction in UHS systems operating in depleted hydrocarbon reservoirs. Using 3,000 high-fidelity simulation scenarios, we examine the influence of 11 key parameters, including reservoir and aquifer depths, wellbore permeability and porosity, initial saturations of water, oil and gas fractions (hydrogen, light, intermediate, and heavy hydrocarbons), reservoir pressure multiplier, and the aquifer-to-reservoir volume ratio, to simulate leakage behavior over a 1,000-year timescale. We train ROMs using a two-step classification-regression approach, achieving R 2 values exceeding 99 % across all targets. These ROMs effectively capture the leakage evolution and identify critical controls of leakage, guiding the design of mitigation strategies. Results indicate that gas leakage occurs in about 27 % of scenarios as early as five years post-operation, reaching volumes of up to 106 ft3. Oil leakage is less frequent (~17 %) and typically begins decades later. Our findings also show that hydrogen often migrates first, owing to its smaller molecular size and higher buoyancy, followed by heavier hydrocarbons. Over time, these heavier components contribute significantly to the total leaked volume, reinforcing the need for targeted monitoring and remediation strategies. Our analysis highlights that deeper storage reservoirs, shallower aquifers, and low-permeability wellbores significantly reduce leakage risks. In conclusion, this work offers a robust framework for risk-informed UHS deployment, supporting energy security through reliable large-scale hydrogen storage while safeguarding environmental integrity.

08 HYDROGEN↗
DOE OSTI2025

Thin-Film Fracture Behavior for Diketopyrrolopyrrole Semiconducting Polymeric Films

Fracture energy, which quantifies a material’s resistance to the propagation of a pre-existing crack, is a key parameter for ensuring the mechanical reliability of stretchable organic electronic devices. However, most existing methods, such as a four-point bending fracture energy, utilized for measuring the fracture energy of semiconducting polymeric thin films are complicated by substrate effects, making it challenging to isolate the intrinsic behavior of the film from interfacial influences. In this study, we employed a pseudo free-standing pure shear method to systematically investigate the cohesive fracture energy of poly(diketopyrrolopyrrole-terthiophene) P(DPP-T)-based thin films to examine the effects of nanoconfinement, side chain length, degree of crystallinity, and strain rates. This method effectively eliminates substrate interference, enabling a direct assessment of the cohesive fracture energy of P(DPP-T) thin films. We found that thinner films and those with lower molecular weights exhibited significantly reduced fracture energies due to diminished chain entanglements. Additionally, films with shorter side chains displayed notably higher fracture energies, which were attributed to an increase in the degree of crystallinity. Finally, slower strain rates led to higher fracture energies, consistent with an enhanced stress relaxation. These insights offer practical guidelines for designing mechanically robust semiconducting polymers, contributing to the advancement of reliable, durable, flexible, and wearable electronic devices.

Zhang, Song [Univ. of Southern Mississippi, Hattie↗
DOE OSTI2025

Defect-Engineered High-Entropy Spinel Oxide@Onion-Like Carbon Catalysts for High-Areal-Energy Rechargeable Zinc–Air Batteries

Rechargeable zinc-air batteries (ReZAB) have emerged as the next-generation batteries with several advantages over the conventional lithium-ion battery. In this work, single nanocrystals of inverse-type high-entropy spinel oxides (HESOx, particle size of 10−12 nm) confined in highly curved defective onion-like carbons (HESOx/OLC AT ) as efficient electrocatalysts for oxygen evolution reaction (OER), oxygen reduction reaction (ORR), and ReZAB, have been synthesized. The HESOx materials were thoroughly characterized using several analytical techniques, including X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), scanning transmission electron microscopy (STEM), Raman, and electron paramagnetic resonance (EPR). HESOx/OLC AT catalyst was tested for ReZAB using literature-recommended parameters that would allow for real technological application. These parameters include a current loading of 10 mA cm −2 and a discharge areal energy density of 35 mWh cm geometric −2 , which maps a Li-ion battery pack-level specific energy of 120 Wh kg pack −1 . HESOx/OLC AT electrocatalysts allowed for continuous discharging and charging at a current loading of 10 mA cm −2 with discharge areal energy densities between 37 and 74 mWh cm geometric −2 , thus outperforming the recommended threshold of 35 mWh cm geometric −2 . Considering that most studies (>90%) hardly meet the recommended threshold for technological application of ReZAB, the present work represents one of the top-performing electrocatalysts for ReZAB. The excellent electrocatalytic properties of defect-rich HESOx/OLC AT toward ORR/OER and ReZAB are governed by the strong electronic modulation arising from d-π hybridization, the availability of multiple catalytic sites for intermediates, and weakened d-band centers of the rate-determining intermediates (i.e., *O adsorption for ORR and *OOH formation for OER) compared to the pristine HESOx. This work introduces an effective approach for the design and synthesis of single nanocrystals of high-entropy electrocatalysts for the development of low-cost, robust, and technologically relevant rechargeable zinc−air batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

SANS quantification of bound water in water-soluble polymers across multiple concentration regimes

Contrast-variation small-angle neutron scattering (CV-SANS) is a widely used technique for quantifying hydration water in soft matter systems, but it is predominantly applied in the dilute regime or for systems with a well-defined structure factor. Here, CV-SANS was used to quantify the number of hydration water molecules associating with three water-soluble polymers with different critical solution temperatures and types of water–solute interactions in dilute, semidilute, and concentrated solution through the exploration of novel methods of data fitting and analysis. Multiple SANS fitting workflows with varying levels of model assumptions were evaluated and compared to give insight into SANS model selection. These fitting pathways ranged from general, model-free algorithms to more standard form and structure factor fitting. In addition, Monte Carlo bootstrapping was evaluated as a method to estimate parameter uncertainty through simulation of technical replicates. The most robust fitting workflow for dilute solutions was found to be form factor fitting without CV-SANS ( i.e. polymer in 100% D 2 O). For semidilute and concentrated solutions, while the model-free approach can be mathematically defined for CV-SANS data, the addition of a structure factor imposes physical constraints on the optimization problem, suggesting that the optimal fitting pathway should include appropriate form and structure factor models. The measured hydration numbers were consistent with the number of tightly bound water molecules associated with each monomer unit, and the concentration dependence of the hydration number was largely governed by the chemistry-specific interactions between water and polymer. Polymers with weaker water–polymer interactions ( i.e. those with fewer hydration water molecules) were found to have more bound water at higher concentrations than those with stronger water–polymer interactions due to the increase in the number of forced water–polymer contacts in the concentrated system. This SANS-based method to count hydration water molecules can be applied to polymers in any concentration regime, which will lead to improved understanding of water–polymer interactions and their impact on materials design.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Creating multifunctional synthetic lichen platforms for sustainable biosynthesis of biofuel precursors

In this project, we were creating a sustainable platform for biofuel production, utilizing carbon-fixing autotrophs to supply oxygen and organic substrates to heterotrophic partners, which in turn produce carbon dioxide to feed the autotrophs. This symbiotic lichen community could lower the input cost, optimize metabolic exchanges and improve the generation of biofuel precursors through multi-omics driven genetic engineering. The cyanobacteria Synechococcus elongatus (S. elongatus) was used as the primary autotroph to provide oxygen and organic substrates, especially sucrose, to a co-culture system. The strain with overexpression of sucrose transporter cscB demonstrated a significant increase in sucrose production under salt stress as what we expected. We also implemented 13C metabolic flux analysis on the sucrose secreting strain S. elongatus cscB-NaCl. Next, transporters proteins like glutamate exporter mscCG from Corynebacterium glutamicum was overexpressed in S. elongatus to improve metabolite exchange. We provided sucrose and glutamate to filamentous fungi to enable their growth and production of biochemicals using substrates from the cyanobacterium. Two fungi ( Aspergillus nidulans and Aspergillus niger ) and two yeast strains ( Rhodotorula toruloides and Lipomyces starkeyi ) served as heterotrophs in this system. They were co-cultured with S.elongatus under different pH condition, and growth on different carbon source to understand the symbiotic system and optimize the parameters for production. The two fungi strained grew much better at pH 7, 8 and 9 with Yeast Nitrogen Base (YNB ) supplementation. Two yeast strains grew well on pH 7 and 8 with YNB supplementation. All four species achieved the highest biomass level when utilizing sucrose as the main carbon sources, which demonstrates a great pairing with S.elongatus . In addition to investigate the metabolite exchange and growth condition of the co-culture system, we incorporated an Aspergillus strain expressing three genes (PAND, BAPAT and HPDH). Expression of these three genes enabled production of 3-hydroxypropionic acid (3HP) with a titer of 3-6 g/L regardless of pH and other nutrient conditions. Moreover, we developed a computational model with different heterotrophic fungi symbiosing with S.elongatus to evaluate the exchange of hundreds of metabolites. The accuracy and sensitivity had been optimized by high-throughput phenotyping assays. This project demonstrated ways to enhance a synthetic lichen platform’s efficiency and scalability to create a robust system for sustainable bioproduct synthesis. The project has enabled us to explore the technological and commercial potential of sustainable lichen co-culture. The project has also trained the next generation of scientists to address sustainability, carbon fixation, and design solutions to contemporary global challenges. The lichen platform represented a novel approach to sustainably producing bioproducts with a reduced carbon footprint and lower costs.

09 BIOMASS FUELS↗
DOE OSTI2021

Towards a Robust Sampling Approach: A Computational Review and Design

As pointed out in several other works, the estimation of the reliability of the electrical grid can not be conducted without the estimation of the stochastic phenomena of electicity demand and electricity production by variable renewable sources. Therefore, sampling procedures have become integral in the design of engineering structures and analysis. Commonly referred to as Monte Carlo uncertainty analysis or integration, the general objective of these procedures is to establish specifics about the uncertainty of an output characteristic of such an engineering system, given uncertainty about its input characteristics. These sampling procedures are applied in a context in which establishing such specifics cannot be performed through other means. Variance-reduction techniques are designed to lessen the variability among estimators to estimate statistics of those output uncertainties. Importance-sampling techniques, on the other hand, are designed to reduce the number of samples needed to estimate a particular statistic—e.g., a tail probability. The combination of these approaches can reduce the computational burden considerably for a particular estimator and statistic. Importance sampling—geared and designed as it is toward improving a particular estimator—suffers, unfortunately, from the unintended consequence of reducing the performance of other estimators in terms of their variance. The objective of this paper is to offer an alternative sampling procedure where this variance does not grow unacceptably large for a suite of estimators. Moreover, it is anticipated that, with additional knowledge of how an engineering output characteristic responds to its input characteristics, tuning parameters of the input’s sampling procedure can be set to improve the output characteristic’s estimation. This work proves the effectiveness of the suggested alternative approach. Such positive outcome will lead to a decrease of the computational burden of performing stochastic optimization of integrated energy systems (e.g., components dispatch, and portfolio composition). In particular, capturing the contribution to the overall system cost of rare and unlikely events and patterns of the electricity demand and production will become less computationally expensive. This is due to the fact that the approach demonstrated here will allow the sampling of those rare occurrences more frequently without misrepresenting their probabilistic impacts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

A model to assess Zircaloy’s mechanical property changes following a transient beyond critical heat flux

Maintaining the integrity of nuclear fuel rods is essential for ensuring public health and safety in nuclear power generation. During reactor operation, this integrity is confirmed by demonstrating compliance with established regulatory acceptance criteria. For moderate-frequency events, such as limiting transients and anticipated operational occurrences (AOOs), the current fuel integrity criterion is based on preventing boiling transition. This criterion assumes that prevention of boiling transition will prevent excessive cladding heating and, thus, fuel failure during normal operations. While conservative, this approach places significant constraints on core design, fuel cycle economics, and a plant’s ability to perform major power uprates, leading to suboptimal fuel utilization and inefficient carbon-free energy production. A more efficient approach could be achieved by revising the failure criterion to a material-specific limit rather than strictly preventing the boiling transition, since boiling transition per se is not a cause of fuel cladding failure. Here, as a result, a new licensing framework based on material properties, termed time-at-temperature (t@T), is needed. This approach would allow for brief periods of post–critical heat flux operation during an AOO without compromising safety. Implementing the t@T licensing strategy requires a robust technical foundation in material properties, which must be established through comprehensive data collection on both unirradiated and irradiated fuel and cladding materials. This foundation would enable the development of a safety basis that ensures safe operation while providing greater flexibility and efficiency for reactor operation. This paper documents a thorough review of the available data to establish a baseline knowledge that can inform the development of cladding mechanical models, as well as identify experimental data gaps that need to be addressed in future research. Machine learning and data informatics were utilized to extract the importance of parameters on the t@T parameter. Industry tools were used to perform baseline analyses to define the relevant transient conditions for data analysis. The subsequent review successfully identified applicable experimental data, as well as sufficient data to evaluate changes in cladding mechanical properties following an AOO transient. Rather than developing new models, this work coupled existing irradiation annealing and recrystallization models to calculate changes in hardness, yield stress, and ultimate tensile stress following an AOO event. The findings from this review were summarized to highlight the experimental data needs required to fill remaining gaps and support the development of future t@T licensing methodologies.

Cladding performance↗
DOE OSTI2024

${\tt Warm}$SPy: a numerical study of cosmological perturbations in warm inflation

We present ${\tt Warm}$SPy, a numerical code in Python designed to solve for the perturbations' equations in warm inflation models and compute the corresponding scalar power spectrum at CMB horizon crossing. In models of warm inflation, a radiation bath of temperature T during inflation induces a dissipation (friction) rate of strength Q ∝ T c /Φ m in the equation of motion for the inflaton field Φ. While for a temperature-independent dissipation rate (c = 0) an analytic expression for the scalar power spectrum exists, in the case of a non-zero value for c the set of equations can only be solved numerically. For c > 0 (c < 0), the coupling between the perturbations in the inflaton field and radiation induces a growing (decaying) mode in the scalar perturbations, generally parameterized by a multiplicative function G(Q) which we refer to as the scalar dissipation function. Using ${\tt Warm}$SPy, we provide an analytic fit for G(Q) for the cases of c = {3,1,-1}, corresponding to three cases that have been realized in physical models. Compared to previous literature results, our fits are more robust and valid over a broader range of dissipation strengths Q ϵ [10 -7 ,10 4 ]. Additionally, for the first time, we numerically assess the stability of the scalar dissipation function against various model parameters, inflationary histories as well as the effects of metric perturbations. As a whole, the results do not depend appreciably on most of the parameters in the analysis, except for the dissipation index c, providing evidence for the universal behaviour of the scalar dissipation function G(Q).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Using dorsal surface for individual identification of dairy calves through 3D deep learning algorithms

Advances in machine learning techniques have allowed the development of computer vision systems (CVS) that can accurately predict several phenotypes of interest for livestock operations. In this context, 3D images taken from a top-down view are particularly useful for estimating body condition score, growth development, and body biometrics in cattle. Frequently, such CVS rely on identification (ID) systems, such as electronic tags, as a way to match animal ID and the predicted phenotype. However, the same 3D images used to predict body weight and other animal biometrics could be adopted for animal recognition as well. Such alternative would optimize CVS to recognize animal ID and monitor growth development simultaneously while leveraging the same hardware infrastructure. Furthermore, this strategy could be used to recognize animals with similar color patterns. Nonetheless, growing animals are continuously changing body shape, which could limit its use as an invariant feature for pattern recognition. Thus, the objectives of this study were: (1) to compare algorithms for different 3D object representations to identify individual animals; and (2) to evaluate how short-term changes in body shape due to animal growth affect the predictive performance of these algorithms. For objective 1, the algorithms were trained (n = 4,558) and tested (n = 1,139) using images from 38 Holstein calves. For objective 2, we designed three different experiments using images (n = 2,347) from five Holstein calves taken over six weeks during their growing period, always training and testing on different weeks. Each experiment evaluated how changing a different parameter of the image capturing procedure affected the predictive ability of the trained algorithms. In the first experiment, we varied the total number of images per animal in the training set; in the second experiment, we varied the number of weeks while keeping a fixed number of images in the training set; and in the third experiment, we skipped weeks between images in the training and test sets. The F 1 score for objective (1) was up to 0.804 when testing with the last frames of each video, and up to 0.959 when using random frames for testing. For objective (2), the F 1 score was up to 0.947 for the first experiment when using 130 images per animal; up to 0.979 for the second experiment when using all five weeks; and up to 0.917 when not skipping weeks between training and testing. In conclusion, these results show that deep learning algorithms can be used to identify individual animals through their dorsal area 3D surfaces, and, from our experiments using calves in their growing period, that they are robust enough to account for changes in body shape and size, making them a promising tool for animal recognition during growth.

3D neural networks↗
DOE OSTI2021

The AGORA High-resolution Galaxy Simulations Comparison Project. III. Cosmological Zoom-in Simulation of a Milky Way–mass Halo

In this work, we present a suite of high-resolution cosmological zoom-in simulations to z = 4 of a 10 12 M⊙ halo at z = 0, obtained using seven contemporary astrophysical simulation codes (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, and Gizmo) widely used in the numerical galaxy formation community. The physics prescriptions for gas cooling and heating and star formation are the same as the ones used in our previous Assembling Galaxies of Resolved Anatomy (AGORA) disk comparison but now account for the effects of cosmological processes such as the expansion of the universe, intergalactic gas inflow, and the cosmic ultraviolet background radiation emitted by massive stars and quasars. In this work, we introduce the most careful comparison yet of galaxy formation simulations run by different code groups, together with a series of four calibration steps each of which is designed to reduce the number of tunable simulation parameters adopted in the final run. In the first two steps, we methodically calibrate the gas physics, such as cooling and heating, in simulations without star formation. In the third step, we seek agreement on the total stellar mass produced with the common star formation prescription used in the AGORA disk comparison, in stellar-feedback-free simulations. In the last calibration step, we activate stellar feedback, where each code group is asked to set the feedback prescription to as close to the most widely used one in its code community as possible, while aiming for convergence in the stellar mass at z = 4 to the values predicted by semiempirical models. After all the participating code groups successfully complete the calibration steps, we achieve a suite of cosmological simulations with similar mass assembly histories down to z = 4. With numerical accuracy that resolves the internal structure of a target halo (≲100 physical pc at z = 4), we find that the codes overall agree well with one another, e.g., in gas and stellar properties, but also show differences, e.g., in circumgalactic medium (CGM) properties. We argue that, if adequately tested in accordance with our proposed calibration steps and common parameters, high-resolution cosmological zoom-in simulations can have robust and reproducible results. New code groups are invited to join and enrich this comparison by generating equivalent models or to test the code's compatibility on their own, by adopting the common initial conditions, the common easy-to-implement physics package, and the proposed calibration steps. Further analyses of the zoom-in simulations presented here will be presented in forthcoming reports from the AGORA Collaboration, including studies of the CGM, simulations by additional codes, and results at lower redshift.

79 ASTRONOMY AND ASTROPHYSICS↗