Abstract
Spectral element methods (SEM), extensions of finite element methods (FEM), have emerged as significant techniques for solving partial differential equations in physics and engineering. SEM can potentially deliver superior accuracy due to the potential superconvergence in nodal solutions for well-shaped tensor-product elements. However, the accuracy of SEM often degrades in complex geometries due to geometric inaccuracies near curved boundaries and the loss of superconvergence with simplicial or non-tensor-product elements. To overcome the first issue, we propose using h- and p-geometric refinement, which refines the mesh near high-curvature regions and increases the degree of geometric basis functions, respectively. We show that when using mixed-element meshes with tensor-product elements in the interior of the domain, curvature-based geometric refinement near boundaries can improve the accuracy of the interior elements by reducing pollution errors and preserving the superconvergence in nodal solutions. To address the second issue, we introduce a post-processing technique using the adaptive extended stencil finite element method (AES-FEM) to recover the accuracy near the curved boundaries. The combination of curvature-based geometric refinement and accurate post-processing offers an effective and easier-to-implement alternative to methods reliant on exact geometries. We demonstrate our techniques by solving the convection–diffusion equation in 2D and show up to two orders of magnitude of improvement in the solution accuracy, even when the elements are poorly shaped near boundaries.
Similar content being viewed by others
Data availability
Data will be made available on request.
References
Babuška I, Miller A (1984) The post-processing approach in the finite element method—part 1: calculation of displacements, stresses and other higher derivatives of the displacements. In. J Numer Methods Eng 20(6):1085–1109
Bao H, Ding D, Chen R (2017) A hybrid spectral-element finite-difference time-domain method for electromagnetic simulation. IEEE Antennas Wirel Propag Lett 16:2244–2248
Bassi F, Rebay S (1997) High-order accurate discontinuous finite element solution of the 2D Euler equations. J Comput Phys 138(2):251–285
Boroomand B, Zienkiewicz O (1997) An improved REP recovery and the effectivity robustness test. Int J Numer Methods Eng 40(17):3247–3277
Boroomand B, Zienkiewicz OC (1997) Recovery by equilibrium in patches (REP). Int J Numer Methods Eng 40(1):137–164
Botti L, Di Pietro DA (2018) Assessment of hybrid high-order methods on curved meshes and comparison with discontinuous Galerkin methods. J Comput Phys 370:58–84
Bramble JH, Schatz AH (1977) Higher order local accuracy by averaging in the finite element method. Math Comp 31(137):94–111
Brenner S, Scott LR (2008) The mathematical theory of finite element methods, 3rd edn. Springer
Brenner SC, Carstensen C (2004) Finite element methods. Comput Mech 1:73–114
Chen C (2005) Orthogonality correction technique in superconvergence analysis. Int J Numer Anal Model 2005(1):31–42
Chen C, Hu S (2013) The highest order superconvergence for bi-k degree rectangular elements at nodes: a proof of 2k-conjecture. Math Comput 82(283):1337–1355
Chen CM (1981) Superconvergence of finite element solutions and their derivatives. Numer Math J Chin Univ 3(2):118–125
Chen CM (1999) Superconvergence for triangular finite elements. Sci China Ser A Math 42(9):917–924
Chen CM, Huang YQ (1995) High accuracy theory of finite elements (in Chinese). Hunan Science and Technique Press, Changsha
Cheung J, Perego M, Bochev P, Gunzburger M (2019) Optimally accurate higher-order finite element methods for polytopial approximations of domains with smooth boundaries. Math Comput 88(319):2187–2219
Ciarlet PG, Raviart P-A (1972) The combined effect of curved boundaries and numerical integration in isoparametric finite element methods. In: A.K. Aziz (ed) The mathematical foundations of the finite element method with applications to partial differential equations. Elsevier, pp 409–474
Conley R, Delaney TJ, Jiao X (2016) Overcoming element quality dependence of finite elements with adaptive extended stencil FEM (AES-FEM). Int J Numer Methods Eng 108(9):1054–1085
Conley R, Delaney TJ, Jiao X (2020) A hybrid method and unified analysis of generalized finite differences and Lagrange finite elements. J Comput Appl Math 376:112862
Dassi F, Mola A, Si H (2014) Curvature-adapted remeshing of CAD surfaces. Proc Eng 82:253–265
Davies R, Morgan K, Hassan O (2009) A high order hybrid finite element method applied to the solution of electromagnetic wave scattering problems in the time domain. Comput Mech 44:321–331
do Carmo M (1976) Differential geometry of curves and surfaces. Prentice-Hall
Dyedov V, Ray N, Einstein D, Jiao X, Tautges T (2015) AHF: array-based half-facet data structure for mixed-dimensional and non-manifold meshes. Eng Comput 31:389–404
Engvall L, Evans JA (2020) Mesh quality metrics for isogeometric Bernstein-Bézier discretizations. Comput Methods Appl Mech Eng 371:113305
Eslami MR (2014) Finite elements methods in mechanics. Springer
Finlayson BA (2013) The method of weighted residuals and variational principles. SIAM
Hughes TJ, Cottrell JA, Bazilevs Y (2005) Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Comput Methods Appl Mech Eng 194(39–41):4135–4195
Jiao X, Wang D, Zha H (2011) Simple and effective variational optimization of surface and volume triangulations. Eng Comput 27:81–94
Karniadakis G, Sherwin S (2005) Spectral/hp element methods for computational fluid dynamics. OUP Oxford
Khan D, Plopski A, Fujimoto Y, Kanbara M, Jabeen G, Zhang YJ, Zhang X, Kato H (2020) Surface remeshing: a systematic literature review of methods and research directions. IEEE Trans Vis Comput Graph 28(3):1680–1713
Komatitsch D, Martin R, Tromp J, Taylor MA, Wingate BA (2001) Wave propagation in 2-D elastic media using a spectral element method with triangles and quadrangles. J Comput Acoust 9(02):703–718
Křížek M, Neittaanmäki P (1987) On superconvergence techniques. Acta Appl Math 9:175–198
Ledger P, Morgan K, Hassan O, Weatherill N (2002) Arbitrary order edge elements for electromagnetic scattering simulations using hybrid meshes and a PML. Int J Numer Methods Eng 55(3):339–358
Lenoir M (1986) Optimal isoparametric finite elements and error estimates for domains involving curved boundaries. SIAM J Numer Anal 23(3):562–580
Li Y, Zhao X, Ray N, Jiao X (2021) Compact feature-aware hermite-style high-order surface reconstruction. Eng Comput 37(1):187–210
Lions P-L et al (1988) On the Schwarz alternating method. i. In: First International Symposium on domain decomposition methods for partial differential equations, volume 1, p. 42. Paris, France
Lloyd S (1982) Least squares quantization in PCM. IEEE Trans Inf Theory 28(2):129–137
Luo X, Shephard MS, Remacle J-F (2001) The influence of geometric approximation on the accuracy of high order methods. Rensselaer SCOREC report, p 1
Maday Y, Rønquist EM (1990) Optimal error analysis of spectral methods with emphasis on non-constant coefficients and deformed geometries. Comput Methods Appl Mech Eng 80(1–3):91–115
Mansouri S, Ebrahimnezhad H (2016) Segmentation-based semi-regular remeshing of 3D models using curvature-adapted subdivision surface fitting. J Vis 19:141–155
Mikhlin S (1951) On the schwarz algorithm. Dokl Akad Nauk SSSR 77(4):569–571
Moxey D, Green M, Sherwin S, Peiró J (2015) An isoparametric approach to high-order curvilinear boundary-layer meshing. Comput Methods Appl Mech Eng 283:636–650
Naga A, Zhang Z (2004) A Posteriori error estimates based on the polynomial preserving recovery. SIAM J Numer Anal 42(4):1780–1800
Patera AT (1984) A spectral element method for fluid dynamics: laminar flow in a channel expansion. J Comput Phys 54(3):468–488
Remacle J-F, Lambrechts J, Toulorge T, Johnen A, Geuzaine C (2014) Optimizing the geometrical accuracy of 2D curvilinear finite element meshes. In: 6th International Conference on advanced computational methods in engineering (ACOMEN 2014)
Ruiz-Gironés E, Sarrate J, Roca X (2016) Generation of curved high-order meshes with optimal quality and geometric accuracy. Proc Eng 163:315–327
Ruiz-Gironés E, Sarrate J, Roca X (2021) Measuring and improving the geometric accuracy of piece-wise polynomial boundary meshes. J Comput Phys 443:110500
Scott MA, Borden MJ, Verhoosel CV, Sederberg TW, Hughes TJ (2011) Isogeometric finite element data structures based on Bézier extraction of T-splines. Int J Numer Methods Eng 88(2):126–156
Sevilla R, Fernández-Méndez S, Huerta A (2008) NURBS-enhanced finite element method (NEFEM). Int J Numer Methods Eng 76(1):56–83
Sevilla R, Fernández-Méndez S, Huerta A (2011) Comparison of high-order curved finite elements. Int J Numer Methods Eng 87(8):719–734
Sevilla R, Rees L, Hassan O (2016) The generation of triangular meshes for NURBS-enhanced FEM. Int J Numer 108(8):941–968
Smith B, Bjorstad P, Gropp W (2004) Domain decomposition: parallel multilevel methods for elliptic partial differential equations. Cambridge University Press
Strang G, Fix G (1973) An analysis of the finite element method. Prentice-Hall
Taylor MA, Wingate BA, Vincent RE (2000) An algorithm for computing Fekete points in the triangle. SIAM J Numer Anal 38(5):1707–1720
The CGAL Project (2023) CGAL user and reference manual, vol 5.5.2. CGAL Editorial Board
The MathWorks Inc (2022) Partial Differential Equation Toolbox version: 9.4 (R2022b)
Toulorge T, Lambrechts J, Remacle J-F (2016) Optimizing the geometrical accuracy of curvilinear meshes. J Comput Phys 310:361–380
Wang D, Clark B, Jiao X (2009) An analysis and comparison of parameterization-based computation of differential quantities for discrete surfaces. Comput Aided Geom Des 26(5):510–527
Zhang Z, Naga A (2005) A new finite element gradient recovery method: superconvergence property. SIAM J Sci Comput 26(4):1192–1213
Zhu Q (1998) A review of two different approaches for superconvergence analysis. Appl Math 43(6):401–411
Zienkiewicz OC, Zhu J (1992) The superconvergent patch recovery (SPR) and adaptive finite element refinement. Comput Methods Appl Mech Eng 101(1–3):207–224
Zlámal M (1978) Superconvergence and reduced integration in the finite element method. Math Comput 32(143):663–685
Zlàmal M (2006) Some superconvergence results in the finite element method. In: Mathematical Aspects of Finite Element Methods: Proceedings of the Conference Held in Rome, December 10–12, 1975, pp 353–362. Springer
Zou X, Lo SB, Sevilla R, Hassan O, Morgan K (2023) Towards a volume mesh generator tailored NEFEM. In: International Meshing Roundtable, Amsterdam, Mar
Zwanenburg P, Nadarajah S (2017) On the necessity of superparametric geometry representation for discontinuous Galerkin methods on domains with curved boundaries. In: 23rd AIAA Computational Fluid Dynamics Conference, p 3946
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors have no relevant financial or non-financial interests to disclose.
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
Jones, J., Conley, R. & Jiao, X. Preserving superconvergence of spectral elements for curved domains via h- and p-geometric refinement. Engineering with Computers 40, 1877–1892 (2024). https://doi.org/10.1007/s00366-023-01895-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00366-023-01895-8