Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular Dynamics”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Simulation studies for surfaces and materials strength

Investigations were carried out in the following three areas. First, a simulation study was conducted to investigate the early stages of a nucleation and growth mechanism taking place on diamond surfaces. The energetically most favorable binding sites were calculated along with the binding energies for up to three carbon atoms deposited on the planes of a diamond crystal. Secondly, calculations were carried out to analyze energies and structural properties for restructured patterns of the surface of a diamond. Two different model functions, developed recently for carbon, were employed in the calculations. Lastly, a crack propagation process for the graphitic basal plane, was investigated using a molecular dynamics technique. Interactions among carbon atoms in the system were calculated using the Tersoff potential energy function.

Halicioglu, Timur↗
NASA NTRS1992

Simulation studies for surfaces and materials strength

Investigations were carried out in two major areas during the last reporting period. Energy- and structure-related properties of small gold clusters deposited on the GaAs(110) surface were investigated using a molecular dynamics procedure. Additionally, a comparative study of the many-body potentials of silicon systems was performed.

Halicioglu, Timur↗
NASA NTRS1992

Neural network error correction for solving coupled ordinary differential equations

A neural network is presented to learn errors generated by a numerical algorithm for solving coupled nonlinear differential equations. The method is based on using a neural network to correctly learn the error generated by, for example, Runge-Kutta on a model molecular dynamics (MD) problem. The neural network programs used in this study were developed by NASA. Comparisons are made for training the neural network using backpropagation and a new method which was found to converge with fewer iterations. The neural net programs, the MD model and the calculations are discussed.

Shelton, R. O.↗
NASA NTRS1991

Surface disordering of Pb(110)

A molecular dynamics simulation of Pb(110), using embedded atom method energy functional and two-body potential derived for the purpose, shows that near 400 K the three surface layers begin to disorder well before the bulk (interior) layers. Transfer of atoms from the first subsurface layer to the surface layer occurs. Disordering propagates beyond the top three layers above 550 K, accompanied by formation of an adlayer and vacancies in the top three layers. Behavior of the two-dimensional layer structure factors indicates that disordering is anisotropic. Simulation results are consistent with experimental observations of surface roughening near 400 K and more extensive surface ordering above 525 K. Results are consistent with simulations for Ni and Al.

Tibbits, P.↗
NASA NTRS1993

Momentum loss in proton-nucleus and nucleus-nucleus collisions

An optical model description, based on multiple scattering theory, of longitudinal momentum loss in proton-nucleus and nucleus-nucleus collisions is presented. The crucial role of the imaginary component of the nucleon-nucleon transition matrix in accounting for longitudinal momentum transfer is demonstrated. Results obtained with this model are compared with Intranuclear Cascade (INC) calculations, as well as with predictions from Vlasov-Uehling-Uhlenbeck (VUU) and quantum molecular dynamics (QMD) simulations. Comparisons are also made with experimental data where available. These indicate that the present model is adequate to account for longitudinal momentum transfer in both proton-nucleus and nucleus-nucleus collisions over a wide range of energies.

Khan, Ferdous↗
NASA NTRS1993

Thermophysical properties of simple liquid metals: A brief review of theory

