Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupling reaction”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Development of Comprehensive Reduced Kinetic Models for Supersonic Reacting Shear Layer Simulations

Large-scale simulations of multi-dimensional unsteady turbulent reacting flows with detailed chemistry and transport can be computationally extremely intensive even on distributed computing architectures. With the development of suitable reduced chemical kinetic models, the number of scalar variables to be integrated can be decreased, leading to a significant reduction in the computational time required for the simulation with limited loss of accuracy in the results. A general MATLAB-based automated mechanism reduction procedure is presented to reduce any complex starting mechanism (detailed or skeletal) with minimal human intervention. Based on the application of the quasi steady-state (QSS) approximation for certain chemical species and on the elimination of the fast reaction rates in the mechanism, several comprehensive reduced models, capable of handling different fuels such as C2H4, CH4 and H2, have been developed and thoroughly tested for several combustion problems (ignition, propagation and extinction) and physical conditions (reactant compositions, temperatures, and pressures). A key feature of the present reduction procedure is the explicit solution of the concentrations of the QSS species, needed for the evaluation of the elementary reaction rates. In contrast, previous approaches relied on an implicit solution due to the strong coupling between QSS species, requiring computationally expensive inner iterations. A novel algorithm, based on the definition of a QSS species coupling matrix, is presented to (i) introduce appropriate truncations to the QSS algebraic relations and (ii) identify the optimal sequence for the explicit solution of the concentration of the QSS species. With the automatic generation of the relevant source code, the resulting reduced models can be readily implemented into numerical codes.

Zambon, A. C.↗
NASA NTRS2016

P21C-2113: Constraining Hesperian Martian PCO2 from Mineral Analysis at Gale Crater

Carbon dioxide is an essential atmospheric component in martian climate models that attempt to reconcile a faint young sun with widespread evidence of liquid water at the planet's surface in the Noachian and Early Hesperian. Current estimates of ancient martian CO levels, derived from global inventories of carbon, and orbital detections of Noachian and Early Hesperian clay mineralbearing terrains indicate CO levels that are unable to support warm and wet conditions. These estimates are subject to various sources of uncertainty however. Mineral and contextual sedimentary environmental data collected by the Mars Science Laboratory rover Curiosity in Gale Crater provide a more direct means of estimating the atmospheric partial pressure of CO (P ) coinciding with a long-lived lake system in Gale crater at approximately 3.5 Ga. Results from a reaction transport model, which simulates mineralogy observed within the Sheepbed member at Yellowknife Bay by coupling mineral equilibria with carbonate precipitation kinetics and rates of sedimentation, indicate atmospheric levels in the 10's mbar range. At such low P levels, climate models are unable to warm Hesperian Mars anywhere near the freezing point of water and other gases are required to raise atmospheric pressure to prevent lakes from boiling away. Thus, lacustrine features of Gale formed in a cold environment by a mechanism yet to be determined, or the climate models lack an essential component that would serve to elevate surface temperatures, at least temporally and/or locally, on Hesperian Mars. Our results also impose restrictions on the potential role of atmospheric CO in inferred warmer conditions of the Noachian.

Bristow, Thomas↗
NASA NTRS2016

Constraining Hesperian martian PCO2 from mineral analysis at Gale crater

Carbon dioxide is an essential atmospheric component in martian climate models that attempt to reconcile a faint young sun with planet-wide evidence of liquid water at the planets surface in the Noachian and Early Hesperian. Current estimates of ancient martian CO2 levels, derived from global inventories of carbon, and orbital detections of Noachian and Early Hesperian clay mineral-bearing terrains indicate CO2 levels that are unable to support warm and wet conditions. These estimates are subject to various sources of uncertainty however. Mineral and contextual sedimentary environmental data collected by the Mars Science Laboratory rover Curiosity in Gale Crater provide a more direct means of estimating the atmospheric partial pressure of CO2 (PCO2) coinciding with a long-lived lake system in Gale crater at ~3.5 Ga. Results from a reaction-transport model, which simulates mineralogy observed within the Sheepbed member at Yellowknife Bay by coupling mineral equilibria with carbonate precipitation kinetics and rates of sedimentation, indicate atmospheric PCO2 levels in the 10’s mbar range. At such low PCO2 levels, climate models are unable to warm Hesperian Mars anywhere near the freezing point of water and other gases are required to raise atmospheric pressure to prevent lakes from boiling away. Thus, lacustrine features of Gale formed in a cold environment by a mechanism yet to be determined, or the climate models still lack an essential component that would serve to elevate surface temperatures, at least temporally and/or locally, on Hesperian Mars. Our results also impose restrictions on the potential role of atmospheric CO2 in inferred warmer conditions of the Noachian.

