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 289 records · Page 16

Ab initio molecular dynamics study of the properties of cerium in liquid sodium at 1000 K temperature

For liquid-sodium-cooled fast nuclear reactor systems, it is crucial to understand the behavior of lanthanides and other potential fission products in liquid sodium or other liquid metal solutions such as liquid cesium-sodium. In this study, we focus on lanthanide behavior in liquid sodium. Using ab initio molecular dynamics, we found that the solubility of cerium in liquid sodium at 1000 K was less than 0.78 at. %, and the diffusion coefficient of cerium in liquid sodium was calculated to be 5.57 × 10{sup −9} m{sup 2}/s. Furthermore, it was found that cerium in small amounts may significantly alter the heat capacity of the liquid sodium system. Our results are consistent with the experimental results for similar materials under similar conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coupling TRACE with nodal neutronics code Ants using the Exterior Communications Interface and VTT's multi-physics driver Cerberus

In this paper a recently developed coupling between system code TRACE and nodal neutronics code Ants is presented. The TRAC/RELAP Advanced Computational Engine (TRACE) is a reactor system code developed by U.S. Nuclear Regulatory Commission (NRC). The code has been designed to analyse loss-of-coolant accidents (LOCAs), operational transients and other accident scenarios in light water reactors. Ants is the reduced order nodal neutronics code of the Kraken framework. The diffusion solution method in Ants is based on the Analytic Function Expansion Nodal method (AFEN) and Flux Expansion Nodal Method (FENM). Rectangular, hexagonal and triangular geometries are supported. Both steady state and burnup problems including micro-depletion can be simulated. The coupling is implemented using VTT's multi-physics driver Cerberus and the Exterior Communications Interface (ECI) which is supplied with TRACE and enables the coupling of TRACE without source code modifications. The main focus in the paper is on describing the TRACE-Ants coupling in detail in order to share to other code developers what we have learned from the ECI based coupling process. The coupled code system is briefly tested by simulating a simple transient in a SMR core initiated by control rod movement. The results of the test calculation look reasonable and more complicated test cases will be modelled in the near future.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Energetic Light Fragment Production Capability in MCNP6

Monte Carlo Methods, Codes, and Applications group within the Computational Physics Division at Los Alamos National Laboratory has led the development of the transport code MCNP6 (Monte Carlo N-Particle transport code, version 6). MCNP6 is a general-purpose, continuous energy, generalized-geometry, time-dependent, Monte-Carlo radiation-transport code designed to track many particle types over broad ranges of energies. It is used around the world in applications ranging from radiation protection and dosimetry, nuclear-reactor design, nuclear criticality safety, detector design and analysis, decontamination and decommissioning, accelerator applications, medical physics, space research, and beyond. At lower energies, the code uses tables of evaluated nuclear data, while for higher energies (> 150 MeV), MCNP6 uses the cascade-exciton model, version 03.03 (CEM03.03), and the Los Alamos quark-gluon string model, version 03.03 (LAQGSM03.03) to model nuclear reactions. Emission of energetic heavy clusters heavier than {sup 4}He from nuclear reactions play a critical role in several applications, including electronics performance in space, human radiation dosages in space or other extreme radiation environments, proton- and hadron-therapy in medical physics, accelerator and shielding applications, and more. None of the available models are able to accurately predict emission of light fragments (LF) from arbitrary reactions. The CEM and LAQGSM event generators in MCNP6 describe quite well the spectra of fragments with sizes up to {sup 4}He across a broad range of target masses and incident energies (up to ∼ 5 GeV for CEM and up to ∼ 1 TeV/A for LAQGSM). However, they do not predict the high-energy tails of LF spectra heavier than {sup 4}He well. Most LF with energies above several tens of MeV are emitted during the precompound stage of a reaction. The 03.03 versions of CEM and LAQGSM do not account for precompound emission of LF larger than {sup 4}He. The goal of this research is to enable MCNP6 to produce high-energy light fragments. These energetic LF may be emitted by our models through three processes: Fermi breakup, preequilibrium, and coalescence. We explore the emission of light fragments through each of these mechanisms and demonstrate an improved agreement with experimental data achieved by extending precompound models to include emission of fragments heavier than {sup 4}He. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Computation of Air Discharge into a Water Pool by Using Volume-of-Fluid Method

