[논문 리뷰] Isogeometric analysis with $C^1$ functions on unstructured quadrilateral meshes

Mario Kapl, Giancarlo Sangalli외 1명|arXiv (Cornell University)|2018. 12. 21.

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

한 줄 요약

이 논문은 분석에 적합한 $G^1$ 다면체 매개변수화를 사용하여 비구조적 사각형 메esh에서 $C^1$-연속적인 이소지오메트릭 스퍼지 공간을 구성하는 방법을 제안한다. 이를 통해 4차 PDE의 최적 갈레르킨 이산화가 가능해지며, 특이점이 있는 정점 근처의 스퍼지 공간을 풍부화시켜 분석에 적합하지 않은 $G^1$ 기하 구조로의 접근을 확장하여 수치 실험에서 최적 수렴 속도를 달성한다.

ABSTRACT

In the context of isogeometric analysis, globally $C^1$ isogeometric spaces over unstructured quadrilateral meshes allow the direct solution of fourth order partial differential equations on complex geometries via their Galerkin discretization. The design of such smooth spaces has been intensively studied in the last five years, in particular for the case of planar domains, and is still task of current research. In this paper, we first give a short survey of the developed methods and especially focus on the approach [26]. There, the construction of a specific $C^1$ isogeometric spline space for the class of so-called analysis-suitable $G^1$ multi-patch parametrizations is presented. This particular class of parameterizations comprises exactly those multi-patch geometries, which ensure the design of $C^1$ spaces with optimal approximation properties, and allows the representation of complex planar multi-patch domains. We present known results in a coherent framework, and also extend the construction to parametrizations that are not analysis-suitable $G^1$ by allowing higher-degree splines in the neighborhood of the extraordinary vertices and edges. Finally, we present numerical tests that illustrate the behavior of the proposed method on representative examples.

연구 동기 및 목표

• 비구조적 평면 다면체 영역에서 직접 4차 PDE의 갈레르킨 해법을 위한 전역적으로 $C^1$-연속적인 이소지오메트릭 공간을 개발하는 것.
• 특이점이 있는 정점이 있는 복잡한 비구조적 사각형 메쉬에서 패치 간면에서 $C^1$-연속적인 스퍼지 함수를 구성하는 데 도전하는 것.
• 분석에 적합한 $G^1$ 매개변수화를 초월하여 고차 스퍼지 함수를 특이점 근처에 적용함으로써 기존 $C^1$ 구성 방법을 일반적인 다면체 기하 구조로 확장하는 것.
• 대표적인 다면체 영역에서 비하모닉 및 관련 PDE를 위한 최적 수렴 속도를 수치적으로 검증하는 것.

제안 방법

• 패치 간면에서의 전역적 $C^1$ 연속성을 확보하기 위해 $C^1$-연속적인 이소지오메트릭 함수와 $G^1$-연속적인 그래프 표면 간의 등가성을 활용한다.
• 최적의 근사 성질을 보장하고 복잡한 평면 영역 표현을 허용하는 분석에 적합한 $G^1$ 다면체 매개변수화의 프레임워크를 적용한다.
• 특이점 정점과 모서리 주변에서 국소적으로 스퍼지 차수를 증가시켜 분석에 적합하지 않은 $G^1$ 기하 구조로의 구성 확장을 수행한다.
• 고차 영역에서 차수를 감소시켜 고차 요소의 수를 최소화하는 보완된 공간 구성 기법을 활용하지만, 이로 인해 비중첩 공간이 발생한다.
• 강제로 딜리클레 및 뉴먼 경계 조건을 $L^2$ 투영을 통해 구현하는 표준 갈레르킨 이산화를 비하모닉 방정식에 적용한다.
• 메쉬 크기 $h = 1/4, 1/8, 1/16, 1/32$ 를 사용한 수치적 균형 연구를 통해 다항식 차수 $p = 3$ 과 $p = 4$ 에 대해 수렴 속도를 검증한다.

실험 결과

연구 질문

• RQ1 특이점이 있는 정점이 있는 비구조적 사각형 메쉬에서 $C^1$-연속적인 이소지오메트릭 공간을 구성하면서도 최적의 근사 성질을 유지할 수 있는 방법은 무엇인가?
• RQ2 분석에 적합한 $G^1$ 매개변수화를 초월하여 일반적인 다면체 기하 구조로의 $C^1$ 구성 방법을 확장하기 위해 필요한 수정 사항은 무엇인가?
• RQ3 비하모닉 방정식과 같은 4차 PDE를 풀이할 때 도출된 $C^1$ 이소지오메트릭 공간의 수렴 행동은 어떠한가?
• RQ4 특이점 정점 근처에서 국소적으로 차수를 증가시키는 것이 이산 공간의 조건수와 근사 품질에 어떤 영향을 미치는가?
• RQ5 차수 감소로 인해 비중첩 공간이 발생하더라도, 방법이 여전히 최적의 수렴 속도를 유지할 수 있는가?

주요 결과

• 시험한 AS-$G^1$ 다면체 기하 구조에서 $p=3$ 과 $p=4$ 모두에 대해 $L^2$-노름에서 $\mathcal{O}(h^{p+1})$, $H^1$-노름에서 $\mathcal{O}(h^p)$, $H^2$-노름에서 $\mathcal{O}(h^{p-1})$ 의 최적 수렴 속도를 달성한다.
• 수치 결과는 확장된 방법을 통해 구성된 $C^1$ 이소지오메트릭 공간이 분석에 적합하지 않은 $G^1$ 매개변수화에서도 여전히 최적의 근사 성질을 유지함을 확인한다.
• 특이점 정점 근처에서 고차 스퍼지 함수를 사용함으로써 기존 방법이 실패하는 복잡한 기하 구조에서도 패치 간면에서 $C^1$ 연속성이 가능해진다.
• 공간 구성 과정으로 인해 고차 영역에서 차수를 감소시켜 비중첩 이산 공간 가족이 발생하지만, 메쉬 균형 수준 간에도 여전히 최적의 수렴이 유지된다.
• 이 방법은 선택된 테스트 케이스에서 정확한 해가 기계 정밀도 수준까지 복원되는 비하모닉 방정식의 갈레르킨 이산화를 성공적으로 가능하게 한다.