CS-041: Flow Measurement Orifice Plate Sizing and Uncertainty analysis

CHALLENGE

The main challenge for this case study is the application of Flownex to:

• Size an orifice plate to be used in conjunction with a flow transmitter to serve as an accurate flow meter for natural gas.

• Evaluate the measurement uncertainty of the flow meter with variations in operating conditions and manufacturing tolerances.

BENEFITS

Flownex is an ideal tool to design gas flow systems, including piping, valving and most other components that are typically found in the oil and gas industry. Not only is Flownex also the ideal tool to design accurate flow measurement orifice plates, but it also has the capability to evaluate the uncertainty of the flow meter in general and the orifice plate in particular with the inevitable variations in conditions and manufacturing tolerances.

SOLUTION

Using the Designer and the Sensitivity Analysis features built into Flownex, the orifice plate can be designed and its operational uncertainty evaluated when functioning in combination with the pressure transducer/flow transmitter.

INTRODUCTION

Flow transmitters use a variety of measurement techniques, the simplest (and cheapest) of which is probably the orifice plate. These orifice plates are often designed in accordance with ISO 5167 or BS 1042. The pressure transducers used by flow transmitters are typically calibrated by mapping the orifice plate minimum to maximum measured pressure difference, to the typical 4-20 mA signals they transmit. At zero flow the transmitter would issue a 4 mA signal. As a rule of thumb, the minimum actual flow would be calibrated to approximately 30% above the minimum signal (i.e. 8.8 mA) whereas the maximum flow would be calibrated to approximately 70% to 80% of full scale (14 to 16 mA).

However, when such an orifice plate is sized, the sizing calculation has to be based on one specific operational point - one pressure, one temperature and one flow rate - resulting in one fixed orifice plate size. In practice, however, the pressure and temperature will not remain constant and neither will the flow rate. Furthermore, the physical orifice plate has to be manufactured and the manufacturing tolerance will also have an influence on the measurement uncertainty of the orifice plate. According to ISO 5167 there is also an uncertainty associated with the calculated discharge coefficient which is directly proportional to the discharge mass flow rate calculated for the orifice plate. This uncertainty is also discussed later in this case study.

Flownex has a very powerful design facility in terms of the Designer feature. The Designer is typically used to “goal seek” for a solution, and in this case study the Designer has been used to find the required orifice plate size to comply with functional requirements. The measurement uncertainty of the orifice plate is evaluated by using the Sensitivity Analysis feature of Flownex. This feature allows for parametric studies where a single independent variable is varied to study its influence on a single or multiple dependent variables. Flownex also offers the ability to perform Monte Carlo analyses where multiple independent variables are randomly varied such that the maximum and minimum combined influence of the independent variables may be observed on the dependent variables.

FLOWNEX MODEL

Figure 1 shows the example network with an orifice plate between two pipes. The link between the upstream pipe and the orifice plate element is used to transfer the pipe internal diameter to the orifice plate element.

For the case at hand, the design information and requirements are as follows:

• Maximum fuel gas flow rate = 1500 Nm3/hr.

• Nominal fuel gas flow rate = 1200 Nm3/hr.

• Minimum fuel gas flow rate = 350 Nm3/hr.

• Maximum temperature = 30°C.

• Nominal temperature = 25°C.

• Minimum temperature = 5°C.

• Supply pressure = 120 kPa-g +3.7 kPa -6.3 kPa.

• Target pressure drop at maximum flow = 5 kPa.

• Target pressure drop at minimum flow = 0.5 kPa.

• Transmitter accuracy = ±0.075% of span (from manufacturer).

• Transmitter repeatability/stability = ±0.01% of span  (from manufacturer).

• Orifice plate diameter manufacturing tolerance = ±0.1 mm.

The nominal operating upstream pressure and temperature are specified at the inlet boundary, whilst the required flow rate is specified as a sink (negative mass source) at the outlet boundary.

SIZING THE ORIFICE PLATE

