Дипломная работа на тему "Кинетические уравнения Власова"

ГлавнаяФизика → Кинетические уравнения Власова




Не нашли то, что вам нужно?
Посмотрите вашу тему в базе готовых дипломных и курсовых работ:

(Результаты откроются в новом окне)

Текст дипломной работы "Кинетические уравнения Власова":


Дипломная робота

Пояснительная записка

«Кинетические уравнения Власова»

Студент группы Иванов И. И.

Руководитель работы Пересечанский В. М.

Заведующий кафедры "Математики"

Певнев В. Я.

2011

Утверждаю

Заведующий кафедрой

математики

Певнев В. Я.

"02" февраля 2011 г.

Задание на дипломную роботу

студенту Иванову Ивану Ивановичу пятого курса

1. Тема роботы: «Кинетические уравнения Власова»

Утверждена приказом

2. Срок сдачи студентом оконченной работы

3. Содержание пояснительной записки (перечень вопросов, которые подлежат рассмотрению): рассмотреть общие понятия кинетических уравнений, рассмотреть и вывести кинетические уравнения Власова, решить и описать одномерную модельную задачу для уравнения Власова

Заказать написание дипломной - rosdiplomnaya.com

Актуальный банк готовых успешно сданных дипломных проектов предлагает вам скачать любые проекты по необходимой вам теме. Безупречное написание дипломных работ по индивидуальным требованиям в Ростове-на-Дону и в других городах России.

Дата выдачи задания

Руководитель работы Пересечанский В. М.

Задание к выполнению принял

Согласовано Утверждаю

Руководитель дипломной работы Заведующий кафедрой

Пересечанский В. М. Певнев В. Я.

2011г. 2011г.

Календарный план дипломной работы

студента Иванова Ивана Ивановича

тема «Кинетические уравнения Власова»

--------------------------------------------------
Содержание работы | Срок исполнения (дата) |

Отметка о выполнении

(дата)

|
---------------------------------------------------------

1. Изучение литературы.

2. Анализ выбранной темы.

3. Обоснование актуальности темы.

4. Вопросы специального 1-го раздела

_-“- 2-го раздела

_-“- 3-го раздела

5. Устранение замечаний консультантов и руководителя.

6. Оформление пояснительной записки.

7. Предоставление работы на кафедру.

8. Предоставление работы на рецензию.

9. Предоставление работы на защитыв ГЭК

|

25.02

10.03

15.03

29.03

05.05

20.05

30.05

11.06

12.06

23.06

30.06

|

Выполнено.

Выполнено.

Выполнено.

Выполнено.

Выполнено.

Выполнено.

Выполнено.

Выполнено.

Выполнено.

Выполнено.

Выполнено.

|
--------------------------------------------------------- --------------------------------------------------

Студент группы: Иванов И. И.

2011 г.

План

Перечень условных сокращений и аббревиатур

Введение

Глава 1 Кинетические уравнения: основные понятия

1.1 Кинетические уравнения типа Больцмана

1.2 Уравнения типа Власова

Глава 2 Уравнение Власова-Максвелла, Власова-Эйнштейна и Власова-Пуассона

2.1 Сдвиг плотности вдоль траекторий динамической системы

2.2 Уравнения геодезических и эволюция функции распределения на римановом многообразии

2.3 Как ведет себя мера риманова пространства при преобразованиях

2.4 Вывод уравнения Власова-Максвелла

2.5 Схема вывода уравнения Власова-Эйнштейна

2.6 Система уравнений Власова-Пуассона для плазмы и электронов

Глава 3 Одномерная модельная задача для уравнения Власова

3.1 Условия

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

3.3 Математическая формализация задачи

3.4 Алгоритм разложения решения системы по параметру ε

3.5 Операторы Власова порядка n

3.6 Общая формула для поправки к полю порядка n

3.7 Классическое и релятивистское решения уравнения Власова

Заключение

Список литературы

Перечень условных сокращений и аббревиатур

ЭМП - Электромагнитное поле

Введение

Кинетические уравнения описывают эволюцию функции распределения F(t, v.x) молекул или других объектов (электронов, ионов, звезд, галактик или галактических скоплений) по скоростям v и пространству х в момент времени t. Это означает, что число частиц в элементе фазового объема dvdx есть F (t, v, x) dvdx.

Простейшее уравнение — уравнение свободного движения:

Рисунок убран из работы и доступен только в оригинальном файле. (1.1)

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

Глава 1 Кинетические уравнения: основные понятия

1.1 Кинетические уравнения типа Больцмана

Первым изученным кинетическим уравнением было уравнение Больцмана. Оно учитывает процессы столкновений добавлением интеграла столкновений в (1.1):

Рисунок убран из работы и доступен только в оригинальном файле. (1.2)

Интеграл столкновений J[F, F] — это квадратичный оператор, учитывающий парные столкновения частиц. Уравнение (1.2) было получено Максвеллом и Больцманом для вывода максвелловского распределения по скоростям, которое тогда только что было использовано для объяснения закона Менделеева – Клапейрона, который будет кратко рассмотрен далее.

Максвелловское распределение связано с одним из первых успехов уравнения Больцмана (1.2) — доказательством Н - теоремы.

Теорема утверждает, что функционал

Рисунок убран из работы и доступен только в оригинальном файле.

для уравнения Больцмана не возрастает: dH/dt <= 0. Этот факт был интерпретирован Больцманом как доказательство возрастания энтропии (Н есть энтропия с обратным знаком), т. е. обоснования 2-го закона термодинамики.

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

Однако уравнение Больцмана писалось Максвеллом для более широких целей. Программа Максвелла состояла в том, чтобы получить уравнения сплошной среды — типа уравнений Навье-Стокса — из уравнения Больцмана и тем самым получить коэффициенты переноса — вязкости и теплопроводности — и их зависимость от межмолекулярного взаимодействия. Ему это удалось для потенциала межмолекулярного взаимодействия U(r) = r -4 (максвелловские молекулы), когда интеграл столкновений сильно упрощается. Достичь аналогичных результатов для других потенциалов не удалось ни Больцману[1], ни Гильберту. однако это сделали Чэпмсн и Энског[2] с помощью специальной схемы теории возмущений (метод Чэпмсна-Энскога). Ставки здесь были очень высоки; такое решение давало бы (и дало: оно предсказало термодиффузию) количественные предсказания в молекулярно-кинетической теории, которая в то время подвергалась критике (в полемику включились не только ученые, например Мах и Авенариус, но и политики, например В. И. Ленин «Материализм и эмпириокритицизм». Чэпмсн и Энског «немного опоздали»: определение разными независимыми способами числа Авогадро с близкими ответами убедило ученых, и страсти улеглись.

В наше время это уравнение со своими следствиями работает в нескольких направлениях. Одно из них — средние слои атмосферы. Высокие слои хорошо описываются уравнением свободного движения (1.1) — газ Кнудсена или свободный газ. Низкие слои — уравнениями газодинамики, которые выводятся из уравнения Больцмана. Сопряжение хотя бы на ЭВМ верхних и низких слоев атмосферы — одна из актуальных задач[3] в связи с летательными аппаратами. Другое направление — химическая кинетика: моделирование смесей. Со всем этим связаны дискретные модели уравнения Больцмана

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

Предельным случаем уравнения Больцмана служит уравнение Ландау, когда наибольший вклад вносит сильное рассеянье вперед. Оно используется для описания плазмы.

Используются также квантовые аналоги уравнения Больцмана — уравнения Улинга-Уленбека. Для этих уравнений стационарными распределениями вместо максвелловского оказываются распределения Ферми-Дирака или Бозе-Эйнштейна.

Таким образом, можно представить иерархию уравнений типа Больцмана в виде следующей схемы:

Рисунок убран из работы и доступен только в оригинальном файле.

Схема 1

Линии с вопросительными знаками означают, что соответствующие уравнения еще, может быть, не выведены (например, приближение Ландау для уравнения Улинга-Уленбека).

1.2 Уравнения типа Власова

Если уравнения типа Больцмана описывают короткодействующие взаимодействия, то уравнения типа Власова описывают дальнодействие.

Уравнения Власова или уравнения самосогласованного поля имеют вид:

Рисунок убран из работы и доступен только в оригинальном файле. (2.1)

Здесь сила f сама есть функционал от функции распределения, а уравнение (2.1) имеет вид уравнения сдвига вдоль характеристик. Простейший вид зависимости силы f от функции распределения соответствует парному потенциалу взаимодействия К(х, у):

Рисунок убран из работы и доступен только в оригинальном файле. (2.2)

Такой вид взаимодействия дает систему уравнений Власова. Обычно говорят о системах уравнений «Власова плюс ещё кого-то» для того, чтобы различать виды взаимодействий. Бывают уравнения Власова-Пуассона, Власова-Максвелла, Власова-Эйнштейна и Власова--Янга--Миллса.

Уравнение Власова-Пуассона бывает двух видов — для гравитации и для плазмы: в обоих случаях (2.2) заменяется на уравнение Пуассона действием оператора Лапласа, при условии, что К (х, у) — фундаментальное решение оператора Лапласа. Таким образом, К сеть потенциал единичного заряда в трехмерном случае, нити — в одномерном случае и плоскости — в двумерном.

Если в гравитационном случае мы заменяем взаимодействие по Ньютону на взаимодействие по Эйнштейну, то получаем уравнение Власова-Эйнштейна.

Если в случае плазмы мы заменяем электростатику на электродинамику, то получаем уравнения Власова-Максвелла. Если у нас сохраняется не заряд, а векторная величина (изотопический заряд или цвет), то вместо электромагнитных 4-иотснциалов мы должны взять матрицы, и получаем уравнения Янга-Миллса. Такие уравнения дают принятую в настоящее время теорию объединенного электрослабого и сильного взаимодействия. Таким образом, все уравнения чипа Власова дают следующую иерархию:

Рисунок убран из работы и доступен только в оригинальном файле.

Схема 2

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

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

Энергетическая подстановка, когда функция распределения зависит только от энергии. В этом случае уравнение (2.1) удовлетворяется, а (2.2) переходит в нелинейное уравнение для потенциала. Это уравнение аналогично уравнениям Бернулли для уравнения Эйлера. И уравнения типа Власова по своей судьбе аналогичны уравнениям Эйлера: их частные случаи стали появляться раньше, чем были написаны сами уравнения Власова. При этом в той же самой энергетической подстановке, выражающей закон сохранения энергии. В приложениях это были плазменный диод (диод Ленгмюра), уравнение Дебая для электролитов и уравнение Лэна-Эмдена в гравитации. В математике такое уравнение еще раньше было изучено в геометрии и называется уравнением Лиувилля. В двумерном случае оно имеет огромную группу симметрии (конформная группа).

Глава 2 Уравнение Власова-Максвелла, Власова-Эйнштейна и Власова-Пуассона

Вторую главу диплома хотелось бы посвятить непосредственно выводу или/и обоснованию системы уравнений Власова-Максвелла. Эта системауравнений выписана А. А. Власовым в работах[4], и широко используется для описания плазмы. Уравнения Власова—Эйнштейна обосновываются аналогично, и я только коротко остановлюсь на них. Под названием уравнений Власова-Максвелла разные исследователи понимают разные уравнения. Наиболее популярно уравнение с нерелятивистской зависимостью скорости от импульса для функции распределения. Важно связать это уравнение с классическим лагранжианом, чтобы, с одной стороны, надежно иметь «правильное» уравнение, а с другой — понимать характер сделанных приближений. Далее, при выводе уравнения Власова-Максвелла будет приведён кратчайший, видимо, путь, связывая с лагранжианом электромагнетизма. Т. к. процесс вывода уравнения Власова-Максвелла является неоднозначным, то перед этим необходимо представить вспомогательные пункты. В 2.1 будет рассмотрено как обосновываются уравнения для функции распределения частиц, сдвигаемых вдоль траекторий произвольной динамической системы хi = Xi(x). А далее изучается уравнение Эйлера-Лафанжа для случая, когда действие есть длина, а также обосновывается выбор функции распределения в переменных х, р (пространство-импульсы).

2.1 Сдвиг плотности вдоль траекторий динамической системы

Рассмотрим произвольную динамическую систему, т. е. систему нелинейных дифференциальных уравнений в k-мерном пространстве:

хi = Xi(x), i = 1,…,k (1.1)

Пусть мы раскидали частицы с какой-то начальной плотностью f(0, x), а в момент времени t эта плотность есть f(t, x), так что число частиц в области G

Рисунок убран из работы и доступен только в оригинальном файле.

Какова эволюция f(t, x)?

Покажем, что соответствующее уравнение имеет вид (по повторяющимся верхним и нижним индексам предполагается суммирование):

Рисунок убран из работы и доступен только в оригинальном файле. (1.2)

Способ 1 . Метод Рисунок убран из работы и доступен только в оригинальном файле.-функций.

Рассмотрим функцию распределения N частиц, сдвигающихся по траекториям этой системы:

Рисунок убран из работы и доступен только в оригинальном файле.

где для каждого l функция xi(t) удовлетворяет уравнениям (1.1). Тогда, дифференцируя по времени, получаем

Рисунок убран из работы и доступен только в оригинальном файле.

С другой стороны, имеем

Рисунок убран из работы и доступен только в оригинальном файле.

Складывая полученные выражения, находим, что уравнения (1.2) для такой функции удовлетворяются. При взятии дивергенции воспользовались формулой

Рисунок убран из работы и доступен только в оригинальном файле.

Для произвольной функции f равенство (1.2) получается переходом к пределу при аппроксимации ее суммой Рисунок убран из работы и доступен только в оригинальном файле.-функций (в слабом смысле).

Способ 2 . Баланс частиц.

Скорость роста частиц в области G есть

Рисунок убран из работы и доступен только в оригинальном файле. (1.3)

Это следует из того, что на малом участке границы ds количество вылетевших за время dt частиц есть f ds dt(X, n), так как все вылетевшие частицы заметают цилиндр с основанием ds и стороной X dt, а поэтому высотой (X, п) dt. Знак минус берется потому, что нормаль — внешняя, и считаются вылетевшие частицы, тогда как слева в (1.3) стоит скорость роста числа частиц в области G. Преобразуя в правой части (1.3) интеграл из поверхностного в объемный по формуле Стокса, мы получим уравнение (1.2), проинтегрированное по области G, а отсюда в силу произвольности G — и само уравнение (1.2).

Перепишем уравнение (1.2) в виде

Рисунок убран из работы и доступен только в оригинальном файле. (1.4)

Если divX = 0, то левая часть (1.4) — это полная производная f(t, x) пo времени.

Вывод. Уравнение для функции распределения частиц, сдвигающихся вдоль траекторий динамической системы (1.1), имеет вид (1.2).

2.2 Уравнения геодезических и эволюция функции распределения на римановом многообразии

Рассмотрим метрику gijdxidxj в пространстве Rn, xРисунок убран из работы и доступен только в оригинальном файле.Rn, gij(x)-n2 функций. Это означает, что длина кривой определяется формулой:

Рисунок убран из работы и доступен только в оригинальном файле. (2.1)

а уравнение геодезических получается из принципа наименьшего действия (принципа наименьшей длины). Если, более обще, действие записывается в виде S = Рисунок убран из работы и доступен только в оригинальном файле.dt, где L—лагранжиан, то уравнения Эйлера-Лагранжа даются варьированием с фиксированными концами траекторий:

Рисунок убран из работы и доступен только в оригинальном файле.

Получаем уравнения Эйлсра-Лагранжа:

Рисунок убран из работы и доступен только в оригинальном файле.

В случае геодезических L = Рисунок убран из работы и доступен только в оригинальном файле. имеем

Рисунок убран из работы и доступен только в оригинальном файле. (2.2)

Функционал длины инвариантен относительно замены t = Рисунок убран из работы и доступен только в оригинальном файле. для любой гладкой функции Рисунок убран из работы и доступен только в оригинальном файле., и то же свойство имеют уравнения (2.2).

Этим свойством иногда распоряжаются так, чтобы максимально упростить уравнения. Выберем[5] в качестве параметра Рисунок убран из работы и доступен только в оригинальном файле.длину линии (интервал, собственное время) s : ds =Рисунок убран из работы и доступен только в оригинальном файле., после деления на ds получим Рисунок убран из работы и доступен только в оригинальном файле. = 1, и уравнения (2.2) превращаются в

Рисунок убран из работы и доступен только в оригинальном файле. (2.3)

Последние совпадают с уравнениями Эйлера-Лагранжа для действия с лагранжианом

Рисунок убран из работы и доступен только в оригинальном файле.

Преобразуем их к виду

Рисунок убран из работы и доступен только в оригинальном файле. (2.4)

Здесь gki — матрица, обратная gij а Рисунок убран из работы и доступен только в оригинальном файле. называются символами Кристоффеля.

Запишем уравнения (1.2) для функции распределения f(х, v, s) по пространству и скоростям (с длиной. s вместо времени) для уравнения (2.4), как это показано в предыдущем параграфе:

Рисунок убран из работы и доступен только в оригинальном файле. (2.5)

Последний член в левой части соответствует тому, что система (2.4) имеет дивергенцию, отличную от нуля. Переход к бездивергентному виду можно осуществить двумя способами.

Способ 1 . Переход к переменным координата-импульс и гамильтонов формализм.

Введем стандартным образом импульсы.[6] Если L = (gijxixj)/2 (этот лагранжиан дает те же уравнения движения, что и (2.1)), то импульсы pi = = dL/dxi= gijxj, а гамильтониан Н = pivi — L = (pipgij)/2. Тогда уравнения (2.3) приобретают гамильтонов вид

Рисунок убран из работы и доступен только в оригинальном файле.

Упражнение

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

Решение

Рисунок убран из работы и доступен только в оригинальном файле.

Получаем уравнения для функции распределения f(s, х,р) по координатам и импульсам (1.2) в виде

Рисунок убран из работы и доступен только в оригинальном файле. (2.6)

Эго уравнение имеет вид

df/ds + {Н, f} = 0, где {Н, f} — скобка Пуассона:

Рисунок убран из работы и доступен только в оригинальном файле.

Способ 2. Переход к инвариантной мере в пространстве координаты-скорости.

Пусть g — определитель матрицы gij. Вместо f в (2.5) введем новую функцию распределения

F{x, v,s) = F(x, v,s)/g.

Упражнение

Показать, что для новой функции распределения уравнение эволюции бездивергентно и имеет вид

Рисунок убран из работы и доступен только в оригинальном файле.

Решение

Воспользуемся операцией дифференцирования определителя. При этом второе слагаемое в (2.5) преобразуется следующим образом:

Рисунок убран из работы и доступен только в оригинальном файле.

В (a) используется тождество

Рисунок убран из работы и доступен только в оригинальном файле.

Для новой функции распределения число частиц записывается в виде

Рисунок убран из работы и доступен только в оригинальном файле.

Поэтому g dxdv есть инвариантная мера: F не растет, т. е. полная производная от неё есть ноль, и поскольку число частиц сохраняется, то мера g dxdv сохраняется тоже.

Вывод. В качестве переменных в функции распределения можно брать импульсы или скорости, а в качестве времени — время или интервал s. Для простоты уравнений брали интервал, который в теории относительности называется собственным временем[7]. Возможность выбрать s в качестве параметра означает синхронизацию собственного времени различных частиц. С этим связан парадокс близнецов. Тот из них, чей интервал (собственное время) меньше, т. е. который «двигался больше», оказывается младше. Поэтому использование s хотя формально и возможно, но делает затруднительным интерпретацию результатов.

2.3 Как ведет себя мера риманова пространства при преобразованиях

Пусть проведена замена координат хк = f (Рисунок убран из работы и доступен только в оригинальном файле.). Как преобразуется при этом метрика? Имеем:

Рисунок убран из работы и доступен только в оригинальном файле.

Поэтому Рисунок убран из работы и доступен только в оригинальном файле. где J — это det (дxi/дРисунок убран из работы и доступен только в оригинальном файле.), Отсюда следует, что так как dx = |J|dРисунок убран из работы и доступен только в оригинальном файле., то Рисунок убран из работы и доступен только в оригинальном файле. = Рисунок убран из работы и доступен только в оригинальном файле., т. е. Рисунок убран из работы и доступен только в оригинальном файле.—инвариант преобразований.

Дифференцируя по параметру, имеем Рисунок убран из работы и доступен только в оригинальном файле., а поэтому dV=|J|dv. Отсюда следует, что g dxdv = Рисунок убран из работы и доступен только в оригинальном файле.Рисунок убран из работы и доступен только в оригинальном файле. - инвариантная мера, где каждый из сомножителей инвариантен при преобразованиях.

Вывод. В качестве переменных функции распределения удобно брать импульсы. В качестве параметра Рисунок убран из работы и доступен только в оригинальном файле. возьмем время, в качестве переменных функции распределения — t (время), х (пространственная координата), р(импульсы): f= f(t, x,p).

2.4 Вывод уравнения Власова-Максвелла

Система уравнений Власова-Максвелла описывает движение частиц в собственном электромагнитном поле. Стартуем с обычного действия для электромагнитного поля[8], действия Власова-Максвелла или Лоренца (по повторяющимся верхним и нижним индексам идет суммирование):

Рисунок убран из работы и доступен только в оригинальном файле.

Рисунок убран из работы и доступен только в оригинальном файле. (4.1)

Рисунок убран из работы и доступен только в оригинальном файле.

где Sр означает действие частиц (particles), Sf — действие полей (fields), Sp-f — действие частиц-полей (particles-fields).

Здесь а означает сорт частиц, отличаемый по массе mа и заряду еa, q нумерует частицы внутри сорта, Рисунок убран из работы и доступен только в оригинальном файле.(q. t) (Рисунок убран из работы и доступен только в оригинальном файле. = 0,1.2,3; q =1,...,Na; a=1..... r) — 4 координаты q-й частицы copтa a, Au(x) — потенциал, Рисунок убран из работы и доступен только в оригинальном файле. — электромагнитные поля, Рисунок убран из работы и доступен только в оригинальном файле. - метрика Минковского: Рисунок убран из работы и доступен только в оригинальном файле., т. е. диагональная матрица с 1 на первом месте и (-1) на остальных. Варьирование проводим специальным способом: сначала получаем движение частицы в поле, потом поля с заданными движениями частиц. Однако для частиц мы перейдем к функциям распределения, что и даст искомую систему уравнений.

1. Варьирование Sp + Sp+f по координатам Рисунок убран из работы и доступен только в оригинальном файле.(q. t)) даст уравнение движения зарядов в поле. Перепишем Рисунок убран из работы и доступен только в оригинальном файле. для метрики Минковского (в дальнейшем греческие индексы Рисунок убран из работы и доступен только в оригинальном файле., Рисунок убран из работы и доступен только в оригинальном файле. пробегают четыре значения: Рисунок убран из работы и доступен только в оригинальном файле. = 0,1,2,3; латинские i, j —три: i = 1,2,3):

