Pele Combustion: Early Frontier Runs and Beyond
This presentation reviews the challenges we encountered porting our combustion codes to the Frontier computer system, work sponsored under the Exascale Computing Program.
codes↗
SEARCH · Search NASA
Results for “Computer programming”
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.
Spherical Radiative Transfer in C++ (SRTC++): A Parallel Monte Carlo Radiative Transfer Model for Titan
We present a new computer program, SRTC++, to solve spatial problems associated with explorations of Saturn’s moon Titan. The program implements a three-dimensional structure well-suited to addressing shortcomings arising from plane-parallel radiative transfer approaches. SRTC++'s design uses parallel processing in an object-oriented, compiled computer language (C++) leading to a flexible and fast architecture. We validate SRTC++ using analytical results, semianalytical radiative transfer expressions, and an existing Titan plane-parallel model. SRTC++ complements existing approaches, addressing spatial problems like near-limb and near-terminator geometries, non-Lambertian surface phase functions (including specular reflections), and surface albedo nonuniformity.
97 MATHEMATICS AND COMPUTING↗
JaqalPaq v.1.1
Jaqal-pup (Jaqal Programming Utilities Project) is a high-level Python metaprogramming interface for Jaqal, a quantum programming language for near-term quantum hardware (see https://arxiv.org/abs/2003.09382 for a description of Jaqal). In addition to parsing and writing out Jaqal programs, it includes the ability to read in and write out quantum programs in other languages and frameworks, e.g. Qiskit (IBM), Quil (Rigetti), Cirq (Google), Q# (Microsoft), ProjectQ (ETH Zurich), and tket (Cambridge Quantum Computing). It implements a gate scheduler, an emulator of a small-scale quantum computer programmed using Jaqal, and tie-ins for other simulation back-ends. SAND2020-12710 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
Russo, Antonio↗
Water Quality for High Performance Computing
To support future generations of high performance computers, the Advanced Simulation and Computing Program has a mission need for 80MW of peak power and cooling. Currently the water for the cooling towers at the Strategic Computing Complex (SCC) are cooled by the LANL waste water treated Sanitary Effluent Recovery Facility (SERF) or potable water from the county. The SCC cooling towers currently operate by cycling water between 2-3 cycles of concentration (COC), which is less than the ideal 6 cycles for maximizing water usage. SERF is the ideal water supplier, however when it’s not available potable water is used. The amount of water that the SCC needs nominally is around 45 Million-Gallons/Year (MGY). The future super computer infrastructure (FSI) will need around 125 MGY nominally. SERF is approaching its production limit for reclaimed water. In order to avoid using 100% potable water for the FSI, the option of using non-potable reuse water from Los Alamos County (LAC) was analyzed. This project scope is to focus on the water quality needed for the new super computer cooling towers. Once water quality needs are established, the ideal treatment process needs to be determined. Due to high silica content in water at Los Alamos, emphasis has been placed on treatment processes to remove silica. SERF’s current and future ability to treat water was assessed, as well as the cost of running a SERF like facility or alternative.
97 MATHEMATICS AND COMPUTING↗
Two new SciDAC institutes promote mathematical tools and software technology for high-performance computing
Bigger is often said to be better, and the newest extreme-scale computers certainly are bigger, with millions of processing units. Moreover, the breadth of science performed on the U.S. Department of Energy (DOE) computing facilities is expanding, with new technology such as artificial intelligence emerging. These advances are exciting, creating new opportunities for scientific discovery; however, they also raise new questions for scientists who want to exploit these advances for tackling more complex problems. Will my simulation code be able to utilize the accelerators in extreme-scale computing systems? Can I take advantage of the deepening memory hierarchy in heterogeneous processors? Is there a way around bottlenecks caused by the widening ratio of peak floating-point operations per second to I/0 bandwidth? How can I manage my huge amounts of data effectively? Can I analyze data in situ, or must I transfer it to offline storage for later analysis? To address such questions, DOE announced that it is providing $57.5 million over the next five years for two multidisciplinary teams — FASTMath and RAPIDS2 — to develop new tools and techniques to harness supercomputers for scientific discovery. The teams, called SciDAC Institutes, are part of the Scientific Discovery through Advanced Computing program.
97 MATHEMATICS AND COMPUTING↗
Design of a Do-It-Yourself Geiger-Muller Counter With Smartphone Mapping Application
Legacy Geiger-Muller (GM) survey meters recovered from fallout shelters have been used by several nuclear scientific societies as part of high school outreach programs. A donated antique instrument helps teachers demonstrate radiological principles, but fails to develop student's electronics skills, generate excitement for nuclear careers, or provide individuals with their own devices to explore the radioactive planet. A simple, affordable GM survey meter built by each student would increase direct engagement while providing hands-on experience with circuit-building, soldering, and computer programming. The inclusion of an affordable single-board computer as a component in the survey meter would enable students to tackle more various computer science and electronics projects, thereby potentially recruiting more students into technology and engineering. This paper details the challenges faced by an interdisciplinary undergraduate team designing an easy-to-assemble smart GM survey meter. Their iterative research, design, and testing process included modification to a basic circuit to enable use of different tube types, component cost reduction, application development, and data communication. The ultimate product of the team's efforts, a survey meter with affordable components and a smartphone application capable of creating radiation maps, is detailed in full.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Molecular-level origin of the carboxylate head group response to divalent metal ion complexation at the air-water interface
We exploit gas-phase cluster ion techniques to provide insight into the local interactions underlying divalent metal ion-driven changes in the spectra of carboxylic acids at the air-water interface. This information is used to clarify the experimental findings of surface sensitive spectroscopies that the CO stretching bands of long chain acids appear at very similar energies when the head group is deprotonated by high sub-phase pH or is exposed to relatively high concentrations of Ca2+ divalent metal ions (e.g., Ca2+ and Mg2+). To this end, we report the evolution of the vibrational spectra of size-selected [Ca2+·RCO2¯]+·(H2O)n=0-14 and RCO2¯·(H2O)n=0-14 cluster ions obtained using infrared photofragmentation mass spectrometry toward the features observed at the air-water interface. The carboxylate with R = -CD2CD3 was chosen in order to capture the longer chain behavior and to suppress band congestion that complicates the spectra of the light isotopologue. Surprisingly, not only does stepwise hydration of the RCO2¯ anion and the [Ca2+·RCO2¯]+ contact ion pair yield solvatochromic responses in opposite directions, but in both cases the responses of the two CO bands (?_s^COO and ?_as^COO) to hydration are opposite to each other! The result is that the two CO bands evolve toward their interfacial and bulk asymptotes from opposite directions. Theoretical simulations of the [Ca2+·RCO2¯]+·(H2O)n clusters indicate that the metal ion remains directly bound to the head group in a contact ion pair motif as the asymmetric CO stretch converges at the bulk value by n = 12, establishing that direct metal complexation can account for the interfacial and bulk behavior. We discuss these effects in the context of a qualitative model that invokes the water network-dependent local electric field along the C-C bond connecting the head group to the tail as the key parameter driving the observed trends. M.A.J. and H.C.A. thank the National Science Foundation through the NSF Center for Aerosol Impacts on Chemistry of the Environment (NSF-CAICE), CHE-1305427 for supporting this work. M.A.J. thanks the Air Force Office of Scientific Research for the cryogenic spectrometer used in the gas phase experiments. K.D.J. thanks the Department of Energy under grant no. DE-FG02-00ER15066. C.J. Mundy and S.M. Kathmann are supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences. M.D. Baer was supported by the U.S. Department of Energy, Office of Basic Energy Sciences, Biomolecular Materials Program at PNNL. PNNL is operated by Battelle Memorial Institute for the U.S. Department of Energy under Contract No. DE-AC05- 76RL01830. Computing resources for the calculation of harmonic spectra were allocated by PNNL’s Institutional Computing program.
Denton, Joanna↗
Design and Testing for a New Thermosyphon Irradiation Vehicle
The High Flux Isotope Reactor (HFIR) at the Oak Ridge National Laboratory (ORNL) requires most materials and all fuel experiments to be placed in a pressure containment vessel to ensure that internal contaminants such as fission products cannot be released into the primary coolant. It also requires that all experiments be capable of withstanding various accident conditions (e.g., loss of coolant) without generating vapor bubbles on the surface of the experiment in the primary coolant. These requirements are intended to artificially increase experiment temperatures by introducing a barrier between the experimental materials and the HFIR coolant, and by reducing heat loads to the HFIR primary coolant, thus ensuring that no boiling can occur. A proposed design for materials irradiation would remove these limitations by providing the required primary containment with an internal cooling flow. This would allow for experiments to be irradiated without concern for coolant contamination (e.g., from cladding failure of advanced fuel pins) or for specimen heat load. This report describes a new materials irradiation experiment design that uses a thermosyphon cooling system to allow experimental materials direct access to a liquid coolant. The new design also increases the range of conditions that can be tested in HFIR. This design will provide a unique capability to validate the performance of current and advanced fuels and materials. Because of limited supporting data for this kind of irradiation vehicle, a test program was initiated to obtain operating data that can be used to (1) qualify the vehicle for operation in HFIR and (2) validate computer models used to perform design- and safety-basis calculations. This report also describes the test facility and experimental data, and it provides a comparison of the experimental data to computer simulations. A total of 51 tests have been completed: four tests with pure steam, 12 tests with argon, and 35 tests with helium. A total of 10 tests were performed at subatmospheric pressure, and four of these were performed with pure steam. One test was conducted at a high power of 92.7 kW, six tests were HFIR startups, and two tests were HFIR loss of offsite power (LOOP). Pressures up to 10 MPa, vapor temperatures up to 583 K (310°C), and heater temperatures above 600 K (327°C) have been reached in these tests. Two computer programs, RELAP5-3D and TRACE, have been used to simulate the tests. The TRACE code has shown good agreement with the test data and has been used to model a variety of tests. This experimental facility has been very useful in demonstrating the viability of this new type of irradiation facility.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Highly interactive, steered scientific workflows on HPC systems: optimizing design solutions
Scientific workflows are becoming increasingly important in high performance computing (HPC) settings, as the feasibility and appeal of many simultaneous heterogeneous tasks increases with increasing hardware capabilities. Currently no HPC-based workflow platform supports a dynamically adaptable workflow with interactive steering and analysis at run-time. Furthermore, for most workflow programs, compute resources are fixed for a given instance, resulting in a possible waste of expensive allocation resources when tasks are spawned and killed. Here we describe the design and testing of a run-time-interactive, adaptable, steered workflow tool capable of executing thousands of parallel tasks without an MPI programming model, using a database management system to facilitate task management through multiple live connections. We find that on the Oak Ridge Leadership Computing Facility pre-exascale Summit supercomputer it is possible to launch and interactively steer workflows with thousands of simultaneous tasks with negligible latency. For the case of particle simulation and analysis tasks that run for minutes to hours, this paradigm offers the prospect of a robust and efficient means to perform simulation-space exploration with on-the-fly analysis and adaptation.
Ossyra, John R↗
A Comparison of Select USDOE and USEPA Risk and Dose Assessment Tools Utilized at Two FUSRAP Projects - 17388
The US Army Corps of Engineers (USACE) utilizes several computer programs and models to conduct risk and dose assessments at Formerly Utilized Sites Remedial Action Program (FUSRAP) projects. In the past, the USACE has primarily relied on the US Department of Energy (USDOE)'s Argonne National Lab (ANL) RESidual Radioactivity (RESRAD) computer code for performing dose and risk assessments involving radiological contamination in onsite soils and the US Environmental Protection Agency (USEPA) Preliminary Remediation Goal (PRG) and Dose Compliance Calculator (DCC) for radiological contamination in ground water. Recent USEPA guidance requires use of USEPA calculators for all media at National Priority List sites, which includes some FUSRAP sites. The USACE recently utilized both RESRAD and the USEPA PRG and DCC calculators in assessments at two FUSRAP projects for all media in order to evaluate differences in the models and determine adjustments to potentially resolve differences. The two projects involve soil and groundwater assessments and address typical FUSRAP contaminants such as uranium, radium, thorium and their progeny. This paper discusses differences between the models, lessons learned, and methods used to address stakeholder and other concerns regarding these assessments. The paper's discussions are useful for anyone performing risk or dose assessments for radiological contamination and especially those using RESRAD and the PRG and DCC calculators. (authors)
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Method and system for semantic metadata compression in a two-tier storage system using copy-on-write
Example embodiments of the present invention relate to a method, a system, and a computer program product for managing metadata. The method includes determining a common characteristic of metadata entries for a plurality of files in a directory of a file system and computing a representative metadata entry for the plurality of metadata entries according to the common characteristic of the metadata entries. The plurality of metadata entries in the file system then may be represented according to the computed representative metadata entry.
97 MATHEMATICS AND COMPUTING↗
On-the fly scheduling of execution of dynamic hardware behaviors
Methods for dynamically executing computer code across multiple disparate processing unit architectures are disclosed. During execution of a first portion of computer code on a first processing unit, it is determined that a first dynamic hardware behavior of a plurality of dynamic hardware behaviors will occur at a subsequent point in time, based on a second dynamic hardware behavior that is occurring. The methods include determining to execute code corresponding to the first dynamic hardware behavior on a second processing unit, rather than the first processing unit, and scheduling computer program code corresponding to the first dynamic hardware behavior to execute on the second processing unit rather than the first processing unit. Upon completion of execution of the computer code corresponding to the first dynamic hardware behavior, a remaining portion of the computer code is scheduled to execute on the first processing unit.
97 MATHEMATICS AND COMPUTING↗
Atomic Gaussian type orbitals and their Fourier transforms via the Rayleigh expansion
Gaussian type orbitals (GTOs), which are one of the types of exponential type orbitals (ETOs), are used usually as basis functions in the multi-center atomic and molecular integrals to better understand physical and chemical properties of matter. In the Fourier transform method (FTM), basis functions have not simplicity to make mathematical operations, but their Fourier transforms are easier to use. In this work, with the help of FTM, Rayleigh expansion and some properties of unnormalized GTOs, we present new mathematical results for the Fourier transform of GTOs in terms of Laguerre polynomials, hypergeometric and Whittaker functions. Physical and analytical properties of GTOs are discussed and some numerical results have been given in a table. Finally, we compare our mathematical results with the other known literature results by using a computer program and details of evaluation are presented.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
A Comparison of GoldSim, RESRAD-Onsite, and NORMALYSA Modeling Tools at Radioactive Waste Sites - 19480
Radioactive waste disposal and associated decision making is a challenging problem around the world. A strong performance assessment model plays a crucial role in making good management and remediation decisions and in optimizing the performance of a site. The International Atomic Energy Agency (IAEA) MODARIA II Working Group 1 (WG1) is comparing the performance of various modeling platforms used at Naturally Occurring Radioactive Material (NORM) and legacy waste sites to support risk-informed decision making. Neptune and Company, Inc. (Neptune) has developed probabilistic performance assessment models using the GoldSim system modeling software at a variety of radioactive waste sites. This paper presents a comparison between GoldSim and two other tools used to evaluate radiation exposure and their associated risks: the NORM And Legacy Site Assessment (NORMALYSA) tool used to assess risk at the Zapadnoe mill tailings site, and the residual radioactivity (RESRAD-Onsite) computer program used to calculate Authorized Limits at the Los Alamos National Laboratory (LANL). The International Atomic Energy Agency (IAEA) supported development of NORMALYSA based off of Ecolego 6.0, in which the user connects several universal modules, representing various environmental transport and exposure pathways, and enters in site-specific data. Further work with NORMALYSA has continued under the MODARIA II WG 1. A radionuclide transport and dose calculation model of the Zapadnoe uranium mill tailings site at the Pridneprovsky Chemical Plant in the Ukraine has been configured by MODARIA II WG1 to test the capabilities of NORMALYSA. Under this test case, contaminant transport and effective dose calculations were performed for two exposure scenarios. The onsite worker is exposed via external ground irradiation and radon and radioactive dust inhalation. The offsite resident defines a residence established 300 years after institutional control has ended, 800 meters downstream from the site with exposure pathways including: ingestion of garden vegetables irrigated by contaminated ground water, ingestion of fish from the Konoplyanka River and agricultural activities in the contaminated garden. GoldSim modeling software is a dynamic probabilistic modeling platform that evaluates uncertainty in a robust and transparent manner. Although GoldSim provides contaminant and radionuclide transport modules to facilitate integrated modeling of radionuclide migration in the environment, these features are not embedded in an existing modeling framework, except as created by the user. Neptune developed a GoldSim model for the two exposure scenarios of the Zapadnoe test case. Neptune's Generic PA Model was used as a starting point for the Zapadnoe site. Input parameters for the GoldSim model include site-specific parameters from the MODARIA Working Group as well as parameters from the Generic PA Model and other PA models created by Neptune. For the other comparison, Version 7.0 of RESRAD-Onsite, developed by Argonne National Laboratory for the U.S. Department of Energy (DOE), was used to calculate Authorized Limits for LANL. Authorized Limits are concentrations of radionuclides in soil that support the evaluation process to clear land for release under different uses. The RESRAD-Onsite computer code was used to calculate Authorized Limits because this program was specifically designed to derive radionuclide soil guidelines corresponding to a specific dose criterion. The scope of the RESRAD dose or risk evaluation is constrained to the transport and exposure pathways that have been included in the computer code. To assist in illustrating the capability of GoldSim modeling software, a GoldSim model was created to replicate the calculation of LANL's Authorized Limits for soil. The main exposure pathways in this model are inhalation of dust and gas, plant ingestion, soil ingestion, and external ground irradiation. All of the processes and outputs of RESRAD-Onsite and NORMALYSA can be created in a GoldSim model, but the value of GoldSim lies in its flexibility for application to any of a wide variety of environmental modeling problems. GoldSim allows development of site-specific inventory, transport mechanisms and risk assessment scenarios. If all potentially complete release, transport, and exposure pathways are addressed in an existing dose and risk assessment computer code, and if the code's capabilities fulfill the project modeling objectives, then expending resources to create a GoldSim model may not be necessary. Site optimization and good decision making about site management are often best when supported with the flexibility of GoldSim as a modeling tool. (authors)
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
The case for free and open source software in research and scholarship
Free and open source software (FOSS) is any computer program released under a licence that grants users rights to run the program for any purpose, to study it, to modify it, and to redistribute it in original or modified form. Here, our aim is to explore the intersection between FOSS and computational reproducibility. We begin by situating FOSS in relation to other ‘open’ initiatives, and specifically open science, open research, and open scholarship. In this context, we argue that anyone who actively contributes to the research process today is a computational researcher, in that they use computers to manage and store information. We then provide a primer to FOSS suitable for anyone concerned with research quality and sustainability—including researchers in any field, as well as support staff, administrators, publishers, funders, and so on. Next, we illustrate how the notions introduced in the primer apply to resources for scientific computing, with reference to the GNU Scientific Library as a case study. We conclude by discussing why the common interpretation of ‘open source’ as ‘open code’ is misplaced, and we use this example to articulate the role of FOSS in research and scholarship today.
97 MATHEMATICS AND COMPUTING↗
Overview of the Tusas Code for Simulation of Dendritic Solidification
The aim of this project is to conduct a parametric investigation into the modeling of two dimensional dendrite solidification, using the phase field model. Specifically, we use the Tusas code, which is for coupled heat and phase-field simulation of dendritic solidification. Dendritic solidification, which may occur in the presence of an unstable solidification interface, results in treelike microstructures that often grow perpendicular to the rest of the growth front. The interface may become unstable if the enthalpy of the solid material is less than that of the liquid material, or if the solute is less soluble in solid than it is in liquid, potentially causing a partition [1]. A key motivation behind this research is that a broadened understanding of phase-field formulation and microstructural developments can be utilized for macroscopic simulations of phase change. This may be directly implemented as a part of the Telluride project at Los Alamos National Laboratory (LANL), through which a computational additive manufacturing simulation tool is being developed, ultimately to become part of the Advanced Simulation and Computing Program within the U.S. Department of Energy [2].
97 MATHEMATICS AND COMPUTING↗
Fundamental Optical Functions and Fine-Structure Parameters of Strontium–Selenide Crystal Transition Bands in the 0–40 eV Range
The spectral complex of fundamental optical functions of strontium-selenide (SrSe) crystal in the range 0–40 eV was determined. Spectra of the imaginary parts of three functions [transverse dielectric permittivity ε{sub 2}(E) and longitudinal volume –Im ε{sup –1} and surface –Im (1 + ε){sup –1} characteristic electron energy losses] of SrSe crystal were deconvoluted into 36 elementary bands. Energies of maxima E{sub i}, half-widths H{sub i}, amplitudes I{sub i}, areas S{sub i}, and oscillator strengths f{sub i} were determined for every band of the three spectra. The main spectral features of the optical functions and 36 elementary bands of SrSe in the range 3–35 eV due to excitons and transverse and longitudinal interband transitions were established. The calculations were based on experimental reflectance spectra and used the Kramers–Kronig integral relations, a modified combined Argand diagram method, and a number of computer programs.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Benchmark problems in aerosol evolution: Comparison of some exact and DSMC results
Aerosols are generated in many normal or accident situations associated with the nuclear enterprise. For a good understanding and modeling of the nuclear source term, for example, good experimental data and computational programs relating to aerosol evolution are needed. In the past several years there has been an effort to explore use of the Direct Simulation Monte Carlo (DSMC) approach for such estimations to improve fidelity of computations to the actual physics and chemistry of the accidents. An integral part of these efforts has been verification and validation of the DSMC technique against other available results wherever possible. Herein this paper explores verification of DSMC against one existing and two new benchmark problems covering condensation, coagulation, deposition and two-component aerosols. The simulations compare well with the exact results, providing further confidence in the use of DSMC.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