Search NASA⌕ Search

Engineering topics

Brown, Jesse M.

Publications and source records attributed to Brown, Jesse M..

At least 19 records

Recognizing Unrecognized Sources of Uncertainty (USU) in Nuclear Data

Historically, time-of-flight (TOF) nuclear cross section measurements on different nuclides are assumed to be uncorrelated if they were recorded in different facilities with identical methods, identical facilities with different methods, and even identical facilities with identical methods. Ideally, measurements of different nuclides would truly be uncorrelated thus providing independent assessments of some cross section. In reality, correlations exist between measurements but are simply assumed to be unimportant. To eliminate these qualitative assumptions, in this paper we make a counter-intuitive suggestion to perform an intentionally correlated measurement of energy-differential fission and capture reactions for nuclides in a single criticality safety benchmark during a single experimental campaign. While this would introduce undesirable correlations, it would fully quantify correlations between datasets, rather than assume that the correlations do not exist.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A new 181 Ta neutron resolved resonance region evaluation

A new 181 Ta neutron resolved resonance region evaluation has been performed from the thermal energy range up to approximately 2.5 keV. The R-matrix SAMMY code was used with the Reich–Moore approximation to evaluate resonance parameters from several experimental data sets. A Monte Carlo approach was used for resonance spin assignments and generating 59 small fictitious resonance levels which were shown to improve the cumulative level, Porter-Thomas, and Wigner distributions as compared to theoretical predictions. Covariance information was also generated for the entire resolved resonance region. Finally, the positive impact of the new evaluation was validated through benchmark calculations which were sensitive to the 181 Ta cross section and showed improvement in the reactivity bias for several benchmark cases.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Generalized Bayesian Framework for Evaluation of Integral Benchmark Experiments

A recently published generalized Bayesian optimization framework has provided a way to retract any or all of the three common assumptions underlying the conventional Generalized Linear Least Squares (GLLS) optimization method based on the concepts introduced in Ref. [2]. These assumptions are: 1. Perfection: The model used for data evaluation and the prior probability distribution function (PDF) of generalized data are perfect. 2. Normality: The prior and posterior PDF are normal. 3. Linearity: The model is linear. In this work we outline how the framework in [1] could be directly adopted for improved evaluation of nuclear criticality integral benchmark experiments (IBEs) by: 1. Removing the first assumption alone by utilizing the concept of imperfections introduced in [1] to enable evaluation in the presence of discrepancies between the data and model or of missing covariance information by a GLLS method that will be seen as a generalization of the conventional GLLS method employed by the TSURFER code, and by 2. Removing the remaining two assumptions by implementing a Markov Chain Monte Carlo method for computation of the posterior PDF in the SAMPLER code, where TSURFER and SAMPLER are the uncertainty quantification (UQ) codes for IBEs in the SCALE code system based on the GLLS and the stochastic method, respectively. The graphic in Figure 1 categorizes the methods discussed in terms of the assumptions that they employ to determine posterior PDFs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

New Measurements to Resolve Discrepancies in Evaluated Model Parameters of 181 Ta

To resolve discrepancies in evaluated cross sections among major nuclear data libraries, energy-differential neutron transmission and radiative capture yield of 181 Ta were measured from 0.15 to 100 keV using multiple sample thicknesses. The new measurements provide resolution such that the resolved resonance region (RRR) can be evaluated up to at least 2.5 keV and the unresolved resonance region can be evaluated up to at least 100 keV. The transmission and capture yield measurements were modeled using resonance parameters from three major libraries to assess the predictive capability of each. It was found that JENDL-5.0 performed best in the RRR. Because of the poor performance of the U.S. ENDF/B evaluation, it is recommended that ENDF/B be reevaluated for 181 Ta.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optimal Chloride Salt Mixture for a Fusion Blanket

