Mech vent people list for CO2 control

Author: mjwitte

src/EnergyPlus/MixedAir.cc

@@ -2761,6 +2761,18 @@ void InitOAController(EnergyPlusData &state, int const OAControllerNum, bool con
                 }
             }
 
+            // Fill People index lists if needed
+            if (vent_mech.SystemOAMethod == DataSizing::SysOAMethod::ProportionalControlDesOcc) {
+                for (int peopleNum = 1; peopleNum <= state.dataHeatBal->TotPeople; ++peopleNum) {
+                    for (auto &thisMechVentZone : vent_mech.VentMechZone) {
+                        if (state.dataHeatBal->People(peopleNum).ZonePtr == thisMechVentZone.zoneNum) {
+                            thisMechVentZone.peopleIndexes.push_back(peopleNum);
+                            break;
+                        }
+                    }
+                }
+            }
+
             // Check to see if any zones on an air loop are not accounted for by a mechanical ventilation object
             for (int AirLoopZoneInfoZoneNum = 1; AirLoopZoneInfoZoneNum <= state.dataAirLoop->AirLoopZoneInfo(AirLoopNum).NumZones;
                  ++AirLoopZoneInfoZoneNum) {
@@ -3608,19 +3620,18 @@ Real64 VentilationMechanicalProps::CalcMechVentController(EnergyPlusData &state,
     // new local variables for DCV
     // Zone OA flow rate based on each calculation method [m3/s]
     std::array<Real64, static_cast<int>(DataSizing::OAFlowCalcMethod::Num)> ZoneOACalc{0.0};
-    Real64 CO2PeopleGeneration = 0; // CO2 generation from people at design level
-    Real64 ZoneOABZ;                // Zone breathing-zone OA flow rate [m3/s]
-    Real64 ZoneOA;                  // Zone OA flow rate [m3/s]
-    Real64 ZoneOAFrac;              // Zone OA fraction (as a fraction of actual supply air flow rate)
-    Real64 SysOAuc;                 // System uncorrected OA flow rate
-    Real64 SysOA;                   // System supply OA volume flow rate [m3/s]
-    Real64 SysEv;                   // System ventilation efficiency
-    Real64 NodeTemp;                // node temperature
-    Real64 NodeHumRat;              // node humidity ratio
-    Real64 ZoneMaxCO2 = 0.0;        // Breathing-zone CO2 concentartion
-    Real64 ZoneMinCO2 = 0.0;        // Minimum CO2 concentration in zone
-    Real64 ZoneOAMin = 0.0;         // Minimum Zone OA flow rate when the zone is unoccupied (i.e. ZoneOAPeople = 0)
-    Real64 ZoneOAMax = 0.0;         // Maximum Zone OA flow rate (ZoneOAPeople + ZoneOACalc[static_cast<int>(DataSizing::OAFlowCalcMethod::PerArea)])
+    Real64 ZoneOABZ;         // Zone breathing-zone OA flow rate [m3/s]
+    Real64 ZoneOA;           // Zone OA flow rate [m3/s]
+    Real64 ZoneOAFrac;       // Zone OA fraction (as a fraction of actual supply air flow rate)
+    Real64 SysOAuc;          // System uncorrected OA flow rate
+    Real64 SysOA;            // System supply OA volume flow rate [m3/s]
+    Real64 SysEv;            // System ventilation efficiency
+    Real64 NodeTemp;         // node temperature
+    Real64 NodeHumRat;       // node humidity ratio
+    Real64 ZoneMaxCO2 = 0.0; // Breathing-zone CO2 concentartion
+    Real64 ZoneMinCO2 = 0.0; // Minimum CO2 concentration in zone
+    Real64 ZoneOAMin = 0.0;  // Minimum Zone OA flow rate when the zone is unoccupied (i.e. ZoneOAPeople = 0)
+    Real64 ZoneOAMax = 0.0;  // Maximum Zone OA flow rate (ZoneOAPeople + ZoneOACalc[static_cast<int>(DataSizing::OAFlowCalcMethod::PerArea)])
     Real64 MechVentOAMassFlow = 0.0;
 
     // Apply mechanical ventilation only when it is available/allowed
@@ -3882,8 +3893,8 @@ Real64 VentilationMechanicalProps::CalcMechVentController(EnergyPlusData &state,
                                         // Calculate zone maximum target CO2 concentration in PPM
                                         if (this->SystemOAMethod == DataSizing::SysOAMethod::ProportionalControlDesOcc) {
                                             // Accumulate CO2 generation from people at design occupancy and current activity level
-                                            for (int PeopleNum = 1; PeopleNum <= state.dataHeatBal->TotPeople; ++PeopleNum) {
-                                                if (state.dataHeatBal->People(PeopleNum).ZonePtr != ZoneNum) continue;
+                                            Real64 CO2PeopleGeneration = 0.0;
+                                            for (int const PeopleNum : thisMechVentZone.peopleIndexes) {
                                                 CO2PeopleGeneration +=
                                                     state.dataHeatBal->People(PeopleNum).NumberOfPeople *
                                                     state.dataHeatBal->People(PeopleNum).CO2RateFactor *

src/EnergyPlus/MixedAir.hh

@@ -284,6 +284,7 @@ namespace MixedAir {
         DataSizing::OAFlowCalcMethod ZoneOAFlowMethod = DataSizing::OAFlowCalcMethod::PerPerson; // OA flow method for each zone
         int ZoneOASchPtr = 0;              // Index to the outdoor air schedule for each zone (from DesignSpecification:OutdoorAir or default)
         Real64 OAPropCtlMinRateSchPtr = 0; // Outdoor design OA flow rate schedule from DesignSpecification:OutdoorAir
+        EPVector<int> peopleIndexes; // List of People objects in this zone (for SystemOAMethod == DataSizing::SysOAMethod::ProportionalControlDesOcc)
     };
 
     struct VentilationMechanicalProps // Derived type for Ventilation:Mechanical data