c fev 06 àj0207 0509
Solution semblable de Lévêque
Nouvelle version de calcul avec FreeFEM++.
Enoncé et corrigé de la PC ENSTA
Fichiers souces: un cas Lévêque
exécuter, puis tracer....:
./FreeFem++-CoCoa neu07.edp
pour Pe =1000 exemple de tracé de ligne isotempérature
tracé de la solution semblable analytique ( notebook mathematica )
tracé de la solution numérique (avec gnuplot):
p[0:3][0:1]'leveq1.txt','leveq10.txt' ,'leveq100.txt','leveq1000.txt'w l,'leveq10000.txt' w l,1-0.538*x
p[0:3][0:1] 'leveq1000.txt',''u 3:4,''u 5:6,'leveq10000.txt',''u 3:4,''u 5:6,''u 7:8,1-0.53*x,'gammainc.txt' w l
tester l'influence du nombre d'Eckert (il varie de 0.1 1 et 10)
