|
|
程序如下：
T1=243.15:303.15;
T2=243.15:303.15;
T3=243.15:303.15;
%汽化 （液态-气态）
es0=101325; %Pa
H1=40640 %J/mol
T0=373.15;
R=8.314 %J/(mol.K) %s水汽的比气体常数
A1=H1./R;
es1=es0*exp(-A1.*(1./T1-1./T0))
plot(T1,es1,'r')
hold on
%固体融化（固态-液态）
%J/(mol.K) %s水汽的比气体常数
es0=610.7 %融化焓
T0=273.15;
H2=6004; %
Vm=(1.963-1.8).*1e-5 %m3/mol
A2=H2./Vm;
%es2=1e6.*(es0+A.*log(T2./T0));
es2=es0-(log(T0./T2))./A2
hold on
plot(T2,es2,'b')
%固体升华（固态-气态）
H3=46004; %升华焓
A3=H3./R;
T0=273.15
es0=610.7
es3=es0./exp(A3.*(1./T3-1./T0))
plot(T3,es3,'g') |
|
回复此楼
+1/436
