Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Adaptive Sampling”

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 1,477 records · Page 82

Explainable AI (XAI)-driven vibration sensing scheme for surface quality monitoring in a smart surface grinding process

Local Interpretable and Model-agnostic Explanation (LIME), an explainable artificial intelligence (XAI) approach is adapted to identify the globally important time-frequency bands for predicting average surface roughness (Ra) in a smart grinding process. The smart grinding setup consisted of a Supertech CNC precision surface grinding machine, instrumented with a Dytran piezoelectric accelerometer attached to the tailstock quill along the tangential direction (Y-axis). For every grinding pass, vibration signatures were captured, and the ground truth surface roughness values were recorded using a Mahr Marsurf M300C portable surface roughness profilometer. The roughness values ranged from 0.06 to 0.14 microns over the complete set of experiments. Time-frequency domain spectrogram frames were extracted for each of the vibration signals collected during the grinding process. Convolutional Neural Networks (CNNs) were modeled to predict the surface roughness based on these spectrogram frames and their image augmentations. The best CNN model was able to predict the roughness values with an overall R2-score of 0.95, training R2-score of 0.99, and testing R2-score of 0.81 with only 80 sets of vibration signals corresponding to 4 experiments with 20 trials each. Although the data size is not large enough to guarantee such performance metrics in real-world scenarios, one can extract statistically consistent explanations underlying the relationships these complex deep learning models capture. Further, the LIME methodology was implemented on the developed surface roughness CNN model to identify the important time-frequency bands (i.e., the superpixels of a spectrogram) influencing the predictions. Based on the identified important regions on the spectrogram frames, the corresponding frequency characteristics were determined that influence the surface roughness predictions. The important frequency range based on LIME results was approximately 11.7 to 19.1 kHz. The power of XAI was demonstrated by cutting down the sampling rate from 160 kHz to 30, 20, 10, and 5 kHz based on the important frequency range and considering Nyquist criteria. Separate CNN models were developed for these ranges by only extracting time-frequency contents below their corresponding Nyquist cut-offs. A proper data acquisition strategy is proposed by comparing the model performances to argue the selection of a sufficient sampling rate to capture the grinding process successfully and robustly.

42 ENGINEERING↗

Reprint of: Explainable AI (XAI)-driven vibration sensing scheme for surface quality monitoring in a smart surface grinding process

Local Interpretable and Model-agnostic Explanation (LIME), an explainable artificial intelligence (XAI) approach is adapted to identify the globally important time-frequency bands for predicting average surface roughness (Ra) in a smart grinding process. The smart grinding setup consisted of a Supertech CNC precision surface grinding machine, instrumented with a Dytran piezoelectric accelerometer attached to the tailstock quill along the tangential direction (Y-axis). For every grinding pass, vibration signatures were captured, and the ground truth surface roughness values were recorded using a Mahr Marsurf M300C portable surface roughness profilometer. The roughness values ranged from 0.06 to 0.14 microns over the complete set of experiments. Time-frequency domain spectrogram frames were extracted for each of the vibration signals collected during the grinding process. Convolutional Neural Networks (CNNs) were modeled to predict the surface roughness based on these spectrogram frames and their image augmentations. The best CNN model was able to predict the roughness values with an overall R2-score of 0.95, training R2-score of 0.99, and testing R2-score of 0.81 with only 80 sets of vibration signals corresponding to 4 experiments with 20 trials each. Although the data size is not large enough to guarantee such performance metrics in real-world scenarios, one can extract statistically consistent explanations underlying the relationships these complex deep learning models capture. Further, the LIME methodology was implemented on the developed surface roughness CNN model to identify the important time-frequency bands (i.e., the superpixels of a spectrogram) influencing the predictions. Based on the identified important regions on the spectrogram frames, the corresponding frequency characteristics were determined that influence the surface roughness predictions. The important frequency range based on LIME results was approximately 11.7 to 19.1 kHz. The power of XAI was demonstrated by cutting down the sampling rate from 160 kHz to 30, 20, 10, and 5 kHz based on the important frequency range and considering Nyquist criteria. Separate CNN models were developed for these ranges by only extracting time-frequency contents below their corresponding Nyquist cut-offs. A proper data acquisition strategy is proposed by comparing the model performances to argue the selection of a sufficient sampling rate to capture the grinding process successfully and robustly.

47 OTHER INSTRUMENTATION↗

Rapid Bayesian High Entropy Alloy Designs Fabricated via Wire Arc Additive Manufacturing

