New $C^0$ interior penalty method for Monge-Ampère equations

Tianyang Chu,Hailong Guo,Zhimin Zhang
2024-08-31
Abstract:Monge-Ampère equation is a prototype second-order fully nonlinear partial differential equation. In this paper, we propose a new idea to design and analyze the $C^0$ interior penalty method to approximation the viscosity solution of the Monge-Ampère equation. The new methods is inspired from the discrete Miranda-Talenti estimate. Based on the vanishing moment representation, we approximate the Monge-Ampère equation by the fourth order semi-linear equation with some additional boundary conditions. We will use the discrete Miranda-Talenti estimates to ensure the well-posedness of the numerical scheme and derive the error estimates.
Numerical Analysis
What problem does this paper attempt to address?
### Problems the paper attempts to solve This paper aims to design and analyze a new C0 interior penalty method to approximately solve the viscosity solutions of the Monge - Ampère equation. The Monge - Ampère equation is a typical second - order fully nonlinear partial differential equation, which appears in many important applications, such as differential geometry, reflector design and optimal transport problems. Due to the fully nonlinear nature of the Monge - Ampère equation, the traditional weak solution theory based on the variational method is not directly applicable. Therefore, new numerical methods need to be proposed to solve this problem. Specifically, the paper makes the following innovations: 1. **New C0 interior penalty method**: Design and analyze a new C0 interior penalty method through discrete Miranda - Talenti estimates. This method uses standard Lagrange elements, which are flexible and easy to implement, and are especially suitable for high - order PDE solvers. 2. **Vanishing moment method**: Use the vanishing moment method to approximate the Monge - Ampère equation as a fourth - order semilinear equation with additional boundary conditions. This method not only simplifies the problem handling but also ensures the goodness of the numerical scheme and error estimates. 3. **Calculation of viscosity solutions**: The new method aims to calculate the viscosity solutions of the Monge - Ampère equation, which requires lower regularity conditions and thus can handle singular solutions. 4. **Error analysis**: Prove the goodness of the numerical method and error estimates through discrete Miranda - Talenti estimates and fixed - point techniques. ### Main contributions - **Theoretical contributions**: Propose a new C0 interior penalty method and prove its goodness and error estimates through strict mathematical analysis. - **Numerical experiments**: Verify the effectiveness of the new method and the correctness of the theoretical results through a series of numerical experiments. ### Numerical experiments The paper verifies the performance of the new method and the theoretical results through multiple numerical experiments. These experiments include tests in two - dimensional and three - dimensional cases, as follows: 1. **Two - dimensional numerical experiments**: - **Numerical example I**: Study the approximation of the perturbed equation (1.2) to the fully nonlinear Monge - Ampère equation (1.1). The results show that as the parameter ϵ decreases, the solution uϵ converges to the exact solution u at a rate of O(ϵ^0.25) in the H2 norm. - **Numerical example II**: Analyze the error convergence rate between uϵ and uϵh when ϵ = 0.01 is fixed. The results show that for quadratic and cubic elements, the optimal O(h^(k - 1)) convergence rates are observed respectively. - **Numerical example III**: Test the ability of the new method to calculate viscosity solutions. The results show that the new method can effectively track viscosity convex solutions. 2. **Three - dimensional numerical experiments**: - **Numerical example IV**: Study the approximation of the perturbed equation (1.2) to the fully nonlinear Monge - Ampère equation (1.1). The results show that, similar to the two - dimensional case, the solution uϵ converges to the exact solution u at a rate of O(ϵ^0.25) in the H2 norm. In summary, by proposing a new C0 interior penalty method, this paper successfully solves the numerical solution problem of the Monge - Ampère equation and verifies the effectiveness of the method through theoretical analysis and numerical experiments.