Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic 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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Excited Electron Dynamics at Semiconductor–Molecule Type-II Heterojunction Interface: First-Principles Dynamics Simulation

Excited electron dynamics at semiconductor–molecule interfaces is ubiquitous in various energy conversion technologies. However, a quantitative understanding of how molecular details influence the quantum dynamics of excited electrons remains a great scientific challenge because of the complex interplay of different processes with various time scales. Here in this paper, we employ first-principles electron dynamics simulations to investigate how molecular features govern the dynamics in a representative interface between the hydrogen-terminated Si(111) surface and a cyanidin molecule. Hot electron transfer to the chemisorbed molecule was observed but was short-lived on the molecule. Interfacial electron transfer to the chemisorbed molecule was found to be largely decoupled from hot electron relaxation within the semiconductor surface. While the hot electron relaxation was found to take place on a time scale of several hundred femtoseconds, the subsequent interfacial electron transfer was slower by an order of magnitude. Finally, at the same time, this secondary process of picosecond electron transfer is comparable in time scale to typical electron trapping into defect states in the energy gap.

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

Molecular Dynamics Simulation of the Structures, Dynamics, and Aggregation of Dissolved Organic Matter

Dissolved organic matter (DOM) plays a significant role in the transport and transformation of pollutants in the aquatic environment. However, the experimental characterization of DOM has been limited mainly to bulk properties, and the molecular-level interactions among various components of DOM remain to be fully characterized. In this study, we use molecular dynamics (MD) simulations to probe the structural properties of model DOM systems at atomic detail. The 200 ns simulations, validated by available experimental data, reveal processes and mechanisms by which chemical species (cations, peptides, lipids, lignin, carbohydrates, and some low-molecular-weight aliphatic and aromatic compounds) aggregate to form complex DOM. The DOM aggregates are dynamic, consisting of a hydrophobic core and amphiphilic exterior. The lipid tails and other hydrophobic fragments form the core, with hydrophilic and amphiphilic groups exposed to water, making DOM accessible to both polar and nonpolar species. Thus, the lipid component acts as a nucleator, whereas cations (especially Ca 2+ ) connect the molecular fragments on the surface by coordinating with the O-containing functional groups of DOM. The structural details revealed here provide new insights including surface accessible atoms, overall assemblage, and interactions among the molecules of DOM for understanding the kinetics and mechanisms through which DOM interacts with metal and other contaminants.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Structure, solvation, and dynamics of Mg{sup 2+}, Ca{sup 2+}, Sr{sup 2+}, and Ba{sup 2+} complexes with 3-hydroxyflavone and perchlorate anion in acetonitrile medium: A molecular dynamics simulation study

Molecular dynamics simulations of complexes of Mg{sup 2+}, Ca{sup 2+}, Sr{sup 2+}, and Ba{sup 2+} with 3-hydroxyflavone (flavonol, 3HF) and ClO {sub 4}{sup −} in acetonitrile were performed. The united atoms force field model was proposed for the 3HF molecule using the results of DFT quantum chemical calculations. 3HF was interpreted as a rigid molecule with two internal degrees of freedom, i.e., rotation of the phenyl ring and of the OH group with respect to the chromone moiety. The interatomic radial distribution functions showed that interaction of the cations with flavonol occurs via the carbonyl group of 3HF and it is accompanied with substitution of one of the acetonitrile molecules in the cations’ first solvation shells. Formation of the cation–3HF complexes does not have significant impact on the rotation of the phenyl ring with respect to the chromone moiety. However, the orientation of the flavonol's OH-group is more sensitive to the interaction with doubly charged cations. When complex with Mg{sup 2+} is formed, the OH-group turns out of the plane of the chromone moiety that leads to rupture of intramolecular H-bond in the ligand molecule. Complexation of Ca{sup 2+}, Sr{sup 2+}, and BaClO {sub 4}{sup +} with 3HF produces two structures with different OH-positions, as in the free flavonol with the intramolecular H-bond and as in the complex with Mg{sup 2+} with disrupted H-bonding. It was shown that additional stabilization of the [MgClO{sub 4}(3HF)]{sup +} and [BaClO{sub 4}(3HF)]{sup +} complexes is determined by strong affinity of perchlorate anion to interact with flavonol via intracomplex hydrogen bond between an oxygen atom of the anion and the hydrogen atom of the 3-hydroxyl group. Noticeable difference in the values of the self-diffusion coefficients for Kt{sup 2+} from one side and ClO {sub 4}{sup −}, 3HF, and AN in the cations’ coordination shell from another side implies quite weak interaction between cation, anion, and ligands in the investigated complexes.

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