Purpose: This project seeks to demonstrate a new high-throughput (rapid) alloy design technique applied to creating new high entropy alloys (HEAs) for extreme environments. High entropy alloys shift the design paradigm from being focused on a single principal element (e.g. nickel-based alloys) to target alloys that include high atomic fractions (X >10%) of multiple elements. These HEA materials can exhibit sluggish diffusion and enhanced corrosion resistance, ideal for potential applications in advanced ultra supercritical (A-USC) steam cycles for power generation. Scope: The addition of multiple elements in high atomic fractions creates an enormous design space that cannot easily be investigated by traditional material design strategies such as designed of experiments (DOE). This project utilizes a Bayesian machine learning algorithm that has been modified to work with calculation of phase diagrams (CALPHAD) software. This Bayesian algorithm reduces manual inputs and increase the likelihood of achieving an optimal solution. Compositional inputs to this algorithm will be assessed using existing material property models for high temperature strength and corrosion resistance. The target for alloy performance will be a 15% (~100 ⁰C) increase in allowable service temperature beyond heat-resistant stainless steels while maintaining or improving alloy cost and corrosion resistance. Haynes 230 was selected as a baseline, which is 57 wt% Ni with 22 wt% Cr 14 wt% W, and 2 wt% Mo as solid solution strengtheners. In addition to rapid design via Bayesian machine learning, the alloys were rapidly fabricated using a multi-wire arc additive manufacturing (mWAAM) technique which allows for precise control of alloy composition and assessing of alloy design “windows” to study composition effects. Build speeds for wire-arc additive processes are among the highest for additive technologies enabling rapid and reliable sample fabrication when compared to conventional methods such as arc button melting. The mWAAM samples will be rapidly characterized via instrumented indentation for room temperature modulus and strength and for elevated temperature strength via hot hardness tests. After being screened with hardness testing, potential alloys will be further evaluated with conventional microscopy techniques including scanning electron microscopy (SEM) and transmission electron microscopy (TEM) to assess agreement with modeling results. The most promising compositions will also be evaluated by printing full sized tensile specimens for mechanical behavior tests at elevated temperatures. Results: Bayesian machine learning of a single performance function was initially used to optimize five performance metrics: 1) single phase stability, 2) yield strength, 3) creep resistance (low diffusion coefficient), 4) freezing range (weldability), and 5) material cost. The single performance function was suboptimal as assumptions had to be made about the results while formulating the optimization. A goal-oriented Bayesian optimization strategy (Hanaoka, 2021) was implemented with CALPHAD for use with the five metrics above. This multi-objective Bayesian optimization (MOBO) enabled the design of NiCrCoFe alloys with V and W additions. A base composition of NiCoCr was selected as Ni provides a stable FCC matrix, Cr aids corrosion/oxidation resistance, and Co is a solid-solutions strengthener that also improves creep by increasing the activation energy. Fe helps reduce diffusion coefficients and cost. Finally, V and W were selected for their reasonable solubility and high atomic misfit to aid in solid solution strengthening. Cracking of the mWAAM specimens was an early issue, and the Easton solidification cracking model (Easton et al., 2014a) was selected for addition to the MOBO function. High performing alloys fabricated by mWAAM included Ni 28 Cr 25 Co 26 Fe 15 V 8 and Ni 62 Cr 18 Co 1 Fe 3 W 15 . It was observed that even after adapting the mWAAM process for W, the W did not fully dissolve. To fully evaluate the Ni 62 Cr 18 Co 1 Fe 3 W 15 composition, a cored wire (80-20 NiCr sheath/powder core) was manufactured and printed via WAAM, and HIP’ing was utilized to homogenize and densify the printed alloy. The V and W alloys produced met metrics 1 (solid solution), 4 (solidification cracking), and 5 (cost). However, an unmodeled mechanism of thermal stress cracking was identified in the WAAM produced materials, perhaps exacerbated by the lack of grain boundary strengthening elements (B, C). Conclusions & Recommendations: A high-throughput (rapid) alloy design technique was applied to designing and manufacturing new high entropy alloys (HEAs) for extreme environments utilizing MOBO and mWAAM. The developed process was rapid and effective in addressing the mechanisms included in the model. The lack of grain boundary strengthening element additions (e.g., B, C) was a simplification that likely produced thermal stress cracking that turned into a large part of the investigation. Additions on the order of 0.005 wt% B and 0.05 wt% C likely would have minimized thermal stress grain boundary cracking. Overall, the high throughput design strategy is promising for rapid design of metrics-driven alloys for advanced ultra supercritical (A-USC) steam cycles for power generation. The MOBO and mWAAM process could be commercialized to accelerate metrics-driven alloy design. In addition, the cored-wire process utilized for scale-up is a promising high-volume process for WAAM alloy development and scale-up.

36 MATERIALS SCIENCE↗

Machine Learning Algorithms for Matching Theories, Simulations, and Observations in Cosmology (Final Project)

The search for simple and elegant equations to explain and predict observable natural phenomena has been a driving force in scientific theory for several centuries. The great successes of Newton's laws, Maxwell's equations, the Schrodinger equation, and the theory of General Relativity have encouraged this perspective, and perhaps a beautiful equation comprising a “theory of everything” remains to be discovered. But now scientists are grappling with systems of ever-increasing complexity, fraught with nonlinearity, and with the quantities of interest often only indirectly observed. Nature may resist a simple description, and the most important discoveries of the next century may be complex theories with countless variables and parameters. The era of big data opens up a promising new approach to scientific discovery in this setting. The predictions of modern theories, even if complex or nonlinear, can be examined through detailed and computationally-intensive simulations that take days or weeks to run and fill petabytes of disk space. The upcoming generation of scientific instruments will provide petabytes of observations at unprecedented resolution and depth. The challenge is to compare data to simulations so as to test theories and identify those that best match the observations. The data are numerous and rich but are subject to noise and various systematic biases. The simulations are large and detailed but are so costly that the effective sample size is small. Moreover, the mapping from data to simulation to the parameters of interest is typically complicated and ill-conditioned. Although these challenges, and the solutions we propose, apply across the sciences, we will ground our efforts by focusing on cosmology. Here, we are truly entering the era of big data in several senses. Ever larger experiments survey larger and larger cosmological volumes, resulting in enormous amounts of pixel data to process and extract the information of cosmological interest. These include large imaging surveys, which collect all the light in broad wavelength ranges over significant fractions of the sky, and spectroscopic surveys, which measure the light as a function of wavelength to map out 3D distances. These experiments are designed to shed light on the big unanswered questions in cosmology. Distilling the petabytes of data into observed quantities for comparison with theory is a daunting task. But these large experiments do not only result in large amounts of data to analyze: if we wish to compare to theories of physics, then we need increasingly extensive predictions for the observed quantities as well, which requires ever-larger cosmological simulations that take substantial -- and at times prohibitive -- computational resources to produce. Our goal is to develop statistical and machine learning methods for using observed and simulated data to advance machine learning with applications cosmology. In particular, we focus our research efforts on the following tasks: The challenge for automated science is that it is computationally impossible and statistically dangerous to consider every possible model in order to find the best one. We will develop Bayesian Optimization based active-learning methods that accelerate both the execution of the simulations and the search for best-fitting parameters. The key idea is to make the simulations adaptive -- across resolutions, time, and parameters -- using the data to search as the simulation runs; Most ML methods operate on simple finite dimensional feature vectors. However, many cosmology and other science applications require ML methods that can operate on more complex objects as inputs or outputs such as functions, distributions, or sets and point clouds. Our goal is to develop efficient ML methods for this problem and demonstrate their applicability in Cosmology, Astrophysics, and other science problems.

79 ASTRONOMY AND ASTROPHYSICS↗

Synthesis of magnetically recyclable ZIF-8@SiO{sub 2}@Fe{sub 3}O{sub 4} catalysts and their catalytic performance for Knoevenagel reaction

Novel magnetic ZIF-8@SiO{sub 2}@Fe{sub 3}O{sub 4} catalysts were synthesized by encapsulating magnetic SiO{sub 2}@Fe{sub 3}O{sub 4} nanoparticles into ZIF-8 through in situ method. The structures of the catalysts were characterized by TEM, SEM, XRD, FT-IR, VSM, N{sub 2} adsorption/desorption and CO{sub 2}-TPD technology. The catalytic activity and recovery properties of the catalysts for the Knoevenagel reaction of p-chlorobenzaldehyde with malononitrile were evaluated. The results showed that the magnetic ZIF-8@SiO{sub 2}@Fe{sub 3}O{sub 4} catalysts had the larger surface areas, the suitable superparamagnetism, and good catalytic activity for Knoevenagel reaction. The conversion of p-chlorobenzaldehyde can reach ~98% and the selectivity of the production can reach ~99% over35.8%ZIF-8@SiO{sub 2}@Fe{sub 3}O{sub 4} (MZC-5) catalyst under the reaction condition of 25 °C and 4 h. The magnetic ZIF-8@SiO{sub 2}@Fe{sub 3}O{sub 4} catalysts also had good substrates adaptation. After reaction, the catalyst can be easily separated from the reaction mixture by an external magnet. The recovery catalyst can be reused five times and the conversion of p-chlorobenzaldehyde can be kept over 90%. - Graphical abstract: Novel magnetically recyclable ZIF-8@SiO{sub 2}@Fe{sub 3}O{sub 4} catalysts were synthesized by encapsulating magnetic SiO{sub 2}@Fe{sub 3}O{sub 4} nanoparticles into ZIF-8 and the as-synthesized catalysts exhibited a good catalytic activity for the Knoevenagel reaction. - Highlights: • A series of novel magnetic ZIF-8@SiO{sub 2}@Fe{sub 3}O{sub 4} catalysts were synthesized. • The catalysts had the larger surface areas and the suitable superparamagnetism. • The catalysts exhibited good catalytic activity for the Knoevenagel reaction. • After reaction the catalyst can be easily separated by an external magnet. • The recovery catalyst can be reused five times and can keep its catalytic activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Population genomic analysis reveals differential evolutionary histories and patterns of diversity across subgenomes and subpopulations of Brassica napus L.

The allotetraploid species Brassica napus L. is a global crop of major economic importance, providing canola oil (seed) and vegetables for human consumption and fodder and meal for livestock feed. Characterizing the genetic diversity present in the extant germplasm pool of B. napus is fundamental to better conserve, manage and utilize the genetic resources of this species. We used sequence-based genotyping to identify and genotype 30,881 SNPs in a diversity panel of 782 B. napus accessions, representing samples of winter and spring growth habits originating from 33 countries across Europe, Asia, and America. We detected strong population structure broadly concordant with growth habit and geography, and identified three major genetic groups: spring (SP), winter Europe (WE), and winter Asia (WA). Subpopulation-specific polymorphism patterns suggest enriched genetic diversity within the WA group and a smaller effective breeding population for the SP group compared to WE. Interestingly, the two subgenomes of B. napus appear to have different geographic origins, with phylogenetic analysis placing WE and WA as basal clades for the other subpopulations in the C and A subgenomes, respectively. Finally, we identified 16 genomic regions where the patterns of diversity differed markedly from the genome-wide average, several of which are suggestive of genomic inversions. The results obtained in this study constitute a valuable resource for worldwide breeding efforts and the genetic dissection and prediction of complex B. napus traits.

58 GEOSCIENCES↗

Discovering the desirable alleles contributing to the lignocellulosic biomass traits in Saccharum germplasm collections for energy cane improvement

1) The success in crop improvement programs depends largely on the extent of genetic variability available. Germplasm collections assembles all the available genetic resources and are critical for long-term crop improvement. This world sugarcane germplasm collection contains enormous genetic variability for various morphological traits, biomass yield components, adaptation and many quality traits, prospectively imbeds a large number of valuable alleles for biofuel traits such as high biomass yield, quantity and quality of lignocelluloses, stress tolerance, and nutrient use efficiency. The germplasm collection is of little value unless it is characterized and utilized for crop improvement. In this project, we phenotypically and genotypically characterized the sugarcane world germplasm collection (The results were published in two papers already and another two papers are to be published). This data will be made available for public to refer to for germplasm unitization specifically in the sugarcane and energy cane breeding programs. In addition, we are identifying the alleles contributing to the biomass traits in sugarcane germplasm. This part of project is very challenging due to the large genome and highly polyploid level of this crop. We firstly established a high throughput sugarcane genotyping pipeline in the genome and bioinformatics era (a paper is published in 2016). We identified and modified a software for genome-wide association analysis of polyploid species. The results of the alleles associated to the biomass traits will be published soon, which will help the scientific community understand the genetic makeup of the biomass components of sugarcane. Molecular breeders can develop markers for marker assisted selection of biomass traits improvement. Further, the development and release of new energy cane cultivars through this project not only improved genetic diversity but also improved dry biomass yields and resistance to diseases. These new cultivars were tested on marginal soils in Florida and showed very promising yield potential that is important for the successful use of energy cane as a dedicated feedstock for lignocellulosic ethanol production. 2) Multiple techniques at different project progress stages were utilized. For example, for the whole world germplasm accession genotyping, a cheap widely used SSR marker genotyping platform was utilized due to the large number of samples (over thousand). But the throughput of this technique is low in generating data points. However, the purpose the genotyping is to form a core collection for further high throughput genotyping. Thus the results from the SSR genotyping was quite good enough to generated the core collection. To genotype the few hundred core collection accessions, an target enrichment sequencing technology was used, which is not only high throughput in generating large number of genotyping data, but also has the candidate genes targeted to genotyping. The data generated would be sufficient in identifying the alleles contributing to the traits of interests. All the techniques used in this project are effective though extensive time was invested specifically for establish the pipeline in the experimental design, data analysis, and different approach comparison. 3) the research can benefit to the public in polyploid genotyping and new and cost efficient genotyping platform development

59 BASIC BIOLOGICAL SCIENCES↗

Angular Distributions of High-Mass Dilepton Production in Hadron Collisions

The SeaQuest experiment is a fixed-target dimuon experiment currently running at the Fermi National Accelerator Laboratory (FNAL). By utilizing the high-intensity, 120 GeV proton beam delivered by the FNAL Main Injector (MI), SeaQuest is able to measure proton-induced Drell-Yan dimuon production off of various nuclear targets in kinematic regions inaccessible to previous similar experiments. A suitably large fraction of the final dataset has been recorded, reconstructed, and analyzed. Very preliminary results from light-sea flavor asymmetry, nuclear dependence, and partonic energy loss analyses have been presented at numerous international conferences. A novel, FPGA-based trigger system has been designed, implemented, and optimized for the SeaQuest experiment. By implementing the trigger decision logic in FPGA firmware, it is more adaptable to changing experimental conditions. Additionally, the peripheral tasks of timing alignment, “trigger matrix” generation, and firmware uploading have been mostly automated, reducing the likelihood of user error in the maintenance and operation of the trigger system. Significant upgrades to hardware and firmware have greatly improved the performance of the trigger system since the 2012 commissioning run of SeaQuest. Four additional v1495 modules were added to facilitate thorough pulser testing of the firmware designs and in-situ pulser tests of all compiled firmware. These pulser tests proved crucial for diagnosing many errors that may have otherwise gone unnoticed. A significant change to the internal clocking of the trigger system eliminated a subtle source of rate-dependent trigger efficiency. With this upgrade, the trigger finally meets the “dead-time free” design specification. Drell-Yan dimuon data have been collected and analyzed for central θCS , with nearly flat acceptance in φCS , in the mass range 5.0 GeV < Mγ* < 10.0 GeV at forward xF with the SeaQuest spectrometer at FNAL. A very preliminary extraction of λ has been performed, and the remaining difficulties in extracting ν have been evaluated. Although the results are not yet publishable, significant progress has been made in developing this very challenging angular distributions analysis. A simple scheme for correcting for the angular acceptances of the spectrometer, trigger, and reconstruction has been developed and demonstrated. A generally applicable correction for the kinematically-dependent, rate-dependent reconstruction efficiency has been developed and applied to all current analyses on SeaQuest data. This rate-dependence correction was the first major hurdle in the path to publication of many preliminary SeaQuest results. The last remaining major correction for all analyses, but especially important for the angular parameter extraction, is the full characterization, rate-dependence correction, and subtraction of the combinatoric background contribution to the reconstructed dimuon sample. Independently, an intuitive, kinematic derivation of the single-event definitions of the Drell-Yan angular parameters has been developed under the assumption of unpolarized annihilating quarks within unpolarized nuclei. At O(αs), where the quarks remain co-planar with the hadrons in the photon rest frame, this kinematic method reproduces the Lam-Tung relation and derives an additional equality for µ2, which is only interpretable for single-event parameters. This method has been extended to the case of quark non- coplanarity, and the coplanar equalities become inequalities. A new equality was discovered, which should be obeyed by single-event parameters even in the case of a non-coplanar quark axis. The non-coplanar parameter relations have been used to derive constraints on the experimentally accessible values of λ and ν. These constraints are compared with existing data and have been found consistent, except in the cases where significant contributions from non-zero Boer-Mulders functions are expected. Finally, the kinematically- derived parameter definitions have been applied to high-precision CMS data. The relative contributions of the qq¯ and qg processes to the Z-boson “Drell-Yan” cross-section have been extracted. Further, an average measure of non-coplanarity, likely caused by O(α2) and higher processes, has been extracted.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Pressure-Driven Fiber-Optic Sensor for Online Corrosion Monitoring

For many industrial applications, corrosion is a life-limiting phenomenon and therefore requires careful consideration and rigorous experimental testing before structural materials can be reliably deployed, particularly in harsh chemical environments. The traditional approach to measuring corrosion requires exposing many samples to the intended environment and then extracting them individually at discrete intervals for postexposure characterization. This approach does not provide a high degree of temporal resolution, nor does it provide any real-time information regarding dynamic changes in corrosion rates. Additionally, developing an online corrosion monitor capable of surviving harsh chemical environments could provide valuable information regarding the structural health of components and changing process conditions that could accelerate corrosion. To this end, a corrosion sensor was developed based on a pressure-driven Fabry-Pérot cavity (FPC). This sensor uses a pressure control system to internally pressurize the FPC formed between the sensor’s housing and a metal-embedded singlemode optical fiber. Simultaneous measurement of the change in FPC length using low-coherence interferometry and the applied pressure enables the calculation of the relative changes in the sensor’s diaphragm thickness due to corrosion on its outer surface. Measurements were made in situ while actively corroding the sensor and were validated against surveillance specimens that were corroded simultaneously. The uncertainties in the measured corrosion were analyzed using error propagation of the uncertainties in the measured pressures and displacements and were found to be <1% of the initial diaphragm thickness.

42 ENGINEERING↗

MPEX AI Digital Twins

All magnetically confined plasma fusion power plant concepts (Tokamak, Spherical Tokamak, Stellarator, Mirror, ...) must exhaust the heat and plasma from the core confinement region to the material walls. The primary channel for this exhaust is through a plasma divertor which directs plasma along open magnetic field lines to a material target. The Material Plasma Exposure eXperiment (MPEX) illustrated in Figure 1, is a high-power, steady-state linear plasma device designed to produce the plasma material interaction (PMI) conditions of the divertor of future magnetic confinement fusion power plants: energy flux 20MW/m 2 , ion fluence 1031/m 2 , pulse duration 106 sec. These goals of plasma exposure in MPEX are well beyond those achieved in magnetic fusion experimental devices. Successfully achieving these high power steady state conditions for long pulses requires operational control of the heating and particle sources and the plasma flux to the walls and target. The MPEX AI Hot Spot Controller, proposed in this project, will help achieve the operational milestones of MPEX. The MPEX device will begin commissioning at the end of FY26. A smaller proto-MPEX was operated for 14,666 plasma discharges and will resume operation in September of 2025 as proto-MPEX-lite, with reduced capability, to test a new window for the Helicon plasma source. The proto-MPEX data has undergone surrogate modeling with machine learning methods (R. Archibald, 2022 IEEE International Conference on Big Data). This proto-MPEX data will be used to begin development of the AI digital twins described in this white paper. The scientific mission of MPEX is to qualify materials of different composition for use in the high energy and plasma flux conditions of a fusion power plant. The materials exposed in MPEX will in some cases be exposed to high neutron fluxes at other ORNL facilities to measure the changes to their PMI properties. The targets exposed in MPEX will be transported under vacuum to a Surface Analysis Station (SAS). The SAS will be equipped with the following diagnostics: Focused Ion Beam (FIB) for trench milling, 100-400 angstrom resolution scanning electron microscope (SEM), surface mapping x-ray spectrometer, high resolution camera, and a future upgrade to a laser induced breakdown spectroscopy quadruple mass spectrometer (LIBS-QMS). The MPEX experiments will generate diverse pre- and post-exposure measurement data of detailed material properties down to the crystal grain level in 3D for post-exposure assessment of PMI damage (e.g. cracking, melting, erosion and redeposition of the material). Physics models for the PMI, and how the material composition and manufacturing impact its performance under high energy plasma exposure, need to be validated with MPEX data to guide the selection of new candidate materials. Our vision for the MPEX AI Digital Twins project is to supply experimental and physics model simulation data to train Artificial Intelligence (AI) models for data processing, analysis, operational control, PMI and materials simulation to maximize the scientific output of the MPEX device. Ultimately, an AI digital twin of MPEX material assessment metrics for tested and synthetic material types with simulated PMI will be trained by the AI Modeling Teams on the experimental and physics simulation data submitted to the American Science Cloud by this project. A purely empirical search for the best material is inefficient given the finite number of samples that can be tested on MPEX. In order to expand the material properties database for training the MPEX Material Assessment AI Digital Twin, and to gain physics understanding of the PMI processes, physics models of the material properties and PMI processes are required. The physics simulations provide detailed simulation data, like impact angles for plasma ions, sputtering yields, transport of the ionized sputtered target material in the plasma, and redeposition locations. This simulation data expands the measurement data for deeper physics understanding. The experimental data is essential to validate the PMI and material structure simulation models. The validated models can then be used to generate new simulation data of MPEX material assessments for synthetic material compositions that have not been exposed in MPEX. These predictive simulations, plus the whole experimental dataset, will be used to train the MPEX Material Assessment AI Digital Twin allowing a rapid generative AI search for new materials with reduced PMI damage by interpolating the domain of the training set. These new optimum materials can be simulated with the physics codes and/or tested in MPEX. The ability of AI neural networks to interpolate multi-dimensional parameter spaces and generate virtual data is exploited for a more efficient search for optimum materials. The advent of the Transformational AI Models Consortium (TAIMC) is an opportunity to engage with state of the art private and public AI developers to achieve the goals of the AI digital twins and AI accelerated physics models proposed in this project. Our partners at ORNL from the Advance Scientific Computing Research (ASCR) organization will collaborate in accelerating the integrated plasma material interaction simulation framework. This simulation framework will provide a platform for generating simulation data across a range of physical fidelities, including hybrid methods that produce multi-fidelity results. This data will be leveraged for AI model development, both for generation of surrogates and the automation of simulation campaigns. A part of the research below will include collaborative efforts with the TAIMC to (i) adapt data storage approaches to ensure AI-readiness, (ii) provide a protypical exemplar to inform and exercise constructed workflows, and (iii) generate and share data, using the TAIMC unified AI data standard, for foundational models that will be trained from multiple sources across the DOE complex. We will also collaborate with the TAIMC, as well as the planned AI modeling teams, to develop approaches for reducing the cost of data generation. These include tailored multi-fidelity approaches as well as fine-tuning strategies to augment general, large-scale foundational models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exploratory Project: Rigid nanostructured organic polymer monolith for in situ collection and analysis of plant metabolites from soil matrices