Flownex uses the BS 1042 formulation to calculate the orifice plate measured pressure difference which is equivalent to the ISO 5167 correlation. Refer to the Flownex Library help on restrictors for the Flownex calculation method. However, the ISO 5167 correlation is more comprehensive and also provides methods to estimate flow rate uncertainties. 

According to ISO 5167 an orifice plate discharge may be calculated as follows:

𝑚̇ is the fluid discharge mass flow rate

𝐶 is the discharge coefficient

𝛽 is the ratio of orifice to pipe internal diameter (d/D )

𝜀 is the expansibility factor

𝑑 is the orifice diameter

∆𝑃 is the measured pressure difference across the orifice plate tappings

𝜌1 is the upstream fluid density

For a given (hypothetically fixed) mass flow rate and a desired measured pressure difference, the orifice plate diameter d is the required unknown. It may be more convenient, therefore, to reorganise (Eq. 1) in terms of the measured pressure difference:

Using the Flownex Designer, the orifice diameter d is varied until the required measured pressure difference ∆P is found for the given (fixed) mass flow rate 𝑚̇.

Since the flow rate is given as a volume flow at normal conditions (NTP), the equivalent mass flow has to be determined since flow rates are specified in terms of mass in Flownex and not volume flows. Unfortunately, like standard conditions (STP), normal conditions are not universally taken at the same pressure and temperature. In this case study, normal conditions are taken to be at 0°C and 101.325 kPa-a. The easiest method for finding the mass flow rate equivalent to the specified NTP volume flow rate is by employing the Designer. But in order to allow Flownex to correctly understand NTP conditions, the ambient conditions in the Flow Solver settings has to be set correctly:

Ensure that the elevation as well as the ambient temperature is set to zero as indicated. Then set up a Designer configuration which searches the upstream pipe “Volume flow based on ambient conditions” for a value of 1500 m³/h. The outlet boundary condition mass flow rate, which will be a negative value, is the independent variable in the search. The configuration for this case is shown in Figure 4. Running the Designer with this configuration set active will quickly find the equivalent mass flow rate and assign it to the outlet boundary.

Alternatively, the inlet boundary pressure and temperature values could also be temporarily set to normal condition values (0°C and 101.325 kPa-a) and the actual volume flow used in the Designer search.

With the mass flow rate known, the Designer is used again, this time to find the required orifice diameter whilst keeping the mass flow rate constant. The orifice plate result “Measured static pressure difference” is the equality constraint whilst the orifice plate input “Orifice diameter” forms the independent variable. Running the Designer with this configuration set as active arrives at the results shown in Figure 1 where an orifice diameter of 64.170 mm is calculated. Obviously, the actual orifice diameter will be chosen to represent a more standard size such as 64.0 mm or 64.5 mm. However, before a final decision can be taken, the orifice plate measured pressure drop has to be checked against the minimum flow case. Once again, the Designer is used to find the minimum flow.

As shown, the measured pressure difference across the orifice plate is only 0.26 kPa which is significantly less than the target value of 0.5 kPa. If accuracy at low flow is important, a pressure transmitter with the capability of accurately measuring a pressure difference range from 0.2 kPa to 5 kPa must be sourced.

Alternatively, if the pressure transmitter can handle pressure differences higher than 5 kPa, the orifice may be sized for the minimum flow case instead and tested against the maximum flow case. Using Designer, an orifice of 57.851 mm would be required to produce a 0.5 kPa pressure difference at 350 Nm³/h as shown in Figure 6. This orifice size, however, would produce a 9.7 kPa pressure difference at maximum flow as shown in Figure 7, so the pressure transmitter must be able to measure the range 0.5 kPa to 10 kPa.

Assuming that a more standard size such as a 57.5 mm orifice plate will be manufactured instead of the calculated 57.851 mm produced by the Flownex Designer, the calculations at minimum, normal and maximum flow must be repeated. The results are given in the following table.

Table 1: Flownex Results for the Actual Orifice Size of 57.5 mm

