УДК 538.9 ББК 22.36 Г 12
Гаврилов А.И.
Кандидат физико-матаиатических наук, доцент кафедры физики Кубанского государственного технологического университета, тел. (861) 255-85-32, e-mail: [email protected] Москаленко Ф.В.
Старший преподаватель кафедры физики Кубанского государственного технологического университета, тел. (861) 255-85-32, e-mail: [email protected] Хлопова Т.П.
Кандидат педагогических наук, доцент кафедры физики Кубанского государственного технологического университета, тел. (861) 255-85-32, e-mail: [email protected]
Траектории частиц стеклообразующих материалов при получении заготовок для световодов парофазным осаждением
(Рецензирована)
Аннотация
В рамках разрабатываемой модели MCVD-технологии производства световодов получены траектории частиц стеклообразующих материалов. Анализируются особенности осаждения и изменения температуры на траекториях при смешанном конвективно-термофоретическом механизме.
Ключевые слова: MCVD-технология, моделирование, траектория частиц конденсата.
Gavrilov A.I.
Candidate of Physics and Mathematics, Associate Professor of Physics Department at Kuban State University of Technology, ph. (861) 255-85-32, e-mail: [email protected] Moskalenko F.V.
Senior Lecturer of Physics Department at Kuban State University of Technology, ph. (861) 255-85-32, email: [email protected]
Khlopova T.P.
Candidate of Pedagogy, Associate Professor of Physics Department at Kuban State University of Technology, ph. (861) 255-85-32, e-mail: [email protected]
Particle trajectories of glass-forming materials upon receipt of glass blanks for optical fiber by vapor-phase deposition
Abstract
The particle trajectories of glass-forming materials have been obtained within the framework of the developed MCVD technique of fiber manufacturing. Precipitation features and temperature changes on the trajectories are analyzed for the mixed convective-thermophoretic mechanism.
Key words: MCVD-technology, modeling, trajectories of condensate particles.
Успех в моделировании производства заготовок для световодов методом осаждения из паровой фазы (MCVD-технология [1]) во многом связан с реалистичностью предположений, положенных в основу анализа теплообмена. В континуальной конвективной модели осаждения обобщения приводят к следующей системе уравнений гидродинамики и теплообмена для парофазной смеси, продуваемой через цилиндрическую вращающуюся кварцевую трубку (рис. 1):
Рис. 1. Постановка задачи теплообмена
Э / ч ди. д / ч
— (гиг)+----------1----(гь7) = 0;
-\ \ г / -Ч V 2 / ’
дг д. дъ
(1)
и
ди V. ди V.
дг г д. г
-+и
диг
д2
1_ дП Рт дг
+ Рт£ • сое.-(Г - Тт)+
+ ,(Д„г - );
и
г дг г д.
ди.
дъ
г д. г 1 дП
+ У'
Рт г -д.
Ли + .и
. 2 Л 2
г д. г
вт - В - - Тт ) +
диъ и. диъ ди
и —-+——- + и
дг г д.
дъ
1 дП _ дТт ^
-----вт -8-Г-С08^^— + УДиъ
Рт дъ дъ
V
дТЁ и. дТЁ
иг—- + ——- + ь7 г дг г д. 2
дТВ
дъ
Рв СР£ г дг
.дії
дг
(2)
(3)
(4)
(5)
Уравнения (2) - (4) учитывают естественную конвекцию, записаны в приближении Буссинеска при обычных упрощениях для длинных труб [1]. Обозначено:
п= Р + Рт - £- г - С0. квазидавление (рт - плотность смеси при среднемассовой температуре);
д = у2= і! +1 -А+ 1 9 2
дг2 г дг г2 д.2
(6)
(7)
Остальные обозначения стандартные для подобных задач [1] (см. также рисунок 1); величины с индексом «£» и «к» относятся к газу и кварцу, соответственно. Граничные условия:
г = го; Т£ = Тк;
Л.А = 4 дТ‘
дг
дг
и = 0; и = 0;
(8)
(9)
г
Подобная задача решалась в упрощенной постановке [1] для постоянной плотности радиального теплового потока на поверхности трубки и в пренебрежении теплоемкостью и теплопроводностью опорной трубки. Однако в реальных технологических условиях эти ограничения не соблюдаются. Будем искать решение системы (1) - (5) методом разложения по малому параметру [2]
*=4°Т, (11)
Яе2
где число Грасгофа (вт - объемная сжимаемость, у£ - кинематическая вязкость)
Ог = г05-вт • £ ^) , (12)
4-Т, -V
число Рейнольдса (и0 - осевая скорость вдали от горелки)
Яе = ^1^. (13)
Плотность радиального теплового потока на поверхности ) находилась из приближенного решения задачи теплообмена без учета естественной конвекции путем сопоставления решения с экспериментально установленным распределением температуры на внешней поверхности трубки. В дальнейшем полученный результат используется в качестве нулевого приближения.
Было получено следующее представление [2]:
4(2)= ^І^), (14)
г0
где ¥(г) - установленная функция, изображенная (для одного из возможных наборов теплофизических параметров) на рисунке 2.
Рис. 2. Определяющие функции нулевого приближения (Яе=4,8)
В качестве нулевого приближения по е принималось: 1) для скорости
и
(0)
0; и
(0)
а- г;
и
(0)
и
Ґ \ 2
г
1 -
V -0 ,
(15)
2) для температуры
Т (0)
Т"" '
(0)
Г (г )(і - 3-2 + 2г3)+ г2 {2- - 3)[(г )іп(і + ё)-ф(г )]-(1 - г )Л¥(г)}. (16)
Здесь Г (г), (р{т.) также известные функции нулевого порядка (см. рисунок 2). Приняты обозначения:
число Прандтля
параметр вращения трубки
- _ . г - — • л -
г ’ - ’ ’
Рг _.^ Рр
N г
а- -
V
(21)
(22)
Поля скорости и температуры с учетом первого приближения представлены на рисунках 3 - 7.
Полученные результаты свидетельствуют о существенном влиянии естественной конвекции как на поле скорости, так и на поле температуры. Учет радиальной и азимутальной составляющих скорости потока приводит к принципиально иной картине осаждения частиц стеклообразующих материалов, взвешенных в потоке газа. Очевидно, что
термофоретический дрейф частиц к поверхности трубки со скоростями ьй ~ 1 — может
с
быть существенным фактором только при рассмотрении непосредственно осаждения. Для частиц удаленных от поверхности термофорез проявляется лишь как тенденция.
Данный вывод подтверждается видом траекторий выделенных малых объемов газа при их движении вдоль трубки, а также характером изменения температуры этих объемов.
Т 7
Л
1 / I
/ 2' г 4
// \
// 3 / / / \ \ \
V
Рис. 3. Зависимость безразмерной температуры от осевой координаты для некоторых г :
1) 7=1; 2) г =0,5; 3) 7=0
Рис. 4. Зависимость безразмерной температуры от радиуса для некоторых значений 2
1) 2 =50; 2) 2 =62; 3) 2 =67; 4) 2 =90
_ 40
1 /\ \ .
- ^ 2
4'
120
1
Рис. 5. Зависимость осевой составляющей безразмерной скорости от безразмерной координаты 2 для некоторых значений 7 и ф : 1) 7=0,5, ф =л; 2) 7=0, ф =0; 3) 7=1, ф =0; 4) 7=0,5, ф =0
Рис. 6. Зависимость радиальной составляющей безразмерной скорости от безразмерного радиуса для некоторых значений ф при 2 =35:
1) ф =ж; 2) ф =3ж/2; 3) ф =ж/2; 4) ф =0
Рис. 7. Зависимость угловой составляющей безразмерной скорости от безразмерного радиуса для некоторых значений ф при г =40: 1) ф =ж/2; 2) ф =ж; 3) ф =0; 4) ф =Зж/2
90
9
Рис. 8. Вид траекторий малых объемов в плоскости " г - ф", заданных различными начальными параметрами г и ф : 1 - траектория частицы, запущенной из точки
г = 0,5 и ф = 270°; 2 - траектория частицы, запущенной из центра трубки во входном сечении
На рисунке 8 можно видеть два несимметричных конвективных вихря в плоскости " г -ф", которые захватывают траектории в зависимости от положения объема во входном сечении трубки. Области действия вихря изменяются при увеличении угловой скорости вращения трубки и можно добиться исчезновения одного из них (в рассмотренных условиях при (й> 10 рад/с). Кроме того, траектории имеют сложную петлеобразную форму и при развертке вдоль оси г трубки (рис. 9), т.е. на некоторых участках наблюдается спутное течение.
Изменения температуры в зоне осаждения носят характер апериодических колебаний, существенно различных для различных траекторий (рис. 10). Указанные особенности неизбежно отразятся на группировании частиц при действии термофоретиче-ских сил. Следовательно, реальное представление о характере осаждения стекломате-риалов, а, значит, и об изменениях в показателе преломления заготовки, можно получить только на основании расчетов по совокупности траекторий частиц.
Рис. 9. Траектории малых объемов
t
Рис. 10. Изменение температуры на траектории в области осаждения стеклообразующих материалов
Примечания:
1. Тепло и массоперенос: теория и практическое приложение: процессы газофазного осаждения при формировании заготовок для оптического волокна: сб. науч. тр. / на-уч.-исслед. НИИ высшего образования; отв. ред. В.Л. Колпащиков, Ю.К. Кривоше-ев, ОТ. Мартыненко [и др.]. Минск: ИТМО АН БССР, 1983. С. 59-64.
2. Гаврилов А.И., Шапошникова Т.Л., Хай-рюзова Е.В. Эффективность осаждения частиц из парогазовой фазы как функция
// -
тика. Компьютер. Образование: тр. Восьмой междунар. конф., Москва-Пущино, 31 января - 5 февраля 2001. М., 2001.
References:
1. Warmth and mass transfer: theory and practical appendix: processes of gas-cycle sedimentation while forming the half-ready products for optical fibre: coll. of scient. works / Scientific research of the Scientific research institute of higher education; ed. by V.L. Kolpashchikov, Yu.K. Krivosheev, O.G. Martynenko [etc.]. Minsk: ITMO AN of the BSSR, 1983. P. 59-64.
2. Gavrilov A.I., Shaposhnikova T.L., Hayryuzova E.V. Efficiency of sedimentation of particles from steam-to-gas phase as a function of thermodynamic parameters // Mathematics. Computer. Education: proceedings of the Eighth international conf., Moscow - Pushchino, January, 31st - February, 5th 2001. M., 2001.