Plant metabolites released from litter leachates and root exudates enable plants to adapt and survive in a wide range of habitats by facilitating resource foraging and plant-organismal interactions, and could influence soil carbon storage. The biological functions of these plant inputs and the organismal interactions they facilitate in soil are strictly governed by their composition and molecular identity. Our current understanding about the molecular identity of exudates is based on physiological studies that are done in soil-less axenic cultures. On the other hand, ecological studies that rely on isotope labeling to track the fluxes of carbon from plants to soil have found the complexities of soil-microbe matrices as an insurmountable barrier to undertake any meaningful molecular level characterization of plant inputs. Although it is constantly advocated to undertake a molecular level identification of the dynamicity of plant metabolites in soils, the complexity of soil system has thus far prevented any such endeavors. We developed polymeric probes through in-situ polymerization of poly(styrene-co-vinylbenzyl chloride-co-divinylbenzene) to elucidate the chemical environment of the soil to which the plant roots are exposed. Hypercrosslinking of the polymeric probes through an in-situ Friedel-Crafts alkylation significantly increased the surface area and the sorption capacity of the probes. Surface functionalization of the probes using a hybrid approach was also attempted. The efficacy of these probes was tested using batch equilibration. Scanning electron microscopy revealed extensive modification of the surface of the probes through hypercrosslinking. The probes exhibited a lower site specific sorption (slope of Freundlich adsorption isotherm close to unity) and percent recovery of the sorbed compounds from the probes were >70, indicating a predominance of reversible sorption. Further we imparted specificity to this copolymer matrix by using molecular-template imprinting approach for the selective capture of less abundant plant metabolites from a crowded soil/litter leachate. Our results suggests that i) the root exudate pattern of plants is highly dependent on the nutrient status of the plant, with greater specificity of root exudation occurring in growing medium with low available form of nutrients, ii) the chemical composition of root exudation is a function of the distance of sampling from the source-roots, with the composition of root exudation being more enriched in polar metabolites farther from the source-roots, iii) Further we demonstrated that the compounds present in real root exudates diffuse farther from the source roots than those in artificial root exudates that are traditionally used. Thus, our project highlights how the soil matrix is instrumental in modifying the chemical composition of root exudates, and highlights that, apart from the plant physiology, the specificity and function of root exudates is also modified by environmental factors.

54 ENVIRONMENTAL SCIENCES↗

PROPER MOTIONS OF THE ARCHES CLUSTER WITH KECK LASER GUIDE STAR ADAPTIVE OPTICS: THE FIRST KINEMATIC MASS MEASUREMENT OF THE ARCHES

We report the first detection of the intrinsic velocity dispersion of the Arches cluster-a young ({approx}2 Myr), massive (10{sup 4} M{sub Sun }) starburst cluster located only 26 pc in projection from the Galactic center. This was accomplished using proper motion measurements within the central 10'' Multiplication-Sign 10'' of the cluster, obtained with the laser guide star adaptive optics system at Keck Observatory over a three-year time baseline (2006-2009). This uniform data set results in proper motion measurements that are improved by a factor {approx}5 over previous measurements from heterogeneous instruments. By careful, simultaneous accounting of the cluster and field contaminant distributions as well as the possible sources of measurement uncertainties, we estimate the internal velocity dispersion to be 0.15 {+-} 0.01 mas yr{sup -1}, which corresponds to 5.4 {+-} 0.4 km s{sup -1} at a distance of 8.4 kpc. Projecting a simple model for the cluster onto the sky to compare with our proper motion data set, in conjunction with surface density data, we estimate the total present-day mass of the cluster to be M(r < 1.0 pc) = 1.5{sup +0.74}{sub -0.60} Multiplication-Sign 10{sup 4} M{sub Sun }. The mass in stars observed within a cylinder of radius R (for comparison to photometric estimates) is found to be M(R < 0.4 pc) = 0.90{sup +0.40}{sub -0.35} Multiplication-Sign 10{sup 4} M{sub Sun} at formal 3{sigma} confidence. This mass measurement is free from assumptions about the mass function of the cluster, and thus may be used to check mass estimates from photometry and simulation. Photometric mass estimates assuming an initially Salpeter mass function ({Gamma}{sub 0} = 1.35, or {Gamma} {approx} 1.0 at present, where dN/d(log M){proportional_to}M{sup {Gamma}}) suggest a total cluster mass M{sub cl} {approx} (4-6) Multiplication-Sign 10{sup 4} M{sub Sun} and projected mass ({approx} 2 {<=} M(R < 0.4 pc) {<=} 3) Multiplication-Sign 10{sup 4} M{sub Sun }. Photometric mass estimates assuming a globally top-heavy or strongly truncated present-day mass function (PDMF; with {Gamma} {approx} 0.6) yield mass estimates closer to M(R < 0.4 pc) {approx} 1-1.2 Multiplication-Sign 10{sup 4} M{sub Sun }. Consequently, our results support a PDMF that is either top-heavy or truncated at low mass, or both. Collateral benefits of our data and analysis include: (1) cluster membership probabilities, which may be used to extract a clean-cluster sample for future photometric work; (2) a refined estimate of the bulk motion of the Arches cluster with respect to the field, which we find to be 172 {+-} 15 km s{sup -1}, which is slightly slower than suggested by previous measurements using one epoch each with the Very Large Telescope and the Keck telescope; and (3) a velocity dispersion estimate for the field itself, which is likely dominated by the inner Galactic bulge and the nuclear disk.

