Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mathematical physics”

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 865 records · Page 48

Impact of Increased Latent Generations on Sensitivity Calculations with SCALE

Analyses of cross section sensitivity data from systems with fissile material allow analysts to associate an importance for each material, nuclide, reaction, and neutron energy by simulating real world criticality scenarios. Although criticality safety validation efforts can be guided by the cross-section sensitivity and uncertainty data generated for a particular system, these calculations can often be computationally expensive and sometimes cumbersome without proper guidance. The TSUNAMI suite within the SCALE code package has several methods for generating sensitivity data, including multigroup and continuous energy (CE) capabilities. The release of SCALE 6.3 has three different CE methods for generating cross section sensitivity data: (1) the Iterated Fission Probability (IFP) method with the KENO Monte Carlo transport solver, (2) the IFP method with the Shift Monte Carlo transport solver, and (3) the Contributon-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method with the KENO Monte Carlo transport solver. Although the CLUTCH method has additional parameters for generating sensitivity data files relative to the IFP method, all three methods use latent generations, which are the generations between an event (i.e., fission) and the assessment of importance based on the asymptotic population of progeny neutrons. Increasing the number of latent generations in a calculation leads to increased discrimination of the sensitivity coefficients but at the cost of the increased uncertainty associated with those generated values. Analysts must balance the accuracy of the sensitivity calculations and its uncertainty with the associated computational cost involved in generating the values. This paper discusses the impact of adjusting the latent generation parameter for a range of sensitivity values and how these changes compare with the direct perturbation values obtained from a change of ±0.5% Δ k in both benchmark and safety application models. Two benchmarks from the International Handbook of Evaluated Criticality Safety Benchmark Experiments and the MPC-32 dual purpose canister for spent nuclear fuel are used for analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

Impact of Thermal Scattering Law on Similarity Assessment in Light-Water or Polyethylene-Moderated Systems

A collaborative effort between Pacific Northwest National Laboratory (PNNL) and Oak Ridge National Laboratory (ORNL) is underway to provide a technical basis and methodology for the criticality safety community to use the sum of fractions (SoF) method for generating limits for mixtures of “selected actinide nuclides” included in the ANSI/ANS 8.15 standard. The PNNL scope in this project is to define a range of mixtures of 233 U, 235 U, and 239 Pu moderated with either light water or polyethylene and to examine the critical masses for these mixtures. The ORNL scope is primarily to provide validation support for these studies. More complete discussion of the project and its validation aspects will be presented at the upcoming International Conference on Nuclear Criticality Safety (ICNC) this Fall in Sendai, Japan. Clear differences in benchmark similarity to application systems as assessed by the integral parameter ck are noted in the validation studies performed as part of this project as a function of moderator. The c k value is a correlation coefficient that represents that amount of shared uncertainty in k eff due to cross sections between two systems. Individual nuclide-reaction contributions between the two systems can be simply summed to arrive at the total c k value. Specifically, the c k values for light-water–moderated solution experiments are higher for a water-moderated application than for a polyethylene-moderated application. This result is neither totally unexpected nor surprising, but the magnitude of the difference was difficult to anticipate. The TSUNAMI sequence, in the SCALE 6.2.4 code package developed by ORNL, was used to generate eigenvalues and reactivity effects with perturbation-theory based approach through sensitivity coefficients for all nuclides in the system with all reactions and energy groups. The TSUNAMI-Indices and Parameters (IP) sequence then uses the sensitivity data generated through TSUNAMI to generate relational parameters (i.e., c k ) to determine the degree of similarity between systems. One detail of the SCALE material and data implementation must be discussed at this point. Several thermal scattering laws (TSLs) are available for 1 H. SCALE uses a different nuclide ID number for each TSL; essentially, each version of 1 H is treated as a unique nuclide. For example, 1 H bound in water ( 1 H-H2O) is assigned the nuclide ID 1001, whereas 1 H bound in polyethylene (h-poly) is assigned the nuclide ID 9001001. The same cross section data are used for all reactions in 1 H, regardless of TSL, except for scattering below the TSL cutoff energy. TSUNAMI-IP treats different nuclide IDs as different nuclides; thus, no uncertainty is shared between 1 H-H2O and h-poly, despite much of the same data, including covariance data, being used for both nuclides. This presents a question: how much of the difference in assessed similarity between water- and polyethylene-moderated systems is due to the moderators, and how much is caused by the treatment of 1 H-H2O and h poly with cross section and covariance data. The extended edits generated by TSUANMI-IP allow for an investigation of this issue specifically, as well as a demonstration of the general techniques available within TSUNAMI to understand the results of the similarity assessment. This paper presents and analyzes the similarity assessment of both water- and polyethylene-moderated systems for a single benchmark: PST-002-001.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Explaining machine-learning models for gamma-ray detection and identification

