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,621 records · Page 90

Global Gas-Phase Oxidation Rates of Select Products from the Fast Pyrolysis of Lignocellulose

The oxidation kinetics for products of fast pyrolysis at low temperatures (<600°C) are not well known. These will be important in effort to model autothermal pyrolysis, which has been recently developed to intensify the process, but which occurs at much lower temperatures than combustion. Furthermore, this study determines global oxidation rates at 400-600°C for three important products of fast pyrolysis: levoglucosan, xylose, and acetic acid. Experiments were performed in a fluidized bed pyrolyzer with the reactor modeled as a series of CSTRs and PFRs to determine reaction rates. Oxidation rates at 500°C for the three model compounds varied by a factor of ten.

09 BIOMASS FUELS↗

Quantifying the Contributions of Aerosol- and Snow-Produced ClNO 2 through Observations and 1D Modeling

Nitryl chloride (ClNO 2 ) is a radical reservoir that forms and accumulates in the nocturnal atmospheric boundary layer influenced by combustion emissions and chloride (e.g., sea salt and road salt). Upon sunrise, ClNO 2 rapidly photolyzes to generate highly reactive chlorine radicals (Cl • ) that affect the air quality by generating secondary air pollutants. Recent studies have shown road salt aerosols and saline snowpack to be sources of ClNO 2 in the wintertime urban environment; however, the quantitative contributions of each chloride source are not known. In this study, we examine the vertically resolved contributions of aerosol particles and saline snowpack as sources of ClNO 2 by using an observationally constrained snow–atmosphere coupled one-dimensional model applied to wintertime Kalamazoo, Michigan, U.S. Model simulations show that ClNO 2 emitted from urban snowpack can be vertically transported throughout the entire atmospheric boundary layer and can be a significant source of ClNO 2 , contributing up to ~60% of the ClNO 2 budget near the surface. Here, modeled snowpack ClNO 2 emission rates were 6 (±7) times higher than the observationally derived emission rates, suggesting that not all snow chloride is available for reaction. ClNO 2 production from both aerosol particles and snow emissions are required to best simulate the observed surface-level ClNO 2 . Using the traditional bulk parameterization for ClNO 2 produced from particles significantly overestimated ClNO 2 due to the assumption of having equivalent dinitrogen pentoxide (N 2 O 5 ) uptake and chloride availability for the entire particle population. In comparison, the chemically resolved surface area-based parameterization slightly underestimated the observations and had uncertainties deriving from ClNO 2 production from residential wood burning particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Biomass Combustion in a Circulating Fluidized Bed Combustor

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has skyrocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing resulting in higher combustion efficiency and the ability to use lower quality fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions. SO 2 emissions are conveniently controlled through the addition of Ca-based sulphur sorbents within the combustor. This report summarizes the current modeling effort on a 50 kW th CFB combustor with a diameter of 10 cm and a height of 5 m; designed, built, and operated at CanmetENERGY in Ottawa, Canada employing the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The MFiX-PIC model parameters for the simulation are tuned against cold-flow experiments from CanmetENERGY using 9 kg of olivine sand as the inert bed material. It is shown that for the relatively coarse fluid meshes and large parcel sizes necessitated by the scale of the simulation, filter size dependent corrections to the drag law must be incorporated to ensure accuracy of the simulation results. The validated cold flow model is extended to simulate reacting flow with torrefied hardwood as the feedstock and to validate the combustion reaction scheme. The species concentrations at the riser outlet are compared against CanmetENERGY’s experiments and show satisfactory agreement. The simulations demonstrate the ability of MFiX-PIC to accurately capture the physics and chemistry of a circulating fluidized bed combustor at bench scales, which can be further extended to pilot- and industrial-scale systems.

20 FOSSIL-FUELED POWER PLANTS↗

Demonstration of a Fiber Optic Regression Probe

