openstudio-standards 0.2.14 → 0.2.15.pre.rc1

data/lib/openstudio-standards/standards/ashrae_90_1/ashrae_90_1_2019/ashrae_90_1_2019.AirLoopHVAC.rb

@@ -1,28 +1,10 @@
 class ASHRAE9012019 < ASHRAE901
   # @!group AirLoopHVAC
 
-  # Apply multizone vav outdoor air method and
-  # adjust multizone VAV damper positions
-  # to achieve a system minimum ventilation effectiveness
-  # of 0.6 per PNNL.  Hard-size the resulting min OA
-  # into the sizing:system object.
-  #
-  # return [Bool] returns true if successful, false if not
-  # @todo move building-type-specific code to Prototype classes
-  def air_loop_hvac_apply_multizone_vav_outdoor_air_sizing(air_loop_hvac)
-    # First time adjustment:
-    # Only applies to multi-zone vav systems
-    # exclusion: for Outpatient: (1) both AHU1 and AHU2 in 'DOE Ref Pre-1980' and 'DOE Ref 1980-2004'
-    # (2) AHU1 in 2004-2019
-    # TODO refactor: move building-type-specific code to Prototype classes
-    if air_loop_hvac_multizone_vav_system?(air_loop_hvac) && !(air_loop_hvac.name.to_s.include? 'Outpatient F1')
-      air_loop_hvac_adjust_minimum_vav_damper_positions(air_loop_hvac)
-    end
-
-    return true
-  end
-
   # Determine the limits for the type of economizer present on the AirLoopHVAC, if any.
+  #
+  # @param air_loop_hvac [OpenStudio::Model::AirLoopHVAC] air loop
+  # @param climate_zone [String] ASHRAE climate zone, e.g. 'ASHRAE 169-2013-4A'
   # @return [Array<Double>] [drybulb_limit_f, enthalpy_limit_btu_per_lb, dewpoint_limit_f]
   def air_loop_hvac_economizer_limits(air_loop_hvac, climate_zone)
     drybulb_limit_f = nil
@@ -31,11 +13,9 @@ class ASHRAE9012019 < ASHRAE901
 
     # Get the OA system and OA controller
     oa_sys = air_loop_hvac.airLoopHVACOutdoorAirSystem
-    if oa_sys.is_initialized
-      oa_sys = oa_sys.get
-    else
-      return [nil, nil, nil] # No OA system
-    end
+    return [nil, nil, nil] unless oa_sys.is_initialized
+
+    oa_sys = oa_sys.get
     oa_control = oa_sys.getControllerOutdoorAir
     economizer_type = oa_control.getEconomizerControlType
     oa_control.resetEconomizerMinimumLimitDryBulbTemperature
@@ -67,17 +47,21 @@ class ASHRAE9012019 < ASHRAE901
 
   # Determine if the system economizer must be integrated or not.
   # All economizers must be integrated in 90.1-2019
+  #
+  # @param air_loop_hvac [OpenStudio::Model::AirLoopHVAC] air loop
+  # @param climate_zone [String] ASHRAE climate zone, e.g. 'ASHRAE 169-2013-4A'
+  # @return [Bool] returns true if required, false if not
   def air_loop_hvac_integrated_economizer_required?(air_loop_hvac, climate_zone)
-    integrated_economizer_required = true
-    return integrated_economizer_required
+    return true
   end
 
   # Check the economizer type currently specified in the ControllerOutdoorAir object on this air loop
   # is acceptable per the standard.
   #
-  # @param (see #economizer_required?)
+  # @param air_loop_hvac [OpenStudio::Model::AirLoopHVAC] air loop
+  # @param climate_zone [String] ASHRAE climate zone, e.g. 'ASHRAE 169-2013-4A'
   # @return [Bool] Returns true if allowable, if the system has no economizer or no OA system.
-  # Returns false if the economizer type is not allowable.
+  #   Returns false if the economizer type is not allowable.
   def air_loop_hvac_economizer_type_allowable?(air_loop_hvac, climate_zone)
     # EnergyPlus economizer types
     # 'NoEconomizer'
@@ -91,11 +75,9 @@ class ASHRAE9012019 < ASHRAE901
 
     # Get the OA system and OA controller
     oa_sys = air_loop_hvac.airLoopHVACOutdoorAirSystem
-    if oa_sys.is_initialized
-      oa_sys = oa_sys.get
-    else
-      return true # No OA system
-    end
+    return true unless oa_sys.is_initialized
+
+    oa_sys = oa_sys.get
     oa_control = oa_sys.getControllerOutdoorAir
     economizer_type = oa_control.getEconomizerControlType
 
@@ -162,13 +144,12 @@ class ASHRAE9012019 < ASHRAE901
   end
 
   # Determine if multizone vav optimization is required.
+  # @note code_sections [90.1-2019_6.5.3.3]
   #
-  # @code_sections [90.1-2019_6.5.3.3]
-  # @param (see #economizer_required?)
-  # @return [Bool] Returns true if required, false if not.
-  # @todo Add exception logic for
-  #   systems with AIA healthcare ventilation requirements
-  #   dual duct systems
+  # @param air_loop_hvac [OpenStudio::Model::AirLoopHVAC] air loop
+  # @param climate_zone [String] ASHRAE climate zone, e.g. 'ASHRAE 169-2013-4A'
+  # @return [Bool] returns true if required, false if not
+  # @todo Add exception logic for systems with AIA healthcare ventilation requirements dual duct systems
   def air_loop_hvac_multizone_vav_optimization_required?(air_loop_hvac, climate_zone)
     multizone_opt_required = false
 
@@ -231,7 +212,7 @@ class ASHRAE9012019 < ASHRAE901
       return multizone_opt_required
     end
 
-    # TODO: Not required for dual-duct systems
+    # @todo Not required for dual-duct systems
     # if self.isDualDuct
     # OpenStudio::logFree(OpenStudio::Info, "openstudio.standards.AirLoopHVAC", "For #{controller_oa.name}: multizone optimization is not applicable because it is a dual duct system")
     # return multizone_opt_required
@@ -244,9 +225,9 @@ class ASHRAE9012019 < ASHRAE901
   end
 
   # Determines the OA flow rates above which an economizer is required.
-  # Two separate rates, one for systems with an economizer and another
-  # for systems without.
-  # are zero for both types.
+  # Two separate rates, one for systems with an economizer and another for systems without.
+  #
+  # @param air_loop_hvac [OpenStudio::Model::AirLoopHVAC] air loop
   # @return [Array<Double>] [min_oa_without_economizer_cfm, min_oa_with_economizer_cfm]
   def air_loop_hvac_demand_control_ventilation_limits(air_loop_hvac)
     min_oa_without_economizer_cfm = 3000
@@ -255,7 +236,10 @@ class ASHRAE9012019 < ASHRAE901
   end
 
   # Determine the air flow and number of story limits for whether motorized OA damper is required.
-  # @return [Array<Double>] [minimum_oa_flow_cfm, maximum_stories]
+  #
+  # @param air_loop_hvac [OpenStudio::Model::AirLoopHVAC] air loop
+  # @param climate_zone [String] ASHRAE climate zone, e.g. 'ASHRAE 169-2013-4A'
+  # @return [Array<Double>] [minimum_oa_flow_cfm, maximum_stories]. If both nil, never required
   def air_loop_hvac_motorized_oa_damper_limits(air_loop_hvac, climate_zone)
     case climate_zone
     when 'ASHRAE 169-2006-0A',
@@ -289,6 +273,8 @@ class ASHRAE9012019 < ASHRAE901
   # Determine the number of stages that should be used as controls for single zone DX systems.
   # 90.1-2019 depends on the cooling capacity of the system.
   #
+  # @param air_loop_hvac [OpenStudio::Model::AirLoopHVAC] air loop
+  # @param climate_zone [String] ASHRAE climate zone, e.g. 'ASHRAE 169-2013-4A'
   # @return [Integer] the number of stages: 0, 1, 2
   def air_loop_hvac_single_zone_controls_num_stages(air_loop_hvac, climate_zone)
     min_clg_cap_btu_per_hr = 65_000
@@ -305,10 +291,11 @@ class ASHRAE9012019 < ASHRAE901
   end
 
   # Determine if the system required supply air temperature (SAT) reset.
-  # For 90.1-2019, SAT reset requirements are based on climate zone.
+  # For 90.1-2019, SAT reset requirements are based on climate zone. More exceptions are added for 90.1 2019 6.5.3.5
   #
-  # @param (see #economizer_required?)
-  # @return [Bool] Returns true if required, false if not.
+  # @param air_loop_hvac [OpenStudio::Model::AirLoopHVAC] air loop
+  # @param climate_zone [String] ASHRAE climate zone, e.g. 'ASHRAE 169-2013-4A'
+  # @return [Bool] returns true if required, false if not
   def air_loop_hvac_supply_air_temperature_reset_required?(air_loop_hvac, climate_zone)
     is_sat_reset_required = false
 
@@ -317,16 +304,61 @@ class ASHRAE9012019 < ASHRAE901
       return is_sat_reset_required
     end
 
+    # check if design outside air is less than 10,000cfm (5000L/s) 90.1 2019 6.5.3.5 Exception 1 and 2
+    design_oa_m3s = nil
+    if air_loop_hvac.sizingSystem.designOutdoorAirFlowRate.is_initialized
+      design_oa_m3s = air_loop_hvac.sizingSystem.designOutdoorAirFlowRate.get
+    elsif air_loop_hvac.sizingSystem.autosizedDesignOutdoorAirFlowRate.is_initialized
+      design_oa_m3s = air_loop_hvac.sizingSystem.autosizedDesignOutdoorAirFlowRate.get
+    else
+      OpenStudio.logFree(OpenStudio::Warn, 'openstudio.standards.AirLoopHVAC', "For #{air_loop_hvac.name} design outdoor air flow rate is not available.")
+    end
+    design_oa_cfm = OpenStudio.convert(design_oa_m3s, 'm^3/s', 'cfm').get
+
+    # check if there is erv 90.1 2019 Exceptions to 6.5.3.5 Exception 3
+    has_erv = air_loop_hvac_energy_recovery?(air_loop_hvac)
+    design_sa_m3s = air_loop_hvac_find_design_supply_air_flow_rate(air_loop_hvac)
+
+    oa_ratio = 0
+    if design_sa_m3s > 0
+      oa_ratio = design_oa_m3s / design_sa_m3s
+    end
+    has_large_oa = (oa_ratio >= 0.8)
+
     case climate_zone
     when 'ASHRAE 169-2006-0A',
          'ASHRAE 169-2006-1A',
-         'ASHRAE 169-2006-2A',
          'ASHRAE 169-2006-3A',
          'ASHRAE 169-2013-0A',
          'ASHRAE 169-2013-1A',
-         'ASHRAE 169-2013-2A',
          'ASHRAE 169-2013-3A'
-      OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.AirLoopHVAC', "For #{air_loop_hvac.name}: Supply air temperature reset is not required per 6.5.3.4 Exception 1, the system is located in climate zone #{climate_zone}.")
+      if design_oa_cfm < 3000
+        is_sat_reset_required = false
+        OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.AirLoopHVAC', "For #{air_loop_hvac.name}: Supply air temperature reset is not required per 6.5.3.5 Exception 1, the system is located in climate zone #{climate_zone}.")
+        return is_sat_reset_required
+      end
+      if has_erv && has_large_oa
+        is_sat_reset_required = false
+        OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.AirLoopHVAC', "For #{air_loop_hvac.name}: Supply air temperature reset is not required per 6.5.3.5 Exception 3, the system is located in climate zone #{climate_zone}.")
+        return is_sat_reset_required
+      end
+      is_sat_reset_required = true
+      OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.AirLoopHVAC', "For #{air_loop_hvac.name}: Supply air temperature reset is required.")
+      return is_sat_reset_required
+    when 'ASHRAE 169-2006-2A',
+         'ASHRAE 169-2013-2A'
+      if design_oa_cfm < 10000
+        is_sat_reset_required = false
+        OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.AirLoopHVAC', "For #{air_loop_hvac.name}: Supply air temperature reset is not required per 6.5.3.5 Exception 2, the system is located in climate zone #{climate_zone}.")
+        return is_sat_reset_required
+      end
+      if has_erv && has_large_oa
+        is_sat_reset_required = false
+        OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.AirLoopHVAC', "For #{air_loop_hvac.name}: Supply air temperature reset is not required per 6.5.3.5 Exception 3, the system is located in climate zone #{climate_zone}.")
+        return is_sat_reset_required
+      end
+      is_sat_reset_required = true
+      OpenStudio.logFree(OpenStudio::Info, 'openstudio.standards.AirLoopHVAC', "For #{air_loop_hvac.name}: Supply air temperature reset is required.")
       return is_sat_reset_required
     when 'ASHRAE 169-2006-0B',
          'ASHRAE 169-2006-1B',
@@ -370,8 +402,11 @@ class ASHRAE9012019 < ASHRAE901
 
   # Determine the airflow limits that govern whether or not an ERV is required.
   # Based on climate zone and % OA, plus the number of operating hours the system has.
-  # @return [Double] the flow rate above which an ERV is required.
-  # if nil, ERV is never required.
+  #
+  # @param air_loop_hvac [OpenStudio::Model::AirLoopHVAC] air loop
+  # @param climate_zone [String] ASHRAE climate zone, e.g. 'ASHRAE 169-2013-4A'
+  # @param pct_oa [Double] percentage of outdoor air
+  # @return [Double] the flow rate above which an ERV is required. if nil, ERV is never required.
   def air_loop_hvac_energy_recovery_ventilator_flow_limit(air_loop_hvac, climate_zone, pct_oa)
     # Calculate the number of system operating hours
     # based on the availability schedule.
@@ -419,11 +454,26 @@ class ASHRAE9012019 < ASHRAE901
         erv_cfm = energy_recovery_limits['greater_than_80_percent_oa']
       end
     else
+      # Check if air loop serves a non-transient dwelling unit,
+      # currently non-transient dwelling units are residential
+      # spaces in the apartment prototypes
+      building_data = model_get_building_climate_zone_and_building_type(air_loop_hvac.model)
+      building_type = building_data['building_type']
+      nontrans_dwel = false
+      if building_type == 'MidriseApartment' || building_type == 'HighriseApartment'
+        air_loop_hvac.thermalZones.each do |zone|
+          next unless thermal_zone_residential?(zone)
+
+          nontrans_dwel = true
+        end
+      end
+
       # Table 6.5.6.1-2, above 8000 hrs
       search_criteria = {
         'template' => template,
         'climate_zone' => climate_zone,
-        'under_8000_hours' => false
+        'under_8000_hours' => false,
+        'nontransient_dwelling' => nontrans_dwel
       }
       energy_recovery_limits = model_find_object(standards_data['energy_recovery'], search_criteria)
       if energy_recovery_limits.nil?
@@ -431,7 +481,11 @@ class ASHRAE9012019 < ASHRAE901
         return nil
       end
       if pct_oa < 0.1
-        erv_cfm = nil
+        if nontrans_dwel
+          erv_cfm = energy_recovery_limits['0_to_10_percent_oa']
+        else
+          erv_cfm = nil
+        end
       elsif pct_oa >= 0.1 && pct_oa < 0.2
         erv_cfm = energy_recovery_limits['10_to_20_percent_oa']
       elsif pct_oa >= 0.2 && pct_oa < 0.3
@@ -453,4 +507,236 @@ class ASHRAE9012019 < ASHRAE901
 
     return erv_cfm
   end
+
+  # Adjust minimum VAV damper positions and set minimum design
+  # system outdoor air flow following ASHRAE Std. 62.1-2019
+  #
+  # @param (see #economizer_required?)
+  # @return [Bool] Returns true if required, false if not.
+  # @todo Add exception logic for systems serving parking garage, warehouse, or multifamily
+  def air_loop_hvac_adjust_minimum_vav_damper_positions(air_loop_hvac)
+    # Do not apply the adjustment to some of the system in
+    # the hospital and outpatient which have their minimum
+    # damper position determined based on AIA 2001 ventilation
+    # requirements
+    if (@instvarbuilding_type == 'Hospital' && (air_loop_hvac.name.to_s.include?('VAV_ER') || air_loop_hvac.name.to_s.include?('VAV_ICU') ||
+                                                air_loop_hvac.name.to_s.include?('VAV_OR') || air_loop_hvac.name.to_s.include?('VAV_LABS') ||
+                                                air_loop_hvac.name.to_s.include?('VAV_PATRMS'))) ||
+       (@instvarbuilding_type == 'Outpatient' && air_loop_hvac.name.to_s.include?('Outpatient F1'))
+
+      return true
+    end
+
+    # Total uncorrected outdoor airflow rate
+    v_ou = 0.0
+    air_loop_hvac.thermalZones.each do |zone|
+      # Vou is the system uncorrected outdoor airflow:
+      # Zone airflow is multiplied by the zone multiplier
+      v_ou += thermal_zone_outdoor_airflow_rate(zone) * zone.multiplier.to_f
+    end
+
+    v_ou_cfm = OpenStudio.convert(v_ou, 'm^3/s', 'cfm').get
+
+    OpenStudio.logFree(OpenStudio::Debug, 'openstudio.standards.AirLoopHVAC', "For #{air_loop_hvac.name}: v_ou = #{v_ou_cfm.round} cfm.")
+
+    # Retrieve the sum of the zone minimum primary airflow
+    vpz_min_sum = air_loop_hvac.autosizeSumMinimumHeatingAirFlowRates
+
+    air_loop_hvac.thermalZones.sort.each do |zone|
+      # Breathing zone airflow rate
+      v_bz = thermal_zone_outdoor_airflow_rate(zone)
+
+      # Zone air distribution, assumed 1 per PNNL
+      e_z = 1.0
+
+      # Zone airflow rate
+      v_oz = v_bz / e_z
+
+      # Primary design airflow rate
+      # max of heating and cooling
+      # design air flow rates
+      v_pz = 0.0
+      clg_dsn_flow = zone.autosizedCoolingDesignAirFlowRate
+      if clg_dsn_flow.is_initialized
+        clg_dsn_flow = clg_dsn_flow.get
+        if clg_dsn_flow > v_pz
+          v_pz = clg_dsn_flow
+        end
+      else
+        OpenStudio.logFree(OpenStudio::Warn, 'openstudio.standards.AirLoopHVAC', "For #{air_loop_hvac.name}: #{zone.name} clg_dsn_flow could not be found.")
+      end
+      htg_dsn_flow = zone.autosizedHeatingDesignAirFlowRate
+      if htg_dsn_flow.is_initialized
+        htg_dsn_flow = htg_dsn_flow.get
+        if htg_dsn_flow > v_pz
+          v_pz = htg_dsn_flow
+        end
+      else
+        OpenStudio.logFree(OpenStudio::Warn, 'openstudio.standards.AirLoopHVAC', "For #{air_loop_hvac.name}: #{zone.name} htg_dsn_flow could not be found.")
+      end
+
+      # Zone ventilation efficiency calculation is computed
+      # on a per zone basis, the zone primary airflow is
+      # adjusted to removed the zone multiplier
+      v_pz /= zone.multiplier.to_f
+
+      # Set minimum damper position
+      air_loop_hvac_set_minimum_damper_position(zone, [0.01, [1.5 * v_oz / v_pz, 1.0].min].max.round(3))
+    end
+
+    # Occupant diversity (D): Ps / sum(Pz)
+    # Current value is based on school prototypes
+    # which are assumed to have the most diversity
+    occ_diver_d = 0.66
+
+    # From ASHRAE Std 62.1-2019 Section 6.2.5.3
+    if occ_diver_d < 0.6
+      e_v = 0.88 * occ_diver_d + 0.22
+    else
+      e_v = 0.75
+    end
+
+    # Total system outdoor intake flow rate
+    v_ot = v_ou / e_v
+    v_ot_cfm = OpenStudio.convert(v_ot, 'm^3/s', 'cfm').get
+
+    # Get maximum OA fraction schedule
+    oa_ctrl = air_loop_hvac.airLoopHVACOutdoorAirSystem.get.getControllerOutdoorAir
+    max_oa_frac_sch = oa_ctrl.maximumFractionofOutdoorAirSchedule
+
+    if !max_oa_frac_sch.is_initialized
+      max_oa_frac_sch = OpenStudio::Model::ScheduleConstant.new(air_loop_hvac.model)
+      max_oa_frac_sch.setName("#{air_loop_hvac.name}_MAX_OA_FRAC")
+      max_oa_frac_sch.setValue(1.0)
+      max_oa_frac_sch_type = 'Schedule:Constant'
+      oa_ctrl.setMaximumFractionofOutdoorAirSchedule(max_oa_frac_sch)
+    else
+      if max_oa_frac_sch.to_ScheduleRuleset.is_initialized
+        max_oa_frac_sch = max_oa_frac_sch.to_ScheduleRuleset.get
+        max_oa_frac_sch_type = 'Schedule:Year'
+      elsif max_oa_frac_sch.to_ScheduleConstant.is_initialized
+        max_oa_frac_sch = max_oa_frac_sch.to_ScheduleConstant.get
+        max_oa_frac_sch_type = 'Schedule:Constant'
+      elsif max_oa_frac_sch.to_ScheduleCompact.is_initialized
+        max_oa_frac_sch = max_oa_frac_sch.to_ScheduleCompact.get
+        max_oa_frac_sch_type = 'Schedule:Compact'
+      end
+    end
+
+    # Add EMS to "cap" the OA calculated by the
+    # Controller:MechanicalVentilation object
+    # to the design v_ot using the maximum OA
+    # fraction schedule
+
+    # Add EMS sensors
+    # OA mass flow calculated by the Controller:MechanicalVentilation
+    air_loop_hvac_name_ems = "EMS_#{air_loop_hvac.name.to_s.gsub(' ', '_')}"
+    oa_vrp_mass_flow = OpenStudio::Model::EnergyManagementSystemSensor.new(air_loop_hvac.model, 'Air System Outdoor Air Mechanical Ventilation Requested Mass Flow Rate')
+    oa_vrp_mass_flow.setKeyName(air_loop_hvac.name.to_s)
+    oa_vrp_mass_flow.setName("#{air_loop_hvac_name_ems}_OA_VRP")
+    # Actual sensed OA mass flow
+    oa_mass_flow = OpenStudio::Model::EnergyManagementSystemSensor.new(air_loop_hvac.model, 'Air System Outdoor Air Mass Flow Rate')
+    oa_mass_flow.setKeyName(air_loop_hvac.name.to_s)
+    oa_mass_flow.setName("#{air_loop_hvac_name_ems}_OA")
+    # Actual sensed volumetric OA flow
+    oa_vol_flow = OpenStudio::Model::EnergyManagementSystemSensor.new(air_loop_hvac.model, 'System Node Standard Density Volume Flow Rate')
+    oa_vol_flow.setKeyName("#{air_loop_hvac.name} Mixed Air Node")
+    oa_vol_flow.setName("#{air_loop_hvac_name_ems}_SUPPLY_FLOW")
+
+    # Add EMS actuator
+    max_oa_fraction = OpenStudio::Model::EnergyManagementSystemActuator.new(max_oa_frac_sch, max_oa_frac_sch_type, 'Schedule Value')
+    max_oa_fraction.setName("#{air_loop_hvac_name_ems}_MAX_OA_FRAC")
+
+    # Add EMS program
+    max_oa_ems_prog = OpenStudio::Model::EnergyManagementSystemProgram.new(air_loop_hvac.model)
+    max_oa_ems_prog.setName("#{air_loop_hvac.name}_MAX_OA_FRAC")
+    max_oa_ems_prog_body = <<-EMS
+    IF #{air_loop_hvac_name_ems}_OA > #{air_loop_hvac_name_ems}_OA_VRP,
+    SET #{air_loop_hvac_name_ems}_MAX_OA_FRAC = NULL,
+    ELSE,
+    IF #{air_loop_hvac_name_ems}_SUPPLY_FLOW > 0,
+    SET #{air_loop_hvac_name_ems}_MAX_OA_FRAC = #{v_ot} / #{air_loop_hvac_name_ems}_SUPPLY_FLOW,
+    ELSE,
+    SET #{air_loop_hvac_name_ems}_MAX_OA_FRAC = NULL,
+    ENDIF,
+    ENDIF
+    EMS
+    max_oa_ems_prog.setBody(max_oa_ems_prog_body)
+
+    max_oa_ems_prog_manager = OpenStudio::Model::EnergyManagementSystemProgramCallingManager.new(air_loop_hvac.model)
+    max_oa_ems_prog_manager.setName("SET_#{air_loop_hvac.name.to_s.gsub(' ', '_')}_MAX_OA_FRAC")
+    max_oa_ems_prog_manager.setCallingPoint('InsideHVACSystemIterationLoop')
+    max_oa_ems_prog_manager.addProgram(max_oa_ems_prog)
+
+    # Hard-size the sizing:system
+    # object with the calculated min OA flow rate
+    sizing_system = air_loop_hvac.sizingSystem
+    sizing_system.setDesignOutdoorAirFlowRate(v_ot)
+    sizing_system.setSystemOutdoorAirMethod('ZoneSum')
+
+    return true
+  end
+
+  # Add occupant standby controls to air loop
+  # When the thermostat schedule is setup or setback
+  # the ventilation is shutoff. Currently this is done
+  # by scheduling air terminal dampers (so load can
+  # still be met) and cycling unitary system fans
+  #
+  # @param air_loop_hvac [OpenStudio::model::AirLoopHVAC] OpenStudio AirLoopHVAC object
+  # @param standby_mode_space [Array] List of all spaces required to have standby mode controls
+  # @return [Boolean] true if sucessful, false otherwise
+  def air_loop_hvac_standby_mode_occupancy_control(air_loop_hvac, standby_mode_spaces)
+    if air_loop_hvac_include_unitary_system?(air_loop_hvac)
+      unitary_system = nil
+      # Get unitary system
+      air_loop_hvac.supplyComponents.each do |comp|
+        if comp.to_AirLoopHVACUnitarySystem.is_initialized
+          unitary_system = comp.to_AirLoopHVACUnitarySystem.get
+        end
+      end
+      return false unless !unitary_system.nil?
+
+      # Set fan operating schedule during assumed occupant standby mode time to 0 so the fan can cycle
+      new_sch = model_set_schedule_value(unitary_system.supplyAirFanOperatingModeSchedule.get, '12' => 0)
+      unitary_system.setSupplyAirFanOperatingModeSchedule(new_sch) unless new_sch.nil?
+    else
+      # Get thermal zones associated with spaces having standby mode occupancy requirements
+      standby_mode_zones = []
+      standby_mode_spaces.sort.each do |space|
+        standby_mode_zones << space.thermalZone.get
+      end
+      # Schedule the MDP of terminals to a low value during occupant standby mode
+      # The intent is to reduce ventilation while still allowing the terminal to
+      # meet loads
+      standby_mode_zones.each do |zone|
+        air_terminal = zone.airLoopHVACTerminal
+        if air_terminal.is_initialized
+          air_terminal = air_terminal.get
+          if air_terminal.to_AirTerminalSingleDuctVAVReheat.is_initialized
+            air_terminal = air_terminal.to_AirTerminalSingleDuctVAVReheat.get
+            if air_terminal.zoneMinimumAirFlowInputMethod == 'Constant' || air_terminal.zoneMinimumAirFlowInputMethod == 'FixedFlow'
+              if air_terminal.zoneMinimumAirFlowInputMethod == 'FixedFlow'
+                mdp_org = air_terminal.fixedMinimumAirFlowRate.get / air_terminal.autosizedMaximumAirFlowRate.get
+                air_terminal.setFixedMinimumAirFlowRate(0)
+              else
+                mdp_org = air_terminal.constantMinimumAirFlowFraction.get
+                air_terminal.setConstantMinimumAirFlowFraction(0)
+              end
+              air_terminal.setZoneMinimumAirFlowInputMethod('Scheduled')
+              air_terminal.setMinimumAirFlowFractionSchedule(model_set_schedule_value(model_add_constant_schedule_ruleset(air_loop_hvac.model, mdp_org, name = "#{air_terminal.name} - MDP", sch_type_limit: 'Fraction'), '12' => 0.1))
+            elsif air_terminal.zoneMinimumAirFlowInputMethod == 'Scheduled'
+              air_terminal.setMinimumAirFlowFractionSchedule(model_set_schedule_value(air_terminal.minimumAirFlowFractionSchedule.get, '12' => 0.1))
+            else
+              OpenStudio.logFree(OpenStudio::Warn, 'openstudio.ashrae_90_1_2019.AirLoopHVAC', "The air terminal associated with #{zone.name} uses a zone minimum air flow input method that is currently not supported so occupant standby controls were not modeled.")
+            end
+          else
+            OpenStudio.logFree(OpenStudio::Warn, 'openstudio.ashrae_90_1_2019.AirLoopHVAC', "The air terminal associated with #{zone.name} isn't of the SingleDuctVAVReheat type so occupant standby controls were not modeled.")
+          end
+        end
+      end
+    end
+
+    return true
+  end
 end

