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

The nano- and meso-scale structure of amorphous calcium carbonate

Abstract Understanding the underlying processes of biomineralization is crucial to a range of disciplines allowing us to quantify the effects of climate change on marine organisms, decipher the details of paleoclimate records and advance the development of biomimetic materials. Many biological minerals form via intermediate amorphous phases, which are hard to characterize due to their transient nature and a lack of long-range order. Here, using Monte Carlo simulations constrained by X-ray and neutron scattering data together with model building, we demonstrate a method for determining the structure of these intermediates with a study of amorphous calcium carbonate (ACC) which is a precursor in the bio-formation of crystalline calcium carbonates. We find that ACC consists of highly ordered anhydrous nano-domains of approx. 2 nm that can be described as nanocrystalline. These nano-domains are held together by an interstitial net-like matrix of water molecules which generate, on the mesoscale, a heterogeneous and gel-like structure of ACC. We probed the structural stability and dynamics of our model on the nanosecond timescale by molecular dynamics simulations. These simulations revealed a gel-like and glassy nature of ACC due to the water molecules and carbonate ions in the interstitial matrix featuring pronounced orientational and translational flexibility. This allows for viscous mobility with diffusion constants four to five orders of magnitude lower than those observed in solutions. Small and ultra-small angle neutron scattering indicates a hierarchically-ordered organization of ACC across length scales that allow us, based on our nano-domain model, to build a comprehensive picture of ACC formation by cluster assembly from solution. This contribution provides a new atomic-scale understanding of ACC and provides a framework for the general exploration of biomineralization and biomimetic processes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

On the Prospect of Chemically Transferable Coarse-Grained Electronic Models for Soft Materials

Electronic coarse-graining (ECG) methods predict quantum-mechanical electronic properties directly from coarse-grained (CG) molecular configurations, enabling electronic predictions at mesoscale length scales. Here, we present a diagnostic assessment of the feasibility of chemically transferable ECG models across a broad polymer-relevant chemical space using all-atom, united-atom, and Martini-scale representations. While high-resolution ECG models achieve near-quantitative accuracy, we show that chemically transferable ECG at the Martini resolution fails because the CG force field does not sample the same configurational distribution of local molecular structure as that underlying the DFT-parameterized ECG model. We demonstrate that our proposed Element-Count-Label (ECL) representation, which augments Martini beads with explicit stoichiometric data, significantly improves chemical generalization across diverse polymer chemistries. However, we find that even with improved chemical resolution, the model cannot recover electronic property distributions that are absent from the configurational space sampled by the CG force field. These results demonstrate that chemically transferable ECG requires future Martini-like force fields to explicitly preserve quantum chemistry–compatible local molecular structure in addition to thermodynamic and structural fidelity.

