Numerical solution of the Dirichlet problem for a Pucci equation in dimension two. Application to homogenization

Abstract The main goal of this article is two fold: (i) To discuss a methodology for the numerical solution of the Dirichlet problem for a Pucci equation in dimension two. (ii) Use the ensuing algorithms to investigate the homogenization properties of the solutions when a coefficient in the Pucci equation oscillates periodically or randomly in space. The solution methodology relies on the combination of a least-squares formulation of the Pucci equation in an appropriate Hilbert space with operator-splitting techniques and mixed finite element approximations. The results of numerical experiments suggest second order accuracy when globally continuous piecewise affine space approximations are used; they also show that the solution of the problem under consideration can be reduced to a sequence of discrete Poisson–Dirichlet problems coupled with one-dimensional optimization problems (one per grid point).

[1]  R. Glowinski,et al.  Numerical methods for the first biharmonic equation and for the two-dimensional Stokes problem , 1977 .

[2]  Finite-element Approximations and Iterative Solutions of a Fourth-order Elliptic Variational Inequality , 1984 .

[3]  Claudianor O. Alves,et al.  Nonvariational elliptic systems , 2002 .

[4]  R. Glowinski.  Finite element methods for incompressible viscous flow , 2003 .

[5]  H. Keller,et al.  Continuation-Conjugate Gradient Methods for the Least Squares Solution of Nonlinear Boundary Value Problems , 1985 .

[6]  R. Glowinski.  Lectures on Numerical Methods for Non-Linear Variational Problems , 1981 .

[7]  L. Reinhart,et al.  On the numerical analysis of the Von Karman equations: Mixed finite element approximation and continuation techniques , 1982 .

[8]  L. Evans.  Periodic homogenisation of certain fully nonlinear partial differential equations , 1992, Proceedings of the Royal Society of Edinburgh: Section A Mathematics.

[9]  R. Glowinski,et al.  Numerical Analysis of Variational Inequalities , 1981 .

[10]  R. Glowinski,et al.  Augmented Lagrangian and Operator-Splitting Methods in Nonlinear Mechanics , 1987 .

[11]  P. G. Ciarlet,et al.  Basic error estimates for elliptic problems , 1991 .

[12]  Roland Glowinski,et al.  On the numerical solution of a two-dimensional Pucci's equation with dirichlet boundary conditions : a least-squares approach , 2005 .

[13]  Panagiotis E. Souganidis,et al.  Homogenization of fully nonlinear, uniformly elliptic and parabolic partial differential equations in stationary ergodic media , 2005 .

[14]  Roland Glowinski,et al.  An approximate factorization/least squares solution method for a mixed finite element approximation of the Cahn-Hilliard equation , 1996 .

[15]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[16]  D. Anderson,et al.  Algorithms for minimization without derivatives , 1974 .

[17]  M. Fortin,et al.  Augmented Lagrangian methods : applications to the numerical solution of boundary-value problems , 1983 .

[18]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[19]  L. Caffarelli,et al.  Fully Nonlinear Elliptic Equations , 1995 .

[20]  Roland Glowinski,et al.  Applications of operator-splitting methods to the direct numerical simulation of particulate and free-surface flows and to the numerical solution of the two-dimensional elliptic Monge-Ampère equation , 2008 .

[21]  Roland Glowinski,et al.  Numerical solution of the two-dimensional elliptic Monge–Ampère equation with Dirichlet boundary conditions: a least-squares approach , 2004 .