SEARCH · Search NASA

Results for “Confined Combustion”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Stacking Multiple Ion Captures in The High Performance Antiproton Trap (HiPAT)

The High performance Antiproton Trap (HiPAT) research project was initiated by the Marshall Space Flight Center's propulsion Research Center to examining the fundamental behavior of low energy antiprotons. Stored antiproton would ultimately be used for experimental demonstration of basic propulsive concepts. Matter-antimatter annihilation produces approximately 10(exp 8) MJ/g nearly 10 orders of magnitude more energy per unit mass than chemical based combustion, hence NASA's interest. To achieve containment, HiPAT utilizes a type of electromagnetic bottle know as a Penning trap positioned within an ultrahigh vacuum test section. Recently, the HiPAT hardware configuration has been enhanced to facilitate the capture of multiple normal matter ion burst. This endeavor is often referred to as "stacking" and used to increasing the number of captured particles. A prior normal matter experimental effort, successfully demonstrated the effectiveness of single burst capture. The stacking process is accomplished by manipulating the electric field generated by the confinement electrodes i.e. adjusting the well potential depth. These potential well values are initially configured to maximize the quantity of captured ions per burst; shallow wells with a depth of 100 volt or less (referenced to the incoming ion beam energy) are typically selected. Once captured, a cooling interval is required to reduce the energy of trapped particles below the lower extent of the "trap door" (or leading electrode) ion emitting potential. This is necessary such that a new burst of hot ions can be introduced while preventing those already inside from escaping. The cooling time is driven by a combination of mechanisms such as synchrotron radiation, background gas scattering, and resistive damping in a time scale on the order of minutes. A potential for reducing this hold period is to actively manipulate the electric field shape, using the power supply control system, to produce a deeper potential well. The trapped ions ride down the morphing well to a lower potential energy, leaving margin to once again cycle the "trap door" capturing a new burst while limiting the number of escaping ions. By adjusting the magnitude and timing of the potential well shaping, subsequent shots from the ion beam can be captured on a time scale shorter then that dictated by the standard inject, capture and cool technique. Currently, experimental tests are under way examining the functionally of this system for stacking multiple ions bursts within the HiPAT system.

Martin, James J.↗

NASA NTRS1998

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗

NASA NTRS1975

Analytical study of mixing and reacting three-dimensional supersonic combustor flow fields

An analytical investigation is presented of mixing and reacting hydrogen jets injected from multiple orifices transverse and parallel to a supersonic airstream. The COMOC computer program, based upon a finite-element solution algorithm, was developed to solve the governing equations for three-dimensional, turbulent, reacting, boundary-region, and confined flow fields. The computational results provide a three-dimensional description of the velocity, temperature, and species-concentration fields downstream of hydrogen injection. Detailed comparisons between cold-flow data and results of the computational analysis have established validity of the turbulent-mixing model based on the elementary mixing-length hypothesis. A method is established to initiate computations for reacting flow fields based upon cold-flow correlations and the appropriate experimental parameters of Mach number, injector spacing, and pressure ratio. Key analytical observations on mixing and combustion efficiency for reacting flows are presented and discussed.

Baker, A. J.↗

NASA NTRS1978

Flame characteristics of a NASA contra swirler module

Results are presented for an investigation of the structure of unconfined propane flames stabilized on a NASA contra-swirl can. General flame characterization, flame visualization by sodium chloride, temperature measurements, combustion products analysis, nitric oxide concentrations, and laser anemometer studies are considered. The results are discussed relative to sodium chloride visualization, temperature distributions, chemical composition of flames, nitric oxide distributions, and velocity distributions. It is shown that unconfined propane-air flames can be stabilized downstream of a NASA contra-swirl can under a wide range of airflow rates and fuel-air ratios. However, on increasing fuel input from the lowest at which flames can be stabilized, a compact recirculation zone-stabilized flame is formed. The effect of confinement on changes in flame shape and structure is determined.

Thompson, D.↗

NASA NTRS1984

Calculation of a hollow-cone liquid spray in a uniform air stream

Fluid dynamic computer codes for the simulation of flows in gas turbine engine combustion systems are being developed. NASA is currently sponsoring a two-phase program for the evaluation of the performance of current codes, taking into account also an improvement of accuracy, if needed. The present investigation forms a part of this program. The numerical technique used includes a Lagrangian spray model for liquid fuels. The spray model, in conjunction with the turbulence model, determines the distribution of fuel in the burning zone of the combustor. The numerical technique was applied to a hollow-cone pressure atomizer spraying water into a coflowing confined airstream.

Sturgess, G. J.↗

