SEARCH · Search NASA

Results for “Intrinsic computation”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Higher-order phonon scattering: advancing the quantum theory of phonon linewidth, thermal conductivity and thermal radiative properties

Phonon scattering plays a central role in the quantum theory of phonon linewidth, which in turn governs important properties including infrared spectra, Raman spectra, lattice thermal conductivity, thermal radiative properties, and also signifi- cantly affects other important processes such as hot electron relaxation. Since Maradudin and Fein’s classic work in 1962, three-phonon scattering had been considered as the dominant intrinsic phonon scattering mechanism and has seen tremendous advances. However, the role of the higher-order four-phonon scattering had been persistently unclear and so was ignored. The tremendous complexity of the formalism and computational challenges stood in the way, prohibiting the direct and quantitative treatment of four-phonon scattering. In 2016, a rigorous four-phonon scattering formalism was developed, and the prediction was realized using empirical potentials. In 2017, the method was extended using first-principles calculated force constants, and the thermal conductivities of boron arsenides (BAs), Si and diamond were predicted. The predictions for BAs were later confirmed by several independent experiments. Four-phonon scattering has since been investigated in a range of materials and established as an important intrinsic scattering mechanism for thermal transport and radiative properties. Specifically, four-phonon scattering is important when the fourth-order scattering potential or phase space becomes relatively large. The former scenario includes: (i) nearly all materials when the temperature is high; (ii) strongly anharmonic (low thermal conductivity) materials, including most rocksalt compounds, halides, hydrides, chalcogenides and oxides. The latter scenario includes: (iii) materials with large acoustic–optical phonon band gaps, such as XY compounds with a large atomic mass ratio between X and Y; (iv) two- dimensional materials with reflection symmetry, such as single-layer graphene, single-layer boron nitride and carbon nanotubes; and (v) phonons with a large density of states, such as optical phonons, which are important for Raman, infrared and thermal radiative properties. Four-phonon scattering is expected to gain broad interest in various technologically important materials for thermoelectrics, thermal barrier coatings, thermal energy storage, phase change, nuclear power, ultra-high temperature ceramics, infrared spectra, Raman spectra, radiative transport, hot electron relaxation and radiative cooling. Four-phonon scattering has been, and will continue to be, established as an important intrinsic phonon scattering mechanism beyond three-phonon scattering. The prediction of four-phonon scattering will transition from a breakthrough to a new routine in the next decade.

Feng, Tianli↗

DOE OSTI2026

Selective Chemical Looping Combustion of Terminal Alkynes in Mixtures with Alkenes

The selective combustion of terminal alkynes in mixtures with alkenes is demonstrated during anaerobic reduction half-cycles on bulk bismuth oxide (Bi 2 O 3 ) as an approach to remove alkynes, which act as inhibitors in olefin polymerization. Bi 2 O 3 combusts phenylacetylene in styrene, 3-methylphenylacetylene in 3-methylstyrene, propyne in propylene, 1-hexyne in 1-hexene, and 1-octyne in 1-octene, with alkyne combustion selectivities exceeding 96%. Near unity reaction orders for hydrocarbon consumption during reduction half-cycles are consistent with combustion pathways initiated by rate-determining initial C–H activation, which drive selective alkyne combustion through intrinsic differences in the first-order rate constants for alkyne and alkene combustion rather than preferential adsorption of alkynes on Bi 2 O 3 surfaces. Computational assessments of initial C–H activation pathways for alkynes and alkenes on (010) α-Bi 2 O 3 surfaces using density functional theory illustrate that heterolytic transition states which form proton-carbanion pairs on Bi–O sites kinetically favor the activation of alkynes rather than alkenes due to differences in C–H bond acidity, and the barrier for heterolytic C–H activation is dictated in part by the sum of the molecular deprotonation energy and the energy to bind an R – carbanion to a Bi site in its transition-state geometry. Finally, these heterolytic reactivity channels during selective chemical looping combustion present novel routes for purifying olefin gas streams containing alkyne impurities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2005

GOES Satellite Data Validation Via Hand-held 4 LED Sun Photometer at Norfolk State University

Sun photometry is a passive means of measuring a quantity of light radiation. The GIFTS- IOMI/GLOBE Water Vapor/Haze Sun photometer contains four light emitting diodes (LEDs), which are used to convert photocurrent to voltage. The intensity of the incoming and outgoing radiation as detected on the Earth s surface can be affected by aerosols and gases in the atmosphere. The focus of this research is primarily on aerosol and water vapor particles that absorb and reemit energy. Two LEDs in the photometer correspond to light scattered at 530 nm (green spectrum) and 620 nm (red spectrum). They collect data pertaining to aerosols that scatter light. The other two LEDs detect the light scattered by water vapor at wavelengths of 820 nm and 920 nm. The water vapor measurements will be compared to data collected by the Geostationary Observation Environmental Satellite (GOES). Before a comparison can be made, the extraterrestrial constant (ET), which is intrinsic to each sun photometer, must be measured. This paper will present determination of the ET constant, from which the aerosol optical thickness (AOT) can be computed for comparison to the GOES satellite to ascertain the reliability of the sun photometer.

Reynolds, Arthur, Jr.↗

DOE OSTI2025

DS-TIDE: Harnessing Dynamical Systems for Efficient Time-Independent Differential Equation Solving

Time-Independent Differential Equations (TIDEs) are central to modeling equilibrium behavior across a wide range of scientific and engineering domains, from electrostatics to porous media flow. Conventional numerical solvers offer reliable solutions but incur significant computational costs due to fine-grained discretization and iterative procedures. Machine learning-based approaches address this by replacing iterative solving processes with one-time inference; however, their sophisticated models require extensive training resources that often exceed those of traditional solvers. Consequently, designing a TIDE solver that achieves high accuracy, broad applicability, and exceptional computational efficiency remains a fundamental challenge. In this paper, we propose DS-TIDE, a novel hardware solver that is inspired by, and subsequently leverages, the intrinsic connection between Dynamical Systems (DS) and Differential Equations (DEs) to efficiently and accurately solve TIDEs. DS-TIDE employs a CMOS-compatible DS-based processor, whose physical states evolve under carefully designed DE-driven dynamics and naturally converge to equilibrium -- the solution of the target TIDE -- within ~1µs on a ~1-watt DS-TIDE processor. To enhance expressivity, DS-TIDE incorporates Heterogeneous Dynamics with Temporal Layering (HDTL), which solves TIDEs through a three-stage DS evolution -- conditioning, solving, and decoding -- each governed by specialized dynamics. The entire evolution process is analogous to an infinitely deep neural network temporally unrolled, offering the system the capability of representing complex equations. Furthermore, DS-TIDE is equipped with an on-device DS-DE Auto-Alignment mechanism that dynamically adapts intrinsic hardware dynamics within milliseconds, effectively aligning the system’s dynamics to diverse target DEs. Experimental results across TIDEs from a wide range of scientific and engineering domains demonstrate that DS-TIDE achieves ~10^3× speedup, ~10^5× energy savings, and competitive or superior accuracy compared to state-of-the-art numerical and ML-based solvers.

Liu, Chuan↗

DOE OSTI2024

Physics-Guided, Physics-Informed, and Physics-Encoded Neural Networks and Operators in Scientific Computing: Fluid and Solid Mechanics

Abstract Advancements in computing power have recently made it possible to utilize machine learning and deep learning to push scientific computing forward in a range of disciplines, such as fluid mechanics, solid mechanics, materials science, etc. The incorporation of neural networks is particularly crucial in this hybridization process. Due to their intrinsic architecture, conventional neural networks cannot be successfully trained and scoped when data are sparse, which is the case in many scientific and engineering domains. Nonetheless, neural networks provide a solid foundation to respect physics-driven or knowledge-based constraints during training. Generally speaking, there are three distinct neural network frameworks to enforce the underlying physics: (i) physics-guided neural networks (PgNNs), (ii) physics-informed neural networks (PiNNs), and (iii) physics-encoded neural networks (PeNNs). These methods provide distinct advantages for accelerating the numerical modeling of complex multiscale multiphysics phenomena. In addition, the recent developments in neural operators (NOs) add another dimension to these new simulation paradigms, especially when the real-time prediction of complex multiphysics systems is required. All these models also come with their own unique drawbacks and limitations that call for further fundamental research. This study aims to present a review of the four neural network frameworks (i.e., PgNNs, PiNNs, PeNNs, and NOs) used in scientific computing research. The state-of-the-art architectures and their applications are reviewed, limitations are discussed, and future research opportunities are presented in terms of improving algorithms, considering causalities, expanding applications, and coupling scientific and deep learning solvers.

