SEARCH · Search NASA

Results for “Dynamical Simulations”

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 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Radiation damage in cubic ZrO 2 and yttria-stabilized zirconia from molecular dynamics simulations

Here, we perform molecular dynamics simulation on cubic ZrO 2 and yttria-stabilized zirconia (YSZ) to elucidate defect cluster formation resulting from radiation damage, and evaluate the impact of Y-dopants. Interstitial clusters composed of split-interstitial building blocks, i.e., Zr-Zr or Y-Zr are formed. Moreover, oxygen vacancies control cation defect migration; in their presence, Zr interstitials aggregate to form split-interstitials whereas in their absence Zr interstitials remain immobile, as isolated single-interstitials. Y-doping prevents interstitial cluster formation due to sequestration of oxygen vacancies.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Analysis of a Li-Ion Nanobattery with Graphite Anode Using Molecular Dynamics Simulations

In this work, molecular dynamics simulations were performed of the initial charging of a Li-ion nanobattery with a graphite anode and lithium hexaflourphosphate (LiPF 6 ) salt dissolved in ethylene carbonate (CO 3 C 2 H 4 ) solvent as the electrolyte solution. The charging was achieved through the application of external electric fields simulating voltage sources. A variety of force fields were combined to simulate the materials of the nanobattery, including the solid electrolyte interphase, metal collectors, and insulator cover. Some of the force field parameters were estimated using ab initio methods and others were taken from the literature. We studied the behavior of Li-ions traveling from cathode to anode through electrolyte solutions of concentrations 1.15 and 3.36 M. Time-dependent variables such as energy, temperature, volume, polarization, and mean square displacement are reported; a few of these variables, as well as others such as current, resistance, current density, conductivity, and resistivity are reported as a function of the external field and charging voltage. A solid electrolyte interphase (SEI) layer was also added to the model to study the mechanism behind the diffusion of the Li-ions through the SEI. As the battery is charged, the depletion of Li atoms in the cathode and their accumulation in the anode follow a linear increase of the polarizability in the solvent, until reaching a saturation point after which the charging of the battery stops, i.e., the energy provided by the external source decays to very low levels. Lastly, the nanobattery model containing the most common materials of a commercial lithium-ion battery is very useful to determine atomistic information that is difficult or too expensive to obtain experimentally; available data shows consistency with our results.

25 ENERGY STORAGE↗

NASA NTRS1975

Dispersion analysis techniques within the space vehicle dynamics simulation program

The Space Vehicle Dynamics Simulation (SVDS) program was evaluated as a dispersion analysis tool. The Linear Error Analysis (LEA) post processor was examined in detail and simulation techniques relative to conducting a dispersion analysis using the SVDS were considered. The LEA processor is a tool for correlating trajectory dispersion data developed by simulating 3 sigma uncertainties as single error source cases. The processor combines trajectory and performance deviations by a root-sum-square (RSS process) and develops a covariance matrix for the deviations. Results are used in dispersion analyses for the baseline reference and orbiter flight test missions. As a part of this study, LEA results were verified as follows: (A) Hand calculating the RSS data and the elements of the covariance matrix for comparison with the LEA processor computed data. (B) Comparing results with previous error analyses. The LEA comparisons and verification are made at main engine cutoff (MECO).

Snow, L. S.↗

NASA NTRS2007

Local Elastic Constants for Epoxy-Nanotube Composites from Molecular Dynamics Simulation

A method from molecular dynamics simulation is developed for determining local elastic constants of an epoxy/nanotube composite. The local values of C11, C33, K12, and K13 elastic constants are calculated for an epoxy/nanotube composite as a function of radial distance from the nanotube. While the results possess a significant amount of statistical uncertainty resulting from both the numerical analysis and the molecular fluctuations during the simulation, the following observations can be made. If the size of the region around the nanotube is increased from shells of 1 to 6 in thickness, then the scatter in the data reduces enough to observe trends. All the elastic constants determined are at a minimum 20 from the center of the nanotube. The C11, C33, and K12 follow similar trends as a function of radial distance from the nanotube. The K13 decreases greater distances from the nanotube and becomes negative which may be a symptom of the statistical averaging.

Frankland, S. J. V.↗

DOE OSTI2015

Classical Molecular Dynamics Simulation of Nuclear Fuel

Molecular dynamics simulation is well suited to study primary damage production by irradiation, defect interactions with fission gas atoms, gas bubble nucleation, grain boundary effects on defect and gas bubble evolution in nuclear fuel, and the resulting changes in thermo-mechanical properties. In these simulations, the forces on the ions are dictated by interaction potentials generated by fitting properties of interest to experimental data. The results obtained from the present generation of potentials are qualitatively similar, but quantitatively different. There is a need to refine existing potentials to provide a better representation of the performance of polycrystalline fuel under a variety of operating conditions, and to develop models that are equipped to handle deviations from stoichiometry. In addition to providing insights into fundamental mechanisms governing the behaviour of nuclear fuel, MD simulations can also provide parameters that can be used as inputs for mesoscale models.

