Rate-based column simulation

feed -> distillate + bottoms (tray efficiency)

Parameters

same params as distillation, plus stage_efficiency [0..1 Murphree, optional — else O'Connell from feed relative volatility] and liquid_viscosity_cp [cP, ~0.2]; N actual trays behave as N*E_o equilibrium stages (nonequilibrium)