Pull to refresh

Что происходит в критическом сечении камеры сгорания ракетного двигателя (ЖРД) и причём тут «скорость звука»?

Reading time 89 min
Views 11K

Газодинамика сверхзвукового сопла Лаваля

Газодинамика сверхзвукового сопла Лаваля.

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

Пришлось выдумать собственный взгляд на атомную физику в целом, чтобы в итоге подойти к пониманию таких казалось бы изученных вопросов, как "Подъёмная сила крыла" и "Газодинамика работы ЖРД".

Теме "Подъёмная сила крыла" посвящены предыдущие семь глав соседней длинной многолетней статьи ( https://habr.com/ru/post/438854/ ), а теперь настала пора соединить её воедино с темой "Модель твёрдого Ядра без электронных оболочек, то есть Ядро равно всему атому по размеру" ( https://habr.com/ru/post/440848/ ).

Именно написание статьи про устройство атома подтолкнуло меня в изучение аэродинамики, чтобы там найти феномены, объяснимые только моей теорией Статического отталкивания в газах (далее СТГ) , что позже привело и к моему пониманию физики работы ЖРД.

Физическое описание работы ЖРД

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

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

Обход противоречий Кинетической теории газов

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

По СТГ молекулы газов висят неподвижно, а их тепловая энергия заключена в сжатом упругом поле отталкивания вокруг них.

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

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

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

Надо ещё особо отметить, что в КТГ на микроуровне не выполняется закон сохранения импульса.

Так в КТГ при столкновении молекул разных масс невозможен перевод части кинетической энергии в тепло при упругом их соударении с сохранением суммарного импульса после столкновения. Ведь тепло - это и есть кинетическая энергия молекул согласно КТГ.

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

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

Из сказанного получается, что исходные допущения КТГ противоречат базовым законам физики-механики в самой своей основе!

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

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

То есть по СТГ возможны варианты разных комбинаций кинетической и тепловой энергии в молекулах газа:

  • летящие в пустоте быстрые и холодных молекулы c сильно растянутым в пустом пространстве слабо напряженным полем отталкивания (именно такие молекулы должны вылетать из идеального ЖРД в вакууме). Замечу, что по СТГ у летящей сквозь космос молекулы газа появляется объективная температура, независимая от системы отсчёта определения её скорости и соответствующей этой скорости кинетической энергии.

  • медленные горячие молекулы, висящие неподвижно в сосуде под давлением газа; -группы горячих молекул, коллективно летящих на высокой скорости внутри сосуда под давлением (спутник в космосе). Подробнее об этом в моей статье «Модель твёрдого Ядра без электронных оболочек, то есть Ядро равно всему атому по размеру» на ХАБР по ссылке и краткая версия с дополнениями о фазовых переходах на Дзен.

Рис.1 Схематическое изображение пограничной зоны контакта разных фаз Жидкость-Газ (Вода-Воздух) в случае Статической Теории существования Газа:
а) Газ при давлении 1 атм ( плотность воздуха 1,2кг/м3) над жидкой водой (расстояние между центрами молекул приблизительно в 10 раз больше, чем в жидкой фазе);
б) Тот же газ над водой при той же температуре, но при давлении 10 атм (плотность 12кг/м3)., при этом расстояние между центрами молекул газа сокращается в 10^0,333=2,15 раз. Такие же расстояния между молекул газа в КС ЖРД, но уже при давлении 100 атм и температуре 3800К.

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

Внедрение СТГ в привычное описание работы ЖРД


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

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

Так в интервал времени дt через сечение площадью S протекает газ со скоростью V и плотностью q.

Тогда расход массы за время дt составит:

дG=дt*V*Sq

В этом выражении произведение дtVS= дW, где дW - это объём элементарных параллельных дисков пространства сопла, в которых содержится постоянное одинаковое количество вещества (см.рис.2).

Толщина этих дисков зависит только от скорости потока V в данном сечении сопла.

Толщина диска дН=дtV

Площадь S этого диска зависит только от геометрии сопла в данном сечении.

Плотность q зависит только от величины дW при постоянном расходе G (кг/с), где:

дG =дt* G

А плотность равна:

q=дG/дW

откуда

дW= дG/ q

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

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

На сводном Графике параметров газа по длине сопла имеется параметр плотности, по которому можно определять объём дW элементарного слоя в каждом сечении (см.рис.2).

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

Существует более точный график из более специализированного учебника по "Теории ЖРД" [1], где графики представлены значительно точнее, а размерные оси оцифрованы (см.рис.3).

Правда в более точном графике (рис.3) отсутствует параметр плотности. Это косвенно демонстрирует наличие неувязок в термодинамических формулах ЖРД, которые при наличии графика плотности стали бы вопиющими.

Рис.2. Классическая обобщённая картинка с графиком из учебника для объяснения принципов работы и газовых состояний рабочей среды в ЖРД, дающая качественное описание с мелкими огрехами по точности, но с дополнительным графиком плотности газа по длине ЖРД.

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

Из величин S и дW можно однозначно получить осевую скорость V в любом сечении ЖРД:


V= дW/(S*дt)= (дG/ q)/(S*дt)= дG/ (q*S*дt)= G/ (q*S)

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

Так получается, что элементарный газовый объём (содержащий одну и туже порцию вещества) в работающем РД имеет непостоянный объём и непостоянную форму при перемещении по оси РД. (см.рис.4)

Непостоянство формы элементарного газового объёма по длине ЖРД (см.рис.4) предполагает изменение соотношения R/дН. Это соотношение R/дН определяет количество единичных молекулярных столбцов в диаметральном сечении к количеству единичных молекулярных дисковых слоёв в объёме дW при равномерном распределении молекул в объёме.

Рис.4. Распределение толщины слоя дW =дt*V одной доли вещества по осевому ходу в прямоточном РД:
А. Плоские параллельные слои с толщиной, пропорциональной скорости на оси (понятное соотнесение с графиком V на рис.1).
Б. Сферические слои с толщиной, пропорциональной скорости на оси. Выпуклости слоёв направлены в сторону доминирующего из действующих на арку усилий сжатия (избыточное давление в КС и ускорение разгона молекул газа в сопле Лаваля)
.

Если считать, что внутри элементарного объёма dW послойное расположение молекул газа не меняется, то геометрическое изменение сечения сопла будет вызывать такое же изменение расстояния между молекулами в единичном слое, при этом расстояние между слоями будет связано с осевой скоростью потока.


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

Подобный скоротечный адиабатный режим осуществим и при разгоне газа до околозвуковой скорости, что и наблюдается вблизи критического сечения в КС ЖРД.

Именно этот вопрос о рекомбинации молекул в единичном объёме dW по длине КС и сопла ЖРД в зависимости от скорости потока будет предметом подробного рассмотрения в данной статье.

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

Рассмотрение поведения отдельной молекулы газа в её персональном объёме по СТГ

Рассмотрение поведения отдельной молекулы газа в её персональном объёме по СТГ.
Изменение скорости – это ускорение.

Ускорение — это следствие действия силы на молекулу газа с некой массой (разница противоположно направленных сил = равнодействующая сила).

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

Таким образом мы ввели неравновесную плотность газа между соседними рядами молекул в привязке к направлению изменения скорости перемещаемого объёма :
-вдоль оси ЖРД будет плотность qv,

-поперёк ЖРД будет плотность qs.

Также мы выявили разницу давлений на границе слоя, инициирующую ускорение движения данного слоя по оси РД Возникает вопрос:

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

Ответом служит известны феномен "Звуковой волны", где звук распространяется в параллельных слоях, меняя послойные расстояния вдоль направления распространения звука, при этом не меняя расстояния между молекулами вдоль фронта волны по слоям (перпендикулярно к направлению распространения звука).

По КТГ мы разделяем по ортогональным направлениям кинетическо-тепловую энергию молекул газа, выделяя только 1/3 на направление перемещения звуковой волны. Именно эта 1/3 тепловой энергии в одном направлении определят ограничение скорости звука. То есть кинетическая энергия на скорости звука должна быть равна 1/3 общей тепловой энергии молекулы газа выраженной в кинетической энергии по одному направлению. Из этого по законам механики следует равенство :

Vа^2*Mr/2=(1/3)*Emr= R*T

Откуда: Vа^2=2*R*T/Mr

Vа=(2*R*T/Mr)^0,5

В классической «кинетической теории газов» вместо множителя 2 в правой части присутствует показатель адиабаты 1,4 ( для двухатомного газа и воздуха 7/5=1,4). То есть наше предположение отличается от существующей теории на величину Ктеор=(2/1,4)^0,5=1,195 Разница в скорости составит 1,2 раза или на 20%. Причём наша расчетная скорость выше «классической» расчётной скорости звука. Тут сразу есть некоторые дополнения. Обычно на бытовом уровне в качестве скорости звука воспринимают скорость движения фронта волны, что и определяется обычно в экспериментах при измерении скорости звука по отставанию звука выстрела от вспышки выстрела. Но в случае передачи плоской волны в упругой среде максимальная скорость отдельной молекулы будет отличаться от средней скорости фронта волны.

При идеально упругом колебательном движении это отличие составляет:

Vмах=1,41*Vср

Где 1,41= 2^0,5.

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

    Так именно на такую величину 1,41 отличается действующее напряжение U= 220В в электрической сети переменного тока от максимальной амплитуды напряжения Uмах=310В в синусоиде (см.рис.5):

    Uмах =220*1,41= 310 В

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

У обывателя на слуху больше величина U=380 В трёхфазного напряжения (см.рис.6), где:

U=380= 220*(3^0,5)= 220*1,73

При этом есть ещё интересные цифры в электротехнической среде :

230*1,73=400В

Цифра 230В и 400В (или 0,4кВ) также часто встречается в технической среде электротехников (см.рис.7).

Рис.5 Синусоида действующего напряжения электрической сети 220В с максимумом на 310В.

Рис.6 Векторное взаимодействие фаз трёхфазного электроснабжения 3х380В.

Рис.7 Распределение напряжения по сети электроснабжения 220В по длине сети с учётом потерь в проводах.

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

Экспериментальные данные в нормальных земных условиях дают скорость звука 331м/с.

Откуда максимальная скорость по направлению будет:

Vа-max=331*1,41=468м/с.

И эта оценка для Vа-max= 468м/с даёт практически точное попадание в полную тепловую энергию газа в результате расчёта по привычной формуле

Vт=(3*R*T/Mr)^0,5=484м/с

с погрешностью всего в 3,5% для молекулярного веса воздуха Mr=0,029кг/моль:

484/468=1,035 или 3,5% разницы.

В ниже приведённой картинке из Википедии (см.рис.8) тепловая скорость молекул воздуха заявлена Vт=464м/с при +20С, что вообще точно совпадает с полученной нами скоростью Vа-max= 468м/с!

Хотя в статье из Википедии не уточняется по какой именно из формул эта скорость получена и с какой принятой молярной массой газов воздуха (воздух- это смесь газов).

Рис.8. Страница из Википедии с набором формул по КТГ для различных тепловых скоростей молекул и таблица с непосредственно цифрами этих скоростей для разных газов.


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

Так что в величине скорости звука в воздухе сомневаться не приходится.

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

Из приведённых расчётов получается , что при звуковых колебаниях на одно направление полностью срабатывается вся тепловая энергия (кинетическая энергия), а не только статистическая доля в 33% по этому направлению.

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

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

Получается, что энергия поперек скорости звука в слое вообще не расходовалась и давление в слое в направлении вдоль слоя не изменилось!

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

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

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

Но тогда за счёт какой энергии происходит разгон в сопле Лаваля до сверхзвуковой скорости после преодоления критического сечения ЖРД?

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

И это предположение подтверждается энергетическим расчётом для газов в КС ЖРД!

Необходимо найти "пропавшую" часть энергии газа, которая покидает КС в огромном количестве (более 70% общей кинетической энергии струи из ЖРД) и срабатывается в сопле Лаваля на почти 2-х кратное ускорение газов.

Сравнение топливного энергетического баланса ЖРД с теоретическим газодинамическим балансом по КТГ.

Сравнение топливного энергетического баланса ЖРД с теоретическим газодинамическим балансом по КТГ.

Расчёт по химической энергии топлива

Расчёт№ 1

Зайдём к вопросу энергобаланса ЖРД с тепло-механической стороны.

Мы знаем, что в ЖРД сжигают 1 кг керосина с теплотой сгорания 42 МДж /кг и температурой газов 3800К, при этом расходуется 2,7 кг кислорода.

В результате получается струя газов с общей массой m=1+2,7=3,7кг и со скоростью истечения 3400м/с при температуре 1800К (см. граф .рис.3.).

Посчитаем баланс теплотворной энергии топлива и кинетической энергии у струи газов из ЖРД.

Ек= V^2*m/2=3400^2*3,7/2= 21 386 000Дж= 21,38 МДж


То есть баланс сошёлся практически идеально, учитывая остаточную энергию температуры струи 1800К от исходной 3800К в КС (см.рис.1-б), что даёт остаток энергии в струe по температуре 1800/3800=0,473 или 47,3%.

При этом кинетическая энергия составляет 52,7% от тепла сгорания топлива, что совпадает с тепловым КПД (2000/3800)=52,67%.

Етоп=21,38/0,527= 40,57 МДж

Разница между 42 и 40,57 МДж составляет всего около 3,5%, что частично списывается на нагрев и испарение топлива (керосин) и окислителя (жидкий кислород). Также присутствует погрешность замера температур и скоростей (как на ЖРД так и на графике рис.3.), что не позволяет получить 100% точность расчётов.

Реакция горения проходит в газообразной фазе обоих компонентов при постепенном испарении керосина и кислорода с поверхности микрокапель распылённых в объёме КС жидких компонентов топлива, при этом энергия на испарение вычитается из энергетического выхлопа реакции окисления.

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

Теплота испарения керосина Q=240кДж/кг=0,24 МДж/кг

Теплота испарения кислорода Q=3,41 кДж/моль

Испарение 2,7кг кислорода О2 с молярной массой 0,032кг/моль потребует затрат энергии

Еисп-о=3,41*2,7/ 0,032=287,7 кДж=0,288 МДж

Что составит

Еисп-топ/Етоп=(0,288+0,24)/42=0,0125 или 1,25% энергии сгоревшего керосина.

Данным расчётом теплоты испарения мы нашли 1,25% энергии сгоревшего топлива из посчитанной величины 3,5% невязки в балансе по расчёту. Оставшийся 2,25 % уходит на тепловые потери сквозь стенки ЖРД и на работу турбонасосного агрегата, нагнетающего топливо в КС под высоким давлением.

То есть при анализе режимов работы ЖРД вполне допустимо использовать привычные формулы механики по Ньютону.

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

Vмах=(Ет2/Мтоп)^0,5=(40 570 0002/3,7)^0,5=4682м/с

Эта цифра Vмах=4682м/с нам ещё пригодится.

Проверка тепло-механического энергобаланса по КТГ.

Пересчитаем тепловую энергию топлива на один моль газа в смеси СО2+Н2О в пропорции 1:1 по объёму (числу молей).

Из чего можно рассчитать среднюю молярную массу смеси газов ( СО2+Н2О):

((12+162)+(12+16)) /2= (44+18)/2=62/2=31 г/моль=0,031кг/моль

Кинетическая энергия струи из ЖРД составит

Ек-mr=21,38*0,031/3,7=0,179МДж/моль


Теперь проведём тот же расчёт для начального состояния и конечного состояния газов по уравнению полной энергии для 1 моль газа:


Eтmr=3*R*T

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

На срезе сопла поток имеет кинетическую энергию от скорости Vсоп=3400м/с и тепловую энергию от температуры 1800К, которые в сумме должны давать энергию сгорания топлива в КС:

Eжрд=Екин+Етmr

Ежрд=Vсоп^2*Mr/2+3*R*Tсоп =3*R*Tкс