In this paper, we review the current theory of the thermophysical properties of simple liquid metals. The emphasis is on thermodynamic properties, but we also briefly discuss the nonequilibrium properties of liquid metals. We begin by defining a 'simple liquid metal' as one in which the valence electrons interact only weakly with the ionic cores, so that the interaction can be treated by perturbation theory. We then write down the equilibrium Hamiltonian of a liquid metal as a sum of five terms: the bare ion-ion interaction, the electron-electron interaction, the bare electron-ion interaction, and the kinetic energies of electrons and ions. Since the electron-ion interaction can be treated by perturbation, the electronic part contributes in two ways to the Helmholtz free energy: it gives a density-dependent term which is independent of the arrangement of ions, and it acts to screen the ion-ion interaction, giving rise to effective ion-ion pair potentials which are density-dependent, in general. After sketching the form of a typical pair potential, we briefly enumerate some methods for calculating the ionic distribution function and hence the Helmholtz free energy of the liquid: monte Carlo simulations, molecular dynamics simulations, and thermodynamic perturbation theory. The final result is a general expression for the Helmholtz free energy of the liquid metal. It can be used to calculate a wide range of thermodynamic properties of simple metal liquids, which we enumerate. They include not only a range of thermodynamic coefficients of both metals and alloys, but also many aspects of the phase diagram, including freezing curves of pure elements and phase diagrams of liquid alloys (including liquidus and solidus curves). We briefly mention some key discoveries resulting from previous applications of this method, and point out that the same methods work for other materials not normally considered to be liquid metals (such as colloidal suspensions, in which the suspended microspheres behave like ions screened by the salt solution in which they are suspended). We conclude with a brief discussion of some non-equilibrium (i.e., transport) properties which can be treated by an extension of these methods. These include electrical resistivity, thermal conductivity, viscosity, atomic self-diffusion coefficients, concentration diffusion coefficients in alloys, surface tension and thermal emissivity. Finally, we briefly mention two methods by which the theory might be extended to non-simple liquid metals: these are empirical techniques (i.e., empirical two- and three-body potentials), and numerical many-body approaches. Both may be potentially applicable to extremely complex systems, such as nonstoichiometric liquid semiconductor alloys.

Stroud, David↗
NASA NTRS1994

Powder agglomeration in a microgravity environment

This is the final report for NASA Grant NAG3-755 entitled 'Powder Agglomeration in a Microgravity Environment.' The research program included both two types of numerical models and two types of experiments. The numerical modeling included the use of Monte Carlo type simulations of agglomerate growth including hydrodynamic screening and molecular dynamics type simulations of the rearrangement of particles within an agglomerate under a gravitational field. Experiments included direct observation of the agglomeration of submicron alumina and indirect observation, using small angle light scattering, of the agglomeration of colloidal silica and aluminum monohydroxide. In the former class of experiments, the powders were constrained to move on a two-dimensional surface oriented to minimize the effect of gravity. In the latter, some experiments involved mixture of suspensions containing particles of opposite charge which resulted in agglomeration on a very short time scale relative to settling under gravity.

Cawley, James D.↗
NASA NTRS1994

Wetting and spreading at the molecular scale

We have studied the microscopic aspects of the spreading of liquid drops on a solid surface by molecular dynamics simulations of coexisting three-phase Lennard-Jones systems of liquid, vapor and solid. We consider both spherically symmetric atoms and chain-like molecules, and a range of interaction strengths. As the attraction between liquid and solid increases we observed a smooth transition in spreading regimes, from partial to complete to terraced wetting. In the terraced case, where distinct monomolecular layers spread with different velocities, the layers are ordered but not solid, with qualitative behavior resembling recent experimental findings, but with interesting differences in the spreading rate.

Koplik, Joel↗
NASA NTRS1994

Storage and recycling of water and carbon dioxide in the earth

The stabilities and properties of water- and carbon-bearing phases in the earth have been determined from phase equilibrium measurements, combined with new data on the equations of state of water, carbon dioxide, carbonates and hydrates. The data have then been used to predict the fate of calcite and hydrous phases in subducting oceanic lithosphere. From the compositions of MORB's one can estimate concentrations of water and carbon of around 200 ppm and 80 ppm respectively in the upper mantle. Lower mantle estimates are very uncertain, but 1900 ppm water and 2000 ppm C are plausible concentrations. Measurements of the density of supercritical water to 3 GPa demonstrate that this phase is less compressible than anticipated from the equations of state of Haar et al. or Saul and Wagner and is closer to predictions based on molecular dynamics simulations. Conversely, fugacity measurements on carbon dioxide to 7 GPa show that this fluid is more compressible than predicted from the MRK equation of state. The results imply that hydrates are relatively more stable and carbonates less stable at pressures greater than 5 GPa than would be predicted from simple extrapolation of the low pressure data. Nevertheless, carbonates remain extremely refractory phases within both the upper and lower mantle.

Wood, Bernard J.↗
NASA NTRS1995

Sensitivity analysis in multipole-accelerated panel methods for potential flow

