Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Повх, И. Л. Магнитная гидродинамика в металлургии

.pdf
Скачиваний:
28
Добавлен:
21.10.2023
Размер:
10.04 Mб
Скачать

решаются при следующих граничных условиях:

иг = щ = и, — 0

при

г = R0,

 

иг , ыф, и, Ф оо

при

г = 0 .

 

Ищем периодические по z решения в виде:

 

иг =

ept и (г)

cos k z ;

 

ыф — ept v (г)

cos k z ;

(187)

uz =

ept w (r)

sin k z ;

 

q =

ept q (r)

cos k z ,

 

где p — некоторая комплексная постоянная,

 

k — волновое число.

в систему уравнений (186),

Подставив решения (187)

получим:

 

 

 

 

v [DD^ — k2-----v = -^3.

 

 

 

v \ D D * - & -

Bt a

v (D# V0) и = 0 ;

(188)

 

v(DD%— k2 -------) w = — k q ,

аиз уравнения неразрывности Z)* ы = — k w ,

где

D =

 

r

;

didD r *—(— =1 / r .

 

Исключив до из уравнений

(188), воздействуя на по­

лученное уравнение оператором Я,

введя a = p R I /v,

k2 — a2/Rо

и используя

 

асимптотическое

'выражение

скорости Уо Для Яыд-<1,

получим следующую задачу:

(D £>,- а2 -

а) (Я Я * - а2) и =

а>Hal a2

 

 

----- ^ (1 -г 2) о,

 

 

 

 

 

 

4 м

 

(189)

 

 

 

 

 

 

 

( Я Я * - а 2 — a — H a2) v =

 

 

(1 — 2 г2) ы.

 

 

Д

 

 

4 v

 

 

В уравнениях

(189) 'все члены имеют

размерность

скорости; ы=и =

0 при г = 0

и г A, D u = 0 при г = 0 , г =

= \.

 

 

 

 

 

 

 

Исключив и из уравнений системы (189) и положив

сг=0, согласно принципу смены устойчивости, придем к одному уравнению шестого порядка для определения v\

100

(D D* — a2)2 (D D* — a2 — Я a2) « =

T„ a2 (1—r2) Q* o, (190)

где

 

 

0,707

1

 

Q *= J Dt V0dr +

J

Dt V0dr

о

0,707

и Q* — ступенчатая аппроксимация завихренности. Видно, что в рассматриваемой задаче собственными

числами являются числа Гартмана, так как число Тэй­ лора в свою очередь есть функция числа Гартмана.

Задала формулируется следующим образом: найти наименьшее положительное число Гартмана при неко­ тором фиксированном Rea, , соответствующее некоторому значению а, определяющему критическую длину стоя­

чей волны (вторичного течения).

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

разложим G = ( D D *—b2)v

по полной

системе

ортого­

нальных функций:

 

 

 

 

G = '2l PJS1 {aJ r),

(62 = a2 +

tfa 2),

(191)

Si (a,- r) = A] Jx (ay r) + Cy lx (a; r) ,

 

aj — корни уравнения:

 

 

 

 

Л (ay)

h

(ay)

 

(192)

J0 (ay)

J0

(ay.)

 

 

 

Обозначив:

«/ ('') = A j J i (ay r),

(r) = Cj ;i К r) >

и учитывая, что

D D* Uj = — a2 Uj, D D* v; = a2 ,

получим следующее решение для a:

101

Величину pj определим из граничных условий:

2 а?

и, (1)

(194)

р,- = —j— -—

---------.

'а ) - h {оа)*

Подставим значение v в уравнение (190), умножив обе части уравнения на rSi(akr) и проинтегрировав по г

от 0 до 1, .получим уравнение для определения критиче­ ского числа Гартмана:

Rl* (<** +

а2)

(а®'+ а2+ Я а*) +

QJk (а2— а2) X

Х'(«/ — а2 — Н а2)

Тма*

P / h (b) -\-

а2 Mjk aj N,k

 

 

= 0,

 

 

 

 

 

 

а)— а4

 

 

 

 

 

 

 

 

(195)

где

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rjk

=

| ги;- (uk + vk)

d r ;

 

 

 

 

 

о

 

 

 

 

 

 

Qik =

l

 

(Ub + vk)

 

 

 

 

j*rv,

dr;

 

 

 

 

 

6

 

 

 

 

 

h (6) =

-J" f

~

 

r 7i (br) (“* +

y*) d r >

 

о

 

 

 

 

 

 

Mjk =

f

(1 — r2) r

(iиу+ vj)

