Помощь в написании студенческих работ
Антистрессовый сервис

Математическое моделирование течений жидкости и газа на основе разрывного метода Галеркина

Диссертация: Математическое моделирование течений жидкости и газа на основе разрывного метода Галеркина
ДиссертацияПомощь в написанииУзнать стоимостьмоей работы

В последнее время наблюдается возрастающий интерес к RKDG-методу в связи с необходимостью численного решения задач динамики многофазных сред, описываемых неконсервативными гиперболическими системами уравнений. В качестве математической модели двухфазных сред наиболее часто используются уравнения Баера-Нунциато, которые в трехмерном случае представляют собой систему одиннадцати уравнений в частных… Читать ещё >

Математическое моделирование течений жидкости и газа на основе разрывного метода Галеркина (реферат, курсовая, диплом, контрольная)

Содержание

  • Глава 1. Применение разрывного метода Галеркина для численного моделирования течений жидкости и газа
    • 1. 1. Численное решение одномерных уравнений с использованием ШШС-метода
      • 1. 1. 1. Схема ШШС-метода для одномерного гиперболического уравнения
      • 1. 1. 2. Решение одномерного квазилинейного уравнения переноса
    • 1. 2. Численное моделирование двумерных течений идеального газа на основе ШШС-метода
      • 1. 2. 1. Схема ШШС-метода для гиперболических уравнений
      • 1. 2. 2. Численное решение уравнений Эйлера
  • Глава 2. Программная реализация и оптимизация алгоритма разрывного метода Галеркина
    • 2. 1. Численное моделирование двумерных течений вязкого теплопроводного газа на основе ШШС-метода
      • 2. 1. 1. Схема ШШС-метода для уравнений конвекции-диффузии
      • 2. 1. 2. Численное решение уравнений Навье-Стокса
    • 2. 2. Оптимизация параметров лимитера в алгоритме ШШС-метода
      • 2. 2. 1. Классический лимитер
      • 2. 2. 2. Алгоритм автоматического выбора параметров лимитера
      • 2. 2. 3. Результаты расчетов
    • 2. 3. Сравнение вычислительных затрат метода конечных объемов и ШШС-метода
    • 2. 4. Параллельный алгоритм ШШС-метода
    • 2. 5. Организация программного комплекса
    • 2. 6. Применение ШШС-метода к расчету аэродинамических характеристик крылового профиля КАСА23 012 с закрылком
      • 2. 6. 1. Постановка задачи
      • 2. 6. 2. Результаты расчетов
  • Глава 3. Численное моделирование двухфазных потоков с использованием разрывного метода Галеркина
    • 3. 1. Математическая модель и задача Римана
      • 3. 1. 1. Одномерные уравнения Баера-Нунциато
      • 3. 1. 2. Точное решение задачи Римана
    • 3. 2. НЬЬС-метод решения задачи Римана .97 ¦
    • 3. 3. Тестирование НЬЬС-метода на локальных задачах Римана
    • 3. 4. НЬЬС-метод решения задачи Римана в методе конечных объемов, ИКБС-методе и РС-методе
      • 3. 4. 1. Численный поток типа НЬЬС
      • 3. 4. 2. Метод конечных объемов
      • 3. 4. 3. ШФС-метод
      • 3. 4. 4. РС-метод
    • 3. 5. Анализ эффективности ШШС-метода с потоком типа НЬЬС по сравнению с другими методами
      • 3. 5. 1. Задачи Римана
      • 3. 5. 2. Задача о взаимодействии двух материалов
  • Выводы

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

Математическое моделирование течений жидкости и газа является важнейшим элементом решения сложных инженерных задач [1−13]. Основным инструментом моделирования таких течений является вычислительный эксперимент [14]. При рассмотрении реальных технических устройств, как правило, приходится работать с очень сложной геометрией, что приводит к необходимости использования в расчетах неструктурированных сеток.

Одним из главных требований, предъявляемых к методам решения задач газовой динамики, является правильность воспроизведения решения в областях, где оно претерпевает сильные изменения во времени и пространстве, в частности, на ударных волнах, волнах разрежения и контактных разрывах [15−17]. Для решения задач газовой динамики часто применяются конечно-объемные схемы типа Годунова [15,18]. Для получения конечно-объемной схемы исходные уравнения интегрируются по элементу (ячейке) сетки. Эта процедура соответствует методу Галеркина-Петрова, в котором в качестве пробных функций используются кусочно-постоянные характеристические функции ячеек сетки. В этом случае определяемые в ходе решения величины являются средними значениями искомого решения на ячейке сетки. Таким образом, решение на каждой ячейке сетки приближается постоянной функцией. При построении простейшего метода конечных объемов межэлементные граничные потоки вычисляются с использованием свойств приближенного решения задачи Римана о распаде разрыва [15], в которой исходными данными служат средние значения численного решения в смежных ячейках. Такая схема обеспечивает только первый порядок аппроксимации по пространству, поэтому для качественной передачи особенностей решения приходится очень сильно измельчать сетку. Этого можно избежать в методах конечных объемов повышенного порядка аппроксимации, в которых межэлементные потоки вычисляются на основе величин, полученных полиномиальной интерполяцией средних значений решения. Однако, как правило, повышение порядка метода связано с расширением шаблона аппроксимации, что может отрицательно сказаться на качестве решения на неструктурированных сетках [18], а также вычислительных затратах и эффективности распараллеливания метода. Поэтому предпочтительным является использование в расчетах численных методов, обладающих высоким порядком аппроксимации и сохраняющих при этом компактность шаблона аппроксимации. Одним из таких методов является метод RKDG (Runge-Kutta discontinuous Galerkin) [19−24], развитию и применению которого и посвящена настоящая работа.

