Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “information criterion”

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 253 records · Page 14

A Fast Implementation of the ISOCLUS Algorithm

Unsupervised clustering is a fundamental tool in numerous image processing and remote sensing applications. For example, unsupervised clustering is often used to obtain vegetation maps of an area of interest. This approach is useful when reliable training data are either scarce or expensive, and when relatively little a priori information about the data is available. Unsupervised clustering methods play a significant role in the pursuit of unsupervised classification. One of the most popular and widely used clustering schemes for remote sensing applications is the ISOCLUS algorithm, which is based on the ISODATA method. The algorithm is given a set of n data points (or samples) in d-dimensional space, an integer k indicating the initial number of clusters, and a number of additional parameters. The general goal is to compute a set of cluster centers in d-space. Although there is no specific optimization criterion, the algorithm is similar in spirit to the well known k-means clustering method in which the objective is to minimize the average squared distance of each point to its nearest center, called the average distortion. One significant feature of ISOCLUS over k-means is that clusters may be merged or split, and so the final number of clusters may be different from the number k supplied as part of the input. This algorithm will be described in later in this paper. The ISOCLUS algorithm can run very slowly, particularly on large data sets. Given its wide use in remote sensing, its efficient computation is an important goal. We have developed a fast implementation of the ISOCLUS algorithm. Our improvement is based on a recent acceleration to the k-means algorithm, the filtering algorithm, by Kanungo et al.. They showed that, by storing the data in a kd-tree, it was possible to significantly reduce the running time of k-means. We have adapted this method for the ISOCLUS algorithm. For technical reasons, which are explained later, it is necessary to make a minor modification to the ISOCLUS specification. We provide empirical evidence, on both synthetic and Landsat image data sets, that our algorithm's performance is essentially the same as that of ISOCLUS, but with significantly lower running times. We show that our algorithm runs from 3 to 30 times faster than a straightforward implementation of ISOCLUS. Our adaptation of the filtering algorithm involves the efficient computation of a number of cluster statistics that are needed for ISOCLUS, but not for k-means.

Memarsadeghi, Nargess↗

An Examination of Two Non-Cooperative Detect and Avoid Well Clear Definitions

NASA’s Unmanned Aircraft Systems Integration into the National Airspace System (UAS in the NAS) project examines the technical barriers associated with the operation of UAS in civil airspace. The present study explored the differential effects of two candidate non-cooperative Detect-and-Avoid Well Clear (DWC)definitions on pilot and system performance in a human-in-the-loop simulation. Active-duty UAS pilots were recruited to maintain DWC against representative Class 4 encounter types with a low size, weight, and power (SWaP) radar declaration range of 3.5 nautical miles (nmi). Objective performance indicated that pilots could consistently maintain DWC against non-cooperative intruders with either DWC candidate, with negligible differences in response times and separation performance against caution and warning-level threats. While losses of DWC were avoided at rates comparable to Phase 1 findings, pilots uploaded their responses to caution-level alerts over 5 seconds faster in the current setup relative to Phase 1. Encounters with faster closure rates were susceptible to shortened caution-level alert durations, especially when employing the DWC criterion with the additional ‘Tau’ (temporal) component. Consequently, caution-level threats frequently elevated to warning-level status (nearly twice as often with theTau candidate). The variable caution alert durations appeared to impact pilots’ coordination with air traffic control (ATC), as ATC approval rates were lower with the ‘Tau’and ‘Disc’ candidates relative to Phase 1 research. Ultimately, the increased alerting time enabled by the Disc candidate deemed it more suitable for any reductions to the assumed radar declaration range requirement, which was re-evaluated in a follow-on study. Findings from this study will inform Phase 2 Minimum Operational Performance Standards (MOPS)development for UAS with alternative surveillance equipment and performance capabilities.

Kevin J Monk↗

Sensing, Spectra and Scaling: What's in Store for Land Observations

Bill Pecora's 1960's vision of the future, using spacecraft-based sensors for mapping the environment and exploring for resources, is being implemented today. New technology has produced better sensors in space such as the Landsat Thematic Mapper (TM) and SPOT, and creative researchers are continuing to find new applications. However, with existing sensors, and those intended for launch in this century, the potential for extracting information from the land surface is far from being exploited. The most recent technology development is imaging spectrometry, the acquisition of images in hundreds of contiguous spectral bands, such that for any pixel a complete reflectance spectrum can be acquired. Experience with Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) has shown that, with proper attention paid to absolute calibration, it is possible to acquire apparent surface reflectance to 5% accuracy without any ground-based measurement. The data reduction incorporates in educated guess of the aerosol scattering, development of a precipitable water vapor map from the data and mapping of cirrus clouds in the 1.38 micrometer band. This is not possible with TM. The pixel size in images of the earth plays and important role in the type and quality of information that can be derived. Less understood is the coupling between spatial and spectral resolution in a sensor. Recent work has shown that in processing the data to derive the relative abundance of materials in a pixel, also known is unmixing, the pixel size is an important parameter. A variance in the relative abundance of materials among the pixels is necessary to be able to derive the endmembers or pure material constituent spectra. In most cases, the 1 km pixel size for the Earth Observing System Moderate Resolution Imaging Spectroradiometer (MODIS) instrument is too large to meet the variance criterion. A pointable high spatial and spectral resolution imaging spectrometer in orbit will be necessary to make the major next step in our understanding of the solid earth surface and its changing face.

Goetz, Alexander F. H.↗

Decode the Workload: Training Deep Learning Models for Efficient Compute Cluster Representation

Monitoring the status of a high throughput computing cluster running computationally intensive production jobs is a crucial yet challenging system administration task due to the complexity of such systems. To this end, we train autoencoders using the Linux kernel CPU metrics of the cluster. Additionally, we explore assisting these models with graph neural networks to share information across threads within a compute node. The models are compared in terms of their ability to: 1) Produce a compressed latent representation that captures the salient features of the input, 2) Detect anomalous activity, and 3) Make distinction between different kinds of jobs run at Jefferson Lab. The goal is to have a robust encoder whose compressed embeddings are used for several downstream tasks. We extend this study further by deploying these models in a human-in-the-loop production-based setting for the anomaly detection task and discuss the associated implementation aspects such as continual learning and the criterion to generate alarms. This study represents a first step in the endeavor towards building self-supervised large-scale foundation models for computing centers.

Mohammed, Ahmed↗

Implications of the weak gravity conjecture for tidal Love numbers of black holes

The weak gravity conjecture indicates that extremal black holes in the low-energy effective field theory should be able to decay. This criterion gives rise to nontrivial constraints on the coefficients of higher-order derivative corrections to gravity. In this paper, we investigate the tidal deformability of neutral black holes due to higher-order derivative corrections. As a proof of concept, we consider a correction of cubic order in the Riemann curvature tensor. The tidal Love numbers of neutral black holes receive leading-order corrections from higher-order derivative terms, since black holes in pure General Relativity have vanishing tidal Love number. We conclude that the interplay between the tidal deformability of black holes and the weak gravity conjecture provides useful information about the effective field theory.

79 ASTRONOMY AND ASTROPHYSICS↗

Blind Modeling Validation Exercises Using the Horizontal Dry Cask Simulator

