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.
159<br />
B.7.4 QUARTA VARIÁVEL DE INTERESSE – MÉDIA DA NORMA ( l 1<br />
):<br />
Para acharmos o erro <strong>de</strong> truncamento e consequent<strong>em</strong>ente o erro <strong>de</strong> discretização do<br />
valor da média da norma l , t<strong>em</strong>os consi<strong>de</strong>rar inicialmente que a equação que nos dá esta<br />
1<br />
variável é,<br />
l<br />
1<br />
<br />
1<br />
N<br />
N<br />
<br />
P1<br />
<br />
P<br />
(B.7.31)<br />
P<br />
Consi<strong>de</strong>rando a Eq.(B.7.31) po<strong>de</strong>mos ver que o erro <strong>de</strong> discretização <strong>de</strong>sta variável<br />
<strong>de</strong>pen<strong>de</strong>rá da solução numérica da variável , ou seja, o erro <strong>de</strong> discretização <strong>de</strong>sta variável<br />
P<br />
será dada pelos erros <strong>de</strong>finidos pela variável I. Desta forma, o que foi comentado para a<br />
variável I, <strong>em</strong> relação aos erros <strong>de</strong> truncamento dos termos advectivos e difusivos, vale<br />
também para a variável IV.<br />
B.8 ANÁLISE A PRIORI DAS ORDENS VERDADEIRAS E ASSINTÓTICA DA<br />
VARIÁVEL I<br />
Iniciar<strong>em</strong>os a análise a priori da variável I com a aproximação CDS-2 para os termos<br />
advectivo e difusivo. Para as <strong>de</strong>mais funções <strong>de</strong> interpolação somente indicar<strong>em</strong>os suas<br />
or<strong>de</strong>ns seguindo o mesmo raciocínio do esqu<strong>em</strong>a CDS-2 a seguir.<br />
Finalmente, após analisarmos os erros <strong>de</strong> truncamento envolvidos na obtenção da<br />
solução numérica com as funções <strong>de</strong> interpolação CDS-2, para a <strong>de</strong>rivada primeira na face<br />
(termo difusivo) e para a variável na face (termo advectivo), po<strong>de</strong>mos calcular o erro <strong>de</strong><br />
discretização total para a primeira variável.<br />
Aplicando à Eq.(A.3.3) os erros analisados para as aproximações e representados<br />
pelas Eqs.(B.6.3) e (B.6.13), obt<strong>em</strong>os a Eq.(B.8.1) como segue:<br />
Pe<br />
i<br />
i<br />
i<br />
(<br />
) ( ) e(<br />
) <br />
(<br />
) ( ) e(<br />
) 0<br />
.<br />
CDS 2 j CDS 2<br />
j CDS 2<br />
j CDS 2<br />
j CDS 2<br />
j CDS 2<br />
j<br />
<br />
(B.8.1)<br />
Aplicando ainda o operador diferencial:<br />
i<br />
.<br />
<br />
0<br />
d( ) Pe<br />
CDS2 CDS2<br />
<br />
(B.8.2)<br />
j<br />
Substituindo a Eq.(B.8.2) <strong>em</strong> (B.8.1), t<strong>em</strong>os:<br />
j