Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational structural dynamics”

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 1,333 records · Page 74

Direct Ab Initio Simulation of the Synthesis of BaZrO 3 and the Microstructure Impacts on Proton Transport

Controlling and predicting the processing-structure-performance relationship in functional materials is a grand challenge in materials science, with important implications for a wide range of emerging applications; a high fidelity understanding of the performance impact of microstructures formed under synthesis conditions is required to develop advanced materials, such as solid-state fuel cells and electrolyzers. Using the ceramic BaZrO 3 as a case study, we directly simulate the synthesis and investigate how proton transport is dictated by microstructures. We develop a framework that couples density functional theory (DFT), machine-learning interatomic potential (MLIP) driven molecular dynamics, and grand canonical Monte Carlo to perform large-scale, microstructure-resolved, atomistic simulations of proton transport in experimentally representative polycrystalline structures. Our fully ab initio approach, using a MLIP as a proxy for DFT, allows us to quantify the competition between two distinct diffusion mechanisms: one associated with grain-boundary regions and another within grains. When the impacts of grain boundaries are taken into account, proton transport exhibits substantial deviation from the bulk oxide limit. This addresses long-standing discrepancies between theory and experiments. Our integrated approach provides atomistic insight into microstructure-dependent proton pathways in BaZrO 3 and establishes a general protocol for predicting processing-structure-performance relationships.

organic↗

Stochastic finite volume method for uncertainty quantification of transient flow in gas pipeline networks

We develop a weakly intrusive framework to simulate the propagation of uncertainty in solutions of generic hyperbolic partial differential equation systems on graph-connected domains with nodal coupling and boundary conditions. The method is based on the Stochastic Finite Volume (SFV) approach and can be applied for uncertainty quantification (UQ) of the dynamical state of fluid flow over actuated transport networks. The numerical scheme has specific advantages for modeling intertemporal uncertainty in time-varying boundary parameters, which cannot be characterized by strict upper and lower (interval) bounds. We describe the scheme for a single pipe, and then formulate the controlled junction Riemann problem (JRP) that enables the extension to general network structures. In conclusion, we demonstrate the method's capabilities and performance characteristics using a standard benchmark test network.

97 MATHEMATICS AND COMPUTING↗

Inverse Methods - Users Manual 5.8

The inverse methods team provides a set of tools for solving inverse problems in structural dynamics and thermal physics, and also sensor placement optimization via Optimal Experimental Design (OED). These methods are used for designing experiments, model calibration, and verification/validation analysis of weapons systems. This document provides a user’s guide to the input for the three apps that are supported for these methods. Details of input specifications, output options, and optimization parameters are included.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sierra/SD – User's Manual – 5.18

Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high-fidelity, validated models used in modal, vibration, static and shock analysis of weapons systems. This document provides a user’s guide to the input for Sierra/SD. Details of input specifications for the different solution types, output options, element types and parameters are included. The appendices contain detailed examples, and instructions for running the software on parallel platforms.

97 MATHEMATICS AND COMPUTING↗

Surface structure of linear nanopores in amorphous silica: Comparison of properties for different pore generation algorithms

In this work, we compare the surface structure of linear nanopores in amorphous silica (a-SiO 2 ) for different versions of “pore drilling” algorithms (where the pores are generated by the removal of atoms from the preformed bulk a-SiO 2 ) and for “cylindrical resist” algorithms (where a-SiO 2 is formed around a cylindrical exclusion region). After adding H to non-bridging O, the former often results in a moderate to high density of surface silanol groups, whereas the latter produces a low density. The silanol surface density for pore drilling can be lowered by a final dehydroxylation step, and that for the cylindrical resist approach can be increased by a final hydroxylation step. In this respect, the two classes of algorithms are complementary. We focus on the characterization of the chemical structure of the pore surface, decomposing the total silanol density into components corresponding to isolated and vicinal mono silanols and geminal silanols. The final dehyroxylation and hydroxylation steps can also be tuned to better align some of these populations with the target experimental values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Machine learning assisted search for Fe–Co–C ternary compounds with high magnetic anisotropy

