Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stefan model”

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 361 records · Page 20

Semimetal–Mott insulator quantum phase transition of the Hubbard model on the honeycomb lattice

We take advantage of recent improvements in the grand canonical hybrid Monte Carlo algorithm, to perform a precision study of the single-particle gap in the hexagonal Hubbard model, with on-site electron-electron interactions. After carefully controlled analyses of the Trotter error, the thermodynamic limit, and finite-size scaling with inverse temperature, we find a critical coupling of U c /κ = 3.834(14) and the critical exponent zν = 1.185(43). Under the assumption that this corresponds to the expected antiferromagnetic Mott transition, we are also able to provide a preliminary estimate β = 1.095(37) for the critical exponent of the order parameter. We consider our findings in view of the SU(2) Gross-Neveu, or chiral Heisenberg, universality class. Here, we also discuss the computational scaling of the hybrid Monte Carlo algorithm, and possible extensions of our work to carbon nanotubes, fullerenes, and topological insulators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Virtual Reality for Shoot/No-Shoot Decision Training in Law Enforcement: A Literature Review and Research Agenda

Virtual reality (VR) can materially improve “shoot / no-shoot” (SNS) training by giving officers realistic, repeatable practice making high-stakes decisions under pressure. Traditional tools—live-fire ranges and video simulators—build basics, but they cannot adapt to each officer in real time or fully mirror the complexity of the field. VR closes that gap by creating immersive scenarios that are safer, more flexible, easier to scale across units, and able to capture objective performance data. SNS decisions are not just about marksmanship; they rely on perception, judgment, memory, and the ability to hold fire when a threat is uncertain. Effective training therefore needs realism, decision complexity, and branching outcomes that reflect the true consequences of choices. These elements strengthen recognition of hostile intent while reducing false positives and building the self-control required in ambiguous situations. VR brings specific advantages: dynamic environments, full-body interaction, and the ability to measure performance with precision—enabling targeted feedback and better transfer of learning to the street. At the same time, responsible deployment must address scenario quality (credible environments and behaviors), lawful decision models, and user wellbeing (appropriate stress levels, comfort, and safety). Sandia’s VIPER Lab is positioned to lead this work. The team combines human-performance science, AI/ML, and VR/AR development with a deep equipment bench (e.g., omnidirectional treadmill, eye-tracking, haptics, multiple HMDs). This ecosystem supports building and validating next-generation SNS training that is immersive, measurable, and trustworthy. Bottom line: Investment in VR-enabled SNS training that blends evidence-based design with careful validation and legal safeguards is expected to pay off in safer, more consistent decision-making and improved community trust, delivered through training that is practical to deploy at scale.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

The FLUKA code: Overview and new developments

The FLUKA Monte Carlo Radiation Transport and Interaction code package is widely used to simulate the interaction of particles with matter in a variety of fields, including high energy physics, space radiation, medical applications, radiation protection and shielding assessments, accelerator studies, astrophysical studies and well logging. This paper gives a brief overview of the FLUKA program and describes recent developments, in particular, improvements in the modelling of particle interactions and transport are described in detail. In addition, an overview of selected applications is given.

Ballarini, Francesca↗

Complete set of bound negative-parity states in the neutron-rich nucleus 18 N

High-resolution γ -ray spectroscopy of N 18 is performed with the Advanced GAmma Tracking Array, following deep-inelastic processes induced by an O 18 beam on a Ta 181 target. Six states are newly identified, which together with the three known excitations exhaust all negative-parity excited states expected in N 18 below the neutron threshold. Spin and parities are proposed for all located states on the basis of decay branchings and comparison with large-scale shell-model calculations performed in the p-sd space, with the YSOX interaction. Of particular interest is the location of the 0 1 - and 1 2 - excitations, which provide strong constrains for cross-shell p-sd matrix elements based on realistic interactions and help to simultaneously reproduce the ground and first-excited states in N 16 and N 18 , for the first time. Overall, understanding the N 18 structure may also have significant impact on neutron-capture cross-section calculations in r-process modeling including light neutron-rich nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Chemical and elemental mapping of spent nuclear fuel sections by soft X-ray spectromicroscopy

