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 163 records · Page 9

NEAMS Workbench and Bison Fuel Performance Remote Application Configuration

The mission of the US Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program is to develop, apply, deploy, and support state-of-the-art predictive modeling and simulation tools for the design and analysis of current and future nuclear energy systems. This is accomplished using computing architectures that range from laptops to leadership class facilities. The NEAMS Workbench [1,2,3] initiative will facilitate the transition from conventional computer codes to high-fidelity ones by providing a common analysis environment and interface for model creation, review, execution, output review, and visualization for integrated codes. The BISON fuel performance code [4] is a finite element method (FEM)-based software being developed by multiple organizations led by the Idaho National Laboratory for the US Department of Energy Office of Nuclear Energy. BISON can simulate a variety of fuel and cladding compositions and geometries in varying dimensionalities, depending upon the problem type and complexity, including one-dimensional (1D), 1.5D, 2D axially symmetric (2D-rz), 2D radial-azimuthal (2D-rθ), or 3D models. The NEAMS Workbench allows execution of the BISON fuel performance code on local and remote compute resources. In this report, the detailed instructions of configuring the NEAMS Workbench for remote execution of the BISON fuel performance code on a portable batch system (PBS)-scheduled compute resource, as well the transfer of associated Bison files, is illustrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The SAS4A/SASSYS-1 Safety Analysis Code System, Version 5

The SAS4A/SASSYS-1 computer code is developed by Argonne National Laboratory for thermal, hydraulic, and neutronic analysis of power and flow transients in liquidmetal- cooled nuclear reactors (LMRs). SAS4A was developed to analyze severe core disruption accidents with coolant boiling and fuel melting and relocation, initiated by a very low probability coincidence of an accident precursor and failure of one or more safety systems. SASSYS-1, originally developed to address loss-of-decay-heat-removal accidents, has evolved into a tool for margin assessment in design basis accident (DBA) analysis and for consequence assessment in beyond-design-basis accident (BDBA) analysis. SAS4A contains detailed, mechanistic models of transient thermal, hydraulic, neutronic, and mechanical phenomena to describe the response of the reactor core, its coolant, fuel elements, and structural members to accident conditions. The core channel models in SAS4A provide the capability to analyze the initial phase of core disruptive accidents, through coolant heat-up and boiling, fuel element failure, and fuel melting and relocation. Originally developed to analyze oxide fuel clad with stainless steel, the models in SAS4A have been extended and specialized to metallic fuel with advanced alloy cladding. SASSYS-1 provides the capability to perform a detailed thermal/hydraulic simulation of the primary and secondary sodium coolant circuits and the balance-ofplant steam/water circuit. These sodium and steam circuit models include component models for heat exchangers, pumps, valves, turbines, and condensers, and thermal/hydraulic models of pipes and plena. SASSYS-1 also contains a plant protection and control system modeling capability, which provides digital representations of reactor, pump, and valve controllers and their response to input signal changes.

97 MATHEMATICS AND COMPUTING↗

Practical methods to improve the development of computational software

The use of computation has become ubiquitous in science and engineering. As the complexity of computer codes has increased, so has the need for robust methods to minimize errors. Past work has show that the number of functional errors is related the number of commands that a code executes. Since the late 1960's, major participants in the field of computation have encouraged the development of best practices for programming to help reduce coder induced error, and this has lead to the emergence of 'software engineering' as a field of study. Best practices for coding and software production have now evolved and become common in the development of commercial software. These same techniques, however, are largely absent from the development of computational codes by research groups. Many of the best practice techniques from the professional software community would be easy for research groups in nuclear science and engineering to adopt. This paper outlines the history of software engineering, as well as issues in modern scientific computation, and recommends practices that should be adopted by individual scientific programmers and university research groups. (authors)

97 MATHEMATICS AND COMPUTING↗

Applications of “Synchrotron Radiation Workshop” Code

The “Synchrotron Radiation Workshop” (SRW) computer code was started in the group of Insertion Devices of the European Synchrotron Radiation Facility in 1997. The idea behind it was to combine in one software library electrodynamics methods for calculating electric fields (and derived characteristics) of synchrotron radiation (SR) generated by relativistic electrons moving in magnetic fields of nearly arbitrary configuration with wave-optics methods for simulating propagation of radiation through optical elements of beamlines. Further, the use of this general approach was initially driven by applications in electron beam diagnostics using SR. However, with time it started to be used for more complicated beamlines of synchrotron light sources, and the use of the code for these applications continues to increase, with the new light sources offering brighter and more coherent x ray beams for wide ranges of experiments in different areas of science.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SURVEY OF COUPLING SCHEMES IN TRADITIONAL COUPLED NEUTRONICS AND THERMAL-HYDRAULICS CODES

The design and the analysis of nuclear power plants (NPPs) require computational codes to predict the behavior of the NPP nuclear components and other systems (i.e., reactor core, primary coolant system, emergency core cooling system, etc.). Coupled calculations are essential to the conduct of deterministic safety assessments. Inasmuch as the physical phenomena that govern the performance of a nuclear reactor are always present simultaneously, ideally computational modeling of a nuclear reactor should include coupled codes that represent all of the active physical phenomena. Such multi-physics codes are under development at several institutions and are expected to become operational in the future. However, in the interim, integrated codes that incorporate modeling capabilities for two to three physical phenomena will remain useful. For example, in the conduct of safety analyses, of paramount importance are codes that couple neutronics and thermal-hydraulics, especially transient codes. Other code systems of importance to safety analyses are those that couple primary system thermal-hydraulics to fission product chemistry, neutronics to fuel performance, containment behavior and structural mechanics to thermal-hydraulics, etc. This paper surveys the methods used traditionally in the coupling of neutronic and thermal-hydraulics codes. The neutron kinetics codes are used for computing the space-time evolution of the neutron flux and, hence, of the power distribution. The thermal-hydraulics codes, which compute mass, momentum and energy transfers, model the coolant flow and the temperature distribution. These codes can be used to compute the neutronic behavior and the thermal-hydraulic states separately. However, the need to account with fidelity for the dynamic feedback between the two sets of properties (via temperature and density effects on the cross section inputs into the neutronics codes) and the requirement to model realistically the transient response of nuclear power plants and to assess the associated emergency systems and procedures imply the necessity of modeling the neutronic and thermal-hydraulics simultaneously within a coupled code system. The focus of this paper is a comparison of the methods by which the coupling between neutron kinetics and thermal-hydraulics treatments has been traditionally achieved in various code systems. As discussed in the last section, the modern approaches to multi-physics code development are beyond the scope of this paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computational Physics Overview [Slides]

Computational physics is an important part of the overall investment in National Security Science at LANL. Computational physics is the study and implementation of numerical analysis to solve problems in physics for which a quantitative theory already exists. Historically, computational physics was the first application of modern computers in science. There are three key elements to computational physics: mathematical models of physical phenomena and conservation equations, computer codes that implement these models, and computer platforms that execute the code instructions and manipulate the data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Demonstrating GPU code portability and scalability for radiative heat transfer computations

High performance computing frameworks utilizing CPUs, Nvidia GPUs, and/or Intel Xeon Phis necessitate portable and scalable solutions for application developers. Nvidia GPUs in particular present numerous portability challenges with a different programming model, additional memory hierarchies, and partitioned execution units among streaming multiprocessors. This work presents modifications to the Uintah asynchronous many-task runtime and the Kokkos portability library to enable one single codebase for complex multiphysics applications to run across different architectures. Scalability and performance results are shown on multiple architectures for a globally coupled radiation heat transfer simulation, ranging from a single node to 16,384 Titan compute nodes.

Computer Science↗

The CCONE Code System and its Application to Nuclear Data Evaluation for Fission and Other Reactions

