How to model defrost cycles in EnergyPlus

edit

Hello, I am quite new to EnergyPlus modeling and am working on a project where I am trying to understand the impact of urban microclimate-adjusted weather files, produced through the Urban Weather Generator tool, on air source heat pump performance. To start, I am using the DOE reference model for a midrise apartment building. In OpenStudio I used the measure add_apszhp_to_each_zone to swap out all existing HVAC systems for air source heat pumps. This works and I am able to run my simulations in EnergyPlus through the run_ep_from_python functions. I am able to retrieve the timeseries output variables I need, being: "Cooling:Electricity J" and "Heating:Electricity J" from the output .csv file. This is great, because I can see how annual hourly energy consumption varies for the same building model across various microclimate weather files (different locations within the same city).

However, to truly model how microclimate influences air source heat pump performance, I want to capture the role of defrost cycling and its impact on energy consumption. The logic here is that during the winter, different microclimate relative humiditiy levels, on top of other localized environmental features, will influence the operation of defrost cycling, impacting the accumulation of frost on the outdoor coil. What I would like to do is plot defrost cycle energy consumption and use throughout the winter, which to my understanding would require the output variables:

• 'Output:Variable,*,Heat Exchanger Defrost Time Fraction,hourly; !- HVAC Average []'
• 'Output:Variable,*,Heating Coil Defrost Electricity Rate,hourly; !- HVAC Average [W]'
• 'Output:Variable,*,Heating Coil Defrost Electricity Energy,hourly; !- HVAC Sum [J]'

which I see listed in my simulation output .rdd file.

Though when I try to actually plot these defrost output time series variables, from the output .eso file, converted to a .csv file, these output variables are not to be found, which I am confused about. For reference, here is how I am defining defrost cycling operation in my IDF for the object Coil:Heating:DX:SingleSpeed:

  Coil:Heating:DX:SingleSpeed,
G N1 Apartment ZN Coil Heating DX Single Speed,  !- Name
Always On Discrete,      !- Availability Schedule Name
Autosize,                !- Gross Rated Heating Capacity {W}
3.1,                     !- Gross Rated Heating COP {W/W}
Autosize,                !- Rated Air Flow Rate {m3/s}
773.3,                   !- 2017 Rated Supply Fan Power Per Volume Flow Rate {W/(m3/s)}
934.4,                   !- 2023 Rated Supply Fan Power Per Volume Flow Rate {W/(m3/s)}
Node 163,                !- Air Inlet Node Name
Node 164,                !- Air Outlet Node Name
Curve Cubic 34,          !- Heating Capacity Function of Temperature Curve Name
Curve Cubic 35,          !- Heating Capacity Function of Flow Fraction Curve Name
Curve Cubic 36,          !- Energy Input Ratio Function of Temperature Curve Name
Curve Quadratic 59,      !- Energy Input Ratio Function of Flow Fraction Curve Name
Curve Quadratic 60,      !- Part Load Fraction Correlation Curve Name
FlatDefrostEIRCurve,     !- Defrost Energy Input Ratio Function of Temperature Curve Name
-8,                      !- Minimum Outdoor Dry-Bulb Temperature for Compressor Operation {C}
,                        !- Outdoor Dry-Bulb Temperature to Turn On Compressor {C}
5,                       !- Maximum Outdoor Dry-Bulb Temperature for Defrost Operation {C ...