zbMATH — the first resource for mathematics

From \(h\) to \(p\) efficiently: implementing finite and spectral/hp element methods to achieve optimal performance for low- and high-order discretisations. (English) Zbl 1194.65138
Summary: The spectral/\(hp\) element method can be considered as bridging the gap between the – traditionally low-order – finite element method on one side and spectral methods on the other side. Consequently, a major challenge which arises in implementing the spectral/\(hp\) element methods is to design algorithms that perform efficiently for both low- and high-order spectral/\(hp\) discretisations, as well as discretisations in the intermediate regime.
In this paper, we explain how the judicious use of different implementation strategies can be employed to achieve high efficiency across a wide range of polynomial orders. Furthermore, based upon this efficient implementation, we analyse which spectral/\(hp\) discretisation (which specific combination of mesh-size \(h\) and polynomial order \(P\)) minimises the computational cost to solve an elliptic problem up to a predefined level of accuracy. We investigate this question for a set of both smooth and non-smooth problems.

65N30 Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs
65N35 Spectral, collocation and related methods for boundary value problems involving PDEs
35J05 Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation
Full Text: DOI
[1] Sherwin, S.J.; Karniadakis, G.E., Tetrahedral hp finite elements: algorithms and flow simulations, J. comput. phys., 124, 14-45, (1996) · Zbl 0847.76038
[2] Hesthaven, J.; Warburton, T.C.E., Nodal high-order methods on unstructured grids: I. time-domain solution of maxwell’s equations, J. comput. phys., 181, 186-221, (2002) · Zbl 1014.78016
[3] Bernard, P.-E.; Remacle, J.-F.; Combien, R.; Legat, V.; Hillewaert, K., High-order discontinuous Galerkin schemes on general 2D manifolds applied to the shallow water equations, J. comput. phys., 228, 6514-6535, (2009) · Zbl 1261.76026
[4] Zienkiewicz, O.C.; Taylor, R.L., The finite element method, (1989), McGraw-Hill New York · Zbl 0991.74002
[5] Hughes, T.J.R., The finite element method, (1987), Prentice-Hall New Jersey
[6] Szabó, B.; Babuška, I., Finite element analysis, (1991), Wiley New York
[7] Karniadakis, G.E.; Sherwin, S.J., Spectral/hp element methods for computational fluid dynamics. numerical mathematics and scientific computation, (2005), Oxford University Press Oxford
[8] Deville, M.O.; Fischer, P.F.; Mund, E.H., High-order methods for incompressible fluid flow. Cambridge monographs on applied and computational mathematics, (2002), Cambridge University Press Cambridge · Zbl 1007.76001
[9] Hesthaven, J.S.; Warburton, T., Nodal discontinuous Galerkin methods: algorithms, analysis, and applications, Springer texts in applied mathematics, vol. 54, (2008), Springer-Verlag New York · Zbl 1134.65068
[10] D. Gottlieb, S.A. Orszag, Numerical Analysis of Spectral Methods: Theory and Applications, CBMS-NSF, Society for Industrial and Applied Mathematics, Philadelphia, 1977. · Zbl 0412.65058
[11] Remacle, J.F.; Flaherty, J.E.; Shepard, M.S., An adaptive discontinuous Galerkin technique with an orthogonal basis applied to compressible flow problems, SIAM rev., 45, 55-73, (2003)
[12] Hughes, T.J.R.; Cottrell, J.A.; Bazilevs, Y., Isogeometric analysis: cad, finite elements, nurbs, exact geometry and mesh refinement, Comput. meth. appl. mech. eng., 194, 39-41, 4135-4195, (2005) · Zbl 1151.74419
[13] Dubiner, Moshe, Spectral methods on triangles and other domains, J. sci. comput., 6, 4, 345-390, (1991) · Zbl 0742.76059
[14] Bagheri, Babak; Ridgway Scott, L.; Zhang, Shangyou, Implementing and using high-order finite element methods, Finite elem. anal. des., 16, 3-4, 175-189, (1994), 186260 · Zbl 0812.76042
[15] R.C. Kirby, M. Knepley, L.R. Scott, Evaluation of the action of finite element operators, Technical Report TR-2004-07, University of Chicago, 2004.
[16] P.E.J. Vos, S.J. Sherwin, R.M. Kirby, From h to p efficiently, Technical Report UUCS-09-007, University of Utah, 2009.
[17] Orszag, S.A., Spectral methods for problems in complex geometries, J. comput. phys., 37, 1, 70-92, (1980) · Zbl 0476.65078
[18] E.M. Rønquist, Optimal spectral element methods for unsteady three-dimensional incompressible Navier-Stokes equations, Ph.D. Thesis, Massachusetts Institute of Technology, Cambridge, MA, 1988.
[19] Fischer, P.; Gottlieb, D., On the optimal number of subdomains for hyperbolic problems on parallel computers, Int. J. supercomputer appl. high perform. comput., 11, 65-76, (1996)
[20] Hesthaven, J.S., A stable penalty method for the compressible navier – stokes equations: ii. one-dimensional domain decomposition schemes, SIAM J. sci. comput., 18, 3, 658-685, (1997) · Zbl 0882.76061
[21] Hesthaven, J.S., A stable penalty method for the compressible navier – stokes equations: iii. multidimensional domain decomposition schemes, SIAM J. sci. comput., 20, 1, 62-93, (1998) · Zbl 0957.76059
[22] Wasberg, C.E.; Gottlieb, D., Optimal decomposition of the domain in spectral methods for wave-like phenomena, SIAM J. sci. comput., 22, 2, 617-632, (2000) · Zbl 0967.65102
[23] Babuška, I.; Suri, M., The p and h – p versions of the finite element method, basic principles and properties, SIAM rev., 36, 4, 578-632, (1994) · Zbl 0813.65118
[24] Seshaiyer, P.; Suri, M., hp submeshing via non-conforming finite element methods, Comput. meth. appl. mech. eng., 189, 3, 1011-1030, (2000) · Zbl 0971.65101
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.