Численное моделирование потока газа в дугогасительном устройстве высоковольтного выключателя тема диссертации и автореферата по ВАК РФ 01.02.05, кандидат физико-математических наук Аверьянова, Светлана Андреевна

• Аверьянова, Светлана Андреевна
• кандидат физико-математических наук
• 2005

Оглавление диссертации кандидат физико-математических наук Аверьянова, Светлана Андреевна

Введение.

Глава 1. Математическая модель.

1.1. Общее описание дугогасительного процесса в высоковольтном выключателе.

1.2. Модель для оценки величины предельно допустимого восстанавливающегося напряжения.

1.2.1. Исходная концепция для оценки величины предельно допустимого восстанавливающегося напряжения.

1.2.2. Модель для расчета напряженности пробоя.

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

1.3.1. Дуга в сопле Лаваля.

1.3.2. Схема течения в дугогасительном устройстве.

1.3.3. Схематизация изменения тока во времени.

1.4. Основные уравнения. Материальные функции.

1.5. Модель турбулентности для стационарной дуги.

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

1.5.2. Исходные положения для формулировки модели.

1.5.3. Окончательная формулировка модели.

1.6. Моделирование турбулентности в случае быстро убывающего тока.ЗЗ

1.6.1. Источник роста кинетической энергии турбулентности в единице объема в дуге для случая быстро уменьшающегося тока.

1.6.2. Модель турбулентности для случая быстро уменьшающегося тока.

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

1.8. Модель излучения.

Глава 2. Используемые численные методы и тестирование расчетных ко

2.1. Переход к криволинейной системе координат.

2.2. Численный метод для решения полной газодинамической системы уравнений.

2.2.1. Факторизованная схема Бима и Уорминга.

2.2.2. Вычисление конвективных членов.

2.3. Численный метод интегрирования параболизованной системы уравнений.

2.3.1. Общие соображения к выбору метода.

2.3.2. Метод Поспелова.

2.3.3. Метод обратной задачи.

2.4. Расчетные сетки.

2.5. Тестирование программного кода (Навье-Стокс).

2.5.1. Невязкое течение в сопле Лаваля.

2.5.2. Вязкое течение в круглой трубе.

2.6. Тестирование программного кода для параболизованной задачи.

2.6.1. Невязкое течение в сопле Лаваля.

2.6.2. Вязкое течение в круглой трубе.

Глава 3. Электрическая дуга в сопле Лаваля.

3.1. Обзор экспериментальных и расчетных исследований дуги в сопле Лаваля.

3.2. Постановка задачи для дуги, горящей в сопле Лаваля.

3.3. Методические расчеты.

3.3.1. Исследование сеточной сходимости.

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

3.4. Результаты расчетов и сравнение с опытом.

3.4.1. Влияние учета излучения.

3.4.2. Определение эмпирических коэффициентов.

3.4.3. Стационарная дуга в спутном сверхзвуковом потоке.

3.5. Угасающая дуга в сопле Лаваля.

3.6. Расчет стационарной дуги в сопле Лаваля, основанный на параболизо-ванной постановке.

3.6.1. Постановка задачи.

3.6.2. Методические расчеты: исследование сеточной сходимости.

3.6.3. Выбор эмпирической константы в модели учета турбулентных пульсаций.

3.7.4. Сопоставление расчетов в полной и параболической постановках.

Глава 4. Дуга в радиально натекающем потоке газа.

4.1. Постановка задачи о расчете восстановления электрической прочности.

4.1.1. Описание экспериментальной установки.

4.1.2. Математическая постановка задачи.

4.2 Методические расчеты.

4.2.1. Тест на чувствительность к шагу по времени.

4.2.2. К вопросу о выборе констант в модели турбулентности.

4.3 Влияние силы тока в дуге на течение в канале.

4.4 Динамическая картина размыкания дуги (результаты расчетов).

4.5 Результаты расчетов процесса восстановления электрической прочности.

