Search NASA⌕ Search

SEARCH · Search NASA

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Nonplanar effects in simulations of laser-driven ejecta microjet experiments

Recent experiments of laser-driven ejecta microjets performed at OMEGA 60 reveal tortuous jets whereby the jets appear to deviate from their initial trajectory as they travel in vacuum. To understand these data, we perform two dimensional numerical simulations, considering different target thicknesses, pressures, and models of the drive conditions. In particular, modeling the finite laser spot size appears essential in reproducing qualitatively the non-planar shock observed in the experiment. Simulations capture jet deflection by accounting for a slight misalignment of the laser pointing with respect to the groove axis along with spatial variation of the laser pulse intensity. The principal physical mechanism appears to be that lateral momentum is imparted by release waves arising from the non-planar drive. The induced off-axis velocity is small in comparison to the jet axial velocity but integrates into a pronounced deflection over the course of the experiment. The analysis of jet axial and lateral mass distributions is found to be reproduced reasonably by the simulations. Simulated radiographs are in qualitative agreement with the experiments, though their differences point to potential shortcomings in modeling strictly three-dimensional experiments using two-dimensional hydrodynamic simulations. The simple analysis is able to explain part of the observed discrepancy in simulated vs experimental jet masses.

42 ENGINEERING↗
DOE OSTI2024

Thermonuclear performance variability near ignition at the National Ignition Facility

We describe our current understanding of the variability and degradation mechanisms observed through a series of five indirectly driven inertial fusion implosions fielded at the National Ignition Facility in the fall of 2021, four of which attempted to reproduce the first experiment to achieve Lawson's criterion for ignition with a thermonuclear yield of 1.35 MJ on August 8, 2021. A large number of absolutely calibrated (imaging, time-resolved, and spectrally resolved) x-ray and neutron diagnostics are fielded on the NIF along multiple lines of sight for each experiment. This allows for a reconstruction of the DT fuel and ablator mix injected into the hotspot around peak burn. We show that nuclear yield variations are well reproduced by numerical modeling when the measured low mode asymmetries and mix mass are included. Furthermore, these observed perturbations during burn are linked to small variations in laser delivery and capsule defects. Stringent specifications are then set to achieve robust ignition with the implosion design studied in this paper.

Divol, L. (ORCID:0000000269699898)↗
DOE OSTI2025

Analysis of the HAE activity in the TJ-II stellarator using a Landau closure model

The aim of this study is to analyze the stability of helical Alfvén eigenmodes (HAEs) in TJ-II discharges and the stabilizing effect of the energetic particles generated by the neutral beam injector (NBI) on pressure gradient-driven modes (PGDMs). HAE and PGDM stability is studied using the linear version of the gyro-fluid code FAR3d and the continuous structure by the STELLGAP code. First, Alfvén eigenmode (AE) and PGDM activity observed in the experiments is reproduced by the simulations, identifying unstable m/n = 4/7 − 2/3 and 7/12 − 5/8 HAEs triggered around ρ = 0.66 showing a frequency of 209 and 204 kHz, respectively, as well as 5/3 PGDM. Next, a parametric study is performed with respect to the thermal ion density and iota profile in the middle-outer plasma region to verify the robustness of the simulation results with respect to the uncertainty of experimental profiles. The analysis confirms that experimental uncertainty does not cause large deviations in the simulation results, showing the destabilization of the same HAEs for all the configurations tested. The simulations also indicate the decay of the 5/3 PGDM growth rate as the energetic particle (EP) population in the plasma increases, consistent with the experiment. Stability analysis of the n = 3, 7, 11, n = 5, 9, 13, n = 6, 10, 14, and n = 8, 12 helical families is performed with respect to the NBI operational regime for different EP energies, β as well as deposition profiles. The most unstable configuration is the radially localized on-axis NBI operation (stiff EP density profile gradients nearby the magnetic axis). Using the simulation model that reproduces the observed Alfvén activity, we extend the study to analyze NBI performance within a theoretical framework. It shows that increasing NBI voltage (which raises EP energy) leads to a degradation in NBI performance for a given power (related to EP β and their density). To achieve better NBI operation, higher voltage must be balanced with lower injection power, ensuring stable AEs while keeping the same EP β.

Plasma confinement↗
DOE OSTI2025

Ponderomotive electron physics captured in a single-fluid extended magnetohydrodynamics model

