From cdf81e08e338b14edbb755b10695a3c6287bbf0c Mon Sep 17 00:00:00 2001 From: Nils Gaukroger Date: Tue, 22 Jun 2021 17:49:54 +0200 Subject: [PATCH] Pitt_Peters input details --- .../__pycache__/utilities_DI.cpython-38.pyc | Bin 18814 -> 20418 bytes A1_Dynamic Inflow/main_uBEMT.py | 8 ++++++++ 2 files changed, 8 insertions(+) diff --git a/A1_Dynamic Inflow/__pycache__/utilities_DI.cpython-38.pyc b/A1_Dynamic Inflow/__pycache__/utilities_DI.cpython-38.pyc index 947efa976c4fb5658e14c1ed6f7d8d2280c92a53..847b543bae785eaf0f16a7133becd5f63dad25e6 100644 GIT binary patch delta 6286 zcmZuVTWlQHb#rHCcW0O6l2S;M;zJr!lBJbR%A{o_(OX+76T^Y0K%zuI$iyFU5(EmEK3o@V5uk0-HV-!+t*5n(Hbx!iBZ!L> zZGYOHbB6mcogw$kx#zylJ?}gJ`x5)X*I4@ZIyzDUJWu`S>!11UQ{PL!EHiC;Tb3Bp zf_o(aFU4ailjV(6Emq z?L-4SjM!+fSGmAG8q6uj`cD#}(cdz1`%B76#?;^^$?3MEIN%6^;b1s5_R(X2hGKRe zezWkqnvtW;58#jtg7XMkR>S7W?dMbfDArWOFCaLKU>Lz61b8GmVJN%Sy2&qM`x1gj z0ECiTSq_sP*KN1#xF#RM9y%YLc@~==LvR(rH2|4p#Fh+r$UIFla-1J~2wGaq$m#1= zb9F6Y%sYIgVzx3Jq|y@?PT^{&5u8Cl0(%^>=IIC{(Vm44WW>cS6{f`b+29N5p?hvZ z(-bsj2%2axLws5QW)1{@k=_Rqej`10q}S*JTAOjNvCrs#nrZFEfH4S^4r9pJ5AU>* zniDd^!M&ZQ`$^BZK?@7C`^oKUXFD_AeTbf4oadb|DAUGiEQy9hOyW|Z3q)$!WKN;v z={9din(VG{FHGn>52SIJP1}(!CZIpk*RkN=yYB7EV!J7ZhlAcdZ?K-={XIRCDlQeb zY|rH7NT)HKbNxod%Lw$=ZKewsnM5#!pmn)V?0vg1heSF@ESgIiSCRSz0^~W9ikKub zC8qM&)mSu?OIE=PWrOQQvZYAxty{hLk&&7Nlj@HG2>L0~-0uT;jx7mG%wX3=zF>$O zk|BBGrsPXotXuG@i}`34OvAGo%4{Z^TyaA_Yep@#ah;-X{zQ zz0hyZe99O3A3V3yZNz+$wveNZX&bRm=%z)xVG>Lkd{%gt6+|OpC{HUJaU;2r@M3=4 zPuvnWVtxz^CZ*pAe!u&#wm;wVYsTh+_1+`3g zFA8SO1yRZn>}Wo%gyt5?+48yU+ojODc9-WhXs z+BPbA&vNWbdfs#R9q4rX00>=s+nLPvGKKYtsw}h7?O*g=Wi<*vD)fhvxmJ#X3Nc1O z_aJt}!58w5%kRTR%yMnrt}{!sSAec1mMb~#I3D+)#itQ`0>NbjEtUq&uzMDQnBkS? zt#Vt|;oKyKe-Qk6e}SC~4i7xaeiQJ4BC7`Bz%%UE!Dj|9vi}agKKRPMabMW%_xJl< zTjIt6d)AkXj_cwVMSe3_8u}!=6a3N8R7dC?*i`)kTMX!L1&SHYLreWwquav6zVB}AHwD+ z0utGN#2QlNCN`V-12YFV5TIyz5y29IAq3>x6rLzhxj+lsW-ZUl&DKl^C2o`xNT#be z0ae1QI84$PT>dCFA46~j0mZ8XSHb2EpS1GKa_-L)o*IO)ZsoYPMma zS-F}rIh}jVG0c*-cSn!8d_##^Y9?m^r@3W)EsVh+Pp5>a#k`zltc9ZC@rQ5|l2Ah~ zAzR`knr56OesVZ;0fCuE_`43Xh#eLJOp{S)BPT+;jO?#rj#^wVqa*FM$oQO7s+#=g z!7qo;v2<|c;Lom-Xi-V^=-PBl#`$M)Z8D%3_~(kN$wt-iLRRIp>jw><`!M=kaQV>T zngUq;*#k%net0PNE&woY!56ofpXe4g6hC?=F<%;z4{;V_r4LbfN%Ru~5WqH4ek^Ja z2!{L-!A}9LjaGrUOslj*fJ?zF1p+Gc*8Fx~=@a}oPVa?)jS+RD!38mo>GHe$bU}pp8lNeC6#f|DZHeu}eZuLr6ZkU%CNO`y$Ru+#{sXex?OiVpCKMLc5 z(xK|(%?xlIS)nu_ausn5j+uBAtI=AuM6skp-Q-0%po}C*&yDTkfjVP1bP_^l{)16Usm^bTl&!E~l6Bm6Gn629&nV$Z{2jk||P!*A+2! z0jY-($OtF~QAoNOym<7U+B^~|Y~f6Z;kIF}HS5GMnJwp{T(wj0gv$KfwOq+^y)fl1 znx0O0G^4S(mW;X_UW|OY$>oc%a2k*Ttiii1}8EjwYa2+>m4NfdT(BW z$@K@ba0-6pBVPx^0}n9REs?j|iU*!?i`|4|z62h!PuP$MlL3AW z#cr`l;mJcA2|vCW^WvL{nlFM2Nk-g=KLbux5cdgN;5tlLnSgG^k8iOJg)kYgq_1oy zx6thrKLtxG;1l>x;!@zQ(B6{3D-!yQ2yCl6&ootm4(d_`iK%7F%P*>wbDmHabk{Ug z$5wUKFkOp7!&~GP*HiP3StwX}%e1{wb=p%auIbYJ?Ci)q>SP3T1fy)vsc=^CC9NgrT z`r&o1lyN)A2zp8Ntg`BsP1u`As#73F&1y|ARcMe}cEDXnQ;lV>=@y))(Hx^b5DM7m!`UmR$YOpp;Upp7i2=ah&-S&xw%pyt!3w1PRG#!)~hp43s+U0TQz;5 zWU2+H0{p16B@>jd@`|l~zyXfdeeVI-V8`_`xL;kxs^zJ!qfV=N-Bv-YaK|9)8dc}v zq}SajK2Q)qmz>F@!o;HP>Yi8QVX`4xVlujVV(5T6uUs~5BZ?>{P-$KGziR~ut5I~V&1Beo zrjwt7#dcq6sVZw0WR$e1$hwJxF@^Kz5nv$aFCq9cf<}~!FgynTG=O#ZQRR)ir8tg3G;xtT-A-fhudn^b;HW3v{cwr1Btf|-( z=?9DQBD8ldc|qa+0|Uolke0sYK-r&jxK*^kHKT&DKiIs#r@P51F?kY(|7!b3_v>N} z(_N-3D)40;ufoL7SAZNPWK?%#7an3?MSzJgs$wgMO(LM{(G+5MInAV_M9+Zjd`mS- zYcy)Ike$AyiP;@et>Y(Z z-@<)K&te!~Kp%6%1XOhk@UfyuLpm27AXxW=>)MX#{5wpzfpNZ6XFE zTW6D@kd54ij`u*MF4nDJ$i2#brGwwzeqnr=)jkijjL09v$zuraB6tO<5{V5&RSv~M z`NLHpZLo-SlmmeC)zaS%8s9D!P}>A_mHY6?qZo1N=v54?+IQw z^I!2^Wc}sf`LlCuHhAmoyTO}h6T!s!f#BV9e>zIaYbv~Xjd>JjT|sae!6bsI;HT%m z0Ut}<3*X-J8W2LcT!L%&F!mh{j$a(Bp-qK}ox*^B3S04fNnbEaVa&5Uxaks$AIBN{ z5MbPiE}$PnjFf{W8YUo^Le}8$q0Hlv{h|$KdbDIc%oH1F+0)vtfk)K9B_RG3qZT!w ztD}*gsjszJ_9APM3RV%mrq1oY24f_oFPotRIW&yZQC!L_p=4O_Jp@3`DZ}+k(_&f| zqOIBt?1@3~aWVM9#5>~K!Zy$TlAXhf`@I-0{AU2x6{v=emMlBNe}qkOWMp0O^o3Ik zrHUDK!UvEH82b-#%HIZGp1REbI{4SAXXF2dgt@@Jbei1_zIEwT9POX~P4KHr=d*tY zXsEcA1=Ntse~itPYkIjGbOio75@Oz}Wv@YShfk4r8Y9-C5kCn^j~uN11CaU85s<>; zC8r4g{HI9y82|`we+aKA=UO(sfnTk?Zrf13+bu2Z+xtNep%J@jI4)uGET>+;IC!SZ z-+;YB5pt2FuYKsM)q`ivErajmz^Ura^G7^RoH~MMMk@d>sLW3yKI7Gj1d8RRkoXw-Li#TGSr}vtV-i Qc-u;PPU)9Lr6zRzKh8+pF8}}l delta 4702 zcmZu#eQ;D)72o^b+ufH9`5*)ml8`LlEafZI0s|fIH(y)zQ|m_OG^7N99&66-J9Xj?Q3Pow4O+ zei!@$dvF0v>51vQ34VDq!eWGF2rUSU0Aw&dVkbeLy4aYVWVKcVY zi#q)t_DRvNuWs5x14Vh+%)*!TO^huTuT|`;jNyb;2usD1Nc-kHpedEE0r>X9_i#+} zH;>|!)d)Kg%2vbfmcr*Fzg6;dkWL&}fv^-|IRYMuPAE0knXq{m_V*&}1CXkl9g!6t zH&bpV?b^H*C3HSIb02oMBJ?8+0K_VMwKQl;=1J8uIA47i_JWpSPB>+&hpmA$Kbf`5 zl@^ie&U?1wYIh^-Kp=xXh^^A;s70pj06{N&UEESn;L z8}}st+#%kpU0k~h`=w)V7AxmoVG9JSo8L{#hn$pW^NjDqpvk#8?(4D+9-N7hT3lo| z!k*d74c5IqbOZ#9cu%(~YC&O07r0%J?a8CcMG8xK+AW$Oa6>*{d7k7@Er=nY<4*;5} zMAD9$Ui8$k6CaMAf;NY5%1n+=IH^R`8@8i7?WK9))AK)I?1(thu#&A2lMM^SR~wr4 zWWWTx5`lc7`t)Ps0C|8>c0Q+YN;{&OVirh!d(s>^W|=#FEP{=9KcLNen4^(ra`HTP$2}*V`nh_c(3<+-*kxxo3RTx%|HTsf}hYA>jhZK+-eimDUJS%puDzD4`l*ToZyc2!;V z+RjIFZ5J4fUl;#eB<7BtWhX~zv;{^Xt_qK}zQuN8(#r~n2ZJJkFT{Rj#_{6A{zy#W z81RXQ>%X(-or%Xf_yXv0O20hL6&IsB^5j?|bYwOjMIhED!d{PEGH4sNiZ>ffm zj6kGCq%ezVij}0yjR*0GUIe7{*fMgz)I4XzmXUZmpRyeEp)2+7ea6AGWhZ6bw3m%| z(Zw>-V>qCxF=CEEWwDcqFp3A&dD?w`M%zyPaL~L2sZKB z?GBSf@d&~`1Ueg~5=N~boyA~yp8^IwZMmh52lqGstEhsm_)q#25( zHBmpwIR7H9O@S4JG(G78!;V^>WTTvRy)_;9BQU!|Y*@Bj{Cj!!{t}_2ag`bXh(O|X z%*$A)RCi=Jku;$adp@~%AO=qx?X&dfkAYyXczMNk<-dV*EB{q^f8{f*ssRQ$!xb5C z5zohBd2pnejAxUkXIpSfD#a_^6X<|+S*Y=~DcwvoT#6!lJ+{sP$O>+cm`=*FxeR7F z)H{A3$|yTvnLH8MZysNb`gr8%4cjMOVh~33S$prQVe} zK+M_)Cv?|{U#?!8*OiH?&Rquw*3HP^;Qsl0Vr(E*#cy|kp=2!OA?6&X+`>bW=g7MJ zDeRW~N-es290dVVqT5st2sa%yQlp3TipWZRZsQ5~jqsy|9UIoN4hn9}Y{~k3<73uK z*6C*hUJay*iNki01CW;!(4Q{n;#=))-AKx^x(~}07|cmHP{w>#K%&SYlTkvvB`}ZUm`%pxJZmu&h3`;qw6rqH(a5HyBB297&ao{YsejBfAzhdIPEo*yj7vywW)Sbi3@}MZ4>}e33#gpmK zfGBRzTlFgbbm8)rB`p6n7z2)K$A>o}yoB%~j;druuv|teLz$1Colg7&PS+4HTlf@! zoRhKn5R{HA-js<_O5DNE;4}xPX+ue;BNL0raF0RvXc_L$m~ptc5hu!Od(4h!u{0S} zt5hsAReplF{zNG*ZEaX1^&TTW46jS6*`ZZ)&~}IW?UV}@Zq)I{uQ73RYm~iSn7{33 zA&9fz?Ch+R`l0D;lVWM-SLTzeh!K&@q$hY1o7f%HZHSY|eIS0~!eh&Ir6gTSmf${}6b5{39mb z*n3?Gu)_EDy~~=gNZtsdCEf&>(7OhrtCCJC#{VRKu)j5Q6(!fi@AvOze-z6PoC&Fr zBK(N>?tyLWCGn2~TSJ%x{B5!BU`N+=XiDA99>c9&{%7n)T-!6o=>Yt%I1uzsIjJ$= zd3gVPwm4(VpYa!Q@!-<@-(Z;k9f8~(N~RO=$KS<~8vwC6{1y0-A-H;guJjAncR6P! zz-hG9Ls{LeX+kq*kKw$q%~x@fLW@Isy8Hv&OM$YanqwQw%`7H*n%Eb_<2_x=o53ux z1{pRCD;+ls8BAqIGUL-#ECcTXrpwoe|Mj#5$w*@Hq39lBK%@eG8bL##+XHdL)7X6- ffy@y=yW=ip3i#wDrVg#IY^?6rn>0nw2de%L5!hpc diff --git a/A1_Dynamic Inflow/main_uBEMT.py b/A1_Dynamic Inflow/main_uBEMT.py index 444f3bd..bd9fe22 100644 --- a/A1_Dynamic Inflow/main_uBEMT.py +++ b/A1_Dynamic Inflow/main_uBEMT.py @@ -45,6 +45,14 @@ CT_step['results'].append(Calc.Results) +#%% A.2 - Sinusoidal change in quasi-steady thrust coefficient + +#Define time array and reduced frequency +time_arr = np.arange(0,10,0.01) +omega = np.arange(0.05,0.301,0.05) + + +