UNIVERSIDADE FEDERAL DO ESPÍRITO SANTO

DEPARTAMENTO DE ENGENHARIA MECÂNICA

PROGRAMA DE PÓS-GRADUAÇÃO EM ENGENHARIA MECÂNICA

ESTEFÂNIA PINTOR CANZIAN

MÉTODO DOS ELEMENTOS DE CONTORNO

VITÓRIA

2019

1
 SUMÁRIO

1. PROBLEMA PROPOSTO ............................................................................................. 3

2. METODOLOGIA DE DESENVOLVIMENTO ............................................................... 3
3. RESULTADOS E CONCLUSÃO................................................................................... 7

2
1. PROBLEMA PROPOSTO
O problema em análise, condução de calor em uma placa, conta com a seguinte
configuração:

Em que as condições de contorno são:

𝑇(𝑥, 0) = 𝑥(1 − 𝑥)
𝑇(1, 𝑦) = 0
𝑇(0, 𝑦) = 0
{ 𝑇(𝑥, 1) = 0

2. METODOLOGIA DE DESENVOLVIMENTO
Para resolver o problema, na Solução Analítica, foi empregado o método da
separação de variáveis. Considerando:

𝑇(𝑥, 𝑦) = 𝑋𝑛 (𝑥 ) 𝑌𝑛 (𝑦)

Para Laplace:

𝑑2𝑇 𝑑2𝑇
+ =0
𝑑𝑥 2 𝑑𝑦 2

Em que:

𝑋𝑛" (𝑥 ) 𝑌(𝑦) = − 𝑌𝑛" (𝑦) 𝑋𝑛 (𝑥) = 0

𝑋𝑛" (𝑥) 𝑌𝑛" (𝑦)
=− = 𝜆2𝑛
𝑋𝑛 (𝑥) 𝑌𝑛 (𝑦)

Sendo 𝜆 uma constante de separação:

3
 𝜆2𝑛 = 0 ⇒ 𝑒
{ 𝜆2𝑛 > 0 ⇒ 𝐴𝑛 cos + 𝐵𝑛 𝑠𝑒𝑛
𝜆2𝑛 < 0 ⇒ 𝐶𝑛 cosh + 𝐷𝑛 senh

Resulta em:

𝑋 " = 𝐴𝑛 cos(𝜆𝑛 𝑥) + 𝐵𝑛 𝑠𝑒𝑛(𝜆𝑛 𝑥 )

{ "𝑛
𝑌𝑛 = 𝐶𝑛 cosh(𝜆𝑛 𝑦) + 𝐷𝑛 𝑠𝑒𝑛ℎ(𝜆𝑛 𝑦)

Considerando as condições de contorno, temos:

Para 𝑥 = 0 → 𝑇(0, 𝑦) = 0:

𝑋𝑛 (0) = 𝐴𝑛 cos(𝜆𝑛 0) + 𝐵𝑛 𝑠𝑒𝑛(𝜆𝑛 0)

𝐴𝑛 = 0

Para 𝑥 = 1 → 𝑇(1, 𝑦) = 0 e com 𝐴𝑛 = 0:

𝑋𝑛 (1) = 0 cos(𝜆𝑛 1) + 𝐵𝑛 𝑠𝑒𝑛(𝜆𝑛 1)

𝐵𝑛 𝑠𝑒𝑛(𝜆𝑛 1) = 0

Com 𝑠𝑒𝑛(𝜆𝑛 ) = 0, onde 𝜆𝑛 = 𝑛𝜋, (𝑛 = 1, 2 … , 𝑛).

Assim:

𝑋𝑛 (𝑥 ) = 𝐵𝑛 𝑠𝑒𝑛(𝑛𝜋𝑥 )

Para 𝑦 = 1 → 𝑇(𝑥, 1) = 0:

𝐶𝑛 cosh(𝜆𝑛 1) + 𝐷𝑛 𝑠𝑒𝑛ℎ(𝜆𝑛 1) = 0

𝐶𝑛 cosh(𝜆𝑛 1) = −𝐷𝑛 𝑠𝑒𝑛ℎ(𝜆𝑛 1)

𝐶𝑛 = −𝐷𝑛 𝑡𝑔ℎ(𝜆𝑛 )

Para 𝑦 = 0 → 𝑇(𝑥, 0) = 𝑥(1 − 𝑥 ):

𝑇(𝑥, 0) = 𝑋𝑛 (𝑥 ) 𝑌𝑛 (0) = 𝑥 (1 − 𝑥 )

𝐵𝑛 𝑠𝑒𝑛(𝜆𝑛 𝑥){[−𝐷𝑛 𝑡𝑔ℎ(𝜆𝑛 ) cosh(𝜆𝑛 𝑦)] + 𝐷𝑛 𝑠𝑒𝑛ℎ(𝜆𝑛 𝑦)} = 𝑥(1 − 𝑥 )

𝐵𝑛 𝑠𝑒𝑛(𝜆𝑛 𝑥){[−𝐷𝑛 𝑡𝑔ℎ(𝜆𝑛 ) cosh(𝜆𝑛 0)] + 𝐷𝑛 𝑠𝑒𝑛ℎ(𝜆𝑛 0)} = 𝑥(1 − 𝑥 )

4
 𝐵𝑛 𝑠𝑒𝑛(𝜆𝑛 𝑥 ){[−𝐷𝑛 𝑡𝑔ℎ(𝜆𝑛 )]} = 𝑥(1 − 𝑥 )

Considerando:

𝐵𝑛 (−𝐷𝑛 ) = 𝐸𝑛

Assim:

𝐵𝑛 𝑠𝑒𝑛(𝜆𝑛 𝑥 ){[−𝐷𝑛 𝑡𝑔ℎ(𝜆𝑛 )]} = 𝑥(1 − 𝑥 )

𝐸𝑛 𝑡𝑔ℎ(𝜆𝑛 ) 𝑠𝑒𝑛(𝜆𝑛 𝑥) = 𝑥(1 − 𝑥 )

Sendo 𝜆𝑛 = 𝑛𝜋 e a função ortogonal 𝑚𝜋:

1 1
∫ 𝑥(1 − 𝑥) 𝑠𝑒𝑛(𝑚𝜋𝑥 ) 𝑑𝑥 = ∫ 𝐸𝑛 𝑠𝑒𝑛(𝑛𝜋𝑥 ) 𝑠𝑒𝑛(𝑚𝜋𝑥) 𝑡𝑔ℎ(𝑛𝜋)𝑑𝑥
0 0
1/2
1 1
∫ 𝑥(1 − 𝑥) 𝑠𝑒𝑛(𝑚𝜋𝑥 ) 𝑑𝑥 = 𝐸𝑛 𝑡𝑔ℎ(𝑛𝜋) ∫ 𝑠𝑒𝑛(𝑛𝜋𝑥 ) 𝑠𝑒𝑛(𝑚𝜋𝑥)𝑑𝑥
0 0

𝑛≠𝑚→0
Em que: {𝑛 = 𝑚 → 1
2

Assim:

1
𝐸𝑛 𝑡𝑔ℎ(𝑛𝜋)
∫ 𝑥(1 − 𝑥) 𝑠𝑒𝑛(𝑚𝜋𝑥 ) 𝑑𝑥 =
0 2

1
2 ∫ 𝑥(1 − 𝑥) 𝑠𝑒𝑛(𝑚𝜋𝑥 ) 𝑑𝑥 = 𝐸𝑛 𝑡𝑔ℎ(𝑛𝜋)
0

Resolvendo a integral por partes, temos que se 𝑛 for par, 𝐸𝑛 𝑡𝑔ℎ(𝑛𝜋) = 0, dessa
forma 𝑇(𝑥, 𝑦) = 0. Do contrário, para 𝑛 ímpar:

𝐵𝑛 𝑠𝑒𝑛(𝜆𝑛 𝑥){[−𝐷𝑛 𝑡𝑔ℎ(𝜆𝑛 ) cosh(𝜆𝑛 𝑦)] + 𝐷𝑛 𝑠𝑒𝑛ℎ(𝜆𝑛 𝑦)} = 𝑥(1 − 𝑥 )

𝑠𝑒𝑛ℎ (𝑛𝜋𝑦)
𝑇(𝑥, 𝑦) = 𝐸𝑛 𝑡𝑔ℎ(𝑛𝜋) 𝑠𝑒𝑛(𝑛𝜋𝑥 ) {cosh (𝑛𝜋𝑦) − }
𝑡𝑔ℎ(𝑛𝜋)
∞
8 𝑠𝑒𝑛 (𝑛𝜋𝑥) 𝑠𝑒𝑛ℎ (𝑛𝜋𝑦)
𝑇(𝑥, 𝑦) = ∑ {cosh (𝑛𝜋𝑦) − }
(𝑛𝜋)3 𝑡𝑔ℎ(𝑛𝜋)
𝑛

5
Para encontrar os fluxos, foi calculada a derivada em relação a x e y:
∞
8 (𝑛𝜋) 𝑐𝑜𝑠 (𝑛𝜋𝑥) 𝑠𝑒𝑛ℎ (𝑛𝜋𝑦)
𝑞𝑥 (𝑥, 𝑦) = ∑ 3
{cosh (𝑛𝜋𝑦) − }
(𝑛𝜋) 𝑡𝑔ℎ(𝑛𝜋)
𝑛

∞
8 𝑐𝑜𝑠 (𝑛𝜋𝑥) 𝑠𝑒𝑛ℎ (𝑛𝜋𝑦)
𝑞𝑥 (𝑥, 𝑦) = ∑ 2
{cosh (𝑛𝜋𝑦) − }
(𝑛𝜋) 𝑡𝑔ℎ(𝑛𝜋)
𝑛

E para y:
∞
8 𝑠𝑒𝑛 (𝑛𝜋𝑥) (𝑛𝜋)𝑐𝑜𝑠ℎ (𝑛𝜋𝑦)
𝑞𝑦 (𝑥, 𝑦) = ∑ {− (𝑛𝜋)senh (𝑛𝜋𝑦) − }
(𝑛𝜋)3 𝑡𝑔ℎ(𝑛𝜋)
𝑛

∞
8 𝑠𝑒𝑛 (𝑛𝜋𝑥) 𝑐𝑜𝑠ℎ (𝑛𝜋𝑦)
𝑞𝑦 (𝑥, 𝑦) = ∑ {− senh (𝑛𝜋𝑦) − }
(𝑛𝜋)2 𝑡𝑔ℎ(𝑛𝜋)
𝑛

O objetivo do trabalho é encontrar o erro entre a solução analítica e a solução

numérica no cálculo do fluxo. Assim, para a solução analítica, foram calculados
no software Matlab as temperaturas [𝑇(𝑥, 𝑦)] e os fluxos [𝑞 (𝑥, 𝑦)] em x e y para
todos os casos, de 16, 40 e 80 elementos (considerando elementos constantes).
Sendo que foi considerado 𝑛 = 101, 𝑛 ímpar.

Para a solução numérica, foram criados arquivos txt para os 3 casos (em anexo),
e utilizando o Fortran foi possível obter os dados de potencial (temperatura) e
derivada do potencial (fluxo), para assim, comparar com os resultados da
solução analítica. Neste caso em estudo a temperatura é prescrita, quanto a
análise do fluxo, que é o que queremos determinar, foi calculado o erro entre a
solução analítica e a numérica através da seguinte equação:

𝑉𝑛𝑢𝑚é𝑟𝑖𝑐𝑜 𝑉𝑎𝑛𝑎𝑙í𝑡𝑖𝑐𝑜
𝐸𝑟𝑟𝑜 = | | 𝑥 100
𝑀𝑎𝑖𝑜𝑟 𝑉𝑎𝑛𝑎𝑙í𝑡𝑖𝑐𝑜

6
 3. RESULTADOS E CONCLUSÃO

Foram criadas tabelas de 16, 40 e 80 elementos com todos os dados obtidos, do

analítico, numérico e o erro da derivada do potencial (fluxo) calculado.

Tabela 1: 16 elementos

Derivada Derivada Erro %

Potencial Potencial
x y Potencial Potencial Derivada
(Numérico) (Analítico)
(Numérico) (Analítico) Potencial
0,1250 0,0000 0,1094 0,10938 0,4037 0,4188 2,108868
0,3750 0,0000 0,2344 0,23437 0,7309 0,7136 2,4311669
0,6250 0,0000 0,2344 0,23437 0,7309 0,7136 2,4311388
0,8750 0,0000 0,1094 0,10938 0,4037 0,4188 2,10884
1,0000 0,1250 0,0000 0,00000 -0,7076 -0,5799 17,901215
1,0000 0,3750 0,0000 0,00000 -0,2257 -0,2478 3,0917053
1,0000 0,6250 0,0000 0,00000 -0,1000 -0,1034 0,4762006
1,0000 0,8750 0,0000 0,00000 -0,0240 -0,0283 0,6005552
0,8750 1,0000 0,0000 0,00000 -0,0238 -0,0269 0,4276118
0,6250 1,0000 0,0000 0,00000 -0,0654 -0,0648 0,0822388
0,3750 1,0000 0,0000 0,00000 -0,0654 -0,0648 0,0822247
0,1250 1,0000 0,0000 0,00000 -0,0238 -0,0269 0,4276076
0,0000 0,8750 0,0000 0,00000 -0,0240 -0,0283 0,6005454
0,0000 0,6250 0,0000 0,00000 -0,1000 -0,1034 0,4762146
0,0000 0,3750 0,0000 0,00000 -0,2257 -0,2478 3,0917053
0,0000 0,1250 0,0000 0,00000 -0,7076 -0,5799 17,901187

Tabela 2: 40 elementos

Erro %
Potencial Potencial Derivada Potencial Derivada Potencial
X y Derivada
(Numérico) (Analítico) (Numérico) (Analítico)
Potencial
0,0500 0,0000 0,0475 0,04750 0,2394 0,2263 1,6987487
0,1500 0,0000 0,1275 0,12750 0,4763 0,4673 1,16573182
0,2500 0,0000 0,1875 0,18750 0,6165 0,6124 0,52359225
0,3500 0,0000 0,2275 0,22750 0,7021 0,6993 0,36764073
0,4500 0,0000 0,2475 0,24750 0,7428 0,7404 0,31097473
0,5500 0,0000 0,2475 0,24750 0,7428 0,7404 0,31094889
0,6500 0,0000 0,2275 0,22750 0,7021 0,6993 0,36769241
0,7500 0,0000 0,1875 0,18750 0,6165 0,6124 0,52357933
0,8500 0,0000 0,1275 0,12750 0,4763 0,4673 1,16579643

7
0,9500 0,0000 0,0475 0,04750 0,2394 0,2263 1,69877454
1,0000 0,0500 0,0000 0,0000 -0,8869 -0,7739 14,6080446
1,0000 0,1500 0,0000 0,0000 -0,5135 -0,5303 2,16938366
1,0000 0,2500 0,0000 0,0000 -0,3730 -0,3762 0,41400801
1,0000 0,3500 0,0000 0,0000 -0,2671 -0,2694 0,28976687
1,0000 0,4500 0,0000 0,0000 -0,1912 -0,1926 0,1791194
1,0000 0,5500 0,0000 0,0000 -0,1354 -0,1363 0,11154435
1,0000 0,6500 0,0000 0,0000 -0,0934 -0,0939 0,06364124
1,0000 0,7500 0,0000 0,0000 -0,0608 -0,0610 0,02767533
1,0000 0,8500 0,0000 0,0000 -0,0344 -0,0343 0,01372972
1,0000 0,9500 0,0000 0,0000 -0,0100 -0,0111 0,14269509
0,9500 1,0000 0,0000 0,0000 -0,0100 -0,0110 0,1234457
0,8500 1,0000 0,0000 0,0000 -0,0323 -0,0319 0,06092123
0,7500 1,0000 0,0000 0,0000 -0,0500 -0,0496 0,05095245
0,6500 1,0000 0,0000 0,0000 -0,0629 -0,0625 0,05255492
0,5500 1,0000 0,0000 0,0000 -0,0697 -0,0693 0,05309547
0,4500 1,0000 0,0000 0,0000 -0,0697 -0,0693 0,05312002
0,3500 1,0000 0,0000 0,0000 -0,0629 -0,0625 0,05256268
0,2500 1,0000 0,0000 0,0000 -0,0500 -0,0496 0,05095762
0,1500 1,0000 0,0000 0,0000 -0,0323 -0,0319 0,06092123
0,0500 1,0000 0,0000 0,0000 -0,0100 -0,0110 0,12345216
0,0000 0,9500 0,0000 0,0000 -0,0100 -0,0111 0,14267894
0,0000 0,8500 0,0000 0,0000 -0,0344 -0,0343 0,01372584
0,0000 0,7500 0,0000 0,0000 -0,0608 -0,0610 0,02768179
0,0000 0,6500 0,0000 0,0000 -0,0934 -0,0939 0,06365675
0,0000 0,5500 0,0000 0,0000 -0,1354 -0,1363 0,11150559
0,0000 0,4500 0,0000 0,0000 -0,1912 -0,1926 0,17913233
0,0000 0,3500 0,0000 0,0000 -0,2671 -0,2694 0,28975395
0,0000 0,2500 0,0000 0,0000 -0,3730 -0,3762 0,41400801
0,0000 0,1500 0,0000 0,0000 -0,5135 -0,5303 2,16937074
0,0000 0,0500 0,0000 0,0000 -0,8869 -0,7739 14,6079412

Tabela 3: 80 elementos

Erro %
Potencial Potencial Derivada Potencial Derivada Potencial
X y Derivada
(Numérico) (Analítico) (Numérica) (Analítico)
Potencial
0,0250 0,0000 0,0244 0,02438 0,2167 0,1351 9,44090
0,0750 0,0000 0,0694 0,06938 0,3062 0,3002 0,69325
0,1250 0,0000 0,1094 0,10938 0,4237 0,4188 0,57078
0,1750 0,0000 0,1444 0,14438 0,5131 0,5102 0,33676

8
0,2250 0,0000 0,1744 0,17438 0,5843 0,5824 0,21573
0,2750 0,0000 0,1994 0,19937 0,6403 0,6388 0,16877
0,3250 0,0000 0,2194 0,21937 0,6833 0,6822 0,12303
0,3750 0,0000 0,2344 0,23437 0,7145 0,7136 0,10922
0,4250 0,0000 0,2444 0,24437 0,7349 0,7341 0,08999
0,4750 0,0000 0,2494 0,24937 0,7450 0,7442 0,08497
0,5250 0,0000 0,2494 0,24937 0,7449 0,7442 0,08177
0,5750 0,0000 0,2444 0,24437 0,7348 0,7341 0,07982
0,6250 0,0000 0,2344 0,23437 0,7143 0,7136 0,09066
0,6750 0,0000 0,2194 0,21937 0,6830 0,6822 0,09383
0,7250 0,0000 0,1994 0,19937 0,6399 0,6388 0,12504
0,7750 0,0000 0,1744 0,17438 0,5837 0,5824 0,14846
0,8250 0,0000 0,1444 0,14438 0,5122 0,5102 0,23239
0,8750 0,0000 0,1094 0,10938 0,4220 0,4188 0,37009
0,9250 0,0000 0,0694 0,06938 0,3074 0,3002 0,83274
0,9750 0,0000 0,0244 0,02438 0,1623 0,1351 3,15290
1,0000 0,0250 0,0000 0,00000 -0,9688 -0,8649 12,02088
1,0000 0,0750 0,0000 0,00000 -0,6849 -0,6994 1,67996
1,0000 0,1250 0,0000 0,00000 -0,5775 -0,5799 0,27819
1,0000 0,1750 0,0000 0,00000 -0,4840 -0,4858 0,20715
1,0000 0,2250 0,0000 0,00000 -0,4081 -0,4093 0,13822
1,0000 0,2750 0,0000 0,00000 -0,3451 -0,3459 0,09951
1,0000 0,3250 0,0000 0,00000 -0,2921 -0,2927 0,07337
1,0000 0,3750 0,0000 0,00000 -0,2473 -0,2478 0,05486
1,0000 0,4250 0,0000 0,00000 -0,2092 -0,2096 0,04066
1,0000 0,4750 0,0000 0,00000 -0,1767 -0,1769 0,02959
1,0000 0,5250 0,0000 0,00000 -0,1486 -0,1488 0,02067
1,0000 0,5750 0,0000 0,00000 -0,1244 -0,1245 0,01303
1,0000 0,6250 0,0000 0,00000 -0,1034 -0,1034 0,00653
1,0000 0,6750 0,0000 0,00000 -0,0849 -0,0849 0,00086
1,0000 0,7250 0,0000 0,00000 -0,0686 -0,0686 0,00448
1,0000 0,7750 0,0000 0,00000 -0,0540 -0,0539 0,00952
1,0000 0,8250 0,0000 0,00000 -0,0407 -0,0406 0,01453
1,0000 0,8750 0,0000 0,00000 -0,0285 -0,0283 0,01967
1,0000 0,9250 0,0000 0,00000 -0,0170 -0,0167 0,03063
1,0000 0,9750 0,0000 0,00000 -0,0052 -0,0055 0,03690
0,9750 1,0000 0,0000 0,00000 -0,0052 -0,0055 0,03297
0,9250 1,0000 0,0000 0,00000 -0,0167 -0,0164 0,03862
0,8750 1,0000 0,0000 0,00000 -0,0271 -0,0269 0,03130
0,8250 1,0000 0,0000 0,00000 -0,0369 -0,0367 0,03004
0,7750 1,0000 0,0000 0,00000 -0,0458 -0,0456 0,02905

9
0,7250 1,0000 0,0000 0,00000 -0,0536 -0,0534 0,02831
0,6750 1,0000 0,0000 0,00000 -0,0601 -0,0598 0,02790
0,6250 1,0000 0,0000 0,00000 -0,0651 -0,0648 0,02765
0,5750 1,0000 0,0000 0,00000 -0,0685 -0,0682 0,02740
0,5250 1,0000 0,0000 0,00000 -0,0702 -0,0700 0,02736
0,4750 1,0000 0,0000 0,00000 -0,0702 -0,0700 0,02736
0,4250 1,0000 0,0000 0,00000 -0,0685 -0,0682 0,02739
0,3750 1,0000 0,0000 0,00000 -0,0651 -0,0648 0,02770
0,3250 1,0000 0,0000 0,00000 -0,0601 -0,0598 0,02821
0,2750 1,0000 0,0000 0,00000 -0,0536 -0,0534 0,02939
0,2250 1,0000 0,0000 0,00000 -0,0459 -0,0456 0,03220
0,1750 1,0000 0,0000 0,00000 -0,0370 -0,0367 0,03377
0,1250 1,0000 0,0000 0,00000 -0,0283 -0,0269 0,16615
0,0750 1,0000 0,0000 0,00000 -0,0119 -0,0164 0,51405
0,0250 1,0000 0,0000 0,00000 0,0175 -0,0055 2,66384
0,0000 0,9750 0,0000 0,00000 -0,0046 -0,0055 0,10826
0,0000 0,9250 0,0000 0,00000 -0,0169 -0,0167 0,02659
0,0000 0,8750 0,0000 0,00000 -0,0284 -0,0283 0,01537
0,0000 0,8250 0,0000 0,00000 -0,0407 -0,0406 0,00834
0,0000 0,7750 0,0000 0,00000 -0,0539 -0,0539 0,00163
0,0000 0,7250 0,0000 0,00000 -0,0685 -0,0686 0,00494
0,0000 0,6750 0,0000 0,00000 -0,0848 -0,0849 0,01169
0,0000 0,6250 0,0000 0,00000 -0,1033 -0,1034 0,01929
0,0000 0,5750 0,0000 0,00000 -0,1243 -0,1245 0,02816
0,0000 0,5250 0,0000 0,00000 -0,1485 -0,1488 0,03844
0,0000 0,4750 0,0000 0,00000 -0,1765 -0,1769 0,05165
0,0000 0,4250 0,0000 0,00000 -0,2090 -0,2096 0,06778
0,0000 0,3750 0,0000 0,00000 -0,2470 -0,2478 0,09009
0,0000 0,3250 0,0000 0,00000 -0,2917 -0,2927 0,12096
0,0000 0,2750 0,0000 0,00000 -0,3445 -0,3459 0,16826
0,0000 0,2250 0,0000 0,00000 -0,4071 -0,4093 0,25569
0,0000 0,1750 0,0000 0,00000 -0,4828 -0,4858 0,35295
0,0000 0,1250 0,0000 0,00000 -0,5590 -0,5799 2,41581
0,0000 0,0750 0,0000 0,00000 -0,7625 -0,6994 7,28931
0,0000 0,0250 0,0000 0,00000 -0,9688 -0,8649 12,02088

A partir da análise dos resultados encontrados nas soluções numérica e analítica

da temperatura (potencial), podemos observar que os valores são iguais. E como
a temperatura é prescrita no problema em análise, esse era o resultado esperado
quanto ao potencial.

10
Além disso, foi calculada a média dos erros para analisar a derivada do potencial
(fluxo):

Tabela 4: Média dos erros

16 el. 40 el. 80 el.

3,38993901 1,12136317 0,73459

Dessa forma, foi possível gerar o seguinte gráfico com a média dos erros (%) no
cálculo do fluxo para 16, 40 e 80 elementos:

Gráfico 1: Erro médio

3,5

2,5
Erro médio

1,5

0,5

0
0 10 20 30 40 50 60 70 80 90
Elementos

A partir do gráfico acima e das tabelas (1), (2) e (3) é possível concluir que os
erros estão relativamente baixos e que quanto maior o número de elementos,
menor é a média dos erros, ou seja, mais preciso é o resultado. Assim, para
obter resultados mais precisos devemos refinar mais, uma vez que o MEC
apresenta a tendência de maior precisão ao aumentar o número de elementos,
como se é esperado de um método de discretização.

11