Journal of Energy

Volume 2014 (2014), Article ID 670236, 6 pages

http://dx.doi.org/10.1155/2014/670236

## Determining Infiltration Rates and Predicting Building Occupancy Using CO_{2} Concentration Curves

Faculty of Engineering and Applied Science, Memorial University of Newfoundland, St. John’s, NL, Canada A1C 5S7

Received 29 August 2013; Accepted 10 December 2013; Published 6 February 2014

Academic Editor: Harry Boyer

Copyright © 2014 P. Parsons. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Demand controlled ventilation (DCV) reduces energy loss by reducing the air exchange flow rate to the minimum required to maintain acceptable indoor air quality (IAQ). DCV commonly uses carbon dioxide (CO_{2}) as a proxy for human activity and increases the ventilation rate once a preset CO_{2} threshold is exceeded. Significant improvements over threshold based ODV strategies are possible if the natural infiltration rate of the building is measured and the occupancy schedule determined by analysing the CO_{2} concentration continuously. These calculated parameters allow mathematical modeling of the ventilated space and the determination of future CO_{2} concentrations and allow prediction of future ventilation demands. The natural infiltration rate and the onset and duration of vacancy periods in a residential dwelling were determined by analysing CO_{2} concentration data. Concentration declines which fit an exponential decay curve with a correlation coefficient >0.90 identified all vacancy periods. The measured natural infiltration rate was found statistically correlated with average wind speed. A dynamic predicted occupancy map was constructed that has the potential to facilitate significant energy savings via deferred ventilation and intelligent cooling and heating strategies.

#### 1. Introduction

The energy required for heating, ventilation, and air conditioning (HVAC) of buildings can be reduced by minimizing uncontrolled air infiltration through leaks in the building envelop. The natural air infiltration rate of older North American residential structures is large enough that IAQ is seldom a problem. However, the continuous air and vapour barrier installed in new construction has necessitated the introduction of mechanical ventilation systems to remove excess moisture and maintain acceptable IAQ. This mechanical ventilation is often provided by a heat or energy recovery ventilator (HRV, ERV) which is sized based on the number and types of rooms and provides a fixed quantity of air exchanges per hour (ACH). This static approach results in excessive ventilation and energy waste during periods of low occupancy.

DCV systems reduce these losses by varying the ventilation rate as needed to maintain a minimum acceptable IAQ. These systems typically use CO_{2} as a proxy for IAQ and trigger additional ventilation once a threshold is exceeded. Unfortunately, using a gas sensor as a demand switch wastes the opportunity to derive additional information from the continuous analysis of CO_{2} concentration, including the natural infiltration of the space being ventilated and its patterns of occupancy. In order to enable intelligent DCV strategies that can reduce energy waste by supplying heating and cooling only when necessary and eliminate overventilation, the characteristics of the building and its usage patterns are required. The natural infiltration rate is required to determine the current CO_{2} generation rate and to predict future CO_{2} concentrations. Decisions to defer ventilation or plan a heating and cooling schedule require a knowledge of upcoming occupancy patterns. While the occupancy schedule could be entered manually by occupants similar to programming a thermostat, the infiltration rate would require a gas diffusion study. The ideal control system would learn by analysing real-time CO_{2} concentration data and require minimal user setup. Currently available control systems for residential use do none of this, using simple on and off controls and a humidistat. Unfortunately, humidity is a poor proxy for human generated pollutants because it varies slowly as weather patterns change and with the seasons except for short periodic spikes driven by bathroom showers. Furthermore, natural infiltration rates are seldom known and the standard blower door tests provide leakage data at 50 Pa which far exceeds normal conditions.

To confirm the assumption that an HVAC control system can learn building occupancy and calculate the natural infiltration rate from an analysis of CO_{2} concentration data, the author equipped a residential dwelling with a modified ventilation system and instrumentation for the continuous recording of CO_{2}. The data was used to generate a weekly occupancy map and the natural infiltration rate was calculated every time the dwelling became vacant.

This paper presents the results of analysing the CO_{2} concentration curves, the generated vacancy map, and the correlation between natural infiltration and wind speeds and its implication for improved HVAC control systems.

#### 2. Background

