Численное моделирование течения жидкости в плоской каверне при больших числах Рейнольдса

В работе обсуждается вопрос существования стационарного решения задачи, в которой исследуется течение несжимаемой вязкой жидкости в двумерной квадратной каверне с подвижной верхней крышкой при числах Рейнольдса Re>10000. На основании обзора литературных источников и анализа численного решения задачи при 10000≤Re≤20000 делаются выводы об условиях и особенностях построения требуемого решения. Достоверность полученных результатов подтверждается хорошим совпадением расчетных параметров вихрей вплоть до четвертого уровня с аналогичными данными других авторов при Re=10000 и Re=20000. Для этих же чисел Рейнольдса приводится подробная структура течения. Демонстрируется эволюция нижних угловых вихрей четвертого уровня для диапазона 15000≤Re≤20000. Показывается, что при Re>10000 течение жидкости в вихрях третьего и четвертого уровней носит выраженный стоксовский характер. Проведено сравнение значений параметров вихрей в углах каверны, полученных численно и аналитически на основании решения Моффатта, описывающего стоксовское течение жидкости около острых углов. Выявлено их вполне удовлетворительное соответствие. Обоснование корректности найденного численно стационарного решения при больших числах Рейнольдса осуществляется с точки зрения его насыщения при кратном возрастании количества расчетных узлов вдоль каждой пространственной координаты, влияния схемной вязкости, а также рассмотрения процесса установления решения в контексте поведения динамических систем.