Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamics 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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Ab initio molecular dynamics (AIMD) simulations of NaCl, UCl 3 molten salts

Ab initio molecular dynamics (AIMD) simulations are used to calculate select thermophysical and thermodynamic properties of NaCl, UCl and NaCl-UCl 3 molten salts. Following established approaches, the AIMD simulations include a model for Van der Waals interactions. The Langreth & Lundqvist (vDW-DF2), DFT-D3, and density-dependent energy correction (DFT-dDsC) dispersion models are first tested for molten NaCl in order to assess their accuracy for density and heat capacity predictions across a range of temperatures. Based on the NaCl results, the vdW-DF2 and DFT-dDsC methods are extended to the UCl 3 system and compared to available experimental data. Next, mixtures of NaCl-UCl 3 are investigated and analyzed with respect to the deviation from ideal solution behavior. For the DFT-dDsC simulations, density deviates by up to 2% from an ideal mixture, with the maximum occurring close to the eutectic composition. The mixing energy also deviates from an ideal solution and exhibits a minimum of about -0.074 to eV per formula unit, again close to the eutectic composition. Finally, the compressibility and species diffusivity of the pure and mixed salt systems are calculated. The diffusivities are slightly reduced in the mixed compared to the pure systems and the compressibilities loosely follow a linear correlation as a function of the UCl 3 composition. The trends observed for mixing properties are rationalized by correlating them to the coordination chemistry, which emphasizes the importance of maintaining the extended network formed by U and Cl ions as the NaCl concentration increases. The concentration at which breakdown of the extended network occurs, roughly coincides with the minimum of both the mixing energy and the deviation from ideal solution behavior for density.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Machine Learning-Based Upscaling of Finite-Size Molecular Dynamics Diffusion Simulations for Binary Fluids

Molecular diffusion coefficients calculated using molecular dynamics (MD) simulations suffer from finite-size (i.e., finite box size and finite particle number) effects. Results from finite-sized MD simulations can be upscaled to infinite simulation size by applying a correction factor. For self-diffusion of single-component fluids, this correction has been well-studied by many researchers including Yeh and Hummer (YH); for binary fluid mixtures, a modified YH correction was recently proposed for correcting MD-predicted Maxwell–Stephan (MS) diffusion rates. In this study we use both empirical and machine learning methods to identify improvements to the finite-size correction factors for both self-diffusion and MS diffusion of binary Lennard-Jones (LJ) fluid mixtures. Using artificial neural networks (ANNs), the error in the corrected LJ fluid diffusion is reduced by an order of magnitude versus existing YH corrections, and the ANN models perform well for mixtures with large dissimilarities in size and interaction energies where the YH correction proves insufficient.

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

Atomistic Simulation of Interfaces in Materials of Solid State Ionics

The possibilities of describing correctly interfaces of different types in solids within a computer experiment using molecular statics simulation, molecular dynamics simulation, and quantum chemical calculations are discussed. Heterophase boundaries of various types, including grain boundaries and solid electrolyte‒solid electrolyte and ionic conductor‒electrode material interfaces, are considered. Specific microstructural features and mechanisms of the ion transport in real heterophase structures (cationic conductor‒metal anode and anionic conductor‒cathode) existing in solid state ionics devices (such as solid-state batteries and fuel cells) are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

On fast simulation of dynamical system with neural vector enhanced numerical solver

The large-scale simulation of dynamical systems is critical in numerous scientific and engineering disciplines. However, traditional numerical solvers are limited by the choice of step sizes when estimating integration, resulting in a trade-off between accuracy and computational efficiency. To address this challenge, we introduce a deep learning-based corrector called Neural Vector (NeurVec), which can compensate for integration errors and enable larger time step sizes in simulations. Our extensive experiments on a variety of complex dynamical system benchmarks demonstrate that NeurVec exhibits remarkable generalization capability on a continuous phase space, even when trained using limited and discrete data. NeurVec significantly accelerates traditional solvers, achieving speeds tens to hundreds of times faster while maintaining high levels of accuracy and stability. Moreover, NeurVec’s simple-yet-effective design, combined with its ease of implementation, has the potential to establish a new paradigm for fast-solving differential equations based on deep learning.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

From Femtoseconds to Gigaseconds: The SolDeg Platform for the Performance Degradation Analysis of Silicon Heterojunction Solar Cells

