Search NASA⌕ Search

SEARCH · Search NASA

Results for “production simulation codes”

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 55 records · Page 3

Solar-Reflective “Cool” Walls: Benefits, Technologies, and Implementation

Raising the albedo (solar reflectance) of a building’s walls reduces unwanted solar heat gain in the cooling season. This saves electricity and lowers peak power demand by decreasing the need for air conditioning. It can also cool the outside air, which can mitigate the urban heat island effect and also improve air quality by slowing the reactions that produce smog. This project quantified the energy savings, peak demand reduction, urban cooling, and air quality improvements attainable from solar-reflective “cool” walls in California; collaborated with industry to assess the performance of existing cool-wall technologies, and to develop innovative cool-wall solutions; and worked with state and federal government agencies, utilities, and industry to create a cool-wall infrastructure, including application guidelines, a product rating program, incentives, and building code credits. Simulations indicate that cool walls provide annual energy savings, peak demand reduction, annual emission reduction, and summer heat island mitigation benefits comparable to those yielded by cool roofs, and are helpful across California and in most of the southern half of the United States (that is, in U.S. climate zones 1—4). Natural exposure trials conducted at three sites in California and another three sites across the United States indicate that cool-wall materials tend to stay clean and reflective. Significant advances were made in novel cool-wall technologies, such as fluorescent cool pigments that expand the color palette for cool-wall products. We prepared guidelines for the climate- and building-appropriate use of cool walls, convened a stakeholder workshop, and created a working group. Ongoing efforts seek to introduce or expand cool-wall provisions in building energy standards, green building programs, and energy efficiency incentive programs, and to develop a cool-wall product rating system.

utility incentives↗

Initial Investigation of the Small Feature Detection Capability of the ORNL Fast Neutron Tomographic Imaging System

A variety of fast neutron imaging simulations and measurements were performed to investigate the spatial resolution and small feature detection capability of a version of the Oak Ridge National Laboratory fast neutron tomographic imaging system. Two image quality indicators were examined: (1) a variety of long steel pins embedded axially in the interior of a polyethylene cylinder, and (2) a variety of short steels pins embedded radially in the surface of a polyethylene cylinder. Measurements demonstrated detectability of about 3 mm. Further data analysis may show improved detection performance and allow a more precise determination of the spatial resolution. Our current simulation code takes only partial consideration of the DT generator neutron production spot size, with neutron production effectively originating at a single point. Simulations showed that features less than 2 mm should be detectable with an improved generator. Our current ING-27 DT generator has a neutron production spot size of about 3.5 mm, which appears to be the limiting factor at present. Generators with spot sizes as small as 1 mm are commercially available. Higher-power, small–spot size associated particle generators are at the prototype stage and are expected to be commercially available soon.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

MPAS-Ocean NESAP Status Report

NESAP performance improvements on MPAS-Ocean have resulted in a 5% to 7% speed-up on each of the examined systems including Cori-KNL, Cori-Haswell, and Edison. These tests were configured to emulate a production workload by using 128 nodes and a high-resolution ocean domain. Overall, the gap between standard and many-core architecture performance has been narrowed, but Cori-KNL remains considerably under-performing relative to Edison. NESAP code alterations affected 600 lines of code, and most of these improvements will benefit other MPAS codes (sea ice, land ice) that are also components within ACME. Modifications are fully tested within MPAS. Testing in ACME across many platforms is underway, and must be completed before the code is merged. In addition, a ten-year production ACME global simulation was conducted on Cori-KNL in late 2016 with the pre-NESAP code in order to test readiness and configurations for scientific studies. Next steps include assessing performance across a range of nodes, threads per node, and ocean resolutions on Cori-KNL.

54 ENVIRONMENTAL SCIENCES↗

