uncertainty of measurement: water vapour ...

8 downloads 0 Views 82KB Size Report
Apr 25, 2016 - The rate at which water vapour is transmitted through a sheet material is important, for example, if the sheet is to be used to damp-proof a ...
Uncertainty of measurement: water vapour transm...

http://metrologyrules.com/uncertainty-of-measur...

Metrology Rules

UNCERTAINTY OF MEASUREMENT: WATER VAPOUR TRANSMISSION OF MATERIALS APRIL 1, 2016 | HANS LIEDBERG | LEAVE A COMMENT

The rate at which water vapour is transmitted through a sheet material is important, for example, if the sheet is to be used to damp-proof a building. The documentary standards ASTM E96 and SANS 952-1 describe methods of evaluating water vapour transmission (WVT), using a desiccant or water to establish a humidity gradient across the sheet. The following article will brie�y describe the method, and then focus on estimation of the uncertainty of the calculated WVT value. The Desiccant Method [E96 clause 4.1] involves sealing the specimen to be tested over the mouth of an impermeable test dish, with desiccant (anhydrous calcium chloride) inside. This dish is placed inside a test chamber or enclosure, where the temperature and humidity are controlled at (38 ± 1) °C and (95 ± 5) %rh [952-1 clause 6.11.1.5], for example. The dish is brie�y removed from the chamber and weighed, periodically, over a period of several days. The mass gain indicates the amount of water vapour moving through the specimen into the desiccant. Typically, three specimens are tested [E96 clause 9.1], with a fourth “dummy” (control) specimen mounted on a dish and weighed, but without any desiccant or water in the dish to maintain a humidity gradient [E96 clause 9.6]. The mass change of the dummy is subtracted from the change in each of the test specimens, to cancel out the effect of variations in temperature and air buoyancy [E96 clause 11.3]. (In practice, vigorous air circulation within the test chamber, required to keep the surrounding air uniform in temperature and humidity [E96 clause 6.2], causes the dummy to gain signi�cant mass, too, because it was sealed at a lower ambient humidity than the 95 %rh maintained in the chamber.) Measurements of specimen mass are performed (perhaps once every 24 hours), until approximately six consecutive changes in mass (corrected for variation in the dummy), plotted versus time, fall “on a reasonably straight line” [952-1 clause 6.11.4.2.3]. The slope of this line (equal to rate of mass change, in grams per hour) is converted to WVT in grams per square metre per 24 hours [952-1 clause 6.11.5], using the measured area of the mouth of each dish. The mean of the WVT values of the three specimens is reported to the customer. Here is a data set, for use in our discussion below: Time

Dummy

(h)

m (g)

dm (g)

m (g)

dm (g)

dm_corr (g)

m (g)

dm (g)

dm_corr (g)

m (g)

dm (g)

dm_corr (g)

0

163.4055

0.0000

197.6220

0.0000

0.0000

200.0488

0.0000

0.0000

190.0098

0.0000

0.0000

24

163.4061

0.0006

197.6260

0.0040

0.0034

200.0545

0.0057

0.0051

190.0167

0.0069

0.0063

48

163.4112

0.0057

197.6346

0.0126

0.0069

200.0615

0.0127

0.0070

190.0238

0.0140

0.0083

72

163.4148

0.0093

197.6408

0.0188

0.0095

200.0685

0.0197

0.0104

190.0326

0.0228

0.0135

120

163.4165

0.0110

197.6462

0.0242

0.0132

200.0753

0.0265

0.0155

190.0392

0.0294

0.0184

144

163.4182

0.0127

197.6554

0.0334

0.0207

200.0847

0.0359

0.0232

190.0471

0.0373

0.0246

Fitted slope (g/h):

1 of 3

Specimen 1

Specimen 2

0.000130

Specimen 3

0.000145

0.000158

25/04/2016 12:14

Uncertainty of measurement: water vapour transm...

http://metrologyrules.com/uncertainty-of-measur...

Area (mm^2):

3117

3117

3117

WVT (g/m^2 /24h):

1.00

1.12

1.21

The slopes above were �tted using the LINEST spreadsheet function. Using the formula =LINEST(measured_y_values, measured_x_values, intercept, stats), with intercept=TRUE (allowing a non-zero y-intercept) and stats=TRUE (displaying not just the �tted values of slope and intercept, but also four additional rows of statistical information), for specimen 1, the following output is obtained: 0.000 130

0.000 11

0.000 013

0.001 09

0.963

0.001 59

103.10

4

0.000 26

