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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within ≲3σ of the mean of all simulations on a given grid (128 3 and 256 3 grid cells), where σ describes the statistical variation due to the flow’s time dependence. They also agree well with a 512 3 reference run. The 128 3 and 256 3 simulations agree within 9% and 4%, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Characterization and dynamics of the peculiar stream Jhelum: A tentative role for the Sagittarius dwarf galaxy

Stellar streams are a promising tool to study the Milky Way’s dark matter subhalo population, as interactions with subhalos are expected to leave visible imprints in the streams in the form of substructure. However, there may be other causes of substructure. Here we studied the kinematics and the unusual morphology of the stellar stream Jhelum. Using a combination of ground-based photometry and Gaia EDR3 astrometry, we characterized the morphology of Jhelum. We combined this new data with radial velocities from the literature to perform orbit integrations of the stream in static Galactic potentials. We also carried out N-body simulations in the presence of the Sagittarius dwarf galaxy. The new data reveal a previously unreported tertiary component in the stream, as well as several gaps and a kink-like feature in its narrow component. We find that for a range of realistic Galactic potentials, no single orbit is able to reproduce Jhelum’s radial velocity data entirely. A generic property of the orbital solutions is that they share a similar orbital plane to Sagittarius and this leads to repeated encounters with the stream. Using N-body simulations that include a massive Sagittarius, we explored its effect on Jhelum, and we show that these encounters are able to qualitatively reproduce the narrow and broad components in Jhelum, as well as create a tertiary component in some cases. We also find evidence that such encounters can result in an apparent increase in the velocity dispersion of the stream by a factor up to four due to overlapping narrow and broad components. Our findings suggest that the Jhelum stream is even more complex than once thought; however, its morphology and kinematics can tentatively be explained via the interactions with Sagittarius. In this scenario, the formation of Jhelum’s narrow and broad components occurs naturally, yet some of the smaller gap-like features remain to be explained.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A giant umbrella-like stellar stream around the tidal ring galaxy NGC 922

Tidal ring galaxies (TRGs) are rarely observed in the local universe due to their intrinsically transient nature. The tidal ring structures are the result of strong interactions between gas-rich stellar disks and smaller galactic systems, and do not last longer than ~500 Myr. Therefore, these are perfect scenarios in which to look for the debris of recently accreted dwarf galactic systems. Our goal is to study the low surface brightness stellar structures around the TRG NGC 922 and to revise the hypothesis of its formation in light of these new data. We present new deep images of the TRG NGC 922 and its surroundings from the DESI Legacy survey data and from our observations with an amateur telescope. These observations are compared with results from high-resolution N-body simulations that were designed to reproduce an alternative formation scenario for this peculiar galaxy. Our new observations unveil that the low surface brightness stellar tidal structures around NGC 922 are much more complex than reported in previous works. In particular, the formerly detected tidal spike-like structure at the north-east of the central galaxy disk is not connected with the dwarf companion galaxy PGC 3080368, which has been suggested as the intruder triggering the ring formation of NGC 922. The deep images reveal that this tidal structure mainly has a fainter giant umbrella-like shape, and thus it was formed from the tidal disruption of a different satellite. Using the broad-band g, r, and z DESI LS images, we measured the photometric properties of this stellar stream, estimating a total absolute magnitude in the r band of M r = -17.0 ± 0.03 mag and a total stellar mass for the stream of between 6.9 and 8.5×10 8 M ⊙ . We performed a set of N-body simulations to reproduce the observed NGC922-intruder interaction, suggesting a new scenario for the formation of its tidal ring from the infall of a gas-rich satellite around 150 Myr ago. Finally, our deep images also reveal a tidal shell around the dwarf galaxy PGC 3080368, a possible fossil of a recent merger with a smaller satellite, which may suggest it is in its first infall towards NGC 922.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Anomalous kaon correlations measured in Pb-Pb collisions at the LHC as evidence for the melting and refreezing of the QCD vacuum

