def main():
house_param = load_config("config2R2C.yml")
days_sim = house_param['timing']['days_sim']
CF = house_param['ventilation']['CF']
Rair_wall, Cwall, Rair_outdoor, Cair = calculateRCOne(house_param)
print(days_sim)
#Loading the radiator and buffervessel parameters
#Heat transfer coefficient of the radiator and het capacity
cpwater = house_param['radiator']['cpwater']
rhowater = house_param['radiator']['rhowater']
Urad = house_param['radiator']['Urad']
Arad = house_param['radiator']['Arad']
volumeRadiator = house_param['radiator']['volume_rad']
UAradiator = Urad * Arad
Crad = cpwater * volumeRadiator * rhowater
#Heat capacity of the buffervessel
volumeBuffervessel = house_param['radiator']['volume_buffervessel']
Cbuffervessel = cpwater * volumeBuffervessel * rhowater
df_nen = nen5060_to_dataframe()
df_irr = run_qsun(df_nen)
print(df_irr.head())
time_sim = df_irr.iloc[0:days_sim * 24, 0].values
Qsolar = (df_irr.total_E * house_param['glass']['E'] +
df_irr.total_SE * house_param['glass']['SE'] +
df_irr.total_S * house_param['glass']['S'] +
df_irr.total_SW * house_param['glass']['SW'] +
df_irr.total_W * house_param['glass']['W'] +
df_irr.total_NW * house_param['glass']['NW'] +
df_irr.total_N * house_param['glass']['N'] +
df_irr.total_NE * house_param['glass']['NE']).values
Qsolar *= house_param['glass']['g_value']
Qsolar_sim = Qsolar[0:days_sim * 24]
Qint = internal_heat_gain(house_param['internal']['Q_day'],
house_param['internal']['delta_Q'],
house_param['internal']['t1'],
house_param['internal']['t2'])
Qinternal_sim = Qint[0:days_sim * 24]
Toutdoor = df_nen.loc[:, 'temperatuur'].values / 10.0 # temperature
T_outdoor_sim = Toutdoor[0:days_sim * 24]
week_day_setpoint = thermostat_sp(
house_param['setpoint']['t1'], house_param['setpoint']['t2'],
house_param['setpoint']['Night_T_SP'],
house_param['setpoint']['Day_T_SP'],
house_param['setpoint']['Flex_T_SP_workday'],
house_param['setpoint']['Wu_time'],
house_param['setpoint']['Work_time'],
house_param['setpoint']['back_home_from_work'])
day_off_setpoint = thermostat_sp(
house_param['setpoint']['t1'], house_param['setpoint']['t2'],
house_param['setpoint']['Night_T_SP'],
house_param['setpoint']['Day_T_SP'],
house_param['setpoint']['Flex_T_SP_dayoff'],
house_param['setpoint']['Wu_time'],
house_param['setpoint']['shopping_time'],
house_param['setpoint']['back_home'])
SP = SP_profile(week_day_setpoint, day_off_setpoint)
SP_sim = SP[0:days_sim * 24]
# solve ODE
data = house_buffervessel(T_outdoor_sim, Qinternal_sim, Qsolar_sim, SP_sim,
time_sim, CF, Rair_outdoor, Rair_wall, Cair,
Cwall, UAradiator, Crad, Cbuffervessel, cpwater)
# plot the results
plt.figure(figsize=(15, 5)) # key-value pair: no spaces
plt.plot(data[4], data[0], label='Tair')
plt.plot(data[4], data[1], label='Twall')
plt.plot(data[4], data[2], label='Treturn')
plt.plot(data[4], data[3], label='Tbuffervessel')
plt.plot(time_sim, SP_sim, label='SP_Temperature')
plt.plot(time_sim, T_outdoor_sim, label='Toutdoor')
plt.legend(loc='best')
plt.show()
def main():
house_param = load_config("config2R2C.yml")
days_sim = house_param['timing']['days_sim']
CF = house_param['ventilation']['CF']
Rair_wall, Cwall, Rair_outdoor, Cair = calculateRCOne(house_param)
print('Simulation days:', days_sim)
df_nen = nen5060_to_dataframe()
df_irr = run_qsun(df_nen)
#df_weeks = read_week('NEN_data')
print(df_irr.head())
time_sim = df_irr.iloc[0:days_sim * 24, 0].values
Qsolar = (df_irr.total_E * house_param['glass']['E'] +
df_irr.total_SE * house_param['glass']['SE'] +
df_irr.total_S * house_param['glass']['S'] +
df_irr.total_SW * house_param['glass']['SW'] +
df_irr.total_W * house_param['glass']['W'] +
df_irr.total_NW * house_param['glass']['NW'] +
df_irr.total_N * house_param['glass']['N'] +
df_irr.total_NE * house_param['glass']['NE']).values
Qsolar *= house_param['glass']['g_value']
Qsolar_sim = Qsolar[0:days_sim * 24]
#print(len(Qsolar_sim))
Qint = internal_heat_gain(house_param['internal']['Q_day'],
house_param['internal']['delta_Q'],
house_param['internal']['t1'],
house_param['internal']['t2'])
Qinternal_sim = Qint[0:days_sim * 24]
Toutdoor = df_nen.loc[:, 'temperatuur'].values / 10.0 # temperature
T_outdoor_sim = Toutdoor[0:days_sim * 24]
#plt.plot(T_outdoor_sim)
week_day_setpoint = thermostat_sp(
house_param['setpoint']['t1'], house_param['setpoint']['t2'],
house_param['setpoint']['Night_T_SP'],
house_param['setpoint']['Day_T_SP'],
house_param['setpoint']['Flex_T_SP_workday'],
house_param['setpoint']['Wu_time'],
house_param['setpoint']['Work_time'],
house_param['setpoint']['back_home_from_work'])
day_off_setpoint = thermostat_sp(
house_param['setpoint']['t1'], house_param['setpoint']['t2'],
house_param['setpoint']['Night_T_SP'],
house_param['setpoint']['Day_T_SP'],
house_param['setpoint']['Flex_T_SP_dayoff'],
house_param['setpoint']['Wu_time'],
house_param['setpoint']['shopping_time'],
house_param['setpoint']['back_home'])
SP = SP_profile(week_day_setpoint, day_off_setpoint)
#SP = temp_sp(house_param['setpoint']['t1'],
# house_param['setpoint']['t2'],
# house_param['setpoint']['Night_T_SP'],
# house_param['setpoint']['Day_T_SP'],
# house_param['setpoint']['Wu_time'],
# house_param['setpoint']['Work_time'],
# house_param['setpoint']['back_home'])
SP_sim = SP[0:days_sim * 24]
# addition NTA8800 house model
# Controller value
kp = house_param['controller']['kp']
# solve ODE
data = house(T_outdoor_sim, Qinternal_sim, Qsolar_sim, SP_sim, time_sim,
CF, Rair_outdoor, Rair_wall, Cair, Cwall, kp)
# plot the results
plt.figure(figsize=(15, 5)) # key-value pair: no spaces
plt.plot(data[0], label='Tair')
plt.plot(data[1], label='Twall')
plt.plot(SP_sim, label='SP_Temperature')
# plt.plot(T_outdoor_sim,label='Toutdoor')
plt.legend(loc='best')
plt.title("Simulation2R2C")
plt.show()
'''
def main():
house_param = load_config("config2R2C.yml")
days_sim = house_param['timing']['days_sim']
CF = house_param['ventilation']['CF']
Rair_wall, Cwall, Rair_outdoor, Cair = calculateRCOne(house_param)
print(days_sim)
df_nen = nen5060_to_dataframe()
df_irr = run_qsun(df_nen)
print(df_irr.head())
time_sim = df_irr.iloc[0:days_sim * 24, 0].values
Qsolar = (df_irr.total_E * house_param['glass']['E'] +
df_irr.total_SE * house_param['glass']['SE'] +
df_irr.total_S * house_param['glass']['S'] +
df_irr.total_SW * house_param['glass']['SW'] +
df_irr.total_W * house_param['glass']['W'] +
df_irr.total_NW * house_param['glass']['NW'] +
df_irr.total_N * house_param['glass']['N'] +
df_irr.total_NE * house_param['glass']['NE']).values
Qsolar *= house_param['glass']['g_value']
Qsolar_sim = Qsolar[0:days_sim * 24]
Qint = internal_heat_gain(house_param['internal']['Q_day'],
house_param['internal']['delta_Q'],
house_param['internal']['t1'],
house_param['internal']['t2'])
Qinternal_sim = Qint[0:days_sim * 24]
Toutdoor = df_nen.loc[:, 'temperatuur'].values / 10.0 # temperature
T_outdoor_sim = Toutdoor[0:days_sim * 24]
week_day_setpoint = thermostat_sp(house_param['setpoint']['t1'],
house_param['setpoint']['t2'],
house_param['setpoint']['Night_T_SP'],
house_param['setpoint']['Day_T_SP'],
house_param['setpoint']['Flex_T_SP_workday'],
house_param['setpoint']['Wu_time'],
house_param['setpoint']['Work_time'],
house_param['setpoint']['back_home_from_work'])
day_off_setpoint = thermostat_sp(house_param['setpoint']['t1'],
house_param['setpoint']['t2'],
house_param['setpoint']['Night_T_SP'],
house_param['setpoint']['Day_T_SP'],
house_param['setpoint']['Flex_T_SP_dayoff'],
house_param['setpoint']['Wu_time'],
house_param['setpoint']['shopping_time'],
house_param['setpoint']['back_home'])
SP = SP_profile(week_day_setpoint, day_off_setpoint)
SP_sim = SP[0:days_sim * 24]
g = GasBoiler(kp=5000, ki=0, kd=0, T_setpoint=19.9, T_node=15, T_amb=10, dead_band=1, P_max=10000, P_min=2000)
# solve ODE
data = house_buffervessel(T_outdoor_sim, Qinternal_sim, Qsolar_sim, SP_sim, time_sim,
CF, Rair_outdoor, Rair_wall, Cair, Cwall, g)
# plot the results
plt.figure(figsize=(15, 5)) # key-value pair: no spaces
plt.plot(data[3], data[0], label='Tair')
plt.plot(data[3], data[1], label='Twall')
# plt.plot(time_sim, SP_sim, label='SP_Temperature')
# plt.plot(time_sim,T_outdoor_sim,label='Toutdoor')
plt.legend(loc='best')
plt.show()
from housemodel.sourcesink.NEN5060 import nen5060_to_dataframe, run_qsun, run_qsun_new
def my_assert(condition, fail_str, suc_str):
assert condition, fail_str
print(suc_str)
df_nen = nen5060_to_dataframe()
Qsolar = np.zeros((8760, 1)) # 8760 rows x 2 cols
az = 0
tlt = 0
logging.info(f"Azimuth {az} and Tilt {tlt}")
df_irr_old = run_qsun(df_nen)
# north_is_zero (niz), check azimuth 0(N), 90(E), 180(S), 270(W)
# in qsun, azimuth = azimuth - 180: north = -180 internally (works!)
df_irr_north_niz = run_qsun_new(df_nen, 0, 90, north_is_zero=True)
Qsolar_north_niz = df_irr_north_niz.total_irr.values
df_irr_east_niz = run_qsun_new(df_nen, 90, 90, north_is_zero=True)
Qsolar_east_niz = (df_irr_east_niz.total_irr).values
df_irr_south_niz = run_qsun_new(df_nen, 180, 90, north_is_zero=True)
Qsolar_south_niz = (df_irr_south_niz.total_irr).values
df_irr_west_niz = run_qsun_new(df_nen, 270, 90, north_is_zero=True)
Qsolar_west_niz = (df_irr_west_niz.total_irr).values