Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Combustion 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 1,027 records · Page 57

Turbulent reacting flow computations including turbulence-chemistry interactions

A two-equation (k-epsilon) turbulence model has been extended to be applicable for compressible reacting flows. A compressibility correction model based on modeling the dilatational terms in the Reynolds stress equations has been used. A turbulence-chemistry interaction model is outlined. In this model, the effects of temperature and species mass concentrations fluctuations on the species mass production rates are decoupled. The effect of temperature fluctuations is modeled via a moment model, and the effect of concentration fluctuations is included using an assumed beta-pdf model. Preliminary results obtained using this model are presented. A two-dimensional reacting mixing layer has been used as a test case. Computations are carried out using the Navier-Stokes solver SPARK using a finite rate chemistry model for hydrogen-air combustion.

Narayan, J. R.↗

Reflight of the Solid Surface Combustion Experiment: Flame Radiation Near Extinction

In flame spreading in quiescent and low-velocity opposing flows, effects of surface reradiation and flame radiation are important in establishing the spread rate and whether the flame, once ignited, survives to steady spread or extinguishes after a time long compared to the ignition event. A reflight of the Solid Surface Combustion Experiment (SSCE), supported by modelling, demonstrates that for thick, flat fuels, the ultimate fate of the flame is extinction rather than steady spread. A mismatch between the thermal scale in the gas, driven by radiation, and the species diffusion scale, driven by mass diffusion, develops such that the high temperature regions of the flame are ultimately located in a region to which oxygen cannot be supplied at a sufficient rate to sustain reaction, and extinction occurs. Results of the experiment conducted on Space Shuttle mission STS 85 on 9 August 1997 are reviewed. For the flat surface geometry, while the hydrodynamic phenomena associated with opposed-flow flame spread may be treated two dimensionally, the radiative effects are three dimensional, and so modelling the radiative processes, with the mismatch in dimensionality, is difficult. The cylindrical geometry at least one long compared to the radius, provides a configuration in which the radiative processes for spread in the axial direction are two dimensional, thus simplifying the modelling. The cylindrical geometry allows for the development of more sophisticated radiative models without the complication of dimensionality concerns, e.g., discrete transfer, which is discussed in detail by Bundy (1998). Additionally, the cylindrical geometry results in a "focussing" of the heat transfer to the surface and may allow for steady spread for radii that for thick fuels of the same half-thickness there is no steady spread.

Altenkirch, R. A.↗

Transient Evolution of a Planar Diffusion Flame Aft of a Translating Flat Plate

The high degree of spatial symmetry of a planar diffusion flame affords great simplifications for experimental and modeling studies of gaseous fuel combustion. Particularly, in a microgravity environment, where buoyancy effects are negligible, an effectively strain-rate-free, vigorous flame may be obtained. Such a flame can also provide long residence times and large length scales for practical probing of flame structures and soot processes. This 2-D numerical study explores the feasibility of establishing such a planar diffusion flame in an enclosed container utilizing a realistic test protocol for a microgravity experiment. Fuel and oxygen mixtures, initially segregated into two half-volumes of a squat rectangular container by a thin separator, are ignited as soon as a flammable mixture is formed in the wake of the separator withdrawn in the centerplane. A triple-flame ensues that propagates behind the trailing edge of the separator. The results of calculations show that the mechanically- and thermally-induced convection decays in about two seconds. The establishment of a planar diffusion flame after this period seems feasible in the central region of the container with sufficient quantities of reactants left over for subsequent studies. An analysis of the flame initiation and formation process suggests how the feasibility of creating such a flame can be further improved.

Gokoglu, Suleyman A.↗

Droplet Combustion in a Slow Convective Flow

The influences of slow convective flow on droplet combustion, particularly in the low Reynolds number regime, have received very little attention in the past. Most studies in the literature are semi-empirical in nature and they were motivated by spray combustion applications in the moderate to high Reynolds number regime. None of the limited number of fundamental theoretical studies applicable to low Reynolds numbers have been verified by rigorous experimental data. Moreover, many unsteady phenomena associated with fluid-dynamic unsteadiness, such as impulsive starting or stopping of a burning droplet, or flow acceleration/deceleration effects, have not been investigated despite their importance in practical applications. In this study we investigate the effects of slow convection on droplet burning dynamics both experimentally and theoretically. The experimental portion of the study involves both ground-based experiments in the drop towers and future flight experiments on board the International Space Station. Heptane and methanol are used as test fuels, and this choice complements the quiescent-environment studies of the Droplet Combustion Experiment (DCE). An analytical model that employs the method of matched asymptotic expansions and uses the ratio of the convective velocity far from the droplet to the Stefan velocity at its surface as the small parameter for expansion has also been developed as a part of this investigation. Results from the ground-based experiments and comparison with the analytical model are presented in this report.

Nayagam, V.↗

Probing High-Temperature Amine Chemistry: Is the Reaction NH 3 + NH 2 ⇌ N 2 H 3 + H 2 Important?

The reaction NH 3 + NH 2 ⇌ N 2 H 3 + H 2 (R1) has been identified as a key step to explain experimental results for pyrolysis and oxidation of ammonia. However, no direct experimental or theoretical evidence for the reaction has been reported. In the present work, the reaction was studied by ab initio theory and by reinterpretation of experimental data. We could not locate a transition state for R1 occurring as a direct process, but alternative mechanisms yielded an upper bound to k 1 of 1.5 × 10 13 exp(–58.9 kcal mol –1 /RT) cm 3 mol –1 s –1 over 1000–2500 K, several orders of magnitude lower than values applied in modeling. Consistent with the theoretical work, re-evaluation of NH 3 pyrolysis data supported a very low value of k 1 . However, this finding opens up a novel unresolved issue. Current kinetic models cannot capture the NH 3 oxidation behavior in a number of laminar flow reactor and jet-stirred reactor experiments without adopting an improbably high value for k 1 . Furthermore, important oxidation steps might be underestimated or missing from mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Initial Observations on the Burning of an Ethanol Droplet in Microgravity

Combustion of liquid ethanol represents an important system both from fundamental and practical points of view, Ethanol is currently being used as an additive to gasoline in order to reduce carbon monoxide and particulate emissions as well as to improve the fuel octane rating. A detailed physical understanding of liquid ethanol combustion is therefore necessary to achieve an optimal performance of such fuel blends in practical conditions. Ethanol is also a relatively simple model compound suitable for investigation of important combustion characteristics typical of more complex fuels. In particular, ethanol has been proposed for studies of sooting behavior during droplet burning. The sooting nature of ethanol has pressure sensitivities similar to that of n-heptane, but shifted to a higher range of pressures (1-3 atm). Additionally, liquid ethanol is miscible with water produced during its combustion forming mixtures with azeotropic behavior, a phenomenon important for understanding multi-component, liquid fuel combustion. In this work, we present initial results obtained in a series of recent space-based experiments and develop a detailed model describing the burning of ethanol droplet in microgravity.

Kazakov, Andrei↗

Analysis of pressure spectra measurements in a ducted combustion system

Combustion noise propagation in an operating ducted liquid fuel combustion system is studied in relation to the development of combustion noise prediction and suppression techniques. The presence of combustor emissions in the duct is proposed as the primary mechanism producing the attenuation and dispersion of combustion noise propagating in an operating liquid fuel combustion system. First, a complex mathematical model for calculating attenuation and dispersion taking into account mass transfer, heat transfer, and viscosity effects due to the presence of liquid fuel droplets or solid soot particles is discussed. Next, a simpler single parameter model for calculating pressure auto-spectra and cross-spectra which takes into account dispersion and attenuation due to heat transfer between solid soot particles and air is developed. Then, auto-spectra and cross-spectra obtained from internal pressure measurements in a combustion system consisting of a J-47 combustor can, a spool piece, and a long duct are presented. Last, analytical results obtained with the single parameter model are compared with the experimental measurements. The single parameter model results are shown to be in excellent agreement with the measurements.

