Solução Numérica do Escoamento Unidimensional de Fluido Incompressível: Equações da Conservação da Massa e da Quantidade de Movimento Linear Problema de Moody Dados de entrada 12 = N: número de volumes de controle (incluindo fictícios) 5.0000000000E+00 = L: comprimento do domínio [m] 1.0000000000E-03 = mi: viscosidade absoluta [Pa.s] 1.0000000000E+03 = ro: massa específica [kg/m3] 2.0000000000E-02 = f: fator de atrito de Darcy [adim.] 1.0000000000E+01 = Uin: velocidade na entrada do duto [m/s] 2.0000000000E-02 = D0: diâmetro de entrada do duto [m] 4.0000000000E-03 = cd: fator de expansão do diâmetro do duto [adim.] 1.0000000000E+00 = dt: passo de tempo [s] 0.0000000000E+00 = beta : coeficiente do termo difusivo (funções de interpolação): 0 = UDS; 1 = CDS 10000 = itmax: número máximo de iterações (geral) 2 = itmaxm: número máximo de iterações (ciclo da massa) Posições, diâmetros e áreas das seções transversais (pontos nodais) Volume i Posição [m] Diâmetro [m] Área da seção transversal [m2] 1 0.000000000000000E+00 2.000000000000000E-02 3.141592653589793E-04 2 2.500000000000000E-01 2.100000000000000E-02 3.463605900582747E-04 3 7.500000000000000E-01 2.300000000000000E-02 4.154756284372501E-04 4 1.250000000000000E+00 2.500000000000000E-02 4.908738521234052E-04 5 1.750000000000000E+00 2.700000000000000E-02 5.725552611167398E-04 6 2.250000000000000E+00 2.900000000000000E-02 6.605198554172541E-04 7 2.750000000000000E+00 3.100000000000000E-02 7.547676350249478E-04 8 3.250000000000000E+00 3.300000000000000E-02 8.552985999398212E-04 9 3.750000000000000E+00 3.500000000000000E-02 9.621127501618743E-04 10 4.250000000000000E+00 3.700000000000001E-02 1.075210085691107E-03 11 4.750000000000000E+00 3.900000000000000E-02 1.194590606527519E-03 12 5.000000000000000E+00 4.000000000000000E-02 1.256637061435917E-03 Coeficientes e termos-fontes (QML) Volume i aW(i) aP(i) aE(i) bP(i) 1 0.000000000000000E+00 1.000000000000000E+00 -1.000000000000000E+00 2.000000000000000E+01 2 3.141593281908323E+00 4.056794018197393E+00 7.602654221687299E-07 1.937764246225734E+00 3 3.141593413855215E+00 4.027581274281614E+00 9.047786842338605E-07 1.977031617607192E+00 4 3.141593558368477E+00 4.012055062831132E+00 1.061858316913350E-06 1.951209764974955E+00 5 3.141593715448110E+00 4.007243332374027E+00 1.231504320207199E-06 1.895933168568996E+00 6 3.141593885094113E+00 4.011589168250133E+00 1.413716694115407E-06 1.845905776846603E+00 7 3.141594067306486E+00 4.024092654975725E+00 1.608495438637974E-06 1.807645663923755E+00 8 3.141594262085231E+00 4.043916938867884E+00 1.815840553774901E-06 1.777078837069465E+00 9 3.141594469430347E+00 4.070348658355687E+00 2.035752039526187E-06 1.751062605936435E+00 10 3.141594689341832E+00 4.102830364696986E+00 2.268229895891832E-06 1.729028748959849E+00 11 3.141594921819689E+00 4.141546522471601E+00 2.513274122871835E-06 1.728999949621754E+00 12 1.000000000000000E+00 1.000000000000000E+00 0.000000000000000E+00 -2.864685201230142E-01 Coeficientes e termos-fontes (MASSA) Volume i aW(i) aP(i) aE(i) bP(i) 1 0.000000000000000E+00 1.000000000000000E+00 0.000000000000000E+00 0.000000000000000E+00 2 1.188946584054528E-07 2.799558377516147E-07 1.610611793461618E-07 4.336808689942018E-19 3 1.610611793461618E-07 3.946898257407335E-07 2.336286463945717E-07 -4.336808689942018E-19 4 2.336286463945717E-07 5.589135917739317E-07 3.252849453793600E-07 0.000000000000000E+00 5 3.252849453793600E-07 7.626664466633988E-07 4.373815012840389E-07 0.000000000000000E+00 6 4.373815012840389E-07 1.007988405764050E-06 5.706069044800115E-07 -4.336808689942018E-19 7 5.706069044800115E-07 1.295696997084491E-06 7.250900926044791E-07 4.336808689942018E-19 8 7.250900926044791E-07 1.625659499025813E-06 9.005694064213338E-07 0.000000000000000E+00 9 9.005694064213338E-07 1.997072872512135E-06 1.096503466090802E-06 0.000000000000000E+00 10 1.096503466090802E-06 2.408234874085009E-06 1.311731407994207E-06 0.000000000000000E+00 11 1.311731407994207E-06 2.812974669630924E-06 1.501243261636717E-06 0.000000000000000E+00 12 0.000000000000000E+00 1.000000000000000E+00 0.000000000000000E+00 0.000000000000000E+00 Soluções numéricas - Velocidades nodais Volume Posição [m] Velocidade [m/s] 1 0.000000000E+00 1.000000000000000E+01 2 2.500000000E-01 8.997798102990291E+00 3 7.500000000E-01 7.509335896088595E+00 4 1.250000000E+00 6.366437261212128E+00 5 1.750000000E+00 5.464279675221650E+00 6 2.250000000E+00 4.739383436188342E+00 7 2.750000000E+00 4.149226141892936E+00 8 3.250000000E+00 3.662852134490656E+00 9 3.750000000E+00 3.257279821018151E+00 10 4.250000000E+00 2.915569652283207E+00 11 4.750000000E+00 2.629101132160193E+00 12 5.000000000E+00 2.485866872098685E+00 Soluções numéricas - Pressões Volume Posição [m] Pressão [Pa] 1 0.000000000E+00 0.000000000000000E+00 2 2.500000000E-01 -5.478725453051807E+02 3 7.500000000E-01 -1.643617635915542E+03 4 1.250000000E+00 -2.555492553674823E+03 5 1.750000000E+00 -3.227124159789862E+03 6 2.250000000E+00 -3.713920928925282E+03 7 2.750000000E+00 -4.076992164591235E+03 8 3.250000000E+00 -4.354634312894516E+03 9 3.750000000E+00 -4.569598339501530E+03 10 4.250000000E+00 -4.737390663607423E+03 11 4.750000000E+00 -4.870197828430370E+03 12 5.000000000E+00 -4.936601410841839E+03 Soluções numéricas - Velocidades nas faces Volume Posição [m] Velocidade [m/s] 1 0.000000000E+00 1.000000000000000E+01 2 5.000000000E-01 8.264462809917356E+00 3 1.000000000E+00 6.944444444444443E+00 4 1.500000000E+00 5.917159763313608E+00 5 2.000000000E+00 5.102040816326531E+00 6 2.500000000E+00 4.444444444444444E+00 7 3.000000000E+00 3.906250000000000E+00 8 3.500000000E+00 3.460207612456747E+00 9 4.000000000E+00 3.086419753086418E+00 10 4.500000000E+00 2.770083102493074E+00 11 5.000000000E+00 2.500000000000000E+00