Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reactor Physics Method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Benchmark Analysis of the HTR-10 with the MAMMOTH Reactor Physics Application

This report documents a set of benchmarks used to validate the Serpent and MAMMOTH neutron transport models of the HTR-10 reactor developed at Idaho National Laboratory. The high-fidelity Serpent Monte Carlo models of the HTR-10 critical and full core configurations include both random, discrete distributions of TRISO particles in the pebbles and random, discrete pebble distributions in the pebble bed core. The Serpent results agree very well with the critical and control rod worth measurements provided in the International Re- actor Physics Experiment Evaluation Project (IRPhEP) report if the ENDF/B-VII.r1 data- set is used. These Serpent models are subsequently used to prepare cross sections and flux tallies in a 10 coarse energy group structure for the MAMMOTH Reactor Physics MOOSE- based application. We find that MAMMOTH can reproduce the Monte Carlo solution for the HTR-10 reactor by using the PJFNK-SPH equivalence method. In all cases studied, the MAMMOTH results are within 120 pcm of the Monte Carlo reference calculation and the maximum errors in neutron absorption and generation rates are within 0.536% and 0.215% from the reference calculation, respectively. A set of benchmark exercises from the International Atomic Energy Agency (IAEA) HTR-10 benchmark in IAEA-TECDOC-1382 were conducted including parts B1 (critical), B2 (temperature coefficient) and B3 (control rod worth). The results we obtain for the critical core are in excellent agreement with the expected value and we confirm that with ENDF/B-VII.r1 data we match the critical eigen- value keff=1 to within 100 pcm. Our temperature coefficient calculation is consistent with the VSOP results generated by the German participants. The Chinese and South African VSOP models produce lower temperature coefficients, a fact that was noted in the original IAEA benchmark, but for which no explanation was provided. Both a full and single control rod worth calculations with Serpent and MAMMOTH agree well with the results from the IAEA benchmark. Finally, and most importantly, this work demonstrates that MAM- MOTH can recover the Monte Carlo high-fidelity solutions to an excellent accuracy using the PJFNK-SPH method. This sets the stage for the ultimate goal of using MAMMOTH in transient, multi-physics scenarios where correction of the cross sections is obtained for steady-state conditions and then used in the transient scenario for which no Monte-Carlo solution is available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Evaluation of Oak Ridge National Laboratory Health Physics Research Reactor Operation Data for Critical Benchmark Creation

The Oak Ridge National Laboratory (ORNL) Health Physics Research Reactor (HPRR) was a research reactor designed and built at ORNL in 1961. The critical assembly used a highly enriched uranium and molybdenum alloy as the fuel and could be operated in steady-state or burst modes. The HPRR has recently been the object of an investigation to create a criticality benchmark. Such benchmarks are very important, as they are used primarily to show the accuracy of newly developed modeling codes and to help experimental validation and reactor licensing. The evaluated experiments considered in this paper were carried out between 1974 and 1986 from various HPRR activities such as steady-state subcritical, steady-state critical, and burst prompt super-critical operations of the reactor for dosimetry, irradiation, or training purposes. By using the HPRR experimental logbook information and the as-built drawings of the critical assembly, a highly detailed model of the HPRR was created with SCALE 6.2.4/KENO-VI, and a first version of a critical benchmark of the HPRR was developed following the International Criticality Safety Benchmark Evaluation Project (ICSBEP) guidelines for thorough description and uncertainty/sensitivity quantification. Unfortunately, in most of the evaluated experiments, the obtained difference between calculated and experimental k eff is around 1,000 pcm, corresponding to a relative error of approximately 1%, beyond the quality standards of the ICSBEP recommending a relative error below 0.1%. Moreover, the derived experimental uncertainty is high, around 4% relative, mainly due to the U-Mo fuel density uncertainty, but also from numerous other factors. For these reasons, the creation of a valuable critical benchmark from HPRR operation data is thus far compromised. In this paper, the different steps of the experiments’ evaluation are summarized, and the reasons for the experimental/calculation discrepancies and potential ways to solve them are explored. This paper also aims to remind us always to exercise considerable care when performing experimental work, and to record all the data possible for potential future uses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

NDCS Group Meeting Work overview: NCSP Task ORNL IP&D-5 [Slides]

