Xreferat.com » Рефераты по физике » Кинетические уравнения Власова

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

Размещено на /


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

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

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


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

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

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

Певнев В.Я.


2011

Утверждаю

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

математики

Певнев В.Я.

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


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


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

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

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

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

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

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

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

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


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

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

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

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 Уравнение Власова-Максвелла, Власова-Эйнштейна и Власова-Пуассона


Вторую главу диплома хотелось бы посвятить непосредственно выводу или/и обоснованию системы уравнений Власова-Максвелла. Эта система уравнений выписана А.А. Власовым в работах1, и широко используется для описания плазмы. Уравнения Власова—Эйнштейна обосновываются аналогично, и я только коротко остановлюсь на них. Под названием уравнений Власова-Максвелла разные исследователи понимают разные уравнения. Наиболее популярно уравнение с нерелятивистской зависимостью скорости от импульса для функции распределения. Важно связать это уравнение с классическим лагранжианом, чтобы, с одной стороны, надежно иметь «правильное» уравнение, а с другой — понимать характер сделанных приближений. Далее, при выводе уравнения Власова-Максвелла будет приведён кратчайший, видимо, путь, связывая с лагранжианом электромагнетизма. Т.к. процесс вывода уравнения Власова-Максвелла является неоднозначным, то перед этим необходимо представить вспомогательные пункты. В 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).

Этим свойством иногда распоряжаются так, чтобы максимально упростить уравнения. Выберем1 в качестве параметра Кинетические уравнения Власовадлину линии (интервал, собственное время) 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 . Переход к переменным координата-импульс и гамильтонов формализм.

Введем стандартным образом импульсы.1 Если 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. Для простоты уравнений брали интервал, который в теории относительности называется собственным временем1. Возможность выбрать s в качестве параметра означает синхронизацию собственного времени различных частиц. С этим связан парадокс близнецов. Тот из них, чей интервал (собственное время) меньше, т.е. который «двигался больше», оказывается младше. Поэтому использование s хотя формально и возможно, но делает затруднительным интерпретацию результатов.


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


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


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


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

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

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


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


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


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

Кинетические уравнения Власова (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 , а где

Если Вам нужна помощь с академической работой (курсовая, контрольная, диплом, реферат и т.д.), обратитесь к нашим специалистам. Более 90000 специалистов готовы Вам помочь.
Бесплатные корректировки и доработки. Бесплатная оценка стоимости работы.

Поможем написать работу на аналогичную тему

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

Похожие рефераты: