Dynamic Simulators and MPC PDF
Dynamic Simulators and MPC PDF
Dynamic simulators
& Model based process control
Jari Lappalainen
VTT
Systems modelling and simulation
Overview
Introduction
Diversity of DS
Dynamic process simulation
MPC
Commercial dynamic simulation tools
Application examples
02/12/2015 2
Dynamic simulation: mimicking time-varying
behaviour of a system by computational models
Multibody systems
Dynamic process
simulation
Source: Wikipedia
Discrete
simulation Source: Wikipedia
System Dynamics
CFD
Courtesy of
Valmet
3
Source: Labome.org
Modelling and simulation at and across several
time and magnitude scales
Global
A process industry simulation point-of-view
phenomena
global
Process unit
operation
e.g. a digester
Basic physical
microscopic
5
Different needs − different methods and models
Rigourous thermal
hydraulics, Accurate
reactions, heat physical
transfer, diffusion,
description
material properties
details
A piece of Large
equipment system
Scope of modelling
Nozzle, reactor, Level Process plant,
pump, header, power plant,
mixer, experimental system dynamics,
set-up, .. …
Time
start-up
shut-down
load changes,
grade changes
7
Example: Kappa profile in continuous pulp digester
02/12/2015 9
Where dynamic behaviour comes from?
Paper mill example
Period Time
10 min 1 min 10 s 1s 100 ms 10 ms 1 ms
Vibration
Slow control Fast control Mechanical
loops loops condition
Surges Pulsation
Ref: Cutshall, K.,
Nature of paper
Stock preparation Short circulation variation, TAPPI Wet
Pulp handling End Operations Short
Dilution circulation Course, TAPPI Press,
Memphis 1997, pp.
371-402.
Headbox and former
process
Fidelity of
physical description Single
Equipment
3D CFD Design and
Optimization
2D CFD
All realistic
Dynamic simulation
controls Piping
Detailed Engineering,
incl. noise hydrodynamics
Control Design,
Consistency Consistency DCS Checkout,
controls dynamics Operator Training,
Troubleshooting
Level Simple tank
controls dynamics
No Conceptual Design
controls Steady-state
Comprehensiveness
of simulation
11
Selecting detail of modelling
Note, the level of details greatly affect the time needed for M&S
02/12/2015 12
Model Predictive Control
13
MPC
02/12/2015 14
MPC
Source: wikipedia
02/12/2015 15
Dynamic simulation
tools
16
Requirements for commercial
Dynamic Simulator
02/12/2015 17
What is needed to build up a successful DS
model
The purpose of the model (the questions !)
Capable M&S tool(s)
Modelling input, e.g. P&I-diagrams, lay-outs, equipment data-sheets,
automation diagrams, system descriptions, ..
And basic understanding of how the process works and is operated
Collaboration with process owner, and preferably support team for
modelling
To get quick answers and feedback for intermediate results
Validation
Measurement data (from the automation system) with the relevant
manipulated, measured and controlled variables, sampling frequency, and
transients
Experts for qualitative validation
02/12/2015 18
Verification and validation
02/12/2015 19
Examples of commercial dynamic simulation tools
02/12/2015 21
Main domains in Apros modelling & simulation
Process
systems Automation
systems
Electrical
systems
Thermal Calculation
hydraulics order
Execution Equivalent
bands circuit
approach
02/12/2015 22
Modelling paradigms
prosess model
02/12/2015 23
Mechanistic models
Mechanistic models combine
02/12/2015 24
Discretisation with respect to time
Explicit
explicit Semi-implicit
semi-implicit Implicit
implicit
methods methods
methods methods
methods
2 25
Dynamic conservation equations (3-eq model)
Dynamic conservation equations for mass, momentum, energy and mass fractions
are solved simultaneously for the entire model network
A Av
v Mass balance 0
t z
+ component balance
02/12/2015 26
Source terms
02/12/2015 27 27
Dynamic conservation equations (6-eq model)
02/12/2015 28
Main principles (6-eq model)
Thermal hydraulics
02/12/2015 29 29
Space discretization
02/12/2015 30 30
Edward’s pipe
02/12/2015 31 31
Edward’s pipe
Thermal hydraulics
02/12/2015 32 32
Edward’s pipe
Thermal hydraulics
02/12/2015 33 33
Edward’s pipe
Void fraction in the middle of the pipe
Thermal hydraulics
02/12/2015 34 34
Pipe flow: Overlapping Grid Approach
flow flow
p pressure pressure pressure
h entalphy entalphy entalphy
ci massfractions massfractions massfractions
dm/dt dm/dt
h dm/dt h dm/dt
ci dm/dt ci dm/dt
Homogeneous Flow Model (3-eq): Same temperature of different phases within the
control volume and same velocity of the phases in the homogeneous flow between
control volumes.
Separate Phase Nonequilibrium Model (6-eq): Different temperature of each phase in
same control volume and different velocities of the separate phase flows between
control volumes.
02/12/2015 35
Heat Diffusion and Transfer
T T T
02/12/2015 36
Matrix Solvers for Large Systems Compared
CPU
time
per Inverse
node of A
Full matrix
Ax=b
Sparse solution
of Ax=b
Number
10 of nodes
The solution time per node is presented as a function of the number of nodes.
The sparse matrix solver seems to be faster than that of a full matrix solver for
dimensions N exceeding 10. The inverse of a matrix takes N times more time
than a single solution of x. The inverse of a sparse matrix is usually a full
matrix and also the multiplication x = A-1b could take much time.
37
Building up a process model
02/12/2015 38
Automation is an essential part of a dynamic
simulation model
APROS MODEL
PROCESS MODEL
AUTOMATION
MODEL
MEASUREMENTS ACTUATORS
CONTROL
&
LOGIC
02/12/2015 39 39
Options to simulate the automation system
02/12/2015 40 40
Simulation model and automation
Measurements &
Actuators
Process Automation
For example
Advanced
Process Automation
Process Control
unit
02/12/2015 41
Calculation principles
By default
Automation time step = thermal hydraulics time step
Calculation order solved automatically
Current time step Current time step decreases =>
(= max time step) also Automation is calculated more frequently
e.g. 200 ms
02/12/2015 43 43
Apros Automation Components – Binary
Pulse_generator
Ab_converter
02/12/2015 44
44
Controller
02/12/2015 45 45
Controller
02/12/2015 46
46
Dynamic simulation can be used for basic
automation engineering too
Efficient
environment to
design control and
logic systems
…and evaluate the
functionality by
simulation
Interfaces with
other process and
automation design
tools
02/12/2015 47
Application examples
48
Various uses of Apros simulation
Performance
Evaluation of Safety analysis analysis
new process
concepts Automation testing
Operator
training
Integrated Process &
Control design Control room Trouble-
design evaluation shooting
Development of
new control strategies
Automation Optimisation
engineering
02/12/2015 49
Various uses and users of Apros simulation
Process and automation engineering, testing, analysis
Fortum, Metso, Alstom, Andritz, Amec Foster Wheeler, Doosan, Eskom, Tecnatom,
TVO, EON, TÜV, nuclear safety authorities, …
Training simulators
Nuclear power plants, conventional power plants, pulp&paper mills, …
Academic research and education
Nuclear, combustion, solar, fuel cells, smart grids, pulp&paper, control engineering, …
VTT, Aalto University, Oulu University of Technology, Lappeenranta University of
Technology, TU Wien, TU Budapest, TU Darmstadt, KTH, TU München, …
02/12/2015 50
Examples of large Apros models:
diagrams 1060
th2 nodes 3540
th6 nodes 1065
valves 2617+347+316+600+144 = 4024
(sv+cv+bv+chv+cv+sf = total)
pumps 301
actuators 2103
analog signals 27272
binary signals 111613
controllers 245 + 18 = 263
es nodes 295
es switches 188
DB size 58 MB (in .gz format)
Introduction
02/12/2015 51 51
Case 1: Process & control concept evaluation
• Example: a CCS capable power plant
concept including
• Interface to ASU
• Supercritical CFB boiler
• Turbine island
• CPU
• Controls
• Also as co-simulation with Aspen Plus
Dynamics
Simulation example:
Mode switch from air
firing to oxy firing and
back
Ref: Lappalainen et al.
International Journal of
Greenhouse Gas Control,
28(2014)11-24
02/12/2015 52
Case 2: Loviisa NPP Automation Renewal
02/12/2015 53
Case 3: Distributed energy production and consumption
Energy systems in the Smart City – Area level Plans & Simulation
54 54
Case 4: Grade Change Development at Board Mill
Mill data
• Stora Enso Imatra mills
• ply board machine
• packaging & graphical boards
• 200 000 tons/year
• width 5.6 m
• speed 200-450 m/min
• basis weight 170-350 g/m2
• grade change freq > 1 per day
55
Case 4: Target of Simulation Studies
02/12/2015 56
Case 4: Automatic Grade Change
headbox slice
openings 1, 2 & 3
steam pressures
57
Case 4: AGC Involves Ramping and Mutual
Timing of Control Variables
GRADE #1 GRADE #2
BASIS WEIGHT
MOISTURE
QUALITY ASH CONTENT
VARIABLES THICKNESS
COLOR
...
MACHINE SPEED
THICK STOCK FLOW(S)
STEAM PRESSURE(S)
MANIPULATED SLICE OPENING(S)
VARIABLES JET-WIRE RATIO(S)
RAW MATERIAL RECIPE
LINEAR LOADS IN WET PRESSING
COATING SETTINGS
REFINER LOADS
10 20 Time / min ...
58
Case 4: Target is to Minimize Off-spec Production
Specification of
GRADE A2
Specification of
GRADE A1
Off-spec
59
Case 4: Tuning Parameters for Each AGC Variable
Start Stop
delay delay
(up, down) Maximum (up, down)
Ramping
Rate
(up, down)
GC Start
60
Case 4: What-if Studies Varying Parameter Values
Procedure:
1. Test performance in simulator
using various parameter values
2. Take parameters into use
61
Case 4: Scope of the Simulator
62
Case 4: Simulator Performance Tested Against
Measurements
SPEED [m/min]
MANIPULATED VARIABLES
360
350
QUALITY VARIABLES
340
330
320
940 290 11
920 280 10
-300 0 300 600 900 1200
900 9
THICK STOCK FLOWS [l/s]
880
140 8
860 130
7
840 120
6
110
820
100 5
800 90
-300 0 300 600 900 1200 4
80
70 3
JET-WIRE-RATIO 2
60
-300 0 300 600 900 1200
1.15
OVEN DRY BASIS WEIGHT [g/m²]
50
1.13 300
40
-300 0 300 600 900 1200 290
1.11
SLICE OPENINGS [mm]
1.09 50 280
270
1.07
45 260
1.05
250
1.03 40
240
1.01
35 230
0.99
220
0.97 30
210
0.95
200
-300 0 300 600 900 1200 25
-300 0 300 600 900 1200
-300 0 300 600 900 1200
63
Case 4: Results
02/12/2015 64
Benefits of dynamic process simulation
02/12/2015 65
Future trends
02/12/2015 66
Summary