TY - JOUR

T1 - An alternating direction method of multipliers for the numerical solution of a fully nonlinear partial differential equation involving the Jacobian determinant

AU - Caboussat, Alexandre

AU - GLOWINSKI, Roland

N1 - Funding Information:
The authors thank Prof. B. Dacorogna (EPFL) for suggesting the investigation of this problem and some test problems of interest, Prof. M. Picasso (EPFL) for helpful comments and discussions, and the anonymous referees for constructive comments. This work was partially supported by the Swiss National Science Foundation (grant 165785) and the National Science Foundation (grant NSF DMS-0913982).
Funding Information:
∗Submitted to the journal’s Methods and Algorithms for Scientific Computing section October 11, 2016; accepted for publication (in revised form) October 17, 2017; published electronically January 5, 2018. http://www.siam.org/journals/sisc/40-1/M109407.html Funding: This work was partially supported by the Swiss National Science Foundation (grant 165785) and the National Science Foundation (grant NSF DMS-0913982). †Haute Ecole de Gestion de Genève (Geneva School of Business Administration), University of Applied Sciences Western Switzerland (HES-SO), 1227 Carouge, Switzerland (alexandre. caboussat@hesge.ch, http://campus.hesge.ch/caboussata/). ‡Department of Mathematics, University of Houston, Houston, TX 77204-3008 (roland@math. uh.edu, http://math.uh.edu/∼roland/), and Hong-Kong Baptist University, Kowloon, Hong-Kong.

PY - 2018

Y1 - 2018

N2 - We consider the Dirichlet problem for a partial differential equation involving the Jacobian determinant in two dimensions of space. The problem consists in finding a vector-valued function such that the determinant of its gradient is given pointwise in a bounded domain, together with essential boundary conditions. This problem was initially considered in Dacorogna and Moser [Ann. Inst. H. Poincaré Anal. Non Linéaire, 7 (1990), pp. 1-26], and several theoretical generalizations have been derived since. In this work, we design a numerical algorithm for the approximation of the solution of such a problem for various kinds of boundary data. The proposed method relies on an augmented Lagrangian algorithm with biharmonic regularization, and low order mixed finite element approximations. An iterative method allows us to decouple the nonlinearity and the differential operators. Numerical experiments show the capabilities of the method for benchmarks and then for more demanding test problems.

AB - We consider the Dirichlet problem for a partial differential equation involving the Jacobian determinant in two dimensions of space. The problem consists in finding a vector-valued function such that the determinant of its gradient is given pointwise in a bounded domain, together with essential boundary conditions. This problem was initially considered in Dacorogna and Moser [Ann. Inst. H. Poincaré Anal. Non Linéaire, 7 (1990), pp. 1-26], and several theoretical generalizations have been derived since. In this work, we design a numerical algorithm for the approximation of the solution of such a problem for various kinds of boundary data. The proposed method relies on an augmented Lagrangian algorithm with biharmonic regularization, and low order mixed finite element approximations. An iterative method allows us to decouple the nonlinearity and the differential operators. Numerical experiments show the capabilities of the method for benchmarks and then for more demanding test problems.

KW - ADMM algorithm

KW - Augmented Lagrangian methods

KW - Biharmonic regularization

KW - Finite element method

KW - Jacobian determinant

KW - Quadratically constrained minimization

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

U2 - 10.1137/16M1094075

DO - 10.1137/16M1094075

M3 - Article

AN - SCOPUS:85044535359

VL - 40

SP - A52-A80

JO - SIAM Journal of Scientific Computing

JF - SIAM Journal of Scientific Computing

SN - 1064-8275

IS - 1

ER -