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 ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
65<br />
10) imprimir e visualizar os resultados referentes a φ<br />
P<br />
( x P<br />
);<br />
11) imprimir e visualizar os resultados finais.<br />
É um processo iterativo que envolve em seu ciclo o cálculo dos coeficientes e termos fontes<br />
<strong>de</strong> todos os volumes <strong>de</strong> controle internos, inteiros e reais, sendo resolvido pelo método direto<br />
TDMA e em seguida, calculado os valores das variáveis <strong>de</strong> interesse.<br />
3.3.1 Condições <strong>de</strong> Contorno Aplicadas Sem Volume Fictício<br />
Conforme a Fig. 2.1 e consi<strong>de</strong>rando os três mo<strong>de</strong>los matemáticos propostos, a malha<br />
utilizada é estruturada e uniforme. Portanto as distâncias entre os centros dos volumes <strong>de</strong><br />
controle são iguais e as distâncias entre as faces também são iguais, implicando que o<br />
domínio é composto por uma geometria <strong>de</strong> nós centrados entre as faces do volume e faces<br />
centradas entre os centros dos volumes. Assim, o cálculo <strong>de</strong>ssa distância uniforme ( ∆ x)<br />
é<br />
realizado por:<br />
L<br />
∆ x =<br />
(3.14)<br />
N<br />
on<strong>de</strong> L é o comprimento do domínio <strong>de</strong> cálculo e N é o número <strong>de</strong> volumes <strong>de</strong> controle.<br />
O refino da malha computacional está <strong>de</strong>finido na seção 2.4 do capítulo 2 e foi<br />
empregado neste trabalho com duas razões: q = 3 para obter os valores da variável em<br />
x = 1/ 2 diretamente no centro dos volumes <strong>de</strong> controle e q = 2 para obter os valores da<br />
variável em x = 1/ 2 fazendo a média aritmética entre dois volumes <strong>de</strong> controle, obtendo seus<br />
valores nas faces dos volumes. Para obter o valor do número dos volumes <strong>de</strong> controle ( N )<br />
com a forma <strong>de</strong> aplicar as condições <strong>de</strong> contorno sem volume fictício foram utilizadas as<br />
seguintes equações:<br />
N = ( 2 ∆x<br />
, 4 ∆x<br />
, 8 ∆x<br />
, K )<br />
(3.15)<br />
N<br />
= ( 3 ∆x<br />
, 9 ∆x<br />
, 27 ∆x<br />
, K )<br />
(3.16)<br />
on<strong>de</strong><br />
∆ x neste caso, indica a quantida<strong>de</strong> <strong>de</strong> volumes e não o cálculo do tamanho <strong>de</strong>sses<br />
volumes que é dado pela Eq. (3.14). A Eq. (3.15) refere-se à razão <strong>de</strong> refino q = 2 que