SEARCH · Search NASA

Results for “core simulation”

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.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Galaxy formation with wave/fuzzy dark matter: The core-halo structure and the solitonic imprint

Dark matter-dominated cores have long been claimed for the well-studied local group dwarf galaxies. More recently, extended stellar halos have been uncovered around several of these dwarfs through deeper imaging and spectroscopy. Such core-halo structures (inner flat core and a characteristic r −3 asymptotic outer halo profile) are not a feature of conventional cold dark matter (CDM). In contrast, smooth and prominent dark matter cores are predicted for wave/fuzzy dark matter (ψDM). The question arises as to what extent the visible stellar profiles should reflect this dark matter core structure. Here we compare cosmological hydrodynamical simulations of CDM, “WDM” (model used as a proxy for ψDM) & ψDM, aiming to predict the stellar profiles for these three DM scenarios. We show that cores surrounded by extended halos are distinguishable for ψDM, where the stellar density is enhanced in the core due to the presence of the relatively dense soliton. Our analysis demonstrates that, in our simulations, a distinctive core-halo structure does not appear in the case of CDM in the DM, gas, or stars. Whereas we do find a core-halo transition for DM, gas, and stars for ψDM, and the scale of this transition is in line with the predicted core radius set by the soliton scale anticipated for the adopted boson mass of 2.5×10 −22 eV. The presence of a core-halo structure in the stellar profile for Galaxy 1 for ψDM is visible for the most massive and the first galaxy to form in the simulation. Clearly, further simulations are needed to establish how strict this possible relationship is between the DM and stellar core-halo profile as a potential observational discriminator. Furthermore, we observe the anticipated asymmetry for ψDM due to the soliton's motion (jumping and random walk), a distinctive characteristic not found in the symmetric distributions of stars in the warm and CDM models.

dark matter↗

DOE OSTI2024

Plasma edge and scrape-off layer turbulence in gyrokinetic simulations of negative triangularity plasmas

Gyrokinetic simulations in the long-wavelength or drift-kinetic limit are carried out of DIII-D inner-wall-limited (IWL) plasmas to investigate the effect of triangularity on edge and scrape-off layer (SOL) turbulence. The effect of neutral interactions and triangularity on plasma blobs is explored due to the impact blobs can have in setting the SOL width or introducing impurities through interactions with plasma-facing components. Seeded blob simulations with neutrals in shaped SOL scenarios demonstrate that increasing elongation, triangularity, or Shafranov shift decreases radial blob velocities, but neutral interactions have a minor effect. Fully turbulent simulations of DIII-D IWL plasmas include both open- and closed-field-line regions. The negative triangularity (NT) simulation has lower average core T e , lower normalized T e fluctuations, and lower fluxes, but a greater number of coherent structures (blobs) identified with increased size and velocity, on average. Density and electron temperature profiles are within a factor of 2 of experimental values. In conclusion, the increased trapped electron particle fraction in NT simulations is consistent with previous studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Demonstration of Monte Carlo Reference Equilibrium Core Solutions for Pebble Bed Reactors Using Shift in kugelpy

This report describes how to produce reference solutions for the equilibrium core search of pebble bed reactors using the pyrates Python library. The pyrates library uses the kugelpy methodology to perform running-in calculations to reach an equilibrium core. This approach relies on full-core Monte Carlo neutron transport calculations and the shifting of fuel pebble inventory through the reactor to simulate pebble movement in the core from one step to the next of a running-in scenario. In addition to the use of Serpent for the Monte Carlo neutron transport calculations, as part of this work, updates for the use of the Shift Monte Carlo code were integrated into the official pyrates GitHub repository. Comparisons of running-in simulations for a generic pebble bed reactor produced by using pyrates with Serpent and Shift are shown. Although consistent trends are shown between the use of the two Monte Carlo codes, the limiting factor in either calculation is the computation time due to the need to perform several hundred full-core neutron transport calculations before an equilibrium core is reached.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Demonstration of Monte Carlo Reference Equilibrium Core Solutions for Pebble Bed Reactors Using Shift in kugelpy

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Multiphysics Running-In Simulations for Pebble-Bed Reactors with Griffin

