AD 121 Reference Guide
AD 121 Reference Guide
Reference Guide
Version: 12.1
August 2003
Copyright © 2003 by Aspen Technology, Inc. All rights reserved.
AspenTech product names that are mentioned in this manual, e.g., Aspen Dynamics™, Aspen Custom Modeler®,
Aspen Plus™, Enterprise Optimization and the aspen leaf logo are trademarks or registered
trademarks of Aspen Technology, Inc., Cambridge, MA.
All other brand and product names are trademarks or registered trademarks of their respective companies.
This manual is intended as a guide to using AspenTech's software. This documentation contains AspenTech
proprietary and confidential information and may not be disclosed, used, or copied without the prior consent of
AspenTech or as set forth in the applicable license agreement. Users are solely responsible for the proper use of the
software and the application of the results obtained.
Although AspenTech has tested the software and reviewed the documentation, the sole warranty for the software
may be found in the applicable license agreement between AspenTech and the user. ASPENTECH MAKES NO
WARRANTY OR REPRESENTATION, EITHER EXPRESSED OR IMPLIED, WITH RESPECT TO THIS
DOCUMENTATION, ITS QUALITY, PERFORMANCE, MERCHANTABILITY, OR FITNESS FOR A
PARTICULAR PURPOSE.
Corporate
Aspen Technology, Inc.
Ten Canal Park
Cambridge, MA 02141-2201
USA
Phone: (1) (617) 949-1000
Fax: (1) (617) 949-1030
URL: http://www.aspentech.com/
Contents
Introduction ix
About This Guide............................................................................................................... ix
Associated Guides.............................................................................................................. ix
What You Need To Use This Guide ....................................................................................x
For More Information ..........................................................................................................x
Online Help..............................................................................................................x
Printing the Portable Document Format Files (PDFs)........................................... xi
World Wide Web ................................................................................................... xi
Technical Support .................................................................................................. xi
1 Streams 1-1
MaterialStream in Dynamic Simulations......................................................................... 1-2
Reverse Flow for Material Stream in Dynamic Simulations ............................... 1-2
Specifying MaterialStream Feeds for Dynamic Simulations .............................. 1-3
Using the Stream Sensor...................................................................................... 1-3
HeatStream in Dynamic Simulations............................................................................... 1-5
Specifying HeatStream Feeds for Dynamic Simulations..................................... 1-6
WorkStream in Dynamic Simulations ............................................................................. 1-6
Specifying WorkStream Feeds for Dynamic Simulations ................................... 1-6
3 Separators 3-1
Flash2 in Dynamic Simulations ....................................................................................... 3-2
Reverse Flow for Flash2 in Dynamic Simulations .............................................. 3-4
Specifying Flash2 for Dynamic Simulations ....................................................... 3-4
Flash3 in Dynamic Simulations ....................................................................................... 3-8
Reverse Flow for Flash3 in Dynamic Simulations ............................................ 3-10
Specifying Flash3 for Dynamic Simulations ..................................................... 3-11
Decanter in Dynamic Simulations ................................................................................. 3-15
Reverse Flow for Decanter in Dynamic Simulations ........................................ 3-17
5 Columns 5-1
Distl in Dynamic Simulations.......................................................................................... 5-2
Specifying Distl for Dynamic Simulations .......................................................... 5-2
References for Distl ......................................................................................................... 5-2
Extract in Dynamic Simulations ...................................................................................... 5-3
Specifying Extract for Dynamic Simulations ...................................................... 5-3
RadFrac in Dynamic Simulations .................................................................................... 5-4
Pressure/Flow Effects for RadFrac ...................................................................... 5-5
Pressure and Level Control for RadFrac.............................................................. 5-8
Specifying RadFrac for Pressure Relief............................................................. 5-13
Specifying RadFrac Overall Heat Transfer Coefficient .................................... 5-13
Specifying RadFrac for Dynamic Simulations .................................................. 5-15
Changing RadFrac Operating Specifications..................................................... 5-20
Simulating Column Startup ............................................................................... 5-22
References for RadFrac ..................................................................................... 5-23
PetroFrac in Dynamic Simulations ................................................................................ 5-24
Pressure/Flow Effects for PetroFrac .................................................................. 5-26
Pressure and Level Control for PetroFrac.......................................................... 5-28
Specifying PetroFrac for Pressure Relief........................................................... 5-32
Specifying PetroFrac Overall Heat Transfer Coefficient .................................. 5-32
Specifying PetroFrac for Dynamic Simulations ................................................ 5-34
References for PetroFrac ................................................................................... 5-39
Column Hydraulic and Pressure Drop Equations .......................................................... 5-40
12 Procedures 12-1
Procedure pChairX ........................................................................................................ 12-1
Mathematical Description.................................................................................. 12-1
Chair Function ................................................................................................... 12-2
Exces Function................................................................................................... 12-2
Usage ............................................................................................................................. 12-2
Liquid Case ........................................................................................................ 12-2
Vapor Case......................................................................................................... 12-3
What is the Availability (Av)?........................................................................... 12-3
Index
Associated Guides
Aspen Dynamics also comes with these other guides:
Operating Hours
Phone Numbers
This chapter describes the unit operation models for mixers and
splitters:
• Mixer
• FSplit
For more details of these methods, see the Aspen Dynamics User
Guide, Chapter 1.
Entrainment in Flash2
In Aspen Dynamics, liquid entrainment in Flash2 is specified
differently than in Aspen Plus, when the dynamics of the flash
vessel are included. This is because the Aspen Plus entrainment
specification assumes steady state, which is not applicable to a
dynamic flash because the liquid flow and vapor flow out of the
vessel are independent of each other, so the level and pressure in
the vessel can vary.
In Aspen Plus, the liquid entrained in the vapor is specified as the
liquid molar flow in the vapor stream over the total of liquid molar
flow, where the total liquid molar flow is given by:
[molar flow through the liquid outlet] plus
[molar flow of liquid that leaves via the vapor outlet]
For an instantaneous flash, Aspen Dynamics uses the same
specification. The variable ENTR is fixed and matches the Aspen
Plus specification.
Reverse Flow for Flash3 is supported in reverse flow simulations when the Vessel
Flash3 in Dynamic type is Vertical or Horizontal. Flash3 is not supported when the
Simulations Vessel type is Instantaneous.
When configuring your simulation for reverse flow, use the
RFlowPhase parameter on the Flash3 Configure form to specify
which phase in the Flash3 will flow out through the feed stream(s)
if the pressure gradients cause the flows in those streams to
reverse. You can select the vapor phase (Vapor), the 1st liquid
phase (Liquid1), the 2nd liquid phase (Liquid2), a homogeneous
liquid phase mixture (Liquid) or a homogeneous three phase
mixture (Mixed) with the same phase fractions as in the Flash3.
The default value is Mixed.
When reverse flow is enabled, material flow may enter or exit
Flash3 via any of its connected streams, depending on the pressure
gradients around the flowsheet.
Note If you specify a pressure drop for Flash3 in Aspen Plus, the
inlet stream pressure will not equal the block pressure and, for the
dynamic simulation, the inlet pressure change will be fixed to the
steady state value. Upon export to Aspen Dynamics, the Pressure
Entrainment in Flash3
In Aspen Dynamics, liquid entrainment in Flash3 is specified
differently than in Aspen Plus, when the dynamics of the flash
vessel are included. This is because the Aspen Plus entrainment
specification assumes steady state, which is not applicable to a
dynamic flash because the liquid flow and vapor flow out of the
vessel are independent of each other, so the level and pressure in
the vessel can vary.
In Aspen Plus, the liquid entrained in the vapor is specified as the
liquid molar flow in the vapor stream over the total of liquid molar
flow, where the total liquid molar flow is given by:
[molar flow through the liquid outlet] plus
[molar flow of liquid that leaves via the vapor outlet]
Where:
K = A constant determined by fitting to steady-
state conditions
∆P = Pressure drop
Rho = Mass density at outlet conditions
Fv_out = Outlet volumetric flow rate
If the steady-state simulation predicts a pressure rise, the pressure
drop equation is not used. Instead the pressure difference will be
specified.
Reverse Flow for Heater is supported in reverse flow simulations. No additional
Heater in Dynamic specifications are required to configure Heater for reverse flow.
Simulations When reverse flow is enabled, Heater is modeled as an
instantaneous mixer followed by a single input, single outlet flow
resistance with heat transfer. Therefore, depending on the local
pressure gradients during simulation, it is possible for some of the
feed streams to be in forward flow and some to be in reverse flow,
so bypassing the flow resistance/heat transfer section of model and
not receiving the Heater duty.
Notes:
• When using Heater in reverse flow simulations, it is
recommended that you do not use Constant Duty for the Heat
Transfer Option. Instead, it is better to use either the Constant
Temperature or LMTD Heat Transfer Option in Heater to
model the heat transfer in a more realistic manner. The reason
for this follows:
With a fixed heat duty, as the flow rate becomes small, the
process temperature is driven above its upper bound towards
infinity. To avoid this, an in-built procedure in Aspen
Dynamics attempts to ramp the fixed duty to zero to keep the
temperature within its bounds. However, when the flow
reverses, the possibility of a discontinuity in the temperature
can lead to convergence problems with this procedure.
• It is recommended that you use Calculated for the Pressure
Drop Option when using Heater in reverse flow simulations. If
you use Fixed for the Pressure Drop Option, the pressure drop
For more details of these methods, see the Aspen Dynamics User
Guide, Chapter 1.
Pressure Heat
drop here transfer
Cold side
Where:
Where:
U_ref = Reference overall heat transfer coefficient
(value fixed by user)
Flow = Hot-side flowrate
Flow_ref = Reference hot-side flowrate (value fixed by
user)
Where:
U_ref is Reference overall heat transfer coefficient
(value fixed by user)
Flow is Cold-side flowrate
Flow_ref is Reference cold-side flowrate (value fixed by
user)
Uexp is exponent (value fixed by user)
The flow can be either on a molar or mass basis determined by the
state of the HTCBasis parameter.
Where:
U_ref is Reference overall heat transfer
coefficient (value fixed by user)
Flow_hot is Hot-side flowrate
Flow_hot_ref is Reference hot-side flowrate (value fixed by
user)
Uexp_hot is Hot-side exponent (value fixed by
user)
Flow_cold is Cold-side flowrate
Flow_cold_ref is Reference cold-side flowrate (value fixed by
user)
Uexp_cold Cold-side exponent (value fixed by
user)
Where:
h is Film heat transfer coefficient
D is Diameter
KW is Tube metal thermal conductivity
Subscripts:
SHELL : Shell-side
TUBE : Tube-side
o : Outer tube
i : Inner tube
and
u exp
Flow
h = h _ ref .
Flow _ ref
If you have data on the dynamics of your plant, you can adjust
these volumes to give a good match to your data.
For dynamic heat exchangers you can choose to model equipment
heat capacity and heat transfer with the environment.
Pressure Heat
drop here transfer
Cold stream
Pressure/Flow Effects for MHeatX
For each stream, the pressure drops are related to the outlet
volumetric flow rates by:
∆P = K ⋅ Rho ⋅ Fv _ out 2
Where:
K = A constant determined by fitting to steady-
state conditions
∆P = Pressure drop
Rho = Mass density at outlet conditions
Fv_out = Outlet volumetric flow rate
If you are using the dynamic model type, all pressure drop is
assumed to take place between the inlet and outlet volumes.
If the steady-state simulation predicts a pressure rise, the pressure
drop equation is not used, and the pressure difference for that
stream will be specified.
Reverse Flow for MHeatX is supported in reverse flow simulations. No additional
MHeatX in Dynamic specifications are required to configure MHeatX for reverse flow.
Simulations When reverse flow is enabled, the flow direction used for the
calculation of the logarithmic mean temperature difference
(LMTD) in each MHeatX stream is adjusted during dynamic
simulation, based on the current hot and cold side flow directions.
For example, if you specify a Countercurrent flow direction in
Aspen Plus, and during a reverse flow simulation the Hot Side
flow reverses, MHeatX will then use a Co-current method for the
calculation of the LMTD. Note that the exchanger configuration
Pressure and Level Pressure and level controllers are automatically added to the
Control for RadFrac flowsheet when appropriate.
For flow-driven simulations:
Reflux Drum
LC
Reflux Drum
LC
Reflux Drum
LC
ILC
Liquid 1 Distillate
Liquid 1 Reflux
LC
Bottoms
Default sump inventory control configuration
Constant Method
U is fixed and may be changed directly by the user
Where:
U_ref is Reference overall heat transfer coefficient
(value fixed by user)
Flow is Actual flowrate
Flow_ref is Reference flowrate (value fixed by user)
Uexp is exponent (value fixed by user)
For the process side the flow can be either on a molar or mass
basis determined by the state of the HTCBasis parameter. The
medium side is always be on a mass basis.
Where:
h is Film heat transfer coefficient
D is Diameter
KW is Tube metal thermal conductivity
Subscripts:
SHELL : Shell-side
TUBE : Tube-side
o : Outer tube
i : Inner tube
and
u exp
Flow
h = h _ ref .
Flow _ ref
For RadFrac you can choose to model equipment heat capacity and
heat transfer with the environment for the Reflux Drum, Sump and
Stages.
The modeling approach lumps the reboiler volume with the sump
volume. You should allow for this when entering the sump size.
For more information about vessel geometry, see the Aspen
Dynamics User Guide, Chapter 1.
If one or more liquid phases are present, specify the initial liquid
volume fraction in the column sump. The default value is 0.5.
If two liquid phases are present and there are separate product
streams for each liquid phase, you can specify the initial liquid 1
volume fraction in the column sump. If you do not specify a value,
then the model calculates the initial liquid 1 volume using the
steady state liquid ratio (beta).
Note You do not need to specify this information for the top and
bottom stages of the column when these correspond to the
condenser and reboiler. However, if there is no condenser or
reboiler, you must complete the specifications for these stages.
Note:
• If you wish to model column startup you should use the simple
tray or packing hydraulics correlations (hyperlink to Selecting
the Method for Hydraulic and Pressure Drop Calculations).
The rigorous correlations are not suitable for conditions where
there is zero flow in the column.
References for Perry, R.H. and Green, D., (1984), Perry’s Chemical Engineers’
RadFrac Handbook, 6th edn, pp 18-10, Eqn (18-19) and pp 18-7, Eqn (18-
3).
Stichlmair, J., Bravo, J.L., and Fair, J.R., (1989), “General Model
for Prediction of Pressure Drop and Capacity of Countercurrent
Gas/Liquid Packed Columns” In: Gas Separation and Purification,
Vol 3, pp 19-28, Eqns (14) and (15).
Where:
K = Constant determined by fitting to steady-state conditions
∆P = Pressure drop
ρ = Mass density at outlet conditions
Fv_out = Outlet volumetric flow rate
Pressure/Flow Effects The pressure in the top stage, whether it is a condenser or a normal
for PetroFrac tray, is always determined from the conditions on that stage, that is
the temperature and the composition.
For other stages:
• Pressure difference between stages is determined using a
pressure drop flow equation.
• Flow rates from stages are related to the liquid holdup on the
stage using a hydraulics equation.
Pressure and Level Pressure and level controllers are automatically added to the
Control for PetroFrac flowsheet when appropriate.
For flow-driven simulations:
Reflux Drum
LC
Reflux Drum
LC
LC
Bottoms
Default sump inventory control configuration
Specifying PetroFrac To specify the overall heat transfer coefficient (U) for the
Overall Heat Transfer condenser
Coefficient • Open the Configure form for the relevant submodel of the
instance of the PetroFrac model on the flowsheet.
The overall heat transfer coefficient can be calculated by these
different methods:
• Constant
• Power law relationship based on process-side
• Power law relationship based on heating/cooling medium-side
• Film coefficient method
The method used is controlled by the HTCOpt parameter.
Constant Method
U is fixed and may be changed directly by the user
Where:
U_ref is Reference overall heat transfer coefficient
(value fixed by user)
Flow is Actual flowrate
Flow_ref is Reference flowrate (value fixed by user)
Uexp is exponent (value fixed by user)
For the process side the flow can be either on a molar or mass
basis determined by the state of the HTCBasis parameter. The
medium side is always be on a mass basis.
Where:
h is Film heat transfer coefficient
D is Diameter
KW is Tube metal thermal conductivity
Subscripts:
SHELL : Shell-side
TUBE : Tube-side
o : Outer tube
i : Inner tube
and
u exp
Flow
h = h _ ref .
Flow _ ref
The modeling approach lumps the reboiler volume with the sump
volume. You should allow for this when entering the sump size.
For more information about vessel geometry, see the Aspen
Dynamics User Guide, Chapter 1.
If a liquid phase is present, specify the initial liquid volume
fraction in the column sump. The default value is 0.5.
Where:
QL = Volumetric liquid flow rate from the stage
KWeir = Weir constant
The value used for the weir constant is the same as that used in the
rigorous tray rating methods.
LWeir = Total weir length
hCrest = Height of the liquid crest over the weir
The total weir length is specified using the ratio of weir length to
column diameter. The default for this ratio is 0.7267.
The liquid crest is the difference between the height of liquid on
the tray and the weir height. The height of liquid is the ratio of the
volume of liquid to the active area of the tray. In turn, the active
area of the tray is specified as a percentage of the tray area. The
default value for this percentage is 90%.
• To simulate a single pass tray, use the default ratio of weir
length to column diameter (that is, 0.7267), and the default
value of active area as a percentage of the tray area (that is,
90%).
• To simulate two, three or four pass trays, adjust these two
default values.
Where:
∆P = Pressure drop per theoretical stage
KO = Orifice constant
ρV = Vapor molar density
QV = Volumetric flow rate of vapor to the packed section
Where:
K dc = A constant
QL = Volumetric flow rate of liquid from the tray for that
pass
Ada = Minimum area under the downcomer apron
The values of the constant, K dc , and the minimum area under the
downcomer apron, Ada , are the same as those used in the rigorous
tray rating calculations in the steady-state simulation.
The downcomer backup is given by the following equation (Perry,
1984):
hb = hWeir + hCrest + htot + hdc
Where:
hWeir = Height of the weir
hCrest = Liquid crest over the weir
htot = Total liquid head loss across the tray
Where:
K BL = A dimensionless constant, equal to 0.555.
uL = Superficial liquid velocity
a = Specific surface area of the packing
g = Acceleration due to gravity
ε = Packing voidage
∆P = Pressure drop across the section
Z = Height of packing in the section
ρL = Liquid density
The specific surface area, a, the packing voidage, ε, and the height
of packing in the section, Z, are all obtained from the steady-state
tray and packing calculations.
Where:
Fl Drain = Draining liquid molar flow rate
Where:
Fv = Vapor flow rate through orifices
Fv Norm = Normal (steady-state) vapor flow rate
through orifices
C DrainPoint = Drain point coefficient, set to 0.1
The two coefficients may be adjusted in the tray model to fit the
draining behavior of real trays.
This chapter first describes the reaction kinetics that are available
in Aspen Dynamics for rate-based reactions.
It then describes the following unit operation models for reactors:
• RStoic
• RYield
• RGibbs
• RCSTR
• RPlug
Adsorption expression
(∑ K ∏ C )i
νj m
j
Option Description
Stirred tank (CSTR) Perfect mixing is assumed.
The reaction equations are applied at the outlet conditions of the
reactor, so that the reactor outlet always satisfies the reaction
specification.
You can use a vertical or horizontal vessel. Liquid and vapor holdup
in the vessel are modeled and the liquid level can vary during a
dynamic run. You can also choose the heating and cooling methods.
The reactions are instantaneous.
Tubular or plug flow The reaction equations are applied at the outlet conditions of the
reactor (PFR) reactor, so that the reactor outlet always satisfies the reaction
specification.
The reactions are instantaneous.
Instantaneous No dynamic effects are modeled.
PFR Geometry
For PFR geometry, specify the length and internal diameter of the
reactor.
For more information about vessel geometry, see the Aspen
Dynamics User Guide, Chapter 1.
For PFR geometry, specify the length and internal diameter of the
reactor.
For more information about vessel geometry, see the Aspen
Dynamics User Guide, Chapter 1.
For PFR, specify the length and internal diameter of the reactor.
For more information about vessel geometry, see the Aspen
Dynamics User Guide, Chapter 1.
You can specify either the vessel diameter or the vessel length.
The other dimension is calculated from the volume, or residence
time specification entered for the steady-state simulation.
Specifying the When simulating a pressure relief system, the RCSTR model
Pressure Relief enables you to specify the following:
RCSTR Model • Reactions
• Geometry
• Holdups
• Heating
• Disengagement
• Catalyst
General
• Reactions on/off
Power Law
• Number of reactions
• Reaction concentration basis
• Reaction phase
• Activation energy
• Pre-exponential factor
• Stoichiometric coefficient for each component
• Temperature exponent
For more information on the equations, see the start of this chapter.
Langmuir-Hinshelwood-Hougen-Watson
• Number of reactions
• Number of reaction terms
• Reaction concentration basis
• Reaction phase
• Activation energy
• Pre-exponential factor
• Stoichiometric coefficient for each component
• Temperature exponent
• Adsorption constant
• Driving force 1 constant
• Driving force 2 constant
• Adsorption power
• Adsorption exponent
• Driving force 1 exponent
• Driving force 2 exponent
For more information on the equations, see the start of this chapter
User Kinetics
• Stoichiometric coefficient for each component
• User constants
• User reaction subroutine name
Method Description
Vapor Only vapor
Liquid Only liquid
Liquid with Breakthrough Liquid until level reaches low limit, then
homogeneous
Homogeneous Fully mixed
DIERS Bubbly From DIERS manual
DIERS Bubbly-Morris As above with corrections from Morris
DIERS Churn-Turbulent From DIERS manual
Tip The key to the above methods is to predict when the liquid
swells to the vent line and when it does, to calculate the quality of
the fluid venting.
Disengagement Equations
The equations for disengagement are taken from Fisher, 1992.
Where:
K = A constant determined by fitting to steady-
state conditions
∆P = Pressure drop
Rho = Mass density at outlet conditions
Fv_out = Outlet volumetric flow rate
If the steady-state simulation predicts a pressure rise, the pressure
drop equation is not used. Instead the pressure difference will be
specified.
When an isentropic turbine is modeled, the discharge pressure is
specified and the power output is calculated.
The discharge flow rate is determined by the material balance.
Where liquid knockout streams are modeled, their flow rates are
also determined by the material balance.
Using Performance When using performance curves for MCompr in a dynamic
Curves for MCompr simulation, note the following points:
head
Head ratio =
operating head
hydraulic efficiency
Efficiency ratio =
operating hydraulic efficiency
This table shows the values of the ratios used to give typical
performance curves:
Point Flow ratio Head ratio Efficiency ratio
1 0.0 1.25 0.2
2 0.1667 1.24 0.4
3 0.3333 1.23 0.595
4 0.5 1.21 0.765
5 0.6667 1.17 0.893
6 0.8333 1.11 0.975
7 1.0 1.0 1.0
8 1.1667 0.82 0.93
9 1.3333 0.49 0.6
10 1.4167 0.0 0.2
1.2
Operating
point
1
0.8
Head ratio
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
Flow ratio
0.8
Efficiency ratio
0.6
0.4
0.2
0
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6
Flow ratio
If the flow through the pump and pressure change across the pump
vary during the dynamic simulation, it is recommended that you
use performance curves.
Model Equations
The flow rate through a control valve under turbulent flow
conditions is given by:
FM w = N 6 F p C v ∆Peff ρ eff (1)
Flow Coefficient
The flow coefficient is a constant of proportionality with units of
L2 that relates that pressure drop driving force across the valve to
the mass flow rate through it. The actual flow coefficient depends
on the valve opening position and is related to the maximum flow
coefficient (at 100% opening) by the valve characteristic equation.
Alternatively, the user may specify the valve C v versus percentage
opening position in a tabular form to which a cubic spline is fitted.
Effective Density
The effective density at the valve inlet is determined from the fluid
phase type For liquid (incompressible) flow the effective density is
simply the actual fluid density. For vapor (compressible) flow the
effective density depends on an expansion factor which accounts
for the change in density of the fluid as it passes from the valve
inlet to the vena contracta and for the change in area of the vena
contracta as the pressure drop is varied.
Multiphase Flow
The valve sizing equations presented by the ISA are based on
single phase liquid or vapor flows. The ISA recommend avoiding
using control valves for mixed (two- and three-) phase service if at
all possible, hence formal sizing equations are not presented for
these cases.
In Valve multi-phase situations are handled in an approximate
manner. The effective choked pressure drop in mixed phase
service is taken as the smaller of the two choked pressure drops
calculated on the basis of all liquid and all vapor phase flows. The
user may also choose the calculation basis of the choked pressure
drop. The effective density is calculated using a phase fraction-
weighted mean value of the effective liquid and vapor densities.
Cavitation Index
The cavitation index K c is a ratio of the total pressure drop across
the valve to the pressure drop from the valve inlet to the vapor
pressure, the latter representing the pressure drop across the valve
Simple Valve Model The shortcut (simple) valve model uses the inlet conditions and a
simple sizing equation to relate the pressure drop cross the valve to
the flow rate through it:
1 2
F 2M w = C v ρ m ∆Pact (2)
2
where:
Pos MAX
Cv = Cv (3)
100
The outlet conditions are then computed from an isenthalpic flash
at the outlet pressure.
The maximum flow coefficient in (3) is fitted to match the steady
state Aspen Plus SM results of its shortcut method using a default
opening position of 50%.
These equations are used to allow Aspen Plus shortcut
(PressureChanger) valves to be exported and configured for use in
Aspen Dynamics pressure driven simulations.
Valve Dynamics Valve includes equations that model the dynamics of a valve
actuator. These equations relate the exposed valve position (i.e. the
valve position that can be manipulated, for example, by a
controller) to the actual valve position (i.e. the valve position in the
model equations) in terms of the dynamic valve position, a time
constant and the dynamic valve position first- and second-order
time derivatives.
The valve dynamics model the effect of:
1. The valve action, direct or reverse
2. First- or second-order or dynamics
3. Velocity limiting
4. Hysteresis effects caused by valve dead band
Actual Valve
Stem
Position
Equations for Valve When the first-order dynamics option is used, the valve stem
Dynamics response is modeled using the following first-order differential
equation:
d (Pos dyn )
Pos = Pos dyn + Tau1.
dt
Where:
Tau1 = First-order time constant
For second-order dynamics, the following second-order differential
equation is used:
d 2 (Pos dyn ) d (Pos dyn )
Pos = Tau 2 . 2
+ 2.Gamma.Tau. + Pos dyn
dt dt
StemMass
Tau =
Hookes
Where:
Pos = Exposed valve position that can be
manipulated (e.g. by a controller)
Continued overleaf…
StrokeTimeDo Time for full scale s 0.01 -> 1E8 1.0
down stroke
StrokeTimeUp Time for full scale s 0.01 -> 1E8 1.0
up stroke
Dband Dead band % 0 -> 100 0.0
If your process models are written to work in time units other than
hours, you will need to change the control model time units. See
Time Units in Control Models, earlier in this chapter.
Valve Closing
position
Opening
0 100
Input
signal
DBand
Note.
• In some circumstances, such as when creating a pressure relief
simulation, you may need to use the Setup button on the
Configure form to initialize an instance of Pipe. The Setup
procedures for Pipe are based on the assumption of forward
flow, so you should avoid using a negative pressure gradient
for the initial starting point.
What Does the The Initialize buttons performs the following tasks:
Initialize Button Do? • Changes the pressure specifications by fixing the inlet and
freeing the outlet pressure
• Initializes all the flashes from the upstream unit
• Changes the initialize parameter to stage1
• Performs an initialization run
• Changes the initialize parameter to stage2
Option Description
Stage1 Assumes all flashes and flows are at inlet condition
Stage2 Assumes all flows are equal
No Normal
Displaying Results To display results for the model, open the required results form.
for the Pipe Model These forms are available for displaying results:
Equations for the The following equations within the Integrate method Pipe model
Pipe Model are for single-phase or homogeneous flow. The Pipe model uses
similar equations for multi-phase, separated flows, with
appropriate adaptations.
Mass Balance
∂M i ∂Fi
+ =0
∂t ∂z
M = Moles
i = Component index
t = Time
F = Molar flow rate
z = Axial position
Energy Balance
∂E ∂
+ ( FH ) − Q = 0
∂t ∂z
t = Time
F = Molar flow rate
z = Axial position
E = Energy
H = Specific enthalpy
Q = Duty
Momentum Balance
Dynamic
∂ ∂ ∂
( Mom) + ( Fm u ) + ( Ap) + FRICT + Aρg sin θ = 0
∂t ∂z ∂z
Instantaneous
∂ ∂
( Fmu ) + ( Ap) + FRICT + Aρg sin θ = 0
∂z ∂z
Comparator
Input1
Output_
Input2
Input_ Output_
If your process models are written to work in time units other than
hours, you will need to change the control model time units. See
Time Units in Control Models, earlier in this chapter.
Discretize
Input_ Output_
If your process models are written to work in time units other than
hours, you will need to change the control model time units. See
Time Units in Control Models, earlier in this chapter.
FeedForward
HiLoSelect
Input 1
Output
Input 2
IAE
Input
IAE
SP
Where:
e = Deviation of the variable from the desired set point
t = Time
IAE = Integral absolute error value
T = Current time
SP for IAE
SP specifies the required value of the input variable. The IAE
element calculates its performance criterion as the integral of the
absolute difference between this required value and the input
value. SP may be specified within the IAE block, or supplied
through an input signal to the block. If SP is supplied through an
input signal, make sure you change its Spec from Fixed to Free on
the block Configure table.
ISE
Input_
ISE
SP
ISE calculates the integral of the squared error between its input
and a set point value.
The model provides a time-integral performance criterion. You can
use ISE to measure how successful a control system has been in
keeping a process variable at its set point (SP) over the entire
dynamic response of a controlled process. SP can be specified
within the block or be an input from another block.
ISE Equation The equation used in the model is the standard integral of the
absolute value of the error form:
T
ISE = ∫
0
e(t ) 2 . dt
Where:
e = Deviation of the variable from the desired set point
t = Time
ISE = Integral absolute error value
T = Current time
SP for ISE
SP specifies the required value of the input variable. The ISE
element calculates its performance criterion as the integral of the
squared difference between this required value and the input value.
SP may be specified within the ISE block, or supplied through an
input signal to the block. If SP is supplied through an input signal,
make sure you change its Spec from Fixed to Free on the block
Configure table.
Lag_1
Input_ Output_
Lag_1 models a first order lag between the input and output
signals.
Lag_1 Equations The Laplace domain transfer function for the first order lag is:
Gain
g (s ) =
Tau.s + 1
Where:
Gain = Steady-state gain (ultimate change in output divided by change in input)
Tau = Time constant in minutes
If your process models are written to work in time units other than
hours, you will need to change the control model time units. See
Time Units in Control Models, earlier in this chapter.
Lead_lag
Input_ Output_
If your process models are written to work in time units other than
hours, you will need to change the control model time units. See
Time Units in Control Models, earlier in this chapter.
Noise
Input_ Output_
Noise generates a Gaussian noise signal. You can use the Noise
model in two ways:
• To add noise to the input signal
• To generate a noisy signal
Noise can be used for testing how well your control system rejects
noise.
Noise Equations The Noise model generates a random value with amplitude that
follows a Gaussian probability distribution. The mean of the noise
is zero, and the magnitude of the noise is determined by the
specified standard deviation. The output of this model is the sum
of the input to the model plus the noise.
The model approximates Gaussian white noise, which is often a
good representation of noise found in process measurement and
control systems.
Probability
−3σ −2σ −σ 0 σ 2σ 3σ
Amplitude
σ = Standard deviation
Configuring Noise To add noise to an input signal, connect to both the Input_ and
Output_ connections. To generate a noisy signal, connect Output_
only, and Fix Input_ to the required mean value of the noisy signal.
Noise has the following configuration parameters:
Property Description Units Valid values Default Value
StdDev Standard deviation – -1E9->1E9 0.0
SPRemote
PV OP
Configuring PID
Use the Configure form to enter parameters for PID.
The form is divided into four tabs for configuring different aspects
of the controller. Each of these is explained below. You will need
to change values on the Tuning and Ranges tab, but the default
values on the Filtering and Other tab are suitable for most
applications.
To help you configure the controller, ensure that you have
connected the Process Variable (PV) and output (OP) connections,
and then use the Initialize Values button on the Configure form.
Name Initialized to
Operator set point Measured Variable
Bias Manipulated Variable
PV range minimum If Measured Variable > 0 0
If Measured Variable < 0 2 x Measured Variable
PV range maximum If Measured Variable > 0 2 x Measured variable
If Measured Variable < 0 0
Output range If Manipulated Variable > 0 0
maximum If Manipulated Variable < 0 2 x Manipulated
Variable
Output range If Manipulated Variable > 0 2 x Manipulated
minimum Variable
If Manipulated Variable < 0 0
If your process models are written to work in time units other than
hours, you will need to change the control model time units. See
Time Units in Control Models, earlier in this chapter.
Manipulated variable
Proportional error
Series Algorithm
The equation used to determine the controller output (OP) is:
1 d (E D )
OP = Bias + Gain.E P +
IntegralTi me ∫ E I .dt .1 + DerivTime
dt
Parallel Algorithm
The equation used to determine the controller output (OP) is:
1 d (E D )
OP = Bias + Gain.E P +
IntegralTime ∫ E I .dt + DerivTime
dt
Note. All values of the variables used in these equations are the
scaled values based on the process variable range minimum and
maximum or the output range minimum and maximum.
PID Controller The PID model includes two controller faceplates that you can use
Faceplates to interact with the controller during a running simulation:
• Full faceplate
• Compact faceplate
The full faceplate is similar to that found on real PID controllers. It
includes three horizontal bars which show the set point (SP),
process variable (PV), and output (OP) as a percentage of range.
To the right are the actual numerical values of SP, PV, and OP in
process units.
Using the ISE and IAE The ISE model and IAE model can be used to give a measure of
Models with a PID how successful a control system is at keeping a process variable
Controller (PV) at its set point (SP). When used with a PID controller, it is
convenient to link the input of the ISE and IAE models directly to
the setpoint and process variable of the controller. For this reason,
the PID setpoint (SP) and process variable (PV) are defined as
output control connections, and you can use control streams to
connect the PID controller and ISE or IAE block, as in the
following table:
Connect this PID controller: To this connection of
the ISE or IAE block:
SP SP
PV Input_
SPRemote
PV OP
Configuring PIDIncr
Use the Configure form to enter parameters for PID.
The form is divided into four tabs for configuring different aspects
of the controller. Each of these is explained below. You will need
to change values on the Tuning and Ranges tab, but the default
values on the Filtering and Other tab are suitable for most
applications.
To help you configure the controller, ensure that you have
connected the process variable (PV) and output (OP) connections,
and then use the Initialize Values button on the Configure form.
Parameter Initialized to
Set point Measured Variable
Initial Output Manipulated Variable
PV range minimum If Measured Variable > 0 0
If Measured Variable < 0 2 x Measured Variable
PV range maximum If Measured Variable > 0 2 x Measured variable
If Measured Variable < 0 0
Output range If Manipulated Variable > 0 0
maximum If Manipulated Variable < 0 2 x Manipulated
Variable
Output range If Manipulated Variable > 0 2 x Manipulated
minimum Variable
If Manipulated Variable < 0 0
If your process models are written to work in time units other than
hours, you will need to change the control model time units. See
Time Units in Control Models, earlier in this chapter.
Filter time constant (PVFilter) has units of time. The default value
is 0.0333 minutes, which is 2 seconds.
Ideal Algorithm
The equation used to determine the controller output (OP) is:
1 d (E D )
OP = Bias + Gain ∗ E P +
IntegralTi me ∫ E I .dt + DerivTime
dt
Series Algorithm
The equation used to determine the controller output (OP) is:
1 d (E D )
OP = Bias + Gain.E P +
IntegralTi me ∫ E I .dt .1 + DerivTime
dt
Parallel Algorithm
The equation used to determine the controller output (OP) is:
1 d (E D )
OP = Bias + Gain.E P +
IntegralTime ∫ E I .dt + DerivTime
dt
Note. All values of the variables used in these equations are the
scaled values based on the process variable range minimum and
maximum or the output range minimum and maximum.
Manipulated variable
Proportional error
PIDIncr Controller The PID model includes two controller faceplates that you can use
Faceplates to interact with the controller during a running simulation:
• Full faceplate
• Compact faceplate
The full faceplate is similar to that found on real PID controllers. It
includes three horizontal bars which show the set point (SP),
process variable (PV), and output (OP) as a percentage of range.
To the right are the actual numerical values of SP, PV, and OP in
process units.
Pressing the fifth button opens the Configure form so that you can
easily change configuration parameters.
The sixth button opens the plot for the controller which shows
values of SP, PV and OP either in process units, or as a percentage
of the range, versus time.
The seventh button opens the Tune form which you can use to
automatically determine tuning parameters for the controller.
To save space the compact faceplate only includes the three
buttons required to change controller mode, and does not include
bars to represent values. Otherwise the behavior is the same as for
the main faceplate.
Using Automatic Before using automatic tuning ensure that the controller Process
Controller Tuning Variable and Output connections are in place. Also ensure that you
have entered the controller input and output ranges and the Set
Point. If you use the Closed Loop ATV test method you must also
have entered the correct controller action.
Next choose your required test method and test settings.
Start a dynamic simulation and open the controller faceplate. Wait
until the controller PV and OP are stable. This is easiest to see if
you plot the values using the plot button. Then press the Start test
button to start the test.
When the test has run press the Finish Test button. For an Open
loop test this is when the PV has stopped changing significantly.
For a closed loop ATV test this is when each cycle looks the same
as the previous one.
The test results will then be calculated and displayed
automatically.
For step tests, in some cases the automatic tuner may find it
difficult to estimate the Time constant and Dead time from the test
results. This can be true if you are using the Discretize model to
simulate the effect of an analyzer in the control loop. For this
reason we recommend you visually check that the values shown
look reasonable. The diagram below shows how the values can be
determined from the step response curve.
54
62.
5
53
52 60
Controller Output
Process Variable
Set Point
51 57.
5
50
55
49
48 52.
5 Dead time Time constant
47
50
46
45
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7
Time Hours
If you would like to change the values used for the calculation of
the tuning parameters type in the values you want to use in place
of those determined by the automatic tuner.
If for any reason you wish to stop a test without calculating the
results press the Cancel test button to abort the test.
To calculate tuning parameters select the Tuning parameters tab.
Select the tuning parameter options and press the Calculate button
to calculate the results. If you wish to update the controller with
the calculated values click the Update controller button.
Tip:
• If your control loop has very fast dynamics, or if you are using
Closed Loop ATV testing on a loop with a very small dead
time, you may need to reduce the Simulation Communication
Interval to ensure that the fast changes in the PV are
adequately captured. You can do this on the Run Options
dialog in the section Time control. You may also need to
reduce the integrator minimum step size. Do this on the
Integrator tab of the Solver Options dialog.
Output_
Upper
bound
Amplitude
Datum
Amplitude
Lower
bound
Time
Upper
bound
Amplitude
Datum
Amplitude
Lower
bound
Time
Ratio
Input1
Output_
Input2
Input_ Output_
SplitRange Equations The following illustration shows the relationship between the input
and output signals:
SteamPtoT
P PT T
Sum evaluates its output as the sum of the two input signals.
Sum Equations The equation used in the Sum model is:
Output _ = Input 1 + Input 2
Transform
Input_ Output_
Superseded Models
Valve_dyn Input_
Output_
StemMass
Tau =
Hookes
Where:
Gamma = Damping coefficient
Tau = Natural period
StemMass = Moving mass of valve stem
Hookes = Hooke’s law constant for spring
Continued overleaf…
StrokeTimeDo Time for full scale s 0.01 -> 1E8 1.0
down stroke
StrokeTimeUp Time for full scale s 0.01 -> 1E8 1.0
If your process models are written to work in time units other than
hours, you will need to change the control model time units. See
Time Units in Control Models, earlier in this chapter.
Valve Closing
position
Opening
0 100
Input
signal
DBand
ε =
1
1−
(1 − x& )
K
S (
1 − x& 1 − 0.05 ⋅ We 0.27 ⋅ Re )
0.05
S=
x& (1 − ε ) ρ fl
(1 − x& ) ε ρg
We =
& x&
m 2 2
d (ρ fl − ρg )
ρg σ *
ρg
& (1 − x& )d
m
Re =
ηfl
0 .7
2 ηfl , 2
f * Exp = 0.0049 x& 2 (1 − x& )
ηg , 2
ΓExp = 1 − σ 2,3 0.25
Subscripts
Exp = Expansion
Con = Contraction
1 = Inlet with fully formed flow
2 = Point of size change
3 = Outlet with fully formed flow
fl = Liquid
g = Gas
“J. Schmidt & L. Friedel – Two-phase Flow Pressure Changes
Across Sudden Expansion in Duct Areas”.
Chem. Eng. Comm; Vols 141-142; pp175-190; 1996.
m3
& 2
1
−
σ 3 .1 *
+ f Conρeff ,3
x&
−
(1 − x& )
2
(
1 − σ 3,1 )
2
( 3 )ρfl
ρeff ,3 ρeff ,1 ε 3 ρg , 3 − ε
1
∆PCon =
1
1 − ΓCon − 1
σ 3,1
(1 − x& ) (1 − x& )
2 2 2
1 x& 2 ε x&SK
+ ρ 1− ε
K
= +
fl
−
ρ eff ρ ε ρ (1 − ε )
g fl 1− ε
K
ρgε (1 − ε )ρfl
ε =
1
1−
(1 − x& )
K
0.50
S 1 − x& (1 + 0.18 ⋅ We 0.12 ⋅ Re )
S=
x& (1 − ε ) ρ fl
(1 − x& ) ε ρg
We =
& x&
m 2 2
d (ρ
fl − ρg )
ρg σ *
ρg
& (1 − x& )d
m
Re =
ηfl
0.8
ηfl
f * Con = 5.2 ⋅ 10 −3 x& 0.1 (1 − x& ) σ 3,1
ηg ,3
ΓCon = 0.77σ 3,1(1 − σ 3,1 0.306 )
Subscripts
Exp = Expansion
Con = Contraction
1 = Inlet with fully formed flow
2 = Point of size change
3 = Outlet with fully formed flow
fl = Liquid
g = Gas
“J. Schmidt & L. Friedel – Two-phase Flow Pressure Changes
Across Sudden Expansion in Duct Areas”.
Chem. Eng. Comm; Vols 141-142; pp175-190; 1996.
∆Prec = ρ
(ν 2
2
−ν1
2
)
2
∆Pnonrec = ρ ν 12 K
2
A
K = 1 − 1
A2
Where:
ν = Velocity (m/s)
A = 2
Cross-sectional area (m )
ρ = Density (kg/m )
3
K = Loss coefficient
= Upstream
1
= Downstream
2
∆Prec = ρ
(ν 2
2
−ν1
2
)
2
∆Pnonrec = ρ ν 22K
2
A
K = 0.45 ( 1 − 2 )
A1
Where:
ν = Velocity (m/s)
A = 2
Cross-sectional area (m )
ρ = Density (kg/m )
3
K = Loss coefficient
= Upstream
1
= Downstream
2
Specifying the Component List and Valid Phases for the PSV
Model
To specify the component list and valid phases, open the Configure
form. On the Type tab, you can specify the following:
For this option You can specify
Component List From user defined list
Valid phases Liquid
Vapor
Vapor-Liquid
Vapor-Liquid-Liquid
100%
Valve
Opening %
Reset
Lift
Primary Lift
Blowdown Overpressure
Reset Set
Pressure Pressure Inlet Pressure
Note Only one of the Full lift, closing and Reset pressures needs
to be specified. The other is calculated to maintain a slope parallel
to the opening characteristics.
By careful selection of these parameters, you can configure many
types of relief valve, for example:
• Pop open, pop close
• Gradual open, gradual close
Incompressible Flow
For all flow methods the following relationship is applied:
V
Flow MASS = C d .Flux
100
Where:
Flux = Mass flux (kg/m2/s)
V = Valve opening (%)
Cd = Discharge coefficient
Simple (incompressible fluid)
Flux = Cv 2.∆P.ρ 0
Cv = Sizing coefficient
∆P = Pressure drop (bar)
ρ0 = Inlet density (kg/m3)
v P β
− 1 = α N − 1
vN P
Where:
N = A, B or C and
IF PA 〈 2.Pexit
PB = 0.5(PA + Pexit )
PC = Pexit
ELSE
PB = 0.75 PA
PC = 0.5 PA
ENDIF
Solve for α and β
2 Pa
G2 = −
v a2 ∫
P0
v .dP
MAX
Pa
η=
Po
Gc2 =
Pa v fa
− C
(
(Vga − V fa ) 1 / n −1 / k )
X0 n
Pg
(S g 0 − S f 0 )
(1 − α 0 ) (1 − η ) + k
ηc
k −1
k
=
α0 c
(k − 1)
1 k
+
2ϕα 2
a (k − 1)
Where:
C pf
( [ (1 − X )
]+ X )
C pg
n=
C
( [ (1 − X ) pf ] + X )
C pg k
1 1
−
1 n k
ϕ = − CPg
n (S g 0 − S f 0 )
1 (1 − X ) v f
=1+
α X vg
Subscripts
0 = Stagnation
c = Critical
a = Throat
g = Vapor
f = Liquid
Schmidt-Westphal Equations
The literature reference for these equations is:
J. Schmidt and F. Westphal (1997). “Praxisbezogens Vorgehen bei
der Auslegung von Sicherheitsventilen und deren
Abblaseleitungen für die Durchströmung mit Dampf/Flüssigkeits-
Gemischen” In: Chemie Ingenieur Technik, Vol. 69, pp. 776–1091.
P ' (T0 )
ηs =
P0
2 ⋅ω
η s,MIN =
1+ 2 ⋅ω
2
v g0 − v f 0 η s ⋅ C Pf 0 ⋅ P0 ⋅ T0 v g 0 − v f 0
ω = x& 0 ⋅ + ⋅
v0 v0 ∆hv 0
ELSEIF : η s > 1
k
2 k −1
ηc =
K + 1
k k
2 k +1
ψ = η − η k
k − 1
ELSE
ω 2 − 2ω + 1 2 η 3
0= η c − 2(ω − 1)η c + ω ⋅ η s ⋅ ln c + ω ⋅ η s − 1
2 ⋅ ηs ⋅ ω ηs 2
η
2 ⋅ (1 − η s ) + 2 ⋅ ω ⋅ η s ⋅ ln s − (ω − 1)(η s − η )
η
ψ =
η
2 ⋅ ω ⋅ s − 1 + 1
η
ENDIF
WHERE
Pexit
ηa =
P0
IF : η a ≤ η c
η = ηc
ELSE
η = ηa
ENDIF
α = α g ⋅ ε 0 + α f ⋅ (1 − ε 0 )
FLUX = α ⋅ ψ
Where:
v = Specific volume (m3/kg)
η = Pressure ratio
P = Pressure (bar)
x& = Mass vapor fraction
Cp = Specific heat at constant pressure (kJ/kg/K)
α = Flow factor
Diener Equations
The literature reference for these equations is:
R.Diener and L.Friedel (1997). “Procedure to calculate the mass
flowrate through certain types of valves, rupture discs, orifices,
nozzles and other flow contractions”, (internal document).
2(1 − x& i ) 2
ε i = 1−
v g ,i
1 − 2 x& i + 1 + 4 x& i (1 − x& i ) − 1
v f ,i
x& i ⋅ v g ,i (1 − x& i ) ⋅ v f ,i
2
v 2 ph,i = +
εi 1− ε i
v hom,i = x& i ⋅ v g ,i + (1 − x& i ) ⋅ v f ,i
(v 2 ph,a − σ *2 ⋅ v 2 ph,in )
FLUX 2(σ * ⋅ v 2 ph,a − σ *2 ⋅ v 2 ph,exit
Pa = Pexit − − ζ form ⋅ v hom, a
1 − (1 − σ )(1 − σ * )
*0.25
2
Leung-Grolmes Equations
The literature reference for these equations is:
J. C. Leung and M. A. Grolmes (1988). “A Generalized correlation
for flashing choked flow of initially subcooled liquid”.
Journal of AIChE, Vol. 34 (4). Pp.688-691.
Ps P P
ηs = ; ηc = c ; η = ;
Po Po Po
2
C pf To Ps υ fgo
ω=
υ fo h fgo
2ω − 1
IF ηs <
2ω
Gc = 2 ρ fo
2
( Po − Ps )
ηc = ηs
ELSE
Pρ
Gc = ηc 0 fo
2
ωη s
2ω 1 2ω − 1
ηc = ηs 1 − 1 −
2ω − 1 η s 2ω
flux = Gc
flux = G
ENDIF
Where:
η = Pressure ratio
ω = intermediate function
Cp = specific heat at constant pressure (kj/kg k)
ρ = Pressure (N/m2)
T = Temperature (k)
υ = Specific volume (m3/kg)
h = Specific enthalpy (kj/kg)
Subscripts
s = Saturation
0 = Stagnation
c = Critical
f = Liquid
g = Vapor
fg = Difference between vapor and liquid phase property
Where:
G = Mass flux at discharge
Cp
k = Isentropic expansion exponent ( )
Cv
ρ0 = Mass density at stagnation conditions
P0 = Stagnation pressure
Pc = Critical pressure
P
η= = Pressure ratio
P0
For choked flow:
k
2 k −1
η = ηc =
k + 1
Where:
Pc
ηc = = Critical pressure ratio
P0
This equation can be used for non-ideal gases provided the
assumptions of constant compressibility and constant k are not
seriously violated.
The literature references for these equations are:
1. “Sizing, Selection and Installation of Pressure-Relieving
Devices in Refineries. Part I- Sizing and Selection” API
Recommended Practice 520, 5th edition, (1990), pp. 30-34.
2. “Guidelines for Pressure Relief and Effluent Handling
Systems”, Center for Chemical Process Safety of the
American Institute of Chemical Engineers(1998), pp. 253-261
3. McCabe W., L., Smith J. C., Harriott, P., “Unit operations of
Chemical Engineering”, (1985), pp. 112-118.
Where:
Visc is viscosity in cP
T is Temperature in K
A, B are specified by the user
∆Pdesign
K= * ρ design
Flow 2 design
Where:
∆P = Actual pressure drop
ρ = Fluid density
K = Flow coefficient
Flow = Volumetric flow rate
∆Pdesign = Design pressure drop
Flow design = Design volumetric flow rate
ρ design = Design mass density
Procedure pChairX
Compute "chair" function and "excess" function:
call (Y1, Y2) = pChairX (YR, X, YMAX, XMINF,
XMAX, DXFL, DXFU, TYPE);
Output Description
Y1 Actual value of Y
Y2 Excess (overflow)
(ramp)
If X > XLOW: Y1 = YR (normal)
Where:
XLOW = XMIN + DXL
Usage
This procedure has 2 output arguments, typically the availability
Av, and an vent flowrate Fl_over or Fv_vent.
// liquid case
Liquid Case call (Av, Fl_over) = pchairx (1.0, Vl, F_max,
V_heel, V, del_Vl, del_Vu, 1);
F_max is 1e6
V_heel is -0.001
del_Vl is -0.005
del_Vu is -0.005
This gives XMIN = V*0.001
// vapor case
Vapor Case call (Av, Fv_vent) = pchairx (1.0, P, F_max,
P_min, P_max, del_Pl, del_Pu, 0);
What is the This is a variable that vary from 0 to 1. The normal case is that it is
Availability (Av)? equal to 1. For a flow driven simulation, instead of having the user
specifying directly the flowrate F, we specify the "required"
flowrate FR. The two variables are related by:
F = Av*FR;
If Av is equal to 1, F is simply equal to FR. For feed streams, Av is
actually fixed to 1. For the case of a liquid stream from a vessel
(i.e. a Flash2 block), Av is related to the liquid level. As the liquid
level becomes too low, Av is ramped down to 0. This means that
regardless of the value of FR, the actual flowrate F is zero, which
prevents trying to remove more liquid from the vessel that is
practically empty.
For the pressure driven case, the principle is similar. The
flow/pressure drop relation is based not on F but on FR. As Av
goes down to zero, the actual pressure difference to drive the flow
has then to increase dramatically, which has the effect of stopping
the flow, thus preventing the removal of further liquid from the
empty vessel.
Lead_lag 9-10
A Multiply 9-12
Noise 9-12
Algorithms 9-24, 9-25, 9-36, 9-37
PID 9-14, 9-28
PRBS 9-45
B
Ratio 9-47
Bursting disk model 10-2 Scale 9-48
SplitRange 9-49
C SteamPtoT 9-51
Catalyst deactivation 6-26, 6-32 Sum 9-52
Catalyst heat transfer 6-36 time units 9-2
Choking 7-13 Transform 9-52
Closed loop tuning 9-26 Valve_dyn 9-53
Column models Cooling options 6-33, 6-34
about 5-1
Distl 5-2 D
Extract 5-3 Dead_time model 9-3
RadFrac 5-4, 5-24 Decanter model 3-15
Comparator model 9-2 DeltaP model 10-38
Compartmentation 5-2 Discretize model 9-3
Compr model 7-2 Disengagement 6-25, 6-32
Compressor types 7-2 Distl model 5-2
Continuous stirred tank reactor (CSTR) 6-17 Draining of trays 5-44
Contraction equations (pipe) 10-8 Dupl model 8-2
Control models 9-1
Comparator 9-2 E
Dead_time 9-3
Edmister approach 5-2
Discretize 9-3
Entrainment in separators
Feedforward 9-4
Flash2 3-7
HiLoSelect 9-6
Flash3 3-14
IAE 9-7
Equipment heat transfer 7-23
ISE 9-8
Decanter model 3-19
Lag_1 9-9
Flash2 model 3-7
W
WorkStream model 1-6
Z
Ziegler-Nichols technique 9-26