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.
52<br />
or<strong>de</strong>m. Consi<strong>de</strong>rando a Eq.(2.58) po<strong>de</strong>-se dizer que a or<strong>de</strong>m do erro é <strong>de</strong> or<strong>de</strong>m p L , também<br />
conhecida por or<strong>de</strong>m assintótica do erro <strong>de</strong> truncamento.<br />
No caso ex<strong>em</strong>plificado acima (UDS) resulta que, para esta aproximação do valor <strong>de</strong><br />
<br />
e<br />
, seu erro <strong>de</strong> truncamento é <strong>de</strong> primeira or<strong>de</strong>m.<br />
Mostrou-se então que o erro <strong>de</strong> truncamento <strong>de</strong> uma função representada por uma<br />
série <strong>de</strong> Taylor t<strong>em</strong> seu erro <strong>de</strong> truncamento <strong>de</strong>finido como a diferença entre seu valor<br />
analítico no ponto “x” e seu valor numérico neste mesmo ponto.<br />
Até agora pô<strong>de</strong>-se ver o erro <strong>de</strong> truncamento como uma operação <strong>de</strong> truncamento<br />
propriamente dita <strong>em</strong> uma série representativa <strong>de</strong> <strong>de</strong>terminada variável, <strong>em</strong> <strong>de</strong>terminado<br />
ponto nodal. É necessário agora ampliar este conceito para o erro <strong>de</strong> truncamento gerado pela<br />
substituição <strong>de</strong> uma ou mais <strong>de</strong>stas funções, representadas por séries <strong>de</strong> Taylor, numa EDP<br />
que está sendo discretizada.<br />
Para isto será consi<strong>de</strong>rada a inserção <strong>de</strong> duas aproximações <strong>em</strong> <strong>de</strong>terminada EDP <strong>em</strong><br />
seu processo <strong>de</strong> discretização. Cada uma <strong>de</strong>stas aproximações irá inserir na equação<br />
discretizada, um erro <strong>de</strong> truncamento local.<br />
Levando <strong>em</strong> consi<strong>de</strong>ração a Eq.(2.6) representativa <strong>de</strong> um processo <strong>de</strong> discretização<br />
ainda não finalizado, e consi<strong>de</strong>rando ainda nesta relação que seu termo fonte seja nulo, t<strong>em</strong>-se<br />
sua redução a,<br />
d<br />
<br />
dx<br />
e<br />
d<br />
<br />
dx<br />
w<br />
0<br />
(2.59)<br />
Adotando agora, como aproximações para os operadores diferencias, a aproximação<br />
por diferenças centrais dadas pelas Eqs.(2.34) e (2.35), para as faces leste e oeste<br />
respectivamente.<br />
Consi<strong>de</strong>rando ainda que cada uma <strong>de</strong>stas aproximações irá gerar um erro <strong>de</strong><br />
truncamento, conforme po<strong>de</strong>-se ver da relação geral dada pela Eq.(2.53), obt<strong>em</strong>-se então, s<strong>em</strong><br />
necessida<strong>de</strong> <strong>de</strong> <strong>de</strong>monstração no momento, as seguintes relações,<br />
<br />
x<br />
e<br />
<br />
<br />
E<br />
<br />
h<br />
P<br />
1<br />
<br />
(2.60)<br />
<br />
x<br />
w<br />
<br />
<br />
P<br />
<br />
h<br />
W<br />
2<br />
<br />
(2.61)