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

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

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

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

Добавлен: 11.04.2024

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

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

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

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

Уравнение индукции выведем, используя уравнения (1)(3)

Максвелла и обобщенный закон Ома (8) (можно использовать

урав­

нение закона Ома с учетом эффекта

Холла

и скольжения

ионов,

но это усложнит рассмотрение вопроса).

 

 

Исключив электрический ток

из этих

уравнений, имеем

r o t # = o-£-j-o-(â X В).

Найдем Е и подставим в уравнение (2), в результате получим искомое уравнение для магнитной индукции:

Ц- =

rot X Ъ) — rot

rot H j

 

или

 

 

 

^ - =

rot(üXß) r o t ( v m r o t ß ) ,

(93)

Отметим, что электрическое поле

Е, существующее в

среде,

не входит в уравнение (93) индукции,

вследствие принятых огра­

ничений.

 

 

 

Если рассматривать случай, когда

ѵ,„ = const, то ѵ,„

можно

вынести за знак операции rot. Дважды примененная к вектору В операция rot может быть представлена в виде

rot (rot В) = grad (div ß) Aß.

С учетом уравнения (4) уравнение для магнитной индукции принимает вид]

^ = r o t ( ö x ß ) + vm Aß.

(94)

Точно такую же структуру имеет обобщенное уравнение Гельмгольца

^ = rot(uxcû) + vAcû

(95)

для гидродинамического вихря скорости со = - у rot ѵ при дви­ жении жидкости, имеющей кинематическую вязкость ѵ, когда: 1) массовые силы fm имеют потенциал; 2) процесс в жидкости баротропный, т. е. когда плотность является функцией давления

Р = / (Р).

108


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

нения (94), будут описывать свойство не только

вектора

индук­

ции В магнитного

поля,

но вектора со гидродинамического

вихря

скорости.

 

 

 

 

 

 

Величина \ т

1/оцт а

является своеобразным

аналогом коэф­

фициента кинематической

вязкости ѵ и поэтому ее часто называют

«магнитной» вязкостью.

 

 

 

 

Уравнения (93) и (94) выявляют причины, вызывающие изме­

нение вектора В магнитного поля

и вектора со вихря в

рассма­

триваемой

точке пространства, а

именно:

 

 

1) первый член правой части этих уравнений отражает изме­

нение поля

В и вихря со в результате движения

жидкости, т. е.

вследствие

конвекции;

 

 

 

 

2) второй член

характеризует

изменение поля

В и вихря со

в результате рассеивания их в пространстве, т. е. вследствие процесса диффузии.

Эти утверждения рассмотрены будут более подробно, а сейчас придадим уравнению индукции вектора В (вихря со) еще одну

форму, используя

векторное

тождество

 

rot (vxB)

= vdivB

ß d i v V + ( ß V ) v — (oV) B.

С учетом того, что div В =

0 и div со =

0, получим:

^ L +

( 0 V ) ß = dF = (Bv)v-Bdivv

+ vmAB;

 

 

 

 

(96)

=

(соѴ) v — со — со d іѵ v - j - vÀco .

 

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

движением жидкости

[два первых

члена правой

части выраже

ний (96)] и диффузией

(третий

член).

 

Запись уравнений индукции

(93)

и вихря (96)

в безразмерном

виде облегчает отдельный анализ процессов диффузии и кон­ векции. Представим каждую величину, входящую в уравнение (93),

например вектор

В,

в виде произведения масштаба

В0 на

безраз­

мерное значение

В,

т. е. В = В0В;

ѵ — ѵ0ѵ;

x =

l0x; y

= l0y\

z = /0 z; t = t0t; vm

= v m 0 v m .

 

 

 

 

Подставим эти выражения в уравнение (96) и учтем, что:

Л,.2

 

Л„22

~

Л,'2

дВ0В

+

 

 

~

Lö(V)

 

 

дх

 

ду

 

dz

 

 

 

 

В0

 

д2Ъ~ .

 

 

 

 

 

 

/„

 

дх2

L

 

 

 

109



Из-под знака производной -щ- выйдет масштабный множитель

-^г-,

а из-под знака

операции rot

{v X В) множитель

.

' о

 

 

 

 

' о

Тогда

будем иметь

в

безразмерном

виде следующее уравнение:

 

^ • ^ - = ^ r o t ( w x î ) + ^ # - А Ѣ

 

 

' о

dt

' о

/ 5

 

или

= r o t ( ' 0 x ß ) + ^ - A 5 ,

где Rem =

Очевидно, что для гидродинамического вихря ю безразмерное уравнение будет иметь вид

 

-ПГ-^

= rot (ихсо) -4- p^-Acû.

 

 

' о « о

dt

 

R e o

 

Числовое соотношение между членами, ответственными за

процессы конвекции

и

диффузии,

определяется масштабными

множителями, стоящими

перед ними: единицей и дробью

или

Если процесс конвекции преобладает над процессом диф­

фузии,

то

 

 

 

 

 

1

»

Re7' Т - е -

^ » 1 "

 

Если же процесс диффузии преобладает над процессом конвек­ ции, то -p-î— ^> 1, т. е. Rem <^ 1. В данном случае характерный безразмерный комплекс

Rem = Т~ = ѵ1^та^

называется магнитным числом Рейнольдса (по аналогии с гидро­ динамическим числом Рейнольдса Re = vl/v).

Рассмотрим отдельно процессы конвекции и диффузии. Для

определенности будем рассматривать вектор В.

Конвекция вектора индукции магнитного поля В наиболее четко проявляется в движущейся среде, когда в ней велико ма­ гнитное число Рейнольдса Rem . Это возможно, в частности, при больших значениях электропроводности о, когда явление конвек­ ции резче проявляется в жидких металлах, чем в ионизованных газах. При полном отсутствии процесса диффузии явление кон­ векции приводит к тому, что магнитное поле перемещается вместе с движущейся электропроводной жидкостью, а если движется магнитное поле, то оно приводит в движение жидкость.

ПО


В рассматриваемом случае (когда число Rem велико) говорят об эффекте «вмороженности» магнитного поля в среду, поэтому

размерное уравнение индукции упрощается

и имеет вид

f - = rot(äx5) .

(97)

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

Рассмотрим на примере вектора В магнитной индукции тео­ рему, имеющую место при больших числах Rem .

Рис. 50. Эффект «вмороженности» магнитного поля в проводящую среду

При движении жидкости с бесконечно большим магнитным числом Рейнольдса поток вектора В магнитной индукции через любую «жидкую» поверхность s не изменяется с течением времени.

В частности, это означает следующее: пусть в начальный мо­ мент времени tx взята какая-либо «жидкая» поверхность s, це­ ликом состоящая из магнитных силовых линий вектора индук­ ции В (рис. 50). Поток J Вп ds в момент времени tx равен нулю

s

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

сохранит

постоянную величину,

если число Rem велико; т. е.

в момент

t2 он будет также равен

нулю, а это возможно, если век­

торы индукции В и /г взаимно перпендикулярны в каждой точке поверхности s. Следовательно, эта поверхность в любой момент времени состоит из магнитных силовых линий, которые оказались как бы «вмороженными» в частицы жидкости, составляющие ее.

Если условие Re,„ > сю не соблюдается, то жидкость при дви­ жении проскальзывает через магнитные силовые линии, или же можно сказать, что магнитное поле диффундирует внутрь объема вещества, выделенного для рассмотрения.

Подобная теорема для гидродинамического вихря со скорости известна как вторая теорема Гельмгольца о вихрях.

111