Peng-Robinson EoS

Calculates the compressibility factor, fugacity coefficient and density of a pure component.

Description

This function allows calculating the compressibility factor, the fugacity coefficient and density of a pure compound for both liquid and vapour phases by means of the Peng-Robinson equation of state (PR EOS). Temperature, pressure, molar weigth, and supercritical properties of the selected compound are needed to be specified as parameters.

It is a Beta version, full version will include multicomponent systems and other equations of state.

 function [Z,fhi,density] = PengRobinson(T,P,Tc,Pc,w,MW,Liquido)
 Parameters: T,P,w,Tc,Pc,w,MW,Liquido
 T: Temperature [=]K
 P: Presure [=] Pa
 Tc: critical temperature [=] K
 Pc: critical presure [=] Pa
 w: accentic factor
 MW: molar weigth [=] kg /mol
 Liquido: if Liquido = 1, liquid fugacity;
           if Liquido = 0 vapor fugacity
 Example:
 [Z fhi density] = PengRobinson(273,2*1.013*1e5,304.21,7.382*1e6,0.225,0.044,1)

This MATLAB function provides an educational tool to predict fluid phase equilibrium and can be easily extended to a multicomponent system.

• MATLAB release: MATLAB 6.5 (R13)

[Farshid Farzaneh]

That's great

[reza laee]

full programm

[Miguel Angel Viascan Chavez]

 

[Ivo Kordic]

Bastante bueno

[Tulio Quijada]

Hola Juan, Mi nombre es Tulio soy Ing Quimico, estoy desarrollando un modulo parecido al tuyo pero con la diferencia que estoy trabajando con compuestos no definidos Pseudo-compuestos(diesel) me gustaria conocer tu programa pero no he podido descargarlo, tambien me gustaria intercambiar ideas contigo.

[Herald Arkenbosch]

Juan, Thanks for the file you shared. I had difficulties getting back ground information on the equation of state; this helps me a lot.

[ibe ugo herbie]

 

[Simon Jönsson]

Excellent easy-to-follow program for educational purpose!

[dharma r maddala]

 

[Miguel D. B.]

Hola Juan, quisiera hacer unas recomendaciones para mejorar el programa:
1. En lugar de usar "roots" sería conveniente utilizar el algoritmo de Ulrich.
2. Utilizar propiedades reducidas para no limitarse un sistema de unidades.
3. También sería bueno utilizar la EOS Peng-Robinson-Strajik-Vera para mejorar los resultados.

[rasoul anwar]

 

[hamed heidarzadeh]

compressibility factor

[Kamran Moradi]

Can solve my article problem , nut not fully.

[mehdi escrochi]

 

[julio gomez]

 

[Cc Cc]

thank you my friend

[Patrick Egan]

Functional and very useful.

[mostafa mafi]

DEar My Friend
If it is possible please send me the multicomponent systems version

[Nand Khatri]

Dear Friend..
Is it possible to give some idea for Binary mixture using general two parameter EOS.
Thanx

[muhammad rizal]

 

[lior kaufman]

 

[Enjou G]

extremely helpful!
how can i use it to determine liquid and vapor molar densities cm^3/mol) and mass densityof water at given T,P?