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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Towards performance portability in the Spark astrophysical magnetohydrodynamics solver in the Flash-X simulation framework

Simulations of core-collapse supernovae, and other astrophysical phenomena, are quintessential extreme-scale computing challenges. For core-collapse supernova simulations to be carried out by the ExaStar project under the Exascale Computing Project umbrella, a robust, efficient, and state-of-the-art magnetohydrodynamics solver is a critical requirement. In Flash-X, the primary software instrument for ExaStar, a new magnetohydrodynamics solver has been designed and implemented from the ground up to achieve accuracy and efficiency for simulations of complex astrophysical flows. This new solver, dubbed Spark, uses high-order spatial reconstruction, Runge-Kutta time integration, and an efficient cell-centered approach to satisfying the divergence-free condition for the magnetic fields. Spark was written to be optimized for data locality in cache hierarchy of CPUs. Since data locality optimizations for cache hierarchy are not directly compatible with those of accelerators, we have taken the approach of using program synthesis to avoid massive amounts of code replication that would be necessary if we were to maintain two different versions of the solver. Our program synthesis relies on a simple key-dictionary approach, implemented in python, that enables us to assemble the version of the solver suitable for the target hardware from code fragments identified by specific keys. In this work, we describe the data locality optimizations of the solver for CPUs and accelerators and the program synthesis tools that enable this portability. We also detail the parallel performance of Spark for both CPUs and accelerators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Electromagnetic total- f algorithm for gyrokinetic particle-in-cell simulations of boundary plasma in XGC

We report the simplified δf mixed-variable/pullback electromagnetic simulation algorithm implemented in XGC for core plasma simulations by Cole et al. [Phys. Plasmas 28, 034501 (2021)] has been generalized to a total- f electromagnetic algorithm that can include, for the first time, the boundary plasma in diverted magnetic geometry with neutral particle recycling, turbulence, and neoclassical physics. The δf mixed-variable/pullback electromagnetic implementation is based on the pioneering work by Kleiber and Mischenko et al. [Kleiber et al., Phys. Plasmas 23, 032501 (2016); Mishchenko et al., Comput. Phys. Commun. 238, 194 (2019)]. An electromagnetic demonstration simulation is performed in a DIII-D-like, H-mode boundary plasma, including a corresponding comparative electrostatic simulation, which confirms that the electromagnetic simulation is necessary for a higher fidelity understanding of the electron particle and heat transport even at the low- β pedestal foot in the vicinity of the magnetic separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Benchmark of Atucha-2 PHWR RELAP5-3D control rod model by Monte Carlo MCNP5 core calculation

Atucha-2 is a Siemens-designed PHWR reactor under construction in the Republic of Argentina. Its geometrical complexity and peculiarities require the adoption of advanced Monte Carlo codes for performing realistic neutronic simulations. Therefore core models of Atucha-2 PHWR were developed using MCNP5. In this work a methodology was set up to collect the flux in the hexagonal mesh by which the Atucha-2 core is represented. The scope of this activity is to evaluate the effect of obliquely inserted control rod on neutron flux in order to validate the RELAP5-3D{sup C}/NESTLE three dimensional neutron kinetic coupled thermal-hydraulic model, applied by GRNSPG/UNIPI for performing selected transients of Chapter 15 FSAR of Atucha-2. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2010

Automatic Tip Vortex Core Profiling for Numerical Flow Simulations of Rotorcraft in Hover

An automated approach is presented that extracts visual and quantitative data from vortex cores produced by Navier-Stokes simulations of rotorcraft in hover mode. This approach extracts contiguous rotor tip vortex-core trajectories, cross-flow velocity profiles, and vortex-core diameter variation with wake age (azimuth angle). This automated approach is faster and more accurate than a conventional manual approach. Moreover, this new approach allows for an efficient way to quantitatively compare vortex-core profiles from different flow simulations, e.g., grid resolution studies, and validate computed results with experimental data