Warthog Framework Update: Depletion Capabilities

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program the US Department of Energy (DOE) Office of Nuclear Energy supports two simulation frameworks to provide advanced simulation techniques for analyzing nuclear reactors. Primarily at Idaho National Laboratory (INL), the fuels product line (FPL) is developing the fuel performance code BISON from Hales et al. [2016], which depends on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework developed by Gaston et al. [2009]. Likewise, at Argonne National Laboratory (ANL), in the reactor product line (RPL), PROTEUS is being developed by Lee et al. [2013] and Nek5000 by Fischer et al. [2015] based on the Simulation-based High-efficiency Advanced Reactor Prototyping (SHARP) framework from Tautges et al. [2013]. MOOSE employs a bottom-up approach to multiphysics coupling, where component codes are built upon a common framework. SHARP employs a top-down approach in which legacy component codes are connected after being separately developed. Each approach has benefits and disadvantages, but this work focuses on having the two approaches work in harmony.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Towards exascale for wind energy simulations

We examine large-eddy-simulation modeling approaches and computational performance of two open-source computational fluid dynamics codes for the simulation of atmospheric boundary layer flows that are of direct relevance to wind energy production. The first code, NekRS, is a high-order, unstructured-grid, spectral element code. The second code, AMR-Wind, is a second-order, block-structured, finite-volume code with adaptive mesh refinement capabilities. The objective of this study is to co-develop these codes in order to improve model fidelity and performance for each. These features will be critical for running ABL-based applications such as wind farm analysis on advanced computing architectures. To this end, we investigate the performance of NekRS and AMR-Wind on the Oak Ridge Leadership Facility supercomputers Summit, using 4 to 800 nodes (24 to 4,800 NVIDIA V100 GPUs), and Crusher, the testbed for the Frontier exascale system, using 18 to 384 Graphics Compute Dies on AMD MI250X GPUs. We compare strong- and weak-scaling capabilities, linear solver performance, and time to solution. We also identify leading inhibitors to parallel scaling.

17 WIND ENERGY↗

Analysis of near-field and far-field aerosol dispersion for microreactors

The current paper presents a simulation-based analysis of aerosol dispersion in the near-field and far-field of a generic, conceptual microreactor operating at pressures close to the ambient pressure; therefore, in the event of an accident that causes radionuclide leakage from the microreactor containment, the radionuclide particles are less likely to travel too far from the reactor, as opposed to conventional reactors. Accordingly, the presented work provides estimates of average and 95-percentile values of the relative effluent concentration. A parametric study is then performed to narrow down the parameters which affect the aerosol dispersion characteristics most significantly. Simulations were performed in the computer code ARCON96, and the parameters found to affect aerosol dispersion characteristics are the atmospheric stability class, and the distance between the release point and the receptor. It is recommended that the computer code RADTRAD be used to calculate actual dosage over distance, using the outputs from ARCON96 as inputs, along with reactor-specific core term inventories. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

On producing CO 2 from subsurface reservoirs: simulations of liquid-gas phase change caused by decompression

Carbon dioxide (CO 2 ) extraction from deep reservoirs is currently important in CO 2 enhanced oil recovery (EOR) and may become more important in the future if interim CO 2 storage becomes common. In late 2014, we were involved in a production test of liquid CO 2 from the Middle Duperow dolostone at Kevin Dome, Montana. The test resulted in lowering the temperature at the well bottom to ~2 °C, and showed that the well and reservoir had very low CO 2 productivity. We have used the CO 2 modeling capabilities of the TOUGH codes to simulate the test and to show that liquid CO 2 in the reservoir changes to gas phase as the pressure is lowered in the well during production testing. The associated phase change and decompression combine to drastically lower the bottom-hole temperature, creating the potential for water ice or CO 2 hydrate to form. By hypothesizing a relatively high-permeability damage zone near the well surrounded by lower permeability reservoir rock, we can match the observed pressure, temperature, and production rate. Moving from the Kevin Dome test to the question of CO 2 extraction from deep reservoirs in general, we carried out a parametric study to investigate the effects of reservoir depth and transmissivity on CO 2 production rate for a prototypical reservoir. Simulations show that depth and high transmissivity favor productivity. Complex phase changes within the ranges of P-T encountered in typical CO 2 production wells affect production rates. The results of our parametric study may be useful for the preliminary feasibility assessment of CCO 2 extraction from deep reservoirs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of bulk composition on production rates of cosmogenic nuclides in meteorites

The bulk chemical composition of meteorites has been suggested as a main factor influencing the production of cosmogenic nuclides. Numerical simulations with Los Alamos Monte Carlo production and transport codes were done for Ne-21/Ne-22 ratios and Ar-38 production rates in meteorites with a wide range of compositions. The calculations show that an enhanced flux of low-energy secondary particles in metal-rich phases is the essential key for the explanation of experimentally observed differences in nuclide production processes in various meteorite classes.

Masarik, Jozef↗

A Two-Stage Procedure Toward the Efficient Implementation of PANS and Other Hybrid Turbulence Models

The main objective of this article is to introduce and to show the implementation of a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for Partial Averaged Navier-Stokes (PANS) and other hybrid models. It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The first step is to solve the unsteady or steady Reynolds Averaged Navier-Stokes (URANS/RANS) equations. From this preprocessing step, the turbulence length-scale field is obtained. This is then used to compute the characteristic length-scale ratio between the turbulence scale and the grid spacing. Based on this ratio, we can assess the finest scale resolution that a given grid for a given flow can support. Along with other additional criteria, we are able to analytically identify the appropriate hybrid solver resolution for different regions of the flow. This procedure removes the grid dependency issue that affects the results produced by different hybrid procedures in solving unsteady flows. The formulation, implementation methodology, and validation example are presented. We implemented this capability in a production Computational Fluid Dynamics (CFD) code, PAB3D, for the simulation of unsteady flows.

Abdol-Hamid, Khaled S.↗

A ROS-based Simulator for Testing the Enhanced Autonomous Navigation of the Mars 2020 Rover

In order to achieve the ambitious objectives of the Mars 2020 (M2020) mission, in particular the ability to autonomously traverse more challenging terrains more efficiently, new surface mobility software was developed for Enhanced Navigation (ENav). That decision was made early in the project, before most of the new surface flight software (FSW) existed, which created a need for a separate framework where the new navigation algorithms could be quickly prototyped and tested, before more realistic FSW-based testbeds became available. The JPL robotics team chose the Robot Operating System [1] (ROS) as the environment in which to test the new ENav algorithms. This made it possible to write the algorithms in the C language required by the FSW, so they could be directly ported over to the flight module later on, while leveraging all the C++ libraries and tools provided by ROS for simulation and testing. The ENav algorithms were developed as a separate C library, and stubs were used to replace any FSW-specific code, such as Event Reporting (EVRs) and data products (DPs). A ROS simulator was developed to generate a rich set of varied 3D terrains representative of the candidate Mars landing sites and simulate the physics of the rover motion, the point cloud perceived by the rover’s stereo vision system, and the new thinking-while-driving (TWD) navigation logic which directs the rover to drive autonomously to user-specified waypoints. To simulate the rover motion and perception, a ROS node was developed that uses a software library called HyperDrive Sim (HDSim), which is a wrapper for the Rover Sequencing and Visualization Program [2] (RSVP). That library provides roverterrain settling, realistic slip modelling, and camera rendering capability based on the rover’s NavCam machine vision models. To simulate the navigation logic, a ROS node was created that initializes and runs the ENav algorithms in a way that mimics the FSW execution, while also providing the capability to load and replay data products, including re-running the recorded inputs through the ENav algorithms for testing. An engineering Graphical User Interface (GUI) was also developed to visualize various elements, such as the rover pose during the drive, the simulated and perceived terrain, the selected local and global paths to the goal, the evaluated candidate paths and the reasons why they were rejected, the keep-in and keep-out zones (KIOZs), etc. Finally, an advanced Monte Carlo (MC) framework that can run many simulations in parallel on the Cloud and automatically generate reports that capture the key ENav performance metrics was developed to evaluate the system in a statisticallymeaningful way. This paper provides an overview of the ROSbased simulator used for testing the M2020 ENav algorithms.

Toupet, Olivier↗

Alternative CAES Technology Using Depleted Unconventional Gas Wells and Subsurface Thermal Energy Storage (GeoCAES)

