[논문 리뷰] Analysis of a high order Trace Finite Element Method for PDEs on level set surfaces

Jörg Grande, Christoph Lehrenfeld외 1명|arXiv (Cornell University)|2016. 11. 03.

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

한 줄 요약

이 논문은 레벨 셋 함수를 기반으로 배경 체적 메esh의 등가형 매핑을 사용하여 암시적으로 정의된 표면에서 PDE를 해결하기 위한 고차원 추적 유한요소법을 제시한다. 이 방법은 최적의 $H^1(/Gamma)$-노름 오차 한계 $\mathcal{O}(h^{k+1})$를 달성하며, 고차원 이산화에서 조건수를 균일하게 유한하게 유지하는 법선 도함수 체적 정규화 기법을 도입한다.

ABSTRACT

We present a new high order finite element method for the discretization of partial differential equations on stationary smooth surfaces which are implicitly described as the zero level of a level set function. The discretization is based on a trace finite element technique. The higher discretization accuracy is obtained by using an isoparametric mapping of the volume mesh, based on the level set function, as introduced in [C. Lehrenfeld, \emph{High order unfitted finite element methods on level set domains using isoparametric mappings}, Comp. Meth. Appl. Mech. Engrg. 2016]. The resulting trace finite element method is easy to implement. We present an error analysis of this method and derive optimal order $H^1(Γ)$-norm error bounds. A second topic of this paper is a unified analysis of several stabilization methods for trace finite element methods. Only a stabilization method which is based on adding an anisotropic diffusion in the volume mesh is able to control the condition number of the stiffness matrix also for the case of higher order discretizations. Results of numerical experiments are included which confirm the theoretical findings on optimal order discretization errors and uniformly bounded condition numbers.

연구 동기 및 목표

• 표면을 암시적으로 정의된 방식으로 사용하는 고차원 유한요소법을 개발하여 표면 PDE를 해결하는 것.
• 레벨 셋 표면에서 고차원 이산화에 대해 $H^1(\Gamma)$-노름에서 최적의 수렴 속도를 달성하는 데 도전하는 것.
• 추적 FEM의 정규화 기법을 분석하고 비교하여 고차원 설정에서 강성 행렬의 조건수에 초점을 맞추는 것.
• 표면이 배경 메쉬를 가로질러도 조건수가 균일하게 유한하게 유지되는 법선 도함수 체적 정규화 기법이 성능을 유지하는지 확인하는 것.
• 복잡한 표면 매개변수화를 피하면서도 고차원 정확도와 견고한 조건수를 동시에 확보하는 실용적이고 쉽게 구현 가능한 방법을 제공하는 것.

제안 방법

• 레벨 셋 함수를 사용하여 배경 체적 메쉬의 등가형 매핑을 적용하여 표면 기하 구조의 고차원 근사치를 확보한다.
• 변형된 배경 메쉬 위에 유한요소 공간을 구성하고, 해를 추적 연산을 통해 표면으로 제한한다.
• 추적 유한요소 형식은 체적 유한요소 공간에서 유도된 기저 함수를 사용한 갈레르킨 이산화를 활용한다.
• 고차원 설정에서 강성 행렬의 조건수를 제어하기 위해 새로운 정규화 기법인 법선 도함수 체적 정규화 기법을 도입한다.
• 기존의 정규화 기법들인 가짜 보상 및 전체 기울기 체적 정규화와의 비교를 가능하게 하는 통합 프레임워크를 사용하여 분석한다.
• 이론적 분석에는 일致성 한계 및 오차 추정이 포함되며, 다항식 확장 연산자와 레벨 셋 함수의 정규성 가정에 기반한 증명이 포함된다.

실험 결과

연구 질문

• RQ1 표면의 암시적 정의와 고차원 레벨 셋 근사치만을 사용하여 고차원 추적 유한요소법을 구성할 수 있는가? (명시적 표면 매개변수화가 필요 없이)
• RQ2 표면이 배경 메쉬를 가로질러도 고차원 추적 FEM 이산화에서 조건수를 균일하게 유한하게 유지하는 정규화 기법은 무엇인가?
• RQ3 제안된 등가형 매핑 기법이 $k$차 유한요소에 대해 최적의 $H^1(\Gamma)$-노름 오차 수렴 속도 $\mathcal{O}(h^{k+1})$를 달성하는가?
• RQ4 법선 도함수 체적 정규화 기법은 조건수 및 정확도 측면에서 다른 정규화 기법과 비교하여 어떻게 성능을 내는가?
• RQ5 복잡한 표면 재구성 또는 매개변수화가 필요 없이 이 방법을 효율적으로 구현할 수 있는가?

주요 결과

• 제안된 고차원 추적 FEM은 레이마 9에서 증명되고 수치 실험으로 뒷받침된 바에 따라 최적의 $H^1(\Gamma)$-노름 오차 한계 $\mathcal{O}(h^{k+1})$를 달성한다.
• 연구된 방법들 중에서 법선 도함수 체적 정규화 기법만이 표면이 배경 메쉬를 가로질러도 고차원 이산화에서 조건수를 균일하게 유한하게 유지한다.
• 이 방법은 고차원 레벨 셋 함수의 표준 유한요소 해석과 일반적인 유한요소 어셈블리만을 사용하므로 쉽게 구현할 수 있다.
• 수치 실험은 최적의 수렴 속도와 균일하게 유한한 조건수를 모두 확인하여 이론적 분석을 검증한다.
• 등가형 매핑 기법은 명시적 표면 매개변수화나 복잡한 표면 재구성 없이도 고차원 기하 구조 근사치를 가능하게 한다.
• 분석 결과, 조각별 평면 근사치로 인한 기하 오차는 등가형 매핑을 통해 제어되며, 이는 최적의 수렴 속도에 기여한다.