Рисунок убран из работы и доступен только в оригинальном файле.

где Lp, — лагранжиан частиц.

Здесь Рисунок убран из работы и доступен только в оригинальном файле. (i = 1,2,3) — трехмерный квадрат скорости, и мы учли, что х° = ct и вынесли с2 из-под корня. Проварьируем это выражение (опуская а):

Рисунок убран из работы и доступен только в оригинальном файле.

Варьируем Sp-f (снова опускаем а):

Рисунок убран из работы и доступен только в оригинальном файле.

Рисунок убран из работы и доступен только в оригинальном файле.

Отсюда из условия Рисунок убран из работы и доступен только в оригинальном файле. = 0 получаем уравнение движения заряженной частицы в поле:

уравнение больцман власов динамический модельный

Рисунок убран из работы и доступен только в оригинальном файле.

где Рисунок убран из работы и доступен только в оригинальном файле.

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

Рисунок убран из работы и доступен только в оригинальном файле.

Обозначая Рисунок убран из работы и доступен только в оригинальном файле. получаем Рисунок убран из работы и доступен только в оригинальном файле. = Рисунок убран из работы и доступен только в оригинальном файле.Отсюда находим уравнение для функции распределения fa(x, p,t) (аналог 1.4):

Рисунок убран из работы и доступен только в оригинальном файле. (4.2)

