В рамках приближении Обербека—Буссинеска рассмотрена термогравитационная 2Б-конвекция в верхней мантии Земли. На границе астеносфера—литосфера принимаются допущения о том, что динамика теплоты определяется ее потоком с ближнего к границе слоя астеносферы, рассеиванием части теплоты вдоль границы и расходами теплоты на плавление вещества литосферы. Численное решение определяющих уравнений выполняется в переменных функция тока — завихренность. Приводится итерационная схема их решения. Обсуждаются вопросы программной реализации аппарата численного моделирования. Показано, что при таких краевых условиях в рассматриваемой системе формируется квазипериодический режим колебаний теплоты. Ключевые слова: моделирование, тепловой поток, конвекция
