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

книги из ГПНТБ / Бесконтактный контроль потока жидких металлов

..pdf
Скачиваний:
9
Добавлен:
23.10.2023
Размер:
19.29 Mб
Скачать

Г Л А В А 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 тл, v01 м/сек, 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

 

 

где Ар — градиент давлений; — расстояние между неподвижными твердыми плоско­

стями.

Средняя скорость определится как

 

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] в этих полях ожи­ дается переход к турбулентному режиму для заданного числа Гартмана М=іі?03ффУст/ру при значениях числа Re меньших, чем при использовании постоянного поля. ВЭфф — эффективное значе­ ние магнитной индукции в зоне жидкого металла. Это в. первую очередь относится к бегущим и пульсирующим в пространстве магнитным полям. Влияние постоянного магнитного поля можно рассматривать только как нижний предел влияния бегущего маг­ нитного поля, поскольку наличие переменных во времени и пространстве электромагнитных сил в последнем случае может привести к дополнительным возмущениям потока [7].

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

Физические свойства чаще всего используемых жидких метал­ лов рассмотрены в следующем параграфе.

§ 3. ФИЗИЧЕСКИЕ СВОЙСТВА Ж И Д К И Х МЕТАЛЛОВ

Без знания основных сведений по физическим параметрам жидких металлов невозможно успешное проектирование датчи­ ков индукционных МГД-устройств.

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

Щелочные металлы применяют в качестве теплоносителей в некоторых типах ядерных реакторов [9—11] и в контурах для на­ учных исследований. Основные физические свойства наиболее распространенных жидких металлов даны в табл. 1 приложения.

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