2.5.5.3. Propagation of error for many variables

2. Measurement Process Characterization
2.5. Uncertainty analysis
2.5.5. Propagation of error considerations

2.5.5.3. Propagation of error for many variables

Example from fluid flow with a nonlinear function

Computing uncertainty for measurands based on more complicated functions can be done using basic propagation of errors principles. For example, suppose we want to compute the uncertainty of the discharge coefficient for fluid flow (Whetstone et al.). The measurement equation is

where

Assuming the variables in the equation are uncorrelated, the squared uncertainty of the discharge coefficient is

(s_Cd)^2 = (dCd/dm)^2*(s_m)^2
+(dCd/dd)^2*(s_d)^2 + (dCd/dD)^2*(s_D)^2 + (dCd/drho)^2*(s_rho)^2
+ (dCd/dDeltaP)^2*(s_DeltaP)^2

and the partial derivatives are the following.

dCd/dm = (sqrt(1-(d/D)^4))/
(K d^2 F sqrt(rho) sqrt(DeltaP))

dCd/dd =
(-2 m d)/(K F D^4 sqrt(rho) sqrt(DeltaP) sqrt(1-(d/D)^4)))
-2 m sqrt(1-(d/D)^4)/(K d^3 F sqrt(rho) sqrt(DeltaP))

dCd/dC =
2 m d^2 / (K F D^5 sqrt(rho) sqrt(DeltaP) sqrt(1-(d/D)^4))

dCd/drho =
-m sqrt(1-(d/D)^4) / (2 K F d^2 rho^(3/2) sqrt(DeltaP))

dCd/dDeltaP =
-m sqrt(1-(d/D)^4) / (2 K F d^2 sqrt(rho) DeltaP^(3/2))

Software can simplify propagation of error

Propagation of error for more complicated functions can be done reliably with software capable of symbolic computations or algebraic representations.

Symbolic computation software can also be used to combine the partial derivatives with the appropriate standard deviations, and then the standard deviation for the discharge coefficient can be evaluated and plotted for specific values of the secondary variables, as shown in the comparison of check standard analysis and propagation of error.

Simplification for dealing with multiplicative variables

Propagation of error for several variables can be simplified considerably for the special case where:

the function, Y, is a simple multiplicative function of secondary variables, and
uncertainty is evaluated as a percentage.

For three variables, X, Z, W, the function

has a standard deviation in absolute units of

s(Y) = SQRT{(Z*W)**2*s(x)**2 + (X*W)**2*s(z)**2 + (X*Z)**2*s(w)**2} =
Y*SQRT{s(x)**2/X**2 + s(z)**2/Z**2 + s(w)**2/W**2}

In percent units, the standard deviation can be written as

s(Y)/Y = SQRT{s(x)**2/X**2 + s(z)**2/Z**2 + s(w)**2/W**2}

if all covariances are negligible. These formulas are easily extended to more than three variables.