Heterojunction Si solar cells exhibit notable performance degradation. Here, we modeled this degradation by electronic defects getting generated by thermal activation across energy barriers over time. To analyze the physics of this degradation, we developed the SolDeg platform to simulate the dynamics of electronic defect generation. First, femtosecond molecular dynamics simulations were performed to create a-Si/c-Si stacks, using the machine learning-based Gaussian approximation potential. Second, we created shocked clusters by a cluster blaster method. Third, the shocked clusters were analyzed to identify which of them supported electronic defects. Fourth, the distribution of energy barriers that control the generation of these electronic defects was determined. Fifth, an accelerated Monte Carlo method was developed to simulate the thermally activated time-dependent defect generation across the barriers. Our main conclusions are as follows. (1) The degradation of a-Si/c-Si heterojunction solar cells via defect generation is controlled by a broad distribution of energy barriers. (2) We developed the SolDeg platform to track the microscopic dynamics of defect generation across this wide barrier distribution and determined the time-dependent defect density N(t) from femtoseconds to gigaseconds, over 24 orders of magnitude in time. (3) We have shown that a stretched exponential analytical form can successfully describe the defect generation N(t) over at least 10 orders of magnitude in time. (4) We found that in relative terms, Voc degrades at a rate of 0.2%/year over the first year, slowing with advancing time. (5) We developed the time correspondence curve to calibrate and validate the accelerated testing of solar cells. We found a compellingly simple scaling relationship between accelerated and normal times t normal ∝ t accel T(accel)/T(normal) . (6) We also carried out experimental studies of defect generation in a-Si:H/c-Si stacks. We found a relatively high degradation rate at early times that slowed considerably at longer time scales.

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

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl 3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.’s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U param eter for the 5$f$ electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl 3 ) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

HPSim, Version 1.0

HPSim is a GPU-accelerated online multi-particle beam dynamics simulation tool for ion linacs. It was originally developed for use on the Los Alamos 800-MeV proton linac. It is a “z-code” that contains typical linac beam transport elements. The linac RF-gap transformation utilizes transit-time-factors to calculate the beam acceleration therein. The space-charge effects are computed using the 2D SCHEFF (Space CHarge EFFect) algorithm, which calculates the radial and longitudinal space charge forces for cylindrically symmetric beam distributions. Other space- charge routines to be incorporated include the 3D PICNIC and a 3D Poisson solver. HPSim can simulate beam dynamics in drift tube linacs (DTLs) and coupled cavity linacs (CCLs). Elliptical superconducting cavity (SC) structures will also be incorporated into the code. The computational core of the code is written in C++ and accelerated using the NVIDIA CUDA technology. Users access the core code, which is wrapped in Python/C APIs, via Pythons scripts that enable ease-of-use and automation of the simulations. The overall linac description including the EPICS PV machine control parameters is kept in an SQLite database that also contains calibration and conversion factors required to transform the machine set points into model values used in the simulation.

Pang, Xiaoying↗
DOE OSTI2016

HPSim, Version 1.0

HPSim is a GPU-accelerated online multi-particle beam dynamics simulation tool for ion linacs. It was originally developed for use on the Los Alamos 800-MeV proton linac. It is a “z-code” that contains typical linac beam transport elements. The linac RF-gap transformation utilizes transit-time-factors to calculate the beam acceleration therein. The space-charge effects are computed using the 2D SCHEFF (Space CHarge EFFect) algorithm, which calculates the radial and longitudinal space charge forces for cylindrically symmetric beam distributions. Other space- charge routines to be incorporated include the 3D PICNIC and a 3D Poisson solver. HPSim can simulate beam dynamics in drift tube linacs (DTLs) and coupled cavity linacs (CCLs). Elliptical superconducting cavity (SC) structures will also be incorporated into the code. The computational core of the code is written in C++ and accelerated using the NVIDIA CUDA technology. Users access the core code, which is wrapped in Python/C APIs, via Pythons scripts that enable ease-of-use and automation of the simulations. The overall linac description including the EPICS PV machine control parameters is kept in an SQLite database that also contains calibration and conversion factors required to transform the machine set points into model values used in the simulation.

Pang, Xiaoying↗
DOE OSTI2023

High-temperature active oxidation of nanocrystalline silicon-carbide: A reactive force-field molecular dynamics study

Flexible woven SiC ceramics are prone to accelerated fiber embrittlement under high temperature oxidation in dynamic oxygen environments. The nanocrystalline structure of the constituent fibers impacts the reaction kinetics and phase transformations during active oxidation. However, fundamental understanding and quantification of grain boundary effects on oxidation behavior in nanocrystalline SiC remain elusive when temperatures exceed 1500 K. This study deploys large-scale molecular dynamics simulations with a reactive force-field to elucidate the complex roles of atomic oxygen reservoir conditions and grain size on oxidation kinetics and the nature of oxides produced in both monocrystalline and nanocrystalline 3C-SiC between 1100 K and 2000 K. The simulations with dynamically replenished oxygen provide good agreement with oxidation kinetics and activation energies for the monocrystalline Si(100) and C(100) orientations published in the available literature. This study reveals that, by contrast, nanocrystalline SiC samples exhibit two distinct oxidation kinetics with a transition point at 1500 K due to surface melting, which is supported by experimental evidence. The introduction of a grain-boundary network produces a two-fold decrease in oxidation activation energies compared to monocrystalline SiC below 1500 K. Above 1500 K, however, the activation energies rise substantially due to the formation of a liquid Si phase at the SiC/Si oxide interface. Further, it is shown that the stability of the interfacial liquid phase is promoted by incoherent grain boundaries in the crystalline SiC. These findings are important for the deployment of nanocrystalline SiC fibers in advanced thermal protection systems for high-temperature applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular Insight into the Effects of Clustering on the Dynamics of Ionomers in Solutions

