Calculo do experimento 2 do laboratorio de PME3033

exp2.jl

######################################################
# Por: Guilherme S Salustiano                        #
# Eng da Computação 2020                             #
######################################################

using Measurements

temperatura_da_agua_graus = 30 ± 0.5
distancia_tubos_pilot_centros = 5.4e-3
distancia_tubos_pilot_parede_centro = distancia_tubos_pilot_centros/2
inclinação_tubos_pilot = pi/6
comprimento_pressao_estatica = 20e-3
altura_pressao_estatica = comprimento_pressao_estatica * sin(inclinação_tubos_pilot)
massa_agua = 219.1 ± 0.1
tempo_agua = 323.84 ± 0.01
densidade_agua = 995.61
viscosidade_cinematica=817e-9
raio_tubo = 0.027
area_tubo = pi*raio_tubo^2
gravidade = 9.8

comprimento_colunas_tubos_pitot = [44e-3, 53e-3, 54e-3, 56e-3, 52e-3, 55e-3, 50e-3, 48e-3, 47e-3, 49e-3]
comprimento_colunas_tubos_pitot = map(x -> x ± 0.1, comprimento_colunas_tubos_pitot) # adicionando a incerteza de 0.1
altura_colunas_tubos_pitot = map(x -> x * sin(inclinação_tubos_pilot), comprimento_colunas_tubos_pitot)
velocidades_tubos_pitot = map(x -> sqrt(2*gravidade*(x-altura_pressao_estatica)), altura_colunas_tubos_pitot)
velocidade_medias = [(velocidades_tubos_pitot[5-x]+velocidades_tubos_pitot[6+x])/2 for x = 0:4 ]

distancia_tubos_centro = [distancia_tubos_pilot_centros/2 + x*distancia_tubos_pilot_centros for x in -5:5]

vazao_massica = massa_agua/tempo_agua
vazao_volumetrica = vazao_massica/densidade_agua
velocidade_media = vazao_volumetrica/area_tubo

numero_reynolds = velocidade_media*raio_tubo*2/viscosidade_cinematica

println("vazão massica: ", vazao_massica)
println("vazão volumetrica: ", vazao_volumetrica)
println("velocidade média: ", velocidade_media)
println("numero de Reynolds: ", numero_reynolds)
println("velocidades: ", )
for v in velocidades_tubos_pitot
  println("\t", v.val, "\t± ", v.err)
end

println("velocidades médias: ", )
for v in velocidade_medias
  println("\t", v.val, "\t± ", v.err)
end

println("(velocidades médias)^2: ", )
for v in map(x -> x^2, velocidade_medias)
  println("\t", v.val, "\t± ", v.err)
end

println("(velocidades médias)^3: ", )
for v in map(x -> x^3, velocidade_medias)
  println("\t", v.val, "\t± ", v.err)
end