PYL: c Août 05 (Lagrée P.-Y.), LMM-Univ PARIS 6, B 162,
4 place Jussieu, 75252 PARIS pyl(a)ccr.jussieu.fr
http://www.lmm.jussieu.fr/~lagree

Résolution numérique de l'équation de la chaleur stationnaire

On veut résoudre l'équation de la chaleur stationnaire 2D: 
     ò2T   ò2T 
     --- + --- = 0
     òx2   òy2
avec comme condition initiale 
   T(x,t=0)=1
et comme conditions aux limites 
   T(0,t)=T(1,t)=0.
avec comme condition initiale 
   T(x,t=0)=1
et comme conditions aux limites 
   T(0,t)=T(1,t)=0.
On fait une relaxation de type Gauss Siedel:
alpha=dy*dy/(2*(dx*dx+dy*dy)); beta=dx*dx/(2*(dx*dx+dy*dy));
Tij=(1-relax)*T[i][j] + relax*(alpha*(T[i+1][j]+T[i-1][j])+beta*(T[i][j+1]+T[i][j-1]));
retour à l'applet Java de la résolution de l'équation de la chaleur, modification des conditions aux limites en direct.
La source pour la solution de l'équation de la chaleur, La source pour l'interface graphique (inspirée de GraphLayout Sun/Apple).