Ionizable groups tethered to polymers enable their many current and potential applications. However, these functionalities drive the formation of physical networks through clustering of the ionic groups, resulting in constrained dynamics of the macromolecules. Understanding the molecular origin of this hindrance remains a critical fundamental question, whose solution will directly impact the processing of ionizable polymers from molecules to viable materials. Here, in this work, using quasielastic neutron scattering accompanied by molecular dynamics simulations, segmental dynamics of slightly sulfonated polystyrene is studied in solutions as the cohesion of the ionic assemblies is tuned. We find that in cyclohexane the ionic assemblies act as centers of confinement, affecting dynamics both on macroscopic lengths and in the vicinity of the ionic assemblies. Addition of a small amount of ethanol affects the packing of the ionizable groups within the assemblies, which in turn enhances the chain dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Learning-Based Load Control to Support Resilient Networked Microgrid Operations

Microgrids have proven to be an effective option for increasing the resiliency of critical end-use loads during extreme events. Building on past operational experiences, some microgrid operators are examining the potential to network microgrids to further improve resiliency. However, the frequency deviations experienced on isolated microgrids during transient events, such as switching operations, step increases in load, and loss of generation, are significantly larger than those typically seen on bulk transmission systems. The larger frequency deviations can cause a loss of inverter-connected assets, resulting in a loss of power to critical end-use loads. This paper presents a method of mitigating the impact of transient events by engaging end-use loads using Grid-Friendly Appliance TM (GFA) controllers. An online, i.e., real-time, device-level algorithm is presented, which adjusts individual GFA controller frequency set-points based on the operational characteristics of each end-use load, and on the changing grid dynamic characteristics. The presented method improves the dynamic stability of the networked microgrid operations while minimizing the interruptions to end-use loads. The presented work is validated with dynamic simulations using a modified version of the IEEE 123-node test system with three microgrids, using the GridLAB-D TM simulation environment.

Radhakrishnan, Nikitha↗
DOE OSTI2016

Coherent dynamic structure factors of strongly coupled plasmas: A generalized hydrodynamic approach

A generalized hydrodynamic fluctuation model is proposed to simplify the calculation of the dynamic structure factor S(ω, k) of non-ideal plasmas using the fluctuation-dissipation theorem. In this model, the kinetic and correlation effects are both included in hydrodynamic coefficients, which are considered as functions of the coupling strength (Γ) and collision parameter (kλ{sub ei}), where λ{sub ei} is the electron-ion mean free path. A particle-particle particle-mesh molecular dynamics simulation code is also developed to simulate the dynamic structure factors, which are used to benchmark the calculation of our model. A good agreement between the two different approaches confirms the reliability of our model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Dynamical Structure Factor of the Three-Dimensional Quantum Spin Liquid Candidate NaCaNi 2 F 7

We study the dynamical structure factor of the spin-1 pyrochlore material NaCaN 2 F 7 , which is well described by a weakly perturbed nearest-neighbour Heisenberg Hamiltonian, Our three approaches—molecular dynamics simulations, stochastic dynamical theory, and linear spin wave theory—reproduce remarkably well the momentum dependence of the experimental inelastic neutron scattering intensity as well as its energy dependence with the exception of the lowest energies. Here, we discuss two surprising aspects and their implications for quantum spin liquids in general: the complete lack of sharp quasiparticle excitations in momentum space and the success of the linear spin wave theory in a regime where it would be expected to fail for several reasons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Computing Free Energies with Fluctuation Relations on Quantum Computers

