The Boston-Sullivan Algorithm

This algorithm is based on the use of newly defined energy and volatility parameters as the primary successive approximation variables, and Broyden s method is used to iterate on these parameters. A brief review of the Boston-Sullivan Algorithm follows. 

The component-material balances are formulated in a manner analogous to that shown in Chap. 2 except for the fact that Boston and Sullivan3 stated these balances in terms of the liquid rates /,-, rather than the vapor rates t ,. Temperatures were computed use of a variation of the method wherein a different base component is used for each plate as suggested by Billingsley.2 Partial molar enthalpies were also expressed in terms of the deviation function Q. 

A formulation of the model equations follows. The S s appearing in the component-material balance for each stage j 

The K value for the base component was taken to be an exponential function of temperature as follows  

The vapor enthalpies were expressed in terms of the enthalpy departure function Q (see Chap. 14) as follows 

---