1 linné flow centre kth mechanics isaac newton institute workshop 8-12 september, 2008 dan...
TRANSCRIPT
![Page 1: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/1.jpg)
1
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Dan Henningson
collaborators
Shervin Bagheri, Espen ÅkervikLuca Brandt, Peter Schmid
Input-output analysis, model reduction and control applied to the Blasius boundary layer - using balanced modes
![Page 2: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/2.jpg)
2
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Linearized Navier-Stokes for Blasius flow
Discrete formulationContinuous formulation
![Page 3: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/3.jpg)
3
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Input-output configuration for linearized N-S
![Page 4: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/4.jpg)
4
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Solution to the complete input-output problem
• Initial value problem: flow stability
• Forced problem: input-output analysis
![Page 5: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/5.jpg)
6
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Input-output operators
• Past inputs to initial state: class of initial conditions possible to generate through chosen forcing
• Initial state to future outputs: possible outputs from initial condition
• Past inputs to future outputs:
![Page 6: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/6.jpg)
7
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Most dangerous inputs, creating the largest outputs
• Eigenmodes of Hankel operator – balanced modes
• Observability Gramian
• Controllability Gramian
![Page 7: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/7.jpg)
8
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Controllability Gramian for GL-equation
• Correlation of actuator impulse response in forward solution
• POD modes:
• Ranks states most easily influenced by input
• Provides a means to measure controllability
Input
![Page 8: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/8.jpg)
9
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Observability Gramian for GL-equation
• Correlation of sensor impulse response in adjoint solution
• Adjoint POD modes:
• Ranks states most easily sensed by output
• Provides a means to measure observability
Output
![Page 9: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/9.jpg)
11
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Balanced modes: eigenvalues of the Hankel operator
• Combine snapshots of direct and adjoint simulation
• Expand modes in snapshots to obtain smaller eigenvalue problem
![Page 10: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/10.jpg)
12
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Snapshots of direct and adjoint solution in Blasius flow
Direct simulation: Adjoint simulation:
![Page 11: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/11.jpg)
13
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Balanced modes for Blasius flow
adjoint forward
![Page 12: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/12.jpg)
15
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Model reduction
• Project dynamics on balanced modes using their biorthogonal adjoints
• Reduced representation of input-output relation, useful in control design
![Page 13: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/13.jpg)
16
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Impulse response
DNS: n=105
ROM: m=50
Disturbance Sensor
Actuator Objective
Disturbance Objective
![Page 14: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/14.jpg)
17
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Frequency response
DNS: n=105
ROM: m=80 m=50 m=2
From all inputs to all outputs
![Page 15: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/15.jpg)
18
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Optimal Feedback Control – LQG
controller
Find an optimal control signal (t) based on the measurements (t) such that in the presence of external disturbances w(t) and measurement noise g(t) the output z(t) is minimized.
Solution: LQG/H2
cost function
KLw z g(noise)
![Page 16: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/16.jpg)
19
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
LQG controller formulation with DNS
• Apply in Navier-Stokes simulation
![Page 17: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/17.jpg)
20
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Performance of controlled system
Noise
Sensor
Actuator
Objective
controller
![Page 18: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/18.jpg)
21
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Performance of controlled system
Noise Sensor Actuator Objective
Control offCheap ControlIntermediate controlExpensive Control
![Page 19: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/19.jpg)
22
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Conclusions
• Input-output formulation ideal for analysis and design of feedback control systems – applied in Blasius flow
• Balanced modes Obtained from snapshots of forward and adjoint solutions Give low order models preserving input-output relationship
between sensors and actuators
• Feedback control of Blasius flowReduced order models with balanced modes used in LQG controlController based on small number of modes works well in DNS
• Outlook: incorporate realistic sensors and actuators in 3D problem and test controllers experimentally
![Page 20: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/20.jpg)
23
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
![Page 21: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/21.jpg)
24
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Message
• Need only snapshots from a Navier-Stokes solver (with adjoint) to perform stability analysis and control design for complex flows
• Main example Blasius, others: GL-equation, jet in cross-flow
![Page 22: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/22.jpg)
25
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Outline
• Introduction with input-output configuration
• Matrix-free methods using Navier-Stokes snapshots
• The initial value problem, global modes and transient growth
• Particular or forced solution and input-output characteristics
• Reduced order models preserving input-output characteristics, balanced truncation
• LQG feedback control based on reduced order model
• Conclusions
![Page 23: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/23.jpg)
26
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Background
• Global modes and transient growthGinzburg-Landau: Cossu & Chomaz (1997); Chomaz (2005)Waterfall problem: Schmid & Henningson (2002)Blasius boundary layer, Ehrenstein & Gallaire (2005); Åkervik et al. (2008)Recirculation bubble: Åkervik et al. (2007); Marquet et al. (2008)
• Matrix-free methods for stability propertiesKrylov-Arnoldi method: Edwards et al. (1994)Stability backward facing step: Barkley et al. (2002)Optimal growth for backward step and pulsatile flow: Barkley et al. (2008)
• Model reduction and feedback control of fluid systemsBalanced truncation: Rowley (2005)Global modes for shallow cavity: Åkervik et al. (2007)Ginzburg-Landau: Bagheri et al. (2008)
• Invited session on Global Instability and Control of Real Flows, Wednesday 8-12, Evergreen 4
![Page 24: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/24.jpg)
27
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
The forced problem: input-output
• Ginzburg-Landau example• Input-output for 2D Blasius configuration• Model reduction
![Page 25: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/25.jpg)
28
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Input-output analysis
• Inputs: Disturbances: roughness, free-stream turbulence, acoustic wavesActuation: blowing/suction, wall motion, forcing
• Outputs: Measurements of pressure, skin friction etc.
• Aim: preserve dynamics of input-output relationship in reduced order model used for control design
![Page 26: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/26.jpg)
29
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Feedback control
• LQG control design using reduced order model• Blasius flow example
![Page 27: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/27.jpg)
30
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
LQG feedback control
Estimator/ Controller
Reduced model of real system/flow
cost function
![Page 28: 1 Linné Flow Centre KTH Mechanics Isaac Newton Institute Workshop 8-12 September, 2008 Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca](https://reader037.vdocument.in/reader037/viewer/2022110402/56649e4a5503460f94b3ddcc/html5/thumbnails/28.jpg)
31
Linné Flow CentreKTH Mechanics
Isaac Newton Institute Workshop 8-12 September, 2008
Riccati equations for control and estimation gains