extWallIndoorSurface,
intWallIndoorSurface,
intGainsConv, intGainsRad, ports,
model,
dt=int(3600/times_per_hour))
# Compute averaged results
T_air_c = T_air - 273.15
T_air_mean = np.array([np.mean(T_air_c[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])
T_air_1 = T_air_mean[0:24]
T_air_10 = T_air_mean[216:240]
T_air_60 = T_air_mean[1416:1440]
# Load reference results
(T_air_ref_1, T_air_ref_10, T_air_ref_60) = tc.load_res("inputs/case05_res.csv")
T_air_ref_1 = T_air_ref_1[:,0]
T_air_ref_10 = T_air_ref_10[:,0]
T_air_ref_60 = T_air_ref_60[:,0]
# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
plt.figure()
ax_top = plt.subplot(211)
plt.plot(res, label="Reference", color="black", linestyle="--")
plt.plot(ref, label="Simulation", color="blue", linestyle="-")
plt.legend()
plt.ylabel("Temperature in degC")
plt.title(title)
times_per_hour))
# Compute averaged results
T_air_c = T_air - 273.15
T_air_mean = np.array([
np.mean(T_air_c[i * times_per_hour:(i + 1) * times_per_hour])
for i in range(24 * 60)
])
T_air_1 = T_air_mean[0:24]
T_air_10 = T_air_mean[216:240]
T_air_60 = T_air_mean[1416:1440]
# Load reference results
(T_air_ref_1, T_air_ref_10,
T_air_ref_60) = tc.load_res("inputs/case03_res.csv")
T_air_ref_1 = T_air_ref_1[:, 0]
T_air_ref_10 = T_air_ref_10[:, 0]
T_air_ref_60 = T_air_ref_60[:, 0]
# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
plt.figure()
ax_top = plt.subplot(211)
plt.plot(res, label="Reference", color="black", linestyle="--")
plt.plot(ref, label="Simulation", color="blue", linestyle="-")
plt.legend()
plt.ylabel("Temperature in degC")
plt.title(title)
extWallIndoorSurface,
intWallIndoorSurface,
intGainsConv, intGainsRad, ports,
model,
dt=int(3600/times_per_hour))
# Compute averaged results
T_air_c = T_air - 273.15
T_air_mean = np.array([np.mean(T_air_c[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])
T_air_1 = T_air_mean[0:24]
T_air_10 = T_air_mean[216:240]
T_air_60 = T_air_mean[1416:1440]
# Load reference results
(T_air_ref_1, T_air_ref_10, T_air_ref_60) = tc.load_res("inputs/case12_res.csv")
T_air_ref_1 = T_air_ref_1[:,0]
T_air_ref_10 = T_air_ref_10[:,0]
T_air_ref_60 = T_air_ref_60[:,0]
# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
plt.figure()
ax_top = plt.subplot(211)
plt.plot(ref, label="Reference", color="black", linestyle="--")
plt.plot(res, label="Simulation", color="blue", linestyle="-")
plt.legend()
plt.ylabel("Temperature in degC")
plt.title(title)
windowIndoorSurface,
extWallIndoorSurface,
intWallIndoorSurface,
intGainsConv, intGainsRad, ports,
model,
dt=int(3600/times_per_hour))
# Compute averaged results
Q_hc_mean = np.array([np.mean(Q_HC[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])
Q_hc_1 = Q_hc_mean[0:24]
Q_hc_10 = Q_hc_mean[216:240]
Q_hc_60 = Q_hc_mean[1416:1440]
# Load reference results
(Q_hc_ref_1, Q_hc_ref_10, Q_hc_ref_60) = tc.load_res("inputs/case07_res.csv")
Q_hc_ref_1 = Q_hc_ref_1[:,0]
Q_hc_ref_10 = Q_hc_ref_10[:,0]
Q_hc_ref_60 = Q_hc_ref_60[:,0]
# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
plt.figure()
ax_top = plt.subplot(211)
plt.plot(ref, label="Reference", color="black", linestyle="--")
plt.plot(res, label="Simulation", color="blue", linestyle="-")
plt.legend()
plt.ylabel("Heat load in W")
plt.title(title)
windowIndoorSurface,
extWallIndoorSurface,
intWallIndoorSurface,
intGainsConv, intGainsRad, ports,
model,
dt=int(3600/times_per_hour))
# Compute averaged results
Q_hc_mean = np.array([np.mean(Q_HC[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])
Q_hc_1 = Q_hc_mean[0:24]
Q_hc_10 = Q_hc_mean[216:240]
Q_hc_60 = Q_hc_mean[1416:1440]
# Load reference results
(Q_hc_ref_1, Q_hc_ref_10, Q_hc_ref_60) = tc.load_res("inputs/case06_res.csv")
Q_hc_ref_1 = -Q_hc_ref_1[:,0]
Q_hc_ref_10 = -Q_hc_ref_10[:,0]
Q_hc_ref_60 = -Q_hc_ref_60[:,0]
# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
plt.figure()
ax_top = plt.subplot(211)
plt.plot(ref, label="Reference", color="black", linestyle="--")
plt.plot(res, label="Simulation", color="blue", linestyle="-")
plt.legend()
plt.ylabel("Heat load in W")
plt.title(title)
Q_iw_mean = np.array([np.mean(Q_iw[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])
Q_ow_mean = np.array([np.mean(Q_ow[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])
Q_hc_1 = Q_hc_mean[0:24] + Q_iw_mean[0:24] + Q_ow_mean[0:24]
Q_hc_10 = Q_hc_mean[216:240] + Q_iw_mean[216:240] + Q_ow_mean[216:240]
Q_hc_60 = Q_hc_mean[1416:1440] + Q_iw_mean[1416:1440] + Q_ow_mean[1416:1440]
T_air_c = T_air - 273.15
T_air_mean = np.array([np.mean(T_air_c[i*times_per_hour:(i+1)*times_per_hour]) for i in range(24*60)])
T_air_1 = T_air_mean[0:24]
T_air_10 = T_air_mean[216:240]
T_air_60 = T_air_mean[1416:1440]
# Load reference results
(load_res_1, load_res_10, load_res_60) = tc.load_res("inputs/case11_res.csv")
Q_hc_ref_1 = load_res_1[:,1]
Q_hc_ref_10 = load_res_10[:,1]
Q_hc_ref_60 = load_res_60[:,1]
T_air_ref_1 = load_res_1[:,0]
T_air_ref_10 = load_res_10[:,0]
T_air_ref_60 = load_res_60[:,0]
# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
plt.figure()
ax_top = plt.subplot(211)
plt.plot(ref, label="Reference", color="black", linestyle="--")
plt.plot(res, label="Simulation", color="blue", linestyle="-")
Q_hc_1 = Q_hc_mean[0:24] + Q_iw_mean[0:24] + Q_ow_mean[0:24]
Q_hc_10 = Q_hc_mean[216:240] + Q_iw_mean[216:240] + Q_ow_mean[216:240]
Q_hc_60 = Q_hc_mean[1416:1440] + Q_iw_mean[1416:1440] + Q_ow_mean[1416:1440]
T_air_c = T_air - 273.15
T_air_mean = np.array([
np.mean(T_air_c[i * times_per_hour:(i + 1) * times_per_hour])
for i in range(24 * 60)
])
T_air_1 = T_air_mean[0:24]
T_air_10 = T_air_mean[216:240]
T_air_60 = T_air_mean[1416:1440]
# Load reference results
(load_res_1, load_res_10, load_res_60) = tc.load_res("inputs/case11_res.csv")
Q_hc_ref_1 = load_res_1[:, 1]
Q_hc_ref_10 = load_res_10[:, 1]
Q_hc_ref_60 = load_res_60[:, 1]
T_air_ref_1 = load_res_1[:, 0]
T_air_ref_10 = load_res_10[:, 0]
T_air_ref_60 = load_res_60[:, 0]
# Plot comparisons
def plot_result(res, ref, title="Results day 1"):
plt.figure()
ax_top = plt.subplot(211)
plt.plot(ref, label="Reference", color="black", linestyle="--")
plt.plot(res, label="Simulation", color="blue", linestyle="-")