Abner J. Salgado,Wujun Zhang
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\to0$, 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.
What problem does this paper attempt to address?
This paper attempts to solve the finite - element approximation problem of Isaacs equations. Specifically, the author proposes and analyzes a two - scale finite - element method for Isaacs equations. The following is a detailed interpretation of this problem:
### Research Background
Fully nonlinear elliptic partial differential equations (PDEs) are widely used in differential geometry, optimal mass transportation, stochastic optimal control, etc. However, the numerical solution methods for such equations are still in the development stage, especially when the finite - element method (FEM) is applied. For fully nonlinear PDEs, the correct concept of solution is the so - called viscosity solution, which is based on the maximum principle rather than the variational principle.
### Research Objectives
The main objective of this paper is to extend the existing two - scale finite - element method for linear elliptic PDEs to fully nonlinear elliptic PDEs, especially Isaacs equations:
\[ F[u](x) := \inf_{\alpha \in A} \sup_{\beta \in B} \left[ A_{\alpha, \beta}(x) : D^2u(x) \right] = f(x) \quad \text{in } \Omega, \quad u = 0 \text{ on } \partial \Omega, \]
where \(\Omega \subset \mathbb{R}^d\) is an open, bounded convex domain, \(f \in C^{0,1}(\overline{\Omega})\), the sets \(A\) and \(B\) are arbitrary finite sets, and the matrix \(A_{\alpha, \beta} \in C^{0,1}(\overline{\Omega}, GL_d(\mathbb{R}))\) is symmetric and uniformly elliptic.
### Main Contributions
1. **Two - scale finite - element method**: The author introduced a two - scale finite - element method. The fine scale is given by the mesh size \(h\), while the coarse scale \(\varepsilon\) is determined by the integro - differential approximation.
2. **Discrete maximum principle**: It is proved that this method satisfies the discrete maximum principle under weakly acute - angled meshes.
3. **Convergence**: Through weak operator consistency, it is proved that the numerical solution converges to the viscosity solution under the conditions that \(\varepsilon, h \to 0\) and \(\varepsilon \gtrsim h^{1/2}|\log h|\).
4. **Convergence rate**: Using the discrete Alexandrov - Bakelman - Pucci estimate, the convergence rate under appropriate smoothness assumptions is derived.
### Summary of Mathematical Formulas
- The form of Isaacs equations is:
\[
F[u](x) := \inf_{\alpha \in A} \sup_{\beta \in B} \left[ A_{\alpha, \beta}(x) : D^2u(x) \right] = f(x) \quad \text{in } \Omega, \quad u = 0 \text{ on } \partial \Omega,
\]
where \(D^2u(x)\) represents the Hessian matrix of \(u(x)\).
- The integro - differential operator \(I_{\alpha, \beta}^\varepsilon[u]\) is defined as:
\[
I_{\alpha, \beta}^\varepsilon[u](x) = \frac{1}{\varepsilon^{d + 2}} \det M_{\alpha, \beta}(x) \int_{\mathbb{R}^d} d_u(x, y) \phi\left( M_{\alpha, \beta}(x)^{-1} \frac{y}{\varepsilon} \right) dy,
\]
where \(d_u(x, y) = u\)