РЕФЕРАТ
Отчет содержит 114 стр., 104 форм., 96 рис., 17 табл., 32 ист., 4 прил.
Ключевые слова: ЛАЗЕРНЫЙ РАЗРЯД, ОПТИЧЕСКИЙ РАЗРЯД,
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ, ДВУХТЕМПЕРАТУРНАЯ МОДЕЛЬ,
РАВНОВЕСНАЯ МОДЕЛЬ, ОДНОМЕРНАЯ МОДЕЛЬ, ВОЗДУШНАЯ
ПЛАЗМА, АРГОНОВАЯ ПЛАЗМА, ПЛАЗМОХИМИЧЕСКИЕ ПРОЦЕССЫ,
УРАВНЕНИЕ
ЭНЕРГИИ,
РАВНОВЕСНОСТИ
СЕЧЕНИЕ
СТОЛКНОВЕНИЙ,
ВЫСОКОСКОРОСТНОЙ
ЭКСПЕРИМЕНТАЛЬНЫЙ
СТЕНД,
ГАЗОВЫЙ
ОПТИЧЕСКАЯ
ОЦЕНКА
ПОТОК,
СХЕМА,
ИМПУЛЬСНЫЙ ЛАЗЕР, НЕПРЕРЫВНЫЙ ЛАЗЕР.
Разрабатываются математические модели непрерывного оптического
разряда, поддерживаемого лазерным излучением. На основе полученных
моделей и результатов предлагается проект стенда для проведения
экспериментов. Результатами работы являются: одномерная равновесная
модель лазерного разряда в высокоскоростном газовом потоке; предложенная
математическая модель неравновесного (двухтемпературного) разряда;
методика экспериментального исследования исследуемой темы, позволяющая
проверить результаты моделирования; оптические схемы установки для
регистрации пробоя и заведения излучения, а также представлена компоновка
всех элементов в лаборатории.
2
СОДЕРЖАНИЕ
ПЕРЕЧЕНЬ СОКРАЩЕНИЙ И ОБОЗНАЧЕНИЙ ..................................... 5
ВВЕДЕНИЕ ..................................................................................................... 6
1
Теоретическая справка ............................................................................... 9
2
Математическая модель равновесного стационарного разряда в
высокоскоростном газовом потоке. ................................................................ 13
2.1
Методика расчета ............................................................................... 14
2.2
Выбор параметров .............................................................................. 16
2.3
Первоначальная оценка равновесности плазмы ................................. 17
2.4
Аппроксимационные формулы........................................................... 19
2.4.1
Коэффициент поглощения ............................................................ 19
2.4.2
Молекулярный вес и энтальпия .................................................... 28
2.4.3
Теплоемкость ................................................................................ 30
2.4.4
Теплопроводность......................................................................... 31
2.5 Одномерная модель .................................................................................. 33
2.6 Коды процедур ......................................................................................... 38
2.6.1
Расчет надежности программы ..................................................... 41
2.6.1.1 Метод Мусса ................................................................................. 41
2.6.1.2 Расчет по методу Мусса ................................................................ 42
2.7
3
Выводы по пункту 2 ........................................................................... 44
Моделирование неравновесного непрерывного оптического разряда в
высокоскоростном газовом потоке в двухтемпературном приближении ...... 45
3.1
Кулоново взаимодействие одиночного электрона с положительным
ионом. .......................................................................................................... 46
3.2
Оценка равновесности ........................................................................ 57
3.3
Математическая модель непрерывного лазерного разряда в
высокоскоростном газовом потоке в двухтемпературном приближении ... 61
3
4
3.3.1
Уравнение энергии тяжелых частиц: ............................................ 61
3.3.2
Уравнение энергии электронов:.................................................... 66
3.4
Плазмохимические процессы ............................................................. 72
3.5
Выводы по пункту 3 ........................................................................... 75
Методика экспериментального исследования .......................................... 76
4.1
Оптическая схема ............................................................................... 77
4.2
Источники лазерного излучения......................................................... 81
4.3
Расчет коллиматоров .......................................................................... 84
4.4
Техническое описание СХПК ............................................................. 88
4.4.1
Две подсистемы СХПК ................................................................. 89
4.4.2
Пневмогидравлическая система СВД ........................................... 91
4.5
4.5.1
Манометры.................................................................................... 93
4.5.2
Преобразователи давления ........................................................... 97
4.6
5
Погрешности измерительных приборов в ПГС .................................. 93
Вывод по пункту 4 ............................................................................ 102
Компоновка экспериментальной установки в лаборатории. .................. 103
ЗАКЛЮЧЕНИЕ .......................................................................................... 107
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ................................. 108
Приложение А ............................................................................................. 111
Приложение Б ............................................................................................. 112
Приложение В ............................................................................................. 113
Приложение Г ............................................................................................. 114
4
ПЕРЕЧЕНЬ СОКРАЩЕНИЙ И ОБОЗНАЧЕНИЙ
В настоящем отчёте о ВКР применяются следующие сокращения и
обозначения:
СВЧ – сверхвысокая частота
FPC – Free Pascal Compiler (бесплатно распространяемый компилятор
языка программирования Pascal)
ЛС – лазерная система
СХПК – система хранения и подачи компонентов
ПГС – пневмогидравлическая система
СВД – система высокого давления
ПГГ – парогазогенератор
УВ – устройство воспламенения
ПЭК ДД – электропневмоклапан двойного действия
5
ВВЕДЕНИЕ
В настоящее время технические характеристики и требования к технике,
которая связана с полетами со сверхзвуковой скоростью возрастают. Способы,
которые дают возможность воздействовать на поток для улучшения
аэродинамических характеристик, являются наиболее значимыми в данной
области.
Актуальность и применяемость работы. Существует множество
способов подвода энергии в скоростном потоке, в большинстве случаев
происходит формирование протяженной области неоднородности – теплового
следа, который является основным элементом воздействия.
Результаты многочисленных исследований показали, что тепловой след
может привести к перестройке ударно-волновой структуры, изменить
аэродинамические моменты и способствовать снижению аэродинамического
сопротивления. Наиболее технологичный и удобный для практического
применения способ воздействия на сверхзвуковой поток - ввод энергии
лазерного излучения. Организуется оптический пробой газа и поглощение
энергии излучения созданной плазмой.
Поддержание оптического разряда в высокоскоростном газовом потоке
позволит заменить применяемые в данный момент конструктивные элементы
(из-за которых образуются в потоке скачки уплотнения) в гиперзвуковых
летательных аппаратах, которые нужны для зажигания топлива, с помощью
локализации плазмы и поддержания ее лазерным излучением появляется
возможность управлять поджогом, а также влиять на сам поток, тем самым
увеличив существующие пределы скоростей.
Сам импульсный разряд и фундаментальная проблема взаимодействия
лазерного излучения на атмосферный воздух, движущийся газ и др. были
изучены в различных научных трудах, которые помогали в освоении данной
темы.
Одним из
направлений исследования фундаментальной задачи
взаимодействия лазерного излучения с атмосферным воздухом является
6
влияние движения газа на образование оптического разряда, поглощение
энергии в нем, взаимодействие оптического разряда с ударной волной [1 -3].
Применение
лазерного
сверхзвукового потока
излучения
для
газа является
воздействия
на
структуру
перспективной задачей [4-5],
исследуемой и в рамках аэрокосмических разработок [6].
Основная цель работы – исследовать возможность непрерывного
лазерного разряда, локализированного в высокоскоростном газовом потоке,
предложить математические модели, а также проект экспериментального
стенда, с помощью которого можно будет применять теорию на практике. В
процессе углубления в данную тему возникали задачи, которые следовало
решить, так основной задачей было предложить работающую математическую
модель плазмы в равновесном (п. 2) и двухтемпературном (п. 3)
приближениях. Данная задача делилась на множество подзадач, решения
которым нужно было найти.
Апробация результатов. Результаты работы были опубликованы в
научных журналах:
−
«Интернаука» 22 (198), часть 2 (ISSN: 2687-0142), стр. 92-95
(DOI:10.32743/26870142.2021.22.198.288416).
Статья:
Определение
траекторий движения одиночного электрона вокруг иона в полностью
ионизированной плазме,
−
«Студенческий вестник» 21 (166), часть 8 (ISSN: 2686-9810),
стр. 99-103. Статья: Одномерная равновесная модель непрерывного
оптического разряда в высокоскоростном газовом потоке,
−
Сборник научных статей по материалам V – международной
научно-практической конференции «Вестник науки», фундаментальные и
прикладные исследования: инноватика в современном мире, К-196-0,
стр. 8-19. Статья: Разработка одномерной модели непрерывного оптического
разряда в высокоскоростном равновесном потоке воздуха.
В настоящей работе были исследованы возможности получения
стационарного оптического разряда в высокоскоростном газовом потоке.
7
Сформулирована математическая модель стационарного разряда (для воздуха
в
условиях
равновесия,
для
аргона
приведена
неравновесная
двухтемпературная модель), приведены оценки параметров лазерного
излучения, необходимого для поддержания стационарного разряда. По
полученным математическим моделям написаны коды процедур для
получения всех значений параметров, нужных для описания данного процесса.
Приведена принципиальная оптическая схема для регистрации теневой
картины. Представлена обстановка эксперимента, предложен подвод
непрерывного излучения.
8
1 Теоретическая справка
Эффект оптического пробоя был обнаружен еще в 1963 году такими
учеными, как П.Мейкер, Р.Терхун и У.Сэвидж [7]. Само обнаружение данного
эффекта стало доступным благодаря созданию лазера с модулированной
добротностью, который дает очень мощный импульс (длительность ±30 нс,
энергия 1 Дж, пиковая мощность 30 МВт). Луч данного рубинового лазера был
сфокусирован линзой, в комнатном воздухе в области фокуса вспыхнула искра
и там образовалась плазма.
В последствии, был теоретически предсказан и получен на опыте в 1970
году непрерывный оптический разряд [8]. По сравнению с традиционными
способами поддержания плазмы с температурой примерно 10000 К при
помощи дугового, индукционного, СВЧ-разрядов, для подвода энергий к
плазме оптическим способом не требуется конструктивных элементов:
электродов, индуктора, волновода. Световая энергия свободно передаётся на
расстояние световым лучом. Это открывает возможность зажигания плазмы
на расстоянии от лазера и в любых, даже труднодоступных местах. Если
продувать холодный газ через горящий непрерывный оптический разряд,
подобно тому, как это делается в дуговых и прочих генераторах непрерывной
плазменной струи - плазмотронах, получается оптический плазмотрон (рис. 1)
Рисунок 1 - Принципиальная схема оптического плазмотрона: 1 - контуры
светового канала; 2 - линза; 3 - плазма; F - точка фокуса; и - газовый
поток; РО - лазерное излучение. [9]
9
Оптический разряд – вид высокочастотного разряда в газах, который
наблюдается в воздухе или других газах, в частотах излучения оптического
диапазона. Инициирование оптического разряда обычно создается мощным
лазерным излучением. Различают два основных вида оптического разряда:
•
Оптический пробой (лазерная искра) – этот эффект может
наблюдаться при фокусировки достаточно мощных лазерных импульсов.
Величина интенсивности порядка 1011 Вт/см2 является необходимой для
получения данного эффекта в воздухе. Происходит бурное нарастание
ионизации ранее не ионизированного газа. [10]
•
Непрерывный оптический разряд – стационарный газовый разряд,
который с помощью лазерного излучения поддерживается в созданной ранее
относительно плотной плазме, то есть существующее ионизированное
состояние газа поддерживается под действием светового излучения. [1 1]
Пороговое поле, при котором наблюдается оптический разряд, зависит
от многих факторов, таких как:
•
Род газа
•
Давление
•
Длительность импульса
•
Размер фокусного пятна
И др.
Пороговая интенсивность снижается при не очень больших давлениях.
Основой для оптического пробоя является эффект развития электронной
лавины, однако при этом затравочные электроны появляются в результате
многофотонной ионизации молекул или атомов газа. При взаимодействии
лазерного излучения с атомом происходит поглощение нескольких (порядка
10—20) фотонов, сопровождающееся вырыванием электрона. Электрон
ускоряется в поле лазерной волны и при соударении с другими атомами
производит их ионизацию, рождая ещё один электрон. Затем уже два
электрона ускоряются полем и при столкновениях с атомами рождают ещё два
10
электрона. Таким образом происходит лавинообразное нарастание числа
свободных электронов. [10]
Для того, чтобы произошёл пробой необходимо, чтобы за время
действия лазерного импульса родилось достаточно большое количество
электронов.
Электронная лавина - нарастающий процесс размножения электронов в
результате ионизации атомов и молекул, как правило, электронным ударом;
является главнейшим элементом электрического пробоя газов. В большинстве
случаев электронная лавина развивается в электрическом или электромагнитном поле, хотя возможно лавинное размножение электронов чисто
тепловой природы, например в ударной волне. [10]
Для создания ионизированного поля (плазмы) стоит упомянуть, что для
процессов ионизации и диссоциации молекул необходимо добиться высокой
энергии, которая достигается путем нагрева среды сфокусированным
лазерным пучком. Для получения пробоя газа (оптического разряда) с
помощью короткого импульса нужно достичь высоких температур в короткий
промежуток времени, при которых возникает процесс образования ионов из
нейтральных атомов и молекул. Положительно заряженный ион образуется,
если электрон в
молекуле
получает
достаточную
энергию
преодоления потенциального
барьера,
равного наименьшей
для
энергии,
необходимой для удаления электрона от свободного атома в его низшем
энергетическом (основном) состоянии на бесконечность. Отрицательно
заряженный ион, наоборот, образуется при захвате дополнительного
электрона атомом с высвобождением энергии.
Классифицировать разрядные процессы можно на 3 категории:
1. Пробой газа, развитие в нем бурной лавинной ионизации под
действием приложенного внешнего поля, превращение неионизованного
вначале газа в плазму. (Оптический разряд)
2. Поддержание полем неравновесной плазмы, в которой температура
электронов, обеспечивающих ионизацию, достаточно высока, а газ атомов,
11
молекул, ионов остается холодным. Обычно это — слабоионизованная плазма
при не слишком высоких давлениях, ниже сотни Торр. Степень ионизации при
этом гораздо меньше равновесной величины, отвечающей температуре
электронов.
3. Поддержание полем равновесной плазмы, в которой температуры
электронов и тяжелых частиц близки друг к другу, а степень ионизации близка
к термодинамически равновесной. Это — так называемая низкотемпературная
плазма с температурами порядка 10000 К, обычно при давлениях порядка
атмосферного. (Непрерывный оптический разряд) [11].
12
2
Математическая модель равновесного стационарного разряда в
высокоскоростном газовом потоке.
В данном пункте рассмотрена равновесная модель непрерывного
оптического разряда. Равномерный стационарный поток воздуха набегает на
неподвижное
фокальное пятно
непрерывного лазерного излучения;
исследуется возможность получения непрерывного оптического разряда в
высокоскоростном потоке в одномерном приближении. Была получена
математическая модель стационарного разряда, при определенных параметрах
лазерного излучения, необходимых для достижения поставленной цели.
Основаниями
построения
равновесной
модели
послужили
первоначальные оценки равновесности, которые показали, что в данном
случае плазму можно считать равновесной.
13
2.1
В
Методика расчета
работе рассмотрена одномерная геометрия
взаимодействия
изобарического потока воздуха с фокальным пятном лазерного излучения. В
данном приближении параметры излучения и потока зависят только от
продольной координаты. Было принято, что только теплопроводность
осуществляет теплоотвод от горячей зоны пробоя. Тогда уравнение энергии
будет иметь следующий вид:
𝑑ℎ
𝑑2 𝑇
𝜌𝑣
= 𝛼𝐼 + 𝜆 2
𝑑𝑥
𝑑𝑥
(2.1)
где 𝜌 – плотность потока; 𝑣 – скорость потока ℎ - энтальпия единицы массы 𝑥
– координата 𝛼 – коэффициент поглощения 𝐼 – интенсивность излучения 𝜆 теплопроводность 𝑇 – температура.
Было принято, что параметры потока, после прохождения лазерного
пятна, в данной одномерной задаче остаются постоянными, следовательно на
выходе из потока были поставлены граничные условия и введено численное
интегрирование навстречу потоку.
𝑥 = 𝑥𝑜𝑢𝑡 : 𝑇 = 𝑇𝑜𝑢𝑡 ;
Физическим
параметром,
𝑑𝑇
=0
𝑑𝑥
определяющий
(2.2)
решение,
является
температура на выходе 𝑇𝑜𝑢𝑡 .
Зная, что 𝑑ℎ = 𝑐𝑝 (𝑇)𝑑𝑇 можно записать уравнение (1) в виде:
𝑑
𝜆𝑌 = 𝜌𝑣𝑐𝑝 𝑌 − 𝛼𝐼
𝑑𝑥
𝑑𝑇
=𝑌
𝑑𝑥
14
(2.3)
Система уравнений (3) интегрируется навстречу потоку, от 𝑥 = 𝑥𝑜𝑢𝑡 до
𝑥 = 𝑥𝑖𝑛 < 𝑥𝑜𝑢𝑡 с «начальными» условиями (2).
В (3) входят параметры воздушной плазмы: 𝑐𝑝 ≡ 𝑑ℎ/𝑑𝑇, 𝛼, 𝜆, 𝜇.
Теплоемкость 𝑐𝑝 для равновесной воздушной плазмы можно определить
дифференцированием энтальпии по температуре при постоянном давлении, а
энтальпия до 20kK вычисляется по процедуре, описанной в [13].
Молекулярный вес 𝜇 входит в уравнение состояния: 𝜌 = 𝜇𝑃/(𝑅0 𝑇) и также для
равновесной плазмы вычисляется в зависимости от температуры и давления
согласно [13].
15
2.2
Выбор параметров
Были рассмотрены следующие варианты параметров набегающего
потока:
Таблица 1 – Выбор параметров
Все значения параметров (давление температура), которые были
изначально взяты, соответствуют стандартной земной атмосфере на
соответствующей высоте. Предполагается, что перед стационарным лазерным
разрядом в набегающем потоке возникает ударная волна, параметры за
которой приведены в таблице с пометкой «^». Рабочий вариант, в данной
работе был взят по варианту 2. Для получения одномерной модели при
параметрах 1 атм, 10 м/с, 300 К и 1 атм, 300 м/с, были вручную изменены
значения параметров p^, v^, T^, после чего получена модель при таких
параметрах.
16
2.3
Первоначальная оценка равновесности плазмы
Если лазерный разряд формируется в фокусе лазерного пучка, то время
пребывания высокоскоростного потока в пределах фокального пятна
достаточно мало. Возникает вопрос: насколько существенны неравновесные
эффекты в такой плазме и насколько правомерно применение модели
равновесной плазмы для описания стационарного лазерного разряда в
высокоскоростном воздушном потоке? Энергия лазерного излучения
воспринимается электронами плазмы и передается ионам через электронионные столкновения. Оценим характерное время установления равновесия
между электронами и ионами. Пусть газокинетический диаметр ионов
составляет 𝑑 =3A. Тогда сечение электрон-ионных столкновений 𝜎 = 𝜋𝑑 2 /
4 =7.1×10 -20 м2. В действительности сечение больше из-за действия
дальнодействующих кулоновых сил, так что приведенную оценку можно
рассматривать как нижнюю. Средняя длина свободного пробега электронов
между столкновениями с ионами: 𝜆 = 1/(𝑛𝜎),
где 𝑛 = 𝑃/(2𝑘𝑇) –
концентрация ионов в полностью однократно ионизированной равновесной
плазме. Для давления 𝑃 =3 атм (вариант 2 таблицы 1), согласно [12],
температура полной однократной ионизации – около 20kK. Для этих условий
𝑛 =5.4×10 23 1/м3, 𝜆 =2.6×10 -5 м. Среднеквадратическая скорость теплового
движения равновесных электронов 𝑣𝑒 = √3𝑘𝑇/𝑚𝑒 =9.5×105 м/с, частота
электрон-ионных столкновений 𝑓𝑒 = 𝑣𝑒 /𝜆 =3.7×1010 Гц. Если диаметр
фокального пятна 𝑑 =40 мкм, то при скорости потока за скачком уплотнения
𝑣 =295 м/с (вариант 2) число электрон-ионных столкновений за то время, пока
газ находится в пределах фокального пятна, составляет 𝑁 = 𝑓𝑒 𝑑/𝑣 =5000. При
каждом столкновении от электрона к иону передается часть энергии,
пропорциональная соотношению масс электрона и иона 𝑚𝑒 /𝑀, а число
столкновений, необходимое для того, чтобы энергия горячих электронов
уменьшилась в 𝑒 раз, составляет 𝑁1/𝑒 = 𝑀/(4𝑚𝑒 ) =6670. Тогда число
Дамкеллера, определяемое как 𝐷𝑎 = 𝑁1/𝑒 /𝑁, т.е. соотношение между числом
17
столкновений, необходимым для установления равновесия, и действительно
имеющим место числом столкновений, равно ≈1.3. Принимая во внимание, что
размер области разряда должен быть больше размера фокального пятна,
получаем, что в описанных условиях плазму стационарного лазерного разряда
можно приближенно считать равновесной. Во всяком случае, для оценок
возможности
осуществления
стационарного
лазерного
разряда
в
высокоскоростном воздушном потоке такой подход можно принять за основу.
18
2.4
Аппроксимационные формулы
2.4.1 Коэффициент поглощения
Расчет коэффициента поглощения производился по точкам из [1 3].
Первоочередно, выбраны данные при частотах 9250 и 9500 см -1.
Для более точной аппроксимации был выбран диапазон давлений, при
котором аппроксимации будут хорошо работать, так в таблице (2) можно
видеть выделенные значения коэффициента поглощения из [13] для частоты
9250 см-1, по которым была сделана аппроксимация. При температуре
T<=10000 К аппроксимационная зависимость была построена в диапазоне
давления от 0,003 до 10 атм, а при температуре T>=10000 К от 0,1 до 10. Таким
образом были получены графики (2-19), на которых есть аппроксимационные
уравнения для коэффициентов поглощения в зависимости от выше указанных
диапазонов давлений при T=const.
После получения формул коэффициента поглощения от давления для
каждой температуры, был определен общий вид уравнения, как формула (2.4)
α(p, T) = A(T) × p B(T)
где
(2.4)
α(p, T) – коэффициент поглощения, зависящий от давления и
температуры, A(T) – коэффициенты, полученные при аппроксимации
коэффициента поглощения по давлению при Т=const, умножающийся на
давление в степени B(T), зависит от температуры, p – давление, B(T) –
коэффициент степени при давлении, которые были так же получены при
аппроксимации коэффициента поглощения по давлению при Т=const .
Далее были аппроксимированы A(T) и B(T), рисунки (20-21). Эти
аппроксимационные формулы используются в основной формуле (2.4).
19
Таблица 2 – значения коэффициента поглощения для частоты 9250 см -1. [13]
Рисунок 2 - α(p) при T=2000 К
Рисунок 3 - α(p) при T=2500 К
Рисунок 4 - α(p) при T=3000 К
Рисунок 5 - α(p) при T=3500 К
Рисунок 6 - α(p) при T=4000 К
Рисунок 7 - α(p) при T=5000 К
20
Рисунок 8 - α(p) при T=6000 К
Рисунок 9 - α(p) при T=7000 К
Рисунок 10 - α(p) при T=8000 К
Рисунок 11 - α(p) при T=9000 К
Рисунок 12 - α(p) при T=10000 К
Рисунок 13 - α(p) при T=11000 К
Рисунок 14 - α(p) при T=12000 К
Рисунок 15 - α(p) при T=13000 К
21
Рисунок 16 - α(p) при T=14000 К
Рисунок 17 - α(p) при T=16000 К
Рисунок 18 - α(p) при T=18000 К
Рисунок 19 - α(p) при T=20000 К
Рисунок 20 – log10A(T) для температур от 2000 до 20000 К
22
Рисунок 21 – B(T) в диапазоне температур от 2000 до 20000 К
Рисунок 22 – Наложенные значения из [13] и работающей аппроксимацией
для данных значений.
Таким образом, мы можем посчитать коэффициент поглощения для
частоты 9250 см-1 в диапазоне давлений от 0,03 до 10 атм при Т<=10000 К и в
диапазоне от 0,1 до 10 атм при Т>=10000К. Выше приведен график (22), где
наглядно показаны значения коэффициента поглощения для фиксированной 1
23
атм по [13], а также значения, полученные с помощью аппроксимационных
формул.
Процесс получения модели расчетов для частоты 9500 см -1 аналогичен,
ниже (23-42) приведены все аппроксимации по данной частоте.
Значения коэффициента поглощения брались в диапазоне давлений от
0,03 до 10 атм при T<12000 K, и при температуре T>=12000 K от 0,3 до 10 атм.
Таблица 3 – значения для коэффициента поглощения для частоты 9500 см -1
[13].
Рисунок 23 - α(p) при T=2000 К
Рисунок 24 - α(p) при T=2500 К
24
Рисунок 25 - α(p) при T=3000 К
Рисунок 26 - α(p) при T=3500 К
Рисунок 27 - α(p) при T=4000 К
Рисунок 28 - α(p) при T=5000 К
Рисунок 29 - α(p) при T=6000 К
Рисунок 30 - α(p) при T=7000 К
25
Рисунок 31 - α(p) при T=8000 К
Рисунок 32 - α(p) при T=9000 К
Рисунок 33 - α(p) при T=10000 К
Рисунок 34 - α(p) при T=11000 К
Рисунок 35 - α(p) при T=12000 К
Рисунок 36 - α(p) при T=13000 К
26
Рисунок 37 - α(p) при T=14000 К
Рисунок 38 - α(p) при T=16000 К
Рисунок 39 - α(p) при T=18000 К
Рисунок 40 - α(p) при T=20000 К
Рисунок 41 – log10(A(T)) для температур от 2000 до 20000 К
27
Рисунок 42 – B(T) для температур от 2000 до 20000 К
Таким образом могут быть вычислены значения коэффициента
поглощения для частот 9250-9500 см-1, давлений 0,03 – 10 атм и температур
2000-20000 К.
Длина волны лазерного излучения равна 1064 нм, исходя из этого
частота равна примерно 9398,5 см -1. Следовательно для получения более
близких к реальности значений стоит выбрать модель для получения
коэффициента поглощения при частоте 9500 см-1.
2.4.2 Молекулярный вес и энтальпия
Молекулярный вес равновесной воздушной плазмы вычисляется
согласно [12]. Считается, что исходный состав воздуха – 21% О2 и 79% N2,
причем влиянием оксидов азота пренебрегается. Это достаточно обосновано,
т.к. оксиды азота существуют в ограниченном диапазоне температур (8-10 kK),
хотя и могут в этом диапазоне влиять на процессы ионизации. Результаты
вычислений приведены на рисунке (43) (на фоне графика, заимствованного из
[12], точками показаны результаты вычислений). Таким образом может быть
вычислен молекулярный вес воздуха в диапазоне температур до 20 kK и
давлений от 0.001 до 1000 атм. Зависимости молекулярного веса от
температуры обнаруживают ряд характерных особенностей, имеющих ясный
физический смысл. Особенно это заметно при низких давлениях. Так, для
P=0.001 атм в диапазоне T=3-4kK имеется «полка», соответствующая полной
диссоциации кислорода O2->2O. При T=6-7kK - такая же особенность,
связанная с диссоциацией азота. Молекулярный вес на этой «полке» ровно в 2
раза меньше исходного веса воздуха, т.к. все двухатомные молекулы
полностью диссоциированы. При T=12-17kK – полная однократная ионизация,
при которой молекулярный вес уменьшается еще вдвое, поскольку масса
электрона пренебрежимо мала по сравнению с массами атомов кислорода и
азота. При дальнейшем повышении температуры молекулярный вес снижается
дальше, что говорит о начале многократной ионизации.
28
Рисунок 43 – Точки среднего молекулярного веса при разных давлениях по
формулам [12], наложенные на график из модели [12]
Рисунок 44 – Точки удельной массовой энтальпии при разных давлениях по
формулам из [12], наложенные на график из модели [12] (красная линия –
Р=3 атм)
29
2.4.3 Теплоемкость
Поскольку
𝑐𝑝 = (𝜕р/𝜕𝑇)𝑝 ,
теплоемкость
можно
вычислить
дифференцированием (45, 46).
Рисунок 45 – Зависимость теплоемкости [Дж/(кг*К)] воздуха от температуры
при Р=1 атм
Рисунок 46 – Зависимость теплоемкости [Дж/(кг*К)] воздуха от температуры
при Р= 3 атм (x=T, y=log10(Cp))
30
Аппроксимационная зависимость теплоемкости от температуры для
давлений Р=1 атм и Р=3 атм приведены на соответствующих рисунках (45, 46).
2.4.4 Теплопроводность
Рисунок 47 – Апроксимационная модель для расчета теплопроводности при 1
и 3 атм (x=T; y= 𝜆)
Вязкость воздуха определяется по модифицированной формуле
Сатерленда до температур 9000 К [14]
𝜇
𝑇 1/2 1 + 111/𝑇∞
1
1.42
(
),
=( )
+
𝜇∞
𝑇∞
1 + 111/𝑇 1 + 0.89𝑐𝐴 /𝑐𝑀 1 + 1.26𝑐𝑀 /𝑐𝐴
(2.5)
сМ и сА – концентрация молекул и атомов в высокотемпературном частично
диссоциированном воздухе, вычисляемый согласно [12]. При давлениях
свыше 0.1 атм и температурах меньше 9000К степень ионизации не превышает
3% - это позволяет не учитывать влияние электронов на теплопроводность.
31
Принимая число Прандтля в этом диапазоне температур постоянным и равным
0.72, получаем:
𝑐𝑝 𝜇∞ 𝑇 1/2 1 + 111/𝑇∞
1
1.42
( )
(
)
𝜆=
+
𝑃𝑟 𝑇∞
1 + 111/𝑇 1 + 0.89𝑐𝐴/𝑐𝑀 1 + 1.26𝑐𝑀 /𝑐𝐴
(2.6)
На рисунке (47) эта зависимость представлена кривой 7.
При температурах, превышающих 8000-9000К, на теплопроводность
значительное влияние оказывает ионизация [15], поскольку при наличии
свободных электронов тепловая энергия переносится в основном именно
электронами. Приведены данные по теплопроводности ионизированной
равновесной воздушной плазмы в диапазоне температур от ~10kK до 20 kK и
давлений от 0.001 до 10 атм. Эти данные показаны на рисунке (47) кривыми 15. Для P=3 атм данные получены интерполяцией (кривая 6). На рисунке (47)
существует
значительный
разрыв
между
данными
вязкости
и
теплопроводности. Данные в данном промежутке отсутствуют, в связи с этим
была искусственно внедрена интерполяционная зависимость. Таким образом
на рисунке (47) показаны зависимости, которые позволяют оценить
теплопроводность воздуха в условиях равновесия для температур до 20000 К,
давлений 1 и 3 атм. Для данных давлений была использована одна
аппроксимационная модель, так как вязкость при данных давлениях
незначительно меняется, так же можно заметить на рисунке (47), что кривые
находятся близко друг к другу, а аппроксимационная зависимость для 1 атм
так же незначительно отличается от зависимости для 3 атм, следовательно в
основной модели мы можем использовать одинаковую зависимость для
давлений 1 атм и 3 атм.
32
2.5
Одномерная модель
Таким образом, имеется вся информация о свойствах воздуха при
давлениях P=3 и P=1 атм и температурах до 20 kK, необходимая и достаточная
для интегрирования уравнений (3).
Приведем численные решения уравнений (3) с граничными условиями
(2), полученные интегрированием назад (против потока) для P=3 и P=1 атм.
Задавая различные значения температуры на выходе потока из области фокуса
лазерного пучка, получаем температуру на входе. Профиль интенсивности
излучения – гауссов:
𝐼
𝑥 − 𝑑/2 2
) )
= 𝑒𝑥𝑝 (− (
𝐼0
𝑑/2
(2.7)
𝑑 – диаметр пятна, 𝐼0 = 4𝑃/𝜋𝑑 2, 𝑃 – мощность непрерывного лазерного
излучения.
Рисунок 48 – Зависимость выходной температуры от входной при
определенных диаметрах пятна для давления 1 атм и скорости 300 м/с
33
Рисунок 49 – Профиль пятна и зависимость температуры от координаты для
давления 1 атм и скорости 300 м/с при диаметре пятна 40 мкм
Рисунок 50 – Зависимость выходной температуры от входной при
определенных диаметрах пятна для давления 1 атм и скорости 10 м/с
34
Рисунок 51 – Профиль пятна и зависимость температуры от координаты для
давления 1 атм и скорости 10 м/с при диаметре пятна 40 мкм
Рисунок 52 – Зависимость выходной температуры от входной при
определенных диаметрах пятна для давления 3 атм и скорости 300 м/с
35
Рисунок 53 – Профиль пятна и зависимость температуры от координаты для
давления 3 атм и скорости 300 м/с при диаметре пятна 40 мкм
На графиках видно, что существуют различные взаимодействия потока
с излучением. Таким образом, при малых температурах коэффициент
поглощения мал, поэтому поток не поглощает излучение и покидает область
фокуса с незначительным приростом температуры. Однако, при задании
выходной температуры, как например для рисунка (53), превышающей ~8000
K, для рисунка (51), ~5000 K, для рисунка (49), ~10000 K, видно, что в ходе
решения входная температура начинает уменьшаться. Это означает, что
существует некоторая максимальная входная температура, которая на данной
ветви решения не может быть превышена. Физический смыл этого можно
понять из того, что в диапазонах выходной температуры от ~18 до ~22 kK на
рисунке (53), от ~7 до ~25 kK на рисунке (51), от ~15 до ~25 kK на рисунке (49)
решений нет. При выходной температуре свыше ~22 kK на рисунке (53), ~25
kK на рисунке (51), ~25 kK на рисунке (49) возникает «горячее» решение, как
видно на рисунках (53,51,49) по кривым 5, данное решение отличается
прогревом среды еще до входа в лазерное пятно – за счет теплопередачи вверх
по потоку.
36
На рисунках (48,50,52) показаны зависимости выходной температуры
𝑇𝑜𝑢𝑡 от входной 𝑇𝑖𝑛 для диаметров пятна 𝑑 от 10 до 320 мкм на рисунках
(48, 52) и от 10 до 10000 мкм на рисунке (50) при постоянной мощности
излучения 𝑃 =1 кВт. Видно, что «зажигание» полностью ионизированной
плазмы разряда имеет пороговый характер. При температурах 𝑇𝑖𝑛 выше
некоторого порога (Для каждого 𝑑 - своего) «холодное» решение перестает
существовать, и температура 𝑇𝑜𝑢𝑡
повышается скачком до значений,
соответствующих полной ионизации. Это поведение можно ассоциировать с
зажиганием разряда. При снижении входной температуры разряд устойчиво
горит, вплоть до 𝑇𝑖𝑛 порядка комнатных температур.
37
2.6
Коды процедур
В открытой среде разработки программного обеспечения – Lazarus, на
языке Object Pascal для компилятора Free Pascal (часто используется
сокращение FPC — Free Pascal Compiler, бесплатно распространяемый
компилятор языка программирования Pascal), были написаны коды программ:
Рисунок 54 – объявление всех аргументов в процедуре, начальные параметры
Рисунок 55 – код процедуры с использованием аппроксимационной модели
для коэффициента поглощения
38
Таким образом, был построен график (56), который показывает, что
аппроксимационные формулы достаточно хорошо описывают значения
коэффициента поглощения (для (56) при 1 атм и температурах от 2000 до
20000 К), а также значения, полученные при запуске кода, совпадают с
аппроксимацией.
Рисунок 56 – log10(α(p,T))
Рисунок 57 – Код процедуры для молекулярного веса и энтальпии с
использованием модели [12]
39
Рисунок 58 – Используя процедуру для вычисления молекулярного веса и
энтальпии вычисляем теплоемкость с помощью дифференцирования
энтальпии
Рисунок 59 – код процедуры для теплопроводности при давлении 1 и 3 атм
40
2.6.1 Расчет надежности программы
2.6.1.1
Метод Мусса
В качестве оценки надежности программ использовалась модель Мусса.
В этой модели надежность программного обеспечения на этапе эксплуатации
оценивается по результатам тестирования.
Т - суммарное время тестирования, n – число отказов, произошедших за
время тестирования.
Тогда по модели Мусса средняя наработка до отказа определяется по
формуле
𝑡 = 𝑡0
𝐶×𝑇
𝑛×𝑡
×𝑒 0
(2.8)
где t 0 – средняя наработка до отказа до начала тестирования, C –
коэффициент, учитывающий уплотнение тестового времени по сравнению с
временем реальной эксплуатации, C = 10.
𝑡0 =
1
𝑓 ×𝐾 ×𝑁
(2.9)
где N – первоначальное число ошибок, его можно оценить с помощью
другой модели на основе статистических данных, полученных при
тестировании; К – коэффициент проявления ошибок. Обычно это значение
изменяется от 1,5*10 -7 до 4*10 -7; f – средняя скорость выполнения одного
оператора программы, равная отношению средней скорости исполнения (A) к
числу операторов (B).
Надежность для периода эксплуатации 100 ч определяется по формуле
100
𝑡
𝑅(𝑡) = 𝑒 −
41
(2.10)
Первоначальное количество ошибок находится по модели Шумана
метода подбора из уравнения
(2.11)
Наименьшее различие значений правой и левой частей этого уравнения
достигается при N=3, следовательно это и есть первоначальное количество
ошибок в программе. [16]
2.6.1.2
Расчет по методу Мусса
В качестве начальных параметров было взято 3 временных промежутка:
Рисунок 61 – количество отказов на
Рисунок 60 – время работы на
соответственных промежутках
разных промежутках
времени
Рисунок 62 – принятые и полученные ранее величины
Тогда
42
Рисунок 63 - средняя наработка до отказа до тестирования и средняя
наработка до отказа
От сюда надежность для периода работы 100 ч равна
Рисунок 64 – Ответ
43
2.7
Выводы по пункту 2
Построена математическая
и
численная
модель
непрерывного
оптического разряда в высокоскоростном газовом потоке. На основании
таблиц [13] и модели [12] получены аппроксимационные формулы для:
коэффициента
поглощения
лазерного
излучения,
теплоемкости,
теплопроводности, молекулярного веса и энтальпии воздушной плазмы при
давлениях 1 и 3 атм и температурах до 20000 К.
Результаты моделирования для давлений 1, 3 атм и скоростей 10, 300 м/с
показывают, что для мощности непрерывного излучения порядка 1 кВт
существуют условия, при которых стационарный оптический разряд в
сверхзвуковом газовом потоке с числами Маха до 5 существует.
В действительности потери энергии из области разряда больше, чем в
описанной одномерной равновесной модели. Это может иметь как
количественные последствия, такие, как повышение порога, так и
качественные – полное исчезновение стационарного «горячего» решения.
Выше была проделана первоначальная оценка равновесности плазмы
(п. 2.3), она является довольно грубой и требует дальнейших уточнений, как
будет видно в (п. 3.2), более корректная оценка показывает, что плазму в
работе равновесной назвать нельзя. В связи с чем, следует построить
неравновесную модель, которая предполагает под собой наличие разницы
температур электронов и тяжелых частиц, этому будет посвящен (п. 3).
44
3
разряда в
Моделирование неравновесного непрерывного оптического
высокоскоростном газовом потоке в
двухтемпературном
приближении
Прошлые
результаты
показали
возможность
локализации
и
поддержания непрерывного лазерного разряда в высокоскоростном газовом
потоке. Данная равновесная модель не учитывает множество параметров,
которые есть в реальности. Низкотемпературная плазма к равновесию не
приходит, как будет видно дальше (п. 3.2), однако для первоначальных оценок
и первых этапов работы она отлично подходит.
В этом пункте работы будет рассмотрена неравновесная аргоновая
плазма, но в двухтемпературном приближении, то есть существуют две
температуры: тяжелых частиц и электронов, отдельно электроны и тяжелые
частицы находятся в термическом равновесии, за счет столкновений
электроны (горячие частицы, температуры до 100 кК) разогревают тяжелые
частицы.
В качестве основного газа предложен аргон, так как масса атома больше
массы атома других газов, входящих в состав воздуха, атомы аргона будут
медленнее разогреваться, если теория покажет возможность поддержания
лазерного разряда в аргоне, то с уверенностью можно сказать, что
поддерживать разряд в воздухе возможно.
Предложена теория для расчета траекторий движения электрона вокруг
иона для определения сечения взаимодействия между частицами, проведена
оценка равновесности, показывающая неравновесность плазмы, приведена
математическая
модель, основываясь на
которой можно
провести
моделирование и узнать возможность поддержания лазерного разряда в
высокоскоростном газовом (в данном случае аргоновом) потоке.
45
3.1
Кулоново
взаимодействие
одиночного
электрона
с
положительным ионом.
Рассмотрим траекторию движения свободного электрона вокруг иона в
полностью ионизированной плазме для определения сечения взаимодействия
между частицами.
Примем во внимание, что масса электрона намного меньше массы иона,
следовательно, в данной траекторной задаче ион будет представлять из себя
стационарный объект.
По второму закону Ньютона:
𝐹⃑ = 𝑚𝑎⃑
(3.1)
Проецируя на оси x и y, получаем:
𝐹 = 𝑚𝑎𝑥
{ 𝑥
𝐹𝑦 = 𝑚𝑎𝑦
(3.2)
𝐹 = |𝐹|𝑐𝑜𝑠𝛽
{ 𝑥
𝐹𝑦 = −|𝐹|𝑠𝑖𝑛𝛽
(3.3)
где |𝐹| – кулоновское взаимодействие между частицами [17]
𝑛𝑒 2
𝑛𝑒2
|𝐹| =
=
4𝜋𝜀0𝑟 2 4𝜋𝜀0((𝑥 − 𝑥𝑓 )2 + 𝑦 2 )
(3.4)
2
𝑟 = (𝑥 − 𝑥𝑓 ) + 𝑦 2 – расстояние между электроном и ионом (x и y –
координаты электрона, 𝑥𝑓 -фокус гиперболы, по траектории которой движется
электрон, в фокусе расположен ион); e – заряд электрона, 𝜀0 – диэлектрическая
постоянная.
46
Тогда
𝑥𝑓 − 𝑥
𝑐𝑜𝑠𝛽 =
2
√((𝑥 − 𝑥𝑓 ) + 𝑦 2 )
𝑦
𝑠𝑖𝑛𝛽 =
(3.5)
2
√((𝑥 − 𝑥𝑓 ) + 𝑦 2 )
{
Используя, систему уравнений (3.5) и выражение (3.4), подставим их в
(3.3),
𝐹𝑥 =
𝐹𝑦 =
{
𝑛𝑒2
2
4𝜋𝜀0 ((𝑥 − 𝑥𝑓 ) + 𝑦 2)
−𝑛𝑒 2
2
4𝜋𝜀0((𝑥 − 𝑥𝑓 ) + 𝑦 2 )
𝑥𝑓 − 𝑥
×
2
√((𝑥 − 𝑥𝑓 ) + 𝑦 2)
𝑦
×
√((𝑥 − 𝑥𝑓 )
2
(3.6)
+ 𝑦2)
Таким образом, получены выражения для x и y компонент силы,
действующей на электрон.
𝐹𝑥 =
𝑛𝑒 2 × (𝑥𝑓 − 𝑥)
2
3
4𝜋𝜀0((𝑥 − 𝑥𝑓 ) + 𝑦 2 )2
−𝑛𝑒2 × 𝑦
𝐹𝑦 =
3
2
2 )2
4𝜋𝜀
((𝑥
−
𝑥
)
+
𝑦
{
0
𝑓
(3.7)
Из (3.2) выведем ускорения:
𝐹𝑥
𝑚
{
𝐹𝑦
𝑎𝑦 =
𝑚
𝑎𝑥 =
(3.8)
47
Из этого следует, что
𝑑2 𝑥
𝑛𝑒2 × (𝑥𝑓 − 𝑥)
=
𝑑𝑡2 4𝜋𝜀 𝑚((𝑥 − 𝑥 )2 + 𝑦 2 )32
0
𝑓
2
2
𝑑 𝑦
−𝑛𝑒 × 𝑦
=
𝑑𝑡2 4𝜋𝜀 𝑚((𝑥 − 𝑥 )2 + 𝑦 2 )32
{
0
𝑓
(3.9)
Так, рассмотрим задачу, при которой электрон подлетает к иону с
постоянным прицельным параметром 𝑃 = 1 м [18]; Прицельный параметр
представляет из себя расстояние между вектором начального направления
электрона (при 𝑉𝑥 ∞, то есть на дальних расстояниях от иона, при которых
электрон не ощущает на себе воздействия) и осью, параллельной направлению
начального вектора электрона и проходящей через фокус гиперболы.
Для дальнейшей работы с полученными выражениями для 2-ых
производных координат х и у их следует привести в безразмерный вид, то есть
1
поделить на прицельный параметр P (× ).
м
𝑥𝑓 𝑥
𝑛𝑒2 × ( 𝑃 − 𝑃) × 𝑃
𝑃×
=
𝑑𝑡2
𝑥 𝑥 2 𝑦2 3
4𝜋𝜀0𝑚(( − 𝑓 ) + )2 × 𝑃3
𝑃 𝑃
𝑃
𝑦
𝑦
𝑑 2 (𝑃)
−𝑛𝑒2 × 𝑃
𝑃×
=
𝑑𝑡2
𝑥 𝑥𝑓 2 𝑦 2 32
4𝜋𝜀
𝑚((
− ) + ) × 𝑃3
0
{
𝑃 𝑃
𝑃
𝑥
𝑑 2 (𝑃)
Так, получаем размерность
𝑛𝑒2
4𝜋𝜀 0 𝑚𝑃3
1
= 2, обозначим ее как:
с
𝑛𝑒2
1
=
4𝜋𝜀0𝑚𝑃3 𝜏 2
и обозначим все безразмерные величины, как:
48
(3.10)
(3.11)
𝑥
𝑥𝑓
𝑦
= 𝑥 ; = 𝑥𝑓 ; = 𝑦;
𝑃
𝑃
𝑃
𝑑2 𝑥 1
(𝑥𝑓 − 𝑥)
=
×
𝑑𝑡2 𝜏 2 ((𝑥 − 𝑥 )2 + 𝑦2 )32
𝑓
2
(3.12)
𝑑 𝑦
1
𝑦
=
−
×
𝑑𝑡2
𝜏 2 ((𝑥 − 𝑥 )2 + 𝑦2 )32
{
𝑓
Переобозначим и приведем в безразмерный вид, поделив на
1
𝜏2
𝑥 = 𝑥 ; 𝑥𝑓 = 𝑥𝑓 ; 𝑦 = 𝑦;
Таким образом получены безразмерные выражения для ускорения (𝑥̈ и
𝑦̈ ), зависящие от координат.
𝑥̈ =
(𝑥𝑓 − 𝑥)
2
3
2
+ 𝑦 )2
((𝑥 − 𝑥𝑓 )
𝑦
𝑦̈ = −
3
2
2
((𝑥 − 𝑥𝑓 ) + 𝑦 )2
{
(3.13)
Рассмотрим уравнение гиперболы [19]
𝑥2 𝑦 2
−
=1
𝑎2 𝑏2
(3.14)
На больших расстояниях от фокуса гиперболы (в ∞) получаем:
𝑥2 𝑦 2
≈
𝑎2 𝑏2
(3.15)
49
Тогда тангенс угла между осью х и асимптотой гиперболы равен:
𝑦
𝑏
sin 𝛼
𝑏 2
sin2 𝛼
= ∓ = tan 𝛼 =
;( ) =
𝑥
𝑎
1 − sin2 𝛼
√1 − sin2 𝛼 𝑎
(3.16)
Найдем синус этого угла
𝑏
2
( )2
𝑏
𝑏2
𝑎
2
2
2
sin 𝛼 = (1 − sin 𝛼 ) ( ) ; sin 𝛼 =
=
𝑏 2 𝑏2 + 𝑎2
𝑎
1+( )
𝑎
(3.17)
Расстояние от начала координат до фокуса гиперболы при b=1:
1
√𝑎2 + 𝑏2
𝑥𝑓 =
=
= √𝑎2 + 1
sin 𝛼
𝑏
(3.18)
Теперь возьмем производную из уравнения гиперболы для того, чтобы
узнать выражения для скоростей на всем диапазоне гиперболы:
2
𝑥𝑥̇
𝑦𝑦̇
−
2
=0
𝑎2
𝑏2
𝑦𝑎2
𝑉𝑥 = 𝑥̇ = 2 𝑦̇
𝑥𝑏
(3.19)
(3.20)
Берем вторую производную и выражаем 𝑦̇ = 𝑉𝑦 . 𝑦̈ и 𝑥̈ подставляем из
полученных ранее выражений.
(𝑥̇ )2 + 𝑥̈ × 𝑥
𝑎2
(𝑦̇ )2 + 𝑦̈ × 𝑦
−
=0
𝑏2
50
(3.21)
Получаем
𝑥̈ × 𝑥𝑏2 − 𝑦̈ × 𝑦𝑎2
𝑉𝑦 = 𝑦̇ = √
𝑦 2 𝑎4
𝑎2 − 2 2
𝑥 𝑏
(3.22)
Тогда, подставив все начальные условия 𝑥0 = 𝑎; 𝑏 = 𝑃 = 1; 𝑦0 = 0; 𝑉𝑥 =
0; данное выражение упроститься:
𝑉𝑦 2 =
𝑥̈
1
=
𝑎 𝑎 × (𝑥𝑓 − 𝑎)2
(3.23)
Таким образом, мы имеем систему
𝑥̈ =
(𝑥𝑓 − 𝑥)
3
2
((𝑥 − 𝑥𝑓 ) + 𝑦 2 )2
𝑦
𝑦̈ = −
3
2
((𝑥 − 𝑥𝑓 ) + 𝑦 2)2
{
(3.24)
С начальными условиями для интегрирования:
𝑥0 = 𝑎; 𝑏 = 𝑃 = 1; 𝑦0 = 0; 𝑉𝑥 = 0; 𝑉𝑦 =
1
;
√𝑎 × (𝑥𝑓 − 𝑎)
Зная 𝑉𝑦 , получаем, что момент
𝑀 = 𝑉∞ 𝑃 = 𝑉𝑦 (𝑥𝑓 − 𝑎) =
51
1
𝑃√ 𝑎
(3.25)
Следовательно, 𝑉∞ , при P=1.
𝑉∞ =
1
√𝑎
(3.26)
Приведем все обратно в размерный вид:
𝑉̅∞ =
1
𝑎
𝑃
; 𝑎̅ = ; 𝑉∗ = ;
𝑃
𝜏
√𝑎̅
𝑉̅∞ =
𝑉∞ 𝑉∞ 𝜏
=
𝑉∗
𝑃
(3.27)
Тогда
𝑉∞ 𝜏
𝑃 𝑃
𝑉∞ 𝜏 2
=√ ; =(
) ;
𝑃
𝑎 𝑎
𝑃
𝑃3
𝑎 = 2 2;
𝑉∞ 𝜏
Подставим
(3.28)
1
𝜏2
𝑛𝑒2
𝑎=
4𝜋𝜀0 𝑚𝑉∞ 2
(3.29)
Исходя из полученных выражений можно варьировать, например
скорость подлета электрона и интегрировать по прицельному параметру.
Верхняя граница интегрирования будет равняться Дебаевскому радиусу [20],
таким образом получаться все возможные гиперболические траектории
движения электрона вокруг иона.
52
Ниже представлены графики траектории движения электрона вокруг
неподвижного иона с различными скоростями подлета электрона (рис. 65-68),
а также разными прицельными параметрами при фиксированной скорости
подлета (рис. 69-72).
Рисунок 65 - Траектория движения
Рисунок 66 - Траектория движения
электрона вокруг иона при
электрона вокруг иона при
параметрах V=5, P=1
параметрах V=10, P=1
Рисунок 67 - Траектория движения
Рисунок 68 - Траектория движения
электрона вокруг иона при
электрона вокруг иона при
параметрах V=25, P=1
параметрах V=50, P=1
53
Рисунок 69 - Траектория движения
Рисунок 70 - Траектория движения
электрона вокруг иона при
электрона вокруг иона при
параметрах V=50, P=0,5
параметрах V=50, P=1
Рисунок 71 - Траектория движения
Рисунок 72 - Траектория движения
электрона вокруг иона при
электрона вокруг иона при
параметрах V=50, P=1,5
параметрах V=50, P=2
Теперь стоит рассмотреть импульс, который электрон придает иону,
после прохождения его по гиперболической траектории.
Скорость электрона до взаимодействия с ионом равна скорости после.
𝑉∞ + = 𝑉∞ −
(3.30)
54
Рисунок 73 – представление треугольника скоростей.
Рассмотрим равнобедренный треугольник, образованный данными
скоростями, основание данного треугольника – скорость, приобретенная
ионом (∆𝑉∞ ), формулу которой следует вывести, угол между 𝑉∞ − и 𝑉∞ + будет
равен
𝜃 = (𝜋 − 2𝛼)
(3.31)
где 𝛼 – половина угла между асимптотами гиперболы.
Тогда проведя биссектрису угла 𝜃 к основанию (∆𝑉∞ ), получим, что 𝑉∞ +
является гипотенузой данного прямоугольного треугольника, от сюда
получаем, что
𝜃 𝜋
= −𝛼
2 2
(3.32)
Углы при основании треугольника будут равны 𝛼, так как
55
90 −
𝜋
−𝛼 =𝛼
2
(3.33)
Далее рассмотрим треугольник, образованный перпендикуляром из угла
основания к стороне начальной скорости (𝑉∞ −). Гипотенузой данного
треугольника является ∆𝑉∞ , катеты: ∆𝑉∞ перп. и ∆𝑉∞ парал., соответственно,
перпендикулярная и параллельная составляющие ∆𝑉∞ .
∆𝑉∞ = √∆𝑉∞ перп . 2 + ∆𝑉∞ парал. 2
(3.34)
Рассмотрим ∆𝑉∞ перп. и ∆𝑉∞ парал .
∆𝑉∞ перп . = 𝑉∞ + × sin 𝜃 = 𝑉∞ + × sin 2𝛼 = 𝑉∞ + × 2 sin 𝛼 cos 𝛼
(3.35)
∆𝑉∞ парал. = 𝑐𝑡𝑔 𝛼 × ∆𝑉∞ перп . = 𝑉∞ + × 2 cos 2 𝛼
(3.36)
Теперь подставим полученные выражения в (3.34)
∆𝑉∞ = √(𝑉∞ + × 2 sin 𝛼 cos 𝛼)2 + (𝑉∞ + × 2 cos 2 𝛼)2
(3.37)
Получаем выражение для скорости, переданной электроном иону в
процессе взаимодействия
∆𝑉∞ = 2 𝑉∞ + × cos 𝛼
где 𝑉∞ + = 𝑉∞ , которую можно выразить из (3.29).
56
(3.38)
Оценка равновесности
3.2
Очень важно оценить равновесность изучаемой плазмы, так как от этого
зависит
дальнейший путь
всех
расчетов и
построений моделей.
Термодинамика изучает системы, которые находятся в состоянии теплового
равновесия, а температура является важнейшей характеристикой данной
системы. Понятие температуры имеет смысл только при наличии хотя бы
частичного равновесия.
Плазма, находящаяся в состоянии полного термодинамического
равновесия, имеет одну определенную температуру. Однако, очень часто
приходится иметь дело с частичным термодинамическим равновесием. Так,
из-за огромной разницы в массе между электронами и ионами, обмен энергией
электронов с ионами происходит гораздо медленней, чем обмен между
частицами близкими по массе. Таким образом, плазма может длительное
время существовать состояние, при котором она характеризуется двумя
температурами: электронной и ионной, каждая из которых находиться в
термодинамическом равновесии. Плазму с равной температурой ионов и
электронов
довольно
сложно
получить,
такую
плазму
называют
изотермической.
Из написанного выше стоит проверить, можно ли изучаемую в данной
работе, плазму, назвать изотермической.
Перед описанием метода проверки стоит сделать некоторые допущения.
В данном методе начальная скорость иона считается равной нулю, то есть ион
до момента соударения с электроном находился в состоянии покоя, а процесс
удара электрона с ионом является абсолютно упругим, следовательно,
деформаций не происходит, и полная кинетическая энергия сохраняется.
Рассмотрим систему
{
̅𝑖
𝑀𝑉𝑖 + 𝑚𝑉𝑒 = 𝑚𝑉̅𝑒 + 𝑀𝑉
2
̅𝑖 2
𝑀𝑉𝑖 2 + 𝑚𝑉𝑒 2 = 𝑚𝑉̅𝑒 + 𝑀𝑉
57
(3.39)
где М – масса иона, m – масса электрона, 𝑉𝑖 – начальная скорость иона, 𝑉𝑒 –
̅𝑖 –
начальная скорость электрона, 𝑉̅𝑒 – скорость электрона после соударения, 𝑉
скорость иона после соударения.
Так как было допущено, что начальная скорость иона равна 0, тогда
система (3.39) примет вид
{
̅𝑖
𝑚𝑉𝑒 = 𝑚𝑉̅𝑒 + 𝑀𝑉
2
̅𝑖 2
𝑚𝑉𝑒 2 = 𝑚𝑉̅𝑒 + 𝑀𝑉
(3.40)
̅𝑖
Для каждого из уравнений в системе выводим 𝑉
𝑚
(𝑉𝑒 − 𝑉̅𝑒 )
𝑀
{
𝑚
̅𝑖2 = (𝑉𝑒 2 − 𝑉̅𝑒 2 )
𝑉
𝑀
̅𝑖 =
𝑉
(3.41)
Теперь приравняем данные уравнения и выведем результат для 𝑉̅𝑒
𝑚 2
𝑚
2
(𝑉𝑒 − 𝑉̅𝑒 ) = ( )2 (𝑉𝑒 − 𝑉̅𝑒 )2
𝑀
𝑀
𝑉̅𝑒 = 𝑉𝑒 (
𝑀−𝑚
)
𝑀+𝑚
(3.42)
(3.43)
̅𝑖
Результат (3.43) нужно подставить в (3.41) к 𝑉
̅𝑖 =
𝑉
𝑚
𝑀−𝑚
𝑚 2 𝑉
)) = 2 ( ) ( 𝑒 𝑚)
𝑉𝑒 (1 − (
𝑀
𝑀+𝑚
𝑀
1+
𝑀
58
(3.44)
Для оценки необходимо знать, достаточно ли таких электрон-ионных
столкновений, в области, где локализован сфокусированный в пучек d = 50
мкм лазерный луч, со скоростью газового потока U = 300 м/с.
Для этого нужно посчитать число Дамкеллера
𝐷𝑎 =
𝜏𝑈
𝑑
(3.45)
𝜏 – время электрон-ионных столкновений
𝜏=
𝐸𝑒
𝑓𝐸𝑖
(3.46)
𝐸𝑒 и 𝐸𝑖 – энергии электрона и иона соответственно, 𝑓 – частота столкновений
𝑓 = 𝜎𝑛𝑒 𝑉𝑒
(3.47)
где 𝜎 – сечение столкновения, которое равняется Дебаевской длине (𝑙𝐷 ), 𝑛𝑒 –
концентрация электронов (𝑛𝑒 =
𝑝
2𝑘𝑇
, где p – давление, k – постоянная
3𝑘𝑇
Больцмана, T – температура), 𝑉𝑒 - скорость электрона (𝑉𝑒 = √ , где 𝑚𝑒 - масса
𝑚
𝑒
электрона).
Дебаевская длина определяется по формуле
𝑞𝑖 2 𝑛𝑖 −1
𝑙𝐷 = (∑
) 2
𝜀0 𝑘𝑇𝑖
(3.48)
𝑖
Здесь индекс i обозначает элемент (элетрон/ион), q – заряд, n –
концентрация, 𝜀0- диэлектрическая постоянная.
Межмолекулярное расстояние равно
59
𝐿=
1
√𝑛
3
(3.49)
При высоких температурах Дебаевская длина приближается к значению
межмолекулярного расстояния.
Энергия иона после удара с подстановкой (3.44) равна
̅𝑖 2
𝑚𝑉
4
𝑚 3 𝑚𝑉𝑒 2
̅
𝐸𝑖 =
=
𝑚 ( ) 2
2
(1 + )2 𝑀
𝑀
(3.50)
Энергия электрона
𝑚𝑉𝑒 2
𝐸𝑒 =
2
(3.51)
Таким образом, подставив выражения (3.48) в (3.47), (3.51, 3.50, 3.47) в
(3.46) и (3.46) в (3.45) получим оценку равновесности.
Используя эти оценки, получилось, что плазму изотермической считать
нельзя. Однако, равновесие по отдельно взятым температурам должно быть,
следовательно в данной задаче необходимо использовать двух - (или более)
температурную модель.
60
3.3
Математическая модель непрерывного лазерного разряда в
высокоскоростном газовом потоке в двухтемпературном приближении
Для двухтемпературного моделирования аргоновой плазмы в первую
очередь нужно рассмотреть уравнение энергии тяжелых частиц, то есть Ar,
Ar *, Ar + (соответственно атом аргона, возбужденный аргон и ион аргона), а
также уравнение энергии электронов, и раскрыть поочередно каждый член.
3.3.1 Уравнение энергии тяжелых частиц [21]:
5
𝛻 ( 𝑘𝐵 𝑇ℎ 𝑛ℎ 𝜐 ) + 𝛻 (−𝜆ℎ 𝛻𝑇ℎ + ∑ ℎ𝑖 𝐽𝑖 ) = 𝑄𝑒𝑙
2
(3.52)
𝑖
где 𝑛ℎ - концентрация (всех частиц) соответствующего элемента h (Ar, Ar*,
Ar+), 𝜐 – конвективная скорость газа, 𝑘𝐵 – постоянная Больцмана, 𝑇ℎ температура тяжелых частиц, 𝜆ℎ - коэффициент теплопроводности тяжелых
5
частиц, ℎ𝑖 = 𝑘𝐵 𝑇ℎ - удельная энтальпия тяжелых частиц, 𝐽𝑖 - массовый
2
диффузионный поток тяжелых частиц, 𝑄𝑒𝑙 - передача энергии электронов за
счет упругих столкновений.
Коэффициент теплопроводности газовой смеси [22]
𝜆ℎ = ∑ (
𝑖
𝑛𝑖
) 𝜆∗𝑖
∑𝑟 𝑛𝑟 𝑀𝑖𝑟
(3.53)
где
𝑚𝑖𝑟 𝜎̅𝑖𝑟
𝑀𝑖𝑟 = √2
𝑚𝑟 𝜎̅𝑖𝑖
(3.54)
𝜆∗𝑖 = 𝑘𝐵 𝑛𝑖 𝜐̅ 𝑖𝑙𝑖
(3.55)
𝜐̅𝑖 - средняя относительная тепловая скорость тяжелых частиц, равная
61
(3.56)
8𝑘𝐵 𝑇ℎ
𝜐̅𝑖 = √
𝜋𝑚𝑖
а 𝑙𝑖 - средняя длина свободного пробега, определяемая выражением
−1
1
𝑚𝑖𝑟 2
𝑙𝑖 = (∑2 (
𝑟
𝑚𝑖
(3.57)
) 𝑛𝑟 𝜎̅𝑖𝑟 )
𝑚𝑖𝑟 – приведенная масса:
𝑚𝑖𝑟 =
𝑚𝑖 𝑚𝑟
𝑚𝑖 + 𝑚𝑟
(3.58)
𝜎̅𝑖𝑟 – среднее сечение упругого столкновения, взятые из баз данных [23] для
Ar-Ar + и [24] для Ar-Ar
Рисунок 74 – Сечения упругих реакций с Ar в зависимости от температуры.
62
Выше приведены сечения упругих столкновений Ar – Ar и Ar – Ar+,
сечение столкновений Ar + – Ar + в данной работе принято равным сечению
упругих столкновений Ar – Ar+, однако этот вопрос в дальнейшем нужно будет
изучить.
Массовый диффузионный поток тяжелых частиц [25]
Диффузионный поток с учетом амбиполярной диффузии определяется
следующим образом:
𝐽ℎ = −
𝑝𝜇 ℎ 𝐷ℎ
𝑝𝜇𝑟 𝐷𝑟
𝐻ℎ + 𝑦ℎ ∑
𝐻
𝑅0 𝑇ℎ
𝑅0 𝑇𝑟 𝑟
(3.59 a)
𝑟
+
1
(𝜇 𝑞 𝜌 𝐷 − 𝑦ℎ ∑ 𝜇𝑟 𝑞𝑟 𝜌𝑟 𝐷𝑟 ) 𝐸𝑎
𝑅0 𝑇ℎ ℎ ℎ ℎ ℎ
𝑟
где p - давление газа, 𝜇ℎ , 𝑇ℎ и 𝐷ℎ - соответственно молекулярная масса,
температура и эффективный коэффициент диффузии частиц h, 𝑦ℎ - массовая
доля, 𝐻ℎ - движущая сила, 𝑞ℎ - заряд, 𝐸𝑎 – амбиполярное электрическое поле,
а 𝑅0 – универсальная газовая постоянная.
Так как в нашей работе плазма является квазинейтральной, тогда 𝐸𝑎 =
0, уравнение в таком случае принимает простой вид:
𝐽ℎ = −
𝑝𝜇 ℎ 𝐷ℎ
𝑝𝜇𝑟 𝐷𝑟
𝐻ℎ + 𝑦ℎ ∑
𝐻
𝑅0 𝑇ℎ
𝑅0 𝑇𝑟 𝑟
(3.59 б)
𝑟
Эффективный коэффициент диффузии 𝐷ℎ
частиц h определяется
следующим образом:
𝐷ℎ =
1 − 𝑥ℎ
𝑥
∑𝑟 𝑟
𝐷ℎ𝑟
(3.60)
63
где 𝑥ℎ - мольная доля частиц h, а 𝐷ℎ𝑟 - коэффициент массовой диффузии
бинарной смеси определяется как:
𝐷ℎ𝑟
(
2𝑘𝐵 𝑇(𝑚ℎ + 𝑚𝑟 ) 0.5
)
=
(1,1) (
𝜋𝑚
𝑚
16𝑛ℎ 𝛺ℎ𝑟
ℎ 𝑟
3
(3.61)
)
где 𝛺ℎ𝑟1,1 - температурный интеграл столкновений из [26]
Рисунок 75 – Интегралы столкновений
(
)
1,1
Значение в источнике для Ar – Ar 𝛺𝐴𝑟𝐴𝑟
= const =1,26E-15 м2, однако есть
сомнения в этих данных, скорее всего значения приведены в см2, тогда график
будет выглядеть правдоподобнее.
Движущая сила определяется выражением
1
𝐻ℎ = 𝛻𝑧ℎ + (𝑧ℎ − 𝑦ℎ )𝛻𝑙𝑛𝑝 − ∑(𝛽ℎ𝑟 𝛻𝑙𝑛 𝑇ℎ − 𝛽𝑟ℎ 𝛻𝑙𝑛𝑇𝑟 )
𝑝
𝑟
64
(3.62)
где 𝑧ℎ =
𝑝ℎ
𝑝
, 𝑝ℎ - парциальное давление частиц h, 𝛽ℎ𝑟 - коэффициент
термодиффузии. между видами h и r [27].
Коэффициент термодиффузии
𝛽ℎ𝑟 =
𝑘𝐵 𝑇𝑝ℎ 𝑝𝑟 𝜏𝑟 𝑚ℎ𝑟
2𝑝𝑚𝑟 2 𝐷ℎ𝑟
(3.63)
τr – обратная частота столкновений, Dhr – по формуле (3.61)
1 1
𝑛𝑖 𝜎ℎ𝑖
𝜏𝑟 =
= (∑
)
𝜈𝑟 2
√𝜋𝛾ℎ𝑖
−1
(3.64)
𝑖
где 𝛾ℎ𝑖 =
𝑚ℎ𝑖
2𝑘𝐵 𝑇ℎ𝑖
, а 𝑇ℎ𝑖 =
(𝑚ℎ 𝑇𝑖 +𝑚𝑖 𝑇ℎ )
(𝑚ℎ +𝑚𝑖 )
Передача энергии электронов за счет упругих столкновений равна [21]:
3𝑚𝑒
𝑄𝑒𝑙 =
𝑛 𝑘 (𝑇 − 𝑇ℎ ) (∑ 𝜈𝑒ℎ )
𝑚ℎ 𝑒 𝐵 𝑒
(3.66)
ℎ
где νeh частота упругих столкновений из [22]
∑ 𝜈𝑒ℎ = ∑ 𝑛ℎ 𝜐̅𝑒 𝜎̅𝑒ℎ
ℎ
ℎ
65
(3.67)
3.3.2 Уравнение энергии электронов [21]:
5
𝛻 ( 𝑛𝑒 𝑘𝐵 𝑇𝑒 𝜐) + 𝛻(−𝜆𝑒 𝛻𝑇𝑒 + ℎ𝑒 𝐽𝑒 ) = 𝑗𝐸 − 𝑄𝑒𝑙 + 𝑄𝑖𝑛 − 𝑄𝑟𝑎𝑑
2
(3.68)
где 𝑛𝑒 – концентрация электронов, 𝑘𝐵 – постоянная Больцмана, 𝑇𝑒 –
температура электронов, 𝜐 – конвективная скорость газа, λ𝑒 – электронная
5
теплопроводность, ℎ𝑒 = 𝑘𝐵 𝑇𝑒
2
– энтальпия электронов, 𝐽𝑒 – массовый
диффузионный поток электронов, 𝑄𝑒𝑙 – передача энергии за счет упругих
столкновений, 𝑄𝑖𝑛 – передача энергии за счет неупругих столкновений, 𝑗𝐸 –
член Джоулевого нагрева, 𝑄𝑟𝑎𝑑 – радиационные потери.
В работе [21] рассматривается плазма, созданная плазменной горелкой
поэтому здесь присутствует член Джоулевого нагрева, в нашей работе он
будет заменен на световой член. Также 𝑄𝑟𝑎𝑑 не учитывается в работе, этот
вопрос нужно подробнее изучить, однако ниже будет некоторая информация
по радиационным потерям, которой можно будет пользоваться в дальнейшем.
Электронная теплопроводность [22]
2,4
𝑘𝐵 2 𝑛𝑒 𝑇𝑒
𝜆𝑒 =
𝜈̅
1 + 𝑒𝑖 𝑚𝑒 𝜈̅𝑒𝐻
√2𝜈̅𝑒𝐻
(3.69)
где 𝑚𝑒 – масса электрона, 𝜈̅ 𝑒𝑖 – частота столкновений электронов с ионами
2
4√2𝜋 𝑚𝑒 3
𝑒2
(
)2 (
𝜈̅ 𝑒𝑖 = 𝑛𝑖
) ln(𝛬)
3
𝑘𝐵 𝑇𝑒
4𝜋𝜀0 𝑚𝑒
(3.70)
про ln(𝛬) будет сказано ниже, 𝜈̅𝑒𝐻 – полная средняя частота столкновений
электронов со всеми тяжелыми частицами определяется аналогично формуле
(3.67),
66
Ниже представлен рисунок (76), 𝜎𝑒ℎ - сечение упругих столкновений
электрона с атомом и возбужденным атомом аргона, данные были взяты из
[23]
Рисунок 76 – сечение упругих столкновений e с Ar и Ar *
Сечение электрон электронных и электрон ионных столкновений
определяется по формуле (3.71) [22]
̅̅̅̅2 ln (Λ)
̅̅̅̅
𝑄𝑒𝑖 = 6𝜋𝑏
0
(3.71)
где
Λ=
𝜆𝐷
𝑏0
(3.72)
𝜆𝐷 – Дебаевская длина, которая равна
67
(3.73)
𝜀0 𝑘𝐵 𝑇
𝜆𝐷 = √
𝑛𝑒 𝑒 2
𝜀0 – диэлектрическая постоянная, а 𝑏0 :
𝑍𝑒 2
𝑏0 =
12𝜋𝜀0 𝑘𝐵 𝑇
(3.74)
здесь e и Ze – заряды двух частиц. Величину 𝑏0 в соответствии с классической
механикой можно интерпретировать как прицельное расстояние для рассеяния
на угол, равный 90 (прицельное расстояние – расстояние между частицей
мишенью и первоначальной траекторией падающей частицы)
Диффузионный
поток
электронов
определяется
аналогично
диффузионному потоку тяжелых частиц (3.59), так же все формулы (3.59-3.64)
используются в расчете для электронов, но с некоторыми изменениями, ниже
будут формулы для диффузионного потока электронов.
𝐽𝑒 = −
𝑝𝜇𝑒 𝐷𝑒
𝑝𝜇 𝑟 𝐷𝑟
𝐻𝑒 + 𝑦ℎ ∑
𝐻
𝑅0 𝑇𝑒
𝑅0 𝑇𝑟 𝑟
(3.75 a)
𝑟
Так как массовая доля 𝑦𝑒 для электронов пренебрежимо мала, данный
член суммы не будет вносить существенных изменений в диффузионный
поток в целом. Тогда уравнение примет вид:
𝐽𝑒 = −
𝑝𝜇𝑒 𝐷𝑒
𝐻
𝑅0 𝑇𝑒 𝑒
Эффективный коэффициент диффузии 𝐷𝑒
следующим образом:
68
(3.75 б)
электронов определяется
𝐷𝑒 =
1 − 𝑥𝑒
𝑥
∑𝑟 𝑟
𝐷𝑒𝑟
(3.76)
где 𝑥𝑒 , 𝑥𝑟 - мольная доля частиц 𝑒 и r, а 𝐷𝑒𝑟 - коэффициент массовой диффузии
бинарной смеси определяется как:
𝐷ℎ𝑟
(
2𝑘𝐵 𝑇𝑒 (𝑚𝑒 + 𝑚𝑟 ) 0.5
)
=
( )(
𝜋𝑚𝑒 𝑚𝑟
16𝑛𝑒 𝛺𝑒𝑟1,1
3
(3.77)
)
𝛺𝑒𝑟1,1 - температурный интеграл столкновений из [26]
Рисунок 77 – Интегралы столкновений
Интеграл столкновений для e-Ar + равен интегралу столкновений e-e,
также как и сечения, а интеграл столкновений e-Ar * равен интегралу
столкновений e-Ar
Движущая сила определяется выражением
1
𝐻𝑒 = 𝛻𝑧𝑒 + 𝑧𝑒 𝛻𝑙𝑛𝑝 − ∑(𝛽𝑒𝑟 𝛻𝑙𝑛 𝑇𝑒 − 𝛽𝑟𝑒 𝛻𝑙𝑛𝑇𝑟 )
𝑝
𝑟
69
(3.78)
где 𝑧𝑒 =
𝑝𝑒
𝑝
, 𝑝𝑒 - парциальное давление электронов, 𝛽𝑒𝑟 - коэффициент
термодиффузии. между e и r [27].
Коэффициент термодиффузии
𝛽𝑒𝑟 =
𝑘𝐵 𝑇𝑝𝑒 𝑝𝑟 𝜏𝑟 𝑚𝑒𝑟
2𝑝𝑚𝑟 2 𝐷𝑒𝑟
(3.79)
τr – обратная частота столкновений, Der – по формуле (3.61)
1 1
𝑛𝑖 𝜎𝑒𝑖
𝜏𝑟 =
= (∑
)
𝜈𝑟 2
√𝜋𝛾𝑒𝑖
−1
(3.80)
𝑖
где 𝑖 – электроны и тяжелые частицы 𝛾𝑒𝑖 =
𝑚𝑒𝑖
2𝑘𝐵 𝑇𝑒𝑖
, а 𝑇𝑒𝑖 =
(𝑚𝑒 𝑇𝑖 +𝑚𝑖 𝑇𝑒 )
(𝑚𝑒 +𝑚𝑖 )
, сечения
определяются по формуле (3.70) и рис. (76)
Во всех формулах индекс r может быть, как тяжелыми частицами, так
как диффузия электронов зависит и от тяжелых частиц, так и электронами.
Передача энергии электронов за счет упругих столкновений равна
формуле (3.66) [21], а передача энергии электронов за счет неупругих
столкновений равна [21]:
𝑄𝑖𝑛 = − ∑ ∆𝜀𝑖 𝑘𝑖 𝑛𝑒 𝑛𝑖
(3.81)
𝑖
где ∆𝜀𝑖 - потеря энергии реакции, 𝑘𝑖 - коэффициенты скорости реакции [21],
𝑛𝑒 , 𝑛𝑖 – концентрация электронов и тяжелых частиц.
Радиационные потери из [32]
𝑄𝑟𝑎𝑑 = 𝑛𝑒 𝑛ℎ 𝐾𝑖 𝜀𝑖
70
(3.82)
где 𝜀𝑖 – энергия в электронвольтах, а 𝐾𝑖 - коэффициент выбросов
𝐾𝑖 = ∫ 𝑓(𝜀 ) 𝜎𝑖 (𝜀)𝜐(𝜀)𝑑𝜀
(3.83)
𝜎𝑖 (𝜀) – эмиссионное сечение, 𝜐(𝜀) – скорость электронов, 𝑓(𝜀 ) – функция,
описывающая электронную энергию
𝑓 (𝜀 ) =
3
0,565𝑇𝑒 2
1
𝜀2
𝜀2
𝑒𝑥𝑝(−0.243 2 )
𝑇𝑒
71
(3.84)
3.4
Плазмохимические процессы
Рассмотрим систему плазмохимических реакций для аргона, ниже
приведена таблица (4).
Реакции
1
Коэффициенты скорости реакции
𝑒 + 𝐴𝑟
4.9 ∗ 10
→ 𝑒 + 𝐴𝑟(4𝑠)
2
𝑒 + 𝐴𝑟 → 2𝑒 + 𝐴𝑟 +
3
𝑒 + 𝐴𝑟(4𝑠) → 2𝑒
+ 𝐴𝑟
4
+
𝑒 + 𝐴𝑟(4𝑠) → 𝑒
+ 𝐴𝑟
5
2𝑒 + 𝐴𝑟 + → 𝑒 + 𝐴𝑟
−15 (
1.27 ∗ 10
𝑇𝑒 [𝑒𝑉
])0.5
−11.65 𝑚3
)[ ]
exp (
𝑇𝑒 [𝑒𝑉 ]
𝑠
−14 (
𝑇𝑒 [𝑒𝑉 ])0.5 exp (
−13 (
1.37 ∗ 10
𝑇𝑒 [𝑒𝑉
4.8 ∗ 10
])0.5
−15.76 𝑚3
)[ ]
𝑇𝑒 [𝑒𝑉 ]
𝑠
−4.11 𝑚3
)[ ]
exp (
𝑇𝑒 [𝑒𝑉 ]
𝑠
𝑚3
0.5 [
[
])
]
𝑇𝑒 𝑒𝑉
𝑠
−16 (
4.5
1
𝑚6
) [ ] для 𝑇𝑒
𝑇𝑒 [𝑒𝑉 ]
𝑠
−39 (
8.75 ∗ 10
≤ 0.276 𝑒𝑉
11.659
1.29 ∗ 10−44 (
+ 2) ×
𝑇𝑒 [𝑒𝑉 ]
4.11 𝑚6
) [ ] для 𝑇𝑒
× exp (
𝑇𝑒 [𝑒𝑉 ] 𝑠
> 0.276 𝑒𝑉
Таблица 4 – химические реакции и их коэффициенты скорости из [21]
Температура приведена в электронвольтах.
Электронвольт – это энергия, приобретаемая частицей, заряд которой
равен заряду электрона e, после ускорения в поле с разностью потенциалов в
1 В, переведем в Кельвин
эВ – e*U
72
3
𝑒𝑈 = 𝑘𝐵 𝑇
2
(3.85)
где 𝑇 – температура, 𝑘𝐵 – постоянная Больцмана
𝑇=
2
2𝑒
𝑒𝑈 | =≫ |𝑈 = 1В| =≫ 𝑇 =
3𝑘𝐵
3𝑘𝐵
(3.86)
От сюда следует, что 1 эВ = 7 729,47 К
Размерность коэффициентов скоростей реакций – м3/с
𝑘 = 𝜎𝜐
(3.87)
где 𝜎 – сечение конкретной реакции, 𝜐 – средняя тепловая скорость
электрона
Рассмотрим реакцию возбуждения аргона:
𝑒 + 𝐴𝑟 → 𝑒 + 𝐴𝑟(4𝑠)
(3.88)
𝑑[𝐴𝑟]
→ −𝑘[𝑒][𝐴𝑟(4𝑠)]
𝑑𝑡
𝑥 м3
1
=
∗ 𝑥 ∗ 𝑥 | =≫ 𝑥 = 3
𝑐
с
м
От сюда видно, что размерность в «штуках» на кубический метр,
переведем в моль на кубический метр, для этого умножим и разделим на число
Авогадро
𝑁𝐴 𝑑[𝐴𝑟]
𝑁𝐴
𝑁𝐴
→ −𝑘[𝑒] [𝐴𝑟(4𝑠)]
| =≫
𝑁𝐴 𝑑𝑡
𝑁𝐴
𝑁𝐴
73
(3.89)
𝑑 [𝐴𝑟]𝑚
м3 1
𝑚 [ ]𝑚 [
𝑚 [
(
)]
]
| =≫
→ −𝑘 𝑒 𝐴𝑟 4𝑠
𝑑𝑡
с моль
где
м3
8𝑘 𝑇
𝑚
] = 𝜎𝜐̅𝑒 𝑁𝐴 = 𝜎𝑁𝐴 √ 𝐵 𝑒
𝑘 = 𝑘𝑁𝐴 [
с ∗ моль
𝜋𝑚𝑒
74
(3.90)
3.5
Выводы по пункту 3
Таким
образом,
была
предложена
двухтемпературного лазерного разряда в
математическая
аргоне.
модель
Основываясь
на
предложенной методике построения математической модели описывается
процесс поддержания неравновесного лазерного разряда в высокоскоростном
газовом потоке. По данной модели, учитывая уравнение неразрывности,
уравнение сохранения массы и все физические величины, могут быть
проведены построения одномерной и трехмерной модели процесса.
Обязательно нужно учесть, что полученные результаты не являются
конечными, модель требует дальнейших оценок и уточнений.
75
4
Методика экспериментального исследования
Для будущих исследований обязательно нужно предложить методику
экспериментального исследования. Основываясь на оценке параметров
лазерного излучения, необходимого для поддержания стационарного разряда
были предложены и посчитаны оптические схемы, разработана выхлопная
труба с оптическими окнами для регистрации и подвода излучения,
перечислены и подобраны все комплектующие, наличие которых обязательно
для проведения экспериментов
76
4.1
Оптическая схема
Предложена полная оптическая схема, для реализации и фиксирования
эксперимента представлена на рисунке (78). Чертеж данной схемы можно
посмотреть в приложении А.
Рисунок 78 – Оптическая схема с визуальным представлением потока
Такое расположение объектов на оптической схеме позволит получать с
помощью сфокусированных в точке излучений непрерывного и импульсного
лазеров, а также регистрировать непрерывный оптический пробой в точке
пересечения сфокусированных лазерных излучений. Ниже представлена
таблица с описанием всех элементов, входящих в данную схему.
77
Таблица 5 – Элементы оптической схемы
Рисунок 79 – ход лучей, построение изображения
После построения хода лучей на был проведен расчет увеличения
изображения и площади, занимаемой изображением на матрице камеры, при
начальном размере 1 мм, для двух объективов с фокусными расстояниями 135
и 50 мм, выбор объектива:
78
Размеры изображения для объектива с фокусным расстоянием 135 мм
больше, таким образом данный объектив будет выбран для фиксирования
происходящего процесса.
Рисунок 80 – Оптическая схема объектива Юпитер -37а 135 мм
80
4.2
Источники лазерного излучения
В качестве мощного непрерывного лазерного источника, с выходной
мощностью 1000 Вт был выбран волоконный иттербиевый лазер ЛС-1. С
помощью данного лазера планируется поддерживать стационарный разряд.
Рисунок 81 – Схема корпуса ЛС-1
81
Рисунок 82 – характеристики ЛС-1
Исходя из характеристик и зная фокусное расстояние фокусирующей
линзы на оптическом окне, а также конечный диаметр пучка, можем найти,
какой диаметр пучка должен быть после прохождения коллиматора:
𝐷=
4 × 𝐹 × 𝐵𝑃𝑃
𝑑п
(4.1)
где D – диаметр пучка, F – фокусное расстояние линзы, F=450 мм, BPP
– качество пучка BPP = 2,5 мм*мрад, dп – диаметр в перетяжке, dп = 50 мкм.
82
Рисунок 83 – модель ЛС-1
Для получения оптического разряда понадобиться импульсный
источник излучения, в качестве которого была выбрана модель квантрона К301В. Для этого квантрона используется резонатор Фабри-Перо (зеркала 100%
и 50% пропускания).
Рисунок 84 – Модель квантрона К-301В.
83
4.3
Расчет коллиматоров
Исходя из требуемых параметров к излучению, а также фокусирующей
линзы на оптическом окне по формуле (4.1), были рассчитаны коллиматоры
для непрерывного волоконного лазера и импульсного лазера.
Таблица 6 – результаты габаритных расчетов коллиматора для волоконного
лазера с расчетными формулами, обозначениями и полученными
значениями.
84
Рисунок 85 – Оптическая схема коллиматора для волоконного лазера
Таблица 7 – результаты габаритных расчетов коллиматора для импульсного
лазера с расчетными формулами, обозначениями и полученными
значениями.
85
Рисунок 86 – Оптическая схема коллиматора для импульсного лазера
Таким образом, зная все диаметры и фокусы линз, были определены и
подобраны все комплектующие, сделаны оптические схемы коллиматоров.
(Приложения Б-В)
FL,
Наименование
Приложение
Производитель
6mm Dia. x -9mm
Б
Edmund Optics
Б
Edmund Optics
Б
Edmund Optics
NIR
II
Coated,
Double-Concave Lens
12.0mm Dia.
30.0mm FL,
Coated,
NIR
x
II
Plano-Convex
Lens
12mm Dia. x 12mm FL, NIR II Coated,
Double-Concave Lens
86
9mm Dia. x -9mm
FL,
NIR
II
Б
Edmund Optics
-
заказная
-
заказная
Coated,
Double-Concave Lens
120.0mm Dia. x
90.0mm FL,
Coated,
NIR
II
Plano-Convex
Lens
120.0mm Dia. x
120.0mm FL, NIR II
Coated,
Plano-Convex
Lens
Таблица 8 – Подобранные линзы
87
4.4
Техническое описание СХПК
Система
обеспечивает
хранение
и
подачу
компонентов
в
парогазогенератор и устройство воспламенения, в агрегаты и элементы СВД
(теплообменник, диффузор, пневмозаслонки и др.), поддерживает требуемые
параметры (расход, давление, температура) компонентов во вр емя пуска,
обеспечивает автоматическое управление и контроль аварийных ситуаций.
Начальные параметры для расчета ПГС, полученные из технического
задания, газодинамических расчетов эжектора и
парогазогенератора,
представлены в таблице (9).
Суммарный расход через ПГГ GΣ
- стартовый режим
9 кг/с
- номинальный режим
6,5 кг/с
Расход основного окислителя
- стартовый режим G1
6,19 кг/с
- номинальный режим G2
4,47 кг/с
Длительность пуска
- стартовый режим t 1
2с
- номинальный режим t 2
13 с
Расход азота на срыв эжектора G3
4 кг/с
Длительность срыва эжектора t 3
3с
Расход на наддув баков основного горючего и
разбавителя и на продувки G4
0.2 кг/с
Длительность наддува и продувок t 4
15 c
Расход окислителя устройства воспламенения G5
0.18 кг/с
Длительность работы устройства воспламенения t 5
3c
Количество пусков с одной заправки n1
4
Объем баллона V1
0.4 м3
Начальное давление в баллоне Р н
200 атм
88
Конечное давление в баллоне Р к
100 атм
Таблица 9 - Начальные параметры для расчета ПГС
4.4.1 Две подсистемы СХПК
4.4.1.1
СХПК газовых компонентов
В ПГС СВД газообразные компоненты (воздух и азот) хранятся в
стандартных цилиндрических баллонах емкостью 400 л. (Р раб=400 ат) и 40 л.
(Рраб=150 ат). На этапе были проведены оценки количества газов [28],
необходимого для проведения цикла испытаний и баллонов для их хранения.
Данные представлены в таблице (10).
Компоненты
Начальная
масса Расходная
компонента в баллонах, масса
кг
Количество
баллонов
компонента, кг
Основной воздух 665
290,62
7 (400 л)
Управляющий
13,3
6,7
2 (40 л)
95
53,5
1 (400 л)
азот
Азот
Таблица 10 - оценки количества газов и баллонов для их хранения
4.4.1.2
СХПК жидких компонентов
Суммарный расход через ПГГ GΣ
- стартовый режим
9 кг/с
- номинальный режим
6,5 кг/с
Расход основного горючего (спирт)
- стартовый режим G1
0,74 кг/с
- номинальный режим G2
0,54 кг/с
89
Расход разбавителя (вода)
- стартовый режим G1
2,0 кг/с
- номинальный режим G2
1,5 кг/с
Длительность пуска
- стартовый режим t 1
2с
- номинальный режим t 2
13 с
Давление наддува баков Р над
70 атм
Плотность горючего ρ г
780 кг/м3
Плотность разбавителя ρ р
1000 кг/м3
Расход
горючего
устройства
воспламенения 0,01 кг/с
(бензин) Gув
Давление наддува бака устройства воспламенения 60 атм
Рув
Плотность горючего устройства воспламенения ρ ув 760 кг/м3
Расход воды в теплообменник Gто
600 л/мин
Количество пусков с одной заправки n1
4
Таблица 11 - Начальные параметры и технические требования к СХПК
В ПГС используется вытеснительная система подачи [29]. Жидкие
компоненты хранятся в баках, рассчитанных на высокое давление. Требуемое
количество жидких компонентов и объемы баков для проведения заданного
непрерывного цикла исследований представлены в таблице (12).
Компоненты
Начальная
масса Расходная
компонента
в масса
баллонах, кг
Количество
баллонов
компонента, кг
Основное горючее 37,4
34
1 (60 л)
Разбавитель
94
1 (124 л)
103,4
90
Горючее
устройства
0,22
0,2
1 (0,4 л)
660
600
1 (800 л)
воспламенения
Вода
теплообменника
Таблица 12 - Требуемое количество жидких компонентов и объемы баков
4.4.2 Пневмогидравлическая система СВД
Принципиальная схема ПГС СВД приведена на рисунке (87), где
выделены следующие основные составляющие системы:
−
система подачи воздуха к ПГГ и другим агрегатам
−
система подачи топлива к ПГГ
−
система подачи топлива к УВ
−
система подачи воды к ПГГ
−
система подачи воды к теплообменнику
−
система подачи азота в диффузор и на срыв эжектора
Пневмогидравлическая система состоит из следующих элементов:
- Воздушный компрессор
- Баллон высокого давления емкостью 400 л. (B1-B8)
- Баллон высокого давления емкостью 40 л. (B9-B10)
- Вентиль (V1-V44)
- Редуктор (R1-R11)
- Электропневмоклапан (EV1-EV16)
- ПЭК ДД (1-22)
- Клапан отсечной (PV1-PV22)
- Клапан обратный (BV1-BV9)
- Клапан предохранительный (SV1-SV3)
- Ресиверы (Rc1-Rc9)
- Дроссельные устройства (T1-9)
91
- Фильтр топливный (F3, F4, F5)
- Фильтр воздушный (F1,F2)
- Манометр (M1-M3)
- Датчик давления (P1-P57)
- Трубопровод (Ø 50; 32; 20; 10; 4 мм)
Liquid Fuel Supply System
V12
V23
Air-Supply System
Rc2 EV2
F1
P2
32
R2
EV12
High pressure air bottles
V13
V22
SV1
P3
EV3
Rc3
B1-B7
R3
EV13
V1
V10
P9
F
V8 3
B9-B10
PV4
10
N2
PH1
BV2
V7
V25
50
P22
N2
V19
V31
F2 V9
V5
FM1
1
2
3
4
5
6
7
8
9
10
:
22
Electrical double-acting valves
PV2 V2
PV3
P8
R9
PV6 BV1
Air compressor
BV9
M1
Rc1
EV1
4
P21
V40
20
V11
Liquid Fuel Igniter Supply System
4
PV21 V4
4
EV16
P56
P12
F4
Distilled
water
PV9
V18
EV14
V15
P13
P50-55
(movable)
BV5
PH2
BV4
Water Supply System
V16 P6
R5
Rc6 EV6
FM3
50
20
SV3 V26
V17
TC7
EV10
Rc7 EV7
P7
R7
PV14
V38
FM2 PV13
TC5
TC8
P45
P44HE
TC2
GG
P18
10
BV7
P26
V44
Ign.
P16
T2
P11
P27
P35
T9
P32 P41-TC4
Mirrors shutter
P43 SD
P30
Rc5
EV5
P5 DEV1
P37-P40
DEV2
(Pitot)
V29
EV8
LC
F5
V28
P36
TC6
P31
P17
V37
P33
R6
EV9
10
PV11
P24
V32
EN
P46
V3 P29
32
T1
TC3
P57
V21 V20
50
BV3
PV1
PV18
P47-49
PV12 V6
PH3
PV7
T4
PV5
Boosters
P25
T3
Rc8
PV19
P10
PV10
V35
P23
10
PUMP
T8
V39
P4
PV17
V41
4
32
R10
V24
V27
32
Rc4
EV4
4
SV2
V33
P1
V43
R11
V42
R4
PV22
EV11
Rc9
EV15
V14
V30
R1
32
Main air
P20
V36
TC1
P34
T6 10
P14 T5
50
PV20
PV15
V34
PV16
P28
P19
T7
Unstart nitrogen
4
P15
20
PV8
R8
BV8
10
BV6
Рисунок 87 - Принципиальная схема ПГС СВД
Таким образом, все компоненты ПГС СВД были рассмотрены. Из всех
элементов можно выделить измерительные датчики: манометр (M1-M3),
датчик давления (P1-P57).
92
4.5
Погрешности измерительных приборов в ПГС
4.5.1 Манометры [30]
Данные датчики давления, используемые в работе предназначены для
контроля давления в трубках, предназначенных для подвода газа (воздуха и
N2) из баллонов в парогазогенератор. Расположение датчиков представлено
на рисунке (88) (М1-М3).
Рисунок 88 - Принципиальная схема ПГС.
Начальный объем 8-ми используемых баллонов - 400 л, с рабочим
давлением 400 атм = 4,053*10^7 Па. Тогда плотность воздуха в баллоне будет
определена по формуле:
93
(4.2)
— плотность воздуха, M — молярная масса (29 г/моль = 0,029 кг/моль для
сухого воздуха), p — абсолютное давление, R = 8,314 — универсальная газовая
постоянная, T — абсолютная температура в кельвинах. Таким образом,
подстановкой получаем:
4,053 ∗ 107 ∗ 0,029
кг
𝜌=
= 482,5 3
8,314 ∗ 293
м
Тогда масса воздуха равна 482,5 * 0,4 = 192 кг
Суммарный расход газа меняется в диапазоне от 6,5 кг/с до 9 кг/с.
Длительность пропуска от 13 до 3 с. Сечение трубки, по которой идет газ из
баллонов 50 см = 0,5 м, длина 10 м, следовательно при расходе 6,5 кг/с
давление за 13 с давление в трубке будет составлять
Рассмотрим ситуацию при расходе 9 кг/с и времени 3 с:
За 3 секунды 1 баллон успевает накачать 27 кг воздуха в трубку, тогда
давление в этой трубке будет составлять:
𝑝=
27 ∗ 8,314 ∗ 293
= 453,6 кПа
0,5 ∗ 10 ∗ 0,029
Тогда для 8 баллонов: 453600 * 8 = 3,63 МПа
Аналогично для ситуации при расходе 6,5 кг/с и времени 13 с:
За 13 секунд 1 баллон успевает накачать 84,5 кг воздуха в трубку, тогда
давление в этой трубке будет составлять:
𝑝=
84,5 ∗ 8,314 ∗ 293
= 1,419 МПа
0,5 ∗ 10 ∗ 0,029
94
Тогда для 8 баллонов: 1419000 * 8 = 11,356 МПа
Таким образом, для данной работы подойдет манометр МТИ-100-ДИИМ16М с верхним пределом 16 МПа, из перечисленных в таблице (13).
Таблица 13 - коды моделей, верхние пределы измерений P B, диапазоны
шкального индикатора и максимальные (испытательные) давления Рисп.
95
Таблица 14 - Пределы допускаемой основной приведенной погрешности
По таблице (14) пределов допускаемой погрешности достаточно
выбрать класс точности D06, то есть базовое исполнение, таким образом
предел допускаемой приведенной погрешности равен +- (0,3+0,3*Р/Рв).
Следовательно, для расхода 9 кг/с приведенная погрешность равна +(0,32%), а при 6,5 кг/с +-(0,51%).
Таким образом, зная приведенную погрешность можно вычислить и
абсолютную по формуле
(4.3)
Тогда ∆𝑥 = 𝑥𝑁 ∗ 𝛾 = {
4540,8 Па = 0,045 атм, при 9
кг
с
57915,6 Па = 0,571 атм, при 6,5
Ниже, на рисунке (89) представлен чертеж манометра
96
и3с
кг
с
и 13 с
Рисунок 89 - Габаритные, присоединительные и монтажные размеры
манометров электронных для точных измерений МТИ-100/М1
4.5.2 Преобразователи давления [31]
Данные датчики нужны для отслеживания давления на выходе. В работе
используются давления от 0,01 атм до 100 атм. В паспорте устройства (таблица
15) была выбрана подходящая конфигурация, которая способна отслеживать
данные значения.
97
Таблица 15 - Основные метрологические характеристики преобразователей
давления измерительных «ЭЛЕМЕР-АИР-30М».
Были выбраны модели TG16 с максимальным пределом 60 Мпа и
минимальным 0,6 МПа для 100 атм и TAH9 для давлений 1 атм с
максимальным значением 250 кПа и минимальным в 6 кПа.
В таблице (16) показана приведенная погрешность данного устройства
98
Таблица 16 - Пределы допускаемой основной приведенной погрешности
аналогового выхода, выраженные в процентах от диапазона измерений.
Для выбранных моделей предел допускаемой приведенной погрешности
равен +-(0,2*РBMAX/PB).
Таким образом, приведенная погрешность на данных устройствах при 1
атм = 0,493%, а при 100 атм = 1,184%.
В паспорте устройства прописаны температурные погрешности
(таблица 17), которые стоит учесть, так, для выбранных моделей она
составляет +-(0,03+0,2*РBMAX/PB), при +10 С: 1 атм = 0,523%, 100 атм = 1,214%.
Дополнительная погрешность АИР-30М, вызванная изменением
температуры окружающего воздуха от нормальной (20+-5 С) до любой
температуры в пределах рабочих температур на каждые 10 С изменения
температуры, γт в %/10 С, не превышает значений, приведенных в таблице
(17).
99
Таблица 17 – Дополнительная температурная погрешность.
Дополнительная погрешность АИР-30М, вызванная воздействием
повышенной влажности, не превышает 0,2 предела допускаемой основной
погрешности
Зная, какие факторы могут влиять на данные модели, можно составить
полную относительную погрешность, для 1 атм (+-5 С и повышенная
влажность) = 0,493 + 0,523 + 0,099 = +- 1,115 %, для 100 атм (+- 5 С и
повышенная влажность) = 1,184 + 1,214 + 0,237 = +- 2,635 %
Таким образом абсолютная погрешность выбранных приборов равна:
при 1 атм: 250кПа ∗ 1,115 = 2787,5 Па = 0,0275 атм
{
при 100 атм: 60 МПа ∗ 2,635 = 1,581 МПа = 15,603 атм
Ниже представлен чертеж с габаритными размерами устройства (90)
100
Рисунок 90 - Габаритные, присоединительные и монтажные размеры
преобразователей давления измерительных «ЭЛЕМЕР-АИР-30М».
101
4.6 Вывод по пункту 4
Основываясь на
параметрах лазерного излучения, разработаны
оптические схемы для подвода излучения импульсного и непрерывного
лазеров,
для
регистрации
непрерывного
оптического
разряда
в
высокоскоростном потоке также разработана оптическая схема. Подобраны
комплектующие экспериментальной установки
102
5
Компоновка экспериментальной установки в лаборатории.
Рисунок 91 – Компоновка элементов оптической схемы в лаборатории.
На рисунке (91) представлено расположение всех комплектующих в
лаборатории. Видно, что элементы вписываются в габариты данного
помещения, таким образом, это место подходит для проведения эксперимента,
и дальнейшего изучения непрерывного оптического разряда на сверхзвуковом
газовом потоке.
103
Рисунок 92 – Расположение элементов на оптическом столе (вид сверху)
Рисунок 93 – Расположение элементов на оптическом столе (вид справа)
104
Рисунок 94 – Расположение элементов на оптическом столе (вид спереди)
Рисунок 95 – Расположение элементов на оптическом столе
Все оптические элементы были размещены на столе согласно
предложенной ранее оптической схеме (78). Ниже представлена выхлопная
труба с оптическими окнами, схема данной установки представлена в
приложении Г.
105
Рисунок 96 – Выхлопная труба в разрезе с изображением хода лучей
Таким образом, в лаборатории были расставлены все элементы по
оптической схеме. Элементы, которые были использованы в постановке:
оптоволоконный иттербиевый лазер лс-1,чиллер, коллиматор, 2 объектива со
стойками, камера, выхлопная труба с четырьмя оптическими окнами, два
оптических окна представляет из себя линзу, которая фокусирует непрерывное
и импульсное излучения в пучок диаметром 50 мкм, в конструкции данных
окон предусмотрена подстройка по всем осям с помощью болтов, движение
обеспечено использованием сильфона.
106
ЗАКЛЮЧЕНИЕ
1.
В работе построена математическая и численная модель
непрерывного оптического разряда в сверхзвуковом газовом потоке в
условиях равновесия. Получены аппроксимационные формулы всех
параметров для данной модели.
2.
Приведены результаты моделирования равновесной модели. Для
мощности непрерывного излучения порядка 1 кВт существуют условия, при
которых стационарный оптический разряд в сверхзвуковом газовом потоке с
числами Маха до 5 существует. При наличии основного непрерывного
лазерного пучка самопроизвольно разряд не возникает, а для его зажигания
необходимо внешнее импульсное энергетическое воздействие.
3.
разряда в
Предложена математическая модель непрерывного оптического
высокоскоростном газовом потоке
в
двухтемпературном
приближении, на основе которой возможно моделирование процесса.
Приведены важные оценки и методики расчетов. Рассмотрена траектория
движения свободного электрона вокруг иона в полностью ионизированной
плазме для определения сечения взаимодействия между частицами.
4.
Предложена
методика
экспериментального
исследования
непрерывного оптического разряда в высокоскоростном потоке, позволяющая
проверить результаты моделирования. Разработана оптическая схема
установки и компоновка всех элементов в лаборатории.
107
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ
1)
В.Н. Зудов, П.К. Третьяков - Взаимодействие оптического разряда
с ударной волной, журнал технической физики, 2018, том 88, вып. 3
2)
Юрьев А.С., Пирогов С.Ю., Рыжов Е.В. Управление обтеканием
тел с использованием подвода лазерной энергии в высокоскоростные потоки
газа. — СПб: ВКА им. А.Ф. Можайского, 2005. — 159 c.
3)
Третьяков П.К., Грачев Г.Н., Иванченко А.И., Крайнев В.Л.,
Пономаренко А.Г., Тищенко В.Н. Стабилизация оптического разряда в
сверхзвуковом потоке аргона // Доклады РАН. — 1994. — Т. 336, № 4. — C.
466-467.
4)
Зеркл Д.К., Шварц С., Мертогул Э., Чжэнь С., Криер Х., Мазумдер
Дж. Исследование лазерной аргоновой плазмы применительно к проблеме
электротермического ракетного двигателя // Ракетная техника и космонавтика.
— 1990. — № 11. — C. 3-13.
5)
Тищенко В.Н., Гулидов А.И. Ускорение лазерной плазмы
оптическим пульсирующим разрядом, движущимся в воздухе с гиперзвуковой
скоростью // Письма в ЖТФ. — 2000. — Т. 26, № 19. — C. 77-83.
6)
Я.В. Ханхасаева, В.Е. Борисов, А.Е. Луцкий - Энергетическое
воздействие на обтекание гиперзвуковых летательных аппаратов
7)
Maker P. D., Terhune R. W., Savage С. М. — In: III Intern. Conference
on Quantum Electronics — Paris, 1963.
8)
Островская Г. В., Зайдель А. Н., Лазерная искра в газах, "УФН",
1973, т. 111, с. 579
9)
Райзер Ю.П. О возможности создания оптического плазмотрона и
необходимой для этого мощности // Письма в ЖЭТФ. — 1970. — Т. 11, № 3.
— C. 195-199.
10)
Райзер Ю. П. — Лазерная искра и распространение разрядов,
Москва, «Наука», 1974.
11)
Райзер Ю. П. — Оптические разряды, "УФН", 1980.
108
12)
Крайко А.Н. Аналитическое представление термодинамических
функций воздуха. Инженерный журнал т.4, вып.3, 1964, стр. 548-550.
13)
И. В. Авилова, Л. М. Биберман, В. С. Воробьев, В. М. Замалин, Г.
А. Кобзев, А. Н. Лагарьков, А. Х. Мнацаканян, Г. Э. Норман - Оптческие
свойства горячего воздуха, под редакцией профессора Л. М. Бибермана, издво «Наука», Москва, 1970.
14)
Краснов Н.Ф. Аэродинамика тел вращения. М.: Машиностроение,
1964. – 572 с.
15)
Ключарев А.Н., Мишаков В.Г., Тимофеев Н.А. Введение в физику
низкотемпературной плазмы.
16)
Н.В. Василенко, В.А. Макаров “Модели оценки надежности
программного обеспечения”, 2004 г.
17)
Сивухин Д. В. – Общий курс физики. — М.: Физматлит; Изд-во
МФТИ, 2004. — Т. III. Электричество. — С. 17. — 656 с. — ISBN 5-9221-02273.
18)
Королев Ю.Д. – Элементарные и кинетические процессы в
газоразрядной
плазме,
издательство
Томского
политехнического
университета, 2008 г.
19)
Бронштейн И.Н. – Гипербола/ Квант. — 1975. — № 3.
20)
Ландау Л.Д., Лифшиц Е.М. – Теоретическая физика, том 10,
Москва, «Наука», главная редакция физико-математической литературы,
1979 г.
21)
Two‑Temperature Chemical Non‑equilibrium Modeling of Argon DC
Arc Plasma Torch, Jiang‑Hong Sun, Su‑Rong Sun, Li‑Hui Zhang, Hai‑Xing Wang
13
July
2020,
Plasma
Chemistry
and
Plasma
Processing,
https://doi.org/10.1007/s11090-020-10108-9
22)
Mitchner M., Kruger Charles H., JR. – Partially ionized gases/
Department of mechanical engineering, Stanford university, New York, 1973 г.
23)
База
данных
сечений
столкновений
–
https://nl.lxcat.net/data/set_databases.php (дата обращения 2.05.2021)
109
URL:
24)
Collision cross sections for identical and non-identical rare-gas atom
pairs for energies from 0.01 eV to 10 keV. A. V. Phelps JILA, University of
Colorado and National Institute of Standards and Technology, Boulder, Colorado
80309-0440 (Dated: December 5, 2004)
25)
Chemical Non‑equilibrium Simulation of Anode Attachmentof an
Argon Transferred, Su‑Rong Sun, Hai‑Xing Wang, Tao Zhu, Anthony B. Murphy,
30
August
2019,
Plasma
Chemistry
and
Plasma
Processing.
https://doi.org/10.1007/s11090-019-10027-4
26)
Программа T&TWinner – URL: http://ttwinner.free.fr/ (дата
обращения 15.04.2021)
27)
Short Communication Simple Approximation For Thermal Diffusion In
Gas Mixtures, John D. Ramshaw, Idaho National Engineering Laboratory, Idaho
Falls, U.S.A., J. Non-Equilib. Thermodyn. Vol. 21 (1996),
28)
Володин В.А. – Конструкция и проектирование ракетных
двигателей – М.: Машиностроение, 1971 г.
29)
Беляев Н.М. – Расчет пневмогидравлических систем ракет – М.:
Машиностроение, 1983 г.
30)
Манометры электронные для точных МТИ-100/М1, МТИ-100/М3,
руководство по эксплуатации НКГЖ.406233.058РЭ
31)
Преобразователи давления измерительные «ЭЛЕМЕР-АИР-30М»,
руководство по эксплуатации НКГЖ.406233.064РЭ
32)
Transition radiation energy loss in inductively coupled argon plasma,
Yuping Ren, J. D. Long and S. Xu, Plasma Sources and Applications Center, NIE,
and Institute of Advanced Studies, Nanyang Technological University, 1 Nanyang
Walk, Singapore 637616, 11 July 2007, doi: 10.1063/1.2749235.
110
Приложение А
1. Чертеж оптической схемы (БГТУ.015113.000)
111
Приложение Б
Закупные линзы коллиматоров:
Производство Edmund Optics
1.
6mm Dia. x -9mm FL, NIR II Coated, Double-Concave Lens
2.
12.0mm Dia. x 30.0mm FL, NIR II Coated, Plano-Convex Lens
3.
12mm Dia. x -12mm FL, NIR II Coated, Double-Concave Lens
4.
9mm Dia. x -9mm FL, NIR II Coated, Double-Concave Lens
112
Приложение В
Оптические схемы коллиматоров для импульсного и непрерывного лазеров:
1. Оптическая схема коллиматора для непрерывного лазера
(БГТУ.015113.001)
2. Оптическая схема коллиматора для импульсного лазера
(БГТУ.015113.002)
113
Приложение Г
1.
Сборочный чертеж выхлопной трубы с оптическими окнами
(БГТУ.015113.003)
114
Отзывы:
Авторизуйтесь, чтобы оставить отзыв