SEARCH · Search NASA

Results for “Astrophysical & cosmological 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.

388 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗

DOE OSTI2023

Dense nuclear matter equation of state from heavy-ion collisions

The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. Furthermore, this white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Development of Novel Combustion Codes for Supercritical CO2 Combustion: Cooperative Research and Development Final Report, CRADA Number CRD-20-17182

The Allam-Fetvedt Cycle, a novel supercritical CO2 (sCO2) power cycle developed by 8 Rivers Capital, LLC, generates low-cost power from fossil fuels while capturing combustion-CO2. This cycle employs high-temperature oxy-combustion using sCO2 as the working fluid, leading to a challenging environment for materials. The most critical component for reliable deployment of this technology is the combustor-turbine interface. In the proposed partnership between NREL and 8 Rivers, we will develop a simulation-based tool to provide key design inputs for defining the materials and geometries required in this critical area.

09 BIOMASS FUELS↗

DOE OSTI2025

High-Fidelity Building Emulator for Integrated Comfort and Energy Analysis using EnergyPlus and Radiance

The growing need for smart, energy-efficient, and occupant-centric buildings has created a demand for advanced control systems that can optimize building operations to balance energy savings, demand flexibility, and comfort. However, current building energy simulation tools, such as EnergyPlus, have limitations that hinder the development and evaluation of these complex control systems. To address this challenge, we introduce a high-fidelity building emulator that dynamically couples EnergyPlus with Radiance for enhanced daylight modeling. The introduced workflow allows researchers and practitioners to rapidly develop and evaluate innovative control solutions. An example study looking at a south-facing office zone revealed up to 67% deviation in predicted light levels, which can significantly impact building assessment.

Yu, Tammie↗

DOE OSTI2026

Probing the atomic dynamics of ultrafast melting with femtosecond electron diffraction

Melting is an every-day phase transition that is determined by thermodynamic parameters like temperature and pressure. In contrast, ultra-fast melting is governed by the microscopic response to a rapid energy input and, thus, can reveal the strength and dynamics of atomic bonds as well as the energy flow rate to the lattice. Accurately describing these processes remains challenging and requires detailed insights into transient states encountered. Here, we present data from femtosecond electron diffraction measurements that capture the structural evolution of copper during the ultrafast solid-to-liquid phase transformations. At absorbed energy densities 2-4 times the melting threshold, melting begins at the surface slightly below the nominal melting point followed by rapid homogeneous melting throughout the volume. Molecular dynamics simulations reproduce these observations and reveal a weak electron-lattice energy transfer rate for the given experimental conditions. Both simulations and experiments show no indications of rapid lattice collapse when its temperature surpasses proposed limits of superheating, providing evidence that the inherent dynamics limits the speed of disordering in ultrafast melting of metals.

FOS: Physical sciences↗

DOE OSTI2026

Integrated Magnetics for Multiport Inductive and Conductive Power Transfer Architecture for Electric Vehicles

This paper proposes a shared magnetic arrangement for an electric-vehicle charger that supports both conductive and inductive charging while preserving independent power control at the two ports. In the proposed implementation, the high-frequency transformer used for the conductive path and the resonant inductor required by the inductive path are realized in an integrated form. By combining the magnetic structure and primary-side converter hardware, the charger can deliver improved hardware utilization and higher power density without a corresponding increase in cost or implementation complexity. The paper describes the magnetic realization and presents finite-element and circuit-level simulation results to confirm minimal interaction between the two charging paths. The study demonstrates that the proposed charger is a viable option for flexible, high-power EV charging systems.

Mukherjee, Subho [ORNL] (ORCID:0009000672297925)↗

DOE OSTI2024

Highly Ordered Eutectic Mesostructures via Template‐Directed Solidification within Thermally Engineered Templates

