Modeling Upper Mantle Convection in the Subduction Zone

A model of upper mantle convection in the subduction zone of a cold lithospheric plate (subduction) into the Earth’s upper strata is developed. The issues of constructing the initial distributions of model variables are discussed. Computational schemes for solving the model equations are given. Calculation of dynamics of mantle convection and reorganization of its structure are performed in the vorticity-current function variables, and dynamics of the plate subduction is calculated on the basis of the smoothed-particle hydrodynamics method (SPH). A series of computational experiments are performed.