Monday, March 26, 2012

PSV SIZING FOR GAS OR VAPOR RELIEF

PSV SIZING FOR GAS OR VAPOR RELIEF
(Single Phase)
Critical Flow Behavior: If a compressible gas is expanded across a nozzle, an orifice, or the end of a pipe, its velocity and specific volume increase with decreasing downstream pressure.  For  a given set of upstream conditions the  mass  rate  of  flow  through  the  nozzle  will  increase until  a  limiting  velocity  is  reached  in  the  nozzle.  It can be shown that the limiting velocity is the velocity of sound in the flowing fluid at that location. The flow rate that corresponds to the limiting velocity is known as the critical flow rate.
Critical flow pressure: The absolute pressure ratio of the pressure at the nozzle exit at sonic velocity (Pcf) to the inlet pressure (P1) is called the critical pressure ratio. Pcf is known as the critical flow pressure.
If the inlet pressure (PSV relieving pressure) is P1
Outlet pressure at sonic velocity is Pcf
Than
Critical Flow Pressure rato =Pcf/P1
Normally critical flow pressure Pcf is 0.40 to 0.60 of inlet pressure P1, it’s a thumb rule.
Example: If inlet relieving pressure is 10 bar than critical flow will occur when the downstream pressure will reach 6 bar to 4 bar.
Note: Under critical flow conditions, the actual pressure at  the  nozzle  exit  of  the  pressure  relief  device  cannot  fall below the critical flow pressure even if a much lower pressure exists downstream(let say atm pressure). At critical flow, the expansion from nozzle pressure to downstream pressure takes place irreversibly with the energy dissipated in turbulence into the surrounding fluid.
Explanation: if the relieving pressure of PSV is 10 bar and critical flow pressure is 4 bar, than minimum pressure at outlet of PSV will be 4 bar. This means that if the back pressure is more than the pressure required at the downstream of the PSV for critical flow condition, at that moment of time you cannot utilize the full relief capacity of PSV (which occurs at critical flow condition). So you need to oversize        (oversize means ratio of actual orifice area to the required area at critical flow condition) the PSV orifice area for compensating the backpressure.

The critical flow pressure ratio in absolute units may be estimated using the ideal gas relationship:
where
Pcf   = critical flow nozzle pressure, in psia,
P1   = upstream relieving pressure, in psia,
k  =  ratio of specific heats for any ideal gas.
The sizing equations for pressure relief devices in vapor or gas service fall into two general categories depending on whether the flow is critical or subcritical.
1.       If the pressure downstream of the nozzle is less than, or equal to, thecritical flow pressure, Pcf, then critical flow will occur, and API 520-I  3.6.2 procedures should be applied.
2.       If the downstream pressure exceeds the critical flow pressure, Pcf, then subcritical flow  will  occur,  and  the API 520-I  procedures  in  3.6.3  or  3.6.4 should be applied.
3.6.2   Sizing for Critical Flow: Pressure  relief  devices  in  gas  or  vapor  service that operate  at  critical  flow  conditions  may  be sized using following Equations, Each of the equations may be used to calculate the effective discharge area, A, required to achieve a required flow rate through a pressure relief device. A  pressure  relief  valve  that  has  an  effective  discharge  area equal to or greater than the calculated value of A is then chosen for the application from API Std 526.

where
A  = required effective discharge area of the device,[in.2] .
W  = required flow through the device, [lb/hr].
C  =  coefficient determined from an expression of the  ratio of the specific heats (k = CP/Cv) of the gas or vapor at inlet relieving conditions. This can be obtained from Figure 32 or Table 8. Where k cannot be determined, it is suggested that a value of C equal to 315 be used.


Kd   = effective coefficient of discharge. For preliminary sizing, use the following values:
= 0.975 when a pressure relief valve is installed with or without a rupture disk in combination,
= 0.62 when a pressure relief valve is not installed and sizing is for a rupture disk.
P1   = upstream relieving pressure, psia. This is the set pressure plus the allowable overpressure plus atmospheric pressure(14.7).
Kb   = capacity correction factor due to back pressure.
This can be obtained from the manufacturer’s literature or estimated for preliminary sizing from Figure30. The back pressure correction factor applies to balanced bellows valves only.
 For conventional and pilot operated valves, use a value for Kb equal to 1.0

