Bonjour,
J'aimerai réaliser un programme en fortran 90 me calculant le température qu'il fait à chaque instant dans une pièce en fonction des conditions aux limites imposées par la température des mur et des chauffages.
Je suis sensé utiliser la méthode Gauss-Seidel pour y arriver.
La formule de la chaleur est un Laplacien de la forme :
dT²/dx² + dT²/dy² = 0
en discrétisant cette équation en 5 points, cela donne :
T(i+1,j) + T(i-1,j) + B² T(i,j+1) + B² T(i,j-1) - C T(i,j) =0
avec, B= DELTAx/DELTAy
et, C= 2( 1+B²)
La pièce est rectangulaire de dimensions 15m en x et 20 m en y. La température des mur est de 5°C et celle des chauffage est de 45°C.
Le 1er chauffage est situé en x=0m et 5m<y<10m
Le second en x=20m et 5m<y<10m
Le 3ème en y=0m et 7m<x<12m
Le 4ème en y=15m et 7m<x<12m
Si qqn pourrait me programmer la méthode, ce serai très sympa.
Merci pour votre aide !
Littletwist
