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 1,567 records · Page 87

PHISICS multi-group transport neutronic capabilities for RELAP5

PHISICS is a neutronic code system currently under development at INL. Its goal is to provide state of the art simulation capability to reactor designers. This paper reports on the effort of coupling this package to the thermal hydraulic system code RELAP5. This will enable full prismatic core and system modeling and the possibility to model coupled (thermal-hydraulics and neutronics) problems with more options for 3D neutron kinetics, compared to the existing diffusion theory neutron kinetics module in RELAP5 (NESTLE). The paper describes the capabilities of the coupling and illustrates them with a set of sample problems. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DEVELOPMENT OF A SIMULATION CODE FOR A COOL-DOWN PROCESS OF THE CRYOGENIC HYDROGEN SYSTEM

Supercritical hydrogen with a pressure of 1.5 MPa and a temperature of 20 K has been selected as a moderator material in an intense spallation neutron source (JSNS), which is one of main experimental facilities in J-PARC. The cryogenic hydrogen system, in which a hydrogen circulation system is cooled by a helium refrigerator with the refrigeration power of 6.45 kW at 15.5 K, has been designed to provide the supercritical hydrogen to the moderator and to remove the nuclear heating generated there. In this study, we have developed a simulation code that predicts temperature behaviors in the hydrogen circulation system during its cool-down process. Cool-down process analyses have been performed, and an operational method for the cool-down process has been studied. The analytical results indicate that the hydrogen circulation system would be able to be cooled down to 18 K within 19 hours.

42 ENGINEERING↗

Comparison of nuclear mass models at nuclearmasses.org

Nuclear masses are crucial in many astrophysics studies, as well as other areas of basic and applied nuclear science. There are now many different global theoretical models of nuclear masses - but choosing the best model for a given application can be challenging. The suite of codes online at nuclearmasses.org greatly facilitates the comparison of nuclear models with evaluated masses as well as the intercomparison of theoretical models. We demonstrate the utility of this suite by comparing the RMS deviations of 13 different theoretical mass models from the AME2003 evaluated masses over mass ranges appropriate for astrophysics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Validation and bias quantification of criticality safety codes for SRNS operations

The validation and bias quantification methods employed by the Criticality Safety Engineering group at Savannah River Nuclear Solutions (SRNS) are discussed. Both SCALE and MCNP are validated against benchmark experiments, where these benchmark models have been divided into subgroups representative of the types of fissile material operations at SRNS. A quantification of the bias is made for each subgroup for each code and the three methods employed for quantification are discussed. The bias information is employed directly in the establishment of the Upper Safety Limit (USL) value of the neutron multiplication factor when using standardized documented SRNS methods. A simple example of validation and bias quantification is given and the procedure for setting the USL using the derived bias information is also presented. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Additional Research Needs to Support the GENII Biosphere Models

In the course of evaluating the current parameter needs for the GENII Version 2 code (Snyder et al. 2013), areas of possible improvement for both the data and the underlying models have been identified. As the data review was implemented, PNNL staff identified areas where the models can be improved both to accommodate the locally significant pathways identified and also to incorporate newer models. The areas are general data needs for the existing models and improved formulations for the pathway models.

54 ENVIRONMENTAL SCIENCES↗

SAM Code Validation on Frictional Pressure Drop through Pebble Beds

The System Analysis Module (SAM) is currently under development at Argonne National Laboratory as a modern system-level modeling and simulation tool for safety analyses of advanced non-light water reactors. This report presents a recent effort to validate the capability of SAM to predict the frictional pressure drop through pebble beds. Selected experimental data were used for code validation, including data from test facilities at Texas A&M University, Missouri University of Science and Technology, and North-West University of South Africa. SAM implements three empirical correlations to predict frictional pressure drop: the classical Ergun correlation; the KTA correlation, which is widely used in high-temperature gas-cooled reactor applications; and the Eisfeld and Schnitzlein correlation, which explicitly considers wall effect. Code validation was performed using all three correlations. For all selected experimental data, the KTA correlation shows the best performance and agrees very well with experimental measurement; the Eisfeld and Schnitzlein correlation, explicitly considering wall effects, shows acceptable accuracy, while there is no evidence that it is better than the KTA correlation; the Ergun correlation, however, over-predicts frictional pressure drop for most selected data points.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A trigonal nodal SP3 method with mesh refinement capabilities - Development and verification

