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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Towards a Unified Setup to Simulate Mid‐Latitude and Tropical Mesoscale Convective Systems at Kilometer‐Scales

Abstract Mesoscale convective systems (MCSs) are the main source of precipitation in the tropics and parts of the mid‐latitudes and are responsible for high‐impact weather worldwide. Studies showed that deficiencies in simulating mid‐latitude MCSs in state‐of‐the‐art climate models can be alleviated by kilometer‐scale models. However, whether these models can also improve tropical MCSs and whether we can find model settings that perform well in both regions is understudied. We take advantage of high‐quality MCS observations collected over the Atmospheric Radiation Measurement (ARM) facilities in the US Southern Great Plains (SGP) and the Amazon basin near Manaus (MAO) to evaluate a perturbed physics ensemble of simulated MCSs with 4 km horizontal grid spacing. A new model evaluation method is developed that enables to distinguish biases stemming from spatiotemporal displacements of MCSs from biases in their reflectivity and cloud shield. Amazon MCSs are similarly well simulated across these evaluation metrics than SGP MCSs despite the challenges anticipated from weaker large‐scale forcing in the tropics. Generally, SGP MCSs are more sensitive to the choice of model microphysics, while Amazon cases are more sensitive to the planetary boundary layer (PBL) scheme. Although our tested model physics combinations had strengths and weaknesses, combinations that performed well for SGP simulations result in worse results in the Amazon basin and vice versa. However, we identified model settings that perform well at both locations, which include the Thompson and Morrison microphysics coupled with the Yonsei University (YSU) PBL scheme and the Thompson scheme coupled with the Mellor‐Yamada‐Janjic PBL scheme.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

On the implementation of dislocation reactions in continuum dislocation dynamics modeling of mesoscale plasticity

The continuum dislocation dynamics framework for mesoscale plasticity is intended to capture the dislocation density evolution and the deformation of crystals when subjected to mechanical loading. It does so by solving a set of transport equations for dislocations concurrently with crystal mechanics equations, with the latter being cast in the form of an eigenstrain problem. Incorporating dislocation reactions in the dislocation transport equations is essential for making such continuum dislocation dynamics predictive. A formulation is proposed to incorporate dislocation reactions in the transport equations of the vector density-based continuum dislocation dynamics. This formulation aims to rigorously enforce dislocation line continuity using the concept of virtual dislocations that close all dislocation loops involved in cross slip, annihilation, and glissile and sessile junction reactions. The addition of virtual dislocations enables us to accurately enforce the divergence free condition upon the numerical solution of the dislocation transport equations for all slip systems individually. A set of tests were performed to illustrate the accuracy of the formulation and the solution of the transport equations within the vector density-based continuum dislocation dynamics. Comparing the results from these tests with an earlier approach in which the divergence free constraint was enforced on the total dislocation density tensor or the sum of two densities when only cross slip is considered shows that the new approach yields highly accurate results. Bulk simulations were performed for a face centered cubic crystal based on the new formulation and the results were compared with discrete dislocation dynamics predictions of the same. The microstructural features obtained from continuum dislocation dynamics were also analyzed with reference to relevant experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Trends in Warm Season Mesoscale Convective Systems Over Asia in 2001–2020

