Решение дифференциальных уравнений эллиптического типа.
       5. Метод установления с использованием схемы переменных направлений.

     Запишем для уравнения (11.2) схему переменных направлений:
 
  (11.8)
Как и схема расщепления (11.7), схема переменных направлений (11.8) является абсолютно устойчивой, что обеспечивает возможность произвольного выбора значения шага итерации. Однако схема переменных направлений (11.8) имеет более высокий (второй) порядок аппроксимации по времени, что позволяет использовать при её решении более грубый шаг итерации, чем при решении схемы расщепления (11.7), для получения результатов с такой же точностью; а это в свою очередь способствует ускорению сходимости итерационного процесса и уменьшению количества итераций:
 

     Каждая из подсхем схемы переменных направлений (11.8) решается с помощью метода прогонки. Коэффициенты, соответствующие уравнению (4.10), имеют вид:

  • для первой подсхемы
     
  • для второй подсхемы
     
    Легко видеть, что для обеих подсхем достаточное условие сходимости прогонки (4.16) выполняется. Алгоритм решения схемы (11.8) аналогичен алгоритму решения схемы расщепления (7.7), (7.8), аппроксимирующей двумерное дифференциальное уравнение параболического типа; методики определения прогоночных коэффициентов и решения на правой границе аналогичны описанным ранее в разделах 4.2.3-4.2.4.
         Также как и в случае метода простой итерации, в качестве нулевой итерации (начального условия, необходимого для решения в связи с введением фиктивной производной по времени) обычно задают свободный член:
     
    Расчёт итераций следует продолжать до тех пор, пока итерационный процесс не сойдётся, т.е. пока не будет выполняться условие (11.3), в разностном представлении соответствующее неравенству:
     
    В качестве итерационного выражения служит прогоночное соотношение (4.11), имеющее вид:
  • для первой подсхемы
  • для второй подсхемы