openstudio-standards 0.3.0 → 0.3.1.rc2

data/lib/openstudio-standards/standards/necb/NECB2020/service_water_heating.rb

@@ -7,7 +7,7 @@ class NECB2020
   #
   # @return [Bool] true if successful, false if not
   #
-  # This was modified in PCF 1630
+  # NECB2020 uses a different procedure calculate gas water heater efficiencies (compared to previous NECB)
   #
   def water_heater_mixed_apply_efficiency(water_heater_mixed)
     # Get the capacity of the water heater
@@ -34,7 +34,7 @@ class NECB2020
       volume_m3 = volume_m3.get
     end
     volume_gal = OpenStudio.convert(volume_m3, 'm^3', 'gal').get
-
+    volume_litre = OpenStudio.convert(volume_m3, 'm^3', 'L').get
     # Get the heater fuel type
     fuel_type = water_heater_mixed.heaterFuelType
     unless fuel_type == 'NaturalGas' || fuel_type == 'Electricity' || fuel_type == 'FuelOilNo2'
@@ -46,21 +46,24 @@ class NECB2020
     # Calculate the energy factor (EF)
     # From PNNL http://www.energycodes.gov/sites/default/files/documents/PrototypeModelEnhancements_2014_0.pdf
     # Appendix A: Service Water Heating
-	# and modified by PCF 1630 as noted below.
+	  # and modified by PCF 1630 as noted below.
+
     water_heater_eff = nil
     ua_btu_per_hr_per_f = nil
     sl_btu_per_hr = nil
+    q_load_btu_per_hr = nil
+    uef = nil
     case fuel_type
     when 'Electricity'
-      volume_l_per_s = volume_m3 * 1000
+      volume_litre_per_s = volume_m3 * 1000
       if capacity_btu_per_hr <= OpenStudio.convert(12, 'kW', 'Btu/hr').get
         # Fixed water heater efficiency per PNNL
         water_heater_eff = 1
         # Calculate the max allowable standby loss (SL)
-        sl_w = if volume_l_per_s < 270
-                 40 + 0.2 * volume_l_per_s # assume bottom inlet
+        sl_w = if volume_litre_per_s < 270
+                 40 + 0.2 * volume_litre_per_s # assume bottom inlet
                else
-                 0.472 * volume_l_per_s - 33.5
+                 0.472 * volume_litre_per_s - 33.5
                  # assume bottom inlet
                end
         sl_btu_per_hr = OpenStudio.convert(sl_w, 'W', 'Btu/hr').get
@@ -68,73 +71,137 @@ class NECB2020
         # Fixed water heater efficiency per PNNL
         water_heater_eff = 1
         # Calculate the max allowable standby loss (SL)   # use this - NECB does not give SL calculation for cap > 12 kW
-        sl_btu_per_hr = 20 + (35 * Math.sqrt(volume_gal))
+        sl_w = 0.3 + 102.2/volume_litre_per_s
+        sl_btu_per_hr = OpenStudio.convert(sl_w, 'W', 'Btu/hr').get
       end
       # Calculate the skin loss coefficient (UA)
       ua_btu_per_hr_per_f = sl_btu_per_hr / 70
     when 'NaturalGas'
-      volume_l = volume_m3 / 1000
-      if capacity_btu_per_hr <= 75_000
-        # Fixed water heater thermal efficiency per PNNL
+      # Performance requirements from NECB2020 Table 6.2.2.1 Gas-fired storage type
+      
+      # Performance requirement based on FHR and volume
+      # Water heater parameters derived using the procedure described by:
+      #   Maguire, J., & Roberts, D. (2020). DERIVING SIMULATION PARAMETERS FOR STORAGE-TYPE WATER HEATERS 
+      #   USING RATINGS DATA PRODUCED FROM THE UNIFORM ENERGY FACTOR TEST PROCEDURE. 2020 Building Performance 
+      #   Analysis Conference and SimBuild co-organized by ASHRAE and IBPSA-USA (pp. 325-331). Chicago: ASHRAE.
+      #   https://www.ashrae.org/file%20library/conferences/specialty%20conferences/2020%20building%20performance/papers/d-bsc20-c039.pdf
+      #
+      #   AND
+      #
+      #   PNNL http://www.energycodes.gov/sites/default/files/documents/PrototypeModelEnhancements_2014_0.pdf
+
+      # Assume fhr = peak demand flow
+      tank_param = auto_size_shw_capacity(model:water_heater_mixed.model, shw_scale: 'NECB_Default')
+      fhr_L_per_hr = tank_param['loop_peak_flow_rate_SI']
+      fhr_L_per_hr = fhr_L_per_hr * 3600000
+      if capacity_w <= 22000 and volume_litre >= 76 and volume_litre < 208
+        if fhr_L_per_hr < 68
+          uef = 0.3456 - 0.00053*volume_litre
+          q_load_btu_per_hr = 5561
+          volume_drawn_gal = 10
+        elsif fhr_L_per_hr >= 68 and fhr_L_per_hr < 193
+          uef = 0.5982 - 0.00050*volume_litre
+          q_load_btu_per_hr = 21131
+          volume_drawn_gal = 38
+        elsif fhr_L_per_hr >= 193 and fhr_L_per_hr < 284
+          uef = 0.6483 - 0.00045*volume_litre
+          q_load_btu_per_hr = 30584
+          volume_drawn_gal = 55
+        elsif fhr_L_per_hr >= 284 
+          uef = 0.6920 - 0.00034*volume_litre
+          q_load_btu_per_hr = 46710
+          volume_drawn_gal = 84
+        end
+
+        # Assume burner efficiency  (PNNL)
         water_heater_eff = 0.82