(и* +

иА)

dr;

 

0

 

 

 

 

 

 

 

yV/ft = у1

jЛ (1 — г2)

r

(Uj vj)

(«ft - f

«ft)

dr .

Критические значения числа Гартмана и волнового числа рассчитаны по ступенчатой аппроксимации завих­ ренности

0,707

1

J

 

 

2* = J

D t V 0d r +

D*V0d r .

(196)

0

0,707

 

 

Зависимость критического числа

Гартмана

от ReU]

показана на рас. 43,

102

4.ЛАМИНАРНО-ВИХРЕВОЕ И ТУРБУЛЕНТНОЕ ТЕЧЕНИЯ

Трехмерное

ламинарное

 

 

 

 

течение, возникающее вслед­

 

 

 

 

ствие

неустойчивости

пер­

 

 

 

 

вичного

чисто

азимутально-

 

 

 

 

ного течения и представляю­

 

 

 

 

щее собой систему меридио­

 

 

 

 

нальных вихрей, накладыва­

 

 

 

 

ющуюся на это

азимуталь­

 

 

 

 

ное течение,

получило

наз­

 

 

 

 

вание ламинарно-вихревого.

 

 

 

 

Теоретический анализ ла­

 

 

 

 

минарно-вихревых МГД-те-

Рнс.

43.

Зависимость

критического

чений

весьма

сложен,

по­

значения

числа Гартмана от вели­

этому мы ограничимся рас­

 

 

чины Re

 

смотрением

течений такого

 

сосуде, помещенном

типа в гомополярнике и длинном

во вращающееся однородное магнитное поле.

полная

В первой из задач [69] численно

решается

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

На рис. 44 показаны профили азимутальной скорос­ ти, построенные по различным решениям. Видно, что вторичные течения, существенно уменьшая азимуталь­ ную скорость, значительно изменяют таклее ее .профиль, делая его более заполненным. Скорости вторичных те­ чений соизмеримы со скоростью основного течения.

Как было показано в предыдущем параграфе, ази­ мутальное течение во .вращающемся магнитном поле ус­

тойчиво лишь при малых значениях

числа

Гартмана,

соответствующих числам Рейнольдса

Re<. 1.

При 'более

высоких скоростях вращения жидкости возникают вто­ ричные вихревые периодические по z течения. Распреде­

ление радиальной компоненты скорости, дающие представление о структуре этих вихрей, может 'быть вы­ числено из второго уравнения системы (189), как толь­ ко будут определены критическое значение числа Гарт­ мана, критическое волновое число и коэффициенты раз­ ложения (19.1).

Возмущение азимутальной скорости, появляющееся при образовании ламинарно-вихревой структуры тече­ ния, можно рассчитать по выражению (193).

ЮЗ

На рис. 45 показаны профили азимутального и ради­ ального возмущений.

Видно, что при возникновении вихрей Тэйлора про­ филь азимутальной скорости становится более плоским и возрастает градиент скорости в пограничном слое, что

0,8

 

 

 

 

 

 

 

 

 

 

 

/

 

 

 

 

 

 

 

 

 

 

 

0,4

 

 

 

 

 

 

 

 

 

 

 

I

01

о,4

\ 0,8 f

 

 

 

 

 

 

 

 

 

3

\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-0,4

 

 

Рнс.

44.

Профили

азимутальной

 

 

Г

2

скорости

течения

жидкости

в

го-

 

 

 

 

мополярнике при 2 = 0:

 

в

 

\

 

 

/ — система,

не

ограниченная

 

 

 

осевом направлении; 2 — течение

в

 

- 0,8

 

 

ограниченной

системе

без

учета

 

 

 

вторичных

перетеканий; 3 — расчет

 

- 1,0

 

 

с

учетом

вторичных

перетеканий

 

 

 

т/ __

/ \ & r i t

 

 

 

 

 

 

 

 

 

 

у о------------ , где }\-—плотность тока

 

 

 

 

на

 

vp

 

электроде радиуса

Рис. 45. Профили азимутального 2

внутреннем

 

 

г,,

/ — высота системы

 

 

и радиального / возмущений ско­

 

 

 

 

 

 

 

 

 

 

 

рости

 

объясняет наблюдаемое в эксперименте резкое увеличе­ ние напряжения трения на стенке.

Профиль радиальной компоненты скорости свиде­ тельствует о возникновении двойного вихревого слоя, причем общая граница этих слоев близка к поверхнос­ ти, на которой обращается в нуль завихренность пер­ вичного азимутального течения, т. е. выполняется усло­ вие неустойчивости Рэлея.