Откуда нам надо получить равенство энергий:

Vкс^2*Mr/2 =3*R*(Tкс-Тсоп)

Посчитаем отдельно левую и правую части по имеющимся данным:

Vкс^2*Mr/2=3400^2*0,031/2=179 180 Дж

3*R*(Tкс-Тсоп)=3*8,3114*(3800-1800)=49 868 Дж

Мы получили отличие энергий в 179180/49868=3,59 раза, что крайне далеко от совпадения с теорией.

Тогда реальный расчётный баланс ЖРД будет выглядеть так:

Vкс^2*Mr/2 =(3*3,59)*R*(Tкс-Ткр)

То есть теперь как-то надо объяснить появление нового коэффициента Кр=3,59*3=10,77 в газовой термодинамике ЖРД.

Проблема где-то есть, и она находится где-то в недрах как самой КТГ, так и в термодинамической модели ЖРД.

Ниже приведена страница учебника по теории ЖРД (см рис.8-9), где в формуле (2.4) вместо коэффициента 3 или 10,77 приводится дробь из показателя адиабаты:

Кр=(к/(к-1))=(1,15/(1,15-1)= 7,667.

Вполне возможно, что для данного ЖРД больше подходит показатель адиабаты к=1,1
что даёт: Кр=(1,1/(1,1-1)=10,8

Вывод по расчёту энергобаланса ЖРД:

Из приведённого расчёта видно, что учебник даёт вполне точное согласование тепло-механического баланса кинетической и тепловой энергий на срезе сопла Лаваля с теплотворной способностью топлива в КС, но с подгоном показателя адиабаты к показателю Кр(к) по топливному энергобалансу.

Правда при к=1,1 и адиабатическом расширении от 100атм при 3800К до 1 атм на сопле Лаваля температура оказывается не опускается до необходимой нам Тсоп=2200К, а оказывается на уровне 2500К, что сильно выше необходимого по термобалансу.
Получается, что для расчёта энергосодержания нам подходит коэффициент к=1,1, а для расчёта адиабатического расширения по конечным точкам температуры от 3800К до 2200К уже нужен к=1,134

Стоит заметить, что для обычного воздуха с показателем адиабаты к=1,4 подгонный коэффициент энергосодержания Кр=1,4/(1,4-1)=3,5.

То есть даже для обычного воздуха разрыв теории ЖРД с другой обычной теорией из учебника КТГ составляет 3,5/3= 1,166 раза или 17%, что довольно много для игнорирования.
Кстати, в учебнике по теории ЖРД (рис.9) ещё и не верно представлены размерности в формуле (2.4), так как в кинетической энергии почему-то потеряли массу одного моля летящего вещества Mr (молярную массу).

Но это уже мои "придирки" к редактору учебника...)))

Рис.9. Расчётные формулы из учебника.

Рис.10. Расчётные формулы из учебника.

Рис.11. Расчётные формулы из учебника.

Доли импульса тяги ЖРД приходящиеся на критическое сечение КС и на сопло Лаваля.

Прикладные аспекты создания РД- это область знаний, где за последние 80 лет развития ракетостроения было измерено и посчитано всё, что только можно измерить и посчитать.

В учебнике теории ЖРД [1] приводятся совершенно в явном виде цифры, которые говорят об отличии силы тяги РД в критическом сечении выше расчётной по давлению в КС на совершенно конкретную величину 22,8-24,8% (см.таб..рис.12).

То есть тяга РД на границе критического сечения без закритической части превышает в 1,234 раза тягу Fкс:

Fкр=1,234Fкс

При этом Fкc рассчитана по формуле:

Fкс=Sкр*Ркс.

Где Fкс- тяга от статического давления дна КС на площади критического сечения, Sкр- площадь критического сечения КС, Ркс- давление в КС на дне (равно давлению наддува и впрыска камеры).

Рис.12. Фрагмент текста учебника с таблицей, где написано (выделено жирном в третьей строчке таблицы 3.1), что тяга РД на границе критического сечения без закритической части превышает в 1,234 раза тягу КС Fкс=Sкр*Ркс.

Бытовая логика заставляет нас считать, что скоростной напор струи из сосуда под давлением должен быть равен давлению в сосуде. И это действительно так для несжимаемой жидкостей!

Но сила реактивной тяги струи вовсе не давлением скоростного напора струи определяется, а импульсом массы выбрасываемой струи!

Сила тяги струи из сосуда но закону Ньютона F=m*a=d(mv)/dt оказывается аж в 2 раза больше, чем произведение площади отверстия на давление скоростного напора:

Fр=Sотв*Ркс

Fтяги=2*

Этот парадокс физики с избытком тяги в отверстии сосуда под давлением достаточно хорошо известен ещё из курса физики-гидравлики. (см.рис.13)

Парадокс о превышении тяги от струи из отверстия над нескомпенсированным давлением жидкости на площади отверстия

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

Это парадокс описывается так:

Через круглое отверстие в стенке сосуда вытекает струя воды. Сила тяги импульса силы этой струи равна:

Fа=m*a=d(mV)/dt=G*V

При этом на площади проекции отверстия действует не скомпенсированная сила:

Fp=S*P=S*q*g*h.

Парадокс в том, что сила тяги импульса струи Fа превышает не скомпенсированную силу от давления жидкости на площади сечения отверстия Fp конкретно в 2 раза:

Fа : Fp=2:1

Выводится это следующим образом.

Массовый секундный расход жидкости G=V*S*q

Подставляем выражение для V и G в формулу для Fа:

Fа=G*V =(V*S*q)*V= S*q *V^2

Скорость струи и давление подчиняется равенствам:

Рh=q*g*h
Pv=V^2*q/2

Где Рh- давление столба воды высотой h с плотностью q при ускорении свободного падения g, Pv- скоростной напор той же воды. Так как в отверстии верно равенство Pv=Ph, то верно равенство:

q*g*h=V^2*q/2

Откуда при сокращении в обеих частях одноимённых величин плотности (q) получаем:

g*h=V^2/2

откуда V^2=2*g*h

Подставляем выражение для V^2 в формулу Fа и получаем:

Fа= S*q V^2=S*q*(2*g*h)=2*(q*S*g*h)

На противоположную от отверстия стенку бака действует гидростатическое давление:

Р=q*g*h.

При этом на площади проекции отверстия действует нескомпенсированная сила гидростатического давления:

Fp=S*P=S*q*g*h.

ИТОГО:

Реактивный импульс силы струи Fa равен :

Fа= 2*(q*S*g*h)

Сила Fp от давления на площади отверстия равна:

Fp=(q*S*g*h)

Откуда: Fа:Fp = 2:1

В рамках гидравлики данный парадокс толком качественным образом не объясняется , сводя всё к расчётам с использованием коэффициентов местных сопротивлений на входе потока жидкости в отверстие в стенке бака.(см.таб.Рис.14)

Рис.13. Схема истечения жидкости из отверстия в стенке сосуда.

Рис.14. Таблица характеристики местного сопротивления на входе в отверстия в резервуаре из-под уровня жидкости.

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

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

В реальной гидравлике существуют сопротивления и потери, поэтому величина отличия двух величин сильно зависит от формы краёв отверстия. Так разная форма захода потока в отверстие снижая эту чёткую двойку до диапазона значений Fa/Fp=1,3-1,96 (см.таб.Рис.13)


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

Максимальный коэффициент потерь 0,5-0,53 при острых кромках на входе в трубу насадка (первые две строчки), с учётом которых получается следующая величина скоростного напора в сечении отверстия:

Рh=0,53*Pv+Pv=(0,53+1)*Pv=1,53*Pv

Откуда: Pv=Рh/1,53=0,65*Ph

То есть потери входа составляют 35% от давления в сосуде, при этом Fa/Fp=1,3 .


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

Pv=Рh/1,02=0,98*Ph

То есть потери входа в коноидальный насадок (см.рис.13,15-16) составляют всего 2% от давления в сосуде.

Именно две последние строчки дают наибольший разрыв в отношении Fa/Fp=1,96

Попытка решить парадокс за счёт учёта местного сопротивления на заужение струи в отверстии (см.рис.13) не является решением парадокса, так как коноидальный и конический заходы дают столь малые потери напора на входе, что обеспечивают практически 2-х кратный разрыв (без 2%) в расчётных значениях "тяги сопла" в явном виде.

Решение парадокса об избытке тяги в струе из бака под давлением

Так откуда берётся избыток тяги в сосуде с дыркой и в ракетном ЖРД?

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

Если отказаться от уравнения Бернулли и принять объяснения с центростремительными ускорениями и реакциями закруглённых стенок на протекающий мимо высокоскоростной изогнутый поток жидкости или газа (см.главу №1 данной статьи) и при согласии на существование разных давлений вдоль и поперёк потока (стакан с вращающимся чаем в гравитационом поле см.рис.17-18), то "чудо" избытком тяги рассеивается как в случае с дырявым баком воды, так и в случае ракетного сопла.


Рис.17. Форма поверхности в стакане с вращающимся чаем.

Рис.18. Разложение сил на поверхности в стакане с вращающимся чаем.


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

Надо помнить, что коноидальный насадок в потоке жидкости не создаёт эффективную форму потока, а скорее её сопровождает.

Таким образом поток жидкости над коноидальным насадком имеет нулевое давление. Возможен даже режим пролёта потока жидкости над коноидальным насадком с тонким воздушным зазором (кавитация).

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

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

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

Пи*(R+r)^2=Пи*2*R^2

Откуда :

r=(2^0,5)R-R=(2^0,5-1)R=0,414*R или приблизительно r=0,2*D


То есть радиус скругления идеального коноидального захода сравнительно мал и составляет всего 20% диаметра основного отверстия.

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

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

Загадка "вертушки Фейнмана"

Удвоенная тяга от реактивной струи при вытекании из отверстия в сосуде под давлением комплекте с формулой подъёмной силы крыла от изгибаемого потока позволяют решить загадку "вертушки Фейнмана" в режиме всасывания.

Рис.19. Вращение газонного разбрызгивателя типа "вертушка Фейнмана" в прямом направлении.

Вращение вертушки-разбрызгивателя в прямом направлении с выбрасыванием воды (см.рис.19) легко объяснимо, так как на конце трубки струя воды создаёт реактивную тягу:

Fт=G*V=(q*S*V)*V=q*S*V^2

Но почему вертушка перестаёт вращаться в режиме всасывания при погружении в воду?

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

Fs=2*S*Pv=2*S*(q*V^2/2)=q*S*V^2=G*V

При этом сила Fs направлена вдоль потока, но против скорости потока внутри трубки.

То есть мы нашли силу, которая должна была бы толкать трубку в обратную сторону!

Так почему трубка вместе с вертушкой не вращается?

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

Рис.20. Схема потоков воды, сил и давлений в "вертушке Фейнмана": А- при прямом потоке (разбрызгивание); Б- при обратном потоке (всасывание).

Расчёт этой силы Fт в качестве "подъёмной силы крыла" производился в первой главе статьи по "аэродинамике".

Итог расчёта таков:

Рт= dH*q*а=dH*q*(V^2/R)

Fт= Рт*B*R

Где dH-толщина изогнутого потока, В-ширина изогнутого потока, R-радиус кривизны загиба трубы. V-cкорость воды в трубе, q- плотность воды в трубе.

Причём для рассматриваемой трубы произведение dH*B=S трубы.

Подставляя первую формулу для Рт во втору для Fт получаем:

Fт= Р*т*B*R=dH*q*(V^2/R)*B*R

Сокращая одноимённые R и подставляя выражение S получаем:

Fт= dH*q(V^2/R)*B*R = dH*В*q*(V^2) = S*q*(V^2)

Так как G=S*q*V , то подставляем в формулу и получаем

Fт= S*q*(V^2) = ( S*q*V)*V = G*V

То есть для реактивной центростремительной силы от потока в трубе при повороте мы получили выражение Fт=G*V , тождественное к формуле реактивной тяги струи из отверстия.

Вывод по вертушке Фейнмана: Оказывается, что реактивная тяга струи создаётся вовсе не на срезе отверстия, из которого струя вылетает, а где-то в недрах трубы на каких-то загибах водяного потока. И это весьма нетривиальный результат!

ИТОГО:

При обратном направлении работы вертушки в режиме всасывания Fт=Fs (по модулю) и направлены противоположно, взаимно компенсируя друг друга, что и не даёт вертушке крутиться.

В прямом направлении вращения вертушки тяга Fs трубки тоже есть (нагнетание в трубку из шланга), но находится она на оси вращения вертушки (где шланг крепится).

Таким образом силы от нагнетания воды в трубку от оси вертушки не создаёт крутящего момента на трубке относительно оси вращения вертушки, так как при нагнетании воды в трубку сила Fs направлена по радиусу от оси вращения вертушки.

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

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

Гидравлические потери по длине трубки не учитываем.

Это позволяет записать Р=0 в трубке при прямом вращении (выбрасывание струи). При этом поперечное динамическое центростремительное давление на изгибе трубки носит локальный характер и не оказывает влияние на течение жидкости вдоль прямой трубы.

Избыточное давление P в вертушке целиком расходуется на первоначальный разгон воды на входе в трубку (см.рис.17-А) , при этом: Р=Рv , где Pv=q*V^2/2.

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

Р= - Рv

Таким образом мы получили чёткое силовое обоснование отсутствие вращения вертушки Фейнмана в режиме всасывания.

Вопрос:

А для чего вообще мы рассматривали "вертушку Фейнмана"?

Ответ:

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

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

Величина прибавки тяги КС ЖРД от коноидально-конического захода на критическое сечение

Разобравшись с механизмом работы коноидального насадка на выходе потока из бака под давлением можно продолжить разбор работы КС ЖРД.

На срезе критического сечения КС (и срезе сопла) начинает действовать "великий интегратор"- закон импульса силы Ньютона:

F=m*a= m*dV/dt= d(m*v)/dt.

То есть на срезе критического сечения КС становится абсолютно всё равно, что происходило с газами внутри самой КС, а тягу определяет однозначное произведение массового расхода на реальную продольную скорость истечения этих газов

Fкр=G*Vкр.

При выстраивании силового баланса в работе ЖРД нужно куда точнее измерять скорость струи в критическом сечении КС, чем это отображается на графиках параметров газа в тракте ЖРД (см.рис.2-3). На критическом сечении графики давления и температуры становятся почти вертикальными, теряя точность в привязке к критическому сечению.

Удивительно то, что в КС ЖРД даже при очень плавном заходе потока на критическое сечение у потока обнаруживается всего 22-25% прибавки силы тяги струи, вместо ожидаемого 2-х кратного роста тяги, как в случае с несжимаемой жидкостью при истечении из коноидального насадка на отверстие в стенке бака под давлением.

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

К тому же скоростной напор струи газов в критическом сечении даже не создаёт запирающего давления, равного давлению в КС!

Скоростной напор Pv на срезе КС составляет всего 62% от давления Ркс на дне КС.

Это следует из табличных данных (см.рис.11), где сила тяги по импульсу силы:

Fa =124% от .

Нам известно, что 24% приращение тяги связано с двукратностью импульса силы Fa по отношению к силе Fр-кр по площади критического сечения от реального динамического давления Pv-кр (скоростного напора струи газов) на срезе КС.

Следовательно избыточная сила тяги по импульсу силы Fa =124% от Fрv-кр означает, что скоростной напор в критическом сечении находится всего на уровне 124/2= 62% от возможного по имеющемуся давлению Ркс в КС.

0,62*Ркс=qкр*Vкр^2/2

Для сохранения силового баланса в КС где-то надо найти потерянные 38% от давления Ркс на выходном конце КС.

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

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

Расчёт скорости газов в критическом сечении КС исходя из известной плотности газов и известного импульса тяги ЖРД в критическом сечении КС.

Расчёт № 2

