Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “USe”

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 631 records · Page 35

High broadband photoconductivity of few-layered MoS 2 field-effect transistors measured using multi-terminal methods: effects of contact resistance

Among the layered two dimensional semiconductors, molybdenum disulfide (MoS 2 ) is considered to be an excellent candidate for applications in optoelectronics and integrated circuits due to its layer-dependent tunable bandgap in the visible region, high ON/OFF current ratio in field-effect transistors (FET) and strong light–matter interaction properties. In this study, using multi-terminal measurements, we report high broadband photocurrent response ( R ) and external quantum efficiency (EQE) of few-atomic layered MoS 2 phototransistors fabricated on a SiO 2 dielectric substrate and encapsulated with a thin transparent polymer film of Cytop. The photocurrent response was measured using a white light source as well as a monochromatic light of wavelength λ = 400 nm–900 nm. We measured responsivity using a 2-terminal configuration as high as R = 1 × 10 3 A W -1 under white light illumination with an optical power P opt = 0.02 nW. The R value increased to 3.5 × 10 3 A W -1 when measured using a 4-terminal configuration. Using monochromatic light on the same device, the measured values of R were 10 3 and 6 × 10 3 A W -1 under illumination of λ = 400 nm when measured using 2- and 4-terminal methods, respectively. The highest EQE values obtained using λ = 400 nm were 10 5 % and 10 6 % measured using 2- and 4-terminal configurations, respectively. The wavelength dependent responsivity decreased from 400 nm to the near-IR region at 900 nm. The observed photoresponse, photocurrent–dark current ratio (PDCR), detectivity as a function of applied gate voltage, optical power, contact resistances and wavelength were measured and are discussed in detail. The observed responsivity is also thoroughly studied as a function of contact resistance of the device.

36 MATERIALS SCIENCE↗

Using active learning to improve quasar identification for the DESI spectra processing pipeline

The Dark Energy Spectroscopic Instrument (DESI) survey uses an automatic spectral classification pipeline to classify spectra. QuasarNET is a convolutional neural network used as part of this pipeline originally trained using data from the Baryon Oscillation Spectroscopic Survey (BOSS). In this paper we implement an active learning algorithm to optimally select spectra to use for training a new version of the QuasarNET weights file using only DESI data, with the goal of improving classification accuracy. This active learning algorithm includes a novel outlier rejection step using a Self-Organizing Map to ensure we label spectra representative of the larger quasar sample observed in DESI. We perform two iterations of the active learning pipeline, assembling a final dataset of 5600 labeled spectra, a small subset of the approximately 1.3 million quasar targets in DESI's Data Release 1. When splitting the spectra into training and validation subsets we achieve similar performance to the previously trained weights file in completeness and purity calculated on the validation dataset but do so with less than one tenth of the amount of training data. The new weights also more consistently classify objects in the same way when used on unlabeled data compared to the old weights file. In the process of improving QuasarNET's classification accuracy we discovered a systemic error in QuasarNET's redshift estimation and used our findings to improve our understanding of QuasarNET's redshifts.

Machine learning↗

Measuring impacts of California agri-environmental programs using field-scale satellite data

In the past decade, California has invested over $\$$200 million in direct grants to growers to support the adoption of agricultural practices that save water and/or improve soil health while also reducing greenhouse gas emissions. Ex-post evaluation of agri-environmental outcomes of these grant programs, however, is limited. We use satellite data to monitor changes in field-level consumptive water use and greenness (i.e. normalized difference vegetation index), a proxy for agricultural productivity, for the most frequently funded crop-types (almonds, grapes, and walnuts) in two California Department of Food and Agriculture programs. Nearly 600 fields receiving funding during the 2014–2022 period were analyzed using two causal inference methods. Fields that received grants to both upgrade irrigation systems and install irrigation water management sensors showed reduced consumptive water use and greenness by an average of 3.5% and 4.2%, respectively (significant at the 10% level). In contrast, we find that the adoption of only irrigation water management sensors, which are designed to inform irrigation scheduling and management, resulted in an average increase of 4.1% and 4.8% in consumptive water use and greenness respectively (significant at the 5% level). We find negligible effects for either consumptive water use or greenness when both pump efficiency upgrades and sensors were implemented. We further find that grants for compost addition and cover cropping led to small greenness increases of 1.7% and 2.8% respectively (significant at the 10% level) and had insignificant effects on consumptive water use. Our analysis of five agri-environmental program interventions reveals that several practice outcomes may be at odds with stated program goals of reducing water use while maintaining or improving agricultural productivity.

