Search NASA⌕ Search

SEARCH · Search NASA

Results for “structural dynamics testing”

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,243 records · Page 69

Diffuse scattering from dynamically compressed single-crystal zirconium following the pressure-induced α → ω phase transition

The prototypical α → ω phase transition in zirconium is an ideal test bed for our understanding of polymorphism under extreme loading conditions. After half a century of study, a consensus had emerged that the transition is realized via one of two distinct displacive mechanisms, depending on the nature of the compression path. However, recent dynamic-compression experiments equipped with diffraction diagnostics performed in the past few years have revealed new transition mechanisms, demonstrating that our understanding of the underlying atomistic dynamics and transition kinetics is in fact far from complete. We present classical molecular dynamics simulations of the α → ω phase transition in single-crystal zirconium shock compressed along the [0001] axis using a machine-learning-class potential. The transition is predicted to proceed primarily via a modified version of the two-stage Usikov-Zilberstein mechanism, whereby the high-pressure ω phase heterogeneously nucleates at boundaries between grains of an intermediate β phase. We further observe the fomentation of atomistic disorder at the junctions between β grains, leading to the formation of highly defective interstitial material between the ω grains. We directly compare synthetic x-ray diffraction patterns generated from our simulations with those obtained using femtosecond diffraction in recent dynamic-compression experiments, and show that the simulations produce the same unique, anisotropic diffuse scattering signal unlike any previously seen from an elemental metal. Our simulations suggest that the diffuse signal arises from a combination of thermal diffuse scattering, nanoparticlelike scattering from residual kinetically stabilized α and β grains, and scattering from interstitial defective structures. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗

Testing $f(R)$ gravity with scale dependent cosmic void velocity profiles

Here, we study the impact of cosmological scale modifications to general relativity on the dynamics of halos within voids by comparing N-body simulations incorporating Hu-Sawicki f(R) gravity, with |f R0 | = 10 -6 and 10 -5 , to those of Λ CDM. By examining the radial velocity statistics within voids classified based on their size and density-profile, as “rising” (R-type) or “shell” (S-type), we find that halo motions in small R-type voids, with effective radius <15 Mpc/h, reveal distinctive differences between f(R) and Λ CDM cosmologies. To understand this observed effect, we study the linear and nonlinear fifth forces, and develop an iterative algorithm to accurately solve the nonlinear fifth force equation. We use this to characterize the Chameleon screening mechanism in voids and contrast the behavior with that observed in gravitationally collapsed objects. The force analysis underscores how smaller R-type voids exhibit the highest ratios of fifth force to Newtonian force, which source distinguishable differences in the velocity profiles and thereby provide rich environments in which to constrain gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Validating Molecular Dynamics Simulations against Experimental Observables in Light of Underlying Conformational Ensembles

Far from the static, idealized conformations deposited into structural databases, proteins are highly dynamic molecules that undergo conformational changes on temporal and spatial scales that may span several orders of magnitude. These conformational changes, often intimately connected to the functional roles that proteins play, may be obscured by traditional biophysical techniques. Over the past 40 years, molecular dynamics (MD) simulations have complemented these techniques by providing the “hidden” atomistic details that underlie protein dynamics. However, there are limitations of the degree to which molecular simulations accurately and quantitatively describe protein motions. Here we show that although four molecular dynamics simulation packages (AMBER, GROMACS, NAMD, and il mm) reproduced a variety of experimental observables for two different proteins (engrailed homeodomain and RNase H) equally well overall at room temperature, there were subtle differences in the underlying conformational distributions and the extent of conformational sampling obtained. This leads to ambiguity about which results are correct, as experiment cannot always provide the necessary detailed information to distinguish between the underlying conformational ensembles. Furthermore, the results with different packages diverged more when considering larger amplitude motion, for example, the thermal unfolding process and conformational states sampled, with some packages failing to allow the protein to unfold at high temperature or providing results at odds with experiment. While most differences between MD simulations performed with different packages are attributed to the force fields themselves, there are many other factors that influence the outcome, including the water model, algorithms that constrain motion, how atomic interactions are handled, and the simulation ensemble employed. Here four different MD packages were tested each using best practices as established by the developers, utilizing three different protein force fields and three different water models. Differences between the simulated protein behavior using two different packages but the same force field, as well as two different packages with different force fields but the same water models and approaches to restraining motion, show how other factors can influence the behavior, and it is incorrect to place all the blame for deviations and errors on force fields or to expect improvements in force fields alone to solve such problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗

A comparison of methods used to unveil the genetic and metabolic pool in the built environment

