Rgularisation dun problme dobstacle bilatral

B. Achachab(*) A. Addou(**), J. Zahi(*)

(*)Laboratoire de Modlisation et calcul conomique, Universit Hassan I, Settat, Maroc. (**)Dpartement de Mathmatiques et Informatique, Universit Mohammed I, Oujda, Maroc. [email protected], [email protected], [email protected]
RSUM. Le but essentiel de ce travail est la rsolution dun problme dobstacle bilatral laide dune mthode de rgularisation . ABSTRACT. The main objective of this work is the resolution of an bilateral obstacle problem by means of the regularization method. MOTS-CLS : rgularisation, inquation variationnelle, obstacle, problme, bilatral. KEYWORDS : regularization, variational inequality, obstacle, problem, bilateral.

Volume 5 - 2006, pages 350 360 - ARIMA

351 - ARIMA - Volume 5 - 2006

1. Introduction

1 /2

Soient un domaine born de R I n de frontire = assez rgulire, g 1 () et 1 ,2 H (), tels que T r(1 ) g T r(2 ) sur .

On considre le problme dinquation variationnelle - dit problme dobstacle bilatral - : Trouver u K = {v H 1 () : v = g sur , 1 v 2 p.p. dans } telle que u(v u) + f, v u 0 v K

[1]

[2]

o f L2 (). Il est bien connu que ce problme admet une solution unique (voir D. Kinderlehrer et G. Stampacchia). On suppose que 1 , 2 L2 (), avec ces hypothses, nous pouvons crire un autre problme dinquation variationnelle sans contraintes, quivalent (1) (2). Ce problme quivalent fait intervenir un terme non differentiable. Notre ide de rgularisation consiste remplacer ce terme par un autre qui est differentiable, on obtient ainsi une famille de problmes dquations aux drives partielles dont la suite des solutions converge vers la solution du problme initial. Par la mthode de dualit par les fonctions conjugues (voir I. Ekeland, R. Temam) on fournit une estimation a-posteriori de lerreur qui est desire pour limplementation pratique de la mthode de rgularisation.

ARIMA - numro spcial TAM TAM'05

Rgularisation d'un problme d'obstacle bilatral - 352

2. Reformulation et rgularisation
2.1. Reformulation de (1)-(2)
Soit un domaine born de IRn , de frontire = assez rgulire. Soient f L2 (), g H 1/2 (), 1 , 2 H 1 (), quon suppose vrier : 1 g 2 p.p. sur . 1 , 2 L2 () On notera : F = f 2 L2 (). On dnit : F = max{F, 0} , F et = 1 2 . Soit le convexe :
+

[3]

= max{0, F }

1 K = {v H g () : 1 v 2 p.p. sur }.

On considre le problme dinquation variationnelle -dit problme dobstacle bilatral- : (P1,2 ) o a(., .) est dni par a(u, v ) =

Trouver u K a(u, v u) + f, v u 0

pour tout v K

u.v,

u, v H 1 ().

Il est bien connu que le problme (P1,2 ) admet une solution unique (voir D. Kinderlehrer et G. Stampacchia). Dans toute la suite on utilisera la mme notation g pour designer un lment de H 1/2 () et un lment de H 1 () dont la trace sur est g , de mme par k = g 1 /, on dsignera aussi un lment k de H 1 () dont la trace sur est k . On va crire le problme dobstacle (P1,2 ) sous une nouvelle forme.

ARIMA - numro spcial TAM TAM'05

353 - ARIMA - Volume 5 - 2006

Thorme 1 : u est solution du problme (P1,2 ) si et seulement si w = u 1 est solution du problme dinquation variationnelle suivant : (P )
1 Trouver w Hk () 1 a(w, v w) + (v ) (w) + (H , v w) 0 pour tout v Hk ()

[4]

o est la fonctionnelle dnie par

1 (v ) = (F , (v + ) ) + (H + , v + ) v Hk (),

et H = F + Preuve :(voir A. Addou, E.B. Mermri)

