openstudio-standards 0.1.13 → 0.1.14

data/lib/openstudio-standards/prototypes/Prototype.utilities.rb

@@ -306,3 +306,178 @@ def adjust_infiltration_to_prototype_building_conditions(initial_infiltration_ra
 
   return adjusted_infiltration_rate_m3_per_s
 end
+
+# Convert biquadratic curves that are a function of temperature
+# from IP (F) to SI (C) or vice-versa.  The curve is of the form
+# z = C1 + C2*x + C3*x^2 + C4*y + C5*y^2 + C6*x*y
+# where C1, C2, ... are the coefficients,
+# x is the first independent variable (in F or C)
+# y is the second independent variable (in F or C)
+# and z is the resulting value
+#
+# @author Scott Horowitz, NREL
+# @param coeffs [Array<Double>] an array of 6 coefficients, in order
+# @return [Array<Double>] the revised coefficients in the new unit system
+def convert_curve_biquadratic(coeffs, ip_to_si=true)
+  if ip_to_si
+    # Convert IP curves to SI curves
+    si_coeffs = []
+    si_coeffs << coeffs[0] + 32.0 * (coeffs[1] + coeffs[3]) + 1024.0 * (coeffs[2] + coeffs[4] + coeffs[5])
+    si_coeffs << 9.0 / 5.0 * coeffs[1] + 576.0 / 5.0 * coeffs[2] + 288.0 / 5.0 * coeffs[5]
+    si_coeffs << 81.0 / 25.0 * coeffs[2]
+    si_coeffs << 9.0 / 5.0 * coeffs[3] + 576.0 / 5.0 * coeffs[4] + 288.0 / 5.0 * coeffs[5]
+    si_coeffs << 81.0 / 25.0 * coeffs[4]
+    si_coeffs << 81.0 / 25.0 * coeffs[5]        
+    return si_coeffs
+  else
+    # Convert SI curves to IP curves
+    ip_coeffs = []
+    ip_coeffs << coeffs[0] - 160.0/9.0 * (coeffs[1] + coeffs[3]) + 25600.0/81.0 * (coeffs[2] + coeffs[4] + coeffs[5])
+    ip_coeffs << 5.0/9.0 * (coeffs[1] - 320.0/9.0 * coeffs[2] - 160.0/9.0 * coeffs[5])
+    ip_coeffs << 25.0/81.0 * coeffs[2]
+    ip_coeffs << 5.0/9.0 * (coeffs[3] - 320.0/9.0 * coeffs[4] - 160.0/9.0 * coeffs[5])
+    ip_coeffs << 25.0/81.0 * coeffs[4]
+    ip_coeffs << 25.0/81.0 * coeffs[5]
+    return ip_coeffs
+  end
+end
+
+# Create a biquadratic curve of the form
+# z = C1 + C2*x + C3*x^2 + C4*y + C5*y^2 + C6*x*y
+#
+# @author Scott Horowitz, NREL
+# @param coeffs [Array<Double>] an array of 6 coefficients, in order
+# @param crv_name [String] the name of the curve
+# @param minX [Double] the minimum value of independent variable X that will be used
+# @param maxX [Double] the maximum value of independent variable X that will be used
+# @param minY [Double] the minimum value of independent variable Y that will be used
+# @param maxY [Double] the maximum value of independent variable Y that will be used
+# @param minOut [Double] the minimum value of dependent variable Z
+# @param maxOut [Double] the maximum value of dependent variable Z
+def create_curve_biquadratic(coeffs, crv_name, minX, maxX, minY, maxY, minOut, maxOut)
+  curve = OpenStudio::Model::CurveBiquadratic.new(self)
+  curve.setName(crv_name)
+  curve.setCoefficient1Constant(coeffs[0])
+  curve.setCoefficient2x(coeffs[1])
+  curve.setCoefficient3xPOW2(coeffs[2])
+  curve.setCoefficient4y(coeffs[3])
+  curve.setCoefficient5yPOW2(coeffs[4])
+  curve.setCoefficient6xTIMESY(coeffs[5])
+  curve.setMinimumValueofx(minX) unless minX.nil?
+  curve.setMaximumValueofx(maxX) unless maxX.nil?
+  curve.setMinimumValueofy(minY) unless minY.nil?
+  curve.setMaximumValueofy(maxY) unless maxY.nil?
+  curve.setMinimumCurveOutput(minOut) unless minOut.nil?
+  curve.setMaximumCurveOutput(maxOut) unless maxOut.nil?
+  return curve
+end
+
+# Create a bicubic curve of the form
+# z = C1 + C2*x + C3*x^2 + C4*y + C5*y^2 + C6*x*y + C7*x^3 + C8*y^3 + C9*x^2*y + C10*x*y^2
+#
+# @author Scott Horowitz, NREL
+# @param coeffs [Array<Double>] an array of 10 coefficients, in order
+# @param crv_name [String] the name of the curve
+# @param minX [Double] the minimum value of independent variable X that will be used
+# @param maxX [Double] the maximum value of independent variable X that will be used
+# @param minY [Double] the minimum value of independent variable Y that will be used
+# @param maxY [Double] the maximum value of independent variable Y that will be used
+# @param minOut [Double] the minimum value of dependent variable Z
+# @param maxOut [Double] the maximum value of dependent variable Z
+def create_curve_bicubic(coeffs, crv_name, minX, maxX, minY, maxY, minOut, maxOut)
+  curve = OpenStudio::Model::CurveBicubic.new(self)
+  curve.setName(crv_name)
+  curve.setCoefficient1Constant(coeffs[0])
+  curve.setCoefficient2x(coeffs[1])
+  curve.setCoefficient3xPOW2(coeffs[2])
+  curve.setCoefficient4y(coeffs[3])
+  curve.setCoefficient5yPOW2(coeffs[4])
+  curve.setCoefficient6xTIMESY(coeffs[5])
+  curve.setCoefficient7xPOW3(coeffs[6])
+  curve.setCoefficient8yPOW3(coeffs[7])
+  curve.setCoefficient9xPOW2TIMESY(coeffs[8])
+  curve.setCoefficient10xTIMESYPOW2(coeffs[9])
+  curve.setMinimumValueofx(minX) unless minX.nil?
+  curve.setMaximumValueofx(maxX) unless maxX.nil?
+  curve.setMinimumValueofy(minY) unless minY.nil?
+  curve.setMaximumValueofy(maxY) unless maxY.nil?
+  curve.setMinimumCurveOutput(minOut) unless minOut.nil?
+  curve.setMaximumCurveOutput(maxOut) unless maxOut.nil?
+  return curve
+end
+
+# Create a quadratic curve of the form
+# z = C1 + C2*x + C3*x^2
+#
+# @author Scott Horowitz, NREL
+# @param coeffs [Array<Double>] an array of 3 coefficients, in order
+# @param crv_name [String] the name of the curve
+# @param minX [Double] the minimum value of independent variable X that will be used
+# @param maxX [Double] the maximum value of independent variable X that will be used
+# @param minOut [Double] the minimum value of dependent variable Z
+# @param maxOut [Double] the maximum value of dependent variable Z
+# @param is_dimensionless [Bool] if true, the X independent variable is considered unitless
+# and the resulting output dependent variable is considered unitless
+def create_curve_quadratic(coeffs, crv_name, minX, maxX, minOut, maxOut, is_dimensionless=false)
+  curve = OpenStudio::Model::CurveQuadratic.new(self)
+  curve.setName(crv_name)
+  curve.setCoefficient1Constant(coeffs[0])
+  curve.setCoefficient2x(coeffs[1])
+  curve.setCoefficient3xPOW2(coeffs[2])
+  curve.setMinimumValueofx(minX) unless minX.nil?
+  curve.setMaximumValueofx(maxX) unless maxX.nil?
+  curve.setMinimumCurveOutput(minOut) unless minOut.nil?
+  curve.setMaximumCurveOutput(maxOut) unless maxOut.nil?
+  if is_dimensionless
+    curve.setInputUnitTypeforX("Dimensionless")
+    curve.setOutputUnitType("Dimensionless")
+  end
+  return curve
+end
+
+# Create a cubic curve of the form
+# z = C1 + C2*x + C3*x^2 + C4*x^3
+#
+# @author Scott Horowitz, NREL
+# @param coeffs [Array<Double>] an array of 4 coefficients, in order
+# @param crv_name [String] the name of the curve
+# @param minX [Double] the minimum value of independent variable X that will be used
+# @param maxX [Double] the maximum value of independent variable X that will be used
+# @param minOut [Double] the minimum value of dependent variable Z
+# @param maxOut [Double] the maximum value of dependent variable Z
+def create_curve_cubic(coeffs, crv_name, minX, maxX, minOut, maxOut)    
+  curve = OpenStudio::Model::CurveCubic.new(self)
+  curve.setName(crv_name)
+  curve.setCoefficient1Constant(coeffs[0])
+  curve.setCoefficient2x(coeffs[1])
+  curve.setCoefficient3xPOW2(coeffs[2])
+  curve.setCoefficient4xPOW3(coeffs[3])
+  curve.setMinimumValueofx(minX) unless minX.nil?
+  curve.setMaximumValueofx(maxX) unless maxX.nil?
+  curve.setMinimumCurveOutput(minOut) unless minOut.nil?
+  curve.setMaximumCurveOutput(maxOut) unless maxOut.nil?
+  return curve
+end
+
+# Create an exponential curve of the form
+# z = C1 + C2*x^C3
+#
+# @author Scott Horowitz, NREL
+# @param coeffs [Array<Double>] an array of 3 coefficients, in order
+# @param crv_name [String] the name of the curve
+# @param minX [Double] the minimum value of independent variable X that will be used
+# @param maxX [Double] the maximum value of independent variable X that will be used
+# @param minOut [Double] the minimum value of dependent variable Z
+# @param maxOut [Double] the maximum value of dependent variable Z
+def create_curve_exponent(coeffs, crv_name, minX, maxX, minOut, maxOut)
+  curve = OpenStudio::Model::CurveExponent.new(self)
+  curve.setName(crv_name)
+  curve.setCoefficient1Constant(coeffs[0])
+  curve.setCoefficient2Constant(coeffs[1])
+  curve.setCoefficient3Constant(coeffs[2])
+  curve.setMinimumValueofx(minX) unless minX.nil?
+  curve.setMaximumValueofx(maxX) unless maxX.nil?
+  curve.setMinimumCurveOutput(minOut) unless minOut.nil?
+  curve.setMaximumCurveOutput(maxOut) unless maxOut.nil?
+  return curve
+end

data/lib/openstudio-standards/standards/Standards.Model.rb

@@ -1244,9 +1244,9 @@ class OpenStudio::Model::Model
         end
 
         # If electric zone heat
-        electric_reheat = false
+        reheat_type = 'Water'
         if zone_heat_fuel == 'Electricity'
-          electric_reheat = true
+          reheat_type = 'Electricity'
         end
 
         # Group zones by story
@@ -1288,7 +1288,7 @@ class OpenStudio::Model::Model
                            0.9,
                            OpenStudio.convert(4.0, 'inH_{2}O', 'Pa').get,
                            nil,
-                           electric_reheat,
+                           reheat_type,
                            nil)
           end
 

data/lib/openstudio-standards/standards/Standards.PlantLoop.rb

@@ -254,7 +254,7 @@ class OpenStudio::Model::PlantLoop
         # Determine the leaving CW temp
         max_leaving_cw_t_f = 85
         leaving_cw_t_10f_approach_f = design_oat_wb_f + 10
-        leaving_cw_t_f = [max_leaving_cw_t_f, leaving_cw_t_10f_approach_f].max
+        leaving_cw_t_f = [max_leaving_cw_t_f, leaving_cw_t_10f_approach_f].min
 
         # Calculate the approach
         approach_r = leaving_cw_t_f - design_oat_wb_f

data/lib/openstudio-standards/standards/Standards.ThermalZone.rb

@@ -42,7 +42,7 @@ class OpenStudio::Model::ThermalZone
       oa_for_people = number_of_people * dsn_oa.outdoorAirFlowperPerson
       oa_for_floor_area = floor_area * dsn_oa.outdoorAirFlowperFloorArea
       oa_rate = dsn_oa.outdoorAirFlowRate
-      oa_for_volume = volume * dsn_oa.outdoorAirFlowAirChangesperHour
+      oa_for_volume = volume * dsn_oa.outdoorAirFlowAirChangesperHour / 3600
 
       # First check if this space uses the Maximum method and other spaces do not
       if dsn_oa.outdoorAirMethod == 'Maximum'