agriculture↗

Past agricultural land use affects multiple facets of ungulate antipredator behavior

Antipredator behavior affects prey fitness, prey demography, and the strength of ecological interactions. Although predator-prey interactions increasingly occur in habitats that experience multiple forms of human-generated disturbance, it is unclear how different forms of disturbance might affect antipredator behavior. Fire is a contemporary disturbance that has dramatic effects on terrestrial habitats. Such habitats may have also experienced past disturbances, like agricultural land use, that leave lasting legacies on habitat structure (e.g., overstory and understory composition). It is unclear how these past and present disturbances affect the use of different antipredator behaviors, like temporal avoidance and vigilance. We examined whether variation in disturbance regimes generates differences in ungulate antipredator behavior by using cameras to measure white-tailed deer vigilance and activity time across 24 longleaf pine woodlands that vary in past land use and contemporary fire regime. Regardless of land-use history, woodlands with high fire frequencies had 4 times less vegetation cover than low-fire woodlands, generating riskier habitats for deer; however, deer responded to fire with different antipredator strategies depending on land-use history. In nonagricultural woodlands, fire affected deer activity time such that activity was nocturnal in low-fire woodlands and crepuscular in high-fire woodlands. In post-agricultural woodlands, fire affected vigilance and not activity time such that deer were more vigilant in high-fire woodlands than in low-fire woodlands. Furthermore, these results suggest that ungulate antipredator behavior may vary spatially depending on past land use and contemporary fire regime, and such disturbances may generate “landscapes of fear” that persist for decades after agricultural use.

60 APPLIED LIFE SCIENCES↗

Conversion of Site-Specific Meteorological Data for use in CAP-88 PC

CAP-88 PC is a commonly used radiological atmospheric dispersion model. This US EPA-approved model is used to demonstrate compliance with atmospheric emission regulations for radionuclides. While the model includes a large library of meteorological data for use across the United States, there are applications when users may wish to use onsite meteorological data as an input to the CAP-88 PC model. Here, in this work, we present a work-around process for preparing and converting onsite data for use in CAP-88 PC. However, the use of local data should provide a more realistic estimate of doses to members of the public in the immediate vicinity of a facility, although the regulatory agency having jurisdiction may not accept the use of local data for compliance. Additionally, the historical meteorological records from 20+ years ago (at a site many kilometers away) might not be representative of current local weather patterns, highlighting another benefit of using local meteorological data.

54 ENVIRONMENTAL SCIENCES↗

Characterizing the nonequilibrium response of FeRh thin films using time-domain thermoreflectance

Time-domain thermoreflectance (TDTR) characterization of FeRh throughout its first-order antiferromagnetic (AF) to ferromagnetic (FM) transition shows that the transient reflectance ΔR⁡(t)/R strongly depends on the magnetic order of the sample. Using TDTR, which uses optical pulses to induce small temperature excursions, we have found that Δ⁢R⁡(t)/R of the AF phase exhibits a large negative response, while the response of the FM phase is positive. This magnetic phase sensitivity has allowed us to study the transient response of both the AF and FM phases to the pump-pulse excitation and the mixed phase of the material. These results are significant since the ultrafast properties of antiferromagnetic materials and mixed antiferromagnetic and ferromagnetic materials are difficult to detect using other conventional techniques. We have found that the AF phase exhibits a strong subpicosecond decaying signal not observed in the FM phase. The magnetic phase dependence of the sign of Δ⁢R⁡(t)/R is qualitatively explained using the results of ab initio density functional theory calculations. Using the two-temperature model, we found that the change in the thermalization time across the transition is caused by differences in both the electronic heat capacity and the electron-phonon coupling factor of the AF and FM phases. The electron-phonon coupling constant in the AF phase is also determined using the two-temperature model conducted using the ntmpy code package. For the FM phase, we provide boundaries for the magnitude of the electron-phonon coupling factor for the FM phase. Finally, these results indicate that TDTR can be used to study the transient properties of magnetic materials that are otherwise challenging to probe.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Weakly Supervised Event Classification Using Imperfect Real-world PMU Data with Scarce Labels