-		
-        # Calculate the minimum Energy Factor (EF) (This was introduced in PCF 1630)
-        if volume_l < 68
-          uef = 0.5982 - 0.0005 * volume_l
-          ef = 1.0005 * uef + 0.0019
-        elsif volume_l >= 68 and volume_l < 193
-          uef = 0.6483 - 0.00045 * volume_l
-          ef = 1.0005 * uef + 0.0019
-        elsif volume_l >= 193 and volume_l < 284
-          uef = 0.692 - 0.00034 * volume_l
-          ef = 1.0005 * uef + 0.0019
+
+        # Estimate recovery efficiency (RE) and UA (Maguire and Robers, 2020)
+        q_load_btu = volume_drawn_gal*8.30074*0.99826*(125-58) #water properties at 91.5F
+        capacity_btu_per_hr = OpenStudio.convert(capacity_w, 'W', 'Btu/hr').get
+        re = water_heater_eff + q_load_btu*(uef-water_heater_eff)/(24*capacity_btu_per_hr*uef)
+        ua_btu_per_hr_per_f = (water_heater_eff-re)*capacity_btu_per_hr/(125-67.5)  
+
+      elsif capacity_w <= 22000 and volume_litre >= 208 and volume_litre < 380
+        if fhr_L_per_hr < 68
+          uef = 0.6470 - 0.00016*volume_litre
+          q_load_btu_per_hr = 5561
+          volume_drawn_gal = 10
+        elsif fhr_L_per_hr >= 68 and fhr_L_per_hr < 193
+          uef = 0.7689 - 0.00013*volume_litre
+          q_load_btu_per_hr = 21131
+          volume_drawn_gal = 38
+        elsif fhr_L_per_hr >= 193 and fhr_L_per_hr < 284
+          uef = 0.7897 - 0.00011*volume_litre
+          q_load_btu_per_hr = 30584
+          volume_drawn_gal = 55
+        elsif fhr_L_per_hr >= 284 
+          uef = 0.8072 - 0.00008*volume_litre
+          q_load_btu_per_hr = 46710
+          volume_drawn_gal = 84
         end
-		
-        # Calculate the Recovery Efficiency (RE)
-        # based on a fixed capacity of 75,000 Btu/hr
-        # and a fixed volume of 40 gallons by solving
-        # this system of equations:
-        # ua = (1/.95-1/re)/(67.5*(24/41094-1/(re*cap)))
-        # 0.82 = (ua*67.5+cap*re)/cap
-        cap = 75_000.0
-        re = (Math.sqrt(6724 * ef**2 * cap**2 + 40_409_100 * ef**2 * cap - 28_080_900 * ef * cap + 29_318_000_625 * ef**2 - 58_636_001_250 * ef + 29_318_000_625) + 82 * ef * cap + 171_225 * ef - 171_225) / (200 * ef * cap)
+
+        # Assume burner  efficiency  (PNNL)
+        water_heater_eff = 0.82
+
+        # Estimate recovery efficiency (RE) and UA (Maguire and Robers, 2020)
+        q_load_btu = volume_drawn_gal*8.30074*0.99826*(125-58) #water properties at 91.5F
+        capacity_btu_per_hr = OpenStudio.convert(capacity_w, 'W', 'Btu/hr').get
+        # Estimate recovery efficiency (RE) and UA (Maguire and Robers, 2020)
+        re = water_heater_eff + q_load_btu*(uef-water_heater_eff)/(24*capacity_btu_per_hr*uef)
+        ua_btu_per_hr_per_f = (water_heater_eff-re)*capacity_btu_per_hr/(125-67.5)      
+
+      elsif capacity_w > 22000 and capacity_w <= 30500 and volume_litre <= 454 
+        # NOTE: volume_litre 454L in this case, refers to manufacturer stated volume.
+        # Assume manufacturer rated volume = actual tank volume (value used in EnergyPlus)
         
-		# Calculate the skin loss coefficient (UA)
-        # based on the actual capacity.
-        ua_btu_per_hr_per_f = (water_heater_eff - re) * capacity_btu_per_hr / 67.5
-		
-      # This capacity band was introduced in PCF 1630
-      elsif capacity_btu_per_hr > 75_000 and capacity_btu_per_hr < 103977 and volume_l < 454
+        uef = 0.8107 - 0.00021*volume_litre
+
+        # Assume burner efficiency  (PNNL)
         water_heater_eff = 0.82