Mesoscale convective systems (MCSs) frequently occur over Asia during the warm season, often producing intense precipitation with associated socioeconomic impacts. Here we reveal significant trends in MCS occurrence frequency and related precipitation in Asia during the warm season (March–September) in 2001–2020, using a tracking method that combines cloud and precipitation criteria with high-resolution satellite data from the Global Precipitation Measurement mission. To examine whether there are differences between MCSs of different scales, both meso-α scales (MαCSs) and meso-β scales (MβCSs), with horizontal scales of 200–2,000 km and 20–200 km, are tracked. The distribution pattern of frequency and related precipitation of both MαCSs and MβCSs are quite similar and manifest positive trends over East Asia (EA) and Northeast Asia, and negative trend over Southeast Asia (SEA). The MCS precipitation trend contributes significantly to total precipitation trend, with MαCSs contributing the most. Our analysis indicates the trend in lower-tropospheric water vapor flux convergence has a similar spatial pattern to the MCS frequency and related precipitation trend. Based on an atmospheric moisture flux decomposition analysis, the water vapor flux convergence trend can largely be explained by the change in horizontal wind convergence, while the specific humidity trend driven largely by temperature change plays a minor role. In conclusion, the trend in wind convergence in EA and SEA is possibly related to the evident trend in the lower-tropospheric anticyclone over the western North Pacific and SEA, which might be due to the relatively stronger warming in the Indian Ocean during the past two decades.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Multimode Metastructures: Novel Hybrid 3D Lattice Topologies

With the rapid proliferation of additive manufacturing and 3D printing technologies, architected cellular solids including truss-like 3D lattice topologies offer the opportunity to program the effective material response through topological design at the mesoscale. The present report summarizes several of the key findings from a 3-year Laboratory Directed Research and Development Program. The program set out to explore novel lattice topologies that can be designed to control, redirect, or dissipate energy from one or multiple insult environments relevant to Sandia missions, including crush, shock/impact, vibration, thermal, etc. In the first 4 sections, we document four novel lattice topologies stemming from this study: coulombic lattices, multi-morphology lattices, interpenetrating lattices, and pore-modified gyroid cellular solids, each with unique properties that had not been achieved by existing cellular/lattice metamaterials. The fifth section explores how unintentional lattice imperfections stemming from the manufacturing process, primarily sur face roughness in the case of laser powder bed fusion, serve to cause stochastic response but that in some cases such as elastic response the stochastic behavior is homogenized through the adoption of lattices. In the sixth section we explore a novel neural network screening process that allows such stocastic variability to be predicted. In the last three sections, we explore considerations of computational design of lattices. Specifically, in section 7 using a novel generative optimization scheme to design novel pareto-optimal lattices for multi-objective environments. In section 8, we use computational design to optimize a metallic lattice structure to absorb impact energy for a 1000 ft/s impact. And in section 9, we develop a modified micromorphic continuum model to solve wave propagation problems in lattices efficiently.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Km‐Scale Simulations of Mesoscale Convective Systems Over South America—A Feature Tracker Intercomparison

Mesoscale convective systems (MCSs) are clusters of thunderstorms that are important in Earth's water and energy cycle. Additionally, they are responsible for extreme events such as large hail, strong winds, and extreme precipitation. Automated object-based analyses that track MCSs have become popular since they allow us to identify and follow MCSs over their entire life cycle in a Lagrangian framework. This rise in popularity was accompanied by an increasing number of MCS tracking algorithms, however, little is known about how sensitive analyses are concerning the MCS tracker formulation. Here, we assess differences between six MCS tracking algorithms on South American MCS characteristics and evaluate MCSs in kilometer-scale simulations with observational-based MCSs over 3 years. All trackers are run with a common set of MCS classification criteria to isolate tracker formulation differences. The tracker formulation substantially impacts MCS characteristics such as frequency, size, duration, and contribution to total precipitation. The evaluation of simulated MCS characteristics is less sensitive to the tracker formulation and all trackers agree that the model can capture MCS characteristics well across different South American climate zones. Dominant sources of uncertainty are the segmentation of cloud systems in space and time and the treatment of how MCSs are linked in time. Our results highlight that comparing MCS analyses that use different tracking algorithms is challenging. We provide general guidelines on how MCS characteristics compare between trackers to facilitate a more robust assessment of MCS statistics in future studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Impact of Momentum Perturbation on Convective Boundary Layer Turbulence