Devanathan, Ram↗

DOE OSTI2024

Dynamical simulation via quantum machine learning with provable generalization

Much attention has been paid to dynamical simulation and quantum machine learning (QML) independently as applications for quantum advantage, while the possibility of using QML to enhance dynamical simulations has not been thoroughly investigated. Here we develop a framework for using QML methods to simulate quantum dynamics on near-term quantum hardware. We use generalization bounds, which bound the error a machine learning model makes on unseen data, to rigorously analyze the training data requirements of an algorithm within this framework. Our algorithm is thus resource efficient in terms of qubit and data requirements. Furthermore, our preliminary numerics for the XY model exhibit efficient scaling with problem size, and we simulate 20 times longer than Trotterization on IBMQ-Bogota. Published by the American Physical Society 2024

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1987

Tether Dynamics Simulation

The proceedings of the conference are presented. The objective was to provide a forum for the discussion of the structure and status of existing computer programs which are used to simulate the dynamics of a variety of tether applications in space. A major topic was different simulation models and the process of validating them. Guidance on future work in these areas was obtained from a panel discussion; the panel was composed of resource and technical managers and dynamic analysts in the tether field. The conclusions of this panel are also presented.

Source record↗

NASA NTRS1989

Controlled multibody dynamics simulation for large space structures

Multibody dynamics discipline, and dynamic simulation in control structure interaction (CSI) design are discussed. The use, capabilities, and architecture of the Large Angle Transient Dynamics (LATDYN) code as a simulation tool are explained. A generic joint body with various types of hinge connections; finite element and element coordinate systems; results of a flexible beam spin-up on a plane; mini-mast deployment; space crane and robotic slewing manipulations; a potential CSI test article; and multibody benchmark experiments are also described.

Housner, J. M.↗

DOE OSTI2016

Multilevel summation with B-spline interpolation for pairwise interactions in molecular dynamics simulations

The multilevel summation method for calculating electrostatic interactions in molecular dynamics simulations constructs an approximation to a pairwise interaction kernel and its gradient, which can be evaluated at a cost that scales linearly with the number of atoms. The method smoothly splits the kernel into a sum of partial kernels of increasing range and decreasing variability with the longer-range parts interpolated from grids of increasing coarseness. Multilevel summation is especially appropriate in the context of dynamics and minimization, because it can produce continuous gradients. This article explores the use of B-splines to increase the accuracy of the multilevel summation method (for nonperiodic boundaries) without incurring additional computation other than a preprocessing step (whose cost also scales linearly). To obtain accurate results efficiently involves technical difficulties, which are overcome by a novel preprocessing algorithm. Numerical experiments demonstrate that the resulting method offers substantial improvements in accuracy and that its performance is competitive with an implementation of the fast multipole method in general and markedly better for Hamiltonian formulations of molecular dynamics. The improvement is great enough to establish multilevel summation as a serious contender for calculating pairwise interactions in molecular dynamics simulations. In particular, the method appears to be uniquely capable for molecular dynamics in two situations, nonperiodic boundary conditions and massively parallel computation, where the fast Fourier transform employed in the particle–mesh Ewald method falls short.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Molecular dynamics simulation studies of 1,3-dimethyl imidazolium nitrate ionic liquid with water

The fundamental understanding of intermolecular interactions of ionic liquids (ILs) with water is essential in predicting IL–water thermodynamic properties. In this study, intermolecular or noncovalent interactions were studied for 1,3-dimethyl imidazolium [DMIM]+ cation and nitrate [NO 3 ] – anion with water, employing quantum mechanics and molecular dynamics simulations. Molecular dynamics simulations were performed using a revised multipolar polarizable force field. The effect of water on ionic liquids was evaluated in terms of thermodynamic and dynamic properties. Thermodynamic properties included liquid densities ρ, excess molar volumes ΔV E , and liquid structures g(r)⁠. Dynamic properties included self-diffusion coefficients D of mixture constituents as a function of water concentration. The density of ionic liquid–water mixtures monotonically decrease with increasing concentration of water. A negative excess volume was obtained for low and high water concentrations, demonstrating strong intermolecular interactions of water with ionic liquid components. Liquid structures of ionic liquid–water mixtures revealed a tendency for anions to interact with cations at shorter intermolecular distances when the water concentration is increased. Diffusion rates were found to increase for all mixture components with increase in the fraction of water. A significant change in the diffusion rate was found at ~0.3 weight fraction of water. However, the water self-diffusion coefficient was dominant at all concentrations. Furthermore, the ratio of water/anion and anion/cation self-diffusion coefficients was found to decrease linearly with increasing concentration of water molecules.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Structured Ionomer Thin Films at Water Interface: Molecular Dynamics Simulation Insight