4.6 Расчеты в условиях реальной геометрии дугогасительных устройств.

Введение диссертации (часть автореферата) на тему «Численное моделирование потока газа в дугогасительном устройстве высоковольтного выключателя»

Газовые высоковольтные выключатели используют как основное силовое коммутационное оборудование в электроэнергетических системах высокого напряжения. Разработка новых и модернизация существующих образцов высоковольтных выключателей, а так же их экспериментальные исследования сопряжены с большими временными и материальными затратами. Поэтому создание математических моделей, позволяющих описывать поведение сильно неизотермического потока газа в проточной части нового высоковольтного выключателя и взаимодействие его с электрической дугой, чрезвычайно важно. Особого внимания требуют процессы, связанные с гашением дуги спутным потоком газа (обычно элегаза — шестифтористой серы, SF6).

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

Суть процесса отключения тока в высоковольтном выключателе состоит в следующем. В момент отключения контакты, находящиеся в дугогаси-тельном устройстве, расходятся, и между ними зажигается мощная электрическая дуга. Одновременно с расхождением контактов начинает двигаться поршень, создающий в дугогасительном устройстве высокоскоростной спут-ный поток газа, который в фазе сильного тока защищает стенки дугогаси-тельного устройства от теплового разрушения. Поскольку отключается переменный ток (с частотой 50 Гц), он быстро уменьшается, а с ним ослабевает и дуга. В момент прохождения тока через ноль начинают конкурировать два процесса: нагрев газа за счет джоулева тепла и отвод тепла от дуги спутным потоком гасящего газа. Если преобладает охлаждение дуги, то газ на некотором участке между контактами становится практически неэлектропроводным. Происходит восстановление электрической прочности газа.

Опыт показывает, что дугогасительные свойства высоковольтного выключателя сильно зависят от формы проточной части дугогасительного устройства, в которой происходит гашение. В канале могут возникать области обратных потоков, в которых задерживается горячий газ, нагретый в фазе сильного тока; динамика газа иногда определяет возникновение зон с пониженной плотностью газа (например, [1]) и другие особенности, которые понижают отключающие способности высоковольтного выключателя.

В элегазовых выключателях (в которых рабочим газом является SF6) восстановление электропрочности происходит в течение очень малого интервала времени (10 — 20мксек) в окрестности момента нуля тока. В течение этого промежутка времени в плазме дуги происходят сложные и пока не полностью изученные физические неравновесные микропроцессы, такие как рекомбинация ионов и электронов в атомы, ассоциации атомов в молекулы, образование так называемых электроотрицательных ионов, ведущее к резкому падению проводимости газа, и другие. В последние несколько лет появились работы, в которых рассматриваются неравновесные микропроцессы [24]. Однако получение численного решения в этом случае даже для простейшей двумерной геометрии проточной части сложно и связано с большими затратами времени и машинных ресурсов. Известно, однако, что, чем больше давление в газе, тем меньше характерные времена этих процессов. В элегазовых выключателях обычно поддерживают давление 6 — 8атм, и характерные времена микропроцессов оказываются малыми - микросекунды и меньше [5]. Эти времена существенно меньше (на порядок и более) характерных времен газодинамических процессов взаимодействия дуги и спутного потока, поэтому в первом приближении для построения модели нестационарного течения в дугогасительном устройстве можно принять гипотезу о локальном термодинамическом равновесии (JTTP) плазмы. Для расчетов ее свойств (material functions, материальных функций) можно использовать соответствующие таблицы, как это делается, например в [6-21].

Большую роль в гашении электрической дуги играет формирование гасящего потока и характер его взаимодействия с дугой, например: распределение скорости и давления в спутном потоке вдоль оси дуги; угол, под которым поток натекает на дугу; турбулизация плазмы дуги при взаимодействии со спутным потоком и так далее. Турбулентность плазмы играет особенно важную роль. Сегодня общепризнано, что отвод тепла от плазмы дуги в период падения тока и прохождения его через ноль происходит за счет именно турбулентной, а не молекулярной теплопроводности [16, 21-25].

