Three-dimensional convection in spherical shells

Abstract

Abstract In this study, the equations of the three-dimensional convective motion of an infinite Prandtl number fluid are solved in spherical geometry, for Rayleigh numbers up to 15 times the critical number. An iterative method is used to find stationary solutions. The spherical parts of the operators are treated using a Galerkin collocation method while the radial and time dependences are expressed using finite difference methods. A systematic search for stationary solutions has led to eight different stream patterns for a low Rayleigh number (1.28 times the critical number). They can be classified as: I) Axisymmetrical solutions, analogous to rolls in plane geometry. II) Solutions which have several ascending plumes within a large area of ascending current, and also several descending plumes within an area of descending current. This type of flow is analogous to bimodal circulation in plane geometry. III) Solutions characterized by isolated ascending (or descending) plumes separated from each other by a closed polyhedral network of descending (or ascending) currents. This type of circulation is called ‘polygonal’ in analogy with hexagonal circulation in plane geometry. The behaviour of each of the eight solutions has been studied by increasing the Rayleigh number up to 15 times the critical number. A trend towards transitions from type (I) and type (II) solutions to type (III) solutions is observed. It is inferred that only the “polygonal” solutions are stable for a Rayleigh number greater than 15 times the critical number.