Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular modelling”

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 433 records · Page 24

Rotational relaxation of molecular hydrogen at moderate temperatures

Using a coupled rotation-vibration-dissociation model the rotational relaxation times for molecular hydrogen as a function of final temperature (500-5000 K), in a hypothetical scenario of sudden compression, are computed. The theoretical model is based on a master equation solver. The bound-bound and bound-free transition rates have been computed using a quasiclassical trajectory method. A review of the available experimental data on the rotational relaxation of hydrogen is presented, with a critical overview of the method of measurements and data reduction, including the sources of errors. These experimental data are then compared with the computed results.

COUPLED ROTATION VIBRATION DIS↗

Extinction of a Free Methanol Droplet in Microgravity

The combustion of isolated, unsupported methanol droplets under microgravity conditions has been studied in the NASA-Lewis 2.2-sec drop tower. The substitution of O2/He mixtures for O2/N2 mixtures results in a substantial increase in the oxygen index required for ignition, the rate of gasification, and the diameter at extinction. As a result, droplets with initial size larger than 1 mm can be produced, deployed, ignited, and burned to extinction, all under microgravity conditions in this facility. A comprehensive time dependent numerical model incorporating detailed multicomponent molecular transport mechanisms and complex elementary chemical kinetics was developed and applied to simulate the entire transient droplet combustion processes (ignition/quasi-steady burning/extinction). Model outputs show good agreement with the experimental measurements, without parameter adjustment. Analysis of modeling results shows that the thermal conductivity of the gas phase under combustion environments increased by a factor of two when the 50-percent O2/50 percent He mixture is used instead of air. The increased thermal conductivity accelerates the heat transfer from the flame front to the droplet surface, resulting in a factor of two increase in the burning rate.

Cho, S. Y.↗

Origin of the Eumetazoa: testing ecological predictions of molecular clocks against the Proterozoic fossil record

Molecular clocks have the potential to shed light on the timing of early metazoan divergences, but differing algorithms and calibration points yield conspicuously discordant results. We argue here that competing molecular clock hypotheses should be testable in the fossil record, on the principle that fundamentally new grades of animal organization will have ecosystem-wide impacts. Using a set of seven nuclear-encoded protein sequences, we demonstrate the paraphyly of Porifera and calculate sponge/eumetazoan and cnidarian/bilaterian divergence times by using both distance [minimum evolution (ME)] and maximum likelihood (ML) molecular clocks; ME brackets the appearance of Eumetazoa between 634 and 604 Ma, whereas ML suggests it was between 867 and 748 Ma. Significantly, the ME, but not the ML, estimate is coincident with a major regime change in the Proterozoic acritarch record, including: (i) disappearance of low-diversity, evolutionarily static, pre-Ediacaran acanthomorphs; (ii) radiation of the high-diversity, short-lived Doushantuo-Pertatataka microbiota; and (iii) an order-of-magnitude increase in evolutionary turnover rate. We interpret this turnover as a consequence of the novel ecological challenges accompanying the evolution of the eumetazoan nervous system and gut. Thus, the more readily preserved microfossil record provides positive evidence for the absence of pre-Ediacaran eumetazoans and strongly supports the veracity, and therefore more general application, of the ME molecular clock.

NASA Discipline Evolutionary Biology↗

Optical properties of grains in molecular clouds and accretion disks

A baseline model of the composition and abundances of grains and gases in molecular cloud cores and accretion disks around young stars is defined by employing: a wide range of astronomical data and theory; the composition of primitive bodies in the solar system; and solar elemental abundances. It is proposed that in the coldest portions of these objects the major grain species include amorphous olivine, amorphous orthopyroxene, volatile and refractory organics, water ice, troilite, and metallic iron. Using a combination of laboratory measurements of optical constants and asymptotic theory, values of the real and imaginary indices of refraction of these grain species over a wavelength range that runs from the vacuum UV to the radio domain are derived. Auxiliary information for these grain species, such as their vaporization temperatures bulk densities, and their fractional abundances by mass. The above information on grain properties in molecular cloud cores and accretion disks is used to estimate the Rosseland mean opacity of the grains in both environments and the IR and microwave opacity of grains in accretion disks.

Pollack, J. B.↗