Kao, David L.↗
DOE OSTI2025

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation↗
DOE OSTI2013

MULTI-DIMENSIONAL MODELS FOR DOUBLE DETONATION IN SUB-CHANDRASEKHAR MASS WHITE DWARFS

Using two-dimensional and three-dimensional simulations, we study the ''robustness'' of the double detonation scenario for Type Ia supernovae, in which a detonation in the helium shell of a carbon-oxygen white dwarf induces a secondary detonation in the underlying core. We find that a helium detonation cannot easily descend into the core unless it commences (artificially) well above the hottest layer calculated for the helium shell in current presupernova models. Compressional waves induced by the sliding helium detonation, however, robustly generate hot spots which trigger a detonation in the core. Our simulations show that this is true even for non-axisymmetric initial conditions. If the helium is ignited at multiple points, then the internal waves can pass through one another or be reflected, but this added complexity does not defeat the generation of the hot spot. The ignition of very low-mass helium shells depends on whether a thermonuclear runaway can simultaneously commence in a sufficiently large region.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

Numerical simulations of thermal-chemical instabilities at the core-mantle boundary

Numerical simulations of thermal-chemical instabilities in the D-double-prime layer at the core-mantle boundary are presented which show that strong lateral heterogeneities in the composition and density fields can be initiated and maintained dynamically if there is continuous replenishment of material from subduced slabs coming from the upper mantle. These chemical instabilities have a tendency to migrate laterally and may help to support core-mantle boundary topography with short and long wavelengths. The thermal-chemical flows produce a relatively stagnant D-double-prime layer with strong lateral and temporal variations in basal heat flux, which gives rise to thermal core-mantle interactions influencing the geodynamo.

Hansen, Ulrich↗
DOE OSTI2014

Three-dimensional NDE of VHTR core components via simulation-based testing. Final report

A next generation, simulation-driven-and-enabled testing platform is developed for the 3D detection and characterization of defects and damage in nuclear graphite and composite structures in Very High Temperature Reactors (VHTRs). The proposed work addresses the critical need for the development of high-fidelity Non-Destructive Examination (NDE) technologies for as-manufactured and replaceable in-service VHTR components. Centered around the novel use of elastic (sonic and ultrasonic) waves, this project deploys a robust, non-iterative inverse solution for the 3D defect reconstruction together with a non-contact, laser-based approach to the measurement of experimental waveforms in VHTR core components. In particular, this research (1) deploys three-dimensional Scanning Laser Doppler Vibrometry (3D SLDV) as a means to accurately and remotely measure 3D displacement waveforms over the accessible surface of a VHTR core component excited by mechanical vibratory source; (2) implements a powerful new inverse technique, based on the concept of Topological Sensitivity (TS), for non-iterative elastic waveform tomography of internal defects - that permits robust 3D detection, reconstruction and characterization of discrete damage (e.g. holes and fractures) in nuclear graphite from limited-aperture NDE measurements; (3) implements state-of-the art computational (finite element) model that caters for accurately simulating elastic wave propagation in 3D blocks of nuclear graphite; (4) integrates the SLDV testing methodology with the TS imaging algorithm into a non-contact, high-fidelity NDE platform for the 3D reconstruction and characterization of defects and damage in VHTR core components; and (5) applies the proposed methodology to VHTR core component samples (both two- and three-dimensional) with a priori induced, discrete damage in the form of holes and fractures. Overall, the newly established SLDV-TS testing platform represents a next-generation NDE tool that surpasses all existing techniques for the 3D ultrasonic imaging of material damage from non-contact, limited-aperture waveform measurements. Outlook. The next stage in the development of this technology includes items such as (a) non-contact generation of mechanical vibrations in VHTR components via thermal expansion created by high-intensity laser; (b) development and incorporation of Synthetic Aperture Focusing Technique (SAFT) for elevating the accuracy of 3D imaging in highly noisy environments with minimal accessible surface; (c) further analytical and computational developments to facilitate the reconstruction of diffuse damage (e.g. microcracks) in nuclear graphite as they lead to the dispersion of elastic waves, (d) concept of model updating for accurate tracking of the evolution of material damage via periodic inspections; (d) adoption of the Bayesian framework to obtain information on the certainty of obtained images; and (e) optimization of the computational scheme toward real-time, model-based imaging of damage in VHTR core components.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Simulating Valence-to-Core X-ray Emission Spectroscopy of Transition Metal Complexes with Time-Dependent Density Functional Theory

