[논문 리뷰] Dual-Primal Isogeometric Tearing and Interconnecting Solvers for large-scale systems of multipatch continuous Galerkin IgA equations

Christoph Hofer, Ulrich Langer|arXiv (Cornell University)|2015. 11. 23.

Advanced Numerical Analysis Techniques참고 문헌 26인용 수 4

한 줄 요약

이 논문은 비균일 기하학적 매핑을 가진 대규모 다패치 연속 갈레르킨 이소지오메트릭 해석(IgA) 문제에 대해 이중-이중 이소지오메트릭 펴내기 및 연결 방법(IETI-DP)을 확장한다. C0 연속성 하에서 패치 경계에서의 강성 스케일링 조건수 예측을 다항로그 조건수 경계로 증명하였으며, 2차원 및 3차원 문제에서 다항수와 점프 계수의 변화에 따라 수치 실험을 통해 강건한 성능을 검증하였다.

ABSTRACT

The dual-primal isogeometric tearing and interconnecting (IETI-DP) method is the adaption of the dual-primal finite element tearing and interconnecting (FETI-DP) method to isogeometric analysis of scalar elliptic boundary value problems like, e.g., diffusion problems with heterogeneous diffusion coefficients. The purpose of this paper is to extent the already existing results on condition number estimates to multi-patch domains, which consist of different geometrical mappings for each patch. Another purpose is to prove a polylogarithmic condition number bound for the preconditioned system with stiffness scaling in case of $C^0$ smoothness across patch interfaces. Numerical experiments validate the presented theory.

연구 동기 및 목표

• 각 패치마다 다른 기하학적 매핑을 가진 다패치 영역에 대해 IgA에서 BDDC 조건수 추정을 기존의 것에서 확장한다.
• 패치 경계에서 C0 연속성 하에서 강성 스케일링을 사용한 IETI-DP 방법에 대해 다항로그 조건수 경계를 유도한다.
• 다양한 다항수와 비균질 확산 계수를 가진 2D 및 3D 영역에서의 수치 실험을 통해 이론적 결과를 검증한다.
• 점프하는 확산 계수 존재 시 IETI-DP 방법이 강성 스케일링과 함께 강건한가를 보여준다.

제안 방법

• B-스플라인 및 NURBS 기저 함수를 활용하여 약한 적합성 갈레르킨 이산화를 위한 FETI-DP 프레임워크를 이소지오메트릭 해석(IgA)에 적응시킨다.
• 기하학적 분할을 통해 상호 겹치지 않는 패치로 영역 분할을 수행하며, 각 패치는 별개의 기하학적 매핑을 통해 매개변수 영역에서 변환된다.
• 정점 평가 및 변의 평균 기반의 원시 제약 조건을 사용하여 패치 경계에서의 연속성을 강제하는 이중-이중 접근법을 구현한다.
• 스체우르 보조 행렬 시스템에서 강성 스케일링을 적용하여 조건수 경계를 향상시키며, H/h에 대해 다항로그 의존성을 달성한다.
• 패치 수준의 해법자로 국소 희소 LU 분해를 사용하고, 전역 경계 문제에 대해 반복적 Krylov 해법자(PCM)를 적용한다.
• G+SMO 라이브러리와 Eigen의 희소 LU를 사용하여 방법을 구현하며, 성능는 벽시계 시간과 반복 횟수로 측정된다.

실험 결과

연구 질문

• RQ1 각 패치마다 다른 기하학적 매핑을 가진 다패치 IgA 영역에 대해 IETI-DP의 조건수 경계를 확장할 수 있는가?
• RQ2 패치 경계에서 C0 연속성 하에서 강성 스케일링이 IETI-DP 방법에 대해 다항로그 조건수 경계를 제공하는가?
• RQ3 다양한 메쉬 세분화(H/h)와 다항수에서 다패치 IgA 문제의 조건수 및 반복 횟수는 어떻게 변화하는가?
• RQ4 패치 경계에서 점프하는 확산 계수 존재 시 IETI-DP 방법이 강성 스케일링과 함께 강건한가?
• RQ5 자유도 수가 증가함에 따라 2D 및 3D 환경에서 IETI-DP 방법의 계산 성능 및 확장성은 어떠한가?

주요 결과

• 강성 스케일링을 적용한 IETI-DP 방법은 C0 연속성 다패치 IgA 문제에서 O(log^2(H/h))의 다항로그 조건수 경계를 달성한다.
• 균일 확산 계수를 가진 2D 문제에서는 다항수 10일 때도 조건수가 3.5 이하로 유지되며, 반복 횟수는 11~15 사이에서 안정되어 있다.
• 점프하는 확산 계수에 대해서도 강건성을 보이며, 다양한 계수 비율에서도 낮은 조건수와 반복 횟수를 유지한다.
• 수치 실험 결과, 조건수는 다항수 증가에 따라 느리게 증가하며, 2D에서는 다항수 10일 때 3.5 이하, 3D에서는 4.0 이하로 유지된다.
• 성능 분석 결과, 실행 시간의 약 90%가 국소 행렬 조립 및 LU 분해 계산에 소요되며, 이는 3D에서 근사 해법자의 필요성을 시사한다.
• 실제로 잘 확장되며, 문제 크기가 증가함에 따라 조건수 및 반복 횟수가 느리게 증가하여 이론적 예측을 확인한다.