SEARCH · Search NASA

Results for “Computational Physics and Simulations”

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.

SourceResearch connectionFrom yearThrough year

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

1,696 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Adsorption, charge transfer and a coverage-driven transition of alkali metals on rutile TiO 2 (110)

The interaction of alkali metals with metal oxide surfaces is central to tuning surface reactivity in heterogeneous catalysis and photocatalysis. Here we present a comprehensive DFT+U study of the adsorption of alkali metals (Li, Na, K, Rb, Cs) on the (110) surface of rutile TiO 2 . At low coverage (θ = 1/8), all alkali metals bind preferentially to bridging oxygen sites with adsorption energies in the range −4.06 to −3.33 eV, transferring nearly one full electron (0.92–0.99 |e|) to the substrate and inducing Ti 4+ → Ti 3+ reduction. The excess charge localizes preferentially at subsurface Ti sites in the form of small polarons. Diffusion barriers indicate facile motion along bridging-oxygen rows, whereas inter-row hopping is strongly hindered. Coverage effects were examined systematically for potassium: adsorption energy and charge transfer per K atom decrease monotonically with increasing θ. Strikingly, a sharp energy discontinuity occurs between θ = 4/8 and θ = 5/8 (ΔE ≈ 1 eV per atom), which we identify as a coverage-driven structural transition arising from steric packing constraints and enhanced K–K electrostatic repulsion once every (1×1) surface cell is occupied. This structural transition perfectly correlates with a dramatic drop in the work function down to an ultra-low minimum of 0.84 eV at θ=5/8, followed by a metallization- driven recovery at higher coverages. Ab initio molecular dynamics simulations confirm zigzag K arrangements at moderate coverage (θ = 1/3), while at high coverage (θ = 2/3) short-range K–K correlations emerge without long-range order. These results provide atomistic insight into the structure–activity relationships underlying alkali promotion effects on oxide-supported catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Experimental investigation of a closed vapour box module for a divertor-like configuration in Magnum-PSI

Efficient management of extreme heat fluxes in the divertor region to extend the lifetime of the components remains a critical challenge for the realization of nuclear fusion-based power plants. Among the alternative concepts explored for the divertor region, the use of liquid metals, particularly lithium, is of interest due its ability to dissipate the incoming plasma heat flux through the vapour shielding effect (VS). In this work, we experimentally investigated a ‘closed’ configuration of a dedicated Vapour Box Module (VBM) in the linear plasma device Magnum-PSI. The goal of the experiments is to simulate the vapour box divertor environment conditions and assess its performance in terms of power mitigation and redistribution and lithium confinement. Initial testing without Li demonstrated the efficacy of a closed VBM structure in inducing detachment via neutral gas accumulation. Apertures which enabled non-condensing gas to be effectively pumped while ensuring lithium condensed on the inner surfaces were therefore added. With a lithium capillary porous structure target used, lithium is directly vaporized by the plasma, forming a dense lithium vapour cloud that interacts with the incoming plasma. This resulted in a significant reduction of the target temperature of at least 48%, together with a temperature locking effect, a phenomenon typically observed in the VS regime. Lithium vapour confinement within the VBM was strongly correlated with the wall temperature. Relatively cold walls promoted Li re-condensation and therefore improved Li confinement, although with the expected trade-off of increased hydrogenic retention on lithium-wetted surfaces. As the wall temperature increased, the confinement efficiency decreased, consistent with reduced Li re-condensation and thermally activated Li–H chemistry and remobilization at the walls. Diagnostic measurements through embedded thermocouples and calorimetry revealed that lithium vaporization and re-condensation processes also playedsignificant roles in plasma power dissipation. The results advance the case for a closed divertor chamber with direct lithium evaporation from the strike-points as a viable method to manage divertor heat fluxes in future fusion reactors.

Romano, Fabio [Dutch Institute for Fundamental Ene↗

DOE OSTI2026

UL 1741SC Conformance Development, Evaluation, and Proposal [Slides]

V2G AC is a type of grid interconnection by mobile battery that transforms an EV into a DER with the potential to stabilize the grid and lower costs for electricity. UL 1741SC is a critical new standard to enable grid interconnection of V2G AC-capable vehicles, by defining requirements for V2G AC-intended EVSEs as a gatekeeper between the grid and the EV to support grid safety, such as through redundant voltage trip (namely oversight). EVSE manufacturers need a clear set of conformance test cases and criteria by which to evaluate their equipment and confirm the compliance with UL 1741SC, so that the EVSEs performs as intended for the grid safety. Hence, NLR identified necessary conformance tests and performed the tests identified to understand the maturity.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2026

In Situ UV–vis Quantification of Chlorobasicity in Binary Molten Chlorides

Optical basicity measurements provide a basis for determining the relative chlorobasicity of a given molten salt composition. This study uses Bi3+ as a spectroscopic probe to investigate optical basicity across three binary molten chloride systems (LiCl-KCl, MgCl2–KCl, and ZnCl2–KCl) through in situ UV–vis spectroscopy measurements of the 1S0→3P1 transition. LiCl-KCl systems exhibit the highest optical basicity with minimal compositional dependence, while both divalent systems show stronger, yet nonlinear, dependence that is reminiscent of acid–base titration curves. The influence of cation electronegativity and polarizing power (Li+ < Mg2+ < Zn2+) directly reflects the basicity curves, with ZnCl2–KCl displaying the steepest gradients and lowest basicity values. Nonlinear dependencies in divalent systems indicate formation of chloride-bridged network structures with coordination states dictated by the amount of KCl in a salt composition, which was further investigated and validated by molecular dynamics simulations of polymeric structures formed in ZnCl2 and MgCl2-based salts using machine learning force fields. The temperature dependence of chlorobasicity in terms of optical basicity was investigated, and a loose, positive correlation between the two was revealed. These findings establish Bi3+ as a sensitive probe for understanding electronic environments in molten chloride systems and provide insights for predicting the chemical behavior of molten salts for use in high-temperature environments.

Buttice, Kailee [Department of Nuclear Engineering↗

DOE OSTI2022

Internet of Things Data Characterization Process: Pattern of Life Behavioral Data Study

The HoneyBee™ TARDIS LDRD team completed a data scoping study that identified the initial processes and procedures to baseline the normal and expected behaviors during operability and interoperability of Internet of Things (IoT) device networks. This research is the initial step in developing a process (or methodology) to inform a much broader information framework incorporating machine learning to determine device pattern-of-life which enables the detection of abnormal IoT behaviors on an individual device, as well as in the context of a larger network.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

From Pollutant Removal to Renewable Energy: MoS2-Enhanced P25-Graphene Photocatalysts for Malathion Degradation and H2 Evolution

The widespread presence of pesticides—especially malathion—in aquatic environments presents a major obstacle to conventional remediation strategies, while the ongoing global energy crisis underscores the urgency of developing renewable energy sources such as hydrogen. In this context, photocatalytic water splitting emerges as a promising approach, though its practical application remains limited by poor charge carrier dynamics and insufficient visible-light utilization. Herein, we report the design and evaluation of a series of TiO2-based ternary nanocomposites comprising commercial P25 TiO2, reduced graphene oxide (rGO), and molybdenum disulfide (MoS2), with MoS2 loadings ranging from 1% to 10% by weight. The photocatalysts were fabricated via a two-step method: hydrothermal integration of rGO into P25 followed by solution-phase self-assembly of exfoliated MoS2 nanosheets. The composites were systematically characterized using X-ray diffraction (XRD), Raman spectroscopy, transmission electron microscopy (TEM), UV-Vis diffuse reflectance spectroscopy (DRS), and photoluminescence (PL) spectroscopy. Photocatalytic activity was assessed through two key applications: the degradation of malathion (20 mg/L) under simulated solar irradiation and hydrogen evolution from water in the presence of sacrificial agents. Quantification was performed using UV-Vis spectroscopy, gas chromatography–mass spectrometry (GC-MS), and thermal conductivity detection (GC-TCD). Results showed that the integration of rGO significantly enhanced surface area and charge mobility, while MoS2 served as an effective co-catalyst, promoting interfacial charge separation and acting as an active site for hydrogen evolution. Nearly complete malathion degradation (~100%) was achieved within two hours, and hydrogen production reached up to 6000 µmol g−1 h−1 under optimal MoS2 loading. Notably, photocatalytic performance declined with higher MoS2 content due to recombination effects. Overall, this work demonstrates the synergistic enhancement provided by rGO and MoS2 in a stable P25-based system and underscores the viability of such ternary nanocomposites for addressing both environmental remediation and sustainable energy conversion challenges.

Chemistry↗

DOE OSTI2026

Neutron detector response modeling in NOvA

Neutrons can present a significant challenge for neutrino experiments in which energy reconstruction is critical. With the ability to escape detection completely and with a weak correlation between their kinetic energy and any eventual energy deposition, it is difficult to fully account for neutrons produced in neutrino interactions. This in turn leads to significant model dependence when evaluating neutron-related systematic uncertainties. The NOvA experiment is a long-baseline neutrino oscillation experiment with a high-statistics sample of antineutrino data collected by its near detector. We report an excess relative to data of simulated neutron candidates with low energy depositions when using standard Geant4 physics lists. The simulation excess is traced to an overabundance of secondary photons produced from interactions of neutrons with kinetic energy greater than \SI{20}{\mega\eV}. Improved agreement with data is obtained by applying the data-driven neutron-on-carbon \menate model for neutrons between \SI{20}{\mega\eV} and ${\sim}$\SI{100}{\mega\eV}. With \menate, the residual oversimulation is more uniform across the calorimetric neutron energy spectrum, suggesting possible overproduction of primary neutrons by the GENIE neutrino interaction generator. These results motivate the adoption of \menate-supplemented Geant4 simulation as the nominal simulation in the production of future \nova simulation.

Abubakar, S.↗

DOE OSTI2026

Boron-Based Neutron Scintillator Screen Characterization with X-Rays and Neutrons

Recent work on boron-based neutron scintillator screens suggests these screens can offer superior performance when compared to commonly used screens. Borated neutron scintillator screens perform well in terms of light output (5-6 times greater than a standard Gadox screen) and detection effi-ciency (larger than standard LiF+ZnS screens). However, previously manu-factured boron-based screens have exhibited non-uniform surface coating and a poor mixture between phosphor and converter particles. The objective of this work was to evaluate newly fabricated scintillator screens to deter-mine if enhanced fabrication methods produced a more homogeneous distribution between neutron converter and scintillation phosphor particles. Uniformity of scintillator material deposition was also inspected. This new iteration of screens appeared more uniform than previous generations with the new coating method improving surface chemistry and scintillator material homogeneity. Additionally, a new methodology for screen characterization, involving the correlation of a neutron image taken with a borated scintillator screen to X-ray computed tomography of that same screen, was demonstrated to elucidate a relationship between scintillator screen thickness and relative light output of the screen under neutron exposure. This method suggested that the ideal thickness of scintillator material was ~150 µm to maximize light output of the screen.

36 - MATERIALS SCIENCE↗

DOE OSTI2026

The Hawaii Meteorology, Energy and Transmission (MET) Toolkit

Reliable long-term resource adequacy and grid planning in Hawaii require precise, multi-decadal meteorological records. This paper introduces the Hawaii Meteorology, Energy and Transmission (MET) Toolkit, a 26-year (2000-2025) high-fidelity atmospheric dataset developed by the National Laboratory of the Rockies (NLR). We present a validation study of the underlying WRF model configurations, comparing the legacy MYNN PBL scheme against an alternative YSU formulation. Using vertical lidar profiles and surface buoy data, our analysis identified a foundational geometric distortion in the legacy NOW-23 Hawaii dataset caused by an incorrect grid projection. When evaluated on a corrected, zero-distortion grid, the YSU scheme demonstrated superior performance in bias and cRMSE compared to the legacy setup. To address these findings, the MET Toolkit has been re-produced as a unified 26-year record using the optimized YSU setup and corrected geometry. This dataset offers the Hawaii power sector a robust, validated, and homogenous reference for future grid resilience and energy integration.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

Design and Modeling of the Charge Readout of a SiMOS Quantum Dot with a Single Electron Transistor and CryoCMOS

Single electron spin qubits trapped in SiMOS quantum dots are a promising technology for scaling to thousand- or million-qubit systems due to their compatibility with mature CMOS manufacturing processes. A readout system that combines a single electron transistor with a custom cryogenic CMOS amplification and digitization chain offers key advantages by avoiding the use of bulky RF components or room-temperature interconnects. We present design techniques and simulation results for an optimized qubit-SET cryoCMOS interface, culminating in the design of the QNDR1 ASIC, the first cryogenic readout ASIC designed under the Quandarum project, which targets the development of a many-channel spin qubit based detector for use in high energy physics.

Quinn, Adam [Fermilab] (ORCID:0009000605056371)↗

DOE OSTI2026

Discovery of magnetic-field-tunable density modulations and spin tilting in a layered altermagnet

Altermagnets recently emerged as a new class of magnetic materials, arising from specific spin crystal symmetries. They exhibit a spin-polarized electronic band structure similar to ferromagnets, yet possess zero net magnetization, promising exotic properties. Here we study a layered triangular lattice altermagnet, cobalt-intercalated NbSe 2 using scanning tunneling microscopy and spectroscopy (STM/S). Spectroscopic-imaging STM and spin-polarized STM reveals emergent 2 a 0 tri-directional charge and spin density modulations on the selenium surface. Density functional theory simulations suggest these modulations reflect the underlying cobalt superstructure. We discover that an out-of-plane magnetic field tunes the modulation amplitudes and the electronic density-of-states in a manner dependent on the field direction and strength. This behavior is attributed to the field-induced tilting of cobalt spins, which can have profound implications on the electronic properties of the altermagnet. Our results provide atomic-scale insights to uncover a magnetic-field tunable altermagnetic band structure, highlight the importance of understanding spin canting in altermagnets.

condensed-matter physics↗

DOE OSTI2026

Modeling Radiolysis and Chemical Reactions during Dry Storage of Aluminum-clad Spent Nuclear Fuel

After aluminum-clad spent nuclear fuel (ASNF) is removed from the reactor, it is initially stored in spent fuel pools, which are specially designed water-filled basins that provide temporary cooling to reduce the temperature of the fuel assemblies and provide radiation shielding. ASNF continues to generate heat due to the radioactive decay of elements within the fuel, which persists for many years post-shutdown as the residual radioactive products decay into more stable elements. During the wet storage period, an oxyhydroxide layer composed of boehmite/bayerite forms on the surfaces of the aluminum cladding from exposure to water in the pools. Road-ready packaging for long-term disposition of the ASNF involves dry storage in helium backfilled DOE standard canisters (DSCs). When the ASNF is removed from water storage and dried, most of the water is removed, but some physisorbed and chemisorbed water remains in the oxyhydroxide layers. This residual water can produce hydrogen when exposed to radiation from the ASNF during dry storage. Predicting hydrogen accumulation over time in the DSCs is critical for long-term storage considerations. Previous modeling efforts have developed coupled computational fluid dynamics (CFD)-chemical models to simulate temperature, pressure, and gas phase concentrations within the DSCs. These models use the thermal field predicted by CFD as input to a radiolysis model for the gas phase and the surface oxyhydroxide layer chemistry. Given the long storage period of the DSCs and the impracticality of long-term experiments, a simulation-based approach is necessary to assess chemical evolution within the canisters. This study advances the development of a modeling framework designed to simulate the chemical evolution of spent fuel canisters. Both thermal and radiation-driven reactions are considered, with radiation kinetics quantified using G-values. Sensitivity analysis identifies key parameters influencing species composition. Reaction pathway diagrams offer insight into dominant species formation routes, enabling more effective comparisons between model predictions and experimental observations, particularly regarding the production of hydrogen. Results show that the model predicts significant hydrogen gas production with minimal oxygen generation, primarily due to hydrogen formation via boehmite pathways. These findings underscore the importance of accurately characterizing surface-bound species and radiolysis kinetics. A deeper understanding of these mechanisms is critical for evaluating the long-term safety of nuclear waste storage.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

DOE OSTI2026

Polarized target nuclear magnetic resonance measurements with deep neural networks

Continuous-wave Nuclear Magnetic Resonance (CW-NMR) operated in constant-current mode has served as a foundational technique for polarization measurement in solid-state dynamically polarized targets within nuclear and high-energy physics experiments for several decades, and it remains an essential tool. Conventional Q-meter-based phase-sensitive detection is critical for precise real-time determination of target polarization during scattering runs. However, the accuracy and reliability of these measurements are frequently compromised by elevated noise levels, baseline drift, and systematic uncertainties arising from signal isolation and fitting, ultimately degrading the overall experimental figure of merit. In this work, we report the first successful application of neural network architectures to continuous-wave NMR polarization metrology. By leveraging advanced machine learning techniques for signal extraction and denoising, we achieve a substantial reduction of fitting uncertainties under a variety of realistic simulated and experimental conditions. These improvements translate directly into more robust real-time (online) polarization monitoring and higher precision in subsequent offline analysis. By reducing analysis-induced uncertainty, the resulting methodology can improve the effective figure of merit for scattering experiments employing dynamically polarized targets and provides a new toolset for NMR-based polarimetry in high-energy and nuclear physics.

Metrology↗

DOE OSTI2024

Circumventing data imbalance in magnetic ground state data for magnetic moment predictions

Abstract Magnetic materials play a crucial role in the transition to more sustainable forms of energy and electric vehicles. There is an anticipated shortage in magnetic materials in the future, and as a result there is an urgent need to discover and design new magnetic materials. Computational magnetic material design using density functional theory is daunting because of the challenge in identifying magnetic ground states from a combinatorially large set of possibilities. Machine learning offers a path forward by enabling efficient surrogate models that can more readily enumerate these states, but there is a dearth of training data available, and what is available tends to be imbalanced with too much non-magnetic data. In this work we show that the discrete and previously tackled data imbalance that exists at the level of the magnetic ordering leads to an imbalanced continuous distribution with many zeros when the data is unraveled at the atomic magnetic moment level, which subsequently leads to models with low accuracy for magnetic properties. We mitigate this by using a two-part model framework. Our scheme is able to classify atoms into magnetic and non-magnetic with an F1 score and Matthew’s correlation coefficient (MCC) of ~91% and then to provide an implicit embedding representation that maps directly onto the magnitude of the magnetic moment with a mean absolute error of 0.1 μ B . Beyond screening for new magnetic materials, we demonstrate an additional practical use case of our scheme: the provision of good initial guesses for magnetic moments in first-principles electronic relaxations. Such initialization is shown to lead to faster convergence to configurations that lie closer to the ground state.

Computer Science↗

DOE OSTI2025

Polaron catastrophe within quantum acoustics

The quantum acoustic framework has recently emerged as a nonperturbative, coherent approach to electron–lattice interactions, uncovering rich physics often obscured by perturbative methods with incoherent scattering events. Here, we model the strongly coupled dynamics of electrons and acoustic lattice vibrations within this framework, representing lattice vibrations as coherent states and electrons as quantum wave packets, in a manner distinctively different from tight-binding or discrete hopping-based approaches. We derive and numerically implement electron backaction on the lattice, providing both visual and quantitative insights into electron wave packet evolution and the formation of acoustic polarons. We investigate polaron binding energies across varying material parameters and compute key observables—including mean square displacement, kinetic energy, potential energy, and vibrational energy—over time. Our findings reveal the conditions that favor polaron formation, which is enhanced by low temperatures, high deformation potential constants, slow sound velocities, and high effective masses. Additionally, we explore the impact of external electric and magnetic fields, showing that while polaron formation remains robust under moderate fields, it is weakly suppressed at higher field strengths. These results deepen our understanding of polaron dynamics and pave the way for future studies into nontrivial transport behavior in quantum materials.

Science & Technology - Other Topics↗

DOE OSTI2026

HTGR Multiphysics Application Drivers FY26 Updates

This report summarizes FY26 progress under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program's high-temperature gas-cooled reactor (HTGR) application driver work, covering a wide range of activities such as code validation and multi-physics code assessment. 1) A detailed SAM model of the High-Temperature Engineering Test Reactor (HTTR) was developed using a unique-block grouping approach, with an extended parallel thermal network method to capture block-to-block conduction and radiation heat transfer, and applied to steady-state simulations of the HTTR 30~MW and 9~MW cases. 2) In another activity, SAM's newly implemented multi-component gas flow model was validated against the Natural convection Shutdown heat removal Test Facility (NSTF) argon ingress experiment, correctly capturing the density-driven suppression and thermal recovery of natural circulation observed when argon is introduced into the air-cooled Reactor Cavity Cooling System (RCCS) loop. 3) For the OECD/NEA High Temperature Test Facility (HTTF) benchmark, we co-led the international benchmark activities as well as the OECD/NEA final benchmark report to be released at the end of this year. 4) Finally, the coupled Griffin-SAM modeling capability for pebble-bed HTGRs was advanced by verifying the Griffin neutronics solution against Serpent Monte Carlo for a realistic non-uniform temperature distribution, resolving several deficiencies in the SAM-to-Griffin temperature transfer scheme, and enabling distinct fuel kernel, moderator, and coolant temperatures for cross section feedback. These new features were demonstrated in a PBR load-following transient.

Lee, Alvin↗

DOE OSTI2024

Chemomechanics in alloy phase stability

We describe a first-principles statistical mechanics method to calculate the free energies of crystalline alloys that depend on temperature, composition, and strain. The approach relies on an extension of the alloy cluster expansion to include an explicit dependence on homogeneous strain in addition to site occupation variables that track the degree of chemical ordering. The method is applied to the Si-Ge binary alloy and is used to calculate free energies that describe phase stability under arbitrary epitaxial constraints. We find that while the incoherent phase diagram (in which coexisting phases are not affected by coherency constraints) hosts a miscibility gap, coherent phase equilibrium predicts ordering and negative enthalpies of mixing. Instead of chemical instability, the chemomechanical free energy exhibits instabilities along directions that couple the composition of the alloy with a volumetric strain order parameter. Furthermore, this has fundamental implications for phase field models of spinodal decomposition as it indicates the importance of gradient energy coefficients that couple gradients in composition with gradients in strain.

Materials Science↗

DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