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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Gas Gun Model and Comparison to Experimental Performance of Pipe Guns Operating with Light Propellant Gases and Large Cryogenic Pellets

Injection of multiple large (~10 to 30 mm diameter) shattered pellets into ITER plasmas is presently part of the scheme planned to mitigate the deleterious effects of disruptions on the vessel components. To help in the design and optimize performance of the pellet injectors for this application, a model referred to as “the gas gun simulator” has been developed and benchmarked against experimental data. The computer code simulator is a Java program that models the gas-dynamics characteristics of a single-stage gas gun. Following a stepwise approach, the code utilizes a variety of input parameters to incrementally simulate and analyze the dynamics of the gun as the projectile is launched down the barrel. Using input data, the model can calculate gun performance based on physical characteristics, such as propellant-gas and fast-valve properties, barrel geometry, and pellet mass. Although the model is fundamentally generic, the present version is configured to accommodate cryogenic pellets composed of H2, D2, Ne, Ar, and mixtures of them and light propellant gases (H2, D2, and He). The pellets are solidified in situ in pipe guns that consist of stainless steel tubes and fast-acting valves that provide the propellant gas for pellet acceleration (to speeds ~200 to 700 m/s). The pellet speed is the key parameter in determining the response time of a shattered pellet system to a plasma disruption event. The calculated speeds from the code simulations of experiments were typically in excellent agreement with the measured values. With the gas gun simulator validated for many test shots and over a wide range of physical and operating parameters, it is a valuable tool for optimization of the injector design, including the fast valve design (orifice size and volume) for any operating pressure (~40 bar expected for the ITER application) and barrel length for any pellet size (mass, diameter, and length). Key design parameters and proposed values for the pellet injectors for the ITER disruption mitigation systems are discussed.

Combs, S. K.↗
DOE OSTI2013

High-Efficacy Modeling of Human Perceptual Similarity Opinions

The computer code relates to human nature studies and more specifically to a novel graphical user interface (GUI) that facilitates high-efficacy collection of perceptual similarity opinions of a user in an effective and intuitive manner. The GUI presents a base image for rating its similarity to seven peripheral images that are displayed simultaneously following a circular layout. The user is asked to report the base image's pairwise similarity to each peripheral image on a fixed scale while preserving the relative ranking among all peripheral images. The collected data are then used to predict the user's subjective opinions regarding the perceptual similarity of images.

Xu, Songhua↗
DOE OSTI2016

An Analytical Methodology to Evaluate Thermal-Hydraulic Performance of Compact Heat Exchangers, Accounting for Heat Loss

High efficiency and compactness, offered by the Printed Circuit Heat Exchanger (PCHE) makes it a strong candidate for use as an Intermediate Heat Exchanger (IHX) and/or Secondary Heat Exchanger (SHX) for Advanced High Temperature Reactor (AHTR) systems and also the Very High Temperature Reactor (VHTR). Idaho National Laboratory (INL) has been developing a computer code for the PCHE design and analysis [1], and will construct an experimental loop to evaluate the thermal hydraulic performance of PCHEs [2]. In order to correctly compare experimental results with the simulations, the extraneous heat losses from the PCHEs must be determined and incorporated. These heat exchangers operate at high temperature and consequently, despite the use of external thermal insulation, experience non-negligible heat losses to the surroundings. The objective of this work is to develop an analytical methodology for quantification of these heat losses based on experimental data to evaluate and predict the performance of the heat exchanger accurately.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

JMOCUP Physics Depletion Calculations for the As-Run AGR-5/6/7 TRISO Particle Experiment in ATR Northeast Flux Trap