This project assessed the technical viability of a process called GeoCAES. The process stores electrical energy by injecting natural gas into shale gas formations using a compressor, storing it, and producing it through an expander to generate electricity. This data submission includes the models of temperature and pressure changes in the wellbore, surface plant equipment (compressor and expander), and the code used in CMG GEM reservoir modeling software to simulate injection and production. Note - the wellbore and surface plant equipment models use the REFPROP Excel Add-in from NIST (linked in submission) to calculate natural gas properties. Note - the reservoir model code requires a license for the Computer Modeling Group (CMG) GEM reservoir modeling software (linked in submission) to run it.

15 GEOTHERMAL ENERGY↗

Deploy production sliding mesh capability with linear solver benchmarking.

Wind applications require the ability to simulate rotating blades. To support this use-case, a novel design-order sliding mesh algorithm has been developed and deployed. The hybrid method combines the control volume finite element methodology (CVFEM) with concepts found within a discontinuous Galerkin (DG) finite element method (FEM) to manage a sliding mesh. The method has been demonstrated to be design-order for the tested polynomial basis (P=1 and P=2) and has been deployed to provide production simulation capability for a Vestas V27 (225 kW) wind turbine. Other stationary and canonical rotating ow simulations are also presented. As the majority of wind-energy applications are driving extensive usage of hybrid meshes, a foundational study that outlines near-wall numerical behavior for a variety of element topologies is presented. Results indicate that the proposed nonlinear stabilization operator (NSO) is an effective stabilization methodology to control Gibbs phenomena at large cell Peclet numbers. The study also provides practical mesh resolution guidelines for future analysis efforts. Application-driven performance and algorithmic improvements have been carried out to increase robustness of the scheme on hybrid production wind energy meshes. Specifically, the Kokkos-based Nalu Kernel construct outlined in the FY17/Q4 ExaWind milestone has been transitioned to the hybrid mesh regime. This code base is exercised within a full V27 production run. Simulation timings for parallel search and custom ghosting are presented. As the low-Mach application space requires implicit matrix solves, the cost of matrix reinitialization has been evaluated on a variety of production meshes. Results indicate that at low element counts, i.e., fewer than 100 million elements, matrix graph initialization and preconditioner setup times are small. However, as mesh sizes increase, e.g., 500 million elements, simulation time associated with \setup-up" costs can increase to nearly 50% of overall simulation time when using the full Tpetra solver stack and nearly 35% when using a mixed Tpetra- Hypre-based solver stack. The report also highlights the project achievement of surpassing the 1 billion element mesh scale for a production V27 hybrid mesh. A detailed timing breakdown is presented that again suggests work to be done in the setup events associated with the linear system. In order to mitigate these initialization costs, several application paths have been explored, all of which are designed to reduce the frequency of matrix reinitialization. Methods such as removing Jacobian entries on the dynamic matrix columns (in concert with increased inner equation iterations), and lagging of Jacobian entries have reduced setup times at the cost of numerical stability. Artificially increasing, or bloating, the matrix stencil to ensure that full Jacobians are included is developed with results suggesting that this methodology is useful in decreasing reinitialization events without loss of matrix contributions. With the above foundational advances in computational capability, the project is well positioned to begin scientific inquiry on a variety of wind-farm physics such as turbine/turbine wake interactions.

97 MATHEMATICS AND COMPUTING↗

Deploy production sliding mesh capability with linear solver benchmarking (ECP Milestone Report, Ver. 1.0)