Soft X-ray spectromicroscopy at the O K -edge, U N 4,5 -edges and Ce M 4,5 -edges has been performed on focused ion beam sections of spent nuclear fuel for the first time, yielding chemical information on the sub-micrometer scale. To analyze these data, a modification to non-negative matrix factorization (NMF) was developed, in which the data are no longer required to be non-negative, but the non-negativity of the spectral components and fit coefficients is largely preserved. The modified NMF method was utilized at the O K -edge to distinguish between two components, one present in the bulk of the sample similar to UO 2 and one present at the interface of the sample which is a hyperstoichiometric UO 2+ x species. The species maps are consistent with a model of a thin layer of UO 2+ x over the entire sample, which is likely explained by oxidation after focused ion beam (FIB) sectioning. In addition to the uranium oxide bulk of the sample, Ce measurements were also performed to investigate the oxidation state of that fission product, which is the subject of considerable interest. Analysis of the Ce spectra shows that Ce is in a predominantly trivalent state, with a possible contribution from tetravalent Ce. Atom probe analysis was performed to provide confirmation of the presence and localization of Ce in the spent fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Unveiling the underlying interactions in Ta 2 NiSe 5 from photoinduced lifetime change

We present a generic procedure for quantifying the interplay of electronic and lattice degrees of freedom in photodoped insulators through a comparative analysis of theoretical many-body simulations and time- and angle-resolved photoemission spectroscopy (TR-ARPES) of the transient response of the candidate excitonic insulator Ta 2 NiSe 5 . Our analysis demonstrates that the electron-electron interactions dominate the electron-phonon ones. In particular, a detailed analysis of the TR-ARPES spectrum enables a clear separation of the dominant broadening (electronic lifetime) effects from the much smaller band-gap renormalization. Theoretical calculations show that the observed strong spectral broadening arises from the electronic scattering of the photoexcited particle-hole pairs and cannot be accounted for in a model in which electron-phonon interactions are dominant. The competing interactions were quantified using the scaling analysis in the weak fluence regime. We demonstrate that the magnitude of the weaker subdominant band-gap renormalization sensitively depends on the distance from the semiconductor/semimetal transition in the high-temperature state, which could explain the apparent contradictions between various TR-ARPES experiments. The analysis presented here indicates that electron-electron interactions play a vital role (albeit not the sole one) in stabilizing the insulating state, and establishes the comparison between lifetime and gap evolution as an important probe of correlated insulators.

36 MATERIALS SCIENCE↗

The Artificial Scientist: in-Transit Machine Learning of Plasma Simulations

Large-scale simulations or scientific experiments produce petabytes of data per run. This poses massive challenges for I/O and storage when scientific analysis workflows are run manually offline. Unsupervised deep learning-based techniques to extract patterns and non-linear relations from these large amounts of data provide a way to build scientific understanding from raw data, reducing the need for manual pre-selection of analysis steps, but require exascale compute and memory to process the full dataset available. In this paper, we demonstrate a heterogeneous streaming workflow in which plasma simulation data is streamed directly to a Machine Learning (ML) application training a model on the simulation data in-transit, completely circumventing the capacity-constrained filesystem bottleneck. This workflow employs openPMD to provide a high level interface to describe scientific data and also uses ADIOS2, to transfer volumes of data that exceed the capabilities of the filesystem. We employ experience replay to avoid catastrophic forgetting in learning from this non-steady state process in a continual manner and adapt it to improve model convergence while learning in-transit. As a proof-of-concept, we approach the ill-posed inverse problem of predicting particle dynamics from radiation in a particle-incell (PIConGPU) simulation of the Kelvin-Helmholtz instability (KHI). We detail hardware-software co-design challenges as we scale PIConGPU to full Frontier, the Top-1 system as of June 2024 Top500 list.

