[논문 리뷰] Convergent semi-Lagrangian methods for the Monge-Ampère equation on unstructured grids

Xiaobing Feng, Max Jensen|arXiv (Cornell University)|2016. 02. 15.

Geometric Analysis and Curvature Flows참고 문헌 2인용 수 6

한 줄 요약

이 논문은 비구조적 삼각형 격자에서 몽체-암페르 방정식에 대해 수렴성 있는 준라그랑주 방법을 제안하며, 이를 해밀턴-자코비-벨리만 방정식으로 재구성함으로써 명시적 볼록성 제약 조건이 필요 없도록 한다. 이 방법은 단조성, 국소 일致성, 균일 안정성 있는 스킴을 통해 유일한 점성 해로의 균일 수렴을 달성하며, 비가속성 및 열화된 경우에도 호워드의 알고리즘을 사용하여 초선형 수렴을 이룬다.

ABSTRACT

This paper is concerned with developing and analyzing convergent semi-Lagrangian methods for the fully nonlinear elliptic Monge-Ampère equation on general triangular grids. This is done by establishing an equivalent (in the viscosity sense) Hamilton-Jacobi-Bellman formulation of the Monge-Ampère equation. A significant benefit of the reformulation is the removal of the convexity constraint from the admissible space as convexity becomes a built-in property of the new formulation. Moreover, this new approach allows one to tap the wealthy numerical methods, such as semi-Lagrangian schemes, for Hamilton-Jacobi-Bellman equations to solve Monge-Ampère type equations. It is proved that the considered numerical methods are monotone, pointwise consistent and uniformly stable. Consequently, its solutions converge uniformly to the unique convex viscosity solution of the Monge-Ampère Dirichlet problem. A super-linearly convergent Howard's algorithm, which is a Newton type method, is utilized as the nonlinear solver to take advantage of the monotonicity of the scheme. Numerical experiments are also presented to gauge the performance of the proposed numerical method and the nonlinear solver.

연구 동기 및 목표

• 일般적인 비구조적 격자에서 몽체-암페르 방정식의 수치적 방법에서 볼록성 조건을 강제하는 문제를 해결하기 위해.
• 이산 공간에서 볼凸성 제약 조건이 필요 없는 완전 비선형 타원형 편미분방정식을 위한 단조성, 일致성, 안정성 있는 수치 스킴을 개발하기 위해.
• 넓은 스텐실 유한차분 스킴과 호환되며, 효율적인 비선형 해법(예: 호워드의 알고리즘)을 적용할 수 있도록 단조성을 확보하기 위해.
• 비구조적 삼각형 메esh에서 몽체-암페르 딜리클레 문제의 유일한 볼록 점성 해로 수치 해의 수렴성을 확립하기 위해.
• 벨르만 재구성에 의해 비볼凸 함수로의 비교 원리를 확장하여, 전통적 프레임워크를 초월한 엄밀한 수렴 분석을 가능하게 하기 위해.

제안 방법

• 점성 의미에서 몽체-암페르 방정식을 등가의 해밀턴-자코비-벨리만(HJB) 형태로 재구성하여, 연산자 구조에 의해 자연스럽게 볼凸성이 포함되도록 한다.
• 유한요소 배경 격자 위에 넓은 스텐실 유한차분 스킴을 구성하며, 삼각형 메쉬와 고차 스텐실을 조합하여 HJB 연산자를 근사한다.
• 최대 원리와 진동을 피하기 위해 준라그랑주 이산화를 사용하여 스킴의 단조성을 확보한다.
• 국소 보간과 방향 도함수를 사용한 준라그랑주 접근을 통해 HJB 연산자를 근사함으로써 국소 일치성을 확보한다.
• 스킴의 단조성을 활용하여 강력한 수렴성을 보장하는 초선형 수렴 성질을 지닌 호워드의 알고리즘(준연속 뉴턴 방법)을 비선형 해법으로 적용한다.
• 이중 격자 접근을 사용한다: 함수 표현을 위한 굵은 유한요소 공간과 HJB 연산자의 정확한 이산화를 위한 넓은 스텐실

실험 결과

연구 질문

• RQ1 비구조적 격자에서 몽체-암페르 방정식을 점성 의미에서 해밀턴-자코비-벨리만 방정식으로 재구성할 수 있는가?
• RQ2 HJB 재구성은 이산 해공간에서 볼凸성 조건을 명시적으로 강제할 필요가 없음을 의미하는가?
• RQ3 비구조적 삼각형 메쉬에서 HJB 재구성에 대해 단조성, 일치성, 균일 안정성 있는 수치 스킴을 구성할 수 있는가?
• RQ4 비가속성 및 열화된 해에 대해 호워드의 알고리즘이 초선형 수렴 성질을 보이며 효과적으로 적용 가능한가?
• RQ5 매끄럽고 비가속성 정확해에 대해 각각 $L^2$, $L^rown$, $H^1$ 노름에서 수치 해의 수렴 속도는 어떠한가?

주요 결과

• 제안된 준라그랑주 스킴은 단조성, 국소 일치성, 균일 안정성을 확보하여 유일한 볼凸 점성 해로의 수렴을 보장한다.
• HJB 재구성은 이산 공간에서 볼凸성 조건을 강제할 필요 없이, 연산자 구조 자체에 의해 볼凸성이 내재되어 있음을 의미한다.
• 사차함수 정확해에 대해 $L^rown$-오차는 $h$에 대해 약 $1.02$의 수렴 차수와 $k$에 대해 약 $1.78$의 수렴 차수를 보이며, 시간 유사 매개변수 $k$에서 초선형 수렴을 나타낸다.
• 비가속성 해 $|x_1|$에 대해 $L^rown$-오차는 $h$에 대해 $0.64$의 수렴 차수와 $k$에 대해 $1.50$의 수렴 차수를 보이며, 정규성 부족에도 불구하고 강건함을 입증한다.
• 호워드의 알고리즘은 모든 메쉬 세밀도와 스텐실 크기에서 중간 정도의 반복 횟수(4–16회)를 보이며, 세밀한 메쉬에서도 효율적인 해를 제공한다.
• 메쉬와 일치하지 않는 특이선이 있는 경우에도 $|x_1|$과 같은 비가속성 해를 성공적으로 계산하며, 열화된 경우의 강건함을 확인한다.