Kc   = combination correction factor for installations with a rupture disk upstream of the pressure relief valve
= 1.0 when a rupture disk is not installed,
= 0.9 when a rupture disk is installed in combination with a pressure relief valve and the combination does not have a published value.
T = relieving temperature of the inlet gas or vapor, R (°F + 460).
Z = compressibility factor for the deviation of the actual gas from a perfect gas, a ratio evaluated at inlet
relieving conditions.
M = molecular weight of the gas or vapor at inlet relieving conditions.
V  = required flow through the device, scfm at 14.7 psia and 60°F.
G  = specific gravity of gas at standard conditions referred to air at standard conditions [normal conditions]. In other words, G = 1.00 for air at 14.7 psia and 60°F.
3.6.3   Sizing for Subcritical Flow, Gas or Vapor (Single Phase):
(Conventional and Pilot-Operated Pressure Relief Valves)
When the ratio of back pressure to inlet pressure exceeds the critical pressure ratio Pcf /P1, the flow through the pressure relief device is subcritical. Following Equations may be used to calculate the required effective discharge area for a conventional pressure relief valve that has its spring setting adjusted to compensate for superimposed back pressure; equations may also be used for sizing a pilot-operated relief valve.
US Customary Units:
where
A  = required effective discharge area of the device,[in.2].
W  = required flow through the device,[lb/hr].
F2   = coefficient of subcritical flow, see Figure 34 for values or use the following equation:

k = ratio of the specific heats.
r = ratio of back pressure to upstream relieving pressure, P2/P1.
Kd   = effective coefficient of discharge. For preliminary sizing, use the following values:
= 0.975 when a pressure relief valve is installed with or without a rupture disk in combination,
= 0.62 when a pressure relief valve is not installed and sizing is for a rupture disk.
Kc   = combination correction factor for installations with a rupture disk upstream of the pressure relief valve.
= 1.0 when a rupture disk is not installed,
= 0.9 when a rupture disk is installed in combination with a pressure relief valve and the combination does not have a published value.
Z = compressibility factor for the deviation of the actual gas from a perfect gas, evaluated at relieving inlet conditions.
T = relieving temperature of the inlet gas or vapor, R (°F + 460).
M = molecular weight of the gas or vapor.

P1   = upstream relieving pressure, psia. This is the set pressure plus the allowable overpressure plus atmospheric pressure (14.7).
P2   = back pressure, psia.
V = required flow through the device, scfm at 14.7 psia and 60°F.
G = specific gravity of gas at standard conditions referred to air at standard conditions                     (normal conditions). In other words, G = 1.00 for air at 14.7 psia and 60°F.
3.6.3.3   Balanced Pressure Relief Valves
Balanced  pressure  relief  valves  should  be  sized  using Equations  3.2  through  3.4  in  paragraph  3.6.2(critical flow equations). 
The back pressure correction factor in this application accounts for flow velocities that are subcritical as well as the tendency for the disc to drop below full lift (the use of subcritical flow equations are appropriate only where full lift is maintained). The back  pressure  correction  factor,  Kb,  for  this  application should be obtained from the manufacturer.

Wednesday, March 21, 2012

More accurate Method to calculate the relieving temperature for Fire case PSV sizing by Shell DEP method.