Итак, центробежная неустойчивость течения во вра­ щающемся магнитном поле проявляется в виде системы вихрей, возникающих в неустойчивом слое 0,707< г^ 1. Эти вихри возбуждают в устойчивом ядре вторую систе­ му вихрей, вращающихся в противоположном направ-' лении.

104

Ламинарно-вихревое течение у криволинейной стен­ ки при определенных значениях числа Рейнольдса ста­ новится неустойчивым, что проявляется в дроблении ме­ ридиональных вихрей и нарушении их пространственной ориентации. Упорядоченная однозначно детерминиро­ ванная вихревая структура, присущая ламинарно-вих­ ревому течению, разрушается, возникает новая вихре­ вая структура со случайным пространственно-времен­ ным распределением скоростей и завихренности — тур­ булентное течение.

Точной количественной теории МГД-турбулентности, и в частности МГД-турбулентности вращающейся жид­ кости, не существует. Немногочисленные работы в этой области выполнены на основе теории размерностей и от­ носятся к полуэмпирическим методам исследования.

Теория турбулентного течения жидкости в гомопо­ лярнике, разработанная Е. П. Ваулиным и В. И. Бабец-

ким [ДО], была экспериментально подтверждена в рабо­ те [61].

Размерный анализ турбулентного течения прово­ дящей жидкости в кольцевом канале с движущейся гра­ ницей при круговом вращающемся магнитном поле, был сделан А. В. Ивановым-Смоленским . и Г. С. Тамояном

[62].

В работе установлено существование трех характер­ ных режимов течения:

1) магнитное поле отсутствует (На 0): f(6/Rl t . Re , т > V2) = 0 ;

2) интенсивность магнитного поля мала ( На ~ 1) —

режим соответствует ламинарно-вихревой структуре:

f № i ,

т% Vй, Re, Н a?/Rе) = 0;

 

3) интенсивность магнитного поля велика

(На^>\):

f ( т°/рУ2, Ha/Re) = О,

 

 

где S — ширина зазора между цилиндрами.

 

экспе­

Достоверность

последней совокупности была

риментально проверена в диапазоне чисел

Рейнольдса

О< .Re^0,42-105 и Гартмана 0 < # а < 7 5 .

при

На^> 1

Зависимость %°!pV2=f(Ha/Re) (рис. 46)

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

105

кривой, что подтверждает правильность выбора крите­ риев для третьего режима течения.

- Явная зависимость коэффициента сопротивления от аргумента HalRe хорошо согласуется с выводами Гар­

риса [22], отметившего сходство между течениями и бе-

Д

 

 

 

 

 

 

 

АД

 

 

 

 

 

o

A

 

 

 

 

 

 

А Д

 

 

 

 

 

<

 

 

 

 

 

 

о

Д

 

 

 

 

 

 

О Ад

 

 

 

 

 

 

.

 

 

 

 

 

 

 

а е * д

По Норгетройду

 

 

 

Ал

 

 

 

 

*

 

Re-0,3-tO5

 

 

 

* 9 # SД

 

e Re-0,175-105

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, 1

° ^ - ^ ~ ~ R e =

Fv/v

*

° ° V

S

'

Re*0,m05 ° №’0,105-10*

 

 

 

 

х Re-0,27;/0^ A Re-0,27-105

 

 

 

 

Re=0,78-105 a Re-0,I35‘I0S

 

 

 

 

®Re--0t8-Ws vRe--Q2M05

0

 

W

20

 

30

Ha/Re■IOl,

 

Рнс. 46. Зависимость коэффициента сопротивления от Ha/Re:

I — теоретическая

кривая

для ламинарного

течения

между

цилиндрами

гущем и постоянном 'магнитных полях при малом сколь­ жении. Сравнение опытных точек с теоретической кри­ вой для ламинарного режима течения показывает, что увеличение интенсивности |бегущего 'магнитного поля не приводит к ламинаризации потока.

5. ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ МГД-ВРАЩЕНИЯ

Теоретические исследования движения проводящей жидкости во вращающемся магнитном поле сравнива­ лись с экспериментальными результатами. Эксперимен­ тальные измерения параметров движущейся жидкости проводились в широком диапазоне чисел Гартмана. Из­ меряли как некоторые интегральные характеристики по­ токов '[26, с. 166; 63], так и локальные [44, с. 167; 67].

106

Изучали такйё .Интегральные характеристики движе­ ния, как момент трения вращающейся жидкости, при­ ложенный к боковой .поверхности 'Цилиндрического со­ суда, и среднюю скорость стационарного и нестационар­ ного вращения свободной поверхности жидкого металла,

