Главная » Просмотр файлов » Диссертация

Диссертация (1150536), страница 7

Файл №1150536 Диссертация (Радиационные эффекты в неравновесной плазме дуговых и тлеющих разрядов) 7 страницаДиссертация (1150536) страница 72019-06-29СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 7)

Подробности вычислений можно найти в Приложении B.3. Для краткости записи, будем использоватьследующие обозначения+Δ+ (,) = +1(,) − + (,),−Δ− (,) = − (,) − +1(,),± (,) → ± , ( ) → ,(2.34) = () → = (0).Выражения для коэффициентов имеют вид⎧(︁)︁(1)(2)⎪⎪, = , + , ,⎪⎪⎪⎪⎡⎤⎪⎪⎪∫︁⎪⎪1 11⎪⎢⎥(1)⎪⎪, = 3/2 · ⎣ √︁− √︁⎦d,⎪⎪(1)(1)++⎪, + +1, + ⎪⎪⎪⎪⎪−1⎪∑︁⎨ (1)++ Δ− + , < : , = +1 +⎪=+1 ⎡⎤⎪⎪⎪+1⎪−1 ∫︁⎪∑︁⎪111⎢⎥(2)⎪⎪√︁√︁=·−⎣⎦d,⎪,3/2 ⎪(2)(2)⎪++⎪=0 ,, + +1,, + ⎪⎪⎪⎪−1−1⎪∑︁∑︁⎪(︀ +)︀⎪(2)+−−−+⎪=Δ+Δ+−⎪+1+1 + − ,,,⎪⎩=+1=+1(2.35)40⎧(︁)︁(3)(4)(5)⎪⎪, = , + , + , ,⎪⎪⎪⎪⎪⎪∫︁⎪⎪11(3)⎪⎪√︁=·d,⎪,3/2⎪+⎪⎪ +1⎪⎪⎪⎪⎪∫︁⎪⎪11⎨ (4), = 3/2 · √︀ − d,=: ⎪0⎪⎤⎡⎪⎪+1⎪−1 ∫︁⎪∑︁⎪111⎥⎢⎪(5)⎪√︁√︁=·−⎪⎦d,⎣,⎪3/2(5)(5)⎪++⎪=0 ,, + +1,, + ⎪⎪⎪⎪−1−1⎪∑︁∑︁⎪(︀ +)︀⎪(5)+−−⎪⎪=Δ+Δ+−+ − − + ,+1+1⎪,,⎩=+1=+1(2.36)⎧)︁(︁(8)(7)(6)⎪⎪, = , + , + , ,⎪⎪⎪⎪⎡⎤⎪⎪+1⎪∫︁⎪⎪1 11⎢⎥⎪(6)⎪− √︁, = 3/2 ·⎪⎣ √︁⎦d,⎪(6)(6)⎪+−⎪, + +1, + +1⎪⎪⎪⎪−1⎪∑︁⎪⎪(6)−−⎪⎪, = Δ+ − − +1,⎪⎪⎪=+1⎪⎡⎤⎪⎪⎪+1⎪−1 ∫︁⎪⎨ (7)111 ∑︁⎢⎥√︁√︁=·−⎣⎦d,,3/2>:(7)(7) =0−−⎪,, + ,, + +1⎪⎪⎪⎪−1⎪∑︁⎪(7)⎪−−⎪,, = Δ+ − − +1,⎪⎪⎪⎪=+1⎪⎡⎤⎪⎪+1⎪∫︁−1⎪⎪1 ∑︁11⎪⎢⎥(8)⎪·=− √︁⎪⎣ √︁⎦d,,⎪3/2⎪(8)(8) =0++⎪,, + +1,, + ⎪⎪⎪⎪−1−1⎪∑︁∑︁⎪(︀ +)︀⎪(8)+−−−+⎪⎪=Δ+Δ+− +1⎪+1 + − .⎩ ,,=+1=+1(2.37)41Переход к равномерной сетке можно осуществить, вынося√︀/ за знак инте-грала и переобозначая ± :± () =(︂)︂1+cos ±2√︃(︂)︂212 − +sin2 .2(2.38)Представляя коэффициент поглощения в виде 0 ( ) = ˜ 0 ( )·0 (0), где 0 (0) - значение коэффициента в центре объёма, можно вынести эффективную вероятностьперехода:eff =.0 (0)Приведенные формулы для неоднородного коэффициента поглощения в бесконечном цилиндре будут использованы далее в Главе 5 при моделировании контрагированного положительного столба и в Главе 3 при моделировании дуговогоразряда.Влияние неоднородности коэффициента поглощения на решение уравненияХолстейна-БиберманаПусть ищется решение системы уравнений в следующем видеeff∑︁ · ( ) = ( ), () = 0 · ( = 0),(2.39)где 0 = 106 [см−1 ].Для того, чтобы наглядно проиллюстрировать влияние неоднородности поглощающих атомов (здесь - нейтральных) на распределение излучающих (в данномслучае - резонансных), можно провести анализ функций Грина.∫︁ Функции Гринапозволяют найти решение уравнения (2.39) в виде () =(,′ ) (′ )d3 ′ .Таким образом, функция Грина соответствует обратной матрице −1 , подразумеваяпереход от дискретных переменных к непрерывным.Рассмотрим несколько случаев неоднородности коэффициента поглощения,вводя для наглядности функции (), нормированные на усредненное по радиусу420 .2 50 .1 5(b )(a )0 .2 00 .1 00 .1 0G ( r , r ')G ( r , r ')0 .0 50 .1 50 .0 51100.

6. 8. 61100. 88r'0.60.60.0 .20 .2r0 . 20 .440 . 20.r0 . 40.8r'0 . 40 .2 5(c )0 .1 50 .1 0G ( r , r ')0 .2 00 .0 510810. 6r'0.60.0 . 4. 840 .20.r0 . 2Рисунок 2.7: Функции Грина оператора переноса излучения для однородного() и неоднородного поглощения при ˜ 0 () = 102 (/)2 /33.34 (), и ˜ 0 () =2410 (/) /3334.33 ().значение avg :˜ 0 () = ()/avg = 1,[︂(︁ )︁2 ]︂12 ˜ 0 () = 1 + 10·,34.33[︂(︁ )︁2 ]︂1˜ 0 () = 1 + 104·.3334.33(2.40)(2.41)(2.42)На Рис.2.7 приведены рассчитанные функции Грина.

В однородном случае функция Грина практически постоянна вдоль линии = ′ , и наблюдается тольконебольшой спад вблизи центра цилиндра и границы . С увеличением неоднородности возникает заметный спад по диагонали, и усиливается асимметрия. Более43-1-21 0-31 0-41 0-52~κ ~ 1 + 1 0 4 ( r / R )02~κ = 10~κ ~ 1 + 1 0 2 ( r / R )0N (r)/N01 001 0~κ = 10N (r)/N (0 )1 01 0-11 0-21 0-31 0-42~κ ~ 1 + 1 0 ( r / R )024(a )0 .0~κ ~ 1 + 1 0 2 ( r / R )0(b )0 .20 .40 .60 .81 .00 .0r/R0 .20 .40 .60 .81 .0r/RРисунок 2.8: Радиальные распределения относительной концентрации резонансных атомов в случае -образного источника возбуждения при наличии неоднородности поглощения.того, наличие неоднородности ведет к большему спаду заселенности излучающегорезонансного уровня на оси разряда (Рис. 2.8a), но выносит больше атомов напериферию разряда (Рис.

2.8b). Здесь 0 = (0)/eff , а (0) - заселенность вточке = 0. Дальнейший анализ влияния неоднородности поглощающих атомовна решение задачи необходимо проводить в рамках конкретных моделей (см. Главу3).2.3Метод трассировки лучейВ предыдущем параграфе был рассмотрен метод дискретизации интегральногоуравнения Холстейна-Бибермана (2.16) на элементарных объёмах, обладающихсимметрией источника. Данный подход весьма удобен для 1D-симметричных моделей, бесконечных по одной из координат и симметричных по другой. Такжеимеются работы, в которых рассчитывались матричные коэффициенты для 2Dсимметричной модели конечного цилиндра с однородным поглощением [172].

Вэтом случае аналитический вывод выражений для коэффициентов являлся довольно громоздким. Представляет интерес метод, который позволит рассматриватьперенос излучения методом конечных объёмов в 3D-источниках произвольнойформы, в том числе, асимметричных.44(x,y,z)ziryρi,k(x’,y’,z’)kxr’Рисунок 2.9: Схема дискретизации на воксельной сетке.2.3.1Схема дискретизации и алгоритм трассировкиПусть произвольная область объёма разбивается на подобъёмов Δ = / , где Δ - элементы равномерной декартовой 3D-сетки (кубические элементы,или воксели). Соответственно, расстояние между точками r и r′ будет определятьсякак√︁|r − r | = = ( − ′ )2 + ( − ′ )2 + ( − ′ )2 · Δ.(2.43)′Здесь ,, имеют смысл целочисленных координат вокселей, а Δ =√3Δ - длинаребра воксела. Можно сопоставить векторам координат (,,), (′ , ′ , ′ ) индексы и для каждой пары элементов:(,,) → ; (′ , ′ , ′ ) → ;√︀, = Δ ( − )2 + ( − )2 + ( − )2, = 1,...,.Данная схема дискретизации проиллюстрирована Рис.

2.9.(2.44)45Для учета всех пересеченных вокселей на пути между элементами и былиспользован быстрый алгоритм трассировки лучей Аманатидеса и Ву [173]. Сутьалгоритма заключается в следующем:1. Задаются начальная и конечная координаты, соответствующие элементам синдексами и .2. Определяются углы наклона луча ∠ = / max ( , , ), ∠ и ∠.3.

Рассчитываются расстояния ( , , ) из исходной точки до граней пересечения с соседними элементами. Минимальная из этих величин определяетрасстояние, которое можно пройти вдоль луча → , при этом оставаясь впределах текущего элемента.4. Затем рассчитываются инкременты для каждой из компонент Δ = Δ/∠,Δ и Δ . К минимальному на данный момент расстоянию min ( , , )прибавляется инкремент, а соответствующая координата увеличивается нафиксированный шаг .5. Производится переход в другой воксел, затем пункты 2 − 5 повторяются,пока не будет достигнута конечная точка в элементе .В результате, на пути движения луча образуется массив пересеченных вокселей(Рис.

2.10). Если рассматривать задачу с однородным поглощением, то достаточнознать и . Для учета неоднородности необходимо сохранять координаты всехпересеченных лучом элементов. Индексы данных элементов записываются в мас,сив , = [, ..., ] и могут быть получены через = 1, , ..., , где - числовокселей в массиве.2.3.2Вычисление матричных коэффициентовПоскольку сетка является декартовой и равномерной, то интеграл по объёму вуравнении (2.16) можно заменить на Δ . Таким образом, основная сложность подхода с симметричными элементарными объёмами, связанная с расчетом пространственных интегралов, устраняется. Для дискретизации ядра () необходиморассмотреть два случая.46J1i,k= iJ2i,kJ3i,kJ4i,kJ5i,kJ6i,kJmi,k= kРисунок 2.10: Алгоритм трассировки лучей на воксельной сетке.

̸= : Здесь ̸= 0, и получаемΔ = ·4∫︁∞0 exp [− ( )]d.2(2.45) = : В этом случае = 0, и в интеграле появляется сингулярность. Онаможет быть устранена путем перехода к сферическим координатам в пределахрассматриваемого элемента, где , → Δ:⎛ ,⎞ ⎤∫︁4 ⎝ − (′ ) ′ 2 ′ ⎠ ⎦= · ⎣1 +e d d4′ 200⎡⎛ Δ⎞ ⎤∞∫︁ ∫︁= · ⎣1 + ⎝ e− () d⎠ d ⎦ .∫︁∞⎡,0Здесь Δ = (3Δ /4)1/3 .0(2.46)47В случае однородного поглощения, итоговые выражения для матричных коэффициентов записываются как,⎧∫︁∞⎪⎪Δ ⎪⎪exp (− Δ)d, ̸= ;⎪⎪2,⎨ 40 ⎡⎤=·∞ ∫︁Δ∫︁⎪⎪⎪⎪1 + ⎣ exp (− )d⎦d, = .⎪⎪⎩0(2.47)0При учете неоднородности интеграл в экспоненте необходимо заменить на суммупо массиву , :,⎛⎞⎧,∫︁∞⎪∑︁⎪Δ () ()⎪⎝⎪exp− ()Δ⎠d, ≠ ;⎪2⎪,⎨ 4=1,0 ⎡⎤=·∫︁∞ ∫︁Δ⎪⎪⎪⎪⎪1 + ⎣ () () exp (− ())d⎦d, = .⎪⎩0(2.48)0Указанные выражения коэффициентов справедливы для произвольного контураспектральной линии.

Используя функции () и (), можно легко показать,что формулы сводятся к более простым в случае лоренцевского контура и егоасимптотики.Для лоренцевского контура(︁ )︁ (︁ )︁00 () = exp −0,22(︁ )︁ [︁ (︁ )︁(︁ )︁]︁01000·exp−−,() =0142 2222(2.49)(2.50)и коэффициенты, соответственно,,⎧(︁ )︁ [︁ (︁ )︁(︁ )︁]︁0 Δ0 ,0 ,0 ,⎪⎪exp−−, ≠ ;01⎨−28222)︂ (︂)︂(︂ ,=·ΔΔ⎪00⎪0, = .⎩exp −22(2.51)48Когда оптическая плотность плазмы в отдельной ячейке достаточно высока(0 Δ ≫ 1), можно получить переход к асимптотической форме: (0 Δ) = √1,0 Δ1(0 Δ) = √ 3.8 0 Δ7(2.52)Используя безразмерный параметр ˜, = , /Δ, коэффициенты можно записать,⎧⎨−0.0313 · (˜, )−7/2 , ≠ ;=√0 Δ ⎩1, = .(2.53)Таким образом, в простейшем случае, полученные выражения сводятся к независимому от величины коэффициента поглощения произведению эффективнойвероятности eff на матрицу переноса, что подтверждает корректность метода.Представленный метод достаточно тривиален с точки зрения математики.

Однако, задача требует выполнения большого числа простых операций - возникаетнеобходимость в применении техник высокопроизводительных вычислений (HighPerformance Computing).49⇒Рисунок 2.11: Дискретизация области расчета путем вокселизации.2.3.3Алгоритм численного решения задачи1. Подготовка трехмерной области источника плазмы, в которой будет решаться задача.

Вычисляется объём исходной непрерывной области , затем осуществляется дискретизация (Рис. 2.11). Для дискретизации на однороднойдекартовой сетке был использована программа для вокселизации, разработанная Моррисом и Солсбери [174]. Пространственное разрешение сеткихарактеризуется максимальным числом вокселей вдоль оси - dim .2. Составляется вектор всех возможных пар элементов и . Каждая пара элементов содержит информацию о координатах (,,) и (′ , ′ , ′ ). Задаютсявходные параметры, такие, как , , , а также 0 - единое значение в случаеоднородного поглощения, или вектор значений при наличии неоднородности.3.

Характеристики

Список файлов диссертации

Свежие статьи
Популярно сейчас
Зачем заказывать выполнение своего задания, если оно уже было выполнено много много раз? Его можно просто купить или даже скачать бесплатно на СтудИзбе. Найдите нужный учебный материал у нас!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6363
Авторов
на СтудИзбе
310
Средний доход
с одного платного файла
Обучение Подробнее