Code of Federal Regulations (alpha)

CFR /  Title 40  /  Part 1065  /  Sec. 1065.640 Flow meter calibration calculations.

This section describes the calculations for calibrating various flow meters. After you calibrate a flow meter using these calculations, use the calculations described in Sec. 1065.642 to calculate flow during an emission test. Paragraph (a) of this section first describes how to convert reference flow meter outputs for use in the calibration equations, which are presented on a molar basis. The remaining paragraphs describe the calibration calculations that are specific to certain types of flow meters.

(a) Reference meter conversions. The calibration equations in this section use molar flow rate, nref, as a reference quantity. If your reference meter outputs a flow rate in a different quantity, such as standard volume rate, Vstdref, actual volume rate, Vactref, or mass rate, mref, convert your reference meter output to a molar flow rate using the following equations, noting that while values for volume rate, mass rate, pressure, temperature, and molar mass may change during an emission test, you should ensure that they are as constant as practical for each individual set point during a flow meter calibration:[GRAPHIC] [TIFF OMITTED] TR28AP14.017 Where: nref = reference molar flow rate.Vstdref = reference volume flow rate, corrected to a standard

pressure and a standard temperature.Vactref = reference volume flow rate at the actual pressure

and temperature of the flow rate.mref = reference mass flow.pstd = standard pressure.pact = actual pressure of the flow rate.Tstd = standard temperature.Tact = actual temperature of the flow rate.R = molar gas constant.Mmix = molar mass of the flow rate. [GRAPHIC] [TIFF OMITTED] TR28AP14.018 [GRAPHIC] [TIFF OMITTED] TR28AP14.019

(3) Perform a least-squares regression of PDP volume pumped per revolution, Vrev, versus PDP slip correction factor, Ks, by calculating slope, a1, and intercept, a0, as described in Sec. 1065.602.

(4) Repeat the procedure in paragraphs (b)(1) through (3) of this section for every speed that you run your PDP.

(5) The following example illustrates these calculations:

Table 1 of Sec. 1065.640--Example of PDP Calibration Data------------------------------------------------------------------------

a1 (m\3\/ a0 (m\3\/

fnPDP (revolution/s) s) revolution)------------------------------------------------------------------------12.6............................................ 0.841 0.05616.5............................................ 0.831 -0.01320.9............................................ 0.809 0.02823.4............................................ 0.788 -0.061------------------------------------------------------------------------

(6) For each speed at which you operate the PDP, use the corresponding slope, a1, and intercept, a0, to calculate flow rate during emission testing as described in Sec. 1065.642.

(c) Venturi governing equations and permissible assumptions. This section describes the governing equations and permissible assumptions for calibrating a venturi and calculating flow using a venturi. Because a subsonic venturi (SSV) and a critical-flow venturi (CFV) both operate similarly, their governing equations are nearly the same, except for the equation describing their pressure ratio, r (i.e., rSSV versus rCFV). These governing equations assume one-dimensional isentropic inviscid compressible flow of an ideal gas. In paragraph (c)(4) of this section, we describe other assumptions that you may make, depending upon how you conduct your emission tests. If we do not allow you to assume that the measured flow is an ideal gas, the governing equations include a first-order correction for the behavior of a real gas; namely, the compressibility factor, Z. If good engineering judgment dictates using a value other than Z=1, you may either use an appropriate equation of state to determine values of Z as a function of measured pressures and temperatures, or you may develop your own calibration equations based on good engineering judgment. Note that the equation for the flow coefficient, Cf, is based on the ideal gas assumption that the isentropic exponent, g, is equal to the ratio of specific heats, Cp/Cv. If good engineering judgment dictates using a real gas isentropic exponent, you may either use an appropriate equation of state to determine values of g as a function of measured pressures and temperatures, or you may develop your own calibration equations based on good engineering judgment. Calculate molar flow rate, as follows:[GRAPHIC] [TIFF OMITTED] TR28AP14.020 Where: Cd = discharge coefficient, as determined in paragraph (c)(1)

of this section.Cf = flow coefficient, as determined in paragraph (c)(2) of

this section.At = venturi throat cross-sectional area.pin = venturi inlet absolute static pressure.Z = compressibility factor.Mmix = molar mass of gas mixture.R = molar gas constant.Tin = venturi inlet absolute temperature.

(1) Using the data collected in Sec. 1065.340, calculate Cd using the following equation:[GRAPHIC] [TIFF OMITTED] TR28AP14.021 Where: nref = a reference molar flow rate.

(2) Determine Cf using one of the following methods:

(i) For CFV flow meters only, determine CfCFV from the following table based on your values for b and g, using linear interpolation to find intermediate values:

Table 2 of Sec. 1065.640--CfCFV Versus b and g for CFV Flow Meters------------------------------------------------------------------------

CfCFV-------------------------------------------------------------------------

gdexh =

gexh = gair =

b 1.385 1.399

