[논문 리뷰] Analysis of interior penalty discontinuous Galerkin methods for the Allen-Cahn equation and the mean curvature flow

Xiaobing Feng, Yukun Li|arXiv (Cornell University)|2013. 10. 28.

Solidification and crystal growth phenomena참고 문헌 5인용 수 4

한 줄 요약

이 논문은 앨런-클로우 방정식과 그의 날카운 계면 극한인 평균 곡률 흐름에 대해 완전히 이산화된 내부 페널티 불연속 갈레르킨(IP-DG) 방법을 개발하고 분석한다. 보다 정교한 비표준 오차 분석 기법을 도입함으로써, 저자들은 오차 한계를 $ε^{-1}$ 에 대해 다항식적으로만 의존하도록 유도하였으며, 지수적 의존성은 배제하여, 영수준 집합이 평균 곡률 흐름으로 수렴할 수 있도록 한다.

ABSTRACT

This paper develops and analyzes two fully discrete interior penalty discontinuous Galerkin (IP-DG) methods for the Allen-Cahn equation, which is a nonlinear singular perturbation of the heat equation and originally arises from phase transition of binary alloys in materials science, and its sharp interface limit (the mean curvature flow) as the perturbation parameter tends to zero. Both fully implicit and energy-splitting time-stepping schemes are proposed. The primary goal of the paper is to derive sharp error bounds which depend on the reciprocal of the perturbation parameter $ε$ (also called "interaction length") only in some lower polynomial order, instead of exponential order, for the proposed IP-DG methods. The derivation is based on a refinement of the nonstandard error analysis technique first introduced in [12]. The centerpiece of this new technique is to establish a spectrum estimate result in totally discontinuous DG finite element spaces with a help of a similar spectrum estimate result in the conforming finite element spaces which was established in [12]. As a nontrivial application of the sharp error estimates, they are used to establish convergence and the rates of convergence of the zero level sets of the fully discrete IP-DG solutions to the classical and generalized mean curvature flow. Numerical experiment results are also presented to gauge the theoretical results and the performance of the proposed fully discrete IP-DG methods.

연구 동기 및 목표

• 특수한 페널티 파rameter에 의존하는 정규화된 영역에서 안정성과 정확성을 유지하는 앨런-클로우 방정식에 대한 완전히 이산화된 IP-DG 스킴을 개발하기 위해.
• $\epsilon^{-1}$ 에 대해 오직 다항식 의존성만 가지는 오차 추정을 유도하여, 작은 $\epsilon$ 에서 수렴 분석을 약화시키는 지수적 증가를 피하기 위해.
• 영수준 집합이 $\epsilon \to 0$ 일 때 평균 곡률 흐름으로 수렴함을 증명하고, 그 수렴 속도를 규명하기 위해.
• 이산 갈레르킨 메쉬를 사용한 이산화 방법에 기반한 확산 계면 방법의 이론적 기초를 마련하기 위해.

제안 방법

• 내부 페널티 안정화를 사용하는 두 가지 완전히 이산화된 IP-DG 스킴을 설정: 하나는 완전히 암시적이고, 다른 하나는 에너지 분할 시간 적분 방식을 사용한다.
• 불연속 갈레르킨 공간에서의 이산 스펙트럼 추정에 기반한 보다 정교한 비표준 오차 분석 기법을 적용한다.
• 기존의 조인트 유한요소 공간에서의 스펙트럼 추정 결과([12] 참조)를 불연속 갈레르킨 설정으로 확장하여, $\epsilon^{-1}$ 의 의존성을 제어한다.
• 이산 에너지 감소 성질 $J^{h}_{\epsilon}(u_{h}^{\ell})$ 를 활용하여 수치적 안정성과 일致성을 검증한다.
• 영수준 집합 $\Gamma_t^\epsilon$ 의 진화를 모니터링하고 평균 곡률 흐름과 비교하기 위해 레벨셋 추적 전략을 활용한다.
• 확산 계면 프로파일과 일치시키기 위해, 초기 자료를 $\tanh(d_0(x)/\sqrt{2}\epsilon)$ 형태로 구성한 수치 실험을 수행한다.

실험 결과

연구 질문

• RQ1 IP-DG 방법이 앨런-클로우 방정식에 대해 $\epsilon^{-1}$ 에 대해 다항식 의존성만 가지는 오차 한계를 달성할 수 있는가? 지수적 의존성은 피한다.
• RQ2 영수준 집합이 $\epsilon \to 0$ 일 때 평균 곡률 흐름으로 수렴하는가?
• RQ3 이산 스펙트럼 추정에 기반한 제안된 오차 분석 기법을 조인트 유한요소 공간에서 불연속 갈레르킨 유한요소 공간으로 확장할 수 있는가?
• RQ4 영수준 집합이 평균 곡률 흐름으로 수렴하는 속도는 얼마이며, $\epsilon$ 에 어떻게 의존하는가?
• RQ5 완전히 암시적 스킴과 에너지 분할 시간 스텝 스킴은 에너지 감소와 정확성을 유지하는 데 어떻게 비교되는가?

주요 결과

• 제안된 IP-DG 방법은 선형 요소를 사용할 경우 $L^\infty(L^2)$ 및 $L^2(H^1)$ 노름에서 최적의 $h$-수렴 속도를 달성하며, 관측된 수렴 차수는 각각 약 2와 1이다.
• 수치 결과는 이산 에너지 $J^{h}_{\epsilon}(u_{h}^{\ell})$ 가 시간이 지남에 따라 단조롭게 감소함을 확인하였으며, 이는 이론적 기대와 일치한다.
• $\epsilon \to 0$ 일 때 수치 해의 영수준 집합이 평균 곡률 흐름으로 수렴함을 보였으며, 시간 스크린샷을 통해 더 큰 $\epsilon$ 에서 더 빠른 진화가 일어남을 확인하였다.
• 오차 추정은 $\epsilon^{-1}$ 에 대해 오직 다항식 순서로 의존하며, 이는 최적의 결과이며 계면 운동의 엄밀한 수렴 분석을 가능하게 한다.
• 이 방법은 비볼록 영역에서도 곡률에 의해 구동되는 복잡한 계면 역학, 예를 들어 곡선과 표면의 진화를 성공적으로 포착하였다.
• 수치 실험에서 관측된 수렴 속도는 이론적 예측과 일치하여 유도된 오차 한계의 날카로움을 검증하였다.