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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Considerations for Femtosecond Laser Electronic Excitation Tagging in High-Speed Flows

Femtosecond laser electronic excitation tagging (FLEET) is an unseeded method for molecular tagging which offers valuable opportunities for measurement of high-speed (transonic, supersonic or hypersonic) flows. The unique nature of high-speed testing demands certain performance from FLEET such as satisfactory signal-to-noise ratio (SNR) at depressed static conditions (i.e., low temperatures, pressures and densities), wide dynamic range for velocity determination (especially single-shot), and measurements with acceptable accuracy and precision. This dissertation strives to evaluate FLEET in those regards and provide strategies to maximize the method's capabilities. A zero-dimensional kinetics model in nitrogen explains FLEET signal changes with pressure/density and/or temperature in terms of plasma-chemical reactions. Poorly known rate coefficients are tuned by comparing model output to measurements, with temporal agreement up to several hundred nanoseconds. Modeling reveals that initial signal peaks at reduced density because of slowed temporal evolution (and decay) of excited populations. Low temperatures enhance signal by enlarging cluster ion populations which contribute to excited species via electron-ion dissociative recombination. A purpose-built free jet facility provides experimental validation of the kinetics model and assesses FLEET velocimetry in low temperature and pressure/density conditions. Signal, lifetime, accuracy and precision results are obtained from unheated subsonic through Mach 4.0 operation of the facility, with best results noted. FLEET measurements of a sweeping jet (SWJ) actuator in compressible operation showcase its advantages in a highly unsteady jet containing subsonic through supersonic velocities. FLEET velocimetry is performed in the device's internal and external flow fields, with the latter compared to hot-wire anemometry. Internal measurements reveal the absence of shockwaves theorized to occur at high pressure ratios. Simultaneous qualitative measurements of compressible jet mixing are shown as a proof-of-concept. Overall, the work demonstrates that previous understanding of SWJ incompressible operation readily extends into the compressible realm. Practical aspects of performing FLEET velocimetry are detailed, along with strategies for improving measurement quality. Determination of a fundamental precision in nitrogen and air is attempted. Experiments show that increasing time delay and/or SNR improves velocimetry precision. A comparison of five camera systems indicates sensors with larger pixels capture higher SNR data and produce more precise results.

FLEET↗
NASA NTRSDate not supplied

Multiscale Analysis of Thermo-Mechanical Behavior of Boron Nitride-Reinforced Epoxy Nanocomposites

The effect of size, shape, morphology, and arrangement of micro constituents of Boron Nitrite (BN) nanoplatelet/epoxy composites on their properties were investigated using a multi scale approach that includes Molecular Dynamics (MD) and micromechanics. The thermo-mechanical properties of the composites were evaluated using molecular theory and the analysis showed that the elastic constants of BN/epoxy composites were not severely affected by temperature. Also, the micromechanical analysis of Generalized Method of Cells (GMC) was utilized at higher length scale to evaluate elastic properties of the composites, for different geometries and arrangements of micro constituents. The predicted results of the analysis showed that the size, aspect ratios, morphology, and the arrangements of inclusions in BN/epoxy nanocomposites all have remarkable effect on the mechanical performances of the material systems.

nanocomposites↗
NASA NTRSDate not supplied

Interfacial thermodynamics of cryogenic fluids: The effect on non-condensable gas on fluid storage

Propellant tanks that contain cryogenic fluids (CFs) in low gravity conditions are usually pressurized with non-condensable (NC) gases for fast extraction. Unfortunately, the presence of NC gases causes CFs to exhibit higher boil-offs compared to pure CF systems. For optimal utilization of the cryogenic fuels, these higher boil-off rates must be minimized. Our goal is to quantify the effects NC gases have on the evaporation and condensation dynamics of CFs to develop strategies that mitigate the higher boil-offs. The hypothesis is that the NC gas accumulates around the CF’s liquid vapor interface in the Knudsen layer, creating a kinetic barrier for both evaporation and condensation. Testing this hypothesis with experimental methods is difficult due to the transient nature of the Knudsen layer, which is only a few nanometers thick. Furthermore, experimental capabilities are limited in low gravity conditions and extremely expensive. Accordingly, we approach this problem from a theoretical perspective. In this study, we employ molecular dynamics (MD) simulations to probe the interfacial mechanisms that affect CF’s evaporation and condensation at varying concentrations of NC gases. Specifically, we use nitrogen (N2) and oxygen (O2) as our CFs and neon (Ne) as our NC gas. Using MD simulations, we show that Ne accumulates at N2 liquid vapor interface across a wide range of Ne concentrations, thereby impeding mass transport of N2. Our simulations allow for direct computation of the molar flux as well as the mass accommodation coefficient (MAC) which can then be used as input parameters to continuum fluid dynamics (CFD) models for optimal storage tank design.

