火箭發動機原理
2023年推力室熱力計算大作業報告
學號 12151171
學生姓名段毓
2023年5月18日
首先設定兩個假定的燃燒溫度t1,t2,用迭代程式設計方法計算出兩種溫度下的燃燒室平衡組分,再由平衡組分計算出兩個總焓,然後線性插值求出tf。
根據tf,計算燃燒室真正的平衡組分、焓值、熵值以及比熱比。
根據壓強比和燃燒室溫度求出出口截面燃燒產物溫度te,並確定相鄰兩個溫度te1和te2。
分別求出噴管出口截面在te1、te2條件下的平衡組分、焓值、熵值,因為是等熵流動過程,所以用熵的「內插法」計算出真正te。
由te求出出口的平衡組分和發動機的理論比衝。
#include <>
#include <>
#define nc 11.2738
#define nh 40.8672
#define no 25.8632
#define nn 5.6773
#define ncl 5.7692
#define ns 2.4177
#define nal 1.8532
#define p 70
#define ip0 -2400
#define r0 8.314
int main(void)
double ip1 = 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");
printf("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);
printf("\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;
for (i=0;i<10;i++)
double ip2 = 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("\nt2=2000k");
printf("\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.",nal2o3,nh2o,nh2,nhcl,no2,nno,ncl,nco2,nco,nn2,nso2,noh,nh);
printf("\nip2=%lfkj/kg\n",ip2);
double tf = (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;
nh = 0,ncl = 0,no2 = 0,noh = 0,nno = 0;
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 - ncnoz為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);
北航數值分析計算實習大作業
目標 使用帶雙步位移的qr分解法求矩陣的全部特徵值,並對其中的每乙個實特徵值求相應的特徵向量。已知 i,j 1,2,10 演算法 1 輸入需要求解的矩陣 首先輸入需要求解的矩陣。即為下述程式中的void creaa 子程式。2 對上述生成的矩陣進行擬上三角化 為了減少求特徵值和特徵向量過程中的計算量...
計算機網路大作業
141021班,141022班 題目一設計乙個宿舍區域網方案,構建乙個多網段的網際網路絡,使得宿舍之間可以通過該區域網進行網路互聯和通訊。要求每個宿舍內使用乙個區域網,宿舍之間採用路由器進行相連。要求 1.每個區域網網段位址可以自行選取,但是要求在每個宿舍的區域網路中網路主機只能接入到8臺以內 通過...
物聯網與雲計算課程大作業方案
大連理工大本科學生連連 空間間跡追蹤系間 院 系,學連件院學 連連,嵌入式 生姓名,學連連連 學號 201392241 完成日期,年月日 2015?12 27 研究背景,在連算機技連高速連展的今天 人連連空連連跡追蹤和智慧型人連作連連越越重連。體來如今大型的連航系連如美的國gps 中的北斗 洲伽利略...