Occupant comfort is affected by temperature, humidity, odours, and CO_{2}. The permissible exposure limit for CO_{2} is 5000 ppm [1] which far exceeds levels that will be encountered in a home. Yawning is not caused by CO_{2} even at concentrations of 30,000 ppm [2, p.384]; however, 20% of visitors are dissatisfied with body odours at CO_{2} concentrations of 1880 ppm [3]. CO_{2} generation is a function of the quantity and size of the occupants and their activity levels. CO_{2} also correlates proportionally with pollutants such as body odour and in a more general manner with cooking, cleaning, or smoking that produce volatile organic compounds (VOC), moisture, and particles. This makes CO_{2} a useful proxy for IAQ and for controlling ventilation systems.

Since CO_{2} is nonreactive in the indoor environment and the dwelling volume is a constant, the total mass of CO_{2} and therefore its concentration can change only if more is added by human activity or if it is diluted via air exchange from outside. The portion of CO_{2} concentration change over time within the dwelling due to human activity is the CO_{2} generation rate divided by the total volume of the ventilated area. Assuming that the indoor air is well mixed and therefore the concentration throughout the building is a constant at any instant, the rate of concentration change due to outgoing exhaust air is the mechanical ventilation rate times the indoor concentration divided by the room volume. Similarly since the outdoor concentration is approximately a constant, the change due to incoming fresh air is the mechanical ventilation rate times the outdoor CO_{2} concentration divided by the room volume. This is expressed as a differential equation:
where is indoor concentration of CO_{2} at time (mass/volume), is mechanical ventilation rate (volume/time), is dwelling volume, is ambient atmospheric CO_{2} concentration (mass/volume), and is CO_{2} generation rate by occupants (mass/time).

Separating the variables and integrating from zero to and from to we obtain:
where is indoor concentration of CO_{2} at time (mass/volume) and is initial indoor CO_{2} concentration (mass/volume).

When all occupants vacate a dwelling, the remaining CO_{2} can be used as a tracer gas to measure the natural infiltration rate. When the house is vacant, is zero and the CO_{2} concentration is a decaying exponential asymptotic to . The beginning of the exponential decay marks the onset of vacancy and the natural infiltration rate of the dwelling is its volume divided by the time constant of the decaying exponential.

The CO_{2} generation rate of twenty adults between 21 and 28 years of age, standing and sitting at a desk, was found to average 0.24 L/min and 0.18 L/min, respectively [4]. Assuming that household activities such as standing at a kitchen counter or watching television require similar effort, the CO_{2} generation rate for three adults, (3), would be 79.1 g/hr. Substituting this into (2) with large and a natural infiltration rate of 80 m^{3}/hr yields (4). Therefore, the maximum CO_{2} concentration in a well-mixed home with leakage of 80 m^{3}/hr and three low activity adults would be 1175 ppm:

#### 3. Methodology

Testing was performed on an occupied two-story, 650 m^{3}, wood-framed home, Figure 1, built to Canadian residential standards in 2003. The air change rate with 50 Pa de-pressurization, , was 2.55 air changes per hour (ACH). The maximum flow rate of the HRV was 385 m^{3}/hr with fresh air ducts into bedrooms and living areas and return ducts in bathrooms and kitchen. The kitchen and bathrooms were equipped with dedicated exhaust fans and operable windows.

The HRV unit was modified for computer control of recirculation and fan speed. A Figaro K30 nondispersive infrared (NDIR) CO_{2} sensor was installed in a central sampling device for differential measurements between the indoor and outdoor air streams. Figure 2 shows the HRV unit with the access door open. The translucent box to the left of the HRV is the air source selector which also contains the NDIR sensor.

The CO_{2} concentration was measured at the return duct every 180 seconds while the HRV recirculated air at 153 m^{3}/hr. Every 2.5 hours recirculation was disabled and fresh air was pumped into the building at 385 m^{3}/hr. Once the fresh air readings stabilized the CO_{2} concentration was recorded and the system reverted to recirculation at 153 m^{3}/hr. All measurements were logged into a database.

A moving window of 200 CO_{2} concentration values, corresponding to 8.3 hours, was curve-fit to a theoretical decay by subtracting the background CO_{2} concentration from each datum and then taking its natural logarithm. The linearised data were then fit to a straight line using the least squares method. The window was advanced one data point at a time until the correlation coefficient , calculated as per (5), peaked at a value ≥0.90 marking the beginning of a vacancy period. The window was further advanced sequentially until dropped below 0.50, marking the end of the vacancy period. Each exponential decay meeting the criteria was then mapped into one-hour time slots to produce an occupancy map. The map was then compared with personal notes from the occupants for verification:
where is seconds elapsed from start of window, is natural log of CO_{2} concentration, and is number of measurements in window period.

A Davis Vantage Pro 2 weather station positioned 4 m above ground level and 5 m from the home was used to record local meteorological data.

#### 4. Results

Sequences of exponential CO_{2} decline were identified and the curves they were fit to, shown as red segments in Figure 3, were superimposed on a graph of actual indoor and outdoor CO_{2} concentrations. These overlays were generated using the equation calculated when the correlation coefficient was maximal and used to visually verify the fit.

Table 1 lists the start and end of potential vacancy periods, correlation with a pure exponential decay, decay time constant, lowest CO_{2}, and whether the identification of vacancy was correct.

The vacancy map in Table 2 shows a weekly calendar with 24 hourly time slots. Each period shows the CO_{2} concentration at the end of the period. Cells with multiple values represent multiple vacancies during the same time slot, over the five-week observation period.

#### 5. Discussion

Setting the threshold for correlation to identified all true vacancies and included three false positives.

Using decays where , the 95% confidence interval for the mean of the decay constant during vacancies is (8.5, 11.5) hours with a standard deviation of 2.3 and mean of 10 hours.

Low wind speeds of 4 kmh yielded an infiltration rate of 650 m^{3}/12.2 = 53 m^{3}/hr. The mean time constant of 10 hours represents an infiltration rate of 650 m^{3}/10 = 65 m^{3}/hr.

The infiltration rate was too low to return the indoor CO_{2} concentration to background levels for short vacancy periods so some vacancies would have been missed had the identification been based on absolute levels.

Table 3 shows a variation in infiltration rate with wind speed. The correlation between decay constants and average wind speed during the measurement period was significant with a Pearson correlation coefficient of −0.833 and a value of 0.001. Figure 4 shows a negative correlation between the time constants and average wind speed which is expected. Wind direction is also a factor in asymmetrical houses as is turbulence [5–7]. The low end of the infiltration rate range can be used to conservatively predict CO_{2} and the end of a vacancy period, although if a wind sensor was installed the correlation could provide a better guess.

The maximum mechanical ventilation rate of the HRV was 384 m^{3}/hr; therefore, natural infiltration was approximately 17% of what the HRV was capable of providing. This natural leakage rate prevented interior CO_{2} levels from exceeding 710 ppm during the observation period with one exception during a gathering of many people where it reached 880 ppm. Three adults engaged in light activity would have been expected to ultimately produce CO_{2} concentrations levels of approximately 650 ppm.

Periodic variation in outside CO_{2} concentrations, Figure 3 lower green line, was correlated with changes in indoor concentrations. The intake and exhaust vents from the HRV were installed as per building code but this resulted in mixing of the air streams that varied with atmospheric conditions. To improve the accuracy of outdoor readings, an extension duct was added to the exhaust vent which greatly reduced mixing. Random mixing made differential measurements inaccurate and resulted in missing most vacancy periods because of lower values. Since the NDIR sensor drift was negligible, the analysis was based solely on absolute indoor concentrations.

During the observation period, CO_{2} concentrations never exceeded 710 ppm with three occupants indicating that the greatest saving from DCV would come from disabling fresh air intake most of the time.

A control system which has learned the occupancy map and average infiltration rate over a period of a few weeks could be implemented as follows.

If the CO_{2} rises above a threshold then look ahead to the next hourly time slot. If the probability of vacancy exceeds 50%, calculate the CO_{2} concentration at the end of the projected vacancy using the average infiltration rate. If this is less than the threshold, then do not introduce fresh air.