Miles, J. H.↗

Propulsive Reaction Control System Model

This software models a propulsive reaction control system (RCS) for guidance, navigation, and control simulation purposes. The model includes the drive electronics, the electromechanical valve dynamics, the combustion dynamics, and thrust. This innovation follows the Mars Science Laboratory entry reaction control system design, and has been created to meet the Mars Science Laboratory (MSL) entry, descent, and landing simulation needs. It has been built to be plug-and-play on multiple MSL testbeds [analysis, Monte Carlo, flight software development, hardware-in-the-loop, and ATLO (assembly, test and launch operations) testbeds]. This RCS model is a C language program. It contains two main functions: the RCS electronics model function that models the RCS FPGA (field-programmable-gate-array) processing and commanding of the RCS valve, and the RCS dynamic model function that models the valve and combustion dynamics. In addition, this software provides support functions to initialize the model states, set parameters, access model telemetry, and access calculated thruster forces.

Brugarolas, Paul↗

Combustion of LOX with H2(sub g) under subcritical, critical, and supercritical conditions and experimental observation of dense spray and mixing of impinging jets

Progress made during the period of February 1 to October 15, 1992 is reported. The overall objective of Task 1 of the investigation is to achieve a better understanding of the combustion processes of liquid oxygen and gaseous hydrogen under subcritical, critical, and supercritical conditions. Specific objectives of the research program are: (1) to determine the evaporation- and burning-rate characteristics of LOX in hydrogen/helium environments under broad ranges of operating conditions; (2) to measure species concentration profile and surface temperature of LOX employing the gas chromatography and fine-wire thermocouples under non-reacting flow situations; (3) to perform a fugacity-based multicomponent thermodynamic phase equilibrium analysis for examining the high-pressure vapor-liquid equilibrium behavior at the liquid surface of LOX; (4) to formulate and solve a theoretical model for simulating the evaporation and combustion processes in a LOX/H2/He system; and (5) to validate the theoretical model with the measured experimental data. Task 2 of the investigation is described. Observation of a like-on-like injector element in the near-injector region performed in the previous phase of this project has identified the existence of a high Reynolds number regime in which the pre-impingement jets are fully turbulent and undergoing surface breakup. The new spray regime, which has not been observed by previous investigators, is characterized by the presence of many fine droplets and the disappearance of the well-defined liquid breakup wave pattern in the post-impingement region. It is speculated that a cavitating region may be present within the orifice so that it could induce strong turbulence, leading to an onset of atomization of the jets prior to impingement. To further investigate the dense spray behavior of the impinging jets in the high Reynolds number region, experiments were conducted using Plexiglas injector components for direct internal flow observation. The main objective is to determine under what conditions a cavitating region would form and whether or not the cavitation is reponsible for the development of the high Reynolds number spray regime. The procedure and major findings of the injector cavitation study are described.

Kuo, Kenneth K.↗

Numerical analysis of soot emissions from gasoline-ethanol and gasoline-butanol blends under gasoline compression ignition conditions

