열역학 상태방정식 계산(매트랩 이용)

소개

화학공학 열역학 과목에서 상태방정식은 여러가지가 있다.
그중 반데르 발스 방정식 Redlich-Kwong ,Soave-Relich-Kwong 3차방정식,비리얼 상태방정식등 여러가지 상태방정식을 매트랩을 이용하여 Saturated liquid Volume 및Saturated vapor volume 등을 구해본다.

목차

없음

본문내용

<매트랩 실행 화면 (eos.m, cubfunV.m, cubfunL.m) >
eos.m
disp(`상태식 계산을 위한 기본 입력 사항`);
T=input(`온도(K) : `);
P=input(`압력(bar) : `);
Tc=input(`임계온도(K) : `);
Pc=input(`임계압력(bar) : `);
w=input(`이심인자 : `);

Tr=T/Tc;
Pr=P/Pc;
R=83.14;

disp(`Ideal gas Equation of State : `)
V1=R*T/P; %기본적인 이상기체방정식을 이용
fprintf(`기체 부피 : %9.4f \n`, V1

disp(`\nVirial Equation of State : `)
B0=0.083-(0.422/Tr^1.6);
B1=0.139-(0.172/Tr^4.2);
Z=1+B0*Pr/Tr+w*B1*Pr/Tr;
V2=Z*R*T/P; %이상기체에서 얻은 부피에 Z값을 곱함
fprintf(`기체 부피 : %9.4f \n`, V2);

disp(`\ncubic state equations: `);
disp(`1) van der waals 2) Redlich-Kwong 3) Soave-Redlich-Kwong 4) Peng-Robinson`);
n=input(`사용할 방정식의 번호를 입력 : `);

%if문을 이용하여 4개의 EOS중 택1
if(n==1) % VDW EOS
alp=1;
sig=0;
epi=0;
ome=0.125;
phi=0.421875;
elseif(n==2) % RK EOS
alp=Tr^(-0.5);
sig=1;
epi=0;
ome=0.08664;
phi=0.42748;
elseif(n==3) % SRK EOS
alp=(1+(0.48+1.574*w-0.176*w*w)*(1-sqrt(Tr)))^2;
sig=1;
epi=0;
ome=0.08664;

참고자료

· 없음