[논문 리뷰] Stabilization of Higher Order Cut Finite Element Methods on Surfaces

Mats G. Larson, Sara Zahedi|arXiv (Cornell University)|2017. 10. 09.

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

한 줄 요약

이 논문은 $[\mathbb{R}^d]$에 내장된 표면에서 타원형 PDE를 해결하기 위한 컷 유한요소법에 대해 새로운 안정화 항을 제안한다. 이 방법은 표면에서의 정규 도함수의 스케일링과 요소 면을 넘어서는 도함수의 점프 제어를 조합하여 고차원 요소를 안정화시키며, 강성 행렬의 조건수 조건이 $O(h^{-2})$로 스케일링되도록 보장함으로써 선형 및 고차원 요소에 대해 에너지 노름과 $L^2$ 노름에서 최적의 오차 추정치를 가능하게 한다. 이는 표면과 일반적인 $n$차원 다양체에 대해 성립한다.

ABSTRACT

We develop and analyze a stabilization term for cut finite element approximations of an elliptic second order partial differential equation on a surface embedded in $\mathbb{R}^d$. The new stabilization term combines properly scaled normal derivatives at the surface together with control of the jump in the normal derivatives across faces and provides control of the variation of the finite element solution on the active three dimensional elements that intersect the surface. We show that the condition number of the stiffness matrix is $O(h^{-2})$, where $h$ is the mesh parameter. The stabilization term works for linear as well as for higher-order elements and the derivation of its stabilizing properties is quite straightforward, which we illustrate by discussing the extension of the analysis to general $n$-dimensional smooth manifolds embedded in $\mathbb{R}^d$, with codimension $d-n$. We also formulate properties of a general stabilization term that are sufficient to prove optimal scaling of the condition number and optimal error estimates in energy- and $L^2$-norm. We finally present numerical studies confirming our theoretical results.

연구 동기 및 목표

• 표면에서 타원형 PDE를 해결할 때 컷 유한요소법이 고차원 요소에 대해 불안정하고 조건수가 나쁜 문제를 해결하기 위해.
• 표면과 교차하는 활성 3차원 요소에서의 해의 변동성을 제어하는 안정화 메커니즘을 개발하기 위해.
• 강성 행렬의 조건수 조건이 요소 차수에 관계없이 최적 스케일링되도록 보장하기 위해.
• 코디멘션 $d-n$을 가진 $[\mathbb{R}^d]$에 내장된 일반적인 $n$차원 매끄러운 다각체로 안정화 프레임워크를 확장하기 위해.
• 에너지 노름과 $L^2$ 노름에서 최적의 오차 추정치를 도출하기 위한 일반 안정화 항에 대한 충분한 조건을 도출하기 위해.

제안 방법

• 표면에서 적절하게 스케일링된 정규 도함수와 요소 면을 넘어서는 정규 도함수의 점프 제어를 조합한 안정화 항을 도입하기 위해.
• 표면과 교차하는 활성 3차원 요소에 이 안정화 항을 적용하여 임베딩 공간 내 해의 변동성을 제어하기 위해.
• 강성 행렬의 조건수에 대한 이론적 경계를 유도하여, 이 조건수가 $O(h^{-2})$로 스케일링됨을 증명하기 위해.
• 유사한 안정화 원칙을 사용하여 $[\mathbb{R}^d]$에 내장된 코디멘션 $d-n$을 가진 $n$차원 매끄러운 다각체로 분석을 일반화하기 위해.
• 에너지 노름과 $L^2$ 노름에서 조건수의 최적 스케일링과 최적 오차 추정치를 보장하기 위한 일반 안정화 항에 대한 충분한 조건을 설정하기 위해.
• 이론적 결과의 타당성을 검증하기 위해 수치 실험을 수행하여 조건수 스케일링과 수렴 속도를 검증하기 위해.

실험 결과

연구 질문

• RQ1 표면에서 컷 유한요소법의 최적 조건수 스케일링을 보장하기 위해 안정화 항을 어떻게 설계할 수 있는가?
• RQ2 어떤 조합의 정규 도함수 항과 점프 페널티가 컷 표면에서 고차원 유한요소 근사치를 안정화시키는가?
• RQ3 코디멘션 $d-n$을 가진 $[\mathbb{R}^d]$에 내장된 일반적인 $n$차원 다각체로 안정화 프레임워크를 확장할 수 있는가?
• RQ4 에너지 노름과 $L^2$ 노름에서 최적의 오차 추정치를 달성하기 위해 안정화 항에 필요한 충분한 조건은 무엇인가?
• RQ5 조건수와 수렴 속도에 대한 이론적 예측이 수치 실험에서 실제로 성립하는가?

주요 결과

• 제안된 안정화 항은 유한요소 공간의 다항식 차수에 관계없이 강성 행렬의 조건수가 $O(h^{-2})$로 스케일링됨을 보장한다.
• 표면에서 선형 및 고차원 유한요소에 대해 에너지 노름에서 최적의 오차 추정치를 달성한다.
• 안정화된 컷 유한요소법에 대해 $L^2$ 노름에서 최적의 오차 추정치가 확립된다.
• 안정화 메커니즘은 $[\mathbb{R}^d]$에 내장된 $n$차원 매끄러운 다각체로 자연스럽게 일반화되며, 최적의 조건수와 수렴성을 유지한다.
• 수치 연구는 이론적 예측을 확인하며, 기대되는 수렴 속도와 조건수 스케일링을 보여준다.
• 안정화 항은 표면과 교차하는 활성 3차원 요소에서 해의 변동성을 효과적으로 제어하여 부적절한 조건수를 방지한다.