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

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

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

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

Добавлен: 11.04.2024

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

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

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

В итоге имеем систему граничных условий

/ ( о , ) = / ( а 2 ) = 0 ;

 

а3

-)-(Fn-FT2)

= [ / ^ ф = — ,

полностью определяющую задачу .

Введем числа Рейиольдса, Гартмама и магнитное число Рей -

иольдса

и определим

их следующим

образом:

 

 

 

 

 

 

 

 

 

 

pVVm

 

 

Vm

 

 

 

 

 

 

Вновь введенные

величины связаны с п а р а м е т р а м и

 

S u p

соот-

 

 

 

 

Re

 

 

 

 

 

 

 

 

 

 

 

 

ношениями

р = - D

 

и Н а 2 = 5 р . П о л а г а я д а л е е / = Re и,

F=RemK

а і = — -g-, а 2 = + 2-,

получим формулировку

задачи

в

виде

 

« " + 4 « + R e u 2 + 4 H a 2 A . ' - 4 R e „ i - H a 2

Л2 + С 3 = 0;

 

 

 

 

А' + и = 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.32)

 

 

 

 

 

 

а/2

 

 

 

 

 

 

 

 

 

 

« = ( ± | ) = 0 ,

-а/2

 

 

 

я . ( - | ) - х ( - - | ) - ± і .

 

/

И ( Ф ) Л р = ± 1 ,

 

З д е с ь С 3

 

с

2

в

двух последних

в ы р а ж е н и я х знак

« + »

соот-

= =р-, а

 

 

 

Re

 

 

 

 

 

 

 

 

 

 

 

 

 

ветствует

источнику

( д и ф ф у з о р ) ,

знак

« —» — стоку

 

(конфузор) .

З а д а ч а

(2.32)

существенно

упрощается,

если

ограничиться

случаем

 

малой проводимости

жидкости

(Rem <g:l). Пренебрегая

в первом

 

уравнении

(2.32)

членом

с коэффициентом

R e m - H a 2

(малым по сравнению с членом

4 H a 2

V )

и заменяя

 

в нем К' на

— и, получим

уравнение

 

 

 

 

 

 

 

 

 

 

u " + ( 4 - H a 3 ) u + R e u 2

+ C 3 = 0,

 

 

 

 

 

 

 

 

 

(2.33)

порядок дифференцирования которого на единицу

меньше ис­

ходной

системы

(2.32). Соответственно

из

граничных

условий

в (2.32)

существенными остаются л и ш ь первые два :

 

 

 

 

 

 

 

 

 

 

а/2

 

 

 

 

 

 

 

 

 

 

и ( ± —

) = 0

и

 

J и ( ф ) г і Ф = ± 1 .

 

 

 

 

 

 

 

(2.34)

^- а / 2


 

З а м е т и м ,

что безындукционное

приближение

(2.33)

м о ж н о

т а к ж е получить, р а з л а г а я функцию

F в р я д по м а л ы м p<gcl

(или

по

ReT O ): F — F0+fiF[+

. . . . К а к

следует

из (2.30),

решением

д л я

FQ

является

Fo = const,

которую

следует

положить

равной

нулю,

если считать, что внешнее радиальное магнитное поле отсутст­

вует, a Fi будет определяться

уравнением

F'\ = — f = Reu.

П о д ­

становка

этих решений в (2.29)

дает

(2.33).

 

 

2.2.1.

Ползущее течение.

П р и

очень

м а л ы х числах

Re в

(2.33) можно пренебречь квадратичным членом, после чего ре­ шение линейного уравнения

м " + ( 4 - Н а 2 ) и + С 3 = 0

 

 

 

 

 

 

 

 

у ж е не представляет труда.

 

 

 

 

 

 

При

Н а 2 < 4

/

= 4 - Н а 2 ,

 

2 я \

 

 

 

 

со2

 

 

 

 

 

 

 

 

\

 

 

а

 

со /

 

а

 

 

 

 

 

 

 

 

 

,

 

 

 

 

cos com — cos со —-

 

or cos со —

 

и = ± с о

 

 

 

 

 

2

С 3

= +

 

2

 

 

.

а

 

 

 

,

 

 

 

а

0

 

 

 

а

 

а

 

 

 

2 sin со — — с о а cos—

 

2 sin со — — с о а cos со —

При Н а 2 = 4

получаем

 

 

 

 

 

 

" = ± 2 а ^ ( а 2 - 4 ф 2 ) '

 

 

С з = ±

^ '

 

 

 

 

 

т. е. профиль скорости, аналогичный профилю

П у а з е й л я

в плос­

кой или круглой трубе в отсутствие магнитного

поля.

 

 

Если

Н а 2 > 4 ,

то решение

приобретает вид

(со2 = Н а 2 — 4):

 

ch со

— ch шф

 

со3

ch со

 

 

« = ± с о

і.

 

 

п

и

 

, С 3 = ±

 

о

и

 

 

а

 

а

 

и

а

а

 

coachco— — 2shco—

 

с о а с п с о — — г

sn со —

 

 

 

 

 

 

 

 

 

 

 

 

(2.35)

П р и м е р воздействия поля

на

распределение

скорости

в д и ф -

фузоре с углом

а =

 

показан на рис. 2.7 [22]. В предельном слу­

чае больших

чисел

Г а р т м а н а

( Н а » 1 ) профиль

скорости

в я д р е

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


множи ­ порядка
Рис. 2.7. Распределе­ ние скорости в стоксовом режиме тече­ ния.

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

ной на рис. 2.6, то в случае

д и ф ф у з о р а мы имеем значительный

отрицательный

градиент

давления

в отличие от течения в отсут­

ствие поля и при числах

Н а 2 < 4 .

 

 

 

К а к

показано

в работе

А. Б .

Цинобера

Иа--0

[24], приведенное

выше

решение (2.35) дл я

 

медленного

(стоксова)

течения

имеет

 

место не только при м а л ы х Re,

но

и при

 

любых

фиксированных

Re и а, если числа

 

Н а достаточно

велики,

причем

решение не­

 

линейной з а д а ч и

(2.33), (2.34)

равномерно

 

по ф стремится

к линейному

решению

(2.35)

на--ют

при Н а - ^ о о , отличаясь

от него

на

величину

 

Re

порядка

Н а 2

2.2.2. Течение в диффузоре при больших

числах

Рейнольдса. П р е ж д е чем перейти

к анализу течения при больших числах Re,

заметим,

что уравнение (2.33) имеет интегрирующий

тель и', та к что его можно свести к уравнению первого

"'2=4З Re[«-

•u3-

2 Re

) .„

Re

RP J

Я

 

 

3 ( 4 - Н а 2

ЗСз

. J C j l

2

интеграл которого, в свою очередь, есть du

У"

(2.36)

Р а з л о ж и м

 

подкоренное

в ы р а ж е н и е

в (2.36) на простые множи ­

тели

Р(и)

= (еі — и) (и—е2)

(и—е3).

Тогда

(2.36) запишется ка к

 

2 Re

 

 

 

du

 

(2.37)

 

 

 

У(еі~и)

 

( ы - е ? )

{и-е3)

 

 

 

 

 

где

корни

полинома еи

е2,

е3 связаны определенными соотноше­

ниями с коэффициентами

полинома

Р(и),

из которых существен­

ным д л я дальнейшего является соотношение

Єі +

Є2 + Є3

=

3 ( Н а З - 4 )

 

 

 

(2.38)

2 Re

 

 

 

 

 

 

 

 

 

 

Т а к к а к на стенках к а н а л а и = 0 , то в промежутке

аа

Т2


д о л ж н а быть

по крайней

мере одна точка, где и' = 0, т. е. из кор­

ней полинома

по меньшей мере один должен

быть вещественным .

Пусть это будет корень

в\. Остальные корни

е2 и е% могут быть

либо комплексно

сопряженными, либо

вещественными.

 

 

О с т а в л я я в стороне

случай

комплексно сопряженных

корней,

рассмотрим в о з м о ж н ы е

решения при больших

Re в

случае

ве­

щественных корней полинома. Пусть

 

 

 

 

 

 

 

 

 

 

 

 

Н а 2 —4

 

Тогда, учитывая

(2.38),

м о ж н о

показать, что

. Под -

 

 

 

 

 

 

 

I

Re

 

 

коренное

в ы р а ж е н и е

в

(2.37)

д о л ж н о

быть

неотрицательным,

т. е. Р(и)

> 0 . Это возможно лишь в тех случаях, когда

« < е 3

или

в2<и<.Єі.

Ограничимся

 

теперь такими течениями, в которых

дви­

жение жидкости осуществляется в одном направлении, т. е. и

везде

имеет один знак . Т а к ка к в диффузор е

скорость и положи ­

тельна, а на стенках

и = 0, то ясно, что д о л ж н о быть либо

а)

0 ^

 

 

 

Н а 2

—4

 

 

 

 

 

 

 

 

 

=С«Ї^<?З,

^ 5

— , либ о

б)