As more complex predictive models are used for gamma-ray spectral analysis, methods are needed to probe and understand their predictions and behavior. Recent work has begun to bring the latest techniques from the field of Explainable Artificial Intelligence (XAI) into the applications of gamma-ray spectroscopy, including the introduction of gradient-based methods like saliency mapping and Gradient-weighted Class Activation Mapping (Grad-CAM), and black box methods like Local Interpretable Model-agnostic Explanations (LIME) and SHapley Additive exPlanations (SHAP). In addition, new sources of synthetic radiological data are becoming available, and these new data sets present opportunities to train models using more data than ever before. In this work, we use a neural network model trained on synthetic NaI(Tl) urban search data to compare some of these explanation methods and identify modifications that need to be applied to adapt the methods to gamma-ray spectral data. We find that the black box methods LIME and SHAP are especially accurate in their results, and recommend SHAP since it requires little hyperparameter tuning. We also propose and demonstrate a technique for generating counterfactual explanations using orthogonal projections of LIME and SHAP explanations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

(U) Transport Corrections Implemented in SENSMG

Truncating the spherical harmonics expansion of the neutron scattering source in the Boltzmann transport equation leads to truncation error. Methods for approximately correcting for this truncation error are called transport corrections. The PARTISN multigroup discrete ordinates neutron transport code has three transport correction options. These options can be specified in the SENSMG multigroup neutron sensitivity code and passed to PARTISN for the neutron transport.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An Analytic Benchmark for the Solution to the Isotopic Fission Spectrum Mixture Problem

In neutron transport calculations, it is common to specify material compositions in terms of constituent isotopes. Material compositions may be described by isotope number densities and associated microscopic cross sections. For general reaction cross sections, the macroscopic cross sections of a composition are simply the summation of the sum of the products of isotope number densities and microscopic cross sections. A notable exception is the mixture of the neutron fission spectrum in fissile material. To demonstrate proper and improper mixture of the neutron fission spectrum, a zero-dimensional test problem is developed. Using the test problem, it is demonstrated that the improper mixing of the fission spectrum results in an eigenvalue error of approximately 65 pcm. Though this may seem small, any error in such a simple calculation is unacceptable. In similar infinite-homogeneous test problems, eigenvalue errors of more than 300 pcm have been observed. It is also shown that the error due to fission spectrum mixing can be exaggerated for computer program verification purposes. It is concluded that the proper mixture of the neutron fission spectrum is essential for accurate neutron transport simulations. The effect of improper neutron fission spectrum mixture is demonstrated and quantified and this test problem may be used for verification purposes in the future.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Online and Offline Monitoring of the ICARUS Detector

In preparation for the launch of the ICARUS -T600 far detector in the Short Baseline Neutrino program, (SBN), it is important to monitor the performance of the detector and quality of the data in real-time. Particularly important to monitor are noise levels of the TPC readout electronics, and the purity of the liquid argon in the TPC, as these can affect the overall sensitivity of the detector to detecting and reconstructing neutrino interactions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

PAGOSA Mesh Block Calculator

A method for automatically generating the “mesh namelist block” portion of a PAGOSA input file, which specifies the boundaries and zone sizes of the computational mesh, has been developed. This application standardizes the mesh block for 2-, and 3-dimensional problems and is computed after the user specifies the extent and resolution of the mesh and the number of processors which will be used to run the problem. The inputs may be provided directly through a series of command-line guided questions, or a separate input file. The mesh block is written to both standard output and to disk.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

(U) LNK3DNT Capability in SENSMG Tested with Transport Correction

The SENSMG multigroup neutron sensitivity code can read materials and geometry from redoin and lnk3dnt files written by the PARTISN multigroup discrete ordinates code or the MCNP6 Monte Carlo code. Reference 5 described an initial version of this capability. Reference 6 introduced two important extensions: 1) Reaction-rate ratios are now allowed as responses with redoin/lnk3dnt inputs; and 2) Sensitivities of reaction-rate ratios to mass densities and interface locations are now calculated (for all problems, not just redoin/lnk3dnt problems). In this report, we document the testing of the redoin/lnk3dnt capability when used with transport correction. PARTISN’s transport corrections are described in Ref. 3, and Refs. 7 and 8 discuss the calculation of cross-section sensitivities when transport corrections are applied. References 7 and 8 verify SENSMG for problems with transport correction using central differences. This report verifies the capability with redoin/lnk3dnt inputs by comparison with regular SENSMG inputs. Unlike Refs. 5 and 6, in this report we use a spherical test problem because the sensitivity of the response with respect to the highest-order scattering is zero in cylinders due to the use of flux moments in the inner products.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Release of NDI Libraries MENDF80 and MT80—Libraries Based on ENDF/B-VIII.0

Several new NDI libraries have been released by the Nuclear Data Team. There are multigroup-energy neutron libraries (see Section II.) and associated production-depletion libraries (see Section IV.). The libraries are based on data from ENDF/B-VIII.0. The libraries will be available with NDI version 2.1.4 when it becomes available in Spring of 2020. These libraries are the libraries recommended for use by the Nuclear Data Team. We encourage all to adopt these libraries at their earliest convenience.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