The U.S. Department of Energy (DOE) established a need to understand the thermal-hydraulic properties of dry storage systems for commercial spent nuclear fuel (SNF) in response to a shift towards the storage of high-burnup (HBU) fuel (> 45 gigawatt days per metric ton of uranium, or GWd/MTU). This shift raises concerns regarding cladding integrity, which faces increased risk at the higher temperatures within spent fuel assemblies present within HBU fuel compared to low-burnup fuel (≤ 45 GWd/MTU). A dry cask simulator (DCS) was built at Sandia National Laboratories (SNL) in Albuquerque, New Mexico to produce validation-quality data that can be used to test the accuracy of the modeling used to predict cladding temperatures. These temperatures are critical to evaluating cladding integrity throughout the storage cycle of commercial spent nuclear fuel. A model validation exercise was previously carried out for the DCS in a vertical configuration. Lessons learned during the previous validation exercise have been applied to a new, blind study using a horizontal dry cask simulator (HDCS). Three modeling institutions – the Nuclear Regulatory Commission (NRC), Pacific Northwest National Laboratory (PNNL), and Empresa Nacional del Uranio, S.A., S.M.E. (ENUSA) – were granted access to the input parameters from the DCS Handbook, SAND2017-13058R, and results from a limited data set from the horizontal BWR dry cask simulator tests reported in the HDCS update report, SAND2019-11688R. With this information, each institution was tasked to calculate peak cladding temperatures and air mass flow rates for ten HDCS test cases. Axial as well as vertical and horizontal transverse temperature profiles were also calculated. These calculations were done using modeling codes (ANSYS/Fluent, STAR-CCM+, or COBRA-SFS), each with their own unique combination of modeling assumptions and boundary conditions. For this validation study, the ten test cases of the horizontal dry cask simulator were defined by three independent variables – fuel assembly decay heat (0.5 kW, 1 kW, 2.5 W, and 5 kW), internal backfill pressure (100 kPa and 800 kPa), and backfill gas (helium and air). The plots provided in Chapter 3 of this report show the axial, vertical, and horizontal temperature profiles obtained from the dry cask simulator experiments in the horizontal configuration and the corresponding models used to describe the thermal-hydraulic behavior of this system. The tables provided in Chapter 3 illustrate the closeness of fit of the model data to the experiment data through root mean square (RMS) calculations of the error in peak cladding temperatures (PCTs), PCT axial locations, axial temperature profiles, vertical and horizontal temperature profiles at two different axial locations, and air mass flow rates for the ten test cases, normalized by the experimental results. The model results are assigned arbitrary model numbers to retain anonymity. Due to the relatively flat axial temperature profiles, small temperature gradients resulted in large deviations of all models’ PCT axial location from the experimental PCT axial location. When the PCT axial location error is excluded in the calculation of the combined RMS of the normalized errors that considers PCT, the temperature profiles, and the air mass flow rates, the model data fits the experimental data to within 5%. When the vault information is excluded, the model data fits the experimental data to within 2.5%. An error analysis was developed further for one model, using the model and experimental uncertainties in each validation parameter to calculate validation uncertainties. The uncertainties for each parameter were used to define quantifiable validation criteria. For this analysis, the model was considered validated for a given comparison metric if the normalized error in that metric divided by the validation uncertainty was less than or equal to 1. When considering the combined RMS of the normalized errors of all metrics divided by their validation uncertainties, the model was found to have satisfied the criterion for model validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Wavefront Control Toolbox for James Webb Space Telescope Testbed

We have developed a Matlab toolbox for wavefront control of optical systems. We have applied this toolbox to the optical models of James Webb Space Telescope (JWST) in general and to the JWST Testbed Telescope (TBT) in particular, implementing both unconstrained and constrained wavefront optimization to correct for possible misalignments present on the segmented primary mirror or the monolithic secondary mirror. The optical models implemented in Zemax optical design program and information is exchanged between Matlab and Zemax via the Dynamic Data Exchange (DDE) interface. The model configuration is managed using the XML protocol. The optimization algorithm uses influence functions for each adjustable degree of freedom of the optical mode. The iterative and non-iterative algorithms have been developed to converge to a local minimum of the root-mean-square (rms) of wavefront error using singular value decomposition technique of the control matrix of influence functions. The toolkit is highly modular and allows the user to choose control strategies for the degrees of freedom to be adjusted on a given iteration and wavefront convergence criterion. As the influence functions are nonlinear over the control parameter space, the toolkit also allows for trade-offs between frequency of updating the local influence functions and execution speed. The functionality of the toolbox and the validity of the underlying algorithms have been verified through extensive simulations.

Shiri, Ron↗

HIRS-AMTS satellite sounding system test - Theoretical and empirical vertical resolving power

