where N j vap/liq are the number of moles in the bulk liquid or vapor in the j th tray of the distillation column. We call the condenser Tray 1 and the reboiler as Tray 40. V j vap/liq,out are the vapor/liquid volumetric flowrate leaving the j th tray, C j vap/liq are concentrations in the bulk liquid or vapor in the j th tray, k liq A are the effective rates of distillation/condensation of each species, T j is the temperature in the j th tray, and k H are Henry's Law constants. S liq 13 is the matrix of stoichiometric coefficients for liquid+liquid reactions, r liq are the reaction rates, V liq is the volume of liquid in each tray. Here we assume all the trays are full, with a diameter of 2.5 m and a height of 0.3 m, i.e., V 2 liq = V 3 liq = ...V 39 liq = (2.5/2) 2 × π × 0.3 = 1.47m 3 . We assume the liquid volume in the condenser and the reboiler are the same, i.e.,