As a central thermodynamic property, free energy enables the calculation of virtually any equilibrium property of a physical system, allowing for the construction of phase diagrams and predictions about transport, chemical reactions, and biological processes. Thus, methods for efficiently computing free energies, which in general is a difficult problem, are of great interest to broad areas of physics and the natural sciences. The majority of techniques for computing free energies target classical systems, leaving the computation of free energies in quantum systems less explored. Recently developed fluctuation relations enable the computation of free energy differences in quantum systems from an ensemble of dynamic simulations. While performing such simulations is exponentially hard on classical computers, quantum computers can efficiently simulate the dynamics of quantum systems. Here, we present an algorithm utilizing a fluctuation relation known as the Jarzynski equality to approximate free energy differences of quantum systems on a quantum computer. In this work we discuss under which conditions our approximation becomes exact, and under which conditions it serves as a strict upper bound. Furthermore, we successfully demonstrate a proof of concept of our algorithm using the transverse field Ising model on a real quantum processor. As quantum hardware continues to improve, we anticipate that our algorithm will enable computation of free energy differences for a wide range of quantum systems useful across the natural sciences.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Load Margin Constrained Moving Target Defense against False Data Injection Attacks

Cyber physical security of power systems with high penetration of renewable generation has attracted attention from researchers. One critical issue is that cyber-physical attacks, disguised as uncertain renewable generation, can target conventional power system state estimation (SE). Moving target defense (MTD) is a promising defense strategy to detect stealthy false data injection (FDI) attacks against SE. However, all existing studies myopically perturb the reactance of transmission lines equipped with distributed flexible AC transmission system (D-FACTS) devices without adequately considering the system voltage stability. Exacerbated by the renewable generation uncertainty, existing MTD may cause voltage instability when the power grid is under stress. To address this issue, we propose a novel MTD framework that explicitly considers system voltage stability by using continuation power flow. We utilize the sensitivity matrix of power injection to line impedance, on which an optimization problem for maximizing load margin is formulated. This framework is validated on the IEEE 14-bus system and the IEEE 118-bus system, in which net load redistribution attacks are launched by sophisticated attackers. Steady-state simulations and dynamic simulations on PSS/E show the effectiveness of the proposed framework in circumventing the voltage instability while maintaining the detection effectiveness of MTD. The impact of the proposed method on attack detection effectiveness is also revealed.

Zhang, Hang↗
DOE OSTI2022

Development of Large Scale Extrusion Deposition for Structural Applications

Large Scale Extrusion Deposition (LSED) is an evolving additive manufacturing (AM) technology that research, such as that taking place at Oak Ridge National Laboratory (ORNL) and companies like Local Motors, are continuing to utilize and develop new applications. A major LSED application of interest has been molds and tooling as it allows for much shorter production time and lower cost. However, interest has been expanding into using LSED for more structural applications due to more frequent use of high performing polymer composites as feedstock. The use of LSED for structural applications is of particular interest to Local Motors as it is currently being used to create commercially viable, energy efficient electric vehicles. LSED offers a unique opportunity when compared to traditional manufacturing methods as it can significantly reduce the number of components necessary, while decreasing embodied energy and carbon emissions. Even with these benefits however, it is important that LSED is properly understood from a structural aspect as this field has not been as heavily researched as tooling. To ensure a high level of safety and repeatability is an essential responsibility of a manufacturer whose products are structural in nature. As such it is important to understand the materials and LSED process to make structural objects that the manufacturer can be confident in. The main goals of this project were to: develop and investigate materials that are of interest for structural LSED applications, further develop and understand the current machines used in LSED and develop simulations tools of LSED and the mechanical properties of the created structure. Material development focused on composite materials that have high mechanical properties and are stable in a variety of environments. The materials were tested to determine their as printed mechanical and thermal properties as these are necessary for simulations. Once the materials were investigated thoroughly, it allowed for simulations to be performed to compare to experimental data with simulations. Materials were also vetted to determine candidates for multi-material printing. The machine development focused around the areas of process monitoring, non-destructive evaluation, and quality control. Finally, the goal of the simulations was to develop a realistic model of printed structures, including in-process simulation and dynamic simulation. The routes to get to some of these goals and the depth in which they were investigated changed throughout the project due to personnel changes and the COVID-19 pandemic. This project resulted in many valuable results such as the development of an nondestructive evaluation (NDE) technique for interlayer defects, proof of simulation for warpage in simple parts, the development and utilization of a profilometer to monitor a print for defects or inconsistencies, thorough investigation of a material used commercially for structural LSED applications and valuable experimental data on the applicability and advantages of multi-material crush structures versus their neat counterparts by creation and testing of samples.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Space life support engineering program

A comprehensive study to develop software to simulate the dynamic operation of water reclamation systems in long-term closed-loop life support systems is being carried out as part of an overall program for the design of systems for a moon station or a Mars voyage. This project is being done in parallel with a similar effort in the Department of Chemistry to develop durable accurate low-cost sensors for monitoring of trace chemical and biological species in recycled water supplies. Aspen-Plus software is being used on a group of high-performance work stations to develop the steady state descriptions for a number of existing technologies. Following completion, a dynamic simulation package will be developed for determining the response of such systems to changes in the metabolic needs of the crew and to upsets in system hardware performance.

Seagrave, Richard C.↗