Nano-Rover Solar Sail Dynamic Simulation
Explore the source record for details and available documents.
Solar Sail ADAMS simulation dynamic↗
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.
Nanorover Solar Sail Dynamic Simulation
The Nanorover Solar Sail is a concept being evaluated for the deployment of small rovers that may be used for extra-terrestial exploration.
solar sail ADAMS simulation dynamic↗
Adaptive variational quantum dynamics simulations with compressed circuits and fewer measurements
The adaptive variational quantum dynamics simulation (AVQDS) method performs real-time evolution of quantum states using automatically generated parametrized quantum circuits that often contain substantially fewer gates than Trotter circuits. Here we report an improved version of the method, which we call AVQDS(T), by porting the tiling efficient trial circuits with rotations implemented simultaneously technique. The algorithm adaptively adds layers of disjoint unitary gates to the ansatz circuit so as to keep the McLachlan distance, a measure of the accuracy of the variational dynamics, below a fixed threshold. Here we perform benchmark noiseless AVQDS(T) simulations of quench dynamics in local spin models and compare with an alternative adaptive variational approach on quantum resource requirement. Quantum dynamical simulations implementing realistic noise channels are also reported. Finally, we propose a way to substantially alleviate the measurement overhead of AVQDS(T) while maintaining high accuracy by synergistically integrating quantum circuit calculations on quantum processing units with classical calculations using, e.g., tensor networks to evaluate the quantum geometric tensor. We showcase that this approach enables AVQDS(T) to deliver more accurate results than simulations using a fixed ansatz of comparable final depth for a significant time duration with fewer quantum resources.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Requirements specification part 1 and analysis of dynamic simulation methods for launch vehicle component level simulation
Requirements specification and analysis of dynamic simulation methods for launch vehicle component level simulation computer program
Source record↗
Translocation of a Polymer Chain across a Nanopore: A Brownian Dynamics Simulation Study
We carried out Brownian dynamics simulation studies of the translocation of single polymer chains across a nanosized pore under the driving of an applied field (chemical potential gradient). The translocation process can be either dominated by the entropic barrier resulted from restricted motion of flexible polymer chains or by applied forces (or chemical gradient across the wall), we focused on the latter case in our studies. Calculation of radius of gyrations at the two opposite sides of the wall shows that the polymer chains are not in equilibrium during the translocation process. Despite this fact, our results show that the one-dimensional diffusion and the nucleation model provide an excellent description of the dependence of average translocation time on the chemical potential gradients, the polymer chain length and the solvent viscosity. In good agreement with experimental results and theoretical predictions, the translocation time distribution of our simple model shows strong non-Gaussian characteristics. It is observed that even for this simple tubelike pore geometry, more than one peak of translocation time distribution can be generated for proper pore diameter and applied field strengths. Both repulsive Weeks-Chandler-Anderson and attractive Lennard-Jones polymer-nanopore interaction were studied, attraction facilitates the translocation process by shortening the total translocation time and dramatically improve the capturing of polymer chain. The width of the translocation time distribution was found to decrease with increasing temperature, increasing field strength, and decreasing pore diameter.
Tian, Pu↗
Ab initio molecular dynamics simulations of AlN responding to low energy particle radiation
Ab initio molecular dynamics simulations of low energy recoil events in wurtzite AlN have been performed to determine threshold displacement energies, defect production and evolution mechanisms, role of partial charge transfer during the process, and the influence of irradiation-induced defects on the properties of AlN. Here, the results show that the threshold displacement energies, E d , along the direction parallel to the basal planes are smaller than those perpendicular to the basal planes. The minimum E d values are determined to be 19 eV and 55 eV for N and Al atom, respectively, which occur along the [$\overline{11}20$] direction. In general, the threshold displacement energies for N are smaller than those for Al atom, indicating the N defects would be dominant under irradiation. The defect production mechanisms have been analyzed. It is found that charge transfer and redistribution for both the primary knock-on atom and the subsequent recoil atoms play a significant role in defect production and evolution. Similar to the trend in oxide materials, there is a nearly linear relationship between E d and the total amount of charge transfer at the potential energy peak in AlN, which provides guidance on the development of charge-transfer interatomic potentials for classic molecular dynamics simulations. Finally, the response behavior of AlN to low energy irradiation is qualitatively investigated. The existence of irradiation-induced defects significantly modifies the electronic structure, and thus affects the magnetic, electronic and optical properties of AlN. In conclusion, these findings further enrich the understanding of defects in the wide bandgap semiconductor of AlN.
36 MATERIALS SCIENCE↗
Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method
Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. Here, the nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Notes on implementation of Coulomb friction in coupled dynamical simulations
A coupled dynamical system is defined as an assembly of rigid/flexible bodies that may be coupled by kinematic connections. The interfaces between bodies are modeled using hinges having 0 to 6 degrees of freedom. The equations of motion are presented for a mechanical system of n flexible bodies in a topological tree configuration. The Lagrange form of the D'Alembert principle was employed to derive the equations. The equations of motion are augmented by the kinematic constraint equations. This augmentation is accomplished via the method of singular value decomposition.
Vandervoort, R. J.↗
Validation of Cryogenic Propellant Tank Filling using Computational Fluid Dynamics Simulation
The Fluid Dynamics Branch at MSFC has positioned itself to support a wide range of customers in need of Cryogenic Fluid Management (CFM) analysis. A computational fluid dynamics (CFD) tool used for all manner of internal and external propulsion applications has been extended and refined to better model cryogenic propellant storage and tanking operations. Through the CFM Portfolio project, several validation activities were initiated. Validation of propellant tank self-pressurization, autogenous pressurization, slosh-induced ullage collapse, and jet-induced mixing all aid in defining model accuracy. The on-going validation effort has enabled confident application of the tool to in-line design and evaluation of CFM hardware and operations. Recent project support included defining the impact of in-space slosh dynamics on reaction control system mass for Space Launch System (SLS) upper stages. Propellant mixing strategies were defined to improve performance of a thermal vent system for a Commercial Lunar Payload Services (CLPS) partner. Design support of in-space maneuvers, tank hardware, and autogenous pressurization operations was also provided through Human Landing System (HLS) collaboration work. The branch has engaged the CFM community to share recent findings and capabilities through several forums including conferences, technical interchange meetings, and workshops. Development and demonstration of CFM modeling capabilities continues in this work on the no-vent fill of propellant tank in micro-gravity to meet the needs of NASA and its industry partners in the endeavor to sustainably reach the Moon and beyond.
CFD↗
Kernel fusion in atomistic spin dynamics simulations on Nvidia GPUs using tensor core
In atomistic spin dynamics simulations, the time cost of constructing the space- and time-displaced pair correlation function in real space increases quadratically as the number of spins N, leading to significant computational effort. The GEMM subroutine can be adopted to accelerate the calculation of the dynamical spin-spin correlation function, but the computational cost of simulating large spin systems (>40000 spins) on CPUs remains expensive. In this work, we perform the simulation on the graphics processing unit (GPU), a hardware solution widely used as an accelerator for scientific computing and deep learning. Here we show that GPUs can accelerate the simulation up to 25-fold compared to multi-core CPUs when using the GEMM subroutine on both. To hide memory latency, we fuse the element-wise operation into the GEMM kernel using CUTLASS that can improve the performance by 26% ~ 33% compared to implementation based on cuBLAS. Furthermore, we perform the on-the-fly calculation in the epilogue of the GEMM subroutine to avoid saving intermediate results on global memory, which makes the large-scale atomistic spin dynamics simulation feasible and affordable.
97 MATHEMATICS AND COMPUTING↗
Method of information entropy for convergence assessment of molecular dynamics simulations
The lack of a reliable method to evaluate the convergence of molecular dynamics simulations has contributed to discrepancies in different areas of molecular dynamics. Here, the method of information entropy is introduced to molecular dynamics for stationarity assessment. The Shannon information entropy formalism is used to monitor the convergence of the atom motion to a steady state in a continuous spatial domain and is also used to assess the stationarity of calculated multidimensional fields such as the temperature field in a discrete spatial domain. It is demonstrated in this work that monitoring the information entropy of the atom position matrix provides a clear indicator of reaching steady state in radiation damage simulations, non-equilibrium molecular dynamics thermal conductivity computations, and simulations of Poiseuille and Couette flow in nanochannels. A main advantage of the present technique is that it is non-local and relies on fundamental quantities available in all molecular dynamics simulations. Unlike monitoring average temperature, the technique is applicable to simulations that conserve total energy such as reverse non-equilibrium molecular dynamics thermal conductivity computations and to simulations where energy dissipates through a boundary as in radiation damage simulations. The method is applied to simulations of iron using the Tersoff/ZBL splined potential, silicon using the Stillinger–Weber potential, and to Lennard–Jones fluid. Its applicability to both solids and fluids shows that the technique has potential for generalization to other areas in molecular dynamics.
74 ATOMIC AND MOLECULAR PHYSICS↗
Efficient fully-coherent quantum signal processing algorithms for real-time dynamics simulation
Simulating the unitary dynamics of a quantum system is a fundamental problem of quantum mechanics, in which quantum computers are believed to have significant advantage over their classical counterparts. One prominent such instance is the simulation of electronic dynamics, which plays an essential role in chemical reactions, non-equilibrium dynamics, and material design. These systems are time-dependent, which requires that the corresponding simulation algorithm can be successfully concatenated with itself over different time intervals to reproduce the overall coherent quantum dynamics of the system. In this paper, we quantify such simulation algorithms by the property of being fully-coherent: the algorithm succeeds with arbitrarily high success probability 1 − δ while only requiring a single copy of the initial state. Here we subsequently develop fully-coherent simulation algorithms based on quantum signal processing (QSP), including a novel algorithm that circumvents the use of amplitude amplification while also achieving a query complexity additive in time t, ln(1/δ), and ln(1/ϵ) for error tolerance ϵ: $Θ‖\mathscr{H}‖|t|+ln(1/ϵ)+ln(1/δ)$. Furthermore, we numerically analyze these algorithms by applying them to the simulation of the spin dynamics of the Heisenberg model and the correlated electronic dynamics of an H2 molecule. Since any electronic Hamiltonian can be mapped to a spin Hamiltonian, our algorithm can efficiently simulate time-dependent ab initio electronic dynamics in the circuit model of quantum computation. Accordingly, it is also our hope that the present work serves as a bridge between QSP-based quantum algorithms and chemical dynamics, stimulating a cross-fertilization between these exciting fields.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Reactive Molecular Dynamics Simulation of Kerogen Thermal Maturation and Cross-Linking Pathways
For this study, molecular dynamics simulations were performed with a ReaxFF reactive force field to investigate bond breaking and bond formation mechanisms during the early thermal maturation of kerogen and potential crosslinking pathways towards a three-dimensional (3D) quasi-infinite molecular network (crosslinked kerogen). Starting with small ensembles of high molecular mass models for immature type I Green River Shale kerogen (Kerogen 1-I), top of the oil window type II kerogen (Kerogen 2-L), and low maturity type III kerogen (Kerogen 3-L), low molecular mass species including H 2 O, H 2 S, C 2 H 4 , C 3 H 6 were produced as the maturities of the remaining kerogens increased. Highly reactive fragments, which are not detected in pyrolysis experiments, were also produced. Further, the crosslinking mechanism in the newly develop polymeric kerogen networks appears to be highly complex and it is influenced by the nature of the kerogen. Covalent —C—S—, —C—O— and —C—C— bonds are the primary crosslinks that structurally binds kerogen monomers together. The trends observed in the thermochemical transformation of kerogen and the kerogen crosslinking pathways are consistent with the theoretical and experimental studies reported in the scientific literature. The reactive force field molecular dynamics simulation provides a potentially valuable alternative approach to the development of a realistic 3D molecular network of kerogen. However, the conversion of kerogen molecules into a 3D molecular network was low.
02 PETROLEUM↗
The selective heating effect of microwave irradiation on a binary mixture of water and polyethylene oxide: a molecular dynamics simulation approach
In this study, we investigate the molecular mechanisms of a microwave-driven selective heating process by performing molecular dynamics simulations for three different systems including pure water, pure polyethylene oxide (PEO), and water-PEO mixed systems in the presence of a microwave with two different intensities of electric field such as 0.001 V/Å -1 and 0.01 V/Å -1 at a frequency of 100 GHz. First, from performing molecular dynamics simulations of CO and CO 2 in the presence of the microwave, it is confirmed that the molecular dipole moment is responsible for the rotational motion induced by the oscillating electric field. Second, by analyzing the MD simulations of the pure water system, we discover that the dipole moment of water exhibits a time lag with respect to the microwave. During the heating process, however, the temperature, kinetic, and potential energies increase synchronously with the oscillating electric field of the microwave, showing that the heating of the water system is caused by the molecular reaction of water to the microwave. Comparing the water-PEO mixed system to the pure water and pure PEO systems, the water-PEO mixed system has a higher heating rate than the pure PEO system but a lower heating rate than the pure water system. Therefore, we conclude that heating the water-PEO mixed system is driven by water molecules selectively activated by microwave irradiation. Furthermore, we also calculate the diffusion coefficients of water molecules and PEO chains by describing their mean square displacements, demonstrating that the diffusion coefficients are increased in the presence of microwaves for both water and PEO in pure and mixed systems. Lastly, during the microwave heating process, the structures of the water-PEO mixed system are altered as a function of the intensity of electric field, which is mainly driven by the response of water molecules.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Ensemble Docking in Drug Discovery: How Many Protein Configurations from Molecular Dynamics Simulations are Needed To Reproduce Known Ligand Binding?
Ensemble docking in drug discovery or chemical biology uses dynamical simulations of target proteins to generate binding site conformations for docking campaigns. We demonstrate that 600 ns molecular dynamics simulations of four G-protein-coupled receptors in their membrane environments generate ensembles of protein configurations that, collectively, are selected by 70–99% of the known ligands of these proteins. Thus, the process of ligand recognition by conformational selection can be reproduced by combining molecular dynamics and docking calculations. Clustering of the molecular dynamics trajectories, however, does not necessarily identify the protein conformations that are most often selected by the ligands.
60 APPLIED LIFE SCIENCES↗
Synthesis of formamide and related organic species in the interstellar medium via chemical dynamics simulations
We show, by means of direct dynamics simulations, how it is possible to define possible reactants and mechanisms leading to the formation of formamide in the interstellar medium. In particular, different ion–molecule reactions in the gas phase were considered: NH{sub 3}OH{sup +}, NH{sub 2}OH{sub 2}{sup +}, H{sub 2}COH{sup +}, and NH{sub 4} {sup +} for the ions and NH{sub 2}OH, H{sub 2}CO, and NH{sub 3} for the partner neutrals. These calculations were combined with high level ab initio calculations to investigate possible further evolution of the products observed. In particular, for formamide, we propose that the NH{sub 2}OH{sub 2}{sup +} + H{sub 2}CO reaction can produce an isomer, NH{sub 2}OCH{sub 2}{sup +}, that, after dissociative recombination, can produce neutral formamide, which was observed in space. The direct dynamics do not pre-impose any reaction pathways and in other reactions, we did not observe the formation of formamide or any possible precursor. On the other hand, we obtained other interesting reactions, like the formation of NH{sub 2}CH{sub 2}{sup +}. Finally, some radiative association processes are proposed. All of the results obtained are discussed in light of the species observed in radioastronomy.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Use of a large jet transport as an inflight dynamic simulator
Modified jet transport used as inflight dynamic simulator for supersonic transport low speed approach and landing characteristics
FLIGHT SIMULATOR↗
Staggered solution procedures for multibody dynamics simulation
The numerical solution procedure for multibody dynamics (MBD) systems is termed a staggered MBD solution procedure that solves the generalized coordinates in a separate module from that for the constraint force. This requires a reformulation of the constraint conditions so that the constraint forces can also be integrated in time. A major advantage of such a partitioned solution procedure is that additional analysis capabilities such as active controller and design optimization modules can be easily interfaced without embedding them into a monolithic program. After introducing the basic equations of motion for MBD system in the second section, Section 3 briefly reviews some constraint handling techniques and introduces the staggered stabilized technique for the solution of the constraint forces as independent variables. The numerical direct time integration of the equations of motion is described in Section 4. As accurate damping treatment is important for the dynamics of space structures, we have employed the central difference method and the mid-point form of the trapezoidal rule since they engender no numerical damping. This is in contrast to the current practice in dynamic simulations of ground vehicles by employing a set of backward difference formulas. First, the equations of motion are partitioned according to the translational and the rotational coordinates. This sets the stage for an efficient treatment of the rotational motions via the singularity-free Euler parameters. The resulting partitioned equations of motion are then integrated via a two-stage explicit stabilized algorithm for updating both the translational coordinates and angular velocities. Once the angular velocities are obtained, the angular orientations are updated via the mid-point implicit formula employing the Euler parameters. When the two algorithms, namely, the two-stage explicit algorithm for the generalized coordinates and the implicit staggered procedure for the constraint Lagrange multipliers, are brought together in a staggered manner, they constitute a staggered explicit-implicit procedure which is summarized in Section 5. Section 6 presents some example problems and discussions concerning several salient features of the staggered MBD solution procedure are offered in Section 7.
Park, K. C.↗