Методика расчёта основных параметров камеры пульсирующего горения Текст научной статьи по специальности «Математика»

УДК 621.438

МЕТОДИКА РАСЧЁТА ОСНОВНЫХ ПАРАМЕТРОВ КАМЕРЫ ПУЛЬСИРУЮЩЕГО ГОРЕНИЯ

© 2012 А. В. Солодовников, Е. Н. Вышегородцев, В. В. Голубятник

Серпуховской военный институт ракетных войск

Рассматривается методика расчёта параметров прямоточной камеры пульсирующего горения (V = const). Приведён алгоритм расчёта.

Камера пульсирующего горения, математическая модель, расчёт, испытания.

Принципиально новым направлением в двигателестроении является разработка силовых установок, работающих на принципе высокоскоростного сгорания топлива ({ -const) и использующих камеру пульсирующего горения (КПГ) с резонансной трубой (РТ), например, пульсирующие воздушно-реактивный [1] и жидкостной ракетный двигатели [2] или турбореактивный двигатель с вибрационной камерой. При этом КПГ потребляет меньше топлива, чем традиционная, при одинаковой их мощности, и в связи с этим фактом ожидаемая эффективность с точки зрения экономии горючего составит 10...30% [1].

Проведённая авторами статьи научно-исследовательская работа [2, 3] за период с 2006 по 2012 годы в Серпуховском военном институте ракетных войск (СВИ РВ) позволила заключить, что одним из возможных вариантов КПГ (рис. 1) является камера сгорания (КС), оборудованная механическим мембранным клапаном, расположенным перпендикулярно к её оси, и резонансной трубой, которая обеспечивает устойчивый пульсирующий режим горения топливной смеси (V = const). Принцип работы камеры пульсирующего горения следующий. Воздух всасывается через мембрану вследствие возникающего в КС перепада давления в конце фазы выхлопа. При каждом заполнении камеры воздухом в неё подаётся горючее для образования топливной смеси заданного состава. Сжатие рабочего тела осуществляется вследствие сгорания топлива. Расположение впускного клапана перпендикулярно к оси камеры сгорания исключает контакт горячего газа с мембраной и, соответственно, повышается надёжность и ресурс всего агрегата [1].

Рис. 1. Модельная КПГ с перпендикулярным расположен нем м ехан ического мембранного клапана

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

В настоящее время разработано несколько различных методов инженерного расчёта КПГ [1], применяемых для проектирования аэрозольных генераторов и пульсирующих воздушно-реактивных двигателей, которые основаны на применении относительных величин для расчёта геометрических размеров агрегатов. Эти данные носят экспериментальный характер и, соответственно, успех в решении задачи определяется научной и технической интуицией конструкторов.

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

дельных элементов КПГ (например, резонансной трубы, форкамер).

Проведённый анализ математических моделей показал, что на данный момент не существует единого подхода к описанию рабочего процесса в КПГ. Так, в работе [4] выявлено пять наиболее распространённых подходов к моделированию рабочего процесса, протекающего в камерах пульсирующего горения: акустический, механический, газодинамический, волновой и вихревой.

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

По мнению авторов, для комплексного описания цикла пульсирующей КС необходимо применить метод расчёта, основанный на волновой теории моделирования рабочего процесса. Данный метод основан на составлении и решении систем уравнений, выражающих основные законы сохранения, которые могут быть записаны как в интегральной, так и в дифференциальной формах [4]. Поскольку описание даже одномерного нестационарного движения газа основывается на дифференциальных уравнениях в частных производных гиперболического типа, не имеющих аналитического решения, то эти системы уравнений, как правило, решаются различными численными методами.

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

надёжность получения сходящегося решения.

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

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

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

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

В-четвёртых, в теории одномерных течений имеется достаточное количество решённых конкретных задач.

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

Ьк І'С

а 1 ► «4 ►

-

і V -ф л Ь V ►

Л

СІфі.

Рис. 2. Расчётная схема КПГ

Эти уравнения имеют следующий вид:

— (р-р) + — (р-и-р) = О,

Ы дх

д

а

д_

а

(р-г)

дс

( и^

ел—

V 2У

др-Р)

а

-т0-/:

(1)

_э

йг

р-и-Р ■

и2

= #•/

У

где /> - плотность газа, кг/м ;

/' - площадь сечения канала, м2; и - скорость движения газа, м/с; е - удельная внутренняя энергия;

Р - давление, Па;

То - среднее по периметру напряжение силы трения;

/ - периметр сечения канала, м; q - количество теплоты, приходящее в единицу времени через единичную площадь поверхности, Дж/м2;

^ - время, с;

х - координата вдоль оси канала, м. Первое уравнение системы (1) - уравнение неразрывности, отражающее закон сохранения массы; второе - уравнение движения, определяющее закон сохранения количества движения; третье - уравнение энергии, представляющее закон сохранения энергии. Система (1) справедлива при следующих допущениях:

- движение рабочего тела одномерное;

- вектор скорости потока в каждой точке поперечного сечения является нормальным по отношению к плоскости сечения;

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

- влиянием массовых сил пренебрегаем;

- стенки канала непроницаемы для рабочего тела;

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

Замкнутое описание математической модели процессов в рабочем теле (газе) предполагает дополнение системы (1) термическим и калорическим уравнениями состояния, которые имеют различный вид в зависимости от модели рабочего тела (иде-

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

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

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

В качестве упрощений вводятся следующие допущения:

- не учитывается наличие трения о стенки канала;

- не учитывается наличие дополнительных источников массы.

С учётом вышесказанного, а также допущений для системы (1), уравнения нестационарного одномерного движения невязкого сжимаемого газа в проточной части КПГ могут быть представлены в виде:

др др д1<’ др ди л Р— +р— + ри— + иЕ— + рр— = О ск ск дх дх дх

ди ди дР р — + ри — + — = О,

с1 дх дх

>

др дР 2

— + и — - а 9/ дх

др др

— + и — \ 9/ дх

= (к - 1 )/з - и ■ q■

(2)

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

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

Для решения системы уравнений (2) целесообразно использовать метод характеристик, алгоритмическая и программная реализация которого была отработана в Серпуховском военном институте ракетных войск при проектировании различных КПГ для ракетных двигателей, энергетических и аэрозольных генераторов.

Методика инженерного расчёта КПГ, предложенная авторами статьи, представлена на рис. 3.

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

Представленная математическая модель позволяет определить искомую длину

резонансной трубы (Ьрт\ основываясь на газодинамическом течении пульсирующего газового потока. Предлагаемый способ определения Ьрт наиболее приемлем для инженерных расчётов, так как моделирование процесса распределения давления в РТ учитывает такие параметры, как массы топлива и воздуха, необходимые для одного цикла работы КПГ. Всё это, в свою очередь, повышает точность расчётов.

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

В некоторых задачах, например расчёт КПГ для аэрозольного генератора, необходимо определение теплового потока по длине камеры сгорания и резонансной трубы. Поэтому для решения приведённой системы уравнений (2) необходимо знать функцию удельного теплового потока С{ =/(х, I), подводимого к газу и обусловленного подводом тепла при горении топливной смеси в КПГ, а также теплообменом при движении рабочего тела в проточной части. Процессы смесеобразования, воспламенения и сгорания топливной смеси в КС являются важнейшими и одновременно наименее изученными явлениями. Динамика рабочего процесса в КС определяет работоспособность проектируемого агрегата, экологичность и эффективность использования топлива.

В основу рабочего процесса камеры пульсирующего горения авторами статьи заложена динамика резонансно-волновых процессов. При этом источником энергии для всего рабочего цикла является давление, создаваемое в КПГ. Амплитуда и длина волны определяется процессами, протекающими в КС. Так как рабочий процесс КПГ имеет циклический характер, то к времени протекания процессов смесеобразования, воспламенения и сгорания предъявляются жёсткие требования, определяемые геометрией проточной части энергоустановок (двигателя).

Исходные данные: Топливо, тяга (мощность), частота пульсаций, а, температура горения.

Вычисляется: масса продуктов сгорания за один период, объём зоны горения топлива У1г длина Ьз и диаметр зоны горения (Із

ли

ь3

Вычисляется: площадь форкамеры диаметр ^ и

длина ф оркамеры Z^м. диаметр объём Уж и длина камеры сгорания 1ке

_____________XX

Вычисляется: длина

сужающегося участка диаметр резонансной трубы ^

XX

Математическим моделированием движения газового потока по резонансной трубе, опр ед еляется длина волны продуктов сгорания и, с о ответственно, длина резонансной трубы Ьрк

______________XX

Итоговые параметры КПГ:

1. Геометрические размеры.

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

3. Распределение тепловых потоков по длине КПГ и резонансной трубе

Х^

£

Рис. 3. Методика расчёта камеры пульсирующего горения (пример определения оптимальной длины резонансной трубы)

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

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

Кроме того, применение сложной модели динамики горения приводит к существенному усложнению математической модели рабочего процесса КПГ и делает её «тяжёлой» даже при использовании возможностей современной вычислительной техники. В ряде работ, например [1], рабочий процесс камеры сгорания рассматривается в упрощённой форме.

Для определения удельного теплового потока с! =/(х, I), необходимого для решения (2), примем следующие допущения:

- теплоподвод осуществляется только в КС, а по длине РТ не учитывается;

- зона теплоподвода в КС совпадает с зоной камеры и не изменяет своего положения до полного выгорания смеси;

- распределение подводимой теплоты по длине зоны теплоподвода равномерное;

- смесь считается однородной;

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

Для описания динамики процесса теплоподвода в зоне тепловыделения используются различные модели динамики сгорания смеси. Например, в работе [4] предложено использовать уравнение К. Неймана, которое было получено исходя из представления реакции горения как бимолекулярной реакции между горючим и кислородом воздуха:

0,21 -Ку.М0.Оц

П

(!-,)•

С \

а

--1

-У’

(3)

где V-

Ц'

у - доля прореагировавшего горючего; Ку - условная константа скорости реак-

Сц - подача топлива за цикл, кг/с;

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

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

т - текущее время на интервале одного цикла, с.

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

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

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

%{?) =-----------------, (4)

Бк-т

где у(т) - доля топлива, прореагировавшего к моменту времени т, получаемая из решения уравнения (3) на каждом временном шаге при численном интегрировании;

Сц - масса топлива, поступившего в КС за один цикл пульсации, кг/с;

(2т ~ теплотворная способность компонентов топлива;

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

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

Следует заметить, что уравнение (3) может быть решено аналитически (в квадратурах). Представим уравнение (3) в виде

Г \

ции;

Ф

ёт

а

--1

'У:

(5)

Ук - объём КС, м ;

где А =

0,21 -к-м -а, у О Ч

VI

Уравнение (5) преобразуем к виду

Ф 2 ( \

— -А-у +А-(а + 1)- у-А-а = 0 (6)

ск

или

Ф 2

— л-а-у + Ь- у + С = 0, (7)

ск

где а = - А;

Ь = А-(а + 1); (8)

с = -Аа.

Уравнение (7) относится к классу уравнений Риккати [4]. Как известно, если коэффициенты а, Ъ, с в уравнении Риккати постоянны, то допускается разделение переменных и общий интеграл уравнения (7) имеет вид:

ф

С\-т = \---^--------- (9)

а-у + Ъ-у + С

Выражение для интеграла от алгебраической функции в (9) равно [4]

Ф 1

I-

■у^+Ь-у+С Ь2-4-а-с

а

2-а-у+Ъ-^Ъ -4-а-с

х1п ----------,

2-а-у+Ь+^Ь -4-а-с если Ь2 - 4ас > 0,

Ф

(10)

I- „

а-у +Ь-у + С V 4 • а • с-1? 2 ■ а ■ х + Ъ

х ап^ . ,

л/ 4 - а ■ с-Ъ если Ь2 - 4ас < 0,

Ф

I-

а-у + Ъ-у + С 2-а-х + Ъ

(П)

(12)

если Ъ - 4ас = 0.

(8):

Определим значение Ъ - 4ас с учётом

Ь2 -4-а-с=А2 -(а+1)2 -4-А2-а=А2-(а-1)2. Поскольку а - 7, то (Ъ2 - 4ас) > 0 и выражение для интеграла (9) имеет вид (10).

Обозначим

(13)

g = л/&2 - 4- а- с = А-(а-\У

Тогда решение (9) с учётом (10) запишется:

1 2■a■y + b-g

g 2■a■y + b + g

(14)

откуда

2 • а■у + Ъ-g 1п---------------= - тg

или

2■a■y + Ъ + g Я

У =

а (х ~С 2 ■ е )

2 а

(15)

Определим постоянную С2 из началь ных условий решения уравнения (3): т = 0, то у = 0.

Тогда из (15) имеем

Ъ + £

Я

а

откуда

О -с2)

= 0

2 а

(16)

(17)

(18)

С. =1--^-,

2 b + g

или с учётом (8)

1

с2 =— .

а

Подставив выражения (18) и (8) в (15), получим окончательное решение уравнения (3) в виде

/ \ а-1

у(т)=а-—------------------ . (19)

1 -А(а-1)т \- — е у '

а

Отметим, что при т —> 0, у —> 0, а при т —> ас, у —> 1, что соответствует физической сущности процесса сгорания топлива и косвенно подтверждает корректность решения (19).

С учётом (19) выражение (4) для удельного теплового потока в КС на временном интервале цикла пульсаций примет вид:

{ \

а-1

а -

л 1 -А(а- 1)т

1- — е у ’

,(20)

где В =

Сч<2тг!т

(21)

Выражение (20) определяет характер изменения удельного теплового потока на временном интервале цикла пульсаций и повторяется периодически с периодом, равным периоду пульсаций Г, т. е. 0 < т < Т. Для интегрирования уравнения (2) необходимо определить характер изменения удельного теплового потока во времени q(t) для всего времени работы КПГ, т. е. при 1>0.

Таким образом, функция q(t) является негармонической периодической функцией времени с периодом Т и удовлетворяет условиям Дирихле, т. е. она ограничена, имеет разрывы первого рода с конечным числом экстремумов на конечном интервале. В этом случае изменение q(t) в течение всего времени работы КПГ можно представить рядом Фурье

?(0= “Г +

00

+ £ аг

П—\

> (22)

где

сое п ■ со ■ Ї + Ъп ■ БІП п ■ СО ■ Ї

а0 =—(23)

то

2 т

ап = — ) q ^ - соэ п- со- тЛ, (24)

ТО1!

2 т

ъп = — 1 % V ) ■ 8ІП- СО ■ т Л, (25)

ТО

со =

2ж

Т

- круговая частота пульсации.

На основе (20), (23-25) получим выражение для коэффициентов Фурье. В форме (22) можно математически описать изменение q в течение всего времени работы КПГ.

Однако для интегрирования уравнений (2) необходимо задать изменение удельного теплового потока не только во времени, но и по длине КПГ. В предположении, что теп-лоподвод осуществляется только в КС, выражение удельного теплового потока как функции координаты х и времени t имеет вид:

Ґ

q(x,t) =

о

+2Х

V п=1

л

■сс&п-со-і+Ь -^тп-со-і

о=0 при х >ЬКС,

Ьк - длина камеры сгорания. Приведённая методика расчёта удельного теплового потока совместно с уравнениями основных параметров КПГ использует вычислительный эксперимент, в результате которого получено распределение теплового потока по её длине и резонансной трубе.

В данном примере с помощью математической модели, учитывающей удельный тепловой поток, уточняется длина РТ.

Исходные данные для вычислительного расчёта представлены в табл. 1.

Таблица 1. Исходные данные для расчёта

Наименование параметров, Значение

размерность

Мощность, КВт ~ 12

Г орючее бензин

Расход горючего, кг/с 0,00035

Расход воздуха, кг/с 0,0044

Температура в КПГ, К ~ 1200

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

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

Расчётные данные представлены в табл. 2.

Таблица 2. Расчётные параметры распределения температуры по длине камеры сгорания и резонансной трубы

Наименование Значение

Расстояние от торца КС, м 0 0,08 0,14 0,8

Температура стенки, °К 980 1173 852 800

Распределение удельного теплового потока по длине РТ показано на рис. 4.

где а=1 при 0 <х <ЬК

С целью подтверждения вычислительных расчётов авторами была подготовлена экспериментальная модель КПГ (рис. 1). На специальном стенде проведено 5 огневых испытаний. Запуск модели осуществлялся при помощи сжатого воздуха, поступающего из ресивера компрессора. Бензин подавался из бака под действием избыточного давления. Модель КПГ запускалась без замечаний. На времени ~ 4 с агрегат выходил на устойчивый пульсирующий режим. Частота пульсации давления в КПГ во всех испытаниях составляла ~ 130 Гц. После осуществления запуска модели КПГ и выхода на режим подача сжатого воздуха из компрессора прекращается. Камера работает в автомодельном режиме.

Система измерения параметров состояла из четырёх датчиков температуры и двух датчиков давления (один измерял параметр в камере сгорания, второй - в РТ). Осреднён-ные значения теплового потока представлены четырьмя точками (значения температуры снимались с датчиков).

В результате проведённых исследований модели [2] были получены экспериментальные данные распределения теплового потока по длине КПГ (на рис. 4 обозначено цифрой 1) при работе на компонентах топлива «воздух + бензин».

Рис. 4. Распределение теплового потока по длине КПГ и резонансной трубы при работе на компонентах топлива воздух + бензин:

1 - экспериментальные данные, 2 - теоретический расчёт, 3 - КПГ с резонансной трубой

В результате испытаний были получены:

1. Основные параметры КПГ на различных видах горючего.

2. Значения температуры газа в различных точках, что позволило построить распределение теплового потока по длине КПГ с резонансной трубой. При этом предложенная методика численного решения задачи позволяет получить удовлетворитель-ные результаты сходимости с экспериментальными данными порядка 10 %.

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

Библиографический список

1. Нечаев, Ю.Н. Пульсирующий детонационный двигатель - это реальность [Текст] / Ю.Н. Нечаев, А.П. Полев, А Н Тарасов // Вестн. Воздушного флота. - 2003. -№ 4, 2003.- С. 72 - 76.

2. Солодовников, А.В. Аппарат пульсирующего горения - основа комбинированного пульсирующего воздушно - ракетного двигателя [Текст] / А.В. Солодовников, Е.Н. Вышегородцев, В.В. Голубятник. - М.: ЦВНИ МО РФ, 2009.

3. Солодовников, А.В. Исследование пульсирующих камер сгорания и их применение в аэрокосмической технике [Текст] /

A.В. Солодовников, Е.Н. Вышегородцев,

B.В. Голубятник // Вестн. СГАУ. - 2009. - № 3(19). - 4.2. - С. 335 - 344.

4. Быченок, В.И. Теплотехника рабочего процесса в аппаратах пульсирующего горения [Текст] / дис. ... д-ра техн. наук: 05.14.04.: защищена 05.08.03 г., утверждена 10.09.04 г. / Быченок Вячеслав Иванович. -Воронеж, 2004. - 350 с. - Библ.: с 283 - 302.

CALCULATION METHODOLOGY OF INTERMITTENT BURNING CHAMBER MAIN PARAMETERS

© 2012 A. V. Solodovnikov, E. N. Vishegorodcev, V. V. Golubyatnik Serpukhov Military Institution of Missile Forces

Calculation methodology of direct-flow intermittent burning chamber main parameters (V = const) is treated in this paper. The calculation algorithm is given.

Intermittent burning chamber, mathematical model, calculation, testing.

Информация об авторах

Солодовников Алексей Витальевич, преподаватель кафедры Серпуховского военного института ракетных войск, подполковник. E-mail: [email protected]. Область научных интересов: термодинамика и тепловые машины.

Вышегородцев Евгений Николаевич, кандидат технических наук, доцент, начальник кафедры Серпуховского военного института ракетных войск, полковник. Область научных интересов: аэродинамика, газовая динамика, конструкция ракетных двигателей.

Голубятник Вячеслав Васильевич, кандидат технических наук, научный консультант Серпуховского военного института ракетных войск. E-mail: slava [email protected]. Область научных интересов: термодинамика и газовая динамика, конструкция ракетных двигателей.

Solodovnikov Aleksey Vitalievich, the chair teacher of Serpukhov military institution of missile forces, lieutenant-colonel. E-mail: [email protected]. Area of research: thermodynamics and heat engines.

Vishegorodcev Evgeny Nikolayevich, the candidate of engineering sciences, senior lecturer, chair chief of Serpukhov military institution of missile forces, colonel. Area of research: engineer-aerodynamics, gas dynamics, rocket engines.

Golubyatnik Vyacheslav Vasilyevich, candidate of engineering sciences, Serpukhov military institution of missile forces. E-mail: slava 6123@rambler. ru. Area of research: thermodynamics and gas dynamics, design of rocket engines.