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.
126<br />
EQUAÇÃO DE BURGERS<br />
A equação <strong>de</strong> Burgers é dada pela Eq. (B.63) e seu termo fonte ( S ) pela Eq. (B.64).<br />
Os coeficientes e termos fontes foram obtidos com as quatro formas <strong>de</strong> aplicar as condições<br />
<strong>de</strong> contorno: sem volume fictício, com volume fictício, com meio-volume e com volume <strong>de</strong><br />
espessura zero. O termo fonte ( S ) foi integrado pela regra do retângulo, numericamente é<br />
dado pela Eq. (B.65).<br />
2<br />
dφ<br />
Re<br />
dx<br />
2<br />
d φ<br />
+ S<br />
dx<br />
=<br />
2<br />
(B.63)<br />
S = Re<br />
2<br />
e<br />
x Re<br />
x Re Re<br />
( 2e<br />
− e −1)<br />
Re<br />
( e −1) 2<br />
(B.64)<br />
S<br />
xP<br />
Re Re<br />
( 2e<br />
− e −1)<br />
Re<br />
( e −1) 2<br />
2 xP<br />
Re<br />
P<br />
= Re e<br />
(B.65)<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.66)<br />
w<br />
a = 2<br />
(B.67)<br />
e<br />
a<br />
p<br />
∗<br />
( φ<br />
∗ + φ )<br />
= 6 + Re ∆x<br />
P E<br />
(B.68)<br />
b<br />
p<br />
∗ 2 ∗ 2 ⎤<br />
( φE<br />
− φ ) ⎥⎦<br />
2<br />
⎡ 2 β<br />
= 2S<br />
P∆x<br />
+ 4φ a<br />
+ Re ∆x⎢2φ<br />
a<br />
−<br />
P<br />
(B.69)<br />
⎣ 2<br />
on<strong>de</strong> o símbolo ∗ representa os valores explícitos da iteração anterior.<br />
• para os volumes <strong>de</strong> controle entre P = 2 a P = N −1: