[논문 리뷰] Finite element approximation of the Isaacs equation

Abner J. Salgado, Wujun Zhang|arXiv (Cornell University)|2015. 12. 30.

Advanced Mathematical Modeling in Engineering참고 문헌 13인용 수 3

한 줄 요약

이 논문은 완전 비선형 타원형 PDE인 이자크스 방정식을 위한 이중 척도 유한요소법을 제안한다. 이는 체계적 크기 $h$에서의 유한요소 이산화와 체계적 크기 $\varepsilon$에서의 적분형 미분 근사(구적법 기반)를 조합함으로써 이루어진다. 이 방법은 약간 예각인 메쉬에서 이산 최대원리(discrete maximum principle)를 만족하며, $C^{2,s}$ 정규성 조건 하에서 $\mathcal{O}(\varepsilon^s + h/\varepsilon^2 |\log h|)$ 수렴 속도로 점근 해에 수렴한다. 이는 비볼록 비선형 PDE에 대해 새로운 이중 척도 FEM 프레임워크를 통해 수렴성과 수렴 속도를 확립한 것이다.

ABSTRACT

We propose and analyze a two-scale finite element method for the Isaacs equation. The fine scale is given by the mesh size $h$ whereas the coarse scale $\varepsilon$ is dictated by an integro-differential approximation of the partial differential equation. We show that the method satisfies the discrete maximum principle provided that the mesh is weakly acute. This, in conjunction with weak operator consistency of the finite element method, allows us to establish convergence of the numerical solution to the viscosity solution as $\varepsilon, h o0$, and $\varepsilon \gtrsim h^{1/2}|\log h|$. In addition, using a discrete Alexandrov Bakelman Pucci estimate we deduce rates of convergence, under suitable smoothness assumptions on the exact solution.

연구 동기 및 목표

• 클래식한 변분 방법이 점근 해 이론에 의해 실패하는 비선형 타원형 PDE, 특히 이자크스 방정식에 대해 수렴하는 유한요소법을 개발하는 것.
• 완전 비선형 PDE에 대한 강력한 유한요소법의 부족을 보완하기 위해 메쉬 기반 FEM과 비국소 적분형 미분 근사의 조합을 통한 이중 척도 접근법을 도입하는 것.
• 이산 최대원리와 약한 연산자 일致성(weak operator consistency)을 활용하여 수렴성과 수렴 속도를 확립하는 것.
• 기존의 볼록/볼록이 아닌 방정식에 대한 수렴 결과를 이자크스 방정식을 모형 문제로 삼아 비볼록/비볼록이 아닌 경우로 확장하는 것.
• 수렴 속도를 유지하면서 효율적인 계산을 가능하게 하는 실용적인, 구적법 기반의 구현을 제공하는 것.

제안 방법

• 이 방법은 이중 척도 분해를 사용한다: 메쉬 크기 $h$에서의 미세 척도 FEM과 체계적 크기 $\varepsilon$에서의 비국소 적분형 미분 연산자 근사. 이는 카프레리와 실베스트레의 헤시안 항 근사에서 유도된다.
• 이자크스 연산자는 라플라스 연산자와 중심 두 번째 차분 연산자 $\mathfrak{d}u(x,y) = u(x+y) - 2u(x) + u(x-y)$를 포함하는 비국소 항의 합으로 재작성된다.
• 비국소 항은 컴팩트 지지, 반지름 대칭 커널 $\varphi$ 위에서의 구적 공식을 사용하여 근사되며, 이는 이차 다항식에 대해 정확한 이산 연산자 $I_{\varepsilon,h}[w](x)$를 유도한다.
• 이산 스킴은 모든 노드 $z \in \mathscr{N}_h$ 에 대해 $\bar{\mathfrak{L}}_h^\varepsilon[\tilde{v}_h^\varepsilon](z) = \frac{\lambda}{2}\Delta_h \tilde{v}_h^\varepsilon(z) + I_{\varepsilon,h}[\tilde{v}_h^\varepsilon](z) = f_z$ 로 정의되며, 일致성과 단조성 보장.
• 메쉬가 약간 예각이면 이 방법은 이산 최대원리를 만족하며, 이는 점근 해로의 수렴에 필수적이다.
• 구현은 효율적이다: 시스템 행렬은 사전에 가중치 $\omega_j \varphi(\xi_j)$ 를 계산하고, 기저 함수를 $z + \varepsilon M(z)\xi_j$ 에서만 그 지지 내에서 평가함으로써 구성된다.

실험 결과

연구 질문

• RQ1 표준 FEM 도구만을 사용하여 비볼록 완전 비선형 PDE인 이자크스 방정식의 점근 해로 수렴하는 이중 척도 유한요소법을 설계할 수 있는가?
• RQ2 제안된 방법이 약간 예각인 메쉬에서 이산 최대원리를 유지하는가? 이는 점근 해 이론에 의한 수렴을 가능하게 하는가?
• RQ3 정확한 해가 $C^{2,s}(\bar{\Omega})$ 에 속할 경우, 이 방법의 수렴 속도는 무엇인가?
• RQ4 비국소 연산자에서의 구적 오차는 전체 수렴 속도에 어떤 영향을 미치는가?
• RQ5 수렴성과 속도 추정을 유지하면서 효율적인 어셈블리가 가능한 실용적인 구현이 가능한가?

주요 결과

• 약간 예각인 메쉬에서 이 방법은 이산 최대원리를 만족하며, 이는 점근 해로의 수렴을 보장하는 데 필수적이다.
• 메쉬 크기 $h$와 체계적 크기 $\varepsilon$ 가 모두 0으로 수렴할 때 수렴성이 확립되며, 조건 $\varepsilon \gtrsim h^{1/2} |\log h|$ 가 필요하여 안정성과 일치성 확보.
• $C^{2,s}$ 정규성 조건 하에서 이 방법은 $L^\infty$ 노름에서 $\mathcal{O}(\varepsilon^s + h/\varepsilon^2 |\log h|)$ 의 수렴 속도를 달성한다.
• 비국소 연산자의 구적 오차는 $C^{2,s}$ 함수에 대해 $\mathcal{O}(\varepsilon^s)$ 로 유계이며, 이는 $x + B_\varepsilon$ 위에서 국소 이차 근사에 기인한다.
• 이산 알렉산드로프-바켈레프-푸크 추정식을 사용하여 수렴 속도를 유도하였으며, 이는 최대원리와 오차 통제를 연결한다.
• 구현은 계산적으로 가능하다: 기저 함수의 평가는 $z + \varepsilon M(z)\xi_j$ 에서만 국소적으로 수행되며, 가중치는 사전 계산되어 있어 효율적인 행렬 어셈블리가 가능하다.