Valence-to-core (VtC) X-ray emission spectroscopy (XES) has emerged as a power- ful technique for the structural characterization of complex organometallic compounds in realistic environments. Since the spectrum represents electronic transitions from the ligand molecular orbitals to the core holes of the metal centers, the approach is more chemically sensitive to the metal-ligand bonding character compared with con- ventional X-ray absorption techniques. In this paper we study how linear-response time-dependent density functional theory (LR-TDDFT) can be harnessed to simulate K-edge VtC X-ray emission spectra reliably. LR-TDDFT allows one to go beyond the single-particle picture that has been extensively used to simulate VtC-XES. We con- sider seven low- and high-spin model complexes involving chromium, manganese and iron transition metal centers. Our results are in good agreement with experiment.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2014

Synthetic observations of the evolution of starless cores in a molecular cloud simulation: Comparisons with JCMT data and predictions for alma

Interpreting the nature of starless cores has been a prominent goal in star formation for many years. In order to characterize the evolutionary stages of these objects, we perform synthetic observations of a numerical simulation of a turbulent molecular cloud. We find that nearly all cores that we detect are associated with filaments and eventually form protostars. We conclude that observed starless cores that appear Jeans unstable are only marginally larger than their respective Jeans masses (within a factor of three). We note single dish observations such as those performed with the James Clerk Maxwell Telescope appear to miss significant core structure on small scales because of beam averaging. Finally, we predict that interferometric observations with Atacama Large Millimeter Array Cycle 1 will resolve the important small-scale structure, which has so far been missed by millimeter wavelength observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

LOW MACH NUMBER MODELING OF CORE CONVECTION IN MASSIVE STARS

This work presents three-dimensional simulations of core convection in a 15 M{sub Sun} star halfway through its main sequence lifetime. To perform the necessary long-time calculations, we use the low Mach number code MAESTRO, with initial conditions taken from a one-dimensional stellar model. We first identify several key factors that the one-dimensional initial model must satisfy to ensure efficient simulation of the convection process. We then use the three-dimensional simulations to examine the effects of two common modeling choices on the resulting convective flow: using a fixed composition approximation and using a reduced domain size. We find that using a fixed composition model actually increases the computational cost relative to using the full multi-species model because the fixed composition system takes longer to reach convection that is in a quasi-static state. Using a reduced (octant rather than full sphere) simulation domain yields flow with statistical properties that are within a factor of two of the full sphere simulation values. Both the octant and full sphere simulations show similar mixing across the convection zone boundary that is consistent with the turbulent entrainment model. However, the global character of the flow is distinctly different in the octant simulation, showing more rapid changes in the large-scale structure of the flow and thus a more isotropic flow on average.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Enhanced Performance Fast Reactors with Engineered Passive Safety Systems

