LiCl物性计算公式MATLAB版

8、热导率
function lamda_LiCl=conductivity(t,x) lamda_H2O=refpropm('L','T',t+273.15,'Q',0,'water'); % lamda_H2O 为水的热导率，调用标准物性计算软件 refpropm 8.0 rho_LiCl=density_licl(t,x); R=(10.8958e-3)+(-11.7882e-3)*x; Xeq=x*rho_LiCl/42.39; lamda_LiCl=lamda_H2O-Xeq*R; end
6、吸收热
function absheat=absheat_licl(t,x) %已知 LiCl 溶液温度 t 和浓度 x 求溶液的吸收热 %吸收热单位是 kJ/kg hdil=differential_licl(t,x);%微分稀释热 p=(pressure_licl(t,x))/100000;%此处 p 的单位化为 bar hv=(0.1638*t-30.44969)*p+(1.88163*t+2500.87);%饱和水蒸气焓，对应溶液温度 hfpie=4.19*t-0.0419;%纯水焓，对应溶液温度 hfgpie=hv-hfpie;%纯水的汽化潜热（汽化焓） ，对应溶液温度 absheat=hdil+hfgpie;%吸收热的计算 end
பைடு நூலகம்
3、扩散系数
function temp=diffusion_licl(t,x) %已知溶液温度 t 和浓度 x。 %扩散系数单位 m2·s-1 R=8.314; A=0.11353e-16; Vc=1/17873; T=t+273.15; Tc=647.096; rouc=17873; T0=T/647.226; alfa=[1.99274064,1.09965342,-0.510839303,-1.75493479,-45.5170352,-674694.45]; beita=[1/3,2/3,5/3,16/3,43/3,110/3]; tem1=0; for i=1:6 tem1=tem1+(alfa(i)*((1-(T/Tc))^beita(i))); end tem2=rouc*(1+tem1); tem3=1/tem2; H=[1.000,0.978197,0.579829,-0.202354]; temp1=0; for i=1:4 temp1=temp1+H(i)*(T0^(-1*(i-1))); end vis0=(T0^0.5)/temp1; tao=1-T0; B=[1.9937718430,1.0985211604,-0.5094492996,-1.7619124270,-44.9005480267,723692.2618632]; beita=[1/3,2/3,5/3,16/3,43/3,110/3]; temp2=0; for i=1:6 temp2=temp2+(B(i)*(tao^beita(i))); end temp3=1+temp2; temp4=(1/T0)-1; temp5=0.5132047+0.2151778*temp2-0.2818107*power(temp2,2)+0.1778064*power(temp2,3)0.0417661*power(temp2,4)+0.3205656*temp4+0.7317883*temp4*temp2; temp6=-1.070786*temp4*power(temp2,2)+0.4605040*temp4*power(temp2,3)0.01578386*temp4*power(temp2,5)+1.241044*temp2*power(temp4,2)1.263184*power(temp2,2)*power(temp4,2)+0.2340379*power(temp2,3)*power(temp4,2)+1.4767 83*temp2*power(temp4,3); temp7=0.4924179*power(temp4,3)*power(temp2,3)+0.1600435*power(temp4,3)*power(temp2,4)-
0.003629481*power(temp4,3)*power(temp2,6)0.7782567*power(temp4,4)+0.1885447*power(temp4,5); temp8=exp(temp3*(temp5+temp6+temp7)); temp9=(55.071*power(10,-6))*vis0*temp8; D0=(R*T*A*((Vc)^(2/3)))/(temp9*tem3); p1=0.52; p2=-4.92; p3=-0.56; temp=D0*(1-((1+(((x^0.5)/p1)^p2))^p3)); end
9、表面张力
function temp=surface_licl(t,x) %已知温度和浓度求表面张力 %单位是 N/m T=t+273.15; T0=T/647.226; sigma0=235.8; b=-0.625; mu=1.256; sigma=[2.757115 -12.011299 14.751818 2.443204 -3.147739]; alpha_h20=sigma0*(1-b*(1-T0))*((1-T0)^mu); temp=alpha_h20*(1+sigma(1)*x+sigma(2)*x*T0+sigma(3)*x*T0*T0+sigma(4)*x*x+sigma(5)*x *x*x); temp=temp/1000; end
5、表面蒸汽压
function p=pressure_licl(t,x) %已知溶液温度 t 和浓度 x %表面蒸汽压单位 Pa tao=1-(t+273.15)/647.14; p=exp((-7.85823*tao+1.83991*tao^1.5-11.7811*tao^3+22.6705*tao^3.515.9393*tao^4+1.77516*tao^7.5)/(1-tao))*22064000*((1-(1+(x/0.362)^(-4.75))^(-0.4)-0.03*exp(200*(x-0.1)^2))*(2-(1+(x/0.28)^4.3)^0.6+((1+(x/0.21)^5.1)^0.49-1)*(t+273.15)/647.14)); end
4、定压比热
function temp=cp_licl(t,x) %已知溶液温度 t 和浓度 x。 %求 LiCl 溶液比热容 单位 kJ·kg-1·K-1 T=t+273.15; theta=(T/228)-1; A=[88.7891,-120.1958,-16.9264,52.4654,0.10826,0.46988]; B=[0,0.02,0.04,0.06,1.8,8]; temp1=0; for i=1:6 temp1=temp1+A(i)*(theta^B(i)); end C=[1.43980,-1.24317,-0.12070,0.12825,0.62934,58.5225,-105.6343,47.7948]; f1=0; if x<=0.31 for i=1:3 f1=f1+(C(i)*(x^i)); end else f1=0.12825+0.62934*x; end D=[0.02,0.04,0.06]; f2=0; for i=1:3 f2=f2+C(i+5)*(theta^D(i)); end temp=temp1*(1-f2*f1); end
7、微分稀释热
function temp=differential_licl(t,x) %微分稀释热单位 kJ/kg %已知 LiCl 溶液温度 t 和浓度 x T=t+273.15; T0=T/647.226; H1=0.845; H2=-1.965; H3=-2.265; H4=0.6; H5=169.105; H6=457.850; x=x/(H4-x); delta_hd0=H5+H6*T0; temp=delta_hd0*((1+((x/H1)^H2))^H3); end
LiCl 物性计算公式 MATLAB 版
所有公式均用 MATLAB 语言编写，输入参数 t 均为摄氏温度，x 代表 LiCl 质量分数的 小数表达，即，x 的范围为 0~1。
1、密度
function temp=density_licl(t,x) %已知溶液温度 t 和浓度 x 求密度，适用范围 x= 0~0.56。 %密度单位 kg·m-3 c0=1.993771843; c1=1.09852116; c2=-0.5094493; c3=-1.761912427; c4=-44.90054803; c5=-723692.2619; tao=1-(t+273.15)/647.14; rouw=322*(1+c0*tao^(1/3)+c1*tao^(2/3)+c2*tao^(5/3)+c3*tao^(16/3)+c4*tao^(43/3)+c5*tao^(1 10/3)); nongdu=x/(1-x); rou1=0.540966; rou2=-0.303792; rou3=0.100791; temp=rouw*(1+rou1*nongdu+rou2*nongdu^2+rou3*nongdu^3); end
2、动力粘度 dynamic viscosity
function temp=viscosity_licl(t,x) %已知溶液温度 t 和浓度 x。 %动力粘度单位 Pa·s T=t+273.15; T0=T/647.226; H=[1.000,0.978197,0.579829,-0.202354]; temp1=0; for i=1:4 temp1=temp1+H(i)*(T0^(-1*(i-1))); end vis0=(T0^0.5)/temp1; tao=1-T0; B=[1.9937718430,1.0985211604,-0.5094492996,-1.7619124270,-44.9005480267,723692.2618632]; beita=[1/3,2/3,5/3,16/3,43/3,110/3]; temp2=0; for i=1:6 temp2=temp2+(B(i)*(tao^beita(i))); end temp3=1+temp2; temp4=(1/T0)-1; temp5=0.5132047+0.2151778*temp2-0.2818107*power(temp2,2)+0.1778064*power(temp2,3)0.0417661*power(temp2,4)+0.3205656*temp4+0.7317883*temp4*temp2; temp6=-1.070786*temp4*power(temp2,2)+0.4605040*temp4*power(temp2,3)0.01578386*temp4*power(temp2,5)+1.241044*temp2*power(temp4,2)1.263184*power(temp2,2)*power(temp4,2)+0.2340379*power(temp2,3)*power(temp4,2)+1.4767 83*temp2*power(temp4,3); temp7=0.4924179*power(temp4,3)*power(temp2,3)+0.1600435*power(temp4,3)*power(temp2,4)0.003629481*power(temp4,3)*power(temp2,6)0.7782567*power(temp4,4)+0.1885447*power(temp4,5); temp8=exp(temp3*(temp5+temp6+temp7)); temp9=(55.071*power(10,-6))*vis0*temp8; N=x/(power(1-x,1/0.6)); temp=temp9*exp(0.090481*power(N,3.6)+1.390262*N+0.675875*(N/T0)0.583517*power(N,2)); end