% to calculate the specific volume v of n butane at 500 k and 18 atmp
function [y] = specvol(V)
Tc = 425.2; % data from Perry's Chemical Handbook critical constant
pc = 37.5; % data from Perry's Chemical Handbook { bar*1*10^5Pa/bar*1 atm/10
T = 500 ;
P = 18 ;
R = 0.082050 ;
aRK = 0.42748*(R*Tc)^2/pc ;
aRK = aRK*(Tc/T)^0.5 ;
bRK = 0.08664*(T*Tc/pc) ;
y = P*v^3- R*T*v^2 + ( aRK - P*bRK^2 - R*T*bRK)*v- aRK*bRK ;
end

 Réponse acceptée

madhan ravi
madhan ravi le 13 Juin 2021
function [y] = specvol(V) % change capital V to v

1 commentaire

madhan ravi
madhan ravi le 13 Juin 2021
Ofcourse if v is a vector then vectorize the last line by changing ^ to .^

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Chemistry dans Centre d'aide et File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by