Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer 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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Benchmarking the New RESRAD-OFFSITE Source Term Model with DUST-MS and GoldSim - 13377

RESRAD-OFFSITE is a computer code developed by Argonne National Laboratory under the sponsorship of U.S. Department of Energy (DOE) and U.S. Nuclear Regulatory Commission (NRC). It is designed on the basis of RESRAD (onsite) code, a computer code designated by DOE and NRC for evaluating soil-contaminated sites for compliance with human health protection requirements pertaining to license termination or environmental remediation. RESRAD-OFFSITE has enhanced capabilities of modeling radionuclide transport to offsite locations and calculating potential radiation exposure to offsite receptors. Recently, a new source term model was incorporated into RESRAD-OFFSITE to enhance its capability further. This new source term model allows simulation of radionuclide releases from different waste forms, in addition to the soil sources originally considered in RESRAD (onsite) and RESRAD-OFFSITE codes. With this new source term model, a variety of applications can be achieved by using RESRAD-OFFSITE, including but not limited to, assessing the performance of radioactive waste disposal facilities. This paper presents the comparison of radionuclide release rates calculated by the new source term model of RESRAD-OFFSITE versus those calculated by DUST-MS and GoldSim, respectively. The focus of comparison is on the release rates of radionuclides from the bottom of the contaminated zone that was assumed to contain radioactive source materials buried in soil. The transport of released contaminants outside of the primary contaminated zone is beyond the scope of this paper. Overall, the agreement between the RESRAD-OFFSITE results and the DUST-MS and GoldSim results is fairly good, with all three codes predicting identical or similar radionuclide release profiles over time. Numerical dispersion in the DUST-MS and GoldSim results was identified as potentially contributing to the disagreement in the release rates. In general, greater discrepancy in the release rates was found for short-lived, fast-moving radionuclides than for long-lived, slow-moving radionuclides. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

SENSMG: First-Order Sensitivities of Neutron Reaction Rates, Reaction-Rate Ratios, Leakage, k eff , and α Using PARTISN

SENSMG is a tool for calculating the first-order sensitivities of reaction-rate ratios, k eff , and α in critical problems and reaction-rate ratios, reaction rates, and leakage in fixed-source problems to multigroup cross sections, isotope densities, material mass densities, and interface locations using the PARTISN multigroup discrete-ordinates code by implementing Generalized Perturbation Theory. SENSMG can be used for one-dimensional spherical and slab (r) and two-dimensional cylindrical (r-z) geometries. For fixed-source (leakage) problems, SENSMG relies on the MISC and/or SOURCES4C codes to compute neutron source rate densities from spontaneous fission and (α,n) sources. SENSMG is a combination of Python and Fortran and was developed under Linux. This computer code abstract describes all user inputs, the input file, and output files. Furthermore, this computer code abstract describes how SENSMG can be modified to support different computer platforms, PARTISN versions, or cross-section availability. Several verification problems are presented in which SENSMG results are compared with MCNP6, SCALE6.2, and direct perturbations (central differences).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

phq: A Fortran code to compute phonon quasiparticle properties and dispersions

We here introduce a Fortran code that computes anharmonic free energy of solids from first-principles based on our phonon quasiparticle approach. In this code, phonon quasiparticle properties, i.e., renormalized phonon frequencies and lifetimes, are extracted from mode-projected velocity auto-correlation functions (VAF) of modes sampled by molecular dynamics (MD) trajectories. Using renormalized frequencies as input, the code next constructs an effective harmonic force constant matrix to calculate anharmonic phonon dispersions over the whole Brillouin zone and thus the anharmonic free energy in the thermodynamic limit ( N → ∞ ). A detailed description of all the input parameters and the subroutines is provided as well. We illustrate the use of this code to compute ab initio temperature-dependent anharmonic phonons of Si in the diamond structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Code Verification and Solution Verification framework in pin-resolved neutron transport code MPACT