The capability to provide localized, real-time monitoring of material regression rates in various applications has the potential to provide a new stream of data for development testing of various components and systems, as well as serving as a monitoring tool in flight applications. These applications include, but are not limited to, the regression of a combusting solid fuel surface, the ablation of the throat in a chemical rocket or the heat shield of an aeroshell, and the monitoring of erosion in long-life plasma thrusters. The rate of regression in the first application is very fast, while the second and third are increasingly slower. A recent fundamental sensor development effort has led to a novel regression, erosion, and ablation sensor technology (REAST). The REAST sensor allows for measurement of real-time surface erosion rates at a discrete surface location. The sensor is optical, using two different, co-located fiber-optics to perform the regression measurement. The disparate optical transmission properties of the two fiber-optics makes it possible to measure the regression rate by monitoring the relative light attenuation through the fibers. As the fibers regress along with the parent material in which they are embedded, the relative light intensities through the two fibers changes, providing a measure of the regression rate. The optical nature of the system makes it relatively easy to use in a variety of harsh, high temperature environments, and it is also unaffected by the presence of electric and magnetic fields. In addition, the sensor could be used to perform optical spectroscopy on the light emitted by a process and collected by fibers, giving localized measurements of various properties. The capability to perform an in-situ measurement of material regression rates is useful in addressing a variety of physical issues in various applications. An in-situ measurement allows for real-time data regarding the erosion rates, providing a quick method for empirically anchoring any analysis geared towards lifetime qualification. Erosion rate data over an operating envelope could also be useful in the modeling detailed physical processes. The sensor has been embedded in many regressing media for the purposes of proof-of-concept testing. A gross demonstration of its capabilities was performed using a sanding wheel to remove layers of metal. A longer-term demonstration measurement involved the placement of the sensor in a brake pad, monitoring the removal of pad material associated with the normal wear-and-tear of driving. It was used to measure the regression rates of the combustable media in small model rocket motors and road flares. Finally, a test was performed using a sand blaster to remove small amounts of material at a time. This test was aimed at demonstrating the unit's present resolution, and is compared with laser profilometry data obtained simultaneously. At the lowest resolution levels, this unit should be useful in locally quantifying the erosion rates of the channel walls in plasma thrusters. .

Korman, Valentin↗

Analytical combustion/emissions research related to the NASA high-speed research program

Increasing the pressure and temperature of the engines of new generation supersonic airliners increases the emissions of nitrogen oxides to a level that would have an adverse impact on the Earth's protective ozone layer. In the process of implementing low emissions combustor technologies, NASA Lewis Research Center has pursued a combustion analysis program to guide combustor design processes, to identify potential concepts of greatest promise, and to optimize them at low cost, with short turn-around time. The approach is to upgrade and apply advanced computer programs for gas turbine applications. Efforts have been made to improve the code capabilities of modeling the physics. Test cases and experiments are used for code validation. To provide insight into the combustion process and combustor design, two-dimensional and three-dimensional codes such as KIVA-II and LeRC 3D have been used. These codes are operational and calculations have been performed to guide low emissions combustion experiments.

Nguyen, H. Lee↗

In-stream measurements of combustion during Mach 5 to 7 tests of the Hypersonic Research Engine (HRE)

Results of in-stream combustion measurements taken during Mach 5 to 7 true simulation testing of the Hypersonic Research Engine/Aerothermodynamic Integration Model (HRE/AIM) are presented. These results, the instrumentation techniques, and configuration changes to the engine installation that were required to test this model are described. In test runs at facility Mach numbers of 5 to 7, an exhaust instrumentation ring which formed an extension of the engine exhaust nozzle shroud provided diagnostic measurements at 10 circumferential locations in the HRE combustor exit plane. The measurements included static and pitot pressures using conventional conical probes, combustion gas temperatures from cooled-gas pyrometer probes, and species concentration from analysis of combustion gas samples. Results showed considerable circumferential variation, indicating that efficiency losses were due to nonuniform fuel distribution or incomplete mixing. Results using the Mach 7 facility nozzle but with Mach 6 temperature simulation, 1590 to 1670 K, showed indications of incomplete combustion. Nitric oxide measurements at the combustor exit peaked at 2000 ppmv for stoichiometric combustion at Mach 6.

