An Analysis of Stability of the Flux Reconstruction Formulation on Quadrilateral Elements for the Linear Advection–Diffusion Equation

An Analysis of Stability of the Flux Reconstruction Formulation on Quadrilateral Elements for the... The Flux Reconstruction (FR) approach to high-order methods is a flexible and robust framework that has proven to be a promising alternative to the traditional Discontinuous Galerkin (DG) schemes on parallel architectures like Graphical Processing Units (GPUs) since it pairs exceptionally well with explicit time-stepping methods. The FR formulation was originally proposed by Huynh (AIAA Pap 2007-4079:1–42, 2007). Vincent et al. (J Sci Comput 47(1):50–72, 2011) later developed a single parameter family of correction functions which provide energy stable schemes under this formulation in 1D. These schemes, known as Vincent–Castonguay–Jameson–Huynh (VCJH) schemes, offer control over properties like stability, dispersion and dissipation through the variation of the VCJH parameter. Classical schemes like nodal-DG and Spectral Difference (SD) can also be recovered under this formulation. Following the development of the FR approach in 1D, Castonguay et al. (J Sci Comput 51(1):224–256, 2012) and Williams et al. (J Comput Phys 250:53–76, 2013) and Williams and Jameson (J Sci Comput 59(3):721–759, 2014) developed correction functions that give rise to energy stable FR formulations for triangles and tetrahedra. For the case of tensor product elements like quadrilaterals and hexahedra however, a simple extension of the 1D FR approach utilizing the 1D VCJH correction functions was possible and has been adopted by several authors (Castonguay in High-order energy stable flux reconstruction schemes for fluid flow simulations on unstructured grids, 2012; Witherden et al. in Comput Fluids 120:173–186, 2015; Comput Phys Commun 185(11):3028–3040, 2014). But whether such an extension of the 1D approach to tensor product elements is stable remained an open question. A direct extension of the 1D stability analysis fails due to certain key difficulties which necessitate the formulation of a norm different from the one utilized for stability analysis in 1D and on simplex elements. We have recently overcome these issues and shown that the VCJH schemes are stable for the linear advection equation on Cartesian meshes for any non-negative value of the VCJH parameter (Sheshadri and Jameson in J Sci Comput 67(2):769–790, 2016; J Sci Comput 67(2):791–794, 2016). In this paper, we have extended the stability analysis to the advection–diffusion equation, demonstrating that the tensor product FR formulation is stable on Cartesian meshes for the advection–diffusion case as well. The analysis in this paper also provides additional insights into the dependence on the VCJH parameter of the diffusion and stability characteristics of these schemes. Several numerical experiments that support the theoretical results are included. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Journal of Scientific Computing Springer Journals

An Analysis of Stability of the Flux Reconstruction Formulation on Quadrilateral Elements for the Linear Advection–Diffusion Equation

Loading next page...
 
/lp/springer_journal/an-analysis-of-stability-of-the-flux-reconstruction-formulation-on-1nmPab9KMK
Publisher
Springer US
Copyright
Copyright © 2017 by Springer Science+Business Media, LLC
Subject
Mathematics; Algorithms; Computational Mathematics and Numerical Analysis; Mathematical and Computational Engineering; Theoretical, Mathematical and Computational Physics
ISSN
0885-7474
eISSN
1573-7691
D.O.I.
10.1007/s10915-017-0513-9
Publisher site
See Article on Publisher Site

Abstract

The Flux Reconstruction (FR) approach to high-order methods is a flexible and robust framework that has proven to be a promising alternative to the traditional Discontinuous Galerkin (DG) schemes on parallel architectures like Graphical Processing Units (GPUs) since it pairs exceptionally well with explicit time-stepping methods. The FR formulation was originally proposed by Huynh (AIAA Pap 2007-4079:1–42, 2007). Vincent et al. (J Sci Comput 47(1):50–72, 2011) later developed a single parameter family of correction functions which provide energy stable schemes under this formulation in 1D. These schemes, known as Vincent–Castonguay–Jameson–Huynh (VCJH) schemes, offer control over properties like stability, dispersion and dissipation through the variation of the VCJH parameter. Classical schemes like nodal-DG and Spectral Difference (SD) can also be recovered under this formulation. Following the development of the FR approach in 1D, Castonguay et al. (J Sci Comput 51(1):224–256, 2012) and Williams et al. (J Comput Phys 250:53–76, 2013) and Williams and Jameson (J Sci Comput 59(3):721–759, 2014) developed correction functions that give rise to energy stable FR formulations for triangles and tetrahedra. For the case of tensor product elements like quadrilaterals and hexahedra however, a simple extension of the 1D FR approach utilizing the 1D VCJH correction functions was possible and has been adopted by several authors (Castonguay in High-order energy stable flux reconstruction schemes for fluid flow simulations on unstructured grids, 2012; Witherden et al. in Comput Fluids 120:173–186, 2015; Comput Phys Commun 185(11):3028–3040, 2014). But whether such an extension of the 1D approach to tensor product elements is stable remained an open question. A direct extension of the 1D stability analysis fails due to certain key difficulties which necessitate the formulation of a norm different from the one utilized for stability analysis in 1D and on simplex elements. We have recently overcome these issues and shown that the VCJH schemes are stable for the linear advection equation on Cartesian meshes for any non-negative value of the VCJH parameter (Sheshadri and Jameson in J Sci Comput 67(2):769–790, 2016; J Sci Comput 67(2):791–794, 2016). In this paper, we have extended the stability analysis to the advection–diffusion equation, demonstrating that the tensor product FR formulation is stable on Cartesian meshes for the advection–diffusion case as well. The analysis in this paper also provides additional insights into the dependence on the VCJH parameter of the diffusion and stability characteristics of these schemes. Several numerical experiments that support the theoretical results are included.

Journal

Journal of Scientific ComputingSpringer Journals

Published: Sep 25, 2017

References

You’re reading a free preview. Subscribe to read the entire article.


DeepDyve is your
personal research library

It’s your single place to instantly
discover and read the research
that matters to you.

Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.

All for just $49/month

Explore the DeepDyve Library

Search

Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly

Organize

Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.

Access

Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.

Your journals are on DeepDyve

Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.

All the latest content is available, no embargo periods.

See the journals in your area

DeepDyve

Freelancer

DeepDyve

Pro

Price

FREE

$49/month
$360/year

Save searches from
Google Scholar,
PubMed

Create lists to
organize your research

Export lists, citations

Read DeepDyve articles

Abstract access only

Unlimited access to over
18 million full-text articles

Print

20 pages / month

PDF Discount

20% off