Методика моделирования процесса течения жидкости в смесителе, состоящем из двух малорасходных центробежных форсунок Текст научной статьи по специальности «Физика»

УДК 621.454:532.5 Б01: 10.18287/2541-7533-2018-17-4-141-154

МЕТОДИКА МОДЕЛИРОВАНИЯ ПРОЦЕССА ТЕЧЕНИЯ ЖИДКОСТИ В СМЕСИТЕЛЕ, СОСТОЯЩЕМ ИЗ ДВУХ МАЛОРАСХОДНЫХ

ЦЕНТРОБЕЖНЫХ ФОРСУНОК

© 2018

Е. В. Семкин начальник группы конструкторского отдела;

АО «Научно-исследовательской институт машиностроения», г. Нижняя Салда; шimash@,Пstш

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

Жидкостный ракетный двигатель малой тяги; смеситель; гидравлический тракт; центробежные форсунки; моделирование течения; методика расчёта.

Цитирование: Семкин Е.В. Методика моделирования процесса течения жидкости в смесителе, состоящем из двух малорасходных центробежных форсунок // Вестник Самарского университета. Аэрокосмическая техника, технологии и машиностроение. 2018. Т. 17, № 4. С. 141-154. DOI: 10.18287/2541-7533-2018-17-4-141-154

Введение

При проектировании смесителей жидкостных ракетных двигателей малой тяги (ЖРДМТ) тягой 55 Н и менее, работающих на самовоспламеняющихся компонентах топлива АТ-НДМГ, ставятся повышенные требования к организации процесса перемешивания компонентов топлива за время жидкофазной индукции [1; 2]. Малые расходные характеристики и повышающиеся требования к динамическим характеристикам таких двигателей (суммарный расход компонентов топлива менее 0,018 кг/с) не позволяют иметь большое количество смесительных элементов. Обычно ограничиваются одним смесительным элементом, который должен обеспечивать равномерное и полное перемешивание компонентов топлива в камере сгорания до заданного массового соотношения за короткое время жидкофазной индукции тжф = 1 • 10 4 с для их максимально

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

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

га. Учёт всех действующих сил на компоненты топлива при их перемешивании является трудоёмкой и сложной задачей. В настоящее время оптимизация работы смесительного элемента в процессе создания ЖРДМТ проводится по результатам проектных расчётов и десятка испытаний на натурных компонентах топлива, что связано с большими финансовыми и временными затратами. Благодаря развитию вычислительной техники в инженерных расчётах в настоящее время появились модели, отражающие работу смесительных элементов и форсунок [2 - 4], реализованные в коммерческих машинных кодах, например Ansys SFX, Ansys Fluent, Flow Vision. Применение моделей работы смесителей и форсунок ещё на проектном этапе работ позволяет провести оптимизацию параметров их работы и существенно сократить финансовые и временные затраты на проведение разработки ЖРДМТ. В статье приведены результаты разработки методики расчёта смесительного элемента ЖРДМТ тягой 10 - 15 Н.

Объект исследования

В качестве объекта исследований принят смеситель камеры перспективного ЖРДМТ тягой 10 - 15 Н разработки АО «НИИМаш». В результате проведённых исследовательских работ было собрано большое количество экспериментальных данных, полученных при проведении гидравлических испытаний смесителя и огневых испытаний ЖРДМТ со смесителем данной конструкции, что делает возможным сравнение результатов расчёта с экспериментами. Смеситель состоит из двух соосно установленных центробежных форсунок. Гидравлический тракт смесителя и присоединённая камера сгорания, используемые в расчёте, приведены на рис. 1.

При перепаде давления 6,8 кгс/см2 расход воды через внутреннюю форсунку составляет 1,85 г/с. Расход воды через внешнюю форсунку составляет 2,54 г/с при перепаде давления 6,1 кгс/см2. Закрутка жидкости во внутренней и внешней форсунках производится в разные стороны. Стенки гидравлического тракта смесителя выполнены с регулярной шероховатостью 0,8 мкм.

Рис. 1. Геометрия расчётной области смесительного элемента двигателя тягой 10-15 Н

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

Постановка задачи

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

Для решения задачи используется двухскоростная модель описания течения двухфазного потока [5]:

d

-=Г {raPa) + V.{raPoUa) = 0, (1)

dt

д

dt (roPoUo) + V.(ro(PoUo®Uo)) =

dt (2)