In the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine were performed to investigate the impact of blending two biofuels, ethanol and n-butanol, with gasoline on the trade-off between combustion phasing and soot emissions under low load conditions. Here, in order to represent market gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF) + ethanol (ETPRF) surrogate (with 20% ethanol by mole; E20) was formulated using a neural network based octane predictor such that the surrogate had the same ethanol content, Research Octane Number (RON) and Octane Sensitivity (S). In addition, a novel skeletal kinetic mechanism for ETPRF and TPRF + n-butanol (BTPRF) blends, incorporating polycyclic aromatic hydrocarbon (PAH) chemistry, was developed. A three-dimensional (3D) engine CFD formulation employing the skeletal mechanism, adaptive mesh refinement (AMR), finite-rate chemistry approach, and hybrid method of moments (HMOM) was adopted to capture the in-cylinder combustion phenomena and soot emissions. The engine CFD model was validated against RD5-87 experimental data for a broad range of start-of-injection (SOI) timings (-21/-27/-36/-45 crank angle degrees (CAD) after top-dead center (aTDC)), with respect to in-cylinder pressure, heat release rate, combustion phasing, and soot emissions. The closed-cycle simulation results were analyzed to elucidate the non-monotonic trend of soot emissions versus SOI timing: SOI-36 > SOI-45 > SOI-21 > SOI-27. Thereafter, the validated CFD model was employed to simulate the combustion of a gasoline-ethanol blend with 45% (by mole) ethanol (E45) and a gasoline-butanol blend with 45% (by mole) n-butanol (B45) under the same operating conditions to study the effects of fuel composition and SOI timing on combustion phasing and soot emissions. The sooting propensity followed the trend: B45 > E20 > E45 at all SOI timings. Overall, it was observed that the autoignition propensity was primarily related to fuel chemistry. On the other hand, sooting propensity showed strong coupling with both fuel chemistry and physical properties, with greater impact of fuel physical properties at advanced SOI timings.

30 DIRECT ENERGY CONVERSION↗

Investigation of Thermal Radiation under Pressurized Oxy-combustion Conditions

Thermal radiation of the gaseous and particle phases in a pilot-scale pressurized oxy-combustor is computationally studied. In particular, the radiation characteristics of gases and particles are estimated by employing the statistical narrow-band model and the large-particle model. It is found that thermal radiation of the particle cloud dominates in the combustor under a furnace temperature of 1500 K and when there is no substantial loss of particles to the walls. Another important observation is that radiation from the gas and particles can be approximately treated as a graybody under these conditions. More specifically, the results on the spectral radiation intensity of a gas comprising 40% (vol) H 2 O and 60% CO 2 show that when the pressure is increased to 15 bar, and the radiation pathlength is 100 cm, the spectral radiation profile of the gas phase approaches that of a blackbody at the respective temperature. In addition, the emissivity of the particulate cloud has been evaluated as a function of the particle concentration and diameter by employing the large-particle model. It is shown that the emissivity grows with the particle concentration but decreases with the particle size for the same mass of the particles. Finally, this outcome of the present study is expected to be used to validate the assumption of the gray-gas model adopted in the numerical simulations of pressurized oxy-combustion.

large-particle model↗

Technique for predicting high-frequency stability characteristics of gaseous-propellant combustors

A technique for predicting the stability characteristics of a gaseous-propellant rocket combustion system is developed based on a model that assumes coupling between the flow through the injector and the oscillating chamber pressure. The theoretical model uses a lumped parameter approach for the flow elements in the injection system plus wave dynamics in the combustion chamber. The injector flow oscillations are coupled to the chamber pressure oscillations with a delay time. Frequency and decay (or growth) rates are calculated for various combustor design and operating parameters to demonstrate the influence of various parameters on stability. Changes in oxidizer design parameters had a much larger influence on stability than a similar change in fuel parameters. A complete description of the computer program used to make these calculations is given in an appendix.

Priem, R. J.↗

Measurements of reactive gaseous rocket injector response factors

The results presented represent the first successful attempt at the measurement of the driving capabilities of coaxial gaseous propellant rocket injectors. The required data have been obtained by employing the modified impedance tube technique with compressed air as the oxidizer and acetylene gas as the fuel. The data describe the frequency dependence of the injector admittances, from which the frequency dependence of the injection response factors can be calculated. The measured injector admittances have been compared with the predictions of the Feiler and Heidmann (1967) analytical model assuming different values for the characteristic combustion time. The values of combustion time which result in a best fit between the measured and predicted data are indicated for different equivalence ratios. It is shown that for the coaxial injector system investigated in this study the characteristic combustion times vary between .7 and 1.2 msec for equivalence ratios in the range of .57 to 1.31. The experimental data clearly show that the tested injector system could indeed drive combustion instabilities over a frequency range that is in qualitative agreement with the predictions of the Feiler and Heidmann model.

