TP2 - Partie 1 - Aides : 3. 1

TP2 - Partie 1 - Aides : Etape 3 : programme principal : 3. 1 - Les points

3. 1 - Les points

Les données sont les mêmes que celles de la partie 2 su TP1 à part la source : $f(x)\; = \; sin(\pi \, \dfrac {x}{L})$

1. la donnée de la source et sa discrétisation :

une fois la discrétisation du domaine faite, la source discrétisée peut être définie directement.

sous scilab : fs = sin( %pi * x / L ) ;

Remarque : ${\it \color{brown}{fs}}$ est un tableau ligne.

2. l’appel du sous-programme de résolution :

comme pour le test précédent :

sous scilab : T=chaleur1DP1plein(nx,x,kc,cldi,clne,fs’) ;

Remarque : fs est un tableau ligne, on passe le tableau colonne fs’.

3. la solution exacte de ${\color{mycolor}-\frac{d}{dx}(k_{c} \frac{dT}{dx})=sin(\pi \, \frac{x}{L}) }$ sur ${\color{mycolor}[0,L]}$ + CL :

${\color{mycolor}\triangleright }$ trouver $A$ , une primitive de $f(x)=sin(\pi \, \dfrac {x}{L})$ ,
${\color{mycolor}\triangleright }$ sur $[0,L/2]$ , résoudre $-k_{c_1} \frac{dT}{dx}=A$ avec $T(0)=To=0$ ,
${\color{mycolor}\triangleright }$ sur $[L/2,L]$ , résoudre $-k_{c_2} \frac{dT}{dx}=A$ avec $-k_{c_2}\frac{dT}{dx}(L)=q_ o$ ,
${\color{mycolor}\triangleright }$ imposer la continuité en $x=\frac{L}{2}$ .