Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core simulation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Determination of the spent fuel decay heat with the VERA core simulator

VERA core simulator capabilities to predict spent fuel decay heat have been explored. The comparison with the Serpent2 and NEWT/TRITON sequence from the SCALE package, performed on the NPP Krško fuel test cases, has enabled independent verification of the successful implementation of the VERA depletion/decay calculation sequence applied on the extended nuclide set. Since the VERA is a 3-D core simulator, the focus in the rest of the paper was on the investigation of the 3-D effects and non-linearity of the soluble boron concentration averaging process, which are difficult to assess with most of available computational tools. A difference in the decay heat prediction between the 3-D and 2-D model approaches was analyzed. In addition, the impact of the Inconel grids on the decay heat was determined. Averaging the concentration of soluble boron over the fuel irradiation period has been found to be a viable approach, since it results in slight decay heat overprediction, which is considered conservative for most applications. Finally, demonstration on the Watts Bar Unit 1 cycle 1 depletion has shown formidable VERA capabilities to accurately predict spent fuel decay heat.

Kromar, Marjan↗
DOE OSTI2012

A Two-Step Approach to Uncertainty Quantification of Core Simulators

For the multiple sources of error introduced into the standard computational regime for simulating reactor cores, rigorous uncertainty analysis methods are available primarily to quantify the effects of cross section uncertainties. Two methods for propagating cross section uncertainties through core simulators are the XSUSA statistical approach and the “two-step” method. The XSUSA approach, which is based on the SUSA code package, is fundamentally a stochastic sampling method. Alternatively, the two-step method utilizes generalized perturbation theory in the first step and stochastic sampling in the second step. The consistency of these two methods in quantifying uncertainties in the multiplication factor and in the core power distribution was examined in the framework of phase I-3 of the OECD Uncertainty Analysis in Modeling benchmark. With the Three Mile Island Unit 1 core as a base model for analysis, the XSUSA and two-step methods were applied with certain limitations, and the results were compared to those produced by other stochastic sampling-based codes. Based on the uncertainty analysis results, conclusions were drawn as to the method that is currently more viable for computing uncertainties in burnup and transient calculations.

Yankov, Artem↗
DOE OSTI2014

An approach for coupled-code multiphysics core simulations from a common input

This study describes an approach for coupled-code multiphysics reactor core simulations that is being developed by the Virtual Environment for Reactor Applications (VERA) project in the Consortium for Advanced Simulation of Light-Water Reactors (CASL). In this approach a user creates a single problem description, called the “VERAIn” common input file, to define and setup the desired coupled-code reactor core simulation. A preprocessing step accepts the VERAIn file and generates a set of fully consistent input files for the different physics codes being coupled. The problem is then solved using a single-executable coupled-code simulation tool applicable to the problem, which is built using VERA infrastructure software tools and the set of physics codes required for the problem of interest. The approach is demonstrated by performing an eigenvalue and power distribution calculation of a typical three-dimensional 17 × 17 assembly with thermal–hydraulic and fuel temperature feedback. All neutronics aspects of the problem (cross-section calculation, neutron transport, power release) are solved using the Insilico code suite and are fully coupled to a thermal–hydraulic analysis calculated by the Cobra-TF (CTF) code. The single-executable coupled-code (Insilico-CTF) simulation tool is created using several VERA tools, including LIME (Lightweight Integrating Multiphysics Environment for coupling codes), DTK (Data Transfer Kit), Trilinos, and TriBITS. Parallel calculations are performed on the Titan supercomputer at Oak Ridge National Laboratory using 1156 cores, and a synopsis of the solution results and code performance is presented. Finally, ongoing development of this approach is also briefly described.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Demonstrate VERA Core Simulator Performance Improvements L2:PHI.P13.03

This report describes the performance improvements made to the VERA Core Simulator (VERA-CS) during FY2016. The development of the VERA Core Simulator has focused on the capability needed to deplete physical reactors and help solve various problems; this capability required the accurate simulation of many operating cycles of a nuclear power plant. The first section of this report introduces two test problems used to assess the run-time performance of VERA-CS using a source dated February 2016. The next section provides a brief overview of the major modifications made to decrease the computational cost. Following the descriptions of the major improvements, the run-time for each improvement is shown. Conclusions on the work are presented, and further follow-on performance improvements are suggested.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Best estimate plus uncertainty analysis of departure from nucleate boiling limiting case with CASL core simulator VERA-CS in response to PWR main steam line break event

