19221

ТЕРМОДИНАМИЧЕСКОЕ РАВНОВЕСИЕ В ГАЗЕ

Лекция

Физика

Лекция 2 ТЕРМОДИНАМИЧЕСКОЕ РАВНОВЕСИЕ В ГАЗЕ Одним из известных подходов к описанию плазмы является ее сопоставление с термодинамической системой. При этом состояние плазмы характеризуется такими величинами как температура энтропия и т.д. В термодинамик...

Русский

2013-07-11

101.5 KB

10 чел.

Лекция 2

ТЕРМОДИНАМИЧЕСКОЕ РАВНОВЕСИЕ В ГАЗЕ

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

       Для модели ПТР требуется выполнение следующих требований: 1) распределения Максвелла по скоростям частиц (ионов и электронов), 2) распределения Больцмана для населенностей уровней частиц, 3) распределения Планка для излучения, 4) распределения Саха для концентрации заряженных частиц. В качестве примера лабораторного источника плазмы, для которого пригодна модель ПТР, можно отнести капиллярный разряд конструкции Подмошенского. В хорошем соответствии с моделью ПТР находится фотосфера Солнца.

       Для модели ЛТР не требуется выполнение необходимых условий во всем объеме плазмы. Подразумевается выполнение следующих распределений в небольшом локальном объеме плазмы: 1) распределения Максвелла по скоростям частиц, 2) распределения Больцмана для населенностей уровней частиц, 3) распределения Саха для концентрации заряженных частиц.  Модель ЛТР широко применяется для разных видов лабораторной и космической плазмы.

       Изложим понятия, необходимые для определения соответствия изучаемой плазмы той или иной модели. Рассмотрим вопрос, связанный с определением оптической толщины плазмы. Для этой цели вводится величина – длина свободного пробега фотона  lп. Ее значение обычно определяется экспериментально, т.к. расчет достаточно сложен. В зависимости от соотношения длины пробега  lп и размеров плазмы  L, плазменная среда может быть оптически тонкая и толстая. Для оптически тонкой плазмы длина пробега фотона должна быть больше размеров плазмы:  lп >L. В случае оптически толстой плазмы реализуется обратная ситуация:  lп <L. Условие оптически толстой плазмы означает, что фотон, возникший во внутренней области плазмы на пути к поверхности, может испытать многократные поглощения, т.е. возникает процесс переизлучения фотона. Анализ контуров спектральных линий в данном случае может дать информацию о сильном искажении контура линии, как, например, об уменьшении интенсивности и провале в центре линии. В случае ПТР плазма должна быть обязательно оптически толстой.

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

       Для количественного описания ионизационного равновесия известным индийским астрофизиком Мегнадом Саха в 1920 г. была получена формула, характеризующая зависимость концентрации плазмы от температуры и энергии ионизации в случае водородной плазмы. При выводе данной формулы предполагается, что плазма достигла состояния термодинамического равновесия. Допустим, что ионизация водорода осуществляется из основного состояния (n=1) в непрерывный спектр (рис.1). Это, конечно, является упрощенным подходом, т.к. не учитывается ионизация из других состояний (n>1). Для вероятности нахождения электрона в состоянии с энергией En считается справедливым распределение Гиббса:

                        

                                                                               

                                                                      Рис.1

                                                                                    

                                                                                                                      

Для отношения вероятностей нахождения электрона в состоянии непрерывного спектра с энергией E2 и в основном состоянии с энергией E1 можно записать следующее соотношение:

              

В данной формуле основное состояние считается невырожденным и его статистический вес равен g1=1. Для нахождения g2 воспользуемся формулой для числа состояний фазового пространства непрерывного спектра:

              

Где s – число степеней свободы. Числитель данной формулы записывается следующим образом:

               ,        ,                   

Запишем g2 с учетом этих формул:

             

При данном рассмотрении предполагается, что для плазмы реализуется распределение Больцмана. Поэтому концентрация частиц в определенном состоянии пропорциональна вероятности, т.е. можно записать следующее отношение:

             

В результате формула Саха для концентрации водородной плазмы записывается в следующем виде:

                    

Ввиду квазинейтральности плазмы (neni) данную формулу часто переписывают, перенося величину  ne из правой части в левую. Тогда в правой части остаются величины, зависящие от Т, т.е. множитель с температурой в степени 3/2 и экспоненциальный множитель, содержащий температуру.

       Рассмотрим основные предположения, которые используются для вывода формулы Саха для ионов. Допускается наличие распределения Больцмана для населенностей иона с зарядностью  z  и с зарядностью  z+1. Предполагается, что прямая ионизация может происходить как из основного состояния иона (n=1), так и из других состояний с большей энергией (n>1). Поэтому в окончательную формулу подставляется статистическая сумма G, содержащая произведения статистических весов отдельных уровней и экспоненциального множителя из распределения Больцмана.

                   

          В результате формула Саха для ионов будет иметь вид:                 

                  

                       

       Зависимость концентрации плазмы от температуры в случае равновесной плазмы, т.е. формула Саха, позволяет получить количественное выражение для такой важной характеристики как степень ионизации плазмы. Формула для степени ионизации плазмы имеет вид:

                     

        - концентрация электронов, - концентрация атомов

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

                 

                   ,         [n]=см-3,    []=эВ

                                

                                                                 Рис.2

                                                                                                                  

