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.
118<br />
APÊNDICE B. Coeficientes e termos fontes das equações governantes<br />
As equações <strong>de</strong> Poisson, advecção-difusão e Burgers foram discretizadas pelo método<br />
dos volumes finitos, utilizando as funções <strong>de</strong> interpolação CDS-2 e o tipo <strong>de</strong> condição <strong>de</strong><br />
contorno <strong>de</strong> Dirichlet, como mostra a Eq. (B.1). Os coeficientes e termos fontes das três<br />
equações foram obtidos a partir da equação algébrica linearizada dada pela Eq. (B.2).<br />
( x = 0 ) φa<br />
e φ ( x L) = φb<br />
φ =<br />
= (B.1)<br />
a φ = a φ + a φ + b<br />
(B.2)<br />
p<br />
P<br />
w<br />
W<br />
e<br />
E<br />
p<br />
EQUAÇÃO DE POISSON<br />
A equação <strong>de</strong> Poisson é dada pela Eq. (B.3) e seu termo fonte ( S ) pela Eq. (B.4). Os<br />
coeficientes e termos fontes foram obtidos com as quatro formas <strong>de</strong> aplicar as condições <strong>de</strong><br />
contorno: sem volume fictício, com volume fictício, com meio-volume e com volume <strong>de</strong><br />
espessura zero. A forma <strong>de</strong> aplicar as condições <strong>de</strong> contorno com volume fictício foi dividida<br />
em duas etapas: primeiro, o termo fonte integrado analiticamente e em seguida, o termo fonte<br />
integrado pela regra do retângulo.<br />
2<br />
d φ<br />
+ S = 0<br />
2<br />
dx<br />
(B.3)<br />
2 xC<br />
( C e )<br />
C<br />
( e −1)<br />
S = −<br />
(B.4)<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: