Article ID Journal Published Year Pages File Type
10356409 Journal of Computational Physics 2005 28 Pages PDF
Abstract
We present a Fourier-spectral element approach for modeling thermal convection in a rotating axisymmetric container. Following the theory detailed in Bernardi et al. [C. Bernardi, M. Dauge, Y. Maday, Spectral Methods for Axisymmetric Domains, Gauthier-Villars, Paris, 1999], a Fourier expansion of the field variables is performed in the periodic direction, and the resulting collection of meridional problems is discretized by means of a parallel spectral element method. A Gauss-Lobatto-Jacobi (0,1) quadrature, which incorporates the cylindrical radius in its weight, is introduced to avoid a potential degeneracy of the discrete set of equations, due to those nodes located on the axis of symmetry. A second-order timestepping scheme is presented, which treats the Coriolis and viscous forces implicitly. Numerical comparisons with analytical and published numerical solutions in spherical and cylindrical geometries are presented which highlight the accuracy of the model and demonstrate its spectral convergence properties.
Related Topics
Physical Sciences and Engineering Computer Science Computer Science Applications
Authors
, , , ,