Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesoscale science”

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 721 records · Page 40

The Thunderstorms With the Greatest Lightning Densities on Earth

The most intense thunderstorms on Earth were surveyed using the comprehensive meteorological instrumentation on the Tropical Rainfall Measuring Mission (TRMM) satellite. Expansive land-based Mesoscale Convective Systems (MCSs) were consistently identified among the Earth's most intense thunderstorms, with their organization into many convective cells spanning a large areal extent permitting exceptional overall flash rates for these storms. In this study, we identify a new class of extreme thunderstorm. Lightning-dense thunderstorms are relatively compact convective storms whose concentrated lightning activity hinders our ability to accurately measure their flash rates. The top storms have a flash rate of one flash spanning many seconds, as there is insufficient separation to distinguish one flash from another. While any particularly active convective cell could be capable of producing high lightning densities, we find that thunderstorms with the greatest lightning densities on Earth are found in maritime thunderstorms that have not been appreciated in prior work due to the inaccurate flash rate measurements. These storms that are mostly found throughout the Gulf of Mexico and east of South Africa (among other coastal and oceanic regions) have measured TRMM proxies for convective intensity that rival the top MCS thunderstorms, but their horizontal and vertical dimensions are small by comparison. Thus, the necessary microphysical elements for electrification processes are more highly concentrated, enabling the observed extreme lightning densities.

54 ENVIRONMENTAL SCIENCES↗

The influence of multiple groups of biological ice nucleating particles on microphysical properties of mixed-phase clouds observed during MC3E

Abstract. A new empirical parameterization (EP) for multiple groups of primary biological aerosol particles (PBAPs) is implemented in the aerosol–cloud model (AC) to investigate their roles as ice nucleating particles (INPs). The EP describes the heterogeneous ice nucleation by (1) fungal spores, (2) bacteria, (3) pollen, (4) detritus of plants, animals, and viruses, and (5) algae. Each group includes fragments from the originally emitted particles. A high-resolution simulation of a midlatitude mesoscale squall line by AC is validated against airborne and ground observations. Sensitivity tests are carried out by varying the initial vertical profiles of the loadings of individual PBAP groups. The resulting changes in warm and ice cloud microphysical parameters are investigated. The changes in warm microphysical parameters, including liquid water content and cloud droplet number concentration, are minimal (<10 %). Overall, PBAPs have little effect on the ice number concentration (<6 %) in the convective region. In the stratiform region, increasing the initial PBAP loadings by a factor of 1000 resulted in less than 40 % change in ice number concentrations. The total ice concentration is mostly controlled by various mechanisms of secondary ice production (SIP). However, when SIP is intentionally shut down in sensitivity tests, increasing the PBAP loading by a factor of 100 has an effect of less than 3 % on the ice phase. Further sensitivity tests revealed that PBAPs have little effect on surface precipitation and on the shortwave and longwave flux (<4 %) for a 100-fold perturbation in PBAPs.

54 ENVIRONMENTAL SCIENCES↗

Sensitivity of mesoscale modeling to urban morphological feature inputs and implications for characterizing urban sustainability

We examine the differences in meteorological output from the Weather Research and Forecasting (WRF) model run at 270 m horizontal resolution using 10 m, 100 m and 1 km resolution 3D neighborhood morphological inputs and with no morphological inputs. We find that the spatial variability in temperature, humidity, and other meteorological variables across the city can vary with the resolution and the coverage of the 3D urban morphological input, and that larger differences occur between simulations run without 3D morphological input and those run with some type of 3D morphology. We also find that the inclusion of input-building-defined roughness length calculations would improve simulation results further. We show that these inputs produce different patterns of heat wave spatial heterogeneity across the city of Washington, DC. These findings suggest that understanding neighborhood level urban sustainability under extreme heat waves, especially for vulnerable neighborhoods, requires attention to the representation of surface terrain in numerical weather models.

54 ENVIRONMENTAL SCIENCES↗

The Hazards Posed by Mesoscale Lightning Megaflashes