This paper studies event classification using imperfect real-world phasor measurement unit (PMU) data with scarce event types (labels). By investigating the real-world PMU data, it is observed that most real-world PMU data's event type is unknown, which makes it challenging to directly use such dataset to build event classifiers as existing classification techniques require high-quality training data with known event type (i.e., label). To address this challenge, a weakly supervised learning based event classification approach is developed, which can use noisy and low-quality PMU data for the training. First, data quality issues are fixed using data preprocessing techniques and then event features are constructed from the PMU data. Using these features, a series of labeling functions are learnt to generate initial estimates of the labels of large amounts of unlabeled PMU data. As the labeling functions are learnt using the same data with scarce labels, the label estimates from the labeling functions can be correlated, noisy, and bias. To enhance these initial estimates, a generative model is developed to characterize the dependencies among the estimated labels, based on which better labels are obtained for training event classifiers. Numerical experiments using the real-world dataset from the Western Interconnection of the U.S. power transmission grid show that the proposed weakly supervised event classifier trained using the dataset with only 5% labeled data can achieve 78.4% classification accuracy.

Liu, Yunchuan↗

Spatially Local Surrogate Modeling of Subgrid-Scale Effects in Idealized Atmospheric Flows: A Deep Learned Approach Using High-Resolution Simulation Data

Abstract We introduce a machine learned surrogate model from high-resolution simulation data to capture the subgrid-scale effects in dry, stratified atmospheric flows. We use deep neural networks (NNs) to model the spatially local state differences between a coarse-resolution simulation and a high-resolution simulation. The setup enables the capture of both dissipative and antidissipative effects in the state differences. The NN model is able to accurately capture the state differences in offline tests outside the training regime. In online tests intended for production use, the NN-coupled coarse simulation has higher accuracy over a significant period of time compared to the coarse-resolution simulation without any correction. We provide evidence of the capability of the NN model to accurately capture high-gradient regions in the flow field. With the accumulation of the errors, the NN-coupled simulation becomes computationally unstable after approximately 90 coarse simulation time steps. Insights gained from these surrogate models further pave the way for formulating stable, complex, physics-based spatially local NN models which are driven by traditional subgrid-scale turbulence closure models. Significance Statement Flows in the atmosphere are highly chaotic and turbulent, comprising flow structures of broad scales. For effective computational modeling of atmospheric flows, the effects of the small- and large-scale structures need to be captured by the simulations. Capturing the small-scale structures requires fine-resolution simulations. Even with the current state-of-the-art supercomputers, it can be prohibitively expensive to simulate these flows when computed for the entire earth over climate time scales. Thus, it is necessary to focus on the larger-scale structures using a coarse-resolution simulation while capturing the effects of the smaller-scale structures using some parameterization (approximation) scheme and incorporating it into the coarse-resolution simulation. We use machine learning to model the effects of the small-scale structures (subgrid-scale effects) in atmospheric flows. Data from a fine-resolution simulation is used to compute the missing subgrid-scale effects in coarse-resolution simulations. We then use machine learning models to approximate these differences between the coarse- and fine-resolution simulations. We see improved accuracy for the coarse-resolution simulations when corrected using these machine learned models.

54 ENVIRONMENTAL SCIENCES↗

Contraceptive Sabotage and Contraceptive Use at the Time of Pregnancy: An Analysis of People with a Recent Live Birth in the United States

Contraceptive sabotage and other forms of intimate partner violence (IPV) can interfere with contraceptive use. We used 2012 to 2015 Pregnancy Risk Assessment Monitoring System data from 8,981 people residing in five states who reported that when they became pregnant, they were not trying to get pregnant. We assessed the relationships between ever experiencing contraceptive sabotage and physical IPV 12 months before pregnancy (both by the current partner) and contraceptive use at the time of pregnancy using multivariable logistic regression. We also assessed the joint associations between physical IPV 12 months before pregnancy and ever experienced contraceptive sabotage with contraceptive use at the time of pregnancy. Few people ever experienced contraceptive sabotage (1.8%; 95% confidence interval [CI]: 1.4, 2.3) or physical IPV 12 months before pregnancy (2.8%; 95% CI: 2.3, 3.3). In models adjusted for age, race/ethnicity, marital status, education, and state of residence, ever experiencing contraceptive sabotage was associated with contraceptive use at the time of pregnancy (adjusted odds ratio [aOR]: 1.73; 95% CI: 1.06, 2.82), but not with physical IPV 12 months before pregnancy (aOR: 0.69; 95% CI: 0.46, 1.02). When examining the joint association, compared to not ever experiencing contraceptive sabotage or physical IPV 12 months before pregnancy, ever experiencing contraceptive sabotage was significantly related to contraceptive use at the time of pregnancy (aOR: 1.72; 95% CI: 1.00, 2.95). However, it was not associated with experiencing physical IPV 12 months before pregnancy (aOR: 0.68; 95% CI: 0.45, 1.04) or with experiencing both contraceptive sabotage and physical IPV 12 months before pregnancy (aOR: 1.21; 95% CI: 0.42, 3.50), compared to not ever experiencing contraceptive sabotage or physical IPV 12 months before pregnancy. Our study highlights that current partner contraceptive sabotage may motivate those not trying to get pregnant to use contraception; however, all people in our sample still experienced a pregnancy.