*Scenario 1. *The CO_{2} concentration is 800 ppm above outdoor levels at 7:30 a.m. Monday morning, and the threshold for ventilation is 600 ppm above background levels. Should mechanical ventilation be used to bring this down by 200 ppm?

The Monday 8 a.m. to 1 p.m. time slot was vacant the last five weeks without exception; the probability of vacancy is estimated at 100%. Anticipating 5 hours of vacancy and time constant of 2.6 hours, CO_{2} should decayed to ppm above outdoor levels. This is negligible and well below the 600 ppm threshold. Therefore, fresh air ventilation at 7:30 is not necessary and will not be required until some time after 1 p.m.

*Scenario 2. *Assume there will be a 1 hour vacancy, and that the probability is 75% that it will occur between 8 a.m to 9 a.m, otherwise conditions are identical to Scenario 1.

Since there is a greater than 50% chance of upcoming vacancy the system expects CO_{2} should decay to ppm above outdoor levels. This is just below the threshold, so ventilation can be deferred.

Had this been a newer home with a tighter building envelop and a time constant of 12 hours, the system would expect CO_{2} to decay to ppm above outdoor levels. In this case, fresh air would be provided until CO_{2} levels returned to a concentration that could be reduced to 600 ppm by natural infiltration during the remainder of the vacancy period. Alternatively, the fan could be run at a lower speed sufficient to augment natural infiltration such that 600 ppm was reached just as occupants were expected to return.

The major energy savings for such a control strategy would come from the following:(1)automatically setting occupied and unoccupied temperatures based on the vacancy probability map;(2)disabling fresh air ventilation when not needed, which in this case, for a threshold of 1000 ppm, would have been all the time;(3)disabling fresh air ventilation before vacancy periods, which for a 9–5 work schedule might eliminate an hour of fresh air each work day depending on the initial concentrations.

#### 6. Conclusion

Periods of vacancy in residential homes were successfully detected using occupant generated CO_{2} as a tracer gas and following the decay in concentration. An accurate occupancy map was generated, indicating the probability and duration of future vacancies. The upper and lower limits of natural infiltration were determined from the time constant of the decay and statistically correlated with wind speed as expected.

The control system for a residential DCV system could be enhanced to use these parameters to optimize energy losses by minimizing excess ventilation while maintaining IAQ. The vacancy map could also be used to automatically program heating and cool systems.

#### Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The author thanks Dr. Tariq Iqbal for research guidance and the Faculty of Engineering and Applied Science, Memorial University of Newfoundland, Technical Services for manufacturing the sampler hardware.

#### References

- “Occupational safety and health standards, toxic and hazardous substances, table z-1 limits for air contaminants 29 cfr 1910. 1000,” 1970.
- R. R. Provine, B. C. Tate, and L. L. Geldmacher, “Yawning: no effect of 3–5% CO
_{2}, 100% O_{2}, and exercise,”*Behavioral and Neural Biology*, vol. 48, no. 3, pp. 382–393, 1987. View at Google Scholar · View at Scopus - “European collaborative action indoor air quality and its impact on man (eca) guidelines for ventilation requirements in buildings,”
*Office for Publications of the European Communities*, no. 11, 1992. - C. Reiff, K. Marlatt, and D. R. Dengel, “Difference in caloric expenditure in sitting versus standing desks,”
*Journal of Physical Activity and Health*, vol. 9, no. 7, pp. 1009–1011, 2012. View at Google Scholar - F. Pan, “Wind-induced internal pressures of buildings with multiple openings,”
*Journal of Engineering Mechanics*, vol. 139, no. 3, pp. 376–385, 2013. View at Publisher · View at Google Scholar - J. H. Oh, G. A. Kopp, and D. R. Inculet, “The UWO contribution to the NIST aerodynamic database for wind loads on low buildings: part 3. Internal pressures,”
*Journal of Wind Engineering and Industrial Aerodynamics*, vol. 95, no. 8, pp. 755–779, 2007. View at Publisher · View at Google Scholar · View at Scopus - G. A. Kopp, J. H. Oh, and D. R. Inculet, “Wind-induced internal pressures in houses,”
*Journal of Structural Engineering*, vol. 134, no. 7, pp. 1129–1138, 2008. View at Publisher · View at Google Scholar · View at Scopus