République algérienne démocratique et populaire
الجمهورية الجزائرية الديمقراطية الشعبية
Ministère de l'Enseignement supérieur et de la Recherche
scientifique
وزارة التعليم العالي والبحث العلمي
Domaine : Science et Technologie
Filière : Génie mécanique
Spécialité : Energétique M1
TP méthode numérique
Réalisé par :
AGHA Chaima
Année universitaire : 2024/2025
Aperçu du Programme
Le programme fourni est écrit en langage C et résout l'équation de Laplace à l'aide de la méthode
des différences finies. L'équation de Laplace, qui est une équation aux dérivées partielles du
second ordre, est largement utilisée en physique et en ingénierie pour modéliser des phénomènes
comme la conduction thermique stationnaire, l'électrostatique et l'écoulement de fluides.
Pour résoudre l'équation différentielle parabolique suivante :
2 2
∂ u ∂ u
2
+ 2 =0
∂x ∂y
Description du Fonctionnement
1. Initialisation et Conditions aux Limites :
o Une grille 2D de taille NX x NY (10x10 dans ce cas) est initialisée avec des
valeurs initiales de 0,0.
o Les conditions aux limites sont définies comme suit :
Bord supérieur (u [1][0]) : valeur fixée à 100,0.
Bord inférieur (u [1] [NX-1]) : valeur fixée à 0,0.
Bords gauche et droit (u[0][j] et u[NX-1] [j]) : valeurs fixées à 0,0.
2. Méthode Itérative :
o Le programme utilise une méthode itérative pour résoudre l'équation. Chaque
point de la grille u[i][j] est mis à jour en prenant la moyenne des points voisins.
o L'itération se poursuit jusqu'à ce que la différence entre deux itérations successives
(diff) soit inférieure à une tolérance spécifiée (TOL), fixée ici à 1e-5.
3. Vérification de la Convergence :
o À chaque itération, le programme vérifie si la solution a convergé. Une fois la
convergence atteinte, il affiche le nombre d'itérations nécessaires.
4. Affichage des Résultats :
o Après convergence, le programme affiche les valeurs finales de la grille u[i][j]
sous forme de tableau.
1. Programme C pour l'équation de Laplace :
2. Compile and Run in DOSBOX :
3. Résultats :
La solution finale montre la répartition de température dans la grille 10x10. La
température maximale (100,0) se trouve au bord supérieur, tandis que la température
diminue progressivement jusqu'au bord inférieur (0,0). Les bords gauche et droit étant
également à 0,0, cela crée un gradient descendant du haut vers le bas.
Conclusion
Le programme permet de résoudre avec succès l'équation de Laplace pour une grille 2D avec des
conditions aux limites fixées. Cependant, le code nécessite des corrections pour les erreurs de
syntaxe et les parties incomplètes. La méthode itérative utilisée est adaptée pour ce type de
problème, et le programme illustre bien les principes fondamentaux des méthodes numériques
pour les équations aux dérivées partielles.