Huber-Krum, Sarah↗

Selected Uses of TSUNAMI in Critical Experiment Design and Analysis

Validation in criticality safety is performed by comparing the results of critical experiments with the calculated results from models of the experiments using the computational method to be validated. Laboratory critical experiments are controlled systems that achieve a k eff of approximately 1 and enable investigation of the parameters at which such a critical condition is achieved. For the critical experiments used in a validation to capture the biases of the materials and neutron energy spectra of interest, those materials must be included in the experiment such that they influence k eff or another observable parameter with statistical significance. This paper discusses the use of sensitivity uncertainty (S/U) methods to develop critical experiments for various purposes. S/U techniques are useful for understanding the underlying components of nuclear data which affect the k eff or another parameter of a given configuration. S/U calculations are most commonly used to compare existing experiments to applications of interest; however, S/U techniques can also be used to identify, optimize, or assess features of proposed experiments so that they can better test specific portions of nuclear data or match an application of interest. The S/U techniques discussed here are from the TSUNAMI code system. The two primary codes discussed in this work are TSUNAMI-3D, which implements the KENO criticality code to calculate the sensitivity of k eff to nuclear data, and TSAR, which calculates the sensitivity of a reactivity difference between two configurations based on their TSUNAMI-3D generated sensitivity profiles. The methods used in these tools are discussed in more detail in the SCALE manual. This paper is one of a series on the development and use of TSUNAMI tools. The other papers address development of TSUNAMI methods and a review of TSUNAMI applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗

Slag Management of Carbon Feedstock Used in Gasification FY17 TCF between Eastman Chemical and NETL (Final Report)

NETL, in collaboration with Eastman Chemical Corp (Kingsport, TN), targeted the joint evaluation and refinement of a slag management program developed by NETL (NETL IP, patent applied for) to control gasifier slag properties formed from ash in individual or combinations of different coal-based carbon feedstock materials and specific slag modifiers used in gasification. Slag forms from mineral or organic-metallic impurities in the carbon feedstock that melt at the high temperatures of gasification, coalescing to form a liquid (slag) that must flow from the gasifier. Slag additives, also known as slag modifiers, are often added to the carbon feedstock to modify the ash chemistry and achieve desired slag flow properties at target temperatures used in gasification. If slag does not flow from a gasifier, its buildup within the gasification chamber will cause either gasifier shutdown to physically remove the slag (consuming time and damaging the protective refractory lining), or the gasification temperature must be increased to lower slag viscosity so it can flow from the gasifier, leading to greater refractory wear/corrosion, along with increased fouling and downstream component wear. Because of these issues, gasification is typically targeted at the lowest temperature possible where slag will flow in a controlled manner (a temperature range called T100 for the lowest temperature, with T250 being the highest temperature where the slag has achieved the desired fluidity). That temperature range also needs to yield a desired syngas composition, highlighting the importance of knowing and being able to control carbon feedstock slag viscosity and the temperature of thermal breakdown of the carbon feedstock. Slag viscosity varies with changes in the gasification chamber temperature, oxygen partial pressure, and slag chemistry. The slag model work with Eastman Chemical targeted predicting and controlling molten coal ash slag behavior at specific gasification conditions. The final goal of the refined model was to develop/refine the NETL slag model so a gasifier operator would be able to predict slag behavior prior to purchasing or using a carbon feedstock and control slag properties once a carbon feedstock is being gasified, allowing a gasifier user to determine what mineral additives or blends of different carbon materials would be needed to achieve desired slag properties at targeted gasification conditions. For a given carbon feedstock and desired gasification condition, the actual gasification cost and slag behavior would be able to be predicted using the slag management program in advance of purchasing or using the carbon feedstock – allowing accurate control of slag ash behavior in real time. As part of the TCF project work with Eastman Chemical, the viscosity of several ash compositions of interest to them would be measured and used to populate the slag model, along with a large database of historical slag temperature, chemistry, and viscosity information they possess – information that would in-crease the predictive accuracy of the slag management model with slag ash chemistry of interest to them. As part of the model evaluation and refinement, the predictive model and its database will be installed on computers at Eastman Chemical and used to manage slag properties in a gasifier during real-time operation of their gasifier. In summary, success of the slag management program at Eastman Chemical would support its use by other gasifier operators to manage slag chemistry and properties over other desired gasification temperature ranges, allowing greater process control and efficiency of their carbon feedstock conversion during gasification.

