Примеры разностных аппроксимаций

Для уравнения (9) неравенство |q| 1 выполняется согласно (11) при всех тогда и только тогда, когда

Примеры разностных аппроксимаций

Информация

Математика и статистика

Другие материалы по предмету

Математика и статистика

Сдать работу со 100% гаранией

Иногда для упрощения записи индексы i и n будем опускать, обозначая

 

 

 

(xi, tn+1)(xi-1, tn+1)(xi, tn+1)(xi+1, tn+1)

 

 

 

(xi-1, tn)(xi, tn)(xi+1, tn)(xi, tn)

 

(xi-1, tn+1)(xi, tn+1)(xi+1, tn+1)(xi, tn+1)

 

 

 

(xi-1, tn)(xi, tn)(xi+1, tn)(xi-1, tn)(xi, tn)(xi+1, tn)

 

(xi, tn-1)

 

 

Чтобы аппроксимировать уравнение (1) в точке (xi, tn), введем шаблон, изображенный на рисунке и состоящий из четырех узлов (xi1, tn), (xi, tn), (xi, tn+1). Производную u/t заменим в точке (xi, tn) разностным отношением ynt, i, а производную 2u/2x второй разностной производной ynxx, i. Правую часть f(x, t) заменим приближенно сеточной функцией ni, в качестве ni можно взять одно из следующих выражений:

 

 

 

В результате получим разносное уравнение

 

(5)

 

которое аппроксимирует исходное дифференциальное уравнение в точке (xi, tn) с первым порядком по и вторым порядком по h при условии, что разность ni f(xi, tn) имеет тот же порядок малости.

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

 

 

(6)

 

 

 

 

Эта схема представляет собой систему линейных алгебраических уравнений с числом уравнений, равным числу неизвестных. Находить решение такой системы следует по слоям. Решение на нулевом слое задано начальными условиями y0i = u0(xi), i = 0, 1,…, N. Если решение yni, i = 0, 1,…, N, на слое n уже найдено, то решение yin+1 на слое n+1 находится по явной формуле

 

(7)

 

а значениядоопределяются из граничных

условий. По этой причине схема (6) называется явной разностной схемой. Несколько позже мы познакомимся и с неявными схемами, в которых для нахождения yin+1 при заданных yin требуется решать систему уравнений.

Погрешность разностной схемы (6) определяется как разность zin = yin u(xi, tn) между решением задачи (6) и решением исходной задачи (1) (3). Подставляя в (6) yin = zin + u(xi, tn), получим уравнение для погрешности

 

 

(8)

 

где погрешность аппроксимации разностной

схемы (6) на решении задачи (1) (3), in = O( + h2). Можно оценить решение zin уравнения (8) через правую часть in и доказать тем самым сходимость разностной схемы (6) с первым порядком по и вторым по h. Однако это исследование мы отложим, а сейчас на примере схемы (6) продемонстрируем один распространенный прием исследования разностных схем с постоянными коэффициентами, называемый методом гармоник. Хотя данный метод не является достаточно обоснованным, в частности не учитывает влияния граничных условий и правых частей, он позволяет легко найти необходимые условия устойчивости и сходимости разностных схем. Покажем, например, что явную схему (6) можно применять лишь при условии 0,5h2, означающем, что шаг по времени надо брать достаточно малым.

Рассмотрим уравнение

 

(9)

 

т.е. однородное уравнение, соответствующее (5). Будем искать частные решения (9), имеющие вид

 

yjn () = qneijh,(10)

 

где i мнимая единица, любое действительное число и q число, подлежащее определению. Подставляя (10) в уравнение (9) и сокращая на eijh, получим

 

 

 

откуда найдем

 

(11)

 