We employ a machine learning (ML)-guided framework to explore rare earth free magnetic materials, specifically focusing on Fe–Co–C ternary compounds for potential use in permanent magnets. Utilizing a specifically trained crystal graph convolutional neural network model, we efficiently screen a vast space of nearly a million substitutional structures to select 620 promising structures for further investigation by first-principles calculation. We predict five low-energy metastable Fe–Co–C compounds with formation energy less than 150 meV/atom above the convex hull. These compounds exhibit high magnetization (Js > 1.0 T) and significant magnetic anisotropy (K1 > 1.0 MJ/m3), making them promising candidates for permanent magnet applications. The phonon calculations indicate these compounds are dynamically stable. Our ML-guided framework demonstrates the utility of rapidly identifying novel materials with tailored magnetic properties.

36 MATERIALS SCIENCE↗

Investigation of Farm-to-Farm Interactions and Blockage Effects from AWAKEN using Large-Scale Numerical Simulations

A large-scale numerical computation of five wind farms was performed as a part of the American Wake Experiment (AWAKEN). This high-fidelity computation used the ExaWind/AMR-Wind LES solver to simulate a 100km × 100km domain containing 541 turbines under unstable atmospheric conditions matching previous measurements. The turbines were represented by Joukowski and OpenFAST coupled actuator disk models. Results of this qualitative comparison illustrate the interactions with wind farms with large scale ABL structures in the flow, as well as the extent of downstream wake penetration in the flow and blockage effects around wind farms.

Cheung, Lawrence C.↗

On the Stark Effect of the O I 777-nm Triplet in Plasma and Laser Fields

The O I 777-nm triplet transition is often used for plasma density diagnostics. It is also employed in nonlinear optics setups for producing quasi-comb structures when pumped by a near-resonant laser field. Here, we apply computer simulations to situations of the radiating atom subjected to the plasma microfields, laser fields, and both perturbations together. Our results, in particular, resolve a controversy related to the spectral line anomalously broadened in some laser-produced plasmas. The importance of using time-dependent density matrix is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Control of Permanent Porosity in Type 3 Porous Liquids via Solvent Clustering

Porous liquids (PLs) are an exciting new class of materials for carbon capture due to their high gas adsorption capacity and ease of industrial implementation. They are composed of sorbent particles suspended in a nonadsorbed solvent, forming a liquid with permanent porosity. While PLs have a vast number of potential compositions based on the number of solvents and sorbent materials available, most of the research has been focused on the selection of the sorbent rather than the solvent. Therefore, PL design criteria on the supramolecular structures of the solvent are explored to create a fundamental understanding of how the solvent enables PL formation for rapid discovery of new PL compositions. Atomistic molecular dynamics simulation of eight solvents with a range of molecular sizes, shapes, and intramolecular bonding was performed, identifying that the shape and size of molecular clusters formed in the solvent are the driving predictor of PL formation rather than the size of the individual solvent molecule. The results demonstrate a significant departure from common approaches to PL formation based on the steric exclusion of solvent molecules from the sorbent via the size of the pore aperture. A modeling and experimental validation study further supports these findings. In conclusion, through this computational material design study, a previously unexplored mechanism in PL formation, solvent–solvent clustering, is identified as a critical factor for the accelerated discovery of liquid phase carbon capture materials.

Carbon capture↗

Emergent disorder and mechanical memory in periodic metamaterials

Ordered mechanical systems typically have one or only a few stable rest configurations, and hence are not considered useful for encoding memory. Multistable and history-dependent responses usually emerge from quenched disorder, for example in amorphous solids or crumpled sheets. In contrast, due to geometric frustration, periodic magnetic systems can create their own disorder and espouse an extensive manifold of quasi-degenerate configurations. Inspired by the topological structure of frustrated artificial spin ices, we introduce an approach to design ordered, periodic mechanical metamaterials that exhibit an extensive set of spatially disordered states. While our design exploits the correspondence between frustration in magnetism and incompatibility in meta-mechanics, our mechanical systems encompass continuous degrees of freedom, and thus generalize their magnetic counterparts. We show how such systems exhibit non-Abelian and history-dependent responses, as their state can depend on the order in which external manipulations were applied. We demonstrate how this richness of the dynamics enables to recognize, from a static measurement of the final state, the sequence of operations that an extended system underwent. Thus, multistability and potential to perform computation emerge from geometric frustration in ordered mechanical lattices that create their own disorder.