Просчитаем данные по КС и критическому сечению по данным из таблицы в учебнике ( рис.12) и формулам из учебника (рис.11).

Рис.11. Расчётные формулы из учебника.

Посчитаем расширение газа от зоны горение топлива в КС до критического сечения через уравнение адиабаты:

Ркс/Ркр=(qкс/qкр)^1,15

qкр=((qкс)^1,15*Ркр/Ркс)^(1/1,15)

qкр=((9,84)^1,15 * 5740000 / 10000000)^(1/1,15)=6,0723кг/м3

То есть данные зависимости из учебника (рис.10) вполне воспроизводимы по расчётным формулам.

Подставляя известное значение давления в КС на уровне Ркс=100атм =10МПа получаем:


Ркр=0,574*100=57,4атм

Температуру в КС считаем 3800К, что даёт температуру в критическом сечении:

Ткр=0,93* 3800 = 3534 К

Плотность газов в КС получаем из молярного объёма газа при давлении 10МПа и температуре 3800К:

Vm=RT/P=8,31* 3800 / 10000000=0,00315 м3/моль

При молярном весе 0,031кг/моль для смеси газов (СО2+Н2О) в равных объёмах плотность составляет :

qкс=0,031/0,00315=9,84 кг/м3
qкр=0,617*qкс=0,617* 9,84= 6,07 кг/м3

Давление в критическом сечении равно:

Ркр=0,574*100=57,4атм

Зная плотность газов в критическом сечении КС можно получит скорость Vкр из уравнения тяги:

1,24*Ркс=qкр*Vкр^2

Откуда:

Vкр=(1,24*Ркс/qкр)^0,5=(1,24*10000000/6,07)^0,5=1429м/с

Итого по расчёту №2:

Мы получили скорость в критическом сечении

Vкр=1429м/с при плотности газов qкр=6,07кг/м3.

Полученная скорость Vкр=1429м/с прекрасно укладывается в предложенную ранее модель с отличием скорости фронта звуковой волны от максимальной скорости молекул во фронте волны.

Так при расчётной по учебнику скорости "звука" в критическом сечении КС:

Vа=(1,15*8,31* 3534 /0,031)^0,5=1043м/с

Vкр-м=Vа* 1,41=1043 * 1,41=1470м/с

То есть Vкр-м отличается от Vкр в 1470/1429=1,029 раза или на 2,9%.

Если посчитать плотность газов в критическом сечении КС по скорости звука Vа=1043м/с, то плотность получится несуразно большой:

qкр=1,24*Ркс/( Vкр^2)=1,24 * 10000000/(1043^2)=11,39 кг/м3

Получилось, что плотность газов в критическом сечении qкр= 11,39 кг/м3 превосходит плотность газов qкс= 9,84 кг/м3 при Ткс=3800, что невозможно при разгоне газов в критическом сечении КС за счёт своего расширения.

Таким образом расчётная "скорость звука" никоим образом не применима для расчёта тяги в критическом сечении КС ЖРД.

Первый крупный смысловой разрыв в теории ЖРД

Расчёт №3

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

Несколькими главами ранее мы посчитали предельную скорость газовой струи при полном переводе тепла в кинетическую энергию.

Для керосин-кислородного топлива это составило:

Vмах= (Ет*2/Мтоп)^0,5 = (40 570 000*2/3,7)^0,5 = 4682 м/с

Правда, в учебниках приводится цифра скорости 3599м/с при равновесном расширении для кислород-керосиновой пары с отношением (3,07:1). Но я не готов принять эту цифру, так как кинетическая энергия такой струи составляет всего 65% от теплотворной способности керосина при сгорании.

Согласно таблице балансов тяги по ЖРД (см.рис.12) на КС приходится чуть менее 50% (точнее 100%*1,234 / 2,5= 49,36%) от конечной теоретической скорости потока из сопла бесконечного расширения при показателе адиабаты к=1,15.

На основании этого струя из КС должна иметь скорость:

Vкс-мах =Vмах*0,4936 = 4682*0,4936 = 2311м/с

Если в реальном ЖРД (см.рис.3) на срезе сопла скорость 3400м/с (в вакууме), то доля тяги КС реального ЖРД составит:

Fкс / Fжрд = 2311/3400 = 0,679 или 68%

А доля энергии срабатываемой на КС составляет:

Ек-кс/Ек-жрд =(Vкс/Vжрд)^2=(2311/3400)^2=0,462 или 46%

Из этого энергобаланса становится понятно почему развитой сопловой частью стали оснащать ЖРД только на крупных орбитальных РН, где по формуле Циолковского борьба за прибавку скорости газов из сопла ЖРД имело первостепенную важность:

Vрн= Vжрд * ln(M/m)

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

Критерием для такого снижения КПД РД тактических ракет (чаще твёрдотопливных, а не жидкостных) является снижение общего габарита ракеты, чтобы помещалось больше боеприпасов в ограниченном объёме отсеков самолётов или транспортных машин, а сами ракеты стоили дешевле в производстве. (см.рис.22-23)

Рис.22. Разрез реактивного неуправляемого снаряда. Под цифрой (6) указан реактивный двигатель, состоящий из следующих элементов: конический заход из КС на критическое сечение, короткое цилиндрическое критическое сечение, коническое сопло. Интересно , что коническое сопло всего в 2 раза шире критического сечение (Fa/Fкр=4-4,5), а параллельная работа нескольких малых РД позволяет сократить длину РД в целом при неизменной его мощности (силе тяги).

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

Пересчитаем параметры ЖРД под новую скорость Vкр=2311м/с.

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

1,24*Ркс=qкр* Vкр^2

Откуда:

qкр = 1,24*Ркс /( Vкр^2 )= 1,24* 10000000 / (2311^2) = 2,321 кг/м3

Из полученной плотности можно посчитать давление в струе газов на срезе КС по уравнению адиабаты:

Ркр = 10/(9,84/2,321)^1,15=1,899 МПа

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

Ткр=1899000*(0,0224*1,38 / 2,31)/8,31=3058К

То есть в КС срабатывается дТкс= 3800-3058=742К

Или Екс/Ежрд=742/(3800-1800)=0,371 или 37% полезной кинетической энергии из ЖРД, при Vсоп=3400м/c.

Опять проявилась расчётная неувязка, если учесть что ранее была получена Екс/Ежрд большей величины, а именно:

Екс/Ежрд= (2311/3400)^2 = 0,46 или 46%

То есть температура газов из КС должна быть ещё ниже, а именно

дТкс=0,462*2000=920К.

Тогда:

Ткр=3800 - 920= 2880К.

Если мы не сомневаемся в полученных значениях плотности газов qкс=9,84 и qкр=2,32кг/м3, то следовательно требуется корректировка величины давления Ркр.

При Ткр=2880 К давление будет

Ркр=2880* 8,31/ (0,0224* 1,38 / 2,31)=1 788 456 Па =1,79 МПа

путём подгона коэффициента адиабаты получаем нужное равенство при к=1,1914 :

Ркр = 10 / (9,84/2,321)^1,1914 =1,789 МПа

Интересно, что нужную нам Ткр=2880К по энергетическому балансу в КС ЖРД можно получить при расчёте температуры по формуле из учебника с показателем адиабаты к=1,639 :

Ткр=Ткс*2/(к+1)=3800* 2 / (1,639+1)= 2880 К

В случае работы ЖРД в плотных слоях атмосферы с давлением на срезе сопла Рсоп=0,1Мпа (1атм)

Скорость газов из сопла составит Vcоп= 3100м/с (см.рис.3).

А баланс тяги в ЖРД составит

Fкс / Fжрд=2311 / 3100 =0,745 или 75%

Это вполне попадает в привычные 60-70% в соотношение для распределения тяги между частями ЖРД при работе с подпором атмосферного давления.

В формуле постоянного массового расхода G плотность и скорость друг с другом жёстко связаны:

G=q*V,

При этом саму величину G определяет сам двигатель в исходных данных, и мы его изменить не можем, не нарушив условие неразрывности потока.

G=qкр*Vкр=6,07 * 1429=8674 (кг/(с*м2))- по формулам КТГ из учебника.

G=qкр*Vкр=2,321 * 2311=5363 ( кг/(с*м2)) - по термическому балансу и законам Ньютона

Разночтения в массовых расходах G по разным методикам подсчёта составляют отношение:

8674/5365=1,617 или около 62%.

F=qкр*Vкр^2=6,07*1429^2=12 395 188 Н/м2 - по формулам КТГ из учебника.

F=qкр*Vкр^2=2,321*2311^2=12 395 813 Н/м2 -по термическому балансу и законам Ньютона

При этом силы тяги Fкр у нас совпали (их мы приняли изначально одинаковыми в качестве исходных данных).

Если бы мы знали реальный параметр G для двигателя, то мы могли бы сделать вывод о том, какая из цифр скоростей Vкр ближе к реальности.

Итого по расчёту №3:

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

Расчёт параметров струи в сопле Лаваля.

Аналогичный расчёт необходимо провести для Сопла Лаваля, где нам известен рост скорости в

Кv=Vсоп / Vкр= 3400 / 2311= 1,47 раза при разгоне от 2311 до 3400м/с ;

  • геометрический рост сечения сопла 18,86 раза ( справочно для РД-107);

  • начальная плотность q=2,321кг/м3 на срезе КС.

  • Тсоп=1800К ( с графика рис.3) В итоге мы получим параметры на срезе сопла: q, P

    Исходя из предыдущих выкладок нам известно что при постоянном расходе по сечению G элементарный объём dW имеет связь с параметром скорости и плотности зависимостью:

    G=dW*q=const

    Так как кратность изменения геометрия слоёв dW нам известна (по росту сечения и росту скорости для РД-107), то можем посчитать параметр q среза сопла:

    qcоп=qкс / ((Vcоп/Vкр)* (Scоп/Sкр)) = 2,321 / (1,47*18,86)=0,084кг/м3

    Исходя из полученной плотности можно получить молярный объём Vm, а из него посчитать давление газов, зная конечную температуру Тсоп=1800К:

    Vm-соп=(0,0224 * 1,38 / 0,084)=0,368 м3

    Р=R*Т / Vm = 8,31 * 1800 / 0,368 = 40 646 Па или 40,6кПа = 0,406 атм.

    Откуда проверочный расчёт по температуре сопла:

    Тсоп=PVm/R=40646* (0,0224 * 1,38 / 0,084) / 8,31=1800К

    Теперь зная конечные и начальные состояния газа можно получит показатель адиабаты для сопла Лаваля в РД-107 со стартовым давлением Ркр=17,9 атм на срезе КС , который при подгоне получается k=1,405:

    Ркр=Рсоп * (qкр/qсоп)^1,1405 =0,0406* (2,321 / 0,084)^1,1405 =1,788 атм.

    При подгоне под давление 17,9 атм в критическом сечении КС получен показатель адиабаты к=1,1405 для сопла Лаваля в РД-107, что уже лучше соответствует табличному показателю к=1,15, чем ранее посчитанный для КС показатель к=1,639.

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

    С фактами и физикой-механикой я спорить не готов, а вот к формулам КТГ в быстротекущих неравновесных процессах у меня доверия нет.

    Так что примем для расширения в сопле Лаваля коэффициент адиабаты k=1,14 при Ркр=17,9 атм и Vкр =2311м/с, а затем посмотрим куда это нас приведёт в дальнейшем.

Происходящее в критическом сечении КС и в сопле Лаваля с точки зрения СТГ

Рассмотрим происходящее в сопле Лаваля на уровне отдельной молекулы газа в её элементарном объёме по теории СТГ.

Для расходования остаточного поперечного давления необходимо, чтобы газ совершил работу в поперечном от скорости потока направлении, то есть вдоль элементарного слоя dW.

Но ведь именно поперечное расширение тонких слоёв и создаёт расширяющееся сопло Лаваля в РД!

При этом темп расширения стенок сопла оказывается ниже скорости потока по оси сопла и ниже скорости звук в этом слое газа (угол менее 40 градусов у стенки сопла к оси РД).

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

В этой аэродинамической аналогии получается, что от кромки сопла Лаваля начинает расходится фронт возмущения звуковой волны поперёк потока газа из критического сечения КС в след за отступающей вбок стенкой сопла ЖРД.

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

Этот фронт разрежения газов от кромки критического сечения КС достигнет оси потока газов из КС, образуя конус плотного газа из КС в перевёрнутой усеченной круглой пирамиде разреженного газа в начальной части сопла Лаваля. (см.рис.24-25).

При этом в объёме газа между стенкой сопла Лаваля и конусом плотного газа из КС возникает пространство постоянного низкого давления, в котором частицы газа имеют дополнительную скорость поперёк сопла со скоростью отступания стенки сопла от частиц в потоке струи из КС:

Vст-соп = Vкр * sin(Асоп)

Где Асоп- угол стенки сопла к оси ЖРД.

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

Суммарная скорость молекул газа на границе конуса возмущения составит ( по теореме Пифагора):

Vкон-кс= ( Vкр^2+Vст-соп^2)^0,5

Где Vст-соп- это скорость отступания стенки сопла в радиальном направлении от потока по оси ЖРД на скорости Vкр.

При этом направление разлёта молекул газа в сопле Лаваля со скоростью V4=Vкон-кс будет под углом близким к углу расхождения стенок конической части сопла к оси ЖРД (см.рис.24-25)

Зона постоянного давления в сопле сохраняется на начальной конусной части сопла до момента полного срабатывания конуса постоянного давления газов из КС на оси ЖРД.

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

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

Рис.24. Схема сечения ЖРД с указанием скоростей в разных сечениях.

Рис.25. Схема распределения элементарных объёмов (плотностей газа) и скоростей молекул газа в ЖРД. На рисунке вектора скорости по сечениям показаны в масштабе, который определяет условное соотношение толщин элементарных объёмов dW по тракту ЖРД и скорости потока по сечениям ЖРД. Номерами на рисунке Б обозначены: 1-жидкое топливо; 2-газ в КС в широком сечении с максимальной температурой; 3-газ КС с максимальной скоростью при кинематическом разгоне ( с постоянной плотностью) перед критическим сечением; 4- газ в критическом сечении с Vкр; 5-газ на срезе сопла с минимальной плотностью при максимально скорости Vсоп.

Рис.26. Разрез одной камеры РД-107 (РН типа "Союз") с рабочим давлением 60атм в КС. Хорошо заметно резкое заужение в критическом сечении КС. В данном случае большая по объёму КС выполнена для обеспечения более спокойного и стабильного режима горения топливной смеси на меньших скоростях движения на большем протяжении КС. Интересно, что длина КС больших и маленьких ЖРД на фото одинаковые, тогда как длина сопел пропорциональна поперечному сечению КС. То есть длина КС ЖРД не масштабируется с изменением мощности, так как определяется режимом горения топлива, а не мощностью ЖРД в целом.

Поведение газов в КС на сужении

Поведение газов в КС на сужении.

В реальной КС ЖРД двигателя РД-107 в конической части КС сечение сокращается по диаметру в 2,5 раза от широкой цилиндрической части к самому узкому критическому сечению. Геометрические оценки даны по фотографии с разреза реального РД-107 (см.рис.26).

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

Изменение в 2,5 раза по диаметру соответствует изменение площадей сечения конуса в 2,5^2=6,25 раза в начале и конце сужения КС.

Также нужно учесть изменение плотности газов при разгоне, которая падает от 9,84 до 2,321кг/м3 (согласно ранее проведённым расчётам), то есть в 9,84/2,321=4,24 раза.

При сохранении постоянного массового расхода G по всем сечениям КС получаем, что скорость в коническом сечении возрастает в отношении

V3 / V1 =6,25*4,24 = 26,5 раза

Таким образом при скорости V3=Vкр=2311 м/с получаем:

V1=2311 / 26,5 = 87м/с