Начальные условиясоответствующие решениям вида (10) (их называют гармониками), ограничены. Если для некоторого множитель q станет по модулю больше единицы, то решение вида (10) будет неограниченно возрастать при n. В этом случае разностное уравнение (9) называется неустойчивым, поскольку нарушается непрерывная зависимость его решения от начальных условий. Если же |q| 1 для всех действительных , то все решения вида (10) ограничены при любом n и разностное уравнение (9) называется устойчивым. В случае неустойчивости найти решение разностной задачи (6) по формулам (7) практически невозможно, так как погрешности (например погрешности округления), внесенные в начальный момент времени, будут неограниченно возрастать при увеличении n. Такие разностные схемы называются неустойчивыми.

Для уравнения (9) неравенство |q| 1 выполняется согласно (11) при всех тогда и только тогда, когда 0,5. Таким образом, использование схемы (6) возможно лишь при выполнении условия 0,5h2. Разностные схемы, устойчивые лишь при некотором ограничении на отношение шагов по пространству и по времени, называются условно устойчивыми. Следовательно, схема (6) возможно устойчива, причем условие устойчивости имеет вид /h2 0,5. Условно устойчивые схемы для уравнений параболического типа используются редко, так как они накладывают слишком сильное ограничение на шаг по времени. Действительно, пусть, например, h = 10-2. Тогда шаг не должен превосходить 0,5 * 10-4, и для того чтобы вычислить решение yjn при t = 1, надо взять число шагов по времени n = -1 2 * 104, т.е. провести не менее 2 * 104 вычислений по формулам (7).

 

3.3. Неявные схемы. Чисто неявной разностной схемой для уравнения теплопроводности теплопроводности (схемой с опережением) называется разностная схема, использующая шаблон (xi, tn), (xi1, tn+1), (xi, tn+1) и имеющая вид

 

(12)

 

 

 

 

Здесь ni = f(xi, tn+1) + O( + h2). Схема имеет первый порядок аппроксимации по и второй по h. Решение системы (12) находится, как и в случае явной схемы, по слоям, начиная с n = 1. Однако, теперь, в отличие от явной схемы, для нахождения yin+1 по известным yin требуется решить систему уравнений

 

(13)

 

где = /h2, Fin = yin + in. Эту систему можно решать методом прогонки, так как условия устойчивости прогонки выполнены.

Для исследования устойчивости разностной схемы (12) будем искать частные решения уравнения

 

 

 

имеющие вид (10). Тогда получим

 

 

 

следовательно, |q| 1 при любых , , h. Таким образом, схема (12) абсолютно устойчива, т.е. устойчива при любых шагах и h. Абсолютная устойчивость является основным условием неявных схем. Теперь уже не надо брать шаг слишком малым, можно взять, например, = h = 10-2. Величина шагов сетки , h определяются теперь необходимой точностью расчета, а не соображениями устойчивости.

Шеститочечной симметричной схемой называется разностная схема

(14)

 

для которой начальные и граничные условия задаются так же, как и в схеме (12). Эта схема использует шеститочечный шаблон, изображенный на рисунке.

Обобщением трех рассмотренных схем является однопараметрическое семейство схем с весами. Зададим произвольный действительный параметр и определим разностную схему

 

(15)

 

 

 

 

При = 0 получим отсюда явную схему, при = 1 чисто неявную схему и при = 0,5 симметричную схему (14). Исследуем погрешность аппроксимации схемы (15) на решении исходной задачи (1) (3). Представим решение задачи (15) в виде yin = u(xi, tn) + zin, где u(xi, tn) точное решение дифференциальной задачи (1) (3). Тогда для погрешности получим систему уравнений

 

(16)

 

i = 1, 2,…, N 1, n = 0, 1,…, K 1,

z0n+1 = zNn+1 = 0, n = 0, 1,…, K 1, zi0 = 0, i = 0, 1,…, N.

 

Сеточная функция in, входящая в правую часть уравнения (16) и равная

 

(17)

называется погрешностью аппроксимации схемы (15) на решении з

Похожие работы

<< < 1 2 3 4 >