79 ASTRONOMY AND ASTROPHYSICS↗

Final Scientific/Technical Report

A variety of calcifying organisms produce a transient or metastable amorphous calcium carbonate (ACC) precursor phase that is assembled and subsequently transformed into a crystalline biomineral, typically calcite or aragonite. The complex shapes, hierarchical structures, and unique physical properties of the biominerals that result from this calcification pathway have stimulated interest in adapting these concepts for the design and creation of bio-inspired functional materials in the laboratory. ACC also forms as a reactive precursor in diverse inorganic systems and is likely to play a much broader role in calcium carbonate formation. Knowledge of the structure, composition, and behavior of this metastable phase is critical for establishing a structural and mechanistic framework for calcium carbonate formation and its role in biogeochemical processes, including carbon cycling. Minor additives, such as magnesium, phosphorus, and organic macromolecules, are known to play important roles in controlling ACC stability, transformation kinetics, and selection of final crystalline polymorph. Molecular water also occurs in many types of ACC and is thought to play a structural role in its stability and transformation behavior. One of the major challenges that remain unresolved is identification of the structural basis for the role of these minor additives and molecular water. The absence of long-range order in ACC, and other amorphous phases, has posed a challenge for study by techniques commonly used for crystalline solids. Preliminary studies in our group show that the combination of two techniques, synchrotron X-ray-based pair distribution function (PDF) analysis and nuclear magnetic resonance (NMR) spectroscopy can provide entirely new insight to structural properties of synthetic ACC over length scales that are most relevant for understanding its transformation properties. Building on preliminary experiments, we propose a systematic study of synthesis, structure, and transformation behavior in abiotic systems. The work will specifically address the influence of phosphate as a minor additive. PDF analysis will utilize total X-ray scattering data collected at synchrotron facilities optimized for this method and will provide direct characterization of the short- and intermediate-range structure of ACC synthesized under controlled conditions. Parallel computational work using reverse Monte Carlo methods will allow structural models to be constructed for a more complete analysis of PDF results. NMR spectroscopy, using a variety of single- and double-resonance techniques, will provide information on H and CO3 components, including dynamical properties, and their relationship to stabilizing additives. PDF and NMR results will be complemented by parallel studies using X-ray absorption and FT-IR spectroscopy to allow direct comparison to previous studies. These techniques will be used to follow the transformation of different ACC samples, with specific additives, to crystalline phases under controlled condition (e.g., relative humidity). This work will provide the structural and mechanistic basis for understanding ACC stability, its transformation behavior, and the factors that govern polymorph selection. This new insight will directly benefit researchers in diverse fields, as well as adding to the framework of knowledge for understanding and controlling calcium carbonate formation in natural and engineered systems.

58 GEOSCIENCES↗

Shed Light in the DaRk LineagES of the Fungal Tree of Life—STRES

The polyphyletic group of black fungi within the Ascomycota (Arthoniomycetes, Dothideomycetes, and Eurotiomycetes) is ubiquitous in natural and anthropogenic habitats. Partly because of their dark, melanin-based pigmentation, black fungi are resistant to stresses including UV- and ionizing-radiation, heat and desiccation, toxic metals, and organic pollutants. Consequently, they are amongst the most stunning extremophiles and poly-extreme-tolerant organisms on Earth. Even though ca. 60 black fungal genomes have been sequenced to date, [mostly in the family Herpotrichiellaceae (Eurotiomycetes)], the class Dothideomycetes that hosts the largest majority of extremophiles has only been sparsely sampled. By sequencing up to 92 species that will become reference genomes, the “Shed light in The daRk lineagES of the fungal tree of life” (STRES) project will cover a broad collection of black fungal diversity spread throughout the Fungal Tree of Life. Interestingly, the STRES project will focus on mostly unsampled genera that display different ecologies and life-styles (e.g., ant- and lichen-associated fungi, rock-inhabiting fungi, etc.). With a resequencing strategy of 10- to 15-fold depth coverage of up to ~550 strains, numerous new reference genomes will be established. To identify metabolites and functional processes, these new genomic resources will be enriched with metabolomics analyses coupled with transcriptomics experiments on selected species under various stress conditions (salinity, dryness, UV radiation, oligotrophy). The data acquired will serve as a reference and foundation for establishing an encyclopedic database for fungal metagenomics as well as the biology, evolution, and ecology of the fungi in extreme environments.

59 BASIC BIOLOGICAL SCIENCES↗

GENPLAT: an Automated Platform for Biomass Enzyme Discovery and Cocktail Optimization

The high cost of enzymes for biomass deconstruction is a major impediment to the economic conversion of lignocellulosic feedstocks to liquid transportation fuels such as ethanol. We have developed an integrated high throughput platform, called GENPLAT, for the discovery and development of novel enzymes and enzyme cocktails for the release of sugars from diverse pretreatment/biomass combinations. GENPLAT comprises four elements: individual pure enzymes, statistical design of experiments, robotic pipeting of biomass slurries and enzymes, and automated colorimeteric determination of released Glc and Xyl. Individual enzymes are produced by expression in Pichia pastoris or Trichoderma reesei, or by chromatographic purification from commercial cocktails or from extracts of novel microorganisms. Simplex lattice (fractional factorial) mixture models are designed using commercial Design of Experiment statistical software. Enzyme mixtures of high complexity are constructed using robotic pipeting into a 96-well format. The measurement of released Glc and Xyl is automated using enzyme-linked colorimetric assays. Optimized enzyme mixtures containing as many as 16 components have been tested on a variety of feedstock and pretreatment combinations. GENPLAT is adaptable to mixtures of pure enzymes, mixtures of commercial products (e.g., Accellerase 1000 and Novozyme 188), extracts of novel microbes, or combinations thereof. To make and test mixtures of ˜10 pure enzymes requires less than 100 μg of each protein and fewer than 100 total reactions, when operated at a final total loading of 15 mg protein/g glucan. We use enzymes from several sources. Enzymes can be purified from natural sources such as fungal cultures (e.g., Aspergillus niger, Cochliobolus carbonum, and Galerina marginata), or they can be made by expression of the encoding genes (obtained from the increasing number of microbial genome sequences) in hosts such as E. coli, Pichia pastoris, or a filamentous fungus such as T. reesei. Proteins can also be purified from commercial enzyme cocktails (e.g., Multifect Xylanase, Novozyme 188). An increasing number of pure enzymes, including glycosyl hydrolases, cell wall-active esterases, proteases, and lyases, are available from commercial sources, e.g., Megazyme, Inc. (www.megazyme.com), NZYTech (www.nzytech.com), and PROZOMIX (www.prozomix.com). Design-Expert software (Stat-Ease, Inc.) is used to create simplex-lattice designs and to analyze responses (in this case, Glc and Xyl release). Mixtures contain 4-20 components, which can vary in proportion between 0 and 100%. Assay points typically include the extreme vertices with a sufficient number of intervening points to generate a valid model. In the terminology of experimental design, most of our studies are "mixture" experiments, meaning that the sum of all components adds to a total fixed protein loading (expressed as mg/g glucan). The number of mixtures in the simplex-lattice depends on both the number of components in the mixture and the degree of polynomial (quadratic or cubic). For example, a 6-component experiment will entail 63 separate reactions with an augmented special cubic model, which can detect three-way interactions, whereas only 23 individual reactions are necessary with an augmented quadratic model. For mixtures containing more than eight components, a quadratic experimental design is more practical, and in our experience such models are usually statistically valid. All enzyme loadings are expressed as a percentage of the final total loading (which for our experiments is typically 15 mg protein/g glucan). For "core" enzymes, the lower percentage limit is set to 5%. This limit was derived from our experience in which yields of Glc and/or Xyl were very low if any core enzyme was present at 0%. Poor models result from too many samples showing very low Glc or Xyl yields. Setting a lower limit in turn determines an upper limit. That is, for a six-component experiment, if the lower limit for each single component is set to 5%, then the upper limit of each single component will be 75%. The lower limits of all other enzymes considered as "accessory" are set to 0%. "Core" and "accessory" are somewhat arbitrary designations and will differ depending on the substrate, but in our studies the core enzymes for release of Glc from corn stover comprise the following enzymes from T. reesei: CBH1 (also known as Cel7A), CBH2 (Cel6A), EG1(Cel7B), BG (βglucosidase), EX3 (endo-β1,4-xylanase, GH10), and BX (β-xylosidase).

09 BIOMASS FUELS↗

Guidance for Monitoring Passive Groundwater Remedies Over Extended Time Scales

Passive remediation can be appropriate where natural processes and actions such as institutional controls mitigate exposure to contaminated groundwater, achieving remedial action objectives and protectiveness of human health and the environment. Monitored natural attenuation (MNA) is a prevalent passive remediation strategy supported by a regulatory framework and monitoring design guidance. However, long-term passive remedies are usually selected in combination with at least one active remedy, such as source removal, in situ treatment, or pump-and-treat, functioning as a complementary method for achieving remediation objectives and meeting the applicable statutory and regulatory requirements. However, MNA and existing monitoring guidance primarily target situations where the remedial action objectives are met within a few decades. When time scales for passive remediation extend to many decades (50 years or more), a corresponding change in monitoring strategy is needed to adapt to the extended time scale. This document provides guidance for implementing an extended-scale monitoring (ESM) approach appropriate for long-duration passive remediation. Extended-scale is defined in this document with respect to time (i.e., a long duration of remediation) and a large enough physical scale such that downstream receptors will not be impacted within the remediation timeframe. ESM applies to slow-moving groundwater contaminant plumes and emphasizes monitoring primarily for potential exposure pathways. For this approach, the primary monitoring objective is to demonstrate that the plume diminishes before reaching the receptor zone or point of compliance and/or a receptor does not receive concentrations above the compliance limit. While the overall objectives of protecting human health and the environment are the same as for plumes where remediation can occur over a shorter time period, the time scales between decisions are longer and the dynamics of plume evolution are slower. To this end, a scenario-based strategy is described for different plume and source conditions, defining a containment and receptor zones. The containment zone is the area where the risk of exposure to groundwater contamination can be mitigated (e.g., through institutional controls) during the remediation time period. The receptor zone is defined as the area where exposure to groundwater contamination cannot be mitigated and compliance concentration standards must be met. Within the containment zone, slow plume migration may occur, leading to concentrations that exceed compliance standards. However, where distance to the receptor zone is large relative to plume migration and attenuation rate, this approach can be protective of the receptor zone. Selection of a long-duration passive remedy needs to be based on sufficient understanding of contaminant sources, hydrogeology, and contaminant plumes. A strong technical basis, supported by predictive analysis, is recommended to substantiate that contamination is expected to stay within the containment zone during the active remediation and attainment phase of the remedy, and diminish to meet compliance standards within the extended timeframe prior to reaching the receptor zone (e.g., many decades or even centuries). The ESM approach is based on verification of plume behavior and not on detailed plume dynamics. Monitoring is conducted to confirm expected behavior with an emphasis on exposure pathways to verify that plumes remain contained in areas where the protectiveness objectives can be met. ESM should not be adopted if there is significant risk of the plume extending beyond the containment zone. Given the slow movement within the containment zone, less frequent sampling is required relative to approaches used for conventional-scale remediation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Final Technical Report to DOE for the Award DE-SC0004601

