External-degrees-of-freedom parameter in the UNIFAC-FV model

Jul 1, 1982 - External-degrees-of-freedom parameter in the UNIFAC-FV model. Reply to comments. T. Oishi, J. M. Prausnitz. Ind. Eng. Chem. Process Des...
1 downloads 11 Views 127KB Size
Ind. Eng. Chem. Process Des. Dev.

treated as an adjustable parameter which depends on both and the particular desired property the (solvent activity, Henry’s constant, x , etc.). It is also important to realize that c1 values cannot be estimated from values of other solvents with similar chemical structure or size. Literature Cited Arai, Y.; Iwal, Y. Ind. Eng. Chem. Process D e s . Dev. 1980, 19, 508.

1982,21, 537

537

Beret, S.; Prausnitz, J. M. AIChE J. 1975, 21, 1123. Gottlieb, M.; Herskowitz, M. Macromolecules 1981, 14, 1468. Olshi, T.; Prausnitz, J. M. Ind. Eng. Chem. Process Des. Dev. 1978, 17, 333.

Chemical Engineering Department Moshe Gottlieb* Ben Gurion University Mordechai Herskowitz Beer Sheva 84120, Israel

Response to “External-Degrees-of-Freedom Parameter In the UNIFAC-FV Model” Sir: Gottlieb and Herskowitz have pointed out what is all-too-often forgotten: the UNIFAC model provides no more than a rough approximation which should only be used to estimate thermodynamic properties when, as unfortunately happens so often, no experimental data at all are available for a particular liquid mixture. For liquid mixtures, Raoult’s law provides a zeroth approximation. UNIFAC provides a first approximation. But even that may have appreciable error. What is remarkable about UNIFAC is not that it sometimes gives poor results. What is remarkable is that so often it does surprisingly well. Reliable vapor-liquid equilibrium data are rare for polymer solutions. Oishi observed that when UNIFAC was applied to such data, using parameters obtained from ordinary liquid-mixture data, calculated solvent activities were consistently low. Since UNIFAC is based on a lattice-type theory of liquid solutions (UNIQUAC),and since Flory and Patterson have shown that for polymer solutions a lattice-type theory must be corrected for appreciable free-volume effects, it appeared reasonable to add to the UNIFAC equation a free-volume contribution based on Flory’s equation of state. Upon doing so, Oishi found that he could significantly reduce the consistent deviation between observed solvent activities and those calculated by (uncorrected) UNIFAC. Upon applying Oishi’s method to experimental data for solutions of poly(dimethylsiloxane), Gottlieb and Herskowitz obtain some appreciable deviations; they propose to remove these deviations by assigning some unrealistic values to the external-degrees-of-freedom parameter cl. This is an unsatisfactory procedure. A possible explanation for the observed deviations may be related to the well-known fact that poly(dimethy1siloxanes) are much more flexible than the polymers considered by Oishi; as shown by Lichtenthaler and Liu (1974), the configurational contribution to the Gibbs energy of poly(dimethylsi1oxane) solutions appears to be somewhat different from that calculated by Staverman’s formula, which is used in UNIFAC. The inadequacy of Staverman’s formula may therefore be responsible, in part, for the results reported by Gottlieb and Herskowitz.

0196-4305/82/1121-0537$01.25/0

However, a more likely explanation follows from uncertainties in the residual UNIFAC group-energy parameters. Calculations are often sensitive to UNIFAC energy parameters. These parameters have been revised several times as new experimental data have become available for ordinary liquid mixtures. Since reliable data are rare for mixtures containing hydrocarbons and low-molecular dimethylsiloxanes, residual UNIFAC group-energy parameters are not known well for the group-group interactions required for UNIFAC calculations as performed by Gottlieb and Herskowitz. It is relevant to note that such parameters have not been reported by Fredenslund and co-workers (Lyngby) nor by Gmehling and co-workers (Dortmund), the two research teams that have concerned themselves most extensively with data reduction to obtain UNIFAC parameters. Indeed, it is not clear what UNIFAC parameters were used by Gottlieb and Herskowitz, nor how they were obtained. It is likely that if the data for poly(dimethylsi1oxane) solutions are examined in toto, and if relatively small adjustments are made in the (in any case uncertain) UNIFAC energy parameters, the calculated and experimental results could be brought into much better agreement. Since publication of the original UNIFAC work (Fredenslund, Jones, and Prausnitz, 1975), numerous articles have suggested modifications, and new tables of parameters have been proposed. These articles underline the approximate nature of UNIFAC which is based on essentially arbitrary choices, concerning what shall be defined as a group and concerning what shall be used as the size of a unit cell. The arbitrariness of these choices, as well as uncertainties in the parameters, can be reduced only through global comparison with extensive and reliable experimental data. Literature Cited Fredenslund, Aa.; Jones, R. L.; Prausnitz, J. M. AIChE J. 1975, 21, 1086. Lichtenthaler, R. N.; Llu, D. D.; Prausnitz, J. M. Ber. Bunsenges. W y s . Chem. 1974, 78, 470.

Department of Chemical Engineering University of California, Berkeley Berkeley, California 94720

0 1982 American

Chemical Society

J. M. Prausnitz