Search NASA⌕ Search

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.

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

1,696 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Adaptation of virtual synchronous generators to dynamic conditions in power grids

Virtual synchronous generators (VSGs) are widely adopted as grid-forming controls for inverter-based resources. However, when grid conditions vary significantly as characterized by changes in short-circuit ratio (SCR) and the reactance-to-resistance (X/R) ratio, fixed-gain designs and the commonly used P–Q decoupling assumption can become inaccurate. Such conditions can degrade transient power performance, leading to oscillations, prolonged settling, and overshoot, particularly in stiff-grid operating points. This paper quantifies how grid strength and impedance-dependent coupling affect the active–reactive power dynamics of a conventional VSG over a broad range of SCR and X/R values. An adaptive VSG tuning framework is then developed by combining (i) a coupling-explicit, impedance-parameterized state-space model to enable systematic controller synthesis, (ii) a full-state-feedback law designed via pole placement to meet prescribed damping and settling-time specifications, and (iii) a physics-informed neural network (PINN)–based online grid-impedance estimator that updates controller gains in real time as grid conditions vary. Offline simulations in MATLAB/Simulink and real-time validation on an OPAL-RT platform show that the proposed method preserves consistent damping and settling behavior with reduced overshoot across wide SCR and X/R ranges, compared with fixed-gain VSG baselines.

Adaptive control↗
DOE OSTI2026

Reduced instability growth and improved radiation trapping with optimized shock timing in double-shell inertial confinement fusion capsules

The double shell is a volume-burn inertial confinement fusion concept consisting of two concentric shells: a low-Z outer shell that collides with and transfers momentum to a high-Z inner shell which compresses and heats the thermonuclear fuel. The increased number of capsule interfaces and severe hydrodynamic instability of the high-density pusher during its acceleration phase provide challenges to the success of the double shell. Two-dimensional radiation-hydrodynamics simulations predict the hydrodynamic instability growth on the outer surface of the pusher can be greatly reduced through appropriate timing of two shocks that cross this interface. One of these shocks, unique to multi-shell designs, arises from radiation-driven ablation of the inner shell ahead of the main shock, the second shock of concern. The shock timing is optimized by increasing the thickness of a low-Z tamper layer exterior to the pusher, resulting in only minimal changes to the implosion timing. Reducing the instability growth on the outer surface of the high-Z pusher can dramatically decrease the modulations that feedthrough to the pusher inner surface, improving the efficacy of radiation trapping in the thermonuclear fuel and increasing the predicted thermonuclear yield by ≳20×.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Correlating Impacts of Injected Fuels on Carbon Emissions in Blast Furnace with Computational Fluid Dynamics Modeling