Тема данной работы возникла благодаря сотрудничеству кафедры «Гидроаэродинамика» СПбГПУ с «Научно-исследовательским институтом высоковольтного аппаратостроения» («НИИВА»). Основная цель работы заключалась в создании математической модели, позволяющей рассчитывать нестационарные поля турбулентного течения в дугогасительном устройстве высоковольтного выключателя при наличии на ее оси электрической дуги переменного тока. Второй целью являлась разработка полуэмпирической модели восстановления электрической прочности газа между разошедшимися контактами, использующей результаты расчетов нестационарных полей течения.

Базовыми составляющими развиваемой на основе ЛТР-приближения модели являются система осредненных по Рейнольдсу уравнений Навье -Стокса, уравнения неразрывности, энергии и закона Ома, а так же таблицы равновесных свойств плазмы рабочего газа. Дополнительные, но весьма важные составляющие - это предлагаемые в работе модели для расчета турбулентной вязкости и потерь тепла, обусловленных излучением.

Как правило, в предыдущих работах рассматривалась дуга, горящая в спутном потоке, текущем параллельно ее оси, и для описания турбулентного взаимодействия дуги и спутного потока использовались первая или вторая модели Прандтля для расчета коэффициента турбулентной вязкости в струйных течениях. Попытки использования k-с модели при расчетах электрических дуг показали [26], что для согласования расчетных и экспериментальных данных необходимо подбирать зависимость эмпирических коэффициентов модели от величины отключаемого тока. В реальных конструкциях высоковольтных выключателей гасящий поток набегает на дугу (осесимметрично) под большим углом к ее оси, а потом разделяется на два потока, текущие вдоль ее оси в разные стороны. В этих условиях применение моделей Прандтля оказывается неправомерным. Таким образом, необходимы дальнейшие разработки моделей турбулентности, которые позволяли бы рассчитывать взаимодействие осесимметричного потока с нестационарной дугой при натекании на нее потока под произвольным углом. При выполнении данного исследования было принято целесообразным сосредоточится на разработке новой алгебраической модели.

При больших значениях силы тока на характеристики электрической дуги и потока, взаимодействующего с ней, большое влияние оказывает лучистый теплообмен [6, 27-29]. Для определения потерь тепла на излучение существуют разные способы [30, 31]. Однако, достаточно аккуратные методы, базирующиеся на учете спектрального состава излучения, очень сложны для реализации и их применение представляет собой отдельную задачу. В работах, посвященных численному моделированию взаимодействия потока с электрической дугой, чаще всего используются приближенные методы [18, 21]. В настоящей работе для расчета потерь тепла на излучение принята относительно простая полу эмпирическая модель, требующая небольшого числа экспериментальных данных для определения входящего в нее эмпирического коэффициента.

В первой главе диссертации дается описание работы высоковольтного выключателя. Представлены схемы потоков, взаимодействующих с электрическими дугами, описаны основные упрощения, принятые в работе и представлены основные уравнения (полная и парабол изо ванная системы уравнений Навье-Стокса) и модели, используемые в работе. Модель для определения предельно допустимого восстанавливающегося напряжения с гипотезой о поведении напряженности пробоя. Алгебраическая модель для коэффициента турбулентной вязкости, основанная на предположении о двух основных механизмах возникновения турбулентности в потоке при наличии электрической дуги на оси канала: сдвиговой неустойчивости и неустойчивости Рэлея-Тейлора. И модель излучения, основанная на приближении лучистой теплопроводности.

В главе 2 представлены два метода численного решения используемых уравнений. Метод Бима и Уорминга для полной газодинамической системы уравнений, с определением конвективных слагаемых по схеме Роу и комбинированный метод для решения параболизованной системы уравнений. Комбинированный метод представляет собой комбинацию метода обратной задачи в области горловины и метода Симуни (с модификацией Поспелова для М>1) для всей остальной области. В обоих случаях проведено тестирование программного кода: представлены результаты расчетов невязкого течения в сопле Лаваля и вязкого течения в круглой трубе.

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