The ponderomotive force, arising from the interaction between electromagnetic waves and plasma, plays a critical role in laser fusion, astrophysical plasmas, and laser diagnostics. Traditionally, modeling this force requires multi-fluid or particle-in-cell simulations due to its strong coupling to electron-scale dynamics. In this work, we demonstrate that a one-fluid, two-temperature extended magnetohydrodynamics (XMHD) model—augmented with a generalized Ohm's law (GOL) including electron inertia—can accurately reproduce key ponderomotive effects. We derive the ponderomotive force within this framework using a phasor-based approach and then validate its nonlinear manifestations through direct numerical simulations in the PERSEUS code, where steepening and density modulation phenomena typically associated with kinetic-scale models are reproduced. These results and prior work establish XMHD as a robust and efficient alternative for modeling nonlinear laser–plasma dynamics, bridging the gap between ideal MHD and fully kinetic approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sparse optimization of two-dimensional terahertz spectroscopy

Two-dimensional terahertz spectroscopy (2DTS) is a low-frequency analog of two-dimensional optical spectroscopy that is rapidly maturing as a probe of a wide variety of condensed matter systems. However, a persistent problem with 2DTS is the long experimental acquisition times, which prevent its broader adoption. A potential solution, requiring no increase in experimental complexity, is signal reconstruction via compressive sensing. In this work, we apply the sparse exponential mode analysis (SEMA) technique to 2DTS of a cuprate superconductor. We benchmark the performance of the algorithm in reconstructing terahertz nonlinearities and find that SEMA reproduces the asymmetric photon echo line shapes at sampling rates as low as 10%, reaching the reconstruction noise floor at sampling rates beyond 20%–30%. The success of SEMA in reproducing such subtle, asymmetric line shapes confirms compressive sensing as a general method to accelerate 2DTS and multidimensional spectroscopies more broadly.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Development of an ultraprecise glue-free bimorph deformable mirror with a length of 460 mm

Lead zirconate titanate (PZT) bimorph mirrors are essential for performing multiple functions, such as beam shaping, wavefront correction, and dynamic focusing with adjustable beam sizes. They typically consist of a mirror substrate with PZT electrodes bonded to the substrate using a thin epoxy film. However, glue-bonded PZT mirrors exhibit low reproducibility in measurement results. One source of this variability is water absorption by the substrate’s epoxy adhesive during wet fabrication processes, such as elastic emission machining (EEM), which leads to swelling and deformation. Additionally, longer deformable bimorph mirrors require larger PZT electrodes, further complicating their design and fabrication, particularly during bonding and ultraprecise metrology. To address these challenges, silver nanoparticles were employed to bond PZT elements to the silicon (Si) mirror, eliminating the need for epoxy glue. Here, in this study, we developed an inorganic-glue-free PZT bimorph mirror with a length of 460 mm and 28 channels. The glue-free bonding method successfully ensured the high reproducibility of measurements after water immersion. This approach resulted in a shape error of only 0.41 nm root mean square after the figure correction process using an ultraprecise EEM process. Scanning acoustic tomography and basic response function tests confirmed the success of the bonding process, showing no critical voids. Furthermore, the shape changes induced by applying voltages to the PZT electrodes closely matched predictions obtained from simulations. These results demonstrate the reliability and precision of the proposed glue-free bonding technique, paving the way for improved performance of new PZT bimorph mirrors.

