Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Multimetallic Layered Composites (MMLCs) for Rapid, Economical Advanced Reactor Deployment (Final Report)

This project focused on the development of multi-metallic layered composites (MMLCs) for advanced fission reactor technologies. There are many instances where one alloy or material simply cannot meet all the demands thrown at it by a reactor system, or cannot allow it to perform as strongly as one would like. Instead of focusing all our effort on developing one perfect alloy, we seek to leverage the design principle of “separation of functionality,” used in many other arenas in design, to boost performance beyond single alloys alone. One illustrative example shows the power of this approach for molten salt-cooled reactors: A three meter tall, three meter diameter reactor vessel made of Incoloy 800 was quoted at $\$$500k in 2018. A Hastelloy N vessel was quoted at $\$$5M. An MMLC vessel, in which a layer of Hastelloy N would be weld-overlaid onto Incoloy 800, was quoted at $\$$700k, and it would achieve the same performance. The potential economic gains of leveraging this approach are therefore substantial. At a minimum, each MMLC would contain one core structural layer and one coolant-facing corrosion-resistant layer. Sometimes, MMLCs required buffer layers, as the structural and corrosion-resistant layers were metallurgically incompatible. In other words, they didn’t always play nice, thus separating layers compatible with both functioned as intermediaries to keep the composite together. However, in doing so we inevitably produce new interfaces, where new issues can arise. Therefore, this project focused on what happens at these interfaces from a combination of high temperatures, irradiation, corrosion, and time. After all, a reactor makes money when it is operating, and outages of any kind erode its economic viability. First, we set out to experimentally prove that MMLCs for at least two advanced reactor systems can be made, today, in US domestic facilities. In this respect we were successful – one MMLC (a Ni-201/Incoloy 800H composite) was successfully made and drawn into two-inch coolant piping. Others were attempted, though new issues relating to cracking in vanadium layers for one and radiation damage performance of the corrosion-resistant layer in another prevented us from moving further in those specific arenas – these are engineering problems which deserve continued focus after this project. Additional experimental work focused on long-term corrosion testing of the outermost layers of the salt-cooled and liquid lead-cooled MMLC concepts, which would then be fed into predictions of how long the MMLCs could last. Next, computational (thermodynamics and atomistic) simulation studies studied how much we expect the interfaces to “blend,” due to the mixing action of neutron irradiation. This eats into both the margin for the structural layer of each MMLC, as dilution from the corrosion-resistant layer into the structural layer would decrease the total load-bearing capacity of an MMLC of finite size. On the other hand, dilution of the corrosion-resistant layer into the structural layer further reduced the margin of corrodible material, reducing the lifetime of the MMLC or necessitating extra thickness to be imparted to the MMLC to meet its functional requirements. Work here focused on irradiation-induced segregation to predict new phases which may embrittle the MMLCs, as well as quantifying irradiation-induced mixing at each interface. The results showed that mixing is expected, but it is both steady and therefore predictable, and not lifetime-limiting for most MMLC concepts – it simply has to be accounted for in calculations of reactor performance when utilizing an MMLC. Then, full-core simulations using the experimentally-derived corrosion data, the computationally discovered irradiation-induced mixing data (partially validated by experiment), and existing, benchmarked core designs for large and small sized reactor concepts (one salt-cooled, one lead-cooled) were conducted to quantify any expansion of reactor operating envelopes achieved by utilizing these MMLCs. This new framework, called REX (Reactor Envelope Expansion), incorporates a combination of core neutronics, thermal hydraulics, and the material performance data derived from this project to see how using an MMLC expands advanced fission reactor operating envelopes. It was discovered that in some cases, MMLC utilization does indeed increase the maximum operating temperatures and cycle lengths of reactor concepts, while in other cases it does not. Finally, our tech-to-market (T2M) strategy was not necessarily to create specific embodiments of MMLCs for immediate sale (because getting into the nuclear market is incredibly slow and laden with regulation, this is a long-term goal), but rather immediate stimulation of US industry using the design approach of MMLCs derived from this project. In this respect we were successful, as one of the PhD students funded on this project co-founded Allium Engineering, Inc., which created a stainless steel / low-alloy steel MMLC to function as chloride corrosion-resistant rebar for embedding into concrete structures. Allium Engineering continues to be successful, having recently opened their first factory as of this writing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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

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

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)]↗