Controlling the structure and dynamics of thin films of ionizable polymers at water interfaces is critical to their many applications. As the chemical diversity within one polymer is increased, controlling the structure and dynamics of the polymer, which is a key to their use, becomes a challenge. Here molecular dynamics simulations (MD) are used to obtain molecular insight into the structure and dynamics of thin films of one such macromolecule at the interface with water. The polymer consists of an ABCBA topology with randomly sulfonated polystyrene (C), tethered symmetrically to flexible poly(ethylene- r -propylene) blocks (B), and end-capped by a poly( t -butylstyrene) block (A). The compositions of the interfacial and bulk regions of thin films of the ABCBA polymers are followed as a function of exposure time to water. We find that interfacial rearrangements take place where buried ionic segments migrate toward the water interface. The hydrophobic blocks collapse and rearrange to minimize their exposure to water. In conclusion, the water that initially drives interfacial reengagements breaks the ionic clusters within the film, forming a dynamic hydrophilic internal network within the hydrophobic segments.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Extended-Term Dynamic Simulations with High Penetrations of Photovoltaic Generation.

The uncontrolled intermittent availability of renewable energy sources makes integration of such devices into today's grid a challenge. Thus, it is imperative that dynamic simulation tools used to analyze power system performance are able to support systems with high amounts of photovoltaic (PV) generation. Additionally, simulation durations expanding beyond minutes into hours must be supported. This report aims to identify the path forward for dynamic simulation tools to accom- modate these needs by characterizing the properties of power systems (with high PV penetration), analyzing how these properties affect dynamic simulation software, and offering solutions for po- tential problems. We present a study of fixed time step, explicit numerical integration schemes that may be more suitable for these goals, based on identified requirements for simulating high PV penetration systems. We also present the alternative of variable time step integration. To help determine the characteristics of systems with high PV generation, we performed small signal sta- bility studies and time domain simulations of two representative systems. Along with feedback from stakeholders and vendors, we identify the current gaps in power system modeling including fast and slow dynamics and propose a new simulation framework to improve our ability to model and simulate longer-term dynamics.

14 SOLAR ENERGY↗

DOE OSTI2014

Hydration dynamics in water clusters via quantum molecular dynamics simulations

We have investigated the hydration dynamics in size selected water clusters with n = 66, 104, 200, 500, and 1000 water molecules using molecular dynamics simulations. To study the most fundamental aspects of relaxation phenomena in clusters, we choose one of the simplest, still realistic, quantum mechanically treated test solute, an excess electron. The project focuses on the time evolution of the clusters following two processes, electron attachment to neutral equilibrated water clusters and electron detachment from an equilibrated water cluster anion. The relaxation dynamics is significantly different in the two processes, most notably restoring the equilibrium final state is less effective after electron attachment. Nevertheless, in both scenarios only minor cluster size dependence is observed. Significantly different relaxation patterns characterize electron detachment for interior and surface state clusters, interior state clusters relaxing significantly faster. This observation may indicate a potential way to distinguish surface state and interior state water cluster anion isomers experimentally. A comparison of equilibrium and non-equilibrium trajectories suggests that linear response theory breaks down for electron attachment at 200 K, but the results converge to reasonable agreement at higher temperatures. Relaxation following electron detachment clearly belongs to the linear regime. Cluster relaxation was also investigated using two different computational models, one preferring cavity type interior states for the excess electron in bulk water, while the other simulating non-cavity structure. While the cavity model predicts appearance of several different hydrated electron isomers in agreement with experiment, the non-cavity model locates only cluster anions with interior excess electron distribution. The present simulations show that surface isomers computed with the cavity predicting potential show similar dynamical behavior to the interior clusters of the non-cavity type model. Relaxation associated with cavity collapse presents, however, unique dynamical signatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

A high-performance implementation of atomistic spin dynamics simulations on x86 CPUs

