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.
Quasi-equilibrium prediction of rate of volatilization /erosion/ of solid tungsten by reaction with gaseous oxygen at high temperature and low pressure
An iterative in vitro selection procedure was used to isolate a new class of catalytic RNAs (ribozymes) from a large pool of random-sequence RNA molecules. These ribozymes ligate two RNA molecules that are aligned on a template by catalyzing the attack of a 3'-hydroxyl on an adjacent 5'-triphosphate - a reaction similar to that employed by the familiar protein enzymes that synthesize RNA. The corresponding uncatalyzed reaction also yields a 3',5'-phosphodiester bond. In vitro evolution of the population of new ribozymes led to improvement of the average ligation activity and the emergence of ribozymes with reaction rates 7 million times faster than the uncatalyzed reaction rate.
Rate coefficients for the reactions of C2H + HCN yields products and C2H + CH3CN yields products have been measured over the temperature range 262-360 K. These experiments represent an ongoing effort to accurately measure reaction rate coefficients of the ethynyl radical, C2H, relevant to planetary atmospheres such as those of Jupiter and Saturn and its satellite Titan. Laser photolysis of C2H2 is used to produce C2H, and transient infrared laser absorption is employed to measure the decay of C2H to obtain the subsequent reaction rates in a transverse flow cell. Rate constants for the reaction C2H + HCN yields products are found to increase significantly with increasing temperature and are measured to be (3.9-6.2) x 10(exp 13) cm(exp 3) molecules(exp -1) s(exp -1) over the temperature range of 297-360 K. The rate constants for the reaction C2H + CH3CN yields products are also found to increase substantially with increasing temperature and are measured to be (1.0-2.1) x 10(exp -12) cm(exp 3) molecules(exp -1) s(exp -1) over the temperature range of 262-360 K. For the reaction C2H + HCN yields products, ab initio calculations of transition state structures are used to infer that the major products form via an addition/elimination pathway. The measured rate constants for the reaction of C2H + HCN yields products are significantly smaller than values currently employed in photochemical models of Titan, which will affect the HC3N distribution.
Theoretical evolutionary sequences of models for stars of 15 and 30 solar masses were computed from the zero-age main sequence to the end of core helium burning. During the earliest stages of core helium depletion, the envelope rapidly expands into the red-supergiant configuration. At 15 solar mass, a blue loop on the H-R diagram ensues if the initial metals abundance, initial helium abundance, or C-12 + alpha particle reaction rate is sufficiently large, or if the 3-alpha reaction rate is sufficiently small. These quantities affect the opacity of the base of the outer convection zone, the mass of the core, and the thermal properties of the core. The blue loop occurs abruptly and fully developed when the critical value of any of these quantities is exceeded, and the effective temperature range and fraction of the lifetime of core helium burning during the slow phase of the blue loop vary surprisingly little. At 30 solar mass no blue loop occurs for any reasonable set of input parameters.
A numerical study for confined, axisymmetrical, turbulent diffusion flames is presented. Local mean gas properties are predicted by solving the appropriate conservation equations in the finite difference form with the corresponding boundary conditions. The k-epsilon two-equation turbulence model is employed to describe the turbulent nature of the flow. A two-step kinetic model is assumed to govern the reaction mechanism. The finite reaction rate is the smaller of an Arrhenius type of reaction rate and a modified version of eddy-breakup model. Reasonable agreement is observed between calculations and measurements, but to obtain better agreement, more work is needed on improvements of the above mathematical models. However, the present numerical study offers an improvement in the analysis and design of the gas turbine combustors.
An S-matrix formalism is developed which can be applied to reactions in which electron screening is important, including three-body reactions and reactions involving weak interactions. The various regimes of the triple-alpha reactions are systematically discussed, and a new nonresonant regime at high densities is identified. Using the S-matrix formalism, an analytic expression is obtained for the screened triple-alpha reactions which is accurate for all temperatures and densities. The results are compared with those of Cameron (1959) and Nomoto et al. (1985), and the latter's expression for the unscreened reaction rate is verified. However, it is shown that the reaction rate in the pycnonuclear regime cannot be obtained from the unscreened rate using a screening factor, and that the results of Nomoto et al. therefore cannot be used in this regime.
The paper presents a detailed theoretical analysis of the process of gas absorption to a thin liquid film adjacent to a horizontal rotating disk. The film is formed by the impingement of a controlled liquid jet at the center of the disk and subsequent radial spreading of liquid along the disk. The chemical reaction between the gas and the liquid film can be expressed as a zero-order homogeneous reaction. The process was modeled by establishing equations for the conservation of mass, momentum, and species concentration and solving them analytically. A scaling analysis was used to determine dominant transport processes. Appropriate boundary conditions were used to solve these equations to develop expressions for the local concentration of gas across the thickness of the film and distributions of film height, bulk concentration, and Sherwood number along the radius of the disk. The partial differential equation for species concentration was solved using the separation of variables technique along with the Duhamel's theorem and the final analytical solution was expressed using confluent hypergeometric functions. Tables for eigenvalues and eigenfunctions are presented for a number of reaction rate constants. A parametric study was performed using Reynolds number, Ekman number, and dimensionless reaction rate as parameters. At all radial locations, Sherwood number increased with Reynolds number (flow rate) as well as Ekman number (rate of rotation). The enhancement of mass transfer due to chemical reaction was found to be small when compared to the case of no reaction (pure absorption), but the enhancement factor was very significant when compared to pure absorption in a stagnant liquid film. The zero-order reaction processes considered in the present investigation included the absorption of oxygen in aqueous alkaline solutions of sodiumdithionite and rhodium complex catalyzed carbonylation of methanol. Present analytical results were compared to previous theoretical results for limiting conditions, and were found to have very good agreement.
Composition changes within the shock layer due to chemical reactions are considered. The Lighthill ideal dissociating gas model was used in an effort to describe the oxygen type molecule. First, the two limiting cases, when the chemical reaction rates are very slow and very fast in comparison to local convective rates, are investigated. Then, the problem is solved for arbitrary chemical reaction rates.
The role of divalent cations in the mechanism of pyrophosphate-activated, template-directed oligoribonucleotide ligation has been investigated. The dependence of the reaction rate on Mg2+ concentration suggests a kinetic scheme in which a Mg2+ ion must bind before ligation can proceed. Mn2+, Ca2+, Sr2+, and Ba2+ can also catalyze the reaction. Although Pb2+ and Zn2+ do not catalyze the reaction in the absence of other divalent ions, they significantly modulate the reaction rate when added in the presence of Mg2+, with Pb2+ stimulating the reaction (up to 65-fold) and Zn2+ inhibiting the reaction. The logarithm of the ligation rate increases linearly, with slope of 0.95, as a function of pH, indicating that the reaction involves a single critical deprotonation step. The ligation rates observed with the different divalent metal ion catalysts (Mn2+ > Mg2+ > Ca2+ > Sr2+ = Ba2+) vary inversely with the pKa values of their bound water molecules. The pH profile and these relative ligation rates suggest a mechanism in which a metal-bound hydroxide ion located near the ligation junction promotes catalysis, most likely by deprotonation of the hydroxl nucleophile. The effects of changing either the leaving group or the attacking hydroxyl, together with the large delta S(++) value for oligonucleotide ligation (about -20 eu), are consistent with an associative transition state.
The paper determines the physical properties - speed, width, and density structure - of conductive burning fronts in degenerate carbon-oxygen (C + O) and oxygen-neon-magnesium (O + Ne + Mg) compositions for a grid of initial densities and compositions. The dependence of the physical properties of the flame on the assumed values of nuclear reaction rates, the nuclear reaction network employed, the thermal conductivity, and the choice of coordinate system are investigated. The occurrence of accretion-induced collapse of a white dwarf is found to be critically dependent on the velocity of the nuclear conductive burning front and the growth rate of hydrodynamic instabilities. Treating the expanding area of the turbulent burning region as a fractal whose tile size is identical to the minimum unstable Rayleigh-Taylor wavelength, it is found, for all reasonable values of the fractal dimension, that for initial C + O or O + Ne + Mg densities above about 9 x 10 exp 9 g/cu cm the white dwarf should collapse to a neutron star.
A fire in a spacecraft poses detrimental consequences and risks mission success in addition to crew safety. This is compounded during long-duration missions when the crew has limited options to recover from a fire. A common spacecraft fire concern is the smoldering of wire insulation, typically made from Polyvinyl chloride (PVC) or Polytetrafluoroethylene (PTFE). This creates acid gases such as Hydrogen Chloride (HCl), Hydrogen Fluoride (HF) and Hydrogen Cyanide (HCN). These poisonous gases are hazardous to the crew. They also interact with common surfaces within the spacecraft more than dominant combustion products such as CO2 and H2O. This makes them more difficult to track for potential fire detection techniques, or for postfire clean-up. It is imperative to be able to understand and predict the fate of these poisonous species in a microgravity environment in order to design a safe vehicle. HCl interacts with a number of materials inside a spacecraft. Primary among these materials is aluminum, which is abundantly used due to its strong and light weight nature. Aluminum has a natural oxide layer that protects it from corrosion but is typically treated to enhance this oxide layer. Among these treatments is a chromate conversion coating (CCC), which provides a thin enough protective oxide layer to still conduct electricity, and a traditional anodized material that has a thicker oxide layer that does not conduct electricity. Nomex is another common material found inside a spacecraft. It is a flame-resistant woven polymer that is related to nylon. This commercially available material is used for cargo storage bags and as a fire barrier. Physics-based models were developed to predict the uptake of HCl by these materials. The ultimate objective of these models is to predict the fate of HCl within the spacecraft so that sensors can be placed in meaningful locations in future missions based on the model predictions. To support these modeling efforts, experiments were performed in a cast acrylic test cell that measured the difference between the inlet and outlet concentration of HCl after inserting a sample rod of the test material. Different uptake capacities were realized for each type of sample tested. A computational fluid dynamics model (CFD) model of the reactor was then constructed that used a one-step global reaction rate with calibratable reaction (or kinetic) constants. These constants were calibrated to match the HCl uptake on the CCC aluminum samples, and the same kinetic constants were then tested for the stock and anodized aluminum samples. Model predictions matched the experimental data for the stock aluminum, and to a much lesser extent, the anodized aluminum. The model was additionally validated at different flow rates, sample surface areas, and inlet concentrations, and showed good agreement for all stock and CCC samples. The model did not accurately predict the HCl uptake in the anodized samples compared to the other two types of aluminum. Adjusting the kinetic constants and transport properties did little to improve the prediction. X-Ray Photoelectron Spectroscopy (XPS) was used to determine that the oxide layer thickness of anodized aluminum is approximately 5,000 nm, compared to 250 nm for CCC and 50 nm for stock. XPS also revealed presence of chlorine further down in the aluminum oxide layer in anodized samples than CCC and stock samples after the samples were saturated with HCl, indicating that accounting for diffusion of HCl into the oxide layer is important for accurate prediction of HCl uptake onto anodized aluminum. Consequently, a multi-scale model was developed and tested. First, a single pore inside the anodized aluminum oxide layer was modeled and is referred to as the pore-scale model. In this model, HCl diffused through the pore and reacted with the aluminum oxide pore wall to create aluminum chloride. The sample was then saturated when the mass transfer resistance through the growing aluminum chloride layer became too large for the HCl to reach the aluminum oxide wall and continue the reaction. This pore-scale model was coupled to the reactor-scale model using a concentration-dependent diffusion coefficient, resulting in much more accurate predictions (approximately half the sum square error of the aforementioned reactor-scale model that produced good agreement for stock and CCC) for a variety of operating conditions. The amount of water vapor or relative humidity (RH) in the flow during a reactor experiment was determined to influence HCl uptake. Experiments were performed to understand the interaction of gaseous HCl with aluminum surfaces in the presence of water vapor. The results show that increasing levels of RH increased the capacity of aluminum to adsorb HCl but decreased the capacity of Nomex to uptake HCl. A series of tests were performed on individual aluminum samples after they had been saturated with a fixed concentration of HCl in dry air conditions with the goal of determining how their HCl uptake capacity changes after various treatments with water relative to the original saturation tests. HCl-saturated aluminum samples subjected to a second dry air flow at the same HCl concentration as the original test had an uptake of 23.5% of the original sample with no treatment in between. Saturated aluminum samples subjected to an air flow with a RH of 90% in between tests had an uptake of 35.6% of the original. Saturated aluminum samples submerged in distilled water for 12 hours in between tests had an uptake of 82.2% of the original sample. Previously saturated aluminum tested with HCl and a 50% RH air flow resulted in similar uptake characteristics in multiple repeated tests. The experiments show the profound effect water vapor has on HCl uptake onto aluminum surfaces. In the samples subjected to water vapor or liquid water, capillary condensation and capillary diffusion alters the transport of HCl significantly. A model was proposed that developed a relationship between RH and the coefficient of HCl diffusion in aluminum chloride. This produced an “S-shaped” curve with diffusion coefficient as a function of RH, with 45% RH represented as the point where the diffusion coefficient is halfway between no water saturation and 100% water saturation in the aluminum chloride product layer. No difference in uptake characteristics for the experiment or model were realized between 50% and 62% RH. The results from the large-scale microgravity experiment, Saffire, are discussed as they pertain to the fate of HCl throughout a spacecraft. HCl was released, both as a standalone event, and in concurrence with the burning of a structured cloth. These events only produced a small response in the far field HCl sensor, while a PMMA burn that did not produce HCl had a significantly greater response. A ground-based large-scale facility was constructed to flow acid gas at the scale and configuration realized in the Saffire experiments. A CFD model of this duct was constructed to test kinetic parameters developed in this work at a larger scale and different geometric configuration and to predict the results of the large-scale facility. The models developed in this work were used to interpret the results of the microgravity tests and lead the discussion on what further experiments and models are needed in order to predict the fate of acid gas in a spacecraft environment. To summarize, the major contributions of this work are as follows: the capacity to uptake HCl, with and without the presence of water vapor, was measured for a variety of real spacecraft surfaces. Several different models (single reactor-scale, multiscale, spacecraft-scale) were developed and with the aid of modeling, the rate of uptake for those surfaces was also predicted and validated. The kinetic parameters determined from the small-scale reactor experiments and models were used to predict large-scale and microgravity tests. Conclusions from this research will be used in the design of spacecraft vehicles and large-scale microgravity fire safety experiments. The models built by this work will aid designers in sensor placement and could be used to predict acid gas transport from fires in partial gravity, as would be seen in Lunar and Martian habitats.
Reaction kinetic measurements on the hydrochlorination of SiCl4 and metallurgical grade (m.g.) silicon metal were made at a wide range of experimental variables. The effect of pressure on the reaction rate was studied at 25 psig, 100 psig, 150 psig and 200 psig, respectively. Results of these experiments show a large pressure effect on the hydrochlorination reaction. As expected, higher pressures produce a higher equilibrium SiHC13 conversion, since the hydrochlorination reaction results in a net volume contraction as product SiHC1 is formed. However, the reaction rate, namely, the rate at which the hydrochlorination reaction reaches its equilibrium SiHC13 conversion, was found to be much faster at low pressures.
The Goddard Space Flight Center (GSFC) two-dimensional model of stratospheric photochemistry and dynamics has been used to calculate the O3 response to stratospheric aircraft (high-speed civil transport (HSCT)) emissions. The sensitivity of the model O3 response was examined for systematic variations of five parameters and two reaction rates over a wide range, expanding on calculations by various modeling groups for the NASA High Speed Research Program and the World Meteorological Organization. In all, 448 model runs were required to test the effects of variations in the latitude, altitude, and magnetitude of the aircraft emissions perturbation, the background chlorine levels, the background sulfate aerosol surface area densities, and the rates of two key reactions. No deviation from previous conclusions concerning the response of O3 to HSCTs was found in this more exhaustive exploration of parameter space. Maximum O3 depletions occur for high-altitude, low altitude HSCT perturbations. Small increases in global total O3 can occur for low-altitude, high-altitude injections. Decreasing aerosol surface area densities and background chlorine levels increases the sensitivity of model O3 to the HSCT perturbations. The location of the aircraft emissions is the most important determinant of the model response. Response to the location of the HSCT emissions is not changed qualitatively by changes in background chlorine and aerosol loading. The response is also not very sensitive to changes in the rates of the reactions NO + HO2 yields NO2 + OH and HO2 + O3 yields OH + 2O2 over the limits of their respective uncertainties. Finally, levels of lower stratospheric HO(sub x) generally decrease when the HSCT perturbation is included, even though there are large increases in H2O due to the perturbation.
The extent, and some of the consequences, of nonequilibrium effects in the inviscid orbiter flow field are examined. Calculations have been made for the reacting flow along selected streamlines in the windward plane of symmetry of a blunted delta shape representative of a shuttle orbiter. An approximate model of the pressure field in the windward plane of symmetry of a blunted delta wing was constructed from wind tunnel measurements of the body surface pressure distribution and approximate calculations of the shock shape and streamline locations. Numerical solutions for the flow along streamlines were obtained for this pressure distribution, and solutions were obtained for both equilibrium and finite-reaction rate, nonequilibrium chemical reactions. Since the streamtube solutions are much easier to compute than fully three-dimensional, reacting flows, results were obtained for a wide range of conditions. The validity of this approach is based on the assumptions that the pressure field is insensitive to nonequilibrium effects and that the flow along streamlines is quasi-one dimensional. Numerical solutions have been used to verify this approach for other blunted slender shapes such as spherically blunted cones.
Chemcial and physical processes which may give rise to isotope fractionation of molecular species in the atmospheres of both earth and other planets are reviewed, along with observations of isotopically substituted molecules in planetary atmospheres. Mechanisms for production of isotope fractionation considered include atmospheric escape and the effect of isotope substitution on equilibrium constants (including those of phase changes), photolysis rates, and chemical reaction rates. The isotopes considered for compounds in the terrestrial atmosphere include D, T, C-13, C-14, N-15, O-18, and S-34. Compounds for which data about isotopic composition in the terrestrial atmosphere are summarized include CO, CO2, O3, N2O, NH3, SO2, H2S, H2O, H, H2, and CH4. Planetary atmospheres discussed include those of Venus, Mars, Jupiter, Saturn, Uranus, and Titan; isotopes reviewed are D, C-13, N-15, and O-18. Suggestions for additional research in the area of isotopically substituted molecules in atmospheres are offered.
The objective was to obtain accurate potential energy surfaces (PES's) for a number of reactions which are important in the H/N/O combustion process. The interest in this is centered around the design of the SCRAM jet engine for the National Aerospace Plane (NASP), which was envisioned as an air-breathing hydrogen-burning vehicle capable of reaching velocities as large as Mach 25. Preliminary studies indicated that the supersonic flow in the combustor region of the scram jet engine required accurate reaction rate data for reactions in the H/N/O system, some of which was not readily available from experiment. The most important class of combustion reactions from the standpoint of the NASP project are radical recombinaton reactions, since these reactions result in most of the heat release in the combustion process. Theoretical characterizations of the potential energy surfaces for these reactions are presented and discussed.
The flame velocities were measured of 20 percent oxygen and 80 percent carbon monoxide mixtures containing either light water or heavy water. The flame velocity increased from 34.5 centimeters per second with no added water to about 104 centimeters per second for a 1.8 percent addition of light water and to 84 centimeters per second for an equal addition of heavy water. The addition of heavy water caused greater increases in flame velocity with equilibrium hydrogen-atom concentration than would be predicted by the Tanford and Pease square-root relation. The ratio of the flame velocity of a mixture containing light water to that of a mixture containing heavy water was found to be 1.4. This value is the same as the ratio of the reaction rate of hydrogen to that of deuterium and oxygen. A ratio of reaction rates of 1.4 would also be required for the square-root law to give the observed ratio of flame-velocity changes.
The objective of this work was to outline and apply an efficient and accurate parametric un-certainty propagation approach to the analysis of convective heating on an Earth entry vehicle concept. The described approach was based on Least Angle Regression used to solve a sparse and underdetermined linear system in the point-collocation non-intrusive polynomial chaos surrogate method. This approach involved an iterative process to computing the non-zero terms of the underlying polynomial chaos model using only enough samples to converge uncertainty interval predictions and Sobol index values based global nonlinear sensitivity estimates. The Earth entry vehicle was analyzed at three points along a representative trajectory for a Mars return mission. 329 sources of uncertainty were identified in the computational fluid dynamics model used to predict the forebody convective heating. These included uncertainty in flow field chemical rates, collision integrals, heats of formation, surface finite rate char model reaction rates, wall roughness height, and the turbulent Schmidt number. Results from this study showed that convective heating uncertainty as high as 50% of the nominal was predicted with only about 50 evaluations of the computational model. This was far fewer than would be required for a sampling-based approach or a full basis polynomial chaos model, which would have required over 50,000 samples. Additionally, results showed that over 90% of the total convective heating uncertainty was due to uncertainty in the N2catalytic rate on the surface, while the remainder of the uncertainty was attributed to the turbulent Schmidt number and the wall roughness uncertainties.