Molecular processes in a high temperature shock layer

Models of the shock layer encountered by an Aeroassisted Orbital Transfer Vehicle require as input accurate cross sections and rate constants for the atomic and molecular processes that characterize the shock radiation. From the estimated atomic and molecular densities in the shock layer and the expected residence time of 1 m/s, it can be expected that electron-ion collision processes will be important in the shock model. Electron capture by molecular ions followed by dissociation, e.g., O2(+) + e(-) yields 0 + 0, can be expected to be of major importance since these processes are known to have high rates (e.g., 10 to the -7th power cu/cm/sec) at room temperature. However, there have been no experimental measurements of dissociative recombination (DR) at temperatures ( 12000K) that are expected to characterize the shock layer. Indeed, even at room temperature, it is often difficult to perform experiments that determine the dependence of the translational energy and quantum yields of the product atoms on the electronic and vibrational state of the reactant molecular ions. Presented are ab initio quantum chemical studies of DR for molecular ions that are likely to be important in the atmospheric shock layer.

Guberman, S. L.↗

Infrared-pumped third-harmonic and sum-frequency generation in diatomic molecules

Using molecular hydrogen as a model system, we present a calculation which demonstrates the feasibility of efficient third-harmonic and sum-frequency generation of IR laser radiation by diatomic molecular gases. A conversion efficiency of 10% without phase matching is achievable with a modest laser intensity of about 20 MW/sq cm. In addition to a two-photon resonance, the efficient conversion depends on a mechanism which makes use of virtual vibronic transitions exclusively. Possible methods for phase matching are also suggested.

She, C. Y.↗

Discrete Rotational Energy for Polyatomic Molecules in Direct Simulation Monte Carlo

The environments experienced by vehicles upon entry into planetary atmospheres generally result in severe aerothermodynamic loading. These flows are characterized by the formation of strong shock waves, behind which high temperatures and non-equilibrium environments are generated, and the accurate prediction of internal energy relaxation and chemical kinetics becomes important. This manuscript focuses on the internal relaxation of the rotational energy of polyatomic molecules. Previous work has described the implementation of a discrete rotational energy model for diatomic molecules (linear rotators). Later, Gimelshein et. al extended this model to include polyatomic molecular internal energies. However, they did not include details on how to generate the particles at equilibrium conditions and centrifugal distortion was not included in their model. The purpose of the present study is to extend the afore-mentioned models to the treatment of discrete rotational energy in polyatomic molecules with the addition of centrifugal distortion. An overview of the background theory needed for the implementation of the model will be given followed by the model itself and verification of the model. For example, sampled versus Boltzmann distributions are compared in Figure 1 for CH3 (oblate rotator – two quantum numbers) and examples of thermal relaxation behavior for linear and non-linear rotators are presented in Figure 2. In addition, comparisons of thermophysical properties will be presented, as well as examples of adiabatic relaxation and application to hypersonic flow.

Rarefied Gas Dynamics↗

The influence of ozone and aerosols on the brightness and color of the twilight zone

The radiance and color of the twilight sky are calculated for single scattered radiation with the use of spherically symmetric models of the earth's atmosphere. Spherical geometry is used throughout the calculations with no plane parallel approximations. Refraction effects are taken into account through fine subdivision of the atmosphere into spherical shells of fixed index of refraction. Shell's law of refraction is used to calculate a direction of travel each time that a photon traverses the interface between layers. Five different models of the atmosphere were used: a pure molecular scattering atmosphere; molecular atmosphere plus ozone absorption; and three models with aerosol concentrations of 1, 3, and 10 times normal together with molecular scattering and ozone absorption. The results of the calculations are shown for various observation positions and local viewing angles in the solar plane for wavelengths in the range of 0.40 microns to 0.75 microns.

Adams, C. N.↗

The influence of ozone and aerosols on the brightness and color of the twilight sky

