TY - JOUR

T1 - A Finite Element/Operator-Splitting Method for the Numerical Solution of the Three Dimensional Monge–Ampère Equation

AU - Liu, Hao

AU - GLOWINSKI, Roland

AU - Leung, Shingyu

AU - Qian, Jianliang

N1 - Funding Information:
The work of R. Glowinski was partially supported by the Hong Kong Kennedy Wong foundation. The work of S. Leung was partially supported by the Hong Kong RGC Grants 16302819 and 16309316. The work of J. Qian was partially supported by NSF.

PY - 2019/12/1

Y1 - 2019/12/1

N2 - In the present article we extend to the three-dimensional elliptic Monge–Ampère equation the method discussed in Glowinski et al. (J Sci Comput 79:1–47, 2019) for the numerical solution of its two-dimensional variant. As in Glowinski et al. (2019) we take advantage of an equivalent divergence formulation of the Monge–Ampère equation, involving the cofactor matrix of the Hessian of the solution. We associate with the above divergence formulation an initial value problem, well suited to time discretization by operator splitting and space approximation by low order mixed finite element methods. An important ingredient of our methodology is forcing the positive semi-definiteness of the approximate Hessian by a hard thresholding eigenvalue projection. The resulting method is robust and easy to implement. It can handle problems with smooth and non-smooth solutions on domains with curved boundary. Using piecewise affine approximations for the solution and its six second-order derivatives, one can achieve second-order convergence rates for problems with smooth solutions.

AB - In the present article we extend to the three-dimensional elliptic Monge–Ampère equation the method discussed in Glowinski et al. (J Sci Comput 79:1–47, 2019) for the numerical solution of its two-dimensional variant. As in Glowinski et al. (2019) we take advantage of an equivalent divergence formulation of the Monge–Ampère equation, involving the cofactor matrix of the Hessian of the solution. We associate with the above divergence formulation an initial value problem, well suited to time discretization by operator splitting and space approximation by low order mixed finite element methods. An important ingredient of our methodology is forcing the positive semi-definiteness of the approximate Hessian by a hard thresholding eigenvalue projection. The resulting method is robust and easy to implement. It can handle problems with smooth and non-smooth solutions on domains with curved boundary. Using piecewise affine approximations for the solution and its six second-order derivatives, one can achieve second-order convergence rates for problems with smooth solutions.

UR - http://www.scopus.com/inward/record.url?scp=85074768798&partnerID=8YFLogxK

U2 - 10.1007/s10915-019-01080-4

DO - 10.1007/s10915-019-01080-4

M3 - Article

AN - SCOPUS:85074768798

VL - 81

SP - 2271

EP - 2302

JO - Journal of Scientific Computing

JF - Journal of Scientific Computing

SN - 0885-7474

IS - 3

ER -