Aspen Plus

Aspen Plus Ammonia Model

Version Number: V7.0 July 2008

Copyright 2008 by Aspen Technology, Inc. All rights reserved. Aspen Plus, Aspen Properties, the aspen leaf logo and Plantelligence and Enterprise Optimization are trademarks or registered trademarks of Aspen Technology, Inc., Burlington, 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. Aspen Technology, Inc. 200 Wheeler Road Burlington, MA 01803-5501 USA Phone: (1) (781) 221-6400 Toll Free: (1) (888) 996-7100 URL: http://www.aspentech.com

Contents
1 Introduction .........................................................................................................2 2 Components .........................................................................................................3 3 Process Description..............................................................................................5 4 Physical Properties...............................................................................................8 5 Chemical Reactions ..............................................................................................9 Reforming Unit .............................................................................................. 9 Carbon Monoxide Conversion ..........................................................................12 Carbon Dioxide Removal ................................................................................13 Methanation Unit...........................................................................................15 Synthesis Unit ..............................................................................................16 6 Simulation Approaches.......................................................................................18 7 Simulation Results .............................................................................................20 8 Conclusions ........................................................................................................21 9 References .........................................................................................................22

1 Introduction

This model simulates an Ammonia Production process using natural gas as a feedstock. This model includes the following features: A set of chemical species and property parameters for this process. Typical process areas including: Desulfurization, Reforming Unit, Carbon Monoxide Conversion, Carbon Dioxide Removal, Methanation Unit, Synthesis Unit, Refrigeration and the main streams connecting these units. Usability features such as an Excel file which allows the user to collect simulation results from the synthesis reactor. Definition of property model parameters with user data.

2 Components

The table below lists the components modeled in the ammonia plant model.
Component ID NH3 H2 N2 CH4 AR CO CO2 H2O O2 C2H6 C3H8 N-BUTANE I-BUTANE N-PENTAN I-PENTAN N-HEXANE N-HEPTAN SULFUR H2S H3O+ OHNH4+ NH2COOHCO3CO3-NH4HCO3S NH4HCO3 Type CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV CONV SOLID CONV Component name AMMONIA HYDROGEN NITROGEN METHANE ARGON CARBON-MONOXIDE CARBON-DIOXIDE WATER OXYGEN ETHANE PROPANE N-BUTANE ISOBUTANE N-PENTANE I-PENTANE N-HEXANE N-HEPTANE SULFUR-8-ATOMIC-GAS HYDROGEN-SULFIDE H3O+ OHNH4+ CARBAMATE HCO3CO3-AMMONIUM-HYDROGEN-CARBONATE AMMONIUM-HYDROGEN-CARBONATE Formula NH3 H2 N2 CH4 Ar CO CO2 H2O O2 C2H6 C3 H 8 C4H10 C4H10 C5H12 C5H12 C6H14 C7H16 S8 H2S H3O+ OHNH4+ NH2COOHCO3CO3-2 NH4HCO3 NH4HCO3

Every component is not included in all the plant sections in order to keep the model as simple as possible. The higher hydrocarbons (C2-C7) are only

included in the reformer section. The CO2 Removal section only includes ions and electrolytes. The synthesis and refrigeration sections only include the components H2, N2, Ar, CH4, NH3 and H2O.

3 Process Description

Ammonia is produced basically from water, air, and energy. The energy source is usually hydrocarbons, thus providing hydrogen as well, but may also be coal or electricity. Steam reforming of light hydrocarbons is the most efficient route, with about 77% of world ammonia capacity being based on natural gas.

Figure 1 Block diagram of the steam/air reforming process.

Ammonia plants using natural gas as a feedstock include the following units:

Natural Gas Desulfurization

The natural gas is delivered as dry gas containing a maximum of 40 ppm by weight of sulfur, which is a poison for the reformer catalyst. The desulfurization unit reduces the sulfur content to about 5 ppm by hydrogenating it to hydrocarbons and hydrogen sulfide and then absorbing the hydrogen sulfide in zinc oxides.