DOE OSTI2023

Validation of a detonation product equation of state for an insensitive high explosive via slab geometry expansion tests

Slab expansion (SLABEX) tests are conducted to validate a process for calibrating the detonation products equation of state (EOS) of a high explosive (HE). The SLABEX tests use rectangular slabs of PBX 9502, a polymer-bonded HE formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene bound with Kel F-800, a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Three PBX 9502 slab thicknesses are examined, each confined symmetrically by two rectangular copper (Cu) plates approximately one-tenth the thickness of the HE slab. For the duration of each experiment, the detonation flow along the central axis of the PBX 9502 slab remains two-dimensional. The lateral flow velocity component of the outer surfaces of the expanding Cu plates is measured, along with the steady axial detonation speed along the central axis of the SLABEX. Hydrodynamic simulations of the Cu plate expansion in the SLABEX geometry, driven by the energy stored in the detonation products by the detonation combustion event, are conducted using a Jones–Wilkins–Lee EOS for the detonation products. This EOS form was recently parameterized for PBX 9502 in the cylinder expansion test geometry using a newly developed calibration technique [Voelkel et al., Combust. Flame 246, 112373 (2022)]. Good agreement between the experiment and prediction is found in each SLABEX test, demonstrating that the detonation product EOS calibration technique produces EOSs that are predictive when applied to other geometries.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

NASA NTRS1977

Analytical study of laser supported combustion waves in hydrogen

A one-dimensional energy equation, with constant pressure and area, was used to model the LSC wave. This equation balances convection, conduction, laser energy absorption, radiation energy loss and radiation energy transport. Solutions of this energy equation were obtained to give profiles of temperature and other properties, as well as the relation between laser intensity and mass flux through the wave. The flow through the LSC wave was then conducted through a variable pressure, variable area streamtube to accelerate it to high speed, with the propulsion application in mind. A numerical method for coupling the LSC wave model to the streamtube flow was developed, and a sample calculation was performed. The result shows that 42% of the laser power has been radiated away by the time the gas reaches the throat. It was concluded that in the radially confined flows of interest for propulsion applications, transverse velocities would be less important than in the unconfined flows where air experiments have been conducted.

Kemp, N. H.↗

NASA NTRS2001

Reusable Solid Rocket Motor Nozzle Joint-4 Thermal Analysis

This study provides for development and test verification of a thermal model used for prediction of joint heating environments, structural temperatures and seal erosions in the Space Shuttle Reusable Solid Rocket Motor (RSRM) Nozzle Joint-4. The heating environments are a result of rapid pressurization of the joint free volume assuming a leak path has occurred in the filler material used for assembly gap close out. Combustion gases flow along the leak path from nozzle environment to joint O-ring gland resulting in local heating to the metal housing and erosion of seal materials. Analysis of this condition was based on usage of the NASA Joint Pressurization Routine (JPR) for environment determination and the Systems Improved Numerical Differencing Analyzer (SINDA) for structural temperature prediction. Model generated temperatures, pressures and seal erosions are compared to hot fire test data for several different leak path situations. Investigated in the hot fire test program were nozzle joint-4 O-ring erosion sensitivities to leak path width in both open and confined joint geometries. Model predictions were in generally good agreement with the test data for the confined leak path cases. Worst case flight predictions are provided using the test-calibrated model. Analysis issues are discussed based on model calibration procedures.

Clayton, J. Louie↗

NASA NTRS2002

Overview of the High Performance Antiproton (HiPAT) Experiment

The annihilation of matter with antimatter represents the highest energy density of any known reaction, producing 10(exp 8) MJ/g, approximately 10 orders of magnitude more energy per unit mass than chemical based combustion. To take the first step towards using this energy for propulsion applications the NASA MSFC Propulsion Research Center (PRC) has initiated a research activity examining the storage of low energy antiprotons. Storage was identified as a key enabling technology since it builds the experience base necessary to understand the handling of antiprotons for virtually all utilization and high-density storage concepts. To address this need, a device referred to as the High Performance Antiproton Trap (HiPAT) is under development at the NASA MSFC PRC. The HiPAT is an electromagnetic system (Penning-Malmberg design) consisting of a 4 Tesla superconductor, a high voltage confinement electrode system (operation up to 20 KV), and an ultra high vacuum test section (operating in the 10(exp -12) torr range). The system was designed to be portable with an ultimate goal of maintaining 10(exp 12) charged particles with a half-life of 18 days. Currently, this system is being experimentally evaluated using normal matter ions which are cheap to produce and relatively easy to handle. These normal ions provide a good indication of overall trap behavior, with the exception of assessing annihilation losses. The ions are produced external to HiPAT using two hydrogen ion sources, with adjustable beam energy and current. Ion are transported in a beam line and controlled through the use of electrostatic optics. These optics serve to both focus and gate the incoming ions, providing microsecond-timed pulses that are dynamically captured by cycling the HiPAT electric containment field like a 'trap door'. The layout of this system more closely simulates the operations expected at an actual antiproton production facility where 'packets' of antiprotons with pulse widths measured in 100's of nanoseconds could be provided. Initial dynamic capture experiments have been performed with both trap and ton source system functioning at approximately 10% of maximum levels, minimizing the potential for extraneous effects. Dynamic trapping techniques demonstrated the successful capture of millions of hydrogen ions with good agreement with the predicted loading based on the timing sequence, trap electric field, and ion beam current. These techniques will be expanded to examine multiple shot capture or 'stacking' to increase the total number of ions stored within HiPAT.

