QUATIONS DIFFRENTIELLES GEN-0135

41.

RSOLUTION DE LQUATION DE LA CHALEUR PAR SPARATION DES VARIABLES

Introduction
( 2 ) est loprateur Laplacien scalaire

41.1.

Lquation de la chaleur est la seconde quation en importance dans les connaissances dun ingnieur:

2u 2u 2u u K = c 2 2u = c 2 2 + 2 + 2 , avec c 2 = x t y z

K est la conductivit thermique, la chaleur spcifique et est la densit volumique du corps. Notez la faon dcrire ci-dessus du Laplacien scalaire dans le systme de coordonnes cartsiennes.
Si on suppose quon travaille avec une poutre infiniment plus longue qupaisse, on veut ngliger la propagation radiale de la chaleur, ou encore, si on sintresse un fil dont la section est constante et de densit volumique homogne, situ selon laxe x , et aprs passage en coordonnes cylindriques :
u 2u = c2 2 t x

Il sagit de lquation de la chaleur selon une seule dimension. Noubliez que u ( x, t ) est la temprature du corps en un point donn de laxe x un instant donn t . Si nous avons les conditions aux frontires suivantes (lorsque que les deux extrmits sont plongs une temprature de zro degrs):
u (0, t ) = 0, t u ( L, t ) = 0, t et lunique condition initiale : u ( x, 0) = f ( x ) Nous avons tout un ensemble (quation donde + conditions) pour trouver une solution particulire. Pour cela, nous adoptons de nouveau la stratgie suivante : 1) Faire en sorte de simplifier lexpression de la forme solution en lexprimant comme tant le produit de deux fonctions indpendantes u ( x, t ) = v( x).w(t ) . Dterminer une solution v( x) qui satisfasse les conditions aux frontires. Dterminer la forme solution w(t ) en consquence des imposs sur v( x) . Exprimer la solution u ( x, t ) sous forme de srie de Fourier pour quelle satisfasse prsent la condition initiale sur t .

2) 3) 4)

- page 1 -

Chapitre 41

quation de la chaleur sparation des variables

41.2.

Fonction variables spares :

Encore une fois, afin de simplifier la recherche de la forme solution, nous posons la solution comme tant le produit de deux fonctions variables spares : u ( x, t ) = v( x).w(t ) ds lors, les drives partielles premire et seconde respectives sont : 2u x
2

2v x 2

.w(t ) = v ''.w

u w = v( x) = v.w ' t t Et en pratiquant la sparation des variables : u 2u = c2 2 t x v.w ' = c 2 v ''.w v '' 1 w ' = v c2 w

Comme il sagit de deux variables x et t qui sont indpendantes, les expressions dans chaque membre de droite et de gauche ne peuvent tre en tout temps gales que si elles sont gales une constante arbitraire k :
v '' 1 w ' = = k. v c2 w

Ce qui nous permet de sparer en deux quations :

v " k .v = 0 2 w ' kc .w = 0
Cela nous ramne la rsolution de deux quations diffrentielles ordinaires, une du premier ordre, et une du second ordre.

41.3.

Satisfaire les conditions aux limites :

Lquation du second ordre en v( x) :

v " k .v = 0
a pour forme solution :

v( x) = ax + b, si k = 0. v( x) = C1e +
kx

+ C2 e

kx

, si k > 0.

v( x) = A cos k x + B sin k x, si k < 0. Or lintroduction des conditions aux limites :

- page 2 -

quations diffrentielles

Chapitre 41

u (0, t ) = v(0) w(t ) = 0, t v(0) = 0 u ( L, t ) = v( L) w(t ) = 0, t v ( L ) = 0

nous dmontre que les deux premires formes ne peuvent exister, puisque : v(0) = a 0 + b b = 0 v( x) 0, x. puis v( L ) = aL + 0 = 0 a = 0 et : v(0) = C1e + puis v( L ) = C1e +
k .0 k .L