This ECAR documents the as-run Jim Sterbentz’s MCNP ORIGEN Coupled Utility Program (JMOCUP) physics depletion calculation and the calculated results for the AGR-5/6/7 irradiation experiment in the Advanced Test Reactor (ATR). AGR 5/6/7 was irradiated for nine power cycles in the northeast flux trap. The depletion calculations were performed to provide input data for a variety of other engineering analyses supporting the AGR-5/6/7 experiment along with post-irradiation characterization of the tri-structural isotropic (TRISO) particle fuel compacts. Detailed full-core MCNP models and Oak Ridge Isotope Generation (ORIGEN2) radionuclide generation models were specifically developed as part of the JMOCUP Monte Carlo depletion calculations. The MCNP ORIGEN2 computer codes were coupled using the well-established JMOCUP utility modules to couple the two codes and run the depletion calculations. The physics calculations were performed in support of the Advanced Gas Reactor (AGR) program.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

RELAP5/MOD2 Overview and Developmental. Assessment Results from TMl-1 Plant Transient Analysis

RELAP5/MOD2 is a new version of the RELAP5 thermal-hydraulic computer code containing improved modeling features that provide a generic capability for pressurized water reactor transient simulation. The objective of this paper is to provide code users with an overview of the code and to report developmental assessment results obtained from a Three Mile Island Unit One plant transient analysis. The assessment shows that the injection of highly sub-cooled water into a high-pressure primary coolant system does not cause unphysical results or pose a problem for RELAP5/MOD2. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Education and training in the field of nuclear instrumentation and measurement: CEA/INSTN (National Institute for Nuclear Sciences and Technologies) strategy to improve and develop new pedagogical tools and methods

Part of the French Alternative Energies and Atomic Energy Commission (CEA), the National Institute for Nuclear Science and Technology (INSTN) is a higher education institution whose mission is to provide students and professionals a high level of scientific and technological qualification in all disciplines related to nuclear energy applications. In this frame, INSTN carries out education and training (E and T) programs in nuclear instrumentation and radioprotection. Its strategy has always been to complete theoretical courses by training courses and laboratory works carried out on an extensive range of training tools that includes a large panel of nuclear instrumentation as well as software applications. Since its creation in 1956, the INSTN has conducted both education and vocational programs on ionizing radiation detection. An extensive range of techniques have commonly been used during practical works with students and employees of companies who need to get the knowledge and specialization in this field. Today, the INSTN is mainly equipped with usual detectors and electronics in large numbers in order to be able to accommodate up to 48 trainees at the same time in two classrooms, with only two trainees for one workstation in order to optimize their learning. In the field of the neutron detection systems, the INSTN has strongly developed its offer taking advantage of the use of research reactors, such as ISIS reactor (700 kW) at Saclay. The implementation of neutron detection systems specific to the courses offers a unique way of observing and analysing the signal coming from neutron detectors, as well as learning how to set the parameters of the detection system in real conditions. Providing the trainees with an extensive overview of each part of the neutron monitoring instrumentation apply to a nuclear reactor, hands-on measurements on the ISIS reactor play a major role in ensuring a practical and comprehensive understanding of the neutron detection system and their integration in the safety system of nuclear reactors. It also gives a solid background for the follow up and the development of the neutron detection systems. Another field of activity of the INSTN is the development of new teaching tools using software applications. We describe hereafter two applications that have recently been developed by the institute, O.S.I.R.I.S and DOSIMEX. The INSTN and the company OREKA have developed an innovative teaching tool named O.S.I.R.I.S. The tool is built on a virtual 3D environment in which users operate in a totally free way. The action is located in a pressurised water reactor in the environment of a steam generator building. The users are students or professionals who want to learn and train on concrete situations on how to protect the workers against radiation. Through the use of this serious game, users must: establish a predictive dose evaluation, implement the principles of radiation protection, perform the dose result of the operation, analyse the gap between the predictive collective dose and the collective dose achieved, and decide about the lessons to provide feedback. Taking into account the increasing constraints related to radiation protection, it is essential to use computer codes to evaluate the dose rates generated by sources of ionizing radiation. Many powerful codes exist but are often relatively complex to implement and are dedicated to experts. In addition, these codes are often 'black boxes' that does not allow the user to understand the underlying physics. The objective of DOSIMEX codes developed by the INSTN is to give, by the means of a VBA code with a simple interface, the ability to easily conduct a wide range of radiological exposure situations: calculation of gamma, beta, neutron equivalent dose rate, calculation of internal contamination and atmospheric transfer.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