Michael Robert DeLyser↗
NASA NTRS1986

Modeling of infrared flux spectra from disk-shaped interstellar dust clouds

Although there is growing observational and theoretical evidence for many disk-shaped objects of astrophysical interest, spherical geometry is assumed in most radiative transfer models. Recently we generalized the quasi-diffusion method developed by Leung (1975, 1976) for spherical geometry to solve the problem of scattering, absorption, and reemission by dust grains in a medium of disk geometry. The method is applicable to a variety of astronomical sources whose dynamics are angular-momentum dominated and hence not accurately treated by spherical geometry, e.g., protoplanetary nebulae, circumstellar disks, bipolar-flow molecular clouds, accretion disks and disk galaxies. Using this technique and realistic grain opacities we construct theoretical models to determine self-consistently the dust temperature distribution and infrared emission from disk-shaped, quiescent dark globules heated externally by the intersteller radiation field. The effects of the following parameters on the temperature structure and the emergent spectrum are studied: grain type (graphite and silicate), optical depth, density inhomogeneity, and degree of disk flattening. The disk models are characterized by an aspect ratio and an optical depth at 0.55 microns. For inhomogeneous models, a gaussian density distribution is assumed such that the ratio of central to surface density is 100. To study the effects of source geometry, we also compare results for models with spherical and disk geometry. In this case both models have the same radius and central optical depth, the disk models having a 1:1 aspect ratio. While the dust temperature distributions in the two cases are very similar, the emergent flux for the disk model depends sensitively on the viewing angle. For clouds which are unresolved, one would expect, since the thermal emission is isotropic in the neighborhood of an emitting grain, and since the emission (in the far infrared) is optically thin, that the emergent flux spectrum should be characteristic only of the dust temperature, and independent of viewing angle. However, because of the lack of complete symmetry and the resulting radiation anisotropy, this is found not to be the case for the disk models. Angle is that it implies a large uncertainty in estimating the radiation dust mass.

Spagna, G. F., Jr.↗
NASA NTRS1993

A fast recursive algorithm for molecular dynamics simulation

The present recursive algorithm for solving molecular systems' dynamical equations of motion employs internal variable models that reduce such simulations' computation time by an order of magnitude, relative to Cartesian models. Extensive use is made of spatial operator methods recently developed for analysis and simulation of the dynamics of multibody systems. A factor-of-450 speedup over the conventional O(N-cubed) algorithm is demonstrated for the case of a polypeptide molecule with 400 residues.

Jain, A.↗
NASA NTRS1983

Remote sensing by infrared heterodyne spectroscopy

The use of infrared heterodyne spectrocopy for the study of planetary atmospheres is discussed. Infrared heterodyne spectroscopy provides a convenient and sensitive method for measuring the true intensity profiles of atmospheric spectral lines. Application of radiative transfer theory to measured lineshapes can then permit the study of molecular abundances, temperatures, total pressures, excitation conditions, and dynamics of the regions of line formation. The theory of formation of atmospheric spectral lines and the retrieval of the information contained in these molecular lines is illustrated. Notable successes of such retrievals from infrared heterodyne measurements on Venus, Mars, Jupiter and the Earth are given. A discussion of developments in infrared heterodyne technology is also presented.

Kostiuk, T.↗
NASA NTRS1983

Remote sensing by IR heterodyne spectroscopy