Учитывая квадратичность кинетической энергии от скорости и 26 кратный рост скорости от V1 до Vкр получаем, что до сужения КС кинетическая энергия газов в 700 раз меньше, чем конечная кинетическая энергия струи в критическом сечении на скорости Vкр (или Ек-s1=0,15% от Ек-кр).

Можно сказать, что при перепаде дТ= 920 К на струе газов в КС при разгоне погрешность в менее 0,15% или дТ=1,4К можно считать несущественной для описания физики процесса.

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

Скорость жидкого топлива в донной части КС

При плотности жидкого керосина 810кг/м3 и жидкого кислорода 1141кг/м3 средняя плотность жидкой топливной смеси в соотношении 1:2,7 составит

qcр=(810+ 1141*2,7) / ,7=1051кг/м3

Тогда в донной части КС, где впрыскивается жидкий кислород и жидкий керосин со средней плотностью около 1051кг/м3, при таком же постоянном массовом расходе G средняя скорость топливной смеси составит всего

Vо=879,84/1051=0,81м/с.

То есть график скорости потока в КС ЖРД будет практически неотличим от нуля в масштабе изображения (см.рис.3), пока газы не дойдут до горячей зоны КС с максимальной температурой Ткс=3800К и до геометрического разгона суживающейся части КС.

Поведение газов в критическом сечении КС ЖРД

Сопряжение конической части КС и конической части начала сопла Лаваля имеет скругление, а потому рассмотреть это участок разгона газов в ЖРД нужно особо внимательно.

Переход между двумя конусными поверхностями КС и сопла Лаваля осуществляется через тороидальную (бублик) поверхность сопряжения (см.рис.24-25).

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

В коноидальном насадке на коротком участке меняется геометрия слоёв газа: от формы сферических слоёв равных скоростей в конической части КС на плоский слой в критическом сечении.

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

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

При этом при разгоне совершается значительная трансформация внутренней энергии из тепловой в кинетическую с падением плотности при расширении газа.

В критическом сечении будет наблюдаться резкий рост скорости согласно постоянному массовому расходу по сечениям:

V3* S3* q3=V1* S1 * q1

Откуда

V3 = V1*S1* q1 / (S3*q3)

Геометрические величины конфигурации КС будут нами сняты с условной геометрии по фото разреза ЖРД РД-107.

При том, что величину G=V1*S1*q1=V0*S0*q0, можно взять из любого сечения КС.

В данном случае нам будет интереснее величина скорости V2, так как скорость V3=2311м/с при q3=2,321кг/м3 на критическом сечении КС мы уже получили из уравнения фактической тяги КС.

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

Для двигателя РД-107 (см.фото.рис.26) угол конуса составляет около 32 градусов (оценочно).

Отношение диаметров Д1/Д2= 2,2 (оценочно), соответственно отношение площадей:

S1/S2= (Д1/Д2)^2= 2,2^2=4,84.

Откуда V2=V1* S1/S2= 87*4,84=421м/с

Изменение кинетичской энергии между сечениями S2 и S3 cоставляет отношение:

Ек-s2/ Ек-s3=(V2 / V3)^2= (421 / 3211)^2=0,0179 или 1,8%

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

Разгон газов в коноидальном заходе КС

Площадь выпуклой поверхности сферического сегмента будет больше площади плоского круга приблизительно в S2/S2-0=1,11 раза.

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

S2/S3=(Д2/Д3)^2=(92/69)^2=1,26^2=1,59

Где Д2/Д3=92/69=1,26 - это оценочный замер с фото РД-107 (см.рис.26)

Тогда площади сферической поверхности крайнего конического сечения S2 будет давать с площадь плоского критического сечения S3 соотношение:

S2/S3=(S2/S2-0)*( S2/S3)=1,11 * 1,59=1,76

Получается, что на коноидальном сопряжении в КС скорость увеличивается в 1,76 раза только за счёт резкого изменения кривизны и диаметра сечений.

Общее увеличение скорости на коноидальном заходе:

Vкр / V2= 3211 / 421= 7,62 раза.

При этом расчёт последовательного геометрического разгона и последующего линейного расширения даёт отношение :

(S2/S3)* (qкс/qкр)=1,76* (9,84 / 2,321) = 7,46

Получились весьма схожие цифры между 7,62 / 7,46= 1,021 или разница 2%.

Таким образом можно смело сказать, что практически весь разгон за счёт изменения плотности

qкс/qкр= 9,84 / 2,321 = 4,84 раза

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

То есть для коноидального захода от конуса КС на цилиндрическое критическое сечение в расчётной модели вполне можно предположить , что сначала газ при постоянной плотности чисто геометрически разгоняется от V2=421м/с до скорости V3=421*1,76=741м/с.

И эта скорость V3=741м/с даёт долю всего 5,3% от конечной кинетической энергии КС:

(V3/ Vкр)^2 = (741 / 3211 )^2 = 0,053 или 5,3%

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

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

Зона разрежения в КС на коноидальном заходе

Прибавка 24% от Fкс=Ркс*S3

накладывается на узкую полоску вокруг критического сечения, где уже на относительно высокой скорости в интервале 421..741м/с действует отрывное течение искривлённого потока газов.

(R+r)^2=1,24* R^2

Откуда r=((1,24)^0,5 - 1)*R = 0,11*R

То есть 24% прибавки тяги обеспечивается кольцом вокруг отверстия шириной всего 5% от диаметра Д3 отверстия S3 в критическом сечении КС.

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

Так температуру сильно меняют значительные тепловые потоки вдоль стенки ЖРД от горячей КС к «холодному» соплу Лаваля.

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

ИТОГО по разгону газов в коноидальном заходе от КС на критическое сечение:

В связи с мгновенным разгоном молекул в критическом сечении возникает вопрос:

Может тогда и критическое сечение не нужно?

Почему бы не делать КС в виде цилиндра с раструбом сопла на сверхзвуковой части?

Собственно ради этого вопроса, по большей мере, и была затеяна эта статья.

Ответ на вопрос состоит из двух частей:

Ответ-Часть 1.

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

Разгон от скорости Vкс3 до Vкр с одновременным падением СРЕДНЕЙ плотности газа в 4,8 раза осуществляется на очень коротком участке коноидального захода перед критическим сечением. Так в критическом сечении происходит почти мгновенное разуплотнение слоёв потока газа по скорости, при этом поперёк потока плотность молекулярных столбцов молекул газа сохраняется неизменной.

Небольшой участок прямой трубы нужен для полного выравнивания скорости по сечению, чтобы в сопло Лаваля поступил уже поток равной скорости.

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

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

Для лучшего понимания механизма почти мгновенного разгона до Vкр от V3 можно представить себе стенку с мелкими отверстиями, равными элементарной площади одной грани единично- молекулярного объёма газа.

Так получится «молекулярное сито» с тонкими трубками для прохода молекул строго по очереди.

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

Гигантское ускорение такого разгона объясняется весьма существенным перепадом давления на открытой в вакуум молекулярной трубе, где на бесконечно малой длине элементарного объёма Wмол одной молекулы происходит резкое срабатывание давления в КС равное Ркс=60-250атм до Р=0 в вакууме. При этом при бесконечно малой массе отдельной молекулы такой резкий скачёк давления создаёт гигантское ускорение.
Удельная молекулярная сила разгона одной молекулы составляют:

Fм=Ркс*Vмол^(2/3)=Ркс*(0,0224м3 / Nа* q/qмол)^(2/3)

где qмол- плотность газа в нормальных условиях.

При этом масса одной молекулы составляет

Mм= Mr*Vмол * q/qмол= Mr* (0,0224м3/Nа) * q / qмол

Подставляя эти две бесконечно малые величины разных порядков в закон Ньютона F=ma получаем:

а=Fм / Mм= (Ркс*(Vмолq/qмол )^(2/3)) / (Mr *( Vмол* q/qмол))

Сокращая одноимённые множители в числителе и знаменателе получаем:

а= Ркс/ ( Mr * (Vмол* q / qмол) ^(1/3) )

Учитывая крайнюю малость величины:

Vмол=0,0224 м3 / Nа=0,0224 / (6,02*10^23) =3,7 / (10^26)м3

То даже кубический корень из неё останется бесконечно малой величиной:

Lмол=Vмол^(1/3)= 3,34 / (10^9)м

Подставляем эту чудовищно малую величину в выражение для ускорения и получаем гигантское до невообразимости ускорение:

а= Ркс/ ( Mr (Vмолq /qмол) ^(1/3))

a=10 000 000 / (0,031*(3,7 / (10^26)*10/1,38)^(1/3))= 5*10^16 м/с2

Разделив эту цифру на величину g=9,81 м/с2 получим ускорение в более привычных величинах:

а=5* 10^15 g= 5000 * 10^12 g

Словами эта цифра звучит так : « Пять тысяч триллионов Ж»

Ответ-Часть 2.

При смене точки зрения оказывается, что это КС расширяется, а не «выхлопная труба» сужается.

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

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

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

В результате в сильно удлиненной КС тепло из зоны горения не сможет попасть назад в сильно отдалившуюся зону испарения.

В итоге произойдёт срыв пламени, горение в КС остановится совсем, КС захлебнётся не газифицированным холодным жидким топливом, и ЖРД потухнет.

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

Инновационная форма КС.

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

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

Подпор плоской стенки полусферической КС вполне возможно осуществлять за счёт охлаждаемой полости с топливом-охладителем, находящимся за стенкой КС под давлением впрыска, равным давлению в КС. (см.рис.27)

Так в ЖРД малой мощности, где относительная длина КС в составе ЖРД очень велика (см.рис.26.), использование более компактной шаровой формы КС может сильно сократить длину ЖРД, а большая площадь полусферического форсуночного дна КС позволит удобнее подводить большое количество трубопроводов к отдельным форсункам.

Так КС с длинной цилиндрической частью с соотношением сторон 3*Д=L можно заменить на полусферическую КС с диаметром 2Д, или на полную сферу диаметром 1,7Д.

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

Но не экономия места главное качество сферической КС!

В сферической КС нет наружных перегретых стенок! Получается что вся наружная поверхность сферической КС закрыта форсунками и подводящими камерами с холодным топливом.

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

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

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

Рис.27. Сравнение габаритов ЖРД с КС разной формы, но одинакового объёма, при одинаковых массовых расходах топлива G и при одинаковом давлении Ркс:
А- цилиндрическая КС с соотношением Длина/диаметр как 3:1.;
Б- Полусферическая КС;
В- сферическая КС с глубоким входом коноидального захода внутрь сферы КС
.

Форма захода на критическое сечение КС

Участок коноидального перехода оказывает максимальный вклад в формирование дополнительной тяги КС от пониженного давления в сужении КС, нескомпенсированного постоянным максимальным давлениям на дно КС.

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

Рис.28. Варианты захода на критическое сечение КС с плавным коноидальным заходом и с ломаной образующей на конус с малым углом конуса.

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

То есть коноидальный заход выгоднее резкого излома разгонных поверхностей в КС.

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

Такие резкие ломанные переходы от КС к расширяющемуся соплу являются часто встречающимся решением, упрощающим конструкцию малых РД.

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

При переходе от прямолинейного течения конусной части КС на коноидальную поверхность об изменении направления потока узнает не всё сечение струи сразу, а только крайние к стенке КС слои газа. При этом распространение информации по сечению от стенки к оси ЖРД пойдёт со скоростью не выше скорости «звука» для данного сечения при данной температуре. Учитывая относительно небольшую ещё скорость на сечении S2 при высокой температуре, то выравнивание скоростей по сечению создаст не очень сильное запаздывание на участке начала коноидального захода. Так в итоге в критическом сечении (любой формы) формируется поток газов с почти постоянной скоростью по критическому сечению КС.

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

Наклон образующей конуса определяется векторным отношением осевой скорости потока V3 в критическом сечении к скорости звука в критическом сечении V3а.

Так как по нашим расчётам скорость газов V3кр в критическом сечении больше скорости звука V3а приблизительно в 2 раза, то конус плотного потока из КС в сопле будет иметь вид острого конуса.

Фактически коноидальное сопряжение КС и сопла носит технологический характер, снижая ударные нагрузки и тепловые нагрузки (перегрев острой кромки излома стенки) при срыве потока газов с острых кромок сопряжений. Плавное скругление в коноидальном сопряжении конусов КС и сопла обеспечивает размывание фронта звуковой волны поперёк оси ЖРД, тем самым обеспечивая более мягкий и устойчивый режим работы ЖРД.

Если же тороидального скругления в сопряжении конусов нет, а конусы граничат по излому образующих (что допустимо), то скачёк скорости в 1,5 раз на сопряжении КС и сопла происходить практически также, но без предварительного выравнивания потока по критическому сечению и с некоторыми резкими возмущениями при отрыве струи на изломе.

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

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

Разгон газов в сопле Лаваля


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

Но если у массы газа уже имеется максимальная возможная продольная скорость по оси РД в критическом сечении КС, то как её увеличивать в сопле Лаваля?

Здесь можно предположить, что работа слоя при расширении как раз и расходуется на ускорение газа.

Как адиабатное расширении тонких слоёв поперёк сопла даст приращения скорости газа в перпендикулярном к расширению направлении?

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

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

А с чего вообще должны образовываться какие-то арочные-сферические слои в сопле РД?
Как выясняется, наличие таких сферических слоёв подразумевается при построении линий тока в КС и в сопле.

В сужении КС линии тока сходятся в пучёк в критическом сечении, а потом расширяются веером в сопле за критическим сечением.

Линии тока- это по сути траектория движения молекул газа в РД.

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

То есть в направлении перпендикулярном к линиям градиента давления изменение давления будет равно нулю, а направление постоянного давления- это линия изобары (линия постоянного давления).

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

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

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

То есть сферический слой по изобаре в РД и есть та самая арка из молекул под давлением, существование которой мы предположили ранее!!!

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

Так при круговом движении потока возникают центростремительные ускорения, перпендикулярные к вектору скорости. При этом центростремительные силы, направленные по радиусу кривизны, не совершают работы.

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

В расширяющемся коническом сопле постоянное по величине давление действует вдоль изогнутой поверхности сферического-арочного слоя, на каждой из молекул которого при этом создаётся перпендикулярная к поверхности слоя сила (см.рис.30-33).

Рис.30. Схема упругих сил в слоях и между слоями на сужении в КС и расширении в сопле ЖРД. (разрез по оси). На схеме показана неравновесная работа упругих сил взаимного отталкивания в ускоряемых газах в различных частях ЖРД (по СТГ). При разгоне увеличиваются расстояния между арочными слоями в начале и в конце заужения КС на схеме. В сопле Лаваля происходит разгон газов в 1,5-2 раза от скорости газов на срезе КС до среза сопла, в такой же степени увеличиваются расстояния между арочными слоями в начале и в сопла на схеме.


Рис.32. Схема сил отталкивания в сферическом слое, создающие разгонную силу на молекулах перпендикулярно слою. (поперечный разрез)


Рис.33. Схема сил отталкивания в сферическом слое, создающие разгонную силу на молекулах перпендикулярно слою.Вид ссверху на шестигранную ячейку молекул в единичном сферическом изобарном слое в ЖРД.(Вид на силы в плоскости слоя в плане)

Именно эта перпендикулярная к изобарному слою сила и разгоняет слой в перпендикулярном к его поверхности направлении.

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

Торможение потоков в КС ЖРД


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

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

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

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

Попытаемся найти этим предположениям подтверждения в существующих учебниках по теории работы ЖРД.

Что-то в работе ЖРД мы смогли посчитать по известным нам законам газовой динамики и свести всё к вполне вменяемому энергетическому балансу по КС от дна до критического сечения.

Осталось только объяснить поддержание избыточного давления в КС непосредственно перед дырой наружу.

Что мешает струе газов разогнаться в пустом отверстии до запирающего давления Рv=Ркс, а не останавливаться на величине 62% от Ркс?

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

