Метод коллокаций и наименьших невязок для трехмерных уравнений Навье-Стокса

Метод коллокаций и наименьших невязок, предложенный ранее для численного решения двумерных уравнений Навье-Стокса, описывающих стационарные течения вязкой несжимаемой жидкости, обобщен на трехмерный случай. В реализованном варианте метода решение ищется в виде разложения по базисным соленоидальным функциям. Для коэффициентов разложения в каждой ячейке расчетной сетки получается переопределенная система линейных алгебраических уравнений, которая решается методом вращений. Для ускорения сходимости итерационного процесса предложен новый алгоритм, основанный на подпространствах Крылова. Результаты верификации метода подтверждают второй порядок сходимости для составляющих вектора скорости. Представлены результаты решения эталонной задачи о течении в кубической каверне с движущейся крышкой для чисел Рейнольдса ${\rm Re} = 100$ и ${\rm Re} = 1000$. Показано, что полученные результаты весьма близки по точности к наиболее точным результатам, полученным другими авторами с помощью различных численных методов высокой точности. Работа выполнена при частичной финансовой поддержке РФФИ (код проекта 13-01-00227).