Search NASA⌕ Search

SEARCH · Search NASA

Results for “RECOMBINATION COEFFICIENT”

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.

244 records · Page 14

Viscous reacting flows with wall slip and catalysis applied to spheres in arc jets and flight

The influence of wall slip and catalytic atom-recombination on the flow field and wall heat flux are calculated for high altitude flight and arc jet flow conditions. Boundary equations, which include velocity slip, temperature jump, and wall catalytic atom recombination, are coupled to the viscous reacting multicomponent Navier-Stokes equations. These equations are solved using a time-dependent finite difference technique applied to spheres in an arc jet flow (Reynolds number of 550) and a high altitude flight case representative of the Space Shuttle Orbiter (Reynolds number of 450). The results indicate that catalysis strongly influences the temperature jump, but not the velocity slip. Slip increases the atom fraction and temperature at both the wall and the flow field. Likewise, the shock stand-off distance, the wall heat flux, and friction coefficient are increased over the nonslip cases. The reacting gas calculations confirm the chemically frozen nature of the shock layer in arc jet flows.

Scott, C. D.↗

Photochemical modeling of Titan's atmosphere

We have developed a new photochemical model of Titan's atmosphere which includes all the important compounds and reactions in spherical geometry from the surface to 1240 km. Compared to the previous model of Yung et al. (1984), the most significant recent change in the reactions used is the updated methane dissociation scheme (Mordaunt et al. 1993). Moreover, the transfer of the solar radiation in the atmosphere and the photolysis rates have been calculated by using a Monte Carlo code. Finally, the eddy diffusion coefficient profile is adjusted in order to fit the mean vertical distribution of HCN retrieved from millimeter groundbased observations of Tanguy et al. (1990); using new values for the boundary flux of atomic nitrogen (Strobel et al. 1992). We have run the model in both steady-state and diurnal modes, with 62 species involved in 249 reactions. There is little difference between diurnal and steady-state results. Overall our results are in a closer agreement with the abundances inferred from the Voyager infrared measurements at the equator than the Yung et al. results. We find that the catalytic scheme for H recombination invoked by Yung et al. only slightly improves the model results and we conclude that this scheme is not essential to fit observations.

Toublanc, D.↗

Comprehensive Laboratory Measurements Resolving the LMM Dielectronic Recombination Satellite Lines in Ne-like Fe XVII Ions

We investigated experimentally and theoretically dielectronic recombination (DR) populating doubly excited configurations 3l3l¢ (LMM) in Fe XVII, the strongest channel for soft X-ray line formation in this ubiquitous species. We used two different electron beam ion traps and two complementary measurement schemes for preparing the Fe XVII samples and evaluating their purity, observing negligible contamination effects. This allowed us to diagnose the electron density in both EBITs. We compared our experimental resonant energies and strengths with those of previous independent work at a storage ring as well as those of configuration interaction, multiconfiguration Dirac–Fock calculations, and many-body perturbation theory. This last approach showed outstanding predictive power in the comparison with the combined independent experimental results. From these we also inferred DR rate coefficients, unveiling discrepancies from those compiled in the OPEN-ADAS and AtomDB databases.

Filipe Grilo↗

Calculations of the displacement damage and short-circuit current degradation in proton irradiated (AlGa)As-GaAs solar cells

A theoretical model for computing the displacement damage defect density and the short-circuit current (I sub sc) degradation in proton-irradiated (AlGa)As-GaAs p-n junction solar cells is presented. Assumptions were made with justification that the radiation induced displacement defects form an effective recombination center which controls the electron and hole lifetimes in the junction space charge region and in the n-GaAs active layer of the irradiated GaAs p-n junction cells. The degradation of I sub sc in the (AlGa)As layer was found to be negligible compared to the total degradation. In order to determine the I sub sc degradation, the displacement defect density, path length, range, reduced energy after penetrating a distance x, and the average number of displacements formed by one proton scattering event were first calculated. The I sub sc degradation was calculated by using the electron capture cross section in the p-diffused layer and the hole capture cross section in the n-base layer as well as the wavelength dependent absorption coefficients. Excellent agreement was found between the researchers calculated values and the measured I sub sc in the proton irradiated GaAs solar cells for proton energies of 100 KeV to 10 MeV and fluences from 10 to the 10th power p/square cm to 10 to the 12th power p/square cm.

Yeh, C. S.↗

Spectral identification/elimination of molecular species in spacecraft glow

