Исследование течения газа в канале при направленном движении потока пара металла методом пробной частицы

Тип работы:
Реферат
Предмет:
Физика


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 81. 29. 14
Н. К. Никулин, О. А. Шемарова
ИССЛЕДОВАНИЕ ТЕЧЕНИЯ ГАЗА В КАНАЛЕ ПРИ НАПРАВЛЕННОМ ДВИЖЕНИИ ПОТОКА ПАРА МЕТАЛЛА МЕТОДОМ ПРОБНОЙ ЧАСТИЦЫ
Для исследования течения газа в канале через поток пара металла предложена математическая модель, сочетающая в себе два принципиально разных подхода, основанных на статистическом методе для моделирования течения газа и на газодинамических методах сплошной среды — для металлического пара.
E-mail: nkn@bmstu. ru
Ключевые слова: вакуум, канал, металлический пар, метод, плотность потока частиц, разреженный газ, расчет, скольжение, скорость, статистическое моделирование, течение.
Объектом исследования является течение газа в канале, где движется поток металлического пара высокой концентрации, на переходном или даже вязкостном режиме. Относительно большая концентрация металлического пара позволяет описать его течение ламинарными законами, в то время как для моделирования течения разреженного газа применяется метод Монте-Карло.
Статистический подход позволяет создать относительно простую математическую модель, с помощью которой можно довольно точно описать процесс течения газа через поток пара металла. Достоинством данного подхода является относительно легко реализуемая возможность введения новых факторов и оценки влияния сопутствующих процессов- в частности, данная методика позволяет учесть поглощение молекул газа металлическим паром и влияние вектора скорости потока пара на газ.
При описании течения металлического пара необходимо принять во внимание, что при уменьшении давления в канале возникает скачок скорости на поверхности стенки, т. е. скорость потока газа в отличие от ламинарного течения, при котором скорость на поверхности стенки считается равной нулю, имеет конечное значение, отличное от нуля. Такое течение газа называют течением со скольжением. Физически такую задачу можно рассматривать как предельный случай течения Куэтта. Математически задача состоит в решении уравнения Больц-мана с соответствующими граничными условиями на стенке.
В качестве расчетной схемы (рис. 1) рассмотрим течение газа в прямом тонком цилиндрическом канале (R ^ L), в котором движется металлический пар. Направление движения пара задается следующим образом: прямой поток (диффузия) — поток пара вместе с газом движется в направлении откачки- встречный поток — противодиффузия.
Рис. 1. Расчетная схема
Геометрическая структура состоит из поверхностей входа, боковой и выхода. Вся поверхность трубы разбита на тонкие кольца шириной Дж. В качестве легкоплавкого металла рассматривается кадмий, исследуемый газ — азот, причем давление пара кадмия намного больше, чем давление азота /pCd ~ 10−5.
Для построения относительно простой (легко модифицируемой и не требующей большой вычислительной мощности) математической модели a priori вводятся следующие допущения: распределение скоростей соответствует закону Максвелла- при взаимодействии молекул газа и стенки коэффициент аккомодации равен единице- соударение молекул рассматривается как упругий удар жестких сфер- влиянием потенциальных полей можно пренебречь- распределение скорости потока пара металла в сечении представляет собой параболический профиль с поправкой на скорость скольжения- влиянием газа на движение пара металла можно пренебречь- для учета сорбирующих свойств пара кадмия и поверхности трубы вводятся коэффициенты K и в, задаваемые как исходные данные.
Математическая модель. При статистическом моделировании течения газа в системе исследуют большое число траекторий движения молекул от момента старта с входного сечения системы до момента
возвращения к сечению входа или выхода либо до момента поглощения молекул поверхностью трубы или паром кадмия. Вероятность перехода молекул через вакуумную систему (коэффициент проводимости) будет определяться как отношение
р = N1,
N '-
где N1 — число молекул, попавших в выходное сечение- N — общее число рассматриваемых молекул.
Аналогично определяются другие молекулярные характеристики. Процедура вычисления конкретных молекулярных характеристик методом Монте-Карло (методом пробной частицы) приведена на рис. 2 в виде структурной схемы для случая вычисления коэффициента проводимости.
Рассматривая элементарное событие и его вероятность, не составляет труда вывести выражения для определения координат старта, траектории частицы и координат точек пересечения траектории с поверхностью [1−3]. Подробнее стоит рассмотреть взаимодействие молекул газа с паром кадмия.
Необходимо найти расстояние, которое пройдет молекула до столкновения. Средняя длина свободного пути молекулы [4] определяется по формуле
Л =
ипа2
где и — число молекул в единице объема- а — эффективное сечение столкновения.
Пусть концентрация кадмия и0, тогда на длине пути & amp-х с поперечным сечением Б содержится и0па2 молекул, сумма поперечных сечений которых? Б = и0па2Б& amp-х, следовательно, вероятность попадания молекулы азота в одну из частиц кадмия составляет
& amp-Б 2
& amp-Р = -- = и0па ах. Б
Далее нужно определить число молекул, которые столкнулись на участке от I до I+& amp-1. Пусть число молекул газа, которые пролетели расстояние I без столкновения, равно N, а число молекул, пролетевших расстояние I + & amp-/, составляет N — dN. Относительное число убывших частиц из потока нерассеянных частиц равно
dN 2
— = -Попа & amp-1 = - л. (1)
После интегрирования, учитывая, что число падающих частиц при начальной координате х = 0 равно N0, получаем
¦1, 2 N = N0 ¦ е- а = N0 ¦ е~1пожа. (2)
Рис. 2. Структурная схема определения вероятности перехода молекул через элемент вакуумной системы методом Монте-Карло
Формулой (2) определяется число молекул, прошедших путь I без столкновений, т. е. вероятность прохождения пути I без столкновений составляет
N (I)
P (l) =
N0
= e л.
Чтобы получить функцию распределения р (1) (плотность вероятности), надо определить вероятность столкновения частицы на участке от I до I + сИ:
I ы -1 с
= = = ^^х = 1е-а
W ы0 Ыо Ыо Л '-
1
СР (/) = = = е А
Л
Итак, плотность вероятности (или функция распределения вероятности) столкновения частицы
1 гст
р (гст) = = е- а, Л
kT
2
+ Д + MN2
где
Л =
(^N2
ПР I --,, , _ ,
2 J V Med
— средняя длина свободного пробега [1].
Расстояние гст, пролетев которое молекула газа столкнется с частицей пара, определяется из выражения
гст
J p (x)dx = о
гст
Г 1 -..
=e аdx =? ^ гст = -Лln (1 — ?),
о
где? — случайное число, равномерно распределенное на участке от 0 до 1, генерируемое ДСЧ.
Таким образом, если расстояние, пройденное молекулой газа, меньше гст, то столкновения не произойдет, если больше — молекула азота столкнется с кадмием, пройдя расстояние гст. В последнем случае вводится поправка на массовую скорость, после чего уточняется траектория молекулы и координаты столкновения.
Давление пара кадмия относительно высокое, порядка 1 Торр, что соответствует переходному режиму течения. Это позволяет описать течение потока пара металла с помощью законов ламинарного течения жидкости. За основу взято течение Пуазейля [2].
Скорость скольжения, как отмечалось ранее, определяется по методике, основанной на рассмотрении предельного случая течения Ку-этта. Считаем, что на расстоянии нескольких длин свободного пробега любые особенности функции распределения, обусловленные влия-
нием стенки, полностью исчезают и справедливо обычное гидродинамическое решение.
Чтобы не усложнять задачу, стоит принять простейшую модель линеаризованного оператора столкновений, а именно линеаризованную модель БГК [5], а также учесть, что в рассматриваемой задаче состояние газа меняется только в направлении х, гидродинамическая скорость имеет только ?-компоненту, а изменениями плотности и температуры можно пренебречь. Таким образом, необходимо решить уравнение Больцмана:
C дф, С) _ v
Cx-^- _ v
дх
п-3/2Cz C'-ze-C ip (x, C'-) d3C- (х, C)
где C — безразмерная собственная скорость- & lt- - потенциал межмолекулярного взаимодействия- г — радиус-вектор молекулы- V — показатель отталкивания межмолекулярного потенциала.
Введем граничные условия на диффузное отражение без изменения температуры, означающее, что в результате столкновения со стенкой молекулы приобретают максвелловское распределение скоростей с нулевой гидродинамической скоростью & lt-(0, С) = 0 при Сх & gt- 0.
Граничные условия задают лишь половину функции распределения. Чтобы завершить постановку задачи, остается потребовать, чтобы функция & lt- менялась при хто по линейному закону, т. е. чтобы гидродинамическая скорость на больших расстояниях от стенки росла линейно относительно расстояния: & lt- ~ vx (х ^ то).
В результате решения данной задачи о течении со скольжением получили соотношение для безразмерной гидродинамической скорости:
сю сю
V,(х) = I е-сХд (х, Сх)& lt-1Сх = х + а+ + ^ Л (Л)е-х/А^Л. (3)
-с 0
Первый член уравнения (3) представляет собой как раз гидродинамическое решение и его роль важна вдали от стенки. Последний член исчезает на расстоянии от стенки порядка нескольких длин свободного пробега, и обычно вблизи стенки он тоже мал. Таким образом, поправка определяется вторым слагаемым.
В потоках разреженных газов эффект скольжения может стать значительным. Его, по крайней мере, можно частично учесть, заменив обычное граничное условие газовой динамиких (0) = 0 граничным условием вида
^(0) = а+ & quot-ТТ ,
dx
x=0
или в размерной форме
vz (0) = а+Л dVz dx
ж=0
где использовано соотношение между длиной свободного пробега Л и коэффициентом вязкости п, вычисленным на основе модели БГК [5]- р — плотность массы. Численное значение коэффициента а+ составляет 1,0161.
В соответствии с принятыми допущениями течение потока пара в трубе может быть описано уравнением [2]
dvM = Др (4)
dp 2Ln '-
где vM — массовая скорость (скорость потока пара кадмия) — р — расстояние от оси трубы.
После интегрирования уравнения (4) с учетом граничного условия (скорость пара на самой поверхности трубки равна скорости скольжения)
dvM
Vm® = a+L-p |P=R
получаем зависимость
VM (p) = 4Ln (R2 — p2) + a+LЩAp-
Максимальная скорость в сечении и скорость на оси трубы определяются как массовый расход газа через единицу поперечного сечения
АР С& gt-2 (5)
VM max = R • (5)
4Ln
По закону Пуазейля расход пара в трубе определяется по формуле
кАр 4 б
V = --- R4, следовательно, поток пара в трубе составляет 8Ln
п пАр г?4 (6)
Q = р° (6)
Пользуясь зависимостями (5) и (6), можно переписать уравнение распределения скоростей потока пара металла в зависимости от расстояния до оси трубы в удобной форме:
Р, r 2Vm max
Vm (Р) = VMma^ 1 — R^) + а+Lp-
R2
где р = /у2 + г2 — расстояние от оси трубы- Я и Ь — радиус и длина трубы.
Теперь можно найти массовую скорость в произвольной точке сечения. Тепловая скорость определяется из выражения для функции
распределения Максвелла модуля скорости в безразмерной форме
4
v© = e c2c2dc, /п
где c = ---безразмерная скорость молекулы- VH = 129
vh v Med
наиболее вероятная скорость.
Тепловая скорость частицы кадмия vH определяется из выражений:
c
4
i = J -^=e-c2c2dc = f (c) — 0
c = f-1 (i) ^ VH = c (i) ¦ 129 M
V Med
Тепловую скорость молекулы азота можно найти аналогично.
Скорость атома кадмия складывается из тепловой (vH) и массовой (vM) составляющих скорости в данной точке сечения. Теперь, когда известны обе составляющие скорости, необходимо определить траекторию. Из треугольников скоростей определяют новые координаты 9'- и ф'- для случаев диффузии и противодиффузии (рис. 3):
-V -" I -«
V = VH + VM- (v'-)2 = vH + V2M — 2vhVm cos (2 + 9) ^
^ v'- =vH + vM — 2vhVm sin 9 —
., vh sin 9 ± vm u,. vh sin 9 ± vm sin 9 =---=& gt-• 9 = arcsin----
^ /п, Д vh sin 9 ± vm /п д
tg Ь -Ф'- =-tg (o -Ф) ^
2 vH sin 9 2
п. fvH sin в ± Vm (n \
^ V = o — arct^ --- tg — - Ф)
2 vH sin в V 2 /у
Выражения, где встречается & quot-±"- или & quot-+"-, применяются в случае диффузии (рис. 3, а), выражения, где встречается & quot--"-, — в случае противодиффузии (рис. 3, б). Под диффузией в данном случае понимается течение разреженного газа и пара металла в направлении откачки, под противодиффузией — течение газа в направлении откачки через встречный поток металлического пара.
Теперь можно определить скорость и траекторию молекулы азота после столкновения, воспользовавшись уравнениями сохранения импульса и энергии в проекциях на оси X и Z:
mN2 v1 sin cos + mCdv2 sin #2 cos =
= mN2 v3 sin #3 cosз + mCdv4 sin #4 cos
mN2 v1 cos + mCdv2 cos #2 = mN2 v3 cos #3 + mCdv4 cos #4-
2 2 2 2 v! X + mCdv2x = vgx + mCdv4x.
2 2 2 2 '-
2 2 2 2 vb + mCdv2z = vjz + mCdv4z
2 2 2 2 '- где vix = vi sin 9i cos viz = vi cos v1 и v3 — скорости молекул азота до и после столкновения с кадмием- v2 и v4 — скорости молекул кадмия до и после столкновения.
Решив данную систему, получим скорость v3 и траекторию молекулы азота после столкновения. Траектория движения молекулы азота после соударения с частицей кадмия определяется углами и #3.
Частицы кадмия после соударения с молекулами азота не рассматриваются, так как азот не влияет на статистическую картину распределения пара кадмия (давление пара кадмия много выше, чем давление газа).
Поглощающие свойства поверхности трубы учитываются c помощью коэффициента захвата в & gt- 0, задаваемого как исходные данные. С использованием ДСЧ генерируется случайное число р. Если в точке столкновения с поверхностью р ^ в, то молекула считается захваченной и фиксируется счетчикомпогл. тр, а при р & gt- в — отраженной. Аналогично с помощью коэффициента K поглощения газа кадмием учитываются и сорбирующие свойства пара кадмия.
Численный эксперимент. Для построения зависимостей плотностей потока падающих ^пад) и поглощенных (v^m) молекул азота на поверхность трубы по ее длине вся поверхность разбивается на кольца, в каждом из которых фиксируется число падающих (^пад) и поглощенных (^погл) молекул:
N
_ 1'- шад
пад = NF'-
_ ^?погл
vin°ra = & quot-NFT'-
где Fi — площадь поверхности ?-го кольца.
Для расчета модели в среде Matlab 7.9.0 была написана программа, с помощью которой проведен численный эксперимент.
Распределения плотностей потоков падающих молекул газа по относительной длине трубы, построенные для обоих случаев (диффузии и противодиффузии) при различных значениях коэффициента прилипания трубы в, представлены на рис. 4.
vm
а
б
Рис. 3. Расчетная схема для уточнения траектории с поправкой на массовую скорость ум в случае диффузии (а) и противодиффузии (б)
-2
V, М
10 000 с
0 10 20 30 40 50 60 70 80 90 LID
Рис. 4. Изменение плотности потока падающих частиц по длине трубы при различных значениях коэффициента прилипания трубы в в случае диффузии и противодиффузии
V, м
9000
8000
? = 0
-3 3
7000
6000
5000
4000
3000
2000
1000
0
2 4
6
8 LID
Рис. 5. Изменение плотности потока падающих частиц по длине трубы при различных значениях коэффициента сорбции К для случая диффузии
Изменение плотности потока падающих частиц по относительной длине трубы при разных значениях коэффициента сорбции К, характеризующего поглощающую способность металлического пара, в случае диффузии приведено на рис. 5. Очевидно, что чем выше сорбирующие свойства пара, тем ниже плотность потока молекул газа.
Влияние на вероятность перехода молекул направления и значения потока пара кадмия ^ показано на рис. 6. Эти данные также позволяют оценить влияние потока кадмия на проводимость системы. В случае, когда потоки газа и пара сонаправлены, наблюдается дополнительный откачной эффект, вероятность перехода молекул увеличивается (рис. 6, а), в случае противодиффузии — уменьшается (рис. 6, б).
Выводы. Из анализа полученных графических зависимостей для рассмотренной математической модели следует, что в случае сонапра-вленного движения газа и пара металла наблюдается увеличение числа частиц в направлении потока диффундирующего газа, что ведет к увеличению вероятности перехода (8% при потоке пара кадмия 1 м3-Па/с и до 40% при больших потоках). Аналогично в случае противодиффузии это ведет к уменьшению вероятности перехода (4% при потоке 10−3 м3-Па/с, с увеличением потока до 1 м3-Па/с падает на 70% и до 99% при больших потоках).
После уточнения значений коэффициентов захвата поверхностью трубы и сорбции паром данная модель может использоваться для определения параметров течения газа в канале (распределение давления по длине канала), определения молекулярных характеристик (например,
Рис. 6. Зависимость вероятности перехода P молекул газа через трубу от потока пара кадмия Q для сонаправленного (а) и встречного (б) движения потоков газа и пара кадмия в случае диффузии (а) и противодиффузии (б)
коэффициента проводимости), в котором присутствуют пары легкоплавких металлов.
СПИСОК ЛИТЕРАТУРЫ
1. Вакуумная техника: Справочник / К. Е. Демихов, Ю. В. Панфилов, Н. К. Никулин и др. — М.: Машиностроение, 2009. — 590 с.
2. Лойцянский Н. Г. Механика жидкости и газа: Учеб. для вузов. — М.: Высш. шк., 2003. — 840 с.
3. СаксаганскийГ. Л. Молекулярные потоки в сложных вакуумных структурах. — М.: Атомиздат, 1980. — 216 с.
4. Яворский Б. М., ДетлафА. А. Справочник по физике. — М.: Наука, 1979. — 942 с.
5. Ферцигер Д ж., Капер Г. Математическая теория процессов переноса в газах. — М.: Мир, 1976.
Статья поступила в редакцию 15. 06. 2011

ПоказатьСвернуть
Заполнить форму текущей работой