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 397 records · Page 22

Comparison of Temperature-Dependent Cross Section Treatments Within the Shift Monte Carlo Radiation Transport Code

Multiple calculation methods have been developed to account for Doppler broadening of cross sections in Monte Carlo radiation transport simulations. This work compares the accuracy and performance of four different broadening methods, including an on-the-fly multipole treatment, applied to three different geometries of increasing complexity, including a nuclear reactor assembly with a realistic temperature profile. The comparisons reveal how the combination of physics implementations in the Shift Monte Carlo code result in slightly different temperature-dependent behavior, and different levels of performance, for a variety of materials and situations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NUMERICAL VERIFICATION OF THE RELAP-7 CORE CHANNEL SINGLE-PHASE MODEL

The RELAP-7 code is the next generation of nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). All the physics in RELAP-7 are fully coupled and the errors resulted from the traditional operator-splitting approach are eliminated. By using 2nd order methods in both time and space and eliminating operator-splitting errors, the numerical error of RELAP-7 can be minimized. Numerical verification is the process to verify the orders of numerical methods. It is an important part of modern verification and validation process. The core channel component in RELAP-7 is designed to simulate coolant flow as well as the conjugated heat transfer between coolant flow and the fuel rod. A special treatment at fuel centerline to avoid numerical singularity for the cylindrical heat conduction in the continuous finite element mesh is discussed. One steady state test case and one fast power up transient test case are utilized for the verification of the core channel model with single-phase flow. Analytical solution for the fuel pin temperature and figures of merit such as peak clad temperature and peak fuel temperature are used to define numerical errors. These cases prove that the mass and energy are well conserved and 2nd order convergence rates for both time and space are achieved in the core channel model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Bridging particle and reactor scales in the simulation of biomass fast pyrolysis by coupling particle resolved simulation and coarse grained CFD-DEM

Efficient utilization of biomass is challenging due to the complexity of the feedstock at molecular, particle, and reactor scales. Direct coupling of particle scale model and reactor scale simulation is computationally infeasible. In this work, the results from a particle resolved simulation were used to calibrate the heat transfer and reaction kinetics of particles. A hybrid drag model was used to consider the different fluidization properties of sands and biomass particles. The computation cost was reduced using a coarse-grained Discrete Element Method (DEM). The simulated conversion rate, char yield, tar yield, and pyrolysis gas yield compare well with experiment. The residence time and axial distributions of biomass were also analyzed. This novel multi-scale method provides an efficient and accurate tool for the modeling of biomass pyrolysis reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-Level Considerations for Access and Access Controls by Design

Nuclear power plants must be, by design and construction, robust structures and difficult to penetrate. Limiting access with difficult-to-penetrate physical barriers is going to be key for staffing reduction. Ideally, for security, the reactors would be sited underground, beneath a massive solid block, too thick to be penetrated by tools or explosives with all communications and power transfer lines also underground and fortified. Having the minimal possible number of access points and methods to completely block access from these points if a threat is detected will greatly help us justify staffing reduction.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Effects of Different Lattice Codes on Results of Sensitivity and Uncertainty Analysis

According to the UAM ('Uncertainty Analysis in Modeling') benchmark, uncertainties in the LWR system calculations at all stages should be determined and quantified. As the neutronics calculations are prerequisite for the predictions of the reactor system, and the uncertainties introduced by the neutronics calculations would be propagated to the subsequent calculations, hence uncertainty analysis has been firstly focused on the neutronics phase to quantify the uncertainties of the reactor-physics predictions in UAM. For the lattice physics, the objective of uncertainty analysis is focused on uncertainty quantifications of the few-group constants; for the core physics, uncertainties of the core steady-state stand-alone neutronics calculations are focused on. Based on the previous researches, the nuclear-data uncertainty has been treated as one of the most significant source of uncertainty and received the most focus. In the previous works, the domestic code UNICORN has the capability of performing sensitivity analysis and uncertainty analysis using the direct numerical perturbation (DNP) method and statistical sampling (SS) method. In this paper, the capability of the UNICORN code has been improved. The domestic lattice code NECP-CACTI has been added in UNICORN to perform the lattice calculations, with DRAGON5.0 and NECP-CACTI available for lattice calculations. For DRAGON5.0 and NECP-CACTI, all the methods and technologies are shared in UNICORN. Therefore, the effects of different lattice codes on the results of sensitivity and uncertainty analysis have been researched and determined in this paper. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Effect of parallel flow on resonant layer responses in high beta plasmas