Bristow, T. F.↗
NASA NTRS1998

Evolution of Interstellar Grains

During the past two decades observations combined with laboratory simulations, have revolutionized our understanding of interstellar ice and dust, the raw materials from which planets, comets and stars form. Most interstellar material is concentrated in large molecular clouds where simple molecules are formed by dust-grain and gas-phase reactions. Gaseous species striking the cold (10K) dust stick, forming an icy grain mantle. This accretion, coupled with UV photolysis, produces a complex chemical mixture containing volatile, non-volatile, and isotopically fractionated species. Ices in molecular clouds contain the very simple molecules H2O, CH3OH, CO, CO2, H2, and perhaps some NH3 and H2CO, as well as more complex species. The evidence for these compounds, as well as carbon-rich materials, will be reviewed and the possible connections with comets and meteorites will be presented in the first part of the talk . The second part of the presentation will focus on interstellar/precometary ice photochemical evolution and the species likely to be found in comets. The chemical composition and photochemical evolution of realistic interstellar/pre-cometary ice analogs will be discussed. Ultraviolet photolysis of these ices produces H2, H2CO, CO2, CO, CH4, HCO, and more complex molecules. When ices representative of interstellar grains and comets are exposed to UV radiation at low temperature a series of moderately complex organic molecules are formed in the ice including: CH3CH2OH (ethanol), HC(=O)NH2 (formamide), CH3C(=O)NH2 (acetamide), and R-C=N (nitriles). Several of these are already known to be in the interstellar medium, and their presence indicates the importance of grain processing. After warming to room temperature an organic residue remains. This is composed primarily of hexamethylenetetramine (HMT, C6H12N4), with lesser amounts of polyoxymethylene-related species (POMs), amides, and ketones. This is in sharp contrast to the organic residues produced by irradiating unrealistic interstellar ice analogs or thermally promoted polymerization-type reactions in unirradiated realistic ice mixtures.

Allamandola, Lou J.↗
NASA NTRS1999

Aerosols and Particulates Workshop Sampling Procedures and Venues Working Group Summary

The Sampling Procedures and Venues Workgroup discussed the potential venues available and issues associated with obtaining measurements. Some of the issues included Incoming Air Quality, Sampling Locations, Probes and Sample Systems. The following is a summary of the discussion of the issues and venues. The influence of inlet air to the measurement of exhaust species, especially trace chemical species, must be considered. Analysis procedures for current engine exhaust emissions regulatory measurements require adjustments for air inlet humidity. As a matter of course in scientific investigations, it is recommended that "background" measurements for any species, particulate or chemical, be performed during inlet air flow before initiation of combustion, if possible, and during the engine test period as feasible and practical. For current regulatory measurements, this would be equivalent to setting the "zero" level for conventional gas analyzers. As a minimum, it is recommended that measurements of the humidity and particulates in the incoming air be taken at the start and end of each test run. Additional measurement points taken during the run are desirable if they can be practically obtained. It was felt that the presence of trace gases in the incoming air is not a significant problem. However, investigators should consider the ambient levels and influences of local air pollution for species of interest. Desired measurement locations depend upon the investigation requirements. A complete investigation of phenomenology of particulate formation and growth requires measurements at a number of locations both within the engine and in the exhaust field downstream of the nozzle exit plane. Desirable locations for both extractive and in situ measurements include: (1) Combustion Zone (Multiple axial locations); (2) Combustor Exit (Multiple radial locations for annular combustors); (3) Turbine Stage (Inlet and exit of the stage); (4) Exit Nozzle (Multiple axial locations downstream of the nozzle). Actual locations with potential for extractive or non-intrusive measurements depend upon the test article and test configuration. Committee members expressed the importance of making investigators aware of various ports that could allow access to various stages of the existing engines. Port locations are engine si)ecific and might allow extractive sampling or innovative hybrid optical-probe access. The turbine stage region was one the most desirable locations for obtaining samples and might be accessed through boroscope ports available in some engine designs. Discussions of probes and sampling systems quickly identified issues dependent on particular measurement quantities. With general consensus, the group recommends SAE procedures for measurements and data analyses of currently regulated exhaust species (CO2, CO, THC, NO(x),) using conventional gas sampling techniques. Special procedures following sound scientific practices must be developed as required for species and/or measurement conditions not covered by SAE standards. Several issues arose concerning short lived radicals and highly reactive species. For conventional sampling, there are concerns of perturbing the sample during extraction, line losses, line-wall reactions, and chemical reactions during the sample transport to the analyzers. Sample lines coated with quartz.or other materials should be investigated for minimization of such effects. The group advocates the development of innovative probe techniques and non-intrusive optical techniques for measurement of short lived radicals and highly reactive species that cannot be sampled accurately otherwise. Two innovative probe concepts were discussed. One concept uses specially designed probes to transfer optical beams to and from a region of flow inaccessible by traditional ports or windows. The probe can perturb the flow field but must have a negligible impact on the region to be optically sampled. Such probes are referred to as hybrid probes and are under development at AEDC for measurement in the high pressure, high temperature of a combustor under development for power generation. The other concept consists of coupling an instrument directly to the probe. The probe would isolate a representative sample stream, freeze chemical reactions and direct the sample into the analyzer portion of the probe. Thus, the measurement would be performed in situ without sample line losses due either to reactions or binding at the wall surfaces. This concept was used to develop a fast, in situ, time-of-flight mass spectrometer measurement system for temporal quantification of NO in the IMPULSE facility at AEDC. Additional work is required in this area to determine the best probe and sampling technique for each species measurement requirement identified by the Trace Chemistry Working Group. A partial list of Venues was used as a baseline for discussion. Additional venues were added to the list and the list was broken out into the following categories: (1)Engines (a) Sea Level Test Stands (b) Altitude Chambers; (2) Annular Combustor Test Stands, (3) Sector Flametube Test Stands, (4) Fundamentals Rigs/Experiments.

