Математическое моделирование движения расплавов теллуритного стекла при получении двухслойных волокон методом двойного тигля Текст научной статьи по специальности «Физика»

Механика

Вестник Нижегородского университета им. Н.И. Лобачевского, 2007, № 1, с. 137- 143

УДК 621.396.22.02:535.8:532.536

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ДВИЖЕНИЯ РАСПЛАВОВ ТЕЛЛУРИТНОГО СТЕКЛА ПРИ ПОЛУЧЕНИИ ДВУХСЛОЙНЫХ ВОЛОКОН МЕТОДОМ ДВОЙНОГО ТИГЛЯ

© 2007 г. М. Ф. Чурбанов 1, Г.Е. Снопатин 1, Р.М. Шапошников 1,

В.В. Шабаров 2

1 Институт химии высокочистых веществ РАН, Нижний Новгород 2 Нижегородский госуниверситет им. Н.И. Лобачевского

vestnik [email protected]

Поступили в редикцию 25.12.2006

Обсуждаются результаты экспериментальных исследований течения расплавов теллуритных стекол. Установлено, что течение этих расплавов описывается моделью вязкопластической жидкости. С позиций гидродинамики рассмотрен начальный этап вытяжки двухслойного световода из двойного тигля. В расчетных исследованиях используется программный комплекс ANSYS CFX.

Постановка задачи

В августе-сентябре 2005 г. в ИХВВ РАН проведены две серии экспериментов по вытягиванию волокон из расплавов теллуритных стекол системы TeO2-WOз. Целью первой серии экспериментов являлось

определение теплофизических и транспортных свойств этих расплавов. Схема

экспериментальной установки и технология экспериментов полностью аналогичны

описанным в [1] экспериментам с расплавами селенида мышьяка.

Целью второй серии экспериментов являлась методическая отработка получения

двухслойных волокон из расплавов TeO2-WO3 методом двойного тигля. Схема

экспериментальной установки показана на рис.

1.

Двойной тигель, состоящий из контейнеров

сердцевины 2 и оболочки 3, выходного цилиндрического канала сердцевины 4 диаметром Асерд и кольцевого канала оболочки 5 диаметром Аобол, установлен в температурной камере 1. Под действием избыточных давлений, подаваемых на входы сердцевины Арсер и оболочки Аробол, расплавы стекол текут в конический канал. После выхода из этого канала формируется свободная двухслойная струя, которая вытягивается из отверстия диаметром Авых, постепенно остывает и переходит в двухслойное волокно, затем направляется на барабан.

Давление в контейнерах сердцевины и оболочки менялось таким образом, чтобы обеспечить доминирующую роль материалов сердцевины или оболочки в выходном двухслойном волокне. Экстремальные значения давления на входах в контейнеры сердцевины и оболочки составляли Артт= 30 кПа и Артах= 200 кПа.

Задача заключается:

- в построении математической модели движения расплава стекол системы TeO2-WO3 с использованием полученных экспериментальных результатов;

- в математическом моделировании течения расплавов при получении двухслойных волокон методом двойного тигля с целью исследования поведения поверхностей раздела «расплав-расплав», «расплав-воздух».

Реологическая модель расплава теллуритных стекол

м3/с

м3/с

Т=390°С

Т=410°С