Lezberg, Erwin A.↗

Study the thermal radiation effects in gas-solid flows with gray and non-gray P1 models implemented in MFiX

Thermal radiation is a dominant mode of heat transfer in combustion/gasification, packed/circulating bed reactors, and energy storage/conversion devices. In hydrocarbon fuel combustion, the high absorption and emission of product gases like CO 2 and H 2 O significantly affect the heat transfer characteristics. It is essential to understand the radiative energy propagation in such applications. In the present study, we developed a framework to implement the P-1 radiation model in MFiX for Eulerian-Eulerian (Two fluid-TFM) and Eulerian-Lagrangian (Discrete element-DEM and Particle in cell-PIC) models. The P-1 is a simplification of the spherical harmonics method. The framework has been verified on various flue gas conditions. The verified solver results are validated against the benchmark results available in the literature. The verified and validated gray and non-gray weighted sum of gray gas models is evaluated for dry and wet flue gas conditions. Further, they are also used to study gas-solid flows in a fluidized bed/riser.

20 FOSSIL-FUELED POWER PLANTS↗

Transient simulation of biomass combustion in a circulating fluidized bed riser

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has sky-rocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing, which results in higher combustion efficiency and the ability to use lower rank fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions, while SO 2 emissions can be conveniently controlled through the addition of Ca-based sulfur sorbents within the combustor. Herin this paper summarizes the modeling effort on a 50 kW th CFB combustor designed, built, and operated at CanmetENERGY in Ottawa, Canada. The numerical model employs the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The MFiX-PIC model parameters for the simulation are tuned against cold-flow experiments from CanmetENERGY using olivine sand as the inert bed material. It is shown that for the relatively coarse fluid meshes and large parcel sizes necessitated by the scale of the simulation, filter size dependent corrections to the drag law must be incorporated to ensure accuracy of the simulation results. The validated cold flow model is extended to simulate reacting flow with torrefied hardwood as the feedstock and to validate the combustion reaction scheme. The species concentrations at the riser outlet are compared against CanmetENERGY’s experiments and show satisfactory agreement. The simulations demonstrate the ability of MFiX-PIC to accurately capture both the physics and chemistry of a CFB combustor at bench scales, which can be further extended to pilot- and industrial-scale systems.

09 BIOMASS FUELS↗

SYSTEMATIC MULTI-INDEX VALIDATIONS OF SIDI ENGINE FLOW FIELD LES COMPUTATIONS USING CRANK ANGLE-RESOLVED PIV MEASUREMENTS

In-cylinder flow fields make significant impacts on the fuel atomization, fuel mixture formation, and combustion process in spark ignition direct injection (SIDI) engines. In recent years, model-based simulation approaches are preferred in regard to investigating the transient in-cylinder flow field characteristics. Most commonly, the simulation models are validated using single representative flow field at a typical crank angle measured by particle image velocimetry (PIV). However, it provides only limited knowledge about the flow field which is highly three-dimensional and of transient nature. In this study, crank angle-resolved Ply measurements are conducted on three distinct planes inside the cylinder to capture the transient process of flow field characteristics which vary with the crank angle. These three planes consist of one tumble plane through the spark plug tip, one tumble plane along two intake ports, and one swirl plane at 30 mm below the cylinder head. Large eddy simulation (LES) is employed for the numerical computations using the CONVERGE codes. On the basis of large datasets for both temporal and spatial domains, a multi-index systematic validation approach is conducted. Crank angle-resolved calculations of global indices and local indices are implemented using the flow fields velocity data obtained from both PIV and LES on select planes. Global indices reveal the trends in similarities of different crank angle degrees and locations, while local indices give the detail comparison results. In summary, with the systematic multi-index validation approach, the crucial crank angle degrees and locations for model verification will be detected. Furthermore, the corresponding critical flow features are analyzed. Practical guideline of flow field validation is proposed.

Crank angle-resolved PIV↗

The Deflagration-to-Detonation Transition in Two Dimensions

