[논문 리뷰] Two-scale method for the Monge-Ampère Equation: Pointwise Error Estimates

Ricardo H. Nochetto, Dimitrios Ntogkas외 1명|arXiv (Cornell University)|2017. 06. 28.

Geometry and complex manifolds참고 문헌 5인용 수 4

한 줄 요약

이 논문은 차원 $d \geq 2$ 에서 몽주-암페르 방정식에 두 스케일 유한요소법을 적용한 경우의 점별 오차 추정을 수립한다. 이는 이산 알렉산드로프 추정을 통해 이산 해가 데이터에 대해 연속적으로 의존한다는 것을 증명함으로써, 힐베르트 공간에서의 $L^\infty$ 수렴 속도를 $h^{2/\beta}$ 로 도출한다. 이는 하올더 또는 소볼레프 정규성을 갖는 고전적 해에 대해 성립하며, 유계이지만 가능성이 떨어지는 헤시안을 갖는 점성 해에 대해서는 수렴 속도가 감소한다.

ABSTRACT

In this paper we continue the analysis of the two-scale method for the Monge-Ampère equation for dimension $d \geq 2$ introduced in [10]. We prove continuous dependence of discrete solutions on data that in turn hinges on a discrete version of the Alexandroff estimate. They are both instrumental to prove pointwise error estimates for classical solutions with Hölder and Sobolev regularity. We also derive convergence rates for viscosity solutions with bounded Hessians which may be piecewise smooth or degenerate.

연구 동기 및 목표

• 차원 $d \geq 2$ 에서 몽주-암페르 방정식에 두 스케일 유한요소법을 적용한 경우의 점별 $L^\infty$ 오차 추정을 수립하는 것.
• 오차 추정을 도출하기 위한 핵심 요소인 이산 해가 데이터에 대해 연속적으로 의존한다는 것을 증명하는 것.
• 하올더 및 소볼레프 정규성을 갖는 고전적 해에 대한 수렴 분석을 확장하고, 유계이지만 가능성이 떨어지는 헤시안을 갖는 점성 해에 대해서도 적용하는 것.
• 비정규성 조건 $f > 0$ 를 제거하고 $f \geq 0$ 인 탈중력 경우의 수렴 속도 분석을 수행하는 것.
• 이전 연구에서 관측된 선형 수렴 속도와 이론적 오차 추정 간의 괴리를 해소하기 위해, $\delta/h \gg 1$ 인 요소가 수렴 속도를 열 劣화시키는 원인임을 규명하는 것.

제안 방법

• 헤시안 행렬식의 최소-최대 표현을 바탕으로 한 정규직교 방향 기저를 활용한 두 스케일 수식을 사용한다.
• 비정규적인 균일 메쉬의 크기 $h$ 에서 연속적인 조각별 선형 유한요소 공간 위에 이산 해를 정의한다.
• 粗모양 스케일 $\delta$ 는 두 번째 순서 방향 차분 $\nabla^2_\delta u(x;v)$ 를 계산하는 스텐실을 결정하며, 헤시안을 근사한다.
• 이산 몽주-암페르 연산자 $T_\varepsilon[u_\varepsilon](x_i)$ 는 정규직교 기저의 이산 집합 $\mathbb{S}^\perp_\theta$ 에서의 최소값으로 정의되며, 두 번째 순서 차분의 양수 및 음수 부분을 사용한다.
• 이산 알렉산드로프 추정을 통해 이산 해가 데이터에 대해 연속적으로 의존함을 확보함으로써 오차 전파 분석이 가능해진다.
• 장애물 함수와 이산 볼록성을 사용하여 보간 오차 및 일致성 오차를 분석하며, $\delta/h$ 스케일 인자에 대한 신중한 제어를 수행한다.

실험 결과

연구 질문

• RQ1 하올더 또는 소볼레프 정규성을 갖는 몽주-암페르 방정식의 고전적 해에 대해 두 스케일 유한요소법의 수렴 속도는 무엇인가?
• RQ2 유계이지만 가능성이 떨어지는 헤시안을 갖는 점성 해에 대해 이 방법은 어떻게 성능을 발휘하는가?
• RQ3 비정규성 조건 $f > 0$ 를 가정하지 않고도 탈중력 경우 $f \geq 0$ 에서 최적의 수렴 속도를 달성할 수 있는가?
• RQ4 왜 이론적 오차 추정은 이전 연구에서 관측된 선형 수렴 속도와 격차를 보이는가?
• RQ5 $\delta/h$ 비율이 수렴 속도를 제한하는 데 어떤 역할을 하는가? 이를 완화할 수 있는가?

주요 결과

• 하올더 정규성을 갖는 고전적 해에 대해, 이 방법은 $h^{2/\beta}$ 의 순서로 $L^\infty$ 수렴 속도를 달성한다. 여기서 $\beta > 2$ 는 정규성 지수에 따라 달라진다.
• 소볼레프 정규성을 갖는 고전적 해에 대해서도, 더 약한 정규성 조건 하에서 동일한 수렴 속도 $h^{2/\beta}$ 를 확보한다.
• 조각별 소볼레프 정규성과 불연속 헤시안을 갖는 점성 해에 대해서는 수렴 속도가 감소하지만 여전히 정량적으로 유계이다.
• 탈중력 경우 $f \geq 0$ 에서는 정규화 절차를 통해 $h^{\frac{1}{d}(1 - 2/\beta)}$ 의 순서로 수렴 속도를 도출한다. 이는 최적성이 아니지만 문헌상 처음으로 제시된 추정이다.
• 이론적 비최적성은 이산 알렉산드로프 접근법에서 필수적인 $\delta/h \gg 1$ 인자에 기인하며, 이는 연속성 의존 추정에 기인한다.
• 적절한 이산화 조건 하에서 두 스케일 방법이 펜그 및 점슨의 준라그랑주 방법과 동일시되며, 이로 인해 오차 추정이 해당 방법에도 확장된다.