Aerospace Technic and Technology

The problem of numerical simulation of the gas flow with the combustion of atomized liquid fuel was solved (the equilibrium combustion model pdf was used along with the partially mixed mixture model) in the annular combustion chamber of a gas turbine engine. Numerical modeling was performed in Ansys Fluent calculation complex. The purpose of the calculations was to simulate the radial and circumferential unevenness of the gas temperature pattern at the outlet of the combustion chamber. As a result of the calculations, it was found that the accuracy of modeling the radial and circumferential unevenness of the gas temperature pattern at the outlet of the combustion chamber is unsatisfactory when using the k–e turbulence model with the initial settings for the Ansys Fluent calculation complex. Moreover, the maximum value of the radial non-uniformity of the gas temperature pattern at the outlet of the combustion chamber exceeded the value obtained in the experiment by 12.61 %, and the maximum value of the circumferential non-uniformity by 12.69 %. To improve the accuracy of modeling the temperature pattern non-uniformity at the outlet of the combustion chamber, a numerical experiment was conducted to study the effect of the degree of turbulent diffusion of gas components on the value of temperature pattern non-uniformity. To reduce the non-uniformity of the temperature pattern at the outlet of the combustion chamber, the degree of turbulent diffusion of gas components was increased with respect to the initial version of the calculation, performed using the k–e model of turbulence with the initial settings for the Ansys Fluent calculation complex, by reducing the turbulent Schmidt number Sc in the turbulence model. For the initial settings of the k–e turbulence model in the Ansys Fluent calculation complex, the turbulent Schmidt number Sc = 0.85. A numerical experiment was performed for the values of Sc = 0.6, Sc = 0.4, and Sc = 0.2. The results of a numerical experiment confirmed the influence of the turbulent Schmidt number Sc on the result of calculating the gas temperature pattern at the outlet of the combustion chamber; as the value of Sc decreases, the level of the circumferential and radial non-uniformities of the gas temperature pattern decreases. However, the degree of reduction of radial and circumferential irregularities with a decrease in Sc is different. Therefore, to ensure high accuracy in calculating both the circumferential and radial non-uniformities of the gas temperature pattern, it was proposed to use a variable value of the turbulent Schmidt number Sc depending on the gas temperature instead of a constant value. The functional dependence of the turbulent Schmidt number Sc on the gas temperature was implemented in the Ansys Fluent calculation complex using the user function (UDF). The results of modeling the gas temperature pattern using the proposed UDF function for the turbulent Schmidt number Sc are in satisfactory agreement with the experimental data for both radial and circumferential non-uniformities of the gas temperature pattern at the outlet of the combustion chamber.

combustion chamber; turbulent Schmidt number; gas temperature pattern; computer simulation; UDF; ANSYS Fluent

Lefevr A. Protsessy v kamerakh sgoraniya GTD [Processes in the combustion chamber of a gas turbine engine]. Moscow, World Publ., 1986. 566 p.

Kostyuk, V. E., Kirilash, E. I., Konoplenko, A. A. Analiz tochnosti chislennykh otsenok perepada davleniya na ZhT, raspredeleniya raskhoda vozdukha po otverstiyam ZhT, koeffitsientov raskhoda otverstii ZhT i neravnomernosti polya temperatur na vykhode KS, dostignutoi drugimi avtorami, a takzhe ispol'zovannykh imi metodicheskikh priemov. Vypolnenie testovykh chislennykh raschetov koeffitsienta raskhoda otverstiya v tonkoi stenke na setkakh s razlichnym razresheniem i analiz ikh tochnosti. Formirovanie vyvodov i gipotez: Rabochie materialy etapa № 1. TS № 65 ot 28.12.2009 [An analysis of the accuracy of numerical estimates of the pressure difference across the flame tube, the distribution of air flow over the openings of the flame tube, the flow coefficients of the openings of the flame tube, and the unevenness of the temperature field at the exit of the combustion chamber, achieved by other authors, as well as the methods used by them. Test numerical calculations of the flow coefficient of the hole in the thin wall on grids with different resolutions and analysis of their accuracy. Formation of conclusions and hypotheses: Working materials of stage No. 1. Technical information No. 65 of 12.28.2009]. Kharkov, 2009. 159 p.

Kostyuk, V. E., Kirilash, E. I., Konoplenko, A. A. Formirovanie chislennykh mode-lei, prigodnykh dlya otsenki effektivnosti kon-struktivnykh meropriyatii, napravlennykh na opti-mizatsiyu kharakteristik KS po neravnomernosti polya temperatur na vykhode KS: Rabochie materia-ly etapa № 5. Dogovor № 203-22/06 ot 28.04.2006 [Formation of numerical models suitable for assessing the effectiveness of structural measures aimed at optimizing the characteristics of the compressor by the uneven temperature field at the output of the compressor: Working materials of stage No. 5. Contract No. 203-22 / 06 of 04/28/2006]. Kharkov, 2009. 39 p.

Averin, S. I., Minaev, A. N., Shvydkii, V. S., Yaroshenko, Yu. G. Mekhanika zhidkosti i gaza [Fluid and gas mechanics]. Moscow, Metallurgija – Metallurgy, 1987. 304 p.

Users Guide FLUENT – Fluent Co., 2005. 2216 p.

He, G., Guo, Y., Hsu, A. T., Brankovic, A., Syed, S., Liu, N. S. The Effect of Schmidt Number on Turbulent Scalar Mixing in a Jet-In-Crossflow. ASME Paper 99-GT-0137. 1999. 12 p.

King, Phil T. CFD Predictions of Isothermal Fuel-Air Mixing in a Radial Swirl Low NO x Combustor using Various RANS Turbulence Models. ASME Paper GT2012-69299, 2012. 11 p.