VERA-CS (Virtual Environment for Reactor Applications, Core Simulator) is a coupled neutron transport and thermal-hydraulics subchannel code under development by the Consortium for Advanced Simulation of Light Water Reactors (CASL). VERA-CS was used to simulate a typical pressurized water reactor (PWR) full core response with 17x17 fuel assemblies for a main steam line break (MSLB) accident scenario with the most reactive rod cluster control assembly stuck out of the core. The accident scenario was initiated at the hot zero power (HZP) at the end of the first fuel cycle with return to power state points that were determined by a system analysis code and the most limiting state point was chosen for core analysis. The best estimate plus uncertainty (BEPU) analysis method was applied using Wilks’ nonparametric statistical approach. In this way, 59 full core simulations were performed to provide the minimum departure from nucleate boiling ratio (MDNBR) at the 95/95 (95% probability with 95% confidence level) tolerance limit. The results show that this typical PWR core remains within MDNBR safety limits for the MSLB accident.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2012

Implementation and validation of the variational nodal expansion method core simulator

The variational nodal expansion method (VNEM) has been developed by the OECD Halden Reactor Project for the purpose of implementation into a full-core simulator. This version of the VNEM does not make use of the diffusion approximation, but instead provides a fast converged solution to the even-parity neutron transport equation in several minutes using nodal methods and variational techniques. This 'method to simulator' implementation has been completed and the VNEM-simulator is currently undergoing continued validation efforts to ensure that it can track more nuclear reactor cycles effectively. This paper presents the results of the most recent development on the 3D VNEM core simulator implementation and validation efforts. The method is briefly described and data is presented and discussed for three cases for the initial criticality tests at the beginning of life, hot stand-by state for the Ringhals-3 PWR. Also presented here are the results for one year of operation of the maiden cycle of this unit. The effective neutron multiplication factor, k{sub eff}, and relative neutron detector readings are compared showing excellent agreement. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Evaporation Of Hanford Waste Treatment Plant Direct Feed Low Activity Waste Effluent Management Facility Core Simulant

The Hanford Waste Treatment and Immobilization Plant (WTP) Low Activity Waste (LAW) vitrification facility will generate an aqueous condensate recycle stream (LAW Melter Off-Gas Condensate, LMOGC) from the off-gas system. The baseline plan for disposition of this stream during full WTP operations is to send it to the WTP Pretreatment Facility, where it will be blended with LAW, concentrated by evaporation, and recycled to the LAW vitrification facility. However, during the Direct Feed LAW (DFLAW) scenario, planned disposition of this stream is to evaporate it in a new evaporator, in the Effluent Management Facility (EMF), and then return it to the LAW melter. It is important to understand the composition of the effluents from the melter and new evaporator, so that the disposition of these streams can be accurately planned and accommodated. Furthermore, alternate disposition of the LMOGC stream would eliminate recycling of problematic components, and would reduce the need for closely integrated operation of the LAW melter and the Pretreatment Facilities. Long-term implementation of this option after WTP start-up would decrease the LAW vitrification mission duration and quantity of glass waste, amongst the other operational complexities such a recycle stream presents. In order to accurately plan for the disposition path, it is key to experimentally determine the fate of contaminants. To do this, testing is needed to accurately account for the buffering chemistry of the components, determine the achievable evaporation end point, identify insoluble solids that form, and determine the distribution of key regulatory-impacting constituents. The LAW Melter Off-Gas Condensate stream will contain components that are volatile at melter temperatures, have limited solubility in the glass waste form, and represent a materials corrosion concern, such as halides and sulfate. Because this stream will recycle within WTP, these components will accumulate in the Melter Condensate stream, exacerbating their impact on the number of LAW glass containers that must be produced. Diverting the stream reduces the halides and sulfates in the recycled Condensate and is a key outcome of this work. This overall program examines the potential treatment and immobilization of this stream to enable alternative disposal. The objective of this task was to demonstrate evaporation of a simulant of the LAW Melter Off-gas Condensate expected during DFLAW operations, in order to predict the composition of the effluents from the EMF evaporator to aid in planning for their disposition. This document describes the results of that test using the core simulant. This simulant formulation is designated as the “core simulant”; other additives will be included for specific testing, such as volatiles for evaporation or hazardous metals for measuring leaching properties of waste forms. The results indicate that the simulant can easily be concentrated via evaporation. During that the pH adjustment step in simulant preparation, ammonium is quickly converted to ammonia, and most of the ammonia was stripped from the simulated waste and partitioned to the condensate. Additionally, it was found that after concentrating (>12x) and cooling that a small amount of LiF and Na 3 (SO 4 )F precipitate out of solution. With the exception of ammonia, analysis of the condensate indicated very low to below detectable levels of many of the constituents in the simulant, yielding very high decontamination factors (DF).

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