The Acoustic Emission signal acquired by the microphones placed in the CABRI test device along the fourteen last R.I.A. experiments: an example of reproducible research in nuclear science

The CABRI facility is an experimental nuclear reactor of the French Atomic Energy Commission. It is located at the Cadarache Research Centre in southern France and it is designed to act as a support to the French nuclear infrastructure. The purpose of the new testing programme termed, 'CABRI International Programme' (CIP) is to study the behaviour of PWR-type fuel rods at high burnup, equipped with an 'advanced' cladding, under Reactivity Initiated Accident (RIA) conditions (such as the scenario of a control rod ejection). Within the framework of this programme, piloted and funded by the French Institute of Nuclear Radioprotection and Safety (IRSN), ten tests are to be conducted with a frequency of two tests per year. The LPRE laboratory of the CEA which is in charge of the Preparation, realisation and breakdown of the test results studies the possibility to set up a new test analysis based on the processing of signals coming from sensors placed within the test equipment. During the experimental phase, the behaviour of the fuel element generates acoustic waves which can be detected by two microphones placed upstream and downstream of the test device. Studies showed the interest to realize temporal and spectral analyses on these signals by showing the existence of signatures which can be correlated with physical phenomena as the rod failure or the test shutdown (i.e. scram). The work presented in this article results from the will to consolidate these studies. Since the main phenomenon to be tracked is the fuel rod failure, the aim would be to highlight specific events which would have been precursors of the rod failure in order to use in the future these signals for further interpretation. Such an antecedent works resumption leads to a better understanding of the experimental needs and constitutes a good initial state to prepare the new very fast digital data acquisition systems. Although all the raw data are accessible in the form of text files, analyses and graphics representations were not straightforward to reproduce from the ancient studies since that, on one hand, people who were in charge of the original work left the laboratory and on the other hand because it is not easy when the time passes, even with our own work, to be able to remember the steps of data manipulations and the exact setup: - During the ancient experiments the use of analog data acquisition systems required to digitize tapes to be able to realize computer treatments. That had had for consequence to lose the initial dating. This one must be correctly edited to do temporal comparisons. - Analyses require functions for calculations whose parameters has to be well-known to reach the same results. We thus wished to manage our workflow in the idea that it can be easily reproducible on all the experiments. The object of the work presented in this article was to put in practice this strong bind between the data, treatments and generation of the document in order not to hesitate to do the iteration principle in action. We do not have to be afraid by the data driven analyses. According to the philosophy of the literate programming, the text of the technical document is woven with the computer code that produces all the printed output as tables, graphs for the study eliminating hence the unrealistic cut and paste. This difficulty is not specific to the nuclear domain. For many years, researchers have been worked out solutions to this mundane issue. And, presently, new technologies and high-level programming languages offer us actual answers. We will firstly present the tools applied in our laboratory to implement this workflow, then we will describe the global perception carried out to continue the study of the Acoustic Emission signals recorded by the two microphones during the fourteen last CABRI R.I.A. test.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

First Industrial Tests of a Drum Monitor Matrix Correction for the Fissile Mass Measurement in Large Volume Historic Metallic Residues with the Differential Die-away Technique

