[논문 리뷰] High-order finite elements on pyramids: approximation spaces, unisolvency and exactness

Nilima Nigam, Joel Phillips|arXiv (Cornell University)|2006. 10. 05.

Advanced Numerical Methods in Computational Mathematics참고 문헌 26인용 수 4

한 줄 요약

이 논문은 유리 기저 함수를 사용하여 피라미드 메esh에서 고차수 유한요소의 가족을 제시하며, $H^1$, $H(\operatorname{curl})$, $H(\operatorname{div})$, 및 $L^2$ 공간에서 유일결정성, 적합성, 정확성을 보장한다. 주요 기여는 호환성, 고차수 근사, 그리고 공명도 다이어그램 성질을 만족하는 구성으로, 하이브리드 타트라헤드론-헥사헤드론-피라미드 메쉬에서 안정적이고 적합한 유한요소 방법을 가능하게 한다.

ABSTRACT

We present a family of high-order finite element approximation spaces on a pyramid, and associated unisolvent degrees of freedom. These spaces consist of rational basis functions. We establish conforming, exactness and polynomial approximation properties.

연구 동기 및 목표

• 하이브리드 타트라헤드론-헥사헤드론 메쉬에서 사용하기 위한 피라미드 요소에서 고차수 적합 유한요소를 개발하기 위해.
• 공유 면과 모서리에서의 자유도를 일치시켜 표준 타트라헤드론 및 헥사헤드론 유한요소와 호환되도록 보장하기 위해.
• $H^1$, $H(\operatorname{curl})$, $H(\operatorname{div})$, 및 $L^2$ 공간에 대해 근사 공간의 유일결정성과 정확성을 확립하기 위해.
• de Rham 복합체의 구조를 이산 설정에서 유지하기 위해 공명도 다이어그램 성질을 만족하기 위해.
• 매달린 노드를 피하기 위해 피라미드 요소가 필요한 복잡한 3차원 기하 구조를 위한 안정적이고 고차수의 유한요소 프레임워크를 제공하기 위해.

제안 방법

• 무한한 기준 피라미드 $\Omega_\infty$ 에서의 유리 변환을 통해 유한 피라미드로의 사상에 의해 근사 공간 $\mathcal{U}^{(s),k}(\Omega)$ 를 구성한다.
• 모든 다항식 차수 $k$ 와 $s=0,1,2,3$ 에 대해 각각 $H^1$, $H(\operatorname{curl})$, $H(\operatorname{div})$, $L^2$ 공간에 대응하는 형상 함수를 유도한다.
• 무한 피라미드에서의 회전 불변 구조를 이용하여 모서리, 면, 체적에 대한 대표 형상 함수를 분류하고 생성한다.
• 이 함수들을 유한 피라미드로의 역추적을 적용하여 불변성을 유지하고 유일결정성을 보장한다.
• 자유도 기반 $\Sigma^{(s),k}$ 를 구성함으로써 유일결정성을 검증하여 형상 함수를 유일하게 결정한다.
• 미분 연산자 $\nabla$, $\nabla\times$, $\nabla\cdot$ 와의 상호작용을 보장함으로써 보간 연산자 $\Pi^{(s)}$ 가 공명도 다이어그램 성질을 만족하도록 한다.

실험 결과

연구 질문

• RQ1 표준 타트라헤드론 및 헥사헤드론 요소와 자유도를 일치시키는 고차수 적합 유한요소를 피라미드 요소에서 구성할 수 있는가?
• RQ2 제안된 유리 기저 함수가 모든 관련 소볼레프 공간에서 유일결정성과 정확성을 보장하는가?
• RQ3 de Rham 복합체에 대해 피라미드 요소에서 이산 설정에서 공명도 다이어그램 성질이 유지되는가?
• RQ4 근사 공간이 임의의 $p$ 에 대해 차수 $p$ 의 모든 다항식을 포함하는 고차수 다항식 완비성에 도달할 수 있는가?
• RQ5 타트라헤드론, 헥사헤드론, 피라미드를 포함하는 하이브리드 메쉬에서 이 구성이 안정적이고 호환되는가?

주요 결과

• 제안된 유한요소 공간 $\mathcal{U}^{(s),k}(\Omega)$ 는 $s=0,1,2,3$ 에 대해 유일결정성과 적합성을 보장하여 이산 문제의 유일한 해를 보장한다.
• 근사 공간은 모든 차수 $k$ 의 다항식을 포함하며, 고차수 근사 성질 P2 를 만족한다.
• 면과 모서리에서의 자유도가 표준 타트라헤드론 및 헥사헤드론 요소와 일치하여 호환성(P1)을 보장한다.
• 공명도 다이어그램 성질(2)이 만족되어 이산 설정에서 de Rham 복합체의 구조가 유지된다.
• 형상 함수는 무한 피라미드에서 유도된 유리 기저 함수를 사용하여 명시적으로 구성되며, 벡터 및 스칼라 장의 정확한 표현을 가능하게 한다.
• 기존의 1차 요소(예: 윌리엄스 요소)를 일반화하고 임의의 차수 $k$ 로 확장한다.