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
Y1 - 2019/12
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 - Journal article
AN - SCOPUS:85074768798
SN - 0885-7474
VL - 81
SP - 2271
EP - 2302
JO - Journal of Scientific Computing
JF - Journal of Scientific Computing
IS - 3
ER -