Dual Presure
Dual Presure
DUAL-PRESSURE PROCESS
EXAMPLE PURPOSE
This example corresponds to a dual-pressure nitric acid production process. It is a rather usual process of industrial
production of nitric acid. The main modules specific to the simulator ProSimPlus HNO3 are used in this simulation:
absorption column of nitrous vapors, nitrous vapors condenser, oxidation reactors, heat exchangers with oxidation
volumes, nitrous vapor compressors, etc.
Reader is reminded that this use case is only an example and should not be used for other purposes. Although this example is based on actual
case it may not be considered as typical nor are the data used always the most accurate available. ProSim shall have no responsibility or
liability for damages arising out of or related to the use of the results of calculations based on this example.
TABLE OF CONTENTS
This example is extracted from [BAD96] and [CLA96], which describe summarily the process. The production of
nitric acid includes three main steps:
Ammonia oxidation
4 +5 →4 +6
4 +3 →2 +6
2 +2 → +3
2 + →2 ↔
| ↔ |
3
+ ↔2 +
2
The mathematical representations of the equilibria and kinetics of the reactions described here above is complex.
Side reactions occur, in particular the formation of nitrogen trioxide:
+ ↔
The dual-pressure process (pressure of the absorption sensibly higher than the one of the catalytic oxidation) lies
on the reaction:
3
+ ↔2 +
2
The main characteristic of this reaction is that each time two molecules of nitric acid are formed there is also
production of one molecule of nitric oxide (NO). This molecule of NO has to be oxidized in NO2, and then absorbed
and so on. These successive oxidations are mainly done in gas phase in oxido-absorption tower where the
reactions of nitric acid production and NO oxidation are performed in parallel:
2 + →2
This process is mainly characterized by the use of two different pressures for the ammonia oxidation and the
absorption of the nitrogen oxides. This ensures a high yield in nitrogen and decreases the consumption of platinum.
The process flowsheet is given in paragraph 1.2. The liquid ammonia is vaporized (E101), then after overheating
(E102), is sent in an air – ammonia mixer (AIR NH3 MIXER). The atmospheric air is compressed
(AIR COMPRESSOR), then divided in two streams: the primary air (HP AIR 2) sent to the air – ammonia mixer
(AIR NH3 MIXER) and the secondary air sent to the bleaching column (BLEACHING COLUMN). The air –
ammonia mixture is sent to the burner (BURNER).
For the reactor modelling, three global reactions are taken onto account:
4 +5 →4 +6
4 +3 →2 +6
2 +2 → +3
The gas after combustion is composed of nitrogen oxides, nitrogen and oxygen. Its sensible heat is recovered in a
series of heat exchangers (E104, E105, E106, E107 and E108). After condensation (CONDENSER 1), a huge
quantity of weak nitric acid is formed and sent to the absorption tower (OXIDO-ABSORPTION COLUMN). The
second compressor (NOx COMPRESSOR) received the gas coming from the first condenser and from the
bleaching column. After cooling in the heat exchanger (E109), these gases enter in the sieve trays oxido-
absorption column (OXIDO-ABSORPTION COLUMN). Process water is introduced at the top of this column and
nitric acid at the whished concentration is withdrawn at its bottom. Then, this acid goes to the bleaching column
(BLEACHING COLUMN). A stripping by the secondary air (HP AIR 4) is done in this equipment. The gas leaving
the absorber C101 is sent to a serie of gas-gas heat exchangers (E103, E105, E107 and E109) and then in an
expander (EXPANDER). The gas leaves the process via a stack. The possible reactor for NOx reduction is not
modeled in this example. This aspect is detailed in the example “PSPH_E01_EN - Mono-pressure process”.
In parallel, steam is produced by energy integration. For that, the water is preheated (E106 and E108), then the
steam is produced in the boiler (BOILER) and overheated (OVERHEATER). A part of this steam is expanded
(STEAM EXPANDER) to bring the mechanical energy needed by the process. This steam is then condensed
(STEAM CONDENSER). From a practical point of view, the two compressors (AIR COMPRESSOR and
NOx COMPRESSOR) are on the same shaft than the two expanders (EXPANDER and STEAM EXPANDER).
Information streams and information stream handler modules are used to model this coupling.
The objective of this process is to produce 1 000 t/d of nitric acid (eq. 100%) at a concentration of 58% mass. At
the steam production level, this plant produces 15 bars steam by energy integration.
The oxygen amount in the tail gas is fixed at 2.5% volume (mol.). The temperature at the outlet of the burner is
fixed to 890°C. The water flow rate for absorption, the ammonia feed flow rate, the air feed flow rate and the ratio
between the primary and the secondary air are automatically adjusted to ensure this production. The amount of
steam expanded is adjusted to balance the power available at the turbines and those necessary for the
compressors.
TG 6
Dual-pressure process
d(TBurn …
Water
E104 FLARE
D(Water)
WATER FEED
E103a
PG 02
D(Air)
OXIDO-ABSORPTION COL…
Process flowsheet
E106a
E105a
PIPE 4 PG 11
E109a
PG 07
NOx COMPRESS… PG 13 PG 14
Q(Air compressor)
PG 15
d(DHN…
ACID 3 Th ACID 4
CONDENSE … ACID 1
Q(NOx compressor)
ACID MIX… Calculator Switch
d(wHN…
PUMP
BLEACHING COLU…
ACID 5
Dual-pressure HP AIR 4
AAA AAA
STORAGE 1 1
Q(E1… Q(E1… Q(E104)
COMPRESSO… w(HN … D(HN…
STEAM 4 C ONDENSAT …
STEAM EXPAND…
Q(Boil… Q(Overhea…
Q(Steam expand…
STEAM 3
EXPANDE…
SPLIT RA…
Q(Expande…
SPEC DEVIATI …
WU 1
STEAM SPLITT…
WU 2 WU 3 STEAM 1 STEAM 2
www.prosim.net
Dual-pressure process
1.3. Components
Components taken into account in the simulation, their chemical formula and CAS numbers are presented in the
following table. Pure components physical properties are extracted from the ProSimPlus HNO3 specific database
(“HNO3”).
For the main part of the process, the “HNO3 specific” thermodynamic model is selected. This model takes into
account the non-ideality of the liquid phase through correlations based on experimental data of partial pressures of
water and nitric acid over aqueous solutions of nitric acid. The perfect gas model is used for the gas phase. A
correlation based on experimental data is used to take into account the excess enthalpies of the water – nitric acid
binary.
For the bleaching column the “Engels (strong acids)” thermodynamic model is used. This model well represents the
NOx solubility and the complexity of the H2O – HNO3 equilibria. The “H*=DH0f, ideal gas, 25°C, 1 atm” enthalpy
basis is selected instead of the default one of this model, to ensure the coherence with the “HNO3 specific” model.
For the cooling water circuit, the “Specific thermodynamic model for water” is used [HAA84].
NO oxidation, NO2 dimerization and nitric acid oxido-absorption reactions are pre-coded. Thus, they don’t have to
be described by the user for modules specific to ProSimPlus HNO3. The reactions to be described are the ones for
the combustion of the ammonia:
4 +5 →4 +6
4 +3 →2 +6
2 +2 → +3
The instantaneous kinetic model is chosen because conversion ratios are specified for each reaction in the
reactors “BURNER” (§ 1.6.1).
The NO2 dimerization equilibrium in gas phase is taken into account in the bleaching column. Indeed this column is
modelled by a non-specific module of ProSimPlus HNO3 (§ 1.6.1).
2 ↔
6891.64
( ) = −21.24366 +
The partial orders are 2 for the NO2 and 1 for the N2O4.
The heats of reaction of these different reactions are computed from the standard enthalpies of formation at 25°C.
Process feeds
Oxidation reactors
These modules model the oxidation of the NO and the dimerization of the NO2 in the gas volumes of pipes.
Cooler/Heater modules
Name Outlet temperature (°C) Pressure drop (bar) Oxidation volume (m3)
E101 Dew temperature 8,00 0,0
E102 080 0.05 0,0
E103a 070 0.05 0,0
E104 440 0.05 0,0
E105a 290 0.05 1.8
E106a 220 0.05 1.8
E107a 180 0.05 1.8
E108a 110 0.05 1.8
E109a 125 0.05 1.8
The following parameters are used for the heat exchangers in which the chemical reactions are taken into account
(the ones with an oxidation volume different from 0):
Name Heat duty (kW) Pressure drop (bar) Oxidation volume (m3)
E103b 0 0.05 0
E105b 0 0.05 0
E107b 0 0.05 0
E109b 0 0.05 0
The heat duties are sent by information streams from the corresponding cooler/heater modules (E103a, E105a,
E107a and E109a).
CONDENSER 1 CONDENSER 2
Tubes length (m) 6
Number of tubes 840 250
Circulation of the vapors Inside the tubes
Inner diameter of the tubes (mm) 25.4
Equivalent diameter (mm) 0
Cooling water
Temperature (°C) 20 15
Flow rate (kg/h) 375 000 271 000
Flow direction Counter-current
Heat transfer coefficients (kcal/h/m2/K)
Oxido-absorption 300 600
Oxidation 0
Pressure drop (bar) 0.05
Reactions
Calculation of the oxidation rate constant Koukolik
Calculation of the dimerization equilibrium constant Koukolik
Calculation of the absorption constant of N2O4 in water Miller (bubble caps)
Number of intermediate points for print 10
Compressors
Mixers
The default parameters are used for the two mixers “AIR NH3 MIXER” and “ACID MIXER3.
“OXIDO-ABSORPTION COLUMN”
The column stages are numbered from top to bottom (plate 1: top plate; plate 30: bottom plate).
Calculator switch
This module is used to change the thermodynamic model from the “HNO3 specific” one used in the previous
modules to the “Engels (strong acids)” used in the bleaching column. The calculation of the liquid enthalpies being
not the same between these two models, the use of a “Calculator switch” module is necessary in order to not
distort the enthalpy balance on the bleaching column. In the case of the switch from the “HNO3 specific” model to
the “Engels (strong acids)” model, and inversely, the “Calculator switch” module is necessary only on liquid streams
(except if they leave the process).
“BLEACHING COLUMN”
Specifications
Specifications Value
Outlet temperature of the ammonia combustion burner (°C) 890
Oxygen amount in the tail gas (% molar) 2.5
Concentration of the produced nitric acid (% mass) 58
Mass flow rate (eq. 100%) of the produced nitric acid (t/j) 1 000
Acting variables
Acting variables
Ammonia feed flow rate
Air feed flow rate
Water feed flow rate of the oxido-absorption column
Splitting ratio between the primary air and the secondary air
ProSimPlus HNO3 determines automatically the calculation sequence and the tear streams. The tear
stream selected by default (“PG 12” outlet stream of the NOx compressor) needs an initialization to reach
the convergence of the simulation. To avoid this, another stream, “PG 15”, is selected directly in the
module “Constraints and Recycles”. Thus, this choice avoids the need to specify an initialization of the
tear stream.
Process feed
“STEAM EXPANDER”
Cooler/Heater modules
Name Outlet temperature (°C) Pressure drop (bar) Oxidation volume (m3)
BOILER Dew temperature 0.05 0
STEAM CONDENSER Bubble temperature 0.05 0
The heat duties of the simple heat exchangers E106b and E108b are transmitted by information streams from the
corresponding cooler/heater modules (E106a and E108a).
The heat duty of the simple heat exchanger “OVERHEATER” is the sum of the heat duties calculated by the
cooler/heater modules “E104” and “BOILER”. The information is transmitted by information stream and the sum is
computed by an “Information stream handler” module (§ 1.7). The parameters of these different items are the
following:
o Information stream from the cooler/heater “E104” to the information stream handler “BOILER HEAT DUTY”
Parameters Value
Information type to be emitted Heat necessary to reach the specified temperature
Information type to be received Input information stream value (In)
o Information stream from the cooler/heater “BOILER” to the information stream handler “BOILER HEAT
DUTY”
Parameters Value
Information type to be emitted Heat necessary to reach the specified temperature
Information type to be received Value of the additive factor (B)
o Information stream from the information stream handler “BOILER HEAT DUTY” to the simple heat
exchanger “OVERHEATER”
Parameters Value
Information type to be emitted Output information stream value (Out)
Information type to be received Heat duty
Parameters Value
A 1
B 0
C 0
Power 1 (Integer value)
In this process, the compressors “AIR COMPRESSOR” and “NOx COMPRESSOR” and the expanders
“EXPANDER” and “STEAM EXPANDER” compose a turbo-expander. The objective is to adjust the steam flow rate
expanded by action on the splitting ratio of the “STEAM SPLITTER” in order to balance the power consumption of
the compressors (“AIR COMPRESSOR” and “NOx COMPRESSOR”) and the power available at the expanders
(“EXPANDER” and “STEAM EXPANDER”). For that, several information stream handlers (§ 1.7) as well as a
constraints and recycles module are used. Several information streams link these different items.
Sum of the power of the compressors “AIR COMPRESSOR” and “NOx COMPRESSOR”
o Information stream from the compressor “NOx COMPRESSOR” to the information stream handler
“COMPRESSORS”
Parameters Value
Information type to be emitted Mechanical power
Information type to be received Input information stream value (In)
o Information stream from the compressor “AIR COMPRESSOR” to the information stream handler
“COMPRESSORS”
Parameters Value
Information type to be emitted Mechanical power
Information type to be received Value of the additive factor (B)
Parameters Value
A 1
B 0
C 0
Power 1 (Integer value)
Parameters Value
Information type to be emitted Mechanical power
Information type to be received Input information stream value (In)
o Information stream from the expander “STEAM EXPANDER” to the information stream handler
“EXPANDERS”
Parameters Value
Information type to be emitted Mechanical power
Information type to be received Value of the additive factor (B)
Parameters Value
A 1
B 0
C 0
Power 1 (Integer value)
Calculation of the deviation between the power needed for the compressors and the power available at the
turbines and action on the splitting ratio of the “STEAM SPLITTER” splitter
o Information stream from the information stream handler “COMPRESSORS” to the information
stream handler “TURBO-EXPANDER”
Parameters Value
Information type to be emitted Output information stream value (Out)
Information type to be received Input information stream value (In)
o Information stream from the information stream handler “EXPANDERS” to the information stream
handler “TURBO-EXPANDER”
Parameters Value
Information type to be emitted Mechanical power
Information type to be received Value of the additive factor (B)
The two values are summed although we wish to compute the deviation between the power needed for
the compressors and the power available at the turbines. In fact, the power needed for a compressor are
negative values, whereas the power available for a turbine are positive values.
Parameters Value
A 1
B 0
C 0
Power 1 (Integer value)
o Information stream from the information stream handler “TURBO-EXPANDER” to the “Constraints
and Recycles 1” constraints and recycles module
Parameters Value
Information type to be emitted Output information stream value (Out)
Information type to be received Automatic
o The default values are kept for the module “Constraints and Recycles 1”
o Information stream from the module “Constraints and Recycles 1” to the splitter “STEAM
SPLITTER”
Parameters Value
Information type to be emitted Automatic
Information type to be received Distribution ratios of output streams
Regarding the information type to be received, well select the stream for which the splitting ratio is
specified and not the stream automatically calculated. The choice is done from the bottom red framed of
the following figure.
To minimize the risk of errors and to correct more easily those eventually committed, it’s advisable to proceed in
four steps to build this example:
1. Simulation of the nitric acid production part without the specifications, the utility production and the coupling
between the compressors and the turbines.
2. Add of the specifications to the previous simulation.
3. Add of the utility production part to the previous simulation.
4. Add of the coupling between the compressors and the turbines to the previous simulation.
In this simulation example, most of the heat exchangers are simulated by decoupling the hot stream and the cold
stream. This way of modeling a two-stream heat exchanger avoids a recycle stream that would penalize the
calculations by decoupling the heat exchanger into two parts.
Several “Information stream handler” modules are used. In ProSimPlus HNO3 this module performs simple
operation on information streams: add a constant (B), subtract a constant (C), multiply its contents (A and P = 1),
divide its contents (A and P = -1) or raised its contents to an entire or real power (P). The result of this operation is
available in the outlet information stream of the “Information stream handler” module:
= ∗( ) + −
In this formula, A, B, C, and P are constants specified by the user. By putting several information streams in an
"Information stream handler" module and playing on their position in the module parameter zone, it’s possible to
use their information to specify the values of the constants A, B, C and P. Thus, operations between several
information streams can be performed. Using several modules of this type, it’s possible to perform more complex
operations. For more complex operation on information streams it’s advised to use a “Windows Script” module.
It’s possible in ProSimPlus HNO3 to use several constraints and recycles modules in one simulation. In this
example for a better readability of the simulation, two constraints and recycles modules are used:
Using several constraints and recycles modules has an effect on the calculation sequence of the module (the list of
the calculations automatically determined by ProSimPlus HNO3) and on the convergence of each cycle. But this
has no effect on the results obtained at convergence.
This possibility could be used in some case for ease of the convergence of complex flowsheet, or, like here, to not
unnecessarily complicate the representation of the process.
2. RESULTS
This document presents only the most relevant stream results. In ProSimPlus HNO3, mass and energy balances
are provided for all streams. Results are also available at the unit operation level (“Report” tab in the configuration
window).
The production of 1 000 t/d of nitric acid eq. 100% at a concentration of 58% mass with this process while having a
temperature of 890°C at the outlet of the ammonia burner and 2.5% molar of oxygen in the tail gas needs:
The NOx amount at the outlet of the oxido-absorption column is 1 070 ppmv eq. NO.
The balance between the consumption of compressors and the power available at the turbines needs to expand
86% of the 15 bar steam produced (966 t/d).
The following curve shows that 15% of the length of the “CONDENSER 1” are not used for the oxido-absorption.
The following curve shows that 33% of the length of the “CONDENSER 2” are not used for the oxido-absorption.
The column stages are numbered from top to bottom (plate 1: top plate; plate 30: bottom plate).
The column stages are numbered from top to bottom (plate 1: top plate; plate 5: bottom plate).
3. REFERENCES
[BAD96] BADOUAL C., “Acide nitrique”, Techniques de l'Ingénieur, J6085 (1996)
[CLA96] CLARKE S.I., MAZZAFRA W.J., “Nitric Acid”, Kirk-Othmer Encyclopedia of Chemical Technology,
4th edition, 17, 80-107 (1996)
[HAA84] HAAR L., GALLAGHER J.S., KELL J. H., “NBS/NRC Steam Tables”, Hemisphere Publishing
Corporation (1984)
[JOU81] JOULIA X., “Contribution au développement d'un programme général de simulation. Application à
l'analyse du fonctionnement d'un atelier de production d'acide nitrique”, Thèse INPT (1981)
[KOU68] KOUKOLIK M., MAREK J., “Mathematical Model of HNO3 Oxidation-Absorption Equipment”, Proc.
Fourth European Symp. on Chem. Reaction Eng. (suppl. Chem. Eng. Sci.), 347-359 (1968)