Здесь Рисунок убран из работы и доступен только в оригинальном файле. Использовано, что Рисунок убран из работы и доступен только в оригинальном файле.

В это уравнение записано для ионов и электронов в следующем виде:

Рисунок убран из работы и доступен только в оригинальном файле. (4.3)

Здесь fi(t, р, х) — функция распределения ионов по пространству и импульсам в момент времени t (i в (4.3) — первая буква слова ion. а не индекс), fе(t, р, х) — функция распределения электронов, ze — заряд иона, (—е) — заряд электрона, [v, B] — векторное произведение. Не выписано выражение v через р, однако часто его берут классическим: vаj = pj/ma, и тогда удобно записать уравнения через функцию распределения f(t, v, х) по скоростям вместо импульса. В записи (4.3) v надо брать различными для электронов и ионов, т. е. (4.3) требует уточнения, где vi, а где vc вместо v, и каковы эти функции, как функции импульса vi(p) и vc(p).

3. Уравнение для полей. Используем функцию распределения вместо плотности. Сначала надо переписать Sp-f через функцию распределения, совершив переход

Рисунок убран из работы и доступен только в оригинальном файле.

после чего Sp-f запишется в виде

Рисунок убран из работы и доступен только в оригинальном файле.

Теперь варьируем по потенциалам Аu(х):