Wind applications require the ability to simulate rotating blades. To support this use-case, a novel design-order sliding mesh algorithm has been developed and deployed. The hybrid method combines the control volume finite element methodology (CVFEM) with concepts found within a discontinuous Galerkin (DG) finite element method (FEM) to manage a sliding mesh. The method has been demonstrated to be design-order for the tested polynomial basis (P=1 and P=2) and has been deployed to provide production simulation capability for a Vestas V27 (225 kW) wind turbine. Other stationary and canonical rotating flow simulations are also presented. As the majority of wind-energy applications are driving extensive usage of hybrid meshes, a foundational study that outlines near-wall numerical behavior for a variety of element topologies is presented. Results indicate that the proposed nonlinear stabilization operator (NSO) is an effective stabilization methodology to control Gibbs phenomena at large cell Peclet numbers. The study also provides practical mesh resolution guidelines for future analysis efforts. Application-driven performance and algorithmic improvements have been carried out to increase robustness of the scheme on hybrid production wind energy meshes. Specifically, the Kokkos-based Nalu Kernel construct outlined in the FY17/Q4 ExaWind milestone has been transitioned to the hybrid mesh regime. This code base is exercised within a full V27 production run. Simulation timings for parallel search and custom ghosting are presented. As the low-Mach application space requires implicit matrix solves, the cost of matrix reinitialization has been evaluated on a variety of production meshes. Results indicate that at low element counts, i.e., fewer than 100 million elements, matrix graph initialization and preconditioner setup times are small. However, as mesh sizes increase, e.g., 500 million elements, simulation time associated with "setup-up" costs can increase to nearly 50% of overall simulation time when using the full Tpetra solver stack and nearly 35% when using a mixed Tpetra- Hypre-based solver stack. The report also highlights the project achievement of surpassing the 1 billion element mesh scale for a production V27 hybrid mesh. A detailed timing breakdown is presented that again suggests work to be done in the setup events associated with the linear system. In order to mitigate these initialization costs, several application paths have been explored, all of which are designed to reduce the frequency of matrix reinitialization. Methods such as removing Jacobian entries on the dynamic matrix columns (in concert with increased inner equation iterations), and lagging of Jacobian entries have reduced setup times at the cost of numerical stability. Artificially increasing, or bloating, the matrix stencil to ensure that full Jacobians are included is developed with results suggesting that this methodology is useful in decreasing reinitialization events without loss of matrix contributions. With the above foundational advances in computational capability, the project is well positioned to begin scientific inquiry on a variety of wind-farm physics such as turbine/turbine wake interactions.

Domino, Stefan P.↗

Electron Thermalization in the Solar Wind and Planetary Plasma Boundaries

The work carried out under this contract attempts a better understanding of whistler wave generation and associated scattering of electrons in the solar wind. This task is accomplished through simulations using a particle-in-cell code and a Vlasov code. In addition, the work is supported by the utilization of a linear kinetic dispersion solver. Previously, we have concentrated on gaining a better understanding of the linear mode properties, and have tested the simulation codes within a known parameter regime. We are now in a new phase in which we implement, execute, and analyze production simulations. This phase is projected to last over several reporting periods, with this being the second cycle. In addition, we have started to research to what extent the evolution of the pertinent instabilities is two-dimensional. We are also continuing our work on the visualization aspects of the simulation results, and on a code version that runs on single-user Alpha-processor based workstations.

Krauss-Varban, Dietmar↗

Deployment of the HFIRCON transport and depletion tool for plutonium-238 production studies

Irradiation of {sup 237}Np-bearing targets in Oak Ridge National Laboratory's (ORNL) High Flux Isotope Reactor (HFIR) results in the efficient production of {sup 238}Pu, which, in the form of heat source PuO{sub 2}, is used as a reliable power source for deep-space and planetary NASA missions. A technology demonstration subproject was initiated at ORNL in 2011 to develop and implement the technology required to establish a {sup 238}Pu supply chain. A systematic progression of NpO{sub 2}/Al cermet (20 vol.% NpO{sub 2}) activities to date has successfully demonstrated target fabrication, irradiation, and chemical recovery processes. Recent program tasks have included the development of the HFIRCON transport and depletion tool for efficient reactor physics analyses and the evaluation of increased NpO{sub 2} loadings (i.e., beyond 20 vol.%) and NpN-based targets. This paper documents the deployment of the HFIRCON code to assess various Np concentrations in NpO{sub 2}- and NpN-based targets in HFIR's inner small vertical experiment facilities. Results indicate that {sup 238}Pu production and quality can be enhanced with increased Np loadings; however, target conversion rates are reduced. The results recorded in this paper, thermal and material balance evaluations, and testing requirement planning will be used to determine whether increased NpO{sub 2} loadings or NpN-based targets will be further considered. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Characterization of W production during ICRF operations: experiments and modeling

