[논문 리뷰] Stabilized CutFEM for the Convection Problem on Surfaces

Erik Burman, Peter Hansbo외 2명|arXiv (Cornell University)|2015. 11. 07.

Advanced Numerical Methods in Computational Mathematics참고 문헌 17인용 수 7

한 줄 요약

이 논문은 배경 테트라헤드론 메esh에서 연속적인 조각별 선형 유한요소를 사용하여 표면 상의 대류 문제를 위한 안정화된 컷 유한요소 방법(CutFEM)을 제시한다. 이 방법은 전체 요소 면에 걸쳐 기울기 점프 안정화를 적용하여 이산 안정성과 선형 시스템의 대수적 안정성을 모두 확보하며, 에너지 노름에서 최적의 $h^{3/2}$ 수렴 속도와 $L^2$ 노름에서 전체 기울기의 $h^{3/4}$ 수렴 속도를 달성한다. 조건수는 $h^{-2}$로 유계이다. 이 접근법은 표면에 맞는 메쉬가 필요 없이 표준 유한요소 분석을 컷 메쉬로 확장한다.

ABSTRACT

We develop a stabilized cut finite element method for the convection problem on a surface based on continuous piecewise linear approximation and gradient jump stabilization terms. The discrete piecewise linear surface cuts through a background mesh consisting of tetrahedra in an arbitrary way and the finite element space consists of piecewise linear continuous functions defined on the background mesh. The variational form involves integrals on the surface and the gradient jump stabilization term is defined on the full faces of the tetrahedra. The stabilization term serves two purposes: first the method is stabilized and secondly the resulting linear system of equations is algebraically stable. We establish stability results that are analogous to the standard meshed flat case and prove $h^{3/2}$ order convergence in the natural norm associated with the method and that the full gradient enjoys $h^{3/4}$ order of convergence in $L^2$. We also show that the condition number of the stiffness matrix is bounded by $h^{-2}$. Finally, our results are verified by numerical examples.

연구 동기 및 목표

• 표면에 맞는 메쉬가 필요 없는 표면 상의 대류 문제를 위한 안정적이고 수렴하는 유한요소 방법을 개발하기 위해.
• 표준 유한요소의 안정성 및 수렴 이론을 적합하지 않은, 임의의 배경 메쉬를 사용하는 컷 유한요소 방법으로 확장하기 위해.
• 이산 수식에서 발생하는 선형 시스템의 대수적 안정성을 확보하기 위해.
• 특히 컷 메쉬에서도 에너지 노름과 전체 기울기에서 최적의 수렴 속도를 달성하기 위해.
• 토러스 표면에서의 수치 실험을 통해 이론적 결과를 검증하기 위해.

제안 방법

• 메쉬는 $\mathbb{R}^3$ 내의 배경 테트라헤드론 메쉬에서 정의된 연속적인 조각별 선형 유한요소 공간을 사용하며, 표면은 레벨셋으로 표현된다.
• 변분 수식은 표면 적분만 포함하며, 해와 시험 함수는 이산 표면 $\Gamma_h$에 제한된다.
• 전체 배경 테트라헤드론 요소 면에 걸쳐 기울기 점프 안정화 항을 추가하여 방법의 안정성과 대수적 안정성을 확보한다.
• 안정성과 일致성을 균형 잡기 위해 안정화 매개변수 $c_F$를 선택하며, 수치 실험에서는 $c_F = 10^{-2}$를 사용한다.
• 표면 미분 연산자를 정의하기 위해 가장 가까운 점 연장 기법을 사용하고, 이산 함수와 연속 함수를 연결하기 위해 리프팅 기법을 활용한다.
• 분석은 복잡도가 유계인 요소 집합을 다루기 위해 커버링 추론을 사용하여 안정성 추정치를 국소화한다.

실험 결과

연구 질문

• RQ1 표면에 맞는 메쉬가 필요 없이, 안정화된 컷유한요소법이 표면 상의 대류 문제에서 최적의 수렴 속도를 달성할 수 있는가?
• RQ2 전체 요소 면에 걸친 기울기 점프 안정화가 컷 유한요소 설정에서 이산 안정성과 대수적 안정성을 모두 보장하는가?
• RQ3 에너지 노름과 $L^2$ 노름에서 전체 기울기의 수렴 차수는 무엇인가?
• RQ4 스티프니스 행렬의 조건수는 메쉬 크기 $h$에 따라 어떻게 변화하는가?
• RQ5 유사한 안정성 메커니즘을 사용하여 표준 유한요소의 분석을 적합하지 않은 컷 메쉬로 확장할 수 있는가?

주요 결과

• 이 방법은 표면에 맞는 메쉬에서의 표준 유한요소 방법과 동일한 최적의 $h^{3/2}$ 수렴 속도를 에너지 노름에서 달성한다.
• 전체 기울기의 $L^2$ 노름 오차는 $h^{3/4}$ 속도로 수렴하며, 삼각형 메쉬에서의 스트림라인 확산 방법보다 우수하다.
• 스티프니스 행렬의 조건수는 $h^{-2}$로 유계이므로 선형 시스템의 양호한 조건성을 보장한다.
• 안정화 항은 자연 에너지 노름에서 이산 안정성을 보장하여 최적의 사전 오차 추정치 유도를 가능하게 한다.
• 수치 실험을 통해 $L^2$ 노름, 에너지 노름, 기울기 오차의 이론적 수렴 속도가 확인된다.
• 조건수 분석 결과는 渐近적 $h^{-2}$ 스케일링을 보이며, 더粗한 메쉬와 더 큰 안정화 매개변수에서는 이격 현상이 관찰된다.