Energy stable flux reconstruction schemes for advection-diffusion problems on triangles.

*(English)*Zbl 1349.65528Summary: The Flux Reconstruction (FR) approach unifies several well-known high-order schemes for unstructured grids, including a collocation-based nodal discontinuous Galerkin (DG) method and all types of Spectral Difference (SD) methods, at least for linear problems. The FR approach also allows for the formulation of new families of schemes. Of particular interest are the energy stable FR schemes, also referred to as the Vincent-Castonguay-Jameson-Huynh (VCJH) schemes, which are an infinite family of high-order schemes parameterized by a single scalar. VCJH schemes are of practical importance because they provide a stable formulation on triangular elements which are often required for numerical simulations over complex geometries. In particular, VCJH schemes are provably stable for linear advection problems on triangles, and include the collocation-based nodal DG scheme on triangles as a special case. Furthermore, certain VCJH schemes have Courant-Friedrichs-Lewy (CFL) limits which are approximately twice those of the collocation-based nodal DG scheme. Thus far, these schemes have been analyzed primarily in the context of pure advection problems on triangles. For the first time, this paper constructs VCJH schemes for advection-diffusion problems on triangles, and proves the stability of these schemes for linear advection-diffusion problems for all orders of accuracy. In addition, this paper uses numerical experiments on triangular grids to verify the stability and accuracy of VCJH schemes for linear advection-diffusion problems and the nonlinear Navier-Stokes equations.

##### MSC:

65M70 | Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs |

65M12 | Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs |

76F10 | Shear flows and turbulence |

##### Keywords:

high-order; discontinuous Galerkin; spectral difference; flux reconstruction; triangles; advection-diffusion
PDF
BibTeX
XML
Cite

\textit{D. M. Williams} et al., J. Comput. Phys. 250, 53--76 (2013; Zbl 1349.65528)

Full Text:
DOI

**OpenURL**

##### References:

[1] | Cockburn, B.; Hou, S.; Shu, C., The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws IV: the multidimensional case, Math. Comput., 54, 545-581, (1990) · Zbl 0695.65066 |

[2] | Bassi, F.; Rebay, S., High-order accurate discontinuous finite element solution of the 2D Euler equations, J. Comput. Phys., 138, 251-285, (1997) · Zbl 0902.76056 |

[3] | Bassi, F.; Rebay, S., A high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier-Stokes equations, J. Comput. Phys., 131, 267-279, (1997) · Zbl 0871.76040 |

[4] | Kopriva, D. A.; Kolias, J. H., A conservative staggered-grid Chebyshev multidomain method for compressible flows, J. Comput. Phys., 125, 244-261, (1996) · Zbl 0847.76069 |

[5] | Liu, Y.; Vinokur, M.; Wang, Z. J., Spectral difference method for unstructured grids I: basic formulation, J. Comput. Phys., 216, 780-801, (2006) · Zbl 1097.65089 |

[6] | H. Huynh, A flux reconstruction approach to high-order schemes including discontinuous Galerkin methods, AIAA P., in: 18th AIAA Computational Fluid Dynamics Conference, Miami, FL, Jun 25-28, 2007. |

[7] | H. Huynh, A reconstruction approach to high-order schemes including discontinuous Galerkin for diffusion, AIAA P., in: 47th AIAA Aerospace Sciences Meeting, Orlando, FL, Jan 5-8, 2009. |

[8] | H. Huynh, High-order methods including discontinuous Galerkin by reconstructions on triangular meshes, AIAA P., in: 49th AIAA Aerospace Sciences Meeting, Orlando, FL, Jan 4-7, 2011. |

[9] | H. Gao, Z.J. Wang, A high-order lifting collocation penalty formulation for the Navier-Stokes equations on 2D mixed grids, AIAA P., in: 19th AIAA Computational Fluid Dynamics Conference, San Antonio, TX, June 22-25, 2009. |

[10] | Wang, Z.; Gao, H., A unifying lifting collocation penalty formulation including the discontinuous Galerkin, spectral volume/difference methods for conservation laws on mixed grids, J. Comput. Phys., 228, 8161-8186, (2009) · Zbl 1173.65343 |

[11] | Yu, M.; Wang, Z., On the connection between the correction and weighting functions in the correction procedure via reconstruction method, J. Scient. Comput., (2012), DOI: 10.1007/s10915-012-9618-3 |

[12] | Haga, T.; Gao, H.; Wang, Z., A high-order unifying discontinuous formulation for the Navier-Stokes equations on 3D mixed grids, Math. Model. Nat. Phenom., 6, 28-56, (2011) · Zbl 1239.76044 |

[13] | T. Haga, H. Gao, Z. Wang, A high-order unifying discontinuous formulation for 3D mixed grids, AIAA P., in: 48th AIAA Aerospace Sciences Meeting, Orlando, FL, Jan 4-7, 2010. · Zbl 1239.76044 |

[14] | Nordström, J.; Forsberg, K.; Adamsson, C.; Eliasson, P., Finite volume methods, unstructured meshes and strict stability for hyperbolic problems, Appl. Numer. Math., 45, 453-473, (2003) · Zbl 1019.65066 |

[15] | Ham, F., Improved scalar transport for unstructured finite volume methods using simplex superposition, annual research briefs, Center for Turbulence Research, NASA-AMES, 347-358, (2008) |

