140 13 23MB
French Pages 529 Year 1984
V. SMIRNOV
COURS de ,. MATHEMATIQUES SUPÉRIEURES tome IV deuxième partie
~DITIONS MIR • MOSCOU
Traduit du rùsse par Djtlali Embarek
,
©
"
"
lfSAaTeXbCTBO «HaYKal) fxaBHaH peAaK~IUI «Ï>I13I1KO-MaTeMaTiPIeCKOH JI1ITepaTypLI, 1981, C I13MeHeHIIHMII
© Traduction française Editions Mir 1984
TABLE DES MATIERES
Préface CHAPITRE
... .. . . . . . . . . . ... .. ... . . . . . . PREMIER. TH~ORIE G~NÉRALE DES ÉQUATIONS D~RIVÉES PARTIELLES • . • . . . •
9
AUX
fO
1-1. Equations du premier ordre • • • • • • • • . • • • • •• 1-1-1. Equations linéaires à deux variables indépendantes (10). 1-1-2. Problème de Cauchy et caractéristiques (13). 1-1-3. Cas d'un nombre quelconque de variables (18). 1-1-4. Exemples (23). 1-1-5. Théorème auxiliaire (24). 1-1-6. Equations non linéaires du premier ordre (28). 1-1-7. Variétés caractéristiques (32). 1-1-8. Méthode de Cauchy (32). 1-1-9. Problème de Cauchy (35). 1-1-10. Unicité de la solution (37). 1-1-11. Cas singulier (40). 1-1-12. Cas d'un nombre quelconque de variables indépendantes (42). 1-1-13. Intégrales complète, générale et singulière (45). 1-1-14. Intégrale complète et problème de Cauchy (47). 1-1-15. Exemples (49). 1-1-16. Cas d'un nombre quelconque de variables (53). 1-1-17. Théorème de J acobi (55). 1-1-18. Systèmes de deux équations du premier ordre (57). 1-1-19. Méthode de Lagrange-Charpie (59). 1-1-20. Systèmes d'équations linéaires (61). 1-1-21. Systèmes complets et jacobiens (63). 1-1-22. Intégration des systèmes complets (65). 1-1-23. Crochet de Poisson (66). 1-1-24. Méthode de Jacobi (69). 1-1-25. Systèmes canoniques (71). 1-1-26. Exemples (72). 1-1-27. Méthode des séries majorantes (73). 1-1-28. Théorème de Kovalevskaïa (76). 1-1-29. Equations d'ordre supérieur (82).
fO
1-2. Equations d'ordre supérieur . . . . . . . . . . . .. 1-2-1. Types d'équations du second ordre (84). 1-2-2. Equations à coefficients constants (86). 1-2-3. Formes normales dans le cas de deux variables indépt'ndantes (88). 1-2-4. Problème de Cauchy (91). 1-2-5. Bandes caractéristiques (94). 1-2-6. Dérivées d'ordre supérieur (96). 1-2-7. Caractéristiques réelles et imaginaires (99). 1-2-8. Théorèmes fondamentaux (100). 1-2-9. Intégrales intermédiaires (102). 1-2-10. Equations de Monge-Ampère (103). 1-2-11. Caractéristiques dans le cas de n variables (104). 1-2-12. Bicaractéristiques (107). 1-2-13. Lien avec un problème variationnel (111). 1-2-14. Propagation de la surface de discontinuité (113). 1-2-15. Discontinuités fortes (116). 1-2-16. Méthode de Riemann (119). 1-2-17. Conditions initiales
84
6
T ABLE DES MATŒRES
caractéristiques (124). 1-2-18. Théorèmes d'existence (125). 1-2-19. Formule d'intégration par parties et formule de Green (130). 1-2-20. Méthode de Volterra (132). 1-2-21. Formule de Sobolev (136). 1-2-22. Formule de Sobolev (suite) (140). 1-2-23. Construction de la fonction cr (142). 1-2-24. Conditions initiales générales (14'i). 1-2-25. Equation des ondes généralisée (149). 1-2-26. Cas d'un nombre quelconque de variables indépendantes (151). 1-2-27. Inégalité énergétique (154). 1-2-28. Théorèmes d'unicité et de dépendance continue des solutions (160). 1-2-29. Cas de l'équation des ondes (163). 1-2-30. Théorème d'immersion dans l'espace des fonctions continues et certaines de ses conséquences (166). 1-2-31. Solutions distributionnelles des équations du second ordre (171). 1-2-32. Sur l'existence et l'unicité des solutions distributionnelles du problème de Cauchy pour l'équation des ondes (177). 1-2-33. Equations de type elliptique (178). 1-3. Systèmes d'équations 1-3-1. Caractéristiques de systèmes d'é~uations (184). 1-3-2. Conditions cinématiques de compatibilite (188). 1-3-3. Condi. tions dyn~m,iq:ue~ O. Par analogie avec l'espace à trois dimensions on se servira de termes géométriques pour étudier l'équation (16). Les énoncés et les démonstrations ne seront pas donnés dans le détail. Nous avons affaire à un espace à (n 1) dimensions de point générique (Xl' . . . , X n , u). Appelons variété à m dimensions de cet espace, l'ensemble des points dont les coordonnées s'expriment en fonction de m paramètres arbitraires:1
+
+ + ... +
+
Xli
= Xli (tH ••• ,
t m );
u = u (tu ••• , t m )
(k =
1, 2, •.. , n).
On admet que m de ces équations sont résolubles par rapport à t l , • . • , t m . Pour m = n on a une variété à n dimensions que l'on appellera surface. Si l'on prend Xl' . . ., X n pour paramètres, on obtiendra l'équation explicite de la surface: u = u (Xl' .•. , x n ). C'est cette forme que doit avoir l'équation de la surface intégrale de l'équation (16). Pour m = 1, la variété à une dimension correspondante s'appelle ligne de l'espace à (n 1) dimensions. Définissons les caractéristiques de l'équation (16) à l'aide du système suivant:
+
dXk
~=ak
) (Xlt ... , Xn , u;
du ----a;s=C (Xl'
... , Xn , U),
(17)
où s est un paramètre auxiliaire. Toute solution de ce système définit une ligne de l'espace à (n 1) dimensions, car il n'existe pas de solution dont tous les Xk et u soient constantes, puisque par hypothèse a~ a; a; > O. Les coordonnées des points de cette ligne s'exprimeront en fonction du paramètre s. Pour construire une surface avec ces lignes, il nous faut prendre une famille de telles lignes dépendant de (n - 1) paramètres arbitraires. On obtient un ensemble de points dépendant de n paramètres. Si une surface régulière u = u (Xl' . . . , X n) est formée par une famille de caractéristiques dépendant de (n - 1) paramètres, alors c'est une surface intégrale de l'équation (16). En effet, en dérivant u {Xl' ••.
+
+ + ... +
1-1-3. CAS D'UN NOMBRE QUELCONQUE DE VARIABLES
19
..., xn) par rapport à s et en se servant des équations (17), on obtient ft
~
du
dS = LJ ux,.,a,.,. k=l
Comme du = c en vertu de la dernière équation (17), il vient l'équada tion (16). Réciproquement, toute surface intégrale régulière u = = U (Xl' .•. , xn) peut être engendrée par une famille de caractéristiques, qui dépend de (n - 1) paramètres. En effet, si l'on connaît la surface intégrale u = u (XH' •• , xn), on peut déterminer les X,.. . à partir du système d'équations
d~1&==a",[xlt •.• ,Xn ,
u(xl
, •••
,xn )]
(k=1, .o.,n),
(18}
ce qui nous donne (n - 1) constantes arbitraires. La constante qui figure par sa somme avec s n'est pas essentielle. En portant la solution du système (18) dans le second membre de u = u (Xl' .•. , xn),. en dérivant par rapport à s et en se servant des équations (16) et (18), on s'assure que u vérifie la dernière des équations (17). On admet comme dans [1-1-2] que u (Xi' ... , xn) et les seconds membres des équations (17) possèdent des dérivées premières continues. Le problème de Cauchy pour l'équation (16) revient à déterminer une surface intégrale contenant une variété à (n - 1) dimensions donnée: Xk = Xk (t l , • • ., t n -1) ; U = U (t l , • t n -1) (k = 1, ..., n), (19) les seconds membres de ces équations étant continus et possédant des dérivées partielles premières continues à l'intérieur d'un do,.. maine D de l'espace à (n - 1) dimensions (t l , . . . , t n - l ). On admet que le rang de la matrice formée avec les dérivées ôô:X k est égal à (n - 1) et qu'à des systèmes différents de valeurs tl (t l , ••• , t n - l ) correspondent des points (Xl' • . • , xn) différents.. Par ailleurs, comme indiqué plus haut, on admet que les coefficients ak (Xl' •.• , Xn' u) et c (Xl' •• 0' Xn, u) possèdent des dérivées premières continues dans un domaine de l'espace, contenant. la variété (19) en son intérieur. En particulier, cette condition du problème de Cauchy peut setraduire par la donnée de la fonction cherchée u comme fonction h à Xl fixe, des autres variables: u IX1=xiO) = cp (x 2 , ••• , Xn)' (20) Ce problème se résout comme dans le cas de deux variables indépendantes. On prend les expressions (19) pour conditions initiale~ 0
2*
.,
20
CH. 1. TH:eORIE DES :eQUATIONS AUX D:eRIV:eES PARTIELLES
dans l'intégration du système (17). On obtient ainsi la solution sous la forme Xl = Xlt (s, Le déterminant ÔXn , · .. , as ÔS ôX2 ÔXn ·.. , at , atl ' . . . . . . • . . . . . . l. . 8x'j lJz n ÔXI , , · .. , at n - l ôt n ôt n- l
7iï , ÔXI , ôt l ÔZI
â=
Dx'j
(22)
l
jouera un rôle essentiel dans la suite. Compte tenu de l'équation (17), on peut le mettre sous la forme al'
â=
a 2,
ÔZI
ÔZ"
at l '
ôtl '
• • • 8Zi , 8t n- 1
.• • • 8z" atn-l'
... , an ÔZn ... , ~ • • • . • .
... ,
(23)
•
ÔXn ~
n-l
Si ce déterminant est non nul sur la variété (19), c'est-à-dire pour $ = 0, alors les n premières équations (21) sont soluble~' par rapport à s, tu ... , t n -1. En portant cette solution dans la dernière des équations (21), on obtient une surface intégrale de l'équation (16). Le problème de Cauchy n'admet pas d'autres solutions dans ce cas. Ceci se démontre exactement de la même façon que pour·deux variables indépendantes. Considérons le cas où la condition initiale est de la forIlle (20); le rôle des paramètres t l , . • • , t n -1 est assumé par X2' .' •• , Xn. Considérons l'équation linéaire et. supposons que le déterIllinant (23) est différent de zéro sur la variété (19). Comme ~=p= 0 pour P =1= q et = 1, on obtient d , al '=1:- o. Un~ diiision par le coefficient al ~ous conduit à l'équation
:;p
Pl
+ a2 (Xl'
••• , Xn )
P2
+ · · · + an (Xl'
• •• , X n) ,
Pn ..:-: ,
= b (Xit ••• , Xn)U+ C (XI~' '. :., Xn). (24)
Supposons que alt, b et c sont continus et ,possècie~t. des ,dérivées partielles premières par rapport à Xi' •• • '~n continùes poùr a~ ~. Xl ~ Ji ~t X 2 , • • • , X n quelconques. Suppbsons en outre que sous 'ces ..conditions les fonctions 'indiquées sont bornées:' 1ali··1 ~ M; ·lbl~'M;''I'c'I~M. - "'.' l, h ' . : ,
1-1-3. CAS D'UN NOMBRE QUELCONQUE DE VARIABLES
2t
Prenons xi pour variable indépendante et mettons le système (17) sous la forme
= ait (xu
dxll d.%1
::i
(k = 2, ... , n),
... ,xn )
= b (Xl' ••• , X n ) U
+ C (Xi!
x n )·
••• ,
(25) (26),
Soit x\O) E[a,~] la valeur initiale de Xl. Intégrons le système (25) avec des conditions initiales quelconques: x" 1X1=Cl(1) =x~O)
De
l'ak
(k=2, ... , n).
1~M il s'ensuit que les solutions
Xk
possèdent des dérivées bornées 1 ::: 1~M, donc
du système (25) Xk
restent bornées
en valeur absolue: 'Xk-X~O) 1~M (~-a). En appliquant la méthode des approximations successives (tome II, [11-3-2]) on s'assure immédiatement que les solutions indiquées (k = 2, ... , n) existent sur l'intervalle a ~ Xl ~ ~ tout entier pour des condition~ initiales arbitraires x1&O) (k = 2, ... , n). On peut dire que la courbe intégrale qui passe par le point A o (x~O), X~O), ••• t x~» arrive au point A (Xl' X 2 ' • • • , x n ) dont les coordonnées sont définies par les formules (27). Le théorème d'unicité nous permet d'affirmer que; si l'on prend le point A pour point initial, alors la courbe intégrale correspondante passera par le point A o• D'où il s'ensuit que les équations (27) sont solubles pour Xk quelconque par rapport à X~O), ••• • . •, x~) et la solution est de la forme I l ' Xl' X , • • • , X ) (k =2, ... , n ) . Xk(0) -- -définissent, quelles que soient a fixe et b = w (a), la solution du .système (107) qui vérifie la condition (106). On peut admettre que les formules (121) définissent quatre des variables x, y, U, p, q en fonction de la cinquième et des trois constantes arbitraires a, b et c = w' (a). L'intégrale générale du système >(107) contient quatre constantes arbitraires. Mais la relation (106) nous dit que la famille .des bandes caractéristiques ne doit dépendre que de trois constantes arbitraires. C'est ce qu'expriment les formules (121). Dans un prochain paragraphe consacré au cas d'un nombre quelconque de variables indépendantes, on s'assurera par des calculs ,directs que les équations (121) donnent bien la solution du système {107). Voyons maintenant s'il est possible de définir une intégrale singulière directement à partir de l'équation différentielle sans le se.cours d'une intégrale complète. Une dérivation de (110) par rap-
l-l-t5. EXEMPLES .
49
port à a et à b nous donne
0; F uCPb + FpCPXb + F qCPub = O. La définition de l'intégrale singulière (113) nous permet d'affirmer que sur la surface intégrale singulière on a : FuCPa
+ Fpxa + F qCPu a = FpCPxa
+ F qCPu a = 0;
FpfPxb
+ F qub = O.
On admettra que le déterminant de ce système homogène par rapport à F p et F q ne s'annule pas sur la surface intégrale singulière, ce qui revient en fait à admettre que les équations (109) sont solubles par rapport à a et b. Ce système nous donne F p = F q = O.
(122)
Donc, l'intégrale singulière peut être obtenue par élimination de p et q entre les trois équations suivantes:
P, q) = 0; F q (x, y. u. Pt q) = O. (123) Les équations (122) indiquent qu'lI est impossible d'appliquer le théorème des fonctions implicites à l'équation (106) relativement à la variable P ou q. Ceci exprime qu'il est impossible d'obtenir l'intégrale singulière en résolvant le problème de Cauchy ainsi que nous l'avons fait dans [1-1-9], en admettant que l'équation était résolue en p (ou q). On aurait pu arriver à ce résultat par une autre voie. Quelle que soit la courbe considérée sur la surface intégrale singulière, le déterminant (96) sera nul le long de cette courbe en vertu des conditions (122). Ceci montre que le problème de Cauchy n'admet pas de solution pour toute courbe de la surface intégrale singulière. F (x, y,
U,
P, q)
=
0; F p (x, Y.
Ut
1·1·15. Exemples. 1. L'équation U
+
= xp
yq
+ f (P,
q)
(124)
est analogue à l'équation de Clairault étudiée au (tome II, [l-t-tin. En remplaçant p et q par a et b on obtient immédiatement l'intégrale complète: u = ax
+ by + f (a,
L'équation u = xp
b).
+ yq+ pq
admet pour intégrale complète:
u=
ax
+ by + ab
En appliquant la méthode indiquée ci-dessus, on obtient l'intégrale singulière: 'U
== -xY.
Si l'on considère une courbe quelconque XD
'-01017
=
cp (t) ; Ye
= 1{> (t) ;
Uo
=
-cp (t) 1{> (t)
(125)
50
CH. 1. TH1!:ORIE DES 1!:QUATIONS AUX D1!:RIV1!:ES PARTIELLES
de cette surface, alors les équations (86)
'l' (t) qo + Poqo + cp (t).'1' (t) = 0, cp' (t) '1' (t) + 'l" (t) cp (t) + cp' (t) Po + '1" (t) qo = admettent les solutions Po = -'l' (t), qo = -cp (t) et -le long de la courbe (125) cp (t) Po +
°
on aura
Fpo=q+cp(t) =:0;
Fq.=p+'P(t)=:O.
L'équation
1
u=xp+yq-T (p2+ q2)
admet pour intégrale singulière 1 u= l (x, a) et q = q>2 (y, a) et l'on obtient l'intégrale complète u=
J
q>d x , a) dx+
J
(j>2
(y, a) dy+b,
.où b est la deuxième constante arbitraire. Appliquée à l'équation
pq-xy=O ou
.E...=.J!.. x q'
(129)
cette méthode nous donne 1 1 u="2 ax2 + 2a y2+b.
Soit à déte rminer la surface intégrale quîpasse par la courbe 1 x
= t;
Y= - ;
t
(130) ~
u = 1.
En portant ces expressions dans (130) et en dérivant par rapport à t, on trouve 1 1 t 1=2' at2 2atl +b; at- at 3 =0.
+
L'élimination de t entre ces équations nous donne b = 0 et nous obtenons une famille de surfaces intégrales à un paramètre. 1 1 u=2" ax 2 + 2a yI. L'env.eloppe de cette famille est la surface intégrale cherchée, soit:. u = xy. Si pour condition initiale on avait pris la courbe x = t; Y = t; u = t~,
. (t3t)
S2
CH. J. THSORJE DES :SQUATIONS AUX :DSRIVSES PARTIEI.J·ES
la méthod~ ;précédente nous aurait conduits aux équations ( ~ a+ ~ -1) ,1I+b=O; 2( ~, a+ ~ -1) , ' 0,
'"
d'où il s'ensuit que a =1,b = 0 et nous n'aurions pas pu trouver une surface intégrale passant par la courbe (131). Il est immédiat de voir qu'on peut prolonger la courbe (131) en une bande caractérist ique en posant p = , etq -:.' t. En effet, les fonctions, x = t; Y = t; u = t~; P = t; q = 1 vérifient l'équation (129) et le système (1.07).
4. Si l'équation ne contient pas les variah' et de la forme F (u, P, q) = 0
indépendantes, t'est-à-di11!
on peut déterminer une intégrale complète en cherchant la solution d, l'équa'tionsous la forme ' u = q> (x + ay), (t32) où
est une constante arbitraire. Voyons à titre d'exemple l'équation
CI
pq -
u
=
O.
(t33)
En effectuant la substitution (132) et en posant CI
[q>' (,)]' -
~
=
x
on trouve
dOLDe
l'intégrale
q> (,) = O.
L'intévation de cette équation différentielle ordinaire nous œmplète de l'équation (t33): (x+ay+b)1 u= 4a ,~
+ G1/,
système (68) s'écrit dans le cas de l'équation (133): ~
IIi"' ~ q;
~
~
et son intégration nous donne x = qoe'' (xo - qo) ; y = poe' ... ..'('
'.
~
d,::;= P; fiS = 2pq;
+ .
i
P
=
+ (Yo -
poe';
q
d,
Po) ;
=
= P; u
=
~
fiï=- q, Poqo~'
+ (Ue- PolJe) ;
qoe'.
Supposons qu'on cherche la, sUrface intégr~le qui p:asse par la droite ft;;; ,
,'Xo
=
t;
1I~ ,
t;
Uo
= t•
.on détermine Po et qo à partir des équations
=
t; Po = 1, a"où Po = J et qo=t., En, portant Po et qo, dans les trois premièreséquati()DS , Prio
(f3St) et en posant e'':'''- ;v~"on obtient les équations paramétriqueS deJ8'SUiface d1erchée en fonction des paramètres ,v ,et t: , , x 7""", Iv; 11 = v; u 7"""tv~,
:w,,' sous
la forme explicite,
.U, ~
xy. .'
I.;I-te~
CAS D'UN NOMBRE QUELCONQUE DE
:vARIAB:J,.ES
53:
'1..1·16. CaS d'un nombre quelconque de variables. On 'appellfJ, de l'équation F (Xit •• 0' Xn, U, Pl' •• 0' Pn) = 0 (134) une solution (135) intégrale complète
de cette équation contenant n constantes arbitr aires as et telle que l'élimination des as entre les équations P. = q>%J& (xi t
Xn , al'
0.0'
000,
an), (k = 1, 2,
00.'
n)
(135t )
et l'équation (135) donne l'équation (134).On admettra que ait sont des fonctions de n - 1 paramètres: ait
=
ait (tt,
0
•
0'
t n -l )
= 1, 2, ... ,
(k
n).
(136)
En portant ces expressions dans (135) et en éliminant les n - 1 paramètres entre les n équations u = Cf> (Xl' Xn , al' ••• , an), , Xn , al' ••• , an)=O (j=1, ... , n-1), 0
Cf>Cj(xh
•• ,
(137)
on obtient l'intégrale générale de l'équation (134). Cette intégrale dépend du choix des n fonctions (136). Passons maintenant au problème de Cauchy. On demande de trouver une surface intégrale de l'équation (134) contenant une variété de n - 1 dimensions donnée: u
=
u (tt, • -
0'
t n -l );
Xit
=
Xit (t lt
0
0
tn-t) (k = 1, 2, .. 0' n).
.,
(138)
Ce problème se résout exactement comme pour le cas de deux variables indépendantes. En portant les expressions (138) dans (135) OD obtient une égalité de la forme
t n -It al' •
an) = O.
(139) En adjoignant à cette égalité les n - 1 relations obtenues par dérivation de (139) par rapport à tlt t n -1' soit '" (t it
••
0'
0
"'Cl = 0 ; """2 = 0;
0
0
.,
.,
"""n-l = 0,
(140) on obtiendra n équations à partir desquelles on peut déterminer ait (k = 1, 2, ... t n) en fonction des paramètres tt, ., 0' t n -ll c'est-à-dire les fonctions (136). En portant ces fonctions dans (137) et en éliminant t lt t n -1 entre les n équations (137), on obtiendra la surface intégrale qui contient la variété (t38). A noter que dans la formule (139) le nombre de paramètres peut être inférieur 0
•• ,
o ••
;
S4
CH. I. THeORIE DES ::eQUATIONS AUX DeRIv:eES PARTIELLES
à n _. 1. Dans ce cas il faut dériver (139) par rapport à ces paramètres. Si l'on fixe les paramè~res tjt les n équations (137) à n 1 variables (u, Xl' ••• , x n ) définissent une courbe dans l'espace à n 1 dimensions de point générique (U, Xl' . . • , x n ). En adjoignant l'équation (135 1) à ces équations, on prolonge cette courbe en une .bande du premier ordre. Cette bande appartient à deux surfaces intégrales: à l'enveloppe qui s'obtient par élimination des paramètres i j entre les équations (137), et à une surface enveloppée. Donc, cette bande est caractéristique, c'est-à-dire vérifie le système de Cauchy (98). Ceci nous permet, si l'on connaît une intégrale complète (135), de trouver la solution du système (98) en fonction de 2n -1 constantes arbitraires. On admettra pour simplifier que an est une fonction de (al' .. ~, an -1)' ces derniers jouant le rôle des paramètres t l , • • . , t n -1' Les formules (137) et (135 1) deviennent:
+
+
ail "', an), (J)a j + cp an bi -:- 0 (j = 1, 2, ..., n - 1), U=(J)(Xl , " ' , X n ,
Ph=(J)Xh(X1,
... , X n , al' ... , an)
(k=1,
(141)
2, ... , n),
où par bi on a désigné la dérivée de an par rapport à ai' Les formules (141) définissent la bande du premier ordre mentionnée. Dans ces formules al' ... , an et bl , • • • , bn -1 sont arbitraires, puisque le choix de la fonction an (al' ... , an -1) l'est. Prouvons formellement que la bande définie par les formules (141) vérifie le système (98). En portant les expressions (141) de U et Ph dans l'équation (134), on obtient une identité en Xh et ah (k = 1, ... , n) dont la dérivation par rapport à as nous donne n
U(J)aj+'
LJ
Ph(J)Xkaj-O
k=1
(j
=
1, 2, ... , n!~l),
n
U(J)an
+Li
k=1
PkCf'xkan
=0.
En multipliant la dernière égalité par bit en l"aJoutant à laprécédente et en se servant de (141), on obtient les n ....... 1 égalités suivantes: l
n ,
.LJ
k= 1 :
•
Ph (1I'tp~ ne s'annule pas dans la partie considérée de D. On a cP~ = Ilq>u;
\j)~ =
121Py,
d'où q>x,!,,I -
lI'tpX!= (fI -
12)
fPy'tpye
(17)
De ces formules il s'ensuit que si le déterminant s'annule en un point, il en sera de même des deux dérivées partielles premières de fP ou 'l'. Il faut donc prendre des solutions de (16 1 ) et (16 2 ) dont les deux dérivées partielles premières ne sont pas simultanément nulles. Les fonctions cp et 'tp vérifient l'équation (15) et en vertu de (9) on a a' = c' = O. De la formule (10) il s'ensuit que b' =1== 0, si bien que l'équation (7) se réduit à la forme (13). Nous avons vu dans [I-1-2l que les solutions des équations (16 1) et (16 2 ) étaient locales, c'est-à-dire qu'elles ne sont distinctes de constantes que dans un domaine qui, en général, est une partie du domaine dans lequel Ih (x, y) sont continûment dérivables, et l'équation (7) ne peut être réduite à une forme normale que dans le domaine indiqué. Cette remarque sur le caractère local de la réduction de l'équation (7) à une forme normale est également valable pour la suite de l'exposé. Passons maintenant à une équation elliptique. Dans ce cas ac - b2 > et les racines de l'équation (14) sont conjuguées complexes. Considérons l 'u~e des équations (16):
°
Ux
=
-b+i vac-b2 a
u y•
Si l'on admet que les coefficients a, b et c sont des fonttions analytiques de x et y et que a =1= 0, on peut trouver une solution de cette équation sous forme d'une fonction analytique [1-1-28]: U = = cp (x, y) i'tp (x, y). De plus
+
ab CPy - .yac=a b b vac-b2 li 'tpy + a 2
CPx = 'l'x = -
'tpy, q>y.
Faisons la substitution (8). En se servant du système en fP et '1' et des formules (9), on trouve b'
= 0;
a'
=
c'
=
(ac - b2)
(cp~
+ 'tp~) ; .
1-2-4. PROBLt1ME DE CAUCHY
91
en divisant par a' on met l'équation sous la forme f)2 u
f)~2
f)2
+ f)11u + ... = o.
(18)
2
La formule (17) devient CPx'1'y-CPy'1'x= -
2 -.1-ac-b2
a
CPy'1'y.
Le problème est donc résolu pour le cas elliptique aussi. La résolution de ce problème dans le domaine tout entier, sous certaines conditions portant sur les coefficients a, b et c, qui ne sont pas considérés comme des fonctions analytiques, est accessible dans le travail de 1. V e k ua, Problème de réduction de formes différentielles elliptiques à la forme canonique et système généralisé de Cauchy-Riemann. DAN SSSR, 1955, 100, nO 2 (en russe); voir également l'ouvrage du même auteur: Fonctions analytiques généralisées, M., Fizmatguiz, 1959 (en russe). Il reste à étudier l'équation parabolique. Dans ce cas l'équation (14) possède une racine double et l'équation (15) nous donne des équations (16 1 ) et (16 2) identiques. Pour cp (x, y), prenons une solution de (16 1) et pour '1' (x, y), une fonction quelconque telle que le jacobien de cp et '1' soit non nul. Le coefficient a' sera nul dans l 'équation transformée en vertu du choix de cp (x, y). De plus, l'équation étant parabolique, on a ac - b2 = 0 et la formule (10) nous dit que b' = O. Donc a' = b' = O. La fonction c' ne peut être identiquement nulle, car le cas échéant on aurait obtenu une équation du premier ordre et la transformation inverse qui fait passer de (~, 11) à (x, y) ne nous aurait pas donné une équation (7) du second ordre. Donc, dans le cas parabolique, on obtient la forme canonique suivante: U1I11
+ ... =
0,
(19)
où les termes non écrits ne contiennent pas de dérivées secondes mais contiennent nécessairement la dérivée première par rapport à ~. 1-2-4. Problème de Cauchy. On a vu au [I-1-2Jl que pour conditions initiales de l'équation du second ordre F (x, y,
U,
p, q, r,
St
t)
=
0
(20)
on peut en particulier prendre la valeur de la fonction u et de sa dérivée U x = P au point x = x o :
u Ix=xo = cp (y);
P Ix=xo = '1' (y).
(21)
Ces conditions seront appelées conditions initiales spéciales. Ces conditions initiales reviennent à se donner la valeur de la fonction inconnue u et de sa dérivée partielle p le long d'une courbe x = x o du plan (x, y). A noter que les valeurs de l'autre dérivée partielle du
92
CH. I. TImORIE DES :eQUATIONS AUX D:eRIV:eES PARTIELLES
premier ordre, q 1%=%0 = cp' (y), s'obtiennent immédiatement à partir de la première condition (21). Donc, en vertu des conditions initiales, on connaîtra la fonction u et ses dérivées partielles premières le long de la courbe x = X o' Les conditions initiales se généralisent sans peine. Soient données dans le plan (x, y) une courbe  sans point double et les valeurs de la fonction inconnue u le long de cette courbe. On connaît ainsi les valeurs prises le long de 'A. par la dérivée de u suivant une direction tangente à Â. Pour déterminer la dérivée première suivant une direction quelconque, on doit disposer encore d'uIle donnée le long de Â, plus exactement, on doit connaître les valeurs prises le long de  par la dérivée de u suivant une direction quelconque distincte d'une direction tangente à Â. La connaissance des dérivées suivant deux directions du plan (x, y) le long de 'A. nous permet de définir la dérivée de u suivant une direction quelconque du plan (x, y) le long de 'A.. Donc, dans le cas considéré, on doit connaître les valeurs prises sur  par la fonction U et par sa dérivée suivant une direction quelconque non tangente à Â. La donnée des valeurs de u le long d'une courbe 'A. du plan (x, y) se traduit par la donnée d'une courbe l dans l'espace (x, y, u). On connaît de plus les valeurs des dérivées partielles p et q le long de Â. En définitive donc, les conditions initiales consistent à se donner une courbe l de l'espace (x, y, u) et la position du plan tangent le long de 1. Paramétriquement ces conditions se traduisent par la donnée de cinq fonctions (22) x (t), Y (t), u (t), P (t), q (t), liées par la relation du
=
p dx
+ q dy.
(23)
La relation (23) revient à exiger que la donnée des deux dérivées partielles p et q le long de 'A. ne contredise pas la donnée de la fonction u le long de 'A., c'est-à-dire que la dérivée de u suivant une ;direction tangente à Â calculée d'après p et q prenne les mêmes valeurs que eelles obtenues à partir de la donnée de u le long de 'A.. Les cinq fonctions (22) liées par la relation (23) définissent une bande dans l'espace (x, y, u) et le problème de Cauchy revient à chercher une surface intégrale de l'équation (20) contenant cette bande. Le problème de Cauchy se pose dans les mêmes termes dans le cas plus général où la fonction inconnue dépend d'un nombre quelconque de variables. Considérons par exemple l'équation différentielle du second ordre dans le cas de trois variables indépendantes: F(xl , x 2 , x a, u, U%l' U%I' u x3 , u x1xi ' U X1X2 ' " .)=0. (24) Les conditions initiales consistent ici à se donner la fonction U et ses dérivées partielles du premier ordre sur une surface S de l'espace (Xl' x 2 , xa). La fonction u étant donnée sur la surface S, pour défi-
1-2-~. PROBL~ME
DE CAUCHY
93
nir toutes seS dérivées partielles du premier ordre le long de S, il suffit de se donner le long de S une dérivée suivant une direction quelconque nOil contenue dans un plan tangent à S. Si la surface S support des conditions initiales est le plan Xl = xl°), on obtient les conditions initiales spéciales: u Ixl=xio> =
3 elles se définissent de façon analogue; en particulier, si S est définie par l'équation explicite où X' alors
x m = q> (Xl' •.. , Xm -1), = (Xl' . . . , x m - l ) parcourt un domaine borné D m - 1 de Rm-l., dS
V/1 +
= -.
mij1 q>i. i=1
(x') dx' ,
1
dx' =dx1
•••
dx m _1 ,
au point x=(x', q> (x'» ES et
1 J/ fi +): ",~, f dS =
(x')
(x') dx',
où j (x') est la valeur de j au point X = (x', q> (x'» E S. Si S est la frontière d'un domaine borné D de l'espace Rm et si S est une surface différentiable, on peut la subdiviser en un nombre fini de parties S ktk = 1, ... , N, définies chacune par une équation explicite exprimant l'une des coordonnées Xl1{. en fonction des autres, et définir l'intégrale) j dS sur S comme la somme des intégrales j dS sur Sk.
J
~
s
La formule d'intégration par parties s'écrit:
) ::i D
v dx =
- ) D
u
::i
dx
+ ) uv cos (n. Xi) dS, S
i
= 1, ... , m.
(107)
1-2-19. FORMULE D'1NTfiGRATION PAR PARTIES ET FORMULE DE GREEN
131
Cette formule est de toute évidence vraie dans le cas oùD est un domaine borné d'un espace euclidien Rm, sa frontière S est unehypersurface régulière et les fonctions u et v de la classe Cl (D) (c'est-à-dire la classe des fonctions continues et continûment dérivables dans D = = DUS). Dans cette formule, cos (n, Xi) représente le cosinus de l'angle de l'axe des Xi et de la normale extérieure n à S. La formule (107) résulte de la formule
~ ;~
dx=
D
~
w cos (n, xi) dx,
D
qui est v~lable quel que soit i = 1, ... , m pour S et w douées des'. propriétés de dérivabilité indiquées ci-de~sus. En effet, si w= uv, on obtient (107). Utilisons la formule (107) pour établir la formule de Green pour un opérateur différentiel linéaire du deuxième ordre. On supposera, que toutes les dérivées qui interviendront dans la suite sont continues dans lJ et que S est régulière. . Soit m
L (u) =
m
. Li ai1~uXiXk + k=1 2J bkuXk + cu,
(108)
t, k=1
et m
L*(v)= ~
(109)
i. k=l
où aik, bk et c sont des fonctions données de x. Considérons l'intégrale ) vL (u) dx,
dx = dX1
•••
dx m ,
D
et transformons-la à l'aide de la formule (107). Ceci nous conduit à la formule de Green
1
[vL (u) - uL* (v)) dx =
D
1
[vP (u) -- uP (v)
+ uvQJ as,
(110}
s
où
(111~
et n désigne la normale extérieure unitaire à S. 9*
132
CH. I. THeORIE DES ÉQUATIONS AUX DÉRIVÉES PARTIELLES
Définissons sur S une direction v que nous appellerons conormale à S. A cet effet, posons ,
N=
/
m
m
V k=l 2j [~ aik cos (n, i=l
Xi)]2
(112)
et définissons la direction v par les formules m
cos (v,
x,J= ~ ~
cos (n,
aik
Xi)
(k= 1, 2, ... , m).
(113)
i=l
La première formule (111) peut s'écrire encore m
au
~
peu) = N L.J
au
cos (v, x h) = N av '
aXk
k=l
et la formule de Green (110) devient en définitive ) [vL(u)-uL*(v)]d,;= )[N(v D
~~-u :~)+uvQJdS.
(114)
8
A noter que si m
~
ih =b aa. aXk i
L.J
(i=l, 2, ... , m),
k=l
alors Q = 0, l'opérateur L* (v) est confondu avec L (v) et l'on peut mettre L (u) sous la forme m
m
L (u)
=
~ aax . ~ i=l
~
aikuXk
+ cu.
k=l
Dans ce cas, l'opérateur L (u) est dit symétrique. Dans le cas général, l'opérateur différentiel L* n'est pas confondu avec L. On l'appelle adjoint de L au sens de Lagrange.
1-2-20. Méthode de Volterra. La résolution du problème de Cauchy pour équations du second ordre dans le cas où le nombre de variables indépendantes est supérieur à deux est bien plus compliquée. Nous avons donné les solutions explicites du problème de Cauchy pour l'équation des ondes avec des conditions initiales supportées par le plan t= 0 (tome II, [VII-1-9l). La méthode qui nous a permis d'obtenir ces solutions ne se généralise malheureusement pas. Dans ce numéro, on expose une autre méthode de résolution du problème de Cauchy pour équations à coefficients constants. Cette méthode, qui est une généralisation de celle de Riemann, repose, de même que cette dernière, sur un usage original de la formule de Green.- Cette méthode
1-2-20. M:eTHODE DE VOLTERRA
133
donne la solution du problème de Cauchy lorsque les conditions initiales sont supportées non seulement par le plan t = 0, mais aussi par des surfaces non caractéristiques. L'esprit de cette méthode est proche de celui des méthodes utilisées pour les équations à coefficients variables. Supposons que S est une hypersurface caractéristique de l'équation L (u) = 0 ou L (u) = f, où f est une fonction donnée. Soit 0) et par suite l'équation (33) définit bien un ellipsoïde. La résolution de l'équation en Â. nous donne .en tout point du corps trois racines strictement positives Â. = P5, P5 étant une fonction homogène de degré 2 par rapport à Pl' P2' Ps. Si l'on divise les deux membres de l'équation (33) par g2, les Pk se transforment en cos a.h, où cos a.h sont les cosinus directeurs de la normale à la surface de l'onde, et la racine p~ se transforme en P2. Donc, en chaque point on obtient trois vitesses possibles de propagation de l'onde dans une direction fixe. Les composantes (hl' h 2, ha) du vecteur de discontinuité se déduisent du système d'équations sans seconds membres qui définit les directions des axes de symétrie de l'ellipsoïde (33). Donc, si la direction est donnée, on aura en chaque point trois vecteurs de discontinuité deux à deux orthogonaux correspondant aux trois vitesses de propagation. Pour obtenir des ondes longitudinales et transversales, il est nécessaire et suffisant que l'un des axes de symétrie de l'ellipsoïde soit dirigé suivant la normale à l'onde correspondante. Sous cette condition, on aura une onde longitudinale et deux transversales. On admet que, la direction étant fixée, l'équation cubique ci-dessus admet trois racines distinctes. Dans le cas d'un milieu isotrope homogène, on a vu qu'une racine était double. Les cosinus directeurs de la normale à l'onde sont proportionnels à Pl' P2' Pa et par suite la condition nécessaire et suffisante revient à dire que pour une certaine racine Â. = PP5 les quantités (hl. h 2, h~) doivent être proportionnelles à (Pit P2' Ps), quel que soit Pk' c'est-à-dire quelle que soit la direction. En remplaçant ces quantités par les quantités proportionnelles h k dans le système homogène en hk , on obtient:
+
+
(api+ c"pl+ b"pl-ppi) Pl {
+
+ (c' + c") PIPi+(b' +b") PIP§ =0,
(c' +c") pip2+(c"pl+bpi+a"pâ--PP5) P2+(a' +a") P2P§=Ü, (b' +blf) pips+(a'
(34)
+ a") P~P3+ (b"pi+alfp~+cpi-PP~) Ps=O.
Si l'on tient compte du fait que quel que soit le choix de Pl' P2 et Ps, on doit obtenir à partir des équations (34) la même valeur pp~, on est conduit aux conditions suivantes pour les coefficients du potentiel élastique A :
=
c
d'où l'on déduit que PPi nales est
=
a
=
b
=
a'
ag~
+ 2a" =
b'
+ 2b" =
c'
+ 2c",
(35)
et la vitesse de propagation des ondes longitudi-
p=".1 V
Il
p'
200
CH. 1.
TH~OFJ.IE
DES~QUATIONS
D~RIV~ES
AUX
PARTIELLES
Les deux autres racines sont généralement distinctes et dépendent du choix de la direction de l'onde, c'est-à-dire du choix de Pk' Les égalités (35) nous fournissent cinq conditions pour les neuf coefficients qui figurent dans l'expression du potentiel éb.stique A. 1-3-7. Ondes électromagnétiques. Considérons les deux premières équations de Maxwell en milieu isotrope: crot H
=
ÂE
+
crot
BEt ,
E
=
-JA-Ht ,
(36)
où E et H sont les champs électrique et magnétique. c, la vitesse de la lumière, Â, le coefficient de conductibilité du milieu, B et ~, la constante diélectrique et la perméabilité magnétique. Les vecteurs E et H sont fonctions des variables indépendantes (xl' x~, Xa, t). En désignant leurs composantes par (el' e2' ea) et (hl' h 2 , ha), on peut mettre l'équation (36) sous la forme
J:. c
J:.. c
ôh t ôt ôh 2 ôt
+ ôes + ôel
ôXa
J:. iJh a + c
_
iJx 2
Dt
ôe 2 =0, ôXa
_~ ôes =0, ôXI
ôe z _ ôXI
(37)
ôet =0, ôX z
les points de suspension remplacent des termes ne contenant pas de dérivées des fonctions ek et h k • Nous avons affaire à un système de six équations du premier ordre à six fonctions. Numérotons ces fonctions comme suit: UI = el;
= e2;
Ua = ea;
= hl;
= hz;
= ha· En formant les expressions (5) et en écrivant l'équation (6), on obtient l'équation du premier ordre suivante pour les surfaces caractéristiques: U2
B
U4
u6
-Po
0
0
0
Pa
-Pz
0
-Po
0
-Pa
0
Pl
0
0
-Po c
c
0 Pa
-P2
B
c
-Pa
B
Pz
P2
-Pl
ua
0 =0.
~
-Po
0
0 0
c
~
0
-Pl
0
-Po
Pl
0
0
0
c
(38)
~Po c
Multiplions les trois premières colonnes de ce déterminant par
J:c
Po- Puis à la
première colonne ajoutons la cinquième multipliée par (-Pa) et la sixième multipliée par P2; à la deuxième colonne, ajoutons la quatrième multipliée par Pa et la sixième multipliée par (-Pl); à la troisième colonne ajoutons la quatrième multipliée par (-P2) et 'la cinquième multipliée par Pl- En développ ant ensuite ce déterminant suivant les éléments de la sixième, cinquième
I~3~7.
20t
ONJ)ES :eLECTROMAGN:eTIQUES
et quatrième ligne, on est conduit à l'équation
9+ pf
PIP2
P2PI
q+ P~
PaPI
PSP2
=0,
où q= e~ P~_ g2. c
(40).
Le développement de ce déterminant nous amène à l'équation (g2=pi+p~+pi),
q2(q+g2)=0
(41}
°
qui se scinde en deux équations. Si l'on égale la somme entre parenthèses à zéro,. on trouve que Po = et l'on obtient une onde stationnaire [1-2-14]. On examinera ultérieurement le second cas où q = 0, c'est-à~dire lorsque
qt c2 p2_g2=0 0 ,
(42\.r
ce qui conduit à une expression notoire de la vitesse de propagation de l'onde: V=
c Vef-L •
(43)·
Etudions maintenant la nature de la discontinuité. Désignons par (al' a 2 , as) les coefficients de discontinuité des dérivées des composantes du champ E, par (~l' ~2' ~s), les mêmes grandeurs relatives au champ H. Introduisons comme toujours les vecteurs de discontinuité a (al' a 2 , as) et li (~l' ~2' ~s). On a [EXh] = Pha, {
(k=O, 1, 2, 3;
[H Xh ] = Phli
X o=
(44)-
t).
Les trois premières équations (18) deviennent ici: (
1
~
~ POal + PS~2 - P2~3 = 0,
+ c
POa 2
+ PI~S
poas
+ P2~1 -
1
l cB
-
PS~l =
0,
PI~2 = 0,
ou, en désignant par n le vecteur unitaire de la normale au front d'onde dirigé dans le sens où WI > 0, -
(45):
eV a=li X n, c
WI
= 0(46).
le ~econd membre représentant le produit vectoriel des vecteurs li et n. De façon· analogue, les trois dernières équations (18) s'écrivent:
~V c
li= -a X n.
(47)-
De ces équations il s'ensuit immédiatement que les vecteurs a et li sont situés dans un plan tangent à l'onde et sont deux à deux orthogonaux. Supposons que devant le front d'onde, c'est-à-dire dans la région où wi > 0, le milieu soit au repos, ce qui s'exprime par la nullité de E et H. Les formules.
:202
CH. I. TH~ORIE DES ~QUATIONS AUX D~RIV:E:ES PARTIELLES
.(44) nous donnent les valeurs des dérivées de E et H sur le front d'onde:
E Xk = -Pk a ; H Xk = -Pk~.
(48)
'Considérons le développement taylorien de E et H limité aux termes contenant !les dérivées premières au voisinage du front d'onde. Comme E et H sont nuls :sur le front d'onde, on obtient, compte tenu de (48), les formules approchées :suivantes: 3
E,.., -a
3
2J
H ,.., -~ ~ Pk (Xk -xj"o»,
Pk (Xk-Xit°»;
k=O
'Où (xilO>, fonction 'ID I
k=O
x~o>tx~°l, X~OI) û.)l
est un point du front d'onde. En développant la en série de Taylor, on peut éc.rire, compte tenu de ce que
(x60>, x1°>,
x~o>,
xâo» = 0 : 3
(Xo,
û.)l
Xl'
X 2, X3)"""
2J
Pk (Xil -xk.° l ),
k=O
,et les formules précédentes deviennent (cf. [1-3-4]) E ,..,
-û.)l
(xo,
Xb X 2 '
x3)
a;
H ,..,
-û.)l
(xo,
Xl' X2' X3)
p.
(49)
-Ces formules approchées sont valables au voisinage de l'onde du côté qui est le siège d'un processus électromagnétique. Dans un milieu homogène anisotrope la quantité 8 ne doit plus être considé!l'ée comme un nombre, mais comme une matrice d'ordre 3. Cette quantité figure dans la relation qui lie le vecteur déplacement électrique au vecteur E (tome II, nV-2-12]). La quantité f.t est un nombre comme avant. Choisissons les axes de -coordonnées de telle sorte que la matrice 8 se réduise à la forme diagonale et soient 8 3 > 8 2 > 8 1 > 0 ses valeurs propres (tome 1111 , [II-2-1], [11-2-2]). Sous -ces conditions, les trois premières équations (37) s'écrivent:
!.L C
oel ot
8S
C
oh 2 oXa
oh a
_
+ ... =0,
OX2
a_ + oh oX
oh l oXs
+ . . . =0 ,
oe s +oh l _oh 2 ot OX2 oXl
+ ... =0,
~ oe 2 C
+
0t
l
-et au lieu de l'équation (39) on obtiendra l'équation
+
ql pi PlP2 PIPa PlP2 q2+ p~ P2PS PlPS
qs
P2Pa
=0,
(50)
+ P~
-où En posant
O. En procédant comme indiqué au (1-3-10) et en effectuant le changement
+
1fa~ u.; VI = Va u,. on ebtient deux équations séparées pour VI et v.: VI
aVI aXI
-ya
=
aVI aX2
+F (x 2) =0;
+ya
aV2 aXI
aV2 aX2
UlI'
(93)
-F (x2) =0.
(94)
Le changement de variables 2't']= --fa Xl +X 2 ,
2s=ytÏ %1 +X2; ramène le système à la forme
(95)
Cherchons la solution de (95) qui vérifie la condition initiale vllx1==o=v.lx1==o=O, Le.
VI
111=i=O;
v2111=~=0.
On déduit à partir de (95) i+11
J2,
,+11
J
F (t) dt;
F (t) dt.
211
Dans les anciennes variables indépendantes, on a 1
l'
F (t) dt;
",= );
Y
F (t) dt,
Y;Xl+ X• - VaxI+X2 et les formules (93) nous permettent de remonter aux solutions Ul' (92) qui vérifient les conditions initiales UI
Ix1==o =
U2
IX1=o =0.
Il est évident que cette solution est unique. Pour a = 0, le système (92) devient aUI _ aXI
au. =0; aX 2
aU2 ax}
+F (X2) = 0,
U2
du système (96)
214
CH. 1. TH:eiORIE DES :eQUATIONS AUX D:eRIV:eiES PARTIELLES
et sa solution qui vérifie les conditions (96) est: UI
= - 2,x~
F' (,x2);
U2
= -,xiF (,x2),
il est entendu que F (,x2) possède une dérivée seconde continue. Traitons enfin le cas où a = - b2 < O. En posant y
b,xl=,x;
,x2=y,
Vt=bUl+
~
) F(t)dt;
(97)
V2 =U2,
e
on écrit le système (92) sous la forme
+ aVI ôy
Ôl'2 =0' aV 2 ôy 'ô,x
ÔVI _ ô,x
=0 •
+
+
On remarque que VI V2i doit être une fonction régulière de z = ,x yi qui en vertu de (96) et (97) doit tendre pour ,x -- 0 vers la fonction réelle . y
~
(98)
) F (t) dt. e
-On peut affirmer que la fonction régulière z doit être prolongée par analyticité .au-delà de la droite x = 0 et par suite doit être analytique sur cette droite même {tome 111 2 , (1-24]). Donc, la fonction (98) et, par suite, la fonction F (y), doivent être analytiques pour y réel. En développant la fonction (98) en série entière de (y - Yo), où Yo est un réel quelconque, soit:
-+ ) y
00
F(t) dt=
~
ak (y-yo)k,
k=O
c
on obtient 00
VI+v 2i= ~ (-i)kak(z-iYo)1~
(z=,x+yi)
k=O
pour les z proches de iyo. Sachant
VI
et V2' on détermine
UI
et U2 grâce à (97).
Chapitre II
. PROBLÈMES AUX LIMITES
11-1. Problèmes aux limites pour une équation différentielle ordinaire 11-1-1. Fonction de Green relative à une équation linéaire du second ordre. Le présent chapitre sera consacré à l'étude des problèmes aux limites aussi bien pour les équations différentielles ordinaires que pour les équations aux dérivées partielles. Nous avons à maintes reprises rencontré de tels problèmes. Nous nous proposons de donner un exposé systématique de cette question. L'application de la méthode de Fourier à la résolution de problèmes aux limites de la physique mathématique nous a conduits à plusieurs reprises au problème aux limites suivant pour une équation différentielle ordinaire du second ordre contenant un paramètre: trouver les valeurs du paramètre  pour lesquelles l'équation sans second membre d
dx
[p (x) y']
+ [Âr (x) -
q (x)] y =
°
(1)
possède une solution non nulle sur un intervalle borné fini [a, b] et vérifiant aux bornes de cet intervalle des conditions aux limites homogènes: O-lY (a)
+ O-zy' (a)
=
°;
~lY (b)
+ ~zY' (b)
= 0,
(2)
et ~k sont des nombres donnés non tous nuls. On admettra que p (x), q (x) et r lx) sont des fonctions continues sur [a, b] et de plus que p (x) ne s'annule en aucun point de cet intervalle et possède une dérivée continue. Désignons la somme des termes de l'équation (1) ne contenant pas le paramètre  par OÙO-k
d
L (y) = Ch [p (x) y']- q (x) y. On appellera comme toujours valeurs propres, les valeurs du paramètre  pour lesquelles le problème homogène possède des solutions non triviales, et fonctions propres, ces solutions non triviales. Ces fonctions sont de toute évidence définies à un facteur multiplicatif constant près. Il est immédiat de voir qu'à toute valeur propre est associée une fonction propre et une seule. En effet, supposons
216
CH. II.
PROB~MES
AUX LIMITES
à l'inverse que l'équation (1) admet pour une valeur propre  deux solutions linéairement indépendantes vérifiant les conditions aux limites (2). L'intégrale générale de l'équation (1) vérifie alors les conditions aux limites (2). Ce qui est absurde, puisqu'on peut trouver une solution de l'équation (1) pour des valeurs initiales y (a) et y' (a) ne satisfaisant pas la première condition (2). Par les transformations élémentaires dont on s'est maintes fois servi (tome 111 2 , [V-8J, lVI-2-3J, lVI-3-2l) on montre que des fonctions propres CPI (x) et CP2 (x) associées à des valeurs propres distinctes sont orthogonales, soit: b
Jr (x) CPl (x) CP2 (x) dx = O. a
_Introduisons pour l'opérateur L (y) une fonction identique à celle qui a défini la déformation statique subie par une corde sous l'action d'une force concentrée (tome IV 1 , (1-1]). Le rôle de l'opérateur L (y) incombait à y". Pour arriver de façon naturelle aux propriétés de la fonction introduite, considérons l'équation d
L (y) = dx [p (x) y'] -q (x) y=
-1 (x)
(3)
et supposons que la fonction 1 (x) est nulle sur l'intervalle la, bJ tout entier sauf un petit intervalle [~ - e, ~ el, où ~ est un point fixe intérieur à [a, blet de plus que
+
~+s
J f (x) dx= 1.
(4)
5- S
Lorsque e tend vers 0, on obtient à la limite l'analogue d'une force concentréa au point x = ~. Sous la condition imposée à f (x) considérons la solution Ye (x) de l'équation (3) qui vérifie les conditions aux limites (2) (on admet que cette solution existe). En intégrant les deux membres de l'équation (3) par rapport à x et en tenant compte de (4), on obtient
p(x)y~(x)
J-.-J. X=S+8
~+8
q(x)y.(x)dx= -1,
ou à la limite lorsque e -+ 0, y'
(~ + 0) -
y' (s - 0) =
-
p
~~)
,
c'est-à-dire que la dérivée y' (x) de la solution y (x) doit subir au point z = ~ un saut égal à - p ~~) • Cette solution dépendra de
217
11-1-1. FONCTION DE GREEN
toute évidence du point de l'intervalle [a, b] qui sera pris pour s~ de sorte qu'elle sera fonction des variables x et On la désignera par G (x, s) et on l'appellera fonction de Green relative à l'opérateur L (y) sous les conditions aux limites (2). Les considérations précédentes nous conduisent à la définition rigoureuse suivante de la fonction de Green: on appelle fonction de Green relative à l'opérateur L (y) avec les conditions aux limites (2) une fonction G (x, s) vérifiant les conditions suivantes: 1) G (x, s) est définie et continue dans un carré k o défini par les inéquations a:::( x, s:::( b; 2) G (x, s), comme fonction de x, possède pour a:::( x < set S< x:::( b des dérivées premières et secondes continues et est solution de l'équation sans second membre L (y) = 0; 3) G (x, s), comme fonction de x, satisfait les conditions aux limites (2) ; 4) la dérivée de G (x, s) par rapport à x que nous désignerons par G' (x, s), subit un saut sur la diagonale du carré ko, c'est-àdire pour x = S, et de plus
s.
r {
G' (s + 0, s) -- G' (s - 0, s) = G'
(s, s+ 0) -
G'
(s, s-
0)
=
p
p
~6)
~6) ,
(5}
•
Les conditions (5) se traduisent par une seule, savoir: à l'approche de tout point x = S de la diagonale, aussi bien d'en haut, c'est-à-dire à partir du domaine S > x, que d'en bas, c'est-à-dire à partir du domaine S< x, la dérivée G' (x, s) doit prendre des valeurs bien définies et la différence de ces valeurs limites doit être égale à p ~s) . Dans chacun de ces domaines, la dérivée seconde par rapport à x s'exprime, en vertu de L (G) = 0, de la manière suivante: p (x) G" (x, s)
=
-p' (x) G' (x, s)
+ q (x) G (x,
s),
et par suite cette dérivée seconde prendra des valeurs bien définies à l'approche de la diagonale d'en haut ou d'en bas. Prouvons maintenant qu'il existe une fonction de Green et une seule satisfaisant toutes les conditions indiquées ci-dessus. On admettra que Îv = n'est pas une valeur propre, c'est-à-dire que l'équation L (y) = ne possède pas de solutions non ~riviales satisfaisant les conditions (2). Nous verrons dans la suite comment il faut modiest une fier la définition de la fonction de Green lorsque Îv = valeur propre. Construisons la solution YI (x) de l'équation L (y) = 0 en prenant pour valeurs initiales YI (a) et Y~ (a) des nombres vérifiant la première condition (2). La solution YI (x) et plus généralement toutes les solutions CIYI (x), où Cl est une constante arbitraire, satisferont la première condition (2). Il est immédiat de voir qu'il n'existe pas d'autres solutions vérifiant la première condition (2). En effet, s'il en existait une, Y (x), on aurait deux équations sans second
° °
°
218
CH. II.
PROB~MES
AUX LIMITES
membre en Ct l et Ct 2 : CtIYI (a)
+ (Z2Y~ (a)
= 0;
(ZiY
(a)
+
Ct 2 y'
(a)
= 0,
et comme on admet naturellement que l'un au moins de ces nombres est non nul, le déterminant de ce système devrait être nul, autrement dit le wronskien des solutions y (x) et YI (x) serait nul pour x = a, donc Y (x) et YI (x) seraient linéairement indépendantes, c'est-à-dire que Y (x) = CYl (x) (tome II, 11-1-1). De façon tout à fait analogue, supposons que C2Y2 (x), où C2 est une constante arbitraire, sont les solutions de l'équation L (y) = qui vérifient la deuxième condition (2). Le théorème d'existence et d'unicité nous dit que les solutions YI (x) et Y2 (x) sont définies sur l'intervalle la, b] et sont linéairement indépendantes. En effet, si elles étaient linéairement dépendantes, alors YI (x) vérifierait les deux conditions aux limites (2) et À = serait valeur propre, ce qui est contraire à l'hypothèse. La fonction G (x; s) est de la forme CIYI (x) pour x ~ et de la forme C 2Y2 (x), pour x ~ Il reste à choisir les constantes Cl et C2 de telle sorte que cette fonction soit continue en x = S et que sa dérivée subisse le saut signalé. Ceci nous conduit aux deux équations suivantes pour la détermination de Cl et C2 : ( CIYI Œ) - C2Y2 (s) = 0, ) 1 (6) l cIY; (s) - C2Y~ (s) = P (6) •
°
°
s,
s.
y;
Le déterminant [Y2 (s) Y~ (S) - YI (S) (S)] de ce système est nOIl nul, puisque les solutions YI (x) et Y2 (x) sont linéairement indépendantes. Donc on obtient des valeurs bien définies pour les constantes Cl et C2. Le wronskien des solutions YI (x) et Y2 (x) s'exprime (tome II, 11-1-11) à l'aide de la formule YI (x) Y~ (x) - Y2 (x) Y~ (x)
= P ~X)
,
où C est une constante non nulle. Quitte à multiplier l'une des solutions, par exemple YI (x), par c, on ramène la relation précédente à la forme p (x) [YI (x) Y~ (x) - Y2 (x) Yi (x)] = 1De cette relation il s'ensuit immédiatement que le système (6) a pour solution: Cl = Y2 (s) et C 2 = YI (S) et la fonction de Green G (x, S) est définie comme suit: G (x, S) = { YI (x) Y2 (s) Y2 (x) YI (S)
(x~ S),
(7)
(x~6)·
On vérifie sans peine qu'elle remplit toutes les quatre conditions. Son unicité résulte directement des raisonnements précédents.
R~DUCTION
11-1-2.
INT~GRALE
A UNE .eQUATION
219
11-1-2. Réduction à une équation intégrale. Considérons l'équation avec second membre L (y)
=
d
- f (x),
(IX [p (x) y'] - q (x) y =
(8)
où f (x) est une fonction donnée, continue sur un intervalle [a, b]. On cherchera la solution de l'équation (8) qui vérifie les conditions aux limites (2). Cette solution est unique. En effet, s'il en existait deux, leur différence satisferait l'équation sans second membre L (y) = 0 et les conditions aux limites (2), c'est-à-dire que 1.. = 0 serait valeur propre. Assurons-nous que la solution unique de l'équation (8) qui satisfait les conditions aux limites (2) est donnée par la formule b
Y (x) = ) G (x, s) f (s) d
s·
(9)
a
L'analogue de cette formule (tome IV l , (1-1]) admet une signification mécanique simple, à savoir que si l'on connaît la déformation statique due à une force ponctuelle, on peut au moyen d'une intégration obtenir la déformation statique due à une force continûment distribuée. Passons à la démonstration du fait que la fonction définie par la formule (9) satisfait (8) et les conditions aux limites (2). Partageons l'intervalle d'intégration en deux, compte tenu du saut subi par la fonction de Green: x
b
y= ) G(x, s)/ (s)ds+ ) G(x, s)f(s)ds. a
x
Une dérivation par rapport à x nous donne x
y' ==
) G' (x,
s) f (6) dx + G (x, x-O) / (x) +
a
b
+ ) G' (x, s) / (s) ds -
G (x, x + 0) f (x)
x
ou encore x
y' =
) G' (x, s) a
b
f (6) ds + ) G' (x, x
b
s) / (s) ds = ) G' (x, s)
f (s) ds, (10)
a
+
puisque la fonction de Green est continue, i. e. G (x, x 0) = = G (x, x - 0). Des formules (9) et (10) et du fait que G (x, s) satisfait les conditions aux limites (2), il s'ensuit aussitôt que la
CH. II. PROBL:E:MES AUX LIMITES
220
fonction (9) vérifie également ces conditions. Pour vérifier l'équation (8) dérivons encore une fois y' par rapport à x. Des transformations peu compliquées nous conduisent à l'expression b
y" = ~ G"(x, s)f(s)ds--!--[G'(x,x-O)-G'(x,x+O)]f(x), a
et de (5) il s'ensuit que b
y" =
JG" (x, s) f (s) d s-
;
~:~
•
(11)
a
En portant dans (8) les expressions (9), (10) et (11) de y, y' et y", on trouve b
~ L(G)f(s)ds-f(x)= -f(x), a
c'est-à-dire que l'équation (8) est bien vérifiée puisque la fonction G (x, s) est solution de l'équation sans second membre L (y) = O. Signalons que les formules précédentes nous disent que la fonction y définie par (9) possède des dérivées premières et secondes continues sur l'intervalle [a, b]. On est donc conduit à la proposition suivante: Si  = n'est pas une valeur propre pour l'équation (8), alors cette équation admet pour toute fonction f (x) définie et continue sur [a, bl une solution unique vérifiant les conditions aux limites (2) et cette solution est donnée par la formule (9). Cette proposition peut encore s'énoncer: Pour toute fonction f ex) continue donnée, la fonction (9) possède des dérivées premières et secondes continues et vérifie l'équation (8) avec les conditions aux limites (2). A noter que si y (x) est une fonction quelconque possédant des dérivées première et seconde continues sur un intervalle fa, blet vérifiant les conditions aux limites (2), alors en la portant dans l'équation (8) on peut construire une fonction f (x) continue. D'après ce qui a été prouvé plus haut, la fonction y (x) s'exprime au moyen de f (x) par la formule (9). Donc, les formules (8) et (9) établissent une correspondance biunivoque entre deux classes de fonctions: la première étant composée des fonctions y (x) possédant sur fa, b 1 des dérivées première et seconde continues et vérifiant les conditions aux limites (2); la seconde, des fonctions f (x) continues sur fa, bl. Le passage de y (x) à f (x) est réalisé par la formule (8), le passage de f (x) à y (x), par la formule (9). On voit donc à la lumière de ce qui vient d'être dit qu'il est possible de ramener le problème aux limites à une équation inté-
°
221
II-f-3. SYM:eTRIE DE LA FONCTION DE GREEN
grale. En effet, en mettant l'équation (1) sous la forme L (y) = -Âr (x) y, on constate immédiatement en vertu des résultats établis que l'équation (1) avec les conditions aux limites (2) est équivalente à l'équation intégrale b
Y (x) =
Â
~ G (x, 6) r (6) y (6) d6.
(12)
a
De façon tout à fait analogue, l'équation avec second membre d dx [p (x) y']
+ [Âr (x) .- q (x)] y =
F (x)
(12 1 )
avec les conditions aux limites (2) équivaut à l'équation intégrale b
y (x) = FI (x)
+ Â ~ G (x, 6) r (6) y (6) d6,
(12 2 )
a
où b
FI (x)
=
-
~ G (x, 6) F (6) d6, a
la solution étant dans les deux cas une fonction continue y (x). 11-1-3. Symétrie de la fonction de Green. La formule (7) définit la fonction de Green non seulement pour x Ela, bl, mais aussi aux bornes x = a et x = b de cet intervalle, c'est-à-dire sur le carré fermé k o : a ~ x, 6~ b, et cette formule entraîne directement la symétrie de la fonction de Green sur k o : G (x, 6)
=
G (6, x).
(13)
Donnons une u utre démonstration de la symétrie de la fonction de Green, une démonstration fondée sur une idée utilisée dans des cas plus généraux. L'identité suivante est de vérification immédiate: d
uL (v) - vL (u) = dx [p (x) (uv' - vu')].
(14)
Dans cette identité, u (x) et v (x) sont deux fonctions arbitraires possédant des dérivées. premières et secondes continues. Posons u = G (x, 61) et v = G (x, 62) dans (14) en admettant pour fixer les idées que 61 < 62" En intégrant sur les intervalles [a, 61], [617 62) et [62' b] et en tenant compte du fait que la fonction de Green est
222
CH. II. PROBL1!:MES AUX LIMITES
solution de l'équation L (y) = 0, on obtient [p (x) (G (x, ~l) G' (x, ~2) -
G (x, ~2) G' (x, ~l»]~_~l = 0,
[p (x) (G (x, ~l) G' (x, ~2) -
G (x, ~2) G' (x, ~l»]~ i~ = 0,
[p (x) (G (x, ~l) G' (x, ~2) -
G (x, ~2) G' (x, ~l»]~=t~2 = O.
En ajoutant ces trois égalités et de par la continuité de la fonction de Green et la discontinuité de sa dérivée première, on trouve G (~l' ~2) - G (~2' ~l) = = [p (x) (G (x, ~l) G' (x, ~2) -
G (x, ~2) G' (x, ~l»]~=~'
(15)
Il est immédiat de vérifier que la différence du second membre s'annule pour x = a et x = b. En effet, la fonction de Green satisfait la première condition (2), soit
alG (a, ~l) alG (a, ~2)
+ a G' (a, + a G' (a, 2
~l) =
2
~2)
0,
= 0,
et comme on admet naturellement que les constantes al et a 2 ne sont pas simultanément nulles, le déterminant de ce système homogène doit être nul, autrement dit, la différence ci-dessus s'annule bien pour x = a. On démontre de même qu'elle s'annule pour x = b et la formule (15) établit la symétrie de la fonction de Green. On peut envisager des conditions aux limites plus générales que les conditions (2), c'est-à-dire des conditions dans lesquelles figurent les valeurs prises par la fonction et sa dérivée première aux bornes de l ' intervalle [a, b]:
alY (a) ~lY (a)
+ a y' (a) + aaY (b) + a,.y' (b) = 0, + ~2Y' (a) + ~aY (b) + ~,.y' (b) = O. 2
Tous les raisonnements précédents, excepté la démonstration de la symétrie de la fonction de Green, restent en vigueur. Pour que la démonstration de la symétrie de la fonction soit valable, il est nécessaire et suffisant que soit remplie la condition suivante: p (b)
al' a 2 R
R
Pl'
P2
=p (a)
a 3 , a.. R
R
P3'
P"
•
Nous glisserons sur la démonstration de ce fait. On vérifie immédiatement que la fonction de Green est symétrique pour des conditions aux limites périodiques: y (a) = y (b); y' (a) = y' (b) si p (a) = = p (b), c'est-à-dire si la fonction p (x) est périodique. Signalons que si les autres coefficients q (x) et r (x) sont périodiques, alors le problème aux limites avec les conditions périodiques mentionnées revient à chercher les valeurs du paramètre 'A pour lesquelles l'équation (1) possède une solution périodique.
11-1-4. VALEURS ET FONCTIONS PROPRES
223
II-1-4. Valeurs et fonctions propres du problème aux limites. Etant donné qu'on a ramené le problème aux limites à une équation intégrale, on peut mettre en œuvre la théorie générale des équations intégrales et établir ainsi des propositions relatives aux valeurs et fonctions propres du problème aux limites. Commençons par le cas r (x) == 1, où l'équation (1) est de la forme
ëfX [p (x) y'] + (1.. - q (x» y = 0, d
(16)
les conditions aux limites étant telles que la fonction de Green est symétrique. L'équation intégrale (12) est une équation à noyau symétrique. Elle admettra des valeurs propres réelles, et les fonctions propres associées aux valeurs propres distinctes seront orthogonales. On a vu au (II-1-11 qu'à toute valeur propre était associée une fonction propre et une seule. Ceci a été prouvé pour des conditions aux limites du type (2). Si les conditions aux limites sont périodiques, à toute valeur propre il peut correspondre deux fonctions propres au plus, car l'équation (16) ne possède que deux solutions linéairement indépendantes. Prouvons encore que le noyau G (x, ~) de l'équation (16) est un noyau complet, c'est-à-dire qu'il n'existe pas de fonction continue j (x) non identiquement nulle qui soit orthogonale au noyau. Si par absurde on suppose qu'une telle fonction existe, c'est-à-dire que b
~
G(x, s)j(s)ds=Ot
a
on constate alors que la fonction (9) d'un côté est identiquement nulle et de l'autre est, en vertu de ce qui a été démontré plus haut, solution de l'équation avec second membre (8), ce qui est impossible. Le noyau étant complet, il existe une infinité de valeurs propres (tome IV l' (1-421). Soient 1.. n (n = 1, 2, ... ) les valeurs propres de l'équation (16), c'est-à-dire du problème aux limites envisagé, et
> O. Les résultats du (tome IV l' [1-44]) nous disent que le problème aux limites (2) pour l'équation (1) se ramène alors à une équation intégrale à noyau symétrique. En particulier, toute fonction satisfaisant les conditions aux limites (2) et possédant des dérivées première et seconde continues sur [a, b), se développe en une série de Fourier régulièrement convergente suivant les fonctions propres du problème: 00
f (x) = ~ cnCPn (x),
(18)
n=1
dont les coefficients sont définis pas les formules b
Cn
=
Jr (x) f (x) CPn (x) dx.
(19)
a
Pour prouver cette proposition, on observera qu'en vertu du (IJ-1-2] on a b
f(x)= -
JG(x,
s)Lff(6))ds·
a
Or, de toute évidence L If (6)] = - V r (6) h Œ), où la fonction h (6) est continue sur· la, bl, puisque r (6) >0. Donc, la fonction V r (x) f (x) se représente par le noyau de l'équation intégrale
II-1-S. SIGNE DES VALEURS PROPRES
225
symétrique : b
V r (x) t (x) = ) G (x, 6) V r (x) r (s) h (s) dS,
(20)
a
et les raisonnements du (tome IV!, [1-44]) nous donnent immédiatement le théorème de décomposition formulé ci-dessus. Comme plus haut on peut prouver la fermeture du noyau et par suite l'existence d'une infinité de valeurs propres. En reprenant les raisonnements du [11-1-1) pour le cas où f (x) possède une dérivée première continue et une dérivée seconde continue par morceaux et en se rappelant que le théorème 2 du (tome IV h [1-31]) est valable pour le cas où une fonction est représentable par un noyau à l'aide d'une fonction h (x) continue par morceaux, on s'assure que le théorème de décom· position est valable pour le cas où la fonction f (x) qui satisfait les conditions aux limites (2) possède une dérivée première continue et une dérivée seconde continue par morceaux. Dans la suite, on signalera les cas où les théorèmes de décomposition restent valables pour une dérivée première continue par morceaux. 11-1-5. Signe des valeurs propres. On admettra pour alléger l'écriture des formules que r (x) == 1. Les raisonnements s'étendent sans peine au cas général. Exhibons tout d'abord une formule exprimant les valeurs propres au moyen des fonctions propres. Soient Ân les valeurs propres et ~ (x) + q (x) q>~ (x)] dx.
(22)
a
Supposons que p (x) > O. Si l'on admet de plus que q (x);;::: 0 sur [a, b], alors la formule (22) nous dit que les valeurs propres sont toutes strictement positives. Supposons maintenant que q (x) est une fonction continue arbitraire et soit m son minimum sur [a, b], c'est-à-dire que q (x);;::: m, x E [a, bl. De la formule précédente, il s'ensuit immédiatement b
Ân~
) P (x)
q>~2 (x) dx+ m~m.
a
Donc, dans le cas envisagé, seul un nombre fini de valeurs propres peuvent être strictement négatives. Supposons maintenant que les conditions aux limites sont de la forme: yi (a) - hly (a) = 0; yi (b) h 2 y (b) = 0, (23) où hl et h 2 sont des constantes strictement positives. Le second terme de la formule (21) est alors strictement positif et l'on s'assure comme plus haut que sous les conditions aux limites (23) et pour q (x) ~ 0 les valeurs propres sont toutes strictement positives. Si toutes les valeurs propres sont strictement positives ou si seul un nombre fini d'entre elles sont strictement négatives, alors est vrai le théorème de Mercer et l'on peut développer le noyau en une série absolument et uniformément convergente:
+
00
G (x, ~) = ~
fPn
(Xi:
n
(~)
•
(24)
n=l
Cette égalité nous permet d'étendre le théorème de décomposition suivant les fonctions propres qui a été prouvé au (11-1-4] à une classe de fonctions plus vaste, plus exactement, supposons que f (x) est continue, possède une dérivée continue sur l'intervalle [a, b] t01,lt entier sauf en un point x = c où elle subit un saut: f' (c 0) - f' (c - 0) = k,
+
et une dérivée seconde continue par morceaux. Supposons par ailleurs que f (x) satisfait comme toujours les conditions aux limites (23). Composons la différence:
f (x) -
k
p (c)
G (x, c),
qui possède une dérivée continue sur l'intervalle la, b] tout entier. Cette -différence est justiciable du théorème de décomposition sui-
II-t-6. EXEMPLES
227
vant les fonctions propres. D'autre part, en vertu de la formule (24) le terme G (x, c) peut être développé suivant les fonctions propres, donc la fonetion f (x) se décompose en une série de Fourier absolument et uniformément convergente suivant les fonctions propres. Ces raisonnements sont valables de toute évidence dans le cas où la dérivée f' (x) présente un nombre fini de sauts dans [a, bl. Ils sont identiques à ceux utilisés au (tome II, [VI-2-8J) pour améliorer la convergence des séries de Fourier. 11-1-6. Exemples. 1. Considérons l'équation y"
avec les conditions aux limites y (0) de Green s'écrit G (x,
s)= {
+ 'Ay = =
y (1)
0
=
O. On a L (y)
s) x
(x ~ ~),
(1-x)s
(s~x).
(1-
=
y" et la fonction
Les valeurs et fonctions propres se définissent sous forme explicite 'An
=
n 2n 2 ;
CPn (x)
= 11 2 sin nnx
(n
=
1, 2, •••),
et les fonctions satisfaisant toutes les conditions du numéro précédent sont justiciables du théorème de décomposition suivant ces fonctions propres. Il est possible d'élargir les conditions d'applicabilité du théorème de décomposition, mais nous glisserons là-dessus. 2. Considérons l'équation de l'exemple précédent avec les conditions aux limites y (0) = y' (1) = O. On a
et les valeurs et fonctions propres sont: 31;2
'An =(2n+ 1)24'""";
-
CPn (x) =-.12 sin (2n+ 1)
n
T
x.
3. Considérons encore l'équation de l'exemple 1 avec les conditions aux limites y (0) = 0; y (1) hy' (1) = O. Composons la fonction de Green. Prenons deux solutions de l'équation y" = 0, l'une vérifiant la première condition aux limites, l'autre, la seconde, soit: YI (x) = x; Y2 (x) = (1 h) - x. Les raisonnements du [1I-1-1J nous· conduisent à la fonction de Green suivante: '
+
+
G (x,
s)= {
1~~~s
x
(x~s),
1+h-x ~ 1+h
(S ~ x).
Les valeurs propres sont toutes strictement positives et en posant  = fJ.2 on s'assure sans peine que les valeurs propres correspondantes fJ. se déduisent à· partir de l'équation tg fJ. + hfJ. = 0 et les fonctions propres seront Cn sin fJ.nx, ' la constante Cn étant déterminée à partir des conditions de normalisation de ces fonctions... 4. L'étude des vibrations d'une membrane circulaire fixée par son contour nous a conduits au problème aux limites suivant: trouver les valeurs du para15*
CH. II. PROBL1l1MES AUX LIMITES
228
mètre  pour lesquelles l'équation
1l'+.!...Y'+ z
(Â- Zlni ) Y=O
(25)
admet une solution vérifiant les conditions aux limites sup 1 y (0) 1 < 00 et y (l) = 0, (*) n étant un entier positif. Ce problème aux limites diffère de ceux étudiés jusqu'à maintenant par le fait que les coefficients de l'équation présentent un pôle en z = et au lieu d'une condition bien définie en x = on exige simplement que la solution soit bornée au voisinage de x = O. La solution de l'équation (25) prend une valeur finie pour x = O. Nous avons rencontré plusieurs fois de tels problèmes aux limites. En multipliant les deux membres de (25) par x, on met l'équation sous la forme usuelle:
°
°
2
d (xy')+ ( ÂX--;n ) y=O, dx
(26)
°
où n > est un nombre donné. La fonction de Green se définit comme précédemment sauf qu'en x = au lieu d'une condition aux limites on exige que la fonction de Green soit finie. L'équation L (y) = est une équation d'Euler (tome II, [II-1-18]) qui admet les solutions linéairement indépendantes x n et x-no La fonction de Green devant être finie en x = 0, on doit la prendre de la forme Clxn sur l'intervalle [0, S]. Sur l'intervalle [S, 1] on doit former une combinaison linéaire des solutions zn et x- n qui s'annule pour x = 1, autrement dit, la fonction de Green doit être de la forme C2 (zn- x- n ). Les constantes Cl et Cl! se déterminent comme toujours à partir des conditions de continuité de la fonction de Green et de saut de sa dérivée première en x = S. On obtient en définitive la formule suivante:
°
°
(x ~ ~), (~~
x).
L'équation avec -second membre L(y) = -f (x) possède une solution unique satisfaisant les conditions aux limites (.), soit 1
y(X)=) G(x,
S)!(S)ds.
o Les raisonnements du [II-1-5] nous disent que les valeurs propres sont toutes strictement positives. En posant  = f.t2 on déterminera les valeurs propres à partir de l'équation transcendante J n (f.t) = 1 et les fonctions propres seront 9>n (x) = cnJn (f.tn x ). Pour n = 0, l'équation 'L (y) = possède deux solutions linéairement indépendantes: YI (x) = 1'et Y2 (x) = ln x. La fonction de Green est donc: -ln x~ (27) G(x,s)= { -lnx, s~x.
°
s,
s,
A noter que la ,formule (9) nous donne ici: 1
y
(x),-
J
x
~
0
G(x, s)!(s)d;=-lnx ) f(S)
1
J
d~- !(s)ln~ dS, x
229
11-1-7. DISTRIBUTION DE GREEN
et l'on vérifie immédiatement que cette fonction satisfait l'équation L (y) = = -1 (x) avec les conditions aux limites (-\10). La forme de l'é~uation (26) indique que r (x) = x. En ramenant le problème aux limites à une equation intégrale on obtient le noyau G (x, 6) V6X qui est continu sur le carré tout entier, y compris en x = 6 = O. Les fonctions propres de cette équation intégrale sont q>n (x) = = cn V; Jo (....nx ) et l'on obtient la série de Fourier absolument et uniformément convergente suivante: 00
G (x,
s) -ys x =
~
Cj)n
(X~:n
(6) •
n=1
Une division par
-V 6x
nous donne
~ c;J0 ( ....n x ) Jo 0 et soit alla partie ·de S définie par (17) _ Montrons que l' on peut choisir dl assez petit pour que pour toute position de M l'on ait au voisinage de No
1l l I-t ~N) dS 1< : .
(45)
al
()n a
1ll I-t~N) dSI~ ll P~
dl;dll,
(46)
a~
al
-OÙ a~
est un disque de centre No et de rayon dl et Pl la longueur du ·projeté MINI de MN sur le plan tangent en No- Supposons que M 'est intérieur à la boule de centre No et de rayon dl- Le point Ml ;appartient alors au disque a; et le disque a~ de centre Ml et de rayon 2dl du plan Œ, 11) contient entièrement le disque de sorte ~que, en vertu de (46),
a;,
1ll
I-t
~N)
dS
1~A ll
d~ld'YJ
2:-r 2d l
=
A
Pl~2dl
al
Jl 0
Pl
d~l dit
=
4Jtd I A.
0
~ : et l'on obtient la majoration (45) quelle que soit la position du point M sur la boule de centre No et de rayon dl' Mettons la fonction (44) sous :la forme u (M) = UI (M) u 2 (M),
in reste à fixer dl de telle sorte que l'on ait 4Jtdi A
+
UI
(M)
=
ll
I-t
~N)
dS ;
I-t (N) dS r
'
al
·la fonction U 2 (M) est continue en No, quant à la fonction U (M), -on démontre qu'elle l'est aussi en procédant exactement comme
II-2-4. DÉRIVEE NORMALE DU POTENTIEL DE SIMPLE COUCHE
283
au [II-2-21 pour la fonction (29). On a donc la proposition suivante: le potentiel de simple couche (44) est défini et continu sur l'espace tout entier. Comme au [II-2-2l on démontre que u (M) -+ 0 lorsque M -+ 00. 11-2-4. Dérivée normale du potentiel de simple couche. Soit no la normale extérieure en un point No de S. Supposons que MEtS et formons la dérivée de la fonction (44) suivant no. Le facteur .!.r étant le seul à dépendre de M, on peut dériver sous le signe d'intégration:
ai-
au (M) = \ \ (N) _ r dS = \ \ (N) cos2 'i' ano J J Il ano J J Il r
dS.
(47)
s Signalons la différence qui existe entre cette intégrale et l'intégrale (19) qui définit le potentiel de double couche. no Dans l'intégrale (19), l'angle • q> est l'angle de r et de n, 1 vecteur unitaire de la normale extérieure en N qui est le point variable d'intégration, alors que dans l'intégrale (47), l'angle 'l' est l'angle de r et de no, vecteur unitaire de la normale extérieure au point fixe No. Dans les deux cas, le vecteur r est de même sens que Fig. 7 S
~
MN (fig. 7). Montrons que l'intégrale (47) existe dans le cas où M est confondu avec le point No. Mettons l'intégrale (47) sous la forme
U
ft (N) cO: 1Po dS =
où r o = NoN et Nous montrerons l'intérieur ou de de (47) tend vers mules
6
U
ft (N) cos (~r no) dS,
(48) -~
l'angle 'l'o = (r o' no) est l'angle de NoN et no. ensuite que lorsque le point M tend vers Node l'extérieur de S suivant une normale, la dérivée des limites définies qui sont données par les for-
(auô~o)L = ~ ~ Il (N) cO~6'i'o dS + 2nll (No), s
(aUa~: 0») e = ) ) ~ (N) co;~'I'o dS - 21[~t (No), s
(49)
284
CH. II. PROBUlMES AUX LIMITES
où comme dans [11-2-2] les indices i et e signifient que M tend vers No de l'intérieur et de l'extérieur de S.
Dans le système de coordonnées locales d' origine No, le vecteur no est de même sens que l'axe N oZ. On désignera comme plus haut les coordonnées de M par (x, y, z) et celles de N dans le système local par (6, 11, ~). En introduisant comme plus haut la partie (Jo de la surface S, on peut écrire l'intégrale (47) sous la forme ) ) f1 (N)
\3 ZdS.
(50)
0'0
Si M
== No,
alors z = 0 et l'intégrale devient
) ) f1 (N)
r~t
dS = ) ) f1 (s,11) O'~
al
rg ~~;(n~) Z) ds d11,
OÙ ~ = ~ (S,11). On s'assure immédiatement que cette intégrale a un sens en vertu de (15) et de l'inégalité r o ;:::: Po. On a donc prouvé l'existence de l'intégrale (47) pour les M situés sur S. Passons à la démonstration des formules (49). Composons la différence de l'intégrale (47) et du potentiel de double couche de même densité f1 (N) :
a~~~)
_ w
(M) =
) ) I-t (N)
cos 'P~cos
qJ
dS.
(51)
s Cette intégrale a un sens si M ~ S ou si M ==. No. Montrons que cette différence reste continue lorsque M traverse la surface S en No. Il suffit à cet effet de prouver comme dans les numéros précédents que cette intégrale étendue à la petite portion (J1 de S, définie par (17), peut être rendue aussi petite que l'on veut en module. On supposera dans les estimations ultérieures que M est situé sur la normale à S en No, c'est-à-dire que x = y = 0 en coordonnées locales. On a alors cos 'P-;cos qJ =
r
__ ~ cos (n, Y)- ~ 3 Z (cos'6'o-1). (52) r cos (n, X) -4 r r
En tenant compte de (15) et de 1 si::::;; Po,
1 11 1::::;; Po;
1
r> Po > 2" r ;
1~ -
z 1::::;; r ::::;; 2po, ~
où Po= Vs 2 + 11 2 est la longueur du projeté de MN sur le plan XN oY, on obtient 1 cos 'P-cos qJ
r2
1
~ ~
b1 2-a.'
Po
où bl est une constante. Compte tenu de (22), il vient
285
II-2-5. DERIVEE NORMALE DU POTENTIEL DE SIMPLE COUCHE
1~ )
~ (N) cos 'P~COs cp dS 1~
)
~ ::~~ ds d11 =
Po:e::;;;d 1
al
0 2Jt dl
= 2Ab1
~ d~o_:
0 donné, on peut exhiber un 11 > 0 ne dépendant pas de la position de No sur S, tel que t (0 (M) - (0 (No) 1 ~ e si 1 MN o 1 ~ 11, M étant situé sur la normale à S en No. Figeons un dl tel que b2d~ ~ : et mettons (0 (M) sous la forme
==
(0
(M)
=
(01
(M)
+
(02
(M),
où {Ù1(M)=)) ~(N)cos'i'~cos: dl et r o >: dl' D'autre part, pour toute position des points N et No sur S, les modules des quantités YI, ~ sont inférieurs au diamètre de la surface S, c'est-à-dire à la distance des points les plus éloignés de S. On a par ailleurs , r - r o 1 ~ 1 z 1 et
s,
..!- __ 1 I-Ir-r0 1 (_1 +_1_+_1_) ~ r~r r~r2 ror 3 -..-: : : : / r~ r 3 -
1z 1
31 z 1 •
dt
'
r3
1z 1
~(j3' 1
et en vertu de la formule (56) 1
[
COS
,!,-cos
(60)
TO
1
Appelons 0'1 la portion de la surface S, découpée par le cylindre circulaire ayant pour axe la normale à S en No et pour base le disque de rayon 2To l ' Décomposons l'intégrale étendue à S en deux, l'une sur 0'1' l'autre sur S"'O'I:' J
1
= ~ ~
~
~
1
(JI
J2=
cos (ro, no) 1 as
1 cos (rI' nd
0
., (61)
\ ~ 1 cos (~r nt>
~"(Jl
L'introduction du produit scalaire nous permet d'écrire: cos (rll' no) rl,nl ro·no cos (ri' nt) r 12
ri
r~
_
rl,nO -
-
ro,no 3
Tg
+ rl,nl -
Tl
rt,no 3
rI
+ru .no (_1__..!-) 3
rI
3'
ro
«)Ù, comme toujours, no et ni sont les vecteurs unitaires de la normale extérieure -en No et NI' De ce qui précède il s'ensuit que cos (ro' no) 1::::::::: I r l,nO-rO,n61 ,/ cos (~j' nI)
+ + Irl,nl-rl,nol + 1ro,no1 1 Tt1 - rt1 1. Tt
T~
-.:::
Tf
(ü2)
Majorons séparément les divers termes: 1 rl,n l -rI"no 1 = 1ri' (ni -
no) 1 ~ rll ni - no 1.
En considérant le triangle de côtés no et ni' on obtient 1 nI - no 1 ~ '6', où '6' .est l'angle de no et nI' La condition (3) nous permet d'écrire: 1 rl,nl - rl,no 1 ~ arlr~,l' 'où a est une constante, Par ailleurs Irl,no -ro,no
1
=
1 (rl-rO),no 1
=
Iro.l,no
1
=
1 ~ll,
,où ~l est une coordonnée de NI dans le système local d'origine No. On a grâce
.à (15)
Irl,n O -ro,no 1 ~ crJ,ta..
Si enfin le point d'intégration N est assez proche de No, alors en vertu de (15) .on a 1 ro'no 1 = 1 ~ 1 ~ cr~+a., Mais, comme pour l'inégalité (59), on peut admet-
D2RIV~E
11-2-6. VALEUR DIRECTE DE LA
NORMALE
289
tre que l'inégalité précédente est valable pour toutes les valeurs de TO' En portant les majorations obtenues dans (62), on trouve cos (rOt no) 1~
cos (ri' nt)
T8-"':::
ri
Clrlr~, 1 + ctr~~f
_
":::::::
~
~
1
+CIT 0
(1 +-+-, fI) ~~ ~~ ~~
+CX Irl-rol
(63)
-3
où Cl = max {a, cl. On déduit à partir du triangle NoNIN que
Tt
+ To.
1
~
TO'
Or, en intégrant !ur (S ",,",0'1) on a ro. 1 ~ ~ et par suite rI ~ ~ • En se servant de ces inégalités et aussi de l'inégalité 1rI - ro 1 ~ ro, 1t on peut mettre (63) sous la forme cos (rI' nI) r12
1 (_ 0, 1 r21
~ c rcx
-...:::
1
cos (ro' no) 1~ r02 -"':::
+~ + r3 1
r1+ar1-CX 0 0, 1 r 3r
+
01
~ C1r~,
1
(
r 1+cx r 1- a 0 0, 1 r 2r 2
+
r 1+a r 1- a
01
1
4
2
4
0
O. 1 ) ~
r3 r01
-...:::
8)
19c1r~
-;:2+-;:2+-;:2+7=2+7=2 = o 0
0
0
0
r2' 0
1
•
En revenant à la deuxième formule (61), on obtient
arr
dS
IJ21~corO',1JJ
ri'
(64)
8"-01
où Co = 19c1' On a considéré un cylindre de rayon de base 2rO,I' Prenons maintenant un cylindre de même axe et de rayon fixe ~. Ce cylindre découpe sur S une portion 0'0 qui est comprise dans 0'1' On a
~ J ~f =
8"01
D~ns
la deuxième intégrale ro
JJ ~: + JJ ~: .
0 0 "al
~ d/3,
rr
JJ
8"-00
donc 9
dS
r~ ~ d2' 1SI,
8"00
où 1Siest l'aire de S. On peut ramener l'intégration sur 0'0""0'1 à une intégration sur le plan tangent en No puis à l'aide des estimations habituelles
( ro ~ Po et cos (n, Z) 2n d/3
JJ ~: ~ J J 0 0"-01
~ -{-) 2po ~~o
obtenir:
de
4n (ln
~
-ln 2ro•
0 2rO• 1
En portant ceci dans (64) on trouve J2~Alr~, 111nrO.ll+B1r~.1' 19-01017
1) .
CH. II. PROBUlMES AUX LIMITES
290
et BI sont des constantes. Si 0 < ration par une majoration de la forme:
OÙ Al
J 2 ~ A2
~
< cx on peut remplacer cette majo-
rg. 1.
Passons à la majoration de JI. On a JI
~ ) ) Icos
(:r
dS
nl)1
+)
al
Icos (~r no)! dS.
(65)
al
Les majorations habituelles nous donnent:
)) ICOS(~r
no)/
dS=))
al
2n 2ro• 1
':3'
dS~C)
)
0
0
al
=
A3r~.
l'
1
où A 3 est une constante. Pour évaluer la première intégrale (65) traçons la
sphère de centre NI et de rayon 4rO.l en remarquant que 4rO.l < ~ . Cette sphère découpe sur S une portion (J2 contenant (JI. La portion (J2 est définie par une équation explicite dans le système de coordonnées locales de centre NI et on peut lui appliquer les estimations habituelles en ramenant l'intégration au plan tangent à S en NI. Le domaine d'intégration est une partie du disque de centre NI et de rayon 4ro.l. En intégrant sur le disque tout entier, on obtient la majoration:
rr
JalJ
1cos(rl;
nI) 1 dS ~
rI
r r 1cos (r~, nt>
J02J
1
rI
dS ~ A4r~ 1. •
En portant toutes les majorations obtenues dans (60), on aura 1 F (NI) - F (No) 1 ~ A (A 2 r Afir~.l)'
e.! +
d'où l'on déduit en définitive la formule (59), où ~ > 0 est strictement inférieur à cx. 11-2-7. Dérivée du potentiel de simple couche suivant une direction quel. conque. Au (11-2-4] nous avons étudié les valeurs limites de la dérivée normale du potentiel de simple couche lorsque M -+ No suivant la normale. Si on impose à la densité fA. (N) des conditions plus fortes que la continuité, on peut prouver que les dérivées suivant une direction fixe convergent et de plus que ces limites ne dépendent pas de la façon dont M -+ No. On admet que la densité vérifie la condition de Lipschitz: 1 J.l. (N 2) - J.l. (NI) 1 ~ Brt2' (66) où rl. 2 = NIN 2' B et 6 sont des constantes strictement positives (6 ~ 1). Soit (NoX, Y, Z) un système de coordonnées locales d'origine No E S. Calculons la dérivée de u (M) suivant l'axe NoX situé dans le plan tangent à S en No. On admettra provisoirement que le point M se trouve sur la normale à S en No. On admettra pour fixer les idées que M est intérieur à S. On a 8u(M) 8x
) ) J.l.(N)
r~
dS
(r=MN).
(67)
S
Introduisons la quantité r' = l!~2
JJ
J.l. (No)
00
r~3
+ "1 +
cos (n, Z) dS = J.l. (N 0)
2
)
Z2
et considérons l'intégrale
J(;2+ ,/+
00
Z2)3/2
d~ d"l
(z =F 0), (68)
11-2-7.
où
(J~
D~RIV~E
291
DU POTENTIEL SUIVANT UNE DIRECTION
est le disque S2+T]2 ~
~
• On a de toute évidence d
J\ \J
~
(SJ+T]2+ Z2)8/2
a'o L'expression (67) devient
0 un nombre donné. Considérons la partie (JI de (Jo définie par S2 1'J2 ~ dl et choisissons dl assez petit pour que l'intégrale
+
) ) I~I 1fJ.(N)~fJ.(No)1 dS al
soit ~ :
en module quelle que soit la position du point M sur la normale
en No. Ceci est possible en vertu de la première majoration (76). Mettons ensuite VI. 1 (M) sous la forme V1, 1
(M)
=)) ~
fJ.
(N)~I-t (No)
al
dS
+
+ II
~ I-t(N)~I-t(No)
dS=V)l,>l (M)+v\:\(M),
0 donné, il existe un 11 tel que ôu(M) -
1
iJ n
(ÔU(N) ) {) n,
1--· ~E
pour
,1 '';;
II-2-fO. SUITES DE FONCTIONS HARMONIQUES
301
-le point M étant sur la normale intérieure à S en N. En fixant e on obtient 1aUa~M) 1~ (B e) pour MN ~ l), d'où 1 U (M 2) U (Ml) 1~ (B e) Ô1 ,2, où Ô1,2 = M l M 2 • De là il s'ensuit que u (M) -+ u (N) lorsque M -+ N suivant la normale en N. On peut écrire par ailleurs
+
+
ô
u (M) - u (N)
=)
au
~~1)
dô 1 ,
o où Ml est un point variable sur la normale en N, ôl = N Ml et Ô = N M, ôl ~ Ô ~ l). De la majoration de la dérivée normale, il vient: 1 u (M) - u (N) 1~ (B e) ô, d'où l'on voit que u (M) -+ -+ U (N) uniformément par rapport à la position de N sur S. En tenant compte de ceci, il est aisé de prouver (11-2-7] que u (M) -+ u (N) quelle que soit la façon dont M -+ N, et que u (M) est continue dans Di. Les raisonnements sont les mêmes pour De. En définitive, si la fonction u (M) admet une dérivée normale régulière elle est continue dans Di (resp. De US). Donc, l'application des formules intégrales ci-dessus est subordonnée à la seule existence des dérivées normales régulières des fonctions u (M) et v (M). Tout ce qui a été dit à propos deD i s'étend au plan. Si l'on considère un domaine illimité dans le plan, on verra plus loin que les choses seront différentes.
+
11-2-10. Suites de fonctions harmoniques. Avant de passer à la résolution des problèmes aux limites pour l'équation de Laplace à l'aide des potentiels de simple et de double couche, on se propose d'établir quelques propriétés des fonctions harmoniques qui viendront compléter celles que l'on connaît déjà. Considérons des suites de fonctions harmoniques ou, ce qui revient au même, des séries de fonctions harmoniques. Toutes les démonstrations seront effectuées pour le plan. Dans l'espace, elles sont exactement les mêmes. Il suffit seulement de remplacer la formule de Poisson par la formule qui donne la solution du problème de Dirichlet pour la sphère. Le thé 0 r ème fondamental des séries uniformément convergentes de fonctions harmoniques présente une très grande ressemblance avec un théorème analogue de la théorie des fonctions d'une variable complexe régulières (tome 111 2 , (1-12]). Si les termes de la série 00
2J
k=f
Uh
(x, y)
(100)
sont des fonctions harmoniques à l'intérieur d'un domaine borné B et continues dans B et si cette série converge uniformément sur le contour
CH. II.
302
PROB~MES
AUX LIMITES
l de B, alors elle converge uniformément dans B et sa somme est une fonction harmonique dans B. Soit donné un nombre e > O. La convergence étant uniforme sur l, il existe un N tel que pour tout n ~ N et tout p > 0 l'on a n+p
1
~
(x, y) I~e «x, y) sont sur l).
Uk
k=n
Cette somme finie de fonctions harmoniques est une fonction harmonique à l'intérieur de B et continue dans B. En vertu de la propriété fondamentale des fonctions harmoniques relative aux extrémums sur le contour (tome II, [V11-3-3l), on peut affirmer que cette inégalité étant réalisée sur le contour l, elle le sera à fortiori en tout point intérieur, autrement dit elle est réalisée dans le domaine B, ce qui exprime la convergence uniforme de la série (100) dans B. Donc, la somme S (x, y) de la série (100) est une fonction continue dans B. Montrons qu'elle est harmonique dans B. Soient Mo un point quelconque de B et ~ 0 un disque de centre Mo et de rayon R tel que ~ 0 c B. Désignons par Sn (x, y) la somme partielle de la série (100). Cette somme est une fonction harmonique et les valeurs qu'elle prend en un point intérieur de ~ 0 s'expriment en fonction de celles qu'elle prend en un point frontière de ~ 0 à l'aide de la formule de Poisson:
où (p, 0 donné il existe un A > 0 tel que 1 u (M) 1~ 8 dès que r;;;::: A. Construisons une boule B 0 centrée en l'origine et de rayon assez grand pour que u (M) soit harmonique à t'extérieur et sur âB o. Construisons une fonction Vl (M') harmonique
CH. II. ·PROBLM4ES AUX LIMITES
312
à l'intérieur de B o et prenant les mêmes valeurs limites sùr ôB o que u (M). Soit 11,1 (M) l'image de la fonction VI (M') par la transformation de Kelvin par rapport à B o. La différence 11, (M) - 11,1 (M) est une fonction harmonique à l'extérieur de B o, nulle sur ôB o, qui tend vers 0 lorsqtier ~ 00. On a vu plus haut qu'une telle fonction est identiquement nulle. Donc, 11, (M) == 11,1 (M). On a vu par ailleurs que 1es .produits rUt
_2 ôu!(M) (M) , r ôx '
r
2 ôu!
(M) ôy
t
r
2
au! (M) ô:
(A 15) ~
sont bornés pour r ~ 00. On remarque donc que la convergence de 11, (M) vers 0 lorsque r ~ 00 entraîne que les produits (115) pour 11, (M) restent bornés lorsque r ~ 00. On dira qu'une fonction u (M), harmonique au point à l'infini, est régulière au point à l'infini si u (M) ~ 0 pour r ~ 00. Si l'on sait seulement que u (M) converge vers une limite finie b, alors on peut dire qu'une telle fonction est égale à la somme de la constante b et d'une fonction harmonique régulière au point à l'infini. Si pour des fonctions 11, (M) et v (,M) harmoniques à l'extérieur d'une sùrface S les produits (115) sont bornés et si ces fonctions possèdent sur S des dérivées normales régulières, alors comme on l'a vu au [II-2-91, elles sont justiciables des formules (94) et (95) dans lesquelles l'intégration est étend ue à la région de l'espace extérieure à S. Le pro b 1 ème e x t é rie u r d e D i r i chI e t consiste à trouver une fonction u (M) harmonique à l'extérieur de S, régulière au point à l'infini, continue à l'extérieur de S et sur S et prenant sur S des valeurs f (N) données à priori. En plaçant l'origine des coordonnées en un point Mo intérieur à S, on ramène par une transformation de Kelvin le problème extérieur de Dirichlet à un problème intérieur pour le domaine image. On démontre l'unicité de la solution du problème extérieur par les raisonnements habituels. L'existence de la solution du problème extérieur se ramène à celle de la solution du problème intérieur, laquelle peut être établie sous des hypothèses générales sur la surface avec des conditions aux limites continues. Signalons la différence qui existe dans les positions du problème extérieur de Dirichlet pour le plan et pour l'espace. Dans le cas plan nous avons donné les conditions aux limites sur la frontière et avons simplement exigé que la fonction tende vers une limite finie lorsque r ~ 00. Dans l'espace, nous imposons cette limite, à savoir nous la supposons égale à o. Il est évident qu'on aurait pu admettre que la fonction u (M) tend vers une limite b pour r ~ 00 et revenir à la position précédente en considérant la différence u (M) - b. Il est immédiatde voir que dans l'espace il ne suffit pas d'exiger que u (M)tende vers une limite finie lorsque r ~ 00 pour que le pro-
II-2-14. UNIClT:e DE LA SOLUTION DU PROBUlME DE NEUMANN
313
blème extériéur soit bien posé. En effet, supposons qu'une certainequantité d'électricité se trouve en équilibre sur une surface conductri-ce S. Ce potentiel électrostatique de simple couche sera égal à une~ constante c sur S et de plus on démontre sans peine que U (Ml est une fonction harmonique à l'extérieur de S et tend vers pourr ~ 00. La constante c sera elle aussi une fonction harmonique à l'extérieur de S et prendra sur S les mêmes valeurs limites, mais elle ne sera pas régulière au point à l'infini au sens de notre définition. Dans le plan, ce raisonnement n'a pas de sens, car le potentiel électrostatique de simple couche sur une courbe l est infini à l'infini. Signalons encore que la fonction u (M) est dite harmonique à l'extérieur de S dans le cas seulement où elle est régulière à l'infini, c'està-dire que certains auteurs incluent dans la définition d'une fonction harmonique à l'extérieur de S la régularité au point à l'infini. Le pro b 1 ème e x t é rie u r deN e u man n consiste à trouver une fonction harmonique à l'extérieur de S, régulièreà l'infini dont la dérivée normale prend des valeurs données sur S. Les valeurs frontières de la dérivée normale ne doivent plus satis-faire ici à la condition (110). La démonstration de cette conditiol1 dans le cas plan ne s'étend pas à l'espace, car l'aire de la sphère de· rayon R est de l'ordre de R2 et l'intégrale de :: sur une sphère derayon assez grand ne doit pas tendre vers lorsque R ~ 00. Si l'on suppose que la solution du problème extérieur de Neumann possède· une dérivée normale régulière, alors l'unicité de la solution résulteimmédiatement de la formule (94). Nous avons effectué le même· raisonnement pour le problème intérieur de Neumann. Signalons en conclusion que les propriétés de la fonction harmonique au voisinage du point à l'infini peuvent être directement acquises à partir du développement de cette fonction suivant desfonctions sphériques au voisinage du point à l'infini.
°
°
11-2-14. Unicité de la solution du problème de Neumann. On sepropose de prouver l'unicité de la solution du problème intérieur de Neumann sans supposer que la dérivée normale est régulière. On considère au préalable un solide de forme spéciale sur lequel on. construit une fonction harmonique douée de propriétés qui seront précisées plus bas. Soit T (ex, k, h) un solide limité par la surface1+a
J=
+ y2)-2-
k (x 2
et le plan z = h, où k, ex et h sont des constantes strictement positi-ves. Soit No (0, 0, 0) le sommet de ce solide. Désignons par a' laI portion de frontière contenue dans le plan z = h et par a" la portion restante. Soit par ailleurs Uo et Ul deux constantes réelles telles queU o < Ul. Con.struisons une fonction w (M) harmonique à l'intérieur-
CH. II. PROBLœiES AUX LIMITES
'314
.(le T (ex, k, h), continue dans T (ex, k, h) et telle que w(M) < Ut sur 0', w (M) ~ U o sur 0" et w (No) = u o. Si r = V x 2 y2 Z2, 8 est l'angle du rayon vecteur de M et de l'axe NoZ, et P n (x), la fonction de Legendre (tome 111 2 , [VI-1-131), alors on sait (tome 111 2 , tVI-1-7l) que pour tout n la fonction rnP n (cos 8) sera harmonique à l'intérieur de T (ex, k, h). Construisons w (M) de la forme w (M) = Y [r cos 8 r 1 +PPl+P (cos 8)] Uo, (117)
+ +
+
+
0 et ~ > 0 sont des constantes qui seront déterminées plus loin. De toute évidence w (No) = U o. Montrons tout d'abord que pour tou t ~ assez proche de 0 P1+rdO) < O. (118) La fonction P n (x) est la somme d'une série hypergéométrique: 1 Pn(x)=F n+1; -n; 1; 1-x 2 ) (Ixl 0 et les k - 2 autres, y,
où y > 0 est un nombre quelconque. On admet qu'il existe un solide T (a, k, h) tangent à S en No et dont tous les points sauf No sont intérieurs à S. Prenons h fixe assez petit pour que l'on ait (119) sur la portion (J" de la surface du solide. Soit u1 la plus petite valeur de u (M) sur (J'. La fonction u (M) n'étant pas constante, on a U o < Ul' En choisissant y assez petit, on peut construire une fonction w (M) douée des propriétés mentionnées plus haut et telle que w (No) = u (No) et w (N) < u (N) sur le reste de la surface de T (a, k, h). En orientant l'axe NoZ dans le sens de la normale intérieure à S en No, on obtient u(M)-u(No) NoM -
ce que nous voulions.
K(M)-.(N o) z
>
w(M)-w(No) z
>
y,
(122)
316
CH. II. PROBL1!:MES AUX LIMITES
Ce théorème entraîne immédiatement l'unicité de la solution du problème de Neumann au sens suivant: Si une fonction u (M) harmonique à l'intérieur de Di est continue dans ÏJ i et (a~~N») i = 0 sur la surface S tout entière, alors u (M) est constante,. Soit No le point de S en lequel u (M) présente son minimum. De (122) il résulte aussitôt que la dérivée suivant la normale en No ne peut tendre vers 0 lorsque M ~ No suivant la normale en No. Si cela était vrai, la formule des accroissements finis nous donnerait u (M)-u (No) ~ 0 J
J
ce qui contredit (122). L'unicité de la solution du problème extérieur de Neumann se prouve exactement de la même manière. La démonstration exhibée a été donnée dans le travail commum de M. K e 1 d Y c h et M. L a v r e n t i e v (DAN SSSR, 1937,. 16, nO 3). La démonstration du théorème d'unicité est élémentaire dans le cas où la surface S est tangente intérieurement à une sphère (S. Z a-: rem b a, UMN, 1946, 1, nOS 3-4). 11-2-15. Résolution des problèmes aux limites dans l'espace. Posons les problèmes intérieurs de Dirichlet et de Neumann pour un domaiBe Di limité par une surface S. On cherchera la solution du problème intérieur de Dirichlet sous la forme d'un potentiel de double couche:
u (M) =
JJ
fA. (N) cos ~~. n) dS,
(123)
8
-+
où r= MN etn est la normale extérieure au point N de S. On cherche la densité fA. (N). En vertu de la première formule (42), le problème intérieur de Dirichlet avec la condition aux limites (124) (U)i 18 - t (N) est équivalent à l'équation intégrale suivante pour la densité fA. (N) : t(N o)=
JJfA.(N)
cos~o.
8
n)
dS+2nfA.(N o)
(ro = NoN).
0
En considérant le noyau K (N . N) o,
= __1_ 2n
cos (ro. n) r~
,
on peut mettre cette équation sous la forme: Jl (No) =
2~ t (No) +
JJIl (N) K (No; N) dS. 8
(125)
11-2-15. ReSOLUTION DES PROBL:I!:MES AUX LIMITES DANS L'ESPACE
31't
Le noyau K (No; N) n'est pas symétrique, car la normale est prise -~ en N et r o = NoN. On obtient le noyau transposé de l'équation adjointe (tome IVu [1-9]) en prenant la normale en No et r o = NN o• Le noyau transposé KI (No; N) est donc défini par la formule: ~
K 1 (N 0 ., N) -- K (N'' 0 N ) --
cos (ro, no) 2nr~
(126)
,
"0 est la normale extérieure en No.
Nous avons changé le signe du noyau, car le sens de r o a été inversé dans KI (No; N) et conservé dans la formule (126). Cherchons la solution du problème intérieur de Neumann avec la condition aux limites OÙ
(127) sous la forme d'un potentiel de simple couche: U
(M) =
JJ s
fi
~N) dB.
(128)
La première formule (49) nous conduit à une équation intégrale équivalente au problème posé:
f (No) =
) ) fi (N) cos (~r no) dB + 211;J.' (No).
s Cette équation peut être mise sous la forme
J.' (No) = 2~ f (No) - ) ) Jt (N) KI (No; N) dB. (129) s De façon analogue, les deuxièmes formules (42) et (49) nous donnent pour les problèmes extérieurs de Dirichlet et de Neumann avec les conditions aux limites (Ue>ls= f (N) ;
( auan(N) ) • = f (N) les équations intégrales
JI (No) = -
2~ f (No) +
) ) JAr (N)
cos
(;r
n)
dB,
s JAr
(N 0) =
-
_1_
2n
f (N 0) + J{" J{" JI (N) s
cos (ro, no) dB 2nr2 , 0
dont les solutions sont cherchées sous la forme (123) (pour le problème de Dirichlet) et (128) (pour le problème de Neumann).· Consi"
CH;;, U.
318
PROB~~ES Â, UX
LIMITES
dérons les équatjOD,S
~ (No) = q:> (No) + Â ,
f1 (No) =
O. Donc, la nécessité de la condition (141) résulte de la continuité de la solution cherchée dans D. Nous avons indiqué plus haut les propriétés des dérivées des solutions du problème de Neumann à l'approche de la surface S. L'étude analogue de la solution du problème de Dirichlet est plus compliquée du fait que cette solution se représente par un potentiel de double couche. Le potentiel de double couche a été étudié par A. Lia pou nov dans le travail mentionné plus haut ainsi que dans Sur le principe fondamental de Neumann dans le problème de Dirichlet (1902). Enumérons les résultats obtenus par A. Liapounov sur le potentiel de double couche et la résolution du problème de Dirichlet. 1. La valeur prise sur S par le potentiel de double couche w (N) de densité continue est une fonction vérifiant la condition de Lipschitz dans un rapport quelconque < 1. 2. Si le potentiel de double couche de densité continue possède une dérivée normale régulière d'un côté de S, alors il possède une dérivée normale régulière de l'autre côté aussi et ces dérivées normales sont égales en tout point de S. 3. Pour que la solution du problème intérieur ou extérieur de Dirichlet qui prend des valeurs f (N) continues sur S, possède une dérivée normale régulière, il est nécessaire et suffisant que le potentiel de double couche de densité f (N) possède une dérivée normale régulière sur S. Ce théorème a été prouvé dans l'hypothèse que le nombre ct qui figure dans la condition (3) soit égal à 1. 4. Soient F (x, y, z) une fonction continue avec ses dérivées premières et secondes, définie dans un voisinage de S, et 1 (N), la valeur prise par F (x, y, z)
II-2-18.
PROBL~MES
AUX LIMITES DANS LE PLAN
325
sur S. Sous ces conditions la dérivée, suivant une direction quelcmique fixe, d'une fonction u (M) harmonique dans Di ou De et prenant les valeurs f (N) sur S est continue dans DioU De US. Ce théorème a été établi' sans admettre que ct = 1. A. Liapounov a prouvé aussi une condition suffisante pour que le potentiel de double couche possède une dérivée normale régulière. Exhibons-la. Soit No un point de S qui est pris pour origine des coordonnées polaires (p, 8) dans le plan tangent à S en No. Les valeurs de la densité l.l (N) au voisinage. de No peuvent être traitées, par projection de N sur le plan tangent, comme une fonction de p et de 8. Posons 2:rt
Jr l.l(p,
1 ~(p)= 2n
8)d8.
o La condition indiquée ci-dessus se ramène à la suivante: il existe deux nombres > 0 et ~ > 0, tels que pour tout point No l'on ait
b
I~ (p)-~ (No) 1 ::;;; bpfH-1.
Ce théorème a été prouvé par Liapounov dans l'hypothèse que ct = 1. D'autres résultats relatifs à la théorie du potentiel newtonien de volume de simple et de double couche, ainsi qu'à l'étude du problème de Dirichlet sont accessibles dans l'ouvrage déjà cité de N. Gunter et dans l'article: Kh. S m 0 1 i t ski, Estimations des dérivées des fonctions fondamentales. DAN SSSR, 1950, 24, nO 2.
11-2-18. Problèmes aux limites dans le plan. Les problèmes de Dirichlet et de Neumann dans le plan se traitent exactement comme dans [11-2-15]. La solution du problème de Dirichlet est cherchée sous la forme d'un potentiel de double couche u (M) = ) ~ (N) cos (~' n) ds,
(144)
l
et celle du problème de Neumann, sous la forme d'un potentiel de simpIe couche: u (M) =
i~ l
(N) ln
+
ds.
(145)
On obtient les équations ,intégrales suivantes pour la densité:
i
~(No)=cp(No)+Â ~(N)K(No; N)ds,
(146)
l
~ (No) =
cp (No)
+ Â ~ ~ (iV) K
l
(No; N) ds,
l
où
(r o= NoN).
(147)
326
CH. II. PROBLeMES AUX LIMITES
L'équation (146) correspond au problème intérieur de Dirichlet f (No) et au problème extérieur de Diripour  = 1 et qJ (No) = ~ :Tt chlet pour  = - 1 et cp (No) = _...!f (No). L'équation (147) défi:Tt nit le problème extérieur de Neumann pour  = 1 et qJ (No) = = _..!. f (No) et le problème intérieur de Neumann pour Â=-1 :Tt . et qJ (No) =..!..f (No). Dans tous ces cas, la fonction :Tt dans la condition aux limites. L'équation (146) peut être mise sous la forme
qJ
(No) figure
lo
P- (so) =
qJ
(so)
+ Â l p- (s) K (so, s) ds,
- -
o
où set So sont les longueurs des arcs LN et LN0 du contour l, mesurées à partir d'un point fixe L dans un sens donné, lo la longueur du contour l. L'équation (147) s'écrit de façon analogue. Les noyaux K (so; s) et KI (so; s) sont continus sous les hypothèses imposées au contour l dans (I1-2-8l. Comme au [11-2-16], le nombre  = 1 n'est pas une valeur propre et  = - 1 est une valeur propre simple. Ceci étant, pour l'équation (146), la fonction propre associée à cette valeur propre est une constante et pour l'équation (147), c'est la densité électrostatique P-o (N) pour laquelle le potentiel de simple couche (145) est constant sur et à l'intérieur de l. Les problèmes intérieurs plans de Dirichlet et de Neumann se résolvent comme dans l'espace. Penchons-nous sur les problèmes extérieurs. La solution du problème extérieur de Neumann est reliée à l'équation (147) pour  = 1. Ceci étant, la fonction cp (No) doit satisfaire à la condition [II-2-12] : l
lo
cp (No) ds o = O.
En intégrant les deux parties de (147) pour  == 1 par rapport à No, on obtient z l1(N o)ds=O, o
J
et par suite la fonction (145), où 11 (N) est la solution de l'équation (147) pour  = 1, est une fonction harmonique régulière à l'infini [11-2-12], donc est solution du problème extérieur de Neumann.
11-2-19.
~QUATION INT~GRALE
POUR FONCTIONS SPH:eRIQUES
Passons au problème extérieur de Dirichlet. Si la condition ) flo (N)
f
327
(N) satisfait
f (N) ds = 0,
l
.alors la substitution de la solution de l'équation (146) pour À = - 1 dans la formule (144) nous donne la solution du problème. Si cette condition n'est pas remplie, on prend alors une constante .a telle que [11-2-16] ) f-lo (N) [f (N)
-al ds=
°
l
et comme plus haut on obtient, grâce à la formule (144), la solution (M) qui prend les valeurs f (N) - a sur l et la somme w (M) a sera la solution cherchée qui prend les valeurs f (N) sur l. L'adjonction de la constante est due au fait que la formule (144) nous donne une fonction harmonique nulle à l'infini, alors que dans le plan la solution du problème extérieur ne s'annule pas à l'infini.
+
:Il)
11-2-19. Equation intégrale pour fonctions sphériques. Considérons l'équation homogène (133) dans le cas d'une sphère unitaire ~ centrée --+ en l'origine. Dans ce cas, no est de même sens que le rayon ON 0 et cos (r o, no) = - r~ de sorte que l'équation homogène (133) sera Ge la forme (148)
On serait arrivé au même résultat en partant de l'équation (132). L'intégrale du second membre représente la valeur au point No (8 0 , CPo) du potentiel de couche sphérique de densité - 4~fl (N) =
= -
Â
4n fl (8, cp) .
Considérons tout d'abord ce potentiel en un point M (p, sr, cp') intérieur à ~. En désignant par r la distance MN et par p, la distance OM, on aura le développement (tome 111 2 , [VI-1-4])
+= ~ 00
Ph (cos 1') ph (p< 1),
(149)
h=O
Pk (x) sont les polynômes de Legendre et l' l'angle des rayons vecteurs OM et ON. Prenons pour fl (8, cp) une fonction sphérique d'ordre
OÙ
~
n:
~
CH. II. PROBUlMES AUX LIMITES
Le développement (149) qui converge uniformément pour p nous permet d'écrire
ll Y
(8, 'fP) +da=
n
2:~t
Y n (8', . (k= 1,2, ••. , n).
m=O i*k
En retranchant la somme p-1
2J
m=O
u m • k (N)
des deux membres de l'égalité précédente, on obtient p-1
n
2J LJ
m=O i=l
um,i(N)=!h(N)-up,h(N)
sur 8 h (k=1, 2, ••• , n). (163)
II-2--23.
M~THODE
DE 'SCHWARZ "(SUITE)
337
Si l'on démontre que pour p -+00 les fonctions u p , k (M) (k = 1, 2, ... , n) tendent uniformément vers 0 dans le dO,maine D" alors de (163) il s'ensuivra que la fonction p-1
n
LJLJ
m=O i=1
U m.
dM}
qui est harmonique à l'intérieur de D et continue -dan's D; nous donnè la solution du problème de Dirichlet pour D, qui prend les valeurs tk (N) sur Sk: 00
U
(NI)
=
LJ
n
~ u m • i (M).
(164)
m=O '1.=1
Etudions maintenant les conditions sous lesquelles les fonctions Up,k (NI) tendent uniformément vers 0 dans le domaine D. Désignons pàr Vk (M) (k = 1, 2, ... , n) une fonction harmonique à l'intérieur de D k et égale à 1 sur Sk. Ceci étant, Vk (M) ~ 0 à l'intérieur de DA et puisque Vk (M) -+ 0 lorsque M -+ 00 il ré'suIte qu'il existe une constante qk: 0 < qk < 1, telle que . Vk (M) ~ qk sur Si pour i =1= k (k = 1, 2, ... , n). (165) Si Wk (M) (k = 1, 2, ... , n) sont des fonctions harmoniques à l'intérieur de D k , continues dans D k et vérifiant la condition 1 Wk (M) 1~ ak sur S k (k = 1, 2, ... , n), (166) où ak sont des constantes, alors akvk (M) - wk(M) seront des fonctions harmoniques à l'intérieur de D k et positives sur Sk' d'où il s'ensuit que akvk (M) -!:!!k (M) ~ 0 dans ï5k~ c'est-à-dire: qlîe Wk (M) ~ ahvh (M) dans D h. Sans modifier la condition (166), peut changer le signe de la fonction harmonique Wk (M) et considérer de la sorte que Wk (M) ~ 0 au point M envisagé. Des raisonnements précédents il s'ensuit donc , 1 Wk (M) 1~ akvk (M), (167)
on
d'où, en vertu de (165), 1 Wh (N) 1~ akqk sur Si pour i =1= k (k - 1,2, .•. , n). (168) Cette inégalité découle donc de (166). Soit a > 0 un nombre tel que 1tk (N) 1~ a pour k = 1, 2, ... , n et q ',max {ql' Q2' ... , qn}. Il est évident que 0 < q < 1. En vertu de (160), il vient 1· UO,k (N) 1~ ~ a sur S h' En vertu de (161) et (168), on a 1 Ul.k (N) 1~ (n - 1) aq sur Sk (k = 1,2, ... , n). En se servant de (162) pour m = 2 et de (168), on obtient 1 U2. k (N) 1 ~ (n - 1)2aq2 sur Sk et, d'une façon générale 1Up.k (N) I~. (n - 1)P aqP sur S k et par suite IUp.k(M) I~ f(n 22-010t7
1) q]P~ (M E D k) (k- 1, 2, ... , n). (169)
338
CH. II.
PROB~MES
AUX LIMITES
Si le nombre de surfaces n = 2, alors il s'ensuit de là que Up,k (M) -+ -+ 0 lorsque p -+ 00 uniformément dans DA et à fortiori uniformément dans D. Si n > 2, on obtient la condition suffisante suivante pour que u p , h (M) -+ 0 : (170) (n - 1) q < 1. Par construction, le nombre q ne dépend pas des conditions aux limites fh (N) et est défini uniquement par le domaine D. On aurait pu exactement de la même manière considérer le cas où le domaine D est un domaine fini borné par une frontière extérieure SI et des frontières intérieures S 2' S 3, • • • , Sn' Ceci étant, nous aurions eu affaire au problème intérieur de Dirichlet pour le domaine Dl limité par 8 1 , et au problème extérieur pour les domaines DA (k = 2, ... , n). La construction indiquée ci-dessus a été proposée par G. G 0 1 0 u z i n e (Matem. Zbornik, 1934, 41, nO 2) (en russe). On s'assure sans peine qu'elle n'est pas valable dans le plan en se donnant des valeurs frontières constantes sur les contours fermés lh' 11-2-24. Fonctions sub- et superharmoniques. La méthode des équations intégrales de résolution du problème de Dirichlet assujettit la frontière du domaine à des conditions relativement fortes. On se propose d'exposer une autre méthode de résolution de ce problème, valable pour des conditions assez générales sur la frontière du domaine et les valeurs frontières. Cette méthode, appelée souvent méthode de balayage, a été proposée par Poincaré ensuite précisée par Perron (O. P e r r 0 n. Eine neue Behandlung der ersten Randwertaufgabe jür du = O. Math. Z., 1923, 18, voir également l'article de 1. P é t r 0 v ski, Méthode de Perron de résolution du problème de Dirichlet. UMN, 1941, 8 et son ouvrage sur les équations aux dérivées partielles) et par Vallée et Poussin. Dans ce numéro, on définit des notions nouvelles qui seront utilisées pour développer la méthode indiquée. Ces notions présentent un intérêt en physique mathématique. On se place dans le plan. Les raisonnements sont exactement les mêmes dans l'espace. La différence dans l'étude des valeurs frontières de la fonction harmonique construite d'après cette méthode est mentionnée à la fin. L'équation de Laplace pour la fonction d'une seule variable indépendante y (x) est y" (x) = O. Son intégrale générale est le binôme y = ax b et son graphique, une droite. Le problème de Dirichlet, c'est-à-dire le problème qui consiste à déterminer la solution de l'équation y" (x) = 0 à l'intérieur de l'intervalle ra, b] prenant des valeurs données aux bornes de cet intervalle, consiste tout simplement à mener une droite par les deux points donnés. Le binôme du premier degré est caractérisé par le fait que sa valeur en un point quelconque x = Xo est égale à la moyenne arithmétique de ses va-
+
II-2-24. FONCTIONS SUB ET SUPERHARMONIQUES
339
+
leurs aux points x = Xo h et x = Xo - h. Considérons maintenant une courbe continue ayant sa concavité tournée vers les y > o. Soit y = y (x) son équation. En tout point de cette courbe, on a 1
Y (x o) 0 tel que f (P) soit inférieure à la valeur moyenne de f (M) sur un cercle de centre P et de rayon p < ô. Si x et y sont les coordonnées de P, cette condition s'écrit alors f(x,
y)~ 2~
2:rt
l f(x+pcos0 (à l'intérieur de B).
domain~
(173 1)
En appliquant la formule de Green au disque K p de centre P (x, y) .contenu à l'intérieur de B et en posant u=f et v=1, on obtient
~ :~ cp
ds=
~~
/1fda,
(174)
Kp
où Cp est la frontière du disque K p • Appliquons encore à la fonction .f la formule (tome II, [VII-3-2])
2~
f (x, y) =
J(f
ô~: r
-ln r
:~
) ds +
cp
2~ ~
J/1f ln r da,
Kp
où r est la distance du point de coordonnées (x, y) au point variable d'intégration. Sur Cp, la direction de n est celle de r, et ds = pdçp. En se servant de (174), on peut mettre la formule précédente sous la forme 2:1:
tex, y)=
2~
)
f(x+pcosçp, y+psinçp)dçp+
2~)) /1fln
;'da.
o Kp Dans le disque K p , on a r/p ~ 1 et, en vertu de (173 1), la dernière formule nous donne l'inégalité (172 1), autrement dit, si la condition (173 1 ) est remplie, la fonction f (M) est une fonction subharmonique à l'intérieur de B. De façon analogue, si t!f ~ 0 (dans B), (173~) alors f (M) est une fonction superharmonique à l'intérieur de B. La définition principale des fonctions sub- et superharmoniques n'implique pas l'existence des dérivées. Les conditions (173 1 ) et (173 2 ) sont identiques aux conditions classiques de convexité et concavité d'une courbe (tome l, II-5-2). Etudions quelques propriétés simples des fonctions sub- et superharmoniques. Supposons que f (M) est une fonction continue dans un domaine fermé et subharmonique à l'intérieur de ce domaine. Soùs cette condition il s'ensuit aussitôt de (172 1) qu'une fonction subharmonique prend sa plus grande valeur sur le contour. De plus, à l'intérieur du domaine elle ne peut avoir un maximum au voisinage duquel elle n'est pas constante. D'une façon analogue, une fonction superharmonique prend sa plus petite valeur sur le contour.
11-2-25. Propositions auxiliaires. Prouvons quelques propositions Sur les fonctions sub- et superharmoniques qui nous seront utiles lors de la résolution du problème de Dirichlet. -
341
II-2-25. PROPOSITIONS. AUXILIAIRES
r ème I. Si ht (M) (k = 1, 2, ... , m) sont des fonctions continues dans un domaine B et subharmoniqucs dans B, alors la foncThé
0
tion
n (M) = 'li> (M) pour tout n. Les suites {'lI>n (M)} peuvent varier d'un point Mo à un autre. Prouvons le Thé 0 r ème. La fonction u (M) est harmonique dans B. Dans la suite, les fonctions seront affectées d'indices supérieurs notés entre parenthèses. Prouvons préalablement le Lem m e. Si P est un point arbitraire fixe de B, il existe alors une suite monotone de fonctions supérieures (J)(1) (M) ~ (J)(2) (M) ~ ... (M
E B)
(178)
telles que (J)(n) (P) ~ u (P). On a vu plus haut qu'il existe une suite {'lI>n (M)} de fonctions supérieures 'lI>n (M) telles que 'lI>n (P) ~ u (P). Posons (J)(n)
(M) = min {'lI>1 (M), 'lI>2 (M), .•. , 'lI>n (M)}.
(179)
Les fonctions (J)(n) (M) sont, nous l'avons vu, des fonctions supérieures continues. Lorsque n croît, le nombre de fonctions 'li> s (M) dont on considère le minimum croît aussi et par suite çp(n) (M) satisfont la condition (178). Le fait que la fonction u (M) est la borne inférieure des fonctions supérieures et la relation (179) entraînent que u (P) ~
,CH. II.PROBUlMES AUX
~
LIMITE~
cp(n) (P) ~ 'Pn (P) d'où cp(n) (P) ~ u (P) Ceci achève la démonstration du lemme.
pui~que
'l'n (P)
~
u (P).
Rem a r que. Soit K un disque de c-entre P contenu dans B. Construisons les fonctions cpW) (M) comme indiqué dans le théorème III'. Les inégalités (178) étant vérifiées sur la frontière du disque K, elles le seront sur le disque K tout entier. A l'extérieur de K, les fonctions cp'lP (M) sont confondues avec cp(n) (M) et par suite vérifient aussi la condition (178): cp 0 tel que dans ~ÔI' l'on ait 'Pl (M) ~ w (No) 28 (M E ~Ô2)' et à fortiori u (M) ~ w (No) 28 (M E B et ~Ô2)' (185) Soit ô = min{ô 1 , ô 2 }. D'après (184) et (185), on a w (No) - 28 ~ u (M) ~ w (No) 2e (M E B et ~ô). (186)
+
+
+
Le nombre 8 étant arbitrairement choisi, il s'ensuit que u (M) ~ ~ w (No) lorsque M ~ No de l'intérieur du domaine, ce qui prouve le théorème. Cette démonstration est valable dans le plan comme dans l'espace. Si w (N) est continue sur l et si la condition 1 est remplie pour chaque point No de l, alors la fonction u (M) est continue dans B et est égale à w (N) sur l. Dé fin i t i o:n. On dit qu'un point frontière No est régulier si pour toute fonction w (N) continue sur l, la fonction u (M) ~ w (No) lorsque M ~ No. Les points frontières ne jouissant pas de cette propriété sont dits points frontières irréguliers. Le théorème prouvé plus haut nous dit que la condition 1 est une condition suffisante de régularité du point No. Indiquons maintenant une condition suffisante simple de nature géométrique pour qu'un point frontière soit régulier dans l'espace. Supposons qu'un point frontière No jouit de la propriété suivante: il existe une boule qui ne contient aucun point de B autre que No. Soient Ml le centre de cette boule et R son rayon. Posons r = M 1M et considérons la fonction 1
1
w(M)~--R r •
CH. II. PROBU:MES AUX LIMITES
348
Cette fonction satisfait visiblement la condition 1 et de plus est harmonique dans B. Plaçons-nous maintenant dans le plan et supposons que la frontière de B est composée d'un nombre fini de courbes fermées simples d'équations x = x (t), y = y (t), où x(t) et y (t) sont des fonction~ de t périodiques continues (fig. 14). Supposons tout d'abord que le point No se trouve sur le contour extérieur lI. Plaçons l'origine des coordonnées en z = 0 et choisissons l'échelle de telle sorte que le domaine B contienne le disque 1 z 1 < 1. Composons la fonction 1
o
F(z)=--
Fig. 14
ln z·
Lorsque z se déplace dans B il ne peut contourner l'origine et F (z) est une fonction continue dans B et régulière dans B et de plus F (0) = O. En posant z = peï\P, on obtient pour la partie réelle de F (z) l'expression ln p
W
(z) = - (ln p)2+cp2 ,
ln p < O. Cette fonction harmonique remplit toutes les conditions ci-dessus. En particulier, à l'extérieur de .~ 8 on a w (z)
>-
InR (ln e)2...1-",2' l '1'0
où CPo est la plus grande valeur de cp dans Ë et R, la distance de l'origine au point de jj le plus éloigné de cette origine. , Supposons maintenant que No est situé sur le contour l2. Choisis~ sons à l'intérieur de l2 un point quelconque ex. et effectuons la transforma tion conforme
z
,
1 =--
z-a.·
Le contour l2 se transforme en un contour extérieur et l'on peut construire une fonction w (M) pour le point considéré par la méthode ci-dessus. En revenant à la variable z, on obtient la fonction cherchée. Donc, si 00 (N) est continue sur le contour l, alors u (M) le sera à l'intérieur et sur l et sera égale à 00 (N) sur l. Supposons maintenant que No est un point de discontinuité de première espèce de 00 (N), autrement dit 00 (N) tend vers des limites différentes selon que N -+ No le long du contour d'un côté ou de l'autre. Désignons ces limites par 001 (No) et 00 2 (No) et supposons que 00 1 (No) < 00 2 (No). En raisonnant comme plus haut, au lieu de
I1-2-27.
~TUDE
DES VALEURS.
FRONTI~RES
349
i
_~
(187) Si (ù (N) est une fonction bornée, c'est-à-dire vérifie les conditions (177), alors on a vu qu'il en est de même de la fonction u (M). Donc u (M) est une fonction harmonique bornée prenant les valeurs frontières (ù (N) en tous les points de continuité de (ù (N). Revenons maintenant à l'espace. On peut construire Une surface relativement simple possédant des points irréguliers. Cette construction a été proposée par Lebesgue et ensuite, indépendamment de lui, par P.S. Uryshon. Le problème des valeurs frontières est étudié plus en détail dans un article de M. K el d Y c h (UMN, 1941, 8). Citons encore un exemple dans le cas du plan. Soit B un disque centré en l'origine des coordonnées et privé de son centre. L'ensemble l des points frontières est constitué du cercle de ce disque et de son centre. Supposons que (ù (N) = 0 sur le cercle et (ù (N) = 1 au centre. Une telle fonction (ù (N) est continue sur l. La fonction harmonique u (M) tend visiblement vers zéro lorsque M tend vers le cercle. Montrons que u (M) ne peut tendre vers 1 lorsque M tend vers le centre. S'il en était ainsi, u (M) serait harmonique dans le disque au cas où elle serait égale à 1 au centre [I1-2-121. Or ceci contredit le théorème de la valeur moyenne d'une fonction harmonique au centre du disque. Donc, l'origine des coordonnées est un point frontière irrégulier. 0 dans le cas considéré. En On démontre aisément que u (M) effet, u (M) est bornée et par suite converge lorsque M tend vers le centre (I1-2-12] et si cette limite est égale à la valeur de u (M) au centre, alors u (M) est partout harmonique dans le disque [II-2-12] et nulle sur le cercle, c'est-à-dire que u (M) O. A noter qu'au lieu de la condition l, on aurait pu poser une condition portant uniquement sur le voisinage du point No et montrer en outre que cette nouvelle condition est identique à la condition 1.
==
==
Con dit ion II. Pour un voisinage ~" du point No, il existe une fonction w" (M) continue sur ~'" superharmonique dans ~1l et telle que w" (No) = 0 et w" (M) > 0 dans ~,,' {No}.
350
CH. II. PROBU;MES AUX LIMITES
On démontre que dans l'espace le point No vérifie la condition II s'il est le sommet d'un cône circulaire dont les points assez proches de No sont situés à l'extérieur de ff (sauf No). Donc, ces points sont réguliers (cf. 1. P é t r 0 v ski, Leçons sur le.séquationsauxdérivée$ partielles, M., Fizmatguiz, 1961 (en russe)). Dans la suite, sauf mention expresse du contraire, on admettra que les contours ou les surfaces limitant les domaines étudiés ont tous leurs points réguliers. Ce sera par exemple le cas des surfaces de Liapounov. Pour ces dernières nous avons construit la solution du problème de Dirichlet à l'aide de la théorie du potentiel et des équations intégrales. Si une fonction continue (ù (N) est définie sur la frontière et si tous les points de cette dernière sont réguliers, alors la fonction harmonique construite u (M) est continue à l'intérieur du domaine et sur la frontière et prend sur cette dernière les valeurs (ù (N). On sait que la fonction u (M) est unique. Si la frontière présente des points irréguliers, alors la fonction harmonique u (M) est bornée à l'intérieur du domaine et prend les valeurs (ù (N) en tous les points frontières réguliers. On démontre qu'il n'existe qu'une seule fonction possédant de telles propriétés. On peut trouver la démonstration de cette proposition et d'autres faits intéressants, dont le critère de régularité de Wiener, dans l'article déjà cité de M. Keldych. 11-2-28. Equation de Laplace dans un espace à n dimensions. Jusqu'ici nous avons étudié l'équation de Laplace dans le plan et dans l'espace. Les résultats acquis se généralisent facilement à un espace à n dimensions dans lequel l'équation prend la forme n
.1u=2J ux.x.=O. i=1
l
t
Citons les principaux résultats concernant les solutions de cette équation. Les fonctions qui possèdent des dérivées premières et secondes continues et qui sont solutions de cette équation sont dites harmoniques. La solution singulière fondamentale est de la forme C rn- 2
(n> 2),
la constante C étant égale à (n- ~) W ,où O>n est l'aire de la n sphère unité dans l'espace à n d'imensions, de sorte que la solution singulière fondamentale s'écrit 1
2).
II-2-29. FONCTION DE' GREEN DE
L'OP~RATEUR
DE LAPLACE
35~
Le volume Un d'une boule de rayon r à n dimensions s'exprim& par la formule (tome II, [III-4-11]) n Un
(2n)2 n 2' r n (n- ) ..• 2
=
pour n paIr,
n+ 1 n-l
Un
_ - -22 2 n n •. n(n-2) ... 1 r pour n ImpaIr,
=
ce qUI, on le vérifie sans peine, peut être mis sous la forme unIque
d'où, en dérivant par rapport à r et en posant r= 1, 2
w
n
(v~n)n
= r (~ ) .
Pour une fonction harmonique dans un domaine D limité par une' surface S on a la formule (tome II, [VII-3-3]) u (M) =
1 ( s
(j)o
(r)
:~
- u
a~on(r)
)
dS
(on n'écrira partout qu'un seul signe d'intégration). La quantité rest la distance du point variable d'intégration à M. Les propriétés. fondamentales des fonctions harmoniques restent en vigueur, y compris le théorème de la moyenne d'une fonction harmonique au centre d'une boule et l'unicité de la solution du problème de Dirichlet. La formule qui donne la solution du problème de Dirichlet pourune boule de rayon R est de la forme u(M)=_1_r (N)
Jf
CJJnR S
(R2_ p2)dS n ,
(188),
(R2+ p2_2Rp cos 8)2
où p est la distance du centre 0 à M, N un point variable de la boule,. e l'angle de ON et OMo La méthode des fonctions supérieures et inférieures 'se généralisesans changement à un espace à n dimensions. La condition de régularité des points frontières établie plus haut est valable aussi en, dimension n. 11-2-29. Fonction de Green de l'opérateur de Laplace. Nous pou-vons définir la fonction de Green pour une équation aux dérivées.
_:
352
CH. II.
PROB~MES
AUX LIMI'l'ES"
partielles de la même manière que nous l'avons fait pour une équation différentielle ordinaire. Commençons par définir la fonction de Green pour l'équation de Laplace dans l'espace avec l'une des conditions aux limites homogènes suivantes:
u/s=o, au
. ~+p (N) uls=O
(189)
.
(190)
(p (N) >0).
Nous pouvons construire la fonction de Green aussi bien pour un domaîne fini Di intérieur à S que pour un domaine illimité De extérieur à S. Commençons par le domaine Di. La fonction de Green G (P; Q) doit être une fonction du couple de points (P ; Q). En tant que fonction de P, la fonction G (P ; Q) doit, dans Di" {Q}, posséder des dérivées premières et secondes continues et être harmonique, et sur la frontière, satisfaire la condition aux limites. En tant que fonction de P, elle doit par ailleurs présenter une singularité au point Q correspondant à une charge (ou une masse) ponctuelle en Q. En tenant compte du facteur 4n figurant dans la formule (tome II, [VII-2-81) il [) ) )
~ ~M)
dTMJ =
- 4rqt (Mo) Mo ED, (r = MoM) , (191)
Di
on définit la fonction de Green pour les conditions (189) ou (190) de la manière suivante: D é fin i t ion. On appelle fonction de Green de l'opérateur de Laplace, correspondant aux conditions aux limites (189) ou (190), une fonctionG (P; Q) vérifiant comme fonction de P, à Q E Di fixe, les conditions suivantes: 1) cette fonction est harmonique dans Di" {Q} ; 2) elle vérifie la condition aux limites (189) ou (190); 3) elle peut être mise sous la forme 1 +g(P,• Q), G(P; Q)=G(x, y, z; ~,11, ~)= 4nr
où r
=
(192)
PQ et g (P ; Q) est une fonction partout harmonique dans D i'
Construire la fonction de Green revient à en trouver la partie régulière g (P ; Q). Les conditions aux limites (189)~et (190) deviennent respectivement pour la fonction g (P; Q): 1
g (N: Q) = - - 4 ,!'l ES, (r= NQ); nr
C1g (~n; Q)L + p (N) g (N ; Q) =
ai.
-
4~ [ a~
+ p ~N)
(193)
] ,
NES. (194)
"~
II-2-29. FONCTION DE GREEN DE
L'OP~RATEUR
DE LAPLACE
353
Donc, la construction de la fonction de Green se ramène à la résolution des problèmes de Dirichlet et mixte pour l'équation de Laplace et l'on peut admettre comme acquise l'existence de la fonction de Green si S est une surface de Liapounov. . Dans le cas du domaine De' il faut inclure dans la définition de la fonction de Green la condition de régularité à l'infini, c'est-à-dire que G (P ; Q) à Q fixe situé à une distance finie, doit tendre vers 0 lorsque P ~ 00. Soit Di un domaine borné de frontière S. Le problème de Dirichlet admet dans Di une solution distributionnelle vérifiant la condition aux limites (193). La formule (192) définit la distribution de Green pour le domaine Di avec la condition aux limites (189). Si No est un point frontière régulier, alors G (P ; Q) ~ 0 pour P ~ ~ No. On peut prouver la réciproque, savoir que si G (P; Q) ~ 0 pour P~No' alors No est un point frontière régulier. Dans le plan, la fonction de Green se définit de façon analogue, sauf qu'il faut remplacer la formule (192) par G (P; Q)
1 1 = -2 ln -r + g (P ; Q) . n
(195)
Les formules (192) et (195) nous disent que G (P; Q) = 00 pour P Q et G (P ; Q) > 0 pour les P assez proches de Q. Le point Q s'appelle pôle de la- fonction de Green. On étudiera la fonction de Green uniquement pour la condition aux limites (189). Montrons que G (P ; Q) est une fonction continue des points P et Q dans Di, pourvu gue ces points ne soient pas confondus. En vertu ~e (192), on peut affirmer que la démonstration de la continuité de G (P; Q) peut être ramenée à celle de la continuité de g (P; Q). Estimons la différence g (P' ; Qf) _ g (P"; Q") ; en ajoutant et en retranchant g (P' ; QIf), on obtient
==
1g
(P'; Q') -
< 1 g (P';
g (P"; Q")
1~
Qf) _ g (P'; Q")
1
+
1
g (P' ; Q") -
g (P"; Q")
1.
La différence g (P' ; Q") - g (P"; Q") est la différence des valeurs de g (P ; Q") aux points P' et P" et cette différence tend visiblement vers 0 lorsque P" ~ P'. La différence g (P' ; Qf) - g (P' ; Q") est la valeur en P' de la fonction harme)ilique g (P; Q') - g (P; Q") 1 qui prend les valeurs frontières -4 (~-~) sur S,. où r' et r" n r r sont les· distances du point variable NES aux points Q' et Q". Si Q" est assez proche de Q', alors la: valeur absolue de la différence est aussi petite que l'on veut lorsque N varie sur S. Or la fonction harmon~gue g (P; Q') - g (P; Q") prend sa valeur minimale et sa valeur maximale sur la frontière S et l'on
(+- ;,,)
23-01017
354
CH. II. PROBL:E:MES AUX LIMITES
peut affirmer que g (P' ; Q') - g (P' ; Q") ~ 0 pour Q" ~ Q'. Ceci prouve la continuité de la fonction g (P; Q), donc de G (P ; Q). La fonction G (P; Q) est strictement positive au voisinage du point Q et nulle sur S, donc elle est strictement positive dans le domaine Di. Ce raisonnement est valable dans l'espace pour le domaine De. Etablissons une inégalité simple pour G (P; Q). La fonction g (P ; Q) prend les valeurs strictement négatives (193) sur S. Donc, g (P; Q) < 0 dans Di et par suite 1
0< G (P; Q) < 4nr
dans
D,
(r= PQ).
(196)
On obtient une estimation analogue pour De. Plaçons-nous maintenant dans le plan. Soit d le diamètre d'un domaine fini B du plan, c'est-à-dire la borne supérieure des distances des couples de points de B. La fonction harmonique g (P ; Q) + 2~ ln ~ prEnd sur l lES valeurs 2~ ln ~ qui sont strictement négatives quelle que soit la pcsition du pôle dans B. 1 1. 1 1 Donc, g (P; Q) +:Et ln ëf< 0, 1. e. g (P; Q) < - 2n ln d dans B. Ce qui nous donne 1
1
1
1
G (P ; Q) < "7) ""n ln -r - -2 n ln -d '
c'est-à-dire 0< G (P; Q) < a ln
~
+b
(dans B),
(197)
où a et b sont des constantes. Les doubles inégalités (196) et (197) nous donnent des encadrements de la fonction de Green qui dépendent de la distance r des points P et Q. 11-2-30. Propriétés de la fonction de Green. Considérons la fonction de Green dans Dien désignant comme plus haut par r la distance du point courant au point Q E Di. Définissons la fonction g (P; Q),
v (P) = {
1 -
4nr '
PEDit PED e-
(198)
Cette fonction est harmonique dans Di et De et nulle à l'infini. Dans De elle possède des dérivées de tout ordre continues dans De U S. On peut traiter v (P) comme la solution dans De du problème de Neumann qui prend les valeurs frontières
f (M)= -
1 418
:n (~),
MES,
(199)
II-2-S0.
PROPRI~T~S
DE LA FONCTION DE GREEN
355
et représenter de la sorte v (P) comme un potentiel de simple couche de densité continue: v(P)= ) ) ~;:V) dS s
(r'=NP).
(200)
Ce potentiel est égal à - 4....!sur S, où r = MQ, autrement dit, il nr prend les mêmes valeurs que g (M; Q). De là il résulte que la formule (200) est valable dans l'espace tout entier, c'est-à-dire que g (P ; Q) =
) ) ~ ~~V) dS
(P EDt),
(201)
s
et par suite g (P ; Q) possède dans Dt des dérivées normales régulières sur S. On peut de toute évidence en dire autant de la fonction G (P; Q). Signalons à propos de la condition aux limites (199) que la fonction
1 T
possède, quel que soit Q E Di' des dérivées de tout ordre non seulement sur S mais dans un voisinage de S. Sur S le second membre de (199) vérifie la condition de Lipschitz 1
f (N 2 )1 - f (NI)
1 ~ arl.2
(rl,2
= N I N 2 ),
et l'on peut affirmer que Il (N) vérifie la condition de Lipschitz [11-2-5] et par suite G (P; Q) possède des dérivées premières continues dans Dt [11-2-7].
Prouvons maintenant que la fonction de Green est symétrique, c'est-à-dire que G (P; Q) = G (Q; P). (202) A noter qu'on vient de démontrer que G (P; Q) possède des dérivées normales régulières sur S. Par ailleurs elle possède des dérivées continues dans Dt' {Q}. Appliquons maintenant la formule ) ) ) (u L\v-v L\u) dT=
Di
JJ(u :: -v :: ) dS B'
aux fonctions u = G (P; QI) et v = G (P; Q2)' le domaine d'intégration Di étant constitué de Dt privé de deux sphères SI et S2 de centres QI et Q2 et de rayon e petit. L'application de cette formule est licite d'après ce qui a été dit plus haut. L'intégrale triple étendue à Di est nulle, puisque la fonction de Green est harmonique en dehors des pôles. L'intégrale étendue à S est nulle en vertu de la condition aux limites (peu importe laquelle), donc ] dS=O. 23·
356
Cff. II. PROBL:ElMES AUX LIMITES
Au point Q2' la fonction G (P; QI) ne présente aucune singularité, quant à la fonction G (P; Q2)' elle tend vers l'infini comme..!.. Le r
°
produit de .!par l'aire 4ne 2 de la sphère tendant vers avec e, B on voit que les seuls termes de la formule à ne pas tendre vers avec e sont ceux qui contiennent la dérivée normale de G (P; Qi), i = 1, ~, au voisinage du point où G = +00. Ces termes seront au nombre de deux, et, en les explicitant, on obtient la somme
4~ ~ ~ G (P ; QI) ~J
a_1_
'a:
dS -
4~ ~ ~ G (P; Q2)
°
a _1_
a:'·
dS
+ 11 = 0,
SI
où 11 -+ 0 avec e, r~ = PQ1' r 2 = PQ2' La formule de Green fait intervenir la 'normale extérieure, donc dans les dernières formules la normale doit être orientée vers l'intérieur des sphères, c'est-à-dire dans le sens contraire aux rayons et l'on a '--'-'
' ; j
4~B2
)l Œ
(P; QI) dS -
S2
4~82 ~ ~ G (P;
Q2) dS + 11 = O.
SI
En appliquant le-théorè:r;ne de la moyenne à ·ces intégrales, on peut écrire , G (P 2; QI) - G (Pl; Q2) 11 = 0,
+
est un point de S2 et pour e -+ 0, on trouve OÙP2
P1~
un point de SI' En passant à la limite
G (Q2; QI) = G (QI; Q2)'
ce qui prouve' que la fonction de Green est symétrique. Il s'ensuit de cette égalité que g (P; Q) est continue en (P, Q) dans Di X Di à l'exception de l'ensemble où; P = Q E S. Dans le cas d'une boule, la fonction de Green s'écrit (tome II, [VII-3-7l)
1 (1
R ) G(P' . , Q)=4n r -pr!'
(203)
où p est la distance de Q au' centre, rI' la distance de P à Q', symétrique de Q par rapport à la sphère, R, le rayon de la boule. En désignant par (x, y, z) et (6, 11, ~) les coordonnées de P et Q, on peut écrire
11-2-31. FONCTION DE GREEN DANS LE PLAN
357
En dérivant (203) par rapport à x par exemple, et comme 1x -
~:ç
1
-'--...:....---:.. ~ 1, rI
on trouve aG(p;
ax
l
Q)I~_1 ~ 4n
(_1 +~) r pri· 2
Les points P et Q étant intérieurs à la sphère, rI> r et G (P; Q) > on obtient
> 0,
aG(p; Q) 1~_1_
1
ax
~
2nr 2
(204)
•
On a des majorations analogues pour les autres dérivées partielles. Supposons que u (M) est la solution du problème intérieur de Dirichlet pour un domaine Di de frontière S avec les conditions aux limites f (N). Si l'on sait que u (M) possède une dérivée normale régulière, alors on peut appliquer la formule (91) en posant v = = g (P; Q). On obtient alors (tome II, [VII-3-7l) u (Q)
= - ) ) f (N)
aG
(:n;
Q) dS N.
(205)
S
A. Lia pou nov a prouvé que cette formule donne la solution du problème de Dirichlet pour toute fonction continue f (N) de la condition aux limites. On lui doit également la première démonstration rigoureuse de la symétrie de la fonction de Green. Ces résultats et d'autres portant sur la théorie du potentiel abordée plus haut se trouvent dans son travail Sur quelques questions liées au problème de Dirichlet, 1898 (en russe). 11-2-31. Fonction de Green dans le plan. L'étude de la fonction
de Green dans le plan diffère sur quelques points de celle effectuée dans l'espace. On envisagera la fonction de Green pour un domaine borné Bide contour l avec la condition aux limites (189) sur l. Définissons comme dans (11-2-30] la fonction v (P) dans le plan:
v (P) =
{
g (P; Q)
à l'intérieur de l,
2n
à l'extérieur de l.
__1_ ln 1r
(206)
Construisons comme au [II-2-30] le votentiel de=simple couche= VI
(P) =
~ ~ (S) ln 1
+-
ds,
(207)
CH.. Ii. PROB~MES AUX LIMITES
358
où r' = PN, N étant un point variable de l. La dérivée normale ( ÔV~~P) ) e prend sur l les valeurs
f (M) =
--
2~ :n ln ~ =
-
cos2~; n), MEL,
(208)
-+ où r est de même sens que MQ, et n, la normale extérieure à l en M. Composons maintenant la fonction
w (P) = ) JI (s) ln ~ ds + 2~ ln ~
(r = PQ, Q est intérieur à l)
l
(209)
harmonique dans Be et telle que ô~ M.
11-2-31. FONCTION DE GREEN DANS LE PLAN
359
Donc, la fonCtion (211) est harmonique dans Be' possède sur lune dérivée normale régulière nulle et tend vers 0 lorsque P -+ 00. Cette fonction est justiciable de la formule
d'où il s'ensuit que w (P) = 0 dans Be, c'est-à-dire que
r
1
J II (s) ln --;:' ds = -
1
2n
ln
r1
dans Be.
1
De là il résulte aussitôt comme dans (I1-2-301 que le potentiel de simple couche (207) est confondu avec la fonction v (P) définie par (206) sur le plan tout entier et l'on peut affirmer que g (P ; Q) possède sur l une dérivée normale régulière. Par ailleurs, comme dans [11-2-30] on peut affirmer que g (P; Q) possède dans Bides dérivées premières continues dans Bi. On démontre que G (P; Q) est symétrique exactement comme dans [II-2-301. Pour un disque de rayon R la fonction de Green s'écrit G (P ; Q) =
2~ ln ';; ,
(212)
les notations étant les mêmes que dans [II-2-30l. Ceci nous conduit aux majorations suivantes: ôG (P; ôx 1
0) 1:::;:=_1_. ~ nr '
ÔG (P; Q) 1:::;:=_1_ ôy ~ nr • I
(213)
La solution du problème de Dirichlet dans Bi est justiciable d'une formule identique à (205). La fonction de Green de l'opérateur de Laplace pour un domaine simplement connexe plan avec la condition aux limites (189) est étroitement liée avec une fonction réalisant une transformation conforme du domaine considéré en le disque 1 w 1 < 1 (tome II1 2 , [11-7]). Soient B un domaine simplement connexe de contour l et Zo = S l'li, un point intérieur à B. Supposons par ailleurs que w = f (z) est une fonction réalisant une transformation conforme de B sur le disque unité, et de plus f (zo) = 0, c'est-à-dire que l'image de Zo est le centre du disque. Si le contour l remplit certaines conditions de régularité, alors la fonction f (z) est continue sur le disque 1 z 1 ~ 1 et transforme le cercle 1 z 1 = 1 en l. La transformation étant injective, f (z) possède une raCIne simple en z = Zo:
+
f (z)
=
(z -
zo) [a o
+ al (z -
zo)
+ ... J
(a o =1= 0).
(214)
360
CH. II. PROBL:E:MES AUX LIMITES
Formons la fonction G (x, y ; ~, '1'1) =
-
1
2n ln 11(z) 1·
(215)
On vérifie immédiatement que cette fonction est la fonction de Green pour le domaine B avec un pôle en (~, '1'1). En effet, ln 1 1 (z) 1 est la partie réelle de ln 1 (z), donc est une fonction harmonique. En vertu de (214), la partie infinie de la fonction (215) sera égale à 1 ln 1 1 1 au point (~, '1'1) et enfin le contour l de B se trans2n Z-Zo forme en le cercle unité, autrement dit Il (z) 1 = 1 sur l, et la fonction (215) est nulle. Soit H (x, y; ~, '1'1) la fonction conjuguée harmonique de (215). On a G -t- il! =
- 2~
ln 1 (z) ,
(216)
et par suite l'on peut exprimer 1 (z) en fonction de G et de H:
f (z) = e- 2Jt(G+Hi). La fonction H est définie à un terme additif constant près, donc au second membre de la dernière formule on aura un facteur constant arbitraire de module 1 qui correspond à une rotation arbitraire du disque 1 w 1 < 1 autour de l'origine. Supposons que le contour l du domaine B est doué de la propriété suivante. L'angle 8 (s) formé par la tangente à l avec une direction quelconque fixe vérifie, comme fonction de la longueur d'arc s, la condition de LipEchitz (217)
où b > 0 et ~ > 0 wnt des constantes. On prouve que dans cette condition la dérivée f' (z) est continue dans B et il existe des constantes m > 0 et .M > 0 telles que m~
1
f'
(z)
1~
(218)
M.
Ces constantes dépendent évidemment du choix du point Zo qui se transforme en l'origine des coordonnées dans le plan de la variable w. Fixons le point Zo et construisons maintenant une transformation conforme de B sur le disque 1 w 1 < < 1, telle que l'antécédent de l'origine soit un point z' intérieur à B. Pour cela, il faut effectuer d'abord une transformation conforme w = f (z) et ensuite appliquer le disque 1 w 1 ~ 1 sur lui-même de telle sorte que l'image du point f (z') soit l'origine des coordonnées. Cette application est une application homographique et en écartant le facteur constant de module 1, on obtient en définitive -2nG(z· z')=Re [ln t
f(z)-f(z') 1-f(z)f(z')
J,
où Re désigne comme toujours la partie réelle et G (z, z'), la fonction de Green pour le domaine B de pôle en z'. En dérivant par rapport à x, où-x est une direc-
361
11-2-32. EXEMPLES
tion quelconque, on trouve -2n iJG (z; z') = Re [ iJx
+ W'J f' (z)
l'
(z) f(z)-f(z')
J
=
1-f(z)f(z')
2
_ Re { f' (z) [1-lf (z') 1 ] [f (z)- f (z')] [1- f (z) f (z')]
}
,
ou, si l'on remplace la partie réelle par le module 2n/ iJG~; z') I~
x
2
II'(z)II1-If(z')1 1
,
If(z)-f(z')II1-f(z)f(z')1
et si l'on tient compte du fait que If (z)1 < 1 et If (z')1 < 1, on trouve iJG (z; z') 1 If' (z)1 11-lf(z')1 2 1 211' (z)1 2n 1 iJx < If (z)-f(z') 1 (1-If(z')I) < If(z)-f(z')I· (218 1 ) Soit z = cp (w) la fonction réciproque de w = f (z). Cette fonction est définie sur le disque 1 w 1 ~ 1. De (218), il s'ensuit que 1 cp' (w) 1 ~ 1/m et l'on obtient w'
1 0, alors en tenant compte de l'encadrement (196) on peut prendre dl assez petit pour que l'intégrale de la formule (222) étendue à Di soit < ~ en valeur absolue pour tout point P intérieur à la sphère Sd1 (No). Lors de l'intégration sur Di, le point Q appartient à Di et le point P est supposé se trouver dans un petit voisinage du point No, par exemple, dans une sphère Sd 1/2 (No). Dans ces con.
II-2-33. FONCTION DE GREEN ET
~QUATION
AVEC SECOND MEMBRE
365
ditions, PQ > ~l et de [II-2-291 il s'ensuit que G (P; Q) est une fonction continue de (P; Q). Donc, dans l'intégrale étendue à Di on peut passer à la limite par rapport à P lorsque P -+- No et cette limite est nulle, car G (P ; Q) satisfait la condition (221). Par conséquent, l'intégrale étendue à Di sera < ~ en module si P est assez proche de No et par suite l'intégrale de la formule (222) sera pla
av. [au]=iroeiCt)(t-:)v [u]=eiCt)(t-:)v', [auJ=eiCt)(t-:) an an' _at ' et en intégrant sur Sp on obtient une intégrale de la forme
~~
e-;k
T
(~~ + ikv) as + ))
;2 e-
ikT
as,
(253)
sp
sp
p. Il est naturel d'exiger que cette expression tende vers 0 lorsque p --+ 00 (ceci correspond à une absence de source de vibrations à l'infini). L'aire élémentaire de la sphère contient le facteur p2 et la condition indiquée ci-dessus sera remplie si l'on impose à v les deux conditions suivantes:
()Ù r =
rv est borné et r (~~
+ ikv) ~ 0
pour r --+ 00, ces deux conditions devant être réalisées quelle que soit l'origine des rayons vecteurs r et uniformément par rapport aux directions de ces rayons vecteurs. Dans la suite on se servira des notations suivantes. Par 0 (,.ex) on désignera une quantité x telle que le rapport :x/ra soit borné lorsque r --+ 00 et par 0 (ra..), une quantité x telle qUé le rapportj x/ra.. --+ 0 lorsque r --+ 00, la convergence devant être uniforme par rapport à la direction du rayon vecteur r et ne pas dépendre du choix de son origine. Les conditions précédentes peuvent s'écrire: v = :J (r- 1) ; (254)
~~ + ikv =
0
(r- 1 ).
(255)
Ces conditions traduisent mathématiquement le principe de radiation dans l'espace. De façon analogue, elles s'écrivent dans le plan: 1
v=O(r-2); 1 av -a;+ ikv=o(r 2).
(256) (257)
376
CH..
n.
PROBL:Il:MES AUXLIMITES
La solution singulière fondamentale vérifiant le principe de radiation dans l'espace aura pour expression (P) =
V
e- ikT r
,
(258)
où r est la distance d'un point 0 fixe au point variable P. En dérivant la solution (258) par rapport à r, on s'assure qu'elle vérifie une condition plus forte que (255), plus exactement, au second membre au lieu de 0 (r- l ) on aura 0 (r- 2 ). Ceci étant, on admet que dans la formule (255), les distances sont mesurées à partir du même point O. Vérifions les formules (254) et (255) en admettant que les distances sont mesurées à partir d'un même point 0 1 et posons 0lP = p. Le fait que pv est borne résulte immédiatement de ce que p/r -+- 1. La formule (255) se vérifie par une simple dérivation de la solution (258) par rapport à p par l'intermédiaire de r. Ceci étant, on a àr
-=ccs'" àp J' où,\, -est l'angle des directions de r et de p; en appliquant la formule qui donne le carré du côté 00 1 dans le triangle OOlP, on obtient -cos,\,
1
=
+0
(r- 2 ).
(259)
Dans le plan, la solution fondamentale qui vérifie le principe de radiation sera la solution H~2) (kr) , où H~2) (z) est la deuxième fonction de Hankel. -Pour vérifier ceci, il suffit de se servir de l'expression asymptotique des fonctions de Hankel et de la formule
~ H~2) (z) =
-
H~2) (z).
(260)
La condition (257) sera remplie sous la forme forte, c'est-à-dire qu'au 3
1
second membre on aura 0 (r -2) au lieu de 0 (r- 2). En multipliant la solution H~2) (kr) par une constante qui ramène la singularité en r = 0 à ln .!, on obtient la solution r
v=
31
2, 1,
H~2) (kr).
(261)
On démontre comme plus haut que le principe de radiation sera aussi vérifié par les solutions H~)
(kr) cos mfl); Hi/;,) (kr) sin mfl) (m
= 1,
2, 3, ...).
(262)
Il-2-38. Théorème d'unicité. Si le principe de radiation a lieu, on peut démontrer le théorème d'unicité suivant: si une fonction v est solution de l'équation (252) à l'extérieur d'un contour fermé l, vérifie le principe de radiation à tinfini et une condition aux limites homogène.,
II-2-38. TH:eOR:eME D'UNICIT:e
par exemple la condition v quement nulle. _ Appliquons la formule
Il = ° ::/l = 0, ou
377
alors elle est identt,.:-
2 au! ) d s J\ J\ (Ul~U2-U2~UI) d't= Jr Jr (au u 1 an -u 2 a;-
(263)-
aBl
Bl
au domaine BI limité intérieurement par le contour l et extérieurement par un cercle ST de centre en un point fixe quelconque et de ray~~ assez grand, et posons U 1 = v et U 2 - li, où v est la conjuguée complexe de v. On admet que v est continue dans BI et possède une dérivée normale régulière. L'intégrale double est nulle en vertu de (252) et l'intégrale prise le long de l est nulle aussi en vertu de la condition aux limites. Reste l'intégrale étendue à ST' contour surlequel n est de même sens que r. La condition (257) nous permet deremplacer
-av
ar = -
1 -ikv .,.l- 0 (r 2). tl .
,
av-
--1
-
ar=ikvrt-o(r 2),
et l'on est conduit à l'égalité 2ik
JIv
1
l2
ds+ ) v·o(r -"2) ds+
Sr
Sr
Jv.o(r-2") ds=O. 1
br
Comme v Vr- et V Vr- sont bornés lorsque r ~ 00, les deux der-niers termes tendent vers et en considérant l'angle polaire cp sur lecercle ST' on obtient
°
2:rt
Vrv 12 dcp~ O.
(264)
) 1
o
Appliquons à présent la formule de Green à la solution v et à la première solution (262). L'intégrale double est nulle comme précédemment et il reste les intégrales prises le long de l et de ST' donc l'intégrale sur ST ne dépend pas de r. Les deux solutions considéréesvérifient le principe de radiation, et de plus les solutions (262) satisfont la condition (257) dans la forme forte au même titre que la solution H~2) (kr). En se servant comme plus haut de la condition (257), on trouve que l'intégrale prise le long de ST tend vers 0, et comme elle ne dépend pas de r, elle est tout simplement nulle, autre-ment dit
Hi/i) (kr)
av Jr 7ir
Sr
dH(2) (kr)
cos mcp dcp - -"':';'~::-r"""""""-';'"
r
J v cos Sr
mcp dcp = O.
CH. II.
PROBL~MES
AUX LIMITES
Si l'on pose lm (r) =
) vcosmcpd (kr), où dm est une constante -aussi. L'équation de fermeture (tome IV l , [1-31) et la formule (264) ~ous disent que pour m fixe et r -+ 00 Cm
Vr H~) (kr)
et
dm
V;: Hg) (kr)
-+
O.
'Ûr, de l'expression asymptotique de Hi:"J (kr) il s'ensuit que VrH~1 (kr) reste supérieur en module à un nombre strictement positif pour les grands r, donc Cm =dm = 0, c'est-à-dire que lm (r) = = gm (r) = 0, d'où il résulte, en vertu de l'équation de fermeture, que v = sur les cercles Sr' Si l est un cercle, alors en assimilant Sr à un cercle concentrique ·à l, on trouve que v à l'extérieur de l, c. q. f. d.
°
== °
°
Dans le cas d'un contour quelconque, les raisonnements précéau voisinage du point à l'infini. On ·dents montrent que v = montrera plus bas [11-2-39] que v (x, y) doit être, comme dans le cas de l'équation de Laplace, une fonction analytique et, en vertu ·du principe du prolongement analytique, la nullité de v au voisinage ·du point à l'infini entraîne sa nullité partout à l'extérieur de l. Le théorème d'unicité se démontre de façon tout à fait analogue . dans l'espace. 11-2-39. Principe de l'amplitude limite et principe de l'absorption limite. Comme au numéro précédent, on peut montrer que la condition de radiation ;permet de sélectionner une solution unique de l'équation dv
+ k~v = -F (P)
(k
>
0),
àéfinie dans l'espace tout entier. On admettra que F (P) est une fonction de point . continûment dérivable définie dans l'espace tout entier et nulle à l'extérieur ..d'un domaine D fini. La solution indiquée est alors donnée par la formule 1 v(P)= 4n
eJr 5Jr -,-, --F (Q) d'LQ ikr
D
(r=PQ).
II-2-39. PRINCIPE DE L'AMPLITUDE LIMITE
379
On est conduit- à cette solution en envisageant le problème non stationnaire des ()scillations forcées sous l'action d'une force périodique. Plus exactement, de la formule de Kirchhoff (tome II, [VII-3-11]) il s'ensuit immédiatement que r:(P)= Hm u(P, t) e- ikt , t-+oo
()ù u (P, t) est la solution de l'équation des ondes li.u -
Utt
=
-F (P) eikt ,
qui vérifie des conditions initiales nulles. Pour cette raison on dit de la solution v (P) qu'elle est 1'« amplitude limite d'une oscillation périodique établie pour les grands t sous l'action d'une force périodique ». Un tel principe de sélection des solutions de l'équation (2641 ) s'appelle principe de l'amplitude limite. Un autre principe de sélection des solutions de l'équation (264 1 ) appelé principe d'absorption limite (cf. V. 1 g n a t 0 vs k i, Ann. Phys., 1905, t8) consiste en ce qui suit: on introduit un paramètre complexe -ie (e > 0) «( absorption ») dans l'équation (2641 ):
+
li.ve (k 2 - ie) ve = -F (P), et l'on prend la solution ve (P) qui tend vers 0 à l'infini (cette solution est uni-
«ue) :
Ve
1 (P)= 4JT.
r
') ~
D
J
-ir(ae-ib e )
r
e
F(Q) dTQ,
()ù a e - ib e = V k 2 - ie (a e > 0, b e > 0), de plus b e -- 0 pour e -- +0. Lorsque e -- +0, la solution ve (P) admet une limite qui est confondue avec la solution v (P) définie par (264 2 ), Les trois principes envisagés ici et au [11-2-37] conduisent à la même solution (264 2), Il semble naturel de s'attendre à ce qu'ils soient applicables à des problèmes plus généraux, par exemple aux problèmes aux limites pour des équations elliptiques dans des domaines illimités. Cependant ces principes trouvent des champs d'application différents. Ainsi, le principe de radiation est favorable au cas où l'équation (2641 ) est considérée dans l'espace tout entier ()u dans un domaine E contenant le point à l'infini (11-2-38]. Si le domaine E est par exemple la bande 0 ~ z ~ 1, alors l'équation (264 1) n'admet aucune solution nulle pour z = et pour z = 1 et satisfaisant les conditions de radiation "Sous la forme (254) et (255) (F (Q) 0). Néanmoins si l'on modifie légèrement ces conditions, le problème admettra une solution unique (cf. A. S v e c h n i k 0 v, Sur le principe de radiation, DAN SSSR, 1950, 73, nO 5). En vertu des deux autres principes applicables ici sans aucun changement, .cet exemple montre que les « conditions de radiation » doivent dépendre de la forme du domaine E à l'infini. Des considérations physiques suggèrent que le principe d'absorption limite n'est pas toujours applicable dans la forme indiquée .ci-dessus si E se rétrécit assez vite à l'infini. Le problème de l'applicabilité des principes formulés ici reste ouvert. Indiquons à ce propos les travaux de F. Rellich (J ahresber. Deutsch. Math. Verein., 53, 57) qui étudient la forme des « conditions de radiation » pour l'équation (2641 ) dans des domaines illimités d'aspect divers, le travail de A. P 0 v z !Il e r, Sur le développement des fonctions suivant les fonctions propres de l'opérateur -li.u + cu. Matem. sb., 1953, 32, nO 1, qui démontre le principe d'absorption limite pour l'équation
°
li.u
*'
+ q (P) u + k
2
u = -F (P)
dans un espace illimité à trois dimensions et enfin l'article de Mme O. Lad y je n s ka ï a, Sur le principe d'amplitude limite. UMN, 1957, 12, nO 3. Cet article est consacré au principe d'amplitude limite pour l'équation ci-dessus.
380
CH. II. PROBUMES AUX LIMITES
Les travaux de Povzner et de Mme Ladyjenskaïa montrent que le principe d'absorption limite possède un plus vaste champ d'application que le principe d'amplitude limite, du moins tel qu'il est formulé plus haut. Plus exactement, les solutions des équations .1ve
+ q (P) ve + (k 2 -
iE) Ve
=
-F (P),
nulles à l'infini, admettent une limite pour E -+ +0 si seulement k 2 n'est pas une valeur propre de l'opérateur .1v q (P) v; la limite Hm u (P, t) e- i kt des t .... +oo solutions du problème non stationnaire correspondant est susceptible de ne pa~ exister si l'opérateur possède une valeur propre au moins. Signalons qu'un nomq (P) v si l'équation .1v bre c s'appelle valeur propre de l'opérateur .1v q (P) v cv = 0 admet une solution non triviale dont le carré du module est intégrable sur l'espace tout entier. Ces principes sont étudié~ ,pour des équations plus générales dans des travaux de D. E i d 0 u s, B. Van ber g et autres. .
+
+
+
+
+
II-2-40. Problèmes aux limites pour l'équation de Helmholtz. La solution (258) de l'équation (252) présente une singularité i. en r r = 0, ce qui nous permet d'élaborer pour cette équation une théorie du potentiel identique à celle du potentiel newtonien pour l'équation de Laplace. En désignant par r la distance du point variable NES au point P, on aura dans l'espace les analogues suivants des potentiels de simple et de double couche
J~ (N) e-;k ., s
v (P) = \"
T
dS,
aan (e- r Jr J ~r(N)
ikT
W
(P)
= -
(265) )
dS,
s
où n est la normale extérieure à S en N. En mettant les noyaux e- ikT 1 e- ikT _1 des potentiels v (P) et w (P) sous la forme - r = -r r on obtient des potentiels déjà étudiés dans lesquels le passage à la limite lorsque P tend vers la surface s'effectue d'après les formules de [11-2-2] et [11-2-4]. Dans l'intégrale restante, le noyau ne présentera plus de singularité pour r = et le passage à la limite est possible sous le signe d'intégration. On obtient ainsi les formules analogues à celles de [11-2-2] et [11-2-4]:
+
,.
°
L= 2a~ L
=
+ .\ J~t (N) iJ~o ( e-:o
kTO
(No)
)
dS,
s
(ro=NoN)
-2n~(No) + JJ ~ (N) a~o ( e-r~ro) dS s
(266)
II-2-40.
PROBL~MES
;LIMITES~VRL'ÉQU~TIONDE
AUX
et Wi
W
= 2n~ (No) -
(No)
) )
s
~
HELMHOLTZ
ro e-:: ) as, :n (e-::ro ) dS,
381
(N) :n (
e(No) = -2n~ (No) - ) ) ~l (N) s
(267)
à noter que le noyau de l'intégrale est dans (266) la valeur prise en No par la dérivée suivant la normale no, et dans (267), la valeur prise en N suivant la normale n; l'intégration est étendue à S comme dans toutes les formules de cette nature. Dans le plan, on obtient les potentiels de simple et de double couche v (P)
= )~
(Pl)
~ H~2) (kr) ds,
l
w(P)
= )~
(N) af)n
[2~ H~2)
(kr)
(268)
Jds,
l
qui sont justiciables de formules analogues aux formules (266) et (267) dans lesquelles il faut remplacer 2n par n. Ces potentiels sont solutions de l'équation (252) et en vertu du choix spécial des noyaux, chaque élément des intégrales écrites et ces intégrales elles-mêmes vérifient le principe de radiation. Considérons le noyau K (N
0'
N
;
k)
f) = an
(e- ikrO
)
= -
2nro
e
-ikrO('k o +1) l r ~nr5
cos m
't'O'
~
OÙ
(8 = 1, 2, ... , n). Par ailleurs
°
n
Tr (AB) = Tr (A' [fJ.l' fJ.2' •.. , fJ.n]) = Lj {A'}ss
fl.,
s=1
d'où il s'ensuit immédiatement que si tous les fls ~ 0, alors Tr (AB)~ ~ 0 et si tous les fls;:::: 0, alors Tr (AB) ;:::: O. Ce qui prouve le lemme. Signalons encore un fait qui nous sera utile pour la suite. Soit u (P) = u (Xl' X 2 , • • • , Xn) une fonction continue réelle définie dans un ouvert D de l'espace de point générique (Xl' X 2 , • • • , xn) admettant dans D des dérivées premières et secondes continues. Supposons que la fonction u (P) atteint son maximum en un point Po E D. Sous ces conditions. la forme symétrique réelle n
. .2J uXi~k (P 0) ~iÇk
t, R=1
ne peut prendre de valeurs strictement positives (tome 111 1 , [11-2-4]), c'est-à-dire que les valeurs caractéristiques de la matrice symétrique réelle Il u xixk (Po) Il ne sont pas toutes strictement positives. De façon analogue, elles ne sont pas toutes strictement négatives au point où u (P) présente un minimum. Considérons l' équation elHptique linéaire n
L (u)
n
= L} aUPx,xk z
i, k=1
+ i=1 .L: biuxz. + cu = f,
(289)
en admettant que aik' bh c et f sont des fonctions continues dans un domaine fini D et que la forme quadratique n
~ aik~i6k
i, k=1
est définie strictement positive dans D. Signalons que la somme (290)
II-2-43.
UNICIT~
DE LA SOLUTION DU PROBUJME DE DIRICHLET
391
est la trace du produit des matrices symétriques réelles Il aik Il et " U xixk II· On admet que la fonction U (P) possède des dérivées premières et secondes dans D. L'unicité de la solution du problème de Dirichlet pour l'équation (289) repose sur le théorème suivant. Thé 0 r ème. Si c sans second membre
0 dans D, on peut choisir ~ de telle sorte que bl~ - au~2 < 0 dans D et ensuite choisir a assez grand pour que a - e-/3xl > 0 dans D. Ceci étant, dans l'équation (295) on a c' < 0 dans D et la fonction w (P) qui est nulle sur S est justiciable du théorème ci-dessus, d'où w (P) 0 dans D. De (294) il s'ensuit qu'il en est de même de v (P). L'hypothèse c ~ 0 est essentielle pour l'unicité de la solution du problème de Dirichlet. Il est aisé d'exhiber un exemple dans lequel l'équation L (v) = 0 avec c > 0 possède une solution nulle sur S mais non identiquement nulle. Un tel exemple nous est fourni par l' équation (297)
==
que nous allons considérer sur le carré 0 ~ Xl ~ n, 0 ~ X2 ~ n. Si k est un entier naturel, alors l'équation (297) admet la solution v = sin kXl sin kx 2 ,
qUI est nulle sur la frontière du carré. On rappelle que l'équation 8. v
+ Âv =
0
(298)
possède une infinité de valeurs propres strictement positives  = = 11,1; 11, = 11,2' ••• , telles que pour 11, = Âh l'équation possède des solutions non identiquement nulles mais nulles sur la frontière· S du domaine considéré. Rem a r que. Les résultats ci-dessus nous permettent d'établir quelques estimations pour les solutions du problème de Dirichlet. Signalons les plus simples d'entre elles. Soit u (P) la solution du problème L (u) = f dans D; u 1s = o. (299) On admet que c < 0 dans D. Désignons par !-t le minimum de et par AI le maximum de 1 f 1 dans D.
1c 1
11-2-44. L'ÉQUATION
t\.1! -
Âv
=
393;
0
Introduisons la fonction v en posant v = u + k, où k est une·constante. On a L (v) = L (u) + ck = f + ck, de sorte que v (P) est la solution du problème L (v) = f + ck dans D; vis = k. (300}< Supposons d'abord que k
=
M•• ~
-
Sous cette condition f
+ ck ~ 0 dans-
D et les solutions de l'équation L (v) = f + ck ne peuvent présenter de minimum strictement négatif dans D. Or v (P) = M > 0 sur S, donc v (P) ~ 0, c'est-à-dire v = u
+
~
~ 0 ou u ~ _
M ~
M • De ~
façon analogue, en posant k = - M , on obtient u ~ ~ • En: !1 !1 définitive on peut affirmer que la solution du problème (299) (si elle· existe' oit vérifier dans D l'inégalité M lul~-. !1
Considérons maintenant le problème
L(u)=O
uls=cp
D;
dans
(301);"
et désignons par N le maximum de 1cp 1sur S. En posant de nouveau v = u k, on est conduit au problème
+
L (v)
=
ck
dans
D;
vis = cp
+ k.
Posons k = N. Vu que l'on a convenu que c < 0 dans D, la fonction v ne peut présenter de minimum strictement négatif dans D. Pour les valeurs frontières de v (P), on a de toute évidence cp N ~ O. D'où, comme plus haut, il résulte que 1l = U N ~ 0, c'est-à-dire-' u ~ -N. De façon analogue, en posant k = -N, on trouve u ~ N, c'est-à-dire 1 u 1 ~ N dans D pour le problème '(301).
+
11-2-44. L'équation
~v
-
+
J.v = O. Considérons l'équation
~v -
Âv =
0,
(303)
0 est un nombre donné et posons le problème intérieur deDirichlet avec la condition aux limites v
Is
=
f
(N).
(304);
Les solutions de l'équation (303) ne peuvent présenter ni de maximums strictement positifs ni de minimums strictement négatifsdans Di [II-2-43], d'où il résulte que le problème posé admet une solution unique. Si la fonction f (N) est telle que -a ~ f (N) ~ b, où a > 0 et, b > 0 sont des nombres arbitraires, alors la solution du problème.
CH. II. PROBUlMES AUX LIMITES
294
0 est donné. Par ailleurs, en vertu de (326), on a pour les assez grands
1
JJJVÂ 'P (P, Â) d't- ~~ I~ ~ , D'
0 et
VÂ s (Â) =
H,
(330)
D-2-U. EXPRESSION ASYMPTOTIQUE DES VALJ!lURS PROPRES
401
alors (331)
la sommation étant étendue aux valeurs de k telles que Âk ~Â. Appliquons ce théorème à la série (328). Dans ce cas Ch = Âj;l et H = :11: et l'on obtient Hm _~
"'-++00
~
y Â. '" k-...-; ::::::Â
;k - 2:
2
ou, ce qui revient au même, (332) où
8 (Â) -+
0 pour
 -+
+ 00.
Si
 = Ân .,
alors (333)
Posons
et désignons le second membre de (332) par-cp (Â). C'est une fonction croissante de  et de plus 0 pour  (Â,,) =
" C"'. LJ k=l
(343)
Dans la formule (342) la sommation est étendue aux valeurs de k pour lesquelles Â", ~ Â. La fonction cp 0..,) est une fonction de Ât eroissante et positive: pour A < Â1 pour Âm~Â< Âm+t.
q> (A) = {O ,
Comme A",
+ A~ 2A
am
(344)
pour Â", ~ Â, il vient
cp(A)~2Â ~ Âkc~Â )."'~~
êt en vertu de (330), on obtient
,
cp (Â) = 0 (V~),
(345)
autrement dit, cp (Â}/Vl est borné pour  -+ 00. On a d'autre part 11-1
~ Â"'C~). = ~
k=l
~t
O'k (
k=l
Âk~Â - Âk+~ +1.) +1.,,0,+-).,
de plus
La formule (334) nous pèrmet d'écrire "
).n
"" c'" '~ LI A",+A = 1&=1 ',:.,1,
cp (x)
(X+A)I
d + 0'" x A,,+l·
'.
Or 'Un ='(P'(l n ) et, de (~45) il s'ensuit que 0'"/(1,, + 1) -.. 0 pour n -+ 00., donc la dernière formule nous donne 00
• (1) =,~
k=1
00
A1&c~A = ~ (x~~)2 "0
dx.
(346)
405
11-2-.6. DSMONSTRATION DU THSOlUl:ME AUXILIAIRE
De (345) il résulte aussitôt que la fonction à intégrer est de l'ordre de x-s/2 lorsque x -+ 00. Pour alléger l'écriture, si '1' (1) = a1b e (Â) Âb , où e (Â) -+ 0 pour 1 -+ 00, on écrira 'l' (Â) ~ aÂb (tome III" [V-151). Prouvons deux lemmes:
+
+
Lem m e I. Si 1 (Â) est définie pour tous les 1 > 0 assez grands~ admet une dérivée continue, ÂI' (Â) croît avec  et 1 (Â) ~ a1q (q > 0), alors l' (1) ~ aq1Q- l . Prouvons d'abord ce lemme pour a = 1 et q = 1. On al (Â) ~  et il nous faut montrer que l' (Â) ~ 1, c'est-à-dire il faut prouver que l' (Â) -+ 1 lorsque 1 -+ 00. Raisonnons par l'absurde. Si l' (Â) ne tend pas vers 1, il existeune suite de valeurs Ân , telle que Ân -+ 00 et l' (Ân ) -+ h, où h =1= 1. Supposons par exemple que h > 1. Soit y > 0 un nombre arbitraire.. Comme ÂI' (1) est une fonction croissante, il vient / (')..,n
+ y n ) -
/ (Â n ) _
y n
= y~"
-
ln+Vln
J
f'
(Â) d1>
Ân~~~n)
Ân+V"'n
l d:
~n
=
j'
~n)
ln (1
+ y) ..
ln
+
Le second membre tend vers !!:..ln (1 y) > 1 si y est assez voisin. y de O. Or de 1 (1) ~ 1, il s'ensuit immédiatement que / ( n +Y n ) - f ( n ) -+ 1. y /1
Cette contradiction prouve le lemme pour a = q = 1. Passons au cas 1 Q général. Posons 11. = 1 Q et fI (11) = ..!. a 1 (11 / ). On a
Donc, ,..,/; (11) est une fonction croissante et l'on peut appliquer le lemme à Il (,...) pour a = q = 1, d'où il vient
Il (,...)
~
1,
Le.
1 l..-1 ..!.. q q ) ~ 1, -l1 f' (l1 aq
et f' (1) ~ aqÂQ-1, ce qui prouve le lemme. Cette démonstration reste en vigueur pour h = 00. Considérons l'intégrale 00
Kp =
J
u
1 2
P+-
----,-,..."........,..2=-
o
(U+1)2p+
du
(p
= 1, 2, ... ).
(347)
. 'CH. II. PROBLt;:MES AUX LIMITES
406
La substitution u = x/(1-x) nous permet d'écrire cette intégralesous la forme (tome 1112 , [111-17)): 1
Kp
=
. 1 P 21 x + (1-x)P- 2dx=
J
(P+.!.)
r (p+.!) r 2 r
(2P+2)
(348}
2
o Lemme II. Soient l-a.
Jo
K p ,l =
u
P+-12
(U+1)2p+2
du; 1 P+u 2
r
00
J
K p ,3=
ô~,
du,
1+a.
-où 0 < a < 1. Alors
.qù
(U+1)2p+2
Kptl~Ô~Kp; K pt2 >(1-ô;)K p ; Kp.3~6;Kp-,
(349)
ô; et ô; qui dépendent de a, tendent vers 0 lorsque
p~ 00.
En appliquant la formule de Stirling (tome III" [111-20]) 1
r (z) = ·V21t zZ-2e- z [1 + 8 (z)]
(8 (z) ~ 0 pour Z~ 00),
au second membre de (348), on obtient 3
K p = V21t 2- 22-21'
(p+~)P+l 2
(p+1)
(p+.!.)P
3 2p+-
2
(1 +e p ),
2
()ù 8 p ~ 0 pour p ~ 00. Le produit de la fraction par vers 1 lorsque p ~ 00, donc
'V p
tend
1
Kp
= Ap- 22-2P (1 + 8~)
(8; ~ 0),
(350)
3
ù A= V2i2- 2. La fonction u/(u+1)2 présente un maximum égal il 1/4 en u=1, d'où
T' (U~1)' 1
K P.1";;'k"
r(U~1)' 1
du< k
P
du,
-où 0 < k < 1/4 et k dépend de a. On obtient donc K
Pt 1 ~AlkP,
(351)
()ù Al = 1t/2, et de façon analogue
K P. 3~AlkP.
(352)
407.
11-2-46. DSMQNSTRATION DU. TH:eOR~ME AUXILIAIRE
On a kP = (1/4- ô)P, où Ô> 0 et dépend de cx. De ce qui précède, il s'ensuit que k P.2 2P p 1/2=(1-4ô)P p1/2-+0 pour p-+oo et en vertu de (350), (351) et (352), on obtien~ les inégalités .(349) pour K P. 1 et K p. 3" On a par ailleurs K P.
2
= K P - K P. 1 - K p. 3> K p - (ô; + ô;) K P'
de là on déduit (349) pour K p. 2' ce qui prouve le lemme. Passons à la démonstration du théo;fème formulé au [11-2-45]. 1) 'après l'hypothèse de ce théorème
r cp (x) J (X+Â.)2
1
00
s (Â) =
dx ~ HÂ
- -
(353)
2.
o Considérons la fonction Â2S (Â) et montrons que sa dérivée est strictement positive et croît avec Â: 00
"2"
2
d dÂ. [l'v S (l'v)] =
r Â-xcp (x) J (X+Â.)8
00
d
Jr
ucp (Â.u) (U+1)3
d
u. o 0 De la dernière expression et du fait que la fonction cp (x) est croissante, il résulte immédiatement que la dérivée du premier membre est strictement positive et croissante. On peut donc appliquer le lemme 1 à la fonction Â2S (Â) et, en vertu de (353), obtenir d dÂ.
x=
2
3 ~ [Â2 S (Â)] ~"2 HÂ2,
d'où 00
-8'(Â)=2~
(354)
On a par ailleurs -
 3s'
r (u+ ucp (Â-u) 1)4
00
:!.
1
d
(Â) ~"2 HÂ2,
-
dÂ
[Â3 S '
(Â)] = 2·3
J
du,
o et l'on peut appliquer le lemme 1 à la fonction - 1..3S ' (1..) : -
d
dÂ.
[1.. 3s' (1..)]
~
1
!
3
2:."2 H)..2.
Une dérivation combinée à (354) nous donne 5
00
S" (")
l'v
= 3.f J\ o
cp (x)
(X+Â.)4
1
3
dx ~ 2'2H1..
- -
2.
(355)
408
CH. II.
PROB~MES
AUX LIMITES
En prolongeant cette procédure, on est conduit à la formule 00
(_1)ffl sm (Â)
=
+ 1) 1 )
(m
(X+~~+2 dx __ {·3 .. ~~m-1) RÂ-
2m+l -2-.
o
(356) Etudions maintenant le comportement de l'intégrale (357) pour de grands Â. On a (x
~
xp ).,)2p+2 -
( 1- x  ) (x+).,)p+1
+
p
-
~
s)
(
LJ
P
.=0
~(~).,).~ (x+).,)P+2+"
où S ) (
= p (p -1) ... (p - s + 1) .
psI
0
\1 p ) = 1 ,
'
et par suite p
J p (Â) =
00
~
( ; ) ( - Â)' ) (X;).,\Xj+2+. dx .
•=0
0
En tenant compte de (356) on obtient 1 p
J (Â) __ RÂ -p- 2 ~ (_1)' ( s) 1·3 ... (2p+2s+1) LJ
p
p
s=o
(p+s+1)12p+,'
(358)
1
Jo (Â) =s (Â) -- RÂ - 2. La première de ces formules peut être mise sous la forme -p_!
J p (Â) -- RÂ
2
1
PaS (-1) (p)
~
vn
r r
~
(p+s+ ) (p+s+2).
(359)
s=o Prouvons maintenant la formule 1
(S
p
../n ~
(_1)8
p)
r r
~
(p+s+ ) 2 K po (p+s+2) =
s=o A cet effet, considérons l'intégrale
rr:
(360)
(361)
409
11-2-46. DEMONSTRATION DU THEOR1lJME AUXILIAIRE
qui par la substitution x= Âu se ramène à la forme 1
p+ u 2 1 , 1 p+ - • (u+ 1)2p+2 du = 00
Lp= Â
2
1
(362)
P+.!. K po Â
0
2
On a p
X) P =.LJ(-1) "" & ( S) Â8 ( x+Â p (X+À)8'8=0
et l'intégrale (361) devient p
Lp =
~
00
(_1)8
(~)
Â8 ) 0
8=0
1
p
1
1
').P+'i
""
.LJ
00-
(_1)8
(Sp) Jr
8=0
u 2
(u+1)P+8+2
d
u.
(363)
0
Le changement u=x/(1-x) nous donne
r J
1
u"2 (U+1)P+8+2 du=
o
rJ x...!.. (1-x) P+8-i-~ dx= r({)r(p+s+{-) 2
r (p+s+2)
0
Portons ceci dans la formule (363) :
V-Lp=
p
~
p:J.. 2').
2
(_1)8
8=0
(~)
r
(p+s+-.!- )
r
(p+s+~
.
En comparant avec (362) on obtient (360) et la formule (359) devient 1
J
P
(Â) ~ 2H Â-P-T K Tt
(364)
P
ou (365) où Tb. dépend de p et de
Â
et ll;l.. -+ 0 pour p fixe et
 -+
00.
CH. II. PROBL:E}MES AUX LIMITES
410
Mettons l'intégrale J p (Â) sous forme d'une somme de quatre termes: (l-a)1I.
(Xx: 0, f-l > 0, f-ll;;;:: 0, f-l2 et ~L3 des nombres quelconques (pas forcément strictement positifs). Posons (u, v) =
l
n
uv dx,
Il u Il = V (u, u), ui: = ~ U~i' i=1
D
1U x 1= Il u
Hv Il Ux 11
2
-Ill" Ux
1111
u 1l-1l311 u 11 2 •
De là il s'ensuit En vertu de (416)
L(u, u»(v-81)lIuxIl2-(1l3+
~~)
lIull 2
(417)
quel que soit 8 1 > O. Supposons que u est une solution distributionnelle de classe W~ (D) du problème (403), de sorte que cette solution est justiciable de l'identité (407). En faisant 11 = u dans (407), on obtient
L (u, u) + Â Il U 11 2 =
.-
~ fu dx.
(418)
D
En se servant de (417), de l'inégalité de Bouniakovski-Schwarz et de (416), on déduit de (418) l'inégalité suivante: (v -- 8 1 )
Il Ux 11 2 + (Â -1l3 - ~tl ) Il U 11 2 ~ ~ Il f" Il u Il ~ë211
11 2+ 48 Il f !1 2, 1
U
2
(419)
où 8 2 > 0 est un nombre arbitraire. Cette inégalité s'appelle première inégalité fondamentale (ou énergétique). Sur (419) on voit que la norme Il U x Il de la solution distributionnelle du problème (403) se majore en fonction de Il fI/et Il u Il. Si (420) alors (419) nous permet d'estimer la norme 1/ U x Il uniquement en fonction de Il f 1/. En effet, prenons par exemple 8 2 = ô/4 et 81 tel que  - Ila - Il~ (481t1 = 3ô/4. Alors par des calculs élémentaires on trouve que 8 1 = vll~ (ôv + Il~) -1, et de (419) on déduit la majoration voulue: (421)
11-2-52. PREMIÈRE IN);:GALIT);: FONDAMENTALE
429
ou, pour abréger, "u
Il (1)
~ Cô
Il fil.
(422)
Grâce à cette majoration on a le théorème d'unicité suivant: Thé 0 r ème 1. Si les coefficients de L satisfont les conditions (413) et si est réalisée la condition (420), alors le problème (403) possède une solution distributionnelle au plus de classe W~ (D) (le domaine D peut- être illimité). En effet, la différence v = u' - u" de deux solutions distributionnelles éventuelles du problème (403) vérifie l'identité (407) avec / = 0, donc l'inégalité (422) avec f = 0, inégalité de laquelle il résulte que v = 0, c'est-à-dire que u' = u". Si le domaine D est borné, on peut affaiblir la condition (420). On se sert à cet effet de l'inégalité
l u2dx~CD l u~ dx D
(423)
D
(l'inégalité (32) du (tome IV1 , [III-7]) valable pour toute fonction o u EW~ (D) et tout domaine borné D. Pour 8 1 E]0, vI il s'ensuit de (419) que ( v-el CD
+ Â - ~l3 -
4l-ti - 8 2 ) Il u el
112~-41 Il f e
11 2 ,
(424)
2
donc, si max
O 0,
(425)
= 61 /2 dans (424), on obtient
Il u 11 ~ 2
612 Il f
2 11
•
(426)
De là ·et de (419) dans laquelle on aura posé par exemple 8 1 = v/2 et 8 2 = 1 on trouve la majoration de la norme complète de u dans W~ (D), soit: Il u 11(1) ~ C Il f Il. (427) La condition (425) est remplie par exemple pour L = â quel que soit  ~ O. La condition (420) au contraire n'est pas satisfaite pour L = ~ et  = O. Si les coefficients sont donnés dans un domaine quelconque Dl et qu'ils satisfont les conditions (413) dans Dl' alors la condition (425) est vérifiée pour tout  fixe (en particulier pour  = 0), pourvu que le domaine D c Dl soit de volume assez petit, car CD -+ lorsque mes D -+ O. On dit pour cette raison que le théorème d 'unici té est valable pour le problème (403) dans des domaines de « volume assez petit ».
°
CH. II. PROBL:E:MES AUX LIMITES
430
Rem a r que. Au (tome IV I , [111-7]), on a montré que CD -+ 0 si le diamètre deD tend vers O. Des raisonnements plus fins montrent que CD est proportionnel à 1 mes D 1 2/ n •
11-2-53. L'espace W;, 0 (D) et la deuxième inégalité fondamentale. Prouvons préalablement l'inégalité
i Ivl dS~q)~(IVI + Ivxl)
(428)
dx.
D
S
Cette inégalité est valable pour toute fonction v E WI (D), c'est-à dire pour toute fonction v E LI (D) possédant des dérivées distributionnelles premières sommables sur D. La frontière S du domaine D doit être assez régulière, par exemple être de classe Cl. L'inégalité (428) est valable pour une classe plus large de domaines: pour des domaines à frontière lipschitzienne. Nous ne donnerons pas leur définition exacte mais de la conclusion déduite plus bas il est aisé d~ voir quelles sont les propriétés de S qui sont suffisantes pour qu~ (428) ait lieu. L'inégalité (428) nous sera utile uniquement pour les fonctions v E Cl (D). C'est pourquoi nous nous bornerons à la démontrer seulement pour de telles v. De ce fait et de la densité d~ Cl (D) dans Wi (D), il s'ensuit que (428) a lieu pour tout v E Wi (D). Considérons un système de coordonnées cartésiennes locales (YI' ... , Yn) ayant pour origine un point quelconque XO ES, et un morceau Sp (XO) de S d'équation
y'
== (YI'
.•. , Yn-l) E B p = {y': 1 y' 1
D E1
,,2 = Be ' l'inégalité (440) devient s
~ \.2" U xx 112~ "L (U) 11 2_ ~ \ 1 (s) dS + C~ ( Il U WU)2. (441) S
Jusqu'ici nous n'avons pas utilisé le fait que u s'annule sur S. Montrons que si l'on fait jouer cette condition, alors l'intégrale ) 1 dS s
peut être ramenée à une forme ne contenant pas les dérivées secondes de u par rapport à x. Ce fait est primordial lors de la déduction de l'inégalité (434). Pour le prouver, considérons un point arbitraire XO E S et introduisons un système de coordonnées cartésiennes locales y: Yh = Chl (Xl - xV, k = 1, ... , n, c'est-à-dire un système dont l'axe des Yn est orienté suivant la normale extérieure n à S en XO et la matrice (Ch l), orthogonale. Soit Yn = (ù (YI' . • ., Yn -1) l' équation de S au voisinage de l'origine des coordonnées Y = (0, ... , 0). Par hypothèse, (ù (YI' ••. , Yn-I) E C 2• La matrice (Ch l) étant orthogonale, on a Xl - xY = ChlYh, l = 1, ... , n. Donc cos (n, Xl) = = Cnz, l = 1, ... , n. Considérons l'expression 1 (s) au point XO et passons dans cette expression aux coordonnées y: 1 (XO) -
= aihajlCmiUYmCpjCqlUYpylnhaihajlCmiuYmc PkCqlUYpylnj
== (bmnb pq -
bmpbqn ) uYmuYpYq'
(442)
où
b pq
=
ajlCpjcql'
p, q
= 1, ... , n.
Utilisons maintenant la condition aux limites U 1s = O. Au VOISInage du point XO dont les coordonnées Y l sont nulles, cette condition devient U (YI' ••• , Yn-l'
(ù
(YI' ••• , Yn-l))
= 0,
et de plus elle est identiquement satisfaite par rapport à YI' ••• • • • , Yn -1 au voisinage de YI = ... = Yn -1 = O. Dérivons cette identité par rapport à Yi et Yh, i, k = 1, , n - 1, et prenons en considération le fait que (ÙYi = 0, i = 1, , n - 1, au point xo. 28*
CH. II.
436
PROBL~MES
Ceci nous donne au point Uv.~
AUX LIMITES
XO:
(443)
= 0,
i, k=1, ... , n-1. Grâce à
:~ /xo=O, i=1, ... , n-1, on a ôu au bnpbqn ) -a- a ô ' (444) n YP Yq 2
1
(xO) =
(bnnb pq -
Pour p = n et q quelconque, ainsi que pour q = n et p quelconque, les termes entre parenthèses dans (444) se simplifient mutuellement, ce qui, combiné avec (443), nous donne pour 1 (XO) l'expression suivante: n-1
1
~
(XO) = -
(bnnb pq -
bnpbqn ) ( ; : ) 2 (OYpY q •
(445)
P. q=1
On admettra que les coordonnées YI' ... , Yn -1 sont choisies dans le plan tangent à S en XO de telle sorte que toutes les dérivées mixtes (OlIpllq' p =1= q, soient nulles en XO (on sait qu'on peut toujours réaliser ceci par une transformation orthogonale des coordonnées YI' ••• , Yn -1)' Alors n-1
l (xO)
= -
~
(bnnbpp-b;p) ( ; :
)2 (OYpYp.
(446)
p=1
En vertu de l'hypothèse SE C2 , on peut exhiber un nombre tel que (OYpYplxo~K,
p=1, ... , n-1,
K~O
(447)
quel que soit XO E S. Si, en particulier, le domaine D est convexe, alors on peut prendre pour K. De la condition (413), il s'ensuit bnnb pp - b~p ~ ~2, et par suite
°
0-
0, si u =f= O.
(455)
Le produit scalaire (u, u) permet de définir un nombre strictement positif appelé norme de l'élément u: plus exactement, Il u Il = = V (u, u). La norme est nulle uniquement pour l'élément nul de H (que nous désignerons simplement par 0). Ce qui précède nous permet d'introduire la notion de convergence: une suite {Uk} k=t converge vers un élément u si lim Il Uk - U Il = O. Une suite {Uk} k-l est k-...oo
°
dite suite de Cauchy (ou suite fondamentale) si Il Uk - Um Il ~ lorsque k et m ~ 00. L'espace H est dit complet si toute suite de Cauchy {Uk}:' 1 converge vers un élément u E H. Enfin, l'espace H est dit séparable s'il existe un ensemble dénombrable d'éléments de H: {Vk} k-1'
440
CH. Il. PROBL1lJMES AUX LIMITES
dense dans H. Ceci signifie que pour tout élément u E H et tout 8 > 0 il existe un élément Vk tel que Il u - Vk Il ~ 8. Si H possède les propriétés énumérées ci-dessus, on l'appelle espace hilbertien complexe séparable complet. Souvent, les propriétés de complétude et de séparabilité sont implicitement incluses dans la notion d'espace hilbertien et dans la suite ne seront pas spécifiées. De même, nous ne les signalerons pas séparément, car on admettra qu'elles sont toutes deux réalisées. Si a, b, ... sont des réels, alors le produit scalaire (u, v) l'est aussi et la conjugaison complexe disparaît des égalités (455). Dans ce cas on parle d'un espace hilbertien réel H. On démontre que dans un espace hilbertien abstrait H, il existe une base orthonormée, c'est-à-dire un système d'éléments {(})k}k' 1 tels que {(})k' (})z} = {hl et tout élément u se représente par une série 00
u =
L
(u,
(})k) (})k
convergente vers u pour la norme de H. De là il
k=l 00
s'ensuit que" u
11
2
=
LJ
1 (u, (})k) 12 •
Il existe une infinité de telles
k=l
bases. Mais en en choisissant une on peut établir un isomorphisme entre deux espaces hilbertiens complexes quelconques Hi et notamment entre H et l'espace hilbertien complexe l2 dont les éléments sont les suites x = (Xl' X 2 , • • • ) de nombres complexes Xk 00
tels que
2:
1 Xk
k=l
/2 < 00. La multiplication par un nombre et
+
l'addition se définissent comme suit: ax = (axI' ax 2 , • • • ), X y = = (Xl YI' X2 Y2' ... ) et le produit scalaire par l'égalité
+
+
00
(X, Y)l2 =
LJ
k=l
XkY k·
L'isomorphisme indiqué se construit comme suit: à un élément u E H on associe la collection ;; = «u, (})I), (u, (})2), ••. ) de ses coefficients par rapport à la base orthonormée {(})k}k' 1 choisie dans H. '" '" Il est immédiat de voir que (u, v) = (u, v) 1 Grâce à cela tous les espaces hilbertiens (il s'agit uniquement des espaces séparables complets) ont la même structure d'un point de vue abstrait et les faits prouvés pour l'un peuvent être reformulés pour l'autre. Les espaces signalés au début de ce numéro sont isomorphes à l'espace l2 et partant l'un à l'autre. Donc, de nombreux résultats établis pour L 2 (D) sont valables pour les autres espaces. Par exemple, l'inégalité de Bouniakovski-Schwarz est un simple corollaire des axiomes (455). Dans le cas d'un espace abstrait H elle devient 1(u, v) 1~ " u Il Il v Il. Cependant nous n'utiliserons pas l'isomorphisme entre les espaces fonctionnels qui nous intéressent, car les faits qu'on en tire sont peu suggestifs et de plus nous n'avons 2
•
11-2-54. QUELQUES
CONSID~RATIONS
SUR LES ESPACES HILBERTIENS
441
pas encore prouvé que les espaces WJ (D) avec l entier sont séparao bles. La séparabilité de ~ (D) et de W;,o (D) résultera des propositions du numéro [II-2-57). Il est aisé de prouver la séparabilité detous les espaces (D) à partir de considérations plus simples et plus générales, ce qui sera fait au tome V, chapitre IV. Dans un espace hilbertien donné H, on peut introduire une nou-· velle structure hilbertienne (c'est-à-dire un autre produit scalaire) de telle sorte que cet espace soit de nouveau complet. Pour cela, il faut que la nouvelle norme Il • 11* correspondant au nouveau produit. scalaire (,)* (c'est-à-dire Il u 11* = V (u, uf*) soit équivalente à la norme initiale Il . Il. L'équivalence des normes Il • 11* et Il . Il signifiequ'il existe deux constantes Cl > 0 et C 2 > 0 telles que pour tous les éléments de Hl' on ait
w1
Cl
Il
u
~
Il
Il
11*
u
~ C2
Il
u
II·
Il est clair que si une suite d 'éléments u k, k = 1, 2, . . . converge· pour la norme Il . Il, alors elle convergera pour la norme Il . 11* et réciproquement. Ceci garantit la complétude de H relativement à la convergence pour la nouvelle norme. Si H est séparable, alors il en est de même du nouvel espace hilbertien engendré par lui. Les opérateurs linéaires sont l'un des principaux objets d'étudedans les espaces hilbertiens. On dit qu'un opérateur A est un opérateur linéaire borné dans H s'il applique H dans H de telle sorteque A (au bv) = aA (u) bA (v) et s'il existe un nombre C >Û' tel que pour tous les éléments de Hl' on ait *)
+
+
Il A
(u)
Il
-< C Il u II·
La borne inférieure de tous les tels C s'appelle norme de l'opérateur A et se note Il A Il. Démontrons une proposition simple qu'en fait nous avons rencontrée à maintes reprises en étudiant des espaces fonctionnels concrets et des opérateurs concrets évoluant dans ces espaces. Lem m e 1. Soit A un opérateur linéaire borné dans un espacehilbertien H et supposons que sa norme Il A Il < 1. Alors l'équation u
=
A (u)
+v
(456}
admet une solution unique dans H quelle que soit v. Pour intégrer l'équation (456), on se sert de la méthode des approximations successives dans sa forme élémentaire, c'est-à-dire: soient UI = v, Un+l = A (un) v, n = 1, 2, ... Il est aisé de voir-
+
*) On écrit souvent Au pour A (u).
CH. II. PROBL:E:MES AUX LIMITES
442
-que " Ul
= Il v /l, Il ~ Il A Il Il Un
/1
Un +1 -
Un
Un -1
Il = Il A
/1 ~ •••
-
O. Il est aisé de *) Au (tome II, [VIJ-1-15], [VII-1-16]) et au (tome III 2 , [VI-1-5] à [VI-1-8], (VI-2-3], (VI-2-11] à [VI-2-13]) on donne la forme explicite de tous les Uh pour de tels D, forme de laquelle il résulte que Uh E Coo {D); la complétude de {uk}~ 1 est prouvée au [11-2-34].
II-2-56. LA R:E:SOLUBILIT:E: AU SENS DE FREDHOLM
447
voir que si u (x) E W:. o (D), alors u (y) = u (x (y)) sera un élément de W;.o (D) et réciproquement. Lorsqu'on passe de x à y, l'opérateur différentiel L (u) s'écrit:
'" ""
L
(u) =
,...,
a (""
OYi
aik
OU)
(y) oYk
""
'"
+ b i (y)
ou
,. ., + '"C (y) u,
0Yi
où , '" aik
bi (y) =
[b
j
(x)
(y)
=
[
a il (x)
°Yi (x) ] 1
OX j
OX j
aYk (x) ] oxz
1 x=x(y) ,
-
X=X(Y)
-r
aZk (x)
,...,
0Yi (x)
f)Yi (x) ] 1
a (a yj
(x) ) 1 dXz x=x(y)
On vérifie sans peine que les coefficients de L oXk
x=x(y)
ay j
,
et C (y) = C (x (y». satisfont les inégalités (413) et (433), mais avec d'autres constantes. Ces constantes nous permettent de choisir un  o assez grand pour ,..., "" que L - ÂoE vérifie la condition (458) dans D, et par suite l' iné,..., galité (459). Si D est une boule, un parallélépipède ou une couronne sphérique, alors de ce qui a été démontré plus haut, il s'ensuit que le problème L (;;) - Âo;;' = (y), ;;, laD = 0, admet une solution ,..., "" unique dans W:. o (D) pour tout f E L 2 (D). En passant dans ce problème à l'ancienne variable x, on s'assure que le problème (471) admet une solution unique dans W;,o (D) pour la valeur de  o retenue et pour toute fonction f E L 2 (D). On a donc prouvé le théorème suivant:
7
Thé 0 r ème 2. Si les coefficients de L de (402) vérifient lesconditions (413) et (433) et si le domaine D est une boule, un parallélépipède ou une couronne sphérique ou est susceptible d'être transforméen l'un de ces derniers par une fonction régulière y = y (x) E C 2 (D),. alors le problème (471) admet une solution unique dans W;,o (D) pour toute fonction f E L 2 (D) et pour  o assez grand. Rem a r que. Ce théorème nous permet de prouver que le problème posé admet une solution unique dans tout domaine D (éventuellement illimité) dont la frontière est une surface de classe C2. 11-2-56. Sur la résolubilité au sens de Fredholm du problème de Dirichlet. Considérons maintenant la résolubilité du problème de Dirichlet L (u) = Âu + f, u 1s = 0, (474) pour  réel arbitraire. Tout ce qui va suivre est valable pour  complexe, mais on se bornera au cas réel. Supposons que les conditions
448
CH. II. PROBLÈMES AUX LIMITES
.du théorème 2 du numéro précédent sont remplies pour L et D. Prenons Âo assez grand pour que LI = L - ÂoE soit justiciable de la -conclusion de ce théorème et ramenons le problème (474) à la forme -équivalente u = (Â - Âo) L -~ (u) L -~ (f). (475)
+
Au [lI-2-55J on a montré que L-~ transforme les éléments de L 2 (D) .en ceux de W 2: 0 (D) et que pour tout f E L 2 (D)
Il
L~l (f) "
::s; C 2 1/ f Il
(2)
(476)
{cf. (461)). En vertu du théorème de Rellich (tome IVI , [lII-7J) il ;s'ensuit de là que L-~ est un opérateur complètement continu dans L 2 (D) (c'est-à-dire comme un opérateur de L 2 (D) dans L 2 (D)). Donc, l'équation (475) est justiciable des théorèmes de Fredholm (tome IVI , [I-28J et tome V, [V-1-14J à [V-1-16J) qui nous disent -que l'équation (475) admet une solution unique dans L 2 (D) pour tous les  à l'exception d'un ensemble dénombrable que l'on désignera par { h } , k = 1, 2, ... Supposons que  =1=  h et que u est la .solution correspondante de l'équation (475). Cette solution est un -élément de W2~O (D), puisque les deux termes du second membre de (475) le sont si f E L 2 (D). Donc, u est une solution du problème (474) appartenant à l'espace W 270 (D). Pour  =  h (et uniquement pour de tels Â), l'équation sans se.cond membre admet les solutions non triviales Uh
=
(Â h -
 o) L-~ (Uh).
(477)
Ces solutions sont aussi des éléments de W 270 (D), car l'appartenance de Uh à L 2 (D) et les propriétés de l'opérateur L-~ entraînent que L-~ (Uh) E W 270 (D). Donc, Uh sont solutions du problème L
(Uh)
=
ÂhUh'
Uh 1S
= O.
(478)
Ces solutions s'appellent fonctions propres de l'opérateur L dans le domaine D pour la condition aux limites du problème de Dirichlet, et Âh , les valeurs propres correspondantes. La théorie des équations linéaires à opérateurs complètement continus nous dit qu'à chaque  h est associé un nombre fini de solutions linéairement indépendantes du problème (478), autrement dit chaque Âh est de multiplicité finie. Pour formuler les autres propositions, il nous faut introduire les opérateurs adjoints de L et de LI agissant dans L 2 (D). L'opérateur A * adjoint d'un opérateur (linéaire) A donné se définit par l'identité (479) (A (u), v) = (u, A * (v)), qui est réalisée pour tous les éléments u sur lesquels est donné A. On admet que le domaine de définition :!lJ (A) de A est l'espace L 2 (D) tout entier si A est borné et un ensemble :!lJ (A) dense dans
1
4i9
II-2-58. LA RaSOLUBILIT&AU SENS DE FREDHOLM
L.iD) si A n'ést pas borné. L'ensemble de définition ~(A *) de l'opérateur A* est composé des éléments v (et d'elix seuls) pour lesquels l'expression (A (u), v) qui est une fonction linéaire de u sur l'ensemble ~(A) peut être mise sous la forme (u, w), où west un élément quelconque de L,(D). L'ensemble ~(A) étant dense dans L,(D), l'élément west unique et il doit être égal à A* (v). Si A est un opérateur borné on démontre que ~(A *) est l'espace Ls(D) tout entier (ou l'espace hilbertien H tout entier si L 2 (D) = == H). Ceci résulte aussitôt du théorème de Riesz sur la forme ~énérale des fonctionnelles linéaires de L,(D) (tome V, [V-1-4]). En effet, pour tout v E L" (D) l'expression (A (u), v) est une fonction continue linéaire u définie sur l'espace L 2 (D) tout entier (c'est-àdire une fonctionnelle linéaire), donc il existe un élément unique w tel que (A (u), v) = (u, w) pour tous les u E L,,(D). Or, ceci signifie que v E ~(A *), w = A* (v) et ~(A *) = L,(D). Si A n'est pas borné, la détermination de 9J(A*) et de la forme explicite de A* est une tâche ardue. Les opérateurs L et LI considérés sont des opérateurs non bornés, définis sur l'ensemble ~(L) = = ~(Ll) = W,~o (D), dense dans L 2(D). On montrera que les opérateurs adjoints L* et LT sont de la forme
L*(u)=-aa (a'k aau (biu)+cu, Li(u) =L* (u)-Âou. (480) Xi xk )--aa Xi (c'est-à-dire sont des opérateurs différentiels adjoints de L et LI respectivement au sens de Lagrange (cf. formule (410» et leurs domaines de définition sont confondus avec W2~O (D). Assurons-nous à cet effet que le problème L~ (u) = f, u Is = 0, (481) où L*(u) est donné par (480),· admet une solution unique dans W,~o lD) pour toute fonction f E L, (D). En effet, l'opérateur L~ remplit toutes les conditions du théorème 2 de [11-2-55], y compris la condition L~ (u, u) ~ cS l Il u Il', cS l > 0 si Âo ~ 1, car Li(u, u)=(-L~(u), u)=
=
J(aikuXkux,~biUUXi-cu2+Âou2) dx=L
1
(u, u).
g
Désignons par (L~) _1 l'opérateur associant à f la solution u du problème (481). Cet opérateur est un opérateur borné dans La(D) mettant L,(D) et W,~o (D) en correspondance bijective. D'autre part, pour toutes fonctions u et v de W,~o (D), on a la relation . (LI (u), v) = (u, L~ (v», (481,.) qui s'obtient par une double intégration par parties du premier membre. En comparant cette égalité avec celle qui a servi à définir l' opé29-01017
450
CB.IL PROBLDlES AUX LIMITES,
rateur A * on voit que les éléments.v de W2~O (D) appartiennent à 9J(A*) = 9J (L~) et A* = L~ est défini par (480). Il reste à véri-' fier que ~ (L~) == W2~O (D), c'est-à-dire à prouver que si pour de~ éléments quelconques v et f de L 2 (D) on a (L 1 (u), v)
= (u, f)
pour toutu E W2~O (D), alors v,E W2~O (D) et f = L~ (v). Déterminons à l'aide de f l'élément w = (L;tif, solution du problème (481) (cette solution est un élément de Wi,o (D» et mettons (u, f) sous la forme (u, L~ (w» = (L 1 (u), w) en vertu de la formule (482). Alors (u, f) = (L 1 (u), w) = (L 1 (u), v), et puisque les éléments L 1 (u) parcourent l'espace L 2 (D) tout entier lorsque u parcourt W2~O (D), il s'ensuit de là que v = w, c'est-à-dire que v appartient bien à W2~O (D) et f = L~ (v). , On a donc démontré que l'opérateur L~, adjoint de l'opérateur non borné L 1 de L 2 (D) dans L 2 (D) et défini sur l'ensemble 3J (L 1 ) = W2~O (D), est également défini sur W2~O (D) et agit Sur les éléments de W 2~ 0 (D) en vertu de (480). L'opérateur L (u) ne différant de L 1 (u) que par un terme additif Âou auquel correspond un opérateur borné, on a 9J (L) = ~ (L 1) = W2~O (D), ~ (L*) = = 9J (L~) = W2~O (D) et l'action de l'opérateur L* sur W2~O (D) est définie par (480). Rem a r que. Quand on dit que l'opérateur différentiel L* est défini par (480), on entend par là la méthode de calcul de L* (u) sur les fonctions u (x). Il faut distinguer cette notion de la notion d'opérateur adjoint d'un opérateur non borné L engendré par un opérateur différentiel L. Dans ce cas il faut indiquer le domaine de définition de l'opérateur différentiel L et trouver le domaine de définition de l'opérateur adjoint. Revenons à l'étude des opérateurs L et L* ainsi que des opérateurs L 1 = L - ÂoE et L! = L* - ÂoE, définis sur W~,o (D) de l'espace L 2 (D). Montrons que l'opérateur adjoint de l'opérateur L-~ de L 2 (D) dans L 2 (D) est l'opérateur borné (L1)-1, c'est-à-dire que
(L-D* = (L1t 1 et 9J«L1)-1) = $«L-D*) = L 2 (D). Ceci résulte immédiatement de faits déjà connus et notamment de ce que l'opérateur L 1 et son adjoint Li mettent en correspondance biunivoque l'espace ~ (L 1 ) = ~ (L~) = W2~O (D) et l'espace L 2 (D) tout entier et vérifient l'identité (481 1) pour tous u, v E W 2~ 0 (D). Lorsque u et v parcourent W2~O (D), les opérateurs L 1 (u) et L 1 (v) parcourent l'espace L 2 (D) tout entier. Désignons L 1 (u) par (C), L~ (v) par 1j), .et mettons l'équation (481 1) sous la forme
(cp,
(L~t11j)
=
(L-~cp, 1j).
11-2-56. LA R:eSOLu:aILIT:e AU SENS DE FREDHOLM
451
Cette égalité 'étant valable pour tous !p, '" E L 2 (D), elle traduit lefait que (L-~)· est égal à l'opérateur borné (L~)-l dans L 2 (D). 'Revenons maintenant au problème (474) et aux problèmes (475)" (477) et (478) qui lui sont reliés. L'opérateur L-i étant complètement continu dans LI (D), ses valeurs propres sont confondues avec celles de son adjoint (L-~)*, c'est-à-dire avec O-k - Âo)' autrement dit l'équation v = ( - Âo) (L-D* (v) (482)\ n'admet des solutions non triviales que pour  =  k , k = 1, 2, . . .' (on rappelle qu'on s'est limité à l'étude des  réels). A chaque Âlt est associé un nombre fini de solutions non triviales Vk de l'équation (482), égal au nombre des solutions linéairement indépendantes del'équation (477) correspondant au même  k • En d'autres termes,. les valeurs propres  k - Âo des opérateurs L-~ et (L-i)* ont mêmemultiplicité. Comme (L-D* = (L~)-\ l'équation (482) équivaut a,a système L~ (v) = ( - Âo) v, vis = 0, qui n'est autre que le système L * (v) = Âv, vis = 0, (483) dont les solutions sont des éléments de W2~O (D), l'opérateur diffé-: rentiel L* étant défini par (480). Donc, le spectre réel des opérateurs L et L* pour la condition aux limites du problème de Dirichlet est au plus dénombrable, chaque valeur propre  k est de multiplicité; finie et les seuls points limites de { k } sont  = + 00. Mais le; théorème 2 du [11-2-55] nous dit que pour  ~ Âo, le problème (478), ne possède pas de solutions non triviales, donc tous les Âk < Âo~ Si nous avions dès le départ envisagé un espace hilbertien complexe L 2 (D) et le problème (474) avec des Âcomplexes, alors en appliquant les théorèmes de Fredholm à l'équation (475) nous aurions déduit que le spectre des opérateurs L et L* est composé d'un nombre, au, plus dénombrable de valeurs propres ayant chacune une multiplicité finie. Des estimations relativement simples du type de celles d& [II-2-52] montrent que le spectre est compris dans une parabole d~ la forme Â"2 = al (a 2 - Â'), où ai sont des nombres réels qui se dé.... terminent sans peine à l'aide de la constante d'ellipticité '\1 et des majorants .... i des coefficients de L, ceci étant, al > 0 etÂ' etl'~ sont les parties réelle et imaginaire de  (Le.  = Â' + iÂ'").. Uneanalyse plus poussée montre que le spectre du problème (474) est: constitué d'un nombre infini de valeurs propres {JI'kl k~l et Re Â.k ~. ~ - 00 pour k ~ 00. Ceci a été étab~i par Carleillan (cf. [II-2~45I) ..\ Au numéro suivant on démontrera cette proposition pour un opé-;. rateur symétrique L, c'est-à-dire lorsque L* = L ou, ce qui est. équivalent, lorsque bi (x) === O. Pour l'instant on se propose de prou-· ver le troisième théorème de Fredholm pour le problème (474)~ 29*
452
CH. ,II. PROBLIlMES AUX LIMITES
Supposons que dans ce problème  est égal à une valeur propre  k • Si ce problème admet une solution u pour f E LI (D) (dans la suite toutes les solutions sont supposées appartenir à W2~O (D)), alors en multipliant la première équation (474) par une fonction arbitraire v E W2~O (D) et en intégrant ensuite sur D, on peut transformer le résultat de la manière suivante:
J
L (u) vdx=
J
uL* (v) dx=Â.
J Jt uvdx+
vdx .
Si v est confondu avec une solution Vk du problème (483) avec =Âk , alors cette égalité devient
Â
=
(484)
) fVkdx=O, D
et par suite elle est nécessaire pour que le problème (474) admette une solution. Montrons que la condition (484) est également suffisante. En effet, le troisième théorème de Fredholm appliqué à une équation pour un opérateur complètement continu L-~ nous dit qu'une condition nécessaire et suffisante pour que l'équation (475) admette une solution dans L'J, (D) est que
O=(L~l(f), vk)=(f, (L~l)* (v,J)=
Â. f Â. k- 0
(f, Vk)'
où Vk est une solution de (482) avec  = Âk • Or le problème (482) est équivalent au problème (483) avec  = Âk et leurs solutions sont confondues, de sorte que cette condition n'est autre que la condition (484). Cette condition est de la même forme pour les Âk complexes. Récapitulons ce que nous venons de démontrer sous la forme du Thé 0 r ème 1. Si les hypothèses du théorème 2 du numéro précédent sont remplies pour L et D, alors le problème (474) admet une solution unique dans wto (D) quelle que soit f E L 2 (D) pour tous les l'réels sauf pour un nombre au plus dénombrable de Âk , k = 1, 2, ... , qui constituent les valeurs propres réelles du spectre de L dans D pour la 'Condition aux limites du problème de Dirichlet. Le problème homogène (478) admet des solutions réelles non triviales pour ces valeurs Âk et pour elles seules et à chaque Âk est associé un nombre fini de solutions linéairrement indépendantes du problème (478). L'ensemble {Âk } est composé de valeurs propres qui sont les valeurs propres réelles du spectre du problème adjoint (483), lesÂkétantde même multiplicité que dans le problème (478). ~es valeurs Âk , k =1, 2, ... peuvent être rangées par ordre' de décroissance et le seul point d'accumulation de {Âk } est·Ze point 1= -' 00. Une condition nécessaire et suffisante pour que le problème (474) (J,l)ec  = ,Â"l possède une solution unique est que soit réalisée, la condi';'
453
11-2-56. LA USOLUBlLIftAU. SBNS DE :PJUmBOLM
tion (484) d'orthogonalité de f à toutes les solutions du problème homogène adjoint (483) avec le même Âk • Sous ces conditions la solution géNil.
nérale du problème (474) est u =
Uo
+ ~CtUkt' où Uo est une solution ~1
particulière du problème (474), Ch des nQmbres arbitraires et Uk t ' i = 1, ... , N k , les solutions du problème (478) avec ÂIl, Ce théorème exprime l'un des énoncés de la résolubilité au sens de Fredholm du problème de Dirichlet (474). Comme indiqué plus haut, le problème (474) s'étudie de façon analogue pour les  complexes. La condition de résolubilité du problème (474) pour ÂIl complexes est de la même forme que (484). Il est utile de remarquer que pour les solutions du problème (474) on a Il u Il 2 ~ c). Il f Il, (485) pour tous les  =1= À k , k = 1, 2, ... , la constante C). dépendant des coefficients de L, du domaine D et du  considéré. Pour les  vérifiant la condition (420), on aurait pu expliciter C). en fonction de caractéristiques relativement simples du domaine D et des para~ mètres ~, t-'t et 'Â. Dans le cas général, lorsqu'on sait seulement que  =1= Âkt' k = t, 2, on peut affirmer l'existence de la constante C,,-, mais sans pouvoir l'expliciter. Il est clair que C,,- tend vers l'infini lorsque  tend vers une valeur propre de {ÂIl}. Signalons qU& la détermination des valeurs propres ÂIl pour un L et un D concrets donne lieu à des calculs assez compliqués. Le théorème t de ce numéro et le théorème 2 du numéro précé~ dent entraînent le
v,
0
•
0'
Thé 0 t ème 2. Si les hypothèses duthéotème 2 du numéro précédent sont remplies pour L et D, alors toute solution distributionnelle du problème (474) de classe ~ (D) est un élément de W2~0 (D). En effet, supposons qu'une fonction u est une solution distributionnelle du problème (474) de classe ~ (D), c'est-à-dire appartient o à ~ (D) et vérifie l'identité (407) pour toute fonction Tl E ~ (D) .. Mettons l'identité (407) sous la forme 0
L (u, 1J) =
+
-
J
(Aou
+ F) 1J dx.
(486)
où F = f ( - Âo) u. Prenons le nombre  o assez grand pour que soit réalisée l'inégalité (420), inégalité qui assure l'existence d'une solution unique du problème (474) avec  - Âo, c'est-à-dire du problème L (v) = Âov + F, v lB -:- 0, (487) la fonction F - f ( - Âo) u étant une fonction connue (puisqueu • est par hypothèse une solution connue qe W~ (D». Il est manifeste
+
454
'CH. II. PROBL~MES AUX LIMITES
que F E L 2 (D). En vertu de (486), la fonction u est une solution d.istributionnelle du problème (487) de classe ~ (D) pour F = f (À - Ào) u. D'autre part, le théorème 2 du numéro précédent nous dît que le problème (487) admet une solution unique dans la classe W 2: 0 (D) (le premier théorème de Fredholm a lieu grâce à (420». En vertu du théorème 2 de (l1-2-52J cette solution doit nécessairement être confondue avec la solution distributionnelle étudiée u o de ,~(D), donc u est u,n élément de W2~O (D), c.q.f.d.
+
+
]1-2-57. Sur le spectre d'un opérateur symétrique. Considérons -des opérateurs différentiels symétriques L de la forme (402), c'està-dire des opérateurs tels que L = L* (l'adjoint L* de L au sens de Lagrange est défini par (480»: L (u) =
f)~i (ai1~ :~
)
+ cu.
(488)
Posons le problème
L (u) = Àu,
u 1s
= 0,
(489)
dans un domaine borné Den admettant que léS coefficients de L et le domaine D remplissent les conditions du théorème 2 de [1I-2-55l. Des résultats prouvés au numéro précédent, il s'ensuit que l'adjoint d'un opérateur non borné L défini sur un ensemble :2J(L) = W2~o(D) dense dans L 2 (D) par (488) est cet opérateur lui-même défini sur le même ensemble W2~O (D). On dit dans ce cas que l'opérateur Lest hermitien (ou auto-adjoint) et on note ceci par L = L*. La forme quadratique correspondant à L s'écrit:
L (u, u)
= -
)
L (u) u dx =
D
)
(aikuxiuXk - cu 2) dx.
D
Comme par hypothèse -c(x)~-. J!3' on a pour ÂO=l-t3 LI (u, U) E = -
)
(L- Ào) (u) dx =
D
=L(u, u)
+Âol/uIl 2 >v )
u~dx>vCÏllluIl2t (490)
D
où CD > 0 est la constante de l'inégalité (423) [11-2-52]. En vertu des résultats de (11-2-55], l'opérateur LI == L o - ÂoE admet un réciproque L~l qui est complètement continu, hermitien de, L 2 (D) dans L 2 (D). Donc, il est justiciable des théorèmes prouvés au (tome IVI , [1-30 à 1-36]). Plus exactement, son spectre est composé d'un nombre au plus dénombrable de valeurs propres réelles que l'on désignera par I1k, k = 1, 2, ... A chaque I1k correspond un nombre fini de solutions linéairement indépendantes de l'équation L;l Uk = J!kUkt (491)
455
II-2-57. SUR LE SPECTRE D'UN OP:eRATEUR SYM:eTRIQUE
°
qui sont des éléments de L 2 (D). L'opérateur L-~ étant inversible, le nombre !..t = n'est pas un point de son spectre (c'est-à-dire qu'aucun !..th n'est nul), donc il existe une infinité de valeurs propres !..tA distinctes. Des propriétés de L-~ il s'ensuit que Uh E W2~O (D) et tout UA est solution du problème 1
L 1 u=-u,
U
fL
18=0,
(492)
pour !..t = !..th et, réciproquement, toute solution du problème (492) appartenant à W2~ 0 (D) est solution du problème (491). Le problème (492) n'est autre que le problème (489) correspondant à Îv = = Îv o !..t _1. De l'inégalité (490), il s'ensuit que les !..th '< et -!..th ~ V-1CV • En effet, pour U = Uk, l'inégalité (490) nous donne
°
+
L l (u h ,
Uh) =
-.
--!-.II Uh I/2~vCÏlll Uh 1/2. fLk
En vertu toujours des propriétés établies au (tome IV1 , (1-30] à '(1-36]) pour les opérateurs symétriques complètement continus, les valeurs propres !..th, k = 1, 2, ... peuvent être supposées numérotées dans leur ordre de croissance. Par ailleurs, dans la suite {!..th} h CX\ il est commode de faire figurer chaque valeur propre en un nombre .d'exemplaires égal à sa multiplicité et d'associer à chacune d'elles une fonction propre Uh normée dans L 2 (D), ces fonctions devant être choisies de manière à être deux à deux orthogonales. D'après ce qui vient d'être dit, on peut a dmettre que
-v_1CV
~ !..tl ~ ft2 ~ •••
0, alors la fonction
°
+
+°
°
°
t
l
w (x, t) = .\ dl' ) G (x,
o
0
ôw
2
~,
a) n
(~,
T)
d~
(31)
est la solution de l'équation 7jt=a
ô2w
ôx 2
+ n (x, t)
qui vérifie des conditions aux limites et initiale homogènes. Tout ce qui vient d'être dit est valable en dimension n. Une démonstration des propositions avancées figure dans le travail mentionné ci-dessus de A. Tikhonov. 11-3-5. Usage de la transformation de Laplace. Nous avons déjà signalée lors de la résolution du système d'équations intégrales (18) qu'on pouvai t user dIa transformation de Laplace. Cette transformation peut être appliquée directement à l'équation différentielle (5). On se servira ici de la transformation unila térale 00
/(s)=) .-stF(t) dt=L 1 (F).
(32)
o Soient données les conditions aux limites (6) et la condition initiale homogène (10). Au lieu de la fonction u (x, t), on cherchera sa transformée de Laplace: 00
les dérivées requises continues et satisfait l'équation (79). Une dérivation de la série de la formule (84) terme à terme par rapport à t nous donne
°
00
t
00
~ bk (t) v k (P) - ~ k=l
ÂhV k
(P) ~ i'kW-t)b k (t') dt'.
h=l
0
La deuxième série (85) étant régulièrement convergente, la sé00
rie
2i
ÂhVh
(P) converge uniformément pour P
EÏ3 et tE [0, Tl. La
k=l 00
somme de la série
2i
h=l
bk (t) vh (P) est égal~ à .il (P ; t), car cette
série converge régulièrement par hypothèse (tome IV l , [1-31]) .. Donc t
00
Ut
(P ; t) =:Tt (P; t) - ~ v h (P) Â h ~ i'kW-t)b k (t') dt',
(86)
0
k=l
Ut (P; t) est continue pour P E B et t E [0, Tl. En remplaçant P par Q dans cette formule, en multipliant ses deux membres par G (P; Q) et en intégrant sur B, on obtient, compte tenu de l'équation intégrale pour Vk (P) :
où
JJG(P; Q)u,(Q; t)dS= B t
00
= SSG (P ; Q) n (Q ; t) dS - 2} B
.
k=l
v k (P) ~ e"'kW-t)b k (t'>.dt', 0
la somme de la dernière série étant égale à u (P; t). Donc u (P ; t) =
- ~ ~ G (P ; Q) ut{Q ; t) dS + ~ ~ G (P ; Q) n (Q ; t) dS. B
B
(87) Comme :Tt (Q; t) possède dans B des dérivées continues, il vient que la dernière intégrale possède des dérivées premières et secondes par rapport à x et y continues dans B et l'opérateur de Laplace de cette intégrale est égal à -:Tt (Q; t).
480
CH. II. PROBLE:MES AUX LIMITES
Mettons maintenant à profit la convergence régulière de la troisième série (85) pour prouver que u (P; t) possède dans B des dérivées premières et secondes continues et satisfait l'équation (79). Posons t
00
w(P; t)= -Ut(P; t)+n(P; t)= ~ Vk(P)Âk~i·k(t'-t)bk(t')dt'. k=l
. 0
La troisième série (85) étant régulièrement convergente, on peut dériver terme à terme par rapport à t la série uniformément convergente ci-dessus: 00
00
Wt (P: t) = ~
(t)
Âhb k
Uk
(P) - ~
k=l
Uk
(P) Â~ ~ i'k(t'-t)b k (t') dt', 0
k=l
les séries intervenant dans cette formule étant uniformément convergentes. En remplaçant P par Q dans la dernière formule, en multipliant ses deux membres par G (P; Q), en intégrant par rapport à Q et en tenant compte de l'équation intégrale pour Uk (P), on ~btient
5~ G (P: Q) wdQ; t) dB = B 00
t
00
= ~ bk (t) U k (P) - ~ U k (P) k=l
k=l
Âk
~ eJokW-t)bk (t') dt', 0
donc, en vertu de (86), il vient
ut< P
; t)
= ) ~ G (P ; Q) wt{ Q ; t) dS, B
d'où il s'ensuit que Ut (P ; Q) possède dans B des dérivées premières continues. La formule (87) montre que U (P; t) possède dans B des dérivées premières et secondes continues et satisfait l'équation !1u (P; t) -
Ut (P; t)
= n (P; t),
ce qui prouve la formule (84). Si l'on ne considère que les solutions distributionnelles de l'équation (79), on peut légitimer cette formule pour des conditions moins astreignantes sur le second membre. Rappelons la définition d'une solution distributionnelle de l'équation (79). Soient D le cylindre mentionné dans [11-3-11, DT la partie de ce dernier limitée supérieurement par le plan t = T. On dit qu'une fonction u (P; t) est une solution distributionnelle de l'équation (79) si toute fonction (J (P; t) possédant dans DT des dérivées premières et secondes continues et
481
II-3-9. SOLUTIONS DE L'11:QUATION DE LA CHALEUR
s'annulant dans un voisinage de BD T est justiciable de la formule
~ ~ ~U(Œxx+ Œyy+Œt)dxdydt= DT
-
J'.\ InŒdxdydt.
(88)
DT
On se bornera à étudier les solutions distributionnelles de la classe C (D T). Supposons que la première série (85) converge régulièrement si P Effet t E [0, Tl. Alors la somme de cette série est égale à n (P; t) et nous avons vu plus haut que la série (84) converge uniformément, de sorte que U E C (D~). Désignons par nn (P; t) la somme partielle de la première série (85) : n
nn (P ; t)
= ~ bk (t) Vh (P), k=1
et par Un (P; t), la somme partielle de la série (84): n
Un (P ; t) = ~ h=l
t Vh
(P) ~ e?-"k 0, d'où il s'ensuit que l'intégrale considérée converge uniformément sur l'intervalle [0, t] par rapport au A B x paramètre x = x' + x "i pour tous les x complexes assez proches de X o et, d'autre part, l'intégrant de cette intégrale est une fonction F1Og • 15 entière de x pour 't E [0, tL De là il résulte que cette intégrale est une fonction holomorphe de x au voisinage de tout point (x, t) intérieur à ABCD (tome 111 2 , [I1I-16]) et en particulier du point M. On peut en dire autant de la deuxième intégrale de la formule (90). Donc, les solutions de l'équation (89) sont des fonctions analytiques de x.
Cette assertion est mise en défaut pour la variable t. En effet, si toute solution u (x, t) de l'équation (89) était une fonction analytique de t, alors les valeurs de u (x, t) sur toute droite d parallèle à l'axe des t et appartenant à la demi-bande de la figure 15 seraient entièrement définies, en vertu du principe du prolongement analytique, par les valeurs que u (x, t) prend sur le segment de d compris dans ABCD. Or ceci est impossible, car les valeurs de u (x, t) dépendent manifestement du procédé qui a servi à prolonger les fonctions 001 (t) et Cù2 (t) données initialement seulement sur les segments AD et BC des droites x = 0 et x = 1. Jusqu'ici on a supposé que f (xo) 0 sur ]0, lL Si f (xo) =;É 0 on peut prolonger cette fonction à un intervalle [a, b] plus large de
==
II-3-9. SOLUTIONS DE V:eQUATION DE LA CHALEUR
483
sorte qu'elle soit nulle en a et b et ensuite la prolonger par 0 au delà de cet intervalle. Composons la différence (~-X)I
b
U-
1
\
2 V nt.
f (~) e-
. 4t
ds.
a
Cette différence est nulle sur le segment AB et elle est justiciable des raisonnements ci-dessus. Il reste à considérer la solution ua (x, t) = 2 ~ nt
(~-X)I
b
Jf (~)
e-
4t
ds.
a
En appliquant le théorème relatif à une intégrale dépendant d'un paramètre (tome II1 2 , [III-16J), on voit que Ua (x, t) sera une fonction régulière de (x, t) au voisinage de tout point situé au-dessus de l'axe t = 0, c'est-à-dire pour t > O. Signalons encore que la formule (90) nous dit que la fonction U possède des dérivées de tout ordre par rapport à t pour x E JO, ZL Il est possible de majorer les dérivées de la solution de l'équation (89) par rapport à t. Supposons que U (x, t) est analytique en x, possède des dérivées de tout ordre par rapport à t au voisinage de x = t = 0 et est impaire en x. On a alors le développement de Maclaurin U
=
t x + 3! x
U1 ( )
Us
(t)
3
+ . . • + (2n U2n+l (t) 2n+1 + + 1)! X •• • ,
(91)
où (n=O, 1,2, ... ).
L'équation (89) nous permet d'écrire U 2n +I
(t)
=
an ( au ) atn ax x=o
=
dnUl (t) dt n •
(92)
Si p > 0 est strictement inférieur au rayon de convergence de la série (91), alors on a l'inégalité (tome 111 2 , [IV-3D u2n+1 (t)
1(2n+ 1)1
1
~
M pn ,
où M > 0 est un nombre arbitraire. De (92) on déduit la majoration suivante pour les dérivées de Ul (t) : dnul (t) I:s:::: M (2n + 1) ! pn • 1 dtn -...;:::: Cette majoration n'assure pas l'analyticité de la fonction Ut (t). Si nous avions la majoration plus forte dnUl (t) I:s:::: M·nl 1 dt n -...;:::: pn , 81*
484
CH. II. PROBLi!.:MES AUX LIMITES
alors la série de Maclaurin de la fonction U1 (t) aurait été convergente et cette fonction, régulière au voisinage de l'origine. II-3-10. Potentiels de simple et de double couche généralisés en dimension un. Au [11-3-2] nous avons résolu un problème aux limites pour l'équation (5) dans une demi-bande limitée inférieurement par la caractéristique t = 0 et latéralement t par les droites x = 0 et x = l. Considérons maintenant un domaine du plan (x, t), limité inférieurement par la caractéristique t = b et latéralement par les courbes li d'équations (fig. 16) 4 1
b'\
x = 0'1 (t) ;
o
x = 0'2 (t) [0'1 (t)
0 assez petit pour que l 'li' (t') - 'li' (to) 1 ~ e pour 1 t' - t o 1 ~ ô, et partageons l'intervalle d'intégration [b, t) en deux parties lb, t o - ô] et [to - ô, t). La fonction exprimée par l'intégrale étendue au premier de ces intervalles est continue en (x o' t o)' il suffit donc de montrer que l'intégrale .
est assez petite quel que soit (x, t) assez proche de (x o, to) ou confondu avec lui. Le module de cette intégrale est ~ à
rJt
e 2a vrn
,
lx-a (t )1 e -
[0'(t')-x]2 4a 2 (1-1')
(t- t')3/2
dt' •
to-6
Supposons que la différence x - cr (t') change de signe k fois au plus, où k > 0 est un entier bien défini, pour t' E [to - Ô, t] et quelle que soit la position du point (x, t) dans un voisinage de (x o' to). Sous ces conditions l'intégrale t
~
lx-a
(t') 1
(t- t')3J2
e
-
[0' {t')-x]2 4a2 (t-I')
dt'
10- 6
se représente par la somme de k intégrales au plus de la forme [0' (t')-x]2
x-a (t') (t- t')3/2
e
4a 2 (t-t')
dt'
487
II-3-10. POTENTIELS DE SIMPLE ET DE DOUBLE COUCHE
qUi diffèrent des intégrales x - a (li+l) 2a Vl - li+1
~
+ 4a
x-a
e- z2 dz
par la quantité
(li)
2a Vi-t.l
-2cr' (t') ---==:::-'- e l/ t - t'
[a (l')-x]2 4a 2 (l-l')
dt' ,
qui est ::::;;; en module à une certaine constante. Donc, l'intégrale étudiée reste bornée lorsque le point (x, t) se trouve dans un voisinage de (x o, t o). Cette intégrale est multipliée par 8, donc la démonstration du fait que le premier terme du second membre de la formule (101) est continu lorsque le point (x, t) traverse l en (x o' t o) s'achève exactement comme au [II-2-2]. En se servant des potentiels mentionnés plus haut, on peut ramener le problème aux limites pour le domaine de la figure 16 à une équation intégrale comme on l'a fait au (11-3-2]. Supposons que les conditions sont les suivantes: u = 0 sur la caractéristique t = b et u = Wi (t) sur li. On cherchera la solution sous la forme 2
t) u ( x"
=
1
~
2ayn L.J
t
Jr
[ai (t')-xJ2
'\Ji
(t')
(t-t')3/2
[x - al' (t')] e
4a 2 (l-l')
dt'.
(102)
i=1 b
Sous ces conditions, quelle que soit 'Pi (t'), l'équation (5) est satisfaite avec la condition aux limites supportée par la caractéristique t = b; quant aux conditions aux limites supportées par li' elles nous donnent le système d'équations intégrales de Volterra suivant pour 'Pi (t):
(103)
-ai(t')]e
488
CH. II. PROBLÈMES AUX LIMITES
COllsidérons les intégrales de la première équation t
r (t'\JI (t') t')3/2
J
[al
[a (t)-a (t')] 1
e-
(t')-al (t)J2 4a 2 (t-t')
1
dt' ,
(104)
dt' .
(105)
b
t
r '\J2 (t') J (t-t')3/2
[a2 (t')-al (t)J2 4a 2 (t- t')
[al (t) -a2 (t')] e
b
Dans l'intégrale (104), la singularité pour t' = t est affaiblie par le numérateur durapport (JI
(t) - (JI (t') (t-t')3/2
exactement comme nous l'avons signalé plus haut. Dans la deuxième intégrale, l'exposant de e tend vers - 0 0 pour t' -+ t et ceci fait entièrement disparaître la singularité. On étudie de façon analogue les intégrales de la deuxième équation (103). Donc, le système (103) possède une solution unique et peut être résolu par la méthode des approximations successives. On peut chercher la solution du problème aux limites posé ci-dessus sous la forme d'une somme de deux potentiels de simple couche: t
2
u (x, t)=
, [[a j(t')-x]2 0 la série (133) converge pour la norme de W~, 0 (B) uniformément en t E [e, 00[, où e > 0 est un nombre arbitraire. En effet, en vertu de (1.32) p
jm.
p
(t)
==" k=m ~ (cp,
Uh) eÎ.ktuk (x)
II~ =
p
2J
k=m
Â~e2Âkt (cp,
Uk)2.
Or pour t>e > 0, les fonctions Â~e2Âkt, k = 1, 2, ... sont ~ à un nombre Ce dépendant seulement de e. Donc pour t>8>0, on a p
im, p (t) ~ Ce 2J
k=m
(cp,
U k )2.
Ceci prouve la convergence de la série (133). De cette convergence, il s'ensuit que la somme U de la série (133) est un élément de W~ 0 (B) dépendant continûment de t > 0 pour la norme de cet espacé. 32-01017
498
CH. II. PROBL:E:MES AUX LIMITES
Une dérivation terme à terme de la série (133) par rapport à t nous donne la série 00
2J
h=1
(134)
'Ah ( 0, où 8 > est un nombre arbitraire. On prouve ceci comme plus haut en prenant en considération le fait que les fonctions l 'Ah Ir /"k t , k = 1, 2, ... , pour t ~ 8 > sont ~ à un nombre C e,r dépendant uniquement de 8 et de r. Une pareille convergence de la série (134) garantit l'appartenance de sa somme à wto (B) pour tous les t > 0, et que cette somme est la dérivée distributionnelle par rapport à t de la somme u (x, t) de la série (133) dans les domaines De, T = B X ]8, T[, 8 > O. De tout ce qui vient d'être dit, il s'ensuit que la somme u (x, t) de la série (133) est une solution distributionnelle du problème (126), (127), (128) dans le domaine DT = B X JO, T[ de la classe mT dont les éléments v (x, t) sont doués des propriétés suivantes: ce sont des éléments de L 2 (B) dépendant continûment de t E [0, Tl pour la norme de L 2 (B) ; ils possèdent une dérivée distributionnelle Vt (x, t) dans DT et de plus v (x, t) et Vt (x, t) sont des éléments de W~. 0 (B) dépendant continûment de t E JO, Tl. La somme u (x, t) de la série (133) est solution de l'équation (126) pour tout t > et pour presque tous les x E B. Elle satisfait la condition aux limites (127) en ce sens que pour t > elle est un élément de W;,o (B). La condition initial e est vérifiée «en moyenne»:
°
°
°
" u (x, t) -
0, cette solution est douée de la régularité requise pour la déduction de l'inégalité énergétique (124), donc elle est justiciable de la majoration ) v 2.dx ~ eC l(t B(t)
e)
)
B(e)
V2
dx
(136)
II-3-t3.
M~THODE
pour tout 8 E JO, t[ et t E [8, TJ. En faisant tendre 8 vers que v 2dx -+' pour 8 -+ 0, on s'assure que v (x, t)
j
499
DE FOURIER
°
°
==
et du fait O. Ce qui
B(e)
prouve le théorème d'unicité. Si l'on soumet {fJ (x) à des conditions subsidiaires, alors la série (133) convergera plus vite, et sa somme sera douée de meilleures o propriétés différentielles. Si, par exemple, {fJ E W~ (B), alors la o série (133) converge pour la norme de W~ (B) uniformément en • t ~ et sa somme sera donc un élément de ~ (B) dépendant continûment de t pour la norme de cet espace pour tous les t ~ O. En effet, en vertu de (131) et (129), on a
°
p
p
\1 ~
({fJ, Uh) i'htU k
II~ = ~ PI/hl
k=m
({fJ,
Uh)2e2"'kt~
k=m
et 00
~l [ 'l', l"" ~;k 1 J=
['l', 'l' l , o
c'est-à-dire que la série (133) converge pour la norme de ~ (B) o uniformément en t E [0, 00[. Montrons que si {fJ E ~ (B), les séries obtenues par une double dérivation terme à terme de la série (133) par rapport à Xh i = 1, ... , n, convergent pour la norme de L 2 (D T), où D T = B X JO, T[ et T est un nombre fini quelconque. Nous avons bien ceci, puisque T
p
~ Il ~ o
({fJ, Uh) e"'ktuk
1/: dt =
h=m p
=
p
T
~ ~ "-l ({fJ, k=m 0
2
U k )2 e "'k
t
dt~ ~ -} l''-hf ({fJ,
U h )2 -+
0
k=m
pour m, p -+ 00. Cette inégalité exprime que la série (134) est convergente pour la norme de L 2 (D). Enfin, le système {Uk}k 1 formant une base dans W~, 0 (B) (11-2-571, il s'ensuit que si {fJ E Wt 0 (B), alors la série (133) converge uniformément en t E [0, oo[ pour la norme de W; (B) et sa somme est un élément de Wt 0 (B) dépendant continûment de t E [0, oo[ pour la norme de W; (B). La série (134) convergera 32*
CH. II. PROBLÈMES AUX LIMITES
500
pour la norme de L 2 (B) uniformément en t E [0, oo[ et sa somme v (x, t) = Ut (x, t) sera un élément de L 2 (B) dépendant continûment de t E [0, oo[ pour la norme de L 2 (B). Nous avons donc prouvé le théorème suivant: Thé 0 r ème 1. Si les coefficients de l'équation (126) et le domaine B remplissent les conditions du théorème 2 de [II-2-55l, c (x) ;:::= ~ et cp E L 2 (B), alors la solution du problème (126), (127), (128) est donnée par la série (133) qui converge uniformément en t E [0, oo[ pour la norme de L 2 (B); pour t > 0, cette série converge uniformément en t E [8, 00[, où 8 > est un nombre arbitraire, pour la norme de w; (B). La série obtenue par une dérivation terme à terme de la série (133) par rapport à t converge uniformément en t E [0, oo[ pour la norme de L 2 (B). La somme de la série (133) est une solution distributionnelle du problème de classe T (T quelconque) et de plus cette
°
°
:m
o
classe est justiciable du théorème d'unicité. Si cp
E W~ (B), alors la 0,'11-
série (133) converge uniformément en t E [0, oo[ pour la norme de W~ (B) et les séries obtenues par une simple dérivation terme à terme par rapport à t et une double dérivation par rapport à x convergent pour la norme de L 2 (D T), DT = B X ]0, TL Enfin si cp E W~, 0 (B), la série (133) converge pour la norme de W: (B), et la série (134), pour la norme de L 2 (B) uniformément en t E [0, 00[.
iConsidérons encore le problème L (u)
=f
(x, t), u
1ST
= 0,
~ u 1 t-o
= 0,
(137)
où L est le même que dans (126) et f E L 2 (D T)' Ce problème est formellement satisfait par la somme de la série t
00
U (x, t)
= ~ ~ f k Ct) i"k(t-T.) d't Uk (x),
(138)
k=l 0
où fk ('t)
= ~
f (x, 't) Uk (x) dx. Montrons que la série (138) et les
B
séries déduites par une simple et une double dérivations terme à terme par rapport à x convergent pour la norme de L 2 (D T)' de sorte que la somme u (x, t) de la série et ses dérivées U X z• et u Xoxlt appartiennent z à L 2 (D T). Pour cela il suffit de s'assurer que T
jp,
m
== ~ o
m
Il:2}
t
~
k=p 0
fk ('t) e"k(t-T.) d'tuk (x)
Il: dt
(139)
501
II-3-t3. M:E:THODE DE FOURIER
tendent vers 0 pour p, m -+ grâce à (132), on a t
T
m
i p , m = ~ ~ Â~ 1 ~ k=p 0 m
IR (l:) /'k(t-'t) dl: 1 dt~ 2
0
T
t
t
~ ~ ~ Â~ (~ f~ (l:) k=p 0
0
t)
t
T
i in
k=p
f (x,
0
dl:) dt~ m
(l:)
dl:
/'k(t-'t)
dt~ ~
i f~ T
(l:) d't,
k=p 0
0
EL 2 (DT) satisfont les égalités
T
T
f2
i"h(t-'t)
0
~ ~ IÂkl et les fonctions
(i
dl:)
i"k(t-'t)
m
lJ
Ce qui est immédiat, puisque
00.
(x, l:) dx dl:
I~ n
=
00
T
00
i
= ~
(l:) dl:
f'k
(l:) dl:.
oB 0 k=l k=l 0 De la majoration de jp,m il s'ensuit également que la série déduite par une dérivation terme à terme de la série (133) par rapport à t converge pour la norme de L 2 (D). Les sommes finies N
UN
(x, t)
= ~
i t
fk
(l:)
i"k(t-'t)
dl: Uk (x)
k=l 0
satisfaisant (pour presque tous les (x, t) de DT) l'équation L
(UN)
=
N
=
où
fN,
fN
(x, t)
=
~
fk
(t)
Uh
(x), et
fN
convergent vers
f
pour
k=l
la norme de L 2 (D T), la somme de la série (138) sera solution de l'équation (137) pour presque tous les (x, t) de DT' Pour que les conditions initiale et aux limites de (137) soient remplies, il faut que la série (138) converge pour la norme de W~ (B) pour tout t ~ O. En effet, les raisonnements qui nous ont servi à majorer jp,m nous donnent 00
I~
t
~
k=l 0
t
~ IÂ k Il
i
k=l
0
00
fk
(l:)
/'h(t-'t)
dl: Uh (x)
Il: =
00
~~
fk
Ct) i"k(t-'t) dl:
h=l 0
2
t
t
i
1~
f'k
(l:) dl: =
i i If
(x, t) 12 dx dt.
0 B
De là il est évident que pour tout t ~ 0, la somme o élément de W~ (B) et Il U (x, t) 111 -+ 0 avec t. Nous venons de prouver le
U
(x, t) est un
502
CH. II. PROBL:RMES AUX LIMITES
Thé 0 r ème 2. Si L et B remplissent les conditions du théorème 2 [II-2-55], f E L 2 (D T), DT = B X ]0, Tl, alors la solution du problème (137) dans DT est donnée par la série (138). Cette série et les séries déduites par une simple dérivation terme à terme par rapport à t et x et une double dérivation par rapport à x convergent pour la norme de L 2 (D T)' La somme de la série (138) satisfait l'équation (137) pour o presque tous les (x, t) E D T' appartient à W~ (B) pour tout t E [0, T] et " u (x, t) III -+ pour t -+ O. Signalons que les solutions des problèmes de Dirichlet et mixte avec des conditions initiale et aux limites pour l'équation L (u) = = f (x, t) sont données par des séries de la forme (133) et (138) dans lesquelles il faut remplacer Uh (x) par les fonctions propres de L correspondant à la condition aux limites considérée.
°
11-3-14. Deuxième inégalité fondamentale et existence de la solution du problème de Dirichlet-Cauchy. Décrivons une méthode de démonstration de l'existence de la solution du problème (113), (127), (128) dans la classe des fonctions W~, 1 (D T), DT = B X ]0, Tl, composée de toutes les fonctions u (x, t) appartenant à L 2 (D T) avec leurs dérivées distributionnelles Ut, U Xt et u xixk ' Cet ensemble peut être traité comme un espace hilbertien complet muni du produit scalaire (u,
v)~~'.J~ = ..\ (uv + UxV x + UxxV xx + UtVt) dx dt
(140)
DT
(on utilise les abréviations de [II-2-52J et [11-2-53]). Désignons la norme de W;,l (DT) par 11·11 ;~rJ~. L'adhérence pour la norme de W;,l (DT) de l'ensemble de toutes les fonctions de C2 (DT)' nulles sur la surface latérale ST du cylindre D T' est un sous-espace de W;,l (D T) qui sera désigné par W~:~ (D T)' On admettra que le domaine B remplit les conditions du théorème 2 de [II-2-55] et que les coefficients de M satisfont les conditions (114), (115), (116) et ôaih (x, t)
1
at
1::;:::: -..;;::~
(141)
Ils'
Montrons que les opérateurs paraboliques M vérifient une inégalité de nature voisine de l' inégalité (434) de [I -2-53]. Considérons à cet effet l'intégrale
~ (M(u»)2dxd1:= ~ Dt
[U'{-L.(U)]2 dx dl",
Dt=BxJO, t[,
Dt
où L (u) ,
ô = -ôXi (
ih
(x, t)
U Xk ) -
bi (x, t)
U.r; i -- C
(x, t) u,
II-3-14. DEUXI:SlME IN:eGALIT:e FONDAMENTALE
503
et U (x, t) est Une fonction arbitraire de C2 (D T), nulle sur ST. Transformons cette intégrale à l'aide de la formule d'intégration par parties (107) de [1-2-19J:
l
[u~ + (L (U))2 + 2aihuXhuxit +
(M (U))2 dx dT: = )
Dt
Dt
+ 2 (biu Xj + cu) u,J dx dT: =
l [u~ +
(L (U))2+.
Dt
+ :1:
(aihUXiuXh) -
Ô;~h
UXiU xh + 2 (biu Xi
+ cu) UT] dx dT:.
(142)
De cette égalité et des conditions (116) et (141), on déduit
l
aihuXiuXh dx+
~ [u~ +. (L (U))2J dx dT:~ ~ Dt
B(t)
+ Cl ~
aihUXjUXk dx+
B(O>
[eui + (1 + e- i ) (ui + u 2)J dx dT:
+~
(M (U))2 dx dT:,
(143)
~
Dt
inégalité dans laquelle la constante Cl ne dépend que de !L2 et !Ls. ~ > 0 est un nombre arbitraire et B (t), B (0) désignent respectivement les bases supérieure et inférieure de Dt. Utilisons maintenant l'inégalité (452) de [II-2-53J. Cette inégalité et la condition (114) nous permettent de déduire de (143) l'inégalité v
~ u~dx+ ~ [u~+ B(t)
:2 (uix+U~+u2)J dxd-r~
Dt
~ f-t ~ U~ dx + C2 ~ [eu~ + (1 + e- i )
(ui + u 2)J dx dT:
+
Dt
B(O)
+ ) (M (u))2dx dT:,
(144)
Dt
où la constante C 2 dépend de v, f-t, f-ti et du domaine B. Faisons e = (2C 2)_1 dans (144). Une majoration grossière de (144) nous donne
~ u~ dx + ~ (2u~ + uix + u~ -+- u2) dx dT:~ B(t)
Dt
~C3[ ~ u~dx+ ~ (u~+u2)dxdT:+ ~O)
Dt
l
Dt
(M(U))2 dx d-rJ.
(145)
·CH. II. PROBUlMES AUX LIMITES
504
L'inégalité
u2dx~ ) u 2dx+) (u~+u2):d't,
) B(t)
(146)
Dt
B(O)
qui se déduit aisément à partir de la formule de Newton-Leibnitz pour u 2 (x, t) à l'aide de l'inégalité de Bouniakovski-Schwarz (cf. (198) [1-2-27]) nous permet de tirer de (145): )
J(ui + u~:c+u~ + u 2) dx d't~
(u~+ u 2) dx +
Dt
B(t)
~c~ [) (U~+U2) dx+ ) (U~+U2) dx dl' + Dt
B(O)
) (M (U»2 dxd't
J.
(147)
Dt
Le lemme de [1-2-27] nous permet de déduire toujours de (145) la majoration désirée )
(u~ + u 2) dx + ) (u~ + u;~ + ui: -t- u 2) dx dl' ~ Dt
B(t)
~ C~eC4t
[i
(U~ + u 2) dx +
B(O)
i
(M (U»2 dx dl'
J
(148)
Dt
(la marche est presque la même que pour (124». L'inégalité (148) n'est autre que la deuxième inégalité énergétique pour les opérateurs paraboliques M pour la condition (127). Cette inégalité a été établie pour toute fonction u (x, t) E C 2 (DT)' nulle sur ST. Montrons que pour de telles fonctions les intégrales () (u~ u 2 ) dx t E
f/2,
+
B(t)
Il u 1I~~'ri~. Considérons à cet effet une régulière, égale à 1 pour tE [~ , T et à 0
E [0, Tl, sont majorées par fonction X (t) ~ pour t E [ 0, ~
i
°
J et utilisons
J
les égali tés
t
(u~ + u 2) X dx =
i i :,; [(ui: +
u 2) xl dx dl' =
0 B('t)
B(t)
=)
(2u x u x -rX + u~X'
+ 2uu-rX + u 2x') dx dl' =
Dt
= ) ( - 2u't~uX Dt
+ uix' + 2uu-rX + u X') dx dT. 2
11-3-14. DEUXUlME INÉGALITÉ FONDAMENTALE
De là il s'ensuit pour tE [
~ ,
TJ :
~ (ui+u 2) dX~C5 ~ (u~+u~x+u~+u2)dxd't'. B(t)
505-
DT
(149)1
J.
On établit une inégalité similaire pour tE [0, ~ La démonstra-tion est la même sauf qu'il faut remplacer la fonction X (t) parune fonction régulière positive égale à 1 Sur [0, ~ ] et à sur-
°
[3{ , TJ.
Grâce à cela la norme
équivalente à la norme "u
11* = max [~ O~t~T
(ui
B(t)
1I·1I~~·~~ de w~: 6(DT) est-
+ u 2) dx f /2 + Il u II~~'~~.
(150).
Ce fait et l'inégalité (148) prouvée pour les u E C 2 (D T), nuls surS T' entraînent la validité de (148) pour tout élément u E W~:6 (D T). Cette inégalité permet de prouver l'unicité de la solution du problème (113), (127), (128) quelles que soient f (x, t) E L 2 (D) et cp (x) E o E w~ (B), en utilisant la méthode de prolongement par rapport au paramètre ([11-2-55] et le théorème 2 de [11-3-131). Plus exactement,. il nous faut considérer les problèmes
==
+
M,; (u) (1- 't') Mo (u) 't'Ml (u) = { ulsT-O, Ult=o-cp, 't'E[O, 1],
f,
(151)1
n
où Mo (u)
= Ut -
2J i=1
UXiX • et Ml (u)
= 111 (u), ainsi que le couple-
t
o
d'espaces hilbertiens ~:~ (DT) et W = L 2 (DT) X W~ (B). L'espaceWest constitué des couples de fonctions {f (x, t) ; cp (x)} et est muni du produit scalaire
(ff;
cp},
ft; ~})w =
l ff
dx dt
+ l (CPxq;x + cp~) dx.
DT
B
Les problèmes (151) peuvent être interprétés comme une familled' équations opératorielles A'( (u) = {f; cp},
't' E [0, il,
où les opérateurs A'( sont définis par A,; (u) = {M'( (u); u It=o},
't' E [0, 1].
(152} (153}
Les opérateurs A,; sont des opérateurs de W~:6 (D T) dans W. L'unicité· de la solution de l'équation (152) pour 't' = et pour {f; cp} E W quelconques est acquise grâce aux théorèmes de (11-3-131. A partir de
°
CH. II. PROBUiMES AUX LIMITES
là on démontre sans peine l'unicité de la solution de tous les problè;mes (152) quelles que soient {f; cp} E W à l'aide des inégalités «149) et de l'inégalité (148), valable pour tous les M.r;, l' E [0, 1] avec une constante C 4 que l'on peut choisir commune pour tous les l' E E [0, il. Nous glisserons sur cette démonstration, car elle est calquée :sur celle du théorème 1 de [II-2-551. Enonçons seulement le résultat :final: Thé 0 r ème. Si les coefficients ,conditions (114), (115), (116) et (141) et 2 de [II-2-55J, alors le problème (113), .unique dans W~', A0(D T), DT =1= B X lE L 2 (D T) et cp E W~ (B).
de M de (113) remplissent les le domaine B, celles du théorème (127), (128) admet une solution ]0, T[ quelles que soient f E
11-3-15. Equations hyperboliques de forme générale. Inégalité 'énergétique pour le problème de Dirichlet-Cauchy. Les numéros suivants de cet ouvrage seront consacrés au problème de DirichletCauchy pour les équations hyperboliques. Considérons l'équation n
n
M(u)= ~ aikuX'Xk+~ biuxi+cu-Utt=f i, k=l
~
i=l
(154)
f(Jui est de la même forme que dans [1-2-27] dans un domaine DT = = B X ]0, T[ d'un espace euclidien Rn+l (x E B c Rn, t E ]0, TO. Les coefficients de l'équation (154) et le second membre f peuvent -dépendre de (x, t). On supposera que ces coefficients sont des fonctions mesurables bornées sur DT et de plus que aik sont dérivables . , ---at' 8aïk 8aih ( , , 1 t d'IS t rl. ;par rapport a" x et t et 1eurs d'erlvees ~ genera emen butionnelles) sont bornées dans DT' Supposons que f E L 2 (DT)' L 'hyperbolicité de l'équation (154) est assurée par la condition n
aih (x, t) ~i~k~ V
2J sr
(155)
i=l
-dans laquelle v > 0 est une constante, ~i des paramètres réels -arbitraires (on convient comme partout que aik = aki)' Posons le problème de Dirichlet-Cauchy pour l'équation (154) dans le domaine .D T' Ce problème consiste à trouver la solution u (x, t) de l'équation (154) dans DT qui satisfait les conditions initiales
ult=o=cp(x),
utlt....o=,p(x),
xEB,
(156)
:et la condition aux limites u
1 ST
= 0,
ST
= S
X [0, T],
(157)
.où S est la frontière de B. Autrement dit, nous cherchons une solution u (x, t) de l'équation (154) dans B c Rn pour t E ]0, T[ vérifiant
II-3-15. EQUATIONS HYPERBOLIQUES DE FORME GENÉRALE
507
les conditions initiales (156) à l'instant t = 0 et s'annulant sur la frontière de B pour t E JO, TL Montrons que ce problème est bien posé, c'est-à-dire qu'il admet une solution au plus et que cette dernière dépend continûment des conditions initiales et du second membre de l'équation. La continuité est comprise au sens des normes intégrales définies par l'inégalité énergétique que nous allons établir. On prouvera l'existence de la solution du problème (154), (156), (157) au numéro suivant en imposant aux coefficients de l'équation (154) des conditions supplémentaires permettant de construire la solution sous la forme d'une série de Fourier. Dans la suite, on aura affaire à des solutions distributionnelles de la classe W~ (D T). Les éléments U (x, t) de cette classe satisfont l'équation (154) pour presque tous (au sens de Lebesgue) les points (x, t) de DT. Ils appartiennent à W~ (B) pour tous les t E [0, Tl et dépendent continûment de t pour la norme de W~ (B) (c'est-à-dire 1/ U (x, t At) - U (x, t) Il "i,B -+0 pour At -+ 0 et t, t At E [0, Tl). La dérivée Ut (x, t) est un élément de L 2 (B) pour tous les t E [0, Tl et dépend continûment de t pour la norme de L 2 (B). Pour cette raison les conditions initiales seront les suivantes: 11 U (x, t) - fP (x) 1/ ~~)B -+ 0 et Il Ut (x, t) - '" (x) 112,B -+ 0 o pour t -+ 0 et la condition aux limites (157) deviendra U (x, t) E W~ (B) pour tous les t E [0, Tl (pour plus de détails sur les espaces W~, voir tome V chap. IV). On supposera naturellement que
+
+
+
o
fP (x) E W~ (B), '" E L 2 (B).
(158)
L'inégalité énergétique pour les solutions du problème (154), (156), (157) (nous omettrons plus bas le terme « distributionnel », car il sera partout question de solutions distributionnelles de la classe W; (D T) sauf, mention expresse du contraire) s'établit de la même manière que l'inégalité énergétique pour les solutions du problème de Cauchy [1-2-271. Sa déduction est même plus simple ici, car les intégrales étendues à la surface l~térale ST sont nulles en vertu de la condition (157). De nombreuses notations et majorations sont empruntées au [1-2-271. Plus exactement, désignons par K (t) l'intégrale K (t)
=
(aikuxiuXh
)
+ un dx,
B(t)
où B (t o) est la section du cylindre DT par le plan t = t o, et par D tt le cylindre B X JO, t l [, en admettant que t i ~ T. Multiplions l'équation (154) par -2u t et intégrons le résultat obtenu sur Dt: - ) 2M (u) Dt
Ut
dx dt = - ) 2fut dx dt. Dt
(159)
508
CH. II. PROBLÈMES AUX LIMITES
Intégrons le premier membre de (159) par parties en procédant comme suit
:t (Ut)2 dx dt =
) 2UttUt dx dt = ) Dt Dt
~\ ll~ dx B(t)
) 2aikUXiXk . Ut dx dt = ) (2aikUxiUtxk + 2 Dt Dt
-
= Jr r_at()
(aikUXiuXk) -
)
ui dx,
B(O)
~;:
UxiUt) dx dt =
éJai"k ] d a t UXiU xk + 2 éJaik éJxk uxiUt x dt =
Dt
Pour établir ces relations on s'est servi du fait que uls T = 0, donc que Ut IS T = O. La relation (159) équivaut donc à la suivante: K (t)
=
K (0) +
) [ éJ;;h
UXiU Xk - 2
~:~h
UxiUt
+
Dt
+ 2biuxiUt + 2cuut -
2fut ] dx dt.
(160)
L'intégrale du second membre de (160) se majore comme au [I-2-27) (le fait que Dt est un cylindre et non pas un cône importe peu). Nous ne reprendrons pas ici ces raisonnements et écrirons seulement l'inégalité énergétique déduite à partir de (160). Cette inégalité (dont l 'homologue est l'inégalité (201) dans [I-2-271) s'écrit )
(ui
+ aikuxiuxk + u 2) dx~
B(t)
~eCt[ ) (u~+aikuX,uXk + u 2 ) dx+ ) j2dx dt]. B(O)
.
(161)
Dt
La constante C dépend dES coefficients de M, plus exactement du nombre v de la condition (155) et des maximums des modun
les des fonctions éJ:~k, ~ (~:ikk - bi ) et c dans le domaine Dte k=l
Les fonctions du second membre de (161) sont connues à partir des conditions du problème, en particulier l'intégrale étendue
509
II-3-16. M:eTHODE DE FOURIER. :eQUATIONS HYPERBOLIQUES
il B (0) est égale à
l ('1'2 +
aikcpXiCPXk +cp2) dx.
B
De (161), il s'ensuit
l
(uï
+ vui + u
2
)
dx~
Blt)
~ eCt
[ )
('1'2 + ~cpi + cp2) dx
B
+ .\ f2 dx dt],
tE [0, Tl, (162)
Dt n
OÙ
~ est la constante de l'inégalité aiksisk ~ ~ ~
i=l
s~ (de sorte que
!.t est un majorant pour aih). L'inégalité (162) s'appelle aussi inégalité énergétique. On en déduit le
Thé
0
r ème
1. Le problème (154), (156), (157) admet une
solution au plus de W~ (D T) qui dépend continûment des conditions initiales et du second membre. En effet, si Ui (x, t), i = 1, 2, sont deux solutions de W~ (D T) correspondant aux conditions initiales CPi (x), 'Pi (x) et aux seconds membres (aux forces) fi (x, t), i - 1, 2, alors leur différence u (x, t) est une solution du même problème de W~ (D T), correspondant aux conditions initiales cp (x) = CPI (x) - CP2 (x), 'P (x) = 'Pl (x) - 'P2 (x) et au second membre f (x, t) = fI (x, t) - f2 (x, t). Donc, u (x, t) est
justiciable de l'inégalité (162), d'où résultent les deux propositions du théorème. 11-3-16. Méthode de Fourier pour les équations hyperboliques. On se propose de prouver l'existence d'une solution du problème de Dirichlet-Cauchy pour les équations hyperboliques de la forme n
i, h=l
+
Supposons que la partie elliptique L (u) = :Jô (aihux) cu UXi h de cette équation et le domaine B remplissent les conditions du théorème 2 [II-2-551, donc que le système de fonctions propres {Uh (X)}~l et de valeurs propres {Î"h} de l'opérateur L pour la condition u 1s = 0 est doué des propriétés décrites dans [II-2-57J. Supposons provisoirement que c (x) ~ O. Ceci exprime que  h < 0 de sorte qu'il est commode de désigner Âh par -~~, en admettant que lth > O. Les solutions u (x, t) = X (x) T (t) de l'équation (163) vérifiant la condition aux limites (157) sont toutes de la forme (ah cos ~ht bit sin ~ht) Uh (x), où ah et b h sont des constantes
+
510
CH. II.
PROBL~MES
AUX LIMITES
arbitraires et Uh (x), la fonction propre associée à la valeur propre Âk
=
-fJ~·
De ce fait, on cherchera la solution du problème (163), (156), (157) sous la forme de la série 00
U (x, t) = ~ (ah cos fJkt h=1
+ bk sin fJht) Uh (x),
(164)
dont on déterminera les coefficients ah et b h à partir des conditions initiales (156). La première de ces conditions nous donne ah = (cp, Uh),
(165)
bA = _1_ ('1', Uk).
(166)
et la deuxième, JLk
Formellement cette série satisfait toutes les conditions de notre problème. On se propose d'étudier sa convergence et en particulier de montrer qu'on peut la dériver terme à terme deux fois par rapport à x et à t. Cette dérivation est nécessaire pour justifier les égalités 00
M (u)
=
M
(2J (ah cos fJh t + bh sin fJh t) Uk (x») = k=1 00
= ~ M (ah cos fJkt -I- bh sin fJht) Uk h=1
(x»
=
0,
(167)
c'est-à-dire s'assurer que la somme de la série (164) est solution de l'équation (163). On montrera que la série (164) et les séries obtenues par une double dérivation terme à terme par rapport à x et t convergent pour la norme de L 2 (B) uniformément en t E [0, 00[. Ceci prouvera que la somme de la série (164) est solution de l'équation (163) pour presque tous les (x, t) de DT (et même plus, pour tous les t E [0, co[ pour presque tous les x E B). Les conditions initiales seront satisfaites au sens suivant:
Il U (x,
t) - cp (x)
IIh2.) B
-+
0,
Il Ut (x,
t) -
'1' (x)
II~~) B -+
°
(168)
pour t -++. O. La condition aux limites (157) sera exprimée par l'appartenance de u (x, t) à W~, 0 (B) pour tous les t ~ O. Tout o ceci aura lieu si cp (x) E 0 (B) et '1' (x) E W~ (B). En effet, on a vu au [II-2-57] que la fonction cp (x) se décompose en la série
wt
00
de Fourier cp (x) =
2J
(cp, uhfuh (x) convergente vers elle pour la_
h=1
norme 11·112 et de plus 00
Il cp Il; = ~ (cp, k=1
00
Uh)2 Jt~
=
2J
h=l
a~fJ~,
(169)
11-3-16.
M~THODE
DE FOURIER. :eQUATIONS HYPERBOLIQUES
5tl
DO
et la fonction '1' (X), en la série de Fourier '\j? = convergeant vers '\j? (x) pour la norme
/1·/11 et
00
/1 '\j? /I~ = ~
k=l
2J ('\j?,
Uk) Uk (x),.
k=1
00
('\j?, Uk)211~
=
2J b~fl:· k=l.
(170).
D'autre part,
Il k=m f (ak cos IIk t + bk sin IIh t ) Uk II~ = p
=
2J
k=m
p
(ak cos IIk t + bk sin IIk t )2 111 ~ 2
2J
h=m
(a: + b:) (-tt,
(171»)
p
Il k=m ~ :t (ah cos IIh t + bh sin (-th t ) Uh Il: = p
=
p
~ (-t~(-ahsin(-tht+bhcos(-tkt)2~t~~2 ~ (a~+b~)111 (172p
k=m
k=m
et p
;;1 (ak cos (-tk t + bk sin J.tk t ) Uk W = k=m
Il ~
p
=
p
~ (ah cos (-tk t + bh sin IIh t )2 (-tt ~ 2 ~ (a~ + b~) (-th.
(1731
h=m h=m En comparant les majorations (171), (172) et (173) à (169) et (170, on s'assure de la validité des propositions concernant la convergence de la série (164) et des séries obtenues par dérivation terme à terme' de (164) par rapport à x et t. Ceci étant, on se rappellera que les normes II· 1/1 et 11·112 sont équivalentes aux normes initiales des espaces~ o W~ (B) et Wt 0 (B) respectivement. Que pour tous les t ~ 0, la somme U (x, t) de la série (164) appartienne à W~,o (B) et Ut (x, t) et Utt (x, t) o respectivement à W~ (B) et L 2 (B) résulte de la convergence de cesséries et du fait que leurs sommes partielles sont des éléments deo W~, 0 (B), W~ (B) et L 2 (B). On a donc prouvé le Thé 0 r ème 1. Supposons que les coefficients de M et le domaineB remplissent les conditions du théorème 2 de [II-2-55J et c (x) ~ o. o Si cp E W~,o (B) et'\j? E W~ (B), alors lasomme U (x, t) de la série (164), dont les coefficients ak et bk sont définis par les formules (165) et (166)" est une solution du problème (163), (156), (157). Cette solution appar-
:512
CH. II.
PROBL~MES
AUX LIMITES
°
.tient à W~. 0 (B) pour t ~ et dépend continûment de t pour la norme de cet espace. Ses dérivées Ut (x, t) et Utt (x, t) sont des éléments de o
W~ (B) et L 2 (B) dépendant continûment de t ~
°
pour les normes de ces espaces. Les séries obtenues par une double dérivation terme à terme par rapport à x et t de la série (164) convergent pour la norme de L 2 (B) .uniformément en t ~ O.
Rem a r que 1. Si l'on écarte la condition c (x) ~ 0, lt théorème reste en vigueur mais il est possible que quelques-unes des -premières valeurs propres Âk soient strictement positives ou nulles. Les termes correspondants seront alors de la forme (akeV"'k i bke- V"'k t ) Uk (x) ou (ak bkt) Uk (x). La convetgence de la série (164) n'est pas affectée par ces termes: en effet, leur somme peut ,être représentée par un seul terme dans (164).
+
+
+
Rem a r que 2. Les problèmes de Cauchy et de Dirichlet.. 'Cauchy pour les équations hyperboliques se résolvent de façon analogue pour t > ou t < O. La série (164) converge comme ci-des:sus pour ·t ~ O. Considérons maintenant l'équation avec second membre
°
M (u) =
f
(x, t),
(174)
'Üù M est donnée par (163). Trouvons la solution qui vérifie les conditions initiales et aux limites homogènes suivantes: u It=o =
Développons à cet effet
0,
Ut
It=o
=
0,
U
IS T = O.
f en une série suivant
{Uk }h' 1: 00
00
t (x,
(175)
t) = ~ (f (x, t) 1
Uh
h=l
(x))
Uh
(x) ===
LJ
fh
(t)
Uh
(x),
h=1
et trouvons les solutions U (x, t) = X (x) T (t) de l'équation M (u) = fh (t) Uh (x) (176' qui satisfont les conditions (175). Posons u (x, t) = De (176), on obtient alors pour Th (t) l'équation
Th
(t)
Uh (x)~
(177)
On sait que la solution qui s'annule avec sa dérivée pour t = la fonction
°
esi
t
Th (t) = - _1_ f-lh
r sin P'h (t-T) th ('t) dT,
Jo
(178) .1
11-3-16.
M~THODE
où comme plus haut
DE FOURIER.
Ill' =
-Â h •
(x, t) = -
~ _1_ \ k-l
513
HYPERBOLIQUES
La somme de la série
t
00
U
~QUATIONS
~h ~
fk
(T) sin Ilh (t -T) dT u k (x)
(179)
satisfait formellement toutes les conditions du problème (174), (175). Pour justifier la formule (179), il faut s'assurer que la série (179) converge comme la série (164) du théorème 1. Prouvons la proposition suivante. Thé 0 r è ID e 2. Supposons que les coefficients de M et le domaine B vérifient les mêmes conditions que dans le théorème 1. Si o f (x, t) E L 2 (D 7), alors la série (179) converge pour la norme de W~ (B) uniformément en t E [0, Tl. Si, de plus, f (x, t) possède une dérivée distributionnelle ft (x, t) EL 2 (DT)' alors la série (179) et lesséries obtenues à partir d'elle par une simple et une double dérivations terme à terme par rapport à x et t convergent pour la norme de L 2 (B) uniformément en t E [0, Tl. Dans le dernier cas la somme u de la série est une solution distributionnelle du problème (174), (175) de la classe (D T) (et même d'une meilleure classe). . Les propositions du théorème résultent des relations et majorations suivantes:
W: p
P
2
~ T k (t) uk (x)
1-
~ Il~T~ (t)~ k=m
k=m
P
t
P
T
~ ~ (~I th (T) 1d't)2 ~T ~ ~ h=m
00
00
(180)
k=m 0
0
et T
n('t) d'tt
T
~ ~ /l('I:) d'l: = ~1 ~ Il ('1:) d'l: =
l
12 (x,
t) dx dt.
(181)
Si de plus ft (x, t) EL 2 (DT)' alors t
'r k (t) =
-
= - !'~
--4\ th ('t) d cos J.tk (t -'t) = ~h ~
[-1 t
1. ('1:) cos P.k (t-'t) d't+ fh (t) -tk (0) cos J.tktJ = t
= -
:1 [)
o
th ('t) (1- cos flk (t - 't» dT
+ th (0)~(1- cos flk t )J.
CH. II. PROBLÈMES AUX LIMITES
5t4
D'où T
T~(t)
(M, t-T) ôx
ôq>(M, t-T) oy
y
x
+
+oq>(M, t-T) ] 0 oz Z ~=t-r = •
(185)
Or, en vertu du théorème d'Euler sur les fonctions homogènes, on a (tome l, [V-1-41) ôq>(M,t-T)(t_ )+oq>(M,t-'t) +oq>(M,t-'t) +oq>(M,t-'t) =0 ot 't OX X oy Y oz z •
On s'assure par la substitution 't = t - r que la relation (185) est réalisée et par suite la formule (184) donne bien une solution de l'équation (183). Cherchons maintenant une solution de l'équation (183) sous la forme spéciale
u-'1'
(+) Y n (e, 1-r,
pour
t~1-r
(200) et 00 ('t) = 0 pour 't ~ O. De cette expression, il s'ensuit comme al' [11-3-17] que la condition (192) est satisfaite pour tout 00 ('t). Il es. immédiat de vérifier que les formules (199) et (200) donnent let solution de l'équation (183) pour n = 0 aussi, pourvu que 00 (T~ possède une dérivée continue. Signalons que la fonction (198) n'e~ pas solution de l'équation (183) pour n = O. La condition aux limites (193) nous conduit à l'équation suivantt t
~
('t) P n ('t --+-1- t) d't = f (t).
00
(201
t-2
Supposons que f (t) admet une dérivée continue et f (0) = t (0) = O. En dérivant l'équation (201) par rapport à t, on trouve 00 (t) +(-1)n+ 1 oo(t-2)t
~
-
00
('t) P~ ('t + 1- t) d't = j' (t)
(n> 1) (2021 ,
t-2
et pour n = 0
/00 (t) -
00
(t -
(202 2)
2) = j' (t).
L'équation (202 1) nous permet de construire 00 ('t) par la méthode des approximations successives. Déterminons d'abord 00 (t) dans l'intervalle [0, 2] à partir de l'équation de Volterra: t
00
(t) -
J ('t) P~ ('t + 1- t) d't= f' (t), 00
o
ensuite
00
(t) dans l'intervalle J2, 4J à partir de l'équation
t
00 (t)-
) oo(t)P~('t+1-t)d't= 2
l 2
= j' (t) + (_1)n 00 (t-2) + co ('t)~P~
('t
+ 1-t) dT .
o
dont le second membre est connu, et ainsi de suite. Portons la fonctior 00 ('t) trouvée dans le second membre de (200).
521
II-3-t8. VIBRATIONS DE L'INT9RIEUR DE LA SPH1l:RE
Utilisons la transformation unilatérale de Laplace pour résoudre l'équation (202 1). Esquissons cette méthode. La formule définitive sera établie plus bas par une autre méthode. Dans l'équation (202 1), représentons l'intégrale par une somme de deux intégrales avec des bornes inférieures nulles, multiplions l~ deux membres par e -st, où S = al 0'2i et al > 0 est un nombre 'Wsez grand et intégrons sur t entre 0 et 00. Posons
+
IH .t~
B
00
00
Q (s)
=
Je-:stw (t) dt,
Je-stf (t) dt;
F (s) =
o
o
("
(203)
;tal
il.ilisons le théorème de convolution (tome IV1 , [1-52]) et les 1:1g,nules
~es formules se déduisent 'sans peine par une intégration direct des premiers membres. Cette méthode nous conduit à l'équation suivante pour Q (s): (~
~
-----
.. /" 2n
V
n
efT (2n+1) e-sJ
s
1 (_
is) Q (s)
= F (s).
n+2'
~
La transformation réciproque de Laplace nous donne co (t) =
- i"::"'( 2n+ 1)
4 ---:::-21t~i-t
O'+ioo
r
J
O'-ioo
1
JI
_ 8
2n
(8) J 1 ( - is) n+e(t+1)SF
ds.
(205)
2
Le nombre réel al est pris assez grand pour que tous les points singuliers de la fonction F (s) soient situés à gauche de la droite d'in tégration. La justification de la possibilité d'appliquer les transformations d,irecte et réciproque de Laplace est facilitée par le fait qu'on a déjà établi l'existence de co (t) par la méthode des approximations suc(~essives et qu'on peut majorer cette fonction pour de grands t en lHt]>osant certaines conditions à f' (t) pour de grands t. En portant j'expression (205) dans (200), en permutant l'ordre d'intégration et 36-01017
522
CH. II. PRDBL:I!:MES AUX LIMITES'
en utilisant la relation facilement démontrable 1
~
J
3t
ePxPn (x) dx=
t ( - ip)
_y--=p=--_,
eiT(2n+t) _n_+_2
V2:rt
(206)
-1
on obtient O'+ioo J
CPn (r, t) =
r
1 _Î
y
T21li
1 ( - irs)
n+-
J.
~
J
(_ is)
3t
F (s) e ds
(n=O, 1, 2, ... ).
n+2
0'-1,00
(207)
Indiquons un moyen plus rapide pour établir la formule (207). En portant (199) dans (183) et en se servant de l'équation pour Y n (e, cp), on obtient l'équation suivante pour q:Jn (r, t): 8 2cpn _ 8t 2 -
8 2 cpn 8r 2
+!
8CPn _
r
n (n+ 1)
8r
r2
q:Jn'
(208)
A cette équaiton il faut adjoindre les conditions q:Jn
It=o =
8~n
CPn Ir=1 =
t=o = 0;
(209)
t (t).
(210)
En multipliant les deux membres de l'équation (208) par e-st , en intégrant sur t entre 0 et 00 et en tenant compte de la condition (209), on obtient pour la fonction 00
X n (r, s)= ~ o
e-stq:Jn
(r, t)dt
(211)
l'équation
d~~n
+; d:; + (-S2- n(~:1) ) Xn=O.
(212)
La transformation de Laplace appliquée à (210) nous donne
= F (s).
X n Ir=1
(213)
La fonction X n (r, 8) doit de plus être finie pour r = O. L'équation (212) se ramène à une équation de Bessel et en se servant de (213) et du fait que X n est finie pour r = 0, on trouve J
X n (r, s.) =
1
y- r
(-irs)
t
n+'2 J
1 ( - is
n+'2
)
F (8).
II-3-19. DISCUSSION DE LA SOLUTION
523
La transformation réciproque de (211) nous conduit ensuite à la formule (207). Les conditions qu'il faut imposer à f (t) pour justifier l'application de la transformation de Laplace et de la formule (207) figurent dans le travail de G. Pet r a che n e, Problèmes dynamiques de la théorie de l'élasticité pour une sphère isotrope. Utch. Zap. LGU, ser. matem. naouk, 1950, nO 21. L'exposé de ce numéro et du suivant est emprunté à ce travail. II-3-19. Discussion de la solution. Etudions la solution
i
o+ioo J lf/l
1 (r, t)=-= , _ .
JI r2m. 0-100
1 (-
J
n+-
.) F (s)
2( 1
irs)
-lS
n+ 2
e~t ds.
(214)
On admettra pour fixer les idées que la fonction f (t) =1= 0 seulement sur un intervalle fini [0, Tl, possède des dérivées première et seconde continues et j (0) = f' (0) = j (T) = f' (T) = ü. Sous ces conditions, une double intégration par parties nous donne T
F(s)=
T
i e-~tf(t)dt= :2 i
e-stj"(t) dt.
(215)
o 0 On admettra que la fonction F (s) est de la forme F (s) = FI (s) F 2 (s) e- sT , (216} où FI (s) et F 2 (s) sont des fractions rationnelles dont le degré du dénominateur est de deux unités au moins supérieur à celui du numérateur. Il est immédiat de vérifier que la fonction F (s) jouira de cette propriété par exemple dans les deux cas suivants:
+
j(t)=t 2 (T-t)2;
j(t)=sin 2
n.:,
tElÜ, T].
On voit sur la formule (215) que F (s) est une fonction entière. La fonction J n +I / 2 (-is) possède des zéros imaginaires que l'on désignera par +ksi (8 = 1, 2, 3, ...), où k s sont les zéros de l'équation J n +I / 2 (k) = (1. (tome 1H z, [V1-2-3]). Utilisons la formule J
1 (z) =
~ [H(1)
n+2 2
+
1 (z) H(2) 1 (z) n+n+2 2
i ]
et les expressions suivantes des fonctions de Hankel:
(217)
524
PROBL~MES
CH. II.
AUX LIMITES
OÙ