Suppression tanks, in-containment refueling water storage tanks (IRWSTs) and filtered containment venting systems (FCVSs) of nuclear power plants use water pools to accommodate high-pressure and high-temperature discharge of non-condensable gases and/or steam. Dynamic behavior of gas bubbles inside a water pool determines pressure and temperature loading to the system and thus we need validated method for applications. The underlying physics relevant to these problems is still a difficult issue due to the complicated nature of gas- liquid interfacial transfers of mass, momentum, and energy in a two-phase flow. Therefore, notwithstanding many studies performed over the past two or three decades, two-fluid theory is still a formidable task and direct numerical simulation (DNS) needs huge computing cost. Pellegrini et al. analyzed steam condensation and chugging in the water pool of the PPOOLEX STB-24 test by using STAR-CCM+ code implemented with two-phase interface model based on the Rayleigh-Taylor instability theory (RTI) which triggers steam depressurization. For the solution of the fine interface, they used Eulerian-Eulerian model for the fluid domain since the volume-of-fluid (VOF) model requires larger mesh points. The RTI resulted in complete steam condensation consistent with the experiment. The fluid tank was open to the atmosphere. They reported the volume fraction distributions in the pool, not the pool temperature or pressure history. For non-condensing gas discharges, however, the VOF is still a good alternative to the Eulerian two-fluid or the DNS approaches. The spatial distribution of each phase is defined in terms of volume fractions and the cell physical properties are obtained from linear combinations of those of constituent phases weighted by each volume fraction. Therefore, conditions where thermal and phase change effects become important are less likely to be simulated realistically. The VOF model is suited for the simulation of flows where each phase constitutes a large structure with a relatively small total contact area between phases. Present work is also interested in this VOF method for application to air discharge into a water pool which occurs in the early operating stages of the said facilities. Our purpose is to validate the VOF method against an existing air discharge experiment of PPOOLEX by using a commercial CFD code, STAR-CCM+ version 10.04, in terms of macroscopic flow properties. Effect of turbulence is also discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Residual Stress Measurements of Cr-coated Cladding

Accident-tolerant fuel concepts for light-water reactor applications have been developed and tested in diverse research programs around the world. Industry teams have developed coated cladding concepts using different coating application methods, such as physical vapor deposition and coldspraying. Recently, variations in thickness and application method were investigated, and their effects on the microstructure and mechanical properties were reported. The results suggest an increase in performance in an unirradiated state, however, it remains unclear what phenomenon drove this performance. One explanation to this performance is associated with the coating process imparting a residual stress. Residual stress analysis was performed on coatings produced with various process parameters to understand how a thin layer of Cr impacts the mechanical properties after coating application and after simulated operating conditions. The gathered information is intended to address safety concerns or support margin identification related to accident-tolerant fuels and to support subsequent modeling and simulation efforts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

RELAP-7 Code Assessment Plan and Requirement Traceability Matrix

The RELAP-7, a safety analysis code for nuclear reactor system, is under development at Idaho National Laboratory (INL). Overall, the code development is directed towards leveraging the advancements in computer science technology, numerical solution methods and physical models over the last decades. Recently, INL has also been putting an effort to establish the code assessment plan, which aims to ensure an improved final product quality through the RELAP-7 development process. The ultimate goal of this plan is to propose a suitable way to systematically assess the wide range of software requirements for RELAP-7, including the software design, user interface, and technical requirements, etc. To this end, we first survey the literature (i.e., international/domestic reports, research articles) addressing the desirable features generally required for advanced nuclear system safety analysis codes. In addition, the V&V (verification and validation) efforts as well as the legacy issues of several recently-developed codes (e.g., RELAP5-3D, TRACE V5.0) are investigated. Lastly, this paper outlines the Requirement Traceability Matrix (RTM) for RELAP-7 which can be used to systematically evaluate and identify the code development process and its present capability.

42 ENGINEERING↗

Thermal scattering law for structure-dependent-Doppler broadening in FLASSH

