книги из ГПНТБ / Повх, И. Л. Магнитная гидродинамика в металлургии
.pdfрешаются при следующих граничных условиях:
иг = щ = и, — 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 = |
|
|
R° |
(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