Nuclear data uncertainty propagation and modeling uncertainty impact evaluation in neutronics core simulation

Uncertainty analysis is a critical requirement in reactor simulation as it is used to quantify the reliability of best-estimate calculation. A comprehensive uncertainty analysis should characterize all sources of uncertainties in a computationally-feasible and scientifically-defendable manner. Here we employ a well-established reduced order modeling (ROM) based uncertainty quantification methodology to propagate uncertainties throughout neutronic calculations. ROM relies on recent advances in randomized data mining techniques applied to large data streams. In our proposed implementation, the nuclear data uncertainties are first propagated from multi-group level through lattice physics calculation to generate few-group parameter uncertainties, described using a vector of mean values and a covariance matrix. Employing an ROM-based compression of the covariance matrix, the few-group uncertainties are then propagated through downstream core simulation in a computationally efficient manner. This straightforward approach, albeit efficient as compared to brute force forward and/or adjoint-based methods, often employs a number of assumptions that have been unquestioned in the literature of neutronic uncertainty analysis. This manuscript argues that these assumptions could introduce another source of uncertainty referred to as modeling uncertainties, whose magnitude needs to be quantified in tandem with nuclear data uncertainties. Thus, our primary goal is to explore the interactions between these two uncertainty sources in order to assess whether modeling uncertainties have an impact on parameter uncertainties. To explore this endeavor, the impact of a number of modeling assumptions on core attributes uncertainties is quantified. The study employs a CANDU reactor model, with Serpent and NEWT as lattice physics solvers and NESTLE-C as core simulator. The modeling assumptions investigated include those related with the uncertainty propagation method employed, e.g., deterministic vs. stochastic, the few-group energy structure employed to represent the cross-sections, the resonance treatment in lattice physics calculation, the reference values for the cross-section, and the number of samples employed to render ROM compression. Results indicate that some of the modeling assumptions could have a non-negligible impact on the core responses propagated uncertainties, highlighting the need for a more comprehensive approach to combine parameter and modeling uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Summary of VERA Core Simulator Performance Improvements

This report describes the performance improvements made to the VERA Core Simulator (VERA CS) during FY2016. At the beginning of the year, a timing study was performed to identify areas of the code where the run-times could be improved. This study identified six broad areas for improvement: • Cross section lookup and CMFD stabilization, • CTF parallel partitioning, • MOC methodology, • Subgroup methodology, • MOC parallel partitioning, and • CMFD Krylov solution techniques. The goal of this work was to reduce the VERA-CS run-times so a complete cycle depletion could be run with 1,000 cores overnight (i.e., the “1,000 core depletion” criterion). This goal was set by our industry partners so they could achieve reasonable turnaround on medium-sized computer clusters. The results presented in this report show a code speed-up by a factor of 5, which successfully satisfies the goal set by industry.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2012

ABWR start-up test analysis using BWR core simulator with three-dimensional direct response matrix method

The ABWR start-up test analysis has been done with the BWR core simulator using the three--dimensional direct response matrix (3D-DRM) method. The Monte Carlo code VMONT made the sub-response matrices for the 3D-DRM method. Each boundary surface was subdivided by 4 x 4 for transverse segments, by 4 for angular segments and by 4 for axial zones in a node. For the calculation speedup, the 3D-DRM code used the divided sub-response matrices data set. The code used the MPI and OpenMP for the parallelized method. The median value is set as the average critical eigenvalues. The changes from the maximum value to the minimum value are 0.34 %{Delta}k with the spectral history method and 0.40 %{Delta}k without it, and the respective standard deviations were 0.12 % and 0.14 %. Using the spectral history method decreased the variation by 0.06 %{Delta}k. The root mean square differences of the axial power distribution were about 6 % between the analysis results and the plant data. Using the currents which converged in the previous exposure step reduced the number of iterations when the CR pattern changed only slightly. The averaged calculation time for each exposure step was about 5 hours on 12 PC Linux cluster servers with Core 2 Quad 3 GHz. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

VERA Core Simulator methodology for pressurized water reactor cycle depletion

