автореферат диссертации по энергетическому, металлургическому и химическому машиностроению, 05.04.12, диссертация на тему:Камеры сгорания газотурбинных двигателей

На правах рукописи

Митрофанов Валерий Александрович

КАМЕРЫ СГОРАНИЯ ГАЗОТУРБИНННЫХ ДВИГАТЕЛЕЙ: МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, МЕТОДОЛОГИЯ РАСЧЕТА, КОНЦЕПЦИЯ ОПТИМАЛЬНОГО ПРОЕКТИРОВАНИЯ.

Специальность 05.04.12. Турбомашины и комбинированные турбоустановки.

автореферат диссертации на соискание ученой степени доктора технических наук.

Санкт-Петербург, 2004г.

Работа выполнена на Федеральном Государственном унитарном предприятии «Завод имени В Я. Климова»

ОФИЦИАЛЬНЫЕ ОППОНЕНТЫ

доктор технических наук, профессор Лапшин Кирилл Леонидович, доктор технических наук, профессор Сударев Анатолий Владимирович, доктор технических наук, профессор Марчуков Евгений Ювенальевич.

ВЕДУЩАЯ ОРГАНИЗАЦИЯ

Филиал ОАО «Силовые машины» «ЛМЗ» в Санкт-Петербурге

Защита состоится «27» апреля 2004 г. в \6_ час.00 мин. на заседании диссертационного совета Д212.229.06 Государственного образовательного учреждения высшего профессионального образования «Санкт-Петербургский государственный политехнический университет» по адресу: Санкт-Петербург, Политехническая ул., д.29, главное здание.

С диссертацией можно ознакомиться в библиотеке ГОУ «СПбГПУ»

Автореферат разослан

Ученый секретарь диссертационного совета, Д212.229.06 д.т.н., профессор

Н.Н. Кортиков

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ.

Актуальность. Современное состояние развития теории камер сгорания свидетельствуют о приоритетном направлении, которое характеризуется совершенствованием расчетных «поверочных» методологий на основе численного решения системы уравнений гидроаэромеханики.

Поэтому, распространенный в настоящее время метод создания камер сгорания можно классифицировать терминологией - «метод проб и ошибок» - проектирование по прототипу. Подобный подход неизбежно связан с большим объемом проектных и расчетно-экспериментальных работ при доводке и совершенствовании характеристик рабочего процесса создаваемой или модернизируемой камеры сгорания.

Таким образом, разработка направления в теории камер сгорания — совершенствование методологии «целенаправленного» проектирования - определение размеров проточной части по заданным характеристикам рабочего процесса, решение задач оптимизации конструкции, выработка концепции проектирования совершенной конструкции камеры сгорания, представляет собой актуальную задачу, имеющую важное научное и практическое значение.

Цели работы. 1.Разработка направления в теории камер сгорания - совершенствование методологии «целенаправленного» проектирования. Суть названной методологии состоит в определении размеров проточной части устройства по заданным величинам характеристик потока сплошной среды на выходе из него и решении задач оптимизации конструкции.

2.Создание интегральной математической модели камеры сгорания газотурбинного двигателя для решения задач оптимального проектирования.

3. Выработка, на базе созданной математической модели, концепции проектирования совершенной конструкции камеры сгорания - т.е. максимально отвечающей всем требованиям, определяемым назначением двигателя.

Научная новизна.

1. На основе преобразований уравнений гидроаэромеханики и тепломассообмена в потоке химически реагирующей смеси теоретически обоснованы гидродинамические и геометрические критерии и параметры камеры сгорания газотурбинного двигателя, определяющие протекание рабочего процесса в ней.

2. Разработана система уравнений для расчета характеристик стационарного процесса в камере сгорания: коэффициента полноты сгорания топлива, индексов эмиссии оксидов азота, монооксидов углерода, несгоревших углеводородов, бензо-(а)-пирена, числа дымности,

коэффициентов окружной и радиальной потерь пол-

библиотека I С.ПетербЯ* /V /V' оэ 1007»**''

ного давления и температурного состояния стенок жаровой трубы.

4. Разработаны математическая модель нестационарных процессов - воспламенения топливовоздушной смеси в различных высотно-климатических условиях, прекращения горения при обеднении топливовоздушной смеси и критерий для оценки стабильности течения газа в диффузоре камеры сгорания, показана связь величины этого критерия со стабильностью окружной неравномерности поля температуры газа.

5. Теоретически обосновано и экспериментально подтверждено существование минимально допустимых значений размеров камеры сгорания — высоты (для кольцевой камеры сгорания) или диаметра жаровой трубы, при которых влияние «холодных» стенок не влияет на процесс горения; длины части жаровой трубы до отверстий подвода воздуха «на смешение» с продуктами сгорания.

6. Обоснованы на основании теоретического анализа и подтверждены многочисленными испытаниями принципы физического моделирования процесса горения в камере сгорания любой схемы и конструкции (прямоточной, противоточной, кольцевой, петлевой, ра-диально-осевой) на ее простейшей цилиндрической модели.

Практическая ценность и внедрение.

- создан алгоритм и программа расчета размеров камеры сгорания для газотурбинных двигателей различного назначения;

- выработаны рекомендации по совершенствованию характеристик модернизируемых камер сгорания;

- определены допустимые пределы для величин геометрических параметров и размеров элементов конструкции камеры сгорания;

- материалы диссертации использованы на ФГУП «Завод имени В.Я.Климова», в частности при создании камер сгорания двигателей: ТВа - 3000, малоразмерного турбовального вертолетного двигателя, высокотемпературного турбореактивного двигателя пятого поколения, приводов стационарных энергоустановок ГТП - 1,5, ГТЭ - 100 и модернизации камер сгорания двигателей: ТВЗ - 117, РД - 33, вспомогательной силовой установки ГТДЭ - 117; материалы диссертации также используются на ОАО «НПО «Сатурн» и в учебном процессе СПбГПУ при подготовке специалистов соответствующего профиля

Публикации. Основные результаты диссертации изложены в 20 печатных работах и трех авторских свидетельствах.

Апробация работы. Результаты диссертационной работы доложены и обсуждены: на ХЫП сессии РАН «Повышение эффективности газотурбинных и парогазовых установок» (Санкт-Петербург, 1996г.); на ХЬУ сессии РАН «Применение газотурбинных и парогазовых

установок в энергетике и на газопроводах» (Санкт-Петербург, 1998г.); на Семинаре по экологически чистым камерам сгорания (Московское представительство фирмы СНЕКМА, Москва 1998г); на Международной конференции «Проблемы промышленной теплотехники» (Киев, 1999г.); на специализированных научно-технических семинарах «Процессы горения, теплообмена и экология тепловых двигателей» (Самара, 1999г.,2000г.); на семинарах кафедры турбинных двигателей и установок, кафедры теоретических основ теплотехники Санкт-Петербургского Государственного политехнического университета; на семинарах научно-технического совета завода им. В.Л. Климова; на Международной научной конференции «Двигатели XXI века», посвященной 70-тилетию ЦИАМ, Москва,2000г.; на Общероссийской конференции «Первые Уткинские чтения», «Военмех», СПб, 2002г.; на II Международной конференции «Проблемы промышленной теплотехники», Киев,2001; на IV Всероссийской научно-технической конференции «Проблемы горения, теплообмена и экология тепловых двигателей», Самара,2002г.; на Международной научно-технической конференции посвященной памяти Генерального конструктора аэрокосмической техники академика Н.Д. Кузнецова, Самара, 2001 г; на L научно-технической сессии РАН по проблемам газовых турбин, СПб, 2003 г; на международном симпозиуме по горению и загрязнению атмосферы, СПб, 2003г.

Личный вклад автора.

1. Проведение теоретических исследований, заключающихся в математических преобразованиях, анализе фундаментальных уравнений механики сплошной среды и тепломассообмена в этой среде.

2. Организация и проведение экспериментальных исследований, обобщение материалов по ним.

3. Разработка системного подхода, математических моделей и концепций, представленных в работе.

Структура диссертации. Диссертация состоит из введения, 11-ти разделов, содержит 298 страниц основного текста, 73 рисунка, списка литературы из 137 наименований и двух приложений.

СОДЕРЖАНИЕ РАБОТЫ.

Во введении сформулирована концепция и основные положения диссертационной работы.

В первом разделе анализируются тенденции развития науки о камере сгорания.

С развитием и совершенствованием поверочных методологий наметилось отставание в развитии методов «непосредственного» или «целенаправленного» определения размеров конструкции для обеспечения требуемых от нее характеристик. Совместное развитие двух направлений - совершенствование поверочных расчетов, в том числе на базе численного интегрирования и совершенствование методологии «целенаправленного» расчета размеров камеры сгорания, качественным образом сократит трудоемкость процесса создания камеры сгорания, сведя ее лишь к доводке, не выходящей за рамки изменений размеров конструкции «на бумаге».

Представленная работа имеет целью развитие второго метода - «целенаправленного» расчета размеров камеры сгорания по заданным характеристикам рабочего процесса в ней на основе интегральной математической модели камеры сгорания газотурбинного двигателя.

Во втором разделе предложена концепция математического моделирования камеры

для решения задач оптимизации конструкции. Математическое моделирование камеры сгорания производится с помощью гидродинамических, геометрических, термодинамических критериев и параметров, установленных па основании преобразований, выполненных автором (представленных в приложении 1), системы дифференциальных уравнений переноса

энергии, концентрации реагирующих веществ, количества движения и уравнения неразрывности (с учетом уравнения состояния газа) в движущемся по каналу переменной геометрии потоке сплошной среды, при наличии процессов химического реагирования в нем:

SV| :SVi =HUWP + diss;SV2 ■= pWP\SVj =-gradP + [V + ^ jgrad divW;

SV4 =ftGrn,Gn,P,r);SV5 =i2(z0,G0,P,T).

где Т -температура газа; С— массовая концентрация конечного продукта реакции (С02,И20, N02); Р,р — давление и плотность газа; — вектор скорости потока газа (ос-

вязкости газа; — динамический коэффициент турбулентной «вязкости» газа; Рг,$с -

число Прандтля и Шмидта; Ргт,8ст -турбулентные числа Прандтля и Шмидта; ¡Уг —массовая скорость химической реакции образования конечного продукта;

где к, — константа соударения молекул, к0 ~Рт\т- показатель степени, характеризующий влияние давления на скорость химической реакции; Е— энергия активации; С^ — массовая

концентрация окислителя в потоке газа; — массовая концентрация исходного реагента {С„Нт, СО, N\, углерод) в потоке газа; Я — универсальная газовая постоянная; п\, Пг — порядок реакции окисления углеводородных топлив, «ь И2=соп81; Ни — количество тепла, выделяющееся при сгорании тошшвовоздушной смеси:

где tj — коэффициент полноты сгорания топлива; Hv -низшая удельная теплота сгорания топлива; Q —вторая вязкость газа; diss- диссипативная функция, значение которой характеризует количество механической энергии жидкости (газа), которое преобразуется вследствие трения во внутреннюю энергию за единицу времени в единице объема сплошной среды.

Функции S^.S^характеризуют процессы фазового перехода (испарения, плавления), зависящие от термодинамических параметров и суммарной площади поверхности частиц, переходящих в другое агрегатное состояние. Величина суммарной площади поверхности частиц зависит от их количества, которое является функцией значений их характерного размера и массового расхода. При протекании процессов химического реагирования исходных компонентов в одной фазе - массовый расход частиц исходного

реагента и окислителя не в газообразпой фазе (капли, твердые частицы); - характер-

реднешюй по Рейнольд су для турбулентного потока); f.lt - коэффициент динамической

ный размер этих частиц; в дальнейшем величины характерных размеров будем обозначать одним символом - х- При записи уравнения (1) положено —— » 0 (что справедливо для дозвуковых газовых потоков в проточной части камер сгорания).

Для нестационарных процессов воспламенения топлива, прекращения его горения и

неустойчивого течения газа по каналу переменной геометрии (например, в диффузоре) - ^ 0, поэтому применена форма записи дифференциальных уравнений со знаками не-

др\у ,

3/

равенств. В результате преобразований с использованием основных положений векторного анализа (понятием о «линиях тока» и выражением для полного дифференциала функции), интегральной формулы Грина, автором получено выражение, характеризующее интенсивность переноса субстанции в потоке сплошной среды при наличии источников и стоков энергии и вещества:

где Ур - величина характерного объема; - величина площади поверхности, охватывающей

характерный объем; —— ~1'р\Р- площадь поперечного сечения капала; 1р —относительная

длина характерного объема; индексы обозначают изменение величины параметра

и ее среднюю величину от выхода из рассматриваемого объема канала ко входу в этот объем;

— величины, характеризующие неравномерность рас-

— V —

пределения параметров потока сплошной среды в рассматриваемом объеме канала (парамет-

ры в числителе представляют собой местные значения); О - массовый расход среды через объем канала; - Яе.

К характерным объемам отнесены:

1) объемы, в которых происходит окисление исходных реагентов до конечных продуктов реакции (объемы реакционных зон);

2) объемы, в которых происходит формирование газодинамических параметров потока воздуха, перед его поступлением в объемы зон химического реагирования; для этих объемов = 0.

Левая часть написанного выражения характеризует скорость переноса потоком сплошной среды энергии — — вещества — количества

движения - Правая часть - скорость преобразования материи из одной формы в

другую с учетом молекулярно-кинетических и турбулентных свойств потока сплошной среды.

Автором выполнен математический анализ модельных дифференциальных уравнений переноса потоком сплошной среды параметров турбулентности (например, кинетической энергии турбулентных пульсаций, диссипации кинетической энергии турбулентных пульсаций, частоты турбулентных энергосодержащих пульсаций - в зависимости от модели турбулентности) и выявлены параметры, определяющие изменение величины динамического коэффициента турбулентной вязкости -

С целью определения величины динамического коэффициента турбулентной «вязкости», которая, в соответствии с представлениями о механизме турбулентного обмена количеством движения, может быть представлена соотношением:

удельная кинетическая энергия турбулентных пульсаций; L - интегральный масштаб турбулентности; к =—(м'2 +у'2 и',у',М>' - проекции пульсационной состав-2

ляющей вектора скорости потока газа на декартовы оси координат (знак, означающий осреднение по времени, опускаем), автором предложена и подтверждена опубликованными экспериментальными данными критериальная модель турбулентности для канала переменной геометрии (на основании интегрирования и преобразований уравнений Рейнольдса для

установившегося турбулентного течения и уравнения Навье-Стокса для нестационарного процесса).

В результате получено выражение для комплекса, входящего в (2):

где с - коэффициент пропорциональности.

6Р V,

Отношения

являются интегральными параметрами, характеризующими

р Г У л тп л тЛ

геометрическую форму канала и определяющими процессы турбулентного обмена количеством движения в нем.

^ 1

Физический смысл критерия ——состоит в том, что он представляет отношение

скоростей микро переноса (молекулярного переноса) к скорости макро переноса (конвективного переноса) количества движения в потоке сплошной среды. С другой стороны:

Комплекс —у имеет размерность [кг/(м3 с)] и характеризует скорость молекулярного переноса субстанции в сплошной среде.

Соответственно, комплексы характеризуют соотношение скоро-

стей молекулярного и конвективного переноса энергии и вещества.

Теоретически и экспериментально установлены факторы, определяющие неравномерность распределения параметров сплошной среды в пространстве, занимаемом этой средой при ее движении в нем. К таким факторам относятся вышеупомянутые гидродинамиче-

ские критерии ~——^"р"'"^—' комплексы ^ , и степень загромождения проточной части канала различными препятствиями. При наличии в канале препятствий для потока, неравномерность распределения параметров определяется процессами'Турбу-лентного смешения разделенных препятствием потоков. Коэффициенты неравномерности распределения скорости, концентрации, температуры, параметров турбулентности по характерным объемам канала переменной геометрии —А у ^Ау А ^ представлены как

функции геометрических критериев и параметров, характеризующих геометрию проточной для движения сплошной среды части канала переменной геометрии.

Третий раздел посвящен изучению и описанию процессов в камере сгорания газотурбинной установки на основе разработанного автором интегрального метода.

Определены характерные размеры, безразмерные и размерные геометрические комплексы, определяющие протекание процессов в камере сгорания газотурбинного двигателя.

Проиллюстрировано представление о комплексе физико-химических явлений, происходящих в камере сгорания, на основе разработанного интегрального подхода к описанию явлений микро и макропроцессов в движущемся, при наличии в нем турбулентных пульсаций, потоке сплошной среды, испытывающей преобразование форм материи из одной в другую (в результате происходящих химических реакций). Определены объемы характерных зон камеры сгорания: объем диффузора; объем части жаровой трубы до отверстий смесителя; весь объем жаровой трубы; объем части жаровой трубы от отверстий смесителя до сечения выхода из нее.

На основании уравнения баланса тепла, выделяемого в процессе химической реакции в объеме жаровой трубы и тепла, передаваемого через боковую поверхность стенок жаровой трубы (к наружному корпусу и от него в окружающее пространство), теоретически обосновано существование минимально допустимого поперечного размера жаровой трубы камеры сгорания (для металлического материала стенок жаровой трубы). С использованием экспериментальных данных определено его численное значение (для кольцевой и одиночной камеры сгорания). Физически, существование минимально допустимого поперечного размера жаровой трубы, объясняется усилением отвода тепла из зоны химической реакции (горения) по мере приближения к ней «холодных» стенок жаровой трубы.

С применением теории турбулентного смешения газовых струй и статистического анализа серийно реализованных конструкций жаровых труб, определена минимально допустимая величина длины части жаровой трубы, на которой должны быть расположены отвер-

стая подвода воздуха для смешения с ним «горячих» продуктов сгорания - т.н. «длина зоны горения».

В четвертом разделе, на базе изложенного материала, представлена в критериально-параметрической форме математическая модель камеры сгорания в виде функции критериев и параметров, определяющих рабочий процесс в ней:

где Тк и Рк - температура и давление на входе в реакционную зону, равны температуре и давлению воздуха на входе в камеру сгорания; Тг —температура газа на выходе из объема реакционной зоны; (X — среднее значение коэффициента избытка воздуха в объеме реакци-

Р Р

онной зоны; —степень расширения диффузора камеры сгорания; -отношение

площади кольцевого канала камеры сгорания к эффективной площади всех отверстий в жа-

ровой трубе (//— коэффициент расхода);

— отношение площади жаровой трубы к

эффективной площади всех отверстий в ней;

-отношение площадей жаровой трубы и

выходного сечения ее; -относительная пропускная способность завихрителя фронтового устройства, представляющая собой отношение эффективной площади завихрителя фронтового устройства к суммарной эффективной площади всех отверстий жаровой тру— иР иР

бы- - отношение эффективной площади отверстий подвода воз-

Е/^г Х/^

духа для охлаждения стенок жаровой трубы к эффективной площади всех отверстий в жаровой трубе (для конструкций с конвективно-пленочным охлаждением стенок);

-интенсивность крутки воздушного потока на выходе из завихрителя, как

функция его геометрической характеристики и представляющая собой отношение тангенциальной и осевой составляющих скорости воздуха в центре кольцевой струи, где

—радиус, проходящий через центры тяжести межлопаточных каналов; —угол наклона лопаток к оси на выходе из завихрителя; — площадь на выходе из межлопаточных каналов; — наружный радиус завихрителя; —количество форсунок,

приходящееся на площадь поперечного сечения жаровой трубы; —относительная координата разделения потоков (топливных, воздушных, газовых).

В конце раздела дан анализ физической сущности критериев, входящих в функцию

(3).

Входящие в функцию (3) параметры характеризуют: скорость протекания химических реакций окисления топлива - вид топлива; скорость конвективного переноса

энергии, вещества и количества движения соотношение скоростей молекулярного и

конвективного переноса количества движения, энергии, вещества -

Ь X Рг Х

интенсивность процессов фазового превращения, химического реагирования, молекулярного тепло-массообмена в движущемся потоке сплошной среды - Р^,Т^,Тг,СС,вид топлива,^; интенсивность процессов турбулентного тепломассообмена -параметры в функции (3), ха-рактеризущие геометрию проточной части камеры сгорания.

Пятый раздел посвящен описанию установок и стендов для испытаний камер сгорания, величин погрешностей определения параметров в процессе проведения экспериментов, методов обработки данных исследований.

В шестом разделе представлены обработанные экспериментальные данные испытаний

в виде функциональной зависимости параметров, характеризующих рабочий процесс в камере сгорания от ее критериев (входящих в функцию (3)):

А =-^втВ\ Дл

1 -п = (Е1„с + 0,232Е1со) Ш}

ж.

288;

xexp

(и V

N.

Ф

V^ )

+ с.

1+ d.-=-%— + d. 'Mar 4

M

ож)

[l + d5W +diW1]x

l + e.

Г . \2

A;

(l + e3 expí^jx

i.+L MFf"

1

\ Zjf" o¡.

ElHC=K(PrJh%

V

\ ж J

xexp

'N

4M

l + e.

/Ч

V>J -

(l + e6exp/.r)x

Г ^ 1

•• , " 0J07

I. +'

V 42>Ъ

£/ =APl

Е1бп=А0Е1НС ; f G \

V

Vir J

1 + 6,

f Рж УьÍ fx )

хехр (b3r + ¿4í«F3)

l+¿ ехр

«

ж J

SN = MPI

( G

V -V

\ 'ж r зг J

l + cu.

-л-CO.

2Ä 41 ^ож)

\2'

X tOjïT + <04(fP)2] ехр^г,

рж J

1 + е.

\2'

8 Д

з J

(l + e explj;

>0,8;

S0,05м;

exp

U

ч2

Х^ОЖ,

Х^ож,

+1

\с13

У^вых J

1 - cu exp

I M

ОЖ J щ

-ís-ái;

где Е1С0. Е1нс, Е1вп,Е1

-индексы эмиссии монооксидов углерода, несгоревших углево-

дородов, бензо- (аг)-пирена, оксидов азота [г/кг топлива], - число дымности [%] (уменьшите отражательной способности бумажного фильтра после пропускания через него пробы отработавших газов); в ,0. —коэффициенты максимальной окружной и радиальной не-

Л12Х ГШИХ

равномерности поля температуры газа на выходе из камеры сгорания; ¿о - теоретическое количество воздуха, необходимое для сжигания 1кг топлива; а у —коэффициент избытка воздуха в части жаровой трубы до отверстий смесителя; —шаг по окружности между форсунками; —длина части жаровой трубы до отверстий смесителя; 1Д —длина диффузора, — высота жаровой трубы; —соответственно диаметр на выходе из разделителя

воздушных потоков («аэрационного» насадка) в завихрителе фронтового устройства и диаметр на выходе из завихрителя (при отсутствии «аэрационного» насадка -

коэффициент потерь полного давления в камере сгорания; Т^г —максимальная величина температуры стенки жаровой трубы; Узг —объем части жаровой трубы до отверстий смеси-

К

теля, — объем жаровой трубы; —объем диффузора; —относительная длина

ср

реакционного объема газосборника жаровой трубы, где =———^яых . ^^ — высота канала в выходном сечении жаровой трубы; длина части жаровой трубы от отверстий смесителя до сечения выхода из нее; — относительная длина камеры закручивания топлив-

ной форсунки (определяет диаметр капель распыленного топлива);

2.0ТБ

доля отбираемо-

го из камеры сгорания воздуха для охлаждения турбины; — скорость потока воздуха во входном сечении в камеру сгорания; - температура газа на выходе из камеры сгора-

- эмпирические коэффициенты (постоянные для определенного вида топлива). Тестирование расчетных зависимостей для описания характеристик стационарного процесса в камере сгорания, подтверждающее их достоверность, иллюстрируется графиками, где произведено сопоставление расчета с экспериментом. Для примера, в автореферате на рис. 1 иллюстрируется точность расчета индексов эмиссии оксидов азота.

Седьмой раздел посвящен нестационарным процессам в камере сгорания: процессам воспламенения тошшвовоздушной смеси, прекращения горения и нестационарному процессу течения газа в диффузоре.

При описании процессов воспламенения и прекращения горения в качестве исходной рассматривается система дифференциальных уравнений сохранения энергии и концентрации горючего в движущемся потоке газа; при описании нестационарного течения газа в диффузоре рассматривается дифференциальное уравнение сохранения количества движения в потоке газа:

для процесса воспламенения:

Я/""

для процесса прекращения горения: _ < о ■

Ы ' Зг

где Сг - массовая концентрация горючего.

При отрыве пограничного слоя от стенок диффузора поток газа вначале ускоряется, поэтому:

В соответствии с выражением (2), условия воспламенения и прекращения горения то-пливовоздушной смеси определяются дисбалансом:

- между скоростью тепловыделения и скоростью отвода тепла потоком газа из реакционного объема; ' - между скоростью расходования горючего в результате химической реакции и скоростью его поступления в реакционный объем с потоком газа (условие непрерывности протекания химической реакции). На основании соотношения (2) и обработки экспериментальных данных по испытаниям моделей и полноразмерных конструкций камер сгорания (в т.ч. пусковых воспламенителей) различных схем и размерностей с кольцевыми коническими, ступенчатыми диффузорами с центральным телом и без центрального тела, криволинейными диффузорами, получено: для процесса воспламенения топливовоздушной смеси в камере сгорания:

для процесса прекращения горения в камере сгорания:

Уог - объем очага, поддерживающего горение топливовоздупшой смеси в камере сгорания (так как объем очага, поддерживающего горение, определяется величиной «тепловой ширины зоны горения», то Уог =^Рк,Т^,Т1,а,вид топлива),тогда:

Уог ^ (л вид топлива)

где Я - коэффициент скорости потока газа, набегающего на очаг горения); Л0ж-коэффици-ент скорости потока воздуха на выходе из стабилизирующего процесс горения устройства (отверстий жаровой трубы); - объем диффузора камеры сгорания; - соответственно площадь входа и выхода из диффузора камеры сгорания (например, для прямоточной камеры

сгорания ^ = = Д ! ^вх - высота канала на входе в диффузор; с1и —ме-

дианный диаметр капель жидкого топлива; ОС ,СС

-соответственно коэффициент из-

бытка воздуха при «бедном» и «богатом» срыве горения.

В разделе приведены данные, иллюстрирующие точность расчета по полученным зависимостям характеристик нестационарных процессов - пределов воспламенения и прекращения горения топливовоздушной смеси в камере сгорания. Для примера, в автореферате представлены графики, иллюстрирующие точность расчета пределов воспламенения топли-вовоздушной смеси в камере сгорания - рис.2, и коэффициента избытка воздуха при «бедном» срыве горения - рис.3. Произведен анализ влияния величин критериев камеры сгорания на качество процесса воспламенения топливовоздушной смеси в ней и условия прекращения горения.

Показана связь значения критерия устойчивости течения газа в диффузоре камеры сгорания с величиной, характеризующей стабильность поля температуры газа на выходе из камеры сгорания.

Восьмой раздел посвящен принципам физического моделирования рабочего процесса в камере сгорания, созданным на основе интегральной математической модели в виде функции (3).

Модель представляла собой индивидуальную камеру сгорания круглого поперечного сечения диаметром, длиной и расстоянием до отверстий смесителя, равными таковым жаровой трубе полноразмерной камеры сгорания: {х^мод =(%\с- Топливная форсунка соответствовала полноразмерной конструкции. Воздух в камеру сгорания подавался либо от заводской сети, либо от баллонов высокого давления. При этом выдерживалось равенство критериев модели и натуры:

{Рх,Ту.,Тг,а,вид топлива)лад = (Рк,Тг.,Тг,а,вид топлива)^

Индекс «мод» обозначает модель, индекс «кс» - полноразмерную камеру сгорания.

Произведенное сопоставление результатов испытаний полноразмерных камер сгорания индивидуальных и кольцевых, как прямоточных, так и противоточных, и их моделей, с одинаковыми для модели и натуры величинами критериев функции демонстрирует адекватность протекания рабочего процесса в модели и полноразмерной камере сгорания.

Применение простейшей модели камеры сгорания, созданной на основе критериальной функции (3), с учетом возможности приведения режимных параметров по представленным выше уравнениям, позволяет исключить проведение автономных дорогостоящих испытаний полноразмеряой конструкции на рабочих режимах для исследования всех ее характеристик, уменьшить объем испытаний и сократить сроки доводки конструкции.

В девятом разделе приводится разработанная система уравнений для расчета всех характеристик камеры сгорания и произведен анализ влияния значений критериев камеры сгорания на характеристики ее рабочего процесса.

Система уравнений является математической моделью характеристик рабочего процесса в камере сгорания и дает возможность совершенствовать (оптимизировать) конструкцию в процессе ее проектирования с учетом ограничений, включающих и нестационарные режимы работы.

Десятый раздел посвящен разработке концепции проектирования камеры сгорания для газотурбинных двигателей трех типов:

- для стационарной энергетической установки;

- для транспортной машины (наземный, воздушный, судовой транспорт);

- для высокоманевренной машины (военная истребительная авиация);

К каждому типу двигателя предъявляются специфические требования. Главным требованием к стационарной энергетической установке является минимальный выброс токсичных компонентов, загрязняющих атмосферу и экономичность; от двигателя для транспортной машины требуется высокая экономичность, низкие выбросы загрязняющих атмосферу веществ и высокая устойчивость работы на переменных режимах работы; от двигателя, при-

меняемого в военной истребительной авиации требуется максимальная устойчивость всех режимов его работы, надежность запуска в различных высотно-климатических условиях и минимальный уровень дымления.

Вышеназванные требования определяют значения конструктивных параметров камеры сгорания для каждого типа двигателя.

В приложении 2 приведен алгоритм расчета размеров камеры сгорания для трех типов двигателей.

В расчете заданы: термодинамические параметры, размеры выходного сечения компрессора - ¥К и выходного сечения камеры сгорания - Iхарактеристика топлива, Т^ ограничения габарита.

Расчет проводится с учетом ограничений (в зависимости от назначения двигателя) по максимальному наружному (габарит двигателя) и минимальному внутреннему (вал ротора) диаметрам камеры сгорания, по длине и по допустимой температуре стенок жаровой трубы.

В результате расчета должны быть получены размеры камеры сгорания, обеспечивающие, при заданных ограничениях, заданные характеристики.

Одиннаддатый раздел иллюстрирует практическое применение разработанной интегральной математической модели камеры сгорания для совершенствования характеристик рабочего процесса в ней.

Адекватность полученных расчетным путем данных результатам экспериментальных исследований, представленная в разделе - рис.4 - 6, дополнительно (с аналогичным сопоставлением в шестом и седьмом разделе расчетных и экспериментальных характеристик рабочего процесса в камерах сгорания) подтверждает достоверность разработанной математической модели, а выработанная концепция проектирования совершенной конструкции камеры сгорания существенно сокращает объем экспериментальных работ при ее доводке.

В приложении 1 представлены математические преобразования, выполненные автором, системы дифференциальных уравнений гидроаэродинамики и тепломассообмена с наличием источников и стоков энергии и вещества. В результате этих преобразований получено выражение (2), которое обосновывает критериально-параметрическую функцию как математическую модель процесса формирования характеристик потока сплошной среды от входа в объем до выхода из объема рассматриваемого канала при наличии в нем источников и стоков энергии и вещества. На основании этого получена функция (3), представляющая собой интегральную математическую модель камеры сгорания.

В приложении 1 также произведено тестирование предложенного автором подхода к преобразованию исходной системы дифференциальных уравнений и обоснованию критери-

ально-параметрической функции на примере обобщения экспериментальных данных Нику-радзе по величине гидравлического сопротивления шероховатых труб, экспериментальных данных по исследованию теплообмена в условиях свободной конвекции, а также на примере математического обоснования с использованием выражения (2) критериальной функции, описывающей процесс теплообмена в условиях вынужденной конвекции.

В приложении 2 представлен алгоритм расчета размеров камеры сгорания, в соответствии с предложенной концепцией оптимального проектирования.

Основные результаты работы.

1. Теоретически обоснованы интегральные гидродинамические и геометрические критерии и параметры, определяющие процесс формирования характеристик движущегося потока газа в камере сгорания с происходящими в нем химическими реакциями окисления (горения) топлива.

2. Разработана интегральная математическая модель камеры сгорания газотурбинного двигателя.

Модель представляет систему математических зависимостей величин, характеризующих протекание рабочего процесса в камере сгорания газотурбинного двигателя от значений ее режимных и конструктивных параметров. Система состоит из уравнений для расчета:

- коэффициента полноты выгорания топлива,

- коэффициентов максимальной окружной и радиальной неравномерности поля температуры газа на выходе из камеры сгорания,

- величины коэффициента потерь полного давления,

- пределов воспламенения топливовоздушной смеси в камере сгорания для различных вы-сотно-климатических условий,

- пределов устойчивого горения топлива в камере сгорания,

- степени токсичности продуктов сгорания (характеризуемой индексами эмиссии загрязняющих атмосферу веществ - оксидов азота, мопооксидов углерода, несгоревших углеводородов, бензо- -пирена, частиц углерода),

- максимальной величины температуры стенок жаровой трубы.

3. На базе разработанной математической модели камеры сгорания найдены оптимальные значения таких параметров ее конструкции как число топливных форсунок, относительная пропускная способность, интенсивность крутки потока воздуха на выходе из завихрителя фронтового устройства, относительная площадь отверстий жаровой трубы, а также парамет-

ров, учитывающих особенности конструкции фронтовых устройств и форсунок. В частности, оптимизация конкретной конструкции камеры сгорания авиационного двигателя по ряду вышеназванных параметров увеличит высотность запуска двигателя более чем на 50% или снизит эмиссию оксидов азота более чем в 2 раза. Оптимальные величины геометрических параметров и критериев положены в основу концепции проектирования камеры сгорания для различных типов двигателей.

4. Обоснована концепция оптимального проектирования камеры сгорания для газотурбинного двигателя различного назначения - энергетическая установка, двигатель транспортной машины, двигатель высокоманевренной транспортной машины (военный истребитель). Суть концепции проектирования заключается в определении размеров камеры сгорания с помощью установленных значений ее критериев и параметров, обеспечивающих такие значения характеристик камеры сгорания, которые в наибольшей степени отвечают предъявляемым к двигателю требованиям.

5. Разработан алгоритм расчета размеров камеры сгорания для каждого из трех типов двигателей.

6. Выполнен анализ системы дифференциальных уравнений переноса субстанции в потоке сплошной среды с наличием источников (стоков) энергии, вещества и силы. Произведено преобразование исходной системы дифференциальных уравнений. Преобразованная форма уравнений позволяет обосновать критерии и параметры, определяющие процесс формирования характеристик потока сплошной среды с источниками (стоками) энергии и вещества в объеме рассматриваемого канала.

7. На основании преобразований уравнений Рейнольдса для установившегося турбулентного течения и уравнения Навье-Стокса для нестационарного процесса, обоснована критериальная модель турбулентности для потока сплошной среды в канале переменной геометрии, каковым является проточная часть камеры сгорания.

8. Теоретически обоснована и определена связь температурного состояния стенок жаровой трубы с ее размерностью и формой. Получена количественная зависимость максимального значения температуры стенок жаровой трубы (с металлическим материалом стенок) от величины размера - высоты жаровой трубы в поперечном сечении.

9. На базе теоретических и экспериментальных исследований определены минимально допустимые размеры - высота поперечного сечения жаровой трубы и длина соны горения» жаровой трубы камеры сгорания.

10. Теоретически обоснованы и экспериментально подтверждены принципы физического моделирования рабочего процесса в камере сгорания. Принцип физического моделирования

основан на равенстве величин разработанных критериев и параметров, определяющих характеристики рабочего процесса в камере сгорания для физической модели и полноразмерной камеры сгорания.

11. Теоретически обоснован безразмерный геометрический критерий, характеризующий стабильность течения газа в диффузоре камеры сгорания. Получена связь величины этого критерия со стабильностью окружной неравномерности поля температуры газа.

12. Разработана математическая модель для расчета характеристик нестационарных процессов: областей воспламенения топливовоздушной смеси в различных высотно-климатических условиях и устойчивости горения при обеднении топливовоздупшой смеси в камере сгорания.

13. Материалы работы использованы при создании камер сгорания двигателя ТВа - 3000, малоразмерного турбовального вертолетного двигателя, высокотемпературного турбореактивного двигателя пятого поколения, приводов стационарной установки ГПТ - 1.5, ГТЭ -100 и при модернизации камер сгорания серийно изготавливаемых турбовальных и турбореактивных двигателей: ТВЗ-117, РД-33, вспомогательной силовой установки ГТДЭ-117.

14. Гидродинамические, геометрические критерии и параметры, полученные в работе, по мнению автора, могут быть применены для изучения и описания процессов, протекающих в жидкостных реактивных двигателях и других тепломассообменных аппаратах.

Основное содержание диссертации опубликовано в следующих работах:

1. Теория и расчет образования, вредных выбросов в камерах сгорания ГТД / ВАМигрофанов, ОАРудаков, НДСаливон и др. // Промышленная теплотехника. 1999.№6.С.34- 38.

2. Организация процессов сжигания топлив с минимальным выбросом NOx и СО /ВАМигрофанов, ОАРудаков, НДСаливон и др.// Тезисы докладов XLV сессии РАН. Санкт- Петербург 30 - 1 июля 1998г.С42.

3. Зависимость эмиссии NOx от конструктивных и режимных параметров камеры сгорания газотурбинного двигателя / ВАМитрофанов, ОАРудаков, НДСаливон и др. // Теплоэнергетика. 1999. №12. С.53 - 56.

4. Обобщенная эмиссионная характеристика ГТД как функция конструктивных и режимных параметров камеры сгорания / ВАМитрофанов, ОАРудаков, НД.Саливон и др,// Теплоэнергетика. 2000.№3. С.63 - 66.

5. Математическая модель процессов образования и распета загрязняющих веществ и опгамизация камер сгорания ГТД/ В.А.Митрофанов, О.А.Рудаков, Н-Д.Саливон и др. // Теплоэнергетика. 2000. №5. С.52 - 55.

6. Камера сгорания ГТД с низким уровнем эмиссии/ ВАМитрофанов, ОАРудаков, НДС&тавон и др.// Тезисы докладов XLIII сессии РАН. 25-26 июня Санкт- Петербург 1996года.

7. Разработка критериев запуска камер сгорания ГТД/ ВАМитрофанов, ОАРудаков, НД.Саливон и дрУ/ «Вестник СГАУ» (сборник научных трудов Самарского государственного аэрокосмического университета). 1999. №2. С. 177 -183.

8. Патент № 2094705 / Камера сгорания авиационного или энергетического газотурбинного двигателя. Сигалов Ю.В., Рудаков ОА, Ефимов Е.В., Гурский С.Э., Митрофанов ВА, Федоров А.М.

9. Патент №2145669 / Способ работы газотурбинного двигателя и огневой блок для сжигания смеси топлив. Митрофанов ВА, Рудаков ОА, Саркисов АА, Саливон Н.Д., Сигалов Ю.В.

10. Патент №2171432 / Камера сгорания с оптимальным числом форсунок. Митрофанов ВА,Саркисов АА, Рудаков ОА, Саливон Н.Д., Сигалов Ю.В.

11. Экологически «чистая» камера сгорания авиационного ГТД / ВАМитрофанов, О АРудаков, НД.Саливон и др.// Тезисы докладов Международной научной конференции «Двигатели XXI века», посвященной 70-тилетию ЦИАМ. Москва.2000г.

12. Расчет устойчивости горения в ГТД / Митрофанов ВА, Рудаков ОА, Саливон Н.Д. и дрУ/ «Вестник СГАУ» (сборник научных трудов Самарского государственного аэрокосмического университета). 2000. №З.С262 - 266.

13. Влияние размерности двигателя на характеристики камеры сгорания/ ВА.Митрофанов, ОАРудаков," Н.Д.Саливон и дрУ/ Тезисы докладов Международной научно-технической конференции, посвященной' памяти Генерального конструктора аэрокосмической техники академика НДКузнецова. Самара.2001г. С.186-187.

14. Метод решения дифференциальных уравнений тепломассообмена и газовой динамики, описывающих процессы горения в камере сгорания/ Митрофанов В.А., Рудаков ОА., Саливон Н.Д. и дрУ/Пром. теплотехника. 2001. Т.23. №6. С. 15-18.

15. Система уравнений для расчета характеристик камеры • сгорания/ Митрофанов ВА, Рудаков ОА, Саливон Н.Д. и др.// Пром. тедлотехника.2001. Т.23. №3. С.83-88.

16. Митрофанов В.А., Рудаков О.А. / Определение параметров турбулентности потока газа в камере сгорания на основе интегрирования и преобразований уравнений Рейнольдса и Навье-Стокса //Тезисы докладов IV Всероссийской научно-технической конференции «Процессы горения теплообмена и экология тепловых двигателей». Самара. 2002г. С. 137-141.

17. Исследование области воспламенения топлива в камере сгорания / Митрофанов ВА, Рудаков ОА, Саливон Н.Д. и др.//Теплоэнергетика. 2003 г. №3. С.47 - 51.

18. Митрофанов В.А., Рудаков О.А., Саркисов А.А. / Интегральная математическая модель камеры сгорания //Теплоэнергетика. 2004. №2.

19. Митрофанов В А, Рудаков О А/Интегральные гидродинамические критерии, характеризующие перенос энергии, вещества и количества движения потока сплошной среды в каналах переменной геометрии при горении // Материалы общероссийской конференции «Первые Уткинские чтения». «Военмех». СПб. 2002.Т2. С.61-62.

20. Митрофанов ВА, Рудаков ОАУ Исследование интегральных критериев движения сплошной среды в канале переменной геометрии с источниками тепла и массы на основе решения дифференциальных уравнений переноса и неразрывности// Сборник «Актуальные вопросы ракетостроения». С- Петербург. 2003. Выпуск 2. С. 132 -135.

21. Митрофанов В А, Рудаков О А / Определение параметров турбулентности потока газа в камере сгорания на основе интегрирования и преобразований уравнений Рейнольдса и Навье-Стокса// «Вестник СГАУ» (сборник научных трудов Самарского государственного аэрокосмического университета).. 2003. №2. С.92-96.

22. Митрофанов ВА., Рудаков ОАУРешение в критериальной форме дифференциальных уравнений переноса для потока сплошной среды в канале переменной геометрии с источниками и стоками энергии и массы // Тезисы докладов 50-й научно-технической секции по проблемам газовых турбин. РАН. Комиссия по газовым турбинам. ЛМЗ. С- Петербург. 2003г. С.20 - 2}.

23. V.Mitrofanov and O.Rudakov. Mathematical simulation and calculation of flame ignition and quenching in gas-turbine engine / Combustion and atmospheric pollution [Edited by G.D. Roy, S.M.Frolov, A.M.Starik]. - Moscow. TORUS PRESS Ltd. 2003. -P.273-276.

г/кг топл. (расч.)

т

Р'

.3®

О 10 20 30 40

г/кг топл. (эксп.)

Точками на графике обозначены: О - ГТД -350; О -ТВ2-117; Л - ТВЗ-117; □ -РД-33; (б) - РД-ЗЗМ; О -ТВ7-117; +-ТВ7-117 (с экспериментальной камерой сгорания);Х - Тва-3000; * - ГТД - 1250; • - Р29-300; + - РД (экспериментальная камера сгорания); и , СЗ - модели камер сгорания.

Двигатели, не участвующие в обобщении: Ь. - ГТН-25 (газовое топливо);

>-ГТУ «Сименс», модель V 94.2 (газовое топливо);0 -ИК-8б;Т - СРб-50;« -модель двухзонной камеры сгорания с «богато-бедным» горением.

Рис.1 Сопоставление расчетных и экспериментальных данных по эмиссии оксидов азота.

а

б 4 2 0

10 30 ' 50 . 70

0[гж*к М3)

1 - механическое распыливание топлива, Н=0км, Тв=Тт=288К, топливо - керосин;

2 - пневматическое распиливание топлива, Н=4,5км, Тв=Тт=223К; топливо - керосин;

3 - механическое распиливание топлива, Н=4км, Тв=Тт=288К, топливо - керосин;

4 - «аэрационное» распыливание топлива, Н=4км, Тв=Тт=223К, топливо - керосин;

5 - пневматическое распыливание топлива, Н=10км, Тв=Тт=223К, топливо - керосин; 5 - пневматическое распыливание топлива, 11=6,5км, Тв=Тт=223К, дизельное топливо;

Ц #. ^, О, X, Ж - результаты испытаний различных вариантов конструкций камер

сгорания для турбовинтовых, турбовапьных, турбореактивных двигателей; _- расчет.

(Н -высота запуска, ТвЛт-температура воздуха и топлива).

Рис. 2 Сопоставление расчетных и экспериментальных данных по воспламенению тошшвовоздушной смеси в камере сгорания.

а"™ (эксп.) 100 80 60 40 20 О

О 20 40 60 . 80 100 а*"* (расч.)

♦ - ГТД-550; ГТЭ-117 (камера сгорания с вращающейся форсункой);О- ГТД-350; Ж - ГТД-1000; ГТД-1250; О - ТВ2-117; Д -ТВЗ-117; Х-ТВ7-117; □ Щ - РД-33 и модификации.

Рпс.З. Сопоставление расчетных и экспериментальных данных по «срыву» горения в камере сгорания (при «обедпении» смеси).

к

У*

Рис. 4. Иллюстрация практического применения разработанной автором интегральной математической модели камеры сгорания для улучшения характеристик запуска газотурбинного двигателя.

--расчет с применением интегральной математической модели;

экспериментальные данные: ® -исходный вариант; О-модернизированный вариант (камера сгорания с новым фронтовым устройством).

(Тв=Тт=288К; Н=0км).

Е1т^(т/кпогт)

Рис.5 Иллюстрация практического применения разработанной автором интегральной математической модели камеры сгорания для улучшения эмиссионных характеристик газотурбинного двигателя.

-- расчет с применением интегральной математической

модели; экспериментальные данные: Д- исходный вариант; А-модернизированный вариант.

отношение частоты вращения вала двигателя к его максимальном; значению в %

Рис.6. Иллюстрация практического применения разработанной автором интегральной математической модели камеры сгорания для улучшения эмиссионных характеристик газотурбинного двигателя.

--расчет с применением интегральной математической модели,

экспериментальные данные: □ - исходный вариант; ■ -модернизированный вариант.

Лицензия ЛР №020593 от 07.0S.97.

Подписано в печать ¿ШОс^ Объем в л.л.

Тираж /00 Заказ 5(3,

Отпечатано с готового оригинал-шкета, предоставленного автором, в типографии Издательства СПбГПУ 195251, Санкт-Петербург, Политехническая ул., 29.

Отпечатано на ризографе ЯЬт-2000 ЕР Постазщик оборудования — фирма "Р-ПРИНТ" Телефон: (812) 110-65-09 Факс: (812) 315-23-04

lü-4216

Введение.

Основные обозначения.

I. Общая характеристика современного состояния проблемы и постановка задачи.

II. Принципы математического моделирования камеры сгорания для расчета характеристик рабочего процесса в ней с помощью предложенной автором критериально-параметрической функции.

III. Применение разработанных интегральных критериев и параметров к изучению и описанию характеристик процессов в камере сгорания газотурбинной установки.

III. I. Принципиальная схема камеры сгорания.

111.2. Характерные размеры и безразмерные геометрические критерии камеры сгорания.

111.3. Объемы реакционных зон камеры сгорания.

111.4. О минимальном поперечном размере камеры сгорания. ф III.5. Экспериментально-теоретическое обоснование выбора величины длины части жаровой е трубы до отверстий смесителя.'.

IV. Интегральный подход к математическому моделированию камеры сгорания.

IV. 1. Представление о рабочем процессе в камере сгорания.

IV.2. Интегральная математическая модель камеры сгорания.

IV.3. Анализ физической сущности критериев, определяющих рабочий процесс в камере сгорания.

Выводы к разделу III и IV.

V. Экспериментальное исследование камеры сгорания.

V.1. Установки и стенды для испытаний камер сгорания.

V.2. Методика проведения испытаний.

V.3. Максимальные величины погрешностей при определении параметров в процессе испытаний камер сгорания.

V.4. Метод обработки экспериментальных данных.

VI. Характеристики камеры сгорания как функции ее критериев.

VI.1. Эмиссия оксидов азота.

VI.2. Эмиссия оксидов углерода.

VI.3. Эмиссия нес горевших углеводородов и бензопирена.

VI.4. Уровень дымления камеры сгорания.

VI.5. Взаимосвязь между единицами, характеризующими уровень токсичности газов на выходе из камеры сгорания.

VI.6. Коэффициент полноты выгорания топлива.

VI.7. Зависимость величин коэффициентов неравномерности поля температуры газа на выходе из камеры сгорания от ее критериев.

VI.8. Зависимость величины коэффициента потерь полного давления в камере сгорания от значений ее критериев.

VI.9. Тепловое состояние стенок жаровой трубы.,.

Выводы к разделу VI.

VII.Oco6eHHocTH применения интегрального подхода к исследованию нестационарных процессов в камере сгорания.

VII. 1. Процессы воспламенения и прекращения горения топливовоздушной смеси в камере сгорания.

VII. 1.1. Процесс воспламенения топливовоздушной смеси в камере сгорания.

VI 1.1.2. Расчет величины медианного диаметра капель распыленного жидкого топлива на режимах запуска двигателя.

VII.2. Процесс прекращения горения при обеднении и обогащении топливовоздушной смеси в камере сгорания газотурбинного двигателя.

VII.3. Критерий стабильного течения потока газа в диффузоре камеры сгорания.

Выводы к разделу VII.

VIII. Физическое моделирование рабочего процесса в камере сгорания.

IX. Система уравнений для расчета характеристик камеры сгорания.

IX. 1.Система математических зависимостей характеристик камеры сгорания от ее критериев.

IX.2. Анализ влияния критериев на характеристики камеры сгорания.

IX.3. Оптимальные значения критериев камеры сгорания.

X. Концепция проектирования камеры сгорания газотурбинного двигателя.

X. 1. Анализ технических требований к газотурбинным двигателям.

X.2.Совершенная конструкция камеры сгорания в зависимости от назначения двигателя.

XI. Практическое применение интегральной математической модели камеры сгорания.

Основные результаты работы.

В данной работе обоснован интегральный подход к изучению и математическому описанию характеристик потока движущейся сплошной среды при наличии развитых турбулентных пульсаций с источниками тепла и вещества по каналу переменной геометрии. Показано, что этот подход имеет важное прикладное значение для решения технических задач, связанных с проектированием аппаратов, в которых происходят молекулярные, конвективные, турбулентные тепло-массообменные процессы (к таким аппаратам, например, относятся камеры сгорания газотурбинных установок, подогреватели котлов и т.п.).Хотя, общепризнанным является тот факт, что основной задачей науки, изучающей закономерности процессов, происходящих в сплошной среде, является определение параметров (характеризующих состояние среды) в любой момент времени и в каждой точке пространства, занимаемого этой средой, однако, такой подход не позволяет решать технические задачи, связанные с проектированием. Параметры среды в любой точке занимаемого ею пространства можно (с какой-либо степенью достоверности) определить при наличии информации о фаничных условия и размерах проточной части канала. Однако такая задача в данной постановке создает известную долю неопределенности для разработчика аппарата, так как конечный результат заранее неизвестен.Для проектирования, с точки зрения автора, рациональной является иная постановка задачи: задав характеристики потока сплошной среды на выходе из канала переменной геометрии, определить их формирующие размеры проточной части канала. При этом необходимо выявить функциональную взаимосвязь параметров потока от критериев, характеризующих размеры пространства, ограничивающие этот поток (размеры проточной части). Универсальность такой взаимосвязи может быть обоснована на базе физических* представлений о происходящем процессе, подтвержденных математически и экспериментально.Автором накоплен большой экспериментальный материал по испытаниям камер сгорания газотурбинных двигателей различных схем и размерностей.Современные тенденции развития теории камер сгорания свидетельствуют о приоритетном направлении поиска решения «основной задачи гидроаэромеханики» применительно к горению в движущемся потоке сплошной среды (газа) методами численного решения системы дифференциальных соотношений, характеризующих состояние среды в сколь угодно малой области пространства, занимаемой этой средой [42, 46, 60, 69, 72, 88, 94, 99]. Несмотря на развитие вычислительной техники, облегчающей процесс расчета и совершен/-.твование численных методов, из-за сложной взаимосвязи протекания микро и макропроцессов, самой сложности этих процессов в камере сгорания, требуются определенного рода допущения, упрощения, эмпирические поправки, что, в конечном итоге, заставляет отходить от концепции «точного» решения. Для совершенствования этих методов, так или иначе, необходимо проведение испытаний, подтверждающих достоверность расчетов для различных конструкций камеры сгорания. Практическое использование таких методов в настоящее время требует нахождения поправочных коэффициентов, вводимых в расчет для конкретной конструкции камеры сгорания [60,72]. Но, такие способы позволяют анализировать тенденции протекания процессов в камере сгорания.Развитие данного направления позволяет предположить, что в обозримом будущем электронные вычислительные машины заменят собой установки и стенды для экспериментальной доводки камер сгорания. Однако, данное направление можно охарактеризовать термином — «поверочная» методология.Несмотря на постоянное совершенствование и развитие теории камеры сгорания, отсутствие математических зависимостей характеристик протекающего в ней процесса от величин интегральных критериев, определяющих геометрию камеры сгорания, заставляет проектировщиков (как в России, так и за рубежом) при ее разработке традиционно выбирать конкретный прототип ранее отработанной конструкции. Из всех интегральных критериев камеры сгорания, наиболее используемых в настоящее время на этапе проектирования, можно назвать «критерий форсирования» [24,46], в который входит конструктивный параметр - величина объема камеры сгорания. Поэтому, метод создания камер сгорания для проектировщика как в случае «физических», так и «электронных» испытаний можно классифицировать терминологией - «метод проб и ошибок» - проектирование по прототипу.Следующим шагом в развитии теории камер сгорания, по мнению автора, явилась бы разработка универсальной математической модели камеры сгорания, как функции ее геометрических, гидродинамических, термодинамических и теплофизических критериев, позволяющей описать все характеристики рабочего процесса в ней. Наличие такой модели обозначит развитие методологии «целенаправленного» проектирования конструкции, суть которой состоит в том, что по заданным характеристикам рабочего процесса в камере сгорания определяются размеры ее проточной части.15 ' Совместное развитие двух вышеупомянутых направлений — «поверочного» метода и методологии «целенаправленного» проектирования, качественным образом сократит трудоемкость процесса создания камеры сгорания, сведя ее лишь к доводке, не выходящей за рамки изменений размеров конструкции «на бумаге».Представленная работа посвящена развитию метода «целенаправленного» проектирования конструкции камеры сгорания.Наличие большого экспериментального материала по испытаниям камер сгорания в сочетании с использованием положений, полученных в результате анализа и математического преобразования фундаментальных уравнений газовой динамики и тепломассообмена в сплошной движущейся среде, создает предпосылки для успешного решения поставленной автором задачи. '• II. Принципы математического моделирования камеры сгорания для расчета характеристик рабочего процесса в ней с помощью предложенной автором критериально-параметрической функции.Для решения поставленной в работе задачи, с точки зрения автора, целесообразно искать систему функциональных связей следующего вида: Zi<(>)f(Xij,Yij); (II.1) где Zi - значение критерия или параметра, интегрально характеризующего протекание рабочего процесса в камере сгорания (например, коэффициент полноты сгорания топлива, индекс эмиссии загрязняющего вещества, коэффициент избытка воздуха при воспламенении топливовоздушной смеси или прекращении ее горения, максимальная величина температуры стенки и т.п.); Xij - значение критериев и параметров, характеризующих режим работы двигателя; Yij - значение критериев и параметров, характеризующих геометрию проточной части камеры сгорания.Знаки неравенств будут использованы в дальнейшем для исследования нестационарных процессов воспламенения топливовоздушной смеси, прекращения ее горения и при исследовании влияния стабильности течения газа в диффузоре камеры сгорания на стабильность величины коэффициента максимальной окружной неравномерности поля температуры газа на входе в сопловой аппарат турбины (знак равенства используется для описания характеристик стационарного процесса).Назовем систему функциональных зависимостей (II. 1), преобразованных к определенным математическим зависимостям, позволяющим рассчитать величины Zi для любой и конструкции и любого режима работы - «интегральной математической моделью камеры сгорания».Х- характерный размер реакционного объема, источников или стоков энергии и вещества; W - интенсивность закрутки потока газа; G^p - отношение массовых расходов примеси (топлива) в струе к расходу газа в струе.По мнению автора, функция этих критериев позволит описать процесс формирования характеристик рабочего процесса в камере сгорания, в которой он определяется не только скоростью протекания химических реакций окисления горючего, но и процессами испарения топлива (для жидкого топлива), процессами смешения топливных и воздушных, горячих и холодных потоков - т.е. тепломассообменными процессами.В следующем разделе, представлен анализ вышеназванных критериев для конструкции камеры сгорания газотурбинного двигателя. III. Применение разработанных интегральных критериев и параметров к изучению и описанию характеристик процессов в камере сгорания газотурбинной установки. III.1. Принципиальная схема камеры сгорания.Принципиальная конструктивная схема камеры сгорания газотурбинного двигателя представлена на рис.3.1. Воздух из компрессора по специально спрофилированным каналам поступает внутрь замкнутого объема, образованного стенками обечайки, называемого «жаровой» или «пламенной» трубой. Внутри этого объема происходит процесс химического реагирования топлива и воздуха (горение). За счет превращения химической энергии топлива в тепло, температура воздуха повышается, который, выходя из камеры сгорания, совершает работу в турбине и выходном устройстве двигателя.Топливо Воздух I • ^ t Газ Рис.3.1. Принципиальная схема камеры сгорания. III.2. Характерные размеры и безразмерные геометрические критерии камеры сгорания.Проанализируем критерии и параметры2", 1р, 'SF^ Ч ^ т Л ,l,W,G„p, входящие в функцию (1.17) (см.Приложение 1 — раздел 1.7) для конструкции проточной части камеры сгорания газотурбинного двигателя.8) dj^ - медианный диаметр капель распыливаемого форсункой топлива.Интенсивность закрутки газового потока в жаровой трубе камеры сгорания будем характеризовать отношением величины тангенциальной W составляющей скорости к осевой - — -^в кольцевой закрученной струе на выходе завихрителя фронтового устройства (рис.3.4).Условие сохранения тангенциальной составляющей скорости в закрученной струе: = const. или: W,R = const, (III.l) где W, - тангенциальная составляющая скорости, R -текущий радиус.Тангенциальная составляющая скорости может быть выражена через скорость на выходе из межлопаточных каналов завихрителя W^^ (рис.3.4): W,=W,sm(p^, (III.2) К где Rj^ - радиус, проходящий через центры тяжести межлопаточных каналов, ^ -угол наклона лопаток к оси на выходе из завихрителя (рисЗ.4).Из уравнения неразрывности: К=^^ (Ш.З) где Fj - объемный расход воздуха через завихритель, Fj^ - площадь на выходе из межлопаточных каналов завихрителя.Объемный расход воздуха через завихритель равен: V,=M^Rs"J^^, (Ш.4) где // - коэффициент расхода через завихритель, R^ -наружный радиус завихрителя (рис.3.4), АР -перепад давления на отверстиях жаровой трубы.Подставим (Ш.З) с учетом (III.4) в (III.2): W, = M—^sm(p-^ . (III.5) TuR.Rj, .Комплекс —^—^sin^- есть геометрическая характеристика завихрителя, обозначаемая А и выведенная таким же образом для центробежной форсунки в [1,4].Геометрическая характеристика завихрителя является конструктивным критерием подобия.Площадь поперечного сечения турбулентной струи с примесью обозначим через F^, где F^ - определим как элемент площади поперечного сечения жаровой трубы, на который приходится расход топлива через форсунку.К объемам реакционных зон камеры сгорания отнесены следующие: - объемы, в которых происходит формирование газодинамических параметров потока воздуха перед поступлением его внутрь объема жаровой трубы; для камеры сгорания таким объемом является объем диффузора - Уд (подробнее об особенностях течения газа в объеме диффузора- см. раздел о нестационарных процессах в камере сгорания) и объем кольцевого канала камеры сгорания; для этих объемов массовая скорость химической реакции —Wp=0; - объем жаровой трубы, где происходит процесс химического реагирования топлива и воздуха с температурой, имеющей величину >1800К; таким объемом является объем части жаровой трубы до последнего ряда отверстий в ее боковых стенках - Узг ', для этого объема характерным является максимальное значение скорости протекания химических реакций, максимальная интенсивность образования NO2 и выделения углерода в местных «переобогащенных» топливом зонах; - весь объем жаровой трубы до сечения входа в сопловой аппарат турбины — Уж', во всем объеме жаровой трубы происходят процессы окисления топлива и продуктов его недожога до конечных продуктов реакции - СО2; Н2О; NOx (iVO^ является продуктом окисления N02)\ объем газосборника жаровой трубы, равный разности величин ¥ж Узг ; в этом объеме происходят процессы окисления углерода, образовавшегося в местных «переобогащенных» топливом зонах Узг и процессы выравнивания температурной неравномерности по длине объема газосборника. ( » IIL4. О минимальном поперечном размере камеры сгорания.В настоящее время проявляется большой интерес, как в авиации, так и в энергетике к двигателям малой размерности, включая микродвигатели. Для двигателей малой размерности характерно применение центробежных компрессоров с противоточной камерой сгорания, расположенной над турбиной.Анализ показывает, что противоточные камеры сгорания малоразмерных двигателей обеспечивают более равномерное поле температуры газа при более низких потерях по сравнению с прямоточными камерами сгорания. Это объясняется большой конфузорностью и длиной газосборника жаровой трубы и малыми значениями приведенной скорости Х^ за компрессором.В камерах сгорания двигателей малой размерности существенно уменьшается высота жаровой трубы, что однозначно приводит к увеличению площади боковой поверхности жаровой трубы 8ж относительно ее объема. Как показывает опыт, увеличение относительной площади поверхности жаровой трубы отрицательно влияет на тепловое состояние стенки жаровой трубы А^ =Т^-Т^ (Г -максимальная температура стенки жаровой трубы) и может W »' приводить к снижению эффективности горения.Как следует из графика, наибольшие трудности возникают с охлаждением малоразмерных, особенно противоточных камер сгорания с большим отношением поверхности жаровой трубы к ее объему.Изменение знака градиента функции А/ = / (изображенной на графике рис. 3.6), начиная с величины —^ = 85м , связано с уменьшением ж W> коэффициента полноты тепловыделения в камере сгорания.Исследования, подтвержденные опытами, включая эксперименты других авторов, показывают, что при этом, из-за снижения температуры газа в реакционной зоне, в связи с интенсивным отводом тепла в стенку, уменьшается скорость химических реакций и увеличивается недожог топлива.Так как максимальная температура стенки жаровой трубы определяется допустимым значением для материала и примерно одинакова для всех камер, то из формулы (III. 19) следует, что для двигателей малой размерности существует ограничение по величине температуры воздуха за компрессором (т.е. по степени повышения давления) из условий теплового состояния стенки.Для камеры сгорания -^^ {h^ )'^.На графике рис. 3.7 показана зависимость высоты жаровой трубы от величины отношения —^, полученная статистической обработкой результатов расчета отношений площади поверхности жаровой трубы к ее объему для различных конструкций кольцевых камер сгорания.Из графика на рис.3.6 и 3.7 следует, что минимально допустимым значением /гж для кольцевой камеры сгорания является {h^)mm 0,045 м .Указанные выше ограничения касаются камер сгорания с металлическим материалом стенок жаровой трубы. Как следует из приведенного анализа, применение керамических материалов, уменьшающих потери тепла из зоны реакции, позволит расширить минимальное ограничение поперечного размера жаровой трубы. I1I.5. Экспериментально-теоретическое обоснование выбора величины длины части жаровой трубы до отверстий смесителя.При проектировании жаровой трубы камеры сгорания важным является правильный выбор расстояния, на котором расположены отверстия для подвода воздуха в зону смешения. С одной стороны, слишком близкое расположение отверстий к фронтовому устройству, т.е. к зоне стабилизации пламени может привести к быстрому «замораживанию» химических реакций окисления топлива. С другой стороны, объем реакционной зоны образования оксидов азота возрастает с увеличением расстояния от торцев форсунок до отверстий подвода воздуха в зону смешения, что приводит к повышению концентрации оксидов азота в продуктах сгорания.Процесс выгорания топлива в жаровой трубе камеры сгорания определяется аэродинамической структурой потока, формируемой конструктивными элементами ее проточной части, например, завихрителем фронтового устройства, создающим зону стабилизации пламени, где происходит наиболее интенсивное протекание химических реакций и процессов тепломассообмена между закрученной струей воздуха, горячими продуктами сгорания (в зоне рециркуляции) и топливом. Аэродинамическую структуру потока закрученной струи воздуха и обратного тока рассмотрим с точки зрения теории смешения турбулентных струй.Соотношение (III.23) показывает, что длину зоны смешения закрученных струй прямого и обратного токов определяют интенсивность крутки W и отношение —-.К На рисунке 3.10 изображена зависимость относительной длины смешения закрученной струи воздуха с газами в зоне обратного тока от интенсивности крутки воздушного потока, рассчитанная по формуле (III.23) для —^ = 0,38 (среднестатистическая величина) и —- = 0,19 (минимальная величина из ^ж реализованных конструкций). Коэффициент А, входяп).ий в формулу (III.23) и характеризующий турбулентность потока, выбран в соответствии с рекомендациями в [52].На рисунке 3.11 изображен график, построенный на основании данных экспериментальных исследований, где по оси ординат отложены значения величины недожога топлива 1 - т; отнесенные к ее минимальному значению. По оси абсцисс - величина относительной длины зоны смешения струй прямого и обратных токов, рассчитанная по (III.23) . Из рисунка видно, что минимальный недожог топлива достигается при — = 0,8 ч-1,0. Величина этой относительной К длины является минимально допустимой для расположения отверстий в боковых стенках жаровой трубы.1 - форсунка, 2 - завихритель, 3 - камера предварительного смешения.Н^{м) 0.3 0,25 0,2 0.15 0.1 0,05 О О 10 20 30 40 50 60 70 80 90 100 S яс яс (м-) Рис 3.7. Зависимость высоты жаровой трубы от величины площади ее относительной поверхности для камеры сгорания газотурбинного двигателя.Рис.3.8. К расчету длины зоны смешения закрученного завихрителем потока воздуха с продуктами сгорания рециркуляционной зоны.Рис.3.9. Схема взаимодействия зон прямого и обратного тока в завихрительном фронтовом устройстве. i w Рис.3.10. Зависимость относительной длины смешения струй прямого и обратного тока от интенсивности крутки потока воздуха на выходе из завихрителя. 1 ^=0,38, 2 —^=0,19 Ж Ж.

1-V 3.5 2.5 1.5 0.5 0.4 0.6 0.8 h ж Рис.3.11. Зависимость величины недожога топлива от относительной длины зоны смешения струй прямого и циркуляционного потока во фронтовом устройстве. IV. Интегральный подход к математическому моделированию камеры егорания. IV.l. Представление о рабочем процессе в камере сгорания.Рабочий процесс в камере сгорания отличается многообразием одновременно протекающих и взаимосвязанных физико-химических превращений в турбулентном потоке газа со сложной аэродинамической структурой. Рисунок 4.1 графически иллюстрирует представления о физических явлениях, происходящих в камере сгорания. Изображенная на рисунке схема является моделью физических процессов в камере сгорания. Модель построена на основе теоретического анализа явлений с помощью разработанного в предыдущих разделах интегрального подхода и полученных в результате критериев и параметров.Критерии и параметры определяют характеристики потока сплошной среды на выходе из рассматриваемого объема, с учетом интенсивности внутренних молекулярных превращений, молекулярных и молярных процессов обмена энергией, веществом, количеством движения, процессов выравнивания неоднородностей среды в пространстве этого объема и скорости переноса этих процессов потоком сплошной среды от входа к выходу из объема камеры сгорания (рис.4.1). IV.2. Интегральная математическая модель камеры сгорания.На основании изложенного в Приложении 1 материала по преобразованию дифференциальных уравнений переноса энергии, концентрации и количества движения в неразрывном потоке сплошной среды с развитыми турбулентными пульсациями при наличии стоков и истоков теплоты и массы и анализа полученных критериев применительно к камере сгорания газотурбинной установки, представленного в разделе III, определены основные критерии и параметры, определяющие формирование характеристик рабочего процесса в камере сгорания.Величины температуры и давления определяют скорость химического реагирования веществ в потоке газа и образования конечного продукта реакции, так как она зависит от скорости движения молекул реагирующих веществ и частоты их соударений. Скорость движения молекул пропорциональна квадратному корню значения температуры, а длина свободного пробега молекул обратно пропорциональна величине давления газа. Таким образом, скорость движения молекул и частота их соударений определяются значениями температуры и давления газа.Физико-химические параметры: вид топлива и коэффициент избытка окислителя определяют выделяющуюся при химическом реагировании окислителя и горючего энергию. От вида топлива также зависит скорость его дробления на капли, испарение, диффузия в движущемся потоке газа, т.е. местные величины концентраций, определяющих местное значение скорости протекания химической реакции.Гидродинамический параметр G/Vp характеризует скорость конвективноГО переноса энергии, вещества и количества движения. Размерность параметра такая же как и у массовой скорости образования или расходования вещества в результате химической реакции в единице, объема. С ростом величины G/Vp возрастает «вынос» из реакционного объема энергии, вещества и количества движения, что приводит к снижению температуры, концентрации конечного продукта реакции и росту величины потерь полного давления в этом объеме.Универсальность параметра подтверждается практикой успешного обобщения экспериментальных данных по характеристикам горения как в России, так и за рубежом с помощью параметров, содержащих G/Vp (критерий Дорошенко, Лонгвелла, Лефевра и др.). " ' F Степень расширения диффузора камеры сгорания ^^ характеризует стабильность течения потока воздуха на входе в реакционный объем жаровой трубы. При большой степени расширения и малой длине диффузора могут возникать отрывы пограничного слоя с последующим вихреобразованием попеременно то от внутренней, то от наружной стенок, в разных местах по окружности кольцевого диффузора. Отрыв пограничного слоя от стенок происходит под воздействием положительного градиента давлений.Такая нестабильность течения воздуха является причиной значительного изменения температуры стенок жаровой трубы и неравномерности поля температуры газа на входе в сопловой аппарат турбины. F Практика показала [39], что при —^ < 2 удается обеспечить стабильное протекание процессов в камере сгорания, независимо от того безотрывное или отрывное течение воздуха в диффузоре (подробнее о процессе нестабильного течения газа в диффузоре - см. раздел о нестационарных процессах в камере сгорания) Критерии ^— и ^— отражают протекание в камере сгорания таких процессов как проникновение струи в сносящий поток газа с учетом угла наклона этой струи к направлению сносящего потока. F Параметр ^—характеризует аэродинамику процесса смешения холодных струй с горячим потоком и определяет снижение местных максимальных температур газа. Для примера рассмотрим процесс втекания одиночной струи воздуха перпендикулярно сносящему Потоку горячего газа (рисунок 4.2).Обозначим: y = h -глубина проникновения струи воздуха в сносящий поток на длине х = 1 (Рис.4.2); VQ - скорость воздуха в отверстиях жаровой трубы; Mf - скорость потока газа в жаровой трубе.Глубина проникновения струи в сносящий поток газа зависит от отношения скоростей потока в струях, которое определяет параметр — , так как ож F V '•ж ^0 Запишем выражение для потока импульса через поперечное сечение втекающей в объем жаровой трубы через ее бо1?овую стенку струи в виде: где v,d- текущая скорость на оси струи и текущий диаметр; х^ - координата, отсчитываемая вдоль оси струи.Примем, что вертикальная составляющая импульса постоянна вдоль струи: v^f{x^ ) COS 0 = const; 2 w \2 const V fM =—^; COS^ где ^-угол между касательной к траектории оси и вертикалью.После интегрирования получаем уравнение траектории струи: Из этого выражения следует, что отклонение оси втекающей струи воздуха определяется ее начальным диаметром и отношением скоростных напоров ри потоков —Y. Величина угла 9 определяет глубину проникновения струи в сносящий поток (Рис.4.2). Увеличение скоростного напора втекающей струи воздуха уменьшает угол в, а следовательно увеличивает глубину ее проникновения в набегающий поток.Проанализируем, как в случае параллельного втекания струи воздуха в горячий поток газа (истечение воздуха через щели охлаждения жаровой трубы) F параметр •=-^— влияет на процессы смешения. Рассмотрим истечение плоской, ограниченной с одной стороны струи воздуха, имеющей начальную скорость UQ В поток горячего газа, движущегося со скоростью щ в ту же сторону, что и воздух (рисунок 4.3).Обозначим:Z)o -начальная ширина струи воздуха; Н -начальная ширина горячего потока газа; b^ -текуш;ая координата зоны смешения струи воздуха с потоком газа, отсчитываемая от плоскости начального разделения потоков (рисунок 4.3); / - координата внешней границы зоны смешения, отсчитываемая вдоль оси х.Параметр W характеризует интенсивность процессов смешения, определяющих уровень местной температуры в потоке газа. Значительно улучшить смесеобразование и смешение потоков на меньшей длине в зоне горения позволяют закрученные струи, формируемые завихрителем фронтового устройства.Течение в закрученной струе можно разделить на область прямого и обратного токов, между которыми происходит интенсивный турбулентный обмен. Область обратных токов несет к фронтовому устройству горячие продукты сгорания и обеспечивает стабилизацию пламени.При увеличении параметра W происходит более интенсивное выравниваif ние концентраций реагентов и температуры в зоне взаимодействия прямого и обратного тока в кольцевой закрученной струе воздуха.Для камеры сгорания газотурбинного двигателя характерно комбинированное использование закрученных (завихритель фронтового устройства), поперечных (отверстия для подачи воздуха в стенках жаровой трубы), параллельных (щели пленочного охлаждения стенок жаровой трубы) струй.Параметр —— характеризует равномерность распределения концентрации топлива по всему сечению жаровой трубы. При бесконечно большом количестве форсунок топливовоздушная смесь становится однородной.Проанализируем влияние отдельных конструктивных особенностей горелочных устройств (воздушный завихритель и топливная форсунка) камеры сгорания на процессы смешения топлива с воздухом в ее реакционных объемах.Значения критериев Вебера, Лапласа, критерия М при 5^^ = const на установившихся режимах работы двигателя для конкретной конструкции форсунки зависят только от термодинамических параметров - Р}^,Тк,Т, а и вида топлива.На основании выведенной автором критериально-параметрической функции и анализа конструктивных схем камер сгорания: 1. Обоснованы критерии и параметры, определяющие рабочий процесс формирования характеристик потока газа в камере сгорания газотурбинного двигателя.2. Произведен анализ физической сущности каждого из критериев и параметров камеры сгорания.3. Представлена схема физического процесса в камере сгорания газотурбинного двигателя, иллюстрирующая представление о происходящих в ней физико-химических явлениях.4. Определены «граничные условия» - минимально допустимая высота жаровой трубы для металлического материала стенок и минимально допустимая длина от завихрителя воздуха фронтового устройства до отверстий поперечного подвода струй воздуха в объем жаровой трубы.Рис.4.1. Модель физического процесса в камере сгорания.Явления пространственного переноса массы, тепла, количества движения в потоке с горением.Парат.тетры кокструицш ВХОД КАМЕРА СГ0РАШ1Я 6F/Fi^;W; Гидродикат.шческий критерии и параг.тйтры ^1 %2G/f;(Sc^r) G/r Тер1<юдпкз1,шческпе и тегатофизическиа парат.тетры К К а Б11Д топлива Процессы пространстненкого турбулентного ст.гешеши ПОТОИОЕ (т.годепь зурбупентиостк для KSJitepH сгорания'! турбулентная гш44|узия1 турбулентная ЕДКОСТЬ турбулентная теплопроводность Дробление капель топшаа Молекулярная диффузия Молекулярная БЯЗВОСТЪ Молекулярная теплопроводность Скорость кга11ИЧйСЕ0Й реакции Испарение Влияние на скорость протекания т.хакрои т.шкропроцессоЕ неравколтерности распределения napaiiteTpoE по o6-bei,iaiit каракгернык зонкаг.теры сгорания Скорость протекания т.такропроцессоЕ Скорость протекания т.тюфопроцессоЕ ВЫХОД Характеристики рабочего процесса Б каг.хере сгорания Скорость переноса т,1атерш{ Рис.4.2. Схема втекания струи воздуха перпендикулярно сносящему потоку.Экспериментальное исследование характеристик рабочего процесса в камере сгорания проводилось на стендах при автономных испытаниях в составе специальных установок для испытаний как моделей, так и полноразмерных камер сгорания.На стендах для автономных испытаний исследовались такие характеристики как неравномерность поля температуры газа на выходе из камеры сгорания, величина коэффициента полноты выгорания топлива в ней, коэффициент восстановления полного давления (при «холодной продувке»), коэффициент избытка воздуха при «бедном» срыве пламени, эмиссионные характеристики и характеристики запуска в различных высотно-климатических условиях (термобарокамера).На рис.5.1 представлена схема стенда с установкой для автономных испытаний камер сгорания и их моделей, на рис.5.2 - схема термобарокамеры, позволяющей имитировать различные высотно-климатические условия с установкой для исследований нестационарных процессов воспламенения топливовоздушной смеси в камере сгорания, Стендовая система включает в себя: - устройство для подвода или отвода (эксгаустер барокамеры) воздуха — компрессор, баллоны высокого давления, при испытаниях моделей также используется воздух от заводской сети; (# - газовый или электрический подогреватель воздуха; - рессивер; - выхлопная шахта.Установка для испытаний камер сгорания состоит из следующих конструктивных элементов: - входного устройства, с плавным входом и участком для измерения параметров потока газа на входе в камеру сгорания, стыкуемого к фланцу модуля камеры сгорания; - испытываемого модуля камеры сгорания; - выходного устройства с установленными в нем (в зависимости от цели проводимых работ) средствами измерений - термопарами, насадками для измерения полного давления, пробоотборниками.На двигательных стендах исследовались такие характеристики камеры сгорания как неравномерность поля температуры на входе в сопловой аппарат турбины и эмиссионные характеристики.В данной работе использованы также данные по измерению эмиссии, которое проводилось в составе объекта. V.2. Методика проведения испытаний.Методика проведения автономных испытаний камеры сгорания обеспечивает определение всех характеристик рабочего процесса в ней.В процессе испытаний производится измерение следующих параметров: - на входе в камеру сгорания - величин давления, температуры, массового расхода топлива, определение массового расхода воздуха; - на выходе из камеры сгорания - величин температуры, давления потока и химического состава газа в различных точках поперечного сечения (многоточечные гребенки термопар, насадки полного давления, пробоотборники, устанавливаемые как стационарно, так и с организацией непрерывного перемещения в плоскости поперечного сечения камеры сгорания).При исследовании нестационарных процессов - запуска и прекращения горения фиксируются параметры на входе в камеру сгорания в момент воспламенения или погасания топливовоздушнои смеси в камере сгорания.Момент воспламенения или погасания определяется по показаниям термопар в выходном сечении камеры сгорания. V.3. Максимальные величины погрешностей при определении параметров в процессе испытаний камеры сгорания.Максимальные погрешности измеряемых параметров в процессе проведения испытаний камеры сгорания представлены в таблице I.Таблица I.Погрешности измерения давления газа Погрешности измерения температуры Погрешности измерения массового расхода жидкого топлива Погрешности определения массового расхода газа Погрешности измерения концентрации в продуктах сгорания конечных и промежуточных продуктов реакции NO^,CO^,H^O,HC,C ±0,5% измеряемого значения ±1% верхнего предела измеряемого значения ±0,5% измеряемого значения ±0,7% расчетного значения в соответствии с ГОСТ 17.22.04-86 V.4. Метод обработки экспериментальных данных.Удобство такого способа заключается в том, что он дает возможность учесть вид функции (линейная, квадратичная, степенная, экспоненциальная и т.п.) влияния каждого фактора с точки зрения его физической сущности на протекание процесса.Для выявления вида функции влияния и степени значимости каждого фактора использованы способы: 1) графического представления экспериментального материала построением зависимостей в координатах: значение характеристики камеры сгорания - [величина критерия, определяющего характеристику] 2) метод «случайного баланса», в котором производится оценка величины изменения характеристики камеры сгорания при изменении каждого из критериев от минимального до максимального значения, реализуемого в ходе эксперимента. # стендовая магистраль подвода воздуха к камере сгорания установка средств для измерения характеристик рабочего процесса на выходе из камеры сгорания, газовый или электрический подогреватель рессивер измерение параметров потока газа на входе в камеру сгорания выхлопная шахта стенда Рис5.1. Схема стенда с установкой для автономных испытаний камер сгорания.Измерение давления и температу1зы возд^оса в термобарокамере.Измерение расхода Подвод воздз'ха через камеру топлива сгорания Объект испытаний Измерение температуры Камера сгорания газа на выходе из камеры сгорания Подвод BosfljTca от т>рбодетандера Измерение давления топливоподачи н температ>фы топлива Рис.5.2. Схема стенда с установкой для проверки запуска камеры сгорания. «' VI. Характеристики камеры сгорания как функции ее критериев. VI.1. Эмиссия оксидов азота. Автором выполнен большой объем исследований камер сгорания с фронтовыми устройствами, охватывающими широкий спектр типовых и перспективных конструкций (завихрительные, аэрационные, микрофакельные, комбинированные) как на газе, так и на жидком топливе, с высоко - и низконапорной подачей топлива. Схема испытанных фронтовых устройств представлена на рис,6.1.Исследования проводились на моделях и полноразмерных камерах сгорания в условиях автономных испытаний, испытаний двигателей на стендах, на летательных аппаратах и на автономных газотурбинных электростанциях. Испытано более 30-и вариантов фронтовых устройств (включая серийные конструкции) и проведено более 300 экспериментов в рабочем диапазоне режимов двигателей: Рк^0,4 - 2,1 МПа; 7^=300 - 820 К; Гг=600 - 1700 К; Т^ - температура газа на выходе из камеры сгорания.Обработкой экспериментального материала получено: для жидкого топлива (керосин, дизельное топливо): Е1,,^=943,6Р, ( G У Т ( Т ехр зг J 288.Размерности в формулах (VI. 1) и (VL2): Давление - МПа; температура - К; площадь - м ;объем - м ; массовый расход — кг /сек.Для образования NOx характерным является объем зоны горения, т.е. область высокой температуры газа. При этом 'длина зоны горения равна расстоянию от среза топливной форсунки до отверстий смесителя, где происходит «замораживание» химических реакций образования NOxТочность расчета по приведенным формулам представлена на рис.6.2, откуда следует удовлетворительная сходимость расчетных и экспериментальных данных. VI.2. Эмиссия оксидов углерода.Индекс эмиссии бензо-(«)-пирена связан с индексом эмиссии несгоревших углеводородов следующей зависимостью: Сопоставление расчетных и экспериментальных данных по эмиссии несгоревших углеводородов представлено на рис.6.3. VI.4. Уровень дымления камеры сгорания. Уровень дымления газов на выходе из камеры сгорания принято оценивать по величине уменьшения отражательной способности бумажного фильтра (в %) при пропускании через него пробы газа определенного объема.Объемом реакционной зоны, определяющей содержание сажи в газах на выходе из камеры сгорания, является объем газосборника жаровой трубы. Насыщение газа кислородом с выгоранием чабтиц углерода происходит в объеме газосборника, равном разности объемов жаровой трубы и зоны горения.При этом функции Рк? Тк и Гг учитывают как химические процессы окисления сажи в газосборнике жаровой трубы, так и физические процессы ее образования в зоне горения. f Точность расчета по приведенной формуле иллюстрируется графиком на рис.6.4. VI.5. Взаимосвязь между единицами, характеризующими уровень токсичности газов на выходе из камеры сгорания.Взаимосвязь размерностей, выражающих количество загрязняющего вещества в миллионных частях массы пробы [ррт] и концентрации в [мг/м^] описывается соотношениями: NO^ [ppm] = — • СО [ppm] = ЯС[ррт] = 2,05 СО [мг/м'] 1,25 НС [мг/м'] ' 0,715 f VI.6. Коэффициент полноты выгорания топлива.После преобразований, учитывая, что GQIGJ мало, выражение для описания связи полноты сгорания топлива с индексами эмиссии СО и НС можно представить в виде: \-т] = {Е1,,+0,2Ъ2Е1,,)\0-\ В такой форме это уравнение приведено в [46]. VI.7. Зависимость величин коэффициентов неравномерности поля температуры газа на выходе из камеры сгорания от ее критериев.Повышенная величина коэффициента окружной неравномерности поля температур газа приводит к перегреву или прогару лопаток соплового аппарата, а эпюра средней температуры по высоте канала определяет работоспособность лопаток турбины.Точность описания экспериментальных данных по измерению максимального значения температуры стенки жаровой трубы для камер сгорания разI личных конструкций, аппроксимирз^ощей их зависимостью, иллюстрируется графиком на рисунке 6.7.Выводы к разделуУ!.Разработана интегральная математическая модель характеристик стационарного процесса в камере сгорания газотурбинного двигателя, позволяющая произвести расчет, в зависимости от режима работы и геометрии проточной части камеры сгорания, значений следующих характеристик рабочего процесса: - коэффициента полноты сгорания топлива; - индексов эмиссии всех составляющих вредных выбросов (оксидов азота, монооксидов углерода, несгоревших углеводородов, бензопирена) и параметра выброса сажи; - коэффициентов окружной и радиальной неравномерности поля температуры газа; - коэффициента потерь полного давления в камере сгорания; - максимальнз^ю величину температуру стенки жаровой трубы камеры сгорания. б) в) г) Рис.6.1. Конструктивные схемы фронтовых устройств.Двигатели, не участвуюндие в обобщении: Ь^ - ГТН-25 (газовое топливо); О - ГТУ «Сименс», модель V 94.2 (газовое топливо); О НК-86; Т - GF6-50; Л' - модель двухзонной камеры сгорания с «богато-бедным» горением.Рис.6.2. Сопоставление расчетных и экспериментальных данных по эмиссии рксидов азота. .t EL CO г/кг.топл. (расч.) ь * * * • 40 80 120 EI. CO г/кг.топл. (эксп.) EI f^C г/кг.топл. (расч.) 0.1 ^ (^ • П . . ^ 0,1 10 100 £:/ ^c г/кг.топл. (эксп.) Рис. 6.3. Сопоставление расчетных и экспериментальных данных по эмиссии оксидов углерода и несгоревших углеводородов, (обозначения точек на графиках соответствуют Рис 6.2) SN(%) расч.ЛГ\ . ни 'jn n - к' л*' ,0 20 40 SN(%) ЭКСП. Рис.6.4. Сопоставление расчетных и экспериментальных данных по числу дымности отработавших газов двигателя. (обозначения точек на графике соответствуют Рис 6.2.) Цр Й W расч.0.3 0.25 0.2 к.у^ у^ р% 0,2 0.25 0.3 ^шах ^КСП. Рис. 6.5. Сопоставление расчетных и экспериментальных данных по определению коэффициента окружной неравномерности поля температуры газа для разных конструкций камер сгорания.2^ \ дСг дСг дСг дСг / ^ , ^ Jd'Cr д'Сг д'Сг dt дх ду dz \^dx ду dz j + dx dx dy dy dz dz где Cj- - концентрация горючего.Процесс воспламенения топливовоздушной смеси в камере сгорания характеризуется резким повышением температуры за короткий промежуток времени и уменьшением концентрации горючего из-за превращения его в # конечные продукты реакции за тот же промежуток времени: dt dt <0 .Условие непрерывности протекания химической реакции при воспламенении топлива, определяющееся превышением скорости поступления горючей смеси с потоком газа в реакционный объем над скоростью ее выгорания в этом объеме выражается неравенством (VIL7).Важнейшей задачей при создании камеры сгорания газотурбинного двигателя является обеспечение надежного воспламенения топливовоздушной смеси в различных высотно-климатических условиях его эксплуатации.Особенно актуально решение проблем запуска для авиационных двигателей, стационарных и транспортных установок, работающих в высокогорных условиях при отрицательной температуре воздуха и топлива.В настоящее время доводка камеры сгорания по обеспечению запуска осуществляется дорогостоящим методом проб и ошибок с привлечением специального оборудования (термобарокамеры) и испытания большого количества вариантов.Для описания процесса воспламенения топливовоздушной смеси в камере сгорания ГТД используем систему неравенств (VILS) и (VII.9) Для процесса запуска, осуществляемого системой непосредственного розжига (от свечи зажигания): J, си' где E^jj = энергия разряда, сосредоточенная в конденсаторе; С - емкость конденсатора; U- напряжение.При нахождении функций, входящих в неравенства, использованы многочисленные экспериментальные материалы по исследованию полноразмерных камер сгорания и моделей на натурных режимах на автономном стенде в термобарокамере.Схема стенда представлена на рисунке 5.2.Объектами испытаний явились камеры сгорания экспериментальных и серийных двигателей: прямоточные (ТВЗ-117, РД-33 и др.), противоточные (ГТД-1250, ТВ7-117, стационарный двигатель для электростанций ГТЭ-1,5 и др.), одиночные (ГТД-350) с завихрительными, аэрационными и микрофакельными фронтовыми устройствами, с вращающейся форсункой.Конструктивные схемы исследованных фронтовых устройств показаны на рис.6.1.Испытания проводились на керосине, дизельном летнем, зимнем и арктическом топливе, без вдува и со вдувом вспомогательного воздуха во второй контур форсунок с перепадом давления 0,02-0,3 МПа. Воспламенение производилось от свечей поверхностного и кумулятивного разряда с частотой разрядных импульсов от 2 до 25Гц, расположенных на расстоянии от торца форсунки 15... 100мм по длине и высоте, с погружением внутрь жаровой трубы на 0...35ММ от ее стенок.Параметр переноса, как видно из анализа неравенств (VII. 12), (VII. 13) и рисунка 7.1, вместе с величиной а определяет область воспламенения топливовоздушной смеси в камере сгорания на запуске двигателя.Неравенство (VII. 12) получено на основе преобразования и решения уравнения энергии в движущемся потоке для нестационарного процесса и отражает условие превышения тепловыделения в результате химических реакций над отводом тепла с продуктами сгорания в движущемся потоке газа.Неравенство (VII. 13) получено решением уравнения концентраций для f нестационарного потока и отражает условие непрерывности протекания химической реакции при запуске, т.е. превышение скорости поступления новых порций горючего в реакционную зону над скоростью выгорания.Величина коэффициента избытка воздуха определяет качество запуска, т.е. температуру газа перед турбиной и надежность воспламенения топливовоздушной смеси в камере сгорания.Запуск без заброса температуры газа на выходе из камеры сгорания происходит при значениях коэффициента избытка воздуха « > 2,5. Если or<2(но не менее а-\\ для предотвращения прогаров концов рабочих лопаток турбины при запуске необходимо обеспечить резкий сброс давления топлива после его воспламенения путем соответствующих мероприятий в системе регулирования двигателя. Анализ неравенств (VII.12), (VII.13) и графика на рисунке 7.1 показывает, что запуск камеры сгорания осуществляется при выполнении двух условий: необходимого (кривая «а») и достаточного (кривая «Ь»).При значениях и максимальном значении df неравенства (VII. 12), ч ^ж J opt (VII.IS) обращаются в тождества (точка пересечения кривых «а» и «Ь» на графике рис.7.1). Так как — пропорционально частоте вращения ротора турбокомпрессора на пусковых режимах, то точка пересечения кривых на рис,7.1 определяет оптимальную частоту вращения ротора, при которой происходит воспламенение топливовоздуп1ной смеси в камере сгорания с минимальной величиной температуры газа перед турбиной.Размерности в (VII. 12) и (VII. 13): давление - в МПа, температура - в К, расход воздуха - в кг/сек, объем - в м"^ , диаметр капли — в метрах, энергия разряда - в джоулях.Неравенства (VII. 12) и (VII. 13) представляют собой математическую модель процесса воспламенения топливовоздушной смеси в камере сгорания в различных высотно - климатических условиях.Достоверность расчета характеристик запуска подтверждается сравнением с опытными данными, что иллюстрируется на рисунке 7.2.Неравенства (VII. 12) и (VII. 13) апробированы в следующем диапазоне режимных и конструктивных параметров: Рк = 0,0227...0,1 МПа (Н=0...12км); Тк = 223...296К; G=0,181...2,6 кг/сек; ci:=1...15; Есв=0,3...8Дж; Vж=0,008...0,0405м^; /zF^ =0,044...0,44; F^//zF^^ =5,61...11,4; число ( топливных форсунок 1...80; миделевая скорость изменялась в пределах 4,7... 12 м/сек.Влияние режимных и конструктивных параметров на область воспламения топливовоздушной смеси в камере сгорания иллюстрируют графики на рисунках 7.3 - 7.9. Из рисунков видно, что повышение температуры И давления воздуха, увеличение параметра — и энергии разряда свечи ' Z/^ож расширяют область воспламенения в камере сгорания, а увеличение диаметра капель распыленного топлива и числа форсунок, приходящихся на единицу площади, наоборот, сужают эту область. Из графика на рис.7.7 видно наличие оптимальной величины относительной пропускной способности завихрителей jiiFy, при которой а максимально.Анализируя неравенства, отметим, что конструктивные параметры /^^и F ^ для обеспечения запуска не могут изменяться в широких пределах, так ож как они имеют оптимальные значения из условий обеспечения минимальной неравномерности поля температур и минимальных выбросов оксидов азота, монооксидов углерода, несгоревших углеводородов и минимального недожога топлива на рабочих режимах (см. раздел VI).Основными параметрами, с помощью которых могут быть достигнуты высокие показатели высотности запуска, являются энергия разряда Е^^ и медианный диаметр капель распыливания топлива d^^. Максимальное значение Есв на серийных отечественных двигателях достигает 6 Джоулей. За рубежом эта величина составляет 16 Джоулей. Повышение энергии разряда свечи является существенным резервом улучшения запуска двигателей.Эффективным средством снижения диаметра капли d^^ является распыливание топлива с помощью воздушного потока. Так, подача воздуха от воздушного стартера во второй контур центробежной топливной форсунки уменьшает диаметр капли в 10 раз при перепаде давления 0,15 МПа и температуре топлива 233К. На рисунке 7.10 иллюстрируются возможности улучшения высотности запуска на примере современного двигателя для истребителя. По оси ординат отложена высота полета, а по оси абсцисс - варианты конструктивного исполнения.Как видно из диаграммы, запуск двигателя без кислорода осуществляется до 8 км. Применение воздушного распыливания повышает высотность запуска до 13,5 км. Дополнительное использование свечи с знергией разряда 16 джоулей повышает высоту запуска до 18 км.Запуск с подводом кислорода реализуется на серийном двигателе до 11 км, а с применением аналогичных мероприятий высота запуска составляет 22 км.При обобщении результатов экспериментальных исследований по воспламенению топливовоздушнои смеси в камере сгорания величина медианного диаметра капель распыленного жидкого топлива d^^ рассчитывалась с помощью следующих зависимостей [45,73,74,75].1) Распыливание топлива центробежной форсункой.Устойчивость горения в камере сгорания на переходных режимах работы двигателя является одной из основных ее характеристик, определяющей диапазон работы двигательной установки в целом.В литературе отсутствуют данные о влиянии геометрии проточной части камеры сгорания на пределы устойчивого горения в ней, характеризуемых максимальной и минимальной величиной коэффициента избытка воздуха при срыве-аг'"^иа^,„.При обеднении или обогащении смеси горение поддерживается в отдельных очагах за счет уровня местной температуры газа.Размер очага будем характеризовать тепловой шириной зоны горения [28,33,95,105]: S —; СрРи„ где Uj^ -нормальная скорость распространения пламени: и и = /(Рд. ,Т^,Т2,а, вид топлива ).Из анализа неравенств (VII. 14), (VII. 15), (VII. 16) и графика на рисунке 7.17 следует, что срыв горения в камере сгорания осуществляется, как и процесс запуска при выполнении двух условий: необходимого (кривая и линия «а» для «бедной» и «богатой» границы ) и достаточного (кривая «Ь»).Неравенства (VII.14), (VII.15), (VII.16) представляют собой математическую модель для описания устойчивости горения в камере сгорания при обеднении и обогащении топливовоздушной смеси.Влияние режимных и конструктивных параметров на фаницы устойчивого горения (по «бедному» срыву) в камере сгорания иллюстрируется графиками на рисунках 7.20 - 7.22.Из рисунков видно, что увеличение температуры (рис.7.20 а) и давления (рис.7.20 б) воздуха на входе в камеру сгорания расширяют границы «бедного» срыва пламени.Увеличение относительной пропускной способности завихрителя, характеризуемого .величиной параметра; /JF^ (рис7.21а,в), увеличение параметра F^l /JF^^ (рис.7.21 б) и уменьшение закрутки воздушного потока на выходе из завихрителя (рис.7.22) приводит к сужению границ области устойчивого горения в камере сгорания.Главной задачей при создании камеры сгорания газотурбинного двигателя является обеспечение стабильности протекающих в ней процессов.Источником нестабильности может явится диффузор камеры сгорания при неправильном его проектировании. Нестабильность течения потока газа в диффузоре приводит к повышению неравномерности поля температуры газа на выходе из камеры сгорания, причем проявление повышенной неравномерности отличается нестационарным характером.Диффузоры, испытанных камер сгорания газотурбинных двигателей представлены на рисунке 7.23: - кольцевые диффузоры с центральным телом — рис.7.23 а); - кольцевые ступенчатые диффузоры (фиксированный «отрыв») с центральным телом- рис.7.23 б); - криволинейные ступенчатые диффузоры - рис.7.23 в); кольцевые диффузоры - рис.7.23 г). в качестве критерия, характеризующего стабильность характеристик рабочего процесса в камере сгорания, которая зависит от стабильности процесса течения газа в диффузоре, примем величину, равную отношению где Абтах- «разброс» значения коэффициента окружной неравномерности поля температуры газа от экземпляра к экземпляру камеры сгорания с определенной геометрией проточной части (т.е одинаковой конструкции); (0тах)ср - Среднестатистическая величина коэффициента окружной неравномерности поля температуры газа, характерная для камеры сгорания определенной конструкции.График на рисунке 7.24 иллюстрирует зависимость нестабильности характеристик (относительного «разброса» величины окружной неравномерности поля температуры газа) камер сгорания серийно изготавливаемых и опытных двигателей от значения конструктивного параметра диффузора, представляющего собой левую часть неравенства (VII. 17), условно обозначенного Кд и характеризующего условие стабильности течения газа в нем.Из графика следует, что для диффузоров камер сгорания, выполненных с учетом условия (VII. 17) «разброс» характеристик на выходе из камеры сгорания минимален и обусловлен только разбросом геометрических размеров проточной части камер сгорания при их изготовлении в пределах полей допусков, оговоренных чертежом.При отрывном течении газа в диффузоре, т.е. при: Кд > 0,374, значение относительного разброса величины коэффициента окружной неравномерности поля температуры газа на выходе из камеры сгорания, определяется величиной критерия Кд - pHC.7i24. ' Нестабильность течения газа в диффузоре камеры сгорания, приводяпдая к нестабильности поля температуры газа на выходе из камеры сгорания, может явиться причиной дефектов конструктивных элементов турбины, проявляющихся в процессе эксплуатации двигателя.Выводы к разделу VTI.1. Разработана система критериально-параметрических зависимостей в виде неравенств, позволяк^щая производить расчет границ области воспламенения топливовоздушной смеси в камере сгорания и устойчивого горения в ней в зависимости от режимных и конструктивных параметров, 2. Произведен анализ критериев и параметров камеры сгорания на положение границ обдасти воспламенения и устойчивого горения топлива в ней.Энергия разряда свечи 6 Джоулей; запуск на земле (О км).Высота запуска 4 км; температура воздуха и топлива 233К. 10 20 30 G К кг/(с м^) ж Рис.7.3.Влияние на границы области воспламенения топливовоздушной смеси в камере сгорания высоты запуска.Энергия разряда свечи 6 Джоулей; температура воздуха и топлива 233К. Энергия разряда свечи 6 (Джоулей; высота запуска 7км; температура воздуха и топлива 23 ЗК. Энергия разряда свечи 6 Джоулей; высота запуска 4км; температура воздуха и топлива 233 К. Энергия разряда свечи 6 Джоулей; высота запуска 4км; температура воздуха и топлива 233К; jjF^ = 0,16.1 - завихритель; 2 - кольцевая закрученйая струя воздуха; 3 - топливная форсунка; 4 - коническая пелена жидкости на выходе из сопла центробежной форсунки. !^ Рис.7.14. Схема взаимодествия воздушного и топливного потока в «аэрационном» фронтовом устройстве.1 - сопло пневматической форсунки; 2 - жидкостная пленка; 3 воздушный поток. i.Рис.7.18. Сопоставление расчетных и экспериментальных данных по «бедному» срыву. \ ПИП / : эксп > ж й у^ X пУ [ > ^ п 02 04 06 [а ) \ Iran /1 расч Различные обозначения точек на графике соответствуют испытанным вариантам пусковых воспламенителей и моделей камер сгорания.6)l-F^/2/zF^^=6,5;2-F^/X//F,^=ll.154а а {^т^) opt -1 .ГТД-1000; Ж - ГТД-1250; Х-ТВ7-117; А-ТВЗ-117; П-РД-33; О-ТВ2-117; О- ГТД-350; • - камера сгорания экспериментального двигателя. ф ф VIII. Физическое моделирование рабочего процесса в камере сгорания.При работе стенда в камеру сгорания подается подогреваемый воздух.Воздух подогревается в газовом или электрическом подогревателе.Регулирование расхода воздуха осуществляется с помощью крана, установленного перед подогревателем.В выходном сечении камеры сгорания устанавливаются, в зависимости от цели проводимых работ, соответствующие средства измерения (пробоотборники, хромель-алюмелевые термопары, гребенка полного давления).Доказательством применимости данного метода моделирования характеристик рабочего процесса камеры сгорания служит сопоставление результатов испытаний полноразмерных прямоточных и противоточных камер сгорания и их моделей, с одинаковыми величинами их критериев.Это сопоставление представлено на графике, изображенном на рисунке 8.2, где по оси ординат отложены величины характеристик рабочего процесса, полученные при испытании моделей камер сгорания, а по оси абсцисс величины характеристик рабочего процесса, полученные при испытаниях натуральных камер сгорания.Применение простейшей модели камеры сгорания, созданной на основе интегральной математической модели в виде критериальной функции, описанной выше, с учетом возможности пересчета режимных параметров по уравнениям разделаVI, полностью исключает необходимость проведения автономных, дорогостоящих испытаний полноразмерной конструкции на рабочих режимах для исследования всех ее характеристик. BbIZ вых Рис.8.1. Схема одногорелочной модели полноразмерной камеры сгорания.П 0.1 0.01 л к X ( о • 0.01 0.1 10 100 Рис.8.2. Сопоставление характеристик рабочего процесса при испытании физической модели и полноразмерной камеры сгорания. I - Характеристики полноразмерных камер сгорания; II - характеристики камер сгорания, полученные при испытании их моделей.Х,Л,А - индексы эмиссии загрязняющих веществ NOx, СО, НС, г/кг топлива (двигатели ТВ7-117, TB3-1I7, ГТЭ-1,25; керосин, дизельное топливо, метан); • число дымности, %(модификации камер сгорания двигателя РД-33; керосин);О,• — коэффициент окружной неравномерности поля температуры газа (двигатели ТВ7-117, РД-33; керосин, дизельное топливо); • - коэффициент потерь полного давления, % (двигатель РД-33);0 — коэффициент избытка воздуха при «бедном» срыве пламени (двигатель РД-33, керосин.). IX. Система уравнений для расчета характеристик камеры сгорания. IX.1. Система математических зависимостей характеристик камеры сгорания от ее критериев.При проектировании камер сгорания решается задача эффективного превращения химической энергии топлива в тепло с обеспечением устойчивости горения, надежного воспламенения топливовоздушной смеси при низких температурах и давлениях газа. Должны быть рассмотрены проблемы экологичности окружающей среды, неравномерности распределения температур газа на выходе из камеры сгорания и потерь полного давления в ее проточной части.Процессы в камере сгорания газотурбинного двигателя отличаются крайней сложностью в связи с одновременным протеканием в ней физикохимических превращений, таких как распыливание топлива, смешение топлива с окислителем и продуктами сгорания, смешение холодных и горячих струй газа в турбулентном потоке со сложной аэродинамической структурой, которая оказывает влияние на протекание химических реакций окисления топлива в реакционных зонах жаровой трубы. Для их описания необходимо привлечение исследований из области химической кинетики, термодинамики, тепломассообмена, механики жидкости и газа.На основе анализа фундаментальных уравнений тепломассообмена для химически реагирующей смеси в потоке газа, автономных испытаний камер сгорания, исследований моделей и обобщения материалов по созданию большого числа серийных двигателей разработан интегральный метод решения этих уравнений с введением критериев подобия камеры сгорания, определяющих протекание и взаимодействие вышеописанных процессов, позволивший математически описать характеристики камеры сгорания.Получены уравнения, позволяющие определить изменение характеристик камеры сгорания в зависимости от режима работы и параметров конструкции: / -г! = (Е1нс + 0,232Е1со)10-'; Е1со=В(Р^-''' V ехр ( Т Л 288/^''' хехр (N Ф (N + с.Ф ж [1 + d^^Fj + d. ifPsfl с \-\-d. ж + d.Система уравнений и неравенств является математической моделью характеристик рабочего процесса в камере сгорания и дает возможность оптимизировать конструкцию в процессе ее автоматизированного проектирования с учетом ограничений, включающих и нестационарные режимы работы.Приведенные зависимости представлены в виде произведений функций режимных и конструктивных параметров, что дает возможность проанализировать влияние этих параметров на характеристики камеры сгорания. IX.2. Анализ влияния критериев на характеристики камер сгорания.Проявление химической кинетики в работе камеры сгорания связано с изменением температуры и давления. Изменение характеристик рабочего процесса в камере сгорания от давления и температуры приведено на рисунке 9.1 и 9.2. По оси ординат отложены отношения текущего значения характеристики к его значению при заданных условиях. С ростом Рк и Тк наблюдается снижение недожога топлива, увеличение выбросов NOx и уменьшение содержания СО и НС, потому что NOx представляет собой конечный продукт реакции (исходные - азот и кислород), а СО и НС — исходные продукты (конечные продукты на выходе из камеры сгорания - СОг и Н2О).Существенное влияние относительного количества форсунок на процессы в камере сгорания иллюстрируется на рисунке 9.3. Здесь и на последующих рис.9.4, 9.5, 9.6 по оси абсцисс отложены отношения текущего значения критерия к его значению; соответствующему минимальному недожогу топлива (максимальной полноте сгорания топлива). По оси ординат - отношение текущего значения характеристики к его значению в точке, соответствующей минимальному недожогу топлива.По физическому смыслу увеличение количества форсунок на единицу 171 I площади означает приближение состояния топливовоздушной смеси к однородной, что определяет протекание приведенных на графике характеристик.Из графика следует, что наивыгоднейщие условия запуска (максимальное значение коэффициента избытка воздз^а на запуске, соответствующее минимальному забросу температуры газа) обеспечивается вполне определенным количеством форсунок.Существует , оптимальное' количество форсунок, обеспечивающее максимальную полноту сгорания топлива (минимальный недожог) при оптимальных значениях выброса NOx и числа дымности. Снижение дымления и выброса NOx путем увеличения количества форсунок выше оптимального значения связано с увеличением недожога топлива.Критерий fiFj определяет долю воздуха, участвующего в формировании топливовоздушной смеси в отдельной горелке, а критерий W характеризует интенсивность этого процесса, рис.9.4 и 9.5. Совместное проявление этих факторов характеризует процессы смешения топливных и воздушных потоков и процессы тепломассообмена между зоной обратных токов и закрученной струей воздуха, сформированных завихрителем в отдельной горелке.Как видно из графиков, существует оптимальное значение критерия 1лР^, при котором процесс протекает с минимальным недожогом топлива, при этом запуск камеры сгорания осуществляется с минимальным забросом температуры газа перед турбиной {а максимален).Из графиков же следует, что обеспечение снижения выброса NOx и дымления приводит к резкому снижению устойчивости рабочего процесса ( а™" уменьшается) и высотности запуска Н, определяемой атмосферными условиями.Из рисунка следует, что интенсивность крутки завихрителя оказывает существенное влияние на устойчивость горения. Как и в случае //F^ имеется значение этого критерия, при котором недожог топлива и уровень дымления минимальны, а уменьшение этого критерия приводит к увеличению выброса NOx.Процессы перемешивания и взаимодействия холодных и горячих турбулентных струй, оказывающих существенное влияние на протекание химических реакций и выравнивание температур в камере сгорания, представляющей собой термогазодинамическую решетку, управляются F .Ж / важнейшим критерием у^ ' р Универсальность критерия и его влияние на характеристики камеры сгорания иллюстрируется графиками на рис.9.6. При оптимальной величине критерия имеет место не только максимальное значение полноты сгорания топлива, но и минимальная неравномерность поля температур перед турбиной.Минимальному значению недожога топлива и минимальной неравномерности температурного поля всегда соответствует вполне определенное значение потерь полного давления. Как снижение, так и увеличение потерь полного давления ухудшает основные характеристики, в том числе выбросы NOxУменьшение критерия приводит к большей устойчивости горения, но снижает высотность запуска.На рисунке показана область устойчивого горения при максимальном обеднении смеси, описанная (IX. 1) и (IX.2). По оси абсцисс отложено отношение коэффициентов скорости в отверстиях жаровой трубы к их величине при максимальном значении коэффициента избытка воздуха при срыве пламени. По оси ординат отложено' отношение коэффициента избытка воздуха при срыве пламени к его максимальному значению. Физический смысл значения коэффициента скорости в отверстиях жаровой трубы состоит в том, что его величина пропорциональна отношению 4vr для очага, / ''^ог поддерживающего горение при обеднении топливовоздушной смеси в камере сгорания.В соответствии с приведенным анализом, сужение границ по левой ветке кривой обусловлено «проскоком» пламени, а по правой ветке - «выносом» пламени из зоны стабилизации.Эмпирическое уравнение (IX.5) описывает тепловое состояние стенки жаровой трубы камеры сгорания.Влияние величины h^ на значение максимальной температуры стенки жаровой трубы получено из уравнения теплового баланса между теплом, отдаваемым стенкам жаровой трубы из ее объема, и теплом, отводимым воздухом, поступающим на охлаждение стенок жаровой трубы.Рассмотренная система уравнений явилась базой для разработки алгоритма и программ проектирования оптимальной конструкции конструкции камеры сгорания.Графическое представление уравнений и анализ физической сущности протекания процессов вооружают проектировщика знаниями о влиянии важнейших конструктивных критериев и термодинамических параметров в их взаимосвязи на изменение характеристик камеры сгорания. IX.3. Оптимальные значения критериев камеры сгорания.Оптимальные значения критериев камеры сгорания положены в концепцию ее проектирования для газотурбинных двигателей различного назначения.1 - г, NO, SN a тех.К > / / {w Hw X,, Рис.9.5. Влияние параметра W на характеристики камеры сгорания.1 -i^;2- NO ^;3- SN ;4-а'"'^.Х.1. Анализ технических требований к газотурбинным двигателям.Можно выделить следующие области использования газотурбинного двигателя: - привод для стационарной энергетической установки; - привод для транспортной машины (наземной, судовой, авиационной); - привод для высокоманевренной машины (истребительная военная авиация).' В соответствии с назначением, предъявляются различные требования к газотурбинным приводам.От стационарной энергетической установки в настоящее время требуется максимальная чистота отработавших газов, т.е. минимальная концентрация в них вредных веществ, обеспечивающая экологические требования к окружающей среде (особенно для энергетических установок, располагаемых в районах населенных пунктов) и высокая экономичность работы.От двигателя для транспортной машины требуется максимальная экономичность, устойчивость работы на переходных режимах и низкая концентрация загрязняющих атмосферу веществ в отработавших газах.Двигатель для высокоманевренной машины должен обладать наилучшими характеристиками запуска в различных высотно- климатических условиях, обеспечивающими его надежность, максимальной устойчивостью работы на переменных режимах и минимальным уровнем дымления.Х.2. Совершенная конструкция камеры сгорания в зависимости от назначения двигателя.Требования, изложенные в предыдущем разделе, определяют конструктивный облик камеры сгорания для двигателя, в зависимости от его назначения.Организация процесса горения в камере сгорания привода энергетической установки должна обеспечить минимальную концентрацию NOx, СО, НС в газах на выходе из нее.Значения критериев камеры сгорания, обеспечивающие минимальный выброс NOx, выбираются из диапазона их изменения на основании анализа отработанных серийно реализованных конструкций камер сгорания и расчетных зависимостей раздела VI: //Гз =0,6; ^ = 2 ; ^ = 1;-^10-' =4,5...5.В соответствии с материалами раздела VI, минимальный выброс NOx ' G может быть обеспечен при максимальной величине параметра — . С другой стороны, для камеры сгорания (Я^ ;^ = 0,15-г 0,3) справедливо соотнощение: V F I ' ЗГ •* Ж^3r По материалам раздела III.5, минимальная величина длины части жаровой трубы до отверстий смесителя: /з^=0,8-Ы,0Л^.Тогда: Гэг ^ ^ h ж ож I Отсюда следует, что для минимизации выбросов оксидов азота, величина параметра Д^ ж должна быть наибольшей, а величины параметров F ^ ,/Zj;j.-наименьшими. Минимально допустимая величина высоты

Основные результаты работы.

1. Теоретически обоснованы интегральные гидродинамические и геометрические критерии и параметры, определяющие процесс формирования характеристик движущегося потока газа в камере сгорания с происходящими в нем химическими реакциями окисления (горения) топлива.

2. Разработана интегральная математическая модель камеры сгорания газотурбинного двигателя.

Модель представляет систему математических зависимостей величин, характеризующих протекание рабочего процесса в камере сгорания газотурбинного двигателя от значений ее режимных и конструктивных параметров. I

Система состоит из уравнений для расчета:

- коэффициента полноты выгорания топлива,

- коэффициентов максимальной окружной и радиальной неравномерности поля температуры газа на выходе из камеры сгорания,

- величины коэффициента потерь полного давления,

- пределов воспламенения топливовоздушной смеси в камере сгорания для различных высотно-климатических условий,

- пределов устойчивого горения топлива в камере сгорания,

- степени токсичности продуктов сгорания (характеризуемой индексами эмиссии загрязняющих атмосферу веществ - оксидов азота, монооксидов углерода, несгоревших углеводородов, бензо- (а)-пирена, частиц углерода),

- максимальной величины температуры стенок жаровой трубы.

3. На базе разработанной математической модели камеры сгорания найдены оптимальные значения таких параметров ее конструкции как число топливных форсунок, относительная пропускная способность, интенсивность крутки потока воздуха на выходе из завихрителя фронтового устройства, относительная площадь отверстий жаровой трубы, а также параметров, учитывающих особенности конструкции фронтовых устройств и форсунок. В частности, оптимизация конкретной конструкции камеры сгорания авиационного двигателя по ряду вышеназванных параметров, увеличит высотность запуска двигателя более чем на 50% или снизит эмиссию оксидов азота более чем в 2 раза. Оптимальные величины геометрических параметров и критериев положены в основу концепции проектирования камеры сгорания для различных типов двигателей.

4. Обоснована концепция оптимального проектирования камеры сгорания для газотурбинного двигателя различного назначения — энергетическая установка, двигатель транспортной машины, двигатель высокоманевренной транспортной машины (военный истребитель). Суть концепции проектирования заключается в определении размеров камеры сгорания с помощью установленных значений ее критериев и параметров, обеспечивающих такие значения характеристик камеры сгорания, которые в наибольшей степени отвечают предъявляемым к двигателю требованиям.

5. Разработан алгоритм расчета размеров камеры сгорания для каждого из трех типов двигателей.

6. Выполнен анализ системы дифференциальных уравнений переноса субстанции в потоке сплошной среды с наличием источников (стоков) энергии, вещества и силы. Произведено преобразование исходной системы дифференциальных уравнений. Преобразованная форма уравнений позволяет обосновать критерии и параметры, определяющие процесс формирования характеристик потока сплошной среды с источниками (стоками) энергии и вещества в объеме рассматриваемого канала.

7. На основании преобразований уравнений Рейнольдса для установившегося турбулентного течения и уравнения Навье-Стокса для нестационарного процесI са, обоснована критериальная модель турбулентности для потока сплошной среды в канале переменной геометрии, каковым является проточная часть камеры сгорания.

8. Теоретически обоснована и определена связь температурного состояния стенок жаровой трубы с ее размерностью и формой. Получена количественная зависимость максимального значения температуры стенок жаровой трубы (с металлическим материалом стенок) от величины размера - высоты жаровой трубы в поперечном сечении.

9. На базе теоретических и экспериментальных исследований определены минимально допустимые размеры - высота поперечного сечения жаровой трубы и длина «зоны горения» жаровой трубы камеры сгорания.

10. Теоретически обоснованы и экспериментально подтверждены принципы физического моделирования рабочего процесса в камере сгорания. Принцип физического моделирования основан на равенстве величин разработанных критериев и параметров, определяющих характеристики рабочего процесса в камере сгорания для физической модели и полноразмерной камеры сгорания.

11. Теоретически обоснован безразмерный геометрический критерий, характеризующий стабильность течения газа в диффузоре камеры сгорания. Получена связь величины этого критерия со стабильностью окружной неравномерности поля температуры газа.

12. Разработана математическая модель для расчета характеристик нестационарных процессов: областей воспламенения топливовоздушной смеси в различных высотно-климатических условиях и устойчивости горения при обеднении топливовоздушной смеси в камере сгорания.

13.Материалы работы использованы при создании камер сгорания двигателя ТВа - 3000, малоразмерного турбовального вертолетного двигателя, высокотемпературного турбореактивного двигателя пятого поколения, приводов стационарной установки ГТП— 1.5, ГТЭ - 100 и при модернизации камер сгорания серийно изготавливаемых турбовальных и турбореактивных двигателей: ТВЗ-117, РД-33, вспомогательной силовой установки ГТДЭ-117.

14. Гидродинамические, геометрические критерии и параметры, полученные в I работе, по мнению' автора, могут быть применены для изучения и описания процессов, протекающих в жидкостных реактивных двигателях и других теп-ломассообменных аппаратах.

1. Абрамович Г.Н. Прикладная газовая динамика. М.: Наука. 1969. 824с.

2. Абрамович Г.Н. Теория турбулентных струй. М.: Наука. - 1986. - 716с.

3. Абрамович Г.Н. Теория центробежной форсунки. М., ЦАГИ, 1944. 114с.

4. Абрамович Г.Н. Турбулентные свободные струи жидкостей и газов. М. -Л.: Госэнергоиздат. - 1948. - 288с. '

5. Авиационные ГТД в наземных установках. С.П.Изотов, В.В.Шашкин, В.М.Капралов и др. Под общ. ред. В.В.Шашкина. Л. Машиностроение. Ленинград, отделение. 1984. 288с. (Экономия топлива и электроэнергии).

6. Адамов В.А. Сжигание мазута в топках котлов. Л.: Недра. 1989. - 304с.

7. Акатнов Н.И., Лавров А.В. О влиянии пульсаций температуры на образование NOx в высокотемпературной среде // Тр. ЛПИ. Сер. Механика иIмашиностроение. — Л. — 1976. — №352. с.46 51.

8. Ахмедов Р.Б. Дутьевые горелочные устройства. М. Недра. 1977. 272с.

9. Ахназарова С.Л., Кафаров В.В. Оптимизация эксперимента в химии и химической технологии. М.: Высш. школа. 1978. 319с.

10. Аэродинамика закрученной струи. Под ред. Р.Б.Ахмедова. М. Энергия. 1977. 240с.

11. Бронштейн И.Н., Семендяев К.А. Справочник по математике для инженеровIи учащихся вузов. Москва. Наука. 1980.

12. Влияние ввода влаги в топку БКЗ 320 - 140 ГМ на выбросы окислов азота/

13. А.Д.Горбаненко, С.П.Титов, В.П.Лукашявичус и др.// Электрические станции. 1984. №5. С.59 60.

14. Тумановский А.Г., Тульский В.Ф. /Влияние впрыска воды на образование окислов азота за камерой сгорания с последовательным вводом воздуха в зону горения // Теплоэнергетика. 1982. №6. С.34 36.

15. Внуков А.К., Алыневский В.Н. Влияние двухступенчатого сжигания и рециркуляции газов на генерацию оксидов азота и серного ангидрида. -Электрические станции. 1980. №7. С. 15 18.

16. Котлер В.Р. /Выбросы оксидов азота при совместном сжигании угля с газом или мазутом // Теплоэнергетика. 1996. №5. — С.47 - 57.

17. Гаврилов А.Ф., Горбаненко А.Д., Туркестанова Е.А. /Влияние влаги, вводимой в горячий воздух, на содержание окислов азота в продуктах сгорания газа и мазута// Теплоэнергетика. 1983. №9. С. 13 15.

18. Методика нормирования эмиссии вредных веществ турбореактивными двигателями дозвуковых самолетов гражданской авиации в зоне аэропорта/ Горбатко А.А., Щербаков В.И., Худяков Е.И. и др. // Труды ЦИАМ. 1982.

19. Горбунов Г.М. Выбор параметров и расчет основных камер ГТД. Учебное пособие. М. Изд. МАИ им. Орджоникидзе. 1972. 229с.

20. Горелочные устройства с пониженным выходом NOx // Механика и энергетика. Выпуск 5. Москва. 1991.

21. ГОСТ 17.22.04.86. Двигатели газотурбинные самолетов гражданской авиации. Нормы и методы определения выбросов загрязняющих веществ.

22. Груздев В.Н. Концентрация токсичных компонентов в продуктах сгорания керосино-воздушных смесей // Изв. ВУЗ. Сер. Авиационная техника. 1978. — №1. — С.49 - 52.

23. Длин A.M. Факторный анализ в производстве. М. «Статистика». 1975.

24. Дорошенко В.Е. О процессе горения в камере ГТД // Тр. ЦИАМ. 1959. №354.

25. Дрегалин А.Ф., Черенков А.С. Общие методы теории высокотемпературных процессов в тепловых двигателях. — М. «Янус — К». 1997. 328с.

26. Зависимость эмиссии NOx от конструктивных и режимных параметров камеры сгорания газотурбинного двигателя /А.А.Саркисов, О.А.Рудаков, Н.Д.Саливон, Ю.В.Сигалов, В.А.Митрофанов // Теплоэнергетика. 1999. №12.I

27. Зельдович Я.Б., Садовников П.Я., Франк-Каменецкий Д.А. Окисление азота при горении. М.: Изд-во АН СССР. 1947.

28. Зуев B.C., Скубачевский JI.C. Камеры сгорания воздушно-реактивных двигателей. М. Изд. МАИ им. Орджоникидзе. 1958.

29. Идельчик И.Е. Справочник по гидравлическим сопротивлениям. М.: Машиностроение. 1975.-559с.

30. Тумановский А.Г. /Испытание кольцевой камеры сгорания авиационного двигателя НК-8 на природном газе // Теплоэнергетика. 1976. №8. - С.60 -64.

31. Исследование области воспламенения топлива в камере сгорания /А.А.Саркисов, О.А.Рудаков, Н.Д.Саливон, В.А.Митрофанов // Теплоэнергетика. №3. 2003.

32. Камеры сгорания авиационных газотурбинных двигателей./ Под редакцией Б.П.Лебедева, Г.Н.Абрамовича и Ю.Ф.Дитякина // Министерство авиационной промышленности СССР. Институт им. П.И.Баранова.

33. Канило П.М. Токсичность ГТД и перспективы применения водорода. Киев: Научная мысль. 1982.

34. Канило П.М., Подгорный А.Н., Христич В.А. Энергетические иэкологические характеристики ГТД при использовании углеводородныхtтоплив и водорода. Киев: Наукова думка. 1987. -224с.

35. Котлер В.Р. Оксиды азота в дымовых газах котлов. М. Энергоатомиздат. 1987.- 144с.

36. Кропп Л.И. Пути сокращения вредных выбросов ТЭС // Теплоэнергетика. -1978. №11. - с.2 — 7.

37. Круглов М.Г., Меднов А.А. Газовая динамика комбинированных двигателей внутреннего сгорания. Учебное пособие. — М.: Машиностроение. 1988. -360с.

38. Кузнецов В.Р. Аналитические методы определения концентрации окиси азота в камерах сгорания газотурбинных двигателей. Новое в зарубежном авиадвигателестроении. 1973. №2.

39. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика: Учебное пособие. TVI.

40. Гидродинамика. М.: Наука. 1988. - 736с.

41. Ластовцев A.M. Гидродинамический расчет вращающихся распылителей. Труды МИХМ т. 11. Москва. 1957.

42. Лефевр А. Процессы в камерах сгорания ГТД. М.: Мир. 1986. 566с.I

43. Лукачев В.П., Белоусов А.Н., Ланский A.M. Влияние конструктивных параметров вихревых газовых горелок на характеристики кольцевой авиационной камеры сгорания. // Горение в потоке. Казань. КАИ. 1982.

44. Малотоксичные камеры сгорания для энергетических ГТУ/ Тумановский А.Г., Гутник М.Н., Соколов К.Ю. // Теплоэнергетика. 1997. №3. С.48 52.

45. Математическая модель процессов образования и расчета загрязняющих веществ и оптимизация камер сгорания ГТД/А.А.Саркисов, О.А.Рудаков, Н.Д.Саливон, Ю.В.Сигалов, В.А.Митрофанов // Теплоэнергетика. 2000. №5. С.52 55.

46. Постников A.M. /Методика обработки результатов измерений выбросов вредных веществ газотурбинными двигателями // Теплоэнергетика. 1992. №9.С.64 66.

47. Камера сгорания ГТД с низким уровнем эмиссии/ Митрофанов В.А., Федоров A.M., Рудаков О.А' ., Сигалов Ю.В.// Тезисы докладов XLIII Научно-технической сессии по проблемам газовых турбин. 25-26 июня Санкт- Петербург 1996года.

48. Рабочий процесс и расчет камер сгорания газотурбинных двигателей /Михайлов А.И., Горбунов Г.М., Борисов В.В., Квасников Л.А., Марков Н.И.

49. Тр. МАИ им. Орджоникидзе. 1959. Вып. 106. 285с.

50. Некоторые вопросы математического описания и оптимизацииIмногофакторных процессов. Г.К.Круг и др. Труды МЭИ. 1963. Вып.51.

51. Тумановский А.Г./ Некоторые особенности образования окислов азота в высокофорсированных камерах сгорания с последовательным вводом воздуха в зону горения/ // Теплоэнергетика 1977. - №12. С.70 - 72.

52. Новое в зарубежном авиадвигателестроении. №1. 1978.

53. Новые газотурбинные установки фирмы Семенз /КВУ./ Новостизарубежной науки и техники. Серия «Двигатели для авиации иtкосмонавтики» №5 6. 1992. С. 15 - 20.

54. Новый газотурбинный двигатель мощностью 110 МВт для стационарных энергетических установок/ Романов В.И., Рудометов С.В., Жирицкий О.Г., Романов В.В. // Теплоэнергетика. 1992. №9. С.27 — 31.

55. Обобщенная эмиссионная характеристика ГТД как функция конструктивных и режимных параметров камеры сгорания /А.А.Саркисов, О.А.Рудаков, Н.Д.Саливон,' Ю.В.Сигалов, В.А.Митрофанов // Теплоэнергетика. 2000. №3.

56. Образование и разложение загрязняющих веществ в пламени: Пер. с англ. //

57. Ред. Н.А.Чигир. М.: Машиностроение. 1981. — 407с.

58. Основы практической теории горенйя: Учебное пособие для вузов/ В.В.Померанцев, К.М.Арефьев, Д.Б.Ахмедов и др. Л.: Энергия. 1973. 264с.

59. Отчет о научно-исследовательской работе «Исследование экологических характеристик двигателей ТВ2-117, ТВЗ-117 и ГТД-350». Центр социальной экспертизы республиканского инновационного фонда СМ РСФСР. 1991.

60. Отчет экспертизы по выбросам изделиями 88, 21, 55Б(77) веществ, загрязняющих окружающую среду. ЭД85.90. ТМКБ «Союз». 1990.

61. Павлов В.А., Штейнер И.Н. Условия оптимизации процессов сжиганияжидкого топлива и газа в энергетических и промышленных установках.: Л.

62. Энергоатомиздат. 1984. 120с.

63. Перемешивание поперечной натекающей на стенку струи с основным потоком./ Новости зарубежной науки и техники. Серия «Двигатели дляавиации и космонавтики» 1992. №5 6. С.32 — 36.

64. Перспективные проекты газовых турбин в Японии./ Новости зарубежной науки и техники. Серия «Двигатели для авиации и космонавтики» 1992. №5 -6. С.32-36. ' |

65. Применение концепции двухзонной камеры сгорания с низким уровнем эмиссии NOx-/ Новости зарубежной науки и техники. Серия «Двигатели для авиации и космонавтики» 1992. №3. С.34 39.

66. Пчелкин Ю.М. Камеры сгорания газотурбинных двигателей. Учебноепособие.- 3-е изд. перераб. и доп. М. Машиностроение. 1984. 280с.

67. Разработка критериев запуска камер сгорания ГТД./ Саркисов А.А., Рудаков1

68. О.А., Саливон Н.Д., Сигалов Ю.В., Митрофанов В.А.// Вестник СГАУ. Сер. Процессы горения, теплообмена и экология тепловых двигателей. Самар. Гос. Аэрокосм. Университет. Самара. 1999. Вып.2. -С.177-183

69. Разработка методов снижения эмиссии загрязняющих веществ стационарными газотурбинными установками на природном газе. Научно-технический отчет 95-001 ЦИАМ. «Эколэн». 1995.

70. Распыливание жидкостей / В.А.Бородин, Ю.Ф.Дитякин, Л.А.Клячко, В.И.Ягодкин - М. 1967.

71. Распыливание жидкостей / Ю.Ф.Дитякин, Л.А.Клячко, Б.В.Новиков, В.И.Ягодкин. М., Машиностроение. 1977. 208с.

72. Распыливание топлива центробежными форсунками турбореактивных двигателей / В.А.Бородин, Ю.Ф.Дитякин, Л.А.Клячко, Н.Н. Струлевич М.1. ЦИАМ. 1952.

73. Расчет устойчивости горения в ГТД / Саркисов А.А., Рудаков О.А., Саливон Н.Д., Сигалов Ю.В., Митрофанов В.А.// Вестник СГАУ. Сер.Процессы горения, теплообмена и экология тепловых двигателей. Самар. Гос. Аэрокосм. Ун-т. Самара. 2000. Вып.З.-С.262-266.

74. Рудаков О.А. Методика расчета оптимальных параметров камеры сгорания высокотемпературного ГТД. Научно-технический сборник. Серия XX. 1978. Выпуск 80. С. 8- 11.

75. Рудаков О.А., Пыховский JI.JI. Расчетная оптимизация камеры сгорания по удельному расходу топлива ГТД /Автоматизация проектирования авиационных двигателей// Тезисы докладов IV отраслевой конференции. Труды ЦИАМ №1095. 1984. С.81 82.

76. Руководящий технический материал авиационной техники. Р.Т.М. 1626-80. Камеры сгорания газотурбинных двигателей. Метод поверочного гидравлического расчета на ЭВМ. 1980. 39с.

77. Салливан Е. Простое уравнение для расчета выбросов NOx из камеры сгорания газотурбинного двигателя, учитывающее загрязнение воздуха./ Труды американского общества инженеров-механиков. Сер. Энергетические машины и установки. 1978. №2. С.1 8.I

78. Сжигание жидкого топлива в кольцевой камере сгорания конструкции H3JI/ А.Г.Тумановский, А.В.Сударев, В.В.Маев и др. // Теплоэнергетика. 1986. №3. С.37 —42.

79. Система уравнений для расчета характеристик камеры сгорания. / Саркисов А.А., Рудаков О.А., Саливон Н.Д., Митрофанов В.А.// Промышленная теплотехника. 2001.-Т.23.- №3.- С.83-88.

80. А.Г.Тумановский, М.Н.Гутник, К.Ю.Соколов /Снижение выбросов оксида углерода и углеводородов на пусковых и переменных режимах энергетических ГТУ при сжигании жидких и газообразных топлив // Электрические станции. 1991. №11. С.2 6.

81. Совершенствование ГТУ для электростанций. / Г.Г.Ольховский, А.Г.Тумановский, А.И.Механиков и др. // Теплоэнергетика. 1991. №6. С.66 -73.

82. Современное состояние развития газотурбинных установок в Европе./ Новости зарубежной науки и техники. Серия «Двигатели для авиации и космонавтики» 1992. №5 6.С.37 - 40.

83. Создание перспективных газовых турбин в США./ Новости зарубежной науки и техники. Серия «Двигатели для авиации и космонавтики» 1992. №5 -6. С.7-11.

84. Соколов К.Ю. Совершенствование камер сгорания ГТУ на основеIматематического моделирования рабочего процесса. // Авторефератдиссертации на соискание ученой степени кандидата технических наук. -Москва: ВТИ. 1989. 23с.

85. Спейшер В.А., Горбаненко А.Д. Повышение эффективности использованияIгаза и мазута в энергетических установках. М. Энергоиздат. 1982. 239с.

86. Сторожук Я.П. Камеры сгорания стационарных газотурбинных и парогазовых установок. — JL: Машиностроение. 1978. 231с.

87. Сторожук Я.П., Ожигов Г.Е. Обзор материалов по газомазутным горелочным устройствам зарубежных фирм. JI. Тр. ЦКТИ. 1975. Вып. 128.

88. Сторожук Я.П., Павлов В.А. Процессы распыливания и смесеобразования при сжигании мазута в топках паровых котлов. — В кн.: Опыт сжигания мазута и газа на электростанциях. М. Энергия. 1968. С.46 71.

89. Сударев А.В. Перспективы создания для стационарной энергетики экологических газотурбинных двигателей // Изв. РАН. Энергетика 1992. №1.

90. Сударев А.В., Маев В.А. Камеры сгорания газотурбинных установок. Интенсификация горения. JL: Недра. 1990. - 274с.

91. Талантов А.В. Горение в потоке. М.: Машиностроение. 1978. 160с.

92. Теория двухконтурных турбореактивных двигателей/ Под ред. д-ра техн. наук проф. С.М.Шляхтенко, д-ра техн. наук проф. В.А.Сосунова. М. Машиностроение. 1979. 432с.

93. Теория и расчет воздушно-реактивных двигателей./ Под ред. С.М.Шляхтенко. Учебное пособие. 2-е изд. перераб. и доп. - М. Машиностроение. 1987. 568с.

94. Теория и расчет образования вредных выбросов в камерах сгорания ГТД/ А.А.Саркисов, О.А.Рудаков, Н.Д.Саливон, Ю.В.Сигалов, В.А.Митрофанов. // Промышленная теплотехника. 1999. Т.21.- №6. - С.34-38.

95. Тишин А.П., Худяков В.А., Артамонов А.К. Исследование возможностей уменьшения концентрации окислов азота при сжигании топлив в теплоэнергоагрегатах./ Калининград. М.О. ЦНИИмаш. 1994. 60с.

96. Тумановский А.Г. Исследование камер сгорания ГТУ и разработка мероприятий по снижению концентраций окислов азота // Автореферат диссертации на соискание ученой степени кандидата технических наук. — Москва: ВТИ. 1971.

97. Тумановский А.Г. Образование окислов азота в камерах сгорания ГТУ при сжигании природного газа // Теория и практика сжигания газа. — Вып.5. Л.: Недра. 1972. - С.330 - 340. t

98. Тумановский А.Г. Предотвращение загрязнения атмосферы ГТУ // Энергетика и охрана окружающей среды. М.: Энергия. 1975. С.223 -234.

99. Турбулентное смешение газовых струй. Под редакцией Г.Н.Абрамовича. Изд-во «Наука». Главная редакция физико-математической литературы. 1974. 272с.

100. Устименко Б.П., Ткацкая О.С. Аэродинамика закрученной струи. В кн.:<

101. Проблемы теплоэнергетики и прикладной теплофизики. Алма-Ата. «Наука». 1970. Вып.6. С.211-216.

102. Физические основы рабочего процесса в камерах сгорания ВРД/

103. Б.В.Раушенбах, С.А.Белый, И.В.Беспалов и др. М.: Машиностроение. 1964.525с.

104. Франк-Каменецкий Д.А. Диффузия и теплопередача в химической кинетике. Изд-во АН СССР. 1947.

105. Христич В.А., Тумановский А.Г. Газотурбинные двигатели и защита окружающей среды. Киев: Техника. 1983. 144с.

106. Шаулов Ю.Х., Лернер М.О. Горение в жидкостных ракетных двигателях./ Государственное научно-техническое издательство ОБОРОНГИЗ. Москва. 1961.

107. Щуровский В.А., Синицын Ю.Н.' Экологические характеристикигазотурбинных агрегатов на переменных режимах. // Газовая промышленность. 1991. №11.

108. Экспериментальное исследование фронтового устройства камерыIсгорания ГТЭ-115 с предварительным образованием гомогенной топливо-воздушной смеси./ А.Г.Тумановский, К.Ю.Соколов, М.Н.Гутник и др. // Теплоэнергетика. 1993. №7. С.42 48.

109. Яворский Б.М., Детлаф А.А. Справочник по физике. М., 1977., 944с.

110. Зысина-Моложен Л.М., Зысин Л.В., Поляк М.В. Теплообмен в турбомашинах.«Машиностроение» ( Ленинград, отд-ние). 1974. 336с.

111. Технический отчет ЦИАМ № 5836. Исследование течения в моделяхдиффузоров камер сгорания. Авторы: Безменов В.Я., Асеев Г.Д., Тарараксина А.П. 1967.

112. Лойцянский Л.Г. Механика жидкости и газа. М.: Наука. 1970. 847с.

113. Верцинский 3., Супрун Т.Т., Эпик Э.Я. Характеристики ламинарно-турбулентного перехода, индуцированного следом одиночного движущегося цилиндра // Промышленная теплотехника. 2001. Т. 23.№ 3. С. 22 30.

114. Appleton J.P., Heywood J.B.< The effects of Imperfect Fuel-air Mixing in a Burner on NO Formation from Nitrogen in the Air and Fuel. 14-th Symposium (international) on Combustion.

115. Bahr D.W., Control and Reduction of Aircraft Turbine Engine. Exhaust Emissions, in W.Cornelius and W.G.Agnew (eds.), Emissions from Continous Combustion Systems. Plenum. New York. 1972. P. 345 372

116. Bowmann C.T. Kinetics off pollutant formation and distinction in combustion. Progress in Energy and Combustion Sc. 1975. Vol.1. №1.

117. Brandon D.B. I.S.A. Journal. 1959. 6. №7.

118. Buchheim R. Influences on Exhaust Emissions from Automotive Gas Turbines, ASME. 1978. Paper 78 GT - 85i

119. Gleason C.C., Rogers D.W., Bahr D.W. Experimental clean combustor program, phase 2. NACA. CR- 13 1971, 1976.

120. Grobman J. Effect of Operating Variables on Pollutant Emissions from Aircraft Turbine Combustors, in W.Cornelius and W.G.Agnew (eds.), Emissions from Continous Combustion Systems. Plenum. New York. 1972. P. 279 303

121. Hey wood J.B. Gas Turbine Combustor Modelling for Calculating Nitric Oxide Emissions, AIAA. Paper 1971. №712.i

122. Lavoie G.A., Heywood J.B., Keck J.C. Experimental and Theoretical Study of Nitric Oxide Formation in Internal Combustion Engines, Combust. Sci. Technol. 1970. Vol.l.P. 313-326

123. Lipfert F.W. Correlation of Gas Turbine Emission data. ASME Publication. 1972. Paper 72-GT-60.

124. Morton H.L., Marshall R.L. Impact of emission regulations on gas turbine combustion design. Canad. Aeronaut. Space, Journal. 1974. Vol.20, №3.

125. Quan F., Marble F.E., Klilgel J.R. Nitric Oxide Formation in Turbulent Diffusion Flames. 14-th Symposium (international) on Combustion.

126. Roberts R., Fiorentino A.J., Diehl L.A. The pollution reduction technologyprogramm for can-annular combustor engines. AIAA. 1976. Paper№76 761.

127. Roberts R., Peduzzi A. Low pollution combustor for GTOL engines. AIAA. 1976. Paper №76-761.

128. Verkamp E.J., Verdouw A.J., Tomlinson J.G. Impact of emission regulations on Future Gas Turbine Engine Combustors. AIAA. 1973. Paper №73 1277.

129. Westenberg A.A. Combust. Sci. Technology. 1971. 59.Vol. 4.1. Патенты.

130. Патент на изобретение № 2094705 /Камера сгорания авиационного или энергетического газотурбинного двигателя. Сигалов Ю.В., Рудаков О.А., Ефимов Е.В., Гурский С.Э., Митрофанов В.А., Федоров A.M.

131. Патент на изобретение № 2145669 /Способ работы газотурбинного двигателя и огневой блок для сжигания смеси топлив. Саркисов А.А., Митрофанов В.А., Рудаков О.А., Саливон Н.Д., Сигалов Ю.В.I

132. Патент на изобретение № 2171432/Камера сгорания с оптимальным числом форсунок. Саркисов А.А., Митрофанов В.А., Рудаков О.А., Саливон Н.Д., Сигалов Ю.В.