Deuterium-tritium fusion reactors cannot operate for a significant period without a closed tritium fuel cycle, according to a recent National Academies of Sciences (NAS) report on bringing fusion reactors to the US electrical grid. This fact places breeder blankets as one of the foundational systems for self-sustained fusion reactor operation. The typical functional requirements for a breeder blanket system include producing tritium, absorbing kinetic energy, transporting thermal energy, and being environmentally attractive. State-of-the-art research on liquid blankets has converged to primarily focus on (LiF) 2 and BeF 2 (FLiBe) molten salts and a metallic eutectic of lead and lithium (PbLi), but “virtually all of the technologies related to the tritium fuel cycle are at a low technological readiness level”. This work sought to explore optimum blanket configurations as it aligned with the Oak Ridge National Laboratory (ORNL) FY 2023 Laboratory Directed Research and Development Program’s research priority of developing and expanding the current understanding of fusion blanket science and technology. This purpose of this work was to address ORNL research priorities and NAS recommendations by investigating novel liquid blanket materials that could provide self-sustaining operation and draw on experience from research on molten salts used for advanced fission reactors, concentrated solar, and thermal energy storage. The hypothesis when proposing this research was that there could be chloride-based blanket designs that can exceed the tritium breeding ratios of (FLiBe) molten salt blankets while reducing the use of Be (FLiBe), avoiding the generation of HF (FLiBe), and minimizing magnetohydrodynamic (MHD)-perturbed flow fields (PbLi). The fastest and most cost-effective path to deploying liquid fusion breeder blankets could be from maximizing the synergistic technological overlap between fusion, fission, concentrated solar, and thermal energy storage industries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Zirconium Nuclear Data Campaign: Measurement of 90 Zr ( n, γ ) Cross Section

The isotopes of Zr with A = [90, 91, 92, 94] make up more than 97% of naturally occurring Zr and are important to many nuclear applications such as nuclear reactors. One of the attractive qualities of naturally occurring Zr isotopes is that they have a low σ γ /σ t ratio at most neutron energies. Thus, they improve the neutron economy in reactors by preferentially scattering neutrons rather than absorbing them. This same quality also presents a challenge to measuring the capture cross section, σ γ , of Zr isotopes. The ENDF/B VIII.0 library has a relative uncertainty of approximately 10–20% for incident neutron energies < 0.1 MeV and an uncertainty greater than 20% for energies > 0.1 MeV for the majority of natural Zr isotopes. This motivated the Nuclear Criticality Safety Program to embark on a campaign to accurately measure and evaluate these Zr isotopes. In this work, we demonstrate energy-dependent neutron capture cross section measurements for the first enriched sample to be measured: 90 Zr.

07 ISOTOPE AND RADIATION SOURCES↗

Bayesian Monte Carlo Evaluation of Imperfect (n, 233 U) Data and Model

Conventional nuclear data evaluation methods using generalized linear least squares make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate to minimize the cost function (even for nonlinear models); the model (e.g., of neutron cross section) and experimental data (including covariance data) are without defect and prior PDFs of parameters and measured data are known perfectly. Neglect of covariance between model parameters and measured data in conventional evaluations contributes to imperfections. These assumptions are inherent to the generalized linear least squares minimization method commonly used for resolved resonance region neutron cross section evaluations but are often not justified due to the presence of non-normal PDFs, nonlinear models (e.g., R-matrix formalism), and inherent imperfections in data and models (e.g., imperfect covariance data). Here, these assumptions are removed in a mathematical framework of Bayes’ theorem, which is implemented using the Metropolis-Hastings Monte Carlo method. Most importantly, new parameters are introduced to parameterize discrepancies between the theoretical model and measured data to quantify judgement about discrepancies or imperfections in a reproducible manner. An evaluation of 233U in the eV region using the ENDF-B/VIII.0 library and transmission data (Guber et al.) is presented, and posterior parameters are compared to those obtained by conventional evaluation methods. This example illustrates the effects of removing the most harmful assumption: that of model-data perfection.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nucleation of cubic boron nitride on boron-doped diamond via plasma enhanced chemical vapor deposition

