The viability of using collocation methods in radius and spherical harmonics in the angular variables to calculate convective flows in full spherical geometry is examined. As a test problem the stability of the conductive state of a self-gravitating fluid sphere subject to rotation and internal heating is considered. A study of the behavior of different radial meshes previously used by several authors in polar coordinates, including or not the origin, is first performed. The presence of spurious modes due to the treatment of the singularity at the origin, to the spherical harmonics truncation, and to the initialization of the eigenvalue solver is shown, and ways to eliminate them are presented. Finally, to show the usefulness of the method, the neutral stability curves at very high Taylor and moderate and small Prandtl numbers are calculated and shown.
CitationSanchez, J., Garcia, F., Net, M. Radial collocation methods for the onset of convection in rotating spheres. "Journal of computational physics", 01 Març 2016, vol. 308, p. 273-288.
All rights reserved. This work is protected by the corresponding intellectual and industrial property rights. Without prejudice to any existing legal exemptions, reproduction, distribution, public communication or transformation of this work are prohibited without permission of the copyright holder. If you wish to make any use of the work not provided for in the law, please contact: firstname.lastname@example.org