Computer Science↗

NASA NTRS1994

A Fortran-90 Based Multiprecision System

The author has developed a new version of his Fortran multiprecision computation system that is based on the Fortran-90 language. With this new approach, a translator program is not required - translation of Fortran code for multiprecision is accomplished by merely utilizing advanced features of Fortran-90, such as derived data types and operator extensions. This approach results in more reliable translation and also permits programmers of multiprecision applications to utilize the full power of the Fortran-90 language. Three multiprecision datatypes are supported in this system: multiprecision integer. real and complex. All the usual Fortran conventions for mixed mode operations are supported, and many of the Fortran intrinsics, such as SIN, EXP and MOD, are supported with multiprecision arguments. This paper also briefly describes an interesting application of this software, wherein new number-theoretic identities have been discovered by means of multiprecision computations.

Bailey, David H.↗

DOE OSTI2022

Computational determination of a primary diffusion mode in γ U-10Mo under irradiation

Low enriched uranium (< 20 % 235 U)-molybdenum (U-Mo) monolithic fuel is the primary candidate for high performance research and test reactors, and is in the process of being qualified to replace highly enriched uranium (≥ 20 % 23 5U) fuel. As part of the qualification process, it is critical to understand and predict the behavior of fission gas bubbles under irradiation, which affects fuel swelling and fuel failure. Mechanistic fuel models are being developed that can both reproduce the existing experimental data for fuel swelling, and be further applied to irradiation conditions beyond the experimental scope. Diffusion of species under irradiation conditions is an important parameter in the mechanistic fuel models; however, no temperature-relevant experimental diffusion data exists. In the present work, radiation-enhanced diffusion coefficients of U, Mo, and Xe in λU-10wt.%Mo were calculated in the temperature range between 300 K and 1400 K via rate-theory models and molecular dynamics simulations with an embedded-atom method interatomic potential for the U-Mo-Xe system. Accordingly, total diffusion coefficients under relevant irradiation conditions are determined using previously obtained intrinsic thermal diffusion and radiation-driven diffusion coefficients, as well as the newly calculated radiation-enhanced diffusion coefficients presented herein. Radiation-enhanced diffusion of U and Mo was dominant in the intermediate temperature range, whereas radiation-enhanced diffusion of Xe did not significantly contribute to total diffusion of Xe at the relevant fission rate densities. Radiation-enhanced diffusion of Xe became faster than both intrinsic thermal diffusion and radiation-driven diffusion at a fission rate density of 5 x 10 22 fissions/m 3 /s, which is higher than the typical fission rate density range in research reactors. The temperature regime that radiation-enhanced diffusion of each element dominated was dependent on the fission rate density. Finally, the total diffusion coefficients of U, Mo, and Xe, updated in this work, will be utilized as parameters in the mechanistic fuel models to help predict the behavior of fission gas bubbles under irradiation more accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

The critical role of intrinsic defects and many-body interactions on the stability of MnBi2Te4

Intrinsic antisite defects pose a major challenge to understanding and predicting the exotic properties of the layered topological magnetic insulator MnBi2Te4 (MBT). In this work, we study the origin of the abundance of intrinsic defects in MBT, including many-body defect–defect interactions and many-body electronic correlations. Until now, ab initio methods have struggled to explain thermodynamic stability and properties influenced by defect behavior in MBT. We model native Mn–Bi antisite defects in MBT at finite temperatures using a cluster expansion that includes defect–defect interactions. To overcome the limitations of conventional density functional theory (DFT), we introduce a hybrid approach that incorporates high-accuracy quantum Monte Carlo (QMC) calculations, introducing missing correlations. This strategy allows for accurate estimation of defect energetics and finite-temperature properties. We compute the configurational free energy, defect concentration, and configurational heat capacity, revealing a second-order order–disorder phase transition near the experimental synthesis temperature. Our study provides the first theoretical insight into the thermodynamics of intrinsic defects in MBT. The negative free energy relative to pristine MBT at synthesis temperatures indicates that Mn–Bi antisite formation is thermodynamically spontaneous. We also present a broadly applicable general framework for correcting low-level theoretical theories using highly accurate many-body corrections from QMC.