Рисунок убран из работы и доступен только в оригинальном файле.

Рисунок убран из работы и доступен только в оригинальном файле.

Полагаем Рисунок убран из работы и доступен только в оригинальном файле. и получаем

Рисунок убран из работы и доступен только в оригинальном файле. (4.4)

Система (4.2), (4.4) и есть система уравнений Власова-Максвелла.

Замечание 1. Уравнения (4.4) являются второй парой уравнений Максвелла, а первая следует из равенств Рисунок убран из работы и доступен только в оригинальном файле. что записывается в эквивалентном виде на языке дифференцирования кососимметрических тензоров Рисунок убран из работы и доступен только в оригинальном файле. Первая пара уравнений Максвелла записывается в виде Рисунок убран из работы и доступен только в оригинальном файле.

Замечание 2. При выводе уравнений Власова-Максвелла по схеме Боголюбова мы должны были бы стартовать с гамильтоновых систем с потенциалами Лиенарта-Вихерта (запаздывающие потенциалы). Для слабого релятивизма соответствующий лагранжиан называется лагранжианом Дарвина.

Замечание 3. Можно таким же способом получить уравнения Вла-сова-Янга-Миллса, заменив в четырех потенциалах Au числа на матрицы.

Вывод, Система уравнений Власова-Максвелла получается при варьировании действия электромагнетизма (действия Лоренца) с переходом к функции распределения. Уравнение для функции распределения получается как уравнение сдвига вдоль траекторий движения частиц.

2.5 Схема вывода уравнения Власова-Эйнштейна

Рассмотрим действие для частицы в гравитационном поле и для поля[9]:

Рисунок убран из работы и доступен только в оригинальном файле. (5.1)

Здесь R — кривизна; вариация но метрике производится переписыванием первого слагаемого в эйлеровых координатах:

Рисунок убран из работы и доступен только в оригинальном файле. (5.2)

При этом, как и в предыдущем параграфе, получается уравнение для полей. Варьируя траектории частиц в Sp в (5.1), получаем уравнение для их движения в гравитационном поле. Уравнение для функции распределения (как и в предыдущем параграфе) сеть уравнение сдвига вдоль характеристик.

