A Non-oscillatory Advection Operator for the Compatible Spectral Element Method

Authors:
M. A. Taylor;A. St. Cyr;A. Fournier
Affiliations:
Sandia National Laboratories, Albuquerque, NM 87185;National Center for Atmospheric Research, Boulder, CO 80303;National Center for Atmospheric Research, Boulder, CO 80303
Venue:
ICCS 2009 Proceedings of the 9th International Conference on Computational Science
Year:
2009

Citing 5
Cited 1

The spectral element method for the shallow water equations on the sphere

Journal of Computational Physics
A 14.6 billion degrees of freedom, 5 teraflops, 2.5 terabyte earthquake simulation on the Earth Simulator

Proceedings of the 2003 ACM/IEEE conference on Supercomputing
High-Resolution Mesh Convergence Properties and Parallel Efficiency of a Spectral Element Atmospheric Dynamical Core

International Journal of High Performance Computing Applications
The NCAR Spectral Element Climate Dynamical Core: Semi-Implicit Eulerian Formulation

Journal of Scientific Computing
A standard test set for numerical approximations to the shallow water equations in spherical geometry

Journal of Computational Physics

CAM-SE: A scalable spectral element dynamical core for the Community Atmosphere Model

International Journal of High Performance Computing Applications

Quantified Score

Hi-index	0.00

Visualization

Abstract

The spectral element method is well known as an efficient way to obtain high-order numerical solutions on unstructured finite element grids. However, the oscillatory nature of the method's advection operator makes it unsuitable for many applications. One popular way to address this problem is with high-order discontinuous-Galerkin methods. In this work, an alternative solution which fits within the continuous Galerkin formulation of the spectral element method is proposed. Making use of a compatible formulation of spectral elements, a natural way to implement conservative non-oscillatory reconstructions for spectral element advection is shown. The reconstructions are local to the element and thus preserve the parallel efficiency of the method. Numerical results from a low-order quasi-monotone reconstruction and a higher-order sign-preserving reconstruction are presented.