------------------------------------------------------------------------0.000............................................. 0.6822 0.68460.400............................................. 0.6857 0.68810.500............................................. 0.6910 0.69340.550............................................. 0.6953 0.69770.600............................................. 0.7011 0.70360.625............................................. 0.7047 0.70720.650............................................. 0.7089 0.71140.675............................................. 0.7137 0.71630.700............................................. 0.7193 0.72190.720............................................. 0.7245 0.72710.740............................................. 0.7303 0.73290.760............................................. 0.7368 0.73950.770............................................. 0.7404 0.74310.780............................................. 0.7442 0.74700.790............................................. 0.7483 0.75110.800............................................. 0.7527 0.75550.810............................................. 0.7573 0.76020.820............................................. 0.7624 0.76520.830............................................. 0.7677 0.77070.840............................................. 0.7735 0.77650.850............................................. 0.7798 0.7828------------------------------------------------------------------------

(ii) For any CFV or SSV flow meter, you may use the following equation to calculate Cf:[GRAPHIC] [TIFF OMITTED] TR28AP14.022 Where: g = isentropic exponent. For an ideal gas, this is the ratio of specific

heats of the gas mixture, Cp/Cv.r = pressure ratio, as determined in paragraph (c)(3) of this section.b = ratio of venturi throat to inlet diameters.

(3) Calculate r as follows:

(i) For SSV systems only, calculate rSSV using the following equation:[GRAPHIC] [TIFF OMITTED] TR28AP14.023 Where: [Delta]pSSV = Differential static pressure; venturi inlet

minus venturi throat.

(ii) For CFV systems only, calculate rCFV iteratively using the following equation: [GRAPHIC] [TIFF OMITTED] TR28AP14.024

(4) You may make any of the following simplifying assumptions of the governing equations, or you may use good engineering judgment to develop more appropriate values for your testing:

(i) For emission testing over the full ranges of raw exhaust, diluted exhaust and dilution air, you may assume that the gas mixture behaves as an ideal gas: Z = 1.

(ii) For the full range of raw exhaust you may assume a constant ratio of specific heats of g = 1.385.

(iii) For the full range of diluted exhaust and air (e.g., calibration air or dilution air), you may assume a constant ratio of specific heats of g = 1.399.

(iv) For the full range of diluted exhaust and air, you may assume the molar mass of the mixture is a function only of the amount of water in the dilution air or calibration air, xH2O, determined as described in Sec. 1065.645, as follows:[GRAPHIC] [TIFF OMITTED] TR28AP14.144 Example: Mair = 28.96559 g/molxH2O = 0.0169 mol/molMH2O = 18.01528 g/molMmix = 28.96559 [middot] (1 - 0.0169) + 18.01528 [middot]

0.0169Mmix = 28.7805 g/mol

(v) For the full range of diluted exhaust and air, you may assume a constant molar mass of the mixture, Mmix, for all calibration and all testing as long as your assumed molar mass differs no more than 1% from the estimated minimum and maximum molar mass during calibration and testing. You may assume this, using good engineering judgment, if you sufficiently control the amount of water in calibration air and in dilution air or if you remove sufficient water from both calibration air and dilution air. The following table gives examples of permissible ranges of dilution air dewpoint versus calibration air dewpoint: Table 3 of Sec. 1065.640--Examples of Dilution Air and Calibration Air

Dewpoints at Which You May Assume a Constant Mmix------------------------------------------------------------------------

assume the for the following

If calibration Tdew ( [deg]C) is following ranges of Tdew (

. . . constant Mmix [deg]C) during

(g/mol) . . . emission tests \a\------------------------------------------------------------------------dry.............................. 28.96559 dry to 18.0................................ 28.89263 dry to 21.5................................ 28.86148 dry to 22.10............................... 28.81911 dry to 24.15............................... 28.76224 dry to 26.20............................... 28.68685 -8 to 28.25............................... 28.58806 12 to 31.30............................... 28.46005 23 to 34.------------------------------------------------------------------------\a\ Range valid for all calibration and emission testing over the

atmospheric pressure range (80.000 to 103.325) kPa.

(5) The following example illustrates the use of the governing equations to calculate the discharge coefficient, Cd of an SSV flow meter at one reference flow meter value. Note that calculating Cd for a CFV flow meter would be similar, except that Cf would be determined from Table 2 of this section or calculated iteratively using values of b and g as described in paragraph (c)(2) of this section. Example: nref = 57.625 mol/sZ = 1Mmix = 28.7805 g/mol = 0.0287805 kg/molR = 8.314472 J/(mol [middot] K) = 8.314472 (m\2\ [middot] kg)/(s\2\

[middot] mol [middot] K)Tin = 298.15 KAt = 0.01824 m\2\pin = 99.132 kPa = 99132.0 Pa = 99132 kg/(m [middot] s\2\)g = 1.399b = 0.8[Delta]p = 2.312 kPa[GRAPHIC] [TIFF OMITTED] TR28AP14.025

(d) SSV calibration. Perform the following steps to calibrate an SSV flow meter:

(1) Calculate the Reynolds number, Re#, for each reference molar flow rate, nref, using the throat diameter of the venturi, dt. Because the dynamic viscosity, m, is needed to compute Re#, you may use your own fluid viscosity model to determine m for your calibration gas (usually air), using good engineering judgment. Alternatively, you may use the Sutherland three-coefficient viscosity model to approximate m, as shown in the following sample calculation for Re#:[GRAPHIC] [TIFF OMITTED] TR28AP14.026 Where: m0 = Sutherland reference viscosity.T0 = Sutherland reference temperature.S = Sutherland constant.

Table 4 of Sec. 1065.640--Sutherland Three-Coefficient Viscosity Model Parameters----------------------------------------------------------------------------------------------------------------

Temp range within

2% error \b\ limit \b\

kg/(m [middot] K K K kPa

s)----------------------------------------------------------------------------------------------------------------Air.................................. 1.716 [middot] 273 111 170 to 1900 <=1800

10-5CO2.................................. 1.370 [middot] 273 222 190 to 1700 <=3600

10-5H2O.................................. 1.12 [middot] 350 1064 360 to 1500 <=10000

10-5O2................................... 1.919 [middot] 273 139 190 to 2000 <=2500

10-5N2................................... 1.663 [middot] 273 107 100 to 1500 <=1600

10-5----------------------------------------------------------------------------------------------------------------\a\ Use tabulated parameters only for the pure gases, as listed. Do not combine parameters in calculations to

calculate viscosities of gas mixtures.\b\ The model results are valid only for ambient conditions in the specified ranges. Example: m0 = 1.716 [middot] 10-5 kg/(m [middot] s)

[GRAPHIC] [TIFF OMITTED] TR28AP14.027

(3) Perform a least-squares regression analysis to determine the best-fit coefficients for the equation and calculate the equation's regression statistic, SEE, according to Sec. 1065.602.

(4) If the equation meets the criterion of SEE <= 0.5% [middot] Cdmax, you may use the equation to determine Cd for emission tests, as described in Sec. 1065.642.

(5) If the equation does not meet the specified statistical criterion, you may use good engineering judgment to omit calibration data points; however you must use at least seven calibration data points to demonstrate that you meet the statistical criterion. This will usually involve narrowing the range of flow rates for a better curve fit.

(6) Take corrective action if the equation does not meet the specified statistical criterion even after omitting calibration data points. For example, select another mathematical expression for the Cd versus Re# equation, check for leaks, or repeat the calibration process. If you must repeat the process, we recommend applying tighter tolerances to measurements and allowing more time for flows to stabilize.

(7) Once you have an equation that meets the specified statistical criterion, you may use the equation only for the corresponding range of Re#.

(e) CFV calibration. Some CFV flow meters consist of a single venturi and some consist of multiple venturis, where different combinations of venturis are used to meter different flow rates. For CFV flow meters that consist of multiple venturis, either calibrate each venturi independently to determine a separate discharge coefficient, Cd, for each venturi, or calibrate each combination of venturis as one venturi. In the case where you calibrate a combination of venturis, use the sum of the active venturi throat areas as At, the square root of the sum of the squares of the active venturi throat diameters as dt, and the ratio of the venturi throat to inlet diameters as the ratio of the square root of the sum of the active venturi throat diameters (dt) to the diameter of the common entrance to all the venturis. (D). To determine the Cd for a single venturi or a single combination of venturis, perform the following steps:

(1) Use the data collected at each calibration set point to calculate an individual Cd for each point using Eq. 1065.640-4.

(2) Calculate the mean and standard deviation of all the Cd values according to Eqs. 1065.602-1 and 1065.602-2.

(3) If the standard deviation of all the Cd values is less than or equal to 0.3% of the mean Cd, use the mean Cd in Eq. 1065.642-4, and use the CFV only up to the highest r measured during calibration using the following equation:[GRAPHIC] [TIFF OMITTED] TR28AP14.028 Where: [Delta]pCFV = Differential static pressure; venturi inlet

minus venturi outlet.

(4) If the standard deviation of all the Cd values exceeds 0.3% of the mean Cd, omit the Cd value corresponding to the data point collected at the highest r measured during calibration.

(5) If the number of remaining data points is less than seven, take corrective action by checking your calibration data or repeating the calibration process. If you repeat the calibration process, we recommend checking for leaks, applying tighter tolerances to measurements and allowing more time for flows to stabilize.

(6) If the number of remaining Cd values is seven or greater, recalculate the mean and standard deviation of the remaining Cd values.

(7) If the standard deviation of the remaining Cd values is less than or equal to 0.3% of the mean of the remaining Cd, use that mean Cd in Eq. 1065.642-4, and use the CFV values only up to the highest r associated with the remaining Cd.

(8) If the standard deviation of the remaining Cd still exceeds 0.3% of the mean of the remaining Cd values, repeat the steps in paragraph (e)(4) through (8) of this section. [79 FR 23785, Apr. 28, 2014]