Université des Sciences et de la Technologie 1ère année Master Génie des Systèmes

Houari Boumediene Industriels et Tertiaires

Faculté de génie électrique TD Techniques d’optimisation et EDP
Département d’électrotechnique

TD n°1 : Equations aux dérivées partielles (méthode de différences finies)

Exercice n°1 :
Soit à résoudre l'équation d’électrostatique dont le potentiel V(x, y) est régi par l'EDP suivante :

 2V  2V
+ = 0; x   0,9 , y   0,12
x 2 y 2
Les conditions aux limites de type Dirichlet sont données par :
V ( x, 0 ) = 0v; V ( x,12 ) = 0v; V ( 9, y ) = 0v; V ( 0, y ) = 50v.

1. Donner L’approximation par différence centrale de cette équation.

2. Si on considère que x = y = 3 , donner le schéma maillage de ce domaine.
3. Ecrire les équations nécessaires à la résolution de ce problème (les mettre sous forme matricielle).
4. Déterminer la solution du problème

Exercice n°2 :
Soit à résoudre l'EDP suivante :
2S 2S
− 3 = S ; x   0,10 , y   0,15
x 2 y 2
Les conditions aux limites de type Dirichlet sont données par :
S ( x, 0 ) = 10; S ( x,15 ) = 20; S ( 0, y ) = 0; S (10, y ) = 0.
1. Donner le type et l'ordre de l’EDP.
2. Donner l’approximation par différence centrale de cette équation.
3. Si on considère que x = y = 5 , donner le schéma de maillage de ce domaine.
4. Ecrire les équations nécessaires à la résolution de ce problème.
5. Déterminer la solution du problème

Exercice n°3 :
Soit l’équation de la chaleur non stationnaire définie par l'EDP suivante :
T '  T
2
−k = 0 ; x   0,   , t   0, 0.4
t x 2

Les conditions aux limites de type Dirichlet sont données par :

T ( x, 0 ) = f ( x ) = x ( − x ) ; T ( 0, t ) = 00 C ; T ( , t ) = 00 C.

1. Si on considère que x = et t = 0.2s , donner le schéma de maillage de ce domaine.
3

1
Mme A.KHEITER Master 1 Pro: Génie des Systèmes Industriels et Tertiaires
 2. Donner la discrétisation de cette équation en utilisant le schéma implicite.
3. Ecrire les équations nécessaires à la résolution de ce problème.
4. Trouver les températures Ti l à t=0.4s.

Exercice n°4 :
Soit à résoudre l'EDP suivante :
f ( x, t )  2 f ( x, t )  4
(1 + t ) − = 0 ; x  0,  , t   0, 0.02
t x 2
 3
Les conditions aux limites sont données par :
4 
f ( x, 0 ) = h ( x ) = sin ( 2 x ) + sin ( 4 x ) ; f ( 0, t ) = f  , t  = 1.
3 
1
1. Si on considère que x = et t = 0.01 , donner le schéma de maillage.
3
2. Donner la discrétisation de cette équation en utilisant le schéma explicite.
3. Ecrire les équations nécessaires à la résolution du problème.
4. Déterminer la solution du problème
Exercice n°5 :
Soit l’équation de Poisson suivante :
2 f 2 f  3
= y 2 (  + x ) ; x  0,  , y   0,1
+
x y  2
2 2

Les conditions aux limites de type Dirichlet sont données par :

3 
f ( x, 0 ) = 5; f ( x,1) = 15; f ( 0, y ) = 10; f  , y  = 0.
2 
1. Donner l’approximation par différence centrale de cette équation en considérant un maillage uniforme (
1
x = y = )
2
2. Déduire l’équation nécessaire à la résolution de ce problème et déterminer la solution du problème (une seule
itération).

Exercice n°6 :
Soit à résoudre l'EDP suivante :
 1  g ( x, t ) −  g ( x, t ) = 0 ; x  0, 3  , t  0; 0.03
2

   
 3t  t x
2
 4 
Les conditions aux limites sont données par :
g ( x, 0 ) = cos ( 2 x ) + sin ( x ) ; g ( 0, t ) = 0 ; g  , t  = 2
3
4 
1
5. Si on considère que x = et t = 0.01 , donner le schéma de maillage.
4
6. En utilisant le schéma explicite, donner la discrétisation de cette équation (dans le but de simplifier l’équation,
t
on prend m = 3t ).
x
2

7. Ecrire les équations nécessaires à la résolution du problème pour t1 , t 2 et t 3 . Déterminer alors g ( x1 , t2 ) ,

g ( x2 , t2 ) et g ( x3 , t3 )

2
Mme A.KHEITER Master 1 Pro: Génie des Systèmes Industriels et Tertiaires
Exercice n°7:
L’équation de diffusion de la température en fonction du temps t, dont la température T(x, y) est régit par l'EDP
suivante :

T  2T
− k ' 2 = 0 ; x   0,12 , t   0,8
t x
1. Donner l’ordre de l’équation. Les conditions aux limites de type Dirichlet sont données par :
T ( x, 0 ) = 0 C ; T (12, t ) = 0 C ; T ( 0, t ) = 50 C. .
2. Donner la discrétisation par schéma explicite de cette équation.
3. Si on considère que x = y = 4 , donner le schéma maillage de ce domaine.
4. Ecrire les équations nécessaires à la résolution de ce problème (les mettre sous forme matricielle).
5. Déterminer la solution du problème
Exercice n°8 :
Soit à résoudre l'EDP suivante :

f ( x, t )  2 f ( x, t )  1
(1 + t ) − = 0 ; x  0,  , t   0, 0.01
t x 2
 3
Les conditions aux limites sont données par :

1 
f ( x, 0 ) = h ( x ) = sin ( 2 x ) + sin ( 4 x ) ; f ( 0, t ) = u  , t  .
3 
1
8. Si on considère que x = et x = 0.01 , donner le schéma de maillage.
9
9. Donner la discrétisation de cette équation en utilisant le schéma explicite.
10. Ecrire les équations nécessaires à la résolution du problème.
11. Déterminer alors f ( x1 , t1 ) et f ( x2 , t1 ) .

Exercice n°9:
Soit l’équation 1D définie par :

2 f
− 2*10−6 f = 75*10−6 x ( 75 − x ) avec 0  x  75
x 2

1. Donner la forme matricielle de ce problème.

2. Trouver les solutions yi aux nœuds xi pour un pas h= 25. On donne f ( 0 ) = f ( 75 ) = 0 .

3
Mme A.KHEITER Master 1 Pro: Génie des Systèmes Industriels et Tertiaires