Mesoscale-to-microscale coupling is an important tool for conducting turbulence-resolving multiscale simulations of realistic atmospheric flows, which are crucial for applications ranging from wind energy to wildfire spread studies. Different techniques are used to facilitate the development of realistic turbulence in the large-eddy simulation (LES) domain while minimizing computational cost. Here, we explore the impact of a simple and computationally efficient Stochastic Cell Perturbation method using momentum perturbation (SCPM-M) to accelerate turbulence generation in boundary-coupled LES simulations using the Weather Research and Forecasting model. We simulate a convective boundary layer (CBL) to characterize the production and dissipation of turbulent kinetic energy (TKE) and the variation of TKE budget terms. Furthermore, we evaluate the impact of applying momentum perturbations of three magnitudes below, up to, and above the CBL on the TKE budget terms. Momentum perturbations greatly reduce the fetch associated with turbulence generation. When applied to half the vertical extent of the boundary layer, momentum perturbations produce an adequate amount of turbulence. However, when applied above the CBL, additional structures are generated at the top of the CBL, near the inversion layer. The magnitudes of the TKE budgets produced by SCPM-M when applied at varying heights and with different perturbation amplitudes are always higher near the surface and inversion layer than those produced by No-SCPM, as are their contributions to the TKE. This study provides a better understanding of how SCPM-M reduces computational costs and how different budget terms contribute to TKE in a boundary-coupled LES simulation.

17 WIND ENERGY↗
DOE OSTI2025

The GFDL‐CM4X Climate Model Hierarchy, Part I: Model Description and Thermal Properties

We present the GFDL‐CM4X (Geophysical Fluid Dynamics Laboratory Climate Model version 4X) coupled climate model hierarchy. The primary application for CM4X is to investigate ocean and sea ice physics as part of a realistic coupled Earth climate model. CM4X utilizes an updated MOM6 (Modular Ocean Model version 6) ocean physics package relative to CM4.0, and there are two members of the hierarchy: one that uses a horizontal grid spacing of 0.25° (referred to as CM4X‐p25) and the other that uses a grid 0.125° (CM4X‐p125). CM4X also refines its atmospheric grid from the nominally 100 km (cubed sphere C96) of CM4.0–50 km (C192). Finally, CM4X simplifies the land model to allow for a more focused study of the role of ocean changes to global mean climate. CM4X‐p125 reaches a global ocean area mean heat flux imbalance of -0.02 W m -2 within $\mathcal{O}$ (150) years in a pre‐industrial simulation, and retains that thermally equilibrated state over the subsequent centuries. This 1850 thermal equilibrium is characterized by roughly 400 ZJ less ocean heat than present‐day, which corresponds to estimates for anthropogenic ocean heat uptake between 1870 and present‐day. CM4X‐p25 approaches its thermal equilibrium only after more than 1000 years, at which time its ocean has roughly 1100 ZJ more heat than its early 21st century ocean initial state. Furthermore, the root‐mean‐square sea surface temperature bias for historical simulations is roughly 20% smaller in CM4X‐p125 relative to CM4X‐p25 (and CM4.0). We offer the mesoscale dominance hypothesis for why CM4X‐p125 shows such favorable thermal equilibration properties.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2005

Lunar and Planetary Science XXXVI, Part 19

