Estimation of uncertainty in tracer gas measurement of air change rates.

Simple and economical measurement of air change rates can be achieved with a passive-type tracer gas doser and sampler. However, this is made more complex by the fact many buildings are not a single fully mixed zone. This means many measurements are required to obtain information on ventilation conditions. In this study, we evaluated the uncertainty of tracer gas measurement of air change rate in n completely mixed zones. A single measurement with one tracer gas could be used to simply estimate the air change rate when n = 2. Accurate air change rates could not be obtained for n ≥ 2 due to a lack of information. However, the proposed method can be used to estimate an air change rate with an accuracy of <33%. Using this method, overestimation of air change rate can be avoided. The proposed estimation method will be useful in practical ventilation measurements.

1. Introduction

Ventilation is an essential and effective countermeasure for reducing chemical pollution in indoor environments. Ventilation conditions are commonly evaluated by measuring the ventilation volume [m3/h], which is the uptake volume of outside air per unit time. An additional indicator is the air change rate [1/h], which is the ratio of the ventilation volume to volume of the indoor space. The air change rate is usually measured using a tracer gas. The main tracer gas methods use concentration decay, continuous emission, and steady concentration. In the concentration decay method, the tracer gas is released as a pulse into an indoor space and the decay of its concentration with time is continuously monitored. The air change rate can be estimated from the decay trend. This method is suitable for short term (up to several hours) monitoring of air change rate. In the continuous emission method, a constant amount of tracer gas is continuously emitted into an indoor space. The steady state concentration of the tracer gas is measured and used to calculate the air change rate. This method is generally used for long-term air change rate measurements. In the steady concentration method, the tracer gas concentration in the ventilated space is kept constant by controlling the emission rate of the gas. In this method, the variation in the ventilation rate with time is monitored. For all of these methods, complete mixing of the ventilated space is required to accurately estimate air change rate. For detailed discussions of these techniques, see Lagus and Persily [1], Serman [2], AIVC [3], or ASHRAE [4].

Passive-type tracer gas doser/sampler systems are a simple and cheap method for measurement of air change rate. Operation of the doser and sampler does not require any electric power, and it is simple enough for residents to operate themselves. Perfluorocarbon (PFC) is a common tracer gas used in these systems. Dietz and Cote [5] and Spengler et al. [6] developed ventilation measurements utilizing PFC. The technique has been researched and applied extensively in North America and Northern Europe [7,8]. In the United States, the Environmental Protection Agency (EPA) recommend methods utilizing PFT-CAT (Method IP4-A) or SF6 (Method IP4-B) for standard measurement of air change rates. However, measurements are more complex when many buildings do not consist of a single fully mixed zone. This means many measurements are required to obtain information on ventilation conditions.

The Conjunction Of Multizone Infiltration Specialists (COMIS) model [9] was developed to estimate ventilation conditions for multiple fully mixed zones. Input parameters for COMIS include air tightness and weather conditions. Yoshino et al. [10] reviewed ventilation measurements for multiple zones using various methods, such as the use of smoke to visualize airflow, and the application of four introduced tracer gases. Okuyama [11] used only one tracer gas to estimate air change rate for multiple zones. In this technique, the tracer gas emission rate was varied and changes in the tracer gas concentration were monitored. The data were used to statistically estimate the air change rate of the multiple zones. However, these methods are still complex, and for practical application simpler and more convenient ventilation measurements are required. Miller et al. [12] and Sherman [13] have discussed the estimation method for ventilation conditions for multi-zone building. Riffat et al. [14] reported the accuracy of single-tracer gas measurement for estimation of air flows between two zones.

In this study, we evaluated the measurement of air change rate for two fully mixed zones from a single measurement with one tracer gas. We discuss the uncertainty in this measurement for n fully mixed zones, and propose a simple method for estimation of air change rate for two fully mixed zones. This provides a simple method for practical application in ventilation studies.

2. Theory

2.1. Uncertainty in the Measurement of Air Change Rate with a Tracer Gas

We first considered air change rate measurements in buildings of n fully mixed zones, where n is a natural number. The ventilation situations in buildings with 1–3 zones are illustrated in Figure 1. Measurements aim to determine the volume of air uptake from the outside into the n zones (air change rate). The uptake volume of air, F, can be represented by the following equation:

Figure 1

Schematic of buildings with 1, 2, or 3 fully mixed zones.

where F0 represents the air uptake volume to each zone. A single tracer gas is emitted at a known rate, E, in zone i. The emitted tracer gas is collected by the passive-type air sampler and the average tracer gas concentration in zone i, C, is measured. The volume of each zone, V, is considered to be a known parameter as it can be measured. On the other hand, the volume of airflow from zone i to zone j, F, cannot be measured. All F have non-negative values. Generally, a gas that does not exist naturally outdoors is selected as the tracer gas, which means its concentration outside (zone 0), C0, is zero. Usually, the measurement time for the tracer gas method is sufficient that the tracer gas concentration is considered to be constant. In each zone, the air volume is also considered to be constant:

This equation reflects the conservation of mass (air) in each zone.

If all fully mixed zones are in contact with each other and the air can freely pass between any two zones, then the number of unknown parameters F is n(n + 1). It should be noted that this provides the maximum number, and the actual number of unknown parameters is reduced compared with this when F = 0 before measurement. If zone i and j is not physically connected, F and F can be assumed as 0. A total of 2n balanced equations can be established, which includes tracer mass balance equations and air mass balance equations for the n fully mixed zones. Therefore, when calculating air change rate the information deficit is equal to n(n − 1) (= n(n + 1) − 2n). This is only equal to 0 when n = 1, and in all other cases we cannot obtain an accurate air change rate. Table 1 shows the maximum information deficit for n fully mixed zones.

Table 1

The maximum information deficit for n fully-mixed zones.

Zone number	Maximum number of unknown parameters (I)	Obtainable information (II)	Maximum information deficit (III) = (I) – (II)
1	2	2	0
2	6	4	2
3	12	6	6
4	20	8	12
n	n(n + 1)	2n	n(n − 1)

To avoid this deficit, we have to supplement n(n − 1) data. For example, preparing and measuring the steady state concentrations of n types of tracer gases for n zones allows calculation of an accurate ventilation rate. Otherwise, with one tracer gas, n ventilation measurement trials are required to calculate the steady state concentration of the tracer gas in n zones. In this case, for each ventilation measurement, the location of the tracer gas doser is varied in n ways. However, these methods require either many tracer gases or many measurements, which means they cannot be practically applied. In addition, the time and money available for air change rate measurements are limited. Consequently, we evaluated the estimation of air change rate from a single measurement using one tracer gas.

2.2. Calculation When n = 1

As already mentioned, F for a single fully mixed zone [Figure 1(a)] can be calculated accurately even from only a single measurement with one tracer gas. The tracer mass balance is given by:

The indoor air mass balance is given by:

where E1 is a known value and C1 is obtained by a tracer gas concentration measurement. F can then be accurately calculated as follows:

2.3. Calculation When n = 2

If building consists of two fully mixed zones [Figure 1(b)], it is more difficult to measure the air uptake F than for n = 1. The tracer mass balance for the two fully mixed zones is as follows:

where C and C are considered to be constant after stabilization of tracer gas concentration. Therefore, we can transform the above equations to:

The air balance is represented by:

There are now six unknown parameters, F01, F02, F10, F12, F20, and F21 in the above four balanced equations, which means F, which is defined by the following equation, cannot be calculated accurately:

A situation such as this usually occurs with n ≥ 2, and the maximum information deficit is equal to n(n − 1). However, practically it is difficult to increase the number of tracer gases used or the number of measurements, which both consequently increase the measurement time. Therefore, it is necessary to estimate the ventilation volume for these zones.

3. Discussion

3.1. Estimation Method for Air Change Rate of Two Fully Mixed Zones from a Single Measurement with One Tracer Gas

For air change rate estimation of two fully mixed zones from limited information, we can utilize information about the volume of each zone, and tracer gas emission rates in each zone. Information about tracer gas concentrations in each zone can be obtained from passive air sampler measurements. As discussed in Section 2.3., four equations can be established for tracer gas mass balance and air balance for the two zones. There are six unknown parameters, F01, F02, F10, F12, F20, and F21, and F cannot be accurately estimated. However, even in this case, we can estimate the air change rate with a practical level of accuracy using the following estimation method. This is due to the symmetry of Equations (8)–(12) and non-negative restrictions for the six unknown parameters.

Here, we use estimated air intake, F, as an indicator of F. F is given by the sum of emission rates divided by the volume-weighted average tracer gas concentration, C̄:

where C̄ is given by:

One of the purposes of the study is to show that F is suitably accurate for practical applications and has a small range of error. The error between F and F can be indicated by the normalized parameter y:

With non-negative restrictions on both parameters, the range of y is at least:

Because F and F are not independent, the range of y is even more limited. A further parameter, x, can be introduced as an indicator of concentration differences between the two zones:

With non-negative restrictions on C and C, the range for x is:

If the steady state concentration of tracer gas in each zone is equal, then the two zones can be treated as one fully mixed zone. So, when x = 0, y is 0, and the ventilation rate can then be accurately calculated. The range for y is mathematically more limited.

The restriction condition is non-negative restriction of six airflow volumes. There are six unknown parameters in the four balanced equations (Equations (8)–(11)), and we can express y as a function of two arbitrary unknown parameters. If F10 and F20 are selected as these two parameters, then y can be expressed as:

If we select F01 and F12, then y can be expressed as:

and for F02 and F21:

It should be noted that only the symmetry of parameters was considered when selecting these three pairs, and other expressions would also give the same results.

In Equation (19), from non-negative restriction of F10 and F20, the range of y as a function of x can be calculated as follows:

In Equation (20), from non-negative restrictions of F01 and F12, we can obtain the range of y as a function of x as follows:

In Equation (21), from non-negative restrictions of F02 and F21, we can obtain the range of y as a function of x as follows:

If non-negative restrictions for all six unknown parameters are considered, then combination of Equations (22)–(27) can accurately express the ranges for y (Appendices I–III). If f1(x) − f4 (x) are determined as follows:

the ranges for y can be expressed as follows:

3.2. Application of the Estimation Method

The estimation method can then be illustrated in application to a specific situation. It should be noted that the kind of tracer gas does not affect the results. This method is based on steady state. Thus, that requires constant airflow rates over a sufficient period of time for the concentrations to stabilize. Here, we attempt to simply measure the air change rate of two fully mixed zones. V1 and V2 are assumed to be 10 m3 and 30 m3, respectively. E1 and E2 are set as 100 μg/h and 200 μg/h, respectively. Figure 2 illustrates this situation. In the case, f1(x) − f4(x) were calculated as:

Figure 2

Schematic of a building consists of two fully mixed zones.

Thus, the range of y is limited as shown in Figure 3. If C1 and C2 are measured to be 30 μg/m3 and 10 μg/m3, respectively, air intake F can be estimated at 20 m3/h using Equations (13) and (14). In this case, x is calculated to be 0.5, so, f1 (0.5) = −0.33, f2 (0.5) = 1, and f3 (0.5) = 0. Thus, from Equation (32):

Figure 3

The range of y (colored area) when V1, V2, E1, and E2 are equal to 10 m3, 30 m3, 100 μg/h, and 200 μg/h, respectively.

So, the estimated F (20 m3/h) is accurate within 33%, and is an underestimate. From Equation 15:

It should be noted that F can be appropriately normalized in accordance with the aim of measurement. If we want to minimize the error rate, we can choose 25 ± 5 m3/h as Fe (error range is ±20%). To minimize the error rate, E1 and E2 should be set at the same value. In this case, the maximum error range of air intake Fe estimated by the methods is always <33%. This error value is equivalent to that reported by Riffat with single tracer gas measurement [14]. Miller et al. [12] reported that airflow rates in two-zone building could be estimated with an accuracy of 8% by two tracer gas decay experiments with the nonlinear least-squares minimization method in controlled conditions. One feature of this method is its simplicity in calculation operations. As described above, the proposed method can estimate air uptake volume with a practical accuracy only by using four arithmetic operations. If this estimation method can be combined with real time monitoring of a tracer gas concentration, we can extemporarily obtain the estimated air uptake without complex calculations. Furthermore, by using this method, overestimation of air change rate can be easily avoided. When aiming to avoiding overestimation of air change rate, the estimation method can give a minimum air change rate (in this case, F,min. = 20 m3/h). Figure 4 summarizes the estimation flow chart for air intake and its error range.

Figure 4

Flow chart for the estimation of air intake and its error range.

4. Conclusions

In this study, we have evaluated the uncertainty in tracer gas measurement of air change rate for n fully mixed zones. We proposed a simple method to estimate the air change rate for two fully mixed zones from a single measurement using one tracer gas. Accurate air change rate could not be obtained for n ≥ 2 due to lack of information. However, the proposed method can provide an estimated air change rate with an accuracy of <33%. Using this method, overestimation of air change rate can be avoided. The proposed estimation method could find practical applications in ventilation studies.