AspenHYSYSDynModeling2006 Ref
AspenHYSYSDynModeling2006 Ref
AspenHYSYSDynModeling2006 Ref
Dynamic Modeling
Controller, Aspen STX Upgrade to TASC, Aspen SULSIM, Aspen Supply Chain Analytics - Demand Management, Aspen
Supply Chain Analytics - Plant Scheduling, Aspen Supply Chain Analytics - S&OP, Aspen Supply Chain Analytics - Supply
Planning, Aspen Supply Chain Connect, Aspen Supply Planner, Aspen Supply Planning - Strategic Analyzer, Aspen Tank
Management, Aspen TASC, Aspen Teams, Aspen TICP, Aspen Transition Manager, Aspen Utilities, Aspen Voice Fulfillment
Management, Aspen Watch, Aspen Water, Aspen Web Fulfillment Management, Aspen XPIMS, Aspen XPRESS Optimizer,
Aspen XPRESS Optimizer for DPO, Aspen Zyqad Development, Aspen Zyqad, aspenONE Product Trading & Blending, SLM,
SLM Commute, SLM Config Wizard, the aspen leaf logo and Plantelligence 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 AspenTechs 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.
Aspen Technology, Inc.
Ten Canal Park
Cambridge, MA 02141-2201
USA
Phone: (617) 949-1000
Fax: (617) 949-1030
Website http://www.aspentech.com
Technical Support
vi
Registered users can also subscribe to our Technical Support eBulletins. These e-Bulletins proactively alert you to important
technical support information such as:
Technical advisories
Product updates and releases
vi
vii
vii
viii
viii
Table of Contents
Technical Support..................................................... v
Online Technical Support Center ............................vi
Phone and E-mail ................................................ vii
1
1.2
1.3
1.4
1.5
1.6
HYSYS Dynamics.............................................1-49
1.7
References .....................................................1-70
2.2
2.3
2.4
Integrator ......................................................2-38
2.5
Event Scheduler..............................................2-45
2.6
2.7
3.2
3.3
3.4
3.5
3.6
References .....................................................3-58
Index.................................................................... I-1
ix
Dynamic Theory
1-1
1 Dynamic Theory
1.1 Introduction................................................................................... 3
1.2 General Concepts ........................................................................... 6
1.2.1 Mathematical Model Classification............................................... 6
1.3 Holdup Model ............................................................................... 13
1.3.1
1.3.2
1.3.3
1.3.4
1.3.5
1.3.6
1.3.7
1-1
1-2
Dynamic Theory
1.7 References....................................................................................70
1-2
Dynamic Theory
1-3
1.1 Introduction
Dynamic simulation can help you to better design, optimize, and
operate your chemical process or refining plant. Chemical plants
are never truly at steady state. Feed and environmental
disturbances, heat exchanger fouling, and catalytic degradation
continuously upset the conditions of a smooth running process.
The transient behaviour of the process system is best studied
using a dynamic simulation tool like HYSYS.
Contact your Aspentech
agent for more
information, or e-mail us
at
[email protected].
1-3
1-4
Introduction
Process optimization
Controller optimization
Safety evaluation
Transitions between operating conditions
Startup/Shutdown conditions
Dynamic Theory
1-5
1-5
1-6
General Concepts
Dynamic Theory
1-7
(1.1)
(1.2)
dY
------ + YY 2 = Kf ( u )
dt
(1.3)
Conservation Relationships
Material Balance
The conservation relationships are the basis of mathematical
modeling in HYSYS. The dynamic mass, component, and energy
balances that are derived in the following section are similar to
the steady state balances with the exception of the
accumulation term in the dynamic balance. It is the
1-7
1-8
General Concepts
(1.4)
d( o V )
----------------- = F i i F o o
dt
(1.5)
where:
Fi = flowrate of the feed entering the tank
i = density of the feed entering the tank
Fo = flowrate of the product exiting the tank
o = density of the product exiting the tank
V = volume of the fluid in the tank
1-8
Dynamic Theory
1-9
Component Balance
Component balances can be written as follows:
Rate of accumulation of component j = Flow of component
j into system - Flow of component j out of system +
Rate of formation of component j by reaction
(1.6)
(1.7)
where:
Cji = concentration of j in the inlet stream
Cjo = concentration of j in the outlet stream
Rj = reaction of rate of the generation of component j
1-9
1-10
General Concepts
Energy Balance
The energy balance is as follows:
Rate of accumulation of total energy = Flow of total energy
into system - Flow of total energy out of system +
Heat added to system across its boundary + Heat
generated by reaction - Work done by system on
surroundings
(1.8)
(1.9)
where:
u = internal energy (energy per unit mass)
k = kinetic energy (energy per unit mass)
= potential energy (energy per unit mass)
V = volume of the fluid
w = shaft work done by system (energy per time)
Po = vessel pressure
Pi = pressure of feed stream
Q = heat added across boundary
Qr = DH rxn r A , heat generated by reaction
1-10
Dynamic Theory
1-11
(1.10)
Solution Method
Implicit Euler Method
Yn+1 is analytically calculated to equal:
tn + 1
Yn + 1 = Yn +
f ( Y ) dt
(1.11)
tn
where:
dY
------ = f ( Y )
dt
(1.12)
1-11
1-12
General Concepts
The figure below shows the integration of f(Y) over time step, h,
using exact integration and the Implicit Euler approximation:
Figure 1.2
Exact Integration
f(Y)
f(Y)
fn + 1
Area
Area
f ( Y )dt
= ( fn + 1 ) h
fn
tn
tn+1
tn
tn+1
Integration Strategy
In HYSYS, dynamic calculations are performed at three different
frequencies:
Volume (pressure-flow)
Energy
Composition
1-12
Dynamic Theory
1-13
1-13
1-14
Holdup Model
1.3.2 Accumulation
The lagged response that is observed in any unit operation is
the result of the accumulation of material, energy, or
composition in the holdup. To predict how the holdup conditions
change over time, a recycle stream is added alongside the feed
streams. For example, the material accumulation in a holdup
can be calculated from:
Material accumulationnew = material flow into system +
material accumulationold (recycle stream) - material
flow out of system
(1.13)
1-14
Dynamic Theory
1-15
Figure 1.3
1-15
1-16
Holdup Model
Efficiencies
In HYSYS, you can indirectly specify the amount of mixing that
occurs between the feed phases and the existing holdup using
feed, recycle, and product efficiencies. These feed efficiency
parameters can be specified on the Efficiencies tab of the unit
operations Advance property view. Click the Advance button on
the Holdup page under the Dynamics tab to open the Advance
property view.
Figure 1.4
1-16
Dynamic Theory
1-17
1-17
1-18
Holdup Model
Nozzles
In HYSYS, you can specify the feed and product nozzle locations
and diameters. These nozzle placement parameters can be
specified in the unit operations Nozzles page under the Rating
tab.
Figure 1.6
1-18
Dynamic Theory
1-19
1-19
1-20
Holdup Model
Energy Balances
Heat is lost (or gained) from the holdup fluid through the wall
and insulation to the surroundings.
Figure 1.7
Dynamic Theory
1-21
(1.14)
k ins
= --------- A ( T wall T ins ) + h ( ins, surr ) A ( T ins T surr )
x ins
(1.15)
where:
A = the heat transfer area
x = the thickness
Cp = the heat capacity
T = the temperature
k = the thermal conductivity
h = the heat transfer coefficient
As shown, both the insulation and wall can store heat. The heat
loss term that is accounted for in the energy balance around the
holdup is h ( fluid, wall ) A ( T fluid T wall ) . If Tfluid is greater than Twall, the
heat is lost to the surroundings. If Tfluid is less than Twall, the heat
is gained from the surroundings.
1-21
1-22
Holdup Model
Simple
Detailed
Simple Model
The Simple model allows you to either specify the heat loss
directly or have the heat loss calculated from specified values:
Overall U value
Ambient Temperature
The heat transfer area, A, and the fluid temperature, Tf, are
calculated by HYSYS. The heat loss is calculated using:
Q = UA(T f T amb )
(1.16)
1-22
Dynamic Theory
1-23
Detailed Model
The Detailed model allows you to specify more detailed heat
transfer parameters. There are three radio buttons in the Heat
Loss Parameters group as described in the table below:
Radio Button
Description
Temperature
Profile
Conduction
fluid
wall
insulation
surroundings
Conductivity of material
Thickness of material
Heat capacity of material
Density of material
Equation (1.14) and (1.15) demonstrate how the
parameters are used by the heat loss model.
Convection
1-23
1-24
Holdup Model
Description
Vessel Level
Calculations
Vessel Pressure
Tray Hydraulics
1-24
Dynamic Theory
1-25
General Tab
This tab provides the same information as shown in the Holdup
page of the Dynamics tab. The accumulation, moles, and
volume of the holdup appear on this tab. The holdup pressure
also appears on this tab.
Figure 1.8
1-25
1-26
Holdup Model
Nozzles Tab
Refer to Section 1.6 HYSYS Dynamics for
more information.
1-26
Dynamic Theory
1-27
Efficiencies Tab
Refer to Section 1.6 HYSYS Dynamics for
more information.
Description
Feed Nozzle
Efficiency
1-27
1-28
Holdup Model
Nozzle Efficiency
Description
Product Nozzle
Efficiency
Recycle Efficiency
Properties Tab
Figure 1.11
Temperature
Pressure
Flow
Molar Fraction of the specific phase in the holdup
Enthalpy
Density
Molecular Weight
1-28
Dynamic Theory
1-29
Compositions Tab
The compositional molar fractions of each phase in the holdup
displays in the Compositions tab.
Figure 1.12
1-30
The holdup model brings the actual feed and product stream
properties to holdup conditions for the volume balance
equations using a rigorous or approximate flash. The pressure
flow solver returns information essential to the holdup model
calculations: the pressure of the holdup or the flow rates of
streams around the holdup.
Dynamic Theory
1-31
Volume Balance
For equipment with holdup, an underlying principle is that the
physical volume of the vessel, and thus, the volume of material
in the vessel at any time remains constant. Therefore, during
calculations in dynamics, the change in volume of the material
inside the vessel is zero:
V = Constant = f ( flow, h, P, T )
dV
------- = 0
dt
(1.17)
(1.18)
where:
V = volume of the vessel
t = time
flow = mass flowrate
h = holdup
P = vessel pressure
T = vessel temperature
(1.19)
1-31
1-32
In the volume balance equation, pressure and flow are the only
two variables to be solved in the matrix. All other values in the
equation are updated after the matrix solves. Each vessel
holdup contributes at least one volume balance equation to the
pressure-flow matrix. When sufficient pressure-flow
specifications are provided by you, any unknown(s) can be
solved whether it be a vessel pressure or one of its flowrates.
The volume balance equation allows you to observe pressure
effects in the vapour holdup due to disturbances in the feed.
Consider a separator whose feed flow suddenly increases.
Assume that the exit streams from the separator are specified
by you and are thus, constant. The vessel pressure would
increase for two reasons:
Resistance Equations
Flows exiting from a holdup are calculated from a volume
balance equation, specified by you, or calculated from a
resistance equation. In general, the resistance equation
calculates flowrates from the pressure differences of the
surrounding nodes. HYSYS contains unit operations such as
valves and heat exchangers which calculate flowrates using
resistance equations. The resistance equations are modeled
after turbulent flow equations and have the form:
Flow = k P
(1.20)
where:
Flow = mass flow rate
k = conductance, which is a constant representing the
reciprocal of resistance to flow
1-32
Dynamic Theory
1-33
(1.21)
where:
Flow = mass flowrate
Cv = conductance, which is a constant representing the
reciprocal of resistance to flow
P1 = upstream pressure
P2 = downstream pressure
For a more detailed
description on the
individual unit operations
and the resistance
equations associated with
them, see the appropriate
unit operation section in
the HYSYS Operations
Guide.
Resistance Term
Valve
Pump
Compressor/
Expander
Heater/Cooler/Heat
Exchanger/Air
Cooler/LNG
1-33
1-34
Unit Operation
Resistance Term
Dynamic Theory
1-35
In the flowsheet shown above, there are eight streams and one
vessel holdup. To fully define the pressure flow matrix, the
pressure and flow for each material stream and the pressure of
each holdup must be solved for. In short, two variables are
required for each material stream and one variable is required
for each holdup:
8 material streams x 2 + 1 vessel holdup x 1 = 17
pressure-flow variables
(1.22)
1-35
1-36
The pressure and flow of material streams are named Pstream name
and Fstream name, respectively. The pressure of the holdup is
named PH.
There are a number of equations which describe the relationship
between the pressures and flows in this network. They are as
shown below:
Pressure-Flow Equation
Description
# of Eqns
dP
---------H- = f ( P, T, holdup, flows )
dt
Separator
Volume Balance equation
PH = P2 = P3 = P5
Valves
Resistance equations
F 2 = K VLV100 P 1 P 2
F 4 = K VLV101 P 3 P 4
F 8 = K VLV102 P 7 P 8
This is the general form of the valve resistance equation.
The actual equations vary according to inlet stream
conditions.
General Flow relations
F1 = F2
F3 = F4
F7 = F8
Since the valves are usually not specified with holdup, this
relation is observed.
Mixer
General Pressure
relation
P5 = P6 = P7
F7 = F5 + F6
13
Dynamic Theory
1-37
Pressure-Flow Specification
Guidelines
The previous section outlined the number of pressure-flow
specifications that are required by the flowsheet in order for the
degrees of freedom to be satisfied. This section presents
possible PF specifications that can be made for the inlet and exit
streams of stand alone operations.
The purpose of this section is to demonstrate the range of
specifications that can be made for different unit operations in
HYSYS. It is hoped that this section provides insight as to what
should and should not be specified for each unit operation.
Valve
Rating information for the valve operation including the valve
type and Cv values can be input on the Sizing page in the
Ratings tab.
The dynamic valve can either be specified as having a set
pressure drop or a pressure flow relation. This option is set on
the Specs page of the Dynamics tab in the valve property view.
1-37
1-38
Heat Exchanger/Cooler/Heater
The dynamic heat exchanger can be specified as having a set
pressure drop or a Overall K-Value (pressure-flow) relation. This
option is set on the Specs page of the Dynamics tab in the heat
exchanger property view:
The heat exchange operations, like the valve, should use the
P-F spec option as much as possible to simulate actual
pressure flow relations in the plant.
1-38
Dynamic Theory
1-39
Separators
Rating information including the volume of the vessel, boot
capacity, and nozzle location can be input on the Sizing and
Nozzles pages in the Ratings tab.
A separator with no valves attached to the inlet and exit streams
requires at most one pressure specification. The other two
streams are specified with flows. A more realistic way to run the
separator is to attach valves to the inlet and exit streams of the
vessel. The boundary streams of the separator with valves
should be specified with pressure.
Condenser/Reboiler
Rating information for the condenser and reboiler including the
vessel volume, boot capacity, and nozzle location can be input
on the Sizing and Nozzles pages of the vessels Ratings tab.
It is highly recommended that the proper equipment be added
to the reflux stream (e.g., pumps, valve, etc.). In all cases, level
control for the condenser should be used to ensure a proper
liquid level.
The Partial Condenser has three exit streams:
1-39
1-40
Separation Columns
For all separation columns, the tray section parameters
including the tray diameter, weir length, weir height, and tray
spacing can be specified on the Sizing page in the Ratings tab of
the Main TS property view.
The basic Absorber column has two inlet and two exit streams.
When used alone, the absorber has four boundary streams and
therefore requires four pressure-flow specifications. A pressure
specification is always required for the liquid product stream
leaving the bottom of the column. A second pressure
specification should be added to the vapour product of the
column, with the two feed streams having flow specifications.
The basic Refluxed absorber column has a single inlet and two
or three exit streams, depending on the condenser
configuration. When used alone, the refluxed ratios has three or
four boundary streams (depending on the condenser) and
requires four or five pressure-flow specifications; generally two
pressure and three flow specifications. A pressure specification
is always required for the liquid product stream leaving the
bottom of the column.
The Reboiled Absorber column has a single inlet and two exit
streams. When used alone, the reboiled absorber has three
boundary streams and therefore requires three pressure-flow
specifications; one pressure and two flow specifications. A
1-40
Dynamic Theory
1-41
Compressor/Expander/Pump
Rating information for the dynamic compressor, expander, and
pump operations can be input on the Curves and Inertia pages
in the Ratings tab.
In general, two specifications should be selected in the
Dynamics Specifications group in the Specs page of the
Dynamics tab in order for these unit operations to fully solve.
You should be aware of specifications which causes
complications or singularity in the pressure flow matrix. Some
examples of such cases are:
Mixer/Tee
The dynamic mixer and tee operations are very similar. It is
recommended that the mixer be specified with the Equalize All
option in Dynamics mode. It is also recommended that the
dynamic tee not use the dynamic splits as specifications. These
1-41
1-42
options are set on the Specs page of the Dynamics tab in their
respective operation property views.
By specifying the dynamic mixer and tee as recommended, the
pressure of the surrounding streams of the unit operation are
equal if static head contributions are not considered. This is a
realistic situation since the pressures of the streams entering
and exiting a mixer or tee must be the same. With the
recommended specifications, flow to and from the tee is
determined by pressures and resistance through the flowsheet.
This is more realistic than using the split fractions which can
also cause complications with regard to flow reversal.
A number of streams can enter or exit a mixer or tee. For stand
alone operations, one stream must be specified with pressure.
The other inlet/exit streams are specified with flow.
1-42
Dynamic Theory
1-43
Steady State
The Steady State mode uses modular operations which are
combined with a non-sequential algorithm. Information is
processed as soon as it is supplied. The results of any
calculation are automatically propagated throughout the
flowsheet, both forwards and backwards.
Material, energy, and composition balances are considered at
the same time. Pressure, flow, temperature, and composition
specifications are considered equally. For example, a columns
overhead flow rate specification is replaced by a composition
specification in the condenser. The column can solve with either
specification.
1-43
1-44
Dynamics
Material, energy and composition balances in Dynamics mode
are not considered at the same time. Material or pressure-flow
balances are solved for at every time step. Energy and
composition balances are defaulted to solve less frequently.
Pressure and flow are calculated simultaneously in a pressureflow matrix. Energy and composition balances are solved in a
modular sequential fashion.
Because the pressure flow solver exclusively considers pressureflow balances in the network, P-F specifications are separate
from temperature and composition specifications. P-F
specifications are input using the one P-F specification per
flowsheet boundary stream rule. Temperature and composition
specifications should be input on every boundary feed stream
entering the flowsheet. Temperature and composition are then
calculated sequentially for each downstream unit operation and
material stream using the holdup model.
Unlike in Steady State mode, information is not processed
immediately after being input. The integrator should be run
after the addition of any unit operation to the flowsheet. Once
the integrator is run, stream conditions for the exit streams of
the added unit operation is calculated.
1-44
Dynamic Theory
1-45
The following table indicates some basic steps you can take to
set up a case in Steady State mode and then switch to
Dynamics mode.
Step
Description
Adding Unit
Operations
Equipment
Sizing
1-45
1-46
Step
Description
Adjusting
Column
Pressure
Logical
Operations
Adding Control
Operations
Enter HYSYS
Dynamic
Environment
1-46
Dynamic Theory
1-47
Step
Description
Adding
Pressure-Flow
Specifications
1-47
1-48
Step
Description
Trouble
Shooting
1-48
Dynamic Theory
1-49
1-49
1-50
HYSYS Dynamics
Dynamic Theory
1-51
This option is useful because you does not see any slowdown due to network problems. HYSYS is assumed,
however, that the license is there and you can make
changes resulting in the case being stopped when the
validation is complete. This option is recommended for
users with network based security.
Check out when in the relevant mode. This option
checks the licenses in a run-time usage mode
(essentially pay-per-use), so the HYSYS Dynamics
license is only checked out while the case is in dynamics
mode, and the HYSYS Dynamics license is returned while
the case is in steady state mode. This option is
recommended for users with token based security.
Dont check out. This option means that licenses are
not checked out. It is recommended for users that do not
have the licenses available. Input requiring the license is
then ignored by HYSYS.
There are two values that are common to each of the three radio
buttons found in the Detailed Heat Loss Model: the Overall Heat
Loss value and the Overall Heat Transfer Area.
1-51
1-52
HYSYS Dynamics
Description
Temperature
Profile
Conduction
Convection
The governing equations relating heat loss from the vessel and
the Detailed heat loss parameters shown here are discussed in
Section 1.3.4 - Heat Loss Model.
1-52
Dynamic Theory
1-53
1.6.2 Nozzles
The Nozzles page of the Rating tab contains information
regarding the elevation and diameter of the nozzles.
Figure 1.17
1-53
1-54
HYSYS Dynamics
The inlet and outlet connection point are all relative to ground
level, which is set by default at 0. All the elevation points can be
user-specified but default settings are supplied by HYSYS as a
starting point. You can adjust the elevations for each nozzle on
the Nozzle page of the Rating tab as desired. To turn off the
elevation display on the PFD, press SHIFT N and the name of
each stream reappears.
The following is a summary of the Nozzle Parameters section:
Nozzle Parameter
Description
Diameter
Elevation (Base)
Elevation (Ground)
Elevation (% of
Height)
Base Elevation
relative to Ground
Diameter
Height
Dynamic Theory
1-55
1-55
1-56
HYSYS Dynamics
Description
Traysection
Elevation Relative to
Ground
Holdup RG
VToAbove
LToBelow
1-56
Dynamic Theory
1-57
Object
Description
Elevation RH
Diameter
1-57
1-58
HYSYS Dynamics
Description
Actuator
Valve
Valve Modes
The valve mode defines the relationship between the desired
actuator position and current actuator position. The desired
actuator position can be set by a PID Controller or Spreadsheet
operation. A controllers output, OP, for instance, is exported to
the desired actuator position.
Depending on the valve mode, the current actuator position can
behave in one of the following three ways:
Instantaneous Mode
First Order Mode
Linear Mode
Instantaneous Mode
In this mode, the actuator moves instantaneously to the desired
actuator position defined by the controller. The equation
defining the relationship is:
Act% = Act
Desired
(1.23)
1-58
Dynamic Theory
1-59
(1.24)
Linear Mode
The actuator can be modeled to move to the desired actuator
position at a constant rate. The actuator moves according to the
following equation (if the desired actuator position is above the
current actuator position):
Act% = ( Actuator Linear Rate )t + Act o %
until Act% = Act Desired %
(1.25)
The linear rate can be specified in the Actuator Linear Rate cell.
Typical stroke times (closure rates) are as follows:
1-59
1-60
HYSYS Dynamics
Valve Stickiness
In reality, the valve does not respond instantaneously to
changes in the actuator. A first order lag can be modeled in the
response of the actual valve position to changes in the actuator
position. The behaviour of the valve percent opening as a
function of the actuator position is shown as follows:
d ( Valve% )
sticky ---------------------------- + Valve% = Act%+ Offset
dt
(1.26)
Fail Modes
Actuators usually have a fail-safe function. If there is a
disruption to the power source driving the valve, the actuator
places the valve in a safe position, either fail open or fail close.
Fail modes can be specified by selecting the corresponding radio
button in the Positions group. The valve can be modeled to fail
by selecting the Actuator has failed checkbox.
(1.27)
1-60
Dynamic Theory
1-61
Positions Group
Various valve position parameters for the actuator and the valve
appear in the Positions group:
Valve Position Parameter
Definition
Minimum
Maximum
Current
1-61
1-62
HYSYS Dynamics
Definition
Desired
Offset
1.6.4 Inertia
The inertia modeling parameters and the frictional or power loss
associated with the rotating equipment in the Pump,
Compressor, and Expander can be specified on the Inertia page
of the Rating tab for these unit operations.
Figure 1.21
The Size Inertia button option is only available in the Pump unit
operation.
Dynamic Theory
1-63
(1.28)
where:
F = molar flow rate
h2 = molar enthalpy of the exit stream
h1 = molar enthalpy of the inlet stream
MW = molecular weight
(1.29)
(1.30)
where:
cf = correction factor
Fm = mass flow rate
H = head (adiabatic or polytropic depending on selection)
= efficiency (adiabatic or polytropic depending on selection)
g = acceleration due to gravity
1-63
1-64
HYSYS Dynamics
(1.31)
where:
M = mass of the impeller and rotating shaft
R = radius of gyration
(1.32)
Ii GRi
where:
Ib = rotational inertia with respect to some base speed
I = rotational inertia of the ith component on the rotating
shaft (with respect to the speed at this gear ratio)
b = base nominal speed
i = nominal speed of the ith component on the rotating shaft
(at this gear ratio)
GRi = gear ratio of the ith component
1-64
Dynamic Theory
1-65
(1.33)
where:
= rotation speed (radians/unit time)
The rotational speed above and for friction loss is in radian
units not revolutions. When the rotating equipment slows
down, the stored rotational energy is either given back to
the fluid or lost to friction.
(1.34)
You can also calculate the inertia by sizing method. The sizing
method uses the following equation to calculation the inertia:
P 0.9556
I = 0.03768 -----3-
N
(1.35)
where:
I = inertia
P = design power
N = design speed
The sizing method for calculating inertia is only available in
the Pump unit operation.
1-65
1-66
HYSYS Dynamics
Friction Loss
The rate of energy lost from mechanical inefficiencies depends
on the frictional power loss factor, ffric, which can be specified in
the Friction Loss group. The frictional work, Pf, can be calculated
as follows:
P f = f fric
(1.36)
c = I f fric
(1.37)
The user can also enter the time constant directly, and the
friction loss factor is calculated from the inertia.
The approach to tuning the friction loss factor for rundown
rate may produce an inaccurate power loss at design
conditions because there is a short-circuit power loss of the
fluid at zero load, which is not currently taken into account1.
Dynamic Theory
1-67
(1.38)
(1.39)
1-67
1-68
HYSYS Dynamics
Description
Enable static
head
contributions
Enable Implicit
Static Head
The Static Head page, found on the Dynamics tab of the Main
Tray Section in the column environment, allows you to choose
the calculation method used to calculate the static head for this
operation. There are four options given on this page:
Dynamic Theory
1-69
For streams:
1-69
1-70
References
1.6.6 Design
The Design page, on the Rating tab of the Pump unit operation,
lets you specify the following variables in the Design Flow
group:
Figure 1.23
1.7 References
1
1-70
Dynamic Tools
2-1
2 Dynamic Tools
2.1 Introduction................................................................................... 3
2.2 Dynamics Assistant........................................................................ 4
2.2.1
2.2.2
2.2.3
2.2.4
2.2.5
2.2.6
2.2.7
2-1
2-2
Dynamic Tools
Theory ..................................................................................45
Event Scheduler Property View .................................................47
Sequence Property View ..........................................................53
Event Property View ................................................................57
Analyzing a Schedule...............................................................70
Running a Schedule.................................................................71
2-2
Dynamic Tools
2-3
2.1 Introduction
Modeling a process in dynamics is a complex endeavour. From
the perspective of defining the model, you must consider
parameters such as vessel holdups, valve sizing, and use of
pressure flow specifications.
To help simplify this process, HYSYS has several dynamic tools
and they are described in the following table:
Tool
Description
Dynamics
Assistant
View Equations
Tool
Integrator
Event Scheduler
2-3
2-4
Dynamics Assistant
Original
Configuration
Configuration and
Pressure Flow Specifications
after using Dynamics Assistant
2-4
Dynamic Tools
2-5
Original Configuration
2-5
2-6
Dynamics Assistant
Figure 2.3
vessel volumes
valve Cvs
k values (for equipment such as heaters, coolers, and
heat exchangers)
2-6
Dynamic Tools
2-7
2-7
2-8
Dynamics Assistant
Pressure Specs
Flow Specs
Uninitialized
Insert Valves
Int. Flow Spec
2-8
Dynamic Tools
2-9
2-9
2-10
Dynamics Assistant
temporary measure.
Figure 2.5
Uninitialized Page
The Uninitialized page contains the list of streams which are not
completely defined. For HYSYS to initialize any streams listed on
this page, you must be in Dynamics mode; uninitialized streams
2-10
Dynamic Tools
2-11
2-11
2-12
Dynamics Assistant
2-12
Dynamic Tools
2-13
valve
compressor
pump
heater
2-13
2-14
Dynamics Assistant
PF versus DP Page
The PF versus DP page lists the unit operations which currently
have a specified pressure drop as the dynamic specification. The
pressure drop option should not be used because this is
physically unrealistic. Material flow is driven by pressure
differences as well as resistances and stops when the pressures
have been equalized.
Figure 2.9
2-14
Dynamic Tools
2-15
LNG Page
The LNG page indicates which LNG exchangers are currently
specified with a pressure drop specification or which LNGs are
missing k values (depending on the dynamic rating method
chosen). Pressure drop specifications should be changed to
either pressure flow equation specifications (k values) or
suitable pressure drop correlations.
Figure 2.10
Valves
Volumes
K values
Valves Page
The Valves page lists the valves that are not sized. The current
conditions for the valve are listed and the calculated valve Cv
based on the pressure drop and percent opening of the valve,
both of which are changed directly on the page.
2-15
2-16
Dynamics Assistant
It is possible to change any of the sizing data for the valve. The
Cv value is updated based on any changes that are made. If you
change Cv, the new value is added to the valve when the Make
Changes button is clicked.
If you do not want HYSYS to size a valve, clear the OK checkbox
for the valve.
2-16
Dynamic Tools
2-17
Volumes Page
The Volumes page lists the unit operations which have unknown
volumes.
Figure 2.12
2-17
2-18
Dynamics Assistant
K Values Page
For more information on
k values or other
pressure flow
parameters, see
Chapter 1.4.2 - Basic
Pressure Flow
Equations.
The k values page lists unit operations for which the k value is
unknown. The information required to calculate the k value is
listed along with the current calculated k value.
Unit which can have k values include:
Heaters
Coolers
Heat exchangers
Air coolers
Valves
Figure 2.13
SS Pressures
Dry Hole Losses
Stream Connection
2-18
Dynamic Tools
2-19
SS Pressures Page
The SS Pressure page identifies tray sections where the total
steady state pressure drop seems to be inconsistent with the
total pressure drop calculated according to the dynamics rating
model.
Figure 2.14
2-20
Dynamics Assistant
The Sizing utility for selected sections button opens the tray
sizing utility view which can be used to minimize the pressure
losses.
2-20
Dynamic Tools
2-21
Misc Specs
Component Splitter
2-21
2-22
Dynamics Assistant
Some examples of the types of changes that are made are listed
below:
2-22
Dynamic Tools
2-23
2-24
Dynamics Assistant
Multiple Connections
Conflicts
Unit Operations
Ignored Opers
Flow direction
2-24
Dynamic Tools
2-25
Conflicts Page
The Conflicts page lists any streams that have their flows
directly controlled by controllers.
Figure 2.20
You can open the property view for a given unit operation by
double-clicking on its name.
2-25
2-26
Dynamics Assistant
2-26
Dynamic Tools
2-27
The Flow direction page is only active when static heads are
disabled.
2-28
Dynamic Tools
2-29
2-29
2-30
2-30
Dynamic Tools
2-31
2-31
2-32
2-32
Dynamic Tools
2-33
The Type cell displays the type of specified equation such as:
Pressure Balance Equation or Flow Balance Equation.
This number
corresponds to
the number
displayed in
the Number of
Variables cell
on the
Summary tab.
2-33
2-34
This number
corresponds
to the
number
displayed in
the User Spec
Vars cell on
the Summary
tab.
2-34
Dynamic Tools
2-35
2-35
2-36
2-36
Dynamic Tools
2-37
2.3.11 Simultaneous
Equations Tab
The Simultaneous Equations (SimulEqns) tab lists the equations
that are solved simultaneously by the integrator.
Figure 2.34
The Type cell displays the type of specified equation such as:
Pressure Balance Equation or Flow Balance Equation.
2-37
2-38
Integrator
2.4 Integrator
The Integrator is used when running a case in Dynamics mode.
You can access the Integrator property view by one of the
following methods:
Figure 2.35
Dynamic Tools
2-39
Once the integrator has executed the number of time steps, the
integrator does not stop, but remains in a holding mode. If
additional time steps are entered, the integrator continues
integration for the given number of time steps.
2-39
2-40
Integrator
Description
Pressure Flow
Solver
Control and
Logical Ops
Energy
Calculations
Composition and
Flash
Calculations
You can apply the execution rate values setting to all operations
by selecting the Use these default periods for all operations
checkbox. Clear the checkbox if you do not want to apply the
execution rate values entered in this group.
You have the option of specifying the composition and energy
balance execution rates per integration time step for individual
dynamic unit operations.
To specify individual execution rates for different unit
operations, it is necessary to add a new Dynamic Equipment
Ops tab to the Workbook.
1. Clear the Use these default periods for all operations
checkbox in the Calculation Execution Rates group of the
Integrator property view.
2. Select the Workbook. Select Workbook from the menu bar,
and then select Setup. This opens the Setup property view.
2-40
Dynamic Tools
2-41
3. Click the Add button in the Workbook Tabs group. The New
Object Type property view automatically opens.
Figure 2.37
2-41
2-42
Integrator
Parameter
Description
Enable implicit
static head
calculation
2-42
Dynamic Tools
2-43
Parameter
Description
Singularity
pressure flow
analysis before
running
Rigorous non
equilibrium mixed
properties
Simultaneously
solve heat transfer
eqns with IOFlash
Model choking of
liquid inside the
valve
2-43
2-44
Integrator
Parameter
Description
Truncate large
volume integration
errors
Reduced recycle
efficiency for small
timesteps
Close component
material and
energy balance
2-44
Dynamic Tools
2-45
2.5.1 Theory
The Event Scheduler property view as shown in Figure 2.42
contains all the Event Schedules in the current HYSYS case.
Each Schedule is comprised of Sequences, which in turn are
made up of Events. An Event must have a Condition (a null
2-45
2-46
Event Scheduler
Schedule 1
Schedule 2
Schedule 3
Schedule N
Sequence A
Event 1
Condition
Action List
Action 1
Action M
Event X
Sequence Z
2-46
Dynamic Tools
2-47
2-47
2-48
Event Scheduler
2-48
Dynamic Tools
2-49
Description
Add
Delete
Copy
Import
Export
Sort
The Legend group displays the status icons, which are shown as
the tree icons in the tree browser, and indicate the current state
of the schedule, sequence, event, or action.
Figure 2.43
The following states are valid:
Schedule: Fully Specified and
Incomplete
Sequence: Complete, Holding, Inactive,
Incomplete and Waiting.
Event: Complete, Fully Specified,
Holding, Inactive, Incomplete, Running,
Time Elapsed and Waiting.
Action: Fully Specified and Incomplete.
2-49
2-50
Event Scheduler
Description
Smart Tree
Trace
Messages
Multi Events
2-50
Dynamic Tools
2-51
Description
Sequence
Run Mode
Status
Event
Waiting For
Pending
Actions
Description
View
Add
Delete
Copy
Import
Export
Sort
2-51
2-52
Event Scheduler
Description
Start
Stop
Resume
Hold
Force
Skip
The Schedule Name field allows you to modify the name of the
schedule.
You can view the status of the selected schedule by clicking on
the Status Panel button although the individual schedule
property view is not used by most users. The individual schedule
property view also allows you to control the schedule by using
the Start, Stop, Resume, and Hold buttons.
Figure 2.45
2-52
Dynamic Tools
2-53
You can also click the View button from the Schedule
Sequences group of the Event Scheduler property view to
view the Sequence property view.
Schedule of Events
Settings
2-53
2-54
Event Scheduler
Description
Event
Specified
Condition
Action List
Jump When
Jump To
2-54
Dynamic Tools
2-55
Button
Description
View
Add
Delete
Copy
Sort
Analyze
Settings Tab
On the Settings tab, you can specify the sequences universal
settings, status window reporting, and the event conditions
default timeout behaviour.
Figure 2.47
2-55
2-56
Event Scheduler
Options
Description
Run Mode
You can select the sequence mode from the dropdown list. There are two options:
One Shot. Executes all its Events in order
then changes its status to Complete.
Continuous. Returns to the first Event after
the last Event has executed in a continuous
loop.
You can also set the Run Mode from the Event
Scheduler property view as shown in Figure 2.42.
You can also specify the default time out behaviour, which
applies to all events unless a specific event overrides this
behaviour.
2-56
Dynamic Tools
2-57
You can also click the View button from the Schedule of
Events tab of the Sequence property view to view the Event
property view.
2-57
2-58
Event Scheduler
The Force and Skip buttons affect the current event which
may not be the event you are viewing.
Condition
Action List
Branching & Time Out Behaviour
Condition Tab
The Condition tab as shown in Figure 2.48 shows the four
possible conditions located within the Wait For group. The
condition the user selects determines:
2-58
Dynamic Tools
2-59
2-59
2-60
Event Scheduler
A Variable to Stabilize
The actions execute after a variable has stabilized.
Figure 2.55
Dynamic Tools
2-61
The name and type of the action selected in the List Of Actions
For This Event group can be changed in the Individual Action
Specification group.
The Name field allows you to change the action name while the
Type drop-down list provides the available action types. After an
2-61
2-62
Event Scheduler
Specify Variable
The Specify Variable action requires an Object and a Value.
Figure 2.57
Start/Stop/Hold/Resume Sequence
The Start, Stop, Hold, and Resume Sequence actions all require
that a Sequence is selected from the drop-down list.
Figure 2.58
2-62
Dynamic Tools
2-63
The available Sequences are only those that are part of the
parent Schedule. You cannot control sequences of other
schedules.
Play Sound
The Play Sound action requires a *.wav audio file, which can be
selected by clicking the Select Audio File button. When a file is
selected, the Test sound button is enabled and the audio file can
be played.
Figure 2.59
Trace Dump
The Trace Dump action requires the selection of a source
variable by clicking the Select Source button.
Figure 2.60
2-63
2-64
Event Scheduler
Save Snapshot
For more information,
refer to Taking a
Snapshot in the Event
Scheduler section from
Section 11.12.3 External Snapshots in
the HYSYS User Guide.
2-64
Dynamic Tools
2-65
Service
Excel
Topic
Sheet1
Command
MyMacro
Stop Integrator
The Stop Integrator action stops the integrator.
Figure 2.64
2-65
2-66
Event Scheduler
Ramp Controller
The Ramp Controller action requires a Controller, a Target Set
Point, and a Ramp Duration as inputs.
Figure 2.65
2-66
Dynamic Tools
2-67
Trace Message
The Trace Message action traces a message to the HYSYS Trace
Window. The user specifies the text of the message in the
Message field.
Figure 2.67
2-67
2-68
Event Scheduler
Description
Never
Always
True
Timeout
False
2-68
Dynamic Tools
2-69
2-69
2-70
Event Scheduler
This property view displays a matrix with a list of all the Events
that are not fully specified. Selecting an Event and doubleclicking in the Event column or clicking the Analyze Event button
opens the Analysis of Event property view.
Figure 2.70
2-70
Dynamic Tools
2-71
2-71
2-72
Control Manager
If the Sequence RunMode is One Shot and the last event in the
list executes, the status changes from Waiting to Complete and
the Sequence is reset.
PIDs Tab
This tab provides a summary of the PID Controllers within the
current simulation. There are three modes displayed:
Dynamic Tools
2-73
There are also three key variables displayed: set point (Sp),
process variable (Pv), and operating target object (Op).
MPC Tab
This tab provides a summary of the MPC Controllers within the
current simulation. There are two modes that appear:
There are also three key variables displayed: set point (Sp),
process variable (Pv), and operating target object (Op).
2-73
2-74
Dynamic Initialization
and off-line.
Correct set up of the cold initialization configuration is not
only necessary for obtaining accurate results, but it can also
prevent movements in the model, and pressure-flow nonconvergence error.
It is recommended that you run the Integrator a few time
steps on any unit operations and streams before dynamic
initialization can have an effect on those operations. After
performing dynamic initialization, if the simulation
progresses through numerous steps without any errors or
significant movement in the process, the initialization is
successful.
Dynamic Tools
2-75
initialization conditions.
There are four buttons in the Dynamic Initialization Manager:
2-75
2-76
Dynamic Initialization
Cold Initialization
Notes
Figure 2.73
2-76
Dynamic Tools
2-77
Objects Page
The Objects page displays a list of objects that have been
selected for this dynamic initialization area. In the Selected
Objects table, the Name, and Type for each selected object is
shown.
Figure 2.74
You can rename
the dynamic
initialization area
in the Name field.
You can add more objects to the list by clicking the Add
Object button.
You can remove objects from the Selected Objects
table, by clicking the checkbox in the Remove column,
or by selecting the row containing the object and clicking
the Remove Object button.
You can remove all of the objects from the Selected
Objects table by click the Remove All button.
2-77
2-78
Dynamic Initialization
2-78
Dynamic Tools
2-79
Object
Description
Fluid Package
drop-down list
View Stream
button
Stream
Properties table
Initialize Objects
button
2-79
2-80
Dynamic Initialization
Configuration Page
The Configuration page allows you to activate the initialization
settings for specific equipment, logical operations, and stream
specifications for this dynamic initialization area.
Figure 2.77
2-80
Dynamic Tools
2-81
Checkbox
Description
Refer to Spreadsheets
section for more
information on Controller,
Digital Point, Selector,
and Transfer Function.
Controller
Initialization
Reset Digital
Point
Reset Transfer
Function
Reset Energy
Streams
Allows you to set the duty and any duty fluid flows to
0.
2-81
2-82
Dynamic Initialization
Checkbox
Reset Stream
Flow Specs
Reset Stream
Pressure Specs
Description
Allows you to reset all stream flow specifications to 0.
Stream flow specifications in HYSYS dynamics are not
a natural way for modeling real life behaviour since
flow usually occurs due to pressure gradients. You
should use the flow specifications with caution. It is
recommended to set the flow specifications to 0 for
cold initialization. It should be noted that when there is
a zero flow specification into a closed valve or flow
path, an absolute zero pressure is produced, and the
corrected pressure message appears in the trace
window.
You can select this checkbox to set the pressure equal
to that of the initialization stream (Atmospheric
pressure is commonly used).
If you insert isolation valves on boundary streams that
have a different pressure than the initialization stream,
then you do not need to select this checkbox. Most
users would prefer to preserve pressure specifications.
Refer to Special
Behaviours and Usage
for Unit Operations
during Dynamic
Initialization for
additional information on
valve behaviour during
dynamic initialization.
Reset Valves
Pump Off
Reset
Compressor/
Expander
Reset Heater/
Cooler specs
2-82
Dynamic Tools
2-83
Valves
When an isolation valve is initialized to be closed, but its
downstream equipment is not in this initialization area, the
valves downstream holdup is still initialized to that of the
initialization stream. As a result, when a valve is cold initialized
without including any of its downstream equipment, the valve
will be closed but its downstream holdup (especially for nonzero volume) will be set to the upstream initialization contents.
The Integrator must be started on a new model before
performing the dynamic initialization. If you are having
difficulties converging the pressure-flow solver or you just
want to quickly cold initialize your new model, just set the
Integrator in manual mode, and start the Integrator. You do
not need to take any steps. Now you can cold initialize.
2-83
2-84
Dynamic Initialization
Spreadsheets
The Spreadsheet operation is often used in control modeling.
However in dynamic initialization, a spreadsheet cell should NOT
be used for the start of control or logic propagation. Instead,
Controller, Selector, Transfer Function, or Digital Control Points
could be used as the operator accessible points for turning on or
off the plant. For instance, the OP state of a Digital Point can be
imported into a spreadsheet cell and then operated upon before
it is sent out to the On/Off switch of a pump. As long as the
spreadsheets do not contain any time dependent behaviour,
they should initialize with the inputs coming from one of the four
mentioned logical operations.
Event Schedules
It is recommended that all the schedules in the Event Scheduler
be reset manually to their OFF state before running the case in
dynamic initialization.
2-84
Dynamic Tools
2-85
2-85
2-86
Dynamic Initialization
2-86
Control Theory
3-1
3 Control Theory
3.1 Introduction................................................................................... 2
3.2 Process Dynamics .......................................................................... 3
3.2.1 Characteristic Parameters of the Process System .......................... 3
3.3 Basic Control................................................................................ 10
3.3.1 Available Control Operations .................................................... 16
3.4 Advanced Control......................................................................... 32
3.4.1 Model Predictive Control.......................................................... 32
3.5 General Guidelines ....................................................................... 39
3.5.1
3.5.2
3.5.3
3.5.4
3-1
3-2
Introduction
3.1 Introduction
Process control on a working level involves the control of
variables such as flowrate, temperature, and pressure in a
continuously-operating plant. Process control in a general sense
attempts to maximize profitability, ensure product quality, and
improve the safety and operating ability of the plant.
While steady state simulation in HYSYS allows the design
engineer to optimize operating conditions in the plant, dynamic
simulation allows you to:
3-2
Control Theory
3-3
3.2.1 Characteristic
Parameters of the Process
System
It is easiest to define a chemical process system using the
general conservation principle which states that:
Rate of accumulation = Input - Output + Internal
Generation
(3.1)
3-3
3-4
Process Dynamics
(3.2)
where:
H = liquid height in the tank
A = cross-sectional area of the tank
Fi = inlet flow rate
Fo = exit flow rate
(3.3)
3-4
Control Theory
3-5
(3.4)
dH
RA ------- + H = RF i
dt
(3.5)
(3.6)
where:
y(t) = output of the system
u(t)= input to the system
K = steady state gain
= time constant of the system
(3.7)
3-5
3-6
Process Dynamics
----t
y ( t ) = MK 1 e
(3.8)
H ( t ) = MR 1 e
t
------AR
(3.9)
3-6
Control Theory
3-7
Process Gain
The process gain is defined as the ratio of the change/deviation
in the process output to the change/deviation in the process
input. The change in the process input is defined in Equation
(3.6) as u(t). The change in the process output is defined as
y(t). The first term in Equation (3.6) is transient and becomes
zero at steady state. Therefore, the gain can be calculated as
shown in the equation below.
y SSnew y SS
- = K
Steady state gain = -----------------------------u SSnew u SS
(3.10)
where:
ySSnew = new steady state y
uSSnew = new steady state U
For this liquid level example, the steady state gain, K, is the
valve resistance, R. Therefore, a step change in the flow into the
tank of magnitude M results in a change in liquid level, H(t), in
the tank equal to MR.
Time Constant
The time constant, , defines the speed of the response. The
response of the system always follow the profile shown in
Figure 3.2. After time units, the response y(t) equals
0.632MK or 63.2% of the ultimate gain. This is always true for
first-order systems without time delays. For this liquid level
example, the time constant is the product of the area of the
tank, A, and the resistance of the exit valve, R.
3-7
3-8
Process Dynamics
Capacity
Definition 1
Capacity can be defined simply as the volume or storage space
of a system. The capacitance of a system dampens the output
causing the response to take time to reach a new steady state.
For electrical systems, the capacity is defined in terms of the
resistance of the system and the time constant of the response:
C = --R
(3.11)
Definition 2
A systems capacity is also defined as its ability to attenuate an
incoming disturbance. Attenuation is defined as:
Response Amplitude out of the systemAttenuation = 1 --------------------------------------------------------------------------------------------Disturbance Amplitude into the system
Attenuation = 1 Amplitude Ratio
(3.12)
3-8
Control Theory
3-9
The input, u(t), to the first order system can be changed from a
step function to a sinusoidal function:
u ( t ) = M sin ( t )
(3.13)
where:
= frequency of the input response
M = amplitude of the input function
e t / + ---------------------------- sin ( t + )
y ( t ) = MK ----------------------( ) 2 + 1
( ) 2 + 1
(3.14)
where:
1
= tan ( )
(3.15)
(3.16)
3-9
3-10
Basic Control
Dead Time
While capacitance is a measure of how fast a system responds
to disturbances, dead time is a measure of the amount of time
that elapses between a disturbance to the system and the
observed response in the system.
Time delays in a system can be significant depending on the
nature of the process and the location of measuring devices
around the process. It is usually the time associated with the
transport of material or energy from one part of the plant to
another that contributes to time delays observed in a system.
The dead time of a process is easily modeled using the Transfer
Function block operation.
3-10
Control Theory
3-11
Terminology
Before reviewing the major control operations that are available
in HYSYS, it is useful to be familiar with the following terms.
Disturbances
A disturbance upsets the process system and causes the output
variables to move from their desired set points. Disturbance
variables cannot be controlled or manipulated by the process
engineer. The control structure should account for all
disturbances that can significantly affect a process. The
disturbances to a process can either be measured or
unmeasured.
3-11
3-12
Basic Control
3-12
Control Theory
3-13
(3.17)
where:
OP(t) = controller output at time t
E(t) = error at time t
Kc = proportional gain of the controller
Ti = integral (reset) time of the controller
Td = derivative (rate) time of the controller
3-13
3-14
Basic Control
(3.18)
where:
SP(t) = set point
PV(t) = measured output process variable
If the PV rises above the SP, the OP, or input to the process,
decreases. If the PV falls below the SP, the OP increases.
For a process with a negative steady state gain, the error should
be set as direct acting:
E ( t ) = PV ( t ) SP ( t )
(3.19)
That is, if the PV rises above the SP, the OP, or input to the
process, increases. If the PV falls below the SP, the OP
decreases.
A typical example of a reverse acting controller is in the
temperature control of a reboiler. In this case, as the
temperature in the vessel rises past the SP, the OP decreases, in
effect closing the valve and hence the flow of heat.
3-14
Control Theory
3-15
Stability
The stability of a system is a very important aspect to consider
when designing control schemes. Most systems can have
oscillatory responses, depending on its controller tuning
parameters. When a process is upset by a bounded disturbance
or bounded change in the input forcing function, the output
typically responds in one of three ways:
Figure 3.6
3-15
3-16
Basic Control
Tuning parameters
Non-linearities in the process
Range and non-linearities in the instruments
Interactions between control loops
Frequency of disturbance
Capacity of process
Noise in measurement of process variables
3-16
Control Theory
3-17
Sensors
Sensors are used to measure process variables. In HYSYS, the
sensing instrument is incorporated directly in the PID Controller
operation. You can choose the range of the sensing instrument
in the Min and Max PV parameters in the controller operation. It
is assumed in HYSYS that the PID controller is perfectly accurate
in its measurement of the process variable.
Valve Mode
Description
Instantaneous
First Order
Linear
3-17
3-18
Basic Control
Description
Linear
0.5
Linear
60
% Cv
40
Equal
Percentage
20
0
0
20
40
60
80
100
% Valve Position
3-18
Control Theory
3-19
Feedback Control
Digital On/Off
Digital On/Off control is one of the most basic forms of
regulatory control. In HYSYS, it is implemented using the Digital
Point operation. An example of On/Off control is a home heating
system. When the thermostat detects that the temperature is
below the set point, the heating element turns on. When the
temperature rises above the set point, the heating element
turns off.
For more information on
the Digital Point
operation in HYSYS, see
Section 5.5 - Digital
Point in the HYSYS
Operations Guide.
(3.20)
(3.21)
3-19
3-20
Basic Control
3-20
Control Theory
3-21
(3.22)
3-21
3-22
Basic Control
3-22
Control Theory
3-23
(3.23)
The integral term serves to bring the error to zero in the control
scheme. The more integral action there is, the slower the
response of the controller. The integral term continuously
moves to eliminate the error. The closed loop response of a
process with PI control and P-only control is shown as follows:
Figure 3.10
3-23
3-24
Basic Control
(3.24)
Control Theory
3-25
3-25
3-26
Basic Control
Cascade Control
Cascade control is a technique that implements a secondary
feed back control loop within a primary feedback control loop.
Cascade control can be used when there are significant
disturbances to the manipulated variable of the primary loop.
A secondary loop is created to control the manipulated variable
of the primary loop. The primary loop then manipulates the set
point of the secondary controller.
Figure 3.12
Control Theory
3-27
3-27
3-28
Basic Control
Feedforward Control
Feedforward control can be used in cases for which feedback
control cannot effectively control a process variable. The main
disadvantage of feedback control is that the controller must wait
until disturbances upset the process before responding.
With feedforward control, the controller can compensate for
disturbances before the process is affected. Cascade control is
useful when measured disturbances significantly affect the input
3-28
Control Theory
3-29
3-29
3-30
Basic Control
(3.25)
where:
Fs = steam flow
= heat of condensation for steam
F = flow of stream Exit
Cp = specific heat of stream Exit
3-30
Control Theory
3-31
From this process model, the desired value of steam flow into
the heat exchanger can be calculated. The flow of steam must
be calculated such that the temperature of stream Exit, T2,
equals the desired temperature, SP. Therefore, Equation
(3.25) becomes:
Cp
F s = ------- F ( SP T 1 )
(3.26)
(3.27)
3-32
Advanced Control
chemicals
food processing
automotive
aerospace
metallurgy
pulp and paper
3-32
Control Theory
3-33
MPC Theory
Currently most model predictive control techniques like Dynamic
Matrix Control (DMC) and Model Algorithmic Control (MAC) are
based on optimization of a quadratic objective function involving
the error between the set point and the predicted outputs. In
these cases, a discrete impulse response model can be used to
derive the objective function.
Let a0, a1, a2,...,aT represent the value of the unit step response
function obtained from a typical open loop process, as shown in
the figure below:
Figure 3.16
3-33
3-34
Advanced Control
cn = cn
(3.28)
c n + 1 = c 0 +
ai mn + 1 i
(3.29)
i=1
Figure 3.17
cn + 1 = c0 +
hi mn + 1 i
(3.30)
i=1
3-34
Control Theory
3-35
cn = c0 +
hi mn i
(3.31)
i=1
c n + 1 = c n +
hi mn + 1 i
(3.32)
i=1
c n* + 1 c n + 1 = cn c n
(3.33)
hi mn + 1 i
(3.34)
i=1
3-35
3-36
Advanced Control
hi mn + j i
(3.35)
i=1
where:
j = 1, 2, ..., V
a1
c n* + 2
a2
a1
a3
a2
c n* + 3
.
.
.
c n* + V
m n
cn + P1
m n + 1
cn + P2
a1
m n + 2
.
.
.
.
.
.
.
.
.
a V a V 1 a V 2 a V U + 1 m n + U 1
cn + P3
.
.
.
cn + PV
(3.36)
where:
ai =
hj
(3.37)
j=1
i
Pi =
Sj
for
i = 1 , 2 , , V
(3.38)
j=1
Si =
h i m n + j i
for
j = 1 , 2 , , V
(3.39)
i = j+1
Control Theory
3-37
E = A m + E'
(3.40)
where:
A = the V U triangular matrix
m = the U 1 vector of future control moves.
and E'
= the closed loop and open loop predictions,
E
respectively, and are defined as follows:
r n + 1 c n* + 1
rn + 2
=
E
En P1
c n* + 2
.
.
.
rn + V cn + V
En P2
E'
=
.
.
.
En PV
(3.41)
3-37
3-38
Advanced Control
(3.42)
m = ( A T A ) A T E' = K c E'
(3.43)
where:
1
( A T A ) A T = pseudo-inverse matrix
K c = matrix of feedback gains (with dimensions V U )
J [ m ] = E T u E + m T y m
(3.44)
where:
u and y are positive-definite weighting matrices for
predicted errors and control moves, respectively. These
matrices allows you to specify different penalties to be
placed on the predicted errors resulting in a better
tuned controller.
m = ( A T u A + y ) A T u E' = K c E'
(3.45)
Control Theory
3-39
Capacity
The ability of a system to attenuate incoming disturbances is a
function of the capacitance of a system and the period of the
disturbances to the system. From Terminology section,
attenuation is defined as:
K
Attenuation = 1 --------------------------( ) 2 + 1
(3.46)
3-39
3-40
General Guidelines
Dead Time
The dead time has no effect on attenuating disturbances to open
loop systems. However, it does have a significant negative
effect on controllability. Dead time in a process system reduces
the amount of gain the controller can implement before
encountering instability. Because the controller is forced to
reduce the gain, the process is less able to attenuate
disturbances than the same process without dead time.
Tight control is possible only if the equivalent dead time in the
loop is small compared to the shortest time constant of a
disturbance with a significant amplitude. It is generally more
effective to reduce the dead time of a process than increase its
capacity.
To reduce dead time:
3-40
Control Theory
3-41
3-41
3-42
General Guidelines
3-42
Control Theory
3-43
Flow Control
Flow in a pipe is typically a fast responding process. The dead
time and capacity associated with a length of pipe is generally
small. It is therefore not unusual for the process to be limited by
the final control element (valve) dynamics. You can easily
incorporate valve dynamics in the HYSYS model by modifying
the valve parameters in the Actuator page of the Dynamics tab.
Tuning a flow loop for PI control is a relatively easy task. For the
flow measurement to track the set point closely, the gain, Kc,
should be set between 0.4 and 0.65 and the integral time, Ti,
should be set between 0.05 and 0.25 minutes.
Since the flow control is fast responding, it can be used
effectively as the secondary controller in a cascade control
structure. The non-linearity in the control loop can cause the
control loop to become unstable at different operating
conditions.
Therefore, the highest process gain should be used to tune the
controller. If a stability limit is reached, the gain should be
decreased, but the integral action should not. Since flow
measurement is naturally noisy, derivative action is not
recommended.
3-43
3-44
General Guidelines
3-44
Control Theory
3-45
Temperature Control
Temperature dynamic responses are generally slow, so PID
control is used. Typically, the controller gain, Kc, should be set
between 2 and 10, the integral time, Ti, should set between 2
and 10 minutes, and the derivative time Td, should be set
between 0 and 5 minutes.
Tuning Methods
An effective means of determining controller tuning parameters
is to bring the closed loop system to the verge of instability. This
is achieved by attaching a P-only controller and increasing the
gain such that the closed loop response cycles with an amplitude
that neither falls nor rises over time. At a systems stability
margins, there are two important system parameters, the
ultimate period and the ultimate gain, which allow the
calculation of the proportional, integral, and derivative gains.
3-45
3-46
General Guidelines
3-46
Control Theory
3-47
Equation
Ultimate Gain
Ultimate Period
4h
K u = -----a
P u = Period taken from limit cycle
Controller Gain
K
K c = ------u3.2
T i = 2.2P u
Equation
Ultimate Gain
3-47
3-48
General Guidelines
Control Theory
3-49
3-49
3-50
General Guidelines
Connections Tab
In the Connections tab, you can specify/select the variable
information entering and exiting the controller.
3-50
Control Theory
Operating
Variable
3-51
Description
Physical Valve
Material Stream
OP ( % )
Flow = ------------------ ( Flow max Flow min ) + Flow min
100
HYSYS varies the flow specification of the material
stream according to the calculated controller output,
OP. (Therefore, a non-realistic situation can arise in
the dynamic case since material flow is not dependent
on the surrounding conditions.)
Energy Stream
OP ( % )
Energy Flow = ------------------ ( Flow max Flow min ) + Flow min
100
If the Utility Fluid option is chosen, you need to specify
the maximum and minimum flow of the utility fluid.
The heat flow is then calculated using the local overall
heat transfer coefficient, the inlet fluid conditions, and
the process conditions.
3-51
3-52
General Guidelines
Parameters Tab
The direction of the controller, the controllers PV range, and the
tuning parameters can be specified in the Parameters tab.
For more information
about whether a
controller is direct or
reverse acting, refer to
Terminology section.
3-52
Control Theory
3-53
From the Strip Charts tab, add a new strip chart by clicking the
Add button and activate the variables to be displayed on the
strip chart. No more than six variables should be selected for
each strip chart to keep it readable.
Figure 3.21
Click on the Strip Chart button in the View group to see the strip
chart. Size as desired and then right-click on the strip chart.
Select Graph Control command from the Object Inspect menu.
3-53
3-54
General Guidelines
There are six tabs, where you can manipulate the strip chart
display features, set the numerical ranges of the strip chart for
each variable, the nature of the lines for each variable, and how
the strip chart updates and plots the data.
Add additional strip charts as desired by going back into the
Databook property view and going to the Strip Charts tab.
3-54
Control Theory
3-55
3-55
3-56
General Guidelines
Control Theory
3-57
Stability
It is shown that the stability of a closed loop process depends on
the controller gain. If the controller gain is increased, the closed
loop response is more likely to become unstable. The controller
gain, Kc, input in the PID Controller operation in HYSYS is a
unitless value defined in Equation (3.47).
PV Range
K c = OP%
-------------------------------------------error
(3.47)
(3.48)
3-57
3-58
References
3.6 References
1
3-58
Index
A
Accumulation 1-14
Advanced Holdup
componsitions tab 1-29
efficiencies tab 1-27
general tab 1-25
nozzles tab 1-26
properties tab 1-28
Advanced Holdup Properties 1-24
Ambient 2-44
ATV Tuning 3-46
B
Balance
See Material Balance, Component Balance and Energy Balance
C
Capacity 3-8, 3-39
Cascade Control 3-26
See also PID Controller
Cold Initialization 2-76
adding objects 2-77
configuring objects 2-79
removing objects 2-77
selecting fluid package 2-79
Component Balance 1-9
Control Strategy 3-48
Controller
available control operations 3-16
choosing correct 3-41
selecting variables 3-50
tuning 3-42
Controller Theory
capacity 3-8
dead time 3-10
process gain 3-7
terminology 3-11
time constant 3-7
Cv
See Valve Flow Coefficient (Cv)
D
Databook 3-53
DDE 2-65
Dead Time 3-10, 3-40
Degrees of Freedom 1-34
Distributed Models 1-6
Dynamic Assistant
general tab 2-7
other specs tab 2-21
pressure flow specs tab 2-14
streams tab 2-8
unknown sizing tab 2-15
user items tab 2-24
Dynamic Initialization 2-73
See also Cold Initialization
Dynamic Simulation
control strategy 3-48
converting steady state models 1-44
degrees of freedom 1-34
differences from steady state 1-43
general concepts 1-6
linear 1-7
non-linear 1-7
theory 1-3, 1-6
Dynamics Assistant 2-4
E
Efficiencies 1-16
Energy Balance 1-10, 1-20
Equation Summary
property view 2-28
Event Scheduler 2-45
F
Face Plate 3-54
Feedback Control 3-12, 3-19
direct acting 3-14
reverse acting 3-14
Feedforward Control 3-28
Flash
non-equilibrium 1-15
Flow Control 3-43
Flowsheet Menu
dynamic initialization 2-73
G
Gas Pressure Control 3-45
H
Heat Loss Model 1-20
detailed 1-23
parameters 1-22
simple 1-22
Holdup Model 1-13
I-1
I-2
advantages of 1-13
assumptions 1-14
I
Implicit Euler Method 1-11
Integration Strategy 1-12
Integrator 2-38, 3-55
K
k Values 1-32
L
Liquid Level Control 3-44
Liquid Pressure Control 3-43
Lumped Models 1-6
M
Material Balance 1-7
N
Nozzles 1-18
O
Open Loop Control 3-11
Operations
general guidelines 1-42
Ordinary Differential Equations 1-7
P
Partial Differential Equations 1-6
PID Controller
ATV tuning 3-46
choosing correct controller 3-41
See also Controller Theory
tuning 3-42, 3-56
Ziegler-Nichols tuning 3-47
Pressure Flow
failed convergence 1-48
model 1-31
specifications 1-34, 1-37
volume balance 1-31
Pressure Flow Solver 1-29
simultaneous solution 1-30
Process Dynamic 3-3
Process Gain 3-7
Proportional Control 3-21
Proportional Integral Control (PI) 3-23
Proportional Integral Derivative Control (PID)
3-24
R
Reactions 1-23
Resistance Equation 1-33
S
Scripts 2-64
Sensors 3-17
Singular Problem 1-48
Solution
Implicit Euler 1-11
Stability 3-15, 3-57
Static Head Contributions 1-19
T
Temperature Control 3-45
Time Constant 3-7
Tuning Methods 3-45
V
Valve
fail-safe function 1-60
modes 1-58
positions 1-61
stickiness 1-58, 1-60
Valve Dynamics 3-17
Valve Flow Coefficient (Cv) 1-33
Valve Type 3-17
Volume Balance 1-31
Z
Ziegler-Nichols Tuning Technique 3-47
I-2