Pachlhofer, Peter↗
NASA NTRS2002

Gas-Phase Reaction Pathways and Rate Coefficients for the Dichlorosilane-Hydrogen and Trichlorosilane-Hydrogen Systems

As part of NASA Ames Research Center's Integrated Process Team on Device/Process Modeling and Nanotechnology our goal is to create/contribute to a gas-phase chemical database for use in modeling microelectronics devices. In particular, we use ab initio methods to determine chemical reaction pathways and to evaluate reaction rate coefficients. Our initial studies concern reactions involved in the dichlorosilane-hydrogen (SiCl2H2--H2) and trichlorosilane-hydrogen (SiCl2H-H2) systems. Reactant, saddle point (transition state), and product geometries and their vibrational harmonic frequencies are determined using the complete-active-space self-consistent-field (CASSCF) electronic structure method with the correlation consistent polarized valence double-zeta basis set (cc-pVDZ). Reaction pathways are constructed by following the imaginary frequency mode of the saddle point to both the reactant and product. Accurate energetics are determined using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations (CCSD(T)) extrapolated to the complete basis set limit. Using the data from the electronic structure calculations, reaction rate coefficients are obtained using conventional and variational transition state and RRKM theories.

Dateo, Christopher E.↗
NASA NTRS1985

Pressure dependence of the absolute rate constant for the reaction Cl + C2H2 from 210-361 K

In recent years, considerable attention has been given to the role of chlorine compounds in the catalytic destruction of stratospheric ozone. However, while some reactions have been studied extensively, the kinetic data for the reaction of Cl with C2H2 is sparse with only three known determinations of the rate constant k3. The reactions involved are Cl + C2H2 yields reversibly ClC2H2(asterisk) (3a) and ClC2H2(asterisk) + M yields ClC2H2 + M (3b). In the present study, flash photolysis coupled with chlorine atomic resonance fluorescence have been employed to determine the pressure and temperature dependence of k3 with the third body M = Ar. Room temperature values are also reported for M = N2. The pressure dependence observed in the experiments confirms the expectation that the reaction involves addition of Cl to the unsaturated C2H2 molecule followed by collisional stabilization of the resulting adduct radical.

Brunning, J.↗
NASA NTRS1982

The global troposphere - Biogeochemical cycles, chemistry, and remote sensing

The chemical composition of the troposphere is controlled by various biogeochemical cycles that couple the atmosphere with the oceans, the solid earth and the biosphere, and by atmospheric photochemical/chemical reactions. These cycles and reactions are discussed and a number of key questions concerning tropospheric composition and chemistry for the carbon, nitrogen, oxygen and sulfur species are identified. Next, various remote sensing techniques and instruments capable of measuring and monitoring tropospheric species from the ground, aircraft and space to address some of these key questions are reviewed. Future thrusts in remote sensing of the troposphere are also considered.

Levine, J. S.↗
NASA NTRS1988

