Приведены результаты численного моделирования теплообмена при обтекании кубических каверн. Исследовано влияние различных параметров ламинарного набегающего потока на устойчивость течения (частота возмущений, толщина пограничного слоя, скорость основного потока). Показано, что интегральные потоки тепла существенным образом зависят от частот возмущений в набегающем потоке, причем частоты, приводящие к интенсификации теплообмена, определяются устойчивостью слоя смешения между каверной и основным потоком. Проведено сопоставление расчетных результатов интегральных характеристик теплообмена с экспериментальными данными. Получено хорошее согласование зависимости числа Нуссельта от числа Рейнольдса с экспериментом.