A computer code system, CCONE, was developed for nuclear data evaluation within the JENDL project. The CCONE code system integrates various nuclear reaction models needed to describe nucleon, light charged nuclei up to alpha-particle and photon induced reactions. The code is written in the C++ programming language using an object-oriented technology. At first, it was applied to neutron-induced reaction data on actinides, which were compiled into JENDL Actinide File 2008 and JENDL-4.0. It has been extensively used in various nuclear data evaluations for both actinide and non-actinide nuclei. The CCONE code has been upgraded to nuclear data evaluation at higher incident energies for neutron-, proton-, and photon-induced reactions. It was also used for estimating β-delayed neutron emission. This paper describes the CCONE code system indicating the concept and design of coding and inputs. Details of the formulation for modelings of the direct, pre-equilibrium and compound reactions are presented. Applications to the nuclear data evaluations such as neutron-induced reactions on actinides and medium-heavy nuclei, high-energy nucleon-induced reactions, photonuclear reaction and β-delayed neutron emission are mentioned.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Validation of MCNP6 Version 1.0 with the ENDF/B-VII.1 Cross Section Library for Plutonium Metals, Oxides, and Solutions on the High Performance Computing Platform Moonlight

This report documents a validation of the MCNP6 Version 1.0 computer code on the high performance computing platform Moonlight, for operations at Los Alamos National Laboratory (LANL) that involve plutonium metals, oxides, and solutions. The validation is conducted using the ENDF/B-VII.1 continuous energy group cross section library at room temperature. The results are for use by nuclear criticality safety personnel in performing analysis and evaluation of various facility activities involving plutonium materials.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Replacement of Legacy Analytical Codes at the Advanced Test Reactor

For each operating cycle of the Advanced Test Reactor (ATR) at Idaho National Laboratory, a Core Safety Assurance Package (CSAP) is necessary to demonstrate compliance with the safety basis approved by the United State Department of Energy (DOE). Certain computer codes are used in CSAP development, most of them developed in-house. This work describes replacement of a large set of these codes and updates previous work. Replacement of legacy codes is necessary due to computer hardware failure but also has generally improved user-friendliness.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Replacement of Legacy Analytical Codes at the Advanced Test Reactor

For each operating cycle of the Advanced Test Reactor (ATR) at Idaho National Laboratory, a Core Safety Assurance Package (CSAP) is necessary to demonstrate compliance with the safety basis approved by the United State Department of Energy (DOE). Certain computer codes are used in CSAP development, most of them developed in-house. This work describes replacement of a large set of these codes and updates previous work. Replacement of legacy codes is necessary due to computer hardware failure but also has generally improved user-friendliness.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A 2D/1D coupling neutron transport method based on the matrix MOC and NEM methods

A new 2D/1D coupling method based on the matrix MOC method (MMOC) and nodal expansion method (NEM) is proposed for solving the three-dimensional heterogeneous neutron transport problem. The MMOC method, used for radial two-dimensional calculation, constructs a response matrix between source and flux with only one sweep and then solves the linear system by using the restarted GMRES algorithm instead of the traditional trajectory sweeping process during within-group iteration for angular flux update. Long characteristics are generated by using the customization of commercial software AutoCAD. A one-dimensional diffusion calculation is carried out in the axial direction by employing the NEM method. The 2D and ID solutions are coupled through the transverse leakage items. The 3D CMFD method is used to ensure the global neutron balance and adjust the different convergence properties of the radial and axial solvers. A computational code is developed based on these theories. Two benchmarks are calculated to verify the coupling method and the code. It is observed that the corresponding numerical results agree well with references, which indicates that the new method is capable of solving the 3D heterogeneous neutron transport problem directly. (authors)

97 MATHEMATICS AND COMPUTING↗

Similarity-transformed equation-of-motion vibrational coupled-cluster theory