Janardan, B. A.↗

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A mathematical model for jet engine combustor pollutant emissions

Mathematical modeling for the description of the origin and disposition of combustion-generated pollutants in gas turbines is presented. A unified model in modular form is proposed which includes kinetics, recirculation, turbulent mixing, multiphase flow effects, swirl and secondary air injection. Subelements of the overall model were applied to data relevant to laboratory reactors and practical combustor configurations. Comparisons between the theory and available data show excellent agreement for basic CO/H2/Air chemical systems. For hydrocarbons the trends are predicted well including higher-than-equilibrium NO levels within the fuel rich regime. Although the need for improved accuracy in fuel rich combustion is indicated, comparisons with actual jet engine data in terms of the effect of combustor-inlet temperature is excellent. In addition, excellent agreement with data is obtained regarding reduced NO emissions with water droplet and steam injection.

Boccio, J. L.↗

Experimental and modeling study of water time histories during H 2 S-N 2 O combustion in a shock tube

Here, with the interest in directly burning sour gas in gas turbines, and the fact that even small amounts of H 2 S or its combustion products can alter combustion characteristics, many research studies have been performed to better understand the combustion chemistry of H 2 S. In the present study, the water formation was followed by laser absorption with N 2 O as an oxidant, instead of O 2 . Nitrous Oxide being essentially consumed via N 2 O (+M) ⇌ N 2 + O (+M), the water formation via the H 2 S + O route can then be probed to further validate the models. Three H 2 S/N 2 O mixtures diluted in 98% Ar were studied to cover the following range of equivalence ratios: 0.5, 1.0, and 2.0, over a wide range of temperatures (1580–1940 K) around atmospheric pressure. A chemical kinetic model was then developed first by validating the base N 2 O kinetics mechanism, then by investigating the experimental results presented herein. The N 2 O kinetics mechanism was updated with recent work on the low- and high-pressure limits for N 2 O decomposition (N 2 O (+M) ⇌ N 2 + O (+M)) as well as a review of rate coefficients for N 2 O + H ⇌ N 2 + OH from the literature. Good agreement is shown for H 2 /N 2 O mixtures. Updates were then made to the H 2 S kinetics mechanism, specifically, an update from the literature on SO 2 + H (+M) ⇌ SO + OH (+M) and an adjustment to SO + SH ⇌ S 2 + OH. Additionally, reactions between SH and N 2 O were determined using W1BD data and transition state theory which required the addition of an NNS sub-mechanism. With these updates, the mechanism provides good agreement with the H 2 S/N 2 O experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Aldehyde cool-flame chemistry explains a missing source of organic acids

Combustion emission is a significant source of organic acids, impacting atmospheric chemistry and climate. Their formation mechanisms, however, remain poorly understood, leading to underestimation in kinetic models. We investigate the cool-flame oxidation of key combustion intermediates—C 1 –C 4 aldehydes and benzaldehyde. Using in-situ synchrotron vacuum ultraviolet photoionization mass spectrometry, we observe the direct conversion of aldehydes to organic acids, a process enhanced by HO 2 radicals. Quantum chemistry calculations reveal that the reaction of RC(O)O 2 with HO 2 on the singlet potential energy surface contributes to organic acids. Incorporating this pathway into a kinetic model significantly improves organic acid prediction. Despite the high-temperature nature of engine combustion, significant spatial and temporal inhomogeneities (e.g., near-wall regions and crevice volumes) lead to localized cool-flame conditions, facilitating organic acid formation and emission. Elucidating the acid formation under cool-flame conditions provides a critical mechanism for accurately modelling anthropogenic organic acid emissions and developing mitigation strategies.

SVUV-PIMS↗