Вывод. Уравнения Власова-Максвелла и Власова-Эйнштейна получаются единообразным способом варьирования соответствующих лагранжианов электромагнетизма и гравитации.

2.6 Система уравнений Власова-Пуассона для плазмы и электронов

Рассмотрим систему уравнений Власова-Максвелла в потенциалах Аu.. Получаем волновую форму релятивистской системы уравнений Власова-Максвелла при условии дuАu = 0 (лоренцова калибровка)[10]:

Рисунок убран из работы и доступен только в оригинальном файле.

Рисунок убран из работы и доступен только в оригинальном файле. (6.1)

Рисунок убран из работы и доступен только в оригинальном файле.

Здесь уравнения Максвелла преобразованы по

Рисунок убран из работы и доступен только в оригинальном файле.

Если мы перейдем к нерелятивистскому пределу, то va = p/ma. Обычно рассматривают функцию f(t, x,v) распределения по скоростям. Получим систему уравнений Власова-Пуассона:

Рисунок убран из работы и доступен только в оригинальном файле.

Рисунок убран из работы и доступен только в оригинальном файле. (6.2)

Рисунок убран из работы и доступен только в оригинальном файле.

Как правило п = 1 (электронная задача) или п = 2 (плазма, состоящая из ионов и электронов).

Глава 3 Одномерная модельная задача для уравнения Власова

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

3.1 Условия

Пусть радиационная генерация электромагнитного поля (далее – ЭМП) идёт в широком потоке быстрых электронов, образующемся при рассеянии в среде ионизирующих частиц. В общем случае ЭМП будет самосогласованным, то есть заметно влиять на динамику быстрых электронов, а решение подобных трёхмерных задач в сложных средах требует зачастую недоступного в настоящее время объема ресурсов ЭВМ. В ряде практически важных случаев параметры потоков ионизирующих частиц и среды оказываются такими, что эффект согласования ЭМП и плотности тока быстрых электронов невелик. Его можно попытаться учесть как малое возмущение при определении значений плотности тока и ЭМП. Другие моменты функции распределения быстрых электронов в задачах радиационной генерации ЭМП как правило интереса не представляют. Этот подход к трёхмерным задачам целесообразно исследовать в рамках более простых дву - и даже одномерных постановок при определении необходимого набора функционалов от функции распределения[11].

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

Целью практической части работы является получение асимптотического разложения решения задачи о торможении потока электронов собственным электрическим полем в одномерной вакуумной модели и преобразование его к виду, который интерпретируется в терминах рядов Бюрмана-Лагранжа. При этом физически очевидно, что влияние поля на электроны в такой задаче заведомо завышено, что гарантирует применимость выводов о сходимости получающихся разложений к решениям для реальных моделей.

Сформулированная проблема имеет универсальный характер: анализ некоторых задач с реакторной тематикой – управление работой реактора, детальный учет эффекта выгорания и др. – приводит к необходимости рассматривать аналогичную ситуацию. Для описания поведения реактора ищется его нейтронное поле N = N(t, , ). Такая информация, разумеется, избыточна, но с помощью неё сравнительно просто находить различные функционалы от уже найденного поля, значения которых и представляют практический интерес. В данном случае интерес представляет величина ЭМП – линейного функционала от электронного поля f = f(t, , ), – которое само входит в систему Власова-Максвелла и, как будет рассмотрено далее, может в некоторых достаточно простых, но имеющих важные следствия, случаях быть найдено, минуя определение функции f(t, , ), нахождение которой в дальнейшем сводится к решению классических задач с линейными уравнениями первого порядка. В этом состоит определённая специфика системы Власова-Максвелла.

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

Модели ЭМП в задачах его радиационной генерации строятся на основе уравнений Максвелла, содержащих роторы электрического и магнитного полей. Это связано с тем, что быстрые электроны сами ионизуют среду, порождая вторичные электроны и ионы низкой энергии. Концентрация вторичных заряженных частиц во много раз превышает концентрацию быстрых электронов. По этой причине вторичные электроны в математических моделях рассматриваются отдельно. Источником ЭМП является сторонний ток быстрых электронов и проводимость слабоионизованного газа вторичных заряженных частиц. Соответствие такой модели закону Гаусса обеспечивается уравнением непрерывности для системы всех заряженных частиц.

Рассмотрим эмиссию монохроматического потока электронов со всей плоскости xOy вдоль оси Oz в вакуум. Функция распределения электронов f(t, z, p) удовлетворяет одномерному уравнению Власова, записанному в импульсной, учитывающей релятивистские эффекты, форме:

Рисунок убран из работы и доступен только в оригинальном файле.

а входящее в него E – самосогласованное электрическое поле – закону Ампера:

Рисунок убран из работы и доступен только в оригинальном файле., Рисунок убран из работы и доступен только в оригинальном файле..

Здесь и далее: с – скорость света;Рисунок убран из работы и доступен только в оригинальном файле.Рисунок убран из работы и доступен только в оригинальном файле. – скорость электрона; m0, re, e, – его масса покоя, классический радиус и заряд соответственно. Плоский поток заряженных частиц порождает одну компоненту плотности электрического тока j(t, z): вдоль оси Oz. В условиях такой симметрии ненулевой является только Εz-компонента ЭМП. Остальные его компоненты равны нулю вследствие уравнений Максвелла и однородных начальных данных: f(0, z, p) = 0, Ε(0, z) = 0.

Функция F(t) описывает интенсивность эмиссии с плоскости, то есть число электронов, вылетающих с единицы площади поверхности в единицу времени. Как и в реальных моделях, для неё выполнено условие F(0) = 0.

3.3 Математическая формализация задачи

