18.07.2014 Views

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 ...

SHOW MORE
SHOW LESS

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!