The use of infrared heterodyne spectroscopy for the study of planetary atmospheres is discussed. Infrared heterodyne spectroscopy provides a convenient and sensitive method for measuring the true intensity profiles of atmospheric spectral lines. Application of radiative transfer theory to measured lineshapes can then permit the study of molecular abundances, temperatures, total pressures, excitation conditions, and dynamics of the regions of line formation. The theory of formation of atmospheric spectral lines and the retrieval of the information contained in these molecular lines is illustrated. Notable successes of such retrievals from infrared heterodyne measurements on Venus, Mars, Jupiter and the earth are given. A discussion of developments in infrared heterodyne technology is also presented. Previously announced in STAR as N83-28551

Kostiuk, T.↗
NASA NTRSDate not supplied

Aerodynamic Implications of Aerocapture Systems for Uranus Orbiters

Exploration of the Ice Giants, Uranus in particular, has been indicated as the highest priority new flagship-class mission by the 2022 Planetary Science Decadal Survey. However, due to it’s location in the solar system, extensive travel time and resources are required for a mission to reach Uranus. Aerocapture offers shorter interplanetary transit times and reduced propellant requirements for such a mission when compared to a traditional propulsive-capture orbiter. Considerations for the aerodynamics of an aerocapture system enabling flagship-class science mission to Uranus are presented. The applicability of a low lift-to-drag entry vehicle aeroshell with flight heritage to a new planet and atmosphere is evaluated. A preliminary aerodatabase is formulated from low-fidelity methods and heritage flight mission data to provide initial estimates of trajectory space to be considered for inserting a spacecraft into orbit at Uranus using aerocapture. Computational fluid dynamics calculations in the free molecular, transitional, and continuum regimes using Direct Simulation Monte Carlo and Navier-Stokes solutions are used to assess the validity of this preliminary aerodatabase. The considered entry vehicle is found to varying levels of agreement with work done by previous flight missions when evaluated for flight in Uranus’ atmosphere. Implications for the aerodynamics of this vehicle as applied in an aerocapture system in flight at Uranus as found by the current work are discussed.

Eli R. Shellabarger↗
NASA NTRSDate not supplied

Aerodynamic Implications of Aerocapture Systems for Uranus Orbiters

Exploration of the Ice Giants, Uranus in particular, has been indicated as the highest priority new flagship-class mission by the 2022 Planetary Science Decadal Survey. However, due to it's location in the solar system, extensive travel time and resources are required for a mission to reach Uranus. Aerocapture offers shorter interplanetary transit times and reduced propellant requirements for such a mission when compared to a traditional propulsive-capture orbiter. Considerations for the aerodynamics of an aerocapture system enabling flagship-class science mission to Uranus are presented. The applicability of a low lift-to-drag entry vehicle aeroshell with flight heritage to a new planet and atmosphere is evaluated. A preliminary aerodatabase is formulated from low-fidelity methods and heritage flight mission data to provide initial estimates of trajectory space to be considered for inserting a spacecraft into orbit at Uranus using aerocapture. Computational fluid dynamics calculations in the free molecular, transitional, and continuum regimes using Direct Simulation Monte Carlo and Navier-Stokes solutions are used to assess the validity of this preliminary aerodatabase. The considered entry vehicle is found to varying levels of agreement with work done by previous flight missions when evaluated for flight in Uranus' atmosphere. Implications for the aerodynamics of this vehicle as applied in an aerocapture system in flight at Uranus as found by the current work are discussed.

Uranus↗
NASA NTRSDate not supplied

Semi-Empirical Interatomic Potential for Large-Scale Molecular Dynamics Simulation of Metal-Oxide Systems

Classical molecular dynamics (MD) simulation can be applied to systems containing billions of atoms during times up to microseconds. However, utilization of a large-scale MD simulation requires reliable but computationally cheap interatomic potentials. In the case of metallic systems, embedded atom method (EAM) and Finnis-Sinclair (FS) potentials are the reasonable choices, but their development for multicomponent alloys is a challenge. Addition of oxygen atoms makes the problem of developing semi-empirical potentials even more difficult because the long-range Coulomb interaction cannot be neglected, and different atoms have different charges which vary during MD relaxation such that a charge equilibration procedure should be applied every MD step. As results researchers tend to avoid simulating metal-oxide systems. In this talk, a new Finnis-Sinclair (FS) Ni-Co-Cr potential will be presented. Special attention will be paid to reproducing of the element melting temperatures and clustering. Next, a charge transfer ionic potential (CTIP) developed to simulate the NiO properties will be presented. This potential utilizes the Ni FS potential and accounts for the Coulomb interaction in the vicinity of oxygen atoms. Fitting of the CTIP parameters and FS potential functions describing the interaction with oxygen atoms to the available experimental and ab initio data will be discussed. Results of MD simulation of interaction of NiO particles with dislocations in Ni will be shown.

