Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reactive Transport Modeling”

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.

At least 91 records · Page 5

Radicals and Reservoirs in the GMI Chemistry and Transport Model: Comparison to Measurements

We have used a three-dimensional chemistry and transport model (CTM), developed under the Global Modeling Initiative (GMI), to carry out two simulations of the composition of the stratosphere under changing halogen loading for 1995 through 2030. The two simulations differ only in that one uses meteorological fields from a general circulation model while the other uses meteorological fields from a data assimilation system. A single year's winds and temperatures are repeated for each 36-year simulation. We compare results from these two simulations with an extensive collection of data from satellite and ground-based measurements for 1993-2000. Comparisons of simulated fields with observations of radical and reservoir species for some of the major ozone-destroying compounds are of similar quality for both simulations. Differences in the upper stratosphere, caused by transport of total reactive nitrogen and methane, impact the balance among the ozone loss processes and the sensitivity of the two simulations to the change in composition.

Douglass, Anne R.↗

Contrast Between 1992 and 1997 High Latitude Spring Haloe Observations of Lower Stratospheric HCl

HCl measurements from HALOE in the northern hemisphere during mid-May 1997 revealed vortex fragments in which the chlorine reservoir partitioning was strongly pushed toward HCl (approx. 90% HCl, approx. 10% ClONO2), similar to partitioning previously observed in the Antarctic vortex region. In contrast, observations of ClONO2 and HCl in the northern polar spring, 1992, and in other years, show these species established the balance typical for gas phase photochemical reactions in this region (approx. 60% HCl, approx. 40% ClONO2). Annually, chlorine reservoirs in the winter lower stratosphere polar vortex are converted to chlorine radicals via heterogeneous reactions on particle surfaces at very cold temperatures (less than about 200 K). As temperatures warm in spring, the heterogeneous processes become insignificant compared with gas phase reactions, and the chlorine reservoirs are reformed. Measurements through the northern winter/spring in 1992 show rapid formation of ClONO2, followed by steady loss of ClONO2 and increasing HCl. Although ClONO2 measurements are not available for 1997, the HCl increase in 1997 is observed to be much more rapid and the eventual HCl mixing ratio is about 50% greater than that of 1992. The observations are examined through comparison with the Goddard three-dimensional chemistry and transport model. This model utilizes winds and temperatures from the Goddard Earth Observing System Data Assimilation System and a complete integration scheme for stratospheric photochemistry. Analysis of the evolution of HCl and ClONO2 shows that the observed difference in the overall rate of HCl formation is explained by the sensitivity of the gas-phase chemistry to the ozone mixing ratio and the temperature. The results show that the model accurately simulates HCl and ClONO2 evolution during these two winters. Model validity is further supported by comparisons with O3 and reactive nitrogen species NO and NO2. This analysis provides a sensitive test of the lower stratospheric chlorine photochemistry, particularly because the analysis considers constituent evolution at a time when the HCl and ClONO2 are far from a photochemical stationary state.

Douglass, A. R.↗

Plume and wake dynamics, mixing, and chemistry behind an HSCT aircraft

The chemical evolution and mixing and vortical motion of a High Speed Civil Transport's engine exhausts must be analyzed in order to track the gas and its speciation as emissions are mixed to atmospheric scales. Attention is presently given to an analytic model of the wake dynamical processes which accounts for the roll-up of the trailing vorticity, its breakup due to the Crow instability, and the subsequent evolution and motion of the reconnected vorticity. The concentrated vorticity is noted to wrap up the buoyant exhaust and suppress its continued mixing and dilution. The species tracked encompass those which could be heterogeneously reactive on the surfaces of the condensed ice particles, and those capable of reacting with exhaust soot particle surfaces to form active contrail and/or cloud condensation nuclei.

Miake-Lye, R. C.↗

Revisions to Photochemical Data for Use in Atmospheric Modeling

Solar and stellar flux incident on an atmosphere can cause molecules to dissociate into highly reactive species and allows for photochemical processes to play a fundamental role in atmospheric chemistry. While models have tried to simulate such processes, they are extremely sensitive to photoabsorption cross-sections and quantum yields: two parameters that are important in determining the photodissociation rate, and hence the lifetime, of atmospheric compounds. Obtaining high-resolution and current data for these parameters is therefore highly desirable. Due to this, database and literature searches for high-quality cross-sections and quantum yields were performed and compiled for KINETICS, a Caltech/JPL Chemical Transport Model that can be used in modeling planetary atmospheres. Furthermore, photodissociation rates determined by running a Titan 1-D model were used to verify the completeness of these latest revisions.

