Untitled

P_atm = 101,325
T_água = 23
 
pv_água = p_sat(Water;T=T_água) "pressão de saturação da água"
rho_água = density(Water;T=T_água;P=P_atm) "densidade da água"
mu = viscosity(Water;T=T_água;P=P_atm) "viscosidade dinâmica da água"
h_vap = (pv_água*1000)/(rho_água*g#) "altura manométrica devido à pressão da água"
H_atm = (P_atm*1000)/(rho_água*g#) "altura manométrica devido à pressão atmosférica"
 
 
{Parâmetros da linha para D_ideal de acordo com Janna}
Q_dot_volume = 3250 "m³/h"
m_dot = rho_água*Q_dot_volume/3600 "kg/s"
n = 1
epsilon = 0,000046
Ro = pi#*epsilon*mu/(4*m_dot)
C_2 = 0,07 "$/kWh"
t = 8640 "h/ano com 5 dias de inatividade"
a = 0,05 "taxa de amortização 1/T em anos"
b = 0,01 "fração do custo de manutenção anual - valor típico"
F = 7 "multiplicador do custo do tubo de 6 a 7"
C_1 = 1000 "constante na curva do ajuste de dados do custo do tubo $/m² valor típico"
eta = 0,82 "eficiência da bomba"
 
 
f_Re = ( 128/(5*pi#^3)*(m_dot^2/mu^5)*(( (4*m_dot)/(pi#*mu) )^n)*( n*(a + b)*(1 + F)*C_1*eta*rho_água^2 )/(C_2*t) )^(1/6)
 
f_atrito = 0,013 "gráfico correspondência"
D_ideal = ( ( 40*f_atrito*m_dot^3*C_2*t )/( n*(a + b)*(1 + F)*C_1*eta*pi#^2*rho_água^2 ) ) ^ ( 1 / (n + 5) )
D_ideal_pol = D_ideal*39,37