Electrochemical Growth Of Crystals In Gels

Nucleation and growth rates readily controlled. Technique developed to grow crystals by controlling rate of transfer of one component into crystallization volume. Method involves electrochemically controlled generation of one of precipitation species, coupled with diffusion barrier. New procedure, developed in connection with formation of lead tin telluride by reaction in gels of metal ions with telluride ions.

Barber, Patrick G.↗
NASA NTRS1996

LES, DNS, and RANS for the Analysis of High-Speed Turbulent Reacting Flows

A filtered density function (FDF) method suitable for chemically reactive flows is developed in the context of large eddy simulation. The advantage of the FDF methodology is its inherent ability to resolve subgrid scales (SGS) scalar correlations that otherwise have to be modeled. Because of the lack of robust models to accurately predict these correlations in turbulent reactive flows, simulations involving turbulent combustion are often met with a degree of skepticism. The FDF methodology avoids the closure problem associated with these terms and treats the reaction in an exact manner. The scalar FDF approach is particularly attractive since it can be coupled with existing hydrodynamic computational fluid dynamics (CFD) codes.

Colucci, P. J.↗
NASA NTRS1995

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO rev. reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X=H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate(DELTA)H(sup o)(sub f) value for HNO using state of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗
NASA NTRSDate not supplied

Multiscale Modeling of the Uptake of Hydrogen Chloride on Anodized Aluminum in Relevance to Spacecraft Fire Safety

Fire safety is a critical component of space missions. An electrical fire within the spacecraft often produces gaseous hydrogen chloride (HCl) from the smoldering of wire insulation. These gases tend to stick to the interior materials of the spacecraft, which is often comprised of treated aluminum. Determining the rate of uptake for HCl by aluminum helps the safe design of the spacecraft to reduce catastrophic mission risk. Our previous studies showed that a one-step global reaction model for adsorption of HCl on aluminum is incapable of predicting the experimentally measured HCl uptake of anodized aluminum, compared to other treatments with thinner oxide layers. X-ray photoelectron spectroscopy revealed that aluminum chloride is found deep in the oxide layer of the anodized aluminum, implying that the transport of HCl within the oxide layer covering the aluminum surface must be considered. This work uses a multiscale modeling approach to predict the uptake of HCl. A pore model is developed wherein the HCl reacts with the aluminum oxide pore walls to create an aluminum chloride layer. The HCl then diffuses through that layer for further reaction, until the diffusion resistance grows too large, and the surface is said to be saturated. The pore model is then coupled with the reactor-scale model to predict overall HCl uptake. Results show that this new multiscale approach predicts the uptake of HCl far more accurately than previously used single-scale models.

Modeling↗
NASA NTRS1995

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO reversible reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HNO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X double bond H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate (Delta)H(sub f, sup o) value for HNO using state-of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +/- 0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗
NASA NTRS1980

Analysis of solid propellant combustion in a closed vessel including secondary reaction

A theory for combustion of solid propellants in a closed vessel is presented allowing for residual exothermic chemical reaction in the bulk of the gas in the vessel. Particular attention is given to propellants exhibiting thick gaseous flame zones such as nitrocellulose, double-base and nitramine propellants. For these, the reaction at high pressures is assumed to involve mainly the oxidation of residual hydrocarbons by NO. It is shown that the direct dynamic coupling between the exothermicity, the molecular weight reduction and the changing pressure can influence the dp/dt-p traces obtained, in a manner not directly related to mass burning rate of the solid. Energy and species conservation equations are derived for the bulk of the vessel in differential form; the system is solved numerically. The results show the effect of extended chemical reaction upon measurable combustion characteristics such as dp/dt-p and burn rate pressure exponent, demonstrating its potential importance in interpretation of closed vessel firing data, depending on the pace of the residual gas phase reactions.

Benreuven, M.↗
NASA NTRS1977

Tropospheric ozone - Coupling transport and photochemistry

The budget of tropospheric ozone is re-examined in the light of lower observed nitrogen oxide densities and revised reaction rates. To estimate the relative importance of transport and photochemistry, model calculations are presented which determine one-dimensional vertical ozone profiles at midlatitudes for equinoctial conditions. These calculations imply that both photochemistry and transport are important in controlling the abundance of tropospheric ozone. It is estimated that photochemical processes produce about 0.6 times 10 to the eleventh and destroy about 1.4 times 10 to the eleventh tropospheric ozone molecules per sq cm per sec, while mixing processes transport about 1.4 times 10 to the eleventh ozone molecules per sq cm per sec from the stratosphere into the troposphere and ground level fluxes destroy about 0.6 times 10 to the eleventh ozone molecules per sq cm per sec. The sensitivity of our model calculations to transport parameters, key rate coefficients, and nitrogen oxide densities are discussed.