Идеи разрывного метода Галеркина, или DG-метода, впервые предложены в работе [25] (Reed, Hill) для численного решения линейного уравнения переноса нейтронов вида где, а б R, а, а — некоторый постоянный вектор. Согласно предложенной в [25] схеме, решение данного уравнения вычислялось последовательно для каждого элемента сетки при условии, что ячейки сетки были расположены вдоль соответствующего характеристического направления.

В дальнейшем Lesaint и Raviart [26] провели подробный анализ DG-метода применительно к задачам данного класса и показали, что порядок сходимости метода при использовании в качестве базисных функций полиномов степени не выше к равен (Ах)к для сеток общего вида и (Ах)к+1 для прямоугольных сеток. Позднее в работе [27] (Johnson, Pitkaranta) получен порядок сходимости (Aa:)fc+½ для произвольных сеток, а в работе [28] (Peterson) доказано, что эта оценка порядка является оптимальной. В работе [29] (Richter) получена оценка порядка сходимости (Ах)к+1 для некоторых неортогональных двумерных сеток.

Положительные результаты применения DG-метода к решению линейных гиперболических задач мотивировали дальнейшие исследования возможности использования метода для решения нелинейных гиперболических уравнений и систем вида с необходимыми начальными и граничными условиями, где и —вектор неизвестных. Пристутствующие в системе нелинейные потоки ^(и) приводят к невозможности рассмотрения ее как совокупности уравнений переноса и поэлементного вычисления решения. Применение первоначальной схемы ОС-метода, предложенной в [25], вело к необходимости решения нелинейной сиои -1- div (an) = /, г=1 стемы уравнений относительно всех неизвестных функций, что делало данную схему крайне неэффективной для сложных гиперболических систем.

Для преодоления указанной трудности СЬауегЛ и Эа1гапо [30] впервые построили явную схему БС-метода для одномерного гиперболического уравнения в дивергентной форме. Для этого использовалась БС-дискретизация исходных уравнений по пространству кусочно-линейными функциями, а получаемая в результате система обыкновенных дифференциальных уравнений относительно коэфициетов аппроксимирующей кусочно-линейной функции решалась явным методом Эйлера. Однако анализ устойчивости построенной схемы показал, что она являлась безусловно неустойчивой при условии постоянства отношения шагов по времени и пространству АЬ/Ах) схема устойчива только если отношение АЬ/ Ах имеет порядок л/Ах, что является слишком жестким ограничением для гиперболических задач.

СЬауе^ и СоскЬигп в работе [31] предложили метод улучшения устойчивости схемы с помощью так называемых ограничителей [32,33]. Таким образом им удалось получить схему, которая относилась к классу ТУБМ (т.е. к классу схем с невозрастающей полной вариацией средних значений) и классу ТУВ (т.е. к классу схем с ограниченной полной вариацией) при фиксированном числе Куранта Ссрь = которое выбиралось из условия Ссрь ^ §• Таким образом гарантировалась сходимость решения к точному. Основным недостатком предложенной схемы являлось то, что она имела только первый порядок аппроксимации по времени при втором порядке аппроксимации по пространству. Ограничитель при этом должен обеспечивать высокий порядок аппроксимации в областях, где решение является гладкой функцией и в то же время подавлять нефизичные осцилляции на разрывах решения.

Этот недостаток исправлен в работе СоскЬигп и ЭИи [34], в которой впервые представлен вариант ШШС-метода. Такая схема сочетает в себе БС аппроксимацию решения по пространству, специальный ТУБ алгоритм метода Рунге-Кутты для решения возникающей системы ОДУ и модифицированную процедуру построения ограничителя. Для результирующей явной схемы с кусочно-линейной пространственной аппроксимацией решения доказана устойчивость при Ссрь ^ | в линейном случае и показано, что схема имеет второй порядок аппроксимации по пространству и по времени, в том числе в точках локального экстремума, а также доказано ТУВМ свойство (ограниченность полной вариации средних значений решения по ячейке). Проведенный численный анализ свидетельствует о сходимости метода к решению, удо-влетворющему условию неубывания энтропии, даже для невыпуклых функций потокаразрывы решения передаются на 2−3 вычислительных ячейках, а решение в областях гладкости и на локальных экстремумах сходится со вторым порядком.