Abstract Resonant layers in a tokamak respond to non-axisymmetric magnetic perturbations by amplifying the mode amplitude and balancing the plasma rotation through magnetic reconnection and force balance, respectively. This resonant response can be characterized by local layer parameters and especially by a single quantity in the linear regime, the so-called inner-layer Δ. The computation of Δ under two-fluid drift-MHD formalism has been progressed by reducing the order of the system in the phase space, where the shielding current is approximated as being only carried by electrons, a posteriori . In this study, we relax the approximation and compute Δ accounted for by the parallel flow associated with the ion shielding current. The posteriori is numerically verified in great agreement with the original SLAYER developed in a previous paper (J.-K. Park 2022 Phys. Plasmas 29 072506). Extending the resonant layer response theory to high β plasmas, our research findings answer two important questions: how the parallel flow influences the resonant layer response and why the parallel flow effect appears in high β plasmas. The complicated plasma compression in high β regime allows the parallel flow response to give rise to the ion shielding current, which not only shifts the zero-crossing condition of the ExB flow but also enhances the field penetration threshold. Technically, the Riccati matrix transformation method is adapted to handle the numerical stiffness due to the increased order of the system. The high fidelity of this numerical method makes use of further extension of the model to higher-order systems to take other physical phenomena into account. This work is envisaged to predict the resonant layer response under high β fusion reactor conditions.

Lee, Yeongsun (ORCID:000000034474416X)↗

UPF and Y-12 Criticality Accident Alarm System (CAAS) Radiation Testing Report

This report documents the experimental conditions and results for the testing of the Mirion Technologies CAAS-3S Criticality Accident Alarm System (CAAS) for the Y-12 National Security Complex (Y-12) and Uranium Processing Facility (UPF). The testing was conducted at the Godiva IV Burst Reactor at the National Criticality Experiments Research Center (NCERC) at the Nevada National Security Site (NNSS) during the week of January 11, 2021. The testing was conducted in accordance with 25774-QL-POA-JR00-00002-VDE-27.0, CAAS-3S Radiation Testing of Cabinet Components Test Plan. The tests subjected the CAAS-3S system to intense and short duration mixed neutron and gamma radiation fields and to high integrated radiation doses. The purpose of the test was to environmentally qualify various CAAS-3S components in a radiation environment. Successful completion of the tests was predicated on the collaborative efforts of several organizations. Sandia National Laboratory (SNL) and Lawrence Livermore National Laboratory (LLNL) provided dosimetry to estimate photon and neutron dose. Los Alamos National Laboratory (LANL) manages the NCERC facility, operated the Godiva IV reactor, and provided field support for CAAS-3S placement and communications interface for data acquisition systems supplied by Mirion Technologies. Throughout the test campaign, Y-12 and UPF personnel provided technical support in determining changes to the test schedule, as needed, interpretation of the results, and comparison between the measured dosimetry data and predicted results. Throughout the test campaign, Mirion personnel provided technical support in the operation of the CAAS-3S system itself.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Parallel Multi-Domain Solution Methodology Applied to Nonlinear Thermal Transport Problems in Nuclear Fuel Pins

This paper describes an efficient and nonlinearly consistent parallel solution methodology for solving coupled nonlinear thermal transport problems that occur in nuclear reactor applications over hundreds of individual 3D physical subdomains. Efficiency is obtained by leveraging knowledge of the physical domains, the physics on individual domains, and the couplings between them for preconditioning within a Jacobian Free Newton Krylov method. Details of the computational infrastructure that enabled this work, namely the open source Advanced Multi-Physics (AMP) package developed by the authors is described. Details of verification and validation experiments, and parallel performance analysis in weak and strong scaling studies demonstrating the achieved efficiency of the algorithm are presented. Furthermore, numerical experiments demonstrate that the preconditioner developed is independent of the number of fuel subdomains in a fuel rod, which is particularly important when simulating different types of fuel rods. Finally, we demonstrate the power of the coupling methodology by considering problems with couplings between surface and volume physics and coupling of nonlinear thermal transport in fuel rods to an external radiation transport code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Parametric study of natural circulation flow in molten salt fuel in molten salt reactor

