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.
Approximate formulas for calculating thermodynamic characteristics and equilibrium constants of chemical reactions
Results are presented from a study of the equilibria, kinetic rates, and the aggregation pathway which leads from a lysozyme monomer crystal to a tetragonal crystal, using dialyzed and recrystallized commercial hen eggwhite lysozyme. Relative light scattering intensity measurements were used to estimate the initial equilibrium constants for undersaturated lysozyme solutions in the tetragonal regime. The K1 value was estimated to be (1-3) x 10 exp 4 L/mol. Estimates of subsequent equilibrium constants depend on the crystal aggregation model chosen or determined. Experimental data suggest that tetragonal lysozyme crystal grows by addition of aggregates preformed in the bulk solution, rather than by monomer addition.
It is shown here that the equilibrium constant for the reversible reaction ClO + O2 - ClO-O2 at stratospheric temperatures must be at least three orders of magnitude less than the current NASA upper limit. The new upper limit greatly diminishes the possible role of ClO-O2 in the chlorine-photosensitized decomposition of O3. Nevertheless, it does not preclude the possibility that it is a significant reservoir of ClO, as well as a possible reactant, at low temperatures characteristic of polar vortices.
Direct determinations of the equilibrium constant for the reaction N2O5 = NO2 + NO3 were carried out by measuring NO2, NO3, and N2O5 using long-path visible and infrared absorption spectroscopy as a function of temperature from 243 to 397 K. The first-order decay rate constant of N2O5 was experimentally measured as a function of temperature. These results are in turn used to derive a value for the rate coefficient for the NO-forming channel in the reaction of NO3 with NO2. The implications of the results for atmospheric chemistry, the thermodynamics of NO3, and for laboratory kinetics studies are discussed.
The computer codes developed provide data to 30000 K for the thermodynamic and transport properties of individual species and reaction rates for the prominent reactions occurring in an 11-species nonequilibrium air model. These properties and the reaction-rate data are computed through the use of curve-fit relations which are functions of temperature (and number density for the equilibrium constant). The curve fits were made using the most accurate data believed available. A detailed review and discussion of the sources and accuracy of the curve-fitted data used herein are given in NASA RP 1232.
A basic, experimental study on the hydrochlorination of silicon tetrachloride and metallurgical grade silicon with hydrogen gas to form trichlorosilane was carried out to greatly expand the range of reaction conditions. The equilibrium constant, K sub p, for the hydrochlorination reaction was measured as a function of temperature, pressure and concentration. The variation of the equilibrium constant as a function of temperature provided the measurement on the heat of reaction, delta H, by the Second Law Method. The value of delta H was measured to give 10.6 Kcal/mole. The equilibrium constant was also studied as a function of concentration. In agreement with the theory, the equilibrium constant remained constant with respect to the varying H2/SiCl4 feed ratios. On the other hand, the effect of pressure on the equilibrium constant was found to be more complex.
Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.
Purpose of research on rough-to-smooth transition of an equilibrium neutral constant stress layer is to develop a model for low-level atmospheric flow over terrains of abruptly changing roughness, such as those occurring near the windward end of a landing strip, and to use the model to derive functions which define the extent of the region affected by the roughness change and allow adequate prediction of wind and shear stress profiles at all points within the region. A model consisting of two bounding logarithmic layers and an intermediate velocity defect layer is assumed, and dimensionless velocity and stress distribution functions which meet all boundary and matching conditions are hypothesized. The functions are used in an asymptotic form of the equation of motion to derive a relation which governs the growth of the internal boundary layer. The growth relation is used to predict variation of surface shear stress.
The entropy of a gas system with the number of particles subject to external control is maximized to derive relations between the thermodynamic variables that obtain at equilibrium. These relations are described in terms of the chemical potential, defined as equivalent partial derivatives of entropy, energy, enthalpy, free energy, or free enthalpy. At equilibrium, the change in total chemical potential must vanish. This fact is used to derive the equilibrium constants for chemical reactions in terms of the partition functions of the species involved in the reaction. Thus the equilibrium constants can be determined accurately, just as other thermodynamic properties, from a knowledge of the energy levels and degeneracies for the gas species involved. These equilibrium constants permit one to calculate the equilibrium concentrations or partial pressures of chemically reacting species that occur in gas mixtures at any given condition of pressure and temperature or volume and temperature.
Light-scattering intensity measurements are a sensitive method for following changes in the hydrodynamic radius of particles in solution. The approach used in this report utilizes the light-scattering intensity ratios of a polydisperse to a monodisperse system. By numerically modeling the process, and fitting the model curves to the data, estimates have been obtained for the dimerization equilibrium constant, the dimer + dimer yields tetramer equilibrium constant, and the association rate constant for the dimerization process.
Hydrogen bonding in deuterochloroform studied using IR spectra and dielectric constant methods, calculating equilibrium constants and heat of formation
To support the development of Controlled Ecological Life Support Systems (CELSS) in the space program, a metabolic simulator has been selected for use in a closed chamber to test functions of the CELSS. This metabolic simulator is a catalytic reactor which oxidizes the methyl acetate to produce carbon dioxide and water vapor. In this project, kinetic studies of catalytic oxidation of methyl acetate were conducted using monolithic and pellet catalysts with 0.5% (by weight) platinum (Pt) on aluminum oxide (Al2O3). The reaction was studied at a pressure of one atmosphere and at temperatures varying from 160 C to 420 C. By-products were identified at the exit of the preheater and reactor. For the kinetic study with the monolithic catalyst, a linear regression method was used to correlate the kinetic data with zero-order, first-order and Langmuir-Hinshelwood models. Results indicate that the first-order model represents the data adequately at low concentrations of methyl acetate. For higher concentrations of methyl acetate, the Langmuir-Hinshelwood model best represents the kinetic data. Both rate constant and adsorption equilibrium constants were estimated from the regression. A Taguchi orthogonal array (L(sub 9)) was used to investigate the effects of temperature, flow rate, and concentration on the catalytic oxidation of methyl acetate. For the monolithic catalyst, temperature exerts the most significant effect, followed by concentration of methyl acetate. For the pellet catalyst, reaction temperature is the most significant factor, followed by gas flow rate and methyl acetate concentration. Concentrations of either carbon dioxide or oxygen were seen to have insignificant effect on the methyl acetate conversion process. Experimental results indicate that the preheater with glass beads can accomplish thermal cracking and catalytic reaction of methyl acetate to produce acetic acid, methanol, methyl formate, and 1-propanol. The concentration of all by-products was measured in ppmv (parts per million by volume). At higher temperatures, greater amounts of these products are produced, as expected. In all cases, methanol was the predominant concentration detected, followed by methyl formate. At temperatures lower than 320 C for the P-type monolithic catalyst, methanol, acetic acid, and acetone were detected, whereas, for the E-type monolithic catalyst, only methanol was detected at 160 C. Both P and E types of the monolithic catalyst were specified with the same substrates (ceramic), washcoat (Al2O3), and promoter (Pt). However, the manufacturing and treatment procedures were quite different. It was therefore concluded that the performance of the E-type monolithic catalyst is superior to that of the P-type for oxidation of methyl acetate. At higher reaction temperatures, e.g., above 420 C, all reactants and byproducts were completely oxidized using these two types of monolithic catalyst to produce carbon dioxide and water vapor. A complex heterogenous catalytic reaction mechanism was proposed to explain the formation of the byproducts (methanol, acetic acid, and methyl formate) as the methyl acetate traveled through the preheater packed with glass beads. The by-product, 1-propanol, may be formed only through a homogeneous reaction, since it is difficult to develop a reasonable sequence of heterogeneous reaction steps to explain its formation. The homogeneous thermal decomposition of methyl acetate to form free radicals was proposed to explain the formation of 1-propanol, and also methanol, in the preheater. A dual-site catalytic reaction mechanism was proposed for the oxidation of methyl acetate over Pt/Al2O3 monolithic catalyst. The dual-site mechanism describes the chemisorption of oxygen molecules as well as a physical adsorption of methyl acetate on the active sites. On the active sites, methyl acetate is oxidized rapidly to form carbon dioxide and water vapor. A rate equation derived from this mechanism gives the Langmuir-Hinshelwood rate formula which has been observed from the experimental data obtained in this project for high methyl acetate concentration (greater than 1000 ppmv) over a monolithic catalyst. If the oxygen concentration is very high and methyl acetate concentration is very low, the reaction rate equation is then reduced to a first-order with respect to methyl acetate concentration. The first-order model has also been observed from the experimental data obtained in this project for low methyl acetate concentration (less than 1000 ppmv).
The hydrochlorination of silicon tetrachloride with hydrogen and metallurgical grade (m.g.) silicon metal, 3 SiCl4 + 2 H2 + Si yields 4 SiHCl3 was shown to be an efficient process to produce trichlorosilane. A research and development program was carried out to study the hydrochlorination reaction over a wide range of reaction conditions. Equilibrium constant and reaction kinetics measurements were made to provide the basis for a theoretical study on the hydrochlorination process. Thermodynamic properties of the hydrochlorination reaction were also measured. The effects of temperature, pressure, and concentration on the equilibrium constant, K sub p, were studied.
The mole fraction x(sub 2) of (H2)2 in equilibrium mixture with H2 under the atmospheric conditions of Jupiter is evaluated from the dimerization equilibrium constant calculated by quantum-chemical treatments and also from the Lennard-Jones potential. The treatments are of an ab initio type with the second and fourth order Moller-Plesset perturbation techniques and a basis set superposition error evaluation. The computed dimerization equilibrium constant is combined with observed height profiles of temperature and pressure. In six treatments considered it is found that the mole fraction decreases with increasing height. Various approximations suggest the dimeric mole fraction at the Jupiter 1 atm pressure level between 0.04 and 1.06%.
Published data are used to calculate equilibrium constants for reactions of the formation of nitric acid hydrates over the temperature range 190 to 205 K. Standard enthalpies of formation and standard entropies are calculated for the tri- and mono-hydrates. These are shown to be in reasonable agreement with earlier calorimetric measurements. The formation of nitric acid trihydrate in the polar stratosphere is discussed in terms of these equilibrium constants.
Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].
Understanding and predicting the effect of epoxide structure on the rate of polymerization in epoxide/CO 2 ring opening copolymerization catalysis is a long-standing challenge. Here, a known highly active Co(III)K(I) catalyst is used to investigate the influences of six different epoxides' binding strengths on their rates of copolymerization. Since calculations and experiments indicate that studying the catalytically relevant Co(III)−epoxide adduct directly is experimentally challenging, epoxide−catalyst binding interactions are quantified using a Co(II)K(I) complex to model the key catalytic intermediate. Epoxide−catalyst coordination is investigated using UV−vis spectroscopy titrations which provide fast and effective determination of association or binding constants. The epoxide−catalyst equilibrium constants show a clear exponential correlation with copolymerization rates and a new catalyst performance linear free energy relationship is revealed. Epoxides exhibiting stronger catalyst binding constants show higher copolymerization rates. The structure−activity correlation is consistent with the polymerization kinetics, mechanism and DFT calculations. Both the methods to investigate epoxide−catalyst coordination and the linear free energy relationship are shown to apply to the series of six epoxides and a second Co(III)K(I) catalyst. These structure− performance relationships are likely applicable to other transition metal catalysts and should expedite future epoxide and catalyst selection to make useful poly(carbonate) materials.
A model is proposed for low-level atmospheric flows over terrains of changing roughness length, such as those found at the windward end of landing strips adjoining rough terrain. The proposed model is used to develop a prediction technique for calculating transition wind and shear-stress profiles in the region following surface roughness discontinuity. The model for the transition region comprises two layers: a logarithmic layer and a buffer layer. The flow is assumed to be steady, two-dimensional, and incompressible, with neutral hydrostatic stability. A diagram is presented for a typical wind profile in the transition region, obtained from the logarithmic and velocity defect profiles using shear stress calculated by relevant equations.