Kelling, Jeffrey [Helmholtz-Zentrum Dresden Rossen↗

Overview of the MOSAiC expedition: Snow and sea ice

Year-round observations of the physical snow and ice properties and processes that govern the ice pack evolution and its interaction with the atmosphere and the ocean were conducted during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition of the research vessel Polarstern in the Arctic Ocean from October 2019 to September 2020. This work was embedded into the interdisciplinary design of the 5 MOSAiC teams, studying the atmosphere, the sea ice, the ocean, the ecosystem, and biogeochemical processes. The overall aim of the snow and sea ice observations during MOSAiC was to characterize the physical properties of the snow and ice cover comprehensively in the central Arctic over an entire annual cycle. This objective was achieved by detailed observations of physical properties and of energy and mass balance of snow and ice. By studying snow and sea ice dynamics over nested spatial scales from centimeters to tens of kilometers, the variability across scales can be considered. On-ice observations of in situ and remote sensing properties of the different surface types over all seasons will help to improve numerical process and climate models and to establish and validate novel satellite remote sensing methods; the linkages to accompanying airborne measurements, satellite observations, and results of numerical models are discussed. We found large spatial variabilities of snow metamorphism and thermal regimes impacting sea ice growth. We conclude that the highly variable snow cover needs to be considered in more detail (in observations, remote sensing, and models) to better understand snow-related feedback processes. The ice pack revealed rapid transformations and motions along the drift in all seasons. The number of coupled ice–ocean interface processes observed in detail are expected to guide upcoming research with respect to the changing Arctic sea ice.

54 ENVIRONMENTAL SCIENCES↗

Nanoscale Confinement of Photo-Injected Electrons at Hybrid Interfaces

A prerequisite for advancing hybrid solar light harvesting systems is a comprehensive understanding of the spatiotemporal dynamics of photoinduced interfacial charge separation. Here, we demonstrate access to this transient charge redistribution for a model hybrid system of nanoporous zinc oxide (ZnO) and ruthenium bipyridyl chromophores. The site-selective probing of the molecular electron donor and semiconductor acceptor by time-resolved X-ray photoemission provides direct insight into the depth distribution of the photoinjected electrons and their interaction with the local band structure on a nanometer length scale. Our results show that these electrons remain localized within less than 6 nm from the interface, due to enhanced downward band bending by the photoinjected charge carriers. This spatial confinement suggests that light-induced charge generation and transport in nanoscale ZnO photocatalytic devices proceeds predominantly within the defect-rich surface region, which may lead to enhanced surface recombination and explain their lower performance compared to titanium dioxide (TiO 2 )-based systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

New avenues for |∆ B | = 2 processes beyond neutron-antineutron oscillations

We explore baryon-number-violating (|∆ B | = 2) processes beyond the well-known neutron-antineutron ($n - \bar{n}$) oscillations, focusing on the $Λ - \bar{Λ}$ system. The presence of a strange quark in the Λ baryon introduces a new set of six-quark operators roughly of the form (uds) 2 , which are different from the (udd) 2 operators responsible for oscillations. Using the Standard Model Effective Field Theory (SMEFT), we classify all dimension-9 operators that cause |∆ B | = 2 transitions and study their UV completions mediated by exotic scalar fields with trilinear interactions. We demonstrate that in these models, oscillations can occur at tree level, with $n - \bar{n}$ mixing potentially appearing at higher loop levels. We employ a chiral effective theory to constrain the effective mass mixing δm Λ , deriving bounds from current experimental limits on $n - \bar{n}$ oscillations and dinucleon decays such as pp → K + K + . These bounds indicate that $Λ - \bar{Λ}$ oscillations probe a complementary parameter space, sensitive to baryon-number violation at scales up to 10 2 − 10 3 TeV. We show that the existing indirect bounds make it challenging to provide a competitive bound on δm Λ at BESIII.

Baryon/Lepton Number Violation↗

Idealized Simulations of the Tropical Climate and Variability in the Single Column Atmosphere Model (SCAM): Radiative-Convective Equilibrium

To explore the interactions among column processes in the Community Atmosphere Model (CAM), the single-column version of CAM (SCAM) is integrated for 1000 days in radiative-convective equilibrium (RCE) with tropical values of boundary conditions, spanning a parameter or configuration space of model physics versions (v5 vs. v6), vertical resolution (standard and 60 levels), sea surface temperature (SST), and some interpretation-driven experiments. The simulated time-mean climate is reasonable, near observations and RCE of a cyclic cloud-resolving model. Updraft detrainment in the deep convection scheme produces distinctive grid-scale structures in humidity and cloud, which also interact with radiative transfer processes. These grid artifacts average out in multi-column RCE results reported elsewhere, illustrating the nuts-and-bolts interpretability that SCAM adds to the hierarchy of model configurations. Multi-day oscillations of precipitation arise from descent of warm convection-capping layers starting near the tropopause, eventually reset by a burst of convective deepening. Experiments reveal how these oscillations depend critically on an internal parameter that controls the number of neutral buoyancy levels allowed for determining cloud top and computing dilute convective available potential energy in the deep convection scheme, and merely modified a little by disabling cloud-base radiation (heating of cloud base). This strong dependence of transient behavior in 1D on this parameter will be tested in the second part of this work, in which SCAM is coupled to a parameterized dynamics of two-dimensional, linearized gravity wave, and in the 3D simulations in future study.

Meteorology & Atmospheric Sciences↗

2025 Advances in NekRS: Supporting improved performance for nuclear applications

This report presents several 2025 advancements in NekRS, a high-fidelity spectral element CFD code developed at Argonne National Laboratory to support the NEAMS thermal-hydraulics program. The forthcoming v25 release consolidates several of these advances, adding new features for portability across heterogeneous GPU architectures, real-time in situ visualization, improved turbulence modeling, and conjugate heat transfer coupling. Over the past year, NekRS has demonstrated strong scalability and performance on DOE’s leading exascale platforms, including Aurora and Frontier, confirming its readiness for some of the largest and most complex simulations attempted to date. These achievements provide a powerful new platform for high-fidelity data generation, which in turn supports the development and validation of advanced closure models critical for reactor safety and design. Significant algorithmic innovations have also been introduced. A new global runtime h-refinement capability simplifies workflows by reducing mesh preparation burdens and enabling coarse-to-fine restarts. Building on this, a novel multigrid strategy was implemented to accelerate pressure and transport solves at scale, addressing long-standing bottlenecks in exascale CFD. Together, these developments improve both the efficiency and accessibility of high-fidelity simulations for reactor-relevant problems. Collectively, these enhancements represent a major step forward in simulation technology, positioning NekRS as a cornerstone of NEAMS efforts to enable accurate, efficient, and scalable high-fidelity analysis of advanced nuclear systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of High-Temperature Emissivity Measurement Capabilities for Molten Salts at PNNL

In this study, an emissivity measurement system for a liquid molten salt was developed on the bases on the Stefan Boltzmann Law. A prototype device operating in ambient atmosphere was constructed an initial result compared well with available literature. The diffuse emissivity of NaCl-KCl system (NaCl, 51mol% NaCl – 49mol% KCl and KCl) was measured from 550 to 850 °C at a wavelength of 1.3 µm. The system appears to be relatively invariant with composition but a clear trend of decreasing emissivity values from ~0.9 at 550 °C to 0.75 at 850 °C was observed. A schematic for an enclosed system for inert atmosphere and actinide salts was drafted. Several of the components for fabrication of this device, improved laser system, furnace, nickel for capsule construction, have been purchased and are ready for assembly. The new laser system will have an improved measurement range of 275-1000 °C. The clam shell furnace has a maximum temperature of 1200 °C and will support the hermetically sealed inert encapsulation vessels. Nickel plate and tubing will create a relatively corrosion resistance capsule for testing that can withstand temperature up to 1000 °C. Emissivity measurements are important for a larger effort to understand radiative heat transfer within molten salts at elevated temperatures (600-1000 °C). Literature review of optical salt properties and radiative heat transfer models indicated that systems within elevated temperature regimes could see upwards of 20% or more of their total heat transferred through photonic mechanisms. Ultimately, future work will measure the total thermal conductivity and break down the molecular and radiative portions. Thermal conductivity, along with other thermal properties of molten salts, both coolant and fuel, are fundamental to modeling and simulations that support MSR design. The results presented here and ongoing work at PNNL will be useful in selecting optimal salt compositions, designing reactor systems, and optimizing heat transfer within a MSR power plant.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An unwanted guest in the electrochemical oxidation of high-voltage Li-ion battery electrolytes: the life of highly reactive protons

Lithium-ion batteries (LIBs) are central to the urgent societal need to decarbonize both transportation and energy storage on the grid. Unfortunately, despite their attractive energy/power density, as well as high coulombic and energy efficiencies, further improvement of this technology – especially their durability – is desperately needed. To support these efforts, our study focuses on fundamental understanding of the decomposition pathways for LIB electrolytes at the cathode–electrolyte interface (CEI), as the nature of these reactions directly controls the extent to which cell capacity and voltage decays in these systems. In this study, we employ electrochemical methods, coupled with product analysis using NMR spectroscopy and mass spectrometry, to determine the decomposition mechanisms in both model and technologically relevant electrolytes. Remarkably, we discovered the electrochemical formation of protons with high chemical activity, comparable to known superacids, at potentials relevant to practical Li-ion batteries. Their reactivity toward every individual component of the CEI provides a unified thermochemical origin for a myriad of side reactions that are commonly associated with the electrochemical reaction. In particular, electrochemically generated protons react with intact EC molecules to form CO 2 and other short and long chain ethers. They also undergo an acid–base reaction with LiPF 6 , to form the weaker acid HF, and with the cathode active material, leaching transition metals into the electrolyte. Collectively, the results of this study all point to the urgent need to either mitigate this proton formation or introduce benign harvesting additives via new electrolyte design strategies.

Ilic, Stefan [Argonne National Laboratory (ANL), A↗

Modeling approaches for addressing unrelaxable bound constraints with unconstrained optimization methods

We explore novel approaches for solving nonlinear optimization problems with unrelaxable bound constraints, which must be satisfied before the objective function can be evaluated. Our method reformulates the unrelaxable bound-constrained problem as an unconstrained optimization problem that is amenable to existing unconstrained optimization methods. The reformulation relies on a domain warping to form a merit function; the choice of the warping determines the level of exactness with which the unconstrained problem can be used to find solutions to the bound-constrained problem, as well as key properties of the unconstrained formulation such as smoothness. We develop theory when the domain warping is a multioutput sigmoidal warping, and we explore the practical elements of applying unconstrained optimization methods to the formulation. Here, we develop an algorithm that exploits the structure of the sigmoidal warping to guarantee that unconstrained optimization algorithms applied to the merit function will find a stationary point to the desired tolerance.

97 MATHEMATICS AND COMPUTING↗

Developing an Automated Uncertainty Quantification Tool to Improve Watershed-Scale Predictions of Water and Nutrient Cycling

Managing the flow of water, nutrients, and contaminants in watersheds is vital to addressing pressing issues related to water scarcity, access to clean drinking water, energy production, resilience to natural and anthropogenic perturbations, and ecological restoration. Decisions about the management of watersheds critically depend on the accuracy with which the flow of water and chemicals through the watershed can be predicted by computer models. Prediction uncertainty can be reduced by matching the model to data, which are collected in the field at great expense. The contribution of watershed characterization data to reducing uncertainty of relevant model predictions can be evaluated in a so-called data-worth analysis, which provides transparent, quantitative metrics about a data set’s value for the support of relevant watershed management objectives. To achieve this goal, we developed a software package that implements the data-worth analysis approach for use with state-of-the-art watershed models. The purpose of the proposed data-worth analysis is to help decision-makers allocate resources for watershed characterization such that the uncertainty in model predictions can be significantly reduced, which leads to better, more effective management decisions. At the same time, watershed characterization costs can be reduced. The specific technical objectives of this SBIR/STTR Phase II project were to develop a framework and associated software toolsets that implement the uncertainty quantification and data-worth analysis approach for use with state-of-the-art watershed models. This goal was achieved by (A) developing a user-friendly, robust software package that is accessible to a wide audience, including watershed managers, policy-makers, and public stakeholders; (B) by demonstrating application of the prototype on several use cases that are representative of complex watershed management challenges spanning a range of scales and that consider different open-source, DOE-based codes and other modeling platforms; and (C) by gathering information about the needs and requirements from potential users to help guide future developments, ensuring that the final product will be commercially viable. The developed software consists of a graphical user interface that guides the user through a sequence of analysis steps, supported by toolsets that leverage state-of-the-art computational simulation-optimization capabilities. A prototype of the software runs on multiple platforms (PC, Mac, multi-processor Linux environment), is linked to diverse watershed simulators (e.g., ECOSYS, TOUGH2, TOUGHREACT, Amanzi-ATS), performs multiple analysis tasks (predictive simulations, sensitivity analysis, uncertainty analysis, automatic parameter estimation, and data-worth analysis, multicomponent geothermometry), and is readily extensible to include external simulators and analysis tools. The software is being commercialized and will be continually updated to address user needs.

58 GEOSCIENCES↗

Probing Active Sites in Cu x Pd y Cluster Catalysts by Machine-Learning-Assisted X-ray Absorption Spectroscopy

Size-selected clusters are important model catalysts because of their narrow size and compositional distributions, as well as enhanced activity and selectivity in many reactions. Still, their structure-activity relationships are, in general, elusive. The main reason is the difficulty in identifying and quantitatively characterizing the catalytic active site in the clusters when it is confined within subnanometric dimensions and under the continuous structural changes the clusters can undergo in reaction conditions. Using machine learning approaches for analysis of the operando X-ray absorption near-edge structure spectra, we obtained accurate speciation of the Cu x Pd y cluster types during the propane oxidation reaction and the structural information about each type. As a result, we elucidated the information about active species and relative roles of Cu and Pd in the clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerated epigenetic aging as a risk factor for chronic obstructive pulmonary disease and decreased lung function in two prospective cohort studies

Chronic obstructive pulmonary disease (COPD) is a frequent diagnosis in older individuals and contributor to global morbidity and mortality. Given the link between lung disease and aging, we need to understand how molecular indicators of aging relate to lung function and disease. Using data from the population-based KORA (Cooperative Health Research in the Region of Augsburg) surveys, we associated baseline epigenetic (DNA methylation) age acceleration with incident COPD and lung function. Models were adjusted for age, sex, smoking, height, weight, and baseline lung disease as appropriate. Associations were replicated in the Normative Aging Study. Of 770 KORA participants, 131 developed incident COPD over 7 years. Baseline accelerated epigenetic aging was significantly associated with incident COPD. The change in age acceleration (follow-up – baseline) was more strongly associated with COPD than baseline aging alone. The association between the change in age acceleration between baseline and follow-up and incident COPD replicated in the Normative Aging Study. Associations with spirometric lung function parameters were weaker than those with COPD, but a meta-analysis of both cohorts provide suggestive evidence of associations. Accelerated epigenetic aging, both baseline measures and changes over time, may be a risk factor for COPD and reduced lung function.

60 APPLIED LIFE SCIENCES↗