The present investigation is concerned with the vertical resolving power of satellite-borne temperature sounding instruments. Information is presented on the capabilities of the High Resolution Infrared Radiation Sounder (HIRS) and a proposed sounding instrument called the Advanced Moisture and Temperature Sounder (AMTS). Two quite different methods for assessing the vertical resolving power of satellite sounders are discussed. The first is the theoretical method of Conrath (1972) which was patterned after the work of Backus and Gilbert (1968) The Backus-Gilbert-Conrath (BGC) approach includes a formalism for deriving a retrieval algorithm for optimizing the vertical resolving power. However, a retrieval algorithm constructed in the BGC optimal fashion is not necessarily optimal as far as actual temperature retrievals are concerned. Thus, an independent criterion for vertical resolving power is discussed. The criterion is based on actual retrievals of signal structure in the temperature field.

Thompson, O. E.↗

The Evolving Chronology of Moon Formation

Defining the age of the Moon has proven to be an elusive task because it requires reliably dating lunar samples using radiometric isotopic systems that record fractionation of parent and daughter elements during events that are petrologically associated with planet formation. Crystallization of the magma ocean is the only event that unambiguously meets this criterion because it probably occurred within tens of millions of years of Moon formation. There are three dateable crystallization products of the magma ocean: mafic mantle cumulates, felsic crustal cumulates, and late-stage crystallization products known as urKREEP (uniform residuum K, rare earth elements, and P). Although ages for these materials in the literature span 200 million years, there is a preponderance of reliable ages around 4.35 billion years recorded in all three lunar rock types. This age is also observed in many secondary crustal rocks, indicating that they were produced contemporaneously (within uncertainty of the ages), possibly during crystallization and overturn of the magma ocean. •The duration of planet formation is key information in understanding the mechanisms by which the terrestrial planets formed. •Ages of the oldest lunar rocks range widely, reflecting either the duration of Moon formation or disturbed ages caused by impact metamorphism. •Ages determined for compositionally distinct crust and mantle materials produced by lunar magma ocean differentiation cluster near 4.35 Gyr. •The repeated occurrence of 4.35 Gyr ages implies that Moon formation occurred late in Solar System history, likely by giant impact into Earth.

79 ASTRONOMY AND ASTROPHYSICS↗

Winds in hot, subluminous stars

Ultraviolet spectra over the wavelength range 1150 to 1750 A of the nuclei of two planetary nebulae are used to derive information concerning their winds. The nucleus of NGC 6826 is a 10th magnitude star with an O3f-type spectrum, an apparent terminal wind velocity of 2000 km/s, and a general level of ionization similar to that in Zeta Pup. The nucleus of Abell 78 is a 13th magnitude star having an OVI-type spectrum, an effective temperature between 100,000 and 150,000 K, and a gravity greater than 100,000 cm/second squared. Its apparent terminal wind velocity is 3400 km/s, and it has a very high general level of ionization. It is suggested that a luminosity-to-mass ratio (in solar units) greater than about 1600 is now an applicable mass-loss criterion for all hot stars, including planetary nuclei.

Heap, S. R.↗

Image Distinguishability Analysis Testing Through Principal Components and Its Application to Hot Spot Scale Invariance

Hot spots are spatial regions of intense energy localization that govern initiation of secondary high explosives. Studies that characterize or compare simulated hot spots are frequently either qualitatively descriptive or resort to quantitative distribution functions that neglect stochastic variations and spatial correlations—effects that are also neglected in common comparison tests like the Kolmogorov–Smirnov test. To this end, we develop an image distinguishability analysis (IDA) test based on principal component (PC) analysis that makes pixel-by-pixel comparisons between small, for example, O(<10), image data sets. The IDA test makes comparisons through a generalized distance metric in the PC space and a test statistic that is derived to calculate mathematical equation-values. Here, we derive a statistical distribution and criticality criterion to determine whether images are distinguishable from established baselines. We apply the IDA test on images generated from molecular dynamics simulations of hot spots from pore collapse in TATB to assess scale invariance in the complex patterns of hot spots that form in a representative high explosive crystal. The IDA test shows that TATB hot spot spatial temperature fields and their derived temperature histograms exhibit scale-invariant features over specific intervals of shock orientation, strength, and initial pore diameter. However, the IDA test also shows that qualitatively different conclusions regarding invariance can be reached depending on whether the hot spot is treated as a spatially correlated field as opposed to a distribution function that lacks spatial information.

organic↗