Бытовая логика не учитывает реальное взаимодействие по отдельным направлениям между отдельными молекулами газа по СТГ.

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

Именно это нам и нужно было доказать!

То есть в критическом сечении КС существует два разных давления :

  1. давление вдоль оси потока Ркр-о близкое к нулю ( равно давлению внешней среды)

  2. давление на стенки критического сечения поперёк потока Ркр-т= 18 атм (18% от максимального давления в КС, равного Ркс=100атм), которое фактически не регистрируемо инструментально из-за высокой кривизны и динамичного ускорения потока в критическом сечении КС. Указанное в расчётных формулах Ркр=0,574*Ркс (см.рис.11) относится не к критическому сечению КС, а к более глубокому сечению в конической части КС с меньшей в два раза скоростью газов, где ещё возможно замерять параметры на стенке, адекватные параметрам по всему сферическому слою данного сечения.

По расчётные формулам из учебника (см.рис.11).

Vа =(1,15* 8,31 * 3534 / 0,031)^0,5=1043 м/с

Полученная скорость звука Vа =1043м/с в целых 2 раза ниже, чем получаемая в реальном ЖРД РД-107 скорость газов согласно замеряемой силе тяги КС и соответствующей ей скорости при заданном постоянном массовом расходе G.

В любом случае остаётся вопросы :

  • Почему газы не могут разгонятся в сужающемся сопле выше некой фиксированной расчётной скорости, называемой "местной скоростью звука"?

  • Какую долю энергии мы срабатываем на достижение расчётной "скорости звука-MAX" в критическом сечении КС РД?

    В итоге пришлось сильно откорректировать расчётные значения.

    Так при принятой Vкр=0,5*Vмах-беск=2311м/с замеряемая температура в критическом сечении составляет Т=2880К (вместо Т=3543К) с поперечным давлением Рт=18 атм (вместо 57 атм), при начальных в КС параметарах Ткс=3800К и Ркс=100атм.

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

Откорректированный график должен быть менее гладким, имея резкие переходы на коноидальном заходе в КС, резкий скачёк на критическом сечении и излом при переходе в сопло Лаваля (см.рис.34). Причём такие резкие трансформации в зоне критического сечения КС будут у графиков всех параметров (V, T , P).

Конусная часть сопла заканчивается при срабатывании центрального конуса газов из КС. На графиках появились горизонтальные полки в зоне постоянного разрежения при срабатывании конуса газов из КС.

Интересно, что после окончания конусного участка КС начинается снижение угла наклона стенок к оси, что означает снижение темпа расширения газов в связи с падением скорости «звука» в слое газа из-за снижения его температур при расширении.

Рис.34. Откорректированный график параметров газов в тракте ЖРД. Для понимания добавлено несколько линий связи разреза ЖРД с графиком для обозначения важных границ отдельных функциональных частей ЖРД. Конусная часть сопла заканчиваяется при срабатывании центрального конуса газов из КС. На графиках появились горизонтальные полки в зоне постоянного разрежения при срабатывании конуса газов из КС.

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

Такое последовательное трёхмерное расширение газов по разным направлением как раз реализовано в типовом ЖРД.

Так в ЖРД сначала происходит расширение по оси КС вдоль скорости, с выработкой продольной части сил упругости поля отталкивания молекул.

Процесс продольного расширения газов за счёт внутренней энергии полностью завершается в критическом сечении КС при достижении "локальной скорости звука" .

Термин "локальная скорости звука" взят в кавычки, так как он является простой натяжкой расчётной формулы на реально полученное значение скорости потока в критическом сечении КС без нормального физичного обоснования данного явления.

Затем в расширяющейся части сопла Лаваля происходит плоское расширение слоёв в поперечных направлениях к скорости (вдоль сферического слоя молекул).

Схема такого неравновесного послойного расширения газов в отдельных частях ЖРД представлена на схеме (см. рис.30-33.)

Численная оценка сил и ускорений в арочных слоях в сопле Лаваля

Численная оценка сил и ускорений в арочных слоях в сопле Лаваля.

Расстояние между молекул в сферическом слое равно

дL=(Vm / Nа)^(1/3)

Где Vm- молярный объём в данных условиях, Nа-число Авогадро (количество молекул в 1 моль вещества).

Угол Ам в сферическом слое крайне мал, но при этом масса молекулы газа так же крайне мала.

Сила Fтс=6*Fс* sin(Ам)

Где множитель 6- это количество ближайших соседних молекул в сферическом слое к одной рассматриваемой молекуле, SinАм- отношения половины растояния между молекул в сферическом слое к радиусу кривизны слоя в этом месте, а сила Fc приближённо может быть расчитана как давление газа Ргаз отнесённое на элементарную грань шестигранной или квадратной призмы элементарного объёма молекулы.

Fc=Ргаз*dL^2=Ргаз*((Vm/Nа)^(1/3))^2=Ргаз*(Vm/Nа)^(2/3)

Величина угла Ам =sinАм при малых углах равна

sinАм=0,5*дL/Rсл=0,5*(Vm/Nа)^(1/3)/Rсл

Ускорение молекулы составит

а=F/m=6* Fс* sinАм/ (Mr/Na)

а= (6*Fс*0,5*(Vm/Nа)^(1/3)/Rсл)/ (Mr/Na)=(3*Fс*Na*(Vm/Nа)^(1/3)) / (Rсл*Mr)

Дробные степени суммируются при перемножении их одинаковых оснований , таким образом показатели степени 1/3 и степени 2/3 при элементарном объёме (Vm/Nа) складываются, становясь единицей:

а=(3*Ргаз*(Vm/Nа)^(2/3)*Na*(Vm/Nа)^(1/3)) / (Rсл*Mr)

а=(3*Ргаз*(Vm/Nа)*Na) / (Rсл*Mr)

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

а=3*Ргаз*(Vm/Nа)*Na / (RслMr)=3*Ргаз*Vm/ (Rсл*Mr)

Итого:

Ускорение молекулы в сферическом слое в сопле Лаваля равно:

а=3*Ргаз*Vm / (Rсл * Mr)

Так как произведение Ргаз*Vm=R*T, то ускорение можно выразить и через одну температуру:

а=3*Ргаз*Vm/(Rсл*Mr)=3*R*T/(Rсл*Mr)

а=3*R*T/(Rсл*Mr)

То есть динамика разгона зависит только от переменных температуры и кривизны слоя.

Кривизна слоя задаётся формой расширения в сопле Лаваля.

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

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

Уже поверхностный анализ показывает, что ускорение будет огромной величиной, так как в числителе присутствует большое число от произведения RT, а в знаменателе находятся малые величины радиуса и молярной массы ( для продуктов сгорания ЖРД Mr=0,031кг/м3), что ещё больше увеличивает конечное ускорение молекул в сферическом слое. Считаем , что величины входящие в формулу равны: Rсл=2*Dкр=0,4м (непосредственно за выходом из КС в сопло)

Т=3500К
Mr=0,031 кг/м3 (СО2+Н2О)

Тогда ускорение равно:

а=3*R*T/ (Rсл*Mr)=3*8,31*3500/ (0,4 * 0,031) =7 036 693 м/с2

При таком ускорении скорость 3400м/с будет достигнута за время:

дt=3400 / 7036693 = 0,00048с

При этом молекулами газа при таком ускорении будет пройдено расстояние:

L=3400* 0,00048 / 2= 0,816м

Цифры вполне согласуются с формой и размерами реальных ЖРД .

В тоже время в самой КС кривизна сферического слоя давит в обратную сторону, то есть внутрь КС, тем самым создавая избыточное статическое давление внутри КС.

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

Малое давление Рторм может существовать при больших избыточных силах Fc (давление поперёк потока), если сферический слой в критическом сечении КС имеет крайне малую кривизну ( Rсл уходит в бесконечность).

В конце коноидального захода на участке прямой трубы в критическом сечении прогиб арок молекулярных слоёв становится нулевым и остаточное давление газов вдоль оси потока мгновенно срабатывается в ускорение газов от Vкс-s3=741м/с до Vкр=2311м/с., мгновенно охлаждая газ на 870К.

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

В то же время в основной части КС скорость газов столь мала, что разгон газов в конической части КС происходит как у несжимаемой жидкости, а подсчёт скорости по сечениям КС производится исходя из постоянного массового расхода при постоянной плотности газа, то есть чисто геометрическими методами.

Резкое расширение газов в сопле Лаваля создаёт давление на стенках сопла, имеющих проекцию на ось ЖРД.

Работа газов поперёк сопла на убегающей от слоя стенке сопла- это и есть полезная работа сопла ЖРД в полёте.

Максимальный КПД ЖРД как движителя ракеты возникает при скорости ракеты равной скорости газов на срезе сопла.

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

Получается, что летящая в атмосфере на скорости 3000м/с ракета оставляет за собой неподвижный столб горячих газов, который даже не расширяется при равенстве давлений на срезе сопла с атмосферой в высотном ЖРД, или незначительно расширяется, если сопло было на первой ступени РН с расчётом на высокое атмосферное давление на уровне моря с достижением скорости 3000м/с в высоких и менее плотных слоях атмосферы.

После остывания этот столб газов за ракетой ещё и сильно сжимается.

Форма сопла Лаваля в ЖРД

Узкое коническое сопло Лаваля - это не единственная возможная форма сопла.

Сопло можно разворачивать на большие углы, в том числе до 180 градусов.

Примером вырожденного сопла может послужить «развёрнутое» или «плоское сопло».

Вспомним, что в вакууме газ расширяется из точки в сферу, увеличивая радиус почти в бесконечность.

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

Таким «идеально-теоретическим» соплом РД будет плоскость.

В одном полупространстве от этого плоского сопла начнёт расширятся облако горячего газа в виде идеального полушария, оказывая строго перпендикулярное давление на плоскость, там где сферический слой газа будет касаться плоскости «сопла». (см.рис.35)

Рис.35. Плоское дисковое сопло для РД. А. Схема переноса слоёв с узкого конического ориентированного сопла на плоское сопло. Б. Схема слоёв расширяющихся газов в дисковом сопле.

Тут возникает только вопрос с КПД такого плоского сопла: Будет ли оно выше, чем у обычного направленного сопла?

КПД плоского сопла будет однозначно ниже, чем у конического сопла Лаваля.

Выражаться это сравнительное КПД плоского и конического сопел будет достаточно простым соотношением: площади сечения шаровой поверхности (площадь круга) к площади полусферы.

Если вспомнить, что площадь сферы: Sсф = 4*R^2 * 3,14

Площадь полусферы: Sсф/2=2*R^2 * 3,14

А площадь круга: Sкр= R^2 * 3,14

То относится КПД плоского и направленного сопел будут как: Sкр : Sсф/2 = 1:2 =0,5

То есть КПД плоского сопла будет почти в 2 раза ниже, чем КПД узкого направленного сопла.

Это значит, что на создание одинаковой тяги потребуется расхода в 2 раза большего количества топлива при в 2 раза меньшем КПД.

Употребление слова «Почти» связано с неустранимой сферичностью потока на срезе сопла (незначительная выпуклость с радиусом сильно больше размера сопла по длине у криволинейных сопел), так что маленькое ухудшения превосходства КПД от 2-х кратности будет всегда.

Кстати, это двукратное усиление КПД при направленном сопле используют при подрыве взрывчатого вещества (далее ВВ) во взрывчатых боеприпасах некоторых типов.

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

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

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

Вернёмся к плоскому соплу.

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

Кстати, шаровой разлёт газа в вакууме с непрерывным приращением скорости отдельных молекул газа- это ещё одна демонстрация недостоверности КТГ.

Так в сильно разреженном шаровом слое молекулы никак не могут друг с другом сталкиваться для взаимного разгона, так как они уже летят в разные стороны на больших расстояниях друг от друга. Но при этом в шаровом слое молекулы продолжают друг от друга отталкиваться продолжая разгон при невероятно высоких степенях расширения (разрежения).

В рамках описанной модели ЖРД по СТГ становится понятным стремление к подъёму давления в КС при сохранении неизменной температуры (скорости струи) на критическом сечении КС. Так при большем давлении на одинаковой скорости потока можно пропустить тот же массовый расход G через сечение меньшей площади. То есть получается более компактная КС, или кратно большая мощность при неизменных размерах. При этом для расширения в сопле Лаваля остается больший диапазон разрежения до уровня внешнего давления подпора атмосферы.

Таким образом, повышение давления уменьшает габариты проточной части КС в ЖРД. Правда радужность картины снижается из- за повышения толщин стенок КС на более высокое давление. Также утяжеляется и становится более мощным турбонасос для пропихивания того же расхода G несжимаемого жидкого топлива в КС с более высоким давлением.

Давление в КС первых баллистических ракетах начиналось с 60 атм ( как в дизельных двигателях) . Теперь давление достигло 250атм (РД-170).

Так по нашему расчёту в критическом сечение поперечное давление составляет 18атм.

Для срабатывания избыточного давления 18 атм до уровня наружного атмосферного давления требуется сопло достаточно небольшого расширения от диаметра критического сечения КС.

У двигателя РД-107 (Ркс=60атм) степень расширения составляла всего 18,86. То есть диаметр увеличивался от критического сечения КС до среза сопла всего в 4,34 раза (см. фотото.рис.26)

Тогда для 250атм (РД-170) потребовалась большая кратность расширения в сопле, достигнув величины 36,8 раза, или рост диаметра сопла в 6,07 раз (рис.36).

Подняв давление в 4 раза от 60атм до 250 атм можно снять в 4 раза большую тягу с КС и Сопла практически неизменного размера.

Интересно заметить, что повышение давления в КС приводит к уменьшению длины КС, так как одинаковое время протекания процесса испарения и горения топлива при большей плотности газов в КС обеспечивается на более коротком пути от форсунок до критического сечения КС (см рис.26,36). Такое сокращение длины КС позитивно сказывается на общем снижении веса ЖРД.

Выше 250 атм поднимать давление уже не выгодно, так как начинается столь значительный рост массы двигателя для обеспечения прочности стенок, что прирост массы уже не компенсируется ростом импульса тяги ЖРД (РД-107 импульс в вакууме 313с при Р=60атм, РД-171 импульс в вакууме 337 с при Р=250атм)

Рис.36. Фото,разрез и чертёж одной камеры РД-180 с рабочим давлением 250атм в КС. Хорошо заметно резкое заужение в критическом сечении. КС выглядит явно "маленькой" по сравнению с большим по размеру соплом с высокой кратностью расширения.

Интересно отметить, что рост мощности отдельной КС ЖРД остановился на уровне тяги 200 тс при давлении в Ркс=250атм ещё 50 лет назад на двигателях типа РД-170.

Это можно считать пределом возможности геометрического увеличения размера КС ЖРД, что связано с достижением критически большого размера отверстия критического сечения КС.

Так РД-107 (тяга 102тс) отличается по тяге от РД-170 (тяга 806тс) в 8 раз, но из-за большего в 4 раза давления в КС (60атм против 250 атм) критическое сечение КС двигателей отличается всего в 2 раза, а диаметры всего в 1,4 раза (корень из 2-х).

При этом в росте линейного размера критического сечения КС ограничивающую роль играют роль два фактора :

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

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

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

    Так по оси КС через критическое сечение проскакивают потоки плотных холодных низкоскоростных газов несгоревшего топлива, резко роняя давление в КС позади себя. Именно эти проскоки недогоревшего топлива через критическое сечение КС и являются основным источником возникновения скачков давления в КС больших ЖРД.

Это наводит на размышления о работоспобности двигателя F-1 c РН Сатурн-5, участвовавшего в лунной программе США.