01 COAL, LIGNITE, AND PEAT↗

KENO V.a Primer: Performing Calculations using SCALE’s Criticality Safety Analysis Sequence (CSAS5) with Fulcrum

The SCALE code system developed at Oak Ridge National Laboratory is widely used and accepted around the world for criticality safety analyses. The well-known KENO V.a three-dimensional Monte Carlo criticality computer code is one of the primary criticality safety analysis tools in SCALE. The KENO V.a primer is designed to help a new user understand and use the SCALE/KENO V.a Monte Carlo code for nuclear criticality safety analyses. It assumes that the user has a college education in a technical field. There is no assumption of familiarity with Monte Carlo codes in general or with SCALE/KENO V.a in particular. The primer is designed to teach by example, with each example illustrating two or three features of SCALE/KENO V.a that are useful in criticality analyses. The primer is based on SCALE 6.2 and 6.3, which includes the Fulcrum graphical user interface (GUI). Each example uses Fulcrum to provide the framework for preparing input data and viewing output results. Starting with a Quickstart section, the primer gives an overview of the basic requirements for SCALE/KENO V.a input and allows the user to quickly run a simple criticality problem with SCALE/KENO V.a. The sections that follow Quickstart include a list of basic objectives at the beginning that identifies the goal of the section and the individual SCALE/KENO V.a features that are covered in detail in the sample problems in that section. Upon completion of the primer, a new user should be comfortable using Fulcrum to set up criticality problems in SCALE/KENO V.a. The primer provides a starting point for the criticality safety analyst who uses SCALE/KENO V.a. Complete descriptions are provided in the SCALE/KENO V.a manual. Although the primer is self-contained, it is intended as a companion volume to the SCALE/KENO V.a training and documentation. The SCALE manual and training schedule are available at https://scale.ornl.gov. The primer provides specific examples of using SCALE/KENO V.a for criticality analyses; the SCALE/KENO V.a manual provides information on the use of SCALE/KENO V.a and all its modules. The primer also contains an appendix with sample input files. In addition, this primer, its errata and sample inputs are also available at https://code.ornl.gov/scale/primers/kenova.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

KENO-VI Primer: Performing Calculations using SCALE’s Criticality Safety Analysis Sequence (CSAS6) with Fulcrum

The SCALE code system developed at Oak Ridge National Laboratory is widely used and accepted around the world for criticality safety analysis. The well-known KENO-VI three-dimensional Monte Carlo criticality computer code is one of the primary criticality safety analysis tools in SCALE. The KENO-VI primer is designed to help a new user understand and use the SCALE/KENO-VI Monte Carlo code for nuclear criticality safety analysis. It assumes that the user has a college education in a technical field. There is no assumption of familiarity with Monte Carlo codes in general or with SCALE/KENO-VI in particular. The primer is designed to teach by example, with each example illustrating two or three features of SCALE/KENO-VI that are useful in criticality analysis. The primer is based on SCALE 6.2 and 6.3, which includes the Fulcrum graphical user interface. Each example uses Fulcrum to provide the framework for preparing input data and viewing output results. Starting with a Quickstart section, the primer gives an overview of the basic requirements for SCALE/KENO-VI input and allows the user to quickly run a simple criticality problem with SCALE/KENO-VI. Each following section begins with a list of basic objectives identifying the goal of the section and the individual SCALE/KENO-VI features covered in detail in the section’s sample problems. Upon completion of the primer, a new user should be comfortable using Fulcrum to set up criticality problems in SCALE/KENO-VI. The primer provides a starting point for the criticality safety analyst who uses SCALE/KENO-VI. Complete descriptions are provided in the SCALE/KENO-VI manual. Although the primer is self-contained, it is intended as a companion volume to the SCALE/KENO-VI training and documentation. The SCALE manual and training schedule are available at https://scale.ornl.gov. The primer provides specific examples of using SCALE/KENO-VI for criticality analysis; the SCALE/KENO-VI manual provides information on the use of SCALE/KENO-VI and all its modules. The primer also contains an appendix with sample input files. In addition, this primer, its errata, and sample inputs are also available at https://code.ornl.gov/scale/primers/kenovi/.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fork Experiments in the Hot Cell Using Spent Fuel Rods for International Nuclear Safeguards

