High‐order cubature rules for tetrahedra
- Faculty of Civil Engineering Cracow University of Technology Cracow Poland
- Department of Civil and Environmental Engineering University of California Davis California
In this paper, we construct new high‐order numerical integration schemes for tetrahedra, with positive weights and integration points that are in the interior of the domain. The construction of cubature rules is a challenging problem, which requires the solution of strongly nonlinear algebraic (moment) equations with side conditions given by affine inequality constraints. We present a robust algorithm based on a sequence of three modified Newton procedures to solve the constrained minimization problem. In the literature, numerical integration rules for the tetrahedron are available up to order p =15. We obtain integration rules for the tetrahedron from p =2 to p =20, which are computed using multiprecision arithmetic. For p ≤15, our approach provides integration rules that have the same or fewer number of integration points than existing rules; for p =16 to p =20, our rules are new. Numerical tests are presented that verify the polynomial‐precision of the cubature rules. Convergence studies are performed for the integration of exponential, rational, weakly singular and trigonometric test functions over tetrahedra with flat and curved faces. In all tests, improvements in accuracy is realized as p is increased, though in some cases nonmonotonic convergence is observed.
- Sponsoring Organization:
- USDOE
- OSTI ID:
- 1596161
- Journal Information:
- International Journal for Numerical Methods in Engineering, Journal Name: International Journal for Numerical Methods in Engineering Journal Issue: 11 Vol. 121; ISSN 0029-5981
- Publisher:
- Wiley Blackwell (John Wiley & Sons)Copyright Statement
- Country of Publication:
- United Kingdom
- Language:
- English
Similar Records
Cubature rules of prescribed merit
Notes on lattice rules.