Файл: Стернин Л.Е. Основы газодинамики двухфазных течений в соплах.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 27.06.2024

Просмотров: 161

Скачиваний: 4

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

вать для широкого использования при расчетах двухфазных те­ чений в соплах с коагуляцией.

В работе [116] проведены также исследования влияния дав­ ления торможения перед соплом, размеров сопла и содержания частиц в газе на их рост и на потери удельного импульса. Полу­ чено, что с ростом давления торможения и содержания частиц в газе коагуляция интенсифицируется, что приводит к интенсивно­ му росту размеров частиц и потерь из-за их отставаний. Из-за увеличения времени пребывания частиц в сопле, например, при увеличении размеров сопла (изменении масштаба) происходит их заметное укрупнение. Однако при этом потери из-за отстава­ ния в целом уменьшаются, так как с ростом размеров сопла уве­ личивается и время релаксации. Основной процесс коагуляции происходит в околокритической части сопла, где отставания и различия в скоростях частиц наибольшие. В расширяющейся части сопла кроме отставаний резко падает и концентрация час­ тиц в единице объема, что существенно уменьшабт коагуляцию, поэтому здесь роста частиц почти не происходит.

Учет дробления при взаимодействии частиц в процессе их коагуляции [6] позволяет более правильно найти размеры час­ тиц в сопле и потери из-за отставания.

На

рис. 2.1 представлены

(при

содержании частиц

в газе

W = 3)

расходные массовые функции

распределения частиц по

размерам gr U g r {r)dr =

в начальном сечении сопла

(/) и

 

/

(пунктиром) нанесена кривая (3),

на срезе сопла

(2). Здесь же

соответствующая расчету без дробления. На рис. 2.2 изображено

изменение среднемассового

радиуса частиц

