Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “function unit testing”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Dark Energy Survey Year 6 Results: Point-Spread Function Modeling

We present the point-spread function (PSF) modeling for weak lensing shear measurement using the full six years of the Dark Energy Survey (DES Y6) data. We review the PSF estimation procedure using the PIFF (PSFs In the Full FOV) software package and describe the key improvements made to PIFF and modeling diagnostics since the DES year three (Y3) analysis: (i) use of external Gaia and infrared photometry catalogs to ensure higher purity of the stellar sample used for model fitting, (ii) addition of color-dependent PSF modeling, the first for any weak lensing analysis, and (iii) inclusion of model diagnostics inspecting fourth-order moments, which can bias weak lensing measurements to a similar degree as second-order modeling errors. Through a comprehensive set of diagnostic tests, we demonstrate the improved accuracy of the Y6 models evident in significantly smaller systematic errors than those of the Y3 analysis, in which all g band data were excluded due to insufficiently accurate PSF models. For the Y6 weak lensing analysis, we include g band photometry data in addition to the riz bands, providing a fourth band for photometric redshift estimation. Looking forward to the next generation of wide-field surveys, we describe several ongoing improvements to PIFF, which will be the default PSF modeling software for weak lensing analyses for the Vera C. Rubin Observatory’s Legacy Survey of Space and Time.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for inverse problems: Determining material properties for heterogeneous systems

Here, we determine spatially varying discontinuous material properties using a domain-decomposition based physics-informed neural networks (PINNs) framework named the Adaptive Interface-PINNs or AdaI-PINNs (Roy et al., 2024). We propose the use of distinct neural networks for the field variables and material properties within each material, utilizing adaptive activation functions. While the neural networks across different materials share the same weights and biases, their activation functions are uniquely tailored using a hyperparameter that influences the slope of the activation function. The proposed framework is tested on several one-dimensional and two-dimensional benchmark examples, and its performance is compared with conventional PINNs and existing domain-decomposition PINNs frameworks, namely, the Multi-domain physics-informed neural network (M-PINN), and the eXtended physics-informed neural networks (XPINNs). The results demonstrate that the proposed approach can determine randomly distributed discontinuous material properties with an L 2 error of $\mathscr{O}$ (10 -3 ) for the material property and the root-mean-square error of $\mathscr{O}$ (10 -3 ) for the primary variable while the other approaches yield errors that are approximately two orders of magnitude larger (that is, $\mathscr{O}$ (10 -1 )). Moreover, the spatial distribution of material properties obtained using the proposed framework is in close agreement with the true distribution, whereas the other approaches fare much worse. Additionally, the proposed approach is approximately 40% faster than its competitors, indicating its potential as a robust alternative for solving inverse problems in heterogeneous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Fermi function and the neutron's lifetime

The traditional Fermi function ansatz for nuclear beta decay describes enhanced perturbative effects in the limit of large nuclear charge Z and/or small electron velocity β. We define and compute the quantum field theory object that replaces this ansatz for neutron beta decay, where neither of these limits hold. We present a new factorization formula that applies in the limit of small electron mass, analyze the components of this formula through two loop order, and resum perturbative corrections that are enhanced by large logarithms. We apply our results to the neutron lifetime, supplying the first two-loop input to the long-distance corrections. Our result can be summarized as τ n x |V ud | 2 [1 + 3λ 2 ] [1 + Δ R ] = $\frac{5263.284(17) s}{1 + 27.04(7) x 10^{-3}}$ with |V ud | the up-down quark mixing parameter, τ n the neutron's lifetime, λ the ratio of axial to vector charge, and Δ R the short-distance matching correction. We find a shift in the long-distance radiative corrections compared to previous work, and discuss implications for extractions of |V ud | and tests of the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Vertical canopy gradients of respiration drive plant carbon budgets and leaf area index

