Fix VRF Heat Pump Total Heating Rate not match the sum of coil heating rate (with piping loss)

doc/input-output-reference/src/overview/group-variable-refrigerant-flow-equipment.tex

@@ -619,7 +619,7 @@ \subsubsection{Outputs}\label{outputs-039}
 
 \paragraph{VRF Heat Pump Total Heating Rate {[}W{]}}\label{vrf-heat-pump-total-heating-rate-w}
 
-This output field is the operating total heating capacity of the variable refrigerant flow heat pump in Watts. The capacity includes any degradation due to defrost mode. This value is calculated for each HVAC system time step being simulated, and the results are averaged for the time step being reported. This value should match the sum of the individual zone terminal unit output variables for Zone VRF Air Terminal Total Heating Rate.
+This output field is the operating total heating capacity of the variable refrigerant flow heat pump in Watts. The capacity includes any degradation due to defrost mode. This value is calculated for each HVAC system time step being simulated, and the results are averaged for the time step being reported. This value should match the sum of the individual zone terminal unit heating coil output variables for Heating Coil Heating Rate plus any piping loss.
 
 \paragraph{VRF Heat Pump Cooling Electricity Rate {[}W{]}}\label{vrf-heat-pump-cooling-electric-power-w}
 
@@ -1153,6 +1153,10 @@ \subsubsection{Outputs}
 
 Note: refer to the rdd file after a simulation for exact output variable names
 
+\paragraph{VRF Heat Pump Total Heating Rate {[}W{]}}\label{vrf-heat-pump-total-heating-rate-w-fluidTCtrl}
+
+This output field is the operating total heating capacity of the variable refrigerant flow heat pump in Watts. The capacity includes any degradation due to defrost mode. This value is calculated for each HVAC system time step being simulated, and the results are averaged for the time step being reported. This value should match the sum of the individual zone terminal unit heating coil output variables for Heating Coil Heating Rate plus any piping loss.
+
 \paragraph{VRF Heat Pump Compressor Rotating Speed {[}rev/min{]}}\label{vrf-heat-pump-compressor-rotating-speed-revmin}
 
 This output only applies for the VRF-FluidTCtrl model. This is the rotating speed of the compressor, which indicates the loading index.

src/EnergyPlus/DXCoils.cc