36 MATERIALS SCIENCE↗

Dynamic and reversible transformations of subnanometre-sized palladium on ceria for efficient methane removal

Reversibly adjusting the active structures of supported metal catalysts in response to dynamic working conditions has long been pursued. Here we report the reaction-environment-modulated transformations of subnanometre-sized Pd on CeO 2 for efficient methane removal, leveraging the reaction environments at different stages of automotive exhaust aftertreatment. During the cold start of vehicles, inactive Pd 1 single atoms are readily transformed into PdO $x$ subnanometre clusters by CO even at room temperature with excess O 2 , resulting in boosted low-temperature CH 4 oxidation. At elevated temperatures, dispersion of PdO $x$ cluster into Pd 1 against metal sintering renders outstanding hydrothermal stability to the catalyst, to be activated during the next vehicle start. Combined experimental and computational studies elucidate the dynamically evolved Pd speciation on CeO 2 at an atomic level. In conclusion, modulating the reversible nature of supported metals helps overcome the long-existing trade-off between low-temperature activity and high-temperature stability, also providing a new paradigm for designing intelligent catalysts that brings single-atom/cluster catalysts closer to real applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-resolution partial differential equations preserved learning framework for spatiotemporal dynamics

Traditional data-driven deep learning models often struggle with high training costs, error accumulation, and poor generalizability in complex physical processes. Physics-informed deep learning (PiDL) addresses these challenges by incorporating physical principles into the model. Most PiDL approaches regularize training by embedding governing equations into the loss function, yet this depends heavily on extensive hyperparameter tuning to weigh each loss term. To this end, we propose to leverage physics prior knowledge by “baking” the discretized governing equations into the neural network architecture via the connection between the partial differential equations (PDE) operators and network structures, resulting in a PDE-preserved neural network (PPNN). This method, embedding discretized PDEs through convolutional residual networks in a multi-resolution setting, largely improves the generalizability and long-term prediction accuracy, outperforming conventional black-box models. The effectiveness and merit of the proposed methods have been demonstrated across various spatiotemporal dynamical systems governed by spatiotemporal PDEs, including reaction-diffusion, Burgers’, and Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗

SAXS-guided unbiased coarse-grained Monte Carlo simulation for identification of self-assembly nanostructures and dimensions

Recent studies have shown that solvated amphiphiles can form nanostructured self-assemblies called dynamic binary complexes (DBCs) in the presence of ions. Since the nanostructures of DBCs are directly related to their viscoelastic properties, it is important to understand how the nanostructures change under different solution conditions. However, it is challenging to obtain a three-dimensional molecular description of these nanostructures by utilizing conventional experimental characterization techniques or thermodynamic models. To this end, we combined the structural data from small angle X-ray scattering (SAXS) experiments and thermodynamic knowledge from coarse-grained Monte Carlo (CGMC) simulations to identify the detailed three-dimensional nanostructure of DBCs. Specifically, unbiased CGMC simulations are performed with SAXS-guided initial conditions, which aids us to sample accurate nanostructures in a computationally efficient fashion. As a result, an elliptical bilayer nanostructure is obtained as the most probable nanostructure of DBCs whose dimensions are validated by scanning electron microscope (SEM) images. Then, utilizing the obtained molecular model of DBCs, here we could also explain the pH tunability of the system. Overall, our results from SAXS-guided unbiased CGMC simulations highlight that using potential energy combined with SAXS data, we can distinguish otherwise degenerate nanostructures resulting from the inherent ambiguity of SAXS patterns.

36 MATERIALS SCIENCE↗

Regulating loops in de Sitter spacetime

Perturbative quantum field theory (QFT) calculations in de Sitter space are riddled with contributions that diverge over time. These contributions often arise from loop integrals, which are notoriously hard to compute in de Sitter. We discuss an approach to evaluate loop integrals that contribute to equal-time correlators of a scalar field theory in a fixed de Sitter background. Our method is based on the Mellin-Barnes representation of correlation functions, which allows us to regulate loop divergences by adjusting the masses of the fields, or by gently deforming the underlying de Sitter spacetime. The resulting expressions have a similar structure as a standard answer from dimensional regularization in flat space QFT. These features of the regulator are illustrated with two examples, worked out in detail. Along the way, we illuminate the physical origin of these divergences and their interpretation with the machinery of the dynamical renormalization group. Our approach regulates the IR divergences of massless and massive particles in the same way. For massless scalars, the loop corrections can be incorporated as systematic improvements to the stochastic inflation framework, allowing for a more precise description of the IR dynamics of such fields in de Sitter. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Sierra/SolidMechanics 4.56 Goodyear User's Guide

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional finite element analysis code for solids and structures subjected to extensive contact and large deformations, encompassing explicit and implicit dynamic as well as quasistatic loading regimes. This document supplements the primary Sierra/SM 4.56 User’s Guide, describing capabilities specific to Goodyear analysis use cases, including additional implicit solver options, material models, finite element formulations, and contact settings.

42 ENGINEERING↗

Developing a complete AI-accelerated workflow for superconductor discovery

The quest to identify new superconducting materials with enhanced properties is hindered by the prohibitive cost of computing electron-phonon spectral functions, severely limiting the materials space that can be explored. Here, we introduce a Bootstrapped Ensemble of Equivariant Graph Neural Networks (BEE-NET), a machine-learning model trained to predict the Eliashberg spectral function and superconducting critical temperature with a mean-absolute-error of 0.87 K relative to DFT-based Allen-Dynes calculations. Intriguingly, BEE-NET achieves a true-negative-rate of 99.4%, enabling highly efficient screening for the rare property of superconductivity. Integrated into a multi-stage, AI-accelerated discovery pipeline that incorporates elemental-substitution strategies and machine-learned interatomic potentials, our workflow reduced over 1.3 million candidate structures to 741 dynamically and thermodynamically stable compounds with DFT-confirmed T c > 5 K. We report the successful synthesis and experimental confirmation of superconductivity in two of these previously unreported compounds. This study establishes a data-driven framework that integrates machine learning, quantum calculations, and experiments to systematically accelerate superconductor discovery.

Gibson, Jason B. [Quantum Formatics, Cambridge, MA↗

Surface reconstructions and electronic structure of metallic delafossite thin films

The growing interest in the growth and study of thin films of low-dimensional metallic delafossites, with the general formula ABO2, is driven by their potential to exhibit electronic and magnetic characteristics that are not accessible in bulk systems. The layered structure of these compounds introduces unique surface states as well as electronic and structural reconstructions, making the investigation of their surface behavior pivotal to understanding their intrinsic electronic structure. In this work, we study the surface phenomena of epitaxially grown PtCoO2, PdCoO2, and PdCrO2 films, utilizing a combination of molecular-beam epitaxy and angle-resolved photoemission spectroscopy. Through precise control of surface termination and treatment, we discover a pronounced 3×3 surface reconstruction in PtCoO2 films and PdCoO2 films, alongside a 2 × 2 surface reconstruction observed in PdCrO2 films. These reconstructions have not been reported in prior studies of delafossites. Furthermore, our computational investigations demonstrate the BO2 surface’s relative stability compared to the A-terminated surface and the significant reduction in surface energy facilitated by the reconstruction of the A-terminated surface. These experimental and theoretical insights illuminate the complex surface dynamics in metallic delafossites, paving the way for future explorations of their distinctive properties in low-dimensional studies.

Materials Science↗

Ensemble Simulation Techniques and Fast Randomized Algorithms

The major goals of the project were to develop and analyze new ensemble simulation techniques, including trajectory stratification and preconditioned MCMC techniques, as well as develop fast numerical linear algebra techniques closely related to ensemble simulation ideas. The trajectory stratification techniques involve simulating in parallel short trajectory fragments of a Markov process confined to a specific region of space‐time and then patching together the statistics gathered to assemble estimates of very general dynamical properties. We have also developed this approach for rare event simulation and extended the techniques to applications requiring a more general framework (such as electronic structure calculations). The preconditioned MCMC techniques involve simulating multiple Markov chains in parallel and then using information from the ensemble to speed the mixing of each individual chain. The fast randomized linear algebra methods are motivated by the diffusion Monte Carlo technique, but are applicable to finding the dominant eigenvalue of (almost) general matrices. For most non‐negative matrices, the schemes result in an error (compared to the power method) that is constant in the dimension of the problem. For more general matrices, we see a very clear sublinear cost trend in computational tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