The topics include: 1) The abundances of Iron-60 in Pyroxene Chondrules from Unequilibrated Ordinary Chondrites; 2) LL-Ordinary Chondrite Impact on the Moon: Results from the 3.9 Ga Impact Melt at the Landing Site of Appolo 17; 3) Evaluation of Chemical Methods for Projectile Identification in Terrestrial and Lunar Impactites; 4) Impact Cratering Experiments in Microgravity Environment; 5) New Achondrites with High-Calcium Pyroxene and Its implication for Igneous Differentiation of Asteroids; 6) Climate History of the Polar Regions of Mars Deduced form Geologic Mapping Results; 7) The crater Production Function for Mars: A-2 Cumulative Power-Law Slope for Pristine Craters Greater than 5 km in Diameter Based on Crater Distribution for Northern Plains Materials; 8) High Resolution Al-26 Chronology: Resolved Time Interval Between Rim and Interior of a Highly Fractionated Compact Type a CAI from Efremovka; 9) Assessing Aqueous Alteration on Mars Using Global Distributions of K and Th; 10) FeNi Metal Grains in LaPaz Mare Basalt Meteorites and Appolo 12 Basalts; 11) Unique Properties of Lunar Soil for In Situ Resource Utilization on the Moon; 12) U-Pb Systematics of Phosphates in Nakhlites; 13) Measurements of Sound Speed in Granular Materials Simulated Regolith; 14) The Effects of Oxygen, Sulphur and Silicon on the Dihedral Angles Between Fe-rich Liquid Metal and Olivine, Ringwoodite and Silicate Perovskite: Implications for Planetary Core Formation; 15) Seismic Shaking Removal of Craters 0.2-0.5 km in Diameter on Asteroid 433 Eros; 16) Focused Ion Beam Microscoopy of ALH84001 Carbonate Disks; 17) Simulating Micro-Gravity in the Laboratory; 18) Mars Atmospheric Sample Return Instrument Development; 19) Combined Remote LIBS and Raman Spectroscopy Measurements; 20) Unusual Radar Backscatter Properties Along the Northern Rim of Imbrium Basin; 21) The Mars Express/NASAS Project at JPL; 22) The Geology of the Viking 2 Lander Site Revisited; 23) An Impact Genesis for Loki Patera? 24) Mars Polar Cap Edges Tracked over 3 Full Mars Years; 25) Elemental Abundance in Presolar SiC: Comparing Grains Separated by Acid Residue and Gently Separation Procedures; 26) First Results from the Descent Imager/Spectral Radiometer (DISR) Experiment on the Huygens Entry Probe of Titan; 27) Minor Element Behavior of Pallasite Olivine: Understanding Pallasite Thermal History and Chronology; 28) Canonical Anorthite in a Grosnaja Forsterite-bearing CAI; 29) Experimental Evidence for Condensation of 'Astrophysical' Carbonate; 30) Distribution and Classification of Multiple Coronae on Venus; 31) Recognition of Rayed Craters on Mars in THEMIS Thermal Infrared Imagery: Implications for Martian Meteorite Source Regions; 32) Geochemical Modeling of Evaporites on Mars: Insight from Meridiani Planum; 33) Hadean Crustal Processes Revealed from Oxygen Isotopes and U-Th-Pb Depth Profiling of Pre-4.0 Ga Detrital Zircons from Western Australia; 34) On Modeling the Seepage of Water into the Martian Subsurface; 35) Martial Gullies and Groundwater: A Series of Unfortunate Exceptions; 36) Olivine and Carbonate Globules in ALH84001: A Terrestrial Analog, and Implications for Water on Mars; 37) A Reevaluation of Mass Movements Within the Valles Marineris Region of Mars Using MOLA and MOC Data; 38) Evidence of Hydrated 109P/Swift-Tuttle Meteoroids from Meteor Spectroscopy; 39) Cr-54 Anomalies in the Solar System: Their Extent and Origin; 40) Reevaluation of the Mn-53-Cr-53 Systematic in the Basaltic Achondrites; 41) Effective Liquid Metal-Silicate Mixing Upon Shock by Power-Law Droplet Size Scaling in Richtmyer-Meshkov Like Perturbations; 42) Post-Impact Deformation of Impact Craters: Towards a Better Understanding Through the Study of Mjolnir Crater; 43) Cutting Silica Aerogel for Particle Extraction; 44) Liquid Hydrocarbons on Titan's Surface? How Cassini ISS Observations Fit into the Story (So Far); and 45) Mesoscale Simulations of Polar Circulations: Late Spring to Late Summe

Source record↗
DOE OSTI2024

Self-Supervised Cloud Classification