2.2. Rgularisation de (P )
Lexpression de la fonctionnelle nous fait penser prendre sous la forme suivante : (v ) =

H + (v ) + F (v + ), ( > 0, destin tendre vers zero).

Or et sont deux fonctions dnies de R I vers R I , drivables et convexes vriant :

0

lim (t) = t+

et

lim (t) = t

Denition 1 : Le problme rgularis de (P ) est : (P )

1 w Hk () 1 a(w , v w ) + (v ) (w ) + H , v w 0 v Hk ()

[5]

le problme (P ) admet une solution unique (voir D. Kinderlehrer et G. Stampacchia ). A prsent on peut proposer les trois exemples suivants de et qui vrient les conditions ci-dessus : t 2 si t 2 t 1 si 0 t (t) = 2 0 si t 0. t
1 t2 2( 2

2 (t) =

si t + ) si 0 t si t 0.

ARIMA - numro spcial TAM TAM'05

Rgularisation d'un problme d'obstacle bilatral - 354

3 (t) = 0

t2 + 2

si t 0 si t 0.

1 (t) =

t2 2

si t 0. si t 0. si t .

2 (t) =

1 t2 2(

si t 0 + ) si t 0 si t . si t 0 si t 0.

3 (t) =

t2 + 2

Pour tablir la convergence de la suite w nous avons besoin nous avons besoin des rsultats suivants (voir R. Glowinski, J.L. Lions, R. Trmolires ) : Lemme 1 Soit V un espace de Hilbert, a : V V R I une forme bilinaire continue et V-elliptique, j : V R I une fonction faiblement continue convexe positive et f une forme linaire, continue dans V . On suppose que j : V R I , ( > 0), une famille de fonctions convexes semicontinues infrieurement (s.c.i.) et positives vriant : j (v ) j (v ) v V, u u faiblement dans V alors on a j (u) lim inf j (u ).
0

[6] [7]

Soient u et u les solutions respectivement des inquations variationnelles suivantes : a(u, v u) + j (v ) j (u) + f, v u 0, v V v V

a(u , v u ) + j (v ) j (u ) + f, v u 0, alors on a u u dans V quand 0. Lemme 2 Soient j (v ) =

(v )dx, j (v ) =

(v )dx

Si (t) (t) uniformment en t, quand 0, alors (6) et (7) sont vris. [8]

ARIMA - numro spcial TAM TAM'05

355 - ARIMA - Volume 5 - 2006

Lemme 3 : On a : | (t) (t)| c1 et | (t) (t)| c2 pour c1 , c2 deux constantes indpendantes de . Proposition 1 : Soit w (resp w) la solution de (P ) (resp (P )), alors on a : (w ) converge vers w dans H 1 (). Preuve : Puisque les , vrient lingalit (9), alors daprs les Lemmes 1 et 2, on a la convergence w w quand 0 dans H 1 (). Prenant v = w (resp v = w) dans linqualit du problme (P ) (resp (P )), on obtient : a(w w , w w ) (w ) (w ) + (w) (w). alors w w Donc w w
2 H 1 () 2 H 1 ()

t R I ,.

[9]

|(w ) (w ) + (w) (w)|. |(w ) (w )| + | (w) (w)|.

Daprs lingalit (6), on a : u u

2 H 1 ()

(F c1 + H + c2 )dx.

Par consquent, on obtient lstimation a priori suivante : u u

H 1 ()

2(F c1 + H + c2 )) 2

dx

ARIMA - numro spcial TAM TAM'05

Rgularisation d'un problme d'obstacle bilatral - 356

3. Estimation a-posteriori de lerreur

Dans cette section nous allons utiliser la dualit par les fonctions conjugues an de dterminer une estimation a-posteriori de lerreur de la solution du problme approch (P ). Nous avons besoin des rsultats prliminaires suivants (voir I. Ekeland, R. Temam ): Soient V et V (resp Y et Y ) deux espaces vectoriels topologiques mis en dualit par la forme bilinaire ., . V (resp ., . Y ). I = R I {, +}, sa fonction conjugue est Soit une fonction de V dans R dnie par (v ) = sup v , v V (v ), v V .
v V