Martin, James J.↗

NASA NTRS1976

Exhaust gas emissions of a vortex breakdown stabilized combustor

Exhaust gas emission data are described for a swirl stabilized continuous combustor. The combustor consists of confined concentric jets with premixed fuel and air in the inner jet and air in the outer jet. Swirl may be induced in both inner and outer jets with the sense of rotation in the same or opposite directions (co-swirl and counter-swirl). The combustor limits NO emissions by lean operation without sacrificing CO and unburned hydrocarbon emission performance, when commercial-grade methane and air fired at one atmosphere without preheat are used. Relative swirl direction and magnitude are found to have significant effects on exhaust gas concentrations, exit temperatures, and combustor efficiencies. Counter-swirl gives a large recirculation zone, a short luminous combustion zone, and large slip velocities in the interjet shear layer. For maximum counter-swirl conditions, the efficiency is low.

Yetter, R. A.↗

DOE OSTI2024

Generation of Large-Volume High-Pressure Plasma by Spatio-Temporal Control of Space Charge

Due to the fundamental limitations of scaling, very little progress has been made towards achieving a large-volume dense non-equilibrium plasma near atmospheric pressures. Commercially available state-of-the art glow-like atmospheric plasma devices for industrial applications are either narrow tubes or wide slits with narrow openings. Often multiple sources are put together in various configurations to process larger surfaces. The traditional methods of exciting electrodes create spatially fixed electric fields. As a result, the space charge at atmospheric pressure tends to be spatially confined resulting in non-uniformity which eventually leads to instability as the discharge is scaled. Theoretical work done under this project showed that it is possible to generate a spatially rotating electric field by exciting a set of electrodes with phase staggered sinusoidal waveforms. The modeling and simulations were done using plasma fluid models. It was shown that such a field can produce a uniform plasma. At the conclusion of the project, experimental proof of concept with an eight-electrode system in various gases (Air, Helium and Argon) was demonstrated. Power measurements and spectral investigation show that the concept can be used to generate a uniform stable plasma. This plasma source has the potential of opening new applications of nonthermal plasma including combustion of carbon-free fuel. The current limitation of the proposed method in scaling to higher volume and pressure is the need for multiple high voltage amplifiers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

SECARB 2020 Semi-Annual Field Review Reports - Anthropogenic Test

The Southern States Energy Board (SSEB) managed the overall Southeast Regional Carbon Sequestration Partnership (SECARB) Phase III project and the Electric Power Research Institute (EPRI) managed all activities associated with the Anthropogenic Test site. Advanced Resources International, Inc. (ARI) assisted EPRI with the management of site activities. EPRI and ARI followed reporting procedures that included monthly and quarterly technical progress reports, as well as periodic reporting on key field activities. Subsequent sections of this report catalog the field activities and are organized in a semi-annual chronology. The SECARB Anthropogenic Test was the largest demonstration of a fully integrated, pulverized coal-fueled carbon capture and storage (CCS) project in the United States as of 2012 and supports a commercial prototype of CO 2 capture; transportation; subsurface storage; and monitoring, verification, and accounting (MVA); and assessment. The demonstration-scale, post-combustion CO 2 capture facility at Alabama Power Company’s Plant Barry diverted flue gas (25MW equivalent) from its Unit #5 coal burning facility and captured the CO 2 using Mitsubishi Heavy Industries (MHI) KM-CDRTM advanced amine technology. The captured CO2 was compressed at Plant Barry and transported by pipeline to the injection location southeast of Citronelle, Alabama. The Citronelle Project’s geologic storage and MVA sites were located on the flanks of the Citronelle Dome, approximately three miles southeast of the city of Citronelle. The injection zone, the Paluxy Formation, a saline formation that occurs at approximately 9,400 feet and overlies the oil production horizon of the Citronelle oilfield, presents a favorable injection reservoir in terms of areal extent and petrophysical characteristics. The confining zone, the basal shale of the overlying Washita-Fredericksburg Formation, is persistent throughout the Citronelle area and possesses the appropriate criteria to act as an effective CO 2 seal.