[16] | Kreiss, H.; Scherer, G., Finite element and finite difference methods for hyperbolic partial differential equations, Math. Aspect. Finite Elem. Part. Differ. Equat., 195-212, (1974), Publ. Math. Res. Center Univ. Wisconsi · Zbl 0355.65085 |

[17] | H. Kreiss, G. Scherer, On the existence of energy estimates for difference approximations for hyperbolic systems, Tech. Rep., Uppsala University, Division of Scientific Computing, 1977. |

[18] | Carpenter, M.; Gottlieb, D.; Abarbanel, S., Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes, J. Comput. Phys., 111, 220-236, (1994) · Zbl 0832.65098 |

[19] | Carpenter, M.; Nordström, J.; Gottlieb, D., A stable and conservative interface treatment of arbitrary spatial accuracy, J. Comput. Phys., 148, 341-365, (1999) · Zbl 0921.65059 |

[20] | Jameson, A., A proof of the stability of the spectral difference method for all orders of accuracy, J. Scient. Comput., 45, 348-358, (2010) · Zbl 1203.65198 |

[21] | Vincent, P. E.; Castonguay, P.; Jameson, A., A new class of high-order energy stable flux reconstruction schemes, J. Scient. Comput., 47, 50-72, (2011) · Zbl 1433.76094 |

[22] | D.M. Williams, P. Castonguay, P.E. Vincent, A. Jameson, An extension of energy stable flux reconstruction to unsteady, non-linear, viscous problems on mixed grids, AIAA P., in: 20th AIAA Computational Fluid Dynamics Conference, Honolulu, Hawaii, June 27-30, 2011. |

[23] | P. Castonguay, High-Order Energy Stable Flux Reconstruction Schemes for Fluid Flow Simulations on Unstructured Grids, Ph.D. thesis, Stanford University, 2012. |

[24] | P. Castonguay, D.M. Williams, P.E. Vincent, A. Jameson, Energy stable flux reconstruction for advection-diffusion problems, Comput. Meth. Appl. Mech. Eng., 2012. · Zbl 1286.65119 |

[25] | Castonguay, P.; Vincent, P. E.; Jameson, A., A new class of high-order energy stable flux reconstruction schemes for conservation laws on triangular grids, J. Scient. Comput., (2011), http://dx.doi.org/10.1007/s10915-012-9618-3 |

[26] | Viviand, H., Conservative forms of gas dynamic equations, La Recherche Aerospatiale, 1, 65-66, (1974) |

[27] | Vinokur, M., Conservation equations of gasdynamics in curvilinear coordinate systems, J. Comput. Phys., 14, 105-125, (1974) · Zbl 0277.76061 |

[28] | J. Hesthaven, T. Warburton, Nodal Discontinuous Galerkin methods: Algorithms, Analysis, and Applications, Springer Verlag, 2007. · Zbl 1134.65068 |

[29] | Cockburn, B.; Shu, C., The local discontinuous Galerkin method for time-dependent convection-diffusion systems, SIAM J. Numer. Anal., 35, 2440-2463, (1998) · Zbl 0927.65118 |

[30] | Peraire, J.; Persson, P., The compact discontinuous Galerkin (CDG) method for elliptic problems, SIAM J. Scient. Comput., 30, 1806-1824, (2009) · Zbl 1167.65436 |

[31] | Arnold, D., An interior penalty finite element method with discontinuous elements, SIAM J. Numer. Anal., 19, 742-760, (1982) · Zbl 0482.65060 |

[32] | F. Bassi, S. Rebay, G. Mariotti, S. Pedinotti, M. Savini, A high-order accurate discontinuous finite element method for inviscid and viscous turbomachinery flows, in: R. Decuypere, G. Dibelius (Eds.), 2nd European Conference on Turbomachinery Fluid Dynamics and Thermodynamics (Antwerpen, Belgium), Technologisch Instituut, 1997, pp. 99-108. |

[33] | Roe, P. L., Approximate Riemann solvers, parameter vectors and difference schemes, J. Comput. Phys., 43, 357-372, (1981) · Zbl 0474.65066 |

[34] | Rusanov, V. V., Calculation of interaction of non-steady shock waves with obstacles, J. Comput. Math. Phys. USSR, 1, 261-279, (1961) |

[35] | Raviart, P.; Thomas, J., A mixed finite element method for second-order elliptic problems, Math. Aspect. Finite Elem. Meth., 292-315, (1977) |

[36] | Dubiner, M., Spectral methods on triangles and other domains, J. Scient. Comput., 6, 345-390, (1991) · Zbl 0742.76059 |

[37] | M.H. Carpenter, C. Kennedy, Fourth-order 2N-storage Runge-Kutta schemes, Technical Report TM 109112, NASA, NASA Langley Research Center, 1994. |

[38] | P. Castonguay, P.E. Vincent, A. Jameson, Application of energy stable flux reconstruction schemes for the Euler equations, AIAA P., in: 49th AIAA Aerospace Sciences Meeting, Orlando, FL, Jan 4-7, 2011. · Zbl 1433.76094 |

[39] | UIUC Applied Aerodynamics Group. UIUC airfoil coordinates database, http://www.ae.illinois.edu/m-selig/ads/coord_database.html, Retrieved April 2013. |

[40] | A. Uranga, P.-O. Persson, M. Drela, J. Peraire, Implicit large eddy simulation of transition to turbulence at low Reynolds numbers using a discontinuous Galerkin method, Int. J. Numer. Meth. Eng., 87 (1-5) (2011) 232-261. · Zbl 1242.76085 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.