= Sm buoy + Faß -raVpa + V • (Oo (VVUo + (VUo))) + Mo,

3

P = ZraPa , (3)

a=1

ß = Z raßa , (4)

3

r

a=1

Z re= 1, (5)

pa= p для всех а - 1, 2 или 3 . (6)

Здесь а - составляющие части двухфазного потока, в решаемой задаче это вода, поступающая в расчётную область через внутреннюю форсунку 1 и внешнюю форсунку 2, и воздух, заполняющий расчётную область на начало расчёта 3;

ra - объёмная доля контрольного объёма расчётной области, занятая фазой а ;

SM buoy - источник импульса силы тяжести в расчётной области, рассчитываемый

по формуле

SM _buoy ={p~pref )g , (7)

где p - текущая плотность в расчётной области, prf - плотность жидкости в расчётной области на начало расчёта, задаваемая пользователем;

Fa/i - сила поверхностного натяжения, действующая на границе раздела фаз, рассчитываемая по формуле

а=1

Faß faß^aß, (8)

ГДе faß = -täaßnaß + V& , (9)

^aß=|V^aß| , (10)

e - сила поверхностного натяжения воды, при контакте с воздухом на границе раздела фаз; naß - вектор нормали на поверхности контакта фаз, направленный от первичной

жидкости к вторичной; kaß - поверхностное искривление поверхности раздела фаз,

определённое как

kaß=V• naß , (11)

Vs - оператор градиента по поверхности контакта фаз.

При проведении решения считается, что обмен массами между фазами протекает равновесно. Член Ma в (2) описывает передачу импульса в процессе взаимодействия фаз потока и представляется как

Ma= Miß . (12)

Здесь Miß учитывает силу сопротивления при взаимодействии двух фаз и определяется из соотношения

MDß = CDPaßAaß\Uß- UjfUß- Ua), (13) где CD - коэффициент сопротивления; pDß - плотность на границе раздела фаз.

Paß= DPa+ rßPß , (14)

Aaß=|Vra| . (15) В системе уравнений (1) - (6) вязкость задаётся как

ßa=ßRa+ßTa , (16)

где /лК - динамическая вязкость; /лТ - турбулентная вязкость по изотропной модели

турбулентности к-м>, при расчёте турбулентность считалась гомогенной.

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

ка. В настоящей работе использована стандартная зависимость Сп от числа Рейнольдса вида [6]:

Сп = — -(1 + 0.15 • Яе0 687), (17)

п Reг V

где Reг. - число Рейнольдса, рассчитанное по относительной скорости фаз потока в расчётной области, г = 1 или 2.

Коэффициент сопротивления в формуле (17) получают отдельно для пар жидкостей, составляющих двухфазный поток в расчётной области: а = 1 и а = 3 (пара «вода-воздух»); а = 2 и а = 3 (пара «вода-воздух»); а = 1 и а = 2 (пара «вода-вода»).

Для пары «вода-воздух» число Рейнольдса в (17) определялось выражением

Яе, = , (18)

^воздух

где Аи - относительная скорость воды и воздуха; ёт - характерный размер течения в паре «вода-воздух»; Уваздух - кинематическая вязкость воздуха.

Для пары «вода-вода» число Рейнольдса в формуле (17) задавалось зависимостью

Кег , (19)

^ воды

где АУ - относительная скорость воды, поступающей в расчётную область через внутреннюю форсунку, и воды, поступающей в расчётную область через внешнюю форсунку; ёп - характерный размер течения в паре «вода-вода»; уводы - кинематическая вязкость воды.

Таким образом, согласно (17), (18) и (19) на коэффициент сопротивления в расчётной области влияют относительные скорости Аи и АУ, входящие в (18) и (19), и характерные размеры течения ёт и ёп. Относительные скорости Аи, АУ вычисляются в процессе решения поля скоростей жидкостей в расчётной области. Значения характерных размеров течения йт и ёп при расчёте остаются постоянными. Критерием правильности выбора величин ёт и ёп является сравнение полученной в результате расчёта и определённой с помощью фотографирования картин течения. В расчётах значения ёт и ёп изменялись в диапазонах:

ёт от 5 до 150 мкм; йп от 5 до 5000 м.

Для сглаживания решения и исключения деления на ноль в (17) коэффициент сопротивления пары «вода-воздух» рассчитывался по алгоритму, показанному на рис. 2.

