Solução Numérica da Eq. de Quantidade de Movimento Linear 1D em Regime Permanente com Propriedades Constantes

Dados de entrada

                10  = N:  número de volumes de controle
  5.0000000000E-02  = R:  raio do duto [m]
 -1.6000000000E+01  = C:  constante (queda de pressão) [Pa/m]
  1.0000000000E-03  = mi: viscosidade dinâmica [Pa.s]
  2.0000000000E-01  = L:  comprimento do duto [m]


Coeficientes e termos-fontes

    Volume    aw(i)                    ae(i)                    ap(i)                    bp(i)
         1    0.000000000000000E+00    5.000000000000000E-03    5.000000000000000E-03    1.000000000000000E-03
         2    5.000000000000000E-03    1.000000000000000E-02    1.500000000000000E-02    3.000000000000000E-03
         3    1.000000000000000E-02    1.500000000000000E-02    2.500000000000000E-02    5.000000000000000E-03
         4    1.500000000000000E-02    2.000000000000000E-02    3.500000000000000E-02    7.000000000000001E-03
         5    2.000000000000000E-02    2.500000000000000E-02    4.500000000000000E-02    8.999999999999999E-03
         6    2.500000000000000E-02    3.000000000000000E-02    5.500000000000000E-02    1.100000000000000E-02
         7    3.000000000000000E-02    3.500000000000000E-02    6.500000000000000E-02    1.300000000000000E-02
         8    3.500000000000000E-02    4.000000000000000E-02    7.500000000000001E-02    1.500000000000000E-02
         9    4.000000000000000E-02    4.500000000000000E-02    8.499999999999999E-02    1.700000000000000E-02
        10    4.500000000000000E-02    0.000000000000000E+00    1.450000000000000E-01    1.900000000000000E-02


Soluções numéricas

    Volume     Posição [m]        u (numérico) [m/s]       u (analítico) [m/s]      Erro numérico [m/s]
         0     0.000000000E+00    1.002500000000000E+01    1.000000000000000E+01   -2.499999999999325E-02
         1     2.500000000E-03    9.999999999999995E+00    9.975000000000001E+00   -2.499999999999325E-02
         2     7.500000000E-03    9.799999999999995E+00    9.775000000000000E+00   -2.499999999999503E-02
         3     1.250000000E-02    9.399999999999993E+00    9.375000000000002E+00   -2.499999999999147E-02
         4     1.750000000E-02    8.799999999999994E+00    8.775000000000000E+00   -2.499999999999325E-02
         5     2.250000000E-02    7.999999999999996E+00    7.975000000000001E+00   -2.499999999999503E-02
         6     2.750000000E-02    6.999999999999996E+00    6.975000000000002E+00   -2.499999999999414E-02
         7     3.250000000E-02    5.799999999999997E+00    5.775000000000000E+00   -2.499999999999680E-02
         8     3.750000000E-02    4.399999999999998E+00    4.375000000000002E+00   -2.499999999999591E-02
         9     4.250000000E-02    2.799999999999999E+00    2.775000000000000E+00   -2.499999999999858E-02
        10     4.750000000E-02    9.999999999999996E-01    9.750000000000009E-01   -2.499999999999869E-02
        11     5.000000000E-02    0.000000000000000E+00    0.000000000000000E+00    0.000000000000000E+00


Velocidade média (retângulo) - res. numérico  [m/s]:        5.049999999999997E+00
Velocidade média - res. analítico [m/s]:                    5.000000000000001E+00
Erro numérico para a velocidade média (retângulo):         -4.999999999999627E-02


Velocidade média (interpolação) - res. numérico  [m/s]:     5.012499999999998E+00
Velocidade média - res. analítico [m/s]:                    5.000000000000001E+00
Erro numérico para a velocidade média (interpolação):      -1.249999999999663E-02


Velocidade máxima - res. numérico [m/s]:                    1.002500000000000E+01
Velocidade máxima - res. analítico [m/s]:                   1.000000000000000E+01


Força do fluido sobre o duto - res. numérico [N]:           2.513274122871834E-02
Força do fluido sobre o duto - res. analítico [N]:          2.513274122871835E-02