0
Research Papers: Gas Turbines: Combustion, Fuels, and Emissions

Large-Eddy Simulation of Turbulent Spray Combustion in a Subscale Aircraft Jet Engine Combustor—Predictions of NO and Soot Concentrations

[+] Author and Article Information
Hideki Moriai

Department of Mechanical
Engineering and Science,
Kyoto University,
Kyoto 615-8540, Japan;
Aerospace Systems,
Mitsubishi Heavy Industries, Ltd.,
Aichi 485-8561, Japan

Ryoichi Kurose

Department of Mechanical
Engineering and Science,
Kyoto University,
Kyoto 615-8540, Japan
e-mail: kurose@mech.kyoto-u.ac.jp

Hiroaki Watanabe

Energy Engineering Research Laboratory,
Central Research Institute of Electric
Power Industry (CRIEPI),
Kanagawa 240-0196, Japan

Yutaka Yano, Satoru Komori

Department of Mechanical
Engineering and Science,
Kyoto University,
Kyoto 615-8540, Japan

Fumiteru Akamatsu

Department of Mechanical Engineering,
Osaka University,
Osaka 565-0871, Japan

Contributed by the Combustion and Fuels Committee of ASME for publication in the Journal of Engineering for Gas Turbines and Power. Manuscript received February 6, 2013; final manuscript received May 2, 2013; published online July 31, 2013. Assoc. Editor: Song-Charng Kong.

J. Eng. Gas Turbines Power 135(9), 091503 (Jul 31, 2013) (10 pages) Paper No: GTP-13-1043; doi: 10.1115/1.4024868 History: Received February 06, 2013; Revised May 02, 2013

Large-eddy simulation (LES) is applied to turbulent spray combustion fields in a subscale (1/2) aircraft jet engine combustor with an air-blast type swirl fuel nozzle and validity is examined by comparing with measurements. In the LES, Jet-A is used as liquid fuel, and individual droplet motion is tracked in a Lagrangian manner with a parcel model. As a turbulent combustion model, the extended flamelet/progress-variable approach, in which heat transfer between droplets and ambient gas including radiation and heat loss from walls can be taken into account, is employed. A detailed chemistry mechanism of Jet-A with 1537 reactions and 274 chemical species is used. The radiative heat transfer is computed by the discrete ordinate (DO) method. The equivalence ratio ranges from 0.91 to 1.29. The comparisons of the predicted droplet velocity and size, gaseous temperature, NO, and soot emissions with the measurements show that the present LES is capable of capturing the general features of the turbulent spray combustion fields in the subscale (1/2) aircraft jet engine combustor.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.

References

