Two-Phase Flashing Flow Methods and Comparisons: J. C. Leung and F. N. Nazario
Two-Phase Flashing Flow Methods and Comparisons: J. C. Leung and F. N. Nazario
Two-Phase Flashing Flow Methods and Comparisons: J. C. Leung and F. N. Nazario
comparisons
Two-phase flashing flow methods used by DIERS, the American Petroleum Institute (API) and
ASME are reviewed and compared. For saturated and low quality inlet conditions, the API
method was found to yield higher theoretical mass flow rates than the DIERS homogeneous
equilibrium flow model (HEM). The degree of disparity can be quantified by two dimensionless
physical property groups. Comparison with flashing flow data for water, refrigerant 12 and a
propane-butane mixture shows good agreement with the HEM prediction, but poor agreement
with the API method. The ASME method only recommends a saturated water flashing flow curve
that is in agreement with the HEM calculation to within 10%.
Two-phase flashing flows are frequently encountered in experimental data. The purpose of this paper is to
both practice and design, ranging from process flow review and evaluate these existing flashing flow proce-
streams through metering devices to pressure relief dures as recommended by the API, ASME4, and
sizing application. The most significant feature of DIERS organizations. This paper, however, will not be
flashing flow is the ease with which choking or critical concerned with the application of various flow/capacity
flow condition is attained. This is a maximum flow coefficients as employed in an actual safety relief valve
situation in which further lowering of the downstream sizing. Rather, the discussion will focus on the theore-
pressure will not lead to further increase in flow. For tical aspect of the flashing flow determination. It is
vapour or gas flow in nozzles or short pipes, the hoped that such a discussion will bring out the
pressure at choking is about half of the upstream differences as well as similarities among these currently
chamber (stagnation) pressure. However, for saturated employed methodologies.
liquid and low quality (vapour mass fraction) inlet
conditions, the pressure at choking is typically only Review of flashing flow methods
l&20% below the upstream pressure. Thus it does not DtEKS recommendation
take much flashing as a result of the depressurization DIERS recommends the HEM as the most appropriate
and fluid acceleration before the flow is ‘choked’, quite flashing flow formulation*. In this model, the flashing
unlike the gas case. two-phase flow mixture is treated much like the
Previously. most of the investigations on flashing classical compressible gas while undergoing an adiaba-
choked flow were related to the so-called ‘loss-of- tic expansion with equal velocities and temperature
coolant’ accident in nuclear reactor safety studies, and (thermodynamic equilibrium) in both phases. Among
therefore water was the fluid being studied most. Much the many other flow models tested in the DIERS
of these earlier works dealt with thermodynamic program, the HEM model yields conservative (low)
non-equilibrium and non-homogeneous (slip) flow estimates of the flow capacity. In flow passages of
behaviour in the attempt to explain the observed greater than 0.1 m in length (frequently encountered in
departure from the reference homogeneous equilib- relief systems), the HEM model in fact provides a
rium flow model (HEM) prediction. In the chemical best-estimate calculation. This observation is consistent
processing and petroleum refining industries, a more with previous investigations5-s.
widely practiced flashing flow method used in relief For flow in frictionless devices, the traditional
valve sizing is based on the API procedure I. The method for evaluating the HEM flow rate would
recently completed DIERS research recommends the require seeking a maximum in the mass velocity
HEM approach?J. which is supported by DIERS expression derived from the first law of
thetmodynamic9, namely,
Received 6 June 1989
Presented at the AlChE Loss Prevention Symposium. 2-6 April
G = Mho - hII”’
1989, Houston, TX, USA V (1)
095c4230/90/020253-08
0 1990 Butterworth &Co. (Publishers1 Ltd
where ho denotes the stagnation enthalpy, and the greater than 8.0) by
specific enthalpy and specific volume are expressed in
terms of quality and the individual phase properties,
G* - 0.9 (44
i.e. h, = hl, + xohf,,, h = hf + xh,, and u =
6
or simply
UJ + XU,,. In addition, the static quality is to be
(1.5
evaluated along a constant entropy line (isentropic
path)
(4b)
x = (S” - St)/SQ (2) A PI recommendation
The above procedures, while not complicated, are API Recommended Practice 521 provides only a brief
rather tedious and would require extensive thermo- description of how to handle two-phase flashing flow
dynamic property tabulation. The latter poses a special through relief devices. That text is repeated here for
problem in fluids without such extensive property completeness:
tables. The DIERS research had suggested a novel way ‘A pressure relief valve handling a liquid at
of carrying out the calculational procedure while vapour-liquid equilibrium or a mixed-phase fluid
reducing the number of flash computations’ will produce flashing with vapour generation as
substantially2.3. More recently, Leungr” proposed, for the fluid moves through the valve. The vapour
the first time, a generalized HEM correlation applic- generation may reduce the effective mass flow
able to all fluids including multicomponent mixtures capacity of the valve and must be taken into
with ideal solution behaviourlJ. The correlation is account. The quantity of flash vapour is common-
based on a single property parameter w given by ly calculated based on adiabatic flashing from the
relieving condition - either to the critical down-
. stream pressure (the critical pressure drop across
Vtl “0 f higo 1
Here all properties are evaluated at the known up- the valve orifice) or to the back pressure,
stream stagnation conditions (with the vapour prop- whichever is higher - and on an orifice area
erties for the case in which the fluid is a saturated liquid calculated for that vapour flow using the same
being based on the composition of the initial vapour pressure drop. An area is also calculated for the
produced when the liquid flashes). This correlation is remaining liquid quantity (after flash) under
reproduced in Figure 1, in which the normalized critical relieving conditions using the total pressure drop
D u”, and the critical pressure (the relieving pressure minus the back pressure).
mass flux. GT = G,/m
ratio, rlc = P,/P,, are correlated solely in terms of w. The orifice selected should have an area equal to
Curve fit equations were given in the original paperr”. or greater than the sum of the individual areas
For the case of saturated liquid inlet, the first term previously calculated’.
in the w definition (which can be approximated by the We could not locate the original paper that
inlet void fraction (Y(,) drops out and the critical mass discussed this procedure and its supporting data,
flux can be approximated (for fluids with o values although numerous recent papers have made reference
to the API flashing flow method1:-r5. In what follows,
we will attempt to give our interpretations of this
model. Our best understanding is that it is a frozen
separated (stratified) two-phase flow formulation with
an inlet quality evaluated based on an adiabatic flash
from the upstream relief condition to the choked
pressure of the vapour phase. Thus it can be characte-
rized as a two-step model, not necessarily rigorous in
the thermodynamic sense.
f1
xcv 1 -X,, -’
w,IA,+- JJ’,/A
and the total flow area, from Equation (A-9)
A =A,+A, (7)
Although these two equations are quite different
in form, they in fact represent the same continuity
statement. Since Equation (7) is stated explicitly in the
API method. the alternate statement, Equation (6).
would also be implied.
Now the individual mass flux terms, W,/A, and 0 400 800 I200 1600 2WO 2400 2000 3200
W,/A f, in Equation (6) can be derived theoretically as Set pressure(psig~
detailed in Appendix B. This appendix also contains Figure 2 ASME flashing flow correlation for saturated water
the rather lengthy conversion of these formulae to the
form used in APJ. These formulae are demonstrated to
(berg)
be equivalent under the following two conditions: 0 50 100 150 200
I I I I iSO
0 The API method in reality permits two different 24 -
pressures to co-exist at the choked plane, with the
E
vapour stream at its own choked pressure as if -409
flowing alone and the liquid stream at the back E
P
pressure. This is therefore a simplifying assump-
tion. - 36%
To summarize then, the API method would Figure3 Comparison of HEM prediction (---_) versus ASME
(---) flashing flow curve for saturated water
evaluate the flashing mass flux according to Equation
(6) with x,, calculated as described earlier, and the
liquid and vapour mass fluxes, WI/A, and W,/A,. pressure ratio is calculated in the ASME method. But
estimated from Equations (B-2) and (B-3) in Appendix in the vicinity of the maximum flow evaluated accord-
B. respectively. ing to Equation (l), the value of G decreases only
gradually with either increasing or decreasing pressure
ASME recommendation ratio. This might help to explain the observed discre-
The ASME Code does not address the issue of flashing pancy if indeed the calculated critical pressure ratio in
two-phase flow except for the case of saturated water. the ASME method does not correspond to the max-
It presents a flashing flow curve for saturated water as imum in G.
reproduced in Figure 2. According to the ASME
publication4, this curve is the theoretical, isentropic Saturated water data
flow value arrived by assuming equilibrium flow and Upon examining choked flow data for saturated water
calculated values for the critical pressure ratio. No from various sources, Moody5 found that equilibrium
recommendation is provided for other fluids. however. flow conditions were closely approached in pipes longer
than 0.1 m. On the other end, pipes approaching zero
length (orifice) tend to have non-flashing Bernoulli-
Comparison
type flow rates. For the first time, Moody demon-
ASME curve and HEM strated that these choked flow data with L > 0.1 m
The ASME saturated water flow capacity curve is first could in fact be predicted quite well simply with the
compared with the published HEM tabulationi6. This HEM. Figure 4 reproduces his comparison of the HEM
comparison, as illustrated in Figure 3, shows reason- with extensive literature data17-*l satisfying this equilib-
able agreement to within 10% throughout the entire rium length criteria. Though Moody’s evaluation of the
pressure range. However, it is not clear how the critical HEM was based on Equation (1) and the steam table,
A
0 I I I I I I z”t
0 2 4 6 e IO 12 I
StoQnotion pressure, PO t MPO)
.- L
IR I
Figure4 Comparison of saturated water choked flow data
-in I6
Y -
against HEM and API (case 1) methods (adapted from MoodyS): E
A, Ref. 17, 4.0 mm diameter; o, Ref. 18, 6.35 mm diameter; n , : l4-
Ref. 19, 7.95 mm diameter; 0, Ref. 20; 173.0 mm diameter; 0,
Ref. 21.13.1 mm diameter 5
x I2 -
2 IO-
his prediction can be reproduced closely by the HEM B
generalized correlation of Figure 1. Also shown in this E 8-
figure are the predictions from the API flashing flow E
2
E 6-
method. Two cases are shown:
4-
Case I. This case employs KD = 1.0 for the vapour
stream and KD = 0.62 for the liquid stream. (Note that 2-
API actually recommends using 0.975 for the vapour I
01 I I I I I
K,, but the results are indistinguishable). 0 0.2 0.4 0.6 0.8 1.0
Stagnation pressye,Pc WPo)
Case 2. This is a modified API method. The choked Figure 5 Comparison of saturated R-12 data against HEM and
pressure for the liquid stream also corresponds to the API (case 1) methods: R-12 data of van den Akker et a/.22, 0,
vapour stream, i.e. the available pressure drop for the stainless steel pipe, o, polycarbonate pipe
Table 1 Comparison with Cp/CI flashing flow data of van den Akker (private communication, 1985). Discharge of 50-5Owt%C&
mixture in a 2 mm dia. x 150 mm length duct
Case 1 Case 2
D 2.150.1 859 9.1 0.144 4800 5300 8655 7203
+ 10% + 80% + 50%
D 2.150.2 916 8.7 0.148 5700 5520 8951 7448
- 3% f 57% + 31%
D 2.150.3 917 8.7 0.148 5600 5525 8960 7456
- 1% + 60% + 33%
D 2.150.4 911 8.8 0.148 5700 5500 8908 7417
- 4% + 56% + 30%
“A discharge coefficient of 0.9 was used to account for the frictional losses in the 75 L/D duct
bBoth cases employed a K, of 0.62 in Equations B-2 or B-8; case 1 assumed AP = PO - Pb while case 2 assumed AP = PO - P,,. After G
was evaluated according to Equation (61. a factor of 0.9 was multiplied to account for the frictional losses
properties. Table 2 illustrates that the predictions quality x,, in Equation (8) can be related to the ‘vapour
based simply on single-component fluid (i.e. either C3 phase’ critical pressure ratio ncV and w, via a modified
or C,) at the average relief pressure and the corres- form of the energy balance (see Appendix C)
ponding component saturation temperature yield for all
intents and purposes the same result. XC” = (1 - n,.,+ (12)
Combining Equations (8), (11). and (12) yields the
Generalized comparison
desired result
To cast the comparison into a more generalized form,
G HEM
we present a closed-form expression relating the HEM -=
flow rate to the API flow rate for the saturated inlet
case. First the mass flux corresponding to the API
method can be written as
G API =
c;G
XC”
+KoG;\/p,p,
l-x,,
1
-’
(8) This equation can be readily evaluated once the
(13)
where G; and G; are the normalized mass fluxes given parameters w, pr/plO, k ( va p our specific heat ratio) and
P,,/Po are specified. In practice, this ratio G HEM/G Art
is governed primarily by two dimensionless property
groups. w and pr/pgo. Figure 6 illustrates how this
quantity G HEM/G API varies with the density ratio
CT= [2(1 - PJP,,)]‘/Z = [2(1 - r7b)]‘/? (10) pf/+, for the conditions of w = 5 or 10, k = 1.2 and
where 77s is the back pressure to stagnation pressure P,/P,, = 0.1. For the two cases shown, i.e.: API
ratio. As for the HEM mass flux, recall method with K, = 0.62 and AP given by the total
pressure drop for the liquid stream; modified API
G HEM = G&M s (11) method with KD = 0.62 for the liquid stream and AP
where G&rM is given by Equation (4a). The ‘critical’ given by (PO - P,“) only, the ratio GueM/GAri is
w G(kg m 2s ‘1
Pll
Test
no. (kPa) Case A Case B Case C Case D Case A Case B Case C Case D Meas.
D2.150.1 859 9.06 9.70 8.12 8.28 5300 5150 5454 5468 4800
+ 10% + 7% + 14% + 14%
D2.150.2 916 8.73 9.39 7.90 7.99 5520 5360 5662 5690 5700
- 3% - 6% - 1% 0%
D2.150.3 917 8.72 9.38 7.89 7.98 5525 5365 5667 5695 5600
- 1% - 4% + 1% + 2%
D2.150.4 911 8.76 9.42 7.91 8.01 5500 5340 5643 5672 5700
- 4% -.6% ~ 1% ~ 1%
Case A, h,,, from rigorous thermodynamic calculation involving C&.,; Case B. h,, from IY,h&, where Y, is vapour mass composition:
Case C. one component assumption based on C,; Case D, one component assumption based on C,
Appendix B: individual mass tlux evaluation where gpm denotes flow rate in gallons per minute.
in API method Again the above equation is in US customary engineer-
ing units (A B is in in2 here). For comparison purposes,
The individual mass flux terms, W,/A, and W,/A I, in
we shall ignore the influence of the discharge coeffi-
Equation (6) can be evaluated separately as follows.
cient K, and the various capacity correction factors,
Under the assumption of frictionless non-interacting
K,, K,, and K,. It is further noted that (1.25 P, - P,,)
flow, the two streams are going to follow isentropic
is simply the available pressure drop A P on the basis of
paths while satisfying the flow area requirement of
25% overpressure for liquid relief sizing. In such cases,
Equation (7). Thus the vapour will expand isentropic-
K, is assigned a value of unity. Now, the Bernoulli
ally in a nozzle according to the well-known equationZh
equation, Equation (B-2), can be rearranged to yield
1 ( Wf/P,)
Af=Kn m
(B-6)
(B-1) and by employing units consistent with the API practice
where a discharge coefficient KD is included to account as shown in the equation below (note that the 6 is
for imperfect nozzles. replaced by 6 here)
Since the liquid is incompressible in the frozen
flow assumption, it will obey Bernoulli’s equation. Its Af(in’)(
li5y;n2)
=
mass flux at the place where the pressure is P will be
1
W, (gpm)(3.785 X 1O-3 m3/gal) kg/m’ ‘I2
- = KD[2pf(Po - P)]‘i” (B-2) sg 2 lOoO-
A, (60 min/sec) ]O
It is widely known that the discharge coefficient K, is
Ko[2AP (psi)(6897~)]“’ “”
about 0.62 for a sharp-edge orifice and this value is
recommended by API as a conservative approach. (B-7)
Theoretically the choked pressure is to be found by
We finally obtain
satisfying the criterion that at the throat dA/dz is zero.
This would yield a choked pressure less than the typical gpmVs&/cc)
A f(in2) = (B-8)
vapour choked pressure given in Step 1. Such a solution 38K&w-.
is, as expected, much more involved. The API method which is identical to Equation (B-5) with K,, K,, and
recommends instead that P in Equation (B-l) for the K, being 1.0, and (1.25P, - Ph) being APon the basis
vapour flow be limited by Equation (5), but P in of 25% overpressure.
Equation (B-2) for the liquid flow be evaluated at the
downstream back pressure P,,. In other words, the total
available pressure drop is allowed in the liquid phase.
Thus for the vapour stream, the classical choked
flow formulaZh is obtained by combining Equations (5)
and (B-l)
Appendix C: approximate solution for flash
quality
(B-3) Based on a constant enthalpy energy balance, the
adiabatic flash quality at a temperature T from an
This is to be compared with Equation (1) in API-520 initially saturated liquid state (x,, = 0) is
[Ref. 27) for gas/vapour relief
C,(To - T)
M X= (C-1)
W, = CKDA~P,, -.L h fE
(B-4)
ZT,,
We invoke the Clapeyron relation
where C is given by
AI-’
-=- h,
,,,[,( k : l )i-y AT (C-2)
or,T
and approximate the properties at the stagnation
and the number 520 contains the proper conversion condition. thus arriving at
when using the US customary engineering units (A g is
x= __-
in in2 here) in Equation (B-4). By noting that pa0 is (C-3)
P,,M,/ZRT,,, Equations (B-3) and (B-4) are there-
and for us,, >> vf
fore equivalent in form.
As for the liquid stream, Equation (B-2) with AP ego
x=xX (C-4)
P = Pb is to be compared with Equation (C-7) in API
RP-520 publication Therefore, at the vapour phase critical pressure
gpm%i
Af= (B-5)
38K,K,K,K,d 1.25P, - P,,