Файл: Щербинин Э.В. Струйные течения вязкой жидкости в магнитном поле.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 152

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

При рассмотрении этих рисунков

прежде

всего о б р а щ а е т на

себя внимание то обстоятельство, что при прохождении

электри­

ческого тока,

направление

которого

совпадает с

направлением

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

радиуса,

жидкость

приходит

в

движение,

харак ­

тером напоминающее

вторичное течение, возбуждаемое

вихре­

 

 

 

 

вой

нитью.

 

Причиной

возникновения

 

 

 

 

движения является вихревой

характер

 

 

 

 

электромагнитной

 

силы,

 

образую ­

 

 

 

 

щейся при

взаимодействии

пропускае­

 

 

 

 

мого

радиального

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

тока

 

 

 

 

и

порожденного

им

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

 

 

 

 

магнитного поля. Нетрудно убедиться,

 

 

 

 

что

по

абсолютной

величине

| /^эм | ~

 

 

 

 

 

1

1

 

/ ' п

ч

 

Л

 

 

 

 

 

 

 

 

~

л3 (1+т)

2 )

\ У 1 + Г | 2

11,

 

а

по

на-

 

 

 

 

 

 

 

 

 

\ у і

+

 

 

 

 

 

 

15 30

45

SO 75

 

правлению

 

она

ортогональна

по­

Рис. 2.21.

Зависимость

ин­

верхностям

r)=const.

В

то

ж е

время

направление

• rot F3 m

совпадает

с

на­

тенсивности

течения от

уг

правлением

 

азимута,

следовательно,

ла конусности воронки.

 

 

 

в

жидкости

 

д о л ж н а возникать

завих­

 

 

 

 

 

 

 

 

 

ренность с осью, направленной по

 

 

 

 

азимуту (см. рис. 2.18).

 

 

 

 

 

По характеру возникновения электромагнитной силы описан­

ное явление сходно с пинч-эффектом

и может быть названо ра­

диальным

пинчем.

Однако

в

отличие,

например,

от

линейного,

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

плотность

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

тока

в

различных

поперечных сече­

ниях

жидкого столба

становится

переменной,

д о л ж н о

возникать

внутреннее

движение

жидкости

в

столбе

аналогично

описан­

ному выше.

 

 

 

 

 

 

 

 

 

В

частном случае

а = 0 решение

(2.118). описывает

«слабый»

смерч

у плоской поверхности,

сопровождаемый коническим раз ­

рядом . В

отсутствие

р а з р я д а

(5 = 0)

это

решение переходит в

найденное

Гольдштиком

[17]1 .

Д р у г о е

частное решение (Re = 0,

а = 0), вытекающее из (2.118), получил Лундквист [30].

Линеаризованное решение задачи (2.118) показывает, что ин­ тенсивность возбуждаемого при пропускании тока движения пропорциональна параметру S, т. е. к в а д р а т у пропускаемого тока. Представляет интерес определение границ изменения

1

Численное решение задачи Гольдштика недавно воспроизведено в ра­

боте

[39].



п а р а м е т р а 5, в которых эта зависимость сохраняется. С этой целью была предпринята попытка численного интегрирования на

Э В М

уравнения

(2.84), в котором

п р а в а я

часть определялась

ре­

шением

(2.115). Результаты

этого

расчета при сс = 0, Re = 0

при­

ведены на рис. 2.22, причем за интенсивность

течения

принима­

лось

значение осевой

скорости

на

оси г = 0,

 

 

 

 

 

 

т. е.

функция

г) (гц|/—я|з).

 

Ка к

видно

из

 

 

 

 

 

 

рисунка,

пропорциональная

зависимость

 