The Molten Salt Reactor (MSR) is one of the most promising system proposed by Generation IV Forum (GIF) for future nuclear reactor systems. Advantages of the MSR are significantly larger compared to other reactor system, and is mainly achieved from its liquid nature of fuel and coolant. Further improvement to this system, which is a natural circulating molten fuel salt inside its tube in the reactor core is proposed, to achieve advantages of reducing and simplifying the MSR design proposed by GIF. Thermal hydraulic analysis on the proposed system was completed using a commercial computation fluid dynamics (CFD) software called FLUENT by ANSYS Inc. An understanding on theory behind this unique natural circulation flow inside the tube caused by fission heat generated in molten fuel salt and tube cooling was briefly introduced. Currently, no commercial CFD software could perfectly simulate natural circulation flow, hence, modeling this flow problem in FLUENT is introduced and analyzed to obtain best simulation results. Results obtained demonstrate the existence of periodical transient nature of flow problem, hence improvements in tube design is proposed based on the analysis on temperature and velocity profile. Results show that the proposed system could operate at up to 750MW core power, given that turbulence are enhanced throughout flow region, and precise molten fuel salt physical properties could be defined. At the request of the authors and the Proceedings Editor the name of the co-author Andrea Cioncolini was corrected from Andrea Coincolini. The same name correction was made in the Acknowledgement section on page 030004-10 and in reference number 4. The updated article was published on 11 May 2015.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A scalable DG solver for the electroneutral Nernst-Planck equations

The robust, scalable simulation of flowing electrochemical systems is increasingly important due to the synergy between intermittent renewable energy and electrochemical technologies such as energy storage and chemical manufacturing. The high Péclet regime of many such applications prevents the use of off-the-shelf discretization methods. In this work, we present a high-order Discontinuous Galerkin scheme for the electroneutral Nernst-Planck equations. The chosen charge conservation formulation allows for the specific treatment of the different physics: upwinding for advection and migration, and interior penalty for diffusion of ionic species as well the electric potential. Similarly, the formulation enables different treatments in the preconditioner: AMG for the potential blocks and ILU-based methods for the advection-dominated concentration blocks. Here we evaluate the convergence rate of the discretization scheme through numerical tests. Strong scaling results for two preconditioning approaches are shown for a large 3D flow-plate reactor example.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laying the Foundation for a Predictive Theory of Electroweak Nuclear Phenomena

Reliable model-independent predictions of nuclear processes such as beta-decay and neutrinonucleus scattering are important for many fundamental science and mission-related investigations; For example, accurate descriptions of as beta-decay are required for modeling the anti-neutrino spectrum of nuclear reactor, whether for trying to detect new physics such as sterile neutrinos, or remotely identifying the purpose of that reactor. For these activities, we need accurate calculations, as well as careful and thorough analyses of theoretical error and uncertainty in the underlying models and the many-body methods and codes used in solving that model. Unfortunately, modern theoretical calculations of these processes are inadequate for this task. At the heart of this inadequacy is the interaction between nucleons and other particles known as nuclear currents. Historically, descriptions of nuclear currents took into account only the simplest of particle interactions. Recently, the development of an effective field theory known as chiral EFT has been used to derive large corrections to these particle interactions. However, the complex nature of these corrections inhibited their use in ab initio methods of predicting nuclear processes. We developed a novel method to address these complexities and applied it in a series of benchmark and demonstrative calculations using state-of-the-art quantum few-nucleon solvers. Specifically, we studied the role of two-nucleon currents in electroweak processes using light nuclei (i.e., nuclei containing fewer than five nucleons) as a workbench to compute these processes and quantify their uncertainties using Bayesian extensions of the EFT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Study of iron structure stability in high temperature molten lead-bismuth eutectic with oxygen injection using molecular dynamics simulation

Corrosion of structural materials in high temperature molten lead-bismuth eutectic is a major problem for design of PbBi cooled reactor. One technique to inhibit corrosion process is to inject oxygen into coolant. In this paper we study and focus on a way of inhibiting the corrosion of iron using molecular dynamics method. For the simulation results we concluded that effective corrosion inhibition of iron may be achieved by injection 0.0532 wt% to 0.1156 wt% oxygen into liquid lead-bismuth. At this oxygen concentration the structure of iron material will be maintained at about 70% in bcc crystal structure during interaction with liquid metal.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Investigation of Xenon Oscillations via Point Kinetics

