CC BY
42
УДК 621.52
Исследование возможностей оптимизации откачных характеристик молекулярновязкостного вакуумного насоса методами статистического моделирования
К.Е. Демихов, Н.К. Никулин, Т.В. Дронова
Рассмотрено влияние основных параметров системы (геометрии от-качной системы и динамики течения газа) на откачные характеристики насоса. Описана методика расчета вероятности прохождения канала частицей методом Монте-Карло с учетом времени движения частицы в канале. Влияние параметров рассмотрено применительно к условиям высокого вакуума.
Ключевые слова: статистический метод, молекула, молекулярновяз-костный насос, вероятность перехода, отношение давлений, быстрота действия, давление, геометрические параметры, форма канала.
Investigation of opportunities to optimize pumping performance of molecular-viscous vacuum pump by statistical simulation methods
K.E. Demikhov, N.K. Nikulin, T.V. Dronova
The paper deals with the influence of system main parameters — the pumping system geometry and the gas flow dynamics — on the pumping characteristics of the pump. The method of calculating the probability for a particle to pass the channel by Monte Carlo method taking into account the time of a particle movement in the channel is described. The influence of the parameters is considered in relation to work in high vacuum.
Keywords: statistical method, molecule, molecular-viscous pump, transition probability, pressure ratio, operating speed, pressure, geometrical parameters, channel shape.
I I ринцип работы механических высоковакуумных насосов основан на передаче дополнительного импульса количества движения молекулам газа движущейся поверхностью при соударении в направлении откачки. Наиболее очевидное использование данного принципа положено в основу работы молекулярного вакуумного насоса (МВН), более поздняя модификация которого — молекулярно-вязкостный вакуумный насос (МВВН), использует тот же принцип, однако из-за конструктивных особенностей позволяет получить иные откачные характеристики.
ДЕМИХОВ Константин Евгеньевич
доктор технических наук, профессор, зав. кафедрой «Вакуумная и компрессорная техника физических установок»
НИКУЛИН Николай Константинович
кандидат технических наук, доцент
ДРОНОВА Татьяна Викторовна
аспирантка кафедры «Вакуумная и компрессорная техника физических установок» (МГТУ им. Н.Э. Баумана) e-mail: nkn@bmstu.ru
ВШзехэшш] выкшшпх ©аведжшй
Рис. 1. Конструктивные схемы сечений каналов МВВН:
прямоугольное; б — трапецеидальное; в — треугольное; г — полукруглое; 1 — ротор; 2
статор
Цель работы — исследование степени и диапазона влияния основных параметров системы на откачные характеристики: максимальную быстроту действия и отношение давлений.
На рисунке 1 показаны конструктивные схемы сечений каналов МВВН: прямоугольное, трапецеидальное, треугольное и полукруглое.
Поскольку в данной статье рассматривается работа МВВН в условиях высокого вакуума с молекулярным режимом течения газа, принимается допущение об отсутствии столкновений между молекулами. При этом их движение зависит только от взаимодействия с поверхностями, ограничивающими траектории. В качестве допущения также принимается модель идеального газа для описания поведения реального газа в рассматриваемых условиях.
Учитывая принятые допущения можно сделать вывод о применимости методов статистического моделирования для расчета откачных характеристик МВВН.
Наиболее зарекомендовавшим себя на практике методом статистического моделирования, ввиду своей наглядности и универсальности, является метод Монте-Карло. Метод основан на моделировании случайных величин и функ-
ций для вычисления производных характеристик системы.
Для модели идеального газа, где взаимодействие между частицами отсутствует, применим следующий алгоритм расчета. Частицы запускаются в систему по очереди, каждая последующая частица запускается после того как закончила полет предыдущая.
Применительно к исследованию возможностей оптимизации МВВН был использован иной алгоритм расчета, основное отличие которого заключается в единовременном испускании частиц с ограниченным временем пребывания в канале Л1_.
вр
Алгоритм представляет более наглядную картину движения потока идеального газа в системе, определяя его основные параметры в целом, а не каждой частицы в отдельности.
Принимаемые допущения [1]:
♦ газ — идеальный;
♦ режим течения — молекулярный;
♦ распределения тепловых скоростей движения молекул газа описывается уравнением Максвелла;
♦ отражение частиц газа от поверхностей подчиняется диффузному закону;
а
♦ поглощение газа на поверхности отсутствует;
♦ процесс рассматривается при постоянной температуре;
♦ коэффициент аккомодации равен 1;
♦ натекания извне и газовыделения отсутствуют.
Алгоритм расчета представлен на рис. 2.
Трехмерная система координат связана с поверхностью входа в канал, начало системы координат установлено в центре этой поверхности так, что ось г направлена вдоль канала, ось у — по нормали к подвижной поверхности.
Задаются геометрические параметры системы а, к, Ькан. Датчиком случайных чисел задаются значения точек старта частицы х0, у0, г0, а также направляющие косинусы движения частицы I, т, п. Вычисляются расстояния от
Рис. 2. Алгоритм расчета вероятности прохождения частицы через заданный канал методом Монте-Карло
точки старта частицы до всех поверхностей канала 51, 52... 5п и находится минимальное расстояние от точки старта до какой-либо поверхности 5т;п. На основе этого вычисляется точка столкновения с поверхностью.
Аналогичным образом происходит дальнейшее движение частицы в канале.
По координате г определяется поверхность выхода частицы из канала. При г < 0 — частица не прошла канал (вернулась в откачиваемый объем). При I > ¿кан, — частица прошла канал (перешла на сторону нагнетания). Отслеживанием времени нахождения частицы в полете I определятся число «пропавших» частиц, которое учитываются при определении вероятности прохождения канала.
Расчет четырех видов каналов выполнен для следующих исходных данных:
♦ газ — азот при температуре 293 К;
♦ длина, ширина, высота канала, угол наклона канала относительно торца у, угол наклона боковых стенок а и отношение и/¥н — переменные значения.
Изменение скорости частицы учитывается только при попадании ее на роторную поверхность. В прямом направлении откачки рассматривается 1 000 000 событий, в обратно — 10 000 000, что обеспечивает необходимую степень точности. Учет обратного потока позволяет вычислить зависимость вероятности прохождения канала £тах и ттах от отношений к/а, ь/к, и/Уя, у.
Математическая модель. 1. Разыгрывание точек старта частицы с поверхности входа в канал:
х0 = а^; у0 = к^; г0 = 0,
где а — безразмерная ширина канала; к — безразмерная высота канала; ^ — случайная величина.
2. Расчет скорости движения частицы по уравнению Максвелла:
^/2 \
^(и)=4п
М
2 пЯТ)
и ехр
Ми2
\
2ЯТ
где и — текущее значение тепловой скорости движения частицы, м/с; М — молекулярная
ВШЗЕХЭШШ] выкшшех ©аведжшй
масса; R = 8 314 — универсальная газовая постоянная, Дж/(кмоль-К); Т — температура газа, К.
3. Расчет направляющих косинусов при старте с плоской поверхности:
3.1) вычисление углов 9 и у:
9 = аге8тЛ/^, у =
3.2) направляющие косинусы [2]:
/ = 8т(9)ео8(у), m = 8т(9)8т(у), п = ео8(9).
4. Вычисление расстояния от точки старта до поверхностей канала и нахождение минимального значения этого расстояния (выбор поверхности столкновения) [2]:
Ах о + Ву0 + Сг, + Б
5 = --
Al + Bm + Cn
где А, В, С, Б — коэффициенты уравнения плоскости Ах + Ву + Сг + Б = 0.
5. Нахождение точек столкновения частицы с поверхностью:
х = /5 + х 0;
у = т5 + у о;
г = п5 + г 0.
6. Расчет времени полета частицы:
г = 5 / V.
Расчетная величина сравнивается с заданным временем нахождения частицы в канале г0. В случае если расчетное время полета частицы г больше времени г0 (г > г0), частица считается «пропавшей».
7. Проверка нахождения частицы в канале. Если г < 0 — частица вернулась в откачиваемый объем, если г > £кан — частица пролетела канал. Если 0 < г< £кан — частица находится в канале. В последнем случае п. 2—6 повторяются, при учете, что в каждом последующем случае поверхностью старта будет поверхность столкновения. Расчет ведется до тех пор, пока четко не
будет определено, пролетела частица канал или вернулась обратно в откачиваемый объем.
8. Расчет ктах и ттах канала с учетом «пропавших» частиц. Ошибку математической модели при использовании метода Монте-Карло [1] можно оценить по формуле
Л < 3д/(1- р)р / N.
Относительная величина ошибки при этом составит
5=л / р < 3^/5—РОлР^.
Последнее уравнение позволяет определить требуемое число испытаний для обеспечения заданной точности расчета
N « 9(1-р) / (р52).
При расчетах погрешность вычислений допускается в пределах 10 % — 5 = 0,1. Для обеспечения данной погрешности необходимое число испытаний составит
N « 9
1-10-
10-4 (од)2
= 10000000 « 10 млн
Исследование возможностей оптимизации.
В работе исследовано влияние на величины основных откачных характеристик насоса kmax и xmax следующих параметров:
♦ отношения высоты к ширине канала h/a;
♦ длины к высоте L/h;
♦ отношения окружной и наивероятной скорости движения частицы Ц/¥н.
Зависимости (для прямоугольного сечения канала a = 4 мм, h = 2 мм) представлены на рис. 3—17.
Выводы
Разработанный алгоритм расчета позволяет смоделировать исходные параметры откачной системы, отвечающей заданным откачным характеристикам, а также исследовать влияние динамических и геометрических параметров насоса на его откачные характеристики.
Проведенное исследование выявило характер зависимости параметров системы, что открывает возможности его дальнейшего использования для целей оптимизации.
0,60 0,50 0,40 0,30 0,20 0,10
0,00
......... --------
i'" ^.
/
¿i
-0,10
o,io
о,зо
0,50
0,70
0,90
1,10
1,30
1,50
Рис. 3. Зависимость ктах от отношения Ц/Уп; длина канала I = 40 мм, угол наклона канала к торцу входного сечения у = 25
— •• — — МВВН — прямоугольник; — • — • — — МВВН — трапеция;
—— — МВВН — полукруг; ------ — МВВН — треугольник
0,01 0,01 0,01 0,00 0,00 о,оо
-------
______
у'
л ¿Я'
-одо
0,10
0,30
0,50
0,70
0,90
1,10
1,30
1,50
Рис. 4. Зависимость Sн от отношения Ц /V,; длина канала I = 40 мм, угол наклона канала к торцу входного сечения у = 25
— •• — — МВВН — прямоугольник; — • — • — — МВВН — трапеция;
—— — МВВН — полукруг; ------ — МВВН — треугольник
600 500 400
зоо
100
/
✓
ш
шг У* У Г *** '
***
Рис. 5. Зависимость т от отношения U /V; длина
max р ' и' ^
канала l = 40 мм, угол наклона канала к торцу входного сечения у = 25
— •• — — МВВН — прямоугольник; — • — • — — МВВН — трапеция;
—— — МВВН — полукруг; ------ — МВВН — треугольник
Рис. 6. Зависимость ктх от отношения Хкан /Акан; окружная скорость и = 300 м/с, угол наклона канала к торцу входного сечения у = 25
— •• — — МВВН — прямоугольник; — • — • — — МВВН — трапеция;
—— — МВВН — полукруг; ------ — МВВН — треугольник
Рис. 7. Зависимость £н от отношения Хкан/йкан; окружная скорость Ц=300 м/с, угол наклона канала к торцу входного сечения у = 25
— •• — — МВВН — прямоугольник; — • — • — — МВВН — трапеция;
—— — МВВН — полукруг; ------ — МВВН — треугольник
юооо
юоо
ю
ti'S ¿ш- /у х"
20 зо ад so
60
Рис. 8. Зависимость ттах от отношения Хкан/кшт; окружная скорость Ц=300 м/с, угол наклона канала к торцу входного сечения у = 25
— •• — — МВВН — прямоугольник; — • — • — — МВВН — трапеция;
—— — МВВН — полукруг; ------ — МВВН — треугольник
ВШзехэшш] выкшшпх ©аведжшй
------
/ /
/
Рис. 9. Зависимость ктах от отношения высоты канала к ширине канала Акан/акан; окружная скорость и = 300 м/с, длина канала I = 40 мм; угол наклона канала к торцу входного сечения у = 25
--— МВВН — прямоугольник;
------ — МВВН — трапеция
Рис. 12. Зависимость k от угла у наклона канала
max J '
к торцу входа; окружная скорость U = 300 м/с, длина канала l = 40 мм:
— •• — — МВВН — прямоугольник; — • — • — — МВВН — трапеция;
—— — МВВН — полукруг; ------ — МВВН — треугольник
Рис. 10. Зависимость Sн от отношения высоты канала к ширине канала ккан/акан; окружная скорость и = 300 м/с, длина канала I = 40 мм; угол наклона канала к торцу входного сечения у = 25
------ — МВВН — прямоугольник;
--— МВВН — трапеция
Рис. 11. Зависимость ттах от отношения высоты канала к ширине канала к /а ; окружная скорость
•г кан' кан' г ^ г
и = 300 м/с, длина канала I = 40 мм; угол наклона канала к торцу входного сечения у = 25
------ — МВВН — прямоугольник;
--— МВВН — трапеция
0,01 0,01 0,01 0,01 0,01 0,01 0,00 0,00 0,00 0,00 о.оо
,—
** __■ ъ
\
t У \
* f
/
St V S
/
10'
20'
30'
ДО'
S0'
60'
70'
Рис. 13. Зависимость Sн от угла у наклона канала к торцу входа; окружная скорость и = 300 м/с, длина канала I = 40 мм:
— •• — — МВВН — прямоугольник; — • — • — — МВВН — трапеция;
—— — МВВН — полукруг; ------ — МВВН — треугольник
2S0
200
150
sf ' "4 \ 4. \
у / / / 7 / vs. N \\
/ / / /
t / -a—* //
10'
20'
30'
40'
50'
60'
70'
Рис. 14. Зависимость ттах от угла у наклона канала к торцу входа; окружная скорость и = 300 м/с, длина канала I = 40 мм:
— •• — — МВВН — прямоугольник; — • — • — — МВВН — трапеция;
—— — МВВН — полукруг; ------ — МВВН — треугольник
Рис. 15. Зависимость kmax от угла а' наклона канала к торцу входа; окружная скорость и = 300 м/с, длина канала ! = 40 мм:
— •• — — МВВН — прямоугольник; — • — • — — МВВН — трапеция;
—— — МВВН — полукруг; ------ — МВВН — треугольник
Рис. 16. Зависимость Sн от угла а' наклона канала к торцу входа; окружная скорость и = 300 м/с, длина канала ! = 40 мм:
—— — МВВН — прямоугольник; — • — • — — МВВН — трапеция;
— •• — — МВВН — полукруг; ----- — МВВН — треугольник
Длина канала практически не влияет на быстроту действия насоса, но существенно влияет
Рис. 17. Зависимость т от угла а' наклона канала
max
к торцу входа; окружная скорость U = 300 м/с, длина канала l = 40 мм:
—— — МВВН — прямоугольник; — • — • — — МВВН — трапеция;
— •• — — МВВН — полукруг; ----- — МВВН — треугольник
на увеличение максимального отношения давлений.
Уменьшение высоты канала по отношению к его ширине приводит к увеличению £тах.
Увеличение окружной скорости канала (или увеличение массы откачиваемого газа) увеличивает ^ах и ^тах.
В результате установлено, что оптимальное отношение высоты канала к его ширине составляет 1/2.
Литература
1. Вакуумная техника: Справочник / К.Е. Демихов, Ю.В. Панфилов, Н.К. Никулин и др.: Под общ. ред. К.Е. Де-михова, Ю.В. Панфилова. М.: Машиностроение, 2009. 590 с.
2. Корн Г., Корн Т. Справочник по математике. М.: Наука, 1977. 832 с.
Статья поступила в редакцию 07.08.2012
