Etape 4 : prog. principal

pour
$L=2, \;  nx=10L+1,\; x_ j=(j-1)h$,   $ h=\frac{L}{nx-1}$,  $k_{c}=1$ et

$q_ o=0$ avec $rc=\rho c_ v=1$, $f(x)=0$ et $T_ o(x)=\left\{ \begin{array}{ll}0& \mbox{si } 0<x<L/2\\ 6&  \mbox{si } L/2 \le x <L\end{array} \right.$.
Choisir différents $\theta $ et faire varier le pas de temps.
Comparer la solution atteinte au bout d’une longue période avec la solution statique. Pour cela,

  1. modifier les données du programme principal du TP2-1 pour prendre en compte celles de cette application,

  2. entrer les paramètres du schéma en temps,

  3. appeler les sous-programmes de résolution statique ($\color{mycolor}\to T_ s$) et instationnaire ($\color{mycolor}\to T_ e$).

  4. tracer les résultats en 3D ($\color{mycolor}x,t,T(x,t)$) et en 2D dynamique.

 
Besoin d’aide sur ces points ?