In the design of an airframe, the effect of changing the geometry on resulting computations is necessary for design optimization. The geometry is defined in terms of a series of design variables, including design variables to define the wing planform, tail, canard, pylon, and nacelle. Design optimization in this research is based on how these design variable affect the potential flow. The potential flow is computed as a function of the geometry and location of a series of panels describing the airframe, which are in turn a function of the design variables. Multipole accelerated panel methods improve the computational complexity of the problem and thus are an attractive approach. To utilize the methods in design optimization, it was necessary to define the appropriate sensitivity derivatives. The overhead incurred from finding the sensitivity derivatives in conjunction with the original computation should be small. This research developed the background for multipole-accelerated panel methods and the framework for finding sensitivity derivatives in the methods. Potential flow panel codes are commonly used for powered-lift aerodynamic predictions for three dimensional geometries. Given an airframe which has been discretized into a series of panels to define the airframe geometry, potential is computed as a function of the influence of all panels on all other panels. This is a computationally intensive problem for which efficient solutions are desired to improve the computational time and to allow greater resolution by use of more panels. One such solution is the use of hierarchical multipole methods which entail approximations of the effects of far-field terms. Hierarchical multipole methods have become prevalent in molecular dynamics and gravitational physics, and have been introduced into the fields of capacitance calculations, computational fluid dynamics, and electromagnetics. The methods utilize multipole expansions to describe the effect of bodies (i.e. particles, astrophysical bodies, panels, etc.) within a sphere on points distant from the sphere, where the influence diminishes as a function of distance. The expansions are exact with infinite series, however, for practical computations, the series are truncated and accuracy is selected based on the number of terms retained in the expansions. A hierarchical tree structure groups bodies together based on proximity to allow definition of multipole expansions for each group. The multipole expansions are then used to compute the effect of the bodies in a group on distant bodies.

Leathrum, James F., Jr.↗
NASA NTRS1995

Fast Multipole Methods for Three-Dimensional N-body Problems

We are developing computational tools for the simulations of three-dimensional flows past bodies undergoing arbitrary motions. High resolution viscous vortex methods have been developed that allow for extended simulations of two-dimensional configurations such as vortex generators. Our objective is to extend this methodology to three dimensions and develop a robust computational scheme for the simulation of such flows. A fundamental issue in the use of vortex methods is the ability of employing efficiently large numbers of computational elements to resolve the large range of scales that exist in complex flows. The traditional cost of the method scales as Omicron (N(sup 2)) as the N computational elements/particles induce velocities at each other, making the method unacceptable for simulations involving more than a few tens of thousands of particles. In the last decade fast methods have been developed that have operation counts of Omicron (N log N) or Omicron (N) (referred to as BH and GR respectively) depending on the details of the algorithm. These methods are based on the observation that the effect of a cluster of particles at a certain distance may be approximated by a finite series expansion. In order to exploit this observation we need to decompose the element population spatially into clusters of particles and build a hierarchy of clusters (a tree data structure) - smaller neighboring clusters combine to form a cluster of the next size up in the hierarchy and so on. This hierarchy of clusters allows one to determine efficiently when the approximation is valid. This algorithm is an N-body solver that appears in many fields of engineering and science. Some examples of its diverse use are in astrophysics, molecular dynamics, micro-magnetics, boundary element simulations of electromagnetic problems, and computer animation. More recently these N-body solvers have been implemented and applied in simulations involving vortex methods. Koumoutsakos and Leonard (1995) implemented the GR scheme in two dimensions for vector computer architectures allowing for simulations of bluff body flows using millions of particles. Winckelmans presented three-dimensional, viscous simulations of interacting vortex rings, using vortons and an implementation of a BH scheme for parallel computer architectures. Bhatt presented a vortex filament method to perform inviscid vortex ring interactions, with an alternative implementation of a BH scheme for a Connection Machine parallel computer architecture.

Koumoutsakos, P.↗
NASA NTRS1996

Electronic Characterization of Defects in Narrow Gap Semiconductors-Comparison of Electronic Energy Levels and Formation Energies in Mercury Cadmium Telluride, Mercury Zinc Telluride, and Mercury Zinc Selenide

We have used a Green's function technique to calculate the energy levels and formation energy of deep defects in the narrow gap semiconductors mercury cadmium telluride (MCT), mercury zinc telluride (MZT) and mercury zinc selenide (MZS). The formation energy is calculated from the difference between the total energy with an impurity cluster and the total energy for the perfect crystal. Substitutional (including antisite), interstitial (self and foreign), and vacancy deep defects are considered. Relaxation effects are calculated (with molecular dynamics). By use of a pseudopotential, we generalize the ideal vacancy model so as to be able to consider relaxation for vacancies. Different charge states are considered and the charged state energy shift (as computed by a modified Haldane-Anderson model) can be twice that due to relaxation. Different charged states for vacancies were not calculated to have much effect on the formation energy. For all cases we find deep defects in the energy gap only for cation site s-like orbitals or anion site p-like orbitals, and for the substitutional case only the latter are appreciably effected by relaxation. For most cases for MCT, MZT, MZS, we consider x (the concentration of Cd or Zn) in the range appropriate for a band gap of 0.1 eV. For defect energy levels, the absolute accuracy of our results is limited, but the precision is good, and hence chemical trends are accurately predicted. For the same reason, defect formation energies are more accurately predicted than energy level position. We attempt, in Appendix B, to calculate vacancy formation energies using relatively simple chemical bonding ideas due to Harrison. However, these results are only marginally accurate for estimating vacancy binding energies. Appendix C lists all written reports and publications produced for the grant. We include abstracts and a complete paper that summarizes our work which is not yet available.

Patterson, James D.↗
NASA NTRS1996

The Micromechanics of the Moving Contact Line

The proposed research is divided into three components concerned with molecular structure, molecular orientation, and continuum averages of discrete systems. In the experimental program, we propose exploring how changes in interfacial molecular structure generate contact line motion. Rather than rely on the electrostatic and electrokinetic fields arising from the molecules themselves, we augment their interactions by an imposed field at the solid/liquid interface. By controling the field, we can manipulate the molecular structure at the solid/liquid interface. In response to controlled changes in molecular structure, we observe the resultant contact line motion. In the analytical portion of the proposed research we seek to formulate a system of equations governing fluid motion which accounts for the orientation of fluid molecules. In preliminary work, we have focused on describing how molecular orientation affects the forces generated at the moving contact line. Ideally, as assumed above, the discrete behavior of molecules can be averaged into a continuum theory. In the numerical portion of the proposed research, we inquire whether the contact line region is, in fact, large enough to possess a well-defined average. Additionally, we ask what types of behavior distinguish discrete systems from continuum systems. Might the smallness of the contact line region, in itself, lead to behavior different from that in the bulk? Taken together, our proposed research seeks to identify and accurately account for some of the molecular dynamics of the moving contact line, and attempts to formulate a description from which one can compute the forces at the moving contact line.

Han, Minsub↗
NASA NTRS1996

Wetting Hysteresis at the Molecular Scale

The motion of a fluid-fluid-solid contact line on a rough surface is well known to display hysteresis in the contact angle vs. velocity relationship. In order to understand the phenomenon at a fundamental microscopic level, we have conducted molecular dynamics computer simulations of a Wilhelmy plate experiment in which a solid surface is dipped into a liquid bath, and the force-velocity characteristics are measured. We directly observe a systematic variation of force and contact angle with velocity, which is single-valued for the case of an atomically smooth solid surface. In the microscopically rough case, however, we find (as intuitively expected) an open hysteresis loop. Further characterization of the interface dynamics is in progress.

Jin, Wei↗
NASA NTRS1996

Computational Materials Research

Computational Materials aims to model and predict thermodynamic, mechanical, and transport properties of polymer matrix composites. This workshop, the second coordinated by NASA Langley, reports progress in measurements and modeling at a number of length scales: atomic, molecular, nano, and continuum. Assembled here are presentations on quantum calculations for force field development, molecular mechanics of interfaces, molecular weight effects on mechanical properties, molecular dynamics applied to poling of polymers for electrets, Monte Carlo simulation of aromatic thermoplastics, thermal pressure coefficients of liquids, ultrasonic elastic constants, group additivity predictions, bulk constitutive models, and viscoplasticity characterization.

Hinkley, Jeffrey A.↗
NASA NTRS1997

Carbon Nanotubules: Building Blocks for Nanometer-Scale Engineering

Proximal probe technology has provided researchers with new ways to investigate and manipulate matter on the nanometer scale. We have studied, through molecular dynamics simulations, using a many-body empirical potential, the indentation of a hydrogen-terminated, diamond (111 ) surface, with a proximal probe tip that consists of an open, hydrogen-terminated, (10,10) carbon nanotubule. The simulations showed that upon indenting 1.8 A, the tubule deforms but returns to its original shape upon retraction. The Young's modulus of the tubule was determined using the predicted Euler buckling force and was found to be comparable to measured and calculated values. In a second series of simulations, an open (10, 10) nanotubule was heated to 4500 K and allowed to close. We find that at this temperature the resulting cap contains numerous imperfections, including some not mentioned previously in the literature.

Sinnott, Susan B.↗
NASA NTRS1995

High-Density Amorphous Ice, the Frost on Interstellar Grains

Most water ice in the universe is in a form which does not occur naturally on Earth and of which only minimal amounts have been made in the laboratory. We have encountered this 'high-density amorphous ice' in electron diffraction experiments of low-temperature (T less than 30 K) vapor-deposited water and have subsequently modeled its structure using molecular dynamics simulations. The characteristic feature of high-density amorphous ice is the presence of 'interstitial' oxygen pair distances between 3 and 4 A. However, we find that the structure is best described as a collapsed lattice of the more familiar low-density amorphous form. These distortions are frozen in at temperatures below 38 K because, we propose, it requires the breaking of one hydrogen bond, on average, per molecule to relieve the strain and to restructure the lattice to that of low-density amorphous ice. Several features of astrophysical ice analogs studied in laboratory experiments are readily explained by the structural transition from high-density amorphous ice into low-density amorphous ice. Changes in the shape of the 3.07 gm water band, trapping efficiency of CO, CO loss, changes in the CO band structure, and the recombination of radicals induced by low-temperature UV photolysis all covary with structural changes that occur in the ice during this amorphous to amorphous transition. While the 3.07 micrometers ice band in various astronomical environments can be modeled with spectra of simple mixtures of amorphous and crystalline forms, the contribution of the high-density amorphous form nearly always dominates.

Jenniskens, P.↗
NASA NTRS1997

Roto-Translational Collision-Induced Absorption of CO2 for the Atmosphere of Venus at Frequencies from 0 to 250 cm(exp -1), at Temperatures from 200 to 800 K

The collision-induced absorption (CIA)of gaseous CO2 is the primary source of far-infrared opacity of the atmosphere of Venus. At the temperatures and densities of the venusian atmosphere, the absorption is due mainly to binary collisions of CO2 molecules. Using a realistic anisotropic intermolecular potential and assuming the absorbing dipole to be due to the electrostatic induction and a quantum overlap, a series of molecular dynamics simulations were performed for the temperature range 200 to 800 K, and the roto-translational (RT) collision-induced absorption spectra at frequencies from 0 to 250 cm(exp -1) were derived. The absorption coefficient in the submillimeter region, used in constituency retrieval studies, decreases more than 10 times in the temperature range 200 to 800 K. On the other hand, the absorption coefficient at 800 K and at the frequency range above 150 cm(exp -1) was found to be almost 10 times higher than at 200 K. Earlier works relied on experimental RT CIA data at a fixed temperature of 300 K. The new, temperature-dependent absorption bands may, when included in the analysis of the atmospheric radiative transfer of the planet, help explain the observed high far-infrared opacity of the lower layers of the atmosphere. To make the results of the simulations readily available for atmospheric abundance and radiative transfer analysis, an analytic model of the roto-translational collision-induced absorption spectral profile, applicable from 200 to 800 K, is being proposed here. The FORTRAN computer code of this newly developed model is available from the authors on request.

Gruszka, Marcin↗