air_physical_.rar_LK 物性_physical_热物性_空气 物性_空气物性

#include<iostream> #include"math.h" using namespace std; const double b10=0.1181193; const double b20=0.265728; const double b30=0.154790; const double b40=0.030323; const double c10=0.0236744; const double c20=0.0186984; const double c30=0; const double c40=0.042724; const double d10=1.55488e-5; const double d20=6.23689e-5; const double beta0=0.65392; const double gama0=0.060167; const double b1r=0.2026579; const double b2r=0.331511; const double b3r=0.027655; const double b4r=0.203488; const double c1r=0.0313385; const double c2r=0.0503618; const double c3r=0.016901; const double c4r=0.041577; const double d1r=4.8736e-5; const double d2r=7.40336e-6; const double betar=1.226; const double gamar=0.03754; double A1=-0.98176,B1=0.227486,C1=-0.037305e-3,D1=0.048302e-6,E1=-0.185243e-10,F1=0.247488e-14,G1=0.124314; double A2=-0.68925,B2=0.253664,C2=-0.014549e-3,D2=0.012544e-6,E2=-0.017106e-10,F2=-0.008239e-14,G2=0.050052; double A3 = 20.786, B3 = 2.825911e-7, C3 = -1.464191e-7, D3 = 1.092131e-8, E3 = -3.661371e-8, F3 = -6.19735, G3 = 179.999; double A4 = 24.99735, B4 = 55.18696, C4 = -33.69137, D4 = 7.948387, E4 = -0.136638, F4 = -403.6075, G4 = 228.2431, H4 = -393.5224; double A5=-2.46342,B5=0.457392,C5=-0.052512e-3,D5=0.064594e-6,E5=-0.202759e-10,F5=0.236310e-14,G5=-0.339830; double omega1=0.021,omega2=0.04,omega3=-0.004,omega4=0.225,omega5=0.344; double Tc1=154.6,Tc2=126.2,Tc3=150.8,Tc4=304.2,Tc5=647.3; double pc1=5.046,pc2=3.394,pc3=4.874,pc4=7.376,pc5=22.05; double Zc1=0.2905-0.085*omega1,Zc2=0.2905-0.085*omega2,Zc3=0.2905-0.085*omega3,Zc4=0.2905-0.085*omega4,Zc5=0.2905-0.085*omega5; double R1=259.8,R2=296.8,R3=208.1,R4=188.9,R5=461.4; void d2crpar(double d,double &omega,double &g1,double &g2,double &g3,double &g4,double &g5,double &R,double &pc,double &Tc) { g1=231.4/(d+1000); g2=755.3/(d+1000); g3=12.8/(d+1000); g4=0.5/(d+1000); g5=d/(d+1000); R=(287048.89+461.4*d)/(d+1000); double M=(8314*d+8314000)/((287048.89+461.4*d)); double X1=7.23125*M/(d+1000),X2=26.95959*M/(d+1000),X3=0.32044*M/(d+1000),X4=0.01136*M/(d+1000),X5=0.05549*M*d/(d+1000); omega=X1*omega1+X2*omega2+X3*omega3+X4*omega4+X5*omega5; double vc1=Zc1*R1*Tc1/pc1/1000000,vc2=Zc2*R2*Tc2/pc2/1000000,vc3=Zc3*R3*Tc3/pc3/1000000,vc4=Zc4*R4*Tc4/pc4/1000000,vc5=Zc5*R5*Tc5/pc5/1000000; double vc=0; double X[5]={X1,X2,X3,X4,X5}; double VC[5][5]={{pow(pow(vc1,1.0/3)+pow(vc1,1.0/3),3),pow(pow(vc1,1.0/3)+pow(vc2,1.0/3),3),pow(pow(vc1,1.0/3)+pow(vc3,1.0/3),3),pow(pow(vc1,1.0/3)+pow(vc4,1.0/3),3),pow(pow(vc1,1.0/3)+pow(vc5,1.0/3),3)}, {pow(pow(vc2,1.0/3)+pow(vc1,1.0/3),3),pow(pow(vc2,1.0/3)+pow(vc2,1.0/3),3),pow(pow(vc2,1.0/3)+pow(vc3,1.0/3),3),pow(pow(vc2,1.0/3)+pow(vc4,1.0/3),3),pow(pow(vc2,1.0/3)+pow(vc5,1.0/3),3)}, {pow(pow(vc3,1.0/3)+pow(vc1,1.0/3),3),pow(pow(vc3,1.0/3)+pow(vc2,1.0/3),3),pow(pow(vc3,1.0/3)+pow(vc3,1.0/3),3),pow(pow(vc3,1.0/3)+pow(vc4,1.0/3),3),pow(pow(vc3,1.0/3)+pow(vc5,1.0/3),3)}, {pow(pow(vc4,1.0/3)+pow(vc1,1.0/3),3),pow(pow(vc4,1.0/3)+pow(vc2,1.0/3),3),pow(pow(vc4,1.0/3)+pow(vc3,1.0/3),3),pow(pow(vc4,1.0/3)+pow(vc4,1.0/3),3),pow(pow(vc4,1.0/3)+pow(vc5,1.0/3),3)}, {pow(pow(vc5,1.0/3)+pow(vc1,1.0/3),3),pow(pow(vc5,1.0/3)+pow(vc2,1.0/3),3),pow(pow(vc5,1.0/3)+pow(vc3,1.0/3),3),pow(pow(vc5,1.0/3)+pow(vc4,1.0/3),3),pow(pow(vc5,1.0/3)+pow(vc5,1.0/3),3)} }; for(int i=0;i<5;i++) { for(int j=0;j<5;j++) { vc+=X[i]*VC[j][i]*X[j]; } } vc=0.125*vc; double Zc=0.2905-0.085*omega; Tc=0; double TC[5][5]={{pow(pow(vc1,1.0/3)+pow(vc1,1.0/3),3)*pow(Tc1*Tc1,0.5),pow(pow(vc1,1.0/3)+pow(vc2,1.0/3),3)*pow(Tc1*Tc2,0.5),pow(pow(vc1,1.0/3)+pow(vc3,1.0/3),3)*pow(Tc1*Tc3,0.5),pow(pow(vc1,1.0/3)+pow(vc4,1.0/3),3)*pow(Tc1*Tc4,0.5),pow(pow(vc1,1.0/3)+pow(vc5,1.0/3),3)*pow(Tc1*Tc5,0.5)}, {pow(pow(vc2,1.0/3)+pow(vc1,1.0/3),3)*pow(Tc2*Tc1,0.5),pow(pow(vc2,1.0/3)+pow(vc2,1.0/3),3)*pow(Tc2*Tc2,0.5),pow(pow(vc2,1.0/3)+pow(vc3,1.0/3),3)*pow(Tc2*Tc3,0.5),pow(pow(vc2,1.0/3)+pow(vc4,1.0/3),3)*pow(Tc2*Tc4,0.5),pow(pow(vc2,1.0/3)+pow(vc5,1.0/3),3)*pow(Tc2*Tc5,0.5)}, {pow(pow(vc3,1.0/3)+pow(vc1,1.0/3),3)*pow(Tc3*Tc1,0.5),pow(pow(vc3,1.0/3)+pow(vc2,1.0/3),3)*pow(Tc3*Tc2,0.5),pow(pow(vc3,1.0/3)+pow(vc3,1.0/3),3)*pow(Tc3*Tc3,0.5),pow(pow(vc3,1.0/3)+pow(vc4,1.0/3),3)*pow(Tc3*Tc4,0.5),pow(pow(vc3,1.0/3)+pow(vc5,1.0/3),3)*pow(Tc3*Tc5,0.5)}, {pow(pow(vc4,1.0/3)+pow(vc1,1.0/3),3)*pow(Tc4*Tc1,0.5),pow(pow(vc4,1.0/3)+pow(vc2,1.0/3),3)*pow(Tc4*Tc2,0.5),pow(pow(vc4,1.0/3)+pow(vc3,1.0/3),3)*pow(Tc4*Tc3,0.5),pow(pow(vc4,1.0/3)+pow(vc4,1.0/3),3)*pow(Tc4*Tc4,0.5),pow(pow(vc4,1.0/3)+pow(vc5,1.0/3),3)*pow(Tc4*Tc5,0.5)}, {pow(pow(vc5,1.0/3)+pow(vc1,1.0/3),3)*pow(Tc5*Tc1,0.5),pow(pow(vc5,1.0/3)+pow(vc2,1.0/3),3)*pow(Tc5*Tc2,0.5),pow(pow(vc5,1.0/3)+pow(vc3,1.0/3),3)*pow(Tc5*Tc3,0.5),pow(pow(vc5,1.0/3)+pow(vc4,1.0/3),3)*pow(Tc5*Tc4,0.5),pow(pow(vc5,1.0/3)+pow(vc5,1.0/3),3)*pow(Tc5*Tc5,0.5)} }; for(int i=0;i<5;i++) { for(int j=0;j<5;j++) { Tc+=X[i]*TC[j][i]*X[j]; } } Tc=Tc/8/vc; pc=Zc*R*Tc/vc/1000000; } void pt2vr(double d,double p,double T,double &vr00,double &vrrr) { double omega=0, g1=0, g2=0, g3=0, g4=0, g5=0,R=0,pc=0,Tc=0; d2crpar(d,omega, g1, g2, g3, g4, g5,R,pc,Tc); double Tr=(T+273.15)/Tc,pr=p/pc; double B0=b10-b20/Tr-b30/pow(Tr,2)-b40/pow(Tr,3); double C0=c10-c20/Tr+c30/pow(Tr,3); double D0=d10+d20/Tr; double Br=b1r-b2r/Tr-b3r/pow(Tr,2)-b4r/pow(Tr,3); double Cr=c1r-c2r/Tr+c3r/pow(Tr,3); double Dr=d1r+d2r/Tr; double Z0,Z1,fZ,fZ1; Z1=1; do { Z0=Z1; fZ=Z0-1-pr/(8*Z0*Tr-pr)+27*pr/Z0/pow((8*Tr),2); fZ1=1+8*pr*Tr/pow((8*Z0*Tr-pr),2)-27*pr/pow(Z0,2)/pow((8*Tr),2); Z1=Z0-fZ/fZ1; } while(fabs(Z1-Z0)>1.0e-5); double vr=Z1*Tr/pr; double vr0,vr1,fvr,fvr1; vr1=vr; do { vr0=vr1; fvr=pr*vr0/Tr-(1+B0/vr0+C0/pow(vr0,2)+D0/pow(vr0,5)+c40/pow(Tr,3)/pow(vr0,2)*(beta0+gama0/pow(vr0,2))*exp(-gama0/pow(vr0,2))); fvr1=pr/Tr+B0/pow(vr0,2)+2*C0/pow(vr0,3)+5*D0/pow(vr0,6)-c40/pow(Tr,3)*exp(-gama0/pow(vr0,2))*((2*beta0-4)*gama0/pow(vr0,5)-2*beta0/pow(vr0,3)+2*pow(gama0,2)/pow(vr0,7)); vr1=vr0-fvr/fvr1; } while(fabs(vr1-vr0)>1.0e-5); vr00=vr1; vr1=vr; do { vr0=vr1; fvr=pr*vr0/Tr-(1+Br/vr0+Cr/pow(vr0,2)+Dr/pow(vr0,5)+c4r/pow(Tr,3)/pow(vr0,2)*(betar+gamar/pow(vr0,2))*exp(-gamar/pow(vr0,2))); fvr1=pr/Tr+Br/pow(vr0,2)+2*Cr/pow(vr0,3)+5*Dr/pow(vr0,6)-c4r/pow(Tr,3)*exp(-gamar/pow(vr0,2))*((2*betar-4)*gamar/pow(vr0,5)-2*betar/pow(vr0,3)+2*pow(gamar,2)/pow(vr0,7)); vr1=vr0-fvr/fvr1; } while(fabs(vr1-vr0)>1.0e-5); vrrr=vr1; } double pt2Z(double d,double p,double T,double &Z00,double &Zrr) { double omega=0, g1=0, g2=0, g3=0, g4=0, g5=0,R=0,pc=0,Tc=0; d2crpar(d,omega, g1, g2, g3, g4, g5,R,pc,Tc); double vr00=0,vrrr=0; pt2vr(d,p,T,vr00,vrrr); double Tr=(T+273.15)/Tc,pr=p/pc; Z00=pr*vr00/Tr; Zrr=pr*vrrr/Tr; return Z00+omega/0.3978*(Zrr-Z00); } void pt2C(double d,double p,double T,double &cp0,double &cp,double &cv) { double omega=0, g1=0, g2=0, g3=0, g4=0, g5=0,R=0,pc=0,Tc=0; d2crpar(d,omega, g1, g2, g3, g4, g5,R,pc,Tc); double vr00=0,vrrr=0; pt2vr(d,p,T,vr00,vrrr); double Tr=(T+273.15)/Tc,pr=p/pc; double B0=b10-b20/Tr-b30/pow(Tr,2)-b40/pow(Tr,3); double C0=c10-c20/Tr+c30/pow(Tr,3); double D0=d10+d20/Tr; double Br=b1r-b2r/Tr-b3r/pow(Tr,2)-b4r/pow(Tr,3); double Cr=c1r-c2r/Tr+c3r/pow(Tr,3); double Dr=d1r+d2r/Tr; double E0=c40/2/pow(Tr,3)/gama0*(beta0+1-(beta0+1+gama0/pow(vr00,2))*exp(-gama0/pow(vr00,2))); double Er=c4r/2/pow(Tr,3)/gamar*(betar+1-(betar+1+gamar/pow(vrrr,2))*exp(-gamar/pow(vrrr,2))); double M0=1/vr00*(1+(b10+b30/pow(Tr,2)+2*b40/pow(Tr,3))/vr00+(c10-2*c30/pow(Tr,3))/2/pow(vr00,2)+d10/5/pow(vr00,5)-2*c40/pow(Tr,3)/pow(vr00,2)*(beta0+gama0/pow(vr00,2))*exp(-gama0/pow(vr00,2))); double N0=