BACKGROUND:A majority of indoor residential microbes originate from humans, pets, and outdoor air and are not adapted to the built environment (BE). Consequently, a large portion of the microbes identified by DNA-based methods are either dead or metabolically inactive. Although many exceptions have been noted, the ribosomal RNA fraction of the sample is more likely to represent either viable or metabolically active cells. We examined methodological variations in sample processing using a defined, mock BE microbial community to better understand the scope of technique-based vs. biological-based differences in both ribosomal transcript (rRNA) and gene (DNA) sequence community analysis. Based on in vitro tests, a protocol was adopted for the analysis of the genetic and metabolic pool (DNA vs. rRNA) of air and surface microbiomes within a residential setting. RESULTS:We observed differences in DNA/RNA co-extraction efficiency for individual microbes, but overall, a greater recovery of rRNA using FastPrep (> 50%). Samples stored with various preservation methods at - 80°C experienced a rapid decline in nucleic acid recovery starting within the first week, although post-extraction rRNA had no significant degradation when treated with RNAStable. We recommend that co-extraction samples be processed as quickly as possible after collection. The in vivo analysis revealed significant differences in the two components (genetic and metabolic pool) in terms of taxonomy, community structure, and microbial association networks. Rare taxa present in the genetic pool showed higher metabolic potential (RNA:DNA ratio), whereas commonly detected taxa of outdoor origins based on DNA sequencing, especially taxa of the Sphingomonadales order, were present in lower relative abundances in the viable community. CONCLUSIONS:Although methodological variations in sample preparations are high, large differences between the DNA and RNA fractions of the total microbial community demonstrate that direct examination of rRNA isolated from a residential BE microbiome has the potential to identify the more likely viable or active portion of the microbial community. In an environment that has primarily dead and metabolically inactive cells, we suggest that the rRNA fraction of BE samples is capable of providing a more ecologically relevant insight into the factors that drive indoor microbial community dynamics.

59 BASIC BIOLOGICAL SCIENCES↗

Software Quality Assurance Plan ANSYS LSDYNA Version 2023R1

ANSYS Inc. develops and markets engineering simulation software and services used in the aerospace, automotive, manufacturing, electronics, biomedical, energy, defense, and many other industries. ANSYS is dedicated to engineering simulation and is the world’s leading software provider. ANSYS was founded in 1970 and is headquartered in Canonsburg, Pennsylvania. ANSYS provides an engineering analysis tool combining structural, thermal, computational fluid dynamics, acoustic and electromagnetic simulation capabilities. ANSYS LS-DYNA is the most used explicit simulation program capable of simulating the response of materials to short periods of severe loading. Its many elements, contact formulations, material models, and other controls can be used to simulate complex models with control over all the details of the problem. ANSYS LS-DYNA has a vast array of capabilities to simulate extreme deformation problems using its explicit solver. Engineers can tackle simulations involving material failure and look at how the failure progresses through a part or through a system. Models with large amounts of parts or surfaces interacting with each other are also easily handled, and the interactions and load passing between complex behaviors are modeled accurately. Using computers with higher numbers of CPU cores can drastically reduce solution times. In addition, many consulting firms and hundreds of universities use ANSYS for analysis, research, and educational purposes. ANSYS is recognized worldwide as one of the most widely used and capable programs of its type. ANSYS has successfully passed over 100 customer quality system audits against American Society of Mechanical Engineers (ASME) NQA-1 and 10 CFR Part 50, Appendix B, since the company was founded, over 60 of which have been since 1997. ANSYS has successfully passed over 100 International Organization for Standardization (ISO) 9001 assessments. ANSYS design analysis software is the first created within a quality system with ISO 9001 certification, which is the internationally accepted quality standard. Product development, testing, maintenance, and support processes also meet the US Nuclear Regulatory Commission’s (NRC’s) quality requirements, as they have for nearly four decades. ANSYS staff perform more than 60,000 software verification tests before releasing each new product. ASME NQA-1-2012 (Subpart 2.7 is specific to software) is the industry- and NRC-accepted approach (consensus standard) for meeting 10 CFR Part 50, Appendix B, requirements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Light Water Sustainability Program: Optimizing Information Automation Using a New Method Based on System-Theoretic Process Analysis