Traditional Doppler broadening methods assume that target nuclei velocities follow a Maxwell-Boltzmann distribution which inherently assumes that the atoms are a free gas or that their velocities are independent of each other. This approximation is reasonable at high target temperatures and high neutron energies. However, nuclei are not independent: the lattice temperature and chemical binding will define the atomic motion. The thermal scattering law (i.e., TSL or S(α,β)) is a material property which describes the chemical binding and temperature response in terms of available momentum and energy states. In the thermal energy range, this TSL will define the thermal scattering cross sections. In the resonance region, the TSL offers a first-principles description of the probability distribution function for the velocity of the target nuclei. Using the TSL in Doppler broadening captures the structure of the material consistently from the thermal range into the resonance cross sections. In this work, both traditional free-gas and precise TSL Doppler broadening methods for resonance cross section evaluation have been implemented within the Full Law Analysis Scattering System Hub (FLASSH). This framework provides a generalized formulation for evaluating low-lying resonance data and streamlines nuclear data processing methods. TSL data generated using ab initio lattice dynamics (AILD) can be used to broaden cross section libraries in ENDF format for high-fidelity input into reactor physics calculations. These capabilities are demonstrated for the lowest absorption resonances of {sup 238}U in UO{sub 2}, UC and UN. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Microstructural and nanomechanical studies of PVD Cr coatings on SiC for LWR fuel cladding applications

The microstructure and nanomechanical properties of pure Cr coatings deposited on SiC with various physical vapor deposition (PVD) methods have been investigated for the mitigation of hydrothermal corrosion of SiC-SiC f (SiC fiber-reinforced SiC matrix) composite fuel cladding in light water reactors (LWR). Cr coatings (4–7 μm thick) were deposited on unbiased SiC substrates using six variants of magnetron sputtering processes: (i) standard (S-DCMS), (ii) pulsed (P-DCMS), (iii) ion-assisted (I-DCMS), and (iv) pulsed ion-assisted direct current magnetron sputtering (PI-DCMS), (v) high-power impulse magnetron sputtering (HiPIMS), and (vi) bipolar HiPIMS (B-HiPIMS). Microstructural characterization and nanoindentation testing were used to evaluate effects of deposition technique on coating microstructure and nanomechanical properties. The application of positive reverse pulses to the sputter target induced coarsening of the columnar grain structure, while ion bombardment of substrates during deposition promoted densification of the coating microstructure. Both effects are exemplified in B-HiPIMS deposition, resulting in a high-density microstructure with compressive residual stress. Material pile-up around nano-indents and plastic work during nanoindentation were analyzed to understand deformation behavior of the coatings. Finally, the results suggest the B-HiPIMS process to be among the promising methods for the deposition of Cr coatings on SiC-SiC f for LWR cladding application.

36 MATERIALS SCIENCE↗

Follow the power—pathways to steady-state tokamak reactors

Abstract Diagramming pathways of dimensionless power is a potent method for extrapolating between operating points on present-day tokamaks and future burning plasma devices. The heat transport power, current drive power, H-mode threshold power and other ‘plasma physics’ powers can be expressed in dimensionally correct (or normalized) form as Pa 3/4 , where P is the power and a is the plasma minor radius, with the relative gyroradius ( ρ * ) dependence ranging from gyro-Bohm-like for transport ( P a 3 / 4 ∝ ρ * − 3 / 2 ), Bohm-like for current drive ( P a 3 / 4 ∝ ρ * − 5 / 2 ) and worse than Bohm-like for H-mode threshold ( P a 3 / 4 ∝ ρ * − 3 ). The D–T fusion power cannot be normalized in the same fashion since it is governed by nuclear physics, but at fixed B T it scales like P a 3 / 4 ∝ ρ * − 9 / 2 . Other ‘mixed physics’ powers can be incorporated into the dimensionless power framework by holding B T fixed in the same manner. Diagramming these dimensionless powers vs ρ * shows how the pathway to a steady-state reactor can be optimized relative to various operational boundaries. Using a steady-state hybrid discharge with β N = 3.2 from DIII-D as the starting point, a multi-parameter optimization finds an attractive pathway to steady-state operation on ITER using 76 MW of current drive power (fusion gain of Q fus = 8), along with a pathway to Q fus = 20 in a JET-sized steady-state reactor with B T = 10 T.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Uncertainty Quantification of Activation Due to Cross Section Data in Neutron Shielding Calculation