Abstract Low-level marine clouds play a pivotal role in Earth’s weather and climate through their interactions with radiation, heat and moisture transport, and the hydrological cycle. These interactions depend on a range of dynamical and microphysical processes that result in a broad diversity of cloud types and spatial structures, and a comprehensive understanding of cloud morphology is critical for continued improvement of our atmospheric modeling and prediction capabilities moving forward. Deep learning has recently accelerated our ability to study clouds using satellite remote sensing, and machine learning classifiers have enabled detailed studies of cloud morphology. A major limitation of deep learning approaches to this problem, however, is the large number of hand-labeled samples that are required for training. This work applies a recently developed self-supervised learning scheme to train a deep convolutional neural network (CNN) to map marine cloud imagery to vector embeddings that capture information about mesoscale cloud morphology and can be used for satellite image classification. The model is evaluated against existing cloud classification datasets and several use cases are demonstrated, including training cloud classifiers with very few labeled samples, interrogation of the CNN’s learned internal feature representations, cross-instrument application, and resilience against sensor calibration drift and changing scene brightness. The self-supervised approach learns meaningful internal representations of cloud structures and achieves comparable classification accuracy to supervised deep learning methods without the expense of creating large hand-annotated training datasets. Significance Statement Marine clouds heavily influence Earth’s weather and climate, and improved understanding of marine clouds is required to improve our atmospheric modeling capabilities and physical understanding of the atmosphere. Recently, deep learning has emerged as a powerful research tool that can be used to identify and study specific marine cloud types in the vast number of images collected by Earth-observing satellites. While powerful, these approaches require hand-labeling of training data, which is prohibitively time intensive. This study evaluates a recently developed self-supervised deep learning method that does not require human-labeled training data for processing images of clouds. We show that the trained algorithm performs competitively with algorithms trained on hand-labeled data for image classification tasks. We also discuss potential downstream uses and demonstrate some exciting features of the approach including application to multiple satellite instruments, resilience against changing image brightness, and its learned internal representations of cloud types. The self-supervised technique removes one of the major hurdles for applying deep learning to very large atmospheric datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Twin nucleation and growth in hexagonal close-packed metals: The role of slip-mediated plasticity on twin embryo formation and evolution

Twinning is a key deformation mechanism in hexagonal close-packed (hcp) metals, which are typified by a lack of easily activated slip systems that can accommodate a general state of loading. Pragmatically, the nucleation and evolution of twin domains occurs concomitantly with slip, such that the eventual twin network is conditioned by both external and internal stresses resulting, among others, from the evolution of dislocations. However our understanding of the interplay between dislocations, and twin nucleation and stability remains limited. This work focuses on elucidating the influence of dislocation-mediated plasticity on the formation, growth, and stability of twin embryos. First, a new mesoscale spectral crystal plasticity-twinning framework is developed and used to quantify the change in the free energy landscape following the nucleation {$10\bar{1}2$} of twins in Mg. Representative twin morphologies are modeled under conditions of limited and profuse slip activity, which are emulative of small- and bulk-scale samples, respectively. Then, the driving traction profiles around twin embryos are investigated via a sharp interface approach to obtain insights into how concurrent slip-mediated plasticity can influence the growth/stabilization of nanometric twins. The driving traction profiles are further utilized to determine the stability of twin embryos post loading. The initial dislocation density in the samples, and within the different domains (i.e., twin vs. parent), is seen to have a significant effect on the twin nucleation stress. Namely, the activation of high levels of concomitant plasticity in the parent grain is seen to significantly drive the formation of nanometric twin nuclei at stresses as low as 250 MPa. Further, the sharp interface analysis reveals that profuse plasticity in the parent grain simultaneously alters the forward and back stresses, such that the magnitude/polarity of the driving tractions become increasingly favorable for nanometric twin growth when slip is active. Finally, prior plasticity in the parent grain is seen to result in favorable driving tractions for nanometric twin growth, even at applied stresses as low as ~ 100 MPa. In conclusion, these results are in stark contrast to a case without any dislocations, wherein applied stresses as high as ~ 500 MPa are necessary to grow the twin domains.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mesoscale Convective Systems Represented in High Resolution E3SMv2 and Impact of New Cloud and Convection Parameterizations