Molten Chloride Thermophysical Properties, Chemical Optimization, and Purification

Next-generation concentrating solar power (CSP) technology requires high-temperature heat-transfer fluids and thermal energy storage that can operate in the temperature range of 500°–720°C. This high temperature demands a new chemistry for the heat-transfer fluids that must have high-temperature thermal and chemical stabilities, high energy density, and low corrosion rate on metallic materials used in the heat exchangers, piping, and thermal energy storage tanks. A ternary MgCl 2 -KCl-NaCl chloride-salt system has been proposed to achieve such goals, given its high-temperature thermal stability. The chloride-salt-based systems have been used by the metallurgy industry as the electrolyte to electrolytically extract metallic magnesium. They have also attracted attentions from the nuclear industry for developing the molten salt fast reactors that can operate at close to atmospheric pressure and at higher temperature than traditional light-water reactors. However, the corrosion properties of the ternary salt are not well known. The corrosion mechanism is believed to be correlated to the presence of MgOHCl impurity at the operation temperature. MgOHCl is fundamentally caused by the hygroscopic nature of the MgCl 2 component. Therefore, the first key objective of this project is to understand the nature of the corrosiveness of the ternary chloride-salt system and its prevention by designing new lab-scale purification processes. Because of its importance, a working group called the Chloride Salt Collective across multiple U.S. Department of Energy national laboratories is established to gain experimental agreement on the corrosion behaviors. In addition, water generation during salt dehydration could lead to highly corrosive HCl/water mixture via a hydrolysis process between water and anhydrous MgCl 2 or MgCl 2 hydrates. MgOHCl will eventually be converted to MgO particulates regardless of the chemical path chosen (i.e., by thermal decomposition, chemical purification with an active metal, or electrochemical purification utilizing the reduction power of electron flows). Large MgO particles with high hardness can wear out bearing and bushing and even clog the delicate channels in the primary heat exchanger. Because of this, a second objective is to understand the chemistry behind salt purification and transfer the scientific understanding to a set of engineering know-how for large-scale salt melter design. Because there is limited information on the thermophysical properties and the sensitivity of these properties around the targeted salt composition, a third objective is to measure thermophysical properties relevant to system and component design and support Topic 1 Liquid Pathway development. For example, composition shift due to salt reaction with water can be a potential issue because it can lead to changes of thermophysical properties, such as (1) melting point (or freeze point), which can jeopardize safe operation of the supercritical CO 2 power cycle proposed for next-generation CSP, and (2) heat capacity and density, which can decrease the overall power-generation capability of the plants. The value of the thermal conductivity of the selected salt composition is also a critical criterion for the solar receiver design.

14 SOLAR ENERGY↗

A Model for the Onset of Vortex Breakdown

A large body of information exists on the breakdown of incompressible streamwise vortices. Less is known about vortex breakdown at high speeds. An interesting example of supersonic vortex breakdown is the breakdown induced by the interaction of vortices with shock waves. The flow in supersonic engine inlets and over high-speed delta wings constitute technologically important examples of this phenomenon, which is termed 'shock-induced vortex breakdown'. In this report, we propose a model to predict the onset of shock-induced vortex breakdown. The proposed model has no adjustable constants, and is compared to both experiment and computation. The model is then extended to consider two other problems: the breakdown of a free compressible vortex, and free incompressible vortex breakdown. The same breakdown criterion is used in all three problems to predict the onset of breakdown. Finally, a new breakdown map is proposed that allows the simultaneous comparison of data from flows ranging from incompressible breakdown to breakdown induced by a shock wave.

Mahesh, K.↗

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗

Investigation of Conjugate Heat Transfer in Turbine Blades and Vanes

