[논문 리뷰] Conforming Discrete Gradgrad-Complexes in Three Dimensions

Jun Hu, Yizhou Liang|arXiv (Cornell University)|2020. 08. 02.

Elasticity and Material Modeling참고 문헌 10인용 수 6

한 줄 요약

이 논문은 3차원에서 처음으로 정확한 이산적 의미에서의 Gradgrad-complex를 구성하는 conforming 유한요소 공간의 가닥을 구축한다. 이는 스칼라 장에 대한 고차 $H^2$-conforming 요소, $H(\operatorname{curl})$-conforming 대칭 텐서 장, 그리고 관련된 발산-conforming 텐서 공간을 개발하여 선형화된 아인슈타인-비앙키 시스템에 대해 대칭성이 강하게 확보된 안정적인 혼합 이산화를 가능하게 한다.

ABSTRACT

In this paper, the first family of conforming discrete three dimensional Gradgrad-complexes consisting of finite element spaces is constructed. These discrete complexes are exact in the sense that the range of each discrete map is the kernel space of the succeeding one. These spaces can be used in the mixed form of the linearized Einstein-Bianchi system.

연구 동기 및 목표

• 3차원에서 정확한 이산 Gradgrad-complex를 이루는 conforming 유한요소 공간의 첫 번째 가닥을 구성하는 것.
• 전기 텐서 장의 대칭성이 강하게 확보된 선형화된 아인슈타인-비앙키 시스템에 대해 안정적인 혼합 유한요소 이산화를 가능하게 하는 것.
• 고정밀도와 정확성을 갖춘 3차원에서의 유한요소 방법으로 미분 복합체 이론을 확장하는 것.
• 엄밀하고 conforming인 프레임워크를 구축하여 향후 보다 단순하고 실용적인 방법의 개발에 기초를 마련하는 것.

제안 방법

• Zenisek와 Zhang의 $C^1$-연속 다항식 기반으로 degree 9 이상의 $H^2$-conforming 유한요소를 사용한다.
• $C^0$-연속 대칭 텐서 공간에 버블 함수를 활용한 강화를 통해 $H(\operatorname{curl})$-conforming 대칭 텐서 유한요소를 구성한다.
• 정점, 변, 면과 같은 하위단순체에서 강화된 정규성 조건을 도입하여 이산 복합체의 정확성을 확보한다.
• 전기장에 대해 대칭성을 강하게 확보하는 혼합 변분형식을 사용하여 선형화된 아인슈타인-비앙키 시스템을 기반으로 한다.
• 크랭크-니콜슨 유사 접근 방식을 사용한 시간 이산화를 포함한 반-이산 및 완전 이산 스킴을 유도한다.
• 에너지 추정과 투영 오차 경계를 통해 안정성과 수렴성을 증명하며, 해의 세 번째 시간 도함수와 $L^2$-노름의 투영을 포함한 항들을 포함한다.

실험 결과

연구 질문

• RQ1 정확한 이산적 의미에서 3차원 Gradgrad-complex에 대한 conforming 유한요소 이산화를 구성할 수 있는가?
• RQ2 대칭 텐서 장에 대한 $H(\operatorname{curl})$-conforming 유한요소는 어떻게 설계되어야 하며, 이는 복합체의 구조를 유지하는가?
• RQ3 이산 복합체의 정확성을 보장하기 위해 필요한 유한요소의 정규성 조건은 무엇인가?
• RQ4 이러한 이산 복합체는 대칭성이 강하게 확보된 선형화된 아인슈타인-비앙키 시스템에 대해 혼합 유한요소 방법을 안정화하는 데 사용될 수 있는가?
• RQ5 최종적으로 도출된 완전 이산 스킴에서 달성 가능한 수렴 속도는 무엇인가?

주요 결과

• 논문은 3차원에서 정확한 부분복합체를 이룬 최초의 conforming 이산 Gradgrad-complex 가닥을 구성한다.
• 이산 복합체는 정확하다: 각 이산 미분 연산자의 범위가 다음 연산자의 핵과 정확히 일치하여 안정성과 일致성을 보장한다.
• 정규성 조건 하에서 에너지 노름에서 최적의 수렴 속도를 달성하며, 공간에서는 $O(h^k)$, 시간에서는 $O(\Delta t^2)$로 스케일링된다.
• 수렴 분석은 투영 오차와 시간 이산화 오차를 고려하며, 해의 세 번째 시간 도함수를 포함한 항들을 포함한다.
• 전기 텐서 장의 대칭성이 $H(\operatorname{curl})$-conforming 대칭 텐서 유한요소를 통해 강하게 확보된다.
• 높은 다항식 차수와 복잡성에도 불구하고, 이 구성은 향후 보다 단순하고 실용적인 방법의 개발을 위한 이론적 기초를 제공한다.