Orifice Flow for Subcooled Flashing Liquids
In a previous Fireside Chat, Orifice Flow Calculation Basis, we provided information on the flow of compressible fluids through sharp-edged orifices. While that topic was primarily geared toward the flow of gases through an orifice (for vapor blowdowns), we had briefly mentioned flow of a two-phase fluid through an orifice and provided a reference to Diener and Schmidt1 (DS Method). We’ve encountered problems implementing this method for flashing liquids (liquids having a saturation or bubblepoint pressure less than the inlet pressure yet greater than the outlet pressure) and wanted to provide an update.
With the flow of two-phase or flashing liquid through an orifice, there are significant non-equilibrium effects, mostly attributed to the time scale associated with vapor nucleation relative to the flow across the restriction. For orifices, in which there isn’t much of a run up distance for the acceleration of the incoming fluid, this results in essentially frozen flow. For us old-timers, we can recall the use of the frozen flow orifice calculations in the API RP 520 Part 1 5th Edition for two-phase flow or flashing liquid flow. The reason those equations were used is that the results actually fit the orifice flow data well – it really wasn’t until DIERS came along and tested nozzles instead of orifices that any significant discrepancy was widely recognized.2 The data for nozzles showed a substantially lower flow rate than what was being predicted by frozen flow, which prompted the implementation of the homogeneous equilibrium method for nozzle flow. This lower flow rate has been attributed to the occurrence of gas formation within the nozzle, especially for nozzles greater than 10 cm (4”) in length.3 It is interesting to note that the data for lowly subcooled or low quality two-phase flow show that there is still likely some thermal non-equilibrium (boiling delay) happening in these cases.3, 4
While working on orifice flow calculations for flashing liquid flow, we ran into difficulties matching the calculation results from the DS Method to experimental results from Richardson, et al,5 prompting us to reconsider the use of the DS Method. Various references state that subcooled flashing liquids flowing through orifices essentially behave like nonflashing liquids, but for us there has always been a question about the effective pressure differential. Using the data of Richardson, et al, we have convinced ourselves that the effective pressure differential is based on the pressure downstream of the orifice (as opposed to the saturation pressure or some other effective pressure, which is what we might use for highly subcooled liquids flowing through a nozzle), regardless of the degree of subcooling of the liquid. In the figure below, we plot the flow data of Richardson for subcooled flashing liquid cases versus the saturated pressure ratio. To generate a useful two-axis plot, we use the dimensionless parameters of Leung6 for the normalized mass flux and saturation pressure ratio. The data tables of Richardson provided fluid composition, inlet pressure (P0), inlet temperatures, orifice diameters, and measured flow rates. The NIST REFPROP program7 was used to generate the inlet densities (ρ0) and saturation pressures (Psat).
Normalized mass flux versus saturation pressure ratio
With this basis, we can obtain flow rates that are very close to the experimental results from Richardson, et al. The average discharge coefficient for subcooled flashing liquid flow from Richardson is 0.594, with a standard deviation of 0.01, which corresponds to the average ASME PTC 19.5-2004 liquid orifice discharge coefficient of 0.596. A comparison of our calculation results to the experimental data is seen in the figure below.
Calculated flow versus experimentally measured flow
Incidentally, the work of Richardson et al was performed for the IP (now Energy Institute) Guidelines for the Safe and Optimum Design of Hydrocarbon Pressure Relief and Blowdown Systems.8
In summary, for subcooled flashing liquids, in which the vapor pressure is greater than the back pressure but less than the incoming pressure, we use the liquid orifice flow calculation with the following parameters, regardless of the degree of subcooling:
- Flow basis = liquid orifice flow 9,10
- Upstream pressure = incoming pressure
- Upstream fluid density = density at the saturation pressure
- Downstream pressure = actual back pressure
- Discharge coefficient = liquid discharge coefficient
[1] Diener R, Schmidt J. “Sizing of Throttling Device for Gas/Liquid Two-Phase Flow Part 2: Control Valves, Orifices and Nozzles”. Process Safety Progress. March 2005; 24 (1): 29-37.
[2] Fisher HG, Forrest HS, Grossel SS, Huff JE, Muller AR, Noronha JA, Shaw DA, Tilley BJ. Emergency Relief System Design Using DIERS Technology – The Design Institute for Emergency Relief Systems (DIERS) Project Manual. AIChE (1992).
[3] Fauske, HK. “Flashing Flows or: Some Practical Guidelines for Emergency Releases”. Plant/Operations Progress. 1985 July; 4(3): 132-134.
[4] Sallet DW, Somers GW. “Flow Capacity and Response of Safety Relief Valves to Saturated Water Flow”. Plant/Operations Progress. 1985 October; 4(4): 207-216.
[5] Richardson SM, Saville G, Fisher A, Meredith AJ, Dix MJ. “Experimental Determination of Two-Phase Flow Rates of Hydrocarbons Through Restrictions”. Trans IChemE, Part B, Process Safety and Environmental Protection, January 2006; 84(B1): 40-53.
[6] Leung JC. “Two-phase flow discharge in nozzles and pipes – a unified approach”. J. Loss Prev. Process Ind. 1990 January; 3: 27-32.
[7] Lemmon, E.W., Huber, M.L., McLinden, M.O. NIST Standard Reference Database 23: Reference Fluid Thermodynamic and Transport Properties-REFPROP, Version 9.1, National Institute of Standards and Technology, Standard Reference Data Program, Gaithersburg, 2013.
[8] Institute of Petroleum. Guidelines for the Safe and Optimum Design of Hydrocarbon Pressure Relief and Blowdown Systems. 2001; London.
[9] Boyce MP. “Transport and Storage of Fluids – Head Meters”. In Perry RH and Green DW. Perry’s Chemical Engineers’ Handbook. 1997; New York, McGraw Hill: 10.11-10-18.
[10] Darby R. Chemical Engineering Fluid Mechanics. 2001; Boca Raton, CRC Press: 304-305.
1 Comments
The selection of sizing method to be used for a two-phase flow application is a decision made by the designer. It is therefore vitally important that the designer has detailed knowledge of the process conditions. The use of an appropriate software tool can also help greatly with the relief valve orifice sizing task.
Post a Comment