Maple R. 4

16 Dez 98

Solução analítica algébrica da temperatura (T) em x=0.5
tendo como parâmetro o número (n) de dígitos usado e
o número de Peclet (Pe).


Expressão genérica em função do Peclet (p):

a := (exp(p*(1/2))-1)/(exp(p)-1);


                              exp(1/2 p) - 1
                         a := --------------
                                exp(p) - 1


Comandos no Maple para Pe=0.1 e n=8:
a := (exp((1/10)*(1/2))-1)/(exp(1/10)-1);
evalf(a,8);

   Pe = 0.1
------------
   n   T
------------
   4   0.4857
   8   0.48750272
  16   0.4875026035157875
  32   0.48750260351578965633849866787034
  50   0.48750260351578965633849866787055286974894192361805


Comandos no Maple para Pe=1 e n=8:
a := (exp((1)*(1/2))-1)/(exp(1)-1);
evalf(a,8);

   Pe = 1
------------
   n   T
------------
   4   0.3778
   8   0.37754069
  16   0.3775406687981454
  32   0.37754066879814543536109943425450
  50   0.37754066879814543536109943425449152124672063469105


Comandos no Maple para Pe=10 e n=8:
a := (exp((10)*(1/2))-1)/(exp(10)-1);
evalf(a,8);

   Pe = 10
------------
   n   T
------------
   4   6.691e-3
   8   6.6928509e-3
  16   6.692850924284854e-3
  32   6.6928509242848555593619803813252e-3
  50   6.6928509242848555593619803813251803937440089072965e-3


Comandos no Maple para Pe=100 e n=8:
a := (exp((100)*(1/2))-1)/(exp(100)-1);
evalf(a,8);

   Pe = 100
--------------
   n    T
--------------
   4    1.929e-22
   8    1.929e-22
  16    1.928749847963918e-22
  32    1.9287498479639177830169708089294e-22
  50    1.9287498479639177830169708089294104911565367534008e-22