Discussion :

[niangcoumba]

slt
Je cherche à résoudre l'équation de diffusion de la chaleur( avec le terme source) dans un champs 2D en utilisant la methode des différence finies (sous matlab)...
Si vous avez quelque chose merci de me le passer svp
Merci

[Jerome Briot]

Quel est la question ?
Où se situe le problème ?

Montre nous ce que tu as déjà commencé à coder par toi-même.

[FR119492]

Salut!
Si tu sais comment fonctionne la méthode des différences finies, il te suffit de la programmer, ce qui n'est pas compliqué.
Jean-Marc Blanc

[autocadNUL]

Lu
En général, l'équation de la chaleur c'est :

Laplacien(T)=S (stationnaire)

soit, en discrétisant (dévellopement de Taylor) au moyen des différences finies :

T(i+1)+T(i-1)-2T(i)=S*(dx)² (1D pas constant)

aux limites, tu as

T(i+1)-T(i)=flux*dx (pour les températures imposées, tu peux prendre un coef de convection h=10^30)

Tu te retrouve donc avec un système linéaire : AT=B
Ou A la matrice de rigidité, tridiagonales, B le vecteur des flux, et T les températures. Il te suffit de résoudre ça avec Gauss-Siedel (très rapide)

Pour la 2D c'est pareil,mais avec plus de terme