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.
70<br />
( φ )<br />
N<br />
1 ⎧ ⎡ P−1<br />
+ φP<br />
⎤⎫<br />
φ = ⎨∆x∑<br />
⎢ ⎥⎬<br />
(3.31)<br />
L ⎩ P=<br />
2 ⎣ 2 ⎦⎭<br />
on<strong>de</strong> φ<br />
P−1<br />
e φ<br />
P<br />
são as soluções numéricas nos volumes oeste e principal.<br />
• Derivada <strong>de</strong> 1ª or<strong>de</strong>m <strong>de</strong> ( )<br />
φ em x = 0 obtida com DDS:<br />
dφ<br />
=<br />
dx<br />
( − )<br />
φ 2<br />
φ 1<br />
∆x<br />
(3.32)<br />
• Derivada <strong>de</strong> 1ª or<strong>de</strong>m <strong>de</strong> ( )<br />
φ em x = 0 obtida com DDS-2:<br />
dφ<br />
=<br />
dx<br />
4φ<br />
2<br />
− 3φ<br />
1<br />
−φ3<br />
2∆x<br />
(3.33)<br />
3.3.4 Condições <strong>de</strong> Contorno Aplicadas Com Volume <strong>de</strong> Espessura Zero<br />
Devido à forma <strong>de</strong> aplicação da condição <strong>de</strong> contorno com volume <strong>de</strong> espessura zero<br />
sobre o contorno, a equação para calcular a distância uniforme ( ∆ x)<br />
, dada pela Eq. (3.14),<br />
recebe uma nova formulação, sendo adaptada para:<br />
L<br />
∆x =<br />
(3.34)<br />
N<br />
( − 2)<br />
O refino da malha foi realizado como na forma <strong>de</strong> aplicar as condições <strong>de</strong> contorno<br />
sem e com volume fictício, utilizando as razões: q = 3 e q = 2 . Para obter o número dos<br />
volumes <strong>de</strong> controle ( N ) com a forma <strong>de</strong> aplicar as condições <strong>de</strong> contorno com volume <strong>de</strong><br />
espessura zero foram utilizadas as seguintes equações:<br />
N = ( 2 ∆x<br />
+ 2, 4 ∆x + 2 , 8 ∆x + 2 , K )<br />
(3.35)<br />
N = ( 3 ∆x<br />
+ 2 , 9 ∆x + 2 , 27 ∆x + 2, K )<br />
(3.36)<br />
on<strong>de</strong> a Eq. (3.35) refere-se a razão <strong>de</strong> refino q = 2 que fornece um número <strong>de</strong> volumes <strong>de</strong><br />
controle pares e a Eq. (3.36) refere-se a razão <strong>de</strong> refino q = 3 que fornece um número <strong>de</strong><br />
volumes <strong>de</strong> controle ímpares.