Program verification in scientific computing encompasses the application of formal and mathematical techniques to a scientific computing code for its credibility, accuracy, and validity. Code Verification identifies bugs and performance issues in the software development stage. Solution Verification assesses the applicability of the code and the accuracy of the solution to problems of interest. Both activities utilize application cases and quantify the error against prescribed acceptance criteria. However, simply executing more application cases does not guarantee stronger or more comprehensive credibility. Here, we establish a verification framework that involves Code Verification and Solution Verification, both of which work together such that the overarching goal of “converge to the correct answer for the intended application” can be reasonably inferred. The application of such a verification framework is demonstrated using the pin-resolved neutron transport code MPACT, where standard unit tests and regression tests are covered, and where the Method of Exact Solutions and the Method of Manufactured Solutions are successfully used. Additionally, the applicability of Method of Manufactured Solutions is extended to the OECD/NEA C5G7 benchmark problems of practical material and geometric configurations. Solution Verification activities are demonstrated on a practical hierarchy of application models of increasing complexity ranging from 2D pin cell problems to 3D assembly problems. The convergence behavior and rate of convergence with respect to each individual variable are studied and provided. This framework can be adapted broadly to other fields involving scientific computing codes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Securely Aggregated Coded Matrix Inversion

Coded computing is a method for mitigating straggling workers in a centralized computing network, by using erasure-coding techniques. Federated learning is a decentralized model for training data distributed across client devices. In this work we propose approximating the inverse of an aggregated data matrix, where the data is generated by clients; similar to the federated learning paradigm, while also being resilient to stragglers. To do so, we propose a coded computing method based on gradient coding. We modify this method so that the coordinator does not access the local data at any point; while the clients access the aggregated matrix in order to complete their tasks. Here, the network we consider is not centrally administrated, and the communications which take place are secure against potential eavesdroppers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Blast-Wave Characteristics Near Site 300

We have developed a procedure for calculating the maximum high-explosive weight limit for use at Site 300 based on the analysis of the blast-wave propagation behavior in the atmosphere. Both analytical and experimental efforts have been involved in the present work, including a literature search, evaluation of other data, measurements of local weather conditions, and implementation of relevant computational programs. Major results of the present effort are • Comparison of the LLNL data with other data reveals discrepancies, suggesting that local effects, such as weather conditions in the downrange regions, may play a role in affecting the wave-propagation behavior. Our measurements of the local weather conditions indeed show that they differ from those at Site 300. • Analyses of other, newer models and data have resulted in the implementation of relevant computer codes, such as the BLASTO code, to be used at Site 300. • A new procedure for predicting the weight limit at Site 300 has been recommended. This new procedure consists of - Site 300 and local valley weather measurements. - The LAP code. - The BLASTO code. The LAP code is a modified SIPS code and computes the acoustic ray trajectories, including the focus ("caustic") locations, under varying local weather conditions. The BLASTO code calculates the peak overpressure vs distance relationship for a given charge weight and weather condition. Both the LAP code and the BLASTO code have been implemented and are operational at Site 300. Future plans include conducting selective tests to verify the validity of the predictive procedure recommended for Site 300.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

3D Coded SUMMA: Communication-Efficient and Robust Parallel Matrix Multiplication

In this paper, we propose a novel fault-tolerant parallel matrix multiplication algorithm called 3D Coded SUMMA that achieves higher failure-tolerance than replication-based schemes for the same amount of redundancy. This work bridges the gap between recent developments in coded computing and fault-tolerance in high-performance computing (HPC). The core idea of coded computing is the same as algorithm-based fault-tolerance (ABFT), which is weaving redundancy in the computation using error-correcting codes. In particular, we show that MatDot codes, an innovative code construction for parallel matrix multiplications, can be integrated into three-dimensional SUMMA (Scalable Universal Matrix Multiplication Algorithm [30]) in a communication-avoiding manner. To tolerate any two node failures, the proposed 3D Coded SUMMA requires ~50% less redundancy than replication, while the overhead in execution time is only about 5–10%.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

SCALE 6.3 Validation: Radiation Shielding

Safe and reliable use of scientific and engineering computer codes requires validation for the types of applications in which they will be used. An example in the nuclear reactor engineering and licensing field is radiation transport employed in shielding analyses. The validity of computer codes for shielding applications is demonstrated in this report for SCALE version 6.3.0. Representative benchmarks corresponding to shielding analyses are selected for the validation study. Typical measurement results analyzed from these benchmarks include neutron fluxes, detector count rates, detector energy response functions, neutron and gamma dose rates, neutron activation rates and activities, neutron leakage fluxes, and skyshine dose rates. Thousands of points of comparison between measurement and calculation are presented in this work. Other than rare outliers typically explained by either a lack of information or large uncertainties in the experiment conditions, material, or dimensions, the Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) radiation transport computer code with built-in variance reduction methods distributed with the SCALE computer code system agrees well with the measurement results. In selected benchmarks, MAVRIC is also compared to Monte Carlo N- Particle® (MCNP® ) 1 calculations. Both computer codes generally agree well within the estimated uncertainties. With the release of SCALE 6.3.0, Shift was integrated as an alternative transport solver in MAVRIC, denoted MAVRIC-Shift. Although the traditional MAVRIC using Monaco was used primarily in this validation study, many results have also been generated using MAVRIC-Shift. Agreement between MAVRIC-Monaco and MAVRIC-Shift is generally very good. The benchmarks presented in this report were obtained from reliable sources such as the International Criticality Safety Benchmark Evaluation Project Handbook, the Shielding Integral Benchmark Archive & Database, and other shielding validation work found in the literature. Additional datapoints and benchmarks will be added to future versions of this report to expand the shielding validation suite.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2017

