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
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
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
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
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)-
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
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