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 37 records · Page 2

Initial full core SMR simulations with NekRS

This document describes the completion of a recent milestone of the ExaSMR program concerning "Full core simulations". The ExaSMR project is developing tools for coupling Monte Carlo (MC) radiation transport solvers to a computational fluid dynamics (CFD) solver. Work is based on the Shift MC, OpenMC MC, and Nek5000/NekRS CFD codes. As part of this milestone, a novel set of pin-resolved CFD full-core simulations have been performed for the first time. These simulations represent a significant increase in capability in what is now possible with CFD on pre-Exascale systems. The simulations have been performed with the spectral element solvers NekRS on the supercomputer Summit. Multiple simulations campaigns have been conducted: (1) LES simulations in bare bundle, (2) RANS simulations in bare bundles, (3) RANS simulations with momentum sources and (4) Conjugate heat transfer calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multidimensional Simulations of Core Convection

The cores of main sequence intermediate- and high-mass stars are convective. Mixing at the radiative–convective boundary, waves excited by the convection, and magnetic fields generated by convective dynamos all influence the main sequence and post-main sequence evolution of these stars. These effects must be understood to accurately model the structure and evolution of intermediate- and high-mass stars. Unfortunately, there are many challenges in simulating core convection due to the wide range of temporal and spatial scales, as well as many important physics effects. In this review, we describe the latest numerical strategies to address these challenges. We then describe the latest state-of-the-art simulations of core convection, summarizing their main findings. These simulations have led to important insights into many of the processes associated with core convection. Two outstanding problems with multidimensional simulations are, 1. it is not always straightforward to extrapolate from simulation parameters to the parameters of real stars; and 2. simulations using different methods sometimes appear to arrive at contradictory results. To address these issues, next generation simulations of core convection must address how their results depend on stellar luminosity, dimensionality, and turbulence intensity. Furthermore, code comparison projects will be essential to establish robust parameterizations that will become the new standard in stellar modeling.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparison of simulation results with sea-level experimental data on 10(14) - 10(16) air shower cores

Simulation predictions for the Leeds 35 sq m horizontal discharge chamber array for proton primaries with a approx. E sup 2.7 spectrum extrapolated from balloon data to 10 to the 16th power eV give power law rho (r)-spectra with constant slope approx. -2 consistent with the experimental data up to the point at which they steepen but overshooting them at higher densities, and at high shower sizes predicted cores which are significantly steeper than those observed. Further comparisons with results for heavy nuclei primaries (up to A = 56) point to the inadequacy of changes in primary composition to account for the observed density spectra and core flattening, and the shower size spectrum together, and point, therefore, to the failure of the scaling interaction model at approx. 10 to the 15th power eV primary energy.

Ash, A. G.↗

A study of the core module simulator floor capability

The floor of the Core Module Simulator (CMS) is required to support various combinations of dead load and live load during the testing process. Even though there is published data on the structural capability of the grating, it is not always evident if the combined loadings with joint loads will cause structural failure.

Foreman, James W.↗

Data accompanying the publication "Sensitivity of grass fires burning in marginal conditions to atmospheric turbulence" in the Journal of Geophysical Research-Atmospheres