A similarity-transformed equation-of-motion vibrational coupled-cluster (STEOM-XVCC) method is introduced as a one-mode theory with an effective vibrational Hamiltonian, which is similarity transformed twice so that its lower-order operators are dressed with higher-order anharmonic effects. The first transformation uses an exponential excitation operator, defining the equation-of-motion vibrational coupled-cluster (EOM-XVCC) method, and the second uses an exponential excitation-deexcitation operator. From diagonalization of this doubly similarity-transformed Hamiltonian in the small one-mode excitation space, the method simultaneously computes accurate anharmonic vibrational frequencies of all fundamentals, which have unique significance in vibrational analyses. We establish a diagrammatic method of deriving the working equations of STEOM-XVCC and prove their connectedness and thus size-consistency as well as the exact equality of its frequencies with the corresponding roots of EOM-XVCC. We furthermore elucidate the similarities and differences between electronic and vibrational STEOM methods and between STEOM-XVCC and vibrational many-body Green’s function theory based on the Dyson equation, which is also an anharmonic one-mode theory. The latter comparison inspires three approximate STEOM-XVCC methods utilizing the common approximations made in the Dyson equation: the diagonal approximation, a perturbative expansion of the Dyson self-energy, and the frequency-independent approximation. The STEOM-XVCC method including up to the simultaneous four-mode excitation operator in a quartic force field and its three approximate variants are formulated and implemented in computer codes with the aid of computer algebra, and they are applied to small test cases with varied degrees of anharmonicit

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computing the cross sections of nuclear reactions with nuclear clusters emission for proton energies between 30 MeV and 2.6 GeV

The cross sections of nuclear reactions involving emission of clusters of light nuclei in proton collisions with a heavy-metal target are computed for incident-proton energies between 30 MeV and 2.6 GeV. The calculation relies on the ALICE/ASH and CASCADE/INPE computer codes. The parameters determining the pre-equilibrium cluster emission are varied in the computation.

74 ATOMIC AND MOLECULAR PHYSICS↗

Monte Carlo simulation of photon buildup factors for shielding materials in radiotherapy x-ray facilities

Purpose: This paper presents the results of a series of calculations to determine buildup factors for ordinary concrete, baryte concrete, lead, steel, and iron in broad beam geometry for photons energies from 0.125 to 25.125 MeV at 0.250 MeV intervals.Methods: Monte Carlo N-particle radiation transport computer code has been used to determine the buildup factors for the studied shielding materials.Results: The computation of the primary broad beams using buildup factors data was done for nine published megavoltage photon beam spectra ranging from 4 to 25 MV in nominal energies, representing linacs made by the three major manufacturers. The first tenth value layer and the equilibrium tenth value layer are calculated from the broad beam transmission for these nine primary megavoltage photon beam spectra.Conclusions: The results, compared with published data, show the ability of these buildup factor data to predict shielding transmission curves for the primary radiation beam. Therefore, the buildup factor data can be combined with primary, scatter, and leakage x-ray spectra to perform computation of broad beam transmission for barriers in radiotherapy shielding x-ray facilities.

61 RADIATION PROTECTION AND DOSIMETRY↗

MCNP ® Code Version 6.3.1 Theory & User Manual

This document acts as a repository of knowledge for the Monte Carlo N-Particle (MCNP) transport computer code. It is maintained alongside the source code and attempts to introduce new users and re-familiarize experienced users with the theory and practices of using the MCNP code for the wide range of particle transport analyses that it is appropriate for. The latest version of the MCNP code, version 6.3.1, provides the Monte Carlo particle transport community with the latest feature developments and bug fixes in the MCNP code. The MCNP code version 6.0 and later is also known as the MCNP6 code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ecological risk assessment software for radioactive waste management facilities - 15464

