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