Alexandre Caboussat, Dimitrios Gourzoulidis, Marco Picasso
{"title":"An adaptive least-squares algorithm for the elliptic Monge–Ampère equation","authors":"Alexandre Caboussat, Dimitrios Gourzoulidis, Marco Picasso","doi":"10.5802/crmeca.222","DOIUrl":null,"url":null,"abstract":"We address the numerical solution of the Dirichlet problem for the two-dimensional elliptic Monge–Ampère equation using a least-squares/relaxation approach. The relaxation algorithm allows the decoupling of the differential operators from the nonlinearities of the equation, within a splitting approach. The approximation relies on mixed low order finite element methods with regularization techniques. In order to account for data singularities in non-smooth cases, we introduce an adaptive mesh refinement technique. The error indicator is based an independent formulation of the Monge–Ampère equation under divergence form, which allows to explicit a residual term. We show that the error is bounded from above by an a posteriori error indicator plus an extra term that remains to be estimated. This indicator is then used within the existing least-squares framework. The results of numerical experiments support the convergence of our relaxation method to a convex classical solution, if such a solution exists. Otherwise they support convergence to a generalized solution in a least-squares sense. Adaptive mesh refinement proves to be efficient, robust, and accurate to tackle test cases with singularities.","PeriodicalId":10566,"journal":{"name":"Comptes Rendus. Chimie","volume":"54 1","pages":"0"},"PeriodicalIF":1.2000,"publicationDate":"2023-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Comptes Rendus. Chimie","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.5802/crmeca.222","RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"CHEMISTRY, MULTIDISCIPLINARY","Score":null,"Total":0}
引用次数: 0
Abstract
We address the numerical solution of the Dirichlet problem for the two-dimensional elliptic Monge–Ampère equation using a least-squares/relaxation approach. The relaxation algorithm allows the decoupling of the differential operators from the nonlinearities of the equation, within a splitting approach. The approximation relies on mixed low order finite element methods with regularization techniques. In order to account for data singularities in non-smooth cases, we introduce an adaptive mesh refinement technique. The error indicator is based an independent formulation of the Monge–Ampère equation under divergence form, which allows to explicit a residual term. We show that the error is bounded from above by an a posteriori error indicator plus an extra term that remains to be estimated. This indicator is then used within the existing least-squares framework. The results of numerical experiments support the convergence of our relaxation method to a convex classical solution, if such a solution exists. Otherwise they support convergence to a generalized solution in a least-squares sense. Adaptive mesh refinement proves to be efficient, robust, and accurate to tackle test cases with singularities.
期刊介绍:
The Comptes Rendus - Chimie are a free-of-charge, open access and peer-reviewed electronic scientific journal publishing original research articles. It is one of seven journals published by the Académie des sciences.
Its objective is to enable researchers to quickly share their work with the international scientific community.
The Comptes Rendus - Chimie also publish journal articles, thematic issues and articles reflecting the history of the Académie des sciences and its current scientific activity.