Студопедия — Разностные схемы для системы модельных уравнений
Студопедия Главная Случайная страница Обратная связь

Разделы: Автомобили Астрономия Биология География Дом и сад Другие языки Другое Информатика История Культура Литература Логика Математика Медицина Металлургия Механика Образование Охрана труда Педагогика Политика Право Психология Религия Риторика Социология Спорт Строительство Технология Туризм Физика Философия Финансы Химия Черчение Экология Экономика Электроника

Разностные схемы для системы модельных уравнений






Среди рассмотренных ранее упрощенных моделей полной системы уравнений Навье-Стокса система уравнений газовой динамики в случае баротропного течения была записана в следующем операторном виде:

, (27)

где .

Рассмотрим наиболее используемые разностные схемы для решения этой системы уравнений.

Одной из наиболее известных разностных схем решения системы (27) является явная схема Маккормака (см.[16]):

(28)

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

Найдем условие устойчивости для системы (28). Для этого преобразуем ее в схему с целыми временными шагами. Из первого уравнения системы (28) находим и, подставляя во второе уравнение, имеем:

или

= .

Исследуем устойчивость схемы (28) для решений вида: .

Подставляя волновой вектор в последнее уравнение, получим:

Теперь можно выписать характеристическое уравнение для системы (28). Обозначая , имеем:

.

Раскрывая определитель, находим корни характеристического уравнения:

После соответствующих преобразований первое уравнение равно: и

Аналогично для второго уравнения: . Тогда – условие устойчивости. Данная схема условно устойчива.

Предложенная схема не является экономичной из-за ограничений на временной шаг. Можно рассмотреть неявные схемы, которые не обладают этим недостатком.

Схема с весами

Для системы уравнений (27) теперь рассмотрим разностную схему с весами

,

имеющую порядок аппроксимации Здесь , при и в остальных случаях. Проведем спектральный анализ устойчивости схемы с весами. В этом случае получаем следующее характеристическое уравнение:

, где .

Схема безусловно устойчива при , что следует из вида характеристического уравнения:

.

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

. Таким образом, данная схема является неэкономичной.

Рассмотрим теперь аппроксимацию частных производных в уравнениях (27) несимметричными разностными операторами. В случае несимметричных схем аппроксимацию конвективных членов следует согласовывать с аппроксимацией члена с давлением. Например чисто неявная схема (здесь разностный оператор выбирается в зависимости от знака скорости )

будет условно устойчива при < .

Для получения безусловно устойчивой схемы необходимо слагаемые, отвечающие градиенту давления, аппроксимировать сопряженным разностным оператором по отношению к конвективным слагаемым, например:

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

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

Рассмотрим расщепление по физическим процессам для системы уравнений: . Представим матрицу в виде суммы двух матриц, где . Матрица содержит слагаемые, связанные с конвективными членами уравнения, а матрица – аналог члена с давлением в уравнении движения и члена вида в уравнении неразрывности. Расщепление матрицы А порождает расщепление дифференциальных операторов на

и как следствие расщепление разностных операторов на

.

Одним из вариантов реализации схемы расщепления есть следующая схема в дробных шагах:

(29)

Рассмотрим более подробно реализацию схемы на дробных шагах, используя для аппроксимации центральные разности.

На первом дробном шаге решается следующая система:

Так как матрица – диагональная, то уравнения можно решать скалярными прогонками независимо друг от друга, а для даже по методу бегущего счета.

На втором дробном шаге решается система:

Из первого уравнения системы выражаем и подставляем во второе уравнение для нахождения параметра , т.е.

Последнее уравнение решается методов пятиточечной прогонки.

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

Для первого уравнения системы (29) в случае симметричной аппроксимации получаем следующее характеристическое уравнение: . Раскрывая определитель, имеем

Аналогично для второго дробного шага получаем

. Все и следовательно их произведение тоже по модулю меньше единицы.

Схема безусловно устойчива.

Схемы приближенной факторизации

Для получения схем второго порядка аппроксимации можно использовать метод приближенной факторизации. Рассмотрим его на примере схемы с весами для уравнения (27):

или в канонической форме: .

Используя расщепление матрицы А по физическим процессам, перепишем уравнение в виде:

.

Факторизация стабилизирующего оператора состоит в следующем:

. Тогда схема вида и называется схемой приближенной факторизации. Она может быть реализована в дробных шагах:

если , то

где .

Данная схема является экономичной и аппроксимирует исходную систему с тем же порядком, что и схема с весами.

Реализация схемы на дробных шагах аналогична ранее рассмотренной схемы.

Схема безусловно устойчива при и корни характеристического уравнения в случае симметричной аппроксимации равны:

.

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







Дата добавления: 2014-11-10; просмотров: 1069. Нарушение авторских прав; Мы поможем в написании вашей работы!



Шрифт зодчего Шрифт зодчего состоит из прописных (заглавных), строчных букв и цифр...

Картограммы и картодиаграммы Картограммы и картодиаграммы применяются для изображения географической характеристики изучаемых явлений...

Практические расчеты на срез и смятие При изучении темы обратите внимание на основные расчетные предпосылки и условности расчета...

Функция спроса населения на данный товар Функция спроса населения на данный товар: Qd=7-Р. Функция предложения: Qs= -5+2Р,где...

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

Билет №7 (1 вопрос) Язык как средство общения и форма существования национальной культуры. Русский литературный язык как нормированная и обработанная форма общенародного языка Важнейшая функция языка - коммуникативная функция, т.е. функция общения Язык представлен в двух своих разновидностях...

Патристика и схоластика как этап в средневековой философии Основной задачей теологии является толкование Священного писания, доказательство существования Бога и формулировка догматов Церкви...

Схема рефлекторной дуги условного слюноотделительного рефлекса При неоднократном сочетании действия предупреждающего сигнала и безусловного пищевого раздражителя формируются...

Уравнение волны. Уравнение плоской гармонической волны. Волновое уравнение. Уравнение сферической волны Уравнением упругой волны называют функцию , которая определяет смещение любой частицы среды с координатами относительно своего положения равновесия в произвольный момент времени t...

Медицинская документация родильного дома Учетные формы родильного дома № 111/у Индивидуальная карта беременной и родильницы № 113/у Обменная карта родильного дома...

Studopedia.info - Студопедия - 2014-2024 год . (0.011 сек.) русская версия | украинская версия