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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

THE ROLE OF NUCLEAR STAR CLUSTERS IN ENHANCING SUPERMASSIVE BLACK HOLE FEEDING RATES DURING GALAXY MERGERS

During galaxy mergers the gas falls to the center, triggers star formation, and feeds the rapid growth of supermassive black holes (SMBHs). SMBHs respond to this fueling by supplying energy back to the ambient gas. Numerical studies suggest that this feedback is necessary to explain why the properties of SMBHs and the formation of bulges are closely related. This intimate link between the SMBH’s mass and the large scale dynamics and luminosity of the host has proven to be a difficult issue to tackle with simulations due to the inability to resolve all the relevant length scales simultaneously. In this paper we simulate SMBH growth at high-resolution with FLASH, accounting for the gravitational focusing effects of nuclear star clusters (NSCs), which appear to be ubiquitous in galactic nuclei. In the simulations, the NSC core is resolved by a minimum cell size of about 0.001 pc or approximately 10{sup −3} of the cluster’s radius. We discuss the conditions required for effective gas funneling to occur, which are mainly dominated by a relationship between NSC velocity dispersion and the local sound speed, and provide a sub-grid prescription for the augmentation of central SMBH accretion rates in the presence of NSCs. For the conditions expected to persist in the centers of merging galaxies, the resultant large central gas densities in NSCs should produce drastically enhanced embedded SMBH accretion rates—up to an order of magnitude increase can be achieved for gas properties resembling those in large-scale galaxy merger simulations. This will naturally result in faster black hole growth rates and higher luminosities than predicted by the commonly used Bondi–Hoyle–Lyttleton accretion formalism.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Chimera D-Series Gravitational Wave Emission Sourced from Neutrino Anisotropy

Gravitational wave data sourced from the time-dependent anisotropic neutrino emission, as well as the time-dependent fluid quadrupole motion, in the Chimera D-Series three-dimensional core collapse supernova simulations. Data from three models initiated from three different progenitors are presented: D9.6-3D, D15-3D, and D25-3D. Please see the README for more information about the data structure and progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Isotope heat source simulator for testing of space power systems

A reliable isotope heat source simulator was designed for use in a Brayton power system. This simulator is composed of an electrically heated tungsten wire which is wound around a boron nitride core and enclosed in a graphite jacket. Simulator testing was performed at the expected operating temperature of the Brayton power system. Endurance testing for 5012 hours was followed by cycling the simulator temperature. The integrity of this simulator was maintained throughout testing. Alumina beads served as a diffusion barrier to prevent interaction between the tungsten heater and boron nitride core. The simulator was designed to maintain a surface temperature of 1311 to 1366 K (1900 to 2000 F) with a power input of approximately 400 watts. The design concept and the materials used in the simulator make possible man different geometries. This flexibility increases its potential use.

Prok, G. M.↗
DOE OSTI2016

CHANGE OF MAGNETIC FIELD-GAS ALIGNMENT AT THE GRAVITY-DRIVEN ALFVÉNIC TRANSITION IN MOLECULAR CLOUDS: IMPLICATIONS FOR DUST POLARIZATION OBSERVATIONS

Diffuse striations in molecular clouds are preferentially aligned with local magnetic fields, whereas dense filaments tend to be perpendicular to them. When and why this transition occurs remain uncertain. To explore the physics behind this transition, we compute the histogram of relative orientation (HRO) between the density gradient and the magnetic field in three-dimensional magnetohydrodynamic (MHD) simulations of prestellar core formation in shock-compressed regions within giant molecular clouds. We find that, in the magnetically dominated (sub-Alfvénic) post-shock region, the gas structure is preferentially aligned with the local magnetic field. For overdense sub-regions with super-Alfvénic gas, their elongation becomes preferentially perpendicular to the local magnetic field. The transition occurs when self-gravitating gas gains enough kinetic energy from the gravitational acceleration to overcome the magnetic support against the cross-field contraction, which results in a power-law increase of the field strength with density. Similar results can be drawn from HROs in projected two-dimensional maps with integrated column densities and synthetic polarized dust emission. We quantitatively analyze our simulated polarization properties, and interpret the reduced polarization fraction at high column densities as the result of increased distortion of magnetic field directions in trans- or super-Alfvénic gas. Furthermore, we introduce measures of the inclination and tangledness of the magnetic field along the line of sight as the controlling factors of the polarization fraction. Observations of the polarization fraction and angle dispersion can therefore be utilized in studying local magnetic field morphology in star-forming regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Exascale Multiphysics Nuclear Reactor Simulations for Advanced Designs

