TP3 - Aides : 3. 3 - Iniitalisation et boucle en temps

TP3 - Aides : Etape 3 : schéma en temps : 3. 3 - Iniitalisation et boucle en temps

3. 3 - Iniitalisation et boucle en temps

Initialiser la 1ère colonne du tableau $\color{mycolor}{T}$ à l’aide des valeurs $\color{mycolor}{To(x_ i)}$ :
ces valeurs sont stockées en colonne dans T : T( : ,1) = To ;

Faire une boucle sur $\color{mycolor}{n}$ (temps) dans laquelle on résout le système linéaire et on stocke la solution $\color{mycolor}{(n+1)}$ dans la colonne $\color{mycolor}{(n+2)}$ du tableau $\color{mycolor}{T}$ :

dans cette boucle (de n=1 à nt-1), il faut :

construire le second membre : F=Fc+B*T(:,n) ;
modifier le second membre pour tenir compte des CL de Dirichlet (FCLdi du TP1)
résoudre le système linéaire $A\, T=LU \, T=F$ : T(:,n+1)=lusolve(A,F)

Sous scilab : écrire en complétant les zones grises du fichier chaleur1DP1temps.sci ,

fin de l’aide sur l’étape 3

$\includegraphics[width=9cm]{./prog-images/chaleur1DP1tempstrou.png}$