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 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Predictive modeling of NSTX discharges with the updated multi-mode anomalous transport module

Abstract The objective of this study is twofold: firstly, to demonstrate the consistency between the anomalous transport results produced by updated Multi-Mode Model (MMM) version 9.0.4 and those obtained through gyrokinetic simulations; and secondly, to showcase MMM’s ability to predict electron and ion temperature profiles in low aspect ratio, high beta NSTX discharges. MMM encompasses a range of transport mechanisms driven by electron and ion temperature gradients, trapped electrons, kinetic ballooning, peeling, microtearing, and drift resistive inertial ballooning modes. These modes within MMM are being verified through corresponding gyrokinetic results. The modes that potentially contribute to ion thermal transport are stable in MMM, aligning with both experimental data and findings from linear CGYRO simulations. The isotope effects on these modes are also studied and higher mass is found to be stabilizing, consistent with the experimental trend. The electron thermal power across the flux surface is computed within MMM and compared to experimental measurements and nonlinear CGYRO simulation results. Specifically, the electron temperature gradient modes (ETGM) within MMM account for 2.0 MW of thermal power, consistent with experimental findings. It is noteworthy that the ETGM model requires approximately 5.0 ms of computation time on a standard desktop, while nonlinear CGYRO simulations necessitate 8.0 h on 8 K cores. MMM proves to be highly computationally efficient, a crucial attribute for various applications, including real-time control, tokamak scenario optimization, and uncertainty quantification of experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Core-edge modeling of gyrokinetic turbulence by coupling the delta-f and total-f models in the XGC code

Whole device gyrokinetic turbulence is modeled using a delta-f model in the core coupled to a total-f model in the edge. This new model is implemented in the gyrokinetic code X-Point Included Gyrokinetic Code (XGC) and works within a single simulation. It employs a canonical Maxwellian background in the core and the usual local Maxwellian in the edge. These different core and edge backgrounds are coupled over an overlap region. A thermal-bath operator is used to keep δf small in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Increased Fidelity and Associated Computational cost of Detailed Integral Experiment Benchmarks [Slides]

It does not seem like the system is significantly more sensitive to diameters of components near the center of the core. Intuitively it is, but was not detectable with simulations run to a Monte Carlo k eff uncertainty of 0.00002. The system is more sensitive to heights of components near the center of the core. Most (if not all) Zeus style benchmarks have perturbed core component heights individually.

42 ENGINEERING↗

DOE OSTI2026

An Exploration of the Equation of State Dependence of Core-collapse Supernova Explosion Outcomes and Signatures

We explore, using a state-of-the-art simulation code in 3D and late-enough times to witness final observables, the dependence of core-collapse supernova explosions on the nuclear equation of state (EOS). Going beyond questions of explodability, we compare final explosion energies, nucleosynthetic yields, recoil kicks, and gravitational-wave and neutrino signatures using the SFHo and DD2 nuclear EOSs for a 9 M ⊙ /solar-metallicity progenitor star. The DD2 EOS is stiffer and has a lower effective nucleon mass. The result is a more extended protoneutron star (PNS) and lower central densities. As a consequence, the mean neutrino energies, final explosion energy, and recoil kick speed are lower. Moreover, the evolution of PNS convection differs between the two EOS models in significant ways. This translates in part into interestingly altered neutrino “light” curves and noticeably altered gravitational-wave signal strengths and frequency characteristics that may be diagnostic. The faster exploding model (SFHo) yields slightly more neutron-rich ejecta and more species with atomic weights between 60 and 90 and a weak r-process. However, this is merely a preliminary study. The next step is a more comprehensive and multiprogenitor set of 3D supernova simulations for various EOSs to late times when the observables have asymptoted. Such a future investigation will have a direct bearing on the neutron star and black hole birth mass functions and the quest toward a fully quantitative theory of supernova observables.

