Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic simulation”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Machine-learning force-field models for dynamical simulations of metallic magnets

We review recent advances in machine-learning (ML) force-field methods for Landau–Lifshitz–Gilbert simulations of itinerant electron magnets, focusing on their scalability and transferability. Built on the principle of locality, a deep neural-network model is developed to efficiently and accurately predict electron-mediated forces governing spin dynamics. Symmetry-aware descriptors constructed through a group-theoretical approach ensure rigorous incorporation of both lattice and spin-rotation symmetries. The framework is demonstrated using the prototypical s-d exchange model widely employed in spintronics. ML-enabled large-scale simulations reveal novel nonequilibrium phenomena, including anomalous coarsening of tetrahedral spin order on the triangular lattice and the freezing of phase-separation dynamics in lightly hole-doped, strong-coupling square-lattice systems. These results establish ML force-field frameworks as scalable, accurate, and versatile tools for modeling nonequilibrium spin dynamics in itinerant magnets.

Artificial neural networks↗
NASA NTRS1995

Simulated dynamic response of a multi-stage compressor with variable molecular weight flow medium

A mathematical model of a multi-stage compressor with variable molecular weight flow medium is derived. The modeled system consists of a five stage, six cylinder, double acting, piston type compressor. Each stage is followed by a water cooled heat exchanger which serves to transfer the heat of compression from the gas. A high molecular weight gas (CFC-12) mixed with air in varying proportions is introduced to the suction of the compressor. Condensation of the heavy gas may occur in the upper stage heat exchangers. The state equations for the system are integrated using the Advanced Continuous Simulation Language (ACSL) for determining the system's dynamic and steady state characteristics under varying operating conditions.

Babcock, Dale A.↗
DOE OSTI2022

Atomic-scale mechanism of carbon nucleation from a deep crustal fluid by replica exchange reactive molecular dynamics simulation

Here we present a mechanistic model of carbon nucleation and growth from a fluid at elevated temperature (T) and pressure conditions, typical of those found in the shallow Earth’s lithosphere. Our model uses a replica exchange reactive molecular dynamics framework in which molecular configurations are swapped between adjacent T replica at regular intervals according to underlying statistical mechanics. This framework allows predicting complex molecular structures and thermodynamics while remaining computationally efficient. Here we simulate the reactivity of an unstable mixture of CO 2 and CH 4 at 1000 K and 1 GPa. We find that the path to thermodynamic equilibrium is initially entropy-driven, producing a diversity of short-lived species, including various alcohols with intermediate carbon oxidation states. Cyclic and polycyclic radicals that are sometimes resonance-stabilized form next and set the stage for carbon nucleation. The carbon exsolution process releases abundant water, is exothermic and starts with the nucleation of a large aggregate of hydrogenated graphene flakes from covalently bonded polycyclic units. The carbon backbone of this nucleus subsequently grows into a hydrogen-depleted fullerene-like structure, before evolving toward a partially bilayered graphene layer. Overall, our results show that the mechanism of graphitic C formation is certainly not bimolecular, and that it may involve a combination of key condensation and radical chain reactions. This will help understand the isotopic, and reactive characteristics of carbon-bearing fluids during their upward transit through the Earth’s mantle and crust. Moreover, the mechanistic insights outlined here present intriguing similarities with the process of soot and interstellar dust formation, which suggests that the widespread distribution of abiotic polyaromatic and graphitic material on Earth and beyond may reflect the prevalence of a fundamental chemical pathway.

58 GEOSCIENCES↗
DOE OSTI2023

Constructing Neural Network Based Models for Simulating Dynamical Systems

Dynamical systems see widespread use in natural sciences like physics, biology, and chemistry, as well as engineering disciplines such as circuit analysis, computational fluid dynamics, and control. For simple systems, the differential equations governing the dynamics can be derived by applying fundamental physical laws. However, for more complex systems, this approach becomes exceedingly difficult. Data-driven modeling is an alternative paradigm that seeks to learn an approximation of the dynamics of a system using observations of the true system. In recent years, there has been an increased interest in applying data-driven modeling techniques to solve a wide range of problems in physics and engineering. Here this article provides a survey of the different ways to construct models of dynamical systems using neural networks. In addition to the basic overview, we review the related literature and outline the most significant challenges from numerical simulations that this modeling paradigm must overcome. Based on the reviewed literature and identified challenges, we provide a discussion on promising research areas.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2024

Experimental Validation of Charged Lunar Dust Dynamics Simulants

The goal of this work is to experimentally verify the electrostatic interaction physics that have been incorporated into the discrete element method (DEM) modeling software package during the FY20 CIF project titled “Charged Particle Dynamics in the Lunar Environment.” Current state-of-the-art (SOA) granular gas dynamics models used to explain rocket plume impingement physics have not taken into account the natural or induced electrostatic environment of the lunar surface, nor the effect of charged regolith grains being present in the plume plasma. This work improved upon the SOA by adding long-range and contact inter-particle electrostatic interactions to a granular mechanics DEM modeling software package. Tribocharging mechanics between spherical particles has been experimentally examined along with the natural electrostatic plasma environments produced via emission from an electron gun and ultraviolet lamp.

Jay Phillips↗
DOE OSTI2021

Validation of Computational Fluid Dynamics Simulations for Biological Performance Assessment in Hydropower units (Final Report)

The biological performance assessment (BioPA) toolset developed by Pacific Northwest National Laboratory (PNNL) estimates the relative biological performance of fish passage at a hydroelectric power turbine unit. The tool is based on the use of computational fluid dynamics (CFD) and fish biological response relationships. The recent release, BioPA-v3, is based on directly computed trajectory and collision of material Lagrangian particles using CFD simulation codes rather than the prior version that relies on Tecplot to compute streamtrace trajectories. Before modifying the toolset, a series of validation tests were performed at the various steps of modification in the toolset. Validation is a critical step of any numerical investigation that reflects the accuracy and reliability of the predicted results. It raises the confidence level of the user to use the modified version of the BioPA toolset. Several test cases were simulated and compared, where available, to observed data. The trajectory and collision of the small spherical and cylindrical particles in a water flume were compared to in-house experiments. The CFD predicted collision rate and flow field compared well with experimental observation for vane array and large cylinder as target bodies. Next, the CFD-predicted flow field and hydraulic performance of a laboratory-scale model of a Francis turbine was also successfully validated. Note that the trajectory of the particles is significantly affected by the flow field in such extreme conditions. In addition to the particle trajectories and flow field, the collision detection method employed in the CFD simulations was also successfully validated. The CFD predicted impact velocity, collision time, velocity, and trajectory of a sphere excellently matched with analytical value for a bouncing ball in the elastic collision. A similar approach was also tested and successfully validated for a collision of sphere with a 45° inclined plane. After successfully validating different cases, the BioPA toolset was modified to use direct output of the CFD prediction and the new version can be used in evaluating biological performance at hydroelectric turbines.

13 HYDRO ENERGY↗
NASA NTRS2019

AladynPi – Adaptive Neural Network Molecular Dynamics Simulation Code with Physically Informed Potential: Computational Materials Mini-Application

This report provides an overview and description of commands used in the Computational Materials mini-application, AladynPi. AladynPi is an extension of a previously released mini-application, Aladyn (https://github.com/nasa/aladyn; Yamakov, V.I., and Glaessgen, E.H., NASA/TM-2018-220104). Aladyn and AladynPi are basic molecular dynamics codes written in FORTRAN 2003, which are designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory method. An input for the ANN is a set of structure coefficients, characterizing the local atomic environment of each atom, for which the atomic energy is obtained in the ANN inference process. In Aladyn, the ANN gives directly the energy of interatomic interactions. In AladynPi, the ANN gives optimized parameters for a predefined empirical function, known as bond-order-potential (BOP). The parameterized BOP function is then used to calculate the energy. AladynPi code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing manycore architectures based on graphic processing units (GPUs). The effort is supported by the High Performance Computing incubator (HPCi) project at NASA Langley Research Center.

Yamakov, Vesselin I.↗
NASA NTRS1999

MRI Measurements and Granular Dynamics Simulation of Segregation of Granular Mixture

A counter intuitive axial segregation phenomenon in a rotating horizontal cylinder has recently captured attention of many researchers in different disciplines. There is a growing consensus that the interplay between the particle dynamics and the evolution of the internal structure during the segregation process must be carefully investigated. Magnetic resonance imaging (MRI) has been used to non-invasively obtain much needed dynamic/static information such as velocity and concentration profiles, and it has proven to be capable of depicting the evolution of segregation processes. Segregation in a rotating cylinder involves two processes: the first is to transport small particles in the radial direction to form a radial core, and the second is to transform the radial core into axially segregated bands. Percolation and/or "stopping" have been proposed as mechanisms for the radial segregation. As to mechanisms for axial band formation, much less is known. The difference in the dynamic angle of repose has been proposed to segregate different components in the axial direction. Recently, Hill and Kakalios have reported that particles mix or demix depending upon the competition between diffusion and preferential drift whose order can be determined by the dynamic angle of repose through the adjustment of the rotation rate. We claim that the dynamic angle of repose could be one of the causes, however, it fails to offer reasonable explanations for certain aspects of the axial migration. For example, we always observe that the radial segregation precedes the axial segregation and small particles migrate in the radial direction to form an axially extended radial core. It then transforms into axially segregated bands. By definition, the effects of the dynamic angle of repose are restricted near the free surface where the flowing layer is present. However, during the process of transforming from the radially segregated core to axially segregated bands, small particles located in the deep core region, which is untouched by the flowing layer, also completely disappear. Usually, the dynamics angle of repose are uniquely defined for individual species to characterize particle properties, and the dynamic angle of repose thus defined provides little information for the dynamic angle of repose of the mixture since the concentration ratio and the internal packing structure do not remain the same during the segregation processes. Under microgravity environment, the dynamics angle of repose argument does not hold since there is simply no flowing layer to influence/determine the preferred directions of segregation. We have thus designed an experiment so that the effects of the dynamic angle of repose can be minimized by filling the cylinder almost completely full. Small particles still formed a radial core and also migrated to form axial bands. As ground based experiments we have designed and conducted both 2D and 3D segregation experiments. The 2D experiments are performed using a thin cylinder (the gap between two end caps is about 5 mm) filled with different combinations of particles. The 3D experiments are conducted with a long cylinder of its length and diameter of 27cm and 7cm, respectively. Results of 2D experiments indicate that different mechanisms govern particle motion in regions near and far from the axis of rotation. Results of 3D experiments indicate that a series of collapses of microstructures of particle packing (micro-collapses) may be responsible for the creation of voids for small particles to migrate through in the axial direction. We have successfully eliminated the dynamic angle of repose as a cause for segregation, however, by almost completely filling the cylinder with the particles, we have lost an opportunity to investigate a possibility of particle "mobility" being a cause for segregation which requires a flowing surface but not the difference in the angle of repose. This is currently being investigated.

Nakagawa, M.↗
DOE OSTI2024

Role of interfaces on the mechanical response of accumulative roll bonded nanometallic laminates investigated via dislocation dynamics simulations

Unraveling the effects of continuous dislocation interactions with interfaces, particularly at the nanometer length scales, is key to a broader understanding of plasticity, to material design and to material certification. To this end, this work proposes a novel discrete dislocation dynamics-based model for dislocation interface interactions tracking the fate of residual dislocation on interfaces. This new approach is used to predict the impact of dislocation/interface reactions on the overall mechanical behavior of accumulative roll bonded nanometallic laminates. The framework considers the dynamic evolution of the interface concurrent with a large network of dislocations, thus, accounting for the local short and long range effects of the dislocations under the external boundary conditions. Specifically, this study focuses on two-phase Fe/Cu nanometallic laminates, and investigates the role of the underlying elastic and plastic contrast of the Fe and the Cu layers on the composite response of the material. Moreover, the role of initial microstructures, resulting from processing is also investigated. Subsequently, the model is used to examine the effect of layer thickness and interface orientation relationship on the residual stresses of the relaxed microstructure. The associated mechanical response of these laminates are compared when loaded under normal direction compression, as well as shear compression. Finally, this work predicts a dominant effect of the layer thickness, as compared to the interface orientation relationship, on the macroscopic response and on the residual stresses of these nanolaminates, while the local dislocation transmission propensity through the interface is significantly influenced by the corresponding orientation relationship.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantum-Based Molecular Dynamics Simulations Using Tensor Cores

Tensor cores, along with tensor processing units, represent a new form of hardware acceleration specifically designed for deep neural network calculations in artificial intelligence applications. Tensor cores provide extraordinary computational speed and energy efficiency but with the caveat that they were designed for tensor contractions (matrix–matrix multiplications) using only low-precision floating-point operations. Despite this perceived limitation, we demonstrate how tensor cores can be applied with high efficiency to the challenging and numerically sensitive problem of quantum-based Born–Oppenheimer molecular dynamics, which requires highly accurate electronic structure optimizations and conservative force evaluations. The interatomic forces are calculated on-the-fly from an electronic structure that is obtained from a generalized deep neural network, where the computational structure naturally takes advantage of the exceptional processing power of the tensor cores and allows for high performance in excess of 100 Tflops on a single Nvidia A100 GPU. Stable molecular dynamics trajectories are generated using the framework of extended Lagrangian Born–Oppenheimer molecular dynamics, which combines computational efficiency with long-term stability, even when using approximate charge relaxations and force evaluations that are limited in accuracy by the numerically noisy conditions caused by the low-precision tensor core floating-point operations. A canonical ensemble simulation scheme is also presented, where the additional numerical noise in the calculated forces is absorbed into a Langevin-like dynamics.

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

Accelerating computational fluid dynamics simulation of post-combustion carbon capture modeling with MeshGraphNets

Packed columns are commonly used in post-combustion processes to capture CO 2 emissions by providing enhanced contact area between a CO 2 -laden gas and CO 2 -absorbing solvent. To study and optimize solvent-based post-combustion carbon capture systems (CCSs), computational fluid dynamics (CFD) can be used to model the liquid–gas countercurrent flow hydrodynamics in these columns and derive key determinants of CO 2 -capture efficiency. However, the large design space of these systems hinders the application of CFD for design optimization due to its high computational cost. In contrast, data-driven modeling approaches can produce fast surrogates to study large-scale physics problems. We build our surrogates using MeshGraphNets (MGN), a graph neural network framework that efficiently learns and produces mesh-based simulations. We apply MGN to a random packed column modeled with over 160K graph nodes and a design space consisting of three key input parameters: solvent surface tension, inlet velocity, and contact angle. Our models can adapt to a wide range of these parameters and accurately predict the complex interactions within the system at rates over 1700 times faster than CFD, affirming its practicality in downstream design optimization tasks. This underscores the robustness and versatility of MGN in modeling complex fluid dynamics for large-scale CCS analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of TATB single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature–pressure–volume equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the two existing nonreactive FFs are more accurate in describing TATB’s physical properties compared to the two variants of the ReaxFF reactive FF considered. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. So this work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