Included in the zipped file you find the prsrk.m routine to calculate molar volume, compressibility factor, fugacity coefficient, Helmholtz energy, entropy, enthalpy, internal energy, and Gibbs energy of vapor and liquid mixtures. The other Matlab functions/scripts are part of a steady-state distillation example as depicted in Fig. 15.6 (pp. 568) of Henley, E. J. and Seader, J. D., "Equilibrium-Stage Separation Operations in Chemical Engineering", 1981, John Wiley & Sons. The files include:
1 ) main_dist.m is the main script containing input data and the call to fmincon.
2 ) feed_dist.m performes BP, DP and flash calculations for the feed conditions.
3 ) model_dist.m is the vectorized distillation column model.
4 ) nonlcon_dist.m is the nonlinear function argument of fmincon.
5 ) parameters_dist.m includes all the parameters needed for the simulation such as tehrmodynamic parameters and pressure drop along the stages.
The results of the simulation in Matlab agree perfectly with those found by simulating the same system in Aspen Plus 2006.