matlab編程天然氣壓縮因子計(jì)算模型-wenkub
2023-06-07 18:26:48 本頁面
【正文】 C3HiC4HnC4HiC5H1nC5H1nC6H1nC7H1nC8H18。 模型描述 已知條件絕對壓力P、熱力學(xué)溫度T、組分?jǐn)?shù)N；各組分的摩爾分?jǐn)?shù)Xi，i = 1～N ；查附表3得到的以下數(shù)據(jù)：58種物質(zhì)的狀態(tài)方程參數(shù)an，bn， ，kn，un，gn，qn，fn，sn，wn ；14種識別組分的特征參數(shù)Mi，Ei，Ki，Gi，Qi，F(xiàn)i，Si，Wi ；14種識別組分的二元交互作用參數(shù)Eij*，Uij，Kij，Gij* 。3 程序代碼function [ Z ] = YSYZ( T,p,x)%計(jì)算天然氣給定組分的壓縮因子% x為天然氣組分，按照CO2 N2 H2 CO CH4 C2H6 C3H8 iC4H10 nC4H10 iC5H12 nC5H12% nC6H14 nC7H16 nC8H18的順序輸入%T 為溫度，單位為K%P為壓力，單位為兆帕N=14。c=[0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,0,0,1,1,1,1,1,1,1,0,1,1,1,1,1,1,1,1,1,0,0,1,1,1,1,1,0,1,1,1,1,0,1,0,1,1,1,1,1,1]。s=[0,0,0,0,0,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]。%特征參數(shù)值M=[,]。K=[,]。F(3)=1。... ,。... ,。... ,。 ,。 ,。 ,。 ,]。... ,。... ,。... ,。 ,。 ,。 ,。 ,]。%參數(shù)Eij、Gij以及第二維利系數(shù)的計(jì)算sum1=0。 for i=1:14 for j=1:14 Eij(i,j)=Ex(i,j)*(E(i)*E(j)^)。 end end B=B+a(n)*ZJCS*sum。 U0=0。 sum2=sum2+x(i)*(E(i)^)。U0=((sum2^2)+2*U0)^。for i=1:(N1) for j=2:N sum2=sum2+x(i)*x(j)*((Kx(i,j)^5)1)*((K(i)*K(j))^)。SUM2=0。 SUM2=SUM2+Cn*(b(i)c(i)*k(i)*(pr^k(i)))*(pr^b(i))*exp(c(i)*(pr^k(i)))。 SUM2=0。endZend4 運(yùn)行結(jié)果對氣樣如下表的天然氣計(jì)算在300K、 氣樣組分氣體組成摩爾分?jǐn)?shù)CO2N2H20CO0CH4C2H6C3H8iC4H10nC4H10iC5H12nC5H12C6H14C7H160C8H180
公司管理相關(guān)推薦
鄂ICP備17016276號-1