We all appreciate the special role Xenon-135 plays in reactor transients. Unless there is sufficient reactivity override, on shutdown a conventional reactor cannot immediately restart because of the buildup of strongly absorbing Xenon from Iodine decay and the absence of neutron flux to continue its elimination. We are also aware of how the flux imbalance in a large reactor core gives rise to Xenon induced power precession around the core in unpredictable ways. In addition, with regard to computational methods, we lack readily accessible numerical algorithms enabling reliable characterization of Xenon oscillations with simple models to explain this fundamentally important phenomenon. Nowhere is this more apparent than in the classroom. In this presentation, we first consider a new 12-element comprehensive decay model to Xenon and its corresponding numerical implementation and verification to establish a predictive model of Xenon buildup in a constant flux. We then follow with a reduced self-consistent flux model in an infinite reactor- a model not found in textbooks. This presentation should appeal to those interested in developing high order numerical methods for reactor transients.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Predicting departure from nucleate boiling with advanced data- and physics-driven approaches

The subcooled and low-quality saturated flow boiling critical heat flux (CHF) corresponding to the departure from nucleate boiling (DNB) is one of the major limiting factors in the design and operation of pressurized water reactors (PWRs). Due to lack of agreement on its physical mechanisms, development of an accurate DNB-type CHF model has been elusive. The various steady-state (SS) predictive tools that have been proposed can be grouped into two categories: (1) data-driven methods—namely best-fit empirical correlations and look-up tables that result in relatively good agreement with specific experimental datasets but often fail to extend beyond their ranges of validity, and (2) the physics-based mechanistic models that rely on reasonable yet limited understanding of the underlying physics and apply constitutive relations to close the conservation equations. On one hand, with recent advances in computational capabilities and optimization techniques, machine learning (ML)-based methods constitute an alternative and a more advanced, data-driven approach. With enough data, they can be particularly useful in engineering fields where the physical phenomena are complex and challenging to model. On the other hand, a state-of-the-art physics-driven approach would not only provide more accurate prediction of SS CHF, but it would also help better understand and model transient applications. The data-driven regression models use deep feed-forward neural networks (NNs) and random forests to cross-validate with 1,865 CHF test cases, covering a wide range of flow conditions and channel geometries. The best-estimate ML-based predictors compare favorably with the widely used look-up table for annulus and plate, and sensitivity analysis has confirmed their effectiveness. The key advantage of ML-based methods is their online extensibility of applicability domain. The proposed physics-driven approach combines key assumptions and parameters in the relatively well-accepted theories of liquid sublayer dryout and near-wall bubble crowding. A more realistic understanding of local mechanisms has been modeled. The new model has been optimized and validated against 1,439 tube data, and it shows considerably improved performance when compared to a recent mechanistic model, demonstrating unbiased close agreement with measurements over a wide range of operating conditions. Future work will extend the improved physics-driven model to non-tube geometry applications and will couple the model with ML via a hybrid approach. The mechanistic SS work will also serve to improve understanding and to model transient CHF scenarios.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A superstructure optimization approach for process synthesis under complex reaction networks

Here in this work, we present mixed integer linear programming methods for the synthesis of processes that involve complex reaction networks. Specifically, we consider the modeling of reactors and interconnecting streams in systems where the composition of the reactor inlet streams can vary substantially, thereby making the determination of the limiting component as well as the calculation of the stream heating/cooling and power requirements challenging. First, towards the modeling of reactors, we develop an extent-based method which detects the limiting reactant of each reaction occurring in parallel with others, based on the inlet flows of the reactants. Second, we develop a computationally tractable method for the calculation of the work and heating/cooling duty needed to condition any stream of a process based on simple calculations that can be performed offline. Finally, we present how the two aforementioned components can be integrated in an optimization model generated based on a process superstructure. We demonstrate the application of the developed methods for the synthesis of a biorefinery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Implementation of a Red-Black SOR CMFD Solver in MPACT