Глава 4 посвящена процессам, происходящим в стационарной и угасающей дугах в каналах, чья форма приближена к форме реального образца дугогасительного устройства высоковольтного выключателя. Особенностью потока в данном случае является то, что его направление существенно меняется по отношению к оси дуги. Описаны экспериментальные данные, которые легли в основу проводимых в данной главе расчетов, подробно описана постановка задачи о дуге, горящей в спутном потоке, радиально набегающем на дугу. На основе представленных экспериментальных данных определены константы в моделях турбулентности (для члена, описывающего неустойчивость Рэлея—Тейлора) и модели о величине предельно допустимого восстанавливающегося напряжения. Представлены результаты сравнения расчетов с опытными данными для восстановления электрической прочности промежутка после нуля тока. Расчеты, сделанные на основании гипотезы о величине предельно допустимого восстанавливающегося напряжения, лежат в пределах разброса опытных данных. В конце главы 4 представлены расчеты стационарного течения в дугогасительном устройстве реальной формы при наличии на оси канала мощной электрической дуги.

В заключении отмечены основные выводы по работе.

Заключение диссертации по теме «Механика жидкости, газа и плазмы», Аверьянова, Светлана Андреевна

Заключение.

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

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

2. Разработана новая алгебраическая модель турбулентности для расчета турбулентной вязкости в дуге, взаимодействующей со спутным потоком в дугогасительных устройствах высоковольтного выключателя при произвольном угле натекания гасящего потока на дугу. Модель учитывает, что турбу-лизация происходит не только за счет сдвиговой неустойчивости, но и за счет неустойчивости Рэлея-Тейлора.

3. Показано, что при быстром убывании тока причиной интенсивного роста 0' кинетической энергии турбулентности в единице объема в дуге является вклад члена уравнения баланса энергии турбулентности, пропорционального полной производной скорости по времени, обычно не учитываемый в известных моделях турбулентности. Это дало основания к использованию гипотезы «замороженности» кинематического коэффициента турбулентной вязкости при расчете угасающей дуги и последующей начальной фазе восстановления электрической прочности газа.

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

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

6. Проведены расчеты взаимодействия стационарных и угасающих дуг в сопле Лаваля. Получено удовлетворительное согласие с экспериментальными данными по изменениям параметров потока на оси дуги. Установлено, что при быстром убывании тока непосредственно за горловиной сопла образуется участок дуги, на котором она остывает быстрее всего, что согласуется с визуализацией угасающих дуг в опытах. Получено согласие расчетов течения в стационарной дуге в полной и параболизованной постановках. Для реализации последней в работе была создана специальная программа для расчетов стационарных дуг, описываемых параболизованной системой уравнений.

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

8. Проведены расчеты взаимодействия стационарной электрической дуги со спутным потоком в дугогасительном устройстве реальной конструкции. Полученные поля согласуются с общепринятыми представлениями о поведении дуги в канале с двухсторонним дутьем.

Список литературы диссертационного исследования кандидат физико-математических наук Аверьянова, Светлана Андреевна, 2005 год

1. Экспериментальный анализ диэлектрических характеристик неоднородного высокотемпературного элегаза в газовом выключателе // Пер. ст. Toshiyuki U., Koichi I., Hiromichi К., Tetsuya N., Katsumi S. из журн.: Т. IEE Japan. 2001. - vol. 121-B, №8. - p.961-966

2. Girard R., Gonzalez J. J., Gleizes A. Modeling of a two-temperature SF6 arc plasma during extinction // J. Phys. D: Appl. Phys. 1999. - v.32. - p. 12291238.

3. Gleizes A., Chervy В., Gonzalez J.J. Calculation of a two-temperature plasma composition: based and application to SF6 // J. Phys. D: Appl. Phys. — 1999. v.32. p.2060-2067.