molecular dynamics simulation↗
NASA NTRS2018

VLTI-GRAVITY Measurements of Cool Evolved Stars I. Variable Photosphere and Extended Atmosphere of the Mira Star R Peg

Context. Dynamic model atmospheres of Mira stars predict variabilities in the photospheric radius and in atmospheric molecular layers which are not yet strongly constrained by observations. Aim. Here we measure the variability of the oxygen-rich Mira star R Peg in near-continuum and molecular bands. Methods. We used near-infrared K-band spectro-interferometry with a spectral resolution of about 4000 obtained at four epochs between post-maximum and minimum visual phases employing the newly available GRAVITY beam combiner at the Very Large Telescope Interferometer (VLTI). Results. Our observations show a continuum radius that is anti-correlated with the visual lightcurve. Uniform disc (UD) angular diameters at a near-continuum wavelength of 2.24 m are steadily increasing with values of 8.7 0.1 mas, 9.4 0.1 mas, 9.8 0.1 mas, and 9.9 0.1 mas at visual phases of 0.15, 0.36, 0.45, 0.53, respectively. UD diameters at a bandpass around 2.05 m, dominated by water vapour, follow the near-continuum variability at larger UD diameters between 10.7 mas and 11.7 mas. UD diameters at the CO 2-0 bandhead, instead, are correlated with the visual lightcurve and anti-correlated with the near-continuum UD diameters, with values between 12.3 mas and 11.7 mas. Conclusions. The observed anti-correlation between continuum radius and visual lightcurve is consistent with an earlier study of the oxygen-rich Mira S Lac, and with recent 1D CODEX dynamic model atmosphere predictions. The amplitude of the variation is comparable to the earlier observations of S Lac, and smaller than predicted by CODEX models. The wavelength-dependent visibility variations at our epochs can be reproduced by a set of CODEX models at model phases between 0.3 and 0.6. The anti-correlation of water vapour and CO contributions at our epochs suggests that these molecules undergo different processes in the extended atmosphere along the stellar cycle. The newly available GRAVITY instrument is suited to conducting longer time series observations, which are needed to provide strong constraints on the model-predicted intra- and inter-cycle variability.

photospheric radius↗
NASA NTRS1993

Calculation of molecular excitation rates

State-to-state collisional excitation rates for interstellar molecules observed by radio astronomers continue to be required to interpret observed line intensities in terms of local temperatures and densities. A problem of particular interest is collisional excitation of water which is important for modeling the observed interstellar masers. In earlier work supported by a different NASA Grant, excitation of water in collisions with He atoms was studied; after many years of successively more refined calculations that problem now seems to be well understood, and discrepancies with earlier experimental data for related (pressure broadening) phenomena are believed to reflect experimental errors. Because of interstellar abundances, excitation by H2, the dominant interstellar species, is much more important than excitation by He, although it has been argued that rates for excitation by these are similar. Under the current grant theoretical study of this problem has begun which is greatly complicated by the additional degrees of freedom which must be included both in determining the interaction potential and also in the molecular scattering calculation. We have now computed the interaction forces for nearly a thousand molecular geometries and are close to having an acceptable global fit to these points which is necessary for the molecular dynamics calculations. Also, extensive modifications have been made to the molecular scattering code, MOLSCAT. These included coding the rotational basis sets and coupling matrix elements required for collisions of an asymmetric top with a linear rotor. A new method for numerical solution of the coupled equations has been incorporated. Because of the long-ranged nature of the water-hydrogen interaction it is necessary to integrate the equations to rather large intermolecular separations, and the integration methods previously available in MOLSCAT are not ideal for such cases. However, the method used by Alexander in his HIBRIDON code is particularly suited for such cases. We have obtained this code and incorporated that part which solves the coupled differential equations as an option in the MOLSCAT program.