В дальнейшем все функции времени предполагаются, если не оговорено противное, продолженными нулем на интервал (t < 0), а функции от – на интервал (z < 0). Аргументы функции (или часть их) опускаются при записи, если это не приводит к недоразумению. Запись производной просто штрихом (или точкой) означает, что она взята по всему аргументу функции, а не по какой-либо его составляющей.

Задача, поставленная в предыдущем пункте, требует предварительного анализа, например, на предмет сокращения числа входящих в основные уравнения параметров – то есть, как минимум, приведения их к безразмерному виду. Проводя эту стандартную ( t= tτ, r= rρ, p= pπ, f = f ̃φ(τ, ρ, π), F = F̃Φ(τ, ρ, π), E(t, r) = EΕ(τ, ρ), J = JJ(τ, ρ), Sext = Qext/Q) операцию, видим, что между масштабными – взволнованными – коэффициентами должны иметь место стандартные же соотношения, дабы безразмерные уравнения не отличались по своей структуре от своих стартовых размерных аналогов. Таковыми являются: связь скорости и импульса (классическая или релятивистская); r = vt, причем, естественнее всегда брать c – за масштаб скоростей, а из t̃ и r̃ выбирать только одну.

В результате приходим к следующей системе соотношений, приводящих исходную систему к полностью безразмерному виду: один из параметров L = r̃, или T = t̃ является свободным, L = cT, P = p̃ = m0v0. Функции, входящие в систему, имеют своими масштабными коэффициентами следующие величины: Q̃ = 2N/(LTP), f ̃ = 2N/(LP), J̃= – |e|N/T, Ẽ = 4π|e|N, а единственный параметр, остающейся сомножителем перед Ε, – ε можно выразить как через начальные данные: ε = 4πreLN, так и через широко используемые (ωплаз)2 = 4πrec2n – плазменную частоту и νист = 1/T – частоту источника: ε = υ0(ωплаз/νист)2.

Приведенная таким образом к безразмерному виду исходная система приобретает вид:

Рисунок убран из работы и доступен только в оригинальном файле. (1)

Рисунок убран из работы и доступен только в оригинальном файле. (2)

Первой задачей анализа системы (1-2) будет получение явных формул для Рисунок убран из работы и доступен только в оригинальном файле.путём разложения его в ряд по степеням ε, что сведёт дальнейшее решение уравнения (1) к решению классического уравнения первого порядка.

3.4 Алгоритм разложения решения системы по параметру ε

Далее, на первом этапе исследования, при получении формул для E(t, z), нам потребуются производные всех порядков от временной компоненты источника F(t). Считаем, что она является действительной аналитической функцией. Зависимость v = v(p) полагаем аналитической по тем же причинам: как классическая, так и квантовая её модели, разумеется, этим свойством обладают, а для построения решения удобнее рассматривать сразу общий случай υ = υ(π) произвольного диффеоморфизма луча π > 0 на луч υ > 0 либо на интервал 0 < υ < υ0.

Поиск начального приближения Рисунок убран из работы и доступен только в оригинальном файле. уравнения (1) приводит к формулам:

Рисунок убран из работы и доступен только в оригинальном файле.

Далее, не обговаривая специально, удобно придерживаться следующих обозначений: χ = τ – ζ/υ, χ0 = τ – ζ/υ0, τ – (ζ – υ(τ – τ̃))/υ0 = χ̃0, χ̃ = χ.

Пусть Рисунок убран из работы и доступен только в оригинальном файле.. Разложив по степеням ε произведение εE φ’ и приравнивая, друг другу коэффициенты при всех последовательных степенях, получаем, как обычно, бесконечную серию уравнений, зацепленных каждое только заодно другое своими правыми частями – последовательными источниками частиц, испытавших данное число взаимодействий (соударений). Начальное уравнение цепочки (с S0 = Sext для φ0) уже выписано. Основным для дальнейшего будет то, что левая часть у всех последующих уравнений одинакова. Правые части их имеют следующий вид: Рисунок убран из работы и доступен только в оригинальном файле.. Тождественность операторов Рисунок убран из работы и доступен только в оригинальном файле., порождающих все уравнения, позволяет следующим образом записать их решения φn, в операторной форме: Рисунок убран из работы и доступен только в оригинальном файле., где Рисунок убран из работы и доступен только в оригинальном файле., а Рисунок убран из работы и доступен только в оригинальном файле. – это оператор сдвига по характеристике (невозмущённого) уравнения переноса: ζ → ζ – υ(τ – τn+1). Далее Εm – это оператор умножения на соответствующую функцию, а Рисунок убран из работы и доступен только в оригинальном файле.; таким образом в развёрнутой записи имеем соотношение Рисунок убран из работы и доступен только в оригинальном файле.. В последней формуле дифференцирования по dπ отмеченного υ = υ(π) НЕ производится.

На этом пути получаются весьма громоздкие явные выражения для поправок φn при малых n.

Рисунок убран из работы и доступен только в оригинальном файле.Из них для J1(τ, χ0) и Ε1(τ, χ0) получаются весьма простые выражения: Подчеркнем, что простота полученных формул есть следствие того, что для данной задачи оператор обращения уравнений Максвелла Рисунок убран из работы и доступен только в оригинальном файле. – это просто интегрирование по dτ от 0 до τ. В результате и все поправки высших порядков выразятся как полиномы от τ с коэффициентами, зависящими только от Ε0(χ0) и её производных. Формула для Ε1(τ,χ0) уже выписана,

Ε2(χ0)=Рисунок убран из работы и доступен только в оригинальном файле., а Ε3(χ0)=

Рисунок убран из работы и доступен только в оригинальном файле.

