A numerical algorithm for simulation of two-component flows in the Earth's mantle

The equations of motion of a two-component isothermal viscous fluid for two cases are considered, when 1) velocity-pressure and 2) stream function - vortex are taken as unknown variables. Such equations can, in particular, be used for the mathematical simulation of some geodynamic processes. The algorithm of numerical solution of the two-dimensional equations is given for the following variables: a stream function, vortex, concentration. An example of the numerical calculation is presented.