Quantification of activation in structure materials around the nuclear reactor through calculation is necessary before start of decommissioning. Then a neutron shielding calculation is carried out to estimate activation of structural materials. Neutron shielding calculation is performed out for spatially large regions around a reactor through deep penetration calculation, calculation result would have significant uncertainty. In order to minimize total amount of structure materials treated as radioactive waste, quantitative estimation on uncertainty of structural activation and its reduction is important. The uncertainties due to calculation arises from that of input data and that of calculation method. Examples of uncertainty of input data is those for nuclear reaction cross section and atom density of the nuclide in the structure. In this study, we focus on the uncertainty of cross section, which would be one of the dominant uncertainties, and try to estimate uncertainty of activation due to cross section data in a neutron shielding calculation. The random sampling method (RS) is used as the uncertainty quantification method. With the RS method, uncertainty estimation can be carried out without modification of computer codes for a neutron shielding calculation. (authors)

36 MATERIALS SCIENCE↗

Evaluation of statistical variation of microstructural properties and temperature effects on creep fracture of Grade 91

This report describes three improvements made to a physically-based model for creep and creep rupture in Grade 91 steel developed as part of the Advanced Reactor Technologies program: 1. Progress on transitioning the model framework to the MOOSE finite element package to improve parallel scalability, increase the physical size of the simulations, and incorporate multiphysics effects into the model. 2. The development of a physically-based method of scaling the model parameters to accurately capture creep in Grade 91 for temperatures in the range of 450° to 500° C. 3. Extending the model to account for microstructural statistical variations, specifically capturing the effect of grain boundary energy on the physical parameters underling the model for grain boundary void nucleation and cavitation. Put together, these improvements result in an accurate, predictive model for the physical response of Grade 91 steel over the expected use temperature range for the material in future liquid metal cooled fast reactors. The model can be used to more accurately predict engineering properties of the material for very long service lives, which could lead to safer, more economical future advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Seventy-Five Years of Nuclear Criticality Safety Documents--A Bibliography

A collection of nuclear criticality documents was begun in the 1970s at the University of California’s Lawrence Livermore Laboratory (now Lawrence Livermore National Laboratory). Documents in this bibliography are dated as early as 1940. The early documents described experimental work that was of vital importance to the safe use of fissile materials at nuclear reactors, the facilities for preparing or transporting nuclear fuel, or its eventual storage. The experiments were of direct value to support safe operations with fissile materials, but they were also needed for validating early physics data and calculational methods and then for validating computational methods used in criticality safety evaluations. These documents provided information for our computerized bibliographic database. The bibliography was originally intended for use by Livermore staff members but criticality safety "pioneers” encouraged us to make the database available to the broader nuclear safety community and even shared copies of their bibliographies and reports. A publication of the bibliography is being prepared and a searchable database is available on the U.S. Department of Energy’s Nuclear Criticality Safety Program’s NCSP website. The NCSP website also includes links from the bibliography to many full-text documents. The database continues to expand--the information in this paper is based on a “snapshot” of its present status. The “vintage” hardcopy document collection is preserved and is being scanned into digital format.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Performing k eff Validation of As-Loaded Criticality Safety Calculations Using UNF-ST&DARDS: Applicable Experiment Selection