This report summarizes the work and results of NEUP Project #15-8251, "Enhanced Performance Fast Reactors with Engineered Passive Safety System," performed from 2015-2018 at UC Berkeley and Argonne National Laboratory. This project concerns the incorporation of a novel engineered passive safety system, called the Autonomous Reactivity Control (ARC) system, into Sodium-cooled Fast Reactors (SFRs) to improve their transient and economic performance. The first part of this project consisted of fleshing out the basic principles and the manufacturing process for an SFR fuel assembly with the ARC system installed, with the details provided in high-quality, step-by-step CAD drawings. Once the detailed design is laid out, high-fidelity simulations were performed to optimize the dimensions of important ARC system components and to quantify the standalone time-response of ARC systems in a typical SFR environment. These high-fidelity simulations were used to create reduced-order models of ARC system response which could be incorporated into the SAS4A/SASSYS-1 whole-core transient models. In order to accurately account for the time-response of ARC systems in whole-core transient simulations, modifications were made to the SAS4A/SASSYS-1 transient code. These new capabilities were tested and incorporated into the production version of the code, and have been used extensively in this analysis. Before performing whole-core transient calculations, the two reference cores selected for this study to cover a wide design space were simulated and characterized. The two reference cores, (1) a medium-sized oxide-fueled Advanced Burner Reactor and (2) a large metal-fueled Breed-and-Burn core, challenge conventional means of achieving passive safety. Because the two cores vary substantially, simulating the ARC system in both allows for conclusions to be drawn about the applicability of ARC systems in different SFRs. Once both cores were fully characterized, transient calculations were performed to understand the complicated interplay between the ARC system and other inherent reactivity feedback mechanisms. By performing parametric studies of the transient behavior, ARC system designs that best improve transient performance in each core were identified. Once transient improvements were established, the transient calculations were extended to improve reactor economics through design modifications. Upon establishing the potential benefits of ARC system inclusion, the study examines the negative neutronic implications of incorporating ARC systems into the reference cores. The study concludes by discussing a number of knowledge gaps that exist which may have important implications for ARC system implementation. Ultimately, the results of this project indicate that ARC systems may provide an effective means to improve transient performance in certain types of systems. The transient simulations in the oxide Advanced Burner Reactor showed that significant margin to coolant boiling and fuel melting could be gained through the use of a properly designed ARC system, thus improving inherent safety. While oscillatory behavior was induced in the transient response for certain combinations of ARC system parameters, the oscillation-free design space that can be effectively used to improve performance was found to be large, and the additional margin was demonstrated to allow economical improvements to the core design. Conversely, difficulties were encountered in the Breed-and-Burn core due to fundamental aspects of the core design. The combination of a particularly strong positive coolant void feedback and the use of metallic fuel allowed for diverging oscillations to manifest in response to certain transient initiators, and in the end it was found that the current ARC system design could not simultaneously enable passive safety in all transients examined. In order for complete passive safety to be achieved, a mechanism for slowing the reverse actuation of ARC systems is necessary. Regardless, the standard ARC system was able to improve the performance of two out of the three examined transients by increasing the margin to boiling in Loss of Heat Sink transients and avoiding boiling in Transient Overpower accidents. Additionally, no significant barriers to implementation were found, although the issue of tritium production should be closely monitored in any core incorporating ARC systems. Still, a number of technology gaps exist, and these are discussed at the conclusion of the report to recommend directions for future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

The Outer Rim Simulation: A Path to Many-core Supercomputers

Herein we describe the Outer Rim cosmological simulation, one of the largest high-resolution N-body simulations performed to date, aimed at promoting science to be carried out with large-scale structure surveys. The simulation covers a volume of (4.225 Gpc) 3 and evolves more than one trillion particles. It was executed on Mira, a BlueGene/Q system at the Argonne Leadership Computing Facility. We discuss some of the computational challenges posed by a system such as Mira, a many-core supercomputer, and how the simulation code, Hardware/Hybrid Accelerated Cosmology Code, has been designed to overcome these challenges. We have carried out a large range of analyses on the simulation data and we report on the results as well as the data products that have been generated. The full data set generated by the simulation totals more than 5 PB of data, making curation and handling of the data a large challenge in itself. The simulation results have been used to generate synthetic catalogs for large-scale structure surveys, including Dark Energy Spectroscopic Instrument and eBOSS, as well as experiments on the cosmic microwave background. A detailed catalog for the data challenges of the Large Synoptic Survey Telescope Dark Energy Science Collaboration has been created as well. We publicly release some of the Outer Rim halo catalogs, downsampled particle information, and lightcone data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Coupled neutronic-thermal-hydraulic simulations of the European SFR core

Within the European SFR (Sodium-cooled Fast Reactor) - Safety Measures Assessment and Research Tools (ESFR-SMART) project, steady-state coupled simulation of the ESFR core has been performed using several core analysis packages, with the objective of quantifying the coupling effect. Focus is on the fuel Doppler effect and coolant expansion effect. Standalone neutronics calculations in TRACE/PARCS (PSI), DYN3D (HZDR) and WIMS (Jacobs) showed superb agreement with the reference Serpent power distribution (root mean square - rms discrepancies of 1.3%, 1.5% and 0.7% respectively). Results for COUNTHER (CIEMAT) were also in reasonable agreement, but with a somewhat higher discrepancy of 3.7%. Temperature distributions from thermal-hydraulic calculations were also compared and are found to be in good agreement. The effect of Doppler and coolant density feedback on core power distribution was predicted by TRACE/PARCS, DYN3D and WIMS to be between 0.4% and 0.8% rms difference in assembly powers. Reactivity coefficients for perturbations in the inlet temperature, flow rate and core power were shown to be negative for these three codes, with values of roughly -0.5 pcm/C. degree, -0.3 pcm/C. degree and -3.5 pcm/C. degree respectively. A preliminary investigation of differential thermal expansion effects indicates that this may have a significant effect on core power distribution of a few %, greater than anticipated a priori and may warrant inclusion in coupled core analysis to ensure the accurate calculation of power distributions. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Condensed phase QM/MM simulations utilizing the exchange core functions to describe exchange repulsions at the QM boundary region

In a recent work, we developed a method [H. Takahashi et al., J. Chem. Phys. 143, 084104 (2015)] referred to as exchange-core function (ECF) approach, to compute exchange repulsion E{sub ex} between solute and solvent in the framework of the quantum mechanical (QM)/molecular mechanical (MM) method. The ECF, represented with a Slater function, plays an essential role in determining E{sub ex} on the basis of the overlap model. In the work of Takahashi et al. [J. Chem. Phys. 143, 084104 (2015)], it was demonstrated that our approach is successful in computing the hydrogen bond energies of minimal QM/MM systems including a cationic QM solute. We provide in this paper the extension of the ECF approach to the free energy calculation in condensed phase QM/MM systems by combining the ECF and the QM/MM-ER approach [H. Takahashi et al., J. Chem. Phys. 121, 3989 (2004)]. By virtue of the theory of solutions in energy representation, the free energy contribution δμ{sub ex} from the exchange repulsion was naturally formulated. We found that the ECF approach in combination with QM/MM-ER gives a substantial improvement on the calculation of the hydration free energy of a hydronium ion. This can be attributed to the fact that the ECF reasonably realizes the contraction of the electron density of the cation due to the deficit of an electron.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2018

High-Resolution Modeling of the Dust and Water Cycles with the NASA Ames Mars Global Climate Model

NASA’s Mars Climate Modeling Center at Ames Research Center is currently undergoing an exciting period of growth in personnel, modeling capabilities, and science productivity. We are transitioning from our legacy Arakawa C-grid finite-difference dynamical core to the NOAA/GFDL cubed-sphere finite-volume dynamical core for simulating the climate of Mars in a global framework. This highly parallelized core is scalable and flexible, which allows for significant improvements in the horizontal and vertical resolutions of our simulations. We have implemented the Ames water ice cloud microphysics package described in Haberle et al. (2018) into this new dynamical core. We will present high-resolution simulations of the dust and water cycles that show that sub-degree horizontal resolution improves the agreement between the vertical distribution of dust and water ice and observations. In particular, both water ice clouds and dust are transported to higher altitudes due to stronger topographic circulations at high resolution. Preliminary results suggest that high-resolution global modeling is needed to properly capture critical features of the dust and water cycles, and thus the current Mars climate.

Kahre, Melinda A.↗