Efficient convolutional sparse coding

Computationally efficient algorithms may be applied for fast dictionary learning solving the convolutional sparse coding problem in the Fourier domain. More specifically, efficient convolutional sparse coding may be derived within an alternating direction method of multipliers (ADMM) framework that utilizes fast Fourier transforms (FFT) to solve the main linear system in the frequency domain. Such algorithms may enable a significant reduction in computational cost over conventional approaches by implementing a linear solver for the most critical and computationally expensive component of the conventional iterative algorithm. The theoretical computational cost of the algorithm may be reduced from O(M.sup.3N) to O(MN log N), where N is the dimensionality of the data and M is the number of elements in the dictionary. This significant improvement in efficiency may greatly increase the range of problems that can practically be addressed via convolutional sparse representations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
DOE OSTI2013

Analysis Approach and Data Package for Mayak Public Doses: Data Supporting the Calculation of Representative Doses to Members of the Public from Airborne Releases of 131 I at the Mayak Production Association from 1948 through 1972

Historical activities at facilities producing nuclear materials for weapons released radioactivity into the air and water. Past studies in the United States have evaluated the release, atmospheric transport and environmental accumulation of 131 I from the nuclear facilities at Hanford in Washington State and the resulting dose to members of the public (Farris et al. 1994). A multi-year dose reconstruction effort (Mokrov et al. 2004) is also being conducted to produce representative dose estimates for members of the public living near Mayak, Russia, from atmospheric releases of 131 I at the facilities of the Mayak Production Association. The approach to calculating individual doses to members of the public from historical releases of airborne 131 I has the following general steps: • Construct estimates of releases 131 I to the air from production facilities. • Model the transport of 131 I in the air and subsequent deposition on the ground and vegetation. • Model the accumulation of 131 I in soil, water and food products (environmental media). • Calculate the dose for an individual by matching the appropriate lifestyle and consumption data for the individual to the concentrations of 131 I in environmental media at their residence location. A number of computer codes were developed to facilitate the study of airborne 131 I emissions at Hanford. Of particular interest is DESCARTES code that modeled accumulation of 131 I in environmental media (Miley et al. 1994). In addition, the CIDER computer code estimated annual doses to individuals (Eslinger et al. 1994) using the equations and parameters specific to Hanford (Snyder et al. 1994). Several of the computer codes developed to model 131 I releases from Hanford are general enough to be used for other facilities. Additional codes have been developed, including the new individual dose code CiderF (Eslinger and Napier 2013), and applied to historical releases of 131 I from Mayak. This document provides a data package that identifies computer code runs and associated input and output files prepared for the purpose of calculating doses to members of the public from atmospheric releases of 131 I at the Mayak Production Association for the time period 1948 through 1972.

Mayak↗
DOE OSTI2020

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↗
DOE OSTI2014

An MCNP6 Version of the Visual Editor - Paper 2

The Visual Editor is a graphical user interface for working with the MCNP computer code developed at Los Alamos National Laboratory. It consists of a graphical user interface written in C++ that is linked directory to the MCNP Fortran code. Work on the original Visual Editor started in 1992 and was first published in 1993. The original version of the Visual Editor was created on Unix/Linux operating systems and was linked to the MCNP computer code. The Visual Editor was adapted to the Windows operating system in 1999. In 2008 the Visual Editor was adapted to work with the MCNPX computer code. This paper discusses the latest upgrade of the Visual Editor to be compatible with the MCNP6 computer code. Two graphical user interfaces have been created, the Visual MCNP6 Plotter and the Visual MCNP6 Editor. The Visual MCNP6 Plotter includes all of the standard plotting features available in the default MCNP6 plotting package including geometry plotting, tally plotting, cross section plotting, mesh tally plotting and weight window mesh plotting. The Visual MNCP6 Editor has additional capabilities to allow for the creation of MCNP input files with graphical tools. The Visual MCNP6 Editor allows user to create and modify the input file within the Editor. Modifications are not allowed in the Visual MCNP6 Plotter. Both graphical user interfaces allow for advanced plotting features including 3D geometry plots and particle track plotting. Complete details on these tools along with their current status can be found at the Visual Editor website (www.mcnpvised.com). (author)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