This paper describes the methodology developed and implemented in the Virtual Environment for Reactor Applications Core Simulator (VERA-CS) to perform high-fidelity, pressurized water reactor (PWR), multicycle, core physics calculations. Depletion of the core with pin-resolved power and nuclide detail is a significant advance in the state of the art for reactor analysis, providing the level of detail necessary to address the problems of the U.S. Department of Energy Nuclear Reactor Simulation Hub, the Consortium for Advanced Simulation of Light Water Reactors (CASL). VERA-CS has three main components: the neutronics solver MPACT, the thermal-hydraulic (T-H) solver COBRA-TF (CTF), and the nuclide transmutation solver ORIGEN. This paper focuses on MPACT and provides an overview of the resonance self-shielding methods, macroscopic-cross-section calculation, two-dimensional/one-dimensional (2-D/1-D) transport, nuclide depletion, T-H feedback, and other supporting methods representing a minimal set of the capabilities needed to simulate high-fidelity models of a commercial nuclear reactor. Results are presented from the simulation of a model of the first cycle of Watts Bar Unit 1. The simulation is within 16 parts per million boron (ppmB) reactivity for all state points compared to cycle measurements, with an average reactivity bias of <5 ppmB for the entire cycle. Comparisons to cycle 1 flux map data are also provided, and the average 2-D root-mean-square (rms) error during cycle 1 is 1.07%. To demonstrate the multicycle capability, a state point at beginning of cycle (BOC) 2 was also simulated and compared to plant data. The comparison of the cycle 2 BOC state has a reactivity difference of +3 ppmB from measurement, and the 2-D rms of the comparison in the flux maps is 1.77%. Lastly, these results provide confidence in VERA-CS’s capability to perform high-fidelity calculations for practical PWR reactor problems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Projected changes to severe thunderstorm environments as a result of twenty-first century warming from RegCM CORDEX-CORE simulations

Hazardous weather related to the occurrence of severe thunderstorms including tornadoes, high-winds, and hail cause significant damage globally to life and property every year. Yet the impact on these storms from a warming climate remains a difficult task due to their transient nature. Here, we investigate the change in large-scale environments in which severe thunderstorms form during twenty-first century warming (RCP2.6 and RCP8.5) in a group of RegCM CORDEX-CORE simulations. Severe potential is measured in terms of Convective Available Potential Energy (CAPE) and vertical wind-shear during the severe seasons in three regions which are known to currently be prone to severe hazards: North America, subtropical South America, and eastern India and Bangladesh. In every region, environments supportive for severe thunderstorms are projected to increase during the warm season months in both the RCP2.6 and RCP8.5 scenarios during the twenty-first century. The number of days supportive for severe thunderstorms increases by several days per season over the vast majority of each region by the end of the century. Analyzing the CAPE and shear trends during the twenty-first century, we find seasonally and regionally specific changes driving the increased severe potential. Twenty-first century surface warming is clearly driving a robust increase in CAPE in all regions, however poleward displacement of vertical shear in the future leads to the displacement of severe environments over North America and South America. The results found here relate that severe impacts in the future cannot be generalized globally, and that regionally specific changes in vertical shear may drive future movement of regions prone to severe weather.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Development and Test of a New Verification Scheme for Transient Core Simulators

Transient calculations in commercial nuclear reactors are performed while typically relying on a time-dependent neutron transport solver or a low-order solver (i.e. diffusion). In order to be licensed, the codes used by the industry need to go through a process of verification and validation, with the verification carried out by comparing the results of simulations to analytical or semi-analytical solutions. Such analytical or semi-analytical solutions can only be obtained if the system to be modelled during the verification process is either fully homogeneous or piecewise homogeneous. This paper reports on the development of a different verification approach that can be applied to fully heterogeneous systems. It relies on the extraction of the point-kinetic response of the reactor (which can be estimated from the results of core simulations) and on its subsequent comparison with its expected analytical form.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Uncertainty quantification and sensitivity analysis with CASL Core Simulator VERA-CS

Uncertainty quantification and sensitivity analysis are important for nuclear reactor safety design and analysis. A 2x2 fuel assembly core design was developed and simulated by the Virtual Environment for Reactor Applications, Core Simulator (VERA-CS) coupled neutronics and thermal-hydraulics code under development by the Consortium for Advanced Simulation of Light Water Reactors (CASL). An approach to uncertainty quantification and sensitivity analysis with VERA-CS was developed and a new toolkit was created to perform uncertainty quantification and sensitivity analysis with fourteen uncertain input parameters. Furthermore, the minimum departure from nucleate boiling ratio (MDNBR), maximum fuel center-line temperature, and maximum outer clad surface temperature were chosen as the selected figures of merit. Pearson, Spearman, and partial correlation coefficients were considered for all of the figures of merit in sensitivity analysis and coolant inlet temperature was consistently the most influential parameter. We used parameters as inputs to the critical heat flux calculation with the W-3 correlation were shown to be the most influential on the MDNBR, maximum fuel center-line temperature, and maximum outer clad surface temperature.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

VERA Core Simulator Methodology for PWR Cycle Depletion

