plt.plot(CWalpha, CWm3, '-', color='#bb0000') # --- (UC) unstarted channel (A8 is sonic) alphamax = 1.1*alphamax # M8 is sonic so M4 is known UCm4 = mf.Mach_Sigma(A3A2/A8A2, .1, gam) # look for subsonic value UCm3low = mf.Mach_Sigma(A3A2, .1, gam) # look for subsonic value alphamin = ray.Ti_Ticri(UCm4, gam)/ray.Ti_Ticri(UCm3low, gam) print " (restart) choked throat for Ti4/Ti0 = %6.3f"%(alphamin) A1A2 = mf.Sigma_Mach(sw.downstream_Mn(M0, gam), gam) print " critical flow for unstart condition if A1/A2 = %6.3f"%(A1A2) print " A1/A0 = %6.3f"%(A1A2*A2A0) UCalpha = np.log10(np.logspace(alphamin, alphamax, npts+1)) UCpi0 = sw.Pi_ratio(M0) UCm3 = ray.SubMach_TiTicri(ray.Ti_Ticri(UCm4, gam)/UCalpha, gam) UCpi3 = UCpi0*np.ones(npts+1) UCpi4 = UCpi3*ray.Pi_Picri(UCm4, gam)/ray.Pi_Picri(UCm3, gam) UCm4 = np.ones(npts+1)*UCm4 plt.subplot(311) plt.plot(UCalpha, UCpi3, '-', color='#990000') plt.subplot(312) plt.plot(UCalpha, UCpi4, '-', color='#990000') plt.subplot(313) plt.plot(UCalpha, UCm3, '-', color='#990000') # --- (RC) unstarted but choked channel (A2 and A8 are sonic) alphamax = alphamin