Griffin, a Multiphysics Object-Oriented Simulation Environment (MOOSE)–based application targeting transient modeling of advanced reactors, has been used recently to model pebble-bed reactors (PBRs). The modeling effort has focused thus far on equilibrium core calculations. A new capability to simulate the running-in phase of PBR operation has been added to Griffin. This work demonstrates the new capability with a coupled multiphysics running-in simulation. Griffin computes power densities in the core at each time step of the running-in simulation and passes these to Pronghorn, which models fluid flow and heat transfer to calculate pebble surface temperatures. These surface temperatures are used along with the power densities in a heat conduction model to compute average fuel and moderator temperatures, which are passed back to Griffin and accounted for with temperature-dependent cross sections. This work also describes a novel methodology for determining appropriate pebble feed rates and control rod positioning during the running-in simulation. Furthermore, the RZ-geometry model used in this work requires minimal computational resources and can be used for optimization and uncertainty studies in future works.

Griffin↗

DOE OSTI2026

Quantum fluctuations in dense plasma simulations

Molecular dynamics (MD) simulations are a powerful tool for modeling warm and hot dense matter. Density functional theory (DFT) MD simulations are often preferred in dense plasmas in order to accurately model quantum electronic structure. However, DFT-MD simulations neglect interaction effects due to fluctuations in excited states. In this work, we present an MD approach that uses excited state method pseudoatoms to run dense plasma simulations with many different core-electron configurations at classical MD speeds. We also allow for transitions between different configurations in our simulations and find that these fluctuations are especially important for highly excited states. Our results suggest that finite configuration lifetimes that are comparable to the inverse ion plasma frequency need to be accounted for in order to accurately model ion distributions in dense plasma simulations. We also demonstrate that excited state fluctuations have a direct impact on ion plasma microfields, generate different plasma microfields for different excitation levels, and thereby induce absorption–emission line shape asymmetries even in steady-state plasmas.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Stochastic Optimization to Find Optimum Beginning-of-Life Core Configuration of Stable Salt Reactor with Online Refueling

A stochastic optimization method has been developed to find an optimum equilibrium cycle core configuration of the waste-burning stable salt reactor, which is a fast-spectrum molten salt reactor with frequent online refueling. An optimum core configuration was determined with the goal of minimizing radial power peaking. Because of the vast number of potential candidate core configurations, stochastic optimization was applied based on simulated annealing and an additional acceleration method, which screened out unpromising core configurations. It has been demonstrated that the developed stochastic optimization method successfully finds the optimal core configuration regardless of the initial guess and outperforms the gradient descent approach. In addition, it has been observed that the use of a so-called out-in core configuration as the initial guess speeds up convergence of the iterative solution more than five times. Based on the searched optimum equilibrium cycle core configuration, new beginning-of-life (BOL) core configurations have been developed. In conclusion, the new BOL core configurations will be used in developing optimum refueling strategies.

Moltex static salt reactor↗

DOE OSTI2024

Neutrino Trapping and Out-of-Equilibrium Effects in Binary Neutron-Star Merger Remnants

We study out-of-thermodynamic-equilibrium effects in neutron-star mergers with 3D general-relativistic neutrino-radiation large-eddy simulations. During mergers, the cores of the neutron stars remain cold (T ~ a few MeV) and out of thermodynamic equilibrium with trapped neutrinos originating from the hot collisional interface between the stars. However, within ~2 to 3 ms matter and neutrinos reach equilibrium everywhere in the remnant massive neutron star. Furthermore, our results show that dissipative effects, such as bulk viscosity, if present, are only active for a short window of time after the merger.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Braiding of Majorana zero modes in vortex cores