This paper describes the methodology developed and implemented in MPACT for performing high-fidelity pressurized water reactor (PWR) multi-cycle core physics calculations. MPACT is being developed primarily for application within the Consortium for the Advanced Simulation of Light Water Reactors (CASL) as one of the main components of the VERA Core Simulator, the others being COBRA-TF and ORIGEN. The methods summarized in this paper include a methodology for performing resonance self-shielding and computing macroscopic cross sections, 2-D/1-D transport, nuclide depletion, thermal-hydraulic feedback, and other supporting methods. These methods represent a minimal set needed to simulate high-fidelity models of a realistic nuclear reactor. Results demonstrating this are presented from the simulation of a realistic model of the first cycle of Watts Bar Unit 1. The simulation, which approximates the cycle operation, is observed to be within 50 ppm boron (ppmB) reactivity for all simulated points in the cycle and approximately 15 ppmB for a consistent statepoint. The verification and validation of the PWR cycle depletion capability in MPACT is the focus of two companion papers.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Pronghorn: A Porous Media Thermal-Hydraulics Core Simulator and its Validation with the SANA Experiments

Pebble bed High Temperature Reactors (HTRs) are characterized by many advantageous design features, such as excellent passive heat removal in accidents, large margins to fuel failure, and online refueling potential. However, a significant challenge in the core modeling of pebble bed reactors is the complex fuel-coolant structure. This paper presents a new porous media simulation code, Pronghorn, that aims to alleviate modeling challenges for pebble bed reactors by providing a fast-running, medium- fidelity core simulator. Pronghorn is intended to accelerate the design and analysis cycle for pebble bed and prismatic HTRs by permitting fast scoping studies and providing boundary conditions for systems-level analysis. Pronghorn is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) using modern software practices and a thorough testing framework. This paper describes the physical models used in Pronghorn and demonstrates Pronghorn’s capability for modeling gas-cooled pebble bed HTRs by presenting simulation results obtained for the German SANA pebble bed decay heat experiments. Within the limitations of the porous media approximation and existing available closure relationships, Pronghorn predicts the SANA experimental pebble temperatures well, expanding the code’s validation base. A brief code-to-code comparison shows a similar level of accuracy with other porous media simulation tools; this enhances Pronghorn’s advantages of an arbitrary equation of state, unstructured mesh capabilities, compressible flow models, high coupling potential to MOOSE fuels performance and systems-level thermal-hydraulics codes, and modern software design.

42 ENGINEERING↗
DOE OSTI2016

Computational and Experimental Examination of Simulated Core Damage and Relocation Dynamics of a BWR Fuel Assembly - 16013

The Japan Atomic Energy Agency (JAEA) plans a large-scale test to evaluate damage and relocation behavior of Boiling Water Reactor (BWR) core materials consisting of fuel rods (zirconia as fuel simulant), channel boxes, control blade assemblies as well as lower core-support structures. Its purpose is to contribute to the understanding of the core material relocation behavior in the event of severe accidents with the BWR design conditions for which existing experimental database is quite limited. Prior to large-scale testing anticipated at JAEA in the future, JAEA desires preliminary investigations to examine melting test pieces, being composed of a mixture of zirconia and metals by a non-transferred thermal plasma heating system. The purpose of such tests is to verify the materials and test piece will be heated by plasma as a primary heat source to the target temperature of ca. 2900 K and to collect data about the relocation behavior of the high-temperature materials. This paper describes the first in a series of tests demonstrating the effective use of thermal plasma technology to recreate the dynamics of a typical fuel core bundle channel assembly section. Results from preliminary computational simulations of simplified fuel bundle assemblies are presented illustrating the effectiveness of a 150 kW, non transferred (NTR) configuration of a thermal plasma argon jet to heat the assembly within a highly-insulated containment vessel. An experimental test program using the computational analyses as a basis and a Phoenix Solutions Co (PSC) commercial model PT150, 150 kW plasma torch is described and the test results are presented illustrating the time versus temperature history of the fuel bundle assembly melting dynamics. Details on plasma operation are discussed including: power requirements, current flow, voltage levels achieved, argon gas flow rates, processing times and standoff distances of the torch discharge from the top of the fuel bundle assembly (achieved with remote control of the axial (vertical) translation of the torch within the containment processing chamber). Photos taken from the top of the processing chamber are presented illustrating the various phases of melting characteristics of the zircaloy-clad, zirconium oxide filled fuel tubes, all within a stainless steel tube bundle channel centered within an insulated zirconia crucible. Based on lessons learned from this initial test campaign, a brief discussion of the requirements and features of a follow-on Phase II test program is provided. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