Mixed Finite Element Methods for the Fully Nonlinear Monge–Ampère Equation Based on the Vanishing Moment Method

2009/01/01 by Xiaobing Feng, Michael Neilan · 12 citations
Mathematics · Engineering · #Nonlinear Partial Differential Equations #Geometric Analysis and Curvature Flows #Advanced Numerical Methods in Computational Mathematics

paper · doi:10.1137/070710378

Abstract

This paper studies mixed finite element approximations of the viscosity solution to the Dirichlet problem for the fully nonlinear Monge–Ampère equation det(D²u⁰)=f (>0) based on the vanishing moment method which was proposed recently by the authors in [X. Feng and M. Neilan, J. Scient. Comp., DOI 10.1007/s10915-008-9221-9, 2008]. In this approach, the second-order fully nonlinear Monge–Ampère equation is approximated by the fourth order quasilinear equation -εΔ² u^ε + detD²u^ε=f. It was proved in [X. Feng, Trans. AMS, submitted] that the solution u^ε converges to the unique convex viscosity solution u⁰ of the Dirichlet problem for the Monge–Ampère equation. This result then opens a door for constructing convergent finite element methods for the fully nonlinear second-order equations, a task which has been impracticable before. The goal of this paper is threefold. First, we develop a family of Hermann–Miyoshi-type mixed finite element methods for approximating the solution u^ε of the regularized fourth-order problem, which computes simultaneously u^ε and the moment tensor σ^ε:=D²u^ε. Second, we derive error estimates, which track explicitly the dependence of the error constants on the parameter ε, for the errors u^ε-u^ε_h and σ⁰-σ_h^ε. Finally, we present a detailed numerical study on the rates of convergence in terms of powers of ε for the error u⁰-u_h^ε and σ^ε-σ_h^ε, and numerically examine what is the “best” mesh size h in relation to ε in order to achieve these rates. Due to the strong nonlinearity of the underlying equation, the standard perturbation argument for error analysis of finite element approximations of nonlinear problems does not work for the problem. To overcome the difficulty, we employ a fixed point technique which strongly relies on the stability of the linearized problem and its mixed finite element approximations.

Citations

• A computational fluid mechanics solution to the Monge-Kantorovich mass transfer problem

Cited by

• Convergent Semi-Lagrangian Methods for the Monge--Ampère Equation on Unstructured Grids
• Fully discrete least-squares splitting scheme for the Monge-Ampère equation: finite element analysis and convergence
• Discrete ABP Estimate and Convergence Rates for Linear Elliptic Equations in Non-divergence Form
• Analysis of a finite element method for second order uniformly elliptic PDEs in non-divergence form
• Low-rank cross approximation of function-valued tensors for reduced-order modeling of parametric PDEs
• Discontinuous Galerkin methods for nonvariational problems
• The Vanishing Moment Method for Fully Nonlinear Second Order Partial Differential Equations: Formulation, Theory, and Numerical Analysis
• Fast Operator-Splitting Methods for Nonlinear Elliptic Equations
• Robust and fast iterative method for the elliptic Monge-Ampère equation
• A sequential least squares method for elliptic equations in non-divergence form
• A Narrow-stencil finite difference method for approximating viscosity solutions of fully nonlinear elliptic partial differential equations with applications to Hamilton-Jacobi-Bellman equations
• Two Numerical Methods for the elliptic Monge-Ampère equation

Related

• Mixed finite element methods for the fully nonlinear Monge-Ampère equation based on the vanishing moment method
• Galerkin Methods for the Fully Nonlinear Monge-Ampère Equation
• The Vanishing Moment Method for Fully Nonlinear Second Order Partial Differential Equations: Formulation, Theory, and Numerical Analysis
• Standard finite elements for the numerical resolution of the elliptic Monge-Ampere equation: mixed methods
• Two-Scale Method for the Monge-Ampère Equation: Convergence to the Viscosity Solution