Data accompanying Jonko, Yedinak, Conley, and Linn (2021): "Sensitivity of grass fires burning in marginal conditions to atmospheric turbulence", JGR-Atmospheres. LA-UR-21-24832. The following files were generated from simulations performed with Los Alamos National Laboratory coupled atmosphere-fire behavior model FIRETEC v1.3: -area_burned.csv contains the area burned per second in m^2 for 45 CORE simulations, 20 INIT simulations and 20 BOUND simulations (85 columns), and 600 seconds of simulation (600 rows). -fire_perimeter.csv containing the fire perimeters in m for 45 CORE simulations, 20 INIT simulations and 20 BOUND simulations (85 columns), and 600 seconds of simulation (600 rows). -consumption_per_second.csv containing consumption in kg/s for 45 CORE simulations, 20 INIT simulations and 20 BOUND simulations (85 columns), and 600 seconds of simulation (600 rows). -rate_of_spread.csv containing the distance of the fire front from the ignition line in m for 45 CORE simulations, 20 INIT simulations and 20 BOUND simulations (85 columns), and 600 seconds of simulation (600 rows). -CORE1_u_10m_300_sec.hdf, INIT3_u_10m_300_sec.hdf, and BOUND3_u_10m_300_sec.hdf contain u wind speeds at 10 m height for 300 seconds (3D: time, x, y) of simulation for simulations CORE1, INIT3, and BOUND3, respectively. -CORE.final.fuel.density.hdf contains final fuel densities in the surface layer of the model (3D: ensemble member, x, y) for all 45 CORE simulations. -CORE{1-45}_3D_u_120_sec.hdf: 45 files containing vertically and horizontally resolved u winds for 120 seconds of simulation (4D: time, x, y, z) for the 45 CORE ensemble members. For additional information about this data, please contact the corresponding author at ajonko@lanl.gov.

54 ENVIRONMENTAL SCIENCES↗

Heat‐Generating Nuclear Waste in Salt: Field Testing and Simulation

Core Ideas A field‐scale experiment and numerical simulations confirm salt backfill behavior. Simulations closely match temperature around and under the piled salt backfill. Results indicate limited dissolution–precipitation reactions around the heat source. Alteration of backfill is unlikely if the drift is allowed to dry before emplacement. Investigations relating to in‐drift disposal of heat‐generating nuclear waste in salt have raised questions about heat–brine interactions in the unsaturated run‐of‐mine (RoM) salt pile used as backfill. These interactions have the potential to change the structure of the RoM salt surrounding the canister, possibly altering long‐term containment of the source. An experiment is in progress at the Waste Isolation Pilot Plant (WIPP), New Mexico, in which a heated canister was placed on the floor of an open drift, covered in a pile of RoM salt, and energized with 1000 W. Temperature in the RoM salt pile had stabilized after about 15 d, allowing evaluation of the heat‐up period of the ongoing experiment. Using a multiphase porous flow simulator that has been modified to handle salt‐specific coupled processes, we examined coupled thermal–hydrological–chemical behavior in the RoM salt pile. Our simulations suggest that for the relatively dry cases examined, porosity changes within RoM salt in a generic salt repository are likely to be minor in the period between waste emplacement and plastic closure of the drift. The primary sensitivity for porosity change is to the early moisture content of the RoM salt used to cover the canister. Secondary influences include moisture availability from the disturbed rock zone (DRZ) surrounding the drift and the capillary pressure ratio between the DRZ and the RoM salt. Early changes in porosity and permeability may be affected by moisture content, but this was not observed in the test. Such changes would be most likely to occur when using damp RoM salt or if waste is emplaced in a drift immediately following opening of the drift before evaporative dewatering of the drift walls occurs.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An Update on Cardinal: Toward Full Core Pebble Simulations

The new DOE-NE Center of Excellence for Thermal-fluids applications in Nuclear Energy inaugurated in April 2018 considers and researches novel new solution strategies for historically challenging flow issues that still plague the current fleet of deployed Light Water nuclear Reactors (LWRs) as well as predicting various fluid flow and fluid related issues with advanced reactor technologies. Our advanced thermal-fluids research and development approach synergistically combines three natural, though overlapping, length and time scales in a hierarchal multi-scale approach to avoid the temptation and pitfalls of attempting to develop a single solve-all algorithm for physical fluid flow problems that will span nine orders of magnitude in spatial and temporal scales. To demonstrate the multi-scale philosophy of the center we focus on Fluoride Cooled High Temperature Reactors (FHRs), and in particular on the Berkley’s PB-FHR Mark-I design. The Fluoride salt cooled High temperature Reactor (FHR) is a class of advanced nuclear reactors that combine the robust coated particle fuel form from high temperature gas cooled reactors, direct reactor auxiliary cooling system (DRACS) passive decay removal of liquid metal fast reactors, and the transparent, high volumetric heat capacitance liquid Fluoride salt working fluids - such as Flibe - from molten salt reactors. This combination of fuel and coolant enables FHRs to operate in a high-temperature low-pressure design space that has beneficial safety and economic implications. The PB-FHR reactor relies on a pebble bed approach and pebble bed reactors are in a sense the poster child for multi-scale analysis. The lower length-scale simulator for pebble reactor cores comprises three physics: neutronics (OpenMC), thermal-fluids (Nek5000/NekRS) and fuel performance (BISON). As part of Center of Excellence ongoing research efforts we have developed Cardinal, a new tool platform for lower length-scale simulation. Cardinal tightly couples all three physics and leverages advances in MOOSE such as the MultiApp system and the concept of MOOSE-wrapped Apps. The present report provides an update on the development of Cardinal with an extension of Cardinal to GPUs. We also perform a first-of-a-kind demonstration simulation on Summit representing a 10x capability increase in terms of pebble count for Cardinal.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Scalable algorithms for domain decomposition of Monte Carlo neutral particle transport simulations on unstructured mesh