Рис. 47. Схема" экспериментальной установки

107

движущегося под действием вращающегося магнитного поля.

Для измерения момента трения вращающейся жид­ кости о стенку сосуда использовалась установка, .пока­ занная на рис. 47 [26, с. 166].

Стеклянный цилиндр 4 со ртутью подвешивался на крутильных весах, состоящих из упругой нити I, цанго­ вых зажимов 2, зеркальца 3, опорной колонны, круго­ вой шкалы радиусом R = 800 мм -и источника света.

Нижняя часть цилиндра 4 была вклеена в переход­ ный цилиндр 6, на который туго насаживался поплавок

7. Эта система устанавливалась по оси расточки статора

5 конструкции Таллинского политехнического института,

имеющего следующие параметры:

Число пар полюсов, р

.......................... 1

60

Номинальный фазный ток / ф, А . .

Диаметр расточки D,

, м м .....................68

Длина

расточки L, м м .............................. 200

Для гашения

крутильных колебаний

упругой нити

поплавок 7 погружался в ванну 8 с вазелиновым мас­

лом. Конструкция установки позволяла подвешивать ци­ линдры различных диаметров. Все детали, находящиеся в магнитном поле, были выполнены из стекла либо орг­ стекла.

Трехфазный статор включался в сеть трехфазного тока 380 В, 50 Гц, через трансформатор ТСШ-4 и регу­ лируемый автотрансформатор РНТ-220-12, которым обеспечивалась плавная регулировка напряжения на об­ мотках статора.

Ток и напряжение измерялись комплектом К-50,

по­

зволяющим производить измерение

режимов

в фазах

без выключения установки.

 

 

 

Напряженность магнитного поля в расточке статора

измерялась катушкой ( w = 2 50 витков, s = 7.3 мм2),

на­

пряжение с которой подавалось на

ламповый

милли­

вольтметр МВЛ-2М.

подсчитывалось

по

Мгновенное значение индукции

формуле

 

 

 

4,44 f w s

Распределение радиальной компоненты магнитного поля по радиусу поперечного сечения статора близко к линейному, а по длине расточки достаточно однородно.

108

Были проведены измерения величины момента тре­ ния жидкого металла о стенку цилиндрического сосуда в диапазоне чисел Гартмана от 0,064 до 32.

Сравнение экспериментальной кривой тгф = /( Я с д) с теоретической, полученной .при рассмотрении одномер­ ного осредненного по координате <р азимутального МГДтечения [85], показало, что для момента трения наблю­ дается удовлетворительное соответствие теоретических и экспериментальных результатов для малых чисел Гар­ тмана (Яад<С 1). При увеличении числа Гартмана до 6 экспериментальные кривые, оставаясь качественно близ­ кими теоретическим, дают заметное количественное рас­ хождение, причем экспериментально измеренный мо­

мент трения выше теоретического. При

более высоких

значениях числа Гартмана зависимость

измеренного в

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

Подобные эксперименты по определению момента трения жидкого металла о стенку сосуда при Яад> 16 и 104 описаны .в работе И. Пуарье. Отмечено, что экспериментально измеренный момент более чем в 3 ра­

за больше

вычисленного теоретически (при R — 10-2

м,

соо=б2,'5 рад/с,

Я =63-Ю 3 А/м измеренный момент

ра­

вен 6-10-5

Н-м;

теоретически рассчитанный момент

ра­

вен 1,9-10-5 Н-м). Опыты показали, что момент трения для больших чисел Рейнольдса и Гартмана не зависит от вязкости используемой жидкости.

Описанная в работе [26, с. 166] установка была усо­ вершенствована с целью измерения момента трения для чисел Гартмана от 0,005 до 0,04. Питание статора было стабилизированным, ток в каждой фазе выравнивался реостатами.

Индукция переменного 'магнитного поля измерялась тарированной катушкой. Постоянная кдгушки составля­ ет 1Г2±1 витков-см3. Напряжение с катушки пода­ валось на ламповый милливольтметр iB3-44.

Сравнение экспериментальной кривой т ,-ф = /(Я а д)

с теоретической, полученной при рассмотрении одномер­ ного осредненного по координате <р азимутального МТД-течения [85], подтвердило, что для момента тре­ ния между ними наблюдается удовлетворительное соот­ ветствие в области |0,005^Яад^'0,01. При увеличении числа Гартмана экспериментальная зависимость прямо­

109

Соседние файлы в папке книги из ГПНТБ