Matlab Methode van Suave Radlich Kwong

Figuur 1: Methode van Suave-Redlich-Kwong artikel deel 1

Figuur 2: Methode van Suave-Radlich-Kwong artikel deel 2

Methode van programmeren:
%%Modelleren opdracht 1 met onderdeel Ammoniak
%Soave-Redlich-Kwong methode
clear all

%%  NH3 Eigenschappen
Tc=405.7+273;
Pc=112.8*10^5;
M= 17,03052;
R=8.3144;
T=240;
%% Soave-Redlich-Kwong
psi = 0.42748;
Ohm = 0.08664;
w   =  0.1
Tr=T/Tc;
Asrk=(1+(0.480+1.574*w-0.0176*(w^2))*(1-sqrt(Tr)))^2
aTr = Asrk
eps = 0;
o = 1;
n = 1;              
%% Startwaardes
s=10000;
V = linspace (0.0008, 0.25, s)*10e-3;

pe=2.3283*10^6;
a = psi * ((aTr*(R)^2 * (Tc)^2)/Pc);
b = Ohm *((R * Tc)/Pc);


for i=1:s

    p(i)= ((R*Tc)./(V(i)-b))-(a./((V(i) + (eps*b)).*(V(i) + (o*b))));  % Druk SRK
    pi(i)=(R*T)/(V(i));              

end

%% Snijpunten

f = @(V) (((R*Tc)./(V-b))-(a./((V + (eps*b)).*(V + (o*b))))-pe);
ve1 = fsolve(f,1);                            %snijpunt 1

f = @(V) (((R*Tc)./(V-b))-(a./((V + (eps*b)).*(V + (o*b))))-pe);
ve2 = fsolve(f,1*10^-3);                      %snijpunt 2

f = @(V) (((R*Tc)./(V-b))-(a./((V + (eps*b)).*(V + (o*b))))-pe);
ve3 = fsolve(f,1*10^-5);                      %snijpunt 3

%% lijn tussen de snijpunten
V2=linspace(ve3,ve2,s);
for i=1:s

    pt(i)= V2(i)-V2(i)+pe;


end
%% Dichtheid
Pl=(M/ve3)/1000;
Pg=(M/ve2)/1000;

%% gegevens weergeven
hold on
plot(V,p)
plot(V,pi, 'r')
plot(ve1,pe, 'x')
plot(ve2,pe, 'x')
plot(ve3,pe, 'x')
plot(V2,pt,'g')
legend('SRK','Ideale gaswet','Isotherm')
axis([0 0.6*10^-2 0 10*10^6])
hold off

fprintf('Dichtheid ammoniak(l) bij  %g Kelvin en %g Pascal is %g kg/m3 . \n', T, pe,Pl);
fprintf('Dichtheid ammoniak(g) bij  %g Kelvin en %g Pascal is %g kg/m3 . \n', T, pe,Pg);

psrk= ((R*Tc)./(ve2-b))-(a./((ve2 + (eps*b)).*(ve2 + (o*b))));
pid=(R*T)/(ve2);
verschil=((pid-psrk)/psrk)*100
xlabel('Volume')
ylabel('Druk')

Door het uitvoeren van deze opdracht heb ik mijzelf verdiept in de chemie en wiskunde. Om tot het eindresultaat te komen heb ik diverse onderzoeken verricht naar mijn stof: Ammoniak, verdiept in de methode van Suave-Radlich-Kwong en hoe dit tot stand is gekomen en lastige weergave in matlab weten te volbrengen. Deze ervaring zal ik wellicht toepassen in het examenjaar van Engineering, echter denk ik dat deze methode te specifiek is voor de doeleinden die ik op dat moment zal aanleren.