При аппроксимации оператора переноса B в уравнении (4) считаем положительными обе компоненты скорости. При иных знаках

изменения очевидны.

Полагаем функциюпостоянной в ячейке (n,m) в том смысле, что именно это значение переносится через правую и верхнюю границы ячейки, то есть

Соответственно

3.3.10. Метод решения дискретного уравнения переноса

Обращение начинаем с левого нижнего угла сетки, так как полагаем u > 0, v > 0. Если знаки u и v другие, то начинаем с других углов.

При нахожденииуже известны (вычислены на предыдущих шагах), а значит, можно определить инаходим из уравнения

Для контроля за точностью решения будем вычислять норму невязки уравнения (3). Итерационный процесс прекращается, когда норма невязки уменьшилась в требуемое количество раз, например, в 1000 или 10000.