Atomistic spin dynamics simulations provide valuable information about the energy spectrum of magnetic materials in different phases, allowing one to identify instabilities and the nature of their excitations. However, the time cost of evaluating the dynamical correlation function S(q,t) increases quadratically as the number of spins N, leading to significant computational effort, making the simulation of large spin systems very challenging. In this work, here, we propose to use a highly optimized general matrix multiply (GEMM) subroutine to calculate the dynamical spin-spin correlation function that can achieve near-optimal hardware utilization. Furthermore, we fuse the element-wise operations in the calculation of S(q,t) into the in-house GEMM kernel, which results in further performance improvements of 44\% - 71\% on several relatively large lattice sizes when compared to the implementation that uses the GEMM subroutine in OpenBLAS, which is the state-of-the-art open source library for Basic Linear Algebra Subroutine (BLAS).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Molecular Dynamics Simulations of Aldol Condensation Catalyzed by Alkylamine-Functionalized Crystalline Silica Surfaces

Molecular dynamics simulations are performed to investigate the cooperatively catalyzed aldol condensation between acetone and 4-nitrobenzaldehyde on alkylamine (or alkylenamine)-grafted silica surfaces, focusing on the mechanism of the catalytic activation of the acetone and 4-nitrobenzaldehyde by the acidic surface silanols followed by the nucleophilic attack of the basic amine functional group toward the activated reactant. From the analysis of the correlations between the catalytically active acid-base sites and reactants, it is concluded that the catalytic cooperativity of the acid-base pair can be affected by two factors: 1) the competition between the silanol and the amine (or enamine) to form a hydrogen bond with a reactant and 2) the flexibility of the alkylamine (or alkylenamine) backbone. Increasing the flexibility of the alkylamine facilitates the nucleophilic attack of the amine on the reactants. From the molecular dynamics simulations, it is found that C 3 propylamine and C 4 butylamine linkers exhibit the highest probability of reaction, which is consistent with the experimental observation that the activity of the aldol reaction on mesoporous silica depends on the length of alkylamine grafted on the silica surface. Furthermore, this simulation work serves as a pioneering study demonstrating how the molecular simulation approach can be successfully employed to investigate the cooperative catalytic activity of such bifunctional acid-base catalysts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1989

The architecture of Newton, a general-purpose dynamics simulator

The architecture for Newton, a general-purpose system for simulating the dynamics of complex physical objects, is described. The system automatically formulates and analyzes equations of motion, and performs automatic modification of this system equations when necessitated by changes in kinematic relationships between objects. Impact and temporary contact are handled, although only using simple models. User-directed influence of simulations is achieved using Newton's module, which can be used to experiment with the control of many-degree-of-freedom articulated objects.

Cremer, James F.↗

NASA NTRS1994

Dynamic Simulation of a Periodic 10 K Sorption Cryocooler

A transient thermal simulation model has been developed to simulate the dynamic performance of a multiple-stage 10 K sorption cryocooler for spacecraft sensor cooling applications that require periodic quick-cooldown (under 2 minutes) , negligible vibration, low power consumption, and long life (5 to 10 years). The model was specifically designed to represent the Brilliant Eyes Ten-Kelvin Sorption Cryocooler Experiment (BETSCE), but it can be adapted to represent other sorption cryocooler systems as well. The model simulates the heat transfer, mass transfer, and thermodynamic processes in the cryostat and the sorbent beds for the entire refrigeration cycle, and includes the transient effects of variable hydrogen supply pressures due to expansion and overflow of hydrogen during the cooldown operation. The paper describes model limitations and simplifying assumptions, with estimates of errors induced by them, and presents comparisons of performance predictions with ground experiments. An important benefit of the model is its ability to predict performance sensitivities to variations of key design and operational parameters. The insights thus obtained are expected to lead to higher efficiencies and lower weights for future designs.

cryocooler design performance sorption simulation ↗

NASA NTRS1993

Fast Auroral Snapshot performance using a multi-body dynamic simulation

This paper examines the complex dynamic interaction between two 2.6 m long stacer booms, four 30 m long flexible wire booms and the attitude control system of the Fast Auroral SnapshoT (FAST) spacecraft. The FAST vehicle will nominally operate as a negative orbit spinner, positioned in a 83 deg inclination, 350 x 4200 km orbit. For this study, a three-axis, non-linear, seven body dynamic simulation is developed using the TREETOPS software package. The significance of this approach is the ability to model each component of the FAST spacecraft as an individual member and connect them together in order to better understand the dynamic coupling between structures and the control system. Both the wire and stacer booms are modeled as separate bodies attached to a rigid central body. The wire booms are oriented perpendicular to the spin axis at right angles relative to each other, whereas the stacer booms are aligned with the spin axis. The analysis consists of a comparison between the simulated in-plane and out-of-plane boom motions with theoretically derived frequencies, and an examination of the dynamic coupling between the control system and boom oscillations. Results show that boom oscillations of up to 0.36 deg are acceptable in order to meet the performance requirements. The dynamic motion is well behaved when the precession coil is operating, however, activation of the spin coil produces an erratic trend in the spin rate which approaches the spin rate requirement.

Zimbelman, Darrell↗