Search NASA⌕ Search

SEARCH · Search NASA

Results for “MOLECULAR DYNAMICS METHOD”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
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 NTRS1973

Study of methods of improving the performance of the Langley Research Center Transonic Dynamics Tunnel (TDT)

A study has been made of possible ways to improve the performance of the Langley Research Center's Transonic Dynamics Tunnel (TDT). The major effort was directed toward obtaining increased dynamic pressure in the Mach number range from 0.8 to 1.2, but methods to increase Mach number capability were also considered. Methods studied for increasing dynamic pressure capability were higher total pressure, auxiliary suction, reducing circuit losses, reduced test medium temperature, smaller test section and higher molecular weight test medium. Increased Mach number methods investigated were nozzle block inserts, variable geometry nozzle, changes in test section wall configuration, and auxiliary suction.

Source record↗
NASA NTRS2000

A New Clinical Instrument for The Early Detection of Cataract Using Dynamic Light Scattering and Corneal Topography

A growing cataract can be detected at the molecular level using the technique of dynamic light scattering (DLS). However, the success of this method in clinical use depends upon the precise control of the scattering volume inside a patient's eye and especially during patient's repeat visits. This is important because the scattering volume (cross-over region between the scattered fight and incident light) inside the eye in a high-quality DLS set-up is very small (few microns in dimension). This precise control holds the key for success in the longitudinal studies of cataract and during anti-cataract drug screening. We have circumvented these problems by fabricating a new DLS fiber optic probe with a working distance of 40 mm and by mounting it inside a cone of a corneal analyzer. This analyzer is frequently used in mapping the corneal topography during PRK (photorefractive keratectomy) and LASIK (laser in situ keratomileusis) procedures in shaping of the cornea to correct myopia. This new instrument and some preliminary clinical tests on one of us (RRA) showing the data reproducibility are described.

Ansari, Rafat R.↗
NASA NTRS2013

Conductance of Ion Channels - Theory vs. Experiment

Transmembrane ion channels mediate a number of essential physiological processes in a cell ranging from regulating osmotic pressure to transmission of neural signals. Kinetics and selectivity of ion transport is of critical importance to a cell and, not surprisingly, it is a subject of numerous experimental and theoretical studies. In this presentation we will analyze in detail computer simulations of two simple channels from fungi - antiamoebin and trichotoxin. Each of these channels is made of an alpha-helical bundle of small, nongenomically synthesized peptides containing a number of rare amino acids and exhibits strong antimicrobial activity. We will focus on calculating ionic conductance defined as the ratio of ionic current through the channel to applied voltage. From molecular dynamics simulations, conductance can be calculated in at least two ways, each involving different approximations. Specifically, the current, given as the number of charges transferred through the channel per unit of time, can be obtained from the number of events in which ions cross the channel during the simulation. This method works well for large currents (high conductance values and/or applied voltages). If the number of crossing events is small, reliable estimates of current are difficult to achieve. Alternatively, conductance can be estimated assuming that ion transport can be well approximated as diffusion in the external potential given by the free energy profile. Then, the current can be calculated by solving the one-dimensional diffusion equation in this external potential and applied voltage (the generalized Nernst-Planck equation). To do so three ingredients are needed: the free energy profile, the position-dependent diffusion coefficient and the diffusive flux of ions into the channel. All these quantities can be obtained from molecular dynamics simulations. An important advantage of this method is that it can be used equally well to estimating large and small currents. In addition, once the free energy profile becomes available the full current-voltage dependence can be readily obtained. For both channels we carried out calculations using both approaches. We also tested the main assumptions underlying the diffusive model, such as uncorrelated nature of individual crossing events and Fickian diffusion. The accuracy and consistency of different methods will be discussed. Finally we will discuss how comparisons between calculated and measured ionic conductance and selectivity of transport can be used for determining structural models of the channels.

Pohorille, Andrew↗
NASA NTRSDate not supplied

Building and Breaking Carbon Composites with REACTER