A Comprehensive Validation Approach Using The RAVEN Code

The RAVEN computer code , developed at the Idaho National Laboratory, is a generic software framework to perform parametric and probabilistic analysis based on the response of complex system codes. RAVEN is a multi-purpose probabilistic and uncertainty quantification platform, capable to communicate with any system code. A natural extension of the RAVEN capabilities is the imple- mentation of an integrated validation methodology, involving several different metrics, that represent an evolution of the methods currently used in the field. The state-of-art vali- dation approaches use neither exploration of the input space through sampling strategies, nor a comprehensive variety of metrics needed to interpret the code responses, with respect experimental data. The RAVEN code allows to address both these lacks. In the following sections, the employed methodology, and its application to the newer developed thermal-hydraulic code RELAP-7, is reported.The validation approach has been applied on an integral effect experiment, representing natu- ral circulation, based on the activities performed by EG&G Idaho. Four different experiment configurations have been considered and nodalized.

PRA↗
DOE OSTI2017

Operator Action Time Evaluation of Emergency Water Injection Strategy for Extreme SGTR Accidents

In terms of severe accidents caused by earthquakes or tsunamis that are beyond expectation, a special safety inspection for operating plants, following the Fukushima accident, has been conducted by the government of Korea to verify that nation's nuclear power plants have been adequately designed to respond to extreme accidents. The inspection results show that Korean nuclear power plants in operation have no imminent risk for an expected maximum potential earthquake or coastal flooding based on the up-to-date investigation. However, there is a need to implement long- and short-term safety improvements in order to secure the level of safety for natural beyond-design-basis events. One of the measures to increase the mitigation capability during extreme accidents is installing the injection flow paths to provide emergency cooling water of extemal sources using fire engines to the steam generators (SGs) or reactor cooling system (RCS). An emergency cooling water system consists of a fixed pipe connected from the RCS or SGs to the outside of the containment. A standby valve is installed on the pipe. Following the occurrence of an extreme accident, movable equipment (for example, a fire truck hose) can be connected to the pipe hole at the opening of the isolation valve. The emergency cooling water injection strategy was not taken into account for the steam generator tube rupture (SGTR) accident management strategy of the OPR-1000. Therefore, it is necessary to develop some guidelines or strategies to cope with an extreme SGTR scenarios using the newly installed injection flow paths and fire engines. As long as the fission products released from the core during severe accidents are kept inside the containment, there may be several ways to reduce the radionuclide release into the environment, for example, by using plant safety features. The containment bypass sequences, however, allow the fission products to be released into the environment or directly to the auxiliary building. The SGTR accidents, which are one of the containment bypass sequences, have shown during various probabilistic safety assessments (PSAs) to be contributors to plant risk, mainly because of their potential for causing a release outside containment. This paper is intended to be an investigation of the characteristics of accident progressions and assessments for the mitigative effectiveness of the external cooling water injection strategies during extreme SGTR accidents in a typical OPR-1000 plant. Analyses were performed using a MAAP computer code version 5.03. The MAAP code is a system-level computer code capable of conducting integral analyses of potential severe accident progressions in nuclear power plants, whose main purpose is to support a level 2 PSA or severe accident management strategy developments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

User's Manual for RESRAD-BUILD Code Version 4: Vol. 2 - User's Guide for RESRAD-BUILD Code Version 4

The RESRAD-BUILD computer code is designed to assess radiological doses to individuals who live or work in a building contaminated with radioactive material. The code is equipped with a user-friendly interface that has many features to facilitate using the computer code and understanding the results. The design of the interface provides various options, from entering data and performing calculations to displaying calculation results. General and context-specific help are available, providing information on editing and viewing the radionuclide database, the definitions of input parameters and their use in the calculations, and the selection of the calculation results for placement into other applications. Two types of sensitivity analysis are supported by the code, i.e., deterministic and probabilistic, that can be used to study the influence of input parameters on the calculation results. This user’s guide provides instructions to help users on how to install the RESRAD-BUILD code, navigate the interface, and use the various features to set up a dose/risk analysis and to view/print the results in text and graphical outputs.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