Mesoscale convective systems (MCSs) play an important role in modulating the global hydrological cycle, general circulation, and radiative energy budget. In this study, we evaluate MCS simulations in the second version of U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SMv2). E3SMv2 atmosphere model (EAMv2) is run at the uniform 0.25? horizontal resolution. We track MCSs consistently in the model and observations using the PyFLEXTRKR algorithm, which defines MCS based on both cloud-top brightness temperature (Tb) and surface precipitation. Results from using Tb only to define MCS, commonly used in previous studies, are also discussed. Furthermore, sensitivity experiments are performed to examine the impact of new cloud and convection parameterizations developed for EAMv3 on simulated MCSs. Our results show that EAMv2 simulated MCS precipitation is largely underestimated in the tropics and contiguous United States. This is mainly attributed to the underestimated precipitation intensity in EAMv2. In contrast, the simulated MCS frequency becomes more comparable to observations if MCSs are defined only based on cloud-top Tb. The Tb-based MCS tracking method, however, includes many cloud systems with very weak precipitation which conflicts with the MCS definition. This result illustrates the importance of accounting for precipitation in evaluating simulated MCSs. We also find that the new physics parameterizations help increase the relative contribution of convective precipitation to total precipitation in the tropics, but the simulated MCS properties are generally not improved. This suggests that simulating MCSs will remain a challenge for the next version of E3SM.

Zhang, Meng↗
DOE OSTI2022

Gap Flows Along the Columbia River Observed During the WFIP2 Field Campaign

Data collected in the Columbia River Gorge and Basin Area during the second Wind Forecasting Improvement Project (WFIP2) was used to study gap flow events in the region. In particular, the relationship between 200 m zonal wind speed and the zonal mesoscale pressure gradient at Troutdale (West of the gorge), Wasco (east of the gorge), and Boardman (farther east of the gorge) sites was examined. An objective criterion was used to identify 169 (Troutdale), 161 (Wasco), and 113 (Boardman) gap flow days that had weak synoptic forcing. The 200 m zonal wind speed exhibited a linear relationship with the zonal pressure gradient at Wasco during both easterly and westerly gap flow days at daily timescales. Such a relationship was only observed during easterly gap flow days at Troutdale and was not observed at Boardman. The zonal wind speed exhibited a diurnal cycle at Wasco and Boardman, with greater changes in the nighttime winds than the daytime winds due to changes in pressure gradient. The root mean square differences (RMSDs) between the observed pressure gradient at Wasco and that simulated by the High Resolution Rapid Refresh (HRRR) model in 3 km (HRRR-3km) and 750 m (HRRR-750m) horizontal resolution were 3.14 hPa and 0.57 hPa, respectively. However, the RMSDs for the 200 m zonal wind speed were 1.40 m s -1 for HRRR-3km and 1.20 m s -1 for HRRR- 750m. Improvement in simulating the zonal pressure gradient accompanied by negligible changes in simulating zonal wind speed due to higher horizontal resolution points to a complex interplay between different components of the model physics when simulating the meteorology. The HRRR-750m resolution model accurately mimicked the observed covariability (phase diagram) between the zonal pressure gradient and the 200 m zonal wind speed during the nighttime, but not during the daytime, suggesting model errors might be related to the representation of boundary layer processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Phase-field simulation of domain size effect on dielectric and piezoelectric responses in K 0.5 Na 0.5 NbO 3 epitaxial thin films with superdomain structures