Ichii, Yoshio [SLAC National Accelerator Laborator↗
DOE OSTI2026

Advanced flip-coil system for magnetic field integral measurements of insertion devices

A novel flip-coil measurement system has been developed for the National Synchrotron Light Source II (NSLS-II) at Brookhaven National Laboratory. This paper describes the design, implementation, and commissioning of the new measurement bench, highlighting its key features, including improved mechanical stability, advanced data acquisition, and enhanced reproducibility. The system enables precise characterization of field integrals and multipole components, ensuring the optimal performance of Insertion Devices (IDs) before installation in the NSLS-II storage ring. The flip-coil system incorporates an innovative approach to minimize mechanical and electrical errors, which significantly improves the reproducibility of measurements. In addition, the system features a state-of-the-art data acquisition system that enables real-time monitoring and analysis, further enhancing the efficiency and accuracy of the measurement process. Furthermore, preliminary tests have demonstrated that the new system meets the stringent requirements for magnetic field characterization of advanced insertion devices, making it an essential tool for future ID commissioning and quality assurance at NSLS-II.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Global energy balance simulations of shattered pellet injection

The dynamics of shattered pellet injection (SPI) shutdowns are simulated using a time-dependent global energy balance model, based on a modification of the KPRAD framework. The new SPI particle source in the model calculates the ablation of individual pellet fragments that enter the plasma as a temporally resolved plume, thus capturing the effects of earlier fragments on the ablation of those that follow, which has a significant impact on the overall assimilation. Despite the reduced physics and the global averaging of all quantities, results from a large number of DIII-D, KSTAR, and JET experiments are well reproduced, including the plasma cooling timescales, particle assimilation, and current quench (CQ) rates. Cooling timescales and CQ rates are in good agreement for pellets containing as little as ∼1% neon, while particle assimilations are most accurate for neon fractions above ∼15% by number of atoms. Below this, the assimilation tends to be overestimated due to the lack of radial transport in the particle balance, which becomes important in the low-Z limit. Predictive simulations of mixed-composition dual-SPI shutdowns in ITER are compared against those with the 3D non-linear magnetohydrodynamic code JOREK and are found to reproduce overall trends observed in the higher-fidelity modeling across a range of injection scenarios. The general success of the model points to the critical role of energy balance in determining SPI particle assimilation and the subsequent disruption dynamics and highlights the value of these simulations for experimental interpretation and for optimizing the deployment of computationally expensive, higher-fidelity models.

energy balance model↗
DOE OSTI2026

Nonlinear electrical transport near the metal–insulator transition in V 4 O 7 thin films

Nonlinear electrical transport associated with correlated electronic states has been widely investigated in transition-metal oxides near metal–insulator transitions. Here, in this study, we investigate nonlinear transport and threshold switching in sputter-deposited V 4 O 7 thin films grown on fused silica substrates. Temperature-dependent transport measurements show a metal–insulator transition near 240 K with negligible thermal hysteresis, defining the temperature scale that governs the nonlinear electrical response. Current–voltage measurements reveal reproducible bipolar threshold switching over a broad temperature range, with the threshold voltage decreasing systematically as the transition temperature is approached. Electrothermal finite-element simulations reproduce the measured switching characteristics and show that switching occurs when localized Joule heating drives a confined region of the device toward the transition temperature. The simulated peak local temperature at threshold lies close to the intrinsic metal–insulator transition temperature of V 4 O 7 , establishing that the nonlinear electrical response originates from electrothermal feedback acting on the strongly temperature-dependent conductivity near the transition. These results identify V 4 O 7 thin films as a model correlated-oxide system in which volatile threshold switching emerges from proximity to an extended metal–insulator transition.

25 ENERGY STORAGE↗
DOE OSTI2017

Mechanics dictate where and how freshwater planarians fission

Asexual freshwater planarians reproduce by tearing themselves into two pieces by a process called binary fission. The resulting head and tail pieces regenerate within about a week, forming two new worms. Understanding this process of ripping oneself into two parts poses a challenging biomechanical problem. Because planarians stop “doing it” at the slightest disturbance, this remained a centuries-old puzzle. In this paper, we focus on Dugesia japonica fission and show that it proceeds in three stages: a local constriction (“waist formation”), pulsation—which increases waist longitudinal stresses—and transverse rupture. We developed a linear mechanical model with a planarian represented by a thin shell. The model fully captures the pulsation dynamics leading to rupture and reproduces empirical time scales and stresses. It asserts that fission execution is a mechanical process. Moreover, we show that the location of waist formation, and thus fission, is determined by physical constraints. Together, our results demonstrate that where and how a planarian rips itself apart during asexual reproduction can be fully explained through biomechanics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Life history variation is maintained by fitness trade-offs and negative frequency-dependent selection

Significance Many populations and species have different strategies for maximizing reproductive success in changing environments. For example, some species of plants reproduce only once in their lifetime (annuals), whereas others can reproduce over multiple years (perennials). Such variation can also occur in fishes, and we examine how this variation in life history strategies is maintained in a Pacific salmonid species ( Oncorhynchus mykiss ). We find that fish that successfully migrate between the ocean and their freshwater spawning grounds multiple times (repeat spawners) have more than twice the lifetime reproductive success of fish that only spawn once (single spawners). We also find a striking pattern of negative frequency-dependent selection in which older, larger females have higher lifetime fitness when they are less abundant.

Christie, Mark R.↗
DOE OSTI2020

Veridical data science

Building and expanding on principles of statistics, machine learning, and scientific inquiry, we propose the predictability, computability, and stability (PCS) framework for veridical data science. Our framework, composed of both a workflow and documentation, aims to provide responsible, reliable, reproducible, and transparent results across the data science life cycle. The PCS workflow uses predictability as a reality check and considers the importance of computation in data collection/storage and algorithm design. It augments predictability and computability with an overarching stability principle. Stability expands on statistical uncertainty considerations to assess how human judgment calls impact data results through data and model/algorithm perturbations. As part of the PCS workflow, we develop PCS inference procedures, namely PCS perturbation intervals and PCS hypothesis testing, to investigate the stability of data results relative to problem formulation, data cleaning, modeling decisions, and interpretations. We illustrate PCS inference through neuroscience and genomics projects of our own and others. Moreover, we demonstrate its favorable performance over existing methods in terms of receiver operating characteristic (ROC) curves in high-dimensional, sparse linear model simulations, including a wide range of misspecified models. Finally, we propose PCS documentation based on R Markdown or Jupyter Notebook, with publicly available, reproducible codes and narratives to back up human choices made throughout an analysis. The PCS workflow and documentation are demonstrated in a genomics case study available on Zenodo.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Interactions between nonfullerene acceptors lead to unstable ternary organic photovoltaic cells

For organic photovoltaic (OPV) devices to achieve consistent performance and long operational lifetimes, organic semiconductors must be processed with precise control over their purity, composition, and structure. This is particularly important for high volume solar cell manufacturing where control of materials quality has a direct impact on yield and cost. Ternary-blend OPVs containing two acceptor–donor–acceptor (A–D–A)-type nonfullerene acceptors (NFAs) and a donor have proven to be an effective strategy to improve solar spectral coverage and reduce energy losses beyond that of binary-blend OPVs. Here, we show that the purity of such a ternary is compromised during blending to form a homogeneously mixed bulk heterojunction thin film. We find that the impurities originate from end-capping C=C/C=C exchange reactions of A–D–A-type NFAs, and that their presence influences both device reproducibility and long-term reliability. The end-capping exchange results in generation of up to four impurity constituents with strong dipolar character that interfere with the photoinduced charge transfer process, leading to reduced charge generation efficiency, morphological instabilities, and an increased vulnerability to photodegradation. As a consequence, the OPV efficiency falls to less than 65% of its initial value within 265 h when exposed to up to 10 suns intensity illumination. We propose potential molecular design strategies critical to enhancing the reproducibility as well as reliability of ternary OPVs by avoiding end-capping reactions.

14 SOLAR ENERGY↗
DOE OSTI2024

Sea-surface temperature pattern effects have slowed global warming and biased warming-based constraints on climate sensitivity

The observed rate of global warming since the 1970s has been proposed as a strong constraint on equilibrium climate sensitivity (ECS) and transient climate response (TCR)—key metrics of the global climate response to greenhouse-gas forcing. Using CMIP5/6 models, we show that the inter-model relationship between warming and these climate sensitivity metrics (the basis for the constraint) arises from a similarity in transient and equilibrium warming patterns within the models, producing an effective climate sensitivity (EffCS) governing recent warming that is comparable to the value of ECS governing long-term warming under CO 2 forcing. However, CMIP5/6 historical simulations do not reproduce observed warming patterns. When driven by observed patterns, even high ECS models produce low EffCS values consistent with the observed global warming rate. The inability of CMIP5/6 models to reproduce observed warming patterns thus results in a bias in the modeled relationship between recent global warming and climate sensitivity. Correcting for this bias means that observed warming is consistent with wide ranges of ECS and TCR extending to higher values than previously recognized. These findings are corroborated by energy balance model simulations and coupled model (CESM1-CAM5) simulations that better replicate observed patterns via tropospheric wind nudging or Antarctic meltwater fluxes. Because CMIP5/6 models fail to simulate observed warming patterns, proposed warming-based constraints on ECS, TCR, and projected global warming are biased low. The results reinforce recent findings that the unique pattern of observed warming has slowed global-mean warming over recent decades and that how the pattern will evolve in the future represents a major source of uncertainty in climate projections.

58 GEOSCIENCES↗
DOE OSTI2026

Data-driven Mori–Zwanzig modeling of Lagrangian particle dynamics in turbulent flows

The dynamics of Lagrangian particles in turbulence play a crucial role in mixing, transport, and dispersion in complex flows. Their trajectories exhibit highly nontrivial statistical behavior, motivating the development of surrogate models that can reproduce these trajectories without incurring the high computational cost of direct numerical simulations of the full Eulerian field. This task is particularly challenging because reduced-order models typically lack access to the full set of interactions with the underlying turbulent field. Novel data-driven machine learning techniques can be powerful in capturing and reproducing complex statistics of the reduced-order/surrogate dynamics. In this work, we show how one can learn a surrogate dynamical system that is able to evolve a turbulent Lagrangian trajectory in a way that is point-wise accurate for short-time predictions (with respect to Kolmogorov time) and stable and statistically accurate at long times. This approach is based on the Mori–Zwanzig formalism, which prescribes a mathematical decomposition of the full dynamical system into resolved dynamics that depend on the current state and the past history of a reduced set of observables, and the unresolved orthogonal dynamics due to unresolved degrees of freedom of the initial state. We show how by training this reduced order model on a point-wise error metric on short time-prediction, we are able to correctly learn the dynamics of Lagrangian turbulence, such that also the long-time statistical behavior is stably recovered at test time. This opens up a range of applications, for example, for the control of active Lagrangian agents in turbulence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Characterizing Tradeoffs in Memory, Accuracy, and Speed for Chemistry Tabulation Techniques

Chemistry tabulation is a common approach in practical simulations of turbulent combustion at engineering scales. Linear interpolants have traditionally been used for accessing precomputed multidimensional tables but suffer from large memory requirements and discontinuous derivatives. Higher-degree interpolants address some of these restrictions but are similarly limited to relatively low-dimensional tabulation. Artificial neural networks (ANNs) can be used to overcome these limitations but cannot guarantee the same accuracy as interpolants and introduce challenges in reproducibility and reliable training. These challenges are enhanced as the physics complexity to be represented within the tabulation increases. Here, we assess the efficiency, accuracy, and memory requirements of Lagrange polynomials, tensor product B-splines, and ANNs as tabulation strategies. We analyze results in the context of nonadiabatic flamelet modeling where higher dimension counts are necessary. While ANNs do not require structuring of data, providing benefits for complex physics representation, interpolation approaches often rely on some structuring of the table. Interpolation using structured table inputs that are not directly related to the variables transported in a simulation can incur additional query costs. This is demonstrated in the present implementation of heat losses. We show that ANNs, despite being difficult to train and reproduce, can be advantageous for high-dimensional, unstructured datasets relevant to nonadiabatic flamelet models. Furthermore we demonstrate that Lagrange polynomials show significant speedup for similar accuracy compared to B-splines.

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

235 U(n, f) Independent fission product yield and isomeric ratio calculated with the statistical Hauser–Feshbach theory

In this work, we have developed a Hauser–Feshbach fission fragment decay model, HF 3 D, which can be applied to the statistical decay of more than 500 primary fission fragment pairs (1,000 nuclides) produced by the neutron induced fission of 235 U. The HF 3 D model numerically integrates fission fragment yield and their initial excitation energy, spin and parity distributions for these 1,000 nuclides. The fission fragment yield Y ( A ) and the total kinetic energy are model inputs, and we estimate them from available experimental data for the 235 U(n th ,f) system. The model parameters in the statistical decay calculation are adjusted to reproduce some fission observables, such as the neutron emission multiplicity v - , its distribution P ( ν ), and the mass dependence v - ( A ). The calculated fission product yield (FPY) and isomeric ratio (IR) are compared with experimental data. The calculated independent FPY Y I ( A ) at the thermal energy reproduces the experimental data well, while the calculated IRs tend to be lower than the Madland–England model prediction. The HF 3 D model is extended to higher incident neutron energies up to the second-chance fission threshold. We demonstrate for the first time that most of the IRs stay constant, although the production of isomeric state changes as the incident energy increases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Application of the coupled-cluster CC( P ; Q ) approaches to the magnesium dimer

Several coupled-cluster approaches with an approximate treatment of connected triply or quadruply excited clusters, including the CCSD(T), CCSD(2) T , and CR-CC(2,3) corrections to CCSD, the active-space CCSDt and CCSDTq methods, and the CC(t;3) and CC(q;4) corrections to CCSDt and CCSDTq derived using the CC(P;Q) framework, and their CCSDT and CCSDTQ parents are applied to the ground-state potential curve and vibrational term values of the magnesium dimer. The correlation-consistent aug-cc-pV(n+d)Z and aug-cc-pwCVnZ (n = T and Q) basis sets are used. Among the noniterative triples corrections to CCSD, the CR-CC(2,3) approach performs the best, but, in analogy to the previously studied beryllium dimer [I. Magoulas et al., J. Phys. Chem. A 122, 1350 (2018)], the CC(t;3) and CC(q;4) methods outperform other employed approaches in reproducing the CCSDT and CCSDTQ data. Composite calculations combining the nearly all-electron CCSDT and valence CC(q;4) and CCSDTQ computations reproduce the experimentally resolved part of the vibrational spectrum of Mg2 to within ~1 cm -1 when the vibrational spacings are examined, while predicting five extra vibrational states near the dissociation threshold, which have not been experimentally resolved. The corresponding binding energies and equilibrium bond lengths are within a few cm -1 and 0.01 A from their experimentally derived values.

74 ATOMIC AND MOLECULAR PHYSICS↗