Understanding the responses, adaptations and feedback mechanisms of biological communities to climate change is critical to project future state of earth and climate systems. Although significant amount of knowledge is available on the feedback responses of aboveground communities to climate change, little is known about the responses of belowground microbial communities due to the challenges in analyzing soil microbial community structure. Thus the goal overall goal of this study is to provide system-level, predictive mechanistic understanding of the temperature sensitivity of soil carbon (C) decomposition to climate warming by using cutting-edge integrated metagenomic technologies. Towards this goal, the following four objectives will be pursued: (i) To determine phylogenetic composition and metabolic diversity of microbial communities in the temperate grassland and tundra ecosystems; (ii) To delineate the responses of microbial community structure, functions and activities to climate change in the temperate grassland and tundra ecosystems; (iii) To determine the temperature sensitivity of microbial respiration in soils with different mixtures of labile versus recalcitrant C, and the underlying microbiological basis for temperature sensitivity of these pools; and (iv) To synthesize all experimental data for revealing microbial control of ecosystem carbon processes in responses to climate change. We have achieved our goals for all four proposed objectives. First, we determined the phylogenetic composition and metabolic diversity of microbial communities in the temperate grassland and tundra ecosystems. For this objective, we have developed a novel phasing amplicon sequencing (PAS) approach for MiSeq sequencing of amplicons. This approach has been used for sequencing various phylogenetic and functional genes related to ecosystem functioning. A comprehensive functional gene array (e.g., GeoChip 5.0) has also been developed and used for soil microbial community analysis in this study. In addition, shot-gun metagenome sequencing along with the above approaches have been used to understand the phylogenetic and functional diversity, composition, and structure of soil microbial communities in both temperature grassland and tundra ecosystems. Second, we determined the response of soil microbial communities to climate warming in both temperate grassland and tundra ecosystems using various methods. Our major findings are: (i) Microorganisms are very rapid to respond to climate warming in the tundra ecosystem, AK, which is vulnerable, too. (ii) Climate warming also significantly shifted the metabolic diversity, composition and structure of microbial communities, and key metabolic pathways related to carbon turnover, such as cellulose degradation (~13%) and CO2 production (~10%), and to nitrogen cycling, including denitrification (~12%) were enriched by warming. (iii) Warming also altered the expression patterns of microbial functional genes important to ecosystem functioning and stability through GeoChip and metatranscriptomic analysis of soil microbial communities at the OK site. Third, we analyzed temperature sensitivity of C decomposition to climate warming for both AK and OK soils through laboratory incubations. Key results include: (i) Alaska tundra soils showed that after one year of incubation, CT in the top 15 cm could be as high as 25% and 15% of the initial soil C content at 25°C and 15°C incubations, respectively. (ii) analysis of 456 incubated soil samples with 16S rRNA gene, ITS and GeoChip hybridization showed that warming shifted the phylogenretic and functional diversity, composition, structure and metabolic potential of soil microbial communities, and at different stages of incubation, key populations and functional genes significantly changed along with soil substrate changes. Functional gene diversity and functional genes for degrading labile C components decrease along incubation when labile C components are exhausting, but the genes related to degrading recalcitrant C increase. These molecular data will be directly used for modeling. Fourth, we have developed novel approaches to integrate and model experimental data to understand microbial control of ecosystem C processes in response to climate change. We compared different methods to calculate Q10 for estimating temperature sensitivity, and new approaches for Q10 calculation and molecular ecological network analysis were also developed. Using those newly developed approaches, our result indicated that Q10s increased with the recalcitrance of C pools, suggesting that longer incubation studies are needed in order to assess the temperature sensitivity of slower C pools, especially at low temperature regimes. This project has been very productive, resulting in 42 papers published or in press, 4 submitted, and 13 in preparation.

54 ENVIRONMENTAL SCIENCES↗

Probability of Pinching Failure in PWR SNF after Long-Term Storage

A probabilistic assessment was performed on pinching failure of PWR SNF rod cladding under accidental SNF cask transport conditions. a set of 3,000 rod pinching scenarios was analyzed, with different geometrical, material, and loading parameters, using a Latin Hypercube Sampling (LHS). The investigated accidents include train derailments leading to cask collisions or cask drops. Pinching failure was related to cladding offset strain, defined as the plastic deformation interval in radial cladding direction normalized with respect to the outer cladding diameter. Finite element (FE) models were used to compute SNF rod deformations and the offset strain demand OSD in each of the assessed pinching scenarios. The numerical models account for the contribution of the degraded fuel pellets to rod pinching stiffness. The offset strain capacity OSC was computed using regression and interpolation models based on ring compression (RCT) test results, and account for hydride-related effects on cladding ductility. The results show a low probability of pinching failure. Hydride-related cladding embrittlement affects the structural integrity of PWR SNF rods after long-term dry storage of up to 300 years. Radial hydride patterns or large amounts of circumferentially oriented hydrides reduce the cladding ductile deformation capability when pinched at low material temperatures. The probability of hydride-related pinching failure of PWR cladding is estimated after long-term dry storage to address: How do PWR SNF rod conditions like hydrogen content, vacuum drying peak cladding hoop stress σ{sub θ}, and cladding material temperature affect cladding ductility? Is hydride-related cladding embrittlement a material degradation mechanism that affects the structural integrity of PWR SNF rods under pinching loading? What are the controlling parameters for deformation demands in PWR SNF rods under pinching loading? Is the pinching load amplitude an important parameter that controls the probability of pinching failure? How do parameters like pellet-cladding gap or fuel pellet degradation affect the pinching loading capacity of PWR SNF rods? Statistical modelling of PWR SNF rod conditions in U.S. inventory: A pseudo rod database was generated based on EIA database. Analysis of train accident scenarios: a probability density function (pdf) of a pinching load distribution P{sub 9m} of a 9-m cask drop is linearly scaled to generate a more representative: distribution P{sub AS}. The new distribution considers train collision velocities, cask drops for on- and off-bridge accidents, and cask impact orientation. LHS generates 3,000 rod pinching scenarios based on nine independent geometrical, material, and loading parameters. Probabilistic FE Analysis: The LHS realizations are used to generate 3,000 FE models. The FE analysis is conducted for P{sub 9m}, and P{sub AS}. The workflow is automatized. The OSC is computed with regression and interpolation models, as a function of the hydrogen content C{sub H}, vacuum drying peak σ{sub θ}, and SNF temperature T{sub FR}: The regression data was obtained in RCTs on PWR cladding. The OSD is derived from the cladding deformations in the FE models. Cladding failure was identified via a comparison of OSC and OSD. A 300 and a 30 years-of-dry-storage duration is investigated: The dry-storage period affects the expected SNF temperature and consequently, the OSC and the OSD. The material temperature LHS population is adapted to account for a SNF temperature variation. The probability of hydride-related, ductile cladding failure during a severe transportation accident is not significantly affected by the pinching load amplitude: The probability is slightly reduced from a range of less than 0.0333% to 0.1667% for P{sub 9m}, to a range of less than 0.0333% to 0.1333% for P{sub AS}, depending on the vacuum drying peak cladding hoop stress modelling assumptions. Further, the SNF dry-storage duration does not significantly affect the probability of ductile pinching failure. Based on available train accidents data, the pinching failure probability for a 5,000 km long transportation route ranges from less than 1.62 E-7 to a value of 6.47 E-7, depending on the vacuum drying peak cladding hoop stress modelling assumptions. The pinching load amplitude is not a controlling parameter for the probability of ductile cladding pinching failure, and does not control the cladding deformations once pellet-cladding contact is established, The pellet-cladding gap width, however, plays a key role. Fuel pellet degradation only affects the deformation capacity, if the pellet stiffness is significantly reduced (more than 90%). Only highly embrittled cladding is likely to fail. The deformation demands on the cladding under pinching loading are small.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