As part of the Consortium for Advanced Simulation of Light Water Reactors (CASL), the Virtual Environment for Reactor Applications (VERA) is being developed to provide high-fidelity multi-physics simulations of nuclear reactor cores. The MPACT code being developed collaboratively by Oak Ridge National Laboratory and the University of Michigan includes the primary deterministic neutron transport solver available within VERA. MPACT employs the 2D/1D method to solve three-dimensional problems, using 2D method of characteristics (MOC) radial transport solvers with 1D NEM-SP3 axial transport solvers, which are coupled using axial and radial transverse leakages. A 3D coarse mesh finite difference (CMFD) solver is used to provide accelerated convergence that is more stable. The Portable, Extensible Tool kit for Scientific Computation (PETSc) generalized minimal residual (GMRES) iterative solver with a block Jacobi preconditioner has been used to solve the CMFD linear systems. This work focused on the addition of a red-black successive over-relaxation (RBSOR) solver to reduce overall runtime by reducing time spent solving the CMFD linear systems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Machine Learning Assisted Safety Modeling and Analysis of Advanced Reactors

With the advances in computational power and numerical methods, analysts can now rely on first-principle simulations to predict ultra-fine details in a variety of applications. Advances in machine learning (ML) have produced algorithms that can now learn high-level abstractions via hierarchical models. This project aims to leverage advances in ML techniques and the available high-resolution simulation data to develop a novel modeling and simulation (M\&S) methodology for reactor safety analysis. While application-agnostic ML techniques are available, complex physics constraints need to be incorporated into ML techniques to build ML-based closures for computationally efficient predictive simulations. This project intends to develop a physics-guided data-driven multi-scale methodology for M\&S of advanced reactors. The project focuses on thermal fluid (T/F) phenomena, which play major roles in advanced reactor safety. Specifically, we propose a data-driven coarse-mesh turbulence model based on local flow features for the transient analysis of thermal mixing and stratification in a sodium-cooled fast reactor (SFR). The model has a coarse-mesh setup to ensure computational efficiency, while it is trained by fine-mesh computational fluid dynamics (CFD) data with Reynolds-averaged Navier-Stokes (RANS) turbulence model to ensure accuracy. Three different neural networks are developed and tested for loss-of-flow transients in the hot pool of SFR, i.e. the densely connected convolutional neural network (DCNN), long-short-term-memory network based on proper orthogonal decomposition (POD-LSTM), and the DCNN informed by LSTM (DCNN-LSTM). The performances of these three neural networks are evaluated based on baseline models. The DCNN-LSTM model has been chosen for further hyperparameter optimization. Furthermore, based on a simplified two-dimensional case, uncertainty quantification (UQ) of the developed ML-based closure are investigated with three methods, i.e. Monte Carlo dropout, deep ensemble, and Bayesian neural network. The developed ML-based turbulent viscosity closure relation based on deep ensemble is then integrated into the system analysis module SAM and serves as a term in the conservation equations. Such a SAM-ML based procedure guarantees that the obtained results are consistent with the physical constraints of the thermal-fluid system. The SAM-ML simulation on the same loss-of-flow transient showed comparable accuracy with the CFD simulation but with a much coarser mesh setup. Last but not least, the ML-based closure improvement with the support of higher-fidelity data from large eddy simulation (LES) is discussed. As a first step towards this direction, a baseline LES simulation is performed to obtain comparable data with RANS results. Based on the early results, future investigation on further improving the ML-based closure is discussed. We believe the developed approach that combines scientific machine learning with nuclear system analysis code can benefit the advanced reactor community as more accurate safety analyses will better characterize reactor safety margins and reduce licensing efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Use Machine Learning to Improve Burnup Measurement in Pebble Bed Reactors

Advanced pebble bed reactor (PBR) designs post new challenges in material control and accountancy (MC&A) because the fuel materials, distributed in many discrete pebbles, are continuously circulated through the reactor core and the refueling path compared to the bulk fuel assembly design in conventional reactors, e.g., light water reactors. In pebble bed reactors, there are hundreds of thousands of fuel pebbles in the reactor core during the normal operation, and the burnup of each pebble is measured when ejected from the core. Accurate burnup measurement is an important step in material control and spent fuel disposition. The measurement is usually based on detection of radiation signatures of fission products accumulated in the pebble fuel over irradiation in the core. Previous research has shown that height of photopeaks of fission products, such as 134 Cs, 137 Cs, 154 Eu, etc., can be used independently or in combination to infer or predict the level of burnup in the fuel. However, it remains challenging to measure such complex sources due to self-shielding effects, strong radiation background and intervening materials. Another operational challenge is the required high throughput of pebbles undergoing burnup measurement, which necessitates limited measurement time and thus impacts quality of measured gamma-ray spectra. Hence, advanced spectral analysis methods are needed to analyze the noisy gamma spectra and predict the burnup values.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