Рассмотрена конвекция вещества верхней мантии Земли, которая в приближении Обербека–Буссинеска обусловлена термогравитационной дифференциацией. В рамках этого приближения выполнено 2-D численное моделирование конвективных течений вещества среды. Уравнение для температуры следует из соотношения баланса энтропии, где вследствие учета в системе переменной вязкости присутствует эффект диссипации энергии. Краевые условия отвечают заданию общепринятой на границе верхней и нижней мантий температуры, а для боковых границ – их теплоизолированность. На границе астеносфера–литосфера приняты допущения о том, что динамика тепла определяется его потоком с ближнего к границе слоя астеносферы, рассеиванием части тепла вдоль границы и расходами тепла на плавление вещества литосферы. Численное решение определяющих уравнений выполнено в переменных функция тока – завихренность. Приведена итерационная схема их решения. Обсуждены вопросы программной реализации аппарата численного моделирования. Показано, что при таких краевых условиях в рассматриваемой системе формируется квазипериодических режим колебаний тепла.
The convection of matter in the Earth's upper mantle is considered, which in the Oberbeck–Boussinesq approximation is due to thermogravitational differentiation. Within the framework of this approximation, a 2-D numerical simulation of convective flows of the medium matter was performed. The equation for temperature follows from the entropy balance relation, where, due to taking into account the variable viscosity in the system, there is an effect of energy dissipation. The boundary conditions correspond to the assignment of the temperature generally accepted at the boundary of the upper and lower mantles, and for the lateral boundaries - their thermal insulation. At the asthenosphere–lithosphere boundary, assumptions were made that the heat dynamics is determined by its flow from the asthenosphere layer closest to the boundary, part of the heat dissipation along the boundary, and heat consumption for melting the lithosphere matter. Numerical solution of the constitutive equations is carried out in variables stream function - vorticity. An iterative scheme for their solution is given. The issues of software implementation of the numerical simulation apparatus are discussed. It is shown that under such boundary conditions, a quasi-periodic regime of heat oscillations is formed in the system under consideration.
The authors declare that there are no conflicts of interest present.