@@ -17483,12 +17483,20 @@ void ControlVRFIUCoil(EnergyPlusData &state,
         if (QCoilSenHeatingLoad > QinSenMin1) {
             // Modulate fan speed to meet room sensible load; SC is not updated
             FanSpdRatioMax = 1.0;
-            auto f = [QCoilSenHeatingLoad, Ts_1, Tin, Garate, BF](Real64 FanSpdRto) {
-                return FanSpdResidualHeat(FanSpdRto, QCoilSenHeatingLoad, Ts_1, Tin, Garate, BF);
+            Tout = Tin + (Ts_1 - Tin) * (1 - BF);
+            Real64 RatedAirMassFlowRate = state.dataDXCoils->DXCoil(CoilIndex).RatedAirMassFlowRate[0];
+            auto f = [QCoilSenHeatingLoad, RatedAirMassFlowRate, Tout, Tin, Win](Real64 FanSpdRto) {
+                return FanSpdResidualHeatUsingH(FanSpdRto, QCoilSenHeatingLoad, RatedAirMassFlowRate, Tout, Tin, Win);
             };
             General::SolveRoot(state, 1.0e-3, MaxIter, SolFla, Ratio1, f, FanSpdRatioMin, FanSpdRatioMax);
             // this will likely cause problems eventually, -1 and -2 mean different things
-            if (SolFla < 0) Ratio1 = FanSpdRatioMax; // over capacity
+            if (SolFla < 0) {
+                if (f(FanSpdRatioMin) <= 0) { // capacity <= demand
+                    Ratio1 = FanSpdRatioMax;  // over capacity
+                } else {                      // capacity > demand even for the minimum fan speed
+                    Ratio1 = FanSpdRatioMin;
+                }
+            }
             FanSpdRatio = Ratio1;
             CoilOnOffRatio = 1.0;
 
@@ -17769,6 +17777,27 @@ Real64 FanSpdResidualHeat(Real64 FanSpdRto, Real64 QCoilSenHeatingLoad, Real64 T
     return (TotCap - ZnSenLoad) / ZnSenLoad;
 }
 
+Real64 FanSpdResidualHeatUsingH(Real64 FanSpdRto, Real64 QCoilSenHeatingLoad, Real64 RatedAirMassFlowRate, Real64 Tout, Real64 Tin, Real64 Win)
+{
+
+    // FUNCTION INFORMATION:
+    //       AUTHOR         Yujie Xu (yujiex)
+    //       DATE WRITTEN   Jul 2024
+    //
+    // PURPOSE OF THIS FUNCTION:
+    //       Calculates residual function (desired zone heating load - actual heating coil capacity)
+    //       This is used to modify the fan speed to adjust the coil heating capacity to match
+    //       the zone heating load. This one uses Hin and Hout difference rather than Tin and Tout difference
+    //       like in FanSpdResidualHeat
+    //
+    Real64 ZnSenLoad = QCoilSenHeatingLoad;
+    // +-100 W minimum zone load?
+    if (std::abs(ZnSenLoad) < 100.0) ZnSenLoad = sign(100.0, ZnSenLoad);
+    Real64 Wout = Win;
+    Real64 TotCap = FanSpdRto * RatedAirMassFlowRate * (PsyHFnTdbW(Tout, Wout) - PsyHFnTdbW(Tin, Win));
+    return (TotCap - ZnSenLoad) / ZnSenLoad;
+}
+
 void SetMSHPDXCoilHeatRecoveryFlag(EnergyPlusData &state, int const DXCoilNum)
 {
 

src/EnergyPlus/DXCoils.hh

@@ -932,6 +932,8 @@ namespace DXCoils {
 
     Real64 FanSpdResidualHeat(Real64 FanSpdRto, Real64 QCoilSenHeatingLoad, Real64 Ts_1, Real64 Tin, Real64 Garate, Real64 BF);
 
+    Real64 FanSpdResidualHeatUsingH(Real64 FanSpdRto, Real64 QCoilSenHeatingLoad, Real64 RatedAirMassFlowRate, Real64 Tout, Real64 Tin, Real64 Win);
+
     void SetMSHPDXCoilHeatRecoveryFlag(EnergyPlusData &state, int const DXCoilNum); // must match coil names for the coil type
 
     void SetDXCoilAirLoopNumber(EnergyPlusData &state, std::string const &CoilName, int const AirLoopNum); // must match coil names for the coil type

src/EnergyPlus/HVACVariableRefrigerantFlow.cc

@@ -11725,12 +11725,14 @@ void VRFCondenserEquipment::CalcVRFCondenser_FluidTCtrl(EnergyPlusData &state)
         this->VRFCondCyclingRatio = CyclingRatio;
 
         Tsuction = this->EvaporatingTemp; // Outdoor unit evaporating temperature
+        this->HeatingCapacityPrev = this->HeatingCapacity;
         this->HeatingCapacity =
             this->CoffEvapCap * this->RatedEvapCapacity * CurveValue(state, this->OUCoolingCAPFT(NumOfCompSpdInput), Tdischarge, Tsuction) +
             this->RatedCompPower * CurveValue(state,
                                               this->OUCoolingPWRFT(NumOfCompSpdInput),
                                               Tdischarge,
                                               Tsuction); // Include the piping loss, at the highest compressor speed
+        this->PipingCorrectionHeatingPrev = this->PipingCorrectionHeating;
         this->PipingCorrectionHeating = TU_HeatingLoad / (TU_HeatingLoad + Pipe_Q_h);
         state.dataHVACVarRefFlow->MaxHeatingCapacity(VRFCond) =
             this->HeatingCapacity; // for report, maximum condensing capacity the system can provide
@@ -12219,7 +12221,14 @@ void VRFCondenserEquipment::CalcVRFCondenser_FluidTCtrl(EnergyPlusData &state)
     }
 
     this->TotalCoolingCapacity = TotalCondCoolingCapacity * CoolingPLR;
-    this->TotalHeatingCapacity = TotalCondHeatingCapacity * HeatingPLR;
+    // adjustment for matching HP heating rate and coil heating rate
+    this->TotalHeatingCapacity = TotalCondHeatingCapacity * HeatingPLR * (this->RatedEvapCapacity / (this->RatedEvapCapacity + Pipe_Q_h));
+    if (this->VRFCondPLR < 1.0) {
+        this->TotalHeatingCapacity = TotalCondHeatingCapacity * HeatingPLR * this->PipingCorrectionHeating;
+    }
+    if (this->TUHeatingLoad / this->PipingCorrectionHeating > TotalCondHeatingCapacity) {
+        this->TotalHeatingCapacity = this->HeatingCapacityPrev * HeatingPLR * this->PipingCorrectionHeatingPrev;
+    }
 
     if (this->MinPLR > 0.0) {
         bool const plrTooLow = this->VRFCondPLR < this->MinPLR;

src/EnergyPlus/HVACVariableRefrigerantFlow.hh

@@ -173,6 +173,7 @@ namespace HVACVariableRefrigerantFlow {
         Real64 OperatingCoolingCOP;            // Operating VRF heat pump cooling COP (W/W)
         Real64 RatedCoolingPower;              // Rated cooling power = Rated Cooling Capacity / Rated COP (W)
         Real64 HeatingCapacity;                // Nominal VRF heat pump heating capacity (W)
+        Real64 HeatingCapacityPrev;            // Nominal VRF heat pump heating capacity (W)
         Real64 HeatingCapacitySizeRatio;       // Ratio of heating to cooling when autosizing
         bool LockHeatingCapacity;              // used in sizing to size VRF heat cap to VRF cool cap
         Real64 TotalHeatingCapacity;           // Nominal VRF heat pump heating capacity (W)
@@ -219,6 +220,7 @@ namespace HVACVariableRefrigerantFlow {
         Real64 PCFHeightHeat;                                             // piping correction factor for height in heating mode
         Real64 EquivPipeLngthHeat;                                        // equivalent piping length for heating
         Real64 PipingCorrectionHeating;                                   // piping correction factor for heating
+        Real64 PipingCorrectionHeatingPrev;                               // piping correction factor for heating
         Real64 CCHeaterPower;                                             // crankcase heater power per compressor (W)
         Real64 CompressorSizeRatio;                                       // ratio of min compressor size to total capacity
         int NumCompressors;                                               // number of compressors in VRF condenser
@@ -394,18 +396,18 @@ namespace HVACVariableRefrigerantFlow {
               WaterCondenserDesignMassFlow(0.0), WaterCondenserMassFlow(0.0), QCondenser(0.0), QCondEnergy(0.0), CondenserSideOutletTemp(0.0),
               SchedPtr(-1), CoolingCapacity(0.0), TotalCoolingCapacity(0.0), CoolingCombinationRatio(1.0), VRFCondPLR(0.0), VRFCondRTF(0.0),
               VRFCondCyclingRatio(0.0), CondenserInletTemp(0.0), CoolingCOP(0.0), OperatingCoolingCOP(0.0), RatedCoolingPower(0.0),
-              HeatingCapacity(0.0), HeatingCapacitySizeRatio(1.0), LockHeatingCapacity(false), TotalHeatingCapacity(0.0),
+              HeatingCapacity(0.0), HeatingCapacityPrev(0.0), HeatingCapacitySizeRatio(1.0), LockHeatingCapacity(false), TotalHeatingCapacity(0.0),
               HeatingCombinationRatio(1.0), HeatingCOP(0.0), OperatingHeatingCOP(0.0), RatedHeatingPower(0.0), MinOATCooling(0.0), MaxOATCooling(0.0),
               MinOATHeating(0.0), MaxOATHeating(0.0), CoolCapFT(0), CoolEIRFT(0), HeatCapFT(0), HeatEIRFT(0), CoolBoundaryCurvePtr(0),
               HeatBoundaryCurvePtr(0), EIRCoolBoundaryCurvePtr(0), CoolEIRFPLR1(0), CoolEIRFPLR2(0), CoolCapFTHi(0), CoolEIRFTHi(0), HeatCapFTHi(0),
               HeatEIRFTHi(0), EIRHeatBoundaryCurvePtr(0), HeatEIRFPLR1(0), HeatEIRFPLR2(0), CoolPLFFPLR(0), HeatPLFFPLR(0), MinPLR(0.0),
               MasterZonePtr(0), MasterZoneTUIndex(0), ThermostatPriority(ThermostatCtrlType::Invalid), SchedPriorityPtr(0), ZoneTUListPtr(0),
               HeatRecoveryUsed(false), VertPipeLngth(0.0), PCFLengthCoolPtr(0), PCFHeightCool(0.0), EquivPipeLngthCool(0.0),
               PipingCorrectionCooling(1.0), PCFLengthHeatPtr(0), PCFHeightHeat(0.0), EquivPipeLngthHeat(0.0), PipingCorrectionHeating(1.0),
-              CCHeaterPower(0.0), CompressorSizeRatio(0.0), NumCompressors(0), MaxOATCCHeater(0.0), DefrostEIRPtr(0), DefrostFraction(0.0),
-              DefrostStrategy(StandardRatings::DefrostStrat::Invalid), DefrostControl(StandardRatings::HPdefrostControl::Invalid),
-              DefrostCapacity(0.0), DefrostPower(0.0), DefrostConsumption(0.0), MaxOATDefrost(0.0),
-              CondenserType(DataHeatBalance::RefrigCondenserType::Invalid), CondenserNodeNum(0), SkipCondenserNodeNumCheck(false),
+              PipingCorrectionHeatingPrev(1.0), CCHeaterPower(0.0), CompressorSizeRatio(0.0), NumCompressors(0), MaxOATCCHeater(0.0),
+              DefrostEIRPtr(0), DefrostFraction(0.0), DefrostStrategy(StandardRatings::DefrostStrat::Invalid),
+              DefrostControl(StandardRatings::HPdefrostControl::Invalid), DefrostCapacity(0.0), DefrostPower(0.0), DefrostConsumption(0.0),
+              MaxOATDefrost(0.0), CondenserType(DataHeatBalance::RefrigCondenserType::Invalid), CondenserNodeNum(0), SkipCondenserNodeNumCheck(false),
               CondenserOutletNodeNum(0), WaterCondVolFlowRate(0.0), EvapCondEffectiveness(0.0), EvapCondAirVolFlowRate(0.0), EvapCondPumpPower(0.0),
               CoolCombRatioPTR(0), HeatCombRatioPTR(0), OperatingMode(0), ElecPower(0.0), ElecCoolingPower(0.0), ElecHeatingPower(0.0),
               CoolElecConsumption(0.0), HeatElecConsumption(0.0), CrankCaseHeaterPower(0.0), CrankCaseHeaterElecConsumption(0.0),

testfiles/US+SF+CZ4A+hp+crawlspace+IECC_2006_VRF.idf

@@ -477,8 +477,8 @@
     0,                       !- No Cooling Supply Air Flow Rate {m3/s}
     0.595,                   !- Heating Supply Air Flow Rate {m3/s}
     0,                       !- No Heating Supply Air Flow Rate {m3/s}
-    autosize,                !- Cooling Outdoor Air Flow Rate {m3/s}
-    autosize,                !- Heating Outdoor Air Flow Rate {m3/s}
+    0,                !- Cooling Outdoor Air Flow Rate {m3/s}
+    0,                !- Heating Outdoor Air Flow Rate {m3/s}
     0,                       !- No Load Outdoor Air Flow Rate {m3/s}
     VRFFanModeSchedule,      !- Supply Air Fan Operating Mode Schedule Name
     drawthrough,             !- Supply Air Fan Placement