Carbon composites have become indispensable for aerospace and other high-performance applications, and a detailed picture of their morphology and failure mechanisms remains difficult to obtain through experiment. REACTER is a versatile computational modeling tool for atomistic molecular dynamics designed to model chemical reactions at the speed and length scales of classical force fields.1 In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites, carbon-fiber reinforced polymers (CFRP) and carbon nanotube (CNT) composites. Carbon fiber core morphologies were created by the method of Desai et al.,2 but using the advanced reaction constraints framework of REACTER, their proposed multistep procedure was reduced to a single uninterrupted molecular dynamics simulation. The carbon fiber filler was embedded into a polymer matrix by simulated in situ polymerization of several thermosetting resins, including bismaleimide and polyarylacetylene, to obtain the final CFRP model. To generate the second class of carbon composite, CNT networks were grown dynamically using the new ‘create atoms’ feature of REACTER, and similarly infiltrated with resin to obtain CNT composites. The resulting models were compared directly to experiment using simulated high resolution transmission electron microscopy and x-ray diffraction. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third-order DFT-based tight-binding (DFTB3) simulations, via a reaction constraint on the total potential energy of the involved atoms.

Molecular Dynamics↗
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 NTRS1975

Dynamic delta method for trace gas analysis

Method has been developed in which measurements are made only over viscous flow range, eliminating fractionation before the molecular leak and problems due to surface elution.

Wood, G. M.↗
NASA NTRS2003

Diagnostics of Carbon Nanotube Formation in a Laser Produced Plume: An Investigation of the Metal Catalyst by Laser Ablation Atomic Fluorescence Spectroscopy

Carbon nanotubes, elongated molecular tubes with diameters of nanometers and lengths in microns, hold great promise for material science. Hopes for super strong light-weight material to be used in spacecraft design is the driving force behind nanotube work at JSC. The molecular nature of these materials requires the appropriate tools for investigation of their structure, properties, and formation. The mechanism of nanotube formation is of particular interest because it may hold keys to controlling the formation of different types of nanotubes and allow them to be produced in much greater quantities at less cost than is currently available. This summer's work involved the interpretation of data taken last summer and analyzed over the academic year. The work involved diagnostic studies of carbon nanotube formation processes occurring in a laser-produced plume. Laser ablation of metal doped graphite to produce a plasma plume in which carbon nanotubes self assemble is one method of making carbon nanotube. The laser ablation method is amenable to applying the techniques of laser spectroscopy, a powerful tool for probing the energies and dynamics of atomic and molecular species. The experimental work performed last summer involved probing one of the metal catalysts, nickel, by laser induced fluorescence. The nickel atom was studied as a function of oven temperature, probe laser wavelength, time after ablation, and position in the laser produced plume. This data along with previously obtained data on carbon was analyzed over the academic year. Interpretations of the data were developed this summer along with discussions of future work. The temperature of the oven in which the target is ablated greatly influences the amount of material ablated and the propagation of the plume. The ablation conditions and the time scale of atomic and molecular lifetimes suggest that initial ablation of the metal doped carbon target results in atomic and small molecular species. The metal atoms survive for several milliseconds while the gaseous carbon atoms and small molecules nucleate more rapidly. Additional experiments and the development of in situ methods for carbon nanotube detection would allow these results to be interpreted from the perspective of carbon nanotube formation.

deBoer, Gary↗
NASA NTRS1987

Numerical solution of the radiation transport equation in disk geometry

An efficient numerical method for solving the problem of radiation transport in a dusty medium with two dimensional (2-D) disk geometry is described. It is a generalization of the one-dimensional quasi-diffusion method in which the transport equation is cast in diffusion form and then solved as a boundary value problem. The method should be applicable to a variety of astronomical sources, the dynamics of which are angular-momentum dominated and hence not accurately treated by spherical geometry, e.g., protoplanetary nebulae, circumstellar disks, interstellar molecular clouds, accretion disks, and disk galaxies. The computational procedure and practical considerations for implementing the method are described in detail. To illustrate the effects of 2-D radiation transport, some model results (dust temperature distributions and IR flux spectra) for externally heated, interstellar dust clouds with spherically symmetric and disk geometry are compared.

Spagna, George F., Jr.↗
NASA NTRS1950

The gas kinetics of very high flight speeds