Так КС F-1 имеет столь малое заужение в критическом сечении (см рис.37-38), что не позволяет надеяться на достижение заявляемого давления 70 атм (сравните с сечением РД-107 с чуть меньшим давлением Ркс=60атм). Скорее такой ЖРД будет работать как довольно мощная факельная горелка, фонтанирующая плохо сгоревшим топливом, что как раз и видно на кадрах его испытания на стенде (см.рис.39).

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

Хорошо известны весьма скептические слова конструктора советских ЖРД Валентина Глушко, когда он узнал про разработку в США однокамерного ЖРД на 800тс тяги: "Пусть попробуют... упрутся в ту же стену, в которую упёрлись мы".

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

Самый мощный на текущий момент ЖРД с тягой 806 тс является именно четырёхкамерный ЖРД Глушко РД-170, а у него тяга одной КС составляла всего 200тс (при Ркс=250 атм), которые и можно считать геометрическим пределом. При пересчёте отношений мощностей и давлений получается, что у F-1 диаметр критического сечения приблизительно в 3,6 раза больше, чем у одной камеры самого мощного РД-170, или в 5 раз больше чем у РД-107 с близким к F-1 давлением в КС.

Разница диаметров критических сечений у РД-170 и РД-107 составляет всего 1,37 раза или +37% от меньшего к большему.

Предположительную схему прорыва холодных плотных газов из КС при чрезмерно широком критическом сечении представлена на схеме (см.рис.41-43)

При этом на малых ЖРД с критическим сечением небольшого диаметра легко достигаются высокие давления в КС, да и работают такие мелкие ЖРД весьма стабильно.

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

Именно к таким поликамерным односопловым ЖРД относятся ЖРД с центральным телом: сопло с центральным диском, сопло с центральным круглым телом (штыревые) или с плоским центральным телом (плоское сопло) (См .рис.21)

Критика от читателей

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

Про твёрдотопливные РД с тягой 1500тс- это хорошее замечание!

Но на него есть ответ, который только добавляет ценности данной статье:

В ракетных ускорителях с твёрдотопливными реактивными двигателями (далее ТРД) присутствует КС огромного размера, в которой размещена и горит вся одновременно огромная топливная шашка. Остановить горение в топливной шашке не в состоянии никакими взрывами никакие пульсации в сопле.

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

Потому данное ограничение с тягой около 200тс для КС с керосин-кислородом не действует на твердотопливные РД.

Возможно, что изменение формы КС (конусная или шаровая) для керосиновых ЖРД и увеличение объёма КС (чтобы горело дольше) позволят сделать более мощные ЖРД на керосин-кислороде, чем текущие 200-220тс на одну камеру.

Рис.37. Разрез камеры ЖРД F-1 от ракеты Сатурн-5. Виден плавный переход с малым заужением от КС к критическому сечению. Также видна малая степень расширения в сопле Лаваля, что говорит о достаточно низком остаточном давлении в струе газов после критического сечения КС.

Рис.38.А. Вид в КС со  стороны сопла.
Рис.38.А. Вид в КС со стороны сопла.
Рис.38.Б. Вид на КС отработанного пуском  двигателя F-1 с оторванным при падении в океан соплом.
Рис.38.Б. Вид на КС отработанного пуском двигателя F-1 с оторванным при падении в океан соплом.

Коментарий к Рис.38. На фото внутренней стороны сопла ЖРД F1 от РН Сатурн-5 виден малый перепад сечения КС на заужение критического сечения. Также на фото поломанного и поднятого с морского дна РД видны сравнительно малые размеры КС по сравнению с большим размером критического сечения КС.

Рис.39. Фото запуска РД F-1 на наземном испытательном стенде.
Рис.39. Фото запуска РД F-1 на наземном испытательном стенде.

Комментарий к рис.39. испытания ЖРД F-1 на стенде. Хорошо видно тёмный поток несгоревшего топлива по краю потока газов из сопла - это выхлоп газогенратора от турбонасоса, который подается внутрь сопла из толстого "рогалика"-коллектора, видного на фото. Несгоревшее газообразное топливо резко воспламеняется от контакта с кислородом воздуха далеко за краем сопла. В глубине струи сквозь внешний дым виден светящийся поток высокоскоростной струи газов из ЖРД, который слишком тонок для данного сопла и не заполняет даже верхню часть сопла, предназначенную специально для эффективной работы в плотных слоях атмосферы на старте.

Рис.40. Момент старта РН Союз с ЖРД типа РД-107.
Рис.40. Момент старта РН Союз с ЖРД типа РД-107.

Коментарий к Рис.40. На фото запуска РН Союз с двигателями РД-107 видно, что пламя из сопла на столько прозрачное, что сквозь струю газов из сопла видна дальняя кромка этого сопла. При этом слой холодного газа из несгоревших паров керосина-охладителя тоже есть, именно это топливо на периферии струи вспыхивает ярким непрозрачным пламенем в кислороде воздуха уже далеко после вылета за края сопла.

Рис.41. Различные варианты масштабирования ЖРД. А- исходный РД малой мощности. Б-линейное масштабирование всех размеров  разгонной части при неизменной длине зоны горения КС. В- искажение  формы КС  и динамики потоков газов в зоне сужения КС при сохранении длины конусной части.
Рис.41. Различные варианты масштабирования ЖРД. А- исходный РД малой мощности. Б-линейное масштабирование всех размеров разгонной части при неизменной длине зоны горения КС. В- искажение формы КС и динамики потоков газов в зоне сужения КС при сохранении длины конусной части.

Комментарий к Рис.41. Схема возникновения прорыва плотных газов из КС через критическое сечение при некорректном масштабировании геометрии КС при повышении мощности ЖРД:

А- стабильно работающий ЖРД с КС малой мощности с длинной переходной часть от КС до критического сечения;

Б- стабильно работающий ЖРД с увеличением диаметра КС в 2 раза ( мощности в 4 раза) и с пропорциональным масштабированием переходной часть от КС до критического сечения. При этом цилиндрическая зона горения топлива имеет постоянную длину по физико-химическим параметрам и не подлежит масштабированию по длине при увеличении мощности ЖРД при увеличении сечения КС;

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

Рис.42. Сравнение двигателей РД-107 с F-1 при пропорциональном масштабировании:
А- вид одного  сопла двигателя РД-107;
Б- Отмасштабированый в 5 раз по линейному  размеру РД-107 и в 25 раз по мощности, при этом зона  горения КС не масштабировалась по длине;
В- реконструкция двигателя F-1 по имеющимся данным по  заявляемой степени расширения сопла 16 и имеющимся  фото.
Рис.42. Сравнение двигателей РД-107 с F-1 при пропорциональном масштабировании: А- вид одного сопла двигателя РД-107; Б- Отмасштабированый в 5 раз по линейному размеру РД-107 и в 25 раз по мощности, при этом зона горения КС не масштабировалась по длине; В- реконструкция двигателя F-1 по имеющимся данным по заявляемой степени расширения сопла 16 и имеющимся фото.
Рис.43. Сравнение занимаемых пространств на РН двигателями различной компоновки:
А- четырёх камерный РД-107 (тяга 102тс) и пакет из 7 разгонных блоков с РД-107 на общую тягу 714тс в сравнении с кругом сопла одного F-1.
Б- Реконструкция двигателя F-1 по доступным данным. 
В- четырёх камерный РД-170 (тяга 806тс) отмасштабированый из КС РД-107 с кратностью расширения сопла 36,7.
Прямоугольнички в критических сечениях и в соплах позволяют визуально оценить кратность расширения внутри  самих ЖРД и кратность масштабирования РД-107 в F-1.
Рис.43. Сравнение занимаемых пространств на РН двигателями различной компоновки: А- четырёх камерный РД-107 (тяга 102тс) и пакет из 7 разгонных блоков с РД-107 на общую тягу 714тс в сравнении с кругом сопла одного F-1. Б- Реконструкция двигателя F-1 по доступным данным. В- четырёх камерный РД-170 (тяга 806тс) отмасштабированый из КС РД-107 с кратностью расширения сопла 36,7. Прямоугольнички в критических сечениях и в соплах позволяют визуально оценить кратность расширения внутри самих ЖРД и кратность масштабирования РД-107 в F-1.

На рис.42 отчётливо видно, что КС F-1 сильно меньше по объёму, чем отмасштабированная из РД-107 версия на равную мощность. Так площадь форсуночного дна у F-1 приблизительно в 3 раза меньше, чем у масштабной версии РД-107. Такое удивительное компактное строение КС позволяет предполагать невозможность полного прогорания топлива в зауженой КС F-1, что должно приводить к выбросу факела с недогоревшим топливом и пониженное давление в КС.

На рисунке (рис.42.Б) в левой части представлена возможная версия формы отмасштабированной КС со сферической формой форсуночного дна и строго конической КС до захода на критическое сечение. Такая форма КС позволяет обеспечить устойчивое структурирование потоков газа до критического сечения, тем самым избавляясь от неравномерности скоростей и неустойчивой динамики потоков газов внутри КС.

На рисунках (рис.43) даны сравнения габаритов двигателей на равную тягу в разном исполнении. Из чего можно сделать заключение о целесообразности борьбы в 1950-х за укрупнение двигателей на низком давлении 60-70 атм в КС. Так на рисунках (рис.43.А-Б) видно, что односопловой F-1 занимал сильно меньше места в плане РН, чем пакет из 7 четырёхкамерных двигателей типа РД-107 на приблизительно ту же развиваемую тягу.

Но сравнение соседних справа картинок (рис.43.Б-В) даёт результат уже не в пользу F-1, так как РД-170 при большей тяге ещё и значительно компактнее при давлении 250 атм в КС. Так же РД-170 (РД-171) всё ещё серийно производится в РФ, причём в разных версиях по числу камер и тяге (двухкамерный РД-181, однокамернный РД-191). А вот производства двигателей F-1 сейчас не существует и покрыто туманом легенд о прошлых победах "Великой Америки".

Ступенчатые процессы расширения в потоке газов из ЖРД с центральным телом

 

(публикация от 31-12-2022)

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

На основании модели ЖРД с применением СТГ возможно объяснить не только происходящее в сопле Лаваля, но и понять происходящее в альтернативных моделях ЖРД с соплами внешнего расширения: штыревые, плоские и тарельчатые.

Эти типы ЖРД отличаются наличием узкого и длинного критического сечений КС, создающего относительно тонкую плоскую струю газов на выходе из КС.

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

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

При рассмотрении сечения кольцевой КС с плоской струёй сечение будет выглядеть очень похоже на сечение круглой КС с соплом Лаваля.

Для сопла Лаваля мы рассматривали половинку разреза, предполагая зеркальное отображение картины относительно оси симметрии. Ну, а  для плоского потока ЖРД с кольцевой КС у похожего разреза вместо  центральной оси появится реальная стенка внешнего сопла (см.рис.44).

Рис.44. Сравнение процессов расширения газов, происходящих в ЖРД с различными конфигурациями критических сечений КС и разными типами внешних сопел: А- расширение газов в ЖРД с соплом Лаваля, Б- расширение газов в ЖРД с соплом с центральным телом (штыревым или тарельчатым).
Рис.44. Сравнение процессов расширения газов, происходящих в ЖРД с различными конфигурациями критических сечений КС и разными типами внешних сопел: А- расширение газов в ЖРД с соплом Лаваля, Б- расширение газов в ЖРД с соплом с центральным телом (штыревым или тарельчатым).

 

Рис.44. Сравнение процессов расширения газов, происходящих в ЖРД с различными конфигурациями критических сечений КС и разными типами внешних сопел: А- расширение газов в ЖРД с соплом Лаваля, Б- расширение газов в ЖРД с соплом с центральным телом (штыревым или тарельчатым).

В ЖРД с центральным телом направление потока газов из критического сечение наклонено вбок относительно оси ЖРД и направления полёта ракеты. Угол наклона может изменяться в широких пределах, достигая 90 градусов в отдельных случаях.

Плоская струя газов из щелевого критического сечения вырывается параллельно участку плоскости наружного сопла, при  этом имея возможность свободно расширяться в пустоту позади ЖРД.

По аналогии с моделью работы традиционного ЖРД с соплом Лаваля можно сказать, что струя газов после критического сечения должна начать  резко расширяться перпендикулярно направлению полёта струи, так как продольное расширение газа уже  завершилось в критическом сечении.

При этом в ЖРД с центральным телом пространство для одностороннего расширение в  бок оказывается ничем не ограничено, что позволяет поверхностным слоям газа на боках струи практически мгновенно срабатывать поперечное давление струи в разгон молекул до максимально возможной скорости, равной остатку энергии давления в этом направлении (см.рис.47-48).

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

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

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

Ширина тарельчатого сопла связана с шириной критического сечения из отношения Vкр/Vа.

То есть длина полного срабатывания плоской струи в поперечное расширение составит длину не более:

 В=Нкр* Vкр/Vа.

Где В- ширина конусного участка сопла, Нкр- толщина щели кольцевого критического сечения КС.

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

Итоговой скоростью струи на таком плоско-тарельчатом сопле составит:

Vcоп1=Vкр*(2^0,5)= Vкр*1,41.

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

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

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

Рис.45. Схема сил в плоско-цилиндрическом слое газов в неравномерно структурированном по трём осям потоке газа из ЖРД с центральным телом.
Рис.45. Схема сил в плоско-цилиндрическом слое газов в неравномерно структурированном по трём осям потоке газа из ЖРД с центральным телом.

Рис.45. Схема сил в плоско-цилиндрическом слое газов в неравномерно структурированном по трём осям потоке газа из ЖРД с центральным телом.

Рис.46. Схема пространственной  неоднородности плотности  распределения  молекул газа в плоском потоке, истекающим из  ЖРД с центральным телом: А-1- штыревое круглое сопло  в поперечном сечении, А-2- штыревое круглое сопло  в продольном сечении, Б-1- поперечное сечение широкой юбки внешнего сопла, Б-2- продольное сечение широкой юбки внешнего сопла
Рис.46. Схема пространственной  неоднородности плотности  распределения  молекул газа в плоском потоке, истекающим из  ЖРД с центральным телом: А-1- штыревое круглое сопло  в поперечном сечении, А-2- штыревое круглое сопло  в продольном сечении, Б-1- поперечное сечение широкой юбки внешнего сопла, Б-2- продольное сечение широкой юбки внешнего сопла

Рис.46. Схема пространственной  неоднородности плотности  распределения  молекул газа в плоском потоке, истекающим из  ЖРД с центральным телом: А-1- штыревое круглое сопло  в поперечном сечении, А-2- штыревое круглое сопло  в продольном сечении, Б-1- поперечное сечение широкой юбки внешнего сопла, Б-2- продольное сечение широкой юбки внешнего сопла.

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

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

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

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

 

Рис. 47. Разгон газов в ЖРД с плоским дисковым центральном телом  в критическом сечении КС:
 А- режим расширения газов в вакууме, Б- режим расширения газов с подпором атмосферы.
Рис. 47. Разгон газов в ЖРД с плоским дисковым центральном телом в критическом сечении КС: А- режим расширения газов в вакууме, Б- режим расширения газов с подпором атмосферы.

Рис. 47. Разгон газов в ЖРД с плоским дисковым центральном телом  в критическом сечении КС: А- режим расширения газов в вакууме, Б- режим расширения газов с подпором атмосферы.

Рис. 48. Разгон газов в ЖРД со штыревым центральном телом  в критическом сечении КС:  
              А- режим расширения газов в вакууме, Б- режим расширения газов с подпором атмосферы.
Рис. 48. Разгон газов в ЖРД со штыревым центральном телом в критическом сечении КС: А- режим расширения газов в вакууме, Б- режим расширения газов с подпором атмосферы.

 

Рис. 48. Разгон газов в ЖРД со штыревым центральном телом  в критическом сечении КС:     А- режим расширения газов в вакууме, Б- режим расширения газов с подпором атмосферы.

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

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

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

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

Динамика разгона молекул газа при прямом разгоне  в направлении градиента от  давления в КС к вакууму на несколько порядков больше, чем динамика разгона перпендикулярно напряжению плоского слоя при малой кривизне самого слоя. Эта разница в динамике разгонов и позволяет  кривизне центрального штыря обгонять кривизну от свободного расширения поперёк потока.

Схожий механизм создания дополнительной тяги на криволинейной юбке ЖРД с центральным диском от недорасширившихся газов с большой остаточной радиальной скоростью в плотных слоях атмосферы.

 

Конечная  предельная скорость истечения при сложении трёх равных ортогональных векторов может составить величину:

Vcоп2=Vкр*3^0,5= Vкр*1,73

При этом предельная  величина Vcоп2=Vкр*1,73 может быть достигнута только при бесконечно большом коэффициенте расширения юбки сопла и в полном вакууме, то есть практически никогда в  реальности.

В тоже время скорость Vcоп1= Vкр*1,41  достигается уже почти полностью на коротком участке тарельчатого сопла с углом около 45 градусов к оси и в условиях атмосферного давления.

Получается, что предельно короткий ЖРД  с центральным диском обеспечивает на крайне компактном сопле тягу на уровне  близком к  1,41/1,71 =0,82 от теоретического максимума или 82 %  от  ЖРД с полноразмерным бесконечным соплом Лаваля.

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

Геометрия односоплового ЖРД неограниченной тяги с линейно-щелевым критическим сечением

(дополнение от 13-01-2023)

Исходя из рассмотренных конфигураций ЖРД с центральным телом различного вида можно сделать предположение, что вполне реально создать односопловой ЖРД неограниченной мощности с единой неразделённой КС.

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

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

В щелевом критическом сечении наименьший  поперечный размер щели будет достаточно мал, то есть много меньше достигнутого предела диаметра критического сечения (как в ЖРД РД-170).

Такая узкая щель из КС будет  обеспечивать формирование стабильных и надёжных тормозных арочно-цилиндрических слоёв высокой кривизны перед линейным критическим сечением.

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

Идея сама по себе не нова, и фигурирует в общедоступных учебниках по проектированию ЖРД [2]. Фрагмент страницы учебника из интернета с описанием такого односоплового мегаЖРД представлена ниже (см.рис.49).

Рис. 49. Страница учебника с описанием ЖРД с   кольцевым  щелевым критическим сечением в составе тороидальной КС.
Рис. 49. Страница учебника с описанием ЖРД с   кольцевым  щелевым критическим сечением в составе тороидальной КС.

Рис. 49. Страница учебника с описанием ЖРД с   кольцевым  щелевым критическим сечением в составе тороидальной КС. На рисунке В)  изображён оптимизированный щелевой мегаЖРД  с тороидальной КС, который практически  аналогичен в разрезе с конфигурацией ЖРД с центральным диском на рисунке Рис.47.

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

 На данный момент уже есть действующий четырёх камерный ЖРД РД-171 с тягой 800тс, при этом ещё и выпускаемый серийно. Четырёх камерный РД-171 по габаритам также более чем в 2 раза короче своего однокамерного соперника F-1 c равной тягой (см.рис.43).

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

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

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

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

В этих экспериментах экономичность вообще мало учитывается, а по стоимости одно устройство  может превосходит стоимость крупных инфраструктурных проектов типа «Крымский мост».

В случае строительства кольцевого мега ЖРД предметом исследования будет не только сам факт создания работоспособного сверх мощного однокамерного  ЖРД, но и возможность изучения на стенде взаимного влияния протяжённых объёмов КС на стабильность  работы ЖРД и формирование устойчивых колебательных систем внутри самого ЖРД.

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

Так в   советской лунной ракете Н-1  на первой ступени было установлено аж 30 шт. ЖРД НК-33 (тяга по 154-171тс). И эта система из трёх десятков ЖРД создавала такой вибрационный фон, что при резонансном наложении вибрационных волн друг на друга отдельные двигатели отключались по превышению давления в топливопроводах от гидроударов.

В тоже время на стабильно летающей ракете «Союз» было всего 5 шт. ЖРД РД-107 (тяга по 102тс), при этом каждый из четырёхкамерных двигателей стоял на отдельном ускорителе, которые в свою очередь объединялись через внешние узлы соединения в единую пакетную 1-ю ступень.

В более поздней мега РН «Энергия» (стартовая масса 2400т) было всего 8 шт. ЖРД, а именно: 4 шт. РД-170 (тяга по 806тс) на паре керосин-кислород и 4шт ЖРД РД-0120 (тяга по 200тс) на паре водород-кислород на центральном блоке РН.

Сборная пакетная конструкция РН «Энергия» имела четыре разгонных блока (пристёгнутых сбоков к центральному блоку), на которых монтировалось всего по одному РД-170. И только на центральном разгонном блоке ставилось параллельно четыре РД-0120 на одну опорную площадку.

Неудачную схему с 33-мя мелкими ЖРД Мерлин-2 (тяга по 200тс на паре метан-кислород) на одной РН сейчас пытаются воплотить американцы на свое сверхтяжелой РН «Старшип» (5тыс.т стартовый вес), а именно на первой ступени «Суперхеви». Но даже современные электронные системы управления не очень хорошо справляются с проблемой взаимного вибрационного влияния отдельных ЖРД в столь сложной многодвигательной системе гигантской РН.

Таким образом при создании стендового макета сверхмощного односоплового щелевого ЖРД было бы возможно определить взаимное влияние работы раздельных КС различного размера и формы, а так же выяснить исходящие импульсно-вибрационные нагрузки от двигателя на РН.

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

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

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

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

Экспериментальная отработка на ЖРД малой размерности эффективности работы КС с щелевым критическим сечением

 (Публикация от 04 февраля 2023г.)

Прежде чем замахиваться на мегаЖРД имеет смысл провести натурный низкозатратный эксперимент на элементной базе уже выпускаемых серийных ЖРД.

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

Что именно нового даст щелевое критическое сечение для работы ЖРД, если КС и сопло останутся малоизменёнными?

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

Вибрация  РН вызывается пульсацией давления в КС ЖРД с такой же частотой.

Именно эта частота пульсации давления в КС создаёт циклическое изменение сечения газовой струи за пределами сопла ЖРД, тем самым являясь мощным источником звуковых колебаний в атмосфере с частотой 5-20Гц. Именно газовая струя из сопла ЖРД своей пульсацией с частотой 5-20Гц и создаёт тот непередаваемый мощный рокот при старте ракеты с космодрома.

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

То есть, для укрупнения ЖРД имеется только один разумный аргумент, а именно:

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

Ту возникают резонный вопрос:

Позволяет ли применение щелевого критического сечения в ЖРД как-то снизить сами вибрации и пульсации давления в КС?

Или кольцевой мега-ЖРД способен обеспечить только механическое замыкание всех источников пульсаций на себя, но не будет их снижать как таковые?

На такие вопросы достаточно просто организовать экспериментальный ответ с низкими затратами.

Так в РФ есть хорошо отработанный за десятилетия выпуска и применения ЖРД РД-107 от РН «Союз». Его вибрационно-акустические характеристики также хорошо изучены. Так же у РД-107 есть рулевые камеры малого размера, которые с низкими затратами можно использовать для экспериментальных переделок (см. фото на рис.26, где рулевая малая КС в нижнем левом углу фото).

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

Для низкозатратного независимого эксперимента в одной из рулевых КС достаточно изменить форму критического сечения с круглого на овальное или прямоугольное с закруглёнными углами, при этом сохранить площадь критического сечения (см.рис.50).

 

Рис.50. Сравнение конфигурации КС и сопловой части ЖРД с круглым и щелевым сечением равной площади: А-исходная конфигурация осесимметричной  камеры ЖРД традиционного исполнения; Б- сравнение отдельных сечений при наложени друг на руга; В- сечение камеры щелевого ЖРД по центральной оси поперёк щелевого критического сечения; Г- сечение камеры щелевого ЖРД по центральной оси вдоль щелевого критического сечения.
Рис.50. Сравнение конфигурации КС и сопловой части ЖРД с круглым и щелевым сечением равной площади: А-исходная конфигурация осесимметричной  камеры ЖРД традиционного исполнения; Б- сравнение отдельных сечений при наложени друг на руга; В- сечение камеры щелевого ЖРД по центральной оси поперёк щелевого критического сечения; Г- сечение камеры щелевого ЖРД по центральной оси вдоль щелевого критического сечения.

Рис.50. Сравнение конфигурации КС и сопловой части ЖРД с круглым и щелевым сечениями равных площадей: А-исходная конфигурация осесимметричной  камеры ЖРД традиционного исполнения; Б- сравнение отдельных сечений при наложении друг на руга; В- сечение камеры щелевого ЖРД по центральной оси поперёк щелевого критического сечения; Г- сечение камеры щелевого ЖРД по центральной оси вдоль щелевого критического сечения.

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

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

Перспективы применения инновационных ЖРД с линейно -щелевым критическим сечением в существующих ракетно-космических программах Роскосмоса

После получения успешных сдвигов характеристике  в эксперименте с малой камерой от РД-107 будет возможно ставить вопрос об полномасштабном эксперименте над уже весьма дорогостоящим двигателем РД-191 с тягой 212 тс при давлении 254атм в КС, являющимся самым мощным односопловым  двигателем на топливе кислород-керосин с исключительно высоким диапазоном дросселирования 27-105% от номинальной тяги.

При анализе современных программ развития ракетно-космической отрасли РФ можно отметить, что вполне очевидна тенденция по замене 4-х соплового ЖРД РД-107 на более мощный односопловой РД-191 на новых РН, таких как:  лёгкой «Союз-2.1в» (масса 160 т) и универсальных разгонных модулях  (УРМ «Ангара» с 1хРД-191) для новых перспективных  лёгких РН «Ангара-1.2» (масса 171 т).

Также на базе УРМ «Ангара» собираются пакетные РН больших грузоподъёмностей:  средняя «Ангара-А3» (3 УРМ, стартовая масса 470т) и тяжёлая «Ангара-А5»(5 УРМ, стартовая масса 800т).  

То есть в привычной пакетной схеме типа РН «Союз» поменяли  4-х камерные РД-107 ( с давлением 60атм и тягой 102тс) на однокамерные  РД-191 (с давлением 254атм и тягой 212тс), а новую линейку более мощных   пакетно-модульных РН назвали «Ангара».

Если удастся в щелевой версии РД-191 поднять тягу за счёт роста давления, или снизить амплитуду  пульсации давления  в КС без изменения тяги при номинальном давлении в КС, то  это буде означать возможность пробития габаритного предела работоспособной КС, причём практически без ограничения по предельной мощности.

На базе нового форсированного РД можно будет создать новую линейку модульных РН «Ангара-Супер» причём с минимальными затратами на переделку конструктива самих УРМ.

Рис.51. Фото общего вида и схема работы односоплового ЖРД РД-191 (тяга 212тс при Ркс=254атм) с одиночным соплом от РД-171. На схеме видно, что криогенный окислитель полностью газифицируется  при выработке генераторного газа для  раскрутки турбонасоса, при этом остаток  жидкого топлива попадает в  КС пройдя по стенкам сопла и КС в качестве охлаждающей жидкости
Рис.51. Фото общего вида и схема работы односоплового ЖРД РД-191 (тяга 212тс при Ркс=254атм) с одиночным соплом от РД-171. На схеме видно, что криогенный окислитель полностью газифицируется  при выработке генераторного газа для  раскрутки турбонасоса, при этом остаток  жидкого топлива попадает в  КС пройдя по стенкам сопла и КС в качестве охлаждающей жидкости

Рис.51. Фото общего вида и схема работы односоплового ЖРД РД-191 (тяга 212тс при Ркс=254атм) с одиночным соплом от РД-171. На схеме видно, что криогенный окислитель полностью газифицируется  при выработке генераторного газа для  раскрутки турбонасоса, при этом остаток  жидкого топлива попадает в  КС пройдя по стенкам сопла и КС в качестве охлаждающей жидкости.

Рис.52. Разрез односоплового ЖРД РД-191 (макет). На разрезе видно, что разница Rкс и Rкр невелика (Rкс/Rкр < 1,5), то есть возможность уплощения критического сечения в  границах существующей геометрии КС сильно меньше, чем у КС РД-107 (где Rкс/Rкр >2,5).Тем не менее, уменьшение критического размера  по ширине щели будет приблизительно 1,7 раза (см.рис.53), что должно значительно снизить амплитуду пульсации давления в КС и повысить их частоту
Рис.52. Разрез односоплового ЖРД РД-191 (макет). На разрезе видно, что разница Rкс и Rкр невелика (Rкс/Rкр < 1,5), то есть возможность уплощения критического сечения в  границах существующей геометрии КС сильно меньше, чем у КС РД-107 (где Rкс/Rкр >2,5).Тем не менее, уменьшение критического размера  по ширине щели будет приблизительно 1,7 раза (см.рис.53), что должно значительно снизить амплитуду пульсации давления в КС и повысить их частоту

Рис.52. Разрез односоплового ЖРД РД-191 (макет). На разрезе видно, что разница Rкс и Rкр невелика (Rкс/Rкр < 1,5), то есть возможность уплощения критического сечения в  границах существующей геометрии КС сильно меньше, чем у КС РД-107 (где Rкс/Rкр >2,5).Тем не менее, уменьшение критического размера  по ширине щели будет приблизительно 1,7 раза (см.рис.53), что должно значительно снизить амплитуду пульсации давления в КС и повысить их частоту.

 

Рис.53. Сравнение конфигурации КС и сопловой части ЖРД РД-191 с круглым и щелевым сечением равной площади: А-исходная конфигурация осесимметричной  камеры ЖРД традиционного исполнения; Б- сравнение отдельных сечений при наложени друг на руга; В- сечение камеры щелевого ЖРД по центральной оси поперёк щелевого критического сечения; Г- сечение камеры щелевого ЖРД по центральной оси вдоль щелевого критического сечения, уменьшение критического размера  по ширине щели будет приблизительно в 1,7 раза от исходного Dкр
Рис.53. Сравнение конфигурации КС и сопловой части ЖРД РД-191 с круглым и щелевым сечением равной площади: А-исходная конфигурация осесимметричной  камеры ЖРД традиционного исполнения; Б- сравнение отдельных сечений при наложени друг на руга; В- сечение камеры щелевого ЖРД по центральной оси поперёк щелевого критического сечения; Г- сечение камеры щелевого ЖРД по центральной оси вдоль щелевого критического сечения, уменьшение критического размера  по ширине щели будет приблизительно в 1,7 раза от исходного Dкр

Рис.53. Сравнение конфигурации КС и сопловой части ЖРД РД-191 с круглым и щелевым сечением равной площади: А-исходная конфигурация осесимметричной  камеры ЖРД традиционного исполнения; Б- сравнение отдельных сечений при наложени друг на руга; В- сечение камеры щелевого ЖРД по центральной оси поперёк щелевого критического сечения; Г- сечение камеры щелевого ЖРД по центральной оси вдоль щелевого критического сечения, уменьшение критического размера  по ширине щели будет приблизительно в 1,7 раза от исходного Dкр.

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

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

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

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

 

Лабораторный эксперимент по проверке тезиса о «звуковом» пределе скорости в критическом сечении «пневматического реактивного двигателя» (ПРД).

 Предлагается поставить лабораторный эксперимент для подтверждения высказанных предположения о более высокой скорости в критическом сечении КС ЖРД, чем «скорость звука» при данной температуре газов.

Эксперимент будет проводиться в нормальных лабораторных условиях ( Т= +22С, Р=1 атм.)

В качестве пневматического реактивного двигателя (ПРД) используется пневматическай установка состоящая из следующих элементов:

- пневматически бак мембранный 24л с предварительным давлением 1,5атм (рессивер),

