Depressurisation: A Practical Guide Updated For Version 3.2 / 2004
Depressurisation: A Practical Guide Updated For Version 3.2 / 2004
Depressurisation: A Practical Guide Updated For Version 3.2 / 2004
This guide has been prepared based upon questions frequently asked regarding the Dynamic
Depressuring utility introduced in HYSYS 3.0. It should provide users with an explanation how to use
the utility and correctly interpret the results. It is divided into three sections:
1.0 Overview
4.0 Appendices
1.0 Overview
Why has the old depressuring utility gone?
The original utility in HYSYS was a pseudo-dynamic calculation based on a series of steady state
calculations. The Dynamic Depressuring utility was introduced in HYSYS 3.0 to allow users to perform
proper time-dependant calculations. A HYSYS Dynamics licence is NOT required to use this new
utility. In version 3.2 onwards you now only have the option to run the new Dynamic Utility.
Make sure the template dyndepressuring.tpl is dated 19/04/2004 or later
The Depressuring utility can be used to simulate the depressurisation of gas, gas-liquid filled vessels,
pipelines and systems with several connected vessels or piping volumes depressuring through a
single valve. References to “vessel” in this guide can also refer to piping or combinations of the two.
• Adiabatic Mode is used to model the blowdown of pressure vessels or piping with no external
heat supplied.
A more in depth discussion of the different methods follows in Section 2.0.
A Depressuring utility can be added to the case by selecting "Tools" "Utilities", highlighting
"Depressuring - Dynamics" and pressing the "Add Utility" Button. Once selected always rename the
Utility to something that is recognisable the next time you open the case. Eg: DP-V1234-Fire
2.1 Connections
How to connect the utility to a stream
On the "Design" tab, "Connections" page, choose the stream that represents the fluid you want to use
as the source for the depressuring. If you have a single vessel, for example, the stream would be the
feed stream into the vessel. Attaching the stream to the utility is accomplished as shown in the view
Ideally, the vessel size will be known and this data can be entered into the appropriate fields on the
form shown above.
If the vessel size is unknown, then the vessel sizing utility in HYSYS can be
used to estimate the required parameters.
The initial liquid volume is normally calculated at the normal liquid level (NLL). HYSYS does not take
the volume of the heads into account so the volume will be the liquid in the cylindrical portion only. If
the feed stream is two-phase, the equilibrium composition of the liquid will be calculated. If an initial
HYSYS does not take account of the heads in a vessel so volumes and areas are calculated as for a
simple cylinder. The total vessel volume is calculated from the diameter and height (or length for a
horizontal vessel). To account for piping or head volume contributions, a small amount can be added
to the height or length of the vessel if the user so requires..
If the conditions of the system at settle out are such that the vapour is superheated, HYSYS will not
allow a liquid inventory. The settle out conditions for mixed sources and volumes are calculated on a
constant enthalpy, volume and mass basis.
Correction Factors allow for adjustments to the amount of metal in contact with the top or bottom of
the vessel. This can also be used to account for additional nozzles, piping, strapping or support
steelwork in close contact with the vessel. HYSYS will use the heat content of this metal when
performing the calculations. This is analogous to adding, for example, ten percent to the vessel mass
to account for fittings.
• Fire Mode
can be used to simulate plant emergency conditions that would occur during a plant fire. Pressure,
temperature and flow profiles are calculated for the application of an external heat source to a
vessel, piping or combination of items. Heat flux into the fluid is user defined using the following
equation:
LiquidVolumetime=t
Q = C1 + C2 × time + C3 (C4 − TVESSEL ) + C5 ×
LiquidVolumetime=0
The Fire equation can also be used to simulate the depressuring of sub-sea pipelines where heat
transfer occurs between seawater and the pipeline. If C3 was equal to UA, C4 was equal to T1 and
C1, C2 and C5 were equal to zero, the above equation would reduce to:
Q = UA(∆T )
uses similar heat flux parameters to those used in Fire mode. Three coefficients: C1, C2 and C3
must be specified. The equation used by HYSYS is an extension to the standard API equation for
heat flux to a liquid containing vessel. A wetted area is required and used to calculate the heat
transfer into the vessel.
The following notes are based on extracts from Guide for Pressure-Relieving and Depressuring
System, API Recommended Practice 521, Fourth Edition, March 1997.
The amount of heat absorbed by a vessel exposed to an open fire is affected by:
The following equations are based on conditions where there is prompt fire fighting and adequate
drainage of flammable materials away from the vessel.
Environmental Factor
Table 5 on Page 17 of API 521 lists F factors for various types of vessels and insulation.
Bare vessel, F = 1.
Earth-covered storage, F = 0.03.
Below-grade storage, F = 0.
For insulated vessels, users should consult the reference and select an F value based on the
insulation conductance for fire exposure conditions.
At present HYSYS does not have the F factor shown in the equation. If you want to
derate the heat input then alter the 21000 or 43.116 figure accordingly.
Note the Initial wetted area variable will only be completed if cases from earlier versions of HYSYS
( pre 3.2) are loaded.
The HYSYS equation is an extension of the standard API equation. Therefore, in field units, C1 will be
21000 multiplied by the environmental factor, F and C2 will 0.82.
(In most cases, C1 will be equal to 21000).
C2
Q = C1 × (WettedAreatime=t )
The surface area wetted by the internal liquid content of the vessel is effective in generating
vapour when the exterior of the vessel is exposed to fire. To determine vapour generation API
recommends that you only take into account that portion of the vessel that is wetted by liquid up to
7.6m (25ft) above the source of the flame. Individual companies may deviate from this figure, so
do check. This usually refers to ground level but it can be any level capable of sustaining a pool
fire. The following table indicates recommended volumes for partially filled vessels. Volumes
higher than 7.6m are normally excluded as are vessel heads protected by support skirts.
The following table is an example showing how the C3 term affects the wetted area calculation. An
initial liquid volume of 6m3 and a wetted area of 500 m2 were given.
C3
Liquid Volume 1 0.75 0.5 0.25 0
Time Volume Ratio Wetted Area Wetted Area Wetted Area Wetted Area Wetted Area
(minutes) (m3) (m2) (m2) (m2) (m2) (m2)
0 6 1.0 500.0 500.0 500.0 500.0 500.0
5 4 0.7 333.3 375.0 416.7 458.3 500.0
10 3 0.5 250.0 312.5 375.0 437.5 500.0
15 2 0.3 166.7 250.0 333.3 416.7 500.0
Therefore if a C3 value of 0 is used, the initial wetted area is used throughout the calculations. This
could represent a worst case scenario. Alternatively, if a C3 value of 1 was used, the volume would
vary proportionally with the liquid volume. This would represent a vertical vessel.
uses the Boltzmann constant to take into account radiation, forced convection, flame temperature
and ambient temperature. The method may be considered as an alternative method to the API
standard.
( ( )
Q = Atotal × ε f × ε v × k (T f + 273.15 ) − (TV + 273.15 ) + outsideU × (Tamb − TV )
4 4
)
where:
A total = total wetted surface area
f = flame emissivity generally ranges from 0.2 to 0.5
(for burning heavy HCs)
v = vessel emissivity generally ranges from 0.5 to 1 (for polished metal)
k = Boltzman constant equals 5.67*10 - 8 W/m2 K4
Tf = flame temperature 1500 K and upwards
Tv = vessel temperature
outside U = convective heat transfer between vessel and air
Tamb = ambient air temp
• Use Spreadsheet is an option that allows the user access to the spreadsheet used by the
depressuring utility. Values can be altered in this spreadsheet and additional equations substituted
for calculation of the heat flux. It is recommended that this option only be used by advanced
users.
Notes
2. Simple: allows the user to either specify the heat loss directly or have it calculated from specified
values
3. Detailed: allows the user to specify a more detailed set of heat loss parameters
For accurate calculations the detailed Heat Loss model is the one that
should be used otherwise the contribution of the metal is completely ignored.
HOWEVER, when using the API Fire mode you should NOT use a heat loss model. The API
Fire equations were developed on the assumption that the equation caters for all effects.
Simple Model
• Heat Transfer Area is the cylindrical area of the vessel with no allowance for head area. This
value is calculated using the vessel dimensions specified on the "Connections" page.
• Using the Simple Heat Loss Model, heat loss from the vessel is calculated using the following
formula:
Q = UA(T fluid − Tambient )
General
The General section allows the user to manipulate Recycle Efficiencies and the ambient temperature.
The default value for all three Recycle Efficiencies is 100%. This means that all material in the vessel
has been flashed together and is in thermodynamic equilibrium. If the Recycle Efficiencies were to be
reduced a portion of the material would by-pass the flash calculation and the vapour and liquid would
no longer instantaneously reach equilibrium. In this case, the phases may have different temperatures.
Unfortunately, there is no single typical number suggested for these parameters. The best option
would be to try various scenarios and observe the results.
Conduction
The Conduction parameters allow the user to manipulate the conductive properties of the wall and
insulation.
The metal wall thickness must always have a finite value (i.e.: it cannot be <empty>). To model a
vessel without insulation, the insulation value thickness should be zero. Users are also required to
enter the specific heat capacity of the material(s), the density of the material(s) and the conductivity of
the material(s).
Convection
The Convection view allows users to manipulate the heat transfer coefficient for inside and outside the
vessel as well as between vapour and liquid material inside the vessel.
To use a set of fixed U values, the "Use Fixed U" option should be selected. If the U values are
unknown, the user can press the "Estimate Coefficients Now" button and have HYSYS determine the
U values. In order to have HYSYS vary the U values throughout the depressuring scenario, select the
"Continually Update U" value.
Correlation Coefficients
This feature gives users the opportunity to manipulate the coefficients used in the heat transfer
correlation. By selecting "Use Specified Constants", the user may manually enter the constants used
in the heat transfer correlations.
m
∆T
h =C×
length
The equation used for the other three correlations is:
m
Nu = C × (Gr × Pr )
Where: Nu = Nusselt Number
Gr = Grashof Number
Pr = Prandtl Number
Notes
The Valve Parameters page allows users to select the type of valves to be used for both vapour and
liquid service. In most cases, either the Fisher or the Relief valve should be used for valve sizing.
Their equations are more advanced than some of the others and can automatically handle choked
conditions. Furthermore, these two valve types support other options that can be accessed through
the valve property view accessible through the Depressuring sub-flowsheet. The seven available valve
types are described in the sections that follow.
Fisher
The Fisher option uses the standard valve option in HYSYS. It allows the user to specify both valve Cv
and percent opening. By pressing the "Size Valve", the valve can be sized for a given flow rate.
PSV hysterysis can be modelled by clicking on the Sub flowsheet template for the depressuring utility
and navigating to the Specs Dynamics Tab as shown.
It is possible to model a depressuring valve using the PSV valve. To do this you will
have to force the relief valve to be open at all times. To do this enter a full open pressure
that is lower than the final expected vessel pressure and a set pressure that is slightly
lower than the full open pressure.
For other valve information as used in the original HYSYS Depressuring Valve Criteria See
Appendix A
“ PV Work Term Contribution" refers approximately to the isentropic efficiency of the process. A
reversible process should have a value of 100% and an isenthalpic process should have a value of
0%. As no processes are fully isentropic nor isenthalpic, this parameter is used in all the different
simulation models to tune the models in order to match conditions observed, and has been requested
by some of our users to use to match the test data they have available.
This parameter is defined in HYSYS as: dH (change in enthalpy) = percentage / 100.0 * dP (change in
pressure) / phase_mole_density. In simple terms you can think of this as the work that the fluid in the
vessel does to expel the material that is leaving.
However for design purpose, that is working without any test data, based on various publications on
the subject, we can propose the following values as a guideline:
Furthermore as you can see, from the way the equation is defined, a higher isentropic efficiency
results in a lower final temperature. Hence, if one is checking that the minimum temperature of the
vessel will not fall below a certain value i.e. for validating the steel alloy grade, then 100% will give the
most conservative result. Also, If one is checking that the final Pressure is below the safety regulatory
limit after 15 minutes, might be safer to makes some checks with lower value such as 87%, to be
more conservative, provided there is no significant heat transfer influence on the phase behaviour
inside the vessel.
Operating pressure refers to the initial vessel pressure. By default, this value is the pressure of the
inlet stream. The time step size refers to the integration step size. It may be a good idea to reduce the
step size if the flow rate is significantly larger than the volume or if the vessel depressurises in a
relatively short amount of time , say 1 minute.
Either the Dynamic Depressuring utility can solve for the final pressure or the Cv/Area required to
achieve a specified final pressure.
The "Calculate Pressure" option uses the specified area/Cv to determine the final pressure.
"Calculate Area" is available for Relief, Supersonic, Subsonic and General valves. "Calculate Cv" is
available for Fisher and Masoneilan valves. The two options differ only in the type of value calculated.
Based on API, it is normal to depressure to 50% of the staring pressure or to 100 psig. Before the
calculations start, the user must specify an initial Cv or area. If the depressuring time is reached
before the final pressure is achieved, then the calculations stop and a new Cv or area is calculated
using the final pressure. The calculations are repeated until the final pressure is reached in the given
amount of depressuring time. The user may specify a maximum number of iterations and a pressure
tolerance to improve convergence. If the user wishes to stop the calculations at any time, the <CTRL>
<BREAK> keys can be used.
Once the utility has run, users can go to the "Performance" "Summary" page to view the results.
• You now only have the option to run the new Dynamic Utility the old quasi dynamic
method has been removed.
• HYSYS does not take the volume of the vessel heads into account so the volume will
be the liquid in the cylindrical portion only. Adjust the vessel size if you wish to allow
for the head volume
• HYSYS defaults the liquid volume to be equal to the volumetric flow of the feed liquid
over one hour. This will be disproportionate to the total vessel volume, a more
realistic hold up time is 4 minutes.
• At present HYSYS does not have the F factor shown in the API521 fire equation. If
you want to derate the heat input then alter the 21000 or 43.116 figure accordingly.
• It is possible to model a depressuring valve using the PSV valve. To do this you will
have to force the relief valve to be open at all times. To do this enter a full open
pressure that is lower than the final expected vessel pressure and a set pressure
that is slightly lower than the full open pressure.
• PV work term
* gas-filled systems, 87% to 98%.
* liquid filled systems 40% to 70%.
A higher efficiency results in a lower final temperature.
• For accurate calculations the detailed Heat Loss model should be used otherwise
the contribution of the metal is completely ignored.
• Make sure you run with a smaller enough time step to capture the peak flow
Supersonic
The supersonic valve equation can be used for modelling systems when no detailed information on the
valve is available. The discharge coefficient (Cd) should be a value between 0.7 and 1. P1 refers to the
upstream pressure and ρ1 the density.
0.5
F = C d × A × (P1 × ρ 1 )
Subsonic
The subsonic valve equation can also be used for modelling systems when no detailed information on
the valve is available but the flow is sub-critical. This can occur when the upstream pressure is less
than twice the backpressure. The discharge coefficient (Cd) should be a value between 0 and 1. The
area (A) should be a value between 0.7 and 1. P1 refers to the upstream pressure and ρ1 the density.
0.5
F = Cd × A ×
(P1 + Pback )× (P1 − Pback ) ρ
1
P1
Pback refers
Back Pressure
It is possible to have the depressuring scenario cycle between pressure build-up and relief. To perform
this analysis, ensure a reasonable pressure differential and increase the number of pressure steps.
Masoneilan
This equation was taken from the Masoneilan catalogue. It can be used for general depressuring
valves to flare. When this option is selected, the user must specify Cv and Cf. The remaining
parameters in the equation are set by the Depressuring utility.
0.5
F = C1 × Cv × C f × Y f × (P1 × ρ1 )
where:
The General valve equation is based on the equation used to calculate critical flow through a nozzle as
1
shown in Perry's Chemical Engineers' Handbook . It should be used when the valve throat area is
known. Note that this equation makes certain limiting assumptions concerning the characteristics of
the orifice.
F = Cd × Av × K term × ( g c × P1 × ρ1 × k )
0.5
Cd = discharge coefficient
Av = throat cross sectional area
k +1
Kterm = 2 2 ( k +1 )
k +1
k = ratio of specific heats (Cp/Cv)
P1 = upstream pressure
ρ1 = upstream density
1
Page 5-14, Equation 5.20 (6th Edition) & Page 10-15, Equation 10.26 (7th Edition)
No Flow
This option indicates that there is no flow through the valve.
Use Spreadsheet
Recommended only for advanced users, this option allows the user to customise a valve equation by
editing the valve spreadsheet found inside the Depressuring sub-flowsheet.
Discharge Coefficient
When the relief, supersonic, subsonic or general valve is selected, the user is required to specify a
discharge coefficient. This correction factor accounts for the vena contracta effect. Values ranging
from 0.6 to 0.7 are typically used. In order to disregard this effect, set the discharge coefficient equal to
1.