Optimal Control of Plane Poiseuille Flow
|
|
- Carmel Andrews
- 5 years ago
- Views:
Transcription
1 Optimal Control of Plane Poiseuille Flow Dr G.Papadakis Division of Engineering, King s College London george.papadakis@kcl.ac.uk AIM Workshop on Flow Control King s College London, 17/18 Sept 29 AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 1/62
2 Contents Background and motivation Basic equations and spectral model Linear state-space model (plant) Initial conditions Synthesis of state feedback controllers, state estimators and output feedback controllers Controller implementation in full FV Navier-Stokes solver Recapitulation and future work AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 2/62
3 Background and motivation Ability to alter flow can have tremendous benefits Reduction of flow unsteadiness can benefit: Ocean ships (consumption 2.1bn barrels of oil pa) Aircrafts (consumption 1.5bn barrels of jet oil pa) Cars, trucks, trains, pipelines etc Increase of flow unsteadiness can benefit: Combustion efficiency Product quality that depends on mixing AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 3/62
4 Active control with feedback (within a linear framework) Application in Poiseuille Flow Wall Plane Poiseuille Flow Wall + Flow Disturbance Wall Normal, y Streamwise, x Actuation Sensing (Controller) AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 4/62
5 Control of Poiseuille flow Simple base flow, geometry, boundary conditions Linearly unstable R>5772, Transition R>1 Steps: Selection of state variables, actuation and sensing Transformation to state-space form Synthesis of linear optimal controllers to minimise Transient Energy Growth Implementation and testing of controllers in a NS solver AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 5/62
6 Perturbation Equations around a base flow Perturbation form of Navier-Stokes and Continuity eq. ( ) u + Ub u + ( u ) u +( u ) U b = p + 1 R 2 u u = Non-linear equations Divergence free velocity field Homogeneous (zero) wall boundary conditions on u Base flow U b parabolic AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 6/62
7 Linearised perturbation equations General form ( ) u + Ub u +( u ) U b = p + 1 R 2 u u = Parallel flows only: Wall-normal Velocity-vorticity v-η form η = u z ( 2 v) t η t w x ( 2 v) + U b x + U b y v z + U b 2 U b y 2 v x 1 R 2 ( 2 v)= η x 1 R 2 η = Boundary conditions at walls: v =, v/ y =, η = AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 7/62
8 Discretisation of v-η form Basis for state-space form: Spectral discretisation Fourier in streamwise and spanwise directions x, z One wavenumber pair (α, β) selected ( U b k 2 2 U b y 2 = j ( 3ṽ α y 2 t (k 2 = α 2 + β 2 ) v = R (ṽ(y, t)e j(αx+βz)) ) k4 ṽ + jrα ) ṽ k2 (jαu b + k2 R t ) η 1 R ) (U b + 2k2 2ṽ jrα y 2 2 η y + jβṽ U b 2 y = η t ( 1 ) 4ṽ jrα y 4 AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 8/62
9 Discretisation of v-η form (cont) Chebyshev polynomials Γ n (y) in wall-normal direction y ṽ(y, t) = N Γ n (y)a v,n (t) Evaluation at collocation pts y =cos(nπ/n),n= N ṽ(y,t) Γ (1)... Γ N (1) a v,. = ṽ(y N,t) Similarly for η } Γ ( 1)... {{ Γ N ( 1) } D a v,n AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 9/62
10 Discretisation of v-η form (cont) Calculation of first order derivatives: ṽ (y, t) = N n= a v,n (t)γ n(y) ṽ (y,t). ṽ (y N,t) = Γ (1)... Γ N (1) } Γ ( 1)... {{ Γ N ( 1) } D1 a v,. a v,n Similarly we define matrices D2, D3, D4 for higher order derivatives. AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 1/62
11 Discretisation of v-η form (cont) Substituting to the v-η equations: L 11 L 12 ȧ = A 11 A 12 a L 21 L 22 A 21 A 22 A 11 = a =(a v,...a v,n a η,...a η,n ) T ) ( αūk 2 αū Ik4 jr A 12 = [], A 21 = βū D, A 22 = D + ) (αū + 2Ik2 ) (αū + Ik2 jr jr D D2 jr L 11 = j ( k 2 D + D2 ), L 12 =[], L 21 =[], L 22 = jd D2 D4 jr AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 11/62
12 Boundary conditions (without control) Lower wall ṽ = ṽ/ y = Similarly for upper wall Similarly for η N n= N n= Γ n ( 1)a v,n (t) = Γ n( 1)a v,n (t) = These b.c. replace the first, second, last, next-to-last rows. Leads to spurious eigenvalues (easily mapped to an arbitrary location in the complex plane) AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 12/62
13 Introduction of wall transpiration We aim to control the flow via wall transpiration. New b.c. Lower wall: ṽ( 1,t)= q l (t) Upper wall: ṽ(1,t)= q u (t) q l (t), q u (t) are unknown functions (to be determined) For both walls: ṽ/ y =,η = AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 13/62
14 Transformation of variables Introduction of variable ṽ h (y, t) that has homogenous b.c. ṽ(y, t) =ṽ h (y, t)+f u (y) q u (t)+f l (y) q l (t) We want ṽ h (±1,t)=ṽ h (±1,t)= This can be achieved if f u (1) = f l ( 1) = 1,f u ( 1) = f l (1) = f u(±1) = f l (±1) = Suitable functions (3rd order polynomials) f u (y) = y3 +3y +2,f l (y) = y3 3y AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 14/62
15 Novel Basis functions To eliminate spurious modes, we introduce novel basis functions Γ DN n (y) and Γ D n We require them to satisfy the conditions: Γ D n (y = ±1) = Γ DN n (y = ±1) = Γ DN n (y = ±1) = ṽ(y, t) = η(y, t) = N n= Γ DN n (y)a v,n (t)+f u (y)ṽ u (t)+f l (y)ṽ l (t) N Γ D n (y)a η,n (t) n= AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 15/62
16 Form of new basis functions Γ D n,odd = Γ n Γ 1 Γ D n,even =Γ n Γ 2 Γ DN n,odd = (Γ n+4 Γ 1 ) (n +3) 2 (Γ n+2 Γ 1 )/(n +1) 2 Γ DN n,even = (Γ n+4 Γ 2 ) ((n +3) 2 1) (Γ n+2 Γ 2 )/ [ (n +1) 2 1 ] 1.5 f u f l Γ D last y.5 1 Γ DN 1 Γ DN Γ D Γ DN last Γ D 1 AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 16/62
17 State-space form L 11 L 12 ȧ = A 11 A 12 a + B 11 B 12 q E 11 E 12 q L 21 L 22 A 21 A 22 B 21 B 22 E 21 E 22 }{{} effect of transpiration where q =( q u, q l ) T and (N +1) 1 vectors are B 11 = B 12 = ) ( αūk 2 αū Ik4 jr ) ( αūk 2 αū Ik4 jr f u + f l + ( ) αū + 2k2 I jr ( ) αū + 2k2 I jr B 21 = βū f u, B 22 = βū f l E 11 = j ( k 2 f u + f ) u, E12 = j ( k 2 f l + f ) l E 21 = (), E 22 =() f u f l AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 17/62
18 State-space form (cont) Invert L ȧ = L 1 Aa + L 1 Bq L 1 E q Input vectors q and q are not independent. Recast system so that q is part of the state variable vector and control input is q. ȧ q = L 1 A L 1 B a q + L 1 E q I New system matrix has the same eigenvalues as L 1 A plus two eigenvalues equal to. AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 18/62
19 Final form of State-Space Model General form required for controller synthesis;- X = AX + BU, Y = CX For our problem: ȧ q }{{} X = ỹ = L 1 A L 1 B }{{} A ( ) C D a }{{} q C }{{} X a q }{{} X + L 1 E I } {{ } B q }{{} U AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 19/62
20 Output in physical space Measurement of shear stresses in x,z directions in the 2 walls. τ yx = 1 ( u R y + v ) x τ yz = 1 ( w R y + v ) z Since v is known it can be subtracted out So the measurement vector can be defined as: y = 1 ( ) T u u w w y y y y R y=+1 y= 1 y=+1 y= 1 AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 2/62
21 Output in Fourier space ỹ = 1 R ũ y ũ y w y w y y=+1 y= 1 y=+1 y= 1 = 1 R j α 2 + β 2 α 2 ṽ y 2 α 2 ṽ y 2 β 2 ṽ y 2 β 2 ṽ y 2 β η y β η y + α η y + α η y y=+1 y= 1 y=+1 y= 1 or ỹ = j α 2 + β 2 C 11 C 12 C 21 C 22 a + j α 2 + β 2 D 11 D 21 q ỹ = ( C D ) a q AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 21/62
22 Final form of State-Space Model: a recap States X : Coefficients of basis a n and Wall-normal velocities at walls q u, q l Input U: rate of change of wall-normal velocities at walls q u, q l Output Y: Wall-shear stress measurements τ u, τ l A has two integrators, associated with transpiration from either wall, free of spurious eigenvalues We are ready to synthesize controllers now! AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 22/62
23 Initial Conditions Worst initial conditions lead to non-normal growth Non-orthogonal system matrix A (Trefethen 1993) Variational method used (Butler and Farrell 1992, Bewley 1998) States X (t) =Ξe Λt χ, Energy E = X (t) T QX (t) Transient Energy Growth θ = E(τ)/E() from θ ( Ξ T QΞ ) ) χ = (e ΛT τ Ξ QΞe Λτ χ Hence X () = Ξχ u() AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 23/62
24 Test Cases Two cases representing early transition 1) α =1,β =,R = 1 Tollmien Schlichting waves One Unstable eigenvalue 2) α =,β =2.44,R = 5 Streamwise vortices Stable, but the largest transient energy growth Both are 2-D cases AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 24/62
25 case 1: Comparison with Orszag(1971) 1 λ i λ i (Orszag) / λ i (Orszag) mode i AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 25/62
26 case 1: Poles and zeros.1 Poles Zeros I R AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 26/62
27 case 1: Unstable ṽ eigenvector Comparison with Thomas(1953) R I Thomas R Thomas I.8.6 v y AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 27/62
28 case 1: ṽ eigenvectors Real Part (solid) and Imaginary Part (dot-dashed) 1 y y 1 R I Mode i AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 28/62
29 case 1: ũ eigenvectors Real Part (solid) and Imaginary Part (dot-dashed) 1 (Scaled by 1/25 c.f. ṽ eigenvectors) y y 1 R I Mode i AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 29/62
30 Transpiration mode (m=2) Re=5 u contours y= y= 1 x= v contours x=6.23 y= y= 1 x= x= Results from linear simulations and solution of NS (in-house code) are identical. AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 3/62
31 Case 1: Observability Current Model Bewley s Data 1 4 κ o Mode i AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 31/62
32 Open loop, case 1, Transient energy From Initial Conditions X worst Scaled to E = N=1 N=2 N=3 N=4 N=5 6 E t AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 32/62
33 Case 2: Poles and zeros.5.4 Poles Zeros I R AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 33/62
34 Case 2: ṽ eigenvectors (imaginary) 1 y I y Mode i AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 34/62
35 Case 2: Observability 1 8 Current Model Bewley s Data 1 6 κ o Mode i AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 35/62
36 Open loop, case 2, Transient energy From Initial Conditions X worst Scaled to E = N=1 N=2 N=3 N=4 N=5 E t AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 36/62
37 Optimal State Feedback (LQR controller) AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 37/62
38 Optimal State Feedback (LQR controller) Given the real system;- X = AX + BU, Y = CX State feedback control signal to minimize;- ( X (t) T QX (t)+u(t) T RU(t) ) dt Given by U = KX where K = R 1 B T X X = X T from ARE A T X + XA XBR 1 B T X + Q = AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 38/62
39 Weighting Matrix Q Choose Q to form energy density E; E = X T QX = 1 V E = 1 8k y= 1 2π/α x= 2π/β z= (k 2 ṽ T ṽ + ṽ T ṽ y y + ηt η)dy (u 2 + v 2 + w 2 ) dz dx dy 2 Curtis-Clenshaw quadrature for integration E = 1 8k 2 X T T T Q Q T +( T/ y) T Q ( T/ y) X X T QX Q contains quadrature weights T converts states to ṽ, η at collocation points AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 39/62
40 Weighting Matrix R Choose R = r 2 I Max energy vs r plotted from linear simulations Large r leads to small control effort (and larger energy) Small r leads to larger control effort (and smaller energy) r =.25 for Case 1 and r = 128 for Case 2. AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 4/62
41 Case 1: Open loop v.s. LQR control N=1 N=2 N=3 N=4 N= N=3 N=4 N=5 N=7 N=1 E E t t without control with LQR control AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 41/62
42 Case 2: Open loop v.s. LQR control N=1 N=2 N=3 N=4 N= N=3 N=4 N=5 N=7 N=1 E t without control E t with LQR control AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 42/62
43 How does the controller affect the flow? u(y, z, t) =u v U b t (Ellingsen and Palm (1975)) (a) t = (b) t = τ/8 (c) t = τ (d) t =4τ Figure from Bewley and Liu (1998) (no control) v component (dashed) magnified by 25 as compared to η (solid) Solid contours: u>, Dashed contours: u< t = τ corresponds to maximum transient energy growth AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 43/62
44 How does the controller affect the flow? (cont) u(y, z, t) =u v U b t (Ellingsen and Palm (1975)) (a) t = (b) t = τ/8 (c) t = τ (d) t =4τ Figure from Bewley and Liu (1998) (with control) v component (dashed) magnified by 25 as compared to η (solid) Solid contours: u>, Dashed contours: u< t = τ corresponds to maximum transient energy growth AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 44/62
45 Optimal State Estimation AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 45/62
46 Optimal State Estimation Given the real estimator;- X est = AX est + BU + L (Y CX est ) The optimal L to minimize;- E { } [X X est ] T [X X est ] Is given by L = Y CV 1 where Y = Y T from ARE AY + Y A T Y C T V 1 CY + W = AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 46/62
47 Weighting Matrices W, V W represents covariance of process noise Transform system states to values at collocation points Choose W =(1 y 1 ) 2 (1 y 2 ) 2, where y 1,y 2 are collocation points Disturbances near the centreline more variable/larger Better estimation than W = I V represents covariance of measurement noise Choose V = αi, α by inspection AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 47/62
48 Case 1: Performance of tuned estimator (open loop) E E, N=1 E Est, N=3 E Est, N=4 E Est, N=5 E Est, N=7 E Est, N= t AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 48/62
49 Case 2: Performance of estimator (open loop) E, N=1 E Est, N=3 E Est, N=4 E Est, N=5 E Est, N=7 E Est, N=1 E t AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 49/62
50 Optimal Output feedback AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 5/62
51 Effect of control: A summary Linear Simulation Case OL Time LQR Time LQG Time 1 unstable AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 51/62
52 Implementation of Controller into a Navier-Stokes solver Off-line: Generate matrices A, B, C,Q Synthesize controller matrix K and estimator matrix L (for output feedback) Store K and A, B, C, L (for output feedback) On start-up: Read matrices in NS solver and store Read initial velocity field On each time-step of the solver: For state feedback Calculate X from velocity field Actuation U = KX For output feedback: FFT on τ u,τ l to get Y =( τ u, τ l ) T Run estimator X est = AX est + BU + L (Y CX est ) Actuation U = KX est Integration of U to get q u and q l Inverse FFT to get q u,q l Use q u,q l to set boundary conditions for next time step AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 52/62
53 Navier-Stokes solver Finite-volume discretisation method Second order in space (CDS), implicit 2 nd order in time PISO algorithm, Collocated grid (approx 1 1) Boundary Conditions: Streamwise - cyclic at α Walls - transpiration set by controller Spanwise - cyclic at β Code modified to solve for the non-linear perturbation u about base flow U b AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 53/62
54 Small Perturbations: Cases 1 and 2, LQR Perturbation energy at t =corresponds to ṽ max =1 4 U cl E C1 = (Case 1) E C2 = (Case 2) 1.2 x Non linear E Linear E 2 x Non linear E Linear E E.6 E t t Case 1 (LQR control) Case 2 (LQR control) AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 54/62
55 Large Perturbations that trigger non-linearities Perturbation energy at t =corresponds to ṽ max =1 2 U cl 1 4 E C1 (Case 1) 1 4 E C2 (Case 2) which is around 1 times the transition threshold of Reddy et al (1998) Non linear E Non linear E Est Linear E Linear E Est.12.1 Non linear E Non linear E Est Linear E Linear E Est.5.8 E.4 E t t Case 1 (Open loop) Case 2 (Open loop) AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 55/62
56 Control of large Perturbations (LQR).12.1 Non linear E Linear E Non linear E Linear E E.6 E t t Case 1 (LQR control) Case 2 (LQR control) AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 56/62
57 Multi-wave number control: Motivation Case 2 x E max = at t=123.9 E max /E(t=)= E Σ i=:2 E i E E 1 E E max =1.929e 5 at t=379.5 E max /E(t=)= E Σ i=:2 E i E E 1 E E E t t Open loop, E() = 1 4 E C2 Open loop, E() = E C2 AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 57/62
58 Nominal vs Multi-wave number control Case 2, E() = 1 4 E C2 x E max = at t=87.9 E max /E(t=)=41.39 E Σ i=:2 E i E 9 8 E max = at t=86.1 E max /E(t=)= E Σ i=:2 E i E.8 E 1 E 2 7 E 1 E 2 E.6 E t t Nominal control Multi-wave number control AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 58/62
59 Comparison of nominal vs Multi-wave number control Energy Growth: A summary max t (E(t))/E(t =) Linear-Sized Small Non-Linear Large Non-linear E C2 1 4 E C2 1 6 E C2 Uncontrolled One wavenumber Control Multiple wavenumber Control AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 59/62
60 Control of oblique vortices (3D simulations) α =1, β =1, ṽ max =1 3 U cl 1.5 x 1 4 n. d. E vs n. d. t 1 n. d. E n. d. t Open loop (red line) vs nominal control (green line) and multi-wave number control (blue line). AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 6/62
61 Recapitulation Development of a linear state-space model for plane Poiseuille flow Synthesis of optimal controllers at (α, β) Synthesis of multi-wave number controllers Implementation of controllers in FV Navier-Stokes code. Small Perturbations Spectral linear and FV CFD results very close. Large Perturbations Trigger non-linearities and generate multi wave numbers. Multi-wave number controllers show benefit over single-wave number controllers AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 61/62
62 The End Any Questions? AIM workshop on Flow Control, King s College London, 17/18 Sept 29 p. 62/62
Optimal Control of Plane Poiseuille Flow
Optimal Control of Plane Poiseuille Flow Workshop on Flow Control, Poitiers 11-14th Oct 04 J. Mckernan, J.F.Whidborne, G.Papadakis Cranfield University, U.K. King s College, London, U.K. Optimal Control
More informationFeedback control of transient energy growth in subcritical plane Poiseuille flow
Feedback control of transient energy growth in subcritical plane Poiseuille flow F. Martinelli, M. Quadrio, J. McKernan, J.F. Whidborne LadHyX, École Polytechnique; Politecnico di Milano; King s College,
More informationTransition to turbulence in plane Poiseuille flow
Proceedings of the 55th Israel Annual Conference on Aerospace Sciences, Tel-Aviv & Haifa, Israel, February 25-26, 2015 ThL2T5.1 Transition to turbulence in plane Poiseuille flow F. Roizner, M. Karp and
More informationFeedback control of transient energy growth in subcritical plane Poiseuille flow
Feedback control of transient energy growth in subcritical plane Poiseuille flow Fulvio Martinelli, Maurizio Quadrio, John McKernan and James F. Whidborne Abstract Subcritical flows may experience large
More informationAdvances in Output Feedback Control of Transient Energy Growth in a Linearized Channel Flow
AIAA SciTech Forum 7-11 January 219, San Diego, California AIAA Scitech 219 Forum 1.2514/6.219-882 Advances in Output Feedback Control of Transient Energy Growth in a Linearized Channel Flow Huaijin Yao
More informationStochastic excitation of streaky boundary layers. Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden
Stochastic excitation of streaky boundary layers Jérôme Hœpffner Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden Boundary layer excited by free-stream turbulence Fully turbulent inflow
More informationFeedback Control of Turbulent Wall Flows
Feedback Control of Turbulent Wall Flows Dipartimento di Ingegneria Aerospaziale Politecnico di Milano Outline Introduction Standard approach Wiener-Hopf approach Conclusions Drag reduction A model problem:
More informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More informationFeedback Control of Transitional Channel Flow using Balanced Proper Orthogonal Decomposition
5th AIAA Theoretical Fluid Mechanics Conference 3-6 June 8, Seattle, Washington AIAA 8-3 Feedback Control of Transitional Channel Flow using Balanced Proper Orthogonal Decomposition Miloš Ilak Clarence
More informationReduced-Order Models for Feedback Control of Transient Energy Growth
9th AIAA Flow Control Conference, 5 9 June, Atlanta, Georgia. -9 Reduced-Order Models for Feedback Control of Transient Energy Growth Aniketh Kalur and Maziar S. Hemati Aerospace Engineering and Mechanics,
More informationEECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators.
Name: SID: EECS C28/ ME C34 Final Wed. Dec. 5, 2 8- am Closed book. Two pages of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 2 2 6 3 4 4 5 6 6 7 8 2 Total
More informationAccurate estimation for the closed-loop robust control using global modes
1/21 G. Tissot, GDR CDD 216. 1/21 Accurate estimation for the closed-loop robust control using global modes application to the Ginzburg-Landau equation Gilles Tissot, Jean-Pierre Raymond Institut de Mathématiques
More information4 Shear flow instabilities
4 Shear flow instabilities Here we consider the linear stability of a uni-directional base flow in a channel u B (y) u B = 0 in the region y [y,y 2. (35) 0 y = y 2 y z x u (y) B y = y In Sec. 4. we derive
More informationSystems Theory and Shear Flow Turbulence Linear Multivariable Systems Theory is alive and well
Systems Theory and Shear Flow Turbulence Linear Multivariable Systems Theory is alive and well Dedicated to J. Boyd Pearson Bassam Bamieh Mechanical & Environmental Engineering University of California
More informationLQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin
LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system
More informationRésonance et contrôle en cavité ouverte
Résonance et contrôle en cavité ouverte Jérôme Hœpffner KTH, Sweden Avec Espen Åkervik, Uwe Ehrenstein, Dan Henningson Outline The flow case Investigation tools resonance Reduced dynamic model for feedback
More informationTransient growth on boundary layer streaks. Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden
Transient growth on boundary layer streaks Jérôme Hœpffner Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden Primary instability: TS vs TG Secondary instability of the streaks Movie of streaks
More informationLecture 2 Supplementary Notes: Derivation of the Phase Equation
Lecture 2 Supplementary Notes: Derivation of the Phase Equation Michael Cross, 25 Derivation from Amplitude Equation Near threshold the phase reduces to the phase of the complex amplitude, and the phase
More informationFLOW-NORDITA Spring School on Turbulent Boundary Layers1
Jonathan F. Morrison, Ati Sharma Department of Aeronautics Imperial College, London & Beverley J. McKeon Graduate Aeronautical Laboratories, California Institute Technology FLOW-NORDITA Spring School on
More informationLecture 9. Introduction to Kalman Filtering. Linear Quadratic Gaussian Control (LQG) G. Hovland 2004
MER42 Advanced Control Lecture 9 Introduction to Kalman Filtering Linear Quadratic Gaussian Control (LQG) G. Hovland 24 Announcement No tutorials on hursday mornings 8-9am I will be present in all practical
More informationControl Systems. Design of State Feedback Control.
Control Systems Design of State Feedback Control chibum@seoultech.ac.kr Outline Design of State feedback control Dominant pole design Symmetric root locus (linear quadratic regulation) 2 Selection of closed-loop
More informationModeling flow statistics using convex optimization
Modeling flow statistics using convex optimization Abstract A method is proposed to estimate the covariance of disturbances to a stable linear system when its state covariance is known and a dynamic model
More informationFeedback Control of Boundary Layer Bypass Transition: Comparison of a numerical study with experiments
Feedback Control of Boundary Layer Bypass Transition: Comparison of a numerical study with experiments Antonios Monokrousos Fredrik Lundell Luca Brandt KTH Mechanics, S-1 44 Stockholm, Sweden δ Ω rms L
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationUnderstanding and Controlling Turbulent Shear Flows
Understanding and Controlling Turbulent Shear Flows Bassam Bamieh Department of Mechanical Engineering University of California, Santa Barbara http://www.engineering.ucsb.edu/ bamieh 23rd Benelux Meeting
More informationD(s) G(s) A control system design definition
R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW DATA SET Data provenance: J. Graham 1, M. Lee 2, N. Malaya 2, R.D. Moser 2, G. Eyink 1 & C. Meneveau 1 Database ingest and Web Services: K. Kanov
More informationCourse Outline. Higher Order Poles: Example. Higher Order Poles. Amme 3500 : System Dynamics & Control. State Space Design. 1 G(s) = s(s + 2)(s +10)
Amme 35 : System Dynamics Control State Space Design Course Outline Week Date Content Assignment Notes 1 1 Mar Introduction 2 8 Mar Frequency Domain Modelling 3 15 Mar Transient Performance and the s-plane
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationThe need for de-aliasing in a Chebyshev pseudo-spectral method
The need for de-aliasing in a Chebyshev pseudo-spectral method Markus Uhlmann Potsdam Institut für Klimafolgenforschung, D-442 Potsdam uhlmann@pik-potsdam.de (March 2) Abstract In the present report, we
More informationRobust Control 5 Nominal Controller Design Continued
Robust Control 5 Nominal Controller Design Continued Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 4/14/2003 Outline he LQR Problem A Generalization to LQR Min-Max
More informationActive Control of Instabilities in Laminar Boundary-Layer Flow { Part II: Use of Sensors and Spectral Controller. Ronald D. Joslin
Active Control of Instabilities in Laminar Boundary-Layer Flow { Part II: Use of Sensors and Spectral Controller Ronald D. Joslin Fluid Mechanics and Acoustics Division, NASA Langley Research Center R.
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW AT Re τ = 5200 DATA SET Data provenance: M. Lee 1 & R. D. Moser 1 Database ingest and Web Services: Z. Wu 2, G. Lemson 2, R. Burns 2, A. Szalay
More informationTRACKING AND DISTURBANCE REJECTION
TRACKING AND DISTURBANCE REJECTION Sadegh Bolouki Lecture slides for ECE 515 University of Illinois, Urbana-Champaign Fall 2016 S. Bolouki (UIUC) 1 / 13 General objective: The output to track a reference
More informationApplication of wall forcing methods in a turbulent channel flow using Incompact3d
Application of wall forcing methods in a turbulent channel flow using Incompact3d S. Khosh Aghdam Department of Mechanical Engineering - University of Sheffield 1 Flow control 2 Drag reduction 3 Maths
More informationLinear stability of plane Poiseuille flow over a generalized Stokes layer
Linear stability of plane Poiseuille flow over a generalized Stokes layer Maurizio Quadrio 1,, Fulvio Martinelli and Peter J Schmid 1 Dip. Ing. Aerospaziale, Politecnico di Milano, Campus Bovisa, I-0156
More informationLecture 2 and 3: Controllability of DT-LTI systems
1 Lecture 2 and 3: Controllability of DT-LTI systems Spring 2013 - EE 194, Advanced Control (Prof Khan) January 23 (Wed) and 28 (Mon), 2013 I LTI SYSTEMS Recall that continuous-time LTI systems can be
More informationExam. 135 minutes, 15 minutes reading time
Exam August 6, 208 Control Systems II (5-0590-00) Dr. Jacopo Tani Exam Exam Duration: 35 minutes, 5 minutes reading time Number of Problems: 35 Number of Points: 47 Permitted aids: 0 pages (5 sheets) A4.
More informationChapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control
Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design
More informationOPTIMAL CONTROL SYSTEMS
SYSTEMS MIN-MAX Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University OUTLINE MIN-MAX CONTROL Problem Definition HJB Equation Example GAME THEORY Differential Games Isaacs
More informationOn using the streamwise traveling waves for variance suppression in channel flows
Proceedings of the 2007 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 2007 WeC19.4 On using the streamwise traveling waves for variance suppression
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationEL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1)
EL 625 Lecture 0 EL 625 Lecture 0 Pole Placement and Observer Design Pole Placement Consider the system ẋ Ax () The solution to this system is x(t) e At x(0) (2) If the eigenvalues of A all lie in the
More information3. Identify and find the general solution of each of the following first order differential equations.
Final Exam MATH 33, Sample Questions. Fall 7. y = Cx 3 3 is the general solution of a differential equation. Find the equation. Answer: y = 3y + 9 xy. y = C x + C x is the general solution of a differential
More informationx 1 x 2. x 1, x 2,..., x n R. x n
WEEK In general terms, our aim in this first part of the course is to use vector space theory to study the geometry of Euclidean space A good knowledge of the subject matter of the Matrix Applications
More informationVectors, matrices, eigenvalues and eigenvectors
Vectors, matrices, eigenvalues and eigenvectors 1 ( ) ( ) ( ) Scaling a vector: 0.5V 2 0.5 2 1 = 0.5 = = 1 0.5 1 0.5 ( ) ( ) ( ) ( ) Adding two vectors: V + W 2 1 2 + 1 3 = + = = 1 3 1 + 3 4 ( ) ( ) a
More informationApplied Differential Equation. November 30, 2012
Applied Differential Equation November 3, Contents 5 System of First Order Linear Equations 5 Introduction and Review of matrices 5 Systems of Linear Algebraic Equations, Linear Independence, Eigenvalues,
More informationEE221A Linear System Theory Final Exam
EE221A Linear System Theory Final Exam Professor C. Tomlin Department of Electrical Engineering and Computer Sciences, UC Berkeley Fall 2016 12/16/16, 8-11am Your answers must be supported by analysis,
More informationControl Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) vs Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
More informationChapter 3. State Feedback - Pole Placement. Motivation
Chapter 3 State Feedback - Pole Placement Motivation Whereas classical control theory is based on output feedback, this course mainly deals with control system design by state feedback. This model-based
More information5. Observer-based Controller Design
EE635 - Control System Theory 5. Observer-based Controller Design Jitkomut Songsiri state feedback pole-placement design regulation and tracking state observer feedback observer design LQR and LQG 5-1
More informationMath 304 Answers to Selected Problems
Math Answers to Selected Problems Section 6.. Find the general solution to each of the following systems. a y y + y y y + y e y y y y y + y f y y + y y y + 6y y y + y Answer: a This is a system of the
More informationOptimal Disturbances in Compressible Boundary Layers Complete Energy Norm Analysis
Optimal Disturbances in Compressible Boundary Layers Complete Energy Norm Analysis Simone Zuccher & Anatoli Tumin University of Arizona, Tucson, AZ, 85721, USA Eli Reshotko Case Western Reserve University,
More informationECEEN 5448 Fall 2011 Homework #4 Solutions
ECEEN 5448 Fall 2 Homework #4 Solutions Professor David G. Meyer Novemeber 29, 2. The state-space realization is A = [ [ ; b = ; c = [ which describes, of course, a free mass (in normalized units) with
More informationEE C128 / ME C134 Final Exam Fall 2014
EE C128 / ME C134 Final Exam Fall 2014 December 19, 2014 Your PRINTED FULL NAME Your STUDENT ID NUMBER Number of additional sheets 1. No computers, no tablets, no connected device (phone etc.) 2. Pocket
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationEnergy Transfer Analysis of Turbulent Plane Couette Flow
Energy Transfer Analysis of Turbulent Plane Couette Flow Satish C. Reddy! and Petros J. Ioannou2 1 Department of Mathematics, Oregon State University, Corvallis, OR 97331 USA, reddy@math.orst.edu 2 Department
More informationControl System Design
ELEC4410 Control System Design Lecture 19: Feedback from Estimated States and Discrete-Time Control Design Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science
More information3. Identify and find the general solution of each of the following first order differential equations.
Final Exam MATH 33, Sample Questions. Fall 6. y = Cx 3 3 is the general solution of a differential equation. Find the equation. Answer: y = 3y + 9 xy. y = C x + C is the general solution of a differential
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationControl Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) and Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
More informationÉCOLE DE PHYSIQUE DES HOUCHES. Institut de Mécanique des Fluides Toulouse CNRS Université de Toulouse, also at
ÉCOLE DE PHYSIQUE DES HOUCHES STABILITY OF FLUID FLOWS Carlo COSSU Institut de Mécanique des Fluides Toulouse CNRS Université de Toulouse, also at Département de Mécanique, École Polytechinique http://www.imft.fr/carlo-cossu/
More informationx 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7
Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)
More informationChapter 4: Higher Order Linear Equations
Chapter 4: Higher Order Linear Equations MATH 351 California State University, Northridge April 7, 2014 MATH 351 (Differential Equations) Ch 4 April 7, 2014 1 / 11 Sec. 4.1: General Theory of nth Order
More informationSection 8.2 : Homogeneous Linear Systems
Section 8.2 : Homogeneous Linear Systems Review: Eigenvalues and Eigenvectors Let A be an n n matrix with constant real components a ij. An eigenvector of A is a nonzero n 1 column vector v such that Av
More informationProblem 1: Solving a linear equation
Math 38 Practice Final Exam ANSWERS Page Problem : Solving a linear equation Given matrix A = 2 2 3 7 4 and vector y = 5 8 9. (a) Solve Ax = y (if the equation is consistent) and write the general solution
More informationZeros and zero dynamics
CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)
More information6 OUTPUT FEEDBACK DESIGN
6 OUTPUT FEEDBACK DESIGN When the whole sate vector is not available for feedback, i.e, we can measure only y = Cx. 6.1 Review of observer design Recall from the first class in linear systems that a simple
More informationOverview of the Numerics of the ECMWF. Atmospheric Forecast Model
Overview of the Numerics of the Atmospheric Forecast Model M. Hortal Seminar 6 Sept 2004 Slide 1 Characteristics of the model Hydrostatic shallow-atmosphere approimation Pressure-based hybrid vertical
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationFlow Instability Matlab Tutorial Bangalore, January 2010 Matthew Juniper, Dan Henningson, Peter Schmid 1
Flow Instability Matlab Tutorial Bangalore, January Matthew Juniper, Dan Henningson, Peter Schmid Simple differentiation matrices Imagine a vector that represents the values of a function y at N equispaced
More informationFull-State Feedback Design for a Multi-Input System
Full-State Feedback Design for a Multi-Input System A. Introduction The open-loop system is described by the following state space model. x(t) = Ax(t)+Bu(t), y(t) =Cx(t)+Du(t) () 4 8.5 A =, B =.5.5, C
More informationLinear State Feedback Controller Design
Assignment For EE5101 - Linear Systems Sem I AY2010/2011 Linear State Feedback Controller Design Phang Swee King A0033585A Email: king@nus.edu.sg NGS/ECE Dept. Faculty of Engineering National University
More informationMath 215/255: Elementary Differential Equations I Harish N Dixit, Department of Mathematics, UBC
Math 215/255: Elementary Differential Equations I Harish N Dixit, Department of Mathematics, UBC First Order Equations Linear Equations y + p(x)y = q(x) Write the equation in the standard form, Calculate
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationProperties of Transformations
6. - 6.4 Properties of Transformations P. Danziger Transformations from R n R m. General Transformations A general transformation maps vectors in R n to vectors in R m. We write T : R n R m to indicate
More informationMath 4B Notes. Written by Victoria Kala SH 6432u Office Hours: T 12:45 1:45pm Last updated 7/24/2016
Math 4B Notes Written by Victoria Kala vtkala@math.ucsb.edu SH 6432u Office Hours: T 2:45 :45pm Last updated 7/24/206 Classification of Differential Equations The order of a differential equation is the
More informationQuaternion Curvatures
Quaternion Curvatures Kurt Nalty September 7, 2011 Abstract Four dimensional curves can be parameterized along pathlength by three scalars using the Frenet-Serret formulas. Using quaternion notation, these
More informationFinal Exam Solutions
EE55: Linear Systems Final Exam SIST, ShanghaiTech Final Exam Solutions Course: Linear Systems Teacher: Prof. Boris Houska Duration: 85min YOUR NAME: (type in English letters) I Introduction This exam
More informationarxiv: v4 [physics.comp-ph] 21 Jan 2019
A spectral/hp element MHD solver Alexander V. Proskurin,Anatoly M. Sagalakov 2 Altai State Technical University, 65638, Russian Federation, Barnaul, Lenin prospect,46, k2@list.ru 2 Altai State University,
More informationLecture 4 Continuous time linear quadratic regulator
EE363 Winter 2008-09 Lecture 4 Continuous time linear quadratic regulator continuous-time LQR problem dynamic programming solution Hamiltonian system and two point boundary value problem infinite horizon
More informationReduced-Order Modeling of Channel Flow Using Traveling POD and Balanced POD
3rd AIAA Flow Control Conference, 5 8 June 26, San Francisco Reduced-Order Modeling of Channel Flow Using Traveling POD and Balanced POD M. Ilak and C. W. Rowley Dept. of Mechanical and Aerospace Engineering,
More informationLinear Algebra Review (Course Notes for Math 308H - Spring 2016)
Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,
More informationProblem Set 5 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationEL2520 Control Theory and Practice
EL2520 Control Theory and Practice Lecture 8: Linear quadratic control Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden Linear quadratic control Allows to compute the controller
More informationDynamics and control of wall-bounded shear flows
Dynamics and control of wall-bounded shear flows Mihailo Jovanović www.umn.edu/ mihailo Center for Turbulence Research, Stanford University; July 13, 2012 M. JOVANOVIĆ 1 Flow control technology: shear-stress
More informationCollocation approximation of the monodromy operator of periodic, linear DDEs
p. Collocation approximation of the monodromy operator of periodic, linear DDEs Ed Bueler 1, Victoria Averina 2, and Eric Butcher 3 13 July, 2004 SIAM Annual Meeting 2004, Portland 1=Dept. of Math. Sci.,
More informationExercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
More information9. Boundary layers. Flow around an arbitrarily-shaped bluff body. Inner flow (strong viscous effects produce vorticity) BL separates
9. Boundary layers Flow around an arbitrarily-shaped bluff body Inner flow (strong viscous effects produce vorticity) BL separates Wake region (vorticity, small viscosity) Boundary layer (BL) Outer flow
More informationState Regulator. Advanced Control. design of controllers using pole placement and LQ design rules
Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state
More informationIterative Learning Control Analysis and Design I
Iterative Learning Control Analysis and Design I Electronics and Computer Science University of Southampton Southampton, SO17 1BJ, UK etar@ecs.soton.ac.uk http://www.ecs.soton.ac.uk/ Contents Basics Representations
More informationFrequency domain analysis of the linearized Navier-Stokes equations
Frequency domain analysis of the linearized Navier-Stokes equations Mihailo R. JovanoviC jrnihailoqengineering.ucsb.edu Bassam Bamieh barnieh@engineering.ucsb.edu Department of Mechanical and Environmental
More informationLinear stability of channel entrance flow
Linear stability of channel entrance flow Damien Biau DIAM, University of Genova, Via Montallegro 1, 16145 Genova, Italy Abstract The spatial stability analysis of two dimensional, steady channel flow
More informationDirect Numerical Simulations of converging-diverging channel flow
Intro Numerical code Results Conclusion Direct Numerical Simulations of converging-diverging channel flow J.-P. Laval (1), M. Marquillie (1) Jean-Philippe.Laval@univ-lille1.fr (1) Laboratoire de Me canique
More informationEstimation of Turbulent Dissipation Rate Using 2D Data in Channel Flows
Proceedings of the 3 rd World Congress on Mechanical, Chemical, and Material Engineering (MCM'17) Rome, Italy June 8 10, 2017 Paper No. HTFF 140 ISSN: 2369-8136 DOI: 10.11159/htff17.140 Estimation of Turbulent
More informationControl Design. Lecture 9: State Feedback and Observers. Two Classes of Control Problems. State Feedback: Problem Formulation
Lecture 9: State Feedback and s [IFAC PB Ch 9] State Feedback s Disturbance Estimation & Integral Action Control Design Many factors to consider, for example: Attenuation of load disturbances Reduction
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More informationNon-homogeneous equations (Sect. 3.6).
Non-homogeneous equations (Sect. 3.6). We study: y + p(t) y + q(t) y = f (t). Method of variation of parameters. Using the method in an example. The proof of the variation of parameter method. Using the
More information