This report describes the interim progress for research supporting the design and optimization of information automation systems for nuclear power plants. Much of the domestic nuclear fleet is currently focused on modernizing technologies and processes, including transitioning toward digitalization in the control room and elsewhere throughout the plant, along with a greater use of automation, artificial intelligence, robotics, and other emerging technologies. While there are significant opportunities to apply these technologies toward greater plant safety, efficiency, and overall cost-effectiveness, optimizing their design and avoiding potential safety and performance risks depends on ensuring that human-performance-related organizational and technical design issues are identified and addressed. This report describes modeling tools and techniques, based on sociotechnical system theory, to support these design goals and their application in the current research effort. The report is intended for senior nuclear energy stakeholders, including regulators, corporate management, and senior plant management. We have developed and employed a method to design an optimized information automation ecosystem (IAE) based on the systems-theoretic constructs underlying sociotechnical systems theory in general and the Systems-Theoretic Accident Modeling and Processes (STAMP) approach in particular. We argue that an IAE can be modeled as an interactive information control system whose behavior can be understood in terms of dynamic control and feedback relationships amongst the system’s technical and organizational components. Up to this point, we have employed a Causal Analysis based on STAMP (CAST) technique to examine a performance- and safety-related incident at an industry partner’s plant that involved the unintentional activation of an emergency diesel generator. This analysis provided insight into the behavior of the plant’s current information control structure within the context of a specific, significant event. Our ongoing analysis is focused on identifying near-term process improvements and longer-term design requirements for an optimized IAE system. The latter analyses will employ a second STAMP-derived technique, System-Theoretic Process Analysis (STPA). STPA is a useful modeling tool for generating and analyzing actual or potential information control structures. Finally, we have begun modeling plantwide organizational relationships and processes. Organizational system modeling will supplement our CAST and STPA findings and provide a basis for mapping out a plantwide information control architecture. CAST analysis findings indicate an important underlying contributor to the incident under investigation, and a significant risk to information automation system performance, was perceived schedule pressure, which exposed weaknesses in interdepartmental coordination between and within responsible plant organizations and challenged the resilience of established plant processes, until a human caused the initiating event. These findings are discussed in terms of their risk to overall system performance and their implications for information automation system resilience and brittleness. We present two preliminary information automation models. The proactive issue resolution model is a test case of an information automation concept with significant near-term potential for application and subsequent reduction in significant plant events. The IAE model is a more general representation of a broader, plantwide information automation system. From our results, we have generated a set of preliminary system-level requirements and safety constraints. These requirements will be further developed over the remainder of our project in collaboration with nuclear industry subject matter experts and specialists in the technical systems under consideration. Additionally, we will continue to pursue the system analyses initiated in the first part of our effort, with a particular emphasis on STPA as the main tool to identify weak or weakening control structures that affect the resilience of organizations and programs. Our intent is to broaden the scope of the analysis from an individual use case to a related set of use cases (e.g., maintenance tasks, compliance tasks) with similar human-system performance challenges. This will enable more generalized findings to refine the Proactive Issue Resolution and IAE models, as well as their system-level requirements and safety constraints. We will use organizational system modeling analyses to supplement STPA findings and model development. We conclude the report with a set of summary recommendations and an initial draft list of system-level requirements and safety constraints for optimized information automation systems.

99 GENERAL AND MISCELLANEOUS↗

Fundamental Studies of the Vibrational, Electronic, and Photophysical Properties of Tetrapyrrolic Architectures

The ability to capture and utilize light in the near-ultraviolet (NUV), visible and near-infrared (NIR-I and NIR-II) spectral regions (i.e., 320–400, 400–700, 700–1000, 1000–1700 nm) is essential for any solar-energy conversion scheme. Nature employs chlorophylls and bacteriochlorophylls in light-harvesting architectures to absorb light in the blue and red/NIR regions. Accessory pigments (carotenoids, bilins) augment absorption of the (bacterio)chlorophylls in the green region. The harvested energy is funneled to a reaction center protein, where charge separation occurs. Subsequent migration of the electron and the hole stabilizes and stores the energy from light via redox chemistry. The long-term objective of the Bocian/Holten&Kirmaier/Lindsey research program under this DOE grant has been to develop tetrapyrrole-based molecular architectures that absorb sunlight, funnel energy and separate charge with high efficiency. Integral to the program has been iterative cycles of design, synthesis and characterization that provided deep insights into the relationships between chemical composition, electronic structure, and key static and dynamic properties (vibrational, redox, photophysical, energy/charge transfer) of tetrapyrrolic systems. Such architectures included monomers, dyads, larger arrays, and complexes with accessory components. The objective was to develop molecular designs and guiding principles to enhance current and future energy-conversion schemes. Molecular arrays targeted to address one or more fundamental questions concerning light harvesting and energy/charge transfer were constructed from analogues of the naturally occurring hemes, chlorophylls and bacteriochlorophylls. Diverse, tunable synthetic building blocks were prepared that spanned the three respective tetrapyrrole families, which are the porphyrins, chlorins and bacteriochlorins. Thus, the research focused on porphyrins as well as synthetic surrogates for chlorophylls (chlorins, 13 1 -oxophorbines and chlorin-imides) and bacteriochlorophylls (bacteriochlorins, bacterio-13 1 -oxophorbines and bacteriochlorin-imides), generically termed hydroporphyrins. Although the three tetrapyrrole classes (porphyrins, chlorins and bacteriochlorins) absorb light strongly in the violet-blue spectral region, the long-wavelength absorption band typically lies in the green-orange, red, and NIR regions, respectively, with increasing intensity. Understanding the spectra, electronic structure, and energy/charge-transfer properties of such tetrapyrrolic macrocycles is of central importance for the rational design of molecular architectures for solar-energy conversion. Our integrated program of molecular design and synthesis coupled with a variety of spectroscopic, electrochemical, and computational studies have probed from first principles how structural and electronic properties of tetrapyrrolic macrocycles dictate spectral properties as well as the rates of ground-state hole/electron transfer and excited-state energy flow in multicomponent architectures. Individual molecules and multicomponent architectures were designed to test ideas of fundamental importance, often requiring the development of new synthetic methodology. The members of the collaborative team had almost daily discussions by phone and/or e-mail concerning design of molecules, flow of compounds between the labs, planning of physical characterization studies, discussing results and analysis and integrating into design of next generation architectures, and the preparation of manuscripts. Furthermore, students and postdocs in the different labs routinely communicated with one another to facilitate the advancement of the research activities. In short, a highly integrated and collaborative research program was well established among the groups. The research effort involved molecular design and synthesis of synthetic molecular architectures by the Lindsey group integrated with physicochemical and photophysical characterization by the Bocian group and the Holten&Kirmaier group (Figure 2). The Bocian group carried out electrochemical, electron paramagnetic resonance (EPR), resonance Raman (RR), and Fourier-transform infrared (FT-IR) studies, as well as density functional theory (DFT) calculations and the time-dependent extension (TDDFT) to gain insight into excited-state properties. The Holten&Kirmaier group carried out static and time-resolved absorption and fluorescence spectroscopy studies and simulated absorption spectra using molecular orbital (MO) energies from DFT as input to the four-orbital model to complement the TDDFT calculations. The combined measurements provided understanding of the vibrational/electronic properties of the individual molecules and the changes that occur upon incorporation into multicomponent architectures. This information underpinned elucidating the mechanisms and timescales of ground-state hole/electron transfer and excited-state energy and charge transfer.

14 SOLAR ENERGY↗

Portable Vibro-acoustic Testing System for In Situ Microstructure Characterization and Metrology

Portable Vibro-acoustic Testing System for In situ Microstructure Characterization and Metrology James A. Smith, Corrie I. Nichol Idaho National Laboratory, Idaho Falls ID Mostafa Fatemi, Mayo Clinic, Rochester MN There is a need in research reactors like the INL to inspect irradiated materials and structures. The goal of this work is to develop a portable scanning infrastructure for a material characterization technique called vibro-acoustography (VA) that has been developed by INL [1] for nuclear applications to characterize fuel, cladding materials, and structures. The proposed VA technology is based on ultrasound and acoustic waves; however, it provides information beyond what is available from the traditional ultrasound techniques and can expand the knowledge on nuclear material characterization and microstructure evolution. VA is a three-dimensional (3D) imaging modality based on ultrasound-stimulated acoustic emission. VA uses the force caused by the beating of two frequencies to generate an acoustic emission signal. Due to absorption or reflection, the energy density in the object at an acoustic focal point changes to produce a “radiation force.” This force locally vibrates the object, which results in an acoustic field that depends on the characteristics of the object at that point. This acoustic field is detected for every point in the object; the resulting data is used to make an image of the object’s mechanical properties. This technology is fundamentally different from traditional ultrasound imaging. More specifically, VA is: (a) based on the “radiation force,” a nonlinear phenomenon in sound propagation, and (b) sensitive to both high and low frequency dynamic responses of the object. The combination of these two response types provides information about characteristics of the object beyond what is possible by traditional ultrasound imaging. High spatial resolution is achieved, since the radiation force is localized to the focal point of the ultrasound beam. Low-frequency information, which is a representative of structural characteristics, is gained through the radiation force mechanism at a low (beat) frequency. This paper will report on the development of a portable scanning that can be set up to characterize materials and component in open water reactors and canals in situ. We will show some initial laboratory results of images generated by vibro-acoustics of surrogate fuel plates and graphite structures and discus the design of the portable system. 1. Smith, J. A., and Kunerth, D. C., “Vibro-acoustic Sensors for Materials Characterization and Related Methods and Systems,” Patent Application 13226280, 09/06/2011.

36 MATERIALS SCIENCE↗

Characterization of an architecture for front-end pixel binning in an integrating pixel array detector