Environmental Risk Assessment (ERA) is a systematic process which can be used to identify, quantify and characterize the risk posed by contaminants and physical stressors in the environment on biological receptors, including human and non-human biota. Human receptors can be addressed through a Human Health Risk Assessment and non-human biota can be addressed through an Ecological Risk Assessment (EcoRA). In 2012 the Canadian Standards Association (CSA) issued a CSA standard N288.6 Environmental Risk Assessment at Class I Nuclear Facilities and Uranium Mines and Mills. The standard applies to radioactive waste management facilities, as well as nuclear power plants. The ERA can also be used retrospectively, to estimate the effect that a contaminant or stressor has already had on the environment; and can be used to assist in determining possible remedial or risk management options to return the quality of the affected environment to a predetermined state. While industry-standard tools for estimating radiological doses to humans are readily available, there is a need for a tool to implement the CSA N288.6 standard for Ecological Risk Assessment. To address this need, AMEC developed a computer code called -Assessment of Impact of Contaminants on Ecological Receptors (AICER)-. The code was developed to be consistent with the CSA standard N288.6, by assessing the impact of radiological contamination on non-human biota for normal operations. The code can also be used to screen for the potential effects of conventional contamination on the environment. In addition, the code is also applicable for abnormal events and accidents with the appropriate modification of some parameters used in the code. The AICER code was qualified in accordance with the requirements of the Canadian quality standard CSA N286.7 -Quality Assurance of Analytical, Scientific, and Design Computer Programs for Nuclear Power Plants-. The functionality of the code is illustrated with screen shots in the paper, following the introduction of the theoretical model for the dose calculation. An inter-comparison with RESRAD-BIOTA and ERICA Assessment tool at the conceptual level was conducted. In addition, the practical application of the code for the assessment of the impacts related to an existing radioactive waste management facility is provided. In particular, analysis was conducted to evaluate radiological doses to non-human biota resulting from a fire accident scenario involving an existing waste storage and processing facility. Analysis included all feasible waterborne and airborne exposure pathways. Doses were evaluated for a range of species representing valued ecological components at various trophic levels and compared with benchmark values for indicator species. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Coupled Neutron and Gamma Heating Calculation Based on VARIANT Nodal Transport Method

Accurate determination of temperature distributions in fuel pins and duct walls is essential for fuel-pin performance and assembly bowing analyses. Currently, the multi-assembly, sub-channel analysis code SE2-ANL is used in thermalhydraulic analysis of sodium cooled fast reactors (SFR). SE2- ANL is interfaced with the DIF3D/GAMSOR/DIF3D procedure for coupled neutron and gamma heating calculations, developed at Argonne National Laboratory (ANL). In this procedure, the neutron flux distribution is first calculated using the finite-difference diffusion theory option of the DIF3D code. Using the neutron flux solution, the GAMSOR code computes the gamma source distribution for the subsequent gamma flux calculation with DIF3D. However, in the current heating calculation scheme, the neutron and gamma flux calculations are limited to six triangular meshes per hexagonal assembly. In addition, the pin power distributions are estimated by assuming that the power distribution within an assembly is separable in the radial and axial directions. In each assembly, the axial power profile is obtained by averaging the power densities of six triangular meshes at each axial mesh interval. The radial power profile is evaluated at a single axial mesh (e.g., at the core mid-plane) by approximating the neutron and gamma flux distributions by a linear or quadratic polynomial of two variables x and y, depending on the polynomial order specified by the user. The coefficients of the linear and quadratic polynomials are determined by the least squares fit of six and 24 triangular mesh fluxes, respectively. Furthermore, the total heat generation rate in an assembly duct wall is simply determined by a user-specified fraction of the total gamma heating within that assembly. With the advancement of computer technology, the nodal transport theory code VARIANT is routinely used in core design analyses of SFR instead of DIF3D. It is therefore necessary to develop a new heating calculation procedure based on VARIANT neutron and gamma transport calculations. In this study, a new, coupled neutron and gamma heating calculation scheme has been developed to determine the pin power distributions more accurately by eliminating the aforementioned limitations of the current procedure. For this, the GAMSOR code has been extended to generate the intra-nodal gamma source distribution to be consistent with the trial functions of VARIANT. In the new heating calculation procedure, the neutron flux distribution is determined first by solving an eigenvalue transport problem using VARIANT. The obtained neutron flux is then used in GAMSOR to calculate the gamma source distribution. Using this gamma source distribution, the gamma flux distribution is determined by solving a fixed-source transport problem using VARIANT again. Using the computed neutron and gamma flux distributions, the power distribution in each fuel pin and duct wall is determined using a new utility code CURVE, which has been developed to evaluate the VARIANT flux solution at each fuel pin position and along the duct mid-wall. For each pin segment contained in a VARIANT node, the axial power shape is represented as a quadratic profile. The performance of the new heating calculation scheme based on VARIANT transport solutions has been tested against the existing method of SE2-ANL, using the ABR-1000 metal core design. The gamma source distribution and the assembly, fuel pin and duct wall power distributions were calculated and compared. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