The alloying element dependence of the local lattice deformation and the elastic properties of Ni{sub 3}Al: A molecular dynamics simulation

Molecular dynamics (MD) together with the modified analytical embedded atom method (MAEAM) is employed to study the alloying elements (Re, Ru, Co, and Ta) dependence of the elastic properties of L1{sub 2}-Ni{sub 3}Al. The investigations indicate that the calculated elastic properties of Ni{sub 3}Al are in reasonable agreement with the previous results. The substituting formation energies of the alloying elements in Ni{sub 3}Al are calculated to determine the site preference. It is found out that Re, Ru, and Ta atoms prefer to occupy the Al sites, and the Co atom prefers to occupy the Ni site. Based on Re, Ru, and Ta substituting the 1st, 2nd, 3rd, and 4th nearest-neighbor atoms, we ascertain that the substituting manners of these alloying elements have a decisive effect on the bulk modulus and the local crystal lattice of Ni{sub 3}Al. Moreover, for Re, the bonding interaction plays a predominant role in the improvement in the bulk modulus of Ni{sub 3}Al, whereas the size effect of Ru and Ta on the improvement in the bulk modulus is more obvious. Finally, the potential reasons of which the alloying elements enhance the bulk modulus are discussed in details.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Mechanisms for hyperthermal polyatomic hydrocarbon modification of PMMA surfaces from molecular dynamics simulations

Classical molecular dynamics simulations are performed to determine the mechanisms by which hyperthermal hydrocarbon polyatomics, which are present in low-energy plasmas, chemically modify polymer surfaces. In particular, C 2 H, CH 3 , and C 3 H 5 are deposited on an amorphous poly (methyl methacrylate) (PMMA) substrate with kinetic energies of 4, 10, 25, and 50 eV and compared to the deposition of H at the same energies. The short-range forces on the atoms are determined using the second generation reactive empirical many-body potential, while the long-range forces are determined using a Lennard-Jones potential. The simulations predict that at all these incident energies, the chemical modification of the PMMA is limited to within a nanometer of the surface. Atoms, fragments, and incident polyatomics are further predicted to chemically attach to specific sites on the PMMA monomers at low energies and to attach to a wider range of sites at higher energies. However, no appreciable cross-linking between polymer chains is predicted to occur. Variation in the penetration depth of the deposited polyatomics or H is correlated to differences in their size and bond saturation. The greatest extent of chemical modification of the PMMA surface slab is achieved for C 2 H deposition with 50 eV of kinetic energy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Predicting 1 H NMR relaxation in Gd 3+ -aqua using molecular dynamics simulations

Atomistic molecular dynamics simulations are used to predict 1 H NMR T 1 relaxation of water from paramagnetic Gd 3+ ions in solution at 25 °C. Simulations of the T 1 relaxivity dispersion function r 1 computed from the Gd 3+ – 1 H dipole–dipole autocorrelation function agree within ≃8% of measurements in the range f 0 ≃ 5 ↔ 500 MHz, without any adjustable parameters in the interpretation of the simulations, and without any relaxation models. Here the simulation results are discussed in the context of the Solomon-Bloembergen-Morgan inner-sphere relaxation model, and the Hwang-Freed outer-sphere relaxation model. Below f 0 ≲ 5 MHz, the simulation overestimates r 1 compared to measurements, which is used to estimate the zero-field electron-spin relaxation time. The simulations show potential for predicting r 1 at high frequencies in chelated Gd 3+ contrast-agents used for clinical MRI.

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

Ab initio molecular dynamics simulations reveal localization and time evolution dynamics of an excess electron in heterogeneous CO{sub 2}–H{sub 2}O systems