Despite its importance for determining global carbon fluxes, leaf respiration remains poorly constrained in land surface models (LSMs). We tested the sensitivity of the Energy Exascale Earth System Model Land Model – Functionally Assembled Terrestrial Ecosystem Simulator (ELM-FATES) to variation in the canopy gradients of leaf maintenance respiration (R dark ). We ran global and point simulations varying the canopy gradient of R dark to explore the impacts on forest structure, composition, and carbon cycling. In global simulations, steeper canopy gradients of R dark lead to increased understory survival and leaf biomass. Leaf area index (LAI) increased up to 77% in tropical regions compared with the default parameterization, improving alignment with remotely sensed benchmarks. Global vegetation carbon varied from 308 Pg C to 449 Pg C across the ensemble. In tropical forest simulations, steeper gradients of R dark had a large impact on successional dynamics. Results show the importance of canopy gradients in leaf traits and fluxes for determining plant carbon budgets and emergent ecosystem properties such as competitive dynamics, LAI, and vegetation carbon. The high-model sensitivity to canopy gradients in R dark highlights the need for more observations of how leaf traits and fluxes vary along light micro-environments to inform critical dynamics in LSMs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Improving Trustworthiness of Data-Driven Power Grid Contingency Analysis With Bayesian Residual Graph Neural Networks

The evolving energy landscape requires novel tools to efficiently perform contingency analysis and reliability assessment of power grids, potentially in real-time. The high computational cost of traditional power flow solvers limits their applicability in practice. Machine learning (ML) surrogates such as deep neural networks (NNs) accelerate power flow solvers computations, enabling high-order contingency analysis and real-time decision-making by learning highly nonlinear functions and integrating grid topology via graph architectures. However, (graph) NNs lack predictive power away from training data and do not provide predictive confidence estimates. Here, we present a Bayesian residual graph NN that integrates knowledge from low-fidelity data via residual training and embeds granular quantification of uncertainties, improving trustworthiness critical for high-consequence decision-making. Applying Bayesian concepts to NNs is challenging due to the high-dimensionality of both the parameter space, complicating derivation of a meaningful prior, and the output space in large grid systems, requiring enhanced techniques to assess the predicted high-dimensional uncertainties. Our contributions include: (1) Deriving a prior for fully connected and graph NNs that leverages low-fidelity data to guide mean predictions and appropriately control prior predictive uncertainty. (2) Integrating this prior within an ensembling with anchoring scheme for efficient approximate posterior inference. (3) Deriving enhanced metrics to assess accuracy of both the mean and uncertainty predictions in high dimensions, appropriately accounting for correlations propagated through graph layers. The resulting Bayesian residual graph NN is tested on a contingency analysis task for 14-bus and 118-bus grids.

24 - POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Urban weather modeling using WRF: linking physical assumptions, code implementation, and observational needs

The Weather Research and Forecasting (WRF) model includes urban schemes that simulate the influence of urban surfaces on the atmosphere using parameterizations for flux, and radiative exchanges. Three core schemes – the Bulk urban parameterization, Single-Layer Urban Canopy Model (SLUCM), and Multi-Layer Urban Canopy Model (MLUCM) – represent increasing levels of complexity. Although the parameterizations within these urban schemes are described in the literature, their specific implementation remains poorly documented, thus slowing down model development efforts. This manuscript provides a roadmap to the three urban schemes in WRF version 4.5.2, presenting equations using the same symbols as in the model code, along with references to code lines, and including graphics and explanations that connect the code to its physical foundations. Our thorough review of the urban parameterizations implemented in WRF version 4.5.2 highlighted a handful of parameters that may introduce discontinuities in simulations: (i) in the SLUCM, a 1 mm h −1 rain rate threshold is employed to switch between two minimum moisture availability parameterizations, thus impacting latent heat flux calculations; (ii) in the SLUCM a threshold is used to partition shortwave radiation into direct and diffuse components; (iii) in all three urban schemes, the bulk Richardson number is employed to select the similarity function, which influences the vertical distribution of heat and momentum. We also identified a highly simplified treatment of the radiative balance on roof surfaces. The implications of these simplifications can be assessed through targeted observations across relevant conditions, including varying precipitation rates, cloud cover, and transitions between stability regimes. Furthermore, the widespread application of the Monin-Obukhov similarity theory in these urban schemes warrants model evaluation under highly stable and unstable conditions and in heterogeneous urban settings with variable land cover and building heights on scales finer than model resolution. To address these challenges, we offer guidance on observational strategies, emphasizing the need for multi-parametric measurements to capture potential compensating biases and multi-height measurements that align with the levels where quantities are diagnostic and prognosed in the model (i.e., the lowest atmospheric level of the WRF model). Finally, our inspection of the code revealed implementation bugs that have now been corrected in WRF versions 4.6.0 and 4.6.1. Sensitivity tests over the Atlanta urban area show that these corrections affect surface temperatures, underscoring the importance of performing rigorous documentation and verification of the implementation of parameterizations in model code.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Neural Network‐Based Methods for Ocean Surface Wave Measurement Using Submarine Distributed Acoustic Sensing (DAS)