On suppose quil existe un oprateur linaire L continu de V dans Y , L L(V, Y ), de transpos L L(Y , V ). Soit J une fonction de V Y dans R I , on considre le problme de minimisation suivant : u V, J (u, Lu) = inf J (v, Lv ) [10]
v V

La fonction conjugue de J est donne par J (y , v ) = sup { v , v

v V,y Y V

+ y , y

J (v, y )}

Thorme 2 : Soient V est un espace de Banach rexif et Y un espace vectoriel norm. Soit J : V Y R I une fonction strictement convexe s.c.i. propre et vriant : (i) u0 V , tel que J (u0 , Lu0 ) < et y J (u0 , y ) est continue en Lu0 . (ii) J (v, Lv ) +, quand v V +, v V . Alors le problme (10) admet une solution unique u, et J (u, Lu) = inf J (v, Lv ) = sup J (y , L y ).
v V y Y

Soit un sous ensemble ouvert de R I n , g : R InR I une fonction de Carathodory n c..d., s R I , x g (x, s) est une fonction mesurable et pour presque tout x , La fonction s g (x, s) est continue. Alors la fonction conjugue de la fonction G(v ) =

g (x, v (x))dx

(suppose dnie sur un espace fonctionnel V ) est G (v ) =

g (x, v (x))dx,

v V

ARIMA - numro spcial TAM TAM'05

357 - ARIMA - Volume 5 - 2006

o g (x, y ) = sup {ys g (x, s)}. sR IN Pour le problme (P ) on prend : V = H 1 (), Y = Y = (L2 ())n L2 () Lv = (v, v ) J (v, Lv ) = H (v ) + G(Lv ) 0 si v = k sur H (v ) = + ailleurs. 2 + + G(y ) = 1 2 |y1 | + H y2 + F (y2 + ) + H y2 dx o y = (y1 , y2 ) avec y1 (L2 ())n et y2 L2 (), de mme on utilise cette notation pour y Y . Daprs le Thorme 1, la fonctionnelle conjugue de J est donne par : J (y , L y ) = H (L y ) + G (y ) donc le problme (P ) peut scrire sous la forme (1.10), cest dire : u H 1 (), J (u, Lu) =
v H 1 ()

inf

J (v, Lv ).

Lemme 4 : Pour y = (y1 , y2 ) avec y1 (L2 ())n et y2 L2 (), on a : 1 2 |y | + ky1 + ky2 dx si divy1 + y2 = 0 dans 2 1 J (y , L y ) = et y2 F , y2 H + ailleurs. [11] Preuve : On a : H (L y ) = = =

sup { Lv, y H (v )}
v H 1 ()

sup
1 () v Hk

(vy1 + vy2 )dx + ky2 )dx +