Lightning megaflashes extending over >100-km distances have been observed by the Geostationary Lightning Mappers (GLMs) on NOAA’s R-series Geostationary Operational Environmental Satellites (GOES). The hazards posed by megaflashes are unclear, however, because of limitations in the GLM data. We address these by reprocessing GOES-16 GLM measurements from 1 January 2018 to 15 January 2020 and integrating them with Earth Networks Global Lightning Network (ENGLN) observations. ENGLN verified 194 880 GLM megaflashes as natural lightning. Of these, 127 479 flashes occurred following the October 2018 GLM software update that standardized GLM timing. Reprocessed GLM/ENGLN lightning maps from these postupdate cases provide a comprehensive view of how individual megaflashes evolve. This megaflash dataset is used to generate statistics that describe their hazards. The average megaflash produces 5–7 cloud-to-ground (CG) strokes that are spread across 40%–50% of the flash extent. As flash extent increases beyond 100 km, megaflashes become concentrated in key hot-spot regions in North and South America while the number of CG and intracloud events per flash and the overall peak current increase. CGs in the larger megaflashes occur over 80% of the flash extent measured by GLM, and the majority contain regions where the megaflash is the only lightning activity in the preceding hour. These statistics demonstrate that there is no safe location below an electrified cloud that is producing megaflashes, and current lightning safety guidance is not always sufficient to mitigate megaflash hazards.

54 ENVIRONMENTAL SCIENCES↗

Investigation of δ zirconium hydride morphology in a single crystal using quantitative phase field simulations supported by experiments

In light water nuclear reactors, waterside corrosion of the cladding material leads to the production of hydrogen, a fraction of which is picked up by the zirconium cladding and precipitates into brittle hydride particles. These nanoscale hydride particles aggregate into mesoscale hydride clusters. The principal stacking direction of the nanoscale hydrides precipitated in the cladding tube changes from circumferential in the absence of applied stress to radial under circumferential applied stress. A quantitative phase field model has been developed to predict the hydride morphology observed experimentally and identify the mechanisms responsible for nanoscale hydride stacking. The model focuses on nanoscale hydride precipitation in a single zirconium grain with a detailed description of the anisotropic elastic contribution. The model predictions concerning the shape, orientation, and stacking behavior of nanoscale hydride are analyzed and compared with experimental observations. The model accurately accounts for the experimentally observed elongated nanoscale hydride shape and the stacking of hydrides along the basal plane of the hexagonal zirconium matrix. When investigating the role of applied stress in hydride morphology, the model challenges some of the mechanisms previously proposed to explain hydride reorientation. Finally, although hydride reorientation has been hypothesized to be caused by a change in nanoscale hydride shape, the current study shows that these mechanisms are unlikely to occur.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Physics and chemistry from parsimonious representations: image analysis via invariant variational autoencoders

Electron, optical, and scanning probe microscopy methods are generating ever increasing volume of image data containing information on atomic and mesoscale structures and functionalities. This necessitates the development of the machine learning methods for discovery of physical and chemical phenomena from the data, such as manifestations of symmetry breaking phenomena in electron and scanning tunneling microscopy images, or variability of the nanoparticles. Variational autoencoders (VAEs) are emerging as a powerful paradigm for the unsupervised data analysis, allowing to disentangle the factors of variability and discover optimal parsimonious representation. Here, we summarize recent developments in VAEs, covering the basic principles and intuition behind the VAEs. The invariant VAEs are introduced as an approach to accommodate scale and translation invariances present in imaging data and separate known factors of variations from the ones to be discovered. We further describe the opportunities enabled by the control over VAE architecture, including conditional, semi-supervised, and joint VAEs. Several case studies of VAE applications for toy models and experimental datasets in Scanning Transmission Electron Microscopy are discussed, emphasizing the deep connection between VAE and basic physical principles. Python codes and datasets discussed in this article are available at https://github.com/saimani5/VAE-tutorials and can be used by researchers as an application guide when applying these to their own datasets.

36 MATERIALS SCIENCE↗

Quantifying Uncertainty in Ice Particle Velocity–Dimension Relationships Using MC3E Observations

Abstract Ice particle terminal fall velocity ( V t ) is fundamental for determining microphysical processes, yet remains extremely challenging to measure. Current theoretical best estimates of V t are functions of Reynolds number. The Reynolds number is related to the Best number, which is a function of ice particle mass, area ratio ( A r ), and maximum dimension ( D max ). These estimates are not conducive for use in most models since model parameterizations often take the form , where ( α , β ) depend on habit and D max . A previously developed framework is used to determine surfaces of equally plausible ( α , β ) coefficients whereby ice particle size/shape distributions are combined with V t best estimates to determine mass- ( V M ) or reflectivity-weighted ( V Z ) velocities that closely match parameterized V M ,SD or V Z ,SD calculated using the ( α , β ) coefficients using two approaches. The first uses surfaces of equally plausible ( a , b ) coefficients describing mass ( M )–dimension relationships (i.e., ) to calculate mass- or reflectivity-weighted velocity from size/shape distributions that are then used to determine ( α , β ) coefficients. The second investigates how uncertainties in A r , D max , and size distribution N ( D ) affect V M or V Z . For seven of nine flight legs flown on 20 and 23 May 2011 during the Mesoscale Continental Convective Clouds Experiment (MC3E), uncertainty from natural parameter variability—namely, the variability in ice particle parameters in similar meteorological conditions—exceeds uncertainties arising from different A r assumptions or D max estimates. The combined uncertainty between A r , D max , and N ( D ) produced smaller variability in ( α , β ) compared to varying M ( D ), demonstrating M ( D ) must be accurately quantified for model fall velocities. Primary sources of uncertainty vary considerably depending on environmental conditions. Significance Statement Ice particle fall velocity is fundamental for numerous processes within clouds, and hence is a critical property that must be accurately represented in weather and climate models. Using aircraft observations of ice particle shapes and sizes obtained in clouds behind midlatitude thunderstorms, this work develops a new framework for estimating ice particle fall velocities and their uncertainty, including quantifying the importance of different uncertainty sources from cloud microphysics measurements. Natural parameter variability contributes the most uncertainty in ice particle fall velocity estimates, although other sources can also be important contributors to uncertainty in certain conditions. Additional work examining ice particle data is needed to further understand how dependent uncertainty in certain ice particle properties are to local environmental conditions.

Meteorology & Atmospheric Sciences↗

The fully coupled regionally refined model of E3SM version 2: overview of the atmosphere, land, and river results

Abstract. This paper provides an overview of the United States (US) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2 (E3SMv2) fully coupled regionally refined model (RRM) and documents the overall atmosphere, land, and river results from the Coupled Model Intercomparison Project 6 (CMIP6) DECK (Diagnosis, Evaluation, and Characterization of Klima) and historical simulations – a first-of-its-kind set of climate production simulations using RRM. The North American (NA) RRM (NARRM) is developed as the high-resolution configuration of E3SMv2 with the primary goal of more explicitly addressing DOE's mission needs regarding impacts to the US energy sector facing Earth system changes. The NARRM features finer horizontal resolution grids centered over NA, consisting of 25→100 km atmosphere and land, a 0.125∘ river-routing model, and 14→60 km ocean and sea ice. By design, the computational cost of NARRM is ∼3× of the uniform low-resolution (LR) model at 100 km but only ∼ 10 %–20 % of a globally uniform high-resolution model at 25 km. A novel hybrid time step strategy for the atmosphere is key for NARRM to achieve improved climate simulation fidelity within the high-resolution patch without sacrificing the overall global performance. The global climate, including climatology, time series, sensitivity, and feedback, is confirmed to be largely identical between NARRM and LR as quantified with typical climate metrics. Over the refined NA area, NARRM is generally superior to LR, including for precipitation and clouds over the contiguous US (CONUS), summertime marine stratocumulus clouds off the coast of California, liquid and ice phase clouds near the North Pole region, extratropical cyclones, and spatial variability in land hydrological processes. The improvements over land are related to the better-resolved topography in NARRM, whereas those over ocean are attributable to the improved air–sea interactions with finer grids for both atmosphere and ocean and sea ice. Some features appear insensitive to the resolution change analyzed here, for instance the diurnal propagation of organized mesoscale convective systems over CONUS and the warm-season land–atmosphere coupling at the southern Great Plains. In summary, our study presents a realistically efficient approach to leverage the fully coupled RRM framework for a standard Earth system model release and high-resolution climate production simulations.

54 ENVIRONMENTAL SCIENCES↗

Electronic 1/f noise as a probe of dimensional effects on spin-glass dynamics

