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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Initial study of notch sensitivity of Grade 91 using mechanisms motivated crystal plasticity finite element method

This report describes the development of a crystal plasticity finite element method (CPFEM) model for the multiaxial creep rupture of Grade 91 steel. The report demonstrates that core blocks in heat pipe reactor designs will see significantly different stress states than those typically found in conventional advanced reactor components where structures can often be represented as thin-walled pressure vessels. Current design practice uses effective stress measures calibrated to biaxial creep test data to account for the effect of stress multiaxiality on rupture life. Because these models rely on biaxial data they may not be accurate for the truly 3D, triaxial states found in heat pipe microreactor core blocks. A physically-based model, like the CPFEM model developed here, may be able to more accurately represent creep rupture under triaxial stresses when compared to conventional, empirical methods. This report describes the development, implementation, and validation of such a model and then applies the model to predict the effect of stress triaxiality on rupture in Grade 91 steel. The model results are used to assess the effective stress measures currently used in high temperature design codes and makes specific recommendations on improving these effective stresses to better account for the effect of 3D stress states on rupture in Grade 91. The validation and then implementation of these recommendations in design codes and standards could lead to safer, more effective heat pipe microreactor core block designs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Criticality and Radiation Shielding Impacts of ENDF/B-VIII.0 [Slides]

Several impacts of ENDF/B-VIII.0 on Radiation Shielding Applications were observed by researchers. Most reaction rates, attenuation ratios, and dose rates are similar for ENDF/B-VII.1 and ENDF/B-VIII.0. A significant reduction in reactor dosimetry results was noted. Investigation continues to reveal underlying data for discrepancies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Validation of Oregon State University High Temperature Test Facility Experiments Using Pronghorn

The OSU High Temperature Test Facility is a quarter-scale diameter, 1/64 scale volume test facility meant to replicate thermophysical phenomena in the prototypical General Atomics Modular High Temperature Gas Reactor. Tests pertaining to conduction cooldown events were performed from 2016-2019, providing a large database by which computational methods that are applicable to different length scales can be validated. One of these codes is Pronghorn, which is a coarse-meshed, porous-based subchannel thermal hydraulics code based on the MOOSE application, with the intention of better capturing the physics of both conduction and convection heat transfer within the OSU HTTF core. The goal of this summer project is to develop the framework by which Pronghorn can perform validation exercises of the HTTF core for benchmarking, by generating a mesh appropriate to the geometry of the HTTF core, developing input decks that accurately capture the initial and boundary conditions, materials, and relevant equations to the physics seen in the HTTF core, and using a postprocessor to compare simulation results to various experimental data. While validation of codes is a multi-year project, a mesh has been generated and tested in Pronghorn that meets mass conservation and basic heat transfer principles. The next step is to accurate depict the fluid inlet and outlet boundary conditions, which will be performed using computational fluid dynamics software.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Investigation of Coincidence Counting for Improving Minimal Detectable Activity of 110m Ag in Single Particle Gamma Analysis

Post-irradiation examination (PIE) of fuel particles from the fourth Advanced Gas Reactor Fuel Development and Qualification (AGR) Program irradiation (AGR-5/6/7) is being performed at Oak Ridge National Laboratory (ORNL). Tristructural isotropic (TRISO)-coated particles and associated compacts for the AGR-5/6/7 experiment fabricated by BWX Technologies Nuclear Operations Group were formed into a graphite matrix compact and irradiated at the Advanced Test Reactor at Idaho National Laboratory. At ORNL, particles are deconsolidated from the graphite matrix compact and individually scanned for emitted gamma rays with the Irradiated Microsphere Gamma Analyzer (IMGA). The IMGA system comprises a single high purity germanium (HPGe) detector, an automated particle handling vacuum system, and an ORTEC DSPEC-50 digital spectrometer for gamma ray analysis. IMGA quantifies gamma ray-emitting fission product inventories of individual TRISO particles, and these inventories can be compared with the measured average inventories per particle and radionuclide inventories predicted by AGR-5/6/7 physics calculations to determine if a particle experienced radionuclide release. Details on IMGA data collection methods can be found in the literature. The TRISO particle’s SiC layer provides structural support, as well as a barrier for fission product release during irradiation or subsequent safety testing. A weakened or compromised SiC layer can be identified by the release of radionuclides, such as 137 Cs, which is detected by IMGA. However, select radionuclides, such as 90 Sr, 110 mAg, and 154 Eu have been shown to migrate through an intact SiC layer. Measurement of the radionuclide 110 mAg is significant as its release has been shown to be particularly sensitive to in-reactor conditions (e.g., temperature) with broad variable particle to-particle release behaviors observed within a single compact. As such, 110 mAg activity is often used for particle selection for comprehensive PIE as bounding 110 mAg retention particles are hypothesized to represent limits in particle behaviors within a compact. As TRISO particle fuel PIE activities continue over time, IMGA measurements of the 110 mAg inventory are eventually hindered because of its relatively short half-life (~250 days). As the fuel ages from its end of irradiation (EOI) date, the measurement uncertainty and minimum detectable activity (MDA) of 110m Ag increase because the detector background continuum begins to dominate. For particles from the second AGR irradiation experiment (AGR-2), the 110m Ag MDA was above 20% of the calculated average particle inventory after approximately five half-lives, and 110m Ag activity was no longer measurable with IMGA after approximately seven half-lives. Therefore, coincidence counting approaches have been explored to determine feasibility of leveraging new approaches to overcome limitations associated with increasing MDA over time.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Discrete element simulation of Pebble Bed Reactors on graphics processing units

Prediction of pebble positions in a Pebble Bed Reactor (PBR) is necessary for both reactor physics and thermal hydraulics simulations as the arrangement of pebbles has a significant impact on the resulting core power, coolant flow, and fuel temperature. Knowledge of pebble movement as the fuel is cycled through the core is also critical for predicting the fuel residence time and subsequently, the fuel burnup. Simulation with the Discrete Element Method (DEM) can provide knowledge of both the fuel packing and the fuel movement during cycling. Previous works that have performed 3D full-core DEM simulation of PBRs have used simplified models that neglect reflector wall features. This work employs a graphics processing unit (GPU)-enabled DEM code, Project Chrono, to analyze the differences in pebble packing and pebble velocities between a simplified smooth PBR reflector and a more realistic reflector that includes circular wall features. Additionally, a sensitivity study is performed on the depth of the wall features to ensure that crystallization is prevented. Project Chrono is also validated for PBR cycling applications using experimental data. It is found that wall features with a depth of at least 0.5 pebble diameters significantly reduce crystallization in the near-wall region, leading to discrepancies in both packing fraction and pebble velocity in this region compared to the simplified reflector models. These discrepancies are found to lead to roughly a 5–10% difference in the prediction of the near-wall porosity and a 10% difference in the prediction of the velocity of pebbles near the wall. As a result of these discrepancies, it is suggested that future DEM simulations of PBRs include wall features to reduce modeling errors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Beyond the effectiveness factor: Multi-step reactions with intraparticle diffusion limitations

Though common in many industrial chemical processes, multi-step reactions occurring within porous catalyst particles provide significant challenges to previously established modeling frameworks. In this work, a system of reaction-diffusion equations is solved analytically via eigenanalysis. The intraparticle concentration profiles are integrated over the particle volume to obtain effectiveness factors for each reacting species, referred to herein as a multi-step effectiveness vector (MEV). The MEV solution is tested for the case of uniform flow past a stationary catalyst particle with a state-of-the-art, multi-step reaction scheme for fluid catalytic cracking (FCC). Outputs from particle-resolved direct numerical simulation (DNS) are compared to outputs from ordinary differentials with effective reaction closures coming from the MEV solution, the classic single-step effectiveness factor (SEF) (R. Aris, 'The mathematical theory of diffusion and reaction in permeable catalysts: The theory of the steady state,' Oxford University Press, 1975), and utilizing the intrinsic rates as-is. When compared to DNS, the MEV solution shows excellent agreement for the formation of all species and quantitatively captures the intraparticle concentration profiles. By contrast, the intrinsic treatment shows an over-prediction for product formation and the SEF treatment shows mixed levels of accuracy that depend upon the reaction coupling. Specifically, the SEF treatment accurately characterizes product formation if the kinetic scheme is weakly coupled but shows significant departure if stronger coupling is present. The computational cost associated with the MEV solution is associated with eigendecomposition of the Thiele matrix, an N x N matrix where N is the number of reacting components. Overall, the method is shown to be relatively easy to implement and may significantly enhance the fidelity of chemical reactor simulations that employ multi-step reaction schemes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Unified thermodynamic model to calculate COP of diverse sorption heat pump cycles: Adsorption, absorption, resorption, and multistep crystalline reactions

A straightforward thermodynamic model is developed in this work to analyze the efficiency limit of diverse sorption systems. A method is presented to quantify the dead thermal mass of heat exchangers. Solid and liquid sorbents based on chemisorption or physical adsorption are accommodated. Four possible single-effect configurations are considered: basic absorption or adsorption (separate desorber, absorber, condenser, and evaporator); separate condenser/evaporator (two identical sorbent-containing reactors with a condenser and a separate direct expansion evaporator); combined condenser/evaporator (one salt-containing reactor with a combined condenser/evaporator module); and resorption (two sorbent-containing reactors, each with a different sorbent). The analytical model was verified against an empirical heat and mass transfer model derived from component experimental results. It was then used to evaluate and determine the optimal design for an ammoniate salt-based solid/gas sorption heat pump for a space heating application. The effects on system performance were evaluated with respect to different working pairs, dead thermal mass factors, and system operating temperatures. The effect of reactor dead mass as well as heat recovery on system performance was also studied for each configuration. Based on the analysis in this work, an ammonia resorption cycle using LiCl/NaBr as the working pair was found to be the most suitable single-effect cycle for space heating applications. Furthermore the maximum cycle heating coefficient of performance for the design conditions was 1.50 with 50% heat recovery and 1.34 without heat recovery.

42 ENGINEERING↗
DOE OSTI2021

A Discussion of Strength Reduction Factor Development for Thermal Aging Effect on Nuclear Structural Alloys

In consideration of structural alloy property deterioration during long-term exposure to elevated temperatures, the yield and ultimate tensile strength reduction factors are provided in the Boiler and Pressure Vessel Code Section III for nuclear reactor component design and operation analysis. Because the Gen IV reactor requirement of 40 ~ 60 years of service life makes it difficult to acquire such long exposure test data for developing the reduction factors, they must be derived from test data with relatively short exposure by predictive methods considered to be reasonably reliable. A novel approach with a physically-based model has recently been proposed for application to development of reduction factors for 9Cr-1Mo-V. In the model, contributors to the tensile strength are first identified and related to definite microstructural features of the alloy, then some physically-based methods are employed to simulate the microstructural evolution, and finally the model is assembled with test-data-calibrated parameters to generate the yield and ultimate tensile strength reduction factors covering elevated temperature exposure for up to 57 years. The approach is undoubtedly a trailblazing development that will, if proven reliable, lead to a paradigm shift in predicting thermal aging behavior of many other alloys. Its debut application to Section III, however, concerns nuclear safety and naturally warrants objective, impartial, and thorough technical scrutiny. In the present paper, the novel and conventional approaches are discussed. Necessary improvements to the novel approach are recommended for its application to nuclear structural component design and analysis, and for its potential expanded use to other alloys.

Ren, Weiju↗
DOE OSTI2022

Assessment and Improvement of Fission Product Transport Predictions of Particle Fuel in BISON

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. This study was motivated by the need to incorporate more physics-based models in BISON in order to foster tri-structural isotropic (TRISO) applications. This document details the integration of new modeling capabilities in BISON, including (1) development of pyrolytic carbon (PyC) and silicon carbide (SiC) layer anisotropic thermal and mass transport capabilities, (2) verification of the mass diffusion solution in TRISO modeling, (3) calibration of fission product diffusivity using Advanced Gas Reactor (AGR) experiments, (4) improved fission product release modeling by developing compact diffusion modeling capabilities, and (5) documentation of accelerated failure analysis on the BISON website. Improvements made to the diffusion models and parameters were documented and validated against AGR-1 and -2 experiment data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1997

A DSMC Study of Low Pressure Argon Discharge

Work toward a self-consistent plasma simulation using the DSMC (Direct Simulation Monte Carlo) method for examination of the flowfields of low-pressure high density plasma reactors is presented. Presently, DSMC simulations for these applications involve either treating the electrons as a fluid or imposing experimentally determined values for the electron number density profile. In either approach, the electrons themselves are not physically simulated. Self-consistent plasma DSMC simulations have been conducted for aerospace applications but at a severe computational cost due in part to the scalar architectures on which the codes were employed. The present work attempts to conduct such simulations at a more reasonable cost using a plasma version of the object-oriented parallel Cornell DSMC code, MONACO, on an IBM SP-2. Due to availability of experimental data, the GEC reference cell is chosen to conduct preliminary investigations. An argon discharge is chosen to conduct preliminary investigations. An argon discharge is examined thus affording a simple chemistry set with eight gas-phase reactions and five species: Ar, Ar(+), Ar(*), Ar(sub 2), and e where Ar(*) is a metastable.

Hash, David B.↗
DOE OSTI2022

