From 46059bce01597cb3d52e4184c309587361c192ea Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?China=20Hagstr=C3=B6m?= Date: Thu, 27 Oct 2022 16:22:45 -0400 Subject: [PATCH] cleaned up code, ran all SSME.py altitudes --- SSME.py | 7 ++--- combustion_chamber.py | 4 --- nozzle.py | 5 --- plot_data.h5 | Bin 511456 -> 511456 bytes rockettests/16000m/16000_altitude.yaml | 8 ++--- rockettests/16250m/16250_altitude.yaml | 8 ++--- rockettests/16500m/16500_altitude.yaml | 8 ++--- rockettests/16750m/16750_altitude.yaml | 8 ++--- rockettests/17000m/17000_altitude.yaml | 8 ++--- rockettests/17250m/17250_altitude.yaml | 8 ++--- rockettests/17500m/17500_altitude.yaml | 8 ++--- rockettests/17750m/17750_altitude.yaml | 8 ++--- rockettests/18000m/18000_altitude.yaml | 8 ++--- rockettests/18250m/18250_altitude.yaml | 8 ++--- rockettests/18500m/18500_altitude.yaml | 8 ++--- rockettests/18750m/18750_altitude.yaml | 8 ++--- rockettests/19000m/19000_altitude.yaml | 8 ++--- rockettests/19250m/19250_altitude.yaml | 8 ++--- rockettests/19500m/19500_altitude.yaml | 8 ++--- rockettests/19750m/19750_altitude.yaml | 8 ++--- rockettests/20000m/20000_altitude.yaml | 8 ++--- rockettests/20250m/20250_altitude.yaml | 8 ++--- rockettests/20500m/20500_altitude.yaml | 8 ++--- rockettests/20750m/20750_altitude.yaml | 8 ++--- rockettests/21000m/21000_altitude.yaml | 8 ++--- rockettests/21250m/21250_altitude.yaml | 8 ++--- rockettests/21500m/21500_altitude.yaml | 8 ++--- rockettests/21750m/21750_altitude.yaml | 8 ++--- rockettests/22000m/22000_altitude.yaml | 8 ++--- rockettests/22250m/22250_altitude.yaml | 8 ++--- rockettests/22500m/22500_altitude.yaml | 8 ++--- rockettests/22750m/22750_altitude.yaml | 8 ++--- rockettests/23000m/23000_altitude.yaml | 8 ++--- rockettests/23250m/23250_altitude.yaml | 8 ++--- rockettests/23500m/23500_altitude.yaml | 8 ++--- rockettests/23750m/23750_altitude.yaml | 8 ++--- rockettests/24000m/24000_altitude.yaml | 8 ++--- rockettests/24250m/24250_altitude.yaml | 8 ++--- rockettests/24500m/24500_altitude.yaml | 8 ++--- rockettests/24750m/24750_altitude.yaml | 8 ++--- shocks.py | 42 +------------------------ 41 files changed, 147 insertions(+), 199 deletions(-) diff --git a/SSME.py b/SSME.py index 30e9eb8..09b168e 100644 --- a/SSME.py +++ b/SSME.py @@ -101,13 +101,13 @@ if h >= 25000: T_atm = (-131.21 + 0.00299*h) + 273.14 #[K] - P_atm = (2.488 * ((T_atm)/ 216.6)**-11.388)*1000 #[Pa] + P_atm = (2.488 * ((T_atm)/ 216.6)**(-11.388))*1000 #[Pa] ambient_T = T_atm ambient_P = P_atm elif 11000 < h < 25000: T_atm = -56.46 + 273.14 #[K] - P_atm = (22.65*10**(1.73-0.000157*h))*1000 + P_atm = (22.65*math.exp(1.73-0.000157*h))*1000 ambient_T = T_atm ambient_P = P_atm @@ -123,7 +123,6 @@ F = 2188080 #[N] thrust at 104.5% RPL P_e = 13798.5 #[Pa] exit pressure of nozzle u_e = (F-((P_e-ambient_P)*A_exit))/(mdot_f+mdot_ox) - print(u_e) #Calculate composition of O2, and N2 X_N2 = 0.78084 @@ -191,7 +190,6 @@ mdot_Noz = mdot_f+mdot_ox L_Noz = (math.sqrt(1/3.1415 *math.sqrt(A_throat*Area_ratio)) - (1/3.1415 *math.sqrt(A_throat)))**2 #[m] - print(L_Noz) #Call nozzle function Noz_states = nozzle(T_Noz1, P_Noz1, comp_Noz1, A_throat, A_exit, L_Noz, mdot_ox, mdot_f) n = len(Noz_states.T)-1 @@ -238,7 +236,6 @@ u = results_u[1] #velocity at exit of nozzle (m/s) gamma = 1.1 - print(u) P1 = P_e #TODO Noz_states.P[n] PLACEHOLDER, need design exit pressure or design alt T1 = Noz_states.T[n] M1 = u/math.sqrt(gamma*P1/Noz_states.density[n]) diff --git a/combustion_chamber.py b/combustion_chamber.py index 78a52f5..692b74f 100644 --- a/combustion_chamber.py +++ b/combustion_chamber.py @@ -2,7 +2,6 @@ import os sys.path.insert(0,"/home/chinahg/GCresearch/cantera/build/python") -#print(sys.path) import cantera as ct ct.add_directory('/user/chinahg') @@ -23,7 +22,6 @@ def combustion_chamber(P_CC0, V_CC, mdot_ox, mdot_f): gasExhaust = ct.Solution('h2o2.yaml') gasExhaust.TPY = 300, P_CC0, 'O2: 6, H2:1' gasExhaust.equilibrate('HP') - #print(gasExhaust.report()) gasFuel = ct.Solution('h2o2.yaml') #ct.Hydrogen() gasFuel.TPX = 20.15, 227527, 'H2:1' #https://science.ksc.nasa.gov/shuttle/technology/sts-newsref/et.html#:~:text=LIQUID%20HYDROGEN%20TANK,-The%20liquid%20hydrogen&text=Its%20operating%20pressure%20range%20is,to%20the%20left%20aft%20umbilical. @@ -70,6 +68,4 @@ def combustion_chamber(P_CC0, V_CC, mdot_ox, mdot_f): A = 0.159 #CC cross sectional area - #print("\nFinal CC Composition: ", gasExhaust.report(),"\n", u) - return(state) diff --git a/nozzle.py b/nozzle.py index f8a3d79..cabaf92 100644 --- a/nozzle.py +++ b/nozzle.py @@ -98,10 +98,6 @@ def nozzle(T_Noz1, P_Noz1, comp_Noz1, A_throat, A_exit, L_Noz, mdot_ox, mdot_f): P_throat = P_t*(2*gamma-1)**(-gamma/(gamma-1)) T_throat = T_t*(1/(2*gamma-1)) - print("\nMACH NUMBER AT NOZZLE INLET: ", M_CC) - #print("\nTHROAT TEMPERATURE: ", T_throat) - #print("\nTHROAT PRESSURE: ", P_throat) - #Call diverging section, return final gas state (end of nozzle) states_final = nozzle_div(T_throat, P_throat, comp_Noz1, A_throat, A_exit, L_Noz, mdot_ox,mdot_f) return(states_final) @@ -155,6 +151,5 @@ def nozzle_div(T_Noz1, P_Noz1, comp_Noz1, A_throat, A_exit, L_Noz, mdot_ox, mdot states.append(gas.state, x=solver.t) i = i+1 - #print(i, solver.t, gas.report()) return(states) \ No newline at end of file diff --git a/plot_data.h5 b/plot_data.h5 index 624ee156245e6ba106c092a7368249cc848e7c79..11193cf0700cc025c9735305a2877f632f5151c6 100644 GIT binary patch delta 11070 zcmZXac|g?F+Q*ySr2w_uxGxCWH86bHl&6cATXNvrFIp{Z$Np;@mbDk`J>5(*bg zMT&GL{fQ_pH&Zi_EN?17E}4@en!VlvWqQBQ^ZX_U=P%CkobUOZ=WNgU4Og-QuVe@A z*>2IdbJ|>uGHq7^>se0CK0dz5fG?H6rn=yOytZ#GzwnXu$36+V{fFm=pDPOO^va0` zR?Wut`LG(__wKJa*2r#cd3EcDDZ47Hj$>Jy%{HsPzA|@?zL8_NpF;dS3QIs|cgp&{GRw<9w;k~)vm6BO~;yLs_d z`L`Z;Yu|C(-N2y0!t%)HO`@2_=m$Cf>S zSZu`^wYH)#PF94us3d7dmg5a;8-u=s)0VR#*AX5G2O>>?_lE>iGQh)u%fl-cuc6tp^huChl5X7jQMi9aV9dWpqMQTN5w+Z7;Vy936yG zCE$Hswxe(id5P~(-9;)- zR^6OnN_skNRrQKfKWsh9PRp3q_1Q=9-ehhYa@yM77yl6eva}EWxvnblNwTQLS_fEl zshmt}S8H3cBnaJV3-3)y^TK-z&I3y7EQAwp;h?Y^&O%pTtaGqjiAAmq9#^;tXsKbYHR7#YD97)??*$ayoK0 zz=55#Ri^|~($j%mdwz5E%cI3icI8pC39|#e$+~McL^k+ZQvR}rcB|5BV8hk@>yjng z&TI{$4yxkr*vhtXk{4gDbi3V!N^RkN;&0N2I2BigF4wU*zq8(QuB zs3%R{AU8D|0$HHbl^kwhw`S~~bZ>DdPmpDe5^Xj6f+6n+*bl$pZT`{)XY`WKKQHW5_wnW3UyM zxQMAZ#+2G$(~g?$)&@U!d-(D9PvRcz#yzMN-E^36KvdaHwaV=t zz(dBTCp?P&33FF3!Bk87l0&KX7jv!BdH|c&K^*4z?fIUI_A_GfgKwvPJH{JvsAfZm z4~FUzIs);{$c3Lh{gnrzvZakFS$Q}_Xad^n?Gjq>Jays=W}b5Dow0tcy_wxM#QcvT{R5u^^J}jic+Nl1gPGT7 z`dscoEtpFaj&2cmDMj;yuwot)@_fIsfabfqQ_WdKfZYY~ekz(H?YSabub-5f_@{#b z{h!M(9o4`a&}~D2^JkBHcfHYWJwIXfe+#qf0wybG-USPKz#XWR{ja2P2`V)!w4S$@ zD-Yk}auqQcqMU(*mZ{8Lt^0;}!`28d9Y+mx99#bpcOh#-cOm>EwxlfYw=AzaMzc|3 z|3v)Nt7}bT5FWPWHJ`ECgHYjVZ?H*F!qhv8e-@+g&8T$9e#{&75gREfAHk>Hol43V znDzQOI-Hw5c)(_Hy&x1S4!;4JF zSEb^ASu$alg_15y^->v8_HwF-)YLa7&!x#T_t;`D1_#(K zJacQUH{>JDg#0-8)uEpc1>{fNZ~wS=v?pZg^kk(&C0NiaZc8kHURZ(+b+yXGDlVnl zvRs*XS27f}sV(l&;)ETp%v`nI1DbR+ms`E>iY4#Ht^J&*ct^g`l43yrkoA3~~K zmaW93&XoHtK0UtEXYAO)J74bh&QNdEN16@I`Omkvuk#1%tdwD~2ikeCYL{g@?n1?7 ziAe{jT$CLJ6c?q73UNnftG{;AV_EibstDrX&OYMKmwI5dNPD_-O+owZm!A2Q8Nc21 zbH(b;-i&KC8)Ce9x@D^)PuZ=5CdW4Y-rs|f{n;7yT~v&IM=2XIAWg^mNbqH&y&n{qr&c3<3-iR8z=;>Y`w z(H7}Ceti4xmTh&_e?duSR<*7e6xwnf6W1f2L zvyM)0LbnYO_UXCfqvF5zS_iMsTsZgskB@_Qp6Lc}3Jg6P75}%SJEH8DCeFdbQ3+it z!$?}H$7S9R9z)!A>>ZkIWa!G7w*uBw;qkR+vCq7Si8VWwBPVd4=Pn$!rqmns4Y9Ei z=*6vXWX%ESXOu`&L2@!I=OHIZP7{OGIqB>iT2kM-4OLOK;4JH+~wlH0~4jZ0+Rc z^}3{cW4dj4^O!>?Z2;z@!X4eO6?tIt(rj5GFfqLmx_3DDm@=wVDpE?BQj6QplQa8+ z`2sn!zc_RGMVKP(>C7e@M<=JIZD-2l$cS;j%=4yn+Ysgb_bz-N2g;;ZPF5xqcu*>S z6?pqhyauxVPKn);uX9F4?1sw`r5v}dx*;Y7f0?t@;lI$@9DuU`aa+P|a_Io=$W;4O z>EEqKb#&40{sZh<`rx0OR!S}#YL&8T$PK{sRb%{Ui#Q9jgha86hK?q*B^&DvT62C# ztIG_A)vV1C(bEQ&ppOTNzYk~d^`Knw_mv8XihK?5X&8EFY=Cf>2*)>tFw&Y&)uYWw zI1E=7z4=3fAb=)GK?ZmLrD?$i7ysn|rj01_cPI5~TLVU~$!4~%BfGz&sIn_kdrDo- zSwzIWxC5oQm*`3`+{XZ?4%A-gGY63FK&Kio&`?mH5+%b7*pR5FPdmz)C`cef-*%MS z(UPH$JId`5oRJ7&qyd-4cr=bPm{}b6+Cg>QSOP=qfYTN{c3kaF;qg?OBu&%)R$a|jOaIoqqh!C&P32*^QZy@80u_xZihE;~VYX2pj@A1q zzIUu#N37v^wq=qx!wIxfx%?42K%wg+8aN8+(rB=m3zd+~5JP*MAKB{ZEi`cCeU8DfNIj2&lo+W_xl=InW_N8dhBJHLv z9YM5QfIxX1SE>k<;R3YJ7!OeT%U?SaL{*KMM8nqvx4l?b&~K`TVEQbsROI7iG7uq# z*p@_jpU!yIh!bik!tlAIIaDC8%||n&-H>swlP0Vx-XdQNb-#?3LJ%2w87__$su(vz`?gmD*52cgZ40|cJS@iG92mUYr;td9>Vs$p2BC<=N-oP zNLVc>A{U?s{-$hwj?RnsohGDSQ7PIZHAweM8FJQj?)$%0DC zO>ja(zEk~51CPQPdjV(kx0R$AC8W#m!X9Zi>AibI8ik!#O`}H}Is6`?3uFo&LYDq= zbW%&Q(awD-Cd^Z`2p%_>f+_75>|!!OkP8r|_;96!CTDz2L8#EqKR?c>IR9`^MZ^2g zdWhdkPNX%0 zh%1%$z0n+1xKl8C6w)Lm_T@r_?~{rB&>5uU01P$31F**d2p7PfWsbk?vD(3N(3i)VT+;#M|_K<1Bv(+kym51v+vh8 zNjV13<}?%5X^i-@-%4}=6ZwfJ9j%e!3nh7#WHMHP70ohM)8LVaSx19MdvfJA8@K;^ z@O~Q!Q4>$k5pg;5HjAG8AhU%ta`0_6VVA~YA@dJ1pAcrOs9)O&0i@lojmsl-n98>y zJ49mfAh(kTQ0nsqP+?Acf7d=yRZ9*KF4gkyCx^Jr--_!FYl`w~$4Rc$KH{8`Ow1`( z)jr~!CZEOuiGJ-UA=Wr7)KYRO68L2T%aW?GRug{&eH`4$ze5!KG4j#k$>Y7z25C24m%c>D1vvB~xKcax z24Nhdu{Am9ijSkd4cOpty;CV&@R$kCTT zs$;N5Lzd!~}$71E#?y{ zx+qydE=4hZqd@Syl#R9DDV9lwFO%%IB~$!1Kw&N+-#Y|47QNC>6z`EFGPF$3uLJ_p zZdmIof*p(gY&80@Seo=UE87TFWLz7eN`Ltm1w~ZV;w&1zrr0>AuHv%Iq9_O7F42+a z?rzS=;g`>8H7wtZlWH7d=>pOmCzSDaXd*3?XVnu9Gs1{*?pXRT+9D&dbg^X8ib0A( zxk`@GG93?l{#<{89Dw&Ehf@1b5g&5+z>X1*r(W=J{XBX^ny!-k?;wdZMbi1-OPV(7 zj0@)Ne$*1D6IL_i`6cS_WD_;aDSC1uMc?95CFK^lU@EagC3$dyJ9fBBToYu-%kKjK zY4^n3`<*ahDvzgHs39ZqgG(~1`8aisUA05>Zo?9XJJ(y^X^7m zsqLl47Ve>{vBgC_$u@e{0zW3mb+@U71-J>gRQxSinkeO0)AJBPq0c`p@Fogrp|^a2 z3ze`JAeMlo%|G?DvtaPZap17Elc@2kjPJtse-Hj>SHgs;eAx&E8Zt89{Usy6(G1{> zib==d6A?Ou)8G@urJ0dXMOt&|{vjmF`AS2dAr{*+&jwPKqlY*MAMualqZVv`3X=)+C5p;M3SE94kKQM0Gd?1Y{4=?sp}QEf&TI?nnpe)TYfktlN(?x zVg9ldO$g5?JVF<-kUXIlb6S2%>yRp$m{b;4W0uhPCt((r66GYJd}%o-k#3 zb9#|-)#Hfi%gj_nZMaF0nK>e&@AXQE%kR6^I_G)q{&m;hXYKF)?6v1}?$zv|tJy(& zw)&VGIqfe0O#9WKT0T>wuKP}m{M3FwC?sfSn@huY4$QLs*mIQWz_B-@FC2+%Kk3wO zHr+=>Y_wMTy}YAguqoeGlv4ZBgzyrZ<3!d*Z@bN0Tg{wd_IC`b=P*Iu;qnjU-@pdG z^lJ}nP-buTOemRYVSn0tH)1ofk&AK?-^g}6YTVpnZVbn1xsDS75P|LT|Iyer(>peo6uKN&Q;LXLIGi&zRiprauJ$vWJ$04-`uS~Wte~thUzq1Q0d>KI|V#W}sUDXB-?pphw&% z>Uevv(7^7WRDHMA6EH=5>;!lrbN2S@uD!OOSM9tx=Byi_8h62(Y(G0R1T{FyhoY-d z6ZLQ9FsXEFiM?t#Jyn5AP1%k^FC&0WosQCx2w+n{9jLZ<86Co=8(*u^@=n38ZRG&^ zZ#wa5&SiT6Ewz+|qBhu{>- zwkJ+;;m48qMgfxrP7^X7dNbI=J#raOb!xILYYUTgj#+oz*~^n`mF`24JCalsv%h+; z?Z}sV_JyzfD_PvTdIcDIscfc&RomyzB}RSpy(c!0zIg~V+OBWA>tT)86hwo z1sRwfsIPRDi+*_Mt-;X4qnxz}o=){pc~I5coA^@~-;MT&e}o45uIjSB=p&D5;b*uY zc5_>gsD1Fyy|z;oNqf(~=mw>ay3K@RFCTi(Uc7}(9(cJ@?N)-*2i>*+FLQuXx$0YT z9q;arF#9YFf45fGu^YG-=e^b>HD`k-$zt(QHRHW2{bua?a1clW+J9+H>f}aZU*Ml$ zuiPIBc}L&@#09@m8m0lNXYQcpO}5mV>d9^R(C7=j68G8ls=YQZ$G6UI4cxdh z&#j}i}mm_y7gJ1y7XNrneJ3Ej!ECN%~MQu zIpDq1<=s4~3?HKUd0WE;OF{LZL8?9Cx8ldkowEJQlx*L9H3T)NLOZ7Gk)PQv{lZjj zFTTdM>y)X)Ui^#LmF8^6`s+ejaPu)`{ktBN2{XP)8n$H}Q+7JxSF5MRlhW`Z$}T%v z4D1TZS#ysa9r3*zC9f%W{zDZwe1Ah<12*t&uWEb1?*wI!`rSoYQ4g_U{P0Jz)%-W7 z^6g~Kb=RpFz6<=n$rv}$Gg;1?2bxr7Duw*=Ca`A7pIa>?! zH6eO<9Pe*J_#l9xr?xaS#2#6K^B(TgX>_fPp1it`tjX7BKD8OVww=E> z2;JbuOYu{wfkq8wDm_SW^EUfbtKGqjs|u5?C2DsFCuoh)-Uu~ekqLBa87+x4P)M(} zL0Q=PcvPSNVyf-Ezv6)Z{MAevyzhhjmD4>*|I~dbvgrr2z8MG77F$ zV?ORbaL2$F8O-_PG?kb7oF}K@L!8}4j}L1F&Y=yvj9D|@jZd`oDc zz5ndOx}Idy#fP_#H)fkxS3Ck)RZ05dSa&iWo4rRPgh5V~c*P6jaSZCGYB5j~UXuqY zuYhuHB0o?$i3lmYbFk-Wy~fe_1X$uk>UET(@tHrF%RLdz@LTOpyxe9PM{`%hn9fU! zKlGSSh7Uz^{>YIZ9S3v06*bGVTe~sSp-U~F$uL^pOzL%1i!+xxGD`?cXE7mn`kc3@ z(}VfgrA5!7Fd{*a4^sgN+kL70dVTtZ^Hob2a>K1hPi^1r32FEc*^go4Hn5XQ*?_19J5|a)^n2|b6IRX;*%q~?xfItBV|Z@w6J;swuvq+F*?T3CciK{-vhyWDd!#nXK?9UZAE0Q4NZ(&N#F9_+eV;8i#w}d3kYn37wr-Pz zKu>DjN7d}GRQi4N<+4yvPt80zZt{M2YHX7$2#vN#b>s`C(XdBieMc@fNEiQNbJW(K zCcIFmR`509g$iE%Lg2x6ADVTGr;jW$FJQbau|++-D%N?-pa1!csq`CiV5|jnpB|c$ z(s8{T9WU>9zLyGof4iruoB#}!JyHc+kwnxLKq8ZNN$rPXf#egKw)rROY@Kx2HBZuS zpDmiuH+mkEHoTfL+5e;`>38Czg5LM3^lR^X^Tr>0ZPRQ2r|+}#o=XU8{S%vpf0s0ekz*Z)jt-jQiCnJx z(3F4hTT$?SAm2FtV))4C-H^HC4?iQw*e&m=vyB*8-!U6Y0y}1jy4yq=!`o(4;%(v7 zux(<9wlH?gV7Rn+a*Uoh9v1&u{A^~;o^p2hTfIG5bsyzdg(-CFndFUNUDC>L@lnCN;`EW2ep&Q;hgU6hH9GF!Xd-i&z!NSx#tlH)m{D_Z*(#*7D##&jG1 z7BjwoJoTGDpY~+bebkHXz}WfYb<4r%|D#t*179~rj%V*EqlPjS*OzU@yJfOXC^KwK zsbxf)v7NWX5QxatdJUjXZ)uFmb^c=>j?J0HNGognQm_5I)?;exKI-_hm@s>8%zhV0 z&mXL)yZPY{e_hp6hk*k%@kkv0*v@*HZN=xfBO}1q5+jq&aYtU^glw22sk_0JlR%M# z4sVvP89I+ooLfJfc6h$|?~L8|y%`Hz9Pq@>6(2RW5hl}+TDRW-_G^x7mxp}sj!kbV z)Y&Nvp_Nae-bTH=Q%)7;z$Anbxb&7nJCsbF-3nV}3L&;aQx8hCrZM)Mw&hwUE5 zr4g#<-TWzx`{JtDq(gN)aX%0rRny+9(l0u@s0(l>HSPFJ=VUir?(ZvyxCZ4Kf%X0F zVUVFZ>hC&7rRu~+KdL&%*w#86qyBD#MQSCM^))9bj5=w|7oMaR(+6Me_~dnU{{wSq z;EjpR3$EHd=DbyWaD=b)d$?dg7(nflHvaeRp??K!Rw$tHJ#c_;d*Lx?T#-CTDMC2G zCTF)x6wWZG*5O+UXPAU@krS@4-NX5>gzub=(-S$Ir+z+ku9@ZuY4}j@*K7FB_8?Gh zd?&5X_ubqmRlQQ&fJR>fS(sDKZ^^%KLc{qDdn9uC?6&+@amlOk&ejIqppgltw9C2L zsJj$VQ2X99=hdcM8brH(H>k_h;Xl|`g;y-qs`!xzxgMeADCMf8?g_^zyP2zvuOTfY zXgbz6TSCd!8?G&|SXR`c3_ct^WeI`JFA!9KuOy;;EiTHG!;zl_(F{io^(}~|E;_8y zlcJwD7E!tLM~1MFmKb6|+o2@KnC2oc$Cwr*+ZMNg)+F0jSTn;R3tMMB(2i6&U8YjI z$+MjEI1hJ+R9n0z>JC!+<)d&<3&M#MA3Ub^u^cVFOy+^mOKuxG^Jh%nURZN()RCBvI#aknol~y)W08@Lk9wSn4=OkP1Ya;ECA5sWC0)gxOwvx(Z> z!Ksz!QoB0{WXyKB!ZslH>>$VvIQZ~v&j$^jl4zWRJr!nqi5k3*A845JXg%C9+jF&x z2MJMEUUcM1<%*-aGGDNgzhnWqU@bo_2~1KKiOKvtL%b1qcXf_=9gG8Q57a>6Mm-UG_j)nsG*>#)y2O&Ljy>!6B>jC zFC@qcTqJ0uUJ5}u*aohxIf-^M?wIWWg#^#i;`WkowTCh?vGU=b8DT|a$)DbR?d&Sj zT4HxALfZ+8s?CbQi3(`Zs~aN<^ww0h?I(%&ej+Ew?N|csjB*31;hhC?!eBtcHXw^$ zAV|)@>G&0XsNr;smUnj2d!!PT{3_Iu!0{L_J4cDCT-_dwAxK?%%h)HS>y3AlPEX>8 za9Br4np)lwYCF7m**nI9c3$=h#O7u1IK{IY?g%GIvzu@}JOxeI2Isib z#EA@k{Z)kGuxYARUSvlVcnO?P3A}8@-KZOO1UfpMaZ=lnEgRFOVqJw_QvJ#*`b(-C zk0s-7O5JnGsAZoy(=gSY}U~<;ai#5yryhu~S!iP1ghrWfp z^o$o8-$N+PO;8TDK{>N2iNaM5q@fVPQ{GNmh$kONTe7WhDJ{J)?|Y&MKCWp$gp+zt z%!77b7^FQ(x=e-lKmlZ8KSA4(Ns0-*I8&wd0w&7KdKc-1V8S&jirNlW8lRcb>}p^h z=tq>gcIpfAj1MJ40F3;CLcF7kXIsO;1g-CdI%&llPIyT0iZ=o!VH=INjU?G##)>x% z<&lsTZ@eU2N-ug|fGVviCB`!$rRAq(cq~-b9j@OM*9%#ZYDy_)nKsYbrYY!v0fjevy%FL;~1FhL{H(OaqUkrL8t+ffX* z5z?VMNiz~r}>G5d=yk`|KqHi5xOP8hKXKcXSeA+r2kq6C^q z3?FZ>_QAeW z#~ZAqL$5P@rI?|)$5By_>v%~(+0JOoGy-c16$5z{FB)K@#i8?aHB3!&Z^ ztWiXYD7Hsg-A?MjwwBSs8;wV~ya*7|FE4qWh!81$Md05})VS#6_Gi zn=wKw-XQ=)Jmp*tb zWW~QJkZ_fIQ#DkOkQ0PU5_v8F_QRN|p-93=vG|~8bitF}K2HdCtmVUkr~#)3e}1S> zqWV58US8ja^z>^}Lmz}M7CS3lo0|H7G#1-=fDfXFToSp{5+ca-ZRLYVB9kUivo@S* z=POef+xk0=Fdwu%*68o{qR8`oSVt*iwEgo=gbG*b??~VwA^rWFB+{M(E9yp9qwUXQ zI7L(2pO;`}#zGLbE_k3n>2bQu(zCyhD(PSiHc4-14}~Be@d4r`w)Hm~!-%&(su)4a z{YClGt5Ak*pxP2hy1#LPm!dm=XYX1Oq!hE6}Ae0jpffBYsS+tnMIRlrbmHbfu zq>ur<0Z#qmv`P|kaav97hjP5#O?k_QPb$~@$W1wYGe6XDgh2D5TDS$2sFtVDRv+wi z12KmA6V(oC{XoIJwhQH88@NONNuqF-GsJ#SAt4v11Co$)gvJRuM;!FQ#R(#OaXKdI hTFP-6uLE%%I74F;DspK$D=PA~=$s@H&jD5_{vVcCH|qcZ diff --git a/rockettests/16000m/16000_altitude.yaml b/rockettests/16000m/16000_altitude.yaml index 47760e1..027fd8d 100644 --- a/rockettests/16000m/16000_altitude.yaml +++ b/rockettests/16000m/16000_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 3741.6934729887184 + - 10362.151115706554 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1682.1653588139836 + - 1845.3774305531829 - Shocks Exit Pressure [Pa]: - - 3741.6934729887184 + - 10362.151115706554 - Shocks Exit Velocity [m/s]: - - 4609.12531520912 + - 4625.755965383679 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/16250m/16250_altitude.yaml b/rockettests/16250m/16250_altitude.yaml index 315f45a..ed0a5b0 100644 --- a/rockettests/16250m/16250_altitude.yaml +++ b/rockettests/16250m/16250_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 3418.3632146588047 + - 9963.315043672133 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1668.4012378582238 + - 1838.804521998011 - Shocks Exit Pressure [Pa]: - - 3418.3632146588056 + - 9963.315043672132 - Shocks Exit Velocity [m/s]: - - 4608.313106641779 + - 4624.754085521447 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/16500m/16500_altitude.yaml b/rockettests/16500m/16500_altitude.yaml index 5bd5850..120d4cd 100644 --- a/rockettests/16500m/16500_altitude.yaml +++ b/rockettests/16500m/16500_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 3122.972833474462 + - 9579.830051792753 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1654.7497402094962 + - 1832.2550249824837 - Shocks Exit Pressure [Pa]: - - 3122.9728334744623 + - 9579.830051792753 - Shocks Exit Velocity [m/s]: - - 4607.571083301473 + - 4623.790767713024 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/16750m/16750_altitude.yaml b/rockettests/16750m/16750_altitude.yaml index 93dadc2..60ad691 100644 --- a/rockettests/16750m/16750_altitude.yaml +++ b/rockettests/16750m/16750_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 2853.1079660570754 + - 9211.10528162192 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1641.209944340787 + - 1825.7288561188311 - Shocks Exit Pressure [Pa]: - - 2853.107966057075 + - 9211.10528162192 - Shocks Exit Velocity [m/s]: - - 4606.893180285303 + - 4622.864527716587 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/17000m/17000_altitude.yaml b/rockettests/17000m/17000_altitude.yaml index d4d9cb0..af694f3 100644 --- a/rockettests/17000m/17000_altitude.yaml +++ b/rockettests/17000m/17000_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 2606.5628809591462 + - 8856.572616676596 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1627.7809362653682 + - 1819.2259323162943 - Shocks Exit Pressure [Pa]: - - 2606.562880959146 + - 8856.572616676598 - Shocks Exit Velocity [m/s]: - - 4606.273856775692 + - 4621.973938418331 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/17250m/17250_altitude.yaml b/rockettests/17250m/17250_altitude.yaml index f8c268a..5bb9c80 100644 --- a/rockettests/17250m/17250_altitude.yaml +++ b/rockettests/17250m/17250_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 2381.322450192244 + - 8515.685807105872 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1614.4618094751027 + - 1812.7461707800717 - Shocks Exit Pressure [Pa]: - - 2381.3224501922446 + - 8515.685807105872 - Shocks Exit Velocity [m/s]: - - 4605.708050752694 + - 4621.117627633633 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/17500m/17500_altitude.yaml b/rockettests/17500m/17500_altitude.yaml index 2ecb80f..a84af72 100644 --- a/rockettests/17500m/17500_altitude.yaml +++ b/rockettests/17500m/17500_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 2175.5456786459445 + - 8187.9196280508995 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1601.2516648792484 + - 1806.2894890102593 - Shocks Exit Pressure [Pa]: - - 2175.5456786459445 + - 8187.9196280508995 - Shocks Exit Velocity [m/s]: - - 4605.191137619748 + - 4620.294275992843 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/17750m/17750_altitude.yaml b/rockettests/17750m/17750_altitude.yaml index 97e64b1..e29efe3 100644 --- a/rockettests/17750m/17750_altitude.yaml +++ b/rockettests/17750m/17750_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 1987.5506567760058 + - 7872.769070399271 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1588.1496107437683 + - 1799.8558048008053 - Shocks Exit Pressure [Pa]: - - 1987.5506567760056 + - 7872.769070399273 - Shocks Exit Velocity [m/s]: - - 4604.718892404689 + - 4619.502614908449 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/18000m/18000_altitude.yaml b/rockettests/18000m/18000_altitude.yaml index 8347c9a..3d28c44 100644 --- a/rockettests/18000m/18000_altitude.yaml +++ b/rockettests/18000m/18000_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 1815.8008135730927 + - 7569.748562687052 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1575.154762631136 + - 1793.4450362384607 - Shocks Exit Pressure [Pa]: - - 1815.8008135730931 + - 7569.748562687054 - Shocks Exit Velocity [m/s]: - - 4604.287455227081 + - 4618.74142462049 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/18250m/18250_altitude.yaml b/rockettests/18250m/18250_altitude.yaml index eb000e9..a00d58c 100644 --- a/rockettests/18250m/18250_altitude.yaml +++ b/rockettests/18250m/18250_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 1658.8923574511375 + - 7278.391222949545 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1562.2662433406292 + - 1787.0571017017364 - Shocks Exit Pressure [Pa]: - - 1658.8923574511375 + - 7278.391222949545 - Shocks Exit Velocity [m/s]: - - 4603.893299749618 + - 4618.009532317199 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/18500m/18500_altitude.yaml b/rockettests/18500m/18500_altitude.yaml index 3f4647b..59e6aa0 100644 --- a/rockettests/18500m/18500_altitude.yaml +++ b/rockettests/18500m/18500_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 1515.5428024038706 + - 6998.248139368094 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1549.4831828491194 + - 1780.6919198598653 - Shocks Exit Pressure [Pa]: - - 1515.5428024038706 + - 6998.248139368094 - Shocks Exit Velocity [m/s]: - - 4603.533204355715 + - 4617.305810327982 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/18750m/18750_altitude.yaml b/rockettests/18750m/18750_altitude.yaml index e7806f4..f59df62 100644 --- a/rockettests/18750m/18750_altitude.yaml +++ b/rockettests/18750m/18750_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 1384.5804856484378 + - 6728.88767860459 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1536.8047182523408 + - 1774.3494096717657 - Shocks Exit Pressure [Pa]: - - 1384.5804856484378 + - 6728.88767860459 - Shocks Exit Velocity [m/s]: - - 4603.204225817742 + - 4616.629174385942 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/19000m/19000_altitude.yaml b/rockettests/19000m/19000_altitude.yaml index a757466..b0e06e3 100644 --- a/rockettests/19000m/19000_altitude.yaml +++ b/rockettests/19000m/19000_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 1264.9349910789215 + - 6469.894820757962 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1524.229993706639 + - 1768.0294903850106 - Shocks Exit Pressure [Pa]: - - 1264.9349910789213 + - 6469.894820757963 - Shocks Exit Velocity [m/s]: - - 4602.903675240647 + - 4615.978581957292 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/19250m/19250_altitude.yaml b/rockettests/19250m/19250_altitude.yaml index 2f483c1..28af96b 100644 --- a/rockettests/19250m/19250_altitude.yaml +++ b/rockettests/19250m/19250_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 1155.6284002561815 + - 6220.87051991799 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1511.758160371201 + - 1761.7320815347982 - Shocks Exit Pressure [Pa]: - - 1155.6284002561815 + - 6220.87051991799 - Shocks Exit Velocity [m/s]: - - 4602.629096084364 + - 4615.353030635051 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/19500m/19500_altitude.yaml b/rockettests/19500m/19500_altitude.yaml index cbf132f..9c15f34 100644 --- a/rockettests/19500m/19500_altitude.yaml +++ b/rockettests/19500m/19500_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 1055.767299424274 + - 5981.431089331222 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1499.3883763507547 + - 1755.4571029429283 - Shocks Exit Pressure [Pa]: - - 1055.767299424274 + - 5981.43108933122 - Shocks Exit Velocity [m/s]: - - 4602.378244085363 + - 4614.751556594588 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/19750m/19750_altitude.yaml b/rockettests/19750m/19750_altitude.yaml index a711d53..edbaabc 100644 --- a/rockettests/19750m/19750_altitude.yaml +++ b/rockettests/19750m/19750_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 964.5354772230672 + - 5751.207610231657 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1487.119806638736 + - 1749.2044747167815 - Shocks Exit Pressure [Pa]: - - 964.5354772230672 + - 5751.207610231656 - Shocks Exit Velocity [m/s]: - - 4602.149068913238 + - 4614.173233108585 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/20000m/20000_altitude.yaml b/rockettests/20000m/20000_altitude.yaml index 43965b6..ed76f74 100644 --- a/rockettests/20000m/20000_altitude.yaml +++ b/rockettests/20000m/20000_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 881.1872534120462 + - 5529.845363425352 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1474.951623060926 + - 1742.974117248302 - Shocks Exit Pressure [Pa]: - - 881.1872534120463 + - 5529.845363425352 - Shocks Exit Velocity [m/s]: - - 4601.939697412403 + - 4613.61716911918 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/20250m/20250_altitude.yaml b/rockettests/20250m/20250_altitude.yaml index dffb559..d204a91 100644 --- a/rockettests/20250m/20250_altitude.yaml +++ b/rockettests/20250m/20250_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 805.0413840778685 + - 5317.003282753191 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1462.8830042195427 + - 1736.765951212983 - Shocks Exit Pressure [Pa]: - - 805.0413840778685 + - 5317.003282753191 - Shocks Exit Velocity [m/s]: - - 4601.748418291913 + - 4613.082507865063 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/20500m/20500_altitude.yaml b/rockettests/20500m/20500_altitude.yaml index 634a4bd..acb5c7b 100644 --- a/rockettests/20500m/20500_altitude.yaml +++ b/rockettests/20500m/20500_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 735.475493510073 + - 5112.353429589685 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1450.913135437796 + - 1730.5798975688583 - Shocks Exit Pressure [Pa]: - - 735.475493510073 + - 5112.353429589685 - Shocks Exit Velocity [m/s]: - - 4601.573668138291 + - 4612.568425561406 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/20750m/20750_altitude.yaml b/rockettests/20750m/20750_altitude.yaml index 5b7a37f..0257f89 100644 --- a/rockettests/20750m/20750_altitude.yaml +++ b/rockettests/20750m/20750_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 671.9209872340721 + - 4915.580487568156 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1439.0412087048935 + - 1724.415877555495 - Shocks Exit Pressure [Pa]: - - 671.9209872340721 + - 4915.580487568156 - Shocks Exit Velocity [m/s]: - - 4601.414018637013 + - 4612.074130130609 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/21000m/21000_altitude.yaml b/rockettests/21000m/21000_altitude.yaml index c425f92..97598d4 100644 --- a/rockettests/21000m/21000_altitude.yaml +++ b/rockettests/21000m/21000_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 613.8584046232762 + - 4726.381276753804 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1427.2664226214958 + - 1718.2738126929912 - Shocks Exit Pressure [Pa]: - - 613.8584046232762 + - 4726.381276753804 - Shocks Exit Velocity [m/s]: - - 4601.2681648982325 + - 4611.598859981903 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/21250m/21250_altitude.yaml b/rockettests/21250m/21250_altitude.yaml index c4e22b4..10439b6 100644 --- a/rockettests/21250m/21250_altitude.yaml +++ b/rockettests/21250m/21250_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 560.813173105074 + - 4544.464286516109 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1415.5879823456196 + - 1712.1536247809754 - Shocks Exit Pressure [Pa]: - - 560.813173105074 + - 4544.4642865161095 - Shocks Exit Velocity [m/s]: - - 4601.134914791303 + - 4611.1418828379165 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/21500m/21500_altitude.yaml b/rockettests/21500m/21500_altitude.yaml index 53c2e97..d20150e 100644 --- a/rockettests/21500m/21500_altitude.yaml +++ b/rockettests/21500m/21500_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 512.3517292578193 + - 4369.549226380816 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1404.0050995389834 + - 1706.0552358976145 - Shocks Exit Pressure [Pa]: - - 512.3517292578193 + - 4369.549226380816 - Shocks Exit Velocity [m/s]: - - 4601.013179200946 + - 4610.702494606414 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/21750m/21750_altitude.yaml b/rockettests/21750m/21750_altitude.yaml index 5d885ce..c04119d 100644 --- a/rockettests/21750m/21750_altitude.yaml +++ b/rockettests/21750m/21750_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 468.07797509473795 + - 4201.366594169518 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1392.516992313792 + - 1699.9785683986163 - Shocks Exit Pressure [Pa]: - - 468.077975094738 + - 4201.366594169518 - Shocks Exit Velocity [m/s]: - - 4600.9019631253905 + - 4610.280018295457 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/22000m/22000_altitude.yaml b/rockettests/22000m/22000_altitude.yaml index db191c2..466ad2a 100644 --- a/rockettests/22000m/22000_altitude.yaml +++ b/rockettests/22000m/22000_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 427.63004057031065 + - 4039.6572607613894 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1381.1228851799538 + - 1693.9235449162456 - Shocks Exit Pressure [Pa]: - - 427.6300405703107 + - 4039.6572607613894 - Shocks Exit Velocity [m/s]: - - 4600.800357543773 + - 4609.873802970324 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/22250m/22250_altitude.yaml b/rockettests/22250m/22250_altitude.yaml index 9fe99c0..1af699e 100644 --- a/rockettests/22250m/22250_altitude.yaml +++ b/rockettests/22250m/22250_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 390.6773258475868 + - 3884.1720708373327 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1369.8220089927354 + - 1687.8900883583372 - Shocks Exit Pressure [Pa]: - - 390.67732584758676 + - 3884.1720708373323 - Shocks Exit Velocity [m/s]: - - 4600.707531986279 + - 4609.483222750571 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/22500m/22500_altitude.yaml b/rockettests/22500m/22500_altitude.yaml index ac6f23e..e34137c 100644 --- a/rockettests/22500m/22500_altitude.yaml +++ b/rockettests/22500m/22500_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 356.91779915149897 + - 3734.671458991359 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1358.6136009008394 + - 1681.8781219073146 - Shocks Exit Pressure [Pa]: - - 356.91779915149897 + - 3734.671458991359 - Shocks Exit Velocity [m/s]: - - 4600.622727746324 + - 4609.107675845709 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/22750m/22750_altitude.yaml b/rockettests/22750m/22750_altitude.yaml index 79b4a26..e917e12 100644 --- a/rockettests/22750m/22750_altitude.yaml +++ b/rockettests/22750m/22750_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 326.0755281222745 + - 3590.925080617721 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1347.4969042949097 + - 1675.887569019211 - Shocks Exit Pressure [Pa]: - - 326.0755281222745 + - 3590.925080617722 - Shocks Exit Velocity [m/s]: - - 4600.5452516792975 + - 4608.74658362798 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/23000m/23000_altitude.yaml b/rockettests/23000m/23000_altitude.yaml index 181959a..3a390dd 100644 --- a/rockettests/23000m/23000_altitude.yaml +++ b/rockettests/23000m/23000_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 297.89842449154213 + - 3452.7114570050903 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1336.4711687564588 + - 1669.9183534226975 - Shocks Exit Pressure [Pa]: - - 297.89842449154213 + - 3452.71145700509 - Shocks Exit Velocity [m/s]: - - 4600.474470537162 + - 4608.399389740834 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/23250m/23250_altitude.yaml b/rockettests/23250m/23250_altitude.yaml index 581c6e2..a38eeaf 100644 --- a/rockettests/23250m/23250_altitude.yaml +++ b/rockettests/23250m/23250_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 272.1561836473213 + - 3319.8176340909617 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1325.5356500072082 + - 1663.9703991181082 - Shocks Exit Pressure [Pa]: - - 272.15618364732126 + - 3319.8176340909613 - Shocks Exit Velocity [m/s]: - - 4600.409805792624 + - 4608.065559241718 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/23500m/23500_altitude.yaml b/rockettests/23500m/23500_altitude.yaml index 3b2435b..e4804fe 100644 --- a/rockettests/23500m/23500_altitude.yaml +++ b/rockettests/23500m/23500_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 248.63840224699624 + - 3192.0388543504805 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1314.6896098588522 + - 1658.043630376476 - Shocks Exit Pressure [Pa]: - - 248.6384022469962 + - 3192.0388543504805 - Shocks Exit Velocity [m/s]: - - 4600.350728910563 + - 4607.744577777856 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/23750m/23750_altitude.yaml b/rockettests/23750m/23750_altitude.yaml index 7b16380..b7c7d94 100644 --- a/rockettests/23750m/23750_altitude.yaml +++ b/rockettests/23750m/23750_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 227.152858492648 + - 3069.1782413141887 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1303.9323161632217 + - 1652.1379717385662 - Shocks Exit Pressure [Pa]: - - 227.15285849264802 + - 3069.178241314188 - Shocks Exit Velocity [m/s]: - - 4600.2967570280625 + - 4607.435950793755 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/24000m/24000_altitude.yaml b/rockettests/24000m/24000_altitude.yaml index 8d4b6b9..6c8389b 100644 --- a/rockettests/24000m/24000_altitude.yaml +++ b/rockettests/24000m/24000_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 207.52394101263303 + - 2951.046496228574 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1293.263042762865 + - 1646.2533480139175 - Shocks Exit Pressure [Pa]: - - 207.52394101263303 + - 2951.0464962285737 - Shocks Exit Velocity [m/s]: - - 4600.247449007752 + - 4607.1392027692145 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/24250m/24250_altitude.yaml b/rockettests/24250m/24250_altitude.yaml index 39e8e31..4fe522a 100644 --- a/rockettests/24250m/24250_altitude.yaml +++ b/rockettests/24250m/24250_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 189.5912135079236 + - 2837.461606392069 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1282.6810694420308 + - 1640.3896842798827 - Shocks Exit Pressure [Pa]: - - 189.59121350792358 + - 2837.461606392069 - Shocks Exit Velocity [m/s]: - - 4600.202401832179 + - 4606.8538764866635 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/24500m/24500_altitude.yaml b/rockettests/24500m/24500_altitude.yaml index f3f35ad..3e7710a 100644 --- a/rockettests/24500m/24500_altitude.yaml +++ b/rockettests/24500m/24500_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 173.2081034313961 + - 2728.248564717109 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1272.1856818780452 + - 1634.5469058806768 - Shocks Exit Pressure [Pa]: - - 173.20810343139613 + - 2728.248564717109 - Shocks Exit Velocity [m/s]: - - 4600.161247309764 + - 4606.579532326695 - H2: - 0.23930995583841458 - H: diff --git a/rockettests/24750m/24750_altitude.yaml b/rockettests/24750m/24750_altitude.yaml index 2a130ba..12dbddd 100644 --- a/rockettests/24750m/24750_altitude.yaml +++ b/rockettests/24750m/24750_altitude.yaml @@ -5,7 +5,7 @@ Oxygen Mole Fraction: - 0.2095 Pressure [Pa]: - - 158.2407039820301 + - 2623.239100086162 Temperature [K]: - 216.67999999999998 - Combustion Chamber Mechanism: @@ -55,11 +55,11 @@ - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 1261.7761715930985 + - 1628.724938426424 - Shocks Exit Pressure [Pa]: - - 158.2407039820301 + - 2623.239100086162 - Shocks Exit Velocity [m/s]: - - 4600.123649065397 + - 4606.315747590723 - H2: - 0.23930995583841458 - H: diff --git a/shocks.py b/shocks.py index d22f0fb..ede19ea 100644 --- a/shocks.py +++ b/shocks.py @@ -16,11 +16,6 @@ import scipy.optimize def shock_calc(M1, P1, T1, P_a): - print("M1 = ", M1, "\n") - print("P1 = ", P1, "\n") - print("T1 = ", T1, "\n") - print("P_a = ", P_a, "\n") - gamma = 1.1 #1-2 expansion fan (underexpanded) @@ -37,11 +32,6 @@ def shock_calc(M1, P1, T1, P_a): P02 = (P2P02/P2)**-1 T02T01 = (P02/P01)**((gamma-1)/gamma) T2 = T1*T2T02*T02T01*T01T1 - - print("After exp fan") - print("M2 = ", M2, "\n") - print("P2 = ", P2, "\n") - print("T2 = ", T2, "\n") #2-3 shock 1 P3 = P_a @@ -60,22 +50,6 @@ def shock_calc(M1, P1, T1, P_a): M3 = M3n/math.sin(beta2rad-theta2rad) #Mach number after 1st shock T3T2 = P3P2*((2*gamma*M2n**2)/(gamma+1))# -(gamma-1))*((gamma-1)*M2n**2 +2))/((gamma+1)**2*M2n**2) #1.29 #T2/T1 from normal shock T3 = T2*T3T2 - - print("After shock") - print("M3 = ", M3, "\n") - print("P3 = ", P3, "\n") - print("T3 = ", T3, "\n") - #### 2-3 ### - #beta2deg = beta2rad*180/math.pi - #M2n = M2*math.sin(beta2rad) - #M3n = math.sqrt(((gamma-1)*M2n**2 +2)/(2*gamma*M2n**2 - (gamma-1))) #normal shock table - - #P3P2 = (2*gamma*M2n**2 -(gamma-1))/(gamma+1) #normal shock table - #T3T2 = ((2*gamma*M2n**2 -(gamma-1))*((gamma-1)*M2n**2 +2))/((gamma+1)**2*M2n**2) #1.2 #normal shock table - #M3 = abs(M3n/math.sin(beta2rad-theta1rad)) - # - #P3 = P3P2*P2 - #T3 = T3T2*T2 #Set new gas state gasPlume4 = ct.Solution('h2o2.yaml') @@ -83,21 +57,7 @@ def shock_calc(M1, P1, T1, P_a): state = ct.SolutionArray(gasPlume4, extra=['M4']) state.append(gasPlume4.state,M4=M3) - - #print("M2:",M2) - #print("P2 [Pa]:",P2) - #print("T2 [K]:",T2) - #print("Theta2 [deg]:",theta1deg, "\n") - # - #print("M3:",M3) - #print("P3 [Pa]:",P3) - #print("T3 [K]:",T3) - #print("Theta3 [deg]:",theta1deg, "\n") - # - #print("M4:",M4) - #print("P4 [Pa]:",P4) - #print("T4 [K]:",T4) - #print("Theta4 [deg]:",theta4deg) + return(state)