Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale structure”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Integrating Ultra-Coarse-Grained Protein Models into Accessible Workflows for Multiscale Molecular Dynamics

To capture protein conformational transitions using molecular dynamics (MD), several simulation resolutions covering different spatial and temporal scales are typically needed. All-atom (AA) simulations provide fine resolution, but are computationally infeasible for large systems over longer durations. Coarse-grained (CG) and ultra-coarse-grained (UCG) models have a lower resolution and computational cost while still being able to conserve essential protein features. Prior work on a Multiscale Machinelearned Modeling Infrastructure (MuMMI) combined both AA and CG simulations to study RAS-RAF protein interactions, leveraging CG models for longer time scales and using AA to investigate unusual conformations in greater detail. However, MuMMI is still resource-intensive, and this study aims to maximize exploration of the protein conformational space while reducing computational cost. In this paper, we build on prior work that integrates UCG models based on heterogeneous elastic network modeling (hENM) into the MuMMI workflow. We demonstrate that UCG models enable accurate sampling of protein conformations, focusing on simulating RAS-RAF protein interactions. Using higher-resolution CG Martini simulation data, we can automatically refine intramolecular interactions in UCG models. We present a scalable Python package that uses fluctuations observed in higher-resolution CG Martini simulations to estimate bond coefficients of the UCG model. We built novel machine learning-based backmapping methods to recover more detailed CG Martini structures from UCG structures, using diffusion models to learn the mapping between scales. Finally, we present UCG-mini-MuMMI, an accessible and less compute-intensive version of MuMMI as a resource for the scientific community. Incorporating UCG models into MD studies is applicable to a broad range of systems and proteins, and our study offers insights into the advantages and limitations of these methods.

Chemical structure↗
DOE OSTI2023

Electrodeposited nickel coatings for exceptional corrosion mitigation in industrial grade molten chloride salts for concentrating solar power

Molten chloride salt eutectics are attractive candidates for use as thermal energy storage media and heat transfer fluids in generation-three concentrating solar thermal power (Gen3 CSP) plants. However, corrosion of alloys in molten chloride salts, especially at high temperatures, is an extremely challenging problem that studies focus on lower temperatures, shorter durations, or analytical grade, and high-purity, salts. To date, there has been no study on corrosion or corrosion mitigation in an industrial-grade salt at a high temperature such as 750 °C. To alleviate this knowledge gap, the study presents new multiscale fractal-textured Ni coatings on various alloy surfaces for effective corrosion mitigation at 750 °C in molten chloride salts. Using the electrodeposition method, durable double-layer textured coatings were formed on stainless steel alloys (SS316, SS310, and SS347) and In800H. The corrosion performance of the coatings is investigated in both analytical-grade purity and, for the first time, practically relevant industrial-grade chloride salts. Ni-coated ferrous alloys showed an exceptionally reduced corrosion rate in the range of 350–480 μm/y in analytical-grade salts, and between 450 and 490 μm/y in purified industrial-grade salts at 750 °C. Ni coatings on ferrous alloys reduced corrosion rates by as much as 70% compared to uncoated surfaces and were comparable to the expensive Ha230 alloy with a high Ni content. As a result, by the use of innovative fractal corrosion mitigation coatings, for the first time, low-cost structural alloys are rendered viable for use with industrial-grade chloride salts, which is profoundly beneficial in practical systems.

14 SOLAR ENERGY↗
DOE OSTI2025

Convergence and Quantum Advantage of Trotterized MERA for Strongly-Correlated Systems

Strongly-correlated quantum many-body systems are difficult to study and simulate classically. We recently proposed a variational quantum eigensolver (VQE) based on the multiscale entanglement renormalization ansatz (MERA) with tensors constrained to certain Trotter circuits. Here, we determine the scaling of computation costs for various critical spin chains which substantiates a polynomial quantum advantage in comparison to classical MERA simulations based on exact energy gradients or variational Monte Carlo. Algorithmic phase diagrams suggest an even greater separation for higher-dimensional systems. Hence, the Trotterized MERA VQE is a promising route for the efficient investigation of strongly-correlated quantum many-body systems on quantum computers. Furthermore, we show how the convergence can be substantially improved by building up the MERA layer by layer in the initialization stage and by scanning through the phase diagram during optimization. For the Trotter circuits being composed of single-qubit and two-qubit rotations, it is experimentally advantageous to have small rotation angles. We find that the average angle amplitude can be reduced considerably with negligible effect on the energy accuracy. Benchmark simulations suggest that the structure of the Trotter circuits for the TMERA tensors is not decisive; in particular, brick-wall circuits and parallel random-pair circuits yield very similar energy accuracies.