+ C2 e C1e

k .0 k .L

= 0 C1 + C2 = 0 C1 = C2 = 0 v( x) 0, x. = 0 C1 = 0

Il ne nous reste alors que la troisime possibilit, lorsque k < 0 :

v(0) = A cos k .0 + B sin k .0 = 0

{A = 0

B=0 puis v( L) = 0.cos k .L + B.sin k .L = 0 ou sin k .L = 0 Nous devons choisir B 0 , si nous ne voulons pas avoir une solution triviale. De sorte que nous devons avoir, pour ce quon appelle le nombre donde p : sin k .L = 0 do lexpression de la solution pour v( x) : vn ( x) = sin n x, n = 1, 2,3,.... L k .L = n

k =

n . L Le coefficient B est omis ici, mais rapparatra plus tard, intgr dans la solution de w(t ) .

Il y a donc toute une famille de fonctions solutions qui satisfont les conditions (aux limites) sur x .

41.4.

Consquence sur la fonction temporelle

En rsolvant lquation dordre 1 pour trouver w(t ) : w ' kc 2 .w = 0 .

n Qui se retrouve tre, cause de la contrainte : k = < 0 sur la solution v( x) : L cn w '+ w=0 L
2 2

La solution gnrale de cette quation est, selon le nombre n :

- page 3 -

Chapitre 41

quation de la chaleur sparation des variables

wn (t ) = Cn e

cn dt L 0

+0

En consquence nous avons une famille de formes solutions un ( x, t ) :

n un ( x, t ) = vn ( x) wn (t ) = Cn sin L

xe

cn t L ,

n = 1, 2,3,....

Et la solution gnrale u g ( x, t ) sera videmment la combinaison linaire de toutes ces formes solutions. Cest

donc une srie trigonomtrique qui reprsente la solution gnrale cette quation aux drives partielles homogne :
n u g ( x, t ) = u n ( x, t ) = Cn sin x e L n =1 n =1

cn t L

41.5.

Solution particulire de lquation de la chaleur

Il sagit prsent de fixer les coefficients Cn laide de lunique condition initiale : u ( x, 0) = f ( x ) . n u p ( x,0) = Cn sin xe L n =1

cn .0 L

= f ( x)

Les coefficients seront donc calculs selon les formules dEuler pour les fonctions priodiques p = 2 L :

f ( x) =

n =1

Cn sin

n 1 x Cn = L L

+L

f ( x) sin

n 2 x.dx = L L

f ( x)sin
0

n x.dx L

f ( x) doit ncessairement avoir un prolongement impair fictif entre L < x < 0 (daprs son expression en srie trigonomtrique), de sorte que f ( x)sin n x constitue une fonction paire. Et alors : L Cn = 2 L

f ( x) sin
0

n x.dx L

EXEMPLE : TEMPRATURE INITIALE EN TRIANGLE LE LONG DE LA POUTRELLE

Soit trouver la temprature le long de la poutrelle tout instant ultrieur lorsque la distribution originelle en temprature est triangulaire, et que les bouts sont maintenus zro degrs.

L x, si 0 < x < 2 f ( x) = L ( L x), si <x<L 2

- page 4 -

quations diffrentielles

Chapitre 41

SOLUTION :
La solution est telle que formule pour une poutrelle de longueur L quelconque :
n u g ( x, t ) = Cn sin xe L n =1

cn t L

avec les coefficients comme suit : Cn = 2 L

f ( x) sin

L/2 L n 2 n n x.dx = x.sin x.dx + ( L x).sin x.dx L L L L L/2 0

Ce calcul avait dj t effectu par parties au chapitre 37, dont il suffit de remplacer
8h n
2 2

2h = 1 . De sorte que : L

Cn =

sin

n 4L n , avec n impair . = 2 2 sin 2 2 n

Et la forme extensive la solution particulire est :

2 2 2 2 c 3c 5 c 7 c 4L 1 x L t 1 3 x L t 1 5 x L t 1 7 x L t + + ... e e e e u p ( x, t ) = 2 2 sin sin sin sin L L L L 1 32 52 72

Cette fois, il nest pas facile dutiliser une identit trigonomtrique comme prcdemment pour mettre en vidence la forme f * ( x ct ) et f * ( x + ct ) prolonges impaires de f ( x) . Voici toutefois la distribution de la temprature le long de la poutre diffrents instants. En lvaluant numriquement avec un ordinateur et au dveloppement limit ci-dessus pour u p ( x, t ) .

- page 5 -

Chapitre 41

quation de la chaleur sparation des variables

EXEMPLE : DISTRIBUTION INITIALE DE LA TEMPRATURE EN FORME DE SINUS

Soit trouver la distribution de la temprature le long dune poutrelle de cuivre long de L = 80 cm , si la

x o C , et que les bouts (conditions aux frontires) sont maintenus 80 0 oC . Combien de temps a-il fallu pour que le point de temprature maximal tombe 50 oC .
distribution initiale est f ( x) = u p ( x, 0) = 100sin Les donnes sur le cuivre sont : Conductivit thermique du cuivre : K = 0,95 cal / cm / sec/ oC . Chaleur spcifique du cuivre : = 0, 092 cal / gm / oC . Densit volumique du cuivre : = 8,92 gm / cm3 .

SOLUTION :
Dabord un croquis sur la distribution initiale en temprature :

partir des donnes, on peut calculer la constante de propagation : c2 = Et que :

2 c 1,158 3,14 -1 = 0, 001785 sec L = 80 2

0,95 K = = 1,158 cm 2 /sec . 0, 092 8,92

Comme la condition initiale est telle que :

n u p ( x,0) = Cn sin L n =1

cn L .0 x e

= f ( x) = 100sin L

x .

On en dduit par simple inspection que : C1 = 1 et Cn = 0, n = 2,3, 4,... et que lexpression particulire de la solution est :

- page 6 -

quations diffrentielles

Chapitre 41

u p ( x, t ) = 100sin L

x e 0,001785t

le point de la tige o la temprature est maximale est en :

x = L 2
Et qualors :

x=

L . 2

L u p ( , t ) = 100 e 0,001785t 2 Le moment t1 cette temprature sera redescendue 50 oC est : L 50 u p ( , t1 ) = 100 e0,001785t1 = 50 oC t1 = ln / ( 0, 001785 ) = 388 sec . 2 100

EXEMPLE : VITESSE DE DCROISSANCE DE LA TEMPRATURE

Soit rsoudre le mme problme que lexemple prcdent, mais dont lexpression de la distribution initiale 3 x o de temprature a t change pour : f ( x) = u p ( x, 0) = 100sin C . Les conditions aux limites sont inchanges. 80

SOLUTION :
Dabord un croquis de la distribution de la distribution initiale en temprature :

Cette foisci, le seul terme de la srie de Fourier qui reprsente la distribution en temprature est :
n u p ( x, 0) = Cn sin xe L n =1

cn .0 L

3 = f ( x ) = 100sin L

le troisime terme, ce qui implique : C3 = 100 et Cn = 0, n = 1, 2, 4,5, 6,... Et la solution particulire est :

- page 7 -

Chapitre 41

quation de la chaleur sparation des variables

3 x e u p ( x, t ) = 100sin L

3c t L

Cette fois-ci le temps t2 de dcroissance de temprature 50 oC sera : L 50 u p ( , t2 ) = 100 e0,01607t2 = 50 oC t2 = ln / ( 0, 01607 ) 43 sec . 6 100

Voir la suite dans le chapitre 42.

- page 8 -

quations diffrentielles

Chapitre 41

EXERCICES QUESTIONS & PROBLMES

- page 9 -

Chapitre 41

quation de la chaleur sparation des variables

SOLUTIONNAIRE

- page 10 -