-        uef = 0.8107 - 0.00021 * volume_l
-        ef = 1.0005 * uef + 0.0019
-        cap = 103977
-        re = (Math.sqrt(6724 * ef**2 * cap**2 + 40_409_100 * ef**2 * cap - 28_080_900 * ef * cap + 29_318_000_625 * ef**2 - 58_636_001_250 * ef + 29_318_000_625) + 82 * ef * cap + 171_225 * ef - 171_225) / (200 * ef * cap)
+
+        # Estimate recovery efficiency (RE) and UA (Maguire and Robers, 2020)
+        capacity_btu_per_hr = OpenStudio.convert(capacity_w, 'W', 'Btu/hr').get
+        if fhr_L_per_hr < 68
+          q_load_btu_per_hr = 5561
+          volume_drawn_gal = 10
+        elsif fhr_L_per_hr >= 68 and fhr_L_per_hr < 193
+          q_load_btu_per_hr = 21131
+          volume_drawn_gal = 38
+        elsif fhr_L_per_hr >= 193 and fhr_L_per_hr < 284
+          q_load_btu_per_hr = 30584
+          volume_drawn_gal = 55
+        elsif fhr_L_per_hr >= 284 
+          q_load_btu_per_hr = 46710
+          volume_drawn_gal = 84
+        end
+        q_load_btu = volume_drawn_gal*8.30074*0.99826*(125-58) #water properties at 91.5F
+
+        # Estimate recovery efficiency (RE) and UA (Maguire and Robers, 2020)
+        re = water_heater_eff + q_load_btu*(uef-water_heater_eff)/(24*capacity_btu_per_hr*uef)
+        ua_btu_per_hr_per_f = (water_heater_eff-re)*capacity_btu_per_hr/(125-67.5)      
         
-		# Calculate the skin loss coefficient (UA)
-        # based on the actual capacity.
-        ua_btu_per_hr_per_f = (water_heater_eff - re) * capacity_btu_per_hr / 67.5
-      else
-        # Thermal efficiency (PCF 1630 update)
+      else # all other water heaters
+        capacity_kw = capacity_w/1000
+        capacity_btu_per_hr = OpenStudio.convert(capacity_w, 'W', 'Btu/hr').get
+        # thermal efficiency (NECB2020)
         et = 0.9
-        sl_w = 0.84 * capacity_btu_per_hr / 3412.412 / 0.234 + 16.57 * (volume_l ** 0.5)
-        sl_btu_per_hr = sl_w * 3.412
-        # Calculate the skin loss coefficient (UA)
-        ua_btu_per_hr_per_f = (sl_btu_per_hr * et) / 70
-        # Calculate water heater efficiency
-        water_heater_eff = (ua_btu_per_hr_per_f * 70 + capacity_btu_per_hr * et) / capacity_btu_per_hr
+        # maximum standby losses
+        sl_w = 0.84*(1.25*capacity_kw + 16.57*(volume_litre**0.5))
+        sl_btu_per_hr = OpenStudio.convert(sl_w, 'W', 'Btu/hr').get
+        ua_btu_per_hr_per_f = sl_btu_per_hr*et / 70
+        water_heater_eff = (ua_btu_per_hr_per_f*70 + capacity_btu_per_hr*et)/capacity_btu_per_hr
       end
     end
 
     # Convert to SI
-    ua_btu_per_hr_per_c = OpenStudio.convert(ua_btu_per_hr_per_f, 'Btu/hr*R', 'W/K').get
-
+    ua_w_per_k = OpenStudio.convert(ua_btu_per_hr_per_f, 'Btu/hr*R', 'W/K').get
     # Set the water heater properties
     # Efficiency
     water_heater_mixed.setHeaterThermalEfficiency(water_heater_eff)
     # Skin loss
-    water_heater_mixed.setOffCycleLossCoefficienttoAmbientTemperature(ua_btu_per_hr_per_c)
-    water_heater_mixed.setOnCycleLossCoefficienttoAmbientTemperature(ua_btu_per_hr_per_c)
+    water_heater_mixed.setOffCycleLossCoefficienttoAmbientTemperature(ua_w_per_k)
+    water_heater_mixed.setOnCycleLossCoefficienttoAmbientTemperature(ua_w_per_k)
     # TODO: Parasitic loss (pilot light)
     # PNNL document says pilot lights were removed, but IDFs
     # still have the on/off cycle parasitic fuel consumptions filled in

data/lib/openstudio-standards/standards/necb/common/btap_data.rb

@@ -51,6 +51,7 @@ class BTAPData
     @btap_data.merge!(energy_peak_data)
     @btap_data.merge!(utility(model))
     @btap_data.merge!(unmet_hours(model))
+    @btap_data.merge!outdoor_air_data(model)
 
     # Data in tables...
     @btap_data.merge!('measures_data_table' => measures_data_table(runner)) unless runner.nil?
@@ -294,7 +295,7 @@ class BTAPData
     outdoor_roofs = BTAP::Geometry::Surfaces.filter_by_surface_types(outdoor_surfaces, 'RoofCeiling')
     outdoor_floors = BTAP::Geometry::Surfaces.filter_by_surface_types(outdoor_surfaces, 'Floor')
     outdoor_subsurfaces = BTAP::Geometry::Surfaces.get_subsurfaces_from_surfaces(outdoor_surfaces)
-    ground_surfaces = BTAP::Geometry::Surfaces.filter_by_boundary_condition(surfaces, 'Ground')
+    ground_surfaces = BTAP::Geometry::Surfaces.filter_by_boundary_condition(surfaces, ['Ground', 'Foundation'])
     ground_walls = BTAP::Geometry::Surfaces.filter_by_surface_types(ground_surfaces, 'Wall')
     ground_roofs = BTAP::Geometry::Surfaces.filter_by_surface_types(ground_surfaces, 'RoofCeiling')
     ground_floors = BTAP::Geometry::Surfaces.filter_by_surface_types(ground_surfaces, 'Floor')
@@ -402,7 +403,7 @@ class BTAPData
   def envelope_exterior_surface_table
     surfaces = @model.getSurfaces.sort
     outdoor_surfaces = BTAP::Geometry::Surfaces.filter_by_boundary_condition(surfaces, 'Outdoors')
