Abstract:
We present an application of the spectral-element method to model axisymmetric flows in rapidly rotating domains. The primitive equations are discretized in space with local tensorized bases of high-order polynomials and in time with a second-order accurate scheme that treats viscous and rotational effects implicitly. We handle the pole problem using a weighted quadrature in elements adjacent to the axis of rotation. The resulting algebraic systems are solved efficiently using preconditioned iterative procedures. We validate our implementation through comparisons with analytic and purely spectral solutions to laminar flows in a spherical shell. This axisymmetric tool is the kernel on which complexity will be added subsequently in the long-term prospect of building a parallel spectral-element based geodynamo model.