Reforming Unit
This unit contains two sections, one is primary reforming, and another is secondary reforming. The desulfurized hydrocarbon feed is reformed to hydrogen and carbon oxides in the presence of steam in the primary reformer, and additionally with hot air in the secondary reformer. The reformed gas contains about 0.3 vol% CH4.

Carbon Monoxide Conversion

In the CO-shift conversion, the major part of the CO contained in the reformed gas is catalytically converted to CO2 in two catalytic stages, the first at high temperature and the second at low temperature.

Carbon Dioxide Removal

The carbon dioxide is removed from the converted gas in the CO2 removal. CO2 is captured by NH3 and it will generate ammonium hydrogen carbonate as byproduct. The purified gas with about 0.1 vol% CO2 is called synthesis gas.

Methanation Unit
Even small quantities of CO (0.1 vol%) and CO2 (0.3 vol%) are poisons for the ammonia synthesis catalyst. Thus they are reacted to methane over a nickel catalyst. The residual content of CO + CO2 is less than 10 ppm.

Synthesis Unit
The synthesis gas is pressurized by a centrifugal compressor to approximately 300 bar and hydrogen and nitrogen are catalytically converted to ammonia.

Tailgas Scrubbing
The ammonia in the purge gas from the ammonia unit is recovered in the tailgas scrubbing unit and fed to a refrigeration unit. The treated purge gas is used as fuel for the primary reformer. The building blocks offered in this package do not model this section of the plant.

Refrigeration
The ammonia gas in the synthesis loop is liquefied by ammonia evaporation in the ammonia chiller and discharged as feed to the urea process and to ammonia storage. The building blocks offered in this package do not model this section of the plant.

4 Physical Properties

Equation-of-state models provide an accurate description of the thermodynamic properties of the high-temperature, high-pressure conditions encountered in ammonia plants. The Redlich-Kwong modification RKS-BM was chosen for this application. In most cases, the model parameters were taken from the Aspen physical property databanks. Special accuracy improvements were made to the RKS-BM model as follows: The pure-component properties of ammonia were modified to improve the vapor pressure and enthalpy of vaporization, which is important for the synthesis section. The binary interaction parameters of H2, N2, Ar and CH4 with ammonia were modified in order to obtain accurate solubilities of these species in ammonia, which is important for the synthesis and refrigeration sections. The binary parameters of the air system (N2, O2, Ar) were modified in order to enable an accurate simulation of the vapor-liquid equilibrium in the nitrogen plant.

Electrolyte NRTL method and RK equation of state are used to compute liquid and vapor properties respectively in the CO2 Removal Unit. The NRTL parameters for electrolyte pairs were regressed against VLE data from Maurer (1988) [1]. NH3 is selected as Henry-components to which Henrys law is applied and the Henrys constants are specified.

5 Chemical Reactions

Reforming Unit
The main reaction occurring in the reformer is the conversion of methane to a mixture of CO, CO2 and H2:

CH 4 + H 2 O CO + 3H 2

(Endothermic) (Exothermic)

(A-1) (A-2)

CO + H 2 O CO 2 + H 2

The kinetics expression for Reaction A-1, the reforming reaction, (Moe and Gerhard, 1965)[2] is as follows:
2 4 K 2 PCH 4 PH 2O PH 2 PCO2 lb mole R = k2 379 hr ft 3

(A-3)

The partial pressures are converted to expressions of molar quantities by assuming ideal gas behavior. The factor 379 is needed to convert the units from SCF to moles.
4 k2 P 3 2 ( H 2 ) ( CO 2 ) 2 K 2 ( CH 4 )( H 2O ) R= 3 P SS 379 SS 2

(A-4)

SS is the total number of moles of mixture per mole of methane fed. Moe and Gerhard arbitrarily set the P3 term to 1.0 in order to correlate their data taken from pressures above atmospheric. Therefore, the final form of the model does not have a P3 term, but is lumped in with the specific rate constant k2 as:

k 2 = Ac

31720 7.912 T + 460 e

(A-5)

The factor Ac, catalyst activity, is used in the model to give a reasonable reactor performance for the methane conversion. It may be adjusted such that reactor performance matches plant data. The equilibrium constant, K2, with units of atm2, is equal to K1K3. Hyman (1967) [3]gives:

 49, 435 K1 = exp + 30.707 atm 2 t + 460

K3 = exp K3

(A-6) (A-7) (A-8)

( = exp (

8,240 t + 460

4.335 below 1,100F 3.765 above 1,100F

7,351.24 t + 460

where, t is the temperature in F. Chemical equilibrium is assumed for the water shift reaction for which the following holds:

K 3=

(H 2 )(CO 2 ) (H 2O)(CO)

(A-9)

The reaction kinetics have been implemented in user kinetics Fortran subroutines of the RPlug model. The Fortran subroutines are REFKIN, DRATE, and KFORMC. Subroutines REFKIN and KFORMC are interface routines necessary to interface the Aspen Plus RPlug model with DRATE, the Fortran kinetics subroutine developed by Mok (1982)[4]. The pressure drop expression adopted is as follows:

dP = Pfact 0.04183 + 0.003292 Z - 0.0000395Z 2 dz

Where, Z: P: Pfact: or: tube length (ft) pressure (atm)

(A-10)

pressure drop factor (a function of catalyst characteristics)

P = Pfact 0.04183Z + 1 0.003292 Z 2 - 1 0.0000395Z 3 2 3

(A-11)

The pressure drop expression has been implemented as a user pressure drop Fortran subroutine of RPlug model. The Fortran subroutine is REFPD. Since the net reaction in the primary reformer is endothermic, the heattransfer rate is critical to the rate of reaction in the tube. The heat-transfer rate is determined by the tube inside and outside heat-transfer coefficients, the flame temperature, and the gas temperature. The heat transfer from the flame to the outside of the tube occurs primarily by radiation. Here an outside heat transfer coefficient is defined for ease of calculation.

Flux = htin (TW T ) = htout T f TW htoutT f + htinT htout + htin

(A-12)

The tube-wall temperature Tw is calculated from the second equality.

TW =

(A-13)

The flux is in units of BTU/hr/sq ft of inside tube area. The htin and htout are the inside and outside heat-transfer coefficients in units of BTU/hr/F/ft2 of inside tube-wall area. Both coefficients vary with gas temperature T. The inside heat-transfer coefficient is an overall physical coefficient determined only by the properties of the fluid. Beck (1962) gave the following correlation for estimating htin:

htin

DP kf

= 0.4 2.58 Re p

) ( Prp )
1/ 3

1/ 3

+ 0.094 Re p

) ( Prp )
0.8

0.8

(A-14)

Hyman (1967) reported that the factor 0.4 is needed to account for the smaller value of htin when ring-shaped catalyst is used instead of the pellets used by Beck. The calculation of the thermal conductivity kf, viscosity, and heat capacity were taken from Mok (1982). The heat-transfer rate outside the tube is assumed to be governed by a radiant heat-transfer mechanism:

Flux =

Dout Din

4 c T f4,r Tw,r = htout T f Tw

(A-15)

where Din and Dout are the inside and outside diameters of the tube, c is the product of the Stefan-Boltzman constant and the effective emissivity, Tf,r and Tw,r are the flame temperature and tube-wall temperature in degrees Rankine. The heat transfer rate was implemented as user heat transfer Fortran subroutines of RPlug model. The Fortran subroutines are REFHT, KFORMC, HTCOEF, TWALL, BPARM, and RKEQ. Subroutines REFHT and KFORMC are interface routines necessary to interface Aspen Plus RPlug model with other Fortran subroutines developed by Mok (1982). HTCOEF computes the heat transfer coefficient. TWALL computes the tube wall temperature. BPARM and RKEQ are physical property routines used to compute transport properties needed in HTCOEF and TWALL. The primary reformer is generally represented by 1-3 RPlug reactors in series, representing the number of firing zones of the unit in the plant. This allows the fuel requirements of each fired zone to be calculated separately. The fuel fed to the reformer contains hydrocarbons higher than methane, whereas Equation (A-1) only represents the reforming of methane. These higher hydrocarbons, however, are rapidly converted to methane and are thus accurately modeled using a simple RSTOIC block prior to the reformer tube with the following reactions:

3C 2 H 6 + H 2 O 5CH 4 + CO 3C 3 H 8 + 2H 2 O 7CH 4 + 2CO 3n - C 4 H10 + 3H 2 O 9CH 4 + 3CO 3i - C 4 H10 + 3H 2 O 9CH 4 + 3CO 3C 5 H12 + 4H 2 O 11CH 4 + 4CO 3C 6 H14 + 5H 2 O 13CH 4 + 5CO

(A-16) (A-17) (A-18) (A-19) (A-20) (A-21)

3C 7 H16 + 6H 2 O 15CH 4 + 6CO

(A-22)

The reformer burners are described by RSTOIC blocks in which the combustion is complete. The burner temperature sets the radiant heattransfer temperature for the reformer tubes; and the heat generated in the burner is equal to the heat absorbed in the reformer tubes. The combination of the RPlug blocks for the reformer tubes and RSTOIC blocks for the burners provides an accurate simulation of the reformer unit.

Carbon Monoxide Conversion

Low-Temperature Shift Reactor
The low-temperature shift reactor is modeled as a plug flow reactor, RPlug. The reaction stoichiometry considered is:

CO + H 2 O CO 2 + H 2
The kinetics expression (Slack, 1974)[5] is as follows:
1/ 2 KF Tref k LT YCOYH 2O 1- K 3 rco = Ac 1 T P + k AYCO + k BYCO2

(A-23)

(A-24)

KF =

( Y CO 2 )( Y H 2 ) ( Y CO )( Y H 2O )
= = = = = = = = catalyst activity exp (3620/T - 4.32126) standard LT catalyst activity in lb-mol/hr/ft3/atm exp (8240/(t+459.7) - 4.33) exp (4580/T - 7.4643) exp (1500/T - 2.623) atm-1 513.13K temperature in K temperature in F atm-1

(A-25)

where: Ac kLT K3 kA kB Tref T t

The reaction kinetics has been implemented in a user kinetics Fortran subroutine under the RPlug model. The Fortran rate subroutine is LTKIN, which calls KFORMC for obtaining the component locations.

High-Temperature Shift Reactor

The high-temperature shift reactor is modeled as a plug flow reactor, RPlug. The reaction stoichiometry considered is:

CO + H 2 O CO 2 + H 2
The kinetics expression (Slack, 1974) is as follows:

(A-26)

CO + H 2 O CO 2 + H 2

(A-27)
Kf K3

rCO = Ac K HT P1/ 2YCO 1

(A-28)

8240 - 4.33 below 1100 F K 3 = exp t + 459.7 KF= (YCO2 )(YH 2 ) (YCO )(YH 2O )
= = = = catalyst activity exp (10.3375 - 5787.62/T) standard HT catalyst activity in lb-mol/hr/ft3/atm1/2 temperature in K temperature in F

(A-29)

(A-30)

where: Ac kHT T t

The reaction kinetics was implemented in a user kinetics Fortran subroutine under the RPlug model. The Fortran rate subroutine is HTKIN, which calls KFORMC for obtaining the component locations.

Carbon Dioxide Removal

The electrolyte solution chemistry has been modeled with a CHEMISTRY model with ID of NH3. Chemical equilibrium is assumed with all the ionic reactions in the CHEMISTRY NH3. In addition, a kinetic REACTION model named NH3 has been created. In NH3, all reactions are assumed to be in chemical equilibrium except the reactions of CO2 with OH- and the reactions of CO2 with NH3.

