CC BY
94
УДК 519.6
В.Л. БУЧАРСКИЙ, Е.М. КАЛИНЧУК
СИММЕТРИЧНЫЕ РАЗНОСТНЫЕ СХЕМЫ МЕТОДА СОВМЕСТНОЙ АППРОКСИМАЦИИ ДЛЯ РЕШЕНИЯ ЛИНЕЙНОГО УРАВНЕНИЯ ПЕРЕНОСА
Анотація. У роботі наведені результати застосування методу сумісної апроксимації для побудови симетричних різницевих схем підвищеного порядку точності на компактному шаблоні. Побудовані різницеві схеми від другого до восьмого порядку апроксимації за часом і простором для одно-вимірного лінійного рівняння переносу. Показані результати тестових розрахунків, які підтверджують теоретичні викладки.
Ключові слова: метод сумісної апроксимації, симетрична схема.
Аннотация. В работе приведены результаты применения метода совместной аппроксимации для построения симметричных разностных схем повышенного порядка точности на компактном шаблоне. Построены разностные схемы от второго до восьмого порядка аппроксимации по времени и пространству для одномерного линейного уравнения переноса. Представлены результаты тестовых расчетов, подтверждающие теоретические выкладки.
Ключевые слова: метод совместной аппроксимации, симметричная схема.
Abstract. The results of using the joint approximation method for designing the high order symmetrical finite difference schemes on a compact stencil are presented. Finite difference schemes from second to eighth order of approximation on time and space for the one-dimensional linear equation of transfer have been built. The results of the test calculations that prove the theoretical outfits are shown.
Keywords: joint approximation method, symmetric scheme.
1. Введение
На сегодняшний день численные методы широко используются в различных областях науки и техники для исследования физических процессов. При этом создается вычислительная модель приближенного решения. Основой таких моделей есть уравнения математической физики, точное решение которых является сложной задачей вследствие многообразия возможных граничных условий. Степень соответствия приближенного решения точному определяется погрешностью аппроксимации, поэтому естественна предпочтительность методов с меньшей погрешностью, то есть, методов более высокого порядка точности.
Ранее автором был предложен новый метод построения компактных разностных схем повышенного порядка точности - метод совместной аппроксимации [1-3]. Приведем кратко идею построения разностных схем методом совместной аппроксимации следуя [3]. В основе этого метода лежит принцип аппроксимации разностной схемы исходного уравнения на пространстве функций, являющихся решением последнего. Следующий из этого принцип эквивалентности частных производных позволяет увеличить число свободных параметров при создании схемы без изменения шаблона. При этом учитывается то, что разностная схема аппроксимирует точно дифференциальное уравнение бесконечного порядка, так называемую гиперболическую форму дифференциального представления разностной схемы (далее - Г-форму дифференциального представления разностной схемы [4]). Все это позволяет повысить порядок локальной аппроксимации разностной схемы.
Разностные аппроксимации производных, входящих в исходное уравнение, строятся методом неопределенных коэффициентов. Однако для определения этих коэффициентов, помимо условий аппроксимации, используется не Г-форма, а параболическая форма дифференциального представления разностной схемы (далее - П-форма). Она получается из Г-формы заменой производных по времени и смешанных производных производными по координате с использованием дифференциальных следствий Г-формы. Как следствие,
© Бучарский В.Л., Калинчук Е.М., 2011
ISSN 1028-9763. Математичні машини і системи, 2011, № 4
невязка разностной схемы содержит меньше условии, подлежащих удовлетворению для повышения порядка аппроксимации. Это позволяет либо повысить порядок аппроксимации на заданном шаблоне по сравнению с традиционными подходами, либо при заданном порядке точности сократить число точек в шаблоне разностной схемы.
Использование предложенного выше подхода к построению разностных схем приводит к аппроксимации не отдельных составляющих исходного уравнения, а уравнения в целом (входящие в него производные аппроксимируются совместно).
Ранее в работах [1, 3] для построения разностных схем использовались несимметричные противопоточные шаблоны. Это не давало возможности проводить расчеты при отрицательных значениях скорости переноса.
2. Применение симметричных разностных схем для решения уравнения переноса методом совместной аппроксимации
Цель данной работы состоит в построении симметричных разностных схем для решения дифференциальных уравнений в частных производных методом совместной аппроксимации и применении предлагаемого подхода к решению простейшего гиперболического уравнения - уравнения переноса при произвольных значениях скорости переноса.
Построим методом совместной аппроксимации разностную схему для простейшего гиперболического уравнения - линейного уравнения переноса:
— + а — = 0, а = const > 0 (1)
dt дх w
на десятиточечном шаблоне, включающем девять точек по пространству и два слоя по времени.
Для выполнения свойства консервативности исходную схему запишем в потоковой форме:
—Г - fn ,а—1 - /”), 1
С—+1 ^п
- + ’>1Л’ + Т4 (а/3 + а2/гП-2 + Оз/м + а4/г" + а5/-1 + = 0 (2)
х 2к к
Здесь А^- = g- - g--1 - оператор конечной разности назад по пространст-
ву, а}., ] = 1,8 - свободные параметры, подлежащие определению.
Раскладывая /++1, /+ в ряд Тейлора в окрестности точки (1п, хг) и подставляя эти
разложения в разностную схему (2), получим ее дифференциальное представление, главный член которого
Л/ Л/
----+ (а + (а1 + а2 + а3 + а4 + а5 + а6 + а7 + а8))— = 0(х, к).
Л Лх
Как видно, для аппроксимации исходного уравнения (1) необходимо, чтобы а1 + а2 + а3 + а4 + а5 + а6 + а7 + а8 = 0 .
Полагая, что
а4 = -а1 -а2 -а3 -а5 -а6 -а7 -а8, (3)
получим Г-форму разностной схемы (2). Ввиду ее громоздкости, ограничимся членами первого порядка по х, к :
— — о о \д2— Л д2/
----+ а-+ h(-3ai — 2а2 — аз + а5 + 2а6 + 3а7 + 4а8) —— +1(----------------—)
dt дх дх 2 dt
22
Далее, заменим производные по времени и смешанные производные через производные по координате, используя дифференциальные следствия Г-формы. Поскольку после подстановки (3) в схеме (2) осталось семь свободных параметров, ограничимся заменой производных до девятого порядка включительно. Результатом будет П-форма девятого дифференциального приближения разностной схемы (2):
Э/ Э/ -Л ,_1 Э 3/
— + а — + > к1 1С,—— = 0.
Эх Эх 1=2 1 Эх1
Выражения для С. достаточно громоздки (например,
а 1
С4 = + ^(-31а1 - 12а2 - 3а3 + а5 + 4а6 + 13а7 + 32а8)о +
а(-3а1 - 2а2 - а3 + а5 + 2а6 + 3а7 + 4а8 )2 о +
а 5 о 4
5
- +
1 3 2 1
—а о (1 - 18а1о + 6(-2а2 -а3 + а5 + 2а6 + 3а7 + 4а8)о + — (26а1 + 7а2 + а3 + 7а7 + 26а8 -6 4
4(а1 + 2а2 + а3 - а5 - 2а6 - 3а7 - 4а8 )(6а1 + 3а2 + а3 + а6 + 3а7 + 6а8 )о),
поэтому в дальнейшем приводить их не будем.
Для повышения порядка малости невязки до восьмого необходимо последовательно
Э ^ /
приравнять нулю коэффициенты С. при ——, к = 2,8. Используя полученную систему ал-
1 Эхк
гебраических уравнений и приравнивая некоторые из ак, к = 2,8 нулю, несложно построить разностные схемы от второго до восьмого порядка точности. Структура этих систем уравнений приведена в табл. 1.
Таблица 1. Структура систем уравнений
Свободные параметры в схемах нечетного порядка точности использовались для симметризации разностных схем.
Устойчивость полученных схем была исследована методом Неймана и методом дифференциального приближения. Все схемы оказались устойчивыми при
|а| т
1.
к
Как видно, полученные разностные схемы позволяют проводить устойчивые расчеты при любых знаках скорости переноса а .
3. Результаты тестовых расчетов
Результаты, полученные выше аналитически, были подтверждены серией тестовых расчетов. Для проверки порядка аппроксимации полученных разностных схем было решено уравнение переноса (1) при а = 1 и а = -1
Порядок точно- сти Система уравнений для определения параметров
2 С2 = 0; а1 = а2 = а3 = а6 = а7 = а8 = 0
3 С2 = С3 = 0; а1 = а2 = а7 = а8 = 0
4 С2 = С3 = С4 = 0; а1 = а2 = а7 = а8 = 0
5 С2 = С3 = С4 = С5 = 0; а1 = а8 = 0
6 С = С = С = С = С = 02 ^3 4 5 ^6 а1 = а8 = 0
7 С2 = С3 = С4 = С5 = С6 = С7 = 0
8 С2 = С3 = С4 = С5 = С6 = С7 = С8 = 0
Э/ + а Э/ = 0,
Эх Эх (4)
/(0,х) = /0(x),хе [-1,1)
с периодическими граничными условиями. Решение было получено для двух начальных условий.
В первой задаче
./аСх) = SІn(px), (5)
во второй
/0( х) = 8т4(рх). (6)
Расчет проводился до X = 2 - до завершения одного цикла. Эти задачи считаются традиционными для исследования сходимости разностных схем.
Во всех задачах полагалось о = 0,5 . Известное точное решение задачи (4)
/ (х, х) = /0( х - ах)
позволило определить погрешности решения е = / -/" в каждой точке по окончании расчета и экспериментально оценить скорость сходимости и, соответственно, порядок аппроксимации. Оценка погрешности решения е проводилась в конечномерных аналогах норм в Ь1, Ь2 С :
N
тА\.к = к > \и,I, тА\ к =
II ¿2
N
2
к>
\ил , и = тахи
11С
г=1^ 1
Здесь и - некоторая сеточная функция, N - число узлов сетки по пространству. Порядок аппроксимации определялся по экстраполяции Ричардсона [4].
Результаты расчетов при различных значениях скоростей переноса давали значение погрешности, отличающееся на величину ошибок машинного округления
2=1
2=1
Рис. 1. Зависимость погрешности численного решения от шага дискретизации
для задачи (4, 5)
Для иллюстрации на рис. 1 и 2 приведены зависимости нормы погрешности численного решения ||иг.|| от шага дискретизации в логарифмических координатах для задач (4, 5) и (4, 6).
Рис. 2. Зависимость погрешности численного решения от шага дискретизации для задачи (4, 6)
Здесь к - угловой коэффициент прямых, аппроксимирующих зависимости нормы погрешности от шага дискретизации. По сути, к есть величина порядка аппроксимации конечноразностной схемы.
Как видно из приведенных результатов, с уменьшением шага дискретизации погрешность приближенного решения экспоненциально уменьшается, а порядок аппроксимации асимптотически приближается к заявленному для обеих задач.
4. Заключение
В данной работе построены симметричные схемы второго, третьего, четвертого, пятого, шестого, седьмого и восьмого порядков точности методом совместной аппроксимации для уравнения переноса. Схемы работают как при положительных, так и при отрицательных значениях параметра a. Результаты, полученные в статье, подтверждаются расчетными данными.
СПИСОК ЛИТЕРАТУРЫ
1. Бучарский В. Л. Метод совместной аппроксимации построения разностных схем для решения уравнений в частных производных / В. Л. Бучарский // Техническая механика. - 2007. - № 1. - С. 50
- 57.
2. Бучарский В.Л. Монотонная разностная схема для уравнений Эйлера / В.Л. Бучарский, В.Ф. Присняков // Проблемы высокотемпературной техники. - Д.: Вид-во ДНУ, 1994. - С. 65 - 71.
3. Бучарский В. Л. Разностная схема метода совместной аппроксимации для решения квазилинейных гиперболических уравнений / В. Л. Бучарский // Проблеми обчислювальної механіки і міцності конструкцій. - Дніпропетровськ: ДНУ, 2008 . - C. 7 - 15.
4. Шокин Ю.И. Метод дифференциального приближения. Приложение к газовой динамике / Ю.И. Шокин, Н.Н. Яненко. - Новосибирск: Наука, 1985. - 364 с.
5. Encyclopedia of Computational Mechanics. Volume 1 Fundamentals / Ed. E. Stein, R. de Borst, T.J. R. Hughes. - WILLEY, 2004. - 798 p.
Стаття надійшла до редакції 27.09.2010