Optimization of an area detector involves compromises between various parameters like frame rate, read noise, dynamic range and pixel size. We have implemented and tested a novel front-end binning design in a photon-integrating hybrid pixel array detector using the MM-PAD- 2.0 pixel architecture. In this architecture, the pixels can be optionally binned in a 2 x 2 pixel configuration using a network of switches to selectively direct the output of 4 sensor pixels to a single amplifier input. Doing this allows a trade-off between frame rate and spatial resolution. Tests show that the binned pixels perform well, but with some degradation on performance as compared to an un-binned pixel. The increased parasitic input capacitance does reduce the signal collected per x-ray as well as increases the noise of the pixel. Here, the increase in noise is, however, less than the factor of 2 increase one would observe for binning in post-processing. Spatial scans across the binned pixels show that no measured signal intensity is lost at the inner binning unit boundaries. In the high flux regime, at a 2 x 2 pixel wide beam spot (FWHM) size, binned mode responds linearly up to a photon flux of 10 7 x-rays/s, and performs comparably with un-binned mode up to a photon flux of 10 8 x-rays/s. While this study demonstrates a proof of concept for front-end binning in integrating detectors, we also identify changes to this early-stage prototype which can further improve the performance of binning pixel structures.

47 OTHER INSTRUMENTATION↗

Using ARM Measurements to Understand and Reduce the Double ITCZ Biases in the Community Atmospheric Model

1. Understanding of the observed variability of ITCZ in the equatorial eastern Pacific. The annual mean precipitation in the eastern Pacific has a maximum zonal band north of the equator in the ITCZ where the maximum SST is located. During the boreal spring (referring to February, March, and April throughout the present paper), because of the accumulated solar radiation heating and oceanic heat transport, a secondary maximum of SST exists in the southeastern equatorial Pacific. Associated with this warm SST is also a seasonal transitional maximum of precipitation in the same region in boreal spring, exhibited as a weak double ITCZ pattern in the equatorial eastern Pacific. This climatological seasonal variation, however, varies greatly from year to year: double ITCZ in the boreal spring occurs in some years but not in other years; when there a single ITCZ, it can appear either north, south or at the equator. Understanding this observed variability is critical to find the ultimate cause of the double ITCZ in climate models. Seasonal variation of ITCZ south of the eastern equatorial Pacific: By analyzing data from satellites, field measurements and atmospheric reanalysis, we have found that in the region where spurious ITCZ in models occurs, there is a “seasonal cloud transition” — from stratocumulus to shallow cumulus and eventually to deep convection —in the South Equatorial Pacific (SEP) from September to April that is similar to the spatial cloud transition from the California coast to the equator. This seasonal transition is associated with increasing sea surface temperature (SST), decreasing lower tropospheric stability and large-scale subsidence. This finding of seasonal cloud transition points to the same source of model errors in the ITCZ simulations as in simulation of stratocumulus-cumulus-deep convection transition. It provides a test for climate models to simulate the relationships between clouds and large-scale atmospheric fields in a region that features a spurious double Inter-tropical Convergence Zone (ITCZ) in most models. This work is recently published in Yu et al. (2016). Interannual variation of ITCZ south of the eastern equatorial Pacific: By analyzing data from satellites, field measurements and atmospheric reanalysis, we have characterized the interannual variation of boreal spring precipitation in the eastern tropical Pacific and found the cause of the observed interannual variability. We have shown that ITCZ in this region can occur as a single ITCZ along the Equator, single ITCZ north of the Equator, single ITCZ south of the Equator, and double ITCZ on both sides of the Equator. We have found that convective instability only plays a secondary role in the ITCZ interannual variability. Instead, the remote impact of the Pacific basin-wide SST on the horizontal gradient of surface pressure and wind convergence is the primary driver of this interannual variability. Results point to the need to include moisture convergence in convection schemes to improve the simulation of precipitation in the eastern tropical Pacific. This result has been recently submitted for publication (Yu and Zhang 2016). 2. Improvement of model parameterizations to reduce the double ITCZ bias We analyzed the current status of climate model performance in simulating precipitation in the equatorial Pacific. We have found that the double ITCZ bias has not been reduced in CMIP5 models relative to CMIP4 models. We have characterized the dynamic structure of the common bias by using precipitation, sea surface temperature, surface winds and sea-level. Results are published in Zhang et al. (2015): Since cumulus convection plays a significant role in the double ITCZ behavior in models, we have used measurements from ARM and other sources to carry out a systematic analysis of the roles of shallow and deep convection in the CAM. We found that in both CAM4 and CAM5, when the intensity of deep convection decreases as a result of parameterization change, the intensity of shallow convection increases, leading to very different changes in precipitation partitions but little change in the total precipitation. The different precipitation partition however can manifest themselves in other measures of model performances including temperature and humidity. This study points to the need to treat model physical parameterizations as integrated system rather than individual components. Results from this study are published in Wang and Zhang (2013). Since shallow convection interacts with the deep convection scheme and surface turbulence to trigger the double ITCZ, we studied methods to improve the shallow convection scheme in climate models. We investigated the bulk budgets of the vertical velocity and its parameterization in convective cores, convective updrafts, and clouds by using large-eddy simulation (LES) of four shallow convection cases including one from ARM. We proposed optimal forms of the Simpson and Wiggert equation to calculate the vertical velocity in bulk mass flux convection schemes for convective cores, convective updrafts, and convective clouds as parameterization schemes. The new scheme is published in Wang and Zhang (2014). By using long-term radar-based ground measurements from ARM, we derived a scale-aware inhomogeneity parameterization of cloud liquid water in climate models. We found a relationship between the inhomogeneity parameter and the model grid size as well as atmospheric stability. This relationship is implemented in the CESM to describe the subgrid-scale cloud inhomogeneity. Relative to the default CESM with the finite-volume dynamic core at 2-degree resolution, the new parameterization leads to smaller cloud inhomogeneity and larger cloud liquid-water path in high latitudes, and the opposite effect in low latitudes, with the regional impact on shortwave cloud radiation effect of up to 10 W/m 2 . This is due to both the smaller (larger) grid size in high (low) latitudes in the longitude-latitude grid setting of CESM and the more stable (unstable) atmosphere. This parameterization is expected lead to more realistic simulation of tropical precipitation in high-resolution models. Results from this study are reported in Xie and Zhang (2015).