Computer models of molecular electronic and vibrational emission intensities were developed. Known radiative emission rates (Einstein coefficients) permit the determination of relative excited state densities from spectral intensities. These codes were applied to the published spectra of glow above shuttle surface and to the Spacelab 1 results of Torr and Torr. The theoretical high-resolution spectra were convolved with the appropriate instrumental slit functions to allow accurate comparison with data. The published spacelab spectrum is complex but N2+ Meinel emission can be clearly identified in the ram spectrum. M2 First Positive emission does not correlate well with observed features, nor does the CN Red System. Spectral overlay comparisons are presented. The spectrum of glow above shuttle surfaces, in contrast to the ISO data, is not highly structured. Diatomic molecular emission was matched to the observed spectral shape. Source excitation mechanisms such as (oxygen atom)-(surface species) reaction product chemiluminescence, surface recombination, or resonance fluorescent re-emission will be discussed for each tentative assignment. These assignments are the necessary first analytical step toward mechanism identification. Different glow mechanisms will occur above surfaces under different orbital conditions.

Green, B. D.↗

Theoretical Studies of Important Processes in Planetary and Comet Atmospheres

Using theoretical quantum chemical calculations, I have successfully described the dissociative recombination (DR) of O2(+) leading to the excited S-1 state of atomic oxygen, the upper state of the well known green line emission. The process is described by O2(+) + e(-) yields O(S-1)) + O(D-1) (1) where e(-) is an electron and the product oxygen atoms are both excited. This process is important in the atmospheres of Venus, Mars and Earth. I have shown in prior calculations that only one repulsive potential curve of O2, f(sup l)Sigma(sub u, sup +) can generate O(S-1) from DR of the lowest vibrational levels of O2(+). However, in the prior results, the calculated quantum yield (i.e. the number of O(S-1) atoms produced for every two product atoms) from the O2(+) v = 0 level was smaller than the laboratory and atmospheric measurements by more than an order of magnitude. Including only direct recombination, the calculated quantum yield for O(S-1) is only 0.0016. In a calculation that accounts for both direct and indirect recombination, the quantum yield is 0.0012. The range of experimentally determined quantum yields is between 0.01 and 0.23. Because of this large difference between the theoretical and experimental quantum yields, it was thought for some time that the ionospheric and laboratory O2(+) must be vibrationally excited since for excited levels, theory gave quantum yields that are similar to experimental yields. It was also suggested that some other process was generating O(S-1) but none could be identified. Under current NASA support, I have found that reaction (1) proceeds via an unusual mechanism. The f(sup 1)Sigma(sub u, sup,+) state does not cross the ion between the turning points of the v = 0 level of the O2(+) ground state. The lack of a favorable crossing leads to a very small calculated DR rate coefficient. However, this mechanism assumes that initial electron capture must occur into the repulsive state that leads to )(S-1). Instead, I have found that initial electron capture occurs mostly into the B(sup 3)Sigma(sup u, sub -) state which crosses the ion between the turning points of the v = 0 ion level and has a large DR rate coefficient. The B state dissociates to O(D-1) and )(P-3). After capture, some of the flux is transferred to the f(sup 1)Sigma(sup u, sub +) state via symmetry mixed intermediate Rydberg states. The neutral Rydberg states are a mixture of Sigma(sub u, sup +)-1 and Sigma(sub u, sup -)-3 symmetry.

Guberman, Steven L.↗

Investigation of N2O Production from 266 and 532 nm Laser Flash Photolysis of O3/N2/O2 Mixtures

Tunable diode laser absorption spectroscopy has been employed to measure the amount of N2O produced from laser flash photolysis of O3/N2/O2 mixtures at 266 and 532 nm. In the 532 nm photolysis experiments very little N2O is observed, thus allowing an upper limit yield of 7 x 10(exp -8) to be established for the process O3 + N2 yield N2O + O2, where O3 is nascent O3 that is newly formed via O(3P(sub J)) + O2 recombination (with vibrational excitation near the dissociation energy of O3). The measured upper limit yield is a factor of approx. 600 smaller than a previous literature value and is approximately a factor of 10 below the threshold for atmospheric importance. In the 266 nm photolysis experiments, significant N2O production is observed and the N2O quantum yield is found to increase linearly with pressure over the range 100 - 900 Torr in air bath gas. The source of N2O in the 266 nm photolysis experiments is believed to be the addition reaction O(1D(sub 2)) + N2 + M yields (k(sub sigma)) N2O + M, although reaction of (very short-lived) electronically excited O3 with N2 cannot be ruled out by the available data. Assuming that all observed N2O comes from the O(1D(sub 2)) + N2 + M reaction, the following expression describes the temperature dependence of k(sub sigma) (in its third-order low-pressure limit) that is consistent with the N2O yield data: k(sub sigma) = (2.8 +/- 0.1) x 10(exp -36)(T/300)(sup -(0-88+0.36)) cm(sup 6) molecule(sup -2)/s, where the uncertainties are 2(sigma) and represent precision only. The accuracy of the reported rate coefficients at the 95% confidence level is estimated to be 30 - 40% depending on the temperature. Model calculations suggest that gas phase processes initiated by ozone absorption of a UV photon represent about 1.4% of the currently estimated global source strength of atmospheric N2O. However, these processes could account for a significant fraction of the oxygen mass-independent enrichment observed in atmospheric N2O, and they appear to be the first suggested photochemical mechanism that is capable of explaining the altitude dependence of the observed mass -independent isotopic signature.