The fissile mass in radioactive waste drums filled with compacted metallic residues (spent fuel hulls and nozzles) produced at AREVA La Hague reprocessing plant is measured by neutron interrogation with the Differential Die-away measurement Technique (DDT. In the next years, old hulls and nozzles mixed with Ion-Exchange Resins will be measured. The ion-exchange resins increase neutron moderation in the matrix, compared to the waste measured in the current process. In this context, the Nuclear Measurement Laboratory (NML) of CEA Cadarache has studied a matrix effect correction method, based on a drum monitor ({sup 3}He proportional counter inside the measurement cavity). A previous study performed with the NML R and D measurement cell PROMETHEE 6 has shown the feasibility of method, and the capability of MCNP simulations to correctly reproduce experimental data and to assess the performances of the proposed correction. A next step of the study has focused on the performance assessment of the method on the industrial station using numerical simulation. A correlation between the prompt calibration coefficient of the {sup 239}Pu signal and the drum monitor signal was established using the MCNPX computer code and a fractional factorial experimental design composed of matrix parameters representative of the variation range of historical waste. Calculations have showed that the method allows the assay of the fissile mass with an uncertainty within a factor of 2, while the matrix effect without correction ranges on 2 decades. In this paper, we present and discuss the first experimental tests on the industrial ACC measurement system. A calculation vs. experiment benchmark has been achieved by performing dedicated calibration measurement with a representative drum and {sup 235}U samples. The preliminary comparison between calculation and experiment shows a satisfactory agreement for the drum monitor. The final objective of this work is to confirm the reliability of the modeling approach and the industrial feasibility of the method, which will be implemented on the industrial station for the measurement of historical wastes. (authors)

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

Use of Monte Carlo simulation to quantify CIPS risk

Crud Induced Power Shift (CIPS), caused by the accumulation of boron compounds in corrosion product deposits on high power fuel assemblies, is an important operational issue for Pressurized Water Reactors (PWRs). The Boron-induced Offset Analysis (BOA) computer code has been developed by EPRI in collaboration with Westinghouse to help utility engineers avoid onset of CIPS. The code, used with an established procedure, predicts the maximum boron mass that will be accumulated in the upcoming cycle for a given core design and plant chemistry and operating program. The calculated boron mass above a given threshold predicts the risk of CIPS to the upcoming cycle. BOA has been successful in correctly predicting the risk for onset or lack of CIPS in many cycles. However, like most codes modeling complex physical phenomena, it is subject to uncertainties in several of the inputs used to describe these phenomena. The current BOA procedure results in a single 'best estimate' prediction, which has been judged to provide engineers with a practical way to assess CIPS risk. As such, however, there is a probability that the actual plant CIPS performance may be underestimated. This paper describes a method to assess the probability (with an associated confidence level) of no CIPS, mild to moderate CIPS, or severe CIPS using BOA for a given core design. This is accomplished through the use of a Monte Carlo simulation. The simulation is applied to several plant examples, and the results are evaluated in terms of the likelihood of mild or severe CIPS. A suggested set of criteria for decision making regarding CIPS risk is presented in the conclusion. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

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

Application of the FW-CADIS variance reduction method to calculate a precise N-flux distribution for the FRJ-2 research reactor - 15603

In order to avoid unnecessary radiation exposure and to minimize the radioactive waste and respectively cost by the decommissioning of nuclear facilities, the radiological characteristics and the activity atlases are required. Therefore, a significant step is the accurate determination of the neutron fluence distribution resulting from the operation of a reactor. For this purpose, computer codes based on the Monte-Carlo method are employed because of the flexibility in the modelling of the complex structures and physical processes. In the scope of the research project founded by the German Federal Ministry of Education and Research (BMBF), and implemented by the Institute of Nuclear Engineering and Technology Transfer (NET) -RWTH Aachen, sophisticated simulation tools are developed, which have the capability to cover the range of demands for activity and dose rate distribution using the example of the German research reactor FRJ-2. Due to the complex geometry of FRJ-2, which consists of different structures and configuration, the MCNP5 Monte Carlo Code is the optimal method for a precise neutron transport simulation. This work presents a method to achieve a precise N-Flux distribution for the FRJ-2 research reactor using hybrid deterministic Monte Carlo method. (authors)

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

Limiting Surface Density Method Application to Large Arrays of Highly Heterogeneous Shipping Packages

This work represents the initial phase of an overall effort to reduce the analytical burdens and increase flexibility of the K-Area Complex which provides stewardship of special nuclear materials. A hand calculation method is being extended to apply to large arrays of highly heterogeneous shipping packages. This method has the potential to quickly determine the impact of perturbations on the array. The objective of the particular work reported herein is twofold. First, the analytical results upon which the Limiting Surface Density (LSD) method was originally based were reproduced using the modern computational code KENO-VI in SCALE 6.1. The second objective was to test the classic LSD method on storage of large arrays of highly heterogeneous shipping packages. Specifically, the objective is to determine if the method is adequate for existing arrays of 9975 shipping packages with fissile metal content. The classic LSD method was derived for simple packages of bare spheres of fissile material centered in an air filled cubic unit cell and placed in a cubic-pitch water-reflected array, the properties of which could then be translated to more practical array and package geometries. The classic method did not address highly heterogeneous packages which are close to neutronically decoupled from other identical packages. Therefore, the material based array constant (called c{sub 2}) as defined by the classic method is determined specifically for the packages in question. (authors)

42 ENGINEERING↗
DOE OSTI2019

Radiological Performance Modeling: Platform Comparison

Zapadnoe Uranium Mill Tailings Site, Pridneprovsky Chemical Plant Ukraine: Two exposure scenarios were utilized to develop performance assessment models in GoldSim and NORMALYSA. The first scenario is an onsite worker exposed to radioactivity through external ground irradiation as well as radon and radioactive dust inhalation. The second being an offsite Resident 300 years after institutional control exposed through the ingestion of garden vegetables irrigated by contaminated ground water, ingestion of fish from Konoplyanka River, and agricultural activities in contaminated garden. NORMALYSA utilizes a modular interface to connect extensive pre-created modules to develop contaminant transport and risk assessment models for safety assessments and performance assessments at NORM and legacy sites. Developed by Facilia AB, the program was tested with the creation of the Zapadnoe model under the International Atomic Energy Agency MODARIA Il Working Group 1 (WG1). Default values are provided for input parameters where site specific data and distributions lack. Probabilistic simulations and inputs are limited to parameters with assigned Gaussian probability functions. Similar results for transient worker inhalation dose were found in comparison with the Zapadnoe GoldSim model. However, the NORMALYSA model has results larger by an order of magnitude for external radiation in the transient worker and for aquatic food ingestion, inhalation, and garden food ingestion to the resident. Performance Assessment Modeling: Radiological performance models promote good management and support in remediation and environmental decision making by optimizing radioactive waste disposal site performance. a performance assessment model of the Zapadnoe uranium mill tailings site and a model for calculating authorized limits at Los Alamos National Laboratory are created to demonstrate and highlight the respective capabilities of three modeling platforms: GoldSim, NORMALYSA, and RESRAD-Onsite. GoldSim is a general-purpose modeling platform with the ability to simulate site-specific FT radionuclide transport and risk assessment using a graphical interface to create and link elements that act as indivisible building blocks of the model. The Zapadnoe performance assessment model created in GoldSim is composed of input parameters collected from site specific parameters from WG1 as well as standard assumptions seen in the Generic PA Model from Neptune and Company. The multivariate Zapadnoe simulation contains both deterministic and probabilistic simulation capabilities with the ability to quantitatively represent the uncertainty and associated risks of the site via the Monte Carlo simulation method. Model parameters, documentation, and assumptions are built directly into the model creating a transparent, reproducible, and more purposeful site specific model. To demonstrate this, an ingestion of water pathway was included in addition to the exposure pathways included in the NORMALYSA model. GoldSim offers a wide variety of probability distributions to choose from as well as many tools to visualize the uncertainty of the models and explore which pathways and radionuclides dominate dose. In contrast to the NORMALYSA model, the GoldSim model uses dynamic simulation to calculate tritium concentrations. Authorized limits support the evaluation of radiological health effects from radionuclides present in bulk contaminated soil in order to assist in the process of clearing land for release under different uses. Calculations of authorized limits in RESRAD-Onsite and GoldSim models support the evaluation for 15 radionuclides using a residential exposure scenario, indoor and outdoor inhalation of dust and vapor phase radionuclides (tritium), incidental soil ingestion, and ingestion grown fruits and Vegetables. RESRAD-Onsite, developed by Argonne National Laboratory is a dedicated computer code for calculating radiation doses and cancer risks to an individual located on top of radioactively contaminated soils. It can also be used to provide guidelines for deriving radionuclide soil levels corresponding to a specific dose criterion. The menu driven user interface facilitates selection of radionuclides and soil concentrations, as well as defines parameter values related to materials, environmental factors, and exposure pathways. The program contains radionuclide-specific data such as dose conversion factors and risk coefficients, as well as soil-water partition coefficients and transfer factors for biotic media. The GoldSim model successfully emulates the RESRAD-Onsite code using identical radiological decay and ingrowth factors, radiation dose conversion factors, and other input parameters. DE EEDE ETS Each model supports dynamic modeling of radionuclide ingrowth and decay over time, and integrate this capability with the transport of radionuclides in the environment via processes such as advection and/or diffusion in fluids, or within porous solid media, The applicability of each program is heavily dependent on project context. In the context of developing Authorized Limits, RESRAD-Onsite may be a more desirable option as it was developed for this particular application. For simple disposal or legacy sites or sites that do not need site-specific model structure, NORMALYSA may be a desirable modeling platform, GoldSim can be used to model any problem for which RESRAD or NORMALYSA is applicable, but its strength lies in its flexibility, analytical, and probabilistic capabilities.

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

Passive Safety Feature of KIPT Neutron Source Facility

Argonne National Laboratory (ANL) of the United States and Kharkov Institute of Physics and Technology (KIPT) of Ukraine have cooperated on the development, design, and construction of a neutron source facility. The facility was constructed at Kharkov, Ukraine and its commissioning process is underway. The facility has an electron accelerator driven subcritical assembly. The electron beam power is 100 kW using 100 MeV electrons. Tungsten or natural uranium is the target material for generating neutrons driving the subcritical assembly. The subcritical assembly is composed of WWR-M2 Russian fuel assemblies with U-235 enrichment of 19.7 wt%, surrounded by beryllium reflector assembles and graphite blocks. The water coolant operates at room temperature during normal operation, when the facility is operating at full power using 100-KW electron beam power. The passive safety features of the facility are discussed in in this study. Monte Carlo computer code MCNPX was utilized in the analyses with ENDF/BVII.0 nuclear data libraries. Negative reactivity temperature feedback was consistently observed, which is important for the facility safety performance. Due to the design of WWR-M2 fuel assemblies, slight temperature increase and the corresponding density decrease in the water coolant would cause relatively large reactivity drop, which offset the reactivity gain by mistakenly loading an additional fuel assembly. The increase of fuel temperature also causes sufficiently large reactivity decrease. This enhances the facility safety performance because fuel temperature increase provides prompt negative reactivity feedback. (authors)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

Particle Image Velocimetry Measurements in a Wire-Wrapped 61-Pin Hexagonal Fuel Bundle

Liquid metal fast reactor fuel assemblies using sodium as a coolant typically utilize a tightly packed triangular lattice of fuel pins enclosed in a hexagonal duct. The fuel pins are typically wrapped with a helical wire spacer. The primary function of the wire spacers is to maintain a gap between neighboring fuel pins. They also mitigate vortex-induced vibration and enhance sub-channel mixing to increase convective heat transfer. Many experiments have been performed to study thermal-hydraulic characteristics, primarily pressure drops, of the wire-wrapped fuel bundles. Novendstern developed a semi-empirical model to predict pressure losses in a hexagonal array of wire-wrapped pins in the turbulent flow regime. Rehme proposed pressure drop correlations for a wide range of geometrical parameters, such as pitch-to-diameter ratio P/D from 1.125 to 1.417, lead-to-wire diameter ratio H / dm from 6 to 45, and number of pins n from 7 to 61. Cheng and Todreas introduced their hydrodynamic models for subchannel friction factors and mixing parameters in a 37-pin fuel bundle, which covered the laminar, transition and turbulent flow regimes. Computational fluid dynamics (CFD) offers an opportunity to validate the above correlations and investigate the flow characteristics in the subchannels. Several numerical studies of the wire-wrapped fuel bundles have been performed and most of them are Reynolds-Averaged Navier-Stokes (RANS) based simulations with various turbulence models. Ahmad and Kim performed RANS simulations of 7- and 19-pin bundles using the k-omega SST turbulent model, while Gajapathy et al. performed RANS simulations on 7-pin bundles using the k-epsilons model. Recently, Fischer et al. carried out large-eddy simulation (LES) calculations on wire-wrapped fuel bundles to obtain higher-resolution numerical solutions of subchannel behavior. The geometrical complexity of the hexagonal wire-wrapped fuel bundle has limited experimental activities to pressure drop measurements. A database of pressure drops and flow-field measurements in a 61-pin wire-wrapped hexagonal fuel bundle is not available to benchmark the existing correlations and validate the CFD calculations. Texas A and M University is conducting isothermal flow experiments in a wire-wrapped 61-pin hexagonal fuel bundle to support the research on advanced nuclear fuel development sponsored by the Department of Energy (DOE). The purpose of these tests is to perform high spatial and temporal resolution measurements of the flow and pressure measurements at different locations in the fuel bundle. The experimental facility employs matched-index-of-refraction (MIR) techniques and laser diagnostic velocity measurement techniques such as particle image velocimetry (PIV) and laser Doppler velocimetry (LDV). The result of this experimental activity will be an experimental database of pressure and flow-field measurements that will be suitable for validating CFD computer codes. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Best-Estimate-Plus-Uncertainty-Informed Deterministic Safety Evaluation

Nuclear power plant (NPP) utilities operate contingent on receiving a favorable safety assessment from their national regulatory and safety authority. Among the many facets of NPP safety addressed is verification of integral design-basis performance - a demonstration of a plant's response to a suite of challenging transient and accident scenarios. Thermal-hydraulic system computer codes like RELAP5, TRAC, CATHARE, and ATHLET are used as the simulation engine as part of a broader evaluation model that accounts for a variety of uncertainties. It is common for these evaluation models to apply either a super-bounding/ deterministic or best-estimate plus uncertainty (BEPU) approach. The trend over the past two decades has been away from the classical deterministic methods towards methods explicitly incorporating parametric uncertainty models. These methods center upon the evaluation of parametric and system response uncertainties, where these uncertainties are derived from the convolution of one or more parametric uncertainty models. Acceptance of BEPU methods has provided the nuclear industry benefits regarding operational flexibility and, in some cases, power output. The advantage of BEPU for new NPPs is less obvious. Generally, these new designs include engineered safety features that practically eliminate fuel damage resulting from design-basis events. To substantiate such claims, applying super-bounding assumptions on uncertainties may be a more expedient demonstration of a plant's safety case than those based on BEPU. In addition, the requirements on design inputs for a deterministic evaluation model are substantially less than that for a BEPU evaluation model. As such, applying a deterministic evaluation model is expected to have the advantage of reducing the duration of the licensing review process. For the current operating fleet of light-water reactors (LWRs), the nuclear safety authorities in most countries have followed the U.S. Nuclear Regulatory Commission's rules for the preparation of safety analysis. These rules are specified in the U.S. Code of Federal Regulations, most notably in, Part 10, Section 50.46, entitled 'acceptance criteria for emergency core cooling systems for light-water nuclear power reactors' and its Appendix K. Notably, Appendix K specifies deterministic requirements for analysis of the loss-of-coolant accident (LOCA), which is generally regarded as the maximum credible accident for LWRs. Projecting 10 CFR 50, Appendix K, to new plants is not obviously applicable as it addresses light water fluid flow and heat transfer phenomena associated with core uncovery. Preventing core uncovery is a requirement for several designs in the class of small modular reactors (SMR). In doing so, the majority of phenomenological uncertainties prescribed in Appendix K are eliminated. Similarly, the impact of other phenomena to effect a change in various safety measures, i.e., figures-of-merit, is expected to be different. The implication is that development of a deterministic evaluation model must be derived with an original technical basis rather than solely rely on Appendix K. While a logical conclusion, the expectation of original bases, regardless of the treatment of uncertainties, is embodied in the U.S. NRC's Regulatory Guide (RG) 1.203, Transient and Accident Analysis. Within the evaluation model development and application process framework of RG 1.203, this paper describes a BEPU-informed deterministic evaluation model for NPPs with large inherent safety margins. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Second-Order Uncertainty Analysis for Radiation Shielding Problems Using MCNP

Standard uncertainty propagation is based on a first-order Taylor expansion. This approximation is usually justified by assuming that 'uncertainties are very small numbers compared with the values of their corresponding quantities'. Furthermore, the second derivatives needed for a second-order analysis are not usually available, and a second-order analysis requires more knowledge regarding the distribution of the random variables. For radiation shielding problems, the MCNP6 Monte Carlo code computes first and second derivatives of a response with respect to nuclide density or reaction cross section. This paper explores the application of these derivatives to second-order uncertainty analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Uncertainty Quantification and Sensitivity Analysis for Coupled VERA-CS and FRAPCON

The current nuclear power plant fleet in the U.S. is facing numerous challenges economically ranging from low natural gas price to the rapidly growing renewable energy sources such as wind and solar power. In order for the nuclear industry to keep the viability of its operating fleet, new innovations have to be introduced to lower the costs to operate the fleet. One area that could contribute to the cost reduction is to lower the fuel cycle cost and improve the plant safety through initiatives such as higher burnup fuel and accident tolerant fuel. At the same time, the nuclear industry has to address future challenges on the horizon such as load following and other flexible operating strategies as well as those associated with new regulations such as 10 CFR 50.69 and 10 CFR50.46c etc. All these require coupled analysis of core neutronics, thermal hydraulics and fuel performance, which necessitates the development of a multi-physics toolkit. Since different models, assumptions, and fidelities have gone into various computer codes, how to consistently propagate uncertainties in a multi-physics environment warrants careful evaluation. The U.S. DOE's Light Water Reactor Sustainability (LWRS) Program initiated a project to develop analytical capabilities to support the industry to address these challenges. The general idea behind the initiative was the development of an integrated multiphysics evaluation model. This integrated evaluation model is called LOTUS which stands for LOCA Toolkit for the U.S. Light Water Reactors. The focus of LOTUS is to establish the automation of the interface between the five disciplines. LOTUS framework ensures that uncertainties can be easily propagated to solve the problems forward and backward (inverse) with large numbers of simulations. The objective of this work is to perform uncertainty quantification (UQ) and sensitivity analysis (SA) on a coupling of CASL's Virtual Environment for Reactor Applications-Core Simulator (VERA-CS) and NRC's fuel performance code FRAPCON. The work presented herein is a subset of the overall LOTUS framework. The figures of merit (FOM) selected for analysis are minimum departure from nucleate boiling ratio (MDNBR), maximum fuel centerline temperature (MFCT), and the gap conductance at peak power (GCPP). All FOM are obtained at the end of the of a depletion case for a half cycle lasting 263 days. The minimum is taken over all spatial location for each pin at the end of the half cycle. This value is calculated in VERACS alone and has been analyzed in previous work by Brown and Zhang. The MFCT is available in both VERA-CS and FRAPCON, but is calculated via different methods (see Model Inconsistencies section for more details). For this work, only the FRAPCON values are reported due to the code's maturity from extensive validation in fuel performance. Previous work includes SA on a standalone nuclear fuel performance code VTT-modified FRAPCON- 3.4 by Ikonen. The GCPP is the gap conductance at the location corresponding to maximum power in the fuel rod. This value is available in FRAPCON alone and highlights the capability of LOTUS to obtain unique fuel performance metrics for each pin in a core simulation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