Hysys 1
Hysys 1
Hysys 1
c om
Contents
1Introduction...............................................................................1
3Process Description....................................................................2
4Physical Properties.....................................................................3
!"imulation #pproach................................................................... $
$"imulation Results......................................................................$
%Conclusions................................................................................&
References.................................................................................1'
1 Introduction
This document describes the steady-state Aspen Plus ® model of the high-pressure synthesis loop of a
urea plant, with a capacity of about 1,100 metric tons of prilled urea per year.
This simulation is based on the Stamicarbon ! " Stripping Process, which is a popular and fast growing
process for manufacturing urea.
The wor# demonstrates the capability of Aspen Plus to rigorously model the urea synthesis process. The
modeling is complicated due to the formation of ammonium carbamate, an intermediate product for
which a special property pac#age has to be de$eloped. This type of model is useful to analy%e the plant
performance and to impro$e plant operation, including&
• Studies of indi$idual pieces of e(uipment with a $iew of increasing their throughput and)or
impro$ing their performance.
• +dentifying bottlenec#s.
• As a basis for optimi%ation study and for de$eloping on-line control system of the plant.
hile this document describes the simulation of the Stamicarbon ! " stripping process, the accurate
results obtained support the applicability of Aspen Plus and the data pac#age to other urea processes.
Components
A small amount of biuret 8 "9!": is produced during the synthesis. +n this wor#, the biuret production
is not considered, but the component can easily be added if re(uired.
The pure component properties of all the components e2cept ammonium carbamate can be retrie$ed
from the Aspen Plus databan#s. Special efforts were made to incorporate pure component properties of
ammonium carbamate in the simulation.
A special contribution of this wor# is the de$elopment of a physical property model to describe the
simultaneous physical-chemical e(uilibrium occurring in the urea-synthesis process.
;etric units are used in this wor#, e2cept that temperature unit is o, pressure unit is #g)cm " 8abs:,
enthalpy flow unit is mmcal)hr, and mole flow unit is #mol)day.
! Process Description
A simulation flowsheet of the synthesis loop is shown in *igure 1.
401FLA!"*'S< #T$I%A0"EPS1&'();9
S05
*01
FLA!"
401 S19
'0"
'0
S1
A01
#&I''
;01
(I)E#
S01
'01
S0<
SP1
A0" S1"
#T$I% EP
The high-pressure loop is operated at around 151 #g)cm " 8abs: and consists of the following #ey pieces
of e(uipment&
%&uipment Purpose
401 3rea 4eactor, where ammonium carbamate is dehydrated to urea
'01 .P. !" Stripper, where the bul# of the unreacted carbamate from the reactor
effluent is decomposed by stripping with !" gas and with heat input.
'0" .P. ondenser, where the gaseous !" and condense and react to form
ammonium carbamate.
'0 .P Scrubber, where the recycled carbamate solution from the downstream low-
pressure section is used to absorb unreacted gases from the reactor.
The model for the thermodynamic properties of the -!"-"!-34'A-A46-"-!" system is based
upon the S4-P!@A4 model within Aspen Plus 8Soa$e, 1?<" Penelou2 and 4au%y, 1?="
Schwart%entruber and 4enon, 1?=?:. The model uses an e(uation of state and is thus suitable for the
high-pressure, high-temperature conditions of urea synthesis. *urther, the model contains e2tensions
that enable an accurate description of the phase and chemical e(uilibria, the density and the other
thermodynamic properties 8e.g., enthalpy: of this system.
e chose our approach to the modeling of the thermodynamic properties after studying the pre$ious
modeling attempts in the literature and analy%ing the a$ailable data. *rB>ac(ues 81?5=:, Cawasumi
81?9", 1?9 and 1?95: and @em#owit% 81?=0: de$eloped chemical and thermodynamic models by
postulating reactions for urea formation and $arious simplified assumptions for the phase non-ideality.
The simplifying assumptions do not permit an accurate and general model for the chemical-
thermodynamic properties. 6ernadis et al. 81?=?: and +sla et al. 81??: de$eloped impro$ed theoretical
models by including ionic species and describing the nonideality of the li(uid phase by a modified
3+D3A model. e belie$e that under the high temperatures 8170 to "00°: and the relati$ely low
water concentrations of urea synthesis, the e2tent of ioni%ation will be small. *urther, modern e(uations
of state such as the S4-P!@A4 model are well suited to the description of the thermodynamic properties
of nonideal systems at high pressures and temperatures. Thus we ha$e chosen to use the S4-P!@A4
model as the physical-property option.
e$elopment of a data pac#age for this system is difficult since most of the data are only a$ailable as
combined physical and chemical e(uilibria. isassociation pressure data are a$ailable for ammonium
carbamate 8Eanac#e, 1?0:, which is the pressure at a specified temperature where ammonia and
carbon dio2ide are in e(uilibrium with the condensed-phase ammonium carbamate. ata are also
a$ailable for the e(uilibrium con$ersion of defined mi2tures of -!"-"! to urea 8Cawasumi, 1?9",
1?9 and 1?95 +noue, 1?<":. *inally, bubble pressures ha$e been measured for defined mi2tures of
-!"-"! at chemical e(uilibrium. The Aspen Plus data analysis capabilities 84S and ATA-*+T:
ha$e been used to obtain a simultaneous good fit of these $aried and comple2 sets of data.
Forlo$s#ii and Cucherya$yi 81?<?: ha$e analy%ed the urea con$ersion data reported by many authors
and ha$e produced an empirical correlation that is recogni%ed to be accurate. *igure " compares the
e(uilibrium con$ersions predicted by the present AspenTech model with the Forlo$s#ii-Cucherya$yi
correlation. The present correlation clearly pro$ides an accurate fit of the data.
*igure presents a comparison between the present AspenTech model and the bubble pressure data
measured by @em#owit% 81?<1, 1?<" and 1?<<:. The model pro$ides an accurate description of the data
and, in particular, pro$ides an accurate description of the minimum in bubble pressure, which is
necessary for an effecti$e description of the urea synthesis process.
;easured and
omparison of alculated
ASP' P@3S 6ubble
'!SPressures
to ata ofin@em#owit%
the -!"
et al.System
81?<2:
00
"00G
"90
:
r 1=0G
a"00
b
8
e
r
u
s
s
e190 170G
r
P
e
lb
b
u100 150G
6
90
0
0.79 0.<0 0.<9 0.=0 0.=9 0.?0 0.?9
;ole *raction in *eed
i"ure !# Measured and Calculated 5u((le Pressures in the N6! C*2 System
There are two main reactions that ta#e place in the urea synthesis process&
The first reaction, which ta#es place in the li(uid phase, con$erts ammonia and carbon dio2ide into
ammonium carbamate. This reaction is highly e2othermic and fast. hemical e(uilibrium is readily
reached under the operating conditions in the reactor. The second reaction also ta#es place in the li(uid
phase and is endothermic. +ts rate is slow and e(uilibrium is usually not reached in the reactor.
A user subroutine, 3S34A.*, was de$eloped to include the reaction #inetics of both reactions.
3S34A.* is used in the reactor simulations. 6oth forward and re$erse reactions were considered. The
Aspen Plus Example Library .
Proprietary Information of AspenTech
Unauthorized duplication or distribution strictly prohibited without prior written permission
http://slide pdf.c om/re a de r/full/188882730-sinte sis-de -ure a -c on-a spe n 6/11
5/20/2018 188882730 Sinte sis de Ure a Con Aspe n - slide pdf.c om
#inetics of 4eaction 1 has been set to be rapid so that e(uilibrium is effecti$ely reached. @imited
literature data were used for the #inetics of 4eaction ".
The reaction #inetics has been formulated to approach the e(uilibrium composition for large residence
times. The e(uilibrium has been described in terms of the fugacity coefficients since an e(uation of
state is used as the thermodynamic model. The e(uilibrium constant for 4eaction 1, in terms of mole
fractions, is written as follows&
T - Temperature
P - Pressure
2 - ;ole fraction $ector
4 - Fas constant
P0 - 4eference pressure 8I 1 atmosphere:
0 0
Gi - +deal-gas Fibbs *ree energy of component i at T, P
φi - *ugacity coefficient of component i at T, P, 2
x CARB
K1 =
"
x NH x CO "
The rates for 4eactions 1 and ", in units of #mol)s)m , are as follows&
=k
xCARB
Rate1 x NH
1
"
+
xCO " −
K1
x x
Rate" = k xCARB − UREA H
"
"O
K "
The two rate e2pressions ha$e been formulated so that they will necessarily reach e(uilibrium at large
residence times. The rate constant for 4eaction 1 8C1: is set to a large $alue so that this reaction is
essentially at e(uilibrium. The rate constant for 4eaction " determines the urea con$ersion in the
reactor. !nly scant information is a$ailable to determine C" and it is usually best to ad>ust its $alue to fit
plant data. A reasonable appro2imation for C " is the following&
100.J107
−
RT
k" = 19. J 10 =
e
)v
L
The !" Stripper, '01, is of the falling film type, which was appro2imated by a 4A*4A 8multistage
distillation: bloc# with 10 stages. eat is supplied to the stages "-? to simulate the heat transfer from the
tubes. The urea solution, falling down on tube walls, is stripped off $olatile by the entering ! " gas
8S0<:. !n each stage the model considers the e(uilibrium of carbamate in the li(uid as well as the /@'
of the mi2ture. ote that the #inetics of carbamate formation is large enough to ensure that chemical
e(uilibrium for the carbamate reaction is reached in each stage of the 4A*4A bloc#.
The 3rea 4eactor 401, <.9" ft in diameter and ?9 ft in length, is modeled with an 4P@3F bloc#. The
#inetics is pro$ided by the user subroutine 3S34A in the 4P@3F bloc#.
+n the reactor, the e2othermic carbamate reaction and the endothermic urea formation reaction are
ta#ing place. The reactor is designed such that its $olume is big enough for the desired urea production.
The .P. Scrubber is modeled using a 9-stage 4A*4A bloc#. eat is ta#en out from the bottom
stage. Similar to the Stripper, '01, on each stage the model considers the e(uilibrium of carbamate in
This simulation is based on a closed-loop flowsheet. The downstream section is appro2imated by using
a S'P bloc# to lin# the Stripper bottom urea solution 8S0?: to the recycled carbamate solution 8S1:.
+ Simulation ;esults
The Aspen Plus run was made using /ersion "07.9. Some of the results are shown below, and a
simulation flowsheet with stream data is shown in *igure 5. The simulation results of this generic model
are reasonable compared with similar plant operations.
3 Conclusions
1. This urea process model has been de$eloped using Aspen Plus /ersion "007.9. This is a rigorous
closed-loop model for the plant while the reco$ery section is appro2imated by using a S'P model.
The carbon dio2ide compression section is not included. *rom the results, it is shown that the S4-
P!@A4 property pac#age used for simulation is appropriate.
". *or further refinement of the model, the following upgrades should be made&
a: The cooling water circuit for '0 and the low pressure steam circuit for '0" should be
implemented. This implementation is useful for energy sa$ing studies.
b: '0" is simulated in this wor# by a 4ST!+ model. owe$er, a 4P@3F model is more suitable
for the simulation of this e(uipment. owe$er, to do so, detailed e(uipment data for '0" are
needed.
c: The stripper '01 is a falling-film type e(uipment. +t embodies /@' e(uilibrium, mass transfer,
reaction, and heat transfer. To rigorously simulate this e(uipment, rate-based calculations
84ateSep: should be used. A special subroutine will be needed to incorporate the mass transfer
limitations of the falling-film. +n this simulation, the 4A*4A model was used with component
efficiencies for , !" and "! as a wor# around.
. The accurate results obtained in the present simulation indicate that Aspen Plus and the physical-
property data pac#age will pro$ide accurate simulations of other urea processes.
;eferences
K1L 6ernadis, ;. ar$oli, F. Santini, ;., M3rea--!"-"! /@' alculations 3sing an '2tended
3+D3A '(uation,N *luid Phase '(uilibria, 9, "0<-"1= 81?=?:.
K"L *rB>ac(ues, ;., MTheoretical 6asis of the +ndustrial Synthesis of 3rea,N him. +nd., 70, ""-9 81?5=:.
KL Forlo$s#ii, .;. Cucherya$yi, /.+., M'(uation for etermination of the '(uilibrium egree of !"
on$ersion uring Synthesis of 3rea,N Oh. Pri#l. Chim., 9, "95=-"991 81?<?:.
K5L +noue, S. Canai, C. !tsu#a, '., M'(uilibrium of 3rea Synthesis. +.N 6ull. hem. Soc. Eapan, 59,
1?-159 81?<":.
K9L +sla, ;.A. +ra%o(ui, .A. Fenoud, .;., MSimulation of a 3rea Synthesis 4eactor. 1.
Thermodynamic *ramewor#,N +' 4esearch, ", "77"-"7<0 81??:.
K7L Eanec#e, '., O. 'lectrochem., 7, 759 81?0:.
K<L Cawasumi, S., M'(uilibrium of the !"--3rea-"! System 3nder igh Temperature,N 6ull.
hem. Soc. Eapan, "9, ""<-"= 81?9":.
K=L Cawasumi, S., M'(uilibrium of the !"--3rea-"! System 3nder igh Temperature and
Pressure. ++. @i(uid-/apor '(uilibrium in the @oading ;ole 4atio of " to !",N 6ull. hem. Soc.
Eapan, "7, "1=-""< 81?9:.
K?L Cawasumi, S., M'(uilibrium of the !"--3rea-"! System 3nder igh Temperature and
Pressure. +/. @i(uid-/apor '(uilibrium in the Presence of '2cess Ammonia or arbon io2ide,N
6ull. hem. Soc. Eapan, "<, "95-"9? 81?95:.
K10L Cotula, '.A., M/apor-@i(uid '(uilibrium of the -!"-"!-3rea System at 'le$ated Pressure,N
E. hem. Technol. 6iotechnol., 1, 10-110 81?=1:
K11L @em#owit%, S.;. Foedegebuur, E. $an den 6erg, P.E., M6ubble-Point ;easurements in the
Ammonia-arbon io2ide System,N E. Appl. hem. 6iotechnol., "1, ""?-"" 81?<1:.
K1"L @em#owit%, S.;. $an 'rp E.. 4e#ers, .;. $an den 6erg, P.E., MPhae '(uilibria in the
Ammonia-arbon io2ide Systems at and Abo$e 3rea Synthesis onditions,N E. Appl. hem.
6iotechnol., 0, =9-101 81?=0:.
K1L @em#owit%, S.;. /et, '. $an den 6erg, P.E., MA Phase ;odel for the Fas-@i(uid '(uilibria in the
Ammonia-arbon io2ide ater-3rea System in hemical '(uilibrium at 3rea Synthesis
onditions. ++. '2perimental /erification,N E. Appl. hem. 6iotechnol., "<, 9-5= 81?<<:.
K15L @em#owit%, S.;. /et, '. $an den 6erg, P.E., MPhase 6eha$ior in the Ammonia-arbon io2ide
System at and Abo$e 3rea Synthesis onditions,N E. Appl. hem. 6iotechnol., "", <"<-<< 81?<":.
K19L Penelou2, A. 4au%y, '., MA onsistent orrection for 4edlich-Cwong-Soa$e /olumes,N *luid
Phase '(uilibria, =, <-" 81?=":.
K17L Schwart%entruber, E. 4enon, ., M'2tension of 3+*A to igh Pressures and Temperatures by
the 3se of a ubic '(uation of State,N +' 4esearch, "=, 105?-1099 81?=?:.
K1<L Soa$e, F., M'(uilibrium onstants from a ;odified 4edlich-Cwong '(uation of State,N hem. 'ng.
Sci., "<, 11?7-1"0 81?<":.