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.

74<br />

I<br />

principal ( )<br />

Para encontrar a or<strong>de</strong>m do erro <strong>de</strong> truncamento a priori para a <strong>de</strong>rivada da variável<br />

φ com o esquema CDS-2, subtrai-se as Eqs. (3.43) e (3.44). Assim, o primeiro<br />

e<br />

termo do segundo membro representa a aproximação do esquema CDS-2 para o termo<br />

difusivo, dado pela Eq. (2.10). Os outros termos da série representam o valor do erro <strong>de</strong><br />

truncamento, dado pela Eq. (3.48).<br />

( φ −φ<br />

)<br />

I E P 1 III 2 1 V 4 1 VII 6<br />

φ<br />

e<br />

= − φe<br />

∆x<br />

− φe<br />

∆x<br />

− φe<br />

∆x<br />

−K (3.47)<br />

∆x<br />

24 1920 322560<br />

I 1 III 2 1 V 4 1 VII 6<br />

[ φ ] = − φ ∆x<br />

− φ ∆x<br />

− φ ∆x<br />

−K<br />

ε τ e<br />

e<br />

e<br />

e<br />

(3.48)<br />

24 1920 322560<br />

on<strong>de</strong> as or<strong>de</strong>ns verda<strong>de</strong>iras do erro <strong>de</strong> truncamento são dadas por<br />

pV<br />

= 2,4,6, K<br />

, sendo a<br />

or<strong>de</strong>m assintótica pL = 2 .<br />

3.4 SOLUÇÃO NUMÉRICA<br />

A proprieda<strong>de</strong> física consi<strong>de</strong>rada para as sete variáveis <strong>de</strong> interesse estudadas pelas<br />

equações <strong>de</strong> Poisson e advecção-difusão foi a temperatura ( T ) e pela equação <strong>de</strong> Burgers foi a<br />

velocida<strong>de</strong> ( u ). A Tab. 3.4 mostra a simbologia atribuída para as variáveis <strong>de</strong> interesse.<br />

Tabela 3.4: Simbologia atribuída às variáveis <strong>de</strong> interesse<br />

VARIÁVEL DE INTERESSE<br />

SÍMBOLO<br />

1) φ em x = 1 2 (valor nodal) T<br />

nod<br />

u<br />

nod<br />

2) φ em x = 1 2 (média arit.) T<br />

med<br />

u<br />

med<br />

3) φ com a regra do retângulo m ret<br />

T ,<br />

u<br />

m , ret<br />

4) φ com a regra do trapézio m trap<br />

T ,<br />

u<br />

m , trap<br />

5) l<br />

1<br />

do erro <strong>de</strong> discretização<br />

E<br />

m<br />

6) dx<br />

dφ<br />

em x = 0 com DDS I<br />

DDS<br />

7) dx<br />

dφ<br />

em x = 0 com DDS-2 I<br />

DDS−2

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

Saved successfully!

Ooh no, something went wrong!