Chemical Transport Model (CTM)↗

A Simulation of Bromoform's Contribution to Stratospheric Bromine

Many chlorinated and brominated compounds that are inert in the troposphere are destroyed in the stratosphere and act as an in-situ source of stratospheric reactive chlorine and bromine. Other halogenated compounds that are reactive in the troposphere might contribute to the stratosphere's halogen budget in two ways. First, like their unreactive companions, rapid convective transport might carry them to the upper troposphere and make them available for subsequent advection by the mean circulation into the stratosphere before they are oxidized or photolyzed. Second, it is more likely that they are destroyed in the troposphere, and the chlorine and bromine that is released might then be transported to the stratosphere. We evaluate the relative influence of these processes on stratospheric bromine in a three-dimensional chemistry and transport model which simulates the distribution of bromoform (CHBr3). CHBr3 is parameterized as a short-lived, ocean-surface source gas whose destruction by photolysis and reaction with hydroxyl (OH) in the troposphere and stratosphere yields inorganic bromine (Br(sub y)). Many of the observed features of CHBr3 are simulated well, and comparisons with observations are used to show that the model represents aspects of transport in the upper troposphere and lower stratosphere that are critical to the evaluation. In particular, the model maintains the observed troposphere-stratosphere distinctness in transport pathways and reproduces the observed seasonal dependence of the mixture of air in the middle- and high-latitude lowermost stratosphere. We estimate that adding CHBr3 to models which already include the long-lived organic brominated compounds (halons and methyl bromide) will increase the simulated stratospheric mass of Br(sub y) by about 15 percent. In-situ stratospheric destruction of CHBr3 produces Br(sub y) in amounts which are comparable to that transported into the stratosphere after photolysis and oxidation of CHBr3 in the troposphere. In our simulations the mass of Br(sub y) produced from the destruction of CHBr3 does not exceed the mass of Br(sub y) produced from the destruction of the long-lived compounds at any level in the stratosphere. However, Br(sub y) from the loss of CHBr3 accounts for approximately one-third of the total Br(sub y) in the lowest kilometer of the stratosphere.

Nielsen, J. Eric↗

Pulse Detonation Engine Modeled

Pulse Detonation Engine Technology is currently being investigated at Glenn for both airbreathing and rocket propulsion applications. The potential for both mechanical simplicity and high efficiency due to the inherent near-constant-volume combustion process, may make Pulse Detonation Engines (PDE's) well suited for a number of mission profiles. Assessment of PDE cycles requires a simulation capability that is both fast and accurate. It should capture the essential physics of the system, yet run at speeds that allow parametric analysis. A quasi-one-dimensional, computational-fluid-dynamics-based simulation has been developed that may meet these requirements. The Euler equations of mass, momentum, and energy have been used along with a single reactive species transport equation, and submodels to account for dominant loss mechanisms (e.g., viscous losses, heat transfer, and valving) to successfully simulate PDE cycles. A high-resolution numerical integration scheme was chosen to capture the discontinuities associated with detonation, and robust boundary condition procedures were incorporated to accommodate flow reversals that may arise during a given cycle. The accompanying graphs compare experimentally measured and computed performance over a range of operating conditions for a particular PDE. Experimental data were supplied by Fred Schauer and Jeff Stutrud from the Air Force Research Laboratory at Wright-Patterson AFB and by Royce Bradley from Innovative Scientific Solutions, Inc. The left graph shows thrust and specific impulse, Isp, as functions of equivalence ratio for a PDE cycle in which the tube is completely filled with a detonable hydrogen/air mixture. The right graph shows thrust and specific impulse as functions of the fraction of the tube that is filled with a stoichiometric mixture of hydrogen and air. For both figures, the operating frequency was 16 Hz. The agreement between measured and computed values is quite good, both in terms of trend and magnitude. The error is under 10 percent everywhere except for the thrust value at an equivalence ratio of 0.8 in the left figure, where it is 14 percent. The simulation results shown were made using 200 numerical cells. Each cycle of the engine, approximately 0.06 sec, required 2.0 min of CPU time on a Sun Ultra2. The simulation is currently being used to analyze existing experiments, design new experiments, and predict performance in propulsion concepts where the PDE is a component (e.g., hybrid engines and combined cycles).

Paxson, Daniel E.↗

Mapping the Potential for Eolian Surface Activity in Grasslands of the High Plains using Landsat Images

There are over 100,000 square kilometers of eolian sand dunes and sand sheets in the High Plains of the central United States. These land-forms may be unstable and may reactivate again as a result of land-use, climate change, or natural climatic variability. The main goal of this thesis was to develop a model that could be used to map an estimate of future dune activity. Multi-temporal calibrated Landsats 5 Thematic Mapper (TM) and 7 Enhanced Thematic Map per Plus (ETM+) NDVI imagery were used in conjunction with the CENTURY vegetation model to correlate vegetation cover to climatic variability. This allows the creation of a predicted vegetation map which, combined with current wind and soil data, was used to create a potential sand transport map for range land in the High Plains under drought conditions.

Gutmann, Ethan Dain↗

Combustion Instabilities Modeled

NASA Lewis Research Center's Advanced Controls and Dynamics Technology Branch is investigating active control strategies to mitigate or eliminate the combustion instabilities prevalent in lean-burning, low-emission combustors. These instabilities result from coupling between the heat-release mechanisms of the burning process and the acoustic flow field of the combustor. Control design and implementation require a simulation capability that is both fast and accurate. It must capture the essential physics of the system, yet be as simple as possible. A quasi-one-dimensional, computational fluid dynamics (CFD) based simulation has been developed which may meet these requirements. The Euler equations of mass, momentum, and energy have been used, along with a single reactive species transport equation to simulate coupled thermoacoustic oscillations. A very simple numerical integration scheme was chosen to reduce computing time. Robust boundary condition procedures were incorporated to simulate various flow conditions (e.g., valves, open ends, and choked inflow) as well as to accommodate flow reversals that may arise during large flow-field oscillations. The accompanying figure shows a sample simulation result. A combustor with an open inlet, a choked outlet, and a large constriction approximately two thirds of the way down the length is shown. The middle plot shows normalized, time-averaged distributions of the relevant flow quantities, and the bottom plot illustrates the acoustic mode shape of the resulting thermoacoustic oscillation. For this simulation, the limit cycle peak-to-peak pressure fluctuations were 13 percent of the mean. The simulation used 100 numerical cells. The total normalized simulation time was 50 units (approximately 15 oscillations), which took 26 sec on a Sun Ultra2.

Paxson, Daniel E.↗

Theoretical study of thermodynamic properties and reaction rates of importance in the high-speed research program

One of the primary goals of NASA's high-speed research program is to determine the feasibility of designing an environmentally safe commercial supersonic transport airplane. The largest environmental concern is focused on the amount of ozone destroying nitrogen oxides (NO(x)) that would be injected into the lower stratosphere during the cruise portion of the flight. The limitations placed on NO(x) emission require more than an order of magnitude reduction over current engine designs. To develop strategies to meet this goal requires first gaining a fundamental understanding of the combustion chemistry. To accurately model the combustor requires a computational fluid dynamics approach that includes both turbulence and chemistry. Since many of the important chemical processes in this regime involve highly reactive radicals, an experimental determination of the required thermodynamic data and rate constants is often very difficult. Unlike experimental approaches, theoretical methods are as applicable to highly reactive species as stable ones. Also our approximation of treating the dynamics classically becomes more accurate with increasing temperature. In this article we review recent progress in generating thermodynamic properties and rate constants that are required to understand NO(x) formation in the combustion process. We also describe our one-dimensional modeling efforts to validate an NH3 combustion reaction mechanism. We have been working in collaboration with researchers at LeRC, to ensure that our theoretical work is focused on the most important thermodynamic quantities and rate constants required in the chemical data base.

Langhoff, Stephen↗

Fire Suppression in Low Gravity Using a Cup Burner

Longer duration missions to the moon, to Mars, and on the International Space Station increase the likelihood of accidental fires. The goal of the present investigation is to: (1) understand the physical and chemical processes of fire suppression in various gravity and O2 levels simulating spacecraft, Mars, and moon missions; (2) provide rigorous testing of numerical models, which include detailed combustion-suppression chemistry and radiation sub-models; and (3) provide basic research results useful for advances in space fire safety technology, including new fire-extinguishing agents and approaches.The structure and extinguishment of enclosed, laminar, methane-air co-flow diffusion flames formed on a cup burner have been studied experimentally and numerically using various fire-extinguishing agents (CO2, N2, He, Ar, CF3H, and Fe(CO)5). The experiments involve both 1g laboratory testing and low-g testing (in drop towers and the KC-135 aircraft). The computation uses a direct numerical simulation with detailed chemistry and radiative heat-loss models. An agent was introduced into a low-speed coflowing oxidizing stream until extinguishment occurred under a fixed minimal fuel velocity, and thus, the extinguishing agent concentrations were determined. The extinguishment of cup-burner flames, which resemble real fires, occurred via a blowoff process (in which the flame base drifted downstream) rather than the global extinction phenomenon typical of counterflow diffusion flames. The computation revealed that the peak reactivity spot (the reaction kernel) formed in the flame base was responsible for attachment and blowoff of the trailing diffusion flame. Furthermore, the buoyancy-induced flame flickering in 1g and thermal and transport properties of the agents affected the flame extinguishment limits.

Takahashi, Fumiaki↗

Fire Suppression in Low Gravity Using a Cup Burner

Longer duration missions to the moon, to Mars, and on the International Space Station increase the likelihood of accidental fires. The goal of the present investigation is to: (1) understand the physical and chemical processes of fire suppression in various gravity and O2 levels simulating spacecraft, Mars, and moon missions; (2) provide rigorous testing of numerical models, which include detailed combustion suppression chemistry and radiation sub-models; and (3) provide basic research results useful for advances in space fire safety technology, including new fire-extinguishing agents and approaches. The structure and extinguishment of enclosed, laminar, methane-air co-flow diffusion flames formed on a cup burner have been studied experimentally and numerically using various fire-extinguishing agents (CO2, N2, He, Ar, CF3H, and Fe(CO)5). The experiments involve both 1g laboratory testing and low-g testing (in drop towers and the KC-135 aircraft). The computation uses a direct numerical simulation with detailed chemistry and radiative heat-loss models. An agent was introduced into a low-speed coflowing oxidizing stream until extinguishment occurred under a fixed minimal fuel velocity, and thus, the extinguishing agent concentrations were determined. The extinguishment of cup-burner flames, which resemble real fires, occurred via a blowoff process (in which the flame base drifted downstream) rather than the global extinction phenomenon typical of counterflow diffusion flames. The computation revealed that the peak reactivity spot (the reaction kernel) formed in the flame base was responsible for attachment and blowoff of the trailing diffusion flame. Furthermore, the buoyancy-induced flame flickering in 1g and thermal and transport properties of the agents affected the flame extinguishment limits.

Takahashi, Fumiaki↗

Formaldehyde Production from Isoprene Oxidation Across NOx Regimes

The chemical link between isoprene and formaldehyde (HCHO) is a strong, non-linear function of NOx (= 27 NO + NO2). This relationship is a linchpin for top-down isoprene emission inventory verification from orbital HCHO column observations. It is also a benchmark for overall mechanism performance with regard to VOC oxidation. Using a comprehensive suite of airborne in situ observations over the Southeast U.S., we quantify HCHO production across the urban-rural spectrum. Analysis of isoprene and its major first-generation oxidation products allows us to define both a "prompt" yield of HCHO (molecules of HCHO produced per molecule of freshly-emitted isoprene) and the background HCHO mixing ratio (from oxidation of longer-lived hydrocarbons). Over the range of observed NOx values (roughly 0.1 - 2 ppbv), the prompt yield increases by a factor of 3 (from 0.3 to 0.9), while background HCHO increases by more than a factor of 2 (from 1.5 to 3.3 ppbv). We apply the same method to evaluate the performance of both a global chemical transport model (AM3) and a measurement-constrained 0-D chemical box model. Both models reproduce the NOx dependence of the prompt HCHO yield, illustrating that models with updated isoprene oxidation mechanisms can adequately capture the link between HCHO and recent isoprene emissions. On the other hand, both models under-estimate background HCHO mixing ratios, suggesting missing HCHO precursors, inadequate representation of later-generation isoprene degradation and/or under-estimated hydroxyl radical concentrations. Moreover, we find that the total organic peroxy radical production rate is essentially independent of NOx, as the increase in oxidizing capacity with NOx is largely balanced by a decrease in VOC reactivity. Thus, the observed NOx dependence of HCHO mainly reflects the changing fate of organic peroxy radicals.

VOC↗

Assimilation of Stratospheric Meteorological and Constituent Observations: A Review

This talk reviews the assimilation of meteorological and constituent observations of the stratosphere. The first efforts to assimilate observations into stratospheric models were during the early 1980s, and a number of research studies followed during the next decade. Since the launch of the Upper Atmospheric Research Satellite (UARS) in 1991, model-assimilated data sets of the stratospheric meteorological state have been routinely available. These assimilated data sets were critical in bringing together observations from the different instruments on UARS as well as linking UARS observations to measurements from other platforms. Using trajectory-mapping techniques, meteorological assimilation analyses are, now, widely used in the analysis of constituent observations and have increased the level of quantitative study of stratospheric chemistry and transport. During the 1990s the use of winds and temperatures from assimilated data sets became standard for offline chemistry and transport modeling. variability in middle latitudes. The transport experiments, however, reveal a set of shortcomings that become obvious as systematic errors are integrated over time. Generally, the tropics are not well represented, mixing between the tropics and middle latitudes is overestimated, and the residual circulation is not accurate. These shortcomings reveal underlying fundamental challenges related to bias and noise. Current studies using model simulation and data assimilation in controlled experimentation are highlighting the issues that must be addressed if assimilated data sets are to be convincingly used to study interannual variability and decadal change. observations. The primary focus has been on stratospheric ozone, but there are efforts that investigate a suite of reactive chemical constituents. Recent progress in ozone assimilation shows the potential of assimilation to contribute to the validation of ozone observations and, ultimately, the retrieval of ozone profiles from space-based radiance measurements. Assimilated data sets provide accurate analyses of synoptic and planetary Scale At the same time, stratospheric assimilation is evolving to include constituent

Rood, Richard B.↗

Combustion synthesis of advanced composite materials

Self-propagating high temperature (combustion) synthesis (SHS), has been investigated as a means of producing both dense and expanded (foamed) ceramic and ceramic-metal composites, ceramic powders and whiskers. Several model exothermic combustion synthesis reactions were used to establish the importance of certain reaction parameters, e.g., stoichiometry, green density, combustion mode, particle size, etc. on the control of the synthesis reaction, product morphology and properties. The use of an in situ liquid infiltration technique and the effect of varying the reactants and their stoichiometry to provide a range of reactant and product species i.e., solids, liquids and gases, with varying physical properties e.g., volatility and thermal conductivity, on the microstructure and morphology of synthesized composite materials is discussed. Conducting the combustion synthesis reaction in a reactive gas environment to take advantage of the synergistic effects of combustion synthesis and vapor phase transport is also examined.

Moore, John J.↗

Suppression Characteristics of Cup-Burner Flames in Low Gravity

The structure and suppression of laminar methane-air co-flow diffusion flames formed on a cup burner have been studied experimentally and numerically using physically acting fire-extinguishing agents (CO2, N2, He, and Ar) in normal earth (lg) and zero gravity (0g). The computation uses a direct numerical simulation with detailed chemistry and radiative heat-loss models. An initial observation of the flame without agent was also made at the NASA Glenn 2.2-Second Drop Tower. An agent was introduced into a low-speed coflowing oxidizing stream by gradually replacing the air until extinguishment occurred under a fixed minimal fuel velocity. The suppression of cup-burner flames, which resemble real fires, occurred via a blowoff process (in which the flame base drifted downstream) rather than the global extinction phenomenon typical of counterflow diffusion flames. The computation revealed that the peak reactivity spot (the reaction kernel) formed in the flame base was responsible for attachment and blowoff phenomena of the trailing diffusion flame. The thermal and transport properties of the agents affected the flame extinguishment limits.

Takahashi, Fumiaki↗

Chlorine Deactivation in the Lower Stratospheric Polar Regions during late Winter: Results from UARS

Recovery from enhanced chlorine conditions in the lower stratospheric polar regions of both hemispheres is investigated using data from the Upper Atmosphere Research Satellite (UARS). Microwave Limb Sounder (MLS) measurements of ClO within the polar vortices are used to infer ClO(sub x) (ClO + 2Cl202) abundances that are then correlated with simultaneous Cryogenic Limb Array Etalon Spectrometer (CLAES) measurements of ClON02 and Halogen Occultation Experiment (HALOE) measurements of HCl obtained starting within 5 days of the end of the MLS and CLAES high-latitude observing periods in each hemisphere. Time series of vortex-averaged mixing ratios are calculated on two potential temperature surfaces (585 K and 465 K) in the lower stratosphere for approximately month-long intervals during late winter: August 17 - September 17, 1992, in the southern hemisphere and February 12 - March 16, 1993, in the northern hemisphere. The observed mixing ratios are adjusted for the effects of vertical transport using diabatic vertical velocities estimated from CLAES tracer data. In the northern hemisphere, the decrease in ClO, is balanced on both surfaces by an increase in ClON02- In the southern hemisphere, continuing polar stratospheric cloud activity prevents ClO from undergoing sustained decline until about September 3. In contrast to the northern hemisphere, there is no significant chemical change in vortex-averaged ClON02 at 465 K, and there is an apparent decrease in ClON02 at 585 K, even after the enhanced ClO abundances have started to recede. Results from the SLIMCAT chemical transport model initialized with UARS data and run with OH + ClO yields HCl + 02 as an 8% channel suggest that the primary recovery product in the south during this time period is not ClON02, but HCl. HALOE HCl mixing ratios are extrapolated back to the time of the MLS and CLAES data. At 585 K, the chlorine budget can be made to balance by extrapolating HCl back to a value of 0.6 parts per billion by volume (ppbv) at the beginning of the study period; at 465 K, the contribution from extrapolated HCl is not sufficient to offset the loss in ClOx, and there is a slight imbalance between the decrease in reactive chlorine and the change in chlorine reservoirs. The difficulty in closing the chlorine budget in the southern hemisphere may arise from complications caused by ongoing activation, incomplete photochemical assumptions, and/or inadequate data quality.

Santee, M. L.↗

Space-Based Sensor Web for Earth Science Applications: An Integrated Architecture for Providing Societal Benefits

There is a significant interest in the Earth Science research and user remote sensing community to substantially increase the number of useful observations relative to the current frequency of collection. The obvious reason for such a push is to improve the temporal, spectral, and spatial coverage of the area(s) under investigation. However, there is little analysis available in terms of the benefits, costs and the optimal set of sensors needed to make the necessary observations. Classic observing system solutions may no longer be applicable because of their point design philosophy. Instead, a new intelligent data collection system paradigm employing both reactive and proactive measurement strategies with adaptability to the dynamics of the phenomena should be developed. This is a complex problem that should be carefully studied and balanced across various boundaries including: science, modeling, applications, and technology. Modeling plays a crucial role in making useful predictions about naturally occurring or human-induced phenomena In particular, modeling can serve to mitigate the potentially deleterious impacts a phenomenon may have on human life, property, and the economy. This is especially significant when one is interested in learning about the dynamics of, for example, the spread of forest fires, regional to large-scale air quality issues, the spread of the harmful invasive species, or the atmospheric transport of volcanic plumes and ash. This paper identifies and examines these challenging issues and presents architectural alternatives for an integrated sensor web to provide observing scenarios driving the requisite dynamic spatial, spectral, and temporal characteristics to address these key application areas. A special emphasis is placed on the observing systems and its operational aspects in serving the multiple users and stakeholders in providing societal benefits. We also address how such systems will take advantage of technological advancement in small spacecraft and emerging information technologies, and how sensor web options may be realized and made affordable. Specialized detector subsystems and precision flying techniques may still require substantial innovation, development time and cost: we have presented the considerations for these issues. Finally, data and information gathering and compression techniques are also briefly described.

Habib, Shahid↗

Modeling of SSME fuel preburner ASI

The Augmented Spark Ignitor (ASI) is a LOX/H2/electrical spark system that functions as an ignition source and sustainer for stable combustion. It is used in the SSME preburner combustor, the SSME main combustion chamber, the J-1 and J-2 engines as well as proposed designs of the Space Transportation Main Engine (STME) main combustor and gas generators. In the SSME it is a long circular cylindrical chamber located along the Main Combustor centerline with a truncated conical dome at the top, which contains two oblique LOX injection ports and two spark plugs offset at 90 degrees. Hydrogen injection is through a number of nearly tangential slots downstream which creates a swirl flow intended to cool the AS1 chamber walls. Past incidents of erosion of the AS1 spark plugs have often led to the need for replacement of these very costly devices. Thus it is desirable to understand the complex reactive flow field within the AS1 both during the initial ignition transient and during the main stage steady state combustion (no sparking).

P Y Liang↗