Bayesian estimation of the hydroxyl radical diffusion coefficient at low temperature and high pressure from atomistic molecular dynamics
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
We have published recently the effect of dissimilar diffusion coefficients on the size of the voltammetric waves from a mixture of two redox-active compounds. Similarly, at the potential range where three redox-active species, decamethylferrocene (dMeFc), ferrocene (Fc) and N-methylphenothiazine (MePTZ), are oxidized simultaneously with rates controlled by linear diffusion, electrogenerated radicals diffusing outwards from the electrode react with the original species diffusing towards the electrode from the bulk; thus, Fc(+) reacts with dMeFc producing Fc and dMeFc(+), while MePTZ(+) reacts both with dMeFc producing MePTZ and dMeFc(+), and with Fc producing MePTZ and Fc(+). These reactions replace dMeFc with Fc at the second plateau, and both dMeFc and Fc with MePTZ at the third plateau. Since the diffusion coefficients of the three species are not equal, the mass-transfer limited currents of the second and the third oxidation wave plateaus change by approx. 10%. Numerical simulations of the experimental voltamograms support this mechanism. Similar results were also obtained for a mixture of four redoxactive compounds. The implications of this non-additive nature of currents on: (a) the use of internal voltammetric standards for quantitative analysis of a mixture of redox-active compounds; and, (b) the half wave potentials (E1/2) of the 2nd, 3rd and 4th waves for qualitative analysis, will be discussed.
At the potential range where both decamethylferrocene (dMeFc) and ferrocene (Fc) are oxidized with rates controlled by linear diffusion, electrogenerated Fc(+) radicals diffusing outwards from the electrode react quantitatively (K23 C=5.8 x 10(exp 8) with dMeFc diffusing towards the electrode and produce Fc and dMeFc. That reaction replaces dMeFc with Fc, whose diffusion coefficient is higher than that of dMeFc(+), and the total mass-transfer limited current from the mixture is increased by approximately 10%. Analogous observations are made when mass-transfer is controlled by convective-diffusion as in RDE voltammetry. Similar results have been obtained with another, and for all practical purposes randomly selected pair of redox-active substances, [Co(bipy)3](2+) and N - methylphenothiazine (MePTZ); reaction of MePTZ(+) with [Co(bipy)3](2+) replaces the latter with MePTZ, which diffuses faster and the current increases by approximately 20%. The experimental voltammograms have been simulated numerically and the role of (a) the rate constant of the homogeneous reaction; (b) the relative concentrations; and, (c) the diffusion coefficients of all species involved have been studied in detail. Importantly, it was also identified that within any given redox system the dependence of the mass-transfer limited current on the bulk concentrations of the redox-active species is expected to be non-linear. These findings are discussed in terms of their electroanalytical implications.
A survey of 15 lines of sight for the CN B2Sigma(+) X2Sigma(+) interstellar absorption lines shows that the CN column density in diffuse interstellar clouds follows the relation log N(CN) is proportional to m log N(H2), where m is approximately equal to 3. This result is reproduced by a reaction network in which CN is produced primarily from C2 by the neutral-neutral reaction C2 + N yields CN + C, and photodissociation is the main destruction pathway for the neutral molecules CH, C2, and CN. The CN radical is the first molecular species observed in diffuse clouds that requires a neutral-neutral reaction for its formation in the gas phase. The network also reproduces the observed ratio N(CN)/N(H2).
The observational constraints on the composition of the interstellar dust are analyzed. The dust in the diffuse interstellar medium consists of a mixture of stardust (amorphous silicates, amorphous carbon, polycyclic aromatic hydrocarbons, and graphite) and interstellar medium dust (organic refractory material). Stardust seems to dominate in the local diffuse interstellar medium. Inside molecular clouds, however, icy grain mantles are also important. The structural differences between crystalline and amorphous materials, which lead to differences in the optical properties, are discussed. The astrophysical consequences are briefly examined. The physical principles of grain surface chemistry are discussed and applied to the formation of molecular hydrogen and icy grain mantles inside dense molecular clouds. Transformation of these icy grain mantles into the organic refractory dust component observed in the diffuse interstellar medium requires ultraviolet sources inside molecular clouds as well as radical diffusion promoted by transient heating of the mantle. The latter process also returns a considerable fraction of the molecules in the grain mantle to the gas phase.
We have combined Selected Area Electron Diffraction (SAED) and cryogenic techniques in an instrumental configuration that allows observing the structure of vapor deposited ice as it evolves during warmup. The ice is deposited in-situ inside an Hitachi H-500 H transmission electron microscope at a base pressure of 1-5 x 10(exp -7) torr on a thin amorphous carbon substrate at 15K or 86K and warmed up at a rate of 1-2 K/min. We find a progression of amorphous forms and well defined amorphous to amorphous transitions. Apart from the well known low-density form of ice, we confirm the presence of a high-density form and find a third amorphous form that coexists with cubic ice. We will report too on the amorphous to crystalline transition and the implications of these results for radical diffusion and gas retention observed in laboratory analog studies of interstellar and cometary ices.
An ab initio study of the low-lying electronic states of C4H has been presented where the species studied has a chi(2)sigma(+) ground state and two low lying pi states. Based on the vertical and adiabatic excitation energies between those states it is suggested that the 4428 A diffuse interstellar band is not carried by C4H. The application of the particle in a box model shows strong coincidences between the strong DIB's and predicted wavelengths of pi-pi transitions in C(2n)H series. Based on those coincidences, it is suggested the C(2n)H species as good candidates for carriers of diffuse interstellar bands.
Diffusion flame stabilization is of essential importance in both Earth-bound combustion systems and spacecraft fire safety. Local extinction, re-ignition, and propagation processes may occur as a result of interactions between the flame zone and vortices or fire-extinguishing agents. By using a computational fluid dynamics code with a detailed chemistry model for methane combustion, the authors have revealed the chemical kinetic structure of the stabilizing region of both jet and flat-plate diffusion flames, predicted the flame stability limit, and proposed diffusion flame attachment and detachment mechanisms in normal and microgravity. Because of the unique geometry of the edge of diffusion flames, radical back-diffusion against the oxygen-rich entrainment dramatically enhanced chain reactions, thus forming a peak reactivity spot, i.e., reaction kernel, responsible for flame holding. The new results have been obtained for the edge diffusion flame propagation and attached flame structure using various C1-C3 hydrocarbons.
The structure, ignition, and stabilization mechanisms for a methanol (CH3OH)-air triple flame are studied using Direct Numerical Simulations (DNS). The methanol (CH3OH)-air triple flame is found to burn with an asymmetric shape due to the different chemical and transport processes characterizing the mixture. The excess fuel, methanol (CH3OH), on the rich premixed flame branch is replaced by more stable fuels CO and H2, which burn at the diffusion flame. On the lean premixed flame side, a higher concentration of O2 leaks through to the diffusion flame. The general structure of the triple point features the contribution of both differential diffusion of radicals and heat. A mixture fraction-temperature phase plane description of the triple flame structure is proposed to highlight some interesting features in partially premixed combustion. The effects of differential diffusion at the triple point add to the contribution of hydrodynamic effects in the stabilization of the triple flame. Differential diffusion effects are measured using two methods: a direct computation using diffusion velocities and an indirect computation based on the difference between the normalized mixture fractions of C and H. The mixture fraction approach does not clearly identify the effects of differential diffusion, in particular at the curved triple point, because of ambiguities in the contribution of carbon and hydrogen atoms' carrying species.
The hydroxyl radical is known to be one of the most important intermediate species in the combustion processes. The hydroxyl radical has also been considered a dominant oxidizer of soot particles in flames. In this investigation the hydroxyl concentration profiles in sooting diffusion flames were measured by the laser-saturated fluorescence (LSF) method. The temperature distributions in the flames were measured by the two-line LSF technique and by thermocouple. In the sooting region the OH fluorescence was too weak to make accurate temperature measurements. The hydroxyl fluorescence profiles for all four flames presented herein show that the OH fluorescence intensities peaked near the flame front. The OH fluorescence intensity dropped sharply toward the dark region of the flame and continued declining to the sooting region. The OH fluorescence profiles also indicate that the OH fluorescence decreased with increasing height in the flames for all flames investigated. Varying the oxidizer composition resulted in a corresponding variation in the maximum OH concentration and the flame temperature. Furthermore, it appears that the maximum OH concentration for each flame increased with increasing flame temperature.
A computer model which simulates, in principle, the chemical changes in the photooxidation of hydrocarbons using as input data a set of elementary reactions, corresponding kinetic rate data and appropriate initial conditions was developed. The Model was refined and exploited to examine more closely the photooxidation and photostabilization of a hydrocarbon polymer. The results lead to the following observations. (1) The time to failure, tau sub f (chosen as the level of 5% C-H bond oxidation which is within the range anticipated for marked change in mechanical properties) varies as the inverse square root of the light intensity. However, tau sub f is almost unaffected by both the photoinitiator type and concentration. (2) The time to failure decreases with the rate of abstraction of C-H by peroxy radicals but increases with the rate of bimolecular radical termination controlled by diffusion. (3) Of the various stabilization mechanisms considered, the trapping of peroxy radicals is distinctly the most effective, although the concommitant decomposition of hydroperoxide is also desirable.
The goal of the present study is to assess numerically the ability of single-step and two-step chemical models to describe the main features encountered during the interaction between a two-dimensional vortex pair and a premixed laminar flame. In the two-step mechanism, the reaction kinetics are represented by a first chain branching reaction A + X yields 2X and a second chain termination reaction X + X yields P. This paper presents the fundamental mechanisms occurring during vortex-flame interactions and the relative impact of the major parameters encountered in turbulent premixed flames and suspected of playing a role in quenching mechanism: (1) Influence of stretch is investigated by analyzing the contribution of curvature and tangential strain on the local structure of the flame. The effect of Lewis number on the flame response to a strained field is analyzed. (2) Radiative heat losses which are suspected to be partially or totally responsible for quenching are also investigated. (3) The effect of the diffusion of the radicals is studied using a two-step mechanism in which an intermediate species is present. The parameters of the two-step mechanism are entirely determined from physical arguments. (4) Precise quantitative comparisons between the DNS and the experimental results of Samaniego et al are performed. These comparisons concern the evolution of the minimum heat release rate found along the flame front during the interaction and the distribution of the heat release rate along the flame front.
In this paper we present the results of a series of numerical experiments done on the propagation and initiation of a detonation wave. The calculations are performed in one-dimension, with considerable grid resolution. Of particular interest are the following questions: (1) the nature of periodic and chaotic instabilities generated by the wave; (2) the influence of the grid resolution on these instabilities; (3) the influence of the 'quality' of the numerical scheme; and (4) the influence of 'noise'. In the calculations, we use a second-order Total Variation Diminishing (TVD) scheme as the basic numerical method, with grid spacings as low as a fraction of a micron. Detonations waves are generated at the closed end of a tube, and allowed to propagate for approximately 20 cm. The required energy for successful initiation of the detonation will be measured for different cases of grid resolution and numerical schemes. A modified version of the TVD scheme has also been devised, which allows for much lower numerical diffusion of the radical species in the exponentially growing region behind the shock. The effect of this modification will be demonstrated. Oscillations in peak pressure and induction length are seen to develop in some cases: the oscillations can go through a sequence of modes, from a regular, high frequency mode to a low frequency mode with period doubling. A chaotic regime can also be obtained. General conclusions on the quality of algorithms will be presented. We will also discuss the performance of a version of the code developed on the IBM SP2 parallel computer.
Diffusion flames are commonly used for industrial burners in furnaces and flares. Oxygen/fuel burners are usually diffusion burners, primarily for safety reasons, to prevent flashback and explosion in a potentially dangerous system. Furthermore, in most fires, condensed materials pyrolyze, vaporize, and burn in air as diffusion flames. As a result of the interaction of a diffusion flame with burner or condensed-fuel surfaces, a quenched space is formed, thus leaving a diffusion flame edge, which plays an important role in flame holding in combustion systems and fire spread through condensed fuels. Despite a long history of jet diffusion flame studies, lifting/blowoff mechanisms have not yet been fully understood, compared to those of premixed flames. In this study, the structure and stability of diffusion flames of gaseous hydrocarbon fuels in coflowing air at normal earth gravity have been investigated experimentally and computationally. Measurements of the critical mean jet velocity (U(sub jc)) of methane, ethane, or propane at lifting or blowoff were made as a function of the coflowing air velocity (U(sub a)) using a tube burner (i.d.: 2.87 mm) (Fig. 1, left). By using a computational fluid dynamics code with 33 species and 112 elementary reaction steps, the internal chemical-kinetic structures of the stabilizing region of methane and propane flames were investigated (Fig. 1, right). A peak reactivity spot, i.e., reaction kernel, is formed in the flame stabilizing region due to back-diffusion of heat and radical species against an oxygen-rich incoming flow, thus holding the trailing diffusion flame. The simulated flame base moved downstream under flow conditions close to the measured stability limit.
Diffusion flames are commonly used for industrial burners in furnaces and flares. Oxygen/fuel burners are usually diffusion burners, primarily for safety reasons, to prevent flashback and explosion in a potentially dangerous system. Furthermore, in most fires, condensed materials pyrolyze, vaporize, and burn in air as diffusion flames. As a result of the interaction of a diffusion flame with burner or condensed-fuel surfaces, a quenched space is formed, thus leaving a diffusion flame edge, which plays an important role in flame holding in combustion systems and fire spread through condensed fuels. Despite a long history of jet diffusion flame studies, lifting/blowoff mechanisms have not yet been fully understood, compared to those of premixed flames. In this study, the structure and stability of diffusion flames of gaseous hydrocarbon fuels in coflowing air at normal earth gravity have been investigated experimentally and computationally. Measurements of the critical mean jet velocity (U(sub jc)) of methane, ethane, or propane at lifting or blowoff were made as a function of the coflowing air velocity (U(sub a)) using a tube burner (i.d.: 2.87 mm). By using a computational fluid dynamics code with 33 species and 112 elementary reaction steps, the internal chemical-kinetic structures of the stabilizing region of methane and propane flames were investigated. A peak reactivity spot, i.e., reaction kernel, is formed in the flame stabilizing region due to back-diffusion of heat and radical species against an oxygen-rich incoming flow, thus holding the trailing diffusion flame. The simulated flame base moved downstream under flow conditions close to the measured stability limit.
Impurity-helium solids are formed when a mixture of impurity and helium gases enters a volume of superfluid helium. Typical choices of impurity gas are hydrogen deuteride, deuterium, nitrogen, neon and argon, or a mixture of these. These solids consist of individual impurity atoms and molecules as well as clusters of impurity atoms and molecules covered with layers of solidified helium. The clusters have an imperfect crystalline structure and diameters ranging up to 90 angstroms, depending somewhat on the choice of impurity. Immediately following formation the clusters aggregate into loosely connected porous solids that are submerged in and completely permeated by the liquid helium. Im-He solids are extremely effective at stabilizing high concentrations of free radicals, which can be introduced by applying a high power RF dis- charge to the impurity gas mixture just before it strikes the super fluid helium. Average concentrations of 10(exp 19) nitrogen atoms/cc and 5 x 10(exp 18) deuterium atoms/cc can be achieved this way. It shows a typical sample formed from a mixture of atomic and molecular hydrogen and deuterium. It shows typical sample formed from atomic and molecular nitrogen. Much of the stability of Im-He solids is attributed to their very large surface area to volume ratio and their permeation by super fluid helium. Heat resulting from a chance meeting and recombination of free radicals is quickly dissipated by the super fluid helium instead of thermally promoting the diffusion of other nearby free radicals.
Diffusion flame stabilization in normal earth gravity (1 g) has long been a fundamental research subject in combustion. Local flame-flow phenomena, including heat and species transport and chemical reactions, around the flame base in the vicinity of condensed surfaces control flame stabilization and fire spreading processes. Therefore, gravity plays an important role in the subject topic because buoyancy induces flow in the flame zone, thus increasing the convective (and diffusive) oxygen transport into the flame zone and, in turn, reaction rates. Recent computations show that a peak reactivity (heat-release or oxygen-consumption rate) spot, or reaction kernel, is formed in the flame base by back-diffusion and reactions of radical species in the incoming oxygen-abundant flow at relatively low temperatures (about 1550 K). Quasi-linear correlations were found between the peak heat-release or oxygen-consumption rate and the velocity at the reaction kernel for cases including both jet and flat-plate diffusion flames in airflow. The reaction kernel provides a stationary ignition source to incoming reactants, sustains combustion, and thus stabilizes the trailing diffusion flame. In a quiescent microgravity environment, no buoyancy-induced flow exits and thus purely diffusive transport controls the reaction rates. Flame stabilization mechanisms in such purely diffusion-controlled regime remain largely unstudied. Therefore, it will be a rigorous test for the reaction kernel correlation if it can be extended toward zero velocity conditions in the purely diffusion-controlled regime. The objectives of this study are to reveal the structure of the flame-stabilizing region of a two-dimensional (2D) laminar jet diffusion flame in microgravity and develop a unified diffusion flame stabilization mechanism. This paper reports the recent progress in the computation and experiment performed in microgravity.
Reconsideration of inner radiation belt processes and proton theories. The ability of trapped particles to diffuse across magnetic shells is shown radically to alter the view of magnetospheric processes and to call for such a reconsideration. Inner belt calculations incorporating radial diffusion processes are performed in a manner parametrizing the radial diffusion coefficient. The latter is then varied so as to obtain agreement with experiment. The diffusion coefficient thus obtained is compared with coefficients derived by other means.