Over the past decade, spin-glass simulations have improved to the point that they now access time- and length-scales comparable to experiments at the mesoscale. A recent series of thin-film field-cooled/zero-field-cooled magnetization (FC/ZFC) experiments demonstrated activated spin dynamics, with a temperature-independent activation energy proportional to the logarithm of the film thickness and with coefficients in remarkable agreement with the simulation. These measurements require the application of small magnetic fields, which has been shown to affect the spin-glass energy landscape. Measurements of the 1/f noise in metallic spin-glasses have been previously shown to be a sensitive probe of the spin dynamics, and the measurements can be made without applying a magnetic field. In this mini-review, we review these techniques and discuss how transport measurements can fit into the current landscape of spin-glass measurements. We compare previous measurements to more recent measurements on similar films, made with ostensibly different cooling protocols, and compare both the previous and recent measurements to the magnetometry. The transport measurements—taken over a wider range of temperature than magnetometry—suggest that the maximum spin-glass energy barrier height is temperature-dependent, not fixed, possibly due to two-dimensional dynamics. We discuss this possibility, along with future measurements, which may be able to resolve this mystery.

1/f noise↗

Quantifying Flood Frequency Associated with Clustered Mesoscale Convective Systems in the United States

Mesoscale convective systems (MCSs) that are clustered in time and space can have a broader impact on flooding because they have larger area coverage than that of individual MCSs. The goal of this study is to understand the flood likelihood associated with MCS clusters. To achieve this, floods in the Storm Events Database in April–August of 2007–17 are matched with clustered MCSs identified from a high-resolution MCS dataset and terrestrial conditions in a land surface dataset over the central-eastern United States. In this work, our analysis indicates that clustered MCSs preferentially occurring in April–June are more effective at producing floods, which also last longer due to the greater rainfall per area and wetter initial soil conditions and, hence, produce greater runoff per area than nonclustered MCSs. Similar increases of flood occurrence with cluster-total rainfall size and wetter soils are also observed for each MCS cluster, especially for the overlapping rainfall areas within each cluster. These areas receive rainfall from multiple MCSs that progressively wet the soils and are therefore associated with higher flood likelihood. This study underscores the importance to understand clustered MCSs to better understand flood risks and their future changes.

54 ENVIRONMENTAL SCIENCES↗

The remarkable inefficiency of stratocumulus

Marine stratocumulus clouds play a central role in Earth's climate system by reflecting incoming solar radiation and exerting a strong cooling effect. Their organization into open and closed mesoscale cellular morphologies can be thought of as an example of bistable dynamics driven by aerosol–cloud interactions and mesoscale processes. From the perspective of non-equilibrium thermodynamics, these structures are an example of a far-from-equilibrium open system that continuously produces and exports entropy. While entropy production has been studied in idealized deep convective systems, it has not yet been quantified for shallow clouds. Here, we compute and decompose the internal entropy production of open- and closed-cell stratocumulus using an ensemble of large-eddy simulations. We show that the overall entropy production of stratocumulus is low, reflecting the limited vertical extent and corresponding reduced ability to utilize the energy fluxes at the system's boundaries. Moist processes dominate the overall irreversibility, which, combined with their low entropy production, leads to a mechanical efficiency about an order of magnitude smaller than in deep convective systems. Although the dominant irreversible processes differ between open- and closed-cell regimes, the distributions of total entropy production largely overlap across the ensemble, limiting the ability to distinguish the dynamics of individual cases based solely on total entropy production.

54 ENVIRONMENTAL SCIENCES↗

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Abstract. Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model configurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

54 ENVIRONMENTAL SCIENCES↗

The Application of a Genetic Algorithm to the Optimization of a Mesoscale Model for Emergency Response