The aerodynamic forces on bodies of arbitrary shape were investigated under conditions such that the mean free path of the air molecule is greater than the dimensions of the body. Air pressures and friction forces were calculated from gas kinetic theory for surfaces facing both toward and away from the air stream at any angle. Air forces for an atmosphere of definite composition (molecular hydrogen) were calculated as a function of the flight velocity. The results indicate that the friction stresses between the air and the body surface are of the same magnitude as the dynamic pressure and as the air pressures normal to the surface. The application of the general method to the specific cases such as thin airfoils and projectiles leads to high drag coefficients and poor glide ratios even for the theoretically best wing sections.

AIR FLOW VELOCITY - MOLECULAR ENERGIES↗
NASA NTRS2017

Validation of Test Methods for Air Leak Rate Verification of Spaceflight Hardware

As deep space exploration continues to be the goal of NASAs human spaceflight program, verification of the performance of spaceflight hardware becomes increasingly critical. Suitable test methods for verifying the leak rate of sealing systems are identified in program qualification testing requirements. One acceptable method for verifying the air leak rate of gas pressure seals is the tracer gas leak detector method. In this method, a tracer gas (commonly helium) leaks past the test seal and is transported to the leak detector where the leak rate is quantified. To predict the air leak rate, a conversion factor of helium-to-air is applied depending on the magnitude of the helium flow rate. The conversion factor is based on either the molecular mass ratio or the ratio of the dynamic viscosities. The current work was aimed at validating this approach for permeation-level leak rates using a series of tests with a silicone elastomer O-ring. An established pressure decay method with constant differential pressure was used to evaluate both the air and helium leak rates of the O-ring under similar temperature and pressure conditions. The results from the pressure decay tests showed, for the elastomer O-ring, that neither the molecular flow nor the viscous flow helium-to-air conversion factors were applicable. Leak rate tests were also performed using nitrogen and argon as the test gas. Molecular mass and viscosity based helium-to-test gas conversion factors were applied, but did not correctly predict the measured leak rates of either gas. To further this study, the effect of pressure boundary conditions was investigated. Often, pressure decay leak rate tests are performed at a differential pressure of 101.3 kPa with atmospheric pressure on the downstream side of the test seal. In space applications, the differential pressure is similar, but with vacuum as the downstream pressure. The same O-ring was tested at four unique differential pressures ranging from 34.5 to 137.9 kPa. Up to six combinations of upstream and downstream pressures for each differential pressure were compared. For a given differential pressure, the various combinations of upstream and downstream dry air pressures did not significantly affect the leak rate. As expected, the leak rate of the O-ring increased with increasing differential pressure. The results suggested that the current leak test pressure conditions, used to verify spacecraft sealing systems with elastomer seals, produce accurate values even though the boundary conditions do not model the space application.

helium leak detector↗
NASA NTRS2014

Modeling Materials: Design for Planetary Entry, Electric Aircraft, and Beyond

NASA missions push the limits of what is possible. The development of high-performance materials must keep pace with the agency's demanding, cutting-edge applications. Researchers at NASA's Ames Research Center are performing multiscale computational modeling to accelerate development times and further the design of next-generation aerospace materials. Multiscale modeling combines several computationally intensive techniques ranging from the atomic level to the macroscale, passing output from one level as input to the next level. These methods are applicable to a wide variety of materials systems. For example: (a) Ultra-high-temperature ceramics for hypersonic aircraft-we utilized the full range of multiscale modeling to characterize thermal protection materials for faster, safer air- and spacecraft, (b) Planetary entry heat shields for space vehicles-we computed thermal and mechanical properties of ablative composites by combining several methods, from atomistic simulations to macroscale computations, (c) Advanced batteries for electric aircraft-we performed large-scale molecular dynamics simulations of advanced electrolytes for ultra-high-energy capacity batteries to enable long-distance electric aircraft service; and (d) Shape-memory alloys for high-efficiency aircraft-we used high-fidelity electronic structure calculations to determine phase diagrams in shape-memory transformations. Advances in high-performance computing have been critical to the development of multiscale materials modeling. We used nearly one million processor hours on NASA's Pleiades supercomputer to characterize electrolytes with a fidelity that would be otherwise impossible. For this and other projects, Pleiades enables us to push the physics and accuracy of our calculations to new levels.

Supercomputing↗
NASA NTRS1986

Electron-Diffraction Analysis of Growth of GaAs

Report describes experiments that used reflection high-energy-electron diffraction (RHEED) to investigate behavior of GaAs surfaces during and after growth by molecular-beam epitaxy (MBE). Experimental results show dynamic RHEED measurements useful both as probes of surface and growth kinetics and as methods for determining and reproducing surface and growth conditions.

