def calc_cool_loads_3for2(bpr, t, tsd): """ Calculation procedure for cooling system loads of AHU, ARU and SCU subsystems of 3for2 system Gabriel Happle, Feb. 2018 :param bpr: building properties row object :param t: time step / hour of year [0..8760] :param tsd: time series data dict :return: dict of rc_model_temperatures """ # get values from tsd m_ve_mech = tsd['m_ve_mech'][t] t_ve_mech_after_hex = tsd['theta_ve_mech'][t] x_ve_mech = tsd['x_ve_mech'][t] t_int_prev = tsd['T_int'][t - 1] x_int_prev = tsd['x_int'][t - 1] # *** # RC MODEL # *** # calculate rc model demand qc_sen_rc_demand, rc_model_temperatures = calc_rc_cooling_demand(bpr=bpr, tsd=tsd, t=t) # *** # AHU # *** # calculate ahu loads ahu_loads = airconditioning_model.central_air_handling_unit_cooling( m_ve_mech, t_ve_mech_after_hex, x_ve_mech, bpr) qc_sen_ahu = ahu_loads['qc_sen_ahu'] qc_lat_ahu = ahu_loads['qc_lat_ahu'] tsd['x_ve_mech'][t] = ahu_loads['x_sup_c_ahu'] # *** # ARU # *** # calculate recirculation unit dehumidification demand tsd['T_int'][t] = rc_model_temperatures[ 'T_int'] # dehumidification load needs zone temperature # NOTE: here we might make some error, as we calculate the moisture set point for the # uncorrected zone air temperature (i.e. no over cooling) g_dhu_demand_aru = latent_loads.calc_dehumidification_moisture_load( bpr, tsd, t) # no sensible demand that controls the ARU qc_sen_demand_aru = 0.0 # calculate ARU system loads with T and x control activated aru_system_loads = airconditioning_model.local_air_recirculation_unit_cooling( qc_sen_demand_aru, g_dhu_demand_aru, t_int_prev, x_int_prev, bpr, t_control=False, x_control=True) g_dhu_aru = aru_system_loads['g_dhu_aru'] qc_lat_aru = aru_system_loads['qc_lat_aru'] qc_sen_aru = aru_system_loads['qc_sen_aru'] # *** # SCU # *** # calculate remaining sensible cooling demand to be met by radiative cooling qc_sen_demand_scu = np.min( [0.0, qc_sen_rc_demand - qc_sen_ahu - qc_sen_aru]) # demand is load qc_sen_scu = qc_sen_demand_scu # *** # ADJUST RC MODEL TEMPERATURE # *** # TODO: check, if it is smaller something went wrong in the calculation qc_sen_total = qc_sen_ahu + qc_sen_aru + qc_sen_scu # update rc model temperatures rc_model_temperatures = rc_model_SIA.calc_rc_model_temperatures_cooling( qc_sen_total, bpr, tsd, t) # *** # ZONE MOISTURE # *** # action on moisture tsd['g_hu_ld'][t] = 0.0 # no humidification tsd['g_dhu_ld'][t] = g_dhu_aru latent_loads.calc_moisture_content_in_zone_local(bpr, tsd, t) # *** # emission losses # *** # emission losses on total sensible load # TODO: check # write to tsd tsd['Qcs_sen_rc'][t] = qc_sen_rc_demand tsd['Qcs_sen_ahu'][t] = qc_sen_ahu tsd['Qcs_sen_aru'][t] = qc_sen_aru tsd['Qcs_sen_scu'][t] = qc_sen_scu tsd['Qcs_lat_ahu'][t] = qc_lat_ahu tsd['Qcs_lat_aru'][t] = qc_lat_aru rc_temperatures_to_tsd(rc_model_temperatures, tsd, t) tsd['Qcs_sen_sys'][ t] = qc_sen_ahu + qc_sen_aru + qc_sen_scu # sum system loads tsd['Qcs_lat_sys'][t] = qc_lat_ahu + qc_lat_aru # air flow tsd['m_ve_rec'][t] = aru_system_loads['ma_sup_cs_aru'] # mass flows to tsd tsd['ma_sup_cs_ahu'][t] = ahu_loads['ma_sup_cs_ahu'] tsd['ta_sup_cs_ahu'][t] = ahu_loads['ta_sup_cs_ahu'] tsd['ta_re_cs_ahu'][t] = ahu_loads['ta_re_cs_ahu'] tsd['ma_sup_cs_aru'][t] = aru_system_loads['ma_sup_cs_aru'] tsd['ta_sup_cs_aru'][t] = aru_system_loads['ta_sup_cs_aru'] tsd['ta_re_cs_aru'][t] = aru_system_loads['ta_re_cs_aru'] q_em_ls_cooling = space_emission_systems.calc_q_em_ls_cooling(bpr, tsd, t) tsd['Qcs_em_ls'][t] = q_em_ls_cooling # system status tsd['sys_status_ahu'][t] = 'On' tsd['sys_status_aru'][t] = 'On:R' tsd['sys_status_sen'][t] = 'On' # the return is only for the input into the detailed thermal reverse calculations for the dashboard graphs return rc_model_temperatures
def calc_rc_cooling_demand(bpr, tsd, t): """ Crank-Nicholson Procedure to calculate heating / cooling demand of buildings following the procedure in 2.3.2 in SIA 2044 / Korrigenda C1 zum Merkblatt SIA 2044:2011 / Korrigenda C2 zum Mekblatt SIA 2044:2011 Special procedures for updating ventilation air AC-heated and AC-cooled buildings Author: Gabriel Happle Date: 01/2017 :param bpr: building properties row object :param tsd: time series data dict :param t: time step / hour of year [0..8760] :return: phi_c_act, rc_model_temperatures """ # following the procedure in 2.3.2 in SIA 2044 / Korrigenda C1 zum Merkblatt SIA 2044:2011 # / Korrigenda C2 zum Mekblatt SIA 2044:2011 # ++++++++++++++++ # CASE 2 - COOLING # ++++++++++++++++ # case for cooling # tsd['system_status'][t] = 'Radiative cooling' # STEP 1 # ****** # calculate temperatures with 0 heating power rc_model_temperatures_0 = rc_model_SIA.calc_rc_model_temperatures_no_heating_cooling( bpr, tsd, t) t_int_0 = rc_model_temperatures_0['T_int'] # CHECK FOR DEMAND if not rc_model_SIA.has_sensible_cooling_demand(t_int_0, tsd, t): # return zero demand rc_model_temperatures = rc_model_temperatures_0 phi_c_act = 0.0 elif rc_model_SIA.has_sensible_cooling_demand(t_int_0, tsd, t): # continue # STEP 2 # ****** # calculate temperatures with 10 W/m2 cooling power phi_hc_10 = 10.0 * bpr.rc_model['Af'] rc_model_temperatures_10 = rc_model_SIA.calc_rc_model_temperatures_cooling( phi_hc_10, bpr, tsd, t) t_int_10 = rc_model_temperatures_10['T_int'] t_int_set = tsd['ta_cs_set'][t] # interpolate heating power # (64) in SIA 2044 / Korrigenda C1 zum Merkblatt SIA 2044:2011 / Korrigenda C2 zum Mekblatt SIA 2044:2011 phi_hc_ul = phi_hc_10 * (t_int_set - t_int_0) / (t_int_10 - t_int_0) # STEP 3 # ****** # check if available power is sufficient phi_c_max = -bpr.hvac['Qcsmax_Wm2'] * bpr.rc_model['Af'] if 0.0 > phi_hc_ul >= phi_c_max: # case heating with phi_hc_ul # calculate temperatures with this power phi_c_act = phi_hc_ul elif 0.0 > phi_hc_ul < phi_c_max: # case heating with max power available # calculate temperatures with this power phi_c_act = phi_c_max else: raise Exception("Unexpected status in 'calc_rc_cooling_demand'") # STEP 4 # ****** rc_model_temperatures = rc_model_SIA.calc_rc_model_temperatures_cooling( phi_c_act, bpr, tsd, t) else: raise Exception("Unexpected status in 'calc_rc_cooling_demand'") return phi_c_act, rc_model_temperatures
def calc_rc_model_demand_heating_cooling(bpr, tsd, t, gv): """ Crank-Nicholson Procedure to calculate heating / cooling demand of buildings following the procedure in 2.3.2 in SIA 2044 / Korrigenda C1 zum Merkblatt SIA 2044:2011 / Korrigenda C2 zum Mekblatt SIA 2044:2011 Special procedures for updating ventilation air AC-heated and AC-cooled buildings Author: Gabriel Happle Date: 01/2017 :param bpr: building properties row object :param tsd: time series data dict :param t: time step / hour of year [0..8760] :param gv: globalvars :return: updates values in tsd """ # following the procedure in 2.3.2 in SIA 2044 / Korrigenda C1 zum Merkblatt SIA 2044:2011 # / Korrigenda C2 zum Mekblatt SIA 2044:2011 # ++++++++++++++++++++++++++++++ # CASE 0 - NO HEATING OR COOLING # ++++++++++++++++++++++++++++++ if not control_heating_cooling_systems.is_active_heating_system(bpr, tsd, t) \ and not control_heating_cooling_systems.is_active_cooling_system(bpr, tsd, t): # STEP 1 # ****** # calculate temperatures rc_model_temperatures = rc_model_SIA.calc_rc_model_temperatures_no_heating_cooling( bpr, tsd, t) # write to tsd tsd['theta_a'][t] = rc_model_temperatures['theta_a'] tsd['theta_m'][t] = rc_model_temperatures['theta_m'] tsd['theta_c'][t] = rc_model_temperatures['theta_c'] tsd['theta_o'][t] = rc_model_temperatures['theta_o'] update_tsd_no_cooling(tsd, t) update_tsd_no_heating(tsd, t) tsd['system_status'][t] = 'systems off' # ++++++++++++++++ # CASE 1 - HEATING # ++++++++++++++++ elif control_heating_cooling_systems.is_active_heating_system(bpr, tsd, t): # case for heating tsd['system_status'][t] = 'Radiative heating' # STEP 1 # ****** # calculate temperatures with 0 heating power rc_model_temperatures_0 = rc_model_SIA.calc_rc_model_temperatures_no_heating_cooling( bpr, tsd, t) theta_a_0 = rc_model_temperatures_0['theta_a'] # STEP 2 # ****** # calculate temperatures with 10 W/m2 heating power phi_hc_10 = 10 * bpr.rc_model['Af'] rc_model_temperatures_10 = rc_model_SIA.calc_rc_model_temperatures_heating( phi_hc_10, bpr, tsd, t) theta_a_10 = rc_model_temperatures_10['theta_a'] theta_a_set = tsd['ta_hs_set'][t] # interpolate heating power # (64) in SIA 2044 / Korrigenda C1 zum Merkblatt SIA 2044:2011 / Korrigenda C2 zum Mekblatt SIA 2044:2011 phi_hc_ul = phi_hc_10 * (theta_a_set - theta_a_0) / (theta_a_10 - theta_a_0) # STEP 3 # ****** # check if available power is sufficient phi_h_max = bpr.hvac['Qhsmax_Wm2'] * bpr.rc_model['Af'] if 0 < phi_hc_ul <= phi_h_max: # case heating with phi_hc_ul # calculate temperatures with this power phi_h_act = phi_hc_ul elif 0 < phi_hc_ul > phi_h_max: # case heating with max power available # calculate temperatures with this power phi_h_act = phi_h_max else: raise # STEP 4 # ****** rc_model_temperatures = rc_model_SIA.calc_rc_model_temperatures_heating( phi_h_act, bpr, tsd, t) # write necessary parameters for AC calculation to tsd tsd['theta_a'][t] = rc_model_temperatures['theta_a'] tsd['theta_m'][t] = rc_model_temperatures['theta_m'] tsd['theta_c'][t] = rc_model_temperatures['theta_c'] tsd['theta_o'][t] = rc_model_temperatures['theta_o'] tsd['Qhs_sen'][t] = phi_h_act tsd['Qhs_sen_sys'][t] = phi_h_act tsd['Qhs_lat_sys'][t] = 0 tsd['Ehs_lat_aux'][t] = 0 tsd['ma_sup_hs'][t] = 0 tsd['Ta_sup_hs'][t] = 0 tsd['Ta_re_hs'][t] = 0 tsd['m_ve_recirculation'][t] = 0 # STEP 5 - latent and sensible heat demand of AC systems # ****** if control_heating_cooling_systems.heating_system_is_ac(bpr): air_con_model_loads_flows_temperatures = airconditioning_model.calc_hvac_heating( tsd, t, gv) tsd['system_status'][t] = 'AC heating' # update temperatures for over heating case if air_con_model_loads_flows_temperatures[ 'q_hs_sen_hvac'] > phi_h_act: phi_h_act_over_heating = air_con_model_loads_flows_temperatures[ 'q_hs_sen_hvac'] rc_model_temperatures = rc_model_SIA.calc_rc_model_temperatures_heating( phi_h_act_over_heating, bpr, tsd, t) # update temperatures tsd['theta_a'][t] = rc_model_temperatures['theta_a'] tsd['theta_m'][t] = rc_model_temperatures['theta_m'] tsd['theta_c'][t] = rc_model_temperatures['theta_c'] tsd['theta_o'][t] = rc_model_temperatures['theta_o'] tsd['system_status'][t] = 'AC over heating' # update AC energy demand tsd['Qhs_sen_sys'][t] = air_con_model_loads_flows_temperatures[ 'q_hs_sen_hvac'] tsd['Qhs_lat_sys'][t] = air_con_model_loads_flows_temperatures[ 'q_hs_lat_hvac'] tsd['ma_sup_hs'][t] = air_con_model_loads_flows_temperatures[ 'ma_sup_hs'] tsd['Ta_sup_hs'][t] = air_con_model_loads_flows_temperatures[ 'ta_sup_hs'] tsd['Ta_re_hs'][t] = air_con_model_loads_flows_temperatures[ 'ta_re_hs'] tsd['Ehs_lat_aux'][t] = air_con_model_loads_flows_temperatures[ 'e_hs_lat_aux'] tsd['m_ve_recirculation'][ t] = air_con_model_loads_flows_temperatures[ 'm_ve_hvac_recirculation'] # STEP 6 - emission system losses # ****** q_em_ls_heating = space_emission_systems.calc_q_em_ls_heating( bpr, tsd, t) # set temperatures to tsd for heating tsd['theta_a'][t] = rc_model_temperatures['theta_a'] tsd['theta_m'][t] = rc_model_temperatures['theta_m'] tsd['theta_c'][t] = rc_model_temperatures['theta_c'] tsd['theta_o'][t] = rc_model_temperatures['theta_o'] tsd['Qhs_lat_sys'][t] = 0 tsd['Qhs_em_ls'][t] = q_em_ls_heating tsd['Qhs_sen'][t] = phi_h_act tsd['Qhsf'][t] = 0 tsd['Qhsf_lat'][t] = 0 update_tsd_no_cooling(tsd, t) # ++++++++++++++++ # CASE 2 - COOLING # ++++++++++++++++ elif control_heating_cooling_systems.is_active_cooling_system(bpr, tsd, t): # case for cooling tsd['system_status'][t] = 'Radiative cooling' # STEP 1 # ****** # calculate temperatures with 0 heating power rc_model_temperatures_0 = rc_model_SIA.calc_rc_model_temperatures_no_heating_cooling( bpr, tsd, t) theta_a_0 = rc_model_temperatures_0['theta_a'] # STEP 2 # ****** # calculate temperatures with 10 W/m2 cooling power phi_hc_10 = 10 * bpr.rc_model['Af'] rc_model_temperatures_10 = rc_model_SIA.calc_rc_model_temperatures_cooling( phi_hc_10, bpr, tsd, t) theta_a_10 = rc_model_temperatures_10['theta_a'] theta_a_set = tsd['ta_cs_set'][t] # interpolate heating power # (64) in SIA 2044 / Korrigenda C1 zum Merkblatt SIA 2044:2011 / Korrigenda C2 zum Mekblatt SIA 2044:2011 phi_hc_ul = phi_hc_10 * (theta_a_set - theta_a_0) / (theta_a_10 - theta_a_0) # STEP 3 # ****** # check if available power is sufficient phi_c_max = -bpr.hvac['Qcsmax_Wm2'] * bpr.rc_model['Af'] if 0 > phi_hc_ul >= phi_c_max: # case heating with phi_hc_ul # calculate temperatures with this power phi_c_act = phi_hc_ul elif 0 > phi_hc_ul < phi_c_max: # case heating with max power available # calculate temperatures with this power phi_c_act = phi_c_max else: raise # STEP 4 # ****** rc_model_temperatures = rc_model_SIA.calc_rc_model_temperatures_cooling( phi_c_act, bpr, tsd, t) # write necessary parameters for AC calculation to tsd tsd['theta_a'][t] = rc_model_temperatures['theta_a'] tsd['theta_m'][t] = rc_model_temperatures['theta_m'] tsd['theta_c'][t] = rc_model_temperatures['theta_c'] tsd['theta_o'][t] = rc_model_temperatures['theta_o'] tsd['Qcs_sen'][t] = phi_c_act tsd['Qcs_sen_sys'][t] = phi_c_act tsd['Qcs_lat_sys'][t] = 0 tsd['ma_sup_cs'][t] = 0 tsd['m_ve_recirculation'][t] = 0 # STEP 5 - latent and sensible heat demand of AC systems # ****** if control_heating_cooling_systems.cooling_system_is_ac(bpr): tsd['system_status'][t] = 'AC cooling' air_con_model_loads_flows_temperatures = airconditioning_model.calc_hvac_cooling( tsd, t, gv) # update temperatures for over cooling case if air_con_model_loads_flows_temperatures[ 'q_cs_sen_hvac'] < phi_c_act: phi_c_act_over_cooling = air_con_model_loads_flows_temperatures[ 'q_cs_sen_hvac'] rc_model_temperatures = rc_model_SIA.calc_rc_model_temperatures_cooling( phi_c_act_over_cooling, bpr, tsd, t) # update temperatures tsd['theta_a'][t] = rc_model_temperatures['theta_a'] tsd['theta_m'][t] = rc_model_temperatures['theta_m'] tsd['theta_c'][t] = rc_model_temperatures['theta_c'] tsd['theta_o'][t] = rc_model_temperatures['theta_o'] tsd['system_status'][t] = 'AC over cooling' # update AC energy demand tsd['Qcs_sen_sys'][t] = air_con_model_loads_flows_temperatures[ 'q_cs_sen_hvac'] tsd['Qcs_lat_sys'][t] = air_con_model_loads_flows_temperatures[ 'q_cs_lat_hvac'] tsd['ma_sup_cs'][t] = air_con_model_loads_flows_temperatures[ 'ma_sup_cs'] tsd['Ta_sup_cs'][t] = air_con_model_loads_flows_temperatures[ 'ta_sup_cs'] tsd['Ta_re_cs'][t] = air_con_model_loads_flows_temperatures[ 'ta_re_cs'] tsd['m_ve_recirculation'][ t] = air_con_model_loads_flows_temperatures[ 'm_ve_hvac_recirculation'] # STEP 6 - emission system losses # ****** q_em_ls_cooling = space_emission_systems.calc_q_em_ls_cooling( bpr, tsd, t) # set temperatures to tsd for heating tsd['theta_a'][t] = rc_model_temperatures['theta_a'] tsd['theta_m'][t] = rc_model_temperatures['theta_m'] tsd['theta_c'][t] = rc_model_temperatures['theta_c'] tsd['theta_o'][t] = rc_model_temperatures['theta_o'] tsd['Qcs'][t] = 0 tsd['Qcs_em_ls'][t] = q_em_ls_cooling tsd['Qcsf'][t] = 0 tsd['Qcsf_lat'][t] = 0 update_tsd_no_heating(tsd, t) return