авторефераты диссертаций БЕСПЛАТНАЯ РОССИЙСКАЯ БИБЛИОТЕКА - WWW.DISLIB.RU

АВТОРЕФЕРАТЫ, ДИССЕРТАЦИИ, МОНОГРАФИИ, НАУЧНЫЕ СТАТЬИ, КНИГИ

 
<< ГЛАВНАЯ
АГРОИНЖЕНЕРИЯ
АСТРОНОМИЯ
БЕЗОПАСНОСТЬ
БИОЛОГИЯ
ЗЕМЛЯ
ИНФОРМАТИКА
ИСКУССТВОВЕДЕНИЕ
ИСТОРИЯ
КУЛЬТУРОЛОГИЯ
МАШИНОСТРОЕНИЕ
МЕДИЦИНА
МЕТАЛЛУРГИЯ
МЕХАНИКА
ПЕДАГОГИКА
ПОЛИТИКА
ПРИБОРОСТРОЕНИЕ
ПРОДОВОЛЬСТВИЕ
ПСИХОЛОГИЯ
РАДИОТЕХНИКА
СЕЛЬСКОЕ ХОЗЯЙСТВО
СОЦИОЛОГИЯ
СТРОИТЕЛЬСТВО
ТЕХНИЧЕСКИЕ НАУКИ
ТРАНСПОРТ
ФАРМАЦЕВТИКА
ФИЗИКА
ФИЗИОЛОГИЯ
ФИЛОЛОГИЯ
ФИЛОСОФИЯ
ХИМИЯ
ЭКОНОМИКА
ЭЛЕКТРОТЕХНИКА
ЭНЕРГЕТИКА
ЮРИСПРУДЕНЦИЯ
ЯЗЫКОЗНАНИЕ
РАЗНОЕ
КОНТАКТЫ

Pages:   || 2 | 3 | 4 |

Численное моделирование сжимаемой турбулентности в проблеме образования и эволюции звёзд

-- [ Страница 1 ] --

Институт прикладной математики им.М.В.Келдыша РАН

На правах рукописи

Устюгов Сергей Дмитриевич

Численное моделирование сжимаемой турбулентности в проблеме образования и эволюции звёзд

05.13.18 – Математическое моделирование, численные методы и комплексы программ

А в т о р е ф е р а т

Диссертации на соискание ученой степени доктора физико-математических наук

Москва – 2012

Работа выполнена в Институте прикладной математики им.М.В.Келдыша РАН

Официальные оппоненты: член-корр. РАН,

Петров И.Б.

доктор физико-математических наук,

Жуков В.Т.

доктор физико-математических наук:

Попов В.Ю.

Ведущая организация: Институт астрономии РАН

Защита состоится «___»____________2012 г. в «_____» час. на заседании диссертационного совета Д002.024.02 в Институте прикладной математики им.М.В.Келдыша РАН по адресу: 125047, Москва, Миусская пл., д.4.

С диссертацией можно ознакомиться в библиотеке Института прикладной математики им.М.В.Келдыша РАН

Автореферат разослан «___»____________2012 г.

Ученый секретарь

диссертационного совета Д002.024.02

доктор физико-математических наук Змитренко Н.В.

Общая характеристика работы

Актуальность темы

В настоящее время считается, что сжимаемая турбулентность является необходимым фактором регулирующим темп образования звёзд в межзвёздных молекулярных облаках. Сверхзвуковая турбулентность создает сильные флуктуации плотности в молекулярном газе, которые сжимаются под действием самогравитации. Этот процесс ведёт, в конечном счете, к образованию звёзд. Турбулентность играет двоякую роль в процессе звездообразования. С одной стороны на больших масштабах турбулентность дает динамическую поддержку молекулярному облаку, с другой стороны она облегчает локальное сжатие вещества. Таким образом, рождение звёзд определяется динамическим поведением вещества внутри родительского облака, которое определяет время и место рождения протозвёзд, процесс их сжатия и роста в результате аккреции вещества из падающей внешней оболочки. Турбулентность в межзвёздной среде является сверхзвуковой, сжимаемой, замагниченной и поэтому отличается от несжимаемой турбулентности, проявляющейся на Земле, например, в атмосфере и океанах. В случае несжимаемой турбулентности диссипация энергии происходит на динамическом масштабе l меньше масштаба, на которых действует вязкость. Каскад энергии идет от крупных масштабов, через взаимодействие рожденных вихрей, к все более и более мелким масштабам, до тех пор, пока не достигает масштаба диссипации. В устойчивом состоянии равное количество энергии переносится от одного масштаба к другому в инерциальном диапазоне. При этом дифференциальный спектр кинетической энергии имеет вид . Турбулентное течение в межзвездной среде отличается от этой картины в следующем. Во-первых, турбулентность является сильно сжимаемой с числом Маха в диапазоне от единицы в теплой диффузной среде вплоть до значения 50 в холодных и плотных молекулярных облаках. Во-вторых, уравнение состояния газа очень мягкое, благодаря действию радиационного охлаждения, с индексом политропы в выражении для давления лежащим в интервале в зависимости от давления и температуры [1]. В-третьих, турбулентность направляется из множества источников, действующих на различных масштабах. При сверхзвуковой турбулентности ударные волны дают дополнительный источник диссипации кинетической энергии. Ударные волны могут переносить энергию между широко разделенными масштабами, удаляя тем самым локальную природу турбулентного каскада энергии характерного для несжимаемой турбулентности. Результирующий спектр кинетической энергии в этом случае приобретает вид . Численные моделирования, также как и наблюдения межзвездной среды показывают, что наклон спектра лежит в диапазоне между Колмогоровским значением и значением (Бюргерс турбулентность) [2].



Физическая природа турбулентности в межзвездной среде еще ясно не определена. В некоторых моделях предполагалось, что источником турбулентности на малых масштабах могут быть звездные ветры и истечения [3], процессы, связанные с взрывом сверхновых [4] или расширением HII областей [5]. Другой привлекательной идеей является возникновение турбулентных течений в местах столкновения крупномасштабных сходящихся течений [6]. Численные моделирования столкновения струй теплого нейтрального газа показали, что этот процесс является достаточным, чтобы поддержать высокую степень турбулентности при образовании облаков [7]. В общем, картина звездообразования при наличии сжимаемой турбулентности выглядит следующим образом. На больших масштабах в молекулярных облаках турбулентность является сильно сверхзвуковой. Это означает возможность образования больших контрастов плотности, которые реально наблюдаются в молекулярных облаках, имеющих среднюю плотность около 100 частиц на один кубический сантиметр и содержащих ядра с плотностью и более. В областях образования скоплений звезд некоторые флуктуации плотности могут превышать критическую массу для начала гравитационного коллапса. Действие турбулентности приводит вначале к образованию кластера таких флуктуаций плотности и затем под влиянием гравитации из некоторого их подмножества образуются звёзды [4]. В этом случае формируется более широкий спектр масс звезд, чем в случае образования одиночных звезд при отсутствии турбулентности. В случае ядер малой массы турбулентность является дозвуковой, и поле скорости становится более когерентным [8]. Такие структуры не подвергаются фрагментации и производят одиночные звезды или двойные системы.

Очень трудно моделировать в лаборатории физические условия, характерные для существования сверхзвуковой турбулентности в молекулярных облаках, и также ввиду недостатка необходимой информация из наблюдений. Таким образом, численное моделирование в настоящее время представляет собой единственный способ исследовать статистику сверхзвуковой турбулентности, которая является ключевой информацией любой успешной статистической теории звездообразования [9]. Долгое время, природа сильно сжимаемой и замагниченной межзвездной турбулентности оставалась слабо разработанной. Численное моделирование дает возможность оценить перенос энергии между масштабами и найти ключевые пространственные корреляции между основными характеристиками турбулентных течений. Это стимулирует развитие аккуратных численных методов удобных для моделирования сжимаемой турбулентности в молекулярных облаках, обладающих низкой численной диссипацией и широким инерционным интервалом в спектре мощности основных физических характеристик. Многие численные методы, которые успешно проходят стандартные МГД тесты, становятся неустойчивыми при решении задачи развития сверхзвуковой МГД турбулентности [10]. Добавление небольшой диссипации в численную схему для обеспечения ее устойчивости приводит, в конечном счете, к деградации турбулентных статистик [11]. Поэтому проблема численной устойчивости становится решающей для моделирования, как сверхзвуковой турбулентности, так и самого процесса звездообразования.

Из наблюдений найдено, что поле скоростей в солнечной фотосфере образовано турбулентными конвективными движениями. Эти движения в свою очередь указывают на существование двух основных масштабов развития конвекции – грануляции и супергрануляции. Наличие промежуточного масштаба мезогрануляции до сих пор вызывает многочисленные споры [12]. Ячейки грануляции размером 1.5 Мм отчетливо наблюдаются в виде яркого клеточного узора, который является отображением вертикального движения газа, определенного из сдвига спектральных линий. Структуры супергрануляции покрывают поверхность покоящегося Солнца с типичным горизонтальным масштабом 30 Мм и средним временем жизни около 1.8 дней. Ячейки супергрануляции легко детектируемы по картам допплеровских скоростей в виде флуктуаций горизонтальной компоненты поля скоростей, среднеквадратичное значение которых равно 360 м/сек.

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

В ходе исследования супергрануляции было предложено несколько теорий ее происхождения. Вначале предполагалось [13], что локальный минимум вертикального адиабатического градиента температуры в глубоких слоях Солнца, вызванный ионизацией He I и He II должен способствовать образованию конвективных ячеек с горизонтальным размером сравнимым с глубиной зон ионизации. С другой точки зрения супергрануляция может представлять собой пространственное упорядочивание большого числа маленьких конвективных ячеек [14]. Применение различных методов локальной гелиосейсмологии и долгие серии наблюдений привели к открытию волновых свойств супергрануляции [15], которое означает, что супергрануляция есть пример волновой конвекции.

Понимание процесса развития турбулентной конвекции в солнечной фотосфере сильно возросло с развитием вычислительных ресурсов. Моделирование солнечной конвекции проводилось в рамках двух подходов. В первом случае использовались идеализированные модели, в которых брались простые политропные модели стратифицированной атмосферы и решались стандартные уравнения несжимаемой и сжимаемой гидродинамики с учетом вязкости, теплопроводности, магнитной диффузии в ограниченной области с искусственными граничными условиями [16]. Во втором случае моделирование проводилось с использованием реалистичных физических условий и с включением важных физических факторов, таких как, например, перенос излучения и реалистичное уравнение состояния [17]. В этом случае результаты численного моделирования могут быть сравнены прямо с наблюдениями. В отличие от идеализированных моделей вычисления в большинстве случаев игнорируют физическую вязкость и основываются только на численной вязкости, применяемой для обеспечения устойчивости вычислений. Поэтому такое моделирование не принимает в расчет все многообразие временных и пространственных масштабов диссипативных процессов. Однако недавние прямые моделирования показали, что зависимость статистических свойств различных турбулентных МГД течений от диссипативных процессов имеет большое значение, например, в проблеме флуктуаций в динамо-механизме генерации магнитных полей [18].

При численном моделировании солнечной супергрануляции применяются два различных геометрических приближения: локальное моделирование в небольшом объеме солнечной поверхности и глобальное моделирование в сферической оболочке. В локальном приближении размер вычислительной области обычно равен 1-3 масштабам супергрануляции. Шаг численной сетки в этом случае не должен быть меньше величины, необходимой для расчета конвекции на масштабе грануляции. Отметим также, что динамика на масштабе супергрануляции сильно ограничена периодическими граничными условиями в горизонтальном направлении. В глобальном сферическом приближении наименьший масштаб сравним с размером супергрануляции, подразумевая, что эволюция масштабов супергрануляции сильно диссипативна. Динамика конвекции на масштабе грануляции в этом приближении полностью исключена из рассмотрения. В целом определенные ограничения каждого типа моделирования еще не позволяют исследовать нелинейную динамику и перенос энергии на масштабе супергрануляции полностью самосогласованно. Тем не менее, численное моделирование даёт полезную информацию о крупномасштабной динамике конвекции и магнитных полях в покоящемся Солнце.





Ранние моделирования в локальном приближении [16,19], не показали наличия масштабов супергрануляции, но было обнаружено существование медленно эволюционирующего поля скоростей на масштабе в пять раз превосходящих размер грануляции. Типичное время корреляции и перенос энергии на этом масштабе были намного больше, чем подобные величины для турбулентной конвекции на масштабе грануляции. В итоге исследований был сделан вывод, что развитие конвекции на масштабе мезогрануляции есть результат коллективного взаимодействия маленьких ячеек грануляции. Реалистичное трехмерное моделирование конвекции [20,21] обнаружило монотонное гладкое увеличение размера конвективных ячеек с глубиной и очень малое изменение спектра мощности скорости соответствующее шкале супергрануляции.