ENRICO is a coupled application developed under the U.S. Department of Energy's Exascale Computing Project (ECP) targeting the modeling of advanced nuclear reactors. It couples radiation transport with heat and fluid simulation, including the high-fidelity, highresolution Monte-Carlo code Shift and the Computational fluid dynamics code NekRS. NekRS is a highly-performant open-source code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on rapidly convergent high-order spectral element discretizations that feature minimal numerical dissipation and dispersion. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Frontier, Nek5000/RS has recently achieved an unprecedented milestone in breaching over 1 billion spectral elements and 350 billion degrees of freedom. Shift has demonstrated the capability to transport upwards of 1 billion particles per second in full core nuclear reactor simulations featuring complete temperature-dependent, continuous-energy physics on Frontier. Shift achieved a weak-scaling efficiency of 97.8% on 8192 nodes of Frontier and calculated 6 reactions in 214,896 fuel pin regions below 1% statistical error yielding first-of-a-kind resolution for a Monte Carlo transport application.

Hamilton, Steven P.↗
DOE OSTI2026

Chimera F-Series Gravitational Wave Emission Sourced from Matter and Neutrino Anisotropy

This dataset contains gravitational wave data sourced from the the time-dependent fluid quadrupole motion as well as neutrino anisotropy, in the Chimera F-Series two-dimensional core collapse supernova simulations. Data from two models initiated from two different progenitors are presented: F15.78 and F15.79. Please see the README for more information about the data structure and progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

SYMPATHETIC MAGNETIC BREAKOUT CORONAL MASS EJECTIONS FROM PSEUDOSTREAMERS

We present high-resolution 2.5D MHD simulation results of magnetic breakout-initiated coronal mass ejections (CMEs) originating from a coronal pseudostreamer configuration. The coronal null point in the magnetic topology of pseudostreamers means that the initiation of consecutive sympathetic eruptions is a natural consequence of the system's evolution. A generic source region energization process-ideal footpoint shearing parallel to the pseudostreamer arcade polarity inversion lines-is all that is necessary to store sufficient magnetic energy to power consecutive CME eruptions given that the pseudostreamer topology enables the breakout initiation mechanism. The second CME occurs because the eruptive flare reconnection of the first CME simultaneously acts as the overlying pre-eruption breakout reconnection for the sympathetic eruption. We examine the details of the magnetic and kinetic energy evolution and the signatures of the overlying null point distortion, current sheet formation, and magnetic breakout reconnection giving rise to the runaway expansion that drives the flare reconnection below the erupting sheared field core. The numerical simulation's spatial resolution and output cadence are sufficient to resolve the formation of magnetic islands during the reconnection process in both the breakout and eruptive flare current sheets. We quantify the flux transfer between the pseudostreamer arcades and show that the eruptive flare reconnection processes flux {approx}10 times faster than the pre-eruption breakout reconnection. We show that the breakout reconnection jets cause bursty, intermittent upflows along the pseudostreamer stalk, as well as downflows in the adjacent pseudostreamer arcade, both of which may be observable as pre-eruption signatures. Finally, we examine the flux rope CME trajectories and show that the breakout current sheet provides a path of least resistance as an imbalance in the surrounding magnetic energy density and results in a non-radial CME deflection early in the eruption.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A VERY DEEP CHANDRA OBSERVATION OF A1795: THE COLD FRONT AND COOLING WAKE

We present a new analysis of very deep Chandra observations of the galaxy cluster A1795. Utilizing nearly 750 ks of net ACIS imaging, we are able to resolve the thermodynamic structure of the intracluster medium (ICM) on length scales of ∼1 kpc near the cool core. We find several previously unresolved structures, including a high pressure feature to the north of the Brightest Cluster Galaxy (BCG) that appears to arise from the bulk motion of A1795's cool core. To the south of the cool core, we find low temperature (∼3 keV), diffuse ICM gas extending for distances of ∼50 kpc spatially coincident with previously identified filaments of Hα emission. Gas at similar temperatures is also detected in adjacent regions without any Hα emission. The X-ray gas coincident with the Hα filament has been measured to be cooling spectroscopically at a rate of ∼1 M{sub ⊙} yr{sup −1}, consistent with measurements of the star formation rate in this region as inferred from ultraviolet (UV) observations, suggesting that the star formation in this filament as inferred by its Hα and UV emission can trace its origin to the rapid cooling of dense, X-ray emitting gas. The Hα filament is not a unique site of cooler ICM, however, as ICM at similar temperatures and even higher metallicities not cospatial with Hα emission is observed just to the west of the Hα filament, suggesting that it may have been uplifted by A1795's central active galaxy. Further simulations of cool core sloshing and active galactic nucleus feedback operating in concert with one another will be necessary to understand how such a dynamic cool core region may have originated and why the Hα emission is so localized with respect to the cool X-ray gas.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2014

Asymmetries in Core-Collapse Supernovae from Maps of Radioactiver 44Ti in Cassiopeia A

Asymmetry is required by most numerical simulations of stellar core-collapse explosions, but the form it takes differs significantly among models. The spatial distribution of radioactive 44Ti, synthesized in an exploding star near the boundary between material falling back onto the collapsing core and that ejected into the surroundingmedium1, directly probes the explosion asymmetries. Cassiopeia A is a young2, nearby3, core-collapse4 remnant from which 44Ti emission has previously been detected5-8 but not imaged. Asymmetries in the explosion have been indirectly inferred from a high ratio of observed 44Ti emission to estimated 56Ni emission9, from optical light echoes10, and from jet-like features seen in the X-ray11 and optical12 ejecta. Here we report spatial maps and spectral properties of the 44Ti in Cassiopeia A. This may explain the unexpected lack of correlation between the 44Ti and iron X-ray emission, the latter being visible only in shock-heated material. The observed spatial distribution rules out symmetric explosions even with a high level of convective mixing, as well as highly asymmetric bipolar explosions resulting from a fast-rotating progenitor. Instead, these observations provide strong evidence for the development of low-mode convective instabilities in core-collapse supernovae.

stellarcore-collapse explosions↗
DOE OSTI2017

Modeling and Preparation for Experimental Testing of Heat Fluxes on W7-X Divertor Scraper Elements

Simulations of heat fluxes to the plasma facing components in the Wendelstein 7-X stellarator will be tested in its next operational phase. The simulations consist of core transport calculations that determine the evolution of the kinetic profiles and the toroidal current, which modifies the fluxes to the divertor, as the magnetic geometry changes. An additional divertor component, the scraper element, was designed to protect the edges of the primary divertor throughout this evolution during certain high-power long-pulse operational scenarios. The effect of unknown parameters of the heat flux calculations, namely, the cross-field thermal diffusivity and the magnetic field structure, is explored. The predicted scaling of the heat flux widths and magnitudes is presented, along with a new method of calculating the 3-D magnetic field structure required to perform the flux calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Theory for Neutron Star and Black Hole Kicks and Induced Spins

Abstract Using 20 long-term 3D core-collapse supernova simulations, we find that lower compactness progenitors that explode quasi-spherically due to the short delay to explosion experience smaller neutron star recoil kicks in the ∼100−200 km s −1 range, while higher compactness progenitors that explode later and more aspherically leave neutron stars with kicks in the ∼300−1000 km s −1 range. In addition, we find that these two classes are correlated with the gravitational mass of the neutron star. This correlation suggests that the survival of binary neutron star systems may in part be due to their lower kick speeds. We also find a correlation between the kick and both the mass dipole of the ejecta and the explosion energy. Furthermore, one channel of black hole birth leaves masses of ∼10 M ⊙ , is not accompanied by a neutrino-driven explosion, and experiences small kicks. A second channel is through a vigorous explosion that leaves behind a black hole with a mass of ∼3.0 M ⊙ kicked to high speeds. We find that the induced spins of nascent neutron stars range from seconds to ∼10 ms, but do not yet see a significant spin/kick correlation for pulsars. We suggest that if an initial spin biases the explosion direction, a spin/kick correlation would be a common byproduct of the neutrino mechanism of core-collapse supernovae. Finally, the induced spin in explosive black hole formation is likely large and in the collapsar range. This new 3D model suite provides a greatly expanded perspective and appears to explain some observed pulsar properties by default.

Astronomy & Astrophysics↗
DOE OSTI2019

Hydrodynamic simulations of disrupted planetary accretion discs inside the core of an AGB star