Moin, P., and Apte, S. V., 2006, “Large-Eddy Simulation of Realistic Gas Turbine Combustors,” AIAA J., 44, pp. 698–708. [CrossRef]
Boileau, M., Pascaud, S., Riber, E., Cuenot, B., Gicquel, L. Y. M., PoinsotT. J., and CazalensM., 2008, “Investigation of Two-Fluid Methods for Large Eddy Simulation of Spray Combustion in Gas Turbines,” Flow, Turbul. Combust., 80, pp. 291–321. [CrossRef]
Patel, N., and Menon, S., 2008, “Simulation of Spray-Turbulence-Flame Interactions in a Lean Direct Injection Combustor,” Combust. Flame, 153, pp. 228–257. [CrossRef]
Ihme, M., and Pitsch, H., 2008, “Modeling of Radiation and Nitric Oxide Formation in Turbulent Nonpremixed Flames Using a Flamelet/Progress Variable Formulation,” Phys. Fluids, 20, p. 055110. [CrossRef]
Fujita, A., Watanabe, H., Kurose, R., and Komori, S., 2013, “Two-Dimensional Direct Numerical Simulation of Spray Flames. Part 1: Effects of Equivalence Ratio, Fuel Droplet Size and Radiation, and Validity of Flamelet Model,” Fuel, 104, pp. 515–525. [CrossRef]
Kitano, T., Nakatani, T., Kurose, R., and Komori, S., 2013, “Two-Dimensional Direct Numerical Simulation of Spray Flames. Part 2: Effects of Ambient Pressure and Lift, and Validity of Flamelet Model,” Fuel, 104, pp. 526–535. [CrossRef]
Pitsch, H., 2009, private communication.
Watanabe, H., Kurose, R., Komori, S., and Pitsch, H., 2008, “Effects of Radiation on Spray Flame Characteristics and Soot Formation,” Combust. Flame, 152, pp. 2–13. [CrossRef]
Fiveland, W. A., 1988, “Three-Dimensional Radiative Heat Transfer Solutions by the Discrete-Ordinates Method,” J. Thermophys., 2, pp. 309–316. [CrossRef]
Kurose, R., Makino, H., Komori, S., Nakamura, M., Akamatsu, F., and Katsuki, M., 2003, “Effects of Outflow From the Surface of a Sphere on Drag, Shear Lift, and Scalar Diffusion,” Phys. Fluids, 15, pp. 2338–2351. [CrossRef]
Nakamura, M., Akamatsu, F., Kurose, R., and Katsuki, M., 2005, “Combustion Mechanism of Liquid Fuel Spray in a Gaseous Flame,” Phys. Fluids, 17, p. 123301. [CrossRef]
Moin, P., Squires, K., Cabot, W., and Lee, S., 1991, “A Dynamic Subgrid-Scale Model for Compressible Turbulence and Scalar Transport,” Phys. Fluids, A3, pp. 2746–2757. [CrossRef]
Pierce, C. D., and Moin, P., 2004, “Progress-Variable Approach for Large-Eddy Simulation of Non-Premixed Turbulent Combustion,” J. Fluid Mech., 504, pp. 73–97. [CrossRef]
Pitsch, H., 1998, A C++ Computer Program for 0-D Combustion and 1-D Laminar Flame Calculation, RWTH Aachen, Aachen, Germany.
Watanabe, H., Kurose, R., Hwang, S. M., and AkamatsuF., 2007, “Characteristics of Flamelets in Spray Flames Formed in a Laminar Counterflow,” Combust. Flame, 148, pp. 234–248. [CrossRef]
Baba, Y., and Kurose, R., 2008, “Analysis and Flamelet Modelling for Spray Combustion,” J. Fluid Mech., 612, pp. 45–79. [CrossRef]
Miller, R. S., Harstad, K., and Bellan, J., 1998, “Evaluation of Equilibrium and Non-Equilibrium Evaporation Models for Many-Droplet Gas-Liquid Flow Simulations,” Int. J. Multiphase Flow, 24, pp. 1025–1055. [CrossRef]
Miller, R. S., and Bellan, J., 1999, “Direct Numerical Simulation of a Confined Three-Dimensional Gas Mixing Layer With One Evaporating Hydrocarbon-Droplet-Laden Stream,” J. Fluid Mech., 384, pp. 293–338. [CrossRef]
Kuo, K. K. Y., 1986, Principles of Combustion, John Wiley and Sons, New York.
Hayashi, J., Watanabe, H., Kurose, R., and Akamatsu, F., 2011, “Effects of Fuel Droplet Size on Soot Formation in Spray Flames Formed in a Laminar Counterflow,” Combust. Flame, 158, pp. 2559–2568. [CrossRef]
Watanabe, H., Tanno, K., Baba, Y., Kurose, R., Komori, S., 2009, “Large-Eddy Simulation of Coal Combustion in a Pulverized Coal Combustion Furnace With a Complex Burner. Turbulence,” Heat and Mass Transfer 6, K.Hanjalić, Y.Nagano, and S.Jakirlić, eds., Begell House Inc., New York, pp. 1027–1030.
Kurose, R., Anami, M., Fujita, A., and Komori, S., 2012 “Numerical Simulation of Flow Past a Heated/Cooled Sphere,” J. Fluid Mech., 692, pp. 332–346. [CrossRef]
Muto, M., Tsubokura, M., and Oshima, N., 2012, “Negative Magnus Lift on a Rotating Sphere at Around the Critical Reynolds Number,” Phys. Fluids, 24, 014102. [CrossRef]
Moriai, H., Wang, P., Hayashi, J., Nakatsuka, N., Akamatsu, F., Kurose, R., and Komori, S., 2009, “Recirculation Behavior in Subscale Gas Turbine Combustor,” Proceedings of the 7th World Conference on Experimental Heat Transfer, Fluid Mechanics and Thermodynamics (ExHFT-7), Krakow, Poland, June 28–July 3, MT-22 (CDROM).
Moriai, H., Hayashi, J., Wang, P., Nakatsuka, N., Akamatsu, F., Kurose, R., and Komori, S., 2011, “Optical Measurement of Spray Combustion Fields in a Sub-Scale Model for Aircraft Gas Turbine Engine Combustor,” J. Gas Turbine Soc. Jpn., 39, pp. 125–130.