4. Gonzalez J. J., Girard R., Gleizes A. Decay and post-arc phases of a SFe arc plasma: a thermal and chemical non-equilibrium model // J. Phys. D: Appl. Phys. -.2000. v.33. - p.2759-2768.

5. Кукеков Г. А. Выключатели переменного тока высокого напряжения // Энергия. Ленинградское отделение. 1972. - 337с.

6. Финкельнбург В., Меккер Г. Электрические дуги и термическая плазма. М.: Изд-во иностр. лит., 1961. - 370с.

7. Грановский В. Л. Электрический ток в газе. Установившийся ток. — М.:Наука, 1971. -544с.

8. Райзер Ю. П. Лазерная искра и распостранение разрядов. — М.: Наука, 1974.-308с.

9. Аронзон Н. А. О теоретическом обосновании принципа минимума напряжения дуги // Электричество. 1958. - №3. - с.56-60

10. Peters Th. Bogenmodelle und Steebeck's Minimumprinzip //Proc. 5th Int. Conf. Ionized Phenom. In Gases. V.l, Munich, Aug.28-Sept.l, 1961 Amsterdam, 1962. - p.885-896.

11. Заруди M. E. Методы расчета характеристик дуги в канале // ТВТ, — 1968 № 1. - с.35-43.

12. Заруди М. Е. Методы расчета столба цилиндрической дуги в канале с учетом излучения // Изв. СО АН СССР. Сер. техн. наук. 1967. — №3, вып.1. - с.8-14.

13. Maecker Н., Stablein Н. The channel model of the cylindrical arc // Proc. 10th Int. Conf. Phen. Ionized Gases, Oxford, Sept/13-18, 1971/ Oxford, 1971. -p.178.

14. Whitman A. M., Cohen I. M. Theory of the constant property arc // J. Appl. Phys. 1973.- v.44, №4 -p.l552-1556.

15. Swanson B. W., Roidt R. M. Boundary layer analysis of an SFe circuit breaker // IEEE Trans. On Power Apparatus and Systems. 1970. - v. PAS-90, №3.-p. 1086-1093.

16. Herman W., Kogelshatz U., Niemeyer L., Ragaller K., Shade E.

17. Experimental and theoretical study of a stationary high-current arc in a supersonic nozzle flow // J. Phys. D: Appl. Phys. 1974. - v.9. - p. 1703-1722.

18. Hermann W., Ragaller K. Theoretical Description of the current interruption in HV Gas blast Breakers // IEEE Trans. On Power Apparatus and Systems. 1977. - v. PAS-96, №5. - pp.1546-1555.

19. Blundell R. E., Fang M. Т. C. The similarity and scaling of radiating arcs burning in a turbulent, axially accelerating gas flow // J. Phys. D: Appl. Phys. — 1997. v.30. - p.628-635.

20. Yan J. D., Hall W. В., Fang M. Т. C. Experimental and theoretical investigation of an enclosed free burning arc in SF6// J. Phys. D: Appl. Phys. -2000. v.33. - p. 1070-1080.

21. Zang J. L., Yan J. D., Fang M. Т. C. Influence of cathode vaporization on arc characteristics in a supersonic nozzle // XIVth Symposium on Physics of Switching Arc, Nove Mesto na Morave, Czech Republic, 10-14 Sept., 2001. — p. 177-180.

22. Fang M. Т. C., Zhang J. F., Newland D.B. Theoretical investigation of a 2kA DC nitrogen arc in a supersonic nozzle // J. Phys. D: Appl. Phys. 1987. -v.20. -p.368-379.

23. Nimeyer L., Ragaller K. Development of turbulence by interaction of gas flow with plasma // Z. Naturforschung. 1973. - v.28a. - heft 8. - p.300-312

24. Бельков E. П. Восстановление электрической прочности искровых промежутков после протекания больших импульсов тока //ЖТФ, — 1974. — t.XLIV, №9. с. 1946-1951.

