Динамика вещества в тесных двойных звездах Дмитрий Валерьевич Бисикало. - презентация

1

2 Динамика вещества в тесных двойных звездах Дмитрий Валерьевич Бисикало

3

4

5 План лекций 1. Общие сведения о тесных двойных звездах. 2. Наблюдательные проявления элементов структуры течения в тесных двойных звездах. 3. Численное моделирование течений в тесных двойных системах.

6 План лекций 1. Общие сведения о тесных двойных звездах. 2. Наблюдательные проявления элементов структуры течения в тесных двойных звездах. 3. Численное моделирование течений в тесных двойных системах.

7 Численное моделирование течений в тесных двойных системах Андрей Георгиевич Жилкин

8 Часть 1 ЗАЧЕМ? Философское отступление

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

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

11 Вычислительная физика и вычислительная астрофизика Вычислительная физика Вычислительная астрофизика?

12 Вычислительная физика и вычислительная астрофизика Вычислительная физика Вычислительная астрофизика

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

14 Особенности вычислительной астрофизики В астрофизике определенную роль играют все четыре типа взаимодействий: сильное, слабое, электромагнитное и гравитационное. Сильное взаимодействие определяет энергетику звезд за счет ядерных реакций. Слабое взаимодействие является существенным для нейтринной астрофизики. Электромагнитное взаимодействие, в первую очередь, связано с космическими магнитными полями. Однако наиболее важную роль в астрофизике играет гравитационное взаимодействие. Гравитация главным образом определяет эволюцию звездных систем и галактик, черных дыр, сжатие межзвездного газа и образование протозвезд. Велика роль гравитации в процессах аккреции, в том числе и в тесных двойных системах.

15 Блок-схема получения результата ЦельМодельМетод РезультатРасчетКод

16 Часть 2 НА ЧЁМ? Вычислительные ресурсы

17

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

19 Вычислительные кластеры МСЦ РАН МВС Пиковая производительность 10.1 TFLOPS. В июне 2006 года занял 70 место в списке Top500 и стал самым мощным российским суперкомпьютером. В его состав входило 574 двухпроцессорных узла на базе процессоров IBM PowerPC 970FX, работающих на частоте 2.2ГГц и объединённых сетью Myrinet.

20 Вычислительные кластеры МСЦ РАН MBC-6000 Предназначен для решения сложных научно-технических задач. Пиковая производительность составляет 1.64 TFLOPS. Общий объем оперативной памяти решающего поля Гбайт. Решающее поле состоит из 256 процессоров Intel Itanium 2.

21 Вычислительные кластеры МСЦ РАН МВС-100K Пиковая производительность TFLOPS. В настоящий момент является одним из самых мощных суперкомпьютеров СНГ. В состав входят 1275 модуля, оснащенных двумя 4/6-ядерными процессорами Intel Xeon (всего ядра). Для объединения узлов в единое решающее поле используется технология Infiniband.

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

23 Методы распараллеливания

24

25

26 Разбиение расчетной области в случае неравномерной сетки.

27 Методы распараллеливания Разбиение расчетной области в случае использования динамически адаптивной встраиваемой сетки.

28 Эффективность распараллеливания Тестирование производительности параллельного численного кода проводится с помощью измерения времени, необходимого для расчета некоторого количества временных шагов. T(N) – время счета программы на N процессорах. Spd определяет ускорение работы параллельной программы по сравнению с однопроцессорным вариантом. Eff демонстрирует влияние накладных расходов (передача данных по сети) на время счета программы.

29 Эффективность распараллеливания

30

31

32

33 Параллельные программы Система передачи сообщений – Message Passing Interface (MPI). Процессы объединяются в группы в рамках своего коммуникатора. Коммуникаторы могут иметь различные топологии. Взаимодействие между процессами происходит в результате обмена сообщениями. Сообщения могут быть различного типа.

34 Параллельные программы

35 Параллельные вычисления

36

37

38 Пока код считает…

