TY - JOUR

T1 - Stabilization of the spectral element method in convection dominated flows by recovery of skew-symmetry

AU - Malm, Johan

AU - Schlatter, Philipp

AU - Fischer, Paul F.

AU - Henningson, Dan S.

N1 - Funding Information:
Acknowledgments The authors gratefully acknowledge funding by VR (The Swedish Research Council) Computer time was provided by SNIC (Swedish National Infrastructure for Computing) with a generous grant by the Knut and Alice Wallenberg (KAW) Foundation. The simulations were run at the Centre for Parallel Computers (PDC) at the Royal Institute of Technology (KTH). The third author was supported by the U.S. Department of Energy under Contract DE-AC02-06CH11357.

PY - 2013/11

Y1 - 2013/11

N2 - We investigate stability properties of the spectral element method for advection dominated incompressible flows. In particular, properties of the widely used convective form of the nonlinear term are studied. We remark that problems which are usually associated with the nonlinearity of the governing Navier-Stokes equations also arise in linear scalar transport problems, which implicates advection rather than nonlinearity as a source of difficulty. Thus, errors arising from insufficient quadrature of the convective term, commonly referred to as 'aliasing errors', destroy the skew-symmetric properties of the convection operator. Recovery of skew-symmetry can be efficiently achieved by the use of over-integration. Moreover, we demonstrate that the stability problems are not simply connected to underresolution. We combine theory with analysis of the linear advection-diffusion equation in 2D and simulations of the incompressible Navier-Stokes equations in 2D of thin shear layers at a very high Reynolds number and in 3D of turbulent and transitional channel flow at moderate Reynolds number. For the Navier-Stokes equations, where the divergence-free constraint needs to be enforced iteratively to a certain accuracy, small divergence errors can be detrimental to the stability of the method and it is therefore advised to use additional stabilization (e.g. so-called filter-based stabilization, spectral vanishing viscosity or entropy viscosity) in order to assure a stable spectral element method.

AB - We investigate stability properties of the spectral element method for advection dominated incompressible flows. In particular, properties of the widely used convective form of the nonlinear term are studied. We remark that problems which are usually associated with the nonlinearity of the governing Navier-Stokes equations also arise in linear scalar transport problems, which implicates advection rather than nonlinearity as a source of difficulty. Thus, errors arising from insufficient quadrature of the convective term, commonly referred to as 'aliasing errors', destroy the skew-symmetric properties of the convection operator. Recovery of skew-symmetry can be efficiently achieved by the use of over-integration. Moreover, we demonstrate that the stability problems are not simply connected to underresolution. We combine theory with analysis of the linear advection-diffusion equation in 2D and simulations of the incompressible Navier-Stokes equations in 2D of thin shear layers at a very high Reynolds number and in 3D of turbulent and transitional channel flow at moderate Reynolds number. For the Navier-Stokes equations, where the divergence-free constraint needs to be enforced iteratively to a certain accuracy, small divergence errors can be detrimental to the stability of the method and it is therefore advised to use additional stabilization (e.g. so-called filter-based stabilization, spectral vanishing viscosity or entropy viscosity) in order to assure a stable spectral element method.

KW - Over-integration

KW - Skew-symmetry

KW - Spectral element method (SEM)

KW - Stability

UR - http://www.scopus.com/inward/record.url?scp=84885471926&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84885471926&partnerID=8YFLogxK

U2 - 10.1007/s10915-013-9704-1

DO - 10.1007/s10915-013-9704-1

M3 - Article

AN - SCOPUS:84885471926

SN - 0885-7474

VL - 57

SP - 254

EP - 277

JO - Journal of Scientific Computing

JF - Journal of Scientific Computing

IS - 2

ER -