TP1 - Partie 2 - Aides : 3. 2 - Le programme principal

TP1 - Partie 2 - Aides : Etape 3 : application : 3. 2 - Le programme principal

3. 2 - Le programme principal

Structure

par rapport à la partie 1, point à ajouter en bleu, points à modifier en cyan

- chargement des sous-programmes : exec(...) ; ...

- affectation des constantes du domaine : $L$ ,

- affectation des constantes physiques : $k1$ , $k2$

- affectation des paramètres de la résolution : $nx$ ,

- discrétisation du domaine : $x$ ,

- discrétisation de $kc$ ,

- définition des CL dirichlet dans $cldi$ ,

- définition des CL de flux dans $clne$ ,

- appel du sous-programme de résolution, de solution $T$ ,

- discrétisation de la solution exacte $Texa$ ,

- tracé de $T$ et $Texa$ pour comparaison.

Enfin
Ecrire ce programme principal et comparer la solution obtenue avec la solution exacte.

Visualisation de la température dans le barreau (légende à droite):

$\includegraphics[width=7cm,height=2.5cm]{./prog-images/tp1-2-couleur.png}$

Comparaison de la température calculée et de la solution exacte : $\includegraphics[width=7cm,height=2.5cm]{./prog-images/tp1-2-courbe.png}$

Remarque : la solution calculée est exacte, c’est une propriété de cette méthode en 1D lorsque le terme du à la source est calculé exactement (nul ici), la conductivité étant constante par élément.

fin de l’aide sur l’étape 3