Ghaffar, Abdul [ORNL] (ORCID:0000000241190168)↗

DOE OSTI2025

Dopant Optimization of Donors in Semiconductor Opening Switches to Eliminate Prepulse

Semiconductor opening switches are solid-state devices capable of delivering nanosecond, hundreds of kilovolts pulses by interrupting kiloamps of current. The interruption of the current occurs in a moderately doped p-region when a high electric field region (HFR) is formed. The HFR occurs because the reverse pumping current cannot be supported by the saturation velocity and majority carrier concentration of the doping level. However, the donor profile also significantly affects the pulse performance. A secondary prepulse occurs if a secondary HFR is formed at the interface of the background n-doping and N+ doping (X n ) . By moving the location of X n deeper into the diode, the effect of the prepulse is reduced. This article investigates the effect of the donor doping profile on the performance metrics of semiconductor opening switches through technology computer-aided design (TCAD) simulations and experimental results. Through a SILVACO TCAD optimization, we designed a P + /p/n - base/n/N + where the intersection of the moderate p-region and intrinsic n-base region (X p ) is at 160 μm and X n is at 220 μm. This profile is fabricated via silicon epitaxy. Experimentally, it is shown that a deep X n (220 μm) compared with a shallow X n (300 μm) reduces the rise time by >5× . In addition, the magnitude of current density during interruption affects the prepulse foot and pulse shape. At lower current densities without the graded donor profile, high peak voltages are not achieved. Comparing the experimental results to the TCAD simulations shows that the model is predictive under high-current densities in the semiconductor opening switch (SOS) regime.

nanosecond pulse power↗

DOE OSTI2020

Charge carrier dynamics in two-dimensional hybrid perovskites: Dion–Jacobson vs. Ruddlesden–Popper phases

Two-dimensional (2D) halide perovskites are promising materials for environmentally stable next-generation optoelectronic device applications. Besides the widely investigated Ruddlesden–Popper (RP) phases, Dion–Jacobson (DJ) phases are attracting considerable attention due to their rapid emergence as efficient solar cell materials. However, there is very little atomistic understanding of the charge carrier dynamics under ambient conditions for these DJ-phases, limiting the possibilities to tune their optoelectronic performances through compositional engineering routes. In this study, by combining nonadiabatic molecular dynamics with time-domain density functional theory methods at room temperature, we compare the dominant non-radiative carrier recombination and dephasing processes in RP and DJ monolayered lead halide perovskites. Our systematic study demonstrates that performance-limiting nonradiative carrier recombination processes greatly depends on the electron–phonon interactions induced by structural fluctuations and instantaneous charge localization in these materials. The stiffer interlayer packing due to the presence of single spacer dications, which separates the lead iodide slabs, reduces the thermal fluctuations in the DJ-phase to a greater extent than that in the RP-phase 2D-perovskites. Specific electronic coupling between the closely spaced lead iodide layers enhances the delocalization of band-edge charge densities in DJ-phase systems. Compared to the RP-phase, reduced inelastic electron–phonon scattering in DJ-phase perovskites significantly limits intrinsic non-radiative recombination processes. The consequent enhancement in the photogenerated charge carrier lifetime makes DJ-phase perovskites potentially suitable for various optoelectronic devices. The computational insights gained from this study allow us to outline a set of robust design principles for DJ-phase perovskites to strategically tune their optoelectronic properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Electronic structure in a rare-earth based nodal-line semimetal candidate PrSbTe

Nodal-line semimetals feature topologically protected band crossings between the bulk valence and conduction bands that extend along a finite dimension in the form of a line or a loop. While ZrSiS and similar materials have attracted extensive research as hosts for the nodal-line semimetallic phase, an alternative avenue has emerged in the form of isostructural rare-earth (RE)-based RESbTe materials. Such systems possess intriguing potentialities for harboring elements of magnetic ordering and electronic correlations owing to the presence of 4f electrons intrinsic to the RE elements. In this study, we have carried out angle-resolved photoemission spectroscopy (ARPES) and thermodynamic measurements in conjunction with first-principles computations on PrSbTe to elucidate its electronic structure and topological characteristics. Magnetic and thermal characterizations indicate the presence of well-localized 4f states with the absence of any discernible phase transition down to 2 K. The ARPES results reveal the presence of gapless Dirac crossings that correspond to a nodal-line along the X-R direction in the three-dimensional Brillouin zone. Furthermore, Dirac crossing that makes up the nodal line, which forms a diamond-shaped nodal plane centered at the center of the Brillouin zone, is also identified within the experimental resolution. In conclusion, this study on the electronic structure of PrSbTe contributes to the understanding of the pivotal role played by spin-orbit coupling in the context of the RESbTe family of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Three Artificial Spintronic Leaky Integrate-and-Fire Neurons