Besides solving the equations of momentum, heat, and moisture transport on the model grid, mesoscale weather models must account for subgrid-scale processes that affect the resolved model variables. These are simulated with model parameterizations, which often rely on values preset by the user. Such “free” model parameters, along with others set to initialize the model, are often poorly constrained, requiring that a user select each from a range of plausible values. Finding the values to optimize any forecasting tool can be accomplished with a search algorithm, and one such process—the genetic algorithm (GA)—has become especially popular. As applied to modeling, GAs represent a Darwinian process: an ensemble of simulations is run with a different set of parameter values for each member, and the members subsequently judged to be most accurate are selected as “parents” who pass their parameters onto a new generation. At the U.S. Department of Energy’s Savannah River Site in South Carolina, we are applying a GA to the Regional Atmospheric Modeling System (RAMS) mesoscale weather model, which supplies input to a model to simulate the dispersion of an airborne contaminant as part of the site’s emergency response preparations. An ensemble of forecasts is run each day, weather data are used to “score” the individual members of the ensemble, and the parameters from the best members are used for the next day’s forecasts. As meteorological conditions change, the parameters change as well, maintaining a model configuration that is best adapted to atmospheric conditions. Significance Statement We wanted to develop a forecasting system by which a weather model is run over the Savannah River Site each day and repeatedly adjusted according to how well it performed the previous day. To run the model, a series of values (parameters) must be set to control how the model will calculate winds, temperatures, and other desired variables. Each day the model was run several times using different combinations of these parameters and later compared with observed meteorological conditions. Parameters that produced the most accurate forecasts were preferentially reused to create the forecasts for the next day. The process was tested for the summer of 2020 and exhibited lower errors than forecasts produced by the model using default values of the parameters.

54 ENVIRONMENTAL SCIENCES↗

Meso- to microscale modeling of atmospheric stability effects on wind turbine wake behavior in complex terrain

Abstract. Terrain-induced flow phenomena modulate wind turbine performance and wake behavior in ways that are not adequately accounted for in typical wind turbine wake and wind plant design models. In this work, we simulate flow over two parallel ridges with a wind turbine on one of the ridges, focusing on conditions observed during the Perdigão field campaign in 2017. Two case studies are selected to be representative of typical flow conditions at the site, including the effects of atmospheric stability: a stable case where a mountain wave occurs (as in ∼ 50 % of the nights observed) and a convective case where a recirculation zone forms in the lee of the ridge with the turbine (as occurred over 50 % of the time with upstream winds normal to the ridgeline). We use the Weather Research and Forecasting Model (WRF), dynamically downscaled from the mesoscale (6.75 km resolution) to microscale large-eddy simulation (LES) at 10 m resolution, where a generalized actuator disk (GAD) wind turbine parameterization is used to simulate turbine wakes. We compare the WRF–LES–GAD model results to data from meteorological towers, lidars, and a tethered lifting system, showing good qualitative and quantitative agreement for both case studies. Significantly, the wind turbine wake shows different amounts of vertical deflection from the terrain and persistence downstream in the two stability regimes. In the stable case, the wake follows the terrain along with the mountain wave and deflects downwards by nearly 100 m below hub height at four rotor diameters downstream. In the convective case, the wake deflects above the recirculation zone over 40 m above hub height at the same downstream distance. Overall, the WRF–LES–GAD model is able to capture the observed behavior of the wind turbine wakes, demonstrating the model's ability to represent wakes over complex terrain for two distinct and representative atmospheric stability classes, and, potentially, to improve wind turbine siting and operation in hilly landscapes.

17 WIND ENERGY↗

Competing Interactions between Mesoscale Length-Scales, Order-Disorder, and Martensitic Transformation in Ferromagnetic Shape Memory Alloys

In the present study, the effects of composition and heat treatments and the resulting microstructural changes on the martensitic transformation and ferromagnetic transition have been investigated in NiCoMnIn metamagnetic shape memory alloys. In this shape memory alloy system, it is observed that upon heat treatments at a wide temperature range, the onset temperature of the martensitic transformation follows a non-monotonic behavior with respect to heat treatment time, and the nature of the non-monotonic behavior is also a function of composition. This behavior cannot be attributed to well-known factors such as precipitation, change in local composition due to the precipitation and/or global degree of order. In this work, a systematic investigation through synthesis, thermal processing and characterization via thermo-physical measurements, transmission electron microscopy and in-situ synchrotron x-ray diffraction experiments has been used to correlate the non-monotonic dependence of the martensitic transformation and ferromagnetic transition temperatures to the evolution of L2 1 domains arising from a order/disorder phase transition. A thermodynamic model for the magneto-structural transition is combined with classical nucleation theory to further ascertain the role of microstructural length-scales on the onset of the martensitic transformation. This work thus provides understanding of the thermodynamic and kinetic factors that can be controlled to tune the coupled magneto-structural transformations in NiCoMnIn metamagnetic shape memory alloys.

36 MATERIALS SCIENCE↗

