Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Example: Volume Measurement Uncertainty Budget

This example demonstrates two distinct approaches for measuring the volume of a cylinder, rigorously applying the principles of uncertainty analysis established by the Guide to the Expression of Uncertainty in Measurement (GUM). This process requires formulating the measurement model, evaluating all input standard uncertainties (Type A and Type B), propagating those uncertainties using sensitivity analysis, and reporting a final result with expanded uncertainty.

Example: Volume Measurement Uncertainty Budget

The core measurement result format must be: Result = best estimate ± uncertainty (confidence level).

Method A: Direct Geometric Measurement

This method calculates volume (VV) based on measurements of the cylinder’s diameter (DD) and length (LL) using high-precision dimensional instruments like a micrometer and caliper. The measurement model is an explicit measurement equation Y=f(Xi)Y=f(X_i).

Step 1: Formulation and Calculation of Best Estimate

The measurement model for the volume of a cylinder is:

V=f(D,L)=π4D2LV = f(D, L) = \frac{\pi}{4} D^2 L

Assumed Values:

Best Estimate:

V=π4(50.00 mm)2(150.00 mm)294,524.3 mm3V = \frac{\pi}{4} (50.00 \text{ mm})^2 (150.00 \text{ mm}) \approx 294,524.3 \text{ mm}^3

Step 2: Evaluating Standard Uncertainties (u(xi)u(x_i))

Input uncertainties are evaluated either statistically (Type A) or non-statistically (Type B). The standard deviation of the mean (sxˉ=sx/Ns_{\bar{x}} = s_x / \sqrt{N}) is used for Type A components. Type B components often derive from rectangular distributions (e.g., resolution limits).

SourceTypeInput Quantity (XiX_i)Standard Uncertainty (u(xi)u(x_i))Distribution
Length RepeatabilityAu(L)Au(L)_A0.015 mm0.015 \text{ mm}Normal
Caliper ResolutionBu(L)B1u(L)_{B1}0.0029 mm0.0029 \text{ mm}Rectangular
Micrometer ResolutionBu(D)B1u(D)_{B1}0.00029 mm0.00029 \text{ mm}Rectangular
Calibration Bias (k=2k=2)Bu(L)B2,u(D)B2u(L)_{B2}, u(D)_{B2}0.005 mm0.005 \text{ mm}Normal

Combined Standard Uncertainty for Inputs (RSS):

u(L)=(0.015)2+(0.0029)2+(0.005)20.016 mmu(L) = \sqrt{(0.015)^2 + (0.0029)^2 + (0.005)^2} \approx 0.016 \text{ mm}
u(D)=(0.00029)2+(0.005)20.005 mmu(D) = \sqrt{(0.00029)^2 + (0.005)^2} \approx 0.005 \text{ mm}

Step 3: Sensitivity Coefficients (CiC_i)

The sensitivity coefficients are the partial derivatives of the model:

CL=VL=π4D21963.5 mm2C_L = \frac{\partial V}{\partial L} = \frac{\pi}{4} D^2 \approx 1963.5 \text{ mm}^2
CD=VD=π2DL11,781.0 mm2C_D = \frac{\partial V}{\partial D} = \frac{\pi}{2} D L \approx 11,781.0 \text{ mm}^2

Step 4: Combining Standard Uncertainties (ucu_c)

Assuming the input errors (DD and LL) are independent (uncorrelated), the combined standard uncertainty is calculated using the Root-Sum-Squares (RSS) method (Law of Propagation of Uncertainty, LPU):

uc(V)=[CLu(L)]2+[CDu(D)]2u_c(V) = \sqrt{\left[C_L \cdot u(L)\right]^2 + \left[C_D \cdot u(D)\right]^2}
uc(V)=(1963.50.016)2+(11,781.00.005)2(31.4)2+(58.9)266.8 mm3u_c(V) = \sqrt{\left(1963.5 \cdot 0.016\right)^2 + \left(11,781.0 \cdot 0.005\right)^2} \approx \sqrt{(31.4)^2 + (58.9)^2} \approx 66.8 \text{ mm}^3
Input VariableUncertainty Contribution (Ciu(xi)C_i \cdot u(x_i)) (mm3\text{mm}^3)Contribution Squared (mm6\text{mm}^6)
Length (LL)31.4985.96
Diameter (DD)58.93469.21
CombinedN/A4455.17\mathbf{4455.17}
uc(V)u_c(V)66.8 mm366.8 \text{ mm}^3N/A

The dominant contributor to the volume uncertainty in this geometric method is the uncertainty in the Diameter (DD), as evidenced by the significantly larger squared contribution term.

Step 5: Expanded Uncertainty and Reporting

The expanded uncertainty (UU) is calculated using a coverage factor (kk). Using k=2k=2 provides a confidence level of approximately 95%.

U=kuc(V)=266.8 mm3134 mm3U = k \cdot u_c(V) = 2 \cdot 66.8 \text{ mm}^3 \approx 134 \text{ mm}^3

Reported Result (Method A):

V=(294,524±134) mm3(k=2,95% confidence)V = (294,524 \pm 134) \text{ mm}^3 \quad (k=2, 95\% \text{ confidence})

Method B: Indirect Gravimetric Measurement (Archimedes)

This method relates the volume (VV) to the measured mass (mm) and known density (ρ\rho) of the cylinder material. The fidelity of the final uncertainty statement depends entirely on establishing a complete measurement equation and meticulously evaluating the uncertainty of all input quantities.

Step 1: Formulation and Calculation of Best Estimate

The measurement model for volume is:

V=f(m,ρ)=mρV = f(m, \rho) = \frac{m}{\rho}

Assumed Values:

Best Estimate:

V=1570.005.33294.55 cm3V = \frac{1570.00}{5.33} \approx 294.55 \text{ cm}^3

Step 2: Evaluating Standard Uncertainties (u(xi)u(x_i))

SourceTypeInput Quantity (XiX_i)Standard Uncertainty (u(xi)u(x_i))Distribution
Mass RepeatabilityAu(m)Au(m)_A0.020 g0.020 \text{ g}Normal
Balance CalibrationBu(m)B1u(m)_{B1}0.025 g0.025 \text{ g}Normal
Density ReferenceBu(ρ)B1u(\rho)_{B1}0.0058 g/cm30.0058 \text{ g/cm}^3Rectangular
Thermal ExpansionBu(ρ)B2u(\rho)_{B2}0.0029 g/cm30.0029 \text{ g/cm}^3Rectangular

Combined Standard Uncertainty for Inputs (RSS):

u(m)=(0.020)2+(0.025)20.032 gu(m) = \sqrt{(0.020)^2 + (0.025)^2} \approx 0.032 \text{ g}

u(ρ)=(0.0058)2+(0.0029)20.0065 g/cm3u(\rho) = \sqrt{(0.0058)^2 + (0.0029)^2} \approx 0.0065 \text{ g/cm}^3

Step 3: Sensitivity Coefficients (CiC_i)

Cm=Vm=1ρ0.1876 cm3/gC_m = \frac{\partial V}{\partial m} = \frac{1}{\rho} \approx 0.1876 \text{ cm}^3/\text{g}
Cρ=Vρ=mρ255.26 cm6/g2C_\rho = \frac{\partial V}{\partial \rho} = -\frac{m}{\rho^2} \approx -55.26 \text{ cm}^6/\text{g}^2

Step 4: Combining Standard Uncertainties (ucu_c)

uc(V)=[Cmu(m)]2+[Cρu(ρ)]2u_c(V) = \sqrt{[C_m \cdot u(m)]^2 + [C_\rho \cdot u(\rho)]^2}
uc(V)=(0.18760.032)2+(55.260.0065)2(0.0060)2+(0.359)20.359 cm3u_c(V) = \sqrt{(0.1876 \cdot 0.032)^2 + (-55.26 \cdot 0.0065)^2} \approx \sqrt{(0.0060)^2 + (-0.359)^2} \approx 0.359 \text{ cm}^3
Input VariableUncertainty Contribution (Ciu(xi)C_i \cdot u(x_i)) (cm3\text{cm}^3)Contribution Squared (cm6\text{cm}^6)
Mass (mm)0.00603.6×1053.6 \times 10^{-5}
Density (ρ\rho)-0.3590.12888
CombinedN/A0.128916\mathbf{0.128916}
uc(V)u_c(V)0.359 cm30.359 \text{ cm}^3N/A

Step 5: Expanded Uncertainty and Reporting

Using k=2k=2 for approximately 95% confidence:

U=kuc(V)=20.359 cm30.72 cm3U = k \cdot u_c(V) = 2 \cdot 0.359 \text{ cm}^3 \approx 0.72 \text{ cm}^3

Reported Result (Method B):

V=(294.55±0.72) cm3(k=2,95% confidence)V = (294.55 \pm 0.72) \text{ cm}^3 \quad (k=2, 95\% \text{ confidence})

Method Comparison and Interpretation

By quantifying the specific contributions of each input quantity via sensitivity coefficients, the uncertainty budget functions as the definitive engineering tool for optimization.

  1. Precision: Converting Method B’s uncertainty to mm3\text{mm}^3 (0.72 cm3=720 mm30.72 \text{ cm}^3 = 720 \text{ mm}^3), we see that Method A (U134 mm3U \approx 134 \text{ mm}^3) is significantly more precise than Method B (U720 mm3U \approx 720 \text{ mm}^3) under these assumptions.

  2. Dominant Errors: In Method B, the density reference (u(ρ)u(\rho)) contributes overwhelmingly to the total uncertainty (0.12888 cm60.12888 \text{ cm}^6) compared to the mass measurement (3.6×105 cm63.6 \times 10^{-5} \text{ cm}^6). To improve Method B, efforts must be focused on increasing the precision of the density value, not the precision of the mass balance.

Would you like to explore another metrological example, such as evaluating uncertainty for a thermocouple calibration system, or dive deeper into the use of tools like the NIST Uncertainty Machine (NUM) for Monte Carlo simulations?