Rusakov, Aleksandr [Princeton University, NJ (Unit↗

DOE OSTI2026

Consolidation and Permeability of the B1 and D1 Gas Hydrate Bearing Sands and Associated Seal Sediments of the Extended-Duration Gas Production Test Site on the Alaska North Slope

Gas hydrate, a solid combination of gas (mostly methane in nature) and water molecules stable at low temperatures and elevated pressures, occurs naturally in marine and permafrost-associated environments. Gas hydrate reservoirs, such as those in the Alaska North Slope, have been considered potential energy resources for gas production. To understand the petrophysical and geo-mechanical characteristics of the reservoir, core samples retrieved from the site of the JOGMEC-DOE-USGS collaborative gas hydrate R&D project have been analyzed in the laboratory for their hydraulic and mechanical properties. This paper focuses on both seal and reservoir samples associated with the B1 and D1 sands, which are evaluated for index properties (including porosity, grain size distribution, liquid and plastic limits, specific surface area, and specific gravity), consolidation, permeability, and water retention. Furthermore, the reservoir core samples were tested with pore-filling, laboratory-grown tetrahydrofuran hydrate, in order to assess reservoir behavior during gas production from hydrates. Under simulated in situ stress conditions, the seal and hydrate-free reservoir cores had a permeability anisotropy ratio of k h /k v = 3.0−5.0, and k h /k v = 2.4−3.0 for the reservoir tetrahydrofuran hydrate-bearing cores. The data suggest that depressurizing the reservoir to induce hydrate dissociation alters the reservoir effective permeability in three ways: permeabilities decrease due to porosity lost (e.g., the initial reservoir thickness can decrease by up to 5% upon 7 MPa depressurization), permeability increases due to the loss of solid hydrate in the pore space, and permeability anisotropy k h /k v decreases in response to the evolving pore-space geometry. We show that given the simulated in situ gas hydrate saturations (i.e., S h = 32% in core 7P-2E and S h = 21% in core 20P-4), gas production from the dissociation of tetrahydrofuran hydrate in the two tested cores results in a net increase in effective permeability and a decrease in k h /k v . This study highlights the importance of investigating seal and reservoir sediments and the impacts of depressurization on the porosity and permeability responses during production.

Geological materials↗

DOE OSTI2025

SASHIMI-SIDM: semi-analytical subhalo modelling for self-interacting dark matter at sub-galactic scales

We combine the semi-analytical structure formation model, SASHIMI, which predicts subhalo populations in collisionless, cold dark matter (CDM), with a parametric model that maps CDM halos to self-interacting dark matter (SIDM) halos. The resulting model, SASHIMI-SIDM, generates SIDM subhalo populations down to sub-galactic mass scales, for an arbitrary input cross section, in minutes. We show that SASHIMI-SIDM agrees with SIDM subhalo populations from high-resolution cosmological zoom-in simulations in resolved regimes. Crucially, we predict that the fraction of core-collapsed subhalos peaks at a mass scale determined by the input SIDM cross section and decreases toward higher halo masses, consistent with the predictions of gravothermal models and cosmological simulations. For the first time, we also show that the core-collapsed fraction decreases toward lower halo masses. While the dependence of the collapse time on mass and concentration implies such behaviour, our semi-analytical approach allows us to quantify and illustrate this trend clearly across the full mass spectrum of subhalos, including for subhalo masses below the resolution limit of any current cosmological SIDM simulation. As a proof of principle, we apply SASHIMI-SIDM to predict the boost to the local dark matter density and annihilation rate from core-collapsed SIDM subhalos, which can be enhanced relative to CDM by an order of magnitude for viable SIDM models. Thus, SASHIMI-SIDM provides an efficient and reliable tool for scanning SIDM parameter space and testing it with astrophysical observations. The code is publicly available at https://github.com/shinichiroando/sashimi-si.

cosmological simulations↗

DOE OSTI2025

Challenges and approaches to interpretive modeling of boundary plasma and neutral transport in a closed, pumped divertor

An experimental discharge from the DIII-D tokamak is modeled using the SOLPS-ITER code suite and compared against measurements in the pumped and relatively closed upper divertor. Uncertainties of boundary plasma simulations are identified by attempting to match code inputs to experimental conditions, including iteratively solving transport coefficients to match upstream experimental profiles using varying quantities of core particle flux, different pumping models, and various assumptions of ion thermal transport. Simulated boundary conditions for particle injection at the core interface are shown to be relevant to the plasma solution at the divertor targets, even if upstream transport is modified so that plasma profiles are comparatively similar, although seperatrix density is not held constant. When upstream plasma profiles are matched to experimental measurements by varying diffusive transport coefficients, using either poloidally symmetric or ballooning structure, the model finds a majority of injected energy being transported radially off the computational domain, in conflict with experimental radiated power measurements and heat flux measurements at the divertor target. Imposing a maximum thermal diffusivity or radially shifting the experimental separatrix location of the fitted profiles to increase power conducted to the targets by increasing the upstream electron temperature does not significantly modify this result. Including a thermalizing plenum volume in the simulation domain is shown to maintain the experimental volumetric pumping rate without knowing the neutral energy distribution incident on the pump duct a priori. By modifying transport parameters to match different assumptions for ion temperature, downstream neutral pressure changes by more than a factor of two, suggesting that attention to ion thermal transport may be a critical parameter for simulations to accurately resolve recycling and neutral transport, particularly in a closed divertor geometry. In addition to quantifying various modeling uncertainties, this work motivates both further experimental study and modeling improvements to improve predictive capabilities.

divertor↗

DOE OSTI2024

On the Geometry of the Near-core Magnetic Field in Massive Stars

It is well-known that the cores of massive stars sustain a stellar dynamo with a complex magnetic field configuration. However, the same cannot be said for the field's strength and geometry at the convective–radiative boundary, which are crucial when performing asteroseismic inference. In this Letter, we present 3D magnetohydrodynamic (MHD) simulations of a 7 M ⊙ mid-main-sequence star, with particular attention given to the convective–radiative boundary in the near-core region. Our simulations reveal that the toroidal magnetic field is significantly stronger than the poloidal field in this region, contrary to recent assumptions. Moreover, the rotational shear layer, also important for asteroseismic inference, is specifically confined within the extent of the Brunt–Väisälä frequency peak. These results, which are based on the inferred properties of HD 43317, have widespread implications for asteroseismic studies of rotation, mixing, and magnetism in stars. While we expect our results to be broadly applicable across stars with similar Brunt–Väisälä frequency profiles and stellar masses, we also expect the MHD parameters (e.g., Re m ) and the initial stellar rotation rate to impact the geometry of the field and differential rotation at the convective–radiative interface.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

CaptureFOAM

CaptureFOAM provides utilities for fast and accurate direct air capture (DAC) contactor simulations, based on the OpenFOAM framework. At the core, CaptureFOAM is a suite of mass transfer and chemistry models built upon the OpenFOAM thin film model, which allows coarser mesh resolution and less constrictive time step requirements when compared to a traditional volume of fluid approach. This allows CaptureFOAM to conduct device-scale simulations for complex geometries including commercial packing materials and novel contactor shapes, matching experimental results in a wide variety of applications.

Kincaid, Kellis [Oak Ridge National Laboratory (OR↗

DOE OSTI2024

Initial Demonstration of New Griffin Technologies for Simulating the Running-In Phase of Pebble Bed Reactors

Griffin is a reactor multiphysics modeling application based on MOOSE (Multiphysics Object-Oriented Simulation Environment) and specifically targeting transient modeling of advanced reactors. Griffin has been used recently to model pebble-bed reactors for the Nuclear Regulatory Commission (NRC) Office of Nuclear Regulatory Research and the Advanced Reactor Technology program. This modeling work has focused thus far on the direct calculation of equilibrium cores. This report documents an initial demonstration of a new running-in simulation capability. The new running-in capability is verified using the existing direct equilibrium core calculation capability. A simplified pebble-bed reactor model is then used to demonstrate the running-in simulation capability. This demonstration shows that Griffin is able to simulate years of operation during the running-in phase efficiently with each depletion step taking only several seconds. Two new technologies are also presented in this report which have been developed in Griffin that will be essential for improved accuracy both of the direct equilibrium core computation and the new running-in simulation capability. The first technology is an online cross section generation capability specifically targeted for pebble-bed reactors. This will improve the accuracy of the depletion calculation as the cross sections are generated at the exact core status. This also avoids the difficult step of pre-generating a separate standalone multigroup cross section set. Secondly, a newly implemented discretization for discontinuous finite element method (DFEM) SN transport in cylindrical (RZ) coordinates, which can be solved efficiently using the existing SN sweep solver, is discussed and some results are shown demonstrating the usefulness of the additional accuracy transport provides over a diffusion approximation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Comparison between Alcator C-Mod ICRF experiments and 3D full wave simulations

Reliable modeling of ion cyclotron range of frequencies (ICRF) antenna performance is essential for interpreting present experiments and guiding the design of future reactors. In this work, a 3D model of the Alcator C-Mod field-aligned antenna is implemented in the Petra-M finite-element framework [S. Shiraiwa et al 2023 Nucl. Fusion 63 026024] and benchmarked against experimental results. Four experimental cases are examined. First, the simulated rectified sheath potentials on the antenna limiters are compared with measurements from a power tapering experiment. Second, proof-of-principle far-field sheath simulations are performed. In scenarios with low single-pass absorption, simulations predict enhanced sheath potentials on a distant poloidal limiter in the far field of the antenna, consistent with experimental observations. Third, the simulated antenna loading during edge localized modes follows experimental trends and appears to be dominated by the density gradient at the pedestal. Finally, Petra-M predicts the unintended excitation of high-$k$ ∥ modes during monopole phasing operation, in agreement with experimental evidence of poor wave coupling and accessibility to the plasma core. Overall, reasonable agreement is found between the simulations and experiments. At the same time, areas of imperfect agreement are identified. These provide important guidance on the limits of the current state-of-the-art modeling, which should be kept in mind when using it as a predictive tool for future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Accelerating laser ray tracing in high fidelity physics simulations of laser melting using squeeze U-net

Laser melting is a core component of the ongoing industrial revolution, dubbed Industry 4.0, as lasers facilitate fast and precise melting and fusion in advanced manufacturing. There is a strong need to optimize the laser process using simulations. However, this has proven challenging as high fidelity simulations are needed for predictive modeling and this is currently prohibitively expensive even when run on hundreds of processors on high performance computers. The challenge is capturing complex physics of laser material interaction, fluid dynamics, thermal physics and material phase transformations at various length and time scales. To close this technological gap, we modified a squeeze U-net to accelerate the laser ray tracing component of such high fidelity models by ~4x–40x while preserving the core physics principle of conservation of energy with 97% accuracy. This approach enables the accurate modeling of global laser energy absorption as a function of local surface temperatures and complex surface topologies, which govern the reflection directions and energy losses of laser rays upon interacting with the material surface.

Computer science↗

DOE OSTI2026

Soil microbiome resilience to short-term (30 days, 90 days) and long-term (1000 days) drought

This dataset contains data used for the paper "Drought duration does not impact soil microbiome resilience". The Related References will be updated with a full citation when available. Increasing global droughts exert large but poorly understood effects on the microbial communities and ecology of soil. Microbial communities generally show resilience and return to pre-drought conditions when short-term droughted soils are rewet; soils exposed to long-term drought, however, often show a lag upon rewetting, after which microbial communities may or may not return to their pre-stressed conditions. Though short-term droughts have been widely studied, long-term drought manipulation experiments remain rare, especially those that compare microbial response to short-term and long-term drought in tandem. We conducted a 1000-day drought simulation in controlled laboratory conditions with soil cores collected from a tidal freshwater ecosystem in Washington state, USA, and subsequently exposed them to rewetting for two weeks. We also included short-term (30-day and 90-day) drought and rewet treatments to directly compare microbial community and organic matter responses across drought durations. We found distinct microbial taxa belonging to Firmicutes and Actinobacteria enriched after the 1000-day drought, but not after the short-term droughts. While we hypothesized that the microbial community would recover from a short-term drought after rewetting to resemble pre-drought conditions, our results revealed community dissimilarities between rewet and pre-drought conditions across all drought durations. These findings suggest unique microbial life history strategies within certain microbial phyla that make them successful colonizers during an extended drought period, and the influence of environmental and physiological context on microbial responses to rewetting. The 16SrRNA gene amplicon dataset contains processed DNA sequences in the form of an ASV table with raw unrarefied read counts and representative sequences in .fasta format as described in the ESS-DIVE amplicon sequence reporting format (https://ess-dive.gitbook.io/amplicon-sequencing-reporting-format/instructions). The Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FTICR-MS) dataset consists of processed files containing presence absence data of molecular formulae and molecular characterization of FTICR resolved peaks. The Nuclear Magnetic Resonance (NMR) dataset contains files relevant to NMR spectra and peaks. A sample key file and a sample metadata file is included for the FTICR/NMR and 16S dataset respectively.

1000-day drought↗

DOE OSTI2024

Coastal-Urban-Rural Atmospheric Gradient Experiment (CoURAGE) Science Plan

Understanding the mechanisms governing the urban atmospheric environment is critical for informing urban populations regarding the impacts of climate change and associated mitigation and adaptation measures. Earth system (climate and weather) models have not yet been adapted to provide accurate predictions of climate and weather variability within cities, nor do they provide well-tested representations of the impacts of urban systems on the atmospheric environment. These limitations are largely due to limited field data available for testing and development of these models. We will deploy the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility’s first Mobile Facility (AMF1) to the mid-Atlantic region surrounding the city of Baltimore for the Coast-Urban-Rural Atmospheric Gradient Experiment (CoURAGE). This deployment will create a four-node regional atmospheric observatory network including Baltimore and its three primary surrounding environments – rural, urban, and bay. CoURAGE investigators will study the interactions among the Earth’s surface, the atmospheric boundary layer, aerosols and atmospheric composition, clouds, radiation, and precipitation at each site, and examine how the spatial gradients across the region interact to create the climate conditions in Baltimore. This study will determine the degree to which Baltimore’s atmospheric environment depends on interactive feedbacks in the atmospheric system and the degree to which conditions in Baltimore depend on the surrounding environment. Some topics of interest include how urban land management exacerbates heat waves, the impact of regional mesoscale winds (nocturnal jet, bay breeze) on urban air pollution and cloud cover, and the impact of the urban heat island and aerosol production on heavy precipitation events. Understanding this integrated coast-urban-rural system quantitatively and with good accuracy and precision is critical to informing climate adaptation and mitigation efforts in the city of Baltimore. The understanding gained should be applicable to many similar coastal, mid-latitude urban centers. Another important objective of CoURAGE is to improve the representation of the climate of coastal cities in Earth systems models (ESMs). CoURAGE investigators will use the observations to test current ESMs, identify weaknesses and work towards improved simulations of this complex environment. The ARM core facility will be deployed in the city of Baltimore, complementing the Baltimore Social-Environmental Collaborative (BSEC), a DOE urban integrated field laboratory (UIFL). Ancillary sites will be deployed to rural Maryland northwest of Baltimore, and to the southern end of Kent Island within Chesapeake Bay. The fourth node will be a long-term atmospheric observatory operated in Beltsville, Maryland by Howard University and the Maryland Department of the Environment. Measurements will be conducted for one year, starting in December of 2024. There will be two intensive operational periods (IOPs), one in summer and one in winter, when the ancillary sites will be enhanced with additional balloon launches, tethered balloon system (TBS) operation, and added atmospheric composition measurements.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Project Presentation: Hardware-Based Demonstration of Temperature Control Functions for Reactor Systems

Autonomous thermal regulation in nuclear reactors remains crucial for maintaining stable operation and ensuring the integrity of fuel. To alleviate public skepticism of the safety of nuclear reactors, demonstrating control over this key factor is pivotal. Utilizing electric heat pads to simulate the heat released in a reactor core, thermocouples for temperature monitoring, and an Arduino micro programmable logic controller (PLC) with an embedded proportional-integral-derivative (PID) algorithm for control, a hardware-based demonstration of a reactor heating system will validate the efficacy of reactor control over this key parameter. This report covers internship project presentation as well as relevant experience with nuclear and proposal for project upgrade.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2026

PIE-Enabled Study of Aqueous Corrosion & Zr Hydriding in Cr-Coated Cladding: M3GV-23PN0101132

Chromium-coated zirconium alloy cladding is under investigation as an accident-tolerant fuel (ATF) concept to extend performance to higher burnups via improved oxidation resistance and reduced hydrogen pickup. However, hydrothermal corrosion behavior and hydrogen transport mechanisms governing in-service hydriding of these coatings remain poorly understood. In a collaborative study between Pacific Northwest National Laboratory (PNNL), Idaho National Laboratory (INL), and the University of Huddersfield, cold spray (CS) and physical vapor deposition (PVD) Cr-coated Optimized ZIRLO™ cladding samples were characterized after exposure to PWR-simulated water chemistry under both in-core (neutron irradiation) and out-of-core (aqueous-only) conditions at the MIT Research Reactor. Multi-scale characterization was performed independently at PNNL and INL to evaluate coating integrity, microstructural evolution, Cr/Zr interface chemistry, and hydride formation. Both laboratories observed a consistent divergence in hydriding behavior: out-of-core CS Cr-coated cladding exhibited elevated hydride concentrations exceeding those of uncoated cladding, whereas in-core CS samples showed substantially suppressed hydriding. In-core specimens also displayed irradiation-specific features, including nanoscale voids within the Cr coating and radiation-induced segregation clusters in the Zr substrate. CS coatings retained an interdiffusion-free Cr/Zr bond with no intermetallic layer, whereas PVD coatings exhibited inferior quality with visible cracks and pores. An automated image-based hydride quantification method systematically overestimated bulk hydrogen content relative to inert gas fusion measurements, underscoring the need for standardized sample preparation and reference standards. Interpretation of hydride nucleation and growth is complicated by the open inner-diameter of the specimens, which provides additional hydrogen pathways. The results motivate further foundational studies to support predictive models of hydrothermal corrosion and hydrogen transport in unirradiated and irradiated Cr-coated cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Self-interacting Dark Matter Interpretation of Crater II

Abstract The satellite galaxy Crater II of the Milky Way is extremely cold and exceptionally diffuse. These unusual properties are challenging to understand in the standard model of cold dark matter. We use controlled N -body simulations to investigate the formation of Crater II in self-interacting dark matter (SIDM), where dark matter particles can scatter and thermalize. Taking the orbit motivated by the measurements from Gaia Early Data Release 3, we show a strong self-interacting cross section per particle mass of 60 cm 2 g −1 is favored for Crater II. The simulated SIDM halo, with a 1 kpc core, leads to both a low stellar velocity dispersion and a large half-light radius for Crater II. These characteristics remain robust regardless of the initial stellar distribution.

Astronomy & Astrophysics↗

DOE OSTI2025

Poster: Hardware-Based Demonstration of Temperature Control Functions for Reactor Systems

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