This work leveraged the rare availability of 25 full-length pressurized water reactor spent fuel rods and 1 irradiated mixed-oxide rod at an Oak Ridge National Laboratory hot cell. This was done to collect measurement data with two Fork detectors to assess the detectors’ capability of verifying operator declaration data and detecting partial defects in spent fuel, which are the two primary goals of international safeguards on spent nuclear fuel. The data can also be used to benchmark the ORIGEN module, which has been adopted in the International Atomic Energy Agency’s (IAEA’s)/European Atomic Energy Community’s (Euratom’s) Integrated Review and Analysis Program to predict the Fork detector count rates in real time. In this project, the authors first calibrated two Fork detectors—a standard one and a modified one—by using known strong neutron and gamma sources. Then, the authors measured all 26 fuel rods at multiple locations along the length. The fuel rods were then assembled into three arrays—2 × 2, 3 × 3, and 5 × 5—by using specially designed support grids to mimic fuel assemblies and measure the arrays with both detectors. For the 5 × 5 array, 4 and 8 fuel rods of the array were replaced in two separate cases with short stainless-steel rods to mimic two partial defect scenarios, and the arrays were measured before and after the replacements. Polyethylene blocks were used in this experiment to mimic water. The results show that the Fork detectors were able to verify operator declarations and detect partial defects in spent fuel, and the authors were the first to demonstrate this through experiments. A discovery was also made that determined the root cause of the nonlinear response to gamma dose in the ion chambers used in IAEA and Euratom’s Fork detectors. After the experiments, both detectors were retrieved from the hot cell for future use. The data collected in this project will be used in a parallel International Nuclear Safeguards Engagement Program (INSEP) project to enhance the safeguards in the Finnish spent fuel encapsulation plant, and the data will be useful to other projects in the future given the increased safeguards needs due to spent fuel transfer and disposal activities worldwide.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Verification of Data-Driven Models of Physical Phenomena using Interpretable Approximation

Machine-learned models, specifically neural networks, are increasingly used as “closures” or “constitutive models” in engineering simulators to represent fine-scale physical phenomena that are too computationally expensive to resolve explicitly. However, these neural net models of unresolved physical phenomena tend to fail unpredictably and are therefore not used in mission-critical simulations. In this report, we describe new methods to authenticate them, i.e., to determine the (physical) information content of their training datasets, qualify the scenarios where they may be used and to verify that the neural net, as trained, adhere to physics theory. We demonstrate these methods with neural net closure of turbulent phenomena used in Reynolds Averaged Navier-Stokes equations. We show the types of turbulent physics extant in our training datasets, and, using a test flow of an impinging jet, identify the exact locations where the neural network would be extrapolating i.e., where it would be used outside the feature-space where it was trained. Using Generalized Linear Mixed Models, we also generate explanations of the neural net (à la Local Interpretable Model agnostic Explanations) at prototypes placed in the training data and compare them with approximate analytical models from turbulence theory. Finally, we verify our findings by reproducing them using two different methods.

42 ENGINEERING↗

Calibration of the Diffusivity Predictions of Centipede Using Approximate Bayesian Computation and Applications in Nyx (Engineering Scale) and Xolotl-MARMOT (Meso-Scale) Simulations