Size effects of mesoscale ferroelectric domains on the macroscopic dielectric and piezoelectric responses in domain-engineered bulk piezocrystals have been extensively studied for more than two decades. However, less is known about the domain size effects in ferroelectric epitaxial thin films, especially for films consisting of low-symmetry ferroelectric phases and exhibits hierarchical superdomain structures. Herein, using phase-field simulations, we systemically evaluate the effective out-of-plane dielectric and piezoelectric coefficients, κ$^{*}_{33}$ and d$^{*}_{33}$, as a function of the domain periods for two types of superdomain structures in ferroelectric K 0.5 Na 0.5 NbO 3 epitaxial thin films. In one type of the superdomain structures, we find that more than 70% increase of κ$^{*}_{33}$ and nearly 20% increase of d$^{*}_{33}$ can be achieved by tuning the domain period by a few tens of nanometers. Dissimilar behaviors are found when the domain period varies along different lateral directions, suggesting anisotropic domain size effects in thin films. By analyzing the local dielectric and piezoelectric responses from each domain variants and domain walls, we reveal that the domain size effect is governed by the variation of out-of-plane polarization inside the domains. Moreover, we also demonstrate enhanced domain size effects by modulating the misfit strains and temperature to approach the polymorphic phase boundaries, suggesting tunability of the size effect by doping and strain engineering. Furthermore, our results reveal multimodal domain size dependence of dielectric and piezoelectric responses in low-symmetry ferroelectric epitaxial thin films, implying that domain size engineering can be used to tune macroscopic properties of thin-film ferroelectrics, similar to their bulk counterparts.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluating Mesoscale Convective Systems Over the US in Conventional and Multiscale Modeling Framework Configurations of E3SMv1

Organized mesoscale convective systems (MCSs) contribute a significant amount of precipitation in the Central and Eastern US during spring and summer, which impacts the availability of freshwater and flooding events. However, current global Earth system models cannot capture MCSs well and misrepresent the statistics of precipitation in the region. In this study, we investigate the representation of MCSs in three configurations of the Energy Exascale Earth System Model (E3SMv1) by tracking individual storms based on outgoing longwave radiation using a new application of TempestExtremes. Our results indicate that conventional parameterizations of convection, implemented in both low (LR; ~150 km) and high (HR; ~25 km) resolution configurations, fail to capture almost all MCS-like events, in-part because they underestimate high-level cloud ice associated with deep convection. On the other hand, the multiscale modeling framework (MMF; cloud-resolving models embedded in each grid-column of ~150 km resolution E3SMv1) configuration represents MCSs and their annual cycle better. Nevertheless, relative to observations, the E3SMv1-MMF spatial distribution of MCSs and associated precipitation is shifted eastward, and the diurnal timing is lagged. A comparison between the large-scale environment in E3SMv1-MMF and ERA5 reanalysis suggests that the biases during the summer in E3SMv1-MMF are associated with biases in low-level humidity and meridional moisture transport within the low-level jet. The fact that conventional parameterizations of convection, even with high-resolution, cannot capture MCSs over the US suggests that methods with explicit representation of kilometer-scale convective organization, such as the MMF, may be necessary for improving the simulation of these convective systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Use of Remote Sensing and In-Situ Observations to Develop and Evaluate Improved Representations of Convection and Clouds for the ACME Model

The overachieving goal of the whole CMDV-MCS project is to improve understanding of warm season continental convection and to develop treatments of convection and microphysics capable of representing mesoscale convective systems (MCSs) features in large-scale models. Our tasks for this project contributing to the overachieving goal include: (1) Improve the ice nucleation formulation for MG2 and P3 cloud microphysics schemes; (2) Improve the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM; and (3) Test the performance of improved ice microphysics in E3SM with observation data. In this project, we have (1) Improved the ice nucleation parameterization for MG2 and P3 in E3SM by implementing two advanced empirical parameterizations with connection to aerosols. The two deterministic heterogeneous ice nucleation parameterizations (i.e., DeMott et al., 2015; Niemand et al., 2012) were merged with the MG2 and P3 cloud microphysics schemes in E3SM. Long-term simulations were conducted to examine the impacts of these new parameterizations on simulated cloud properties; (2) Improved the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM. We evaluated the double Gaussian PDF of vertical velocity simulated by the Cloud Layers Unified By Binormals (CLUBB) and the sub-column vertical velocity sampled from the Subgrid Importance Latin Hypercube Sampler (SILHS) in E3SM. We introduced the vertical velocity variance induced by topographic gravity waves for ice nucleation and droplet activation; and (3) Tested the performance of improved ice microphysics in E3SM with observation data. We tested the new treatments of ice nucleation in the single column model (SCM) mode for the stratiform mixed-phase clouds observed during 9-10 October 2004 in the DOE ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) and for the convective clouds observed on 20 May 2011 in the Midlatitude Continental Convective Clouds Experiment (MC3E). Modeled ice nucleating particles (INPs) concentrations were compared against observations collected around the globe.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A protocol and analysis of year-long simulations of global storm-resolving models and beyond

