Underlying mathematical relations for efficient COMET response expansion coefficient and core calculations
Not Available
Engineering topics
Publications and source records attributed to Rahnema, Farzad.
Not Available
In the mid-1970s and early 1980s, almost one hundred critical and approach-to-critical experiments using water-moderated plutonium metal arrays were performed at the Rocky Flats Critical Mass Laboratory. The experiments were performed in two programs: 2x2xN arrays conducted in the 1970s and 3x3x3 arrays conducted in the 1980s. Both sets of experiments used the same double-canned plutonium parts and were water moderated. The major differences between the two programs were the support structure and tank and the administrative limits on criticality. This paper summarizes the critical data and results of both sets of experiments. In conclusion, these experiments are then analyzed in a companion paper.
As part of the nuclear data evaluation and validation cycle, the ENDF/B-VIII.0 cross-section library released in 2018 requires testing to determine areas of improvement and deterioration. Previous work by the authors investigated the performance of 16 O, 56 Fe, and 63,65 Cu cross sections, with this study acting as an extension of the prior work. In addition to the isotopes and nuclear criticality safety benchmarks of interest to the prior work, benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook were selected for their k eff sensitivity to 1 H, C, 58,60 Ni, 182,183,184,186 W, 235,238U, or 239 Pu cross sections and were modeled in the SCALE code system maintained by Oak Ridge National Laboratory. In total, 253 benchmark configurations were selected for their sensitivities and modeled using SCALE 6.2.4 Criticality Safety Analysis Sequences (CSAS) continuous-energy Monte Carlo k eff calculations. This collection includes and expands upon the 99 benchmarks in the prior work. The AMPX-processed ENDF/B-VIII.0 library was decomposed into individual ENDF/B-VIII.0 datum libraries for each isotope of interest. Doing so allowed for the individual substitution of an ENDF/B-VIII.0 cross section in the place of ENDF/B-VII.1, determining isotope-specific effects of ENDF/B-VIII.0 relative to ENDF/B-VII.1. Full library calculations with entirely ENDF/B-VII.1 data or entirely ENDF/B-VIII.0 data were also executed. As a measure of performance, the average relative deviation was determined as the ratio of the deviation between calculated and experimental keff to the propagated calculational and experimental uncertainty. With calculated full library and isotope-specific ENDF/B-VIII.0 k eff ’s, an optimized combination of data libraries was estimated and confirmed with SCALE calculations. This showed that reverting 239 Pu, 58 Ni, 16 O, and 65 Cu cross sections to ENDF/B-VII.1 resulted in improved performance relative to the full ENDF/B-VIII.0 library. Across all 253 benchmarks, the average relative deviation was 1.29σ for the full ENDF/B-VII.1 library, 1.17σ for the full ENDF/B-VIII.0 library, and 0.97σ for the optimized combination. The reversion of 239 Pu, 58 Ni, 16 O, and 65 Cu cross sections to ENDF/B-VII.1 in the 99 benchmarks of the prior work resulted in further improved experimental agreement compared to the previously reported improvement from 16 O and 65 Cu alone. Therefore, it is suggested that applications with significant sensitivities to 239 Pu, 58 Ni, 16 O, and 65 Cu consider their choice of nuclear data library.
Not provided.
Recently completed cross-section evaluations sponsored in part by the Nuclear Criticality Safety Program were incorporated into the 2018 release of the ENDF/B-VIII.0 cross-section library. Evaluated isotopes of interest to the nuclear data and criticality safety community include 16 O, 56 Fe, and 63,65 Cu. For performance validation, benchmark models defined in the International Criticality Safety Benchmark Evaluation Project Handbook were selected based on energy-integrated k eff sensitivities to total cross sections of interest and compared with experimental values. Of the 102 benchmark configurations that were utilized, 63 are sensitive to 16 O, 32 sensitive to 63,65 Cu, and 25 sensitive to 56 Fe. Selected benchmarks were modeled in SCALE 6.2.3 Criticality Safety Analysis Sequence (CSAS) continuous-energy Monte Carlo k eff calculations with ENDF/B-VII.1, with a hybrid ENDF/B-VII.1 with ENDF/B-VIII.0 data substituted for individual isotopes of interest, and with ENDF/B-VIII.0. ENDF/B-VIII.0 showed improved agreement with experimental k eff for 56 Fe, 63 Cu, elemental copper, and full library substitution while producing lessened agreement for 16 O and 65 Cu. Furthermore, with full library and isotope-specific ENDF/B-VIII.0 performance, a best-case ENDF library was formed by excluding underperforming isotopes’ ENDF/B-VIII.0 data, reverting 16 O and 65 Cu cross sections to ENDF/B-VII.1. This resulted in the average relative deviation between calculated and experimental data improving from 1.45σ for the ENDF/B-VIII.0 library to 1.32σ for the best-case library, relative to benchmark uncertainty.
This paper details and implements a framework for evaluating thermal neutron scattering cross sections that provide S(α,β) data and covariance data for hydrogen in light water. This methodology involves perturbing model parameters of molecular dynamics potentials and fitting the simulation results to experimental data. The framework is general and can be applied to any material or simulation method. The fit is made using the Unified Monte Carlo method to experimentally measure double-differential scattering cross sections of light water at the Spallation Neutron Source at Oak Ridge National Laboratory. Mean values and covariance data were generated for model parameters, phonon density of states, double-differential cross sections, and total scattering cross sections. These posterior parameter values were very similar to their prior values with a maximum relative error of 0.54%. This falls within in the Unified Monte Carlo–calculated uncertainties on the order of 2.7%. Additionally, posterior double-differential cross sections agree favorably with ENDF/B-VIII.0 cross sections. The new thermal scattering law was tested by comparing it against benchmarks from the International Criticality Safety Benchmark Evaluation Project Handbook, which showed a slight improvement over the ENDF/B-VIII.0 library. Additionally, the covariance matrix of the phonon density of states was validated to confirm that the spread of k eff from the density of states used to generate the covariance matrix was similar to the spread of k eff from the density of states of the sampled covariance matrix.
Not provided.
Some of the radiation transport solvers that are based on expansion function methods in hexagonal geometry generally use projection techniques that transform the spatial dependence of the angular flux in hexagonal geometry to Cartesian geometry by mapping and then expanding in terms of expansion functions (e.g., Legendre polynomials). However, this approach usually requires a relatively higherorder expansion in space for accurate representation of the spatial distribution of the flux. To overcome this limitation, a set of consistent spatial expansion functions have been developed in hexagonal geometry in this paper. These functions satisfy the orthogonality condition on hexagonal surfaces and guarantee the conservation of currents after expansion regardless of the truncation order. The latter is an important property for preservation of neutron balance. The new spatial expansion method has been implemented into the coarse mesh radiation transport (COMET) code and tested in both the single assembly and the full core Advanced High Temperature Reactor (AHTR) benchmark problems. The eigenvalues and assembly-averaged and stripe-wise fission density distributions predicted by the new method with 2nd order spatial expansion functions were found to exhibit similar accuracy as those computed by the current mapping method with 4th order spatial expansion functions. This indicates that the angular fluxes in AHTRs can be predominantly represented by a low order expansion if the new spatial expansion functions are used. As a result, to maintain the same accuracy, the computational efficiency of COMET is improved significantly (2–5 times) by using the new spatial expansion method.
This document constitutes progress made on this project by the Georgia Tech led Integrated Approach to Fluoride High Temperature Reactor Technology and Licensing Challenges (FHR-IRP). This is the first draft of the project final report. Chapter 3 is focused on Tritium management. Tritium management in Fluoride-salt High-temperature Reactors (FHRs) presents a major operational concern and a significant challenge. Tritium is generated via multiple pathways including various interactions in fuel and structural materials, as well as coolant (primarily interaction with lithium and beryllium), leading to a complex evolution and spatial distribution of tritium generation rate. FHR generates orders of magnitude more tritium than in LWRs, and has potential for a significant release, via permeation through structural materials (due to high temperature), and through heat exchanger(s). Therefore, accurate estimate is necessary. For that purpose relevant tritium production pathways were identified. A methodology to obtain space-time dependent tritium generation rate was developed and implemented, employing SCALE and MCNP codes. Detailed models of the reference Advanced High Temperature Reactor (AHTR) design were developed. A multi-parametric metrics was introduced to capture relevant tritium production indicators, including the instantaneous tritium production rate, total tritium production (neglecting subsequent decay), and accumulated tritium production (accounting for the decay). The obtained results will inform tritium management, including the sizing of the components, and to gain insight on the impact of the core design parameters on tritium generation rates. Research effort on Liquid Salt Impurity Removal, Redox and Corrosion Control is discussed in Chapter 4. Section 4.1 is an introduction chapter. Based on the current available knowledge and operation experience on previous MSRs and other types of nuclear reactors, the report identifies the major impurities that are related to either material compatibility, salt properties or radioactive source terms, and the potential contamination sources of these impurities, their transport in the primary coolant loop and available cleanup methodologies for a particular group of impurity. In section 4.2, the available purification methods were discussed, including the liquid/liquid extraction using liquid Bi-Li which has been applied in MSBR and electrochemical purification which is developed based on Pyroprocessing. In this section, the new multi-stage and multi-salt-component was developed and discussed. Section 4.3 discussed the redox potential control for salt chemistry and corrosion mitigation, including the control theory and the soluble/soluble buffer control method. The section also developed method on how to analyze soluble/soluble redox reaction in a molten salt system and designed an electrochemical sensor for redox potential control. Corrosion models were also reported in this section. Section 4.4 focused on how to measure the fundamental properties in a molten salt. As an example, the section focused on La fluoride in the molten salt. The section developed methods on how analyze electrochemical measurement data to obtain the fundamental properties including the redox potential, diffusion coefficient and exchange current density. All these fundamental information is important for salt purification system development. Section 4.5 focused on measurements of corrosion products in the salts. Because of the different behaviors between Cr and La in the salt, new analysis method for Cr fluorides in the salt was developed. Some results for Fe fluoride in the salt and corrosion measurements were also reported. Section 4.6 demonstrated the separation of oxide from the salt. This section also measured the solubility of La 2 O 3 in the salt as function of temperature. Section 4.7 focused on the separation of La and Ce from the molten salt. These two elements were selected because La was found to be separated by using inert electrode, while for separation of Ce, reactive electrode must be selected. Section 4.8 demonstrated the separation iodine ions from the molten salt. Two kinds of working electrode materials (metal W and graphite) were investigated. Section 4.9 focused on the fundamental data of fission products and actinides in liquid Bi. The chapter reviewed the available data, then based a CALPHAD model, the available data were broadened to the concentration-dependent cases.
Not Available
Not Available
On the path to deployment of any reactor, modeling tool verification and validation is a key step. Fluoride-Salt-Cooled High Temperature Reactors (FHR) pose challenges to neutronics modeling and simulation tools due to several design features of the reactor type. This paper presents a categorized list of phenomena that pose challenges to FHR modeling using current neutronics tools. These phenomena are presented in four categories: Fundamental Cross Section Data, Material Composition, Computational Methodology, and General Depletion. Here, a short path forward for these phenomena is also presented. In addition, a discussion of the resulting gaps in current codes is presented.
Here, this study identifies the phenomena that impact the accuracy of the multiphysics analysis of fluoride salt-cooled high temperature reactors (FHR) and categorizes the level of coupling needed to model these phenomena with sufficient accuracy in the context of three scenarios, namely, normal operation, station blackout, and simultaneous withdrawal of all control rods. Such a study is a necessary step prior to developing or extending existing multiphysics modeling and simulation tools. The identified multiphysics phenomena are categorized into the groups requiring ‘‘tight” and ‘‘loose” coupling. ‘‘Tight” coupling refers to phenomena that require iterative feedback between multiple physics modules/methods whereas ‘‘loose” coupling refers to phenomena that require only data from other physics that can be pre-calculated and shared in tabular or functional form. For the station blackout scenario, it was determined that no identified phenomena required tight coupling. For the rod withdrawal and normal operation scenarios, it is found that there is a large number of phenomena that require tight coupling between neutronics and thermal hydraulics.
With the rise in interest of GEN-IV reactor systems, newer, and more accurate thermal scattering nuclear data are needed. Both GEN-IV reactors and current light water reactors applying for license extensions require high fidelity cross sections and uncertainties to better thermal margins and thereby maximizing energy production and improving economics and safety. Thermal moderator data play a key role in both reactor physics and nuclear criticality safety analyses. Due to the lack of uncertainties data for thermal scattering material, there is no way to quantify the effects of thermal scattering uncertainties in quantities of interest (e.g., k eff ) in reactor or criticality safety systems.
Explore the source record for details and available documents.
Not Available
One new method for constructing probability tables in the unresolved resonance region (URR) has been developed. This new methodology is an extensive modification of the single-level Breit-Wigner (SLBW) pseudo-resonance pair sequence method commonly used to generate probability tables in the URR. The new method uses a Monte Carlo process to generate many pseudo-resonance sequences by first sampling the average resonance parameter data in the URR and then converting the sampled resonance parameters to the more robust R-matrix limited (RML) format. Furthermore, for each sampled set of pseudo-resonance sequences, the temperature-dependent cross sections are reconstructed on a small grid around the energy of reference using the Reich-Moore formalism and the Leal-Hwang Doppler broadening methodology. We then use the effective cross sections calculated at the energies of reference to construct probability tables in the URR. The RML cross-section reconstruction algorithm has been rigorously tested for a variety of isotopes, including 16 O, 19 F, 35 Cl, 56 Fe, 63 Cu, and 65 Cu. The new URR method also produced normalized cross-section factor probability tables for 238 U that were found to be in agreement with current standards. The modified 238 U probability tables were shown to produce results in excellent agreement with several standard benchmarks, including the IEU-MET-FAST-007 (BIG TEN), IEU-MET-FAST-003, and IEU-COMP-FAST-004 benchmarks.