Далее естественно было предположить, что и общая формула для поправки к Ε(τ,χ0) порядка n будет иметь аналогичный вид:

Εn(τ,χ0) = Рисунок убран из работы и доступен только в оригинальном файле.,

где Рисунок убран из работы и доступен только в оригинальном файле. – полином степени k от Рисунок убран из работы и доступен только в оригинальном файле., на что указывает показатель степени υ0 в знаменателях его коэффициентов. То, что получатся именно полиномы, а не мономы, как при малых n, угадывается при анализе характера упрощений в полученных формулах при переходе от функции распределения к току электронов; уже при n = 4 в коэффициент при Рисунок убран из работы и доступен только в оригинальном файле., войдёт сумма A Рисунок убран из работы и доступен только в оригинальном файле.+ B Рисунок убран из работы и доступен только в оригинальном файле. с наперед неизвестными значениями A и B.

3.5 Операторы Власова порядка n

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

Записав с их помощью начальный отрезок разложения поля Ε в ряд по ε, получаем следующее выражение полного поля Ε через невозмущённое поле Ε0:

Рисунок убран из работы и доступен только в оригинальном файле.

Связь скобок – а в них соответственно по 2(n-1) слагаемых – в данной формуле с введенными в предыдущем разделе величинами Εq очевидна, но их внутренняя структура нетривиальна. Например, последнюю из них – Ε3 (и, как очевидно, все последующие) разбить на q компонент можно как минимум двумя способами с различными смысловыми интерпретациями слагаемых.

Во-первых, Εq = Рисунок убран из работы и доступен только в оригинальном файле., где индекс m соответствует номеру φm от которого берётся производная по импульсу. В этом случае получаем:

Рисунок убран из работы и доступен только в оригинальном файле.

что означает различение взаимодействий частиц с полями различных порядков (индексов). Во-вторых, то же разложение пишется в симметричной форме:

Рисунок убран из работы и доступен только в оригинальном файле., (3)

где учитывается только общее число взаимодействий частиц с полем.

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

3.6 Общая формула для поправки к полю порядка n

Вторая форма записи разложения (3) позволит найти эту общую формулу. Рассмотрим операторный вид общей её формы:

Рисунок убран из работы и доступен только в оригинальном файле., (4)

и аналогично для φ, но без оператора Рисунок убран из работы и доступен только в оригинальном файле., решающего систему Максвелла.

Заметим, что в каждой внутренней сумме Рисунок убран из работы и доступен только в оригинальном файле. слагаемых, поскольку операторы Власова различных порядков не коммутируют. То есть, φn является суммой δ-функции и её n первых производных с соответствующими коэффициентами (в которые входят Εν для ν=1,…,(n - 1) ). Такое представление одной из компонент решения задачи – φ, очевидно, плохо тем же, чем плохо представление решения произвольного интегрального уравнения через определители Фредгольма: своей необозримостью. Для того чтобы избегнуть подобной ситуации, надо упростить формулу (4) так, чтобы далее привести задачу определения φ к стандартной.

При упрощении формулы (4) коэффициенты, при последовательных степенях τk, зависящие от υ и её производных по dπ оказываются удовлетворяющими рекуррентному соотношению: Рисунок убран из работы и доступен только в оригинальном файле. = Рисунок убран из работы и доступен только в оригинальном файле.Рисунок убран из работы и доступен только в оригинальном файле. + (n–1+k) Рисунок убран из работы и доступен только в оригинальном файле.Рисунок убран из работы и доступен только в оригинальном файле.. Здесь оператор Рисунок убран из работы и доступен только в оригинальном файле. = Рисунок убран из работы и доступен только в оригинальном файле.Рисунок убран из работы и доступен только в оригинальном файле.+ Рисунок убран из работы и доступен только в оригинальном файле.’’Рисунок убран из работы и доступен только в оригинальном файле. +… .

С учетом вычисленных заранее результатов по поправкам малых порядков это означает, что получившийся для описания поправки порядка n многочлен с символической переменной Рисунок убран из работы и доступен только в оригинальном файле. совпадает с полиномом Рисунок убран из работы и доступен только в оригинальном файле. для производных от обратных функций. То, что порядок его смещен на единицу, вызвано необходимостью снятия одного дифференцирования по dτ (в Рисунок убран из работы и доступен только в оригинальном файле.). По той же причине внутренняя его переменная – это не входящая в уравнение Власова (3) скорость υ, но интеграл от неё по импульсу, то есть энергия: ω = ω(υ). Поскольку упомянутые полиномы Рисунок убран из работы и доступен только в оригинальном файле. стандартно выражаются через имеющие явную формулу полиномы Белла, приходим в итоге к одному из вариантов формулы Бруно.

Суммируя полученные при выводе формул (3-4) результаты, приходим к завершающей формуле (5) для записи решения системы (1-2) в замкнутом виде. Здесь и далее, как и ранее Рисунок убран из работы и доступен 
<p>Здесь опубликована для ознакомления часть дипломной работы Если у вас нет возможности самостоятельно написать дипломную - закажите её написание опытному автору»


Просмотров: 524

Другие дипломные работы по специальности "Физика":

Электроснабжение завода продольно-строгальных станков

Смотреть работу >>

Математическое моделирование пластической деформации кристаллов

Смотреть работу >>

Электроснабжение фермы КРС на 800 голов в ОАО "Петелино" Ялуторовского района Тюменской области с обеспечением нормативных условий надежности

Смотреть работу >>

Электроснабжение судоремонтного завода

Смотреть работу >>

Повышение надежности электроснабжения потребителей н. п. Орлово Армизонского района Тюменской области с выбором оборудования на ПС 110/10 кВ "Орлово"

Смотреть работу >>