18.07.2014 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!