The Deflagration-to-Detonation Transition (DDT) in one-dimensional porous explosive, where combustion in an explosive transitions to detonation, can be described by the following model. This simplified model proceeds in five steps, as follows: 1) Ignition of the explosive, surface burning. 2) Convective burning, with the flame front penetrating through the porous network of the explosive. This proceeds until the pressure grows high enough to result in choked flow in the pores restricting the convective burn. 3) The choked flow results in the formation of a high-density compact of explosive. This compact is driven into undisturbed material by the pressure of the burning explosive. 4) The compression of the undisturbed porous explosive by the compact leads to the ignition of a compressive burn. This builds in pressure until a supported shock forms. 5) The shock builds in pressure until detonation occurs. See Figures 1 and 2 for an overview of the proceeding steps. It has been assumed in the past that the same mechanism which drives DDT in the one dimensional case will apply to both two and three dimensions. However, this has never been experimentally verified and it is possible that an entirely new mechanism drives DDT higher dimensions. In order to test this a series of two dimensional DDT tests were performed to help provide empirical evidence of the mechanism in higher dimensions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Sensitivity Studies of Premixed Flame Response to Transverse, High Frequency Disturbances

This study analyzes the response of premixed flames to high frequency transverse acoustic disturbances. A significant body of experiment, modeling and simulation efforts has been performed to understand the response of the spatially integrated heat release of premixed flames to low-frequency disturbances. These studies quantified this sensitivity through the Flame Transfer Function (FTF), as the global heat release is the key quantity influencing thermoacoustic stability for acoustically compact flames. Recent efforts have expanded these models to consider high-frequency acoustic disturbances, and the sensitivity of premixed flames through the Rayleigh integral (RI). The full Rayleigh Integral, as opposed to simply the spatially integrated heat release, must be considered in these cases where the flame is not acoustically compact. Prior modeling studies considered several effects - stability of axisymmetric flames centered in a cylindrical combustion chamber and the effect of flame offset from the combustor axis and the effect of mean flame asymmetries in the flame centered in the combustor. These studies showed that the flame's stability was sensitive to all these parameters, as well as to the transverse mode in question. In this paper, we quantify these different parametric sensitivities of premixed flame stability in the presence of high-frequency transverse acoustic modes - specifically sensitivity to flame length and nozzle location within the acoustic mode shape for 2 different transverse modes. The sensitivity of the RI to these control parameters provides an improved understanding of which flame parameters are critical to improve the overall thermoacoustic stability of the system.

Acharya, Vishal S.↗

Further analytical study of hybrid rocket combustion

Analytical studies of the transient and steady-state combustion processes in a hybrid rocket system are discussed. The particular system chosen consists of a gaseous oxidizer flowing within a tube of solid fuel, resulting in a heterogeneous combustion. Finite rate chemical kinetics with appropriate reaction mechanisms were incorporated in the model. A temperature dependent Arrhenius type fuel surface regression rate equation was chosen for the current study. The governing mathematical equations employed for the reacting gas phase and for the solid phase are the general, two-dimensional, time-dependent conservation equations in a cylindrical coordinate system. Keeping the simplifying assumptions to a minimum, these basic equations were programmed for numerical computation, using two implicit finite-difference schemes, the Lax-Wendroff scheme for the gas phase, and, the Crank-Nicolson scheme for the solid phase.

Hung, W. S. Y.↗

Predicted exhaust emissions from a methanol and jet fueled gas turbine combustor

A computer model of a gas turbine combustor has been used to predict the kinetic combustion and pollutant formation processes for methanol and simulated jet fuel. Use of the kinetic reaction mechanisms has also allowed a study of ignition delay and flammability limit of these two fuels. The NOX emissions for methanol were predicted to be from 69 to 92% lower than those for jet fuel at the same equivalence ratio which is in agreement with experimentally observed results. The high heat of vaporization of methanol lowers both the combustor inlet mixture temperatures and the final combustion temperatures. The lower combustion temperatures lead to low NOX emissions while the lower inlet mixture temperatures increase methanol's ignition delay. This increase in ignition delay dictates the lean flammability limit of methanol to be 0.8, while jet fuel is shown to combust at 0.4.