Miao, Qiang [Duke Quantum Center, Duke University,↗
DOE OSTI2021

Distribution of Bound and Free Water in Anatomical Fractions of Pine Residues and Corn Stover as a Function of Biological Degradation

Biomass quality is influenced by water’s abundance, distribution, and status in relation to other chemical species within the polymer matrix. Water interacts with polymers that make up the cell walls, and these interactions govern the physical and chemical changes that occur during the storage and preprocessing of biomass feedstocks. Time-domain nuclear magnetic resonance (TD-NMR) was employed to explore variations in the physical constraints of water within the lignocellulosic microstructure in distinct anatomical fractions of biomass and as a function of biological degradation. The Carr–Purcell–Meiboom–Gill sequence, when combined with knowledge of the chemical composition and physical structure of pine residues and corn stover anatomical fractions, gives an accurate measurement of the bound and free water. In this work, the impacts of storage and biological degradation were investigated to elucidate changes in the status and distribution of water within distinct plant tissues. We also investigate how degradation during storage affects water interactions in different pine residues (e.g., bark, branch, and needle) and corn stover (e.g., cob, leaf, and stalk) anatomical fractions using transverse relaxation times (T2). As demonstrated herein, TD-NMR provides quantitative data on lignocellulosic biomass–water interactions within anatomical fractions, which can further aid in the investigation of preprocessing effects on feedstock quality. Our findings suggest that biological heating enhances biomass–water interactions at the cellular and macromolecular scale. In addition, analysis of three-dimensional scanning electron microscopy reconstructions indicates that surface roughness wavelengths align with microscale roughness, suggesting that pine forestry residue and corn stover particles have primarily hydrophobic exterior surfaces. Furthermore, this study offers multiscale insights into understanding the microstructure, wettability, and chemical environment that dictate diffusion, enzyme access, and recalcitrance of lignocellulosic biomass.

09 BIOMASS FUELS↗
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 OSTI2023

Study on the Carbon Nanostructures for Nanosized Si Electrodes

Silicon (Si) has been considered as a next-generation anode material due to natural abundance, low operating potential (<0.5 V vs. Li/Li+), and high theoretical specific capacity of 4200 mAh g-1.[1] However, the electrochemical alloying reaction of Si involves large volume changes of 400% during lithiation and delithiation, causing cracking and pulverization of Si.[1] In addition, solid electrolyte interface (SEI) of Si anode experiences constant changes due to unstable SEI reactivity.[2] Considerable efforts have been made to design nanostructured Si materials to address the issues because nanostructuring can relieve the mechanical strain.[3] However, when compared to micrometer-sized active materials having sufficient electrolyte pathway through interparticle space, the nanoparticle (NP)-based electrodes tend to be densely packed. As mass loading is higher and electrode is thicker, ion transport issue can become more severe in the densely packed electrodes.[4] In this work, we will present studies on engineered porosity impacts performance in Si NP electrodes by altering carbon nanostructures. We explore the extreme limit of effectively non-porous electrodes using quasi-spherical, Si NPs, which result in densely packed electrodes when slurry is prepared with conventional carbon black. We engineer porosity using carbon nanostructures including multi-walled carbon nanotubes and carbon nanorods in place of conventional carbon black to create pore structure in Si NP-based electrodes. These experiments provide a correlation between mass loading, porosity and silicon utilization in Si NP-based electrodes. References: [1] C. K. Chan, H. Peng, G. Liu, K. Mcllwrath, X. F. Zhang, R. A. Huggins, Y. Cui, "High-performance lithium battery anodes using silicon nanowires", Nat. Nanotechnol., 2008, 3, 31, [2] J. D. McBrayer, M.-T. F. Rodrigues, M. C. Schulze, D. P. Abraham, C. A. Apblett, I. Bloom, G. M. Carroll, A. M. Coclasure, C. Fang, K. L. Harrison, G. Liu, S. D. Minteer, N. R. Neale, G. M. Veith, C. S. Johnson, J. T. Vaughey, A. K. Burrel, B. Cunningham, "Calendar aging of silicon-containing batteries", Nat. Energy, 2021, 6, 1164, [3] H. Wu, G. Zheng, N. Liu, T. J. Carney, Y. Yang, Y. Cui, "Engineering Empty space between Si Nanoparticles for Lithium-Ion Battery Anodes", Nano Lett., 2012, 12, 904, [4] X. Zhang, Z. Ju, Y. Zhu, K. J. Takeuchi, E. S. Takeuchi, A. C. Marschilok, G. Yu, "Multiscale Understanding and Architecture Design of High Energy/Power Lithium-Ion Battery Electrodes", Adv. Energy Mater., 2021, 11, 2000808

carbon nanostructure↗
DOE OSTI2024

Initial Development of Fusion Magnet Simulation Capabilities for Performance and Safety Evaluation Using the MOOSE Framework

Fusion energy holds the promise of being a transformative technology as a carbon-neutral, sustainable source of energy. Whole device modeling and the development of fusion digital twins will be increasingly important for emerging fusion device concepts at both national laboratories and within the commercial fusion industry. However, meeting the challenge of whole device modeling of fusion energy devices requires robust, multiphysics, multiscale modeling and simulation technologies capable of running on large-scale supercomputers. Detailed analysis of individual systems at-scale is also required to ensure safe and efficient operation as well as provide the safety basis for future device designs and licensing activities. In a tokamak, toroidal and poloidal magnets confine and shape the fusion plasma to promote the fusion reaction. High plasma temperatures and high magnetic field requirements in modern design concepts (leading to high amounts of energy stored within each magnet) impose electrical, thermal, and mechanical loads on the magnet components, which in turn impacts the safety considerations of the magnet and their supporting systems. Idaho National Laboratory (INL) has a history of working in this space, including development and benchmarking of the Magnetic System Circuitry Analysis Program (MSCAP) and Magnet Arcing (MAGARC) codes to study magnet quench events; notably, MAGARC was used to study quenching during the ITER Engineering Design Activity. However, these legacy codes and capabilities are not parallel and scalable, and new tools are required for future advances in this area, which leads to the INL-developed Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Developed originally for fission reactor systems under United States Department of Energy, Office of Nuclear Energy modeling and simulation programs, the MOOSE framework has been well-suited to multiscale, multiphysics modeling and simulation needs for nuclear systems. The framework is open-source, well-tested, under continuous development and deployment, and developed to a Nuclear Quality Assurance, Level 1 software quality standard. MOOSE has also been used in the fusion space previously in several projects: INL’s Tritium Migration Analysis Program, Version 8 (TMAP8) for tritium migration, UK Atomic Energy Authority’s A Unified Resource for OpenMC (fusion) Reactor Applications (AURORA) code for fusion thermo-mechanical and neutronics analysis, and Argonne National Laboratory’s Cardinal for high-fidelity computational fluid dynamics and neutronics. However, to model superconducting magnets, several MOOSE enhancements are required: additions to the current MOOSE electromagnetic capabilities, new material libraries for superconductors of interest (such as YBCO), as well as fusion-specific models for thermo-mechanics. This talk will discuss initial development activities to build these capabilities in MOOSE, focusing on initial validation and benchmarking activities. Proposed coupling workflows and future work to support the simulation of fusion magnets and magnet structural assemblies for performance and safety evaluation in MOOSE will also be discussed.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Ion Size Effects on the Thermodynamic, Kinetic, and Mechanical Properties during Ion Exchange in Solid-State Electrolytes

Ion exchange offers a pathway to impose residual compressive stresses to mitigate the electro-chemo-mechanical cracking of solid-state electrolytes such as lithium lanthanum zirconium oxide. This study uses a coupled multiscale framework (integrating density functional theory (DFT), molecular dynamics (MD), and continuum modeling) to examine how exchange ion size influences stress, diffusion, fracture toughness, and electronic properties. Larger isovalent ions (Na + , Ag + , K + ) were exchanged with Li + , with DFT confirming their preference for octahedral 96h sites and a linear relationship between ion size and chemical free expansion coefficient. MD simulations reveal stress and concentration effects on exchange ion diffusivity at elevated temperatures, with Na + and Ag + maintaining favorable mobility while K + showing concentration-dependent clustering. Continuum modeling predicts the range of fracture strength improvements and the required ion exchange concentration profile. It was shown that a 5% surface exchange concentration can induce ∼0.6 GPa of surface compressive stress using Na + and ∼1.0 GPa of surface compressive stress using Ag + . On the other hand, larger ion exchange species may penalize Li + diffusivity by increasing the activation volume and activation energy. Interestingly, Na + has a negligible penalty on Li-ion diffusivity. The room temperature Li + ion diffusivity is reduced by ∼40% with Ag + ion exchange. Electronic band structure analysis shows no size-dependent change in the bandgap, though Ag + introduces localized defect states near the valence band maximum. This study highlights ion size as a key factor in optimizing LLZO properties, offering a framework to improve the solid-state battery performance.

Jagad, Harsh D. [Brown Univ., Providence, RI (Unit↗
DOE OSTI2022

The Helicity Sign of Flux Transfer Event Flux Ropes and Its Relationship to the Guide Field and Hall Physics in Magnetic Reconnection at the Magnetopause

Flux Transfer Events (FTEs) are transient magnetic flux ropes typically found at the Earth's magnetopause on the dayside. While it is known that FTEs are generated by magnetic reconnection, it remains unclear how the details of magnetic reconnection controls their properties. A recent study showed that the helicity sign of FTEs positively correlates with the east-west (B y ) component of the Interplanetary Magnetic Field (IMF). With data from the Cluster and Magnetospheric Multiscale missions, we performed a statistical study of 166 quasi force-free FTEs. We focus on their helicity sign and possible association with upstream solar wind conditions and local magnetic reconnection properties. Using both in situ data and magnetic shear modeling, we find that FTEs whose helicity sign corresponds to the IMF B y are associated with moderate magnetic shears while those that do not correspond to the IMF B y are associated with higher magnetic shears. While uncertainty in IMF propagation to the magnetopause may lead to randomness in the determination of the flux rope core field and helicity, we rather propose that for small IMF B y , which corresponds to high shear and low guide field, the Hall pattern of magnetic reconnection determines the FTE core field and helicity sign. In that context we explain how the temporal sequence of multiple X-line formation and the reconnection rate are important in determining the flux rope helicity sign. This work highlights a fundamental connection between kinetic processes at work in magnetic reconnection and the macroscale structure of FTEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Multimodal Nanoscale Mapping of Local Structure and CO 2 Adsorption in Metal–Organic Frameworks

Diamine functionalization of the metal−organic framework Mg 2 (dobpdc) (dobpdc 4− = 4,4′-dioxidobiphenyl-3,3′-dicarboxylate) significantly enhances its selectivity for CO 2 capture from flue gases and air. The structure and CO 2 capacity of such materials are typically assessed using bulk techniques that rely on averaging signal over large ensembles of unit cells, obscuring local heterogeneities, such as variations in CO 2 occupancy across individual nanocrystals. To resolve this limitation, we demonstrate a multimodal, nanoscale characterization of Mg 2 (dobpdc) appended with 1,3-diaminopropane. By employing recently developed characterization techniques at progressively smaller length scales, we uncover insights from correspondingly smaller populations of unit cells. First, we use parallel-beam 3D electron diffraction (3D ED) to identify a prominent expansion in lattice parameters upon desorption of CO 2 , as observed at the level of single nanocrystals. Second, we use convergent-probe 4D scanning transmission electron microscopy (4D-STEM) to quantify associated differences in lattice strain as a function of gas loading and diamine appending. These measurements sample small subvolumes within individual nanocrystals. Finally, we apply infrared scattering scanning near-field optical microscopy (IR s- SNOM) to confirm variable CO 2 chemisorption across adsorption sites at the surface of single nanocrystals. This multimodal, multiscale approach allows us to map heterogeneity within individual nanocrystals. Collectively, these findings emphasize the importance of local, nanoscale characterization of metal−organic frameworks in revealing previously unresolvable features that impact their performance.

Karstens, Sarah L. [University of California, Berk↗
DOE OSTI2022

Lagrangian conditional statistics and flow topology in edge plasma turbulence

Lagrangian statistics and particle transport in edge plasma turbulence are investigated using the Hasegawa–Wakatani model and its modified version. The latter shows the emergence of pronounced zonal flows. Different values of the adiabaticity parameter are considered. The main goal is to characterize the role of coherent structures, i.e., vortices and zonal flows, and their impact on the Lagrangian statistics of particles. Computationally intensive long time simulations following ensembles of test particles over hundreds of eddy turnover times are considered in statistically stationary turbulent flows. The flow topology is characterized using the Lagrangian Okubo–Weiss criterion in order to split the flow into topologically different domains. In elliptic and hyperbolic regions, the probability density functions (PDFs) of the residence time have self-similar algebraic decaying tails. However, in the intermediate regions, the PDFs exhibit exponentially decaying tails. Topologically conditioned PDFs of the Lagrangian velocity, and acceleration and density fluctuations are likewise computed. The differences between the classical Hasegawa–Wakatani system and its modified version are assessed, and the role of zonal flows is highlighted. Here, the density flux spectrum, which characterizes the contributions of different length scales, is studied, and its inertial scaling is found to be in agreement with predictions based on dimensional arguments. Analyzing the angular change of particle tracers at different time scales, corresponding to coarse grained curvature, completes the study, and these multiscale geometric statistics quantify the directional properties of the particle motion in different flow regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Particle Acceleration by Magnetic Reconnection in Geospace

Particles are accelerated to very high, non-thermal energies during explosive energy-release phenomena in space, solar, and astrophysical plasma environments. While it has been established that magnetic reconnection plays an important role in the dynamics of Earth’s magnetosphere, it remains unclear how magnetic reconnection can further explain particle acceleration to non-thermal energies. Here we review recent progress in our understanding of particle acceleration by magnetic reconnection in Earth’s magnetosphere. With improved resolutions, recent spacecraft missions have enabled detailed studies of particle acceleration at various structures such as the diffusion region, separatrix, jets, magnetic islands (flux ropes), and dipolarization front. With the guiding-center approximation of particle motion, many studies have discussed the relative importance of the parallel electric field as well as the Fermi and betatron effects. However, in order to fully understand the particle acceleration mechanism and further compare with particle acceleration in solar and astrophysical plasma environments, there is a need for further investigation of, for example, energy partition and the precise role of turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Revealing 3D Morphological and Chemical Evolution Mechanisms of Metals in Molten Salt by Multimodal Microscopy

Growing interest in molten salts as effective high-temperature heat-transfer fluids for sustainable energy systems drives a critical need to fundamentally understand the interactions between metals and molten salts. This work utilizes the multimodal microscopy methods of synchrotron X-ray nanotomography and electron microscopy to investigate the 3D morphological and chemical evolution of two-model systems, pure nickel metal and Ni-20Cr binary alloy, in a representative molten salt (KCl-MgCl2 50–50 mol %, 800 °C). In both systems, unexpected shell-like structures formed because of the presence of more noble tungsten, suggesting a potential route of using Ni–W alloys for enhanced molten-salt corrosion resistance. The binary alloy Ni-20Cr developed a bicontinuous porous structure, reassembling functional porous metals manufactured by dealloying. This work elucidates better mechanistic understanding of corrosion in molten salts, which can contribute to the design of more reliable alloys for molten salt applications including next-generation nuclear and solar power plants and opens the possibility of using molten salts to fabricate functional porous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Enhanced Planar Antenna Efficiency Through Magnetic Thin-Films

This work proposes to use magnetic material as the substrate of planar antennas to overcome the platform effect caused by the conducting ground plane. The upper bound of the radiation efficiency of an electric-current-driven low-profile antenna is theoretically derived, which is inversely proportional to the Gilbert damping factor of the magnetic material. Meanwhile, the improvement of radiation due to the use of magnetic material is demonstrated by a three-dimensional (3D) multiphysics and multiscale time-domain model. The simulation results match the theoretical derivation, showing 25% radiation efficiency from a planar antenna backed by a FeGaB thin film with 2.56 m thickness. Furthermore, for conductive ferromagnetic materials, it is shown that the eddy current loss can be well suppressed by laminating the thin film into multiple layers. Here, the radiation efficiency of the modeled antenna with a conductive ferromagnetic substrate is improved from 2.2% to 11.8% by dividing the substrate into 10 layers, with a ferromagnetic material fill factor of 93%.

42 ENGINEERING↗
DOE OSTI2024

Investigating Mechanisms Driving Differences in the Characteristics of Precipitation in the E3SM Multiscale Modeling Framework With 2D Versus 3D Cloud Resolving Model Configurations

Abstract In this study, we compare the Energy Exascale Earth Systems Model (E3SM) multiscale modeling framework (MMF) with the cloud resolving model (CRM) configured in two (2dMMF) and three (3dMMF) dimensions. We explore how CRM dimensionality impacts the representation of mean and extreme precipitation characteristics. Our results show that tropical mean precipitation patterns are better represented in 3dMMF compared to observations (Integrated Multi‐satellitE Retrivals for GPM and Global Precipitation Climatology Project One Degree Daily products), while 2dMMF better captures extreme precipitation intensity, with systematic land‐ocean differences in precipitation and cloud‐associated variables. These differences are attributed to the co‐occurrence of CRM throttling (i.e., suppressed convection in due to smaller numbers of CRM columns and domain size) and dilution (i.e., 3‐D cloud circulations with increased entrainment and lower precipitation efficiency) effects. Overall, throttling results in more low‐level humidity in 2dMMF and dilution contributes to more high clouds with less precipitation efficiency in 3dMMF. Since throttling occurs more strongly over the ocean than land, the 3dMMF tends to have less cloud liquid and precipitation over the ocean and more cloud ice and precipitation over land. These results may serve as a guide for choosing the CRM structure to reduce precipitation and cloud‐related biases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Development of Machine-Learned Interatomic Potentials to Predict Structure, Transport, and Reactivity in Platinum-Based Fuel Cells

Machine-learned interatomic potentials (MLIPs) have rapidly progressed in accuracy, speed, and data efficiency in recent years. However, training robust MLIPs in multicomponent systems remains a challenge. In this work, we train an MLIP to describe hydrated Nafion ionomers and platinum catalysts, which are important components of fuel cells, by constructing a diverse training set to describe the bulk polymer and interfacial catalyst–polymer interactions well. We use our trained MLIP to study the properties of the platinum–Nafion system, including polymer structure, proton mobility in a bulk Nafion polymer and near a platinum-Nafion interface, and reactions near and far from the interface, finding excellent results for structure and reactions contained within our training set. Transport seems to be well described, with both vehicular transport and Grotthuss hopping captured, although converged calculations of diffusivities were not computed because they require calculations of tens of nanoseconds that are challenging with current state-of-the-art MLIPs. The combined insights that this model provides can be leveraged to optimize fuel cell performance, and the approach can be applied to other chemical processes and devices where structure, transport, and reactivity all contribute to the overall observed performance.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Multiscale hierarchical and heterogeneous mechanical response of additively manufactured novel Al alloy investigated by high-resolution nanoindentation mapping

Smart alloying and microstructural engineering mitigate challenges associated with laser-powder bed fusion additive manufacturing (L-PBFAM). A novel Al–Ni–Ti–Zr alloy utilized grain refinement by heterogeneous nucleation and eutectic solidification to achieve superior performance-printability synergy. Conventional mechanical testing cannot delineate complex micromechanics of such alloys. This study combined multiscale nanomechanical and microstructural mapping to illustrate mechanical signatures associated with hierarchical heat distribution and rapid solidification of L-PBFAM. The disproportionate hardening effect imparted by Al 3 (Ti,Zr) precipitates in the pool boundaries and the semi-solid zone was successfully demonstrated. Nanomechanical response associated with heterogeneity in particle volume fraction and coherency across melt pool was interpreted from nanoindentation force–displacement curves. The hardness map effectively delineated the weakest and strongest sections in the pool with microscopic accuracy. The presented approach serves as a high throughput methodology to establish the chemistry-processing-microstructure-properties correlation of newly designed alloys for L-PBFAM.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic modeling of LiF microstructure ionic conductivity and its influence on nucleation and plating

We report the formation and degradation of the solid electrolyte interphase (SEI) and its underlying transport properties play an essential role in the overall performance of lithium-ion batteries. This paper presents classical molecular dynamics studies on polycrystalline inorganic lithium fluoride (LiF) layers to model and predict the SEI transport properties. The ionic conductivity is obtained from the lithium-ion diffusivity in polycrystalline structures of LiF using the Nernst-Einstein relation. The predicted molecular dynamics data are used in a continuum scale phase-field model to evaluate the plating kinetics under fast charging conditions. The analysis emphasizes that the SEI ionic conductivity properties impact the plating dynamics, where SEI's low ion conductivity value is prone to large plating and subsequent capacity degradation. The combination of atomic and continuum scale studies shown herein lays a foundation to tune in SEI transport properties to decrease the amount of lithium plating and improve the performance of fast-charging batteries.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