54 ENVIRONMENTAL SCIENCES↗

Localizing Holes as Polarons and Predicting Band Gaps, Defect Levels, and Delithiation Energies of Solid-State Materials with a Local Exchange-Correlation Functional

Here, we assess the performance of seven exchange–correlation functionals (some with and some without a Hubbard U correction) for their ability (i) to predict band gaps of silicon, diamond, and Li-ion battery cathode materials, (ii) to localize hole polarons and predict delithiation energies in Li-ion battery cathode materials, and (iii) to predict transition levels of charge carriers of doped silicon and diamond. Both local and hybrid exchange–correlation functionals were tested. The local functionals tend to underestimate band gaps and delocalize polarons. The hybrid functionals very often give a good description of both properties, but they may not be practical for calculations involving large unit cells, large ensembles, or dynamics, and therefore a local functional with a Hubbard U correction is often used (giving the method called DFT+U), where the value of a parameter U is adjusted according to the system and the property being investigated. Keeping in mind the importance of computational cost and the undesirability of having to adjust an empirical parameter for each system or property of interest, we recently developed a local functional, namely HLE17, to try to accurately predict band gaps and excitation energies, and we validated it using mostly main-group solids and molecules. Here we test the performance of HLE17 for its ability to predict band gaps and localize polarons in other solid-state materials, and we compare its performance to that of popular local functionals (PBE, PBEsol, and TPSS), a range-separated hybrid functional with screened exchange (HSE06), and DFT+U. We find that HLE17 predicts more accurate band gaps than other local functionals and can localize holes as polarons, which other local functionals usually fail to do, and for a number of cases it is comparable in performance quality to Hubbard-corrected functionals without the need for system-specific parametrization and to hybrid functionals without the high cost. Because HLE17 does not predict accurate lattice constants, we use the single-point method of quantum chemistry, where the geometry is optimized with one functional and the band gap is calculated with HLE17, or we perform calculations with the lattice constants obtained by TPSS and both the fractional intracell coordinates and the electronic structure obtained by HLE17 (a new method denoted HLE17\\TPSS). In particular, we performed calculations by HLE17//TPSS, HLE17//HSE06, HLE17//DFT+U, and HLE17\\TPSS, and these methods usually agree well with each other and give values similar to experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mid-Latitude Circulation and Extremes in a Changing Climate. Final Technical Report for DOE Award DE-SC0012599

