CC BY
16
УДК 519.6
БСТ: 10.14529/ mmp230209
СЕТОЧНО-ХАРАКТЕРИСТИЧЕСКАЯ РАЗНОСТНАЯ СХЕМА ДЛЯ РЕШЕНИЯ УРАВНЕНИЯ ХОПФА НА ОСНОВЕ ДВУХ РАЗЛИЧНЫХ ДИВЕРГЕНТНЫХ ФОРМ
В.Е. Карпов1, А.И. Лобанов1
1 Московский физико-технический институт (национальный исследовательский университет), г. Москва, Российская Федерация
Построено новое двухпараметрическое семейство разностных схем для численного решения уравнения Хопфа. Исходная задача заменялась задачей для системы двух дифференциальных уравнений на основе различных дивергентных форм уравнения Хопфа. Потоковые члены выражались в виде линейных комбинаций переменных, входящих в разные дивергентные формы. В отличие от большинства работ, использующих методы неопределенных коэффициентов для построения разностных схем, при таком подходе неопределенные коэффициенты возникают при формулировке дифференциальной задачи. Система уравнений сохраняет гиперболический тип при любых значениях параметров. Для численной реализации за основу выбрана известная сеточно-характеристическая схема в инвариантах Римана, которая в случае линейного уравнения с постоянными коэффициентами переходит в схему Лакса - Вендроффа. Проведены расчеты двух тестовых задач - об эволюции гладкого начального условия и формировании разрывного решения и о распространении «ударной волны». По результатам тестовых расчетов подобраны коэффициенты экстраполяции, позволяющие получить хорошее согласие с точным решением. Исследовался апостериорный порядок сходимости к предельной функции для разрывных решений. При удачно подобранных коэффициентах экстраполяции он незначительно превышает единицу в момент градиентной катастрофы. При распространении сильного разрыва на больших временах порядок сходимости падает до 0,76. Остается открытым вопрос о постановке оптимизационной задачи, позволяющей выбирать коэффициенты экстраполяции наилучшим образом, возможно, в зависимости от локальных свойств решения. Также открытым пока остается вопрос о создании гибридных разностных схем с переменными коэффициентами экстраполяции в зависимости от гладкости решения.
Ключевые слова: уравнение Хопфа; дивергентная форма; неопределенные коэффициенты; схема Лакса - Вендроффа; разрывное 'решение.
Введение
Уравнения и системы уравнений гиперболического типа встречаются в большом количестве приложений. Характерной чертой решений задач для систем уравнений гиперболического типа является наличие слабых (обобщенных) решений, допускающих как разрывы самих функций, так и разрывы производных. Ввиду большого количества постановок задач для уравнений гиперболического типа не сформировался единый подход к их решению численными методами. Подробный обзор конечно-разностных методов решения задач для уравнений гиперболического типа приведен
Для линейного уравнения переноса С.К. Годуновым установлена теорема, запрещающая существование монотонных конечно-разностных схем с порядком аппроксимации выше первого [2]. Хотя аналогов теоремы Годунова для произвольных уравнений и систем гиперболического типа не установлено, практически те же эффекты
в [1].
проявляются и при реализации численных методов для других систем уравнений гиперболического типа. Под монотонностью здесь понимается отсутствие нефизичных осцилляций численного решения. Немонотонность для схем высокого порядка аппроксимации проявляется, например, при решении уравнения Хопфа методом прямых [3].
Несмотря на большое число численных методов решения систем уравнений гиперболического типа, проблема создания разностных схем, минимизирующих как дисси-пативную, так и дисперсионную составляющие погрешности, по-прежнему актуальна. Известны следующие подходы к ее решению.
Использование дифференциальных следствий самого уравнения. По-видимому, исторически первый способ повышения порядка аппроксимации. В 1962 году на его основе Р.П. Федоренко построен класс гибридных схем, использующих дифференциальные продолжения в области гладкости решения (то есть там, где производные определены) [4]. Дифференциальные продолжения использованы К.М. Магомедовым и А.С. Холодовым при построении класса сеточно-характеристических численных методов [5]. Использование дифференциальных продолжений также привело к созданию класса схем повышенного порядка аппроксимации на минимальном шаблоне («компактные разностные схемы>) [6].
Другой способ повышения порядка аппроксимации конечно-разностных схем предложен Б.В. Роговым. Он заключается в одновременном использовании в качестве расчетных переменных не только самой функции, но и интегральных средних ее значений по отрезку. Схемы, построенные подобным образом, получили название «бикомпактные схемы> [7,8].
В оригинальных работах по бикомпактным схемам такие интегральные средние относились к тем же узлам сетки, что и основные переменные. В принципе можно определить переменные бикомпактной схемы на разнесенных сетках.
Третий подход к повышению пространственного порядка аппроксимации конечно-разностной схемы при попытке уменьшить и диссипативную, и дисперсионную составляющую погрешности предложен сравнительно недавно. В работе В.М. Го-ловизнина и Б.Н. Четверушкина [9] построена новая консервативная балансово-характеристическая схема. При этом для интерполяции «консервативных> переменных используются потоковые величины, определяемые на разнесенной сетке.
Отметим, что в настоящее время все упомянутые подходы к повышению порядка аппроксимации разностных схем активно развиваются. Число публикаций за последние годы по данной тематике велико. Кроме того, создание класса бикомпактных схем по своим идеям близко к идеям разрывного метода Галеркина [10] (где интегральное среднее по отрезку есть один из способов приближения потоковых величин). Идея введения разнесенных сеток для консервативных и потоковых переменных для повышения точности разностных схем активно применяется для построения схем в задачах с несколькими пространственными переменными. Такой подход в принципе можно считать общим как для классов бикомпактных и балансово-характеристических схем, так и для вариантов разрывного метода Галеркина.
Все упомянутые подходы объединяет общая идея - для построения разностной схемы использовать не только само уравнение, но и некоторые его следствия.
В данной работе авторами предлагается новый подход для «расширения> дифференциальной задачи. Как известно, для некоторых уравнений гиперболического типа существует не единственная дивергентная форма представления уравнения. Например, для квазилинейного уравнения Хопфа существует счетное число таких форм [11]. Предлагается для построения разностных схем использовать несколько дивергентных форм, объединяя их в систему уравнений, и строить разностную схему уже для системы уравнений гиперболического типа.
1. Постановка задачи
Рассмотрим следующую дифференциальную задачу для уравнения Хопфа:
ди ди
Наряду с характеристической формой записи уравнения Хопфа (1) будем рассматривать дивергентную форму записи
ди 1 д , 9Л
— + -— (и2) = 0. 2
дЬ 2 дх v ; у '
Эта дивергентная форма соответствует закону сохранения импульса.
Выберем теперь другую дивергентную форму записи уравнения (1). Так как решение (1) отлично от тождественного нуля, то его можно переписать в виде
1 ди ди и дЬ дх
Введем новую переменную V = 1п | и |, тогда последнее равенство будет иметь вид
ди ди п ,оЛ
т. + Ш = а <3)
Это - следствие исходного уравнения (1). Обратное преобразование, переводящее решение (3) в решение (1) - и = sign(u)ev.
Для уравнения (1) потоковое слагаемое представим в виде линейной комбинации потоков, вычисленных с использованием функций и и V:
| + + (4)
где в > 0.
Если теперь рассматривать систему уравнений (3), (4) и решать ее численно, то сеточная функция и будет представлять собой сеточную проекцию решения уравнения Хопфа. Такая форма системы уравнений будет самой простой, но, к сожалению, численные результаты в этом случае получаются неудовлетворительными.
Для получения нового семейства разностных схем введем параметр а > 0 и представим потоковое слагаемое в (3) по аналогии с (4). Получим систему нелинейных уравнений в частных производных
^ + + = о,
dv д
— + — (аи + (1 - a) sign (и) е") = 0. dt дх
Начальные и граничные условия для системы будут сформулированы ниже при описании результатов расчетов тестовых задач.
В систему нелинейных дифференциальных уравнений входят два неопределенных коэффициента. В отличие от большинства работ, использующих методы неопределенных коэффициентов для построения разностных схем, при таком подходе неопределенные коэффициенты возникают при формулировке дифференциальной задачи.
Отметим, что использование методов неопределенных коэффициентов широко использовалось при построении сеточно-характеристических схем [5]. Пространства неопределенных коэффициентов для построения линейных разностных схем впервые были введены в [11]. В [12] на основе анализа в пространствах неопределенных коэффициентов построены гибридные неосциллирующие схемы высокого порядка аппроксимации для нелинейных систем уравнений гиперболического типа.
Известно [13], что (1) допускает запись в виде бесконечного числа дивергентных форм. Форма (1) и дивергентная форма, учитывающая закон сохранения энергии для уравнения Хопфа, использовалась для анализа семейства разностных схем в пространстве неопределенных коэффициентов для решения уравнения Бюргерса с преобладающей конвекцией [14].
В ходе эволюции даже гладкого начального возмущения при решении уравнения Хопфа могут сформироваться разрывы [15]. Классическое решение уравнения перестает существовать, существует только слабое (обобщенное) решение. Актуальным является проблема построения разностных схем, позволяющих правильно воспроизводить явление градиентной катастрофы, правильно воспроизводить скорость фронта и давать минимальные осцилляции при возникновении разрывов решения.
2. Сеточно-характеристическая схема второго порядка аппроксимации по пространству
В качестве рабочей схемы строится вариант сеточно-характеристической схемы [5], имеющий второй порядок аппроксимации по пространственной переменной. Несмотря на то, что схема фиксирована, она погружена в двумерное пространство неопределенных коэффициентов. В случае, когда 0 < а, в < 1, выполняется интерполяция при вычислении пространственных производных для обеих переменных, в противном случае - экстраполяция хотя бы для одной переменной.
Для построения разностной схемы перепишем систему в характеристической форме
д_(и\ (Ци (1 — /3) е2"" \д_(и\_ г ,
дг\у) + \а (1-а) ^п (иу) дх ^)
Собственные числа матрицы системы (5) будут
1
А- = 2
А
2
[вu + (1 - a) sign(u)ev] +
\Ри + (1 - a) sign(u)ev] - л/D
2
здесь D — дискриминант возникающего квадратного уравнения
D = [вu + (1 - a) sign(u)ev]2 - 4 [в (1 - а) Ы ev - а (1 - в) e2v] .
Для точного решения уравнения Хопфа должно выполняться равенство e2v = u2, тогда, после упрощения, получим следующую оценку дискриминанта:
D = (в - 1 - a)2 u2. (6)
Из (6) следует, что собственные числа всегда действительные, тогда система (5) имеет гиперболический тип при любых значениях параметров а и в •
1
Соответствующие левые собственные векторы для собственных чисел матрицы системы (5) выписаны ниже:
1 = (а, А1 — ви),
2 = (а, А2 — ви) •
При переходе в базис из левых собственных векторов система уравнений расщепляется. При этом возникают новые расчетные переменные - инварианты Римана. Для системы (5) они имеют вид:
11 = аи + (А1 — в и) V,
12 = аи + (А2 — в и) V.
Обратный переход от инвариантов Римана к естественным переменным задачи и и V осуществляется пересчетом
и\_ —1 /А2 — в и в и — А Л (11
VI а,л/Г> \ а ) \12
Введем некоторые ограничения на выбор параметров а и в. Для собственных чисел матрицы (5) с учетом (6) получаем оценки
= \ [(/3 + (1 — а)) и + \/3 — 1 — а\ |м|] ,
А2 = ^ [(/3 + (1 — а)) и — — 1 — а\ \и\].
Потребуем сначала, чтобы при положительном значении и оба собственных числа были неотрицательными. Как легко показать, для этого требуется выполнение условия в > а• Потребуем далее, чтобы при отрицательном значении и первое собственное число было неотрицательным, а второе - неположительным. Для выполнения этого требования должно выполняться неравенство в ^ а•
Построим теперь двухпараметрическое семейство сеточно-характеристических схем. В расчетной области введем равномерную расчетную сетку по времени £ и по
пространству х: _
Г = пт, хк = х0 + кН, к = О, N.
Значения и и V при п = 0 задаются начальными условиями. Значения и и V при к = 0(или при к = N, в зависимости от постановки задачи) задаются граничными условиями. Для расчетов будем использовать разностную схему второго порядка аппроксимации по пространственным переменным. В качестве основы для построения схемы выбран вариант сеточно-характеристической схемы [5]. В случае линейного уравнения переноса с постоянными коэффициентами она переходит в известную схему Лакса - Вендроффа. Разностные уравнения для инвариантов Римана выглядят следующим образом. Если ит > 0, тогда
с1 = (1 - ы:) (11)1+мпт (ьг^ - \ (ыпт)2 тпт+1 - 2 (1гТт+(ьТт-х)
= (1 - ) №)ш + (^2Тт (ЬТт_1 - \ ((^2)™ )2 (№)"+1 - 2 (/2)" + (/2)"
здесь введены обозначения (01 и (02)^ для локальных значений чисел Куранта для
си (-2):
инвариантов
(rr \п Т („уть (al)m = -Г-, ( 2jr
h h
Если < 0, тогда
(h)T = (1 - № (/1): + (<7хС - \ ((aiTJ2 ((h)nm+l - 2 (hZ + (hZ.,) ,
= (1 + (02)™) (h)nm - (02)™ - ^ ((^2)™ )2 - 2 (/2)™ + (/2)™ _i) •
Значения и и V естественных переменных на (п + 1)-м временном слое восстанавливаются
«у+1 = 1 ((А2 - - (Л1 - //Лга+1
,ут К«А2): - мо:] I <
Для устойчивости разностной схемы потребуем выполнения условия Куранта
н
Т<
max(|Ai| , |Л2|):'
3. Результаты тестовых расчетов
С помощью численной реализации описанной разностной схемы проведены численные расчеты нескольких тестовых задач. Результаты показывают, что лучшие результаты разностная схема дает при использовании экстраполяции в дифференциальной задаче а > 1,в > 1- При а = 1, в =1 уравнения в системе расщепляются, для решения первого уравнения системы предложенная разностная схема переходит в традиционный вариант схемы Лакса - Вендроффа.
Тестовая задача 1. В качестве начальных условий выбираем гладкую функцию
х
и(0,х) = А-сЬ"2 -.
Известно, что в ходе эволюции данного начального возмущения теряется гладкость [15], в решении образуется сильный разрыв.
Начальные условия для функции V очевидны, и(0,х) = 1пА — 21п(сЬ Для определенности положим А = 1. В качестве граничных условий выбираем условия периодичности на отрезке хЕ [-5, 5]
и (¿, -5) = и (¿, 5).
На рис. 1 приведены график точного решения задачи, построенного методом характеристик (пунктирная линия), и графики численного решения на сетке с 3200 ячейками по пространственной переменной и значении числа Куранта 0,4 при значениях параметров а =1, 01 и в =1,1 (сплошная линия) и а = 1,1 и в = 1,15 (штриховая линия) для времени Т = 2, 598 (момент градиентной катастрофы). Пройдено 2079 шагов по времени. Видно, что совпадение графиков удовлетворительное, но численное решение несколько опережает точное из-за того, что численная погрешность носит преимущественно дисперсионный характер.
Т = 2.598
х
Рис. 1. Сравнение численного решения при значениях параметров 1,01 и 1,1 (сплошная линия) и 1,1 и 1,15 (штриховая линия) с точным решением дифференциальной задачи, построенным методом характеристик, (пунктирная линия) для времени Т = 2,598 (момент градиентной катастрофы)
Для количественной оценки порядка сходимости к предельной сеточной функции при утрате гладкости решения вычислялась норма разности численных решений на различных сетках (аналог нормы Ь2) при фиксированном времени
егг 2
\
м
^У] « - «3200^,Ьп))2. (7)
i=0
Расчет проводился на последовательности сгущающихся сеток с числом ячеек 100, 200, 400, 800, 1600 при максимальном значении числа Куранта 0,4. Здесь число Куранта вычислялось по начальным данным, к = ^-тах ((А^ , |А2|)°- Зависимость нормы
погрешности от шага по пространственной координате приведена на рис. 2 в двойных логарифмических координатах. Апостериорный порядок сходимости определялся тангенсом угла наклона прямой, построенной по расчетным данным методом наименьших квадратов. Для момента градиентной катастрофы он оказался равным 1,02 для а =1, 01 и в = 1,1 (сплошная линия) и 1,05 для а = 1,1 и в = 1,15 (штриховая линия).
На рис. 3 показана зависимость нормы погрешности в зависимости от шага по времени при фиксированном шаге по пространственной переменной (от числа Куранта). Расчеты проводились для значений числа Куранта (су) 0,4, 0,2, 0,1, 0,05, график приведен в двойных логарифмических координатах. Для рассматриваемой разностной схемы погрешность, связанная с изменением шага по времени, существенно меньше погрешности, определяемой пространственной дискретизацией, поэтому схема практически нечувствительна к изменению шага по времени (числа Куранта).
Аналогичные оценки проведены для времени 5,0 (после 4000 шагов по времени) На рис. 4 приведены график точного решения задачи, построенного методом характеристик (пунктирная линия) и графики численного решения на сетке с 3200 ячейками по пространственной переменной и значении числа Куранта 0,4 при значениях параметров а = 1, 01 и в =1,1 (сплошная линия) и а = 1,1 и в = 1,15 (штриховая линия). Видно, что совпадение графиков удовлетворительное, но теперь наметилось
T = 2.598
1.01, 1.10 +
1.10, 1.15 □
0.01 0.1 h
Рис. 2. Оценка порядка сходимости схемы для момента градиентной катастрофы
T = 2.598
1.01, 1.10 -1-
1.10, 1.15 —н —
Рис. 3. Оценка зависимости погрешности разностной схемы от шага по времени (числа Куранта) для момента градиентной катастрофы
cv
некоторое отставание фронта численного решения от точного. Это отставание можно объяснить тем, что для численного решения нарушаются имеющиеся в задаче законы сохранения (построенная схема неконсервативна). Так как в основу построения разностной схемы положена немонотонная схема Лакса - Вендроффа для каждого инварианта, то можно ожидать, что построенная схема немонотонна. Рис. 4 иллюстрирует появление нефизичных осцилляций в окрестности разрыва.
т = 5.0
X
Рис. 4. Сравнение численного решения при значениях параметров 1,01 и 1,1 (сплошная линия) и 1,1 и 1,15 (штриховая линия) для времени Т = 5,0 (пройдено 4000 шагов по времени) с точным, построенным методом характеристик (пунктирная линия)
Для количественной оценки порядка сходимости при времени 5,0 также использована норма (7). Расчет проводился на последовательности сгущающихся сеток с числом ячеек 100, 200, 400, 800, 1600 при фиксированном числе Куранта 0,4. Зависимость нормы погрешности от шага по пространственной координате приведена на
Т = 5.0
Т = 5.0
1.01, 1.10 +
1.10, 1.15 □
1.01, 1.10 1.10, 1.15
0.1 СУ
Рис. 5. Оценка порядка сходимости схемы для момента Т = 5,0
Рис. 6. Оценка погрешности схемы для момента Т = 5,0 при изменении шага по времени (числа Куранта)
рис. 5 в двойных логарифмических координатах. Апостериорный порядок сходимости определялся тангенсом угла наклона прямой, построенной по расчетным данным методом наименьших квадратов. Порядок сходимости разностной схемы оказался равным 0,76 для а = 1, 01 и в =1,1 (сплошная линия) и 0,78 для а = 1,1 и в =1,15 (штриховая линия).
На рис. 6 показана зависимость нормы погрешности в зависимости от шага по времени (от числа Куранта) при фиксированном шаге по пространственной переменной для времени Т = 5, 0. Расчеты проводились для чисел Куранта (су) 0,4, 0,2, 0,1, 0,05. график приведен в двойных логарифмических координатах. Для рассматриваемой разностной схемы погрешность, связанная с изменением шага по времени, существенно меньше погрешности, определяемой пространственной дискретизацией, поэтому схема и для этого момента практически нечувствительна к изменению шага по времени (числа Куранта).
Тестовая задача 2. Разные авторы проводят тестирование своих разностных схем на отличных друг от друга постановках задач. В работе [16] приведены результаты тестирования схемы (разрывный метод Галеркина) на точном разрывном решении уравнения Хопфа:
и (х, £)
х + 1
¿ + 2 х
.ГГТ'
, x<xf (£) х > Xf (£)
, где (¿) = л/а + 2) а + 1) + I + 1. (8)
Такого рода решение моделирует распространение ударных волн в газе. Используем это решение в качестве тестового. В качестве начальных данных бралось значение точного решения при £ = 0, в качестве граничного значения на каждом временном слое при х = 0 бралось значение точного решения в этой точке при соответствующем моменте времени. Предложенная разностная схема использована для воспроизведения данного разрывного решения на сетке, совпадающей с использованной в [16] (400 и 800 ячеек по пространству). Расчеты велись с числом Куранта 0,4. Наилучшее совпадение с точным решением достигнуто при значениях коэффициентов
а = в = 1, 0365.
и
T = 5.0
х
Рис. 7. Сравнение численного решения для Ь = 0,5 и сеток с числом ячеек 400 (штриховая линия) и 800 (сплошная линия) по пространственной переменной с точным решением (пунктирная линия)
т = 5.0
0.1 г
400 cells------|
0.09 - 800 cells................|
0.08 - I
I
0 07 " |
0.06 - I
r i L 0.05 - 1
<U I
0.04 - f
0.03 - ;j
0.02 - (!
0.01 - .j
r:__
0 -----1-----1-----^-----I'll_I
0 0.2 0.4 0.6 0.8 1
x
Рис. 8. Погрешность приближенного решения при времени 0,5 для сеток с 400 ячейками (штриховая линия, 500 шагов по времени) и 800 ячейками (пунктирная линия, 1000 шагов по времени)
Сравнение численного и точного решения для времени 0,5 приведено на рис. 7. Видно, что на данном тесте качество приближения заведомо не хуже, чем для разрывного метода Галеркина [16].
На рис. 8 приведена погрешность численного решения для времени 0,5 (500 и 1000 шагов по времени для сеток 400 и 800 ячеек по пространству соответственно). В области гладкости функции погрешность незначительна. Она возрастает в окрестности разрыва решения. Несмотря на то, что при сгущении сетки имеется стремление к точному решению, амплитуда графика погрешности в окрестности разрыва остается примерно постоянной.
Заключение
Для повышения априорного порядка аппроксимации разностных схем и уменьшения нефизичных осцилляций на разрывах при построении разностных схем обычно применяются некоторые дополнительные условия, представляющие собой следствия исходного дифференциального уравнения. Это могут быть, например, дифференциальные продолжения или интегральные средние на сетке. При использовании таких подходов можно применить метод анализа разностных схем в пространстве неопределенных коэффициентов для построения разностной схемы.
В предлагаемом подходе введение в рассмотрение неопределенных коэффициентов предполагается на этапе формулировки дифференциальной задачи. Мы рассматриваем проблему построения разностной схемы для квазилинейного уравнения гиперболического типа (уравнения Хопфа) как проблему построения схемы для системы уравнений, причем для потоковых переменных выполняется экстраполяция на основе использования двух дивергентных форм. Тестовые расчеты показали, что при удачном выборе коэффициентов экстраполяции разностная схема показывает неплохие результаты на разрывных решениях уравнения Хопфа.
При проведении численных расчетов неопределенные коэффициенты подбирались вручную. Открытым остается вопрос о постановке оптимизационной задачи,
позволяющей выбирать коэффициенты экстраполяции наилучшим образом, возможно, в зависимости от локальных свойств решения. Также открытым пока остается вопрос о создании гибридных разностных схем (с переменными коэффициентами экстраполяции и разными аппроксимационными свойствами в зависимости от гладкости решения).
Литература
1. Куликовский, А.Г. Математические вопросы численного решения гиперболических систем уравнений / А.Г. Куликовский, Н.В. Погорелов, А.Ю. Семенов. - М.: Физматлит, 2012.
2. Годунов, С.К. Разностный метод численного расчета разрывных решений уравнений гидродинамики / С.К. Годунов // Математический сборник. - 1959. - Т. 47 (89), № 3. -С. 271-306.
3. Гужев, Д.С. Уравнение Бюргерса - тест для численных методов / Д.С. Гужев, Н.Н. Ка-литкин // Математическое моделирование. - 1995. - Т. 7, № 4. - С. 99-127.
4. Федоренко, Р.П. Применение разностных схем высокой точности для численного решения гиперболических уравнений /Р.П. Федоренко // Журнал вычислительной математики и математической физики. - 1962. - Т. 2, № 6. - С. 1122-1128.
5. Магомедов, К.М. Сеточно-характеристические численные методы / К.М. Магомедов, А.С. Холодов. - М.: Юрайт, 2017.
6. Толстых, А.И. Компактные и мультиоператорные аппроксимации высокой точности для уравнений в частных производных / А.И. Толстых. - М.: Наука, 2015.
7. Рогов, Б.В. Монотонные бикомпактные схемы для линейного уравнения переноса / Б.В. Рогов, М.Н. Михайловская // Математическое моделирование. - 2011. - Т. 23, № 6. - С. 98-110.
8. Аристова, Е.Н. Бикомпактные схемы для неоднородного линейного уравнения переноса / Е.Н. Аристова, Д.Ф. Байдин, Б.В. Рогов // Математическое моделирование. - 2013. -Т. 25, № 5. - С. 55-66.
9. Головизнин, В.М. Алгоритмы нового поколения в вычислительной гидродинамике / В.М. Головизнин, Б.Н. Четверушкин // Журнал вычислительной математики и математической физики. - 2018. - Т. 58, № 8. - С. 20-29.
10. Cockburn, B. An Introduction to the Discontinuous Galerkin Method for Convection-Dominated Problems / B. Cockburn // Lecture Notes in Mathematics. - Berlin: SpringerVerlag, 1998. - P. 151-268.
11. Холодов, А.С. О построении разностных схем с положительной аппроксимацией для уравнений гиперболического типа / А.С. Холодов // Журнал вычислительной математики и математической физики. - 1978. - Т. 18, № 6. - С. 1476-1492.
12. Холодов, А.С. О критериях монотонности разностных схем для уравнений гиперболического типа / А.С. Холодов, Я.А. Холодов // Журнал вычислительной математики и математической физики. - 2006. - Т. 46, № 9. - С. 1638-1667.
13. Федоренко, Р.П Введение в вычислительную физику / Р.П. Федоренко. - Долгопрудный: Интеллект, 2009.
14. Basharov, I.V. On the Finite Difference Schemes for Burgers Equation Solution / I.V. Basharov, A.I. Lobanov // Smart Modelling for Engineering Systems. - Singapore, 2021. - V. 215, P. 151-167.
15. Уизем, Дж.Б. Линейные и нелинейные волны / Дж.Б. Уизем. - М.: Мир, 1977.
16. Криксин, Ю.А. Построение точных решений некоторых уравнений гиперболического типа, содержащих разрыв, распространяющийся по неоднородному фону /Ю.А. Криксин, П.А. Кучугов, М.Е. Ладонкина, О.А. Неклюдова, В.Ф. Тишкин // Препринты Института прикладной математики им. М.В. Келдыша. - М., 2018. - № 17. - 14 c.
Владимир Ефимович Карпов, кандидат физико-математических наук, доцент, кафедра информатики и вычислительной математики, Московский физико-технический институт (национальный исследовательский университет) (г. Москва, Российская Федерация), [email protected].
Алексей Иванович Лобанов, доктор физико-математических наук, профессор, кафедра вычислительной физики, Московский физико-технический институт (национальный исследовательский университет) (г. Москва, Российская Федерация), а1ехеу. i. lobanov@gmail. сот.
Поступила в редакцию 19 февраля 2023 г.
MSC 35L60, 65M06 DOI: 10.14529/mmp230209
GRID-CHARACTERISTIC DIFFERENCE SCHEME FOR SOLVING THE HOPF EQUATION BASED ON TWO DIFFERENT DIVERGENT FORMS
V.E. Karpov1, A.I. Lobanov1
:Moscow Institute of Physics and Technology (National Research University), Moscow,
Russian Federation
E-mails: [email protected], [email protected]
A new two-parameter family of difference schemes for the numerical solution of the Hopf equation is constructed. The original problem was replaced by a problem for a system of two differential equations based on various divergent forms of the Hopf equation. The flux terms were expressed as linear combinations of the variables included in the different divergent forms. In contrast to most works that use uncertain coefficient methods to construct difference schemes, in this approach uncertain coefficients arise in the formulation of the differential problem. The system of equations retains the hyperbolic type for any parameter values. For the numerical implementation, the well-known grid-characteristic scheme in Riemann invariants is chosen as the basis, which in the case of a linear equation with constant coefficients passes into a Lax-Wendroff scheme. Calculations of two test problems -on the evolution of a smooth initial condition and the formation of a discontinuous solution and on the propagation of a "shock wave" - have been performed. Based on the results of the test calculations, we selected extrapolation coefficients that allow us to obtain a good agreement with the exact solution. The a posteriori order of convergence to the limit function for discontinuous solutions was investigated. When the extrapolation coefficients are well chosen, it insignificantly exceeds 1 at the moment of a gradient catastrophe. The order of convergence decreases to 0,76 when the strong discontinuity propagates at large times. The question of setting an optimization problem that allows one to choose the extrapolation coefficients in the best way, possibly depending on the local properties of the solution, remains open. The question of creating hybrid difference schemes with variable extrapolation coefficients depending on the smoothness of the solution also remains open.
Keywords: Hopf equation; divergent form; undefined coefficients; Lax-Wendroff scheme; discontinuous solution.
References
1. Kulikovskii A.G., Pogorelov N.V., Semenov A.Yu. Matematicheskie voprosy chislennogo resheniya giperbolicheskih system uravneniy [Mathematical Aspects of Numerical Solution of Hyperbolic Systems]. Moscow, Fizmatlit, 2012. (in Russian)
2. Godunov S.K. [A Difference Scheme for Numerical Solution of Discontinuous Solution of Hydrodynamic Equations]. Matematicheskii Sbornik, 1959, vol. 89, no. 3, pp. 271-306. (in Russian)
3. Guzhev D.S., Kalitkin N.N. [Burgers Equation is a Test for Numerical Methods]. Mathematical Models and Computer Simulation. 1995, vol. 7, no. 4, pp. 99-127. (in Russian)
4. Fedorenko R.P. The Application of Difference Schemes of High Accuracy to the Numerical Solution of Hyperbolic Equations. USSR Computational Mathematics and Mathematical Physics, 1963, vol. 2, no. 6, pp. 1355-1365. DOI: 10.1016/0041-5553(63)90351-3
5. Magomedov K.M., Kholodov A.S. Setochno-harakteristicheskie chislennye metody [Grid-and-Characteristics Numerical Methods]. Moscow, Yurait, 2017. (in Russian)
6. Tolstykh A.I. Kompaktnye i multioperatornye approximatsii vysokoy tochnosti dlya uravneniy v chastnyh proizvodnyh [Compact and High-Precision Multi-Operator Approximations for Partial Differential Equations]. Moscow, Nauka, 2015. (in Russian)
7. Rogov B.V., Mikhailovskaya M.N. The Monotonic Bicompact Schemes for a Linear Transfer Equation. Mathematical Models and Computer Simulations, 2012, vol. 4, no. 1, pp. 92-100. DOI: 10.1134/S2070048212010103
8. Aristova E.N., Baydin D.F., Rogov B.V. Bicompact Scheme for Linear Inhomogeneous Transport Equation. Mathematical Models and Computer Simulations, 2013, vol. 5, no. 6, pp. 586-594. DOI: 10.1134/S2070048213060033
9. Goloviznin V.M., Chetverushkin B.N. New Generation Algorithms for Computational Fluid Dynamics. Computational Mathematics and Mathematical Physics, 2018, vol. 58, no. 8, pp. 1217-1225. DOI: 10.1134/S0965542518080079
10. Cockburn B. An Introduction to the Discontinuous Galerkin Method for Convection-Dominated Problems. Lecture Notes in Mathematics. Berlin, Springer-Verlag, 1998, pp. 151-268. DOI: 10.1007/BFb0096353
11. Kholodov A.S. Construction of Difference Schemes with Positive Approximation for Hyperbolic Equations. USSR Computational Mathematics and Mathematical Physics, 1978, vol. 18, no. 6, pp. 116-132. DOI: 10.1016/0041-5553(78)90141-6
12. Kholodov A.S., Kholodov Ya.A. Monotonicity Criteria for Difference Schemes Designed for Hyperbolic Equations. Computational Mathematics and Mathematical Physics, 2006, vol. 46, no. 9, pp. 1560-1588. DOI: 10.1134/S0965542506090089
13. Fedorenko R.P. Vvedenie v vychislitel'nuyu fiziku [Introduction to Computational Physics]. Dolgoprudny, Intellekt, 2009. (in Russian)
14. Basharov I.V., Lobanov A.I. On the Finite Difference Schemes for Burgers Equation Solution. Smart Modelling for Engineering Systems, Singapore, 2021, vol. 215, pp. 151-167. DOI: 10.1007/978-981-33-4619-2_12
15. Whitham G.B. Linear and Nonlinear Waves. New York, John Wiley and Sons, 1974.
16. Kriksin Yu.A., Kuchugov P.A., Ladonrina M.E., Neklyudova O.A., Tishkin V.F. [Construction of Exact Solutions of Some Equations of Hyperbolic Type Containing a Discontinuity Propagating over an Inhomogeneous Background]. Keldysh Institute Preprints, Moscow, 2018, no. 17, 14 p. DOI: 10.20948/prepr-2018-17. (in Russian)
Received February 19, 2023
