08.ЛекцДеулРасчПровМетодомМК

Download Report

Transcript 08.ЛекцДеулРасчПровМетодомМК

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Основы Вакуумной Техники
проф.д.т.н. Деулин Е.А
Лекция 8
Расчёт проводимости трубопроводов
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Вопросы к билетам по представленной лекции:
1.-Формулы проводимости цилиндрического трубопровода при различных
режимах течения газа
2.-Связь проводимости трубопровода со степенью вакуума.
3.- Проводимость диафрагмы и проводимость трубопровода сложной формы
4.- Физические и технические основы расчёта проводимости методом
Монте- Карло
5.-Методика расчёта проводимости сложного трубопровода методом МонтеКарло
6.-Преимущества расчёта проводимости методом Монте- Карло перед
«геометрическими» методами расчёта по формулам проводимости
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Основное уравнение ВакТехники :
1
S0

1
U

1
SH
или
S0
SHU
SH  U
связывает параметры трёх основных компонентов вакуумной системы:
быстроту действия насоса, проводимость трубопровода и быстроту откачки
реципиента, поэтому расчёт проводимости трубопроводов актуален при
расчёте вакуумной системы
1 – насос;
2 – вакуумопровод;
3–реципиент
откачиваемый объём ).
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Проводимость параллельно соединенных трубопроводов:
UΣ=U1+U2+…Ui ;
Проводимость последовательно соединенных трубопроводов:
1
1

U
U1

1
 ... 
U2
1
Ui
 U 
1
1
1

U1
 ... 
U2
1
;
Ui
Проводимость отверстия (диафрагмы):
U
 м3
 117  S 
 с
U Д  3 6 ,4
Т
М А

;

S – площадь диафрагмы [м2];
 м 
 л
V1  117 

11
,
7

2
 с  см 2

с  м 
3
 Vi  A

;

U  36 , 4  S 
3
T м 

;
M  с 
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Формулы для расчёта проводимости цилиндрического
трубопровода
1) Вязкостный (ламинарный) режим:
d – диаметр вакуумопровода [м],
l – длина вакуумопровода [м],
P – давление [Па].
Для воздуха при t=20C => d>10·l
3
P1  P2  м 

;
l
2
 с 
d
U B  1360
4
2) Молекулярно-вязкостный режим:
d
U МB  1360
l
4
P1  P2
2
1  1, 95 d
P1  P2
 м3 
2
 121
;
P1  P2  с 
l
1  2,4 d
2
d
3
Для воздуха при t=20C
=> d>10·l
3) Молекулярный режим:
U М  38 ,1
d
l
3
3
T м 

;
M  с 
UМ
3
3
d м 
 121

;
l  с 
Для воздуха при
t=20C => d>10·l
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Проводимость длинного цилиндрического
вакуумопровода как функция геометрических
параметров и давления (режима течения газа)
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Проводимость вакуумо- провода, учитывает режим течения, характер явления
переноса (вязкости) и, тем самым, степень вакуума
Для трубы:
1) Вязкостный (ламинарный) режим (низкий вакуум):
U B  1360
3
P1  P2  м 

;
l
2
 с 
d
4
d – диаметр вакуумопровода
[м],
l – длина вакуумопровода [м],
P – давление [Па].
Для воздуха при t=20C => l
>10· d
2) Молекулярно-вязкостный режим (средний вакуум):
U МB  1360
d
4
P1  P2
l
2
1  1, 95 d
P1  P2
 м3 
2
 121
;
P1  P2  с 
l
1  2,4 d
2
d
3
Для воздуха при t=20C
=> l >10· d
3) Молекулярный режим (высокий вакуум):
UМ
3
3
d м 
 121

;
l  с 
Для воздуха при t=20C =>
l >10· d
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Связь проводимости вакуумопровода с представлением о
степени вакуума
Высокий
вакуум
Средний
вакуум
Низкий
вакуум
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Формулы расчёта проводимости трубопровода сложной формы
трубопровод длинной (L>lOd) может быть представлен цилиндрической трубой постоянного
диаметра его проводимость может быть определена (для воздуха М=29, Т=293 К, молекулярный
3
d
режим течения газа):
U  121
L
, м3с-1
(1)
где:d- расчетный диаметр трубопровода окончательной откачки, м;
расчетная длина трубопровода окончательной откачки, м.
Наличие клапана или затвора с таким же диаметром прохода учитывается увеличением длины
L на величина (2-6) dсобственно вакуумпровод окончательной откачки, который на расчетных
схемах может быть представлен трубопроводом постоянного сечения, если элементы конструкций
клапана, откачного гнезда, ловушки, золотника не диафрагмируют этот трубопровод. В общем
виде суммарная проводимость трубопровода состоящего из последовательно соединенных
1
U 
m
проводимостью Ui:
1
U
м3с-1
, (2)
i 1
i
где: m- количество последовательно соединенных участков трубопровода, учитываемых при
расчете;
Суммарная проводимость трубопровода состоящего из параллельно соединенных участков
(например, в золотниковых машинах):
U  U
,3с-1 (3)
где: n- количество параллельно соединенных участков. Проводимость тонкой диафрагмы
(длина L0 определяется только ее площадью А:
Т
n
i
i 1
U Д  3 6 ,4
М А
 Vi  A
где: Т- температура газа, К; М- молекулярный вес газа, Кмоль; Vi- объем газа, ударяющегося о
единицу поверхности в единицу времени, Vi=ll7 м3м-2с-1 (при М = 29, Т = 293 К).
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Расчёт проводимости трубопровода сложной формы
Проводимость трубопровода произвольной формы может быть найдена как:
U=UДK
(4),
где:UД- проводимость входного сечения рассматриваемого трубопровода.
Чтобы рассчитать проводимость участков следующих за "диафрагмой"), необходимо
ползоваться коэффициентом Клаузинга:
K   N  - N обр

N

где:N- суммарное число молекул, вошедших в трубопровод через впускное сечение;
Nобр- число "обратных" молекул, отраженных от стенок и вернувшихся через
впускное сечение.
При откачке в молекулярно-вязкостном режиме течения газа проводимость
трубопровода может определяться:
UМВ=0,9UМ+UВ
(5)
где:UМВ, UМ, UВ- расчетная проводимость трубопровода в молекулярно-вяэкостном, молекулярном, вязкостном режимах,
соответственно.
Таблицы формул для расчёта проводимости трубопроводов для газа, М=29, Т=293 К
Приведены на последующих слайдах
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Проводимость отверстия (диафрагмы) произвольной формы:
U
 м3
 117  S 
 с

;

Проводимость трубопровода в виде щели
U Щ  97  K 2
l/b
ab
2
l
3
T м 

;
M  с 
0,1
0,2
0,4
0,8
1
2
5
10
>10
0,036
0,068
0,13
0,22
0,26
0,4
0,67
0,94
3/8·ln
(l/b)
a>>b,
l- длина
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Формулы расчёта проводимости трубопроводов различного сечения
Форма сечения
трудопровода
Проводимость
длиной l, м
круг диаметром d
равносторонний
треугольник со
стороной а
Молукулярный режим
течения газа
3
U  121 d l
U 
прямоугольник со
сторонами а,b (ab)
U  1, 36  10  Pcp
3
121 d 1  d 2 
2
d
1
 d2 
