Hachana Dounia
Hachana Dounia
Hachana Dounia
DÉPARTEMENT DE MATHÉMATIQUES
MASTER en Mathématiques
Option : Analyse
Par
HACHANA DOUNIA
Titre :
,
Résolution des EDP dordre 2 linéaire :
(Méthode de séparation des variables)
Membres du Comité d’Examen :
Dr. BERBICHE Mohamed UMKB Président
juin 2021
Dédicace
i
REMERCIEMENTS
En premier lieu, je remercie mes parents car grâce à eux que j’ai pu continuer
mes études à un niveau si avancé.
Je remercie mon encadreur GUIDADD Derradji qui m’a guidé dans mon projet
et m’aider à trouver des solutions pour avancer. Ainsi que les professeurs,
BERBICHE MOHAMED et SOLTANI Sihame qui m’ont accepté de
présider les jurys de soutenance.
En …n, je remercie toutes les personnes, amis, qui directement ou indirectement ont
contribué à la réalisation de ce travail.
ii
Table des matières
Remerciements ii
Introduction 1
2.3 Application . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
iii
Table des matières
3 Examples 29
Bibliographie 39
iv
Table des …gures
v
Liste des tableaux
vi
Introduction
1
Introduction
2
Chapitre 1
Une équation aux dérivées partielles (EDP) est une relation reliant une fonction in-
connue de plusieurs variables u à ses dérivées partielles. Les EDPs se trouvent dans
les applications de la, physique, l’ingénierie, la biologie, l’économie, etc. En e¤et,
dans ces domaines, les phénomènes se modélisent souvent par des systèmes mathé-
matiques impliquant des EDPs. Les di¤érents processus du phénomène se décrivent
en déterminant une relation entre u et ses dérivées partielles.
L’exemple le plus simple est lorsque la fonction u dépend uniquement d’une seule
variable. Alors cette relation est décrite simplement par ce qu’on appelle une EDO.
3
Chapitre 1. Généralités sur les équations ordinaires et les équations aux dérivées
partielles
où x est dans R (ou par fois dans un intervalle de R) et y ! (y0 ; :::; yn ) est dans
Rn+1 .
Exemple 1.1.1 Le mouvement d’un objet sur une droite peut être décrit par l’équa-
tion
00
u (x) = g(u(x))
soit u = u (x; y; :::) une fonction de plusieurs variables indépendantes au nombre …ni.
Une EDP pour la fonction u est une relation qui lie :
4
Chapitre 1. Généralités sur les équations ordinaires et les équations aux dérivées
partielles
@u @u @u
F x; y; :::u; ; ; ; ::: =0 (1.1)
@x @y @x2
@u @u @u
F x; y; :::u; ; ; ::: =0
@x @y @x2
.est satisfaite pour x; y; ::: appartenant à une certaine région de l’éspace des variable
independantes.
Dé…nition 1.1.2 une équation aux dérivées partielles(EDP) est une relation fonc-
tionnelle entre une fonction inconnue u de plusieurs variable (u = u(x1 ; :::xn )) est
ses dérivées (x1 ; :::; xn ) 2 D,où D désigne un ouvert de Rn ; (n 2):
Dé…nition 1.1.3 on appelle ordre d’une EDP l’ordre le plus éléve des dérivées par-
tielles intervenant dans l’EDP.
Exemple 1.1.2
@u @u
+ =0 1er ordre:
@x @y
@u @u
+ =0 1er ordre linéaire.
@x @y
@u @u
+ + sin u = 0 1er ordre non-linéaire.
@x @y
@u @2u
+u 2 =0 2eme ordre non-linéaire.
@x @y
5
Chapitre 1. Généralités sur les équations ordinaires et les équations aux dérivées
partielles
9
>
u1 solution =
u1 + u2 est solution.
u2 solution >
;
Equation de transport :
@u
@t
+c @u
@x
=0 où u (x; t) est utilisée pour modéliser la polution de l’air, la dispersion
des colorants ou même le ‡ux de tra…c.
Equation de Burgers :
@u
@t
+ u @u
@x
= 0 où u (x; t) est issue de l’étude de la dynamique des gazs.
Equation de la chaleur :
@u @2u @2u @2u
@t
k @x2
+ @y 2
+ @z 2
= 0 où u (x; y; z; t) est utilisée dans l’étude de la conduc-
tion thermique
Equation d’Euler-bernoulli :
@2u 4
@t2
+ c4 @@xu4 = 0 où u(x,t) est utilisée dans la theorie des poutres dans la prochaine
sectionon illustre la maniére de dérivation de quelques EDPs citées ci-dessus en
utilisant des lois physique
6
Chapitre 1. Généralités sur les équations ordinaires et les équations aux dérivées
partielles
@u @u
A(x; y) + B(x; y) + C(x; y)u = D(x; y) (1.2)
@x @y
est la forme la plus générale pour une EDP linéaire 1er ordre.
@u
Exemple 1.1.4 @x
+ y @u
@y
= 0 avec u = u(x; y) est une EDP 1er ordre, linéaire
homogène.
8
>
< du = @u @u
@x
dx + @y
dy
>
: = ( ydx + dy) @u
@y
u(x; y) = f (ye x )
où f 2 C 1 (R)
Dé…nition 1.1.5 une équation aux dérivées partielles du 2eme ordre de deux va-
riables x et y est une équation de la forme :
dont un cas particulier important est celui ou A,B,C,D et F sont des constantes.
7
Chapitre 1. Généralités sur les équations ordinaires et les équations aux dérivées
partielles
Equation de la place :
En dimension n, on cherche u = u(x1 ; :::xn ) tel que
8
>
< u = f dans un ouvert :
>
: u = y dans @ :
Equation des ondes :
pour un domaine
8 Rn ; on cherche u = u(x1 ; :::xN ; t) telque
>
> @2u
u = f dans [0; T ]
>
> @t2
>
>
>
< u = g sur [0; T ]
>
> u(x; t = 0) = u0 (x)
>
>
>
>
>
: @u
@t
(x; t = 0) = v0 (x)
1. Conditions initiales :
Si u est une fonction de (x; t) 2 R R?+ on donne u (x; t0 ) = 0 (x) ou
D u (x; t0 ) = 0 (x), on parle aussi des conditions de Cauchy.
Dé…nition 1.1.6 Une condition aux limites est une contrainte sur les valeurs que
prennent les solutions des EDP sur une frontiére.Ces conditions imposent une valeur
de u ou de ses derivées au bord du domaine
8
Chapitre 1. Généralités sur les équations ordinaires et les équations aux dérivées
partielles
Exemple 1.1.6 8
>
< u = f sur .
>
: u = d sur @ .
La variable dépendante n’est pas connue sur la frontiére mais si dérivée est bien
dé…nie.
Exemple 1.1.7 8
>
< u = f sur .
>
: @u
= g sur (n normale à @ ).
@n
cette condition est composée des deux première condition imposée sur la frontière.
8
>
>
>
> u = f sur .
<
u = D sur @ 1.
>
>
>
>
: @u
= g sur @ 2
@n
9
Chapitre 1. Généralités sur les équations ordinaires et les équations aux dérivées
partielles
où A, B et C sont constants de x et y.
les équations aux dérivées partielles du 2eme ordre sont classées comme équation du
type elleptique ou du type hyperbolique ou encore du type parabolique, selon la
valeur du determinant : B 2 4AC:
2 2 2
On dit alors que, l’équation A @@xu2 + B @x@y
@ u
+ C @@yu2 + D(x; y; u; @u ; @u ) @u = 0 sur un
@x @y @x
@2u 2
1. @y 2
C 2 @@xu2 = 0 avec C > 0
B2 4AC > 0 ainsi l’équation des ondes est hyperbolique.
@u 2
2. @x
d @@xu2 = 0 ainsi l’équation de la di¤ussion et parabolique.
@2u @2u
3. @x2
+ @y 2
=0
B2 4AC = 4 < 0 ainsi l’équation de la place est elliptique.
2 @2u
4. y @@xu2 @y 2
= 0 e·quation tricomi.
y > 0 )l’EDP est hyperbolique.
y = 0 )l’EDP est parabolique.
y < 0 )l’EDP est elleptique.
Equation elliptique.
10
Chapitre 1. Généralités sur les équations ordinaires et les équations aux dérivées
partielles
le problème elliptique est celui fourni par l’équation de la place (ou de poisson)
soumise à des conditions aux limites, par exemple de Direchlet :
8
>
< u = f dans .
>
: u = u0 sur .
11
Chapitre 2
12
Chapitre 2. Méthode de séparation des variable de fourier
8
>
> 2 2
>
> a1 (x) @@xu2 + b1 (t) @@t2u + a2 (x) @u @x
+ b2 (t) @u @t
+ (a3 (x) + b3 (t)) u (x; t) = h (x; t) ; (x; t) 2 I J
<
u (x; 0) = f (x) ; x 2 I condition initiale , condition aux limites
>
>
>
>
: ou a1 ; b1 ; a2 ; b2 ; a3 ; b3 ; h (x; t) ; f (x) sont des fonctions données.
Conditions aux limites : On distingue di¤érents types des conditions aux limites
(C:L)
IConditions de Dirichlet :
u (a; t) = u (b; t) = 0
IConditions de Neumann :
ux (a; t) = ux (b; t) = 0
IConditions périodiques :
1. on recherche les solutions séparées de (") ces solutions sont de la forme spéciale
13
Chapitre 2. Méthode de séparation des variable de fourier
2. On résout les problémes (p1 ) et (p2 ), les solutions de (p1 ) nous permettent de
construire une base hilbertiènne (Xi )i2N on désigne par (Ti )i2N les solutions de
(p2 ):
4. dans la derniére étape,on calcule les coe¢ cient de cette série et étudie sa conver-
gence.
Principe de superposition
Si u1 et u2 satisfont une EDP linéaire homogène, alors une combinaison linéaire
arbitraire c1 u1 + c2 u2 ; c1 ; c2 2 R satisfait également la même équation.
Dans la suite on va rappeler quelques notions essentielles sur les séries de fourier
qui seront utilisées dans le reste de ce chapitre.
Exemple 2.2.1 On va voir que les fonctions cos et sin peuvent former une base or-
n o
thonormale. Pour cela, l’ensemble des fonctions p12 ; cos(x) sin(x) cos(2 ) sin(2x)
p ; p ; p ; p ; :::; cos(nx)
p ; sin(nx)
p ; :::
un système orthogonale de fonction sur l’intervalle [ ; ] pour la fonction de poids
14
Chapitre 2. Méthode de séparation des variable de fourier
si m; n 0, on a : 8
>
>
>
> 0 si m 6= n
R R <
cos((m n)x) cos((m+n)x)
sin(mx) cos(nx)dx = 2
dx= si m = n
>
>
>
>
: 2 si m = n = 0
si m 0 et n 1, on a :
R R cos((m+n)x) cos((m n)x)
cos(mx) sin(nx)dx = 2
dx =0
Le but de cette section est d’ecrire une fonction f (x) continue par morceaux et 2
périodique sous la forme
a0 X
1
f (x) = + (an cos (nx) + bn sin (nx))
2 n=1
Dé…nition 2.2.3 Les co¢ cients de fourier associées à la fonction f sont données
par :
1
R 1
R 1
R
a0 = f (x)dx; an = f (x) cos(nx)dx; bn = f (x) sin(nx)dx si n 1
15
Chapitre 2. Méthode de séparation des variable de fourier
et
1
R 2( 1)n+1
bn = (x + x2 ) sin(nx)dx = n
Remarque 2.2.1 1. Si f est une fonction paire, alors sa série de fourie est :
Z
a0 X
1
2
+ an cos(nx) avec an = f (x) cos(nx)dx pour n 0:
2 n=1
X
1 Z
2
bn sin(nx) avec bn = f (x) sin(nx)dx pour n 0:
n=1
2.3 Application
16
Chapitre 2. Méthode de séparation des variable de fourier
@u @2u
k = 0; 0 < x < L; t > 0; (2.1)
@t @x2
et la condition initiale
17
Chapitre 2. Méthode de séparation des variable de fourier
Par di¤érentiation de (2.4) par rapport à t et deux fois par rapport à x et par
substitution dans (2.1), obtient
on peut réecrire
T 0 (t) kX 00 (x)
= ; 0 < x < L; t > 0: (2.5)
T (t) X(x)
Puisque x et t sont des variables indépendantes, cette relation implique qu’il existe
une constante ( appellée constante de séparation) telle que
T 0 (t) kX 00 (x)
= = ; 0 < x < L; t > 0: (2.6)
T (t) X(x)
8
>
< u (0; t0 ) = X (0) T (t0 ) = 0
=) X (0) = X (L) = 0
>
: u (L; t ) = X (L) T (t ) = 0
0 0
8
>
< @2X
@x
+ X = 0; 0 x L
(2.7)
>
: X (0) = X (L) = 0
Et
dT
+ T = 0; t 0 (2.8)
dt
18
Chapitre 2. Méthode de séparation des variable de fourier
triviale de (2.7) est appelée fonction propre avec la valeur propre en distingne 3
cas.
2
Cas 01 : = < 0 alors :
8
>
< + =0
>
: exp( L) + exp( L) = 0
Cas 02 : Si = 0; on obtient:
X (x) = + x
19
Chapitre 2. Méthode de séparation des variable de fourier
8
>
< =0
>
: sin( L) = 0
L=n ; = (n L)2 ; n 2 Z
Il en résulte que
n = (n )2
Xn (x) = n sin(n x L)
sont les fonctions caractéristiques du probléme (2,7) comme sin( x) = sin(x) pour
tout x 2 R, il su¢ t donc de considérer
À la …n de cette étape, on peut considérer qu’on a bien construit une base hilbertiènne
(Xi )i2N
20
Chapitre 2. Méthode de séparation des variable de fourier
X
N
Un (x; t) = n sin(n x L) exp( k(n L)2 )t
n=1
de solution séparées est également une solution de l’équation de la chaleur qui satisfait
les condition aux limites de direchlet
X
N
f (x) = n sin(n x L)
n=1
c’est à dire qu’il s’agit d’une combinaison linéaire des fonctions propres ensuite, une
solution au problème de chaleur (2.1)-(2.3) est donnée par :
X
N
Un (x; t) = n sin(n x L) exp( k(n L))2 t
n=1
L’idée brillante de fourier était possible de représenter une fonction arbitraire f satis-
faisant les conditions aux limites (2.2) comme une combinaison linéaire in…nie unique
des fonctions propres sin(n x L): En d’autres termes, il est possible de trouver des
constantes n telles que
X
1
f (x) = n sin(n x L)
n=1
une telle série est appelée une série (ou extension) de fourier (généralisée) de la
21
Chapitre 2. Méthode de séparation des variable de fourier
X
1
u(x; t) = n sin(n x L) exp( k(n L)2 )t
n=1
n .Remarquons
8
Z Z >
<
L X
1 L 0; m 6= n
sin(m x L)f (x)dx = n sin(m x L) sin(n x L)dx =
0 0 >
: L 2; m = n
n=1
RL Z L
0
sin(m x L)f (x)dx 2
n = RL = sin(n x L)f (x)dx
sin2 (m x L)ds L 0
0
ou
>
>
>
> RL
: n = 2
sin(m x L)f (x)dx
L 0
22
Chapitre 2. Méthode de séparation des variable de fourier
@u @2u
k = 0; 0 < x < L; t > 0: (2.9)
@t @x2
ux (0; t) = ux (L; t) = 0; t 0: (2.10)
où f est une condition initiale donnée est k est une constante positive. A…n de rendre
(2.10) consistant avec (2.11) ou suppose la condition de compatibilité : fx (0) =
fx (L) = 0
XTt = KXxx T
Tt Xxx
=K = (2.13)
T X
est une constante réelle. Comme on cherche des solutions ne s’annule pas identi-
23
Chapitre 2. Méthode de séparation des variable de fourier
8
>
< X 00 + X = 0,0 < x < L
(2.14)
>
: X 0 (0) = X 0 (L) = 0
et
fT 0 + KT = 0; t > 0 (2.15)
Cas 01 :
2
= < 0; alors X(x) = exp( x) + exp( x) et X0(x) = [ exp( x) +
exp( x)]
Cas 02 :
si = 0; X(x) = + x et X0(x) =
24
Chapitre 2. Méthode de séparation des variable de fourier
8
>
< X 0 (0) = 0
) =0
>
: X 0 (L) = 0
Donc = 0 est une valeur propre avec la fonction propre X(x) = (constante).
Cas 03 :
2
Si = > 0, alors X(x) = cos( x)+ sin( x) et X 0 (x) = [ sin( x) + cos( x)]
) = 0 et sin( x) = 0
comme sin( l) = 0; alors sin( L)= 0: Conséquemment L = n et = (n nL)2
avec n 2 Z on obtient n = ( nL )2 et Xn (x) = n cos( nLx )
Parceque cos (x) = cos ( x)pour tout x 2 R, et on prend en considération la valeur
n 2 n x
n =( ) et Xn (x) = n cos( ):n2N
L L
On passe maintenant à l’équation (2.8) dont la solution générale est donnée par
T (t) = n exp( k(n L)2 )t; n 2 N
un (x; t) = n cos(n x L) : n 2 N
25
Chapitre 2. Méthode de séparation des variable de fourier
0
X
N
f (x) = + n cos(n x L)
2 n=1
0
X
1
u(x; t) = + n cos(n x L) exp( k(n L)2 )t
2 n=1
On sait que si m; n 0: On a
8
>
>
Z Z >
> 0; m 6= n
L L
cos((m n)x) + cos((m + n)x)dx <
cos(mx) cos(nx)dx = = L 2; m = 0
0 0 2 >
>
>
>
: L; m = n = 0
Pour m=0 multipliant f par cos (0 x L)et on integre sur (0; L) on obtient ;
Z L Z L X1 Z L
0 0L
cos(0 x L)f (x)dx = cos(0 x L)dx+ n cos(0 x L) cos(n x L)dx = +0
0 2 0 n=1 0 2
RL 2
RL
Car 0 cos(0 x L) cos(n x L)dx = 0 puisque n 6= (m = 0) ; alors : 0 = L 0
cos(0 x L)f (x)dx
R
2 L
L 0
f (x)dx
26
Chapitre 2. Méthode de séparation des variable de fourier
Z L Z L X
1 Z L
0
cos(m x L)f (x)dx = cos(m x L)dx + n cos(m x L) cos(n x L)dx
0 2 0 n=1 0
Z L
0L
= [sin (0) sin (m )] + m cos(m x L) cos(n x L)dx
2m 0
L
= m
2
ce qui donne
R
2 L
m = L 0 cos(m x L)f (x)dx
RL P RL
et f (x)) = L1 0 f (x)dx + L2 1
n=1 cos(n x L) 0 cos(n x L)f (x) dx
Z Z
2X
L 1 L
1 2
u(x; t) = f (x)dx+ cos(n x L) exp( k(n L) )t cos(n x L)f (x) dx
L 0 L n=1 0
Exemple 2.3.1
Z 0
2 2 1 2 1 3
0 = x + x2 dx = x + x
0 2 3 0
Si n 0 alors :
Z 0 Z 0
2 2 1 2 1
n = x+x cos (nx) dx = x+x sin (nx) 0
(1 + 2x) sin (nx) dx
0 n n 0
1
= f[ (1 + 2x) cos (nx)]0 g
n2
1
= 2 [ (1 + 2x) cos (nx) + 1]
n
1
= 2 [(1 + 2 ) ( 1)n 1]
n
27
Chapitre 2. Méthode de séparation des variable de fourier
Donc
1 2 X
1
1
f (x) = + 2
+ 2
((1 + 2 ) ( 1)n 1) cos(n x L)
2 3 n=1
n
X
1
1
u(x; t) = 0 + 2
((1 + 2 ) ( 1)n 1) cos(n x L) exp( k(n L)2 )t
n=1
n
28
Chapitre 3
Examples
8
>
> @2u 2
>
> @t2
= x2 @@xu2 + 3x @u
@x
<
u(1; t) = u(e; t) = 0 (3.1)
>
>
>
>
: u(x; 0) = 0; @u (x; 0) = f (x)
@t
où f est une fonction continue sur ]1,e[. Résoudre ce problème (3.1) par la
méthode de séparation des variables, en recherchant la solution sous la forme
P
u(x; t) = k2N Xk (x)Tk (t):
00 00 0
XT = x2 X T + 3xX T
00 0 00
x2 X +3xX T
Ce qui donne X
= T
= 2R
29
Chapitre 3. Examples
X(1) = X(e) = 0
0
dz 0
0 Z
X = Z =
dx x
00
2 2 0
00 dz 0 dz 00 Z Z
X = 2Z + Z = 2+ 2
dx dx x x
00 0 00 0
x2 X + 3xX X = Z + 2Z Z
4.
Z (0) = Z (1) = 0
1
Z (0) = A + B = 0 et Z (1) = e Ae! + Be !
=0
30
Chapitre 3. Examples
z
r= 1 i!:Ce qui conduit à Z (z) = e (A cos (!z) + B sin (!z)) :ou A et B
sont deux constantes réelles. Les conditions aux limites imposent alors :
1
Z (0) = A = 0 et Z (1) = e (B sin [!]) = 0
Pour obtenir une solution non triviale on doit donc avoir sin (!) = 0 ce qui
implique ! = k :k 2 N Il en resulte :
z
Zk (z) = Bk e sin (k z)
et
ln x
Xk (x) = Zk (z) = Zk (ln x) = Bk e sin (k ln x)
p
Les vecteurs sont normalisés en choisissant Bk = 2: On injectant la forme
P
k2N Xk Tk dans l’éqution on obtient :
X X 00 0
X
Xk Tk = x2 Xk + 3xXk Tk = k Xk Tk
k2NF k2NF k2NF
X 00
Tk k Tk Xk = 0
k2NF
00 00
Tk k Tk = Tk + 1 + (k )2 Tk = 0
de solution generale :
q 0 q
Tk (0)
Tk (1) = Tk (0) cos t 1 + (k )2 +q sin t 1 + (k )2
1 + (k )2
31
Chapitre 3. Examples
@u X 0
(x; 0) = Tk (0) Xk (x) = f (x)
@t F k2N
Finalement :
q
X sin (k ln x) sin t 1 + (k )2
u(x; t) = ck
k2N
x
ou Re
1
sin (k ln x) f (y) dy sin (t)
u(x; t)ck = 2 q
1 + (k )2
Exemple 3.0.3 Résoudre le probléme suivant par la méthode de séparation des va-
riables : 8
>
> 2
>
> a @@xu2 + bu + @u
@t
= 0; (x; t) 2]0; `[ R+
<
u(0; t) = u0 e t ; u(`; t) = u0 e t
>
>
>
>
: limt!+1 u(x; t) = 0
2
où a: b: : : u0 et ` des réels striquetement positifs tel que a > b`2 :
t x x t
U (x; t) = u0 e 1 + e
` `
@2 @ @2U @U
a 2
+b + = a 2
+ bU + (3.2)
@x @t @x @t
t x x t
= u0 (b )e 1 + (b ) e
` `
32
Chapitre 3. Examples
(0; t) = (`; t) = 0
lim = lim u U =0
t!+1 t!+1
00 00
T X
b = a =
T X
q
les conditions aux limites conduit à Xn = 2` sin n x
`
: On injecte en…n la série
P
n2N Xn Tn dans l’EDP
X an2 2
0 x x
t t
b Tn + Tn = u0 (b )e 1 + (b ) e
k2N
` ` `
On projette alors le second membre sur la base : Il faut donc calculer les projections
x
de x ! 1 et x ! `
r Z ` p
2 x 2`
sin(n )dx = (1 ( 1)n )
` 0 ` n
p Z ` p
2 x 2`
p sin(n )dx = ( 1)n+1
` ` 0 ` n
p
an2 2
0 2` t x x t
b Tn + Tn = u0 (b )e 1 + (b ) e
`2 n ` `
an2 2
Par hypothèse : b `2
< 0, ce qui donnerait une exponentielle croissante en
solution de problème homogène, incompatible avec l’hypothèse de solution nulle a
33
Chapitre 3. Examples
l’in…ni. On cherche donc une solution particulière qui tend vers 0 en l’in…ni et donc :
p !
an2 2
0 2` (b ) t n+1 (b ) t
Tn (t) = b Tn +Tn = u0 an2 2 e + ( 1) an2 2 e
`2 n b `2
b `2
ce qui conduit à
!
X 2 (b ) (b ) x
(x; t) = u0 an2 2 e t
+ ( 1)n+1 an2 2 e t
sin n :
n2N
n b `2
b `2
`
Etant donné que le second membre est C 1 , on a convergence uniforme des séries
spatiales. Par ailleurs les jTn j sont bornés et décroissants (pour n croissant) il y
a donc convergence uniforme de la série de la fonction de (x; t) sur tout segment
temporel.
Exemple 3.0.4 Soit la température u(x; t) dans un tige de longueur ` qui est parfai-
tement isolée thermiquement incluant les extrémités à x = 0 et x = `: Aucune chaleur
ne s’échappe des extrémités de la tige. Mathématiquement cette condition signi…e que
@u @u
(0; t) = 0; (`; t) = 0 8t > 0:
@x @x
supposons que la température initiale est u(x; 0) = f (x) pour tout x 2 [0; `] : En
d’autres mots, la température u(x; t) est une solution du problème suivant :
@u @u2
= c2 2
@t @x
34
Chapitre 3. Examples
@u @u
(0; t) = 0; (`; t) = 0; 8t > 0: et u(x; 0) = f (x) 8x 2 [0; `] :
@x @x
ou (3.3)
>
>
>
> R R`
: a0 = 2 ` f (x)dx et an = 2
f (x) cos( n ` x )dx:
` 0 ` 0
@u @u2
= c2 2 où u = u(x; t); 0 x `; t > 0
@t @x
`; t > 0
@u 2
@t
= c2 @u
@x2
où u = u(x; t); 0 x
@u @u
(0; t) = 0 et (`; t) = 0 pour tout t 0
@x @x
1. Nous cherchons à déterminer des solutions de ce dernier problème qui sont de la forme
0 00
u(x; t) = X(x)T (t): En substituant dans L’EDP , nous obtenons XT = x2 X T ou
35
Chapitre 3. Examples
00 0
X estla dérivée seconde de X par rapport à x et T est la dérivée de T par rapport
à t. En divisant les deux cotés de l’équation par c2 XT,nousobtenons
0 00
1 T X
c2 T
= X
Le terme de gauche est une fonction de t et celui de droite une fonction de x. Pour
que ceci soit possible, il faut que ces deux termes soient constants. Donc
0 00
1T X
2
= = où est une constante
c T X
Nous pouvons aussi tenir compte des conditions au bord. Nous obtenons ainsi :
@u 0 0
(0; t) = 0 ) X (0)T (t) = 0 ) X (0) = 0 et
@x
@u 0 0
(`; t) = 0 ) X (`)T (t) = 0 ) X (`) = 0
@x
8
>
< X 00 0 0
X = 0 avec X (0) = 0 et X (`) = 0
>
: T
0
c2 T = 0
2
Il nous faut maintenant considérer les di¤érents cas de si 0; disons = avec
00
0; alors la solution générale de X X = 0 est X(x) = A exp( x)+B exp( x)
0
parce que X (x) = A exp( x) B exp( x) et on considérant les deux conditions
0 0
X (0) = 0 et X (`) = 0 , nous obtenons le système d’équations
36
Chapitre 3. Examples
0 10 1 0 1
B CB A C B 0 C
@ A@ A=@ A
exp( `) exp( `) B 0
0 1
B C 2 2
det @ A= (exp( `) exp( `)) = 2 sinh( `) 6= 0
exp( `) exp( `)
Solution 3.0.2
0 10 1 0 1
B 0 CB A C B 0 C
@ A@ A=@ A
sin( `) cos( `) B 0
n 2 n
= et X (x) = A cos x ou n 2 N; n 1
` `
37
Chapitre 3. Examples
00 0
En conclusion pour que l’équation X X = 0 avec les conditions X (0) = 0
0 n 2
et X (`) = 0 ait une solution non nulle, il faut que soit = `
et X (x) =
n
A cos `
x .
0
Si nous considérons maintenant l’équation T c2 T = 0 pour ces di¤érents : Dans
0
le cas où = 0; alors la solution générale T = 0 et T (t) D où D est une constante.
n 2
Dans le cas ou = n = `
ou n 2 N; n 1: Alors la solution générale de
0 cn 2
T c2 T = 0 est T (t) = D exp `
t : où D est une constante.
Nous pouvons conclure que si = 0 alors u(x; t) = X(x)T (t) AD est une solu-
n 2
tion du problème intermédiaire ( ),alors que si = n = `
; alors u(x; t) =
n cn 2
X(x)T (t) = A cos `
x exp `
t est aussi une solution du problème inter-
médiaire ( ) : Parceque l’équation de la chaleur est linéaire ; nous pouvons additionner
ces solutions. Donc la solution formelle du problème intermédiaire ( ) est
a0 X
1
n x cn 2
u(x; t) = + an cos( ) exp( ( ) t)
2 n=1
` `
a0 X
1
n x
u(x; 0) = f (x) = + an cos( )
2 n=1
`
C’est-à-dire que les coe¢ cients an pour n 0 sont les co¢ cients de la série de fou-
rier paire de f(x). En d’autre mots,
Z `
2 n x
an = f (x) cos( )dx pour tout n 0
` 0 `
f (x) = 3 x2 2x3
38
Bibliographie
Z
2
a0 = (3 x2 2x3 )dx = 3
` 0
et, si n 1 alors
Z
2
an = (3 x2 2x3 ) cos(nx)dx:
0
(1 + ( 1)n+1 24
= :
n4
3 X1
(1 + ( 1)n+1 )24
u(x; t) = + 4
cos(nx) exp( n2 t)
2 n=1
n
39
Bibliographie
[2] Bezard, M. (1996). Equations aux dérivées partielles. Ecole nationale supérieure
des techniques avancées.
·
[3] KELLECHE, A. (2019). Equations de la physique mathématique.
[5] Vladmir Arnold. (1997). Leçons sur les équations aux dérivées partielles.
[6] Zitouni, S., & Zennir, K. (2018). Equations aux Dérivées Partielles de nature
physique : EDPs de nature physique Editions universitaires européennes.
40
41
Annexe B : Abréviations et
Notations
@ : la frontière de .
C1 : classe C 1
EDO : Equation di¤érentielle ordinaire.
sh : Sinus hyperbolique.
42
Résumé
Nous avons traité dans ce mémoire les équations
aux dérivées linéaires de second ordre. Nous avons
discuté sur la méthode de séparation des variables
pour obtenir la solution.
Abstract
In this work, we have dealt with second order
linear derivative equation we discussed about the
method of separating variables to get the solution.
ملخص
في هذه المذكرة تطرقنا إلى المعادالت ذات مشتقات جزئية الخطية
من الدرجة الثانية حيث تطرقنا إلى كيفية إيجاد الحل بطريقة فصل
.المتغيرات