(ky1

sup
1 () v H0

(vy1 + vy2 )dx

(ky1 + ky2 )dx

si divy1 + y2 = 0 dans

ailleurs

ARIMA - numro spcial TAM TAM'05

Rgularisation d'un problme d'obstacle bilatral - 358

et G (y ) = = sup { y , y G(y )}
y Y

1 + sup (y1 y 1 y2 y2 |y1 |2 H + y2 F (y2 + ) H y2 )dx 2 y Y 1 2 | )dx si y2 F et y2 H + ( |y1 = 2 ailleurs. Alors

J (y , L y ) =

1 2 |y | + ky1 + ky2 dx 2 1

si divy1 + y2 = 0 dans et y2 F , y 2 H + ailleurs.

Lemme 5 : Pour w et w solution de (P ) et (P ) respectivement, on a : 1 2 (w w)

2 L2 ()

1 1 2 + |w |2 +H + w +F (w + ) +H w + |y1 | +ky1 +ky2 dx. 2 2

y = (y1 , y2 ) Y , avec divy1 + y2 = 0 et y2 F , y2 H + dans .

Preuve : On a : J (w , Lw ) J (w, Lw) = 1 1 + |w |2 |w|2 + H + w H + w+ + F (w + ) 2 2 F (w + ) + H (w w)dx.

Utilisons linquation de (P) avec v = w , on obtient : J (w , Lw ) J (w, Lw) 1 2 (w w)

2 L2 ()

Appliquons le thorme 1 et utilisons le lemme 1, alors on a : J (w , Lw )J (w, Lw)

1 1 2 + |w |2 +H + w +F (w + ) +H w + |y1 | +ky1 +ky2 dx. 2 2

y =

( y1 , y2 )

Y , avec divy1 + y2 = 0 et y2 F , y2 H + dans .

ARIMA - numro spcial TAM TAM'05

359 - ARIMA - Volume 5 - 2006

Proposition 2 : Si w et w sont les solutions de (P ) et (P ) respectivement alors on a lestimation a-posteriori suivante : 1 2 (w w)

2 L2 () + H + w + F (w + ) w (H + (w )+ F (w + )).

[12]

Preuve : Puisque est diffrentiable linquation de (P ) est quivalente : 1 w Hk () : a(w , v ) +

(H + (w ) + F (w + ) + H )vdx = 0

[13]

do w vrie le problme de Dirichlet suivant : w + H + (w ) + F (w + ) + H = 0 dans w = k sur . Si on pose

y1 = w et y2 = (H + (w ) + F (w + ) + H )w .

alors on a :
divy1 + y2 = 0 et y2 F , y2 H +

Donc daprs le lemme 5 on a une estimation a-posteriori de lerreur :

1 2

(w w)

2 L2 ()

+ w (w k ) + H + w + F (w + ) + H w

+k (H + (w ) + F (w + ) + H ). [14] Si on choisit v = w H 1 () dans (13) on aura : a(w , w ) +

(H + (w ) + F (w + ) + H )w dx = 0

do w (w k ) +

(H + (w ) + F (w + ) + H )(w k ) = 0

Lestimation a-posteriori de lerreur (14), devient : 1 2 (w w)

2 L2 () + H + w + F (w + ) w (H + (w )+ F (w + )).

ARIMA - numro spcial TAM TAM'05

Rgularisation d'un problme d'obstacle bilatral - 360

Pour le choix (1) et (2) de et , lestimation a-posteriori (12) de lerreur est : 1 2 (u u)

2 L2 ()

H + w (1
[0w ]

w )+

F (w + )(1
[w + 0]

w + ).

Pour le choix (3) de et , lestimation a-posteriori (12) de lerreur est : 1 2 (u u)

2 L2 () [w0 ]

H + w (1

w )+ 2 + 2 w

F (w + )(1
[w + 0]

w + (w + )2 + 2

).

Il est facile de voir que le second membre de cette estimation dans les trois cas tend vers 0 lorsque 0. Acknowledgments The authors would like to thank the anonymous referees for interesting remarks. This work was supported in part by the Volkswagen Foundation Grant number I/79315.

4. Bibliographie
[1] A. A DDOU , E.B. M ERMRI, Sur une mthode de rsolution dun problme dobstacle. , Math-Recherche & Applications, N 2, (2000), pp. 59-69, [2] I. E KLAND , R. T EMAM, Analyse convexe et problmes variationnels. , Gauthier-Villars, (1974), [3] R. G LOWINSKI , J.L. L IONS , R. T RMOLIRES, , Numerical analysis of variatinal inequalities. , North-Holland, Amsterdam, (1981), [4] W. H AN , J. Z HOU, The regularization method for an obstacle problem. , Numer. Math. 69, p. 155-166, (1994). [5] D. K INDERLEHRER ET G. S TAMPACCHIA, An introduction to variational inequalities and their applications. , Academic Press, (1980).

ARIMA - numro spcial TAM TAM'05