def ConfigB(): f = pylab.figure(figsize=(4, 2)) ax = f.add_axes((0, 0, 0.5, 1)) geo = ps.geoVals() ps.plotScrollSet(7 * 3.141592654 / 4, discOn=True, lw=1, geo=geo, axis=ax) pylab.gca().set_xlim((-0.017, 0.010)) pylab.gca().set_ylim((-0.021, 0.006)) ax = f.add_axes((0.5, 0, 0.5, 1)) geo = ps.geoVals() ps.setDiscGeo(geo, "2Arc", 0.001) geo.disc_x0 = -0.00269898 geo.disc_y0 = -0.00348883 geo.disc_R = 0.01187952 ps.plotScrollSet(7 * 3.141592654 / 4, discOn=True, lw=1, geo=geo, axis=ax) pylab.gca().set_xlim((-0.017, 0.010)) a = -0.004 pylab.gca().set_ylim((-0.021 - a, 0.006 - a)) pylab.savefig("ConfigBDischargeBlockage.pdf") pylab.savefig("ConfigBDischargeBlockage.eps") pylab.savefig("ConfigBDischargeBlockage.png", dpi=400) pylab.close() path = "ConfigA\\Run 0001" A0 = np.loadtxt(path + "/" + "Adisc.csv", delimiter=",") path = "ConfigB\\Run 0001" A1 = np.loadtxt(path + "/" + "Adisc.csv", delimiter=",") path = "ConfigB\\Run 0005" A2 = np.loadtxt(path + "/" + "Adisc.csv", delimiter=",") path = "ConfigB\\Run 0009" A3 = np.loadtxt(path + "/" + "Adisc.csv", delimiter=",") pi = np.pi axes_size = (0.18, 0.20, 0.77, 0.74) fig = pylab.figure(figsize=(3.5, 2.5)) pylab.axes(axes_size) pylab.plot(A0[:, 0], A0[:, 1] * 1000 ** 2, label="$p_s$") ## pylab.plot(A1[:,0],A1[:,1]*1000**2,label='$p_s$') pylab.plot(A2[:, 0], A2[:, 1] * 1000 ** 2, label="$p_s$") ## pylab.plot(A3[:,0],A3[:,1]*1000**2,label='$p_s$') pylab.text(0.9 * 3.14, 120, "Baseline", ha="center", va="bottom") pylab.text(0.8 * 3.14, 372, "Larger port", ha="center", va="bottom", rotation=7) pylab.xticks([0, pi / 2, pi, 3 * pi / 2, 2 * pi], ("0", "$\pi$/2", "$\pi$", "3$\pi$/2", "2$\pi$")) pylab.xlim([0, 2 * pi]) pylab.ylim([0, 500]) pylab.xlabel(r"$\theta $ [radians]") pylab.ylabel(r"Discharge Port Area [mm$^2$]") pylab.savefig("ConfigBAdisc.pdf") pylab.savefig("ConfigBAdisc.eps") pylab.savefig("ConfigBAdisc.png", dpi=600) pylab.close() path = "ConfigB\\Run 0005" pi = np.pi axes_size = (0.15, 0.18, 0.8, 0.77) theta = np.loadtxt(path + "/" + "theta.csv", delimiter=",") p = np.loadtxt(path + "/" + "p.csv", delimiter=",") p = np.where(p < 1e-10, np.nan, p) # Plot data fig = pylab.figure(num=1, figsize=(3.5, 2.5)) pylab.axes(axes_size) ##pylab.subplot(131) pylab.plot(theta, p[1, :], label="$p_s$") pylab.plot(theta, p[9, :], label="$p_{c1-1}$") pylab.plot(theta, p[2, :], label="$p_{d1}$") pylab.plot(theta, p[6, :], label="$p_{ddd}$") pylab.plot(theta, p[5, :], label="$p_{dd}$") pylab.plot(theta, p[7, :], "k--", label="$p_{dd}$") pylab.plot(theta, p[8, :], "k--") pylab.text(3 * pi / 2.0, 500 + 10, "Suction", ha="center", va="bottom", size=8) pylab.text(pi / 2.0, 1850 - 10, "Discharge", ha="center", va="top", size=8) pylab.text(7 * pi / 4.0, 440, "$s_1$", ha="center", va="top", size=10, color="b") pylab.text(pi, 830, "$c_1$", ha="left", va="top", size=10, color="g") pylab.text(3 * pi / 2, 1150, "$d_1$", ha="left", va="top", size=10, color="r") pylab.text(3 * pi / 2 + 0.2, 1800, "$dd$", ha="center", va="top", size=10, color="m") pylab.plot(np.r_[1.36 * pi, 1.36 * pi], np.r_[900, 2200], "k--") pylab.text(1.36 * pi - 0.05, 1450, "Discharge Angle", rotation=90, size=7, ha="right", va="center") th = 1.75 * pi pylab.plot(np.r_[th, th], np.r_[900, 2200], "k--") pylab.text(th + 0.05, 1000, "Merge", rotation=90, size=7, ha="left", va="center") pylab.text(3 * pi / 4, 1950, "ddd", color="c", ha="center", va="center") pylab.xticks([0, pi / 2, pi, 3 * pi / 2, 2 * pi], ("0", "$\pi$/2", "$\pi$", "3$\pi$/2", "2$\pi$")) pylab.ylim([300, 2200]) pylab.xlim([0, 2 * pi]) pylab.xlabel(r"$\theta $ [radians]") pylab.ylabel(r"Pressure [kPa]") pylab.savefig("ConfigBPressure.pdf") pylab.savefig("ConfigBPressure.eps") pylab.savefig("ConfigBPressure.png", dpi=600) pylab.close()
def ConfigA(): r = prep("ConfigA\\Results.csv") ## print r.dtype.names w = 3.5 h = 3.5 d = 0.15 f = pylab.figure(figsize=(3.5, 3.5)) ax = f.add_axes((0.1, 0.1, 0.8, 0.8)) data = np.r_[r["leakageflank"], r["leakageradial"], r["suction"], r["discharge"], r["mechanical"]] labels = ["Leakage\n(flank)", "Leakage\n(radial)", "Suction", "Discharge", "Mechanical"] (path, labelText, pcts) = ax.pie(data, labels=labels, autopct="%0.1f%%", colors=["yellow", "g", "r", "c", "m"]) [t.set_size(8) for t in labelText] [p.set_size(8) for p in pcts] [p.set_linewidth(0.5) for p in path] f.savefig("ConfigALossPie.pdf") f.savefig("ConfigALossPie.eps", set_bbox="tight") f.savefig("ConfigALossPie.png", dpi=600) pylab.close() path = "ConfigA\\Run 0001" pi = np.pi axes_size = (0.15, 0.18, 0.8, 0.77) theta = np.loadtxt(path + "/" + "theta.csv", delimiter=",") p = np.loadtxt(path + "/" + "p.csv", delimiter=",") p = np.where(p < 1e-10, np.nan, p) # Plot data fig = pylab.figure(num=1, figsize=(3.5, 2.5)) pylab.axes(axes_size) ##pylab.subplot(131) pylab.plot(theta, p[1, :], label="$p_s$") pylab.plot(theta, p[9, :], label="$p_{c1-1}$") pylab.plot(theta, p[2, :], label="$p_{d1}$") pylab.plot(theta, p[6, :], label="$p_{ddd}$") pylab.plot(theta, p[5, :], label="$p_{dd}$") pylab.plot(theta, p[7, :], "k--", label="$p_{dd}$") pylab.plot(theta, p[8, :], "k--") pylab.text(3 * pi / 2.0, 500 + 10, "Suction", ha="center", va="bottom", size=8) pylab.text(pi / 2.0, 1850 - 10, "Discharge", ha="center", va="top", size=8) pylab.text(7 * pi / 4.0, 440, "$s_1$", ha="center", va="top", size=10, color="b") pylab.text(pi, 830, "$c_1$", ha="left", va="top", size=10, color="g") pylab.text(3 * pi / 2, 1150, "$d_1$", ha="left", va="top", size=10, color="r") pylab.text(3 * pi / 2, 1700, "$dd$", ha="center", va="top", size=10, color="m") pylab.plot(np.r_[1.36 * pi, 1.36 * pi], np.r_[900, 2200], "k--") pylab.text(1.36 * pi - 0.05, 1450, "Discharge Angle", rotation=90, size=7, ha="right", va="center") pylab.plot(np.r_[1.72 * pi, 1.72 * pi], np.r_[900, 2200], "k--") pylab.text(1.72 * pi + 0.05, 1000, "Merge", rotation=90, size=7, ha="left", va="center") pylab.text(3 * pi / 4, 2050, "ddd", color="c", ha="center", va="center") pylab.xticks([0, pi / 2, pi, 3 * pi / 2, 2 * pi], ("0", "$\pi$/2", "$\pi$", "3$\pi$/2", "2$\pi$")) pylab.xlim([0, 2 * pi]) pylab.ylim([300, 2200]) pylab.xlabel(r"$\theta $ [radians]") pylab.ylabel(r"Pressure [kPa]") pylab.savefig("ConfigAPressure.pdf") pylab.savefig("ConfigAPressure.eps") pylab.savefig("ConfigAPressure.png", dpi=600) pylab.close() geo = ps.geoVals() ps.plotScrollSet(theta=3 * pi / 2.0, discOn=True, lw=1) pylab.savefig("ConfigAScrollSet.pdf") pylab.savefig("ConfigAScrollSet.eps") pylab.savefig("ConfigAScrollSet.png", dpi=600) pylab.close()