-    ground_surfaces = BTAP::Geometry::Surfaces.filter_by_boundary_condition(surfaces, 'Ground')
+    ground_surfaces = BTAP::Geometry::Surfaces.filter_by_boundary_condition(surfaces, ['Ground', 'Foundation'])
     exterior_opaque_surfaces = outdoor_surfaces + ground_surfaces
     # outdoor_surfaces.each { |surface| puts surface.name}
     # get surface table from sql
@@ -1408,6 +1409,7 @@ class BTAPData
      'energy_eui_natural_gas_gj_per_m_sq',
      'energy_eui_pumps_gj_per_m_sq',
      'energy_eui_total_gj_per_m_sq',
+     'energy_eui_heat recovery_gj_per_m_sq',
      'energy_eui_water systems_gj_per_m_sq'].each { |end_use| data[end_use] = 0.0 }
     # Get E+ End use table from sql
     table = get_sql_table_to_json(@model, 'AnnualBuildingUtilityPerformanceSummary', 'Entire Facility', 'End Uses')['table']
@@ -1495,6 +1497,126 @@ class BTAPData
     return service_water_heating
   end
 
+  # The below method (outdoor_air_data extract a couple of outputs related to outdoor air from the .html output file)
+  def outdoor_air_data(model)
+    # Store outdoor air data
+    outdoor_air_data = {}
+    #===============================================================================================================
+    airloops_total_outdoor_air_mechanical_ventilation_m3 = 0.0
+    airloops_total_outdoor_air_natural_ventilation_m3 = 0.0
+    zones_total_outdoor_air_mechanical_ventilation_m3 = 0.0
+    zones_total_outdoor_air_natural_ventilation_m3 = 0.0
+    zones_total_outdoor_air_infiltration_m3 = 0.0
+    #===============================================================================================================
+    # Total outdoor air by airLoop
+    model.getAirLoopHVACs.sort.each do |air_loop|
+      air_loop_name = air_loop.name.get.upcase
+
+      # Mechanical ventilation of all airloops
+      command = "SELECT Value
+               FROM TabularDataWithStrings
+               WHERE ReportName='OutdoorAirDetails'
+               AND ReportForString='Entire Facility'
+               AND TableName='Total Outdoor Air by AirLoop'
+               AND RowName='#{air_loop_name}'
+               AND ColumnName='Mechanical Ventilation'
+               AND Units='m3'"
+      airloops_total_outdoor_air_mechanical_ventilation_m3 += @sqlite_file.get.execAndReturnFirstDouble(command).to_f
+
+      # Natural ventilation of all airloops
+      command = "SELECT Value
+               FROM TabularDataWithStrings
+               WHERE ReportName='OutdoorAirDetails'
+               AND ReportForString='Entire Facility'
+               AND TableName='Total Outdoor Air by AirLoop'
+               AND RowName='#{air_loop_name}'
+               AND ColumnName='Natural Ventilation'
+               AND Units='m3'"
+      airloops_total_outdoor_air_natural_ventilation_m3 += @sqlite_file.get.execAndReturnFirstDouble(command).to_f
+
+    end
+
+    # Not-normalized mechanical/natural.
+    outdoor_air_data['airloops_total_outdoor_air_mechanical_ventilation_m3'] = airloops_total_outdoor_air_mechanical_ventilation_m3
+    outdoor_air_data['airloops_total_outdoor_air_natural_ventilation_m3'] = airloops_total_outdoor_air_natural_ventilation_m3
+
+    # Normalized mechanical/natural: ACH (air changes per hour)
+    outdoor_air_data['airloops_total_outdoor_air_mechanical_ventilation_ach_1_per_hr'] = airloops_total_outdoor_air_mechanical_ventilation_m3 / (@btap_data['bldg_volume_m_cu'] * 365 * 24)
+    outdoor_air_data['airloops_total_outdoor_air_natural_ventilation_ach_1_per_hr'] = airloops_total_outdoor_air_natural_ventilation_m3 / (@btap_data['bldg_volume_m_cu'] * 365 * 24)
+
+    # Normalized mechanical/natural: normalized by conditioned floor area
+    outdoor_air_data['airloops_total_outdoor_air_mechanical_ventilation_flow_per_conditioned_floor_area_m3_per_s_m2'] = airloops_total_outdoor_air_mechanical_ventilation_m3 / (@conditioned_floor_area_m_sq * 365 * 24 * 3600)
+    outdoor_air_data['airloops_total_outdoor_air_natural_ventilation_flow_per_conditioned_floor_area_m3_per_s_m2'] = airloops_total_outdoor_air_natural_ventilation_m3 / (@conditioned_floor_area_m_sq * 365 * 24 * 3600)
+
+    # Normalized mechanical/natural: normalized by exterior area
+    outdoor_air_data['airloops_total_outdoor_air_mechanical_ventilation_flow_per_exterior_area_m3_per_s_m2'] = airloops_total_outdoor_air_mechanical_ventilation_m3 / (@btap_data['bldg_exterior_area_m_sq'] * 365 * 24 * 3600)
+    outdoor_air_data['airloops_total_outdoor_air_natural_ventilation_flow_per_exterior_area_m3_per_s_m2'] = airloops_total_outdoor_air_natural_ventilation_m3 / (@btap_data['bldg_exterior_area_m_sq'] * 365 * 24 * 3600)
+
+    #===============================================================================================================
+    # Total outdoor air by zone
+    total_outdoor_air_mechanical_ventilation_zones_m3 = 0.0
+    model.getThermalZones.sort.each do |zone|
+      zone_name = zone.name.get.upcase
+
+      # Mechanical ventilation of all zones
+      command = "SELECT Value
+               FROM TabularDataWithStrings
+               WHERE ReportName='OutdoorAirDetails'
+               AND ReportForString='Entire Facility'
+               AND TableName='Total Outdoor Air by Zone'
+               AND RowName='#{zone_name}'
+               AND ColumnName='Mechanical Ventilation'
+               AND Units='m3'"
+      zones_total_outdoor_air_mechanical_ventilation_m3 += @sqlite_file.get.execAndReturnFirstDouble(command).to_f
+
+      # Natural ventilation of all zones
+      command = "SELECT Value
+               FROM TabularDataWithStrings
+               WHERE ReportName='OutdoorAirDetails'
+               AND ReportForString='Entire Facility'
+               AND TableName='Total Outdoor Air by Zone'
+               AND RowName='#{zone_name}'
+               AND ColumnName='Natural Ventilation'
+               AND Units='m3'"
+      zones_total_outdoor_air_natural_ventilation_m3 += @sqlite_file.get.execAndReturnFirstDouble(command).to_f
+
+      # Infiltration of all zones
+      command = "SELECT Value
+               FROM TabularDataWithStrings
+               WHERE ReportName='OutdoorAirDetails'
+               AND ReportForString='Entire Facility'
+               AND TableName='Total Outdoor Air by Zone'
+               AND RowName='#{zone_name}'
+               AND ColumnName='Infiltration'
+               AND Units='m3'"
+      zones_total_outdoor_air_infiltration_m3 += @sqlite_file.get.execAndReturnFirstDouble(command).to_f
+
+    end
+
+    # Not-normalized mechanical/natural/infiltration.
+    outdoor_air_data['zones_total_outdoor_air_mechanical_ventilation_m3'] = zones_total_outdoor_air_mechanical_ventilation_m3
+    outdoor_air_data['zones_total_outdoor_air_natural_ventilation_m3'] = zones_total_outdoor_air_natural_ventilation_m3
+    outdoor_air_data['zones_total_outdoor_air_infiltration_m3'] = zones_total_outdoor_air_infiltration_m3
+
+    # Normalized mechanical/natural/infiltration: ACH (air changes per hour)
+    outdoor_air_data['zones_total_outdoor_air_mechanical_ventilation_ach_1_per_hr'] = zones_total_outdoor_air_mechanical_ventilation_m3 / (@btap_data['bldg_volume_m_cu'] * 365 * 24)
+    outdoor_air_data['zones_total_outdoor_air_natural_ventilation_ach_1_per_hr'] = zones_total_outdoor_air_natural_ventilation_m3 / (@btap_data['bldg_volume_m_cu'] * 365 * 24)
+    outdoor_air_data['zones_total_outdoor_air_infiltration_ach_1_per_hr'] = zones_total_outdoor_air_infiltration_m3 / (@btap_data['bldg_volume_m_cu'] * 365 * 24)
+
+    # Normalized mechanical/natural/infiltration: normalized by conditioned floor area
+    outdoor_air_data['zones_total_outdoor_air_mechanical_ventilation_flow_per_conditioned_floor_area_m3_per_s_m2'] = zones_total_outdoor_air_mechanical_ventilation_m3 / (@conditioned_floor_area_m_sq * 365 * 24 * 3600)
+    outdoor_air_data['zones_total_outdoor_air_natural_ventilation_flow_per_conditioned_floor_area_m3_per_s_m2'] = zones_total_outdoor_air_natural_ventilation_m3 / (@conditioned_floor_area_m_sq * 365 * 24 * 3600)
+    outdoor_air_data['zones_total_outdoor_air_infiltration_flow_per_conditioned_floor_area_m3_per_s_m2'] = zones_total_outdoor_air_infiltration_m3 / (@conditioned_floor_area_m_sq * 365 * 24 * 3600)
+
+    # Normalized mechanical/natural/infiltration: normalized by exterior area
+    outdoor_air_data['zones_total_outdoor_air_mechanical_ventilation_flow_per_exterior_area_m3_per_s_m2'] = zones_total_outdoor_air_mechanical_ventilation_m3 / (@btap_data['bldg_exterior_area_m_sq'] * 365 * 24 * 3600)
+    outdoor_air_data['zones_total_outdoor_air_natural_ventilation_flow_per_exterior_area_m3_per_s_m2'] = zones_total_outdoor_air_natural_ventilation_m3 / (@btap_data['bldg_exterior_area_m_sq'] * 365 * 24 * 3600)
+    outdoor_air_data['zones_total_outdoor_air_infiltration_flow_per_exterior_area_m3_per_s_m2'] = zones_total_outdoor_air_infiltration_m3 / (@btap_data['bldg_exterior_area_m_sq'] * 365 * 24 * 3600)
+    #===============================================================================================================
+
+    return outdoor_air_data
+  end
+
   def sql_data_tables(model)
     puts 'Getting SQL Data into json...'
     start = Time.now

