[논문 리뷰] A stabilized trace finite element method for partial differential equations on evolving surfaces

Christoph Lehrenfeld, Maxim A. Olshanskii외 1명|arXiv (Cornell University)|2017. 09. 21.

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

한 줄 요약

이 논문은 고정 배경 메esh와 시간 이산 유한 차분법을 사용하여 변화하는 표면에서 편미분방정식을 해결하기 위한 안정화된 트레이스 유한요소법(TraceFEM)을 제안한다. 이 방법은 트레이스 공식을 안정화하고 해를 체적 근처로 연장함으로써 최적 수렴성과 강건한 조건수를 확보하여 표면 위치에 관계없이 최적 오차 추정과 조건수 경계를 달성한다.

ABSTRACT

In this paper, we study a numerical method for the solution of partial differential equations on evolving surfaces. The numerical method is built on the stabilized trace finite element method (TraceFEM) for the spatial discretization and finite differences for the time discretization. The TraceFEM uses a stationary background mesh, which can be chosen independent of time and the position of the surface. The stabilization ensures well-conditioning of the algebraic systems and defines a regular extension of the solution from the surface to its volumetric neighborhood. Having such an extension is essential for the numerical method to be well-defined. The paper proves numerical stability and optimal order error estimates for the case of simplicial background meshes and finite element spaces of order $m\ge1$. For the algebraic condition numbers of the resulting systems we prove estimates, which are independent of the position of the interface. The method allows that the surface and its evolution are given implicitly with the help of an indicator function. Results of numerical experiments for a set of 2D evolving surfaces are provided.

연구 동기 및 목표

• 표면 재구성 또는 공간-시간 적분이 필요 없이 시간에 따라 변화하는 표면에서 PDE를 해결하기 위한 강건하고 비일치 유한요소 방법을 개발하기 위해.
• 변화하는 표면에서 고차 유한요소에 대해 수치적 안정성과 최적 수렴 속도를 확보하기 위해.
• 표면의 위치나 기하학적 형태에 관계없이 엄밀한 오차 추정과 조건수 경계를 제공하기 위해.
• 정확도와 안정성을 유지하면서도 음성 표면 기술을 위한 레벨셋 방법을 사용할 수 있도록 하기 위해.

제안 방법

• 방법은 정적 배경 메쉬 위에서 안정화된 트레이스 유한요소 공식을 사용하며, 각 시간 단계에서 해의 트레이스를 변화하는 표면에 제한한다.
• 조건수 제어와 표면에서 그 주변으로의 체적 연장 정의를 위해 안정화 항을 추가한다.
• 공간-시간 메쉬나 공간-시간 다양체 재구성의 필요 없이 표준 유한 차분법으로 시간 도함수를 근사한다.
• 변화하는 표면를 음성적으로 기술하기 위해 레벨셋 함수를 사용하여 기하학적 유연성을 확보한다.
• 핵심 요소는 배경 영역으로의 해의 너비 밴드 연장으로, 잘 정의된 문제와 안정성을 보장한다.
• 분석은 기하학적 사상과 트레이스 부등식을 기반으로 하여 확장된 영역에서 오차를 경계하고 최적 수렴 속도를 도출한다.

실험 결과

연구 질문

• RQ1 표면 재구성 또는 공간-시간 메쉬가 필요 없이 변화하는 표면에서 안정적이고 최적 수렴성을 확보하는 유한요소 방법을 개발할 수 있는가?
• RQ2 표면 위치에 관계없이 잘 조절된 대수적 시스템을 보장하기 위해 안정화를 어떻게 설계할 수 있는가?
• RQ3 변화하는 표면에서 고차 유한요소를 트레이스 유한요소 프레임워크 내에서 어떤 오차 추정을 유도할 수 있는가?
• RQ4 유도된 선형 시스템의 조건수를 표면의 위치나 곡률에 관계없이 경계로 둘 수 있는가?
• RQ5 레벨셋 함수를 통해 음성적으로 정의된 표면일 때, 방법의 수치적 성능은 어떻게 되는가?

주요 결과

• 삼각형 배경 메쉬에서 유한요소 공간의 차수 $ m \geq 1 $ 에 대해 최적 순서 오차 추정을 달성한다.
• 대수적 시스템의 조건수가 표면의 위치나 기하학적 형태에 관계없이 경계지어지며, 강건한 조건수를 보장한다.
• 공간-시간 다양체의 재구성 없이도 레벨셋 함수를 통해 음성 표면 기술이 가능하며, 명시적 표면 추적도 필요로 하지 않는다.
• 수치 실험을 통해 2차원 변화 표면 집합에서 최적 수렴 속도를 확인하여 이론적 결과를 검증한다.
• 분석을 통해 오차 추정의 법선 도함수 항이 안정화를 통해 제어되며, 최종 오차 추정은 $ \delta_n $, $ h $, 및 $ \Delta t $ 에 의존함을 규명하였다.
• 표준 가정 하에 방법은 안정적이고 수렴성이 입증되었으며, 최종 오차 추정은 메쉬 크기와 시간 간격에 대해 최적 비율로 스케일링된다.