CEA/DAM - CEA/DEN - LANL - IAEA - IPHC Collaboration Meeting on Nuclear Data, Nuclear Reaction Theories, and Machine Learning (Summary Report)

Dupuis hosted our biennial CEA/NNSA collaboration meeting, and Kawano and Lovell of LANL and Tamagno of CEA/DEN (Cadarache) visited CEA/DAM to discuss current issues and further developments in nuclear reaction modeling, nuclear data evaluation, and relevant programs of modern computer science, especially the machine learning technology. The discussions included topics of our common interest such as neutron inelastic scattering, capture, and fission, with particular emphasis on the actinide nuclear data evaluations. We report briefly the summary of our discussions. We also had a half-day video meeting with Capote of IAEA and Kerveno of IPHC to discuss nuclear reactions on actinides.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Annular Core Research Reactor (ACRR) Pulse Curve Characterization

Reactor pulse characterization at the Annular Core Research Reactor (ACRR) at Sandia Technical Area V (TA-V) is commonly done through photo conductive detection (PCD) and calorimeter detectors. Each of these offer a mode of analyzing a digital signal with different advantages/methods for determination of integrated dose or temporal based metrology. This report outlines a method and code that takes the millions of data points from such detectors and delivers a characteristic pulse trendline through two main methods: digital signal filtration and machine learning, in particular, Support Vector Machines (SVMs). Each method's endpoint is to deliver a characteristic curve for the many bucket environments of ACRR while considering other points of interest including delayed gamma removal for prompt dose metrics. This work draws and adds on previous work detailing the delayed gamma fraction contributions from CINDER simulations of the ACRR. Results from this project show a method to determine characteristic curves in a way that has previously been limited by data set size.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Theoretical Review of the MFIX Fluid and Two-Fluid Models

MFIX (Multiphase Flow with Interphase eXchanges) is a general purpose code that can be used for describing the hydrodynamics, chemical reactions, and heat transfer of dense or dilute multiphase fluid flows. MFIX calculations give detailed information on pressure, temperature, composition, and velocity distributions in the defined system. This report provides an overview of MFIX’s expressions for single phase flow and continuous dispersed multiphase flow (i.e., the two-fluid model). In addition to the two fluid model (MFIX-TFM), MFIX contains a discrete element model (MFIX-DEM) and a multi-phase particle in cell model (MFIX-PIC) that are not reviewed within this document.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SCALE Lattice Physics Code Assessments of Accident Tolerant Fuel

This report highlights accident-tolerant fuel (ATF) code assessment activities performed under Project NRCHQ- 60-17-T-0017, Lattice Physics Enhancements and Assessment. ATF covers a broad range of advanced fuel and clad designs for light water reactors (LWRs) to enhance performance under several accident conditions. Several ATF concepts are anticipated to be deployed as lead test rods (LTRs) or lead test assemblies (LTAs) within the next five years. The purpose of this work is to assess the predictive capabilities of NRC neutronics codes that underpin various licensing calculations. ATF designs use different fuel and clad materials compared to standard UO 2 fuel and zirconium-alloy claddings (hereafter UO 2 -Zry). These new materials and geometrical designs need to be assessed to quantify the impact of nuclear data uncertainties on quantities of interest (QOIs) in licensing calculations and the impact of modeling approximations which may be valid for UO 2 -Zry but not for ATF. This report outlines a systematic approach for ATF neutronics code assessment which includes sensitivity and uncertainty (S/U) analysis of nuclear data, identification of experimental benchmark and gaps for code validation, investigation of modeling approximations, and code-to-code comparisons of calculated QOIs against high-fidelity reference continuous energy (CE) Monte Carlo (MC) calculations. This report focuses on the assessment of the SCALE/Polaris lattice physics code for reactor safety analysis. Polaris lattice physics calculations generate few group (FG) cross sections for PARCS full-core calculations. (Full-core analysis, spent fuel analysis, and severe accident analysis will be performed in future work.) The selected ATF concepts for this report include Cr 2 O 3 and Al 2 O 3 -Cr 2 O 3 -doped UO 2 fuel, U 3 Si 2 fuel, FeCrAl cladding, SiC cladding, and Cr-coated cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

COMET/ZEUS Activation Foils Report for Irradiations IER-241 and IER-261

The goal of this work was to determine the internal neutron spectrum of the COMET/ZEUS assembly using activation foils. Similar measurements were performed on the COMET assembly when it was located at the Los Alamos Critical Experiments Facility (LACEF) at the Los Alamos National Laboratory (LANL). The present measurements were performed after the COMET assembly was moved to the National Criticality Experiments Research Center (NCERC) located in the Device Assembly Facility (DAF) at the Nevada National Security Site (NNSS). The December 2012 COMET/ZEUS irradiation was the first in a series of experiments, and the COMET/ZEUS assembly was selected because of operational readiness and authorization to run specific targets. The August 2013 irradiation was the second in the series and involved only small changes in the configuration. This report describes the unfolding of the activation foils to determine the internal neutron spectra.

07 ISOTOPE AND RADIATION SOURCES↗