from pyXSteam.XSteam import XSteam
steamTable = XSteam(XSteam.UNIT_SYSTEM_FLS) # ft/lb/sec/°F/psi/btu
def flash(mc, Phc, Pf):
hc = steamTable.h_pt(Phc+14.696, steamTable.tsat_p(Phc + 14.696) - 1)
hfs = steamTable.h_pt(Pf+14.696, steamTable.tsat_p(Pf + 14.696) + 1)
hfc = steamTable.h_pt(Pf+14.696, steamTable.tsat_p(Pf + 14.696) - 1)
return mc * (hc - hfc) / (hfs - hfc)
print("flash steam flow = ", flash(100, 120, 60)) # 100 lb/hr, 120 psig, 60 psig