[논문 리뷰] Discrete Compactness for p-Version of Tetrahedral Edge Elements

Ralf Hiptmair|ArXiv.org|2009. 01. 07.

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

한 줄 요약

이 논문은 $p$-버전의 정합체 원소가 $\boldsymbol{H}(\boldsymbol{curl}, \Omega)$에 속하는 삼각형 모서리 유한요소에 대해 이산 콤���트성 성질을 증명하며, 마이켈슨 고유값 문제에 대한 갈레르킨 근사의 스펙트럼 정확성을 확립한다. 안정적인 벡터 포텐셜 구성과 공통 작용소 기반의 보간을 활용하여, $p \to \infty$일 때 이산 수렴성 하에 있는 발산 자유 부분공간의 유계수열이 $L^2(\Omega)$에서 수렴함을 보이며, 고유값과 고유함수의 점근적 수렴성을 확인한다.

ABSTRACT

We consider the first family of $\Hcurl$-conforming Nedéléc finite elements on tetrahedral meshes. Spectral approximation ($p$-version) is achieved by keeping the mesh fixed and raising the polynomial degree $p$ uniformly in all mesh cells. We prove that the associated subspaces of discretely weakly divergence free piecewise polynomial vector fields enjoy a long conjectured discrete compactness property as $p o\infty$. This permits us to conclude asymptotic spectral correctness of spectral Galerkin finite element approximations of Maxwell eigenvalue problems.

연구 동기 및 목표

• 일반적인 삼각형 메esh에서 $p$-버전 네델렉 모서리 유한요소에 대한 이산 콤팩트성 성질을 확립하기.
• 마이켈슨 고유값 문제에 대한 $p$-버전 유한요소 방법의 스펙트럼 정확성에 관한 오랜 동안의 추측을 해결하기.
• 이론적 기초를 제공하여 전자기학에서 $p$-버전 스펙트럼 갈레르킨 방법의 고유값과 고유함수 수렴성을 엄밀히 증명하기.
• $h$-버전 또는 2차원 경우에 국한된 이전 결과들과는 달리, 3차원 $hp$-유한요소에 대한 이산 콤팩트성 이론을 확장하기.

제안 방법

• 고정 메쉬와 증가하는 다항식 차수 $p$를 사용한 $p$-버전 유한요소 방법을 적용하여 스펙트럼 근사를 달성한다.
• 안정적인 $H^1$-벡터 포텐셜을 생성하기 위해 부드러운 파oincaré 사상 구성 기법을 적용하여 이산 헬름홀츠 분해에서의 안정성을 확보한다.
• 이산 컬 연산자의 커널에서의 근사 오차를 제어하기 위해 공통 작용소 기반의 보간 연산자를 활용한다.
• 이산 벡터장의 $L^2$-직교 헬름홀츠 분해를 통해 비소용성 및 조화 성분으로 분해한다.
• 이𝑟rotational 부분과 조화 부분을 별도로 분석함으로써 이산 해가 $L^2(\Omega)$에서 수렴함을 증명한다.
• 스트랑의 보조정리와 코hom로지 공간 추정을 활용하여 이산 커널 근사 오차를 제한하며, $s > 1/2$인 경우 조화 장의 $H^s$-정규성을 이용한다.

실험 결과

연구 질문

• RQ1 일반적인 삼각형 메쉬에서 $p$-버전의 삼각형 모서리 유한요소가 $p \to \infty$일 때 이산 콤팩트성 성질을 만족하는가?
• RQ2 마이켈슨 고유값 문제에 대한 스펙트럼 갈레르킨 방법이 $p$-버전 모서리 유한요소를 사용하여 스펙트럼적으로 정확함을 증명할 수 있는가?
• RQ3 이산 발산 자유 부분공간이 $L^2(\Omega)$에서 $L^2$-노름에 대해 콤팩트한가?
• RQ4 이산 벡터장의 $L^2$-직교 분해를 사용하여 고유함수의 $L^2(\Omega)$ 수렴을 증명할 수 있는가?

주요 결과

• 이산 콤팩트성 성질은 $\boldsymbol{H}(\operatorname{curl}, \Omega)$ 내에서 $p$-버전의 삼각형 모서리 유한요소에 대해 성립하며, 이는 이산 발산 자유 부분공간의 모든 유계수열이 $L^2(\Omega)$에서 수렴하는 부분수열을 가짐을 의미한다.
• 이산 고유함수의 극한은 연속적인 발산 자유 공간 $\boldsymbol{H}(\operatorname{curl}, \Omega) \cap \boldsymbol{H}_0(\operatorname{div}_{\boldsymbol{\epsilon}}0, \Omega)$에 속하며, 이는 스펙트럼 정확성을 보장한다.
• 스트랑의 보조정리와 $L^2$-오차 추정을 통해 이산 해의 비소용성 부분이 극한의 $L^2$-투영으로 수렴함을 증명한다.
• 조화 부분은 $H^s$-정규성과 $p$-로버스트한 투영 연산자의 근사 성질 덕분에 연속된 조화 장으로 수렴한다.
• 이산 커널 근사 오차는 $\mathcal{O}(\log^{3/2}p / p^s)$로 감소하며, $s > 1/2$인 경우 스펙트럼 수렴을 확인한다.
• 결과적으로, 3차원에서 마이켈슨 고유값 문제에 대한 $p$-버전 모서리 유한요소 방법의 점근적 최적성과 스펙트럼 정확성이 확인된다.