OBJECTIVE: More accurate Method to calculate the relieving temperature for Fire case PSV sizing by Shell DEP method.
(This method is more accurate because in this method we consider the actual volumetric flow for relieving condition; this method gives somewhat less value of relieving temperature as compared to the Method 1 which is based on mass flow)
PFD: PFD developed in Aspen HYSYS
2.     Methodology: Following are the steps to be followed for developing the model in HYSYS.
a.     Open HYSYS and define the stream component, Select fluid package then go to simulation environment and pick a material stream from HYSYS utility panel.
b.     Define the stream conditions Temperature as Operating temperature and operating Pressure give some mass flow rate let say 100 kg/hr.
c.      Once stream (20) is conversed now put a Flash drum (V-105) downstream of material stream (20).
d.     Now define two outlet stream 21 & 22, 21 stream as vapor outlet stream and 22 as liquid outlet stream
e.     Define the vessel volume, vessel volume includes vessel volume on which PSV is mounted and piping volume which comes under fire zone, if piping volume is not known add 10% additional volume in the vessel volume.
f.       Now define the liquid level in vessel V-105 by adjusting the percentage, and note down the liquid volume (VL) and gas volume (VG) in the vessel.
g.     Split stream 21 into stream 23 & 24 and stream 22 into 26 & 27.
h.     Now adjust stream 24 mass flow rate into actual volume flow rate, select user supplied option in Adjuster connection option and specified target value is VG   (m3/hr).
i.       Define some flow rate in stream number 23, HYSYS will adjust stream 24 according to the specified target value, if you are not able to adjust properly than you should see the adjuster parameter and define tolerance, step size minimum value, maximum value and maximum number of iterations than reset adjuster it will adjust.
j.       Now adjust stream 26 mass flow rate into actual volume flow rate, select user supplied option in Adjuster connection option and specified target value is VL (m3/hr).
k.     Define some flow rate in stream number 27, HYSYS will adjust stream 26 according to the specified target value, if you are not able to adjust properly than you should see the adjuster parameter and define tolerance, step size minimum value, maximum value and maximum number of iterations than reset adjuster it will adjust.
l.       Now put a MIX-102 and mix stream 24, 26 and get the final mixed stream 28 from the mixture (MIX-102) outlet.
m.  Now define another stream 29, keeping mass balance same as stream 28 and pressure as Relieving Pressure.
n.     Put a Adjuster for stream 29 to adjust the temperature of stream 29 in such a way that mass density of stream number 29 is equal to the mass density of stream 28.
o.     After adjusting the temperature of stream 29.
p.     Now note down the temperature of stream 29 (Relieving Temperature) which is at relieving pressure.

Result:
1.      The Relieving Temperature at relieving Pressure is 191.8 deg C and on mass basis it was 194.6 deg C in method1.
Conclusion: Both the methods (1 & 2) is correct and according to SHELL philosophy acceptable to Shell but when you are going to analyze a Flare network, if you have some existing PSV size for Fire case at that time in Flare-net you only need to get the relieving temperature to get the value of rated capacity (Calculated by Flarenet) of Existing known size PSV at that time you can prefer METHOD1 and if you are going to size new PSV for Fire case at that time you can use METHOD1.

To calculate the Relieving temperature & Latent Heat of Vaporization for Fire case PSV sizing by Shell DEP method.

OBJECTIVE: To calculate the Relieving temperature & Latent Heat of Vaporization for Fire case PSV sizing by Shell DEP method.
(This method is based on mass flow it gives somewhat higher relieving temperature)
PFD: PFD developed in Aspen HYSYS
1.     Methodology: Following are the steps to be followed for developing the model in HYSYS.
a.     Open HYSYS and define the stream component, Select fluid package then go to simulation environment and pick a material stream from HYSYS utility panel.
b.     Define the stream conditions Temperature as Operating temperature and Pressure as Relieving Pressure give some mass flow rate let say 100 kg/hr.
c.      Once stream (14) is conversed now put a heater (E-100) downstream of material stream (14).
d.     Now define the vapor fraction of stream (17) for 5% flashing by mole if the composition basis is mole % and after that go to heater parameter and define pressure drop across heater zero, now E-100 and stream 17 is conversed.
e.     Now put an Adjuster on stream 17 to adjust vapor fraction 5% by mole basis to convert into 5% weight basis.
f.       Now place a Flash drum (V-104) for flashing stream no 17.
g.     Define two outlet streams 15 & 16 for V-104, 15 for vapor outlet and 16 for liquid outlet from the flash drum.
h.     If the stream 16 vapor fraction is zero then automatically stream 16 is at its bubble point, therefore got into the properties of stream 16 and note down the mass enthalpy (-1969KJ/Kg).
i.       Now place a heater E-102 downstream of stream16 to bring the stream at its dew point.
j.       Define a heater E-102 outlet stream 18, now define the pressure drop across heater zero and stream vapor fraction 0.00001.
k.     Now stream heater E-102 and stream 18 is conversed and stream 18 is at Dew point.
l.       Note down the temperature of stream 18, which your relieving temperature (194.8 deg C) at relieving pressure.
m.  Now got to the stream 18 properties note the value of mass enthalpy (-1610 KJ/Kg).
n.     Subtract the liquid mass enthalpy (-1969 KJ/kg) from vapor mass enthalpy (-1610KJ/Kg), you will get the latent heat of vaporization 358.6 KJ/Kg.

Result:
1.     The Relieving Temperature at relieving Pressure is 194.6 deg C.
2.     Latent Heat of Vaporization is 358.6 KJ/Kg.