Two new data-driven models for estimating ocean surface waves from distributed acoustic sensing (DAS) submarine cable strain rate are developed using supervised machine learning on a 10-day data set collected offshore of Oliktok Point, Alaska. The new models were trained on target data from seafloor pressure moorings at three sites spaced evenly along 27.1 km of cable and were benchmarked against an empirical transfer function method previously used to estimate waves from DAS. A model which uses convolutional neural networks to transform 2-km frequency-wavenumber strain spectra to seafloor pressure spectra outperforms the benchmark in wave height prediction (RMSE of 0.15 vs. 0.41 m) and period prediction (0.29 vs. 0.37 s) when evaluated on a held-out test data set. When applied to a DAS data set collected on the same cable 2 years prior, the CNN-based model maintained similar significant wave height performance (RMSE = 0.23 m) relative to available satellite altimetry data. A two-hidden-layer, fully connected neural network which transforms 1-D strain spectra to seafloor pressure spectra also outperforms the benchmark in wave height prediction (RMSE of 0.19 vs. 0.41 m), but does not generalize as well to the prior data. Regression-based machine learning is useful for estimating waves from DAS data when the pressure-strain relationship varies temporally and spatially across different wave conditions. Models can be applied to DAS data to measure waves with higher spatial resolution and longer temporal coverage than traditional methods, which often measure waves only at a single point.