01 COAL, LIGNITE, AND PEAT↗

NASA NTRSDate not supplied

Physics and Modeling of Quasi-Steady and Forced Shock Trains

The dual-mode scramjet is a propulsive concept for contemporary hypersonic vehicles promising flight cruise speeds upwards of Mach 5. Despite decades of development and a recent resurgence in interest, such engines are not routinely deployed for commercial or military purposes. Instead, they have been confined to the realm of ground-test research and limited one-off flight test missions. One of the key stumbling blocks on the critical path to deployment of a dual-mode scramjet is the optimal design and stable operation of the isolator, a typically constant area duct that separates the scramjet inlet from the combustor. During ramjet mode at Mach numbers from 3 to 5 and during the transition to scramjet mode at higher Mach numbers, the isolator contains the pseudo-shock. Necessary to meet the back-pressure in the combustor, the pseudo-shock is typically divided into a shock dominated shock train region followed by a turbulence driven mixing region. Fluctuations in combustor pressure due to vehicle maneuvering, throttling, or combustion instabilities drive the shock train to shift streamwise within the isolator. In response to sharp pressure fluctuations, the shock train is unable to stabilize within the isolator, and the scramjet unstarts by ejecting the shock train out the inlet. This work addresses three key gaps in understanding of the isolator pseudo-shock. The phenomenon is studied in three experimental facilities to investigate the quasi-steady and unsteady fluid dynamics of the shock train in fueled and un-fueled test conditions. The strong link between boundary layer morphology and shock train structure is studied in detail and used as a tool to evaluate the state-of-the-art definitions and models of the pseudo-shock. Identifying the limitations of these models, a new process driven flux-conserved pseudo-shock model is proposed. This provides a computationally efficient model grounded in the fluid dynamic mechanisms of the pseudo-shock with better predictive capability than the current state-of-the-art across all facilities tested. The unsteady behavior of the shock train is further investigated to identify and understand two distinct forced dynamic modes: mechanically induced, low-rate forcing and high-rate combustion driven unsteadiness that sheds upstream propagating shocks. The proposed process driven model is found to be effective in the the limiting case of high-rate combustion driven unsteadiness by evaluating the model in a shock stationary reference frame. Further analysis identifies the approximately rigid chain motion of the shock train structure to mechanically induced forcing. The time delay between the onset of forced dynamics is found to saturate to a constant value and maximum shock train speed along the isolator with respect to the perturbation rate of change. The result of this work provides dual-mode scramjet designers with improved tools and physical understanding to develop more efficient, robust, and capable propulsion systems for hypersonic vehicles.

JAXA↗

NASA NTRS1988

Confined turbulent fluid-particle flow modeling using multiple-realization particle trajectory schemes

A multiple-realization particle trajectory scheme has been developed and applied to the numerical prediction of confined turbulent fluid-particle flows. The example flows investigated include the vertical pipe upflow experimental data of Tsuji et al. and the experimental data of Leavitt for a coaxial jet flow, comprising a particle-laden central jet and a clean annular jet, into a large recirculation chamber. The results obtained from the numerical scheme agree well with the experimental data, lending confidence to the modeling approach. The multiple-realization particle trajectory turbulent flow modeling scheme is believed to be a more elegant and accurate approach to the extension of single-particle hydrodynamics to dilute multi-particle systems than the more commonly employed two-fluid modeling approach. It is also better able to incorporate additional force items such as lift, virtual mass and Bassett history terms directly into the particle equation of motion as appropriate. This makes it a suitable candidate for particle migration studies and an extension to situations involving liquid particulate phases with possible propulsion applications, such as in spray combustion, follows naturally.

Adeniji-Fashola, A. A.↗

DOE OSTI2025

2D CFD of lean premixed hydrogen–air flame quenching under locomotive engine conditions

