Gauss-Galerkin quadrature rules for quadratic and cubic spline spaces and their application to isogeometric analysis

被引:59
作者
Barton, Michael [1 ,2 ]
Calo, Victor Manuel [1 ,3 ,4 ]
机构
[1] King Abdullah Univ Sci & Technol, Ctr Numer Porous Media, Thuwal 239556900, Saudi Arabia
[2] BCAM, Alameda Mazarredo 14, Bilbao 48009, Basque Country, Spain
[3] Curtin Univ, Fac Sci & Engn, Western Australian Sch Mines, CSIRO Professorial Chair Computat Geosci, Kent St, Perth, WA 6102, Australia
[4] CSIRO, Mineral Resources, Kensington, WA 6102, Australia
基金
欧盟地平线“2020”;
关键词
Optimal quadrature rules; Galerkin method; Gaussian quadrature; B-splines; Isogeometric analysis; Homotopy continuation for quadrature; POLYNOMIAL SYSTEMS; SIMULATION; FLOW;
D O I
10.1016/j.cad.2016.07.003
中图分类号
TP31 [计算机软件];
学科分类号
081202 ; 0835 ;
摘要
We introduce Gaussian quadrature rules for spline spaces that are frequently used in Galerkin discretizations to build mass and stiffness matrices. By definition, these spaces are of even degrees. The optimal quadrature rules we recently derived (Barton and Cabo, 2016) act on spaces of the smallest odd degrees and, therefore, are still slightly sub-optimal. In this work, we derive optimal rules directly for even-degree spaces and therefore further improve our recent result. We use optimal quadrature rules for spaces over two elements as elementary building blocks and use recursively the homotopy continuation concept described in Barton and Cabo (2016) to derive optimal rules for arbitrary admissible numbers of elements. We demonstrate the proposed methodology on relevant examples, where we derive optimal rules for various even-degree spline spaces. We also discuss convergence of our rules to their asymptotic counterparts, these are the analogues of the midpoint rule of Hughes et al. (2010), that are exact and optimal for infinite domains. (C) 2016 Elsevier Ltd. All rights reserved.
引用
收藏
页码:57 / 67
页数:11
相关论文
共 33 条
[1]   Explicit Gaussian quadrature rules for C1 cubic splines with symmetrically stretched knot sequences [J].
Ait-Haddou, Rachid ;
Barton, Michael ;
Calo, Victor Manuel .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2015, 290 :543-552
[2]   A simple algorithm for obtaining nearly optimal quadrature rules for NURBS-based isogeometric analysis [J].
Auricchio, F. ;
Calabro, F. ;
Hughes, T. J. R. ;
Reali, A. ;
Sangalli, G. .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2012, 249 :15-27
[3]   Observations of boundary-layer wind-tunnel flow over isolated ridges of varying steepness and roughness [J].
Ayotte, KW ;
Hughes, DE .
BOUNDARY-LAYER METEOROLOGY, 2004, 112 (03) :525-556
[4]  
Barton M, 2015, GAUSSIAN QUADR UNPUB
[5]   Optimal quadrature rules for odd-degree spline spaces and their application to tensor-product-based isogeometric analysis [J].
Barton, Michael ;
Manuel Calo, Victor .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2016, 305 :217-240
[6]   Gaussian quadrature for splines via homotopy continuation: Rules for C2 cubic splines [J].
Barton, Michael ;
Calo, Victor Manuel .
JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS, 2016, 296 :709-723
[7]   Topologically guaranteed univariate solutions of underconstrained polynomial systems via no-loop and single-component tests [J].
Barton, Michael ;
Elber, Gershon ;
Hanniel, Iddo .
COMPUTER-AIDED DESIGN, 2011, 43 (08) :1035-1044
[8]   Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows [J].
Bazilevs, Y. ;
Calo, V. M. ;
Cottrell, J. A. ;
Hughes, T. J. R. ;
Reali, A. ;
Scovazzi, G. .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2007, 197 (1-4) :173-201
[9]   DynEarthSol2D: An efficient unstructured finite element method to study long-term tectonic deformation [J].
Choi, E. ;
Tan, E. ;
Lavier, L. L. ;
Calo, V. M. .
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 2013, 118 (05) :2429-2444
[10]  
Cortes A. M. A., 2015, J COMPUT SCI