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)





PYL