The presentation discusses the primary goal of the NCSP Task IP&D-5 which is to use available data from HPRR operation to create a benchmark report for inclusion in the ICSBEP, as a Criticality Accident Alarm System (CAAS) shielding benchmark. The presentation elaborates on the HPRR operation and how it relates to the KENO VI model, sulfur pellets measurements, and calculation methodology between KENO and MAVRIC FW-CADIS. The presentation concludes by stating that some work is still needed to evaluate uncertainties and to explain the discrepancy computation/experiment.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Criticality Experiments to Reduce Compensating Errors in Plutonium Nuclear Data

Compensating errors between nuclear data observables in a library can adversely impact application simulations. The primary goal of the EUCLID project (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) is to reduce compensating errors in nuclear data. A new criticality experiment, described in this work, was designed with the specific target nuclear data of 239 Pu fission, inelastic scattering, elastic scattering, capture, nu-bar, and prompt fission neutron spectrum (PFNS). This work will focus on the design and execution of the EUCLID experiment, performed on the Planet vertical lift critical assembly machine at the National Criticality Experiments Research Center (NCERC). The criticality experiment includes two different configurations with very different geometries: one is cube-like to minimize neutron leakage while the other is slab-like to maximize leakage. Having these two widely varying configurations allows the scattering sensitivities of 239 Pu to the neutron multiplication factor to be greatly changed while minimally impacting the other cross section sensitivities. Both configurations utilize the Pu ZPPR (Zero Power Physics Reactor) plates as fuel. The experiments were designed using a D-Optimality criteria, which is an optimization method minimizing the log-determinant of the adjusted nuclear data covariance for the target reactions. These experiments include not only inference of k eff , as done in all critical benchmark experiments, but several other responses as well, such as neutron multiplication measurements and reaction rate ratios. After analysis of the measured data is complete, adjustment of nuclear data will be performed to assess whether the new experimental data successfully reduced compensating errors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Ceramic nuclear fuel fracture modeling with the extended finite element method

Ceramic fuel pellets used in nuclear light water reactors experience significant fracture due to the high thermal gradients experienced under normal operating conditions. This has important effects on the performance of the fuel system. Because of this, a realistic, physically based fracture modeling capability is essential to predict fuel behavior in a wide variety of normal and off-normal conditions. The extended finite element method (X-FEM) is a powerful method to represent arbitrary propagating discrete cracks in finite element models, and has many characteristics that make it attractive for nuclear fuel performance analysis. This paper describes the implementation of X-FEM in a multiphysics fuel performance code and presents applications of that capability. These applications include several thermal mechanics fracture benchmark problems, which demonstrate the accuracy of this approach. It also includes application of this capability to study nuclear fuel fracture, both on stationary and propagating cracks. The study on stationary cracks shows the effects that varying the number of radial cracks and the length of those cracks have on the energy release rate. Here, the propagating crack case demonstrates random initiation and subsequent propagation of interacting thermally induced cracks during an initial ramp to full power with fresh fuel.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Variance Estimation in Monte Carlo Eigenvalue Simulations Using Spectral Analysis Method

The Monte Carlo (MC) method is the 'gold standard' for performing nuclear reactor calculations and criticality safety analyses because of its ability to handle complex geometries and physics. Reactor and criticality safety analyses require solving eigenvalue problems in complex nuclear systems. Once a stationary fission source distribution is obtained in MC simulations, the sample mean of many stationary cycles is calculated. The MC simulation sample mean is typically referred to as a tally. Variance, or standard deviation, of the sample mean/tally is needed to indicate the level of statistical uncertainty in the simulation and to understand the convergence of the sample mean. Current MC codes typically use sample variance to estimate the statistical uncertainty of the simulation and assume that the MC stationary cycles are independent. However, it is known that cycles are correlated, and estimators of the variance that ignore these correlations systematically underestimate the variance. In this paper, we propose a novel method in the Fourier frequency domain to estimate the variance of the sample mean in MC eigenvalue calculations. Results indicate that the proposed estimator can be unbiased compared to the real variance. The proposed estimator requires minimal additional computing effort.

42 ENGINEERING↗
DOE OSTI2025

CONFLUX: A standardized framework to calculate reactor antineutrino flux

