Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reproducibility”

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

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

At least 1,405 records · Page 78

"Reducing Detailed Vehicle Energy Dynamics to Physics-Like Models"

The energy demand of vehicles, particularly in unsteady drive cycles, is affected by complex dynamics internal to the engine and other powertrain components. Yet, in many applications, particularly macroscopic traffic flow modeling and optimization, structurally simple approximations to the complex vehicle dynamics are needed that nevertheless reproduce the correct effective energy behavior. This work presents a systematic model reduction pipeline that starts from complex vehicle models based on the Autonomie software and derives a hierarchy of simplified models that are fast to evaluate, easy to disseminate in open-source frameworks, and compatible with optimization frameworks. The pipeline, based on a virtual chassis dynamometer and subsequent approximation strategies, is reproducible and is applied to six different vehicle classes to produce concrete explicit energy models that represent an average vehicle in each class and leverage the accuracy and validation work of the Autonomie software.

Khoudari, Nour↗

GeoThermalCloud for EGS – An Open-source, User-friendly, Scalable AI Workflow for Modeling Enhanced Geothermal Systems

Enhanced Geothermal Systems (EGS) offer a vast potential to expand the use of geothermal energy. Heat is extracted from this engineered system by injecting relatively cold water into subsurface fractures, which are in contact with hot dry rock, and brought back to surface through production wells. Creating EGS requires improving the natural permeability of hot crystalline rocks. In this short conference paper, we present a reproducible workflow for modeling EGS. Our workflow called the GeoThermalCloud (GTC) for EGS, leverages recent advances in machine learning, deep learning, and high-performance computing. This GTC framework is currently being made open-source, user-friendly, and reproducible through python scripts as well as Google Colab/Jupyter Notebooks. This GTC for EGS modeling scripts are made available at https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS and will constantly be updated to cater for geothermal community. Current GTC framework provides scripts to train deep learning (DL) models for techno-economics and data worth analysis. The Geothermal Design Tool (https://github.com/GeoDesignTool/GeoDT.git), a fast and simplified multi-physics solver, is used to develop a database for training DL models. This short paper provides details on the scripts to curate, process, and train DL models. The scripts can easily be modified to train on databases generated by other popular open-source simulators such as PFLOTRAN, STOMP, TOUGH, and GEOSX or commercial software such as ResFrac and COMSOL.

15 GEOTHERMAL ENERGY↗

Measurement of Nuclear Dependence in Inclusive Antineutrino Scattering

One of the fundamental questions in physics is the matter-antimatter asymmetry of the Universe, explored through charge-parity violation searches in neutrino oscillation experiments. Next-generation oscillation experiments like DUNE need to constrain the systematic uncertainties arising from our understanding of (anti)neutrino-nucleus scattering to within a few percent to achieve their physics goals. MINERvA, a dedicated (anti)neutrino-nucleus scattering experiment, plays a crucial role in constraining these uncertainties. This thesis presents MINERvA's first high-statistics direct measurement of nuclear dependence in charged-current inclusive antineutrino scattering on carbon, hydrocarbon, iron, and lead as a function of antimuon transverse momentum, $p_{\text{T}}$, and Bjorken~$x$. The analysis utilises interactions with a mean antineutrino energy of approximately~6~GeV, with reconstructed antimuons having a scattering angle of less than 17$^\circ$ relative to the antin eutrino beam and an antimuon energy of 2--20 GeV. The measured per-nucleon differential cross-sections are reported with a precision of 7--9\%, while the cross-section ratios of carbon, iron, and lead to hydrocarbon have uncertainties of 5\% or less. The cross-sections for iron and lead indicate strong suppression at low $p_{\text{T}}$ and Bjorken~$x$, and an enhancement at high $p_{\text{T}}$. These effects are observed to be more pronounced with the increasing size of the target nucleus and are not reproduced by the underlying simulation prediction. Comparisons to alternative models used in current (anti)neutrino interaction generators show some improvements in modelling over the base prediction model, yet they are still unable to fully reproduce the observed nuclear dependence in this analysis. Importantly, the analysis provides a direct test of nuclear effects in inclusive antineutrino scattering, with major contributions from resonant pion production, deep inelastic scattering, and the transition region between these channels, which will be significant in DUNE. This measurement also represents one of the largest antineutrino datasets in this energy regime analysed to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Initial Development of a new RELAP5-35 Model of HTTF

Previous validation studies of RELAP5-3D against HTTF data demonstrated an ability to reproduce trends in temperature but an inability to reproduce the transient temperature rise. This is hypothesized to be the result of modeling decisions made when the initial RELAP5-3D model was develop prior to the experiments. To test this hypothesis, a new RELAP5-3D model is being developed. This presentation shows results from full power steady-state and in a pressurized conduction cooldown accident using the new model and compares them to results from the legacy model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Comparison of Results between the Legacy and Refined RELAP5-3D Models of the High Temperature Test Facility in Exercises 1 and 2 of the HTTF Benchmark

Work conducted in FY23 identified that RELAP5-3D was capable of reproducing trends in HTTF data during experiment PG-27 but was incapable of reproducing measured values. The primary cause of this discrepancy between RELAP5-3D results and experimental data was hypothesized to be a distortion in power density that was introduced by the radial nodalization of the model. We further hypothesized that a new model would provide better results when compared to the experiments PG-27 and PG-29. Work this FY developed a new model that is better capable of capturing local heat generation rates and contains a representation of each 1/6 azimuthal sector of the core. We used this model to develop a new set of solutions to Exercises 1 and 2 of Problems 2 and 3 in the benchmark. In this report, we present the first comprehensive comparison of the results between the two models. We see that in Exercise 1A, which is common between problems 2 and 3, the results are similar, though the results from the new model show greater detail than those from the legacy model. In Problem 2 Exercise 1B and Problem 3 Exercise 1B, we see that heat removal is slower in the new model than the legacy model. Problem 3 Exercise 1C shows temperatures that are lower in most places in the new model than the legacy model, but the area with active heat generation has higher block temperatures in the new model than the legacy model. Problem 3 Exercise 1D further shows that long-term heat removal is lower in the new model. Problem 2 Exercise 1C demonstrated that the new model observes higher temperatures in the core regions than the legacy model, justifying the need to preserve the power density in HTTF. The validation of PG-27 and PG-29 also demonstrated the improved temperature agreement in the core regions, particularly with a calibrated model that implements an effective thermal conductivity for the core material. Overall, PG-27 models show reasonable to excellent agreement for steady-state temperatures and minimal to reasonable agreement for transients. PG-29 models showed minimal to insufficient agreement with the data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An R shiny graphical user interface for highprecision mass spectrometric data analysis

• There is currently a lack of software that meets the needs for the analysis of raw data produced by modern isotope ratio mass spectrometers for both R&D and routine use at SRNL and other US national labs • Needs to accommodate multiple isotope systems, instruments, and manufacturers • Include modern statistical methods and handling/visualization of uncertainty • Flexible software with transparent (no “black box”) and reproducible methods • This project is inspired by existing discipline-specific data analysis software (e.g., Tripoli1 , ET_Redux2, IsoplotR3) used in the geochemical community • Our goal is to build an open source data analysis software package that focuses on flexibility, transparency, and reproducibility

LABONE, ELIZABETH↗

Enhancing Cluster Identification in Atom Probe Tomography Data Using Transfer Learning

Atom Probe Tomography (APT) is a powerful technique for visualizing the atomic-scale distribution of solutes in materials, but quantitative cluster analysis of APT datasets remains a challenge due to the need for subjective parameter selection in clustering algorithms. While distance-based and density-based methods such as HDBSCAN are widely used, their performance is highly sensitive to user-defined parameters, which undermines reproducibility and accuracy. This study proposes an image-based, deep learning-aided workflow for automating parameter selection and cluster detection in APT data analysis. By projecting 3D APT point clouds onto 2D planes, we leverage pretrained convolutional neural networks (ConvNeXt-Tiny and ResNet-50) through transfer learning to predict the number of clusters present in synthetic datasets. The output is used to guide K-means clustering and estimate HDBSCAN parameters, specifically minimum cluster size and minimum sample points. This approach reduces reliance on manual parameter tuning, improving consistency and scalability. The methodology demonstrates the feasibility of using image-based deep learning for interpreting complex spatial patterns in APT data, enabling faster and more objective analysis. The complete workflow and code are made publicly available to support reproducibility and future research.

Density-based clustering↗

LeWRON: Agentic Analysis of Electroweak Phase Transitions

The electroweak phase transition (EWPT) is a central topic in particle physics and cosmology, connecting collider phenomenology, baryogenesis, and gravitational-wave observatories. Its analysis requires a technically demanding, convention-sensitive, and model-dependent pipeline, from constructing the finite-temperature effective potential to tracking thermal histories, computing bubble nucleation rates, and predicting gravitational-wave spectra. We present LeWRON (Learning ElectroWeak phase tRansitiON), an agentic framework that orchestrates this pipeline starting from an input Lagrangian. LeWRON combines audited toolbox construction with an Explorer module that uses the generated model-specific code for further analysis, including scans and plots. Intermediate analytic outputs are checked by auditor agents and stored as structured artifacts, enabling reproducible human inspection and downstream use through both a command-line interface and a public Python API. The framework supports a reproduction mode, which infers conventions from the literature and reproduces published results, and a discovery mode, which guides users through structured checkpoints for new models. We demonstrate LeWRON across representative beyond-the-Standard-Model scenarios and release the code on GitHub.

Wang, Isaac R. [Fermilab] (ORCID:000000030789218X)↗

DEM Modeling and Validation of Pebble Bed Packing Using Chrono::GPU

Accurate prediction of pebble packing structure is important for pebble bed reactors because the spatial distribution of void fraction directly affects coolant flow, pressure drop, heat transfer, and neutronic behavior. However, experimentally validated DEM studies that directly evaluate local void-fraction structure in reactor-relevant pebble beds remain limited. In this work, the pebble bed experiment conducted at Missouri University of Science and Technology is simulated using the graphics processing unit (GPU)-based discrete element method (DEM) code Chrono::GPU. The study focuses on evaluating the ability of Chrono::GPU to reproduce the packing arrangement and void-fraction distribution of a randomly packed spherical pebble bed. The DEM results are first verified against established radial void-fraction correlations, including the Mueller and Vortmeyer-Schuster models, to assess the predicted bulk porosity, near-wall behavior, and oscillatory packing structure. The simulation is then verified against reference DEM data and validated against gamma-ray computed tomography (CT) experimental data at three axial locations. The Chrono::GPU results reproduce the main features of the experimental packing, including the high void fraction near the wall, the first near-wall trough, and the damped oscillatory radial profile caused by wall-induced ordering. Quantitative comparison with DEM data and the CT-based radial profiles shows good agreement, with mean absolute errors on the order of 0.07 and root-mean-square errors below 0.09 for the averaged profiles. These results demonstrate that Chrono::GPU can accurately capture the void-fraction structure of spherical pebble beds and provides a reliable DEM framework for future pebble bed reactor packing, recycling, and thermal-hydraulic studies.

97 - MATHEMATICS AND COMPUTING↗

Machine learning for slow-extraction uniformity at the Fermilab Delivery Ring

The Muon-to-Electron Conversion Experiment (Mu2e) requires uniform slow extraction from the Fermilab Delivery Ring during each 43 ms spill. This project analyzed 13,737 spills from seven parquet files collected during March and June 2026 commissioning. Each spill contained 430 samples at 0.1 ms resolution and was summarized by the Spill Duty Factor. Null Action and active proportional-integral-derivative conditions established operational baselines. Fixed-ramp sessions were segmented to evaluate reproducibility, drift, and spill-to-spill spread, while a deliberate quadrupole-ramp step identified the beam-response peak and trough. Median profiles were reproducible within fixed conditions but drifted over multi-hour periods. The resulting baselines and response features support future physics-based and machine-learning approaches to real-time spill regulation.

Prescott, Matthew J. [Purdue U., West Lafayette; F↗

Modeling of a Four-Stage Linear Ionization Cooling Channel for a Muon Collider in G4Beamline

A previous study of an eight-stage rectilinear ionization cooling channel in the ICOOL software demonstrated a five-order-of-magnitude reduction in a muon beam’s 6D emittance. In this study, we look to compare the ways ICOOL and Muons, Inc.’s g4Beamline software model ionization cooling by comparing their modeling of the first four stages to this optimized cooling channel constructed in ICOOL. We begin by identifying the parameters used to construct the optimized ionization cooling channel in ICOOL. We then reconstruct this beam in g4Beamline with identical parameter specifications and simulate the cooling of an identical input beam. Finally, we compare the two simulations based on their beam transmission, longitudinal emittance, and transverse emittance along the channel length. Through this process, we demonstrate that G4Beamline accurately reproduces transverse cooling results but predicts systematically different longitudinal emittance evolution while maintaining similar overall cooling performance, reproducing a 97.9% reduction in 6D emittance over four stages.

Keeler, Dominic [Purdue U., West Lafayette] (ORCID↗

Emergent Degradation Phenomena Demonstrated on Resilient, Flexible, and Scalable Integrated Photoelectrochemical Cells

Photoelectrochemical (PEC) water splitting provides a pathway to generate sustainable clean fuels using the two most abundant resources on Earth: sunlight and water. Currently, most of the successful models of PEC cells are still fabricated on small scales near 1 cm 2 , which largely limits the mass deployment of solar-fuel production. Here, the scale-up to 8 cm 2 of an integrated PEC (IPEC) device is demonstrated and its performance compared to a 1 cm 2 IPEC cell, using state-of-the-art iridium and platinum catalysts with III–V photoabsorbers. The initial photocurrents at 1 sun are 8 and 7 mA cm -2 with degradation rates of 0.60 and 0.47 mA cm -2 day -1 , during unbiased operation for the 1 and 8 cm 2 devices, respectively. Additionally, evaluating under outdoor and indoor conditions at two U.S. National Laboratories reveals similar results, evidencing the reproducibility of this design's performance. Furthermore, the emerging degradation mechanisms during scale-up are investigated and the knowledge gained from this work will provide feedback to the broader community, since PEC device durability is a limiting factor in its potential future deployment.

25 ENERGY STORAGE↗

Water Splitting: Emergent Degradation Phenomena Demonstrated on Resilient, Flexible, and Scalable Integrated Photoelectrochemical Cells (Adv. Energy Mater. 48/2020)

We report Photoelectrochemical (PEC) water splitting provides a pathway to generate sustainable clean fuels using the two most abundant resources on Earth: sunlight and water. Currently, most of the successful models of PEC cells are still fabricated on small scales near 1 cm 2 , which largely limits the mass deployment of solar-fuel production. Here, the scale-up to 8 cm 2 of an integrated PEC (IPEC) device is demonstrated and its performance compared to a 1 cm 2 IPEC cell, using state-of-the-art iridium and platinum catalysts with III-V photoabsorbers. The initial photocurrents at 1 sun were 8 and 7 mA cm -2 with degradation rates of 0.60 and 0.47 mA cm -2 day -1 , during unbiased operation for the 1 and 8 cm 2 devices, respectively. Evaluating under outdoor and indoor conditions at two US National Laboratories revealed similar results, evidencing the reproducibility of this design’s performance. Furthermore, the emerging degradation mechanisms during scale-up are investigated and the knowledge gained from this work will provide feedback to the broader community, since PEC device durability is a limiting factor in its potential future deployment.

25 ENERGY STORAGE↗

A variational multiscale immersed meshfree method for heterogeneous materials

Abstract We introduce an immersed meshfree formulation for modeling heterogeneous materials with flexible non-body-fitted discretizations, approximations, and quadrature rules. The interfacial compatibility condition is imposed by a volumetric constraint, which avoids a tedious contour integral for complex material geometry. The proposed immersed approach is formulated under a variational multiscale based formulation, termed the variational multiscale immersed method (VMIM). Under this framework, the solution approximation on either the foreground or the background can be decoupled into coarse-scale and fine-scale in the variational equations, where the fine-scale approximation represents a correction to the residual of the coarse-scale equations. The resulting fine-scale solution leads to a residual-based stabilization in the VMIM discrete equations. The employment of reproducing kernel (RK) approximation for the coarse- and fine-scale variables allows arbitrary order of continuity in the approximation, which is particularly advantageous for modeling heterogeneous materials. The effectiveness of VMIM is demonstrated with several numerical examples, showing accuracy, stability, and discretization efficiency of the proposed method.

36 MATERIALS SCIENCE↗

Validating automated resonance evaluation with synthetic data

The integrity and precision of nuclear data are crucial for a broad spectrum of applications, from national security and nuclear reactor design to medical diagnostics, where the associated uncertainties can significantly impact outcomes. A substantial portion of uncertainty in nuclear data originates from the subjective biases in the evaluation process, a crucial phase in the nuclear data production pipeline. Recent advancements indicate that automation of certain routines can mitigate these biases, thereby standardizing the evaluation process and enhancing reproducibility. This research aims to provide a methodology, framework, and metrics for the validation of automated nuclear data evaluation software leveraging high-quality synthetic data that closely mimic real experimental observables. An introduced error metric provides a scale and intuitive measure of the evaluation quality by quantifying the estimate’s accuracy and performance across the specified energy range. Synthetic data provides access to experimental observables and underlying resonance parameters, enabling comparison of different evaluations. The methodology is demonstrated using Ta-181 isotope data in the resolved resonance region. The Automated Resonance Identification Subroutine (ARIS), which operates without prior resonance information, was used to test and showcase the framework’s capabilities utilizing the proposed error metrics. The results demonstrate the effectiveness of the proposed approach and framework for optimizing software parameters and testing hypotheses through “what-if” controlled experiments, such as modifying assumptions about experimental conditions or average resonance parameters.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Improving the performance of PM6 donor polymer by random ternary copolymerization of BDD and DTBT segments

Developing high-performance wide bandgap polymer donors to match the rapidly growing non-fullerene small molecule acceptors (NF-SMAs) is of great importance to further improve the photovoltaic performances of organic solar cells (OSCs). PM6 and D18 are two leading wide bandgap donor polymers, but their performances are highly sensitive to the molecular weights, which have distinct impact on the morphology and thereby the charge generation and recombination processes. Here, in this work, we reported the random ternary copolymerization of PM6 and D18 segments to develop random terpolymers with finely tuned photoelectronic properties and optimized blend morphology for high-performance OSCs. The resulting random terpolymers (PMD-5, PMD-10, PMD-15 and PMD-20) exhibited finely tuned absorptions, energy levels and blend morphology. The optimal morphology of PMD-15:L8-BO contributed to the fast exciton diffusion and dissociation, efficient charge extraction with suppressed recombination. Thus, the power conversion efficiency was improved from 17.22% of the control PM6 devices and 17.47% of the control D18 devices to 18.12% of PMD-15 based devices, which is the top performance of random terpolymers. Moreover, the PMD-15 exhibited much higher batch-to-batch reproducibility at a broad molecular weight range than that of PM6 and D18, showing promising for practical applications.

36 MATERIALS SCIENCE↗

A close look onto structural models and primary ligands of metallo-β-lactamases

β-Lactamases are hydrolytic enzymes capable of opening the β-lactam ring of antibiotics such as penicillin, thus endowing the bacteria that produce them with antibiotic resistance. Of particular medical concern are metallo-β-lactamases (MBLs), with an active site built around coordinated Zn cations. MBLs are pan-reactive enzymes that can break down almost all classes of β-lactams, including such last-resort antibiotics as carbapenems. They are not only broad-spectrum-reactive but are often plasmid-borne (e.g., the New Delhi enzyme, NDM), and can spread horizontally even among unrelated bacteria. Acquired MBLs are encoded by mobile genetic elements, which often include other resistance genes, making the microbiological situation particularly alarming. There is an urgent need to develop MBL inhibitors in order to rescue our antibiotic armory. A number of such efforts have been undertaken, most notably using the 3D structures of various MBLs as drug-design targets. Structure-guided drug discovery depends on the quality of the structures that are collected in the Protein Data Bank (PDB) and on the consistency of the information in dedicated β-lactamase databases. Here, we conducted a careful review of the crystal structures of class B β-lactamases, concluding that the quality of these structures varies widely, especially in the regions where small molecules interact with the macromolecules. In a number of examples the interpretation of the bound ligands (e.g., inhibitors, substrate/product analogs) is doubtful or even incorrect, and it appears that in some cases the modeling of ligands was not supported by electron density. For ten MBL structures, alternative interpretations of the original diffraction data could be proposed and the new models have been deposited in the PDB. In four cases, these models, prepared jointly with the authors of the original depositions, superseded the previous deposits. This review emphasizes the importance of critical assessment of structural models describing key drug design targets at the level of the raw experimental data. Since the structures reviewed here are the basis for ongoing design of new MBL inhibitors, it is important to identify and correct the problems with ambiguous crystallographic interpretations, thus enhancing reproducibility in this highly medically relevant area.

60 APPLIED LIFE SCIENCES↗

Liouville operators over the Hardy space

The role of Liouville operators in the study of dynamical systems through the use of occupation measures has been an active area of research in control theory over the past decade. This article investigates Liouville operators over the Hardy space, which encode complex ordinary differential equations in an operator over a reproducing kernel Hilbert space.

97 MATHEMATICS AND COMPUTING↗