##############################
# chem.py                    #
# Didier Descamps 08/03/2019 #
##############################

import sys
from math import pi

if(len(sys.argv)!=3):
    sys.exit("il faut entrer dans l'ordre la puissance (kW)\
 et la température (°C) maximales :") 
P=float(sys.argv[1])    
t=float(sys.argv[2])
h=1 

## débit volume à tirer 
c=1.0	    # propane : 1 m3(n)/kWh(pci) ; 1,07 pour le gaz naturel
qf=c*P*(t+273.0)/3600.0/273.0
print("débit des fumées=""{0:3.5f}"" m3/s".format(qf)) 

## pression motrice pour 1 m de hauteur 
rho0=1.2695 # 1.2345 pour le gaz nat    
pmot=rho0*(1.0-273.0/(t+273.0))*9.81*h 

## calcul direct de a 
lamda=0.04 # facteur de pertes regime laminaire, Re~5500  
rhof=rho0*273.0/(t+273.0)  
a5=lamda*h*rhof*qf*qf/2.0/pmot 
a=a5**0.2 
print("pour le propane. (si gaz naturel, ajouter 3%)")
print("section carrée : a=" "{0:3.1f}"" mm".format(a*1000.0))
d=a*((16.0/pi**2)**0.2)  
print("section circulaire : dia=" "{0:3.1f}"" mm".format(d*1000.0))
print("vitesse=""{0:3.1f}"" m/s".format(qf/a/a))