Давление на выходе из сопла слишком высокое?

Ради забавы я занимался простым ракетостроением и хотел рассчитать характеристики двигателя спуска лунного модуля (LMDE). LMDE использует Aerozine-50 (UDMH) и тетроксид диазота ( Н 2 О 4 ). Я хотел рассчитать выходное давление сопла. Формула для расчета скорости выхлопа на выходе:

В е "=" Т р М 2 γ γ 1 ( 1 ( п е п ) ( γ 1 ) / γ )
где:

  • В е скорость выхлопа. В этом случае он равен 3225 м / с .

Рассчитано по Ф "=" м ˙ В е или

В е "=" Ф м ˙
где Ф является 11 , 965 Н и массовый расход 3,71 к г / с при тяге 25 %.

  • Т - абсолютная температура газа на входе. В этом случае равно 294,216 К . Получено отсюда . Посмотрите это, найдя «Номинальная температура топлива на входе в форсунку».
  • р универсальная газовая постоянная, которая 8314,5 Дж / ( к м о л · К )
  • М - вес молекулярного газа топлива. В данном случае это 20.58 к г / к м о л . Получено отсюда . Обратите внимание, что LM имеет соотношение смеси 1,6, поэтому используйте строку с маркировкой 1,6. Кроме того, давление в камере LM составляет 120 фунтов на квадратный дюйм (около 8 атм). Этот источник также имеет коэффициент удельной теплоемкости.
  • γ - коэффициент изоэнтропического расширения, также известный как коэффициент удельной теплоемкости. В данном случае это 1,232 . (Получено из источника выше).
  • п е давление на выходе из сопла (в Паскалях).
  • п давление газа на входе. В данном случае это 3 , 010 , 000 п а (или 437 фунтов на квадратный дюйм). Получено отсюда .

пытаюсь найти выходное давление( п е ), поэтому я изменил формулу на:

п е "=" п ( 1 В е 2 М ( γ 1 ) Т р 2 γ ) γ γ 1

Замените все значения. Числа и единицы 3.01 × 10 6 п а с с а л с , 3225 м / с , 294,261 К , 8314,5 Дж / к м о л К , 20.58 к г / к м о л , и γ "=" 1,232 является соотношением.

п е "=" 3.01 × 10 6 ( 1 ( 3225 ) 2 ( 20.58 ) ( 1,232 1 ) ( 294,261 ) ( 8314,5 ) ( 2 ) ( 1,232 ) ) 1,232 1,232 1

если я подставляю все значения в уравнение, я получаю 1,105 × 10 11 п а . Около 16,8 миллионов фунтов на квадратный дюйм!

Другой способ рассчитать выходное давление с помощью Ф "=" м ˙ В е + ( п е п а ) А или:

п е "=" Ф м ˙ В е А + п а

А площадь среза сопла ( 1,9 м 2 ) и п а это атмосферное давление, которое 0 п а потому что LMDE находится в вакууме. Результат 0,13158 п а или 0,00001885 п с я .

Вопрос: Что здесь происходит? Для первого уравнения я получаю такое нереальное высокое давление на выходе из сопла. Но во втором уравнении я получаю такое низкое давление на выходе (почти вакуум). Разве эти значения не должны быть примерно одинаковыми?