- клапан для заполнения рессивера  сжатым воздухом от компрессора до 8атм,

- манометр для контроля давления в рессивере,

- запорные шаровые кран (1/2”) на линии заполнения,

- запорные шаровые кран (3/4”) на линии выпуска воздуха из рессивера,

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

Пневматический реактивный двигатель (ПРД)
Пневматический реактивный двигатель (ПРД)

Рис.54. Общий вид лабораторного  «пневматического реактивного двигателя» (ПРД).

Определение тяги ПРД производится с помощью электронных настольных электронных весов с точностью измерения 1гр и цифровой индикацией.

Время засекается секундомером.

Эксперимент по определению тяги ПРД  имеет следующие этапы:

1.       Наполнение рессивера сжатым воздухом до заданного избыточного давления (4-8 атм)

2.       Установка ПРД над весами соплом перпендикулярно весам.

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

 

Лабораторный стенд с ПРД
Лабораторный стенд с ПРД

 Рис.55. Лабораторная установка для измерения тяги ПРД.

Результаты эксперимента №1.

 При предварительном наполнении ПРД до 4атм и конечном давлении 1,5атм (отсечка по мембране бака) весы показывают от 309гр в начале  до 190гр в конце, время истечения около 10с, при диаметре воздушного отверсти D=2,5мм.

Оценка результатов эксперимента №1.

Зная полученую тягу, размер отверстия и давление в рессивере мы можем рассчитать параметры воздушной струи.

Так предположим, что струя имеет скорость Vcт=Vа=330 м/с.

Тогда мы можем рассчитать плотность струи qст  из формулы импульса силы:

F=q*S*V^2

где q –плотность струи в критическом сечении отверсти, S- площадь сечения отверстия, V- скорость воздуха в в критическом сечении отверстия.

 Откуда:

q=F/(S*V^2)

q=3,03/(0,0000049*330^2)=5,67 кг/м3

Учитывая, что плотность воздуха 1,2кг/м3 при нормальных условиях, то при избыточном давлении 4 атм плотность воздуха в ресивере составит:

q.рес=1,2* (4+1)=6 кг/м3

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

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

При адиабатическом разгоне газа преобретается энергия

Езв=(Mr*Va^2)/2= (0,029* 330^2)/2= 1579 дж/моль

То есть именно такую энергию должен потерять газ при расширении от состояния №1:

Р1=(4+1)=5атм, Т1= 273+22=295К, V1= 22,4/5=4,48л

До состояния №2:

Р2, Т2, V2

При это выполняется закон адиабатического расширения:

Р1/Р2= (V2/V1)^1,4

Где 1,4- это показатель адиабаты для воздуха.

 

Так как воздух в отверсти расширяется до уровня внешнено давления, то Р2=1атм.

Тогда можно рассчитать :

5/1=(V2/V1)^1,4

Откуда V2/V1=5^(1/1.4)=3,157

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

V2/V1=q1/q2

q2= q1/(V2/V1)

q2=6/3,157=1,9кг/м3

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

Как мы уже договорились ранее, в критическом сечении срабатывается только одна треть энергии при расширения по осевому направлению.

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

 n=(V2/V1)^(1/3) =3,157^(1/3) = 1,467

 от полного разуплотнения до конечного равновесного состояния всего объёма по трём направлениям.

Итого: в критическом сечении плотность газов составит

qкр=q1/ ((V2/V1)/n) =6/(3,157/1,467)=2,79кг/м3

 

Из полученой плотности мы можем рассчитать скорость струи газа Vст в отверстие  из формулы импульса силы:

F=q*S*V^2

Откуда скорость газов в отверстии:

Vст=(F/(q*S))^0,5 =(3,03/(2,79*0,0000049))^0,5 =471м/с

где q –плотность струи в критическом сечении отверсти, S- площадь сечения отверстия, V- скорость воздуха в в критическом сечении отверстия.

То есть замеряемая сила тяга ПРД по законам термодинамики должна достигаться при скорости воздуха V=471м/с, что  в 1,42 раза или на 42% больше скорости звука в воздухе.

Таким образом, мы достаточно точно попали в ранее выдвинутое предположение, что скорость молекул в звуковой волне полностью срабатывает свою энергию по одному направлению в максимальную скорость, которая в 1,41= (2^0,5) раз больше средней скорости фронта звука.

То есть простой и низкозатратный эксперимент способен доказать несостоятельность общепринятых тезисов как КТГ так и теории работы ЖРД.

При этом подтверждается выдвинутая версия неравновесного однонаправленного расширения газов в критическом сечении КС по модели СТГ.

Сравним энергии газа на скорости : Vа=330м/с при Езв=1579Дж/моль

с альтернативной скоростью: Vк=Vа*1,42=471м/с при Ек=3183Дж/моль

Осталось понять с чем эти цифры сравнивать?

Из ранее приведённых формул полной энергии газа из учебник м ы знаем, что на каждое из 3-х направлений пространства приходится 1/3 полной энергии.

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

Егаз= 3*Ек= 3*3183=9549Дж/моль

Также из формулы Егаз= Кр1*R*Т можем получить величину корректировочного коэффициента Кр:

Кр1= Егаз/(R*Т)=9549/(8,31*295)=3,895

Эта величина близка к величине коэффициента энергосодержания Кр из учебника, связанного с показателем адиабаты газа к=1,4 формулой Кр=к/(к-1)=1,4/(1,4-1)=3,5

Весьма близкие получаются значения.

При этом Кр = Кр1 при показателе адиабаты к=1,345.

Тот же расчёт для Езв даст в два раза меньшую энергию, что приведёт к изменению

Кр-зв=(Кр1)/2=1,948, что достигается при совершенно фантастическом показателе адиабаты к=2,06.

ИТОГО: полученная из расчёт по экспериментальным данным скорость истечения сжатого воздуха из отверстия, равная Vк=Vа*1,42=471м/с вполне удовлетворяет предложеной теоретической модели о срабатывании всей энергии газа на одно направление в максимальную скорость молекул в звуковой волне, равную Vмах=Va*1,41.

Также, полученное значение скорости Vк=Vа*1,42=471м/с даёт близкие значения к расчётным формулам из учебник, определяющим величину коэффициента Кр общего энергосодержания газа исходя из величины адиабаты газа.

Тут возникает вопрос о правомочности последовательного расширения по отдельным направлениям  элементарного объёма газа, чтобы в итоге результат был равен одноактному расширению по всем направлениям?

Сойдётся ли ответ  от произведения при 3-х последовательных расширениях с единократным расширением до конечного объёма?

Это легко считается из балансового уравнения адиабаты:

Р1/Р2= (V2/V1)^1,4

Если считать объём V2=V1*W, то конечный объём V2  мы выразим тремя этапами расширения одинаковой кратности:

 V21= V1*(W^(1/3))

V22= V21*(W^(1/3))

V23= V22*(W^(1/3))

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

V23=V1*(W^(1/3)) *(W^(1/3))*(W^(1/3))=V1*(W^(1/3))^3=V1*W= V2

То есть последовательное 3-х кратное расширение на величину

n=(W^(1/3))

даёт в итоге равенство с однократным расширением на величину W, что и следовало доказать.

 

Также можно подставить трёхкратное расширение и в уравнение адиабаты.

Так возведя каждый из трёх множитель в степень К=1,4 в итоге получаем ровно тот же результат на выходе после перемножения показателя  степеней:

(V2-3/V1)^1,4=((W^(1/3)) ^1,4) *((W^(1/3)) ^1,4)*((W^(1/3)) ^1,4)=

=((W^(1,4/3)) *((W^(1,4/3))*((W^(1,4/3)) ^1,4)=

=(W^(1,4/3))^3=(W^(1,4*3/3))=

=W^(1,4)= W^1,4

 

Скорость газов в сопле при изменении давления в ресивере ПРД

Эксперимент №2

Исходные тезисы теории ЖРД  утверждают, что тяга КС определяется прямой зависимостью от давления в КС при постоянстве температуры:

F= Pкс*Sкр*Ккс

Где Ркс-давление в КС; Sкр- площадь критического сечения; Ккс- это эмпирический коээфициент конкретной модели ЖРД с конкретными специфичными свойствами газа (разные коэффициенты адиабаты газа и формы КС), при этом значение Ккс лежит в диапазоне:  1<Ккс< 2. (тут значение 2 соответствует несжимаемой жидкости с коноидальным заходом в отверстие в стенке бака под давлением и полным 100% срабатыванием давления в баке в кинетическую энергию струи)

То есть при одинаковой постоянной температуре газа в КС (ресивере ПРД) тяга линейно растёт от роста давления.

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

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

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

График тяги ПРД с различными соплами.
График тяги ПРД с различными соплами.

Рис.56. Графики продувки насадков различных диаметров и фор в ПРД: 1- отверстие Ф3мм ( две серии замеров разными цветами), 2- отверстие Ф2,5мм ( две серии замеров разными цветами, черный цвет с фаской-соплом), 3-отверстие Ф2мм ( две серии замеров разными цветами), 4- отверстие Ф2мм с соплом Лаваля (на перерасширение струи), форма графика весьма замысловата  и далека от линейной. Тонкие линии даны для ориентации к линейной опраксимации графиков.

Влияние конического насадка «сопла Лаваля» на разгон струи.

Эксперимент №3

В предыдущих этапах эксперимента не было проверено влияние расширения струи в сопле Лаваля на импульс струи газов.

Для развития эксперимента было установлено конусное сопло с большой кратностью расширения:

Sсоп/Sкр >10.

В результате продувки ПРД с коническим соплом появилось несколько интересных результатов:

1.       Показания весов резко снизились по отношению к этому  же отверстию без насадка.

2.       На давлении около 5 атм ещё и появился резкий свист.

3.       Время полного опорожнения ресивера  оказалось одинаковым  с соплом или без него.

Все три результата были неожиданными.

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

Но я нигде не встречал информацию, что у перерасширенной струи снижается собственный импульс!

Кстати, тормозящие перерасширенные сопла Лаваля широко применяют в стрелковом оружии (пулемёты Дегтярёва, КПВТ, Утёс и др., а также  автоматы АКСУ и др.), но чаще  их называют «пламегасители», хотя  они в значительной мере являются ещё и «дульными тормозами» (см.рис.57)

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

3. Ранее приведённый расчёт плотности истекающей струи из ПРД показал, что для выравнивания  поперечного давления до атмосферного нужно расширение сопла всего-то в 2 раза от площади сечения отверстия.

Такое расширение достигается соплом Лаваля в форме небольшой конической фаской с диаметром внешней кромки Dвнеш=1,5*Dотв.

Правда устройство такой фаски не повлияло на замеряемый на весах импульс самой струи (черная кривая график 2 на рис.56), что вполне укладывается в модель расширения конуса сверхзвуковой струи за пределами КС в начальной части сопла Лаваля (рассмотрено в главе про ЖРД). Но прибавка тяги при этом на ПРД (ЖРД) возникать должна.

Для выявления положительного  или отрицательного влияния  сопла разных степеней уширения нужно измерять непосредственно реакцию струи на сам ПРД.

Для проверки изменения тяги самого ПРД необходимо поменять схему эксперимента.

В новом варианте эксперимента будет взвешиваться на весах  сам ПРД, а не ударное воздействие его струи на весы.

Результаты эксперимента с непосредственным взвешивании самого ПРД будут представлены позднее в новой главе.

 

Пулемёт НСВ "Утёс" с раструбным дульным тормозом-пламегасителем.
Пулемёт НСВ "Утёс" с раструбным дульным тормозом-пламегасителем.

Рис.57. Современный  пулемёт НСВ-«Утёс» с раструбным  дульным тормозом-пламегасителем на конце ствола.

Эксперимент №4

Эксперимент с непосредственным взвешиванием тяги ПРД.

(добавлено 17-02-2023)

Представляется новая схема экспериментальной установки (см.рис.58) для проверки результатов предыдущего эксперимента со взвешиванием реакции струи на весах (см.рис.55-56)

Эксперимент с непосредственным взвешиванием тяги ПРД.
Эксперимент с непосредственным взвешиванием тяги ПРД.

Рис.58. Схема непосредственного замера силы тяги струи ПРД с различными насадками.

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

Результат такой продувки нескольких вариантов насадков представлены на графике (см.рис.59)

Графики тяги ПРД с разными насадками
Графики тяги ПРД с разными насадками

Рис.59. Графики продувки насадков различных диаметров и фор в ПРД:

1-4-старый эксперимент (см.рис.56),

5- совпадающие графики (4шт) от продувки отверстия Ф1,5мм  с гладким соплом Лаваля и без него( два эксперимента по две серии замеров разными цветами). Линии графиков практически совпадают, что означает высокую точность эксперимента (воспроизводимость) и полное отсутствие влияния гладкого сопла на тягу ПРД в данном случае.

6- отверстие Ф1,5мм с гибким шлангом и взвешиванием тяги из отверстия ПРД непосредственно на весах ( 5 серий замеров серым цветом). Схема замера даёт неустойчивые результаты из-за несовершенства метода и сильных влияний промежуточного шлага в переходных режимах на коротких  интервалах выпуска воздуха из ресивера.

7- отверстие Ф1,5мм с ребристым соплом на гибком шланге и взвешиванием тяги из отверстия ПРД непосредственно на весах (2 серий замеров голубым цветом).Схема замера даёт неустойчивые результаты из-за несовершенства метода и сильных влияний промежуточного шлага в переходных режимах на коротких  интервалах выпуска воздуха из ресивера.

8- отверстия Ф1,5мм  с ребристым соплом Лаваля ( два замеров красным и зелёным цветами) при взвешивание струи (1-я схема эксперимента). Данное совпадение графиков двух серий замеров в схеме эксперимента со взвешиванием реакции струи на весах говорит о высокой точность и воспроизводимости результата, что в данном случае достигается применением малого сечения выходного отверстия, тем самым обеспечивается достаточно длительное истечение струи при почти постоянном давлении в ресивере.

Тонкие линии даны для ориентации к линейной аппроксимации графиков.

 

Выводы из Эксперимента №4

Данный эксперимент убедительно доказывает несколько положений в теории работы ПРД:

1.       Происходящее в сопле Лаваля различных конфигураций не влияет на истечение струи газов из критического сечения ПРД (время опорожнения ресивера одинаково с соплами и без).

2.       Гладкое сопло Лаваля не создаёт обнаружимой дополнительной тяги на ПРД в диапазоне данных испытательных давлений в ресивере (2,5-7атм), что подтверждается практически совпадающими графиками экспериментов с разными насадками (см.серию графиков №5 на рис.59).

3.       Ребристое сопло Лаваля создаёт сильное торможение потока газа из ПРД в диапазоне данных испытательных давлений в ресивере (2,5-7 атм), то есть использование таких ребристых конусов вполне возможно в качестве «дульного тормоза» на стрелковом оружие (см.рис.57). В качестве ребристого сопла Лаваля в эксперименте использован металлический конус от шариковой ручки, где на внутренней поверхности имеются кольцевые ступеньки и резьба поперёк потока.

P.S.

Уважаемые читатели!

Если вы добрались до этих строк, то этот весьма длинный и непростой текст вам был чем-то интересен.

Свои замечания и комментарии вы можете обращать мне лично по е-почте: imoninpgd@gmail.com

С уважением,

Монин Илья Алексеевич, ктн.

Библиография:

  1. Яскин, А.В.
    Теория устройства ракетных двигателей: учебное пособие /
    А.В. Яскин; Алт. гос. техн. ун-т, БТИ. Бийск: Изд-во Алт. гос.
    техн. ун-та, 2013 262 с.

Ссылка на ПДФ-файл учебника в интернете.

  1. Учебное пособие по конструкции ЖРД из интернета.

    https://mydocx.ru/3-37968.html

Tags:
Hubs:
+2
Comments 52
Comments Comments 52

Articles