Fission gas evolution and release in UO 2 nuclear fuel are important fuel performance metrics and occur in several distinct stages: 1) nucleation, growth and resolution of intra-granular bubbles, 2) diffusion to grain boundaries and 3) nucleation and growth of bubbles at grain boundaries, which eventually form a connected network (percolation) enabling release of gas from grain boundaries through connections to triple junctions, grain edges or free surfaces. The NE-SciDAC project is developing several computational tools to model this problem, which are connected in a hierarchical multi-scale framework. The information transfer in the multi-scale framework is a critical step that, in addition to best-estimates, should include uncertainty quantification. Despite taking a first-principles multi-scale approach, there is a need to perform parameter calibration to ensure consistency with available experimental data. In the present study, uncertainty quantification (UQ) and parameter calibration is demonstrated for one of the lower length scale codes in the multi-scale framework (Centipede) and then the results, including instances of the propagated uncertainties, are used in other codes within the framework, specifically Nyx and Xolotl-MARMOT. We calibrated the model parameters in Centipede, a computer code used to predict diffusivities of uranium (U) and xenon (Xe) in the context of the simulation of fission gas in uranium oxide (UO 2 ) nuclear fuel. The Centipede code depends on 183 parameters, all of which are subject to uncertainty. The three data sets used in our calibration effort are taken from the literature. This data is available as a set of measurements, including measurement errors. Our goal is to calibrate a statistical model that predicts both the value of the measurement and the uncertainty associated with the measurement. We perform a Bayesian calibration of the model parameters using a dedicated approximate Bayesian computation (ABC) likelihood function. To avoid excessive computational costs, we replace the expensive Centipede simulation code by a higher-order surrogate model, constructed using only the 9 most important parameters. These important parameters are identified by a preliminary global sensitivity analysis (GSA) study. Among the important parameters are T0 (the temperature at which UO 2 is perfectly stoichiometric) and Hf_pO2 (the temperature dependence of the oxygen (O) partial pressure) that should be considered as operating conditions to be estimated along with the other parameters. We consider two different cases: one where we define one set of these operating conditions for all data sets, and one where we define distinct operating condition parameters for each data set. The Xe diffusivities predicted by the latter case show distinct features that could not be observed in the former. Next, we use the diffusivity predictions by Centipede as input to Nyx, a reduced order fuel performance code focused on gas behavior alone, in order to estimate quantities associated with inter-granular bubble formation at conditions specified by the experiments. Finally, the diffusivities obtained from the calibrated Centipede runs were used in coupled Xolotl-MARMOT simulations of intra- and inter-granular gas evolution. The results are compared to simulations using the baseline diffusivities from Turnbull et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Basic Research Needs in The Science of Scientific Software Development and Use: Investment in Software is Investment in Science

Increasingly powerful and affordable computing has revolutionized scientific and scholarly discovery across a broad range of fields. Computing relies on software, which has been rapidly growing in scope, diversity, and complexity. At the same time, the methods, processes, and tools used to produce and utilize this essential software are often ad hoc, and the study and improvement of them are often done without the benefit of direct funding or prioritization. Consequently, concerns are growing about the productivity of the developers and users of scientific software, its sustainability, and the trustworthiness of the results that it produces. Increased investment, especially in the characterization and improvement of how scientific software is developed and used, is important for sustaining and improving the impact of software as the scope and complexity of scientific efforts expand. Without this investment, we face the risk of diminishing returns on our software investments because the demands for increased functionality, usability, reliability, and more will not be sufficiently met. The US Department of Energy Office of Science (DOE/SC) is at the forefront of modern software-enabled scientific discovery across numerous areas of computational, experimental, and observational science, including major investments in national user facilities that support these activities. For many years, DOE/SC software investments have provided tremendous value to the scientific community. We want to continue and further improve the value of DOE/SC software efforts by using a scientific approach to understanding and improving how scientific software is developed and used. In December 2021, the DOE/SC Office of Advanced Scientific Computing Research (ASCR) convened a workshop on basic research needs for the Science of Scientific-Software Development and Use (SSSDU). Through keynote presentations, lightning talks, and breakout groups, which built on insights from 124 pre-workshop position papers, participants discussed the current practice of software development, maintenance, evolution, and use, and considered how the scientific method could be used to examine these practices and develop more evidence-based approaches to enhance the impact of software and computing on all areas of science. Workshop participants identified three priority research directions (PRDs) and three important crosscutting themes that center on the following overarching insight: Software has become an essential part of modern science, impacting discoveries, policy, and technological development. To maintain and improve confidence in science delivered via software, we must improve the processes and tools that help us create and use software, and this enhancement requires a deep understanding of the diverse array of teams and individuals doing the work.

97 MATHEMATICS AND COMPUTING↗