Pse Aspen
Pse Aspen
Speaker: 曾安紘
[email protected]
PSE Laboratory
Department of Chemical Engineering
National Taiwan University
(化工館RM309, 02-33663068)
Week 1~2
Outline
• Part 1 : Introduction
• Part 2 : Startup
• Part 3 : Properties analysis
3-1: Basic Property analysis
3-2: System without built-in thermodynamic parameters
3-3: System with User Defined Component
3-4: Validation of thermodynamic parameters
• Part 4 : Running Simulation in Aspen Plus
(simple units: Mixer, Pump, valve, flash, heat exchanger )
2
Introduction to Aspen Plus
Part 1: Introduction
3
What is Aspen Plus
4
Ref: http://www.aspentech.com/products/aspen-plus.cfm
What Aspen Plus provides
5
Ref: Aspen Plus® Product Brochure
More Detailed
• Properties analysis
– Properties of pure component and mixtures (Enthalpy, density,
viscosity, heat capacity,…etc)
– Phase equilibrium (VLE, VLLE, azeotrope calculation…etc)
– Parameters estimation for properties models (UNIFAC method for
binary parameters, Joback method for boiling points…etc)
– Data regression from experimental deta
• Process simulation
– pump, compressor, valve, tank, heat exchanger, CSTR, PFR, distillation
column, extraction column, absorber, filter, crystallizer…etc
6
What course Aspen Plus can be employed for
7
Lesson Objectives
8
Introduction to Aspen Plus
9
Start with Aspen Plus
11
Properties
Help
12
More Information
Help for
Commands for Controlling Simulations
13
Properties -- Setup
14
Properties -- Input components
15
Properties -- Methods
16
Property Method Selection—General Rule
Example 1:
water - benzene
Example 2:
benzene - toluene
17
Typical Activity Coefficient Models
Non-Randon-Two Liquid Model (NRTL)
Uniquac Model
Unifac Model
Typical Equation of States
NRTL
20
NRTL – Binary Parameters
21
Access Properties Models and Parameters
22
Review Databank Data
Including:
Ideal gas heat of formation at 298.15 K
Ideal gas Gibbs free energy of formation at
298.15 K
Heat of vaporization at TB
Normal boiling point
Standard liquid volume at 60°F
….
23
Pure Component
Temperature-Dependent Properties
24
Properties – Tools
27
Properties – Data Source
28
Properties – Run Mode
29
Properties – Analysis
30
Simulation
Help
Model Palette
(View│Model Palette)
31
Simulation -- Setup – Specification
Input mode
32
Basic Input---Summary
33
Introduction to Aspen Plus
34
Overview of Property Analysis
Use this form To generate
Tables and plots of pure component properties as a function of temperature
Pure
and pressure
Binary Txy, Pxy, or Gibbs energy of mixing curves for a binary system
Residue Residue curve maps
Ternary maps showing phase envelope, tie lines, and azeotropes of ternary
Ternary
systems
This feature locates all the azeotropes that exist among a specified set of
Azeotrope
components.
Ternary diagrams in Aspen Distillation Synthesis feature: Azeotropes,
Ternary Maps Distillation boundary, Residue curves or distillation curves, Isovolatility curves,
Tie lines, Vapor curve, Boiling point
35
Find Components (1/5)
36
Find Components (2/5)
37
Find Components (3/5)
38
Find Components (4/5)
39
Find Components (5/5)
40
Methods (1/2) –
Select Thermodynamic Model
Select NRTL
41
Methods (2/2) –
Check Binary Parameter for V8.4-V9
Properties
Methods
Parameters
NRTL-1
42
Methods (2/2) –
Check Binary Parameter for V10
Properties Analysis – Pure Component (1/5)
44
Available Properties (2/5)
5. Specify pressure
Add “N-butyl-acetate”
46
Example1: Calculation Results of CP (4/5)
47
Example1: Calculation Results of CP (5/5)
Data results
48
Properties Analysis – Binary Components (1/7)
Binary Component Properties Analysis (2/7)
3. Select
compositions basis
1. Select
analysis type (Txy)
2. Select
two component
6. Specify pressure
3. Select
compositions basis
4. Specify
composition range
Data results
55
Example: Generate XY plot (6/7)
Click “Input”
Example: Generate XY plot (7/7)
Properties Analysis – Ternary (1/6)
(add one new components)
Properties Analysis – Ternary (2/6)
(Check NRTL binary parameter)
1. Select
three component
2. Select
property method
5. Specify pressure
3. Select
phase (VLE, LLE…)
4. Specify
number of tie line 6. Specify temperature
(if LLE is slected)
Data results
63
Property Analysis – Conceptual Design
64
Conceptual Design -- Azeotrope Analysis (1/3)
Conceptual Design -- Azeotrope Analysis (2/3)
1. Select components (at least two) 2. Specify pressure
71
A 3-component system
Run Estimation
Property Estimation Mode
81
User Defined components
• For those components which are not included in
Aspen Plus databank.
• Information needed:
--Molecular Structure (Must)
--Molecular Weight (optional)
--Saturation vapor pressure (optional)
--Boiling point, Melting point…..etc (optional)
Example: Building Trisilylamine
Find
Trisilylamine
Example: Building Trisilylamine
No Trisilylamine
in the databanks
Example: Building Trisilylamine
• Check NIST chemistry webbook for those information
you need to build a user defined component.
(http://webbook.nist.gov/chemistry/)
Searching
Method
Example: Building Trisilylamine
86
Example: Building Trisilylamine
87
Example: Building Trisilylamine
88
Example: Building Trisilylamine
Example: Building Trisilylamine
90
Example: Building Trisilylamine
91
Example: Building Trisilylamine
92
Example: Building Trisilylamine
93
Example: Building Trisilylamine
95
Example: Building Trisilylamine
96
Example: Building Trisilylamine
97
Example: Building Trisilylamine
98
Example: Building Trisilylamine
99
Hand-On Exercises
• Which kind of model should be used?
(EOS or activity coefficient model?)
• Draw the Txy diagram for these two components under 1 atm.
• Draw the X-Y diagram for these two components under 1 atm.
Introduction to Aspen Plus
101
Finding Experimental Data in Aspen Plus
<Goal>
Compare the experimental VLE data and
the predicted VLE data by NRTL model.
Comparing the VLE with Built-in Parameter
bp: -24.8 oC
bp: 100 oC
118
Find NIST Data
119
Draw the Figure of Experimental Figure
120
Check Binary Interaction Parameters
121
Use Binary to Draw the Figure from Model
122
Combine the Two Figures
123
Hand-On Exercises
• Find the Binary VLE Data of Ethanol and Methanol at
T=313.15, 333.15, and 353.15 (K).
Simple Units
(Mixer, Pump, valve, flash, heat exchanger)
125
Example 1: Calculate the mixing
properties of two stream
1 4
3
2
Mixer Pump
1 2 3 4
Mole Flow kmol/hr
WATER 10 0 ? ?
BUOH 0 9 ? ?
BUAC 0 6 ? ?
Total Flow kmol/hr 10 15 ? ?
Temperature C 50 80 ? ?
Pressure bar 1 1 1 10
Enthalpy kcal/mol ? ? ? ?
Entropy cal/mol-K ? ? ? ?
Density kg/cum ? ? ? ?
126
Simulation -- Setup – Specification
Select Steady-State
Select Simulation
127
Reveal Model Library
128
Adding a Mixer
129
Adding Material Streams
130
Adding Material Streams (cont’d)
131
Adding Material Streams (cont’d)
132
Adding Material Streams (cont’d)
133
Reconnecting Material Streams (Feed Stream)
134
Reconnecting Material Streams
(Product Stream)
135
Change Material Streams Names
136
Specifying Feed Condition
137
Specifying Feed Condition (cont’d)
1 2
138
Specifying Input of Mixer
139
Specifying Input of Mixer (cont’d)
140
Run Simulation
141
Status of Simulation Results
Message Means
The run has completed normally, and results are
Results Available
present.
Results for the run are present. Warning
messages were generated during the
Results with Warnings
calculations. View the Control Panel or History
for messages.
Results for the run are present. Error messages
Results with Errors were generated during the calculations. View the
Control Panel or History for messages.
Results for the run are present, but you have
changed the input since the results were
Input Changed
generated. The results may be inconsistent with
the current input.
142
Stream Results
143
Heat and Material Balance Table
Stream ID 1 2 3
From B1
To B1 B1
Phase LIQUID LIQUID LIQUID
Substream: MIXED
Mole Flow kmol/hr
WATER 10.00000 0.0 10.00000
BUOH 0.0 9.000000 9.000000
BUAC 0.0 6.000000 6.000000
Total Flow kmol/hr 10.00000 15.00000 25.00000
Total Flow kg/hr 180.1528 1364.066 1544.218
Total Flow l/min 3.096993 29.00351 31.22363
Temperature C 50.00000 80.00000 71.21621
Pressure bar 1.000000 1.000000 1.000000
Vapor Frac 0.0 0.0 0.0
Liquid Frac 1.000000 1.000000 1.000000
Pull down the list and select Solid Frac 0.0 0.0 0.0
Enthalpy cal/mol -67813.56 -94336.22 -83727.15
“Full” to show more properties Enthalpy cal/gm -3764.224 -1037.372 -1355.494
results. Enthalpy Gcal/hr -.6781356 -1.415043 -2.093179
Entropy cal/mol-K -37.52634 -134.8439 -95.61121
Entropy cal/gm-K -2.083028 -1.482816 -1.547890
Density mol/cc .0538156 8.61965E-3 .0133445
Density gm/cc .9695038 .7838510 .8242788
Average MW 18.01528 90.93771 61.76874
Enthalpy and Entropy Liq Vol 60F l/min 3.008333 26.95642 29.96476
144
Change Units of Calculation Results
145
Setup – Report Options
146
Stream Results with Format of Mole Fraction
147
Add Pump Block
148
Add A Material Stream
149
Connect Streams
150
Pump – Specification
3. Efficiencies (Default: 1)
151
Run Simulation
152
Block Results (Pump)
153
154
Streams Results
155
Calculation Results
(Mass and Energy Balances)
1 4
3
2
Mixer Pump
1 2 3 4
Mole Flow kmol/hr
WATER 10 0 10 10
BUOH 0 9 9 9
BUAC 0 6 6 6
Total Flow kmol/hr 10 15 25 25
Temperature C 50 80 71.22 72.33
Pressure bar 1 1 1 10
Enthalpy kcal/mol -67.81 -94.34 -83.73 -83.67
Entropy cal/mol-K -37.53 -134.84 -95.61 -95.45
Density kg/cum 969.50 783.85 824.28 823.01
156
Exercise
1 5 6
2 4
3 Mixer Pump
1 2 3 4 5 6
Mole Flow kmol/hr
Water 10 0 0 ? ? ?
Ethanol 0 5 0 ? ? ?
Methanol 0 0 15 ? ? ?
Total Flow kmol/hr 10 15 15 ? ? ?
Temperature C 50 70 40 ? ? ?
Pressure bar 1 1 1 1 4 2
Enthalpy kcal/mol ? ? ? ? ? ?
Entropy cal/mol-K ? ? ? ? ? ?
Density kmol/cum ? ? ? ? ? ?
120
T-x
T-y
115
T ( C)
T=105 C 110
o
Saturated Feed P=1atm
105
P=1atm
F=100 kmol/hr 100
zwater=0.5 0.0 0.2 0.4 0.6 0.8 1.0
zHAc=0.5 xWater and yWater
6. Specify pressure
2. Select two
component
4. Specify
composition range
8. click “Run Analysis” to generate the results
Calculation Result of T-xy
Calculation Result of T-xy
Data results
167
Generate xy plot
Generate xy plot (cont’d)
Flash Separation
120
T-x
T-y
115
T ( C)
T=105 C 110
o
Saturated Feed P=1atm
105
P=1atm
F=100 kmol/hr 100
zwater=0.5 0.0 0.2 0.4 0.6 0.8 1.0
zHAc=0.5 xWater and yWater
Saturated Feed
(Vapor fraction=0)
P=1atm
F=100 kmol/hr
zwater=0.5
zHAc=0.5
Block Input: Flash2
Flash2: Specification
T=105 C
P=1atm
Required Input Complete
177
Stream Results
Stream Results (cont’d)
42.657 kmol/hr
zwater=0.591
zHAc=0.409
T=105 C
Saturated Feed P=1atm
P=1atm
F=100 kmol/hr
zwater=0.5 57.343 kmol/hr
zHAc=0.5 zwater=0.432
zHAc=0.568
EXAMPLE 3: HEAT EXCHANGE (SIMPLE HEATER)
Q: Calculate the enthalpy required to heat the mixture
from 5.2 atm, 50 ⁰C to 5 atm 250 ⁰C using Peng-
Rob model.
180
HEAT EXCHANGE (SIMPLE HEATER)
Using Aspen Plus…
181
HEAT EXCHANGE (SIMPLE HEATER)
182
HEAT EXCHANGE (SIMPLE HEATER)
183
HEAT EXCHANGE (SIMPLE HEATER)
184
HEAT EXCHANGE (SIMPLE HEATER)
185
HEAT EXCHANGE (SIMPLE HEATER)
186
HEAT EXCHANGE (SIMPLE HEATER)
3. Run
187
HEAT EXCHANGE (SIMPLE HEATER)
188
HEAT EXCHANGE (SIMPLE HEATER)
189
HEAT EXCHANGE (SIMPLE HEATER)
190
HEAT EXCHANGE (SIMPLE HEATER)
191
CALCULATE UTILITY (SIMPLE HEATER)
192
CALCULATE UTILITY (SIMPLE HEATER)
193
CALCULATE UTILITY (SIMPLE HEATER)
194
CALCULATE UTILITY (SIMPLE HEATER)
195
EXAMPLE 4: HEAT EXCHANGE (HEATX)
Hot stream
Inlet T:200℃、Inlet P:0.4 MPa
Flow rate:10000kg/hr
composition:benzene 50%, styrene
20%,EB 20%,water 10%
Cold stream
Inlet T:20℃、inlet P:1.0 MPa
Flow rate:60000 kg/hr。
197
Thermodynamic Model – NRTL
198
Thermodynamic Model – NRTL
199
ADD BLOCK: HEATX
Feeds Conditions
CLD-IN
Feeds Conditions
HOT-IN
202
BLOCK INPUT
Check result
204
Check result
205
EXAMPLE 5: USING HEATER FOR 2 STREAMS
HEAT EXCHANGE
207
Determine Pressure Drop
208
Connecting Heaters
209
Results of Heat Exchange
• Heat exchange amount will be exactly the same, while the positive and
negative value correspond to endothermic and exothermic process.
• No automatic examination for temperature crossing problem using this
technique, so please be careful.
210
EXAMPLE 6: MULTI COMPRESSOR
(MCOMPR)
33 oC 80 oC
1 bar 120 bar
212
MCOMPR - FLOWSHEET
213
SETTINGS FOR MCOMPR
Totally 3 compressors
214
SETTINGS FOR MCOMPR
215
SETTINGS FOR MCOMPR
218
RESULTS OF MCOMPR
219
EXAMPLE 7: MULTI HX (MHEATX)
221
Thermodynamic Model – PR
222
Thermodynamic Model – PR
223
ADD BLOCK: MHEATX
Feeds Conditions
HOT-IN
225
Feeds Conditions
COLDIN1
COLDIN2
BLOCK INPUT
Temperature
approach should be
Extra zone because larger than 5 oC
of different cold
stream temperature
230
Thanks for your attention!
PSE Laboratory
Department of Chemical Engineering
National Taiwan University
(化工館 RM307)