[数据结构与算法]Matlab与热工——水的超临界压力下的热力性质

? ? ? ?对于超临界压力水，采用MATLAB 软件调用美国国家标准与技术研究院的 REFPROP 9.0 软件来获取计算过程中涉及到的热物性参数。计算并分析水的传热系数、密度、动力黏度、比热容的变化情况。在临界压力以上，不同压力时都存在一个定压比热容Cp最大的点，该点称为拟临界点，它对应的温度称为拟临界温度。拟临界区内物性随温度变化剧烈。

1、matlab代码

clear;clc
T1=640;
T2=675;
Tstep=0.1;

i=0;
for p=23000:2000:25000
    i=i+1;
    j=0;
    t=0;
    for T=T1:Tstep:T2
        j=j+1;
        t(j)=T;
        rho(i,j)=refpropm('D','T',T,'P',p,'water');
        cp(i,j)=refpropm('C','T',T,'P',p,'water');
        vis(i,j)=refpropm('V','T',T,'P',p,'water');
        tcon(i,j)=refpropm('L','T',T,'P',p,'water');
    end
end


figure(2)
[AX,H1,H2]=plotyy(t,cp(1,:)./1000,t,tcon(1,:));
set(H2,'linestyle','-.')
set(H1,'linewidth',2)
set(H2,'linewidth',2)
set(AX(1),'ytick',[0:50:300])
set(AX(2),'ytick',[0:0.1:0.6])
set(gca, 'LineWidth',1.5)
set(get(AX(1),'Ylabel'),'string','Cp/kJ·kg^{-1}·K^{-1}')
set(get(AX(2),'Ylabel'),'string','k_f/W·m^{-1}·K^{-1}')
hold on
plot(t,rho(1,:)./3,'m--','linewidth',2)
hold on
plot(t,(10e6)*vis(1,:)./3,'k:','linewidth',2)
xlabel('\fontname{Times New Roman}T/K')

? ? ? ?因为作者菌还不会用matlab画多Y轴图，所以这里采用了双Y轴图加hold on叠加画图的方式，同时对它们纵坐标做了调整，才使得四个曲线在同一张图上出现（PS：会的大神可以教教我~~）

2、加工

? ? ? ?同样采用了比较简单的方法，用PPT加工出另外两个坐标，最后得到的成图如下：

?