Isothermal flash calculation using Soave equation of state

Multicomponent flash calculations using an equation of state and classical mixing rules.

Description

You only have to run main01.m, is easy to change the composition and operating conditions. If you want to add new components just add the critical propeties and acentric factor. Also is possible to use binary interaction parameters differents to zero.

• MATLAB release: MATLAB 7.2 (R2006a)

[Jeremy Kozdon]

This seems to be correct. I would suggest though, since it is a driver file, to remove the global statements.

[Jeremy Kozdon]

Ran into some NAN problems with the following inputs.

Tc = [304.2 190.6 630.1];
Pc = [72.9 45.4 30.17];
w = [0.228 0.008 0.4489];
T0 = 450;
P0 = 300;
%feed composition
xy01 = [0 0 1];