Process Simulation First Steps
Process Simulation First Steps
Maurizio Fermeglia
Maurizio.fermeglia@units.it
Department of Engineering & Architecture
University of Trieste
Agenda
Process simulation goals and definitions
Material and energy balances
Degrees of freedom analysis
Process simulation fundamentals,
Process simulation in the life cycle of chemical processes
Trends of process simulation
Users guidelines
Limitations and points of attention in process simulation
Conclusions and perspectives
System Graphics
Manager Reports
Equipment sizing
Utilities Forms
Cost estimation
Database
(Physical
Economic and Properties,
Environmental Costs, TD
Evaluation models, etc)
Splitter
Solve
model
equations
Solve model
equations –
AE solver
Solve model
equations –
DAE solver
Splitter
Stream
data
Operating
data
Splitter
Unit
models
Method of
solution
IGCC Plant
IGCC Flowsheet
Recycle with decomposition
stream
Tear
stream
Process Simulation – Maurizio Fermeglia Trieste, 28 September, 2021 - slide 19
Material and energy Balances
The general balance equation
𝑁𝑖 = 𝑁𝑣 − 𝑁𝑟
𝑁𝑖 = independent variables
Inherent relations (specific of the system)
𝑁𝑣 = total variables Mass balance relations
𝑁𝑟 = independent equations Energy balance relations
Phase equilibria relations
Chemical equilibria relations
𝑵_𝒊>𝟎: "problem is underspecified and additional variables must be specified in order to determine
the remaining variables"
𝑵_𝒊=𝟎: problem can be solved
𝑵_𝒊<𝟎: problem is overdetermined with redundant and possibly inconsistent relations
Process Simulation – Maurizio Fermeglia Trieste, 28 September, 2021 - slide 21
Stream
𝑁𝑣 = 𝑇, 𝑃, 𝑚,ሶ 𝐻𝑠𝑡𝑟𝑒𝑎𝑚 , 𝑥𝑖 = 1 + 1 + 1 + 1 + 𝑐 − 1 = 𝑐 + 3
𝑁𝑟 = 1
Enthalpy is a function of
composition, T and P
or
𝑁𝑣 = 𝑇, 𝑃, 𝑚,ሶ 𝐻𝑠𝑡𝑟𝑒𝑎𝑚 , 𝑥𝑖 = 1 + 1 + 1 + 1 + 𝑐 = 𝑐 + 4
𝑁𝑟 = 1 + 1
𝑐
𝑥𝑖 = 1
𝑖 𝑁𝑖 = 𝑐 + 2
Enthalpy is a function of
composition, T and P e.g.
c-1: components
1: temperature
1: pressure
1: mass flow
𝑁𝑟 = 2 + 2 ∗ 𝑐 − 1 + 1 + 1 = 2𝑐 + 2
Constant
Inherent relations
composition among
(constant T and P Mass Balance Energy Balance
inlet and 2 outlet
for inlet and outlets)
streams
𝑁𝑖 = 𝑐 + 5
e.g.
c+2: definition of inlet stream
1: temperature
1: pressure
1: any composition or heat
Process Simulation – Maurizio Fermeglia Trieste, 28 September, 2021 - slide 23
Mixer
𝑁𝑣 = 3 ∗ 𝑐 + 2 + 1 = 3𝑐 + 7
Number of
Heat
streams
exchanged
𝑁𝑟 = 𝑐 + 1
𝑁𝑖 = 2𝑐 + 6
e.g.
2*(c+2): definition of inlet streams
1: temperature
1: pressure
𝑁𝑟 = 𝑐 + 1
𝑁𝑖 = 𝑐 + 4
e.g.
(c+2): definition of inlet stream
1: temperature
1: pressure or heat
𝑁𝑟 = 𝑐 + 𝑐 + 1 + 2 = 2𝑐 + 3
T and P constant among
Mass balance Phase equilibrium Energy Balance outlet streams
𝑁𝑖 = 𝑐 + 4
e.g.
(c+2): definition of inlet stream
1: temperature
1: pressure or heat
Equilibrium or
𝑁𝑟 = 𝑐 + 1 + 𝑛 conversion fixed
𝑁𝑖 = 𝑐 + 4 𝑁𝑖 = 𝑐 + 5
If reactions are If reactions are not
specified specified
e.g. e.g.
(c+2): definition of inlet stream (c+2): definition of inlet stream
1: temperature 1: temperature
1: pressure or heat 1: pressure or heat
1: equilibrium constant
𝑁𝑖 = 2𝑐 + 6 q
Ln+1
Vn
e.g.
2(c+2): definition of inlet
streams
1: temperature
1: pressure Ln
Vn-1
𝑁𝑖 = 3𝑐 + 8 q Ln+1
Vn
e.g.
3(c+2): definition of inlet
streams
1: temperature
1: pressure Ln
Vn-1
e.g.
2(c+2): definition of inlet
streams
1: heat exchanged
1: pressure Ln Vn-1
1: mass flow of stream S
𝑁𝑟 = 2 ∗ 𝑐 + 2 ∗ (𝑛 − 1)
1
𝑁𝑖 = 2𝑛 + 2𝑐 + 5
Streams 2 3 4 5 4 3 Total = 21
𝑁𝑟 = 9 ∗ 𝑐 + 2
𝑁𝑖 = 2 ∗ (𝑛 + 𝑚) + 𝑐 + 13
e.g.
c+2: definition of feed
n,m: P ∀ stage
n,m: Heat loss ∀ stage
2: number of stages n,m
1: Reflux ratio
2: LK,HK recovery
6: T,P condenser, reboiler, splitter
ΔH
∆𝑄2 = ∆𝑄𝑡𝑜𝑡 MIXER
3
Process Simulation – Maurizio Fermeglia Trieste, 28 September, 2021 - slide 33
DoF of a Cyclic Process
Gas-phase dehydrogenation of isobutane to isobutene
◼ Reaction: A → B + C C4H10 → C4H8 + H2
◼ Pure isobutane (A) is mixed adiabatically with a recycle containing 90% mole
isobutane and the balance isobutene (B)
◼ Single pass isobutane conversion is 35%
◼ All hydrogen (C) and 10% of the isobutane (as well as some isobutene) is
sent to another part of the plant
◼ Qr and Qs are the net rate of heat transfer in the reactor and separator
Na2 (mol A/s)
Nb2 (mol B/s)
Na1 (mol A/s) Na3 (mol A/s)
Nc2 (mol C/s)
P-1 Nb1 (mol B/s) Nb3 (mol B/s)
100 mol A/s
Nc3 (mol C/s)
20 °C 90 °C
T1 °C
30 °C
E-2 Separator
E-3 Pump
E-1 Reactor
N4 (mol/s)
0.9 (mol A /mol)
0.10 (mol B/mol)
85°C
Overall process
◼ 7 local DoF (1+3+3); 7 ties (Na1, Nb1, Na2, Nb2, Nc2, N4, T1 counted twice) → 0 DoF
E-2 Separator
E-3 Pump
E-1 Reactor
N4 (mol/s)
0.9 (mol A /mol)
0.10 (mol B/mol)
85°C
E-2 Separator
E-3 Pump
E-1 Reactor
N4 (mol/s)
0.9 (mol A /mol)
0.10 (mol B/mol)
85°C
Process R&D
Process
Plant Retrofits
Synthesis
PROCESS SIMULATION
Process
Process &
Improvement -
Equipment
Analysis
Design
Process Simulation – Maurizio Fermeglia Simulation Trieste, 28 September, 2021 - slide 40
The “Life Cycle” of a process
Process
Definition
LCA Software
Process
Optimization
Control Control
System Design system and Sustainability
Process Evaluation
Integration
Dynamic Data
simulation for reconciliation
training and model
Start Up of the maintenance
Plant
process Operation
Evaluate Generate
Ideas Ideas
Does it work?
Set objective
Process Process
targets
Quantify Simulation Synthesis
«real» impact
How to change
the system?
Propose process
configuration &
Operating design (Pinch
Plant technology, train
of columns, …)
Goals definition
Process
& scope
design
Material and
energy
balances
interpretation
Results
Inventory
analysis
Product
design
Multiscale
modeling for
material
Environmental
properties
impact analysis