![](/user_photo/_userpic.png)
книги из ГПНТБ / Бесконтактный контроль потока жидких металлов
..pdfГ Л А В А II
Д В И Ж Е Н ИЕ Ж И Д К О Г О МЕТАЛЛ А В МАГНИТНО М ПОЛЕ
§ 1. МГД-УРАВНЕНИЯ ДЛЯ МЕДЛЕННО ДВИЖУЩИХСЯ ПРОВОДЯЩИХ СРЕД
Прикладные вопросы магнитной гидродинамики — науки о законах взаимодействия магнитного поля и электропроводящей среды — связаны с разработкой и применением различных видов МГД-устройств. К таким устройствам относятся МГД-насосы для перекачивания жидких металлов, МГД-генераторы, расходо меры для измерения расхода электропроводящих жидкостей, плазменные двигатели, дозаторы жидких металлов и т. п.
В МГД-устройствах одновременно протекают как электро магнитные, или электродинамические, так и гидродинамические процессы, причем эти процессы взаимосвязаны.
Электромагнитные процессы описываются уравнениями Макс велла, гидродинамические — уравнением Навье-Стокса. Взаимо связь этих двух процессов обусловливается тем, что в уравнении Навье-Стокса имеется дополнительный член, выражающий объемную электромагнитную силу, действующую на элемент жидкости. Кроме того, в уравнении энергии при решении МГДзадач должно быть учтено выделение джоулева тепла в объеме контролируемой среды.
Жидкий металл можно рассматривать как изотропную вяз кую несжимаемую жидкость [1]. Физические параметры жидких металлов р (плотность), т) (динамическая вязкость), а (электро проводность) считаются не зависящими от напряженностей элек трического и магнитного полей. Для жидких металлов и электро литов практически >д. = |ло, где ц.о = 4я • Ю - 7 гн/м.
Сравнительно большая величина электропроводности жидких металлов позволяет пренебречь в них токами смещения. Рассмат риваемые среды электрически нейтральны, т. е. в них нет объем ных электрических зарядов, поэтому конвективные электрические токи отсутствуют.
Детальное рассмотрение физических процессов является пред метом специальных курсов магнитной гидродинамики. Нас же интересуют лишь основные соотношения и принципиальный под ход к их рассмотрению.
Для вязкой несжимаемой проводящей жидкости с учетом указанных предположений система уравнений магнитной гидро динамики записывается следующим образом [1, 2]:
dv |
1 |
|
1 |
|
— + ( v V ) v = - — V p + v A v + — [ j B ] ; |
(2.1) |
|||
дВ |
|
1 |
ДВ; |
(2.2) |
—— =rot[vB] + |
ЦоО |
|||
Оl |
|
|
|
|
j = — |
rot B = a(E + [vB]); |
(2.3) |
||
1Ю- |
|
|
|
|
|
r o t E = - | L ; |
|
( 2 . 4 ) |
|
|
divv = 0; |
|
|
(2.5) |
|
divB = 0. |
|
|
(2.6) |
Для жидких металлов ввиду небольшой разности температур в зоне измерений можно пренебречь температурными изменени
ями физических параметров. |
|
|
|
|
Запишем уравнение для энергии в указанном |
случае |
|
||
|
РСР (~^r+vVT |
) = ^ + — + № Т |
р , |
(2.7) |
где j2/a |
— мощность, затрачиваемая на джоулево тепло; |
|
||
Wrp |
•— член, характеризующий работу сил трения. |
|
Граничные условия для решения этих уравнений в общем виде являются простым перечнем граничных условий электродина мики, гидродинамики и теории теплопередачи.
В магнитной гидродинамике жидких металлов можно выде лить два основных аспекта [3]:
1) |
изучение влияния |
движущегося металла на |
магнитное |
поле; |
' |
|
|
2) |
изучение влияния |
магнитного поля и связанных |
с ним то |
ков и сил на движение жидкого металла.
Для качественного анализа магнитогидродинамических про цессов, а также численного решения задач на ЭВМ целесооб разно привести уравнения гидро- и электродинамики к безраз мерному виду. При этом уменьшается число переменных и ре зультат может быть обобщен по отношению к целому классу физически подобных явлений с одинаковыми безразмерными па раметрами.
Введем |
следующие безразмерные величины: v = vtVo\ t = |
tJo; |
p=ptpv02; |
B = BtB0; E = E*v0B0; j=j*av0B0; T=T*T0; V = - n i . |
|
AO
При подстановке этих величин в систему уравнений (2.1)— (2.7) получим указанные уравнения в критериальной форме (все величины безразмерные, поэтому * опущена):
I |
dv |
|
у ) v = - |
Vp+—— |
Ї М |
2 |
(2.8) |
St |
at |
( v |
Av + — - [j B]; |
||||
|
|
Rem |
Re |
|
|||
|
|
1 |
дВ |
1 |
|
|
|
|
|
St |
dt |
Re„ - AB + rot [vB]; |
(2.9) |
||
|
|
j = - ^ — r o t B = ( E + [ v B ] ) ; |
|
(2.10) |
|||
|
|
|
|
1 |
c5B |
|
|
|
|
|
divv = divB = 0; |
|
(2.12) |
В теории индукционных расходомеров основными величинами являются индукция магнитного поля и скорость жидкого ме талла. Знание распределения магнитного поля в рабочем зазоре датчика и распределения скоростей по сечению измерительного участка канала позволяет определить основные параметры рас ходомеров. Поэтому основными критериальными параметрами в теории расходомеров являются числа Рейнольдса Re, Rem и число Гартмана М.
Безразмерные критерии и некоторые их комбинации позво ляют оценить степень влияния магнитогидродинамических эф фектов на течение проводящей жидкости, характер переходного процесса при нестационарном движении жидкого металла в ка нале (например, при дозировании), а также влияние джоулева тепла на теплообмен. Последнее может быть существенным при выполнении условия М R e / в ^ 1/Re Рг, что непосредственно сле дует из (2.13).
В табл. 2.1 приведены значения основных критериальных па раметров для жидких металлов при условии -г?0=0,01 м, В0 =
= 0,01 тл, v0— 1 м/сек, f = 100 гц.
§ 2. ОСНОВНЫЕ БЕЗРАЗМЕРНЫЕ КРИТЕРИИ И ИХ ФИЗИЧЕСКИЙ СМЫСЛ
Гидродинамическое число Рейнольдса. Основные критерии подобия можно получить в результате простого анализа, одно временно раскрывающего физический смысл каждого критерия.
|
|
|
|
Численные значения |
критериев |
для, некоторых |
проводящих |
|
|
|
|
|
|
|
|
Ж и д к а я |
|
|
|
Критерій"! |
|
Bi |
Na |
Na56K-l4 |
P b - B i |
|
|
|
|
|
|
|
Т= 300 °С |
|
|
|
|
ч |
|
6,5 • 104 |
2,5610» |
2,8 • 1О4 |
6- 10» |
|
|
|
|
9,55I0~3 |
7,55- Ю-2 |
|
|
|
|
Re |
m |
0 |
- |
2 |
1,07-10- |
2 |
||
|
= u.o-/?t>o |
|
|
2,32 • 10" |
|
|||
M = / ? 0 B o [ / j l |
|
214 |
1320 |
785 |
222 |
|
||
Р г = |
IlfJL. |
|
1,39 • 10~2 |
0,594-10-2 |
1,176- Ю - 2 |
2,3 • Ю-2 |
||
|
|
У. |
|
|
|
|
|
|
fi |
|
R02p2cpT0 |
|
6,2- 1012 |
8,57 • 1012 |
8,421012 |
5,271012 |
|
|
|
Ч2 |
|
|
||||
|
|
|
|
|
|
|
|
|
|
W |
|
0,703 |
68 |
22 |
0,823. |
|
|
|
|
|
|
|
|
|
||
М* Re |
|
4,8- 10-" |
0,52 • Ю-2 |
20,5 • 10-" |
5,61 • Ю-4 |
|||
|
0 |
|
|
|
|
|
|
|
|
|
1 |
|
1,1 • ю - 3 |
6,58- Ю-3 |
3,04 • Ю - 3 |
0,723 • Ю-3 |
|
Re Рг |
|
|
|
|
|
|
||
d> = \ioaa>Ro2 |
|
6,2 • Ю-2 |
0,48 |
0,147 |
6,8 • Ю - 2 |
Рассмотрим, какие силы действуют в жидкости при МГД-тече- нии. В реальной жидкости всегда действуют силы вязкости. В со ответствии с законом Ньютона величина действующих сил вязко
сти будет порядка |
v |
FB~pv-^-. |
(2.14) |
|
Ко |
Другой действующей силой является сила инерции, которая, как нетрудно в этом убедиться, пропорциональна
|
|
F |
p U ° 2 |
|
(2.15) |
|
|
|
До |
|
|
Эти две силы действуют |
при движении |
жидкости |
в отсутствие |
||
электромагнитного поля. |
|
|
|
|
|
Критерий Яе = -~^= 0 |
0 |
— гидродинамическое |
число Рей- |
||
нольдса — равен |
* в |
V |
|
|
|
отношению инерционных сил, действующих в |
|||||
жидкости, к силам |
вязкого |
трения. Этот |
безразмерный показа |
тель получил широкое распространение в гидро- и газодинамике
|
|
|
|
|
|
Т а б л и ц а 2.1 |
|
жидкостей (./?о = 0,01 |
м; Ва—\ ТА; ti0 = 1 м/сек; |
f = 100 гц) |
|
|
|||
среда |
|
|
|
|
|
|
|
Li |
Zn |
|
А1 |
|
H N 0 3 |
NaCI |
Морская |
|
|
(с=25%) |
(с=20%) |
вода |
|||
|
|
|
|
|
|||
Г=500° С |
Г=600° |
с |
Г=700° С |
|
Т=20° С |
|
|
1,36- Ю4 |
3,04 • 104 |
8,2 • 104 |
8,75 • 104 |
9,6-103 |
7,52 • 103 |
9,28103 |
|
3,76 • Ю-3 1,86- Ю - 3 |
4,8 • Ю-3 |
1,33- ю - 2 |
9,67 • Ю - 7 |
2,5 • 10~7 |
0,5 • Ю - 7 |
||
920 |
357 |
|
390 |
262 |
2,54 |
1,15 |
0,6 |
|
|
|
|
2,44 Ю - 2 |
7,23 |
9,6 |
7,25 |
|
|
|
|
10,6- Ю1 2 |
2,96- Ю1 2 |
2,08 • 1012 |
3,37 • 10'2 |
62,23 |
4,18 |
|
1,85 |
0,784 |
0,672 • Ю - 3 0,176- Ю - 3 |
0,388 • Ю-4 |
|
|
|
|
|
5,66 • Ю - 4 |
20,9 • 10~9 |
4,78- 10-э |
ю-9 |
|
|
|
|
0,476 • Ю - 3 |
0,144 • 10-" 0,138- Ю-4 |
0,148- Ю-4 |
|
0,236 |
0,216 |
0,376 |
8 • Ю - 2 |
|
|
|
и является критерием турбулентности. При значении Re^ReKp, когда отношение сил инерции к силам вязкости не превосходит некоторого критического значения, наблюдается упорядоченный режим течения вязкой жидкости (газа), характеризующийся от сутствием перемешивания между соседними слоями жидкости. При достаточно малых скоростях течения случайно возникающие в потоке возмущения гаснут, не вызывая заметного перемешива ния слоев жидкости. С увеличением же скорости в потоке созда ются условия, при которых ламинарное течение перестает быть устойчивым и под влиянием случайных возмущений переходит в турбулентное.
Наличие условий, при которых ламинарное течение перестает быть устойчивым, зависит от гидродинамического числа Рейнольдса. Дл я каждого вида^течения существует такое число ReKp , называемое нижним критическим числом Рейнольдса, когдапри любом Re<ReK p ламинарное течение является устойчивым. При значениях Re>ReK p ламинарное течение при возникновении воз мущений переходит в турбулентное. Например, для течения в
5 — 2939
круглой цилиндрической трубе R e K p « 1100. Это означает, что для воды (v=10~6 м2/сек) устойчивое ламинарное течение со скоро стью v0—\ м/сек возможно лишь в трубках диаметром не более 2,2 мм.
Вследствие вязкого трения жидкость в трубе течет цилиндри ческими слоями, движение которых по мере приближения к стен кам трубы все более замедляется. Распределение скорости по се чению канала можно представить равенством .
° = Ф - ( £ 0 П ' |
І 2 Л 6 ) |
где г — расстояние от оси; и0 = 2/з"ои- При турбулентном течении элементы жидкости совершают не
упорядоченные, неустановившиеся движения по сложным траек ториям. Это приводит к интенсивному перемешиванию слоев жид кости. Профиль усредненной скорости турбулентного течения от личается от параболического профиля ламинарного течения более быстрым возрастанием скорости у стенок и меньшей кри визной в центральной части течения. Для турбулентного течения вязкой несжимаемой жидкости в круглой трубе имеется ряд полуэмпирическнх формул [4]. Например,
0 = 0 м (/ і - |I -г І )V і , (2.17)
где m = 0,14-0,2.
Магнитное число Рейнольдса. Взаимодействие магнитного поля и проводящей среды в случае стационарного течения опре деляется из основных уравнений магнитной гидродинамики, из которых следует, что
rot[vB]+ |
— |
at |
(2.18) |
|
IIQO |
|
|
Известно, что уравнения типа |
|
|
|
!j£=DAu |
(2.19) |
описывают процесс диффузии, где D — коэффициент диффузии (например, при тепловых процессах D — коэффициент темпера туропроводности) .
Уравнение вида
— =rot[vu] |
(2.20) |
может быть приведено в случае v = const к уравнению тепловой конвекции
дТ |
|
— + v grad Т = 0. |
(2.21) |
Таким образом, уравнение (2.18) описывает взаимосвязь двух процессов: первый член определяет конвекцию магнитного поля, второй характеризует диффузию магнитного поля.
Анализ уравнения (2.18), записанного в безразмерном виде, показывает, что критерием преобладания конвекции магнитного поля над его диффузией является магнитное число Рейнольдса
Rem=noavoR0.
Известно, что, например, в вязком потоке вязкость вызывает диффузию вихря вопреки конвекции и обычное число Рейнольдса
R e = ° 0 ^ 0 определяет степень преобладания конвекции над диф-
V
фузией вихря.
В магнитной гидродинамике проводимость может приводить к преобладанию конвекции над диффузией магнитного поля. Дей ствительно, поскольку коэффициент диффузии магнитного поля
п1
D=—, то с увеличением проводимости среды конвекция должна
преобладать.
В индукционных расходомерах жидких металлов определяю щую роль играет диффузия магнитного поля. В этом случае маг нитное число Рейнольдса определяет относительную величину возмущений приложенного магнитного ПОЛЯ .
Легко оценить порядок величины малых возмущений поля. По закону Ома индуцированный в движущейся жидкости ток в об щем случае имеет величину порядка avoB0. Индуцированное поле, согласно соотношению uj„ = rotBH , по порядку величины равно BH=\io<yv0BoRo- Следовательно, отношение индуцированного маг нитного поля Ва к приложенному магнитному полю В0 имеет по рядок
\xoOV0R0='Rem.
Таким образом, значение магнитного числа Rem позволяет оце нить порядок величины индуцированного магнитного поля (вто ричного магнитного поля), а тем самым и порядок величины ско ростного сигнала при заданной частоте, геометрии приемного витка и величине приложенного магнитного поля.
Число Гартмана. Исходя из уравнений магнитной гидродина мики, можно заключить, что магнитное поле оказывает сущест венное влияние на характер течения электропроводящей среды.
При этом важную роль играет взаимное пространственное распо ложение векторов магнитной индукции и скорости проводящей среды. Из общих законов электродинамики следует, что за счет движения токи в электропроводящей среде наводятся только в том случае, если имеется составляющая магнитного поля, пер пендикулярная к направлению перемещения потока.
Наиболее наглядно влияние магнитного поля на течение элек тропроводящей жидкости можно показать на примере задачи Гартмана: найти распределение скоростей в несжимаемой вязкой проводящей жидкости, стационарно движущейся в пространстве между двумя плоскопараллельными твердыми плоскостями; в направлении, препендикулярном к плоскостям, приложено одно
родное внешнее магнитное поле |
В0. |
|
|
||
Как известно [2], решением |
указанной задачи является функ |
||||
ция |
|
|
|
|
|
|
„ . c h M - c h M - |
У |
|
||
v = —Ap' |
М |
\ |
|
Ь |
(2.22) |
|
/ |
shM |
|
||
\ о б о 2 |
|
|
где Ар — градиент давлений; 2Ь — расстояние между неподвижными твердыми плоско
стями.
Средняя скорость определится как |
|
||
v0=-Ap- |
M - t h M |
(2.23) |
|
обо2 th М |
|||
|
откуда вместо (2.22) получаем
chM —chM-
v = v0 M c h M - s h M (2.24)
М/СЕК |
M=0 |
oil |
— 37 |
|
|
|
97,5 |
0J9 |
. 235 |
Oil |
|
0,(5' |
(Б »,ии |
(2 |
Рис. 2.1. Экспериментальные профили скорости течения рту ти в магнитогидродинамическом канале.
На рис. 2.1 показано распределение скоростей ртути при течении в по стоянном поперечном магнитном поле при различных напряженностях последнего. С увеличением на пряженности магнитного поля про исходит уплощение профиля ско рости. Средняя часть потока может считаться однородной и изменение скорости происходит в узком погра ничном слое. Физическое объясне ние может быть следующим: как уменьшение вязкости при уменьше-
нии силы сдвига между слоями выравнивает их скорости, так и увеличение о и В0, вызывая относительное повышение мощности jB0v, развиваемой слоями с большой скоростью, в силу закона о минимуме энергии, способствует выравниванию энергий слоев. Одновременно с выравниванием энергий слоев происходит и вы равнивание их скоростей [5].
Таким образом, критерием влияния магнитного поля на тече
ние жидкости оказывается отношение |
|
|
В0Ь » ] / — |
= М , |
(2.25) |
' ц |
|
|
характеризующее степень деформации профиля скорости. При
малых значениях числа Гартмана М (М<с1) |
получаем v = vM(l — |
|
— y2/b2) — ламинарное течение |
в соответствии с результатами |
|
обычной гидродинамики. Если |
1, то |
|
v = vu(\-e-(b-\v\)№) |
(2.26) |
|
и |
|
|
и 0 = - А р |
= - . |
(2.27) |
Следовательно, увеличение магнитного поля делает профиль ско ростей более плоским на большой части сечения и уменьшает среднюю скорость движения при заданном градиенте давления. При достаточно больших значениях М образуются пограничные слои жидкости с экспоненциальным характером убывания ско рости.
В гидродинамике для оценки толщины пограничного слоя су ществует понятие о толщине вытеснения, т. е. расстоянии, на ко тором можно усреднить отклонения скоростей вблизи стенки от скорости основного потока v0 [6]:
|
Ьі |
|
б = — |
[ (v0-v)dy=~^M. |
(2.28) |
v0J |
4 u |
M |
Для достаточно больших значений М, при которых t h M » l , от носительная толщина пограничного слоя при ламинарном тече нии в магнитном поле имеет порядок М - 1 :
8,= А = М - ' . |
(2.29) |
Ь\ |
|
Таким образом, число Гартмана М — это обратная величина от носительной толщины пограничного слоя в сильном магнитном поле. Более наглядную трактовку числа Гартмана можно дать.
если рассматривать действие сил на объем жидкости, движу щейся в магнитном поле. В движущемся объеме жидкости
разность |
потенциалов, наведенная в нем магнитным полем, |
А ф ~ BQVQRO. |
Разность потенциалов обусловливает электрический |
ток в жидкости ,
/— OBQVQ.
Всвою очередь, ток взаимодействует с создавшим его магнитным
полем В0, и на элементарный объем жидкости действует электро магнитная сила, порядок величины которой
F 8 M = / f i ~ o B o 4 . |
(2.30) |
С другой стороны, при течении жидкости существенную роль иг рают силы вязкости
Итак, для оценки роли влияния электромагнитной силы на те чение электропроводящей жидкости получаем следующий кри терий:
F3JFB = B<?Ro2o/pv.
Нетрудно увидеть, что квадратный корень из этой величины есть не что иное, как число Гартмана
' ' Г в |
' PV |
Таким образом, число Гартмана характеризует отношение элек тромагнитных сил, действующих в объеме электропроводящей жидкости, к силам вязкости, действующим в том же объеме жидкости.
Если при течении жидкости существенны также силы инер ции, то для оценки роли электромагнитной силы получаем
Fun pv0
Этот критерий, называемый числом Стюарта, равен
М 2
Г . - — . .(2.32)
Для жидких металлов практически М^>1. При этом мерой вли яния магнитного поля на течение проводящей жидкости, как это. следует из (2.8), является параметр FS[Q].
Если М > - 1 , a M2 /Re<§;l, электромагнитные силы малы по сравнению с силами инерции и магнитное поле практически не влияет иа течение жидкости. -
При |
М » 1 и M 2 / R e » l магнитное поле |
оказывает существен |
|
ное влияние на течение электропроводящей |
жидкости. |
||
При |
движении проводящей жидкости |
в |
постоянном магнит |
ном поле характерно увеличение устойчивости потока. Наложе нием магнитного поля в принципе может быть ламинаризовано большинство потоков жидкого металла, встречающихся в прак тике [5]. Более того, переход к турбулентному режиму носит бо лее плавный характер и появляется значительная переходная область.
Поскольку индукционные МГД-устройства в основном ис пользуют переменные магнитные поля, благодаря наличию воз мущающего действия вихревого поля сил (7, 8] в этих полях ожи дается переход к турбулентному режиму для заданного числа Гартмана М=іі?0-В3ффУст/ру при значениях числа Re меньших, чем при использовании постоянного поля. ВЭфф — эффективное значе ние магнитной индукции в зоне жидкого металла. Это в. первую очередь относится к бегущим и пульсирующим в пространстве магнитным полям. Влияние постоянного магнитного поля можно рассматривать только как нижний предел влияния бегущего маг нитного поля, поскольку наличие переменных во времени и пространстве электромагнитных сил в последнем случае может привести к дополнительным возмущениям потока [7].
Для вычисления безразмерных критериев необходимо распо лагать физическими характеристиками жидких металлов, пара метрами приложенных электромагнитных полей, условиями тече ния металла.
Физические свойства чаще всего используемых жидких метал лов рассмотрены в следующем параграфе.
§ 3. ФИЗИЧЕСКИЕ СВОЙСТВА Ж И Д К И Х МЕТАЛЛОВ
Без знания основных сведений по физическим параметрам жидких металлов невозможно успешное проектирование датчи ков индукционных МГД-устройств.
Металлургическая промышленность и литейное производство имеют дело с разнообразными жидкими металлами и их спла вами.
Щелочные металлы применяют в качестве теплоносителей в некоторых типах ядерных реакторов [9—11] и в контурах для на учных исследований. Основные физические свойства наиболее распространенных жидких металлов даны в табл. 1 приложения.