data/lib/openstudio-standards/standards/necb/common/btap_datapoint.rb

@@ -117,6 +117,8 @@ class BTAPDatapoint
                                        epw_file: @options[:epw_file],
                                        sizing_run_dir: File.join(@dp_temp_folder, 'sizing_folder'),
                                        primary_heating_fuel: @options[:primary_heating_fuel],
+                                       necb_reference_hp: @options[:necb_reference_hp],
+                                       necb_reference_hp_supp_fuel: @options[:necb_reference_hp_supp_fuel],
                                        dcv_type: @options[:dcv_type], # Four options: @options[: (1) 'NECB_Default', (2) 'No DCV', (3) 'Occupancy-based DCV' , (4) 'CO2-based DCV'
                                        lights_type: @options[:lights_type], # Two options: @options[: (1) 'NECB_Default', (2) 'LED'
                                        lights_scale: @options[:lights_scale],
@@ -167,7 +169,9 @@ class BTAPDatapoint
                                        output_meters: @options[:output_meters],
                                        airloop_economizer_type: @options[:airloop_economizer_type],
                                        shw_scale: @options[:shw_scale],
-                                       baseline_system_zones_map_option: @options[:baseline_system_zones_map_option])
+                                       baseline_system_zones_map_option: @options[:baseline_system_zones_map_option],
+                                       tbd_option: @options[:tbd_option]
+                                       )
       end
 
       # Save model to to disk.
