×

zbMATH — the first resource for mathematics

High-order triangle-based discontinuous Galerkin methods for hyperbolic equations on a rotating sphere. (English) Zbl 1089.65096
Summary: High-order triangle-based discontinuous Galerkin (DG) methods for hyperbolic equations on a rotating sphere are presented. The DG method can be characterized as the fusion of finite elements with finite volumes. This DG formulation uses high-order Lagrange polynomials on the triangle using nodal sets up to 15th order. The finite element-type area integrals are evaluated using order 2\(N\) Gauss cubature rules. This leads to a full mass matrix which, unlike for continuous Galerkin (CG) methods such as the spectral element (SE) method presented by F. X. Giraldo and T. Warburton [ibid. 207, No. 1, 129–150 (2005; Zbl 1177.86002)], is small, local and efficient to invert.
Two types of finite volume-type flux integrals are studied: a set based on Gauss-Lobatto quadrature points (order 2\(N-1)\) and a set based on Gauss quadrature points (order 2\(N\)). Furthermore, we explore conservation and advection forms as well as strong and weak forms. Seven test cases are used to compare the different methods including some with scale contractions and shock waves. All three strong forms performed extremely well with the strong conservation form with \(2N\) integration being the most accurate of the four DG methods studied. The strong advection form with \(2N\) integration performed extremely well even for flows with shock waves. The strong conservation form with 2\(N - 1\) integration yielded results almost as good as those with 2\(N\) while being less expensive. All the DG methods performed better than the SE method for almost all the test cases, especially for those with strong discontinuities. Finally, the DG methods required less computing time than the SE method due to the local nature of the mass matrix.