Monte Carlo simulations are the gold standard in particle transport simulations, for shielding and nuclear reactor simulations Full core simulations can require terabytes of memory for tracking pin-resolved (up to 50k pins), nuclide (200 tracked), axial (20-200 zones) reaction rates Domain decomposition, splitting the calculation and the memory storage, among several nodes is challenging as it involves communication of particles across domain boundaries We created a new MOOSE application, MaCaw, to study this problem Algorithms developed for ray tracing in MOOSE were adapted for simulating neutral particle transport Physics and tally routines in OpenMC, dynamically linked, are called from the application

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Full core 3D simulation of the BEAVRS benchmark with OpenMOC

The Method of Characteristics (MOC) has seen wide interest in reactor physics because of its accuracy and efficiency in computing lattice physics problems. While most of its use has been in solving 2D problems, there has been recent interest in extending MOC to 3D in order to more accurately calculate 3D power distributions in LWRs. While the method is naturally extensible to 3D, it presents significant computational difficulties. OpenMOC, an open source MOC simulator, has been extended from 2D to 3D capabilities. Simulation results are presented in this study for the full core 3D simulation of the BEAVRS benchmark, adequately resolved in space and angle.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Laser anemometry measurements of natural circulation flow in a scale model PWR reactor system

The natural circulation of a single phase fluid in a scale model of a pressurized water reactor system during a postulated grade core accident is analyzed. The fluids utilized were water and SF6. The design of the reactor model and the similitude requirements are described. Four LDA tests were conducted: water with 28 kW of heat in the simulated core, with and without the participation of simulated steam generators; water with 28 kW of heat in the simulated core, with the participation of simulated steam generators and with cold upflow of 12 lbm/min from the lower plenum; and SF6 with 0.9 kW of heat in the simulated core and without the participation of the simulated steam generators. For the water tests, the velocity of the water in the center of the core increases with vertical height and continues to increase in the upper plenum. For SF6, it is observed that the velocities are an order of magnitude higher than those of water; however, the velocity patterns are similar.

Kadambi, J. R.↗

An Assessment of Southern Ocean Water Masses and Sea Ice During 1988-2007 in a Suite of Interannual CORE-II Simulations

We characterise the representation of the Southern Ocean water mass structure and sea ice within a suite of 15 global ocean-ice models run with the Coordinated Ocean-ice Reference Experiment Phase II (CORE-II) protocol. The main focus is the representation of the present (1988-2007) mode and intermediate waters, thus framing an analysis of winter and summer mixed layer depths; temperature, salinity, and potential vorticity structure; and temporal variability of sea ice distributions. We also consider the interannual variability over the same 20 year period. Comparisons are made between models as well as to observation-based analyses where available. The CORE-II models exhibit several biases relative to Southern Ocean observations, including an underestimation of the model mean mixed layer depths of mode and intermediate water masses in March (associated with greater ocean surface heat gain), and an overestimation in September (associated with greater high latitude ocean heat loss and a more northward winter sea-ice extent). In addition, the models have cold and fresh/warm and salty water column biases centred near 50 deg S. Over the 1988-2007 period, the CORE-II models consistently simulate spatially variable trends in sea-ice concentration, surface freshwater fluxes, mixed layer depths, and 200-700 m ocean heat content. In particular, sea-ice coverage around most of the Antarctic continental shelf is reduced, leading to a cooling and freshening of the near surface waters. The shoaling of the mixed layer is associated with increased surface buoyancy gain, except in the Pacific where sea ice is also influential. The models are in disagreement, despite the common CORE-II atmospheric state, in their spatial pattern of the 20-year trends in the mixed layer depth and sea-ice.

