ПОЛЯРИЗАЦИЯ (POLARIZATION) |
Главная | История исследования | Роль поляризации плазмы |
Энергобаланс аргоновой дуги | Каналовая модель |
КАНАЛОВАЯ МОДЕЛЬ ЭЛЕКТРИЧЕСКОЙ ДУГИВведениеНа предыдущей странице "Роль поляризации..." показано, как анализ экспериментальных данных, привел к выводу, что неизлучательные потери энергии из канала электрического сильноточного квазистационарного разряда определяются в основном поляризационными процессами в плазме канала. Этот результат явился новым и совершенно неожиданным, впервые он был неявным образом использован в работе [1, 2] для анализа энергобаланса электрических квазистационарных разрядов как с преобладанием излучательных потерь энергии из разряда, так и без таковых. Здесь упомянутый анализ выполнен более подробно, чем в работе [2], для разрядов без преобладания излучательных потерь энергии. Типичными представителями таких разрядов являются электрические дуги, горящие при давлениях порядка атмосферного или выше. Учет поляризационных процессов позволяет по-новому решить задачи, касающиеся как зависимости температуры канала дуги от внешних параметров, так и расчета вольтамперных характеристик дуг (как свободно горящих, так и стабилизированных стенками). Чтобы наглядно увидеть разницу в решении названных задач с помощью нового подхода, учитывающего поляризацию плазмы канала, и традиционного, корни которого уходят к тридцатым годам прошлого века, напомним сначала суть последнего. 1. Классическая каналовая модель.1.1 Уравнение Эленбааса-Хеллера (1934 г.)Уже в тридцатых годах прошлого столетия на основе теоретических и экспериментальных работ сформировалось представление, что в столбе дуг высокого давления с достаточно хорошим приближением реализуется состояние термической плазмы. Под этим термином подразумевается, что все основные процессы в такой плазме являются однозначными функциями температуры, одинаковой для всех плазменных частиц. На основе этого Эленбаасом и Хеллером была предложена теория столба дуги, которая послужила основой для последующих теорий. Проведем элементарный вывод основного уравнения этой теории, описывающего баланс энергии в плазме разряда. Предварительно оговорим исходные условия. Анализируем однородный по длине цилиндрический столб дуги в продольном электрическом поле, напряженность которого Е. Считаем, что дуга горит в неподвижном газе, находящемся при постоянном давлении в охлаждаемой трубке радиуса Rt. Рассматриваем режимы, при которых излучательные потери энергии из столба малы, а потери энергии из столба определяются в основном процессами теплопроводности в столбе дуги. Как известно, тепловой поток в среде
описывается выражением Здесь J - удельный тепловой поток, λ - коэффициент теплопроводности, Т - температура. В силу допущенной однородности столба по длине и его цилиндрической симметрии все параметры вещества в столбе являются функциями только текущего радиуса столба r. В частности, градиент температуры в этом случае представляет собой производную от температуры по радиусу, то есть в столбе реализуются лишь радиальные потоки тепла. Чтобы записать дифференциальное
уравнение, описывающее связь между параметрами столба, проанализируем
энергобаланс тонкого цилиндрического слоя столба радиуса r и
толщины Δr. Для упрощения записи выкладок рассмотрим энергобаланс
цилиндрического слоя в расчете на единицу его длины. Очевидно, что
энергобаланс слоя складывается из следующих частей: входящего теплового
потока Q+ через внутреннюю стенку слоя радиуса r;
тепловой энергии QE, генерируемой в слое протеканием
тока, и выходящего теплового потока Q- через наружную
стенку слоя радиуса r+Δr. В итоге баланс запишется в виде: Раскроем выражения, входящие в энергобаланс. Входящий тепловой поток представим в виде: Мощность тепловой энергии QE , генерируемой в цилиндрическом слое за счет протекания тока, запишем как: где σ(r) - удельная электропроводность среды в слое радиуса r. Величину Q- выразим следующим образом: Подставляя
выражения (3-5) в уравнение энергобаланса (2), прийдем после
элементарны преобразований полученного выражения к дифференциальному
уравнению, которое, обычно, записывается в виде: Это и есть известное уравнение Элеенбаса-Хеллера. Граничные условия для него: на оси столба (r=0) производная dT/dr=0 вследствие симметрии, при r = Rt температура T = Tt , где Tt - температура стенок трубки. 1.2 Каналовая модель дуги (1932 г.). Принцип минимума Штеенбека. Решить уравнение (6) в общем виде аналитически не удалось до настоящего времени. Однако разрабатывались различные приближенных методы его решения [3]. При постоянной силе тока дуги I и заданной температуре стенок трубки Tt это выражение связывает между собой три переменные величины T, R и E. Чтобы получить еще одно независимое выражение, связывающее эти величины, проинтегрируем уравнение (6) в бестоковой зоне каналовой дуги, учитывая, что в этой зоне σ=0. Интегрирование показывает, что величина
rJ(r) постоянна в этой зоне, что само по себе очевидно, кроме того,
она равна: Это
следует из того, что электрическая мощность, выделяемая в канале, равна
IE, поэтому тепловой поток на границе R равен
IT/2π. Эта
функция обладает тем свойством, что производная от нее, как легко видеть,
равна тепловому потоку J, взятому с обратным знаком: Из сравнения выражений (8) и (10) вытекает легко интегрируемое дифференциальное уравнение: Интегрирование обеих частей этого выражения от Rt до R приводит к следующему уравнению, связывающему величины I, T и R: Таким образом, мы имеем два уравнения (7) и (12), которые при заданной силе тока I и заданных условиях окружающей среды (T=Tt при r=Rt) связывают три переменные T, r и E. Исключая из этих выражений одну из переменных, например, переменную T, получим бесконечную последовательность режимов горения дуги выражаемую функциональной зависимостью E (r). Чтобы из
указанной последовательности режимов можно было выбрать тот, который
фактически реализуется, используется допущение, известное как минимальный
принцип Штеенбека. Согласно этому допущению при заданной силе тока и
фиксированных условиях окружающей среды реализуется тот режим, при котором
напряженность электрического поля минимальна. Другими словами при этом
режиме производная Принцип Штеенбека был не раз проверен при исследовании дуговых разрядов. При этом были получены результаты, удивительно хорошо согласующиеся с экспериментом [4]. Однако вопрос обоснования этого принципа не перестает волновать исследователей до настоящего времени. Подводя итоги,
можно сказать , что каналовая модель позволяет рассчитать температуру и
напряженность поля в дуге. Однако эти расчеты сложны, поэтому
имеются расхождения между данными разных авторов. Ведь при расчете
теплопроводности λ(Т) плазмы приходится учитывать большое
количество процессов, поскольку она складывается из величин
теплопроводности молекул (в случае молекулярных газов), атомов, ионов,
электронов; кроме того, нужно учесть также вклад от диффузии энергий
ионизации и диссоциации.
2. Модернизированная каналовая модель, учитывающая поляризацию плазмы.2.1 Энергобаланс однородного гипотетического канала дуги с учетом поляризации. Плазма представляет собой ионизованный квазинейтральный газ.
Квазинейтральность означает, что число положительных и отрицательных зарядов в плазме
одинаково. Дальнодействующее кулоновское взаимодействие заряженных частиц приводит
к целому ряду качественных особенностей плазмы, отличающему ее от обычного газа,
в котором нет заряженных частиц. Это дает основание считать плазму особым, "четвертым"
состоянием вещества. Таким образом, вещество, согласно сказанному, может находиться в
твердой, жидкой, газообразной и плазменной фазе. Причем переход вещества из одной фазы в другую происходит с ростом температуры. Это уравнение энергобаланса, записанное для зоны проводимости каналовой модели с учетом поляризации плазмы, является третьим уравнением, связывающим между собой величины T, E и R каналовой модели. Подробно это уравнение раскроем и проанализируем чуть позже, а сначала обсудим общие вопросы, связанные с этим уравнением. Установление третьего уравнения позволяет отказаться от принципа
минимума Штеенбека, который, вообще говоря, никак не следует из классической каналовой модели. В свое время это вызвало многочисленные дискуссии [4], которые здесь рассматривать не будем, а обсудим принцип с новых позиций, учитывающих поляризацию плазмы канала.
2.2 Аналитическое выражение для температуры однородного канала.
Найдем аналитическое выражение для модельной температуры канала дуги через ее внешние параметры. С этой целью раскроем уравнение энергобаланса для центральной зоны модели каналовой дуги (14). Чтобы избежать работы с громоздкими выражениями, выкладки проведем для каналовых дуг, горящих в атмосфере газов, состоящих только из одного химического элемента (аргон, азот, кислород и т.д.). Рассмотрение ограничиваем режимами, при которых плазма канала состоит из атомов, однократно заряженных ионов и электронов с концентрациями, равными, соответственно, no, n1 и ne, причем в силу квазинейтральности плазмы ne=n1. Напомним, что, по определению, температура T в центральной зоне модели постоянна. Выкладки проводим для канала единичной длины. Величина излучаемой каналом дуг мощности W для широкого диапазона режимов достаточно незначительна, поэтому при анализе энергобаланса дуг не слишком высокого давления ею пренебрегают. Мы этого делать не будем, но, так как аналитическое выражения для этих потерь не известно, то учтем эти потери формально. Учет потерь выполним с помощью введения множителя δ согласно соотношения: Перепишем уравнение энергобаланса (14) с учетом формул (15-16) в виде, удобном для проведения дальнейших преобразований: Выражение для мощности поляризационных потерь (МПП), которое на этой странице обозначено буквой Õ вместо использовавшейся на других страницах буквы Q, сразу запишем в окончательном виде, вытекающем из общей формулы (2) на странице "История исследования..." при раскрытии ее с помощью выражений (9-15) на той же странице: Символы K1 и α в этой формуле обозначают комбинации физических параметров, причем численные значения всех комбинаций на данной странице сайта даны в системе МКСА: Записанное выражение учитывает то обстоятельство, что порции энергии, затрачиваемые на удаление за пределы канала атомов, ионов и электронов, различаются между собой множителем порядка двойки. В более раннем варианте работы [2] этим различием пренебрегалось. Выразим удельную электропроводность σ плазмы, входящую в энергобаланс (15), через другие параметры плазмы канала. Согласно общеизвестной формуле ее связь с другими параметрами плазмы описывается выражением: Здесь e - заряд электрона, λe - длина свободного пробега электронов, Qeo и Qe1 -эффективные сечения столкновений электронов с атомами и ионами [4]. Через величину Qef обозначена комбинация величин:
При записи (23) учтено, что в слабо ионизованной плазме величина n1no-1 практически равна степени ионизации плазмы η. Запишем рабочее выражение для σ. С этой целью выразим концентрацию электронов в плазменном канале через концентрацию атомов и температуру плазмы с помощью уравнения Саха [6]: Здесь Z1 и Zo - статистические суммы ионов и атомов, h - постоянная Планка, Eo - потенциал ионизации атомов. Для записи зависимости температуры канала в удобном виде используем вместо величины концентрации атомов в канале no величину давления p, при котором горит дуга. Давление связано с концентрацией законом Дальтона; с учетом малой степени ионизации плазмы в дугах, а также полной диссоциации молекул в канале в случае молекулярного плазмообразующего газа можно приближенно записать: Подставив выражения (18) и (25) в (17), учтя (26) и извлекши квадратный корень из обеих частей записанного выражения, сгруппируем множители, в него входящие, следующим образом. В одной части до знака равенства оставим множители, в которые явным образом входит температура T, а в другой части разместим все остальные множители. Для удобства работы полученное выражение разобъем на два отдельных. Для этого каждую его часть приравняем одному и тому же выражению B(T). В результате получим искомую зависимость, связывающую T c I, R, p и параметрами вещества, в следующем виде: Через F в (29) обозначена комбинация: Для полноты картины запишем в явном виде выражение для величины Fo, используя зависимости (20) и (23): Обсудим приемы работы с полученной зависимостью (27-31). Из рассмотрения выражения (28) следует, что величина B(T) зависит только от температуры, но реально это не совсем точно, поскольку энергия ионизации атомов в плазме равна не Eo, а величине Eo-χo, которая уже зависит и от давления плазмы p. Однако, вклад этой поправки незначителен и ею можно пренебречь.
Для дуги, горящей при заданном давлении среды p, левая часть выражения (32) рассчитывается аналитически. Этот расчет удается полностью выполнить для значений давления p, при которых излучательные потери из канала дуги малы (δ<< 1). Затем строится вспомогательный график зависимости температуры T от логарифма lg(BH-1p1/4) левой части выражения (32).
3. Сопоставление модельной зависимости для температуры с экспериментальными данными по аргоновым дугам.3.1. Сравнение с экспериментальной зависимостью. Чтобы записать формулу (27) для конкретного плазмообразующего газа нужно знать сечения столкновения электронов с атомами и ионами этого газа, необходимые для расчета зависимости H(T). Литературные сведения о величинах этих сечений, как будет проиллюстрировано ниже на примере аргона, характеризуются большим разбросом. Поэтому решено было оценить величину H(T), исходя из известных экспериментальных зависимостей. В качестве одной такой зависимости была проанализирована зависимость для средней электронной температуры в аргоновых каналовых дугах атмосферного давления, из [7]. Ниже приведено отсканированное изображение этой зависимости. Сравнение модельной зависимости (27) с эмпирической (33) показало, что, во-первых, величину H, в первом приближении, можно считать постоянной. Во-вторых, из анализа графика зависимости T=lgB(T), изображенного выше, следует, что для обеспечения погрешности расчета температуры около 3% необходимо, чтобы величина Qef, была известна с точностью 26% в области температур около 12000K и 47% в области температур около 7000K. В результате компромиса было выбрано постоянное значение H=1,62· 10-4. Логарифм этого значения lgH=-3,79 [2]. При таком выборе значения H результаты определения величины температуры T по модельной зависимости (27) практически совпадают с расчетными значениями, даваемыми эмпирической зависимостью (33), в области температур T от 7000 до 11500K, причем модельная зависимость оказывается работоспособной в значительно более широком диапазоне режимов [1], чем (33). В частности, модельная зависимость отражает зависимость температуры от давления, чего не делает зависимость (33). 3.2. Оценка эффективных значений величин сечений Qeo и Qe1 на основе экспериментального значения величины lgH. Чтобы рассчитать модельное значение величины H по формуле (29), необходимо знать зависимости для сечений столкновения электронов с ионами и атомами Qe1 и Qeo от температуры. Однако эти зависимости точно не известны. Теоретические расчеты сечений, выполненные разными авторами, существенно различаются. Известные расчеты сечений приведены на сводном графике из [3]. Поэтому представляет интерес оценка величины сечений, исходя из экспериментального значения lgH=-3,79.
Согласно этой формуле экспериментальному значению lgH=-3,79 соответствует величина Fo= 8,55· 10-18. Исходя из выражения (30) для Fo=Qefα-1, оценим величину сечений при различных значениях температуры. Оценки сечений представлены в таблице. В первых столбцах ее приведены термодинамические значения величин no, ne, τ, η и α в плазме канала при различных значениях температуры , а в двух последних столбцах приведены оценки величин Qeo и Qe1. Оценки проведены в предположении, что при T=7000K величина Qeо=αFo из-за малой степени ионизации плазмы, а при температурах T=8000-12000K величина Qe1≈η-1αFo.
Результаты таблицы показывают, что оценки дают в несколько раз большие значения величины сечений Qe1 по сравнению с известными расчетами. Как это объяснить, пока не ясно. На полученной зависимости этого сечения от температуры четко просматривается влияние эффекта Рамзауэра. 3.3. Вольтамперная характеристика каналовой аргоновой дуги.Как уже упоминалось, что, если канал дуги ограничен охлаждаемыми стенками, а величина тока велика, то радиус плазменного канала практичеки совпадает с радиусом ограничивающего канала. В этом случае формула (27) позволяет рассчитать величину средней температуры канала T. Если известна зависимость удельной электропроводности плазмы от температуры, то для упомянутых режимов можно легко рассчитать величину напряженности поля в дуге, то есть вольтамперную характеристику канала дуги. Однако выполнение этого простого расчета упирается в отсутствие надежных сведений о электропроводности плазмы, так как напряженность электрического поля очень чувствительна к этой величине. На следующем рисунке представлены известные экспериментальные данные и расчетные зависимости электропроводности аргоновой плазмы, выполненные разными авторами.
|