[논문 리뷰] Optimal quadrature rules for isogeometric analysis

Michael Bartoň, Victor M. Calo|arXiv (Cornell University)|2015. 11. 12.

Advanced Numerical Analysis Techniques참고 문헌 51인용 수 6

한 줄 요약

이 논문은 등속도 이sovolumetric 분석에서 사용되는 스퍼린 스페이스에 대한 최적의 가우시안 구적법 규칙을 유도하기 위해 호모토피 연속 기반 방법을 제안한다. 이는 유한 도메인에서 정확한 적분을 보장하면서도 최소한의 구적점 수를 확보한다. 접근 방식은 다항식에 대한 고전적 가우시안 구적법에서 시작하여, 목표로 하는 스퍼린 스페이스에 맞게 노드 벡터를 변형함으로써, 양의 가중치와 기계 정밀도 수준의 정확도를 갖는 규칙을 도출한다. 특히 저연속성 스퍼린 스퍼린의 경우, 점차적으로 점점 더 점근적 패턴에 수렴하는 경향이 뚜렷하다.

ABSTRACT

We introduce optimal quadrature rules for spline spaces that are frequently used in Galerkin discretizations to build mass and stiffness matrices. Using the homotopy continuation concept [6] that transforms optimal quadrature rules from source spaces to target spaces, we derive optimal rules for splines generated above finite domains. Starting with the classical Gaussian quadrature for polynomials, which is an optimal rule for a discontinuous odd-degree space, we derive rules for target spaces of higher continuity. We further show how the homotopy methodology handles cases where the source and target rules require different numbers of optimal quadrature points. We demonstrate it by deriving optimal rules for various odd-degree spline spaces, particularly with non-uniform knot sequences and non-uniform multiplicities. We also discuss convergence of our rules to their asymptotic counterparts, that is, the analogues of the midpoint rule of Hughes et al. [34], that are exact and optimal for infinite domains. For spaces of low continuities, we numerically show that the derived rules quickly converge to their asymptotic counterparts as the weights and nodes of a few boundary elements differ from the asymptotic values.

연구 동기 및 목표

• 유한 도메인에서 정확한 적분을 보장하면서도 최소한의 구적점 수를 갖는 스퍼린 스페이스에 대한 최적의 구적법 규칙을 개발하기 위해.
• 비균일하고 고연속성 스퍼린 스페이스에 대한 최적 규칙을 유도하는 데 있어 닫힌 형태 또는 수치적으로 안정적인 방법의 부족을 해결하기 위해.
• 양의 가중치와 높은 수치 정확도를 갖는 규칙을 제공함으로써 효율적이고 안정적인 갈레르킨 이산화를 가능하게 하기 위해.
• 유도된 규칙이 점점 더 점근적 대응체로 수렴하는 정도를 조사하기 위해, 특히 무한 도메인의 해석에 대해.
• 기존의 $ C^1 $ 스퍼린에 대한 재귀적 방법을 통합된 호모토피 프레임워크를 사용해 임의의 차수와 연속성으로 일반화하기 위해.

제안 방법

• 기존 스페이스(예: 불연속 다항식)에서 목표 스퍼린 스페이스로의 최적 구적법 규칙의 해 경로를 추적하기 위해 호모토피 연속을 활용한다.
• 모멘트 조건에서 유도된 다항식 방정식계의 영으로서 구적법 규칙을 표현함으로써, 스퍼린 스페이스 전역에서 정확한 적분을 보장한다.
• 기본 알려진 구성(예: $ C^1 $에 대한 이중 노드)에서 목표 구성(예: 단일 노드)으로 노드 벡터를 연속적으로 변형함으로써, 방정식계의 대수적 구조를 유지한다.
• 기존 스페이스와 목표 스페이스에서 필요한 구적점 수가 다를 경우, 호모토피 추적 중에 시스템을 동적으로 조정함으로써 이를 처리한다.
• 수치 연속 해법기를 활용하여 변형된 시스템의 근(구적 노드 및 가중치)를 계산함으로써 고정밀도와 양의 가중치를 확보한다.
• 임의의 스퍼린 함수에 대해 테스트하여 결과를 검증하였으며, 정확한 값과 16자리 소수까지 일치함을 확인하였다.

실험 결과

연구 질문

• RQ1 유한 도메인에서 임의의 차수와 연속성을 갖는 스퍼린 스페이스에 대해 최적의 구적법 규칙을 체계적으로 유도할 수 있는 방법은 무엇인가?
• RQ2 호모토피 연속 방법은 기존 스페이스(예: $ C^1 $ 세차 스퍼린)에서 다른 노드 다중도 및 연속성을 갖는 목표 스페이스로 최적 규칙을 효과적으로 전이시킬 수 있는가?
• RQ3 요소 수가 증가함에 따라 유도된 구적법 규칙이 점근적 대응체(예: 무한 도메인의 중점 규칙)로 수렴하는 속도는 어떠한가?
• RQ4 유도된 규칙는 부적절하거나 음수 가중치를 갖는 구적법보다 얼마나 더 양의 가중치와 수치적 안정성을 유지하는가?
• RQ5 기존 스페이스와 목표 스페이스에서 필요한 구적점 수가 다를 경우, 이 방법은 이를 처리할 수 있는가?

주요 결과

• 호모토피 연속 방법은 비균일 노드 벡터와 고연속성 스퍼린을 포함한 다양한 스퍼린 스페이스에 대해 최적의 가우시안 구적법 규칙을 성공적으로 도출하였으며, 최소한의 구적점 수를 확보하였다.
• $ C^0 $ 및 $ C^1 $ 스퍼린의 경우, 유도된 규칙는 빠르게 점근적 대응체로 수렴하며, 첫 번째 몇 개의 경계 요소만이 기계 정밀도를 초월하여 반복 패턴과 다를 뿐이다.
• 모든 유도된 규칙는 양의 가중치를 갖으며, 음수 가중치를 갖는 방법보다 수치적 안정성이 향상된다.
• 이 방법은 기계 정밀도 정확도를 달성하였다: 임의의 스퍼린 함수에 대한 테스트 적분 결과는 정확한 값과 16자리 소수까지 일치하였다.
• $ C^1 $ 오목 스퍼린과 $ C^1 $ 세차 스퍼린에 대해 유도된 규칙는 명시적인 대수적 구조를 보이며, 반복적 해법기를 사용하지 않고도 닫힌 형태로 계산이 가능하다.
• 표 10은 $ C^0 $ 및 $ C^1 $ 스퍼린의 경우, 점근적 규칙에서 다를 수 있는 경계 요소가 1~4개 뿐임을 보여주며, 이는 효율적 구현에 적합한 준반복적 구조임을 시사한다.