25. Джонс Дж. Влияние турбулентности на отключение тока // Отключение токов в сетях высокого напряжения // под ред. Рагаллера К. //пер. с англ. М.:Энергиздат, 1981 - с. 268-297.

26. Herman W., Kogelshatz U., Niemeyer L., Ragaller K., Shade E.1.vestigation on the Physical Phenomena around current zero in HV gas Blast Breakers // IEEE Trans. On Power Apparatus and Systems, 1976. - v. PAS-95, №4.-p. 1165-1176.

27. Yan J. D., Nuttall К. I., Fang M. Т. C. A comparative study of turbulence models for SF6 arcs in a supersonic nozzle // J. Phys. D: Appl. Phys. 1999 - v.32. - p.1401-1406.

28. Иевлев B.M. Турбулентное движение высокотемпературных сплошных сред М: «Наука», 1975. - 256с.

29. Бай Ши-и Динамика излучающего газа-М.: Мир, 1968. -323с.

30. Зельдович Я. Б., Райзер Ю. П. Физика ударных волн и высокотемпературных гидродинамических явлений М.: Наука, 1966, — 686с.

31. Низкотемпературная плазма. Т.1 Теория столба электрической дуги // под ред. Энгельшт В. С., Урюков А. Б., Новосибирск: Наука, 1990 — 376с.

32. Aubrecht V., Gross В. Net emission coefficient of radiation in SF6 arc plasmas // J. Phys. D: Appl. Phys. 1994. - v.27. - p.95-100

33. Frost L.S., Liebermann R.W. Composition and transport Properties of SF6 and their use in a simplified enthalpy flow arc model // Proceedings of the IEEE. — 1971. v.59, №4. — p.474-485

34. Betancourt S. Untersuchungen zur dielektrischen Wiederverfestigung magnetusch beblsener Lichtbogen in ruhendem SF6 bei Ausschaltstromen von 500 bis 2500 A // Diss. DI. Aachen. 1988.

35. Kesserling F., Koppelmann B. Das Schaltproblem der Hochspannungstechnik//Arc. f. El., 1936. - Bd.30., - h.l .

36. Лойцянский Л. Г. Механика жидкости и газа М.: Наука, 1987. -736с.

37. Moll H-R. Untersuchungen ueber die dielektrische Verfestigung bei SF6 Leistungsschaltern // Diss.arb.ETH Zuerich. 1983.

38. Лапин, Ю.В., Стрелец M. X. Внутренние течения газовых смесей. -М.: Наука, 1989.-366.

39. Монин А.С., Яглом A.M. Статистическая гидромеханика, т.1 // СПб: Гидрометеоиздат, 1992. 693с.

40. Fang М. Т. С., Zhuang Q., Guo X. J. Current zero behavior of an SFe gas-blast arc. Part II: turbulent flow // J. Phys. D: Appl. Phys. 1994. - v.21. - p.74-83.

41. Chervy В., Gleizez A., Razafinimanana M. Thermodynamic properties and transport coefficient in SF6-Cu mixtures at temperature of 300-30000K and pressures of 0.1-lMpa//J. Phys. D: Appl. Phys. 1994.- v.27. - p. 1139-1206

42. Варгафтик H. Б. Справочник по теплофизическим свойствам газов и жидкостей// М: Наука, 1972. 708с.

43. Nimeyer L., Ragaller К. Development of turbulence by interaction of gas flow with plasma // Z. Naturforschung. 1973. - v.28a. - heft 8. - p.300-312

44. Шнейдер M. H. О турбулентном охлаждении послеискрового канала в газе // РАН, 1993, Препринт №522. 26с.

45. Laufer J The structure of turbulence in fully developed pipe flow //NACA Tec. Rep. №1174. - 1954. - 12p.

46. Хинце И.О. Турбулентность //M.: ГИФМЛ. 1936. - 680с.