MSC:
65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
35L65 Hyperbolic conservation laws
76B15 Water waves, gravity waves; dispersion and scattering, nonlinear interaction
76M12 Finite volume methods applied to problems in fluid mechanics
76M10 Finite element methods applied to problems in fluid mechanics
Software:
chammp; HE-E1GODF
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] Ata, R.; Soulaimani, A., A stabilized SPH method for inviscid shallow water flows, International journal for numerical methods in fluids, 47, 139-159, (2005) · Zbl 1065.76162
[2] Cools, R.; Rabinowitz, P., Monomial cubature rules Since stroud: a compilation, Journal of computational and applied mathematics, 48, 309-326, (1993) · Zbl 0799.65027
[3] Cools, R., Monomial cubature rules Since stroud: a compilation - part 2, Journal of computational and applied mathematics, 112, 21-27, (1999) · Zbl 0954.65021
[4] Diamessis, P.J.; Domaradzki, J.A.; Hesthaven, J.S., A spectral multidomain penalty method model for the simulation of high Reynolds number localized incompressible stratified turbulence, Journal of computational physics, 202, 298-322, (2005) · Zbl 1061.76054
[5] Dubiner, M., Spectral methods on triangles and other domains, Journal of scientific computing, 6, 345-390, (1991) · Zbl 0742.76059
[6] Dupont, F.; Lin, C.A., The adaptive spectral element method and comparisons with more traditional formulations for Ocean modeling, Journal of atmospheric and oceanic technology, 21, 135-147, (2004)
[7] Eskilsson, C.; Sherwin, S.J., A triangular spectral/hp discontinuous Galerkin method for modelling 2D shallow water equations, International journal for numerical methods in fluids, 45, 605-623, (2004) · Zbl 1085.76544
[8] Friedrich, O., Weighted essentially non-oscillatory schemes for the interpolation of Mean values on unstructured grids, Journal of computational physics, 144, 194-212, (1998) · Zbl 1392.76048
[9] F.X. Giraldo, A space marching adaptive remeshing technique applied to the 3D Euler equations for supersonic flow, PhD thesis, University of Virginia, 1995.
[10] F.X. Giraldo, A finite volume high resolution 2D Euler solver with adaptive grid generation on high performance computers, in: M. Morandi-Cecchi, K. Morgan, J. Periaux, B.A. Schrefler, O.C. Zienkiewicz (Eds.), Proceedings of the 9th International Conference on Finite Element in Fluids, Venice, Italy, 1995.
[11] Giraldo, F.X., A spectral element shallow water model on spherical geodesic grids, International journal for numerical methods in fluids, 35, 869-901, (2001) · Zbl 1030.76045
[12] Giraldo, F.X.; Hesthaven, J.S.; Warburton, T., Nodal high-order discontinuous Galerkin methods for the spherical shallow water equations, Journal of computational physics, 181, 499-525, (2002) · Zbl 1178.76268
[13] Giraldo, F.X.; Perot, J.B.; Fischer, P.F., A spectral element semi-Lagrangian (SESL) method for the spherical shallow water equations, Journal of computational physics, 190, 623-650, (2003) · Zbl 1076.76058
[14] Giraldo, F.X., Semi-implicit time-integrators for a scalable spectral element atmospheric model, Quarterly journal of the royal meteorological society, 131, 2431-2454, (2005)
[15] Giraldo, F.X.; Warburton, T., A nodal triangle-based spectral element method for the shallow water equations on the sphere, Journal of computational physics, 207, 129-150, (2005) · Zbl 1177.86002
[16] Hesthaven, J.S., From electrostatics to almost optimal nodal sets for polynomial interpolation in a simplex, SIAM journal on numerical analysis, 35, 655-676, (1998) · Zbl 0933.41004
[17] Hu, C.Q.; Shu, C.W., Weighted essentially non-oscillatory schemes on triangular meshes, Journal of computational physics, 150, 97-127, (1999) · Zbl 0926.65090
[18] Hundsdorer, W.; Jaffré, J., Implicit-explicit time stepping with spatial discontinuous finite elements, Applied numerical mathematics, 45, 231-254, (2003) · Zbl 1029.65115
[19] Karniadakis, G.E.; Israeli, M.; Orszag, S.A., High-order splitting methods for the incompressible Navier-Stokes equations, Journal of computational physics, 97, 414-443, (1991) · Zbl 0738.76050
[20] Koornwinder, T., Two-variable analogues of the classical orthogonal polynomials, () · Zbl 0326.33002
[21] Lyness, J.; Cools, R., A survey of numerical cubature over triangles, Applied mathematics, 48, 127-150, (1994) · Zbl 0820.41026
[22] Proriol, J., Sur une famille de polynomes á deux variables orthogonaux dans un triangle, C.R. acadamic science Paris, 257, 2459, (1957) · Zbl 0080.05204
[23] Rossmanith, J.A.; Bale, D.S.; LeVeque, R.J., A wave propagation algorithm for hyperbolic systems on curved manifolds, Journal of computational physics, 199, 631-662, (2004) · Zbl 1126.76350
[24] Sanchez, J.; Net, M.; García-Archilla, B.; Simó, C., Newton-Krylov continuation of periodic orbits for Navier-Stokes flows, Journal of computational physics, 201, 13-33, (2004) · Zbl 1153.76345
[25] Shu, C.W., High-order finite difference and finite volume WENO schemes and discontinuous Galerkin methods for CFD, International journal of computational fluid dynamics, 17, 107-118, (2003) · Zbl 1034.76044
[26] Stroud, A.H., Approximate calculation of multiple integrals, (1971), Prentice-Hall Publishing New Jersey · Zbl 0379.65013
[27] Taylor, M.A.; Wingate, B.A.; Vincent, R.E., An algorithm for computing Fekete points in the triangle, SIAM journal on numerical analysis, 38, 1707-1720, (2000) · Zbl 0986.65017
[28] Toro, E.F., Riemann solvers and numerical methods for fluid dynamics, A practical introduction, (1999), Springer New York · Zbl 0923.76004
[29] Wang, Z.J.; Zhang, L.; Liu, Y., Spectral (finite) volume method for conservation laws on unstructured grids IV: extension to two-dimensional systems, Journal of computational physics, 194, 716-741, (2004) · Zbl 1039.65072
[30] Williamson, D.L.; Drake, J.B.; Hack, J.J.; Jakob, R.; Swarztrauber, P.N., A standard test set for numerical approximations to the shallow water equations in spherical geometry, Journal of computational physics, 102, 211-224, (1992) · Zbl 0756.76060
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.