ОЦЕНКА ПРОЕКТНЫХ ПАРАМЕТРОВ ДВИГАТЕЛЬНОЙ УСТАНОВКИ КОСМИЧЕСКОГО АППАРАТА С УЧЕТОМ РЕЗУЛЬТАТА ГАЗОДИНАМИЧЕСКОГО МОДЕЛИРОВАНИЯ ТЕЧЕНИЯ ГРАНУЛИРОВАННОГО ТОПЛИВА
ОЦЕНКА ПРОЕКТНЫХ ПАРАМЕТРОВ ДВИГАТЕЛЬНОЙ УСТАНОВКИ КОСМИЧЕСКОГО АППАРАТА С УЧЕТОМ РЕЗУЛЬТАТА ГАЗОДИНАМИЧЕСКОГО МОДЕЛИРОВАНИЯ ТЕЧЕНИЯ ГРАНУЛИРОВАННОГО ТОПЛИВА
Аннотация
При проектировании перспективных энергоустановок, применяющих порошкообразные и гранулированные топлива, проводится численное моделирование неравновесного истечения высококонцентрированной порошково-газовой среды из отверстия форсунки. Для безопасной отработки системы подачи проектируется экспериментальная установка, исключающая горение, тем самым нет необходимости использовать реальный состав топлива. Вместо него используется инертный порошок, идентичный по своим физико-механическим свойствам к реальному составу хлористого калия (KCl). В работе представлены результаты исследования дисперсных характеристик измельченного KCl. Полученные данные использовались при численном моделировании течения гранулированного топлива в форсунке системы подачи двигательной установки. Представлена методика оценки проектных параметров двигательной установки, а на её основе выполнен расчет проектных параметров двигательной установки космического аппарата.
1. Введение
Двигательные установки (ДУ) на гранулированном топливе (ГТ)
, , , — перспективные энергоустановки, использующие порошкообразный высокоэнергетический материал в качестве топлива, обладающие следующими свойствами: многократное включение-выключение, многоразовое использование.На стадии формирования облика разрабатываемого изделия необходимо оценить проектные параметры (давления, расходы, габариты изделия). Двигатели, использующие порошкообразные топлива активно изучаются
, , , , однако сведений о методике оценки проектных параметров ДУ, применяющих гранулированное топливо, практически нет.Надёжность безотказной работы ДУ на ГТ зависит во многом от устойчивости функционирования системы подачи (СП) топливного наполнителя. Для создания высокоэффективных и безопасных энергоустановок система подачи требует критического режима истечения среды из форсунки, при котором процессы, происходящие в камере сгорания (КС) не влияют на процессы, протекающие в системе подачи. Критическое истечение топлива из форсунки достигается при критическом отношении давления в КС и давления на входе в форсунку, т.е., параметры системы подачи напрямую взаимосвязаны с параметрами КС.
Для численного моделирования критического истечения высококонцентрированной порошково-газовой среды из форсунки в настоящее время предложена одномерная математическая модель
. Как известно критическое отношение давлений на входе форсунку и на выходе из нее (Таким образом, условие стабильной работы системы подачи ДУ на ГТ:
где
Из условия (1) следует, что при заданном давлении
Для обеспечения устойчивой работы системы подачи ГТ необходимо провести расчет критических параметров процесса истечения ГТ из форсунки, определяя зависимость
При заданном давлении
При проектировании системы подачи ДУ, зная давление в КС, необходимо определиться с диапазоном давлений перед форсункой
Расходные характеристики СП зависят от выбранной геометрии форсунки (диаметр критического сечения
Численное моделирование предполагает использование дисперсных характеристик рассматриваемого порошкообразного материала (порозность (объем, занимаемый газом в среде порошка), размер частиц), которые могут быть установлены экспериментально.
Для верификации математической модели, представленной в , необходимо разработать экспериментальную установку для отработки процессов подачи ГТ в КС. Изучение процессов подачи не требует организации процессов воспламенения ГТ, поэтому для проведения безопасной отработки предлагается вместо гранулированного топлива использовать порошок, идентичный по своим физико-механическим свойствам к реальному ГТ.
В качестве порошка-«имитатора» выбран измельченный хлорид калия (KCl) . Его физико-механические характеристики близки к предложенному составу ГТ (плотность состава ГТ: ~1900 кг/м3; плотность KCl: 1984 кг/м3). Для повышения представительности расчетов, необходимо уточнение таких характеристик исходного порошка, как дисперсность и порозность.
Данные численного моделирования необходимы для оценки проектных параметров ДУ на ГТ. В качестве примера, полезно рассчитать проектные параметры ДУ на ГТ, применяющейся для космического аппарата (КА), разработка которой представляет актуальную научно-прикладную задачу , .
Целью представляемого исследования является поиск комплексного подхода в логичной последовательности для формирования методики оценки расчётных проектных параметров двигательной установки на гранулированном топливе для космического аппарата. Для достижения указанной цели необходимо решить следующие задачи:
1. Определить дисперсность и порозность исходного KCl, порошка имитатора гранулированного топлива.
2. Определить диапазоны давлений перед форсункой устойчивой работы системы подачи.
3. Изучить влияние геометрии форсунки на параметры критического течения.
4. Сформировать рекомендации для проектных параметров системы подачи (давления в системе, геометрия форсунки).
2. Методы и принципы исследования
2.1. Гранулометрический анализ
Гранулометрический анализ измельченного порошка хлорида калия проводился на экспериментальной лабораторной установке Mastersizer 2000 компании Malver Instruments Ltd. Реальный вид прибора представлен на рисунке 1. Диапазон измерения диаметров частиц порошка составляет от 0,1 мкм до 2000 мкм.

Рисунок 1 - Лабораторный прибор Mastersizer 2000

Рисунок 2 - Гранулометрический анализ пробы
Для исследования порозности хлорида калия использовалось следующее оборудование:
· стакан мерный цилиндрический 500…1000 мл и ценой деления 5 мл;
· точные электронные весы с погрешностью измерения 0,5 г;
· киянка резиновая;
· емкость для хранения порошка объемом не менее 1000…1500 мл.
Параметры порозности и уплотняемости определялись при частоте колебаний 3…4 Гц, амплитуде 5…10 мм и длительности колебаний 2…4 мин у искомого KCl по ГОСТ 4568-95 (100…1000 мкм). Истинная плотность хлорида калия принята
Исследование порошка заключается в определении коэффициентов порозности и способности к уплотняемости путем многократного определения насыпной и уплотненной плотностей материалов. Проводились не менее десяти экспериментов по уплотняемости каждого из порошков, определялись статистические параметры исследования (математическое ожидание (МО), среднеквадратическое отклонение (СКО), критерий Граббса, случайная, неисключённая систематическая погрешности (НСП) и суммарная), приведены порозности и уплотняемости исследуемых порошков с доверительными границами погрешности оценки измеряемой величины. Ниже представлен алгоритм измерения порозности и уплотняемости порошков:
1) проверить отсутствие посторонних веществ в мерном стакане и емкости для хранения порошка, при необходимости промыть водопроводной водой и высушить;
2) в емкость для хранения порошка объемом 1000…1500 мл насыпать один анализируемый порошка;
3) взвесить мерный стакан на электронных весах, зафиксировать результат;
4) из емкости во весь измеряемый объем мерного стакана насыпать анализируемый порошок;
5) мерный стакан с порошком взвесить на электронных весах, из полученного значения вычесть массу стакана, записать результат;
6) произвести уплотнение порошка — на мерный стакан в течение 2…4 мин приложить с помощью резиновой киянки вибро-воздействие, выполняя удары по внешней поверхности емкости;
7) нагрузку прикладывать до тех пор, пока изменение объема не будет менее половины цены деления стакана после 1 мин приложения нагрузки;
8) записать конечный объем, полученный в результате уплотнения порошка;
9) высыпать порошок из мерного стакана в емкость для хранения;
10) рассчитать начальную (до уплотнения) и конечную (после уплотнения) плотности порошка;
11) рассчитать начальную и конечную порозности порошка;
12) повторить пп. 4) – 11) не менее 9 раз;
13) определить статистические параметры исследований;
14) определить доверительные границы погрешности оценки полученных величин.
В п. 3.2 представлен подробный анализ полученных результатов проведенного исследования по определению порозности и уплотняемости порошков.
В ходе исследования порозности и уплотняемости порошков были проведены серии экспериментов по замеру массы, начального и конечного объемов анализируемых порошков, рассчитаны насыпная и уплотненная плотности материалов, определены коэффициенты порозности до и после уплотнения.
Проведено исследование порозности и уплотняемости измельченного хлорида калия. Результаты представлены в таблицах 1–2. Расчет статистических параметров выполнялся при доверительной вероятности
Таблица 1 - Исследование порозности и уплотняемости измельченного KCl
Номер эксперимента | Vнас, мл | Vупл, мл | m, кг |
1 | 1000 | 920 | 1069 |
2 | 1000 | 920 | 1043 |
3 | 1000 | 900 | 1062 |
4 | 990 | 870 | 1054 |
5 | 1000 | 880 | 1053 |
6 | 1000 | 900 | 1078 |
7 | 1000 | 895 | 1065 |
8 | 1000 | 895 | 1079 |
9 | 1005 | 890 | 1087 |
10 | 990 | 870 | 1067 |
xср, Acp | 999 | 894 | 1066 |
S, S(A) | 4,74 | 17,61 | 13,43 |
Sx_{ср} | 1,50 | 5,57 | 4,25 |
ε, ε(P) | 3,39 | 12,59 | 9,60 |
ΘΣ, Θ(P) | 2,50 | 2,50 | 0,50 |
SΘ | 1,44 | 1,44 | 0,29 |
K, (Θ(P))/(S(A)) | 2,00 | 2,15 | 2,23 |
Δ, Δ(P) | 4 | 12 | 9 |
Примечание: прямое измерение параметров
Таблица 2 - Исследование порозности и уплотняемости измельченного KCl
Номер эксперимента | ρнас, г/см3 | ρупл, г/см3 | εнас | εупл | ν |
1 | 1,069 | 1,162 | 0,461 | 0,414 | 0,080 |
2 | 1,043 | 1,134 | 0,474 | 0,429 | 0,080 |
3 | 1,062 | 1,180 | 0,465 | 0,405 | 0,100 |
4 | 1,065 | 1,211 | 0,463 | 0,389 | 0,121 |
5 | 1,053 | 1,197 | 0,469 | 0,397 | 0,120 |
6 | 1,078 | 1,198 | 0,457 | 0,396 | 0,100 |
7 | 1,065 | 1,190 | 0,463 | 0,400 | 0,105 |
8 | 1,079 | 1,206 | 0,456 | 0,392 | 0,105 |
9 | 1,082 | 1,221 | 0,455 | 0,384 | 0,114 |
10 | 1,078 | 1,226 | 0,457 | 0,382 | 0,121 |
xср, Aср | 1,067 | 1,192 | 0,462 | 0,399 | 0,105 |
S, S(A) | 0,01437 | 0,02787 | 0,00542 | 0,01052 | 0,01814 |
Sx_{cp} | - | - | - | - | - |
ε, ε(P) | 0,03250 | 0,06304 | 0,01227 | 0,02379 | 0,04103 |
ΘΣ, Θ(P) | 0,00299 | 0,00372 | 0,00124 | 0,00154 | 0,00370 |
SΘ | - | - | - | - | - |
K, (Θ(P))/(S(A)) | 0,20812 | 0,13341 | 0,22894 | 0,14676 | 0,20380 |
Δ, Δ(P) | 0,033 | 0,063 | 0,012 | 0,024 | 0,041 |
Примечание: косвенное измерение
Таким образом, насыпная порозность хлорида калия составила:
Порозность хлорида калия после уплотнения:
Уплотняемость хлорида калия составляет:
Уплотненная плотность хлорида калия рассчитана по формуле:
Таким образом, уплотненная плотность для следующих расчетов принимается равной 1190,4 кг/м3.
2.3. Численное моделирование течения порошково-газовой среды
В предыдущем разделе определены параметры измельченного порошка KCl, использующегося в качестве имитатора ГТ. Исходные данные для расчета представлены в таблице 3.
Таблица 3 - Исходные данные для численного моделирования истечения KCl с азотом из форсунки
Параметр | Значение |
Давление на входе в форсунку, p0, МПа | 0,5…4,0 |
Радиус критического сечения, rcr, мм | 1,0…2,0 |
Длина участка, z, мм | 21,0 |
Угол конусности форсунки, φ, ° | 24,0 |
Минимальная (предельная) порозность, ε0 | 0,40 |
Средний диаметр по поверхности частиц, d3,2, мкм | 22,24 |
Плотность материала, T, кг/м3 | 1984 |
Температура газа, К | 298,15 |
Газовая постоянная, Rg, Дж/(кг·К) | 296,50 |
Динамическая вязкость, η, Па·с | 17,78·10-6 |
Толщина вытеснения, δ, м | 0,0001 |
Для проведения испытаний для разрабатываемой экспериментальной установки изготовлена партия форсунок, геометрия которых представлена в таблице 2.
Как известно, на критические параметры сильно влияют параметры геометрии форсунки. Изучено влияние диаметра критического сечения при постоянном угле конусности форсунки. Для форсунок с разными диаметрами критического сечения в результате численного моделирования с использованием математической модели получены следующие зависимости от давления на входе в форсунку (рис. 3–6):
· Критического отношения давлений.
· Отношения расходов фаз.
· Критического удельного расхода порошка;
· Неравновесной критической изотермической скорости порошка.
· Координаты начала процесса перехода предельно-уплотненного течения порошка в псевдоожиженное.
Получены следующие зависимости от диаметра критического сечения при постоянном угле конусности форсунки (рис. 7–8):
· Критического удельного расхода порошка.
· Отношения расходов фаза.

Рисунок 3 - Критическое отношение давлений в зависимости от давления на входе в форсунку для разных диаметров критического сечения

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

Рисунок 5 - Отношение расходов фаз в зависимости от давления на входе в форсунку для разных диаметров критического сечения

Рисунок 6 - Неравновесная изотермическая критическая скорость порошка в зависимости от давления на входе в форсунку для разных диаметров критического сечения

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

Рисунок 8 - Отношение расходов фаз в зависимости от диаметра критического сечения при постоянном угле конусности форсунки
Из рисунка 3 также следует, что с повышением давления в системе подачи уменьшаются потери давления в критическом сечении форсунки. До
3. Основные результаты
3.1. Предлагаемые рекомендации
1. При подаче топлива в камеру сгорания должно сохраняться критическое отношение давлений между камерой сгорания и системой подачи (т.е., должно выполняться условие
2. Используя полученную зависимость, можно рассчитать давление, реализуемое в камере сгорания. Колебания давления вверх в КС могут привести к прекращению критического режима истечения порошково-газовой среды из форсунки, что, в свою очередь, может привести к возгоранию топливного бака.
3. При выборе давления в системе необходимо указать коэффициент безопасности, минимальное значение которого может уточняться по результатам огневых испытаний. Можно рассчитать давление в камере сгорания по формуле (4):
где
На рисунке 9 представлен график для определения давления в системе подачи и камере сгорания.
Варьируя давление в камере сгорания, с учетом коэффициента безопасности, по рисунку 9 выбирается минимальное давление в системе подачи, необходимое для достижения критического режима течения топлива из форсунки. Возможен обратный путь: задаваясь давлением в системе подачи определяется максимальное давление КС, при котором сохраняется критический режим истечения топлива из форсунки.

Рисунок 9 - График для определения давления в системе подачи и камере сгорания
Методически предварительная оценка проектных параметров двигательной установки при критическом течении гранулированного топлива в системе подачи предусматривает оценочный расчет ДУ на ГТ для космического аппарата. Исходные данные для расчёта представлены в таблице 4.
Таблица 4 - Исходные данные для проектного расчёта
Параметр | Значение |
Тяга, Н, не менее | 392 |
Время работы двигателя в непрерывном режиме, τ, с, не менее | 1000 |
Удельный импульс ДУ в пустоте, м/с | 2885,5 |
Насыпная плотность топлива, ρt, кг/м3 | 1190,4 |
Номинальное давление в КС, МПа | 1 |
Количество форсунок | 2 |
Коэффициент безопасности системы подачи, Kb | 1,20 |
Коэффициент запаса топлива, Kzb | 1,10 |
Коэффициент запаса вытеснительного газа, Kzg | 1,25 |
Другие исходные данные для расчета – это результаты численного моделирования критического истечения порошково-газовой среды с имитатором гранулированного топлива, а именно:
· Зависимость критического отношения давлений (рис. 3) от давления на входе в форсунку.
· Зависимость критического удельного расхода (рис. 4) от давления на входе в форсунку.
· Зависимость относительного расхода газа (рис. 5) от давления на входе в форсунку.
Давление в камере сгорания задано в исходных данных, согласно которому определяется минимальное давление перед форсункой по графику, представленному на рис. 9 с учетом закладываемого коэффициента безопасности (учитывает возможное повышение давления в КС). Минимальное давление перед форсункой
Кривые на рис.7 предопределяют радиус критического сечения форсунки. Для упрощения расчета задаётся радиус критического сечения 1 мм и далее по графику, представленному на рисунке 4, определяется критический удельный расход гранулированного топлива
По формуле (4) вычисляется необходимая масса топлива:
По формуле (5) необходимая масса вытеснительного газа составит:
Таблица 5 - Результаты расчета проектных параметров
Наименование параметра | Значение |
Масса топлива, кг | 284,984 |
Масса вытеснительного газа, кг | 7,459 |
Требуемое минимальное давление перед форсункой, МПа | 3,662 |
Давление в камере сгорания, МПа | 1 |
Время работы двигателя в непрерывном режиме работы, с | 1000 |
Тяга на непрерывном режиме работы, Н | 747,565 |
Количество форсунок | 2 |
Радиус критического сечения форсунки, м | 0,001 |
Угол, φ, ° | 24,0 |
Радиус входного сечения форсунки, м | 0,005 |
Номинальный массовый расход топлива, кг/с | 0,259 |
Номинальный массовый расход вытеснительного газа, кг/с | 0,0054 |
Удельный импульс тяги двигателя в пустоте, м/с | 2885,5 |
Объем топлива с учетом пор, м3 | 4,18 |
Таким образом, в таблице 5 представлены параметры системы подачи, которая удовлетворяет заданные требования исходных данных: обеспечивает необходимые тяговые характеристики двигателя при сохранении устойчивой работы.
4. Обсуждение
Очевидным является то, что для выполнения заданных требований можно подобрать форсунку с другой геометрией, которая обеспечила бы меньшие значения расходных характеристик. Тем не менее, одной из рассматриваемых задач было предложить методику оценки проектных параметров системы подачи. Чтобы оптимизировать значения проектных параметров необходимо более детально изучить влияние геометрии форсунки на газодинамические параметры критического истечения порошково-газовой среды из форсунки, что позволит выбрать оптимальный вариант геометрии форсунки.
5. Заключение
1. Для проведения экспериментальных исследований и безопасной отработки системы подачи гранулированного топлива предложено использование инертного альтернативного порошка, имитирующего физико-механические свойства предполагаемого состава топлива. В качестве такого порошка выбран измельченный (безопасный) хлорид калия. Методом лазерной дифракции определен уровень дисперсности (средний поверхностный диаметр частиц продукта).
2. Представлены результаты численного моделирования критического истечения хлорида калия в среде инертного газа (азота) из форсунки. Для форсунок с разными диаметрами критического сечения получены следующие зависимости от давления на входе в форсунку:
· Критического отношения давлений.
· Отношения расходов фаз.
· Критического удельного расхода порошка.
· Неравновесной критической изотермической скорости порошка.
· Координаты начала процесса перехода предельно-уплотненного течения порошка в псевдоожиженное состояние.
Получены следующие зависимости от диаметра критического сечения при постоянном угле конусности форсунки:
· Критического удельного расхода порошка.
· Отношения расходов фаза.
3. Приведена методика оценки проектных параметров двигательной установки на гранулированном топливе. Выполнен оценочный расчет проектных параметров двигательной установки на гранулированном топливе космического аппарата.
