From d0b630cd3e4b7e4e037fb29e82b2e3725aa00bea Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?China=20Hagstr=C3=B6m?= Date: Mon, 18 Jul 2022 10:10:39 -0400 Subject: [PATCH] all fixed, scaling plots --- SSME.py | 26 +- UpstreamPlotting.ipynb | 14 +- nozzle.py | 9 +- plot_data.h5 | Bin 511456 -> 511456 bytes plumeSSMEJP.ipynb | 699 ++++--------------------- rockettests/16000m/16000_altitude.yaml | 42 +- rockettests/mult_alt/NOxMassAlt.png | Bin 0 -> 17303 bytes shocks.py | 8 +- 8 files changed, 149 insertions(+), 649 deletions(-) create mode 100644 rockettests/mult_alt/NOxMassAlt.png diff --git a/SSME.py b/SSME.py index 61e6e35..06d678b 100644 --- a/SSME.py +++ b/SSME.py @@ -97,6 +97,7 @@ h = altitudes[g] #altitude [m] ####################################################################### #P_atm = 101325*(1 - 2.25577*(10**-5) * h)**5.25588 + 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] @@ -112,6 +113,17 @@ else: print("ERROR: OUT OF ALTITUDE RANGE") + #calculate exit velocity + #Nozzle Geometry + Area_ratio = 19.8 + A_throat = 0.0599 #[m] + A_exit = A_throat*Area_ratio #[m] + + 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 X_O2 = 0.2095 @@ -177,11 +189,8 @@ comp_Noz1 = state.X[n] mdot_Noz = mdot_f+mdot_ox - #Nozzle Geometry - L_Noz = 0.01 #[m] - A_throat = 0.0599 #[m] - A_exit = A_throat*35 #[m] - + 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 @@ -190,7 +199,7 @@ results_T[1] = Noz_states.T[n] results_P[1] = Noz_states.P[n] results_X[1,:] = Noz_states.X[n] - results_u[1] = (mdot_f+mdot_ox)/(Noz_states.density[n]*A_exit) + results_u[1] = u_e #[m/s] from F = m_dot*u_e (rocket equation) #SAVE NEW STATES TO YAML #Load YAML file to append new data @@ -212,11 +221,12 @@ #Area function A_throat = 0.0599 + r_in = math.sqrt(A_throat)/3.1415 dAdx = np.zeros(len(Noz_states.x)-1) A = np.zeros(len(Noz_states.x)-1) t = 0 for t in range(n): - A[t] = 3.1415*(np.sqrt(Noz_states.x[t])+A_throat)**2 + A[t] = 3.1415*(np.sqrt(t)+r_in)**2 ############################################################################# ### SHOCKS/EXPANSION ### @@ -228,7 +238,7 @@ u = results_u[1] #velocity at exit of nozzle (m/s) gamma = 1.1 print(u) - P1 = 13789.5 #Noz_states.P[n] PLACEHOLDER, need design exit pressure or design alt + P1 = P_e #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]) P2 = P_atm #Pa diff --git a/UpstreamPlotting.ipynb b/UpstreamPlotting.ipynb index f2c249f..79ce4c4 100644 --- a/UpstreamPlotting.ipynb +++ b/UpstreamPlotting.ipynb @@ -68,6 +68,7 @@ } ], "source": [ + "#COMBUSTION CHAMBER\n", "#PLOT TEMPERATURE\n", "T_lit = np.full((len(state.T),1),3588.7)\n", "OxInj_T = np.full((len(state.T),1),145.9)\n", @@ -226,7 +227,18 @@ "id": "9969d84b", "metadata": {}, "outputs": [], - "source": [] + "source": [ + "#NOZZLE\n", + "#velocity vs isentropic\n", + "\n", + "#species vs isentropic\n", + "\n", + "#pressure vs isentropic\n", + "\n", + "#temp vs isentropic\n", + "\n", + "#" + ] } ], "metadata": { diff --git a/nozzle.py b/nozzle.py index f571da5..8d23cba 100644 --- a/nozzle.py +++ b/nozzle.py @@ -30,12 +30,13 @@ def __call__(self, t, y): nsp = 53 #number of species in mechanism if t == 0: - dAdx = 0.1888181/math.sqrt(t+0.001) + 3.1415 + dAdx = (3.1415*math.sqrt(t+0.001) + 0.244745)/math.sqrt(t+0.001) else: - dAdx = 0.1888181/math.sqrt(t) + 3.1415 + dAdx = (3.1415*math.sqrt(t) + 0.244745)/math.sqrt(t) mdot = 67.35 + 404.79 A_in = 0.0599 + r_in = math.sqrt(A_in)/3.1415 # State vector is [T, Y_1, Y_2, ... Y_K] self.gas.TDY = y[1], y[0], y[2:nsp+2] @@ -46,7 +47,7 @@ def __call__(self, t, y): #converging #create new function to find dAdx and A etc - A = 3.1415*(np.sqrt(t)+A_in)**2 + A = 3.1415*(np.sqrt(t)+r_in)**2 MW_mix = self.gas.mean_molecular_weight Ru = ct.gas_constant @@ -97,7 +98,7 @@ 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("\nMACH NUMBER AT NOZZLE INLET: ", M_CC) #print("\nTHROAT TEMPERATURE: ", T_throat) #print("\nTHROAT PRESSURE: ", P_throat) diff --git a/plot_data.h5 b/plot_data.h5 index 3045c5ed3836d3136eac4a8855a3ac163099e1cb..5a3af62ea727382a41de6667282095c85fb96e82 100644 GIT binary patch delta 23751 zcmcJXd0bD~|Hpm$v{<6-6k4p=8Oz90Ok>DUwm~6TQ;~fOg_#K1K6F)x?8`_#bjcEi zFpVjDC6mdnD9vC}>Gyh{bMMEe$Nm2KyMH;abMO6lKhOJp?tQ=SIkzVY=hx3VrxhucU3T1=J;SN_Pc9|# zn>;twf929-e&XpzTSEJEG!Z}EYO_GMr=y+zlB=gi|5vpm8tca+4%M;OoOJaoxHW5H z>&LEQr_q{N6Co(bq?Uf5wO&J3>dR|3<_){r)mcMlS36H?`ar(Y)@GA-OYf*nk)|PoE4EW01)za(ngMCWgpfkgN=L^Rn?QT%zBWdFWwUCD5Ef?&5s6Rkfx z0`5J4({o0_y=ds{5=(X8csua6=NQv9xuXmm7$L&)CchV0_Z2JTE z$Hf^(S{{>eSuj_f;z*VX9U7SD>ujfTv)R0z#t!Iwrj;>1AHNG)?XLi}3j`<{e7*)b zvbOBkC$qkSjtae=lF>s68e!R~!T#ttH9*h#r}sNumf(8Ff5phw5fwn)7t5An)-K`g z5T{aIaIga-hz%~If2B#g3|>72MxPFGdd&(1G6@RgldBwATL!lOm+w{sILl$myN&=Z+P~Lj;!FduT||UrB_{0#2Wsn+H{cdU7&(2ib>}v*P92-I zn=xz8;f&a=@GKg9#tM$CErU*o`U;x5zjwW;^zDi}ly$=RoNsExssS3@a%7b=r3tRT zwyZq-mWK=~&DcE-%+vzUeBGls-e!|_uY|C#jE+kacE7Ee_A*p?Ep42B@lg*G!;N>u zc++vab||q{v|e1fi_Z=78UOT}VxD-RIMeWT)OP5q_hFi-;Q% zSOYow>cx&Z3vbM|{n*Wve)t3`_-L**zgJ)j&O{WIo%~SnG zz5JZu8eTX&_P|-0P}sQHG(c!wMdL`8Qcf&5SLH(R$k*fxgQ{s97(`Gt4MLjWLf&bz zg@Ee3U)b2!BwBQ;p@Ea1fz_j?23|&sye!1k!Ecd}+$>nv(|q-FK;{6K%CX8(m^8vW zcvX!pYK+39XR|E#U=;M(d7{G(c@)Uc!jdKiTl0jqc3cSYdn`(YqxOIWhLC2jDPuf% z9EpC88a$8k0KKs}L8HZ>b6eqB5;C-BYa?izYnhP|GnCM!rwbEnhHO#;dhMC$dlLZN z;_QpCYs)J@Cw1nwiaDJ%2ya*|6r40G74o`(n3Ppl5EBd$AL>TLMKlnNOyb1Z&Kg8G zNvjv;r9#u*xCY9=y{V58_whwdbxV3Fxp(#-aA-iY4Qg>8H;37c=U@TkZlT|!(@zo*+R)`endo( zh`K2b!<4g5Gftzej+RRV`J7q^@+ul(_a-CW>^+tZS}#@dme$(Peth}QYVf{lXkRuR zya{U-FY56wOU_EzGpuS9s9>tjMa&m$qd^71GEkvS3|Ewb3hq1jc|y<*xSh@kI)y;5 zvP*+!cO6K|LdAxA$?=f<;Vt~P-{|T6e!|5>>%mI$gT4D`|7yQV4f2eDc?aA;{`Y*_ zK^JdVB=0C zV4Ec=mn=A%x+ymmGanJP6C~A>;+`A|a7czeEEhjVLPo#z)Cm24rKBHgUQwdoZHpgD zb1Pp1B;C%o`I|U0n>VJLV1TjuvRAdTAO{J4qZfa(N$5=NU zLv9);cB*0mXh`Pm+-mp*60)~FtDAhy$~ARo?HxExIS}tI?;ifZe~}sk@%}@){b>xu z2o~r`b-Aqq64ZPlwH_FWiYKY{C1usXq>KTvB())-_lAR!q+)ubCFQ5~ zI33D>y~)7{d-LT)=XuUb>`Tc`L*j=lPy;*VRo#GCVB0NyP`=Np0`_MoU=suKQ^Yti z@Z7X5XAqE^ma%T6rP1Q+b|#qIZa5#6Cz#PAM`*$v!Hz82i&Ws`h!|K6&&dvB{bw#aAK95qlaPP|yK7^s~W1)bhq zlvyEJSuzZGWH_79!;M7r2rhWxd2S?7ah{Mf67FWCl%tNK2pXToyMCB0djSNf*F`2E#5L6yq`L!R0pz`|j%LYF*3<0fP2&vP6Mhpeh zD8k8XmYh!K41%3Mh|$YlhKgqcD8NVK^qiT5j25k52;a=XuQ0(EFofN6jgSMKZd+ej zt3>{3^ZMUs)Cp1pdC!UQH^%^ZUcu`I4bRKS>{eEFDbQdmPfBfT)_@W;7?W?f!FpijY{%Sf@gRw@BVCZpc>$l zuDJa-2;ld=uiI6-rvkY9CWEagJ2y)LYZ^lE7NCMch}{Bjvj(d=x0>LQJr(^ye;`*G z7)9ieVtKTQQLVkbW7qGhcU4l38u;GsYTD9_HlJ=z77Jp-M7sTL+v@>+ZCK9Ipm zHmRfcb0H@0r|_n3DmqgS0+_%G4g#1ptdWI>AExjQ1Y!I!3U9QSmrYLKS4hY+dRBsw zm;J9q4m(>_$-laT@NM47X=?Ci_l-&a2K*no7_!2n0!+ml%T8w&S22Lj-Qx|3$QMY+A^R$hqq`vC)AVYtBk{Z<3brv4919fiH zoYxP!R-pFE<+k!_GvFzx2?*6D(o(d9_-8<7jdJs7ao`J@)GKIQ?#78@^2pgxqSagD zIsMn(C&|r+Dwig`cF*{siT7SL5T`9FZvGFJCZ}v{Hu`L|ysX3Lbo56slFjs>KJv*( zB>l&Qkd~z_6&@9X9_I-KkF=;aKXWCb-b942#`VcB<3!(bVkSj*@LH%mnUGx_s?~nY z!jBL7bhD1LavX#ozpcD}hbn8*CqDY-(3*5f?WS&nE38QiZ>mYQvU1Z>>p+G}v2I>&T4@pB8{l7*n+R~U81=Av_$3lD^`-ho)XltCRchB%NnJ3e=|6eT_o(6h zkc7+@!Juy1_T;}GLk!;Mm8Lz|U@Tt>>FL0g=;NS8gj0WrBbZr(RhrF+c^jDV4hhU? zG3G6;!AwG8j%;g$Y(He+#~GGN_%ETQsf~%_$@3E~B zE)PE_7xyX+?cA^7MpgD63^32lvjp$x$WP8|r^vi~f*%D|HyVJ;GfJ8{+*A%gn14)o*}+|EOzN<` zG0_Jx0ORJ|UN~!zJOC_Bai43j70Kv4Nw8#eKA4FDTjPY_1zd<&2Ei@-1dRZSOieI_ zHYvJ=pInB^NyyyxAx7Lk5BaCCl3vN({Nj<#Ekk#z!TrYKQ`JS_zWGzbmn-@hxY;Vm zY5;Lt^M#l-l3)uV;amtg;qWyhr6O~kHqca?NFjX;YklFH2pm#$*`hVd!6T0p2U+Y? zLbgBkINNWsDqB%L`MwWy2c%zvx6KY3wxaOq{Cf;M0xRK>JRxfb7h?Gi3U2CV*=Q{9 z7COX&mk5!FMvFIhQ*d7ZVZmMsZnT(_`yIsTBxEO_IArAH(Us=*?y_3Rf1%@qmOVzq zsNv)(i{g$n0skqJj2Ai03{IvsY?}mbtCWbIgk#BamWb-5Lr6rkg+ZtIAraXq!uu)> z!hN{;_8D?BBtu_$4%d*7(F@NTp+6tjD>-IS1)@dhh2AJdciHrOa3jmjl&Lx~>Gj3@Q;HHYJN^j9Kwescb<)#UAJcy{t_ z+l%{w*y`Si<)vQ?OPOqGFY7I+$Yz?EkWAmJ}KFG*!a{rPB z+-@qc7fOgxt z)NJuB11Ov0Z(QINveL@c6whk3bX03r>jy?+p-p!}o9`s$@g0b%o3DPjk+=V)iw97Q zqs2TuaS(o!gzV|`AB{XcecP8_>5r8I@cgLb&nKN#835}~trmX9hNKj~m$Te1%LBmD zQs?VkXAlMgP#=NbT@{uKhr4C!sh$2!} z7X3fC-Bhx7o?UD?YPl-=4;wZtZJ!PH4}z(k!|4j_OQsua#qy><7t+qV01y*mS^zvw zJZyDyhA9d!4GOP~v&i9SG3ea6xR!(r?K$5Fx@6+-BeG8^p@Yu44xN5emFq1FfO(!?r{1n&9%8i(>DY$)(3yZ9 zGS9S2uU`z%tkygSG(w0!N3l&-vqJL%?-p_{fS7tIv==4QWNw7BW3Lb}S`4`88h}a2 zlX=2*BQO8isb^foWF@fq{ljKqR9n zAN_n=CZLUWclZ&}?yf0TtsUWbw(4;+gF+ipdReynkrCzb1KDo3Je8CeGD4^PFleJ1 zv3w%(MNV%}-pN?;{^$$?C3{F$eFh+y%JZ+6 z<28l06Hq(klc&*Qo?cmi(@DslE-W+}1Mf@kx3%l79E0n>35(s@tFm!8vheOl+Bm#m zQgT8`#f`&W#oSgATkH_puv#iHOTbSmz68ju!Ai_BLhcxT?{ttEh|yrgSF~o1OvuC? z&CI^0-!Jyo`~J9-lDK^HrvTSR5o&lnUNkSX2JzZ=v3A}+8N7}PGeWWx#bpg{G+*&Z zwK8*{Lwv;})tW1@z}gJORyP_5dTlhct{I{ml3|-Pz%P)Hv1c_j!d?^h*2JTp65FP& zfA5q!>(#*C;+OY=@9Bn>^lN0e}UB8 z42gGdNLFWNwkC8)absIEqr|&7%w~m^R*Bhh>zrm&>#MRuJbPr7$uEH!oB5aL02djP zzcI7v1{hcgpX3RxyK^C8ArsL}-308_gLexdJpfFEyteFVhTwJsp-Uf%ZnT)2xAev7 zBqTRSerM$6x6v+h3JR24N0)?s9(=m#S~a|UBDDXlyTDJq*Yd*sZw&b5?g(yfChE!` zC7an=$6$CHN2Jnn2v;K3IuP6LG$>S64k!PTVpZixTu(v{ZP!so&`me}7yd3=3Ei_> zr#*vvsInpc_T8utmjGRL(Y;o4S{a~uRp|pF7%Nqk6C}l}N?%Ep`@+|B60a&Jk*|ls z*U3}K*Jv^BTK?cBA#;ZY7;!(af4J20cO`e5?LlF;My*l9*{x06+&>QPpC13=-?*{N z&8GC5xst7H6>}aqi3gP>vQV@HmjxhZ4Ypx+Ax-B|=%=r_#Mev7-H>AJa^vaz)W)aM zKv9Vp-mv(@!{w{hK-|`7$>g2nTGOXq4{OSZd^)#Y14dZMrn7T67t*I(Ye7sr?0r$d zI!W=*xz=+fZd^|>J|1WBqH+^aLo#Q_>F|3bWM^A%GjjHXD3jE^`;;T$@?myft>dbE z0O2ymd&hc=#P*iI-;OnxM}j$f$xefW79oz7V08+_1jCBT(?lGA6VlEQF=AO__o1p}L(#6vA&~#P4bXo9l_`l--d@TNL?vUp) zFrVViSs*4;@l6*V(hLF5Zm(?K&9}@3KL#ug39*mian@kfWe(lO7t_VLN{pb_MyA-d7^s9! zaB!w}HXO`vN*9PycLunhY&{&!shYV8*gC(^L|y%rrBXFsRHl=SQS%vlixX zcNpBAQ`H>qMvHHot%0+ckhjfRnSZT6`dVE{>eyXLek8lI)1Q-8s&P*nWOi+M7mz!B z-o7aKs!YyqjG1#CI9ii0`)4Jjr& znj7Wd>sMzcj%uzXd?^ao(|=Rt`=v@NXDw_A!UI-ABln+`3Hgk6ZwnS!i9p0WZf7L? z+jAjqY;TV6Cb0Mp;AIU~V0I!8>(I}(%GT^kz>qA!uXM-NB;)|^*u%)pwUcMximt63 zgp{FA%g;EeG6-2s=hv~tAY5v`q0!^xhC$#9nSBwH+*s6_{kY{)t@#~4B-NU3vSKTm z-_rm@!^QCfD8$iX>?S|r7f8tLvj!WnN2FXlU#wBG>vJpr<-Sdouelzz)wHkz`-)}z z-#F}(*?FxwlG`d`>nuW?tY&rQXx=S!7!8-BSE@6|5b0)}2{@gEOnk!E zi1>5%rSInbkMGw8UTakHqq!=?zCH^N()VlqzFGTZZpH7{Z2Y;c(nf^n9*X)wRVA-gv=Zo zYs7rR=i;ld^GfD!Kl=?iU9eOQ_q*M$v->HSNB14OSX?YK^MdoBWGgE;4}lYgOWS6T zNXq31s0q{Ms5vUHbHvX+(2S=jhh)^lPU4qHaOzvaHNE-QxBq)v*^iiUR7t)2d%L}j zi-Xmm&UDL2%mj6!NpxnfAembFK=VA<2!wA_xF9LMP2r-XSlQ{FEGxEu_7Xk3uh50; zOU-LE0+3=(it)pHa{X(KXGSY&gZ3r5&Zwpe?byWbPU)Z>_TuBsKPMY#*-LT%bhw7B z#HX5wYqa1bl$>{Xw~k&Y9lZ-q0>hV@_{&E)g;kvQ%{?m{mT>A{>|P~2{C{Biv{*c- zg(@GQx}S>sCk0p^-|E}V@RG6Uk2(o&p2Pf`pu(T_K;nTGU&^z61rq9ac}29;&o+z> zdc*sLpf@0-ZvM(YpJs77B;Q_EKrV%3LA|mF{|5MuD6OlzwH-( z^{+n-pWjN`%d9ONBx8}4)#R2-Sy?T9NGPafQAXzoc=?ZSnd(?z7C#0o+hW$x0v?AF zJ!~zEGVq#LLfqz!Qev*VUogAuz(O^=Yf((*X z{66$o4#K@-zoxx-qRKbN_E*CP?Z+T2t>I_w+sQBpe0!NQ97-VkzTHz&{Jz~wQvAN% zn=29P0EqECx*znV7)Oh7Z*c=R35h$hzY+Ft=Nqlt*iDIjDDCZ$hYqS>*T0ysb0@HO zcUxbtS4$b2Z!a4H7H(>u_S^mzy6Sc8m}V;hp8TMY>qF9~Gg6q~ckn+uC3*Z*A9GlDV=Q0+86pC*!Tg-*D!EA|i$mJ9M z@>4JuqDL?Vd8M^@7>1Cy(PG}7xDr20LiTq0Dx)#T9Csn{y_s?hrgnd2-qli-F}PH0 z+GY*L;M>mA4%V(Bj{!?Xo!1*|#s0DgE~KyPH-Hlu!p05oH*2uJ>9&bt+yc}dI-*94 zQO8B%ViGd7^>!m_%OC#O_x6MGo76SG)*li-b*>toPdW2;)e2C5%vw3Rae0ML6qoEa z*ous74;Lc-rW+z02tryMh>3-zWBV-ds;~wP1pAvAhbX`y8FlUvTuefy_B>`pJ!56| zGR+Gmwb?HXH26t`7c1vm}&5*7cO8Htj@A9l+Y z@@lkgMz+W70^vi6LCMAs?%eL_mhB!WQ71iU*ekWIDyUnZ_IWW6sCiHC tyc(7#qtb&LpEKP60a=O9neePei!TJ&rt*I2bEZEeh{{g^HP>uiq delta 23751 zcmcJXXIvE5*T=KFG%HvV3;si~1Z>y~sGue)SRmMmx|XOxRE)+3D%eF-kX1)P1+kZi zan-SeU_-?Qkr)dY43;PmQ4%#$a&-b2lXU>^>hldGH4-=e@ z?6NgHADwNgthL_s(8ppR zm!84q#daw=@Al1qi{B~U^lrJ{yx{Vxb>8OJ;+@coKHeU^!cQGt-N$>u?vBf9zgX?P z_{Z0-#_@l7zv&USZc^cu9`&aMZQM5Om}y5FC-YgK0G;`dY6o-dF`*+{;X!pYgWw zvkKN{{4U=7u3lqmuiTlhR&il3Ub(=n) z%t`NkUh=icY3)=~bT=2t(jrqrciv80S8h7fgSYc~RIua5@3xrd`7|*Y%4rlfB%6DF z3BQsczn1iIWo?;vY`+h@dB0zJezCQhclyb}eLMT@(UP~yrRZ1LyJTL$m}xcDbS=Pz zOf?q=Afm~BisI5h$^Korrn1581i@*QF5Y}<2;6HJlWoo$4)>yAayn6NYX6N3Z%gh? z8ubD9&X;34zPqO@zvRPd;8Q)F_S#LZ^& zb}~C)GES~we8%rm4TdTJb(#uLH28cCb!BbYuMcN@0KGle;jY_aHE2xrw)J<%Cu<2h zYjA^%?CBqY`puUu#jFkE?FjdBQ+T)wBS;Ldpns)FyAWP=0He_(*{oZPKqf(fd~~TR zYs~r@6>Xjyqm+z^iM(t((_O}*?1uann zi`E$)#K@?E_XXWl)4_)>WU7#9!kzfcpIqc!)RmiZ@|iGI=2x&YDoge8<}HuG+XlR8 ze}NYb##_74m9=Hw$f6ItvjS&St>XXgSD-g+0V#?P#ciqM- zI%mrSJ{VR*ciwlKim;B`!0D+MOj;5?Yg*$%pV=~@5VPJYL}*?`=Sr4RPAqs<>;NPRmg!jeaD0K9@&RM6$e zU!JHH5vwJ1>$9!?Sv~r)BDjUq?YOOCUOOGa8&(U&N3F|+!VVxNW#19R1VhC8IuWri z4MZb-viQBH4iQe${;5s5(6~FUfwFMN_E2zdeD)%H&G+#t?rnYg?d#WMotE5_H|VF| z{q-a6(mvoO%bl9*i-cO-&Zx`4c-7 z&2C?vWo8fh?QzSev$v=)Lsp23BBG+T#H=xFN|p0pU6C^mHnTmy0Vk{!jQ6cEnMZ+z z5n!GNjg}6D4g+-f7=B1d90P~b6;u?8UB}YI9tKR~1d4IASd1T?go{bYG4A}G!q=DD z^i1-ayIC~|EuOVJH^ePU%R%V&b+PX0NqG>cFchL^0F4k*W>Aci)u;=FGc$R&kT(;= z)Jq{g?*|I;(NHK2#h&xX%cSUTRuTqi5;FAI1q$f6|HM^`=Ljn3t6NrQ-Y8q4C3OF} z5tAP^&yYP&xu$S!8SuzhB%4Cca!DZD6iQdf7`t+XvX%UZh$0bn(@P9fUbJqs&M*Nj z@3rJ}Y9Yw0=!B$L1@FHGC ziv-7bP(iQ^RA`aF6{Vnp-&TH}5V{p^r*lG&A<(OA*WuY63`tq2*kC6)9+E%2h5vRd zp5D8f*Y5Dg%#AAYz1@2l{%E~aOY+TirqBJS%2nCZj6CQth{)84!xZ1d5b@5cNC-ay zY*I0g^yNarPm)Ue3Eb38*{C?uOcQ(tU44R_jTVdbTc_}&B;;84Iiv8sA?I@M(NILk~0$NQ4D!Qpciy{ERSP=;-;+>F?oR$mU zyx{$mpjNPMDlf$bFKK2k#aYNrlSTK}1P&?Yy-_~6-A;LIY8|m&1=;>;$nLtQ!?c7P z7Zl&V%RvhyP3nzP^sdb2m@FI2WS!oH6x@T8IDSAC(oWQ8a6b@QRC$ZerT93!r zY8*punk=@hq6ai2^LBnU`~nHt+W|H7AF}d2t1pdt8y%|}h~jfegMXhgPs@RL-@34E zVYTa$w^^W{{*>D)AVDn>u6+hZqT)$vJxSTu*HPhYD-Uia$VsE*@)&=R>| ztKXZm=RZRB>|wFBQ1J8SLSV^hh!3cNLCGmf&z@X~C_RypZv!6Qeu`cDQj9}MkL>z> z`VR$no9?G)B>fnzLfuenREcr>Of6CO%t)G@SpSB63uBM$ut4CE;cP+=R*3i^T=2s4 z%ut}>JRxr=+|5WSM;%TPz7wR0LF8?;n72zt<6083x5pY4-X5I#?_Ia@DAf?`+%)xc zg#LRihoI(&mg5giu?zvvPbUM77>XxTgp=7UISu9vg3}ZbqnEu56;FjwfFHo=dD96Q zEn2@2+RwzVFu@lvgrr#tdILdkNX_FTjdVG$!&@QUC^L}P>=DftmXMa6k;DaAJ_szqk5z!mI}1paWP0z4_< z$>g$(6)N!iL5n&+{&uRC-~$GJInK>+Q+W^fi?!H_vU9y8u%;n|#{m@-LSh`e%^Ixg z+@QxJTSR~8b>u3(TZkM|tc=$yYVEYK{fZYRe+Fk$?Sn*x}vxWn*@TmP1-oJn_MC|OS7 zUG-<_kVg~iT%lSjKO}mjQhZ;;XrRzc7T2F6H$yV?UKeo<2^sxJx&l4t`NTDukxNwQ zdfmz^`}R!K65Z-?jadD)kI)<5oDN_*FRE z<2E@QEk?cME>0&QQ@^~Ypl;UjkN#<$7OSW)1pR*3e0hSF)bql^bDP>^DyMZ&KDU)u zn<0-tO+ctNk(Qz*q&x;PYgF3Ai@r~3QvXWhlASD$EF@<`iS~ah=X9LMtbT7joW#9_IO+5VkSj*@Y;HVD-*J-BMpWRS$Os=pA9kp&zkh+0bT6YY}Is4 z+TpIvp)nt=NeczlBwJazS>515hD+sUO-XszG*qx|UT)eO5a18sUzD2&aI_e8;Ai+H z5;FDKdJ5{j%|zp0Bh}t79@+Sx!Y4bl^nTqFE#1Z@SiH|GO=qydSiThE<-(O{bTJ^p zsXxLM%&fsG%_hYB7|eKw1ZK1tb6j&UlaQDDz`=d_$2ro7hI!lVbMCxL|)Gok)uxc~7y8&?y zNtw8dPVn(j@CKO1Ek9$Y=GA@Ou2+K?P4^xomB$V3^+x996Z}9wFp{a6{S0_c69cVl z3fKGqNQxFVg~I-lvL674qnj^>3^X7Qw;PL#0}R$y^zu{y@hKs8r~>#=O-?H zK582%4**M3{AO8fMKbz>Bv>*!8_Yz3t#Ly594SF5?3{c>Ra7m?evU(r3`T_XA*}U9#S%EA=pAQJ-6M47bkqBNQL?Rk5W+YK?R|R45P6}?c zn3H?$#pxtuCm-3TaPq{VGlD0zSM#T~8T&=oAqiSK`Qp_#C20>VPNp^NhCgvzr9|{7 z97~q7MAR%DLL!|C_*t`x(=eA1O=q2zG6EQs-i8wWA1f#ItnW(FV@IkH z+dup{+y2|lT6(s=PyN2vr&ti#+p_zAgNkgXNnP-Q3tncH6@#7*QHVEuB^_dAreg^| zBqAS0h&RDl@Rg4@1P#f29rzB{kdS?Swp8Kk!)M>vWmTwmB4r$_>h;w^O?O^BE{hmY zE#5K^JRdb!xiAp6HOI?JQhfEpT2c|#Ru$|dzWQNng&1#6_wA}yh;g(ScU%oCS0*Ix zP4-qFiVwG~A^&~8Lyi3;*SO$x;wCNqAGxVU)iD_|HqS}xfsE`X_lMQzc8iJitq|ep zmJ-s2{3K~r(~+p#;2NqmU7W1oaZ>cqzU6{nAR&9)r?CQk{}1=WR`ph+zX)_5Q{|eb zd(Z5+E5(D`tzBoU`cKqb&_0>_dTiGF4O;rx`TT|*&D`Y8SNy?U;0>>kmG(YXcvhoj zqFS?BFEA1dZMqX$^p=!=ZxB;AU;X%!yxpBH?ng0>7W4G@{`gH2vZu2KC_O#t&vo~f zs|Vo8LDywRJvAMGgm*gM*x~X3u*a9@P>Zct{TRlDh}91SIS{0kj}hP}#W&rJB<$Wa z3be5KY7{vdlCdWn@e?Fuc875a_OMi^4fmR=+1r2ruKMsrCM}(B^G}ZfuUCA;9u{n| z70a7bxR7?}}^}ojR{ZD)ZeR?6cRm@vR0S~LCMa>8hld_KlF~P7!&Ba#8#si@8 zh9`@sR*9u5_RArjL#1!qe zbKJm{h)NUDP5o@Ao!2Jb&r8is6x_qem$dc~Pr#(u+Q&8w9XlkieQZ~_xv*_htKW3$ zLAV{|zf$;Ny_RmS);uEeYTb`E8v5d8r_#UJC7bX*8K9;OxZn1dJ8~E=dP433)P)y{z!^xq91c z4C|%_w!OVS`_HbL0$=*q_3zA?GBACMwU*h*HrizYk&LE%wDeb*fHvCQ;zvZMY%8o< z2gC7f)g$9Bg*K%0vMeb_LAj(feNn+(^`_a>yO9(7_K(prmhbwd_7HYhC|M$2^)Y~8 zD&J%G1VDsAUwH^Eo^r)&&A<40Lg-&`G@Zj2J)To&&j6}N5qTOd=IJHHIGu#->6b4R zV=yPd=91HF^%!Ja5axf`TGNfgGwbv(Zk;8M0WUM(aa%=fu|sIXYN^C513#&hG9a@C zD={kwd5$src&Igq(O|^q4c6FdrXX%(ZT%tr&Z$17mUS&P@h5R_Lwp*o)za%_{&Q-t ze=ieLg}GL;6UAjsZZu!2Kw*Tv=#%#A?W|K+vG~FTY z>u!wP79?Zx#Jgc9z`#oQq)=$ynF|pMnTT%cCg7m1yjzIq3Sc7SwWXIef;$R?-aRO~ z(PD0n>xt7zNN(QLTjAxrg}Jt&kJNj^!&j6JJl=Ss#)?sbhfH!NzP=62s4y5h$0Ra!cG?T^XlV)JEgcJD5jCE3bWF@FFj z@u0Fq7K)bOJqN_B!8Xj!rRlsC{q!}L_9~5^e^i5aNA4_I zsU>3W-XBd~_8)J9Z@vnQu#!z@&nPaWPq|iun0VOxqL64w@z1%|a3#j9p%_obS-hx> zC2C0K>|_&ukA&=Ohm8tncbI>p@3sK-NO-^gzOeRTO+SF>V2Hf6+Dje@=IpR-7F$tN zZs$TODiZ-qh?$9SIPtKea)&kY(EV_5%r1&>v>0^$UR+B;h7Q=TfR5Cy{M(~Q{S`z> zCh^|=q{~+J-seY}orpWJic3*8UeS}<^4rH>Nf5~;3TQ2RE zy&@_8C0B;5*p}H0n%>6=l0N6UNpTJ(W#L)BDu8!XE3Wrx6E$#Nhw1_Qf;IipWlNV{ z|2FPC1Q(N#gWUNq#Q+TOdhNITtonnOI}In>%{`&%2QMMPDMx;PCJzAro~syW+*Yh) zzT`qodr2{lZXxd#i0L@1F8`qzKM%#e=Xy(SCPmNel6Qb6Aw!QXQ$VlZ{$lf}(Q4>Z zuD0pLRWyaRue0++y*ly-AMCGv2-gfYF3ea~W?D&tm6_Hy5(br-HvEXRZPwNX?!E$d z=T)_VyV2s?W^3XsCgg3i_BJ2tk4-gpy6rDjf8BF1x1IaV2}`uRr|sz7zTKm4GC8}E zSFH~`LO`8~a3`yY+Or|=7K{yTDp)Vuq2OqPS$z{xz?L(sxY)qokYbvvjUoqcY<2G1 zyv=IDr=oB<>w>1=FLh`(d3nc4G9jPQel5WQD-np8$E_4%N^35}nASE3Zvsp4055B> z0<$f7_%{04R@qt|2^f+E`1#JbnuHwSZMrJl+{Z2I?cGl5K}i4lQRPW@O$TAksnUMc zVl0Ee7czSyCb_YwHG6T(rCPH$KP1(fU&@NDXnsWlkOLQ|^rH|*i?Qnm;1@{9>@x-` z*cXj?pZfc4_19eHe7isVHfs7c*P@^o1ICFmI}7pqL%FRYw$38N$!b<-j^N#b#|XF_ zz0yX7kwlyi;^I+6j20vAH5R9nkcp3sQxGqhm3Vp4Y&EfOdZV%dHkuOO4t}(`L*7Tk z4X1EhrHu+9@F-bM)P7T?L!nbaO_-t6Z1BeH8Tk3PnX@RyA*E;a%h?Ld2LXo%KKN9P z85=b|&Eef5E$?01WBX2BV)M~wTtd(yz>wwmHnT_wi2x^4?EB0}N%8M97fTA?XXXkC zOJv19(p*X&F9vbM3i3Ew%;UGB@arUGkNZR`h9We{_S5$#)SvCO*7eGaE!On2J*QCt zJ^Hn>38(oSzm?%xE#@+S|cQLS~Lk zR505gTN7L~Le1Q1*>{7E7cbD#{RNI0sY?gT%)H>-E7{5l&VAs7;nKF*1CsJS0BXYY zK4^o=tBm;B2bw7q<&caz@F;$X1gE|sq?>I%y#3cLc5Y}cQ-8ad^p(@j|Gf*>l)B#W ztakHdYJqQ4_!(>j!nY};N{Vk&_(f8z?EEH8R&4+5S$cTu$4a%QdpeB(q?mU}`S9)( z*?IP;L^W;bt`k1fs%c7_8+64qc)Nv`y%e941=o<3_*4^djTW4QlJgetHqi^EgW2FD zFnp<5!A@Zn=WUw+JIfMImB4pjmGo9)O*b~1-@lorAD{+KJ2mq0c^QlTsFP6e1m@G$ zR{UuXBpzt-r98)HAfbNmXGBZ=Y{O{(0^To#7J!ht`78e-n#DCC`S!A6aw#MW>Ln%k zKS;<^`Q<-~seIn4v(+@m)Y8elD9UN57ZI!l{Re+YL z`1Z2*l47^+=t{0c>@mZH{uE;sYOuu}3@wIz-p1CI2^qVMo$ZI@+x$ns9otwnw$;OP zLni#7>F2jyjCQ>Os>&Z_O54jEY+WQ{k(JfrmP=V#ZGK27u5DXE=LmTDk8hdk+F}+v z0+wwtYhVkHLy7*5wu%hAB>Qmm=FV!&=-b8LSL~UqrFVObyJ_#<%7V$am^A?z*-TTn zpeYxkeKT9cICWEE($I}}3*K(vMX&T>y}K@C-eH=*Wf1s%+Y=5Y5PskGk`%vhcas#qZ+GWP#5w?C z+>P!BJt@Y~V%%|Gf}4cIy{V4^+pF#Hz`J$S*n2boK5)lHQ|$7jFhis7Wo*8^Y!Fzu zsfC6veQZrN>N+vaRssU}K_M*wKF3+o{<1*wd23)t3?*!|81}6Zz$PKP+b2lD{def> zZS#ZIsXpIY8M>`~bnt8~WB#6Dc6o<5nOoXhHXdv+R$2v^z?GOc!4@}s>aYA3%&fuw zY{n#-*B$75wy!LN;v7=S!Xu|ENM|Jt_cWeUe_6ix&$S_GqNZP#=d65wQ#@)RWwX0r z4%|YPV{;k)T*h3GLh%h|^SO{Vn1xA)yx-z4KZSE4`iE1HdjZuxg1n6u^Y-{9_*oLN zx3iWi#vuJcVz&*B>M@wq`I$|}>Y9$h?q!a{cHNT4fTg0IYb>^6f7x0tq_6AOffE=) z%sTj+HQ3*Di=`Ov32J{6QKQAElQ-dF5;C>JW(D<;+tUxZ&Re63@l`+88mX1#Fhr*Zq5Grk2!BD;8cj?eNiu=BJNwTcuTi<8Uug@xPgI zLQ?p{Zn;9C$Q54>IBAPLW`n@U_L!v-K9m@mrhIU>za5`Fy^9)k!tV{bU2CZ+YDWDR tt)H}(QR&qgK4-cD0= x_max && u == 0#desired xlim\n", + " x_lim = g\n", + " u = 1\n", + " end\n", + " end\n", + "\n", " heatmap((x[2:x_lim]), y, clim = (0,Inf),\n", " var[1:y_len,2:x_lim], colorbar_title=clabel, size = (700, 500), dpi=300, c=colormap)\n", " xlabel!(xlabel)\n", @@ -463,7 +472,7 @@ " A[1, 2] = 1/(Δψ[1])\n", "\n", " # Dirichlet conditions at y -> infinity #change\n", - " A[n,n] = 1.0 \n", + " A[n,n] = 1.0\n", "\n", " # Definition for convenience\n", " Pr = ambient.Pr\n", @@ -528,608 +537,15 @@ }, { "cell_type": "code", - "execution_count": 14, + "execution_count": null, "id": "e9323b4e", "metadata": {}, - "outputs": [ - { - "name": "stdout", - "output_type": "stream", - "text": [ - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in cp/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 16.406463656499994\n", - "\tValue computed using high-temperature polynomial: 4200978274.3966804\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in h/RT detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 11.881653851683332\n", - "\tValue computed using high-temperature polynomial: 1050244576.379208\n", - "\n", - "CanteraWarning: NasaPoly2::validate: \n", - "For species N2O5, discontinuity in s/R detected at Tmid = 1000.0\n", - "\tValue computed using low-temperature polynomial: 59.91566118948061\n", - "\tValue computed using high-temperature polynomial: 1400326145.9123998\n", - "\n", - "(80, 58)" - ] - } - ], + "outputs": [], "source": [ "n_species = 58\n", "upper = 40000 #[m]\n", "lower = 16000 #[m]\n", - "space = convert(Int, (upper-lower)/250+1)\n", + "space = convert(Int, (upper-lower)/1000+1) #250\n", "h = Int.(LinRange(16000, 40000, space))\n", "g = 1\n", "\n", @@ -1250,7 +666,7 @@ " \n", " # println(\"started solver :)\")\n", " \n", - " for i=1:n-1 #x \n", + " for i=1:n-1 #x\n", " for j=1:n_species #species\n", " #calculate f0 at half step 0.5*Δϕ (x)\n", " χ_h0[:,j] = solve_exhaust_flow_χ(u[:,i], T[:,i], ambient, n, 0.5*Δϕ[i], Δψ, χ[:,i,j], i, j)\n", @@ -1281,17 +697,72 @@ " χ_save[m,:,:,:] = χ\n", " \n", " #println(\"Splitting complete\")\n", - " \n", + " print(m,\" of \", length(h), \" altitudes done!\\n\")\n", "end\n", "print(\"done computing! 💕\")" ] }, { "cell_type": "code", - "execution_count": null, + "execution_count": 37, "id": "d011e69b", "metadata": {}, - "outputs": [], + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n", + "┌ Warning: Duplicated knots were deduplicated. Use Interpolations.deduplicate_knots!(knots) explicitly to avoid this warning.\n", + "│ k1 = [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.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, 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.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0]\n", + "└ @ Interpolations /home/chinahg/.julia/packages/Interpolations/Glp9h/src/gridded/gridded.jl:77\n" + ] + } + ], "source": [ "#PLOTTING ONLY IN THIS NOTEBOOK\n", "h = Int.(LinRange(16000, 40000, space)) #make an array and put in loop for all alts\n", @@ -1316,13 +787,14 @@ " xx, yy, u_g, T_g, χ_gNO = regrid_solution(x, y, u, T, χ[:,:,36], 0.01);\n", "\n", " #PLOT 2D MAPS AND SAVE\n", - " plot_H2O = plot_heatmap(xx, yy, χ_gH2O, \"Axial distance, m\", \"Radial distance, m\",\"Concentration\", :magma)\n", + " x_max = 20 #[m]\n", + " plot_H2O = plot_heatmap(xx, yy, χ_gH2O, \"Axial distance, m\", \"Radial distance, m\",\"Concentration\", :magma, x_max)\n", " savefig(\"/home/chinahg/GCresearch/rocketemissions/rockettests/\"*h_string*\"m/X_H2O.png\")\n", - " plot_NO2 = plot_heatmap(xx, yy, χ_gNO2, \"Axial distance, m\", \"Radial distance, m\",\"Concentration\", :magma)\n", + " plot_NO2 = plot_heatmap(xx, yy, χ_gNO2, \"Axial distance, m\", \"Radial distance, m\",\"Concentration\", :magma, x_max)\n", " savefig(\"/home/chinahg/GCresearch/rocketemissions/rockettests/\"*h_string*\"m/X_NO2.png\")\n", - " plot_N2O = plot_heatmap(xx, yy, χ_gN2O, \"Axial distance, m\", \"Radial distance, m\",\"Concentration\", :magma)\n", + " plot_N2O = plot_heatmap(xx, yy, χ_gN2O, \"Axial distance, m\", \"Radial distance, m\",\"Concentration\", :magma, x_max)\n", " savefig(\"/home/chinahg/GCresearch/rocketemissions/rockettests/\"*h_string*\"m/X_N2O.png\")\n", - " plot_NO = plot_heatmap(xx, yy, χ_gNO, \"Axial distance, m\", \"Radial distance, m\",\"Concentration\", :magma)\n", + " plot_NO = plot_heatmap(xx, yy, χ_gNO, \"Axial distance, m\", \"Radial distance, m\",\"Concentration\", :magma, x_max)\n", " savefig(\"/home/chinahg/GCresearch/rocketemissions/rockettests/\"*h_string*\"m/X_NO.png\")\n", "\n", " #CALCULATE NO EI FOR ALTITUDE\n", @@ -1367,14 +839,18 @@ "\n", " #ppm to kg\n", " #Volume for each \"ring\" so we can turn ppm to kg :)\n", - " x_lim = 90\n", + " x_lim = n\n", + " y_lim = s\n", " volume = zeros(s,x_lim)\n", " NOx_mass = zeros(length(h),x_lim)\n", - " for i = 2:x_lim, j = 2:80\n", - " volume[j,i] = 0.27*3.1415*(xx[i]-xx[i-1])*(yy[j]^2 - yy[j-1]^2)\n", + "\n", + " display(xx[n]-xx[n-1])\n", + "\n", + " for i = 2:x_lim, j = 2:y_lim\n", + " volume[j,i] = 0.25*3.1415*(xx[i]-xx[i-1])*(yy[j]^2 - yy[j-1]^2)\n", " NOx_mass[m,i] += (χ_gNO[j,i]+χ_gNO2[j,i]+χ_gN2O[j,i])*volume[j,i]\n", " end\n", - "\n", + " #display(NOx_mass)\n", " plotNOxMass = plot(xx[1:x_lim], NOx_mass[m,:], ylabel = \"Mass NOx\", xlabel = \"Distance from nozzle exit [m]\", label = \"NOx [kg]\", lw = 3)\n", " savefig(\"/home/chinahg/GCresearch/rocketemissions/rockettests/\"*h_string*\"m/NOxMass.png\")\n", "\n", @@ -1390,6 +866,7 @@ " \n", "end\n", "\n", + "display(NOx_mass_alt)\n", "plotNOxMassAlt = plot(NOx_mass_alt, h, ylabel = \"Altitude [m]\", xlabel = \"Mass NOx [kg]\", label = \"NOx Mass [kg]\", lw = 3)\n", "savefig(\"/home/chinahg/GCresearch/rocketemissions/rockettests/mult_alt/NOxMassAlt.png\")\n", "\n", diff --git a/rockettests/16000m/16000_altitude.yaml b/rockettests/16000m/16000_altitude.yaml index d5db250..4f2e63f 100644 --- a/rockettests/16000m/16000_altitude.yaml +++ b/rockettests/16000m/16000_altitude.yaml @@ -33,46 +33,46 @@ - Nozzle Mechanism: - gri30 - Nozzle Exit Temperature [K]: - - 2835.0991405223144 + - 1807.9617523502745 - Nozzle Exit Pressure [Pa]: - - 1616724.2886052236 + - 89174.22494486334 - H2: - - 0.23997922000883562 + - 0.23930995999292012 - H: - - 0.014909323208890678 + - 0.0063865509611238745 - O: - - 0.0005116983515446788 + - 2.7962317493488678e-05 - O2: - - 0.0006879208258059928 + - 4.725534126669615e-05 - OH: - - 0.013149399432864632 + - 0.0013516335073095248 - H2O: - - 0.7307582552725053 + - 0.7528764318026328 - HO2: - - 2.9706399010189216e-06 + - 1.2348342691158568e-07 - H2O2: - - 1.2122596521394767e-06 + - 8.259382660355577e-08 - Shocks Mechanism: - NONE - Shocks Exit Temperature [K]: - - 2443.0010319568387 + - 1682.16524035398 - Shocks Exit Pressure [Pa]: - - 3741.6934729887184 + - 3741.693472988718 - Shocks Exit Velocity [m/s]: - - 235.93049978732773 + - 4609.12531520912 - H2: - - 0.23997922000883562 + - 0.23930995999292012 - H: - - 0.014909323208890678 + - 0.0063865509611238745 - O: - - 0.0005116983515446788 + - 2.7962317493488678e-05 - O2: - - 0.0006879208258059928 + - 4.725534126669615e-05 - OH: - - 0.013149399432864632 + - 0.0013516335073095248 - H2O: - - 0.7307582552725053 + - 0.7528764318026328 - HO2: - - 2.9706399010189216e-06 + - 1.2348342691158568e-07 - H2O2: - - 1.2122596521394767e-06 + - 8.259382660355577e-08 diff --git a/rockettests/mult_alt/NOxMassAlt.png b/rockettests/mult_alt/NOxMassAlt.png new file mode 100644 index 0000000000000000000000000000000000000000..1de756b015f8b5da1a8d6017710d8b2564751193 GIT binary patch literal 17303 zcmdsf2{hGh->%Z2fp$^KtP+)>BJ(UHnWxN^L}oH%jy94EX+X(5j}3+p3Q-XmQe;et zZOA-jIM?oZ-uFGkVJ*dAEt6Yre;2eB|?J48KLZM8`FLF_kl1bi44uMUi&nrbGd zI7cl*MM04h$hLukLRXG}A7ryrt)sa9i9n>Fcx6Q;OF?n|*9viVc_<*0SoTtv$rcJ3i50oiwm*gTKht!O^agy9`ITxNi2Cwlz1uKAcQJ zA$?(DD!An8&!ykLE8V9L2?-UCci!bx3){!9U|{f&@NIZl#CvImno+Q;%q>49<k45hjoaRP5u;B3`SJUJ}_Jp%axjH&7j^tN%7Tc>)9Y1!Ajd+lQ zBRC8jyHZtE#V!{}knJqE(weLoj&)pKo^B!%Xld{6M!d!2>YSsfH73d?-oAa7R+aC7%sCB>$>}C3RuLy+XMeVy zi%a3{+uIAd4;;|a(yIPBot>36G_;@i_Wk>opUtP0lx)8i*|ue9;3>YAx&5m1aZEjV zOglLxB?a;N>C>k-qdO>igL=w6i)}hIXX3+}4$*B*O-ZTnVj|4W47YmDkL#a1CouNp z_J;n|mEVhAt9a`cc09vTK0eb*AMaeP%g)Zu_@*My&!0c7t*x(?7#SOP`}?miRrB$w z7~1T!G&4N*wfJ#tY*lr2cXzj8fn~EawM)>m=g+Mp)izV~9+iw|-M0@bQ$NVNf%G|3 zi}TR<<;J(~P6!GL_Vo0S{=BB9rjTn{Sy_ckU%q}-R8-We$2~kRF6Jbj)7Rgt}@3CVOZJAnWX=zfc3kPUf_fSw&T)ExrZ0P)1#njaI zQeKxJ#m?qnnxU^))T`_|(~HF!OMJ8$Itw&gRp~ z%D<4At1FM|e`IB6f9>edkn(A%TSI0SV^4>MhIoxjiGhJ9R+QFJ^cfyRI6o>2%vYTs z@4j*4#v75ea|Q;trpKdi{E5pv-CP3>4vzl4{d(D2>8e3A%vQY>-UwaVxRdXbJVQ|< z+XDX7e6#Kw?&~NjE+evj{P^+y{rj08w}16k>gnpHsYPaFWhL1i>qwN_tgo+sJSZ32 znUSHct(_v2nVET*hv$A&l>LHcf!Z^~EPi#oy9|%Wj|Dz*#Ohm3@Hd-leXsnwx{S7M z+eT``qFj)GfPik!Iho?2;b9Fu|DCd^9%kpyA9KVn(Q=YTMuLKZLPSOp zq$dRgd{$=lN;!mu&&kP|jdd2PoLfdU+^~M#{hY&x55Io>`l20fb7`pg0Eu+<9l3>h z8PCklqI$+gMd@X)w4^GxrmIPL%naQ?aS*@e@cr!_aLrO)z#ItwQ2ol@RMueahCN7GBM~;2U0P_@Ey}?WKsaIFBaZB4k?_D$$}1?`S`60BGa4q542+EMSQ3ukPRYqt(=dy2 zn^wHk$<~cODW|f%4jo>Sq>%Tdbl84)mQ}(~sW4-jZ*uBKy_oY5fj~HKV{=wW(tV0p zDzYUP74G-k*a70%m=hfP_g`#Il8=jv3!&$$@L3V{=C{y5#$UK_!OqSOf8y)O$w}?9 z$jBDpj_upG6J`7H)TPeDEsw?A6i%H2bhu<;p{%S-Pfvffr0d?jd*z<E?Y0(!`&UdHsAVl9v)F!y0fDL_1gYdPkDV5UvWW! zxZ^isSXd)sQBGENY5rG!LIQ$*=h36tn9#F8fwr(Bc<>lQ>!}RvuI}I(Z)fgZ`9@HUKk}{$YzuUXz!? zCr-3)yI$i^#PJHug(zBiMa>e6VrE2-J0ENv{_wj|7Pk$OCd( z=j5DKQVM1%0+w{qvUqiUot2dpk>Km=OM3kjh!^3T9CydU%8G%RnR(y7%$%H;R-du} zr@sIAp`)X-viLKeiN@l~lMAb>tJGcu+y%+X0|fx>10~x0^5r7?Udcm;2C;i>Pivo= z)d_YtcuRIUm$AGUp3__3TF-UX-WXrfYwGD07ZuIT&BYzNAY0D7;N^oDjfvS8UYG22 zZ|Rc&y5IZzM^TQ&oW7r?9N>3uwH@8LonqEz$WmdC~vPw}*zx2-|iRpr#7C?WUlZ(CHI&TNmlNyf8U4 zu`2*%P0!ERlL!=xw5?p8QvQnwZFBvLK$%uBmE8c+H*03;M_~vF2>}=Yt}!t( zS|)wU9?Zk-KV;SNBqAoueO5d)H0~rO_9H2{d3dZ{ULdXT;X_e?`*t1AIM`FNvaAMb zLy@QmtFh(9X>?9T@BE~sBorkHuLX5|efP##5ob>wolj{hkG>5Igx$`t=Fy`^ckb-mx$`*c#i6Aac&OZqR^i7`fpAetNy!r@PJm^-ZD`2()Q^m8Zf>q! z1Lcq_L<$a|r&m-|$ji&G{GO-Z8W39X89f3WpNW-q4cIEAX!C{*MlT#LfI9#M<>lp- zd(0vM@rn_$H>nF_zT=Ozpt`e*^GZcTYU+f04 z4uogZ&ykUlHGp_zJ_qqOvzRO5u)EAH-6PMS@Qk9O1wLt}h*K@j=*julSfq@ofPjFd zr6mv{!Q;mfbRe<6=YOpYr^i#V174{qD?cwPT0QsLxXhJY0+3m=&wx=;`ayPhcyx4i zS8En29RCB~dG<_GMn)#Kv85&S*RSKA&(30T(FXvo%3MZt-Nrzs(o`R3B`5RSbqUkb z(q^Ti;VUR85D7FvX7={1tgLDq`+-?e$#Gv;)OlR?b*<{Y=vlCKM@dmuy+%K_tsOwwU(ZIi!dkVRD(0H?4L7zr%$-Hs1y}iA& zv$HOom6Ldolamu@kXiH|sR;0cax={-1%)ssk*T+L4vdeFQ(IgB@dVWfW05#RRzry{`O`}bMUPtcZ-W`Gp3iN_=*;}F~HAQ2TCsR9w2D($OL1E;bJI4hD z0Y?Qs`3$#kcw|5|xM^C6hc_~c6_yx&=C|daY7k9)6CWR+ebu^{QZIli(&F?_pfxb6 znb#$NK~+_Wn`;B-PZZ#AqqU+$w|o zuHoTeQ?GyUBK=Jkr}_U1Qpi*diz)S-HLu&=k>di8G;igy5$}C}a)z|#8JBq!qlu0+ zCkHM3dh^2xl<3hVE;J*?J>Axzoe8}EK?wW5S@8dDRqp@rW*mo8zkT~QI%;!Vi2wNU zJNNIKryWF~zFQQcQi~Pt9~d|=c)2NFN=Qh^;;$d{->R_etyn?pVC=5#?tc8#18+Y_ zl9IuO3=a;r8>Qm!r+X@8gs)BXAsIuSe*i{>@Urjd7L;D6RccwDKQG6k;r3icvsv=n z>3HMuKiz@$uLfFL&ykBDZWb06(~DlsI;h>d4sXY=?jr4|0%{K;5}n7pI%;Z&^|_Th z(Ei*K;v%wha=K8~ZtXt0htGt25k;QjU}l^d5zJPc#U5NVGJg{;QAOG!i<*N%)OU^z z(x%Kx-5RtGK)w$@oS=JvmYmYd%+6L$Ifg4JWk^xExVVIdhH|Pr)ONeC&UNXY>iH5%%gTByeSP2GI}F+eF=$s_utsWXs;a6giZDH|k>$7nf@tv* zExAh{KYknyVSZ-jo|nwoM^e5OQBj{~M?OC)gDAUjr0wc6NRGF0Bj_XFvoCbwN4ON= zcO5hrNDX^tHzeM=^&T`7@*BsZ)m|M9ji;C1MczNEC+G%fHL*K$ZEru&dAeGB+~db6 zoJJ$6U14Ei7iq|4$iN^gF_E2vWA>>e)&@lY1_S{r>K5Ixsn; z)RE8mkPN;wH*3ht2M{1cUA%Zv*52gYIW7(kdE+)NK0Ybe(Q}H5btombDgyHJ^W)=t z|6&#*AeY!d6cx>Ee7wUVc8wsLLToUlM}gf#M<;OUEup#Ak8<_(^9$D2Vpj(nw$L6T$b!#<=2bclD9E0^a3L`{Ik}Sjl9K3qcbAi% z``iSy7m~S?)3d_DSty5eyLWG4cXf4z{xdr_wgtXJQhp@ddTsPj*eTlZ#RO}9TKu_H<9|KCpZK;|W(ygvTtG-@?F+pD7u~nAMy7vG6cpJA4kcZOUe>Pk#=dkO);xbcLHMh% zhzN*`U?-21)EdY>lqhJkC_9;jMJXvn)x#Vdskd%zWeL~Q)ZD++ffVYg^i|_Q>>0hh zHVGO6j0JK#(w=L*&^x030cDdB_{274pW=uGtN&Z{YOL%67@^U%j3!S#O9%R&(k9@s(j1-Q$d!zO_PU%PgVAiHbV zE>RcKDXJ^T0$ds5Y*CTmCsk&)gT$AM5_#B=}l8Xl(?Eo8A~gxGjBALZES6=jTiNl@oC8Jcu8+i@ucMBrpIE6-7VnS z7C;3kPN0#bwoJLe?bt;P4Gn!z*v>-hEo8i-{g%hHVmaP-b?4TAnb-0=cI<#yyZfjv zLAEhYVuVC01t-8yV+fi?M&v$>x(cDt6oQ+vF&XwD(ILgrAT$cCKHWK>cul+q(yyf3 z15AR-X6*@cSZqK+l@pj1jY8502g2#E-hB1yG^A4X&phCZ1vVX#r4!|Sd zB#o`E{YXJc&?}!fM7T+yCOr@99Vsqj~_n# z?&HUZik%Pf)F3#HCAhAa8oC`M8I5EG38+hCN@a#p1zZa?oaO&a5)(6oBGKNitE8k9 zd7r!p(sJQLckkXE=_xO(IfHLRy!5?sOib)I+I*;8-_X$8OA@%uE`)|jsq?$Ky7Xph z2aA2^lX1(8M1-i)>C*`haGUASD6vXzZf-!icqj;tqv*)!JRh_HGGLLLr7522hA5PmsR{)L`)1Y z$ybb9ka1l`05TQ)tmV&d$HsO+&H<@C`pya;aBzF&RkFJAsZ!DM!UYU3On+qoN*@{_ zFKour4F}v1z{^ZM#xswE>><7`E-jtO`UldX=P}UEmS-a@EG=noy4UWDV&0FU^%YDsOLZ zvP_XZ35I__G3-*sA6%L#MFl<7P+eU;`)NW#AD+J<3irw;67R4GP867ycIfI~2-G)h zXPlI;q_}u2xg4PbTx|Ywl!piOQWs^F!-y64k-SPveDmfFS_qLK86SWDzK)U6*?dhc ztw+L+-ym*q7_;GBck!+nN-98z)ujiU^1RAP$S6-FY2f|6%VG|3J z17fbDY(#*FjT<&#z|xYvb=x*fmBqal9F*umH(sS6K=bP{&F$`9`5Etfxqtyrx8=7Q z83*GX^_ma`-$K_S(s;;gs)bo?2em8DV z-km)*4glo4G8+{ZMniTU+_WShofj4s1WI2Ed507NQIh8Z=bl~!cm87*@Y;=wfHT+C z)d`s$-u1-l(-BA>E-tfwBW5fXs`!Hky4B_G;w;s($A*ws5Uq{TK_E?E{fvh<-R8&_ ze;OeUQDxh<`VSv=r~9H^!8;%$;};O{=9M31YO1wq$r@S|;Cm|{Xf1x2@dE-Y?hR5g z9{#E@lVcC>(W4`{877*#PoH!OR?s`q29BRNF*Px94HHfvq1;?s#J-)~ctZ?7v@7n$ zKEQTn>+g*x5Q|v+Cz!W>&g%}OX2kqte`5mVfe~O!3~yy*ULz{W()!P7>)GI&Hy3_H znTD;ZDJgCFgt5#LW@x(F+8Z`*gn*(NTZ|`pUQn=|iV7%!&Xzo#Ldk~s0b#E8_C9{E zg$d%H*T;|Qo}Qj5Y-FH%^NFOFOoGB z`=DlpOEx_(J6czV;62L>vD?W=QxjOl zX3Hz#iGY8hR>^o%6O$eyW&7%t8E_@cT%aUjYSjMe(_=wfhN${(4Cth#uP^i;0WpA% z$OEE4D@$4l;776u2xzORHNAfAM|uzu@#=>Ya-51Cv<0%^t=qS4dMhN@Q($Yr91ghT z)925X*g~)blpy8M+y@HC}*SY z-&aAfsj7x~h+y)7r}_N#`OB9^XU?3774iGaPe2y52B|jqzkh$mPnu2DV&leu(z-#b zMKwa$-s9xiC~7w3F|YtPSJy!YVW9lgWsfcTrOug9Sh{<)LFO>7-$_p&6&-yJ=xd;s zT0){cHkJu!urxd@3|y3^kcmKGqucs>b_7xPGGH4$>SGUvn||}Dfk~A&ljKgGk(W2e zvuU#ko#moCa`<5(;LN}QJbl)|qq2IHo$YYiSVA#5&A|D*^9?-ax(+T!ddU_BjqG+Drc>%^t ze?m|?uz&yMrnr@QzDmp;g_g&ucc5@elNWR!{~)L!pUF#P#Y@M5!0hbo#Ka5d&V8}C z?TN|m#%+7hEt=cf?%6GaDdHwG&70%h5K^ix&kGAN7DmlJM6Qo-YHLehyioVf81T1f zjNP(z>%$Y~w{F~!1GS>u_~pwN2jST0=)JOUQ7^FGO5_`u$ls(|U*femg)*Djdx6aA zFe1eO(6RB=ySH!q9E35dJg^!h>Up>HDs^}r?bZNL7NjB;WB5~cl%L;gP}5A?wdH9w zG*3LIp}xLK^f5n{tB8{MU%lz3c1Tuk?%M`FtILZwcODoW8v_SXyuJ4+Rx!)F{s9lv zX;>04(baAL$Gk}<`oRO&k#-&w&%c{2_f`Po z2@Fc|4U6;sC}a6lA^sC6mOw3-%EEwl7F_n!sZ&?=0N6f+-(kzaJiXNU3igtRn>&z( z*#jOik1~W0rtB_5AIXw5VQYZ)0I~!i*BZ0h(NRcr2OGblJXC>m*R(Su0>V&B|XbaqlP|c{=@rgs|- zn*y4ClBYV2jdq|mGuQv;Srj~NQ0MU{w2Fg7u>R%M-@m|$4iMq9(AL#OiTBPq4}Sna zD(*y+Z}axO=k4uH?VDO!G|kO_G{rMm68+D*6D&}B}TybL9P;qPZvtTSF*JV6Exhp8XEz#wDzxG%Q7_CwK0N|OCm1@%f<`STj%vx0(|I^EY%{~SLh zKXmnop{@cJfq8?+=!LqXVtRYl^0F6#&PrcJWekM~+nf}t{Toq~X%1f&HeO6bB!EbS za)h}(8K?Ud+nvu)c=zrdh-ybi$H`7DFo~rfd!qk%{r>LF>m}oRiLLd0tlZ$jd&RFk zfhiCE1yUYJTynk@A7*y|-JlR!S}@$D|3-VVw6y%lhuyQav4P0m)!oggDr2?0DG+!w z;n;2>>N55M0oPk{)uOzz391DKJzDj+5$qa~U}iReDX-`&v}x5sTG<7!isE87z*;;L zNIC{=6@l9lXb3I(M@9eg#9y#jqqBQu z3jKx8v?Arfz`%ffuQjZw+P?uXp{$T4z~tnsOEXs{M~a_6$B3IY-2!Hch3cyGe~C75 zGGKdYX}E&xIT>CJOm0{&(@p7v%0M0e<#0vQ-mHK2pBwxsWY6;7z0j-W2kp1>%TvLD z-lGNaX=h%iF@yb8NWEx4KmYe#)c;mF{;#y9|I(xV|LoO$cETiw6clpX-J1)5hvHUL zlD!}aeXtmH8DH9A>Tii@zz5?uZX;kuFtVr!v-wYikln9sZ35(9jEvkxHa+Ng#baWL zP(?F(6AvH0GJd4O5Sb3Gmi$|YyuY{R(BFQ^Modf~Y$Gcf#umT7qlIRTz2Bn! zMiwJX)B(&duZ+u9$btO8JzG>%^uhS;d}dM-xdtKn|NJc>A>qlBj|~lHV;ceSj~zdL z<;4#FCoCB1AWiK57}}M|^`kW3pA9`uxB9s3f9hzdh3UP*YQ0|{$Rvt0e5A>lnODla z<0QR3el{}04-E?38j2bQ=@2y5MLW6(lWwi*n7);__wOo(?9*@7@JhRN?%A0_t*lIy|CU4IxLV$~n9m08K^R#z#j+ zvQtxCAe)skV-N>qxU#Z>=fIztwGbVhfgmC4;sx4|jf%`&delEScp1G3znPio4hhx3 z8WjQ3(k5yKfeTd|P#`f2(;S2K6}ATF<#>pfXI@U&OXI>Ff85kg2-x2*& z*Tl2D8={Ystz*;Ett~A{l`Y?-2|n5&&_zX12%oBRi8~Le11Vd(^i5s^YKLh?HgN*e zRID;)vhcrv>cS;}YAl=R07o6}b65Tfqr>!|3N!O#wD3gr7g#-vJKIRo#{9}Wb+xq+ zibI#FSrj+BZUR9_H%-A3?P>{<9$X_-<-w_YxS2C8`4tf)?%EwX4cGvo{H> z5scK|fQA2Z0}P&>(io}PFXdn2VTWs4F zf_l9*po&P`KusNrl-U+gg~dQKqS~m;xNwvBmWtKyLCadEtzGIDY^ zR&_uU_1Fj#*ThNYSFQ}FDo0_;;!rmoF5)&G8x*v;5K-ZXXe#waqkyslVTdeSQ)BoF zi@5HuK!>?eeF)4}i)O9L@TZ~xM2Cm(=#B(v*|`%ITIJg8<;!F#Q%T9&#umTp>@09{ z4HOyoo;J9K8^e6Gt*H%vRaI8@octCL{p%*%GIE)ll9JM@1Tq6WF%X6dMTY>z{1nLw zLb~WBbMxVmk=9ToGhps6?&Ys+qGW?BM2f;Ry*5l+VPG1fG)6NEBOyb^4&hXchu3Z2!(ASR{3Z!0Eyztzh_3JBcsMOGTGEf!T zl=tn6K`X>4=L4u`Yb&kqe&0#4br3{ZvI>EHberQA_=?79?A)G}e^q5s!27Iu#D73q{7QSkZ zD4$+YI4If=aEngfHic(`uQXoVz3S)m7ItT$rK8k^H+?)A<5FR3QF|y5>{DV^M?>l* z0_Q=EMJ&5}c$_dEa~m}4)tp>Y%G~BjWKG2sUut&|u9=V^t)`Z*ySdpI=2vmj290Lw z>8U9k*2vb)J*YD)B)a-ClqOBLFn0r_jTQ4^Erp_y@DMOdn4dL5D$yjBs*L&~}d8bHJeaW3FlSk8LUb$#duzyKHnnCJi zW@Uk$X^epd{cMUyX@$&!2%Lh4;PQuusj%Bxf9B1Sb9OF(#Q|+Q9e%s?^e9aIZQ}cn z&MhxqJ$bT_j8M-#kp}I+0E=vf*SvJ;V4)F)WH?Mj%5gHl;1Le-m#$--_;5HOwK!Ad zcp#2VOiyd6tM3)CjEIkiaQz-1z|6#yK?B;0&OqMK7mfoB4GrYFD4QtF_N2I&RoJ1g zHk1+WN*rQ&S1QF%p>YB!Y&V#)~`IadsXFn`ys1e8lkg^OiQhis8p1(*pP+APjsY@j5yyA z2mS@tT}9mTg1M&dCB!|L6VS`h7iS8;qPvkxraKb~AJ%(na&qU&?j=A&w2u7DOw7-i zPg>E9>#65Mq@$;!`^X2?0xfs0@AVDHTI(n`sfWVnjkbwh1^sCI_%VUu=#e839C0k{ zQ<AA;$Vd`pTHHzy5S)}3%qj7qZj=6;Jx&w-9$N)&bI zqxb|b-}1@|N`a(V0VpR-rf9$Jo}L3G@~@cwd^)HysK2}WB5LW^&d!!Gz#TQ9LbM@n z?t5oB04Cb1IR+aaGt$#5_Bx93@u|&}gaf7GaK^i)KZg$?a)bV3fKq<`wLj=Ic}%CP zJICJx_;&mewh;RC0U@{aKG`pl0TCl?1iD=rQX4=0~%tP{n#QXIcCF86RjMMPwR0r!H4!kFap z7$=IIY^B6uu}c!yzq* z4G|vi6?u>cVA)z}XLokHa}Y5+g)C+Zl~wye)5DX>rl#?$=Fflvp=85-2_#O-2{$TE z5`~AW0aM!bR;U!>AesW2Dvruw62isxw!S{3Vv_s%CGt~5W0&Dz1vbu%LSkay$6M3G zE#pB7Duib)d;^1(s&;}1+VUZ3p57z;P?zChp42-Ty-Il#yzC`<6 zS3*;jW;<8mRf1!Sh3ww^tKp7LPB>iE8_A)p5XvAg#T+2Z%$>qQsUHHDdAAHsV8SB{ zGt~6p2S$Vg>bY*dNjIF@YLCHRVqwG5)zQJYz6{1M?E(CliHT__Q`m)cFF5#Jdf4`f z0FZ5!rZ;co_4M3;Sm0EH1BhA5$E#ZsS~ppUx3W4?NLlm42`BsHbDuqZ+EHjtA5}j) zG4Z6RC{3;n5r5{474!=9FyMMj7U~XU|%Q9)Ugr4H@L-2PvZdx`CBV-Rx62?4{QBgp)EiJpof%7#2I3h}*hdWE+SbU54-n5tOKM1ZI%i>fZEfIc zsA(Qb&3%^44&y(6{7634JF6=@_8l!VSP*u%UF4JM2Rav^7Qu|XG=`&A z-dy3=g+eh z-Z))I3<`P_6XPh%gc!h7O;}6}K*rRW2>olPEKuplN(zOk>X|QJJ{+dFzQs#Y3Xn1- zr3O6S*)3uUsub`OOSlCvC4&#dBh+FG z9mA*(5_G7E;#ww7?4D8!UEC#`sE#$&)0;;Ug}TD!09eA}QtXH5L0FhuShN-=d%_k4 z-DRRq^cb}#k>cYetbfL|He2hg?#u=jN>1hZ&NOYVe>FwIRZUn6n|5U{vZ|y(-eF`#8 aoAWK)aZkKWD-kv#7`WwCPvyy;_x~SA!xdNn literal 0 HcmV?d00001 diff --git a/shocks.py b/shocks.py index 4ff6ae6..612d5fe 100644 --- a/shocks.py +++ b/shocks.py @@ -27,8 +27,8 @@ def shock_calc(M1, P1, T1, P_a): nuM1 = 0 nuM2 = nuM1 + math.sqrt((gamma+1)/(gamma-1)) *math.atan(math.sqrt(((gamma-1)/(gamma+1))*(M1**2 -1))) - math.atan(math.sqrt(M1**2 -1)) #7.56 #Prandtl-Meyer function M2 = sp.optimize.newton(Prandtl_Meyer_Mach, 1.5, args=(nuM2, gamma)) #Prandtl-Meyer - P2P1 = ((1+(gamma-1)*0.5*(M1**2))/(1+(gamma-1)*0.5*(M2**2)))**(gamma/(gamma-1)) - P2 = P2P1*P1 + #P2P1 = P2/P1 #((1+(gamma-1)*0.5*(M1**2))/(1+(gamma-1)*0.5*(M2**2)))**(gamma/(gamma-1)) + P2 = P_a T2T02 = (1 + (gamma-1)/2 * M2**2)**-1 #T4/T04 from isentropic relations T01T1 = (1 + (gamma-1)/2 * M1**2) #T03/T3 from isentropic relations P1P01 = (T01T1**-1)**(gamma/(gamma-1)) @@ -58,8 +58,8 @@ def shock_calc(M1, P1, T1, P_a): M3 = M3n/math.sin(beta2rad-theta2rad) #Mach number after 1st shock - T3T2 = ((2*gamma*M2n**2 -(gamma-1))*((gamma-1)*M2n**2 +2))/((gamma+1)**2*M2n**2) #1.29 #T2/T1 from normal shock - T3 = T2*T3T2 + 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")