Multi physics simulation on Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS) by using Cardinal

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended durations or under accident condition. As one of the most significant components of a Reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Ageing degradation mechanisms (e.g. thermo-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Modeling In-Reactor Temperature Gradients in Lithium Aluminate Ceramic Pellets

In reactor thermal performance of lithium aluminate ceramic pellets are calculated to determine operating temperatures if internal components shift during irradiation. A 2D conduction model was setup to model this scenario with additional physics, including thermal expansion of components to reactor temperatures, creep down of the cladding due to differential pressure between the inside of the rod and reactor coolant system, and net radiation exchange between components. The non-linear coupling of the governing equations presented convergence challenges due to the very different effects from a small temperature change on the radiation vs radial and circumferential conduction equations. Initial temperature guesses for the multidimensional solution were provided from conservative 1D models simulating bounding scenarios. Convergence was greatly improved by taking small steps towards the root solutions and averaging with initial guesses. The final method robustly converged all 132 points at an axial plane after a few dozen iterations. The temperature results show that the additional physics have a secondary effect as the dominant heat transfer mechanism is radial conduction across the gas gap followed by circumferential conduction around components. Circumferential conduction in the lithium aluminate pellets is enhanced by the relatively high thermal conductivity of this ceramic.

Carstens, Nathan A.↗
DOE OSTI2021

A novel implicit hybrid machine learning model and its application for reinforcement learning

A novel methodology to develop implicit hybrid models is presented. PyTorch is used to integrate physics-based equations with machine learning models. Automatic differentiation of the hybrid model is leveraged to solve the implicit equations. Iterative solving enables gradient based updates to the machine learning model. The novel methodology is compared to an explicit hybrid approach on a continuously stirred tank reactor (CSTR). The novel method results in a lower modelling error. Both hybrid models effectively train with noisy data. To test the implicit hybrid model, it is employed as a reinforcement learning (RL) training model. The RL algorithm trained on the hybrid model outperforms real time optimization of the CSTR and performs nearly as well as RL trained directly on the CSTR and a traditional gradient based approach. Training RL directly on the CSTR requires over 60,000 system interactions compared to 6000 historical data points for hybrid model development.

42 ENGINEERING↗
DOE OSTI2025

Extending the Nuclide Inventory Validation Basis for High-Burnup Fuel with New Radiochemical Assay Data

Efforts are underway at Oak Ridge National Laboratory to improve the nuclide inventory validation basis for spent nuclear fuel at high burnups. Recently conducted radiochemical assay experiments provided new measurement data for nine samples of fuel irradiated in a pressurized water reactor, with estimated sample burnups in the 30 to 70 GWd/t range. This type of destructive assay data is essential for validating computational methods, tools, and nuclear data applied in nuclear safety analyses and for improving our understanding of the bias and uncertainty in code predictions. The measurement data include key actinides and fission products that span a gamut of needs and interests for nuclear science and engineering applications in criticality safety, reactor physics, nuclide inventory, decay heat, and radiation shielding. The SCALE 6.3 code system with ENDF/B-VII.1 cross-section libraries was used to simulate the irradiation histories of the measured fuel samples. The calculated nuclide concentrations are compared to corresponding measurement data. The significance of the comparisons is discussed, emphasizing how the addition of the new measurement data fills gaps in the validation basis at high burnups and contributes to the decrease in bias and uncertainty for predicted nuclide concentrations. The discussion addresses the effect of the sample burnup used in the simulation—which is based on reactor operator records or on calibration to measured data for burnup indicator fission products—on the validation results.

Nuclide inventory↗
DOE OSTI2023

Advanced Reactor Cyber Analysis and Development Environment (ARCADE) for System-Level Design Analysis

Cybersecurity is a persistent concern to the safety and security of Nuclear Power Plants (NPPs), but has lacked data-driven, evidence-based research. Rigorous cybersecurity analysis is critical for the licensing of advanced reactors using a performance-based approach. One tool that enables cybersecurity analysis is modeling and simulation. The nuclear industry makes extensive use of modeling and simulation throughout the decision process but lacks a method to incorporate cybersecurity analysis with existing models. To meet this need, the Advanced Reactor Cyber Analysis and Development Environment (ARCADE) was developed. ARCADE is a suite of publicly available tools that can be used to develop emulations of industrial control system devices and networks and integrate those emulations with physics simulators. This integration of cyber emulations and physics models enables rigorous cyber-physical analysis of cyber-attacks on NPP systems. This report provides an overview of key considerations for using ARCADE with existing physics models and demonstrates ARCADE’s capabilities for cybersecurity analysis. Using a model of the Small Modular Advanced High Temperature Reactor (SmAHTR), ARCADE was able to determine the sensitivity of the primary heat exchangers (PHX) to coordinated cyber-attacks. The analysis determined that while the PHX’s failures cause disruption to the reactor, they did not cause any safety limits to be exceeded because of the plant design, including passive safety features. Further development of ARCADE will enable rigorous, repeatable, and automated cyber-physical analysis of advanced reactor control systems. These efforts will also help reduce regulatory uncertainty by presenting similar types of cybersecurity analyses in a common format, driving standard approaches and reporting.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

High-burnup isotopics bias and uncertainty estimations for LANCR02

Calculated nuclide concentrations are compared against destructive radiochemical assays of modern, high-burnup BWR fuel to validate the LANCR02 lattice physics code. Good agreement is observed between different lattice physics methodologies - including the LANCR02 implementation of Method of Characteristics neutron and gamma transport solver and the TGBLA06 methodology - and against measurements. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