Issue |
ESAIM: COCV
Volume 19, Number 3, July-September 2013
|
|
---|---|---|
Page(s) | 780 - 810 | |
DOI | https://doi.org/10.1051/cocv/2012033 | |
Published online | 03 June 2013 |
A least-squares method for the numerical solution of the Dirichlet problem for the elliptic monge − ampère equation in dimension two∗
1
Haute École de Gestion/Geneva School of Business
Administration, Genève, Switzerland
alexandre.caboussat@hesge.ch
2
University of Houston, Department of Mathematics,
4800 Calhoun Rd, Houston, 77204-3008
Texas,
USA
caboussat@math.uh.edu
3
University of Houston, Department of Mathematics,
4800 Calhoun Rd, Houston, 77204-3008
Texas,
USA
roland@math.uh.edu
4
Rice University, Department of Computational and Applied
Mathematics, MS
134, Houston,
77251-1892
Texas,
USA
sorensen@rice.edu
Received:
31
January
2011
Revised:
21
August
2012
We address in this article the computation of the convex solutions of the Dirichlet problem for the real elliptic Monge − Ampère equation for general convex domains in two dimensions. The method we discuss combines a least-squares formulation with a relaxation method. This approach leads to a sequence of Poisson − Dirichlet problems and another sequence of low dimensional algebraic eigenvalue problems of a new type. Mixed finite element approximations with a smoothing procedure are used for the computer implementation of our least-squares/relaxation methodology. Domains with curved boundaries are easily accommodated. Numerical experiments show the convergence of the computed solutions to their continuous counterparts when such solutions exist. On the other hand, when classical solutions do not exist, our methodology produces solutions in a least-squares sense.
Résumé
Nous étudions, dans cet article, une méthode numérique, pour le calcul des solutions convexes du problème de Dirichlet pour l’équation de Monge − Ampère elliptique, dans des domaines bi-dimensionnel convexes. Une méthode de moindres carrés est couplée à un algorithme de relaxation, conduisant à la résolution d’une suite de problèmes de Poisson − Dirichlet, et d’une suite de problèmes de valeurs propres de petite dimension d’un type nouveau. Une approximation par éléments finis mixtes, couplée à une méthode de régularisation, est utilisée pour implémenter la méthode de moindres-carrés/relaxation ci-dessus, de sorte que les domaines avec frontière courbe sont traités facilement. Des expériences numériques montrent la convergence des solutions calculées vers la solution convexe du problème continu, lorsqu’une telle solution existe. Par ailleurs, si le problème n’a pas de solution classique, notre méthodologie fournit des solutions au sens des moindres carrés.
Mathematics Subject Classification: 65N30 / 65K10 / 65F30 / 49M15 / 49K20
Key words: Monge − Ampère equation / least-squares method / biharmonic problem / conjugate gradient method / quadratic constraint minimization / mixed finite element methods
© EDP Sciences, SMAI, 2013
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.