An adaptive method for the numerical solution of a 2D Monge-Ampère equation

Caboussat, Alexandre (Haute école de gestion de Genève, HES-SO // Haute Ecole Spécialisée de Suisse Occidentale) ; Gourzoulidis, Dimitrios (Haute école de gestion de Genève, HES-SO // Haute Ecole Spécialisée de Suisse Occidentale ; EcolePolytechnique Fédérale de Lausanne, Switzerland) ; Picasso, Marco (EcolePolytechnique Fédérale de Lausanne, Switzerland)

Parabolic fully nonlinear equations may be found in various applications,for instance in optimal portfolio management strategy. We focus here on a canonical parabolic Monge-Ampère equationin two space dimensions. A numerical method has beeninvestigated in[1]. The goal is toextend themethodology by coupling atime steppingsemi-implicit methodthat relies on a conservative formulation of the Monge-Ampère equationwith mesh adaptation.The parabolic Monge-Ampère equation can be expressed as astronglynonlinear, heat-type, parabolic equation, where the nonlinear diffusion function is expressed asa function ofthe cofactor matrix of the Hessianmatrixof the solution. We linearize this diffusion operator and advocatea semi-implicit time-stepping algorithm.In particular, we use the time-evolutive equation to reach a stationarysolutioncorresponding to a solutionof the elliptic Monge-Ampère equation. A loworder, piecewise linear,finite element method is used for spacediscretization, together with a mixed approach for the approximation of the second derivatives. The error is bounded above by an error indicator plus an extra term that can be disregarded in special cases. A mesh adaptivity strategy based on these estimates is then implemented within thetime-stepping algorithmfor the nonlinear equation. Numerical experiments exhibit appropriate convergence orders and arobust behavior. Adaptive mesh refinement proves to be efficient and accurate to tackle test cases with singularities. Inparticular, we consider equations with exact solutionswith singularitieson the boundary of the domain, or with right-hand sides involving Dirac functions.

Note: Due to the COVID-19 outbreak, the 10th International Conference on Adaptative Modeling and Simulation venue in Gothenburg was cancelled. The proceedings of the online conference are however published according to the original schedule.

Conference Type:
published full paper
Economie et Services
HEG - Genève
CRAG - Centre de Recherche Appliquée en Gestion
Gothenburg, Sweden, 21-23 June 2021
Gothenburg, Sweden
21-23 June 2021
8 p.
Published in:
Proceedings of the 10th International Conference on Adaptive Modeling and Simulation (ADMOS 2021)
Appears in Collection:

 Record created 2021-07-12, last modified 2021-07-12

Download fulltext

Rate this document:

Rate this document:
(Not yet reviewed)