modeling coefficient of performance curves for ASHP

edit

Hi,

I am a university student working on the electrification of heating systems. I want to model the coefficient of performance curves for ASHPs, and I was looking at how ResStock does it. Each of the air source heat pumps uses the corresponding performance curve for EIR based on how many stages it has. Then it is scaled based on some data on COP values.

Here is an excerpt from the code. Could anyone help me understand how to use these coefficients--what would be the variables associated with them or point me towards relevant documentation that might aid my understanding? I could not find their source:

require_relative "constants" require_relative "geometry" require_relative "util" require_relative "unit_conversions" require_relative "psychrometrics" require_relative "schedules"

class HVAC def self.apply_central_ac_1speed(model, unit, runner, seer, eers, shrs, fan_power_rated, fan_power_installed, crankcase_capacity, crankcase_temp, eer_capacity_derates, capacity, dse, frac_cool_load_served)

return true if frac_cool_load_served <= 0

num_speeds = 1

# Performance curves
# NOTE: These coefficients are in IP UNITS
curves_in_ip = true
cOOL_CAP_FT_SPEC = [[3.670270705, -0.098652414, 0.000955906, 0.006552414, -0.0000156, -0.000131877]]
cOOL_EIR_FT_SPEC = [[-3.302695861, 0.137871531, -0.001056996, -0.012573945, 0.000214638, -0.000145054]]
cOOL_CAP_FFLOW_SPEC = [[0.718605468, 0.410099989, -0.128705457]]
cOOL_EIR_FFLOW_SPEC = [[1.32299905, -0.477711207, 0.154712157]]

capacity_ratios = [1.0]
fan_speed_ratios = [1.0]

# Cooling Coil
rated_airflow_rate = 386.1 # cfm
cfms_ton_rated = calc_cfms_ton_rated(rated_airflow_rate, fan_speed_ratios, capacity_ratios)
cooling_eirs = calc_cooling_eirs(num_speeds, eers, fan_power_rated)
shrs_rated_gross = calc_shrs_rated_gross(num_speeds, shrs, fan_power_rated, cfms_ton_rated)
cOOL_CLOSS_FPLR_SPEC = [calc_plr_coefficients_cooling(num_speeds, seer)]

obj_name = Constants.ObjectNameCentralAirConditioner(unit.name.to_s)

thermal_zones = Geometry.get_thermal_zones_from_spaces(unit.spaces)

control_slave_zones_hash = get_control_and_slave_zones(thermal_zones)
control_slave_zones_hash.each do |control_zone, slave_zones|
  # _processCurvesDXCooling

  clg_coil_stage_data = calc_coil_stage_data_cooling(model, capacity, (0...num_speeds).to_a, cooling_eirs, shrs_rated_gross, cOOL_CAP_FT_SPEC, cOOL_EIR_FT_SPEC, cOOL_CLOSS_FPLR_SPEC, cOOL_CAP_FFLOW_SPEC, cOOL_EIR_FFLOW_SPEC, curves_in_ip, dse)


  # _processSystemFan

  fan = OpenStudio::Model::FanOnOff.new(model, model.alwaysOnDiscreteSchedule)
  fan_eff = 0.75 # Overall Efficiency of the Fan, Motor and Drive
  fan.setName(obj_name + " clg supply fan")
  fan.setEndUseSubcategory(obj_name + " clg supply fan")
  fan.setFanEfficiency(fan_eff)
  fan.setPressureRise(calculate_fan_pressure_rise(fan_eff, fan_power_installed / dse))
  fan.setMotorEfficiency(1.0)
  fan.setMotorInAirstreamFraction(1.0)

  # _processSystemAir

  air_loop_unitary = OpenStudio::Model::AirLoopHVACUnitarySystem.new(model)
  air_loop_unitary.setName(obj_name + " unitary system")
  air_loop_unitary.setAvailabilitySchedule(model.alwaysOnDiscreteSchedule)
  air_loop_unitary.setCoolingCoil(clg_coil)
  air_loop_unitary.setSupplyAirFlowRateDuringHeatingOperation(0.0)
  air_loop_unitary.setSupplyFan(fan)
  air_loop_unitary.setFanPlacement("BlowThrough")
  air_loop_unitary.setSupplyAirFanOperatingModeSchedule(model.alwaysOffDiscreteSchedule)
  air_loop_unitary.setMaximumSupplyAirTemperature(UnitConversions.convert(120.0, "F", "C"))
  air_loop_unitary.setSupplyAirFlowRateWhenNoCoolingorHeatingisRequired(0)

  air_loop = OpenStudio::Model::AirLoopHVAC.new(model)
  air_loop.setName(obj_name + " asys")
  air_supply_inlet_node = air_loop.supplyInletNode
  air_supply_outlet_node = air_loop.supplyOutletNode
  air_demand_inlet_node = air_loop.demandInletNode
  air_demand_outlet_node = air_loop.demandOutletNode

  air_loop_unitary.addToNode(air_supply_inlet_node)

  runner.registerInfo("Added '#{fan.name}' to '#{air_loop_unitary.name}' of '#{air_loop.name}'")
  runner.registerInfo("Added '#{clg_coil.name}' to '#{air_loop_unitary.name}' of '#{air_loop.name}'")

  air_loop_unitary.setControllingZoneorThermostatLocation(control_zone)

  # _processSystemDemandSideAir
  # Demand Side

  # Supply Air
  zone_splitter = air_loop.zoneSplitter
  zone_splitter.setName(obj_name + " zone splitter")

  zone_mixer = air_loop.zoneMixer
  zone_mixer.setName(obj_name + " zone mixer")

  diffuser_living = OpenStudio::Model::AirTerminalSingleDuctUncontrolled.new(model, model.alwaysOnDiscreteSchedule)
  diffuser_living.setName(obj_name + " #{control_zone.name} direct air")
  air_loop.multiAddBranchForZone(control_zone, diffuser_living)

  air_loop.multiAddBranchForZone(control_zone)
  runner.registerInfo("Added '#{air_loop.name}' to '#{control_zone.name}' of #{unit.name}")

  prioritize_zone_hvac(model, runner, control_zone)

  slave_zones.each do |slave_zone|
    diffuser_fbsmt = OpenStudio::Model::AirTerminalSingleDuctUncontrolled.new(model, model.alwaysOnDiscreteSchedule)
    diffuser_fbsmt.setName(obj_name + " #{slave_zone.name} direct air")
    air_loop.multiAddBranchForZone(slave_zone, diffuser_fbsmt)

    air_loop.multiAddBranchForZone(slave_zone)
    runner.registerInfo("Added '#{air_loop.name}' to '#{slave_zone.name}' of ...