The neutronics model of the nodal reactor dynamics code DYN3D developed for 3D analyses of steady states and transients in Light-Water Reactors has been extended by a simplified P{sub 3} (SP{sub 3}) neutron transport option - to overcome the limitations of the diffusion approach at regions with significant anisotropy effects. To provide a method being applicable to reactors with hexagonal fuel assemblies and to furthermore allow flexible mesh refinement, the nodal SP{sub 3} method has been developed on the basis of a flux expansion in triangular-z geometry. In this paper, the derivation of the trigonal SP{sub 3} method is presented in a condensed form and a verification of the methodology on quasi-pin level is performed by means of two single-assembly test examples. The corresponding pin-wise few-group cross sections were obtained by the deterministic lattice code HELIOS. The power distributions were calculated using both the trigonal DYN3D diffusion and SP{sub 3} solver and compared to the HELIOS reference solutions. Close to regions with non-negligible flux gradients, e.g., caused by the presence of a strong absorbing material, the power distribution calculated by DYN3D-SP{sub 3} shows a significant improvement in comparison to the diffusion method. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Comparison of a NuScale SMR conceptual core design using CASMO5/simulate5 and MCNP5

A key issue during the initial start-ups of new Small Modular Reactors (SMRs) is the lack of operational data for reactor model validation. To help better understand the accuracy of the reactor analysis codes CASMO5 and SIMULATE5, higher order comparisons to MCNP5 have been performed. These comparisons are for an initial core conceptual design of the NuScale reactor. The data have been evaluated at Hot Zero Power (HZP) conditions. Comparisons of core reactivity, fuel temperature coefficient (FTC), and moderator temperature coefficients (MTC) have been performed. Comparison results show good agreement between CASMO5/SIMULATE5 and MCNP5 for the conceptual initial core design. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

TRACE/PARCS analysis of the OECD/NEA Oskarshamn-2 BWR stability benchmark

On February 25, 1999, the Oskarshamn-2 NPP experienced a stability event which culminated in diverging power oscillations with a decay ratio of about 1.4. The event was successfully modeled by the TRACE/PARCS coupled code system, and further analysis of the event is described in this paper. The results show very good agreement with the plant data, capturing the entire behavior of the transient including the onset of instability, growth of the oscillations (decay ratio) and oscillation frequency. This provides confidence in the prediction of other parameters which are not available from the plant records. The event provides coupled code validation for a challenging BWR stability event, which involves the accurate simulation of neutron kinetics (NK), thermal-hydraulics (TH), and TH/NK. coupling. The success of this work has demonstrated the ability of the 3-D coupled systems code TRACE/PARCS to capture the complex behavior of BWR stability events. The problem was released as an international OECD/NEA benchmark, and it is the first benchmark based on measured plant data for a stability event with a DR greater than one. Interested participants are invited to contact authors for more information. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The experiment filtering tool in the MACSENS workbench - 14341