We report that due to their non-volatility and intrinsic current integration capabilities, spintronic devices that rely on domain wall (DW) motion through a free ferromagnetic track have garnered significant interest in the field of neuromorphic computing. Although a number of such devices have already been proposed, they require the use of external circuitry to implement several important neuronal behaviors. As such, they are likely to result in either a decrease in energy efficiency, an increase in fabrication complexity, or even both. To resolve this issue, we have proposed three individual neurons that are capable of performing these functionalities without the use of an external circuitry. To implement leaking, the first neuron uses a dipolar coupling field, the second uses an anisotropy gradient, and the third uses shape variations of the DW track.

neural network crossbar↗

DOE OSTI2021

Early Drug Discovery and Development of Novel Cancer Therapeutics Targeting DNA Polymerase Eta (POLH)

Polymerase eta (or Pol η or POLH) is a specialized DNA polymerase that is able to bypass certain blocking lesions, such as those generated by ultraviolet radiation (UVR) or cisplatin, and is deployed to replication foci for translesion synthesis as part of the DNA damage response (DDR). Inherited defects in the gene encoding POLH (a.k.a., XPV) are associated with the rare, sun-sensitive, cancer-prone disorder, xeroderma pigmentosum, owing to the enzyme’s ability to accurately bypass UVR-induced thymine dimers. In standard-of-care cancer therapies involving platinum-based clinical agents, e.g., cisplatin or oxaliplatin, POLH can bypass platinum-DNA adducts, negating benefits of the treatment and enabling drug resistance. POLH inhibition can sensitize cells to platinum-based chemotherapies, and the polymerase has also been implicated in resistance to nucleoside analogs, such as gemcitabine. POLH overexpression has been linked to the development of chemoresistance in several cancers, including lung, ovarian, and bladder. Co-inhibition of POLH and the ATR serine/threonine kinase, another DDR protein, causes synthetic lethality in a range of cancers, reinforcing that POLH is an emerging target for the development of novel oncology therapeutics. Using a fragment-based drug discovery approach in combination with an optimized crystallization screen, we have solved the first X-ray crystal structures of small novel drug-like compounds, i.e., fragments, bound to POLH, as starting points for the design of POLH inhibitors. The intrinsic molecular resolution afforded by the method can be quickly exploited in fragment growth and elaboration as well as analog scoping and scaffold hopping using medicinal and computational chemistry to advance hits to lead. An initial small round of medicinal chemistry has resulted in inhibitors with a range of functional activity in an in vitro biochemical assay, leading to the rapid identification of an inhibitor to advance to subsequent rounds of chemistry to generate a lead compound. Importantly, our chemical matter is different from the traditional nucleoside analog-based approaches for targeting DNA polymerases.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2024

Computational Modeling of CO2 Capture By Novel PIM-RU in a Fluidized Bed Riser

The interest in carbon capture, storage, and utilization (CCUS) has increased significantly in the past few decades as it can help mitigate the threat of global warming caused by substantial increase in CO2 emissions due to anthropogenic activities. Although several CO2 capture technologies have been developed, porous solid sorbents, which adsorb CO2 by physisorption, are considered promising candidates for post-combustion CO2 capture because of the easier recovery of adsorbed CO2 and high material stability. Some prominent types of porous solid sorbents are metal-organic frameworks (MOFs), Zeolite, mesoporous silica, and polymer-based sorbents (e.g., polymers with intrinsic microporosity or PIM). Oak Ridge National Laboratory (ORNL) recently developed a novel PIM-based sorbent, referred to here as PIM-RU. The main objective of this work is to investigate the CO2 capture performance of this sorbent using computational fluid dynamics (CFD). While process configuration and reactor design are open questions, a fluidized bed riser was selected as the contactor type for this study.

