M2 GSI Mercredi 3 octobre 2007
Correction : séparation de variables
Durée : 2h00
Exercice 1
Une barre de longueur L, initialement à température nulle, est chauffée par
une source de chaleur P . Alors que ses deux extrêmités sont maintenues à
température nulle. On cherche à déterminer l’évolution de température dans la
barre, solution de
∂2u
∂u
∂t − λ ∂x2 = P
pour (x, t) ∈]0, L[×]0, +∞[
u(x, 0) = 0 pour x ∈]0, L[
u(0, t) = u(L, t) = 0 pour t ∈]0, +∞[
où λ est une constante thermique donnée.
1. Déterminer la température à l’équilibre, notée v(x).
Nous cherchons une température u(x, t) = v(x) indépendante du temps, ce
qui remplacé dans l’équation donne : −λv 00 (x) = P avec comme condition
aux limites v(0) = V (L) = 0. On intègre deux fois et on trouve v(x) =
P
2λ x(L − x).
2. On pose w(x, t) = u(x, t) − v(x). Montrer que w vérifie la même équation
que u avec un second membre nul, les mêmes conditions en x = 0 et x = L
mais une condition initiale non nulle connue.
∂2w ∂2u 00
On a ∂w ∂u
∂t − λ ∂x2 = ∂t − λ ∂x2 + λv = 0. D’autre part u et v étant nulles en
x = 0 et x = L, il en est de même pour w. Enfin la condition initiale sur
P
w est w(x, 0) = u(x, 0) − v(x) = −v(x) = − 2λ x(L − x).
3. Appliquer la méthode de séparation des variables pour résoudre l’équation
en w, en fonction d’une condition initiale w0 (x).
Nous avons donc mis l’équation sous la forme de celle de l’exercice 1 de
la feuille 2, mais avec une condition de température nulle imposée au
lieu d’une condition de flux nul. On applique maintenant la méthode de
séparation des variables.
Etape 1. On cherche une solution sous la forme w(x, t) = f (t)g(x). On
trouve comme équations différentielles
f 0 (t)g(x) − λf (t)g 00 (x) = 0
soit
f 0 (t) g 00 (x)
=λ .
f (t) g(x)
M2 GSI Mercredi 3 octobre 2007
Les deux membres de l’équation ci-dessus étant exprimés en des variables
indépendantes, les deux expressions sont en fait constantes et égales à une
constante α ∈ R. Donc
α
f 0 (t) = αf (t) et g 00 (x) − g(x) = 0.
λ
D’autre part d’après les conditions aux limites on a g(0) = g(L) = 0.
Comme on l’a vu en TD, pour qu’il existe g non identiquement nul il faut
choisir α < 0, de sorte que l’équation en g ait des solutions en sinus et
cosinus qui peuvent s’annuler deux fois. Posons donc αλ = −ω 2 , ce qui
donne comme solution
g(x) = A sin ωx + B cos ωx.
Comme g(0) = 0 on a B = 0 et g(L) = 0 donne A sin ωL = 0. Si on veut
une solution non identiquement nulle, il faut que A 6= 0, donc sin ωL = 0.
∗
Donc ω = kπL , k ∈ N . Cela donne comme solution
kπx
gk (x) = ak sin , ak ∈ R.
L
L’équation en f quant à elle est une équation du premier ordre qui s’intègre
facilement, en
kπ 2
f (t) = b e−λ( L ) t ,
k k b ∈ R. k
Finalement nous obtenons comme solution à variables séparables
−λ( kπ kπx
L ) t
2
wk (t) = ck e , sinck ∈ R.
L
Notons que toute combinaison de ces solutions et aussi solution donc on a
trouvé en général comme solution une fonction de la forme
∞
X 2
L ) t
−λ( kπ kπx
w(x, t) = ck e sin , ck ∈ R.
L
k=1
Il faut maintenant réaliser la condition initiale w(x, 0) = w0 (x). Pour cela
on cherche les coefficients ck tels que
∞
X kπx
ck sin = w0 (x),
L
k=1
On a déjà fait ce calcul en multipliant par l’un des sinus et en intégrant
sur [0, L] on trouve
2 L
Z
kπx
ck = w0 (x) sin dx,
L 0 L
ce qui termine la résolution du problème en w.
M2 GSI Mercredi 3 octobre 2007
4. w0 étant connue, calculer explicitement en fonction de k et des paramètres
du problème, les coefficients intervenant dans le développement en série de
la solution.
P
On a w0 (x) = − 2λ x(L − x) donc on peut calculer les coefficients ck ex-
plicitement en intégrant deux fois par parties l’intégrale
Z L
P kπx
ck = x(L − x) sin dx.
λL 0 L
On trouve
L Z L
P L kπx P kπx
ck = −x(L − x) cos + (L − 2x) cos dx
λL kπ L 0 λkπ 0 L
où le terme entre crochets est nul. En intégrant à nouveau par parties le
terme entre crochet comporte un sinus qui s’annule et on trouve
Z L
2P L kπx
ck = sin dx.
λ(kπ)2 0 L
Enfin
2P L2
ck = 3
(1 − (−1)k )
λ(kπ)
ce qui donne une expression explicite des coefficients ck , et en particulier
ceux de rang pair sont nuls.
5. En déduire l’expression du développement en série de u. Pouvait-on appli-
quer la méthode de séparation des variables directement sur l’équation de
départ ? Justifiez votre réponse.
P
On a u(x, t) = w(x, t) + 2λ x(L − x) ce qui donne une expression de u à
partir de celle de w. On ne pouvait pas appliquer directement la méthode
à l’équation de départ, car sans second membre nul, on ne peut pas séparer
les variables x et t comme on l’a fait.