[논문 리뷰] Construction of analysis-suitable $G^1$ planar multi-patch parameterizations

Mario Kapl, Giancarlo Sangalli외 1명|arXiv (Cornell University)|2017. 06. 10.

Advanced Numerical Analysis Techniques참고 문헌 24인용 수 3

한 줄 요약

이 논문은 주어진 초기 C0 다중패치 기하구조에서 평면 영역을 위한 분석에 적합한 G1(AS-G1) 다중패치 매개변수화를 구성하는 방법을 제시한다. 선형 제약 조건을 가진 이차 최적화 문제를 해결하여, 경계, 정점 및 정점에서의 1차 도함수를 유지하면서도 원래 기하구조에 가능한 한 가까운 AS-G1 매개변수화를 생성한다. 수치 실험에서는 최적의 h수렴성을 입증한다.

ABSTRACT

Isogeometric analysis allows to define shape functions of global $C^{1}$ continuity (or of higher continuity) over multi-patch geometries. The construction of such $C^{1}$-smooth isogeometric functions is a non-trivial task and requires particular multi-patch parameterizations, so-called analysis-suitable $G^{1}$ (in short, AS-$G^{1}$) parameterizations, to ensure that the resulting $C^{1}$ isogeometric spaces possess optimal approximation properties, cf. [7]. In this work, we show through examples that it is possible to construct AS-$G^{1}$ multi-patch parameterizations of planar domains, given their boundary. More precisely, given a generic multi-patch geometry, we generate an AS-$G^{1}$ multi-patch parameterization possessing the same boundary, the same vertices and the same first derivatives at the vertices, and which is as close as possible to this initial geometry. Our algorithm is based on a quadratic optimization problem with linear side constraints. Numerical tests also confirm that $C^{1}$ isogeometric spaces over AS-$G^{1}$ multi-patch parameterized domains converge optimally under mesh refinement, while for generic parameterizations the convergence order is severely reduced.

연구 동기 및 목표

• 평면 영역을 위한 실용적인 분석에 적합한 G1(AS-G1) 다중패치 매개변수화를 구성하는 방법을 개발하기 위해.
• 결과 매개변수화가 초기 C0 다중패치 기하구조의 경계, 정점 및 정점에서의 1차 도함수를 유지하도록 보장하기 위해.
• AS-G1 제약 조건을 만족하면서도 초기 기하구조에 가능한 한 가까운 AS-G1 매개변수화를 생성하기 위해.
• AS-G1 매개변수화를 사용하는 C1 이소지오메트릭 공간이 최적의 h수렴 속도를 달성하는지 수치적으로 검증하기 위해.
• 소수의 변형에 민감할 수 있음에도 불구하고 AS-G1 기하구조 구성의 유연성과 실현 가능성을 입증하기 위해.

제안 방법

• AS-G1 매개변수화의 구성 문제를 선형 등식 제약 조건을 가진 이차 최적화 문제로 수식화하기 위해.
• 초기 C0 다중패치 매개변수화를 기준으로 삼아 경계, 정점 및 정점에서의 1차 도함수를 유지하기 위해.
• 패치 간의 표면 연속성 조건에서 유도된 AS-G1 제약 조건을 도입하기 위해.
• 제약 조건의 민감성으로 인해 정확도를 유지하기 위해 기호 계산을 사용하여 최적화 문제를 수치적으로 해결하기 위해.
• 패치 단위로 매개변수화를 구성하여 AS-G1 조건을 통해 경계를 넘어서도 C1 연속성을 확보하기 위해.
• 향후 작업에서 조건 수를 개선하고 민감도를 감소시키기 위해 국소 지지도가 작은 기저를 활용하기 위해.

실험 결과

연구 질문

• RQ1 주어진 경계만으로 임의의 평면 다중패치 영역에 대해 AS-G1 다중패치 매개변수화를 구성할 수 있는가?
• RQ2 AS-G1 제약 조건은 초기 기하구조에 비해 결과 매개변수화의 품질과 근접도를 어느 정도 제한하는가?
• RQ3 AS-G1 매개변수화의 사용이 이소지오메트릭 분석에서 최적의 h수렴을 보장하는가?
• RQ4 AS-G1 제약 조건은 소수의 변형에 얼마나 민감한가? 이는 기저 설계를 통해 완화될 수 있는가?
• RQ5 제안된 방법은 체적 또는 표면 기반 다중패치 영역으로 확장될 수 있는가?

주요 결과

• 제안된 방법은 초기 C0 기하구조와 시각적으로 구분되지 않지만, 모든 연속성 제약 조건을 만족하는 AS-G1 다중패치 매개변수화를 성공적으로 생성한다.
• 수치 실험 결과, AS-G1 매개변수화를 기반으로 한 C1 이소지오메트릭 공간이 L2-노름에서 순서 O(h⁴)의 최적 수렴 속도를 확보함을 확인한다.
• 반면, 일반적인 C0 매개변수화는 L2-노름에서 최대 O(h¹/²)의 극도로 낮은 수렴 속도를 보이며, L∞-노름에서는 수렴이 전혀 이루어지지 않는다.
• AS-G1 제약 조건은 고품질의 매개변수화 구축을 방해하지 않으며, 기하구조 설계에서 상당한 유연성을 보여준다.
• 여러 테스트 케이스(자동차, 퍼즐 조각, 변형된 원 등)를 통해 초기 기하구조가 AS-G1가 아니더라도 이 방법이 강건하고 효과적임을 입증한다.
• AS-G1 제약 조건이 수치적 변형에 매우 민감하므로 기호 계산이 필수적임을 확인한다.