Aspen Plus Urea Synthesis Loop Model
Aspen Plus Urea Synthesis Loop Model
Aspen Plus Urea Synthesis Loop Model
Aspen Plus®, Aspen Properties®, the aspen leaf logo and Plantelligence and Enterprise Optimization are trademarks
or registered trademarks of Aspen Technology, Inc., Bedford, MA.
All other brand and product names are trademarks or registered trademarks of their respective companies.
This document is intended as a guide to using AspenTech's software. This documentation contains AspenTech
proprietary and confidential information and may not be disclosed, used, or copied without the prior consent of
AspenTech or as set forth in the applicable license agreement. Users are solely responsible for the proper use of
the software and the application of the results obtained.
Although AspenTech has tested the software and reviewed the documentation, the sole warranty for the software
may be found in the applicable license agreement between AspenTech and the user. ASPENTECH MAKES NO
WARRANTY OR REPRESENTATION, EITHER EXPRESSED OR IMPLIED, WITH RESPECT TO THIS DOCUMENTATION,
ITS QUALITY, PERFORMANCE, MERCHANTABILITY, OR FITNESS FOR A PARTICULAR PURPOSE.
3 Process Description..............................................................................................3
4 Physical Properties...............................................................................................5
8 Conclusions ........................................................................................................13
9 References .........................................................................................................14
Contents iii
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 CO2 Stripping Process, which is a
popular and fast growing process for manufacturing urea.
The work 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
package has to be developed. This type of model is useful to analyze the plant
performance and to improve plant operation, including:
Energy saving studies to improve economics of the plant.
Studies of individual pieces of equipment with a view of increasing their
throughput and/or improving their performance.
Flowsheet modification for better plant operation.
Identifying bottlenecks.
As a basis for optimization study and for developing on-line control system
of the plant.
While this document describes the simulation of the Stamicarbon CO2
stripping process, the accurate results obtained support the applicability of
Aspen Plus and the data package to other urea processes.
1 Introduction 1
2 Components and Units of
Measurement
F01
FLASH2
R01 S15
E02
E03
S13
A01
RGIBBS
M01
MIXER
S01
E01
S07
SP1
A02 S12
RSTOIC SEP
3 Process Description 3
The process is described as follows:
Feed CO2 gas (S07) is fed to the CO2 stripper, E01, to strip the urea
solution coming from the reactor. In the stripper, ammonium carbamate
decomposes, liberating more NH3 and CO2 to be stripped out. Heat is
supplied on the shell side of tubes by condensing 285 psig steam while the
urea solution falls inside the tubes countercurrently down past the rising
CO2 stripping gas. The outlet liquid solution from the stripper is rich in
urea and goes to the downstream section for urea purification.
In the adiabatic urea reactor, R01, an aqueous solution of NH3 and CO2
(much in the form of ammonium carbamate) and vapors flow upward
through 8 stages of reactor volume to minimize back-mixing and provide
enough residence time for urea formation. Remaining gases condense
and carbamate decomposes in the reactor to provide heat for the slightly
endothermic reaction of carbamate to urea. The urea solution (S06)
overflows from the top of the reactor back to the Stripper E01, while the
unreacted gases (S05) pass out the top of the reactor.
These unreacted gases pass to the Scrubber, E03, where recycled
carbamate solution from the Evaporation/Recirculation section (S13) is
passed over the top of a packed bed and fills the tube side section of
vertical tubes. The gases rise up through the tubes and pass up through
the packed section before leaving the top of the vessel. The contact with
the carbamate solution absorbs the unreacted NH3 and CO2, while the
inert gases of O2, N2, and others (S15) vent out from the top. On the
bottom, recirculated cooling water flowing inside tubes removes the heat
of absorption from the carbamate solution. Carbamate solution (S17)
overflows out of the vessel.
The solution (S17) together with the top vapor stream from the Stripper
(S08) is fed to the Carbamate Condenser, E02, through the use of an
ejector, where the ammonia feed (S01) serves as the pumping fluid.
Ammonium carbamate forms in this condenser. The vapor-liquid mixture
(S03) falls through tubes and the heat of reaction is removed by
generation of 50 psig steam on the shell side. The mixture goes to the
bottom of the reactor for urea production.
The bottom stream (S09) from the Stripper E01 is sent to the downstream
section to recover urea. The recycled stream after recovering the urea
(S13) is back to the H.P. Scrubber to complete the loop.
4 3 Process Description
4 Physical Properties
4 Physical Properties 5
predicted by the present AspenTech model with the Gorlovskii-Kucheryavyi
correlation. The present correlation clearly provides an accurate fit of the
data.
6 4 Physical Properties
Figure 3 presents a comparison between the present AspenTech model and
the bubble pressure data measured by Lemkowitz (1971, 1972 and 1977).
The model provides an accurate description of the data and, in particular,
provides an accurate description of the minimum in bubble pressure, which is
necessary for an effective description of the urea synthesis process.
200°C
250
Bubble Pressure (bar)
200 180°C
150 160°C
100 140°C
50
0
0.65 0.70 0.75 0.80 0.85 0.90 0.95
Mole Fraction NH3 in Feed
4 Physical Properties 7
5 Chemical Reactions
Kinetics
There are two main reactions that take place in the urea synthesis process:
(1) 2NH3 + CO2 CARB
(2) CARB UREA + H2O
The first reaction, which takes place in the liquid phase, converts ammonia
and carbon dioxide into ammonium carbamate. This reaction is highly
exothermic and fast. Chemical equilibrium is readily reached under the
operating conditions in the reactor. The second reaction also takes place in
the liquid phase and is endothermic. Its rate is slow and equilibrium is
usually not reached in the reactor.
A user subroutine, USURA.F, was developed to include the reaction kinetics of
both reactions. USURA.F is used in the reactor simulations. Both forward and
reverse reactions were considered. The kinetics of Reaction 1 has been set to
be rapid so that equilibrium is effectively reached. Limited literature data
were used for the kinetics of Reaction 2.
The reaction kinetics has been formulated to approach the equilibrium
composition for large residence times. The equilibrium has been described in
terms of the fugacity coefficients since an equation of state is used as the
thermodynamic model. The equilibrium constant for Reaction 1, in terms of
mole fractions, is written as follows:
K1 exp
G0
CARB 2GNH 3 GCO 2
0 0
P NH
22
3 CO 2
0
CARB
RT P
Where:
T = Temperature
P = Pressure
x = Mole fraction vector
R = Gas constant
P0 = Reference pressure (= 1 atmosphere)
2 x
Rate1 k1 xNH x
3 CO 2
CARB
K1
xUREA xH 2O
Rate2 k2 xCARB
K2
The two rate expressions have been formulated so that they will necessarily
reach equilibrium at large residence times. The rate constant for Reaction 1
(K1) is set to a large value so that this reaction is essentially at equilibrium.
The rate constant for Reaction 2 determines the urea conversion in the
reactor. Only scant information is available to determine K2 and it is usually
best to adjust its value to fit plant data. A reasonable approximation for K2 is
the following:
100.*106
RT
k2 15. 108 e /V L
Where R=8314.3 and VL is the molar volume of the liquid.
The CO2 Stripper, E01, is of the falling film type, which was approximated by
a RADFRAC (multistage distillation) block with 10 stages. Heat is supplied to
the stages 2-9 to simulate the heat transfer from the tubes. The urea
solution, falling down on tube walls, is stripped off volatile NH3 by the
entering CO2 gas (S07). On each stage the model considers the equilibrium
of carbamate in the liquid as well as the VLE of the mixture. Note that the
kinetics of carbamate formation is large enough to ensure that chemical
equilibrium for the carbamate reaction is reached in each stage of the
RADFRAC block.
The Urea Reactor R01, 7.52 ft in diameter and 95 ft in length, is modeled with
an RPLUG block. The kinetics is provided by the user subroutine USURA in
the RPLUG block.
In the reactor, the exothermic carbamate reaction and the endothermic urea
formation reaction are taking place. The reactor is designed such that its
volume is big enough for the desired urea production.
The H.P. Scrubber is modeled using a 5-stage RADFRAC block. Heat is taken
out from the bottom stage. Similar to the Stripper, E01, on each stage the
model considers the equilibrium of carbamate in the liquid as well as the VLE
of the mixture.
The E02 H.P. CO2 Condenser is modeled with an RSTOIC block. A Design
Specification paragraph is included to monitor the specified reactor outlet
temperature of 183oC by adjusting the extent of the carbamate reaction in
E02.
This simulation is based on a closed-loop flowsheet. The downstream section
is approximated by using a SEP block to link the Stripper bottom urea solution
(S09) to the recycled carbamate solution (S13).
10 6 Simulation Approach
7 Simulation Results
The Aspen Plus run was made using Version 206.5. Some of the results are
shown below, and a simulation flowsheet with stream data is shown in Figure
4. The simulation results of this generic model are reasonable compared with
similar plant operations.
7 Simulation Results 11
Key Process Simulation results:
Equipment Variable Value Unit
12 7 Simulation Results
8 Conclusions
This urea process model has been developed using Aspen Plus Version
2006.5. This is a rigorous closed-loop model for the plant while the
recovery section is approximated by using a SEP model. The carbon
dioxide compression section is not included. From the results, it is shown
that the SR-POLAR property package used for simulation is appropriate.
The model may be used as a guide for understanding the process and the
economics, and also as a starting point for more sophisticated models for
plant designing and process equipment specifying.
For further refinement of the model, the following upgrades should be
made:
o The cooling water circuit for E03 and the low pressure steam circuit for
E02 should be implemented. This implementation is useful for energy
saving studies.
o E02 is simulated in this work by a RSTOIC model. However, a RPLUG
model is more suitable for the simulation of this equipment. However,
to do so, detailed equipment data for E02 are needed.
o The stripper E01 is a falling-film type equipment. It embodies VLE
equilibrium, mass transfer, reaction, and heat transfer. To rigorously
simulate this equipment, rate-based calculations (RateSep) should be
used. A special subroutine will be needed to incorporate the mass
transfer limitations of the falling-film. In this simulation, the RADFRAC
model was used with component efficiencies for NH3, CO2 and H2O as a
work around.
The accurate results obtained in the present simulation indicate that
Aspen Plus and the physical property data package will provide accurate
simulations of other urea processes.
8 Conclusions 13
9 References
14 9 References
[13] Lemkowitz, S.M.; Vet, E.; van den Berg, P.J., “A Phase Model for the
Gas-Liquid Equilibria in the Ammonia-Carbon Dioxide Water-Urea System in
Chemical Equilibrium at Urea Synthesis Conditions. II. Experimental
Verification,” J. Appl. Chem. Biotechnol., 27, 335-348 (1977).
[14] Lemkowitz, S.M.; Vet, E.; van den Berg, P.J., “Phase Behavior in the
Ammonia-Carbon Dioxide System at and Above Urea Synthesis Conditions,” J.
Appl. Chem. Biotechnol., 22, 727-737 (1972).
[15] Peneloux, A.; Rauzy, E., “A Consistent Correction for Redlich-Kwong-
Soave Volumes,” Fluid Phase Equilibria, 8, 7-23 (1982).
[16] Schwartzentruber, J.; Renon, H., “Extension of UNIFAC to High Pressures
and Temperatures by the Use of a Cubic Equation of State,” I&EC Research,
28, 1049-1055 (1989).
[17] Soave, G., “Equilibrium Constants from a Modified Redlich-Kwong
Equation of State,” Chem. Eng. Sci., 27, 1196-1203 (1972).
9 References 15