fabiana de fátima giacomini - PG-Mec Programa de Pós-Graduação ...
fabiana de fátima giacomini - PG-Mec Programa de Pós-Graduação ...
fabiana de fátima giacomini - PG-Mec Programa de Pós-Graduação ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
123<br />
EQUAÇÃO DE ADVECÇÃO-DIFUSÃO<br />
A equação <strong>de</strong> advecção-difusão é dada pela Eq. (B.42). Os coeficientes e termos<br />
fontes foram obtidos com as quatro formas <strong>de</strong> aplicar as condições <strong>de</strong> contorno: sem volume<br />
fictício, com volume fictício, com meio-volume e com volume <strong>de</strong> espessura zero.<br />
d<br />
Pe<br />
dx<br />
d<br />
dx<br />
2<br />
φ φ =<br />
2<br />
(B.42)<br />
Condições <strong>de</strong> Contorno Aplicadas Sem Volume Fictício<br />
Os coeficientes e termos fontes para a forma <strong>de</strong> aplicar as condições <strong>de</strong> contorno sem<br />
volume fictício são dados por:<br />
• para o volume <strong>de</strong> controle P = 1:<br />
a = 0<br />
(B.43)<br />
w<br />
a e<br />
= 1 −<br />
Pe<br />
(B.44)<br />
∆x<br />
2<br />
a p<br />
=<br />
Pe 3<br />
+<br />
2 ∆x<br />
(B.45)<br />
b<br />
p<br />
⎛ 2 ⎞<br />
= ⎜ Pe + ⎟φ<br />
a<br />
⎝ ∆x<br />
⎠<br />
(B.46)<br />
• para os volumes <strong>de</strong> controle entre P = 2 a P = N −1:<br />
1 Pe<br />
a w<br />
= +<br />
(B.47)<br />
∆x<br />
2<br />
a e<br />
= 1 −<br />
Pe<br />
(B.48)<br />
∆x<br />
2<br />
a = a + a<br />
(B.49)<br />
p<br />
w<br />
e<br />
b = 0<br />
(B.50)<br />
p