SEARCH · Engineering Papers

Results for “Multiscale systems”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2021

First Principles Simulations of Electrified Interfaces in Electrochemistry

This chapter discusses some of the recent advances made in the first principles modeling of electrochemical catalysts. It also discusses the key development, namely the ability to explicitly treat the effects of surface electrification due to electrochemical processes and applied voltages in a computationally efficient manner. The chapter introduces the thermodynamics and statistical mechanics of electrified metal-solution interfaces. It then discusses the structure of the electrode-electrolyte interface and the effects of applied voltages, followed by an overview of the first-principles model and a motivating example. The chapter also provides a brief summary of classical thermodynamics and describes the basics of thermodynamic detour. It considers the thermodynamics of macroscopic systems and introduces several useful fundamental relations. Here, the chapter also considers a system to exist in a certain macrostate, which is a particular thermodynamic state specified by a set of fixed properties such as constant particle number, constant volume, and constant temperature.

36 MATERIALS SCIENCE↗

NASA NTRS2009

Multiscale Failure Analysis of Laminated Composite Panels Subjected to Blast Loading Using FEAMAC/Explicit

This preliminary report demonstrates the capabilities of the recently developed software implementation that links the Generalized Method of Cells to explicit finite element analysis by extending a previous development which tied the generalized method of cells to implicit finite elements. The multiscale framework, which uses explicit finite elements at the global-scale and the generalized method of cells at the microscale is detailed. This implementation is suitable for both dynamic mechanics problems and static problems exhibiting drastic and sudden changes in material properties, which often encounter convergence issues with commercial implicit solvers. Progressive failure analysis of stiffened and un-stiffened fiber-reinforced laminates subjected to normal blast pressure loads was performed and is used to demonstrate the capabilities of this framework. The focus of this report is to document the development of the software implementation; thus, no comparison between the results of the models and experimental data is drawn. However, the validity of the results are assessed qualitatively through the observation of failure paths, stress contours, and the distribution of system energies.

Pineda, Evan J.↗

DOE OSTI2024

United States Multi-Sector Dynamics land use and land cover base maps to support Human-Earth System Modeling

Datasets are land use and land cover (LULC) rasterized base maps at 30-m resolution for the conterminous United States (CONUS) for the years 2008, 2011, 2016, and 2019. Separate base maps are provided where LULC classifications are thematically congruent with Community Land Model (CLM), Land Use Harmonization (LUH2), and Global Change Analysis Model (GCAM), and a detailed decomposition of all combined land classes into a Multisector Dynamics (MSD) LULC product. Base maps were developed using empirically derived satellite (National Land Cover Dataset, MODIS) and combined observation datasets (Crop Data Layer, Protected Areas Database) and represent the most up-to-date accurate information on LULC in the CONUS. The four datasets encompass four different landcover classification systems: MSD Layers - The raw landcover classes obtained from reclassifying NLCD and USDA Crop data layers into a respective landcover class GCAM Layers - The MSD classes mosaiced, reclassified, and combined into the respective GCAM landcover classes CLM Layers - Similar process to GCAM layers, but mosaiced, reclassified, and combined MSD layers to their respective PFT classes LUH2 Layers - Similar process to both GCAM and CLM Layers, but mosaiced, reclassified and combined the MSD layers to align with the respective states

Food↗

DOE OSTI2025

High‐Resolution National‐Scale Water Modeling Is Enhanced by Multiscale Differentiable Physics‐Informed Machine Learning