ABSTRACT Volume complete sky surveys provide evidence for a binary origin for the formation of isolated white dwarfs with magnetic fields in excess of a MegaGauss. Interestingly, not a single high-field magnetic white dwarf has been found in a detached system, suggesting that if the progenitors are indeed binaries, the companion must be removed or merge during formation. An origin scenario consistent with observations involves the engulfment, inspiral, and subsequent tidal disruption of a low-mass companion in the interior of a giant star during a common envelope phase. Material from the shredded companion forms a cold accretion disc embedded in the hot ambient around the proto-white dwarf. Entrainment of hot material may evaporate the disc before it can sufficiently amplify the magnetic field, which typically requires at least a few orbits of the disc. Using three-dimensional hydrodynamic simulations of accretion discs with masses between 1 and 10 times the mass of Jupiter inside the core of an Asymptotic Giant Branch star, we find that the discs survive for at least 10 orbits (and likely for 100 orbits), sufficient for strong magnetic fields to develop.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Graph-Based Representations and Applications to Process Simulation

Rapid and robust convergence of a process flowsheet is critical to enable large-scale simulations that address core scientific questions related to process design, optimization, and sustainability. However, due to the highly coupled and nonlinear nature of chemical processes, efficiently solving a flowsheet remains a challenge. In this work, we show that graph representations of the underlying physical phenomena in unit operations may help identify potential avenues to systematically reformulate the network of equations and enable more robust topology-based convergence of flowsheets. To this end, we developed graph abstractions of the governing equations of vapor-liquid and liquid-liquid equilibrium separation equipment. These graph abstractions consist of a mesh of interconnected variable nodes and equation nodes that are systematically generated through PhenomeNode, a new open-source library in Python developed in this study. We show that partitioning the graph into separate mass, energy, and equilibrium subgraphs can help decouple nonlinearities and guide decomposition algorithms. By employing the graph abstraction on an industrial separation process for separating glacial acetic acid from water, we implemented a new block decomposition scheme in BioSTEAM and demonstrated that this can accelerate convergence over a traditional sequential modular approach.

Distillation↗
NASA NTRS2022

New Constraints on the Dark Matter Density Profiles of Dwarf Galaxies from Proper Motions of Globular Cluster Streams

The central density profiles in dwarf galaxy halos depend strongly on the nature of dark matter (DM). Recently, in Malhan et al. we employed N-body simulations to show that the cuspy cold DM subhalos predicted by cosmological simulations can be differentiated from cored subhalos using the properties of accreted globular cluster (GC) streams since these GCs experience tidal stripping within their parent halos prior to accretion onto the Milky Way. We previously found that clusters that are accreted within cuspy subhalos produce streams with larger physical widths and higher dispersions in line-of-sight velocity and angular momentum than streams that are accreted within cored subhalos. Here, we use the same suite of simulations to demonstrate that the dispersion in the tangential velocities of streams (σ 𝜈 Tan) is also sensitive to the central DM density profiles of their parent dwarfs and GCs that they were accreted from; cuspy subhalos produce streams with larger σ 𝜈 Tan than those accreted inside cored subhalos. Using Gaia EDR3 observations of multiple GC streams we compare their (σ 𝜈 Tan) values with simulations. The measured (σ 𝜈 Tan) values are consistent with both an “in situ” origin and with accretion inside cored subhalos of M ∼ 10(8–9) Me (or very low-mass cuspy subhalos of mass ∼10(8) Mꙩ ). Despite the large current uncertainties in σ 𝜈 Tan, we find a low probability that any of the progenitor GCs were accreted from cuspy subhalos of M ≳ 10 (9) Mꙩ. The uncertainties on Gaia tangential velocity measurements are expected to decrease in future and will allow for stronger constraints on subhalo DM density profiles.

Khyati Malhan↗
DOE OSTI2022

Deterministic modelling of the SPERT IV reactor transients using the multi-physics capability of WIMS

The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

New Constraints on the Dark Matter Density Profiles of Dwarf Galaxies from Proper Motions of Globular Cluster Streams

The central density profiles in dwarf galaxy halos depend strongly on the nature of dark matter (DM). Recently, in Malhan et al. we employed N-body simulations to show that the cuspy cold DM subhalos predicted by cosmological simulations can be differentiated from cored subhalos using the properties of accreted globular cluster (GC) streams since these GCs experience tidal stripping within their parent halos prior to accretion onto the Milky Way. We previously found that clusters that are accreted within cuspy subhalos produce streams with larger physical widths and higher dispersions in line-of-sight velocity and angular momentum than streams that are accreted within cored subhalos. Here, we use the same suite of simulations to demonstrate that the dispersion in the tangential velocities of streams (${\sigma }_{{v}_{\mathrm{Tan}}}$) is also sensitive to the central DM density profiles of their parent dwarfs and GCs that they were accreted from; cuspy subhalos produce streams with larger ${\sigma }_{{v}_{\mathrm{Tan}}}$ than those accreted inside cored subhalos. Using Gaia EDR3 observations of multiple GC streams we compare their ${\sigma }_{{v}_{\mathrm{Tan}}}$ values with simulations. The measured ${\sigma }_{{v}_{\mathrm{Tan}}}$ values are consistent with both an “in situ” origin and with accretion inside cored subhalos of M ~ 10 8–9 M ⊙ (or very low-mass cuspy subhalos of mass ~10 8 M ⊙ ). Despite the large current uncertainties in ${\sigma }_{{v}_{\mathrm{Tan}}}$, we find a low probability that any of the progenitor GCs were accreted from cuspy subhalos of M ≳ 10 9 M ⊙ . The uncertainties on Gaia tangential velocity measurements are expected to decrease in future and will allow for stronger constraints on subhalo DM density profiles.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2018