м3/(

Рвх, кПа

м3/с

3.00Е-08-

Рвх, кПа

Т=400°С

Т=430°С

Рвх, кПа

40 60 рвх, кПа

Рис. 2

4.00Е-09

4.00Е-09

3.00Е-09

3.00Е-09

2.00Е-09

2.00Е-09

1.00Е-09

1.00Е-09

□.00Е+00

0.00Е+00

0

0

6.00Е-09

2.00Е-08

4.00Е-09

1.00Е-08

2.00Е-09

0.00Е+00

0.00Е+0С

0

40

80

20

Результаты опытов первой серии (рис. 2) свидетельствуют о неньютоновском характере течения расплавов стекол TeO2-WO3. При ньютоновском течении зависимость объемного расхода среды Q от давления на входе в фильеру рвх должна носить линейный характер; расход при этом равен нулю только при отсутствии перепада давлений на входе и выходе фильеры. В опытах течение расплава TeO2-WO3 не обнаруживается даже при существенных избыточных давлениях на входе в фильеру (Дрвх » 100 кПа). Такое поведение сплошной текучей среды типично для «пластических» жидкостей, в которых наряду с вязкостью присутствуют также и пластические свойства [2].

В теоретической гидродинамике изучение течений «пластических» жидкостей проводят обычно с использованием реологической модели Бингама-Шведова. В численном моделировании течение таких жидкостей исследуется либо с использованием модели «бивязкости», либо с использованием более общей модели псевдопластической жидкости [3], согласно которой связь между напряжениями т и скоростями деформации в случае течения жидкости в длинной круговой цилиндрической трубе имеет вид:

ёг

^ ,-7ТЛ \п 1

ёУу ёг

(1)

где А и п < 1 почти постоянны в широких интервалах напряжений и скоростей деформаций, но зависят от температуры, т' ~ кажущийся коэффициент вязкости, убывающий с ростом скорости деформации, г - полярная координата.

Фактически модель псевдопластической жидкости является обобщением как линейных моделей вязкопластической жидкости, так и ньютоновской модели жидкости (п = 1 в (1)). Псевдопластическая жидкость течет при любых т, однако с уменьшением касательного напряжения т свойство текучести проявляется в меньшей степени (кажущаяся вязкость т’ ® ¥) и при малых т может экспериментально не обнаруживаться. Поэтому для теоретического анализа течения расплавов TeO2-WO3 принята модель псевдопластической жидкости. Константы модели А и п должны определяться по результатам экспериментов для каждого температурного режима.

Определение реологических констант в гидродинамической модели расплава теллуритных стекол

Авторами в [1] показано, что при медленном течении вязкопластической жидкости контейнер фильеры практически без искажений передает давление на входе в фильеру в ее канал. Это обстоятельство используется для определения реологических констант нелинейной модели вязкопластической жидкости по результатам экспериментов.

С позиций нелинейной модели вязкопластической жидкости

(псевдопластической жидкости) получено выражение для распределения осевых скоростей Уу при течении жидкости в круговом цилиндрическом канале радиуса Я вследствие осевого градиента давлений dpldy:

Ґ л , \11 п

Уу ) =

п

П + 1

}( п+1)1 п

- Г

1 dp 2А dy

(п+1) I п )

X

X ( Я(

Объемный расход при псевдопластической жидкости цилиндрическом канале равен

0 = \0>Уу 2р г^ =

3п +1

\1/ п

Я (3п+1)| п (3)

Рис. 2

жидкости в среде БоЛтаи Ро^га^айоп была создана программа, позволяющая определять А и и по гипотезе (1) исходя из минимума рассогласования расчетных результатов (3) с результатом и экспериментов. При этом полагалось, что на выходе из фильеры давление равно атмосферному.

В таблице представлены полученные результаты по величинам входных параметров модели А и п для различных температур.

Таблица

Т, °С п А, Пасп

390 0,17 5650

400 0,22 5600

410 0,73 3600

430 0,48 1850

(2)

течении

круговом

1 ёр

2 А ёу

ч ^

;тном случае п = 1 из (2) следует й закон Пуазейля. При 0 < п < 1 расход более чувствителен к изменению радиуса трубы и градиента давления в сравнении с течением Пуазейля.

Для идентификации реологических параметров А и п модели псевдопластической

Результаты расчетов по объемному расходу Q расплава TeO2-WO3 в зависимости от давления р на входе в цилиндрический канал фильеры для Т = 390°С показаны на рис. 3. Расчет объемного расхода по модели псевдопластической жидкости вполне

удовлетворительно согласуется с полученными экспериментальными результатами.

Зависимость распределения осевых

скоростей по радиусу канала отличается от распределения Пуазейля и состоит из двух участков: отрезка почти горизонтальной прямой и параболы; на некотором расстоянии от оси канала эти участки склеиваются. Аналогичное распределение скоростей имеет место при расчете, базирующемся на модели Бингама -Шведова.

Математическое моделирование течения расплава стекол системы TeO2-WOз в нижней части двойного тигля и при выходе

О, и /с

/

<

X

Р, Па

эксперимент, Т = 390 С, расчет по (3) Т = 390°С

из конического канала

Течение вязкопластической жидкости изучалось с использованием программного пакета АКБУБ СБХ, ориентированного на решение гидрогазодинамических задач. Идеология пакета базируется на методе конечных объемов [4]; пакет позволяет

моделировать течения нелинейных жидкостей, в том числе с границами раздела жидких сред (свободными границами).

Рассматривается течение двух расплавов системы TeO2-WO3 в нижней части двойного тигля с их последующим выходом из конического канала и начальным этапом

вытягивания. Расчетная область (рис. 4)

имитирует нижнюю часть цилиндрического канала сердцевины ОБСБ, нижнюю часть кольцевого канала оболочки ДЕЕК, конический канал двойного тигля ООБКЬ и, для моделирования выхода струй, часть

пространства в объеме кругового цилиндра АОМЫ высоты АО = 0,01 м радиуса АЫ = 0,006 м.

Г еометрические размеры элементов

расчетной модели тигля принимались равными соответствующим геометрическим размерам двойного тигля, используемого при проведении второй серии опытов: БС = 0,003 м, ЕЕ = 0,005 м, БО = 0,04 м, высота конического канала ОО = 0.01 м, радиус выходного сечения

конического канала ОЬ = 0.004 м, толщина разделительного стекла между каналами сердцевины и оболочки СЕ = 0,001 м. Радиус

скругления в месте сопряжения каналов сердцевины и оболочки принят равным 0,0005 м.

Уравнения осесимметричного течения

расплавов и воздуха имеют вид:

р і (дУ^ + У п дУ^ + дУ^) =

' дґ дг ду

=_э^+^ (2т ) +

дг дг дг

(4)

д дУуі дУ ■ и- дУуі у.

+ду (2т (-^+д^))+—(-дгу- _ —),

ду дг дг г дг г

дУу, дУу

дґ

дУ,

р. (* + У ___у^ + У ) =

( - гі дг уі ду ^

др, д , дуу, дУгі

= _ р,8 + ^Т (иі Ьг1 + ^))■

ду

дг дг де

д дУу, и, ґдУуі дУгі

+ ^Г Ря-у-) + (^ _(5)

ду ду г дг ду

дУуі+Уі+дУі=0.

дг г ду

(6)

В уравнениях количества движения (4), (5) и неразрывности (6) индекс г = 1 соответствует расплаву сердцевины, г = 2 - расплаву

оболочки, г = 3 — воздуху; V„■, -уг — компоненты соответственно осевой и радиальной скорости, рг — плотность, рг — давление, g = 9,81 м/с — ускорение свободного падения, т — вязкость.

Вязкость расплавов TeO2-WO3 является нелинейной функцией скоростей деформации следующего вида

Л = А,Н"—\ (7)

где Нг — инвариант тензора скоростей деформации,

1 д- V ■ 2 V ■ дУу. 2 Нг = (т[(-^ — ~)2 + (~ ^^Т^)2 + (8)

6 дг г

г ду

ґдУуі дУп. . 1 ґдУу1 дУп. ^Ц2

+ (— ----------—)2] + — (—— + ——)2)112.

ду дг 4 дг ду

Характеристики расплавов соответствуют температуре Т» 400°С и приняты следующими:

- для сердцевины

п1 = 0,22, А: = 5600 Пас022, р = 6000 кг/м3;

- для оболочки

п2 = 0,25, А2 = 5600 Па с025, р2 = 6000 кг/м3. Вязкость и плотность воздуха заданы и постоянны.

Решение задачи разделяется на два этапа: течение расплавов до выходной границы АЫ расчетной области и их течение при моделировании начальной процедуры вытяжки.

В начальный момент времени расплавы занимают верхнюю часть цилиндрического и кольцевого каналов, в их нижней части, а также в коническом канале и в области АОМЫ находится воздух. Скорости во всех точках расчетной области равны 0, действующие давления — гидростатические. Далее, на входы БС и ЕЕ подаются некоторые избыточные давления, начинается движение расплавов с их последующим взаимодействием в коническом канале между собой и воздухом.

Граничными условиями задачи (АБ — ось симметрии) являются условия прилипания на твердых стенках СБЕ и ЕКЬМЫ: —у = —г = 0, а также заданные избыточные давления на входах в каналы сердцевины БС и оболочки ЕЕ; на выходной границе АЫ также задается значение избыточного давления.

На не известных до решения задачи поверхностях раздела сред ставятся кинематическое граничное условие (условие непрерывности скоростей при переходе через границу раздела) и динамическое граничное условие, состоящее в равенстве напряжений, действующих на элементарную площадку границы раздела сред.

В численной реализации определение мгновенного положения границ раздела сред осуществляется методом объемного слежения (VOF), согласно которому интегрирование уравнений движения сред осуществляется с учетом перемещений частиц — пассивных маркеров [5, 6]. Пассивные маркеры движутся со скоростью жидкости и указывают вид жидкости, находящейся в рассматриваемый момент времени в каждом конечном объеме расчетной области.

Расчеты проводились на различных сетках. Оптимальным с точки зрения точности и времени счета является сетка с равномерным разбиением по г и сгущением по у к ОЬ, то есть к выходу из конического канала. Далее приведен ряд результатов решения рассматриваемой задачи на сетке, состоящей примерно из 55000 конечных объемов. До начала процедуры вытяжки — первого этапа расчета — граничные условия в терминах избыточного давления имели вид:

р = 200 кПа на БС,

р = 120 кПа на ЕЕ, (9)

р = 0 на AN.

После достижения струями расплавов сердцевины и оболочки границы AN начинается моделирование процедуры вытяжки. В рассматриваемом случае одновременно снижается давление на входе в канал сердцевины ВС до величины р = 50 кПа, повышается давление на входе в канал оболочки ЕЕ до значения р = 150 кПа. Таким образом, в данной задаче моделируется получение двухслойного волокна, в котором доминирует материал оболочки.

Граничные условия на AN модифицируются. На части выходной границы AN, прилегающей к А и имеющей в радиальном направлении протяженность, равную требуемому радиусу двухслойного волокна, задается некоторая отрицательная величина избыточного давления. На другой части AN задается значение осевой скорости (другие компоненты скорости — нулевые). В процессе вытяжки значения давления и осевой скорости постепенно увеличиваются по абсолютной величине до значений, при которых через AN будет проходить струя (волокно) вытягиваемого радиуса. Промежуточные результаты расчета контролируются, поскольку резкие изменения давления и осевой скорости на AN приводят к разрыву сплошности струи, что соответствует, очевидно, обрыву волокна и, математически, к потере устойчивости счета.

На рис. 5-8 показаны «компьютерные кадры» течения расплавов и воздуха до начала процедуры вытяжки. В начальном положении (ґ = 0) расплавы занимали 95% цилиндрического и кольцевого каналов; предварительные расчеты показали, что уменьшение этой величины не приводит к изменению последующего движения расплавов и воздуха в коническом канале.

Приведенные расчетные материалы показывают, что взаимодействие расплавов и воздуха носит достаточно сложный характер. На нижней границе межканального стекла-разделителя (зона Б на рис. 4) в начальные моменты времени образуется воздушная полость, которая смыкается с атмосферой через узкую воздушную горловину. Эта горловина существенно влияет на движение расплавов в коническом канале и ниже его; с развитием процесса воздух из полости уходит через горловину, а сама горловина переносится вниз по потоку за выходную границу расчетной области.

В «кадре», приведенном на рис. 8, горловина уже находится за границей расчетной области.

-ЇМ;

Рис. 5

Рис. 7

Рис. 6

Рис. 8

/'

/

/ / |:ы1 -ІЇ.І; і / / 1 1 У.Ш.

■ "ГЛТ.4: 1 і.:::іЧУл

— ір-и II II. 1 “рйЯПф--

1 ііьсі лі -оі; 1 і |миі. ІІ ^.1

9 ■ / Рис. 10

Характерно, что в ряде случае в физическом Отметим еще следующий полученный эксперименте в зоне О также наблюдается расчетный результат. Если нижняя часть стекла

воздух, который может уноситься с расплавами - разделителя каналов имеет острые кромки, то

в процессе вытяжки. по результатам математического

моделирования получить двухслойное волокно

не удается: между расплавами, в том числе и при их вытяжке, наблюдается воздушная прослойка.

На рис. 9, 10 представлены результаты математического моделирования по начальной процедуре вытягивания двухслойного волокна. Отметим, что моделирование собственно

получения двухслойного волокна в настоящих расчетных исследованиях не рассматривалось.

Вследствие резкого падения давления на входе в канал сердцевины и относительно больших скоростей расплава оболочки

поверхность раздела расплавов весьма быстро, примерно за 1,5 с, смещается в верхней части конического канала к оси симметрии (в сторону канала сердцевины). В нижней части смещение поверхности раздела происходит медленнее, что обусловлено, по-видимому,

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

сформировалась, струя же оболочки только начинает формироваться).

Выводы

1. С позиций модели вязкопластической жидкости исследовано течение расплавов стекол системы TeO2-WOз в различных каналах.

2. В исследованном диапазоне температур проведена идентификация построенной математической модели по результатам опытов. Движение рассматриваемой сплошной среды — расплавов теллуритных стекол — описывается моделью вязкопластической жидкости.

3. Для исследования течения вязкопластических жидкостей с поверхностями раздела сред использован программный комплекс СБХ. Программный комплекс удачно схватывает основные моменты исследуемого явления, и, по-видимому, может использоваться для исследования аналогичных и более сложных задач химической технологии. В связи с этим следует привести перечень задач математического моделирования, решение

которых будет способствовать более глубокому пониманию гидродинамических особенностей течения расплавов теллуритных и халькогенидных стекол, отработке

технологических процедур и, в конечном счете, изготовлению световодов высокого качества. Это:

— численное исследование влияния

месторасположения входа контейнера оболочки в оболочку с целью получения двумерного потока на выходе канала оболочки;

— дальнейшая математическая отработка процедуры формирования двухслойного волокна;

— численное исследование течения

расплавов теллуритных и халькогенидных стекол в неоднородных температурных полях, моделирование «загустевающих» струй;

— аналитическое и численное исследование устойчивости течения струй расплавов по отношению к внешним возмущениям, моделирование течений в каналах различной геометрии в условиях пульсаций температуры и давления;

— построение математической модели

течения расплава стекол с кристаллами наноразмеров и ее численная реализация.

Работа выполнена при финансовой поддержке РФФИ (грант 05—08—137а).

Список литературы

1. Математическое моделирование движения сплошной среды при вытягивании световода через фильеру / М.Ф. Чурбанов [и др.] // Вестник ННГУ. Серия Механика / Н. Новгород: Изд-во ННГУ. — 2004. - Вып. 1(6).— С. 104-115.

2. Лойцянский, Л.Г. Механика жидкости и газа/ Л.Г. Лойцянский. — М.: Наука, 1971.

3. Литвинов, В.Г. Движение нелинейно-вязкой жидкости / В.Г. Литвинов. — М.: Наука,1982.

4. Флетчер, К. Вычислительные методы в динамике жидкостей. Т. 1, 2 / К. Флетчер. — М.: Мир,1991.

5. The MAC method / J.E. Welch [et. al]. - Los Alamos Scientific Lab. Rept. № LA-3425. — Los Alamos: 1966.

6. Harlow, F.H. Numerical study of large amplitude free surface motion / F.H. Harlow, J.E. Welch // Phys. Fluids. - 1966. - 9. P. 842—85.

MATHEMATICAL MODELING OF MOLTEN TELLURITE GLASS FLOWS DURING FABRICATION OF TWO-LAYER FIBERS BY THE DOUBLE-CRUCIBLE TECHNIQUE

M.F. Churbanov, G.E. Snopatin, R.M. Shaposhnikov, V. V. Shabarov

We discuss the results of experimental studies of molten tellurite glass flows. It is found that the molten flows are described by the model of a viscoplastic fluid. The initial stage of drawing a two-layer optical fiber from a double crucible is analyzed within the framework of the hydrodynamical approach. The ANSYS CFX