📄Работа №207261
Тема: Применение усреднённых по Фавру уравнений Навье-Стокса для моделирования перемешивания, обусловленного неустойчивостью Рихтмайера-Мешкова
Тип работы
Дипломные работы, ВКР
Предмет
Школьная математика
Объем:
100 листов
Год:
2020
Просмотров:
41
4240
руб.
🛒 Купить работу
Не подходит эта работа?
Закажите новую по вашим требованиям
Узнать цену на написание
Закажите новую по вашим требованиям
Представленный материал является образцом учебного исследования, примером структуры и содержания учебного исследования по заявленной теме. Размещён исключительно в информационных и ознакомительных целях.
Workspay.ru оказывает информационные услуги по сбору, обработке и структурированию материалов в соответствии с требованиями заказчика.
Размещение материала не означает публикацию произведения впервые и не предполагает передачу исключительных авторских прав третьим лицам.
Материал не предназначен для дословной сдачи в образовательные организации и требует самостоятельной переработки с соблюдением законодательства Российской Федерации об авторском праве и принципов академической добросовестности.
Авторские права на исходные материалы принадлежат их законным правообладателям. В случае возникновения вопросов, связанных с размещённым материалом, просим направить обращение через форму обратной связи.
Настоящий учебно-методический информационный материал размещён в ознакомительных и исследовательских целях и представляет собой пример учебного исследования. Не является готовым научным трудом и требует самостоятельной переработки.
📋 Содержание
ВВЕДЕНИЕ 6
1 Обзор литературы 8
1.1 Гидродинамические неустойчивости и турбулентность в физике
высоких плотностей энергии 8
1.1.1 Астрофизические исследования 8
1.1.2 Перемешивание 10
1.2 Неустойчивость Рихтмайера-Мешкова и перемешивание 10
1.3 Обзор экспериментальных и модельных исследований
неустойчивости Рихтмайера-Мешкова 12
1.3.1 Эксперименты 13
1.3.2 Численное моделирование 16
2 к-sмодель турбулентности и численные методы 20
2.1 Введение 20
2.2 Уравнения k-sмодели турбулентности 21
2.2.1 Усреднение по Рейнольдсу и Фавру 21
2.2.2 Усреднение уравнений Навье-Стокса 25
2.2.3 Замыкание усреднённой системы уравнений 28
2.3 Модификация Авраменко 35
2.4 Применение численных методов 36
2.4.1. Газодинамический этап 37
2.4.2 Диффузионно-турбулентный этап 50
2.5 Начальные условия для неустойчивости Рихтмайера-Мешкова 60
3 Численное моделирование 62
3.1 Введение 62
3.2 Тестирование газодинамических методов 63
3.3 Моделирование эспериментов Леинова и др 74
3.3.1 Моделирование эксперимента №1 75
3.3.2 Моделирование эксперимента №2 76
3.3.3 Моделирование эксперимента №3 77
3.3.4 Моделирование эксперимента №4 78
3.3.5 Моделирование эксперимента №5 79
3.3.6 Моделирование эксперимента №6 80
3.4 Моделирование экспериментов Веттера и Стуртеванта 81
3.4.1 Моделирование эксперимента №1 81
3.4.2 Моделирование эксперимента №2 83
3.4.3 Моделирование эксперимента №3 84
3.5 Моделирование эксперимента Погги и др 86
4 Исследования на сходимость 88
4.1 Введение 88
4.2 Классическая k-sмодель 88
4.2.1 Сходимость в эксперименте Леинова и др. №1 88
4.2.2 Сходимость в эксперименте Веттера и Стуртеванта №1 89
4.2.3 Сходимость в эксперименте Погги и др 89
4.2.4 Выводы 90
4.3 Модификация Авраменко 91
4.3.1 Сходимость в эксперименте Леинова и др. №1 91
4.3.2 Сходимость в эксперименте Веттера и Стуртеванта №1 91
4.3.3 Сходимость в эксперименте Погги и др 92
4.3.4 Выводы 93
ЗАКЛЮЧЕНИЕ 94
БИБЛИОГРАФИЧЕСКИЙ СПИСОК 96
Приложения должны быть в работе, но в данный момент отсутствуют
1 Обзор литературы 8
1.1 Гидродинамические неустойчивости и турбулентность в физике
высоких плотностей энергии 8
1.1.1 Астрофизические исследования 8
1.1.2 Перемешивание 10
1.2 Неустойчивость Рихтмайера-Мешкова и перемешивание 10
1.3 Обзор экспериментальных и модельных исследований
неустойчивости Рихтмайера-Мешкова 12
1.3.1 Эксперименты 13
1.3.2 Численное моделирование 16
2 к-sмодель турбулентности и численные методы 20
2.1 Введение 20
2.2 Уравнения k-sмодели турбулентности 21
2.2.1 Усреднение по Рейнольдсу и Фавру 21
2.2.2 Усреднение уравнений Навье-Стокса 25
2.2.3 Замыкание усреднённой системы уравнений 28
2.3 Модификация Авраменко 35
2.4 Применение численных методов 36
2.4.1. Газодинамический этап 37
2.4.2 Диффузионно-турбулентный этап 50
2.5 Начальные условия для неустойчивости Рихтмайера-Мешкова 60
3 Численное моделирование 62
3.1 Введение 62
3.2 Тестирование газодинамических методов 63
3.3 Моделирование эспериментов Леинова и др 74
3.3.1 Моделирование эксперимента №1 75
3.3.2 Моделирование эксперимента №2 76
3.3.3 Моделирование эксперимента №3 77
3.3.4 Моделирование эксперимента №4 78
3.3.5 Моделирование эксперимента №5 79
3.3.6 Моделирование эксперимента №6 80
3.4 Моделирование экспериментов Веттера и Стуртеванта 81
3.4.1 Моделирование эксперимента №1 81
3.4.2 Моделирование эксперимента №2 83
3.4.3 Моделирование эксперимента №3 84
3.5 Моделирование эксперимента Погги и др 86
4 Исследования на сходимость 88
4.1 Введение 88
4.2 Классическая k-sмодель 88
4.2.1 Сходимость в эксперименте Леинова и др. №1 88
4.2.2 Сходимость в эксперименте Веттера и Стуртеванта №1 89
4.2.3 Сходимость в эксперименте Погги и др 89
4.2.4 Выводы 90
4.3 Модификация Авраменко 91
4.3.1 Сходимость в эксперименте Леинова и др. №1 91
4.3.2 Сходимость в эксперименте Веттера и Стуртеванта №1 91
4.3.3 Сходимость в эксперименте Погги и др 92
4.3.4 Выводы 93
ЗАКЛЮЧЕНИЕ 94
БИБЛИОГРАФИЧЕСКИЙ СПИСОК 96
Приложения должны быть в работе, но в данный момент отсутствуют
📖 Введение
Перемешивание, вызванное гидродинамической неустойчивостью, привлекает большое внимание исследователей во многих областях науки и техники. Интерес к неустойчивости Рихтмайера-Мешкова важен для улучшения современных методов прогнозирования турбулентности в сложных быстропротекающих процессах физики высоких плотностей энергии. В данной работе проводится исследование ударно-турбулентного перемешивания двух веществ, вызванное неустойчивостью Рихтмайера- Мешкова с ударной волной (падающей изначально и отражённой), с использованием многокомпонентной системы уравнений Навье-Стокса, усредненной по Фавру, включающей молекулярный перенос в смеси, по k- sмодели турбулентности. Уравнения модели решаются при помощи конечно-объёмных методов.
Отличием от исследований [1-5] в первую очередь является использование для расчётов схемы первого порядка, в то время как ранее проводились расчёты по схемам третьего порядка и выше. Схемы высших порядков являются более сложными в плане реализации, а также, в силу их сложности, расчёт по ним занимает больше времени, однако, стоит отметить более высокую точность таких схем. Тем не менее, поскольку расчётов на схемах первого порядка не проводилось, то отсутствуют и исследования, оправдывающие использование более сложных схем для решения системы уравнений k-sмодели турбулентности. В данной работе хоть и не будет проведено сравнительного анализа решения на схемах разного порядка, но будет проведено сравнение результатов расчёта с экспериментальными данными.
Из всего вышесказанного можно отметить главную цель настоящей работы, а именно создание программы моделирования перемешивания, с помощью предложенной k-sмодели с использованием численной схемы первого порядка точности. В связи с этим были поставлены следующие задачи:
• Реализовать численное решение уравнений k-sмодели турбулентности
• Провести моделирование перемешивания, обусловленного неустойчивостью Рихтмайера-Мешкова
• Сравнить полученные результаты с экспериментальными данными
• Оценить работу предложенной k-sмодели с использованием численной схемы первого порядка точности для моделирования перемешивания, обусловленного неустойчивостью Рихтмайера- Мешкова
Отличием от исследований [1-5] в первую очередь является использование для расчётов схемы первого порядка, в то время как ранее проводились расчёты по схемам третьего порядка и выше. Схемы высших порядков являются более сложными в плане реализации, а также, в силу их сложности, расчёт по ним занимает больше времени, однако, стоит отметить более высокую точность таких схем. Тем не менее, поскольку расчётов на схемах первого порядка не проводилось, то отсутствуют и исследования, оправдывающие использование более сложных схем для решения системы уравнений k-sмодели турбулентности. В данной работе хоть и не будет проведено сравнительного анализа решения на схемах разного порядка, но будет проведено сравнение результатов расчёта с экспериментальными данными.
Из всего вышесказанного можно отметить главную цель настоящей работы, а именно создание программы моделирования перемешивания, с помощью предложенной k-sмодели с использованием численной схемы первого порядка точности. В связи с этим были поставлены следующие задачи:
• Реализовать численное решение уравнений k-sмодели турбулентности
• Провести моделирование перемешивания, обусловленного неустойчивостью Рихтмайера-Мешкова
• Сравнить полученные результаты с экспериментальными данными
• Оценить работу предложенной k-sмодели с использованием численной схемы первого порядка точности для моделирования перемешивания, обусловленного неустойчивостью Рихтмайера- Мешкова
✅ Заключение
В рамках данной работы была успешно создана программа моделирования перемешивания с помощью предложенной k-sмодели с использованием численной схемы первого порядка точности. Реализовано численное решение уравнений k-sмодели турбулентности, проведено моделирование перемешивания, обусловленного неустойчивостью Рихтмайера-Мешкова, полученные результаты сравнивались с экспериментальными данными. Произведена оценка работы предложенной k-sмодели с использованием численной схемы первого порядка точности для моделирования перемешивания, обусловленного неустойчивостью Рихтмайера-Мешкова.
Проводимое в данной работе моделирование перемешивания, вызванного неустойчивостью Рихтмайера-Мешкова при помощи усреднённых по Фавру уравнений Навье-Стокса показало, что численные результаты согласуются с экспериментальными данными.
Для решения усреднённых уравнений было принято решение разделить процесс вычисления на два этапа. На каждом временном слое сначала решались уравнения обычной газовой динамики с помощью схемы Хэнкока с реконструктором MUSCL в примитивных переменных и решателем HLLC, а после уравнения диффузионных для диффузионных потоков и источников турбулентности при помощи неявной схемы первого порядка. Несмотря на второй порядок точности для первого этапа, вследствие как раз таки разделения на два этапа и использования на втором схемы первого порядка, результирующий порядок схемы не может быть выше первого. Однако, не смотря на это, были получены весьма неплохие результаты для всех экспериментов, кроме эксперимента Веттера и Стуртеванта (Vetter and Sturtevant) с числом Маха падающей ударной волны 1,24.
Описанные выше результаты были получены благодаря подбору параметров модели, причём расчётная сетка тоже являлась таким параметром, поскольку при использовании предложенной k-sмодели наблюдалось отсутствие численной сходимости при повторном прохождении ударной волны через зону турбулентного перемешивания, что объясняется неправильной работой данной модели при наличии ударных волн. Данная проблема была решена при помощи модификации Авраменко, однако, это привело к необходимости повторного подбора параметров модели, который уже не проводился в рамках данной работы.
Проводимое в данной работе моделирование перемешивания, вызванного неустойчивостью Рихтмайера-Мешкова при помощи усреднённых по Фавру уравнений Навье-Стокса показало, что численные результаты согласуются с экспериментальными данными.
Для решения усреднённых уравнений было принято решение разделить процесс вычисления на два этапа. На каждом временном слое сначала решались уравнения обычной газовой динамики с помощью схемы Хэнкока с реконструктором MUSCL в примитивных переменных и решателем HLLC, а после уравнения диффузионных для диффузионных потоков и источников турбулентности при помощи неявной схемы первого порядка. Несмотря на второй порядок точности для первого этапа, вследствие как раз таки разделения на два этапа и использования на втором схемы первого порядка, результирующий порядок схемы не может быть выше первого. Однако, не смотря на это, были получены весьма неплохие результаты для всех экспериментов, кроме эксперимента Веттера и Стуртеванта (Vetter and Sturtevant) с числом Маха падающей ударной волны 1,24.
Описанные выше результаты были получены благодаря подбору параметров модели, причём расчётная сетка тоже являлась таким параметром, поскольку при использовании предложенной k-sмодели наблюдалось отсутствие численной сходимости при повторном прохождении ударной волны через зону турбулентного перемешивания, что объясняется неправильной работой данной модели при наличии ударных волн. Данная проблема была решена при помощи модификации Авраменко, однако, это привело к необходимости повторного подбора параметров модели, который уже не проводился в рамках данной работы.
ИЛИ
Поиск аналога
📕 Список литературы
🛒 Оформить заказ
⚡ Работу высылаем в течении 5 минут после оплаты.