Chameides, W. L.↗
NASA NTRS1998

An Atmospheric General Circulation Model with Chemistry for the CRAY T3E: Design, Performance Optimization and Coupling to an Ocean Model

The design, implementation and performance optimization on the CRAY T3E of an atmospheric general circulation model (AGCM) which includes the transport of, and chemical reactions among, an arbitrary number of constituents is reviewed. The parallel implementation is based on a two-dimensional (longitude and latitude) data domain decomposition. Initial optimization efforts centered on minimizing the impact of substantial static and weakly-dynamic load imbalances among processors through load redistribution schemes. Recent optimization efforts have centered on single-node optimization. Strategies employed include loop unrolling, both manually and through the compiler, the use of an optimized assembler-code library for special function calls, and restructuring of parts of the code to improve data locality. Data exchanges and synchronizations involved in coupling different data-distributed models can account for a significant fraction of the running time. Therefore, the required scattering and gathering of data must be optimized. In systems such as the T3E, there is much more aggregate bandwidth in the total system than in any particular processor. This suggests a distributed design. The design and implementation of a such distributed 'Data Broker' as a means to efficiently couple the components of our climate system model is described.

Farrara, John D.↗
NASA NTRS1996

Interstellar/Precometary Organic Material and the Photochemical Evolution of Complex Organics

During the past two decades ground-, air-, and space-based infrared spectroscopic observations, combined with realistic laboratory simulations, have revolutionized our understanding of interstellar ice and dust, the raw materials from which planets, comets and stars form. Most interstellar material is concentrated in Large molecular clouds where simple molecules are formed by dust grain and gas phase reactions. Gaseous species striking the cold (10 K) dust will stick, forming an icy grain mantle. This accretion, coupled with energetic particle bombardment and UV photolysis, will produce a complex chemical mixture containing volatile, non-volatile, and isotopically fractionated species. Ices in molecular clouds contain the very simple molecules H2O, CH3OH, CO, CO2, H2, and perhaps some NH3 and H2CO, as well as more complex species including nitriles and ketones or esters. The evidence for these compounds as well as carbon rich materials such as polycyclic aromatic hydrocarbons (PAHs), microdiamonds, and amorphous carbon will be reviewed and the possible connections with comets and meteorites will be presented in the first part of the talk. The second part of the presentation will focus on interstellar/precometary ice photochemical evolution. The chemical composition and photochemical evolution of realistic interstellar/pre-cometary ice analogs containing methanol will be discussed. ultraviolet photolysis of these ices produces H2, H2CO, CO2, CO, CH4, HCO, and more complex molecules. Infrared spectroscopy, H-1 and C-13 nuclear magnetic resonance (NMR) spectroscopy, and gas chromatography-mass spectrometry demonstrate that when ices representative of interstellar grains and comets are exposed to UV radiation at low temperature a series of moderately complex organic molecules are formed in the ice including: CH3CH2OH (ethanol), HC(=O)NH2 (formamide), CH3C(=O)NH2 (acetamide), and R-C(integral)N (nitriles). Several of these are already known to be in the interstellar medium, and their presence indicates the importance of grain processing. After warming to room temperature what remains is an organic residue composed primarily of Hexamethylenetetramine (HMT, C6H12N4), with lesser amounts of polyoxymethylene related species (POMs), amides, and ketones. This is in sharp contrast to the organic residues produced by irradiating ices which do not contain methanol (unrealistic interstellar ice analogs) or thermally promoted polymerization-type reactions in unirradiated realistic ice mixtures. Here HMT is only a minor product in a residue dominated by a mixture of polyoxymethylene related species. The implications, for infrared astronomy and astrochemistry, of high concentrations of HMT in interstellar and cometary ices may be profound. The ultraviolet photolysis of HMT frozen in H20 ice produces the "XCN" band observed in the spectra of protostellar objects and laboratory ices, as well as carbon oxides and other nitriles. Thus, HMT may be a precursor of XCN in protostellar objects and a source of CN and CO in the tail of comets. Also, HMT is known to hydrolyze under acidic conditions to yield ammonia and formaldehyde as well as amino acids. Thus, HMT may have been a source of organic material delivered to the early earth by comets.

Allamandola, Lou J.↗
NASA NTRS2007

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.↗