[논문 리뷰] $\mathcal{P}_m$ Interior Penalty Nonconforming Finite Element Methods for $2m$-th Order PDEs in $\mathbb{R}^n$

Shuonan Wu, Jinchao Xu|arXiv (Cornell University)|2017. 10. 20.

Advanced Numerical Methods in Computational Mathematics참고 문헌 1인용 수 10

한 줄 요약

이 논문은 $\mathbb{R}^n$에서 $2m$-차수 PDE에 대해 최소한의 $\pi_m$ 내부 페널티 비일치 유한요소 방법을 제안하며, $\mathcal{P}_m$ 형상 함수와 내부 페널티를 사용하여 $m > n$일 때 수렴성을 확보한다. 이 방법은 $\mathcal{O}(1)$ 페널티 파ameter를 사용하여 에너지 노름에서 준최적 오차 추정을 달성하며, 일치 요소가 필요 없이 고차수 PDE의 강건하고 단순하며 유연한 이산화를 가능하게 한다.

ABSTRACT

In general $n$-dimensional simplicial meshes, we propose a family of interior penalty nonconforming finite element methods for $2m$-th order partial differential equations, where $m \geq 0$ and $n \geq 1$. For this family of nonconforming finite elements, the shape function space consists of polynomials with a degree not greater than $m$, making it minimal. This family of finite element spaces exhibits natural inclusion properties, analogous to those in the corresponding Sobolev spaces in the continuous case. By applying interior penalty to the bilinear form, we establish quasi-optimal error estimates in the energy norm. Due to the weak continuity of the nonconforming finite element spaces, the interior penalty terms in the bilinear form take a simple form, and an interesting property is that the penalty parameter needs only to be a positive constant of $\mathcal{O}(1)$. These theoretical results are further validated by numerical tests.

연구 동기 및 목표

• 임의의 차원 $n$에서 고차수 PDE($2m$-차수)에 대한 간단하고 최소한이며 민첩한 유한요소 방법을 개발하는 것.
• $m \geq 2$ 또는 $n \geq 2$일 때 $C^{m-1}$ 일치 요소를 구성하는 데 어려움이 있으며, 이는 복잡하고 비실용적인 요소로 이어지므로 이를 해결하는 것.
• $m \leq n$의 제약을 초월하여 모를리-왕-서 비일치 요소를 내부 페널티를 도입하여 $m > n$일 경우로 확장하는 것.
• 일치 유한요소의 존재를 요구하지 않고도 수렴성과 준최적 오차 추정을 보장하며, 최소한의 정규성 가정 하에 수행하는 것.

제안 방법

• 비일치 유한요소 공간에 대해 최소한의 $\mathcal{P}_m$ 다항형상 함수를 사용하여 단순성과 국소 효율성을 확보한다.
• 요소 면에서 내부 페널티 항을 구현하여 방법을 안정화시키고 $m > n$일 때 일致성 오차를 제어하며, 페널티 파am터는 $\mathcal{O}(1)$로 설정한다.
• 미분의 약한 연속성을 차수 $m-1$까지 유지하기 위해 자유도를 정의한다. 이는 연속 문제와의 일致성을 보장한다.
• 깨진 $H^m$ 내적과 도함수의 점프를 포함한 면적분을 조합한 이차형식을 정의하며, 페널티 항은 $h_F^{-2m+1}$ 비례로 스케일링된다.
• 대칭성과 강건한 양의 정부호성을 유지하면서 약한 연속성을 페널티 항을 통해 강제하는 변분 형식을 적용한다.
• 지역 다항식 차수 $m$을 갖는 비일치 유한요소 공간 $V_h$를 사용하며, 오차 분석을 위해 깨진 소볼레프 노름과 점프 연산자를 정의한다.

실험 결과

연구 질문

• RQ1 요소 수명이 $m > n$일 때 $\mathcal{P}_m$ 다항식에 기반한 최소한의 비일치 유한요소 방법을 $\mathbb{R}^n$에서 $2m$-차수 PDE에 대해 구성할 수 있는가?
• RQ2 일치 요소가 존재하지 않을 경우 수렴성과 최적 오차 추정을 보장하기 위해 내부 페널티 항을 어떻게 설계할 수 있는가?
• RQ3 이러한 비일치 방법에서 에너지 노름에서 준최적 수렴률을 달성하기 위해 필요한 최소 정규성 조건은 무엇인가?
• RQ4 메쉬 크기와의 의존 없이 $\mathcal{O}(1)$ 상수로 유지할 수 있는 페널티 파am터를 확보할 수 있는가? 이는 조건 수의 강건성을 보장한다.
• RQ5 해의 정규성이 떨어지는 비凸 영역에서도 이 방법이 최적 수렴률을 유지하는가?

주요 결과

• 일치 유한요소가 존재한다는 가정 하에, 에너지 노름에서 준최적 오차 추정 $\|u - u_h\|_h \lesssim h^s |u|_{m+s}$를 달성하며, $s = \min\{1,t\}$이다.
• 해가 $u \in H^{m+t}(\Omega)$에 속할 경우, 깨진 $H^m$ 노름에서 수렴 차수 $h^{\min(1,t)}$를 달성하며, 이는 최적 수렴률과 일치한다.
• 수치 실험 결과는 단위 정사각형에서 매끄러운 해에 대해 $H^3$ 노름에서 선형 수렴을 확인하였으며, $m=3,n=2$일 때 관측된 수렴 차수는 약 $1.00$으로, 정리 4.9를 검증한다.
• 특이 해를 갖는 L자형 도메인($u \in H^{3+1/2}$)에서, 깨진 $H^3$ 노름에서 수렴 차수 약 $0.5$를 달성하였으며, 이는 이론적 예측인 $h^{1/2}$와 일치한다.
• 페널티 파am터 $\eta$는 $\mathcal{O}(1)$ 상수로 선택 가능하며(예: $\eta = 1$), 조건 수의 강건성과 실용적 효율성을 보장한다.
• 이 방법은 강성 행렬의 대칭성과 양의 정부호성을 유지하므로 선형 시스템의 효율적 해법이 가능하다.