Figures

Grahic Jump Location
Fig. 1

Subscale (1/2) rectangular combustor with air-blast type swirl fuel nozzle; (a) combustor, (b) swirl fuel nozzle

Grahic Jump Location
Fig. 2

Computational domain, grids, and nozzle configuration for LES; (a) computational domain, (b) grids in burner, (c) nozzle configuration

Grahic Jump Location
Fig. 3

Initial droplet size distributions of injected fuel droplets; (a) φ = 0.91, (b) φ = 1.09, (c) φ = 1.29

Grahic Jump Location
Fig. 4

Isosurface of instantaneous gaseous temperature (T = 1800 K) and instantaneous droplet distribution in case of ϕ = 1.09

Grahic Jump Location
Fig. 5

Distributions of instantaneous gaseous temperature T and isosurface of droplet volume fraction on the x–y and y–z planes in case of ϕ = 1.09. The color range of black to white corresponds to 300–2000 (K) and white dots show the fuel droplets. (a) xy plane at z = 0; (b) yz plane at z = 5 mm.

Grahic Jump Location
Fig. 6

Comparisons of radial profiles of predicted and measured time-averaged streamwise droplet velocities u¯d at three streamwise locations in case of ϕ = 1.09, together with predicted time-averaged streamwise flow velocities u¯; (a) x = 4.0 mm, (b) x = 8.0 mm, (c) x = 12.0 mm

Grahic Jump Location
Fig. 7

Comparisons of radial profiles of predicted and measured Sauter mean diameters D32 at three streamwise locations in case of ϕ = 1.09

Grahic Jump Location
Fig. 8

Radial profiles of predicted time-averaged azimuthal droplet velocities and flow velocities, u¯θ,d and u¯θ, at three streamwise locations in case of ϕ = 1.09; (a) x = 4.0 mm, (b) x = 8.0 mm, (c) x = 12.0 mm

Grahic Jump Location
Fig. 9

Axial profiles of predicted time-averaged gaseous temperature T¯ and NO mass fraction Y¯NO in case of ϕ = 1.09, together with the measured values at the center of the combustor exit; (a) gaseous temperature, (b) NO mass fraction

Grahic Jump Location
Fig. 10

Distributions of predicted instantaneous streamwise flow velocity u, gaseous temperature T, and mass fractions of evaporated fuel, oxygen O2, carbon dioxide CO2, and carbon monoxide CO, Yfuel, YO2, YCO2, and YCO, on the x–y plane in case of ϕ = 0.91. The color ranges of black to white correspond to −100–100 (m/s) for u, 300–2000 (K) for T, 0 – 1.0 × 10−4 (-) for fuel, 0–0.3 (-) for O2, 0–0.2 (-) for CO2, and 0–0.2 (-) for CO, respectively. (a) u, (b) T, (c) fuel, (d) O2, (e) CO2, (f) CO.