39 Часть 3 ЧТО? Математические модели

40 Типы вычислительных задач С математической точки зрения встречаются такие задачи: 1. Обыкновенные дифференциальные уравнения. 2. Дифференциальные уравнения в частных производных. 3. Интегральные уравнения. 4. Интегро-дифференциальные уравнения. В любом случае различные методы решения состоят в том, что исходная задача сводится к решению некоторой системы линейных алгебраических уравнений. Эта система затем решается на компьютере.

41 Метод расщепления

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

43 Уравнения газовой динамики Уравнение непрерывности. Уравнение движения. Уравнение энергии. Уравнение состояния. Внешние силы.

44 Уравнения магнитной газодинамики Уравнение индукции.

45 Учет других процессов Магнитное поле компактного объекта. В численной модели удобнее вычислять не полное магнитное поле, а только поле b, индуцированное токами в плазме.

46 Учет других процессов Диффузия магнитного поля. Магнитное поле Турбулентный вихрь Токовые слои Два оборота 1. Диссипация токов в турбулентных вихрях аккреционного диска. Магнитное поле Турбулентный вихрь Токовые слои Два оборота Экваториальная плоскость диска К звезде Силовые трубки магнитного поля Условие равновесия трубки 2. Плавучесть силовых трубок тороидального магнитного поля в аккреционном диске.

47 Учет других процессов Процессы нагрева и охлаждения.

48 Часть 4 КАК? Методы расчета

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

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

51 Дискретизация Пример дискретизации одномерной расчетной области.

52 Дискретный анализ Фурье При дискретизации часть информации неминуемо теряется. Возникает некоторая неопределенность в значениях величин. Два примера дискретизации величины u(x). Набор значений в узлах u i дает неполное представление исходной величины u(x). Но чем больше узлов, тем лучше представление.

53 Дискретный анализ Фурье Сеточные функции u i разлагаются в конечные суммы Фурье: Таким образом, при дискретизации теряется вся коротковолновая (большие значения k) информация о решении. Набор значений сеточной функции u i описывает только длинноволновую (длина волны > 2 /N) часть информации о исходной величине u(x). Чем больше N, тем более короткие длины волн можно учесть. В пределе при N, дискретное Фурье-разложение переходит в Фурье- представление для непрерывных функций. Дискретный Фурье-анализ можно использовать для более тонкого анализа разностных схем.

54 Расчетные сетки Все расчетные сетки по типу логических связей между узлами можно разделить на два типа: структурированные и неструктурированные. Структурированные расчетные сетки определяются упорядоченной структурой узлов с явно выраженными сеточными направлениями. В общем случае эти направления можно трактовать как оси некоторой криволинейной системы координат. Ячейки такой сетки являются топологическими прямоугольниками (двумерный случай) или параллелепипедами (трехмерный случай). Узлы такой сетки упорядочены по индексами:

55 Расчетные сетки В криволинейной системе координат (X,Y), определяющейся сеточными направлениями (индексами), структурированная сетка является равномерной по всем координатам.

56 Расчетные сетки Неструктурированные сетки определяются простым набором узлов Логическая связь между узлами сетки задается произвольным образом. Как правило для этого используются треугольные ячейки. Для этого можно использовать алгоритм триангуляции Делоне.

57 Расчетные сетки

58

59 Адаптация сеток Адаптивными называются сетки, каким-либо образом приспосабливающиеся к особенностям течения. Они позволяют без увеличения вычислительных затрат повышать точность расчета. В основе всех методов построения адаптивных сеток положен принцип оптимального распределения узлов сетки в расчетной области. Это распределение узлов должно учитывать взаимное расположение и скорости отдельных подобластей с какими-либо особенностями течения: большие градиенты, сильные разрывы, межфазные границы и т.п. Адаптивные сетки можно разделить на два типа: 1) Геометрически адаптивные 2) Динамически адаптивные

