1IPGP - Institut de Physique du Globe de Paris (IPGP, 1 rue Jussieu, 75238 Paris cedex 05 ; Université Paris Diderot, Bât. Lamarck A case postale 7011, 75205 Paris CEDEX 13 - France)
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.
https://hal-insu.archives-ouvertes.fr/insu-03600288 Contributor : Nathalie PothierConnect in order to contact the contributor Submitted on : Monday, March 7, 2022 - 2:57:49 PM Last modification on : Wednesday, March 9, 2022 - 3:36:00 AM
A. Fournier, H. -P. Bunge, R. Hollerbach, J. -P. Vilotte. Application of the spectral-element method to the axisymmetric Navier-Stokes equation. Geophysical Journal International, Oxford University Press (OUP), 2004, 156, pp.682-700. ⟨10.1111/j.1365-246X.2004.02149.x⟩. ⟨insu-03600288⟩