47. Klebanoff P.S. "Characteristics of turbulence in boundary layers with zero presser gradient" //NACA Report. 1955. - №1247. - 19p.

48. Бельков E. П. Охлаждение газа и восстановление электрической прочности полсе протекания искрового разряда //ЖТФ, 1969. - T.XLI, №8. -с.1678-1681.

49. Greig J. R., Pachacek R. E., Raleigh M. Channel cooling in turbulent convective mixing // Phys. Fluids. 1985. - v.28, №8. - p.2357-2364

50. Кабанов С. H., Маслова Л.И., Тархова Т.И., Трухин В.А., Юров

51. В.Т. Динамика остывания лазерной искры в воздухе // ЖТФ, 1990. - т.60, №6. - с.37-40.

52. Lee S.C., Harsha Р.Т. Use of turbulent kinetic energy in free mixing studies // AIAA, 1970. -т.8. - с. 1026-1032.

53. Турбулентность // под ред. Фрост У., Моулден Т. М: Мир, 1983. -536с.

54. Зельдович Я. Б. К теории горения не перемешанных газов // ЖТФ. -1949. -т.19, №10. -с.1199-1210.

55. Вулис Л. А. К вопросу о роли пульсаций температуры в турбулентном горении // Тр. III Всесоюзного совещания по теории газов М.: Изд-во АН СССР. - 1959.- с.86-90.

56. Турбулентные течения реагирующих газов // под ред. Либби П.А. — М.: Мир, 1983.-325с.

57. Кузнецов В. Р., Сабельников В. А. Турбулентность и горение М.: Наука, 1986.-286с.

58. Акатнов Н. И., Лавров В. А. О влиянии турбулентных пульсаций температуры и концентрации на физико-химические процессы в высокотемпературной турбулентной газовой струе // ТВТ. — 1978. — т. 16, №5. с.1005—1011.

59. Хожателев М. Б., Ярин Л. П. Экспериментальное исследование влияния пульсаций температуры на электропроводность низкотемпературной плазмы // ТВТ. 1967. -т.5, №6. - с.941-948

60. Либрович В. Б., Лисицын В. И. Взаимодействие пульсаций потока с химической реакцией турбулентного пламени // Институт проблем механики АН СССР, препринт №57, Москва, 1975. 43с.

61. Ловке Л. Перенос излучения в дугах отключения // Отключение токов в сетях высокого напряжения / под. ред. Рагаллера К. / пер. с англ. — М.: Энергиздат, 1981 с. 268-297.

62. Gleizes A., Rahmani В, Gonzalez J. J., Liani В Calculation of net emission coefficient in N2, SF6 and SF6 N2 arc plasmas // J. Phys. D: Appl. Phys.- 1991. -v.24. p. 1300-1309

63. Liebermann R.W., Lowke J. J. Radiation emission coefficients for sulfur hexafluoride arc plasmas // J. Quant. Spectrosc. Radiat. Transfer- 1976. v. 16.- p.253-264.

64. Дресвин С. В. и др. Физика и техника низкотемпературной плазмы. -М.: Атомиздат, 1972. 352с.

65. Адрианов В.Н. Основы радиационного и сложного теплообмена М.: Энергия, 1972.-256с.

66. Назаренко И. П., Паневин И. Г. Расчет характеристик стабилизированных дуг с учетом переноса излучения и отрыва температур // Моделирование и методы расчета физико-химических процессов в низкотемпературной плазме М.: Наука, 1974. - с.86-107

67. Lowke J. J., Capriotti E. R. Calculation of temperature profiles of high pressure electric arcs using the diffusion approximation for radiation transfer // JQSRN. 1969. - v.9. - p.207-236

68. Lowke J. J. Characteristic of radiation-dominate electric arc // J. Appl. Phys. 1970. - v.41, №6. - p. 2588-2600

69. Флетчер К. Вычислительные методы в динамике жидкостей, т.2 // М.: Мир, 1991.-552с.