riz по соплу (х =

= х//, где I—длина сопла)

при W = 3.

Верхняя

(пунктирная)

кривая соответствует коагуляции без

учета дробления, а ниж­

няя — с дроблением. Как видно из рис. 2.2,

на

срезе сопла при

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

Па рис. 2.3 показано изменение потерь удельного импульса, обусловленных отставанием частиц в зависимости от W. Пунк­ тирная кривая соответствует коагуляции без учета дробления, а сплошная — коагуляции с учетом дробления. При большом со­ держании частиц в газе потери могут различаться в несколько раз. (Потери при отсутствии учета взаимодействия частиц меж­ ду собой находятся в пределах 1%).

78


дттм'

Рис. 2.1. Изменение функций распреде­ ления частиц в сопле

W M

 

X

\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

X

\

 

 

 

 

 

 

 

 

!

 

 

 

 

 

 

 

 

/

 

 

 

 

 

 

 

 

/

/

 

 

 

 

 

 

 

 

У

 

 

 

 

 

 

 

 

.t

 

 

 

 

 

 

 

 

—/"

 

 

 

 

 

 

 

 

 

-2,5 -2,0

-1,5 -1,0 -0,5 Lffx

О

1

Ч

6

8 W

Рис. 2.2.

Изменение

средне-

 

Рис. 2.3. Кривые зависи-

массового размера

частиц

 

мости

потерь

от

содер-

 

вдоль сопла

 

 

жания частиц в газе

79


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

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

В качестве примеров, иллюстрирующих замену полидисперсного потока жидких частиц (при наличии эффектов коагуляции и дробления) эквивалентным монодисперсным потоком, были проведены расчеты потерь удельного импульса для обоих случа­

ев при W = 3 и 6. Получено,

что при наличии коагуляции при

W = 3 потери составляют около 4,8%

при

среднемассовом диа­

метре частиц в критическом

сечении

(т.

е. где w = a ) — 3,1р.

Для монодисперсного потока с этим же размером частиц потери составляют около 5,6%• При W = 6 среднемассовый диаметр час­ тиц оказывается равным 5,3р, а фактические потери — около 9,5%. Для монодисперсного потока — около 11%. Таким обра­ зом, для приведенных примеров увеличение потерь при эквива­ лентном монодисперсном потоке с размером частиц, равным их среднемассовому размеру в критическом сечении, составляет около 15%.

Результаты серии расчетов коагуляции в соплах приведены в справочнике [115], где в сжатом виде изложена также и общая теория коагуляции.

Исследования,

проведенные недавно А. М. Подвысоцким и

А. А. Шрайбером,

взаимодействия свободно падающих крупных

капель с движущимися в горизонтальной плоскости мелкими кап­ лями привели к формуле для Фц, заметно отличающейся от фор­ мулы (2. 31),

Фц = 1—0,247 Re°'434r 0’133 (п/г})°'273,

(2. 32)

что может быть объяснено армирующим влиянием подвески (ее стабилизирующим воздействием на мишень при низких значени­ ях параметра ReBr 0’4 и воздействием, стимулирующим дробле­ ние, при высоких значениях этого параметра). Потери удельного импульса, определенные по формуле (2.32), меньше (до 1—2%) потерь, вычисленных по формуле (2.31).

В заключение отметим, что расчеты, проведенные этими же авторами показали, что учет произвольного полидисперсного со­ става осколков с произвольными начальными скоростями и тем­ пературами может привести к уменьшению потерь удельного им­ пульса в пределах 0,5%.


Г л а в а III

ЭЛЕМЕНТЫ КИНЕТИКИ КОНДЕНСАЦИИ

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

Конденсация может иметь место в том случае, когда парци­ альное давление пара одного из компонентов смеси в процессе расширения становится большим его давления насыщения. При этом, как правило, пересыщение s (отношение парциального давления к давлению насыщения над плоской поверхностью пе­ рехода), при котором конденсация становится заметной, намно­ го больше единицы. Например, для водяного пара s ä 4-^8, а для паров ртути пересыщение может превышать 2000.

Явление конденсации еще в прошлом веке исследовал Виль­ сон, который установил, что в чистом воздухе, насыщенном во­ дяным паром, при достаточно быстром расширении самопроиз­ вольно появляется облако, состоящее из очень мелких водяных капелек. При первоначальной температуре 20°С облако появля­ лось только в том случае, когда степень расширения превышала приблизительно 1,37 и когда создавалось почти восьмикратное пересыщение. Камера, в которой Вильсон изучал поведение сме­ си паров воды с различными газами, впоследствии стала назы­ ваться его именем.

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

81

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

Интересной иллюстрацией гетерогенной конденсации являет­ ся конденсация паров одного вещества на ядрах другого. Так, при наличии в потоке газа даже небольшого количества другого ве­ щества с более высокими температурой и давлением насыщения в первую очередь происходит гомогенная конденсация этого ве­ щества с образованием большого количества зародышей, на ко­ торых в дальнейшем конденсируется основной газ. При этом пе­ ресыщение газа почти или полностью отсутствует. В качестве примера можно привести конденсацию азота при наличии в нем около 0,05% водяного пара [24]. Эта конденсация оказывается практически равновесной.

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

Фундаментальный теоретический анализ конденсации в соп­ лах был сделан Осватичем в 1942 г. [156], который первым ввел в аэродинамику теорию образования зародышей. Осватич соста­ вил уравнения, описывающие возникновение и рост частиц при движении двухфазной среды в сопле с учетом передачи тепла от частиц к газу, и решил их численным методом. Полученные ре­ зультаты свидетельствуют о росте статического давления в зоне конденсации и хорошо совпадают с экспериментальными данны­ ми Бинни и Вудса. Основные уравнения, упрощенно описываю­ щие гомогенную кондѳндсацию в рост капель при движении сме­ си воздуха и паров воды в сопле аэродинамической трубы, приве­ дены в работе Малюжинца [87].

Теоретическое описание течений с конденсацией, проведенное в работах [87 и 156], основано на использовании выражения для скорости образования зародышей. Вывод такого выражения яв­ ляется одной из наиболее сложных проблем кинетики конденса­ ции. Классическая теория образования зародышей (нуклеации), основанная на жидкокапельной модели определения концентра­ ций групп молекул, была предложена Фольмером и Вебером

[164]и развита в работах Беккера и Дёринга [139], Зельдовича

[56]и Френкеля [123]. Несмотря на некоторые недостатки жидко-

82


капельной теории, теоретическое описание процесса конденсации в камере Вильсона, основанное на этой теории, согласуется с экспериментом [122].

Джёрмати [145] на основании классической теории нуклеации предложил упрощенный аналитический метод оценки место­ положения точки, в которой конденсацию можно считать закон­ ченной, и точки Вильсона (в этой точке разница между факти­ ческим количеством конденсированного вещества и таковым для равновесного процесса является максимальной; вблизи точки Вильсона имеет место максимальное пересыщение). В работе Джёрмати и Мейера [147] представлены результаты экспери­ ментального исследования движения сухого и влажного водяно­ го пара в трех плоских соплах, характеризующихся различными градиентами давления и имеющих различные размеры. Наблю­ далось выпадение конденсата и повышение давления в зоне конденсации. Изменение давления вниз по потоку от зоны кон­ денсации происходило эквидистантно изменению давления для сухого пара. При расчете процесса нуклеации по классической теории имело место хорошее соответствие экспериментальных данных с результатами расчетов при условии увеличения поверх­ ностного натяжения на 13%. Такое завышение коэффициента поверхностного натяжения для совпадения теоретических и экс­ периментальных данных делал Осватич [156], отмечая, что для быстродвижущейся среды нельзя использовать поверхностное натяжение, измеренное в состоянии покоя.

Хилл [150] рассмотрел экспериментальные данные по скоро­ сти роста зародышей после нуклеации и сделал вывод, что ско­ рость роста для малых капель не зависит от их радиуса. Это поз­ волило ему описать процесс конденсации тремя обыкновенными дифференциальными уравнениями. На основании анализа теплоII массобмена капли с внешней средой Хилл вывел уравнение для скорости роста капли и уравнение теплового баланса, из которого определяется температура капли. В эти уравнения бы­ ли введены коэффициенты конденсации ак и термической акко­ модации ß. Был проведен расчет расширения паров воды в соп­ ле в одномерном приближении. Хорошее совпадение с экспери­ ментальными данными Вегенера и Пауринга [165] имело место при ß = 1 и ак = 0,04.

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

83