7 / 18
ISSN 0236-3941. Вестник МГТУ им. Н.Э. Баумана. Сер. «Машиностроение». 2016. № 3
43
В оригинальной версии модели [17, 18] диаметр лигамента для ко-
ротковолнового режима неустойчивости определяется аналогично
длинноволновому по формуле
16 .
l
s
h
d
K
В более поздней версии, реализованной в ANSYS FLUENT, для
коротковолновой неустойчивости диаметр лигамента связывается ко-
эффициентом пропорциональности
L
C
с длиной волны неустойчиво-
сти
2
.
s
K
Этот коэффициент пропорциональности назван постоян-
ной лигамента [16]. Подход к моделированию изменяется вследствие
того, что при коротковолновой форме неустойчивости пленки срыв ли-
гаментов и капель происходит с гребешков, образовывающихся на
волнах, идущих по пленке.
Таким образом, для коротковолновой неустойчивости диаметр ли-
гамента можно определить с помощью следующей формулы [16]:
2 ,
L
l
s
C d
K
где
L
C
— постоянная лигамента (исследуемый параметр модели).
Диаметр результирующих капель определяют по зависимости
1/6
0
1,88 1 3Oh ,
l
d
d
где
l
d
— диаметр лигамента;
Oh
— число Онезорге.
Диаметр результирующих капель затем используется как харак-
терный линейный размер в распределении Розина — Раммлера, с по-
мощью которого определяется диаметр капель на выходе:
1
exp
,
q
q
q
qd
d
f d
D
D
где
q
— параметр распределения, равный по умолчанию
2
0
.
ln2
q
d D
Таким образом, для получения распределения капель по диамет-
ру и скорости с помощью используемой модели первичного распада
необходимо задать следующие входные параметры: физические
свойства жидкой и газообразной фаз, угол распыливания и скорость
жидкой пленки на выходе из сопла форсунки. Угол распыла и ско-
рость пленки могут быть получены из гидравлических проливок
форсунки или путем численного моделирования течения жидкости в
форсунке. Приведенные уравнения позволяют провести численное