Template-directed self-assembly of solidifying eutectics results in emergence of unique microstructures due to diffusion constraints and thermal gradients imposed by the template. Here, the importance of selecting the template material based on its conductivity to control heat transfer between the template and the solidifying eutectic, and thus the thermal gradients near the solidification front, is demonstrated. Simulations elucidate the relationship between the thermal properties of the eutectic and template and the resultant microstructure. The overarching finding is that templates with low thermal conductivities are generally advantageous for forming highly organized microstructures. When electrochemically porosified silicon pillars (thermal conductivity < 0.3 Wm −1 K −1 ) are used as the template into which an AgCl-KCl eutectic is solidified, 99% of the unit cells in the solidified structure exhibit the same pattern. In contrast, when higher thermal conductivity crystalline silicon pillars (≈100 Wm −1 K −1 ) are utilized, the expected pattern is only present in 50% of the unit cells. The thermally engineered template results in mesostructures with tunable optical properties and reflectances nearly identical to the simulated reflectances of perfect structures, indicating highly ordered patterns are formed over large areas. This work highlights the importance of controlling heat flows in template-directed self-assembly of eutectics.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Quantifying System Strength From Grid-Forming Resources Using Frequency Scan Approach: Preprint

Current industry practices for quantifying the system strength contribution from grid-forming (GFM) resources to ensure stability of power systems dominated by inverter-based resources (IBRs) are primarily based on iterative electromagnetic transient (EMT) time-domain simulation studies. While feasible, these approaches are resource-intensive, lack scalability and intuition, and might not evaluate the system strength contribution over the entire frequency range of interest. This paper introduces a novel, frequency-domain approach to quantify system strength support provided by a GFM resource using frequency scans. The proposed method uses transfer functions from the grid voltage magnitude (V) and phase (?), respectively, to the reactive (Q) and active power (P) output of a GFM resource for quantifying its contribution to system strength. These transfer functions provide a direct measure of the ability of a GFM resource to behave as a stiff voltage source behind a reactance over a specified frequency range, enabling robust quantification of its system strength contribution. The key innovation of this work is the development of a frequency domain system strength metric called the dynamic short-circuit ratio (dSCR) that is suitable for IBR-dominated power systems and is directly related with the familiar short circuit ratio (SCR) metric. The new metric, dSCR, enables the assessment of system strength contributions from both synchronous machines and converter-based GFM resources using a unified benchmark, which is not possible with the traditional SCR metric. The paper also demonstrates how impedance scans could identify if an unstable condition observed during weak grid conditions is a result of the lack active or reactive power support or both. By leveraging the proposed frequency-domain dSCR metric for quantifying system strength contribution from GFM IBRs, the paper demonstrates targeted mitigation strategies for weak grid instabilities without resorting to repeated, time-consuming time-domain simulations. The result is a scalable and efficient approach to remediate stability challenges in power systems with high shares of IBRs and accelerating the integration of GFM technologies for system strength support in power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

Complex-Valued Intermolecular Coupling Enables Directional Exciton Transport in Excitonic Circuits

Molecular systems capable of directing the flow of excitons are key to the development and optimization of optoelectronic materials. The transport of excitons across multiple molecules is governed by an intermolecular electronic coupling network. In this article, we consider the effects of complex-valued intermolecular electronic coupling on exciton transport. Here, we present a molecular motif capable of generating complex-valued coupling under excitation with circularly polarized light. We use theoretical modeling and simulation to illustrate how complex coupling can be leveraged to drive the rotational flux of excitons in cyclic molecular networks and direct the exciton population in branched molecular networks.

chromophores↗

DOE OSTI2025

Technology & Design Innovations to Maximize the Reduction Effect on DCFC Unit Cost Economics (Max‐REDUCE)

With the high cost and reliability issues of DCFC infrastructure the project specifically addressed the needs of innovating on the AC to DC power conversion. We proved the integration of a new power module into existing charger units, and that theoretical educational work on single stage energy conversion is commercially viable, and a reduced cost compared to existing methods. Our A sample shown 3 PH 480VAC 60hz to 350VDC 8.5kW power transfer, and our software and simulation proved efficiency greater than 97% from 20kW to > 60kW, over the DC output voltage range 200 – 1200V. The project was stopped before we could implement B sample and execute full power testing, validation and integration into a demonstration charger (budget period 3).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

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

Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling

Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .

Close-clearance impellers↗

DOE OSTI2025

