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 -