В работе [35] (Cockburn, Shu) данный подход формально обобщен для получения RKDG схем высокого порядка аппроксимации скалярных законов сохранения. При построении RKDG-метода порядка (к +1) использована кусочная пространственная аппроксимация решения полиномами степени не выше к, TVD вариант явного метода Рунге-Кутты порядка (к + 1) для дискретизации по времени и обобщенная процедура монотонизации решения. Схемы такого типа применены в работе [36] (Cockburn, Lin, Shu) при решении одномерных систем, а обобщение RKDG-метода на случай многомерных систем уравнений впервые рассмотрено в работе [37] (Cockburn, Hou,.

Shu). Преимущество использования RKDG-метода при решении многомерных систем состоит в том, что он позволяет работать со сложной геометрией расчетной области и использовать при вычислениях неструктурированные сетки для пространственной DG-дискретизации, при этом временная TVD дискретизация остается такой же, как и в случае одномерных уравнений и систем. Единственной проблемой при разработке схемы RKDG для многомерных задач остается построение ограничителя, так как в этой ситуации условия устойчивости становятся более жесткими по сравнению с одномерным случаем. В отличие от одномерных задач, для которых возможно создание TVD схем высокого порядка, порядок многомерных TVD схем на неструктурированных сетках не превышает первого [18,38]. Поэтому использование любого TVD ограничителя неизбежно будет приводить к потере порядка результирующей схемы. В [37] разработан вариант обобщенного ограничителя, основанный на требовании выполнения локального принципа максимума для численных схем. В работах [19−24,39−41] область применения RKDG-метода расширена для решения двухи трехмерных задач газовой динамики.

В последнее время наблюдается возрастающий интерес к RKDG-методу в связи с необходимостью численного решения задач динамики многофазных сред, описываемых неконсервативными гиперболическими системами уравнений. В качестве математической модели двухфазных сред наиболее часто используются уравнения Баера-Нунциато, которые в трехмерном случае представляют собой систему одиннадцати уравнений в частных производных первого порядка и описывают динамику смеси, состоящей, как правило, из твердой и газообразной фазы. Такая модель впервые предложена Ваег и Nunziato [42] при описании динамики гранулированных горючих материалов в газообразных продуктах горения. Отличительной чертой уравнений Баера-Нунциато является то, что они предполагают наличие двух векторов скорости и двух давлений (в соответствии с числом фаз). Другие распространенные модели динамики многофазных сред имеют вид, сходный с уравнениями Баера-Нунциато [43−47], поэтому с практической точки зрения целесообразно исследовать эти уравнения, сохраняя при этом некоторую степень общности. Математические особенности одномерных однородных уравнений Баера-Нунциато исследованы в [48] (Embid, Ваег). Эти уравнения являются гиперболическими, известна их характеристическая структура. Однако уравнения не могут быть записаны в форме законов сохранения (дивергентной форме): в консервативных переменных уравнения сохраняют неконсервативные члены, что делает невозможным применение классических условий Гю-гонио на разрывах решения. Решение неконсервативных уравнений и систем долгое время оставалось сложной проблемой как с теоретической, так и с вычислительной точек зрения [49−54].

В настоящей работе преимущественно исследована задача Римана для уравнений Баера-Нунциато [47,55,56], разработан численный метод ее решения и численный поток типа HLLC, который впоследствии использован в методе конечных объемов, RKDG-методе и PC-методе (англ. path-conservative) [43,44,46,56−61].

Цель и задачи исследования

Целью работы является развитие RKDG-метода, его оптимизация, разработка и применение основанного на нем программного комплекса для математического моделирования процессов газовой динамики и динамики двухфазных сред.

Для достижения поставленной цели потребовалось решение следующих основных задач.

1. Реализация алгоритма RKDG-метода для решения одномерного квазилинейного уравнения переноса pi создание основы программного комплекса.

Проверка реализованного алгоритма и программного комплекса с использованием как известных в литературе, так и специально созданных тестовых задач.

2. Реализация алгоритма ШШС-метода для решения двумерных уравнений Эйлера и создание программного комплекса для численного решения задач идеальной газовой динамики. Проверка реализованного алгоритма и программного комплекса на примере известных тестовых задач и сравнение различных вариантов ШШС-метода с другими известными методами.

3. Реализация алгоритма ШШС-метода для решения двумерных уравнений Навье-Стокса для сжимаемого газа и расширение программного комплекса для учета вязкости и теплопроводности среды. Тестирование реализованного алгоритма на задачах обтекания различных профилей и сравнение полученных в результате вычислительного эксперимента аэродинамических коэффициентов с экспериментальными данными.

4. Оптимизация параметров алгоритма монотонизации численного решения и анализ эффективности монотонизации как элемента ШШС-метода.

5. Разработка и программная реализация параллельного алгоритма ШШС-метода, анализ его эффективности на различных типах вычислительных систем.

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

Методы исследования. Основным методом исследования задач, поставленных в диссертационной работе, является вычислительный эксперимент.

Достоверность и обоснованность полученных результатов гарантируется строгостью используемого математического аппарата и подтверждается сравнением результатов численного моделирования с известными экспериI ментальными данными, а также данными вычислительных экспериментов, выполненных известными численными методами.

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