data/lib/openstudio-standards/standards/ashrae_90_1/ashrae_90_1_2019/ashrae_90_1_2019.AirTerminalSingleDuctVAVReheat.rb

@@ -2,10 +2,11 @@ class ASHRAE9012019 < ASHRAE901
   # @!group AirTerminalSingleDuctVAVReheat
 
   # Specifies the minimum damper position for VAV dampers.
-  # For terminals with hot water heat and DDC, the minimum is 20%,
-  # otherwise the minimum is 30%.
+  # For terminals with hot water heat and DDC, the minimum is 20%, otherwise the minimum is 30%.
   #
+  # @param air_terminal_single_duct_vav_reheat [OpenStudio::Model::AirTerminalSingleDuctVAVReheat] the air terminal object
   # @param has_ddc [Bool] whether or not there is DDC control of the VAV terminal in question
+  # @return [Double] minimum damper position
   def air_terminal_single_duct_vav_reheat_minimum_damper_position(air_terminal_single_duct_vav_reheat, has_ddc = false)
     min_damper_position = nil
     case air_terminal_single_duct_vav_reheat_reheat_type(air_terminal_single_duct_vav_reheat)

data/lib/openstudio-standards/standards/ashrae_90_1/ashrae_90_1_2019/ashrae_90_1_2019.CoolingTower.rb

@@ -3,17 +3,14 @@
 module ASHRAE9012019CoolingTower
   # @!group CoolingTower
 
-  # Above this point, centrifugal fan cooling towers must meet the limits
-  # of propeller or axial cooling towers instead.
-  # 90.1 6.5.5.3 Limit on Centrifugal Fan Open Circuit Cooling Towers.
-  # is 1,100 gallons per minute.
+  # Above this point, centrifugal fan cooling towers must meet the limits of propeller or axial cooling towers instead.
+  # @note code_sections 6.5.5.3 limit on centrifugal fan open=circuit cooling towers is 1,100 gallons per minute.
   #
-  # @param cooling_tower [OpenStudio::Model::CoolingTowerSingleSpeed,
-  # OpenStudio::Model::CoolingTowerTwoSpeed,
-  # OpenStudio::Model::CoolingTowerVariableSpeed] the cooling tower
+  # @param cooling_tower [OpenStudio::Model::StraightComponent] cooling tower object, allowable types:
+  #   CoolingTowerSingleSpeed, CoolingTowerTwoSpeed, CoolingTowerVariableSpeed
   # @return [Double] the limit, in gallons per minute.  Return nil for no limit.
   def cooling_tower_apply_minimum_power_per_flow_gpm_limit(cooling_tower)
-    gpm_limit = 1100
+    gpm_limit = 1100.0
     return gpm_limit
   end
 end

data/lib/openstudio-standards/standards/ashrae_90_1/ashrae_90_1_2019/ashrae_90_1_2019.CoolingTowerVariableSpeed.rb

@@ -5,12 +5,11 @@ class ASHRAE9012019 < ASHRAE901
 
   # Apply the efficiency, plus Multicell heat rejection with VSD
   #
-  # @param cooling_tower_variable_speed [OpenStudio::Model::CoolingTowerVariableSpeed] the cooling tower
+  # @param cooling_tower_variable_speed [OpenStudio::Model::CoolingTowerVariableSpeed] variable speed cooling tower
+  # @return [Bool] returns true if successful, false if not
   def cooling_tower_variable_speed_apply_efficiency_and_curves(cooling_tower_variable_speed)
     cooling_tower_apply_minimum_power_per_flow(cooling_tower_variable_speed)
-
     cooling_tower_variable_speed.setCellControl('MaximalCell')
-
     return true
   end
 end

data/lib/openstudio-standards/standards/ashrae_90_1/ashrae_90_1_2019/ashrae_90_1_2019.FanVariableVolume.rb

@@ -5,7 +5,7 @@ class ASHRAE9012019 < ASHRAE901
   # Per 90.1-2019, table 6.5.3.2.1: the fan motor size for chiller-water
   # and evaporative cooling is 0.25 hp as of 1/1/2014 instead of 5 hp
   #