A. Chemistry ID: NH3

1 2 3 4 5 Equilibrium Equilibrium Equilibrium Equilibrium Equilibrium

2H 2 O H 3 O + + OH
CO 2 + 2H 2 O H 3 O + + HCO 3 HCO 3 + H 2 O H 3 O + + CO 3 2

(A-31) (A-32) (A-33) (A-34) (A-35)

NH 3 + H 2 O NH + + OH 4
H 2 NCOO + H 2 O NH 3 + HCO 3

6 7

Salt Dissociation

+ NH 4 HCO 3 (S) NH 4 + HCO 3

(A-36) (A-37)

+ NH 4 HCO 3 NH 4 + HCO 3

B. Reaction ID: REA-NH3

1 2 3 4 5 6 7 8 EQUIL EQUIL EQUIL KINETIC KINETIC KINETIC KINETIC Salt
+ NH 3 + H 2 O NH 4 + OH

(A-38) (A-39) (A-40) (A-41) (A-42) (A-43) (A-44) (A-45)

2H 2 O H 3 O + + OH
HCO 3 + H 2 O H 3 O + + CO 3 2 CO 2 + OH HCO3

HCO3 CO 2 + OH

NH 3 + CO 2 + H 2O H 2 NCOO- + H 3O + H 2 NCOO- + H 3O + NH 3 + H 2O + CO 2
+ NH 4 HCO 3 (S) NH 4 + HCO 3

The equilibrium expressions for the reactions are taken from the work of Miles et al.[6]. In addition, the power law expressions are used for the ratecontrolled reactions (reactions 4-7 in REA-NH3) and the general power law expression is:

E 1 1 N ai n r = k (T T0 ) exp C i R T T0 i =1
If T0 is not specified, the reduced power law expression is used:

(A-46)

r = kT n exp (

E N ai ) C i RT i =1

(A-47)

In this file, the reduced expressions are used. In equation A-47, the concentration basis is Molarity, the factor n is zero, k and E are given in Table 2. The kinetic parameters for reactions 4-7 in Table 3 are derived from the work of Pinsent[7].

Table 1. Parameters k and E in Equation A-47

Reaction No. 4 5 6 7

k
4.32e+13 2.38e+17 1.35e+11 4.75e+20

E (cal/mol)
13249 29451 11585 16529

The built-in Keq expression (Equation 3) is used for the salt precipitation reaction of NH4HCO3,

ln K eq = A + B T + C ln (T ) + DT

(A-48)

The parameters A, B and C were regressed against SLE data from Trypuc et al. [8].

Table 3. Parameters A, B and C in Equation A-48

Reaction NH4HCO3 salt precipitation A -914.00821 B 38648.2117 C 136.174996

Methanation Unit
The methanation reactor, which removes trace amounts of CO from H2-rich synthesis feed mixtures, is modeled as a plug flow reactor by RPlug. Two reactions occur in the methanation reactor:

CO + 3H 2 CH 4 + H2O CO 2 + 4H 2 CH 4 + 2H 2 O

(A-49) (A-50)

Yadav and Rinker (1993)[9] developed empirical correlations for the methanation of CO (Equation A-49) as part of a laboratory data program. Their correlation, based upon a proposed Langmuir-Hinshelwood mechanism is as follows:

r=

0.5 APH 2 PCO 0.5 1 + BPCO + CPH 2

(A-51)

where r is the reaction rate in mol/g.s, P refers to partial pressures in kPa, and A, B and C are empirical constants determined at each temperature. Yadav and Rinker (1993) covered a broad range of CO concentrations, but for methanation we are mainly interested in very low CO levels. At low CO concentrations, Equation A-38 simplifies to the following form:

r=

A PCO
0.5 C 2 PH 2

(A-52)

The data reported by Yadav and Rinker at 503, 513 and 529 K were fitted to a simple Arrhenius equation as follows:

r = 0.314e

1 1 1,300 T 513

P CO 0.5 PH 2

(A-53)