Estupinan, E. G.↗

Development of a Detailed Surface Chemistry Framework in DSMC

Many of the current direct simulation Monte Carlo (DSMC) codes still employ only simple surface catalysis models. These include only basic mechanisms such as dissociation, recombination, and exchange reactions, without any provision for adsorption and finite rate kinetics. Incorporating finite rate chemistry at the surface is increasingly becoming a necessity for various applications such as high speed re-entry flows over thermal protection systems (TPS), micro-electro-mechanical systems (MEMS), surface catalysis, etc. In the recent years, relatively few works have examined finite-rate surface reaction modeling using the DSMC method.In this work, a generalized finite-rate surface chemistry framework incorporating a comprehensive list of reaction mechanisms is developed and implemented into the DSMC solver SPARTA. The various mechanisms include adsorption, desorption, Langmuir-Hinshelwood (LH), Eley-Rideal (ER), Collision Induced (CI), condensation, sublimation, etc. The approach is to stochastically model the various competing reactions occurring on a set of active sites. Both gas-surface (e.g., ER, CI) and pure-surface (e.g., LH, desorption) reaction mechanisms are incorporated. The reaction mechanisms could also be catalytic or surface altering based on the participation of the bulk-phase species (e.g., bulk carbon atoms). Marschall and MacLean developed a general formulation in which multiple phases and surface sites are used and we adopt a similar convention in the current work. Microscopic parameters of reaction probabilities (for gas-surface reactions) and frequencies (for pure-surface reactions) that are require for DSMC are computed from the surface properties and macroscopic parameters such as rate constants, sticking coefficients, etc. The energy and angular distributions of the products are decided based on the reaction type and input parameters. Thus, the user has the capability to model various surface reactions via user-specified reaction rate constants, surface properties and parameters.

Gas↗

Analytical studies of hypersonic viscous dissociated flows

This project primarily dealt with integral boundary-layer solution techniques that are directly applicable to the problem of determining aerodynamic heating rates of hypersonic vehicles like X-33 in the vicinity of stagnation points, windward centerlines, and swept-wing leading edges. The analyses include effects of finite-rate gas chemistry across the boundary layer and finite-rate catalysis of atom recombination at the surface. A new approach for combining the insight afforded by integral boundary-layer analysis with comprehensive (and expensive) computational fluid dynamic (CFD) flowfield solutions of the thin-layer Navier-Stokes equations was developed. The approach extracts CFD derived quantities at the wall and at the boundary layer edge for inclusion in a post-processing boundary-layer analysis. The post-processed data base allows a designer at a workstation to ask and answer the following questions: (1) How much does the heating change if one uses a thermal protection system (TPS) with different catalytic properties than was used in the original CFD solution? (2) How does the heating change when one moves the interface of two different TPS materials with different catalytic efficiencies for the purpose of reducing vehicle weight and expense? The answer to the second question is particularly critical, because abrupt changes from low catalytic efficiency to high catalytic efficiency can lead to localized increase in heating which exceeds the usually conservative estimate provided by a fully catalytic wall assumption. A secondary issue that was addressed involves the prediction of heating levels in the vicinity of sharp corners that are transverse to or aligned with the flow. An example of the first case is heating at the edge of the COMET reentry module. An example of the second case is heating along the side edge of a deflected body flap on an SSV. The difficulty of putting grids in the vicinity of such corners with continuously varying metric coefficients causes problems in CFD predictions. A preliminary theory for prediction that says the heating at the corner is X percent of the heating N boundary-layer thicknesses inboard was developed. This will prove useful to analytically evaluate the possible benefits of rounding the edges of these configurations and defining how much rounding is sufficient.

Inger, George R.↗

Aerothermal Analysis of the Rocket Lab Venus Probe Heatshield

