В работе рассматриваются два метода численного решения уравнений движения несжимаемой жидкости.

Одним из методов является метод искусственной сжимаемости. Суть метода заключается во введении дополнительного члена (1/β)*(∂p/∂τ) (производная по псевдо времени от давления) в уравнение неразрывности, где β = const – коэффициент искусственной сжимаемости. Для аппроксимации невязких потоков в полученной системе используется противопотоковая схема 3-го порядка с расщеплением матрицы Якоби невязкого потока на сумму матриц A⁺ и A^– , обладающих неотрицательными и неположительными собственными значениями.

Второй рассмотренный метод - метод coupled, в котором для аппроксимации невязких потоков используется противопотоковая схема, заимствованная из метода SIMPLE, где предварительно вычисляется поток массы через грань ячейки. Для расчета потока массы используется Рхи—Чоу интерполяция.

В результате линеаризации схем обоих методов получена СЛАУ Ax = b, где A – это разреженная матрица с блоками 4х4. В работе рассматривается два варианта решения данной СЛАУ - приближенная LU - факторизация неявного оператора, а также решение без факторизации с использованием итерационного метода GMRES из библиотеки MKL. Произведено сравнение методов по скорости сходимости и численному решению на задачах невязкого и вязкого обтекания цилиндра.

Семинар будет проведен в смешанном формате. Очное заседание пройдет в конференц-зале ФИЦ ИВТ (к.513), дистанционно к семинару можно будет подключиться по ссылке https://vcs-6.ict.nsc.ru/rooms/grz-ayy-7ne/join

Для студентов и аспирантов, находящихся в Новосибирске и не имеющих причин для дистанционного подключения, очное присутствие обязательно.