Problèmes elliptiques et paraboliques non variationnels
Problèmes elliptiques et paraboliques non variationnels
9 janvier 1997
UMPA-ENS Lyon
46 Allée d’Italie
69364 Lyon Cedex 07
France
Table des Matières
Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2 Outils de base . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.1 Résultats d’intégration . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.2 Espaces de Sobolev et analyse fonctionnelle . . . . . . . . . . . . . . . 20
10 Espaces de Lorentz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
10.1 Espaces de Lorentz . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
10.2 Quelques fonctions . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
Bibliographie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 155
Introduction
Les équations aux dérivées partielles de types elliptiques et paraboliques ont été
étudiées depuis longtemps et le cas de problèmes linéaires à données H −1 (elliptique)
ou L2 (0, T ; H −1 ) (parabolique), cadre (( variationnel )), est bien connu. Il y a existence
et unicité dans les (( bons espaces )) et les équations sont formulées de la façon suivante
(où Ω est un ouvert borné de RN ) : pour une équation elliptique avec conditions aux
limites de Dirichlet et f ∈ H −1 (Ω), on cherche (Lax-Milgram [32]) u vérifiant
1
u ∈ H0 (Ω)
Z
(A∇u)∇ϕ = hf, ϕi ∀ϕ ∈ H01 (Ω)
Ω
où A est une matrice bornée et coercitive, et pour une équation parabolique avec
u0 ∈ L2 (Ω) et f ∈ L2 (0, T ; H −1 (Ω)), on cherche (Lions-Magenes [35]) u vérifiant
Les hypothèses sur u entraı̂nent que u ∈ C(0, T ; L2 (Ω)), la condition initiale a donc
bien un sens.
Si l’on considère le problème elliptique dans le cas où les données ne sont plus
H −1 , cette formulation n’est plus adaptée. Pour f ∈ M (Ω) = (C(Ω))0 , espace de
mesures, Stampacchia [48] a proposé en 1965 une méthode donnant dans le cas linéaire
existence et unicité des solutions d’une équation elliptique (nous la précisons dans
[43]). Cette méthode utilise la dualité (et un résultat de régularité) et conduit à la
formulation suivante
\
W01,q (Ω)
u ∈
q < NN
−1
(1)
Z
t
− hdiv(A ∇ϕ), ui = ϕ df
Ω
∀ϕ ∈ H01 (Ω) ∩ C(Ω) tel que div(At ∇ϕ) ∈ r > N W −1,r (Ω).
S
8 Introduction
Comme W01,r (Ω) ⊂ C(Ω), et ϕ ∈ W01,r (Ω) implique div(At ∇ϕ) ∈ W −1,r (Ω), les
fonctions test de (2) peuvent aussi être choisies commes fonctions test pour (1), mais
la réciproque est fausse. Aussi (2) n’assure pas l’unicité des solutions comme nous
allons le voir, alors que (1) l’assure.
0
Pour f ∈ W −1,p (Ω) et p > 1, c’est-a-dire, ici encore, dans le cadre (( variationnel )),
les équations non-linéaires de la forme
− div(a(x, u, ∇u)) = f
où a(x, s, ξ) défini un opérateur elliptique ont été étudiées par Leray et Lions (1965)
(voir [33]), qui ont montré qu’il y avait existence d’une solution u vérifiant
1,p
u ∈ W0 (Ω)
(3)
Z
a(x, u, ∇u)∇ϕ = hf, ϕi ∀ϕ ∈ W01,p (Ω)
Ω
0 0
et pour le cas parabolique (Lions [34]) avec u0 ∈ L2 (Ω) et f ∈ Lp (0, T ; W −1,p (Ω)) il
existe u vérifiant
0 0
u ∈ Lp (0, T ; W01,p (Ω)), ut ∈ Lp (0, T ; W −1,p (Ω))
Z T Z TZ Z T
hut , ϕi + a(t, x, u, ∇u)∇ϕ = hf, ϕi ∀ϕ ∈ Lp (0, T ; W01,p (Ω))
0 0 Ω 0
u(0) = u0
[41]). La formulation que vérifie la solution est proche de celle des distributions : pour
l’équation elliptique dans le cas de Dirichlet homogène, et f ∈ M (Ω)
\
u ∈ W01,q (Ω)
q < NN
(p−1)
−1
Z Z [
W01,r (Ω)
a(x, u, ∇u)∇ϕ = ϕ df ∀ϕ ∈
Ω Ω r>N
Dans le cas où a ne dépend pas de u, l’existence d’une solution a d’abord été
donnée dans [11] (et plus rapidement dans [15]) avec une hypothèse technique sur a
afin d’avoir la convergence presque partout de ∇un vers ∇u. Cette condition a été
levée dans [12], et le résultat a été montré pour a(x, u, ∇u) dans le cas elliptique (voir
[26, 25, 41]). Dans le cas parabolique avec a(t, x, ∇u), ceci a été fait par Dall’Aglio
et Orsina [24], puis pour a(t, x, u, ∇u) dans [10].
Nous avons montré, dans[43], que dans le cas elliptique linéaire la méthode par
approximation utilisée pour montrer l’existence conduit à la même solution que celle
de Stampacchia. Cependant la formulation (2) est trop faible pour assurer l’unicité
comme nous le prouvons (voir [43]) avec le contre-exemple de Serrin [46], qui permet
d’obtenir, dans le cas linéaire, une solution non nulle pour un problème à données
nulles.
Pour remédier à cela, dans le cas elliptique avec conditions aux limites homogènes
de Dirichlet, où les données sont L1 (Ω) (ce qui est bien sûr plus restrictif que M (Ω)),
trois approches ont été utilisées.
Dall’Aglio [23] a montré que, même pour un problème non-linéaire, la méthode
par approximation conduit à une unique solution appelée SOLA (Solution Obtenue
comme Limite d’Approximations).
Bénilan, Boccardo, Gallouët, Gariepy, Pierre et Vazquez [4] définissent la notion
de solution entropique (ici p > 2 − 1/N pour simplifier) :
1,1
Tk (u) ∈ W01,p (Ω) ∀k > 0
u ∈ W (Ω),
Z Z
a(x, ∇u)∇Tk (u − ϕ) ≤ f Tk (u − ϕ) ∀ϕ ∈ W01,p (Ω) ∩ L∞ (Ω).
Ω Ω
et l’unicité de solutions entropiques a été montrée pour des conditions aux limites
homogènes plus générales et en particulier de Neumann et Fourier par Andreu, Mazón,
Segura de León et Toledo [3].
Enfin, Lions et Murat [36, 39] ont introduit une autre notion : celle de solution
renormalisée (suite aux solutions du même nom, dues à Di Perna et Lions, pour
l’équation de Boltzmann). Elle vérifie (ici aussi, p > 2 − 1/N pour simplifier) :
u ∈ W 1,1 (Ω), Tk (u) ∈ W01,p (Ω) ∀k > 0
Z
lim |∇u|p = 0 ∀k > 0
h→+∞
h≤|u|≤k+h
Z Z Z
0
S(u) a(x, ∇u)∇ϕ + S (u) ϕ a(x, ∇u)∇u = f S(u) ϕ
Ω Ω Ω
∀ϕ ∈ W01,p (Ω) ∩ L∞ (Ω)
Dans le cas parabolique et p > 2−1/(N +1), nous avons, de façon analogue, défini
la solution entropique d’une équation parabolique (voir [44]) :
u ∈ C(0, T ; L1 (Ω)) ∩ L1 (0, T ; W 1,1 (Ω)), Tk (u) ∈ Lp (0, T ; W01,p (Ω)) ∀k > 0,
Z Z Z T
Θk (u − ϕ)(T ) − Θk (u0 − ϕ(0)) − hϕt , Tk (u − ϕ)i
Ω
Ω
Z TZ
0
Z
+ a(t, x, ∇u)∇Tk (u − ϕ) ≤ f Tk (u − ϕ)
0 Ω Ω
∀ϕ ∈ Lp (0, T ; W01,p (Ω)) ∩ L∞ (]0, T [×Ω) ∩ C(0, T ; L1 (Ω))
0 0
tel que ϕt ∈ Lp (0, T ; W −1,p (Ω)).
Blanchard et Murat [8] ont aussi défini une solution renormalisée qui vérifie
u ∈ C(0, T ; L1 (Ω)) ∩ L1 (0, T ; W 1,1 (Ω)), Tk (u) ∈ Lp (0, T ; W01,p (Ω)),
Z TZ
|∇u|p = 0 ∀k > 0,
lim
h→+∞ 0 h≤|u|≤k+h
u(0) = u0 ,
Z T Z TZ
S 0 (u) a(t, x, ∇u)∇ϕ
− hS(u), ϕt i +
0 0 Ω
Z TZ Z TZ
00
f S 0 (u) ϕ
+ S (u) ϕ a(t, x, ∇u)∇u = ∀ϕ ∈ D(Ω)
0 Ω 0 Ω
Dans le cas où le second membre est une mesure quelconque, on ne sait plus, en
général, démontrer que la notion de SOLA assure l’unicité, et les solutions entropiques
et renormalisées ne sont plus définies. Pour les équations elliptiques, de nouvelles
définitions équivalentes ont été proposées par Dal Maso, Murat, Orsina et Prignet
[22] : elles généralisent les trois notions précédentes, mais ne conduisent pas (encore)
à l’unicité.
Une mesure f ∈ M (Ω) peut être décomposée en f = µ0 + λ+ − λ− où µ0 est une
0
mesure ne chargeant pas les ensembles de p-capacités nulles (µ0 ∈ L1 (Ω)+W −1,p (Ω))
et λ+ et λ− les parties positives et negatives de f − µ0 . Voici une des formulations
vérifiée par u :
Z Z Z Z
a(x, ∇u)∇w dx = w dµ0 + w dλ − w−∞ dλ−
+∞ +
Ω Ω Ω Ω
pour tout w ∈ W01,p (Ω) ∩ L∞ (Ω) tel qu’il existe k > 0, w+∞ et w−∞ ∈ W 1,r (Ω), avec
r > N , tels que
(
w = w+∞ capp -quasi partout sur {u > k},
w = w−∞ capp -quasi partout sur {u < −k}.
Cet ouvrage est constitué de chapitres pouvant être lus indépendamment les uns
des autres. Les chapitres 1 et 2 sont des rappels concernant les équations elliptiques
et paraboliques dans un cadre (( classique )), l’intégration et les espaces de Sobolev.
La première partie est consacrée aux méthodes de dualité pour des problèmes
elliptiques linéaires : le chapitre 3 concerne les équations à seconds membres W −1,s (Ω)
pour s < 2, et le chapitre 4 les équations à seconds membres mesure.
La deuxième partie fait intervenir des équations elliptiques dans une formulation
au sens des distributions. Le chapitre 5 donne l’existence pour des équations non-
linéaires, tandis que le chapitre 6 prouve que même dans le cas linéaire il n’y a pas
unicité.
La troisième partie est consacrée aux équations à seconds membres L1 : le chapitre
7 montre la continuité des solutions entropiques d’une équation elliptique par rapport
à un paramètre (et redémontre leur unicité). Et le chapitre 8 prouve l’existence et
l’unicité des solutions entropiques d’équations paraboliques.
12 Introduction
Enfin la quatrième partie donne des propriétés des solutions des équations
elliptiques à seconds membres mesure (dont on connaı̂t déjà l’existence et la non-
unicité). Dans le chapitre 9, on étudie les singularités de la solution pour un second
membre somme d’un dirac et d’une fonction L1 . Puis le chapitre 10 montre que les
espaces de Lorentz, bien que plus fins que les espaces de Sobolev, ne permettent
pas de distinguer les solutions à second membre dirac de celles à second membre
L1 . Le chapitre 11 est consacré à l’étude d’une formulation plus générale que les
formulations entropiques et renormalisées, mais ne donnant pas l’unicité lorsque le
second membre est un dirac. Lorsque le second membre est L1 , il y a convergence
forte des troncatures des approximations, nous montrons dans le chapitre 12 que ce
n’est plus nécessairement le cas pour une mesure quelconque. Enfin, le chapitre 13
propose une généralisation des notions de SOLA, solutions entropiques, et solutions
renormalisées pour une mesure quelconque au second membre.
Les chapitres 4 et 6 sont les deux parties (en français) de l’article [43] publié à
Rendiconti di Matematica.
Le chapitre 5 est à paraı̂tre aux Annales de Toulouse [41].
Le chapitre 8 est paru (en anglais) à Nonlinear Analysis TMA [44].
Le chapitre 7 est un article en préparation [42].
Le chapitre 13 est un article en préparation [22] écrit en collaboration avec G. Dal
Maso, F. Murat et L. Orsina, une note CRAS est aussi en préparation.
1. Cadre classique pour les équations elliptiques et
paraboliques
2
P
où ∆ = i=1,...,N ∂xi (par la suite, les sommes en i ne seront plus notées). Les
problèmes modèles pour les équations non linéaires font intervenir le p-laplacien
∆p (u) = div(|∇u|p−2 ∇u)
− div(A∇u) = f dans Ω
u = 0 sur ∂Ω
où A est une matrice de taille N . Si l’on note aij les coefficients de A, l’équation
s’écrit −∂xi (aij ∂xj u) = f. Nous supposons que aij ∈ L∞ (Ω) et satisfait la condition
d’ellipticité (ou coercitivité) : il existe α > 0 tel que
X X
∀(ξi ) ∈ RN ξi aij ξj ≥ α ξi ξj .
i,j i,j
14 Chap. 1 Cadre classique pour les équations . . .
et
u(0) = u0 dans L2 (Ω),
qui a bien un sens car u ∈ L2 (0, T ; H01 (Ω)) et ut ∈ L2 (0, T ; H −1 (Ω)) implique
u ∈ C([0, T ], L2 (Ω)).
La solution est obtenue par une méthode de Galerkin, c’est-à-dire comme limite de
solutions appartenant à des espaces de dimensions finies.
Nous supposons que a satisfait les hypothèses suivantes (qui font de A(u) un
opérateur de Leray-Lions) : a est une fonction de Carathéodory, c’est-à-dire que
• a(x, s, ξ) : RN × R × RN → RN est mesurable en x ∈ RN pour tout s ∈ R
et ξ ∈ RN et continue en ξ ∈ RN et s ∈ R pour presque tout x ∈ RN . Nous
noterons a(x, u, ∇u) = a(x, u(x), ∇u(x)).
et a vérifie aussi des conditions de coercitivité, monotonie stricte et croissance : il
existe p vérifiant 1 < p ≤ N
• Il existe α > 0 tel que pour tout s et ξ et presque tout x on ait
Ces hypothèses sont classiques pour l’étude des opérateurs non linéaires sous forme
divergentielle.
0
avec u0 dans L2 (Ω) et f dans Lp (]0, T [; W −1,p (Ω)) où a est une fonction de
Carathéodory satisfaisant des conditions de coercitivité, de monotonie stricte et
de croissance du type de celles de Leray-Lions, définissant un opérateur sur
Lp (]0, T [; W01,p (Ω)).
La solution u, obtenue par Lions [34] (ici aussi grâce à une méthode de Galerkin),
vérifie
0 0
u ∈ Lp (]0, T [; W01,p (Ω)) et ut ∈ Lp (]0, T [; W −1,p (Ω))
Z T Z TZ Z T
p 1,p
∀ϕ ∈ L (]0, T [; W0 (Ω)) hut , ϕi + a(t, x, u, ∇u)∇ϕ = hf, ϕi
0 0 Ω 0
u(0) = u0 ,
ici aussi cette dernière égalité a bien un sens car u ∈ C([0, T ]; L2 (Ω)), ce qui est une
conséquence de la régularité de u et ut .
Donnons les hypothèses sur a, qui font de A(u) = − div(a(t, x, u, ∇u)) un
opérateur de Leray-Lions : a est une fonction de Carathéodory, c’est-à-dire que
• a(t, x, s, ξ) : R × RN × RN → RN est mesurable en t ∈ R et x ∈ RN pour
tout s ∈ R et ξ ∈ RN et continue en s ∈ R et ξ ∈ RN pour presque tout
t ∈ R et x ∈ RN . Nous noterons a(t, x, u, ∇u) = a(t, x, u(t, x), ∇u(t, x)).
et a vérifie aussi des conditions de coercitivité, monotonie stricte et croissance : il
existe p vérifiant 1 < p ≤ N et
• Il existe α > 0 tel que pour tout s et ξ et presque tout t et x on ait
a(t, x, s, ξ)ξ ≥ α|ξ|p
• Pour tout s, ξ et η et presque tout t et x on a
[a(t, x, s, ξ) − a(t, x, s, η)](ξ − η) > 0 pour ξ 6= η
0
• Il existe b(t, x) ∈ Lp (]0, T [×Ω), (où p0 = p/(p − 1)) et β > 0 tels que pour
tout s et ξ et presque tout t et x on ait
|a(t, x, s, ξ)| ≤ β(b(t, x) + |s|p−1 + |ξ|p−1 )
Ces hypothèses sont classiques pour l’étude des opérateurs non linéaires sous forme
divergentielle.
1.4 Troncature
et sa primitive Θk : R → R+
Z x
Θk (x) = Tk (s) ds.
0
1.5 Mesure 17
La fonction Tk sert, par exemple, à montrer le principe du maximum faible pour les
équations elliptiques, elle va être très utilisée dans la suite.
Comme cette fonction est lipschitzienne, le théorème de Stampacchia affirme que
pour u une fonction W 1,p (Ω), avec p ≥ 1, on a Tk (u) ∈ W 1,p (Ω) et
1.5 Mesure
Nous reprenons les énoncés de Kavian [31], les résultats n’y sont pas démontrés
(sauf le théorème de Vitali), les démonstrations peuvent être trouvées dans Rudin
[45]. Dans la suite Ω sera un ouvert borné de RN et nous utiliserons la mesure de
Lesbegue, mais ces résultats restent vrais dans un cadre plus général.
Lemme de Fatou : Soit (fn ) une suite de fonctions mesurables positives. Alors
Z Z
lim fn dx ≤ lim fn dx.
Ω n→∞ n→∞ Ω
Conséquence : Si (un ) est bornée dans W 1,p (Ω) avec p > 1 (et donc convergeant
dans W 1,p (Ω) faible à une sous suite près) et si (∇un ) converge presque partout alors
(un ) converge (à une sous suite près) dans W 1,q (Ω) (fort) ∀q < p.
En effet, d’après le théorème de Vitali, il suffit de montrer que (|∇un |q ) est équi-
intégrable or pour E borelien
Z Z q/p Z 1−q/p
q p
|∇un | ≤ |∇un | 1 ≤ C mes(E)1−q/p .
E Ω E
Lemme : Soit f une fonction mesurable strictement positive. Alors pour tout
ε > 0 il existe δ > 0 tel que pour tout A ⊂ Ω mesurable on ait
Z
f dx < δ =⇒ mes(A) < ε.
A
Lemme : Soit (fn ) une suite de fonctions mesurables de Ω dans R. Si (fn ) est
de Cauchy en mesure, c’est-à-dire si
∀ε > 0, ∀η > 0, ∃n0 ∈ N tel que p, q ≥ n0 ⇒ mes({x ∈ Ω; |fp (x) − fq (x)| > ε}) ≤ η,
Notons Ak = {x ∈ Ω; |fnk+1 (x) − fnk (x)| > 1/2k }, et soit Bp = k≥p Ak alors
S
mes(Bp ) ≤ 1/2p−1 . Sur Bpc , son complémentaire, la suite (fnk ) est de Cauchy et
converge donc uniformément et simplement vers gp . Pour p0 ≥ p on a Bpc ⊂ Bpc0 donc
gp0 /Bpc = gp si bien qu’il existe g tel que g/Bpc = gp . Ainsi (fnk ) converge simplement
vers g sur p∈N Bpc . Et ( p∈N Bpc )c = p∈N Bp qui est de mesure nulle, donc (fnk )
S S T
converge presque partout vers g.
Nous reprenons ici aussi certains énoncés de Kavian [31] et de Brezis [19], pour
une présentation plus complète des espaces de Sobolev, on pourra voir Adams [1].
Par la suite, Ω est un ouvert borné de RN .
2.2 Espaces de Sobolev et analyse fonctionnelle 21
et pour p = ∞, on note
de norme
kukL∞ = sup |u|.
Ω
et pour p = ∞
kukW 1,∞ = max(kukL∞ , k∇ukL∞ ).
L’espace D(Ω) désignant l’ensemble des fonction de classe C ∞ (Ω) à support compact
dans Ω, on note pour 1 ≤ p < ∞
W 1,p (Ω)
W01,p (Ω) = D(Ω) .
Comme pour 1 ≤ p < ∞, l’espace D(Ω) est par définition dense dans W01,p (Ω), on peut
identifier le dual de W01,p (Ω) à un sous-espace de l’espace des distributions D0 (Ω) (le
0
dual de Lp (Ω) étant identifié à Lp (Ω) où p0 = p/(p − 1)). On note :
0
W −1,p (Ω) = (W01,p (Ω))0 .
Dans le cas particulier où p = 2, W 1,2 (Ω), W01,2 (Ω) et W −1,2 (Ω) sont respectivement
notés H 1 (Ω), H01 (Ω) et H −1 (Ω).
Théorème : Soit E un espace de Banach séparable et soit (fn ) une suite bornée
dans E 0 . Alors il existe une sous suite (fnk ) qui converge pour la topologie faible ∗ de
E 0.
Ce théorème est souvent utilisé pour E = L1 (Ω), qui est bien séparable (mais pas
réflexif), comme E 0 = L∞ (Ω), on peut donc extraire d’une suite bornée dans L∞ (Ω),
une sous suite qui converge dans L∞ (Ω) faible ∗.
22 Chap. 2 Outils de base
Théorème : Soit E un espace de Banach réflexif et soit (fn ) une suite bornée
dans E. Alors il existe une sous suite (fnk ) qui converge pour la topologie faible de
E.
Ce théorème permet d’extraire des sous suites convergentes faiblement de suites
bornées dans Lp (Ω) pour 1 < p < ∞, puisque Lp (Ω) est alors réflexif. Il ne s’applique
donc ni à L1 (Ω), ni à L∞ (Ω) (pour lequel, en revanche, le théorème précédent
s’applique). Il en est de même pour les W01,p (Ω).
kukLp ≤ Ck∇ukLp .
23
3. Relation entre regularité et existence du
problème adjoint
On s’intéresse ici à des équations elliptiques linéaires à seconds membres W −1,s avec
s < 2, dont on montre l’existence de solutions par dualité.
Si f est plus régulière, la solution peut être plus régulière, par exemple grâce au
théorème de Meyers [6]. Aussi introduisons nous
r0 = sup{r1 ; f ∈ W −1,r (Ω) ⇒ u ∈ W01,r (Ω) ∀ 2 ≤ r ≤ r1 }.
Par construction, r0 ≥ 2 et le théorème de Meyers dit, en fait, que r0 > 2. Ceci par
dualité va permettre de résoudre (3.1) pour des f moins réguliers, en particulier pour
f ∈ W −1,q (Ω) avec q < 2 : aussi de façon analogue à la définition de r0 définit-on
q0 = inf{q1 ; pour f ∈ W −1,q (Ω),
(3.1) admet une unique solution dans W01,q (Ω) ∀q1 ≤ q ≤ 2}.
Théorème : Régularité des solutions et existence d’une solution sont duales, i.e.
q0 = r 0 0 .
26 Chap. 3 Relation entre regularité et existence du problème adjoint
Démonstration : Soit 1 < s < ∞ tel que ∀ g ∈ W −1,s (Ω) il existe un unique
v ∈ W01,s (Ω) qui vérifie
Z
0
(A∗ ∇v)∇ϕ = hg, ϕiW −1,s ,W 1,s0 ∀ ϕ ∈ W01,s (Ω).
0
Ω
Notons Gs : W −1,s (Ω) → W01,s (Ω) tel que Gs (g) = v, on notera que Gs est bijectif et
Z
0
(3.2) (A∗ ∇Gr (g))∇ϕ = hg, ϕiW −1,s ,W 1,s0 ∀ ϕ ∈ W01,s (Ω),
0
Ω
De plus, − div(A∗ ∇) est continue de W01,s (Ω) dans W −1,s (Ω) et bijective d’inverse Gs
aussi grâce au théorème de l’application ouverte, Gs est continue de W −1,s (Ω) dans
W01,s (Ω), à l’aide du théorème de Riesz, on peut donc définir Gs ∗ son adjoint, tel que
0 0 0
Gs ∗ : W −1,s (Ω) → W01,s (Ω) et par définition pour f ∈ W −1,s (Ω) on a u = Gs ∗ (f ) si
(3.3) ∀ g ∈ W −1,s (Ω) hu, giW 1,s0 ,W −1,s = hf, Gs (g)iW −1,s0 ,W 1,s
0
on a avec (3.3)
Z
∀v ∈ W01,s (Ω) (A∇v)∇u = hf, viW −1,s0 ,W 1,s ,
Ω
0 0
pour f ∈ W −1,s (Ω), ainsi a-t-on trouvé une solution de (3.1) dans W01,s (Ω), et elle
est unique : en effet, il suffit de montrer que pour f = 0 on a u = 0. Or grâce à la
formule de Green
∀ g ∈ W −1,s (Ω) hu, giW 1,s0 ,W −1,s = 0
0
1,s0 0
et g décrit W −1,s (Ω) le dual de W0 (Ω) donc u = 0. Donc pour f ∈ W −1,s (Ω), (3.1)
0
admet une unique solution dans W01,s (Ω) si et seulement si (il suffit d’échanger les
rôles de A et A∗ ) (3.1) écrit avec A∗ admet pour f ∈ W −1,s (Ω) une unique solution
dans W01,s (Ω) donc q0 = r0 0 .
3.2 r0 peut être fini et même très proche de 2 27
Considérons r > r0 tel qu’il existe g ∈ W −1,r (Ω) n’admettant pas de solutions
dans W01,r (Ω), donc
− div(A∗ ∇W01,r (Ω)) 6= W −1,r (Ω).
Or si E ⊂ F ,
E = F ⇐⇒ E fermé et E dense dans F
donc ici, on a soit − div(A∗ ∇W01,r (Ω)) non fermé, soit − div(A∗ ∇W01,r (Ω)) non dense.
Théorème : − div(A∗ ∇W01,r (Ω)) fermé est équivalent à la propriété : (3.1) admet
0
une solution pour tout f ∈ W −1,r (Ω).
Démonstration : D’après la théorie des opérateurs (voir Brezis [19], par exemple),
0
− div(A∗ ∇W01,r (Ω)) fermé équivaut à − div(A∇W01,r (Ω)) fermé, or si f ∈ H −1 (Ω),
0
il existe une solution de (3.1) dans H01 (Ω) et donc dans W01,r (Ω), donc H −1 (Ω) ⊂
0 0 0
− div(A∇W01,r (Ω)) et H −1 (Ω) est dense dans W −1,r (Ω) donc − div(A∇W01,r (Ω))
0 0
aussi. Ainsi − div(A∇W01,r (Ω)) est fermé et dense dans W −1,r (Ω), donc on a
0 0
− div(A∇W01,r (Ω)) = W −1,r (Ω) et il y a donc existence d’une solution de
0
(3.1). Réciproquement si (3.1) admet une solution pour tout f ∈ W01,r (Ω), alors
0
− div(A∇W01,r (Ω)) est fermé, donc − div(A∗ ∇W01,r (Ω)) aussi.
Théorème : − div(A∗ ∇W01,r (Ω)) dense est équivalent à la propriété : (3.1) admet
0
au plus une solution pour tout f ∈ W −1,r (Ω).
et hv, gi = 1, ainsi v est non nul et vérifie, grâce à la formule de Green, (3.1) avec
0
f = 0, et v peut être prolongé grâce à Hahn-Banach à W −1,r (Ω) donc v ∈ W −1,r (Ω).
Comme 0 est aussi solution, il n’y a pas unicité.
Remarque : Le contre exemple de Serrin montre qu’il existe des équations pour
lesquelles il n’y a pas unicité, et donc telles que − div(A∗ ∇W01,r (Ω)) n’est pas dense.
Il existe aussi des A∗ tels que − div(A∗ ∇W01,r (Ω)) ne soit pas fermé.
4. Solution par dualité de Stampacchia
On montre ici l’existence et l’unicité des solutions pour des équations elliptiques
linéraires à seconds membres mesure, à l’aide d’un argument de dualité. Nous
précisons la méthode de Stampacchia [48].
pour α > 0 et f ∈ M (Ω) = (C(Ω))0 espace de mesures qui est le dual de l’espace des
fonctions continues sur Ω muni de sa norme habituelle.
Ce problème admet, pour f ∈ H −1 (Ω), une unique solution variationnelle, qui
appartient à H01 (Ω) et vérifie
Z
1
∀ v ∈ H0 (Ω) (A∇u)∇v dx = hf, viH −1 ,H01
Ω
Pour Tf 6∈ H −1 (Ω), mais f ∈ M (Ω) on ne trouve plus de solutions dans H01 (Ω), mais
dans q < N W01,q (Ω), ce qui conduit à une formulation plus faible, puisque pour que
N −1
la première intégrale ait un sens il faut alors que v ∈ p > N W 1,p (Ω) :
S
[ Z Z
(4.2) ∀v ∈ W01,p (Ω) (A∇u)∇v dx = v df.
p>N Ω Ω
Nous avons W01,p (Ω) ⊂ C(Ω) pour p > N aussi le membre de droite a un sens.
30 Chap. 4 Solution par dualité de Stampacchia
Cette formulation étant plus faible que la formulation variationnelle elle n’assure
plus l’unicité, comme le montre le contre-exemple de Serrin [46] que nous présentons
au chapitre 6.
L’existence de solutions vérifiant cette formulation a été obtenue de plusieurs
façons, nous allons nous intéresser à deux types de solutions : les solutions obtenues
par dualité et les solutions obtenues par approximation. Le premier est dû à Stam-
pacchia [48], les solutions vérifient une formulation plus forte qui assure l’unicité mais
n’est applicable qu’à un problème linéaire et le second à Boccardo et Gallouët [11],
il est applicable à un problème non linéaire mais n’assure pas l’unicité, aussi pour
préciser la formulation, des inégalités d’entropie ont été introduites [4].
Pour p1 > p2 , on a W −1,p1 (Ω) ⊂ W −1,p2 (Ω) etSon peut vérifier aisément que
Gp2 /W −1,p1 (Ω) = Gp1 , on peut donc définir G de p > N W −1,p (Ω) dans C(Ω) par
G/W −1,p (Ω) = Gp et G∗ son adjoint, opérateur de M (Ω) dans q < N W01,q (Ω), par
T
N −1
∀p>N ∀ g ∈ W −1,p (Ω) hG∗ µ, giW 1,p0 ,W −1,p = hµ, GgiM (Ω),C(Ω)
0
A g fixé l’application µ → hµ, GgiM (Ω),C(Ω) est continue de M (Ω), muni de la topologie
faible ∗, dans R donc G∗ est continu au sens faible suivant
Soit p > N et g ∈ W −1,p (Ω) alors µ → hG∗ µ, giW 1,p0 ,W −1,p est continue de
0
M (Ω), muni de la topologie faible ∗, dans R
Nous pouvons aussi utiliser le cadre des espaces localement convexes (voir Conway
[21]). W −1,p (Ω) est un espace normé donc localement convexe et si p1 > p2 nous
4.2 La théorie de Stampacchia 31
avons W −1,p1 (Ω) ⊂ W −1,p2 (Ω) avec inclusion continue, nous pouvons donc munir
−1,p
(Ω) de la topologie limite inductive de celles des W −1,p (Ω). Pour cette
S
p>N W
topologie, un opérateur est continu sur p > N W −1,p (Ω) si et seulement s’il est continu
S
sur chaque W −1,p (Ω) pour p > N , donc G est continu de p > N W −1,p (Ω) dans C(Ω) et
S
( p > N W −1,p (Ω))0 = q < N W01,q (Ω). Nous pouvons alors munir q < N W01,q (Ω)
S T T
N −1 N −1
de la topologie faible ∗ et on vérifie aisément que la continuité faible de G∗ définie
précédemment est équivalente à celle de G∗ de M (Ω) dans q < N W01,q (Ω) pour les
T
N −1
topologies faibles ∗.
Comme g ∈ p > N W −1,p (Ω), Gg ∈ C(Ω), f ∈ M (Ω) et u ∈ q < N W01,q (Ω) (car Ω
S T
N −1
borné et u ∈ H01 (Ω)), nous obtenons
∀p>N ∀ g ∈ W −1,p (Ω) hu, giW 1,p0 ,W −1,p = hf, GgiM (Ω),C(Ω)
0
donc pour f ∈ H −1 (Ω)∩M (Ω) nous avons u = G∗ f , ce qui nous conduit à la définition
suivante
le dual de W01,q (Ω) (avec 1/q + 1/q 0 = 1) donc u est unique. Cette formulation assure
donc bien l’unicité.
1
u ∈ L (Ω)
(4.7)
Z Z
∗ ∗
∀v ∈ H01 (Ω) ∩ C(Ω) tel que L v ∈ C(Ω) u L v dx = v df.
Ω Ω
S
−1,p
Cependant G p>N W (Ω) n’est pas connu, ce qui rend peu claire cette formu-
lation, aussi allons nous montrer quelques inclusions.
Nous avons vu que Gg est, par définition, la solution
S variationnelle de L∗ v = g
1 −1,p
d’où Gg ∈ H0 (Ω) et que de plus Gg ∈ C(Ω) donc G( p > N W (Ω)) ⊂ H01 (Ω) ∩
C(Ω). Soit v ∈ p > N W01,p (Ω) alors L∗ v ∈ p > N W −1,p (Ω) et puisque p > N ,
S S
v ∈ H01 (Ω) donc v = G(L∗ v) et donc p > N W01,p (Ω) ⊂ G( p > N W −1,p (Ω)), ainsi
S S
avons nous
[ S
W01,p (Ω) ⊂ G p>N W −1,p
(Ω) ⊂ H01 (Ω) ∩ C(Ω),
p>N
donc
S
−1,p
= {v ∈ H01 (Ω) ∩ C(Ω) tel que L∗ v ∈ p > N W −1,p (Ω)}
S
G p>N W (Ω)
Avant de montrer l’équivalence entre (4.5) et (4.7), notons que (4.7) a été
introduite par Stampacchia pour avoir une formulation simple c’est-à-dire ne faisant
intervenir que des intégrales.
Soit g ∈ C(Ω) ⊂ p > N W −1,p (Ω) alors il existe v ∈ H01 (Ω)∩C(Ω) tel que g = L∗ v
S
la formulation (4.7) peut donc s’écrire, pour f = 0,
Z
∀g ∈ C(Ω) u g dx = 0
Ω
1
ce qui implique u = 0 puisque u ∈ L (Ω). Ainsi (4.7) assure l’unicité de u.
Pour montrer l’équivalence de (4.5) et (4.7) il suffit donc de montrer que
la
T solution1,qde (4.5) est solution de (4.7). La solution u de (4.5) vérifie u ∈
1
q< N W 0 (Ω) donc u ∈ L (Ω) et
N −1
La solution définie par Boccardo et Gallouët, dans [11], est la limite de solutions
obtenues pour des fonctions f plus régulières. Soit f ∈ M (Ω) et soit hfn i ∈ H −1 (Ω) ∩
L1 (Ω) une suite convergeant vers f au sens des distributions, avec kfn kL1 (Ω) ≤
kf kM (Ω) ; cette suite existe par densité de H −1 (Ω) ∩ L1 (Ω) dans M (Ω) pour la
topologie faible ∗ de M (Ω). Soit un la solution variationnelle de Lun = fn , alors
kun kW 1,q (Ω) < C, où C ne dépend que de q, L, Ω et kf kM (Ω) , pour tout q tel que
0
1 ≤ q < NN−1 (voir [11]).
Il existe alors un u ∈ W01,q (Ω) et une sous-suite notée, à nouveau, (un ) tels que
un * u dans W01,q (Ω) faible, donc Lun → Lu au sens des distributions, et donc
Lu = f dans D0 (Ω), c’est-à-dire que
Z Z
∀ ϕ ∈ D(Ω) (A∇u)∇ϕ dx = ϕ df
Ω Ω
34 Chap. 4 Solution par dualité de Stampacchia
et donc par densité de D(Ω) dans p > N W01,p (Ω), ( Ω v df a un sens puisque p > N
S R
Comme p > N W01,p (Ω) ⊂ G( p > N W −1,p (Ω)) (et l’inclusion peut être stricte), (4.10)
S S
0
est donc plus faible que (4.8). Et pour tout q fixé, L∗ v ne décrit plus W −1,q (Ω), le
0
dual de W01,q (Ω), pour q < NN−1 , mais un de ses sous espaces L∗ (W01,q (Ω)), l’unicité
donc n’est plus assurée.
Ainsi la formulation au sens des distributions est plus faible que celle de Stam-
pacchia et n’assure pas l’unicité (voir le chapitre 6) ; et la solution de Stampacchia
vérifie (4.8) et donc vérifie aussi (4.9).
soit, avec la formule de Green, hun , L viH01 ,H −1 = hfn , viH01 ,H −1 , or G( p > N W −1,p (Ω))
∗
S
Remarque : si aij ∈ W 1,∞ (Ω) alors un lemme de Meyers dit que, pour
g S∈ W −1,p (Ω), la solution de L∗ v = g appartient à W01,p (Ω), dans ce cas
G( p > N W −1,p (Ω)) = p > N W01,p (Ω) et l’équation Lu = f admet une unique so-
S
lution vérifiant (4.9), ce qui est plus fort que l’unicité des solutions de Boccardo-
Gallouët ou de Stampacchia. Ceci est classique, il suffit en fait que aij ∈ C 0,α (Ω).
Part II
35
5. Existence de solutions de problèmes elliptiques
avec second membre mesure et conditions aux
limites non homogènes
Donnons les hypothèses sur a, opérateur de Leray-Lions, nous les appellerons (H)
par la suite : a est une fonction de Carathéodory, c’est-à-dire que
• a(x, s, ξ) : RN × R × RN → RN est mesurable en x ∈ RN pour tout s ∈ R
et ξ ∈ RN et continue en ξ ∈ RN et s ∈ R pour presque tout x ∈ RN . Nous
noterons a(x, u, ∇u) = a(x, u(x), ∇u(x)).
et A vérifie aussi des conditions de coercitivité, monotonie et croissance : il existe p
vérifiant 2 − 1/N < p ≤ N
• Il existe α > 0 tel que pour tout s et ξ et presque tout x on ait
a(x, s, ξ)ξ ≥ α|ξ|p
• Pour tout s, ξ et η et presque tout x on a
[a(x, s, ξ) − a(x, s, η)](ξ − η) > 0 pour ξ 6= η
38 Chap. 5 Existence de solutions de problèmes elliptiques
0
• Il existe b(x) ∈ Lp (Ω) (p0 = p/(p − 1)) et β > 0 tels que pour tout s et ξ et
presque tout x on ait
Ces hypothèses sont classiques (sauf p > 2 − 1/N ) pour l’étude des opérateurs non
linéaires sous forme divergentielle (voir Leray-Lions [33]).
Pour k > 0, nous utiliserons la fonction (( tronquante )) Tk : R → R telle que
Tk (x) = min(k, max(x, −k)). Cette fonction étant lipschitzienne, si u ∈ W 1,p (Ω) on
a Tk (u) ∈ W 1,p (Ω) ∩ L∞ (Ω), cette dernière fonction pourra donc servir de fonction
test.
Nous noterons C toute constante indépendante de n, indice des suites qui
interviendrons ci-après.
avec f ∈ M (Ω) et g ∈ M (∂Ω), où n désigne, ici, la normale à ∂Ω. On résoudra (5.1)
au sens faible suivant : on cherche
\
u∈ W 1,q (Ω)
(p−1)N
q< N −1
tel que
[ Z Z Z
1,r
∀ϕ ∈ W (Ω) a(x, u, ∇u)∇ϕ = ϕ dg + ϕ df.
r>N Ω ∂Ω Ω
où l’on note encore ϕ, la trace de ϕ sur ∂Ω, ce que nous ferons désormais.
Bien entendu, si ce problème admet une solution, le choix ϕ ≡ 1 impose
Z Z
(5.2) dg + df = 0,
∂Ω Ω
ce que nous allons supposer. Ce problème admet alors, en particulier, une infinité de
solutions : plus précisément, il existe une solution de moyenne C, pour toute constante
C. Aussi cherche-t-on les solutions à moyenne donnée.
Théorème : Soient f ∈ M (Ω) (ne chargeant pas le bord), g ∈ M (∂Ω), vérifiant
(5.2), Ω un ouvert borné connexe, régulier de RN , A vérifiant les hypothèses (H) et
e ∈ R, alors le problème (5.1) admet une solution de moyenne u
soit u e.
5.1 Conditions aux limites de Neumann 39
0
Démonstration : Soient (fn ) ∈ W −1,p (Ω) ∩ L1 (Ω) et (gn ) ∈ (W 1−1/p,p (∂Ω))0 ∩
L1 (∂Ω) tels que kfn kL1 ≤ kf kM (Ω) , kgn kL1 ≤ kgkM (∂Ω) et fn * f dans M (Ω) faible
∗ et gn * g dans M (∂Ω) faible ∗ tels que
Z Z
gn + fn = 0.
∂Ω Ω
Étape 1 : Montrons que un est borné dans W 1,q (Ω) pour q <(p − 1)N/(N − 1),
ce qui donnera l’existence de u ∈ W 1,q (Ω) tel qu’une sous-suite de un converge vers
u dans W 1,q (Ω) faible. Rx
Pour m > 0 donné, soit ψm (x) = m 0 dt/(t + 1)m+1 pour x ≥ 0 et ψm (x) =
0
−ψm (−x) pour x ≤ 0 on a |ψm | ≤ 1, |ψm | ≤ m et ψm continue. Alors ψm (un ) ∈
1,p ∞
W (Ω) ∩ L (Ω), il est donc possible de choisir ϕ = ψm (un ) dans (5.3) qui devient
Z Z Z
a(x, un , ∇un )∇ψm (un ) = ψm (un )gn + ψm (un )fn
Ω ∂Ω Ω
0
or ∇ψm (un ) = ψm (un )∇un donc par coercitivité
Z
0
α |∇un |p ψm (un ) ≤ kfn kL1 + kgn kL1
Ω
0
et ψm (x) = m/(|x| + 1)m+1 donc
|∇un |p
Z
α m ≤ C,
Ω (|un | + 1)m+1
grâce à l’inégalité de Hölder nous obtenons pour q < p
q/p Z (p−q)/p
|∇un |p
Z Z
q (m+1)q/(p−q)
|∇un | ≤ m+1
(|un | + 1) .
Ω Ω (|un | + 1) Ω
Choisissons alors m > 0 tel que (m + 1)q/(p − q) < q ∗ (1/q ∗ = 1/q − 1/N ) ce qui est
possible pour q <(p − 1)N/(N − 1) alors
∗
(|un | + 1)(m+1)q/(p−q) ≤ ε|un |q + C(ε)
or la moyenne de un est u
e donc l’inégalité de Poincaré-Sobolev nous donne
kun − u
ekLq∗ ≤ Ck∇un kLq
soit
∗
u|(mes(Ω))1/q
kun kLq∗ ≤ Ck∇un kLq + |e
d’où grâce à (5.4)
Z 1/q∗ Z (p−q)/pq Z (p−q)/pq
q∗ (p−q)/pq q∗ 1 q∗
|un | ≤ Cε |un | +C ≤ |un | +C
Ω Ω 2 Ω
∗
pour ε assez petit. OrR (p − qq)/pq ≤ 1/q pour p ≤ N donc ku 1,q
n kLq∗ ≤ C et
donc, grâce à (5.4), Ω |∇un | ≤ C. Ainsi un est bornée dans W (Ω) pour tout
q <(p − 1)N/(N − 1), il est donc possible d’extraire de un une sous-suite convergeant
faiblement dans W 1,q (Ω).
{|(∇un − ∇u)(x)| ≥ δ} ⊂
{|un | ≥ k} ∪ {|u| ≥ k} ∪ {|∇un | ≥ k} ∪ {|∇u| ≥ k} ∪ {|un − u| ≥ η}
∪ {|(∇un − ∇u)(x)| ≥ δ, |un | ≤ k, |u| ≤ k, |∇un | ≤ k, |∇u| ≤ k, |un − u| ≤ η}
donc Z
1 C
mes(A1 ) ≤ |un | ≤ ≤ε
k Ω k
pour k assez grand, puisque un est borné dans W 1,q (Ω) pour q <(p − 1)N/(N − 1) et
donc dans L1 (Ω). De même ∇un est borné dans L1 (Ω) et u ∈ W 1,q (Ω) ⊂ W 1,1 (Ω),
donc les mesures de A2 , A3 et A4 sont majorées de la même façon. Fixons k tel que
chacune des mesures soit plus petite que ε.
5.1 Conditions aux limites de Neumann 41
or un converge vers u dans L1 (Ω), grâce au théorème de Rellich, donc pour η donné
il existe n1 tel que pour n ≥ n1 on ait
mes(A5 ) ≤ ε.
Z Z
γ≤ [a(x, un , ∇un ) − a(x, un , ∇Tk (u))](∇un − ∇Tk (u))1{|un −Tk (u)|≤η}
A6 A6
car sur A6 , |un − Tk (u)| = |un − u| ≤ η, de plus le terme intégré est positif et
∇Tη (un − Tk (u)) = (∇un − ∇Tk (u))1{|un −Tk (u)|≤η} , nous avons donc
Z Z
γ ≤ [a(x, un , ∇un ) − a(x, un , ∇Tk (u))]∇Tη (un − Tk (u))
A6
ZΩ Z
≤ a(x, un , ∇un )∇Tη (un − Tk (u)) − a(x, un , ∇Tk (u))∇Tη (un − Tk (u))
Ω Ω
et si l’on choisit ϕ = Tη (un − Tk (u)) ∈ W 1,p (Ω) ∩ L∞ (Ω) dans (5.3) on obtient
Z
a(x, un , ∇un )∇Tη (un − Tk (u)) ≤ η (kfn kL1 + kgn kL1 ) ≤ η C
Ω
Soit h > 0, choisissons maintenant ϕ = Th (un ) ∈ W 1,p (Ω) ∩ L∞ (Ω) dans (5.3), alors
Z
a(x, un , ∇un )∇Th (un ) ≤ h (kfn kL1 + kgn kL1 ) ≤ h C
Ω
42 Chap. 5 Existence de solutions de problèmes elliptiques
et par coercitivité
Z Z
a(x, un , ∇un )∇Th (un ) = a(x, un , ∇un )∇un 1|un |≤h
Ω Ω Z Z
≥ α |∇un | 1|un |≤h = α |∇Th (un )|p
p
Ω Ω
donc Th (un ) est borné dans W 1,p (Ω) pour tout h. Aussi pour h = k + η, Tk+η (un )
tend, à une sous-suite près, vers Tk+η (u) faiblement dans W 1,p (Ω) et donc forte-
ment dans Lp (Ω) et presque partout quand n → +∞. Donc a(x, Tk+η (un ), ∇Tk (u))
0
converge presque partout et |a(x, Tk+η (un ), ∇Tk (u))|p est equi-intégrable, aussi
0
grâce au théorème de Vitali, a(x, Tk+η (un ), ∇Tk (u)) converge fortement dans Lp (Ω)
vers a(x, Tk+η (u), ∇Tk (u)). Et comme 1{|Tk+η (un )−Tk (u)|≤η} tend presque partout vers
1{|Tk+η (u)−Tk (u)|≤η} et est borné, le théorème de convergence dominée de Lebesgue
donne la convergence forte de
a(x, Tk+η (un ), ∇Tk (u))1{|Tk+η (un )−Tk (u)|≤η} → a(x, Tk+η (u), ∇Tk (u))1{|Tk+η (u)−Tk (u)|≤η}
0
dans Lp (Ω).
De plus comme Tk+η (un ) converge faiblement dans W 1,p (Ω), ∇(Tk+η (un ) − Tk (u))
converge faiblement vers ∇(Tk+η (u) − Tk (u)) dans Lp (Ω) pour n → +∞ donc
Z
lim a(x, Tk+η (un ), ∇Tk (u))∇Tη (Tk+η (un ) − Tk (u))
n→+∞ Ω
Z
= a(x, Tk+η (u), ∇Tk (u))∇Tη (Tk+η (u) − Tk (u))
Ω
Étape 3 : Montrons que u est solution de (5.1). Soit ϕ ∈ W 1,r (Ω) avec r > N
alors ϕ ∈ W 1,p (Ω) ∩ L∞ (Ω), donc par définition de un
Z Z Z
a(x, un , ∇un )∇ϕ = ϕ gn + ϕ fn .
Ω ∂Ω Ω
ainsi |a(x, un , ∇un )−a(x, u, ∇u)|s est équi-intégrable et converge presque partout vers
0 donc converge dans L1 (Ω), grâce au théorème de Vitali, et donc a(x, un , ∇un ) →
a(x, u, ∇u) dans Ls (Ω) pour s < N/(N − 1), ainsi
Z Z
a(x, un , ∇un )∇ϕ → a(x, u, ∇u)∇ϕ
Ω Ω
et donc
[ Z Z Z
1,r
∀ϕ ∈ W (Ω) a(x, u, ∇u)∇ϕ = ϕ dg + ϕ df.
r>N Ω ∂Ω Ω
Remarquons que ∇un est borné dans W 1,q (Ω) pour q <(p−1)N/(N −1) aussi de la
même façon que ci-dessus, nous pouvons montrer que |∇(un − u)|q est équi-intégrable
et converge presque partout et donc, grâce au théorème de Vitali, que ∇un → ∇u
dans Lq (Ω) pour q <(p − 1)N/(N − 1).
44 Chap. 5 Existence de solutions de problèmes elliptiques
tel que
[ Z Z Z Z
1,r
∀ϕ ∈ W (Ω) a(x, u, ∇u)∇ϕ + λ uϕ = ϕ dg + ϕ df.
r>N Ω ∂Ω ∂Ω Ω
et il existe C2 tel que pour tout u ∈ W 1,q (Ω) (1 < q < N ) on ait
Z q/q∗ Z q Z
q∗ q
|u| ≤ C2 |u| + |∇u|
Ω ∂Ω Ω
R R R R
or Ω
u|q = |e
|e u|q |Ω| et |e
u|q |Ω|q ≤ ( Ω |u|)q donc Ω |e
u|q ≤ C( Ω |u|)q d’où
Z Z q Z
q
|u| ≤ C |u| + |∇u|q
Ω Ω Ω
Étape 1 : Montrons que un est borné dans W 1,q (Ω) pour q <(p − 1)N/(N − 1).
La méthode est la même que pour les conditions de Neumann,R mais l’inégalité de
Poincaré-Sobolev que nous venons de démontrer fait intervenir ∂Ω |un |, il nous faut
donc une estimation de cette intégrale. Comme Tk (un ) ∈ W 1,p (Ω) ∩ L∞ (Ω) nous
pouvons choisir ϕ = Tk (un ) dans (5.6), ce qui nous donne
Z Z Z Z
0
a(x, un , ∇un )∇un Tk (un ) + λ un Tk (un ) = Tk (un ) gn + Tk (un ) fn
Ω ∂Ω ∂Ω Ω
soit Z
1
λ un Tk (un ) ≤ C
∂Ω k
46 Chap. 5 Existence de solutions de problèmes elliptiques
ε0
Z
a(x, un , ∇Tk (u))∇Tη (un − Tk (u)) ≤
Ω 2
Étape 3 : Comme pour les conditions aux limites R de Neumann nous R avons,
1,r
pour
R ϕ ∈ W (Ω) avec r > N , les convergences
R de Ω
a(x, u n , ∇u n )∇ϕ, Ω
ϕ fn et
1,q
∂Ω
ϕ gn , il reste donc à montrer celle de ∂Ω un ϕ. Comme un → u dans W (Ω) pour
q <(p − 1)N/(N − 1) la trace de un sur ∂Ω converge vers celle de u dans W 1−1/q,q (∂Ω)
et donc dans Lq (∂Ω), comme ϕ ∈ W 1,r (Ω) pour r > N , ϕ ∈ C(Ω) donc ϕ ∈ L∞ (∂Ω)
et nous obtenons Z Z
un ϕ → uϕ
∂Ω ∂Ω
[ Z Z
∀ϕ ∈ W01,r (Ω) a(x, u, ∇u)∇ϕ = ϕ df.
r>N Ω Ω
48 Chap. 5 Existence de solutions de problèmes elliptiques
et un − g ∈ W01,p (Ω).
Étape 1 : Nous montrons, là encore, que un est borné dans W 1,q (Ω) pour
q <(p−1)N/(N −1). La démonstration ne peut pas être ici faite avec ϕ = ψm (un −e
g ),
aussi allons nous adopter une variante de la démonstration correspondant aux
conditions aux limites de Neumann.
Soit k ∈ N, choisissons ϕ = Tk (un − ge) ∈ W01,p ∩ L∞ (Ω) alors
Z Z
a(x, un , ∇un )∇(un − ge)1{|un −eg|≤k} = Tk (un − ge) fn
Ω Ω
donc
Z Z Z
a(x, un , ∇un )∇un 1{|un −eg|≤k} = Tk (un − ge) fn + a(x, un , ∇un )∇e
g 1{|un −eg|≤k}
Ω Ω Ω
car (1{|un −eg|≤k} )p/(p−1) = 1{|un −eg|≤k} . L’inégalité de Poincaré nous donne, appliquée à
Tk (un − ge) ∈ W01,p (Ω),
Z Z Z
p p
|un − ge| 1{|un −eg|≤k} ≤ |Tk (un − ge)| ≤ C |∇(un − ge)|p 1{|un −eg|≤k}
Ω Ω Ω
si bien que
Z Z Z
p p p
|un | 1{|un −eg|≤k} ≤ C g | + |∇e
(|e g| ) + C |∇un |p 1{|un −eg|≤k}
Ω Ω Ω
5.3 Conditions aux limites de Dirichlet 49
alors
Z
α |∇un |p 1{|un −eg|≤k} ≤ kkf kM (Ω) + βCke g kW 1,p )p−1 ]
g kW 1,p [kbkW 1,p0 + (ke
Ω
Z (p−1)/p
p
g kW 1,p
+2βCke |∇un | 1{|un −eg|≤k}
Ω
donc Z Z (p−1)/p
1 p 1 p
|∇un | 1{|un −eg|≤k} ≤C +C |∇un | 1{|un −eg|≤k}
k Ω k Ω
et (p − 1)/p < 1 donc il existe C indépendant de n et de k tel que
Z
|∇un |p 1{|un −eg|≤k} ≤ C k
Ω
R
or Ω
g |p 1{|un −eg|≤k} ≤ (k∇e
|∇e g kLp )p donc
Z
(5.10) |∇(un − ge)|p 1{|un −eg|≤k} ≤ C k
Ω
et donc
k−1 Z
X
(5.11) |∇(un − ge)|p 1{l≤|un −eg|≤l+1} ≤ C k
l=0 Ω
posons Z
1
ak = et bk = |∇(un − ge)|p 1{k≤|un −eg|≤k+1}
(1 + k)m+1 Ω
P Pk−1
ce qui nous amène à considérer k∈N ak bk , notons Bk = l=0 bl pour k ≥ 1 et
B0 = 0 alors pour N ∈ N
N
X N
X N
X N
X
ak b k = ak (Bk+1 − Bk ) = ak Bk+1 − ak Bk
k=0 k=0 k=0 k=0
N
X N
X
= ak Bk+1 − ak+1 Bk+1 + aN +1 BN +1
k=0 k=0
XN
= (ak − ak+1 ) Bk+1 + aN +1 BN +1
k=0
50 Chap. 5 Existence de solutions de problèmes elliptiques
et pour q <(p − 1)N/(N − 1), on conclut de la même façon que précédemment que
un − ge est borné dans W01,q (Ω) et donc que un est borné dans W 1,q (Ω).
|∇vn |q ≤ C 0 donc (vn ) est bornée dans W01,q (Ω) pour tout
R
pour tout k ∈ N alors Ω
q <(p − 1)N/(N − 1).
or nous avons montré que Th (un − ge) est borné dans W01,p (Ω) donc ∇(Tk+η (un − ge) −
Tk (u − ge)) converge faiblement dans Lp (Ω) vers ∇(Tk+η (u − ge) − Tk (u − ge)) et nous
pouvons faire le même raisonnement que celui qui a été fait pour les conditions de
Neumann, avec un − ge : 1{|Tk+η (un −eg)−Tk (u−eg)|≤η} tend pp vers 1{|Tk+η (u−eg)−Tk (u−eg)|≤η}
0
et a(x, Tk+η (un − ge) + ge, ∇Tk (u)) converge fortement dans Lp (Ω) vers a(x, Tk+η (u −
ge) + ge, ∇Tk (u)) pour n → +∞. Enfin ∇Tη (Tk+η (u − ge) − Tk (u − ge)) → 0 pp pour
η → 0, on peut donc choisir η < ε0 /2C tel que pour n assez grand
ε0
Z
a(x, un , ∇un )∇Tη (un − ge − Tk (u − ge)) ≤
Ω 2
et
ε0
Z
a(x, un , ∇Tk (u))∇Tη (un − ge − Tk (u − ge)) ≤
Ω 2
donc de même que pour le problème de Neumann on a mes(A6 ) ≤ ε pour n assez
grand, ce qui achève la preuve de la convergence.
pour α > 0 et f ∈ M (Ω) = (C(Ω))0 espace de mesures qui est le dual de l’espace des
fonctions continues sur Ω muni de sa norme habituelle.
Ce problème admet, pour f ∈ H −1 (Ω), une unique solution variationnelle, qui
appartient à H01 (Ω) et vérifie
Z
1
∀ v ∈ H0 (Ω) (A∇u)∇v dx = hf, viH −1 ,H01
Ω
Pour Tf 6∈ H −1 (Ω), mais f ∈ M (Ω) on ne trouve plus de solutions dans H01 (Ω), mais
dans q < N W01,q (Ω), ce qui conduit à une formulation plus faible, puisque pour que
N −1
la première intégrale ait un sens il faut alors que v ∈ p > N W 1,p (Ω) :
S
[ Z Z
1,p
(6.2) ∀v ∈ W0 (Ω) (A∇u)∇v dx = v df.
p>N Ω Ω
54 Chap. 6 Non unicité des solutions . . .
Nous avons W01,p (Ω) ⊂ C(Ω) pour p > N aussi le membre de droite a un sens.
Cette formulation étant plus faible que la formulation variationnelle elle n’assure
plus l’unicité, comme le montre le contre-exemple de Serrin [46] que nous allons
présenter.
Pour N ≥ 2, nous avons étudié, au chapitre 4 deux constructions d’une solution
de (6.1) au sens de (6.2), cette solution appartient à W01,q (Ω), ∀q < NN−1 . Nous allons
montrer que, pour N > 2 et Ω = BRN (la boule unité de RN ), il n’ y a pas, dans cet
espace, unicité des solutions vérifiant (6.2).
L’idée est d’adapter le contre-exemple de Serrin [46] (que nous présentons dans
le paragraphe 6.2) pour construire une solution u de l’équation de (6.1) pour f = 0
telle que u ∈ W 1,q (BRN ) pour tout q < NN−1 et u 6∈ H 1 (Ω) mais dont la trace sur
1
S N −1 la sphère unité de RN est dans H 2 (S N −1 ). Nous construisons cette solution au
paragraphe 6.3 et montrons sa régularité au paragraphe 6.4.
Grâce à la régularité de u sur S N −1 , nous pouvons considérer le problème de
Dirichlet suivant
− div(A∇v) = 0 dans BRN
(6.3)
v = u sur S N −1
qui admet une unique solution (variationnelle) v dans H 1 (BRN ). Comme u 6∈
H 1 (BRN ), on a donc construit deux solutions distinctes de (6.3) qui appartiennent à
1,q
T
q< N W (BRN ). Soit w = u − v, par linéarité, w vérifie
N −1
W01,q (BRN ) et
T
c’est-à-dire que w ∈ q< N
N −1
[ Z
∀ϕ ∈ W01,p (BRN ) (A∇w)∇ϕ = 0
p>N B RN
(c’est la formulation (6.2)) et donc, toujours par linéarité, (6.1) admet au moins deux
solutions dans W 1,q (BRN ) pour tout q < NN−1 : celle de [11], notée u1 , et u1 + w.
Remarque : Cette adaptation n’est donc possible que pour N > 2. Ceci est normal
puisque l’on sait, grâce au lemme de Meyers [37], que
T pour N = 2 le problème de
Dirichlet homogène admet une unique solution dans q < N W01,q (Ω) (voir [29]).
N −1
car ∂xi (aij ∂xj u) = 0 dans r > η puisque u est solution classique de (6.4) dans
{x ∈ RN | r > η}, pour tout η > 0. De plus, les calculs conduisent à
I I
xi 1 x1
ϕ aij ∂xj u dσ = ϕ r−1−ε dσ
r=η r ε r=η r
et comme ϕ est C 1
avec
sup |f (x1 , . . . , xN )| = O(η)
r=η
donc
I I I
x1 x1 x1 −1−ε
ϕ r−1−ε dσ = ϕ(0, . . . , xN ) r−1−ε dσ + f (x1 , . . . , xN ) r dσ.
r=η r r=η r r=η r
Le premier terme est nul car c’est l’intégrale d’une fonction impaire en x1 sur un
domaine symétrique en x1 . Pour le second passons en coordonnées (( cylindriques ))
avec dσ = rdθ dz,
I I
x1 x1
f (x) 2+ε dσ = f (x) 2+ε r g(θ, z) dθ dz
r=η r r=η r
I
1
≤ ε O(η) |g(θ, z)| dθ dz
η r=η
= O(η 1−ε ).
si bien que
Z Z
1−ε
aij ∂xj u ∂xi ϕ dx = O(η ) d’où aij ∂xj u ∂xi ϕ dx = 0,
r>η Ω
6.4 Régularité de u
N
u ∈ W 1,q (Ω) ∀q< .
N −1
6.4 Régularité de u 57
Afin de pouvoir construire une solution v dans H 1 (Ω) vérifiant v = u sur ∂Ω,
1
il faut montrer que u ∈ H 2 (∂Ω). Pour que Ω soit de bord régulier nous choisissons
Ω = BRN , la boule unité de RN si bien que ∂Ω = S N −1 la sphère unité de RN .
L’ouvert Ω = BRN ainsi choisi a la régularité C m pour tout m ∈ N et une
1
frontière bornée, aussi pouvons-nous définir H 2 (S N −1 ).
1
Propriété : La fonction u appartient à H 2 (S N −1 ) pour ε < N 1−1 .
où y 0 = (y1 , . . . , yN −1 ), alors u ∈ W s,p (∂Ω) si θj u ∈ W s,p (RN −1 ). Cette définition est
indépendante du choix des (Uj ), (Ψj ) et (ωj ).
et si v ne dépend que de x1 et x2 on a
Z
N −3
I(v) ≤ M (2γ) |v(x1 , x2 )|β dx1 dx2
[−γ,γ]2
2
or nous avons vu que u et ux ∈ Lβ (Ω2 ) pour β < 1+ε donc I(u) < +∞ et I(ux ) < +∞,
0 0
c’est-à-dire que u(Ψj (y , 0)) et ux (Ψj (y , 0)) ∈ L (BRN −1 ) or ωj est C 1 donc θj u ∈
β
2
W 1,β (BRN −1 ) pour β < 1+ε .
1 1
Afin de montrer que u ∈ H 2 (S N −1 ), il reste à montrer que θj u ∈ H 2 (RN −1 ) pour
cela, utilisons un théorème d’injection de Sobolev (voir [1]) :
1 1 1 1
W 1,q (RN −1 ) ⊂ W 2 ,p (RN −1 ) pour = −
p q 2(N − 1)
1
comme θj u ∈ W 1,β (RN −1 ) pour tout β < 1+ε 2
on déduit que θj u ∈ W 2 ,p (RN −1 ) pour
−2 1 −2
tout p < 2/( NN −1
+ ε). Donc θj u ∈ H 2 (RN −1 ) si 2/( N
N −1
+ ε) > 2 ce qui est le cas pour
1
1 N −1 1
ε < N −1 . Ainsi u ∈ H 2 (S ) pour ε < N −1 (la démonstration montre un peu plus
1,β N −1 2
c’est-à-dire que u ∈ W (S ) pour tout β < 1+ε ).
6.5 Étude des tronqués de u 59
Afin d’obtenir l’unicité des solutions de (6.1) dans le cas f ∈ L1 , il est demandé
dans [4], pour des conditions homogènes de Dirichlet, que la solution u vérifie
Tk (u) ∈ H01 (Ω) pour tout k > 0 et que u vérifie une inégalité supplémentaire. Lorsque
la condition aux limites n’est plus homogènes, on montre dans [41] et au chapitre 5
que Tk (u − ϕ) ∈ H01 (Ω) pour k > 0 et ϕ un relèvement H 1 (Ω) de la condition aux
limites.
Nous allons donc montrer que le contre-exemple de Serrin, noté u au paragraphe
6.3, u vérifie Tk (u) 6∈ H 1 (BRN ) pour tout k ≥ 0 et qu’il existe ϕ ∈ H 1 (Ω) valant u
sur le bord et tel que Tk (u − ϕ) 6∈ H01 (Ω). Pour cela montrons que l’intégrale
Z Z
2
k∇Tk (u)k = k∇(u)k2 1{|u|≤k}
BRN BRN
Passons en polaire, nous obtenons u = cos θ/ρε donc nous avons |u| ≤ k pour
ρε ≥ | cos θ|/k soit ρ ≥ (| cos θ|/k)1/ε que nous noterons ρ(θ), d’où
√
Z Z 2πZ 1/ 2
2
k∇uk 1{|u|≤k} dx1 dx2 = k∇uk2 ρ dρ dθ
x1 2 +x2 2 ≤ 12 0 ρ(θ)
or
1
k∇uk2 = 2
(Ax1 4 + Bx1 2 x2 2 + Cx2 4 )
(x1 + x2 2 )3+ε
où A, B et C dépendent de ε. Considérons chacun des trois termes
√
2πZ 1/ 2 2π
x1 4 cos4 θ
Z Z Z
dx1 dx2 = dρ dθ = K + cos2 θ dθ < +∞
x1 2 +x2 2 ≤ 1
2 (x1 2 + x2 2 )3+ε 0 ρ(θ) ρ1+2ε 0
|u|≤k
et de même
2π
x1 2 x2 2
Z Z
dx1 dx2 = K 0 + sin2 θ dθ < +∞
1
x1 2 +x2 2 ≤ 2 (x1 2 + x2 2 )3+ε 0
|u|≤k
60 Chap. 6 Non unicité des solutions . . .
et
2π
x2 4 sin4 θ
Z Z
dx1 dx2 = K 00 + dθ
1
x1 2 +x2 2 ≤ 2 (x1 2 + x2 2 )3+ε 0 cos2 θ
|u|≤k
π π
or cos θ ∼ − θ en donc la dernière intégrale est divergente, aussi
2 2
Z
k∇uk2 1{|u|≤k} dx1 dx2 = +∞
x1 2 +x2 2 ≤ 12
rayon 1/2, on obtient un relèvement (encore noté ϕ) nul dans cette même boule, le
calculs ci-dessus donnent donc, de même Tk (u − ϕ) = Tk (u) 6∈ H01 (B1/2 ) et donc
Tk (u − ϕ) 6∈ BRN .
Part III
61
7. Dépendance continue par rapport à un
paramètre des solutions de problème elliptique à
seconds membres L1
0
où h., .i est le crochet de dualité entre W01,p (Ω) et W −1,p (Ω).
Lorsque f ∈ M (Ω) = (C(Ω))0 , espace de mesures, l’existence d’une solution
vérifiant \
u ∈ W01,q (Ω)
q < (p−1)N
(N −1)
(7.3) Z Z [
W01,r (Ω)
a(x, α, ∇u α )∇ϕ = ϕ df ∀ϕ ∈
Ω Ω r>N
a été montrée par Boccardo et Gallouët [11, 12]. Cependant il n’y a pas unicité des
solutions au sens ci-dessus comme le montre le contre-exemple adapté de celui de
Serrin (voir [46, 43] et chapitre 6). Aussi lorsque f appartient seulement à L1 (Ω)
faut-il considérer les solutions au sens suivant (Tk fonction (( tronquante )) est définie
64 Chap. 7 Dépendance continue par rapport à un paramètre
par la suite)
\
u ∈ W01,q (Ω) Tk (u) ∈ W01,p (Ω) ∀k > 0
q < (p−1)N
(7.4) Z (N −1)
Z
a(x, α, ∇uα )∇Tk (uα − ϕ) ≤ f Tk (uα − ϕ) ∀ϕ ∈ L∞ (Ω) ∩ W01,p (Ω),
Ω Ω
dont l’unicité a été montrée par Bénilan, Boccardo, Gallouët, Gariepy, Pierre et
Vazquez [4].
Donnons les hypothèses sur a, opérateur de Leray-Lions, nous les appellerons (H)
par la suite : a est une fonction de Carathéodory, c’est-à-dire que
• a(x, α, ξ) : RN × R × RN → RN est mesurable en x ∈ RN pour tout α ∈ R
et ξ ∈ RN et continue en ξ ∈ RN et s ∈ R pour presque tout x ∈ RN . Nous
noterons a(x, α, ∇uα ) = a(x, α, ∇uα (x)).
et a vérifie aussi des conditions de coercitivité, monotonie et croissance : il existe p
vérifiant 2 − 1/N < p ≤ N
• Il existe γ > 0 tel que pour tout α et ξ et presque tout x on ait
Ces hypothèses sont classiques (sauf p > 2 − 1/N ) pour l’étude des opérateurs non
linéaires sous forme divergentielle (voir Leray-Lions [33]).
Nous ajoutons une forme un peu plus forte de continuité par rapport à α :
0
7.1 Cadre variationnel : f ∈ W −1,p (Ω)
et
hf, uα i ≤ kf kW −1,p0 kuα kW 1,p
0
donc
1
(kuα kW 1,p )p−1 ≤ kf kW −1,p0
0 γ
c’est-à-dire que (uα ) est bornée dans W01,p (Ω) indépendamment de α.
Choisissons ϕ = uα − uα0 dans (7.2) écrit pour α et α0 , alors après soustraction
Z
[a(x, α0 , ∇uα ) − a(x, α0 , ∇uα0 )]∇(uα − uα0 )
Ω
Z
= [a(x, α0 , ∇uα ) − a(x, α, ∇uα )]∇(uα − uα0 )
Ω
vérifie Z
lim gα = 0
α→α0 Ω
et est positif donc (gα ) converge vers 0 dans L1 (Ω). Une sous-suite encore notée (gα )
converge donc presque partout vers 0, donc pour presque tout x fixé, par stricte
monotonie, ∇uα (x) → ∇uα0 (x) (par l’absurde grâce à la continuité de a par rapport
à ξ et α). Ainsi (∇uα ) converge presque partout vers ∇uα0 .
De plus (ua ) est bornée dans W01,p (Ω) donc (par unicité de la limite) converge
faiblement vers uα0 dans W01,p (Ω) lorsque α → α0 donc hf, uα i → hf, uα0 i et donc
Z Z
a(x, α, ∇uα )∇uα → a(x, α0 , ∇uα0 )∇uα0 .
Ω Ω
Comme a(x, α, ∇uα )∇uα ≥ 0, ceci implique la convergence de a(x, α, ∇uα )∇uα dans
L1 (Ω) or
γ|∇uα |p ≤ a(x, α, ∇uα )∇uα
donc la suite (γ|∇uα |p ) est dominée par une suite qui converge dans L1 et est donc
équi-intégrable et (|∇uα − ∇uα0 |p ) aussi, comme cette dernière suite converge de plus
presque partout, le théorème de Vitali entraı̂ne que ∇uα → ∇uα0 dans Lp (Ω) (fort).
Ainsi (ua ) converge fortement vers uα0 dans W01,p (Ω).
Par unicité de la limite uα0 , on a en fait la convergence de toute la suite (uα ) de
départ (par l’absurde).
Remarque 2 : Nous venons de montrer que (uα ) converge fortement vers uα0
dans W01,p (Ω), et en particulier converge faiblement donc l’opérateur div(a(x, α, ξ)ξ)
H-converge vers div(a(x, α, ξ)ξ).
pour tout ϕ ∈ L∞ (Ω) ∩ W01,p (Ω). L’unicité est ici essentielle : s’il n’y avait pas unicité,
la continuité n’aurait aucun sens. Aussi ne peut-on pas considérer le cas f mesure,
pour lequel la formulation au sens des distributions n’assure pas l’unicité (voir ci-
dessus).
Démonstration : Ici aussi on peut montrer que (uα ) est borné dans W01,q (Ω)
pour q < N/(N − 1) indépendamment de α, voir [11, 12].
De façon analogue à ce que nous venons de faire, choisissons ϕ = Th (uα0 ) dans
l’équation en uα et ϕ = Th (uα ) dans celle en uα0 (ceci reprend partiellement [4]), et
faisons la somme
Z
[a(x, α, ∇uα ) − a(x, α, ∇uα0 )]∇Tk (uα − Th (uα0 ))
Ω
Z Z
≤ − a(x, α, ∇uα0 )∇Tk (uα − Th (uα0 )) − a(x, α0 , ∇uα0 )∇Tk (uα0 − Th (uα ))
Z Ω Ω
Z
C= f [Tk (uα − uα0 ) + Tk (uα0 − uα )]
|uα |≤h
|uα0 |≤h
Z
+ f [Tk (uα − Th (uα0 )) + Tk (uα0 − Th (uα ))]
|uα |≥h ou |uα0 |≥h
et Z
A2 = a(x, α, ∇uα0 )∇Tk (uα − Th (uα0 ))
|uα0 |≥h
alors
Z
[a(x, α, ∇uα ) − a(x, α, ∇uα0 )]∇Tk (uα − Th (uα0 )) ≤ A − A1 + A2 .
|uα0 |≤h
Enfin, soit Z
B1 = a(x, α, ∇uα0 )∇Tk (uα − Th (uα0 )),
|uα0 |≥h
Z
B2 = a(x, α, ∇uα0 )∇Tk (uα − Th (uα0 )),
|uα |≥h
|uα0 |≤h
Z
B3 = a(x, α0 , ∇uα0 )∇Tk (uα0 − Th (uα )),
|uα |≥h
Z
B4 = a(x, α0 , ∇uα0 )∇Tk (uα0 − Th (uα )).
|uα0 |≥h
|uα |≤h
68 Chap. 7 Dépendance continue par rapport à un paramètre
Z
B≤ [a(x, α, ∇uα0 ) − a(x, α0 , ∇uα0 )]∇Tk (uα0 − uα ) − B2 − B4
|uα |≤h
|uα0 |≤h
Donc
Z
[a(x, α, ∇uα ) − a(x, α, ∇uα0 )]∇Tk (uα − Th (uα0 )) ≤ C − B2 − B4
|uα0 |≤h
Z
+ [a(x, α, ∇uα0 ) − a(x, α0 , ∇uα0 )]∇Tk (uα0 − uα )
|uα |≤h
|uα0 |≤h
or |[a(x, α, ∇uα0 )−a(x, α0 , ∇uα0 )]∇(uα0 −uα )| ≤ c(α−α0 )(b(x)+|∇uα0 |p−1 )|∇(uα0 −
uα )|
Z
[a(x, α, ∇uα ) − a(x, α, ∇uα0 )]∇Tk (uα − Th (uα0 )) ≤ C − B2 − B4
|uα0 |≤h
Z 1/p0 Z !1/p
p0
+c(α − α0 ) |b(x)| + |∇uα |p |∇Tk (uα − uα0 )|p
|uα0 |≤h
|uα |≤h |uα |≤h
Pour cela, fixons δ > 0 et ε > 0, et remarquons que pour k > 0 et h > 0 nous avons
{x; |(∇uα − ∇uα0 )(x)| ≥ δ} ⊂
{|uα0 | ≥ h} ∪ {|∇uα | ≥ h} ∪ {|∇uα0 | ≥ h} ∪ {|uα − uα0 | ≥ k}
∪ {|(∇uα − ∇uα0 )| ≥ δ, |uα0 | ≤ h, |∇uα | ≤ h, |∇uα0 | ≤ h, |uα − uα0 | ≤ k}
7.2 Cadre non variationnel : f ∈ L1 (Ω) 69
donc Z
1 C
mes(D1 ) ≤ |uα0 | ≤ ≤ε
h Ω h
pour h assez grand, puisque uα0 est borné dans Lq (]0, T [×Ω) pour q <(p−1)N/(N −1)
et donc dans L1 (Ω). Fixons donc h tel que mes(D1 ) ≤ ε.
Majorons mes(D2 ) et mes(D3 ), nous avons
Z Z
|∇uα | ≥ |∇uα | ≥ h mes(D2 )
Ω D1
donc Z
1 C
mes(D2 ) ≤ |∇uα | ≤ ≤ε
h Ω h
pour h assez grand, puisque (∇uα ) est borné dans Lq (]0, T [×Ω) pour q <(p−1)N/(N −
1) et donc dans L1 (Ω). Fixons donc h tel que mes(D2 ) ≤ ε et mes(D3 ) ≤ ε.
Majorons maintenant mes(D4 ), nous avons
Z Z
|(uα − uα0 )| ≥ |(uα − uα0 )| ≥ k mes(D4 )
Ω D3
Comme (uα ) est borné dans W 1,q (Ω) et donc dans L1 (Ω) indépendamment de α on a
Z
2C ≥ |(uα − uα0 )| ≥ k mes(D4 )
Ω
γ et de D5 , nous avons
Z Z
γ≤ [a(t, x, ∇uα ) − a(t, x, ∇uα0 )](∇uα − ∇uα0 )1{|uα −uα0 |≤η}
D5 D5
Or on a vu que
Z
[a(x, α, ∇uα ) − a(x, α, ∇uα0 )]∇Tk (uα − Th (uα0 )) ≤ C − B2 − B4
|uα0 |≤h
Z 1/p0 Z !1/p
p0 p p
+c(α − α0 ) |b(x)| + |∇uα | |∇Tk (uα − uα0 )|
|uα0 |≤h
|uα |≤h |uα |≤h
R
et D5 ⊂ {|uα0 | ≤ h} donc ceci majore D5 γ.
Soit v la limite de uα presque partout grâce à Rellich, k a déjà été fixé, nous
allons choisir un nouvel h plus petit que le précédent (majoration de D1 et D2 ) tel
que Z
|f | ≤ ε00
|v|≥h ou |uα0 |≥h
donc
Z Z
p0 0
|a(x, α, ∇uα0 )| ≤ (b(x) + |∇uα0 |p−1 )p
h−k≤|uα0 |≤h h−k≤|uα0 |≤h
Z
0
≤ (b(x) + (|f |/γ)1/(p−1) )p
h−k≤|uα0 |
Z !1/p0 Z 1/p
p0 p
≤ |a(x, α, ∇uα0 )| |∇uα |
h−k≤|uα0 |≤h h≤|uα |≤h+k
car
{|uα | ≥ h, |uα0 | ≤ h, |uα − Th (uα0 )| ≤ k} ⊂ {h − k ≤ |uα0 | ≤ h}
{|uα | ≥ h, |uα0 | ≤ h, |uα − Th (uα0 )| ≤ k} ⊂ {h ≤ |uα | ≤ h + k}
Ainsi pour h assez grand |B5 | ≤ Cε00 et de même pour B6 , |B7 |, B8 . Pour ce h fixé,
on a Z
0
|b(x)|p + |∇uα |p ≤ Ch
|uα |≤h
et Z
|∇Tk (uα − uα0 )|p ≤ Ch
|uα0 |≤h
|uα |≤h
et pour |α − α0 | ≤ α1 on a aussi
Z Z
|f | − |f | ≤ ε00
|uα |≥h ou |uα0 |≥h |v|≥h ou |uα0 |≥h
Les mesures des quatre premiers ensembles sont majorées de la même façon que ci-
dessus, celle du cinquième est petite pour η fixé et α − α0 petit car uα tend vers u
presque partout grâce au théorème de Rellich, il reste donc à fixer η et à majorer la
mesure du sixièmeR ensemble, que nous noterons D6 . Pour cela, et comme ci dessus il
suffit de majorer D6 γ :
Z Z
γ≤ [a(x, α, ∇uα ) − a(x, α, ∇u)]∇(uα − u).
D6 D6
et pour η petit Z Z
f Tη (uα − Th (u)) ≤ η |f | ≤ ε,
Ω Ω
ainsi Z
[a(x, α, ∇uα )∇(uα − u) ≤ ε.
D6
donc Z
a(x, α, ∇Th (u))∇(Th (uα ) − Th (u)) ≤ ε
D6
Ce chapitre est la version française d’un article à paraı̂tre à Nonlinear Analysis TMA
[43].
La formulation classique des problèmes paraboliques dans le cas où les données
sont dans L1 n’assure pas l’unicité des solutions aussi nous présentons ici une
formulation appelée entropique permettant d’obtenir existence et unicité.
avec u0 dans L1 (Ω) et f dans L1 (]0, T [×Ω) où Ω est un ouvert borné de RN et a une
fonction de Carathéodory satisfaisant des conditions de coercitivité, de monotonie
et de croissance du type de celles de Leray-Lions, définissant un opérateur sur
Lp (]0, T [; W01,p (Ω)).
Boccardo et Gallouët [11] ont montré, dans le cas plus général où u0 ∈ M (Ω),
f ∈ M ([0, T ] × Ω) (pour O un ouvert M (O) = (C(O))0 est le dual de l’ensemble des
fonctions continues sur O muni de sa norme habituelle, c’est un espace de mesures),
et p > 2 − 1/(N + 1), qu’il existait une solution au sens suivant
\
u∈ Lq (]0, T [, W01,q (Ω))
(p−1)(N +1)+1
q< N +1
et Z TZ Z Z TZ Z TZ
(8.2) − u ϕt + ϕ(0) du0 + a(t, x, ∇u)∇ϕ = ϕ df
0 Ω Ω 0 Ω 0 Ω
− div(A∇v) = 0 dans Ω
au sens de (8.2).
Afin de remédier à la non-unicité, une formulation entropique est proposée, celle-
ci est proche de celle qui a été introduite pour le cas elliptique dans [4]. Dans le cas
où a ne dépend pas de t, l’existence et l’unicité de la solution entropique ont été
montrées, en utilisant la théorie des semi-groupes, par Andreu, Mazón, Segura de
León et Toledo [2].
Une autre formulation définissant des solutions dites renormalisées et assurant
aussi l’unicité est donnée par Blanchard et Murat [8]. Ces deux formulations con-
duisent à la même solution.
Donnons les hypothèses sur a, nous les appellerons (H) par la suite : a est une
fonction de Carathéodory, c’est-à-dire que
• a(t, x, ξ) : R × RN × RN → RN est mesurable en t ∈ R et x ∈ RN pour tout
ξ ∈ RN et continue en ξ ∈ RN pour presque tout t ∈ R et x ∈ RN . Nous
noterons a(t, x, ∇u) = a(t, x, ∇u(t, x)).
et a vérifie aussi des conditions de coercitivité, monotonie et croissance : il existe p
vérifiant 2 − 1/(N + 1) < p ≤ N et
8.1 Position du problème 77
Ces hypothèses sont classiques (sauf p > 2 − 1/(N + 1)) pour l’étude des opérateurs
non linéaires sous forme divergentielle (voir Leray et Lions [33]).
et sa primitive Θk : R → R+
Z x
Θk (x) = Tk (s) ds.
0
qui implique Z T Z Z
hvt , Tk (v)i = Θk (v(T )) − Θk (v(0)) .
0 Ω Ω
0
h., .i représente la dualité entre W −1,p (Ω) et W01,p (Ω). Il en sera de même par la suite
sauf mention contraire.
Définition : Pour f ∈ L1 (]0, T [×Ω), u0 ∈ L1 (Ω) et Ω ouvert borné de RN , on
appelle solution entropique de (8.1) une fonction u ∈ C([0, T ]; L1 (Ω)) telle que pour
tout k > 0 on ait Tk (u) ∈ Lp (]0, T [; W01,p (Ω)) et qui vérifie
Z Z Z T
Θk (u − ϕ)(T ) − Θk (u0 − ϕ(0)) + hϕt , Tk (u − ϕ)i
Ω Ω 0
Z TZ Z TZ
(8.3) + a(t, x, ∇u)∇Tk (u − ϕ) ≤ f Tk (u − ϕ)
0 Ω 0 Ω
pour tout k > 0 et ϕ ∈ Lp (]0, T [; W01,p (Ω)) ∩ L∞ (]0, T [×Ω) ∩ C([0, T ]; L1 (Ω)) tel que
0 0
ϕt ∈ Lp (]0, T [; W −1,p (Ω)).
78 Chap. 8 Solutions entropiques de problèmes paraboliques
Nous rappelons tout d’abord comment est obtenue la solution de [11] et nous
précisons sa régularité et, dans les deux paragraphes suivants, nous montrons
l’existence puis l’unicité de la solution entropique.
Plus précisément, ici f ∈ L1 (]0, T [×Ω) ⊂ M ([0, T ]×Ω) soit alors (fn ) une suite de
0 0
fonctions de L1 (]0, T [×Ω) ∩ Lp (]0, T [; W −1,p (Ω)) telle que kfn kL1 ≤ kf kL1 et fn → f
dans L1 (]0, T [×Ω) et soit un0 une suite de L2 (Ω) ⊂ L1 (Ω) telle que kun0 kL1 ≤ ku0 kL1
et un0 → u0 dans L1 (Ω). Soit un la solution (( classique )) (voir Lions [34]) de
et que un vérifie
Z T Z TZ Z T
(8.4) hunt , ϕi + a(t, x, ∇un )∇ϕ = hfn , ϕi
0 0 Ω 0
pour tout ϕ ∈ Lp (]0, T [; W01,p (Ω)) et que un (0) = un0 . Il est montré dans [11] que (un )
est borne dans Lq (]0, T [, W01,q (Ω)) pour tout q < (p−1)(N +1)+1
N +1
et que Tk (un ) est borné
p 1,p
dans L (]0, T [, W0 (Ω)) pour tout k > 0 aussi T (à une sous suite près) il existe u tel
que (un ) converge faiblement vers u dans q <[(p−1)(N +1)+1]/(N +1) L (]0, T [, W01,q (Ω))
q
et Tk (un ) converge faiblement vers Tk (u) dans Lp (]0, T [, W01,p (Ω)) pour tout k > 0.
Pour montrer que u est une solution de (8.1), il nous reste à montrer la propriété
suivante.
∀δ > 0, ∀ε > 0, ∃n0 tel que ∀n, m ≥ n0 mes{(t, x); |(∇un − ∇um )(t, x)| ≥ δ} ≤ ε.
Pour cela, fixons δ > 0 et ε > 0, et remarquons que pour k > 0 et η > 0 nous avons
donc Z TZ
1 C
mes(A1 ) ≤ |∇un | ≤ ≤ε
k 0 Ω k
pour k assez grand, puisque (∇un ) est borné dans Lq (]0, T [×Ω) pour q <[(p − 1)(N +
1) + 1]/(N + 1) et donc dans L1 (]0, T [×Ω). Fixons donc k tel que mes(A1 ) ≤ ε et
mes(A2 ) ≤ ε (pour tous n, m ∈ N).
Majorons maintenant mes(A3 ), nous avons
Z TZ Z
|(un − um )| ≥ |(un − um )| ≥ η mes(A3 )
0 Ω A3
il suffit donc de montrer que (un ) est de Cauchy dans L1 (]0, T [×Ω). Comme (un ) est
borné dans Lq (]0, T [; W01,q (Ω)) pour q <[(p − 1)(N + 1) + 1]/(N + 1), div(a(t, x, ∇un ))
est borné dans Lq (]0, T [; W −1,q (Ω)) pour tout q <[(p − 1)(N + 1) + 1]/(p − 1)(N + 1).
De plus L1 (Ω) ⊂ W −1,q (Ω) pour tout q < N/(N − 1), et comme p > 2 − 1/(N + 1), on
a [(p − 1)(N + 1) + 1]/(p − 1)(N + 1) < N/(N − 1) donc (fn ) qui est, par hypothèse,
borné dans L1 (]0, T [×Ω), l’est aussi dans L1 (]0, T [; W −1,q (Ω)) pour q <[(p − 1)(N +
1) + 1]/(p − 1)(N + 1). Or unt = fn + div(a(t, x, ∇un )) donc (unt ) est borné dans
L1 (]0, T [; W −1,q (Ω)) pour tout q <[(p − 1)(N + 1) + 1]/(p − 1)(N + 1). Ainsi
donc, d’après un lemme de compacité du type de celui d’Aubin (voir, par exemple,
Simon [47]), (un ) converge (à une sous-suite près) dans Lq (]0, T [; Lq (Ω)) pour tout
q <[(p − 1)(N + 1) + 1]/(p − 1)(N + 1) et donc dans L1 (]0, T [×Ω). Ainsi la suite (un )
80 Chap. 8 Solutions entropiques de problèmes paraboliques
est convergente (à une sous-suite près) dans L1 (]0, T [×Ω) et est donc de Cauchy dans
L1 (]0, T [×Ω) donc pour η donné, il existe n0 tel que pour n, m ≥ n0 on ait
mes(A3 ) ≤ ε.
Il reste donc à majorer mes(A4 ), et à choisir η. D’après la monotonie de a, nous
avons [a(t, x, ξ1 ) − a(t, x, ξ2 )](ξ1 − ξ2 ) > 0 pour ξ1 − ξ2 6= 0 or l’ensemble des (ξ1 , ξ2 )
tels que |ξ1 | ≤ k, |ξ2 | ≤ k et |ξ1 − ξ2 | ≥ δ est compact et a est continue en ξ pour
presque tout t et x, donc [a(t, x, ξ1 ) − a(t, x, ξ2 )](ξ1 − ξ2 ) atteint sur ce compact son
minimum que nous noterons γ(t, x), qui vérifie γ(t, x) > 0 pp. De plus grâce à un
résultat d’intégration (voir chapitre 2), il existe ε0 > 0, tel que pour tout ensemble
mesurable A ⊂ ]0, T [×Ω
Z
(8.5) γ ≤ ε0 =⇒ mes(A) ≤ ε
A
γ et de A4 , nous avons
Z Z
γ≤ [a(t, x, ∇un ) − a(t, x, ∇um )](∇un − ∇um )1{|un −um |≤η} ,
A4 A4
de plus le terme intégré est positif et ∇Tη (un − um ) = (∇un − ∇um )1{|un −um |≤η} , nous
avons donc
Z Z TZ
γ≤ [a(t, x, ∇un ) − a(t, x, ∇um )]∇Tη (un − um )
A4 0 Ω
et si l’on choisit ϕ = Tη (un − um ) ∈ Lp (]0, T [; W 1,p (Ω)) ∩ L∞ (]0, T [×Ω) qui vérifie
0 0
Tη (un − um )t ∈ Lp (]0, T [; W −1,p (Ω)) dans l’équation (8.4) écrite avec un puis avec
um on obtient
Z T Z TZ
h(un − um )t , Tη (un − um )i + [a(t, x, ∇un ) − a(t, x, ∇um )]∇Tη (un − um )
0 0 Ω
Z TZ
= (fn − fm ) Tη (un − um )
0 Ω
le premier terme est positif (Θη (x) ≥ 0), et Θη (x) ≤ η|x| donc
Z TZ
[a(t, x, ∇un ) − a(t, x, ∇um )]∇Tη (un − um )
0 Ω
Z TZ Z
≤η |fn − fm | + η |un0 − um
0 | ≤ 2 η (kf kL1 + ku0 kL1 )
0 Ω Ω
8.2 La solution obtenue par approximation 81
Ainsi η étant fixé, la majoration de mes(A3 ) nous donne n0 tel que pour tout
n, m ≥ n0 on ait mes(A3 ) ≤ ε et donc
mes({|(∇un − ∇um )(x)| ≥ δ}) ≤ 4ε,
la convergence en mesure de (∇un ) vers ∇u est donc démontrée, ainsi que la propriété
(après extraction d’une sous-suite).
Ceci achève de prouver que u est solution de (8.2). Par la suite nous avons besoin
de la propriété suivante.
soit Z Z Z
[Θ1 (un − um )](t) ≤ Θ1 (un0 − um
0 ) +T |fn − fm |
Ω Ω Ω
or
Θ1 (x) ≤ |x|
donc Z Z Z
[Θ1 (un − um )](t) ≤ |un0 − um
0 | +T |fn − fm |
Ω Ω Ω
pour tout t ≤ T . Notons an,m le second membre, aussi
|un − um |(t)
Z Z Z
2
|un − um | (t) + ≤ [Θ1 (un − um )](t) ≤ an,m
|un −um | < 1 |un −um | > 1 2 Ω
et Z Z Z
|un − um |(t) = |un − um |(t) + |un − um |(t)
Ω |un −um | < 1 |un −um | > 1
Z 12
1
2
≤ |un − um | (t) mes(Ω) 2 + 2 an,m
|un −um | < 1
1 1
≤ (2 mes(Ω)) 2 an,m 2 + 2 an,m .
Comme (fn ) et (un ) convergent dans L1 on a an,m → 0 pour m et n → ∞ donc
(un ) est de Cauchy dans C([0, T ]; L1 (Ω)), aussi u ∈ C([0, T ]; L1 (Ω)) et ∀ t ≤ T on a
un (t) → u(t) dans L1 (Ω). Ce qui achève la démonstration de cette propriété.
si bien que ψ ∈ Lp (]0, T [; W01,p (Ω)). Ce choix de ψ dans (8.7) est donc possible et,
nous obtenons
Z T Z TZ Z TZ
n
hut , Tk (un − ϕ)i + a(t, x, ∇un )∇Tk (un − ϕ) = fn Tk (un − ϕ)
0 0 Ω 0 Ω
8.3 Existence d’une solution entropique 83
ce qui conduit à
Z Z Z T
[Θk (un − ϕ)](T ) − [Θk (un − ϕ)](0) + hϕt , Tk (un − ϕ)i
Ω Ω 0
Z TZ Z TZ
(8.8) + a(t, x, ∇un )∇Tk (un − ϕ) = fn Tk (un − ϕ)
0 Ω 0 Ω
Comme Tk+M (un ) est borné dans Lp (]0, T [, W01,p (Ω)), a(t, x, ∇Tk+M (un )) est borné
0
dans Lp (]0, T [×Ω) et donc converge faiblement. De plus ∇Tk+M (un ) → ∇Tk+M (u)
pp, donc a(t, x, ∇Tk+M (un )) → a(t, x, ∇Tk+M (u)) pp (car ∇Tk+M (un ) → ∇Tk+M (u)
0
presque partout) et dans Lp (]0, T [×Ω) faible. Et comme, grâce à la convergence
dominée, ∇ϕ 1{|Tk+M (un )−ϕ|≤k} converge dans Lp (]0, T [×Ω), on a
Z TZ
a(t, x, ∇Tk+M (un ))∇ϕ 1{|Tk+M (un )−ϕ|≤k}
0 Ω
Z TZ
→ a(t, x, ∇Tk+M (u))∇ϕ 1{|Tk+M (u)−ϕ|≤k}
0 Ω
l’autre terme étant positif, le lemme de Fatou nous donne (nous avons montré la
convergence presque partout de (∇un ) ce qui entraı̂ne celle de ∇Tk+M (un ))
Z TZ
a(t, x, ∇Tk+M (u))∇Tk+M (u)1{|Tk+M (u)−ϕ|≤k}
0 Ω
Z TZ
≤ lim a(t, x, ∇Tk+M (un ))∇Tk+M (un ) 1{|Tk+M (un )−ϕ|≤k}
0 Ω
Nous allons maintenant faire tendre successivement ε vers 0, n vers ∞ puis h vers
∞. Commençons par ε → 0. Notons Aε1 à Aε7 ces sept termes, nous obtenons alors
et donc
Aε1 + Aε2 + Aε4 + Aε5 ≤ Aε6 + Aε7 + |Aε3 |.
comme |Θk (v − Thε (un ))(T )| ≤ k[|v| + |Th (un )|](T ), (Thε )0 (x) ≤ Th0 (x) et |∇Tk (v −
Thε (un ))| ≤ [|∇Tk+h (v)| + |∇Th (un )], les quatre termes du membre de gauche et les
86 Chap. 8 Solutions entropiques de problèmes paraboliques
|un0 |1{|un0 |≥h−ε} et |(Thε )00 (x)| = (Rhε )00 (x) donc
R R
et Rε (un ) ≤
Ω h 0 Ω
Z TZ Z TZ
|(Thε )00 (un )| a(t, x, ∇un )∇un = (Rhε )00 (un ) a(t, x, ∇un )∇un
0 Ω 0 Ω
Z TZ Z
≤ |fn |1{|un |≥h−ε} + |un0 |1{|un0 |≥h−ε} .
0 Ω Ω
Donc Z TZ Z
|Aε3 | ≤k |fn |1{|un |≥h−ε} + |un0 |1{|un0 |≥h−ε}
0 Ω Ω
R TR R
qui converge vers k( 0 Ω |fn |1{|un |≥h} + Ω |un0 |1{|un0 |≥h} ) lorsque ε → 0, par convergence
dominée.
Il reste à observer, pour Aε4 , que
et nous avons montré que ∇un → ∇u presque partout donc ∇Th (un ) → ∇Th (u)
presque partout donc grâce au lemme de Fatou nous avons
Z TZ
[a(t, x, ∇v) − a(t, x, ∇Th (u))]∇Tk (v − Th (u))
0 Ω
Z TZ
≤ lim [a(t, x, ∇v) − a(t, x, ∇Th (un ))]∇Tk (v − Th (un ))
0 Ω
et les autres termes convergent (à une sous-suite près, pour avoir la domination de
(fn ) et (un0 ) (voir chapitre 2)) grâce au théorème de convergence dominée de Lebesgue
quand n → ∞ donc
Z Z
Θk (v − Th (u))(T ) − Θk (u0 − Th (u0 ))
Ω Ω
Z TZ
+ [a(t, x, ∇v) − a(t, x, ∇Th (u))]∇Tk (v − Th (u))
0 Ω
Z TZ Z TZ Z
0
≤ f (1 − (Th ) (u))Tk (v − Th (u)) + k |f |1{|u|≥h} + |u0 |1{|u0 |≥h}
0 Ω 0 Ω Ω
donc Z Z TZ
Θk (v − u)(T ) + [a(t, x, ∇v) − a(t, x, ∇v)]∇Tk (v − u) ≤ 0
Ω 0 Ω
89
9. Etude des singularités de la solution d’une
équation à second membre mesure
−∆w = δ
(9.1)
−∆v = f ∈ L1 (Ω)
vérifiant des conditions de Dirichlet homogène sur ∂Ω. On pourrait supposer qu’en
0, w l’emporte sur v (i.e. qu’elle explose plus au sens w − v → +∞ pour x → 0), il
n’en est rien comme le montre l’exemple suivant qui vérifie
Propriétés : Soit Ω = B1 (RN ), la boule unité de RN , soient w et v solutions de
(9.1) alors
w − v 6→ +∞ en 0, Tk (w − v) non continu en 0,
Z Z
1 1
lim (w − v) = +∞, lim Tk (w − v) = k,
ε→0 |Bε | Ω ε→0 |Bε | Ω
9.1 Construction de w et v
u0
00
− u + (N − 1) =g
r
d’où
(θ − 1)e
g=
rN (− ln(r/e))θ
92 Chap. 9 Etude des singularités de la solution . . .
d’où g ∈ L1 (Ω).
Soit ϕn radiale régulière telle que ϕn (r) = 1 pour r ≤ 1/4n2 et 0 pour r ≥ 1/3n2 .
Alors u ϕn est solution de
donc f ∈ L1 (Ω). De plus les supports des vn sont disjoints donc v(xn ) = ∞, et
v(yn ) = 0, ainsi v n’est pas continue en 0 (puisque xn et yn → 0 lorsque n → ∞).
9.2 Propriétés de w − v
−∆(w − v) = δ − f
En effet Z Z Z
1 1 1
(w − v) = w− v
|Bε | Bε |Bε | Bε |Bε | Bε
et d’autre part Z X 1 Z X 1 Z
v= 2
vn = 2
vn
Bε n≥1
n B ε
n B ε n≥1/ε
donc
Z Z
1 1 X 1
v≤ u .
|Bε | Bε |Bε | Bε n2
n≥1/ε
Rr Rr
Comme u(r) = u0 + u(1) =
1 1
u0 (car u(1) = 0) et (− ln(r/e)) ≥ 1 pour r ≤ 1
donc Z r r
−1
Z
0 1 1
u(r) = u ≤ dr = −1
1 1 rN −1 N − 2 rN −2
ainsi
1
u≤ w
(N − 2)CN
donc Z
1 1 1
u≤ −1
|Bε | Bε N −2 2εN −2
si bien que
Z X
1 1 1 1 ε 1
v≤ −1 ≤ −1
|Bε | Bε N −2 2εN −2 n2 N − 2 2εN −2
n≥1/ε
94 Chap. 9 Etude des singularités de la solution . . .
Comme w − v ne peut être négative que dans le support de v qui est la réunion des
translatés des B1/3n2 (c’est-à-dire le support des vn ) et qu’en dehors, w − v ≥ k si ε
est assez petit, on a
Z Z Z
Tk (w − v) ≥ S
k+ S −k
Bε Bε − n≥1/ε (B1/3n2 +xn ) n≥1/ε (B1/3n2 +xn )
X X 1 N DN ε2N −1
|B1/3n2 | = DN ≤
3n2 3N (2N − 1)
n≥1/ε n≥1/ε
donc Z
Tk (w − v) ≥ k DN (εN − Cε2N −1 ) ≥ k DN εN (1 − CεN −1 )
Bε
d’où Z
1
Tk (w − v) ≥ k(1 − CεN −1 )
|Bε | Bε
et Tk (w − v) ≤ k donc Z
1
Tk (w − v) ≤ k
|Bε | Bε
et finalement Z
1
lim Tk (w − v) = k.
ε→0 |Bε | Bε
Les espaces de Lebesgue ou de Sobolev sont insuffisants pour distinguer les solutions
des équations elliptiques à second membre L1 et mesure : en effet, il existe u et v
solutions de
(10.1) −∆u = δ
(10.2) −∆v = f ∈ L1
N
dans Ω la boule unité de RT , satisfaisant des conditions homogènes de Dirichlet sur
1,q 1,N/(N −1)
∂Ω, et telles que u et v ∈ q < N W0 (Ω), mais u et v 6∈ W0 (Ω). La solution
N −1
u est la fonction (bien connue) de Green, et la fonction v est due à Orsina [40]. Ces
deux fonctions appartiennent aussi à un même espace de Macinkiewicz.
Comme ces espaces ne permettent pas de distinguer ces deux fonctions, il pourrait
être intéressant de considérer une famille plus fine d’espaces et en particulier les
espaces de Lorentz. Nous allons montrer que u et v n’appartiennent pas aux mêmes
espaces, mais que l’on peut construire w solution de
(10.3) −∆w = g ∈ L1
appartenant aux mêmes espaces que u. Les espaces de Lorentz sont donc, eux aussi,
insuffisants.
Nous utilisons la définition de Ziemer [50] des espaces de Lorentz, pour définir f ∗
nous reprenons aussi la présentation de Kavian [31].
où d’une part, mesN {f > t} = mesN {x; f (x) > t} et mesN désigne la mesure de
Lebesgue dans Ω, et d’autre part, mes1 {f ∗ > t} = mes1 {s; f ∗ (s) > t} et mes1 est
la mesure de Lebesgue sur ]0, |Ω|[. On prolonge f ∗ par 0 en dehors de ]0, |Ω|[. Cette
fonction est appelée rearrangement décroissant de f . On peut choisir
c’est-à-dire que si mes{f > s} est bijective, f ∗ (s) = (mesN {f > s})−1 . Grâce au
théorème de Fubini,
Z Z Z ∞ Z ∞Z Z ∞
p
f dx = 1{f p > t} dt dx = 1{f > t1/p } dx dt = mesN {f > t1/p } dt
Ω Ω 0 0 Ω 0
et de même Z ∞ Z ∞
∗ p
(f ) dt = mes1 {f ∗ > t1/p } dt,
0 0
∗
et comme mesN {f > t} = mes1 {f > t} on a finalement
Z Z ∞
p
f dx = (f ∗ )p dt.
Ω 0
En particulier, et selon le calcul précédent, pour 1 < p < +∞ on a Lp,p (Ω) = Lp (Ω).
Rappelons maintenant la définition des espaces de Marcinkiewicz (voir [5])
On a M p (Ω) = Lp,∞ (Ω). En effet notons ϕf (t) = mesN {|f | > t}, alors
ϕf (|f |∗ (t)) ≤ t et |f |∗ (ϕf (t)) ≤ t pour t ≥ 0. On déduit donc aisément que si
|f |∗ (t) ≤ Ct−1/p alors ϕf (s) ≤ C p /sp donc Lp,∞ (Ω) ⊂ M p (Ω), et si ϕf (s) ≤ C/sp
alors |f |∗ (t) ≤ C 1/p t−1/p donc M p (Ω) ⊂ Lp,∞ (Ω).
10.2 Quelques fonctions 97
on peut vérifier aisément que ∇u ∈ Lq (Ω) pour tout q < N/(N − 1) mais que
∇u 6∈ LN/N −1 (Ω). Déterminons à quel espace de Lorentz appartient ∇u. On a
|∇u|(r) = (N − 2)CN /rN −1 donc pour t ≤ |Ω|
C
|∇u|∗ (t) = |∇u|((t/DN )1/N ) = N −1
t N
donc
Z ∞ Z |Ω| q Z |Ω|
dt 1 C dt q q(N −1)
1/p
(t ∗
|∇u| (t))q
= t p
N −1 =C tp− N
−1
dt
0 t 0 t N t 0
qui est fini si et seulement si q/p − q(N − 1)/N − 1 > −1 et donc si p < N/(N − 1)
et q quelconque. De plus t1/p |∇u|∗ (t) = Ct1/p−(N −1)/N qui est borné au voisinage de
0 pour p ≤ N/(N − 1). Ainsi
N
|∇u| ∈ Lp,q (Ω), 1<p< , 1 < q < +∞,
N −1
N
|∇u| ∈ Lp,∞ (Ω), 1<p ≤ ,
N −1
mais N
|∇u| 6∈ L N −1 ,q (Ω), 1 < q < +∞.
Étudions maintenant la solution de (10.2) qui s’écrit, pour v radiale,
v0
00
− v + (N − 1) = f,
r
C
|∇v|∗ (t) = N −1
1
t N (1 − N
ln( DtN ))θ−1
donc d’après des calculs analogues aux précédents, on a, pour 1 < p < N/(N − 1) et
1 < q < ∞,
Z |Ω|
dt
(t1/p |∇v|∗ (t))q < +∞,
0 t
R |Ω| 1/p
et 0 (t |∇v|∗ (t))q /t dt = +∞ si p > N/(N − 1). Mais pour p = N/(N − 1),
Z |Ω| Z |Ω|
1/p ∗ dtq C dt
(t |∇v| (t)) = 1 < +∞
0 t 0 t(1 − N
ln( DtN ))q(θ−1)
si q(θ − 1) > 1, or θ − 1 <(N − 1)/N donc l’intégrale est finie pour N/(N −
1) < 1/(θ − 1) < q < +∞ et p = N/(N − 1). Enfin t1/p |∇v|∗ (t) = Ct1/p−(N −1)/N /(1 −
ln(t/DN )/N )q(θ−1) qui est borné au voisinage de 0 pour p ≤ N/(N − 1), donc
∇v ∈ Lp,∞ (Ω). Ainsi v vérifie, comme u
N
|∇v| ∈ Lp,q (Ω), 1<p< , 1 < q < +∞,
N −1
N
|∇v| ∈ Lp,∞ (Ω), 1<p ≤ .
N −1
Mais on a, de plus,
N N 1
|∇v| ∈ L N −1 ,q (Ω), < < q < +∞.
N − 1 (θ − 1)
w0
00
− w + (N − 1) =g
r
d’où
1
w0 (r) =
rN −1 (ln(2 − ln r))1+α
donc
C
|∇w|∗ (t) = N −1
1
t N [ln(2 − N
ln( DtN ))]1+α
10.2 Quelques fonctions 99
pour tout q < +∞. Enfin t1/p |∇w|∗ (t) = t1/p−(N −1)/N /[ln(2 − N1 ln(t/DN ))]1+α qui est
fini au voisinage de 0 pour p ≤ N/(N − 1). Ainsi
N
|∇w| ∈ Lp,q (Ω), 1<p< , 1 < q < +∞,
N −1
N
|∇w| ∈ Lp,∞ (Ω), 1<p ≤ ,
N −1
mais N
|∇w| 6∈ L N −1 ,q (Ω), 1 < q < +∞.
Donc u et w appartient aux même espaces de Lorentz, qui sont donc insuffisants
pour distinguer les solutions à second membre L1 et δ.
11. Classes de fonctions test pour problèmes
elliptiques à second membre mesure
et
pour tout w ∈ W 1,p (Ω) ∩ L∞ (Ω) tel que ∇w = 0 sur {x; |u(x)| ≥ m} pour un m ≥ 0.
Comme on peut choisir w = Tk (u−ϕ) ci-dessus, cette dernière formulation conduit
aussi à l’unicité. L’équivalence de ces formulations dit, en quelque sorte, que Tk (u−ϕ)
et S(u)ϕ sont (( denses )) dans les w.
11.2 Équivalence
Montrons (d’après [38]) que les solutions entropiques sont solutions w. Soit
w ∈ W01,p (Ω) ∩ L∞ (Ω) tel qu’il existe m tel que ∇w = 0 sur {x; |u(x)| ≥ m}. Alors
soit ϕ = Tn (u) − w, on a bien ϕ ∈ W01,p (Ω) ∩ L∞ (Ω) que l’on peut utiliser dans (11.2),
ce qui donne
Z Z
a(x, ∇u)∇Tk (u − Tn (u) + w) ≤ f Tk (u − Tn (u) + w)
Ω Ω
et donc
Z Z Z
a(x, ∇u)∇(u−Tn (u))+ a(x, ∇u)∇w ≤ f Tk (u−Tn (u)+w)
|u−Tn (u)+w|≤k |u−Tn (u)+w|≤k Ω
et si k ≥ kwk∞ on a donc
Z Z
a(x, ∇u)∇w → a(x, ∇Tm (u))∇w
|u−Tn (u)+w|≤k Ω
R
Reste donc le premier terme |u−Tn (u)+w|≤k a(x, ∇u)∇(u − Tn (u)) mais si |u| ≤ n alors
u − Tn (u) = 0, et si |u| ≥ n alors ∇Tn (u) = 0 donc
Z Z
a(x, ∇u)∇(u − Tn (u)) = a(x, ∇u)∇u 1|u|≥n
|u−Tn (u)+w|≤k |u−Tn (u)+w|≤k
donc Z Z
(11.4) a(x, ∇Tm (u))∇w ≤ f Tk (w)
Ω Ω
Bien entendu le choix de w = Tk (u−ϕ) dans (11.4) conduit à (11.2) donc l’unicité
des solutions entropiques nous donne celles des solutions de (11.4).
104 Chap. 11 Classes de fonctions test pour problèmes elliptiques. . .
Nous allons maintenant montrer que les solutions renormalisées sont solutions de
(11.4) (voir [39, 16] avec des modifications mineures).
Choisissons ϕ = w et S = Sn dans (11.3), avec Sn régulière telle que 0 ≤ Sn ≤ 1,
S(x) = 0 si |x| ≥ n + 1, S(x) = 1 si |x| ≤ n, et radiale linéaire par morceaux. Alors
Z Z Z
0
Sn (u)a(x, ∇u)∇w + Sn (u)w a(x, ∇u)∇u ≤ f Sn (u)w
Ω Ω Ω
et Z Z
w a(x, ∇u)∇u ≤ Ckwk∞ |∇u|p
n≤|u|≤n+1 n≤|u|≤n+1
donc Z
Sn 0 (u)w a(x, ∇u)∇u → 0
Ω
et finalement Z Z
a(x, ∇u)∇w = f w.
Ω Ω
11.3 Non unicité 105
Nous allons montrer que, bien que contenant plus de fonctions test que les
deux autres formulations, cette troisième formulation ne conduit pas à l’unicité des
solutions quand f est une mesure.
En fait Boccardo, Gallouët et Orsina [14] ont montré qu’il y a unicité si f ne charge
pas les ensembles de capacité nulle, c’est-à dire, en particulier, en dimension N ≥ 3, si
f ne charge pas des objets de dimension 0 ou 1 : les points, les segments. . . Autrement
dit, pour montrer que la solution n’est pas unique on peut choisir f chargeant des
points ou des segments, c’est ce que nous allons faire. En dimension 2, il faudrait
que f ne charge que des points et le contre exemple ci dessous ne marcherait pas
(les limites en +∞ seraient finies), en fait il y a unicité : voir [29]. En dimension 1,
M (Ω) ⊂ H −1 (Ω) il y a donc unicité.
Soit Ω = BRN , pour N ≥ 3 la boule unité de RN . Et soit
A = {−1/2 ≤ x1 ≤ 1/2, x2 = · · · = xn = 0}
1
G(x, y) = CN + H(x, y)
|y − x|N −2
N −2
Soit f (z) la primitive de 1/(1 + z 2 ) 2 nulle en 0, qui est donc impaire, alors
1 1
− x1
1 1 2 2
+ x1
B(x) = [f (z1 ) − f (z2 )] = N −3 f +f
ρN −3 ρ ρ ρ
Ainsi, il existe ε0 tel que |u(x)| ≥ k sur {x; x1 ∈ [−1/2, 1/2] et ρ2 ≤ ε0 } ∪ {x; x1 ≥
1/2 et ρ2 + (x1 − 1/2)2 ≤ ε0 } ∪ {x; x1 ≤ −1/2 et ρ2 + (x1 + 1/2)2 ≤ ε0 }. Soit Aε0 cet
ensemble, et soit w tel que ∇w = 0 là où |u| ≥ k, donc ∇w = 0 sur Aε0 . Soit ε00 < ε0
on a bien sûr Aε00 ⊂ Aε0 et même d(Aε00 , Aε0 c ) > 0, on peut donc construire une suite
wn ∈ W01,r (Ω) pour r > N , telle que wn → w dans H01 (Ω), ∇wn = 0 sur Aε00 pour tout
n et wn (0) → w(0).
Comme wn est régulière, ∇wn = 0 sur Aε00 , et wn = const = wn (0) sur Aε00 ⊃ A,
Z Z Z
∇u ∇wn = ∇u ∇wn = wn dδA = wn (0)
Aε00 c Ω Ω
et il existe h tel que |u(x)| ≤ h sur Aε00 c (voir ci-dessus avec ε4 = ε00 ), donc Th (u) = u
sur Aε00 c donc
Z Z Z Z
∇u ∇wn = ∇u ∇wn = ∇Th (u)∇wn = ∇Th (u)∇wn
Ω Aε00 c Aε00 c Ω
lorsque n → +∞ (car |u| ≤ h sur Aε0 c ⊂ Aε00 c ). Or wn (0) → w(0) donc pour tout
w ∈ H01 (Ω) ∩ L∞ (Ω) tel que ∇w = 0 lorsque |u| ≥ k
Z Z
∇u ∇w = w(0) = w dδ.
Ω Ω
Ainsi u est solution w de −∆u = δ alors que −∆u = δA au sens des distributions.
−∆v = δ
bien-sûr
v(x) = G(x, 0) = 1/|x|N −2 + H(x, 0)
108 Chap. 11 Classes de fonctions test pour problèmes elliptiques. . .
donc pour k donné, il existe ε5 tel que |v(x)| ≥ k lorsque |x| ≤ ε5 et pour ε6 donné,
il existe h tel que |v(x)| ≤ h lorsque |x| ≥ ε6 , aussi, comme pour u, on peut montrer
que pour tout w ∈ H01 (Ω) ∩ L∞ (Ω) tel que ∇w = 0 lorsque |v| ≥ k
Z Z
∇v ∇w = w(0) = w dδ
Ω Ω
donc v est solution w de −∆v = δ ainsi que −∆v = δ au sens des distributions.
Il y a donc deux solutions différentes pour la même équation, la formulation est
donc insuffisante.
12. Convergence faible des tronqués des solutions
problèmes elliptiques à second membre mesure
Boccardo et Gallouët [11, 12] ont montré que, pour Ω ouvert borné de RN , le problème
elliptique de Dirichlet homogène
admet une solution pour f ∈ M (Ω) = (C(Ω))0 . Cette solution est obtenue par
approximation, c’est-à-dire qu’elle est la limite de (un ), où un est solution de
u0
00
− u + (N − 1) =f
r
avec r = |x| norme euclidienne dans RN . La solution générale sans second membre
est 1/rN −2 et la solution particulière associée à 1 est −r2 /(2N ).
110 Chap. 12 Convergence faible des tronqués. . .
Soient 1 1
1 2N 3N
an = , bn = , cn = ,
n n n
1 1 1
4N 5N 6N
dn = , en = , fn = ,
n n n
si bien que
1
fn N − en N = en N − dn N = dn N − cn N = cn N − bn N = bn N − an N = .
nN
Et soient hn = nN et gn tels que
hn sur 0 ≤ r ≤ an ,
0 sur an ≤ r ≤ cn ,
gn =
− hn sur cn ≤ r ≤ en ,
hn sur en ≤ r ≤ fn .
dès que Z
gn Tk (un ) → 0
Ω
ce qui est bien le cas par convergence dominée si gn → 0 dans L1 (Ω). Ceci montre
notamment que Z
|∇Tk (un )|2 → 0
Ω
1
dans le cas g ∈ L (Ω).
R Effectivement dans l’exemple que nous venons de construire on peut montrer que
g T (u ) ne tend pas vers 0 (ici gn converge seulement dans M (Ω) faible ∗, aussi
Ω n k n
le théorème de convergence dominée ne permet pas de conclure) :
Z Z in Z fn
gn Tk (un ) = gn Tk (un ) + gn Tk (un )
Ω 0 in
Z an Z en Z fn
= hn Tk (un ) − Tk (un ) + Tk (un )
0 cn en
112 Chap. 12 Convergence faible des tronqués. . .
R en R fn
or hn cn
1 = 2 CN et 0 ≥ hn en
Tk (un ) ≥ −k CN donc
Z
gn Tk (un ) ≥ 3 k CN − k CN = 2 k CN
Ω
Ce chapitre est une partie d’un article en préparation écrit en collaboration avec G.
Dal Maso, F. Murat et L. Orsina [22].
13.1.1 Assumptions
for almost every x in Ω and for every ξ in RN , where β > 0 is given, and b is a non
negative function in Lp (Ω) ;
for almost every x in Ω and for every ξ, ξ 0 in RN , ξ 6= ξ 0 . Observe that (13.1.2), and
the continuity of a with respect to ξ, imply
(see for example [33] and [34]). Here and in the following, we will denote by h·, ·i
0
the duality between W −1,p (Ω) and W01,p (Ω). Furthermore, thanks to (13.1.4), such a
solution is unique.
We define Mb (Ω) as the set of measures with bounded total variation. If p > N ,
0
then Mb (Ω) is contained in W −1,p (Ω), so that both existence and uniqueness of
solutions for (13.0.1) follow from the result quoted above. This explains the restriction
p ≤ N that we have imposed.
13.1 Hypothèses et énoncés des résultats 115
In Section 3, Propositions 13.2.4 and 13.2.5, we will recall a result of [14] that
states that every measure µ in Mb (Ω) can be decomposed as follows :
µ = µ0 + λ+ − λ− = f − div (g) + λ+ − λ− , (13.1.7)
where µ0 is a measure in M0 (Ω), which is the set of measures of Mb (Ω) that do
not charge the subsets of Ω of zero p-capacity (see Section 3 for the definition of
p-capacity), and so can be decomposed as f − div (g), with f in L1 (Ω) and g in
0
(Lp (Ω))N (see Proposition 13.2.5), while λ+ and λ− are two non negative measures
in Mb (Ω) which are concentrated on two subsets E + and E − , respectively, of zero
p-capacity.
We define Tk (s) = max(−k, min(k, s)) for s in R, and for k a positive real number
(see (13.1.26) in Subsection 13.1.5).
(a) the function u is almost everywhere finite, and is such that Tk (u) belongs to
W01,p (Ω) for every k > 0 ;
(b) ∇u, defined in Definition 13.1.1 is such that
|∇u|p−1 belongs to Lq (Ω), for every q < NN−1 ; (13.1.8)
(c) for every function w in W01,p (Ω) ∩ L∞ (Ω) such that there exists k > 0 and w+∞
and w−∞ in W 1,r (Ω), with r > N , such that
(
w = w+∞ capp -quasi everywhere on the set {u > k},
(13.1.9)
w = w−∞ capp -quasi everywhere on the set {u < −k},
116 Chap. 13 Une généralisation des notions de solutions . . .
we have
Z Z Z Z
a(x, ∇u) · ∇w dx = w dµ0 + w dλ − w−∞ dλ− .
+∞ +
(13.1.10)
Ω Ω Ω Ω
Remark 13.1.3 Every term in (13.1.10) has a meaning. Indeed, the integral on the
left hand side can be splitted as
Z Z Z
a(x, ∇u) · ∇w dx + a(x, ∇u) · ∇w dx + a(x, ∇u) · ∇w dx ,
{u < −k} {|u|≤k} {u > k}
where all three terms are well defined : actually, since |∇u|p−1 belongs to Lq (Ω) for
every q < NN−1 , hypothesis (13.1.3) implies that a(x, ∇u) belongs to (Lq (Ω))N , for
every q < NN−1 ; on the other hand
Z Z
a(x, ∇u) · ∇w dx = a(x, ∇u) · ∇w−∞ dx < +∞ ,
{u < −k} {u < −k}
since w−∞ belongs to W 1,r (Ω) with r > N and a(x, ∇u) belongs to (Lq (Ω))N for every
q < NN−1 = N 0 ; the same thing is true for the term
Z
a(x, ∇u) · ∇w dx .
{u > k}
Finally, for the middle term we have, since a(x, 0) = 0 (see (13.1.5)),
Z Z
a(x, ∇u) · ∇w dx = a(x, ∇Tk (u)) · ∇w dx < +∞ ,
{|u|≤k} Ω
since Tk (u) belongs to W01,p (Ω), which implies by (13.1.3) that a(x, ∇Tk (u)) belongs
0
to (Lp (Ω))N , while w belongs to W01,p (Ω).
As for the right hand side, the term
Z Z
w dλ − w−∞ dλ−
+∞ +
Ω Ω
is obviously well defined, since both w+∞ and w−∞ are continuous and bounded on
Ω, while the term Z
w dµ0
Ω
is well defined since w belongs to W01,p (Ω) ∩ L∞ (Ω) and so to L∞ (Ω, dµ0 ) (see
Proposition 13.2.3), and therefore to L1 (Ω, dµ0 ) since µ0 belongs to M0 (Ω).
Example 13.1.4 Observe that we did not assume that the function u belongs to
some space Lr (Ω), but only that u is Lebesgue measurable. Indeed it is possible that
the function u does not belong to L1loc (Ω) : this is the case in the following example.
13.1 Hypothèses et énoncés des résultats 117
Note that u belongs to L1loc (Ω) if and only if γ < N , i.e., if p > N2N
+1
.
Let us show that u is a renormalized solution of the equation
−div (|∇u|p−2 ∇u) = ωN δ0 in Ω,
(13.1.12)
u=0 on ∂Ω,
where δ0 is the Dirac mass concentrated at the origin. Indeed, Tk (u) belongs to
W01,p (Ω) (and is actually Lipschitz continuous and zero on the boundary of Ω), and,
defining ∇u as in Definition 13.1.1, we have
x
∇u = − for every 1 < p ≤ N ,
|x|γ+2
which implies that |∇u|p−1 = |x|1−N belongs to Lq (Ω) for every q < NN−1 . Thus, (i) is
satisfied.
For what concerns (c), consider a function w in W01,p (Ω) ∩ L∞ (Ω) which belongs
to C 1 in a neighbourhood of the origin (so to meet the requirements of Definition
13.1.2). Using the fact that x·∇w
|x|N
belongs to L1 (Ω) (see Remark 13.1.3), integrating by
parts on B1 (0)\Bε (0), using the fact that div ( |x|xN ) = 0 in the sense of distributions
on RN \{0}, and finally using the continuity of w at the origin, we have
x · ∇w
Z Z
p−2
|∇u| ∇u · ∇w dx = lim+ N
dx
Ω ε→0
ZB1 (0)\Bε (0) |x|
x x
= lim+ w N · dσ
ε→0 ∂Bε (0)Z |x| |x|
1
= lim+ N −1 w dσ = ωN w(0)
ε→0 ε ∂Bε (0)
Z
= ωN w+∞ dδ0 .
Ω
2N
In this example u does not belong to L1 (Ω) if p ≤ N +1
, and ∇u is not the
distributional gradient of u.
Remark 13.1.5 Note that a renormalized solution u, which is a priori defined only
almost everywhere in Ω, is actually defined capp -quasi everywhere (that is to say,
except on a set of zero p-capacity), since Tk (u) belongs to W01,p (Ω) for every k > 0,
and since the functions of W01,p (Ω) admit a representative which is defined capp -quasi
everywhere (see Section 3).
118 Chap. 13 Une généralisation des notions de solutions . . .
Remark 13.1.6 Since the choice of test functions in (c) of Definition 13.1.2 is rather
complicated, one may wonder whether the class of admissible test functions is not
empty. This is indeed the case, since for example any function in Cc∞ (Ω) is admissible.
Thus, if u is a renormalized solution of (13.0.1) in the sense of Definition 13.1.2, and
if p > 2 − N1 , then it is also a solution in the sense of distributions.
There are more admissible functions, built after u, such as Tk (u) (choosing
w +∞
≡ k and w−∞ ≡ −k).
If ϕ belongs to W01,p (Ω) ∩ L∞ (Ω), then it is possible to choose in (13.1.10) the
function w = Tk (u − ϕ) ; indeed, this function belongs to W01,p (Ω) ∩ L∞ (Ω), and
we can choose w+∞ ≡ k and w−∞ ≡ −k, on the sets {u > k + kϕkL∞ (Ω) } and
{u < −k − kϕkL∞ (Ω) } respectively. Thus, a renormalized solution of (13.0.1) turns
out to be an entropy solution of (13.0.1) in the sense defined in [4]. Hence, if the
measure µ does not charge the sets of zero p-capacity, there exists at most one
renormalized solution of (13.0.1), due to the uniqueness result of [14]. Note however
that the definition of entropy solution with datum a Dirac mass given in [14], Remark
3.4, did not imply that an entropy solution is a distributional solutions. In contrast,
we have seen that a renormalized solution as defined in Definition 13.1.2 is also a
distributional solution, which rules out the counterexample to uniqueness given in
[14], Remark 3.4.
and Z Z
1
lim a(x, ∇u) · ∇u ϕ dx = ϕ dλ− . (13.1.14)
n→+∞ n {−2n < u≤−n} Ω
(e) For every h in W 1,∞ (R) with compact support in R, and for every ϕ in
W01,p (Ω) ∩ L∞ (Ω), we have
Z Z Z
0
a(x, ∇u) · ∇u h (u) ϕ dx + a(x, ∇u) · ∇ϕ h(u) dx = h(u) ϕ dµ0 .
Ω Ω Ω
(13.1.15)
13.1 Hypothèses et énoncés des résultats 119
Remark 13.1.8 Observe that every term in (13.1.15) has a meaning : indeed, in the
right hand side h(u) ϕ belongs to W01,p (Ω) ∩ L∞ (Ω), and is hence in L∞ (Ω, dµ0 ) by
Proposition 13.2.3. On the other hand, since supp (h) ⊆ [−M, M ] for some M > 0,
the left hand side as to be understood as
Z Z
0
a(x, ∇TM (u)) · ∇TM (u) h (u) ϕ dx + a(x, ∇TM (u)) · ∇ϕ h(u) dx ,
Ω Ω
where both terms are finite by (13.1.3) since both ϕ and TM (u) belong to W01,p (Ω) ∩
L∞ (Ω).
Remark 13.1.9 Definition 13.1.7 is similar to the0 definition of renormalized solution
given in [36] or [39] if µ belongs to L1 (Ω) + W −1,p (Ω) (that is to say, after the result
of [14] cited in Proposition 13.2.4 below, to M0 (Ω)). Indeed in these papers, the
definition of renormalized solution included (e) as well as
Z
1
lim a(x, ∇u) · ∇u dx = 0 , (13.1.16)
n→+∞ n {n≤|u| < 2n}
which coincides with (13.1.13) and (13.1.14) in the case λ+ = λ− = 0 and ϕ ≡ 1. Here
(13.1.13) and (13.1.14) replace (13.1.16), and specify the behaviour of the energy of
u on the set where u is very large.
Conditions (13.1.13) and (13.1.14) can be removed if we change the class of
admissible test functions.
Definition 13.1.10 Assume that a satisfies (13.1.2)–(13.1.4), and let µ be a measure
in Mb (Ω), which is decomposed as µ0 + λ+ − λ− . A measurable function u is a
renormalized solution of (13.0.1) if u satisfies (a) and (b) of Definition 13.1.2, and if
the following holds :
(f) for every h in W 1,∞ (R) such that h0 has compact support in R, and for every
ϕ in W01,r (Ω), with r > N ,
Z Z
0
a(x, ∇u) · ∇u h (u) ϕ dx + a(x, ∇u) · ∇ϕ h(u) dx
Ω Z Z Ω Z (13.1.17)
+∞ + −∞ −
= h(u) ϕ dµ0 + h ϕ dλ − h ϕ dλ .
Ω Ω Ω
+∞ −∞
Here h and h are the limits of h(s) at +∞ and −∞ respectively (which
exist, since h is constant for |s| large enough).
Remark 13.1.11 As in (13.1.15), every term in (13.1.17) is well defined : this is clear
for the right hand side since h(u) ϕ belongs to L∞ (Ω, dµ0 ) (see Proposition 13.2.3),
and thus to L1 (Ω, dµ0 ). Since supp (h0 ) ⊆ [−M1 , M1 ], for some M1 , the left hand side
has to be understood as
Z Z
0
a(x, ∇TM1 (u)) · ∇TM1 (u) h (u) ϕ dx + a(x, ∇u) · ∇ϕ h(u) dx ,
Ω Ω
where both terms are finite : the first one by (13.1.3) since TM1 (u) belongs to W01,p (Ω),
and the second one since a(x, ∇u) belongs to (Lq (Ω))N for every q < NN−1 due to the
hypotheses on |∇u|p−1 and to (13.1.3).
120 Chap. 13 Une généralisation des notions de solutions . . .
Remark 13.1.12 The conditions (13.1.13) and (13.1.14) that were required to hold
in Definition 13.1.7 are now “embedded” in the formulation (13.1.17) : see Step 3 of
the proof of Theorem 13.1.16 in Section 6, where these conditions will be obtained
from (13.1.17).
Remark 13.1.13 We remark that if the function h in both Definition 13.1.7 and
Definition 13.1.10 is such that h(0) = 0, then the test functions ϕ can be respectively
chosen in W 1,p (Ω) ∩ L∞ (Ω) and in W 1,r (Ω), with r > N .
λ+ +
k → λ , λ− −
k → λ in the weak∗ topology of measures.
As usual, every term in (13.1.19) is well defined due to the regularity of Tk (u) and
0
ϕ, and since measures in M0 (Ω) are in W −1,p (Ω) + L1 (Ω) (see Proposition 13.2.5).
Remark 13.1.15 Some comments on these new definitions are in order. Definition
13.1.7 and 13.1.10 state (we refer explicitly to Definition 13.1.10) that the renormal-
ized solution u is in some sense equal to +∞ on the sets charged by λ+ , and to −∞
on the sets charged by λ− ; this is clearly expressed by the presence of the two terms
h+∞ and h−∞ in (13.1.17). Thus, the presence of a measure which is concentrated
on a set of zero p-capacity leads to an unbounded solution that we can suppose to
be unbounded “in a stronger way” than the solutions corresponding to measures in
M0 (Ω) : this fact is very clear if one considers solutions corresponding to a Dirac
mass and to a function in L1 (Ω).
Definition 13.1.14 says something about the equation solved by Tk (u). Indeed,
(13.1.19) can be reformulated (in the sense of distributions, that is, choosing ϕ in
Cc∞ (Ω)) as follows :
−
−div (a(x, ∇Tk (u))) = µ0 {|u| ≤ k} + λ+
k − λk .
13.1 Hypothèses et énoncés des résultats 121
This formulation explains also the required conditions (13.1.18) and the convergence
−
of λ+k and λk . Observe that a renormalized solution in the sense of Definition 13.1.14
−
is a solution obtained by approximation (since the sequence µ0 {|u| ≤ k} + λ+ k − λk
converges to µ = µ0 + λ+ − λ− in the weak∗ topology of measures) with the property
that the truncates are strongly convergent (it is indeed easy to see that Th (Tk (u))
converges strongly to Th (u) in W01,p (Ω) as k tends to infinity).
The equivalence of the four definitions will be proved in Section 6.
in the sense specified by (13.1.6). We then have the following result, which implies
Theorem 13.1.16.
122 Chap. 13 Une généralisation des notions de solutions . . .
Remark 13.1.18 We explicitly remark that we require that the two sequences (λ+ )ε
and (λ− )ε are sequences of non negative functions. If we do not make this requirement,
the result of Theorem 13.1.17 (as well as that of Theorem 13.1.20) may not hold, see
examples 13.5.2 and 13.5.3 in Section 13.5.
Remark 13.1.19 A measure in Mb (Ω) is not the not most general possible datum
0
for (13.0.1). Indeed, there exist elements in W −1,p (Ω) which are not measures. This
implies that a more general datum is of the kind
µ − div (F ) ,
0
with µ in Mb (Ω) and F in (Lp (Ω))N . However, the new term −div (F ) does not give
any problem in the proof of the results, since it can be treated as the term −div (gε ).
Thus, we will restrict ourselves to the case of a datum µ belonging to Mb (Ω), the
case of a datum µ − div (F ) being analogous.
Let us just explicitly state that every of the present paper holds if µ in Mb (Ω)
is replaced by
0
µ − div (F ) with µ in Mb (Ω) and F in (Lp (Ω))N .
The main tool of the proof of Theorem 13.1.16 will be the following strong
convergence result for the truncates of uε , which is interesting on its own, and specifies
the sense in which uε converges to u.
This subsection contains some notation, the definition of some of the “objects”
that we will use in the proof of the results, and some already known results on the
sequence uε of solutions of (13.1.24).
We will use the following functions of one real variable, which may depend on
one or more non negative real parameters such as k and n :
−k
-
k s
−k
0 if s < −n − k,
x+n+k
if −n − k ≤ s < −n,
k
Hn,k (s) = 1 if |s| ≤ n, (13.1.27)
n+k−x
if n < s ≤ n + k,
k
0 if s > n + k.
Hn,k (s)
6
1
Q
Q
Q -
−n−k −n n n+k s
We will also use the function k − Tk (s) (and its companion k + Tk (s)) :
k−Tk (s)
6
2k
@
@
@
@
@ -
−k k s
In the whole present paper we will denote by c a generic constant, which can
vary from line to line, but which depends on N , p, Ω, α, β and b (which appear in
Subsection 2.1), and never on another parameter (such as ε, δ, η, n).
Moreover, if η, δ and ε are positive real numbers, and n belongs to N, we will
denote by ω(η, δ, n, ε) any quantity such that
If the order in which the limits are taken will be different, we will change the order of
appearance of the symbols, from the last limit to be taken, to the first : for example,
ω(η, n, δ, ε) is any quantity whose absolute value converges to zero after taking the
limits on ε, δ, n and η successively. If the quantity we consider does not depend on one
among η, δ, n and ε, we will omit the dependence from the corresponding variable :
as an example, ω(η, ε) is a quantity such that
Finally, we will denote (for example) by ωη,δ (n, ε) a quantity that depends on η, δ,
n, ε and is such that
lim lim+ |ωη,δ (n, ε)| = 0 ,
n→+∞ ε→0
Let now µε = fε − div (gε ) + (λ+ )ε − (λ− )ε be such that fε , gε , (λ+ )ε and (λ− )ε
satisfy hypotheses (13.1.20)–(13.1.23), and let uε be the solution of (13.1.24) in the
sense (13.1.6).
Then the sequence uε has the following properties : there exists a positive constant
c, independent on k, n and ε, such that, for every k > 0, for every n ≥ 0, and for
every ε > 0, we have
Z
1
a(x, ∇Tk (uε )) · ∇Tk (uε ) dx ≤ c , (13.1.29)
k Ω
Z
1
a(x, ∇Tk (uε )) · ∇Tk (uε ) dx ≤ c , (13.1.30)
k {n≤uε < n+k}
Z
1
a(x, ∇Tk (uε )) · ∇Tk (uε ) dx ≤ c . (13.1.31)
k {−n−k < uε ≤−n}
From (13.1.29) it follows, by a result in [4] (see also [36]), that
This implies, again by a result in [4] (see also [36]), that there exists a subsequence
of uε , still denoted by uε , and a measurable function u such that uε converges
almost everywhere to u, and Tk (uε ) converges weakly in W01,p (Ω) to Tk (u) for every
k > 0 ; moreover, u is almost everywhere finite, and there exists a positive constant
c, independent on k and n, such that
Z
1
a(x, ∇Tk (u)) · ∇Tk (u) dx ≤ c ∀k > 0 ,
k Ω
Z
1
a(x, ∇Tk (u)) · ∇Tk (u) dx ≤ c ∀k > 0, ∀n ≥ 0 , (13.1.33)
k {n≤u < n+k}
Z
1
a(x, ∇Tk (u)) · ∇Tk (u) dx ≤ c ∀k > 0, ∀n ≥ 0 , (13.1.34)
k {−n−k < u≤−n}
|∇u|p−1 belongs to Lq (Ω), for every q < NN−1 .
Observe that, in this latter case, ∇u is the “approximate gradient” of u, and not (in
general) its distributional gradient (in (13.1.33) and (13.1.34) it is the distributional
gradient of Tn+k (u)). Furthermore, by a result of [12],
Thus, by (13.1.3), and by the boundedness of |∇uε |p−1 in Lq (Ω), for every q < NN−1 ,
it follows that
a(x, ∇uε ) → a(x, ∇u) strongly in (Lq (Ω))N , for every q < NN−1 . (13.1.35)
In this Section, we give some results about bounded measures on Ω, and define
a fairly general (and suitable for our purposes) way to approximate a measure in
Mb (Ω).
Before giving the results, we recall the definition of p-capacity.
where χK is the characteristic function of K ; we will use the convention that inf Ø =
+∞. The p-capacity of any open subset A of Ω is then defined by :
capp (A, Ω) = sup capp (K, Ω), K compact, K ⊂ A ,
We remark that, once the p-capacity has been defined as before, then, for every
Borelian subset of Ω, Z
p
capp (B, Ω) = inf |∇u| dx ,
Ω
where the infimum is taken over all functions u in W01,p (Ω) such that u = 1 capp -quasi
everywhere on B, and u ≥ 0 capp -quasi everywhere on Ω. In the preceding assertion
we have taken for u its capp -quasi continuous representative.
Definition 13.2.2 We define Mb (Ω) as the set of measures on Ω with bounded total
variation.
We define M0 (Ω) as the set of measures in Mb (Ω) which are absolutely con-
tinuous with respect to the p-capacity, that is, µ belongs to M0 (Ω) if and only if
µ(B) = 0 for every Borelian set B such that capp (B, Ω) = 0.
Proof. Since v belongs to W01,p (Ω), then v has a capp -quasi continuous representa-
tive, which we continue to denote by v. This representative is defined up to a set of
zero p-capacity, so that it is defined µ0 -almost everywhere. Furthermore, since there
exists a non negative constant k such that |v| ≤ k almost everywhere on Ω, we have
that |v| ≤ k capp -quasi everywhere on Ω, and thus µ0 -almost everywhere on Ω. Thus,
v belongs to L∞ (Ω, µ0 ).
We recall that if µ is a measure in Mb (Ω), and E is a borelian subset of Ω, then
the measure µ E is defined by µ E(B) = µ(E ∩ B), for every Borelian subset B
of Ω.
We begin with the following decomposition result.
Proposition 13.2.4 Let µ be a measure in Mb (Ω). Then there exists a unique pair
(µ0 , λ) of measures in Mb (Ω) such that
µ = µ0 + λ ,
E+ ∩ E− = Ø , E+ ∪ E− = E ,
and
λ+ = λ E + , λ− = λ E − .
13.2 Approximation des mesures 127
Thus, we have
capp (E + , Ω) = capp (E − , Ω) = 0 .
Let now δ be a fixed positive number. Due to the regularity of the measures λ+
and λ− , there exist two compact subsets of Ω, Kδ+ and Kδ− , such that
Kδ+ ⊂ E + , Kδ− ⊂ E − ,
Let us explicitly remark that it may happen that Uδ+ ∩E − 6= Ø, and/or Uδ− ∩E + 6= Ø.
Finally, since capp (E + , Ω) = 0, and capp (E − , Ω) = 0, we have
Thus, there exist two functions ψδ+ and ψδ− , and two open sets A+ −
δ and Aδ , with the
following properties :
Kδ+ ⊂ A+δ , Kδ− ⊂ A−δ , (13.2.37)
ψδ+ ∈ Cc∞ (Uδ+ ) , ψδ− ∈ Cc∞ (Uδ− ) , (13.2.38)
0 ≤ ψδ+ ≤ 1 , 0 ≤ ψδ− ≤ 1 , (13.2.39)
ψδ+ ≡ 1 on A+
δ , ψδ− ≡ 1 on A− δ , (13.2.40)
Z Z
|∇ψδ+ |p dx ≤ δ , |∇ψδ− |p dx ≤ δ . (13.2.41)
Uδ+ Uδ−
)From now on, we will consider ψδ+ and ψδ− as functions in Cc∞ (Ω), setting ψδ+ ≡ 0
and ψδ− ≡ 0 on, respectively, Ω\Uδ+ and Ω\Uδ− .
In order to deal with µ0 , we recall the following decomposition result.
µ = µ0 + λ+ − λ− = f − div (g) + λ+ − λ− .
Let now δ and η be two fixed positive real numbers ; since 1 − ψδ+ ψη+ belongs to
C ∞ (Ω), and is identically zero on Kδ+ ∩ Kη+ , with 0 ≤ 1 − ψδ+ ψη+ ≤ 1, we have, using
13.3 Près de E 129
so that Z
(1 − ψδ+ ψη+ ) (λ+ )ε dx = ω(η, δ, ε) . (13.2.45)
Ω
Reasoning in the same way, we get
Z
(1 − ψδ− ψη− ) (λ− )ε dx = ω(η, δ, ε) . (13.2.46)
Ω
13.3 Near E
In this Section we study the behaviour of the approximate solutions uε near the
set E where the measure λ is concentrated ; here the meaning of “near E” will be
specified through the functions ψδ+ and ψδ− (see the previous Section for the definition
of λ, E, ψδ+ and ψδ− ). We prove that, in some sense, the sequence uε of solutions of
(13.1.24) tends to +∞ on a neighbourhood of E + , and to −∞ on a neighbourhood
of E − ; this will reflect on the behaviour of the gradients of Tk (uε ).
We recall here, for the convenience of the reader, that uε is the unique solution
in W01,p (Ω) of
−div (a(x, ∇uε )) = fε − div (gε ) + (λ+ )ε − (λ− )ε in Ω,
(13.3.47)
u =0 ε on ∂Ω,
in the sense (13.1.6). and that fε , gε , (λ+ )ε and (λ− )ε satisfy (13.1.20)–(13.1.23). We
will also suppose that we have already extracted a subsequence, still denoted by uε ,
such that (13.1.29)–(13.1.35) hold true.
Our first result is the following.
Lemma 13.3.1 Let fε , gε , (λ+ )ε and (λ− )ε be sequences of functions satisfying
(13.1.20)–(13.1.23). Let uε be the solution of (13.3.47), and suppose to have extracted
form uε a subsequence, still denoted by uε , such that (13.1.29)–(13.1.35) hold. Let
η be a positive real number, and let ϕ+ and ϕ− be two non negative functions in
W 1,∞ (Ω) such that
Z Z
0≤ +
ϕ− dλ ≤ η , 0≤ ϕ+ dλ− ≤ η . (13.3.48)
Ω Ω
130 Chap. 13 Une généralisation des notions de solutions . . .
We then have
1Z
n
a(x, ∇uε ) · ∇uε ϕ− dx ≤ ωη (n, ε) + η ,
{n≤uε < 2n}
Z (13.3.49)
ϕ− (λ− )ε dx ≤ ωη (n, ε) + η ,
{uε > 2n}
and 1Z
n
a(x, ∇uε ) · ∇uε ϕ+ dx ≤ ωη (n, ε) + η ,
{−2n < uε ≤−n}
Z (13.3.50)
ϕ+ (λ+ )ε dx ≤ ωη (n, ε) + η .
{uε < −2n}
Remark 13.3.2 We are going to comment on the results of the previous lemma. We
will discuss only (13.3.49), since the same comments can be made for (13.3.50).
The first inequality of (13.3.49) says, in some sense, that the energy of uε on the
set {n ≤ uε < 2n}, once divided by n, vanishes as ε tends to zero, and then n tends
to infinity, on the set where λ+ is not concentrated. In contrast, note that
Z
1
a(x, ∇uε ) · ∇uε dx ,
n {n≤uε < 2n}
is bounded with respect to both n and ε (see (13.1.30) with k = n and use the
coercivity (13.1.2)), but does not converge to zero as ε tends to zero and then n tends
to infinity (see (13.5.108), in Section 13.5).
The second part of (13.3.49) describes a similar fact, for what concerns (λ− )ε :
even if we use a test function which is concentrated where (λ− )ε is concentrated
(indeed, no hypotheses are required on the behaviour of ϕ− with respect to λ− ),
the fact that we restrict our attention to the set where uε is larger than 2n yields a
quantity that converges to zero.
Proof of Lemma 13.3.1. Let βn+ (s) = Bn,n (s+ ), where Bn,n is defined in (13.1.28) ;
if we consider the sequence βn+ (uε ), we then have
Z Z
+ p 1 c
|∇βn (uε )| dx = p |∇uε |p dx ≤ p−1 , (13.3.51)
Ω n {n≤uε < 2n} n
βn+ (uε ) → βn+ (u) almost everywhere and weakly∗ in L∞ (Ω). (13.3.52)
Since βn+ (0) = 0, and since βn+ (s) is bounded by 1, (13.3.54) implies that
and that
βn+ (u) → 0 weakly in W01,p (Ω). (13.3.56)
We now choose βn+ (uε ) ϕ− as test function in (13.3.47) (this choice is admissible
since the function belongs to W01,p (Ω)). We obtain
Z
1
a(x, ∇uε ) · ∇uε ϕ− dx (A)
n {n≤uε < 2n}
Z
+ a(x, ∇uε ) · ∇ϕ− βn+ (uε ) dx (B)
Ω
Z
= fε βn+ (uε ) ϕ− dx (C)
Ω
− hdiv (gε ), βn+ (uε ) ϕ− i (D)
Z
+ βn+ (uε ) ϕ− (λ+ )ε dx (E)
ZΩ
− βn+ (uε ) ϕ− (λ− )ε dx . (F)
Ω
Recalling that a(x, ∇uε ) converges to a(x, ∇u) strongly in (Lq (Ω))N for every q < NN−1
by (13.1.35), that ϕ− belongs to W 1,∞ (Ω), that (13.3.52) holds, and using (13.3.55),
we obtain
Z
(B) = a(x, ∇u) · ∇ϕ− βn+ (u) dx + ωn (ε) = ω(n, ε) . (13.3.57)
Ω
Furthermore,
(D) = −hdiv (g), βn+ (u) ϕ− i + ωn (ε) = ω(n, ε) , (13.3.59)
0
due to the strong convergence (13.1.21) of −div (gε ) to −div (g) in W −1,p (Ω) and to
(13.3.53) (for the first limit), and to (13.3.56) (for the second). Finally, since βn+ (uε )
is non negative, and ϕ− is continuous,
Z Z
+
(E) ≤ ϕ− (λ )ε dx = ϕ− dλ+ + ω(ε) ≤ η + ω(ε) , (13.3.60)
Ω Ω
132 Chap. 13 Une généralisation des notions de solutions . . .
where we have used (13.3.48) in the last passage. Thus, observing that
Z
−(F) ≥ ϕ− (λ− )ε dx , (13.3.61)
{uε > 2n}
Lemma 13.3.3 Let k be a positive real number. Let fε , gε , (λ+ )ε and (λ− )ε
be sequences of functions satisfying (13.1.20)–(13.1.23). Let uε be the solution of
(13.3.47), and suppose to have extracted form uε a subsequence, still denoted by uε ,
such that (13.1.29)–(13.1.35) hold. Let ψδ+ and ψδ− , as well as ψη+ and ψη− , be functions
which satisfy (13.2.38)–(13.2.41). Then the following holds
Z
a(x, ∇Tk (uε )) · ∇Tk (uε ) ψδ+ ψη+ dx = ω(η, δ, ε) ,
Z Ω (13.3.62)
+ + +
(k − Tk (uε )) ψδ ψη (λ )ε dx = ω(η, n, δ, ε) .
{−n≤uε ≤k}
Z
a(x, ∇Tk (uε )) · ∇Tk (uε ) ψδ− ψη− dx = ω(η, δ, ε) ,
Ω
Z (13.3.63)
(k + Tk (uε )) ψδ− ψη− −
(λ )ε dx = ω(η, n, δ, ε) .
{−k≤uε ≤n}
Remark 13.3.4 As in Remark 13.3.2, some comments are in order. The first result
in (13.3.62) can be seen as a result giving some properties of Tk (uε ) and Tk (u) near the
set E + . Indeed, using the almost everywhere convergence of a(x, ∇Tk (uε )) · ∇Tk (uε )
to a(x, ∇Tk (u))·∇Tk (u) and Fatou lemma (since, by the coercivity condition (13.1.2),
a(x, ∇Tk (uε )) · ∇Tk (uε ) is non negative), we obtain
Z
a(x, ∇Tk (u)) · ∇Tk (u) ψδ+ ψη+ dx = ω(η, δ) ,
Ω
for every k > 0. This means that u is very large near the set E + , so that Tk (u) is
equal to k, and thus its gradient is zero near E + .
The second inequality of (13.3.62) states the same fact in terms of λ+ : the set
where u is smaller than k, has little (and, actually, zero) measure with respect to
λ+ . Since k is arbitrary, this means, grosso modo, that u is λ+ almost everywhere
positively infinite on E + .
The same remarks can be made on (13.3.63).
13.3 Près de E 133
Proof of Lemma 13.3.3. Let k > 0 be fixed, and let n in N be such that n > k.
Let hn (s) = Hn,n (s), where Hn,n is defined in (13.1.27). Reasoning as in the proof of
Lemma 13.3.1 (that is, using again (13.1.30)), and observing that hn (uε ) is bounded
by 1, we get that for n fixed
hn (uε ) → hn (u) almost everywhere and weakly∗ in L∞ (Ω), (13.3.64)
hn (uε ) → hn (u) weakly in W 1,p (Ω). (13.3.65)
We choose as test function in (13.3.47)
(k − Tk (uε )) hn (uε ) ψδ+ ψη+ ,
+ a(x, ∇uε ) · ∇uε h0n (uε )(k − Tk (uε )) ψδ+ ψη+ dx (B)
ZΩ
+ a(x, ∇uε ) · ∇ψδ+ hn (uε )(k − Tk (uε )) ψη+ dx (C)
ZΩ
+ a(x, ∇uε ) · ∇ψη+ hn (uε )(k − Tk (uε )) ψδ+ dx (D)
Ω
Z
= fε (k − Tk (uε )) hn (uε ) ψδ+ ψη+ dx (E)
Ω
− hdiv (gε ), (k − Tk (uε )) hn (uε ) ψδ+ ψη+ i (F)
Z
+ (k − Tk (uε )) hn (uε ) ψδ+ ψη+ (λ+ )ε dx (G)
ZΩ
− (k − Tk (uε )) hn (uε ) ψδ+ ψη+ (λ− )ε dx . (H)
Ω
Since n is larger than k, then k − Tk (uε ) = 2k on the set {−2n < uε ≤ −n}, and
k − Tk (uε ) = 0 on the set {n ≤ uε < 2n}, we have
Z
1
(B) = a(x, ∇uε ) · ∇uε (k − Tk (uε )) ψδ+ ψη+ dx
n {−2n
Z < uε ≤−n}
1
− a(x, ∇uε ) · ∇uε (k − Tk (uε )) ψδ+ ψη+ dx .
n
Z {n≤uε < 2n}
2k
= a(x, ∇uε ) · ∇uε ψδ+ ψη+ dx ,
n {−2n < uε ≤−n}
and so, since the integrand functions are non negative, and ψδ+ ≤ 1, we have
Z
2k
(B) ≤ a(x, ∇uε ) · ∇uε ψη+ dx .
n {−2n < uε ≤−n}
Thus, by (13.3.50), which we can apply since ϕ+ = ψη+ is such that (13.3.48) holds
thanks to (13.2.44),
(B) = ω(η, n, ε) . (13.3.66)
134 Chap. 13 Une généralisation des notions de solutions . . .
Furthermore, since for k fixed k −Tk (uε ) converges to k −Tk (u) in the weak∗ topology
of L∞ (Ω), and since (13.3.64) holds, we have, since supp (hn ) = [−2n, 2n],
Z
(C) = a(x, ∇T2n (u)) · ∇ψδ+ hn (u)(k − Tk (u)) ψη+ dx + ωη,n,δ (ε) = ωη,n (δ, ε) ,
Ω
(13.3.67)
where the last statement is due to the fact that ψδ+ converges strongly to zero
0
in W01,p (Ω) (see (13.2.41)), that a(x, ∇T2n (u)) hn (u) belongs to (Lp (Ω))N , while
(k − Tk (u)) ψη+ belongs to L∞ (Ω). Similarly, we have
Z
(D) = a(x, ∇T2n (u)) · ∇ψη+ hn (u)(k − Tk (u)) ψδ+ dx + ωη,n,δ (ε) = ωη,n (δ, ε) .
Ω
(13.3.68)
since fε converges weakly to f in L1 (Ω), since (13.3.64) holds, and since k − Tk (uε )
converges to k − Tk (u) weakly∗ in L∞ (Ω) and almost everywhere ; the second limit
is performed using again the fact that ψδ+ converges to zero in the weak∗ topology of
L∞ (Ω), while the term f (k − Tk (u)) hn (u) ψη+ belongs to L1 (Ω).
We then have
and Z
(G) ≥ (k − Tk (uε )) ψδ+ ψη+ (λ+ )ε dx ,
{−n≤uε ≤k}
and using the corresponding properties of ψδ− , ψη− , (λ+ )ε and (λ− )ε .
Lemma 13.4.1 Let k be a positive real number. Let fε , gε , (λ+ )ε and (λ− )ε
be sequences of functions satisfying (13.1.20)–(13.1.23). Let uε be the solution of
(13.3.47), and suppose to have extracted form uε a subsequence, still denoted by uε ,
such that (13.1.29)–(13.1.35) hold. Let ψδ+ and ψδ− , as well as ψη+ and ψη− , be functions
which satisfy (13.2.38)–(13.2.41), and define, for δ > 0 and η > 0
Remark 13.4.2 The meaning of (13.4.74) is, roughly speaking, that a(x, ∇Tk (uε )) ·
∇Tk (uε ) strongly converges to a(x, ∇Tk (u)) · ∇Tk (u) in L1 (Ω) if we stay away from
E.
We split the proof into various lemmas. We begin with the following result.
Proof. We choose (1 − Φδ,η ) Tk (uε ) as test function in (13.3.47) (this can be done
since this function belongs to W01,p (Ω)) ; we get
Z
a(x, ∇Tk (uε )) · ∇Tk (uε ) (1 − Φδ,η ) dx (A)
Ω Z
− a(x, ∇uε ) · ∇Φδ,η Tk (uε ) dx (B)
Ω
Z
= fε (1 − Φδ,η ) Tk (uε ) dx (C)
Ω
− hdiv (gε ), (1 − Φδ,η ) Tk (uε )i (D)
Z
+ (1 − Φδ,η ) Tk (uε ) (λ+ )ε dx (E)
ZΩ
− (1 − Φδ,η ) Tk (uε ) (λ− )ε dx . (F)
Ω
Since Φδ,η belongs to Cc∞ (Ω), since a(x, ∇uε ) converges to a(x, ∇u) strongly in
(Lq (Ω))N for every q < NN−1 by (13.1.35), and since Tk (uε ) converges weakly∗ in L∞ (Ω)
and almost everywhere to Tk (u), we have
Z
(B) = a(x, ∇u) · ∇Φδ,η Tk (u) dx + ωη,δ (ε) . (13.4.76)
Ω
Since for η and δ fixed, (1 − Φδ,η ) Tk (uε ) converges to (1 − Φδ,η ) Tk (u) weakly∗ in
L∞ (Ω) and almost everywhere in Ω, while fε converges weakly to f in L1 (Ω), we get,
by Lemma 13.1.21,
Z
(C) = f (1 − Φδ,η ) Tk (u) dx + ωη,δ (ε) . (13.4.77)
Ω
Since for η and δ fixed, (1 − Φδ,η ) Tk (uε ) converges to (1 − Φδ,η ) Tk (u) weakly in
0
W01,p (Ω), and since −div (gε ) converges to −div (g) strongly in W −1,p (Ω), we obtain
(D) = −hdiv (g), (1 − Φδ,η ) Tk (u)i + ωη,δ (ε) . (13.4.78)
For (E), we have, using the tight convergence of (λ+ )ε to λ+ , the definition
(13.4.73) of Φδ,η , and (13.2.45) and (13.2.43),
Z Z
+
|(E)| ≤ k (1 − Φδ,η ) (λ )ε dx = k (1 − Φδ,η ) dλ+ + ωη,δ (ε)
ZΩ Z Ω
≤ k (1 − ψδ+ ψη+ ) dλ+ + k ψδ− ψη− dλ+ + ωη,δ (ε)
ZΩ Z Ω (13.4.79)
+ + − +
≤ k (1 − ψη ) dλ + k ψη dλ + ωη (δ, ε)
Ω Ω
= ω(η, δ, ε) .
For every term above we can apply the result of Lemma 13.3.1. Indeed, if we define
ϕ− = 1 − ψδ+ ψη+ , we have, by (13.2.45),
Z
ϕ− dλ+ ≤ δ + η ,
Ω
Proof. Let k > 0 be fixed, let n ∈ N be such that n > k, and define hn (s) = Hn,n (s)
as in the proof of Lemma 13.3.3 ; we thus have that hn (uε ) satisfies (13.3.64) and
(13.3.65). Moreover, using the definition of hn and (13.3.64) and (13.3.65) with k = n,
we have
Since a(x, ∇uε ) hn (uε ) = a(x, ∇T2n (uε )) hn (uε ), the sequence a(x, ∇uε ) hn (uε )
0
converges to a(x, ∇T2n (u)) hn (u) weakly in (Lp (Ω))N as ε tends to zero. Thus, we
have, since n > k,
Z
(A) = a(x, ∇T2n (u)) · ∇Tk (u) (1 − Φδ,η ) hn (u) dx + ωη,δ,n (ε)
ZΩ (13.4.85)
= a(x, ∇Tk (u)) · ∇Tk (u) (1 − Φδ,η ) dx + ωη,δ (n, ε) .
Ω
13.5 Démonstration des résultats 139
Moreover, we have
Z
(B) = − a(x, ∇u) · ∇Φδ,η Tk (u) hn (u) dx + ωη,δ,n (ε)
ZΩ (13.4.86)
= − a(x, ∇u) · ∇Φδ,η Tk (u) dx + ωη,δ (n, ε) ,
Ω
where the first statement holds since a(x, ∇uε ) converges strongly to a(x, ∇u) in
(Lq (Ω))N , for every q < NN−1 , Φδ,η belongs to Cc∞ (Ω), and hn (uε ) satisfies (13.3.64) ;
the second one is due to (13.4.83). Furthermore, by the result of Lemma 13.4.4,
Z
k
|(C)| ≤ a(x, ∇uε ) · ∇uε (1 − Φδ,η ) dx = ω(η, δ, n, ε) , (13.4.87)
n {n≤|uε | < 2n}
As for the right hand side, we have
Z
(D) = f Tk (u) (1 − Φδ,η ) hn (u) dx + ωη,δ,n (ε)
ZΩ (13.4.88)
= f Tk (u) (1 − Φδ,η ) dx + ωη,δ (n, ε) ,
Ω
by (13.3.64), and by the weak convergence of fε to f weakly in L1 (Ω) (we have applied
again Lemma 13.1.21) ; the second statement is due to (13.4.83). Moreover,
where the first statement holds since −div (gε ) converges to −div (g) strongly in
0
W −1,p (Ω), while hn (uε ) satisfies (13.3.65), while the second one is due to the fact
that hn (u) satisfies (13.4.84). The two remaining terms are estimated as follows using
(13.2.45) and (13.2.43) :
Z
|(F)| ≤ k (1 − Φδ,η ) (λ+ )ε dx ≤ ω(η, δ, ε) . (13.4.90)
Ω
Analogously, we get
|(G)| ≤ ω(η, δ, ε) . (13.4.91)
Putting together (13.4.85)–(13.4.91), we have proved (13.4.82), since the right hand
isde of (13.4.82) does not depend on ε.
Proof of Lemma 13.4.1. It is enough to put together the estimates (13.4.75) and
(13.4.82) in order to obtain (13.4.74).
Proof of Theorem 13.1.20. Let η and δ be fixed positive real numbers, and
let Φδ,η = 1 − ψδ+ ψη+ − ψδ− ψη− . Let uε be the sequence of solutions of (13.3.47),
and suppose to have extracted from it a subsequence, still denoted by uε , such that
(13.1.29)–(13.1.35) hold. We have
Z
[a(x, ∇Tk (uε )) · ∇Tk (uε ) − a(x, ∇Tk (u)) · ∇Tk (u)] dx (A)
Z Ω
= [a(x, ∇Tk (uε )) · ∇Tk (uε ) − a(x, ∇Tk (u)) · ∇Tk (u)] (1 − Φδ,η ) dx (B)
Ω
Z
+ a(x, ∇Tk (uε )) · ∇Tk (uε ) (ψδ+ ψη+ + ψδ− ψη− ) dx (C)
ZΩ
− a(x, ∇Tk (u)) · ∇Tk (u) (ψδ+ ψη+ + ψδ− ψη− ) dx . (D)
Ω
Since a(x, ∇Tk (uε )) · ∇Tk (uε ) is a sequence of non negative functions that converges
almost everywhere to a(x, ∇Tk (u)) · ∇Tk (u), (13.5.95) implies that
a(x, ∇Tk (uε )) · ∇Tk (uε ) → a(x, ∇Tk (u)) · ∇Tk (u) strongly in L1 (Ω).
But then, by (13.1.2), we have
α |∇Tk (uε )|p ≤ a(x, ∇Tk (uε )) · ∇Tk (uε ) ,
so that, since ∇Tk (uε ) converges almost everywhere, a generalized version of the
Lebesgue theorem yields that ∇Tk (uε ) is strongly compact in (Lp (Ω))N ; since its
almost everywhere limit is ∇Tk (u), we obtain that
∇Tk (uε ) → ∇Tk (uε ) strongly in (Lp (Ω))N ,
and this concludes the proof.
13.5 Démonstration des résultats 141
Remark 13.5.1 We want to remark that if the sequences (λ+ )ε and (λ− )ε are not
made of non negative functions, the result may not be true. This is the case of the
following two examples.
Example 13.5.2 In [20] (see also [17], where this example has been quoted in
relation with strong convergence of truncations) it is considered a sequence uε of
functions in H01 (Ω) that solve the following problem
−∆uε = µε − λε in Ω,
u =0 on ∂Ω,
ε
where µε and λε are two sequences of non negative elements in H −1 (Ω), which both
converge to the Lebesgue measure on Ω, the first in the weak∗ topology of measures,
and the second weakly in H −1 (Ω) (we are thus outside our framework, since we require
the term in H −1 (Ω) to be strongly convergent). Since 0 ≤ uε ≤ 1, and since uε does
not converge strongly in H01 (Ω), we thus have that T1 (uε ) does not converge strongly
in H01 (Ω), and so the result of Theorem 13.1.20 does not hold.
One may think that the result fails since we are approximating a measure which
is not concentrated on a set of zero capacity, but this is not true, as the following
example shows.
Example 13.5.3 Let N = 4, and let Ω = B1 (0) = {x ∈ R4 : |x| < 1} ; consider the
sequence uε of solutions of
−∆uε = (λ+ )ε − (λ− )ε in Ω,
u =0 on ∂Ω,
ε
where
1 1
(λ+ )ε = (χBε (0) + χB √4 6 ε (0)\B √4 5 ε (0) ) , (λ− )ε = χB √4 (0)\B √4 (0) .
ε4 ε4 5ε 3ε
It is easily seen that both (λ+ )ε and (λ− )ε converge in the weak∗ topology of measures
to the Dirac’s delta (times a constant) concentrated in the origin. Since the solution
uε is radial, it can be explicitly calculated ; it is then easy to see that uε ≡ 0 on
B1 (0)\B √ 4
6 ε (0), so that uε converges to zero almost everywhere in Ω. Furthermore,
since (λ )ε +(λ− )ε is bounded in L1 (Ω), then uε is bounded in W01,q (Ω), for every q < 43
+
(see [12]). This implies (by Rellich theorem) that uε is strongly convergent in L1 (Ω)
(for example). Since its almost everywhere limit is zero, then uε converges strongly
to zero in L1 (Ω). Let us study the behaviour of T√k (uε ). Using
√ again the explicit form
4 4
of uε , it is easy to see that, for ρ = |x| between 2 ε and 3 ε, we have
1 1 c
uε (ρ) = − −1 ,
8 ρ2 ε 2
142 Chap. 13 Une généralisation des notions de solutions . . .
√
where c is a positive constant, independent on ε, such√that uε ( 4 2 ε) is positive (and
actually diverges to +∞ as ε tends to zero), while uε ( 4 3 ε) is negative (and actually
diverges to −∞ as ε tends to zero). Thus,√if k is√fixed, and ε is small enough, the set
−
{|uε | ≤ k} contains a subset (ρ+ 4
ε , ρε ) of ( 2 ε,
4
3 ε), where
1 c 1 c
+ 2
= 2 + 8k + 1 , + 2
= 2 − 8k + 1 .
(ρε ) ε (ρε ) ε
so that Tk (uε ) does not converge to Tk (u) = 0 strongly in H01 (Ω). Hence, the result
of Theorem 13.1.20 does not hold. We remark that we have (again by explicit
calculations) that
Z Z
1 2 1
|∇uε | dx = 1 + ω(n, ε) = |∇uε |2 dx .
n {n≤uε < 2n} n {−2n < uε ≤−n}
We give now the proof of Theorem 13.1.16. The proof will be splitted in four
steps. We begin by proving that the function u which is the limit of the subsequence
of uε given by Theorem 13.1.20 satisfies (13.1.15), and then that it is a renormalized
solution in the sense of Definition 13.1.10. The third step is the proof that the limit
function u satisfies (13.1.13) and (13.1.14), so that, together with Step 1, we have
that u is a renormalized solution of (13.0.1) in the sense of Definition 13.1.7 ; we
finally prove that u is a renormalized solution in the sense of Definition 13.1.2.
Step 1 : obtaining (13.1.15). Let h be a function in W 1,∞ (R) such that h has
compact support in R, and let ϕ be a function in W01,p (Ω) ∩ L∞ (Ω). We consider the
sequence of functions h(uε ). If supp (h) ⊆ [−M, M ], we have
Since TM (uε ) is strongly convergent in (Lp (Ω))N by the result of Theorem 13.1.20,
and since ∇h(uε ) is almost everywhere convergent to ∇h(u) (due to the continuity of
h and the almost everywhere convergence of uε and ∇uε to u and ∇u respectively),
we have that ∇h(uε ) is strongly convergent to ∇h(u) in (Lp (Ω))N by a generalized
version of the Lebesgue theorem. Moreover, since h is bounded, we also have that
h(uε ) is bounded in Lr (Ω) for every r ; since it is almost everywhere convergent, we
have that it is strongly convergent to h(u) in Lp (Ω). Thus, h(uε ) converges to h(u)
strongly in W 1,p (Ω). Furthermore, the boundedness of h, and the almost everywhere
convergence of h(uε ), imply that h(uε ) converges to h(u) in the weak∗ topology of
L∞ (Ω).
13.5 Démonstration des résultats 143
since we have the strong convergence of ∇h(uε ) and ∇TM (uε ) in (Lp (Ω))N , and since
ϕ belongs to L∞ (Ω). Furthermore
Z
(B) = a(x, ∇u) · ∇ϕ h(u) dx + ω(ε) , (13.5.97)
Ω
since a(x, ∇uε ) h(uε ) = a(x, ∇TM (uε )) h(uε ), and so it is strongly convergent in
0
(Lp (Ω))N to a(x, ∇u) h(u) as ε tends to zero ; the result then follows since ϕ belongs
to W01,p (Ω). As for the right hand side, we have
Z
(C) = f h(u) ϕ dx + ω(ε) , (13.5.98)
Ω
by (13.2.45). For the second term, we begin remarking that if h belongs to W 1,∞ (R),
and has compact support in R, then there exists a positive real number M1 such that
|h(t)| ≤ M1 − TM1 (t) ∀t ∈ R .
Thus
Z
kϕkL∞ (Ω) (λ+ )ε |h(uε )| ψδ+ ψη+ dx
{|uε |≤MZ}
(13.5.101)
≤ kϕkL∞ (Ω) (λ+ )ε (M1 − TM1 (uε )) ψδ+ ψη+ dx = ω(η, δ, ε) ,
{|uε |≤M }
for every h in W (R) with compact support in R, and for every ϕ in W01,p (Ω) ∩
1,∞
and every term, but (B), (E) and (F), can be treated exactly as in Step 1. We have
Z
(B) = a(x, ∇u) · ∇ϕ h(u) dx + ω(ε) ,
Ω
Since a(x, ∇uε ) converges strongly to a(x, ∇u) in (Lq (Ω))N , for every q < NN−1 by
(13.1.35), since ϕ is in W01,r (Ω) for r > N , and since h(uε ) is weakly∗ convergent in
L∞ (Ω) to h(u). As for the other terms, we have
Z Z
+ +∞
(E) = (λ )ε h ϕ dx + (λ+ )ε (h(uε ) − h+∞ ) ϕ dx .
Ω Ω
Z
+∞
h ϕ dλ+ . (13.5.104)
Ω
As for the second term, we observe that h(uε )−h+∞ is zero on the set uε > M2 , where
M2 is such that supp (h0 ) ⊆ [−M2 , M2 ], and so
Z
(λ+ )ε (h(uε ) − h+∞ ) ϕ dx
Ω Z
≤ kϕkL∞ (Ω) (λ+ )ε |h(uε ) − h+∞ | dx
{uε ≤M2 } Z
(13.5.105)
≤ 2kϕkL∞ (Ω) khkL∞ (R) (λ+ )ε (1 − ψδ+ ψη+ ) dx
Z {uε ≤M2 }
As in Step 1, the first term of the right hand side is such that
Z
2kϕkL∞ (Ω) khkL∞ (R) (λ+ )ε (1 − ψδ+ ψη+ ) dx = ω(η, δ, ε) , (13.5.106)
{uε ≤M2 }
thanks to (13.2.45). Moreover, since there exists a positive real number M3 such that
the second term in the right hand side of (13.5.105) can be estimated as follows (we
146 Chap. 13 Une généralisation des notions de solutions . . .
= ω(η, n, ε) .
To obtain the result, we have used (13.3.50) (to deal with the integral on the set
{uε < −n}), and (13.3.62) (to deal with the integral on the set {−n ≤ uε ≤ M2 }).
Using (13.5.106) and (13.5.107) in (13.5.105), we get
Z
(λ+ )ε (h(uε ) − h+∞ ) ϕ dx = ω(ε) .
Ω
for every h in W 1,∞ (R) such that h0 has compact support in R, and for every ϕ
in W01,r (Ω) for r > N , which is exactly (13.1.17) if we identify, as usual, µ0 with
f − div (g).
Step 3 : obtaining (13.1.13) and (13.1.14). We will prove a more general
version of both (13.1.13) and (13.1.14). Let αn and βn be sequences of positive real
numbers such that
lim βn = +∞ , βn > αn ∀n ∈ N .
n→+∞
13.5 Démonstration des résultats 147
To prove this fact, we choose Bαn ,βn −αn (u+ ) ϕ as test function in (13.1.17), where
Bαn ,βn −αn is as in (13.1.28). Such a function is an admissible test function in (13.1.17)
since the support of the derivative of Bαn ,βn −αn (s+ ) is [αn , βn ], and Bαn ,βn +αn (0) = 0.
We obtain
Z
a(x, ∇u) · ∇ϕ Bαn ,βn −αn (u+ ) dx (A)
Ω Z
1
+ a(x, ∇u) · ∇u ϕ dx (B)
βn − αn {αn ≤u < βn }
Z
= f ϕ Bαn ,βn −αn (u+ ) dx (C)
Ω
− hdiv (g), ϕ Bαn ,βn −αn (u+ )i (D)
Z
+ ϕ dλ+ . (E)
Ω
It is clear that (13.5.108) will follow from the above identity if we prove that
This is true, since Bαn ,βn −αn (u+ ) converges to zero in the weak∗ topology of L∞ (Ω)
(and this takes into account |(A)| and |(C)|), and in the weak topology of W01,p (Ω),
since it is bounded in W01,p (Ω) and converges almost everywhere to zero as n tends
to infinity (and this gives the result for |(D)|).
Repeating the same steps as before, we have
Z Z
1
a(x, ∇u) · ∇u ϕ dx = ϕ dλ− + ω(n) . (13.5.109)
βn − αn {−βn < u≤−αn } Ω
Analogously,
1
a(x, ∇u) · ∇u χ{−βn < u≤−αn } → λ−
βn − αn (13.5.111)
in the weak∗ topology of measure.
148 Chap. 13 Une généralisation des notions de solutions . . .
Step 4 : obtaining (13.1.10). Let w be a function in W01,p (Ω) ∩ L∞ (Ω) such that
there exists k > 0 and w+∞ and w−∞ in W 1, r(Ω) for r > N such that
(
w = w+∞ capp -quasi everywhere on the set {u > k},
(13.5.112)
w = w−∞ capp -quasi everywhere on the set {u < −k}.
We choose as test function in (13.1.15) Hn,1 (u) w, where Hn,1 (s) has been defined in
(13.1.27). We remark that this choice is possible since Hn,1 (s) belongs to W 1,∞ (R)
and has compact support in R. We obtain :
Z
a(x, ∇u) · ∇w Hn,1 (u) dx (A)
Ω
Z
+ a(x, ∇u) · ∇u w dx (B)
{−n−1 < u≤−n}
Z
− a(x, ∇u) · ∇u w dx (C)
{n≤u < n+1}
Z
= f w Hn,1 (u) dx (D)
Ω
− hdiv (g), w Hn,1 (u)i . (E)
Since Hn,1 (u) converges to 1 weakly in W 1,p (Ω) and weakly∗ in L∞ (Ω) as n tends to
infinity, we have Z
(A) = a(x, ∇u) · ∇w dx + ω(n) , (13.5.113)
Ω
and Z
(D) + (E) = f w dx − hdiv (g), wi + ω(n) . (13.5.114)
Ω
On the other hand we have, if n > k (and this is not a restriction since n will tend to
infinity), Z
(B) = a(x, ∇u) · ∇u w−∞ dx ,
{−n−1 < u≤−n}
and Z
(C) = − a(x, ∇u) · ∇u w+∞ dx .
{n≤u < n+1}
that is, (13.1.10) holds true once we identify µ0 with f − div (g). This concludes Step
4 and the proof of the theorem.
13.6 Équivalence entre les définitions 149
(i) implies (ii). If u satisfies (i), and h and ϕ are as in Definition 13.1.7, then it is
possible to choose w = h(u) ϕ as test function in (13.1.10), setting w+∞ = w−∞ = 0
and k = M , where M is such that supp (h) ⊆ [−M, M ]. We thus obtain (13.1.15).
In order to prove that (13.1.13) holds true, we choose w = Bn,n (u+ ) ϕ, where Bn,n is
as in (13.1.28), and ϕ belongs to W 1,r (Ω) for r > N . We obtain, since w+∞ = 1 and
w−∞ = 0,
Z Z
1
a(x, ∇u) · ∇u ϕ dx + a(x, ∇u) · ∇ϕ Bn,n (u+ ) dx
n {n≤u < 2n}Z Z Ω
= Bn,n (u+ ) ϕ dµ0 + ϕ dλ+ .
Ω Ω
and that Z
Bn,n (u+ ) ϕ dµ0 = ω(n) ,
Ω
so that Z Z
1
a(x, ∇u) · ∇u ϕ dx = ϕ dλ+ + ω(n) ,
n {n≤u < 2n} Ω
(ii) implies (iii). If h is as in Definition 13.1.10, let M be such that supp (h0 ) ⊆
[−M, M ]. Let n be greater than M , and define the following function
0 if s < −2n,
s + 2n
h−∞
if −2n ≤ s < −n,
n
hn (s) = h(s) if |s| ≤ n,
2n − s
h+∞
if n < s ≤ 2n,
n
0 if s > 2n.
Since Hk−δ,2δ (u) converges to χ{|u|≤k} in the weak∗ topology of L∞ (Ω) as δ tends to
0
zero, and since a(x, ∇u) Hk−δ,2δ (u) = a(x, ∇Tk+δ (u)) Hk−δ,2δ (u) belongs to (Lp (Ω))N ,
we thus have Z
(C) = a(x, ∇Tk (u)) · ∇ϕ dx + ω(δ) . (13.6.116)
Ω
Reasoning in the same way, we have
Z
(D) = ϕ χ{|u|≤k} dµ0 + ω(δ) . (13.6.117)
Ω
Putting together the latter two identities, (13.6.116) and (13.6.117), we have, for
every ϕ in W01,r (Ω) for r > N ,
Z Z Z Z
+ −
ϕ dλk − ϕ dλk = a(x, ∇Tk (u)) · ∇ϕ dx − ϕ χ{|u|≤k} dµ0 , (13.6.118)
Ω Ω Ω Ω
which coincides with (13.1.19) if we restrict the class of test functions to W01,r (Ω)
−
for r > N . In order to prove that λ+
k and λk satisfy all the required properties, we
152 Chap. 13 Une généralisation des notions de solutions . . .
−
begin observing that since ρ+δ,k and ρδ,k are sequences of non negative functions, then
−
both λ+k and λk are non negative measures. Furthermore, due to the fact that (for
example) the support of ρ+δ,k is contained in the set {k − η ≤ u < k + η} if η > δ, then
+
the support of λk is contained in {k − η ≤ u < k + η} for every η > 0, and so (13.1.18)
holds true. Furthermore, (13.6.118) implies that
−
λ+
k − λk = −div (a(x, ∇Tk (u)) + µ0 {|u| ≤ k} ,
−
so that (by the regularity of Tk (u) and µ0 ), we have that λ+ k − λk belongs to
0
W −1,p (Ω) + L1 (Ω), that is, it belongs to M0 (Ω) (see Proposition 13.2.5). This
implies that we can extend the set of admissible test functions from Cc1 (Ω) to
W01,p (Ω) ∩ L∞ (Ω). Let now K be a compact subset of Ω with zero p-capacity. Then
there exists a sequence ϕδ of Cc∞ (Ω) functions such that 0 ≤ ϕδ ≤ 1 on Ω, ϕδ ≡ 1 on
K, and ϕδ converges to zero strongly in W01,p (Ω), almost everywhere and ∗-weakly in
L∞ (Ω). We fix h < k and choose Th (u+ ) ϕδ as test function in (13.6.118). We obtain,
−
recalling the support properties of λ+ k and λk ,
Z Z
a(x, ∇Tk (u)) · ∇Th (u ) ϕδ dx + a(x, ∇Tk (u)) · ∇ϕδ Th (u+ ) dx
+
Ω Z ΩZ
= + χ
ϕδ Th (u ) {|u|≤k} dµ0 + h ϕδ dλ+
k .
Ω Ω
Since λ+
k and ϕδ are positive, we have, recalling that ϕδ ≡ 1 on K,
Z
+
0 ≤ λk (K) ≤ ϕδ dλ+
k .
Ω
so that λ+k (K) = 0 for every subset K of zero p-capacity. This implies that λk is
+
zero on every borelian set of zero p-capacity ; hence, in view of Proposition 13.2.5, it
belongs to M0 (Ω). The same calculations (taking Th (u− ) ϕδ as test function) yield
the same result for λ− k.
Observe that λ+k and λ− +
k do not depend on the subsequence extracted from ρδ,k and
ρ−
δ,k respectively, so that the whole sequences are convergent in the weak∗ topology
of measures.
To prove the weak∗ convergence of λ+ + +
k to λ , we fix k > 0, and choose Bk,k (u ) ϕ
in (13.1.19) written for 2k + 1, with ϕ in W01,r (Ω) for r > N . We obtain
Z
Bk,k (u+ )ϕ dλ+
2k+1 (A)
Ω Z
Since supp (λ+ 2k+1 ) ⊆ {2k ≤ u ≤ 2k + 2} (for example), and since on this set
Bk,k (u) ≡ 1, we have Z
(A) = ϕ dλ+
2k+1 . (13.6.119)
Ω
On the other hand, since supp (λ−
2k+1 ) ⊆ {−2k − 2 ≤ u ≤ −2k}, and since on this set
+
Bk,k (u ) ≡ 0, we have
(B) = 0 . (13.6.120)
Since Bk,k (u+ ) converges to 0 as k tends to infinity, both almost everywhere and in
the weak∗ topology of L∞ (Ω), and since a(x, ∇Tk (u)) converges to a(x, ∇u) strongly
in (Lq (Ω))N for every q < NN−1 , we have, being ϕ in W01,r (Ω) for r > N ,
(D) = ω(k) . (13.6.121)
The same ideas imply that
(E) = ω(k) . (13.6.122)
Finally, we have, by (13.1.13),
Z Z
1
(C) = a(x, ∇u) · ∇u ϕ dx = ϕ dλ+ + ω(k) .
k {k≤u < 2k} Ω
Passing to the limit on k, which is possible on every term due to the regularity of u
and w, we have
Z Z
a(x, ∇Tk (u)) · ∇w dx = a(x, ∇u) · ∇w dx + ω(k) .
Ω Ω
The second term is fixed, and thus gives no problem ; as for the third (the fourth
being identical), if k > η we have, recalling the support property of λ+
k,
Z Z Z
+ +∞ +
w dλk = w dλk = w+∞ dλ+ + ω(k) ,
Ω Ω Ω
since w+∞ is continuous. Putting together the four terms, we obtain (13.1.10).
[5] Ph. Bénilan, H. Brezis, M.G. Crandall, A semilinear equation in L1 , Ann. Scuola
Norm. Sup. Pisa. Cl. Sci. 2 (1975), 523–555.
[12] L. Boccardo, T. Gallouët, Nonlinear elliptic equations with right hand side
measures, Comm. P.D.E. 17 (1992), 641–655.
[32] P.D. Lax, A.N. Milgram, Parabolic equations, Contributions to the theory of
partial differential equations (L. Bers, S. Bochner, F. John, eds.), Annals of
mathematics studies, vol. 33, Princeton University Press, 1954, 167–190.
[33] J. Leray, J.-L. Lions, Quelques résultats de Višik sur les problèmes elliptiques
non linéaires par les méthodes de Minty-Browder, Bull. Soc. Math. France 93
(1965), 97–107.
[34] J.-L. Lions, Quelques méthodes de résolution des problèmes aux limites non
linéaires, Dunod et Gauthier-Villars, 1969.
[35] J.-L. Lions, E. Magenes, Problèmes aux limites non homogènes et applications,
Dunod, 1968.
[36] P.-L. Lions, F. Murat, Sur les solutions renormalisées d’équations elliptiques non
linéaires, en préparation.
[37] N.G. Meyers, An Lp -estimate for the gradient of solutions of second order elliptic
divergence equations, Ann. Scuola Norm. Sup. Pisa. Cl. Sci. 17 (1963), 189–206.
[40] L. Orsina, Solvability of linear and semilinear eigenvalue problems with L1 data,
Rend. Sem. Mat. Univ. Padova 90 (1993), 207–238.
[41] A. Prignet, Conditions aux limites non homogènes pour des problèmes elliptiques
avec second membre mesure, à paraı̂tre à Ann. Fac. Sciences de Toulouse.
[47] J. Simon, Compact sets in the space Lp (0, T ; B), Ann. Matematica Pura Appli-
cata (1987), 65–96.
[49] J.L. Vazquez, Entropy solutions and the uniqueness problem for nonlinear
second-order elliptic equations, Nonlinear partial differential equations (A. Ben
Kirane, J.P. Gossez, eds.), vol. 343, Addison-Wesley Longman, 1996, 179–203.