Cubic boron nitride (c-BN), with a small 1.4% lattice mismatch with diamond, presents a heterostructure with multiple opportunities for electronic device applications. However, the formation of c-BN/diamond heterostructures has been limited by the tendency to form hexagonal BN at the interface. In this study, c-BN has been deposited on free standing polycrystalline and single crystal boron-doped diamond substrates via electron cyclotron resonance plasma enhanced chemical vapor deposition (ECR-PECVD), employing fluorine chemistry. In situ x-ray photoelectron spectroscopy (XPS) is used to characterize the nucleation and growth of boron nitride (BN) films as a function of hydrogen gas flow rates during deposition. The PECVD growth rate of BN was found to increase with increased hydrogen gas flow. In the absence of hydrogen gas flow, the BN layer was reduced in thickness or etched. The XPS results show that an excess of hydrogen gas significantly increases the percent of sp2 bonding, characteristic of hexagonal BN (h-BN), particularly during initial layer growth. Reducing the hydrogen flow, such that hydrogen gas is the limiting reactant, minimizes the sp2 bonding during the nucleation of BN. TEM results indicate the partial coverage of the diamond with thin epitaxial islands of c-BN. The limited hydrogen reaction is found to be a favorable growth environment for c-BN on boron-doped diamond.

Physics↗

Total thermal neutron cross section measurements of yttrium hydride from 0.0005 - 3 eV

Yttrium hydride serves as a neutron moderator material that enables compact, high temperature nuclear reactors. However, in order to accurately design and simulate a nuclear system relying upon yttrium hydride, the fundamental nuclear data of yttrium hydride must be well understood. Thermal neutron scattering law (TSL) evaluations represent an important aspect of nuclear data as thermal scattering can drastically alter the neutron multiplication factor of a system. Therefore, to support evaluation and validation of thermal neutron scattering for yttrium hydride, researchers at Rensselaer Polytechnic Institute (RPI) performed total thermal neutron cross section measurements for YH 1.68 and YH 1.85 over the energy range of 0.0005 - 3 eV. Further, these measurements represent the first total cross section measurements for yttrium hydride that encompass the entire thermal region. Comparisons were made against the ENDF-B/VIII.0, Zerkle & Holmes and Oak Ridge National Laboratory TSL evaluations, where generally good agreement was found.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Phenomenological R -Matrix parameterization of direct, doorway, and compound nuclear reactions [Abstract]

Although formal expressions for scattering matrix accounting for direct, doorway, and compound nuclear (CN) resonant reactions have been derived several decades ago in both the transition ( T -)matrix formalism and the reactance ( K -)matrix formalism, the absence of corresponding expressions in phenomenological R -matrix formalism has limited the application of the latter to CN resonant reactions only. We remove this limitation by parameterizing direct, doorway, and CN resonant reactions in a phenomenological R -matrix scattering matrix, and provide a parameterization for a corresponding Reich-Moore approximation of eliminated capture channels. Direct reactions induce (previously neglected) mixing among the incoming or outgoing R -matrix channel wave functions, parameterized by real and orthonormal channel-rotation matrix, M , whereby the original scattering matrix U is transformed into M T UM . Any real and orthonormal matrix, M , can be equivalently expressed as e η , where η is a real and skew-symmetric 2 rotation-generating matrix that subsequently yields a more intuitive parameterization of eliminated direct capture reactions in Reich-Moore approximation. A phenomenological R -matrix parameterization of doorway reactions is inferred by equating the expression for reactance ( K -)matrix, given in terms of Brune’s alternative R -matrix parameterization, to a corresponding expression derived using Feshbach’s projection operator formalism. Assuming that all doorway states, just like CN states, are confined within spheres defined by R -matrix channel radii, a new R -matrix-like term induced by doorway states is gleaned, wherein each doorway state is parameterized by its energy, width, and the strength of its coupling to each CN state. Since a Reich-Moore approximation for retained-channel scattering matrix ought to approximate the effect of eliminated capture channels taking place via direct, doorway, or CN reactions, each of the three kinds of reactions contributing to the capture entails a corresponding Reich-Moore parameterization in a first-order approximation: direct contribution is parameterized by introducing finite diagonal elements of a retained-channel rotation-generating matrix, doorway contribution is parameterized by doorway capture widths, while CN contribution is parameterized by conventional Reich-Moore capture widths. We will present evidence of direct and doorway reactions observed in recent measurements of resolved resonance cross sections at the Gaerttner LINAC Center at Rensselaer Polytechnic Institute, and will outline a path for implementing this new R -matrix parameterization into the SAMMY nuclear data evaluation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of ( n, γ ) cross section of 90 Zr [Abstract]