Here, we demonstrate the successful simulation of $\sqrt{𝑍-}$, $\sqrt{𝑋-}$, and 𝑋-quantum gates using Majorana zero modes (MZMs) that emerge in magnetic vortices located in topological superconductors. We compute the transition probabilities and geometric phase differences accounting for the full many-body dynamics and show that qubit states can be read out by fusing the vortex core MZMs and measuring the resulting charge density. We visualize the gate processes using the time- and energy-dependent nonequilibrium local density of states. Our results demonstrate the feasibility of employing vortex core MZMs for the realization of fault-tolerant topological quantum computing.

Majorana bound states↗

DOE OSTI2024

Validation of NSFsim as a Grad-Shafranov equilibrium solver at DIII-D

Plasma shape is a significant factor that must be considered for any Fusion Pilot Plant (FPP) as it has significant consequences for plasma stability and core confinement. A new simulator, NSFsim, has been developed based on a historically successful code, DINA [1], offering tools to simulate both transport and plasma shape. Specifically, NSFsim is a free boundary equilibrium and transport solver and has been configured to match the properties of the DIII-D tokamak. This paper is focused on validating the Grad-Shafranov (GS) solver of NSFsim by analyzing its ability to recreate the plasma shape, the poloidal flux distribution, and the measurements of the simulated diagnostic signals originating from flux loops and magnetic probes in DIII-D. Five different plasma shapes are simulated to show the robustness of NSFsim to different plasma conditions; these shapes are Lower Single Null (LSN), Upper Single Null (USN), Double Null (DN), Inner Wall Limited (IWL), and Negative Triangularity (NT). The NSFsim results are compared against real measured signals, magnetic profile fits from EFIT [2], and another plasma equilibrium simulator, GSevolve [3]. EFIT reconstructions of shots are readily available at DIII-D, but GSevolve was manually ran by us to provide simulation data to compare against.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Kinetic Model of HoxEFU reduction by NADH [SWR-26-087]

This repository is used to release code generated for manuscripts on the Photosynthetic Energy Transduction core program. This code simulates the reduction of HoxEFU by NADH. The electro transfer rate constants for the simulation are specified in the .csv files. The two .csv files correspond tot he two models described in Dawson et al. Cell. Rep. Phys. Sci. 2026. The code utilizes a chemical master equation, a set of differential equations, defining the time evolution of the oxidation and reduction kinetics of NAD+, NADH, a FMN flavin, and a set of iron sulfur clusters. The kinetics of HoxEFU reduction by NADH are evaluated by numerical integration of the chemical master equation using a variable-time-step Runge-Kutta algorithm.

