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

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

doi:10.1016/S1631-073X(03)00025-6

Numerical Analysis/Calculus of Variations

A penalty/Newton/conjugate gradient method for the solution of obstacle problems
[Sur une méthode de pénalité/Newton et gradient conjugué pour la résolution de problèmes d'obstacles]

Présenté par : Philippe G. Ciarlet

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

¹ University of Houston, Department of Mathematics, Houston, TX 77204-3476, USA

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

Résumé
Abstract

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 $K = {v ∣ v \in H_{0}^{1} (Ω)$ , v⩾ψ p.p. sur $Ω}$ . 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 $K = {v ∣ v \in H_{0}^{1} (Ω)$ , v⩾ψ a.e. on $Ω} .$ 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.

Reçu le : 2002-11-18
Accepté le : 2002-11-19
Publié le : 2003-03-01

DOI : 10.1016/S1631-073X(03)00025-6

Affiliations des auteurs :

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

¹ University of Houston, Department of Mathematics, Houston, TX 77204-3476, USA

@article{CRMATH_2003__336_5_435_0,
     author = {Roland Glowinski and Yuri A. Kuznetsov and Tsorng-Whay Pan},
     title = {A {penalty/Newton/conjugate} gradient method for the solution of obstacle problems},
     journal = {Comptes Rendus. Math\'ematique},
     pages = {435--440},
     publisher = {Elsevier},
     volume = {336},
     number = {5},
     year = {2003},
     doi = {10.1016/S1631-073X(03)00025-6},
     language = {en},
}

TY  - JOUR
AU  - Roland Glowinski
AU  - Yuri A. Kuznetsov
AU  - Tsorng-Whay Pan
TI  - A penalty/Newton/conjugate gradient method for the solution of obstacle problems
JO  - Comptes Rendus. Mathématique
PY  - 2003
SP  - 435
EP  - 440
VL  - 336
IS  - 5
PB  - Elsevier
DO  - 10.1016/S1631-073X(03)00025-6
LA  - en
ID  - CRMATH_2003__336_5_435_0
ER  -

%0 Journal Article
%A Roland Glowinski
%A Yuri A. Kuznetsov
%A Tsorng-Whay Pan
%T A penalty/Newton/conjugate gradient method for the solution of obstacle problems
%J Comptes Rendus. Mathématique
%D 2003
%P 435-440
%V 336
%N 5
%I Elsevier
%R 10.1016/S1631-073X(03)00025-6
%G en
%F CRMATH_2003__336_5_435_0

Roland Glowinski; Yuri A. Kuznetsov; Tsorng-Whay Pan. A penalty/Newton/conjugate gradient method for the solution of obstacle problems. Comptes Rendus. Mathématique, Volume 336 (2003) no. 5, pp. 435-440. doi : 10.1016/S1631-073X(03)00025-6. https://comptes-rendus.academie-sciences.fr/mathematique/articles/10.1016/S1631-073X(03)00025-6/

Bibliographie
Cité par

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

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

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

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

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

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

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

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

Cité par Sources :

Commentaires - Politique