Для уравнений Эйлера, описывающих динамику идеального газа, проведен детальный сравнительный анализ ЯКОС-метода с другими известными методами, такими как метод конечных объемов с численными потоками году-новского типа, а именно: потоками Куранта-Изаксона-Риса (КИР), Лакса-Фридрихса и потоками типа Хартена-Лакса-ван Лира (НЬЬ и НЪЪС) [18]. ШФС-метод является одним из обобщений конечно-объемных методов на случай кусочно-полиномиальной аппроксимации решения. Численные потоки выбираются аналогично тем, которые используются в конечно-объемных методах. Шаблон аппроксимации в ШШО-методе состоит из ячейки и ее ближайших соседей, а порядок аппроксимации зависит лишь от порядка полинома базисной функции ячейки. Это дает возможность локально изменять порядок аппроксимации и адаптировать ее к особенностям решения. Г1КВС-метод позволяет работать с геометрически сложными областями, достигая при этом высокой точности. Благодаря узкому шаблону ШФС-метода возможно использование неструктурированных сеток. В методе достаточно просто реализуются параллельные вычисления.

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

Исследована эффективность применения ШФС-метода для решения системы уравнений динамики вязкого теплопроводного газа. В этом случае градиенты неизвестных величин рассматриваются как независимые переменные, в результате чего система сводится к системе дифференциальных уравнений в частных производных первого порядка, после чего применяется стандартная схема ККБО-метода. Невязкие потоки через границу ячейки аппроксимируются так же, как и в случае идеального газа, а для вязких потоков используется центральный численный поток.

Разработана модификация ШШС-метода, позволяющая автоматически оптимизировать параметры монотонизатора решения в процессе расчетов, адаптируя их к локальным особенностям решения. Проведено сравнение результатов, полученных с использованием предложенного алгоритма монотонизации и стандартного алгоритма, использующего фиксированные параметры ограничителя.

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

В работе изучена возможность расширения области применения RKDG-метода для решения неконсервативных гиперболических систем. В качестве примера исследована система уравнений Баера-Нунциато, описывающая движения двухфазной среды без учета фазовых переходов. Разработан численный поток типа HLLC, основанный на аппроксимации решения задачи Рима-на для уравнений Баера-Нунциато, который затем использован при построении численной схемы метода конечных объемов, RKDG-метода и РС-метода. Результирующие схемы проверены на специально подобранных тестовых задачах, сделаны соответствующие выводы о работоспособности RKDG-метода по сравнению с другими методами решения неконсервативных систем.

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

На защиту выносятся следующие положения.

1. Оптимизированный алгоритм RKDG-метода с возможностью пространственно-временной адаптации монотонизатора к особенностям решения.

2. Алгоритм численного решения задачи Римана для уравнений Баера-Нунциато. Применение численного потока на его основе в методах конечных объемов, PC и RKDG.

3. Параллельный алгоритм RKDG-метода.

4. Применение программного комплекса на основе RKDG-метода для расчета течений жидкости и газа, определения аэродинамических характеристик профилей и моделирования динамики двухфазных сред.

Апробация работы. Результаты диссертационной работы апробированы на ХУ1-Й Школе-семинаре молодых ученых и специалистов под руководством академика РАН А. И. Леонтьева (Санкт-Петербург, 2007), Международной конференции «Дифференциальные уравнения. Функциональные пространства. Теория приближений» (Новосибирск, 2008), XII Международной конференции «Современные проблемы механики сплошной среды» (Ростов-на-Дону, 2008), 5-й Всероссийской конференции «Необратимые процессы в природе и технике» (Москва, 2009), Международной конференции по вычислительному теплои массообмену (Гуанчжоу, Китай, 2009), 17-й Международной конференции по математической физике (Прага, Чехия, 2009), 6-й Международной конференции по вычислительной газовой динамике (Санкт-Петербург, 2010).

Публикации. Основные результаты диссертационной работы опубликованы в 2 препринтах [62,63], 6 научных статьях [24,60,64−67], в том числе в 5 статьях из Перечня рецензируемых ведущих научных журналов и изданий [24,60,64−66], и 10 тезисах и докладах конференций [68−77].

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

Структура и объём работы. Диссертационная работа состоит из введения, трех глав, общих выводов и списка литературы. Работа изложена на 141 страницах, содержит 96 иллюстраций и 24 таблицы.

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

включает 110 наименований.

Выводы.

По результатам проведенных в работе исследований могут быть сделаны следующие выводы.

1. Разработан оптимизированный алгоритм ИКОС-метода, позволяющий проводить процедуру монотонизации с автоматическим выбором параметров лимитера в процессе решения задачи с возможностью пространственно-временной адаптации лимитера к особенностям решения. Разработан алгоритм численного решения задачи Римана для уравнений Баера-Нунциато, на его основе создан алгоритм вычисления потока, примененный в качестве базового элемента в методе конечных объемов, РС-методе и ШФС-методе.

2. Создан программный комплекс на основе ККБС-метода, позволяющий проводить расчеты течений как идеального, так и вязкого теплопроводного газа или жидкости, а также моделировать течения двухфазных сред с помощью метода конечных объемов, РС-метода и ШШС-метода. Разработан параллельный алгоритм ШШС-метода и проанализирована эффективность распараллеливания на различных вычислительных комплексах.