For successfully heating plasma with waves in the ion cyclotron range of frequencies (ICRFs), mitigating impurity production is just as crucial as maximizing power coupling, especially in high-Z environments. ICRF can effectively deposit energy on ions, modify turbulence-driven transport, and enhance fusion reaction efficiency, but only when its power coupling has minimal impact on impurity production. To do so, one must rely on a toroidal array of at least three active elements excited with appropriate phasing and power ratio to reduce the currents induced on the antenna frame below levels critical for physical sputtering. In contrast to classic two-strap antennas, which are optimized for dipole phasing with equal power on both straps, three-strap antennas in ASDEX Upgrade (AUG)—but also four-strap antennas in JET, Alcator C-Mod, SPARC and ITER—offer the possibility to act also on the power ratio between the central and outer straps. With optimal settings, impurity production can be reduced substantially, making the ICRF compatible with the high-Z wall. This paper explores the characteristics of the AUG three-strap antennas in terms of impurity production, as well as the key role of plasma composition in this process. Numerical simulations were performed using SSWICH and Petra-M (finite element codes) to quantify impurity production and compare with experimental results. Energies of ions falling on antenna limiters (measured with probes) are well predicted by both codes. These tools are then used to further describe the source of the impurity, namely the gross erosion of tungsten from an ICRF antenna, for different plasma mixtures. Results are also compared to spectroscopy data. Ultimately, we show that deleterious effects of the ICRF on plasma surface interactions will be weaker in plasmas containing larger fractions of highly ionized heavier low-Z impurity, which is typically relevant for experiments relying on impurity seeding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

LArSoft: Toolkit for Simulation, Reconstruction and Analysis of Liquid Argon TPC Neutrino Detectors

LArSoft is a set of detector-independent software tools for the simulation, reconstruction and analysis of data from liquid argon (LAr) neutrino experiments The common features of LAr time projection chambers (TPCs) enable sharing of algorithm code across detectors of very dierent size and conguration. LArSoft is currently used in production simulation and reconstruction by the ArgoNeuT, DUNE, LArIAT, MicroBooNE, and SBND experiments. The software suite oers a wide selection of algorithms and utilities, including those for associated photo-detectors and the handling of auxiliary detectors outside the TPCs. Available algorithms cover the full range of simulation and reconstruction, from raw waveforms to high- level reconstructed objects, event topologies and classication. The common code within LArSoft is contributed by adopting experiments, which also provide detector-specic geometry descriptions, and code for the treatment of electronic signals. LArSoft is also a collaboration of experiments, Fermilab and associated software projects which cooperate in setting requirements, priorities, and schedules. In this talk, we outline the general architecture of the software and the interaction with external libraries and detector-specic code. We also describe the dynamics of LArSoft software development between the contributing experiments, the projects supporting the software infrastructure LArSoft relies on, and the core LArSoft support project.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

How numerical simulations helped to achieve breakeven on the NIF

The inertial confinement fusion program relies upon detailed simulations with inertial confinement fusion (ICF) codes to design targets and to interpret the experimental results. These simulations treat as much physics from essential principles as is practical, including laser deposition, cross beam energy transfer, x-ray production and transport, nonlocal thermal equilibrium kinetics, thermal transport, hydrodynamic instabilities, thermonuclear burn, and transport of reaction products. Improvements in radiation hydrodynamic code capabilities and vast increases in computing power have enabled more realistic, accurate 3D simulations that treat all known asymmetry sources. We describe how numerical simulations helped to guide the program, assess the impediments to breakeven, and optimize every aspect of target design. A preshot simulation of the first National Ignition Facility experiment that surpassed breakeven predicted an increased yield that matches the experimental result, within the preshot predicted uncertainty, with a target gain of 1.5. We will cover the key developments in Lawrence Livermore National Laboratory ICF codes that enabled these simulations and give specific examples of how they helped to guide the program.

Marinak, M. M. (ORCID:0009000331127696)↗