The isotopes of Zr with A = [90, 91, 92, 94] make up more than 97% of naturally occurring Zr and are important to many nuclear applications such as nuclear reactors. One of the attractive qualities of naturally occurring Zr isotopes is that they have a low σγ/σt ratio at most neutron energies, such that they improve the neutron economy in reactors by preferentially scattering neutrons rather than absorbing them. This same quality also presents a challenge to measuring the capture cross section, σγ, of Zr isotopes. The ENDF/B-VIII.0 library has a relative uncertainty of approximately 10-20% for incident neutron energies < 0.1 MeV, and uncertainty greater than 20% for energies > 0.1 MeV for the majority of natural Zr isotopes. This motivated the Nuclear Criticality Safety Program (NCSP) to embark on a campaign to accurately measure and evaluate these isotopes of Zr. Here we will show energy-dependent neutron capture cross section measurements for the first enriched sample to be measured: 90 Zr. The measurements of isotopically enriched samples are being carried out at the Geel Electron Linear Accelerator (GELINA) facility of the Joint Research Center - Geel (JRC-Geel) of the European Union. As isotopic enrichment is a costly process we are careful not to activate any of the samples, as this may hinder future radiation-sensitive measurements. The activation analysis is presented in a report by Brown et al.. Once we were satisfied that the Zr samples would not be activated by the measurements, the 90 Zr sample was fabricated at Oak Ridge National Laboratory (ORNL) and shipped to GELINA. The dimensions of the cylindrical sample are approximately 0.12 cm thick and a radius of 2.5 cm. Since 90 Zr is not chemically reactive to air, bare metallic samples were employed. The sample was measured at a flight path (FP) length of 60 m, using four C 6 D 6 detectors on FP14. The final paper will include experimental details and measured cross section data for 90 Zr compared to current evaluated nuclear data libraries.

07 ISOTOPE AND RADIATION SOURCES↗

Modernization efforts for the R-Matrix code SAMMY [Slides]

This presentation details modernization efforts for the R-Matrix code SAMMY. It provides an overview of SAMMY, modernization and maintenance goals, modernization accomplishments, covariance data, SAMMY API development, as well as future plans.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modernization efforts for the R -Matrix code SAMMY [Abstract]

