Методы решения нестационарных задач

Обратимся к рассмотренным выше аппроксимациям нестационарных задач. Так, для краевой задачи

Методы решения нестационарных задач - student2.ru в Методы решения нестационарных задач - student2.ru , Методы решения нестационарных задач - student2.ru ,

Методы решения нестационарных задач - student2.ru , Методы решения нестационарных задач - student2.ru в Методы решения нестационарных задач - student2.ru

явная разностная схема имеет вид

Методы решения нестационарных задач - student2.ru , Методы решения нестационарных задач - student2.ru , (2.47)

Методы решения нестационарных задач - student2.ru .

В чем состоит метод решения системы сеточных уравнений (2.47) – в проведении вычислений по простым рекуррентным формулам.

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

Методы решения нестационарных задач - student2.ru , Методы решения нестационарных задач - student2.ru . (2.48)

Таким образом, для реализации шага неявной схемы следует решить систему сеточных уравнений вида

Методы решения нестационарных задач - student2.ru . (2.49)

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

Если область, для которой решается задача, является прямоугольником, то существуют специфические весьма эффективные методы решения системы (2.49): прямые методы – быстрое дискретное преобразование Фурье, метод циклической редукции; итерационные методы – метод переменных направлений с оптимальным выбором параметров, метод Р.П. Федоренко, или многосеточный метод (multigrid) и др.

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

Продемонстрируем построение этих схем на примере следующей задачи:

Методы решения нестационарных задач - student2.ru в Методы решения нестационарных задач - student2.ru ,

Методы решения нестационарных задач - student2.ru в Методы решения нестационарных задач - student2.ru .

Предположим, что оператор Методы решения нестационарных задач - student2.ru не зависит от времени. Здесь мы считаем, что оператор Методы решения нестационарных задач - student2.ru рассматривается на функциях, удовлетворяющих краевым условиям.

Относительно оператора Методы решения нестационарных задач - student2.ru предположим, что

Методы решения нестационарных задач - student2.ru .

Предположим также, что решение систем сеточных уравнений, матрицами которых являются соответствующие сеточные операторы Методы решения нестационарных задач - student2.ru и Методы решения нестационарных задач - student2.ru , легко осуществляется и не требует больших вычислительных затрат.

В этом случае удобно рассмотреть следующую схему:

Методы решения нестационарных задач - student2.ru ,

Методы решения нестационарных задач - student2.ru , Методы решения нестационарных задач - student2.ru .

Это схема переменных направлений, или продольно-поперечная схема. Происхождение названия схемы становится понятным, если применить ее к решению задачи

Методы решения нестационарных задач - student2.ru в Методы решения нестационарных задач - student2.ru ,

Методы решения нестационарных задач - student2.ru на Методы решения нестационарных задач - student2.ru ,

Методы решения нестационарных задач - student2.ru в Методы решения нестационарных задач - student2.ru при Методы решения нестационарных задач - student2.ru .

Для данной задачи схема примет вид

Методы решения нестационарных задач - student2.ru , Методы решения нестационарных задач - student2.ru ,

Методы решения нестационарных задач - student2.ru , Методы решения нестационарных задач - student2.ru .

В первом из этих уравнений схема является явной по направлению переменной Методы решения нестационарных задач - student2.ru и неявной по направлению переменной Методы решения нестационарных задач - student2.ru , а во втором уравнении – наоборот.

Еще одним примером широко используемых схем для решения эволюционных задач являются схемы, которые получаются так называемым методом расщепления.

Приведем пример схемы, полученной методом расщепления:

Методы решения нестационарных задач - student2.ru ,

Методы решения нестационарных задач - student2.ru , Методы решения нестационарных задач - student2.ru .

Для рассмотренного примера, где Методы решения нестационарных задач - student2.ru , имеем

Методы решения нестационарных задач - student2.ru ,

Методы решения нестационарных задач - student2.ru , Методы решения нестационарных задач - student2.ru .

У этой схемы много названий: локально-одномерная, дробных шагов, покоординатного расщепления.

Обе рассмотренные схемы при Методы решения нестационарных задач - student2.ru и Методы решения нестационарных задач - student2.ru абсолютно устойчивы.

Наши рекомендации