Представим графические зависимости степени ионизации от температуры, полученные с помощью данной формулы для цезия, водорода и гелия (рис.2). Для расчетов использовалась концентрация газа равная n=1016 см-3. Самым легко ионизуемым газом является цезий, у которого полная ионизация (1) наступает практически при Т1 эВ. Для водорода полная ионизация осуществляется при Т1,8 эВ, а для гелия при Т3 эВ. Для большей концентрации  n=1017 см-3 полная ионизация водорода  1 согласно расчету наступает даже при Т0,16 эВ=1850 К. Приведенные расчеты температур являются несколько завышенными относительно реальных экспериментов, т.к. не учитывается влияние примесей, существующих в плазме и приводящих к уменьшению реальных значений температуры.                                                                              

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

                     

                     

Скорость ионизации ([Qi ]= (част/cсм3)) имеет выражение:

                           

Где  ki  константа ионизации, na – концентрация атомов, ne – концентрация электронов.                   

Скорость рекомбинации ([Qr]= (част/cсм3)) запишется в виде:

                     

Где  kr константа рекомбинации, ni – концентрация ионов.

В стационарном состоянии, когда реализуется равновесие между процессами ионизации и рекомбинации, данные скорости можно приравнять:

                                         

В результате отношение концентраций ионов к концентрации атомов равно отношению констант ионизации  и рекомбинации, что составляет формулу Эльверта:                      

                                           

Следует заметить, что данные константы имеют зависимость от температуры    и.

 


n=2

n=1

1, g1

E2, g2

Eи

4

3

2

1

T, эВ

1,0

0

3

2

1

0,5

     n=1016 см-3

1) Cs (E1=3,9 эВ)

2) H (E2=13,6 эВ)

3) He (E3=24,5 эВ)


 

А также другие работы, которые могут Вас заинтересовать

41968. Дослідження стійкості ланки другого порядку 114.05 KB
  Для лінійних систем автоматичного керування, які описуються характеристичним рівнянням виду a0pn+a1pn-1+…+an-1p+an=0 стійкість не залежить від величини і вигляду збурення і визначається коренями характеристичного рівняння, яке залежить від параметрів системи Для зручності зафіксуємо L C та змінюватимемо R withinttrns; urovnenie:=TTpp2xiTp1; h:=k p urovnenie; l:=invlplcehpt; sol:=solveurovneniep: sol[1];sol[2]; Аперіодичний процес Вибираємо L=50мГн.05;C:=2010^6;R:=250;T:=sqrtLC;xi:=RsqrtC L 2;k:=1;p1:=sol[1];p2:=sol[2];задання параметрів для даного виду процесу l:=invlplcehpt;розрахунок зворотнього перетворення Лапласа plotlt=0.05;C:=2010^6;R:=100;T:=sqrtLC;xi:=RsqrtC L 2;k:=1;p1:=sol[1];p2:=sol[2]; l:=invlplcehpt:...
41969. ДОСЛIДЖЕННЯ ВЕКТОРНИХ ПЛОТТЕРІВ 78.68 KB
  Все рассматриваемые здесь команды находятся в основной части языка HPGL 2. Первыми идут команды ини рйализации для установки размера изображения и другие параметры после них следуют команды для прорисовки линий фигур и трок символов а также одна или две команды для завершения процесса. Некоторые команды имеющие числовые аргументы требуют целых значений в то время как другие команды допускают наличие чисел с десятичной точкой. Некоторые команды передают результаты обратно хосткомпьютеру: например 01 сообщает идентификацию модели...
41970. Функції введення/виведення printf(), scanf(). Лінійні обчислювальні процеси 14.14 KB
  Обладнання: ПКПО Borlnd C Хід роботи Вивчити теоретичні відомості Ознайомитися з форматом функцій printf і scnf.h void min { long ; double b; unsigned c; flot d; cout Вводите n ; cin b c d; cout Long n ; printf 16.5d ; printf n double n ; printf 16.
41971. Обчислювальний процес, що розгалужується, з різними логічними умовами: оператор if... else, умовна операція (?:), оператор switch, оператор break, оператор goto 23.72 KB
  else умовна операція : оператор switch оператор brek оператор goto Ціль роботи: Вивчити реалізацію в мові ветвящихся обчислювальних процесів . Навчитися писати програми використовуючи оператори: розгалуження if.else переключення switch у парі з оператором brek оператор переходу goto тернарную умовну операцію .
41972. Розробка програм з циклічними обчислювальними процесами 44.87 KB
  Розробка програм з циклічними обчислювальними процесами Ціль роботи: Вивчити написання програм мовою С, використовуючи ітераційні (циклічні) методи, освоїти основні оператори, що підтримують роботу з циклами (for, while, do... while). Навчитися писати програми, використовуючи дані оператори.
41973. ПОБУДОВА ОПТИМАЛЬНОГО НЕРІВНОМІРНОГО КОДУ ЗА МЕТОДИКОЮ ХАФФМАНА 53.47 KB
  0 проводиться перехід до побудови дерева коду за допомогою проміжних вузлів. 161 00074 3 В 893 00412 21 Х 156 00072 11 Л 745 00344 29 Ю 148 00068 16 Р 699 00322 22 Ц 126 00058 №п п Символ ni pi №п п Символ ni pi 12 М 656 00303 25 Щ 108 00050 10 К 574 00265 24 Ш 60 00028 5 Д 507 00234 28 Э 59 00027 26 Ы 467 00215 20 Ф 30 00014 19 У 399 00184 8 З 4 00002 Дерево коду за методикою Хаффмана: Визначаємо ентропію джерела за формулою: Визначаємо максимальний ступінь стиснення інформації: Середня довжина кодової комбінації:...
41976. Изучение методики процедурного программирования в СУБД 903.17 KB
  Изучение управленческих конструкций IFEndIF и IIF. Изучение конструкций построения циклов For EndFOR. Изучение управленческих конструкций IFEndIF и IIF.