The R-Matrix code SAMMY is a widely used nuclear data evaluation code focused on the resolved range, which includes corrections for experimental effects. The code is still mostly written in Fortran 77, and uses a memory management system suitable for the time of its initial writing (1984). A modernization effort is under way to bring the code in-line with modern software development practices. A continuous-integration testing framework was added, automating the large existing set of test cases. It is run on every commit. The memory management was updated to current standard practices suitable for modern software analysis tools. The code can be obtained from https://code.ornl.gov/RNSD/SAMMY. The resonance parameters and covariance information are now stored in C++ objects shared by SAMMY and AMPX, the processing code that generates nuclear data libraries for SCALE. This allows for easier maintenance and access to the resonance parameters inside and outside of SAMMY. This feature is already used by accessing and changing parameters in memory in the Bayesian Monte Carlo Evaluation Framework for Cross Sections Nuclear Data and Integral Benchmark Experiments project, Further plans include the switch to the ENDF reading and writing routines in AMPX, as these routines are more robust, easier to maintain, and support more features. Of note here is support for the new GNDS format. Previously it wasn’t easy to share the full covariance matrix for evaluations containing more than one isotope due to limitations on the ENDF format; this is now supported in GNDS. The data are currently available in a binary SAMMY format and can be exported to GNDS to make them more widely available and sharable. The next step will be to use the same resonance processing code at 0K in AMPX and SAMMY as one of the available Reich-Moore R-Matrix formalism. The first step toward this goal is to isolate the reconstruction into a module that takes resonance parameters as its input and does not depend on SAMMY global parameters. This goal has been achieved and it should now be possible to more easily change the resonance formalism and add enhancements as the Phenomenological R-Matrix parameterization of direct, doorway, and compound nuclear reactions discussed elsewhere on this conference. This concerted modernization and enhancement effort provides multiple advantages to the nuclear data community. It will allow parameter optimization using enhanced formalisms, including experimental effects, that better match complex experimental data. Then those evaluated parameters can immediately be passed off to AMPX to be reconstructed with the exact same cross section model and be put into a data library for subsequent testing using SCALE and the Valid Benchmark suite or other suitable benchmark suites.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Bayesian Monte Carlo Evaluation Framework for Imperfect Data and Models [Abstract]

Nuclear data evaluation methods conventionally make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate for minimization of a cost function (even for non-linear models); and that both the model (of, e.g., neutron cross section) and experimental data (including their covariance data) are perfect. These assumptions are inherent to the well-known generalized linear least squares (GLLS) minimization method commonly used for evaluations of resolved resonance region (RRR) neutron cross sections. However, these assumptions are often not justified due to the presence of non-normal PDFs, non-linear models (e.g. R -matrix formalism), and inherent imperfections in data and models (e.g. discrepant data sets, discrepancies between the previous evaluation and newly measured data, or imperfect covariance data). We remove the said assumptions in a mathematical framework of Bayes’ theorem, and implement it using the Metropolis-Hastings Monte Carlo method. Parameters of a new kind are introduced to parameterize inherent imperfections, e.g. , any discrepancies between the theoretical model and measured data. These new parameters enable evaluators to quantify their expert judgement about any discrepancies or imperfections in a reproducible manner. We demonstrate the framework with an ongoing evaluation of 233 U in the eV region using the ENDF-B/VIII library and transmission data measured by Guber, et al. , and compare the posterior parameters to those obtained by conventional evaluation methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Bayesian Optimization Framework for Imperfect Data or Models

Conventional Bayesian optimization methods implicitly assume that the data and model being optimized are “perfect.” This assumption leads to inaccurate posterior probability distribution functions (PDFs) when applied to “imperfect” data or models. The new Bayesian optimization framework presented in this report provides a way to parameterize the effect of imperfections usually encountered in a prior PDF of generalized data or a model on the posterior PDF. The effects of imperfections are parameterized by a set of constraints imposed on the posterior expectation values of deviations between the data and the model and on their covariance matrix elements. A particular set of values for these constraints conveys an evaluator’s best estimate of the effect of imperfections on the corresponding posterior expectation values. When a prior PDF of generalized data is assumed to be normal, an expression for a posterior PDF satisfying an arbitrary set of constraints is derived analytically for linear models. An analogous iterative algorithm is given for nonlinear models. The corresponding posterior PDF should be used to estimate any posterior expectation values in the presence of imperfections parameterized by that set of constraints. A posterior PDF of a conventional Bayesian optimization method is recovered analytically when all evaluator-specified constraints are set to zero (i.e., in the absence of any imperfections). The analytical expressions derived in this report for normal PDFs and linear models were verified numerically by a Metropolis–Hastings Monte Carlo method. The methods presented herein could be applied to any kind of data or models, including differential cross-section data or integral benchmark experiments.

97 MATHEMATICS AND COMPUTING↗