Engineering Papers⌕ Search

Engineering topics

Hamilton, Steven P

Publications and source records attributed to Hamilton, Steven P.

A Cast of Thousands: How the IDEAS Productivity Project Has Advanced Software Productivity and Sustainability

Computational and data-enabled science and engineering are revolutionizing advances throughout science and society, at all scales of computing. For example, teams in the U.S. Department of Energy’s Exascale Computing Project have been tackling new frontiers in modeling, simulation, and analysis by exploiting unprecedented exascale computing capabilities—building an advanced software ecosystem that supports next-generation applications and addresses disruptive changes in computer architectures. However, concerns are growing about the productivity of the developers of scientific software. Members of the Interoperable Design of Extreme-scale Application Software project serve as catalysts to address these challenges through fostering software communities, incubating and curating methodologies and resources, and disseminating knowledge to advance developer productivity and software sustainability. This article discusses how these synergistic activities are advancing scientific discovery—mitigating technical risks by building a firmer foundation for reproducible, sustainable science at all scales of computing, from laptops to clusters to exascale and beyond.

97 MATHEMATICS AND COMPUTING↗

Shift: A Massively Parallel Monte Carlo Radiation Transport Package

This paper discusses the massively-parallel Monte Carlo radiation transport package, Shift, developed at Oak Ridge National Laboratory. It reviews the capabilities, implementation, and parallel performance of this code package. Scaling results demonstrate very good strong and weak scaling behavior of the implemented algorithms. Benchmark results from various reactor problems show that Shift results compare well to other contemporary Monte Carlo codes and experimental results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Three-dimensional discrete ordinates reactor assembly calculations on GPUs

In this paper we describe and demonstrate a discrete ordinates sweep algorithm on GPUs. This sweep algorithm is nested within a multilevel comunication-based decomposition based on energy. We demonstrated the effectiveness of this algorithm on detailed three-dimensional critical experiments and PWR lattice problems. For these problems we show improvement factors of 4 6 over conventional communication-based, CPU-only sweeps. These sweep kernel speedups resulted in a factor of 2 total time-to-solution improvement.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Analysis of Anderson Acceleration on a Simplified Neutronics/Thermal Hydraulics System

ABSTRACT A standard method for solving coupled multiphysics problems in light water reactors is Picard iteration, which sequentially alternates between solving single physics applications. This solution approach is appealing due to simplicity of implementation and the ability to leverage existing software packages to accurately solve single physics applications. However, there are several drawbacks in the convergence behavior of this method; namely slow convergence and the necessity of heuristically chosen damping factors to achieve convergence in many cases. Anderson acceleration is a method that has been seen to be more robust and fast converging than Picard iteration for many problems, without significantly higher cost per iteration or complexity of implementation, though its effectiveness in the context of multiphysics coupling is not well explored. In this work, we develop a one-dimensional model simulating the coupling between the neutron distribution and fuel and coolant properties in a single fuel pin. We show that this model generally captures the convergence issues noted in Picard iterations which couple high-fidelity physics codes. We then use this model to gauge potential improvements with regard to rate of convergence and robustness from utilizing Anderson acceleration as an alternative to Picard iteration.

Toth, Alex↗

A Validation Study of Pin Heat Transfer for UO2 Fuel Based on the IFA-432 Experiments

