×

Non-uniform order mixed FEM approximation: implementation, post-processing, computable error bound and adaptivity. (English) Zbl 1243.65135

Summary: The present work provides a straightforward and focused set of tools and corresponding theoretical support for the implementation of an adaptive high order finite element code with guaranteed error control for the approximation of elliptic problems in mixed form. The work contains: details of the discretisation using non-uniform order mixed finite elements of arbitrarily high order; a new local post-processing scheme for the primary variable; the use of the post-processing scheme in the derivation of new, fully computable bounds for the error in the flux variable; and, an \(hp\)-adaptive refinement strategy based on the a posteriori error estimator.
Numerical examples are presented illustrating the results obtained when the procedure is applied to a challenging problem involving a ten-pole electric motor with singularities arising from both geometric features and discontinuities in material properties. The procedure is shown to be capable of producing high accuracy numerical approximations with relatively modest numbers of unknowns.

MSC:

65N30 Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs
35J25 Boundary value problems for second-order elliptic equations
65N15 Error bounds for boundary value problems involving PDEs
35Q60 PDEs in connection with optics and electromagnetic theory
35Q35 PDEs in connection with fluid mechanics
Full Text: DOI

References:

[1] Schwab, C., p- and hp-Finite Element Methods, Numerical Mathematics and Scientific Computation (1998), The Clarendon Press, Oxford University Press: The Clarendon Press, Oxford University Press New York, ISBN 0-19-850390-3. Theory and applications in solid and fluid mechanics · Zbl 0910.73003
[2] Guo, B.; Babuka, I., The \(h-p\) version of the finite element method. Part 1 the basic approximation results, Comput. Mech., 1, 21-41 (1986), 10.1007/BF00298636; <http://www.dx.doi.org/10.1007/BF00298636> · Zbl 0634.73058
[3] Guo, B.; Babuka, I., The \(h-p\) version of the finite element method. Part 2 general results and applications, Comput. Mech., 1, 203-220 (1986), 10.1007/BF00272624; <http://www.dx.doi.org/10.1007/BF00272624> · Zbl 0634.73059
[4] Brezzi, F.; Fortin, M., Mixed and Hybrid Finite Element Methods (1991), Springer-Verlag New York, Inc: Springer-Verlag New York, Inc NY, USA, ISBN 0-387-97582-9 · Zbl 0788.73002
[5] Ainsworth, M.; Pinchedez, K., hp-approximation theory for BDFM and RT finite elements on quadrilaterals, SIAM J. Numer. Anal., 40, 6, 2047-2068 (2002), (electronic) (2003). doi:10.1137/S0036142901391128. URL <http://www.dx.doi.org/10.1137/S0036142901391128> · Zbl 1035.65119
[6] Demkowicz, L., Computing with hp-adaptive finite elements, Chapman & Hall/CRC Applied Mathematics and Nonlinear Science Series, vol. 1 (2007), Chapman & Hall/CRC: Chapman & Hall/CRC Boca Raton, FL, ISBN 978-1-58488-671-6; 1-58488-671-4. One and two dimensional elliptic and Maxwell problems, With 1 CD-ROM (UNIX) · Zbl 1111.65103
[7] Demkowicz, L.; Kurtz, J.; Pardo, D.; Paszyński, M.; Rachowicz, W.; Zdunek, A., Computing with hp-adaptive finite elements, Chapman & Hall/CRC Applied Mathematics and Nonlinear Science Series, vol. 2 (2008), Chapman & Hall: Chapman & Hall CRC, Boca Raton, FL, ISBN 978-1-58488-672-3; 1-58488-672-2. Frontiers: three dimensional elliptic and Maxwell problems with applications · Zbl 1148.65001
[8] Arnold, D. N.; Brezzi, F., Mixed and nonconforming finite element methods: implementation, postprocessing and error estimates, RAIRO Modél Math Anal Numér, 19, 1, 7-32 (1985) · Zbl 0567.65078
[9] Stenberg, R., Postprocessing schemes for some mixed finite elements, RAIRO Modél Math Anal Numér, 25, 1, 151-167 (1991) · Zbl 0717.65081
[10] Alonso, A., Error estimators for a mixed method, Numer. Math., 74, 4, 385-395 (1996), doi:10.1007/s002110050222. <http://www.dx.doi.org/10.1007/s002110050222> · Zbl 0866.65068
[11] Lovadina, C.; Stenberg, R., Energy norm a posteriori error estimates for mixed finite element methods, Math. Comput., 75, 256, 1659-1674 (2006), (electronic). doi:10.1090/S0025-5718-06-01872-2. URL <http://www.dx.doi.org/10.1090/S0025-5718-06-01872-2> · Zbl 1119.65110
[12] Braess, D.; Verfürth, R., A posteriori error estimators for the Raviart-Thomas element, SIAM J. Numer. Anal., 33, 6, 2431-2444 (1996), doi:10.1137/S0036142994264079. URL <http://www.dx.doi.org/10.1137/S0036142994264079> · Zbl 0866.65071
[13] Carstensen, C., A posteriori error estimate for the mixed finite element method, Math. Comput., 66, 218, 465-476 (1997), doi:10.1090/S0025-5718-97-00837-5. URL <http://www.dx.doi.org/10.1090/S0025-5718-97-00837-5> · Zbl 0864.65068
[14] Vohralı´k, M., Unified primal formulation-based a priori and a posteriori error analysis of mixed finite element methods, Math. Comput., 79, 272, 2001-2032 (2010), doi:10.1090/S0025-5718-2010-02375-0. URL <http://www.dx.doi.org/10.1090/S0025-5718-2010-02375-0> · Zbl 1201.65200
[15] Vohralı´k, M., A posteriori error estimates for lowest-order mixed finite element discretizations of convection-diffusion-reaction equations, SIAM J. Numer. Anal., 45, 4, 1570-1599 (2007), (electronic). doi:10.1137/060653184. URL <http://www.dx.doi.org/10.1137/060653184> · Zbl 1151.65084
[16] Ainsworth, M., A posteriori error estimation for lowest order Raviart-Thomas mixed finite elements, SIAM J. Sci. Comput., 30, 1, 189-204 (2007/08), doi:10.1137/06067331X. URL <http://www.dx.doi.org/10.1137/06067331X> · Zbl 1159.65353
[17] Ainsworth, M.; Coyle, J., Hierarchic finite element bases on unstructured tetrahedral meshes, Internat. J. Numer. Methods Eng., 58, 14, 2103-2130 (2003), doi:10.1002/nme.847. URL <http://www.dx.doi.org/10.1002/nme.847> · Zbl 1042.65088
[18] Marsden, J. E.; Hughes, T. J.R., Mathematical Foundations of Elasticity (1994), Dover Publications Inc.: Dover Publications Inc. New York, ISBN 0-486-67865-2. Corrected reprint of the 1983 original
[19] M. Ainsworth, B. Senior, Aspects of an adaptive hp-finite element method: adaptive strategy, conforming approximation and efficient solvers, Comput. Methods Appl. Mech. Eng. 1997; 150(1-4): 65-87. (doi:10.1016/S0045-7825(97)00101-1. Symposium on Advances in Computational Mechanics, vol. 2 (Austin, TX, 1997); URL <http://www.dx.doi.org/10.1016/S0045-7825(97)00101-1>; M. Ainsworth, B. Senior, Aspects of an adaptive hp-finite element method: adaptive strategy, conforming approximation and efficient solvers, Comput. Methods Appl. Mech. Eng. 1997; 150(1-4): 65-87. (doi:10.1016/S0045-7825(97)00101-1. Symposium on Advances in Computational Mechanics, vol. 2 (Austin, TX, 1997); URL <http://www.dx.doi.org/10.1016/S0045-7825(97)00101-1> · Zbl 0906.73057
[20] Kron, G., A set of principles to interconnect the solutions of physical systems, J. Appl. Phys., 24, 965-980 (1953) · Zbl 0051.19405
[21] Marini, L. D., An inexpensive method for the evaluation of the solution of the lowest order Raviart-Thomas mixed method, SIAM J. Numer. Anal., 22, 3, 493-496 (1985), doi:10.1137/0722029. URL <http://www.dx.doi.org/10.1137/0722029> · Zbl 0573.65082
[22] Stenberg, R., Some new families of finite elements for the Stokes equations, Numer. Math., 56, 8, 827-838 (1990), doi:10.1007/BF01405291. URL <http://www.dx.doi.org/10.1007/BF01405291> · Zbl 0708.76088
[23] Ainsworth, M., A posteriori error estimation for discontinuous Galerkin finite element approximation, SIAM J. Numer. Anal., 45, 4, 1777-1798 (2007), (electronic). doi:10.1137/060665993. URL <http://www.dx.doi.org/10.1137/060665993> · Zbl 1151.65083
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. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.