We report on work carried out to develop a 3-D coupled Finite Volume/BEM-based temperature forward/flux back (TFFB) coupling algorithm to solve the conjugate heat transfer (CHT) which arises naturally in analysis of systems exposed to a convective environment. Here, heat conduction within a structure is coupled to heat transfer to the external fluid which is convecting heat into or out of the solid structure. There are two basic approaches to solving coupled fluid structural systems. The first is a direct coupling where the solution of the different fields is solved simultaneously in one large set of equations. The second approach is a loose coupling strategy where each set of field equations is solved to provide boundary conditions for the other. The equations are solved in turn until an iterated convergence criterion is met at the fluid-solid interface. The loose coupling strategy is particularly attractive when coupling auxiliary field equations to computational fluid dynamics codes. We adopt the latter method in which the BEM is used to solve heat conduction inside a structure which is exposed to a convective field which in turn is resolved by solving the NASA Glenn compressible Navier-Stokes finite volume code Glenn-HT. The BEM code features constant and bi-linear discontinuous elements and an ILU-preconditioned GMRES iterative solver for the resulting non-symmetric algebraic set arising in the conduction solution. Interface of flux and temperature is enforced at the solid/fluid interface, and a radial-basis function scheme is used to interpolated information between the CFD and BEM surface grids. Additionally, relaxation is implemented in passing the fluxes from the conduction solution to the fluid solution. Results from a simple test example are reported.

Kassab, A. J.↗

Estimation of proportions in mixed pixels through their region characterization

A region of mixed pixels can be characterized through the probability density function of proportions of classes in the pixels. Using information from the spectral vectors of a given set of pixels from the mixed pixel region, expressions are developed for obtaining the maximum likelihood estimates of the parameters of probability density functions of proportions. The proportions of classes in the mixed pixels can then be estimated. If the mixed pixels contain objects of two classes, the computation can be reduced by transforming the spectral vectors using a transformation matrix that simultaneously diagonalizes the covariance matrices of the two classes. If the proportions of the classes of a set of mixed pixels from the region are given, then expressions are developed for obtaining the estmates of the parameters of the probability density function of the proportions of mixed pixels. Development of these expressions is based on the criterion of the minimum sum of squares of errors. Experimental results from the processing of remotely sensed agricultural multispectral imagery data are presented.

Chittineni, C. B.↗

In-trail dynamics of multiple CDTI-equipped aircraft queues

One of the potential problems of in-trail self-spacing with a Cockpit Display of Traffic Information (CDTI) is whether dynamic oscillations would occur in a queue of aircraft flying an approach, similar to the ""accordion'' effect seem with the queue of automobiles in stop-and-go traffic. In order to gain some insight into this potential problem, a brief experiment was conducted with the Transport Systems Research Vehicle (TSRV) ground-based simulator equipped with CDTI which presented the position of other aircraft in the area. Three simulation sessions were conducted wherein queues of up to nine aircraft were built, each one self-spacing on the preceding aircraft. The aircraft crews were rotated to ensure that the pilots had no prior knowledge of the lead aircraft behavior they would be following. Two different spacing criteria were employed: a constant time predictor criterion and a constant time delay criterion. The experiment failed to uncover any dynamic oscillatory tendencies in queues of seven to nine aircraft.

Kelly, J. R.↗

Evaluation of Mechanistic and Empirical Models against Existing FFRD and LOCA Experimental Databases

The desire of the nuclear industry to improve the economics of existing nuclear power plants has necessitated research into the potential of a phenomenon known as fuel fragmentation, relocation, and dispersal (FFRD). This phenomenon is possible during a loss-of-coolant accident (LOCA) at relatively high burnup. The Nuclear Energy Advanced Modeling and Simulation program has been developing simulation capabilities for FFRD and LOCA in the BISON fuel performance code for multiple years. This year, the effort has been focused on evaluating new lower length scale informed pulverization thresholds as well as updating and adding new empirical models for various phenomena. Models added or updated this year include a preliminary transient fission gas release model, new high-temperature Zircaloy creep models, a Zircaloy rupture opening area model, and a temperature-dependent emissivity during radiation from the fuel rod to the surrounding atmosphere. The models are verified through implementation tests to demonstrate code correctness after addition to BISON. The models are then assessed against a subset of the existing BISON validation suite for LOCA and FFRD cases. Two new cases, Studsvik Rods 192 and 193, have been added. The results indicate that the inclusion of a bubble pressure evolution model in the bubbles in the high-burnup structure has the largest impact compared to the 3D fracture criterion on reducing the calculated amount of pulverized fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