SEARCH · Engineering Papers
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.
Dissociational rate constants and vibrational relaxation times modifications, considering coupling of vibrational and dissociational nonequilibrium of expanding gas in nozzle
Recent work in the aerothermodynamics community has focused on the development of reduced-order models for thermo-chemical non-equilibrium which avoid the restrictive assumptions of multi-temperature models and the prohibitive cost associated with State-to-State (StS) models. In the present work, this is accomplished by lumping energy states together and assuming groups of states are roughly in equilibrium. As a result, the non-equilibrium behavior of a gas can be captured at a reduced computational cost from a full StS simulation. In this work, we present a hybrid grouping model for studying energy transfer and dissociation in a mixture of nitrogen molecules due to N2-N2 reactions. This is accomplished by making use of a grouping strategy informed by data from the N2-N StS kinetic data. However, due to the massive computational cost associated with constructing StS data for the N2-N2 system, the kinetic data for the hybrid grouping model are calculated using the quasi-classical trajectory (QCT) method by sampling states for trajectory within the groups. This general framework is called the Maximum-Entropy Quasi-Classical Trajectory (ME-QCT) method. The primary challenge associated with this method is that rates for reverse grouped reactions cannot be obtained through detailed balance at a group level, due to the variation of group internal temperatures. To construct the full model for N2-N2 grouped kinetics using the ME-QCT method, detailed balance is invoked at the microscopic level, allowing for the calculation of the full kinetic data from QCT. Results will be presented using the full ME-QCT model for the N2-N2 system in an isothermal and isochoric reactor simulation. In addition, simple CFD test cases for a one-dimensional standing shock and for a quasi-one-dimensional nozzle will be used for demonstration of the ME-QCT method. This method allows for the calculation of non-equilibrium behavior for the N2-N2 system without the prohibitive cost of a full StS simulation. Moreover, it enables the construction of a unified model for the dissociating and recombining non-equilibrium flows.
This work presents a general framework for model reduction of non-equilibrium energy transfer and dissociation processes. The multi-group maximum-entropy method is coupled with the quasi-classical trajectory method to directly construct a reduced order model for chemical non-equilibrium. Kinetic data is calculated by applying detailed balance at a microscopic level, overcoming the limitation of calculating recombination kinetic data. This approach enables the construction of a reduced order model for kinetics which bypasses the need to compute state-to-state kinetic data. This physics based reduced order model ensures that an equilibrium distribution is reached given infinite time, while allowing for non-equilibrium distributions during the relaxation and dissociation processes. A proof-of-concept test case demonstrates the applicability of this model by comparison with state-to-state kinetic data for the N2(Chi (sup 1) Sigma ((sub g), (sup plus sign)) minus N ((exp 4) S (sub u)) system with excellent agreement.
Comparisons are made between potential energy surfaces (PES) for N2 + N and N2 + N2 collisions and between rate coefficients for N2 dissociation that were computed using the quasiclassical trajectory method (QCT) on these PESs. For N2 + N we compare the Laganas empirical LEPS surface with one from NASA Ames Research Center based on ab initio quantum chemistry calculations. For N2 + N2 we compare two ab initio PESs (from NASA Ames and from the University of Minnesota). These use different methods for computing the ground state electronic energy for N4, but give similar results. Thermal N2 dissociation rate coefficients, for the 10,000K-30,000K temperature range, have been computed using each PES and the results are in excellent agreement. Quasi-stationary state (QSS) rate coefficients using both PESs have been computed at these temperatures using the Direct Molecular Simulation of Schwartentruber and coworkers. The QSS rate coefficients are up to a factor of 5 lower than the thermal ones and the thermal and QSS values bracket the results of shock-tube experiments. We conclude that the combination of ab initio quantum chemistry PESs and QCT calculations provides an attractive approach for the determination of accurate high-temperature rate coefficients for use in aerothermodynamics modeling.
Grain boundaries in FCC metals with low stacking-fault energy can form in dissociated configurations of stacking faults. Perhaps the most studied example is 9R stacking at boundaries near Σ3{112}, where the distribution of stacking faults is related to the emission of Shockley partial dislocations. Here, we combine atomic-resolution electron microscopy, atomistic simulations, and dislocation theory to demonstrate that boundaries vicinal to Σ33a support the stabilization of dissociated 9R stacking within a narrow range of inclinations. This boundary is interesting since its misorientation (20.05°) lies in the medium-angle regime, just past the upper misorientation limit for low-angle boundaries, motivating questions for how best to describe it in terms of dislocations. Our HAADF-STEM observations of thin film bicrystals, supported by atomistic modeling, reveal that this inclination dependence arises from specific geometric constraints on the arrangement of Shockley partial dislocations at the interface. Quantification of stacking-fault distributions across multiple boundaries indicates that the density and spacing of faults closely follow the ideal 9R motif, with subtle variations reflecting the complex energy landscape of these boundaries. Through energetic analysis, we establish the presence of competing metastable states enabled by variations in stacking sequences, emphasizing the significant role of crystallographic geometry. We generalize our analysis as a function of misorientation, showing how 9R at the Σ33a boundary is related to previous observations and calculations of HCP at a 29.7° boundary. This study provides a crystallographically grounded framework connecting dislocation structures, stacking-fault distributions, and metastability at grain boundaries in FCC metals.
Knowledge of how trace amounts of more reactive metals influence the oxidation rate and mechanism of Cu surfaces is essential for developing strategies to optimize the performance of Cu-based catalysts. We find that the addition of 1% Rh to Cu(111) increases the initial O 2 dissociation rate by approximately 9-fold. CO poisoning experiments reveal that single Rh atoms activate O 2 and facilitate the spillover of atomic oxygen to Cu sites. Scanning tunneling microscopy (STM) and in situ X-ray photoelectron spectroscopy (XPS) support this mechanism, showing enhanced surface oxygen near Rh atoms. Here, a density functional theory (DFT)-based model demonstrates that Rh binds the O 2 precursor 0.15 eV more strongly than Cu(111) and lowers the O 2 dissociation barrier by 0.02 eV. Both single-crystal and nanoparticle experiments show that at low oxygen pressures, Rh enhances Cu oxidation, whereas at higher pressures, it inhibits deeper oxidation, as evidenced by in situ ultraviolet-visible (UV-vis) spectra.
Bipolar membranes (BPMs) create counteracting spatial gradients of pH and electrostatic potential in electrochemical systems, enabling applications in pH regulation, electrocatalysis, and separations. At the polarized junction of a BPM the water dissociation (WD, 2H2O ⇌ H3O+ + OH-) reaction can be driven, but it remains poorly understood. In this Perspective, we integrate molecular insights from bulk-water autoionization and the associated field effects with continuum descriptions of BPM electrostatics and experimental WD kinetic analyses to describe possible mechanisms of voltage-driven WD. Pristine BPM junctions highlight both the limits of primarily electric-field-driven WD and the practical challenges of junction stability at extreme reverse bias. Introducing heterogeneous catalyst layers, commonly metal oxides and graphene oxides, accelerates WD by orders of magnitude through hypothesized coupled effects in which surface acid-base functionality and high-density hydroxyl sites mediate proton-transfer steps, and catalyst mobile electronic/ionic charges redistribute the junction electric potential drop to shape the local electric fields and reactive microenvironments. Kinetic analyses suggest two regimes of heterogeneous WD mechanism, including field-driven ordering of interfacial water and a Second-Wien-Effect dissociation-barrier lowering. We conclude by defining the key unknown variables (local pH, electrostatic potential, catalyst charge state and relationships among mechanisms) and outlining experimental and multiscale modeling strategies needed for predictive WD catalysis and for controlling related ion-transfer reactions.
Modeling chemical reactions with quantum chemical methods is challenging when the electronic structure varies significantly throughout the reaction and when electronic excited states are involved. Multireference methods, such as complete active space self-consistent field (CASSCF), can handle these multiconfigurational situations. However, even if the size of the needed active space is affordable, in many cases, the active space does not change consistently from reactant to product, causing discontinuities in the potential energy surface. The localized active space SCF (LASSCF) is a cheaper alternative to CASSCF for strongly correlated systems with weakly correlated fragments. The method is used for the first time to study a chemical reaction, namely the bond dissociation of a mono-, di-, and triphenylsulfonium cation. LASSCF calculations generate smooth potential energy scans more easily than the corresponding, more computationally expensive CASSCF calculations while predicting similar bond dissociation energies. Furthermore, our calculations suggest a homolytic bond cleavage for di- and triphenylsulfonium and a heterolytic pathway for monophenylsulfonium.
Electrochemical reactivity is known to be dictated by the structure and composition of the electrocatalyst–electrolyte interface. Here, we show that optically generated electric fields at this interface can influence electrochemical reactivity insofar as to completely switch reaction selectivity. We study an electrocatalyst composed of gold–copper alloy nanoparticles known to be active toward the reduction of CO 2 to CO. However, under the action of highly localized electric fields generated by plasmonic excitation of the gold–copper alloy nanoparticles, water splitting becomes favored at the expense of CO 2 reduction. Real-time time-dependent density functional tight binding calculations indicate that optically generated electric fields promote transient-hole-transfer-driven dissociation of the O─H bond of water preferentially over transient-electron-driven dissociation of the C─O bond of CO 2 . These results highlight the potential of optically generated electric fields for modulating pathways, switching reactivity on/off, and even directing outcomes.
Values of the bond dissociation energies, ionization potentials, and electron affinities that were taken from the literature are presented in tables for some monatomic, diatomic, and polyatomic molecules which are found in many high-temperature chemical reactions including combustion reactions. Much of the information came from literature published after 1950 which either reported experimental and theoretical energy values or gave a review of previous literature on the subject. In some cases values for the bond dissociation energies were calculated from recently published heats of formation.
Integral Hellmann-Feynman formula for dissociation energies in H, Li and LiH molecules, discussing transition densities and dissociation energy
Ground state dissociation energies for diatomic halogens from vibrational spacings near dissociation limit
The dissociation of molecular hydrogen by helium was studied at 5500 deg K using the quasiclassical approximation. The temperature-dependent reaction cross sections were very dependent on the initial quantum states of the molecule. A rate constant was estimated for the dissociation reaction, and agreement between theory and experiment was satisfactory.
In recent studies of the formation of molecules in interstellar clouds, Solomon and Klemperer (1972) and Frisch (1972) were able to fit the observed ratios of CH+ and CH densities only with the assumption that the rate of dissociative recombination is extremely small. It is known that the rate for most diatomic molecules is relatively large, with the exception of H2+ and He2+ for which it is anomalously small. The present work explores the possibility that CH+ also has a low recombination rate. However, preliminary calculations yield a large rate of the dissociative recombination process, seemingly necessitating considerable revision of some recent models of interstellar clouds.
Investigation of the dissociation kinetics of undiluted carbon monoxide over the 5,600 to 12,000 K temperature range. Data are presented that have been obtained as time-resolved pressure measurements on the end wall of a shock tube and radiation emission of a C2 Swan system (0-0 band) behind incident shock waves. The decomposition of CO is complex and includes a chain with C2 as an intermediate species. The dissociation rate for the overall process has been found to be independent of the proportions of the collision partners M = CO, C, and O. The rate constant found is on the average about 10 times that previously measured with argon as the collision partner.
Measurements of night-time 6300 A airglow intensities at the Arecibo Observatory have been compared with dissociative recombination calculations based on electron densities derived from simultaneous incoherent backscatter measurements. The agreement indicates that the nightglow can be fully accounted for by dissociative recombination. The comparisons are examined to determine the importance of quenching, heavy ions, ionization above the F-layer peak, and the temperature parameter of the model atmosphere. Comparable fits between the observed and calculated intensities are found for several available model atmospheres. The least-squares fitting process, used to make the comparisons, produces comparable fits over a wide range of combinations of neutral densities and of reaction constants. Yet, the fitting places constraints upon the possible combinations; these constraints indicate that the latest laboratory chemical constants and densities extrapolated to a base altitude are mutually consistent.
Measurements were made of the rotational profiles of specific bands of the CO fourth-positive group (4PG). The CO 4PG bands were excited by electron impact dissociative excitation of CO2. The results are applicable to analysis of the Mariner observations of the CO 4PG in the dayglow of Mars. The results indicate that dissociative excitation of CO2 by electron impact leads to CO(A 1Pi) fragments with a rotational distribution that is highly nonthermal. The parent CO2 temperature was about 300 K in the experiment, while the fragment CO(A 1Pi) showed emission band profiles consistent with a rotational temperature greater than about 1500 K. Laboratory measurement of the reduced transmission of the hot bands by thermal CO appears to be the most direct way of determining the column density responsible for the CO(v',0) absorption of Mars.
Most physical applications of the swirling flow, defined as a vortex superimposed on an axial flow in the nozzle, involve high temperatures and the possibility of real gas effects. The generalized one-dimensional swirling flow in a converging-diverging nozzle is analyzed for equilibrium and frozen dissociation using the ideal dissociating gas model. Numerical results are provided to illustrate the major effects and to compare with results obtained for a perfect gas with constant ratio of specific heats. It is found that, even in the case of real gases, perfect gas calculations can give a good estimate of the reduction in mass flow due to swirl.