Search NASA⌕ Search

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.

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

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS2004

High Speed Networking and Large-scale Simulation in Geodynamics

Large-scale numerical simulation has been one of the most important approaches for understanding global geodynamical processes. In this approach, peta-scale floating point operations (pflops) are often required to carry out a single physically-meaningful numerical experiment. For example, to model convective flow in the Earth's core and generation of the geomagnetic field (geodynamo), simulation for one magnetic free-decay time (approximately 15000 years) with a modest resolution of 150 in three spatial dimensions would require approximately 0.2 pflops. If such a numerical model is used to predict geomagnetic secular variation over decades and longer, with e.g. an ensemble Kalman filter assimilation approach, approximately 30 (and perhaps more) independent simulations of similar scales would be needed for one data assimilation analysis. Obviously, such a simulation would require an enormous computing resource that exceeds the capacity of a single facility currently available at our disposal. One solution is to utilize a very fast network (e.g. 10Gb optical networks) and available middleware (e.g. Globus Toolkit) to allocate available but often heterogeneous resources for such large-scale computing efforts. At NASA GSFC, we are experimenting with such an approach by networking several clusters for geomagnetic data assimilation research. We shall present our initial testing results in the meeting.

Kuang, Weijia↗
NASA NTRS1979

High temperature UF6 RF plasma experiments applicable to uranium plasma core reactors

An investigation was conducted using a 1.2 MW RF induction heater facility to aid in developing the technology necessary for designing a self critical fissioning uranium plasma core reactor. Pure, high temperature uranium hexafluoride (UF6) was injected into an argon fluid mechanically confined, steady state, RF heated plasma while employing different exhaust systems and diagnostic techniques to simulate and investigate some potential characteristics of uranium plasma core nuclear reactors. The development of techniques and equipment for fluid mechanical confinement of RF heated uranium plasmas with a high density of uranium vapor within the plasma, while simultaneously minimizing deposition of uranium and uranium compounds on the test chamber peripheral wall, endwall surfaces, and primary exhaust ducts, is discussed. The material tests and handling techniques suitable for use with high temperature, high pressure, gaseous UF6 are described and the development of complementary diagnostic instrumentation and measurement techniques to characterize the uranium plasma, effluent exhaust gases, and residue deposited on the test chamber and exhaust system components is reported.

Roman, W. C.↗
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↗
NASA NTRS1996

Research In Nonlinear Flight Control for Tiltrotor Aircraft Operating in the Terminal Area

The research during the first year of the effort focused on the implementation of the recently developed combination of neural net work adaptive control and feedback linearization. At the core of this research is the comprehensive simulation code Generic Tiltrotor Simulator (GTRS) of the XV-15 tilt rotor aircraft. For this research the GTRS code has been ported to a Fortran environment for use on PC. The emphasis of the research is on terminal area approach procedures, including conversion from aircraft to helicopter configuration. This report focuses on the longitudinal control which is the more challenging case for augmentation. Therefore, an attitude command attitude hold (ACAH) control augmentation is considered which is typically used for the pitch channel during approach procedures. To evaluate the performance of the neural network adaptive control architecture it was necessary to develop a set of low order pilot models capable of performing such tasks as, follow desired altitude profiles, follow desired speed profiles, operate on both sides of powercurve, convert, including flaps as well as mastangle changes, operate with different stability and control augmentation system (SCAS) modes. The pilot models are divided in two sets, one for the backside of the powercurve and one for the frontside. These two sets are linearly blended with speed. The mastangle is also scheduled with speed. Different aspects of the proposed architecture for the neural network (NNW) augmented model inversion were also demonstrated. The demonstration involved implementation of a NNW architecture using linearized models from GTRS, including rotor states, to represent the XV-15 at various operating points. The dynamics used for the model inversion were based on the XV-15 operating at 30 Kts, with residualized rotor dynamics, and not including cross coupling between translational and rotational states. The neural network demonstrated ACAH control under various circumstances. Future efforts will include the implementation into the Fortran environment of GTRS, including pilot modeling and NNW augmentation for the lateral channels. These efforts should lead to the development of architectures that will provide for fully automated approach, using similar strategies.

Calise, A. J.↗