In view of the important implications of excess electrons (EEs) interacting with CO{sub 2}–H{sub 2}O clusters in many fields, using ab initio molecular dynamics simulation technique, we reveal the structures and dynamics of an EE associated with its localization and subsequent time evolution in heterogeneous CO{sub 2}–H{sub 2}O mixed media. Our results indicate that although hydration can increase the electron-binding ability of a CO{sub 2} molecule, it only plays an assisting role. Instead, it is the bending vibrations that play the major role in localizing the EE. Due to enhanced attraction of CO{sub 2}, an EE can stably reside in the empty, low-lying π{sup *} orbital of a CO{sub 2} molecule via a localization process arising from its initial binding state. The localization is completed within a few tens of femtoseconds. After EE trapping, the ∠OCO angle of the core CO{sub 2}{sup −} oscillates in the range of 127°∼142°, with an oscillation period of about 48 fs. The corresponding vertical detachment energy of the EE is about 4.0 eV, which indicates extreme stability of such a CO{sub 2}-bound solvated EE in [CO{sub 2}(H{sub 2}O){sub n}]{sup −} systems. Interestingly, hydration occurs not only on the O atoms of the core CO{sub 2}{sup −} through formation of O⋯H–O H–bond(s), but also on the C atom, through formation of a C⋯H–O H–bond. In the latter binding mode, the EE cloud exhibits considerable penetration to the solvent water molecules, and its IR characteristic peak is relatively red-shifted compared with the former. Hydration on the C site can increase the EE distribution at the C atom and thus reduce the C⋯H distance in the C⋯H–O H–bonds, and vice versa. The number of water molecules associated with the CO{sub 2}{sup −} anion in the first hydration shell is about 4∼7. No dimer-core (C{sub 2}O{sub 4}{sup −}) and core-switching were observed in the double CO{sub 2} aqueous media. This work provides molecular dynamics insights into the localization and time evolution dynamics of an EE in heterogeneous CO{sub 2}–H{sub 2}O media.

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

Photodissociation dynamics of CH{sub 3}C(O)SH in argon matrix: A QM/MM nonadiabatic dynamics simulation

In this work, we have first employed the combined quantum mechanics/molecular mechanics (QM/MM) method to study the photodissociation mechanism of thioacetic acid CH{sub 3}C(O)SH in the S{sub 1}, T{sub 1}, and S{sub 0} states in argon matrix. CH{sub 3}C(O)SH is treated quantum mechanically using the complete active space self-consistent field and complete active space second-order perturbation theory methods; argon matrix is described classically using Lennard-Jones potentials. We find that the C-S bond fission is predominant due to its small barriers of ca. 3.0 and 1.0 kcal/mol in the S{sub 1} and T{sub 1} states. It completely suppresses the nearby C—C bond fission. After the bond fission, the S{sub 1} radical pair of CH{sub 3}CO and SH can decay to the S{sub 0} and T{sub 1} states via internal conversion and intersystem crossing, respectively. In the S{sub 0} state, the radical pair can either recombine to form CH{sub 3}C(O)SH or proceed to form molecular products of CH{sub 2}CO and H{sub 2}S. We have further employed our recently developed QM/MM generalized trajectory-based surface-hopping method to simulate the photodissociation dynamics of CH{sub 3}C(O)SH. In 1 ps dynamics simulation, 56% trajectories stay at the Franck-Condon region; the S{sub 1} C—S bond fission takes place in the remaining 44% trajectories. Among all nonadiabatic transitions, the S{sub 1} → S{sub 0} internal conversion is major (55%) but the S{sub 1} → T{sub 1} intersystem crossing is still comparable and cannot be ignored, which accounts for 28%. Finally, we have found a radical channel generating the molecular products of CH{sub 2}CO and H{sub 2}S, which is complementary to the concerted molecular channel. The present work sets the stage for simulating photodissociation dynamics of similar thio-carbonyl systems in matrix.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Tether Dynamics Simulation Workshop summary

The Tether Dynamics Simulation Workshop focused on the efforts of the Tether Applications Simulation Working Group (TASWG) to catalog various dynamics simulations, document environmental models, and provide a set of verified simulation results. The workshop consisted of reports on these activities and the presentation of papers on current dynamics topics of general interest. This paper documents the work performed at the Workshop and the results and recommendations for further work.

Rupp, Charles C.↗
DOE OSTI2015

Molecular dynamics simulations of large macromolecular complexes

Connecting dynamics to structural data from diverse experimental sources, molecular dynamics simulations permit the exploration of biological phenomena in unparalleled detail. Advances in simulations are moving the atomic resolution descriptions of biological systems into the million-to-billion atom regime, in which numerous cell functions reside. In this opinion, we review the progress, driven by large-scale molecular dynamics simulations, in the study of viruses, ribosomes, bioenergetic systems, and other diverse applications. Finally, these examples highlight the utility of molecular dynamics simulations in the critical task of relating atomic detail to the function of supramolecular complexes, a task that cannot be achieved by smaller-scale simulations or existing experimental approaches alone.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Study of ice cluster impacts on amorphous silica using the ReaxFF reactive force field molecular dynamics simulation method