This document provides an overview of the aerothermodynamic analyses performed by the Aerothermodynamics Branch at NASA Langley Research Center for the Rocket Lab Venus Probe (RLVP). In addition to defining the baseline heating environment to the heatshield, this document pursues the experimental validation of key physical models at RLVP-relevant conditions. This experimental validation analysis, which captures the model form uncertainty, is used as one of two primary components of the margin assessment, where the other component is the parametric uncertainty. These model form (experimental) and parametric uncertainty components are used to construct a spatial and time varying margin for the heating to the RLVP heatshield. The margin is evaluated as the sum of the parametric and model form uncertainty components. The model form uncertainty is defined as the difference between the RLVP-relevant measurements and their simulations, using the upper limit uncertainty bounds for both the measurements and simulations in the comparisons. The differences in the dominant physics in the stagnation region and flank lead to the separate RLVP-relevant measurements for assessing the model form uncertainty in these two regions. These regions are addressed as follows: Stagnation Region Heating Environment: For the high-temperature stagnation-region, both the radiative heating and impact of blowing on convective heating are significant, while the impacts of turbulence and roughness are negligible. Coupled radiation and ablation LAURA/HARA solutions with ray-tracing provide the radiative heating over the entire vehicle, including the contributions from the Venus atmosphere and ablation species. Non-ablating LAURA simulations provide the convective heating. During the material-response computation typically used for TPS sizing, this non-ablating convective heating is corrected for the impact of ablation using the blowing correction. Coupled ablation LAURA simulations that capture finite-rate sur-face processes show that this blowing correction may be non-conservative over most of the heatshield. This non-conservatism is due to hydrogen recombination in the finite-rate surface model, which tends to increase the coupled ablation convective heating to near the non-ablating values, therefore making any reduction in the non-ablating value through the blowing correction non-conservative. This non-conservatism due to H catalysis is captured in the parametric component of the margin. The best available ground-test measurements that capture the impact of blowing on stagnation region convective heating, at RLVP-relevant conditions, indicate that the current blowing reduction model is non-conservative by up to 20% at RLVP-relevant blowing rates (the coupled blowing simulations were also non-conservative). Because of the relatively low velocity of the ground tests and the non-Venus atmospheric chemistry, these measurements do not capture the chemistry and therefore do not inform the uncertainty due to H catalysis. However, they do capture the fluid mechanics of blowing. The non-conservatism of the blowing correction implied by these measurements is covered by the model form component of the margin, which leads to total margin values over 50%. For the radiative heating, the shock-tube informed bias approach suggests a model form uncertainty of roughly 20%, while the parametric uncertainty analysis suggests values over 100%. The combined stagnation-point radiation margin of over 100% leads to peak margined radiative heating values of over300 W/cm2, which remains small relative to the peak margined convective heating of nearly 2000 W/cm2. Based on this analysis, at the stagnation point, the peak margined heat rate is 2203 W/cm2 and the margined total heat load is 31.5 kJ/cm2 for the current nominal trajectory. Flank Heating Environment: The forebody flank (and near-shoulder) heating environment is dominated by the impact of turbulence, roughness augmentation, and ablation on the convective heating. An extensive collection of ground test measurements with RLVP-relevant turbulence and roughness is studied to show that the maximum difference between the simulated and measured convective heating is 5%. However, with the exception of the Holden measurements from the 1980s, these measurements do not include roughness elements extending into the supersonic region of the boundary layer, which is likely to occur for RLVP (due to the 45 degree sphere-cone geometry). The interaction between the supersonic flow and roughness could cause convective heating augmentation fundamentally different than for locally subsonic flow. Although these Holden measurements are consistent with the other measurements considered, another path was pursued to assure that the rough-ness height extending into supersonic flow does not fundamentally change the roughness augmentation. This additional path was a computational effort to resolve the roughness elements in the CFD grid, so that the interaction be-tween the roughness elements and locally supersonic flow may be simulated in detail. This roughness-resolved CFD simulation is feasible because of the RLVP forebody TPS’s patterned roughness, which may be approximated analytically, and because of the axisymmetric nominal flow field, which allows a narrow surface region to be simulated and therefore make the computational expense feasible. These grid-resolved roughness simulations, which are performed at actual RLVP flight conditions, result in heating augmentation values that are below the design approach for roughness augmentation. This provides evidence that the design approach for RLVP roughness augmentation is sufficient. Based on this analysis, at this flank or near-shoulder location, the peak margined total heat rate is 2088 W/cm2and the margined total heat load is 26.0 kJ/cm2for the current nominal trajectory. Heat flux, shear, pressure and heat transfer coefficient at the RLVP stagnation point and near shoulder location are evaluated for the entire trajectory, and curved fit to a functional form of F=AρB∞UC∞. These simplified relationships for the nominal and margined aerothermal environments are referred to as aerothermal indicators, and presented at the end of this document.

Christopher O Johnston↗