Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Molecular modeling of interactions in electronic nose sensors for environmental monitoring

We report a study aimed at understanding analyte interactions with sensors made from polymer-carbon black composite films. The sensors are used in an Electronic Nose (ENose) which is used for monitoring the breathing air quality in human habitats. The model mimics the experimental conditions of the composite film deposition and formation and was developed using molecular modeling and simulation tools. The Dreiding 2.21 Force Field was used for the polymer and analyte molecules while graphite parameters were assigned to the carbon black atoms. The polymer considered for this work is methyl vinyl ether / maleic acid copolymer. The target analytes include both inorganic (NH3) and organic (methanol) types of compound. Results indicate different composite-analyte interaction behavior.

electronic noise polymer-carbon composite molecula↗
NASA NTRSDate not supplied

Effect of Aromaticity on Soot Formation: A Reactive Molecular Dynamics Study

This study investigates the formation of soot from n-dodecane and its blend with m-xylene under pyrolytic and low-oxygen conditions with reactive molecular dynamics (MD) simulations. ReaxFF force field is used to simulate the interactions of fuel and oxygen molecules at engine-relevant conditions (2500-3500 K and 70-75 bar). Fuel molecules with carbon atom density of 0.03 g/cm 3 are dispersed into a cube and constant temperature simulations are carried out with a Nosé-Hoover thermostat. Pyrolysis of n-dodecane and the chemical pathway towards the formation of aromatic hydrocarbons are investigated. The initial pyrolysis rate of n-dodecane is observed to be slightly higher with the addition of aromatic compounds and oxygen. While the first aromatics formed from the decomposed n-dodecane are phenyl radicals at 2500 K, decomposition of the aliphatic fuel compound is observed to be necessary to build a carbon cluster, and eventually soot. Morphological properties of the created soot particles are also analyzed. The maturity, size and aromatics content of the soot particle are slightly affected with the addition of m-xylene and oxygen at 3500 K within the 3 ns simulation time possibly due the formation of larger and more condensed aromatic rings, resulting in larger soot molecules. The outcomes of this study provide important insights into the formation mechanism of soot and its morphological characteristics depending on varying fuel concentrations at pyrolysis and low oxygen engine conditions. Results of this analysis will further contribute to the understanding of contrail formation on generated soot particles originated from the combustion of different chemical components with particle chemistry matching companion experiments.

soot formation↗
NASA NTRS2014

Computational and Experimental Investigation of Li-doped Ionic Liquid Electrolytes: [pyr14][tfsi], [pyr13][fsi], and [EMIM][BF4]

We employ molecular dynamics (MD) simulation and experiment to investigate the structure, thermodynamics, and transport of N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N -methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li-salt mole fraction (0.05 xLi+ 0.33) and temperature (298 K T 393 K). Structurally, Li+ is shown to be solvated by three anion neigh- bors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi+ we find the presence of lithium aggregates. Pulsed field gradient spin-echo NMR measurements of diffusion and electrochemical impedance spectroscopy measurements of ionic conductivity are made for the neat ionic liquids as well as 0.5 molal solutions of Li-salt in the ionic liquids. Bulk ionic liquid properties (density, diffusion, viscosity, and ionic conductivity) are obtained with MD and show excellent agreement with experiment. While the diffusion exhibits a systematic decrease with increasing xLi+, the contribution of Li+ to ionic conductivity increases until reach- ing a saturation doping level of xLi+ 0.10. Comparatively, the Li+ conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1-0.3 mScm. Our transport results also demonstrate the necessity of long MD simulation runs ( 200 ns) required to converge transport properties at room T. The differences in Li+ transport are reflected in the residence times of Li+ with the anions (Li), which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, to comment on the relative kinetics of Li+ transport in each liquid, we find that while the net motion of Li+ with its solvation shell (vehicular) significantly contributes to net diffusion in all liquids, the importance of transport through anion exchange (hopping) increases at high xLi+ and in liquids with large anions.

electrolytes↗
NASA NTRS2006

Orion Aerodynamics for Hypersonic Free Molecular to Continuum Conditions

Numerical simulations are performed for the Orion Crew Module, previously known as the Crew Exploration Vehicle (CEV) Command Module, to characterize its aerodynamics during the high altitude portion of its reentry into the Earth's atmosphere, that is, from free molecular to continuum hypersonic conditions. The focus is on flow conditions similar to those that the Orion Crew Module would experience during a return from the International Space Station. The bulk of the calculations are performed with two direct simulation Monte Carlo (DSMC) codes, and these data are anchored with results from both free molecular and Navier-Stokes calculations. Results for aerodynamic forces and moments are presented that demonstrate their sensitivity to rarefaction, that is, for free molecular to continuum conditions (Knudsen numbers of 111 to 0.0003). Also included are aerodynamic data as a function of angle of attack for different levels of rarefaction and results that demonstrate the aerodynamic sensitivity of the Orion CM to a range of reentry velocities (7.6 to 15 km/s).