We study the dynamics of the collisions between amorphous silica structures and amorphous and crystal ice clusters with impact velocities of 1 km/s, 4 km/s, and 7 km/s using the ReaxFF reactive molecular dynamics simulation method. The initial ice clusters consist of 150 water molecules for the amorphous ice cluster and 128 water molecules for the crystal ice cluster. The ice clusters are collided on the surface of amorphous fully oxidized and suboxide silica. These simulations show that at 1 km/s impact velocities, all the ice clusters accumulate on the surface and at 4 km/s and 7 km/s impact velocities, some of the ice cluster molecules bounce back from the surface. At 4 km/s and 7 km/s impact velocities, few of the water molecules dissociations are observed. The effect of the second ice cluster impacts on the surfaces which are fully covered with ice, on the mass loss/accumulation is studied. These studies show that at 1 km/s impacts, the entire ice cluster accumulates on the surface at both first and second ice impacts. At higher impact velocities, some ice molecules which after the first ice impacts have been attached to the surface will separate from the surface after the second ice impacts at 7 km/s impact velocity. For the 4 km/s ice cluster impact, ice accumulation is observed for the crystal ice cluster impacts and ice separation is observed for the amorphous ice impacts. Observing the temperatures of the ice clusters during the collisions indicates that the possibility of electron excitement at impact velocities less than 10 km/s is minimal and ReaxFF reactive molecular dynamics simulation can predict the chemistry of these hypervelocity impacts. However, at impact velocities close to 10 km/s the average temperature of the impacting ice clusters increase to about 2000 K, with individual molecules occasionally reaching temperatures of over 8000 K and thus it will be prudent to consider the concept of electron excitation at these higher impact velocities, which goes beyond the current ReaxFF ability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Study on the thermal resistance in secondary particles chain of silica aerogel by molecular dynamics simulation

In this article, molecular dynamics simulation was performed to study the heat transport in secondary particles chain of silica aerogel. The two adjacent particles as the basic heat transport unit were modelled to characterize the heat transfer through the calculation of thermal resistance and vibrational density of states (VDOS). The total thermal resistance of two contact particles was predicted by non-equilibrium molecular dynamics simulations (NEMD). The defects were formed by deleting atoms in the system randomly first and performing heating and quenching process afterwards to achieve the DLCA (diffusive limited cluster-cluster aggregation) process. This kind of treatment showed a very reasonable prediction of thermal conductivity for the silica aerogels compared with the experimental values. The heat transport was great suppressed as the contact length increased or defect concentration increased. The constrain effect of heat transport was much significant when contact length fraction was in the small range (<0.5) or the defect concentration is in the high range (>0.5). Also, as the contact length increased, the role of joint thermal resistance played in the constraint of heat transport was increasing. However, the defect concentration did not affect the share of joint thermal resistance as the contact length did. VDOS of the system was calculated by numerical method to characterize the heat transport from atomic vibration view. The smaller contact length and greater defect concentration primarily affected the longitudinal acoustic modes, which ultimately influenced the heat transport between the adjacent particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2010

Two-Speed Gearbox Dynamic Simulation Predictions and Test Validation

Dynamic simulations and experimental validation tests were performed on a two-stage, two-speed gearbox as part of the drive system research activities of the NASA Fundamental Aeronautics Subsonics Rotary Wing Project. The gearbox was driven by two electromagnetic motors and had two electromagnetic, multi-disk clutches to control output speed. A dynamic model of the system was created which included a direct current electric motor with proportional-integral-derivative (PID) speed control, a two-speed gearbox with dual electromagnetically actuated clutches, and an eddy current dynamometer. A six degree-of-freedom model of the gearbox accounted for the system torsional dynamics and included gear, clutch, shaft, and load inertias as well as shaft flexibilities and a dry clutch stick-slip friction model. Experimental validation tests were performed on the gearbox in the NASA Glenn gear noise test facility. Gearbox output speed and torque as well as drive motor speed and current were compared to those from the analytical predictions. The experiments correlate very well with the predictions, thus validating the dynamic simulation methodologies.

Lewicki, David G.↗
DOE OSTI2015

Density relaxation and particle motion characteristics in a non-ionic deep eutectic solvent (acetamide + urea): Time-resolved fluorescence measurements and all-atom molecular dynamics simulations

