Аннотация:
Разработан и реализован бессеточный метод для трехмерного численного решения системы нестационарных уравнений Навье–Стокса. Метод основан на дискретизации области конечным набором распределенных вычислительных узлов. Для повышения точности применяется комбинированная аппроксимация пространственных производных: для конвективных потоков используется полиномиальный метод наименьших квадратов (Polynomial Least Squares, PLS), для вязких — разложение в ряд Тейлора (Taylor Least Squares, TLS). Ключевой особенностью, устраняющей асимметрию при расчете обтекания осесимметричных тел, является преобразование ортонормированной системы координат для каждой пары узлов при вычислении конвективных потоков. Реконструкция векторов состояния по схеме MUSCL и векторов градиента обеспечивает второй порядок точности по пространству. Интегрирование по времени выполняется явным методом Рунге–Кутты. Программная реализация на C++ с использованием OpenCL позволяет проводить расчеты на графических процессорах. Метод апробирован на задаче сверхзвукового обтекания сферы; результаты демонстрируют хорошее согласование с эталонными данными, а отклонение конвективного теплового потока при увеличении числа узлов до $2.5 \cdot 10^7$ не превышает 2 %.