Grahic Jump Location
Fig. 11

Distributions of predicted instantaneous streamwise flow velocity u, gaseous temperature T, and mass fractions of evaporated fuel, oxygen O2, carbon dioxide CO2, and carbon monoxide CO, Yfuel, YO2, YCO2, and YCO, on the x–y plane in case of ϕ = 1.09. The color ranges are as shown in Fig. 10. (a) u, (b) T, (c) fuel, (d) O2, (e) CO2, (f) CO.

Grahic Jump Location
Fig. 12

Distributions of predicted instantaneous streamwise flow velocity u, gaseous temperature T, and mass fractions of evaporated fuel, oxygen O2, carbon dioxide CO2, and carbon monoxide CO, Yfuel, YO2, YCO2, and YCO, on the x–y plane in case of ϕ = 1.29. The color ranges are as shown in Fig. 10. (a) u, (b) T, (c) fuel, (d) O2, (e) CO2, (f) CO.

Grahic Jump Location
Fig. 13

Distributions of predicted instantaneous mass fraction of hydroxyl radical OH, YOH, on the x–y plane in cases of ϕ = 0.91, 1.09, and 1.29. The color range of black to white corresponds to 0–5.0 × 10−3 (-). (a) φ = 0.91, (b) φ = 1.09, (c) φ = 1.29.

Grahic Jump Location
Fig. 14

Distributions of predicted instantaneous mixture fraction Z and product mass fraction C on the x–y plane in case of ϕ = 0.91. The color ranges of black to white correspond to 0–0.2 (-) for Z and 0–0.4 (-) for C, respectively. (a) Z, (b) C.

Grahic Jump Location
Fig. 15

Distributions of predicted instantaneous mixture fraction Z and product mass fraction C on the x–y plane in case of ϕ = 1.09. The color ranges are as shown in Fig. 14. (a) Z, (b) C.

Grahic Jump Location
Fig. 16

Distributions of predicted instantaneous mixture fraction Z and product mass fraction C on the x–y plane in cases of ϕ = 1.29. The color ranges are as shown in Fig. 14. (a) Z, (b) C.

Grahic Jump Location
Fig. 17

Distributions of predicted instantaneous mass fraction of nitrogen oxide NO, YNO, on the x–y plane in cases of ϕ = 0.91, 1.09, and 1.29. The color range of black to white corresponds to 0–6.0 × 10−5 (-). (a) φ = 0.91, (b) φ = 1.09, (c) φ = 1.29.

Grahic Jump Location
Fig. 18

Distributions of predicted instantaneous soot radiation energy emission, Esoot, and mass fraction of acetylene C2H2, YC2H2, on the x–y plane in case of ϕ = 0.91, together with direct photographs taken in the experiments. The color range of black to white in C2H2 corresponds to 0–2.0 × 10−2 (-). (a) Esoot with gaseous temperature, (b) C2H2, (c) experiment.

Grahic Jump Location
Fig. 19

Distributions of predicted instantaneous soot radiation energy emission, Esoot, and mass fraction of acetylene C2H2, YC2H2, on the x–y plane in cases of ϕ = 1.09, together with direct photographs taken in the experiments. The color range in C2H2 is as in Fig. 18. (a) Esoot with gaseous temperature, (b) C2H2, (c) experiment.

Grahic Jump Location
Fig. 20

Distributions of predicted instantaneous soot radiation energy emission, Esoot, and mass fraction of acetylene C2H2, YC2H2, on the x–y plane in cases of ϕ = 1.29, together with direct photographs taken in the experiments. The color range in C2H2 is as in Fig. 18. (a) Esoot with gaseous temperature, (b) C2H2, (c) experiment.

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In