Temperature dependent relaxation dynamics, particle motion characteristics, and heterogeneity aspects of deep eutectic solvents (DESs) made of acetamide (CH{sub 3}CONH{sub 2}) and urea (NH{sub 2}CONH{sub 2}) have been investigated by employing time-resolved fluorescence measurements and all-atom molecular dynamics simulations. Three different compositions (f) for the mixture [fCH{sub 3}CONH{sub 2} + (1 − f)NH{sub 2}CONH{sub 2}] have been studied in a temperature range of 328-353 K which is ∼120-145 K above the measured glass transition temperatures (∼207 K) of these DESs but much lower than the individual melting temperature of either of the constituents. Steady state fluorescence emission measurements using probe solutes with sharply different lifetimes do not indicate any dependence on excitation wavelength in these metastable molten systems. Time-resolved fluorescence anisotropy measurements reveal near-hydrodynamic coupling between medium viscosity and rotation of a dissolved dipolar solute. Stokes shift dynamics have been found to be too fast to be detected by the time-resolution (∼70 ps) employed, suggesting extremely rapid medium polarization relaxation. All-atom simulations reveal Gaussian distribution for particle displacements and van Hove correlations, and significant overlap between non-Gaussian (α{sub 2}) and new non-Gaussian (γ) heterogeneity parameters. In addition, no stretched exponential relaxations have been detected in the simulated wavenumber dependent acetamide dynamic structure factors. All these results are in sharp contrast to earlier observations for ionic deep eutectics with acetamide [Guchhait et al., J. Chem. Phys. 140, 104514 (2014)] and suggest a fundamental difference in interaction and dynamics between ionic and non-ionic deep eutectic solvent systems.

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

Constant-depth circuits for dynamic simulations of materials on quantum computers

Abstract Dynamic simulation of materials is a promising application for near-term quantum computers. Current algorithms for Hamiltonian simulation, however, produce circuits that grow in depth with increasing simulation time, limiting feasible simulations to short-time dynamics. Here, we present a method for generating circuits that are constant in depth with increasing simulation time for a specific subset of one-dimensional (1D) materials Hamiltonians, thereby enabling simulations out to arbitrarily long times. Furthermore, by removing the effective limit on the number of feasibly simulatable time-steps, the constant-depth circuits enable Trotter error to be made negligibly small by allowing simulations to be broken into arbitrarily many time-steps. For an N -spin system, the constant-depth circuit contains only $\mathcal {O}(N^{2})$ O ( N 2 ) CNOT gates. Such compact circuits enable us to successfully execute long-time dynamic simulation of ubiquitous models, such as the transverse field Ising and XY models, on current quantum hardware for systems of up to 5 qubits without the need for complex error mitigation techniques. Aside from enabling long-time dynamic simulations with minimal Trotter error for a specific subset of 1D Hamiltonians, our constant-depth circuits can advance materials simulations on quantum computers more broadly in a number of indirect ways.

Materials simulation↗
NASA NTRS1990

Reduced order component models for flexible multibody dynamics simulations

Many flexible multibody dynamics simulation codes require some form of component description that properly characterizes the dynamic behavior of the system. A model reduction procedure for producing low order component models for flexible multibody simulation is described. Referred to as projection and assembly, the method is a Rayleigh-Ritz approach that uses partitions of the system modal matrix as component Ritz transformation matrices. It is shown that the projection and assembly method yields a reduced system model that preserves a specified set of the full order system modes. Unlike classical component mode synthesis methods, the exactness of the method described is obtained at the expense of having to compute the full order system modes. The paper provides a comprehensive description of the method, a proof of exactness, and numerical results demonstrating the method's effectiveness.

Tsuha, Walter S.↗
NASA NTRSDate not supplied

Computational Fluid Dynamics Simulations of the Transonic Dynamics Tunnel Airstream Oscillator System

This paper presents a Computational Fluid Dynamics (CFD) model of the flow in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT) with the Airstream Oscillator System (AOS) in operation. The TDT is a continuous-flow, closed circuit wind tunnel with a 16- by 16-foot slotted test section with cropped corners. The tunnel was originally built as the 19-ft Pressure Tunnel in 1938, but it was converted to the current transonic tunnel in the 1950s, with capabilities to use either air or heavy gas (R-134a) as the test medium. The TDT was also fitted with an AOS which can be used to create a gust field in the test section. To date, no computational analyses of the TDT involving the AOS have been performed. In this study, experimental data acquired of an oscillating airstream in the tunnel will be used to calibrate the computational analyses. A significant motivation of this work is to attempt to compare computational gust velocities to those recorded in the TDT literature. In addition to a validation of this model with experimental data, it may also be possible to supplement the experimental data with computational data. The experimental data is rather sparse and at only a few Mach numbers. Computational results may be able to expand the AOS data set. Another motivation for this work is that the Integrated Adaptive Wing Technology Maturation (IAWTM) semi-span model of the high-aspect-ratio CRM with 10 trailing edge control surfaces is currently being fabricated and will be delivered to the NASA Langley Transonic Dynamics Tunnel (TDT) for testing in late 2020. Among tests to be conducted, gust load alleviation (GLA) will be demonstrated at transonic conditions using the AOS.

Computational Fluid Dynamics↗