-  # @param fan_variable_volume [OpenStudio::Model::FanVariableVolume] the fan
+  # @param fan_variable_volume [OpenStudio::Model::FanVariableVolume] variable volume fan object
   # @return [Double] the limit, in horsepower. Return nil for no limit by default.
   def fan_variable_volume_part_load_fan_power_limitation_hp_limit(fan_variable_volume)
     hp_limit = case fan_variable_volume_cooling_system_type(fan_variable_volume)
@@ -24,7 +24,7 @@ class ASHRAE9012019 < ASHRAE901
   # Per 90.1-2019, table 6.5.3.2.1: the cooling capacity threshold is 75000
   # instead of 110000 as of 1/1/2014
   #
-  # @param fan_variable_volume [OpenStudio::Model::FanVariableVolume] the fan
+  # @param fan_variable_volume [OpenStudio::Model::FanVariableVolume] variable volume fan object
   # @return [Double] the limit, in Btu/hr. Return nil for no limit by default.
   def fan_variable_volume_part_load_fan_power_limitation_capacity_limit(fan_variable_volume)
     cap_limit_btu_per_hr = case fan_variable_volume_cooling_system_type(fan_variable_volume)

data/lib/openstudio-standards/standards/ashrae_90_1/ashrae_90_1_2019/ashrae_90_1_2019.Space.rb

@@ -96,15 +96,17 @@ class ASHRAE9012019 < ASHRAE901
     return [req_top_ctrl, req_pri_ctrl, req_sec_ctrl]
   end
 
-  # Determine the fraction controlled by each sensor and which
-  # window each sensor should go near.
+  # Determine the fraction controlled by each sensor and which window each sensor should go near.
   #
-  # @param space [OpenStudio::Model::Space] the space with the daylighting
+  # @param space [OpenStudio::Model::Space] space object
+  # @param areas [Hash] a hash of daylighted areas
   # @param sorted_windows [Hash] a hash of windows, sorted by priority
   # @param sorted_skylights [Hash] a hash of skylights, sorted by priority
   # @param req_top_ctrl [Bool] if toplighting controls are required
   # @param req_pri_ctrl [Bool] if primary sidelighting controls are required
   # @param req_sec_ctrl [Bool] if secondary sidelighting controls are required
+  # @return [Array] array of 4 items
+  #   [sensor 1 fraction, sensor 2 fraction, sensor 1 window, sensor 2 window]
   def space_daylighting_fractions_and_windows(space,
                                               areas,
                                               sorted_windows,
@@ -160,7 +162,8 @@ class ASHRAE9012019 < ASHRAE901
           sensor_1_window = sorted_skylights[0]
         else
           # Sensor 1 controls toplighted area
-          sensor_1_frac = areas['toplighted_area'] / space_area_m2
+          num_sensors = 2
+          sensor_1_frac = areas['toplighted_area'] / space_area_m2 / num_sensors
           sensor_1_window = sorted_skylights[0]
           sensor_2_frac = sensor_1_frac
           sensor_2_window = sensor_1_window
@@ -178,7 +181,7 @@ class ASHRAE9012019 < ASHRAE901
     return [sensor_1_frac, sensor_2_frac, sensor_1_window, sensor_2_window]
   end
 
-  # Determine the base infiltration rate at 75 PA.
+  # Determine the base infiltration rate at 75 Pa.
   #
   # @return [Double] the baseline infiltration rate, in cfm/ft^2
   # defaults to no infiltration.
@@ -186,4 +189,21 @@ class ASHRAE9012019 < ASHRAE901
     basic_infil_rate_cfm_per_ft2 = 1.0
     return basic_infil_rate_cfm_per_ft2
   end
+
+  # Provide the type of daylighting control type
+  #
+  # @param [OpenStudio::Model::Space] OpenStudio Space object
+  # return [String] daylighting control type
+  def space_daylighting_control_type(space)
+    return 'ContinuousOff'
+  end
+
+  # Provide the minimum input power fraction for continuous
+  # dimming daylighting control
+  #
+  # @param [OpenStudio::Model::Space] OpenStudio Space object
+  # return [Float] daylighting minimum input power fraction
+  def space_daylighting_minimum_input_power_fraction(space)
+    return 0.2
+  end
 end

data/lib/openstudio-standards/standards/ashrae_90_1/ashrae_90_1_2019/ashrae_90_1_2019.ThermalZone.rb

@@ -24,7 +24,8 @@ class ASHRAE9012019 < ASHRAE901
         case bldg_type
         when 'Retail', 'StripMall', 'SuperMarket'
           occ_type = 'retail'
-          # when 'SomeBuildingType' # TODO add publicassembly building types
+          # when 'SomeBuildingType'
+          # @todo add publicassembly building types
           # occ_type = 'publicassembly'
         end
       end