The MACSENS (outils Moret d'Aide a la Caracterisation et aux etudes de SENSibilite) workbench is dedicated to the advanced validation work for criticality studies. The aim of this kind of tool is to determine the experiments from a validation database being the 'closest' to an application case and, from there, to evaluate a calculation bias with a related uncertainty. To achieve this goal, the MACSENS workbench relies on different kind of information stored in various databases depending on the origin of the data. Two main kinds of data can be identified. The first one concerns the data used to describe the experiments of the validation database (e.g. the inputs needed to perform a calculation). These values characterize the chemical and geometrical models. They can be considered as constant, as they do not vary once evaluated, and some of them can be from a non-numerical nature (reflector material for example...). The second kind of data gathers the results obtained with the use of a simulation code, such as reaction rates, fluxes and leakages, sensitivities, k{sub eff}, etc. These values depend on the version of the code and are purely numerical. In the early 2000's, regarding the increasing number of experiments evaluated in the ICSBEP handbook, the OECD/NEA (Organization for Economic Co-operation and Development/Nuclear Energy Agency), under the ICSBEP (International Criticality Safety Benchmark Evaluation Project) working group direction, decided to store such information in a database and to provide a software, called DICE (Database for the International Criticality safety benchmark Evaluation project), allowing users to search easily criticality experiments of interest. Up to now, more than 4400 critical or subcritical configurations have been compiled by ICSBEP, described and stored in the DICE database. Over time, a validation work has been done on the data, with the participation of IRSN, and a quality assurance policy has been defined. For these reasons, the DICE database is highly reliable. These features were fundamental for the MACSENS developers, who were looking for a database containing descriptive data of the experiments included in the MORET code validation database. The paper describes the solution chosen by IRSN to take advantage of the DICE database content and to overcome the DICE limits. It also explains the quality assurance policy applied to guaranty the reliability of the stored data. It will also discuss the use of the experiment-filtering module of the MACSENS workbench to select a set of experiments and to examine numerical parameter trends for an application case. (authors)

42 ENGINEERING↗

Recent Results on Coupling Fast Waves to High Performance Plasmas on DIII-D

Fast Waves (FWs) at 60 MHz and 90 MHz are used in DIII-D for central electron heating and current drive. Coupling of FWs to high-performance discharges is limited by low antenna loading in these regimes. To extend the application of high-power FWs to such regimes, methods of increasing the antenna loading in these regimes are needed. A systematic study of loading enhancement techniques has been carried out in DIII-D, including reduction of the antenna/plasma distance, gas puffing into the far scrape-off layer (SOL), and control of other parameters that affect the particle balance in the far SOL. Quantitative understanding of the physics of the loading resistance and its dependence on edge density profiles is demonstrated. The core FW heating efficiency appeared to be {approx}100% in the Advanced Inductive regime, consistent with the high first-pass direct electron absorption of {approx}75% that is predicted by the ray-tracing code GENRAY in this high electron beta regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Diffusive deposition of aerosols in Phebus containment during FPT-2 test

At present the lumped-parameter codes is the main tool to investigate the complex response of the containment of Nuclear Power Plant in case of an accident. Continuous development and validation of the codes is required to perform realistic investigation of the processes that determine the possible source term of radioactive products to the environment. Validation of the codes is based on the comparison of the calculated results with the measurements performed in experimental facilities. The most extensive experimental program to investigate fission product release from the molten fuel, transport through the cooling circuit and deposition in the containment is performed in PHEBUS test facility. Test FPT-2 performed in this facility is considered for analysis of processes taking place in containment. Earlier performed investigations using COCOSYS code showed that the code could be successfully used for analysis of thermal-hydraulic processes and deposition of aerosols, but there was also noticed that diffusive deposition on the vertical walls does not fit well with the measured results. In the CPA module of ASTEC code there is implemented different model for diffusive deposition, therefore the PHEBUS containment model was transferred from COCOSYS code to ASTEC-CPA to investigate the influence of the diffusive deposition modelling. Analysis was performed using PHEBUS containment model of 16 nodes. The calculated thermal-hydraulic parameters are in good agreement with measured results, which gives basis for realistic simulation of aerosol transport and deposition processes. Performed investigations showed that diffusive deposition model has influence on the aerosol deposition distribution on different surfaces in the test facility. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Uncertainties propagation in the framework of a Rod Ejection Accident modeling based on a multi-physics approach

The control of uncertainties in the field of reactor physics and their propagation in best-estimate modeling are a major issue in safety analysis. In this framework, the CEA develops a methodology to perform multi-physics simulations including uncertainties analysis. The present paper aims to present and apply this methodology for the analysis of an accidental situation such as REA (Rod Ejection Accident). This accident is characterized by a strong interaction between the different areas of the reactor physics (neutronic, fuel thermal and thermal hydraulic). The modeling is performed with CRONOS2 code. The uncertainties analysis has been conducted with the URANIE platform developed by the CEA: For each identified response from the modeling (output) and considering a set of key parameters with their uncertainties (input), a surrogate model in the form of a neural network has been produced. The set of neural networks is then used to carry out a sensitivity analysis which consists on a global variance analysis with the determination of the Sobol indices for all responses. The sensitivity indices are obtained for the input parameters by an approach based on the use of polynomial chaos. The present exercise helped to develop a methodological flow scheme, to consolidate the use of URANIE tool in the framework of parallel calculations. Finally, the use of polynomial chaos allowed computing high order sensitivity indices and thus highlighting and classifying the influence of identified uncertainties on each response of the analysis (single and interaction effects). (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

PHISICS toolkit: Multi-reactor transmutation analysis utility - MRTAU

The principal idea of this paper is to present the new capabilities available in the PHISICS toolkit, connected with the implementation of the depletion code MRTAU, a generic depletion/decay/burn-up code developed at the Idaho National Laboratory. It is programmed in a modular structure and modern FORTRAN 95/2003. The code tracks the time evolution of the isotopic concentration of a given material accounting for nuclear reactions happening in presence of neutron flux and also due to natural decay. MRTAU has two different methods to perform the depletion calculation, in order to let the user choose the best on with respect to his needs. Both the methodologies and some significant results are reported in this paper. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The Data Transfer Kit: A geometric rendezvous-based tool for multiphysics data transfer

The Data Transfer Kit (DTK) is a software library designed to provide parallel data transfer services for arbitrary physics components based on the concept of geometric rendezvous. The rendezvous algorithm provides a means to geometrically correlate two geometric domains that may be arbitrarily decomposed in a parallel simulation. By repartitioning both domains such that they have the same geometric domain on each parallel process, efficient and load balanced search operations and data transfer can be performed at a desirable algorithmic time complexity with low communication overhead relative to other types of mapping algorithms. With the increased development efforts in multiphysics simulation and other multiple mesh and geometry problems, generating parallel topology maps for transferring fields and other data between geometric domains is a common operation. The algorithms used to generate parallel topology maps based on the concept of geometric rendezvous as implemented in DTK are described with an example using a conjugate heat transfer calculation and thermal coupling with a neutronics code. In addition, we provide the results of initial scaling studies performed on the Jaguar Cray XK6 system at Oak Ridge National Laboratory for a worse-case-scenario problem in terms of algorithmic complexity that shows good scaling on 0(1 x 104) cores for topology map generation and excellent scaling on 0(1 x 105) cores for the data transfer operation with meshes of O(1 x 109) elements. (authors)

97 MATHEMATICS AND COMPUTING↗

Monte Carlo calculations of k{sub Q}, the beam quality conversion factor

Purpose: To use EGSnrc Monte Carlo simulations to directly calculate beam quality conversion factors, k{sub Q}, for 32 cylindrical ionization chambers over a range of beam qualities and to quantify the effect of systematic uncertainties on Monte Carlo calculations of k{sub Q}. These factors are required to use the TG-51 or TRS-398 clinical dosimetry protocols for calibrating external radiotherapy beams. Methods: Ionization chambers are modeled either from blueprints or manufacturers' user's manuals. The dose-to-air in the chamber is calculated using the EGSnrc user-code egs{sub c}hamber using 11 different tabulated clinical photon spectra for the incident beams. The dose to a small volume of water is also calculated in the absence of the chamber at the midpoint of the chamber on its central axis. Using a simple equation, k{sub Q} is calculated from these quantities under the assumption that W/e is constant with energy and compared to TG-51 protocol and measured values. Results: Polynomial fits to the Monte Carlo calculated k{sub Q} factors as a function of beam quality expressed as %dd(10){sub x} and TPR{sub 10}{sup 20} are given for each ionization chamber. Differences are explained between Monte Carlo calculated values and values from the TG-51 protocol or calculated using the computer program used for TG-51 calculations. Systematic uncertainties in calculated k{sub Q} values are analyzed and amount to a maximum of one standard deviation uncertainty of 0.99% if one assumes that photon cross-section uncertainties are uncorrelated and 0.63% if they are assumed correlated. The largest components of the uncertainty are the constancy of W/e and the uncertainty in the cross-section for photons in water. Conclusions: It is now possible to calculate k{sub Q} directly using Monte Carlo simulations. Monte Carlo calculations for most ionization chambers give results which are comparable to TG-51 values. Discrepancies can be explained using individual Monte Carlo calculations of various correction factors which are more accurate than previously used values. For small ionization chambers with central electrodes composed of high-Z materials, the effect of the central electrode is much larger than that for the aluminum electrodes in Farmer chambers.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Mathematical framework for the analysis of dynamic stochastic systems with the RAVEN code

RAVEN (Reactor Analysis and Virtual control Environment) is a software code under development at Idaho National Laboratory aimed at performing probabilistic risk assessment and uncertainty quantification using RELAP-7, for which it acts also as a simulation controller. In this paper we will present the equations characterizing a dynamic stochastic system and we will then discuss the behavior of each stochastic term and how it is accounted for in the RAVEN software design. Moreover we will present preliminary results of the implementation. (authors)

97 MATHEMATICS AND COMPUTING↗

Optical boundary reconstruction of tokamak plasmas for feedback control of plasma position and shape

A new diagnostic is developed to reconstruct the plasma boundary using visible wavelength images. Exploiting the plasma's edge localized and toroidally symmetric emission profile, a new coordinate transform is presented to reconstruct the plasma boundary from a poloidal view image. The plasma boundary reconstruction is implemented in MATLAB and applied to camera images of Mega-Ampere Spherical Tokamak discharges. The optically reconstructed plasma boundaries are compared to magnetic reconstructions from the offline reconstruction code EFIT, showing very good qualitative and quantitative agreement. Average errors are within 2 cm and correlation is high. In the current software implementation, plasma boundary reconstruction from a single image takes 3 ms. The applicability and system requirements of the new optical boundary reconstruction, called OFIT, for use in both feedback control of plasma position and shape and in offline reconstruction tools are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