[논문 리뷰] Trace Finite Element Methods for PDEs on Surfaces

Maxim A. Olshanskii, Arnold Reusken|arXiv (Cornell University)|2016. 11. 30.

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

한 줄 요약

이 논문은 표면 피팅 메esh가 필요 없이 표면 상의 PDE를 해결하기 위한 트레이스 유한요소법(TraceFEM)을 소개한다. 배경 타트라에드럴 메쉬 상의 부피 유한요소 함수의 트레이스를 사용하여 표면 PDE를 해결하며, $P_1$ 요소에서 최적의 $O(h^2)$ 수렴을 달성하고 고차원 및 변화하는 표면으로 확장되어 안정성과 정확도가 뛰어나다.

ABSTRACT

In this paper we consider a class of unfitted finite element methods for discretization of partial differential equations on surfaces. In this class of methods known as the Trace Finite Element Method (TraceFEM), restrictions or traces of background surface-independent finite element functions are used to approximate the solution of a PDE on a surface. We treat equations on steady and time-dependent (evolving) surfaces. Higher order TraceFEM is explained in detail. We review the error analysis and algebraic properties of the method. The paper navigates through the known variants of the TraceFEM and the literature on the subject.

연구 동기 및 목표

• 표면 피팅 삼각형 메쉬가 필요 없이도 표면 상의 PDE를 안정적으로 해결할 수 있는 메쉬 독립적인 유한요소 방법을 개발한다.
• 표면 기하구조가 시간에 따라 변화하는 변화하는 표면 상의 PDE 해결 과제를 다룬다.
• 고정된 배경 메쉬 상에서 트레이스 기반 유한요소 공간을 통해 고차 정확도와 효율적인 계산을 가능하게 한다.
• 표면이 배경 메쉬를 임의의 방식으로 가로질러도 안정성과 최적 수렴이 유지되도록 보장한다.
• 시간에 따라 변화하는 문제 및 연립 문제로의 방법 확장을 포함하여, 이동-확산 및 레벨셋 기반 변화 표면에 대해 적용한다.

제안 방법

• 표면 $\Gamma$를 포함하는 배경 타트라에드럴 메쉬 $Ω \subset \mathbb{R}^3$를 사용하며, 표면 기하구조와는 독립적이다.
• 연속된 $P_1$ (또는 고차원) 함수의 부피 공간에서 유한요소 공간을 음영된 표면 $\Gamma_h$ 상의 트레이스로 정의한다.
• 표면 $\Gamma$를 $C^2$ 레벨셋 함수 $\phi$ 의 영점 수준집합으로 근사하며, $\Gamma_h = \{x \in \Omega : \phi_h(x) = 0\}$ 를 정의한다. 여기서 $\phi_h$ 는 $\phi$ 의 $P_1$ 보간이다.
• 표면 $\Gamma_h$ 상에서 PDE의 약한 형태를 트레이스 공간 $V_h^\Gamma$ 를 사용하여 이산화함으로써, 다각형 근사 $\Gamma_h$ 상에서 갈레르킨 형식을 도출한다.
• 배경 메쉬에 의해 유도된 표면 삼각형 메쉬 $\mathcal{F}_h$ 상에서 수치적 적분을 적용하며, 이는 $O(h^2)$ 표면 근사가 가능하다.
• 시간에 따라 변화하는 문제의 경우, TraceFEM을 시간 이산화 기법(예: 암시적 오일러, BDF2)과 결합하고, 해의 확장을 위해 가장 가까운 점 또는 빠른 진행 방법을 사용한다.

실험 결과

연구 질문

• RQ1 표면 피팅 메쉬가 필요 없고, 표면-메쉬 교차가 임의의 방식이어도 안정적인 표면 PDE 유한요소 방법을 개발할 수 있는가?
• RQ2 음영된 표면 상에서 $P_1$ 및 고차원 요소에 대해 TraceFEM의 최적 수렴 속도는 무엇인가?
• RQ3 표면의 위상적 변화가 발생할 경우, 시간에 따라 변화하는 표면으로서의 방법을 어떻게 확장할 수 있는가?
• RQ4 반라그랑주형 또는 특성 기반 시간 이산화와 결합했을 때 TraceFEM의 안정성과 수렴 성질은 어떠한가?
• RQ5 공간-시간에 걸쳐 효율적인 해 확장을 위해, 유한요소 방법을 초과하는 하이퍼볼릭 해법(예: 빠른 진행 방법)과 모듈러하게 결합할 수 있는가?

주요 결과

• $P_1$ TraceFEM은 정적 표면 상의 라플라스-벨트라미 방정식에서 $L^2$ 및 에너지 노름에서 최적의 $O(h^2)$ 수렴을 달성한다.
• 오직 배경 메쉬의 형태 정규성에 의존하므로, 표면-메쉬 교차의 임의성에 대해 강건함을 입증한다.
• 시간에 따라 변화하는 표면의 경우, 암시적 시간 스텝과 해의 좁은 밴드로의 적절한 확장을 통해 안정성과 수렴성을 유지한다.
• BDF2 시간 이산화와 $P_1$ TraceFEM을 결합한 수치 실험은 변화하는 표면에서 공간-시간 기준으로 2차 수렴을 보였다.
• 위상적 변화를 겪는 표면을 성공적으로 처리하였으며, 인용된 연구의 수치 결과에서 이를 입증하였다.
• TraceFEM과 하이퍼볼릭 해법(예: 빠른 진행 방법)을 조합한 하이브리드 접근법은 시간 단계 간에 모듈러하고 효율적이며 정확한 해 확장을 가능하게 하였다.