From h to p Efficiently: Selecting the Optimal Spectral/hp Discretisation in Three Dimensions

There is a growing interest in high-order finite and spectral/hp element methods using continuous and discontinuous Galerkin formulations. In this paper we investigate the effect of h- and p-type refinement on the relationship between runtime performance and solution accuracy. The broad spectrum of possible domain discretisations makes establishing a performance-optimal selection non-trivial. Through comparing the runtime of different implementations for evaluating operators over the space of discretisations with a desired solution tolerance, we demonstrate how the optimal discretisation and operator implementation may be selected for a specified problem. Furthermore, this demonstrates the need for codes to support both low- and high-order discretisa- tions.

[1]  D. Gottlieb,et al.  Numerical analysis of spectral methods : theory and applications , 1977 .

[2]  O. C. Zienkiewicz,et al.  The Finite Element Method: Its Basis and Fundamentals , 2005 .

[3]  George Em Karniadakis,et al.  TetrahedralhpFinite Elements , 1996 .

[4]  Spencer J. Sherwin,et al.  Hierarchical hp finite elements in hybrid domains , 1997 .

[5]  Long Chen.  FINITE ELEMENT METHOD , 2013 .

[6]  Sherwin,et al.  Tetrahedral hp Finite Elements : Algorithms and Flow Simulations , 1996 .

[7]  S. Orszag.  Spectral methods for problems in complex geometries , 1980 .

[8]  Moshe Dubiner.  Spectral methods on triangles and other domains , 1991 .

[9]  Jean-François Remacle,et al.  High-order discontinuous Galerkin schemes on general 2D manifolds applied to the shallow water equations , 2009, J. Comput. Phys..

[10]  Barry F. Smith,et al.  Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations , 1996 .

[11]  U. Lee.  Spectral Element Method in Structural Dynamics , 2009 .

[12]  G. Karniadakis,et al.  Spectral/hp Element Methods for Computational Fluid Dynamics , 2005 .

[13]  Robert Michael Kirby,et al.  From h to p efficiently: Implementing finite and spectral/hp element methods to achieve optimal performance for low- and high-order discretisations , 2010, J. Comput. Phys..

[14]  Steven A. Orszag,et al.  Spectral Methods for Problems in Complex Geometrics , 1979 .

[15]  J. Z. Zhu,et al.  The finite element method , 1977 .

[16]  A. Patera.  A spectral element method for fluid dynamics: Laminar flow in a channel expansion , 1984 .

[17]  S. Sherwin,et al.  From h to p efficiently: Strategy selection for operator evaluation on hexahedral and tetrahedral elements , 2011 .

[18]  J. Hesthaven,et al.  Nodal high-order methods on unstructured grids , 2002 .

[19]  Joel Ferziger,et al.  Higher Order Methods for Incompressible Fluid Flow: by Deville, Fischer and Mund, Cambridge University Press, 499 pp. , 2003 .

[20]  P. Fischer,et al.  High-Order Methods for Incompressible Fluid Flow , 2002 .

[21]  Spencer J. Sherwin,et al.  From h to p Efficiently: Implementing finite and spectral/hp element discretisations to achieve optimal performance at low and high order approximations. , 2009 .