18.07.2014 Views

verificação de funções de interpolação em advecção-difusão 1d ...

verificação de funções de interpolação em advecção-difusão 1d ...

verificação de funções de interpolação em advecção-difusão 1d ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

141<br />

Para o <strong>de</strong>senvolvimento do trabalho, consi<strong>de</strong>rar<strong>em</strong>os para este caso que a equação<br />

geral do erro <strong>de</strong> truncamento dada pela Eq.(B.4.4), terá aqui seus coeficientes dados por “e”,<br />

ficando então:<br />

2 4 6<br />

( ) e . h e . h e . h ...<br />

(B.6.17)<br />

j CDS 2 1 2 3<br />

<br />

Comparando as Eqs.(B.6.16) e (B.6.17), po<strong>de</strong>mos <strong>de</strong>finir os coeficientes como:<br />

e<br />

1<br />

ii<br />

j<br />

<br />

<br />

8<br />

e<br />

2<br />

iv<br />

j<br />

<br />

<br />

384<br />

e<br />

3<br />

vi<br />

j<br />

; etc. (B.6.18)<br />

46080<br />

Sendo as or<strong>de</strong>ns <strong>de</strong> seus erros as seguintes:<br />

Or<strong>de</strong>ns verda<strong>de</strong>iras p 2,4,6....<br />

(B.6.19)<br />

V<br />

Or<strong>de</strong>m assintótica p 2<br />

(B.6.20)<br />

L<br />

B.6.2 UDS / CDS-2<br />

Devido ao fato <strong>de</strong> já termos <strong>de</strong>scrito no ít<strong>em</strong> B.6.1 o processo <strong>de</strong> aproximação e<br />

análise a priori da aproximação do termo difusivo com CDS-2, suprimir<strong>em</strong>os aqui esta<br />

<strong>de</strong>monstração, relacionando resultados do ít<strong>em</strong> já mencionado quando necessário.<br />

Para achar o valor da variável " " na face usando o UDS, simplesmente<br />

reorganizar<strong>em</strong>os a Eq. (B.5.5) <strong>de</strong>ixando na seguinte forma:<br />

2<br />

3<br />

4<br />

5<br />

i h ii h iii h iv h v h<br />

<br />

j<br />

<br />

j 1<br />

<br />

j<br />

. <br />

j<br />

. <br />

j<br />

. <br />

j<br />

. <br />

j<br />

. ....<br />

(B.6.21)<br />

2 8 48 384 3840<br />

Consi<strong>de</strong>rar<strong>em</strong>os o índice “j” para a face leste, “j+1” é o valor do nó “P”. Se for a<br />

oeste, ter<strong>em</strong>os “E”, s<strong>em</strong>pre <strong>de</strong> acordo com a Fig.B.2.1. L<strong>em</strong>brando ainda que, in<strong>de</strong>pen<strong>de</strong>nte<br />

<strong>de</strong> que face consi<strong>de</strong>rarmos, o erro <strong>de</strong> truncamento é o mesmo. Substituindo agora a Eq.(B.4.1)<br />

<strong>em</strong> (B.6.21), t<strong>em</strong>os,<br />

2<br />

3<br />

4<br />

5<br />

<br />

i h ii h iii h iv h v h <br />

<br />

<br />

<br />

<br />

<br />

<br />

j<br />

<br />

j 1<br />

E<br />

j1<br />

j<br />

.<br />

j<br />

.<br />

j<br />

.<br />

j<br />

.<br />

j<br />

. ..... (B.6.22)<br />

2 8 48 384 3840 <br />

Que po<strong>de</strong> ser reescrita da forma:<br />

(B.6.23)<br />

j<br />

( <br />

j<br />

)<br />

UDS<br />

( <br />

j<br />

)<br />

UDS<br />

e ( <br />

j<br />

)<br />

UDS

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

Saved successfully!

Ooh no, something went wrong!