Measurements of the dynamical correlations between neutral and charged kaons in central Pb-Pb collisions at $\sqrt{S_{NN}}$ = 2.76 TeV by the ALICE Collaboration display anomalous behavior relative to conventional heavy-ion collision simulators. We consider other conventional statistical models, none of which can reproduce the magnitude and centrality dependence of the correlations. The data can be reproduced by coherent emission from domains which grow in number and volume with increasing centrality. We study the dynamical evolution of the strange quark condensate and show that the energy released during the expansion and cooling of the system may be sufficient to explain the anomaly.

Kapusta, Joseph↗
DOE OSTI2026

Automated ICRF heating surrogate modeling via machine learning

This work introduces automated machine learning workflows that address critical bottlenecks in surrogate model development for Ion Cyclotron Range of Frequencies (ICRF) heating applications. The automated framework includes data analysis tools that transform raw datasets into actionable insights in seconds, replacing weeks of manual exploratory effort and ensuring consistent, reproducible dataset characterization. By integrating advanced hyperparameter optimization (HPO) methods including Bayesian optimization via BoTorch and Tree-structured Parzen Estimators (TPE), the framework significantly reduces model development time from weeks to hours, decreasing computational cost and required expertise, while enabling high-accuracy surrogate models. Compared to traditional hyperparameter scanning (HPS) techniques such as methodical, randomized, and grid searches, HPO methods achieve superior convergence and predictive performance, even when compared to already well-tuned reference models. On NSTX High Harmonic Fast Wave (HHFW) heating datasets, both Random Forest Regressor (RFR) and neural network surrogates demonstrate improved accuracy, achieving R 2 values beyond 0.97 and 0.98, respectively. The results show that while HPO gains are modest for robust architectures like RFR, they become essential for more sensitive models such as neural networks, highlighting the trade-offs across optimization strategies. Through automated workflows that eliminate manual hyperparameter tuning and require minimal ML expertise, this work enables widespread adoption of high-fidelity surrogate models across the fusion community for real-time plasma control, uncertainty quantification, rapid experimental scenario development, and integrated system optimization.

Sanchez-Villar, Alvaro [Princeton Plasma Physics L↗
DOE OSTI2025

Re-evaluating the prompt fission neutron spectrum of spontaneously fissioning 252 Cf

The prompt fission neutron spectrum (PFNS) of spontaneously fissioning 252 Cf is a Neutron Data Standards observable. Nearly all fission spectra of actinides were measured relative to it, using efficiencies derived from it, or analyzed with simulations validated by it. The current Standards evaluation was published by W. Mannhart in 1987. It could not be updated because the evaluation input, experimental mean values and covariances, were lost. First, we attempt to reproduce it. However, Mannhart’s evaluation can only be reproduced within its one-σ uncertainties as some of its aspects (e.g., experimental covariances, rejected data points) remain unknown. Therefore, a new evaluation is presented: We revisit all existing experimental 252 Cf(sf) PFNS data, including those published after the release of the current Standards evaluation, and re-estimate associated covariances. The newly evaluated 252 Cf(sf) PFNS differs distinctly from Mannhart’s below 300 keV and extends it to lower and higher outgoing neutron energies (500 eV–25 MeV). The new evaluated uncertainties are larger from 3–9 MeV and smaller otherwise. Spectrum averaged cross sections of importance to the International Reactor Dosimetry and Fusion File community calculated with the new spectrum are close to those calculated with Mannhart’s evaluation and agree with experimental values well within their uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Automated Approach for Construction of Long-Term, Data-Intensive Watershed Models

Watershed models such as the Hydrological Simulation Program in FORTRAN (HSPF) are frequently used to analyze large-scale water quantity and water quality issues. The construction of HSPF models is a difficult and time-consuming process, because it requires compilation of extensive quantities of data into formatted files followed by a simultaneous fitting of many uncertain parameters. Because of these complications, studies are often limited to calibration periods of only a few years even when analyzing long time-scale issues such as climate and land-use changes that are difficult or impossible to reproduce. High-level, open source programming languages provide an environment for automating the extraction and processing of various sources of data in addition to the calibration process required by HSPF. Recently developed software can be used to build HSPF input files, run simulations, and postprocess simulation results using the Python programming language. The combination of tools in Python, public data sets on the Internet, and software extensions enables rapid development of long-term, reproducible, and sophisticated models for new hydrologic insight. Herein, the utility of these tools is illustrated by developing an automated 30-year HSPF model with a Nash-Sutcliffe efficiency of 0.88 for monthly flows and 0.75 for daily flows in a simulation time of approximately 3h. The integration of HSPF with a high-level programming language creates opportunities to more rigorously explore model assumptions, calibration methods, and alternative data sets.

Lampert, Dave J↗
DOE OSTI2018

Reactive flow modeling of the polymer bonded explosive LX-17 double shock experiments

Overdriven double shock experiments provide a measurement of the properties of the reaction product states of the 1-3-5-triamino-2-4-6trinitrobenzene-based explosive LX-17. These experiments used two flyer materials mounted on the end of a projectile to send an initial shock through the LX-17, followed by a second shock of a higher magnitude into the detonation products. The experimental results are compared to 2D reactive flow modeling. A reactive flow model that describes only the kinetics of the LX-17 decomposition fails to accurately reproduce the decay of the first shock or the curvature or strength of the second shock. A new model is proposed in which the carbon condensate produced in the reaction zone is controlled by a kinetic rate. This allows the carbon condensate to be initially out of chemical equilibrium with the product gas. Finally, this new model reproduces the initial detonation peak and decay and matches the curvature of the second shock; however, it still over-predicts the strength of the second shock.

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

Using the maximum entropy distribution to describe electrons in reconnecting current sheets

Particle distributions in weakly collisional environments such as the magnetosphere have been observed to show deviations from the Maxwellian distribution. These can often be reproduced in kinetic simulations, but fluid models, which are used in global simulations of the magnetosphere, do not necessarily capture any of this. We apply the maximum entropy fluid closure of Levermore, which leads to well posed moment equations, to reconstruct particle distributions from a kinetic simulation in a reconnection region. Our results show that without information other than the moments, the model can reproduce the general structure of the distributions but not all of the finer details. In conclusion, the advantages of the closure over the traditional Grad closure are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Detonation on a tabletop: Nitromethane with high time and space resolution

Nitromethane (NM), a flammable liquid, has been a model system for the shock-to-detonation transition in homogeneous condensed-phase explosives for over 50 years, but we do not understand the fast processes at the molecular scale in the detonation front at the molecular scale. That is largely because prior studies triggered detonations in bomb-sized charges with input shock durations and times-to detonation that were typically microseconds, which made it impossible to observe the faster processes in real time. Here, we studied NM shocked with 4 ns duration input pulses using a tabletop apparatus with laser-launched flyer plates and arrays of tiny disposable optical cuvettes, where the pressure and temperature were probed in real time (1 ns) with photon Doppler velocimetry, optical pyrometry, and high-speed video. Using a 4 ns shock with an input pressure close to the von Neumann spike pressure of 19 GPa, we achieved the minimum time-to-detonation, about 12 ns, where the time-to-detonation is controlled by fundamental molecular processes. We demonstrated the reproducibility of our detonations and showed that they had the same properties as in bomb-sized charges: our detonation velocity, von Neumann spike and Chapman-Jouguet pressures, temperatures, and reaction zone lengths were the same as in bomb-sized charges. Being able to trigger realistic reproducible detonations from a short pulse makes it possible to investigate molecular and fluid dynamics in the detonation by measuring transient responses in real time. Lastly, we found that it took 6 ns for the temperature to reach 3430 K. The high pressure was observed at about 8 ns, when there was a volume explosion to nearly twice the von Neumann spike pressure before settling down to a steady detonation .

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

Gyrokinetic GENE simulations of DIII-D near-edge L-mode plasmas

We present gyrokinetic simulations with the GENE code addressing the near-edge region of an L-mode plasma in the DIII-D tokamak. At radial position ρ = 0.80, simulations with the ion temperature gradient (ITG) increased by 40% above the nominal value give electron and ion heat fluxes that are in simultaneous agreement with the experiment. This gradient increase is consistent with the combined statistical and systematic uncertainty σ of the charge exchange recombination spectroscopy measurements at the 1.6σ level. Multiscale simulations are carried out with a realistic mass ratio and geometry for the first time in the near-edge. These multiscale simulations suggest that the highly unstable ion temperature gradient (ITG) modes of the flux-matched ion-scale simulations suppress electron-scale transport, such that ion-scale simulations are sufficient at this location. At radial position ρ = 0.90, nonlinear simulations show a hybrid state of ITG and trapped electron modes, which was not expected from linear simulations. The nonlinear simulations reproduce the total experimental heat flux with the inclusion of E × B shear effects and an increase in the electron temperature gradient by ∼23%. This gradient increase is compatible with the combined statistical and systematic uncertainty of the Thomson scattering data at the 1.3σ level. These results are consistent with previous findings that gyrokinetic simulations are able to reproduce the experimental heat fluxes by varying input parameters close to their experimental uncertainties, pushing the validation frontier closer to the edge region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Communication: Nucleation rates of supersaturated aqueous NaCl using a polarizable force field

In this work, we use molecular dynamics simulations with a polarizable force field, namely, the modified AH/BK3 model [J. Kolafa, J. Chem. Phys. 145, 204509 (2016)], in combination with the forward flux sampling technique, to calculate the rates of homogeneous nucleation of NaCl from supersaturated aqueous solutions at 298 K and 1 bar. A non-polarizable model that reproduces the experimental equilibrium solubility {AH/TIP4P-2005 of Benavides et al. [J. Chem. Phys. 147, 104501 (2017)]} is also used for comparison. Nucleation rates calculated from the polarizable force field are found to be in good agreement with experimental measurements, while the non-polarizable model severely underestimates the nucleation rates. These results, in combination with our earlier study of a different non-polarizable force field [H. Jiang et al., J. Chem. Phys. 148, 044505 (2018)], lead to the conclusion that nucleation rates are sensitive to the details of force fields, and a good representation of nucleation rates may not be feasible using available non-polarizable force fields, even if these reproduce the equilibrium salt solubility. Finally, inclusion of polarization could be important for an accurate prediction of nucleation rates in salt solutions.

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

New interaction potentials for alkali and alkaline-earth aluminosilicate glasses

Here, we apply a recently developed optimization scheme to obtain effective potentials for alkali and alkaline-earth aluminosilicate glasses that contain lithium, sodium, potassium, or calcium as modifiers. As input data for the optimization, we used the radial distribution functions of the liquid at high temperature generated by means of ab initio molecular dynamics simulations and density and elastic modulus of glass at room temperature from experiments. The new interaction potentials are able to reproduce reliably the structure and various mechanical and vibrational properties over a wide range of compositions for binary silicates. Previously, we tested these potentials for various ternary systems and find that they are transferable and can be mixed, thus allowing us to reproduce and predict the structure and properties of multicomponent glasses.

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

The laser shock station in the dynamic compression sector. I

The Laser Shock Station in the Dynamic Compression Sector [Advanced Photon Source (APS), Argonne National Laboratory] links a laser-driven shock compression platform with high energy x-ray pulses from the APS to achieve in-situ, time-resolved x-ray measurements (diffraction and imaging) in materials subjected to well-characterized, high stress, short duration planar shock waves. Additionally, the combination of this station and the other DCS experimental stations provides a unique and versatile experimental facility to study condensed state phenomena by undertaking real-time, multi-scale measurements in materials subjected to a wide range of shock amplitudes (to above ~350 GPa) and time-durations (~10 ns to 1 μs). The Laser Shock Station uses a 100 joule, 5-17 ns, 351 nm frequency tripled Nd:glass laser with programmable pulse shaping and focal profile smoothing for maximum precision. The laser is capable of a shot every 30 minutes. The interaction chamber contains multiple ports for simultaneous diagnostics, a sample holder that allows sequential exposure for 14 samples without breaking vacuum, the capability to vary the angle between the x-ray beam and laser beam by 135°, and the ability to translate to select one of two types of x-ray beams. X-ray data can be taken with a temporal resolution of ~90 psec. Highly reproducible, planar shocks can be propagated in samples to undertake precise and well-characterized experiments. In a series of ten shots, the absolute variation in the time the shock breaks out at the back of the sample is less than 500ps and the variation in peak particle velocity at the sample/window interface is 4.3%. Lastly, this paper describes the entire DCS Laser Shock Station, including sample fabrication and diagnostics, as well as representative experimental results from shock compressed tantalum to demonstrate the excellent quality of x-ray diffraction data in solid tantalum and to establish the stress for tantalum melting.

47 OTHER INSTRUMENTATION↗
DOE OSTI2019

Dissipative particle dynamics for directed self-assembly of block copolymers

The dissipative particle dynamics (DPD) simulation method has been shown to be a promising tool to study self-assembly of soft matter systems. In particular, it has been used to study block copolymer (BCP) self-assembly. However, previous parameterizations of this model are not able to capture most of the rich phase behaviors of BCPs in thin films nor in directed self-assembly (chemoepitaxy or graphoepitaxy). Here, we extend the applicability of the DPD method for BCPs to make it applicable to thin films and directed self-assembly. Our new reparameterization not only is able to reproduce the bulk phase behavior but also manages to predict thin film structures obtained experimentally from chemoepitaxy or graphoepitaxy. A number of different complex structures, such as bilayer nanomeshes, 90° bend structures, circular cylinders/lamellae and Frank-Kasper phases directed by trenches, and post arrays or chemically patterned substrates, have all been reproduced in this work. This reparameterized DPD model should serves as a powerful tool to predict BCP self-assembly, especially in some complex systems where it is difficult to implement self-consistent field theory.

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

Photoelectron spectroscopy of [Mo6X14]2− dianions (X = Cl–I)

Photoelectron spectroscopy and theoretical investigations have been performed to systematically probe the intrinsic electronic properties of [Mo6X14]2− (X = halogen). All three PE spectra of gaseous [Mo6X14]2− (X = Cl, Br, I) dianions, which were generated by electrospray ionization, exhibit multiple resolved peaks in the recorded binding energy range. Theoretical investigations on the orbital structure and charge distribution were performed to support interpretation of the observed spectra and were further extended onto [Mo6F14]2−, a dianion that was not available for the experimental study. The measured adiabatic (ADE) and vertical detachment energies (VDE) for X = Cl–I were well reproduced by density functional theory calculations (accuracy ∼0.1 eV). Corresponding ADE/VDE values for the dianions were found to be 1.48/2.13 (calc.) and 2.30/2.65, 2.30/2.62, and 2.20/2.42 eV (all expt.) for X = F, Cl, Br, and I, respectively, showing an interesting buckled trend of electron binding energy (EBE) along the halogen series, i.e., EBE (F) ≪ EBE (Cl) ∼ EBE (Br) > EBE (I). Molecular orbital analyses indicate different mixing of metal and halogen atomic orbitals, which is strongly dependent on the nature of X, and suggest that the most loosely bound electrons are detached mainly from the metal core for X = F and Cl, but from halide ligands for X = Br and I. The repulsive Coulomb barrier (RCB), estimated from the photon energy dependent spectra, decreases with increasing halogen size, from 1.8 eV for X = Cl to 1.6 eV for X = I. Electrostatic potential modeling confirms the experimental RCB values and predicts that the most favorable electron detaching pathway should lie via the face-bridging halide ligands.

Li, Ren-Zhong↗
DOE OSTI2020

New interaction potentials for borate glasses with mixed network formers

Here, we adapt and apply a recently developed optimization scheme used to obtain effective potentials for aluminosilicate glasses to include the network former boron into the interaction parameter set. As input data for the optimization, we used the radial distribution functions of the liquid at high temperature generated by ab initio molecular dynamics simulations, and density, coordination, and elastic modulus of glass at room temperature from experiments. The new interaction potentials are shown to reliably reproduce the structure, coordination, and mechanical properties over a wide range of compositions for binary alkali borates. Furthermore, the transferability of these new interaction parameters allows mixing to reliably reproduce the properties of various boroaluminate and borosilicate glasses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