Proper modeling of the methanation must include the effect of the backward reaction as equilibrium is approached. Hence, the complete version of Equation A-53 is as follows:

r = Ac 0.314e

1 1 1,300 T 513

P yCH 4 yH 2O 0.5 yCO 3 2 yH 2 P K CO PH 2

(A-54)

where, P is the pressure in kPa, y is the component mole fraction, KCO is the equilibrium constant for Reaction A-49 and Ac is the catalyst activity factor.

K CO = e (

38.4523 + 26,270 / T )

(A-55)

Kinetic data are not available for the methanation of CO2. We use the same kinetic expression for CO2 as for CO, recognizing the weakness of this approach, but also recognizing that the concentration of CO2 in the process stream is very small. The rate of Reaction A-50 is given as:

r = Ac 0.314e K CO2 = e (

1 1 1,300 T 513

2 P yCH 4 yH 2O 0.5 yCO2 4 2 yH 2 P K CO2 PH 2

(A-56)

33.923 + 21,621/ T )

(A-57)

Synthesis Unit
The ammonia synthesis converter beds are modeled as a plug flow reactor by RPlug. The reaction stoichiometry considered is:

N 2 +3H 2 2NH 3
The kinetics expression (Nielsen, 1968)[10] is as follows:

(A-58)

a2 2 AK aN keq a 3A H kgmole 3 R = Ac 2 m hr aA 1 + K a W aH
where: Ac w = = = catalyst activity .654 1.523

(A-59)

and where keq is the equilibrium constant, AK, the specific rate constant, and Ka, the adsorption equilibrium constant.

log10 keq = 2.691122 log10 (T ) 5.519265 105 T + 1.848863 107 T 2 + 2001.6 + 2.6899 T
(A-60)

AK = 3.945 1010 exp 5622 T K a = 2.94 104 exp

(A-61) (A-62)

( 12104 ) T