Adelman, H. G.↗

Numerical optimization of gasdynamic lasers

A sophisticated numerical optimization procedure has been coupled to a detailed theoretical model of the carbon dioxide gasdynamic laser (GDL) for the purpose of systematically evaluating the performance of these systems. The integrated model was used to predict the performance of a wide variety of combustion-driven gasdynamic lasers and was then used to maximize the theoretical available specific energy of these lasers. Up to six design variables and three nonlinear constraints were used in these optimizations, although the method is not limited to this number. Results are presented to show the performance of typical lasers before and after optimization. Improvements approaching a factor of two in the available specific energy are found, suggesting that numerical optimization techniques are valuable in improving the design of gasdynamic lasers.

Saunders, R. C., III↗

A comprehensive model to determine the effects of temperature and species fluctuations on reaction rates in turbulent reacting flows

A computationally-viable model describing the interaction between fluid-mechanical turbulence and finite-rate combustion reactions, principally in high-speed flows was developed. Chemical kinetic mechanisms, complete and global, were developed describing the finite rate reaction of fuels of interest to NASA. These fuels included principally hydrogen and silane, although a limited amount of work involved hydrocarbon fuels as well.

Chinitz, W.↗

Radiation Effects on Flow Characteristics in Combustion Chambers

A JANNAF sponsored workshop was held to discuss the importance and role of radiative heat transfer in rocket combustion chambers. The potential impact of radiative transfer on hardware design, reliability, and performance was discussed. The current state of radiative transfer prediction capability in CFD modeling was reviewed and concluded to be substantially lacking in both the physical models used and the radiative property data available. There is a clear need to begin to establish a data base for making radiation calculations in rocket combustion chambers. A natural starting point for this effort would be the NASA thermochemical equilibrium code (CEC).

Brewster, M. Q.↗

Free-jet Tests of a 1.1-inch-diameter Supersonic Ram-jet Engine

Results are presented of free-jet tests of a 1.1-inch-diameter hydrogen-burning ram-jet engine over a Mach number range from 1.42 to 2.28 and a Reynolds number range from 6.01 times 10 to the 6th power to 15.78 times 10 to the 6th power. Tests demonstrated the reliability and wide operating range of the engine and showed its suitability for installation on wind-tunnel models of airplane and missile configurations. A comparison of engine operation with combustion-chamber lengths of 3.33 and 1.51 engine diameters was made at a Mach number of 2.06. A maximum test thrust coefficient of 0.905 was obtained at fuel-air ratio of 0.034 and a Mach number of 2.06 for the engine with the 3.33-diameter combustion chamber.

ENGINES, RAM-JET↗

Computer program for calculation of complex chemical equilibrium compositions and applications. Part 1: Analysis

This report presents the latest in a number of versions of chemical equilibrium and applications programs developed at the NASA Lewis Research Center over more than 40 years. These programs have changed over the years to include additional features and improved calculation techniques and to take advantage of constantly improving computer capabilities. The minimization-of-free-energy approach to chemical equilibrium calculations has been used in all versions of the program since 1967. The two principal purposes of this report are presented in two parts. The first purpose, which is accomplished here in part 1, is to present in detail a number of topics of general interest in complex equilibrium calculations. These topics include mathematical analyses and techniques for obtaining chemical equilibrium; formulas for obtaining thermodynamic and transport mixture properties and thermodynamic derivatives; criteria for inclusion of condensed phases; calculations at a triple point; inclusion of ionized species; and various applications, such as constant-pressure or constant-volume combustion, rocket performance based on either a finite- or infinite-chamber-area model, shock wave calculations, and Chapman-Jouguet detonations. The second purpose of this report, to facilitate the use of the computer code, is accomplished in part 2, entitled 'Users Manual and Program Description'. Various aspects of the computer code are discussed, and a number of examples are given to illustrate its versatility.

Gordon, Sanford↗