We construct a nonlinear least-squares finite element method for computing the smooth convex solutions of the Dirichlet boundary value problem of the Monge-Ampère equation on strictly convex smooth domains in
- NSF-PAR ID:
- 10260020
- Date Published:
- Journal Name:
- Numerische Mathematik
- ISSN:
- 0029-599X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
. It is based on an isoparametric finite element space with exotic degrees of freedom that can enforce the convexity of the approximate solutions. A priori anda posteriori error estimates together with corroborating numerical results are presented. -
We consider the numerical solution of a fourth‐order total variation flow problem representing surface relaxation below the roughening temperature. Based on a regularization and scaling of the nonlinear fourth‐order parabolic equation, we perform an implicit discretization in time and a C0Interior Penalty Discontinuous Galerkin (C0IPDG) discretization in space. The C0IPDG approximation can be derived from a mixed formulation involving numerical flux functions where an appropriate choice of the flux functions allows to eliminate the discrete dual variable. The fully discrete problem can be interpreted as a parameter dependent nonlinear system with the discrete time as a parameter. It is solved by a predictor corrector continuation strategy featuring an adaptive choice of the time step sizes. A documentation of numerical results is provided illustrating the performance of the C0IPDG method and the predictor corrector continuation strategy. The existence and uniqueness of a solution of the C0IPDG method will be shown in the second part of this paper.
-
We prove the existence and uniqueness of a solution of a C0Interior Penalty Discontinuous Galerkin (C0IPDG) method for the numerical solution of a fourth‐order total variation flow problem that has been developed in part I of the paper. The proof relies on a nonlinear version of the Lax‐Milgram Lemma. It requires to establish that the nonlinear operator associated with the C0IPDG approximation is Lipschitz continuous and strongly monotone on bounded sets of the underlying finite element space.
-
In this paper, the convergence of a homotopy method (1.1) for solving the steady state problem of Burgers’ equation is considered. When ν is fixed, we prove that the solution of (1.1) converges to the unique steady state solution as epsilon → 0, which is independent of the initial conditions. Numerical examples are presented to confirm this conclusion by using the continuous finite element method. In contrast, when ν = epsilon → 0, numerically we show that steady state solutions obtained by (1.1) indeed depend on initial conditions.more » « less
-
We design and analyze a C0 interior penalty method for the approximation of classical solutions of the Dirichlet boundary value problem of the Monge–Ampère equation on convex polygonal domains. The method is based on an enhanced cubic Lagrange finite element that enables the enforcement of the convexity of the approximate solutions. Numerical results that corroborate the a priori and a posteriori error estimates are presented. It is also observed from numerical experiments that this method can capture certain weak solutions.more » « less