THE EFFECT OF AMBIENT DISSOCIATION AND SPECIES DIFFUSION ON NONEQUILIBRIUM SHOCK-LAYERS
Effects of ambient dissociation and species diffusion for nonequilibrium shock layers
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.
Effects of ambient dissociation and species diffusion for nonequilibrium shock layers
Blunt body flow studies analyzing the effect of ambient dissociation and species diffusion on nonequilibrium shock layers
This paper presents recent theoretical and experimental results from an ongoing research program that deals with reduced-gravity combustion of bi-component droplets initially in the mm size range or larger. The primary objectives of the research are to study the effects of droplet internal flows, thermal and solutal Marangoni stresses, and species volatility differences on liquid species transport and overall combustion phenomena (e.g., gas-phase unsteadiness, burning rates, sooting, radiation, and extinction). The research program utilizes a reduced-gravity environment so that buoyancy effects are rendered negligible. Use of large droplets also facilitates visualization of droplet internal flows, which is important to this research. This program is a continuation of extensive ground based experimental and theoretical research on bi-component droplet combustion which has been ongoing for several years. The focal point of this research program is a flight experiment (Bi-Component Droplet Combustion Experiment, BCDCE). This flight experiment is still under development. However, supporting ground studies have been performed, and preliminary data have been obtained from flight experiments (Fiber Supported Droplet Combustion Experiment, FSDC-1 and FSDC-2). These flight experiments were performed during the STS-73/USML-2 and STS-94/MSL- I missions. In the experiments, droplets composed of low- and high-volatility species are burned. The low-volatility components are initially present in small amounts. As combustion of a droplet proceeds, the liquid surface mass fraction of the low-volatility component will increase with time, resulting in a sudden and temporary decrease in droplet burning rates as the droplet rapidly heats to temperatures close to the boiling point of the low-volatility component. This decrease in burning rates causes a sudden and temporary contraction of the flame. The decrease in burning rates and the flame contraction can be observed experimentally. Measurements of burning rates as well as the onset time for flame contraction allow effective liquid-phase species diffusivities to be calculated, e.g., using asymptotic theory. A goal of the research is to relate effective liquid species diffusivities to droplet internal flow characteristics. Droplet internal flows will be visualized in future flight and ground-based experiments.
This research deals with reduced-gravity combustion of bi-component droplets initially in the mm size range or larger. The primary objectives of the research are to study the effects of droplet internal flows, thermal and solutal Marangoni stresses, and species volatility differences on liquid species transport and overall combustion phenomena (e.g., gas-phase unsteadiness, burning rates, sooting, radiation, and extinction). The research program utilizes a reduced-gravity environment so that buoyancy effects are rendered negligible. Use of large droplets also facilitates visualization of droplet internal flows, which is important for this research. In the experiments, droplets composed of low- and high-volatility species are burned. The low-volatility components are initially present in small amounts. As combustion of a droplet proceeds, the liquid surface mass fraction of the low-volatility component will increase with time, resulting in a sudden and temporary decrease in droplet burning rates as the droplet rapidly heats to temperatures close to the boiling point of the low-volatility component. This decrease in burning rates causes a sudden and temporary contraction of the flame. The decrease in burning rates and the flame contraction can be observed experimentally. Measurements of burning rates as well as the onset time for flame contraction allow effective liquid-phase species diffusivities to be calculated, e.g., using asymptotic theory. It is planned that droplet internal flows will be visualized in flight and ground-based experiments. In this way, effective liquid species diffusivities can be related to droplet internal flow characteristics. This program is a continuation of extensive ground based experimental and theoretical research on bi-component droplet combustion that has been ongoing for several years. The focal point of this program is a flight experiment (Bi-Component Droplet Combustion Experiment, BCDCE). This flight experiment is under development. However, supporting ground-based studies have been performed. This ground-based research is summarized below.
This research deals with reduced-gravity combustion of bi-component droplets initially in the mm size range or larger. The primary objectives of the research are to study the effects of droplet internal flows, thermal and solutal Marangoni stresses, and species volatility differences on liquid species transport and overall combustion phenomena (e.g., gas-phase unsteadiness, burning rates, sooting, radiation, and extinction). The research program utilizes a reduced gravity environment so that buoyancy effects are rendered negligible. Use of large droplets also facilitates visualization of droplet internal flows, which is important for this research. In the experiments, droplets composed of low- and high-volatility species are burned. The low-volatility components are initially present in small amounts. As combustion of a droplet proceeds, the liquid surface mass fraction of the low-volatility component will increase with time, resulting in a sudden and temporary decrease in droplet burning rates as the droplet rapidly heats to temperatures close to the boiling point of the low-volatility component. This decrease in burning rates causes a sudden and temporary contraction of the flame. The decrease in burning rates and the flame contraction can be observed experimentally. Measurements of burning rates as well as the onset time for flame contraction allow effective liquid-phase species diffusivities to be calculated, e.g., using asymptotic theory. It is planned that droplet internal flows will be visualized in flight and ground-based experiments. In this way, effective liquid species diffusivities can be related to droplet internal flow characteristics. This program is a continuation of extensive ground-based experimental and theoretical research on bi-component droplet combustion that has been ongoing for several years. The focal point of this program is a flight experiment (Bi-Component Droplet Combustion Experiment, BCDCE). This flight experiment is under development. However, supporting ground-based studies have been performed. Some of the most recent ground-based research is summarized.
This research deals with reduced-gravity combustion of bi-component droplets initially in the mm size range or larger. The primary objectives of the research are to study the effects of droplet internal flows, thermal and solutal Marangoni stresses, and species volatility differences on liquid species transport and overall combustion phenomena (e.g., gas-phase unsteadiness, burning rates, sooting, radiation, and extinction). The research program utilizes a reduced-gravity environment so that buoyancy effects are rendered negligible. Use of large droplets also facilitates visualization of droplet internal flows, which is important for this research. In the experiments, droplets composed of low- and high-volatility species are burned. The low-volatility components are initially present in small amounts. As combustion of a droplet proceeds, the liquid surface mass fraction of the low-volatility component will increase with time, resulting in a sudden and temporary decrease in droplet burning rates as the droplet rapidly heats to temperatures close to the boiling point of the low-volatility component. This decrease in burning rates causes a sudden and temporary contraction of the flame. The decrease in burning rates and the flame contraction can be observed experimentally. Measurements of burning rates as well as the onset time for flame contraction allow effective liquid-phase species diffusivities to be calculated, e.g., using asymptotic theory. It is planned that droplet internal flows will be visualized in future flight and ground-based experiments. In this way, effective liquid species diffusivities can be related to droplet internal flow characteristics. This program is a continuation of extensive ground based experimental and theoretical research on bi-component droplet combustion that has been ongoing for several years. The focal point of this program is a flight experiment (Bi-Component Droplet Combustion Experiment, BCDCE). This flight experiment is under development. However, supporting studies have been performed. Because of space limitations, only some of the research performed over the last two years (since the 5th Microgravity Combustion Workshop) is summarized here.
Cool flame quenching distances are generally presumed to be larger than those associated with hot flames, because the quenching distance scales with the inverse of the flame propagation speed, and cool flame propagation speeds are often times slower than those associated with hot flames. To date, this presumption has never been put to a rigorous test, because unstirred, non-isothermal cool flame studies on Earth are complicated by natural convection. Moreover, the critical Peclet number (Pe) for quenching of cool flames has never been established and may not be the same as that associated with wall quenching due to conduction heat loss in hot flames, Pe approx. = 40-60. The objectives of this ground-based study are to: (1) better understand the role of conduction heat loss and species diffusion on cool flame quenching (i.e., Lewis number effects), (2) determine cool flame quenching distances (i.e, critical Peclet number, Pe) for different experimental parameters and vessel surface pretreatments, and (3) understand the mechanisms that govern the quenching distances in premixtures that support cool flames as well as hot flames induced by spark-ignition. Objective (3) poses a unique fire safety hazard if conditions exist where cool flame quenching distances are smaller than those associated with hot flames. For example, a significant, yet unexplored risk, can occur if a multi-stage ignition (a cool flame that transitions to a hot flame) occurs in a vessel size that is smaller than that associated with the hot quenching distance. To accomplish the above objectives, a variety of hydrocarbon-air mixtures will be tested in a static reactor at elevated temperature in the laboratory (1g). In addition, reactions with chemical induction times that are sufficiently short will be tested aboard NASA's KC-135 microgravity (mu-g) aircraft. The mu-g results will be compared to a numerical model that includes species diffusion, heat conduction, and a skeletal kinetic mechanism, following the work on diffusion-controlled cool flames by Fairlie et,al., 2000.
A low-gravity environment, in space or ground-based facilities such as drop towers, provides a unique setting for study of combustion mechanisms. Understanding the physical phenomena controlling the ignition and spread of flames in microgravity has importance for space safety as well as better characterization of dynamical and chemical combustion processes which are normally masked by buoyancy and other gravity-related effects. Even the use of so-called 'limiting cases' or the construction of 1-D or 2-D models and experiments fail to make the analysis of combustion simultaneously simple and accurate. Ideally, to bridge the gap between chemistry and fluid mechanics in microgravity combustion, species concentrations and temperature profiles are needed throughout the flame. However, restrictions associated with performing measurements in reduced gravity, especially size and weight considerations, have generally limited microgravity combustion studies to the capture of flame emissions on film or video laser Schlieren imaging and (intrusive) temperature measurements using thermocouples. Given the development of detailed theoretical models, more sophisticated studies are needed to provide the kind of quantitative data necessary to characterize the properties of microgravity combustion processes as well as provide accurate feedback to improve the predictive capabilities of the computational models. While there have been a myriad of fluid mechanical visualization studies in microgravity combustion, little experimental work has been completed to obtain reactant and product concentrations within a microgravity flame. This is largely due to the fact that traditional sampling methods (quenching microprobes using GC and/or mass spec analysis) are too heavy, slow, and cumbersome for microgravity experiments. Non-intrusive optical spectroscopic techniques have - up until now - also required excessively bulky, power hungry equipment. However, with the advent of near-IR diode lasers, the possibility now exists to obtain reactant and product concentrations and temperatures non-intrusively in microgravity combustion studies. Over the past ten years, Southwest Sciences has focused its research on the high sensitivity, quantitative detection of gas phase species using diode lasers. Our research approach combines three innovations in an experimental system resulting in a new capability for nonintrusive measurement of major combustion species. FM spectroscopy or high frequency Wavelength Modulation Spectroscopy (WMS) have recently been applied to sensitive absorption measurements at Southwest Sciences and in other laboratories using GaAlAs or InGaAsP diode lasers in the visible or near-infrared as well as lead-salt lasers in the mid-infrared spectral region. Because these lasers exhibit essentially no source noise at the high detection frequencies employed with this technique, the achievement of sensitivity approaching the detector shot noise limit is possible.
Ionizing radiation produces several radiolytic species such as.OH, e-aq, and H. when interacting with biological matter. Following their creation, radiolytic species diffuse and chemically react with biological molecules such as DNA. Despite years of research, many questions on the DNA damage by ionizing radiation remains, notably on the indirect effect, i.e. the damage resulting from the reactions of the radiolytic species with DNA. To simulate DNA damage by ionizing radiation, we are developing a step-by-step radiation chemistry code that is based on the Green's functions of the diffusion equation (GFDE), which is able to follow the trajectories of all particles and their reactions with time. In the recent years, simulations based on the GFDE have been used extensively in biochemistry, notably to simulate biochemical networks in time and space and are often used as the "gold standard" to validate diffusion-reaction theories. The exact GFDE for partially diffusion-controlled reactions is difficult to use because of its complex form. Therefore, the radial Green's function, which is much simpler, is often used. Hence, much effort has been devoted to the sampling of the radial Green's functions, for which we have developed a sampling algorithm This algorithm only yields the inter-particle distance vector length after a time step; the sampling of the deviation angle of the inter-particle vector is not taken into consideration. In this work, we show that the radial distribution is predicted by the exact radial Green's function. We also use a technique developed by Clifford et al. to generate the inter-particle vector deviation angles, knowing the inter-particle vector length before and after a time step. The results are compared with those predicted by the exact GFDE and by the analytical angular functions for free diffusion. This first step in the creation of the radiation chemistry code should help the understanding of the contribution of the indirect effect in the formation of DNA damage and double-strand breaks.
Water diffusion experiments were carried out by dehydrating rhyolitic obsidian from Valles Caldera (New Mexico, USA) at 510-980 degrees C. The starting glass wafers contained approximately 0.114 wt% total water, lower than any glasses previously investigated for water diffusion. Weight loss due to dehydration was measured as a function of experiment duration, which permits determination of mean bulk water diffusivity, mean Dw. These diffusivities are in the range of 2.6 to 18 X 10(-14) m2/s and can be fit with the following Arrhenius equation: ln mean Dw (m2/s) = -(25.10 +/- 1.29) - (46,480 +/- 11,400) (J/mol) / RT, except for two replicate runs at 510 degrees C which give mean Dw values much lower than that defined by the above equation. When interpreted according to a model of water speciation in which molecular H2O is the diffusing species with concentration-independent diffusivity while OH units do not contribute to the transport but react to provide H2O, the data (except for the 510 degrees C data) are in agreement with extrapolation from previous results and hence extend the previous data base and provide a test of the applicability of the model to very low water contents. Mean bulk water diffusivities are about two orders of magnitude less than molecular H2O diffusivities because the fraction of molecular H2O out of total water is very small at 0.114 wt% total water and less. The 510 degrees C experimental results can be interpreted as due to slow kinetics of OH to H2O interconversion at low temperatures.
Melt convection, along with species diffusion and segregation on the solidification interface are the primary factors responsible for species redistribution during HgCdTe crystal growth from the melt. As no direct information about convection velocity is available, numerical modeling is a logical approach to estimate convection. Furthermore influence of microgravity level, double-diffusion and material properties should be taken into account. In the present study, HgCdTe is considered as a binary alloy with melting temperature available from a phase diagram. The numerical model of convection and solidification of binary alloy is based on the general equations of heat and mass transfer in two-dimensional region. Mathematical modeling of binary alloy solidification is still a challenging numericial problem. A Rigorous mathematical approach to this problem is available only when convection is not considered at all. The proposed numerical model was developed using the finite element code FIDAP. In the present study, the numerical model is used to consider thermal, solutal convection and a double diffusion source of mass transport.
Mixing efficiency is obtained for a strut fuel injector at hypervelocity flow conditions by using large-eddy simulations (LES), Reynolds-averaged simulations (RAS), and experiments. The injector and flow conditions have been previously investigated by using RAS and experiments as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). Because the fidelity of LES is a strong function of the grid, the mixing efficiency is obtained on two grids, the coarser of which is a factor of two coarser in each of the three dimensions with respect to the fine grid. The RAS uses the two-equation linear eddy viscosity and diffusivity modeling of Menter. In RAS, the species diffusivity model exhibits a strong dependence on the turbulent Schmidt number, which is often adjusted until some metric of engineering interest, such as the mixing efficiency, matches the experimental data. In the absence of experimental data, scale-resolving simulations, such as LES, have been proposed as surrogates for experiments that could provide the data needed to “calibrate” the turbulent Schmidt number in the RAS models. This approach is followed because LES requires significantly more computational resources (CPU, data storage, and time) than RAS, making it prohibitive for use in many engineering applications and specifically for parameter exploration or optimization. Here we examine the mixing efficiency obtained from several RAS with different values of the turbulent Schmidt number, and compare the results with those obtained from the LES and experiments. In addition, the least squares fitting approach was used to demonstrate how to obtain an estimate for the turbulent Schmidt number from LES analytically. These estimates were then used together with prior knowledge about RAS model sensitivity to select a turbulence model that was expected to best match the LES data.
Mixing efficiency is obtained for a strut fuel injector at hypervelocity flow conditions by using large-eddy simulations (LES), Reynolds-averaged simulations (RAS), and experiments. The injector and flow conditions have been previously investigated by using RAS and experiments as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). Because the fidelity of LES is a strong function of the grid, the mixing efficiency is obtained on two grids, the coarser of which is a factor of two coarser in each of the three dimensions with respect to the fine grid. The RAS uses the two-equation linear eddy viscosity and diffusivity modeling of Menter. In RAS, the species diffusivity model exhibits a strong dependence on the turbulent Schmidt number, which is often adjusted until some metric of engineering interest, such as the mixing efficiency, matches the experimental data. In the absence of experimental data, scale-resolving simulations, such as LES, have been proposed as surrogates for experiments that could provide the data needed to “calibrate” the turbulent Schmidt number in the RAS models. This approach is followed because LES requires significantly more computational resources (CPU, data storage, and time) than RAS, making it prohibitive for use in many engineering applications and specifically for parameter exploration or optimization. Here we examine the mixing efficiency obtained from several RAS with different values of the turbulent Schmidt number, and compare the results with those obtained from the LES and experiments. In addition, the least squares fitting approach was used to demonstrate how to obtain an estimate for the turbulent Schmidt number from LES analytically. These estimates were then used together with prior knowledge about RAS model sensitivity to select a turbulence model that was expected to best match the LES data.
Multicomponent gas diffusion is reviewed with particular emphasis on gas flows near solid boundaries-the so-called Kramers-Kistemaker effect. The aim is to derive an appropriate momentum boundary condition which governs many gaseous species diffusing together. The many species' generalization of the traditional single gas condition, either as slip or stick (no-slip), is not obvious, particularly for technologically important cases of lower gas pressures and very dissimilar molecular weight gases. No convincing theoretical case exists for why two gases should interact with solid boundaries equally but in opposite flow directions, such that the total gas flow exactly vanishes. ln this way, the multicomponent no-slip boundary requires careful treatment The approaches discussed here generally adopt a microscopic model for gas-solid contact. The method has the advantage that the mathematics remain tractable and hence experimentally testable. Two new proposals are put forward, the first building in some molecular collision physics, the second drawing on a detailed view of surface diffusion which does not unphysically extrapolate bulk gas properties to govern the adsorbed molecules. The outcome is a better accounting of previously anomalous experiments. Models predict novel slip conditions appearing even for the case of equal molecular weight components. These approaches become particularly significant in view of a conceptual contradiction found to arise in previous derivations of the appropriate boundary conditions. The analogous case of three gases, one of which is uniformly distributed and hence non-diffusing, presents a further refinement which gives unexpected flow reversals near solid boundaries. This case is investigated alone and for aggregating gas species near their condensation point. In addition to predicting new physics, this investigation carries practical implications for controlling vapor diffusion in the growth of crystals used in medical diagnosis (e.g. mercuric iodide) and semiconductors.
The Cool Flame Experiment aims to address the role of diffusive transport on the structure and the stability of gas-phase, non-isothermal, hydrocarbon oxidation reactions, cool flames and auto-ignition fronts in an unstirred, static reactor. These reactions cannot be studied on Earth where natural convection due to self-heating during the course of slow reaction dominates diffusive transport and produces spatio-temporal variations in the thermal and thus species concentration profiles. On Earth, reactions with associated Rayleigh numbers (Ra) less than the critical Ra for onset of convection (Ra(sub cr) approx. 600) cannot be achieved in laboratory-scale vessels for conditions representative of nearly all low-temperature reactions. In fact, the Ra at 1g ranges from 10(exp 4) - 10(exp 5) (or larger), while at reduced-gravity, these values can be reduced two to six orders of magnitude (below Ra(sub cr)), depending on the reduced-gravity test facility. Currently, laboratory (1g) and NASA s KC-135 reduced-gravity (g) aircraft studies are being conducted in parallel with the development of a detailed chemical kinetic model that includes thermal and species diffusion. Select experiments have also been conducted at partial gravity (Martian, 0.3gearth) aboard the KC-135 aircraft. This paper discusses these preliminary results for propane-oxygen premixtures in the low to intermediate temperature range (310- 350 C) at reduced-gravity.
Premixed stoichiometric H2/O2/N2 flames propagating in two-dimensional turbulence were studied using direct numerical simulation (simulations in which all fluid and thermochemical scales are fully resolved) including realistic chemical kinetics and molecular transport. Results are compared with earlier zero-chemistry (flame sheet) and one-step chemistry simulations. Consistent with the simpler models, the turbulent flame with realistic chemistry aligns preferentially with extensive strain rates in the tangent plane and flame curvature probability density functions are close to symmetric with near-zero means. By contrast to simple-chemistry results with non-unity Lewis numbers (ratio of thermal to species diffusivity), local flame structure does not correlate with curvature but rather with tangential strain rate. Turbulent straining results in substantial thinning of the flame relative to the steady unstrained laminar case. Heat release and H2O2 contours remain thin and connected ('flamelet-like') while species including H-atom and OH are more diffuse. Peak OH concentration occurs well behind the peak heat-release zone. The feasibility of incorporating realistic chemistry into full turbulence simulations to address issues such as pollutant formation in hydrocarbon-air flames is suggested.
X-ray photoelectron spectroscopy (XPS) has been used to study the nitridation time and temperature dependence of the nitrogen distribution in thermally nitrided SiO2 films. The XPS data show that the maximum nitrogen concentration near the (SiO(x)N(y)/Si interface is initially at the interface, but moves 20-25 A away from the interface with increasing nitridation time. Computer modeling of the kinetic processes involved is carried out and reveals a mechanism in which diffusing species, initially consisting primarily of nitrogen, react with the substrate, followed by formation of the oxygen-rich oxynitride due to reaction of the diffusing oxygen displaced by the slower nitridation of the SiO2. The data are consistent with this mechanism provided the influence of the interfacial strain on the nitridation and oxidation kinetics is explicitly accounted for.