70. Бим Р. М., Уорминг Р. Ф. Неявная факторизованная разностная схема для уравнений Навье-Стокса течений сжимаемого газа // РТК. 1978. - т. 16, №4.-с. 145-156

71. Чакраварти С.И. Жем К.-Й. Расчет трехмерных сверхзвуковых течений с дозвуковыми зонами на основе уравнений Эйлера // Аэрокосмическая техника. 1987. - №11. - с.22-35.

72. Roe P. L., Approximate Reimann Solvers, Parameter Vectors and Difference Scheme //J. Comput. Phys. 1981. - v.43. - pp.357-372

73. Hirsch C. Numerical Computation of Internal and External flows, v.2 // Wiley, 1990.-736c.

74. Kermani M. J., Plett E. G. Roe Scheme in Generalized Coordinates; Part I Formulations // 39th AIAA Aerospace Science Meeting and Exhibit 8-11 January 2001, Reno, NV, AIAA 2001-0086, 1 lp.

75. Симуни JI. M. Численное решение задачи о неизотермическом движении вязкой жидкости в плоской трубе // Инженерно-физический журнал. 1966. -т.Ю, №1. - с.86-94.

76. Манина М.П., Поспелов В. А., Ходжиев С. Решение прямой задачи сопла Лаваля в приближении «узкого канала» методом установления // Гидроаэродинамика. Сборник научных трудов. Л.: изд. ЛПИ, -1983. — с.26— 30

77. Пирумов У. Г. Обратная задача теории сопла // М.: Машиностроение, 1988,-240с.

78. Belhaouari J. В., Gonzalez J. J., Gleizes A. Simulation of a decaying SF6 arc plasma: hydrodynamic and kinetic coupling study // J. Phys. D: Appl. Phys. -1998. v.31. - p.1219-1232.

79. Frind G., Rich J.A. Recovery speed of axial flow gas blast interrupter: dependence on pressure and dl/dt for air and SF6 // IEEE Trans, on Power Apparatus and Systems. 1974. - v.PAS-93., №5 - pp. 1657-1665.

80. Perkins J.F., Frost L.S. Effect of nozzle parameters on SF6arc interruption // IEEE Trans. On Power Apparatus and Systems, 1973. - v. PAS-92, №3. - p. 961-970.

81. Басин А. С., Ипатьева О. С., Попов В. Н. Моделирование температурных полей в трубчатом электроде плазмотрона отнестационарного воздействия пятна дуги // Теплофизика и аэромеханика. -1998. т.5, №4. - с.583-592

82. Moll R., Schade Е. Dielectric Recovery of Axially Blown SF6-Arcs // 6 Int. Conferenece on Gas Discharge and Their Applications, 1950. - p. 13-16.

83. Бэтчелор Дж. Введение в динамику жидкости и газа // М.: Мир, 1973, -758с.

84. Вывод зависимости /?V = С,/дг824Р

85. Далее заметим, что ^(^У получим требуемую формулу.

86. Моль с импульсом (ри) 1 из области 1, будучи перенесен в слой с радиусом г,создаст отрицательную пульсацию импульса:ри\-(ри) = {р'и\ (А.6)

87. Введем дополнительное обозначение: ри=С\4р=Дг). Тогда:r¥\ . 1 а*/pu\-(pu) = f(r-l)-f(r) = Рри\ (ри) = f(r + /) - Дг) =rdf кдг1 +1 д2/ ,22 дг21. А.7) (А.8)

88. Складывая правые стороны (А.5), (А.6) и (А.7), (А.8) получим:pV^ + C/y^fr'2or

89. Таким образом, учтя вид функции Дг), окончательно имеем:1 дг

90. Матрицы R\ L и D (к разделу 2.2).1. Вид матриц для А =dF 8U1. D =а{и + с) О О О1. О аи О О1. О О аи О1. О 0 0 а(и-с)1. R =1 и+с v h1 и О О0 0 101 и — с v h1.=