Abstract The National Water Model (NWM) is a key tool for flood forecasting, planning, and water management. Key challenges facing the NWM include calibration and parameter regionalization when confronted with big data. We present two novel versions of high‐resolution (∼37 km 2 ) differentiable models (a type of hybrid model): one with implicit, unit‐hydrograph‐style routing and another with explicit Muskingum‐Cunge routing in the river network. The former predicts streamflow at basin outlets whereas the latter presents a discretized product that seamlessly covers rivers in the conterminous United States (CONUS). Both versions use neural networks to provide a multiscale parameterization and process‐based equations to provide a structural backbone, which were trained simultaneously (“end‐to‐end”) on 2,807 basins across the CONUS and evaluated on 4,997 basins. Both versions show great potential to elevate future NWM performance for extensively calibrated as well as ungauged sites: the median daily Nash‐Sutcliffe efficiency of all 4,997 basins is improved to around 0.68 from 0.48 of NWM3.0. As they resolve spatial heterogeneity, both versions greatly improved simulations in the western CONUS and also in the Prairie Pothole Region, a long‐standing modeling challenge. The Muskingum‐Cunge version further improved performance for basins >10,000 km 2 . Overall, our results show how neural‐network‐based parameterizations can improve NWM performance for providing operational flood predictions while maintaining interpretability and multivariate outputs. The modeling system supports the Basic Model Interface (BMI), which allows seamless integration with the next‐generation NWM. We also provide a CONUS‐scale hydrologic data set for further evaluation and use.

Song, Yalan [Civil and Environmental Engineering T↗

DOE OSTI2024

Radiation-Induced Plutonium Redox Chemistry

Plutonium plays a key role in global actinide research and nuclear fuel cycle technologies, and yet, our fundamental understanding of its inherent radiation-induced chemical behavior is limited. These radiation-induced processes cannot simply be switched off, as they are as fundamentally inherent to plutonium as the impact of relativistic effects on its f-electrons. In less chemically complex actinide systems, such as aqueous solutions of neptunium and americium, , radiolysis products play a significant role in the redox cycling of their oxidation states. However, plutonium's multiple, coexisting, and chemically active oxidation states, which comprise of bare ions and dioxo cations, provide additional redox pathways that complicate radiation-induced processes. Oxidation state control is critical for the manipulation of plutonium, especially in used nuclear fuel reprocessing technologies, where oxidation specific states are successfully extracted, and others rejected. Consequently, mechanistically understanding the behavior of plutonium’s multiple oxidation states in the presence of intense ionizing radiation fields is essential for predicting the behavior of this element under multiple conditions that support the development and innovation of nuclear fuel cycle technologies. Here, we present recent advances in our understanding of plutonium radiation chemistry, including the first-ever multiscale model for predicting gamma radiation-induced plutonium redox chemistry, and new chemical kinetics for the reaction of plutonium and its complexes of tributyl phosphate (TBP), N,N-di-(2-ethylhexyl)butyramide (DEHBA), and N,N-di-(2-ethylhexyl)isobutyramide (DEHiBA) with transients radiolysis products, a measured using electron pulse radiolysis.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2025

Development and Validation of MALAMUTE model for Electric Field Assisted Sintering of Structural Materials

Fusion power plant designs feature extreme material performance requirements for structural material candidates. In addition to conventional alloys, more advanced composites and oxide dispersion strengthened (ODS) alloys are being explored, however, achieving the desired microstructures to maximize performance using traditional manufacturing methods can be challenging. The advanced manufacturing (AM) electric field-assisted sintering (EFAS) technique offers improved control over the final microstructure through higher heating and cooling rates and moderate pressures. Modeling and simulation tools show promise in elucidating the process-structure-property-performance (PSPP) correlation for AM-produced parts, including the EFAS process. An inherently multiscale process, the EFAS technique aligns well with the multiscale modeling capability of the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE)[cite]. We present here an electro-thermo-mechanical approach to modeling the EFAS process using the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) code. Prediction of the field and gradient distributions across the EFAS tooling is required to accurately describe the conditions for the lower-scale microstructural evolution models. In this work we present the MALAMUTE model developed to predict the electrical potential, temperature, and mechanical stress distribution across the EFAS graphite tooling and part at the larger engineering-scale. Validation of the MALAMUTE engineering-scale model is completed using data from experimental densification and pre-densified runs of iron powder via EFAS at 1000oC. These runs were conducted using a Thermal Technology DCS-5 EFAS system. Data collected during the experiment runs include the direct current (DC) supplied to the graphite tooling, the temperature of the graphite tooling as measured with a pyrometer, and the force applied to the top of the graphite tooling stack, and the data were recorded every 10 seconds. Our validation approach used the current and force data from the EFAS run as boundary condition inputs to the MAMALUTE simulation; the temperature data were used to evaluate the MALAMUTE EFAS model prediction. Results of the MALAMUTE simulations are employed to connect the external pyrometer temperature measurement to the temperature profile across the part undergoing consolidation. We investigate the impact of material property variation and mesh deformation on the temperature profile as predicted by MALAMUTE. We conclude by highlighting projects where the MALAMUTE EFAS modeling and simulation capabilities will be used to assist experimental design.