@@ -388,7 +392,16 @@ class BTAPDatapoint
         FROM ReportDataDictionary
         WHERE ReportDataDictionaryIndex == #{rdd_index}
       "
-      key_value = model.sqlFile.get.execAndReturnFirstString(query).get
+
+      # In some cases KeyValue has a value and sometimes it does not.  In some cases KeyValue is null.  If the command
+      # below is run and KeyValue is null then the command fails and returns an error.  The fix below assumes that if
+      # the command below fails it is because KeyValue is null.  In that case the "key_value" variable is set to a
+      # blank.
+      begin
+        key_value = model.sqlFile.get.execAndReturnFirstString(query).get
+      rescue StandardError => bang
+        key_value = ""
+      end
 
       #Get Units
       query = "

data/lib/openstudio-standards/standards/necb/common/necb_reference_runs.csv

@@ -1,4 +1,4 @@
-name,_id,run_start_time,run_end_time,status,status_message,btap_standard_building_type_geometry.building_type,btap_results.bldg_standards_template,btap_standard_systems.primary_heating_fuel,btap_standard_weather.epw_file,btap_results.bldg_conditioned_floor_area_m_sq,btap_results.bldg_exterior_area_m_sq,btap_results.bldg_fdwr,btap_results.bldg_name,btap_results.bldg_srr,btap_results.bldg_standards_building_template,btap_results.bldg_standards_building_type,btap_results.bldg_standards_number_of_above_ground_stories,btap_results.bldg_standards_number_of_stories,btap_results.bldg_surface_to_volume_ratio,btap_results.bldg_volume_m_cu,btap_results.cost_equipment_envelope_total_cost_per_m_sq,btap_results.cost_equipment_heating_and_cooling_total_cost_per_m_sq,btap_results.cost_equipment_lighting_total_cost_per_m_sq,btap_results.cost_equipment_shw_total_cost_per_m_sq,btap_results.cost_equipment_total_cost_per_m_sq,btap_results.cost_equipment_ventilation_total_cost_per_m_sq,btap_results.cost_rs_means_city,btap_results.cost_rs_means_prov,btap_results.cost_utility_neb_electricity_cost_per_m_sq,btap_results.cost_utility_neb_natural_gas_cost_per_m_sq,btap_results.cost_utility_neb_oil_cost_per_m_sq,btap_results.cost_utility_neb_total_cost_per_m_sq,btap_results.energy_eui_additional_fuel_gj_per_m_sq,btap_results.energy_eui_cooling_gj_per_m_sq,btap_results.energy_eui_district_cooling_gj_per_m_sq,btap_results.energy_eui_district_heating_gj_per_m_sq,btap_results.energy_eui_electricity_gj_per_m_sq,btap_results.energy_eui_fans_gj_per_m_sq,btap_results.energy_eui_heating_gj_per_m_sq,btap_results.energy_eui_interior_equipment_gj_per_m_sq,btap_results.energy_eui_interior_lighting_gj_per_m_sq,btap_results.energy_eui_natural_gas_gj_per_m_sq,btap_results.energy_eui_pumps_gj_per_m_sq,btap_results.energy_eui_total_gj_per_m_sq,btap_results.energy_eui_water_systems_gj_per_m_sq,btap_results.energy_peak_electric_w_per_m_sq,btap_results.energy_peak_natural_gas_w_per_m_sq,btap_results.energy_principal_heating_source,btap_results.envelope_ground_floors_average_conductance_w_per_m_sq_k,btap_results.envelope_ground_roofs_average_conductance_w_per_m_sq_k,btap_results.envelope_ground_walls_average_conductance_w_per_m_sq_k,btap_results.envelope_outdoor_doors_average_conductance_w_per_m_sq_k,btap_results.envelope_outdoor_floors_average_conductance_w_per_m_sq_k,btap_results.envelope_outdoor_overhead_doors_average_conductance_w_per_m_sq_k,btap_results.envelope_outdoor_roofs_average_conductance_w_per_m_sq_k,btap_results.envelope_outdoor_walls_average_conductance_w_per_m_sq_k,btap_results.envelope_outdoor_windows_average_conductance_w_per_m_sq_k,btap_results.envelope_skylights_average_conductance_w_per_m_sq_k,btap_results.location_cdd,btap_results.location_city,btap_results.location_country,btap_results.location_hdd,btap_results.location_latitude,btap_results.location_longitude,btap_results.location_state_province_region,btap_results.location_weather_file,btap_results.location_zone,btap_results.measures_data_table,btap_results.shw_additional_fuel_per_year,btap_results.shw_electricity_per_day,btap_results.shw_electricity_per_day_per_occupant,btap_results.shw_electricity_per_year,btap_results.shw_natural_gas_per_year,btap_results.shw_total_nominal_occupancy,btap_results.shw_water_m_cu_per_day,btap_results.shw_water_m_cu_per_day_per_occupant,btap_results.shw_water_m_cu_per_year,btap_results.simulation_btap_data_version,btap_results.simulation_date,btap_results.simulation_os_standards_revision,btap_results.simulation_os_standards_version,btap_results.npv_total_per_m_sq
+name,_id,run_start_time,run_end_time,status,status_message,btap_standard_building_type_geometry.building_type,btap_results.bldg_standards_template,btap_standard_systems.primary_heating_fuel,btap_standard_weather.epw_file,btap_results.bldg_conditioned_floor_area_m_sq,btap_results.bldg_exterior_area_m_sq,btap_results.bldg_fdwr,btap_results.bldg_name,btap_results.bldg_srr,btap_results.bldg_standards_building_template,btap_results.bldg_standards_building_type,btap_results.bldg_standards_number_of_above_ground_stories,btap_results.bldg_standards_number_of_stories,btap_results.bldg_surface_to_volume_ratio,btap_results.bldg_volume_m_cu,btap_results.cost_equipment_envelope_total_cost_per_m_sq,btap_results.cost_equipment_heating_and_cooling_total_cost_per_m_sq,btap_results.cost_equipment_lighting_total_cost_per_m_sq,btap_results.cost_equipment_shw_total_cost_per_m_sq,btap_results.cost_equipment_total_cost_per_m_sq,btap_results.cost_equipment_ventilation_total_cost_per_m_sq,btap_results.cost_rs_means_city,btap_results.cost_rs_means_prov,btap_results.cost_utility_neb_electricity_cost_per_m_sq,btap_results.cost_utility_neb_natural_gas_cost_per_m_sq,btap_results.cost_utility_neb_oil_cost_per_m_sq,btap_results.cost_utility_neb_total_cost_per_m_sq,btap_results.energy_eui_additional_fuel_gj_per_m_sq,btap_results.energy_eui_cooling_gj_per_m_sq,btap_results.energy_eui_district_cooling_gj_per_m_sq,btap_results.energy_eui_district_heating_gj_per_m_sq,btap_results.energy_eui_electricity_gj_per_m_sq,btap_results.energy_eui_fans_gj_per_m_sq,btap_results.energy_eui_heating_gj_per_m_sq,btap_results.energy_eui_interior_equipment_gj_per_m_sq,btap_results.energy_eui_interior_lighting_gj_per_m_sq,btap_results.energy_eui_natural_gas_gj_per_m_sq,btap_results.energy_eui_pumps_gj_per_m_sq,btap_results.energy_eui_total_gj_per_m_sq,btap_results.energy_eui_water_systems_gj_per_m_sq,btap_results.energy_peak_electric_w_per_m_sq,btap_results.energy_peak_natural_gas_w_per_m_sq,btap_results.energy_principal_heating_source,btap_results.envelope_ground_floors_average_conductance_w_per_m_sq_k,btap_results.envelope_ground_roofs_average_conductance_w_per_m_sq_k,btap_results.envelope_ground_walls_average_conductance_w_per_m_sq_k,btap_results.envelope_outdoor_doors_average_conductance_w_per_m_sq_k,btap_results.envelope_outdoor_floors_average_conductance_w_per_m_sq_k,btap_results.envelope_outdoor_overhead_doors_average_conductance_w_per_m_sq_k,btap_results.envelope_outdoor_roofs_average_conductance_w_per_m_sq_k,btap_results.envelope_outdoor_walls_average_conductance_w_per_m_sq_k,btap_results.envelope_outdoor_windows_average_conductance_w_per_m_sq_k,btap_results.envelope_skylights_average_conductance_w_per_m_sq_k,btap_results.location_cdd,btap_results.location_city,btap_results.location_country,btap_results.location_hdd,btap_results.location_latitude,btap_results.location_longitude,btap_results.location_state_province_region,btap_results.location_weather_file,btap_results.location_zone,btap_results.measures_data_table,btap_results.shw_additional_fuel_per_year,btap_results.shw_electricity_per_day,btap_results.shw_electricity_per_day_per_occupant,btap_results.shw_electricity_per_year,btap_results.shw_natural_gas_per_year,btap_results.shw_total_nominal_occupancy,btap_results.shw_water_m_cu_per_day,btap_results.shw_water_m_cu_per_day_per_occupant,btap_results.shw_water_m_cu_per_year,btap_results.simulation_btap_data_version,btap_results.simulation_date,btap_results.simulation_os_standards_revision,btap_results.simulation_os_standards_version,btap_results.npv_total_per_m_sq,airloops_total_outdoor_air_mechanical_ventilation_m3,airloops_total_outdoor_air_natural_ventilation_ach_1_per_hr,airloops_total_outdoor_air_mechanical_ventilation_flow_per_conditioned_floor_area_m3_per_s_m2
 DOE Autogenerated 139,02044732-a9d8-4e84-9a37-ee894242bd76,2019-11-22T21:02:22.820Z,2019-11-22T21:04:34.271Z,completed,completed normal,Warehouse,NECB2011,Electricity,CAN_QC_Sept-Iles.AP.718110_CWEC2016.epw,4835.13,7095.242119931756,3.4,Warehouse,1.5,,Warehouse,1,1,0.1873385454726663,37873.90417722118,146.96869991086072,9.340462407422343,52.72710351117757,3.392070120141547,266.74079290525793,54.31245902385251,DALHOUSIE,NEW BRUNSWICK,7.27358941745103,0.0,0.0,7.27358941745103,0.0,0.0021736747512476396,0.0,0.0,0.35343000084796067,0.03717790421353717,0.20311553153689765,0.017571399321217835,0.08201227267932817,0.0,,0.35343000084796067,0.011379218345732174,86.32869643629023,NA,Electricity,0.757,"","",2.2,"","",0.162,0.21,2.2,2.2,14.0,Sept-Iles AP,CAN,6200.0,50.223,-66.266,QC,CAN_QC_Sept-Iles.AP.718110_CWEC2016.epw,"","[
   {
     ""measure_name"": ""btap_standard_building_type_geometry"",

data/lib/openstudio-standards/utilities/simulation.rb

@@ -103,7 +103,7 @@ Standard.class_eval do
 
       # 'touch' the weather file - for some odd reason this fixes the simulation not running issue we had on openstudio-server.
       # Removed for until further investigation completed.
-      # FileUtils.touch("#{run_dir}/#{epw_name}")
+      #FileUtils.touch("#{run_dir}/#{epw_name}")
 
       cli_path = OpenStudio.getOpenStudioCLI
       cmd = "\"#{cli_path}\" run -w \"#{osw_path}\""

data/lib/openstudio-standards/version.rb

@@ -13,5 +13,5 @@ module OpenstudioStandards
     end
     return 'git-not-found-on-this-system'
   end
-  VERSION = '0.3.0'.freeze
+  VERSION = '0.3.1.rc2'.freeze
 end

data/lib/openstudio-standards/weather/Weather.stat_file.rb

@@ -46,6 +46,8 @@ module EnergyPlus
     attr_accessor :typical_winter_dry_week
     attr_accessor :typical_autumn_week
     attr_accessor :typical_spring_week
+    attr_accessor :monthly_undis_ground_temps_0p5m
+    attr_accessor :monthly_undis_ground_temps_4p0m
 
     def initialize(path)
       @path = Pathname.new(path)
@@ -73,6 +75,8 @@ module EnergyPlus
       @typical_winter_dry_week = []
       @typical_autumn_week = []
       @typical_spring_week = []
+      @monthly_undis_ground_temps_0p5m = []
+      @monthly_undis_ground_temps_4p0m = []
       @data = []
       init
     end
@@ -313,7 +317,7 @@ module EnergyPlus
         # use regex to get the temperatures
         regex = /Daily Avg(.*)\n/
         match_data = text.match(regex)
-        if match_data.nil?
+       if match_data.nil?
           puts "Can't find outdoor air temps"
           raise
         else
@@ -331,6 +335,34 @@ module EnergyPlus
           # puts "#{@monthly_dry_bulb}"
         end
 
+        # use regex to get undisturbed ground temps at 0.5 m depth
+        regex = /Monthly.*Calculated.*undisturbed*.*Ground.*Temperatures.*\n.*Jan.*Feb.*Mar.*Apr.*May.*Jun.*Jul.*Aug.*Sep.*Oct.*Nov.*Dec.*\n(.*)\n(.*)\n(.*)/
+        match_data = text.match(regex)
+        if match_data.nil?
+          puts "Can't find undisturbed ground temperatures"
+          raise
+        else
+          # first match is undisturbed ground temperature at 0.5 m and 4.0 m depth
+          monthly_undis_ground_temps_0p5m = match_data[1].strip.split(/\s+/)
+          monthly_undis_ground_temps_0p5m.shift
+          monthly_undis_ground_temps_0p5m.shift
+          monthly_undis_ground_temps_4p0m = match_data[3].strip.split(/\s+/)
+          monthly_undis_ground_temps_4p0m.shift
+          monthly_undis_ground_temps_4p0m.shift
+          # have to be 12 months
+          if monthly_undis_ground_temps_0p5m.size != 12
+            puts "Can't find undisturbed ground temps at 0.5 m"
+            raise
+          end
+          if monthly_undis_ground_temps_4p0m.size != 12
+            puts "Can't find undisturbed ground temps at 4.0 m"
+            raise
+          end
+          # insert as numbers
+          monthly_undis_ground_temps_0p5m.each { |temp| @monthly_undis_ground_temps_0p5m << temp.to_f }
+          monthly_undis_ground_temps_4p0m.each { |temp| @monthly_undis_ground_temps_4p0m << temp.to_f }
+        end
+
         # now we are valid
         @valid = true
 

data/lib/openstudio-standards.rb

@@ -33,6 +33,7 @@ module OpenstudioStandards
   # Standards classes
   require_relative "#{stds}/standard"
   # NECB2011 Code
+  require_relative "#{stds}/necb/NECB2011/system_fuels"
   require_relative "#{stds}/necb/NECB2011/necb_2011"
   require_relative "#{stds}/necb/NECB2011/building_envelope"
   require_relative "#{stds}/necb/NECB2011/lighting"