[논문 리뷰] Discrete Hessian complexes in three dimensions

Long Chen, Xuehai Huang|arXiv (Cornell University)|2020. 12. 20.

Black Holes and Theoretical Physics참고 문헌 37인용 수 5

한 줄 요약

이 논문은 다면체 메esh에서 다항식 텐서 공간 분해를 이용해 3차원에서 일치하는 가상 요소 헤시안 복합체를 구성함으로써 선형화된 아인슈타인-비앙키 시스템의 최적 순서 이산화를 가능하게 한다. 이 방법은 H2-일치하는 가상 요소와 H(div)-일치하는 유한 요소를 코즐 유형의 분해를 통해 결합하여, 적절히 선택된 자유도와 투영 기반 연산자를 통해 안정성과 계산 가능성을 확보하며, 에너지 노름에서 최적 수렴 속도를 달성한다.

ABSTRACT

A family of conforming virtual element Hessian complexes on tetrahedral meshes are constructed based on decompositions of polynomial tensor spaces. They are applied to discretize the linearized time-independent Einstein-Bianchi system with optimal order convergence.

연구 동기 및 목표

• 삼각형 메쉬에서 최소 자유도를 가진 3차원에서 안정적이고 일치하는 이산 헤시안 복합체를 개발하기 위해.
• 가상 요소 방법을 이용해 선형화된 시간에 의존하지 않는 아인슈타인-비앙키 시스템의 최적 순서 이산화를 가능하게 하기 위해.
• 기존의 유한 요소 헤시안 복합체에서의 고차 다항식 복잡성을 가상 요소 기법을 활용해 극복하기 위해.
• 새로운 다항식 텐서 공간 분해를 통해 이산 복합체의 잘 정의됨과 계산 가능성을 보장하기 위해.
• 핵심 구조(예: H(curl)-일치성, L2-투영의 계산 가능성)를 유지하면서도 최적 수렴을 확보하는 프레임워크를 수립하기 위해.

제안 방법

• 다각형에서의 다항식 텐서 공간에 대해 코즐 유형의 분해를 사용하여 이산 헤시안 복합체를 구성: Pk(K; S) = ∇²Pk+2(K) ⊕sym(Pk−1(K; T) × x).
• 가상 요소 공간 Σ(K) = ∇²W(K) ⊕sym(V(K) × x)를 정의하며, 여기서 W(K)는 H2-일치하는 가상 요소 공간이고, V(K) = Pk−1(K; T)는 H(div)-일치하는 유한 요소 공간이다.
• 투영 기반 접근법을 구현: QKkσ는 내부 모멘텀을 통해 계산 가능하며, eQKk+2σ는 ΠKσ로부터 재구성되어 L2-투영의 계산 가능성을 확보한다.
• 이중선형 형식 bK와 cK를 사용하여 혼합 방법의 이산 이중선형 형식을 정의하며, 이산 인프-서프 조건을 통해 안정성을 확보한다.
• 변분 형식 (91)–(92)를 사용한 일치하는 혼합 가상 요소 방법을 적용하며, 해는 국소 형식의 전역 조립을 통해 계산된다.
• 일致하고 안정적인 혼합 유한 요소 형식에 기반한 혼합 방법을 활용하며, 이는 이산 인프-서프 조건과 투영 기반 일致성 추정에 의존한다.

실험 결과

연구 질문

• RQ1 기존의 유한 요소 복합체보다 자유도가 적은 가상 요소를 사용하여 3차원 삼각형 메쉬에서 일치하는 이산 헤시안 복합체를 구성할 수 있는가?
• RQ2 완전한 경계 조건 특성화 없이도 다면체 영역에서 비균일 경계 조건을 가진 국소 H2 문제의 잘 정의됨을 어떻게 보장할 수 있는가?
• RQ3 어떤 다항식 텐서 공간 분해가 헤시안 복합체에 대해 안정적이고 일치하는 가상 요소 공간을 구성하는 데 기여하는가?
• RQ4 이 가상 요소 프레임워크를 통해 선형화된 아인슈타인-비앙키 시스템에 대해 최적 수렴 속도를 달성할 수 있는가?
• RQ5 형상 함수의 명시적 지식 없이도 가상 요소 공간의 L2-투영을 다항식 공간으로 계산 가능하게 할 수 있는가?

주요 결과

• 정규성 가정 하에, 제안된 가상 요소 헤시안 복합체는 헤시안의 H(curl)-노름과 잠재함수의 H2-노름에서 최적 순서 hk의 수렴 속도를 달성한다.
• 새로운 안정성 증명을 통해 이산 인프-서프 조건이 확립되어 혼합 형식의 잘 정의됨이 보장된다.
• 일치성 추정 (96)–(98)은 투영 기반 연산자가 관련 노름에서 각각 hk 및 hk+1의 오차를 유도함을 보여주며, 이는 상한이 존재함을 의미한다.
• 가상 요소 공간의 L2-투영이 다항식 공간 Pk(S)로의 투영이 내부 모멘텀과 재구성 연산자만을 사용하여 계산 가능함을 보장한다.
• 수렴 속도 ∥E − Eh∥H(curl) + ∥σ − σh∥2 ≲ hk(|E|k+1 + ∥σ∥k+2 + |f|k)가 증명되어 부드러운 해에 대해 최적 순서임을 확인한다.
• 기존의 유한 요소 복합체에서 요구하는 고차 다항식 기저 함수(예: P9)의 필요성을 피함으로써 자유도를 크게 감소시킨다.