Electrically Heated Testing of the Kilowatt Reactor Using Stirling TechnologY (KRUSTY) Experiment Using a Depleted Uranium Core

The Kilopower project aims to develop and demonstrate scalable fission-based power technology for systems capable of delivering 1 to 10 kW of electric power with a specific power ranging from 2.5 to 6.5 W/kg. This technology could enable high-power science missions or could be used to provide surface power for manned missions to the Moon or Mars. NASA has partnered with the U.S. Department of Energy's National Nuclear Security Administration, Los Alamos National Laboratory, Nevada National Security Site (NNSS), and Y−12 National Security Complex to develop and test a prototypic reactor and power system using existing facilities and infrastructure. This technology demonstration, referred to as the "Kilowatt Reactor Using Stirling TechnologY (KRUSTY)," will undergo nuclear ground testing by the end of calendar year (CY) 2017 at the NNSS. The 1-kWe variation of the Kilopower system was chosen for the KRUSTY demonstration. The concept for the 1-kWe flight system consists of a 4 kWt highly enriched uranium-molybdenum reactor operating at 800 degC coupled to sodium heat pipes. The heat pipes deliver heat to the hot ends of eight 125-W Stirling convertors producing a net electrical output of 1 kW. Waste heat is rejected using titanium-water heat pipes coupled to carbon composite radiator panels. The KRUSTY test, based on this design, uses a prototypic highly enriched uranium-molybdenum core coupled to prototypic sodium heat pipes. The heat pipes transfer heat to two Advanced Stirling Convertors (ASC−E2s) and six thermal simulators, which simulate the thermal draw of full-scale power conversion units. Thermal simulators and Stirling engines are gas cooled. The most recent project milestone was the completion of nonnuclear system-level testing using an electrically heated depleted uranium (DU) (nonfissioning) reactor core simulator at the NASA Glenn Research Center. System-level testing has validated performance predictions and has demonstrated system-level operation and control in a test configuration that replicates the one to be used at the Device Assembly Facility (DAF) at the NNSS. Fabrication, assembly, and testing of the DU core has allowed for higher fidelity system-level testing at Glenn, and has validated the fabrication methods to be used on the highly enriched uranium core that will supply heat for the DAF KRUSTY demonstration.

Briggs, Maxwell H.↗
DOE OSTI2026

Alternating salt and freshwater floods of coastal soils impact soil structure, hydraulic properties, and oxygen dynamics

Coastal soils are increasingly impacted by hydrologic intensification in the form of rising sea level and flooding from storm surges and precipitation. Alternating saltwater (SW) and freshwater (FW) exposure has the potential to disperse colloids, which can lead to disintegration of soil structure, clogging of pore spaces, and reduction in ecologically and biogeochemically important functions like infiltration or gas exchange. To investigate how hydrologic intensification affects soil structure and oxygen dynamics, we conducted a series of laboratory-based flood simulations. Intact soil cores of the A and B horizons (22 total) from the toe slope of an upland coastal forest along the western shore of Chesapeake Bay. We subjected the cores to 24 h of saturation with either FW or alternating brackish SW and FW, with a 24-h draining event in between flood events. Oxygen diffusion into soil during draining was reduced by up to 30% for soils flooded alternating SW-FW compared to soils flooded only with FW. We attributed the reduced oxygen diffusion to clogging of smaller pores by colloids, with colloid redistribution observed as an increase in specific surface area down the core profile of up to 59%. After three SW-FW floods (six floods total), there were significant changes in pore size distribution, significant redistribution of colloids, and the A horizon became sodic. We concluded that a small number of SW flooding events can induce a measurable change in soil physical properties that directly impacts the biogeochemical dynamics.

Rod, Kenton A. [Pacific Northwest National Laborat↗