fabiana de fátima giacomini - PG-Mec Programa de Pós-Graduação ...

94
Tabela 5.3: Malha, tamanho e número de volumes de controle empregados com a forma de
aplicar as condições de contorno com volume de espessura zero
M
∆ x nodal
N
1 0,33333 5
2 0,11111 11
3 0,03704 29
4 0,01235 83
5 0,00412 245
6 0,00137 731
7 4,57247E-4 2.189
8 1,52416E-4 6.563
9 5,08053E-5 19.685
10 1,69351E-5 59.051
11 5,64503E-6 177.149
12 1,88168E-6 531.443
13 6,27225E-7 1.594.325
14 2,09075E-7 4.782.971
15 6,96917E-8 14.348.909
16 2,32306E-8 43.046.723
0,1
0,01
1E-3
1E-4
Módulo do Erro Numérico
1E-5
1E-7
1E-9
1E-11
1E-13
1E-15
T nod
(Volume Fictício)
T nod
(Sem Volume Fictício)
T nod
(Meio Volume no Contorno)
T nod
(Volume Espessura Zero)
Módulo do Erro Numérico
1E-6
1E-8
1E-10
1E-12
1E-14
1E-16
T m, ret
(Volume Fictício)
T m, ret
(Sem Volume Fictício)
T m, ret
(Meio Volume no Contorno)
T m, ret
(Volume Espessura Zero)
1E-8 1E-7 1E-6 1E-5 1E-4 1E-3 0,01 0,1 1
∆x
1E-8 1E-7 1E-6 1E-5 1E-4 1E-3 0,01 0,1 1
∆x
Figura 5.1: Comparação do erro da variável
T entre as quatro formas
de interesse ( )
nod
de aplicar as condições de contorno
Figura 5.2: Comparação do erro da variável
T ,
entre as quatro formas
de interesse ( )
m ret
de aplicar as condições de contorno
Pela análise das Figs. 5.1 a 5.4, percebe-se que as curvas do erro de discretização
tiveram um decaimento linear, além de estarem qualitativamente iguais nas Figs. 5.1, 5.2 e
5.3. Fazendo um estudo preliminar dos coeficientes, esperou-se que as formas de aplicar as
condições de contorno sem e com volume fictício e com volume de espessura zero tivessem
os mesmos valores dos erros de discretização, pois conforme mostra o apêndice D, os valores
resultantes dos coeficientes foram iguais. Isso não ocorre com a forma de aplicar as condições