Разве это не м ˙ 3,71 кг/с при тяге 25%. Я также использовал 11 965 Н, что составляет примерно 25% от полной мощности LM.
Сделанный. Виноват. Это совершенно вылетело из головы, пока вы не упомянули 15 кг/с.
Можете ли вы убедиться, что р должна быть универсальной газовой постоянной в первом уравнении? Я думаю, что есть конфликтующие соглашения , например: р "=" р ты М где р ты универсальная газовая постоянная и М молекулярная масса:
Хм. Может, все-таки попробовать, вдруг Википедия есть Википедия?
Мне кажется, вы игнорируете горение. Вы используете комнатную температуру как температуру застоя газа. Температура застоя в камере сгорания должна быть намного больше. 2500 градусов K дает правильную скорость выхлопа при условии, что все остальные значения верны, а давление на выходе равно нулю.
Посмотрите на температуру адиабатического пламени в той же ссылке, из которой вы получили молекулярную массу. По крайней мере, поэтому первое уравнение не проверялось. Я позволю вам решить переработанную версию и сообщить, стало ли лучше.
@OrganicMarble Кажется, стало намного лучше. Замена 2500 градусов К дала результат около 283 Па (0,041 фунта на кв. дюйм). Это вроде бы правильно, но опять же. Для последнего уравнения я получил 0,13 Па (0,000019 фунтов на квадратный дюйм).
Прохладный! Получение хорошего значения для Tc довольно сложно и для этого уровня вычислений чревато предположениями. Что произойдет, если вы используете временный файл с того веб-сайта, с которого вы получили другие свойства? Они смотрели выше и выше лучше для вашего ответа. Я только что вытащил 2500 из своего фундамента.
Для проверки рассчитайте Tc, используя первое уравнение и требуемое значение Pe.
@OrganicMarble Но слишком высокий уровень усугубляет ситуацию. Например, я использовал 3000 К, и это было около 40 000 Па. Сейчас немного поздно, и я продолжу это завтра.
Еще одна вещь, прежде чем я назову это ночью. P в первом уравнении должно быть давлением в камере сгорания. Я думаю, что это было около 100 фунтов на квадратный дюйм, но, возможно, меньше при 25% дроссельной заслонки. Не уверен, откуда вы взяли 437. Помните, что выход из камеры сгорания — это вход в сопло де Лаваля, и именно эти свойства вы должны использовать в уравнении скорости на выходе. Давление в камере и температура в камере.
Я решил первое уравнение для Tc, учитывая Pc 120 фунтов на квадратный дюйм / 82700 Па (хотя я не уверен, что такое Pc при дроссельной заслонке 25%), и получил 2558 R. Сейчас я собираюсь написать ответ.
Эм, 2558 К. Я просто не чувствую единиц СИ, мой плохой.

Ответы (1)

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

Самый важный вывод заключается в том, что уравнение скорости на выходе Де Лаваля (найдено в Википедии, а также уравнение 3-14 в Саттоне, 4-е издание) требует, чтобы вы использовали свойства камеры сгорания в качестве входных данных в уравнение. Это включает в себя повышение температуры из-за сгорания.

введите описание изображения здесь

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

Я также слышал, что есть различные онлайн-калькуляторы, которые могут это сделать, но у меня нет опыта работы с ними. В комментарии к этому ответу упоминается « propep ». Комментарий к этому ответу на этот вопрос. Количественные графики v, T, p и положения от камеры через сопло до окружающей среды для нескольких канонических современных двигателей? упоминает " cpropep "

Рассчитать температуру сгорания и давление сгорания. Это ключ. Большое спасибо. Кажется, я рассчитывал давление топлива, поступающего в камеру. Вычислить температуру горения действительно сложно. К счастью, я нашел этот источник , в котором указаны температура и давление сгорания. (увидеть его можно, найдя «TRW TR-201»). Итак, это около 703 265 Па и 2 973 К. Я забыл о расчете на 25%, но вместо этого на 100%. Таким образом, это около 47 860 Н и массовый расход 15 кг/с. В е 3190,7 м/с.
(Продолжение комментария выше (ограничение слов)). Используя приведенные выше значения, это около 143,76 Паскалей (0,021 фунта на кв. дюйм). Это кажется правильным ответом. Можно ли с уверенностью предположить, что при тяге 25% давление на выходе равно давлению на выходе при полной тяге, деленному на 4? Я проголосовал за, но не могли бы вы объяснить, почему второе уравнение (при полной тяге, то есть 47,86 кН и 15 кг / с) дает около 6,67 Па (0,00096 фунтов на квадратный дюйм). Второй результат примерно в 22 раза ниже. Я бы не считал это "приблизительным".
@StarMan Не делая много работы (помните, что я ленив), я думаю, что эти два уравнения имеют разные предположения. Уравнение Ve предполагает изоэнтропический поток и, возможно, другие вещи, а уравнение тяги не делает никаких предположений о потоке.
Я понимаю. Итак, второе уравнение тяги для нахождения п е чем-то похож на п "=" Ф / А как его получение количества и силы выбрасываемого вещества и деление его на площадь.
Попробуйте рассчитать разницу в тяге между двумя давлениями в выходной плоскости. Вероятно, не так много. Радуйтесь, что они так хорошо согласны со всеми предположениями, заложенными в эти простые уравнения! Примечание. Саттон говорит, что постоянное отношение удельных теплоемкостей является еще одним предположением.
Истинный. Разница около 260 Н. В любом случае, большое спасибо.