推力室热力计算大作业文档格式.docx
- 文档编号:22659747
- 上传时间:2023-02-05
- 格式:DOCX
- 页数:17
- 大小:188.10KB
推力室热力计算大作业文档格式.docx
《推力室热力计算大作业文档格式.docx》由会员分享,可在线阅读,更多相关《推力室热力计算大作业文档格式.docx(17页珍藏版)》请在冰豆网上搜索。
#include<
stdio.h>
math.h>
#defineNC11.2738
#defineNH40.8672
#defineNO25.8632
#defineNN5.6773
#defineNCl5.7692
#defineNS2.4177
#defineNAl1.8532
#definep70
#defineIp0-2400
#defineR08.314
intmain(void)
{
//T=2800k
inti;
doublenSO2=NS,nAl2O3=NAl/2;
/*第一次近似计算*/
doubleKp=7.0080,KpH2O=0.2233e-1,KpH2Oz=0.2091e-1,KpH2=0.6649e-2,KpHCl=4.5983e-3,KpNO=0.8786e-2;
doublehAl2O3=-1317.62,hH2O=-127.652,hH2=81.5682,hHCl=-8.20359,hO2=90.35295,hNO=178.2121,hCl=176.07,hCO2=-252.515,hCO=-24.3009,hN2=85.56553,hSO2=-158.738,hOH=124.7157,hH=270.2106;
doublenH=0,nCl=0,nO2=0,nOH=0,nNO=0;
//令浓度较小组分摩尔数为0
doublenN2,nHCl,NOz,A,B,nCO2,nCO,nH2O,nH2,ng;
for(i=0;
i<
10;
i++)
{
nN2=(NN-nNO)/2;
nHCl=NCl-nCl;
NOz=NO-2*nSO2-3*nAl2O3-2*nO2-nOH-nNO;
A=NOz-NC;
//NOz为NO'
的意思
B=(NH-nHCl-nOH-nH)/2-A;
/*计算nCO2*/
nCO2=0.5*(sqrt((Kp*B+NOz)*(Kp*B+NOz)+4*(Kp-1)*A*NC)-Kp*B-NOz)/(Kp-1);
/*代入(4-57)*/
nCO=NC-nCO2;
nH2O=A-nCO2;
nH2=B+nCO2;
/*计算ng*/
ng=nCO2+nH2O+nCO+nH2+nN2+nHCl+nSO2+nO2+nOH+nNO+nH+nCl;
/*第二次近似计算*/
nO2=(nH2O/nH2)*(nH2O/nH2)*KpH2O*KpH2O*pow(p/ng,-1);
nOH=nH2O*KpH2Oz*pow(p/ng,-1)/sqrt(nH2);
nH=sqrt(nH2*KpH2)*pow(p/ng,-0.5);
nCl=nHCl*KpHCl*pow(p/ng,-1)/nH;
nNO=sqrt(nN2*nO2*KpNO);
}
doubleIp1=nAl2O3*hAl2O3+nH2O*hH2O+nH2*hH2+nHCl*hHCl+nO2*hO2+nNO*hNO+nCl*hCl+nCO2*hCO2+nCO*hCO+nN2*hN2+nSO2*hSO2+nOH*hOH+nH*hH;
printf("
T1=2800K\n"
);
nAl2O3=%lfmol/Kg,\nnH2O=%lfmol/Kg,\nnH2=%lfmol/Kg,\nnHCl=%lfmol/Kg,\nnO2=%lfmol/Kg,\nnNO=%lfmol/Kg,\nnCl=%lfmol/Kg,\nnCO2=%lfmol/Kg,\nnCO=%lfmol/Kg,\nnN2=%lfmol/Kg,\nnSO2=%lfmol/Kg,\nnOH=%lfmol/Kg,\nnH=%lfmol/Kg."
nAl2O3,nH2O,nH2,nHCl,nO2,nNO,nCl,nCO2,nCO,nN2,nSO2,nOH,nH);
\nIp1=%lfkJ/kg\n"
Ip1);
/*----------------------------------------我是华丽丽的分割线-----------------------------*/
//T=2000k
Kp=4.7410;
KpH2O=0.2892e-3,KpH2Oz=0.1378e-3,KpH2=0.2631e-5,KpHCl=1.9697e-6,KpNO=0.3926e-3;
hAl2O3=-1460.81,hH2O=-169.323,hH2=53.1347,hHCl=-37.2656,hO2=59.38296,hNO=148.5028,hCl=159.096,hCO2=-301.93,hCO=-53.7004,hN2=56.33483,hSO2=-205.823,hOH=96.07448,hH=253.5697;
nH=0,nCl=0,nO2=0,nOH=0,nNO=0;
doubleIp2=nAl2O3*hAl2O3+nH2O*hH2O+nH2*hH2+nHCl*hHCl+nO2*hO2+nNO*hNO+nCl*hCl+nCO2*hCO2+nCO*hCO+nN2*hN2+nSO2*hSO2+nOH*hOH+nH*hH;
\nT2=2000K"
\nnAl2O3=%lfmol/Kg,\nnH2O=%lfmol/Kg,\nnH2=%lfmol/Kg,\nnHCl=%lfmol/Kg,\nnO2=%lfmol/Kg,\nnNO=%lfmol/Kg,\nnCl=%lfmol/Kg,\nnCO2=%lfmol/Kg,\nnCO=%lfmol/Kg,\nnN2=%lfmol/Kg,\nnSO2=%lfmol/Kg,\nnOH=%lfmol/Kg,\nnH=%lfmol/Kg."
\nIp2=%lfkJ/kg\n"
Ip2);
doubleTf=(Ip0-Ip2)*(2800-2000)/(Ip1-Ip2)+2000;
printf("
\n内插法求得燃烧室温度Tf=%lfK"
Tf);
/*----------------------------------------计算出Tf=2400K-----------------------------*/
//Tf=2400K^_^
Kp=6.0390;
KpH2O=0.3634e-2,KpH2Oz=0.2573e-2,KpH2=0.2516e-3,KpHCl=1.8055e-4,KpNO=0.2410e-2;
hAl2O3=-1378.49,hH2O=-148.809,hH2=67.12096,hHCl=-22.8753,hO2=74.68854,hNO=163.2842,hCl=167.6069,hCO2=-277.412,hCO=-39.0831,hN2=70.85828,hSO2=-182.387,hOH=110.2098,hH=261.8902;
doubleIpf=nAl2O3*hAl2O3+nH2O*hH2O+nH2*hH2+nHCl*hHCl+nO2*hO2+nNO*hNO+nCl*hCl+nCO2*hCO2+nCO*hCO+nN2*hN2+nSO2*hSO2+nOH*hOH+nH*hH;
\n燃烧室平衡组分\nnAl2O3=%lfmol/Kg,\nnH2O=%lfmol/Kg,\nnH2=%lfmol/Kg,\nnHCl=%lfmol/Kg,\nnO2=%lfmol/Kg,\nnNO=%lfmol/Kg,\nnCl=%lfmol/Kg,\nnCO2=%lfmol/Kg,\nnCO=%lfmol/Kg,\nnN2=%lfmol/Kg,\nnSO2=%lfmol/Kg,\nnOH=%lfmol/Kg,\nnH=%lfmol/Kg."
\n总焓\nIpf=%lfkJ/kg\n"
Ipf);
//根据2400K计算的总焓和所给Ip0相差不大
/*----------------------------------------计算总熵-----------------------------*/
doubles0Al2O3=310.181,s0H2O=273.967,s0H2=194.833,s0HCl=252.735,s0O2=275.851,s0NO=279.948,s0Cl=211.43,s0CO2=320.652,s0CO=265.813,s0N2=258.775,s0SO2=356.278,s0OH=248.849,s0H=158.075;
doublesAl2O3=s0Al2O3,sH2O=s0H2O-R0*log(p*nH2O/ng),sH2=s0H2-R0*log(p*nH2/ng);
doublesHCl=s0HCl-R0*log(p*nHCl/ng),sO2=s0O2-R0*log(p*nO2/ng),sNO=s0NO-R0*log(p*nNO/ng);
doublesCl=s0Cl-R0*log(p*nCl/ng),sCO2=s0CO2-R0*log(p*nCO2/ng),sCO=s0CO-R0*log(p*nCO/ng);
doublesN2=s0N2-R0*log(p*nN2/ng),sSO2=s0SO2-R0*log(p*nSO2/ng),sOH=s0OH-R0*log(p*nOH/ng),sH=s0H-R0*log(p*nH/ng);
doubles=nAl2O3*sAl2O3+nH2O*sH2O+nH2*sH2+nHCl*sHCl+nO2*sO2+nNO*sNO+nCl*sCl+nCO2*sCO2+nCO*sCO+nN2*sN2+nSO2*sSO2+nOH*sOH+nH*sH;
总熵\ns=%lfJ/Kg.K\n"
s);
/*----------------------------------------计算比热比k-----------------------------*/
doublecpAl2O3=143.9089,cpH2O=52.18258,cpH2=35.57987,cpHCl=36.3641,cpO2=38.72138,cpNO=37.15544,cpCl=21.20883,cpCO2=61.82268,cpCO=36.77861,cpN2=36.56507,cpSO2=58.88178,cpOH=35.84491,cpH=20.80102;
doublecp=nAl2O3*cpAl2O3+nH2O*cpH2O+nH2*cpH2+nHCl*cpHCl+nO2*cpO2+nNO*cpNO+nCl*cpCl+nCO2*cpCO2+nCO*cpCO+nN2*cpN2+nSO2*cpSO2+nOH*cpOH+nH*cpH;
doublem=1/(ng+nAl2O3);
doubleR=R0/m;
doublek=cp/(cp-R);
比热比\nk=%lf\n"
k);
/*----------------------------------------计算燃烧室理论特征速度c*-----------------------------*/
doublecx=sqrt(R*Tf)/0.6562;
燃烧室理论特征速度\nc*=%lfm/s\n"
cx);
/*----------------------------------------计算出口截面温度Te-----------------------------*/
doubleTe=pow(1.0000000/p,(k-1.000000000)/k)*2400;
出口截面温度\nTe=%lfK\n"
Te);
/*----------------------------------------计算1000K和1100K下燃烧产物平衡组分-----------------------------*/
//1000K
Kp=0.7254;
\nT=1000K时\n出口截面燃烧产物平衡组分\nnAl2O3=%lfmol/Kg,\nnH2O=%lfmol/Kg,\nnH2=%lfmol/Kg,\nnHCl=%lfmol/Kg,\nnCO2=%lfmol/Kg,\nnCO=%lfmol/Kg,\nnN2=%lfmol/Kg,\nnSO2=%lfmol/Kg."
nAl2O3,nH2O,nH2,nHCl,nCO2,nCO,nN2,nSO2);
/*计算焓*/
hAl2O3=-1593.1,hH2O=-215.653,hH2=20.83535,hHCl=-71.1597,hCO2=-368.551,hCO=-88.7816,hN2=21.62586,hSO2=-262.681;
doubleIpe1=nAl2O3*hAl2O3+nH2O*hH2O+nH2*hH2+nHCl*hHCl+nCO2*hCO2+nCO*hCO+nN2*hN2+nSO2*hSO2;
\n总焓\nIpe1=%lfkJ/kg"
Ipe1);
/*计算熵*/
s0Al2O3=182.387,s0H2O=232.8726,s0H2=166.2406,s0HCl=222.8405,s0CO2=269.3539,s0CO=234.9577,s0N2=228.2292,s0SO2=306.3339;
sAl2O3=s0Al2O3,sH2O=s0H2O-R0*log(nH2O/ng),sH2=s0H2-R0*log(nH2/ng),sHCl=s0HCl-R0*log(nHCl/ng);
sCO2=s0CO2-R0*log(nCO2/ng),sCO=s0CO-R0*log(nCO/ng),sN2=s0N2-R0*log(nN2/ng),sSO2=s0SO2-R0*log(nSO2/ng);
doubles1=nAl2O3*sAl2O3+nH2O*sH2O+nH2*sH2+nHCl*sHCl+nCO2*sCO2+nCO*sCO+nN2*sN2+nSO2*sSO2;
\n总熵\ns1=%lfJ/Kg.K\n"
s1);
//1100K
Kp=1.056;
\n\nT=1100K时\n出口截面燃烧产物平衡组分\nnAl2O3=%lfmol/Kg,\nnH2O=%lfmol/Kg,\nnH2=%lfmol/Kg,\nnHCl=%lfmol/Kg,\nnCO2=%lfmol/Kg,\nnCO=%lfmol/Kg,\nnN2=%lfmol/Kg,\nnSO2=%lfmol/Kg."
nAl2O3,nH2O,nH2,nHCl,nCO
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 推力 热力 计算 作业