For the purposes of this case study, it is assumed that a pressure transducer and transmitter with measurements ranging from 0.5 kPa to 10 kPa are available. This concludes the sizing phase of the flow meter and leaves the uncertainty analysis to be performed.

ESTIMATION OF THE FLOW METER MEASUREMENT UNCERTAINTY 

When performing an uncertainty analysis the real value of a property is rarely known since it is deduced from the measured value and estimates of the uncertainties of measurement. For a given mass flow rate 𝑚̇  in (Eq. 2), ∆P may be calculated to provide the theoretical (real) orifice plate pressure difference value, however, several of the terms in the equation contain uncertainties. Therefore, the actual measured pressure difference would also contain some uncertainties which will be interpreted as uncertainties in the flow rate 𝑚̇ by the transmitter.

In this case study, the approach to uncertainty analysis is that the real values (flow rates) have been nominated at the start and are therefore known. The maximum flow rate is to be 1500 Nm3/hr, the nominal flow rate 1200 Nm3/hr and the minimum flow rate is 350 Nm3/hr. What is uncertain is the accuracy with which the orifice plate and transmitter combination is able to measure these real flow rates. Therefore in this case study, the aim is to determine the difference between the real and measured values. This difference is the uncertainty of measurement.

The case study focuses on two ways to determine this uncertainty:

• A mathematical approach where the governing equations are analyzed and an overall flow measurement uncertainty established. This approach demonstrates how quickly pure mathematical modeling will grow beyond the scope of most design engineers.

• Performing a Monte Carlo analysis using Flownex’s Sensitivity Analysis feature will demonstrate the simplicity of conducting complex uncertainty analyses which are still within the capability of any engineer equipped with Flownex.
   

Mathematical Uncertainty Estimation

With a known flow rate the theoretically measured pressure difference may be determined from (Eq. 2). However, the actual installed orifice plate flow meter may measure slightly different pressure differences to what theory suggests because of the following:

• The gas supply pressure may change.

• The gas supply temperature may change.

• The orifice plate has a manufacturing tolerance which will influence measured results.

• According to ISO 5167 there is some uncertainty in (Eq. 2).

• The gas composition may change.

• The pressure transmitter has less than perfect accuracy, stability and repeatability.

The influence of all of these variables combined or in isolation on the flow measurement uncertainty of the orifice plate may be estimated from first principles using the well-known root sum squared approach:

where: 

  is the combined resultant measured pressure difference uncertainty of all the independent uncertainties

      is the pressure-related uncertainty  

   is the temperature-related uncertainty  

   is the orifice plate manufacturing tolerance-related uncertainty

      is the ISO 5167 discharge coefficient uncertainty   (Eq. 2) 

    is the gas composition-related uncertainty  

    is the pressure transducer and transmitter-related uncertainty  

To be able to calculate each of the terms above, (Eq. 2) has to be evaluated in terms of changes of each of the terms in (Eq. 3). However, as can be seen, (Eq. 2) does not contain any direct pressure or temperature terms, but these uncertainties may be accounted for indirectly via the density. Similarly, changes in the gas composition could also be reflected as changes in density. The transducer/transmitter uncertainties are given by the manufacturer and can be readily added to the above equation. Orifice plate diameter tolerances are also known.

The measured pressure difference uncertainty 𝛿(∆𝑃) caused by each of these variables 𝜑 is established by determining the rate of change of the measured pressure difference with changes in each variable 𝜕(∆𝑃)⁄𝜕𝜑 and multiplying it by the uncertainty of that variable 𝛿𝜑. For the case at hand the resultant uncertainty then becomes:

Density Uncertainty

From (Eq. 2), the relationship between the measured pressure difference and each of the variables with uncertainties in question are as follows:

The term 𝛿𝜌 must be calculated from the influence of the pressure and temperature on the density. This may be done from the ideal gas equation of state:

where: 

Z is the ideal gas compressibility factor 

M is the molar mass of the gas 

is the universal gas constant 

Therefore: 

and thus 

Similarly: 

Assume that one could simply represent the expected changes in gas composition in terms of density changes by calculating the density from changes in the gas molar mass:

From the ideal gas law above: 

Typical natural gas has a molar mass between 17 kg/kmol and 20 kg/kmol. For this specific case, the natural gas composition used has a molar mass of 17.55 kg/kmol.  It would be reasonable to assume that the gas molar mass could vary between 17.35 kg/kmol and 17.75 kg/kmol, i.e.:

Combining (Eq. 6), (Eq. 7) and (Eq. 8):

The square of the product of (Eq. 5) and (Eq. 9) represents the first term in (Eq. 4).

Orifice Diameter Uncertainty

Since β = d / D

The term relates 𝛿𝑑 to the uncertainty in the orifice diameter:

The square of the product of (Eq. 10) and 𝛿𝑑 represents the second term in (Eq. 4).

Orifice Plate ISO Uncertainty

According to ISO 5167 an uncertainty exists for the discharge coefficient C in (Eq. 2). Therefore:

The discharge coefficient uncertainty as given by ISO 5167 for orifice plate beta ratios between 0.6 and 0.75 is:

For the case at hand the beta ratio (orifice to pipe internal diameter) is:

Since the diameter is larger than 71.12 mm and the Reynolds Number (Re_D) is larger than 10 000, the additional uncertainty correction terms specified by ISO 5167 are not applicable. Therefore:

Note however that (Eq. 13) may not be valid for other beta ratios, and that other uncertainty terms may be applicable at other Reynolds numbers.

The square of the product of (Eq. 7) and 𝛿𝐶 represents the third term in (Eq. 4).

Transducer and Transmitter Uncertainty

The transducer/transmitter uncertainty is given by the manufacturer directly in terms of a measured pressure difference uncertainty:

Total Combined Uncertainty

The above uncertainty terms must now be substituted into (Eq. 4):

However, this equation is still incomplete. Referring to (Eq. 2), the expansibility coefficient ε is a function of the upstream and downstream pressures P₁ and P₂ and the upstream density ρ₁ and should also have been included in the first two terms of (Eq. 16):

where

P₁ and P₂ are the upstream and downstream pressures respectively

k is the specific heat ratio

There are other factors too that have been omitted in the above discussion. It should be clear by now that even for this relatively simple problem, this approach would quickly escalate beyond the mathematical capabilities of most engineers in industry, including this author. Fortunately Flownex has the capability to perform the same uncertainty analysis numerically in a way that is easy to use for anyone familiar with Flownex.

Estimation of the Flow Meter Measurement Uncertainty using Flownex’s Sensitivity Analysis Feature

Flownex’s Sensitivity Analysis feature may be used to perform uncertainty analyses following two slightly different approaches:

• A series of parametric analyses may be performed, each analysis varying only one parameter such as those represented by (Eq. 3), and then combining their overall influence by using (Eq. 3) or (Eq. 16).

• Alternatively, all fluid flow related variables may be combined into a single Monte Carlo analysis, resulting in a single calculated uncertainty for these variables.

PARAMETRIC STUDY

As an example of the first approach, a parametric study of the influence of temperature alone on the flow measurement is demonstrated. Since the mass flow is kept constant, the density and therefore the actual volume flow rate has to vary in concurrence with the temperature. Variations in the volume flow rate will cause variations in measured pressure difference which the flow transmitter will translate as changes in flow rate. As shown in Figure 8, the temperature is varied between 5°C and 30°C with 25°C being the nominal value.

Running this configuration in the Sensitivity Analysis tool at the maximum flow rate produces the following results:

The orifice plate is supposed to measure 10.0331 kPa as discussed earlier, however, the values range from 9.319 kPa at 5°C to 10.215 kPa at 30°C which represents a total variation of:

Comparing this result to the mathematical approach given by (Eq. 5):

