Abstract
Conventional wisdom suggests that high-order finite-difference methods are more efficient than high-order discontinuous spectral-element methods on smooth meshes, but less efficient as the mesh becomes increasingly distorted because of a significant loss of accuracy on such meshes. This paper investigates the influence of mesh distortion on the relative efficiency of different implementations of generalized summation-by-parts (GSBP) methods, with emphasis on comparing finite-difference and discontinuous spectral-element approaches. These include discretizations built using classical finite-difference SBP operators, with and without optimized boundary closures, as well as both Legendre-Gauss and Legendre-Gauss-Lobatto operators. The traditionally finite-difference operators are also applied as discontinuous spectral-element operators by selecting a fixed number of nodes per element and performing mesh refinement by increasing the number of elements rather than the number of mesh nodes. Using the linear convection equation and nonlinear Euler equations as models, solutions are obtained on meshes with different types and severity of distortion. Contrary to expectation, the results show that finite-difference implementations are no more sensitive to mesh distortion than discontinuous spectral-element implementations, maintaining their relative efficiency in most cases. The results also show that the operators of Mattsson et al. (J Comput Phys 264:91–111, 2014) with optimized boundary operators are often the most efficient for a given implementation strategy (finite-difference or discontinuous spectral-element). While their accuracy as finite-difference operators might be expected, their superior accuracy to LG and LGL nodal distributions when implemented as discontinuous spectral-element operators is not well known.
We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.
Data Availability
The authors declare that the data supporting the findings of this study are available within the paper and its supplementary files in graphical form. Raw data can be provided from the corresponding author on reasonable request.
Notes
We expect that three-dimensional results will exhibit similar trends.
For first derivative operators, the order and polynomial degree of the operator are the same.
The values presented here are rounded to 3 decimal places for the sake of space. Please refer to [14] for higher precision values of the operator’s coefficients.
Figures show the square root of the number of nodes; accounting for the simulations being in two dimensions. This is to enable reference convergence slopes to be added to the figures.
The term cell is used here to indicate the quadrilateral formed by 4 nodes in a mesh with indices (i, j), \((i+1,j)\), \((i+1,j+1)\), and \((i,j+1)\) - not a spectral-element.
References
Craig Penner, D.A.: Development and Investigation of Accurate High-Order Generalized Summation-by-Parts Discretizations for Computational Fluid Dynamics. University of Toronto, Toronto (2022)
Craig Penner, D.A., Zingg, D.W.: Accurate high-order tensor-product generalized summation-by-parts discretizations of hyperbolic conservation laws: General curved domains and functional superconvergence. J. Sci. Comput. 93(2), 36 (2022)
Diener, P., Dorband, E.N., Schnetter, E., Tiglio, M.: Optimized high-order derivative and dissipation operators satisfying summation by parts, and applications in three-dimensional multi-block evolutions. J. Sci. Comput. 32(1), 109–145 (2007)
Del Rey Fernández, D.C., Boom, P.D., Zingg, D.W.: A generalized framework for nodal first derivative summation-by-parts operators. J. Comput. Phys. 266(1), 214–239 (2014)
Del Rey Fernández, D.C., Hicken, J.E., Zingg, D.W.: Review of summation-by-parts operators with simultaneous approximation terms for the numerical solution of partial differential equations. Comput. Fluids 95(22), 171–196 (2014)
Del Rey Fernández, D.C., Hicken, J.E., Zingg, D.W.: Simultaneous approximation terms for multidimensional summation-by-parts operators. J. Sci. Comput. 75(1), 83–110 (2018)
Gassner, G.J.: A skew-symmetric discontinuous Galerkin spectral element discretization and its relation to SBP-SAT finite difference methods. SIAM J. Sci. Comput. 35(3), A1233–A1253 (2013)
Hesthaven, J.S., Warburton, T.: Nodal Discontinuous Galerkin Methods Algorithms, Analysis, and Applications. Springer, Cham (2008)
Hicken, J.E., Del Rey Fernández, D.C., Zingg, D.W.: Multi-dimensional summation-by-parts operators: general theory and application to simplex elements. SIAM J. Sci. Comput. 38(4), A1935–A1958 (2016)
Hicken, J.E., Zingg, D.W.: Parallel Newton–Krylov solver for the Euler equations discretized using simultaneous approximation terms. AIAA J. 46(11), 2773–2786 (2008)
Hindenlang, F.: Mesh Curving Techniques for High Order Parallel Simulations on Unstructured Meshes. University of Stuttgart, Stuttgart (2014)
Kopriva, D.A., Gassner, G.J.: An energy stable discontinuous Galerkin spectral element discretization for variable coefficient advection problems. SIAM J. Sci. Comput. 4(36), A2076–A2099 (2014)
Linders, V., Lundquist, T., Nordström, J.: On the order of accuracy of finite difference operators on diagonal norm based summation-by-parts form. SIAM J. Numer. Anal. 56(2), 1048–1063 (2018)
Mattsson, K., Almquist, M., Carpenter, M.H.: Optimal diagonal-norm SBP operators. J. Comput. Phys. 264(1), 91–111 (2014)
Mattsson, K., Almquist, M., van der Weide, E.: Boundary optimized diagonal-norm SBP operators. J. Comput. Phys. 374, 1261–1266 (2018)
Nordström, J.: Conservative finite difference formulations, variable coefficients, energy estimates and artificial dissipation. J. Sci. Comput. 29(3), 375–404 (2006)
Olsson, P., Oliger, J.: Energy and Maximum Norm Estimates for Nonlinear Conservation Laws, pp. 94–101. The Research Institute of Advanced Computer Science, Mountain View (1994)
Shu, C.W.: High-order finite difference and finite volume WENO schemes and discontinuous Galerkin methods for CFD. Int. J. Comput. Fluid Dyn. 17(2), 107–118 (2003)
Svärd, M., Nordström, J.: Review of summation-by-parts schemes for initial-boundary-value-problems. J. Comput. Phys. 268(1), 17–38 (2014)
Funding
The authors have not disclosed any funding.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors have no conflict of interest to declare that are relevant to the content of this article.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Boom, P.D., Del Rey Fernández, D.C. & Zingg, D.W. Relative Efficiency of Finite-Difference and Discontinuous Spectral-Element Summation-by-Parts Methods on Distorted Meshes. J Sci Comput 102, 28 (2025). https://doi.org/10.1007/s10915-024-02745-5
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10915-024-02745-5