%% 热加速影响, acceleration flow
%% Ack = qw*bb/(G*Cpb*Reb^0.625)*(uw/ub)*(roub/rouw)^0.5;
%% Acj = bb*qw*d/(lamata*Reb^1.63*Prb);
function Accelerationk = Acceleration1(P,G,qn,Tw,Tf,din)
% P ; 试验段压力,MPa
% G : 质量流速, kg/m2s
% qn : 试验段热流密度, W/m2;
% Tw : 壁面温度, C;
% Tf : 流体温度,C;
% din : 试验段内径, m;
roub = refpropm('D','T',Tf+273.15,'P',P*1000,'CO2'); % 主流体温度,kg/m3;
rouw = refpropm('D','T',Tw+273.15,'P',P*1000,'CO2'); % 壁面流体温度,kg/m3
Prb = refpropm('^','T',Tf+237.15,'P',P*1000,'CO2'); % 普朗特数
lamtab = refpropm('L','T',Tf+237.15,'P',P*1000,'CO2'); % 流体导热系数,w/mK
ub = refpropm('V','T',Tf+273.15,'P',P*1000,'CO2') ; % 主流动力粘性系数 [Pa*s]
uw = refpropm('V','T',Tw+273.15,'P',P*1000,'CO2') ; % 壁面流体动力粘性系数 [Pa*S]
betab = refpropm('B','T',Tf+273.15,'P',P*1000,'CO2'); % 体积膨胀系数, 单位 1/K;
Cpb = refpropm('C','T',Tf+273.15,'P',P*1000,'CO2') ; % 流体比热值, 单位 J/(kgK)
Reb = G*din/ub; % 雷诺数
Accelerationk = ((betab*qn)/(G*Cpb*Reb^0.625))*(uw/ub)*(roub/rouw)^0.5; % Kim 热加速无量纲数
评论1