
книги из ГПНТБ / Бальчитис А.А. Емкостная подобласть индукционных процессов преобразования потоков энергии
.pdfобусловленных движением сред с конечной проводимостью в электрическом поле. Процессы оказываются более сложными, чем первоначально предпола гались. Поэтому X . Альфвен делает вывод: „ ... несмотря на всю свою эле гантность, теория первого подхода (т.е. метода магнитной гидродинамики без учета электрогидродинамических явлений . — А . Б . ) почти ничего не может поделать с реальностью ..
МГД-приближение рассматривает лишь индуктивную подобласть процессов электромагнитного („фарадеевского") преобразования потоков энер гии, для которой характерны большие проводимости рабочего тела, малые скорости движения (при больших скоростях движения уменьшается темпе ратура, ионизованность и проводимость рабочего тела), низкие частоты (<о£-»0), большой удельный вес, габариты и значительные потери (наличие магнитопроводов и обмоток).
Преобразовательная техника развивается в направлении освоения более интенсивных процессов, больших скоростей движения рабочего тела, высоких частот и напряжений, преобразователей с хорошими весо-габаритными харак теристиками и малыми потерями. Этим требованиям в значительно большей ме ре, чем индуктивные, отвечают емкостные индукционные методы преобразова ния потоков энергии. Поэтому возникает необходимость рассмотрения пове дения слабопроводящих сред в электрическом поле, т.е. электрогидрогазодинамической подобласти процессов электромагнитной гидродинамики.
Исследования в области электрогидрогазодинамики - науки, предме том которой является описание поведения неэлектропроводящих жидкостей или газов (диэлектрических сред) в электрическом поле, немногочисленны [102-112 и др . 2 ] .
Работы эти не систематизированы и относятся, в основном, к области исследования конвекционных, электрокинетических и электростатических явлений. Меткая общая характеристика состояния исследований в области ЭГД-явлений дана Г. А. Остроумовым: „ . . . эта отрасль науки находится только в начальной стадии своего развития. Это заключение вытекает из того, что в несметном количестве опубликованных работв этой области накопи лось слишком много противоречивого эмпиризма и произвольных частных допущений и рассуждений при отсутствии отчетливых объясняющих идей" [106].
1 X . Альфвен. Физика плазмы, космические исследования и происхождение солнечной системы. У Ф Н , т. 104, вып. 4, 1971, 529-537.
2 Исторический обзор и обширная библиография по исследованиям в области электро гидродинамических явлений даны в работе: W. F . Picard, In Progress in dielectrics. Academic Press, New York, 1965.
100
Из последних работ, посвященных проблемам ЭГД-конвекционных тече ний, следует отметить работу Г. Н . Копылова [107], в которой дается анали тическое решение задачи ЭГД-течения жидкости со слабой униполярной проводимостью в плоской трубе, а т а к ж е выясняется характер изменения параметров течения под действием внешнего электрического поля различной напряженности и направленности.
Попытка обобщения исследований в области конвекционной электро газодинамики предпринята в монографии И. Б . Рубашева и Ю. С. Бортникова [228]. В этой работе анализируется общая система уравнений конвекцион ной электрогазодинамики, рассматриваются частные случаи и приводятся решения некоторых технических устройств. Теоретические данные сравни ваются с экспериментальными. В числе основных направлений технического использования ЭГД-конвекционных процессов указаны: электрогазодинами ческая технология, преобразование энергии, управление свойствами среды, диагностика (определение физических характеристик сред с высокой концен
трацией |
носителей электрических |
зарядов). |
|
||
В конвекционной |
ЭГД |
рассматриваются квазистатические |
электричес |
||
кие поля, и основные |
уравнения представляют собой законы электростатики |
||||
(табл. |
3.1). Существенной |
чертой |
математической модели, |
описывающей |
ЭГД-конвекционные процессы, является равенство нулю вихря вектора на пряженности электрического поля Ё ( у х | = 0 ) .
Конвекционная электрогидродинамика обычно рассматривает движение жидкости с электрическим взаимодействием на границе раздела фаз, т.е.
исследует действие поверхностных |
электрических сил. |
Следовательно, конвекционной подобластью электрогидрогазодинами- |
|
ческих течений может быть названа |
„ . . . т а к а я система, в которой определя |
ющими являются электростатические и газодинамические (или гидродинами ческие — А.Б.) силы" [228, стр. 17].
Следует заметить, что при частичном заполнении рабочего зазора емкост ной машины диэлектриком наблюдается оседание на поверхность поляризо ванного диэлектрика свободных электрических зарядов, уменьшающих напря женность результирующего электрического поля и ухудшающих рабочие
характеристики преобразователя. Подобные явления не могут |
наблюдаться |
в индуктивных машинах (отсутствие свободных магнитных зарядов). |
|
Конвекционные явления в емкостных машинах исследованы эксперимен |
|
тально. Опубликованы некоторые теоретические соображения |
[92 и др.]. |
Эти явления должны учитываться при разработке емкостных машин с частич ным заполнением рабочего зазора диэлектриком.
101
Т а б л и ц а |
3.1 |
|
|
|
|
|
|
Основные дифференциальные уравнения, формулы преобразования и |
граничные условия |
||||||
для квазистатических |
электрических полей, согласно [92] |
|
|
||||
|
Уравнения |
|
Формулы преобразования |
Граничные условия1 |
|||
у х £ = 0 |
|
|
Е' = Е |
|
л х [ £ ] = 0 |
||
VD = qe |
|
|
|
|
й •[!>] = |
а |
|
V 8 + - ^ |
Че = 0 |
|
q'e=qe |
|
|
да'- |
|
|
|
|
|
|
|
|
|
Ъ = |
е0-Ё+Р |
|
8'= 8 — qe |
• v |
= п v [qe] |
- dt |
|
|
|
|
|
Р'=р |
|
|
|
В таблице 3.1 кроме общепринятых использованы следующие обозначения: |
|||||||
а — поверхностная |
плотность |
электрического |
заряда; |
|
|||
К — поверхностная |
плотность |
тока. |
|
|
|||
Квазистатические уравнения конвекционной электрогидродинамики ин |
|||||||
вариантны |
относительно |
преобразований Галилея. |
|
|
|||
В |
отличие от |
конвективной |
индукционная |
электрогидрогазодинамика |
рассматривает поведение слабопроводящих жидкостей и газов не в квази статических, а переменных (динамических) электрических полях, и основными уравнениями являются не уравнения квазиэлектростатики, а уравнения электродинамики (уравнения Максвелла). Отличительной чертой уравнений индукционной электрогидрогазодинамики является неравенство нулю вихря вектора напряженности электрического поля £ (у х £ # 0 ) , использование преобразований Лоренца и исследование действия не поверхностных, а объем ных сил электрического взаимодействия.
Индукционную электрогидрогазодинамику можно рассматривать как раздел классической электродинамики, связанный с исследованием электро механических и электродинамических эффектов в слабопроводящих жидкос тях и газах, вызванных действием сил электрического индукционного взаимо
действия. С другой |
стороны, электрогидрогазодинамику можно классифици |
|||||
ровать как |
часть |
механики сплошных |
сред, |
когда исследование поведения |
||
жидкостей или газов ведется с учетом |
сил |
электрического индукционного |
||||
взаимодействия. |
|
|
|
|
||
1 |
[A] = AW |
— AW |
— скачок величины А |
на |
поверхности. |
|
2 |
V s — поверхностный оператор набла. |
|
|
|
202
Исследование электрогидрогазодинамических индукционных течений важ но не только в связи с перспективой разработки новых эффективных ЭГДИ - и ЭГазДИ-методов преобразования потоков энергии, но и в связи с потребнос тями важных отраслей промышленности, интенсивно осваивающих прогрес сивные процессы электронно-ионной технологии. „Возможное применение электрогазодинамических процессов имеет очень широкие границы — от гене рирования энергии на Земле до снабжения электроэнергией космических ап паратов" [113].
3.2. Основные уравнения индукционной электрогидрогазодинамики
Уравнения индукционной электрогидрогазодинамики могут быть полу чены из более общих уравнений электромагнитной гидрогазодинамики — науки, рассматривающей движение сред с любой заданной проводимостью в электрических и магнитных полях.
Электрогидрогазодинамическое приближение будет хорошо выполняться в случае большого сопротивления среды, т.е. когда токи и соответственно энергия магнитного поля пренебрежимо малы по сравнению с энергией элек трического поля (se0 • E2>\x\i0 • Н2).
Относительная диэлектрическая проницаемость газа близка к единице, поэтому во всех уравнениях злектрогазодинамического приближения и в
некоторых |
случаях электрогидродинамического приближения |
можно |
при |
||||
нять s a |
1. |
|
|
|
|
|
|
Система |
уравнений |
электрогидрогазодинамического |
приближения ма |
||||
тематически |
в значительной степени отличается от известных уравнений |
маг |
|||||
нитной |
гидродинамики |
наличием членов, учитывающих |
токи |
конвекции. |
Когда ж е токи конвекции отсутствуют, эти уравнения математически стано вятся аналогичными уравнениям магнитной гидродинамики и могут быть названы уравнениями электрогидрогазодинамического индукционного при ближения [114].
В уравнениях электрогидрогазодинамического индукционного приближе ния уравнения Максвелла, которые могут быть преобразованы в дифферен циальные уравнения в частных производных (аналогично уравнениям Л а п л а с а , Пуассона или Гельмгольца), рассматриваются совместно с нелинейными дифференциальными уравнениями в частных производных Навье —Стокса, включающих дополнительный член, учитывающий электромагнитные объем ные силы или — совместно с уравнениями газодинамики.
Совместное решение указанных уравнений определяет распределение то ков емкостной индукции в канале и электромагнитные объемные силы, влия -
103
ющие на характер движения элемента рабочего тела, закономерность течения, р я д электромеханических и электрогидродинамических эффектов, а также интегральные величины (напор, мощность, потери напора и мощности, элект рические потери и др.).
Полученные при этом математические соотношения служат основой для расчета и проектирования ЭГДИ-преобразователей с жидкообразным или газообразным рабочим телом, а также при исследовании режимов их работы.
Решение уравнений электрогидрогазодинамического индукционного при ближения в общем виде затруднительно и составляет предмет математической физики. Однако в частных случаях элементарного характера могут быть полу чены простые выражения, удовлетворяющие решению технических задач, связанных с рассмотрением процессов преобразования потоков энергии.
Если течение будет турбулентным, а электромагнитные поля — перемен ными, то решения уравнений электрогидродинамического индукционного приближения аналитически получены быть не могут. В этом случае объемные силы, действующие на элемент рабочего тела, являются переменными и носят нерегулярный характер; решение соответствующих уравнений возможно только при условии введения некоторых упрощающих допущений и использования электродинамического приближения.
При упрощенном рассмотрении нестационарных явлений в уравнении Макс велла (1.28) целесообразно пренебречь реактивной составляющей плотности магнитного тока, полагая
$M)=j<x>[HL0-H=0. |
|
|
|
|
||
Чтобы оценить,когда |
можно пренебречь членом ja>y.[L0 |
• Н, сопоставим |
этот |
|||
член с активной |
составляющей плотности магнитного |
тока |
|
|||
У м - я . |
|
|
|
|
||
Предположим, |
что Н изменяется во времени синусоидально, тогда |
для |
||||
отношения Ь(г) |
к 8( о ) получим |
равенство |
|
|
||
М |
М |
|
•/ |
L |
|
|
[8|,в ) ]« |
|
Гы |
Тм |
|
|
|
или, с учетом (3.9), полагая к^-к^х l 1 , \ixl, - зависимость
1 Произведение к(/ • к\ имеет размерность ж - 2 [см. (3.9)].
104
Величина уе |
для большинства технических задач может быть принята рав |
||
ной Ю - 1 1 ом-1, |
м-1, |
[ а 0 = 1 . 2 6 - Ю - 6 гн-м-1. |
Следовательно, |
№1* - ; 1 |
0 - " . с о . |
|
|
[ 8 L r > l m |
|
|
|
т.е. даже при весьма высоких частотах реактивной составляющей магнитного тока можно пренебречь.
Аналогичное сопоставление величин активной составляющей объемной силы магнитного —/э(°> (1.48) и электрического индукционного взаимодействия —/эл (1-50) позволяет установить условия, при которых можно пренебречь си
лой магнитного индукционного взаимодействия |
/<£) . |
||
Активная составляющая объемной силы магнитного индукционного взаи |
|||
модействия (1.48) может быть |
представлена равенством |
||
\¥еа)хВ\хуе |
(E+vB)B. |
|
|
Если Е — величина |
малая, |
то это равенство |
упрощается |
\Ь{е] xB\xye-v |
ВК |
|
|
Аналогично для активной |
составляющей объемной силы электрического |
индукционного взаимодействия (1.50) можно получить приближенную зависи мость (kbi-kftx 1 м~2):
Отношение активных составляющих объемных сил магнитного и электри ческого индукционного взаимодействия равно
| 8 с а ) х Д : _ Те • vB* _ 2 у.(х0
Подставляя в это отношение ранее принятые значения |
( у е = l O - 1 ^ - 1 - ^ " 1 , |
|||||
( х 0 = 1 , 2 6 - 1 0 - « г и - д е - 1 , |
г0 = 8,85 • Ю - ^ ф - . и " 1 |
, |
е = (г=1), |
|||
получаем |
|
|
|
|
|
|
И ^ х я ! _ П 0 _ 1 Х ч 2 |
1 , 2 6 - ю - " |
, 0 _ 1 7 |
|
|
||
j7jx8<M ">| |
' |
8,85 - 10 - " |
|
|
|
|
т. е. — величину пренебрежимо |
малую . |
|
|
|||
Приведенные соображения показывают, что при рассмотрении технических |
||||||
задач в уравнениях |
электрогидрогазодинамики можно пренебречь реактивной |
|||||
составляющей комплекса |
плотности |
магнитного тока |
и объемной силой ма |
|||
гнитного индукционного |
взаимодействия. |
|
|
105
Таким образом, при рассмотрении технических задач электрогидрогазодинамики можно пользоваться упрощенной моделью электромагнитного поля, описываемого приближенными уравнениями Максвелла
|
у х Я = |
дЪdt |
|
|
(3.1) |
|
|
у х £ = |
- 8 . |
|
|
(3.2) |
|
|
|
|
|
|
|
(3.3) |
где ум — |
удельная |
магнитная |
проводимость (1.29). |
|
||
Из |
уравнения |
(3.3) с учетом (3.2) |
получается: |
|
||
|
|
|
|
|
|
(3.4) |
Подставляя (3.4) в (3.1), после |
преобразования получаем |
|||||
|
дЕ - = V х (v x £ ) ~ V x |
v e(V x E ) , |
(3.5) |
|||
|
dt |
|
|
|
|
|
где |
|
|
|
|
|
(3.6) |
|
Ve = (Ум • £0 |
• s) 1 |
|
|
||
— коэффициент |
электродинамической |
вязкости или |
диффузии. |
Необходимо заметить, что в классической электродинамике закон Ома — зависимость напряженности магнитного поля от напряженно^, ги электричес
кого поля обычно принято выражать равенством |
|
8^ = (у х Н) = ув-£. |
(3.7) |
В данной работе, с целью получения уравнений индукционной |
электро |
гидродинамики, математически аналогичных известным уравнениям |
магнит |
ной гидродинамики, закон Ома записывается в другой форме, а именно [235]:
ам 1 ) = ( у х Ё ) = у м - Я . |
|
(3.8) |
Нетрудно показать, что |
между зависимостями |
(3.7) и (3.8) нет принци |
пиального различия. Действительно, применяя операцию rot и теорему Стокса, уравнения (3.7) и (3.8) могут быть преобразованы следующим образом:
б1«) = ( у х Я )
у х 8 ] а > = [ у х ( у х Я ) ]
Т е ( у х £ ) = [ у х ( у х Я ) ]
Je-S^ |
/ |
dsM |
= (f> ( у х Я ) dL |
|
|
f |
dsM |
|
i ( a ) S_ |
||
o1a ) = y e - S |
ФdlM |
§<a) = ( v x £ )
V x Й а ) = [ y x (v x Ё)]
y„ (V x Я ) = [у x (у x Ё)]
У м - S ^ / dse = j) |
(yxE)dle |
°м — Ум °e
ФdL
106
Следовательно, для перехода от обычно используемой формулы закона Ома (3.7) к (3.8), применяемой в уравнениях электрогидрогазодинамики, необ ходимо произвести замену удельных проводимостей согласно равенству
Ye=Y-1-ky-kM |
, |
(3.9) |
|
где k(f) = ~- |
коэффициент формы (геометрии) электрического поля, |
м~1; |
\dse
S
фdlM
£ ( м ) _ ь_ коэффициент формы магнитного поля, м~1.
\dsM
S
Таким образом, зависимость (3.8) является физическим аналогом закона Ома.
Принимая во внимание векторное тождество
V x ( V x £ ) = v ( v £ ) - V 2 £
и учитывая, что в рассматриваемом |
случае у £ = 0 , |
равенство |
(3.5) |
преобразу |
|||||||||||
ется |
так |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
J |
= V x ( s x £ ) |
+ v e . V 2 |
£ . |
|
|
|
|
(3.10) |
|||||
|
Уравнения (3.5) и (3.10) представляют |
собой |
дифференциальное |
урав |
|||||||||||
нение индукции электрического поля в пространстве и во времени |
[115, |
116]. |
|||||||||||||
|
Используя |
векторное |
уравнение |
|
|
|
|
|
|||||||
|
|
V x ( » x £ ) = (Ёу) |
|
v - |
(z>y) Ё+ v (V-S) - |
Ё (yv), |
|
|
|
||||||
можно |
получить еще одну форму уравнения |
(3.5) |
|
|
|
|
|||||||||
|
|
д? |
- |
- |
- |
- |
- |
|
— |
|
|
|
|
|
|
или |
|
~ |
+ (су) Е-(£у)» |
|
+ Е(v») |
= ve • Д £ |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
§--(Ev)t |
|
+ E(Vv) |
= ve-A£, |
|
|
|
|
(3,10') |
|||||
которое |
сокращенно |
записывается |
так: |
|
|
|
|
|
|||||||
|
|
h e l m £ = v e - A t f , |
|
|
|
|
|
|
|
|
|
(З.Ю') |
|||
где helm — векторный линейный [оператор — гельмгольциан, |
определяемый |
||||||||||||||
для |
любого |
векторного |
поля |
А |
формулой |
|
|
|
|
|
helm А = dA— (A--TJ)V - + A (yv).
Уравнения |
электро гидр о газодинамики |
кроме дифференциального урав |
||||||
нения индукции (3.5) имеют в своем составе следующие известные |
уравнения: |
|||||||
|
•§f |
Pm + V ( p m - ^ ) = 0, |
|
|
|
(3.12) |
||
|
pm[~dt |
w + (oV)o]= -VP |
+ riAv + (rl' +1 |
-/jj v ( V » ) + / , |
(3.13) |
|||
|
P» |
S-=^Pm' |
Рт + У ^ - У ^ О + Ф + б г + а - |
(3.14) |
||||
Здесь / |
— сторонние |
объемные |
электродинамические силы; |
|
||||
е |
— внутренняя энергия |
единицы |
массы |
газа; |
|
|||
X |
— |
теплопроводность; |
|
|
|
|
Ф— диссипативная функция;
QT |
— теплота, получаемая |
единицей объема газа; |
|
|
||||||
Qd |
— член, характеризующий диффузию |
газа. |
|
|
||||||
Суммарные объемные силы индукционного взаимодействия |
определяются |
|||||||||
зависимостью |
|
|
|
|
|
|
|
|
|
|
/ = / э л + Л м == (V * £ ) х £ + (у х Я ) х В = |
|
|
|
|||||||
= |
[(его)"1 (у х D) х D + (щ^оГ1 (у х В) х В]. |
|
|
(3.15) |
||||||
Учитывая, кроме того, векторные тождества |
|
|
|
|||||||
(у |
xD)xD |
= {D\j)D-D(yD) |
|
= {D^y)D-~ |
yDz, |
|
|
|||
(ухВ)хВ |
|
= (Ву)В-В(уВ) |
|
= (Ву)В-~ |
уВ2, |
|
|
|||
равенство |
(3.15) |
преобразуется |
так: |
|
|
|
|
|||
VP = ( ^ о Г 1 |
D + (и-^о)"1 CBV) В - |
|
|
|
||||||
- |
j l ^ y ^ ^ |
+ M - ^ |
s 1 |
] . |
|
|
|
(3.16) |
||
Следовательно, |
|
|
|
|
|
|
|
|||
V |
[р + у |
(еео)-1 |
X»2 + \ |
( № о |
) - 1 Я 2 ] = (е s0 )" * ( 5 у ) Я + Ы |
" 1 (Б V) 5 |
||||
и уравнение движения |
(3.13) может быть представлено равенством |
|
||||||||
Р» [ J |
+ ( 5 V) о] = - V [Р + ( ^ о ) - 1 |
D* + (jx^Lo)-1 2?2 ] + |
|
|
||||||
+ riAv |
+ (rl' + j |
т)) v ( V ^ ) + ( ^ 0 |
) - 1 ( ^ V ) ^ + ( № o ) - 1 ( S v ) 5 . |
(3.17) |
108
|
Теплота, получаемая единицей объема |
газа из электромагнитного |
поля |
||||||||||
в электрогазодинамическом |
приближении, |
определяется |
зависимостью |
||||||||||
|
QT |
= $i?-H=(Vx£)a |
|
[ e e 0 . v e ( v x £ ) - e e 0 ( » x £ ) ] , |
|
(3.18) |
|||||||
где |
= - |
(V х |
Е)а — активная |
составляющая вектора плотности магнит |
|||||||||
|
|
|
|
ного тока; |
|
|
|
|
|
|
|||
H=(vxD) |
— *{u' (V х-Ё)— |
из |
совместного рассмотрения (3.2) и (3.3). |
||||||||||
|
С учетом |
(3.18) |
уравнение |
энергии |
(3.14) |
записывается так (членами Ф |
|||||||
и Qd |
пренебрегаем): |
|
|
|
|
|
|
|
|
|
|
||
|
Рт ^=P?~l |
| - p , „ + V ( ^ V r ) + ( v x £ ) [ss0 • ve (у х Ё) |
- |
|
|||||||||
|
-ze0(vxE)]. |
|
|
|
|
|
|
|
|
(3.19) |
|||
|
Используя понятие энтальпии торможения, уравнение энергии в электро- |
||||||||||||
газодинамическом приближении можно представить равенством |
|
||||||||||||
|
|
dh |
dp |
,, |
_,. |
|
|
|
|
|
|
|
|
|
+ |
(у х Ё) [es0 • ve (у х Е) - |
ге0 (о х £ )], |
|
|
(3.20) |
|||||||
где |
h = Ср |
- Т+ |
^ V2— энтальпия |
торможения; |
|
|
|
|
|||||
|
Ср — теплоемкость при |
постоянном |
давлении. |
|
|
|
|||||||
|
Следовательно, |
основными |
|
уравнениями |
индукционной |
электрогидро- |
|||||||
газодинамики |
являются зависимости (3.5), (3.11), (3.12), |
(3.17) |
и (3.19) [114]. |
||||||||||
|
Нетрудно заметить математическую аналогию рассматриваемых |
урав |
|||||||||||
нений и уравнений |
магнитной |
гидродинамики. Это обстоятельство в |
значи |
тельной степени облегчает решение задач электрогидрогазодинамических индукционных течений.
Если в задачах электрогидрогазодинамических индукционных течений возникает необходимость рассматривать поведение заряженных сред в попереч ном электрическом поле, то в уравнениях индукционной электрогидрогазодина мики будут появляться дополнительные члены, обусловленные конвекцион ными явлениями. Теория электрогидрогазодинамических конвекционных те чений в настоящее время достаточно полно разработана [см. 92 и др.].
0.0. Безразмерная форма уравнений индукционной электрогидрогазодинамики
Чтобы определить наличие того или иного специфического электрогазо динамического явления, а также сравнить течения различного масштаба, целе сообразно провести оценку относительных величин членов, входящих в урав-
109