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

101
Outro ponto a observar é que as curvas das FI`s, principalmente UDS e QUICK4, são
bastante próximas, resultado que não aparecia na Fig.4.5, onde há uma distinção bastante
clara no comportamento de cada curva das FI`s.
Fazendo agora o número de Pe=0,1 (Fig. 4.18), pode-se ver, analogamente ao já
analisado anteriormente, que as FI`s atingem as mesmas ordens de Pe=0,01 e Pe=5, porém,
chegam nesta ordem após o quinto passo de refinamento, ou seja, após refinar a malha inicial
pela quinta vez.
Módulo do Erro de Discretização
0,1
1E -6
1E -11
1E -16
1E -21
1E -26
E h-U D S
E m er-U D S
E h-C D S
E m er-C D S
E h-Q U IC K 4
E m er-Q U IC K 4
1E -31
1E -8 1E-7 1E -6 1E -5 1E -4 1E -3 0,01 0,1 1
h
Figura 4.18 - Gráfico “ módulos dos erros de discretização (com Pe=0,1) x h” com variação das FI’s para
a variável T c
Este comportamento continua ocorrendo quando tem-se Pe=1,0 (Fig. 4.19) e
Pe=10,0 (Fig. 4.20), sendo que, para atingir as mesmas magnitudes de erros ditadas
anteriormente, teve-se que refinar as malhas 7 e 8 vezes respectivamente.
Vale ainda notar, que à medida que aumenta o número de Peclet as curvas das FI`s
que pareciam ter comportamentos, e inclinações muito próximas até agora, em Pe=10 já
apresentam comportamentos bastante característicos, tendo inclinações diferentes.
E finalmente, analisando o comportamento das FI`s para o Pe=100, a Fig. 4.21
mostra comportamento bastante distinto do que foi visto anteriormente, principalmente no que
diz respeito às magnitudes de erro obtidas. Ao contrário das magnitudes atingidas até aqui,
que se aproximavam de 10 -30 , neste caso específico conseguiu-se valores de 10 -44 , ou seja,