@@ -1453,6 +1453,10 @@ class OpenStudio::Model::ThermalZone
       else
         sch_name = space_type_properties['exhaust_schedule']
         exhaust_schedule = model.add_schedule(sch_name)
+        unless exhaust_schedule
+          OpenStudio.logFree(OpenStudio::Warn, 'openstudio.Standards.ThermalZone', "Could not find an exhaust schedule called #{sch_name}, exhaust fans will run continuously.")
+          exhaust_schedule = model.alwaysOnDiscreteSchedule
+        end
       end
 
       # add exhaust fans

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

@@ -0,0 +1,174 @@
+# Adds sql queries to the model
+class OpenStudio::Model::Model
+
+  def sql_file
+
+    # Ensure that the model has a sql file associated with it
+    if self.sqlFile.empty?
+      OpenStudio::logFree(OpenStudio::Error, 'openstudio.model.Model', 'Failed to retrieve data because the sql file containing results is missing.')
+      return false
+    end
+
+    return self.sqlFile.get
+  end
+
+  def annual_occupied_unmet_hours
+
+    sql = self.sql_file
+
+    # setup the queries
+    heating_setpoint_unmet_query = "SELECT Value
+                                    FROM TabularDataWithStrings
+                                    WHERE ReportName='SystemSummary'
+                                    AND ReportForString='Entire Facility'
+                                    AND TableName='Time Setpoint Not Met'
+                                    AND RowName = 'Facility'
+                                    AND ColumnName='During Occupied Heating'"
+
+    cooling_setpoint_unmet_query = "SELECT Value
+                                    FROM TabularDataWithStrings
+                                    WHERE ReportName='SystemSummary'
+                                    AND ReportForString='Entire Facility'
+                                    AND TableName='Time Setpoint Not Met'
+                                    AND RowName = 'Facility'
+                                    AND ColumnName='During Occupied Cooling'"
+    
+    # get the info
+    heating_setpoint_unmet = sql.execAndReturnFirstDouble(heating_setpoint_unmet_query)
+    cooling_setpoint_unmet = sql.execAndReturnFirstDouble(cooling_setpoint_unmet_query)
+    
+    # make sure all the data are availalbe
+    if heating_setpoint_unmet.empty? || cooling_setpoint_unmet.empty?
+      OpenStudio::logFree(OpenStudio::Error, 'openstudio.model.Model', 'Could not get unmet hours information.')
+      return false
+    end
+    
+    #aggregate heating and cooling hrs
+    heating_or_cooling_setpoint_unmet = heating_setpoint_unmet.get + cooling_setpoint_unmet.get    
+ 
+    return heating_or_cooling_setpoint_unmet
+  end
+
+  def annual_occupied_unmet_heating_hours
+
+    sql = self.sql_file
+
+    # setup the queries
+    heating_setpoint_unmet_query = "SELECT Value
+                                    FROM TabularDataWithStrings
+                                    WHERE ReportName='SystemSummary'
+                                    AND ReportForString='Entire Facility'
+                                    AND TableName='Time Setpoint Not Met'
+                                    AND RowName = 'Facility'
+                                    AND ColumnName='During Occupied Heating'"
+
+    # get the info
+    heating_setpoint_unmet = sql.execAndReturnFirstDouble(heating_setpoint_unmet_query)
+    
+    # make sure all the data are availalbe
+    if heating_setpoint_unmet.empty?
+      OpenStudio::logFree(OpenStudio::Error, 'openstudio.model.Model', 'Could not get unmet heating hours information.')
+      return false
+    end
+
+    return heating_setpoint_unmet.get
+  end
+
+  def annual_occupied_unmet_cooling_hours
+
+    sql = self.sql_file
+
+    # setup the queries
+    cooling_setpoint_unmet_query = "SELECT Value
+                                    FROM TabularDataWithStrings
+                                    WHERE ReportName='SystemSummary'
+                                    AND ReportForString='Entire Facility'
+                                    AND TableName='Time Setpoint Not Met'
+                                    AND RowName = 'Facility'
+                                    AND ColumnName='During Occupied Cooling'"
+
+    # get the info
+    cooling_setpoint_unmet = sql.execAndReturnFirstDouble(cooling_setpoint_unmet_query)
+    
+    # make sure all the data are availalbe
+    if cooling_setpoint_unmet.empty?
+      OpenStudio::logFree(OpenStudio::Error, 'openstudio.model.Model', 'Could not get unmet cooling hours information.')
+      return false
+    end
+
+    return cooling_setpoint_unmet.get
+  end
+
+  def annual_eui_kbtu_per_ft2
+
+    sql = self.sql_file
+
+    building = self.getBuilding
+    
+    # make sure all required data are available
+    if sql.totalSiteEnergy.empty?
+      OpenStudio::logFree(OpenStudio::Error, 'openstudio.model.Model', 'Site energy data unavailable.')
+      return false
+    end
+    
+    total_site_energy_kBtu = OpenStudio::convert(sql.totalSiteEnergy.get, "GJ", "kBtu").get
+  
+    floor_area_ft2 = OpenStudio::convert(building.floorArea, "m^2", "ft^2").get
+    
+    site_eui_kbtu_per_ft2 = total_site_energy_kBtu / floor_area_ft2
+
+    return site_eui_kbtu_per_ft2
+  end
+
+  def net_conditioned_floor_area
+
+    sql = self.sql_file
+
+    # setup the queries
+    area_query = "SELECT Value
+                  FROM TabularDataWithStrings
+                  WHERE ReportName='AnnualBuildingUtilityPerformanceSummary'
+                  AND ReportForString='Entire Facility'
+                  AND TableName='Building Area'
+                  AND RowName = 'Net Conditioned Building Area'
+                  AND ColumnName='Area'"
+
+    # get the info
+    area_m2 = sql.execAndReturnFirstDouble(area_query)
+    
+    # make sure all the data are availalbe
+    if area_m2.empty?
+      OpenStudio::logFree(OpenStudio::Error, 'openstudio.model.Model', 'Could not get conditioned area information.')
+      return false
+    end
+
+    return area_m2.get
+  end
+
+  def annual_energy_by_fuel_and_enduse(fuel_type, end_use)
+
+    sql = self.sql_file
+
+    # setup the queries
+    query = "SELECT Value
+             FROM TabularDataWithStrings
+             WHERE ReportName='AnnualBuildingUtilityPerformanceSummary'
+             AND ReportForString='Entire Facility'
+             AND TableName='End Uses'
+             AND RowName = '#{end_use}'
+             AND ColumnName='#{fuel_type}'"
+
+    # get the info
+    energy_gj = sql.execAndReturnFirstDouble(query)
+    
+    # make sure all the data are availalbe
+    if energy_gj.empty?
+      OpenStudio::logFree(OpenStudio::Error, 'openstudio.model.Model', "Could not get energy for #{fuel_type} #{end_use}.")
+      return 0.0
+    end
+
+    return energy_gj.get
+
+  end
+
+end

data/lib/openstudio-standards/version.rb

@@ -1,3 +1,3 @@
 module OpenstudioStandards
-  VERSION = '0.1.13'.freeze
+  VERSION = '0.1.14'.freeze
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: openstudio-standards
 version: !ruby/object:Gem::Version
-  version: 0.1.13
+  version: 0.1.14
 platform: ruby
 authors:
 - Andrew Parker
@@ -21,7 +21,7 @@ authors:
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2017-03-24 00:00:00.000000000 Z
+date: 2017-06-15 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: nokogiri
@@ -733,6 +733,7 @@ files:
 - lib/openstudio-standards/utilities/hash.rb
 - lib/openstudio-standards/utilities/logging.rb
 - lib/openstudio-standards/utilities/simulation.rb
+- lib/openstudio-standards/utilities/sqlfile.rb
 - lib/openstudio-standards/version.rb
 - lib/openstudio-standards/weather/Weather.Model.rb
 - lib/openstudio-standards/weather/Weather.stat_file.rb