Three-dimensional numerical simulation of the inverse problem of thermal convection using the quasi-reversibility method

Inverse (time-reverse) simulation of three-dimensional thermoconvective flows is considered for a highly viscous incompressible fluid with temperature-dependent density and viscosity. The model of the fluid dynamics is described by the Stokes equations, the incompressibility and heat balance equations subject to the appropriate initial and boundary conditions. To solve the problem backward in time, the quasi-reversibility method is applied to the heat balance equation. The numerical solution is based on the introduction of a two-component vector potential for the velocity of the medium, on the application of the finite element method with a special tricubic spline basis for computing this potential, and on the application of the splitting method and the method of characteristics for computing the temperature. The numerical algorithm is designed to be executed on parallel computers. The proposed numerical algorithm is used to reconstruct the evolution of diapiric structures in the Earth's upper mantle. The computational efficiency of the algorithm is analyzed on the basis of the appropriate functionals of residuals.