Eigensolution analysis of spectral/hp continuous Galerkin approximations to advection-diffusion problems: insights into spectral vanishing viscosity
File(s)YJCPH_6288.pdf (5.38 MB) 1-s2.0-S0021999115008256-main.pdf (4.44 MB)
Accepted version
Published version
Author(s)
Moura, RC
Sherwin, SJ
Peiro, J
Type
Journal Article
Abstract
This study addresses linear dispersion–diffusion analysis for the spectral/hp continuous
Galerkin (CG) formulation in one dimension. First, numerical dispersion and diffusion
curves are obtained for the advection–diffusion problem and the role of multiple
eigencurves peculiar to spectral/hp methods is discussed. From the eigencurves’ behaviour,
we observe that CG might feature potentially undesirable non-smooth dispersion/diffusion
characteristics for under-resolved simulations of problems strongly dominated by either
convection or diffusion. Subsequently, the linear advection equation augmented with
spectral vanishing viscosity (SVV) is analysed. Dispersion and diffusion characteristics of
CG with SVV-based stabilization are verified to display similar non-smooth features in
flow regions where convection is much stronger than dissipation or vice-versa, owing to
a dependency of the standard SVV operator on a local Péclet number. First a modification
is proposed to the traditional SVV scaling that enforces a globally constant Péclet number
so as to avoid the previous issues. In addition, a new SVV kernel function is suggested
and shown to provide a more regular behaviour for the eigencurves along with a
consistent increase in resolution power for higher-order discretizations, as measured by
the extent of the wavenumber range where numerical errors are negligible. The dissipation
characteristics of CG with the SVV modifications suggested are then verified to be broadly
equivalent to those obtained through upwinding in the discontinuous Galerkin (DG)
scheme. Nevertheless, for the kernel function proposed, the full upwind DG scheme is
found to have a slightly higher resolution power for the same dissipation levels. These
results show that improved CG-SVV characteristics can be pursued via different kernel
functions with the aid of optimization algorithms.
Galerkin (CG) formulation in one dimension. First, numerical dispersion and diffusion
curves are obtained for the advection–diffusion problem and the role of multiple
eigencurves peculiar to spectral/hp methods is discussed. From the eigencurves’ behaviour,
we observe that CG might feature potentially undesirable non-smooth dispersion/diffusion
characteristics for under-resolved simulations of problems strongly dominated by either
convection or diffusion. Subsequently, the linear advection equation augmented with
spectral vanishing viscosity (SVV) is analysed. Dispersion and diffusion characteristics of
CG with SVV-based stabilization are verified to display similar non-smooth features in
flow regions where convection is much stronger than dissipation or vice-versa, owing to
a dependency of the standard SVV operator on a local Péclet number. First a modification
is proposed to the traditional SVV scaling that enforces a globally constant Péclet number
so as to avoid the previous issues. In addition, a new SVV kernel function is suggested
and shown to provide a more regular behaviour for the eigencurves along with a
consistent increase in resolution power for higher-order discretizations, as measured by
the extent of the wavenumber range where numerical errors are negligible. The dissipation
characteristics of CG with the SVV modifications suggested are then verified to be broadly
equivalent to those obtained through upwinding in the discontinuous Galerkin (DG)
scheme. Nevertheless, for the kernel function proposed, the full upwind DG scheme is
found to have a slightly higher resolution power for the same dissipation levels. These
results show that improved CG-SVV characteristics can be pursued via different kernel
functions with the aid of optimization algorithms.
Date Issued
2015-12-11
Date Acceptance
2015-12-05
Citation
Journal of Computational Physics, 2015, 307, pp.401-422
ISSN
1090-2716
Publisher
Elsevier
Start Page
401
End Page
422
Journal / Book Title
Journal of Computational Physics
Volume
307
Copyright Statement
© 2015 The Authors. Published by Elsevier Inc. This is an open access article under the CC
BY license (http://creativecommons.org/licenses/by/4.0/)
BY license (http://creativecommons.org/licenses/by/4.0/)
License URL
Sponsor
Engineering & Physical Science Research Council (E
Grant Number
EP/L000407/1
Subjects
Applied Mathematics
01 Mathematical Sciences
02 Physical Sciences
09 Engineering
Publication Status
Published