Challenges in Continuous In-Field Critical Current Testing of High-Temperature Superconducting Tapes: Thermal and Mechanical Perspectives

High-temperature superconductors (HTS) are essential for ultra-high-field applications requiring exceptional current-carrying capacity under extreme conditions. However, systematic characterization of critical current in long-length conductors remains challenging due to complex thermal, electromag netic, and mechanical interactions during continuous testing. This study reports the development of a continuous in-field magnetization testing system for position-dependent critical current measurement in HTS tapes at 20 K under 7.5 T fields applied normal to the tape plane, enabling identification of performance-limiting regions that could compromise magnet stability. Here, the system addresses two fundamental challenges inherent to cryogenic reel to-reel testing. First, thermal management requires continuous cooling of a moving conductor to 20 K, achieved through liquid nitrogen precooling combined with a 100 W@20 K Gifford McMahon cryocooler. Second, screening currents in high fields generate Lorentz forces that induce twisting, bowing, and potential delamination. To mitigate these risks, we propose mechanical reinforcement and active current density suppression strategies. Numerical simulations using the stream function formulation reveal four primary failure modes: frictional heating at guide interfaces, unstable equilibria causing deformation, transverse current-induced stresses at guide transitions, and unsupported forces in vertical spans. Our mitigation strategies include PTFE coated guides to minimize friction, spring-loaded stabilization mechanisms to maintain tape alignment, controlled pre-heating using the liquid nitrogen thermal jacket to suppress critical current at stress points, and optimized guide positioning to minimize force accumulation. The experimental system is nearing completion, with testing planned to commence within two months. Preliminary validation at 65 K under 0.5 T demonstrates strong correlation between simulation-predicted mechanical instabilities and observed critical current variations during conductor tran sitions through the measurement region. These findings establish a robust foundation for quality assurance protocols essential to next-generation superconducting magnet applications.

Chen, Siwei [Princeton Plasma Physics Laboratory (↗

DOE OSTI2026

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗

DOE OSTI2025

Topological prethermal strong zero modes on superconducting processors

Abstract Symmetry-protected topological phases 1–4 cannot be described by any local order parameter and are beyond the conventional symmetry-breaking model 5 . They are characterized by topological boundary modes that remain stable under symmetry respecting perturbations 1–4,6–8 . In clean, gapped systems without disorder, the stability of these edge modes is restricted to the zero-temperature manifold; at finite temperatures, interactions with mobile thermal excitations lead to their decay 9–11 . Here we report the observation of a distinct type of topological edge mode 12–14 , which is protected by emergent symmetries and persists across the entire spectrum, in an array of 100 programmable superconducting qubits. Through digital quantum simulation of a one-dimensional disorder-free stabilizer Hamiltonian, we observe robust long-lived topological edge modes over up to 30 cycles for a wide range of initial states. We show that the interaction between these edge modes and bulk excitations can be suppressed by dimerizing the stabilizer strength, leading to an emergent U(1) × U(1) symmetry in the prethermal regime of the system. Furthermore, we exploit these topological edge modes as logical qubits and prepare a logical Bell state, which exhibits persistent coherence, despite the system being disorder-free and at finite temperature. Our results establish a viable digital simulation approach 15–18 to experimentally study topological matter at finite temperature and demonstrate a potential route to construct long-lived, robust boundary qubits in disorder-free systems.

Science & Technology - Other Topics↗

DOE OSTI2026

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Two-mode bosonic quoctit for high energy physics

In this work, we study a two-mode bosonic encoding of a quoctit inside a non-Abelian group-structured constellation of coherent states. This work is motivated by the importance of non-Abelian symmetry in particle physics and the desire to have transversal non-Abelian logical gates. We use the previously developed 2 T constellation of states used to encode a so-called 2 T qutrit. The fidelity of the 2 T quoctit is benchmarked against other bosonic qudits for different noise models and find it compare favorably when power constraints are considered. This paves the way for the construction of higher-dimensional qudits (e.g., a quicosotetrit with 2 T group structure) in bosonic systems with practical applications in quantum simulations of particle physics.

Kürkçüoglu, Doga Murat [Fermilab] (ORCID:000000031↗