- •3. Разностные методы решения уравнений газовой динамики.
- •3.1 Особенности численного решения задач газовой динамики.
- •3.2 Способы описания газодинамических течений и построение разностных схем.
- •3.3 Однородные разностные схемы. Схемы с псевдовязкостью.
- •3.3.1. Способы единообразного описания газодинамических течений.
- •3.4. Разностная схема Неймана-Рихтмайера - «крест» для системы уравнений газовой динамики с вязкостью.
3.2 Способы описания газодинамических течений и построение разностных схем.
Характер применяемых схем интегрирования уравнений газовой динамики существенно зависит от способа описания течения. Выделим три наиболее важных из них.
1-й
способ.
Область
плоскости
,
в которой рассматривается движение,
разбивается сильными и слабыми разрывами
на отдельные области
гладкого течения, в которых удовлетворяются
уравнения газовой динамики, в то время
как на самих разрывах удовлетворяются
условия совместности. В
таком рассмотрении обобщенное решение
есть совокупность гладких решений,
определяемых в областях
и примыкающих друг к другу через линии
разрывов с соблюдением условий
совместности. При
таком описании возникает необходимость
численного интегрирования уравнений
газовой динамики в областях
с выполнением условий примыкания на
линиях разрыва.
Наиболее известным разностным методом, соответствующим этому способу описания, является метод характеристик. Действительно, среди линий раздела мы имеем слабые разрывы и контактные границы, являющиеся характеристиками, что и делает удобной характеристическую разностную схему.
Полная детальность описания течения, составляющая положительную черту метода характеристик, затрудняет одновременно его реализацию на ЭВМ из-за сложной логики расчета особенностей и построения фронта расчета.
Разумеется, что метод характеристик (мы его подробно рассмотрим в последующих лекциях) не является единственным разностным методом, который может быть применен в рамках детального описания течения газа.
2-й способ. Обобщенное решение определяется интегральными законами сохранения в эйлеровых или лагранжевых координатах. Такое описание является единообразным, поскольку как уравнения газовой динамики, так и условия совместности являются следствиями законов сохранения.
Разностные схемы, соответствующие второму способу описания, получаются единообразной аппроксимацией законов сохранения независимо от характера течения и поэтому носят название однородных схем или схем сквозного счета.
3-й способ. Обобщенное решение определяется как предел классического решения некоторой системы квазилинейных параболических уравнений с малыми параметрами при старших производных.
Если
(3.1)
есть исходная система уравнений газовой динамики, записанная в виде законов сохранения, то соответствующая параболическая система имеет вид:
(3.2)
Здесь
-
вектор-функция, описывающая течение,
,
- векторные функции
от векторного аргумента
,
- квадратная матрица,
- малый параметр.
Матрица
должна быть подобрана таким образом,
чтобы решение
системы (3.2) обладало
достаточной гладкостью и при
приближалось в
каком-то смысле к решению системы (3.1).
Разностные схемы, основанные на третьем способе рассмотрения, также имеет характер сквозного счета. В некоторых случаях второй и третий подходы приводят к одинаковым схемам.
Сделаем несколько замечаний по поводу граничных условий в задачах газовой динамике.
Постановка граничных условий зависит как от характера рассматриваемой задачи, так и от способа описания течения и в первую очередь от принятой системы координат. Если перечислить все достоинства и недостатки аппроксимации граничных условий в разных системах описания – эйлеровых и лагранжевых, то становится понятным, что универсального и наилучшего способа разностного решения задач газовой динамики нет. В реальных расчетах приходится приспосабливаться к конкретной задаче и применять различные способы.