Moss, James N.↗
NASA NTRS1974

Turbulent mixing of multiple co-axial gaseous fuel jets in a supersonic airstream

An experimental and analytical investigation of a strut-mounted, four-nozzle, downstream-facing, gaseous fuel-injector assembly was conducted in a Mach 4 airstream at 154 psia and 520 R. Helium was used as the injectant in order to simulate a low molecular weight fuel, and the interjet spacing (S/D) was the main parameter varied. The principal data are in the form of helium concentration profiles at seven axial stations, Mach number distributions at three axial stations, and Schlieren photographs of the flow field at different interjet spacings. An approximate analysis was developed based upon a linearization in the Von Mises plane, Crocco integrals for the temperature and concentrations fields, and an eddy viscosity model. Good agreement with the data was achieved. It was found that, at these conditions, interjet spacing has a significant effect on mixing only up to S/D of approximately 2.5, with the concentration of the injectant at a given point becoming smaller as S/D is increased.

Lorber, A. K.↗
NASA NTRS1977

Silicon material task. Part 3: Low-cost silicon solar array project

The feasibility of a process for carbon reduction of low impurity silica in a plasma heat source was investigated to produce low-cost solar-grade silicon. Theoretical aspects of the reaction chemistry were studied with the aid of a computer program using iterative free energy minimization. These calculations indicate a threshold temperature exists at 2400 K below which no silicon is formed. The computer simulation technique of molecular dynamics was used to study the quenching of product species.

Roques, R. A.↗
NASA NTRS1984

Line emission from H II blister models

Numerical techniques to calculate the thermal and geometric properties of line emission from H II 'blister' regions are presented. It is assumed that the density distributions of the H II regions are a function of two dimensions, with rotational symmetry specifying the shape in three-dimensions. The thermal and ionization equilibrium equations of the problem are solved by spherical modeling, and a spherical sector approximation is used to simplify the three-dimensional treatment of diffuse ionizing radiation. The global properties of H II 'blister' regions near the edges of a molecular cloud are simulated by means of the geometry/density distribution, and the results are compared with observational data. It is shown that there is a monotonic increase of peak surface brightness from the i = 0 deg (pole-on) observational position to the i = 90 deg (edge-on) position. The enhancement of the line peak intensity from the edge-on to the pole-on positions is found to depend on the density, stratification, ionization, and electron temperature weighting. It is found that as i increases, the position of peak line brightness of the lower excitation species is displaced to the high-density side of the high excitation species.

Rubin, R. H.↗
NASA NTRS1985

Laplace-transform technique for deriving thermodynamic equations from the classical microcanonical ensemble

A direct and convenient method is presented for deriving expressions which equate any thermodynamic state function to averages of specific dynamical functions and their fluctuations over the classical microcanonical distribution. Specific expressions are obtained for a variety of thermodynamic quantities. The effect of various entropy definitions on the results are assessed, and the latter are compared to previous work in the literature. The derived formulas are applied to the analysis of molecular-dynamics computer simulations.

Pearson, E. M.↗
NASA NTRS1989

Particle kinetic simulation of high altitude hypervelocity flight

A new model for simulating the transfer of energy between the translational and rotational modes was derived for a homogeneous gas of diatomic molecules. The model was developed specifically for use in discrete particle simulation methods where molecular motion and inter-molecular collisions are treated at the molecular level. In such methods it is normal to assume a constant rotational collision number for the entire flowfield. The new model differs in that a temperature dependence is introduced which has been predicted by theory and observed in experiment. The new model is applied to the relaxation of rotational temperature, and is found to produce significant differences in comparison with the model normally employed at both high and low temperatures. Calculations have also been performed for a Mach 7 normal shock wave. Large differences in the solutions are again observed, with the new model offering an improved correspondence to the available experimental data.

Boyd, Iain D.↗
NASA NTRS1990

Rotational-translational energy transfer in rarefied nonequilibrium flows

A new model for simulating the transfer of energy between the translational and rotational modes is derived for a homogeneous gas of diatomic molecules. The model has been developed specifically for use in discrete particle simulation methods where molecular motion and intermolecular collisions are treated at the molecular level. A temperature dependence is introduced which has been predicted by theory and observed in experiment. The new model is applied to the relaxation of rotational temperature, and is found to produce significant differences in comparison with the model normally employed at both high and low temperatures. Calculations have also been performed for a Mach 7 normal shock wave.

Boyd, Iain D.↗
NASA NTRS1992

Edge technique - Theory and application to the lidar measurement of atmospheric wind

The paper describes the theory of the edge technique, a powerful method for the detection and measurement of small frequency shifts. It can be employed with a lidar to obtain range-resolved measurements of wind with high accuracy and high vertical resolution. The technique can be applied to measure wind with a lidar by using either the aerosol or molecular backscattered signal. Simulations for a ground-based lidar at 1.06 micron using reasonable instrumental parameters show an accuracy of the vector components of the wind which is better than 0.5 m/s from the ground to an altitude of 20 km for a 100-m vertical resolution and a 100-shot average.

Korb, C. L.↗
NASA NTRS1993

Isomerization reaction dynamics and equilibrium at the liquid-vapor interface of water. A molecular-dynamics study

The gauche-trans isomerization reaction of 1,2-dichloroethane at the liquid-vapor interface of water is studied using molecular-dynamics computer simulations. The solvent bulk and surface effects on the torsional potential of mean force and on barrier recrossing dynamics are computed. The isomerization reaction involves a large change in the electric dipole moment, and as a result the trans/gauche ratio is considerably affected by the transition from the bulk solvent to the surface. Reactive flux correlation function calculations of the reaction rate reveal that deviation from the transition-state theory due to barrier recrossing is greater at the surface than in the bulk water. This suggests that the system exhibits non-Rice-Ramsperger-Kassel-Marcus behavior due to the weak solvent-solute coupling at the water liquid-vapor interface.

Benjamin, Ilan↗
NASA NTRS1991

Molecular Dynamics study of Pb overlayer on Cu(100)

Isothermal-isobaric Molecular Dynamics (MD) simulation of a submonolayer Pb film in c(2x2) ordered structure adsorbed on a Cu(100) substrate showed retention of order to high T. The Embedded Atom Method (EAM) calculated the energy of atoms of overlayer and substrate. The time-averaged squared modulus of the two dimensional structure factor for the Pb overlayer measured the order of the overlayer. The results are for increasing T only, and require verification by simulated cooling.

Karimi, M.↗
NASA NTRS1992

Structure of turbulent non-premixed flames modeled with two-step chemistry

Direct numerical simulations of turbulent diffusion flames modeled with finite-rate, two-step chemistry, A + B yields I, A + I yields P, were carried out. A detailed analysis of the turbulent flame structure reveals the complex nature of the penetration of various reactive species across two reaction zones in mixture fraction space. Due to this two zone structure, these flames were found to be robust, resisting extinction over the parameter ranges investigated. As in single-step computations, mixture fraction dissipation rate and the mixture fraction were found to be statistically correlated. Simulations involving unequal molecular diffusivities suggest that the small scale mixing process and, hence, the turbulent flame structure is sensitive to the Schmidt number.

Chen, J. H.↗
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 NTRS1995

Molecular dynamics studies of simple membrane-water interfaces: Structure and functions in the beginnings of cellular life

Molecular dynamics computer simulations of the structure and functions of a simple membrane are performed in order to examine whether membranes provide an environment capable of promoting protobiological evolution. Our model membrane is composed of glycerol 1-monooleate. It is found that the bilayer surface fluctuates in time and space, occasionally creating thinning defects in the membrane. These defects are essential for passive transport of simple ions across membranes because they reduce the Born barrier to this process by approximately 40%. Negative ions are transferred across the bilayer more readily than positive ions due to favorable interactions with the electric field at the membrane-water interface. Passive transport of neutral molecules is, in general, more complex than predicted by the solubility-diffusion model. In particular, molecules which exhibit sufficient hydrophilicity and lipophilicity concentrate near membrane surfaces and experience 'interfacial resistance' to transport. The membrane-water interface forms an environment suitable for heterogeneous catalysis. Several possible mechanisms leading to an increase of reaction rates at the interface are discussed. We conclude that vesicles have many properties that make them very good candidates for earliest protocells. Some potentially fruitful directions of experimental and theoretical research on this subject are proposed.

Pohorille, Andrew↗
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 NTRS1998

Magellan Aerodynamic Characteristics During the Termination Experiment Including Thruster Plume-Free Stream Interaction

Results are presented on the aerodynamic characteristics of the Magellan spacecraft during the October 1994 Termination Experiment, including the effects of the thruster engine exhaust plumes upon the molecular free stream around the spacecraft and upon the aerodynamics coefficients. As Magellan passed through the Venusian atmosphere, the solar arrays were turned in opposite directions relative to the free stream creating a torque on the spacecraft. The spacecraft control system was programmed to counter the effects of this torque with attitude control engines to maintain an inertially fixed attitude. The orientation and reaction engine telemetry returned from Magellan are used to create a model of the aerodynamic torques. Geometric models of the Magellan spacecraft are analyzed with the aid of both free molecular and Direct Simulation Monte Carlo codes. The simulated aerodynamic torques determined are compared to the measured torques. The Direct Simulation Monte Carlo method is also used to model the attitude engine exhaust plumes, the free stream disturbance caused by these plumes, and the resulting torques acting on the spacecraft compared to no-exhaust plume cases. The effect of the exhaust plumes was found to be sufficiently large that thrust reversal is possible.

Cestero, Francisco J.↗