Even though these results differ by only 19 Pa, it should be noted that small variations in the expansibility coefficient and discharge coefficient will also result from changes in temperature since the volume flow rate will be changing. These have not been accounted for in the above equation. This exercise must now be repeated for each variable associated with an uncertainty and for each flow rate of interest.

MONTE CARLO ANALYSIS

However, all of these influences will occur simultaneously, and this is where Flownex’s Monte Carlo Analysis feature comes into its own. For the case at hand, the Monte Carlo configuration setup is as follows:

As before, the temperature is varied between 5°C and 30°C, the pressure is varied between 215 kPa-a and 225 kPa-a and the orifice diameter manufacturing tolerance is ±0.1 mm. For the purposes of demonstration, the Monte Carlo runs have been set to 1000 and could be set higher. The inputs are summarized in the following table:

Table 2: Flownex Inputs for the Monte Carlo Analysis

Even though the uncertainties 𝑢𝑃, 𝑢𝑇 and 𝑢𝑑  are supplied as three separate inputs, Flownex will calculate a single combined measured pressure difference uncertainty 𝑢𝑃𝑇𝑑. Note that the transmitter and discharge coefficient uncertainties cannot be added to the Flownex Monte Carlo run as they are not flow property related. Hence they must be added mathematically at a later stage. Also note that the uncertainty in gas composition has been ignored for the sake of simplicity.

Running this sensitivity analysis configuration allows Flownex to vary all three of these variables simultaneously, searching for the combination that would result in the maximum variation. The results are shown in the following figure:

These results have been copied into Excel and the minimum and maximum measured pressure difference values obtained. For this case, the maximum combined uncertainty for pressure, temperature and orifice diameter is:

The other uncertainties and must be evaluated mathematically and then added to using (Eq. 3):

From (Eq. 12) and (Eq. 14):

As given by the manufacturer in Table 2:

Therefore:

According to (Eq. 1) the uncertainty in measured pressure difference relates to an uncertainty in mass flow rate as follows:

Therefore, for the maximum flow case:

Similarly, a Monte Carlo analysis must also be performed at the minimum and nominal flow cases. The results are summarized in the following table:

Table 3: Flownex Results for the Monte Carlo Analysis

As shown in the results, the mass flow uncertainties are approximately 3.7%. If this level of accuracy is unacceptable, temperature and pressure correction may be implemented by adding pressure and temperature sensors to the upstream flow and adding appropriate logic to the transmitter to compensate for changes in the density.

Furthermore, orifice diameter manufacturing tolerances may also be improved. The results show that the bulk of the uncertainty in this case originates from the pressure, temperature and orifice diameter uncertainty 𝑢𝑃𝑇𝑑 – refer to Figure 12 below. A single parametric study on the influence of the manufacturing tolerance at maximum flow highlights that the contribution of the orifice tolerance to the overall uncertainty is approximately 14.8% of the total measured pressure difference uncertainty 𝑢∆𝑃. It can be concluded therefore that the bulk of the uncertainty originates from pressure and temperature related uncertainties.

Figure 12 shows the comparative uncertainties as a function of flow rate. As expected, the uncertainties increase with increasing flow rate. The transducer and transmitter related uncertainties are negligible in comparison, and the discharge coefficient uncertainty only contributes approximately 18% of the total. The bulk of the uncertainties are due to changes in pressure, temperature and the orifice plate diameter tolerance.

SUMMERY

For the case studied, it is shown that the bulk of the uncertainty originates from pressure and temperature variations. If the calculated mass flow uncertainty of 3.7% is not acceptable, pressure and temperature correction has to be implemented in the transmitter. The orifice plate manufacturing tolerance contributes less to the overall uncertainty but is relatively simple to improve.

Not only is Flownex capable of designing large, sophisticated flow networks, but it is also able to focus on a single element of design such as a simple orifice plate and then analyse that element in incredible depth. This versatility of Flownex is unparalleled amongst similar products and must make Flownex an indispensable software tool amongst process design engineers.

DOWNLOAD PDF