[논문 리뷰] Almost-$C^1$ splines: Biquadratic splines on unstructured quadrilateral meshes and their application to fourth order problems

Thomas Takacs, Deepesh Toshniwal|arXiv (Cornell University)|2022. 01. 27.

Advanced Numerical Analysis Techniques인용 수 4

한 줄 요약

이 논문은 비정규격 사각형 메esh에서 거의-$C^1$ 이차 스플라인을 도입하여, 네 번째 차수 PDE에 대한 $H^2$-비일치하는 등각유한요소 해석을 가능하게 한다. 이 방법은 모든 변과 꼭짓점에서 거의 $C^1$ 연속성을 확보하며, 임의의 메쉬 레이아웃을 지원하고 최적의 근사 성질을 가진 잘 조절된, B-spline 유사 기저를 제공한다.

ABSTRACT

Isogeometric Analysis generalizes classical finite element analysis and intends to integrate it with the field of Computer-Aided Design. A central problem in achieving this objective is the reconstruction of analysis-suitable models from Computer-Aided Design models, which is in general a non-trivial and time-consuming task. In this article, we present a novel spline construction, that enables model reconstruction as well as simulation of high-order PDEs on the reconstructed models. The proposed almost-$C^1$ are biquadratic splines on fully unstructured quadrilateral meshes (without restrictions on placements or number of extraordinary vertices). They are $C^1$ smooth almost everywhere, that is, at all vertices and across most edges, and in addition almost (i.e. approximately) $C^1$ smooth across all other edges. Thus, the splines form $H^2$-nonconforming analysis-suitable discretization spaces. This is the lowest-degree unstructured spline construction that can be used to solve fourth-order problems. The associated spline basis is non-singular and has several B-spline-like properties (e.g., partition of unity, non-negativity, local support), the almost-$C^1$ splines are described in an explicit Bézier-extraction-based framework that can be easily implemented. Numerical tests suggest that the basis is well-conditioned and exhibits optimal approximation behavior.

연구 동기 및 목표

• 기하학적 제약 없이 비정규격 사각형 메쉬에서 분석에 적합한 등각유한요소 해석을 가능하게 하기 위해.
• 고차수 PDE를 위한 CAD 데이터로부터 분석에 적합한 모델을 재구성하는 데 도전하는 문제를 해결하기 위해.
• 네 번째 차수 문제에 적합한 $H^2$-비일치 방법을 지원하는 저차수 스플라인 구축을 개발하기 위해.
• 스플라인 기저가 비특이적이며 국소적으로 지지되며, 분할 단위 및 비음성과 같은 B-spline 유사 성질을 갖도록 보장하기 위해.
• Bézier 추출을 통한 실용적이고 명시적인 구현 프레임워크를 제공하기 위해.

제안 방법

• 완전히 비정규격 사각형 메쉬, 특히 특이 꼭짓점이 있는 메쉬에서 이차 스플라인을 구축하기 위해.
• 모든 꼭짓점과 대부분의 변에서 $C^1$ 연속성을 강제하고, 나머지 변에서는 거의-$C^1$ 연속성을 확보하기 위해.
• 스플라인 기저를 명시적으로 표현하고 효율적인 구현을 가능하게 하기 위해 Bézier 추출 기반 프레임워크를 사용하기 위해.
• 기저가 비음성, 국소적 지지 및 분할 단위와 같은 핵심 성질을 만족하도록 보장하기 위해.
• 스플라인 공간을 $H^2$-비일치로 설계하여 네 번째 차수 PDE에 적합하게 하기 위해.
• 최적 수렴 속도와 잘 조절된 시스템을 확인하기 위해 수치적으로 방법을 검증하기 위해.

실험 결과

연구 질문

• RQ1 비정규격 사각형 메쉬에서 이차 스플라인은 네 번째 차수 문제에 대해 거의-$C^1$ 연속성을 달성하면서도 잘 조절된 기저를 유지할 수 있는가?
• RQ2 제안된 스플라인 구축 방식은 고차수 PDE에 대해 근사 정확도와 조건수 측면에서 어떻게 성능을 발휘하는가?
• RQ3 기하학적 제약 없이 특이 꼭짓점이 있는 임의의 사각형 메쉬를 얼마나 잘 처리할 수 있는가?
• RQ4 Bézier 추출 프레임워크는 스플라인 기저의 효율적이고 안정적인 구현을 가능하게 하는가?
• RQ5 다양한 메쉬 유형에서 네 번째 차수 문제에 대해 이 방법의 수렴 행동은 어떠한가?

주요 결과

• 제안된 거의-$C^1$ 이차 스플라인은 비정규격 사각형 메쉬에서 네 번째 차수 PDE에 대해 최적의 근사 순서를 달성한다.
• 스플라인 기저는 비특이적이며 수치 실험에서 잘 조절된 행동을 보인다.
• 이 방법은 특이 꼭짓점이 있는 임의의 비정규격 메쉬를 지원한다.
• 기저는 비음성, 국소적 지지 및 분할 단위와 같은 B-spline 유사 성질을 유지한다.
• 수치 실험을 통해 최적 수렴 속도와 다양한 시험 문제에서 강건한 조건수를 확인하였다.
• Bézier 추출 프레임워크는 스플라인 기저의 간편하고 효율적인 구현을 가능하게 하였다.