36 - MATERIALS SCIENCE↗

DOE OSTI2020

DNA-assembled superconducting 3D nanoscale architectures

Studies of nanoscale superconducting structures have revealed various physical phenomena and led to the development of a wide range of applications. Most of these studies concentrated on one- and two-dimensional structures due to the lack of approaches for creation of fully engineered three-dimensional (3D) nanostructures. Here, we present a ‘bottom-up’ method to create 3D superconducting nanostructures with prescribed multiscale organization using DNA-based self-assembly methods. We assemble 3D DNA superlattices from octahedral DNA frames with incorporated nanoparticles, through connecting frames at their vertices, which result in cubic superlattices with a 48 nm unit cell. The superconductive superlattice is formed by converting a DNA superlattice first into highly-structured 3D silica scaffold, to turn it from a soft and liquid-environment dependent macromolecular construction into a solid structure, following by its coating with superconducting niobium (Nb). Through low-temperature electrical characterization we demonstrate that this process creates 3D arrays of Josephson junctions. This approach may be utilized in development of a variety of applications such as 3D Superconducting Quantum interference Devices (SQUIDs) for measurement of the magnetic field vector, highly sensitive Superconducting Quantum Interference Filters (SQIFs), and parametric amplifiers for quantum information systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

NASA NTRS2023

Foreshock Ion Motion Across Discontinuities: Formation of Foreshock Transients

In the ion foreshock, hot flow anomalies (HFAs) and foreshock bubbles (FBs) are two types of foreshock transients that have the strongest fluctuations, which can disturb the magnetosphere-ionosphere system and increase shock acceleration efficiency. They form due to interaction between the foreshock ions and solar wind discontinuities: the direction of the foreshock ion-driven current and whether it decreases or increases the magnetic field strength behind the discontinuity determine whether the transient's formation can be promoted or suppressed. Thus, to predict the HFA and FB formation and forecast their space weather effects, it is necessary to predict the foreshock ion-driven current direction. In this study, we derive analytical equations of foreshock ion velocities within discontinuities to estimate foreshock ion-driven current direction, which provides a quantitative criterion of HFA and FB formation. To validate the criterion, we use Acceleration Reconnection Turbulence & Electrodynamics of Moon's Interaction with the Sun to observe pristine solar wind discontinuities and calculate discontinuity parameters. We use Magnetospheric Multiscale to observe the foreshock ion motion around the discontinuities and show that the data support our model. This study is another step toward a predictive model of HFA and FB formation so that we can forecast their space weather effects at Earth using solar wind observations at lunar orbit or L1.

Terry Z. Liu↗

DOE OSTI2023

On a fully-implicit VMS-stabilized FE formulation for low Mach number compressible resistive MHD with application to MCF

This study presents the development and evaluation of a fully-implicit variational multiscale (VMS) stabilized unstructured finite element (FE) formulation for compressible magnetohydrodynamics (MHD) model, at low Mach number regime. The model describes the dynamics of a compressible conducting fluid in the low Mach number limit in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics phenomena. The specific applications that motivate this study are macroscopic simulations of the longer time-scale stability and disruptions of magnetic confinement fusion (MCF) devices, specifically the ITER tokamak. The discussion considers the development of the VMS FE representation, the structure of the stabilizing terms that deal with significant convective flows, the stabilization of the nearly incompressible response of the fluid flow, and the stabilization of the constraint that enforces the solenoidal involution on the magnetic field. The nonlinear discretized system is solved with scalable preconditioned Newton–Krylov iterative methods, which employs a multiphysics block preconditioning method based on approximate block factorizations and Schur complements. The study presents an evaluation of the VMS method on a 2D cartesian tearing mode instability, and illustrates the scalability of the solvers on MCF relevant problems. A set of results are also presented for longer time-scale stability and disruptions for the ITER tokamak. These include a vertical displacement event (VDE), and a (1,1) internal kink mode. Here, the formulation is demonstrated to be scalable and also reasonably robust with respect to the Lundquist number scaling.

