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