EVALUATION OF THE DESIGN PARAMETERS OF THE SPACECRAFT'S PROPULSION SYSTEM TAKING INTO ACCOUNT THE RESULTS OF GAS-DYNAMIC MODELLING OF GRANULATED FUEL FLOW

Research article
DOI:
https://doi.org/10.60797/IRJ.2026.165.71
EDN:
RRQKUM
Suggested:
19.08.2025
Accepted:
06.03.2026
Published:
17.03.2026
Issue: № 3 (165), 2026
Rightholder: authors. License: Attribution 4.0 International (CC BY 4.0)
14
1
XML
PDF

Abstract

When designing perspective power facilities that use powdered and granular fuels, numerical modelling of the non-equilibrium flow of a highly concentrated powder-gas medium from the nozzle hole is carried out. To ensure the safe operation of the delivery system, an experimental setup is designed that excludes combustion, thereby eliminating the need to use real fuel. Instead, an inert powder is used that is identical in its physical and mechanical properties to the real composition of potassium chloride (KCl). The paper presents the results of a study of the dispersion characteristics of crushed KCl. The acquired data were used in numerical modelling of the flow of granular fuel in the nozzle of the propulsion system. A method for evaluating the design parameters of the propulsion system is presented, and based on it, the design parameters of the spacecraft's propulsion system are calculated.

1. Введение

Двигательные установки (ДУ) на гранулированном топливе (ГТ)

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

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

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

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

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

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

Таким образом, условие стабильной работы системы подачи ДУ на ГТ:

(1)

где

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

Из условия (1) следует, что при заданном давлении

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

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

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

При заданном давлении

, вычисленное значение
будет являться минимально допустимым отношением давлений между КС и СП ГТ.

(2)

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

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

Расходные характеристики СП зависят от выбранной геометрии форсунки (диаметр критического сечения

, угол к продольной оси
.). Таким образом, необходимо провести расчеты критических параметров СП с заданными параметрами гранулированного топлива, в зависимости от давления перед форсункой и ее геометрии. Для упрощения задачи предусмотрен угол конусности форсунки постоянным
= 22,5°.

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

Для верификации математической модели, представленной в

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

В качестве порошка-«имитатора» выбран измельченный хлорид калия (KCl)

. Его физико-механические характеристики близки к предложенному составу ГТ (плотность состава ГТ: ~1900 кг/м3; плотность KCl: 1984 кг/м3). Для повышения представительности расчетов, необходимо уточнение таких характеристик исходного порошка, как дисперсность и порозность.

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

,
.

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

1. Определить дисперсность и порозность исходного KCl, порошка имитатора гранулированного топлива.

2. Определить диапазоны давлений перед форсункой устойчивой работы системы подачи.

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

4. Сформировать рекомендации для проектных параметров системы подачи (давления в системе, геометрия форсунки).

2. Методы и принципы исследования

2.1. Гранулометрический анализ

Гранулометрический анализ измельченного порошка хлорида калия проводился на экспериментальной лабораторной установке Mastersizer 2000 компании Malver Instruments Ltd. Реальный вид прибора представлен на рисунке 1. Диапазон измерения диаметров частиц порошка составляет от 0,1 мкм до 2000 мкм.

Лабораторный прибор Mastersizer 2000

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

На рис. 2 представлены результаты гранулометрического анализа измельченного порошка хлорида калия на установке Mastersizer 2000. Таким образом, средний поверхностный диаметр измельченного хлорида калия, используемый при расчете критических параметров по математической модели, описанной в
, составляет 22,24 мкм.
Гранулометрический анализ пробы

Рисунок 2 - Гранулометрический анализ пробы

2.2. Анализ порозности и уплотняемости порошковых материалов

Для исследования порозности хлорида калия использовалось следующее оборудование:

· стакан мерный цилиндрический 500…1000 мл и ценой деления 5 мл;

· точные электронные весы с погрешностью измерения 0,5 г;

· киянка резиновая;

· емкость для хранения порошка объемом не менее 1000…1500 мл.

Параметры порозности и уплотняемости определялись при частоте колебаний 3…4 Гц, амплитуде 5…10 мм и длительности колебаний 2…4 мин у искомого KCl по ГОСТ 4568-95 (100…1000 мкм). Истинная плотность хлорида калия принята

= 1984 кг/м3.

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

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. Расчет статистических параметров выполнялся при доверительной вероятности

= 0,95 по данным всех десяти экспериментов, так как максимальное и минимальное значения прошли проверку на исключение грубых погрешностей по критерию Граббса. Коэффициент Стьюдента
= 2,262 при доверительной вероятности
= 0,95 и степени свободы
= 9.

Таблица 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

Примечание: косвенное измерение

Таким образом, насыпная порозность хлорида калия составила:

= 0,462±0,012.

Порозность хлорида калия после уплотнения:

= 0,399±0,024.

Уплотняемость хлорида калия составляет:

= (10,5±4,1)%.

Уплотненная плотность хлорида калия рассчитана по формуле:

(3)

Таким образом, уплотненная плотность для следующих расчетов принимается равной 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,5 МПа разуплотнение потока порошка происходит вблизи критического сечения (
≤ 0,2 мм), далее координата разуплотнения стремительно изменяется, потери давления резко уменьшаются. Этот факт необходимо учитывать при выборе давления в камере сгорания.

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

3.1. Предлагаемые рекомендации

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

<
).

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

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

(4)

где

коэффициент безопасности.

На рисунке 9 представлен график для определения давления в системе подачи и камере сгорания.

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

График для определения давления в системе подачи и камере сгорания

Рисунок 9 - График для определения давления в системе подачи и камере сгорания

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

Методически предварительная оценка проектных параметров двигательной установки при критическом течении гранулированного топлива в системе подачи предусматривает оценочный расчет ДУ на ГТ для космического аппарата. Исходные данные для расчёта представлены в таблице 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 с учетом закладываемого коэффициента безопасности (учитывает возможное повышение давления в КС). Минимальное давление перед форсункой

, которое обеспечит критический режим истечения топлива из форсунки, составит 3,662 МПа.

Кривые на рис.7 предопределяют радиус критического сечения форсунки. Для упрощения расчета задаётся радиус критического сечения 1 мм и далее по графику, представленному на рисунке 4, определяется критический удельный расход гранулированного топлива

, который составит 41233,27 кг/(с·м2). Разделив это значение на площадь критического сечения и умножив его на количество форсунок, получается номинальный расход топлива (на непрерывном режиме), который составит 0,259 кг/с. Пользуясь графиком на рисунке 5 можно определить номинальный расход газа на непрерывном режиме составит 0,0054 кг/с.

По формуле (4) вычисляется необходимая масса топлива:

(4)

По формуле (5) необходимая масса вытеснительного газа составит:

(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. Приведена методика оценки проектных параметров двигательной установки на гранулированном топливе. Выполнен оценочный расчет проектных параметров двигательной установки на гранулированном топливе космического аппарата.

Article metrics

Views:14
Downloads:1
Views
Total:
Views:14