We propose a protocol to evaluate and analyze year-long simulations of global storm-resolving models (GSRMs). The proposed protocol complements an earlier 40-day simulation protocol under the DYAMOND (DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains) project to allow the analysis of the seasonal cycle and associated climatic relevant phenomena. This intercomparison aims to reveal how GSRMs, which can simulate mesoscale convective systems (MCSs) in the global domain, reproduce atmospheric large-scale structures related to convection beyond month-long simulations. The intercomparison for one-year simulations is conducted by either atmosphere-only models or atmosphere–ocean coupled models with atmospheric horizontal mesh sizes less than 5 km. We recommend the continuous four seasons from March 2020 to February 2021 as a target period for the intercomparison but with options for many groups to join more flexibly. The output variables are collected at 0.25° resolution, and archives of a small set of native grid variables are encouraged to analyze tropical cyclones and MCSs. Through the proposed global storm-resolving simulation, we will evaluate the climatological distributions of the atmospheric large-scale circulations, such as the Intertropical Convergence Zone (ITCZ), monsoon, midlatitude jets, their time evolution, and the upscale impacts on them. We present sample analyses from a one-year simulation using the 3.5 km mesh Nonhydrostatic Icosahedral Atmospheric Model (NICAM), revealing the realistic zonal contrast of tropical precipitation, no double ITCZ structure, the reasonable midlatitude jet position and intensity but a weak bias of storm track activities, and a warm bias over the Eurasia during boreal winter. We also clarify the cross-scale interaction, such as the effects of cold pools on mean precipitation over the Maritime Continent through the precipitation diurnal cycle and the effects of resolved gravity waves on midlatitude mean flows. The proposed one-year simulation protocol is referred to as the “Sendai Protocol.” This protocol is not unique or definite for evaluating GSRMs; we prospect a hierarchical set of experiments from short-term to multi-year simulations as GSRM intercomparisons.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Kinetics and Thermodynamics of Sr Permeation in CeO 2 -Based Barrier Layers for Solid-Oxide Electrolyzer Cells

Solid-oxide electrolyzer cells (SOECs) convert steam to hydrogen efficiently at high temperatures. However, during operation, the diffusion of cations or impurities through the cells due to electrode degradation can cause unwanted secondary phases to form, which may degrade device performance. Here, in this study, we use atomistic and mesoscale simulations coupled with experimental analysis to study the diffusion of Sr through the Gd-doped CeO 2 (GDC) barrier layer used to protect the yttria-stabilized zirconia (YSZ) electrolyte in SOECs. From our atomistic calculations, we find Sr diffusion to be negligibly slow in bulk GDC; however, surface diffusion is much more favorable. Subsequent mesoscale simulations show that Sr diffusion is activated when the porosity of GDC exceeds ∼10% and significantly exceeds diffusion in bulk and grain boundary regions. We also find that SrO-based species can accumulate at GDC surfaces; however, SrO aggregation and coarsening will be limited by the large lattice mismatch between GDC and SrO. Energy-dispersive X-ray spectroscopy (EDS) and electron diffraction confirm that Sr can accumulate within GDC pores and form disperse Sr-containing secondary phases. Altogether, Sr diffusion in dense GDC is unlikely to give rise to thick SrO layers, which would severely limit device performance. The formation of Sr-containing secondary phases can largely be avoided by restricting the porosity of the GDC layer as much as possible.

36 MATERIALS SCIENCE↗