0.000 01

The �rst row gives �tted slope and y-intercept, the second row standard errors (uncertainties at coverage factor k=1) of the �tted slope and intercept, the third row correlation coef�cient R^2 and standard error in the y-values, the fourth row the Fisher F-statistic and the degrees of freedom (number of data points minus number of �tted parameters). We will discuss the use of some of these statistical measures below. (The �fth row is not of interest to this discussion.) How may we estimate the uncertainty of measurement (UoM) of this WVT value? First, consider that the basic task is to determine the change in mass over the change in time, or

. The absolute

accuracy of the balance used to weigh the dishes is not important, nor are the absolute times indicated by the watch: only the changes are of interest. How can we estimate the accuracy of ? We could evaluate the sensitivity of the balance [OIML R 76-1 clause T.4.1] to small changes in a load around 200 g (in the above example). However, there is an easier way, which also takes into account any effects that vary randomly from one balance reading to the next. As we have six data points to determine only two unknowns (slope and y-intercept of the straight line), we have what is called an “over-determined system”, for which we �nd the “best �t” line by the method of least squares. This line typically does not pass through any of the data points, but is “as close as possible” to all (in other words, the best compromise to all available data). The differences between the measured data points and the �tted line give a quantitative estimate of random errors in the measured values. These differences are conveniently combined in the standard error of the �tted slope, namely, 0.000 013 g/h for specimen 1 above. How do we expand this standard (k=1) uncertainty, to reach a level of con�dence of approximately 95%? For a very large data set, we would simply multiply u(k=1) by the coverage factor k=2. However, as our data set is small, we should enlarge the coverage factor somewhat, to compensate for our limited knowledge. The t-distribution of Student tells us what this k-value should be: the spreadsheet function =TINV(1-0.95, 4) gives the value of Student’s t-distribution for 4 degrees of freedom and 95% level of con�dence, namely, 2.8. Multiplying u(k=1) = 0.000 013 g/h by 2.8, we obtain the expanded uncertainty 0.000 036 g/h. (This is what U(k=2) would have been, if we had a very large data set, or, in other words, almost in�nite degreees of freedom.) So, the slope for specimen 1 is (0.000 130 ± 0.000 036) g/h. Or, as a relative uncertainty,

, or 27%. Note that this is signi�cantly worse than the

balance sensitivity of 1% of total mass gain (which would translate to 0.000 001 g/h) required in the test methods [E96 clause 6.3, 952-1 clause 6.11.4.2.2], indicating that balance sensitivity is, in this case, negligible compared to other factors causing “noise” in the mass readings. It highlights the importance of having more data points than unknowns and performing a least squares �t, as this uncertainty component would otherwise be grossly underestimated. It also shows how the standard error in the slope is far more useful than the correlation coef�cient R^2, in quantifying uncertainties. Now, how do we estimate the uncertainty in

2 of 3

? (In the previous paragraph, it seems we already obtained a

25/04/2016 12:14

Uncertainty of measurement: water vapour transm...

http://metrologyrules.com/uncertainty-of-measur...

“complete” uncertainty for the slope. However, the �tting procedure we used assumes no uncertainty in the x-values, only uncertainty in the y-values, so we had better look at the uncertainty in the time intervals, too.) The documentary standards require time intervals to be measured to an accuracy of 1% (for example, a 24 hour interval to an accuracy of 15 minutes), so,

[E96 clause 11.3, 952-1 clause

6.11.4.2.2]. We can see that the relative uncertainty in

is far smaller than that in the �tted slope, so we

expect the uncertainty in slope to completely dominate the �nal, combined, uncertainty. The formula for WVT is

. When parameters are combined by simple multiplication or

division, the relative uncertainties may be added simply in quadrature. So, . As

is typically smaller than 0.01, it, like

can be neglected, so that the �nal relative uncertainty is

,

, or 27%. In other words, U(k=2) =

1.00 g/m^2/24h * 0.27 = 0.27 g/m^2/24h. The above is the UoM for specimen 1. Specimens 2 and 3 have similar uncertainties (0.29 g/m^2/24h and 0.25 g/m^2/24h, respectively). Now, what is the uncertainty of the mean of the three specimens? If the three uncertainties are similar in magnitude, the uncertainty of the mean is . If the uncertainties are quite different, the weighted mean (giving more weight to those specimens with smaller uncertainties) could be used: . The uncertainty of the weighted mean is

, which also works out to be

0.16 g/m^2/24h, in the above example.

3 of 3

25/04/2016 12:14