42 ENGINEERING↗

DOE OSTI2022

Bridging Nano- and Microscale X-Ray Tomography for Battery Research by Leveraging Artificial Intelligence

X-ray computed tomography (CT) is a non-destructive imaging technique in which contrast originates from the materials' absorption coefficient. The recent development of laboratory nanoscale CT (nano-CT) systems has pushed the spatial resolution for battery material imaging to voxel sizes of 50 nm, a limit previously achievable only with synchrotron facilities. Given the non-destructive nature of CT, in situ and operando studies have emerged as powerful methods to quantify morphological parameters, such as tortuosity factor, porosity, surface area and volume expansion, during battery operation or cycling. Combined with artificial intelligence and machine learning analysis techniques, nano-CT has enabled the development of predictive models to analyse the impact of the electrode microstructure on cell performances or the influence of material heterogeneities on electrochemical responses. In this Review, we discuss the role of X-ray CT and nano-CT experimentation in the battery field, discuss the incorporation of artificial intelligence and machine learning analyses and provide a perspective on how the combination of multiscale CT imaging techniques can expand the development of predictive multiscale battery behavioural models.

diagnostics↗

DOE OSTI2020

Multi-Scale Modeling Framework for Mercury Biogeochemistry

Multi-Scale modeling of mercury (Hg) geochemical speciation and reactions has been performed by integrating atomistic quantum chemical calculations with continuum scale speciation models. Major progress has been made in the improvement of quantum chemical models to calculate critical thermodynamic data for Hg complexes in aquatic environments. Rapid and reliable quantum chemical approaches have been developed for calculating acid dissociation constants (pK a ) and stability constants (log K), with calculated mean unsigned errors of 0.5 and 1.5 log units, respectively for ligand molecules and Hg complexes. At the continuum scale, systematic analysis of uncertainty propagation in mercury (Hg) speciation modeling has been conducted and was used to identify environmental conditions under which thermodynamic constant uncertainties are significant and recommended to be accounted for. The integrated framework for multi-scale modeling of mercury geochemistry is open to the research community through the web-based multiscale modeling aqueous speciation resource, AQUA-MER. The improved quantum chemical approaches for thermodynamic constant calculations are accessible through AQUA-MER and can be used to provide the missing constants in the continuum scale speciation calculations. In addition to low molecular mass Hg complex speciation, modeling natural aquatic environments also involve the transport of high molecular weight dissolved organic matter (DOM) in reactive flows simultaneously with equilibrium and kinetic reactions. To this end, atomistic MD simulations were performed to capture the details of aggregation, mechanisms and distribution of functional groups in DOM at the molecular level. The elemental composition and calculated bulk properties of the DOM models are in close agreement with experimental measurements. A travel-time based reactive transport model in the hyporheic zone of stream corridors was established for the multicomponent Hg-DOM-S system and implemented through PFLOTRAN.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

Traditional Kinetic Monte Carlo (KMC) approaches, rooted in Gillespie’s stochastic simulation algorithm, become computationally demanding in systems with a large range of timescales. The goal of this work is to propose and study new adaptive lattice-KMC time integration strategies for spatially non-uniform systems. To that end, two novel adaptive tau-leaping methods and their corresponding time integration strategies are developed based on the idea of the “n-fold” direct KMC method. These strategies allow for the simultaneous execution of multiple reactions, advancing time by adaptively selected coarse increments. We present numerical experiments comparing the proposed methods with existing approaches in a catalytic surface kinetics application involving ammonia decomposition.

Bimolecular reactions↗

NASA NTRS2014

MMS Observatory TV Results Contamination Summary