Lewis, B. F.↗
NASA NTRS1984

A variational method for the calculation of dynamic polarizabilities and two-photon transition moments - The dressed molecule approach

The solution of the time-dependent Schroedinger equation of the molecule + radiation field system is analyzed. A quantized radiation field is used. The relationship between the oscillatory wave function and the dressed state wave function is established, the oscillatory wave function being the solution if the radiation field is turned on at t = 0 and the dressed state being the stationary solution satisfying the boundary condition only if the field is present at t = negative infinity. In general, the oscillatory wave function can be expressed using a complete set of dressed states. However, in the presence of a nonresonant radiation field, the system is well represented by a single dressed state. Molecular properties such as dynamic polarizabilities and two-photon transition moments can be deduced from the dressed state wave function instead of the oscillatory wave function as in previous methods. Because of its stationary character, the dressed state is more amenable to approximate calculations. A CI method is developed for this purpose. The CI matrix is simple in form and can be readily constructed using existing computer codes. The present method can also be adapted to calculate other optical properties.

Huo, W. M.↗
NASA NTRS1997

Low-Density Aerodynamics of the Stardust Sample Return Capsule

The aerodynamics of the Stardust Sample Return Capsule are analyzed in the low- density, transitional flow regime using free-molecular, Direct Simulation Monte Carlo, Navier-Stokes, and Newtonian methods to provide inputs for constructing a transitional flow bridging relation. The accuracy of this bridging relation in reconstructing the aero- dynamic coefficients given by the more exact methods is presented for a range of flight conditions and vehicle attitudes. There is good agreement between the various prediction methods, and a simple sine-squared bridging relation is shown to provide a reasonably good description of the axial force, normal force, and pitching moment over a range of Knudsen numbers from 0.001 to 10. The predictions show a static instability of the Star- dust capsule in the free-molecular regime that persists well into the transitional flow. The addition of a thin disk to the base of the capsule is shown to remove this static instability. However, the extremely high entry velocity of 12.6 km/s for the proposed trajectory introduces difficult design issues for incorporating this disk caused by the high aerothermal loads that occur even under relatively rarefied conditions.

Wilmoth, Richard G.↗
NASA NTRS1992

Estimation of the viscosities of silicate liquids at high pressure from measurements of oxygen diffusivity

The dynamics and evolution of a magma ocean depend significantly on the structure and physical properties of silicate liquids at high pressure. Viscosity is particularly important because of its effect on the dynamics of convection and crystal setting. Traditionally, the viscosities of silicate liquids have been measured at high pressure by falling sphere viscometry but, due to some limitations of this technique, few measurements have been made at pressures exceeding 2.5 GPa. An alternative approach, which has previously been proven up to 2 GPa, is to use high pressure measurements of oxygen self-diffusivity (D) to estimate viscosity (eta) from the Eyring relationship eta = kT/D lambda (where k is the Boltzmann constant, T is absolute temperature, and lambda is the diffusive jump distance). We have performed preliminary experiments on Na2Si4O9 liquid in order to explore the feasibility of this method up to pressures in excess of 10 GPa. Diffusion couples were prepared from glass starting materials with one half of the sample enriched in O-18 and the other half containing the natural abundance (0.2 wt percent). Diffusion experiments were performed at 1600 to 1800 C at pressures in the range of 2.5 to 10 GPa for times up to 6 minutes using a 1200 ton multianvil apparatus. Oxygen diffusivities were derived from the resulting O-18 concentration profiles, which were analyzed using an ion microprobe. At 1800 C, oxygen diffusivities increase continuously from 1 x 10(exp -10) m(exp -2) s(exp -1) at 2.5 Gpa to 5 x 10(exp -10) m(exp -2) s(exp -1) at 10 GPa. From the Eyring relationship, the viscosity of Na2SiO9 liquid is predicted to decrease by approx. 0.7 log units over this pressure range at 1800 C. These trends, which can be related to changes in the structure of the liquid at high pressure, are in agreement with results of molecular dynamics calculations. These results suggest that it should be possible to estimate the viscosities of silicate liquids up to at least 15 GPa and 2200 C by this method.

Rubie, D. C.↗