книги из ГПНТБ / Повх, И. Л. Магнитная гидродинамика в металлургии
.pdfВ работе [57 с. 315—324] решена задача о лами нарном течении вязкой проводящей жидкости со сво бодной поверхностью по наклонной плоскости в одно родных скрещенных электрическом и магнитном полях.
Решение выполнено при условии, что скорость жид кости и индукция наведенного магнитного поля зависят только от координаты у, приложенное электрическое
поле имеет только z-составляющую, градиент давления
= 0 , поскольку на свободной поверхности давление
везде постоянно и равно ро-
В этом случае система уравнений имеет следующий
вид: |
, |
. |
|
j-. г\\ |
|
г> d By |
|
|||
d2 v, |
|
1 |
(88) |
|||||||
f - = — |
(Р g sin a — a Ez B ) ---- В |
dy |
||||||||
dy2 |
|
d2Br = — B |
|
dvx |
|
|
||||
|
1 |
|
|
|
(89) |
|||||
a ji |
dy2 |
|
|
dy |
|
|
|
|||
|
|
|
|
p g cos a -j- a Bx Ez . |
(90) |
|||||
Граничные условия: |
|
|
|
|
|
|
|
|
||
vx = 0 при y = 0 ; |
|
dy |
= 0 |
при у = h; Bx = 0 при у —0; |
||||||
|
|
|
|
|
|
|
|
|
||
d Bx |
0 |
при у = h. |
|
|
|
|||||
|
dy |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
||
Вводя новые переменные |
|
|
|
|
|
|
|
|||
v = vx + |
|
B r |
|
w = vx — |
|
By- |
(91) |
|||
fi у a T) |
p |
y o i| |
||||||||
|
|
|
|
|
||||||
и безразмерные величины |
|
|
|
|
|
|
||||
v = ^ L |
v ; |
w ■ |
|
|
— |
|
У |
(92) |
||
nh2 |
W, y = -*-, |
|||||||||
nh2 |
' |
|
|
|
|
h |
|
|||
где |
|
(p g sin a — a EzB) , |
|
|
|
|||||
n = |
|
|
|
|||||||
получим систему |
двух |
линейных |
уравнении |
второго |
||||||
порядка, решение которой дает: |
|
|
|
|
|
|||||
— _ ch Н а — ch Н а (1 — у) |
|
|
(93) |
|||||||
°х ~~ |
|
|
НИ2 |
|
|
’ |
|
|||
|
|
|
|
|
|
|||||
д _ sh Я a — sh Я а (1— у) — Н а у |
(94) |
|||||||||
х |
|
|
|
Т а 2 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|||
Некоторые результаты теоретических и эксперимен тальных исследований турбулентного течения .проводя-
щей жидкое™ в канале со свободной поверхностью при наличии скрещенных электрического и магнитного полей рассмотрены в работах [24, 25].
В экспериментах |
использовали 30%-ный раствор |
H2SO4. Магнитная индукция в зазоре электромагнита |
|
составляла до 1,3 Т; |
плотность тока менялась в пределах |
О—1,4-103 А/м2. МГД-.поток имел следующие параметры:
R e = (404-60) • Ю3; |
# а = |
164-20; Яа2/Д е= .(54-8) .(Щ-3. |
Исследовали два |
случая |
ориентации электромагнитной |
силы: по направлению силы тяжести и по направлению к свободной поверхности потока.
Установлено, что ориентация электромагнитной силы оказывает значительное влияние на характер деформа ции профиля скорости. В случае, когда электромагнит ная сила совпадает по направлению с силой тяжести, происходит подавление вертикальных и усиление гори зонтальных пульсаций скорости.
Особый интерес представляет исследование влияния скрещенных электрического и магнитного полей на гид равлические характеристики многофазных потоков.
1 |
Рис. 10. Схема экспериментальной МГД-установки: |
||||
— гидродинамические |
решетки; |
2 — канал; |
3 — электроды; |
||
4 |
— устройство |
для послойного отбора проб; |
5 — перегород |
||
ка; 6 — сливной |
бак; |
7 — мерные |
емкости; 8 — насос; 9 — |
||
|
|
|
манометр |
|
|
Экспериментальные исследования профиля осредненной скорости и концентрации двухфазного потока в кана ле МГД-сепаратора, проведенные при различных значе ниях электромагнитной силы, крупности, плотности и концентрации частиц и средней скорости потока, описа ны в работах [26, 27].
Схема экспериментальной установки приведена на рис. 10. В качестве рабочей жидкости использован 22 %-
51
ный раствор |
NaCl (плотность при |
20°С составляет |
1,16-103 иг/м3), |
в качестве дисперсной |
фазы — уголь со |
средним диаметром 0,1—0,3 мм, 0,3—0,7 мм и 1,0—2,5 мм (плотность равна 1,25-103 юг/м3) и кварцевый песок средним диаметром 0,1—0,3 мм (плотность при комнат ной температуре 2 ,6 - 1 0 3 кг/м3).
Требуемая средняя концентрация смеси (0—3% по объему) создавалась засыпкой необходимого количества твердого материала непосредственно в сливной бак ус тановки.
Экспериментами установлено, что наличие взвешен ных частиц в потоке оказывает существенное влияние на профиль средней скорости. На .рис. 11 показаны профили
„ |
— |
и |
„ — с |
|
относительных скоростей |
и = |
— и концентрации |
с = — |
|
|
|
уср |
|
Сср |
Рис. |
11. |
Профили средних скоростей (а) |
н концентраций |
||
(б) |
по |
высоте потока (уголь |
d ср0,3*т*0,7 |
мм: v ср = |
|
=0,7 м/с) |
при уэм> |
Н/м3: |
|
|
|
|
|
1 — 0; |
2 — 0,55-104: |
3 — 10* Н/м3 |
|
(с — расходная концентрация частиц в слое, сер — сред няя расходная концентрация частиц) частиц угля во взвесенесущем потоке в обычных условиях и при нали чии окрещенных электрического и магнитного полей. Как видно, на деформацию профиля скорости значительно влияет перераспределение взвешенных частиц по высоте потока, вызванное присутствием электромагнитных сил в
52
электролите. Это влияние проявляется в подтормаживании слоев жидкости, в которых имеет место наибольшая концентрация частиц.
•*
Рнс. 12. Влияние средней скорости |
Рнс. |
13. |
Распределение |
частиц раз |
|||||||
потока |
на профиль |
концентрации |
личной |
плотности |
(иср= |
*'^ |
У‘‘С' |
||||
(уголь |
dср = 0 ,3 -f 0,7 |
мм; у эм = |
|||||||||
d |
=0,1—0,3 |
мм) при |
V |
, |
Н/м3: |
||||||
= 10* |
Н/м3) при t>cp , м/с: |
ср |
(песок); |
v |
|
эм |
|
3— |
|||
|
/—0,7; 2—1,2; |
3—1,5 |
/—0 |
2—2-10* |
(песок); |
||||||
|
|
|
2-10* (уголь) |
|
|
|
|||||
Распределение взвешен ных в потоке частиц по вы соте зависит не только от величины электромагнитной силы, но и от скорости по тока, плотности и размеров частиц.
Из приведенных на рис. 1 2 кривых можно сделать
вывод, что увеличение сред ней скорости потока спо собствует выравниванию профилей осредненных кон центраций частиц. Такое явление, по-видимому, мо жет быть объяснено ростом интенсивности турбулент ных пульсаций с увеличе нием средней скорости по тока. Кроме того, опре-
/7 J-z
У
0 ----------------------------- -
0,50 х |
!.П |
С |
|
|
O'- |
|
|
Рнс, |
14. Распределение частиц раз |
||
ных |
размеров |
(уголь |
при и = |
= 1,2 м/с и уэм =Ю* Н/м3) при
d ср , мм: /—0,3—0,7; 2 — 1,0—2,5
53
деленное влияние оказывает уменьшение продолжи тельности движения взвешенных частиц через рабочую зону вследствие электромагнитного утяжеления.
Это предположение подтверждается зависимостями, изображенными на рис. 13. Из приведенных кривых так же следует, что, хотя подъемная сила достаточна для всплывания частиц песка, в статическом состоянии пес чинки не успевают перераспределиться по сечению по тока. Более легкие частицы угля всплывают значительно быстрее.
■Существенное влияние на поведение взвешенных ча стиц в потоке оказывает их размер. Кривые на рис. 14 показывают, что с увеличением размеров частиц нерав номерность распределения их по сечению потока значи тельно возрастает. Такое явление можно объяснить увеличением 'гидравлической крупности частиц, т. е. уменьшением влияния вязкости несущей среды.
3. ТЕЧЕНИЯ В БЕГУЩЕМ МАГНИТНОМ ПОЛЕ
Создание бегущего магнитного поля [28]. Бегущее ■магнитное ноле создается системой проводников, кото рые определенным образом включены в сеть многофаз ного переменного тока. При этом токи, протекающие по каждым двум соседним проводникам, смещены по фазе. Магнитное поле, наведенное такой системой токов, сину соидально изменяется во времени. Вдоль поверхности ин дуктора поле этих токов изменяется в пространстве так же по закону, близкому к синусоидальному. Длина полу волны этой синусоиды называется полюсным шагом, или полюсным делением х. В расчетах часто также использу
ют величину а = — . Более подробно механизм образо
вания бегущего магнитного поля рассмотрен в рабо те [29].
Для увеличения создаваемого магнитного поля обмот ку укладывают в пазы магнитолровода, набираемого из листов трансформаторной стали.
Существуют индукторы плоского и цилиндрического типа. Этими индукторами создаются соответственно плоские и цилиндрические поля. В металлургии могут быть применены устройства как с односторонним распо ложением индуктора (например, электромагнитный же лоб), так и с двусторонним (индукционные насосы, элек
54
тромагнитные устройства для обработки отливок в про цессе затвердевания и т. д.). При двустороннем располо жении, как правило, индукторы имеют одинаковый по люсный шаг и одинаковую токовую нагрузку. Подключе нием обмоток к разным фазам или перемещением одного индуктора относительно другого можно добиться опреде ленного угла сдвига между токами в индукторах (ф). В МГД-усдройствах транспортного назначения осуществля ется так называемое согласное включение сторон индук тора. В этом случае токи в противолежащих проводни ках текут в одном направлении (ф = 0 ).
При перемешивании жидкого металла может быть ис пользована другая схема включения сторон индуктора, когда токи протекают в противоположных направлени ях (ф=180°). Такая схема включения называется встреч ным включением.
В зависимости от назначения МГД-устройства можно обеспечить большую или меньшую неравномерность рас пределения электромагнитных сил, а следовательно, раз личную интенсивность перемешивания жидкого металла в канале.
Задачу о стационарном ламинарном течении жидкого
металла в канале плоского индукционного насоса |
(ши |
рина канала много больше его высоты) с учетом |
зату |
хания электромагнитного поля рассмотрел Н. М. |
Охре- |
менко [30]. При решении предполагалось, что существу ет лишь х — составляющая вектора скорости, у — состав
ляющая плотности тока и две составляющие магнитной индукции — Вх и Вг (рис. 15).
Распределение окорости в этом случае может быть определено как
ch (Р На) — ch (Р На
U = U„ |
ch (Р Н а) — 1 |
X |
|
||
|
|
|
|||
X |
ch 2 ср z — 1 |
1— cos 2 ф z |
ГсЬ2ф£> — 1 |
|
|
с2 — ф2 |
|
с1 -f- ф2 |
[ с2 — ф2 |
|
|
|
|
|
|||
|
1— cos 2 ф Ь] |
ch ^Р Н а |
— J — 1 |
(95) |
|
|
|
|
|||
|
с2+ ф2 |
J |
ch (Р На) — I |
||
|
|
||||
|
Р На |
о г |
|
|
|
где с = --- -■ ; Vf='2f%—скорость магнитного поля; |
|
||||
Umax— скорость в центре потока;
55
Р = ' / c h 2<f 6 + c<«2 ,t6 ; Н а ~ У Т У Т ; i|) = a/i;
m = “jTy" l / " 1 + l^ 1 + ео;
" = W / - > + / > + 3 ;
Ha рис. 16 показаны рассчитанные профили скорости для различных чисел Гартмана. Как видно, неравномер ность в распределении электромагнитных сил приводит к появлению точек перегиба в профиле скорости. Наличие
Рис. 15. Схема |
индукционной |
Рис. 16. Профили скорости в Оегу- |
МГД*машнны с двусторонним ин- |
щем магнитном поле при разлнч- |
|
дуктором |
ных числах Гартмана |
|
этого перегиба дает основание предполагать, что устой чивость ламинарного течения в бегущем магнитном поле значительно ниже, чем устойчивость ламинарного тече ния проводящей жидкости в постоянном магнитном поле.
Выражение профиля скорости для индукционного на соса, работающего в режиме перемешивания или элек тромагнитного стопора, имеет следующий вид:
U |
с2V/ I ch 2 ф z— |
1 |
1— cos2трz |
|
|
|
|||
2 |
\ |
с2— ф2 |
|
с2-j-i|)2 |
|
|
|
||
_ 'ch 2 ф 6— |
1 |
1— cos2 г|)b1 ch (Р Я a z/b) — |
1 |
||||||
с2— ф2 |
|
c2 |
+ i|52 |
J |
ch (Р Я а) — 1 |
|
|||
_ 7 |
(sh 2ф 6/2ф Ь) — |
1 |
1— sin 2 op6/2г|>Ь\ |
|
|
||||
Д |
с2— ф2 |
|
|
c2+ il>2’ |
J |
|
|
||
_/ch 2ф Ь— |
1 __ |
1— cos i|)b \ |
sh (Р Я а) |
|
1 |
^ |
|||
V с2— ф2 |
с2+ яр3 / ch (р Я а)— 1J |
|
|||||||
56
ch (P Ha) — ch |
(|3 Ha z/b) |
|
|
(96) |
||
X -ch ф Ha) — sh (P На) |
ф H a)~l |
|
|
|||
В случае развития ламинарного течения проводящей |
||||||
жидкости для осредненной скорости имеем |
|
|
||||
|
|
|
z |
00 |
|
|
|
[ На — |
|
|
|||
Охгп = и о |
V |
ь } |
0 |
|
|
|
ch Н а — На:f-1 sh На |
2 х |
|
|
|||
ср |
|
|
||||
|
|
|
|
т=\ |
|
|
cos ут —cos ут |
( |
На? + ут . . х |
|
|
||
X (Н О2 + Ут) |
ехр |
) |
. (97) |
|||
COS Ут |
V |
Re |
Ь |
1 |
||
где «о — значение |
скорости в |
начальном |
сечении |
(ио= |
||
= oon st); |
|
|
|
|
|
|
у-т— положительные корни уравнения tg z — z.
Как видно, профиль осредиенной окорости не зависит от частоты. Легко показать, что на 'бесконечном-удале
нии от входа выражение |
(63) |
имеет |
вид, полученный |
|||||||
Гартмапом. |
|
|
|
|
|
|
|
|
||
Длина начального участка |
|
|
|
|
||||||
^ _ |
|
b Re |
j |
200 |
(ch На — sh Н а/Н a) |
(cos Ya — 1) |
mg\ |
|||
|
На-- \ - у\ |
|
|
(//a 3 + Yi) (ch Н а — 1) |
cosy! |
|
||||
где yi=4,493 — наименьший отличный от |
нуля |
корень |
||||||||
Из |
|
|
уравнения tg z — z. |
|
|
|
||||
представленного |
на |
рис. |
|
|
|
|||||
17 графика видно, что длина на |
|
|
|
|||||||
чального участка L/L0 уменьша |
|
|
|
|||||||
ется |
с |
увеличением |
магнитного |
|
|
|
||||
поля |
(здесь L0— длина |
началь |
|
|
|
|||||
ного |
участка в отсутствие |
маг |
|
|
|
|||||
нитного поля). |
1 |
и |
ламинарном |
|
|
|
||||
При |
Я а » |
|
|
|
||||||
течении LmbRejHa. |
|
|
|
|
|
|
||||
Турбулентное течение [30]. Те |
|
|
|
|||||||
чения жидких |
|
металлов |
в прак |
|
|
|
||||
тических условиях носят разви
тый турбулентный характер, по |
Рис. |
17. Влияние |
бегущего |
|
магнитного поля |
на |
длину |
||
этому изучение основных законо |
|
начального участка |
||
мерностей этого течения в бегу |
|
прикладное |
зна |
|
щем магнитном поле имеет большое |
||||
чение. Поскольку решение уравнения |
Навье — Стокса |
|||
с учетом турбулентных пульсаций |
скорости и электро |
|||
57
магнитного поля невозможно, возрастает роль экспе риментальных методов исследования. Бегущее магнит ное поле оказывает более сильное влияние «а 'коэффи циент сопротивления, чем постоянное магнитное поле.
На основании опытных данных Н. М. Охремешко предложена следующая зависимость для коэффициента сопротивления в бегущем магнитном поле:
Re (1 — s)
где %о — коэффициент сопротивления при отсутствии маг
нитного поля;
koc — коэффициент ослабления, учитывающий влияние
краевых эффектов и индуцированных в стен ках канала токов;
s — скольжение.
'При значениях числа Стюарта S /= 0 -M 0 - 2 коэффи циенты d равны: Ci=Л2 ,0 ; С2 = 4 ,0 , с3 =П,5.
Коэффициент соответствует начальному участку, на котором происходит основная потеря давления. Коэф фициент с2 относится к потерям давления для активной зоны канала, с3— для зоны стабилизированного течения.
При гидравлически гладком плоском канале и отсут ствии поперечного краевого эффекта имеем
'■ = > • • + <100>
Экспериментально установлено, что под влиянием бе гущего магнитного поля профиль скорости по высоте ка нала становится более плоским, и градиент скорости у стенок канала возрастает . При этом указанные явления проявляются сильнее с увеличением числа Стюарта и в значительной степени зависят от степени проявления по перечного краевого эффекта.
По ширине канала профиль скорости приобретает М.- образную форму, аналогичную кривой распределения ин дукции магнитного поля fio ширине канала.
Незначительное смещение канала от оси индукторов приводит к существенной деформации профиля скорости. Это явление служит причиной возникновения вторичных течений и потерь давления, развиваемого насосом.
58
4. УСТОЙЧИВОСТЬ МГД-ТЕЧЕНИЙ В КАНАЛАХ
Влияние поперечного магнитного поля [31]. Одеон из задач магнитной гидродинамики является задача о влия
нии |
магнитного |
поля |
на |
* |
|
|
устойчивость ламинарного |
те- |
/ |
И1 |
|||
чения. Исследования устойчи |
|
|
||||
вости основываются на методе |
|
|
||||
малых возмущений. |
|
|
|
\В |
||
Рассмотрим плоскопарал |
|
Рис. 18. Течение проводящей |
||||
лельное |
установившееся |
те |
|
|||
|
жидкости между двумя парал |
|||||
чение |
проводящей |
жидкости |
|
лельными пластинами в попе |
||
(рис. |
18) |
при следующих |
ус- |
|
речном магнитном поле |
|
|
|
|||||
ловиях:
vx = vx (у) ; vy = vz == 0 ; р = Агх -f- Л + С (у) ;
Ву = В; Вх = Вх (у); В, = 0,
где А\, Л2 и В — постоянные величины; Вх и С — функции только от у.
При наложении на основное течение малых возмуще
ний составляющие скорости, давления и магнитной |
ин |
||||||||
дукции имеют вид: |
; |
Vy-=v'y ; |
о. - |
о '; |
j |
|
|||
|
Vx = |
vx (у) + v |
|
||||||
|
p = AlX + A2 + C (у) + Р'\ |
|
|
j |
(101) |
||||
|
Вх = |
|
Вх .(у) -j- Ьх ; |
Ву — В Аг by] |
Bz = bz , |
j |
|
||
где vx |
, v'y |
, |
p |
, bx, by, bz величины |
возмуще |
||||
|
|
|
|
|
|
ний, которые предпо |
|||
|
|
|
|
|
|
лагаются малыми по |
|||
|
|
|
|
|
|
сравнению с величи |
|||
|
|
|
|
|
|
нами скорости, |
дав |
||
|
|
|
|
|
|
ления и индукции ос |
|||
|
|
|
|
|
|
новного течения. |
|
||
Подставляя выражения (101) в уравнения (1,1), (12), |
|||||||||
(16) и |
(18) |
|
с учетам уравнения |
(13) |
и сохраняя только |
||||
линейные члены, получим следующие уравнения для воз мущений:
dvx |
dvx |
, |
dvx |
1 |
dp' |
|
|
vx --------1—V |
------— -------------- ■— |
|
|||||
|
dt |
dx |
u |
dу |
p |
3dx- |
|
1 |
„ |
dbx |
В |
dbx |
9BX |
|
|
+ p |
Bx |
X |
dy + by |
dy |
|
|
|
|X |
dx |
|
|
|
|||
|
|
1 |
|
(.Bx bx + Bby); |
( 102) |
||
|
|
P j i |
d x |
||||
|
|
|
|
|
|
||
59