""""

 

 

 

между

интенсивностью

и

 

параметром

S

 

 

 

 

 

 

имеет

место

примерно

до 5 ~ 1 0 . С

увеличе­

 

 

 

 

 

 

нием S интенсивность течения резко воз­

 

 

 

 

 

 

растает,

более

того,

при

5 ^ 1 5 0

 

решение

 

 

 

 

 

 

становится

неограниченным

 

(этот

ж е

ре­

 

 

 

 

 

 

з у л ь т а т

получен

и в работе [31]).

 

 

 

 

 

 

 

 

 

 

В реальных условиях параметр S может

 

 

 

 

 

 

достигать значений порядка 106,

 

та к

что

 

 

 

 

 

 

расчет реальных устройств по методике,

 

 

 

 

 

 

описанной

выше,

оказывается

непригод­

Рис.

2.22.

Расчет ин­

ным. По-видимому, при больших

 

S

наве­

 

тенсивности

течения

денное движением магнитное поле стано­

без

учета

влияния ин­

вится

столь

существенным,

что

пренебре­

дуцированных

токов.

гать

им

в

уравнении

 

(2.84)

у ж е

нельзя.

 

 

 

 

 

 

Действительно, пусть в первом приближе ­

 

 

 

 

 

 

нии i}>»S. Тогда, согласно

(2.89), L i т а к ж е можно

принять

про­

порциональным

S.

В свою

очередь,

в

правой

части

(2.84)

член,

с о д е р ж а щ и й

L]L0,

возрастает,

ка к 5 2 р , т. е. у ж е квадратично по S.

Т а к и м образом,

д а ж е

при

очень м а л ы х

р, но

при

таких

S,

что

« S p ~ l , влияние

индуцируемых

токов

становится

по

порядку

ве­

личины равным влиянию пропускаемого тока. Влияние

ж е

инду­

цируемых токов, ка к будет

 

показано

ниже, сводится

к подтор-

м а ж и в а н и ю

течения,

поэтому

м о ж е т

оказаться,

что-ограничен­

ное решение существует, если решать совместно уравнения

(2.84)

и (2.89).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Влияние индуцированных токов определим на частном при­

мере /4 = 0, В=\.

В этом

случае магнитное поле создается

током,

проходящим по оси симметрии, и индуцированные токи стано­

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

 

В = \, из (2.118)

 

Действительно,

полагая

/ 4 = 0 ,

(для прос­

тоты выберем <х = 0)

имеем

в ы р а ж е н и е

 

Re2

 

 

Arsh г)

 

^ = Ф о = — т г - [ ( 1 - 2 І П 2 )

Ц

 

 

 

 

 

 

~~2~

У і + г , 2

l/1+rv2

 

+ 2У1 + ті2

Arsh-п—-л In ( 1 + т у 2 ) ]

(2.120)

6 — 2274


в

которое

не

входят

члены,

пропорциональные 5, что

связано

с

потенциальным

характером

магнитной

силы. Вихревая

часть

электромагнитной

силы

в

 

этом

случае

учитывается

членом L \

в

разложении

L — L0+^Ll+...

 

 

 

,

она

может

быть

определена

после

решения

уравнения

 

(2.89),

в котором

нужно

п о л о ж и т ь

То = 0,

L Q = 1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, определяющей системой уравнений при ма­

лых р будут служить уравнение

(2.84), в котором член,

пропор­

циональный

5,

следует

исключить,

и

(2.89).

Решение

системы

(2.84), (2.89) представим в

 

виде

р я д а

по

степеням малого

пара ­

метра

Н а 2 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ф ^ ф о + Н а ^ ^

 

 

£ [ = ! , < , + H a 2 L u .

 

 

 

 

 

 

 

Нетрудно убедиться, что решением д л я

г|)0

является

в ы р а ж е н и е

(2.120), а д л я определения

x\i{

и

£ ю с л у ж а т

уравнения

 

 

(1 + ц2)ЧЛ

+ гіфі =

2г,У 1 +

г,2

/

-

^

= r

L j o d t j - 2 ( 1 + 2г|2 )

/ T ) L , 0 A I ;

 

 

 

 

 

 

 

 

 

о

Я + П 2

 

 

 

 

о

(2.121)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( l + T) 2 )L"„>+3r,L' 1 ( ) =

- 2 ^ о

 

 

 

 

 

 

 

 

 

 

.(2.122)

с условиями

г|)і (0)

и

(2.116).

 

 

 

 

 

 

 

 

 

 

 

 

Решением

(2.122)

является

функция

 

 

 

 

 

 

 

 

 

 

 

 

 

(1—2

In 2)

 

( 4 2

+

l n ( l + i i 2

) )

+ ( l + 2 M 2 ) A r s h T 1

+

 

 

Re2

 

 

 

 

 

 

 

 

+

Arshr|

In

( l + i i 2 )

 

 

 

 

 

 

(iW)

 

 

Г| + Г| ІП (1

+Г|2 )

 

( l + ч2 ) 3/2

J

У І + г , 2

 

 

 

 

 

dv\ +

 

 

 

 

1 + T 1 2

 

 

 

CAr\

 

 

 

 

 

 

 

 

 

У і + т і 2

 

 

 

 

 

 

 

 

где C4 находится из условия

L i 0 ( ° o )

= 0 .

 

 

 

 

Численные значения

функции

Ью

отыскивались

графическим

интегрированием, в связи с чем

вместо

бесконечного

верхнего

предела принималось

значение

т) = 4,3.

В

таком

случае С 4 =

= 0,061026.

Численные

значения

Ью использовались затем д л я

вычисления

правой части

(2.121),

после

чего это

уравнение


р е ш а л о с ь численно методом Рунге — Кутта . Результаты решения

приведены на рис. 2.23, та м ж е показана функция -=-4 . Кє

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

Рис. 2.23. Характеристики

Рис. 2.24.

Схема диффузора, обра-

течения,

возникающего

под

зованного

плоскими стенками,

влиянием

магнитного

поля

 

 

индуцированных токов.

 

 

 

2.3.6.

Течение с линейным источником в круговом конусе. З а ­

д а ч а , к

рассмотрению

которой мы здесь приступаем,

в ы б р а н а

потому, что она непосредственно связана

с изучением

про­

странственного течения в диффузоре, схема

которого

приведена

на рис. 2.24, а. Л ю б о й

реальный диффузор,

образованный

плос­

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

нок, в общем случае

расположенных под некоторым углом

друг

к другу.

В частности,

если дв е

последние

стенки о б р а з у ю т

схо­

д я щ у ю с я

трубу, то мы имеем

дело с ситуацией, часто встречаю­

щейся

при конструировании

переходных

участков

рабочих

ка­

н а л о в МГД - насосов [32, 33]: течение в них расходится

в одной

плоскости и сходится в другой.

 

 

 

 

 

Некоторые аспекты анализа

плоского

течения

(при

наличии

л и ш ь двух расходящихся стенок

(см. рис. 2.24,6))

— проблема

Г а м е л я

— были рассмотрены

в п. 2.2. Здесь мы поставим проти­

воположную задачу: попытаемся выяснить роль стенок, ограни­

чивающих течение по координате z, в отсутствие двух

осталь­

ных. К р о м е того, будем считать, что расстояние между

этими

стенками достаточно велико, чтобы можно было пренебречь их взаимным влиянием, а сами стенки представляют собой поверх-