Шпоры по вычмату / 019

Неявная разностная схема для уравнения Пуассона.

Рассмотрим уравнение Пуассона для двух переменных , в прямоугольной области с размерами сторон а, b

10.1

Обычно для уравнения Пуассона задаются граничные условия, т.е. значения функции "" на сторонах прямоугольника считаем известными. Нужно определить значения внутри прямоугольника, см. рис. 10.1.

Пусть шаг дискретизации по обеим координатам одинаков и равен . Для вторых производных выберем формулу численного дифференцирования (4.4), что дает разностную схему

10.2

Этой разностной схеме соответствует пятиточечный шаблон рис.10.1, который называют ''крестом''.

Рис. 12.1. Шаблон "крест" для уравнения Пуассона.

В какое бы место рассматриваемой области, см. рис. 11.1, мы не поместили этот шаблон, уравнение (10.2) будет содержать три или более неизвестных. Это не позволяет явно вычислять значения функции "" на следующем слое, как это выполнялось для явной схемы (11.8). Поэтому разностная схема (10.2) называется неявной.

При трех переменных , , «крест» имеет еще одно направление перпендикулярно плоскости. Чтобы найти неизвестные значения в случае неявной схемы, необходимо записать уравнение (10.2) для каждого внутреннего узла рассматриваемой сетки. В результате получим систему линейных уравнений, в которой количество неизвестных равно количеству уравнений. Например, если имеем внутренних узлов по и внутренних узлов по y, то количество неизвестных и уравнений равно , т.е. обычно очень велико.

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

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

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

На рис.10.2 для случая двух независимых переменных показано как зависит погрешность решения от шага дискретизации по одной переменной, а шаг дискретизации по другой переменной является параметром для каждой кривой. Возрастание погрешности при уменьшении шага для явных схем соответствует неустойчивости. Отметим, что в случае неявных схем при возможно проявление погрешностей округления при численном дифференцировании, рассмотренных в лекции 4, но столь малые шаги при решении уравнений в частных производных не используются.

Рис. 12.2 Зависимость погрешности от шагов дискретизации для явных и неявных разностных схем.

Контроль точности

Пусть при произвольно выбранных шагах дискретизации и получено решение уравнения в частных производных в узлах сетки рис. 11.1. Для контроля, как и в лекции 7 при решении ОДУ, нужно выполнить повторный расчет с другими шагами или другим методом. Обычно расчет повторяют с шагами , и сравнивают два решения в узлах первой (более крупной) сетки. За оценку погрешности можно взять значение

,

10.3

где номера узлов первой сетки, – первое решение, – второе решение в тех же точках .

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