Sea ice↗

Development of Preliminary VERA-CS Crud Induced Localized Corrosion Modeling Capability

The Consortium for Advanced Simulation of Light Water Reactors (CASL) is developing multiphysics core-simulator software for light water reactors (LWRs) known as VERA-CS in order to improve the state of the art in modeling and simulation of challenge problems that are limiting to the nuclear industry. One of these challenge problems includes fuel rod crud deposition, which can lead to crud-induced power shift (CIPS) and crud-induced localized corrosion (CILC). This paper documents work that was performed to develop a preliminary CILC-modeling capability in VERA-CS in support of the crud challenge problem. The CILC capabilities were developed by coupling VERA-CS to the CASL-developed Cicada code, which provides 1D and 3D clad conduction and oxide growth modeling tools, as well as coupling to the CASL-developed MAMBA code, which is used for modeling clad crud deposition. An approach called rod thermal-hydraulic reconstruction (ROTHCON) was developed and integrated into VERA-CS. This allows the modeler to capture spacer-grid turbulence and heat transfer effects in the CTF subchannel code so that the spatial resolution of crud and oxide rod surface growth could be better resolved. After implementing these capabilities, several assessments were performed to ensure that the capabilities function as expected, and a pin-resolved quarter-core simulation was run as a demonstration.

Salko Jr, Robert↗

Verification of Reduced Order Modeling Based Uncertainty/Sensitivity Estimator (ROMUSE)

This article presents a number of verification case studies for a recently developed sensitivity/uncertainty code package. The code package, ROMUSE (Reduced Order Modeling based Uncertainty/Sensitivity Estimator) is an effort to provide an analysis tool to be used in conjunction with reactor core simulators, in particular the Virtual Environment for Reactor Applications (VERA) core simulator. ROMUSE has been written in C++ and is currently capable of performing various types of parameter perturbations and associated sensitivity analysis, uncertainty quantification, surrogate model construction and subspace analysis. The current version 2.0 has the capability to interface with the Design Analysis Kit for Optimization and Terascale Applications (DAKOTA) code, which gives ROMUSE access to the various algorithms implemented within DAKOTA, most importantly model calibration. The verification study is performed via two basic problems and two reactor physics models. The first problem is used to verify the ROMUSE single physics gradient-based range finding algorithm capability using an abstract quadratic model. The second problem is the Brusselator problem, which is a coupled problem representative of multi-physics problems. This problem is used to test the capability of constructing surrogates via ROMUSE-DAKOTA. Finally, light water reactor pin cell and sodium-cooled fast reactor fuel assembly problems are simulated via SCALE 6.1 to test ROMUSE capability for uncertainty quantification and sensitivity analysis purposes.

97 MATHEMATICS AND COMPUTING↗

Gyrokinetic Particle Simulations of Core Turbulence in Fusion Plasmas

In this project, we have further developed GTC gyrokinetic particle simulations to study plasma confinement properties in the stellarators and tokamaks with 3D magnetic fields, where both neoclassical and turbulent transport are self-consistently treated on an equal footing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Drive turbine system for 20-inch turbofan simulators. 2: Core turbine design

Drive-turbines for a given set of 20-inch turbo-fan simulators are described. The simulators had both single-stage and two-stage fans that had design pressure ratios as low as 1.25 and as high as 3.0. The desired objective of the study was to be able to drive all of the single-stage fans with one core turbine and to drive all of the two-stage fans with this same core turbine in combination with a duct turbine. The core turbine is described. Included are the design operating conditions, design velocity diagram and a power-speed envelope determined by an off-design performance procedure. Also discussed is the adaption and scaling of an existing turbine design to this particular application.

Whitney, W. J.↗

An Assessment of Antarctic Circumpolar Current and Southern Ocean Meridional Overturning Circulation During 1958-2007 in a Suite of Interannual CORE-2 Simulations

In the framework of the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II), we present an analysis of the representation of the Antarctic Circumpolar Current (ACC) and Southern Ocean Meridional Overturning Circulation (MOC) in a suite of seventeen global ocean-sea ice models. We focus on the mean, variability and trends of both the ACC and MOC over the 19582007 period, and discuss their relationship with the surface forcing. We aim to quantify the degree of eddy saturation and eddy compensation in the models participating in CORE-II, and compare our results with available observations, previous fine-resolution numerical studies and theoretical constraints. Most models show weak ACC transport sensitivity to changes in forcing during the past five decades, and they can be considered to be in an eddy saturated regime. Larger contrasts arise when considering MOC trends, with a majority of models exhibiting significant strengthening of the MOC during the late 20th and early 21st century. Only a few models show a relatively small sensitivity to forcing changes, responding with an intensified eddy-induced circulation that provides some degree of eddy compensation, while still showing considerable decadal trends. Both ACC and MOC interannual variability are largely controlled by the Southern Annular Mode (SAM). Based on these results, models are clustered into two groups. Models with constant or two-dimensional (horizontal) specification of the eddy-induced advection coefficient show larger ocean interior decadal trends, larger ACC transport decadal trends and no eddy compensation in the MOC. Eddy-permitting models or models with a three-dimensional time varying show smaller changes in isopycnal slopes and associated ACC trends, and partial eddy compensation. As previously argued, a constant in time or space is responsible for a poor representation of mesoscale eddy effects and cannot properly simulate the sensitivity of the ACC and MOC to changing surface forcing. Evidence is given for a larger sensitivity of the MOC as compared to the ACC transport, even when approaching eddy saturation. Future process studies designed for disentangling the role of momentum and buoyancy forcing in driving the ACC and MOC are proposed.

Southern Ocean dynamics↗

Development of Monte Carlo Capability for Orion Parachute Simulations

Parachute test programs employ Monte Carlo simulation techniques to plan testing and make critical decisions related to parachute loads, rate-of-descent, or other parameters. This paper describes the development and use of a MATLAB-based Monte Carlo tool for three parachute drop test simulations currently used by NASA. The Decelerator System Simulation (DSS) is a legacy 6 Degree-of-Freedom (DOF) simulation used to predict parachute loads and descent trajectories. The Decelerator System Simulation Application (DSSA) is a 6-DOF simulation that is well suited for modeling aircraft extraction and descent of pallet-like test vehicles. The Drop Test Vehicle Simulation (DTVSim) is a 2-DOF trajectory simulation that is convenient for quick turn-around analysis tasks. These three tools have significantly different software architectures and do not share common input files or output data structures. Separate Monte Carlo tools were initially developed for each simulation. A recently-developed simulation output structure enables the use of the more sophisticated DSSA Monte Carlo tool with any of the core-simulations. The task of configuring the inputs for the nominal simulation is left to the existing tools. Once the nominal simulation is configured, the Monte Carlo tool perturbs the input set according to dispersion rules created by the analyst. These rules define the statistical distribution and parameters to be applied to each simulation input. Individual dispersed parameters are combined to create a dispersed set of simulation inputs. The Monte Carlo tool repeatedly executes the core-simulation with the dispersed inputs and stores the results for analysis. The analyst may define conditions on one or more output parameters at which to collect data slices. The tool provides a versatile interface for reviewing output of large Monte Carlo data sets while preserving the capability for detailed examination of individual dispersed trajectories. The Monte Carlo tool described in this paper has proven useful in planning several Crew Exploration Vehicle parachute tests.

Moore, James W.↗