[논문 리뷰] A finite element method for the Monge-Ampère equation with transport boundary conditions

Ellya L. Kawecki, Omar Lakkis외 1명|arXiv (Cornell University)|2018. 07. 10.

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

한 줄 요약

이 논문은 최적 운반 경계 조건을 갖는 몽헤-아플랑르 방정식을 해결하기 위한 비변분 유한요소법(NVFEM)을 제시한다. 뉴턴-랩슨 선형화와 기울기/헤시안 복원 기법을 통해 최적 수렴성을 달성한다. 이 방법은 조각다항식 근사에서 L2 및 H1 노름에서 경험적으로 최적 수렴률을 달성하며, 광학 및 메쉬 이동 분야의 수치 실험을 통해 검증된다.

ABSTRACT

We address the numerical solution via Galerkin type methods of the Monge-Ampère equation with transport boundary conditions arising in optimal mass transport, geometric optics and computational mesh or grid movement techniques. This fully nonlinear elliptic problem admits a linearisation via a Newton-Raphson iteration, which leads to an oblique derivative boundary value problem for elliptic equations in nondivergence form. We discretise these by employing the nonvariational finite element method, which lead to empirically observed optimal convergence rates, provided recovery techinques are used to approximate the gradient and the Hessian of the unknown functions. We provide extensive numerical testing to illustrate the strengths of our approach and the potential applications in optics and mesh movement.

연구 동기 및 목표

• 완전 비선형 몽헤-아플랑르 방정식에 대해 최적 운반 경계 조건을 갖는 갈레르킨 형식의 유한요소법을 개발한다.
• 몽헤-아플랑르 PDE의 뉴턴-랩슨 선형화에서 발생하는 비발산 형식의 타원형 방정식을 해결하는 데 도전한다.
• 기울기 및 헤시안 복원 기법을 사용하여 비변분 유한요소를 통해 강력하고 정확한 수치해를 확보한다.
• 수치 실험을 통해 기하광학 및 적응형 메쉬 이동 분야에의 적용 가능성을 입증한다.
• FEniCS와 같은 표준 유한요소 패키지를 사용하여 계산 효율성이 높고 쉽게 구현 가능한 방법을 제공한다.

제안 방법

• 이 방법은 완전 비선형 몽헤-아플랑르 방정식을 비발산 형식의 타원형 문제로 선형화하기 위해 뉴턴-랩슨 반복을 활용하며, 기울기 경계 조건이 있는 문제로 변환한다.
• 선형화된 문제를 이산화하기 위해 비변분 유한요소법(NVFEM)을 사용하여 약한 형태를 회피하고 비정규 메쉬에서도 일치 근사를 가능하게 한다.
• 해의 기울기와 헤시안을 복원하기 위해 기울기 및 헤시안 복원 기법을 적용하며, 이는 운반 지도 ∇u를 근사하는 데 필수적이다.
• 해와 그 도함수를 위한 조각다항식 유한요소 공간(P1, P2, P3)을 사용하며, 복원 기법이 수렴률을 향상시킨다.
• 경계 조건은 볼록성 가정 하에 t(∂Ω) = ∂Υ 와 동치인 운반 조건 ∇u(∂Ω) = ∂Υ 를 통해 구현된다.
• 이 방법은 FEniCS 유한요소 패키지를 사용하여 구현되어 있어 간편한 적응 및 고차원 확장이 가능하다.

실험 결과

연구 질문

• RQ1 갈레르킨 형식의 유한요소법이 최적 운반 경계 조건을 갖는 몽헤-아플랑르 방정식에서 최적 수렴률을 달성할 수 있는가?
• RQ2 기울기 및 헤시안 복원 기법이 비변분 유한요소 체계에서 최적 수렴률을 복원하는 데 얼마나 효과적인가?
• RQ3 NVFEM이 최적 운반에서 유래하는 기울기 경계 조건을 갖는 비선형 타원형 문제에 성공적으로 적용될 수 있는가?
• RQ4 이 방법은 기하광학 및 메쉬 이동과 같은 실용적 응용 분야에서 성능이 뛰어나게 작동하는가?
• RQ5 이 방법은 표준 유한요소 소프트웨어를 사용하여 강력하고 쉽게 구현 가능한가?

주요 결과

• 기울기 복원 기법을 적용한 P1, P2, P3 유한요소에서 L2 및 H1 노름에서 최적 수렴률을 달성하며, H1 및 L2 오차에 대해 관측된 수렴 차수는 약 1.8–1.9이다.
• 기울기 복원 없이 P1 요소를 사용할 경우 수렴률이 최적 수준 이하이거나 수렴하지 않음을 확인하여, 복원 기법의 필수성을 입증한다.
• P2 및 P3 요소의 경우, 헤시안 및 기울기 근사에 대해 관측된 수렴 차수는 각각 1.85 및 1.93로 최적에 가까운 값을 보였다.
• 모니터 함수 기반 메쉬 생성에 성공적으로 적용되었으며, 가스파르 몽게 초상화 예제를 통해 균일 메쉬가 목표 밀도를 따라 왜곡되는 것을 확인했다.
• 복원된 기울기를 통해 운반 지도 ∇u의 정확한 근사를 가능하게 하여 광학 및 메쉬 적응 분야에서 중요한 응용에 기여한다.
• 계산 프레임워크는 공개되어 있어 재현성과 고차원 등각요소로의 확장에 유리하다.