Nuclear fission reactors are abundant sources of antineutrinos for neutrino physics experiments. The flux and spectrum of antineutrinos emitted by a reactor can indicate its activity and composition, suggesting potential applications of neutrino measurements beyond fundamental scientific studies that may be valuable to society. The utility of reactor antineutrinos for applications and fundamental science is dependent on the availability of precise predictions of these emissions. For example, in the last decade, disagreements between reactor antineutrino measurements and models have inspired revision of reactor antineutrino calculations and standard nuclear databases as well as searches for new fundamental particles not predicted by the Standard Model of particle physics. Past predictions and descriptions of the methods used to generate them are documented to varying degrees in the literature, with different modeling teams incorporating a range of methods, input data, and assumptions. The resulting difficulty in accessing or reproducing past models and reconciling results from differing approaches complicates the future study and application of reactor antineutrinos. The CONFLUX (Calculation Of Neutrino FLUX) software framework is a neutrino prediction tool built with the goal of simplifying, standardizing, and democratizing the process of reactor antineutrino flux calculations. CONFLUX includes three primary methods for calculating the antineutrino emissions of nuclear reactors or individual beta decays that incorporate common nuclear data and beta decay theory. The software is prepackaged with the current nuclear databases, including ENDF.B/VIII, JEFF-3.3, and ENSDF, and it includes the capability to predict time-dependent reactor emissions, adjust nuclear database or beta decay inputs/assumptions, and propagate related sources of uncertainty. Here, this paper describes the CONFLUX software structure, details the methods used for flux and spectrum calculations, and provides examples of potential use cases.

