section 8.

Solution de l’équation de la chaleur

Dans cette section, on donnera la solution de l’équation de la chaleur
avec
1 des conditions aux limites homogènes.
2 des conditions aux limites non homogènes.
Équation de la chaleur Problème de Dirichlet Problème de Neumann

1. Équation de la chaleur avec C.L homogènes

L’équation de la chaleur avec des conditions aux limites homogènes est

α2 uxx = ut , pour 0 ≤ x ≤ L, t > 0, (1a)

u(0,t) = 0 et u(L,t) = 0, pour t > 0, (1b)
u(x ,0) = f (x ), pour 0 ≤ x ≤ L. (1c)

Étape 1 : Séparation des variables : en posant u(x ,t) = F (x )G(t), alors G(t) vérifie
l’équation du premier ordre

G ′ (t) − α2 λG(t) = 0,

et F (x ) satisfait au problème avec conditions aux limites homogènes

F ′′ (x ) − λF (x ) = 0, F (0) = 0 et F (L) = 0.

Étape 2 : Calcul des valeurs et fonctions propres : la résolution du problème avec

conditions aux limites en F (x ) donne les valeurs et fonctions propres
2
nπ nπx
  
λn = − , Fn (x ) = cn sin , n = 1,2, · · · .
L L
Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 2 / 17
Équation de la chaleur Problème de Dirichlet Problème de Neumann

Équation de la chaleur avec C.L homogènes (suite)

Étape 3 : Calcul des solutions de l’EDO en G(t)
 nπ 2 αnπ )2 t