Flynn, George↗
NASA NTRS1999

Multiple Point Dynamic Gas Density Measurements Using Molecular Rayleigh Scattering

A nonintrusive technique for measuring dynamic gas density properties is described. Molecular Rayleigh scattering is used to measure the time-history of gas density simultaneously at eight spatial locations at a 50 kHz sampling rate. The data are analyzed using the Welch method of modified periodograms to reduce measurement uncertainty. Cross-correlations, power spectral density functions, cross-spectral density functions, and coherence functions may be obtained from the data. The technique is demonstrated using low speed co-flowing jets with a heated inner jet.

Seasholtz, Richard↗
NASA NTRS1989

Determinations of molecular weight and molecular weight distribution of high polymers by the rheological properties

Several methods are reviewed by which the molecular weight (MW) and the molecular weight distribution (MWD) of polymeric material were determined from the rheological properties. A poly(arylene ether) polymer with six different molecular weights was used in this investigation. Experimentally measured MW and MWD were conducted by GPC/LALLS (gel permeation chromatography/low angle laser light scattering), and the rheological properties of the melts were measured by a Rheometric System Four rheometer. It was found that qualitative information of the MW and MWD of these polymers could be derived from the viscoelastic properties, with the methods proposed by Zeichner and Patel, and by Dormier et al., by shifting the master curves of the dynamic storage modulus, G', and the loss modulus, G'', along the frequency axis. Efforts were also made to calculate quantitative profiles of MW and MWD for these polymers from their rheological properties. The technique recently proposed by Wu was evaluated. It was found that satisfactory results could only be obtained for polymers with single modal distribution in the molecular weight.

Huang, J. Y.↗
NASA NTRS2020

Molecular Dynamics Investigation of the Structural and Mechanical Properties of Off-Stoichiometric Epoxy Resins

We carried out molecular dynamics (MD) simulations to measure the mechanical properties of various off-stoichiometric polymers regarding amine to epoxy ratios (r) and to understand the stiffness of the polymers in terms of their structures. The aerospace-grade API-60 epoxy resin is used as an adhesive bond for assembling large-scale composite structures via the co-curing-ply bonding method. This method will produce a reliable and certifiable composite joint without additional fasteners. Calculated Young's modulus was measured from the uniaxial tension simulation with several high strain rates, and the experimental modulus was estimated by extrapolating the simulation results. We found that the stiffness was associated with molecular packing caused by chemical cross-linking. We also found that the number of network clusters gradually decreased as the ratio approached r = 1.0, which made the tighter cluster and the system much stiffer with an increase in the molecular weight and the degree of cross-linking. Structural properties such as Rg, MSD were measured to figure out the degree of stiffness with respect to the r.

Jang, Chang Woon↗
NASA NTRS1971

The molecular dynamics of atmospheric reaction

Detailed information about the chemistry of the upper atmosphere took the form of quantitative data concerning the rate of reaction into specified states of product vibration, rotation and translation for exothermic reaction, as well as concerning the rate of reaction from specified states of reagent vibration, rotation and translation for endothermic reaction. The techniques used were variants on the infrared chemiluminescence method. Emphasis was placed on reactions that formed, and that removed, vibrationally-excited hydroxyl radicals. Fundamental studies were also performed on exothermic reactions involving hydrogen halides.

Polanyi, J. C.↗
NASA NTRS1999

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗
NASA NTRS2000

Modeling Molecules

The molecule modeling method known as Multibody Order (N) Dynamics, or MBO(N)D, was developed by Moldyn, Inc. at Goddard Space Flight Center through funding provided by the SBIR program. The software can model the dynamics of molecules through technology which stimulates low-frequency molecular motions and properties, such as movements among a molecule's constituent parts. With MBO(N)D, a molecule is substructured into a set of interconnected rigid and flexible bodies. These bodies replace the computation burden of mapping individual atoms. Moldyn's technology cuts computation time while increasing accuracy. The MBO(N)D technology is available as Insight II 97.0 from Molecular Simulations, Inc. Currently the technology is used to account for forces on spacecraft parts and to perform molecular analyses for pharmaceutical purposes. It permits the solution of molecular dynamics problems on a moderate workstation, as opposed to on a supercomputer.

Source record↗