Characteristics of Precipitation and Mesoscale Convective Systems Over the Peruvian Central Andes in Multi 5‐Year Convection‐Permitting Simulations

Abstract Using the Weather Research and Forecasting model with two planetary boundary layer schemes, ACM2 and MYNN, convection‐permitting model (CPM) regional climate simulations were conducted for a 6‐year period, including a one‐year spin‐up period, at a 15‐km grid spacing covering entire South America and a nested convection‐permitting 3‐km grid spacing covering the Peruvian central Andes region. These two CPM simulations along with a 4‐km simulation covering South America produced by National Center for Atmospheric Research (NCAR), three gridded precipitation products, and rain gauge data in Peru and Brazil, are used to document the characteristics of precipitation and mesoscale convective systems (MCSs) in the Peruvian central Andes region. Results show that all km‐scale simulations generally capture the spatiotemporal patterns of precipitation and MCSs at both seasonal and diurnal scales, although biases exist in aspects such as precipitation intensity and MCS frequency, size, propagation speed, and associated precipitation intensity. The 3‐km simulation using MYNN scheme generally outperforms the other simulations in capturing seasonal and diurnal precipitation over the mountain, while both it and the 4‐km simulation demonstrate superior performance in the western Amazon Basin, based on the comparison to the gridded precipitation products and gauge data. Dynamic factors, primarily low‐level jet and terrain‐induced uplift, are the key drivers for precipitation and MCS genesis along the east slope of the Andes, while thermodynamic factors control the precipitation and MCS activity in the western Amazon Basin and over elevated mountainous regions. The study suggests model improvements and better model configurations for future regional climate projections.

Meteorology & Atmospheric Sciences↗

Evaluation and Uncertainty of Radioxenon Transport with a Mesoscale Model after the February 2013 Underground Test in North Korea

Abstract The transport of radioxenon released from the February 2013 underground nuclear weapons test in North Korea was analyzed at two receptors—one at the Comprehensive Nuclear Test Ban Treaty site Rn58 in Russia (400 km downwind) and a second at Rn38 in Japan (1000 km downwind). Transport was modeled with two ensembles of mesoscale simulations, one generated with varying initial and lateral boundary conditions taken from the Global Forecasting System uncertainty ensemble, and a second created from different parameterizations and surface conditions. The wind variability was similar for the two ensembles and consistent with observations at 925 mb (1 mb = 1 hPa) but not at the surface. Biases in calculated surface winds and the radioxenon concentration in the ensembles were attributed mainly to poor simulation of the sea breeze at both locations and mountain lee affects at Rn38 in Japan. These wind regimes affected the timing of the surface radioxenon plume at Rn58 and its duration at Rn38. Surface wind variability induced by terrain and land–sea contrast (the sea breeze) also had a significant effect on the surface winds and plume dynamics, including blocking of flow approaching elevated terrain near Vladivostok and the west side of Japan. Increased plume uncertainty was seen at night because of surface wind variability. Measured surface chemical variability was larger than found in the first European Tracer Experiment in central Europe. The study found that horizontal model resolution contributes to uncertainty but not as much as vertical resolution, boundary layer parameterizations, and assimilation of surface meteorological data near the receptor.

54 ENVIRONMENTAL SCIENCES↗

Microstructural-Level Fuel Performance Modeling of U-Mo Monolithic Fuel

As the physics that governs the microstructural evolution of nuclear fuel span various time and spatial scales, to fully understand the fuel behavior inevitably involves atomic to mesoscale resolution that can be difficult to determine experimentally. Microstructural-level modeling and simulations can be used to develop physics-based materials models that can provide physical understanding to inform fabrication process control, as well as a valuable feedback mechanism between post-irradiation examination (PIE) results and fabrication parameters. In accordance with the program schedule, the primary goals of the microstructure modeling effort are to: 1. Address critical microstructural questions and provide practical guidance to the fabricator via the fuel product specification; 2. Provide mechanistic inputs for the existing fuel performance code to improve its descriptive and predictive capability at the macroscopic scale. In fiscal year (FY)-21, the work scope consisted of six main facets: (1) the effect of carbides on fuel performance; (2) gas diffusivity in different phases; (3) integration of microstructural fuel performance modeling; (4) property degradation; (5) irradiation creep; and (6) historical analysis of microstructure data. Brief summaries of each are included below.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