Нужна помощь с методом МКО для решения уравнения Навье-Стокса

Заранее прошу прощения я не математик, но род действий вынуждает с ней сталкиваться.
Я решил уравнения Навье-Стокса для вязкой не сжимаемой жидкости при помощи МКР. Однако из-за сложной конфигурации области, в которой происходит движение жидкости, решение является не устойчивым. Уменьшил шаг по времени однако это не помогло. Хочу решить данную задачу методом МКО на шахматной сетке но не могу разобраться как это сделать. Буду благодарен если кто-нибудь скинет пример реализации на языке паскаль или с++ или даст ссылки на толковую литературу. Буду рад любой помощи!