The radiance and color of the twilight sky are calculated for single scattered radiation with the use of spherically symmetric models of the earth's atmosphere. Spherical geometry is used throughout the calculations with no plane-parallel approximations. Refraction effects are taken into account through fine subdivision of the atmosphere into spherical shells of fixed index of refraction. Snell's law of refraction is used to calculate a new direction of travel each time that a photon traverses the interface between layers. Five different models of the atmosphere were used: a pure molecular scattering atmosphere; molecular atmosphere plus ozone absorption; and three models with aerosol concentrations of one, three, and ten times normal together with molecular scattering and ozone absorption. The results of the calculations are shown for various observation positions and local viewing angles in the solar plane for wavelengths in the range from 0.40 to 0.75 micron.

Adams, C. N.↗

Development of a Computational Chemical Vapor Deposition Model: Applications to Indium Nitride and Dicyanovinylaniline

A computational chemical vapor deposition (CVD) model is presented, that couples chemical reaction mechanisms with fluid dynamic simulations for vapor deposition experiments. The chemical properties of the systems under investigation are evaluated using quantum, molecular and statistical mechanics models. The fluid dynamic computations are performed using the CFD-ACE program, which can simulate multispecies transport, heat and mass transfer, gas phase chemistry, chemistry of adsorbed species, pulsed reactant flow and variable gravity conditions. Two experimental setups are being studied, in order to fabricate films of: (a) indium nitride (InN) from the gas or surface phase reaction of trimethylindium and ammonia; and (b) 4-(1,1)dicyanovinyl-dimethylaminoaniline (DCVA) by vapor deposition. Modeling of these setups requires knowledge of three groups of properties: thermodynamic properties (heat capacity), transport properties (diffusion, viscosity, and thermal conductivity), and kinetic properties (rate constants for all possible elementary chemical reactions). These properties are evaluated using computational methods whenever experimental data is not available for the species or for the elementary reactions. The chemical vapor deposition model is applied to InN and DCVA. Several possible InN mechanisms are proposed and analyzed. The CVD model simulations of InN show that the deposition rate of InN is more efficient when pulsing chemistry is used under conditions of high pressure and microgravity. An analysis of the chemical properties of DCVA show that DCVA dimers may form under certain conditions of physical vapor transport. CVD simulations of the DCVA system suggest that deposition of the DCVA dimer may play a small role in the film and crystal growth processes.

Cardelino, Carlos↗

Highlights of Transient Plume Impingement Model Validation and Applications

This paper describes highlights of an ongoing validation effort conducted to assess the viability of applying a set of analytic point source transient free molecule equations to model behavior ranging from molecular effusion to rocket plumes. The validation effort includes encouraging comparisons to both steady and transient studies involving experimental data and direct simulation Monte Carlo results. Finally, this model is applied to describe features of two exotic transient scenarios involving NASA Goddard Space Flight Center satellite programs.

Woronowicz, Michael↗

The sugar model: catalysis by amines and amino acid products

Ammonia and amines (including amino acids) were shown to catalyze the formation of sugars from formaldehyde and glycolaldehyde, and the subsequent conversion of sugars to carbonylcontaining products under the conditions studied (pH 5.5 and 50 degrees C). Sterically unhindered primary amines were better catalysts than ammonia, secondary amines, and sterically hindered primary amines (i.e. alpha-aminoisobutyric acid). Reactions catalyzed by primary amines initially consumed formaldehyde and glycolaldehyde about 15-20 times faster than an uncatalyzed control reaction. The amine-catalyzed reactions yielded aldotriose (glyceraldehyde), ketotriose (dihydroxyacetone), aldotetroses (erythrose and threose), ketotetrose (erythrulose), pyruvaldehyde, acetaldehyde, glyoxal, pyruvate, glyoxylate, and several unindentified carbonyl products. The concentrations of the carbonyl products, except pyruvate and ketotetrose, initially increased and then declined during the reaction, indicating their ultimate conversion to other products (like larger sugars or pyruvate). The uncatalyzed control reaction yielded no pyruvate or glyoxylate, and only trace amounts of pyruvaldehyde, acetaldehyde and glyoxal. In the presence of 15 mM catalytic primary amine, such as alanine, the rates of triose and pyruvaldehyde of synthesis were about 15-times and 1200-times faster, respectively, than the uncatalyzed reaction. Since previous studies established that alanine is synthesized from glycolaldehyde and formaldehyde via pyruvaldehyde as its direct precursor, the demonstration that the alanine catalyzes the conversion of glycolaldehyde and formaldehyde to pyruvaldehyde indicates that this synthetic pathway is capable of autocatalysis. The relevance of this synthetic process, named the Sugar Model, to the origin of life is discussed.

Acetaldehyde/analogs & derivatives/chemistry↗

Shock propagation and the generation of magnetohydrodynamic wave fields in inhomogeneous molecular clouds

We develop a simple one-dimensional model for the interaction of a steady, thin, planar shock wave with a nonrigid cloud which may be in motion relative to the surrounding medium, and we apply the model to shocks impinging on, and propagating through, molecular clouds. Both 'adiabatic' (gamma = 5/3) and radiative (gamma = 1) shocks are considered and we allow for the presence of a uniform magnetic field directed either parallel or perpendicular to the shock normal. The former field orientation is equivalent to the hydrodynamic case, and the latter involves only fast MHD shocks. We focus on the manner in which such shocks can generate internal kinetic motions in the cloud on a range of size and density scales through the direct acceleration of cores and clumps by shocks transmitted into them and through the generation of an MHD wavefield via the reflection of the incident shock at clump boundaries. We find that stronger incident Mach numbers and smaller density contrasts lead to more efficient cloud acceleration, as do isothermal intercloud shocks and small intercloud magnetic field strengths. The acceleration efficiency is insensitive to the adiabatic index and the magnetic field strength in the cloud itself. For typical parameter choices, the direct acceleration of clouds and clumps by strong shocks is found to be substantial and could at least in part account for their observed velocity dispersions. If the shocks are moderately weak, the final velocity of the cloud is linearly related to its initial velocity, with higher acceleration giving shallower slopes (i.e., final velocity distributions which are less sensitive to the initial distribution). Compared to the kinetic energy of the postshock cloud, the energy given to the wavefield at each encounter is small, and the heating of the interclump medium by the dissipation of this wavefield is found to be insufficient to balance the cooling rate in the cloud as a whole (although it may be important in particular regions), even if this medium is warm, unless it is also extremely tenuous (n approximately less than 0.1/cu cm). Nevertheless, the correction for the velocity imparted to the cloud leads to a substantial increase in the critical incident Mach numner for wave emission over that reported by Spitzer for the rigid case. The implications of our model for shock-induced star formation are discussed briefly.

Miesch, Mark S.↗

New Insights into the Chemical Composition of Five Oort Cloud Comets After Re-analysis of their Infrared Spectra

We present revised results for the main molecular species in five Oort Cloud (OC) comets observed with NIRSPEC at the Keck Observatory between 1999 and 2012 (C/1999 S4 (LINEAR), C/2001 A2 (LINEAR), C/2007 W1 (Boattini), C/2012 F6 (Lemmon), and C/2012 S1 (ISON)). The re-evaluation of these data shows the improvement of results in some of the datasets, in particular for comets observed and analyzed before the advent of new and revised fluorescence models and terrestrial retrieval methods introduced since 2011. We observe significant improvements in the resulting rotational temperatures and production rates for all species, and in mixing ratios of minor species (relative to water). The re-analysis also allowed us to quantify species not analyzed previously, mostly due to the lack of molecular models (e.g., ammonia and formaldehyde). We note, however, that the improvement of these revised values is less substantial for comets observed (and/or analyzed) since 2010.

Manuela Lippi↗

A 2000 Solar Mass Rotating Molecular Disk Around NGC 6334A

We present millimeter and centimeter wave spectroscopic observations of the H II region NGC 6334A. We have mapped the source in several transitions of CO, CS, and NH3. The molecular emission shows a distinct flattened structure in the east-west direction. This structure is probably a thick molecular disk or torus (2.2 x 0.9 pc) responsible for the bipolarity of the near-infrared (NIR) and radio continuum emission which extends in two "lobes" to the north and south of the shell-like H II region. The molecular disk is rotating from west to east (omega approximately equals 2.4 km/s.pc) about an axis approximately parallel to the radio and NIR emission lobes. By assuming virial equilibrium, we find that the molecular disk contains approximately 2000 solar mass. Single-component gas excitation model calculations show that the molecular gas in the disk is warmer and denser (T(sub k) approximately equals 60 K, n approximately equals 3000/cc) than the gas to the north and south (T(sub k) approximately equals 50 K, n approximately equals 400/cc). High resolution (approximately 5 sec) NH3 (3, 3) images of NGC 6334A reveal several small (approximately 0.1 pc) clumps, one of which lies southwest of the radio continuum shell, and is spatially coincident with a near-infrared source, IRS 20. A second NH3 clump is coincident with an H2O maser and the center of a molecular outflow. The dense gas tracers, CS J = 5 approaches 4 and 7 approaches 6, peak near IRS 20 and the H2O maser, not at NGC 6334A. IRS 20 has a substantial far-infrared (FIR) luminosity L(sub FIR) approximately 10(exp 5) solar luminosity, which indicates the presence of an O 7.5 star but has no detected radio continuum (F(sub 6 cm) < 0.02 Jy). The combination of dense gas, a large FIR luminosity and a lack of radio continuum can best be explained if IRS 20 is a protostar. A third clump of NH3 emission lies to the west of IRS 20 but is not associated with any other molecular or continuum features. The star formation activity in the region has moved west of NGC 6334A to IRS 20 and the H2O maser position. We suggest that NGC 6334A, IRS 20, and the H2O maser spot are part of a "protocluster" of stars which is condensing from the massive molecular disk. The similarity between the structure around NGC 6334A and other large (r approximately 1 pc), massive (M approximately 10(exp 3) solar mass), rotating disks (K3-50A and G10.6-0.4) suggests that this may be a common mechanism by which open clusters form.

Kraemer, Kathleen E.↗

Effect of molecular anisotropy on the intensity and degree of polarization of light scattered from model atmospheres

Computations of the intensity, flux, degree of polarization, and the positions of neutral points are presented for models of the terrestrial gaseous and hazy atmospheres by incorporating the molecular anisotropy due to air in the Rayleigh scattering optical thickness and phase matrix. Molecular anisotropy causes significant changes in the intensity, flux and the degree of polarization of the scattered light. The positions of neutral points do not change significantly. When the Rayleigh scattering optical thickness is kept constant and the molecular anisotropy factor is included only in the Rayleigh phase matrix, the flux does not change and the intensity and positions of neutron points change by a small amount. The changes in the degree of polarization are still significant.

Bahethi, O. P.↗

Scalable Variable Charge Molecular Dynamics Simulations of Metal-Oxide Systems

Interfaces between metals and oxides are an important feature in many technologically relevant materials, e.g., oxidation of metal surfaces, oxide-dispersion strengthened (ODS) alloys, dielectric components, and thermal barrier coatings among others. Experimental studies of such interfaces are challenging since the majority are buried within the bulk, making computational modeling an attractive alternative. Molecular dynamics (MD) simulations operate at the length scales relevant to many interface-mediated mechanisms, but the requisite interatomic potentials for metal-oxide systems require computationally expensive variable charge schemes to account for the disparate bonding types, thus often limiting their effectiveness. Here we introduce several improvements to the charge transfer interatomic potential (CTIP) model which enable greater computational efficiency for large scale MD simulations. Then, using a new CTIP parametrization for the Ni-O system, we demonstrate its capabilities to capture critical atomic scale mechanisms associated with metal-oxide interfaces. Long time scale simulations (>10 ns) are used to investigate high temperature oxidation and oxide precipitation from the melt, and large length scale simulations (> 1 million atoms) are used to study the interaction of dislocations with oxide particles. We have implemented the new CTIP model in the widely used, open-source MD code LAMMPS.

Gabriel Plummer↗

Results of the SPAS-01 RCS plume impingement test

The Orbiter RCS flowfield is modeled in the free molecular region and an impingement analysis is performed to characterize on-orbit plume interactions with nearby bodies, such as the SPAS-01 satellite. The flow in the far field is treated with a source flow model at any point 40 ft from the nozzle. Impingement is calculated according to Newtonian impact theory. The flow density in the far field is expressed in terms of a conservation of energy equation, with account taken of the expansion of the gas, cooling, and the angular distribution. The SPAS-01 platform gathered data at eight different positions around the Orbiter to test the model predictions. Excellent agreement was obtained in predictions of the torque and rotation rate induced on the SPAS by firing the thrusters.

Lazaron, M. P.↗