A penalty/Newton/conjugate gradient method for the solution of obstacle problems

Glowinski, Roland; Kuznetsov, Yuri A.; Pan, Tsorng-Whay

Comptes Rendus. Mathématique, Tome 336 (2003) no. 5, pp. 435-440.

Motivé par la recherche des solutions non négatives d'un système d'équations eiconales, avec conditions aux limites de Dirichlet, on étudie dans cette Note une méthode pour la résolution numérique de problèmes d'inéquations variationnelles paraboliques pour des ensembles convexes du type $\mathrm{K}=\{\mathrm{v}\mid \mathrm{v}\in {\mathrm{H}}_0^1\left(\Omega \right)$, v⩾ψ p.p. sur $\Omega \}$. La méthode numérique combine pénalité et algorithme de Newton, les problèmes linéarisés étant résolus par un algorithme de gradient conjugué qui demande à chaque iteration la résolution d'un problème linéaire pour un analogue discret de l'opérateur elliptique I−μΔ avec μ>0. Les essais numériques montrent que la méthode ainsi obtenue a de bonnes propriétés de convergence, même pour des petites valeurs du paramètre de pénalité.

Motivated by the search for non-negative solutions of a system of Eikonal equations with Dirichlet boundary conditions, we discuss in this Note a method for the numerical solution of parabolic variational inequality problems for convex sets such as $\mathrm{K}=\{\mathrm{v}\mid \mathrm{v}\in {\mathrm{H}}_0^1\left(\Omega \right)$, v⩾ψ a.e. on $\Omega \}.$ The numerical methodology combines penalty and Newton's method, the linearized problems being solved by a conjugate gradient algorithm requiring at each iteration the solution of a linear problem for a discrete analogue of the elliptic operator I−μΔ. Numerical experiments show that the resulting method has good convergence properties, even for small values of the penalty parameter.

[1] Carlsson, H.; Glowinski, R. Vibrations of Euler–Bernoulli beams with pointwise obstacles (Gatignol, R.; Soubbaramayer, eds.), Advances in Kinetic Theory and Continuum Mechanics, Springer-Verlag, Berlin, 1991, pp. 261-275

[2] Ciarlet, P.G. The Finite Element Method for Elliptic Problems, North-Holland, Amsterdam, 1978

[3] Ciarlet, P.G. Basic error estimates for elliptic problems (Ciarlet, P.G.; Lions, J.-L., eds.), Handbook of Numerical Analysis, Vol. II, North-Holland, Amsterdam, 1991, pp. 17-352

[4] Dacorogna, B.; Glowinski, R.; Pan, T.-W. Numerical solution of a system of Eikonal equations, C. R. Acad. Sci. Paris, Sér. I, Volume 336 (2003)

[5] Dean, E.J.; Glowinski, R.; Kuo, Y.M.; Nasser, M.G. On the discretization of some second order in time differential equations. Applications to nonlinear wave problems (Balakrishnan, A.V., ed.), Computational Techniques in Identification and Control of Flexible Flight Structures, Optimization Software Inc., Los Angeles, 1990, pp. 199-246

[6] Duvaut, G.; Lions, J.-L. Inequalities in Mechanics and Physics, Springer-Verlag, Berlin, 1976

[7] Glowinski, R. Numerical Methods for Nonlinear Variational Problems, Springer-Verlag, New York, 1984

[8] Glowinski, R.; Lions, J.-L.; Tremolières, R. Numerical Analysis of Variational Inequalities, North-Holland, Amsterdam, 1981