def OBJECTIVE_Cratio(x):
Run = solveCycle(
"N2",
"Zerol",
"Zerol",
ps,
x[0],
x[1],
TL,
298,
pr,
0.9,
eta,
eta,
eta,
eta,
0.0,
1.0,
1.0,
use_etac_fit=True,
)
return 1 / Run.Cycle.etaII
def ConfigE():
Tc = 233
Th = 298
eta = 0.8
reCalc = False
N = 50
prv = np.linspace(1.5, 4.5, N)
psv = np.linspace(300, 700, N)
(PS, PR) = np.meshgrid(psv, prv)
ETAII = np.zeros_like(PS)
CAP = np.zeros_like(PS)
ETAC = np.zeros_like(PS)
XL = np.zeros_like(PS)
for TL in [223.0, 275.0]:
if reCalc == True:
for i in range(len(prv)):
pr = prv[i]
for j in range(len(psv)):
ps = psv[j]
eta = 0.70
def OBJECTIVE_Cratio(x):
Run = solveCycle(
"N2",
"Zerol",
"Zerol",
ps,
x[0],
x[1],
TL,
298,
pr,
0.9,
eta,
eta,
eta,
eta,
0.0,
1.0,
1.0,
use_etac_fit=True,
)
return 1 / Run.Cycle.etaII
x0 = [10.0, 6.0]
b = [(1.0, 15.0), (1.0, 10.0)]
(x, f, d) = fmin_l_bfgs_b(OBJECTIVE_Cratio, np.array(x0), bounds=b, approx_grad=True)
Run = solveCycle(
"N2",
"Zerol",
"Zerol",
ps,
x[0],
x[1],
TL,
298,
pr,
0.9,
eta,
eta,
eta,
eta,
0.0,
1.0,
1.0,
use_etac_fit=True,
)
print x, Run.Comp.xL, Run.Exp.xL
print ps, pr, Run.Cycle.etaII * 100, Run.ColdHX.Q, Run.Comp.eta_c, Run.Comp.T_in
PS[i, j] = ps
PR[i, j] = pr
ETAII[i, j] = Run.Cycle.etaII
CAP[i, j] = Run.ColdHX.Q
ETAC[i, j] = Run.Comp.eta_c
XL[i, j] = Run.Comp.xL
foldPath = os.path.join("ConfigE Joined", "%0.0f" % (TL))
np.savetxt(os.path.join(foldPath, "PS.csv"), PS)
np.savetxt(os.path.join(foldPath, "PR.csv"), PR)
np.savetxt(os.path.join(foldPath, "ETAII.csv"), ETAII)
np.savetxt(os.path.join(foldPath, "CAP.csv"), CAP)
np.savetxt(os.path.join(foldPath, "ETAC.csv"), ETAC)
np.savetxt(os.path.join(foldPath, "XL.csv"), XL)
else:
foldPath = os.path.join("ConfigE Joined", "%0.0f" % (TL))
PS = np.loadtxt(os.path.join(foldPath, "PS.csv"))
PR = np.loadtxt(os.path.join(foldPath, "PR.csv"))
ETAII = np.loadtxt(os.path.join(foldPath, "ETAII.csv"))
CAP = np.loadtxt(os.path.join(foldPath, "CAP.csv"))
ETAC = np.loadtxt(os.path.join(foldPath, "ETAC.csv"))
XL = np.loadtxt(os.path.join(foldPath, "XL.csv"))
levelsCoarse = []
levelsFine = []
if TL == 275:
levelsCoarse = [5, 6, 7]
levelsFine = [5.5, 6.5, 7.5]
elif TL == 223:
levelsCoarse = [12, 15, 18]
levelsFine = [13.5, 16.5]
fig = pylab.figure(figsize=(3.5, 2.5))
ax = fig.add_axes((0.15, 0.15, 0.8, 0.8))
CS = ax.contourf(PS, PR, CAP, cmap=pylab.get_cmap("Blues"))
cb = fig.colorbar(CS)
cb.set_label("Capacity [kW]")
CS = ax.contour(PS, PR, ETAII * 100, colors="k", levels=levelsCoarse, linewidths=0.8)
ax.clabel(CS, fontsize=8, inline=1, fmt="%1.1f %%", inline_spacing=1)
CS = ax.contour(PS, PR, ETAII * 100, colors="k", levels=levelsFine, linewidths=0.2, linestyles="dashed")
ax.clabel(CS, fontsize=8, inline=1, fmt="%1.1f %%", inline_spacing=1, manual=True)
ax.set_xlabel("$p_s$ [kPa]")
ax.set_ylabel("$p_{ratio}$ [-]")
fig.savefig("ConfigE" + "%0.0f" % (TL) + ".pdf")
fig.savefig("ConfigE" + "%0.0f" % (TL) + ".eps")
fig.savefig("ConfigE" + "%0.0f" % (TL) + ".png", dpi=400)
pylab.close()
## fig=pylab.figure(figsize=(3.5,2.5))
## ax=fig.add_axes((0.15,0.15,0.8,0.8))
## CS=ax.contour(PS,PR,ETAC*100,colors='k')
## ax.clabel(CS,fontsize=10,inline=1,fmt='%1.1f %%')
## ax.set_xlabel('$p_s$ [kPa]')
## ax.set_ylabel('$p_{ratio}$ [-]')
## fig.savefig('ConfigEetac'+'%0.0f'%(TL)+'.pdf')
## fig.savefig('ConfigEetac'+'%0.0f'%(TL)+'.eps')
## fig.savefig('ConfigEetac'+'%0.0f'%(TL)+'.png',dpi=400)
## pylab.close()
fig = pylab.figure(figsize=(3.5, 2.5))
ax = fig.add_axes((0.15, 0.15, 0.8, 0.8))
CS = ax.contour(PS, PR, XL * 100, colors="k", levels=[65, 70, 75, 80], linewidths=1)
ax.clabel(CS, fontsize=8, inline=1, fmt="%1.1f %%")
CS = ax.contour(PS, PR, XL * 100, colors="k", levels=[67.5, 72.5, 77.5], linewidths=0.2, linestyles="dashed")
ax.clabel(CS, fontsize=8, inline=1, fmt="%1.1f %%")
ax.set_xlabel("$p_s$ [kPa]")
ax.set_ylabel("$p_{ratio}$ [-]")
fig.savefig("ConfigExLopt" + "%0.0f" % (TL) + ".pdf")
fig.savefig("ConfigExLopt" + "%0.0f" % (TL) + ".eps")
fig.savefig("ConfigExLopt" + "%0.0f" % (TL) + ".png", dpi=400)
pylab.close()
