р
ур
у
д (р р ),
д
единицу объёма (Вт/м°).
Рис. 1 схематично иллюстрирует расчётную область и постановку 
нижней и верхней поверхностях задаётся избыточная по отношению к 
уровню температура Т = () и условие прилипания, на боковых границах задаё
периодичности (как для скорости, так и для температуры).
Рис. 1. Общий вид расчетной области
Решаемая система уравнений в безразмерной форме относительн
скорости и, давления р и избыточной температуры / представлена ниже (1).
И -и= 0
ди
д—[+и' Йи= — 7р+РГ [72Ц+РГ’ЧЗТ2 
ЭТ уи.рт- Т-
дЁ
Во всех расчетах в данной задаче принималось число Прандтля Рг = у/а
Рэлея в данном исследовании определялось в соответствии с работой (Со!и5
Рое!, 2016) как Ка = 2ВЕа/(рСра’у). В расчетах использовались следующи
числа Рэлея: Ва = 10°; 10% 10'°. Указанные значения числа Рэлея соо
турбулентному режиму течения. В расчётах использовались модели турбулент
КАМЗ, а именно: модель рейнольдсовых напряжений ГЁК (Гаппет, 1975), 
$5УТ, модель Спаларта-Аллмараса.
Вычислительные аспекты. Для проведения численных расчетов испо
программное обеспечение с открытым исходным кодом ОреоЕОАМ у11. В дан
использовалась структурированная сетка с прямоугольными ячейками со сгу
холодным стенкам. Размеры сетки составляли 50 х 50 х 100 узлов, со сг
нижней и верхней стенкам.
Для повышения устойчивости задачи и обеспечения сходимости, 
решатель /Лшй, использующийся для моделирования такого рода задач, был
изменения для уравнения энергии и приведения его к виду (1).
Для определения времени установления в центре расчетной обл
установлена точка мониторинга температуры, по которой при появлении х
пульсаций можно говорить о развитом режиме конвекционного течения.
289