Dahl, Peter [National Laboratory of the Rockies (N↗

DOE OSTI2024

Accelerating CO2 Storage Site Characterization through a New Understanding of Favorable Formation Properties and the Impact of Core-Scale Heterogeneities

CO2 sequestration in deep geologic formations can permanently reduce atmospheric CO2 emissions and help to abate climate change. Target formations must undergo a time- and resource-intensive site evaluation process, assessing storage capacity, environmental safety, and suitability for CO2 trapping via reactive transport models based on data from a limited number of core samples. As such, simulations are often simplified and omit heterogeneities in formation properties that may be significant but are not well understood. To facilitate more rapid site assessment, this work first defines the aquifer properties of favorable storage formations through the analysis of promising and active storage sites. Data show quartz is the most prevalent formation mineral with carbonate minerals, highly reactive with injected CO2, present in over 75% of formations. Porosity and permeability data are highly clustered at 10–30% and 10–1000 mD. Field-scale reactive transport simulations are then constructed and used to analyze CO2 trapping efficiency. The models consider porosity and carbonate mineral heterogeneity as well as the impacts of typical temperature gradients. Simulated sequestration efficiencies are compared to results from a comparable homogenous model to understand the implications of aquifer non-uniformities. The results show a lower sequestration efficiency in the homogeneous model during the injection phase. During the post-injection phase, the homogenization of porosity and carbonate mineralogy results in a higher sequestration efficiency. Incorporating the temperature gradient also increases the sequestration efficiency. Importantly, the maximum deviation between the homogeneous and heterogeneous simulations at the end of the 50-year study period is only ~10%. Larger impacts may be incurred for properties outside the defined, promising ranges suggested here.

Environmental Sciences & Ecology↗

DOE OSTI2026

Mergers Matter: Gravothermal Collapse in Dwarf Halos with Self-Interacting Dark Matter

Self-Interacting Dark Matter (SIDM) models with large cross sections at relative velocities below $\sim100\,{\rm km \, s}^{-1}$ can be tested with dwarf galaxy observations. We analyze six dark-matter-only zoom-in $\sim10^{10}\,{\rm M}_\odot$ halos with diverse assembly histories, adopting a cross section over mass of $σ/m = 70\,cm^2 \, g^{-1}$. We find that mergers inject orbital kinetic energy into the halo, altering the heat transport and the gravothermal evolution of the core. Three of the six halos -- those with the most quiescent merger histories -- show clear signs of core collapse in these simulations. Halos with sustained mergers do not collapse. Furthermore, merger-induced heat transport drives two non-collapsing halos to central densities well below the predictions of the gravothermal fluid model. These findings suggest a novel mechanism for producing dark-matter-deficient galaxies and expanding the diversity of rotation curves beyond what halo concentration alone predicts. Merger histories are thus essential for understanding central density distributions of dwarf galaxy halos in SIDM.

Silverman, Maya [UC, Irvine; Bohr Inst.] (ORCID:00↗

DOE OSTI2024

Effects of stratification on overshooting and waves atop the convective core of M ⊙ main-sequence stars

As a massive star evolves along the main sequence, its core contracts, leaving behind a stable stratification in helium. We simulate two-dimensional convection in the core at three different stages of evolution of a $5\,\mathrm{ M}_{\odot }$ star, with three different stratifications in helium atop the core. We study the propagation of internal gravity waves in the stably stratified envelope, along with the overshooting length of convective plumes above the convective boundary. We find that the stratification in helium in evolved stars hinders radial motions and effectively shields the radiative envelope against plume penetration. This prevents convective overshooting from being an efficient mixing process in the radiative envelope. In addition, internal gravity waves are less excited in evolved models compared to the zero-age-main-sequence model, and are also more damped in the stratified region above the core. As a result, the wave power is several orders of magnitude lower in mid- and terminal-main-sequence models compared to zero-age-main-sequence stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Simulating Muon Capture with MARLEY

Core-collapse supernovae are intense sources of tens-of-MeV neutrinos. Each explosion releases about 1058 neutrinos in a burst lasting tens of seconds. As neutrinos are weakly interacting, they escape more easily compared to photons, allowing them to be observed before the light signal. Understanding low-energy neutrino interactions with argon will be critical for future supernova measurements with the Deep Underground Neutrino Experiment (DUNE). Interpreting future DUNE supernova neutrino data will require a high-quality model of O(10 MeV) neutrino interactions, but the dominant charged-current cross section has never been measured at these low energies for an argon target. We exploit the theoretical similarities between muon capture, which has been measured, and antineutrino scattering to test MARLEY (Model of Argon Reaction Low Energy Yields), the simulation code used by DUNE to model supernova neutrino reactions. I implemented C++ code to add a muon capture mode to MARLEY and compared the predictions of the new simulation to measured isotope branching ratios of the residual nuclei. The preliminary branching ratios predicted by MARLEY running in muon capture mode demonstrate promising baseline performance when compared to experimental data and are expected to improve with further model adjustments.

Wong, Baker Louie↗

DOE OSTI2026

Galaxy Size and Rotation Curve Diversity in ΛCDM with Baryons

The observed rotation curves of dwarf galaxies exhibit significant diversity at fixed halo mass, challenging galaxy formation within the cold dark matter (CDM) model. Previous cosmological galaxy formation simulations with baryonic physics fail to reproduce the full diversity of rotation curves, suggesting that there is a flaw in baryonic feedback models, observational bias, or that an alternative to CDM must be invoked. In this work, we use the Marvelous Massive Dwarf zoom-in simulations, a suite of high-resolution dwarf simulations with M 200 ∼ 10 10 –10 11 M ⊙ and M * ∼ 10 7 –10 9 M ⊙ , designed to target the mass range where the galaxy rotation curve diversity is maximized, i.e., between and 100 km s −1 . We add to this a set of low-mass galaxies from the Marvel Dwarf Zoom Volumes to extend the galaxy mass range to lower values. Our fiducial star formation and feedback models produce simulated dwarfs with a broader range of rotation curve shapes, similar to observations. These simulations both create dark matter cores via baryonic feedback, reproducing the slower-rising rotation curves, while also allowing for compact galaxies and steeply rising rotation curves. Our simulated dwarfs also reproduce the observed size–M * relation, including scatter, producing both extended and compact dwarfs for the first time in simulated field dwarfs. However, the slowly rising and high baryon mass fraction, as well as the steeply rising and low baryon mass fraction, remain missing. We explore star formation and feedback models and conclude that previous simulations may have had feedback that was too strong to produce compact dwarfs.

Cruz, Akaxia [Flatiron Institute, New York, NY (Un↗

DOE OSTI2026

Simulation to a Newborn Supernova Remnant from a Low-mass Iron Core Star

Supernova remnant observations show a high degree of asymmetry, mixing, and inhomogeneity. These asymmetries are seeded during the early seconds of the explosion and are further enhanced and modified as the shock and ejecta move through the stellar progenitor and into the circumstellar medium. We present simulations of a 9.6 M⊙ zero-metallicity progenitor initialized after shock revival and evolved for several years when the ejecta is in the circumstellar medium. A suite of 1D and 2D simulations examines the effects of neutron-star wind and radioactive decay heating. In 1D, decay heating forms a low-density bubble that suppresses the reverse shock. While in 2D, the heating is localized to metal-rich pockets, inflating them and compressing the surrounding material into dense shells. In 3D, the neutron-star wind and decay heating modify the plume morphology, producing more large-scale structures. The extended plume morphology leads to an asymmetrical shock breakout. After breakout, the leading plumes cannot keep up with the shock front, resulting in deceleration and fragmentation by the reverse shock while retaining the large-scale asymmetry. The projected ejecta morphology and velocities are strongly viewing angle dependent. The relatively uniform metal-rich distribution does not resemble the strongly inhomogeneous ejecta structure of Cas A. The 160-isotope decay network shows that 24.4% of the radioactive heating comes from decay chains other than the canonical 56Ni chain. The low explosion energy, low 56Ni yield, and Ni/Fe ratio greater than unity suggest an observational signature similar to an electron capture supernova.

Neopane, Sudarshan [University of Tennessee (UT)]↗

DOE OSTI2026

Enabling Dynamic Probabilistic Risk Assessment of Physical Security Using EMRALD and MAAP (Presentation)

The optimization of physical security in nuclear power plants requires sophisticated methodologies that integrate operator actions and plant behavior through advanced simulation tools. Idaho National Laboratory has developed the Modeling and Analysis for Safety and Security using the Dynamic EMRALD Framework (MASS-DEF) methodology, an approach that integrates force-on-force simulations, dynamic probabilistic risk assessment, and thermal hydraulics modeling to enhance security planning while reducing costs. A reduced order model for thermal hydraulic simulations performed by the Modular Accident Analysis Program (MAAP) was developed to evaluate reactor core behavior during attack scenarios. MAAP simulations are computationally intensive and must be run in a secure environment, complicating analysis and validation. By pre-computed scenario outcomes for a small number of modified parameters, the reduced order model significantly decreases the computational cost and enables offsite review of the results.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