Method heatingAndCooling

src/isomodel/SimModel.cpp:1484–1634

  */
  }
  void SimModel::heatingAndCooling(const Vector& v_E_sol, const Vector& v_Th_avg, const Vector& v_Hve_ht, const Vector& v_Tc_avg,
                                   const Vector& v_Hve_cl, double tau, double H_tr, double phi_I_tot, double frac_hrs_wk_day, Vector& v_Qfan_tot,
                                   Vector& v_Qneed_ht, Vector& v_Qneed_cl, double& Qneed_ht_yr, double& Qneed_cl_yr) const {
    Vector temp = mult(megasecondsInMonth, phi_I_tot, 12);
    Vector v_tot_mo_ht_gain = sum(temp, v_E_sol);

    double a_H0 = 1;
    double tau_H0 = 15;
    double a_H = a_H0 + tau / tau_H0;

    Vector v_QT_ht = mult(mult(dif(v_Th_avg, location->weather()->mdbt()), megasecondsInMonth), H_tr);
    Vector v_QV_ht = mult(mult(mult(v_Hve_ht, structure->floorArea()), dif(v_Th_avg, location->weather()->mdbt())), megasecondsInMonth);
    Vector v_Qtot_ht = sum(v_QT_ht, v_QV_ht);
    /*
  %% Heating and Cooling Needs

%total monthly heat gains (MJ)

v_tot_mo_ht_gain = phi_I_tot*v_Msec_ina_mo + v_E_sol;  % total_heat_gain = total internal + total solar in MJ/m2.

% compute the heating need including thermal mass effects
% NOTE: the building heat thermal time constant, tau, was calculated in the section
% on interior temperature
a_H0=1; % a_H_0 = reference dimensionless parameter
tau_H0=15; % tau_H_0 = reference time constant
a_H = a_H0+tau/tau_H0;  %a_H_building heating dimensionless constant

v_QT_ht = H_tr.*(v_Th_avg-v_mdbt).*v_Msec_ina_mo;  %QT = transmission loss (MJ)
v_QV_ht = v_Hve_ht*In.cond_flr_area.*(v_Th_avg-v_mdbt).*v_Msec_ina_mo; % QV in MJ
v_Qtot_ht = v_QT_ht+v_QV_ht ; %QL_total total heat loss in MJ
*/
    Vector v_gamma_H_ht = div(v_tot_mo_ht_gain, sum(v_Qtot_ht, std::numeric_limits<double>::min()));
    Vector v_eta_g_H(12);
    for (size_t i = 0; i < v_eta_g_H.size(); i++) {
      v_eta_g_H[i] = v_gamma_H_ht(i) > 0 ? (1 - std::pow(v_gamma_H_ht[i], a_H)) / (1 - std::pow(v_gamma_H_ht[i], (a_H + 1)))
                                         : 1 / (v_gamma_H_ht(i) + std::numeric_limits<double>::min());
    }
    v_Qneed_ht = dif(v_Qtot_ht, mult(v_eta_g_H, v_tot_mo_ht_gain));
    Qneed_ht_yr = sum(v_Qneed_ht);

    /*
% compute the ratio of heat gain to heating loss, gamma_H
v_gamma_H_ht = v_tot_mo_ht_gain./(v_Qtot_ht+eps);  %gamma_H = QG/QL  - eps added to avoid divide by zero problem

% for each month, check the heat gain ratio and set the heating gain
% utilization factor, eta_g_H accordingly
v_eta_g_H=zeros(12,1);
for I=1:12
    if v_gamma_H_ht(I)>0
        v_eta_g_H(I) = (1-v_gamma_H_ht(I).^a_H)./(1-v_gamma_H_ht(I).^(a_H+1));  % eta_g_H heating gain utilization factor
    else
        v_eta_g_H(I) = 1./(v_gamma_H_ht(I)+eps);
    end
end
v_Qneed_ht = v_Qtot_ht - v_eta_g_H.*v_tot_mo_ht_gain; %QNH = QL,H - eta_G_H.*Q_G_H

Qneed_ht_yr = sum(v_Qneed_ht);
   */