Mid-latitude extreme weather events are responsible for a large part of climate-related damage. Yet large uncertainties remain in climate model projections of heat waves, droughts, and heavy rain/snow events on regional scales, limiting our ability to effectively use these projections for climate adaptation and mitigation. These uncertainties can be attributed to both the lack of spatial resolution in the models, and to the lack of a dynamical understanding of these extremes. The objectives of this proposal were to first to develop empirical and theoretical understanding of how fine-scale features of the atmospheric circulation that are responsible for extreme events (such as blocking, atmospheric rivers, Rossby wave breaking) are related to more easily quantifiable larger scale circulation patterns, and on the basis of this develop specific metrics that can be applied to a range of atmospheric and coupled atmosphere-ocean models to diagnose their prediction of extremes. These metrics were used to understand what model resolution is necessary to describe the probability of extreme events realistically, and will be applied to CMIP5 models and high resolution time-slice projections of future climate to better understand the potential for the shifts and intensification of extremes over North America in a changing climate. The approach was to relate the fine-scale features to the large scales in current climate simulations, seasonal re-forecasts, and climate change projections in a very wide range of models, including the atmospheric and coupled models of ECMWF over a range of horizontal resolutions (125 to 10 km), aqua-planet configuration of the Model for Prediction Across Scales and High Order Method Modeling Environments (resolutions ranging from 240 km – 7.5 km) with various physics suites, and selected CMIP5 model simulations. The large scale circulation will be quantified both on the basis of the well-tested preferred circulation regime approach, and very recently developed measures, the finite amplitude Wave Activity and its spectrum. The fine scale structures related to extremes were diagnosed following the latest approaches in the literature. The goal is to use the large scale measures as indicators of the probability of occurrence of the finer scale structures, and hence extreme events. These indicators will then be applied to the CMIP5 models and time-slice projections of a future climate. The expected outcomes of this research were: (i) a unified methodology and set of metrics that will enable modelers to assess the probability of extreme events in various current and future climate simulations, (ii) an estimate of the model resolution needed to capture extreme events in a realistic way, and (iii) a better understanding of the extreme events over North America in a changing climate.

54 ENVIRONMENTAL SCIENCES↗

Mid-Latitude Circulation and Extremes in a Changing Climate

Mid-latitude extreme weather events are responsible for a large part of climate-related damage. Yet large uncertainties remain in climate model projections of heat waves, droughts, and heavy rain/snow events on regional scales, limiting our ability to effectively use these projections for climate adaptation and mitigation. These uncertainties can be attributed to both the lack of spatial resolution in the models, and to the lack of a dynamical understanding of these extremes. The objectives of this proposal were to first to develop empirical and theoretical understanding of how fine-scale features of the atmospheric circulation that are responsible for extreme events (such as blocking, atmospheric rivers, Rossby wave breaking) are related to more easily quantifiable larger scale circulation patterns, and on the basis of this develop specific metrics that can be applied to a range of atmospheric and coupled atmosphere-ocean models to diagnose their prediction of extremes. These metrics were used to understand what model resolution is necessary to describe the probability of extreme events realistically, and will be applied to CMIP5 models and high resolution time-slice projections of future climate to better understand the potential for the shifts and intensification of extremes over North America in a changing climate. The approach was to relate the fine-scale features to the large scales in current climate simulations, seasonal re-forecasts, and climate change projections in a very wide range of models, including the atmospheric and coupled models of ECMWF over a range of horizontal resolutions (125 to 10 km), aqua-planet configuration of the Model for Prediction Across Scales and High Order Method Modeling Environments (resolutions ranging from 240 km – 7.5 km) with various physics suites, and selected CMIP5 model simulations. The large scale circulation will be quantified both on the basis of the well-tested preferred circulation regime approach, and very recently developed measures, the finite amplitude Wave Activity and its spectrum. The fine scale structures related to extremes were diagnosed following the latest approaches in the literature. The goal is to use the large scale measures as indicators of the probability of occurrence of the finer scale structures, and hence extreme events. These indicators will then be applied to the CMIP5 models and time-slice projections of a future climate. The expected outcomes of this research were: (i) a unified methodology and set of metrics that will enable modelers to assess the probability of extreme events in various current and future climate simulations, (ii) an estimate of the model resolution needed to capture extreme events in a realistic way, and (iii) a better understanding of the extreme events over North America in a changing climate.

54 ENVIRONMENTAL SCIENCES↗

High burn-up spent nuclear fuel vibration integrity study - 15134

The Oak Ridge National Laboratory (ORNL) has developed the cyclic integrated reversible-bending fatigue tester (CIRFT) approach to successfully demonstrate the controllable fatigue fracture of high burnup (HBU) spent nuclear fuel (SNF) in a normal vibration mode. CIRFT enables examination of the underlying mechanisms of SNF system dynamic performance. Due to the inhomogeneous composite structure of the SNF system, the detailed mechanisms of pellet-pellet and pellet-clad interactions and the stress concentration effects at the pellet-pellet interface cannot be readily obtained from a CIRFT system measurement. Therefore, finite element analyses (FEAs) are being used to complement the experimental work by estimating values of local stresses and strains from the measured global moment-curvature. These results and their implications will be used for further investigation. Major findings of CIRFT on the HBU SNF are as follows: - The degree of bonding between pellets, as well as the bonding between pellets and clad in the SNF system plays an important role in the vibration performance of SNF. - Fuel structure contributes to SNF system stiffness. - There are significant variations in stress and curvature of SNF systems during vibration cycles resulting from segment pellet and clad interactions. - SNF failure initiates at the pellet-pellet interface region and appears to be spontaneous. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantification of spatial and seasonal trends in the atmosphere and construction of statistical models for infrasonic propagation