The general method for performing validation of as loaded criticality safety calculations using UNF ST&DARDS is presented in a paper by Clarity, which includes a description of the UNF-ST&DARDS system. Proof-of-principle analyses were performed in the summer of 2019 for MPC-32 dual purpose canisters (DPCs) containing pressurized water reactor (PWR) fuel assemblies. Summaries of these results are presented in this and a companion paper for this conference. The current paper describes the TSUNAMI-IP calculations performed to select applicable experiments for validation of 11 MPC-32 DPCs. The companion paper discusses the TSUNAMI-3D calculations used to generate sensitivity data to support the experiment selections discussed here. Experiment selection is based on the sensitivity/uncertainty (S/U) methods used to validate criticality safety calculations of as-loaded DPCs containing pressurized water reactor (PWR) spent nuclear fuel (SNF). This process has been demonstrated and is summarized in this paper. The approach is similar to that used in NUREG/CR-7109, which provides an approach for validation of PWR burnup credit (BUC), including major and minor actinides and major fission products. The premise of S/U-based validation is that applicable experiments—those having a similar bias to a given application system—will have similar sensitivities for each isotope and reaction in the two systems. It is assumed that cross sections with larger uncertainties are more likely to contain data errors which contribute to the bias. The integral index c k thus propagates the system sensitivities with the nuclear covariance data to calculate a correlation coefficient representing the similarity of the two systems. In this work, a c k value of 0.8 or higher is interpreted as identifying an experiment with sufficient similarity for use in validation. This paper presents a brief summary of the characteristics of the 11 MPC-32 DPCs used in the proof of-principle analysis for as-loaded criticality safety calculation validation and an overview of the critical experiment suite with which each of these DPC models was compared. A summary and discussion of c k results is also presented, followed by conclusions and a discussion of future work to be performed for validation of UNF ST&DARDS as-loaded criticality safety calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Improved Variable-Transformation Algorithm for the Hexagonal Variational Nodal Method

Hexagonal Variational Nodal Method (VNM) has been widely employed in fast reactor core calculation codes. It provides the weak form solution of the second order transport equation. It employs a direct three-dimensional expansion instead of the transverse integration technique, with the singular term caused by transverse integration over the hexagonal nodal avoided. However, as the increasing Pn expansion order in angle, the number of moments increases dramatically, which significantly increases the requirements in both memory and computational time. As a consequence, Yang proposed the Variable-Transformation Algorithm. It firstly used the reverse Cuthill-McKee (RCMK) reordering method to transform the response matrices into block-diagonal. Then, for the unreducible matrices, orthogonal transformation was applied, after which the matrices were sparse enough to be reduced by using the RCMK reordering once more. In this summary, based on a newly developed hexagonal VNM code NECP-VioletHex, the Variable Transformation Algorithm was investigated and improved by adopting more constraints to the basis functions to further take advantage of the symmetric properties of the hexagonal-z nodes. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Parallel Performance of the Time Dependent Transport Code TDKENO with Applications to TREAT Simulations

Renewed interest in high fidelity simulation of excursion events has prompted the improvement of codes that solve the time-dependent Boltzmann transport equation. Often these codes make approximations to the transport equation in order to achieve results in a reasonable amount of time. One method, the improved quasi-static (IQS), makes few approximations compared to adiabatic, diffusion, and point kinetics. This is the methodology employed by the code TDKENO. One downside to this approach is a portion of the calculation (flux shape) is done with a Monte Carlo code and must be called upon several times. To minimize run time the flux shape calculation is solved with a modified version of KENO-VI from Scale 6.2. This version of KENO runs in parallel across hundreds of nodes, thus improving performance. The parallel capabilities of KENO allows TDKENO to simulate complex transient experiments with a large number of histories. These improvements are applied to the simulation of TREAT calibration experiments. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Survey of Modeling and Simulation Techniques for Advanced Manufacturing Technologies Volume I – Predicting Initial Microstructures

This report summarizes the current state of modeling and simulation methods for predicting the initial structure and properties of material in components assembled using advanced manufacturing technologies (AMTs). The report is the first volume in a two-volume series. This first volume focuses on predicting initial microstructures of AM material. The second volume discusses predicting material properties given the initial microstructure. The focus is on technologies of particular relevance to the design and manufacture of nuclear reactor components. The purpose of the report is to help develop the technical knowledge base to support regulatory decisions that will be needed to assess nuclear components manufactured with AMTs as these components are installed at nuclear power plants (NPPs). The report develops a list of AMTs of particular interest to the NRC and summarizes the key microstructural features and corresponding processing parameters relevant to each technology. Further sections describe available physically-based and data-driven prediction methods and survey widely available software tools. The report concludes with a summary of gaps that may be of particular interest to the NRC when evaluating modeling and simulation methods for AMTs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