Davis, Jacob R. [Univ. of Washington, Seattle, WA ↗
DOE OSTI2025

Extended molecular eigenmodes treatment of dipole–dipole NMR relaxation in real fluids

Traditional models of NMR relaxation fail to account for the complex, multi-exponential behavior of the autocorrelation function in realistic systems characterized by soft-interactions and molecules that are chemically and physically complex. Here, in this study, we describe the relative diffusion of the spin dipoles by means of a Fokker–Planck equation that includes an interaction potential of mean force to account for the response of the physical/chemical environment around the dipoles. By numerically solving the Fokker–Planck equation for the diffusion propagator, we estimate dipole–dipole NMR relaxation for like- and unlike-spin systems via its eigenmode solution. We test the model against molecular simulations of diffusing dipoles with harmonic potentials and also validate using experimental longitudinal relaxation data from real systems, including Gd(III)–aqua and Gd(III)–DO3A–butrol complexes, the latter being an important MRI contrast agent. Using this novel approach, we predict both the inner- and outer-shell contributions to the relaxivity rates with excellent accuracy at frequencies relevant to MRI. We also show that, under the appropriate assumptions, our framework naturally recovers the Bloembergen–Purcell–Pound, the Solomon–Bloembergen–Morgan, and the Hwang–Freed models. Our implementation is general and publicly available for application to a broad range of systems.

Pinheiro dos Santos, Thiago J. [Rice Univ., Housto↗
DOE OSTI2025

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗
DOE OSTI2026

Uniaxial tensile mechanical properties from a scratch on ultra-fine grained tungsten

We demonstrate that surface-based mechanical testing can be used to quantify a wide range of mechanical properties determined using traditional uniaxial tension experiments. These include mapping of dynamic stress–strain behaviour, which can be used as a screening tool for uniaxial tensile ductility (uniform elongation), fracture toughness (K C ), which is shown to be in good agreement with reported plane-strain K IC , as well work-hardening exponents (n), activation volumes (V) energies (Q) and strain-rate sensitivity (m) exponents as a function of temperature. Advantages of the surface-based approach are the ability to work with significantly smaller volumes and accommodate processing defects like pores, voids and cracks. These methods are applied to ultra-fine grained tungsten (average grain size of 650 nm), as it is an ideal candidate for alloy development efforts in advanced energy systems. Tungsten is a useful baseline as it exhibits a near-ambient temperature (approx. 200°C) ductile-to-brittle transition temperature (DBTT), enabling validation efforts for the proposed methods, and extraordinary material properties, including the highest melting temperature of any metal and complex temperature and strain-rate dependencies linked to transitions in dislocation mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The NREL Sensor Laboratory Detection of Hydrogen Emissions

The development of a functional hydrogen detection system is a multifaceted process that integrates hardware, deployments strategies, and analytics which can be supported by the NREL Sensor Laboratory: 1. Support of the design, validation and optimization of sensing prototypes; 2. Guide optimized sensing element development, including control electronics; 3. Laboratory testing to validate/optimize metrological performance (measurement range, detection limit, etc.); 4. Provide test sites for field deployments representative of real-world scenarios with controlled hydrogen releases; 5. Develop sensor placement and operation guidance; 6. Provide guidance on electronics to accommodate facility integration; 7. Electrical safety designs to allow for operation within restricted zones; 8. Integration into facility monitoring and control systems; 9. Guide incorporation of cyber security elements to protect facilities from malicious attacks; 10. Modeling and application of advanced analytics to detect and quantify emissions; 11. Higher Order dispersion models to guide sensor placement for reliable detection; 12. Advanced analytics for improved metrological performances, and to inform inverse modeling; 13. Market support and commercialization (national and international markets); 14. Commercial deployments in H2@SCALE markets (e.g., HUBs and other large-scale hydrogen markets); and 15. Leverage off international collaborations/partnerships (e.g., NREL is on the advisory board for the European initiative "pre-Normative Research on Hydrogen Releases Assessment"-NHyRA).

08 HYDROGEN↗
DOE OSTI2019

SECARB Project Overview Book

In 2003, the United States Department of Energy (DOE) issued awards initiating seven Regional Carbon Sequestration Partnerships (RCSP) spanning the United States and portions of Canada. Managed by DOE’s National Energy Technology Laboratory (NETL), the RCSPs currently represent more than 400 state agencies, universities, and private companies, spanning 43 states, three Native American Organizations, and four Canadian provinces. The Southern States Energy Board (SSEB) received an award on October 1, 2003, which established the Board’s overall management and administration of the Southeast Regional Carbon Sequestration Partnership, or SECARB. The geographical region currently includes 13 states and a network of more than 100 stakeholders. The Nation’s leading scientists, university researchers, national laboratories, industrial representatives, environmental organizations, and many others have taken a proactive approach at identifying and characterizing the most promising options for technology deployment and carbon dioxide (CO2) utilization and storage (CCUS) for the Southeast. The results obtained during these injections will be important to the future commercialization of CCUS technologies and the continued use of fossil fuels as a significant energy source in a manner that is environmentally responsible. Each of the six SECARB project locations has been locally coordinated by its own field team. The field teams assume responsibility for the technical scope of work, local education and outreach, permitting, monitoring, verification, and accounting (MVA), and maintaining the validation test’s schedule and budget. Each team contributes new information to the continued characterization of the region. In addition, a task is dedicated to integrating field data and filling gaps in regional characterization data sets. Data and tools developed in the Continued Characterization task are incorporated into a relational database and geographic information system (GIS). All field tests, the continued characterization project, and the cross-cutting functions are designed to support the DOE roadmap by validating technologies and identifying locations throughout the region that could support full-scale CCUS deployment opportunities. Detailed fact sheets for the SECARB projects are available on the SECARB website at www.secarbon.org (see “Projects” tab).

RCSP↗
DOE OSTI2025

Dual-season thermal energy storage-integrated heat pump system for load shifting: experiments, modeling and nationwide analysis

Here, this paper proposes and validates a dual-season thermal energy storage-integrated heat pump (TES-HP) system that shifts building electric load in both cooling and heating modes using a single phase-change TES. The design uses the TES as a heat sink for the condenser in summer and a heat source for the evaporator in winter, thereby reducing on-peak electric demand while ensuring compatible with existing air-distribution systems. A 14-kW TES-HP prototype with six operating modes and novel refrigerant charge management strategies was experimentally tested under varied conditions. Subsequently, data-driven polynomial performance curves were developed and validated against experiments, then coupled with a phase change material (PCM) model and a calibrated building model within a rule-based controller. Single-day case studies demonstrated effective on-peak demand reduction while maintaining thermal comfort. During cooling, hourly power was lowered by up to 1.5 kW; in heating, average hourly power decreased from over 3.5 kW to 1.8 kW, compared to an HP-only system, achieving 64.5 % load shifting, and the need for resistance heating was eliminated. Seasonal simulations showed typical on-peak electricity savings of 25–35 % in cooling and 40–65 % in heating, with the largest benefits on the hottest and coldest days. Extended response-surface analysis and nationwide mapping quantified load shifting as functions of ambient and TES temperatures, revealing a seasonal trade-off in TES phase-change temperature. These results demonstrate the TES-HP system as a practical and scalable solution for grid-interactive buildings that reduces on-peak demand and electricity use while maintaining thermal comfort and enhancing grid flexibility and reliability.

Grid-interactive buildings↗
DOE OSTI2025

Neutron irradiation & thermomechanical experiment (NITE) - design

For the reliable long-term operation of fusion power plants, it is crucial to understand and predict the lifetime of materials in use. These materials include all structural and functional materials utilized at the first wall, blanket, magnets, and shielding. The key challenge is, that the harsh environment including high heat fluxes, high thermal stress and stress cycling, neutron irradiation, and sputtering on such materials should not be viewed separately. Currently, the synergistic loads cannot be evaluated experimentally because of the lack of adequate facilities. The purpose of that work is to design a synergetic Neutron Irradiation and Thermomechanical Experiment (NITE) for fusion materials. This design will leverage the existing Advanced-Test-Reactor (ATR), a fission reactor at the Idaho National Laboratory. We also acknowledge that with existing fission reactors the exact fusion condition cannot be created, and the limitations are critically discussed. The combination of neutron irradiation with a high heat flux is the focus. This is realized with an irradiation capsule design that includes a TRISO fueled region inside the capsule to enable a steady-state heat flux on one side of the specimen. In conclusion, the experimental design modeling showed that steady-state heat fluxes of 2.4 MW/m 2 with a thermal gradient of above 250°C can be achieved in a 5 mm thick specimen.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

StreamGen: Connecting Populations of Streams and Shells to Their Host Galaxies

In this work, we study how the abundance and dynamics of populations of disrupting satellite galaxies change systematically as a function of host galaxy properties. We apply a theoretical model of the phase-mixing process to classify intact satellite galaxies and stellar streamlike and shell-like debris in ∼1500 Milky Way–mass systems generated by a semi-analytic galaxy formation code, SatGen. In particular, we test the effect of host galaxy halo mass, disk mass, ratio of disk scale height to length, and stellar feedback model on disrupting satellite populations. We find that the counts of tidal debris are consistent across all host galaxy models, within a given host mass range, and that all models can have streamlike debris on low-energy orbits, consistent with that observed around the Milky Way. However, we find a preference for streamlike debris on lower-energy orbits in models with a thicker (lower-density) host disk or on higher-energy orbits in models with a more massive host disk. Importantly, we observe significant halo-to-halo variance across all models. These results highlight the importance of simulating and observing large samples of Milky Way–mass galaxies and accounting for variations in host properties when using disrupting satellites in studies of near-field cosmology.

dark matter↗
DOE OSTI2025

Population Balance Models for Catalytic Depolymerization: From Elementary Steps to Multiphase Reactors

Here, the ongoing accumulation of plastic waste in landfills and in the environment is driving research on chemical processes and catalysts to recycle polymers. Traditional modeling strategies are not applicable to these processes because they involve too many reactants and intermediates, one for each molecular weight and each functionalization. To model the kinetics, we have developed population balance models (PBMs) that account for macromolecular reactants in the bulk and macromolecular catalytic intermediates. These PBMs couple to each other through polymer adsorption and desorption models and to traditional rate equations for small molecule products and co-reactants (like hydrogen or ethylene). The models, in combination with experimental data, are being used in many ways: (i) to test mechanistic hypotheses, (ii) to extract rate parameters, (iii) to quantitatively compare catalyst activities, (iv) to account for mass transfer and vapor–liquid partitioning in two-phase reactors, and (v) to design novel support architectures and catalysts that mimic the processive action of natural depolymerization enzymes. Some key theoretical advances allow PBMs to be constructed from elementary rates and mechanisms, as opposed to traditional formulations with pseudoelementary rate parameters invoked as fitting parameters. We discuss ways to build these models “bottom-up” from first-principles calculations and ways to extract model parameters from “top down” analyses of rate data. The combination provides a quantitative bridge between first-principles calculations and the kinetics of complex macromolecular transformations for polymer upcycling and beyond.

Manis, Lela K. [University of Illinois at Urbana-C↗
DOE OSTI2022

The ICON-A model for direct QBO simulations on GPUs (version icon-cscs:baf28a514)

Abstract. Classical numerical models for the global atmosphere, as used for numerical weather forecasting or climate research, have been developed for conventional central processing unit (CPU) architectures. This hinders the employment of such models on current top-performing supercomputers, which achieve their computing power with hybrid architectures, mostly using graphics processing units (GPUs). Thus also scientific applications of such models are restricted to the lesser computer power of CPUs. Here we present the development of a GPU-enabled version of the ICON atmosphere model (ICON-A), motivated by a research project on the quasi-biennial oscillation (QBO), a global-scale wind oscillation in the equatorial stratosphere that depends on a broad spectrum of atmospheric waves, which originates from tropical deep convection. Resolving the relevant scales, from a few kilometers to the size of the globe, is a formidable computational problem, which can only be realized now on top-performing supercomputers. This motivated porting ICON-A, in the specific configuration needed for the research project, in a first step to the GPU architecture of the Piz Daint computer at the Swiss National Supercomputing Centre and in a second step to the JUWELS Booster computer at the Forschungszentrum Jülich. On Piz Daint, the ported code achieves a single-node GPU vs. CPU speedup factor of 6.4 and allows for global experiments at a horizontal resolution of 5 km on 1024 computing nodes with 1 GPU per node with a turnover of 48 simulated days per day. On JUWELS Booster, the more modern hardware in combination with an upgraded code base allows for simulations at the same resolution on 128 computing nodes with 4 GPUs per node and a turnover of 133 simulated days per day. Additionally, the code still remains functional on CPUs, as is demonstrated by additional experiments on the Levante compute system at the German Climate Computing Center. While the application shows good weak scaling over the tested 16-fold increase in grid size and node count, making also higher resolved global simulations possible, the strong scaling on GPUs is relatively poor, which limits the options to increase turnover with more nodes. Initial experiments demonstrate that the ICON-A model can simulate downward-propagating QBO jets, which are driven by wave–mean flow interaction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Statistical Uncertainty of Inhalation Dose Coefficients in Consequence Management: Propagated Dose Uncertainty in ICRP 66 Human Respiratory Tract Model

Reference inhalation dose models rely on deterministic biokinetics and reference computational phantoms, limiting their applicability to the variability present in population-specific exposures encountered in emergency response scenarios. Here, this study introduces REDCAL, a Python-based computational framework developed to propagate uncertainty in inhalation dose coefficients using the International Commission on Radiological Protection (ICRP) Publication 66 Human Respiratory Tract Model. REDCAL integrates ICRP deposition and clearance models, systemic biokinetics, and governing physics principles, and leverages Sandia National Laboratories’ Dakota toolkit for uncertainty quantification via Latin Hypercube Sampling. REDCAL was validated against DCAL, with biokinetic retention results differing by less than 1% and effective dose coefficients by less than 2% across all tested radionuclides. Stochastic sampling introduced variability in dose coefficients, with geometric standard deviations (GSD) in committed effective dose coefficients (CEDC) ranging from 1.0 to 1.5, based on lognormal distribution fits. Analysis demonstrated that variations in the activity median aerodynamic diameter (AMAD) notably influenced the computed CEDC values. Smaller particles (<1 µm) increased doses by 20–30% due to deeper lung deposition and prolonged retention for alpha emitting radionuclides, such as 241 Am and 239 Pu. Radionuclides with fast clearance, such as 133 I, demonstrated a dose reduction exceeding 50%, as AMAD increased beyond 5 µm due to upper airway deposition and rapid mucociliary clearance. The greatest GSD among the radionuclides reported in this study was for 241 Am. In most cases, the largest GSDs in the CEDC were associated with larger particle sizes, an expected outcome, as ICRP Publication 66 defines GSD in particle size as a function of AMAD, resulting in an extended tail of the lognormal distribution. The findings support improved inhalation dose assessments and enhance consequence management strategies for the U.S. Federal Radiological Monitoring and Assessment Center by quantifying uncertainty in dose coefficients and strengthening decision-making for emergency response scenarios.

Biokinetic Modeling↗