The IFA-432 (Integrated Fuel Assessment) experiments from the International Fuel Performance Experiments (IFPE) database were designed to study the effects of gap size, fuel density, and fuel densification on fuel centerline temperature in light-water-reactor fuel. An evaluation of nuclear fuel pin heat transfer in the FRAPCON-3.4 and Exnihilo codes for uranium dioxide (UO$_2$) fuel systems was performed, with a focus on the densification stage (2.2 \unitfrac{GWd}{MT(UO$_{2}$)}). In addition, sensitivity studies were performed to evaluate the effect of the radial power shape and approximations to the geometry to account for the thermocouple hole. The analysis demonstrated excellent agreement for rods 1, 2, 3, and 5 (varying gap thicknesses and density with traditional fuel), demonstrating the accuracy of the codes and their underlying material models for traditional fuel. For rod 6, which contained unstable fuel that densified an order of magnitude more than traditional, stable fuel, the magnitude of densification was over-predicted and the temperatures were outside of the experimental uncertainty. The radial power shape within the fuel was shown to significantly impact the predicted centerline temperatures, whereas modeling the fuel at the thermocouple location as either annular or solid was relatively negligible. This has provided an initial benchmarking of the pin heat transfer capability of Exnihilo for UO$_2$ fuel with respect to a well-validated nuclear fuel performance code.

Phillippe, Aaron M↗

A Validation Study of Pin Heat Transfer for MOX Fuel Based on the IFA-597 Experiments

Abstract The IFA-597 (Integrated Fuel Assessment) experiments from the International Fuel Performance Experiments (IFPE) database were designed to study the thermal behavior of mixed oxide (MOX) fuel and the effects of an annulus on fission gas release in light-water-reactor fuel. An evaluation of nuclear fuel pin heat transfer in the FRAPCON-3.4 and Exnihilo codes for MOX fuel systems was performed, with a focus on the first 20 time steps ( 6 GWd/MT(iHM)) for explicit comparison between the codes. In addition, sensitivity studies were performed to evaluate the effect of the radial power shape and approximations to the geometry to account for the thermocouple hole, dish, and chamfer. The analysis demonstrated relative agreement for both solid (rod 1) and annular (rod 2) fuel in the experiment, demonstrating the accuracy of the codes and their underlying material models for MOX fuel, while also revealing a small energy loss artifact in how gap conductance is currently handled in Exnihilo for chamfered fuel pellets. The within-pellet power shape was shown to significantly impact the predicted centerline temperatures. This has provided an initial benchmarking of the pin heat transfer capability of Exnihilo for MOX fuel with respect to a well-validated nuclear fuel performance code.

Phillippe, Aaron M↗

Watts Bar Unit 1 Cycle Zero Power Physics Tests Analysis with VERA-CS

The Consortium for Advanced Simulation of Light Water Reactors (CASL) is developing a collection of methods and software products known as VERA, the Virtual Environment for Reactor Applications, including a core simulation capability called VERA-CS. A key milestone for this endeavor is to validate VERA against measurements from operating nuclear power reactors. The first step in validation against plant data is to determine the ability of VERA to accurately simulate the initial startup physics tests for Watts Bar Nuclear Power Station, Unit 1 (WBN1) cycle 1. VERA-CS calculations were performed with the Insilico code developed at ORNL using cross section processing from the SCALE system and the transport capabilities within the Denovo transport code using the SPN method. The calculations were performed with ENDF/B-VII.0 cross sections in 252 groups (collapsed to 23 groups for the 3D transport solution). The key results of the comparison of calculations with measurements include initial criticality, control rod worth critical configurations, control rod worth, differential boron worth, and isothermal temperature reactivity coefficient (ITC). The VERA results for these parameters show good agreement with measurements, with the exception of the ITC, which requires additional investigation. Results are also compared to those obtained with Monte Carlo methods and a current industry core simulator.

CASL↗

Multiphysics Simulations for LWR Analysis

Accurate prediction of the neutron and temperature distributions within an operating nuclear reactor requires the solution of multiple coupled physics equations. In a light water reactor (LWR), there is a very strong coupling between the power distribution (governed by the radiation transport equation) and the temperature and density distributions (governed by a thermal diusion equation in combination with a fluid flow model). This study aims to begin to quantify the impact of such feedback mechanisms as well as identify numerical diculties associated with such multiphysics problems. A description of the multiphysics model and current solution strategy within the Exnihilo code package will be presented, along with numerical results for a representative 3 3 PWR mini-assembly.

multi-physics↗