Несколько локальных моделирований развития конвекции проведено с учетом магнитного поля. Их можно разделить на два класса, по способу введения магнитного поля в начальные условия моделирования. В первом случае исследователи вводили среднее магнитное поле заданной величины и конфигурации в начальные условия, во втором магнитное поле генерировалось в ходе моделирования, стартуя с небольшого начального возмущения магнитного поля. Моделирование со слабым начальным магнитным полем [22] показало, что магнитное поле стремилось организоваться в сеть на масштабах мезогрануляции, также как в спокойном Солнце. В отсутствии начального среднего магнитного поля [23], результаты показали, что мелкомасштабное хаотическое магнитное поле, созданное динамо механизмом, организуется в крупномасштабные магнитные структуры на масштабе супергрануляции. В противоположном пределе сильного среднего магнитного поля, было обнаружено, что основной масштаб турбулентной конвекции сильно уменьшен магнитным натяжением.

Трехмерное моделирование локальной МГД конвекции с реалистичной физикой на масштабе супергрануляции было впервые проведено в работах [24]. Вычисления проводились в расчетной области шириной 60~Мм и глубиной 20~Мм. В начальный момент времени в области задавалось однородное вертикальное магнитное поле величиной 50 Гс. Результаты моделирования показали образование крупномасштабной структуры супергрануляции, концентрацию магнитного поля в расширенных областях размером около 5 Мм и в слоях расположенных по краям ячеек супергрануляции. После этого в работе [25] было проведено МГД моделирование конвекции в области 48 Мм шириной и 20 Мм глубиной с магнитным полем, введенным начально на нижней границе расчетной области, в форме однородного горизонтального потока. Хотя время моделирования только соответствовало времени жизни супергрануляции, результаты работы показали, что выметание магнитных элементов к границам структур супергрануляции приводит к образованию магнитной сети, во многом схожей с той, что была получена ранее в моделировании [24].

Таким образом, численное моделирование супергрануляции становится очень активной областью исследований. Большинство недавних, крупномасштабных моделей успешно воспроизводят наблюдения фотосферной, солнечной магнитоконвекции на масштабе от грануляции до супергрануляции. Одним из главных достижений является очень хорошее совпадение результатов моделирования с данными, полученными при обработке средствами гелиосейсмологии реальных наблюдений [26,27].

Цели диссертации:

1. Разработка численного метода высокого порядка с низкой численной диссипацией с хорошими спектральными характеристиками для моделирования сжимаемой МГД турбулентности для задачи звездообразования.

2. Исследование сверхзвуковой МГД турбулентности в задаче образования звёзд. Получение универсальных спектральных характеристик форсированной турбулентности для изотермического и мультифазного случая.

3. Исследование локальной солнечной супергрануляции с реалистичной физикой с учетом магнитного поля. Определение характерных масштабов развития солнечной конвекции.

Научная новизна:

1. Для гиперболических законов сохранения предложена разностная схема повышенного порядка точности с низкой численной диссипацией. Схема позволяет с высоким разрешением моделировать турбулентные сверхзвуковые МГД течения без возникновения нефизических осцилляций на основных МГД разрывах и искажения спектральных характеристик турбулентности.

2. Построена новая физико-математическая модель сжимаемой МГД турбулентности действующей при образовании иерархических крупномасштабных структур в диффузной межзвездной среде и в областях звездообразования в плотных молекулярных облаках. В основе модели лежит представление о турбулентности как механизме преобразования кинетической энергии поступающей от внешних источников на больших масштабах при одновременном действии случайной вынуждающей силы. В рамках этой модели впервые в мировой практике показано, что закон 4/3 для несжимаемой МГД турбулентности может быть расширен для сверхзвуковых турбулентных течений.

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

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

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

Научная и практическая ценность:



Pages:   || 2 | 3 | 4 |
 

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







 
© 2013 www.dislib.ru - «Авторефераты диссертаций - бесплатно»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.