Рис. 2. Алгоритм расчёта коэффициента сопротивления Св :

и

X ,У ,1

и.

X ,У ,1

и

X ,У ,1

- проекции скорости воды, поступающей в расчётную область

через внутреннюю и наружную форсунки, и воздуха на координатные оси X, У, 1; АиХГ,1 - проекции относительной скорости воды на координатные оси X, У, 1;

Ли, Аи2 - относительные скорости воды и воздуха в расчётной области

На рис. 3 показан алгоритм определения коэффициента сопротивления для пары «вода-вода».

Рис. 3. Алгоритм расчёта коэффициента сопротивления : А Vх,У,1 - проекции относительной скорости воды на координатные оси X, У, 1; А V - относительная скорость воды в расчётной области

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

В расчётной области была сконструирована специальная расчётная сетка [2] с количеством расчётных узлов порядка 10 000 000. Расчёт проводился в лаборатории суперкомпьютерного моделирования ЮУрГУ.

Результаты расчёта

На рис. 4 - 7 приведены результаты расчёта распыла воды смесителем при изменении ёт от 20 до 150 мкм с фиксированным значением ёп = 20 мкм (формулы (18) и

(19)).

Уо1ите Ргасйсо 150киг(ас» 2 Пдиг* 1 нн 1 ОООе-ООО

Усяите Ргасмоп |$(н,шГк* 2 Ндиге 2 I 000е»000

О

а б

Рис. 4. Решение на момент времени А/ = 14,3 мс от начала расчёта, йт =30 мкм: а - структура течения у корня факела распыла и на выходе из гидравлического тракта смесителя; б - структура свободного факела распыла в объёме расчётной области, прилегающей к смесителю. Объёмная доля воды - 0,1

б

Рис. 5. Решение на момент времени А/ = 10 мс от начала расчёта, йт = 50 мкм: а - структура течения у корня факела распыла и на выходе из гидравлического тракта смесителя; б - структура свободного факела распыла в объёме расчётной области, прилегающей к смесителю. Объёмная доля воды - 0,1

а

УсЯите Ргесйоп

СапЧоиг 1 Г|дигп 1 I г 1 000е+000

Уотте Ргасооп Ндиг*2 000е»000

б

Рис. 6. Решение на момент времени Аt = 15 мс от начала расчёта, =100 мкм: а - структура течения у корня факела распыла и на выходе из гидравлического тракта смесителя; б - структура свободного факела распыла в объёме расчётной области, прилегающей к смесителю. Объёмная доля воды - 0,1

а

¥№1ег2 \ЛЯите Рхайюп СаШош 2 Ваига 1

1 ооое*оао

0 ОООе-СЮ |

8 000е-001 МЮОе-ОСЛ 6 ОООе-ОСН 5 000е-001 4 000е-001 3 000е-001

2 ОООвЛЮ1 ■ 1 000в-001

1.ооа*-а15

2 Усяите Ргасиол

2 Гдигя 2 таг 1 СЮОе*ОйО

У,

а б

Рис. 7. Решение на момент времени Аt = 6,9 мс от начала расчёта, йт = 150 мкм: а - структура течения у корня факела распыла и на выходе из гидравлического тракта смесителя; б - структура свободного факела распыла в объёме расчётной области, прилегающей к смесителю. Объёмная доля воды - 0,1

На рис. 8 - 12 приведены результаты расчётов распыла воды смесителем при изменении в диапазоне от 600 мкм до 5000 м (формула (19)).

тЛ&2 \Zo1ume ЯгасГюо СогИоиг 2 Р|дчг9 1 — 1 000е»000 в 000е-001 8 000е-001 7 ООСЫЮ1 6 000е-001 ■ 5 000е-001 4 000^-001 3 000е-001 2 000е-001 1 000*ТО1

^ ихмв-015

улПеП УсянтеРгасмп |5тигГк* I Пдига 2 «г 1 ОООе+ООО

5 000е-001

а б

Рис. 8. Решение на момент времени А/ = 15,0 мс от начала расчёта, йт =20 мкм, йп =600 мкм: а - структура течения у корня факела распыла и на выходе гидравлического тракта смесителя; б - структура свободного факела распыла в объёме расчётной области, прилегающей к смесителю. Объёмная доля воды - 0,1

•"¡Лъ<2 Уо1ите Ргас(юп СогИоиг 2 Р^чт®1

I—г 1 000е»000

в 000е-001 8 [Ю0е-[Ю1 ?.000е-001 6 0006-001 5 000е-001 4 0009-001 3 000в-001 2 ОООе-<Ю1 1 000*-001 1.000*4)15

«»1ег2 Уедите РгаеОап 19в«игГк* 2 Г чзчг» г 1 ООО*.000

а б

Рис. 9. Решение на момент времени А/ = 15,0 мс от начала расчёта, =50 мкм, =5,0 м: а - структура течения у корня факела распыла и на выходе из гидравлического тракта смесителя; б - структура свободного факела распыла в объёме расчётной области, прилегающей к смесителю. Объёмная доля воды - 0,1

*Чйег1 УЫите Ргайюо СогИоиг 1 Пдига 1 ■ I СЙОе'ООО в 000е-001 8 000е-001 7.000ЫЮ1 6 000е-001

5 000е-001 4.0009-001 3 000е-001 2 ОООМЮ1 1 000*ТО1 1.000*-А 15

•лга!&|1 Уо1ите Ргас{(оп (вокиг^яс* 1 Пдига 2 ■ ОООегООО

а б

Рис. 10. Решение на момент времени А/ = 25,0 мс от начала расчёта, йт=50 мкм, йп=5,0 м: а - структура течения у корня факела распыла и на выходе из гидравлического тракта смесителя; б - структура свободного факела распыла в объёме расчётной области, прилегающей к смесителю. Объёмная доля воды - 0,1

1 0СЮв-ДО1 1.000*-615

/га1е<2 УЬкпн РгасЬоп

|$д*игГк* 2 Пдш 2 СЮОе-ООО

б

Рис. 11. Решение на момент времени At = 24.0 мс от начала расчёта. йт=50 м. йп=5000.0 м: а - структура течения у корня факела распыла и на выходе из гидравлического тракта смесителя; б - структура свободного факела распыла в объёме расчётной области. прилегающей к смесителю. Объёмная доля воды - 0.1

а

у*га1ес1 "/спите Ргасиюп

СопШиг 1 Пдит 1 г-г 1 000е»000 3 ОООе-ОО I

а ооо*-оо1 7 ОО Ое-001 а 000е-001

5 000е-001 4.0ОСЫЭО1 3 000в-001 I 2 000*401 т 1 ОООе-О01 ® 1.000»-015

мы

Аагег2 Усяигле Ргасиап

1хо»игГвсв 2 Ндига 2 ОООегООО

б

Рис. 12. Решение на момент времени At = 40.0 мс от начала расчёта. йт=50 мкм. йп =5000.0 м: а - структура течения у корня факела распыла и на выходе из гидравлического тракта смесителя; б - структура свободного факела распыла в объёме расчётной области. прилегающей к смесителю. Объёмная доля воды - 0.1

а

На рис. 4, а; 5, а; 6, а; 7, а; 8, а; 9, а; 10, а; 11, а; 12, а изменение объёмной доли воды, поступающей через внешнюю и внутреннюю форсунки смесителя, поставлено в соответствие с изменением цвета на рисунках. По изменению цвета в факеле распыла можно судить об изменении объёмной доли воды в конкретной точке расчётной области, а также о степени перемешивания между собой воды, подаваемой через внешнюю и внутреннюю форсунки. Так, зелёный цвет в факеле распыла смесителя свидетельствует о том, что в данной точке расчётной области находится половина объёмной доли воды, поступившей из внутренней и внешней форсунок, т.е. произошло полное перемешивание.

Значения расходов воды через внешнюю и внутреннюю форсунки, установившихся в процессе расчётов (рис. 4 - 12), приведены в таблице.

Таблица. Значения расходов воды через внешнюю и внутреннюю форсунки

рисунок № Gвнутр, г/с G внешн /с Точность, %

внутренняя форсунка внешняя форсунка

4 1,800 2,661 2,7 4,8

5 1,797 2,659 2,9 4,7

6 1,797 2,663 2,9 4,8

7 1,800 2,697 2,7 6,2

8 1,797 2,658 2,9 4,6

9 1,800 2,662 2,7 4,8

10 1,797 2,649 2,9 4,3

11 1,797 2,648 2,9 4,3

12 1,797 2,677 2,9 5,4

В таблице приняты обозначения: Gвнутр _ секундный массовый расход воды через внутреннюю форсунку смесителя; Gвнешн - секундный массовый расход воды через внешнюю форсунку смесителя.

Анализ результатов расчётов

Из данных, представленных на рис. 4 - 7, следует, что увеличение параметра dm

от 20 до 150 мкм приводит к изменениям в основном в структуре течения у корня факела распыла и на выходе из гидравлического тракта смесителя. С переходом к величине dm от 100 до 150 мкм пропадает зона «слипания» воды внешней и внутренней форсунок смесителя в корне факела распыла у сопла внешней форсунки.

Кроме того, в диапазоне dm от 20 до 150 мкм выявлено некоторое изменение

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

С ростом dn от 600 мкм до 5000 м в диапазоне dm от 20 до 50 мкм (рис. 8 - 12)

зона смешения воды внешней и внутренней форсунок сначала приближается к внешней форсунке (рис. 8, а; 9, а), а затем проникает в зону «слипания» плёнки воды внешней и внутренней форсунок (рис. 10, а; 11, а; 12, а). Как следует из данных на этих рисунках, объёмные доли воды внешней и внутренней форсунок уже на диаметре среза сопла внешней форсунки составляют 0,5, что свидетельствует о завершении перемешивания воды в указанной точке расчётной области. Зоны смешения воды внешней и внутренних форсунок на рис. 10, а; 11, а; 12, а, полученные при значениях dn = 5,0 м, dn =5000

м и dm =50 мкм отличаются малыми размерами и высокой интенсивностью перемешивания. Когда dn возрастает в вышеуказанном диапазоне, вслед за перемещением слоя

смешения воды место распада плёночного течения факела распыла (рис. 9, б, 10, б, 11, б, 12, б) перемещается ближе к срезу сопла внешней центробежной форсунки. Данные на рис. 9, 10 соответствуют одним и тем же значениям параметров dn =5 м и

dm =50 мкм на различных промежутках времени от начала расчёта At = 15,0 мс (рис. 9) и At = 25,0 мс (рис. 10). Изменения в структуре течения у корня факела распыла и на выходе из гидравлического тракта смесителя на рис. 9, а и 10, а, а также в структуре свободного факела распыла (рис. 9, б и 10, б) свидетельствуют о том, что картина тече-

ния факела распыла в расчётной области за время At = 25,0 мс не установилась. Данные на рис. 11, 12 показаны при одних и тех же значениях параметров =5000 м и

ёт =50 мкм, но для различных промежутков времени от начала расчёта At = 24,0 мс

(рис. 11) и At = 40,0 мс (рис. 12). Изменения в структуре течения у корня факела распыла и на выходе из гидравлического тракта смесителя на рис. 11, а и 12, а, а также в структуре свободного факела распыла (рис. 11, б и 12, б) свидетельствуют о том, что картина течения в расчётной области за время At = 40,0 мс не установилась.

Результаты расчёта наиболее близки к картине течения в гидравлических испытаниях при параметрах ёт =50 мкм и ёп =5000 м, входящих в формулы (18) и (19). На

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

ь

Рис. 13. Сравнение результатов расчёта и фотографирования при гидравлических испытаниях, объёмная доля воды - 0,1

Выводы

При моделировании картины течения факела распыла в смесителе ЖРДМТ тягой 10-15 Н с помощью двухскоростной модели течения двухфазной среды и методики определения коэффициента сопротивления Са, изложенной в настоящей работе, удаётся получать различные картины течения факела распыла воды. Исследовано влияние изменения характерных размеров в диапазонах ёт от 20 до 150 мкм и ёп от 20 мкм до

5000 м. Путём подбора значений ёт = 50 мкм и ёп = 5000 м получено хорошее совпадение расчётной картины течения с результатами гидравлического испытания. В процессе расчётных работ подтверждена работоспособность и возможность использования алгоритмов вычисления коэффициента сопротивления С0, показанных на рис. 2, 3.

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

1. Продолжить исследование влияния характерных размеров йт и йп на параметры течения в факеле распыла.

2. Получить картину течения факела распыла, близкую к установившейся, при увеличении времени расчёта до At = 50 мс и более.

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

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

1. Андреев Ю.З. Исследование зависимостей характеристик ЖРДМТ тягой 50...400 Н на топливе АТ+НДМГ от основных параметров двухкомпонентной соосной центробежной форсунки и струйных форсунок завесы. Дис. ... канд. технических наук. Нижняя Салда, 2004. 181 с.

2. Сёмкин Е.В. Исследование процессов истечения жидкости из малорасходных центробежных форсунок жидкостного ракетного двигателя малой тяги // Вестник Самарского университета. Аэрокосмическая техника, технологии и машиностроение. 2016. Т. 15, № 4. С. 150-161. DOI: 10.18287/2541-7533-2016-15-4-150-161

3. Туманова Е. Использование CAE системы Flow Vision для исследования взаимодействия потоков жидкости в центробежно-струйной форсунке // САПР и графика. 2005. № 9. https://sapr.ru/article/14574

4. Сербин С.И., Вилкул С.В. Моделирование процесса распыливания жидкого топлива центробежной форсункой // Авиационно-космическая техника и технология. 2008. № 7 (54). С. 95-98.

5. ANSYS CFX-Solver, Release 10.0: Theory. ANSYS Europe Ltd, 2005. 266 p.

6. Стернин Л.Е., Шрайбер А.А. Многофазные течения газа с частицами. М.: Машиностроение, 1994. 320 с.

NUMERICAL SIMULATION OF LOW LIQUID FLOW IN A MIXER CONSISTING OF TWO CENTRIFUGAL SWIRLERS

© 2018

E. V. Semkin Leader of Design Department Team;

Research and Development Institute of Mechanical Engineering,

Nizhnyaya Salda, Russian Federation;

niimash@list.ru

The paper describes a procedure of numerical simulation of operating fluid (water) discharge into the air environment through the hydraulic paths of a mixer composed of two low-flow coaxial swirl injectors. The method is based on a two-velocity model of a two-phase liquid flow and determining the interface interaction resistance coefficient versus the Reynolds number in the course of solution. The Reynolds number is calculated from the relative velocity of the liquid components forming the two-phase flow. The paper investigates variations in the mixer spray pattern versus two characteristic dimensions reckoned among the key parameters to calculate the coefficient of interface interaction resistance. Algorithms of calculating the coefficient of interface interaction resistance are proposed. The results of modeling liquid mixing and film flow breakdown within and beyond the mixer hydraulic paths for different values of the characteristic dimensions are presented. It is shown that we can achieve the conformity of the calculation results obtained by using the proposed method with the cold flow data by selecting the values of characteristic dimensions with reference to which interface interaction resistance coefficients are determined. Further works are projected.

Numerical simulation method; mixer; swirl injector; flow simulation; low-thrust liquid rocket engine; hydraulic path.

Citation: Semkin E.V. Numerical simulation of low liquid flow in a mixer consisting of two centrifugal swirlers. Vestnik of Samara University. Aerospace and Mechanical Engineering. 2018. V. 17, no. 4. P. 141-154. DOI: 10.18287/25417533-2018-17-4-141-154

References

1. Andreev Yu.Z. Issledovanie zavisimostey kharakteristik ZhRDMT tyagoy 50...400 N na toplive AT+NDMG ot osnovnykhparametrov dvukhkomponentnoy soosnoy tsentrobezhnoy forsunki i struynykh forsunok zavesy. Dis. ... kand. tekhn. nauk [Dependence of 50...400 N NTO+UDMH liquid rocket engine characteristics on the basic parameters of a bipropellant coaxial swirl injector and jet injectors of the cooling film. Cand. Sci. (Eng.) dissertation]. Nyzhnyaya Salda, 2004. 181 p.

2. Semkin Ye.V. Research of liquid outflow from low-emission swirl injectors of a low-thrust liquid rocket engine. Vestnik of Samara University. Aerospace and Mechanical Engineering. 2016. V. 15, no. 4. P. 150-161. DOI: 10.18287/2541-7533-2016-15-4-150-161. (In Russ.)

3. Tumanova E. CAE-system Flow Vision for investigation of liquid flow interaction within swirl-spray injectors. SAPR I Grafika. 2005. No. 9. Available at: https: //sapr.ru/article/14574. (In Russ.)

4. Serbin S.I., Vilkul S.V. Liquid fuel atomization modeling by pressure-swirl atomizer. Aerospace Technic and Technology. 2008. No. 7 (54). P. 95-98. (In Russ.)

5. ANSYS CFX-Solver, Release 10.0: Theory. ANSYS Europe Ltd, 2005. 266 p.

6. Sternin L.E., Shrayber A.A. Mnogofaznye techeniya gaza s chastitsami [Multiphase gas particle flows]. Moscow: Mashinostroenie Publ., 1994. 320 p.