Hydrogen (H 2 ) is a promising fuel for reducing emissions in heavy-duty internal combustion engines (ICEs), but its low quenching distance increases the risk of flame propagation into narrow crevice regions, such as piston-liner gaps. This work uses detailed CFD simulations with finite-rate chemistry to investigate premixed H 2 –air flame quenching in a two-dimensional (2D) region consistent with the piston-liner gap of a diesel ICE. Model accuracy was assessed by comparison with experimentally measured quenching distances in an annular stepwise diverging tube (ASDT). A parametric study was conducted to assess the influence of crevice width (0.05–1.18 mm), pressure (50–150 bar), unburned gas temperature (431–573 K), and equivalence ratio ( Φ= 0.3–0.6). Results show that the critical Péclet number for flame survival is no greater than 3.25, consistent with prior literature, even under ultra-lean and high-pressure conditions (Φ ≤ 0.3, > 50 bar). Additionally, flames with Péclet numbers exceeding 6.15 exhibited front wrinkling, suggesting the onset of velocity-driven instabilities and enhanced flame robustness. These findings help define thresholds for flame quenching in confined geometries and support the safe design of H 2 fueled ICEs.

08 HYDROGEN↗

NASA NTRS1999

Methane Pyrolysis and Disposing Off Resulting Carbon

Sabatier/Electrolysis (S/E) is a leading process for producing methane and oxygen for application to Mars ISPP. One significant problem with this process is that it produces an excess of methane for combustion with the amount of oxygen that is produced. Therefore, one must discard roughly half of the methane to obtain the proper stoichiometric methane/oxygen mixture for ascent from Mars. This is a waste of hydrogen, which must be brought from Earth and is difficult to transport to Mars and store on Mars. To reduce the problem of transporting hydrogen to Mars, the S/E process can be augmented by another process which reduces overall hydrogen requirement. Three conceptual approaches for doing this are (i) recover hydrogen from the excess methane produced by the S/E process, (ii) convert the methane to a higher hydrocarbon or other organic with a lower H/C ratio than methane, and (iii) use a separate process (such as zirconia or reverse water gas shift reaction) to produce additional oxygen, thus utilizing all the methane produced by the Sabatier process. We report our results here on recovering hydrogen from the excess methane using pyrolysis of methane. Pyrolysis has the advantage that it produces almost pure hydrogen, and any unreacted methane can pass through the S/E process reactor. It has the disadvantage that disposing of the carbon produced by pyrolysis presents difficulties. The goals of a research program on recovery of hydrogen from methane are (in descending priority order): 1) Study the kinetics of pyrolysis to arrive at a pyrolysis reactor design that produces high yields in a confined volume at the lowest possible operating temperature; 2) Study the kinetics of carbon burnoff to determine whether high yields can be obtained in a confined volume at acceptable operating temperatures; and 3) Investigate catalytic techniques for depositing carbon as a fine soot which can be physically separated from the reactor. In the JPL program, we have made significant measurements in regard to goal 1, cursory measurements in regard to goal 2, and would plan to pursue goal 3 if additional resources are secured.

Sharma, P. K.↗

DOE OSTI2021

Experimental and theoretical analysis of carbon driven detonation waves in a heterogeneously premixed Rotating Detonation Engine

Coal dust explosions can be hazardous; however, they can also generate a significant rise in stagnation pressure if adequately harnessed. Rotating detonation combustors seek to take advantage of the stagnation pressure rise phenomenon in a more sustained and controlled manner via confinement to a physical annulus, leading to increased overall thermodynamic efficiency. Here this investigation presents an analysis of detonations fueled by Carbon Black, a solid particulate consisting of virtually pure carbon molecules and lean Hydrogen-Air mixtures. It is realized that with the addition of Carbon Black, an increase of lean mixture detonability and detonation velocities extending the operating limit over that of a pure hydrogen-air mixture is experienced. For all testing conditions, the total equivalence ratio is held at φ = 1, while the fuel mixture's carbon mass fraction is increased from 0 to 0.7 while the hydrogen is decreased. Detonation wave velocities are extracted from high-speed imaging through applying a Discrete Fourier Transform algorithm to determine changes to the wave speed as Carbon Black particles are introduced. As a result, due to the addition of Carbon Black as an auxiliary fuel source, detonations were formed instead of deflagrations in operating conditions where one would expect deflagrations at the same hydrogen-air equivalence ratios without Carbon Black addition. The detonation formation provides evidence that the coal particles are reacting within the detonation wave in a large enough capacity to support a detonation wave within the annulus. Furthermore, the wave speed is shown to increase with the additional of carbon particles. At a constant global equivalence ratio, the detonation wave velocities were found to decrease with hydrogen's incremental replacement with coal particles. Whereby, through a theoretical comparison of the heat of combustion as computed from the experimentally derived detonation wave velocities, a linear relationship of the two was shown to exist. Therefore, the heat of combustion has the potential to describe an operational limit to sustaining a detonation wave.

42 ENGINEERING↗

NASA NTRS1995

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