An Overview of the State-of-the-Art Reactor Consequence Uncertainty Assessment Accident Progression Insights

The U.S. Nuclear Regulatory Commission (NRC) with Sandia National Laboratories (Sandia) have completed three uncertainty analyses (UAs) as part of the State-of-the-Art Reactor Consequence Analyses (SOARCA) program. The SOARCA UAs included an integrated evaluation of uncertainty in accident progression, radiological release, and offsite health consequence projections. The UA for Peach Bottom, a boiling-water reactor (BWR) with a Mark I containment located in the State of Pennsylvania, analyzed the unmitigated long-term station blackout SOARCA scenario. The UA for Sequoyah, a 4-loop Westinghouse pressurized-water reactor (PWR) located in the State of Tennessee, analyzed the unmitigated short-term station blackout SOARCA scenario, with a focus on issues unique to the ice condenser containment and the potential for early containment failure due to hydrogen deflagration. The UA for Surry, a 3-loop Westinghouse PWR with a sub-atmospheric large dry containment located in the State of Virginia, analyzed the unmitigated short-term station blackout SOARCA scenario including the potential for thermally-induced steam-generator tube rupture. These three UAs are currently documented in three NUREG/CR reports. This report provides input to planned NRC documentation on the insights and findings from the SOARCA UA program. The purpose of the summary report is to provide a useful reference for regulatory applications that require the evaluation of offsite consequence risk from beyond design basis event severe accidents. This report focuses on the accident progression and source term insights developed from the MELCOR analyses. MELCOR is the NRC's best-estimate, severe accident computer code used in the SOARCA UAs. In anticipation of the SOARCA UA insights work, NRC and Sandia benchmarked the response of the Peach Bottom model to selected reference calculations from the Peach Bottom SOARCA UA. Peach Bottom was the first SOARCA UA performed and was completed in 2015 using the MELCOR 1.8.6 code. The PWR SOARCA UAs evolved the original methodology and utilized the updated MELCOR 2.2 computer code. The Peach Bottom model has been systematically updated for other NRC research efforts and has been updated to MELCOR 2.2. computer code. The findings from the new reference calculations using the updated model with the MELCOR 2.2 code are also integrated into the report. A second objective is an assessment of the applicability of the results to the other nuclear reactors in the U.S. As the key findings are reviewed, judgments are presented on the applicability of the results to other U.S. nuclear power plants. An important objective of the SOARCA program relied on high- fidelity plant-specific modeling. However, the nature of the insights and conclusions allowed judgements to be made on the applicability of the various insights to the same general classification of plant (i.e., BWR or PWR) or the entire fleet of plants. Finally, the results from the SOARCA UA accident progression calculations contain a wealth of information not previously documented in the NUREG/CRs. This report includes new but related information that can be used to benchmark past or support future regulatory decisions related to severe accidents. The new work includes a benchmark of the NUREG-1465 licensing source term definitions, the variability of key accident progression events and timing to radionuclide release, and an improved understanding of the timing and source terms from consequential steam generator tube ruptures. iii ACKNOWLEDGEMENTS The Sandia authors gratefully acknowledge the significant technical and programmatic contributions from the NRC SOARCA team which are reflected throughout the report. Dr. Tina Ghosh has been involved throughout the SOARCA UAs, providing the primary managerial and technical oversight. The long lists of NRC and Sandia contributors from the SOARCA UAs are cited in the three NUREG/CRs and are also gratefully acknowledged by the small team of authors compiling the results of their efforts. Significant technical contributions, advice, and reviews were provided by Dr. Hossein Esmaili, Dr. Alfred Hathaway, and Dr. Edward Fuller (retired) of the NRC. Dr. Randal Gauntt (retired), Mr. Patrick Mattie, Mr. Joseph Jones (retired), and Dr. Doug Osborn from Sandia are recognized as the SOARCA UA managers guiding the past efforts. There is a comparable list of project managers at the NRC including Ms. Patricia Santiago, Dr. Salman Haq, and Mr. Jon Barr. Sadly, we have lost Mr. Charlie Tinkler and Mr. Robert Prato, who were important contributors to the original SOARCA project. Finally, Mr. Kyle Ross and Mr. Mark Leonard have also retired but were significant technical contributors. Mr. Kyle Ross was the technical lead on all three SOARCA UAs and the original pressurized water reactor SOARCA study. Mr. Leonard was the technical lead on the original boiling water reactor SOARCA study and a key contributor to the first Peach Bottom SOARCA UA. iv

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