Zhang, Xianyi [Lawrence Livermore National Laborat↗
DOE OSTI2017

VUSAT - A VERA-CS Uncertainty and Sensitivity Analysis Toolkit

VERA-CS (the Virtual Environment for Reactor Applications, Core Simulator) is a high-fidelity multiphysics core simulator under development by the Consortium for Advanced Simulation of Light Water Reactors (CASL). VERA-CS includes three main components: MPACT for reactor physics and neutron transport, COBRA-TF for thermal-hydraulics, and ORIGEN for isotopic depletion. VERA-CS has been applied to an array of problems including core physics analysis, and the modeling for all of the fuel cycles of Watts Bar Nuclear Unit 1 (WBN1). MPACT is a 3D pin-resolved reactor transport code developed by Oak Ridge National Laboratory (ORNL) and the University of Michigan. COBRA-TF (CTF) is a three-field two-phase flow subchannel analysis code developed by ORNL and North Carolina State University. CTF solves the mass, momentum, and energy equations for the liquid, vapor, and droplet phases for each subchannel in the core and captures the axial flow within each channel as well as cross-flow between channels. Direct coupling between neutronics (MPACT) and thermal-hydraulics (CTF) in VERA-CS provides temperature and density feedback. CTF includes a number of models important for LWR safety and design analysis including flow regime dependent two-phase wall heat transfer, inter-phase heat transfer and drag, droplet breakup, and quench-front tracking. ORIGEN is developed by ORNL and released through the SCALE package. ORIGEN is directly integrated in MPACT through its API to calculate the changes in nuclides for every depletable cross-section region. Uncertainties exist in safety calculations with computer codes for nuclear reactor modeling and simulation, and they come from code limitations, scaling inaccuracies embedded in the experimental data, and uncertainties associated with the state of the reactor at the initiation of the transient. Consequently, for applications of high fidelity calculations using VERA-CS, proper uncertainty quantification (UQ) and sensitivity analysis (SA) of the computing models of reactor neutronics and thermal-hydraulics are important to comprehensively understand the results. New methods in UQ and SA have to be developed in multi-physics calculations using VERA-CS. Toward this end, a software package called VUSAT which stands for VERA-CS Uncertainty and Sensitivity Analysis Toolkit has been developed to facilitate the applications of VERA-CS within the context of best estimate plus uncertainty (BEPU) analysis framework. In this work, the VUSAT application to analyze departure from nucleate boiling (DNB) is presented.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Utilizing Physics-Informed Synthetic Data to Train a Digital Twin for Predicting Reactor Operations

Understanding techniques to strengthen the nuclear nonproliferation regime is crucial in reducing the creation of nuclear weaponry on the basis of advancements in the nuclear energy industry. Prior to construction of a nuclear power plant, it is necessary to understand the proliferation potential of the plant’s reactor. Digital twins serve as a unique solution to recognizing reactor behavior indicative of nuclear proliferation. The following research conducted serves as a validation to training a digital twin on synthetic data fabricated via means of reactor physics simulations based on parameters of Idaho State University’s AGN-201 reactor. The synthetic data is utilized to train a long short-term memory (LSTM) recurrent neural network model. The accuracy of the predicted data is measured against real operational data to verify the reliability of the synthetic data creation methods and if these methods should be used in the future to inform inspectors of a reactor’s proliferation capabilities.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

A Second Moment Method for k -Eigenvalue Acceleration with Continuous Diffusion and Discontinuous Transport Discretizations

The second moment method is a linear acceleration technique that couples the transport equation to a diffusion equation with transport-dependent additive closures. The resulting low-order diffusion equation can be discretized independent of the transport discretization, unlike diffusion synthetic acceleration, and is symmetric positive definite, unlike quasidiffusion. While this method has been shown to be comparable to quasidiffusion in iterative performance for fixed source and time-dependent problems, it is largely unexplored as an eigenvalue problem acceleration scheme due to the belief that the resulting inhomogeneous source makes the problem ill posed. Recently, a preliminary feasibility study was performed on the second moment method for eigenvalue problems. The results suggested comparable performance to quasidiffusion and more robust performance than diffusion synthetic acceleration. This work extends the initial study to more realistic reactor problems using state-of-the-art discretization techniques. Finally, the results in this paper show that the second moment method is more computationally efficient than its alternatives on complex reactor problems with unstructured meshes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Cryogenic hydrogen fuel for controlled inertial confinement fusion (formation of reactor-scale cryogenic targets)

In inertial fusion energy research, considerable attention has recently been focused on low-cost fabrication of a large number of targets by developing a specialized layering module of repeatable operation. The targets must be free-standing, or unmounted. Therefore, the development of a target factory for inertial confinement fusion (ICF) is based on methods that can ensure a cost-effective target production with high repeatability. Minimization of the amount of tritium (i.e., minimization of time and space at all production stages) is a necessary condition as well. Additionally, the cryogenic hydrogen fuel inside the targets must have a structure (ultrafine layers—the grain size should be scaled back to the nanometer range) that supports the fuel layer survivability under target injection and transport through the reactor chamber. To meet the above requirements, significant progress has been made at the Lebedev Physical Institute (LPI) in the technology developed on the basis of rapid fuel layering inside moving free-standing targets (FST), also referred to as the FST layering method. Owing to the research carried out at LPI, unique experience has been gained in the development of the FST-layering module for target fabrication with an ultrafine fuel layer, including a reactor- scale target design. This experience can be used for the development of the next-generation FST-layering module for construction of a prototype of a target factory for power laser facilities and inertial fusion power plants.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Design Optimization of a Criticality Experiment for the Molten Chloride Reactor Experiment Facility

Neutronics simulations of Molten Chloride Fast Reactors have quantifiable biases that arise from nuclear data, modeling choices, or numerical methods. The multiphysics nature of molten salt reactors makes it challenging to disentangle neutronics modeling biases from biases originating from other physical phenomena. In comparison to a mock-up reactor, criticality experiments can specifically assess the neutronics modeling bias while limiting multiphysics effects. The criticality experiment must be neutronically representative of the full-scale reactor to be valuable. Here, in this paper, we describe the design of a criticality experiment to validate only the neutronics of TerraPower’s Molten Chloride Reactor Experiment (MCRE) and its criticality safety upset scenarios. The proposed experiment uses different chlorine-containing materials to maximize its similarity to the MCRE. The design process uses a constrained Bayesian optimization algorithm to investigate different objective functions that use covariance information for 35 Cl nuclear data. The experiments could reduce the nuclear data–induced uncertainty in k eff of the MCRE from 2161 to 886 pcm. They would also increase the upper subcritical limit of the MCRE criticality safety upset scenario from 0.94101 to 0.94476 when using the WHISPER analysis framework.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Multigroup Cross-section Generation in MCNP6.3 [Slides]

This presentation states that in comparison to the NJOY-produced multigroup cross sections, the MCNP-produced multigroup cross sections are generally consistent. Statistical uncertainties, however, are challenging and the unresolved resonance region may be looked at in the future. It also discusses how the SPM and LCS options were compared to each other for internal consistency. Additionally, some reactor pin-cell-like problems were used to compare to multigroup capabilities in other Monte Carlo codes (e.g., Serpent, OpenMC).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

SCALE 6.3.2 User Manual

SCALE is a comprehensive modeling and simulation suite for nuclear safety analysis and design developed and maintained by Oak Ridge National Laboratory under contract with the U.S. Nuclear Regulatory Commission, U.S. Department of Energy, and the National Nuclear Security Administration to perform reactor physics, criticality safety, radiation shielding, and spent fuel characterization for nuclear facilities and transportation/storage package designs. Visit the SCALE website for additional information.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Performing Numerical Analysis of Cybersecurity Options Using Dynamic Risk Analysis Tool EMRALD

Cyberattacks can have many different attack paths, durations, and goals. There are also many different mitigation options involving hardware, software, and/or humans. Considering a cyber threat should involve defense-in-depth methods and a quantitative or numerical evaluation of overall effectiveness against dynamic, time-dependent attacks to make cost and risk-informed decisions. Typical cyberattack modeling methods only provide a qualitative evaluation. The main areas of cybersecurity are confidentiality, integrity, and availability. For companies with cyber-physical systems such as advanced nuclear reactors, cyber-related safety is a requirement set by North American Electric Reliability and the U.S. Nuclear Regulatory Commission. They are also concerned about availability or reliability as a business case. As cyber threats are evolving to a business-for-hire structure, more attacks may focus on disrupting business success and reliability, causing financial and economic stability risk. Companies want to know business reliability and recovery from those threats, and that requires modeling physical behavior of the targets. Dynamic-state-based and Markov-based modeling provides a method for better cyber scenario modeling with different tools having issues such as state-base explosion. Dynamic modeling enables time and conditional features not found in other numerical evaluation methods. EMRALD (Event Modeling Risk Assessment using Lined Diagrams) is a dynamic risk analysis modeling and simulation tool and has features that reduce modeling issues. It has been used to model different time-dependent events including plant behavior and operator procedures. As a general modeling tool, EMRALD can also be used to model cyberattack scenarios with varying mitigation options and quantify effectiveness, producing numerical data for risk-informed decisions. This paper uses EMRALD to demonstrate that dynamic numerical risk analysis can be used for cyber threat modeling to provide insights for design decision-making and optimize defense strategies. Keywords: cyber modeling; cyber-physical systems; numerical cyber modeling

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Verification of High-Fidelity Neutronics Code PROTEUS for C5G7 Benchmark Problems

Under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the 3D whole core transport code PROTEUS has been developed at Argonne National Laboratory to perform the detailed 3D transport calculation for various reactor types including irregular or complex geometric reactors. PROTEUS includes two distinct transport solvers: SN and MOC. Among them, the MOC solver based on the extruded geometry, thus named MOCEX, can produce the faithful 3-D transport solutions by employing the rigorous and robust 2D/3D MOC formulation. This formulation overcomes the drawbacks of a full 3D MOC method requiring enormous memory and computational resources and resolves the instability issues of the 2D/1D method as well. The solution of MOCEX converges to that of a Monte Carlo code as the discretization of angle, space and energy is refined, unlike the 2D/1D method. As a verification and validation (V and V) effort of PROTEUS-MOCEX, this paper presents verification test results using the well-known 3D C5G7 benchmark problems by solving and comparing eigenvalues and pin powers with MCNP solutions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Scalable Gas Diffusion Electrode Fabrication for Electrochemical CO 2 Reduction Using Physical Vapor Deposition Methods

Electrochemical CO 2 reduction (ECR) promises the replacement of fossil fuels as the source of feedstock chemicals and seasonal storage of renewable energy. While much progress has been made in catalyst development and electrochemical reactor design, few studies have addressed the effect of catalyst integration on device performance. Using a microfluidic gas diffusion electrolyzer, we systematically studied the effect of thickness and the morphology of electron beam (EB) and magnetron-sputtered (MS) Cu catalyst coatings on ECR performance. We observed that EB-Cu outperforms MS-Cu in current density, selectivity, and energy efficiency, with 400 nm thick catalyst coatings performing the best. Additionally, the superior performance of EB-Cu catalysts is assigned to their faceted surface morphology and sharper Cu/gas diffusion layer interface, which increases their hydrophobicity. Tests in a large-scale zero-gap electrolyzer yielded similar product selectivity distributions with an ethylene Faradaic efficiency of 39% at 200 mA/cm 2 , demonstrating the scalability for industrial ECR applications.

36 MATERIALS SCIENCE↗