60 Адаптация сеток Динамически адаптивные сетки делятся на а) Адаптивно-подвижные б) Адаптивно-встраиваимые В первом случае количество узлов сетки остается постоянным, но меняется их взаимное расположение. Во втором случае изменяется число узлов сетки.

61 Наша сетка Структура расчетной сетки в экваториальной плоскости (xy) двойной системы.

62 Наша сетка Структура расчетной сетки в вертикальной плоскости (xz) двойной системы.

63 Теорема о сходимости Исходная задачаРазностная задача Согласованность Решение Сходимость Корректность Устойчивость

64 Анализ устойчивости Определение устойчивости формулируется таким образом, чтобы при наличии согласованности получить необходимое и достаточное условие для сходимости. Однако в отличие от согласованности и сходимости, понятие устойчивости не использует информацию об исходной задаче, а является чисто внутренним свойством разностной схемы. Поэтому ей уделяется большое внимание при анализе разностных задач. Основные методы анализа устойчивости схем. 1. Метод дискретных возмущений. 2. Метод дифференциальных приближений. 3. Спектральный метод. 4. Энергетический метод.

65 Диффузия и дисперсия Спектральный метод, основанный на дискретном анализе Фурье, позволяет провести классификацию ошибок разностной схемы. Наиболее важные амплитудные и фазовые ошибки связаны с явлением диффузии и дисперсии на разностной сетке. Диффузия приводит к размытию профилей, а дисперсия – к появлению нефизических осцилляций в численном решении.

66 Классификация схем Классификация схем для системы уравнений газовой динамики определяется способом построения обобщенного решения. 1. Обобщенное решение строится как предел гладких решений некоторой системы уравнений параболического типа. К этому классу можно отнести схемы с искусственной вязкостью и схему Лакса-Вендроффа. 2. Обобщенное решение строится путем сшивки гладких решений на поверхностях разрывов с помощью условий Гюгонио. Наиболее известный такой подход – метод характеристик. 3. Обобщенное решение определяется как решение системы интегральных законов сохранения. К этому классу относятся все годуновские методы.

67 Метод Годунова В 1959 году С.К. Годунов предложил новый метод для численного решения задач одномерной газовой динамики. В этом методе используется процедура решения задачи Римана о распаде произвольного разрыва. Опишем основную идею метода. Рассмотрим базовую схему, полученную с помощью метода конечных объемов: Здесь для удобства ячейки нумеруются целыми индексами, а узлы – полуцелыми. Поскольку мы все равно не знаем распределение функции u(x,t) внутри ячеек, будем считать эту функцию кусочно-постоянной. Тогда ее значение в ячейке совпадает со средним значением u i.

68 Метод Годунова На границах ячеек имеют место произвольные разрывы. Поэтому в следующий момент должен произойти их распад. Произвольные разрывы на границах ячеек.

69 Метод Годунова В результате распада разрыва на границе ячейки функция u(x,t) принимает некоторое значение u * (x i+1/2,t). В силу автомодельности задачи Римана, все бегущие волны являются центрированными. Поэтому величина u * (x i+1/2,t) не зависит от времени.

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

71 Метод Лакса-Фридрихса Запишем соотношения Гюгонио на левой и правой волнах: Мы получили два уравнения для двух неизвестных u * и F *. Решая их, находим:

72 Повышающие поправки Схема называется TVD-схемой (total variation diminishing, схема не увеличивающая полную вариацию решения), если

73 Повышающие поправки Общая идея построение повышающих TVD-поправок заключается в следующем. 1.Выбирается монотонная базовая схема первого порядка аппроксимации. 2.Путем добавки антидиффузионных членов базовая схема преобразуется в схему с повышенным порядком аппроксимации. 3.Строится определенное правило ограничения коэффициентов антидиффузии для обеспечения монотонности. Ограничения потоков производится на основе принципа TVD. Поэтому окончательная схема является TVD-схемой.

74 Повышающие поправки

75

76

77

78

79

80 Визуализация расчетов

81