Kidder, Katherine M [Department of Chemistry; Univ↗
DOE OSTI2021

Domain Wall Patterning and Giant Response Functions in Ferrimagnetic Spinels

Abstract The manipulation of mesoscale domain wall phenomena has emerged as a powerful strategy for designing ferroelectric responses in functional devices, but its full potential is not yet realized in the field of magnetism. This work shows a direct connection between magnetic response functions in mechanically strained samples of Mn 3 O 4 and MnV 2 O 4 and stripe‐like patternings of the bulk magnetization which appear below known magnetostructural transitions. Building off previous magnetic force microscopy data, a small‐angle neutron scattering is used to show that these patterns represent distinctive magnetic phenomena which extend throughout the bulk of two separate materials, and further are controllable via applied magnetic field and mechanical stress. These results are unambiguously connected to the anomalously large magnetoelastic and magnetodielectric response functions reported for these materials, by performing susceptibility measurements on the same crystals and directly correlating local and macroscopic data.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Long‐Term Large‐Scale Atmospheric Forcing Data From Three‐Dimensional Constrained Variational Analysis for the ARM SGP Site

Here, this study presents a long‐term three‐dimensional large‐scale forcing data set (VARANAL3D) derived from the three‐dimensional constrained variational analysis (3DCVA) method at the Atmospheric Radiation Measurement (ARM) program Southern Great Plains (SGP) site from 2004 to 2018. Building on the same input data sets as the conventional continuous forcing data set (VARANAL), VARANAL3D maintains overall consistency in domain‐averaged fields while introducing spatial variability, offering critical insights into the influence of mesoscale synoptic systems on cloud‐related processes. Evaluations are conducted across four cloud and precipitation regimes: Clear‐sky, Shallow‐clouds, Afternoon‐precipitation, and Nocturnal‐precipitation, presenting high consistency of the domain‐mean forcing data sets while emphasizing the role of subdomain forcing variability particularly in precipitating regimes. Single column model (SCM) simulations demonstrate that subdomain VARANAL3D forcing improves cloud and precipitation representation, with the ensemble outperforming domain‐mean forcing in three cloudy and precipitating regimes. Overall, these results highlight VARANAL3D's value for investigating the impacts of spatial variability of large‐scale forcing on atmospheric processes. The VARANAL3D data set provides new opportunities for evaluating model physics, advancing the development of scale‐aware parameterizations and deepening our understanding of cloud and precipitation dynamics.

Environmental sciences↗
DOE OSTI2025

Cross slip of extended dislocations in face-centered cubic metals through phase-field modeling

Cross slip is a dislocation mechanism that significantly impacts the mechanical behavior of engineering alloys. Here, in this work, we advance a 3D phase-field dislocation dynamics (PFDD) mesoscale technique to simulate cross slip across a broad range of face-centered cubic (FCC) metals. The formulation incorporates elastic anisotropy, an FCC numerical grid, and a high-fidelity representation of the entire γ -surface from density functional theory for eight FCC metals and no adjustable parameters or rules. The relaxed core structures under zero stress for all metals are predicted to extend in plane. The analytical model for stacking fault width agrees well with the PFDD result under the assumption of elastic isotropy but overestimates it under elastic anisotropy, when the degree of anisotropy is large. The dynamic simulations are designed to elucidate the material parameters that influence the propensity for cross slip. Whether cross slip occurs under a non-Schmid stress or to bypass a hard obstacle, the critical stress to cross slip scales strongly with the anisotropic energy coefficient for a screw dislocation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Parameterization of the Elevated Convection With a Unified Convection Scheme (UNICON) and Its Impacts on the Diurnal Cycle of Precipitation

To improve the simulation of nocturnal precipitation, we develop a parameterization for the elevated convection that is launched from the level of maximum grid-mean moist static energy, when grid-mean vertical flow is upward at the launching interface. The parameterized elevated convection is forced by both cold pool-driven and partially resolved external mesoscale organized flows. Properties of the external mesoscale flow are estimated from grid-mean vertical velocity and three dimensional advection tendencies of temperature and moisture. The new parameterization is implemented into a unified convection scheme (UNICON) and tested for both the single-column case at the Southern Great Plain (SGP) site in US and in global simulations. At the SGP site, the parameterized elevated convection strengthens nocturnal convection, increases nocturnal precipitation, and better simulates the observed diurnal cycle of precipitation. It appears that the elevated and surface-based convections interact with each other by stabilizing the atmospheric column, which affects subsequent convection. Global simulation shows that the elevated convection mostly occurs over the continents during the night, and also over the oceanic mid-latitude warm air advection and storm track regions during summer. Without degrading global mean climate, the elevated convection improves the simulation of nocturnal precipitation in the mid-US but simulates somewhat strong nocturnal precipitation in northeastern Asia. It seems that a key to simulate observed nocturnal precipitation is to appropriately parameterize the impacts of external organized flow on the elevated convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Compound Mesoscale Convective Systems and Low‐Pressure Systems in Tropical Monsoon Regions: Assessing Their Meteorology and Precipitation

Mesoscale convective systems (MCS) and low-pressure systems (LPS) are both strongly associated with precipitation across the regions where they occur, particularly within global monsoon systems; however, their co-occurrence and its relationship to precipitation have not been systematically examined. Here, we use LPS and MCS trackers to detect compound MCS and LPS events in five monsoon regions and assess the association of this co-occurrence with anomalies of winds, precipitation, and other atmospheric variables. Additionally, we investigate the spatial distribution of precipitating MCS and LPS events. Our results show that most (∼60%) MCS and LPS co-occurrences are located in the lower latitudes, where they contribute up to 40% of annual precipitation. We find that compound events generally produce more extreme precipitation than MCS-only or LPS-only events. Furthermore, our assessment of the synoptic and mesoscale composites reveals that the underlying dynamics of compound events exhibit anomalously positive convective available potential energy and an anomalously low-pressure within the location of the event. In terms of the synoptic environment of the features, the compound MCS and LPS events are associated with inverted troughs in three out of five monsoon locations assessed.

Environmental sciences↗
DOE OSTI2025

Emergent topological polarization textures in relaxor ferroelectrics

Relaxor ferroelectrics underpin high-performance actuators and sensors, yet the nature of polar heterogeneities driving their broadband dielectric response remains debated. Using a unified, multimodal structural refinement framework— simultaneously fitting complementary X-ray and neutron total scattering, X-ray absorption spectra, and diffuse scattering—we reconstruct 3D mesoscale polarization maps in the classic relaxor system PbMg 1/3 Nb 2/3 O 3 –PbTiO 3 . We uncover self-organized swirling polarization textures with half-skyrmion (meron) vortices, challenging models of independent polar nanoregions. These textures, characterized by smooth changes in the polarization direction, originate from overlapping volumes in which the projections of locally correlated polarization vectors onto each volume’s long axis share the same sign. Vortex cores correlate strongly with local charge and strain gradients imposed by compositional heterogeneities. In this work, our results suggest that chemical disorder, acting via depolarizing and strain fields, stabilizes topological vortex textures of the polarization field, offering a route for engineering new dielectric and ferroelectric functionalities.

Eremenko, Maksim [National Institute of Standards ↗
DOE OSTI2020

Large anomalous Nernst and inverse spin-Hall effects in epitaxial thin films of kagome semimetal Mn 3 Ge

Synthesis of crystallographically well-defined thin films of topological materials is important for unraveling their mesoscale quantum properties and for device applications. Mn 3 Ge , an antiferromagnetic Weyl semimetal with a chiral magnetic structure on a kagome lattice, is expected to have enhanced Berry curvature around Weyl nodes near the Fermi energy, leading to large anomalous Hall/Nernst effects and a large spin-Hall effect. Using magnetron sputtering, we have grown epitaxial thin films of hexagonal D 0 19 Mn 3 Ge that are flat and continuous. Large anomalous Nernst and inverse spin-Hall effects are observed in thermoelectric and spin-pumping devices. The anomalous Nernst signal in our Mn 3 Ge films is estimated to be 0.1 μV/K and is comparable to that in ferromagnetic Fe, despite Mn 3 Ge having a weak magnetization of ~ 3.5 m μ B / Mn at room temperature. In this work, the spin-mixing conductance is 90.5 nm – 2 at the Py / Mn 3 Ge interface, and the spin-Hall angle in Mn 3 Ge is estimated to be about eight times of that in Pt.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Changing Characteristics of Tropical Extreme Precipitation–Cloud Regimes in Warmer Climates

In this study, we investigated the changing characteristics of climatic scale (monthly) tropical extreme precipitation in warming climates using the Energy Exascale Earth System Model (E3SM). The results are from Atmospheric Model Intercomparison Project (AMIP)-type simulations driven by (a) a control experiment with the present-day sea surface temperature (SST) and CO 2 concentration, (b) P4K, the same as in (a) but with a uniform increase of 4K in the SST globally, and (c) the same as in (a), but with an imposed SST and CO 2 concentration from the outputs of the coupled E3SM forced by a 4xCO 2 concentration. We found that as the surface warmed under P4K and 4xCO 2 , both convective and stratiform rain increased. Importantly, there was an increasing fractional contribution of stratiform rain as a function of the precipitation intensity, with the most extreme but rare events occurring preferentially over land more than the ocean, and more so under 4xCO 2 than P4K. Extreme precipitation was facilitated by increased precipitation efficiency, reflecting accelerated rates of recycling of precipitation cloud water (both liquid and ice phases) in regions with colder anvil cloud tops. Changes in the vertical profiles of clouds, condensation heating, and vertical motions indicate increasing precipitation–cloud–circulation organization from the control and P4K to 4xCO 2 . The results suggest that large-scale ocean warming, that is, P4K, was the primary cause contributing to an organization structure resembling the well-known mesoscale convective system (MCS), with increased extreme precipitation on shorter (hourly to daily) time scales. Additional 4xCO 2 atmospheric radiative heating and dynamically consistent anomalous SST further amplified the MCS organization under P4K. Analyses of the surface moist static energy distribution show that increases in the surface moisture (temperature) under P4K and 4xCO 2 was the key driver leading to enhanced convective instability over tropical ocean (land). However, a fast and large increase in the land surface temperature and lack of available local moisture resulted in a strong reduction in the land surface relative humidity, reflecting severe drying and enhanced convective inhibition (CIN). It is argued that very extreme and rare “record-breaking” precipitation events found over land under P4K, and more so under 4xCO 2 , are likely due to the delayed onset of deep convection, that is, the longer the suppression of deep convection by CIN, the more severe the extreme precipitation when it eventually occurs, due to the release of a large amount of stored surplus convective available potential energy in the lower troposphere during prolonged CIN.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

Hurricane Data Analysis Tool

In order to facilitate Earth science data access, the NASA Goddard Earth Sciences Data Information Services Center (GES DISC) has developed a web prototype, the Hurricane Data Analysis Tool (HDAT; URL: http://disc.gsfc.nasa.gov/HDAT), to allow users to conduct online visualization and analysis of several remote sensing and model datasets for educational activities and studies of tropical cyclones and other weather phenomena. With a web browser and few mouse clicks, users can have a full access to terabytes of data and generate 2-D or time-series plots and animation without downloading any software and data. HDAT includes data from the NASA Tropical Rainfall Measuring Mission (TRMM), the NASA Quick Scatterometer(QuikSCAT) and NECP Reanalysis, and the NCEP/CPC half-hourly, 4-km Global (60 N - 60 S) IR Dataset. The GES DISC archives TRMM data. The daily global rainfall product derived from the 3-hourly multi-satellite precipitation product (3B42 V6) is available in HDAT. The TRMM Microwave Imager (TMI) sea surface temperature from the Remote Sensing Systems is in HDAT as well. The NASA QuikSCAT ocean surface wind and the NCEP Reanalysis provide ocean surface and atmospheric conditions, respectively. The global merged IR product, also known as, the NCEP/CPC half-hourly, 4-km Global (60 N -60 S) IR Dataset, is one of TRMM ancillary datasets. They are globally-merged pixel-resolution IR brightness temperature data (equivalent blackbody temperatures), merged from all available geostationary satellites (GOES-8/10, METEOSAT-7/5 & GMS). The GES DISC has collected over 10 years of the data beginning from February of 2000. This high temporal resolution (every 30 minutes) dataset not only provides additional background information to TRMM and other satellite missions, but also allows observing a wide range of meteorological phenomena from space, such as, hurricanes, typhoons, tropical cyclones, mesoscale convection system, etc. Basic functions include selection of area of interest and time, single imagery, overlay of two different products, animation,a time skip capability and different image size outputs. Users can save an animation as a file (animated gif) and import it in other presentation software, such as, Microsoft PowerPoint. Since the tool can directly access the real data, more features and functionality can be added in the future.

Liu, Zhong↗
DOE OSTI2026

A mechanistic, multiscale model for predicting Pd penetration in TRISO fuels using BISON

TRistructural ISOtropic (TRISO) particles use silicon carbide (SiC) as the primary structural member and barrier against metallic fission product (FP) release. palladiums (PDs), produced by fission in the fuel kernel, can diffuse to and chemically interact with the SiC layer, degrading its structural integrity and ability to contain radioactive FPs. Existing temperature-dependent correlations for predicting Pd penetration rely on experimental data with significant scatter due to varying conditions, potentially complicating ongoing fuel qualification and licensing efforts for advanced reactors that would subject TRISO fuels to operating conditions outside of those examined in the experiments. A mechanistic model of Pd production, transport, and reaction is developed in this work to better understand and predict PD attack of SiC in TRISO particles. molecular dynamicss (MDs) simulations are utilized to calculate the Pd diffusivity in SiC grain bulk and grain boundaries. A mesoscale phase-field diffusion model, informed by the MD diffusivities, is used to develop a reduced order model (ROM) for the effect of SiC microstructure on the Pd penetration rate. The engineering scale BISON model calculates Pd production and transport, utilizing the ROM to predict penetration rates consistent with experimental data. This novel mechanistic ROM captures the effect of temperature, microstructure, and irradiation history on Pd penetration. In conclusion, these new capabilities are expected to support ongoing qualification and licensing efforts associated with near-term TRISO-fueled reactor applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Distinguishing Mesoscale Polar Order (Unidirectional vs Bidirectional) of Cellulose Microfibrils in Plant Cell Walls Using Sum Frequency Generation Spectroscopy

Cellulose in plant cell walls are synthesized as crystalline microfibrils with diameters of 3–4 nm and lengths of around 1–10 μm. These microfibrils are known to be the backbone of cell walls, and their multiscale three-dimensional organization plays a critical role in cell wall functions including plant growth and recalcitrance to degradation. The mesoscale organization of microfibrils over a 1–100 nm range in cell walls is challenging to resolve because most characterization techniques investigating this length scale suffer from low spatial resolution, sample preparation artifacts, or inaccessibility of specific cell types. Here, we report a sum frequency generation (SFG) study determining the mesoscale polarity of cellulose microfibrils in intact plant cell walls. SFG is a nonlinear optical spectroscopy technique sensitive to the molecular-to-mesoscale order of noncentrosymmetric domains in amorphous matrices. However, the quantitative theoretical model to unravel the effect of polarity in packing of noncentrosymmetric domains on SFG spectral features has remained unresolved. In this work, we show how the phase synchronization principle of the SFG process is used to predict the relative intensities of vibrational modes with different polar angles from the noncentrosymmetric domain axis. Applying this model calculation for the first time and employing SFG microscopy, we found that cellulose microfibrils in certain xylem cell walls are deposited unidirectionally (or biased in one direction) instead of the bidirectional polarity which was believed to be dominant in plant cell walls from volume-averaged characterizations of macroscopic samples. Here, with this advancement in SFG analysis, one can now determine the relative polarity of noncentrosymmetric domains such as crystalline biopolymers interspersed in amorphous polymer matrices, which will open opportunities to study new questions that have not been conceived in the past.

36 MATERIALS SCIENCE↗
DOE OSTI2025

SULI Intern Final Report: Computationally Investigating Hydrogen Thermo-Diffusion in Yttrium Hydride Using Multiscale Methods

The renaissance of nuclear energy has arrived, heralding an age of abundant inexpensive clean energy, and renewed space exploration. In nuclear-powered spacecraft and microreactors, safety and size are of utmost importance. Yttrium Hydride (YHx) is being researched for its utility as a neutron moderator in nuclear reactors; the hydrogen in YHx slows down neutrons, enabling a continuous nuclear reaction in the reactor. This has the benefit of allowing reactors to be more safe, compact, and efficient. The goal of this effort is to computationally predict the coefficient of temperature-dependent hydrogen diffusion within YHx, the Soret coefficient. This parameter is essential for determining the safe operating modes of YHx moderators. Zirconium Hydride (ZrHx) is used in the Training, Research, Isotopes, General Atomics (TRIGA) reactor, is the reference material for these calculations. In this work, nanoscale atomic modeling in the Vienna Ab initio Simulation Package (VASP) is combined with the mesoscale finite element phase-field module in the Multiphysics Object-Oriented Simulation Environment (MOOSE); this culminates in a new multiscale computational method to simulate Soret diffusion of hydrogen in YHx. This data is useful for predicting experimental outcomes. This workflow involves convergence testing followed by static, Nudged Elastic Band (NEB), Quasi-Harmonic Approximation (QHA), and Molecular Dynamics (MD) calculations - linked with phase field simulation. NEB simulates hydrogen migration, while QHA and MD predict temperature-dependent properties. The static calculations align with literature, and preliminary NEB and QHA calculations yield accurate results. Once the atomic calculations are complete, we will incorporate Electron Backscatter Diffraction (EBSD) images and VASP-generated parameters into the phase field module to simulate intra- and intergranular transport of hydrogen in ZrHx and YHx. Future research will extend our approach to fuel-moderator materials systems such as Uranium-Yttrium Hydride (U-YHx). This work contributes to the development of advanced nuclear energy solutions for space travel.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Mesoscale structural gradients in human tooth enamel

The outstanding mechanical and chemical properties of dental enamel emerge from its complex hierarchical architecture. An accurate, detailed multiscale model of the structure and composition of enamel is important for understanding lesion formation in tooth decay (dental caries), enamel development (amelogenesis) and associated pathologies (e.g., amelogenesis imperfecta or molar hypomineralization), and minimally invasive dentistry. Although features at length scales smaller than 100 nm (individual crystallites) and greater than 50 µm (multiple rods) are well understood, competing field of view and sampling considerations have hindered exploration of mesoscale features, i.e., at the level of single enamel rods and the interrod enamel (1 to 10 µm). Here, we combine synchrotron X-ray diffraction at submicrometer resolution, analysis of crystallite orientation distribution, and unsupervised machine learning to show that crystallographic parameters differ between rod head and rod tail/interrod enamel. This variation strongly suggests that crystallites in different microarchitectural domains also differ in their composition. Thus, we use a dilute linear model to predict the concentrations of minority ions in hydroxylapatite (Mg 2+ and CO 3 2- /Na + ) that plausibly explain the observed lattice parameter variations. While differences within samples are highly significant and of similar magnitude, absolute values and the sign of the effect for some crystallographic parameters show interindividual variation that warrants further investigation. In conclusion, by revealing additional complexity at the rod/interrod level of human enamel and leaving open the possibility of modulation across larger length scales, these results inform future investigations into mechanisms governing amelogenesis and introduce another feature to consider when modeling the mechanical and chemical performance of enamel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Nematic Order, Plasmonic Switching and Self-Patterning of Colloidal Gold Bipyramids

Dispersing inorganic colloidal nanoparticles within nematic liquid crystals provides a versatile platform both for forming new soft matter phases and for predefining physical behavior through mesoscale molecular-colloidal self-organization. However, owing to formation of particle-induced singular defects and complex elasticity-mediated interactions, this approach has been implemented mainly just for colloidal nanorods and nanoplatelets, limiting its potential technological utility. Here, orientationally ordered nematic colloidal dispersions are reported of pentagonal gold bipyramids that exhibit narrow but controlled polarization-dependent surface plasmon resonance spectra and facile electric switching. Bipyramids tend to orient with their C 5 rotation symmetry axes along the nematic director, exhibiting spatially homogeneous density within aligned samples. Topological solitons, like heliknotons, allow for spatial reorganization of these nanoparticles according to elastic free energy density within their micrometer-scale structures. With the nanoparticle orientations slaved to the nematic director and being switched by low voltages ≈1 V within a fraction of a second, these plasmonic composite materials are of interest for technological uses like color filters and plasmonic polarizers, as well as may lead to the development of unusual nematic phases, like pentatic liquid crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Improving Aerosol Radiative Forcing and Climate in E3SM: Impacts of New Cloud Microphysics and Improved Wet Removal Treatments

Abstract Numerous Earth system models exhibit excessive aerosol effective forcing at the top of the atmosphere (TOA), including the Department of Energy's Energy Exascale Earth System Model (E3SM). Here, in the context of the E3SM version 3 effort, the predicted particle property (P3) stratiform cloud microphysics scheme and an enhanced deep convection parameterization suite (ZM_plus) are implemented into E3SM. The ZM_plus includes a convective cloud microphysics scheme, a multi‐scale coherent structure parameterization for mesoscale convective systems, and a revised cloud base mass flux formulation considering impacts of the large‐scale environment. The P3 scheme improved cloud and radiation particularly over the Northern Hemisphere and the frequency of heavy precipitation over the tropics, and the ZM_plus improved clouds in the tropics. P3 decreases aerosol effective forcing by 0.15 W m −2 , while the ZM_plus increases it by 0.27 W m −2 , resulting from excessive direct (0.31 W m −2 ) and indirect forcing (−1.79 W m −2 ). The excessive aerosol forcings are due to aerosol overestimation associated with insufficient aerosol wet removal. By improving the physical treatments in the aerosol wet removal, we effectively mitigate anthropogenic aerosol overestimation and thus attenuate direct (0.09 W m −2 ) and indirect aerosol forcing (−1.52 W m −2 ). Adjustment to primary organic matter hygroscopicity reduces direct and indirect forcing to more reasonable values: −0.13 W m −2 and −1.31 W m −2 , respectively. On climatology, improved aerosol treatments mitigate overestimation of aerosol optical depth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Accelerating Multiscale Materials Modeling with Machine Learning

The focus of this project is to accelerate and transform the workflow of multiscale materials modeling by developing an integrated toolchain seamlessly combining DFT, SNAP, LAMMPS, (shown in Figure 1-1) and a machine-learning (ML) model that will more efficiently extract information from a smaller set of first-principles calculations. Our ML model enables us to accelerate first-principles data generation by interpolating existing high fidelity data, and extend the simulation scale by extrapolating high fidelity data (10 2 atoms) to the mesoscale (10 4 atoms). It encodes the underlying physics of atomic interactions on the microscopic scale by adapting a variety of ML techniques such as deep neural networks (DNNs), and graph neural networks (GNNs). We developed a new surrogate model for density functional theory using deep neural networks. The developed ML surrogate is demonstrated in a workflow to generate accurate band energies, total energies, and density of the 298K and 933K Aluminum systems. Furthermore, the models can be used to predict the quantities of interest for systems with more number of atoms than the training data set. We have demonstrated that the ML model can be used to compute the quantities of interest for systems with 100,000 Al atoms. When compared with 2000 Al system the new surrogate model is as accurate as DFT, but three orders of magnitude faster. We also explored optimal experimental design techniques to choose the training data and novel Graph Neural Networks to train on smaller data sets. These are promising methods that need to be explored in the future.

36 MATERIALS SCIENCE↗