Gn′ (t) + α Gn (t) = 0 ⇒ Gn (t) = dn e −( L , n = 1,2, . . . .
L

Étape 4 : Principe de superposition : on cherche la solution de (1) sous la forme

∞ ∞
X X αnπ )2 t nπx
u(x ,t) = Dn Fn (x )Gn (t) = Dn e −( L sin( ).
L
n=1 n=1

Étape 5 : Calcul des coefficients Dn de sorte que la condition initiale (1c) soit satisfaite
∞ Z L
X nπx 2 nπx
u(x ,0) = Dn sin( ) = f (x ) ⇒ Dn = f (x ) sin( )dx ,n =
L L 0
L
n=1
1,2, . . .

La solution de l’équation de la chaleur avec conditions aux limites homogènes est

∞  Z L

X 2 nπx αnπ )2 t nπx
u(x ,t) = f (x ) sin( )dx e −( L sin( )
L 0
L L
n=1

Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 3 / 17

Équation de la chaleur Problème de Dirichlet Problème de Neumann

2. Équation de la chaleur avec C.L non homogènes

L’équation de la chaleur avec des conditions aux limites non homogènes est
2

 α uxx = ut , pour 0 ≤ x ≤ L, t > 0,
u(0,t) = T1 et u(L,t) = T2 , pour t > 0, (2)
u(x ,0) = f (x ), pour 0 ≤ x ≤ L,


la solution du problème (2) est

∞
X αnπ )2 t nπx x
u(x ,t) = Dn e −( L sin( ) + (T2 − T1 ) + T1 ,
L L
n=1
Z Lh
2 x nπx
i  
avec Dn = f (x ) − (T2 − T1 ) − T1 sin dx .
L 0
L L
x
1 où uperm. (x ,t) = (T2 − T1 ) + T1 est la solution stationnaire (ou stable ou
L
permanente) qui vérifie le problème vxx = 0,v (0) = T1 ,v (L) = T2 et
∞
X αnπ )2 t nπx
2 utrans. (x ,t) = Dn e −( L sin( ) est la solution transitoire qui vérifie
L
n=1
l’équation de la chaleur avec des conditions aux limites homogènes (1) et la
x
condition initiale u(x ,0) = f (x ) − (T2 − T1 ) − T1
L
Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 4 / 17
Équation de la chaleur Problème de Dirichlet Problème de Neumann

Exemple 1 :
On considère le problème suivant :

∂v ∂2v
− = 0, pour 0 < x < 30 et t ≥ 0,
∂t ∂x 2
∂v ∂v (3)
(0,t) = 0 et (30,t) = 0, pour t > 0,
∂x ∂x
v (x ,0) = 25[u5 (x ) − u10 (x )] pour 0 ≤ x ≤ 30,

où ua (x ) = u(x − a) est la fonction unité échelon (ou la fonction de Heaviside) au

point a.
(a) On pose v (x ,t) = F (x )G(t). En utilisant la méthode de séparation des
variables, montrer que la fonction F (x ) satisfait au problème de fonctions
propres
F ′′ (x ) − λF (x ) = 0, F ′ (0) = 0 et F ′ (30) = 0 (4)
et la fonction G(t) satisfait à l’équation différentielle G ′ (t) − λG(t) = 0, où
λ est une constante réelle de séparation.
(b) Résoudre le problème de fonctions propres (4).
On demande une solution complète où chaque étape est justifiée.
(c) En vous servant des résultats obtenus en (a) et (b), trouver la solution du
problème (3).
Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 5 / 17
Équation de la chaleur Problème de Dirichlet Problème de Neumann

Exemple 2 :
On considère le problème suivant :
∂c ∂2c
= α2 2 − β 2 c(x ,t), pour 0 < x < a et t ≥ 0, (5a)
∂t ∂x
c(0,t) = 0 et c(a,t) = 0, pour t > 0, (5b)
c(x ,0) = c0 pour 0 ≤ x ≤ a, (5c)

où α,β,a et c0 sont des constantes. Ce problème décrit la diffusion en régime

transitoire de la concentration c(x ,t) d’une certaine substance dans un milieu
donné. La concentration initiale de la substance est donnée par c(x ,0) = c0
(a) Soit c(x ,t) = F (x )G(t). En utilisant la méthode de séparation des variables,
montrer que la fonction F (x ) satisfait l’équation différentielle
α2 F ′′ (x ) − (λ + β 2 )F (x ) = 0, et la fonction G(t) satisfait l’équation
différentielle G ′ (t) − λG(t) = 0, où λ est une constante de séparation.
(b) Déterminer les conditions limites associées à l’équation différentielle satisfaite
par F (x ) et résoudre le problème de fonctions et valeurs propres associé à la
fonction F (x ).
(c) En vous servant des résultats obtenus en (a) et (b), trouver la solution du
problème (5). On demande une solution complète où chaque étape est
justifiée.

Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 6 / 17

Équation de la chaleur Problème de Dirichlet Problème de Neumann

Exemple 3 :
On considère le problème suivant

1 ∂2p ∂p ∂p

 − = , pour 0 < x < 1 et t > 0;
2 ∂x 2 ∂x ∂t

p(0,t) = 0 et p(1,t) = 0 pour t > 0 (6)

p(x ,0) = 1 pour 0 < x < 1.


(a) Soit p(x ,t) = F (x )G(t), en utilisant la méthode de séparation des variables,
montrer que la fonction F (x ) satisfait au problème de fonctions propres
F ′′ (x ) − 2F ′ (x ) − 2λF (x ) = 0, F (0) = 0 et F (1) = 0, et la fonction G(t)
satisfait à l’équation différentielle G ′ (t) − λG(t) = 0, où λ est une constante
de séparation.
(b) Vérifier que les solutions du problèmede fonctions  propres en F sont données
x n2 π 2 + 1
par Fn (x ) = cn e sin(nπx ) et λn = − , où les cn sont des
2
constantes réelles pour n = 1,2,3, . . . .
(c) En vous servant des résultats obtenus en (a) et (b), trouver la solution du
problème (6).
∞
X 2nπ[(−1)n+1 e −1 + 1]
Indice : 1 = bn e x sin(nπx ), où bn = ,n ≥ 1.
n2 π 2 + 1
n=1

Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 7 / 17

 section 8.8

L’équation de Laplace
Dans cette section, nous étudierons
1 Le problème de Dirichlet homogène sur 3 cotés.
2 Les conditions aux limites de Dirichlet générales.
3 Le problème de Neumann.
Équation de la chaleur Problème de Dirichlet Problème de Neumann

1. Le problème de Dirichlet homogène sur 3 cotés

1 L’équation de la chaleur en dimension deux s’écrit ut = α2 (uxx + uyy )
2 À l’état stationnaire, ut = 0, on obtient l’équation de Laplace

uxx + uyy = 0.

3 Le problème de Dirichlet sur un rectangle avec des conditions homogènes

sur trois cotés est
y
b u(x, b) = 0 (a, b)

u(a, y) = f (y)
u(0, y) = 0

∆u = uxx + uyy = 0

0 x
u(x, 0) = 0 a

uxx + uyy = 0, pour 0 ≤ x ≤ a, 0 ≤ y ≤ b, (7a)

u(x ,0) = 0 et u(x ,b) = 0, pour 0 < x < a, (7b)
u(0,y ) = 0, pour 0 ≤ y ≤ b, (7c)
u(a,y ) = f (y ), pour 0 ≤ y ≤ b. (7d)
Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 9 / 17
Équation de la chaleur Problème de Dirichlet Problème de Neumann

La solution du problème de Dirichlet homogène

Étape 1 : Séparation des variables : en posant u(x ,y ) = F (x )G(y ), alors G(y ) vérifie
le problème avec conditions aux limites homogènes

G ′′ (y ) + λG(y ) = 0, G(0) = 0 et G(b) = 0,

et F (x ) vérifie l’équation différentielle ordinaire

F ′′ (x ) − λF (x ) = 0, F (0) = 0

Étape 2 : Calcul des valeurs et fonctions propres : la résolution du problème avec

conditions aux limites en G(y ) donne les valeurs et fonctions propres
2
nπ nπy
  
λn = , Gn (y ) = cn sin , n ≥ 1.
b b

Étape 3 : Calcul des solutions de l’EDO en F (x )

2
nπ nπx
  
F ′′ (x ) − F (x ) = 0, F (0) = 0 ⇒ Fn (x ) = dn sinh , n ≥ 1.
b b

Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 10 / 17

Équation de la chaleur Problème de Dirichlet Problème de Neumann

La solution du problème de Dirichlet homogène (suite)

Étape 4 : Principe de superposition : on cherche la solution de (7) sous la forme
∞ ∞
nπx nπy
X X    
u(x ,y ) = Dn Fn (x )Gn (y ) = Dn sinh sin
b b
n=1 n=1

Étape 5 : Calcul de Cn et Dn pour que la condition aux limites (7d) soit satisfaite
∞
nπa nπy
X    
u(a,y ) = Dn sinh sin = f (y ) ⇒
b b
n=1
Z b
2
 nπy 
Dn = nπa

f (y ) sin dy , n = 1,2, . . . .
b sinh b 0
b

La solution du problème de Dirichlet est alors

∞
" #
Z b
2 nπy nπx nπy
X      
u(x ,y ) = nπa

f (y ) sin dy sinh sin .
n=1
b sinh b 0
b b b

Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 11 / 17

Équation de la chaleur Problème de Dirichlet Problème de Neumann

2. Les conditions aux limites de Dirichlet générales

La solution du problème de Dirichlet ∆u = uxx + uyy = 0, 0 ≤ x ≤ a, 0 ≤ y ≤ b
avec les conditions aux limites générales
y u2(x,b) = f2 (x)

u2(0,y) = 0

u2 (a,y) = 0
∂ 2 u2 ∂ 2 u2
∂x2
+ ∂y 2
=0

x
u2(x,a) = 0
y y y
u3 (x,b) = 0 b u(x,b) = f2(x) (a,b) u4(x,b) = 0

u4(a,y) = g2 (y)
u(0,y) = g1 (y)

u(a,y) = g2 (y)
u3(0,y) = g1 (y)

u1(a,y) = 0

u4(0,y) = 0
∂ 2 u3 ∂ 2 u3 ∂2u ∂2u ∂ 2 u4 ∂ 2 u4
∂x2 + ∂y 2 =0 ∂x2 + ∂y 2 =0 ∂x2 + ∂y 2 =0

x (0,0) x x
u3(x,0) = 0 u(x,0) = f1 (x) a u4 (x,a) = 0
y u1(x,b) = 0
u1(0,y) = 0

u1(a,y) = 0
∂ 2 u1 ∂ 2 u1
∂x2 + ∂y 2 =0

x
u1(x,0) = f1(x)

u(x ,0) = f1 (x ), u(x ,b) = f2 (x ), 0 < x < a,

u(0,y ) = g1 (y ), u(a,y ) = g2 (y ), 0 ≤ y ≤ b.
4
est par superposition, u(x ,y ) = ui (x ,y ) où ui , i = 1,2,3,4, résout un problème
P
i =1
de Dirichlet avec conditions aux limites homogènes sur 3 côtés.
Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 12 / 17
Équation de la chaleur Problème de Dirichlet Problème de Neumann

3. Le problème de Neumann
Les conditions de Neumann consistent à imposer des conditions sur la dérivée
∂u
normale = ∂~n u = ∇u · ~n aux frontières du domaine.
∂~n
y
1. Pour un rectangle, ces conditions sont
uy (x, b) = 0
b (a, b)
∂~n u(0,y ) = ∇u · ~n(0,y ) = −ux (0,y ),0 < y < b,
ux (a, y) = f (y)
−ux (0, y) = 0

∆u = uxx + uyy = 0 ∂~n u(a,y ) = ∇u · ~n(a,y ) = ux (a,y ),0 < y < b,

0 x
∂~n u(x ,0) = ∇u · ~n(x ,0) = −uy (x ,0),0 < x < a,
−uy (x, 0) = 0 a
∂~n u(x ,b) = ∇u · ~n(x ,b) = uy (x ,b),0 < x < a.

2. Le problème de Neumann sur un rectangle avec des conditions homogènes

sur trois des cotés est

∆u = uxx + uyy = 0, pour 0 ≤ x ≤ a, 0 ≤ y ≤ b, (8a)

uy (x ,0) = 0 et uy (x ,b) = 0, pour 0 < x < a, (8b)
ux (0,y ) = 0, pour 0 ≤ y ≤ b, (8c)
ux (a,y ) = f (y ), pour 0 ≤ y ≤ b. (8d)

Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 13 / 17

Équation de la chaleur Problème de Dirichlet Problème de Neumann

Problème de Neumann sur un rectangle

Étape 1 : Séparation des variables : en posant u(x ,y ) = F (x )G(y ), alors G(y ) vérifie
le problème avec conditions aux limites homogènes

G ′′ (y ) + λG(y ) = 0, G ′ (0) = 0 et G ′ (b) = 0,

et F (x ) vérifie l’équation différentielle ordinaire

F ′′ (x ) − λF (x ) = 0, F ′ (0) = 0

Étape 2 : Calcul des valeurs et fonctions propres : la résolution du problème avec

conditions aux limites en G(y ) donne les valeurs et fonctions propres
2
nπ nπy
  
λn = , Gn (y ) = cn cos , n = 0,1,2 . . . .
b b

Étape 3 : Calcul des solutions de l’EDO en F (x )

2
nπ nπx
  
F ′′ (x )− F (x ) = 0, F ′ (0) = 0 ⇒ Fn (x ) = dn cosh , n = 0,1, . . . .
b b

Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 14 / 17

Équation de la chaleur Problème de Dirichlet Problème de Neumann

La solution du problème de Neumann sur un rectangle (suite)

Étape 4 : Principe de superposition : on cherche la solution de (8) sous la forme
∞ ∞
nπx nπy
X X    
u(x ,y ) = Dn Fn (x )Gn (y ) = D0 + Dn cosh cos
b b
n=0 n=1

Étape 5 : Calcul des coefficients Dn pour que la condition initiale (8d) soit satisfaite
∞
nπ nπa nπy
X    
ux (a,y ) = Dn sinh cos = f (y ) ⇒
b b b
n=1
Z b
2
 nπy 
Dn = nπa

f (y ) cos dy , n = 1,2, . . . .
nπ sinh b 0
b

La solution du problème de Neumann est donc

∞
" #
Z b
2 nπy nπx nπy
X      
u(x ,y ) = D0 + nπa

f (y ) cos dy cosh cos .
n=1
nπ sinh b 0
b b b

Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 15 / 17

Exemple 4 :
La distribution de la température u(x ,y ) dans une plaque carrée d’un matériau
composite est donnée par l’équation aux dérivées partielles

∂2u ∂2u ∂u
+ +2 = 0, pour 0 < x < π, 0 < y < π. (9)
∂x 2 ∂y 2 ∂y
La distribution de température satisfait aux conditions suivantes :
( ∂u
u(x ,0) = 0 et (x ,π) = f (x ) pour 0 ≤ x ≤ π,
∂y
u(0,y ) = 0 et u(π,y ) = 0 pour 0 ≤ y ≤ π,
où f (x ) est une fonction donnée.
(a) Soit u(x ,y ) = F (x )G(y ), en utilisant la méthode de séparation des variables,
montrer que la fonction F (x ) satisfait au problème de fonctions propres
F ′′ (x ) − λF (x ) = 0, F (0) = 0 et F (π) = 0, la fonction G(y ) satisfait à
l’équation différentielle G ′′ (y ) + 2G ′ (y ) + λG(y ) = 0, et à la condition limite
homogène G(0) = 0, où λ est une constante réelle.
(b) Vérifier que les solutions du problème de fonctions propres en F sont données
par Fn (x ) = Cn sin(nx ) et λn = −n2 , où les Cn sont des constantes réelles
pour n = 1,2,3, . . . ..
(c) En vous servant des résultats obtenus en (a) et (b), trouver la distribution de
la température u(x ,y ).
On demande une solution complète où chaque étape est justifiée.
Équation de la chaleur Problème de Dirichlet Problème de Neumann

Exemple 5 :
Soit u(x ,y ), la la solution de l’équation de Laplace dans la plaque rectangulaire
0 < x < 2,0 < y < 3 qui satisfait aux conditions aux limites d
∂u ∂u
(
(0,y ) = 0, (2,y ) = 0, si 0 < y < 3;
∂x ∂x (10)
u(x ,0) = 0, u(x ,3) = f (x ) = 2 − x si 0 < x < 2.

(a) Soit u(x ,y ) = F (x )G(y ).

En utilisant la méthode de séparation des variables, montrer que la fonction
F (x ) satisfait au problème de fonctions propres

F ′′ (x ) − λF (x ) = 0, F ′ (0) = 0 et F ′ (2) = 0, (11)

et que G ′′ (y ) + λG(y ) = 0, G(0) = 0, où λ est une constante réelle.

(b) En vous servant des résultats obtenus en (a), trouver u(x ,y ), la solution de
l’équation de Laplace qui satisfait aux conditions aux limites (11).
Note : Les solutions au problème de fonctions propres (11) sont données par
2
nπx nπ
  
Fn (x ) = cn cos et λn = − ,
2 2
où cn sont des constantes réelles pour n = 0,1,2,3, . . . .
Z. Coulibaly GCH2535 - H2018 - Gr. 1 : www.mgi.polymtl.ca/zoumana.coulibaly 24 janvier 2018 17 / 17