A major challenge for steelmaking is the reduction of CO 2 emissions. In this regard, the blast furnace (BF) is critical due to the high associated CO 2 levels. This investigation assesses the impact of tuyere‐injected fuels on BF CO 2 emissions. Specifically, computational fluid dynamics results obtained previously at Purdue University Northwest are analyzed to obtain CO 2 emissions when natural gas (NG), syngas, hydrogen, or hydrogen/NG are injected. CO 2 emissions are compared with those produced when 95 kg of NG/thm is injected. Among these scenarios, the largest CO 2 reduction occurs when 102 kg of syngas/thm (COG feedstock #1) is injected at 973 K, reducing CO 2 by 190.6 kg thm −1 . The largest CO 2 reduction obtained with NG occurs when 130 kg thm −1 is injected at 600 K, reducing emissions by 65 kg thm −1 . H 2 injection also reduces CO 2 , but requires careful adjusting to reach stable operation. For instance, injecting 35 kg of H 2 /thm reduces CO 2 by 52 kg thm −1 . Increasing gaseous injection rates can significantly reduce CO 2 emissions, with fuel preheating providing an addendum, but high injection rates can lead to unstable operation. Furthermore, results show a correlation between CO 2 emissions and average temperature of shaft region for multiple fuels and injection conditions.

Metallurgy & Metallurgical Engineering↗
DOE OSTI2026

Validating a Dynamic PWR Safety and Security Model?

Nuclear power plants (NPPs) are assessed for safety and security using separate models that cannot capture how an attacker's decisions and a plant's response unfold together in real time, leaving regulators and operators without a complete picture of true plant vulnerability. Traditional probabilistic risk assessment (PRA) methods treat adversarial events as fixed initiators with predetermined outcomes, and are structurally incapable of representing the time-dependent interplay between physical security events, safety system response, and operator mitigative actions. At Idaho National Laboratory (INL), I contributed to the development and validation of Modeling and Analysis for Safety and Security using the Dynamic EMRALD Framework (MASS-DEF). Where static PRA relies on event-tree logic that cannot evolve mid-scenario, MASS-DEF couples a time-dependent dynamic PRA tool EMRALD (Event Modeling Risk Assessment using Linked Diagrams) with attack simulation software, allowing attacker behavior, plant system states, and operator actions to interact across time. My work focused on validating a general Pressurized Water Reactor (PWR) model. I traced model logic against PWR plant to identified errors in logic and confirm accuracy. I then built and tested attack scenarios against a general PWR model to verify that the model produced expected outcomes across all logical pathways. I also contributed a section to a related technical paper applying the same EMRALD platform to radiation dose modeling. Results show that MASS-DEF can quantitatively demonstrate that many plants exceed their regulatory security thresholds. This demonstrated margin provides a technically defensible basis for reducing the number of guards without compromising regulatory compliance. Physical security costs represent roughly 10% of annual operating budgets, making such reductions directly meaningful to INL's mission of sustaining existing commercial NPPs. This internship strengthened my understanding of nuclear systems, probabilistic modeling, and technical writing, and has solidified my pursuit of a career at a national laboratory.

98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL↗
DOE OSTI2026

Frontier Job-Centric Telemetry Dataset

Comprehensive analysis of high-performance computing (HPC) systems requires linking workload execution to system behavior. This kind of analysis is vital for diagnosing performance issues, managing capacity, detecting anomalous workloads, and understanding how applications interact with system hardware. This job-centric telemetry dataset unifies scheduler job records with node-level measurements, enabling direct association between workloads and their corresponding power, thermal, and performance characteristics. It contains sanitized, scheduler related metadata for 152,400 individual jobs that ran on the Frontier supercomputer and ended on selected days throughout 2024 and 2025, a subpopulation of ~6.8% of the total number of allocated jobs with non-zero run time on the system over that same period. Each is linked with files that contain telemetry time series records of the power utilization and temperature behavior of its allocated nodes and their processors during the run time of the job. Where available, a portion of the job files also contain network performance time series. Jobs are sampled from select days that reflect normal levels of user activity and possess job size distributions with large numbers of leadership class jobs (>20% of Frontier nodes). Jobs in this dataset attempt to best represent successful user workflows.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Enhancing Coherence Limits in Superconducting Quantum Systems for Computing and Sensing

This talk will highlight recent efforts at the SQMS Center to develop qudit-based quantum computing architectures using superconducting three-dimensional (3D) cavities, as well as the use of these ultra-coherent cavities for quantum sensing. I will present systematic studies of materials and devices aimed at identifying and mitigating the dominant sources of decoherence—including two-level systems (TLS), quasiparticles, and other noise mechanisms—in both transmons and 3D cavities. These investigations include microwave loss characterization of niobium, tantalum, aluminum, their native oxides, and substrate materials such as silicon and sapphire. By combining measurements on qubits and cavities, we disentangle subsystem-specific loss mechanisms and establish a hierarchy of mitigation strategies, leading to transmon coherence times exceeding one millisecond. I will also discuss studies of quasiparticle dynamics, including quasiparticle bursts observed in qubits operated both above ground and at the Gran Sasso underground laboratory, and the observation that applied magnetic fields can suppress temporal T₁ fluctuations. Building on these advances, we demonstrate a record-coherence two-cell cavity-qudit system with coherence times exceeding 20 milliseconds. Leveraging tunable sideband interactions together with error-resilient protocols, including measurement-based error correction and post-selection, we achieve high-fidelity quantum state control, including the preparation of Fock states up to N=20 with fidelities above 95% and the generation of high-fidelity two-mode entangled states. Finally, I will discuss how these ultra-coherent quantum systems are enabling emerging quantum sensing applications, including searches for dark matter and gravitational waves.

Roy, Tanay [Fermilab] (ORCID:000000019442862X)↗
DOE OSTI2024

Computational modeling of coupled mechanical damage and electrochemistry in ternary oxide composite electrodes

Performance degradation of ternary layered oxide cathodes largely originates from their loss of structural integrity in cyclic usage. Mechanical damage, such as intergranular fracture of the active particles, is not only a mechanical cleavage process but also interferes with electrochemical kinetics such as infiltration of liquid electrolyte, surface corrosion of the constituent primary particles, and may eventually isolate the primary grains from the electron conducting network. Here, in this work, we develop a computational framework that integrates electrochemistry of a LiNi x Mn y Co 1−x−y O 2 (NMC) composite cathode with mechanical damage of the active particles. To fully examine the intricate chemomechanical behavior of the electrode, we evaluate the effects of the anisotropic material properties, the influence of mechanical potential on Li transport, and the concurrent intergranular fracture and electrolyte penetration along the grain boundaries upon multiple cycles. Electrolyte infiltration benefits capacity retention but aggravates further mechanical damage by corrosion. Structural failure mostly occurs in the first charging due to the anisotropic mechanical strain between the primary grains, while the resulting damage remains stable in the later few cycles. The results are consistent with experimental observations and the integration of electrochemistry and mechanical failure enables a step further understanding of the complex mechanism of battery degradation.

Battery degradation↗
DOE OSTI2024

Fluorophilic Sigma(σ)‐Lock Self‐Healable Copolymers

Although F-Containing molecules and macromolecules are often used in molecular biology to increase the binding with Lewis acidic groups by introducing favorable C−F dipoles, there is virtually no experimental evidence and limited understanding of the nature of these interactions, especially their role in synthetic polymeric materials. These studies elucidate the molecular origin of inter- and intra-Chain interactions responsible for self-healing of F-Containing copolymers composed of pentafluorostyrene and n-butyl acrylate units (p(PFS/nBA). Guided by dynamic surface oscillating force (SOF) and spectroscopic measurements supported by molecular dynamics (MD) simulations, these studies show that the reformation of σ-σ orbitals in −C−F of PFS and CH 3 CH 2 − of nBA units enables the recovery of entropic energy via fluorophilic-σ-lock van der Waals forces when PFS/nBA molar ratios are ~50/50. The strength of these interactions determined experimentally for self-healable PFS/nBA compositions is in the order ~0.3 kcal/mol which primarily comes from fluorophilic-σ-lock (~70 %) contributions. These interactions are significantly diminished for non-self-healable counterparts. Strongly polarized −C−F σ orbitals create lateral dipolar forces enhancing the affinity towards −C−H orbitals, facilitating energetically favorable interactions. Entropic recovery driven by non-Covalent bonding offers a valuable tool in designing materials with unique functionalities, particularly self-healable batteries and energy storage devices.

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

Self‐Propelling Macroscale Sheets Powered by Enzyme Pumps

Nanoscale enzymes anchored to surfaces act as chemical pumps by converting chemical energy released from enzymatic reactions into spontaneous fluid flow that propels entrained nano‐ and microparticles. Enzymatic pumps are biocompatible, highly selective, and display unique substrate specificity. Utilizing these pumps to trigger self‐propelled motion on the macroscale has, however, constituted a significant challenge and thus prevented their adaptation in macroscopic fluidic devices and soft robotics. Using experiments and simulations, we herein show that enzymatic pumps can drive centimeter‐scale polymer sheets along directed linear paths and rotational trajectories. In these studies, the sheets are confined to the air/water interface. With the addition of appropriate substrate, the asymmetric enzymatic coating on the sheets induces chemically driven, buoyancy flows that controllably propel the sheet's motion on the air/water interface. The directionality and speed of the motion can be tailored by changing the pattern of the enzymatic coating, type of enzyme, and nature and concentration of the substrate. This work highlights the utility of biocompatible enzymes for generating motion in macroscale fluidic devices and robotics and indicates their potential utility for in vivo applications.

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

Preliminary Assessment for the Electric Load Shifting Potential of Integrating Thermal Energy Storage with Heat Pumps in Residential Buildings in Texas of United States

The widespread adoption of electric-driven heat pumps for heating and cooling is expected to significantly increase electric demand. Cooling electric demand will result in an increase in peak hours electric loads, placing additional strain on the grid during peak hours. This challenge, combined with the current rapidly increasing demand from data centers, exacerbates the electric demand duck curve problem. Integrating thermal energy storage (TES) with heat pump can shift electric use for heating and cooling from peak to off-peak hours of the electric grid, which can help flatten the daily electric demand profile of a building. This paper presents a novel design that integrates heat pump with TES (HP-TES), which uses phase change materials (PCM). Heat pump charges TES by melting or freezing PCM during off-peak hours when there is no thermal demand from the building. TES is then discharged (i.e., by freezing or melting PCM) during peak hours to provide a more favorable heat source or heat sink for the heat pump to meet the thermal demand of the building with lower electricity use than conventional air-source heat pumps. Computer simulations were developed to predict the performance of HP-TES applied to a typical single-family house in the US. The building-level simulation results were scaled up to preliminarily assess the aggregated impacts of deploying HP-TES across all single-family houses in Texas of the United States, including reduction of peak demand of the electric grid.

Anees, Fady [ORNL]↗
DOE OSTI2026

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Block copolymer self-assembly derived mesoporous magnetic materials with three-dimensionally (3D) co-continuous gyroid nanostructure

Magnetic nanomaterials are gaining interest for their many applications in technological areas from information science and computing to next-generation quantum energy materials. While magnetic materials have historically been nanostructured through techniques such as lithography and molecular beam epitaxy, there has recently been growing interest in using soft matter self-assembly. In this work, a triblock terpolymer, poly(isoprene-block-styrene-block-ethylene oxide) (ISO), is used as a structure directing agent for aluminosilicate sol nanoparticles and magnetic material precursors to generate organic–inorganic bulk hybrid films with co-continuous morphology. After thermal processing into mesoporous materials, results from a combination of small angle X-ray scattering (SAXS) and scanning electron microscopy (SEM) are consistent with the double gyroid morphology. Nitrogen sorption measurements reveal a type IV isotherm with H1 hysteresis, and yield a specific surface area of around 200 m 2 g −1 and an average pore size of 23 nm. The magnetization of the mesostructured material as a function of applied field shows magnetic hysteresis and coercivity at 300 K and 10 K. Comparison of magnetic measurements between the mesoporous gyroid and an unstructured bulk magnetic material, derived from the identical inorganic precursors, reveals the structured material exhibits a coercivity of 250 Oe, opposed to 148 Oe for the unstructured at 10 K, and presence of remnant magnetic moment not conventionally found in bulk hematite; both of these properties are attributed to the mesostructure. This scalable route to mesoporous magnetic materials with co-continuous morphologies from block copolymer self-assembly may provide a pathway to advanced magnetic nanomaterials with a range of potential applications.

Chemistry↗
DOE OSTI2026

Optimal experimental design using eigenvalue-based criteria with Pyomo.DoE

New developments in automated optimal experimental design within the PSE+ software ecosystem. Advancements in user experience (to reduce the time taken to perform optimal experiment design) and computational capabilities (allowing more diverse experimental design) are shown with an example relevant to critical minerals and materials. Also, a small tutorial on science-based optimal experimental design and novel contributions therein are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Computational discovery of a tetragonal MBene phase: diamond-shaped boron lattices, viable MAB precursors, and selective CO 2 reduction

MBenes are a class of two-dimensional transition-metal borides derived from layered MAB phases. Here, we use first-principles calculations to predict a new tetragonal MB (1 : 1) MBene phase featuring a diamond-shaped B–B lattice and systematically assess its phase stability, synthetic accessibility, and catalytic performance towards the CO 2 reduction reaction (CO 2 RR). Four tetragonal configurations are benchmarked against known hexagonal and orthorhombic phases, and a multi-tier screening identifies six robust members (CrB, FeB, MoB, WB, IrB, and PtB). Exploration of the I4/mmm MAB space reveals 11 viable precursors, including the experimentally synthesized Ir 2 ZnB 2 , which validates our screening approach. CO 2 RR free-energy diagrams show that CrB favors CH 3 OH while FeB, MoB, WB, IrB, and PtB preferentially yield HCOOH. The corresponding limiting potentials are competitive with, and in some cases superior to, those reported for state-of-the-art orthorhombic and hexagonal MBenes. These findings expand the MBene landscape by establishing a tetragonal phase with structural, synthetic, and catalytic promise.

Lu, Linguo [University of Puerto Rico, San Juan, P↗
DOE OSTI2025

Demonstration Trials of AI/ML Edge+Cloud Suite (CRADA Final Report)

PACE AI and LBNL partnered under this CRADA to test and evaluate the PACE5 edge node prototype, an AI/ML edge and cloud-based suite, at FLEXLAB.The objective of the test was to evaluate the PACE5 edge node prototype's ability to perform demand shed and take to dynamic price signals, and to demonstrate advanced fault detection and microgrid monitoring capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Fingerprinting Uranium Oxides with Electron Energy Loss Spectroscopy Supported by Theoretical Computations

Uranium oxides occur in a variety of phases that differ in crystal structure and uranium oxidation states. Electron energy loss spectroscopy (EELS) is one of the few techniques that has sufficient spatial resolution and sensitivity to electronic structure to distinguish amongst phases at the nanoscale. However, beam-sensitive materials such as uranium oxides are subject to spectral modification due to interactions with the electron beam. Therefore, theory support is essential to reliably exclude the impact of beam damage and generate true reference datasets. Here we use a comparison of theoretical and experimental spectra to probe the impact of beam damage on O K-edge and U N-edge (N6,7 and N4,5) EELS spectra of various single-valent and mixed-valence uranium oxide bulk phases. Using a low-dose experimental set-up, we show that the O K-edge theoretical spectra are in excellent agreement with experiment for both peak positions and relative intensities of respective peaks. In contrast, U N-edge features are less distinguishing due to the partially localized nature of the U 5f orbitals and overlapping multiplet and spin–orbit coupling effects. This work demonstrates that O K-edge EELS is sufficiently diagnostic to distinguish a wide range of uranium oxides and that the experimental approach used here minimizes beam damage and allows valence state discrimination across the U(IV), U(V) and U(VI) series. When combined with imaging modes available in electron mi-croscopy, the work enables detailed investigation and characterization of uranium redox transformations at the nanoscale.

Carbone, Jacopo↗
DOE OSTI2026

CIEPAT for Photovoltaic System Resilience

The Cyber-Informed Engineering Photovoltaic Analysis Tool (CIEPAT) was developed in collaboration with the U.S. Department of Energy's Office of Cybersecurity, Energy Security, and Emergency Response (CESER). This tool is an energy source subcomponent integrated into the CIEMAT ecosystem and is developed to enhance the security and resilience of Photovoltaic installations by incorporating Cyber-Informed Engineering (CIE) principles into the deployment of PV systems.

14 SOLAR ENERGY↗
DOE OSTI2026

Assembly bias and local Primordial non-Gaussianity from DESI DR1 quasars

The analysis of the large-scale clustering of quasars (QSO) observed by the Dark Energy Spectroscopic Instrument (DESI) represents a promising avenue for constraining local Primordial non-Gaussianity (PNG), parameterized by f NL . The signal to be constrained is the scale-dependent bias induced in the 2-point clustering of the considered tracer sample. The resulting constraints on f NL , however, are fully degenerate with the local PNG bias parameter b ϕ , dependent on the assembly bias parameter p. Using IllustrisTNG hydrodynamical simulations, we select a QSO sample reflecting the selection criteria and properties of DESI QSOs, and provide a robust prior for p, and thus for b ϕ , building on the findings of Fondi et al. 2024. We find a distribution with mean p̅ ≃ 1.4 with weak redshift dependence, stable to selection noise and consistent with the expected recent merger history typical of quasar-hosting halos. By comparing with the CAMELS simulations we demonstrate that this prior is robust to astrophysical assumptions and cosmic variance. Finally, applying this prior to the DESI DR1 dataset, we derive updated constraints on local PNG, obtaining f NL = -3.3±9.2.

cosmological parameters from LSS↗