[논문 리뷰] On optimal finite element schemes for biharmonic equation

Shuo Zhang|arXiv (Cornell University)|2018. 05. 10.

Advanced Numerical Methods in Computational Mathematics참고 문헌 16인용 수 5

한 줄 요약

이 논문은 평면 비라미틱 방정식에 대해 조각별 입방 및 사차 다항식을 사용하는 두 개의 비일관성 있는 유한요소 스킴을 구축하며, 일반적인 형태로 정규화된 삼각형 메esh에서 에너지 노름에서 최적 수렴 속도 $\mathcal{O}(h^2)$ 와 $\mathcal{O}(h^3)$ 를 달성한다. 이 스킴들은 이산화된 스토크스 복합체를 기반으로 하며, $\mathcal{C}^1$-일관성 요소의 한계를 극복하기 위해 표면 연속성 조건을 통해 일致성을 확보함으로써 특수한 격자 구조가 필요 없이 최적의 근사가 가능하도록 한다.

ABSTRACT

In this paper, two nonconforming finite element schemes that use piecewise cubic and piecewise quartic polynomials respectively are constructed for the planar biharmonic equation with optimal convergence rates on general shape-regular triangulations. Therefore, it is proved that optimal finite element schemes of arbitrary order for planar biharmonic equation can be constructed on general shape-regular triangulations.

연구 동기 및 목표

• 평면 비라미틱 방정식에 대해 다항식 차수 $k=3$ 와 $k=4$ 에서 최적의 유한요소 스킴의 격차를 메우기 위해.
• 특수 격자 가정 없이 일반적인 형태로 정규화된 삼각형 메쉬에서 최적 수렴 속도를 달성하는 스킴을 개발하기 위해.
• 지나친 연속성 조건으로 인해 Ciarlet의 삼중체로 정의되지 않는 유한요소 공간을 만들기 위해.
• 일반적인 삼각형 메쉬에서 임의의 순서의 최적 스킴을 비라미틱 방정식에 대해 만들 수 있음을 보여주기 위해.
• 지역 기저 함수가 없는 상황에서 안정성과 일치성을 보장하기 위해 이산화된 스토크스 복합체를 사용하는 프레임워크를 수립하기 위해.

제안 방법

• 조각별 입방 및 사차 다항식을 사용하여 각각 $B_h^3$ 와 $B_h^4$ 라는 두 개의 비일관성 유한요소 공간을 구축하며, 셀의 변에서 연속성 조건을 적용한다.
• 표면 연속성 조건을 도입: $\int_e \llbracket w_h \rrbracket p_e = 0$ for $p_e \in P_{k-3}(e)$ 와 $\int_e q_e \llbracket \partial_n w_h \rrbracket = 0$ for $q_e \in P_{k-2}(e)$ 를 통해 일치성과 최적 근사가 보장되도록 한다.
• 지역 기저 함수의 직접 구축을 대체하기 위해 이산화된 스토크스 복합체를 구조적 프레임워크로 활용하여 유한요소 공간을 정의하고 분석한다.
• 포isson 해법, 스토크스 유사 시스템, 최종 헤시안 복구를 포함하는 3단계 해법 절차를 활용하여 갈레르킨 해와 일치하는 스킴을 확보한다.
• 일반 격자에서 Fortin-Soulie와 Crouziex-Falk 쌍의 안정성을 활용하여 이산 시스템의 최적 수렴성과 잘 정의됨을 보장한다.
• 스킴이 효율적으로 구현되고 최적으로 해석될 수 있으며, 내부 페널티를 연속성 조건으로 대체함으로써 연속형 갈레르킨 방법 유사한 구현이 가능함을 보여준다.

실험 결과

연구 질문

• RQ1 일반적인 형태로 정규화된 삼각형 메쉬에서 평면 비라미틱 방정식에 대해 임의의 순서의 최적 유한요소 스킴을 만들 수 있는가?
• RQ2 조각별 다항식 차수 $k=3$ 와 $k=4$ 를 갖는 비일관성 유한요소 공간이 특수 격자 구조 없이 최적 수렴 속도를 달성할 수 있는가?
• RQ3 지역 기저 함수가 제공되지 않는 상황에서 지나친 연속성 조건을 갖는 유한요소 공간은 어떻게 구축하고 분석할 수 있는가?
• RQ4 이산화된 스토크스 복합체는 비일관성 비라미틱 스킴의 최적 근사와 안정성에 어떤 역할을 하는가?
• RQ5 제안된 스킴과 연속형 갈레르킨 방법 사이에는 구현 및 渐近적 행동 측면에서 구조적 연결 고리가 존재하는가?

주요 결과

• 제안된 $B_h^3$ 와 $B_h^4$ 유한요소 공간은 일반적인 형태로 정규화된 삼각형 메쉬에서 에너지 노름에서 최적 수렴 속도 $\mathcal{O}(h^2)$ 와 $\mathcal{O}(h^3)$ 를 달성한다.
• 일반 격자에서 Fortin-Soulie와 Crouziex-Falk 스토크스 쌍의 기초 안정성 덕분에 스킴은 안정적이고 잘 정의되어 있다.
• 지나친 연속성 조건으로 인해 Ciarlet의 삼중체로 정의되지 않으며, 이에 따라 이산화된 스토크스 복합체를 통한 간접적 구성이 필요하다.
• 수렴 분석은 오차 추정을 넘어서 유한요소 공간의 근사 성질을 증명하는 데에도 사용되며, 표준 보간 기반 접근과는 다릅니다.
• 스킴은 효율적으로 구현되고 최적으로 해석될 수 있으며, 내부 페널티를 연속성 조건으로 대체함으로써 연속형 갈레르킨 방법 유사한 구현이 가능하다.
• 모든 $k \geq 2$ 에 대해 $B_h^k$ 의 가족은 모를리 요소를 일반화하며, 임의의 $k \geq 2$ 와 고차원으로 확장 가능할 수 있어 더 넓은 이론적 프레임워크를 시사한다.