aN, aH, aA, in Equation (are the activities of nitrogen, hydrogen and ammonia. and w are parameters. Nielsen also gave the following formula for evaluating the activities:

P Ai Ci ai = xi P exp RT Bi RT 3 + Ai0.5 xi Ai0.5 T

where: R P T xi
i H2 N2 NH3 Ar CH4

(A-63)

= = = =
Ai

Gas constant (=0.0826) Pressure in atm Temperature in K mole fraction of component i

Bi .02096 .05046 .03415 Ci 5.04 x 102 4.20 x 104 4.77 x 106

The values of A's, B's and C's are listed below:

0.1975 1.3445 2.393 1.2907 2.2769

The reaction kinetics has been implemented by a user kinetics Fortran subroutine of RPlug model named SYNKIN. Subroutine SYNKIN calls RNH3, the Fortran kinetics subroutine developed by Mok (1981). Subroutine KFORMC is used to find the component integer variables from the Aspen Plus system.

6 Simulation Approaches

Aspen Plus Unit Operation Blocks Used in Ammonia production Model

Unit Operation Natural Gas Desulfurization Reforming Aspen Plus "Block" Sep2 + RStoic RStoic + RPlug Comments / Specifications The sulfur is hydrogenated to hydrogen sulfide and then removed The reactions are defined in the subroutines using RPlug. The reformed gas contains about 0.3 vol% CH4. CO converts to CO2 in two catalytic stages, the first at high temperature and the second at low temperature. CO2 is captured by NH3. The residual content of CO + CO2 is less than 10 ppm. Reacting with RPlug and high pressure. Ammonia gas in the synthesis loop is liquefied

Carbon Monoxide Conversion Carbon Dioxide Removal Methanation Synthesis Refrigeration

RPlug + FSplit + Heater

Flash2 RPlug + Sep2 RPlug + FSplit + Heater Valve + HeatX

Design Spec Used in Ammonia production Model

Spec Name F18 Spec (Target) Mass flow of stream S-18 is 44100 kg/hr, the tolerance is 0.1. Manipulated Variables The heat duty of Heater E003.

Calculator Used in Ammonia production Model

Name RESULTS Purpose Collect simulation results from the synthesis reactor.

Kinetics Model Implementation

A summary of the Fortran subroutines developed for the kinetics reactors is presented in Table A-1. Table A-2 shows the input parameters to the main subroutines.

Table A-1. Summary of Kinetics Fortran Subroutines

Primary Reformer High Temperature Shift HTKIN KFORMC Low Temperature Shift LTKIN KFORMC Methanator Synthesis Converter SYNKIN KFORMC RNH3

Kinetics Pressure Drop Heat Transfer

REFKIN KFORMC DRATE REFPD REFHT KFORMC HTCOEF TWALL BPARM RKEQ

MTHKIN

Table A-2. Parameter Input to Kinetics Fortran Subroutines

Main User Subroutine Name Parameter Default Value

REFKIN REFPD REFHT

R1 Catalyst activity R1 Pressure drop factor R1 Flame temperature R2 Maximum tube temperature R3 Radiative Heat transfer coefficient R1 - Catalyst activity R1 - Catalyst activity R1 - Catalyst activity R1 - Catalyst activity

1.0 1.0 975C 1200C 1.9*10-8 kcal/hr/m2 1.0 1.0 1.0 1.0

HTKIN LTKIN MTHKIN SYNKIN

7 Simulation Results

Key Stream Simulation Results

Flowsheet Variable Value Unit

Feed

Process Air Natural Gas Process Stream Fuel to Furnace 1 Fuel to Furnace 2 Combustion Air BFW

53000 32000 102041 10184 10854 178220 170180 3101

scmh scmh kg/hr scmh scmh scmh kg/hr kmol/hr

Product

NH3 (Liquid phase)

8 Conclusions

The Ammonia model provides a useful description of the process. The simulation has been developed using many of the capabilities of Aspen Plus including unit operation models, physical property methods, models and data, and flowsheeting capabilities like convergence design specs. The model provides rigorous mass and energy balance information for ammonia production, and can be used to support the conceptual process design. This model also can be used in developing reaction simulations with RPlug. The model is meant to be used as a guide for modeling ammonia production. Users may use it as a starting point for more sophisticated models.

9 References

[1] Goppert U., Maurer G., Vapor-Liquid Equilibrium in Aqueous Solutions of Ammonia and Carbon Dioxide at Temperatures Between 333 and 393K and Pressure up to 7 MPa, Fluid Phase Equilibrium, Vol. 41,153 (1988) [2] Moe, J.H. and Gerhard, E.R., Preprint 36d, 56th National Meeting, AIChE, May, 1965. [3] Hyman, M.H., Computational Study of a Refinery Hydrogen Plant, H.S. Thesis, Univ. of Calif., Berkeley, 1967. [4] Mok, L-F, Sensitivity Study of Energy Consumption in Ammonia Plant Operation, H.S. Thesis, Univ. of California, Berkeley, 1982. [5] Slack, A.V. and G.R. James, ed., Ammonia, Marcel Dekker, Inc., New York, 1974. [6] Miles D.H., Wilson G.M., "Vapor liquid Equilibrium Data for Design of Sour Water Strippers", Annual Report to the API for 1974, October 1975. [7] Pinsent B. R., Pearson L., Roughton F. J., The Kinetics of Combination of Carbon Dioxide with Ammonia, Trans. Faraday Soc., Vol. 52, 1594 (1956) [8] Trypuc M., Kielkowska U., "Solubility in the NH4HCO3 + NH4VO3 + H2O System", J. Chem. Eng. Data, Vol. 41, 1005(1996) [9] Yadav, R. and R.G. Rinker, Steady-State Methanation Kinetics Over a Ni/Al2O3 Catalyst, Can. J. of Chem. Eng., 71, 202 (1993). [10] Nielsen, A., An Investigation on Promoted Iron Catalysts for the Synthesis of Ammonia, 3rd ed., Jul. Gjellerups Forlag, 1968.