The Magnetospheric Multiscale (MMS) mission is a constellation of 4 observatories designed to investigate the fundamental plasma physics of reconnection in the Earth's magnetosphere. The various instrument suites measure electric and magnetic fields, energetic particles, and plasma composition. Each spacecraft has undergone extensive environmental testing to prepare it for its minimum 2 year mission. In this paper, we report on the extensive thermal vacuum testing campaign. The testing was performed at the Naval Research Laboratory utilizing the "Big Blue" vacuum chamber. A total of ten thermal vacuum tests were performed, including two chamber certifications, three dry runs, and five tests of the individual MMS observatories. During the test, the observatories were enclosed in a thermal enclosure known as the "hamster cage". The enclosure allowed for a detailed thermal control of various observatory zone, but at the same time, imposed additional contamination and system performance requirements. The environment inside the enclosure and the vacuum chamber was actively monitored by several QCMs, RGA, and up to 18 ion gauges. Each spacecraft underwent a bakeout phase, which was followed by 4 thermal cycles. Unique aspects of the TV campaign included slow pump downs with a partial represses, thruster firings, Helium identification, and monitoring pressure spikes with ion gauges. Selected data from these TV tests is presented along with lessons learned.

Pressure↗

DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2009

Electrical Resistance of SiC/SiC Ceramic Matrix Composites for Damage Detection and Life-Prediction

Ceramic matrix composites (CMC) are suitable for high temperature structural applications such as turbine airfoils and hypersonic thermal protection systems due to their low density high thermal conductivity. The employment of these materials in such applications is limited by the ability to accurately monitor and predict damage evolution. Current nondestructive methods such as ultrasound, x-ray, and thermal imaging are limited in their ability to quantify small scale, transverse, in-plane, matrix cracks developed over long-time creep and fatigue conditions. CMC is a multifunctional material in which the damage is coupled with the material s electrical resistance, providing the possibility of real-time information about the damage state through monitoring of resistance. Here, resistance measurement of SiC/SiC composites under mechanical load at both room temperature monotonic and high temperature creep conditions, coupled with a modal acoustic emission technique, can relate the effects of temperature, strain, matrix cracks, fiber breaks, and oxidation to the change in electrical resistance. A multiscale model can in turn be developed for life prediction of in-service composites, based on electrical resistance methods. Results of tensile mechanical testing of SiC/SiC composites at room and high temperatures will be discussed. Data relating electrical resistivity to composite constituent content, fiber architecture, temperature, matrix crack formation, and oxidation will be explained, along with progress in modeling such properties.

Smith, Craig↗

NASA NTRS2008

A New Concurrent Multiscale Methodology for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Yamakov, Vesselin↗

NASA NTRS2016

Initial Satellite Formation Flight Results from the Magnetospheric Multiscale Mission

This paper will describe the results that have been obtained to date concerning MMS formation flying. The MMS spacecraft spin at a rate of 3.1 RPM, with spin axis roughly aligned with Ecliptic North. Several booms are used to deploy instruments: two 5 m magnetometer booms in the spin plane, two rigid booms of length 12.5 m along the positive and negative spin axes, and four flexible wire booms of length 60 m in the spin plane. Minimizing flexible motion of the wire booms requires that reorientation of the spacecraft spin axis be kept to a minimum: this is limited to attitude maneuvers to counteract the effects of gravity-gradient and apparent solar motion. Orbital maneuvers must therefore be carried out in essentially the nominal science attitude. These burns make use of a set of monopropellant hydrazine thrusters: two (of thrust 4.5 N) along the spin axis in each direction, and eight (of thrust 18 N) in the spin plane; the latter are pulsed at the spin rate to produce a net delta-v. An on-board accelerometer-based controller is used to accurately generate a commanded delta-v. Navigation makes use of a weak-signal GPS-based system: this allows signals to be received even when MMS is flying above the GPS orbits, producing a highly accurate determination of the four MMS orbits. This data is downlinked to the MMS Mission Operations Center (MOC) and used by the MOC Flight Dynamics Operations Area (FDOA) for maneuver design. These commands are then uplinked to the spacecraft and executed autonomously using the controller, with the ground monitoring the burns in real time.

satellite↗

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↗