0 ^ w ^ e b

е 2

^ 0 , е і > 0 .

 

 

Случай

а)

2, І

 

 

 

 

 

 

 

 

 

 

возможен лишь

при условии

Н а 2 > 4 .

Будем

счи­

тать течение симметричным. Тогда

значение ы = е 3

будет

дости­

гаться в точке к р = 0 ,

а (2.36)

перепишется

в виде

 

 

 

-і/

2 Re

Г

du

 

 

 

а

 

 

 

 

 

 

ф у~—=

/

,

 

О ^ Ф ^ — .

 

 

 

 

(2.39)

'

3

{

уР{и)

 

 

 

2

 

 

 

 

v

'

Значения

корней полинома

определятся

из условия

(2.38)

и ус ­

ловий (2.34), которые приобретут вид

 

 

 

 

 

 

 

ез

 

 

 

 

 

е3

 

 

 

 

 

'

6

{

VP (и)

*

6

 

/

1/P(U)

 

 

 

 

Введем в интегралы

(2.39)

и (2.40)

вместо

и новую

переменную

положив

 

 

 

 

 

 

 

 

 

 

 

 

M=e 8 - (ei - <?8)ctg*it> .

 

 

 

 

 

 

 

(2.41)

Значение

« = е 3

будет

Достигаться

при ip = -g-, а значение

и = 0 —

при

і|>=гЬо, где -шо определяется

из c t g 2 ф 0

go

. Соответст-

=

венно условия

(2.40)

перепишутся ка к

 

Є\ — Єз

 

 

 

 

 

 

 

,

 

 

 

л/2

 

 

 

 

 

 

 

 

 

a І

/

[

;

(2.42)


§ 2. ПЛОСКИЕ И ОСЕСИММЕТРИЧНЫЕ ТЕЧЕНИЯ

 

і

47

 

 

 

 

 

 

 

 

я/2

 

я/2

 

 

'

6

y e i - e 9 J y i - ^ s i n ^

 

^

 

 

 

 

- 2 - j / e j - е 3 ctg tpo y 1 - Л 2 sin 2

яро

,

(2.43)

г д е

A2=_fiz5_<i.

 

 

 

 

 

 

е і - е 3

 

 

 

 

 

Большие значения чисел Re приводят к тому, что правые

части формул (2.42), (2.43) становятся

очень большими. В свою

очередь,

эт о означает, что параметр

k д о л ж е н мало

отличаться

от единицы, в противном случае эллиптические интегралы и, со­

ответственно, правые части

(2.42), (2.43)

имели

бы конечное

зна­

чение. Та к ка к

 

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

е25;

кроме

того, это

позволяет

в правой

части

(2.43)

пренебречь вторым

и третьим

слагаемыми,

которые

в отличие

от первого

имеют конечное

зна­

чение. Сравнивая теперь (2.42) и (2.43), получаем

 

 

 

1

 

л/2

 

 

 

 

 

 

 

 

 

 

 

е 3

и

/

 

* "

-

У

 

- Є

 

«

(2.44)

« - L

 

К Є ( Є 1

З )

 

a

 

f y i - £ 2 s i n 2 o b

 

 

'

6

 

2

 

 

 

 

 

 

Фо

'

 

т

 

 

 

 

 

 

 

 

 

Найдем

значения

остальных

корней

полинома. Д л я этой

цели

представим

 

 

 

 

 

 

 

 

 

 

 

 

 

Я/2

 

 

 

Я/2

 

 

 

-фо

 

dip

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— J і

— /в

J

y i - / f e 2 s i n 2 \ p ~ I

y i - A : 2 s i n 2 T | )

/

 

 

 

 

Vl —/г2 sin2 г|з

 

При &~1 будут иметь место приближенные равенства [25]:

т і -

 

/

Фо

^

, n І + sin ч|ю

 

4

Г

 

 

 

 

 

 

 

 

Уі — k2

 

J

COS \Jj

COS1|)0

 

 

= l n ( V f i z £ L +

y M ,

 

подстановка

которых в (2.44)

приводит к соотношению

 

е . - е з

( У g ' ~ g 3 , y g ' ) 2

6

 

 

 

Є з

 

' ' Є з

(2.45)