[논문 리뷰] Discontinuous Galerkin Isogeometric Analysis on Non-matching Segmentation: Error Estimates and Efficient Solvers

Christoph Hofer, Ulrich Langer외 1명|arXiv (Cornell University)|2016. 10. 12.

Advanced Numerical Analysis Techniques참고 문헌 1인용 수 5

한 줄 요약

이 논문은 간격과 겹침 영역을 가진 비일치하는 다패치 기하 구조에서 타원형 확산 문제를 위한 비연속 갈레르킨 이소지오메트릭 해석(dG IgA) 방법을 제안한다. 표준 수치적 플럭스가 실패할 경우, 테일러 기반 플럭스 근사치를 도입하여 패치별 문제를 결합함으로써 최적 수렴 속도와 뛰어난 성능을 달성하며, IETI-DP 솔버를 통해 효율적으로 작동한다. 이는 큰 확산 계수의 점프가 존재하는 경우에도 성능이 떨어지지 않는다.

ABSTRACT

The Isogeometric Analysis (IgA) of boundary value problems in complex domains often requires a decomposition of the computational domain into patches such that each of which can be parametrized by the so-called geometrical mapping. In this paper, we develop discontinuous Galerkin (dG) IgA techniques for solving elliptic diffusion problems on decompositions that can include non-matching parametrizations of the interfaces, i.e., the interfaces of the adjacent patches may be not identical. The lack of the exact parametrization of the patches leads to the creation of gap and overlapping regions between the patches. This does not allow the immediate use of the classical numerical fluxes that are known in the literature. The unknown normal fluxes of the solution on the non-matching interfaces are approximated by Taylor expansions using the values of the solution computed on the boundary of the patches These approximations are used in order to build up the numerical fluxes of the final dG IgA scheme and to couple the local patch-wise discrete problems. The resulting linear systems are solved by using efficient domainecomposition methods based on the tearing and interconnecting technology. We present numerical results of a series of test problems that validate the theoretical estimates presented.

연구 동기 및 목표

• 비일치하는 매개변수화로 인해 간격과 겹침이 발생하는 계산 영역을 패치로 분할한 타원형 확산 문제를 해결하는 데 도전한다.
• 분할 침해(세그멘테이션 범죄)로 인해 표준 수치적 플럭스를 적용할 수 없는 상황에서 안정적이고 정확한 dG IgA 스킴을 개발한다.
• 기하 구조의 간격과 겹침 영역에서 발생하는 근사 오차에도 불구하고 최적 수렴 속도를 확보한다.
• 특히 IETI-DP 솔버를 활용하여 유도된 선형 시스템을 효율적으로 해결할 수 있도록 보장한다.
• 정확한 해가 알려진 2D 및 3D 문제에 대해 다양한 확산 계수를 가진 경우에 메서드를 검증한다.

제안 방법

• 각 패치의 내부 해 값 기반으로 간격 및 겹침 인터페이스에서의 법선 플럭스를 테일러 전개를 이용해 근사한다.
• 이러한 테일러 기반 플럭스 근사치를 통합하여 비연속 갈레르킨 이소지오메트릭 해석(dG IgA) 스킴을 위한 새로운 수치적 플럭스를 구성함으로써, 비일치 인터페이스에서 약한 연속성을 강제한다.
• 고차원 근사와 기하 구조의 근사를 위해 B-spline 또는 NURBS 기반 이소지오메트릭 이산화를 사용한다.
• 전역 선형 시스템을 효율적으로 해결하기 위해 정점, 변, 면 평균 기반의 원천 변수를 사용하는 IETI-DP 도메인 분할 방법을 적용한다.
• 솔버의 강건성을 향상시키고 조건수를 감소시키기 위해 조건자와 강성 스케일링을 조건자에 적용한다.
• 메esh 크기 $ h $ 와 간격/겹침 영역의 최대 너비에 의존하는 dG-노름 기반 사전 오차 추정치를 유도한다.

실험 결과

연구 질문

• RQ1 패치 간 인터페이스 매개변수화가 비일치하여 간격과 겹침이 발생할 경우, dG IgA에서 최적 수렴 속도를 달성할 수 있는가?
• RQ2 정확한 인터페이스 일치가 없을 경우 수치적 플럭스는 어떻게 구성할 수 있으며, 플럭스 근사치의 영향은 정확도에 어떤 영향을 미치는가?
• RQ3 복잡한 기하 구조를 가진 비일치 패치 분할에서 dG IgA 시스템을 해결할 때 IETI-DP 솔버가 강건한 성능을 유지하는가?
• RQ4 큰 확산 계수의 점프가 이러한 비일치 설정에서 선형 시스템의 수렴성과 조건수에 어떤 영향을 미치는가?
• RQ5 원천 변수 선택 및 스케일링 전략이 이러한 문제에 대한 IETI-DP 솔버 성능에서 어떤 역할을 하는가?

주요 결과

• 제안된 dG IgA 방법은 수치 실험을 통해 이론적 기대와 일치하는 최적 수렴 속도를 보이며, dG-노름에서 수렴 차수를 달성한다.
• 조건자된 IETI-DP 시스템의 조건수는 메쉬 세분화에 관계없이 작고 거의 일정하게 유지되어 강건한 솔버 성능을 나타낸다.
• 240,424개의 자유도를 가진 큰 시스템에서도 콘jugate gradient 반복 수가 낮고 안정적으로 유지되며(예: 7–16회 반복), 성능이 떨어지지 않는다.
• 비균질한 확산 계수를 가진 3D 예제에서 조건자 스케일링이 강성 스케일링보다 더 우수한 성능을 보였다.
• 패치 인터페이스에서 큰 확산 계수의 점프가 존재하더라도 정확도와 솔버 효율성이 유지된다.
• 수치 결과는 이론적 오차 추정치를 확인하며, 메쉬를 세분화하고 간격/겹침 너비를 줄임에 따라 수렴 속도가 이론적 기대값에 가까워지는 것을 보여준다.