[논문 리뷰] Error analysis for a finite difference scheme for axisymmetric mean curvature flow of genus-0 surfaces

Klaus Deckelnick, Robert Nürnberg|arXiv (Cornell University)|2020. 10. 19.

Geometric Analysis and Curvature Flows참고 문헌 16인용 수 4

한 줄 요약

이 논문은 고도수 0 표면의 축대칭 평균 곡률 흐름을 위한 유한차분 스킴을 제시하며, 도는 축에서의 비가역성에 대해 신중한 경계 조건 처리를 통해 문제를 해결한다. 이는 이산 L2 및 H1 노름에서 최적의 오차 한계를 확립하며, 수치 수렴 실험과 자가수축체 및 비통합 표면의 시뮬레이션을 통해 검증된다.

ABSTRACT

We consider a finite difference approximation of mean curvature flow for axisymmetric surfaces of genus zero. A careful treatment of the degeneracy at the axis of rotation for the one dimensional partial differential equation for a parameterization of the generating curve allows us to prove error bounds with respect to discrete $L^2$- and $H^1$-norms for a fully discrete approximation. The theoretical results are confirmed with the help of numerical convergence experiments. We also present numerical simulations for some genus-0 surfaces, including for a non-embedded self-shrinker for mean curvature flow.

연구 동기 및 목표

• 고도수 0 표면의 축대칭 평균 곡률 흐름을 위한 안정적이고 정확한 유한차분 스킴을 개발한다. 이 표면들은 축과 직각으로 만나는 개방 곡선이다.
• 도는 축에서 반경 좌표가 0이 되어 표준 이산화가 실패하는, 제어 방정식의 비가역성 문제를 해결한다.
• 경계에서의 특이성에도 불구하고, 완전히 이산화된 유한차분 스킴에 대해 이산 L2 및 H1 노름에서 엄밀한 오차 추정을 수립한다.
• 부드러운 초기 데이터와 특이성 있는 초기 데이터에 대한 수치 수렴 실험을 통해 이론적 오차 한계를 검증한다.
• 비통합 표면, 특히 콘 특이성이 있는 비통합 자가수축체에 대해도 스킴의 강건성을 보여준다.

제안 방법

• DeTurck의 기법을 사용하여 엄격한 쌍곡성(parabolicity)을 확보하기 위해, 생성 곡선을 반경 및 축 방향 좌표에서 일변수 편미분방정식(PDE)으로 공식화한다.
• PDE를 발산형으로 변환하여 자연스러운 변분 공식화와 안정적인 유한차분 이산화를 가능하게 한다.
• 시간에 대해 반명시적 역유클리드 시간 이산화를 적용하고, 공간에 대해 조각별 선형 유한차분을 사용하여 안정성을 확보한다.
• 축에서의 특이성을 처리하기 위해, 속도 및 곡률 항의 일致한 경계 조건을 유도하기 위해 르호프탈의 정리를 적용한다.
• 좌표 변환과 축 근처의 메esh 정밀화를 통해 정확도를 유지하고 메쉬의 교차를 방지한다.
• 균일한 공간 및 시간 그리드를 사용하는 완전히 이산화된 스킴을 구현하고, 수치 실험을 통해 수렴성을 검증한다.

실험 결과

연구 질문

• RQ1 축대칭 평균 곡률 흐름에 대해 고도수 0 표면의 경우, 축에서의 비가역성으로 인해 최적의 오차 한계를 유한차분 스킴에 대해 확립할 수 있는가?
• RQ2 도는 축에서의 PDE 특이성을 수치적으로 처리하여 정확도와 안정성을 유지할 수 있는가?
• RQ3 비통합 표면과 90도 접촉 각 조건을 위반하는 초기 데이터에 대해 유한차분 스킴이 수렴성과 안정성을 유지하는가?
• RQ4 자기유사 수축 해, 예를 들어 콘 특이성이 있는 비통합 자가수축체를 정확하게 시뮬레이션할 수 있는가?
• RQ5 스킴의 수렴 속도는 이산 L2 및 H1 노름에서 얼마이며, 이론적 예측과 일치하는가?

주요 결과

• 축에서의 비가역성에도 불구하고, 완전히 이산화된 유한차분 스킴에 대해 이산 L2 및 H1 노름에서 최적의 오차 한계가 증명된다.
• 수치 실험을 통해 스킴은 공간에서 2차 수렴, 시간에서 1차 수렴을 달성하며, J = 512 및 Δt = 10−4 조건에서 확인되었다.
• 표면적은 시간에 따라 선형적으로 감소하며, 소멸 시간은 약 1.0로 예상과 일치한다. 이는 자가수축체의 이론적 기대와 부합한다.
• 콘 특이성이 있는 비통합 초기 데이터를 성공적으로 처리한다: 외부 콘은 빠르게 매끄럽아지고, 내부 콘은 수축하는 반원형으로 진화한다.
• 초기 조건에서 위반되었더라도, 수치해는 한계에서 90도 접촉 각 조건을 유지한다. 이는 강건성과 자가수정 기능을 보여준다.
• 이론적 분석은 르호프탈의 정리를 통해 고도수 0 표면로 확장되었으며, 이전 방법이 실패한 영역에서 오차 제어가 가능해졌다.