SUMMARY Infrasonic waves are influenced by variations in the density, pressure and temperature as well as the ambient winds. Modelling infrasonic propagation can be challenging due to the dynamic nature of the atmosphere as well as the sparseness of measurements which result in variability and notable uncertainty. A framework is presented to quantify spatial and seasonal trends in atmospheric structure via analysis of the effective sound speed profile and identification of temporal trends in the middle atmospheric waveguide produced by the circumpolar vortex winds. Seasonal definitions identifying typical atmospheric structures during the summer, winter and spring/fall transition periods are identified using atmospheric data from 2010 through 2020. Seasonal trend analysis is conducted for a number of locations across the contiguous United States to quantify spatial variations in atmospheric structure that impact infrasonic propagation. For each season and location, empirical orthogonal function analysis is used to reduce the historical archive of atmospheric data into a smaller representative set that can be analysed using numerical tools more efficiently. Infrasonic ray tracing and finite-frequency modal propagation analyses are applied to construct propagation path geometry and transmission loss statistics which are useful in localization and yield estimation for infrasonic sources, respectively. An example application is detailed in which transmission loss statistics are combined with an explosive source model and noise statistics to quantify the capability of a network to detect nearby sources.

58 GEOSCIENCES↗

Regional climate change: consensus, discrepancies, and ways forward

Climate change has emerged across many regions. Some observed regional climate changes, such as amplified Arctic warming and land-sea warming contrasts have been predicted by climate models. However, many other observed regional changes, such as changes in tropical sea surface temperature and monsoon rainfall are not well simulated by climate model ensembles even when taking into account natural internal variability and structural uncertainties in the response of models to anthropogenic radiative forcing. This suggests climate model predictions may not fully reflect what our future will look like. The discrepancies between models and observations are not well understood due to several real and apparent puzzles and limitations such as the “signal-to-noise paradox” and real-world record-shattering extremes falling outside of the possible range predicted by models. Addressing these discrepancies, puzzles and limitations is essential, because understanding and reliably predicting regional climate change is necessary in order to communicate effectively about the underlying drivers of change, provide reliable information to stakeholders, enable societies to adapt, and increase resilience and reduce vulnerability. The challenges of achieving this are greater in the Global South, especially because of the lack of observational data over long time periods and a lack of scientific focus on Global South climate change. To address discrepancies between observations and models, it is important to prioritize resources for understanding regional climate predictions and analyzing where and why models and observations disagree via testing hypotheses of drivers of biases using observations and models. Gaps in understanding can be discovered and filled by exploiting new tools, such as artificial intelligence/machine learning, high-resolution models, new modeling experiments in the model hierarchy, better quantification of forcing, and new observations. Conscious efforts are needed toward creating opportunities that allow regional experts, particularly those from the Global South, to take the lead in regional climate research. This includes co-learning in technical aspects of analyzing simulations and in the physics and dynamics of regional climate change. Finally, improved methods of regional climate communication are needed, which account for the underlying uncertainties, in order to provide reliable and actionable information to stakeholders and the media.

54 ENVIRONMENTAL SCIENCES↗

The Aerogel Cerenkov detector for the SHMS magnetic spectrometer in Hall C at Jefferson Lab

Hadronic reactions producing strange quarks such as exclusive or semi-inclusive kaon production, play an important role in studies of hadron structure and the dynamics that bind the most basic elements of nuclear physics. The small-angle capability of the new Super High Momentum Spectrometer (SHMS) in Hall C, coupled with its high momentum reach - up to the anticipated 11-GeV beam energy in Hall C - and coincidence capability with the well-understood High Momentum Spectrometer, will allow for probes of such hadron structure involving strangeness down to the smallest distance scales to date. To cleanly select the kaons, a threshold aerogel Cerenkov detector has been constructed for the SHMS. The detector consists of an aerogel tray followed by a diffusion box. Four trays for aerogel of nominal refractive indices of n =1.030, 1.020, 1.015 and 1.011 were constructed. The tray combination will allow for identification of kaons from 1 GeV/c up to 7.2 GeV/c, reaching 10 -2 proton and 10 -3 pion rejection, with kaon detection efficiency better than 95%. The diffusion box of the detector is equipped with 14 five-inch diameter photomultiplier tubes. Its interior walls are covered with Gore diffusive reflector, which is superior to the commonly used Millipore paper and improved the detector performance by 35%. The inner surface of the two aerogel trays with higher refractive index is covered with Millipore paper, however, those two trays with lower aerogel refractive index are again covered with Gore diffusive reflector for higher performance. The measured mean number of photoelectrons in saturation is ~12 for n =1.030, ~8 for n =1.020, ~10 for n =1.015, and ~5.5 for n =1.011. Here the design details, the results of component characterization, and initial performance tests and optimization of the detector are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