Aziz, Hossain↗

DOE OSTI2020

Subtle metastability of the layered magnetic topological insulator MnBi 2 Te 4 from weak interactions

Abstract Layered quantum materials can host interesting properties, including magnetic and topological, for which enormous computational predictions have been done. Their thermodynamic stability is much less visited computationally, which however determines the existence of materials and can be used to guide experimental synthesis. MnBi 2 Te 4 is one of such layered quantum materials that was predicted to be an intrinsic antiferromagnetic topological insulator, and later experimentally realized but in a thermodynamically metastable state. Here, using a combined first-principles-based approach that considers lattice, charge, and spin degrees of freedom, we investigate the metastability of MnBi 2 Te 4 by calculating the Helmholtz free energy for the reaction Bi 2 Te 3 + MnTe → MnBi 2 Te 4 . We identify a temperature range (~500–873 K) in which the compound is stable with respect to the competing binary phases, consistent with experimental observation. We validate the predictions by comparing the calculated specific heats contributed from different degrees of freedom with experimental results. Our findings indicate that the degrees of freedom responsible for the van der Waals interaction, lattice vibration, magnetic coupling, and nontrivial band topology in MnBi 2 Te 4 not only enable emergent phenomena but also play a crucial role in determining its thermodynamic stability. This conclusion lays the foundation for the future computational material synthesis of novel layered systems.

36 MATERIALS SCIENCE↗

DOE OSTI2021

GeantV: Results from the Prototype of Concurrent Vector Particle Transport Simulation in HEP

Full detector simulation was among the largest CPU consumers in all CERN experiment software stacks for the first two runs of the Large Hadron Collider. In the early 2010s, it was projected that simulation demands would scale linearly with increasing luminosity, with only partial compensation from increasing computing resources. The extension of fast simulation approaches to cover more use cases that represent a larger fraction of the simulation budget is only part of the solution, because of intrinsic precision limitations. The remainder corresponds to speeding up the simulation software by several factors, which is not achievable by just applying simple optimizations to the current code base. In this context, the GeantV R&D project was launched, aiming to redesign the legacy particle transport code in order to benefit from features of fine-grained parallelism, including vectorization and increased locality of both instruction and data. This paper provides an extensive presentation of the results and achievements of this R&D project, as well as the conclusions and lessons learned from the beta version prototype.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Mapping 3D grain and precipitate structure during in situ mechanical testing of open-cell metal foam using micro-computed tomography and high-energy X-ray diffraction microscopy

Open-cell metal foams are ultra-low-density cellular metals with complex hierarchical structures that span bulk, cell, ligament, and sub-ligament scales and give rise to desirable properties such as high strength-to-weight ratio and excellent energy absorption. Although literature suggests that intrinsic material structures at sub-ligament length scales (e.g., grains and precipitates) play an important role in mechanical behavior of open-cell metal foams, there are very few experimental measurements of such structures in three dimensions and for meaningful volumes of foam. This study seeks to map and track the three-dimensional (3D) grain and precipitate structures of an intact volume of open-cell aluminum foam by advancing microstructural characterization techniques that leverage X-ray micro-computed tomography (μCT) and far-field high-energy X-ray diffraction microscopy (FF-HEDM). A 6%-relative-density aluminum foam sample was mechanically tested in compression while μCT and FF-HEDM measurements were collected at interrupted loading states at beamline 1-ID of the Advanced Photon Source. Further, a new scanning strategy and reconstruction algorithm were established to enable characterization of a foam volume with diameter approximately four times wider than the nominal width of the X-ray beam. The result is a set of maps that detail both the 3D grain and precipitate structures throughout the foam volume at successive strain steps. A novel grain tracking procedure was developed to track individual grains within the foam volume by accounting for the large rigid-body motions that individual ligaments can undergo during mechanical loading. The ability to track grains and precipitate structures in three dimensions throughout large bulk deformation of ultra-low-density polycrystalline materials enables new possibilities for validating numerical models and investigating local failure mechanisms. Furthermore, the methods and procedures developed in this study could be applied to other ultra-low-density structures, such as additively manufactured lattices.

36 MATERIALS SCIENCE↗

DOE OSTI2024

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