Application of supercomputers to 3-D mantle convection.

Current generation vector machines are providing for the first time the computing power needed to treat planetary mantle convection in a fully three-dimensional fashion. A numerical technique known as multigrid has been implemented in spherical geometry using a hierarchy of meshes constructed from the regular icosahedron to yield a highly efficient three-dimensional compressible Eulerian finite element hydrodynamics formulation. The paper describes the numerical method and presents convection solutions for the mantles of both the earth and the Moon. In the case of the Earth, the convection pattern is characterized by upwelling in narrow circular plumes originating at the core-mantle boundary and by downwelling in sheets or slabs derived from the cold upper boundary layer. The preferred number of plumes appears to be on the order of six or seven. For the Moon, the numerical results indicate that development of a predominately L = 2 pattern in later lunar history is a plausible explanation for the present large second-degree non-hydrostatic component in the lunar figure.