l
a b
2
U  308 
среднее давление в трубопроводе, Па)
3
U  48 ,1 a l
кольцо
наруж. диам. d1
внутр. диам. d2
Вязкостный режим течения газа. (Рср-
U  299  Pcp
4
l
4
l


2
2

 d 1  d 2   Pcp
4
4
3
U  1, 36  10   d 1  d 2 

l
 d 1
 

ln
d




2
2
a  b . l
a
d
U  865
fPcp ab
l
3
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Формулы предыдущих слайдов показывают необходимость создания
универсального метода для расчёта проводимости сложных трубопроводов
- метода пробных испытаний (Метода Монте- Карло)
Моделирование поведения молекул подчиняется Закону Кнудсена:
dN отр
d
 N N  cos 
dN отр
d
 c  cos 
Молекулы при соударении с поверхностью
задерживаются на ней на время τ, поэтому угол
отражения не зависит от угла падения
В начале расчёта определяется проводимость впускной диафрагмы Uo
рассчитываеиого нами трубопровода
 л 
U 0  Ai  V i  11 , 7 A 
;
2 
 см с 
где А [см2] – площадь диафрагмы
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Формулы предыдущих слайдов показывают необходимость создания
универсального метода для расчёта проводимости сложных трубопроводов
- метода пробных испытаний (Метода Монте- Карло)
Моделирование поведения молекул подчиняется Закону Кнудсена:
dN отр
d
dN отр
 N N  cos 
d
 c  cos 
Молекулы при соударении с поверхностью
задерживаются на ней на время τ, поэтому угол
отражения не зависящий от угла падения является
случайной величиной, подчиняющейся указанному
«Кнудсеновскому» закону распределения и
определяемой генератором случайных чисел
Пример чисел, выдаваемых генератором случайных чисел, подчиняющихся
«Кнудсеновскому» закону распределения, и характеризующих угол
отражения от 0 до 180 0 с округлением значения до 100 (т.е. чисел от 1 до 17 с
функцией распределения по закону cos φ )
12
10
13
11
7
6
6
15
8
7
15
2
9
11
16
14
11
12
6
3
10
15
14
7
9
6
12
8
8
12
5
9
11
14
9
10
7
3
4
11
14
3
5
15
8
4
12
1
13
8
16
2
13
10
14
4
17
6
7
8
7
5
8
6
7
9
4
12
9
12
6
9
16
13
3
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Виды угловых распределений молекул по скоростям:
Слева- равномерное распределение; Справа- распределение Кнудсена
90
120
90
1
60
120
0.8
0.8
0.6
150
30
0.6
150
0.4
180
0
0
210
330
240
30
0.4
0.2
f ( )
0.2
0f (  )
0
180
0
0
210
330
300
270

1
60
240
300
270

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Методика расчёта проводимости сложного вакуумопровода
(«прямопролётного» клапана) методом Монте-Карло
Порядок расчёта проводимости
сложного трубопровода:
U U0 K  K 
N ПР
U  U 0  (1 
 1
N ОБР
N
N
N ОБР
).
N
Общее количество молекул равномерно распределенных по входной
площади:n
n
A
N 
A
i 1
i
 N  N i;
i 1
N i  N ПР
i
A
Где i- количество элементарных площадок на которые
разбиваем входную диафрагму.
;
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Вид в разрезе «углового» клапана КРУТ DY=40
Рассмотрим расчёт проводимости этого клапана
методом МК в плоской ( 2D ) системе координат
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Исходная геометрическая схема для расчёта проводимости этого клапана
КРУТ DY40
методом МК в 2D системе координат
( т.е. без учёта 3х мерного движения молекул)
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Схема, поясняющая необходимость изменения
размеров на плоском чертеже (2D схема) для учёта
3х координатного движения молекул.
Истинное соотношение размеров
S1
A1

S2
 D1  4
2

A2
D 2  4
2

D1
D2
2
2
Соотношение размеров при графическом расчете на
плоскости в 2D системе с учётом 3х координатного
движения молекул
S 2 РАСЧ  S 1
A2
A1
 D1
D2
D1
2
2
Пример распределения чисел от 1 до 17 с функцией распределения по закону cos φ
выдаваемых генератором случайных чисел (см. следующие слайды)
12
10
13
11
7
6
6
15
8
7
15
2
9
11
16
14
11
12
6
3
10
15
14
7
9
6
12
8
8
12
5
9
11
14
9
10
7
3
4
11
14
3
5
15
8
4
12
1
13
8
16
2
13
10
14
4
17
6
7
8
7
5
8
6
7
9
4
12
9
12
6
9
16
13
3
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Геометрическая схема клапана КРУТ DY40 для расчёта его проводимости
методом МК в 2D системе координат (На базе исходной с учётом 3х
мерного движения молекул)
1 запуск 10 молекул
«Впускная
диафрагма»
клапана
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Геометрическая схема клапана КРУТ DY40 для расчёта его проводимости
методом МК в 2D системе координат (На базе исходной с учётом 3х мерного
движения молекул) 2 запуск 10 молекул
«Впускная
диафрагма»
клапана
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Геометрическая схема клапана КРУТ DY40 для расчёта его проводимости
методом МК в 2D системе координат (На базе исходной с учётом 3х
мерного движения молекул) 3 запуск 30 молекул
«обратные»
моле кулы
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Геометрическая схема клапана КРУТ DY40 для расчёта его проводимости
методом МК в 2D системе координат (На базе исходной с учётом 3х
мерного движения молекул) 4 запуск 10 молекул
Видны
траектории
«прямых» и
«обратных»
молекул
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Далее рассмотрен расчёт проводимости нанозазора поляризационного клапана
методом Монте-Карло
(Схема работы клапана создаваемого на каф. МТ-11 представлена ниже)
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Пример расчёта проводимости нанозазора поляризационного клапана
методом Монте-Карло На рис. представлены элементы модели, использованные
для расчёта суммарной проводимости плоского уплотнения:
слева- сегмент кольца, покрытый сканами; справа-трубопровод, эквивалентный
сегменту кольца (сектор кольца уплотнения, с наложенными на него сканами)
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Элементы процесса расчёта проводимости нанозазора поляризационного клапана
методом Монте-Карло
В вакуумной технике для расчета проводимости сложных разветвленных
трубопроводов часто пользуются методом электрических аналогий. При
последовательном соединении n элементов вакуумной системы с известными
проводимостями Ui общая проводимость системы
(1)
n
1

U посл

i 1
1
Ui
При параллельном соединении n элементов вакуумной системы ее общая
проводимость U равна сумме проводимостей всех элементов .
n
U пар 
U
(2)
i
i 1
Для n последовательно соединенных элементов вакуумной системы формула
принимает вид:
(3)
n
1
P

1
P
  n  1
i 1
0
0i
где
P0i – вероятность прохождения молекулы через i-й элемент вакуумной системы
(коэффициент Клаузинга);
P0 – суммарная вероятность прохождения частицы через n последовательно
соединенных элементов вакуумной системы
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Вид сканов поверхностей элементов уплотнения клапана :
а) кремния;
б) алюминия;
в) меди
а)
1.969x1.969mkmx249nm
б)
1.969x1.969mkmx487.5nm
в) 5x5mkm x320nm
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Компьютерная реконструкция нанозазора в уплотнении
между контактирующими поверхностями по сканам поверхностей :
тарели (нижняя поверхность)
и седла (верхняя поверхность).
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Результаты расчёта проводимости нанозазора поляризационного
клапана методом Монте-Карло
Представлено изменение проводимости нанозазора U как функция изменения
напряжения u электрического поля при различных значениях нормальной силы
F: 1 – 476 Н; 2 – 493 Н; 3 – 526 Н; 4 – 594 Н; 5 – 610 Н; 6 – 627 Н; 7 – 644 Н
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ ИМ. Н.Э. БАУМАНА
Пример расчёта проводимости нанозазора поляризационного клапана методом
Монте-Карло выполненный М.В.Косинским
(разработка группы «Прецизионный вакуумный привод» каф МТ-11)