3. Продемонстрирована работоспособность ЯКБС-метода и основанного на нем программного комплекса на задачах моделирования течений газа в каналах сложной формы, а также задачах определения аэродинамических характеристик различных профилей, включая крыловые профили с закрылками. Эффективность созданного алгоритма вычисления потока типа НЬЬС при его использовании в ККБС-методе решения двухфазных уравнений Бае-ра-Нунциато проверена на специально подобранных задачах Римана и задаче взаимодействия материалов.

Показать весь текст

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

  1. В. С., Кувыркин Г. Н. Математические модели термомеханики. М.: Физматлит, 2002. 168 с.
  2. Лойцянский J1. Г. Механика жидкости и газа. М.: Дрофа, 2003. 840 с.
  3. Г. Г. Газовая динамика. М.: Наука, 1988. 424 с.
  4. П. К. Отрывные течения.: В 3 т. М.: Мир, 1972. Т. 1. 1972. 299 с. Т. 2. 1973. 279 с. Т. 3. 1973. 333 с.
  5. Н. Ф., Кошевой В. Н., Калугин В. Т. Аэродинамика отрывных течений. М.: Высшая школа, 1988. 351 с.
  6. Н. Ф., Кошевой В. Н., Данилов А. Н. и др. Прикладная аэродинамика. М.: Высшая школа, 1974. 731 с.
  7. Н. Ф., Кошевой В. Н., Захарченко В. Ф., Данилов А. И. Основы прикладной аэродинамики. Кн. 2. Обтекание тел вязкой жидкостью. Рулевые устройства. М.: Высшая школа, 1991. 351 с.
  8. M. Е. Техническая газодинамика (основы газодинамики турбин). M.-JI.: Энегроиздат, 1953. 544 с.
  9. JI. Гидроаэромеханика. Ижевск: НИЦ «Регулярная и хаотическая динамика», 2000. 576 с.
  10. Л. Д., Лифшиц Е. М. Теоретическая физика: Учебное пособие.: В 10 т. T. VI. Гидродинамика. М.: Наука, 1986. 736 с.
  11. Дж. Введение в динамику жидкости. М.: Мир, 1973. 779 с.
  12. Patankar S. V. Numerical heat transfer and fluid flow. Hemisphere Publishing Corporation, 1980. 196 p.
  13. A. A., Гулин A. В. Численные методы. M.: Наука, 1989. 432 с.
  14. Того Е. F. Riemann Solvers and Numerical Methods for Fluid Dynamics. Berlin: Springer, 2009. 724 p.
  15. . Л., Яненко H. H. Системы квазилинейных уравнений и их приложения к газовой динамике. М.: Наука, 1978. 687 с.
  16. Дж. Линейные и нелинейные волны. М.: Мир, 1977. 624 с.
  17. А. Г., Погорелов Н. В., Семенов А. Ю. Математическиевопросы численного решения гиперболических систем уравнений. М.: Физматлит, 2001. 608 с.
  18. Cockburn В., Shu С. W. Runge-Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems // J. Sci. Сотр. 2001. V. 3. P. 173−261.
  19. Bassi F., Rebay S. A high-order accurate discontinuous finite element solution of the 2D Euler equations // J. Сотр. Phys. 1997. V. 138. P. 251 285.
  20. Bassi F., Rebay S. A high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier-Stokes equations // J. Сотр. Phys. 1997. V. 131. P. 267−279.
  21. Dolejsi V. On the discontinuous Galerkin method for the numerical solution of the Navier-Stokes equations // Int. J. Numer. Meth. Fluids. 2004. V. 45. P. 1083−1106.
  22. Van der Vegt J. J. W., Van der Ven H. Discontinuous Galerkin finite element method with anisotropic local grid refinement for inviscid compressible flows // J. Сотр. Phys. 1998. V. 141. P. 46−77.
  23. M. П., Савенков E. В., Токарева С. А. Решение задач газовой динамики с ударными волнами RKDG-методом // Математическое моделирование. 2008. Т. 20, № 11. С. 55−66.
  24. W. Н., Hill Т. R. Triangular mesh methods for the neutron transport equation // Los Alamos Schientific Laboratory Report LA-UR-73−479. 1973.
  25. LeSaint P., Raviart P. A. On a finite element method for solving the neutron transport equation // Mathematical aspects of finite elements in partial differential equations (C. de Boor, Ed.). Academic Press. 1974. P. 89−145.
  26. Johnson C., Pitkaranta J. An analysis of the discontinuous Galerkin method for a scalar hyperbolic equation // Math. Сотр. 1986. V. 46. P. 1−26.
  27. Peterson T. A note on the convergence of the discontinuous Galerkin method for a scalar hyperbolic equation // SIAM J. Numer. Anal. 1991. V. 28. P. 133 140.
  28. Richter G. R. An optimal-order error estimate for the discontinuous Galerkin method // Math. Сотр. 1988. V. 50. P. 50−75.
  29. Chavent G., Salzano G. A finite element method for the Id water flooding problem with gravity //J. Сотр. Phys. 1982. V. 45. P. 307−344.
  30. Chavent G., Cockburn B. The local projection ^^-discontinuous Galerkinfinite element method for scalar conservation laws // M2AN. 1989. V. 23. P. 565−592.
  31. Van Leer B. Towards the ultimate conservation difference scheme, ii // J. Сотр. Phys. 1974. V. 14. P. 361−376.
  32. А. С. Численные методы решения уравнений и систем гиперболического типа // Энциклопедия низкотемпературной плазмы. Часть 2. 2008. С. 220−235.
  33. Cockburn В., Shu С. W. The Runge-Kutta local projection-discontinuous Galerkin method for scalar conservation laws // M2AN. 1991. V. 25. P. 337 361.
  34. Cockburn В., Shu C. W. Tvb Runge-Kutta local projection discontinuous Galerkin finite element method for scalar conservation laws ii: General framework // Math. Сотр. 1989. V. 52. P. 411−435.
  35. Cockburn В., Lin S. Y., Shu C. W. Tvb Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws iii: One-dimensional systems // J. Сотр. Phys. 1989. V. 84. P. 90−113.
  36. Cockburn В., Hou S., Shu C. W. Tvb Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws iv: The multidimensional case // Math. Сотр. 1990. V. 54. P. 545−581.
  37. Goodman J., LeVeque R. On the accuracy of stable schemes for 2d scalar conservation laws // Math. Сотр. 1985. V. 45. P. 15−21.
  38. Cockburn В., Shu C. W. The p*-Rkdg method for two-dimensional Euler equations of gas dynamics // ICASE Report. 1991. № 91−32.
  39. Cockburn В., Shu C. W. The Runge-Kutta discontinuous Galerkin finite element method for conservation laws v: Multidimensional systems // J. Сотр. Phys. 1998. V. 141. P. 199−224.
  40. С. В. Разработка и исследование численных схем высокого порядка точности для решения уравнений газовой динамики на неструктурированных сетках : Дис.. д. ф.-м. н. М.: 2007. 127 с.
  41. Baer М. R., Nunziato J. W. A Two-Phase Mixture Theory for the Deflagration-to-Detonation Transition (DDT) in Reactive Granular Materials //J. Multiphase Flow. 1986. V. 12. P. 861−889.
  42. Saurel R., Abgrall R. A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows //J. Сотр. Phys. 1999. V. 150. P. 425−467.
  43. Saurel R., Abgrall R. Discrete Equations for Physical and Numerical Compressible Multiphase Mixtures // J. Comp. Phys. 2003. V. 186. P. 361 396.
  44. Stewart H. B., Wendroff B. Two-Phase Flow: Models and Methods // J. Comp. Phys. 1984. V. 56. P. 363−409.
  45. Romenski E., Resnyanski E. D., Toro E. F. Conservative Hyperbolic Formulation for Compressible Two-Phase Flow with Different Phase Pressures and Temperatures // Quarterly of Applied Mathematics. 2007. V. 65. P. 259−279.
  46. Deledicque V., Papalexandris M. V. An Exact Riemann Solver for Compressible Two-Phase Flow Models Containing Non-Conservative Products //J. Comp. Phys. 2007. V. 222. P. 217−245.
  47. Embid P., Baer M. Mathematical Analysis of a Two-Phase Continuum Mixture Theory // Continuum Mech. Thermodyn. 1992. V. 4. P. 279−312.
  48. Toro E. F. Riemann-Problem Based Techniques for Computing Reactive Two-Phase Flows // Lecture Notes in Physics. 1989. V. 351. P. 472−481.
  49. Toumi I. An Upwind Numerical Method for Two-Fluid Two-Phase Flow Models // Nuclear Sci. Eng. 1996. V. 123. P. 147−168.
  50. Saurel R., Abgrall R. A Simple Method for Compressible Multifluid Flows // SIAM J. Sci. Comput. 1999. V. 21, № 3. P. 1115−1145.
  51. Pailliere H., Kumbaro A., Bestion D. et al. Advanced Three-Dimensional Two-Phase Flow Simulation Tools for Application to Reactor Safety (ASTAR) // Nuclear Engineering and Design. 2005. V. 235. P. 379−400.
  52. Chang C. H., Liou M. S. A Robust and Accurate Approach to Computing Compressible Multiphase Flow: Stratified Flow Model and AUSM±up Scheme // J. Comp. Phys. 2007. V. 225. P. 840−873.
  53. Saurel R., Le Metayer O., Massoni J., Gavrilyuk S. Shock jump relations for multiphase mixtures with stiff mechanical relaxation // Shock Waves. 2007. V. 16. P. 209−232.
  54. Andrianov N., Warnecke G. The Riemann Problem for the Baer-Nunziato Model of Two-Phase Flows // J. Comp. Phys. 2004. V. 195. P. 434−464.
  55. Schwendeman D. W., Wahle C. W., Kapila A. K. The Riemann Problem and a High-Resolution Godunov Method for a Model of Compressible Two-Phase Flow //J. Comp. Phys. 2006. V. 212. P. 490−526.
  56. Rhebergen S., Bokhove О., Van der Vegt J. J. W. Discontinuous Galerkin Finite Element Methods for Hyperbolic Nonconservative Partial Differential Equations //J. Comput. Phys. 2008. V. 227. P. 1887−1922.
  57. Pares С. Numerical Methods for Nonconservative Hyperbolic Systems: a theoretical framework // SIAM J. Numer. Anal. 2006. V. 44. P. 300−321.
  58. Luz Munoz-Ruiz M., Pares С. Godunov Method for Nonconservative Hyperbolic Systems // Modelisation mathematique et analyse numerique.2007. V. 41. P. 169−185.
  59. S.A., Того E. F. HLLC-type Riemann solver for the Baer-Nunziato equations of compressible two-phase flow //J. Сотр. Phys. 2010. V. 229. P. 3573−3604.
  60. M., Hidalgo A., Castro M., Pares С., Того E. FORCE Schemes on Unstructured Meshes II: Nonconservative Hyperbolic Systems // Сотр. Meth. Appl. Mech. Eng. 2010. V. 199. P. 625−647
  61. M. П., Савенков E. Б., Токарева С. А. Применение разрывного метода Галеркина для численного решения квазилинейного уравнения переноса. М., 2005. 34 с. (Препринт ИПМ им. М. В. Келдыша РАН, № 105).
  62. Применение RKDG метода для численного решения задач газовой динамики / С. А. Токарева и др. М., 2006. 30 с. (Препринт ИПМ им. М. В. Келдыша РАН, № 52).
  63. Galanin M. P., Tokareva S. A. The RKDG method and its application for the numerical solution of gas dynamics problems // Heat Transfer Research.2008. V. 39, № 2. P. 123−132.
  64. M. P., Savenkov E. В., Tokareva S. A. Solving gas dynamics problems with shock waves using the Runge-Kutta discontinuous Galerkin method // Mathematical Models and Computer Simulations. 2009. V. 5, № 1. P. 635−645.
  65. И. К., Токарева С. А. Сравнение эффективности параллельных алгоритмов решения задач газовой динамики на разных вычислительных комплексах // Вестник МГТУ им. Н. Э. Баумана Естественные науки. 2009. № 1. С. 90−97.
  66. Tokareva S. A. A problem-independent slope limiting algorithm for the Runge-Kutta discontinuous Galerkin method // Computational Methods in
  67. Applied Mathematics. 2010. V. 10, № 3. P. 326−342.
  68. С. А. Применение разрывного метода Галеркина для численного решения задач газовой динамики // Студенческий научный Вестник: Тез. докл. общеуниверситетской научно-технической конференции «Студенческая весна—2006». М., 2006. Т. 3. С. 145.
  69. С. А. RKDG-метод для численного решения задач динамики вязкого газа // Студенческий научный Вестник: Тез. докл. общеуниверситетской научно-технической конференции «Студенческая весна — 2007». М., 2007. Т. 4. С. 86−87.
  70. С. А. Параллельный алгоритм RKDG-метода для расчета заполнения газом канала ускорителя // Студенческий научный Вестник: тез. докл. общеуниверситетской научно-технической конференции «Студенческая весна —2008». М., 2008. Т. 6., ч. 2. С. 25−26.
  71. С. А. RKDG-метод и его применение для численного решения задач газовой динамики // Необратимые процессы в природе и технике: Труды пятой Всероссийской конференции. Москва, 2009. Ч. 2. С. 93−96.
  72. И. К., Токарева С. А. Параллельный алгоритм RKDG-метода для решения задач газовой динамики / / Современные проблемы механики сплошной среды. Ростов-на-Дону, 2008. С. 150−154.
  73. Tokareva S. A. Runge-Kutta Discontinuous Galerkin method for simulation of flow around obstacles // Proceedings of the 6th International Conference of Computational Heat and Mass Transfer. Guangzhou, China, 2009. P. Ill— 116.
  74. Tokareva S. A. Numerical simulation of fluid flow with Discontinuous
  75. Galerkin methods // Proceedings of 16th ICMP, Prague, Czech Republic, 2009.
  76. S. А., Того E. F. HLLC-type Riemann solver for the Baer-Nunziato equations of compressible two-phase flow // Proceedings of the 6th ICCFD, St. Petersburg, Russia, 2010. P. 238−239.
  77. Т. А., Галанин M. П. Нелинейная монотонизация схемы К. И. Бабенко для численного решения квазилинейного уравнения переноса. М., 2003. 35 с. (Препринт ИПМ им. М. В. Келдыша РАН, № 62).
  78. Shu С. W. TVB uniformly high order schemes for conservation laws // Math. Comp. 1987. V. 49. P. 105−121.
  79. Shu C. W., Osher S. Efficient implementation of essentially non-oscillatory shock-capturing schemes // J. Comput. Phys. 1988. V. 77. P. 439−471.
  80. Shu C. W., Osher S. Efficient implementation of essentially non-oscillatory shock-capturing schemes, ii //J. Comput. Phys. 1989. V. 83. P. 32−78.
  81. Shu C. W. TVD time discretizations // SIAM J. Sei. Stat. Comput. 1988. V. 9. P. 1073−1084.
  82. А. А., Попов Ю. П. Разностные методы решения задач газовой динамики. М.: Едиториал УРСС, 2004. 424 с.
  83. О. А. Математическое моделирование двойных звездных систем : Дис.. д. ф.-м. н. М.: 1999. 323 с.
  84. Woodward P., Colella P. The numerical simulation of two-dimensional fluid flow with strong shocks // J. Comp. Phys. 1984. V. 54. P. 115−173.
  85. О. Метод конечных элементов в технике. М.: Мир, 1975. 541 с.
  86. Li В. Q. Discontinuous finite elements in fluid dynamics and heat transfer. Berlin: Springer, 2006. 578 p.
  87. Bassi F., De Bartolo С., Hartmann R., Nigro A. A discontinuous Galerkin method for inviscid low Mach number flows //J. Comp. Phys. 2009. V. 228. P. 3996−4001.
  88. Th. von Karman. Uber den Mechanismus des Widerstandes, den ein bewegter Korper in einer Flussigkeit erfahrt // Gottinger Nachrichten, mathematischphysikalische Klasse. 1912. P. 547−556.
  89. Zahm A. F. Flow and drag formulas for simple quadrics // Aerodynamical Laboratory, U. S. Navy, 1927. Report № 253.
  90. Случановская 3. П. Распределение давления на поверхности прямоугольного, трехгранного и полукруглого цилиндров и их аэродинамические коэффициенты //Тр. Инс-та механики МГУ, № 24 / Под ред. С. М. Горлина. М.: Изд-во МГУ, 1973. С. 52−60.
  91. Novak М., Tanaka Н. Effect of Turbulence on Galloping Instability // Journal of the Engineering Mechanics Division, ASCE. 1984. V. 100. P. 27−47.
  92. Novak M. Aeroelastic Galloping of Prismatic Bodies // Journal of the Engineering Mechanics Division, ASCE. 1969. V. 95. P. 115−142.
  93. Novak M. Galloping Oscillations of Prismatic Structures // Journal of the Engineering Mechanics Division, ASCE. 1972. V. 98. P. 27−46.
  94. Novak M., Davenport A.G. Aeroelastic Instability of Prisms in Turbulent Flow // Journal of the Engineering Mechanics Division, ASCE. 1970. V. 96. P. 17−39
  95. Parkinson G. V., Smith J. D. The square prism as an aeroelastic nonlinear oscillator // Quart. J. Mech. and Appl. Math. 1964. V. 17. P. 255.
  96. Lyn D. A., Einav S., Rodi W., Park J .H. A laser-Doppler velocimetry study of ensemble-averaged characteristics of the turbulent wake of a square cylinder // J. Fluid Mech. 1995. V. 304. P. 285−319.
  97. Dolejsi V. Semi-implicit interior penalty discontinuous Galerkin methods for viscous compressible flows // Commun. in Comput. Phys. 2008. V. 4. P. 231 274.
  98. Ван Дайк M. Альбом течений жидкости и газа. М.: Мир, 1986. 184 с.
  99. В. Д. Параллельное программирование в MPI. Новосибирск: Изд-во ИВМиМГ СО РАН, 2002. 215 с.
  100. В. П. Теория и практика параллельных вычислений. М.: БИНОМ. Лаборатория знаний, 2007. 423 с.
  101. С. А., Посыпкин М. А. Технологии параллельного программирования. М.: Форум Инфра-М, 2008. 208 с.
  102. Wenzinger С. J., Harris Т. A. Wind-tunnel investigation of an N.A.C.A. 23 102 airfoil with various arrangements of slotted flaps // National Advisory Committee for Aeronautics, 1939. Report № 664. vvi.p.
  103. Того E. F., Spruce M., Speares W. Restoration of the Contact Surface in the HLL-Riemann Solver // Technical Report CoA-9204, Department of Aerospace Science, College of Aeronautics, Cranfield Institute of Technology, 1. UK, 1992.
  104. Toro E. F., Spruce M., Speares W. Restoration of the Contact Surface in the HLL-Riemann Solver // Shock Waves 1994. V. 4. R 25−34.
  105. Batten P., Clarke N., Lambert C., Causon D. On the Choice of Wave Speeds in the HLLC Riemann Solver // SIAM J. Sci. and Stat. Comp. 1997. V. 18. P. 1553−1570.
  106. Batten P., Leschziner M. A., Golberg U. C. Average-State Jacobians and Implicit Methods for Compressible Viscous and Turbulent Flows // J. Comput. Phys. 1997. V. 137. P. 38−78.
  107. Volpert A. I. The Space BV and Quasilinear Equations // Math. USSR Sbornik. 1967. V. 73. P. 225−267.
  108. Dal Maso G., LeFloch P. G., Murat F. Definition and Weak Stability of Nonconservative Products //J. Math. Pures Appl. 1995. V. 74. P. 483.
  109. Toumi I. A Weak Formulation of Roe’s Approximate Riemann Solver // J. Comput. Phys. 1992. V. 102. P. 360−373.
Заполнить форму текущей работой

Сессия? Спокойно!