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.
40<br />
Existem várias funções <strong>de</strong> interpolação <strong>de</strong>scritas e utilizadas em CFD. Neste trabalho<br />
foram utilizadas as funções <strong>de</strong> interpolação CDS-2 (Central Differencing Scheme)<br />
(FERZIGER e PERIC, 2002) <strong>de</strong> 2ª or<strong>de</strong>m tanto para os termos difusivos quanto para os<br />
termos advectivos e correção adiada para os termos não-lineares (TANNEHILL et al., 1997).<br />
As funções DDS (Downstream Differencing Scheme) e DDS-2 <strong>de</strong> 1ª e 2ª or<strong>de</strong>m para a<br />
<strong>de</strong>rivada da proprieda<strong>de</strong> ( )<br />
φ em x = 0 e as funções <strong>de</strong> integração pela regra do retângulo e<br />
pela regra do trapézio para a média da variável <strong>de</strong> interesse.<br />
2.5.1 Funções <strong>de</strong> Interpolação para Termos Difusivos<br />
O esquema CDS-2 utiliza uma interpolação linear, conduzindo a erros <strong>de</strong> 2ª or<strong>de</strong>m.<br />
Consi<strong>de</strong>rando a disposição da malha na Fig. 2.1, as aproximações das <strong>de</strong>rivadas da<br />
proprieda<strong>de</strong> ( φ ) para as faces ( w ) e ( e ) do volume <strong>de</strong> controle ( )<br />
P , são:<br />
dφ<br />
dx<br />
e<br />
φE<br />
−φP<br />
=<br />
∆xe<br />
(2.10)<br />
dφ<br />
dx<br />
w<br />
φP<br />
−φW<br />
=<br />
∆xw<br />
(2.11)<br />
on<strong>de</strong> a proprieda<strong>de</strong> ( φ ) é avaliada no centro dos volumes <strong>de</strong> controle e as distâncias ( x)<br />
∆ são<br />
calculadas entre as faces.<br />
Os esquemas DDS e DDS-2 empregam funções <strong>de</strong> interpolação linear, conduzindo a<br />
erros <strong>de</strong> 1ª e 2ª or<strong>de</strong>ns, respectivamente. As aproximações para a <strong>de</strong>rivada da variável <strong>de</strong><br />
interesse ( )<br />
φ em = 0<br />
x para a face ( e ) do volume <strong>de</strong> controle ( )<br />
P , são respectivamente:<br />
( −φ<br />
)<br />
dφ 2 φP<br />
=<br />
dx ∆xe<br />
e<br />
c<br />
(2.12)<br />
dφ<br />
dx<br />
e<br />
9φ<br />
P<br />
− φE<br />
− 8φ<br />
c<br />
=<br />
3∆xe<br />
(2.13)<br />
on<strong>de</strong> ( φ<br />
c<br />
) é a condição <strong>de</strong> contorno conhecida.