УДК 531.51 УЧЕТ МОЛЕКУЛЯРНОЙ ВЯЗКОСТИ ПРИ ПРЯМОМ 3D ЧИСЛЕННОМ МОДЕЛИРОВАНИИ ГРАВИТАЦИОННОГО ТУРБУЛЕНТНОГО ПЕРЕМЕШИВАНИЯ А. Л. <...> Cтадник, В. П. Стаценко, Ю. В. Янилкин РФЯЦ-ВНИИЭФ С помощью газодинамической 3D методики ТРЭК исследуeтся влияние молекулярной вязкости на развитие гравитационного турбулентного перемешивания на плоской границе раздела двух несжимаемых жидкостей (газов). <...> Турбулентное перемешивание при неустойчивости Рэлея – Тейлора представляет собой классическую задачу. <...> Имеется достаточно большое количество экспериментальных данных (например, [1, 2]) для данной задачи по законам роста зоны турбулентного перемешивания (ЗТП), которые однако дают заметно различающиеся друг от друга значения коэффициента роста зоны на участке, который предполагается автомодельным (0,04 < < 0,35). <...> В наших работах [3–5] было показано, что в 3D расчетах по коду ТРЭК без учета молекулярной вязкости, как и в экспериментах, имеет место начальная неавтомодельная стадия развития турбулентного перемешивания. <...> Ширина соответствующего участка зависит от начальных данных, используемой разностной схемы и количества счетных ячеек в расчете, что указывает на значительное влияние схемной вязкости на решение. <...> В работе [7] на основе анализа экспериментальных данных [6, 8] показано, что на начальной стадии развития процесса перемешивания существенна роль молекулярной вязкости. <...> В настоящей работе приводятся результаты 3D расчетов с учетом молекулярной вязкости, для моделирования которой в комплексе ТРЭК была разработана методика решения уравнений Навье – Стокса для сжимаемой среды. <...> С использованием этой методики была выполнена серия расчетов на счетной сетке с числом счетных ячеек N = 2·106. <...> Результаты расчетов, проведенных по разработанной методике, сравниваются с аналогичными результатами без учета вязкости, что позволяет оценить как ее влияние на решение, так и величину схемных эффектов. <...> Ускорение тяжести gz = –1 – g направлено <...>