3. Схема «крест».

Это наиболее простая и довольно точная однородная разностная схема счета газодинамики. Ее шаблзн приведен на рис. 98; значения радиусов приписываются узлам сетки, значения скорости — границам пространственных интервалов на полуцелых слоях, а значения плотности, давления и внутренней энергии — серединам интервалов на целых слоях.

Построение схемы напоминает акустический «крест». Для простоты записи выберем равномерные по массе и времени шаги и t и аппроксимируем систему следующими разностными уравнениями:

Рис. 98.

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

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

Таким образом, вязкое давление (65) принимает вид

где — скорость звука. Выражение (67) написано для плоского случая; но обычно им пользуются при любой симметрии задачи.

Аппроксимация. Из вида шаблона на рис. 98 и симметричного написания схемы (66) нетрудно заметить, что на течениях без сжатий, когда псевдовязкость (67) обращается в нуль, схема «крест» имеет локальную аппроксимацию

На течениях со сжатиями (в том числе — с ударными волнами) псевдовязкость отлична от нуля. Правда, квадратичный член в (67а) имеет величину но линейный член имеет величину и, тем самым, ухудшает порядок аппроксимации. Кроме того, вязкие члены записываются не вполне симметрично по времени. В итоге аппроксимация ухудшается до

Нахождение разностного решения. Схема (66) — явная; вычисления по ней проводятся следующим образом. Пусть все величины на исходном слое известны. Тогда из разностного уравнения импульса (66а) находим во всех интервалах; затем из второго уравнения (66б) определяем а из уравнения (66в) - .

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

Замечание 1. Уравнение энергии в (66) можно сделать яным, используя в нем только значение с исходного слоя:

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

Устойчивость схемы можно исследовать методом разделения переменных, линеаризируя схему и замораживая коэффициенты. Громоздкие выкладки приводят к условию устойчивости типа Куранта.

Например, на гладких течениях с нулевой вязкостью схема устойчива при

Для идеального газа и условие (69) принимает вид где есть адиабатическая скорость звука. На течениях с ненулевой вязкостью ограничение на шаг несколько более сильное; при квадратичной вязкости условие устойчивости принимает вид

где — скачок скорости на ударной волне. Хотя это исследование не является строгим, тем не менее данное условие устойчивости хорошо подтверждается на практике.

Рис. 99.

Таким образом, «крест» — условно устойчивая схема. Отметим любопытное обстоятельство. Для расчета гладких течений вязкость не нужна. А если рассчитать без вязкости ударную волну (выбирая небольшое удовлетворяющее условию (70)), то получим «разболтку», изображенную на рис. 99. Этот расчет устойчив, поскольку амплитуда колебаний не возрастает со временем. Но сходимости к физически правильному решению при нет, так как на разрыве потеряна аппроксимация.

Сходимость газодинамической схемы «крест» не доказана. Однако эта схема успешно используется в расчетах примерно с 1950 г. и проверена на многих трудных задачах с известными точными решениями. При стремлении шагов к нулю наблюдалась сходимость к правильному решению, если шаги удовлетворяли условию устойчивости.

Замечание 2. Схема (66) неконсервативна; однако ее дисбаланс стремится к нулю при

Замечание 3. Газодинамические задачи с очень тонкими слоями особенно трудны для расчета. В самом деле, если , то для вычисления с удовлетворительной точностью по формуле (66в) надо знать радиусы с очень высокой точностью, сравнимой с ошибками округления на ЭВМ. В подобных задачах иногда приходится вести расчет с двойным числом знаков или специально видоизменять разностную схему.