On numerical integration in isogeometric subdivision methods for PDEs on surfaces

详细信息    查看全文

• 作者：Bert Jü ; ttler^b ; ^c ; Angelos Mantzaflaris^c ; Ricardo Perl^a ; ^{ricardo.perl@ins.uni-bonn.de" class="auth_mail" title="E-mail the corresponding author} ; Martin Rumpf^a
• 关键词：Subdivision methods ; Isogeometric analysis ; PDEs on surfaces ; Numerical quadrature
• 刊名：Computer Methods in Applied Mechanics and Engineering
• 出版年：2016
• 出版时间：15 April 2016
• 年：2016
• 卷：302
• 期：Complete
• 页码：131-146
• 全文大小：2019 K

文摘

Subdivision surfaces offer great flexibility in capturing irregular topologies combined with higher order smoothness. For instance, Loop and Catmull–Clark subdivision schemes provide C²$C^2$ smoothness everywhere except at extraordinary vertices, where the generated surfaces are C¹$C^1$ smooth. The combination of flexibility and smoothness leads to their frequent use in geometric modeling and makes them an ideal candidate for performing isogeometric analysis of higher order problems with irregular topologies, e.g. thin shell problems. In the neighborhood of extraordinary vertices, however, the resulting surfaces (and the functions defined on them) are non-polynomial. Thus numerical integration based on quadrature has to be performed carefully to preserve the expected higher order consistency.

This paper presents a detailed case study of different quadrature schemes for isogeometric discretizations of partial differential equations on closed surfaces with Loop’s subdivision scheme. As model problems we consider elliptic equations with the Laplace–Beltrami and the surface bi-Laplacian operator as well as the associated eigenvalue problem for the Laplace–Beltrami operator. A particular emphasis is on the robustness of the approach in the vicinity of extraordinary vertices. Based on a series of numerical experiments, different quadrature schemes are compared. Mid-edge quadrature, which can easily be implemented via a lookup table, turns out to be a preferable choice due to its robustness and efficiency.