Process Simulation and Control Using Aspen
Process Simulation and Control Using Aspen
PROCESS SIMULATION
AND CONTROL USING
METHANOl
BUTENES
RDCOLUMN CCS
AMIYA K. JANA
This
Thispage
pageintentionally
intentionallyleft
leftblank
blank
Rs. 295.00
@ 2009 by PHI Learning Pnvate Limited, New Delhi. All rights reserved. No part of this book may
be reproduced In any form, by mimeograph or any other means, without permission in writing from
the publisher.
ISBN-978-81-203-3659-9
The export rights of this book are vested solely with the publisher.
Published by Asoke K. Ghosh, PHI Learning Private Limited, M-97, Connaught Circus,
New Delhi-110001 and Printed by Jay Print Pack Private Limited, New Delhi-110015.
Rs. 295.00
@ 2009 by PHI Learning Pnvate Limited, New Delhi. All rights reserved. No part of this book may
be reproduced In any form, by mimeograph or any other means, without permission in writing from
the publisher.
ISBN-978-81-203-3659-9
The export rights of this book are vested solely with the publisher.
Published by Asoke K. Ghosh, PHI Learning Private Limited, M-97, Connaught Circus,
New Delhi-110001 and Printed by Jay Print Pack Private Limited, New Delhi-110015.
This
Thispage
pageintentionally
intentionallyleft
leftblank
blank
Rs. 295.00
@ 2009 by PHI Learning Pnvate Limited, New Delhi. All rights reserved. No part of this book may
be reproduced In any form, by mimeograph or any other means, without permission in writing from
the publisher.
ISBN-978-81-203-3659-9
The export rights of this book are vested solely with the publisher.
Published by Asoke K. Ghosh, PHI Learning Private Limited, M-97, Connaught Circus,
New Delhi-110001 and Printed by Jay Print Pack Private Limited, New Delhi-110015.
This
Thispage
pageintentionally
intentionallyleft
leftblank
blank
Rs. 295.00
@ 2009 by PHI Learning Pnvate Limited, New Delhi. All rights reserved. No part of this book may
be reproduced In any form, by mimeograph or any other means, without permission in writing from
the publisher.
ISBN-978-81-203-3659-9
The export rights of this book are vested solely with the publisher.
Published by Asoke K. Ghosh, PHI Learning Private Limited, M-97, Connaught Circus,
New Delhi-110001 and Printed by Jay Print Pack Private Limited, New Delhi-110015.
Contents
Preface vii
Acknowledgements ix
Reference 53
Reference 106
v
Copyrighted maierlal
VI CONTENTS
f
Part II Chemical Plant Simulation using Aspen Plus
4 . Aspen Plus Simulation of Chemical Plants 189-226
4 1 TntrnHnrtion
4 . 2 Aspen Plus Simulation of a Distillation Train 189
4 . 3 Aspen Plus Simulation of a Vinyl Chloride Monomer (VCM)
Production Unit 203
Summary and Conclusions 220
Prnhlpms ; , -
220
References 226
Index 315-317
Copyrlghled maierlal
r
Preface
"
The future success of the chemical process industries mostly depends on the ability to
design and operate complex, highly interconnected plants that are profitable and that
"
meet quality, safety, environmental and other standards To achieve this goal, the software
.
tools for process simulation and optimization are increasingly being used in industry.
By developing a computer program, it may be manageable to solve a model structure
of a chemical process with a small number of equations. But as the complexity of a plant
integrated with several process units increases, the solution becomes a challenge. Under
this circumstance, in recent years, we motivate to use the process flowsheet simulator to
solve the problems faster and more reliably. In this book, the Aspen software package
has been used for steady state simulation, process optimization, dynamics and closed-
loop control.
To improve the design, operability, safety, and productivity of a chemical process
with minimizing capital and operating costs, the engineers concerned must have a solid
knowledge of the process behaviour. The process dynamics can be predicted by solving
the mathematical model equations. Within a short time period, this can be achieved
quite accurately and eficiently by using Aspen lowsheet simulator. This software tool is
f
not only useful for plant simulation but can also automatically generate several control
structures, suitable for the used process flow diagram. In addition, the control parameters,
including the constraints imposed on the controlled as well as manipulated variables.
are also provided by Aspen to start the simulation run. However, we have the option to
modify or even replace them.
This well organized book is divided into three parts. Part I (Steady State Simulation
and Optimization using Aspen Plus ) includes three chapters. Chapter 1 presents the
introductory concepts with solving the lash chambers. The computation of bubble point
f
and dew point temperatures is also focused. Chapters 2 and 3 are devoted to simulation
of several reactor models and separating column models, respectively.
Part II (Chemical Plant Simulation using Aspen Plus ) consists of only one chapter
(Chapter 4). It addresses the steady state simulation of large chemical plants. Several
individual processes are interconnected to form the chemical plants. The Aspen Plus
simulator is used in both Part I and Part II.
vii
Copyrighted maierlal
viii PREFACE
The Aspen Dynamics package is employed in Part III (Dynamics and Control using
Aspen Dynamics ) that comprises Chapters 5 and 6. Chapter 5 is concerned with the
dynamics and control of low-driven chemical processes. In the closed-loop control study ,
f
the servo as well as regulatory tests have been conducted. Dynamics and control of
pressure-driven processes have been discussed in Chapter 6.
The target readers for this book are undergraduate and postgraduate students of
chemical engineering. It will be also helpful to research scientists and practising engineers.
Amiya K. -Jana
Copyrighted maierlal
Acknowledgements
rial
Part I
Copyrigf
This
Thispage
pageintentionally
intentionallyleft
leftblank
blank
Rs. 295.00
@ 2009 by PHI Learning Pnvate Limited, New Delhi. All rights reserved. No part of this book may
be reproduced In any form, by mimeograph or any other means, without permission in writing from
the publisher.
ISBN-978-81-203-3659-9
The export rights of this book are vested solely with the publisher.
Published by Asoke K. Ghosh, PHI Learning Private Limited, M-97, Connaught Circus,
New Delhi-110001 and Printed by Jay Print Pack Private Limited, New Delhi-110015.
CHAPTER
i
Introduction and Stepwise
Aspen Plus Simulation:
Flash Drum Examples
11
. ASPEN: AN INTRODUCTION
Copynghled material
4 PROCESS SIMULATION AND CONTROL USING ASPEN
a single unit or a chemical plant, required Fortran codes are built-in in the Aspen
simulator. Additionally, we can also use our own subroutine in the Aspen package.
The Aspen simulation package has a large experimental databank for
thermodynamic and physical parameters. Therefore, we need to give limited input data
for solving even a process plant having a large number of units with avoiding human
errors and spending a minimum time.
Aspen simulator has been developed for the simulation of a wide variety of
processes, such as chemical and petrochemical, petroleum refining, polymer, and coal-
based processes. Previously, this flowsheet simulator was used with limited
applications. Nowadays, different Aspen packages are available for simulations with
promising performance. Briefly, some of them are presented below.
Aspen Plus-This process simulation tool is mainly used for steady state simulation of
chemicals, petrochemicals and petroleum industries. It is also used for performance
monitoring, design, optimization and business planning.
Aspen Dynamics-This powerful tool is extensively used for dynamics study and closed-
loop control of several process industries. Remember that Aspen Dynamics is integrated
with Aspen Plus.
Aspen BatchCAD-This simulator is typically used for batch processing, reactions and
distillations. It allows us to derive reaction and kinetic information from experimental
data to create a process simulation.
Aspen Chromatography-This is a dynamic simulation software package used for both
batch chromatography and chromatographic simulated moving bed processes.
Aspen Properties-It is useful for thermophysical properties calculation.
Aspen Polymers Plus-It is a modelling tool for steady state and dynamic simulation,
and optimization of polymer processes. This package is available within Aspen Plus or
Aspen Properties rather than via an external menu.
Aspen HYSYS-This process modelling package is typically used for steady state
simulation, performance monitoring, design, optimization and business planning for
petroleum refining, and oil and gas industries.
It is clear that Aspen simulates the performance of the designed process. A solid
understanding of the underlying chemical engineering principles is needed to supply
reasonable values of input parameters and to analyze the results obtained. For example, a
user must have good idea of the distillation column behaviour before attempting to use
Aspen for simulating that column. In addition to the process flow diagram, required input
information to simulate a process are: setup, components properties, streams and blocks.
,
f
icon on our desktop, the Aspen Plus Startup dialog appears. There are three choices
and we can create our work from scratch using a Blank Simulation, start from a
Template or Open an Existing Simulation. Let us select the Blank Simulation option
and click OK (see Figure 1.1).
MM
MM 'Ml I I-
FIGURE 1.1
The simulation engine of Aspen Plus is independent rom its Graphical User
f
Interface (GUI). We can create our simulations using the GUI at one computer and run
them connecting to the simulation engine at another computer. Here, we will use the
'
simulation engine at Local PC'. Default values are OK.
Hit OK in the Connect to Engine dialog (Figure 1.2). Notice that this step is specific
to the installation.
The next screen shows a blank Process Flowsheet Window. The irst step in
f
defined as a blueprint of a plant or part of it. It includes all input streams, unit
operations, streams that interconnect the unit operations and the output streams.
Several process units are listed by category at the bottom of the main window in a
toolbar known as the Model Library. If we want to know about a model, we can use the
Help menu from the menu bar. In the following, different useful items are highlighted
briefly (Figure 1.3).
Copyrighted material
6 PROCESS SIMULATION AND CONTROL USING ASPEN
Connect to Engine
Liter Into
Node name:
Uset name
Password
Working dfedory:
OK Exit Help
FIGURE 1.2
A*<>r rv l u» s*iitiil-('
r|ttRt..|:>|.>l rr
Al/lniAiAioj-MMBSF ZlF
Next button
Status bar
FIGURE 1.3
Copyrighted material
INTRODUCTION AND STKPWISK ASPEN PI.US SIMULATION 7
To develop a lowsheet, irst choose a unit operation available in the Model Library.
f
Proprietary models can also be included in the lowsheet window using User Models
f
option. Excel workbook or Fortran subroutine is required to define the user model. In
the subsequent step, using Material STREAMS icon, connect the inlet and outlet streams
with the process. A process is called as a block in Aspen terminology. Notice that clicking
on Material STREAMS, when we move the cursor into the lowsheet area red and blue
f
arrows appear around the model block. These arrows indicate places to attach streams
to the block. Red arrows indicate required streams and blue arrows are optional.
When the lowsheet is completed, the status message changes from Flowsheet Not
f
Complete to Required Input Incomplete. After providing all required input data using
input forms, the status bar shows Required Input Complete and then only the simulation
results are obtained. In the Data Browsery we have to enter information at locations
where there are red semicircles. When one has finished a section, a blue checkmark
appears. In subsection 1.3.2. a simple problem has been solved, presenting a detailed
stepwise simulation procedure in Aspen Plus. In addition, three more problems have
also been discussed with their solution approaches subsequently.
13
.
STEPWISE ASPEN PLUS SIMULATION OF FLASH DRUMS
In the Model Library, there are ive built-in separators. A brief description of these
f
two-phase (liquid-liquid) system. This separator includes two outlet liquid streams along
with inlet stream(s). It can be used as the separation columns. If there is any tendency
of vapour formation with two liquid phases, it is recommended to use Flash3 instead of
Decanter.
Sep 1: It is a multi-outlet component separator since two or more outlet streams can
be produced rom this process unit. It can be used as the component separation columns.
f
132 . .
Simulation of a Flash Drum
Problem statement
A 100 kmol/hr feed consisting of 10, 20, 30, and 40 mole% of propane, c-butane,
r
n-pentane and n-hexane, respectively, enters a lash chamber at 15 psia and 50oF.
,
f
The lash drum (Flash2) is shown in Figure 1.4 and it operates at 100 psia and 200oF.
f
Applying the SYSOP0 property method, compute the composition of the exit streams.
3 -
FLASH
f
Simulation approach
From the desktop, select Start button followed by Programs, AspenTech, Aspen
Engineering Suite, Aspen Plus Version and Aspen Plus User Interface. Then choose
Template option in the Aspen Plus Startup dialog (Figure 1.5).
I 1- l-MHM*
FIGURE 1.5
As the next window appears after hitting OK in the above screen, select General
with English Units (Figure 1.6).
Copyrighted material
INTRODUCTION AND STEPVV1SE ASPEN PLUSIM SIMULATION 9
-Hi 1
1 #
;1
L -
'
. I
i.-
I i -
FIGURE 1.6
Then click OK. Again, hit OK when the Aspen Plus engine window pops up and
subsequently, proceed to create the lowsheet.
f
Creating flowsheet
Select the Separators tab from the Model Library toolbar. As discussed earlier, there
are ive built-in models. Among them, select Flash2 and place this model in the window.
f
Now the Process Flowsheet Window includes the lash drum as shown in Figure 1.7. By
f
'
nia*lHl mU JM ??1
-
a-i-m * -ai-o "d 3 I l-l SI Hi bl'
r
- I --i
0 9 «=>. 8 - C . 1
FIGURE 1.7
Copyrlghled
10 PROCESS SIMULATION AND CONTROL USING ASPEN1
To add the input and output streams with the block, click on Streams section (lower
left-hand comer). There are three different stream categories (Material, Heat and Work),
as shown in Figure 1.8.
O,
, l , 1 Ma I J--
XQ.o Q lr
- -
FIGURE 1.8
Block Bl includes three red arrows and one blue arrow as we approach the block
after selecting the Material STREAMS icon. Now we need to connect the streams with
the lash chamber using red arrows and the blue arrow is optional. The connection
f
-
i- - rl ...iil il a ! 1
rmfT -1 "| LV I .(Bit
-
( - 11 iwl
-
- - I
. III MM .- . I -.
FIGURE 1.9
Copyrighled mateial
r
INTRODUCTION AND STFPWISK ASPEN PLUS SIMULATION 11
Clicking on Material STREAMS, move the mouse pointer over the red arrow at the
inlet of the lash chamber. Click once when the arrow is highlighted and move the
f
cursor so that the stream is in the position we want. Then click once more. We should
see a stream labelled 1 entering the drum as a feed stream. Next, click the red arrow
coming out at the bottom of the unit and drag the stream away and click. This stream
is marked as 2. The same approach has been followed to add the product stream at the
top as Stream 3. Now the lowsheet looks like Figure 1.10. Note that in the present
f
case, only the red arrows have been utilized.
.. . ,
-
a
0
.
>
-
Of.
.
<o-e-a.o.ir-
FIGURE 1.10
We can rename the stream(s) and block(s). To do that highlight the object we want
to rename and click the right mouse button. Select Rename Block and then give a new
name, as shown in Figure 1.11 for Block Bl.
a
-
«
r
0 %
' '
0 O-P-f c
.
a
FIGURE 1.11
Copynghied material
12 PROCESS SIMULATION AND CONTROL USING ASPEN
Alternatively, highlight the object, press Ctrl + M on the keyboard, change the
name, and finally hit Enter or OK. After renaming Stream 1 to F, Stream 2 to L,
Stream 3 to V and Block Bl to FLASH, the lowsheet inally resembles Figure 1.12.
f
~
- -
c-Q-
0
a-=
Si . , S
FIGURE 1.12
In order to inspect completeness for the entire process lowsheet, look at the status
f
indicator. If the message includes Flowsheet Not Complete, click on Material STREAMS.
If any red arrow(s) still exists in the lowsheet window, it indicates that the process is
f
not precisely connected with the stream(s). Then we need to try again for proper
connection. To ind out why the connectivity is not complete, hit the Next button on the
f
Data Browser toolbar. However, if we made a mistake and want to remove a stream
(or block) from the lowsheet, highlight it. right click on it. hit Delete Stream (or Delete
f
Anyway, suppose that the lowsheet connectivity is complete. Accordingly, the status
f
We have defined the unit operation to be simulated and set up the streams into
and out of the process. Next we need to enter the rest of the information using several
input forms required to complete the simulation. Within Aspen Plus, the easiest way to
ind the next step is to use one of the followings:
f
1 .
click the Next button
2 ind Next in the Tools menu
f
Copynghied material
INTRODUCTION AND STKPWISK ASPEN PLUS SIMULATION 13
r|nf? . .
l ..|..h nr .! -wi i - M
3
a-c
o
- m
(mu, im ml '
FIGURE 1.13
Configuring settings
As we click OiC on the message. Aspen Plus opens the Data Browser window containing
the Data Browser menu tree and Setup/Specifications/Global sheet.
Alternatively, clicking on Solver Settings and then choosing Setup /Specifications in
the left pane of the Data Browser window, we can also obtain this screen (Figure 1.14).
-JUS.'
; I*
.
. >i . ->
rr.Fi F
OQ-o-O-it-
FIGURE 1.14
14 PROCESS SIMUIvVTION AND CONTROL USING ASPEN
Although optional, it is a good practice to ill up the above form for our project giving
f
the Title (Flash Calculations) and keeping the other items unchanged (Figure 1.15).
. !
3af* I 3 ri-i - »i ji .1 H
-
."
y
*
-
(0-eo.o-1
FIGURE 1.15
In the next step (Figure 1.16), we may provide the Aspen Plus accounting information
(required at some installations). In this regard, a sample copy is given with the followings:
User name: AKJANA
Account number: 1
Project ID: ANYTHING
Project name: AS YOU WISH
\ r i-i i-f si
.
iO-Oo.Q.I.
m -
FIGURE 1.16
Copyrighted material
INTRODUCTION AND STEPWISE ASPEN PLUS SIMULATION 15
We may wish to have streams results summarized with mole ractions or some other basis
f
that is not set by default. For this, we can use the Report Options under Setup folder. In the
subsequent step, select Stream sheet and then choose Mole raction basis,
f
... - rJtW
. . --
g
'
""
t -
IZZi U-.-J7-- i
-
* '
i
*
.
(O-eo-e-T-
FIGURE 1.17
As illed out, the form shown in Figure 1.17, inal results related to all inlet and
f
product streams will be shown additionally in terms of mole raction. Remember that all
f
values in the inal results sheet should be given in the British unit as chosen it previously.
f
Specifying components
Clicking on Next button or double-clicking on Components in the column at the left side
and then selecting Specifications, we get the following opening screen (Figure 1.18).
iff i ijLJH .
.
(0-8-o.o.ir.
. * -
FIGURE 1.18
Copynghi
16 PROCESS SIMULATION AND CONTROL USING ASPEN1"
essentially an alias for a component It is enough to enter the formulas or names of the .
components as their IDs Based on these component IDs, Aspen Plus fills out the Type .
,
Component name and Formula columns But sometimes Aspen Plus does not find an .
exact match in its library. Like in the present simulation, we have the following screen ,
(Figure 1.19) after inserting chemical formulas of the components in the Component ID
column.
3513
I Toolt Run Plot Ltrarv . rxWv Help
"
_
L _
r i iL 3 I
-
--
i NoncOTvenbonal | Dalabanki |
Q Srxiaoor Oot-ons Dehne component t
Q StreanOas* Type Component name Fo.mula
.
Jj StXHtreans Convenhonal PflOPANE 3H9
- un*j-Set»
N'C MIO Cunvonbonal
N-C6HU Conventonal
© Spectficanona
'
I Assay/Bend
ught-End Preoert**
-
Jj P«ro Character natwn
Pieudocorrpcrientl
AfW-Corrpj
""
I Merry CiJmp*
"
l Ccrro-Coixrt
Propertes
StT«f"t
'
| Peacfcoro
Tind | EfcoWaaid j U e. Defied | Reade
* Conv Oofins
-
_
21 EOConvOpBora Componen ID. II daia arc 10 be foliieved Ifcim dalobanks. enici Shai Componert Hanw c FwmUa See Help
O sab*
j
"
-
Heai Etclianflet! | Coturr j Reacloit | Prenue Changers Manpiiaioo | Sobdi | U;ei Modets
"
Sep
-c. Sep2
FotMefcj preMFI . C:\ ..aFolde(5\A!penP1ui 11,
FIGURE 1.19
Obviously, only for Component ID C3H8, Aspen Plus provided the Component name
(PROPANE) and Formula (C3H8). This simulator does not recognize the other three
components by their IDs. Therefore, we have to search in the following way
(Figure 1.20) to obtain their names and formulas. Click on a component ID (say, N-C4H10),
then hit Find button.
Now, we have to give a hint with Component name or formula (butane) and then
hit Enter or Find now button (Figure 1.21). Apart from component name or formula,
we can also search a component by giving component class or molecular weight (range)
or boiling point (range) or CAS (Chemical Abstracts Service) number. Click on Advanced
button in the following screen to get these options.
INTRODUCTION AND STKPWISE ASPEN PLUS SIMULATION 17
i r- u.ivi»rT
»j <<||*i -| »| Qlral h»|
7 1
HmmaftmU
. / .'r -.r,,
Careen r
>i cm}: DMNMeMI
M«HU DBIWVWBM
*lf-Con«»
"
| M«fryCa«p«
it
_
j immmm
I COTvOcttins
-
tOCirr.0»ilicr»
3>!>!ic:- BBS wowre >:...<-» ..r .j -
irw, i«Hei
O M<4>
U k4
UmriKitrmx Sapaiataii | Han
Hm> Eichtnsan | COni
-
CH
m 6 o
mo. cm n
-
» st-t
FIGURE 1.20
MWflll III I I
f
I F--IHH>nr i i -i 'Ipi i m\ .i-i
i 35Ji
WMM
' Mirr? 1X414
' I'tM* |M.| fiiM
tOD O* *MM2
.
*
J 1 'ttrVM r tuhtn 114 2n -VJ J-' SrM O-' '
C4M»4 runrn
PURE 11
J
. I
Si K
FIGURE 1.21
18 PROCESS SIMULATION AND CONTROL USING ASPEN
- u m m: vr
f rai-«-l«»|««i|-4| *»-| »"l
,
.) i"! -I vj
"
ttl
I r- I " i I I M -leal : ! !
let >.-
Si - ~
j s- I
n tt-
FIGURE 1.22
option. Notice that if we make a mistake adding a component, right click on the row
and then hit Delete Row or Clear.
data for product streams L and V because those data are asked to compute in the present
problem (see Figure 1.24).
This property method assumes ideal behaviour for vapour as well as liquid phase.
C ll
INTRODUCTION AND STEI'WISK ASPEN PLUS SIMULATION 19
cina
Tiers r
"
i
0 samii (Ham
UVUM .
-
par-
-
AFU
r
I t4 -I - I . - |M
Co
f> . FBI
P j mi«D»
-
a-
HO-e-o-i-it. !
FIGURE 1.23
Ha 'ssH I
0]t*lMI _
rmr i~i-..t>-rv
f5~
Im«7V= 31
nns Dt
f, .rilll
, ri.ttn it:
-
'
I JIU-*"- I'M- . *
h o e czd- @ - it.
FIGURE 1.24
U3SE
Toob Ron Piol Lfciaiy Wmdow Help
~
did -J a M
D - I I 'I -isil I lai alS*l
~
i :r .- u>i"i-
gp-n=-3
O K>OE5I<iN
0 tMCPMAL
(#> TXPOftT -
i
O VIE
* .
ilj Advanced
JQ
-
Lifl >=
-
Input
EO varial
IS FLASH
| Be
Conv Op«noj
i EO Conv Option*
O Setup
DMOBasK
49 DMOAdv
Input CompteK
FIGURE 1.25
Now the Status message (Required Input Complete) implies that all necessary
information have been inserted adequately. Moreover, all the icons on the left are blue.
It reveals that all the menus are completely filled out. If any menu is still red, carefully
enter the required information to make it blue.
simulation, press OK on the message. We can also perform the simulation selecting
Run from the Run pulldown menu or using shortcut key F5 .
3
£=1
8 ! 7.1 CarrvOpllam
33 to Conv Option*
FIGURE 1.26
The Control Panel, as shown in Figure 1.27, shows the progress of the simulation.
It presents all warnings, errors, and status messages.
jNIRODUCTlON AND STEPW1SE ASPEN PLUS SIMULATION 4 21
,
QhrjAj*i j-j an .| ihi .j M
I 1"! _=J 3?) H -iroh L_jih
3 I
- "
3 r "
3 r
. loch:
Command Lr» |
AI bkK+» h«v» bean .
STREAMS
0
FU>»K3
6 -ciD
Fl<nH3 D«canl- Sup S»p2
FIGURE 1.27
Viewing results
Hitting ATex button and then clicking OK, the Run Status screen appears irst (see Figure 1.28).
f
yil l .i.l.lJIII«.II..IIHIII.I.IMItMIIIIH.HI II.Wl'ltlll.Ml.llltHW
-
I I Ffe Edt VKm Data Tools Rvxi Mot Lbtoty Window htetp
ItflHI -I v| daHal
-
3 m I _iJ_iMi_L
B Ru-i Slatut 3 sQg r
S i Streams
QU RaMiU Swranarv
-
Run Statu*
Streams
Convergence
(0-9 o 8 .
FIGURE 1.28
From the Data Browser, choose Results Summary /Streams and get the following
1 29).
screen that includes the inal results of the given problem (see Figure .
f
We can
Save the work by choosing File/Save As/...from the menu list on the top.
name the ile whatever we want. Note that an Aspen Plus Backup ile (*.bkp) takes
f
f
much less space than a normal Aspen Plus Documents ile ( .apw).
f
22 PROCESS SIMULATION AND CONTROL USING ASPEN
I I
» J/l Block*
£1 fo. * r
"
3 5l<»amT»blf[
r rjj 1 - 3
50 0 2000 200 0
i f. nr i ion oo ion on
f
Mas* Ftow b/h. 15906 41* 13312.698 2593 716
v. l:...- Flow culler 1039 561 382.439 3008 065
(V Mixw pMto! SoiMralof* { Heal Enchangon | Column* | Re»cto.. \ P-eume Chongeij \ Mo puMw* | So§*. | Um. Models )
STREAMS
HO 0 cd ' FVwh2
" ""
Flaih3 D>came> 3ep Sep2
For H*te, press Fi C.\- .g Pol<tou\Ajper. PK» 11,1 ! NUM i - .
.
.. Av,4,>:-
*
Start}} Aspen Pkn - Simulatl_
FIGURE 1.29
If we click on Stream Table button, the results table takes a place in the Process
Flowsheet Window, as shown in Figure 1.30.
Fie Edt View Data Tocfc Run Ffevaheet Librvy Whdow Help
F L V
EnlhJpy MMBtu/hi -
16583 -
12.499 -
2236
Hole Flw
Mok Trie
C3K8 0 100
.
0053 0 301
HX4HI0 0 200
.
0 169
.
0329
Mm/Spitlan Sflprntms { Heat Eicchangeit { Cdum | Reactori | PrMtue Chmgeii j Mmpdalai | Soldi j Use. Models j
D-»
<0-8-o
-
C:V.oF<*lefs\A»penMu»n.l
?1 1
FIGURE 1.30
INTRODUCTION AND STKPWISK ASPKN PI.US,M SIMUI.ATION 23
shown in Figure 1.30, along with the input summary in the inal report on the present
f
project (see Figure 1.31).
linput swimtrf crvaccd bv upen Plus "el. U.l tt tiiMtiS rrf jun a, 2007 ~
; Dlr»ctory CtSproarur 11 TBc\Aspanrai:n\WDrklng Polaei ' j' iveft Plus 11.1 tllnnm*
mMPuis
DPLUS RCSULTS-ON
IN-UNir» lii.
DGKRIPriON '
General Sl*u1al1e*< w<th English unl s :
F, dsI, Ib/hf, lEf«ol/»». oiu/hr, eirft/hr,
Propariy Haihooi wona
eln* M*l» for Incur: NOll
"
CJm8 C3h8 /
N-camo caxio-z /
N-cenW CftHH-l
"lOWSHEET
bicc> flash ih-e aut-v l
PROPERTIES SY5OP0
- - plash Flash;
kabah rtwp- ao. "sr.-ic-j .
FIGURE 1.31
Version/ Flash.rep). Subsequently, we may open the saved report ile (Flash.rep) going
f
through My Computer with using a program, such as the Microsoft Office Word or
WordPad or Notepad. A report ile for the present problem is opened below.
f
(i) FAR 52.227-14, Alt. Ill, (ii) FAR 52.227-19, (iii) DEARS
252.227-7013(c)(l)(ii), or (iv) THE ACCOMPANYING LICENSE AGREEMENT,
AS APPLICABLE. FOR PURPOSES OF THE FAR, THIS SOFTWARE SHALL BE DEEMED
TO BE "UNPUBLISHED" AND LICENSED WITH DISCLOSURE PROHIBITIONS.
CONTRACTOR/SUBCONTRACTOR: ASPEN TECHNOLOGY, INC. TEN CANAL PARK,
CAMBRIDGE, MA 02141.
TABLE OF CONTENTS
FLOWSHEET SECTION 2
FLOWSHEET CONNECTIVITY BY STREAMS 2
FLOWSHEET CONNECTIVITY BY BLOCKS 2
COMPUTATIONAL SEQUENCE 2
OVERALL FLOWSHEET BALANCE 2
STREAM SECTION 5
F L V 5
LOCATED IN:
PDF SIZE USED FOR INPUT TRANSLATION:
NUMBER OF FILE RECORDS (PSIZE) = 0
NUMBER OF IN-CORE RECORDS - 256
PSIZE NEEDED FOR SIMULATION - 256
DESCRIPTION
COMPUTATIONAL SEQUENCE
FLASH
TOTAL BALANCE
MOLE( LBMOL/HR) 220.462 220.462 0.000000E+00
MASS(LB/HR) 15906.4 15906.4 -0.782159E-11
ENTHALPY(BTU/HR) -0.165833E+08 -0.147349E+08-0.111463
COMPONENTS
***
MASS AND ENERGY BALANCE ***
INPUT DATA
***
RESULTS ***
N-C4H10 0 20000
.
0 16926
.
0 32874
.
1 9422
.
N-C5H12 0 30000
.
0 31602
.
0 23290
.
0 73697
.
N-C6H14 0 40000
.
0 46260
.
0 13781
.
0 29790
.
STREAM ID L V
FROM : FLASH FLASH
TO FLASH
SUBSTREAM: MIXED
PHASE: MIXED LIQUID VAPOR
COMPONENTS: LBMOL/HR
C3H8 22.0462 9 2754
.
12.7709
N-C4H10 44 0925
.
30.1237 13.9688
N-C5H12 66 1387
.
56.2424 9 8963
.
N-C6H14 88 1849
.
82.3291 5 8558
.
N-C4H10 0 2000
.
0 1693
.
0 3287
.
N-C5H12 0 3000
.
0 3160
.
0 2329
.
N-C6H14 0 4000
.
0 4626
.
0 1378
.
TOTAL FLOW:
LBMOL/HR 220.4623 177.9706 42.4917
LB/HR 1.5906+04 1 3313+04
. 2593.7158
CUFT/HR 1839.5613 382.4385 3008.0650
STATE VARIABLES:
TEMP F 50.0000 200.0000 200.0000
PRES PSI 15.0000 100.0000 100.0000
VFRAC 1.8002-02 0 0 .
1 0000
.
S FRAC 0.0 00 .
0 0
.
28 PROCESS SIMULATION AND CONTROL USING ASPEN1
ENTHALPY:
BTU/LBMOL -7.5221+04 -7.0232+04 -5.2612+04
BTU/LB -1042.5543 -938.9019 -861.9118
BTU/HR -1.6583+07 -1.2499+07 -2.2356+06
ENTROPY:
BTU/LBMOL-R -130.1235 -123.3349 -87.8846
1 4398
BTU/LB-R -1.8035 -1.6488 -
DENSITY:
LBMOL/CUFT 0.1198 0.4654 1.4126 02 -
BLOCK STATUS
**************
********************************************************
* *
*
Calculations were completed normally *
* *
*
All Unit Operation blocks were completed normally *
* *
*
All streams were lashed normally *
f
* *
************************************************************************:!:;!=
13
. .
3 Computation of Bubble Point Temperature
Problem statement
Compute the bubble point temperature at 18 bar of the following hydrocarbon mixture
(see Table 1.1) using the RK-Soave property method.
TABLE 1.1
c2 0 1
.
C3 0 15
.
i-Ci 0 1
.
n-Ci 02
.
i-C5 0 25
.
n-C5 0 15
.
Assume the mixture inlet temperature of 250C, pressure of 5 bar and low rate of
f
120 kmol/hr.
S,MULA' noN 29
Simulation approach
After starting the Aspen Plus simulator, the Aspen Plus Stnrt v i
Among the three choices, select Template option and then S
,.,
e F Tl 3
L L J.-i..'i- I iM BlMtt i ~| S!| -j j jj g j
H !i j
FIGURE 1.32
FIGURE 1.33
In the next ,
press OK in the Connect to Engine dialog. Once Aspen Plus connects to
the simulation engine, we are ready to begin entering the process system.
30 PROCESS SIMULATION AND CONTROL USING ASPEN
Creating flowsheet
Using the Flash2 separator available in the equipment Model Library, develop the
following process low diagram (see Figure 1.34) in the Flowsheet Window by connecting
f
the input and output streams with the lash drum. Recall that red arrows are required
f
ports and blue arrows are optional ports. To continue the simulation, we need to click
either on Next button or Solver Settings as discussed earlier. Note that whenever we
have doubts on what to do next, the simplest way is to click the Next button.
0
o
mm 1
o-e-oi-ir-
_
2£
FIGURE 1.34
Configuring settings
From the Data Browser, choose Setup I Specifications. The Title of the present problem
is given as 'Bubble Point Calculations'. Other items in the following sheet remain
untouched (see Figure 1.35). However, we can also change those items (e.g., Units of
measurement. Input mode, etc).
ij
, u m
it »
« "'E E3
FIGURE 1.35
INTKODUCTION AND STHPWISE ASFKN PLUSIM SIMULVTION 31
In the next, the Aspen Plus accounting information are given (see Figure 1.36).
_ .
.
r
-
t.-'l(.
.
I i Us*-**,
11 -
-
< O Q <=>. @ . 4 .
KM a IV- II I MM »»»»» !r.i-».
i
FIGURE 1.36
Specifying components
Click on TVex button or choose Components /Specifications in the list on the left. Then
define all components and obtain the following window (see Figure 1.37).
_
1 J
-
I-I"I> i j . i i xapji i iw 'H
tgj <«irM »i cuari n.i
r
/ -. ii.- MMM
9 BmiMM VMM
r -
.
L 11 d i it .jf
*
I
'
X . 4 i
.. .
l.w..-->., s-m
1
l
H 0 8 cd. g - f .
t
FIGURE 1.37
Copyrlqhted material
32 PROCESS SIMULATION AND CONTROL USING ASPEN
.
=1 3 JLi Si Mi bl
-
8
i 3;
F-3
. Q-S-o-'g-'ii
D
FIGURE 1.38
Hitting ATex/ button twice, we have the following picture (see Figure 1.39). The binary
parameters are tabulated below. When we close this window or cbck OK on the message.
it implies that we approve the parameter values. However, we have the opportunity to
edit or enter the parameter values in the table. In blank spaces of the table, zeros are
there. It does not reveal that the ideal mixture assumption is used because many
thermodynamic models predict non-ideal behaviour with parameter values of zero.
T£msxS\zi zl 2 '-I H 21 61
3 .ifLdB&teMMI)
n
m
TTD
-3
=w
:
3
MIX
« MM I *
FIGURE 1.39
INTRODUCTION AND STEPWISR ASI-KNJ>LU
sim 33
Specifying stream information
Click OK. Alternatively, use the Data Browser menu tree to navigate to the Streams/1/
Input/Specifications sheet. Then insert all specifications for Stream 1 as shown in Figure 1 40
J &1
J . 1 1,, I*
~
n 1 1 i 1 igi
la El »l aUl
JO
tcxnpojitior.
3 71 n
pr
(5
1 tdvaoced
f
rr s;; -
J
~
r i Rpioftt
& Setup
L55s?P Bos-:
[i hWs/Sphleu Ssp falais j He Esdw ers j Columns | Reaclw: [ Piessine Changers j Manipulators : Solids j UferModefi
Matenal
FIGURE 1.40
form, enter the required inputs (Pressure = 18 bar and Vapour fraction = 0) for block
BUBBLE (see Figure 1.41).
-
al>l
3 *i I *! «iEi
"
si - r
/Speatifotnni FlathOcdoni ' Er**rrr«nl
i Pr<*8rUw
CJ 2
~J 1
-
1
.
_
)
9 ***
-
a-»
STREAMS
0 0 Qfl
Ri ? fWJ Ete i S<c S Bi
FIGURE 1.41
34 PROCESS SIMULATION AND CONTROL USING ASPEN
Running the simulation
Press Next button and then hit OK to run the simulation The following Control Panel .
laillUJLIIIlllBIl
i f** t-t Vwi- Data Toofi Bun Lfciafy WVidwv/ Meto
JDlugB] atfij J-j ral-rlatl-l-qi l «ii [=5 _
d .
IjjJ -J SI afl e-i
f
i 3 - i i -i tim mi
41 - I "I JiJ S i l <| jjB] | IMI
ss 3 r i r
j
i - r-. - and Po:iiotve SckXi
"
Command Lr* f
| Sepaislota | Heat ExchM oeic [ Cohjnmi | Rrac'txt | Prectuie Changed SoMt j Uiei Model: |
STREAMS
'
o o
flathZ Flail-.3 Decade.
e - it -
Sep Sepg
C:\...oFoW«"\A«>enPtuj 1
FIGURE 1.42
Viewing results
Clearly, Figure 1.42 includes the Status message: Results Available. As the simulation
calculations completed, click on Solver Settings and then double-chck on Blocks to obtain
the following screen (see Figure 1.43).
-
-
"
J _
l l-l-PT
-
j j b I 3 ±tti iiJfXi 32>J jaal n i
33 Set-*
O Setup
i SOU. | U-Moa* |
-
CM
<o 9 e i
1 n A?
. FWi3 Dac«4«> Sap -
| ,Md«fl« .»-.. .
FIGURE 1.43
INTKODUCTION AND STKPWISK ASI'KN PLUS SIMULATION 35
JaflHI Ml *1
mi
WMllwilfc ii»»i»y
NM1 »»»
IB 3(v«««iP»*Jl
ra f O Cor- OBban*
V
O 1
tmt
| SOU. | UnMaM |
.
j 0 - 6 -o- f. r.
111,1 ' MM *r
FIGURE 1.44
From the results sheet, we obtain the bubble point temperature = 42.75411960C.
13
. 4
. Computation of Dew Point Temperature
Problem statement
Compute the dew point temperature at 1.5 bar of the hydrocarbon mixture, shown in
Table 1.2, using the RK-Soavc property method.
TABLE 1.2
C2 0 1
.
Ca 0 15
.
<-c4 0 1
.
n-CA 0 2
.
M3a 0 25
.
"-
(>
,
0 15
.
Assume the mixture inlet temperature of 250C, pressure of 5 bar and low rate of
f
120 kmol/hr.
36 PROCESS SIMULATION AND CONTROL USING ASPEN
Simulation approach
As we start Aspen Plus rom the Start menu or by double-clicking the Aspen Plus icon on
f
our desktop theAspe?i Plus Startup dialog appears (see Figure 1.45). Select Template option
,
.
Id *J 1PJ M _j.
l LJ...l-:i.::.l JAI
_.
"
-
I/I I J J_J_J_:J..J gj
-
J
3!
1 empWis
ft? Start] j £
FIGURE 1.45
As Aspen Plus presents the window after clicking OK as shown Figure 1.45, choose
General with Metric Units. Then press OK (see Figure 1.46).
MB
i i i
Flowsheei
Daneel
FIGURE 1.46
INTRODUCTION AND STEPWISE ASPEN PLUS SIMULATION 37
Subsequently, dick OK when the Aspen Plus engine window pops up.
Creating flowsheet
In the next, we obtain a blank Process Flowsheet Window. Then we start to develop
the process lowsheet by adding the Flash2 separator from the Model Library
f
toolbar and joining the inlet and product streams by the help of Material STREAMS
(Figure 1.47).
-
H c-
0 St*
-
CD c
± if
FIGURE 1.47
Now the process low diagram is complete. The Status bar in the bottom right of
f
the above window (see Figure 1.47) reveals Required Input Incomplete indicating that
input data are required to continue the simulation.
Configuring settings
Hitting Next button and then clicking OK, we get the setup input form. The present
'
problem is titled as Dew Point Calculations' (see Figure 1.48).
In Figure 1.49, the Aspen Plus accounting information are provided.
Specifying components
Here we have to enter all the components we are using in the simulation. In the list on
the left, choose Components /Specifications and ill up the table following the procedure
f
Copyrighted malenal
38 PROCESS SIMULATION AND CONTROL USING ASPEN1
,J9J »i
J U
- LT _l_LJ__rv 3 I I _lL
v ...
S«up
</Gh*»l| OeKtw»ion | ilccojnlh.0 ] DiagrwHiei |
© SprtlfttaMoo,
© 5(re*T,GI«S
Urn'-. d( rfi?B;ijHMr*rJ
Ur«j-S«ti
MET I f - 'j 3
Oi>pm (en Inc ji mode 3
.
»J
.
J J 3
Blocks
-
I
CorVrwoente
L
STREAMS S
0 o 8 I Scp2
tilCTh2 Flath3 Decaniei Sap _
FIGURE 1.48
Mi
Fie E* View Data Tools Run Plot Lfcrary Window Help
r
|0 Specfeno
h j/j Setup /GlobalI ./Descnption /Arc Diapnoslici
Specifications
SiroJatron Options : Aspen Plus accounitng ioifnation
r
bfe Subsbeams
j U set name: jAKJANA
Accouil luffnber:
S 1 3 Units-Sets
Q Custom Unrts I PtqedID ANYTHING
0 Report Options j Project name: jAS YDU UKEi
*:
| Components
+ Properties
I Streams
_
iJ Blocks
Si Reactions
+ Convergence
+ FtowsheeSng Options
* Mate'Analysts Tools
*
.
ifl Cor/iguration
Q] RtsJts Sunmary
Project n .
Input Canvtete
pT MMea/SpiKm SefMHtais | HutEictangen | Cokmu | Reactat | Pteuute Changeu j ManpuWcs | Stfcc | Us« Modeh
f
H 0 -0 -0
C < Flash? FImKI DmcmIm
FIGURE 1.49
INTRODUCTION AND STEPWISE ASPEN PLUS SIMULATION 39
i--rr - i-!>-i' it
» «
.
Jj I 111
t r
-
Q*
.
<o-8o.i.ii:
FIGURE 1.50
=1
IMXooc Cfcc* »i c*iiBMr opto* SHI
r . -.
-I
i
-
1 .i
4 'P- .
FIGURE 1.51
Copyrighted material
40 PROCESS SIMULATION AND CONTROL USING ASPEN
Entering all required information one obtains the screen as shown in Figure 1.52 , .
II
-
WO[>l ili
<:< 'i
ptseu i t«um [Met. »]:
fin f hr*"" .] ik
BrTH>l
. its til
_
j Cfctt/ » runt
, .~J>
_
a .
Micw illii
"
-
jO 8 0-8 C
* «-«| I AvenPfcn-SknutotL-
FIGURE 1.52
on Blocks and then select DEW. Filling up the input form, we have Figure 1.53.
Btali
r
: ,. uh.|! rv I 3 /-\ 1 -i -Igil I Ji
3
"
3 FlMhOpooru I Em
30
3n
" -
"
ANWIJ-1
ANCMIJ-I
3
BnSBIJ-1
0XTrI>l
'
I Group
Input CcntfM*
I*
STREAMS
<0 Q cd
FIGURE 1.53
INTRODUCTION AND STEPWISK ASPEN PLUS 3!MUU\TION 41
-
j r-H'-hrr
- II . t .! -1031 I - ! !
HI 3
3 ,
mi r . .1:
IIKMM
D -
o-e-oi-it-
-I*. MM'I (Will
FIGURE 1.54
Viewing results
First click on Solver Settings. From the Data Browser, choose Blocks/DEW/Results
(see Figure 1.55) to get the dew point temperature = 22.19453840C.
i' r ui>.i.rf
-
i* I- *
.MM
-
I
a -.
MM*
MM*
JVM -
hum
.
(O-e-o-i-it-
1 im*f n u t- w ' i
FIGURE 1.55
42 PROCESS SIMULATION AND CONTROL USING ASPEN
13
. .
5 T-xy and P-xy Diagrams of a Binary Mixture
Problem statement
into a lash chamber (Flash2) with a flow rate of 120 kmol/hr at 3 bar and 250C
f
Simulation approach
As usual, start Aspen Plus and select Template. Click OK to get the next screen and
choose General with Metric Units. Then again hit OK. In the subsequent step, click OK
in the Connect to Engine window to obtain a blank Process Flowsheet Window.
Creating flowsheet
From the equipment Model Library at the bottom of the Aspen Plus process lowsheet
f
window, select the Separators tab and insert the Flash2 separator. Then connect the
separation unit with the incoming and outgoing streams. The complete process is shown
in Figure 1.56.
CD o
0
-
o
STfSAMS
9-o
FIGURE 1.56
Configuring settings
After clicking on Solver Settings, select Setup /Specifications in the list on the left. The
Title of the present problem is given as TXY and PXY Diagrams'. Subsequently, the
'
Aspen Plus accounting information are also provided [see Figures 1.57(a) and (b)].
INTKOIHTTION AND STKI'WISK ASl'liN I'l.l'S ' SIMULATION 43
r 1 1 i .1 ..b; : . IP «|J»I
-
'-
: 5! r-
CO . j .
FIGURE 1.57(a)
213
I -IE' - ;« 1
1
gjg j i 3 <iF-i] >>j aal -I
?3
3 <o'»-
-
CH
oeo.@.
«t Mi(.i,tra5»«-i
FIGURE 1.57(b)
Specifying components
Hitting Next button and defining the components (ethanol and water) in the input
form, one obtains Figure 1.58.
For this problem, choose the Wilson model by scrolling down (see Figure 1.59).
44 PROCESS SIMULATION AND CONTROL USING ASPEN1
iL !
r j .Lj.inr
__
L
- - , _
l L
1© Specfrahon:
'
v
Selertionj PtMeum | NoneonvtrMnal | Oat nki {
ComponnHID
_
J Sot rr s
iTHANOl
© SpetKkattww
| iitM) .WwiO
'
lltr-ComOS
f
'
) HerryCofroi
.
_
J P nwt
Strrars
21 EC Ccti. 0«JOftt
_
1J
fi< MotB Ephieri Sopaialod j HetfEfihsrflw; | Colutnu | Rwcfct j PremMChanel | ManpJsiai j SoHi | UteiModeU |
f
Matewl
STREAMS 1 lwW Sep
f
FIGURE 1.58
"
3
"
3M t r
© 5wj(aC»i Opt mm Gtobalj Flow!heel Seciooi j Heleie.-ced j
Streain Class
-
j Substra4fn» ' c . meihodi i moddi
w-y rr lhod fwiUSDN
_
*
Ljfl Un*s-S«s P.oc«i two P
Vl
©
Components
Report OcKnns Base method
Henry component
| WILSON
i 1 Modty piopeity models
1 Asia//Blend
=3
~
"
3
.
_
lj petfo Cha< actof cation Wale, solubtty j 3
P ss-jdoeoinporients
ACw-Comp* Electic le cal ula'lor.oplnni
| I Menry Comp* "
3
UWIFAC (Vo
__
j Comp-Ofoupi
.
_
y ConxHJsts
Ptopertes
SpeclTKatton*
.
_
V) Pf«5<rrty Method*
* Estmalion
.
_
J Molecular Structure
» Paramateri
'
i Data
*
33 Prop-Set)
. Advanced
Input Complete
it' MMKt/5p(tl«M S«pai«foi> j Heat Exchange" | Cokjmnl | Reactors | Piohuio Chonoelj i ManipuWoit | SoWi | LUei Model: |
f
FIGURE 1.59
Once the base property method has been selected and we click the Next button, a
window pops up asking whether to continue to the next step or to modify the properties
(see Figure 1.60).
INTRODUCTION AND STEPWiSK ASPEN PLUS SIMULATION 45
OK Cancel
FIGURE 1.60
f
rate and components mole fraction, one obtains Figure 1.61 as shown.
J6
I -
.
- --r- :
it
3
m
If -
."HI
ft -nvm i iia"
m m» wa«
nr
-
a -
o-e-o-i-ir-
FIGURE 1.61
Copyrighled material
46 PROCESS SIMULATION AND CONTROL USING ASPEN
ffTIwioL
F3
J BSg I
P jT H il 3 s*»h
FIGURE 1.62
We must note that this option can be used to generate T-xy, P-xy or Gibbs energy of
'
mixing diagrams. Select Txy for the present problem. We aim to do an analysis on the
mixture of ethanol and water; so select these components accordingly. The user has the
'
option of specifying, which component will be used for the x-axis (which component s
mole raction will be diagrammed). The default is whichever component is indicated as
f
component 1. Make sure that we are creating the diagram for the mole raction of ethanol.
Entering required information Figure 1.62 takes the following form (see Figure 1.63).
, f
[WATER
His-1
|l
FIGURE 1.63
Click on Go and get the T-xy plot at a constant pressure (1 013 . bar) as shown in
Figure 1.64. Although the Status bar shows Required Input Incomplete ,
but there is no
problem to get the plot based on the given information.
INTRODUCTION AND STEPWISE ASPEN PLUS1 M SIMUL.\TION 47
FIGURE 1.64
It should be noted that if we move the T-xy plot slightly or close it, we ind
f
Figure 1.65 having a databank. Some of these values have been used to make the
plot (Figure 1.64).
n3K|fc!»|-qM!!H 3i
-
i i i .i .m- i m\
Mil* I MMI
rwrrc
Wtfm
: : -
<o e-o-e-a;
«Mp»iM.i lM«4Ml-
FIGURE 1.65
Copyrighted material
46 PROCESS SIMULATION AND CONTROL USING ASPEN
.
r
"
-
:;
- i Mf n: j i i 'i mo; i i"i SM
. .
REAMS
j tRLAMb
f-
FIGURE 1.62
We must note that this option can be used to generate T-xy, P-xy or Gibbs energy of
mixing diagrams. Select Txy' for the present problem. We aim to do an analysis on the
'
mixture of ethanol and water; so select these components accordingly. The user has the
'
option of specifying, which component will be used for the x-axis (which component s
mole fraction will be diagrammed). The default is whichever component is indicated as
component 1. Make sure that we are creating the diagram for the mole fraction of ethanol.
Entering required information. Figure 1.62 takes the following form (see Figure 1.63).
M il SI M SI
2 |W*TER 3
fETHANOL
3
(\oflm<y cmms
-
o-
FIGURE 1.63
Click on Go and get the T-xy plot at a constant pressure (1 013 . bar) as shown in
Figure 1.64. Although the Status bar shows Required Input Incomplete ,
but there is no
problem to get the plot based on the given information.
INTRODUCTION AND STEPWISE ASPEN PLUS SIMULATION 49
Clicking on Go button, we have the following P-xy plot |see Figure 1.68(a)|
at a constant temperature (90oC) and respective databank produced (Figure 1.68(b)|.
I I- I-I.-IkU-
L *
-
i
1 1 1 DM
- .
. ..
.
11
FIGURE 1.68(a)
B I ...
11 n?
TWB ! 11 sarm
A ij " nBsw
"
Wii
"
B .
'
ams IWB T7ZB5 ITFii
r
moo n
B [B TIPl fSiiTE IfifiBT
f
« nw na! 61 1
.
HD-*
HO-Q-o-i-iI-
Ikvil .H«i) Ohm Sh> !»!
IU MM
FIGURE 1.68(b)
Copyrighled material
50 PROCESS SIMULATION AND CONTROL USING ASPEN
Notice that the plot window can be edited by ight clicking on that window and
r
selecting Properties In the properties window, the user can modify the title axis scale
.
,
,
font colour of the plot, etc. Alternatively, double-click on the different elements of the
,
plot and modify them as we like to improve the presentation and clarity .
In this chapter a brief introduction of the Aspen simulator is presented first. It is well
,
recognized that the Aspen software is an extremely powerful simulation tool in which , ,
a large number of parameter values are stored in the databank and the calculations are
pre-programmed. At the preliminary stage of this software course, this chapter may
help to accustom with several items and stepwise simulation procedures. Here four ,
simple problems (flash calculation, bubble point calculation, dew point calculation and
T-xy as well as P-xy plot generation) have been solved showing all simulation steps .
PROBLEMS |
1 . 1 A liquid mixture, consisting of 60 mole% benzene and 40 mole% toluene, is fed
with a low rate of 100 kmol/hr at 3 bar and 250C to a flash chamber (Flash2)
f
operated at 1.2 atm and 100oC Applying the SYSOP0 method, compute the
.
10 mole% o-xylene, is lashed at 1 atm and 110oC. The feed mixture with a low
f
f
rate of 100 kmol/hr enters the lash drum (Flash2) at 1 atm and 80oC Using the .
f
TABLE 1.3
i-C4 12
n -C4(LK) 448
i-C 5(HK) 36
Ce 23
C7 39.1
272.2
c9 31
876.3
The lash chamber (Flash2) operates at 180oF and 80 psia. Applying the SYSOP0
f
1 . 4 Find the bubble point and dew point temperatures of a mixture of 0.4 mole fraction
toluene and 0.6 mole fraction rso-butanol at 101.3 kPa. Assume ideal mixture
and inlet temperature of 50oC, pressure of 1.5 atm, and low rate of 100 kmol/hr.
f
1 5 Find the bubble point and dew point temperatures and corresponding vapour
.
f
100 kmol/hr enters at 50oC and 1 atm. Consider ideality in both liquid and vapour
phases.
1 6 Compute the bubble point and dew point temperatures of a solution of
.
n-C4 0 25
.
n-C5 04
.
Ce 03
.
The ideal solution with a low rate of 100 kmol/hr enters at 50oC and 1 atm.
f
1 . 7 Calculate the bubble point pressure at 40oC of the following hydrocarbon stream
(see Table 1.5).
TABLE 1.6
c2 0 1 .
Ca 0 15 .
i-C4 0 1 .
n-C4 02 .
i-Cs 0 15 .
n-C5 0 15 .
c6 0 1 .
Use the SRK thermodynamic model and consider the inlet temperature of 30oC,
pressure of 4.5 bar and low rate of 100 kmol/hr.
f
1 .
8 A binary mixture, consisting of 50 mole% ethanol and 50 mole% 1-propanol, is
fed to a lash drum (Flash2) with a low rate of 120 kmol/hr at 3.5 bar and 30oC.
f
into a decanter model run at 341.1 K and 308.9 kPa. To identify the second
liquid phase, select n-pentane as a key component (see Table 1.6).
52 PROCESS SIMULATION AND CONTROL USING ASPEN
TABLE 1.6
Applying the NRTL property method simulate the decanter block to compute
,
1 .
10 A ternary mixture having the following low rates is fed to a separator (Sep2) at
f
50oC and 5 bar (see Table 1.7).
TABLE 1.7
n -pentane 33.623
ethanol 0 476
.
water 3 705
.
To solve the present problem using Aspen Plus, the following specifications are
provided along with a T/F ratio of 0.905478 (see Table 1.8 and Figure 1.69).
TABLE 1.8
ethanol 09 .
B -O
Applying the SRK property method, simulate the lowsheet, shown in Figure 1.69,
f
SOLIDS base property method and simulate the dryer model (Flash2) to compute
the recovery of water in the top product.
INTRODUCTION AND STKPWISE ASPEN PLUS SIMULATION 53
Wet
Temperature = 75DC
Pressure = 1 aim AiROur;
Flow rates
S(02 = 800 Ib/hr
H20 = 5 Ib/hr WET
Air
Temperature = 200oC
AIR
0
Pressure = 1 atm
f
REFERENCE
2 1
.
BUILT-IN REACTOR MODELS
In the Aspen Plus model library, seven built-in reactor models are available. They
are RStoic, RYield, REquil, RGibbs, RCSTR, RPlug and RBatch. The stoichiometric
reactor, RStoic, is used when the stoichiometry is known but the reaction kinetics is
either unknown or unimportant. The yield reactor, RYield, is employed in those cases
where both the reactions-kinetics and stoichiometry-are unknown but the product
yields Eire known to us. For single-phase chemical equilibrium or simultaneous phase
and chemical equilibrium calculations, we choose either REquil or RGibbs. REquil model
solves stoichiometric chemical and phase equilibrium equations. On the other hand,
RGibbs solves its model by minimizing Gibbs free energy, subject to atom balance
constraints. RCSTR, RPlug and RBatch are rigorous models of continuous stirred tank
reactor (CSTR), plug low reactor (PER) and batch (or semi-batch) reactor respectively.
,
f
Eor these three reactor models, kinetics is known. RPlug and RBatch handle rate-
based kinetic reactions, whereas RCSTR simultaneously handles equilibrium and rate-
based reactions. It should be noted that the rigorous models in Aspen Plus can use
built-in Power law or Langmuir-Hinshelwood-Hougen-Watson (LHHW) or user defined
kinetics. The user can define the reaction kinetics in Fortran subroutine or in excel
worksheet.
One of the most important things to remember when using a computer simulation
program, in any application, is that incorrect input data or programming can lead to
solutions that are correct based on the program's specifications but unrealistic with
" "
regard to real-life applications. For this reason, a good knowledge is must on the reaction
engineering. In the following, we will simulate several reactor models using the Aspen
Plus software package. Apart from these solved examples, interested reader may
simulate the reactor models given in the exercise at the end of this chapter.
54
ASPEN PLUS SIMULATION OF REACTOR MODELS 55
22
.
ASPEN PLUS SIMULATION OF A RStolc MODEL
Problem statement
Pure ethylbenzene enters the RStoic reactor with a low rate of 100 kmol/hr at 260oC
f
and 1.5 bar. The reactor operates at 250oC and 1.2 bar. We can use the fractional
conversion of ethylbenzene equals 0.8. Using the Peng-Robinson thermodynamic method,
simulate the reactor model.
Simulation approach
As we start Aspen Plus from the Start menu or by double-clicking the Aspen Plus icon
on our desktop, first the Aspen Plus Startup dialog appears (see Figure 2.1). Choose
Template option and then click OK.
_
J
FIGURE 2.1
As the next window pops up (see Figure 2.2), select General with Metric Units and
hit OK button.
Copyrighted materia
56 4 PROCESS SIMULATION AND CONTROL USING ASPEN
jzj
I M I I I lAl I I I- I
.
j .j--jc-r; ] f.-S- -.r 3 C j n-V; j «' if!: VV.
FIGURE 2.2
Here we use the simulation engine at 'Local PC. Click OK when the Connect to
Engine dialog is displayed (see Figure 2.3). Note that this step is specific to the installation .
Connect to Engine
User Info
Node name :
User name:
Password:
Working directory:
OK Exit Help
FIGURE 2.3
Creating flowsheet
We are now ready to develop the process low diagram. Select the Reactors tab from
f
the Model Library toolbar, then choose RStoic icon and inally place this unit in the
f
blank Process Flowsheet Window. In order to connect the feed and efluent streams
f
MODELS 57
two red arrows and one blue arrow. Remember that red aarr0WSfare 're(luired
f
rrow
blue arrows are optional ports. ts and
Click once on the starting point, expand the feed li ne and click a~Hn tv f a
stream is labelled as 1. Addmg the outlet stream to the reactort tJXwa WW
-
n
y' UIiaiiy
we make the image as shown in Figure 2.4.
I .lal I Ml
= 03-
-
Q a
In
i . i . S -O-M-i
o
-
a Ri astt. tb pfvjj
FIGURE 2.4
» l«IVl
-
'
. r- C«J Kf! Pin tr-Kl«- LI'-TV iWoc,-. i
'
c*
f
Eh-
at rsms acs'R
FIGURE 2.5
Obviously, the Sia s md/cator in the bottom right of the mam window h changed
the message rom Flowsheet Not Complete to Required Input /ncom ff . fsimulation.
lete the ation
to enter th* remaining data using several input forms required to comp
f
58 PROCESS SIMULATION AND CONTROL USING ASPEN
Configuring settings
Hitting Next icon and clicking OK on the message sheet displayed we get the setup input ,
form. First the title of the present problem is given as 'Simulation of the RStoic Reactor'
In the next, the Aspen Plus accounting information (required at some installations)
are provided.
User name: AKJANA
Account number: 5
Project ID: ANYTHING
Project name: YOUR CHOICE
Finally, select Report Options under Setup folder choose 'Mole' as well as 'Mass' ,
fraction item under Stream tab (see Figure 2.6(a) (b) and (c)). ,
MM ±S
_
i r- i - i- i jv
_
-
i « i iai
UMsi
[jjttiEjjft
L- - .1
J .
lU -
I- S . S . § -Q-M-O-
B.BM Bi u.
.
'-.C--- KC TIi PFtjj Rfem.
FIGURE 2.6(a)
FIGURE 2.6(b)
ASPEN PLUS SIMULATION OF REACTOR MODEI S 59
Mil
: r-i-hi r» ! .|gi
.
i ip' h-i
,
Dm dm
r _ utM
-
it
I i i I M>l Umomm I
-
O
-
'
tifc waw «
f
FIGURE 2.6(c)
Specifying components
In the Data Browser window, choose Components /Specifications to obtain the component
input form. Now ill out the table for three components, ethylbenzene, styrene and
f
hydrogen (see Figure 2.7). If Aspen Plus does not recognize the components by their
IDs as defined by the user, use the Find button to search them. Select the components
from the lists and then Add them. A detailed procedure is presented in Chapter 1.
I?! 1 -1-
i "" TH III
sr-l© 8 18 0IIU
FIGURE 2.7
fd materic
60 PROCESS SIMULATION AND CONTROL USING ASPEN
property input form. Use the Peng-Robinson thermodynamic package by selecting PENG -
ol lBj _
J_ J w] KW«>|<H m -1 H JpJjJ J
r
"
3
""
"
3
J * ew Proc*li«t 1 3
f
3
"
U **-
3M| #
'
FIGURE 2.8
DZSMSSEGSSSD
: i* Hot Utorr Wrdo* H-fc
f
A|>Mdiedtio<n| FWiOwonTl
|/M1XED
i]
o'***
t
"
Ware 2mu«
3
~
"
3
p3" 3
(i a ****** 3
1 'wr-S-
0
D
1
J -
' ur
ttBUH,
m
eS
'
i q
-1 - Bli. - BO f""
su
SSiS :
r
. ,.JCT»J .
FIGURE 2.9
ASPEN PLUS SIMULATION OF REACTOR MODELS 61
information
Specifying block
From the Data Browser, select Blocks/REACTOR. Specifying operating conditions for
the reactor model, the form looks like Figure 2.10.
FIGURE 2.10
to choose the reactants and products using the dropdown list input the stoichiometric ,
coefBcients and specify the fractional conversion In the Aspen Plus simulator, coefficients .
should be negative for reactants and positive for products (see Figure 2 11). .
** b* bo "e*
> '-'
RiACTQR
FIGURE 2.11
62 PROCESS SIMULATION AND CONTROL USING ASPEN
required input information have been inserted by the user. There are a few ways to
run the simulation. We could select either the Next button in the toolbar which will tell
us that all of the required inputs are complete and ask if we would like to run the
simulation. We can also run the simulation by selecting the Run button in the toolbar
(this is the button with a block arrow pointing to the right). Alternatively, we can go to
Run on the menu bar and select 'Run' (F5).
MM.|8W«'!i ,l|Hllir
DMll
I M ill
Elfb Imeicbah
""
I Comp-GroLps
' Con-.p-Lis's
*
1 1 Cperty Methods
S tstrfi ficn At tequfed npd u ocmpHe Y j can rui the MnuMlon nitw. w
t
':
Jj Moiecua- Sbuctm iiu can erttr more input To er4er more f-pj. Bated Cared th«n
> p ParameJers seled t e ooUont yoj mM tnyn Ihe Dais poldOAT-, menu
D a Rui ir-e sirxilatiwi now?
S Advanced
&reanS
_
Jfl :
Inpu C«nplete
[H " Mnwii/SpWer; Hea<Ev.-. i9Pt; J.,, Rucloi Chsnga, | M<n>M>t ( 5c«> j U»Mo*b |
-
CH
'
STREAMS RStdc RYieW
F-
r H«o press F1
, ,
Stall *
'
Aww.RaocDdr | « Awr.Mcd I
-
.
Boot
_
FIGURE 2.12
Viewing results
As we click OK on the above message the Control Panel appears showing the progress ,
of the simulation. After the simulation is run and converged we notice that the Results ,
Summary tab on the Data Browser window has a blue checkmark Clicking on that tab .
will open up the Run Status. If the simulation has converged it should state ,
"
Calculations were completed normally" (see Figure 2 13). .
Pressing Next button and then OK, we get the Run Status screen In the subsequent .
step, select Results Summary /Streams in the list on the left and obtain the final results
(see Figure 2.14). Save the work done by choosing File/Save As/...in the menu list on
the top.
If we click on Stream Table knob just above the results table, the results are recorded
in the Process Flowsheet Window, as shown in Figure 2.15
ASPEN PLUS SIMULATION OF REACTOR MODELS 63
'
>k r« An [Mi Tot \r Um <.>«. MO
- - -
-
CH
s ' w*. mw< n* w> ii>*j
FIGURE 2.13
-
T I M -I -lei
I I
- "
3» 3 "-'-I
J
-
a)
55T:
-
t« " <w
"
is*" .
, i i.
M», iMm
l>ii l«U
Man'*
m
cna
inpp BUB
incMA
m- @ . i . e u m u
' irw« wto* MMt ac» S
"- I r
FIGURE 2.14
64 PROCESS SIMULATION AND CONTROL USING ASPEN
Ffc Edt Vfew 0«a Into tjn nowheri Lfc fy VAk w Het
f
'
-
lup I , IT _
LiiE| | |a|
STREAMS 1 BSinc
, S 0 U 31U
BEoii HGMis RCSTB BPItg RBaM.
\ s FoWen JJswn Ru» H 1 HUM lfloAi Artahie
f
FIGURE 2.15
CBSES
Fie £* Forw* >Atw
input Sugary created by Aspen Plus K«1. 11.1 at 12:U:CM Thu jul 5, 300?
Oirecrory C: Proqr-5R Pi les'AspenTech .norfcing Pol ders'.Aspen Plus 11.1 Fllep
'
"
&
f
'
fASS-DENSITV-
'
kg.'CUH" W)LE-£NTHALP- kcal,'noV t,
t-ASS-EWTM .P- kcal/kg' HE*T-MMkcal t'OLE-CONC-'mol .
T
'
&
POBOP-bar
DESCRIPTIOH "
General SlHllailoi) mith Metric units :
C, bar, kg/hr, knclhr. MMKcal/hr, c\m/hr.
property Method: Mone
PBOPERTIFS PENG-ROB
5THCAH F
SUBSTBEAf KIXEO TCHP-J60. PRE5-1.S MOLE-FLOW-100.
W>LE-FMC ETHYL-01 1.
'
B i I vjnwi-* |- la»«Jtol | lto.»,-s || -WEME1 :« jpCittU
FIGURE 2.16
- y wkjusu,jO t
f 65
DO
f
instructions as presented in Chapter 1 P eSent Problem, follow the
f
.
2.
3 ASPEN PLUS SIMULATION OF A RCSTR MODEL
Problem statement
f
the liquid-phase reaction, the inlet streams have the specifications shown in Table 2 1
, . .
TABLE 2.1
Fake reaction kinetics data for the Arrhenius law are given as:
Pre-exponential factor = 5 x 105 m3/kmol s
Activation energy = 20,000 Btu/lbmol
[CJ basis = Molarity
Use the SYSOP0 base property method in the simulation. The reaction is irst-order in f
aniline and hydrogen. The reaction rate constant is defined with respect to aniline.
Simulate the CSTR model and compute the component mole fractions in both the liquid
as well as vapour product.
Simulation approach
Start with the General with Metric Units Template, as shown in Figures 2.17(a) and (b).
Click OK in the above screen. When the Connect to Engine dialog appears, again
hit OK knob to obtain a blank Process Flowsheet Window.
Creating flowsheet
Select the Reactors tab from the Model Litwy
models available. Among them, choose RCSTR P cetheit process
t
in tn low magr
them,
Adding inlet and product streams and renaming
f
Q|a|B| JJ J_J
_
nMfel I 1 :1 si 21 __1_L.J ni M M ®l
A1 ] c 8lor+. SmuWen
r OMUsnE.ulr.lSim.jl.j'i-
"
FIGURE 2.17(a)
MMtajJ-V arvtr
nitpwi wi mmi ,
SpNtft/Chmic*
SUrti
'
FIGURE 2.17(b)
ASPEN PLUS1" SIMUIATION OF REACTOR MODELS 67
o|rf|y|giai _
lg|g rj|twi<H;j 3 I ("l jajgj
-
I}-
FIGURE 2.18
Configuring settings
Hit Afe button and then OK and get the setup input form. The present project is titled
as '
Simulation of the RCSTR Reactor' and the accounting information are given as
'
- II -». w . t -.-.
f
FIGURE 2.19(a)
68 PROCESS SIMULATION AND CONTROL USING ASPEN1 M
0 Spiicfcii
l li«< MBW
.
.
jfl IM-SHi
O CuHsflUnli Rovci ID
f
kfUCdRfMi
STREAMS
'
RSioc RYwId
. 11 -y-Bo
REgnl RGMw RCSTfl RWjg REafch
O * $3 17 1'.
FIGURE 2.19(b)
In the subsequent step, choose Setup/Report Options / Stream rom the Data Browser f
'
window and select Mole' as well as 'Mass' fraction basis (see Figure 2.20).
.
Jfl Ml S«t»
P MtJa P Mcta TFF [gENJ T]
! r Uau P Mm |S Standard fa0cdm>i
P S.>- .:abh t
P Componerti t h (wo o-. 01 H itDon
-
t
StfltW
(Bill
BV fvuc RE.M- RGte. RCS1R RPI m j,
f
1 " -
FIGURE 2 20 .
ASPEN PLUS SIMULATION OP REACTOR MODELS -f 69
Specifying components
The example reaction system includes three components. They are aniline, hydrogen
and cyclohexylamine. Defining all these species in the component input form, one obtains
Figure 2.21.
V nt Eik 4n eu To* FU. Pla Uh
f
3 Mdiilfs-3ij bj rl
~
AMIUNE C6H7I11
WyMOGEN K1T1R0G H
Ffesctons
Eire V/cw) UtwCMnd Rtttdei:
' ""'""
D in
I I Sotd. | U«>M«Mt t
Ml
RSac Brttj tEqai RGfcb) lCStft RFtifl Rflaieh
r
FIGURE 2.21
i
I 3
urvac
_
j F rm
I*
FIGURE 2.22
70 PROCESS SIMULATION AND CONTROL USING ASPEN
J
-
3 i«*f
_
j PiAiW
Strunu
fj EOVar-ittai
FIGURE 2.23(a)
'
10
3 Owerti
i Jy MIXED
~
3
3 :iu*f..
-
3
-
O UMFACQtsun 3
« Zj EMMbn
ra; ToW IT
BCSIB «fl, m»
FIGURE 2.23(b)
d « « |.if.|.iu«-.| ne i
9 - .
i -
ff .
j .
.
J
- 1,--.
I
1
F 3
~
J
Si-
__-__
iir
ITXUK
r |®- 9 . S . 9 Q U O Mm fJte >-. -'« »-
FIGURE 2.24
Product streams have been defined with their phases (see Figure 2.25).
Ifflll
71
I.
li:
.
I
-
0 -p uj
I1XJM
m 0 . 8 . o y JE D IMMI Ptmt hm. unit TVl ggg
FIGURE 2.25
Press Afexf button or click on Reactions and get the window (as shown in Figure 2.26).
Copyrighted material
72 PROCESS SIMULATION AND CONTROL USING ASPEN"
I UtatCSTRflCCTn i
f
«b Ed* ««» DKB Tooit tr Pta IJfewy Wirvis Hdp
f
e|ej rgklaKKM "I ! I"l J JJ J ®l
f
_
JJ
-
i 1 i HT I leal; I M Hi
, Setup
MvanMd
a h
.
11 L
Bocta
-
C5TR
e s-up
(J EOVsnai>«
O EOhpu -
O Sp«c GfWJpt
Pott
leschon E3 ID
Peaili
-
ii
[1 " MiMMiyS{«ter9 | S«p«aU»i ] He*E hangefi j Coluwij Heoclou j Preiwe Changeii ] MfloipUaloft | Sf** | UiwH&Wi |
-
a->
Mated \
STREAMS REeril RGtb! iCSTR RPbg
f
'
FIGURE 2.26
Right click on Available reaction sets, hit New button, then either accept default
name R-l or give a name as we want for the reaction set and inally click on OK f .
Subsequently, select POWERLAW in the Enter Type list and hit OK to get the screen
as shown in Figure 2.27.
MHl
r .-l-i- PT -.1 M- I Ml jW
.
/Spccft hm j/S««atm/Ba ljont] PSD j CwvwMAm [
£j Pwwt/Henccs S dect (sacbcn tw to t« ndudsd r nwdel
LJ
"
1 Mdecwer Su-me
J : -
_ j Data
i V) tco-S*:
f
Jfl H
'
i o l
csrn
O '
9 Sp«Gto*i
H SbeenRMub ,
<0 s s s o O
LJ.
FIGURE 2.27
ASPEN PLUS SIMULATION OF REACTOR MODELS 73
S .
FIGURE 2.28
form, we need to enter the stoichiometric coefficient as well as exponent for all
components. The exponents represent the order of the reaction with respect to each
component. Note that there are two types of reactions [kinetic (rate-controlled reactions)
and equilibrium] permitted under Power law reaction ID .
M. | |
j
-
hB- 0 . i . 0 1 1
iifitw4i f.4M Bf»J Km wifl
FIGURE 2.29
74 PROCESS SIMULATION AND CONTROL USING ASPENT
respect to aniline. Accordingly, we may use the following information to specify the
reaction (see Table 2.2).
TABLE 2.2
aniline -
1 1
3 1
hydrogen -
1 0
cyclohexylamine
Recall that in Aspen Plus terminology, coefficients must be negative for reactants
and positive for products. As we ill up the form, it looks like Figure 2.30.
f
illtiliiii ESS
'
"
ANILINE 1 j ; CYCLO-01 1
-
IYDR0GEN .
3
j *
* 1 i
Ctote
Bock,
y Reactiom
r J Chemolry
B Peacuons Edt Delete
t R-I
f
1 Convefgcnce
fj Rowaheetng Onions
1
.r
Reojrad tipul hcowMe
IT Mam pKen | Sepaators | HeatEndiangen | Cokams Haachm | PtenueOiaven ] SoUt UnModeb
KWariel
STREAMS RSIoic RYieU
1.0 .y-U-U-
HMj RStb. RCStfl BPIm Rieldi
"
ReuMtnO
« « b3
FIGURE 2.30
If we do not specify the exponent for a species, Aspen Plus takes a default value of
zero. In Figure 2.31, the resulting relation is displayed in the stoichiometry sheet.
In the subsequent step (see Figure 2.32), we move on to Kinetic tab.
ASPEN PLUS SIMULATION OF REACTOR MODELS 75
PHPI Liasigl -
.
if?:.
3aft l"-" JSldilJP BiiJfllalfil
-
.
j am
I» j 0-.
- 1 <
) - .IW I I I Ihmt
<o- y i Q on o
"
KiMWiingwr r ~ » !
FIGURE 2.31
Irl |x|
IM .Q C3WA\*\<M H
* "I I"! -I vl -I 9|
r .
"
' jfl
- I 3
9
9 .
t into mn*&*B**n*
f
t .
ta
3
.
jfl .
.
-
*
P '
* jfl
FIGURE 2.32
76 PROCESS SIMULATION AND CONTROL USING SPEN]
basis. Accordingly, the
'
'Molarity
As directed in the problem statement, we use
Power law is expressed as:
n
E n 1
r= k exp (2.1)
[T0; R
where r is the rate of reaction, K the reaction rate constant (kinetic factor in Aspen
Plus terminology), k the pre-exponential or frequency factor, T the temperature m degree
K Tn the datum temperature in degree K, n the temperature exponent S the activation
energy R the universal gas constant, C the molarity in kmol/m a the concentration ,
where,
E
K = kTn exp (2.3)
RT
In most of our simple cases, the reaction rate constant is represented by the Arrhenius
law, that is
E N
K - k exp (2.4)
RT)
Notice that when the Arrhenius formula is used we put zero for n and nothing for T0
,
in the Aspen Plus window. Also, the units of the pre-exponential factor are identical to
those of the rate constant and vary depending on the order of the reaction As we .
know, the dimensions of the rate constant for an nth order reaction are:
(time)-1 (concentration)1-'1
Next come back to the problem The kinetic data are required to provide in the above
.
sheet. Here we use the Arrhenius law to represent the reaction rate constant. It is
important to mention that the pre-exponential factor must be specified in SI unit. For
the example CSTR problem ,
the pre-exponential factor and activation energy are given
as 5 x 105 m3/kmol s and 20 000 Btu/lbmol respectively (see Figure 2.33).
,
mam
.. ._ - -J »«. Pte Un* OMn m,
f LMMMto twite
s csm
3
ll ; US
a-i
*
Kdlarai."
stream
farHsfe.pnMn
t
FIGURE 2.33
:
Ffc Dm Taofai Run Lfesry Wirdiw KHp
JJ 1 1- i,JV
,
-Hal
-
iLWjilSla)
5
@ CSTPJ
fo K o .. .
.-
FIGURE 2.34
Viewing results
In the next select Solver Settings, choose igsuto Summary/Sf ms in the list on the
f
left and finally get the results shown in Figure 2.35 in a tabulated form.
78 PROCESS SIMULATION AND CONTROL USING ASPEN1
I f I I i PT ! .leal I - Ml tM
" I i I i i I I i
"
3 '-" l
il il -
am nil 1
1000 0541
0«5
0J30 0 001
MUM nooo
-
sm DOM MPPM
0 98)
R&tac RVWd
i -1 .QMi-O'
SEgJ HQtei RCSIR RPfaa RftWi
FIGURE 2.35
The combination of two benzene molecules forms one molecule of diphenyl and one of
hydrogen (Fogler, 2005). The elementary reversible vapour-phase reaction occurs in a
plug low reactor (PER).
f
1250oF and 15 Psi. The data for the Arrhenius law are given below .
Simulation approach
Select Aspen Plus User Interface. When the Aspen Plus window pops up, choose
Template and click on OK (see Figure 2.36).
i -
...
.
...
iwmmmlt mm
FIGURE 2.36
In the next step (see Figure 2.37), select General with English Units and hit OK button.
1 V-
I
-
- -
FIGURE 2.37
Copyrighted material
TM
80 PROCESS SIMULATION AND CONTROL USING ASPEN
Creating flowsheet
In the Model Library, select the Reactors tab. Expanding the RPlug icon, the following
screen is obtained (see Figure 2.38).
li,-1?:- ?-- IM
-
Uj _
SIftEAMS
jS's - s - § o
'
RStoc lY»fc) W»J RCte RCSTR RBtfd<
f
FIGURE 2.38
Inserting the left bottom symbol in the Process Flowsheet Window adding the feed and ,
product streams, and renaming the block as well as streams, finally we see Figure 2.39.
Be £* *> &M ro* tj> Uonn WnSo* H*
f
>|[T><rr| h~o
,
I* -****** | f«M». t hmI- mw | c*-« iu««« I rM..1,o, _
i 1 '
S'W
iS- SSI Gj
-
FIGURE 2 39 .
ASPEN PLUS SIMULATION OF REACTOR MODELS 81
Configuring settings
At this moment, we are sure that the process low diagram is drawn correctly . The Status
f
message directs us to provide the input information. Hitting Next knob and clicking on
OK, we obtain a form for setup specifications. First we input the Title of the present
nroject (Simulation of the RPlug Model), followed by the accounting information
(AKJANA/7/ANYTHING/AS YOU WANT) and Report Options [see Figures 2.40(a) to (c)] .
f
ISrolWoneilheHPVjgMocW
Vdd|*MMC |
-
o->
SIKAMS '
HSteic tVM
i s u -= u
myt RG|te. RC?tR Rptq RBtuh
t
FIGURE 2.40(A)
UaTSil
> nt Mm OKk TMIp An W L±>»v WWo* H*p
arsi
aiobdj-/Deicnmn >/Acciwnlina| 0>agr>o«(«ci {
SIRLWi
hB- 1 -1 - 8 Q
l5ia R.'* RfrMl RGfaU; W Iff Hf''-.
O
f
"
U
gMdiM-AiMf But " I
__ '
.
FIGURE 2.40(b)
82 PROCESS SIMULATION AND CONTROL USING ASPEN
r m«i
» r SM
»O-S-0 y
FIGURE 2.40(c)
Specifying components
From the Data Browser, select Specifications under the Components folder. As we provide
the chemical formula of the components in the Component ID column, the other columns
of the table are automatically filled up (see Figure 2.41).
< Fit E* Htw D«i Tat* tji RsT Ihwf fntotr Hife
f
IMM FvmU,
cia<io
_
j «r<rCwwi
SET
"
r* rg '. Bin
k«W iooxi Id* -J
d
rvJ.cnttm
'
FIGURE 2.41
tmum
3 1 I
is I
I
3
d
"~
r.
I 3
. » r
FIGURE 2.42
_
i r _ .
IF
{y MIXED 3
State vsmUm
" " "
3 3r 3
|12SJ |f 3
1-5 |p. d H2
UH**C i3rtu»
Toid lwr (m.,- 3
f
|0 02 jbmot/N 3
Miinii
"
Tdat IT
Hoi 'jmvUf t.
i-1 . § . § u-i j
'
FIGURE 2.43
84 PROCESS SIMULATION AND CONTROL USING ASPEN
i Wl tl P*l Lb try Wn
f
f
Jill IMMM
'
r-.i:i i nr
*
-
J Ntw QuwclWMldn
'
3
J/j Prcf n«
-
J
> -a f
aock-
_
J Readx
_
J C -wssxe
I Sold* { UnrMixM. j
BE** ftGMn
FIGURE 2.44
Open the Configuration sheet and enter the reactor dimensions in the next form
(see Figure 2.45).
_
j . /- itarj
36
DwtmUt 06 . -J
S _
j PAttwn
?»cp(rtif Veered
a :
a §is q-u
'
RStot Rrail lE<M Wite' BCS'R Bp
f
HUM
FIGURE 2.45
ASPEN PLUS SIMULATION OF REACTOR MODELS 85
T the subsequent step, we define a reaction set for the simulation. The default name
R-l has been accepted. Then select Power law kinetics and obtain the picture, shown in
Figure 2.46.
' r
M *" D*« ** "rw*
I 22
'
i)
SIR£>M$
H®-1 j . I y
' HStac fiW RCqU RGtti RCSTR iB
a BB*J>
f
FIGURE 2.46
and (b)) for reaction number 1 (2C6H6 -> C12H10 + H2). Since the reaction rate constants
are defined with respect to benzene we convert the stoichiometric coefficient of benzene
,
to unity for both the reactions Obviously the reactions are second-order.
.
,
Rmmm
I1-I'
v R 1
ii
-
o -
-
0 S 0 1J
-
a .a # 1
FIGURE 2.47(a)
86 PROCESS SIMULATION AND CONTROL USING ASPEN
HI »l
l ! |h| -i ~ij j
05
J _2=J
j i
. :
f ue HB**
'
FIGURE 2.47(b)
Aspen Plus uses a default value of zero. As the message on the screen, shown in
Figure 2.47(b) reveals, it is true that the forward reaction rate does not depend on the
product components. After completing the irst reaction, select 'New' rom the Reaction
f
f
No. list. Enter '2; for the reverse reaction QHe C H + H2) and click OK (see Figure 2 48). .
l -
.
. ff alaltfrfi
-
.1 ilBl: I si
f
«|-
r a
1
Oeate a nm Redcton No .
PR
R-t
""
'-
-
--- II *p»f\«-a»i «
FIGURE 2 48 .
433
1 n-i 1.1 nr -
i ,ieii
. i mi *m
i; i
71
J _i5Lj
FIGURE 2.49
Hit A exf knob and obtain two stoichiometric relations as shown in Figure 2 50. .
-
.
3>>J qLJniJ
HmNo Stuctimttry
: Kn«c
_
j MHnnd
, I Sehdt I Ui*M«Jrt )
u E .11
,.
61 bio's
*
c
FIGURE 2.50
In the simulation of the present problem we use partial pressure basis (applicable ,
for vapour only) and therefore, the Power law expression has the following form:
,
( f >
n
E ri 1
r = k exp (2.5)
R To ,
where , P represents the partial pressure (N/m2). If fo is not specified, the above equation
18 replaced by:
88 PROCESS SIMULATION AND CONTROL USING ASPEN
For the prescribed reactions, values of the pre-exponential factor and activation energy
are provided in the two forms, shown in Figures 51(a) and (b). To apply the Arrhenius
2 .
law, we put zero for temperature exponent n and left the box, allotted for datum
temperature T0, empty.
i»f.i
I r mi r» ! .isi; I - IB!
as
±1
"
1 d
E §
ill
* ai F
* a ?
Si Bacfa
i PR
0 R-1
FIGURE 2.51(a)
». Ea »«, 0«, r i, a .
.
i '
r u>i-«i» rr 'i-.joii
-
"
a i
KiMtel«daNUT/T>>|"*'(E<n|m'l/TB|
i a *******
I Mill
FIGURE 2 51(b) .
ASPKN PLUS SIMULATION OF REACTOR MODELS 89
_
i r-i I ! f»
_
'
-i-igi 1 w aisd
(0 9 S 8 O = UM t<< Of*
FIGURE 2.52
(a) Viewing results: Click on Solver Settings knob, choose Results Summary/
Streams in the column at the left side and inally obtain the results for all
f
" J
a S d
IUWi 1 -m
- -
tm-1
" "
TW no-
4 1-
- -
*m sub-
-
aaraocc- ve- im
(stcsss: -
rwm 1
L
ROT Ml
Mr
-
oi
TWI
M
I - I M»l
inuA mm mm mik --
-I * I-
FIGURE 2.53
C ll
90 PROCESS SIMULATION AND CONTROL USING ASPEN
(b) Producing a plot of mole fraction vs length: Use the Data Browser menu
tree to navigate to the Blocks IPFRI Profiles sheet (see Figure 2.54).
MJi HillLlim-WPMlMli -. ..l»1.T71.
»in: «.
Fte &*l
Wi V*t
V*» C#»
C#» Toe*
roe* H i
Hji PW
pw L*f«7
Ltmv > "«
.
ifl - P.OC..I Sbe«. I
it Pt >wt«
t aa 9»um
f
-
e v.
-
a pf
_
pn F z
Utt Sutra.
far-
bt !S iHo
f] nwdb IS TZZf&VS
"
; 15 r?55 sSiTFW
.
4
s !5 00001 u*?
-
m
a 9fM- p« Z' 5
A m RNdm lb i55S
-
QniMgra i* 15
IS last 1 [ri!DK4IIt
f
g LSiQFM,
it. j
FIGURE 2.54
In the next, select Plot Wizard from the Plot pulldown menu Alternatively, press .
a S5
I
: PlOCCUtilMnKtXEflM
1 a t
f
9 EM IJ
if/ t** Ocw
jlE
E
t Fa* ;
f
ii
24 IE
ir
-
9 a EOCor-Ortcm ;
fj LSSOPBu -
L J J
51 REAMS
1
1 HSbe
i . i y=o
HVMJ REcU ROttx
FIGURE 2 55 .
Click on Next icon and get a variety of plots (see Figure 2 56). .
ASI'KN PLUS SIMULATION OF REACTOR MODELS 91
MM
pi-eniajaaisi =i r. -i ht a3f j
t
---
I
a tfmm
i
9 -
? - m N
1
-
n 1
17 f
R
-
n
ri
n 1
(
< H
1
J I i
iTmao -
Nftj mfc < w
- " -
ind* .. -_;=_
f
FIGURE 2.56
Among the available options, select one plot type that is titled as 'Composition' and
press Next button (see Figure 2.57).
r-l |..l'fT
-
- 'i-lci 1
fi ita.l
3 i±d «JP-3a -''ail -i -
.
a - i r
-
io -
-
j
-
m -
a Bin
3S
v- I i-*
-
3 '-
4
'
{© 9 . i 0 Q -O-
'
w > mm "tj~ mm m> mm J
FIGURE 2.57
Again click on Next and get the form, shown in Figure 2.58.
92 PROCESS SIMULATION AND CONTROL USING ASPEN1
V) - , PlOCBU SUUM j I I
i
.
iJ PHI f''V.I fi
,. -1
I
5r .
i
.
i {
CIS-
t
ri t
!
_
Cvitl 'Sack ir i
-
CH
STREWS RYaM
si y
TlSitd mt*» SCS?fl BWug
u BflWch
1« j- 9 -B- «M
FIGURE 2.58
Check whether the information displayed in the window, shown in Figure 2.58 are ,
ok or not. Hitting Finish knob. Figure 2.59 is obtained by plotting 'reactor molar
'
composition (mole fraction) as ordinate against 'reactor length' (in) as abscissa.
111 *
'
.
" ' 8M,
FIGURE 2.59
Note that the plot window can be edited by right clicking on that window and
selecting Properties In the properties window
.
,
the user can modify the title, axis scale,
font and colour of the plot Alternatively, double-click on the different elements of the
.
plot and modify them as we like to improve the presentation and clarity.
ASPEN PLUS SIMULATION OK KKACTOR MOOEI 93
2.
5 ASPEN PLUS SIMULATION OF A RPlug MODEL USING LHHW
KINETICS
Problem statement
In acetic anhydride manufacturing, the cracking of acetone produces ketene and methane
according to the following irreversible vapour-phase reaction:
CH3COCH3 -> CH2CO + CH4
acetone ketene methane
This reaction is irst-order with respect to acetone. Pure acetone feed with a low
f
f
rate of 130 kmol/hr enters a PFR at 7250C and 1.5 atm. The kinetic data for the Aspen
Plus simulation are given below.
k = 1.1 s"1
E = 28.5 x 107 J/kmol
n=0
T0 = 980 K
The unit of pre-exponential factor clearly indicates the |C 1 basis. To use the Langmuir-
,
Hinshelwood-Hougen-Watson (LHHW) kinetic model, set zero for all coeficients under Term 1
f
and that for all coeficients except A under Term 2. Take a very large negative value for
f
coeficient A. The sample adiabatic PFR is 3 m in length and 0.6 m in diameter. Applying
f
the SYSOP0 base method, compute the component mole fraction in the product stream.
Simulation approach
As we select Aspen Plus User Interface, first the Aspen Plus Startup window appears,
as shown in Figure 2.60. Choose Template option and press OK.
1 -I
**mmm*mH MM
FIGURE 2.60
94 PROCESS SIMULATION AND CONTROL USING ASPEN
In the next, select General with Metric Units and again hit OK button (see Figure 2.61) .
pea
f
<*-Sxar
'
Penmen
Mair>
'
"
11 'C*
'
FIGURE 2.61
Creating flowsheet
From the Model Library toolbar we have selected RPlug reactor and developed the
,
s,flt
H8- S . 8 - <SSL QU U
zzz rsr izf
M awif- ~ -
FIGURE 2.62
ASPEN PLUS"1 SIMULATION OK REACTOR MODELS 95
Configuring settings
In the list on the left, choose Setup /Specifications. For the present problem, we wish to
'
give the Title as Simulation of the PFR'. and accounting information as 'AKJANA/8/
ANYTHING/AS WE LIKE'. In addition, choose 'Mole' and 'Mass' fraction basis for the
streams under Report Options [see Figures 2.63(a), (b) and (c)l.
' i LU.
-
Ml o . § 6 onu
m»t «>>.
FIGURE 2.63(a)
FIGURE 2.63(b)
Gopyngt-
96 PROCESS SIMULATION AND CONTROL USING ASPEN1
ti-n» o hi NAKtad W
t
FkMbM hi--- 'i
P Hde PM* TFf, IGEN M.
3
T Mm. P
r
f
FIGURE 2.63(c)
Specifying components
Select Specifications under Components folder in the Data Browser window As we .
out the Component ID column Aspen Plus provides the rest of the information in ,
1 f -
i i-i- r» jiAm \ m
"
3 M£ i J nJ -3 »l Qj -.1 «*!
3 S-L*
O SfamOM*
Tim
tCEIO-JE
SottTSe )3<roi
itENE KEIENE bHJO
$ a«pm)<i seths -pn
»
|1
i Bk«k>
stficwi wiac f .
-
8 . 1 -y-lE-U
Pfcu
'
ns tu ncsTp
f
FIGURE 2.64
JtUI
l r Ll_L_F
_ _
-iCI I ! !
9 "w-cwr.
I 3
j» -
. la
mr'iiir
.
(0- 0 I t»*
: I jn U- ia ).
FIGURE 2.65
in Figure 2.66.
I.UH
-
;
I
, ,
3 i--- 3
I '
-
.
r
I--... --
I '- -
3
figs?
I'* 3
g M
f|7
. »
-
1 i iT
I- 0 . i . 8 OMU
»»«« gjfc «ani ggi
FIGURE 2.66
Copyrighted material
98 PROCESS SIMULATION AND CONTROL USING ASPEN
i ffc Ed» We* Oto Ta* Rn Rot lisrary VAnJo* He* 313
OUIHI I I itelal i?l rsK NKiH »'| 3 _
UjiJ _J -il SI
_
J AM H
.
_
J y c
-
9MM
-
Jfl F
O BUM*:
3;-* Ctro'
~
| | = 41111)
Mta
O EOrw -J
STREAMS RVaM
.
8 i US IJ
BEqut RGhta BCSTfl SPlup HBatet
f
. .
FIGURE 2.67
In the Configuration sheet reactor length as well as diameter are given (see Figure 2.68).
,
I 3
"
36 a SM
3i
t
StMrtl "
.
J6 f 3
rj f
& tdta
f
PnAn
f EOVw-b*.
a EOmpJ -
0
D
SWMlttf I HtHEttl
D-* iisi (
-
FIGURE 2.68
In the next we define a reaction set for the present simulation The default name
, .
R-l is ok Then select LHHW kinetics and obtain the screen exhibited in Figure 2.69.
.
,
ASPEN PUJSTM SIMULATION OF REACTOR MODELS 99
1 W-l-K-l..
-
a .
(@- 0 g g u a u
1 . - WW -G- Hi-' WS.
FIGURE 2.69
i r.:i..i-u rr . mi-*ifc|
d «- 1 d
*.Ml' -
t 2 roK-. i
.
I*
-
in«Mt
Q -
-
<@ S 8 § Q»0 *a>. 'mh gi»
FIGURE 2.70
100 PROCESS SIMULATION ANnjWQlOLUSING ASPEN'
Hitting on Next and clicking Kinetic button, we get Aeldn ics input form. A littl e
description is given below to understand the use of LHHW kmetxc model m Aspen
simulator.
The LHHW rate expression is represented by:
(kinetic factor) (driving force)
r = (2 1)
(adsorption expression) .
The kinetic factor (reaction rate constant) has the following form:
'
n
E \ 1 >
K = k exp (2.8)
R k
T To)
If Tq is ignored, Eq. (2.3) replaces the above expression. Note that all the notations
used in Eq. (2.8) have been defined earlier.
The driving force is expressed by:
f N A
n c?
concentration term C used in the above discussion is dependent on the [CJ basis Say .
for example when [CJ basis is selected as molarity the concentration term represents
,
,
the component molar concentration (kmol/m3); similarly when [CJ basis is partial
pressure, the concentration term represents the component partial pressure (N/m2).
Providing required data we have the illed kinetic sheet shown in Figure 2.71.
,
f
reactants,
set acetone exponent to 1. Similarly for products set ketene and methane ,
exponents to 0 Also enter zero for all four driving force constants as mentioned in the
.
In the subsequent step (see Figure 2.74), select Term 2' rom the pulldown Enter
f
term menu Since the given reaction is irst-order with respect to acetone
.
and there is
f
no second term enter zero for all exponents and coefficients Owing to
,
the method .
Aspen Plus uses to specify a reaction, we should insert a very large negative value for
coetticient A (say ,
-106) to make Term 2 essentially zero [see Eq (2 9)1 Finally click ,
on Next icon .
n- .
ASPEN PLUS SIMULATION OF REACTOR MODELS 101
-
L T -
I I 'i r-i -I -lei I !«!
3
-~ .
.
j ii w
, 3
j9
O 'm<iaii '
-
i
*
Zj i am i(t/T«f
lUlllll cuM.
:
f
J
.
Zj
a *-
.
j f
. .
- i
D -.
Ml
-
1IWJ'
<5 l-B
FIGURE 2.71
I Mk | - --
FIGURE 2.72
102 PROCESS SIMULATION AND CONTROL USING ASPEN
1
j
-
[-EETEEXSC
Idlbsw I
j
_
'.D- o-xr'
-
I
. p eis ep
f
.
1
BEoJ AG**!
FIGURE 2.73
Rwcaigthsai jVapm
3
Hi
ErteHtrm [l«rm2
.
a
_
j
t(*clartr
Expowii
0
) .
lJ Q
u
_
) Data
a
.
SbMM
F
. p
.
21 - i f/fT
PFS - Ml T SeeHfle .
t
lj Omttry
G
j Ci ytw tyw »fival«n enerw 'a t«
.
Jn poww law wpittWin
j
_
u l»
Mom
STRCAHS
BE(M RSfca RCSTR
f
FUin HUM
FIGURE 2.74
comer of the window shown in Figure 2 74. Subsequently run the simulation and obtain .
,
_
j r i
_
- -
i-i'
'
nr -i \ m -sw
i --»- tii
»**"' I
'
II I -*« -- I . *"l »->
£r <@ 6 S 0 O H U
IIIMH HlK «* Mm ".iril A* Mar
FIGURE 2.75
Viewing results
Pressing Solver Settings knob and selecting i?esw/ s Summary /Streams, we obtain the
final results as reported in Figure 2.76.
i i-liisialiil:
"
1 1 I -i
1 3
r[.- , .... el
_
fami ST
____ -
iffan ur
r
rzw
'
i
BTTiW I'M
Tivi
'
ii
FIGURE 2.76
Copynghied material
104 PROCESS SIMULATION AND CONTROL USING ASPF.N
PROBLEMS |
2 . 1 Ethyl acetate is produced in an esterification reaction between acetic acid and
ethyl alcohol.
acetic acid + ethyl alcohol <-> ethyl acetate + water
The feed mixture, consisting of 52.5 mole% acetic acid, 45 mole% ethyl alcohol
and 2.5 mole% water, enters the RCSTR model with a low rate of 400 kmol/hr at
f
750C and 1.1 atm. The reactor operates at 70oC and 1 atm. Both the reactions
are first-order with respect to each of the reactants (i.e., overall second-order). For
these liquid-phase reactions, the kinetic data for the Arrhenius law are given below:
Forward reaction: k = 2.0 x 108 m3/kmol - s
E= 6 0 x 107 J/kmol
.
2 . 2 Repeat the above problem replacing RCSTR model by RStoic model with 80%
conversion of ethyl alcohol.
2 3 Simulate the reactor (Problem 2.1) for the case of an RGibbs model.
.
and 5 mole% acetone, is introduced in a CSTR at 10 atm and 1250C and a low
f
TABLE 2.3
ethane 0 301
.
5 314
.
acetone 0 581
.
19.94
Simulate the RYield reactor and compare the results (mole fractions in the
product) with those obtained for Problem 2.4.
2 6 As stated in Problem 2.1, the reaction between acetic acid and ethanol gives
.
data for the LHHW model are given below (Fogler, 2005).
k = 0.072 s"1
£ = 82 x 103 cal/mol
Tq = 1000 K
|C,] basis = Molarity
The reactor length is 3 m and diameter is 0.8 m. Using the SYSOP0
thermodynamic model, simulate the reactor.
2 8 Repeat the above problem replacing the PFR by a stoichiometric reactor with
.
'
-
rk = K
26586
K= exp 22.8- mol/lit s . atm15
T
25589
K' = exp 19.62- mol/lit . s
Here, -rA is the rate of disappearance of acetone (A), -rk the rate of disappearance
'
of ketene ik), PA the partial pressure of A, and K and K the reaction rate
constants. Pure acetone feed with a flow rate of 130 kmol/hr enters the reactor
at 7250C and 1.5 atm. The reactor with a volume of 1 . 4 m3 operates at 700oC
and 1.5 atm. Applying the SYSOPO base method compute the component mole
,
REFERENCE |
Fogler ,
H. Scott (2005), Elements of Chemical Reaction Engineering Prentice-Hall of India
,
31
.
BUILT-IN DISTILLATION MODELS
An Aspen simulation package has nine built-in unit operation models for the separating
column. In the Aspen terminology, these packages are named as DSTWU, Distl, RadFrac.
Extract. MultiFrac, SCFrac, PetroFrac, RateFrac and BatchFrac. Under these categories,
several model configurations are available. Note that Extract model is used for liquid-
liquid extraction. Among the built-in column models, DSTWU, Distl and SCFrac
represent the shortcut distillation and the rest of the distillation models perform igorous
r
calculations.
DSTWU model uses Winn-Underwood-Gilliland method for a single-feed two-product
fractionating column having either a partial or total condenser. It estimates minimum
number of stages using Winn method and minimum reflux ratio using Underwood
method. Moreover, it determines the actual reflux ratio for the specified number of
stages or the actual number of stages for the specified reflux ratio, depending on which
is entered using Gilliland correlation. It also calculates the optimal feed tray and reboiler
as well as condenser duty. Remember that this model assumes constant molar overflow
and relative volatilities.
Distl model includes a single feed and two products, and assumes constant molar
overflow and relative volatilities. It uses Edmister approach to calculate product
composition. We need to specify a number of stages, e.g. feed location, reflux ratio,
pressure profile and distillate to feed iD/F) ratio. Actually, when all the data are
provided, we can use this column model to verify the product results.
RadFrac is a rigorous fractionating column model that can handle any number of
feeds as well as side draws. It has a wide variety of appUcations, such as absorption,
stripping, ordinary distillation, extractive and azeotropic distillation, reactive distillation, etc.
MultiFrac is usually employed for any number of fractionating columns and any
number of connections between the columns or within the columns. It has the ability to
simulate the distillation columns integrated with flash towers, feed furnaces, side
107
Copyrighted material
108 PROCESS SIMUKATION AND CONTROL USING ASPEN
strippers, pumparrounds, etc. This rigorous column model can be used as an alternative
of PetroFrac, especially when the configuration is beyond the capabilities of PetroFrac .
It is used to model refinery columns, such as atmospheric distillation unit (ADU) and
vacuum distillation unit (VDU).
PetroFrac is commonly employed to fractionate a petroleum feed. This rigorous model
simulates the refinery columns, such as ADU, VDU, luidized-bed catalytic cracking (FCC)
f
fractionator, etc., equipped with a feed furnace, side strippers, pumparounds and so on.
RateFrac is a rate-based nonequilibrium column model employed to simulate all
types of vapour-liquid separation operations, such as absorption, desorption and
distillation. It simulates single and interlinked columns with tray type as well as packed
type arrangement.
BatchFrac is a rigorous model used for simulating the batch distillation columns. It
also includes the reactions occurred in any stage of the separator. BatchFrac model
does not consider column hydraulics, and there is negligible vapour holdup and constant
liquid holdup.
It is worthy to mention that for detailed information regarding any built-in Aspen
Plus model, select that model icon in the Model Library toolbar and press Fl.
In this chapter, we will simulate different distillation models, including a petroleum
refining column, using the Aspen Plus software. Moreover, an absorption column will
be analyzed. In addition to the steady state simulation the process optimization will
,
32 1
. . Simulation of a DSTWU Model
Problem statement
A feed stream, consisting of 60 mole% ethane and 40 mole% ethylene enters a DSTWU
,
column having a low rate of 200 Ibmol/hr at 750F and 15 psia This feed is required to
.
f
component in the distillate and 99 9% of the heavy key component in the bottoms. The
.
sample process operates at 300 psia with zero tray-to-tray pressure drop The pressure
.
in the reboiler as well as condenser is also 300 psia In the simulation, use total
.
30 theoretical stages (including condenser and reboiler) and a total condenser Applying .
the RK-Soave property method simulate the column and calculate the minimum reflux
,
ratio, actual reflux ratio minimum number of stages actual number of stages, and
,
,
feed location.
Simulation approach
From the desktop select Start button and then click on Programs, AspenTech, Aspen
,
,
Engineering Suite Aspen Plus Version and Aspen Plus User Interface. Then choose
,
Template option in the Aspen Plus Startup dialog and hit OK (see Figure 3 1). .
ASPEN PLUS SIMULATION OF DISTILLATION MODKUS 109
I ' l-l-l I- 1 1 I
-
I ! 1-1
FIGURE 3.1
Select General with English Units as the next window appears (see Figure 3.2).
.
4./
MM
-
Hi . -
...
mm ...
__
l
FIGURE 3.2
C aterial
110 PROCESS SIMULATION AND CONTROL USING ASPEN
Again press OK to see the Connect to Engine dialog (see Figure 3.3). Here we choose
T,ocal PC by scrolling down. Hit OK knob and move on to develop the process low diagram
f
.
Connect to Engine a
Server type:
User Info
Node name:
User name:
Password:
Working directory:
( OK 1 Exit Help
FIGURE 3.3
Creating flowsheet
As we select Columns tab in the bottom Model Library toolbar (Figure 3 4), Aspen Plus .
FIGURE 3.4
ASPEN PLUS SIMULATION OK DISTILLATION MODELS 111
In the next, select DSTWU icon to represent the short-cut distillation process.
Once we have selected the icon, place the icon on the lowsheet by clicking with the
f
cross-hair somewhere on the lowsheet background. When inished, click on K | symbol
f
or ight-click on the lowsheet background. By default, the column is named as Bl
r
f
(see Figure 3.5).
i\n 'amiami I
Hi tM Dn <oaa 'hr nrann lw -«
UJ -
iT
_
-
CH
"SAW. ' DIIMI Out "*l<m 1M MtfMi IW l*. ..
i c-.i C- a'aMAcwi a IM AM ru- MMC*
FIGURE 3.5
In the screen, shown in Figure 3.5. only the block is displayed; there are no incoming
and outgoing streams connected with the block. Therefore, the Status message in the
bottom right of the window includes Flowsheet Not Complete. Interestingly, after
connecting all required streams with the unit, this message sometime may also be
retained. This happens because of improper flowsheet connectivity.
To add a single feed stream and two product outlets (distillate and bottom), click on
Material STREAMS tab in the lower left-hand corner. As we move the cursor
(a crosshair) onto the process lowsheet, suddenly three red arrows and one blue arrow
f
appear around the block. These arrows indicate places to attach streams to the block.
As we know, red arrows are required ports and blue arrows are optional ports. Click
once on the connection point between the feed stream and the DSTWU block, enlarge
the feed line and finally click again. By default, this stream is labelled as 1. In the
similar fashion, we can add the two product streams, namely 2 and 3, to the distillation
unit (see Figure 3.6).
Copyrighted malarial
112 PROCESS SIMULATION AND CONTROL USING ASPEN
0 -
&
Ul . J
ilf Mewt/SpUeu 1 Stpaatai | HwlEMhangwt Criumn* j ReKloit | PrwawOiangsi | MsripuWwi | SoW« | UMtMwWt j
f
STREAMS ' DSTWU Dntl Rrfisc EntisO Mutftw; SCFiac PeOoFi Ratrf.ac BWchFiac
Book flOflCfcrJ-
r
FIGURE 3.6
fte EiJI «ew OKa Tocti f**i fte«h«i Uxsy Wndm. H-t
0 O
E -
STREAMS
'
DSVM Rrf(»e E**d MtAffc SCfttc PMcfiae RjuF.ic Eatctfr
.. , - ~r IR Wfc. . nR*!!! MUM TW*o
'
Vn co- I Qw i W l- O il W l UN T * || A -a .
f
FIGURE 3.7
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 113
Now the Status bar in the window, shown in Figure 3.7, says Required Input
Incomplete indicating that the lowsheet is complete and the input specifications are
f
required to provide using available input forms for running the Aspen simulator.
Configuring settings
Recall that within the Aspen simulation software, the simplest way to find the next
step is to use one of the following equivalent commands:
0 -
c-Q-
3 I
| *,0-H-» Wot I tJO t ttoM I <]OmH»1 Mnw |aj<i»Otto»»« || AvxfV S~ Ql'f. t*9t
FIGURE 3.8
Hitting OK on the above message, we obtain the setup input form. Alternatively,
select Solver Settings knob and choose Setup /Specifications in the list on the left
(see Figure 3.9).
Although optional, it is a good practice to ill out the above form with a title and to
f
t
L r i-i
_ _
rv 'I -ipi i w
I "'M'.
'
SintMC WIMAJ Cxi Brf.K MiJtfue 5Ct'«c Pk. i tttrfii Brf.»f«:
f
FIGURE 3.9
_
J W«w
i (Mi
" -
I 3
= 3 EOCsn.C«n
'*
,
7 '*"' I ""I 1 -TO- I Mrow<w I Sold: t UtvUaMt (
FIGURE 3 10 .
In the next (see Figure 3 11) the Aspen Plus accounting information (required
.
.
l\ I r",
9 -
'
3f
a) -.
a mi
5!
) MB I - - - I
FIGURE 3.11
We may wish to have stream results summarized with mole fractions and/or some
other basis that is not set by default. For this, we can use Report Options under Setup
folder. In the subsequent step, open Stream sheet and then choose 'Mole' fraction basis.
In this regard, a sample copy is shown in Figure 3.12. although this is not essential for
the present problem.
i-d 3-J-itiiJa.ii
0 -»
1 .«*>*>
FIGURE 3.12
116 PROCESS SIMULATION AND CONTRQLUSING ASPEN
Specifying components
Use the Data Browser menu tree to navigate to the Components/Specifications/
Selection sheet (see Figure 3.13).
. :
f
. .
1 J
.
5 Drill-id!
Uyr-End Pmpoti
rHudocariiMW
ti aeehi
iO II SiiKi ait to tr rMneved Itom dsiatw* J. erte< C<m(rt<*fK Haw « f om-ia SteHdp
O ,
[if Mtw SpUen | S«p»*« | HealEttJiWBen Cohmnt } lwcto« 1 P-essueO owt | MwpuWWi | Sate | UraM«W. j
f
U *
- PatftFi ? Rahjiac BalchFrae
r
Sr«*,B«.ft
-
r:siHa-Ai>i,ivini1
FIGURE 3.13
In the window, shown in Figure 3.13, the table has four columns; they are under
the headings oi Component ID, Type, Component name and Formula. Among them, the
Type is a specification of how an Aspen software calculates the thermodynamic
properties. For fluid processing of organic chemicals, it is generally suitable to use
Conventional optiom Remember that component ID column should be illed out by the
f
chemical formulas or names of the components as their IDs On the basis of these .
component IDs, Aspen Plus may spontaneously ill up the Type Component name and
f T haPPen' * that AsPen Plu« to ind an eXaCt
,
f
T86 fj!?
f
in
IDs
(see lfir COmPonent hane and ethylene, as thefr
3.14). ThefootherT
6 0f
(see Figure three columns have been automatically filled out. i
.
Set RK-Soave
property method by scrolling down (see Figure 3.15).
ASPEN PLUS SIMULATION OK DISTILLATION MODELS 117
..i-j ;
3 .
FIGURE 3.14
3S3
Plata l I wi Qb3Mslllid5d 3
I r-l I..|-f7 .: .id ! '« -It
zsjti D»ld -3l alig|g|
» BS aJ
9 -
i ! r
-
FIGURE 3.15
Copyrighted material
118 PROCESS SIMULATION ANnCO TROLJ-JSING ASPEN'
i r - I .i-l rv . j J-Igl M.
mm m
.
J PltKTMM 21
-
a) enjr,h«»caa r f THWIE
|pmtu«
-
-
:
RK£6'J1 ;
l
nrtxu i
ri UNKK Owe-
rj um(KO j>(
m i
1
111' if
FIGURE 3.16
Filling out the form, shown in Figure 3.16, with the data given in the problem
statement, one obtains the data, shown in Figure 3.17.
.
He Ed! Vc« ata Took e> FM Ltrary M*km -i*
t
3Mi
.
J/) sT»*r. M«hoC ] EOOpbora |
Conmotnon
-
g , - .-
| Mole f-*: 3r
h Corrconan
2
-
|Pini«o IS
RXSBU-1
[is 04
RKTKUI
Q E*anM*To-is
1J «.
.
_
j «.w
lew [T
ll _
-
o ~*
SIB6W6 : bStWU
i
J- «
.
! »2l
FIGURE 3 17 .
FIGURE 3.18
Under Column specifications option, here we enter the number of stages that is 30.
It is fairly true that we can alternatively specify the reflux ratio when the number of
stages is asked to compute. Note that ethylene is the light key and naturally ethane is
the heavy key. As mentioned in the problem statement, recovery of the light key
component in the distillate (= moles of light key in the distillate/moles of light key in
the feed) is 0.996 and recovery of the heavy key component in the distillate (= moles of
heavy key in the distillate/moles of heavy key in the feed) is 0.001. In addition, the
pressure of the total condenser and reboiler is given as 300 psia. Entering all these
information, one obtains the result, shown in Figure 3.19.
i f -
l-l t' -
1 I
f
,,
--
|
-o -
FIGURE 3.19
CopynghlGd material
120 PROCESS SIMULATION AND CONTROL USING ASPKN'
iT TTJ
f
f
KiHictoipooert
.
s*.
Com [ETMYLEHE
Smm
t
flBMy: [0 0C1
E0V»Mblw
_
33 Dsmni
a -
Q Be-* C-pua-i
Mb
STRWMS
FIGURE 3.20
Click OK on the above message and obtain the Control Panel window that shows
the progress of the simulation (see Figure 3 21). .
-
EH
HUM
it -. I
f
FIGURE 3.21
ASPEN PLUS SIMU1.ATI0N OF DISTILLATION MgggUj 121
Hitting Next followed by OK, we have the Run Status screen (see Figure 3.22).
_
i r -i
_
rr .i.ipi i ibi
-
.
f
HMtf
mutt WIWW (M fMF(« iMMO *fik SOik i*rfi«
f
HH-'Ifci '*. .
FIGURE 3.22
Viewing results
In the next, select Blocks/DSTWU/Results rom the Data Browser. In the following
f
Save the work by choosing File I Save As /... in the menu list on the top. We can name the
ile whatever we like. Remember that a backup ile (*.bkp) takes much less space than
f
Copyrtghtod material
122 PROCESS SIMULATION AND CONTROL USING ASPEN
dim
-
-
i 'j-i
.
i a SH
r
-
1S?497652
»3*J12.;3 Br.il.
-
ifM IM|NMM f
0399
HE IP
sti| s
'
_
I Oa j-WWtW Awcn PIl» - Static Q G9«
FIGURE 3.23
'
lalxt
i Edi Font \kB '
bescfiiPTiON -
Central simulation with Eoallih units :
F. psl, Ib/hr Ibool/hr, Btu/hr. coft/hr.
,
PSOP-SDUSCES PUHEll
C0KPOMEKTS
ETHANE C2H5 /
ElKfLEKE C2H>
PROPERTIES Pk-SOAVE
PROP-OATA RirSKD-l
IH-W.ITS ENC
PROP-LIST BKSKI3
BPVAL ETHANE ETHYLENE . OlOOMOfriJO
BPVAL ETMYLENE ETHANE .0100000000
iTREW F
S085TRE»t fIXEO TEKP"' 5. PRES-1S t«0LE-FLOW-200 .
.
-
»± | <1A
.
FIGURE 3.24
summary, stream information, etc., select Export rom the File pulldown menu. Then
f
Plus Version/ DSTWU.rep). In the next, open the saved report file (DSTWU.rep) goingd
through My Computer and inally using a program, such as the Microsoft Office Wor
f
or WordPad or Notepad. For the present problem, the inal report is shown below.
f
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 123
(i) FAR 52.227-14. Alt. Ill, (ii) FAR 52.227-19. (iii) DEARS
252.227-7013(cMl)(ii). or (iv) THE ACCOMPANYING LICENSE AGREEMENT,
AS APPLICABLE. FOR PURPOSES OF THE FAR, THIS SOFTWARE SHALL BE DEEMED
TO BE "UNPUBLISHED" AND LICENSED WITH DISCLOSURE PROHIBITIONS.
CONTRACTOR/SUBCONTRACTOR; ASPEN TECHNOLOGY. INC. TEN CANAL PARK.
CAMBRIDGE. MA 02141.
TABLE OF CONTENTS
FLOWSHEET SECTION 2
FLOWSHEET CONNECTIVITY BY STREAMS 2
FLOWSHEET CONNECTIVITY BY BLOCKS.. 2
COMPUTATIONAL SEQUENCE 2
OVERALL FLOWSHEET BALANCE 2
STREAM SECTION 5
EOF 5
LOCATED IN:
DESCRIPTION
BTU/HR, CUFT/HR. PROPERTY METHOD: NONE FLOW BASIS FOR INPUT: MOLE
STREAM REPORT COMPOSITION: MOLE FLOW
COMPUTATIONAL SEQUENCE
TOTAL BALANCE
MOLE( LBMOIVHR) 200.000 200.000 0 000000E+00
.
COMPONENTS
INLET STREAM: F
CONDENSER OUTLET: D
REBOILER OUTLET: B
PROPERTY OPTION SET: RK-SOAVE STANDARD RKS EQUATION OF STATE
BDF
STREAM ID B D F
FROM: DSTWU DSTWU
TO : DSTWU
SUBSTREAM: MIXED
PHASE: LIQUID LIQUID VAPOUR
COMPONENTS: LBMOL/HR
ETHANE 119.8800 0 1200
. 120.0000
ETHYLENE 0 3200
. 79.6800 80.0000
COMPONENTS: MOLE FRAC
ETHANE 0 9973
. 1 5038-03
. 0 6000
.
ETHYLENE 2 6622-03
. 0 9985
. 0 4000
.
TOTAL FLOW:
LBMOL/HR 120.2000 79.8000 200.0000
LB/HR 3613.7256 2238.9320 5852.6576
CUFT/HR 140.3489 82.0590 7 5963+04
.
STATE VARIABLES:
TEMP (F) 20.4654 -
18.3114 75.0000
PRES (PSI) 300.0000 300.0000 15.0000
VFRAC 00 00
.
. 1 0000
.
LFRAC 1 0000
. 1 0000
. 00
.
SFRAC 00 00
.
. 00
.
ENTHALPY:
BTU/LBMOL -
4 1532+04
.
1 6983+04
.
-
1 2638+04
.
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 127
BTU/LB -
1381.4403 605.3231 -
431.8608
BTU/HR - 4 9921+06
.
1 3553+06
.
-
2 5275+06
.
ENTROPY:
BTU/T.RMOL-R -
58.6713 -
30.5758 -
28.8269
BTU/LB-R -
1 9515
.
-
1 0898
.
-
0 9851
.
DENSITY:
LBMOiyCUFT 0 8564
.
0 9725
.
2 6329-03
.
BLOCK STATUS
* *
*
Calculations were completed normally *
*
All Unit Operation blocks were completed normally *
« *
*
All streams were lashed normally *
f
« #
322
. . Simulation of a RadFrac Model
Problem statement
We will continue the above problem with few modifications. A hydrocarbon stream,
consisting of 60 mole% ethane and 40 mole% ethylene, enters a RadFrac column having
a low rate of 200 Ibmol/hr at 750F and 15 psia. The distillation process that has total
f
68 theoretical stages (including condenser and reboiler) and a total condenser operates
at 300 psia with zero pressure drop throughout. The distillate rate, reflux ratio and
feed tray location are given as 79.8 Ibmol/hr, 8.75 (mole basis) and 41 (above-stage),
respectively. Consider the RK-Soave property method.
(a) Simulate the column and compute the compositions of top as well as bottom
products.
(b) Is there any discrepancy in product compositions obtained rom RadFrac and
f
Simulation approach
(a) Start with the General with English Units Template, as shown in Figures 3.25(a)
and 3.25(b).
128 PROCESS SIMULATION AND CONTROL USING ASPEN
I I
FIGURE 3.25(a)
Click OiiTin the screen, shown in Figure 3.25(b). When the Connect to Engine dialog
pops up, again press OK button to obtain a blank Process Flowsheet Window.
FIGURE 3.25(b)
Creating flowsheet
Among the built-in columns in the Model Library of Aspen Simulator select RadFrac
,
and place it in the lowsheet window Connecting feed, distillate and bottom product
f
r|gf7-| .|..|' pr .
1 -tCi i iw IH
t 13 0
- J r- -
<
FIGURE 3.26
Configuring settings
In the subsequent step, simply hit Afet button followed by OK to open a setup input
r
form. These two windows, shown in Figures 3.27(a) and (b), include the Global and
Accounting information for the present project.
. % fa M» Dm tM» * imh n m n*
9 i"-
- . I - 11 - I
9
i ii
-
3 >- w.
0 -
FIGURE 3.27(a)
Copyrighted material
130 PROCESS SIMULATION AND CONTROL USING ASPEN
Skv
iwi-s«i
f> 5V-C6AH
Sttwrnf
v a
*
Flo»*e«bftg Opbon*
t1 a MoM ,1-! Took
f
-
4 Gfl
i 2Jil
FIGURE 3.27(b)
; H«MtebancUWinsD««nMpwl
Ffaubwit FiMjonbMa ShWDfamM -
Wo«tWu.iartwiw*
MOM
SIRCMC <i-t-o:-§-iv-(i..#.c-r-
' DSTMl
DSTMJ ag Arf E*ao. SCf >1>rf l .ac e IK
f
FIGURE 3.28
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 131
Specifying components
In the list on the left, choose Components /Specifications to define the components.
Using the component names, ethane and ethylene, as their IDs, we obtain the illed
f
table as shown in Figure 3.29.
r-i-i-i
-
i .w
-
i -igi *m
3a !«
Mn*
;lrM*N£ cats
. IHYUNC
>
O si
o ,.
-
FIGURE 3.29
Soave base method to compute the physical properties (see Figure 3 30). .
' Fit CM
0|<*|H| - ;
>I>'.|--|T »MBJ|v
.
mrai J «J|a.
3
B«i«»a(tw6 fflr.SOAVt
.
- - r 1 J
OwoCTylO j
IJ .
It
FIGURE 3.30
132 PROCESS SIMULATION AND CONTROL USING ASPEN
3 |m* dr
it ,
(200 jbrt*.
1j
FIGURE 3.31
a «*
H*o4erHcuv*i
7]F 1 ?98 (trrotA,
O Dwt i
Ha
-
FIGURE 3.32
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 133
Under Setup subfolder, the filled Streams sheet looks like Figure 3.33.
i 'to
. 1 « «!«.;
> | | MB I M»» |
it r
-
FIGURE 3.33
In the next, simply input 300 psi under Stage 1/Condenser pressure. Aspen simulator
assumes that the column operates isobarically if no additional pressure information is
provided (see Figure 3.34).
IB I' W tl*)
i:.,ir.ir.ii.0.ii'.fi..#.s .j-.
FIGURE 3.34
f
I l"l "li-dail 4;|
l i-.i T ai Jfilp
~
_.
iJ -
-
>C«lcaJ.«ttQn» t»«in . . .
Coif/«cg*nfi» iearttlOK*:
f
01 KL IL Sce/Tol
ill Eo.oas
a i j U.m
J 13 7.oc«a
4 1 « i.73*fr
5
.
.
i : i o.i6a»7i-ai
lli>ra»l
STREAMS PSTWU Dirf Radfw , lAaci MUtfxic; SCFuw fahoFrac Ratrffac Batetfrac __
_
FIGURE 3.35
Viewing results
Click on Solver Settings followed by Results Summary and Streams, we have the table,
shown in Figure 3.36, accompanying the results of all individual streams. Save the
work in a folder as a ile.
f
3 Sii.i»I«e|
J CcnvOotcm
1 i i
" 1
OS* Um 0517
"
- cx.magttgicffTr-wi
FIGURE 3.36
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 135
f
of the present problem (see Figure 3.37).
ire«.i uiu-y itHin bv up*- »lia .>. It.I « ll.'S'l* Sal ).i 11, 00'
Olr.ciory I Pro m (tin -Oin! .s't-t .ic.r Hut II. 1
tiHu« ei-» -
rmkM C!"
W'OITl Ultll'l
1 . ..
FIGURE 3.37
TABLE 3.1
ethane 0 996
.
0 004
.
ethylene 0 004
.
0 99G
.
TABLE 3.2
Component B D
ethane 0 997
.
0 002
.
ethylene 0 003
.
0 998
.
From Tables 3.1 and 3.2, it is obvious that there is a little difference between the
product compositions. However, the main reason behind this fact is that the RadFrac
performs rigorous calculations, whereas the DSTWU is a shortcut model. Another
possibility is the round-off error associated in the reflux ratio and feed tray position.
Copyrighled malarial
'
136 PROCESS SIMULATION AND CONTROL USING ASPEN
3 .
3 ASPEN PLUS SIMULATION OF THE MULTICOMPONENT
DISTILLATION COLUMNS
33
.
1
.
Simulation of a RadFrac Model
Problem statement
C3 5
/-C 15
4
n-C4 20
' C5-
25
A?-C5 35
Simulation approach
(a) As we start Aspen Plus from the Start menu or by double-clicking the Aspen
Plus icon on our desktop the Aspen Plus Startup dialog appears (see Figure 3 39).
,
.
FIGURE 3.39
VSI'KN I'll 'S SIMl'LATION OF DISTOIATIOM MHHKl.S 137
As Aspen Plus presents the window after clicking OK in Figure 3.39, choose General
with English Units. Then hit OK (see Figure 3.40).
FIGURE 3.40
Click OK when the Aspen Plus engine window is displayed (see Figure 3.40).
Remember that this step is specific to the installation.
Creating flowsheet
At present, we have a blank Process Flowsheet Window. So, we start to develop the
process low diagram by adding a RadFrac column from the Model Library toolbar and
f
drawing the inlet and product streams by the help of Material STREAMS.
Now the process lowsheet is complete. The Status bar in the bottom ight of the
f
screen, shown in Figure 3.41, displays a message of Required Input Incomplete indicating
that input data are required to enter to continue the simulation.
Configuring settings
Hitting Next knob and then clicking OK, we get the setup input form. In Figures 3.42(a)
'
and (b), the Title of the problem ( Simulation of a Multicomponent Column') followed
by the Aspen Plus accounting information (AKJANA/ll/ANY ID/FINE) are provided.
Include the additional items in Report Options/Stream sheet under Setup folder
(see Figure 3.43).
C aterial
PROCESS SIMULATION AND CONTROL USING ASPEN
-
CM
STREAMS ' DSm) Diiii Kvfiec E«liacl Muffiac SCFw PelroFiK Hattfims BUchFrac
'
FIGURE 3.41
|a<xl
-
O RepalCptom
d
Conmneti
OuQtitanAt |ENG
Prao«t« CONVEN Si
d
rj Rearforj d
Uu...... '
< 1 J
d
r UwfreswAtcataJabm
_
J RaKliSurray
. SetHat,
FIGURE 3.42(a)
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 139
f
_
I f - l -I-1 fV -1 -iC! i m\ MM
o--
1-1
FIGURE 3.42(b)
.
i -iff] i 1*1 -iti
- I
5
r .
'
FIGURE 3.43
Specifying components
In the left pane of the Data Browser window, select Components /Specifications. Filling
out the Component ID column, we obtain the table as shown in Figure 3.44.
Copynghied material
140 PROCESS SIMULATION AND CONTROL USING ASPEN
it SmJalicn 0*:<-i
f
PROPANE Convert onal fflOPWIE
SOBUfll S06UTAHE C4H10-2
CuBjn Uh*l
0 Papc"
wet* MEmtlBUIiJIt. 9UZ-3
JPCItTWIE
Prtcaties
FIGURE 3.44
Specifications in the column at the left side (Figure 3 45). From the Property method .
pulldown menu, select PENG-ROB. This equation of state model is chosen for
thermodynamic property predictions .
das
Ffe E* W Pa W
0|e?|B|
I W
si «lStl
Jt Specficabom
tw-waioi method j -
f
htufaoampemrti
Qj H-f»)rC<mpj
; . unif«:qw*>! : ChcntbylD j
O GLOBAL
1 (3 Enmaw
1 -ij PramMwi
F
-
CH
STREAMS Pflefnc P fwt B#etf.K
t
FIGURE 3.45
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 141
. ....
Ml II- P=!J 3
I2S I-
~
3
f» k 3 ,
& -
3
I 3
FIGURE 3.46
r£Z3_
3=
'
: 3
. mam
ea
t
FIGURE 3.47
142 4 PROCESS SIMULATION AND CONTROL USING ASPEN
In the subsequent step, specify the feed tray location in the Streams sheet as shown
in Figure 3.48.
f
I
CJ
.
31 P«p 3«.
Q
flow
i rj dv
20 Um4
S .
t . .
9 , '
0- RADFRAC
Hi j j lndxra
f
* Cj T3<*
FIGURE 3.48
-1*1 xl
5* pfboi
: f--3
Cj b**<**t
Ci 'jruJACCw
Cj uwff«&o
t cj «~*
Cj
E ( J Vt. j xr-
* CJ .
cj w
_ CJ earfiMngCkUn
f
'
!-
f*M* Mfr
FIGURE 3.49
the condenser and reboiler) we have the screen, shown in Figure 3.50 irst. ,
f
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 143
f
wal _U a|t1«>hNIH!£l 2l jiiiiii J ©i
r.>i..i..|. fT -
I .W I - W jgfel
'
o
Method
« Set*
Jj
0
:j
9
0
p llhMWHcu
Trr. Smg
»
0
0 Gsmaigra
®
0
T' 4
"
1
FIGURE 3.50 t
Press the knob to open the Vapour-Liquid sheet (see Figure 3 51). .
a '
a
Slartrg Endng
*
-
RSOfRAC
O a.
_
j P*> Sir.)
© Eta««
5-
/Cx- iele
flT <**tmm I I He»L£»-J«wn 1 *-". | Re«a». | l> i«ei>w | Meroiae,. ) Sol* 1 UwHoA*
-
a-»
MeMTiel '
. . -' 5 SEi . n. .... Dj.kT.
FIGURE 3.51
Assume the rectifying along with the stripping zone as Section 1 and fill up the
table , shown in Figure 3.52.
144 -f PROCESS SIMULATION AND CONTROL USING ASPEN1
mmmamm
RU> &H \Aew Dats Tec*
i r-i y i nr am i m
.
-
J <
'
joc*w. St«V9
Ugp
1 19
PAOFRAC
* 1
Mm CM
C vientm hi
Pock Sa j
EMMm
Report
User Suooj-
LiJ_
f -r.
STREAM?
NUM o
Mobetuot* It
FIGURE 3.52
±Jjij
< Loading SlmXat cn Jngtn* lS:61r33 i
D RADFRAC
IHTOUUTTflH
eitJAfly PUJutfOTiRS paru can* set d roa hoo«l SEParrD
mie muivn from scf tasli . xuu kms - lsptste
« esM
1 2195
.
c ;j*4t
Jsasll
"
V J
STREAMS
'
OSIVAJ
OSTWU Dfl« Rrfwc MJf« iCfiK FWrfnc R ft*: gaafwe
J'j jT«-tfc, i- *
FIGURE 3.53
ASPEN PLUS SINfULATION OF DISTILLATION MODELS 145
Viewing results
Click on Solver Settings knob and then choose Results Summary /Streams to obtain the
product compositions (see Figure 3.54).
rl 3'-f 13 "--l
9
_
J
IM Mi
mmf ~
" "
Hi 1 T» TBH
"
«.>* TOB 1 l» . "
J
»*».
-
sin 1
"
rtai-'
-
1|
1 !
Will ue in
TariT ' «« IM
lb i«»
i.«
"
aoni IJH- TIB TM
ntn 1 . a l<» III
Tf«5B 1
|
-
FIGURE 3.54
It is a good habit to save the work done at least at this moment. If we wish to see
the tabulated results with the process low diagram in a single sheet, simply hit Stream
f
Table button just above the results table (see Figure 3.55).
«-C3-
FIGURE 3.55
146 PROCESS SIMULATION AND CONTROL USING ASPEN
imut Su«Hry' cr«««d by Ajp«r Clul R«l. 11-1 it VMM Sun )ul IS, 2007
; Dlr«tory c froqriM fil« .AipefiT«ch\**ork1ng Folders Aspen Plus 11.1
. Fileniw C userj ijjn* *ppOjt« local T*<Bp -*pt>a cn,
"
tiuC
"
'
OtStRIPTION
C»n»ral Slaulation wtth English units ;
f, psl. Ib-tv. Itaol/hr. Btu hr. cuft/hr.
property method: none
fIoh basis for Input: Nola
Strcan report cowposlclon: t»ol« fl«
'
AOUEOUS / 50LIOS / INOBGANIC , A
PROP-SOURCES POREll .
aflUEOUS ,' SOLIDS / INORGANIC
CWPOWMTS
PROPAXE ClHS ,'
ISO»U-0l C4HI0-Z ,'
N-«VT-01 CAMlO-l /
2-t«ET-01 CiMl2-2 /
N-PEN-&1 C5H12-1
rSHCET
BLOCK RAOFR.AC IN-F OVT-OV B
PROPERTIES PEMG-ROB
. ol #
'.
f> " o* t QwpW M I JO IW l lhAJT M- [ wnfta 5 } > AtWwAort* || -APF7EAJ > 1 17a>
t
t
FIGURE 3.56
(b) First, choose Blocks /RADFRAC /Profiles in the column at the left side
Accordingly, we have the stage-wise data as shown in Figure 3 57. .
mm J >>j am w
3fl *«- rmi I c 1
Vapofto.
F oawnd
r *j
Tiit
1 asoU 1 "
B
--
o
t
t
3 aj i
'
9TIJ3K?SI 20 i 171 385017
T- t c I 195 3*a»6 20 j
f
20 i
liJnuiiJlii i
;oo 3ra 20 i
201354309 ia i
9 anwsssi }»
lit
n 121) ;: -
.
.
- i
llJ ||
:iT:;nr-si
In the next, select Plot Wizard rom the Plot dropdown menu or press Ctrl + Alt + W
f
on the keyboard to get Figure 3.58.
I ... '.
1H|_U*1«*) nKW*haMH n I |h| .131.31 _
j
71 .1
. «-
Cj owl-
| am
r
a 'mm
ijte>
!tDM5
I Until-
1
15"
-
CH
FIGURE 3.58
Click on Next button in the Plot Wizard Step 1 dialog and get a variety of plot types ,
i"l pi a)|
ift
~
.
Q?1!* OowfMto fV**! Mva
'
H4W
,
i
( smvm*
J 51 OK
_
J T.w Sot,
naa
j
_
wFodo FkMntln CCCOTHI CQCQS H) Hv4m > BM
J
- s- j seisa
«2J!4
Wti*1
1 Join*
1
60 - .
JP""
Pi"
-
CH
'
m& tS» so.. w« ac.-A»»H.MI NU«
f
f
iQVjre 1740
FIGURE 3.59
148 PROCESS SIMULATION AND CONTROL USING ASPEN
Select the plot type under the heading of Temp and press Finish button to obtain a
'
plot of Temperature (0F) vs. 'Stage' (see Figure 3.60).
rinwii-ii
I I J s < i io u £ 15 i« k w tr tt is-a
: it
FIGURE 3.60
Recall that the above plot window can be edited by ight clicking on that window
r
and selecting Properties. Then the user can easily modify the title, axis scale, font and
colour of the plot.
332
. .
Simulation of a PetroFrac Model
Problem statement
An artificial petroleum refining column (PRC), shown in Figure 3.61, consists of a feed
furnace and a distillation tower. The tower has two pumparound circuits, a partial
condenser and three side strippers. The furnace (single stage lash type) operates at
f
25 psia and provides a fractional overflash of 40% (StdVol basis) in the tower. The
outlet stream of the furnace goes to the tower on Stage 22. The tower has 26 stages
with a Murphree stage eficiency equal to 90%. A steam stream, STEAM, is introduced
f
at the bottom of the ractionator (26th stage with on-stage convention). There are another
f
three steam streams, STM1, STM2 and STM3, used in the side strippers. The condenser
runs at 15.7 psia with a pressure drop of 5 psi. The tower pressure drop is equal to
4 psi. The distillate rate is 10000 bbl/day and the distillate vapour fraction in the
condenser is 0.2 (StdVol basis).
Copyrighted malarial
ASPEN PLUS SIMULATION OF EHSTILLATIQM MODELS 149
LIGHTS <
WATER
STMl
sir,-
STM2
BOT
A hydrocarbon mixture with the following component-wise low rates enters the
f
furnace at 1170F and 44.7 psia (see Table 3.3).
TABLE 3.3
H2O 250
In Table 3.4, two pumparound circuits and three side strippers are specified.
TABLE 3.4
Loeatum Specifications
Pumparound Draw stage Return stage Flow rate Heat duty
(drawoff type) (bbl/day) (MMBtu/hr)
1 (partial) 8 6 49000 -
40 (for cooling)
2 (partial) 1 12 1000 -
17 (for cooling)
Location
Four steam streams used in the column model are described in Table 3.5.
TABLE 3.5
Specifications 1
Steam stream Location Temperature (0F) Pressure (psia) Flow rate (lb/hr)
STEAM Main tower 350 50 11500
STM1 SID1 stripper 350 50 4000
STM2 SID2 stripper 350 50 1500
STM3 SID3 stripper 350 50 1000
Considering the 'BK10' base method under 'REFINERT process type, simulate the
PetroFrac column and report the low rates (bbl/day) of all product streams.
f
Simulation approach
Select Aspen Plus User Interface. When the Aspen Plus Startup dialog appears ,
I I I I I I It
A+ f Blank Smiiabon
A+ Temptale
"
OpendnEwiingSimulatior
0«
I
FIGURE 3.62
As the next window pops up (see Figure 3.63), select Petroleum with English Units
and press OK knob .
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 151
mm
'Mm.
1
linn v-
.
..
L-
EJ J-J
.
. ..
FIGURE 3.63
Click OK when the Connect to Engine dialog appears. The next screen presents a
blank process lowsheet.
f
Creating flowsheet
Select the Columns tab from the Model Library toolbar. As we expand the PetroFrac
block icon, a variety of models is displayed as shown in Figure 3.64. Select a model
icon and press Fl to know more about that.
J? if" IT
[1 it Jb Jt
& ir if
ir .
# lb
t ir Op
f
tr ffc* c
J 1 -
4
. | turn* t
i> B p a .
# I .
-
W 11 1
FIGURE 3.64
152 PROCKSS SIMULATION AND CONTROL USING ASPEN
FIGURE 3.65
Configuring settings
Click Next to continue the simulation (see Figure 3.66). In the Title field, enter
'
t.
r
k
<UbM
-
.
<i.i-?-ir-Q-ii'-i?--i
.
»»No M - -o- Ut- W»
i.-r
W.
FIGURE 3.66
Copyrlghtf
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 153
In the form, shown in Figure 3.67, the Aspen Plus accounting information are included.
I'M MI | i l"!
1 1 F
-
9 -i*--- I 1 ,1'5-l-Mi-' 3 »\ Q| .i -1
-
9 M Mjri '-c
s
»1 9i<>>
. I"- 1
FIGURE 3.67
Specifying components
In the subsequent step, use the Data Browser menu tree to navigate to the Components /
Specifications sheet. Filling out the component input form, we have Figure 3.68.
1 7
-
3
T
,
- -
it! T KJETTH tSruT
-
r? f
Lj i-r
FIGURE 3.68
154 PROCESS SIMULATION AND CONTROL USING ASPEN
f
K-10 method).
Ea Tc<* tr HO th*> .v v t
f
NBI I I a>|e| ¥?! nlt'lftKKI I >»l ~l i I "I -I *\ I J
I f i i I fV .|:-liS|-. I .:M
1 >hl -3idPy!!il
-
"
ficwtw |Bfni;tFf.
I 3
3
< i-j :
1 aj c v-lu
'
FIGURE 3.69
D|tfiB| _
iJ e|jgJ njKjfti i-g -i ».i { |h| Miigj .} _
j
Pn zi bJiiif
-
3 »J QN h
n
|
n
b
iei
cs nm
!3J no
i 0.1
w,
1 ED-**
rv
-
D -
FIGURE 3.70
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 155
As we hit Next icon, an input form for Stream STEAM opens up. After illing out, it
f
looks like Figure 3.71 .
M-i rv -
1*0 |r
. F/traBy
£ Mol(Qior Stioam
i: FtranKn
ict
d ia 1
d H
11500
: g ST
E S
S LGHTS
3 SE>1
3 O SlOS
* a
"
Total NlSOO
'
itiputCMeM*
fV MMd iAw. } SwoWort 1 HutEttharpei' Columw j Rwawt j h w O* | M*v niE | S<*t> | U«Modes |
'
STREAMS DSIWU
ntnjl r. .1 Raifitc EntlKl MJtf . SCfiac t)_-.r...
Pt fwFrac njj
Bat
-
efiae
._
Bwctjj..
b
fiac
;
NUVI cj/Ki i Tote
O
- S . * AwnWM t'
FIGURE 3.71
In Figures 3,72(a) (b) and (c), three illed input forms are shown for STM1, STM2
,
f
JMI mi
dSfb GPET.r-j Ol l fZi T] Q|.M N>J
"
9
| EOOptm )
j'Aijean n*n*
: I
J to
_
_
y t- tr-, Minodi
C1
.
_
j d C2
ISO [p.. d -O
_
J 0*3 C4
d - ICI
- 1
cs
NC5
1 J H20
fEEL
tniii
kM Cans,..
i fM-I-8-If'(S'#'C'4
"
-
'
SWUM MTMI
'
MM* (renFI
FIGURE 3.72(a)
156 PROCKSS SIMULATION AND CONTROL USING ASPEN
O tn Ua J9J»J
S3 2J 1 SI
_
LiJ pr 1 1 1 jaj «iM
») LJ| | N.|
-
1
EOfJi eini
ApecMicotioo<| FtHhOtfHIni ) 1
l/j Cm<p<mIi BkfcibMmnaiq) fTMIXuT 3 J
31 [mmTfIcm 31"1"'
\m |i zJ 1 '.(.-'/r 1 V**
u,,
.
MdMculitSniaM r "
zl -
DA
id
'
Mam 1
ItfalltoM |SUwl id
| |l.H/dv zj cs
»j Steaon
. BIT
. ns 1. J ; m
_
,
.
JH FEED
. LIGHTS
f _
j sin
.
_
, sm
. srew
.
S s '
-
_
stm:
O kvol
RnuU
i .
Duc.
Rwftw. c ,_i
Exlud u. c.u
MiAFiac crc.v
SCFisc c .rX,
Pe'.cfi f4v«Fi»: 84ictfi«
Aspen Plus -S
FIGURE 3.72(b)
at:
Rs £c* We* Dae Tc s Fiwi FU : lAray VNVid&B
"
'
"
3M 5 [ g rg \ \ u j | n>|
~
SLtsdeamr me jyMIXED 3 I
| Frcc et
~
3 K'**-
|3» 1- Vdue
:t
|p, ,. d
|SD If. zi
--
3fl C4
ISidvcJ d riC4
901
1 d CS
£ J BIS
1 J
-
HA-i
i
a FEED
S ij LN5HTS
s _
j sm
g i SID2
._
j 5103
_
> jfl ST = AM
J! STMi
s STM 3 lad fuccT
© lt*xJ
P,
ED A
Llj-_
-
a-»
Rrfrac EJiaci MJf c SCF.« (W.k RaFw
FIGURE 3.72(c)
"
j .
TiU-J.J 1'
_
!l:'Mii,
Sehjoteta'i
o
P«rN*V o, Ligud
* a f-" Uiidtihu«E
3\-ur
S EOT
UGHTS
K BDI
i jj so;
g £j 903
>: M STEAM
STVl
M STM2
WATER
1
J
Ki vjij hGOFWw- ~
4 Main colunn 81
"
Mod
' MTi.ni iv-.o . '.r. ri-e- n. f... K. t-- nJ. ii-...
SWE-MS I OSrwU DitB BrfrK E-trwt HMik SCfitc P fix. (WeF.se B<(j*F(»(,
fs-.net, wean
'
FIGURE 3.73
As we press Next icon, the Blocks /PRC/Setup /Streams sheet appears as shown in
Figure 3.74.
'
'
Fie Edt W C*s T>-A PW Uxsy Wrec* Hefe -
f
jtHGRETRC <<j|Ai 3 j
5 JCcriifjiter, ystiasnej Steam |vn«*tue|wCariMnnn)vFijin«ca{
feed me*"!
31** Conwibon
5 TEAM 26
FEED 22 Funaco
26
ttSHTS Vapd
Fieewatei
DIS
Jfl STBI
.
nw
i, as*!
FIGURE 3.74
n Figure 3.75, the pressure sheet shows the condenser pressure along with the top
as well as bottom stage pressure of the distillation tower.
As given in the problem statement enter 0.2 in the Distillate vapour fraction field ,
rl«l-l
9 -
F P 3
-
F P 3
-
17-
- 3 'no r
.
_j imi
< o «
*
a S'
.
zj m
-
3 iTm
Himdl ' in
.
-I « - - o-
'
FIGURE 3.75
I ] 3
T _
UXTI
w
w:
itfl
FIGURE 3.76
In the next step (see Figure 3.77), the feed furnace is specified by selecting the type
of furnace and giving the values of pressure and ractional overflash.
f
Copyrighted material
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 159
tUiiLH
MM
.
_
!3 J J
_
- zi
| SoUt | UiaMotJtt }
'
.
R£*MS
.
L/ST-mj _
D*L R«*mc E*«» Mrffue StTuc PeboFiK R ef>*e £!5n«C
FIGURE 3.77
In the left pane of the Data Browser window select Blocks I PRC I Efficiencies and ,
Sertgr. Enlnq
>
j4
fir1
FIGURE 3.78
The three windows shown in Figure 3.79(a), (b) and (c), specify the side strippers
,
I IB] I SI MM
_
] Mi 3 ]0pUlrtalF«di | Uqud Retail j P'i
J
_
id Seho
'
ll 3 r 3
a Stj m "
_
J 60T a
-
J DiS
-
FEED
- "
TT-
-
ICHTS
S'ipgng iearr, f 3
r
3
~
SIDl
_
l SlMm lo bo«om pioduct tslw |Wion»l| f
Jj SID! "
j SO] r 3
jfl STEW
STM1 Flow jptcrfcahwi
5TH2
fiofioftTpoduct 1
a STHi
'
- -
EC
_
j , .
PBC
O Scbc
.
ii 3-;
.
il 53
O Healer Co
4ir*«ol:fa!)o«,nctbng(ebo Reboto itprttdtK m Ihe bdlofn sage
-
,1
jl Mawi/SplPen | Sepaitfoi) | He* Ewhanowj Colurani ) RwidDn | Pi«a«Q»ano«: i Man xiatai | SoWi | Ui«MocW: |
-
CM
« "
.
.
I; 13:01
FIGURE 3.79(a)
_ _
h i - IBI iSlti
as'
* forced
| OpbanalFeedt j Rerun
, rj ot Setup
t
U 015
.
Jfl fEEt
NUTlb8(0(tl«(|ei [* 3 3
Siropeiixoduol 1 102 0v«<hMdretuK |n 3
.
j 5DI
-
.
H 5:02 Slnpomgirwdun
: a SiO) I* Stiipoiij *te*fi |STM2 jj
.
STEAM
.
5TM1
Sieamtobollampiodi 1 It .
.
Ji sTMi r Rebcisfduy
.
!r,«
S 23 WATER OpboMt tebo<e. heal ihwmt
ij Brxki t " (. produd 3 3
i|
Srrf/oJ 15000
'
I P i»»wid» -1
-
St r.
*
Si
*
A S-2
> .
S-3
O ttw* Stwca
f
r
0 -
rj CoodMer Hcurv
iwifiii iiHi W
KMancdu dw. JSt
3SA tinM) Sb»«r S2
HutEichanoeii Cokanu | Raadai | Pnuuactrogni I MawMon | SdUs | UmMcxMi
FIGURE 3.79(b)
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 161
as-3
-
a GOT
! U ns
FEED
3
n l "
3 tJv*f>»adi«t«n fie -J
f
-
u UCHTS
. Zj SOI
S«DJ
nB
SOS
5TBW
STM1
1 iV) snc
r Rttolodur)' .
r
*
- j STM3
.
j WATtP
_
n
- FX. 4
1) Jj .....
STREAMS DSTWU
'
fvHtki fn .
FIGURE 3.79(c)
Although the Status bar says Required Input Complete, we have to specify the two
pumparound circuits connected with the main fractionator. Select Blocks/PRC/
Pumparounds in the list on the left. Click on New as the object manager appears. We
may accept the default Pumparound ID T-l Then specify the first pumparound circuit
'
.
Ql -l H
_
J ra» «
-
1) Sdukw
_ j
J W.srtrt
OHRltfMM
_
-
J «*
lj LXirTS
1, 5103
U s-o: 31-ur
STEAM
STM1
jfl snc
ino
j MHTSH
taw
$ fK
Jj "MV
HBEAMS
-
hfin: ' OStWU
OSIWU Dail if.*: ErtM HJ*.* SCfi*c P«ofr*c H«rf'*c
f
r
FIGURE 3.80
"'
162 PROCESS SIMULATION AND CONTHOL USING ASI'KN
f
|»'|H| |
i r n I F .
i
. .. p rr-T n U -1 q -.| ul 1 n.|
-
J i
-
J -
| :
.
f" -
hm*lm
7
'
-
I t:
.
_
j w;
.
V) REE
_
_
j SiOl (..m 1... . d
.
.
1
,
"0-
L-
Sfll 1 -.-. d
f
J Beck*
Jfl ' T - OS
-
Jfl PI
h
_
i
_
J
FIGURE 3.81
...
J rj_iUV 11 Im; i ibi i i
.
Fr nxtcj input i(Mift«Hiau
1 1 O.tlMl
FIGURE 3.82
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 163
Viewing results
From the Data Browser, choose Results Summary /Streams and obtain the table, shown
in Figure 3.83, that includes the low rates (bbl/day) of all product streams. Save the
f
work done.
l.;,ft-i
-
o-
FIGURE 3.83
To obtain the input information (see Figure 3.84), select Input Summary from the
View pulldown menu.
ITU UMilMly of 1 1
.-
Cl CM
ut c*<ao ;
Ml . ««0-|
"Co
FIGURE 3.84
Copyrighted material
164 PROCESS SIMULATION AND CONTROL USING ASPEN
34
.
SIMULATION AND ANALYSIS OF AN ABSORPTION COLUMN
Problem statement
A hydrocarbon vapour enters an absorption column below the bottom stage and the
absorbent enters above the top stage. The column operates at 75 psia with no pressure
drop and it has four equilibrium stages. The absorber is specified in Figure 3.85.
GAS-PDT
Absorbent
Pure n-C10
Temperature = 90oF
Pressure = 75 psia
Flow rate = 1000 Ibmol/hr
ABSORBENT
oAo-rttu
Gas Feed
Temperature = 90oF
Pressure = 75 psia
n-C 150
5
(a) Simulate the absorber model (ABSBR2 under RadFrac) and compute the product
compositions.
(b) Perform the sensitivity analysis by examining the effect of absorbent low rate
f
(c) Compute the absorbent low rate to keep 15 mole% of C, in the gas product
f
(GAS-PDT). 3
Simulation approach
(a) Double-click Aspen Plus User Interface icon on the desktop. When Aspen Plus
window pops up, select General with English Units Template as shown m
Figures 3.86(a) and (b).
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 165
; .
b
,
_£l a AHM -< r-J 3 ! I I Jgl-J _!
1 r OMn»E«F«Sai>H>
0 J.<* *h«MW«U»0
V
..
FIGURE 3.86(a)
1F
FIGURE 3.86(b)
Click OK when the Connect to Engine dialog is displayed and proceed to develop
the process low diagram .
f
Creating flowsheet
Select the Columns tab from the bottom toolbar Among the available RadFrac models,.
select ABSBR2 and then place it on the lowsheet by clicking with the cross hair
f
somewhere
the inlet andon the lowsheet background Right-click to de-select the block. Connecting
.
f
outlet streams and changing the all default labels we have Figure d.»/. ,
166 PROCESS SIMULATION AND CONTROL USING ASPEN
F dl
- Vto* tJ«8 Tccto Btr low wl liw/ Wn<*»«
O 1 «8 'H
-
O 1 :-AS'Ee: f
STREAMS S OSTWU CwU Ratfiac E»»[»c> MtlaP.ac SCfiv: PfUcfrac Raiefiae Batctfrac
; Slwl] S J " Botk j 4] Aapen_ Mo«by Mct | - Cha r3 Hritirft j jUwtdTwro Hwcao )[ Aapen PVj« . Soul « j
FIGURE 3.87
Configuring settings
In the subsequent step, hit Next symbol and ill up the three setup input forms as shown
f
_
L r i-i-I- PT
_
-
5**
© Spectftcj«k>n«
Q SmitMlcnOiMni
ISmJabon ot «i Abtoiba
"
R«#iiVPO 3
OuWtBttlt |ENG r "
DM C
FIGURE 3.88(a)
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 167
* -
FIGURE 3.88(b)
38
Tig-
IS
sir
-
:
FIGURE 3.88(c)
Specifying components
From the Data Browser, choose Components/Specifications. In the input form, shown
in Figure 3.89, all components are defined.
Copyrighted malarial
168 PROCKSS SIMULATION AND CONTROL USING ASPEN11
& .
! -m I W W
'
I "
5!ES
ii
-
11 LI
-
D -
FIGURE 3.89
IF
. > --
-
D-
W W >»«.» M«» -il> «-»»
FIGURE 3.90
' r i i i- rv i ie:
-
m aia
i
"
3 nr
IS 3 .
,
IV
H3-
FIGURE 3.91(a)
3j
7laLJlllrr-3»)aU;
"
I E0O«Mr« |
: 1 EC?*- >]
I .
tea
r L
"
'
> fi '
"
£J T-r1
FIGURE 3.91(b)
«i -fe!
SfEf
P 3 ,
j\ r
-
Li ir
FIGURE 3.92
170 PROCESS SIMULATION AND CONTROL USING ASPEN
Select the Streams tab to specify stream location. Under Convention, there are two
feeding options: On-Stage and Above-Stage. In the present problem, the top stage is
the first stage and the bottom stage is the fourth one. Therefore the absorbent is fed
above Stage 1 and the gas feed is introduced above Stage 5 (see Figure 3.93).
-
He &k Vfe* CW* Tec* Fin fa u , Wi-do- Up
Mai -I..[ **m *?! oMfel M *'l :2l_J_liiJjdElj*J.£iJ
i r 1.14- nr -i i iai
PtoimiI >/!:cnjcnni| 1 I
sws. Com*
WOfiBEH
5
© HMttr, Fto-iatm
_
J MfMffd } Urn* 1W 6*jb Flow :. U*.
| .
FIGURE 3.93
In the next step (see Figure 3 94), select Pressure tab to specify the pressure profile
.
across the absorption column In this case, the column is operated isobarically at
.
75 psia. Under Top stage / Condenser pressure enter 75 psi. Aspen software assumes ,
7]
fop ' Ccmkraii pftlfJI*
* l (t*Sf EEC
f
E AflSOftSER
[V I - I . ..MT-UJ-W: 1 W J
FIGURE 3.94
ASPEN PLUS"' SIMULATION OF DISTILLATION MODELS 171
ci»1H| *U M*} _
aMai H'"I n . I M .| i?i m £
r
-
J T-l l-l' F*
-
'I .IBI I laai «M
mjihTiiiiTiii liiriSij
FIGURE 3.95
Viewing results
Choosing Results Summary /Streams in the left pane of the Data Browser window ,
we
get the results as shown in Figure 3.96.
loHal am
_
Lr_ia_ J rr _
.
i
la .
1 *LiJ «)W
~
3 aaad
I
j|r d
'
Jll .
d
a aria HBB
ei Mm
»> ...
NC4 tuoiyi liH IMHB
a-1 *l)4s
r
mM OlOD
|
tcio 1000 np oca .
1
-
a-
'
_V ~, 26... r-j .., iESfad
sifiEtw' e '' j Dr.« ri.j,3, Ml<, '* SSfigB f,v'r'a- SftS
f
f
f -. kM- j.
I a-w a i e ss i iaajgLg
t
FIGURE 3.96
172 PROCESS SIMULATION AND CONTROL USING ASPEN
(b) In the sensitivity analysis we will manipulate the absorbent low rate and ,
f
examine its effect on the exiting propane concentration In the column at the .
left side, double-click on Model Analysis Tools folder and then select Sensitivity .
As the Object manager is displayed, choose New. On the next window shown in ,
Figure 3.97, Aspen prompts us for an ID. Enter 'C3' as ID, and click OK .
- t* idl .tea
f
1 MM Wl nMfcl-si-JI-i el r I L'iiliJJSll -1
_
i r L m rv
_ _
d |B|
1
.
_ , w
-
:.-
.
_
]
I C mrtn
I C-.**
i j Sohh | Utti
FIGURE 3.97
In the next step (see Figure 3.98), select New under Define tab. Then we are
'
prompted to enter a variable name. Enter C3' and press OK. Subsequently, the following
information are required to provide:
II I I ! . fT . 1 .lEI ,
1 |B|
1 |
StaM l&i-iTOI
- > -
J .
'
FIGURE 3.98
ASPEN PLUS*" SIMULATION OF DISTILLATION MODELS 173
Type: Mole-Frac
Stream: GAS-PDT
Substream: MIXED
Component: C3
Hit Next and select the Vary tab (see Figure 3.99) The manipulated variable iis
specified with the following data:
Type: Mole-Flow
Stream: ABSORBEN
Substream: MIXED
Component: NC10
Overall range
Lower: 500
Upper: 1500
Increment: 50
> Ht aa Vfem Csa Too» fu\ Pa {tar, >Vhd» >Mi 1*1 "i
-
i i»l.li:M;ii A
_
L r _
* 1>>I ! 1 *»1
/Dorr yvacj |«T«ulaM| TaUgr OMMbn I Oaard |
.
_
j On, Ogmrm VsxMtruitei 3
Tor
Cjcpona 3 Low |M0
Upo» [Two
Ina fsT
Hercxljbrti
u«e Tea*
Lm1 I
-
-
Train,
Ure2 |
U C3 LntJ P
'
» - -
MM
.
_
j D»Ff
23 r tttji
MM
UKJAHS ' CrtTVU FVlrfnc lKrfnc BWrt/'K
f
FIGURE 3.99
In the subsequent step (see Figure 3 100), . select the Tabulate tab. This screen is
used by Aspen to set up tables. Insert T under Column No. Then right click on the
adjacent cell under Tabulated variable or expression. Select Variable List and drag and
drop the variable name (C3) into the cell. We may also directly type 03' in the cell. '
Then run the simulation and get the screen, shown in Figure 3.101.
174 PROCESS SlMULVriON AND CONTROL USING ASPEN
'
» (. <« . Mi Im* * uam mmm M*
-
10 > Cfton*
'> -
et " .
ft . .- u~
4 "
s
-
D -
'
DSTWU CMI Mrfe wjftK KTm hmfm *ml~ Mo o
t
< w -m Ik
FIGURE 3.100
1 .
I'ICI I 111 !?1»!J
AM
ir
- Qc>
> nn
I :mii
In Tm.
t trv
I Mil
t Mai
lil/TU
'» «»«
K ua
MtM
i uut
'
MIMU Om IW<k Umo >w i« IOai DMTw taXw
FIGURE 3.101
\SPEN PLUS SIMULATION OF DISTIUJVTION MODELS 175
5 -
FIGURE 3.102
In order to represent the results graphically, highlight a column in the table and
select X-Axis Variable (Ctrl + Alt + X) from the Plot pulldown menu. By the similar
way. select Y Axis Variable (Ctrl + Alt + Y) for the next column. Then select Display
Plot (Ctrl + Alt + P) from the Plot menu and obtain Figure 3.103.
J T
-
I-1 l- iv i-g I i»< siW
I --. I - l - I
J
FIGURE 3 103
176 PROCKSS SIMULATION AND CONTROL USING ASPEN
(c) In the left pane of the Data Browser window {see Figure 3.104), open
Flowsheeting Options folder and then select Design Spec. We need to provide
this design spec a name in the same manner that we did for the sensitivity
analysis. Press New, enter DSC3* and click on OK. '
38
2i£B_uaas -'-i ~ _u-j.ii2a
'
3 ~--
5=.
o
-
FIGURE 3.104
Select ihTettJ under Define tab. Then enter 'CS' as a variable name and press OK. In
the next step (see Figure 3.105), the following information are required to input:
Type: Mole-Frac
Stream: GAS-PDT
Substream: MIXED
Component: C3
.
.
K-MJ-B-M-iM-f-
FIGURE 3.105
Copyrlghiod material
ASPEN PLUS1M SIMULATION OF DISTILLATION MODELS 177
In the subsequent step (see Figure 3.106), select the Spec tab. Design specification
data are noted below:
Spec: C3
Target: 0.15
Tolerance: 0.001
;Ti=E4Mi DSTVU Dag R»Jf>: Erftaci M frac SCF»c PihoF-JC flyrf.K cj'ctfix
f
FIGURE 3.106
Finally (see Figure 3 107), select the Vary tab and enter the following information:
.
> 4. Lfc Ar- Cm« to-Ji Hpi Ptoi Ubr, AWtow »H>
_
j r iii nr
_
imi i ibi i i
' X
.
U
' -
i
_
/-,
6
3 "
-
- J r
0 r r
0
Q ! M6
'
>
r-s.
_
i >
-
?
:
1
| HmMH I r>MMf'CM«M | Mr iW I SOW: | UwMaW- |
»gj ' I a
FIGURE 3.107
178 PROCESS SIMULATION AND CONTROL USING ASPKN"
Type: Mole-Flow
Stream: ABSORBEN
Substream: MIXED
Component: NC10
Manipulated variable limits
Lower: 500
Upper: 1500
As we run the simulation, we get the screen, shown in Figure 3.108.
el w| Mii*MaJi£)!id r -
I I "I 311 <M
j lai
K. K. XL Ici Ir.
1 t 1
I 1 1 C.1K31
ML IL ttt/To)
i i » c.aui*
3TH£*«S ! OSTWU Dan Brfrac Um* UtffiK SOHk ttfwe RaUfwc Badftic
f
f
r
ail * J " . »t
'
FIGURE 3.108
As we choose Streams subfolder under Results Summary folder in the list on the left
(see Figure 3.109), we get the absorbent low rate of 1179.467 Ibmol/hr to keep 15 mole
f
of C3 in the gas product. This answer we can also obtain from the sensitivity plot.
It is well known that Aspen Plus is capable to optimize a function Here, we will continue .
the above absorption problem (Section 3 4) for optimization. In the present study, we .
wish to maximize C3 mole fraction in the gas product (GAS-PDT) with respect to
absorbent inlet temperature (lower limit = 50oF and upper limit = 300oF) .
Simulation approach
First solve part (a) of the previous absorption problem It means, ill up the input
f
forms for setup components, properties, streams and blocks. In the next, simulate tne
,
t -r
4
« -
.
'
)
M
"
TM H
lit
KB
r
r
,
_
1 1
-
- Si.- - JUBJ -
-
i>'«
IT iMT 1
.
sna-fSBRI- - tc
U wra-; wim
M .» i .«<-'.
FIGURE 3.109
In the column at the left side (see Figure 3.110), choose Model Analysis Tools/
Optimization. As the Object manager is displayed, hit New button and accept the default
ID 0-1' Press OK and then New. Entenng variable name 'CS'. again click OK Provide
the following information to maximize C3 mole fraction in the gas product.
I'M*
N t)
FIGURE 3.110
180 PROCESS SIMULATION ANT) CONTROL USING ASPEN1
Hit Next knob twice and get the screen, shown in Figure 3.111.
T-Birry m.v.-.T'-iB
_
LilJ Jl]35l £vj:
-
ii
2J
dp
a
-
iJ
At o i
ir Hitut/SpEttm j
f
FIGURE 3.111
Right-click in the empty cell with selecting Maximize option. Then select Variable
List and drag and drop the variable name (C3) into the cell (see Figure 3 112). We can .
CbiKtive tundnn -
HQ OfWlc*
s 2u eSci
[5) Rma
S ij 01
U .. ..
i'r
FIGURE 3.112
'
In the subsequent step select the Vary tab. Under Variable number, as we choose
,
New', automatically the number T appears FillinR out the form, we have the window .
as shown in Figure 3 .
113.
ASPEN PLUS'" SIMULATION OF DISTILLATION MODELS 181
!5)| I \£\ j l
21 Cw
_
j rw-m
fi nw aw i
f
D -
FIGURE 3.113
Pressing iVex symbol and running the simulation we get the answer (see Figure 3 114). , .
The maximum C3 mole raction of 0.259 is obtained at absorbent inlet temperature of 179 80F . .
f
J 1 i i fir i bi
-
-
. 1 ibi
j 'Jr-.rr.
3 T i
>,
K1.«pu 75 00
noco 10*
foosaoo HH
USC0 £»-
wv.- lUHjU
"
. t
f
Wfa>
'
'
c? mm 1?. *
Ti PIVB
170(00 i»S3
FIGURE 3.114
has been presented Several separating columns, including a petroleum refining column
.
and an absorber have been simulated using Aspen Plus. The process optimization has
,
also been discussed with an example. The present study covers both the binary as well
as multicomponent systems Interested readers may try to simulate the models given
.
in the exercise .
182 PROCESS SIMULATION AND CONTROL USING ASPKN
PROBLEMS|
3 . 1 A feed mixture, consisting of 60 mole% ethanol and 40 niole% water, is to be
separated by using a DSTWU model having a low rate of 100 kmol/hr at 40oC
f
and 1 atm so as to recover at least 85% of the light key component in the liquid
distillate and 80% of the heavy key component in the bottoms. The column
operates at 1 atm with no pressure drop throughout. In the simulation, consider
the reflux ratio of 1.5 and a total condenser. Applying the Wilson property method,
simulate the column and ind out the minimum number of stages, actual number
of stages, and feed position. f
3 . 2 A feed stream, consisting of 50 mole% ethane and 50 mole% ethylene. enters a
Distl column having a low rate of 200 Ibmol/hr at 750F and 15 psia. This separator
f
runs at 300 psia with no tray-to-tray pressure drop. The pressure in the reboiler
as well as condenser is also 300 psia. The feed enters the model at 6th stage and
the column has total 15 theoretical stages (including condenser and reboiler)
and a total condenser. If the reflux ratio is 7 and the distillate to feed ratio is
08
compute the mole fraction of ethane in both the product streams with applying
. .
is 2. A side product (vapour) is withdrawn from 14th stage. Applying the Soave-
Redlich-Kwong (SRK) property method, simulate the column model and report
the product compositions.
Feed
Temperature = 110nF
Pressure = 175 kPa
Dj O
Feed stage = 10 (above stage)
Component Flow rate Sj cC.
(Ibmol/hr)
benzene 250
toluene 80
diphenyl 10 B\ $
3 .
4 A reboiled stripper is to be employed to remove mainly propane and lighter
components from a feed stream, shown in Figure 3.116. It has total 6 stages
(including condenser and reboiler) and no condenser. The bottoms rate is
100 Ibmol/hr and the column pressure is 150 psia throughout. Using the Peng-
Robinson thermodynamic method, simulate the RadFrac model (STRIP2) and
ind out the product compositions.
f
Copyrighted material
ASPEN PLUS SIMULATION OF DISTILLATION MODELS 183
Feed
Temperature = 40oF D -
Pressure = 300 psia
Feed stage = 1 (above stage)
c2 75
c3 150
n-C 175
4
n-C5 60
n-C 35
s
Applying the UNIQUAC property method, simulate the extraction model (ICON1)
and note down the product compositions.
Feed
Temperature = 30oC
Pressure = 1 atm
Feed stage = 1
Solvent - SOLVENT-
RAFFINAT
'
Temperature - 30DC
Pressure = 1 atm
Feed stage = 5
Copyrlghled malarial
184 PROCESS SIMULATION AND CONTROL USING ASPEN
3 .
7 An artificial petroleum refining column (PRC) shown in Figure 3.118 consists of
a feed furnace and a fractionation tower. The tower includes one pumparound
circuit, a partial condenser and one side stripper. The furnace (single stage flash
type) operates at 20 psia and provides a fractional overflash of 50% (StdVol basis)
in the tower. The outlet stream of the furnace enters the tower on stage 18 .
The column has total 20 stages. A steam stream, STEAM, is fed at the bottom
of the fractionator (20th stage with on-stage convention). There is another steam
stream, STEM1, used in the side stripper. The condenser runs at 15 psia with a
pressure drop of 5 psi. The tower pressure drop is equal to 5 psi. The distillate
rate is 12000 bbl/day and the distillate vapour fraction in the condenser is 0.25
(StdVol basis). The liquid product, SID1, is withdrawn from 5th stage with a
flow rate of 2000 bbl/day.
A hydrocarbon mixture with the given component-wise low rates (Table 3.6)
f
enters the furnace at 120oF and 45 psia.
LIGHTS
WATER
IS o
-
SID1
STEM1 -O
FEED
SID2 C>
STEM
BOT O
TABLE 3.6
1800
n-C 7500
4
30000
1-0, 42000
nrCt 250
H 0 250
2
ASHEN PLUS SIMULATION OF DISTILLATION MODELS 185
The pomp around circuit (for cooling) and the side stnpper are specified with the
following information (see Table 3.7).
TABLE 3.7
Location Specifications
i
I (partial) 8 6 40000 20
Location
Two steam streams, used in the column model, are described in Table 3.8.
TABLE 3
Specifications
Steam stream Location Temperature (8F) Pressure (psia) Flow rate Ob/hr)
STEAM Main tower 350 50 12000
Selecting the PENG-ROB base method under RE FINE RV process type simulate ,
the model using a PetroFrac column and report the flow rates (bbl/day > of all
product streams.
This
Thispage
pageintentionally
intentionallyleft
leftblank
blank
Rs. 295.00
@ 2009 by PHI Learning Pnvate Limited, New Delhi. All rights reserved. No part of this book may
be reproduced In any form, by mimeograph or any other means, without permission in writing from
the publisher.
ISBN-978-81-203-3659-9
The export rights of this book are vested solely with the publisher.
Published by Asoke K. Ghosh, PHI Learning Private Limited, M-97, Connaught Circus,
New Delhi-110001 and Printed by Jay Print Pack Private Limited, New Delhi-110015.
Part II
Rs. 295.00
@ 2009 by PHI Learning Pnvate Limited, New Delhi. All rights reserved. No part of this book may
be reproduced In any form, by mimeograph or any other means, without permission in writing from
the publisher.
ISBN-978-81-203-3659-9
The export rights of this book are vested solely with the publisher.
Published by Asoke K. Ghosh, PHI Learning Private Limited, M-97, Connaught Circus,
New Delhi-110001 and Printed by Jay Print Pack Private Limited, New Delhi-110015.
Aspen Plus Simulation of
Chemical Plants
4 1 .
INTRODUCTION
In the last three chapters, we have studied in detail the simulation of individual
processes, such as flash drum, dryer, chemical reactor, distillation column including
petroleum refining process, absorber, stripper and liquid-liquid extraction unit, using
Aspen Plus software. Here, by a 'chemical plant' we mean a chemical process
integrated with several single process units. The chemical process industries usually
include flash chamber, mixer, splitter, heat exchanger, pump, compressor, reactor,
fractionator, ilter and so on. It is easy to simulate even a large chemical plant by the
f
Problem statement
A hydrocarbon stream H is supplied at 50C and 2.5 atm. The pump Pi discharges the
feed F at 10 atm. In Table 4.1 the component-wise low rates are tabulated for stream H.
f
The schematic representation of the complete process integrated with a pump and
ive DSTWU column models (Cl, C2, C3, C4 and C5) is shown in Figure 4.1.
f
189
Copyrk
190 PROCESS SIMULATION AND CONTROL USING ASPEN
TABLE 4.1
10
35
50
130
200
180
200
n-C.
5
pi C1 C2 C3 C4 cs
For Aspen Plus simulation of the distillation train, required information are given
in Table 4.2.
TABLE 4.2
All distillation models have total 20 theoretical stages (including condenser and
reboiler) and a total condenser. For the light key (LK) and heavy key (HK), we expect
99.9% and 0.1% recovery, respectively, in the distillate of all columns. Using the Peng-
Robinson property method, simulate the distillation train and report the compositions
of all distillation products.
Simulation approach
From the desktop, select Start button followed by Programs, AspenTech, Aspen
Engineering Suite, Aspen Plus Version and finally Aspen Plus User Interface. Then
choose Template option in the Aspen Plus Startup dialog (see Figure 4.2).
Copyrighled malarial
\SI'KN PWB SIMULATION OP CHKMK \l PLANTS 191
si
In ro..-
AM n i.tu.
1 1"
-
.' I 1 «
FIGURE 4.2
As wo hit OK button, the following window appears (sec Figure 13). Based on the
units used in the problem statement we select General with Metric Uliits,
,
in-i
.
1»
1- . l-.
.
I J I"-
_
I
'
FIGURE 4 3
192 PROCESS SIMULATION AND CONTROL USING ASPEN
'
Press OK and obtain the Connect to Engine dialog. Select Local PC as Server type
and click OK. Actually, this step is specific to our installation (see Figure 4.4).
Connect to Engine
Server type:
User Info
Node name:
User name;
Password;
Working directory:
OK ] ExB Help
FIGURE 4.4
Creating flowsheet
The next screen represents a Process Flowsheet Window. Add a pump by selecting the
Pressure Changers tab from the Model Library toolbar. Moreover, in the library, select
the Columns tab and then choose DSTWU model to include five such columns
consecutively on the flowsheet. Notice that to incorporate a block click on the ,
appropriate icon and then place the block on the process lowsheet by clicking with the
f
cross hairs somewhere on the lowsheet background. Right click to de-select the block.
f
Now we need to interconnect the blocks and add the inlet as well as outlet streams.
Select Material STREAMS on the left of the toolbar at the bottom In the next, as we .
move the cursor to the process lowsheet window several red and blue arrows appear
,
f
around the blocks. The red arrows indicate required streams and the blue arrows are
optional. In the previous chapters, we have learned how to connect the feed and product
streams with a single block.
Let us observe Figure 4.5 to know how to interconnect the two blocks by a stream.
Here, irst we wish to interconnect the pump PI with the column Cl using the feed
f
stream F. Right-click with highlighting feed block select Reconnect Destination and
,
then move the cursor to click on an arrow that is fed to the column Cl.
m o
-
QD-o
SOD-*
r ]-o
PI Cl
FIGURE 4.5
ASPEN PLUS SIMU1.ATION OF CHEMICAL PLANTS 193
B1 c>
C1
FIGURE 4.6
By the same way, interconnect remaining blocks. Renaming all blocks as well as
incoming and outgoing streams, finally we have the screen shown in Figure 4.7. To
rename a particular stream (or block), irst select it, then right-click, next select Rename
f
Stream (or Rename Block) and inally enter the appropriate name.
f
H _J iU _l J
kl
5","*,s Oirwu o ai
, SCfuc PMiofioc BMefiai Boictfuc
C |iFold».mo»IVilH HUM fw»llr»J
FIGURE 4.7
The status indicator in the bottom right of the window, shown in Figure 4.7, says
Required Input Incomplete indicating that the process flowsheet is complete and input
data are required to enter for running the simulation.
194 PROCESS SIMULATION AND CONTROL USING ASPEN
Configuring settings
As we hit Next icon and then click on OK, the following window pops up (see Figure 4.8).
Remember that in the Data Browser, we need to enter information using data input
forms at locations where there are red semicircles. As we inish a section, a blue
f
checkmark appears.
» r- -i'iv
-
life ,
FIGURE 4.8
V-H
-
id* ~ ' . -
figur: m 9
The next window (see Figure 4.10) includes the Aspen Plus accounting information ,
-
o- i S *
MM -
FIGURE 4.10
If we want the streams results summary sheet to display mole fractions select ,
Report Options under Setup folder to the left Under the Stream tab select 'Mole' as ,
n fT7
-
- .
i i
-
*- . -
-
»
u>U>i
.
2}' I -
n .
- .« <
FIGURE 4.11
196 PROCESS SIMULATION AND CONTROL USING ASPEN
Specifying components
In the subsequent step, use the Data Browser menu tree to navigate to the Components/
Specifications sheet. It is shown in Chapter 1 how to define components in the component
input form. Here, we have this table as shown in Figure 4.12.
g :.r.: Fopui
C2
CJ
C*
IU
cs
.
*
!S
*
_
j wy tm
1: - . .
_
J O pOa**
.
i i r
id
5 -.
_
| > M<MVnt«« -J
FIGURE 4.12
plant, set PENG-ROB base property method by scrolling down (Figure 4 13). .
I 3
J DM
.
3
PT
FIGURE 4.13
ASPEN PLUS SIMULATION OF CHEMICAL PLANTS -f 197
Properties under a particular model of Blocks folder in the list on the left and then
choose the suitable property method.
3
"
3 j. a. -j],. ,. j
2j twi- m*.
FIGURE 4.14
.
_
1 "
: ; £
it:
6 :
FIGURE 4.15
198 PROCKSS SIMULATION AND CONTROL USING ASPEN
Subsequently, the filled input forms are shown in Figures 4.16(a), (b), (c) and (d)
for other four DSTWU columns.
-
LT-
ap u**
t
neb*.
(a
Ct«p fo
1 F
-
Rww jo O
T
C3
O St*: G<t<-C-.
r>n»m
] .JtdiJ Trim
FIGURE 4.16(a)
J iiiJ JZi J
_
_
iJ~_
|0 IrcU 3Mt I HJilJ -3 U
_
j V- n CdaimxtOebor* j Convagervt
.
) .an d Ptuiite
i pi
a Wumbe.
- tjf
j
Reba» fT
_
.
U B3 r Reftw .aha
*
6-
-
_
1
.
U K
I J W
Ccrrp fic<
_
.
0 02 3]
_1 M
'
-
_
j W
j _, D£
Reow (0 CO)
S M Cl
jfi : .
tco- It* cm
f
RMI
t
P»u
I I Sob* I UMModrt t
k
Km-
Fv H* cm F"
FIGURE 4.16(b)
ASPEN PLUS SIMULATION OF CHEMICAL PLANTS 199
i _
J «
(7 ll**c»*
*
Zl r i »>-
S H
-
J '
-
J -
-
1 «
as
FIGURE 4.16(c)
Pta Jrav
fiimii
f
"
CondEnia f
u «
_
J w
_
J .
-
- I
FIGURE 4 16(d) .
Click on Afet and specify the pump (PI) outlet by providing the discharge pressure
r
FU Jtnry Wndov. I Ms
HA
L TJ i i
-
J2J
"
r Preci.ie "iciMie
r Pt«tiu>«(Mo
"
I I 3
) <Jfti«mr* (kthnjs condhoni
Saiw
ij»t SutrciJine
£0 MM
Sw&wjm -1
FIGURE 4.17
The status bar in the window, shown in Figure 4.17, includes a message of Required
Input Complete; it reveals that to run the simulator, sufficient data have been provided.
But there is no such restriction that we cannot specify the process with more input
information. Again, as we click on Next, Aspen Plus shows a message under the heading
of Required Input Complete as shown in Figure 4.18.
J T -
1-1. I pr H jsn l -iai
D>38
(? Ptmp
0 EWSK
j Rtcnc
"
"
3
3
P Result
_
J
I W ci
-
ai c;
j a o
.
ai «
.
a pi
EOCorwOpten
-
Q -
iM
. 9 -m *
FIGURE 4.18
ASPEN PLUS"' SIMULATION OF CHEMICAL PLANTS 201
Notice that if there are no red semicircles in the left it can be said that the data ,
.
1
B a
B O
B w
Z '
93 us states; ~. Jt mu tjxll huz - jjtkt:
n a = c» n
z. =i3»i arc:
fi" -
FIGURE 4.19
Viewing results
Choose Results Summary /Streams in the column at the left side and obtain the
compositions of all distillation products as shown in Figure 4.20.
We may save the work by choosing File/Save As/...using the menu list on the top.
W< tan give a name of the ile whatever we like. Note that if we click on Stream Table,
f
the results summary table is incorporated in the Process Flowsheet Window, as shown
in Figure 4 21.
.
all necessary information on the solved Aspen Plus problem including given process ,
3 >bJ.'<IP~~3»l-i g|
C- « 1 J .
~
J Cv CMn
Wit T a!
o . B -
am
"
TW
"
nrw
9 - l
Tia HH
"
i.r.T
e - - - -
n TTEB 1 m mid
" 1 '
e tUi sm 1 TOTB 1 nuif 1m muff
- ~ ~ " -
B nn- TW Mi
" '-
RD UB ' Wi 1 UK TW
"
-
.
4 2j
TS ~
a TBS
"
713 us 1 TB 1
r
S
" ____ ,
BB TRC -
IMl
"
irK. iv -.
tmufi *M ««
FIGURE 4.20
i MIBlJBI
-
1
.r
lOmm* I " I Ho I
FIGURE 4.21
Copyhghied malarial
ASPEN PLUS SIMULATION OF CHEMICAL PLANTS 203
voilM-CM otir*-!
SOLIDS 1MMCAMIC
CI C2M I
( 1 ClHt
1C4 C4H10-2 i
' ' C4N10-1 /
IC1 CSH12-; '
ici cvai-i /
MCt CftfU-1
LCMMKT
1LOW C» I«l-e4 O'JT-O) Bl
.LOCK C4 IWt} O'.n-.oi M
ukk ci in-c; out-03 el
FIGURE 4.22
43
. ASPEN PLUS SIMULATION OF A VINYL CHLORIDE MONOMER
(VCM) PRODUCTION UNIT
Problem statement
The process low diagram for Aspen Plus simulation of the vinyl chloride monomer
f
manufacturing plant is shown in Figure 4 23. The lowsheet has been developed based .
f
O-fcmi
O-feu
BB 66 B7
F10
[purge]-o
Pure ethylene, stored as a gas at 70oF and 1000 psia, with a flow rate of 20 tons/hr ,
and pure chlorine, stored as a liquid at 70oF and 150 psia, with a flow rate of 50 tons/hr
enter the mixer block Bl operated at 2 atm. The mixer outlet Fl then goes to the
reactor B2 run at 363 K and 1.5 atm. In this stoichiometric reactor (RStoic), the following
chlorination reaction occurs with 98% conversion of ethylene to 1, 2-dichloroethane:
C2H4 + Cl2 -> C2H4C12
ethylene chlorine dichloroethane
In the next, mixer B3 operated at 1.4 atm allows the mixing of the recycled stream
F12 with the reactor product F2. The outlet stream F3 is then condensed fully to liquid
phase in block B4 at 298 K before being pumped to an evaporator. The pump B5 has
discharged the liquid at 26 atm. The evaporator B6 performs the phase change operation
and then the vapour temperature is increased in the same unit to 515 K. In the
subsequent step, stream F6 is introduced in the reactor B7 (RStoic) in which the
following pyrolysis reaction occurs:
C2H4C12 -> C2H3C1 + HC1
dichloroethane VCM hydrogen chloride
The dichloroethane is converted to VCM and it takes place spontaneously at 773 K
and 25 atm with 65% conversion. To reduce carbon deposition in the heat exchanger,
the hot vapour stream leaving the reactor is quenched in block B8 yielding a saturated
vapour stream at 443 K. Quencher effluent stream F8 is condensed to liquid phase in
block B9 at 279 K and then fed to a DSTWU column B10 as stream F9. In the next ,
Stream F10 is introduced in another DSTWU column Bll. The first column mainly
separates HC1 from other components, while the second column purifies VCM from the
rests. Both the distillation columns have 10 theoretical stages (including condenser
and reboiler) and a total condenser along with the specifications shown in Table 4.3.
,
TABLE 4.3
Finally block B12 (FSplit) splits stream Fll to ensure the recycling of 99.999% of
Fll as F12 stream to mixer B3. A purge stream is introduced to prevent accumulation
of unreacted components.
Using the POLYSRK property method simulate the complete plant to compute the
,
Simulation approach
To start Aspen Plus package select Aspen Plus User Interface under Programs. When
,
the Aspen Plus window pops up choose Template and click on OK. In the next, select
,
Polymers with Metric Units (see Figure 4.24) and press OK button.
ASPEN PLUS SIMULATION OF CHEMICAL PLANTS 205
- . . r .ii.-i.*. Ti 3 . i i .rj 3 a
i : ; : I -i
-.
-
:
HZ]
FIGURE 4.24
Click OA" when the Aspen Plus engine window appears to obtain a blank Process
Flowsheet Window.
Creating flowsheet
We can develop the process low diagram (see Figure 4.25) by incorporating the following
f
ED
413
? 4f
FIGURE 4.25
206 PROCESS SIMULATION AND CONTROL USING ASPEN
f
several input forms one by one.
Configuring settings
After creating the lowsheet for the VCM manufacturing unit hit Next button followed ,
f
by OK to open the Setup /Specifications / Global sheet. In the following the irst screen ,
f
,
shown in Figure 4.26(a), includes the Title of the present project as Simulation of a
VCM Production Unit' and the next screen displayed in Figure 4.26(b), shows the ,
© SjtuIblco Cf
|SiT«iabon oi a VCM PieducUon Urd
itm&iB Gfcbal jetting,
_
O irouldars fwiT 3
Input mods
v Sf>ec/>awrn
Sueam das:
31
Flovbuc 3
it Jj Pew Oiar»a*fo
3
r Um see
£1 Cvre-ljdi
it; O Parw
5Tft£*MS
ssft.
FIGURE 4.26(a)
ASPEN PLUS SIMUI.ATION OF CHEMICAL PLANTS 207
Qi BI _
LJ *j«J .] r5h-|ftl»N|w| mj rj J_iilj _J J
I r -i-1-i Fir -
t lei, I i«l MhJ
Mi*
FIGURE 4.26(b)
We wish to have streams results summarized with mass fraction basis that is not
set by default. Accordingly, we choose Mass' fraction basis in the Report Options/
'
DIQ*IHI
'
J .
1 - M|
11 >id illF-Z]»J iJ _
P Hbi
I .1
FIGURE 4 27
208 PROCESS SIMULATION AND CONTROL USING ASPEN
Specifying components
The components that are involved in the monomer manufacturing process are ethylene
(C2H4), chlorine (CI2), 1,2-dichloroethane (C2H4C12), vinyl chloride (C2H:iCl} and hydrogen
chloride (HC1). In order to get a blank component input form, choose Components/
Speciftcatiojis in the left pane of the Data Browser window. Defining all species in the
Selection sheet, we have Figure 4.28.
510*101 _
Tii.J J
JT _1_M.
_
Setjr
Selection] Pelidftum | Noraiorivenliona! j-/Dalabanl s |
Soedicatona
0 Smjiaion Options Define comp«*nli
Q Slrean- Class Componerit ID Type Comconent name Foirtiula
I t-'enf - Comps
iMFJCGrtupi
'
1 Conp-Groops
.
J/] Comp-Ltsts
_
j Fo ert
PfOpefti«
| Stocks
'
I Corv-ergence
'
1 Flowsheeling Oplons
'
z1
if Miyeu/Splittets | Sepaialoi Heat Exchargeti ] Columns ) Reactots ) Ptestute Changei: ] Manpulatots j SoWs ) Usei Modeit ]
f
Matei.at
STREAMS ' »m FSpS SSpK
ForHe*j.pn»j Fl C\. gFtAlfla\Aip«nPliji 11 1 NUM fiML«<! tTua netwowK
'.
ami s - jjChapter4 1*00 | Oiapler I - Mero | T) Uidii Teimi I* | Aspen f\j» - VCM [| Aspen Plus « 9 C . '656
FIGURE 4.28
IM tp ... .' i ~ . .
f
CirftUl .-I i c|g| ni-«-|fcl%l-3M ml I M "I ?l -I _
Li
"
in fHMM
f
a
o
IT
VtBUKi
- -
. - M
FIGURE 4.29
|i*lB|
-
J L
"
3 "! M -1 ""1
10 Ci . I
|-«i.pr. 3 pw=-3P -3
IB 1 - f
fS5
ED
- "
iPin
.
-
d
1
.
_
j
. "
-
U "
.
3 U
5 "
-
"
. 3 .-
"
hi -
iTKMI 1 o rh. im
t
9
-
FIGURE 4.30(a)
ASPEN PLUS SlMULVnON OF CHEMICAL PLANTS 209
IM tp ... .' i ~ . .
f
CirftUl .-I i c|g| ni-«-|fcl%l-3M ml I M "I ?l -I _
Li
"
in fHMM
f
a
o
IT
VtBUKi
- -
. - M
FIGURE 4.29
|i*lB|
-
J L
"
3 "! M -1 ""1
10 Ci . I
|-«i.pr. 3 pw=-3P -3
IB 1 - f
fS5
ED
- "
iPin
.
-
d
1
.
_
j
. "
-
U "
.
3 U
5 "
-
"
. 3 .-
"
hi -
iTKMI 1 o rh. im
t
9
-
FIGURE 4.30(a)
'
t
t
_
i rj _
*
Pnpoty Mdhodi
Estimalcfi
d
Mo cJat Sttxicl'jie 1" F .
J ComponenI Va>je
t j Pafanrten cm
_
_
j Mt* iso P« d 1 J-DOl
PnapSel)
C2HXL
*
Jj y danced
J/J Streams Total How Ma;-. d HCL
-
C2H« 1 kgAr d
© IhdJ
RajuJs
d EO Vanables
-
(a CL2
0 Inpul
Rtsuls
Cf EO Vanables
.
_
j D1
.
_
j 02
.
J
Total IST
_
.
_
J "
t . , F3
t
S F4
. 1 F!
jl " Mixeit/Splitteis | Sepiialois | Heal Exchinaeti Cokjmu | Reactoi: [ Ptessuie Changei; | MarjpOaioii | Sold; j U;eiMc«Jd.-
Material
STREAMS 1 FSpia SS|*
FocHeto orevin C:\ gFoWere'-AspoiPlusIt 1
: NUM Rt»ju«ed ir<x£ TTC0ncJ«»t
3*1 S i Chaptera Aipen Plus - Smjali
FIGURE 4.30(b)
two reactors, four heat exchangers one pump, two distillation columns and one splitter. ,
Although discussed during the Aspen Plus simulation of different single process units
in the preceding chapters, we must remember the following points when we ill up the
f
To simulate a mixer model, at least provide the pressure data and valid phases.
In the simulation of the reactor model coefficients should be negative for ,
In the Vapour fraction ield of a heater model put 0 to indicate bubble point
f
and 1 to indicate dew point. For subcooled liquid and superheated vapour, use
Flash specifications.
The windows, shown in Figures 4.31(a) to 4.31(n) display the block-wise information
using the input forms.
ASPEN PLUS SIMULATION OF CHEMICAL PLANTS T
f* e* d*»
'
** * *****
ol lBj _
U M§) 1 I2iil*l5lll*d J£l 1 - I l"l ! gl l i$l
i I I JT I -imi LJll
-
jj
0
... -
B oi
' '
B u
I '"' d
VMphn. j 3
5 "
1i "
I 3
_
J "
H «
n 11 .
1
il a "
a "
-
a
& a fio
A a *
me
id Boot
0 H"
BDHxJ
Si v>'M>f ! plHMl i(53555B O fW pri- jre J value i 0 PretH/eiioo'vAe . Q G>uo» tnli Oudf (neuuefoc'Ae:
r
M- . FSpM SSj*
*
'
a*l| B«* | .
FIGURE 4.31(a)
a|tf|y| | I hi .tifssi il
i-nr mei. lai im
.
i sir i..e. -i oi i .
1 »l ol >l H
IwRMCtnni | Coitwum j Heal ol Re lun | 5«lectf*y } PSD | Componef An
.
U » Oc«>«cgcardl<or.i
.
a no
F n 3F
"
t a Fn 71
-
U fs 3
i, t*
f
a ei
O SoKGraup«
id
Bock Octon,
'Si* w»
FIGURE 4.31(b)
1"
212 PROCKSS SIMULATION AND CONTROL USING ASPKN
Compwienl
C2H4 i C2H4d2
12 i *
3
3
H Heck Ooi"-.!!!
Spec Group)
Si p
[i Hiwit/Spttlef* j iepatalow ] HealE«hsrwi | Column) | Beattot) ] Ptftsauro Chanijers j ManpjIalMi [ SoWs | UwMwtei?
FIGURE 4.31(c)
He b* Das
_
LJ~- I -1 i PT -I isi I ; m
3 jb [MET <<|
~~~
>>J _
Q vj N>]
~
VnaihOptwiwj
-
Mwo vecJicaMnt
VaWphaier 3
TtmpnaJue eiSmala Convtrjance paiamwlwi
~
3 Mainum toohons
E 3
-
Ji E2
Site CMC
J B3
amte
t
youwrheofenure AtwoUe i OjW p.enue K value > 0 Ptaw/e *«, rW 0 Ga yrj; Oj!*. p,. -* lo- 3h
'
MCorwca
STREAMS
PorHato-MiFl
Hea a. Heat.- MHaaC
-
Hwar.
~
HXFV* HrRCvlST J
4 UtoT jj Japw-J ctCTT:»3<
O
FIGURE 4.31(d)
ASPEN PLUS SIMULATION OF CHEMICAL PLANTS 213
0 c.-"
-
I DM
=4
'
"
_
j
O EO«>a
::
you Who mota. *ac«v ti*clon D 0 loc bubbte pwrJ 1 Olw (tew oo.nl Fm utorooled kjari o> wpcheaied vax* m l&Tpeia»u)fl sryj
iJ
wkfimitm
FIGURE 4.31(e)
:> t- .,
i
« «« ...
-
.
1 .
, ,
jin.
1 1
-
0
'
>
'y
0
.
J if
FIGURE 4.31(f)
214 * PROCESS SIMULATION AND CONTROL USING ASPEN
Raail
EOVa
EOW U 3 [sis r
Sc«ct
Psns |v*poi Irac'ion
S!>ew ReiJit
Selup
O P if ance Cur
B«k Oftior*
R«Jia
(f KVawbtes
51 EO ina/.
'
; ( Fens
~
f ] Stream ResJts
s ja be
© InpUl
"
: i Halves
© Cvrwrit
B*>d( Opteha
ResJls
J BDVafcK*
© EOkBU!
.
Reactws j Pie-jyjie Char eis ] Mawpife!«; ] SoWs j User Modei:
FIGURE 4.31(g)
Ems
'
Spec Gt Lps
J J Roil jTempeialuie
,
VardphasM
'.
] Hcurfes
'
© Bbqk Ostfcns
Pesiu
(f EO VanaWes
EC
© £pec Grmjps
Q Setup
3 optics
ResJia
d EO y*i*i«
l SpccGnjjw
f
£
STREAMS.
.
{>.< %} Mm FS,* JsJ
J| or ty. veil | - J «
FIGURE 4.31(h)
ASPEN PLUS SIMULATION OF CHEMICAL PLANTS 215
I r J- IT
(7-3
MOT j :«*«fl,
HO.
«
*
"
I 3
'
STREAMS Mm, FSct SSpt
an* on.Fl C-gFoMws axrPusIl 1 HUM
f
" 1 . . | - jChootw -'-o; j jdaotgl Ffeo j j Ussft Twrs W |[ Aapen Plus o. Fus '.CM | « 3C:&7
t
FIGURE 4.31(i)
IS
r c-... ._
r3r
V) !-
'J
1
'>
O
5 i
-
1
0 1
.
-
; -
i ii
"
[» Mao./Sclii.,, | j
j K A err-,
ilRUMS '
FIGURE 4.31(j)
216 PROCESS SIMULATION AND CONTROL USING ASPEN
-a
.
t:
-
Took Ru' PM Lbrv> Wlr
hJ J J
j r L _
_
IT l lal I lal ggj
FtothrpecAcMione
EObca
0 Sf«C Cir<XM
Vapo> risc'nn
.
|5 1 J
-
I p<,'1*
"
O 3]
O "
0 Bo* Octsrti
(f EOVautMt
_
J -
"
«i you type (UtWWKrroeiafuo S«H*
i Pwi
"
i
[i Hbeii/SpMtwt | .
-
cvy:.. | HeatEicch)ng»: | Ct-Unrj j ReflCioi)
STREAMS FSpM
"
FIGURE 4.31(k)
Fie Edl 4ew Dots Toott F>in Plot U«*v Wmdcw Hdp
us
i r i -1 i - pt i m i
Hanres *|
i
Resotr ~
1" 3 '>«fcme'
"
bd mill 3
Spec Graj-i
Pom
d Stwi. -mi
Tola* c vrferee.
o kpl 1
locov |0999
i
f
a Sec. Ccfcint
Cwm) |C2hj:l
I £0 Rectw (OOOI
rm
a Stream B«u»!
-
in
O
I9
sc.
Spec Gnue*
| ScMi ) LiMrHodM |
""MS
SIR£AMS ' M-t ts,*
FIGURE 4.31(1)
ASPKN PLUSIM SIMUI.VriON OF CHEMIC-VI PL\NTS 217
I >I = |h| I J
i r -i i i nr - i -m \ m
'
t-- i- F
f
1
O Si*. 3 .
f
« - I
1
Me
FIGURE 4.31(m)
i gr>i
..
'
Hi h* (te "t» fVr ty, ttai Widn >
1
re
FIGURE 4 31(n)
'
f
O BMkOpupni 1
SVt«r<
O So*: * 13 D B9SS
5 N" * -
URGE
V) B11
'
Rj in if.- rcn'
a
Rw-ij 01t 1 1
t*
f
SbMM BenA
lonnt/SplilUK |
FIGURE 4.32
As we hit O-K button on the message the Control Panel window appears as displayed ,
gaBSBEsai
*) £ Cm IceM umj WrDo* m» 4uS
J £
'
» .» %w,cj - lilSSEi
FIGURE 4.33
ASPEN PLUS SIMULATION OF CHEMICAL PLANTS 219
Viewing results
Choose Results Summary /Streams in the column at the left side and rearrange the table to
get the results in the form as shown in Figure 4.34 Save the work positively at this moment. .
_ _
r -i i I'-f -
i HOI
jaTS i »M«Hr- ..|.|_.|n.|
il 1
. .
'
Taw- >fa»gi
O
W - TTOoo EM I nsr i-1 t&nti
-
B
1
' ' -
io)A«tf j lb. 4« 4
1
(Mt
-
SJjSTS-
-
37 3 ;3 : i
i axns ob (SB? ? 00 00
| BU wn
"
51 H
»»r . "IT 1583
IMHM * 004!% KUOEBE )<UCCQ nStWi i
1 <«U< I UVMMW ]
FIGURE 4.34
sntew CiiM
."
n -ncw cl2 y>. ctorn'tr*
FIGURE 4.35
220 PROCESS SIMULATION AND CONTROL USING ASPEN
approach to simulate a lowsheet is that after developing the process low diagram in the
f
f
lowsheet window of Aspen Plus, we can simply use Next button for data entry. As we
f
receive the message of Required Inpu t Completey we can move on to run the simulation.
In the next two chapters, we will study the process dynamics and closed-loop control of
low-driven as well as pressure-driven processes using Aspen Dynamics package.
f
PROBLEMS|
4 . 1 A hydrocarbon stream with component-wise low rates, shown in Table 4.4, enters
f
the isentropic compressor at 120oF and 1 atm. The compressor has discharged the
vapour stream at 3 atm.
TABLE 4.4
10
95
150
n -C4 25
/i-C3 10
n-C 100
6
The complete process lowsheet for lashing and stripping operation is shown in
f
Figure 4.36. The lash drum (Flash2) runs at 1250F and 2.8 atm. The stripper
f
(STRIP2) has total 6 stages (including condenser and reboiler) and bottoms to
feed ratio (mole basis) is 0.8. The feed stream to the stripper is introduced above
the top stage and the pressure throughout the column is 2 atm.
V24
I 0 1
CopyHghlod material
ASPEN PLUS SIMUIJUION OF CHKMICAL PLANTS 221
Using the UNIQUAC property method, simulate the plant to compute the product
compositions and flow rates.
4
.
2 A ternary mixture, as shown in Table 4.5, is fed as stream H at 100oF and 290
psia to a pump Pi employed to increase 20 psi pressure.
TABU 4.5
500
300
"
r 11, 10
The stripper (STRIP2) has total 100 stages (including condenser and reboilen
with a reboiler duty of 107 Btu/hr Stream F enters above 70th stage and Stream
R) mien above 1st stage. The top stage pressure of the stripper is 280 psia
with a stage pressure drop of 0.5 psi The intercolumn pump P2 has increased
25 psi pressure The RECT column has total 120 stages (including condenser
and reboileri with a reflux ratio (mole basis) of 10 and a bottoms to feed ratio
'mole basis) of 0.6. Stream Dl enters below 120th stage. In the simulation,
consider condenser pressure of 275 psia with a pressure drop of 5 psi and a
stage pressure drop of 0.1 psi (see Figure 4.37).
0
-
0
FIGURE 4 37
.
A flowsheet of a propylene-propane mixture separation process
distillate
222 PROCESS SIMULATION AND CONTROL USING ASPEN
FA
S-o
C>-I ANILINE
El
PUMP
CD-
-
CSTR
COMPRESS
The reactor model (RCSTR) operates at 580 psia and 2480F, and its volume is
1200 t3 (75% liquid). For the liquid-phase reaction, the inlet streams have the
f
Both pump and compressor (isentropic) have discharged the fluids at 585 psia.
Data for the Arrhenius law are given as:
Pre-exponential factor = 5x 105 m3/kmol . s
Activation energy = 20,000 Btu/lbmol
[CJ basis = Molarity
Use the SYSOPO base property method in the simulation. The reaction is irst-
f
order in aniline and hydrogen. The reaction rate constant is defined with respect
to aniline. Simulate the process and compute the component mole fractions in
the liquid product and the vent stream.
4 . 4 The process flow diagram for an azeotropic distillation process is shown in
Figure 4.39. The technique involves separating close boiling components
(acetic acid and water) by adding a third component (vinyl acetate), called an
Copyrighted material
ASPEN PLUS SIMULATION OF CHEMICAL PLANTS 223
entrainer, to form a minimum boiling azeotrope which carries the water overhead
and leaves dry product (acetic acid) in the bottom. The overhead vapour is
condensed and then separated in the decanter into two liquid phases: the organic
phase and aqueous phase.
DECANTER
1 VA-RICHh
DIST
VA
1 FEEDf
HW-RICHh
AA
RADFRAC
A feed stream, namely FEED, enters above 15th stage of the azeotropic distillation
column at 330oF and 90 psia in addition to the low rates, shown in Table 4.7.
f
TABLE 4.7
The entrainer, VA (vinyl acetate), with a low rate of 455 Ibmol/hr enters above
f
12th stage of the column at 200oF and 100 psia. The azeotropic column (RadFrac)
has the following specifications:
Number of stages (including condenser and reboiler): 55
Condenser type: total
Valid phases: vapour-liquid-liquid
Reflux ratio (mole basis): 4
Bottoms rate: 2700 Ibmol/hr
Condenser pressure: 66 psia
Column pressure drop: 12 psi
Key component in the second liquid-phase: water
Stages to be tested for two liquid-phases: 1 to 55
The specifications for the decanter model are noted below:
Pressure: 50 psia
Temperature: 110oC
Key component in the second liquid-phase: water
Using the NRTL RK thermodynamic model simulate the process to compute the
-
4 . 5 A hydrocarbon stream H is at 50C and 2.5 atm. The pump has discharged the
liquid feed F at 5 atm. The component-wise low rates are shown in Table 4.8 for
f
stream H.
TABLE 4.8
35
C3 50
i-C, 130
n-C 200
4
'-
c5 180
n-C5 200
n-C 5
6
DRl DR2
CR2
CRT
BR2|-0
BRI
G1
PUMP
h
C2
C3
B3 | C1 B2}<>
TABLE 4.9(a)
TABLE 4.9(b)
CR2 Total 05 .
15
the use of an isentropic compressor Bl. Stream H has component-wise flow rates,
shown m Table 4.10.
TABLE 4.10
mtrogen 100
hydrogen 300
ammonia 0
carbon dioxide 1
A low diagram for the ammonia process (Finlayson 2006) is shown in Figure 4.41.
,
f
B '
m-<
H-
0 H"
ED-0
Bl B2 B3 84
f
block B3 to 900oF at 4000 psia. Note that the reactor (RGibbs) B4 runs at 900oF
and 3970 psia. In the next, the reactor outlet F4 is cooled in a heat exchanger
B5 operated at 80oF and 3970 psia. The flash drum (FIash2) B6 produces Streams
Bl and F6 at 80oF and 3970 psia. In the subsequent step, Stream F6 enters the
splitter (FSplit) B7 and 0.01% of it is used as purge. Finally, an isentropic
compressor B8 has discharged Stream F8 to the mixer block B2 at 4000 psia.
Using the NRTL thermodynamic model and the Newton's iteration method (from
the Data Browser, choose Convergence/Conu Options), simulate the ammonia
process to compute the component-wise low rates and compositions of all streams.
f
REFERENCES |
Finlayson, B.A. (2006), Introduction to Chemical Engineering Computing, 1st ed., Wiley
Interscience, New Jersey.
Seider, W.D., J.D. Seider and D.R. Lewin (1998), Process Design Principles: Synthesis,
Analysis, and Evaluation, 1st ed., John Wiley & Sons, New York.
Copyrighted maierlal
Part III
51
.
INTRODUCTION
used or process design and control. This powerful simulator can automatically initialize
f
the dynamic simulation using the steady state esults of the Aspen Plus simulation .
r
Interestingly, when the ile containing the lowsheet is opened in Aspen Dynamics the ,
f
default control structures are already installed on some loops Usually, level, pressure
.
and temperature controllers are included where appropriate However these default ,
control schemes can be modified or even eplaced with other suitable control loops
r
available in Aspen Dynamic- package Note that there is a scope to include some
additional controllers for the used process Moreover this simulation tool provides a
,
items that should be taken care of For example the size of all equipments must be
,
specified and the control structures must be devised For steady state simulation using
Aspen Plus the size of the equipment is not needed, except for eactors. On the other
,
r
hand ,
for dynamic simulation using Aspen Dynamics, the inventories of material
contained in all the piece* of equipment affect the dynamic response Therefore, the
.
driven dynamic simulation Pres ure-driven simulations include pumps and compressors
-
.
where needed to provide the required pressure drop for material low Control valves
f
must be installed where needed and their pressure drops selected For flow-driven
,
229
230 PROCESS SIMULATION AND CONTROL USING ASPEN
In the present chapter, we wish to study the dynamics and control of the low-
f
driven processes. For this intention, we choose a reactor (RCSTR) as well as a distillation
column (RadFrac) example rom the model library of Aspen simulator.
f
52 .
DYNAMICS AND CONTROL OF A CONTINUOUS STIRRED TANK
REACTOR (CSTR)
Problem statement
Ethyl acetate is produced in an esteriication reaction between acetic acid and ethyl alcohol.
f
acetic acid + ethyl alcohol ethyl acetate + water
A feed mixture, consisting of 52.5 mole% acetic acid, 45 mole% ethyl alcohol and
2 5 mole% water, enters the RCSTR model with a low rate of 400 kmol/hr at 750C and
.
f
1 .
1 atm. The reactor, as shown in Figure 5.1, operates at 70oC and 1 atm.
Both the reactions are first-order with respect to each of the reactants {i.e., overall
second-order). For these liquid-phase reactions, the kinetic data for the Arrhenius law
are given below:
Forward reaction: A = 2.0 x 108 m3/kmol s
S = 6.0 x 107 J/kmol
Reverse reaction: k = 5.0 x 107 m3/kmol . s
£ = 6.0 x 107 J/kmol
Diameter: 0,45711 m
Volume: 0.15 m3
(a) Simulate the reactor model using the SYSOP0 thermodynamic model to compute
the product compositions.
Copyrighied malerial
)YNAMICS AND CONTROL OF PI.OW-DHrVKN PROCKSSKS 231
(b) Report the default controllers tuning parameters and control actions used, and
constraints imposed on variables.
(c) Investigate the servo performance of the default liquid level and temperature
control algorithms and discuss the effect of loop interaction.
(d) Show the regulatory behaviour of both the controllers in presence of disturbance
in feed temperature.
Simulation approach
(a) To open the Aspen Plus Startup dialog box. click the desktop Start button, then
point to Programs, AspenTech, Aspen Engineering Suite, Aspen Plus Version
'
and then click the Aspen Plus User Interface. Let s select the option with
Template and then click OK (see Figure 5.2).
o-tt rt m mots; am -
bHft mm 1 ; :_ tm
FIGURE 5.2
As the next window appears (see Figure 5.3), it is appropriate to select General
with Metric Units and hit OK button.
Here we use the simulation engine at 'Local PC When the Connect to Engine dialog
pops up (see Figure 5.4), press OK. Note that this step is specific to the installation.
Creating flowsheet
The process low diagram, shown in Figure 5.5, includes a reactor, namely RCSTR,
f
Copyrlghiod material
232 PROCESS SIMULATION AND CONTROL USING ASPEN
Nil Zl ,
iJ_l
_
J J
I I 1 1
r
.
.
3*3pe IPE S!fe»T Property kit Fhafnieceul cals hI Caenerai SmJalion vrth
fBark Smutatcn L Fhamaceuticala Mettle iJnts
r
t
MMKca(/W cum/nr
i3 Chemicals w<h Heine Unrts Polymer wrO". Mibfe
j BMMMM »**> BfSlW> UnJs f. P,TT>fr«*atu / /.-eh Property Method Nooe
i S cl VM **h MeMc Lints £jj P>n)fr«(arur3y /. (h
FtewbwfoTB xi Mole
i jGas Process iMi English UrAa So!'*;-fh EngWi I
r
yScWswlh fAetrcU. Strearr report
,
G*n*rai wth English IJhta La Specia*/ QmiOk ccmpoabon Hole flow
"
£2 Specialty ChemV;a*s
d m alkjgy EntfishUnfts
dumelallurgy Mih Metnc Unfis Run Type
1 2] [Rwsheel 3
111
Fot Help prws Fl
,
FIGURE 5.3
Connect to Engine
User Info
Node name:
User name:
PassiAiord:
Working directory:
OK Exit Help
FIGURE 5.4
PVNAMICS AND CONTROL OF FLO\V DRl\-EN PROCESSES 233
M rn,-» 1
.
« I I ---
-
- -;
e i i i ysa
'
-
-
.
i
FIGURE 5.5
Configuring settings
Hitting Next button, we get Global sheet of the Specifications form under Setup folder
in the left pane of the Data Browser window. Enter the Title of the present problem-
Dynamic Simulation of a CSTR'. change the Input mode from 'Steady-State' to "Dynamic"
"
and leave the remaining items at their defaults. The window looks like Figure 5.6.
? ).!
,
81 V
.,
* a
jLr
~
i
3 I 1 !-! -'in: q >!i.
4
-
ia
Km "ttm "CT on
js- b
FIGURE 5.6
234 PROCESS SIMULATION AND CONTROL USING ASPEN
In the next window, as shown in Figure 5.7, the Aspen Plus accounting information
required at some installations are provided.
Hi _l 3 J J
_
j r _
I I I [»" -
I M I M
Accouxlrurbar
© Mm >-
O R»MOttor# PrcrKtraw fr-OUH WISH
i -
lnoMtCflinpi<»
[V Hm>iA(«m j $etur*» j --v . I C(*ftm RMEton | PreiMeChange | Marv ai j Sokb | UwMoiw |
'
FIGURE 5.7
In the subsequent step (see Figure 5.8), select Stream sheet with opening the Report
Options form under Setup folder and include Mole fraction item.
.
tfe E* V1e» 0*4 '<xk Rn tM Ufmy WMoh -.
f
1 hi .JSlal _J
i r-1 i i nr i -m \ |a! «|».|
-
3 »| j ,| H.|
P Gtf«<a(e»ii»xlsrJ*Mn-feKirt
am'o br rcUM r lDcmti itpvt
t
iWhua Fi«cMntt*m
P Malt P Mole TIT fGEfTM
0 RVMlOIMm
r Urn* r Mm.
.
_
l -r. - - -.
.
Jj
-
CH
- s -
FIGURE 5.8
DVN VMICS AND CONTROL OF FLOW-DRIVEN PROCESSES 235
Specifying components
In the Data Browser window, choose Components /Specifications to obtain the component
input form. Filling out the table with the components (acetic acid, ethanol, ethyl acetate and
wateri involved in the present reaction system, the screen looks like Figure 5.9.
MM - . -
E-sC
aun
i 1
I
I*
, ,,.
.
* - -« MUM
FIGURE 5.9
'
-
KW -U MSI *) nKI*M-aH 1 3)
9 «-.*«
. 2
FIGURE 5.10
236 PROCESS SIMULATION AND CONTROL USING ASPEN
oj -
u l-w
-
-
3 ****
* 1
-
- 3
3 (.-. j| 3
T 3
('- * z
h -
i= 3 -
.
3
pr (n -%\
I J
I "in
FIGURE 5.11
-
o -
# l-s-g-Q-m-g-
-
. j
'
T*4m ! . » Iff ..-*-. |.;a .l»-. |<j w> r
FIGURE 5.12
DYNAMICS AND CONTROL OF FLOW-DRIVEN PROCESSES 237
In the next step (Figure 5.13), select RCSTR I Dynamic! Vessel sheet under Blocks
folder and enter the reactor geometry data.
LMf ***** -**
JUIBI U _
51 nMi>H<H 3 I H .-J JJ -1 J
i r i i i fv i lei I Pi 3>M
. i
FIGURE 5.13
The forward reaction as well as the backward reaction is represented with their
stoichiometric coeficients and exponents in two sheets shown in Figures 5.14(a) and (b). ,
f
EAC
_
l ***
gff Vftrtc:
ll
-
CH-
.
1-0 LIZ U
FIGURE 5.14(a)
238 PROCESS SIMULATION AND CONTROL USING ASPEN
f
_
i r i i i nr i m I issi
71 f. .1 I .1 I .<lr:- 7i »i nl «i uf\
-
UMFAC Grows;
Reaction No fyT
Readmit Pto<k«<t
1 V
-
J 0"
Prop Sett
3 -"EEC
j PRODUCT
_
3 Becks
.
J/j 3CSTR
_
j Ow stry
-
Resctcr:
w R-l
J t eti e
ll -
FIGURE 5.14(b)
The power law data for both the reactions provided in the problem statement are
entered in the two Kinetic sheets shown in Figures 5.15(a) and (b)
ar -i MM <<||AI 3 »l o\s
~
/Sto-Jwnwy *Kinebc| : - |
'
; | He«yG)mt>s
UfUF Grauo*
}l) AA- EAL-) EAC . VATER d
iwc Cte, | . -J
J/i Zcrvij&
f
;
i "
k. 200000000
[ W
.
.
. a P«* en i E::
_
l 0** To-
I
(qi«w.
> ) FEED
£ PRODUCT
.
1 SCSTR
y Reacscoi
_
J Ow-wy-
-
nn miiii
SIRtAHS
J- 90
FIGURE 5.15(a)
DYNAMICS AND CONTROL OF FLOW-DRIVEN PROCESSES 239
-.- 1 3
I I
r
.
«
st e
-'
FIGURE 5.15(b)
The status indicator in the above window reveals by the message Required Input
Complete that no more input specifications are required to run the simulation .
dialog box appears Hitting OK on the message, we are displayed the Control Panel
.
where the simulation messages during the run are recorded (see Figure 5.16).
- ?!
- MM
FIGURE 5.16
240 PROCESS SIMULATION AND CONTROL USING ASPEN
Setup
1 I I 1
Omponeras
Streams
3 foumi r
" ~
as Bocks
Displw 3 Stie«mTatile|
Reacts
3fl Ccr.v«gefK«
'
1 Conv Options
41 d1 d
£0 Conv Options Volume Flow cumyhi 24.497 24 01S
O S«up Enthalpy MMkcal/hi .
35035 35 951
o 0M0 Bsslc
o OMOAdv Mole Flow km*hi
o lSSQP Base M 2ia0CB 87110
0 LSSQP Ad*
EAl 180.000 57110
Tear
EAC 122,890
1 Convei ence
Conv Ortef WATER mooo 132.890
Q| Sequence id* fan
Rowsheetng Ctouons
AA 0 525 0 218
Model Analvsis Too .
Convetgence
Jdll
Results Avaisse
Mstets/SpittM! | Seoatatois | Heat Enchangeis \ Cokmra Hoactoti | Pressuie Chatlseis ] Minipulatois ) Sofds | Usei Models |
| j£]Oaiaer5-McmsollW..| 4]Cha(»er2-Wt!roseilW...|| Aspen flu. - SIimSI.. Adohe olal Pntesai |« iS?t}); 1545
t
FIGURE 5.17
(b) Exporting dynamic simulation: In the subsequent stage (see Figure 5.18),
we wish to carry out the simulation dynamically. Accordingly, at this moment,
we have to follow the sequential steps noted below:
Click on Export from the pulldown File menu or simply press Ctrl+E on the
keyboard.
Open the Drive and then Folder where we want to save the work as a ile.
f
Choose 'Flow Driven Dyn Simulation (*.dynf & *dyn.appdf)' rom the options
f
the same folder within which the exported dynamic simulation ile is saved. Originally
f
many iles are saved along with the backup or dynamic ile. Anyway, we are now ready
f
3 ±ljili<JF; »l jlalH
. I I I
.
1 , 1
-
9 Ljii=v.
Tt*-
CJ Omwrngnt **, MW
Can***
EAL
WATER
SrniatTi r irTi i -o/r «Caned |
Jjj
8 . S -O'll-O-
cvi
| 0>»l»5-rfco...| *]am2-«*n,. \.*]Mll<m,-lic..\\ <»«,PIU.-C , '; eJu t |'« igsj
t
f
FIGURE 5.18
Suite, then Aspen Dynamics Version and then click on Aspen Dynamics a blank dynamic ,
D H SQ W |Sl»dySl.lc j-J
3
r «-Jt c ioos j a, k it
Jaw
.
FIGURE 5.19
242 PROCESS SIMULATION AND CONTROL USING ASPEN
f
press Ctrl+O on the keyboard. In the Open dialog box, locate the drive, then folder and
inally the ile 'Ch5 52 . RCSTR' (see Figure 5.20).
f
_ _
ij is H *»Q IS
UIMl
r »-
-
SSmJatn
> Ft OinMITKt
£l Hi
5 *«» ICM .
U f
.
. ITS | CWn |
Opin
B i -r
FIGURE 5.20
inserted level (LCI) and temperature (TC2) controllers appears (see Figure 5.21).
-
Li
Tt MM [MM :"
u
'
i
t- - - - .
-Ill
f
. .
. ..
..
Li I
m
lad *>.
in
FIGURE 5.21
DYNAMICS AND CONTROL OF FLOW DRIven PROCESSES
-
243
Details of the two control loops, to be used inally, are given below
f
.
Loop 1
Controller: LCI
Type of controller: proportional (P) only
Controlled variable: reactor liquid level
Manipulated variable: product flow rate
Controller action: direct
Loop 2
Controller: TC2
Type of Controller: proportional integral (PI)
Controlled variable: reactor temperature
Manipulated variable: heat duty (cooling operation)
Controller action: reverse
Note that the direct acting control system increases the output signal as the input
signal to the controller increases. On the other hand, as the input signal to the control
structure increases, the output signal from the controller must decrease for the case of
reverse acting control strategy. The direct acting control law has negative gain and
increase/increase (or decrease/decrease) term is commonly used to represent it For the .
reverse action, increase/decrease (or decrease/increase) term is used and controller gain
has positive sign.
The reactor flowsheet includes two (LCI and TC2) single-input/single-output (SISO)
control loops. Therefore we can say that this is a multi-input/multi-output (MIMO) or
,
the set point is represented by SP and the controller output or control variable or
manipulated variable is abbreviated by OP .
For the example CSTR system level and temperature controllers are automatically
,
implemented when the Aspen Dynamics simulation is created The default values for
.
SP , PV and OP are computed from the steady state simulation. To achieve better closed-
loop process response the Aspen-generated control structures can be modified or even
,
replaced by the suitable control schemes available in the control library of Aspen
software .
In addition, the default values for controller tuning parameters, such as gain,
integral time derivative time and so on, can also be changed.
,
,
Most of the control strategies are easily tuned by simply using heuristics. As
suggested by Luyben (2004) all liquid levels should use P-only controllers with a gain
,
of 2. All flow controllers should use a gain of 0.5 and an integral time of 0.3 minute
'
also enable iltering with a ilter time of 0.1 minute). The author also mentioned that
f
the default values in Aspen Dynamics for most pressure controllers seem to work
reasonably well But temperature controllers often need some adjustments.
.
,
244 PROCKSS SIMULATION AND CONTROL USING ASPEN
Figure 5.22, is same with the value of process variable. To show the results table of
loop I, highlight the controller block LCI, press the ight mouse button, go to Forms
r
and then select Results.
'
-
i .mi m i -- w*
ffes
El &
m £
11 ..1 in
FIGURE 5.22
We can have the same information in a faceplate, shown in Figure 5.23, simply by
double-clicking on the block LCI. But as a difference, the units are not mentioned here
with the values of SP. PV and OP.
.
. ;i J
-
.
W ,1V- 3
>
--
IP
mar
urn
FIGURE 5.23
Gopyngt-
DYNAMICS AND CONTROL OF FLOW DRIVEN PROCESSES -
245
Similarly, we have the results table, shown in Figure 5 24 for the temperature loop 2. .
,
E3B i«ra
-
-
frol-
r-O-
HS
FIGURE 5.24
r,.-:n .
O j ct
-
FIGURE 5.25
246 PROCESS SIMULATION AND CONTROL USING ASPEN"1
Alternatively, to obtain the Configure dialog box, irst double-click on the controller
f
block LCI and then click on Configure symbol (yellow colour) in the faceplate as shown
in Figure 5.26.
C SsD-
U 1
.
3
«liBtn»ted becAutA t . hsd residuals ouer le-OUS
Ml
FIGURE 5.26
By the similar way, we obtain the tuning data sheet, shown in Figure 5.27, for the
temperature controller TC2.
o i
u
a» F f.A
(MMM |a FT-3
£3
r Dmo>
d5.
FIGURE 5.27
DYNAMICS AND CONTROL OF FLOW-DRIVEN PROCESSES 247
Note that the default Operator set point value is the steady state value of the process
variable (PV). The reactor liquid level is the PV for loop 1 and reactor temperature for
loop 2. Bias signal is the output rom the controller when the error (= SP-PV) is zero.
f
From the results tables shown earlier, it is obvious that the error is zero for both loops.
Therefore, Aspen Dynamics has set the value of OP as the bias value.
The proportional integral (PI) control methodology is automatically installed with
default values for the controller gain (= 10 %/%), integral time (= 60000 minutes) and
derivative time (= 0 minute) to monitor the reactor level. However, as mentioned
previously, the proportional-only controller with a gain of 2 is sufficient to effectively
control the liquid level. Remember that to make the integral action inactive, we can
use a very large value, for example 105 minutes (even the default value of 6 x 104
minutes may also be accepted), for the integral term. For loop 1, the controller action
should be Direct' as set by default (see Figure 5.28).
'
Baa
Ffc MM Took Af) Wltttm H*
|DynBfnic ]
"
U 7f 3 ii t i*
f
BLOCKSfTCn Co mv
r : i
Bw |21083 6171
Gat
jiooooo d
L
m m ConOofc* «aon -
-
-
- -
-
t cj
f
f
Integrftting tro. 62 43 to 62 44
Integrating tr » 62 44 to 6? 45
Mograt ing fro* 62 45 to 62 46
Step 3074 Tiw 6 2504«+001 step size- 5 OOOOe-002 step facto 1 5000e+000 accepted
j
wm -.
FIGURE 5.28
iDynamlc ii m * 93 O ci rj a
f
r ft -jt G,.t| 3 iL Si i 'i
f
5
Ope>«ioi ;a>pan p 1
Mi
G«n h
tmegral dm --
w-lPRoSucrt-H .
,
.
['' . v- tme I
C.>**i Chwi
r Died
Revwte
IfWiataeVabei | Htto
m
alftulation ready for solution 3
equations were not eliMinated because they had residuals ever le-OOS
k total ol 109 equations (29 S'-i) were eliminated
Simulation has 363 variables 260 equations and 968 non-zeros
lAjpenOrn-c.
FIGURE 5.29
D B # Gi ©i ? jDynamic I) M «
Q SniMoh
a Bi.ocKsnrncorfKi
Piocott voiabie
"
Ranjsmawwm (l 8280M (m
o
TC3
Ct>0>n btT»y
.us S- Ktr
Ootore
M ty«m>ce«a
I UCh e'
FIGURE 5.30
DYNAMICS AND CONTROL OF FLOW DRIVEN PROCESSES -
249
As shown in Figure 5.30, the default ranges for both the process and output variables
are too large (± 100% of the steady state values) It may be practical to consider the .
following constraints.
Process variable
Range minimum: 0.6855 m (25% subtracted from steady state value of PV)
Range maximum: 1.1425 m (25% added with steady state value of PV)
Output
Range minimum: 15812.7 kg/hr (25% subtracted from steady state value of OP)
Range maximum: 26354.5 kg/hr (25% added with steady state value of OP)
Entering these upper and lower bounds, we have the window, shown in Figure 5.31,
for the level controller.
Ul«l
ftv Ut» Wnd?* He«
f
r t 7t eiloos J ,w ii 5V s-: ;
t
[M
h
jl M253G h PBODUCT)-C|>
Oulf/
S firty. BUMMM
|1SBI27
91
l-wtions <"-ie nol elikinst«d Ivrrsuse they had residuals over le-OOS
h total ol 109 equations '29 S;| vete .Uaiaat«d
StTOlation has H'i v«ri*&lM 260 aqtrnttoM end 989 iior-:eros
Vital iV«nw:*00(
FIGURE 5.31
Again the typical ranges for the temperature control loop are noted here.
Process variable
Range minimum: 52 50C (25% subtracted)
.
Output
mention that the negative value of heat duty reveals the cooling operation (heat removal).
250 PROCESS SIMULATION AND CONTROL USING ASPEN
o *y .a m w |n,..mic
, .
'
9
Editing Siaulacion
Validation coBplete
FIGURE 5.32
Both the control algorithms are completely specified above. In the next, the controller
performance will be examined in terms of set point tracking (servo) and disturbance
rejection (regulatory).
(c) Starting the Run: Before running the program, we must be accustomed with
some frequently used items of the toolbar as described in Figure 5.33.
. IMJSI
FIGURE 5.33
DYNAMICS AND CONTROL OF FLOW DRIVEN PROCESSES -
251
We wish to carry out the simulation for a certain time, say 5 hours To ix up this .
f
time period, select Pause At from the Run pulldown menu or simply press Ctrl+F5 on
the keyboard. Then select Pause at time, type 5 in the ield or whatever we want and
f
click on OK (see Figure 5.34).
yrr&ntmsmfmsssi : .
,,. 1.1
a uai
,
, s - MA
c; |W)-
1 « 1 : O -
tC5tR
f
S i 1
.
i ;;;.-,;"',:::/;;:,:'m;.;:'::;f iv,i:
,
J
FIGURE 5.34
next, press on Configure and Plot symbols in the faceplate. Alternatively, to open the
faceplate Configure dialog box and ResultsPlot dialog box, first select LCI block, then
.
choose Forms and subsequently press one-by-one on faceplate Configure and ResultsPlot, ,
k* h s a © w [n n 3 B9t> esaspa
< {MB
n mm MUUI
ma : -
.
ja_!J
nMm I fmg, | <
! r.j J
-
>.,,
pi KJS» ,
|r~ 1"
l
~
1
-
,
1 «, 1
I ogji I at.
t
FIGURE 5.35
252 4- PROCESS SIMULATION AND CONTROL USING ASPEN
First make sure that all the items in the Configure dialog box and faceplate are
correct. In order to execute the dynamic closed-loop simulation click on Run button in ,
the toolbar. During the simulation run, give a step change in the set point value of
reactor liquid level from 0.914029 to 1.1 metre at time = 1 5 hours. Typing the new set .
point value in the faceplate, press Enter button on the keyboard so that the Operator
set point value in the Configure dialog box also changes automatically to 1 1 . meter
Note that the new set point must be within the specified ranges of PV In Figure 5 36 ,
the servo performance of the level controller is depicted for 5 hours as selected earlier
Obviously, the plot also includes the manipulated input profile .
EC
DlSySa IE -V | Dynamic 3
* ifl)_l
T tt Tl' ©id 10 05 it
! Rowsheet
Coopoient Lots
} Dagncsl'cs
*
SS
t'POC I
IT
o JDJ*]
t
|21083.61 RCSIR
: lrJe.;tdlme: limtjo |
m
I" l d
-
If5 11
'
Cpndolei adkm
: Direct
-
( Reveise
||
:o» 4 97
-
obi 4
Paused
FIGURE 5.36
Figure 5.36 represents an excellent set point tracking performance of the level
controller (P-only). Obviously, the LCI provides process responses with almost no
deviation from the desired set point value and with very fast approach to reach the
target liquid level.
Notice that the above plot can be edited by right clicking on that plot and selecting
Properties option or by clicking on that plot and pressing Alt+Enter on the keyboard.
In the properties window, user can modify the title, axis scale, font and colour of the
plot. Alternatively, double-click on the different elements of the plot and modify them
as we like to improve the clarity and overall presentation.
Now, we will discuss the interaction of two control loops. When we introduce a set
point step change in the reactor liquid level, the LCI scheme attempts to compensatell
for the changes through the manipulation of the efluent low rate. This, in turn, wi
f
t
disturb the reactor temperature and loop 2 will compensate by manipulating the hea
DYNAMICS AND CONTROL OF FLOW DRIVEN PROCESSES
-
253
removal of the CSTR appropriately. Thus we can say that loop 1 affects loop 2. In Figure 5.37,
Aspen Dynamics window demonstrates the loop interaction under the same set point
step change (0.914029 to 1.1 metre at time = 1.5 hours) as considered previously .
PS
~
rKjTj F= E],. ;- %
c;
o
.
BEDS
1 - 3
J _JSS_J
FIGURE 5.37
'
ehj
f
FIGURE 5.38
Now we wish to conduct the servo performance study for the TC2 controller with two
consecutive set point step (pulse input) changes in reactor temperature (70 -) 750C at
time = 1 2 hours and then 75 -> 70oC at time = 3 hours).
.
254 PROCESS SIMULATION AND CONTROL USING AS PEN
Clearly, the proportional integral controller with default tuning parameters values
shows a high-quality temperature tracking performance. As stated if the performance ,
of any controller is not satisfactory, we have the option to tune the parameters simply
by trial-and-error method.
If we introduce a set point change in the reactor temperature the TC2 controller ,
takes necessary action with adjusting the heat duty to compensate for the changes .
But interestingly, the liquid level remains undisturbed. Figure 5.38 confirms this fact .
At this point we can conclude that loop 1 affects loop 2, but loop 2 does not affect loop 1 .
To change the feed temperature twice as prescribed above, irst we need to open
f
the feed data sheet by double-clicking on the FEED block in the process flowsheet
(see Figure 5.39).
-
Ch55.2 BC5rRd»rf
_
Alllms
jDynamcs o -
; SyiieriLbrjry
,
Lkie Iradton
Solver Componert
tXlOTrl LMs
f
OUgnogucs GMmIc
ntegrating fra» 4 96 to 4 97 1
Integrating fro« 4 97 to 4 98
Integrating fro« 4 99 to 4 99
Integrating fro* 4 99 to S
Step 361 Time- 5 0422e+0O0. step size- S 0000e-002. step tactor* 1 S000»*000 accepted
Run complete at 16 45 29
il I
FIGURE 5.39
In the subsequent step, run the program with Initialization run mode. As it is
inished, go back to Dynamic mode. Then, open the plot sheets for both the controllers.
f
The regulatory behaviour is illustrated in Figure 5.40 giving changes in feed temperature
DYNAMICS AND CONTROL OF FLOW-DRIVEN PROCESSES 255
are not
in the feed data sheet. For brevity, the faceplate and configure dialog box
included in the Aspen Dynamics window, shown in Figure 5.40.
2 a ts 7 I Dynamic
3
-
r u loo? s
Ot* :iC£3=
-
-
trx*'*-** C fi
U
| _
Ml
in
53 a
T mpeiatuie loop
-
i
; So
1
3
0 05 1 1= 2 25 3 35 4 45 5
Time Hems 000e+000. accepted
FIGURE 5.40
It is obvious that the reactor liquid level remains unchanged with a change in feed
temperature since there is no interaction involved On the other hand, the reactor .
he column has total 25 theoretical stages (including condenser and reboiler) and
operates at a pressure in the reflux drum of 18 psia and reboiler of 21 psia. The
ow rate is 285 Ibmol/hr and reflux ratio is 2 2 (mole basis) .
.
256 PROCESS SIMULATION AND CONTROL USING ASPEN
Feed Specifications
Component Mole %
BOTTOM
benzene 45
toluene 55
In Table 5.1, the reflux drum and the base of the column (the 'sump' in Aspen
terminology) are specified. It is fair to use an aspect ratio (length to diameter ratio) of
2 (Luyben, 2004).
TABLE 5.1
Item Vessel type Head type Height / Length (ft) Diameter (ft)
Sump -
elliptical 5 25.
The column diameter is 5 ft. Use default values for other tray hydraulic parameters
(e.g., tray spacing, weir height and weir length to column diameter ratio). Consider log-
mean temperature difference (LMTD) assumptions for the total condenser. Actually the ,
LMTD is calculated using the temperatures of process luid and coolant In the simulation. .
f
assume constant reboiler heat duty and apply the UNIFAC base property method.
(a)Simulate the column model to obtain the products mole fractions.
(b)Keeping the default level and pressure control algorithms unaltered, inspect
the servo as well as regulatory performance of a proportional integral (PI1
controller that is required to insert to control the benzene composition in the
distillate by manipulating the reflux low rate.
f
Simulation approach
(a) Select Aspen Plus User Interface and when the Aspen Plus window pops up.
choose Template and press OK. In the subsequent step, select General with
English Units and hit OK button. To open the process lowsheet window, click
f
Creating flowsheet
From the Model Library toolbar, select the Columns tab. Place the RadFrac model on
the lowsheet window and add the feed as well as two product streams. Renaming all
f
.
imaimAJmm lAdftiauugj*) ji.ojjiIzJII
o: 6 -i
'
twin TmT
FIGURE 5,42
Configuring settings
Ah we hit Nt'xf followed by OK button, the ROtUp input Corni appcarH (h«m* Kitfun! .43).
The diHtillation problem Is tilled oh: 'Cl0S6d-l00p Performanct; of a DihI illation ('olumn'
Iniportuntly, une tin- 'Dynanne' input noxle
mltiMI
lit I:: | n
NiMH
I
(
(
FIGURE 6 43
258 PROCESS SIMULATION AND CONTKOI, IISINC ASPEN
Figure 5.44 includes the Aspen I'lus (iccon/ilin infornuition We can fill up the
.
'
-
i '1
t t vW o**
f
mm jj ai LiinsMibd ad j ..ijj-j -j ai i i
O Ma
o - J-.
f> 1 " .
!(. [Ilflll,
I - MCM
fToirTiTTT
I
_
J M.>i*t 4» M 1><»
Uvain . I -l|a. n/ ' Hi'W.O | ' } rhjiHg fi tVio«imVcrt [[ Aapcn fits »?t>8 111'
FIGURE 5.44
We like to see the composition of all incoming and outgoing streams in mole raction
f
basis in the inal results table. Accordingly, we use Stream sheet under the Report
f
i MHl iJM«l _
BMSLSfcalgl "'I "'I I I "I "I ~l g
""
JJ I I I FV . _
LLM_Lia|
"
FlipdrtOplta 351
bnnil | rbwihWI | BM ./SIiiimi| Piorniy [ ADA |
0 S*-*. ",
'
, w. «..
.
O -
I! . P W .Hum.
. J 1
FIGURE 5.45
DYNAMICS \NI) CONTUOI, OK I'l.OW DRIVKN ('UOCICSSKS 259
Specifying components
'
From (In- Data Brow.scr, hoIoc! ( oniponcnta/Speeipcctions to open the componont input
'
lorm In ill'' lahlc. shown in l- i mc S 1 (> Ihc Ivvo species are dclincd .
Qp
him nl
ni.urm rrniioii mo
Itarimi Ophor*
I Ml/ tC Qrm)«
ml
Mi
FIGURE 5.46
in Aspen simulation a property method originall} Includes several models for calculating
,
the physical properties For the distillation example, set the UNIPAC base method b>
scrolling down (see Figure 5.47).
Specifying stream Information
,"
next, (.pen Streama IFEED IInput ISpecifications sheel Entering the given
Values lor all State variables and teed eompo ion Ihe slream mpnl lorm looks like
.
Figure 6 'IH
Configuration sheet and then insorl the required datn (see Figure 5.49)
260 . PROCESS SIMULATION AND CONTROL USING ASPEN
f
JjiJ J J
_
LJ~- I I I Jje
HU I J±hJill -ZliiJuLJgij
Relerenced |
0 Slream Oae* Piowxtyrwthocbt moded
B | S<i)*ww
PtocMt IVP*
Picwtymethod juNiFftr Tj
Q Cutlom OMi B«e meitwd jufJIFAC r Modlv P<ope<iy rmdett
O Repwt Opicna Heniy corrvcuwitt [
'
Fiw-watm method 3
~
P»e<rfoc«npOfi«rt9
3
Etecddyte c ciisbon option:
; O MrCorpi Chamtiiy ID r
17 . : .
_
Potjma
B l Property Method!
f
.
_
j EHi Un
9 CD Wotecd» Slftjcti«
a
LJNIFAC wlh AedKrh-Kwof.g elation d tlate and Henty't law
Qt | Prap-S«ts
£ r j Adveoced
hpot Complete
[it* MbMi/Splilleri j Sepstfon | HhI Exchangees Cohaint { Readou | Piettue Changeft | Manptiatori | Sokh | UteiModrt |
STREAMS DSTVAJ Di'.H FladFiac E act Mijfrac SCf<ac FelioFisc fWeFrac BaichFrac
C \ d Fitfcn'Asien dd IM HUM
FIGURE 5.47
-
fi- ?d( [Ula T04I1 Run not tibrary Wride)v< HH'
$ JJ Peiro Owaclefuabon
Lghl-End Prcpofto
3 |M*f 3
I d
'
Pteudocorvowta 225
CotrpGrem Vofcje
SEN2ENE 015
HonyCcrrps d TOLUENE 055
UNIFAC Gtohm
I p.. d
) Comp-Groupi
> Compel
d
_
B _
) Potywi lanlltow [Mote
-
,
i Data
_
J MM*
-
Ji Pr -Sett
$ MvarvceC
.
_
| BOTTOM !oUI F
-
Jl) FEED
o
EOVMMu
STREAMS 1 DSTVA)
FoH* PMF1
Book
FIGURE 5.48
DYNAMICS AND CONTROL OK KLOW-DKiVEN PROCESSES 261
cz:
'-
a
Mr
a
3
9
-
ap>
1->- d J
FIGURE 5.49
In the subsequent step (see Figure 5.50), fill up Streams sheet with informing t\
location 113th tray (above stage)!.
1 r - i t
- - : -
nr
j Kim
w
-
I .» -
, »«
fTTTTB
FIGURE 5.50
262 PROCESS SIMULATION AND CONTROL USINKi ASPEN
tttt
O hpu 1 e
iesuts
f
25 21
*
_
j tOP *
(j Bocks
8 S«tup
Q Design Specs
D Vaiy
S Heaierj Co era
lJ Putnoarxjxis
D*c»Tlefs
o BfoencJes
® Reeciavs
D Condenser Hc n.'
Reboter Hcijr\-es
_
J Tray Sung
.
D Tray ReU-g
CJ Peck Sa ig
Peck Relng
0 Praoerties
«3 -
1L 1 w
:6 Stages lf*Jtrcr*s*j-.
MKm/SpUleis Sep«alots j Heat Exchanges Cduran* } Reactors Piessute ChangH! j Man atojs j SofcU ) Uiei Models j
HZH1
STREAMS
(C-I-C-S-M-jfr-C-cT-
' DSTWU SCftac PeltoFrac R eFiac
"
SatchFfdc
'
FIGURE 5.51
set constant. Here the temperature approach represents the difference between the
process temperature and the coolant outlet temperature at the initial steady state
Note that among the heat transfer specifications the coolant inlet temperature and ,
temperature approach may vary during a dynamic simulation whereas the specific ,
heat capacity of the coolant is ixed during a dynamic run (see Figure 5.52).
f
For the reboiler , it is fair to use constant heat duty computed in the Aspen Plus
simulation. However , the reboiler duty may be changed at dynamic state either by
manually or automatically with employing a controller (see Figure 5.53).
Entering geometry data for reflux drum and sump
The reflux drum and sump are specified in Figures 5.54(a) and (b) with their given
geometry data. The information on vessel orientation, head type, length (or height)
and diameter are used to compute the vessel holdup .
DYNAMICS AND CONTROL OF FLOW DmVKN PROCKSSP S -
263
j|a;|B| I I v\ r;|Tlal».|<|w| n i
i r -i i i rv i .isi J2J
] ±LI iijri 3 uu Nj
Sim
p
"
I «np«amB appiMch
9 *«mCoolen
CJ (
_
]
« : M
_
j
a Reboier Moves
ii
_
j
Pat* Siing
i faw Rang
i
_
o
o
o Reoort
o ijm- SutraArm
0 O iarrjc Eaueme
Bock Opuon*
J 1
Mwr»/SpfcJlM« j Sep«arori | HMTExchangen | C(*in>> j RMctai | FtenusChangen | M. ScWi I UwHcmMi J
-
CM
'
MUM
5TRt»S FSfil SSpB
Book
FIGURE 5.52
1 nr 1 leal i ibi
"
-
'
I tmwwj-di
f
S Cleans
'
_
) x<JerMr Mar*
-
O 'JMr SoOTOJUTW
0 a«*C«bor»
1 *L
FIGURE 5.53
264 4- PROCESS SIMULATION AND CONTROL USING ASPEN
In jJIZJ
l 1 1 mm
3 ±l_JiiJR 321] nil
"
3 Efb |m
2 PrccMei
, /Condense, j/Reboiier /Reltmi Dmml Sump j
.
Jfl Steams ~
l Vaiy
i Lengtlk 5 fl z
9
_
Heater CocJeo
?ii«>anxx>ds : Diwnetei 25 t z
CJ
f
.
i Oacartefi
5
_
Iniialspdcificabon
s Reactone
Cordansaf H v
jj RebalerHcu-wss
_
i Tra)' Stw>3
a TrajRing
EJ Pack Sons
_
j PactRalrtg
9 ?ro(j«rt:«
0 Eawnatei
9 Conveyance
9 Repor.
9 Uiw SiijfDU;nes
w Oynanac
'
iTangen' lc- tange'-tl ertgUi J Heod Type-EBclic oi Hemisphenc-al Acfuattenglh it Head Tvp*=FU
t
25 5t»5£S t; Rmpifotnd f
Mwen/SpliHeft | Sepaisios j H»at£xchangeit j Cokimne | Rescloit { Piessue CKsngef: | ManipUators I Solicft j UserModeis {
Mattel
'
Mom
STREAMS J Mbw FSpil
FSpB SSpW
For Help. pre«i FT
1 $ Aspen Oyna j
_
. psty Tta . j j-j spxxK hg . j Calajbto; <m
FIGURE 5.54(a)
S3ESSESS
s % EdR Mew Oa i Too* Pin lat Ubwy
f
3 -I I'M lEi J Si
i r jjj_fv
_ .
"
V1
_
oc-rt es 3
J/j areams
ioda
ft
-
RMFRK : Headiippe
Totatiqi*dvoU?ehac!wn [
~.
FIGURE 5.54(b)
DYNAMICS AND CONTROL OF FLOW-DRIVEN PROCESSES 265
"
f : fV -
I - imi I lai 1 1
3 iilFi 3 ill d l \
~
3
.
jj Stan
-
Mi
_
j «y .
d
« :4 0 726££ 90
Oecanfm »
_
j Retokr Haxves
lj Tufting
:
_
l i
0 BtockOptofa
- . - - V
|rt>J
SeG«*tn | Hea EKChaqot | Cckm | ReKUfl | PieumChanpn | M««i>lois { SotcH | day Mod* |
-
D->
Mam*
1 .1 pr-u j.p t.
_
Caiaiacr
FIGURE 5.55
_ _
window .
266 PROCESS SIMULATION AND CONTROL USING ASPEN
Setup
1 1
_
P pel s
31 9mamt -
3 F-m*;
Bseki
Sue-am Table I
Rnctane
z\ d
-
Com Optons
Zl
-
ConvOder 4cfcFlowlhcnoi/t.
-
I Sequence BENZENE 1 036 270.000
.
268,964
j RoA-shee!r>; Options
TOLUENE
_
Jj EO Cc guratan oleFtac
9 Pijn Status
TOLUENE 0 S97
.
0 550
.
0 056
Convef ence
Roula Avoleble
"
plj Uureit/Spllten ] Separatort ) Heal Ewhangett ] CoUww ) Reacton | Piewuie OwrgBrs | MampUatois | Srtdt ) U«e( Models |
Man
STREAMS '1 Mi»«. M « FSpa
r rJJ SSl*
rA
|4)Ch8i»5 W..,|| Aipen Phm... apm.D) .,- | apen.IVia-.. } apocfeheat:. [ CefcuMof | « Q P'tVi 11 53
f
FIGURE 5.56
_ _
window appears (see Figure 5.57) accompanying with the closed-loop process low
f
diagram. The flowsheet actually includes the three default control schemes LCI, PC2
and LC3 to monitor the reflux drum liquid level, top stage pressure and column base
liquid level, respectively.
In the present discussion, we do not want to change anything of the three
automatically inserted control strategies. All data, including timing parameters, ranges,
bias values and controller actions, remain untouched. A little detail of these control
structures is given below.
Loop 1
Controller: LCI
Type of controUer: P-only (since integral time is very large (60000 minutes))
Controlled variable: liquid level in the reflux drum
Manipulated variable: distillate low rate
f
»g [Dynnmlt J h« .IS'ByiUEiafl
f
inM
-
t
LCI
r>0 .
C<rtenta d Snmiaaon
£5 9 »o--
BOTTOM fO
LlJ
IBjnsilcal 0.(10
'
Aspen fVjs Mty8 [>i_
I A »n IVvrtc-
FIGURE 5.57
Loop 2
Controller: PC2
Type of Controller: PI
Controlled variable: top stage pressure
Manipulated variable: coolant inflow rate
Controller action: direct
Loop 3{
Controller: LC3
Type of controller: P-only
Controlled variable: liquid level in the column base
Manipulated variable: bottoms low rate
f
r tt o p ,
ar s _
Bp
J
FIGURE 5.58
Again hit expand button next to the ControlModels icon Then select PID controller .
drag it to the low diagram and drop the control block near to the top product stream .
f
Uae l*.
lO i£i
Si*?0-
no
-
TwT-o
1
.
Id V
§1 - t
-Q: h tgnm K>
9
-
iai»i
PS Rnuu 3
FIGURE 5.59
wiih holding the ControlSignal icon over a port, the name of that P
.
*ame *
on the
Anyway, move the pointer and release the mouse buttonselect 0fe g wo
the dastillate compos.tio
OutputSignal originated rom TOP (stream) block. To
f
CONTROL OF FLOW
DRTV rPPn cc 269
of benzene as controUed variable, choose 'STREAMSCTOP") ZnCBENyFTWv
raction' by description in the Select the Control Variable 4 bo Fi 5 6of
f
W lOyti.Mc
siREMnrniPi fmanaa
SIBEAMSI'TOP'l fiTtrtTOi U
si a
.
.
SIRDMSTTOPtP
STREaMSriOPlZmrtBEftZE Moslhacten
STSfAMSlTDP
"
'
IQUJE
Ucli oal
t
t
FIGURE 5.60
As we press 0/!l button, the cursor becomes a solid black arrow representing the
input signal to the controller. To transmit this signal to the CCT block, connect the
black arrow with a port marked InputSignal. Since this signal conveys the process
variable (PV) information to the CCT controller select 'CCT.PV by name with 'Process ,
variable
'
by description in the Select the Control Variable dialog box (see Figure 5.61).
-
Ulxl
0«B'<BQi 1% X? |Dyn.mlc
.
ri
.
FIGURE 5.61
270 PROCESS SIMULATION AND CONTROL USING ASPEN
Hit OK button and obtain Figure 5.62. Obviously the CCT controller is partially ,
configured. To complete the top composition loop, the controller output should be
connected with the manipulated variable to pass on the signal .
s FtwrilMl Window
r Utt G..d |0 05 t!
) Modd)
( Parametef Types
>C2
:
| Pol TflJW
i Pre5»ureCnsr»geRi 0 LCI
J PMH . rc- Bl
| Prcceduw
) F actora
| Seoaratora
I Slreani Type* CH FEED [
J Cwviedion LC3
.
>o
CortolStgnal >o- o r
r
ConteniiQl ContiptSignal
-
| B0TT0M]</
aart| p , m " Bopk | -jjChap-er 5 Mi | Jrj UaeMTgms | Adobe Acroba | Goo e Wnd | A*>en Ptu: C Dm
FIGURE 5.62
the blue outgoing arrow marked OutputSignal from the CCT block. As Select the
Control Variable dialog box appears (see Figure 5.63), choose 'CCT.OP' by name and
press OK.
Immediately, a solid black arrow representing the controller output signal is
automatically generated. Move the mouse pointer to reflux stream and make a
connection to InputSignal2 port. To use the reflux low rate as control variable, select
f
'
BLOCKSC'RADFRAC). Reflux.FmR' in the dialog box and click OK (see Figure 5.64).
Now the binary distillation column is coupled with four control schemes, LCI,
PC2, LC3 and CCT, and the closed-loop process looks like Figure 5.65. The subsequent
discussion includes the modification of different tuning properties of the CCT
controller.
DYNAMICS AND CONTROL OK FLOW DRIVKn PROCESSES -
271
HKBB
OwtB W | Dynamic J
M l
r UTtlMln *J Li \i r %
BIBB
Piocmi rmbto
3
Edit tnq SimuIsIton
L Uj . i
FIGURE 5.63
0 * H .[! (B V |Ovn.mlr
IS, , . 'i . 1.. - V..- , * - iDl
_
3 r it 7f o-d|uir, i; wit » s
yi||
i
M oed tW*
t
FIGURE 5.64
272 PROCESS SIMULATION AND CONTROL USING ASPEN
6* a »Q l W |Uyn.mlt J >
f
r « Orflofi ,
f
,
>o
J Stmam Tipca
OH FEED h
>o
J >0-- v m
Editing 5i»ulstior.
FIGURE 5.65
then hit Configure symbol in the faceplate to open the Configure dialog box (see Figure 5.66).
3
* |r n -it c..i|o"-, jj « n
r LCI
CwT-o
Tuwnp | Raryii jf ffamg J Ottw
g>-ffnoT-
OteianMnMlai
y
TtmepaianeieK : H BOTTOM K)
Bm |o
.* V
irOjaJim |» 1- d
1- d
4 .-
-
»
FIGURE 5.66
DYNAMICS AND CONTROL OF FLOW-DRIVEN PROCESSES 273
Obviously, some of the default values set by Aspen Dynamics are not acceptable.
For example, the operator set point value of process variable (benzene composition in
distillate) should not be greater than L Secondly, the CCT controller action must be
Reverse'. In addition, the value of control variable (reflux flow rate) at steady state is
'
process variable, control variable, bias and ranges. Note that the controller action is
changed only manually. It is wise to initialize the values by the help of Aspen Dynamics .
For this, press Initialize Values button in the Configure dialog box and use 'Reverse'
controller action. It is obvious in the window, shown in Figure 5.67, that the values of
SP, PV and OP in the faceplate change automatically to their steady state values. If
this approach fails to initialize the simulation of controller model with the steady state
data, check and replace, if necessary, the values of PV and OP with their steady state
values by double clicking on signal transmission lines (input to the controller and output
from the controller).
Alltto 3
r Gndi-JCB zl 2* H
r
acci Mill
.
O LCI
10 3437
10 9437
|49472 na
< -{TfLCL.j
LCI
if
Turong ] Ranges ] nam T >o-
1
1 BOTTOM K>
0|wdlo«»t|)«* 10 343 33 IfcmoMbmol
Tuwigpawieleit
|o
Conbofct Kljon
r Dieo
nables 1494
» loot DynmcaOOO
J
It | JlJ Ouplg 5 >4ct»<I W j AJOOe oif friM \ tm>on t» *»l»-0« jj '' «
t
FIGURE 5.67
Db Tat* Kr XVhoa* t*
r
- W -; . (ft |lnlllBll?nll<in
jo g*-'.' >C
BBSiwi
Q Hc
>o-
| H bohom f<'
P Opts range
0><pki
Ready
5.11 r
'
FIGURE 5.68
However, here we use the typical variable ranges, shown in Figure 5.69 .
»ff |rnitializMion « i S1 © 0 [3
f
3
3
3 SrUMon
% iaizi
>o -
L
-Q H 60MOM [<>
nwpwjifiiiii [T
.
__g groups 664 lo 670 of type Explicit Linear
Solwing group c'l sire S type Linear
f
fun coaplete at 11 16 14
J
..
/is
FIGURE 5.69
DYNAMICS AND CONTROL OF FLOW DRIVEN PROCESSES
-
275
Process variable
Range minimum: 0.85
Range maximum: 10
Output
Range minimum: 10000 Ib/hr
Range maximum: 120000 Ib/hr
It is important to mention that it is a good idea to carry out Initialization as well
as Dynamic run after each new change in the control scheme so that any error in
controller installation can be detected individually.
We have now completed all required control specifications for the top composition
loop In the ongoing study, we prefer to conduct the simulation experiment to observe
the designed controller performance continuously for 5 hours. As done for the previous
CSTR problem, similarly either simply press Ctrl+F5 on the keyboard or select Pause
At rom the Run menu and put 5 hours as Pause at time.
f
In the next, we will inspect the CCT controller performance first dealing with the
servo problem followed by the regulatory problem.
appears. Then open the Configure as well as ResultsPlot dialog box The second one is
.
basically a blank graph sheet that presents the variations of process variable set point
,
and faceplate are correct. In the next hit Run button to start the dynamic simulation. The
,
plots, shown in Figure 5 70, illustrate the servo behaviour of the PI control algorithm with
a step increase (0 9437 - 0 97 at time = 1.51 hours) followed by a step decrease
.
(0 97 - 0 9 at time = 3 hours) in the set point value of the distillate composition of benzene.
.
gain of 10 %/% and integral time of 10 minutes. These values have been chosen based
on a pulse input test in the distillate composition of benzene and using the trial-and-
error approach It should be kept in mind that the objective at this point is not to come
up with the best control structure or the optimum controller tuning. We only need a
control scheme and tunings that provide a reasonably good tracking performance to
drive the simulation to a new steady state.
Remember that to edit the plots, shown in Figure 5 70, double-click on different
elements of the plots and modify them as we like.
Viewing regulatory performance of CCT
In order to investigate the regulatory performance of the CCT controller, we give a step
input change in the feed pressure (21 -» 23 psia) at time = 1.48 hours and that in the feed
temperature (225 -» 230oF) at time = 3 hours. The PI controller tuning set provides good
disturbance rejection performance (see Figure 5.71) although the tuning parameter values
.
gain and integral time) have been chosen based on a pulse set point input change.
276 PROCESS SIMULATION AND CONTROL USfNG AST'EN
r a it Mi-nr 2} £ i s
»o
jtOTTOMfO
FIGURE 5.70
3 S H « (i V? lOyn.mk
t Eyewntcrari
Tyl-o
--
>o
HQ- -. ccr
'
33 31
'
il
3
Z«-B£KZE«-J
5 05 1 If 2 2£ 3 3S < «5 5
Time Hoiks
1
j
FIGURE 5.71
I.:IM
.
:
1
s
.
41 . .
! m
.
FIGURE 5.72
"
r :: f - >. »- »
1 55 i; 7T i: '».»«
'w. i .
-
f v ..
c MM
-
FIGURE 5.73
278 . PROCESS SIMULATION AND CONTROL USING ASPEN
EC MS*}
-
9 -
=inj l
"
_
r Tt (0 05 3 i£ i « W 1* 3
LCI
>o-
%j am {33-0
SP| 100025
PV; I00O25 LCJ >a-j
--
BOTTOM K>
Tunuvj | Ranges | Ffenr ] Olha
11
Cortiolw action -
Drad > £
i
r Reverte S?
1
-
If r
0 0 5 1 1 5 2 2 5 3 3 = 4 45
,
5 i 2 2.5 3 13
Time Hours Time Hours
TEE
eed
t
FIGURE 5.74
Clearly, the CCB controller shows satisfactory set point tracking performance
against a pulse input change It is observed rom Figure 5.74 that owing to strong
.
f
interaction between the two composition loops of the distillation column the set point ,
changes in bottom loop affect the top product composition Similarly, when any set .
point change is introduced in the top composition loop, the bottom product composition
will also be affected.
DYNAMICS IND CONTROL OF FLOW-DRIN KN I'HOrKSSKS 279
PROBLEMS|
5 1 A feed mixture of benzene and toluene is fed to a flash drum (Flash2). The
.
separator operates at 1.2 atm and 100oC For dynamic simulation, required feed
.
Feed
Temperature = 25°C
Pressure = 3 bar
FLASH
Flow rate = 100 kmol/hr
>o
PI
Component Mole fraction
benzene 06
.
toluene 04
.
'
a) Use the SYSOP0 property method to compute the amounts of liquid and
vapour products and their compositions .
.
b) As shown in Figure 5 75, employ a PI control scheme to monitor the
.
(c) Show the closed-loop servo performance with +10% and then -10% step
changes in the flash temperature .
5 2 A vapour mixture of toluene, methane and hydrogen is heated using a shell and
.
tube heat exchanger (HeatX) The superheated steam is used as a heating medium.
.
Temperature = 2780F I
HOT-OUT
Pressure = 500 psia
Cold Stream Out
Component Flow rate
(kmol/hr) ICOLD-INf j cold-out hoi Temperature = HOOT
Pressure = 498 psia
toluene 200
Dead time pi
methane 2300
HOT-INK >AT
hydrogen 1000 >o
ii
Hot Stream In
Temperature = 1160oF
Pressure = 14.7 psia
Flow rate = 5110 kmol/hr
(a) Simulate the heat exchanger model using the shortcut method, counter-
current low direction and NRTL-RK property method.
f
(b) Include a PI control structure to observe the set point (cold stream outlet
temperature) tracking performance and the manipulated input (steam inflow
rate) profile. In the closed-loop simulation experiment, assume that the
temperature sensor takes 1 minute time (dead time) to measure the controlled
variable. Report the used tuning properties.
(c) Examine the regulatory performance by introducing + 10% and subsequently
-
in Figure 5.77. The low rates for both the disturbance or wild stream (Fw) and
f
the manipulated stream (FE) are measured, and the measured ratio, R = FE/Fw m ,
is calculated. The output of the ratio element is sent to a ratio controller (PI) that
compares the calculated ratio Rm to the desired ratio Rd (set point) and adjusts
the manipulated low rate accordingly.
f
Input 2
Ratio
>o
PI
> | POT >
Input 1
Hint: Double-click on Input 1 transmission line and ill up Tables 5.3(a) and (b).
f
TABLE 5.3(a)
Value Spec
*,
>STREAMS("E ) Fcn("ETHANOL")
.
100.0 Free
Value Spec
w
>STREAMS("W ) Fcn("WATER").
150.0 Free
Use Initialize Values button and incorporate the given tuning properties before
running the program.
282 PROCESS SIMULATION AND CONTROL USING ASI'RN""
5 , 5 A reboiled stripper is used to remove mainly propane and lighter species from a
feed stream, shown in Figure 5.78. It has total 6 stages (including condenser
and reboiler) and no condenser.
The bottoms rate is 100 Ibmol/hr and the column top stage pressure is 150 psia
with a column pressure drop of 8 psi. The diameter of the stripper (Stages 1 to
5) is 6.5 ft. The reboiler heat duty is assumed constant, although it changes at
dynamic state. The sump has elliptical head with a height of 5 ft and diameter
of 2.5 ft.
For the closed-loop simulation, use the following data:
Dead time = 2 minutes
Magnitude of noise (standard deviation) = 0.01 Ibmol/lbmol
Proportional gain of PI = 1 %/%
Integral time of PI = 20 minutes
Controller action = Reverse
PCI
Feed
Temperature = 40oF
Pressure = 160 psia
Feed stage = 1 (above stage)
Dead time Noise Pi
Component Flow rate
(Ibmol/hr) >o
c, 60
c2 75
C3 150
n-C4 175
n-C 60
5
n-C 35
8
composition control scheme (PI) coupling with a 'Dead time' and 'Noise _
closed-loop data and execute the dynamic simulations to test the developed
composition controller performance .
5 .
6 Ethylene is produced by cracking of ethane in a stoichiometric reactor. The
irreversible elementary vapour-phase reaction is given as .
C2H6 -i C2H4 + H2
ethane ethylene hydrogen
Pure ethane feed enters the reactor model shown in Figure 5.79, with a low
,
f
rate of 750 kmol/hr at 800oC and 5 5 atm. The reactor operates at inlet
.
pi
> M >o
f
U) Using the SYSOPO thermodynamic method, simulate the reactor model.
(b) Develop a control loop as configured in the low diagram to maintain the
f
desired reactor temperature by the adjustment of heat duty. Considering the
measurement lag of 1 minute, inspect the servo as well as regulatory control
performance. Report the tuning properties used to achieve a satisfactory
closed-loop performance.
5 . 7 A binary feed mixture consisting of methylcyclohexane fMCH) and toluene is
introduced above tray number 14 of a RadFrac distillation model, shown in
Figure 5.80.
O 1 phenol [
O 1 FEED h
tray number 7 of the column The two input streams have the following
.
TABLE 5.5
FEED 105 14
.
181.44 kmol/hr 0 5/0.5
.
(MCH/toluene)
The column has 22 theoretical stages (including condenser and reboiler) with a
total condenser The distillate rate and reflux ratio are given as 90.72 kmol/hr
.
TABLE 5.6
Sump -
elliptical 15
. 0 75
.
The column diameter and tray spacing are given as 2 m and 0.6 m respectively,
,
(a) Simulate the distillation column using the UNIFAC property method to
compute the composition of MCH in the distillate and that of phenol in the
bottom product.
(b) In addition to the default level and pressure controllers, insert a PID structure
to control the MCH composition in the top product by manipulating the flow
rate of PHENOL stream.
(c) Produce the plots to show the closed-loop control responses, and report the
tuning parameters, control actions and operating ranges for controlled as
well as manipulated variables used.
REFERENCES|
Luyben, W.L., (2004), "Use of Dynamic Simulation to Converge Complex Process
Flowsheets", Chemical Engineering Education, pp. 142-149.
Seborg, D.E., T.F. Edgar and D.A. Mellichamp, (2003), Process Dynamics and Control,
2nd ed., John Wiley & Sons, Inc.
CHAPTER
6
Dynamics and Control of
Pressure-driven Processes
61 .
INTRODUCTION
are inserted where needed, to provide the required pressure drop for material low.
,
f
Control valves are installed where needed, and their pressure drops selected. For good
,
control , the pressure drop across a control valve should be greater than 0.1 bar. The
luid that lows through a valve should normally be liquid-only or vapour-only because
f
It should be pointed out that for a pressure-driven case, we must not insert a valve
in the suction of a pump or at the discharge of a compressor (compressor speed or its
equivalent compressor work is manipulated). The control valves are positioned on the
luid streams such that the controllers can manipulate the valve positions.
f
The simple low-driven dynamic simulations have been discussed in detail in the
f
285
286 PROCESS SIMULATION AND CONTROL USING ASPEN
Problem statement
The methyl tertiary butyl ether (MTBE) column configuration (Jacobs and Krishna,
1993) chosen for the simulation is shown in Figure 6.1.
Temperature = 320 K
Pressure = 1 aim
Flow rate = 711.30 kmol/hr
Jy tCl-fpisT-Q
Feed stage = 10 (above-stage) 0| METHANOL Ft
CV2
PUMP
/so-butene (IB) 0 36
.
n -butene (NB) 0 64
.
'
9244(M
-
.
IB
Forward rate: rf= 3 67 x 1012 exp
.
RT xMeOH ,
-
134454> VMTB!'
Backward rate: r,, = 2 67 . x 1017 exp
RT
DYNAMICS AND CONTROL OF PRESSURE-DRIVEN PROCESSES 287
f
640.8 kmol/hr. In the MTBE synthesis process, it is desirable to obtain a bottom product
containing high-purity MTBE and a distillate containing high-purity NB. In Table 6.1
the reflux drum and the sump (the next-to-last stage in the column) are specified.
TABLE 6.1
The pump delivers the liquid stream POUT at 11.7 atm. The compressor (isentropic)
has discharged the vapour feed FV at 11.5 atm. The three control valves (adiabatic
lash) CV1, CV2 and CV3 have the outlet pressures of 11.5 atm, 10.8 atm and 11.3 atm
f
(b) develop the control configurations to achieve the desired product purity under
disturbance input.
Simulation approach
(a) Start the Aspen program by double-clicking the Aspen Plus User Interface icon
on the desktop. Then select Template option and press Oif (see Figure 6.2).
: aM
FIGURE 6.2
Copyrighted material
1
288 PROCESS SIMULATION AND CONTROL USING ASPEN"
We choose General with Metric Units option and hit OK button (see Figure 6.3).
016*11 _
L I _!_) *d U-i-lfcl I M 3 I I I J jU J J
l l I I l I HI J J J
_ _ ,
J U ali
mi in -
I
.
>.. I
0*. I ,.J h
1
r -
P.*
FIGURE 6.3
When the Connect to Engine window appears, use the default Server type (Local PC).
Creating flowsheet
The process low diagram includes a feed pump a feed compressor, a distillation column ,
f
and three control valves. The complete process lowsheet drawn in an Aspen window
f
should somewhat resemble the one shown in Figure 6 4 Recall that Aspen has a tool . .
in the toolbar that automatically takes the user through the required data input in a
stepwise fashion. The blue Next button does this.
-
o-
" if - *x
FIGURE 6.4
DYNAMICS AND CONTROL OF PRKSSIIRE-DRI\T N PROCESSES ,
289
Configuring settings
At the beginning of data entry, ill up Global sheet followed by Accounting sheet under
f
Specifications of Setup folder. Moreover, select 'Mole' fraction along with 'Std.liq.volume'
low basis in Stream sheet under Report Options [see Figures 6.5(a) ,
(b) and (c)].
f
41
& -
Sim
Pi bw
HUB .
|
MMH .
fMM,
r 3 EC
ir S
FIGURE 6.5(a)
:.
>-B| *ifi »] atmnn-at .! 3 I 1"! jiJ 3
-
T i-l h r* .1 -IE; : II
- a*-=- = -
aa-2-
FIGURE 6.5(b)
290 PROCESS SIMULATION AND CONTROL USING ASPEN
j r
_ _
I I -I 1 ,
;
| .imi I M MgJ
nafhi >bJiiiF-a Jj
0
P lr<Ul«il>«tmd«.CTvtlci.i
f
STREAMS Pmp Conif MCow fV»
FaH«te . (wun gFoUeiiVwenButllV HUM FUntiwIlrpj r cWc
FIGURE 6.5(c)
Specifying components
The components involved in the example system are MeOH (CH40) IB (C4H8-5), NB ,
(C4H8-1) and MTBE (C5H120-D2). Within the parentheses, the chemical formulas
used in Aspen terminology are mentioned (see Figxire 6.6).
I r M I IT I W! I [a) -«l»«l
3MS
h EQWmfMi
NB i-BUIEHE C4HB'1
MIBE
.
_ J PwoCh.
.
J Con
. - . PA
J Si 35! I u'*
i »tr»iiiiin
FIGURE 6.6
DYNAMICS AND CONTROL OF PRESSURE-DRIVEN PROCESSES 291
1 d 1 J
1 3
'1 J 1 J
1 d 1 3
I J
1 J
1 r
FIGURE 6.7
a r i j-j-i ..if; 3 ») qU I
9 S3HS 1I
wo
I
-
ass
.
.
1
FIGURE 6.8(a)
292 PROCESS SIMULATION AND CONTROL USING ASPEN
i r m i nr i ei i m mi
-
DgrfcR FH »| |..|[Tr
) mum* i i
HPS 3
la !.
J r: -.
"lur
1
il Pi9
mibe
i
'
SJIENIS
1 E0v«
.
_
j
j Fl
_
*
li FV
-
ICIHMBI
Res*;
FIGURE 6.8(b)
_
i r- I -1 I- IT
_
I lal I la
afbRTr H JulF 321) QLJ!!!!
-
fir-.
q hm 3
J tUIENES
1
S PI
J POUI
_
COMffltSi
Q Salup
f COVmUk
i tO Ira'
f
. ; MM
FIGURE 6.9(a)
DYNAMICS AND CONTROL OF PRESSURE-DRIVEN PROCESSES
ir LUJV _
'I'N Mai gH
0
-
a, n
A. F\
-
ij M(t>«nOL
.
ai n
.
ii w
0
a
a
a
a Bo.* O(*ora
3 MmnMtatm p
s t'J in"
Spk BMW
_
J
Pm
w
Iflpul
Btae« Ophorj
tP .»rtNei
9 IQlfVJ
a snear O-.ji
FIGURE 6.9(b)
»i Dial n.i
| - 'MM
ni»pMii|tH f -3 * I-
r c w**, I 1 1 r I
[23 ECVwwim
: : .
1107
FIGURE 6.9(c)
294 PROCESS SIMl'UVnON AND CONTROL USING ASPEN
d (BfacK-VlrVJ.rllnpui 0
* E* C«* Tadb Rrf, Plo i±1Mr wnx*. h*
"
© .:>-:
_
J
2] Kea*:
H EO-M.,
- ,
£jc* ?Ow»ii
_
.HIM
FIGURE 6.9(d)
l r _ _
L _
LLJV 'I -isi
a -'
i
.
-r - y
m r
RiMP
f :
O Sim v _J
£TRDMt
J
FIGURE 6.10
DYNAMICS AND CONTROL OF PRESSURE-DRIVEN PROCESSES 295
In the list on the left, choose Blocks I RDCOLUMN I Setup to ill up Configuration
f
sheet (see Figure 6.11).
'
Wal _LJ «m *) fJMfcl'.KM -I -il - I M -I 71 -I _J
i T I I I'
-
I - Wl I Ml -
3 3|.. ji :. -5j
1
I
-
D-
FIGURE 6.11
Streams sheet defines both the feed streams (see Figure 6.12), FL and FV, as well
as product streams, PI and P2.
.«B ftH" Ml
It
- ia Ese
K J . -»...- >. " I MUM
FIGURE 6.12
p l|yi|yTff!l ?lWTfffF!ffBIHTi]f!Bl;I K | l]
E* Vwm OVs Took Run PM Lb«v Window Help
-
I9J 'I
_
i r i i -i f"
_
i 'M i ibi
] kfH l«" " d 1 I aU id .
3 (|-/C0M»»ni« i
_
J Pw, |
-
«3
B«nl.
I* Slag* 2 pfw-.me
r
_
j PM
S'iMnRtnJii PieiM* lor f*tl ol CoUm (oprwiBll
a n»
.
r SUflep-eiuwAop |
P Column pt 'Me <*op ffllT
O IJi«Siix»i«
H.nii
0 EO Irvul
0 Spec Grtt
Pa»i
1 RDCOLUMN
Q Srtup
lei Scwci r
"
J Varv
7 5l*jte OPurcwo***!)
Sep«alo(s [ H»al E'-changen ] Coktfrrtt | Re-scloit Pie*ii«e Channeli j MaMpUaloc. | ScW: | UswMotW; |
CV 0Foltas\AipennAlt.1
| -UCKy 6 Mno-.oliVoftl |[ Atpen Pkrt - Simutali..
FIGURE 6.13
Filling out Specifications and Holdups sheets we have two windows as shown in ,
- if**!
fit EcW Rin PW Lfcraiy Wrrfow Help
i i i i isi
i*7 Haft !HtTr8.sd <<||ii zi»li=JUH
/Specd»GaHon»| HoWjpc ] Re-Jdence T/nm j Coovoiw |
RsjJ: Reacfior. rwe:
(§ EOVanafcist Endng Reactor. 10 leaCton Chem:)iylti
f
aape JSOe
9 SpccGicmm 4 Rl
_
j PW:
Q S 'team RejJt *
ja pump
O Se«»
0 Peiiair.M e Cu
0 Um Suboine
0 bkK Opborc
Hen*.
(9 COVmMei
0 EOlnpti
Q Sc«c G-cmk
Seeen.Ren*i
) RDCOLUMN
9 Sew
OeivSoect
V«y
HeefenCoctoi
PtJTOWOUKll
Decertefi
Enoercc
'
STREAMS caw ypgi v** _
1
C1- BFob»: iO«PV* 1'
FIGURE 6.14(a)
DYNAMICS AND CONTROL OK PRESSURE DRIVEN PROCESSES -
297
- i - -- nu Todi nw .)
"
-
LT -I t-'h f 1 I >im!
' i lai
J POX.
l
_
5- i
*
1
IS
Pw Pt-ir.
Pka 11 1 HUM
FIGURE 6.14(b)
Select Pack Rating under RDCOLUMN of Blocks folder. Creating a new ID, T,
and specifying the packing section as well as packing characteristics, we obtain
Figure 6.15.
.
f » Vmm Da* Ictb Hkfi FV» Lt»v Wnfan H4>
Eg 1
1 3
"
3 ShWWw
XDCOUJMN
'
J
O «>U«(. v*-
T
FIGURE 6.15
298 PROCESS SIMULATION AND CONTROL USING ASPEN
f
200 (see Figure 6.16).
SJ MM
Pwlwm*xe Curve; .
atic jVAIgMhmj Jtuttvcai] Diaywiici |
0 Usci Subroulnt
Bbrt Optwni
BeaJi
1
« E0 Input
_
| Pom
SireamBesuM
ij iDCOLUMW IrJbakzatoiraelhdd
f
o Setup
DesgnSpeci
D Vary UgwMgwd c m* tpKir.g metl d
9 Healed Coofei:
D Punpaiixridi
CJ Dec»nie.j
S Ethaenaei
o itadorn
f
i Conden-ti Haivei
_
Rebaln Hcurvei
_ J IraySnIno
_
1 Tray Rarng
_
J
-
a Pack Pal
<J 1
® EtfrWei
ConvDisence
2
"
HUtod
STBEflMS
FffHsfc. p>G
FIGURE 6.16
In the next, dick Dynamic tab wader Blocks/RDCOLUMN The design specifications .
of the reflux drum and sump are reported in Figures 6 17(a) and (b). .
-
l»lxl
'
Fir EA View Data Took Run PW Lbmy WMm H«t>
I Dynamic
PetJi 1
EOVmabtei
O EOlnpul 3
© SpecGnxo:
RDCOLUMN
Total kwl»«k4MfrKiBn P
Condorao Hcu
BetwteiHcuve;
CJ TiayRama
SIREAMS ! PW
isimmsh
"
M * .;j
- . row. i?. ; 1619
FIGURE 6.17(a)
DYNAMICS AND CONTROL OF PRESSURE-DRIVEN PROCESSES 299
0 SMC Orn*.
j
1
_
I Cfa L St»ei
_
j CortJe HM-«i
4 PMt i*";
f
I 1
0 Cw-waenu
"
5 IfaprflMH
(THKMIS
-
S -
FIGURE 6.17(b)
- 5>
Fh £* V*- luco Rir> Ltto Wi<*>. K-t
j - :sec;
-
-
_
>
2 D~
[T W T|W| i
C - )f -m3e:Au» Pi= II
S'-i r-
FIGURE 6.18
300 PROCESS SIMULATION AND CONTROL USINd ASRKN1
Hit Next icon to open the Reactions folder. For the forward reaction (Reaction No 1) .
and the backward reaction (Reaction No. 2), the stoichiometric coefficients and
exponents are defined under Kinetic' Reaction type in the two sheets as shown in
'
MBj _U *J -'IfeM M 3 1 M _1 ±J J J
ProA-tli
CetKotn 1 t**rm*
MECH 1 l
ISI »
-
FtMoicra
f " Hi SeiBa> | StoWOf | H** ttfhwigsn [ CtWn j RwCCi- Pi«b.b OianiiBn | Mjr xWa! | Soldi 1 UaaMod-h
-
o
"
FIGURE 6.19(a)
RWMtt
rt. .
r.-.,
MICE 1
J UMfl
t HHi t*e»< I bWltoi I " f Kl«ro». | DAim | PmOo. Pimm* OunoMi | UitvJ**, | ioW, ] u«UuM<
f
FIGURE 6.19(b)
DYNAMICS AND CONTROL OF PRESSURE DRIVEN PROCESSES -
301
The Power law kinetdc data for both the reactions are provided in Figures 6 20(a) and (b). .
_
l r LL i
_ _ - - _
fV ! .mi 1 mi
3 :
us 5"'
o«flto.M1[iB-.p»*ii 1 mi - - .
FIGURE 6.20(a)
|jLi sial -J
,
I r I M- fj I l l
"
- 3
I en
It U»M«<PmwL«»
O Bkt* On
MTBE MEOH. iB
>
UAod
.
jo win
-
5!
FIGURE 6.20(b)
i r i i i i m i w 1 1
a 3 j LJ ajfzi 3xJ Nil
Jfl Set* I
'V Gompooenli
i
J/J PlOpMlK)
"
MkM
3 Slre rn r bte|
Re«ctioni
i/\ C onveigence
1
1" Zi 1' zi I Zl 1 Zi {
t
_
j C iv Options Volume Flow ciWhi 5922 949 146 743 97 974 146 743 97 971 29 549
-
ED Lotiv Opl<on:
40136 4053 10 396 41 770 .
10 396 41 770 40150
© DM0 Baiic Mete Flow kmcWn
Q DM0 Adv MEOH 711,300 100 007 0 605 100 607 0605 711 300
© LSSQPBawc
IB 707 465 93 059 4 517 93 059 4 517
f/ LSSQPAA*
Tea. NB 1257,715 1231 290 26 425 1231 290 26 425
_
j Conveiaence MTBE 0635 609 253 0 635
. 609 253
Conv Order
' itcteFcac
i Sequence
j FlomhBeting Options
_
MEOH 1000 0071 344 FfM 0071 944 PPM 1000
_
j ModelAnalysrs Tooli IB 0 360 0 065 01X17 0 065 0 007
/\ EG C&nliguration
'
Minut/SplmEiii | Sep««loi! | Heat EKchans«n ] Cokmm | ReMloii | Pressuie Changeis | MampuWoi! | Soidt | Usei Modds
-
a i K A
Matenal
'
STREAMS 1
Mw
ki
FSpB
r:,L.
SSpH
cc-o
f
_
ChaplBii | -j] l>apla 6 MicKBOll Wold || Atpni Plm - E«_Pb1.., , '» Aipon Dyiamcs EK_Pb | « _
1126
FIGURE 6.21
level controllers (LCI and LC3) and one pressure controller (PC2). Each of these
controllers has an operator set point (SP), a process variable (PV), also known as
controlled variable, and a controller output (OP), also called as manipulated variable,
whose values are obtained from the Aspen Plus simulation. These control structures
also have their own tuning parameters, and so on, suggested by Aspen Dynamics.
However, there is a scope to modify (or remove) the controller and its related items.
The Aspen generated control loops defined below should be used in the closed-loop
study of the prescribed catalytic distillation column.
J SAND QNTROL OF PRF.SSURE-Drivkn PRnrPgg,Q t 303
HSJ-Jl-(»>.
1 nr- <nm
teTH
J
1 ,
-,.,
FIGURE 6.22
Loop 1
Controller: LCI
Type of controller: proportional (P)-only (since reset time is very large)
Controlled variable, liquid level in the reflux drum
Manipulated variable: distillate (DIS) flow rate (percentage opening of valve CV2)
Controller action: direct
Use all default data, except proportional gain of 2 (suggested by Luyben, 2004)
Loop 2
Controller: PC2
Type of controller: proportional integral (PI)
Controlled variable: top stage pressure
Manipulated variable: condenser heat removal
Controller action: reverse
Use all default data (suggested by Luyben, 2004)
the condenser heat removal and P denotes the pressure to be controlled. Assuming
direct control action the controlle r equation can be rewr itten for Aspen Dynamics
Qr = - 47 48 - Kc {PSP - P), where 47.48 is the bias signal (Vr.s' "s
,
**
e from stea dy state position,
.
-'
gn indicates heat removal (cooling operation). If we mov
'
< is dear that when pressure (P) increases, the error (PS/. P) va lue becomes negahv..
d ultimately
.
,
the neRative vain., ofQc decreases. Originally, the negative value should
304 PROCESS SIMULATION AND CONTROL USING ASPEN
increase because if pressure increases, there is a need to increase the heat removal
rate. Therefore, our assumption is wrong and it should be reverse action in Aspen
Dynamics.)
Loop 3
Controller: LC3
f
Controller action: direct
Use all default data, except proportional gain of 2 (suggested by Luyben 2004) ,
For a chemical reaction with two reactants , the type of lowsheet depends on whether
f
we want to operate the catalytic distillation column with no-excess of either reactant or
excess reactant (Kaymak and Luyben, 2005). For a double-feed RD column if there is ,
any imbalance in the inflow of the two reactants ('excess reactant' case) the product
,
purity drops. This is because one of the reactants becomes excess and exits with the
product stream, and this stream would have to be further processed to purify the product
and recover the reactant for recycle. Obviously the 'excess reactant' lowsheet requires
,
f
at least two separating columns and is therefore more expensive. However it is easier
,
to control. On the other hand, the 'no-excess reactant' lowsheet has better steady state
f
economics but presents challenging control problems because of the need to precisely
balance the stoichiometry of the reaction.
Several control structures used to maintain the correct stoichiometric ratio of the
reactants have been proposed by researchers (e g Al-Arfaj and Luyben, 2000; 2002;
. .,
Wang et al., 2003). To meet this control objective the controller requires some type of
,
feedback of information rom within the process to indicate the accumulation or depletion
f
of at least one of the reactants. This can simply be done by the use of an internal
composition controller by manipulating the low rate of one of the resh feeds. There
f
are also other efficient control techniques (e.g. cascade control, inferential control)
,
reported for stoichiometric balancing (Wang et al. 2003). However, it is not practical
,
to simply ratio the two feed streams as has been proposed in some of the literature
,
papers. Flow measurement inaccuracies and feed composition changes doom to failure
DYNAMICS AND CONTROL OF PRESSURE-DRIVEN PROCESSES 305
any ratio controller that does not somehow incorporate information about compositions
inside the reactive system and feed this information back to adjust fresh feed.
For the concerned distillation column the methanol composition is controlled on
,
10* stage by the adjustment of the methanol fresh feed. The butene feed rate is low
f
controlled. It is worthy to mention that manipulating the methanol feed to control an
internal methanol composition is preferred when the butene feed coming rom the
f
upstream units is not ree to be adjusted. If this is not the case then alternatively the
f ,
f
diagram, place the control block near to CVl block and rename it as CC4. In the next,
expand Stream Types and use ControlSignal icon to complete the CC4 configuration,
shown in Figure 6.23. Chapter 5 presents a detail of how to configure a control structure
in Aspen Dynamics.
i a
Urn . -
. - B-
FIGURE 6.23
A little detail of the composition control loop for methanol feed is demonstrated
below .
Loop 4
Controller CC4
.
Type of controller: PI
Controlled variable: liquid phase mole fraction of MeOH on Stage 10
Manipulated variable: resh methanol (FL) low rate (percentage opening of v alve
f
CVl)
Controller action: reverse
306 PROCESS SIMULATION AND CONTROL USING ASPLN
Before executing the simulation run, it is customary to have a look on the data
sheet. For this, double-click on CC4 control block and then press Configure knob in the
faceplate to open the Configure dialog box. As mentioned in Chapter 5 it is wise to ,
click on Initialize Values button. Still one doubt is there: is the value of process variable
(PV) displayed same with the steady state liquid phase concentration of MeOH on
Stage 10 obtained in the Aspen Plus simulation? Be sure about it choose Blocks/ ,
RDCOLUMN/Profiles with opening the Aspen Plus simulation file. Then select 'Liquid'
in the View ield in Compositions sheet and obtain the table shown in Figure 6 24, . with
f
liquid mole fraction of MeOH on 10th stage of 0.04886022. This value is identical with
that of PV in the Configure dialog box.
MBI : 1: 1 m|»l *l
ii
Jfl PUMP
RDCOLUMN
MEOH 0 MB HTBE
)a;o**55i
-
a*
STREAMS
0 w* ?y3>Wl«
.
I[lr53
FIGURE 6.24
The controller CC4 is tuned by trial-and-error approach and the parameter values
have been chosen as:
RD column.
Top section
In addition to the LCI and PC2 control structures, the distillate composition can be
controlled by manipulating the reflux low rate. In an alternative approach, along withf
f
the pressure control (PC2), we can control the reflux drum level by the manipulation o
the reflux rate and the distillate low rate can be adjusted by a ratio control law to give
f
a constant reflux ratio. In the present case, the former control scheme has been
incorporated for performance study.
DYNAMICS AND CONTROL OF PRESSURE-DRIVKN PROCESSES 307
Bottom section
In the bottom section of a distillation column, it is a common practice that either the
bottom product purity or the tray temperature near the bottom of the column, which
has a strong correlation with the product purity, is controlled at its desired value by
the manipulation of the reboiler heat duty. For the sample process, we have implemented
a composition control structure for product quality control.
As the CC4 control block has been connected, similarly we can incorporate the
other control structures discussed above with the distillation flowsheet. The window,
shown in Figure 6.25, includes a closed-loop scheme in which the MTBE purity is
controlled in the bottoms by adjusting the reboiler heat input and the methanol impurity
in the top is controlled by manipulating the reflux flow rate. As stated earlier, the
concentration of methanol on the reactive stage it is being fed to (Stage 10) is measured
and controlled by the manipulation of the fresh methanol feed rate. The butene flow
rate is low-controlled. The liquid levels in the reflux drum and the base of the column
f
are maintained by the distillate flow rate and the bottoms flow rate, respectively. The
condenser heat removal is manipulated to control the column pressure. All of the
structures are single-input/single-output (SISO) structures with PI controllers (P-only
on levels).
c 0 a *q
-
-
M
m 5
JC»r.Lt»*v TV)-HM*
m m
r«uv(.-,
.
m
ll .
I
FIGURE 6.25
Loop 5
Controller: FC5
Type of controller: PI
308 PROCESS SIMULATION AND CONTROL USING ASPEN
f
Manipulated variable: brake power (shaft power or brake power of motor or engine
required to drive a compressor)
Controller action: reverse
Loop 6
Controller: CC6
Type of controller: PI
Controlled variable: MTBE mole fraction in the bottoms
Manipulated variable: reboiler heat input
Controller action: reverse
Loop 7
Controller: CC7
Type of controller: PI
Controlled variable: MeOH mole fraction in the distillate
Manipulated variable: reflux rate (mass low)
f
Now the flowsheet is ready for closed-loop performance study Start the program
.
(F7) rom the Run menu or press 'Re-start Simulation' button on the Run Control toolbar.
f
block or stream, then right-click to point Forms and inally select the item that we
f
want to access.
DYNAMICS AND CONTROL OF PRESSURE-DRIVEN PROCESSES 309
EBB
-
A -
1 r-
A
»5fl
i
jf5
ill 1 4--
FIGURE 6.26
DT2 and DT3) connected with three composition controllers (CC4 CC6 and CC7). ,
Ot- [ZtK- F
o5 3
4
x-i
- -
mf
Lii
13 iJ
FIGURE 6.27
310 PROCESS SIMULATION AND CONTROL USING ASPEN
incorporate a dead time for a measured variable say methanol mole fraction on ,
Stage 10, highlight DTI block right-click on the block, point to Forms and then select
,
Configure to open the configure table. In the Value cell enter 0.25 min as a sensor ,
dead time. Follow the same approach for other two dead time blocks .
Here, we have used the proportional gain of 1 %/% and integral time of 20 min or
f
all composition controllers. The effects of disturbance in butene feed temperature have
been depicted in Figure 6.28.
mm
Fte View Took Wtxjow Heb
f
SfnUahon
r tt Tf Gi«i|oo5 j i; a* v»
_
' K
-
I phut !-:-»&-
o<- gl«
arfcad x se r-
i-r*
U>V>* | jjDlWS-MiCCToll.. I A«»»Pte-t..W. | y <iM»teot.inole«,,||k7 fapm D «o. ... « #r)0 IS26
FIGURE 6.28
Initially a step decrease (76.85 -> 650C at time = 8 hours) and subsequently a step
increase (65 -> 760C at time = 15 hours) have been considered in the simulation study.
The developed closed-loop process lowsheet responds satisfactorily under load variable
f
processes. Here, a case study has been conducted on a MTBE catalytic distillation
column using the pressure-driven dynamics. The complete process low diagram includes
f
a distillation column, a feed compressor, a feed pump and three control valves. In the
MTBE synthesis process, a bottom product containing high-purity MTBE and a top
product enriched with n-butene are obtained. To maintain the MTBE purity in the
bottoms stream, several control structures have been configured with the lowsheet in
f
DYNAMICS AND CONTROL OF PKKSSURK-DRIVKN PROCKSSKS .*J 1 1
Aspen Dynamics. All of the structures are SISO schemes with PI controllers (P-only
on levels). The controllers have been tuned by simply using heuristics. The proposed
closed-loop process provides satisfactory results under disturbance input and
measurement lag.
PROBLEMS
6 . 1 A binary mixture of ethanol and l-propanol enters a flash drum (Flash2) The
feed specifications are shown in Figure 6.29 with the process low diagram .
f
Liquid mixture
(UQ-MIX) CV2
Temperature = 90X
Pressure = 1,4 bar
Flow rate = 120 kmol/hr cCH liq-mix f»B-[fgiE
CV1
Component Mol fract
ethanol 06
1-propanol 04 .
-
(pF]->t'i-|pdt-uq1-o
CV3
The lash chamber operates at 90oC and 1.2 bar. The vertically placed separator
f
with a length of 2 m and diameter of 1 m has elliptical head type. All the control
valves have a pressure drop of 0.2 bar. Applying the RK-Soave thermodynamic
model as a base property method,
(a) simulate the lowsheet to obtain the product compositions,
f
(b) design the two control schemes to maintain the pressure and liquid level in
the lash chamber, and
f
irreversible reaction:
An isentropic compressor discharges the FEED stream that enters the RStoic
reactor at 2 bar The reactor runs at 260oC and 2 bar The control valve involves
.
08
.
.
'
b) observe the closed loop process response employing the flow controllers.
-
312 PROCESS SIMULATION AND CONTROL USING ASPECT
Pure ethylbenzene
Temperature = 260oC
Pressure = 1 bar
Flow rate = 100 kmol/hr
"
M ! [pptI-
o
cv
-
|feed|-1
COMPRESS RSTOIC
FEED FL P1 -
u
F1
CV1 CV2
PUMP
> <
F2
COMPRESS RCSTR
the liquid-phase reaction the inlet streams Fl and F2, are specified in Table 6.2.
, ,
TABLE 6.2
ibi configure the control schemes for maintaining the liquid level pressure and ,
REFERENCES |
Al-Arfaj. M A. and W L Luyben (2000) "Comparison of Alternative Control Structures
.
for an Ideal Two-product Reactive Distillation Column Ind. Eng. Chem. Res., 39,
"
pp 3298-3307.
Al-Arfaj. M A and W L. Luyben (2002) "Control Study of Ethyl fert Butyl Ether Reactive
,
Jacobs. R. and R Krishna (1993) "Multiple Solutions in Reactive Distillation for Methyl
.
.
,
m the Kinetic Regime for the Synthesis of n Butvl Acetate Ind Eng. Chem Re* . "
42. pp B182-5194.
This
Thispage
pageintentionally
intentionallyleft
leftblank
blank
Rs. 295.00
@ 2009 by PHI Learning Pnvate Limited, New Delhi. All rights reserved. No part of this book may
be reproduced In any form, by mimeograph or any other means, without permission in writing from
the publisher.
ISBN-978-81-203-3659-9
The export rights of this book are vested solely with the publisher.
Published by Asoke K. Ghosh, PHI Learning Private Limited, M-97, Connaught Circus,
New Delhi-110001 and Printed by Jay Print Pack Private Limited, New Delhi-110015.
Index
'
t'tnlvtir dialillation 28fi ,
flbemCad ,
. 1 rixpOlll'lltN, 2il7
316 INDKX
Flash 2, 3, 7 Peng-Robinson 60 ,
Pre-exponential factor 65 ,
HETP, 287
Property method 18, 32. 39 ,
HYSYSTM, 3
Pumparound circuits 149 ,
Report options, 15
MTBE column, 286 Requil, 54
Material STREAMS, 7 Results plot dialog box, 251
Measurement lags, 309, 310 Reverse acting control, 243
MELLAPAK, 287 RGibbs, 54
Methane, 93 RK-Soave, 28. 32
Model library, 5 RPlug. 54, 78
Molarity, 76 RStoic, 54, 55
Multi-input/multi-output, 243 Run status, 62
MultiFrac, 107 RYield, 54
Sep 1, 2, 7
Object manager, 179 Separators, 42
Operator set point, 247 Servo performance 275 ,
Single-inputysingle-output, 243
Pause at time, 251 Solver settings, 13
PENG-ROB, 140 SRK, 52
INDEX 317
f
Stream information. 18. 33
Stream table, 22
Styrene, 55 Vapour raction 210
,
f
SULZER, 287 Variable number 180 .
As Ihe complexilv of a plant integrated with several process units increases solving Ihe model structure with a large equation
,
set becomes a challenging task. To overcome this situation, various process flowsheet simulators are used. This book describes
the simulation, optimisation, dynamics and closed-loop control of a wide variety of chemical processes using the most popular
'"
commercial flowsheet simulator Aspen .
The book presents the Aspen simulation of a large variety of chemical units, including flash drum, continuous stirred tank reactor
(CSTR), plug flow reactor (PFR), petroleum refining column, heat exchanger, absorption lower, reactive dislittation, disiillation
train, and monomer production unit. It also discusses the dynamics and control of flow-driven as well as pressure-driven chemical
processes using Ihe Aspen Dynamics package.
KEY FEATURES
Acquaints Ihe students with the simulation of large chemical plants with several single process units.
*
Includes a large number of worked out examples ittustrated in step ay-step format for easy understanding of the concepts.
f
This book is suitable for the undergraduate and postgraduate students of chemical engineering. It will also be helpful to research
scientists and practising engineers.
THE AUTHOR
Amiya K. Jana received his B.E. degree in chemical engineering in 1998 from Jadavpur University, M.Tech. in chemical engineering
in 2000 from IIT Kharagpur, and Ph.D. in chemical engineering in 2004 from IIT Kharagpur.
Presently. Or. Jana is Assistant Professor at IIT Kharagpur. His areas of research include control system process intensification,
,
and modelling and simulation. He is also the author of ChemiesJ Process Mode/ting and Computer Simukuon published by
PHI learning.
lSBN:')7fl-flWD3-3l.S1-,1
Rs. 295.00
www.phlndia.com