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 ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
157<br />
E o erro <strong>de</strong> truncamento dado por:<br />
N ''<br />
'''<br />
iv<br />
1 <br />
p<br />
<br />
2 p<br />
<br />
4 p 6<br />
( ) ( . h . h . h ...)<br />
(B.7.17)<br />
N 24 1920 322560<br />
p1<br />
T<strong>em</strong>os então que o erro <strong>de</strong> truncamento toma a forma:<br />
''<br />
'''<br />
IV<br />
<br />
2<br />
4<br />
6<br />
( ) . h . h . h .....<br />
(B.7.18)<br />
24 1920 322560<br />
E <strong>de</strong>ixando na forma genérica:<br />
2 4 6<br />
( ) i . h i . h i . h ...<br />
(B.7.19)<br />
1 2 3<br />
<br />
Ou seja, o erro <strong>de</strong> truncamento da solução numérica para a média da variável obtida<br />
com a regra do retângulo, t<strong>em</strong> como or<strong>de</strong>m assintótica e or<strong>de</strong>ns verda<strong>de</strong>iras os valores:<br />
Or<strong>de</strong>ns verda<strong>de</strong>iras p 2,4,6....<br />
(B.7.20)<br />
V<br />
Or<strong>de</strong>m assintótica p 2<br />
(B.7.21)<br />
L<br />
B.7.3 TERCEIRA VARIÁVEL DE INTERESSE<br />
<br />
i<br />
x1<br />
Com a terceira variável <strong>de</strong> interesse ocorre o mesmo fenômeno ocorrido com a<br />
segunda, ou seja, ela sofrerá o impacto do erro <strong>de</strong> truncamento da função <strong>de</strong> interpolação<br />
utilizada na obtenção da solução numérica, assim como o erro <strong>de</strong> truncamento gerado pela<br />
relação numérica utilizada para obtenção da variável, que neste caso está representada pela<br />
Eq.(B.3.3).<br />
Para acharmos o erro <strong>de</strong> truncamento, e consequent<strong>em</strong>ente o erro <strong>de</strong> discretização da<br />
aproximação, para obtenção da <strong>de</strong>rivada primeira da variável <strong>em</strong> x 1, t<strong>em</strong>os que multiplicar<br />
a Eq.(B.5.5) por (-9) e somar o resultado à Eq.(B.5.6). Consi<strong>de</strong>rar<strong>em</strong>os para este caso os<br />
seguintes características:<br />
j e j 1<br />
P j 2 W<br />
Resultando então:<br />
i 8.<br />
<br />
e<br />
9. <br />
P<br />
W<br />
3 iii 2 72 iv 3 234 v 4<br />
<br />
e<br />
<br />
. <br />
e<br />
. h . <br />
e<br />
. h . <br />
e.<br />
h .... (B.7.22)<br />
3. h 8 384 3840