Optimal rotary control of the cylinder wake using POD reduced order model p. 1/30
|
|
- Jessie Richards
- 5 years ago
- Views:
Transcription
1 Optimal rotary control of the cylinder wake using POD reduced order model Michel Bergmann, Laurent Cordier & Jean-Pierre Brancher Laboratoire d Énergétique et de Mécanique Théorique et Appliquée UMR 7563 (CNRS - INPL - UHP) ENSEM - 2, avenue de la Forêt de Haye BP Vandoeuvre Cedex, France Optimal rotary control of the cylinder wake using POD reduced order model p. 1/30
2 Outline I - Flow configuration and numerical methods II - Optimal control III - Proper Orthogonal Decomposition (POD) IV - Reduced Order Model of the cylinder wake (ROM) V - Optimal control formulation applied to the ROM VI - Results of POD ROM VII - General observations VIII - Nelder Mead Simplex method Conclusions and perspectives Optimal rotary control of the cylinder wake using POD reduced order model p. 2/30
3 PSfrag replacements I - Configuration and numerical method Two dimensional flow around a circular cylinder at R e = 200 Viscous, incompressible and Newtonian fluid Cylinder oscillation with a tangential velocity γ(t) Γ e U y γ(t) D Γ c x Γ sup Γ s Γ inf Fractional steps method in time Finite Elements Method (FEM) in space Numerical code written by M.Braza (IMFT-EMT2) & D.Ruiz (ENSEEIHT) Optimal rotary control of the cylinder wake using POD reduced order model p. 3/30
4 I - Configuration and numerical method Iso pressure at t = 100. Iso vorticity at t = placements CD, CL C D C L time units Aerodynamics coefficients. Authors S t C D Braza et al. (1986) Henderson et al. (1997) He et al. (2000) this study Strouhal number and drag coefficient. Optimal rotary control of the cylinder wake using POD reduced order model p. 4/30
5 II - Optimal control Definition Mathematical method allowing to determine without empiricism a control law starting from the optimization of a cost functional. State equation F(φ, c) = 0 ; (Navier-Stokes + C.I. + C.L.) Control variables c ; (Blowing/suction, design parameters...) Cost functional J (φ, c). (Drag, lift,...) Find a control law c and state variable φ such that the cost functional J (φ, c) reach an extremum under the constrain F(φ, c) = 0. Optimal rotary control of the cylinder wake using POD reduced order model p. 5/30
6 II - Optimal control Lagrange multipliers Constrained optimization unconstrained optimization Introduction of Lagrange multipliers ξ. Lagrange functional : L(φ, c, ξ) = J (φ, c) < F(φ, c), ξ > Force L to be stationary look for δl = 0 : δl = L L L δφ + δc + φ c ξ δξ = 0 Suppose φ, c and ξ independant each other : L L L δφ = δc = φ c ξ δξ = 0 Optimal rotary control of the cylinder wake using POD reduced order model p. 6/30
7 II - Optimal control optimal system State equation ( L ξ δξ = 0) : F(φ, c) = 0 Co-sate equation ( L δφ = 0) : φ ( ) F ξ = φ Optimality condition ( L δc = 0) : c ( ) J = c ( J φ ( F c ) ) ξ Expensive method in CPU time and storage memory for large system! Ensure only a local (generally not global) minimum Optimal rotary control of the cylinder wake using POD reduced order model p. 7/30
8 II - Optimal control Reduced Order Model (ROM) "without an inexpensive method for reducing the cost of flow computation, it is unlikely that the solution of optimization problems involving the three dimensional unsteady Navier-Stokes system will become routine" Initialization M. Gunzburger, 2000 Recourse to detailed model (TRPOD) High fidelity model f(x), grad f(x) Approximation model a(x), grad a(x) x Optimization Optimization on simplified model Optimal rotary control of the cylinder wake using POD reduced order model p. 8/30
9 II - Proper Orthogonal Decomposition (POD) Introduced in fluid mechanics (turbulence context) by Lumley (1967). Look for a realization φ(x) which is closer, in an average sense, to the realizations u(x). (X = (x, t) D = Ω R + ) φ(x) solution of the problem : max φ (u, φ 2 s.t. φ 2 = 1. Snapshots method, Sirovich (1987) : C(t, t )a (n) (t ) dt = λ (n) a (n) (t). T Optimal convergence L 2 norm (energy) of φ(x) Dynamical order reduction is possible. Decomposition of the velocity field : u(x, t) = N P OD i=1 a (i) (t)φ (i) (x). Optimal rotary control of the cylinder wake using POD reduced order model p. 9/30
10 III - Reduced Order Model of the cylinder wake (ROM) Galerkin s projection of NSE on the POD basis : ( φ (i), u t ) + (u )u = (φ (i), p + 1Re ) u. Integration by parts (Green s formula) leads : ( φ (i), u t ) + (u )u = ( p, φ (i)) 1 Re [p φ (i) ] + 1 Re [( u)φ(i) ]. ( ) ( φ (i) ) T, u with [a] = Γ a n dγ and (A, B) = Ω A : B dω = i, j Ω A ij B ji dω. Optimal rotary control of the cylinder wake using POD reduced order model p. 10/30
11 III - Reduced Order Model of the cylinder wake (ROM) Velocity decomposition with N P OD modes : u(x, t) = u m (x) + γ(t) u c (x) + N P OD k=1 a (k) (t)φ (k) (x). Reduced order dynamical system where only N gal ( N P OD ) modes are retained (state s equation) : d a (i) (t) d t =A i + N gal j=1 + D i d γ d t + B ij a (j) (t) + E i + a (i) (0) = (u(x, 0), φ (i) (x)). N gal j=1 N gal j=1 N gal k=1 C ijk a (j) (t)a (k) (t) F ij a (j) (t) γ + G i γ 2 A i, B ij, C ijk, D i, E i, F ij and G i depend on φ, u m, u c and Re. Optimal rotary control of the cylinder wake using POD reduced order model p. 11/30
12 IV - Reduced Order Model of the cylinder wake Stabilization Integration and (optimal) stabilization of the POD ROM for γ = A sin(2πs t t), A = 2 et S t = 0.5. POD reconstruction errors temporal modes amplification ments a (n) time units Temporal evolution of the first 6 POD temporal modes. Causes : Extraction of large energetic eddies Dissipation takes place in small eddies Solution : Optimal addition of artificial viscosity on each POD mode projection (Navier Stokes) prediction before stabilisation (POD ROM) prediction after stabilisation (POD ROM). Optimal rotary control of the cylinder wake using POD reduced order model p. 12/30
13 IV - Reduced Order Model of the cylinder wake Stabilization s j d b a (n) 10-1 PSfrag replacements a (n) a (n) proj POD modes index Comparison of energetic spectrum POD modes index Comparison of absolute errors. Good agreements between POD and DNS spectrum Reduced reconstruction error between predicted (POD) and projected (DNS) modes Validation of the POD ROM Optimal rotary control of the cylinder wake using POD reduced order model p. 13/30
14 IV - Optimal control formulation based on reduced order model Objective functional : J (a, γ(t)) = T 0 J(a, γ(t)) dt = T 0 N gal a (i)2 + α 2 γ(t)2 dt. i=1 α : regularization parameter (penalization). Co-state s equations : d ξ (i) N (t) gal = B ji + γ F ji + dt j=1 ξ (i) (T ) = 0. N gal k=1 (C jik + C jki ) a (k) ξ (j) (t) 2a (i) Optimality condition (gradient) : N gal dξ (i) N gal δγ(t) = D i + E i + dt i=1 i=1 N gal j=1 F ij a (j) + 2G i γ(t) ξ (i) + αγ. Optimal rotary control of the cylinder wake using POD reduced order model p. 14/30
15 IV - Optimal control formulation based on reduced order model γ (0) (t) done ; for n = 0, 1, 2,... and while a convergence criterium is not satisfied, do : 1. From t = 0 to t = T solve the state s equations with γ (n) (t) ; state s variables a (n) (t) 2. From t = T to t = 0 solve the co-state s equations with a (n) (t) ; co-state s variables ξ (n) (t) 3. Solve the optimality condition with a (n) (t) and ξ (n) (t) ; objective gradient δγ (n) (t) 4. New control law γ (n+1) (t) = γ (n) (t) + ω (n) δγ (n) (t) End do. Optimal rotary control of the cylinder wake using POD reduced order model p. 15/30
16 V - Closed loop results Generalities No reactualization of the POD basis. The energetic representativity is a priori different to the dynamical one : possible inconvenient for control, a POD dynamical system represents a priori only the dynamics (and its vicinity) used to build the low dynamical model. Construction of a POD basis representative of a large range of dynamics : excitation of a great number of degrees of freedom scanning γ(t) in amplitudes and frequencies. Optimal rotary control of the cylinder wake using POD reduced order model p. 16/30
17 V - Closed loop results Excitation γe s -3-4 PSfrag replacements time units S t γ e (t) = A 1 sin(2πs t 1 t) sin(2πs t 2 t A 2 sin(2πs t 3 t)) with A 1 = 4, A 2 = 18, S t 1 = 1/120, S t 2 = 1/3 and S t 3 = 1/60. 0 amplitudes 4 and Fourier analysis 0 frequencies 0.65 γ e initial control law in the iterative process. Optimal rotary control of the cylinder wake using POD reduced order model p. 17/30
18 V - Closed loop results Energy s γ = γ = γ e 10-5 PSfrag replacements c λ (i) POD modes index λ (i) relative Ec POD modes index Stationary cylinder γ = 0 : 2 modes out of 100 are sufficient to restitute 97% of the kinetic energy. Controlled cylinder γ = γ e : 40 modes out of 100 are then necessary to restitute 97% of the kinetic energy Robustness evolution with dynamical evolutions. Optimal rotary control of the cylinder wake using POD reduced order model p. 18/30
19 V - Closed loop results Optimal control γopt s -2 PSfrag replacements time units S t Reduction of the wake instationarity. γ opt A sin(2πs t t) with A = 2.2 and S t = 0.53 J (γ e ) = 9.81 = J (γ opt ) = The control is optimal for the reduced order model based on POD. Is it also optimal for the Navier Stokes model? Optimal rotary control of the cylinder wake using POD reduced order model p. 19/30
20 V - Closed loop results Comparison of wakes organization No mathematical proof concerning the Navier Stokes optimality. no control γ = 0 Isocontours of vorticity ω z. optimal control γ = γ opt no control : γ = 0 Asymmetrical flow. Large and energetic eddies. optimal control : γ = γ opt Symmetrization of the (near) wake. Smaller and lower energetic eddies. Optimal rotary control of the cylinder wake using POD reduced order model p. 20/30
21 V - Closed loop results Aerodynamics coefficients CD CL s 1.1 PSfrag replacements time units time units Very consequent drag reduction : C D 0 = 1.40 for γ = 0 et C D = 1.04 for γ = γ opt C D /C D 0 = 0.74 more than 25%. Decrease of the lift amplitude : C L = 0.68 for γ = 0 et C L = 0.13 for γ = γ opt. Optimal rotary control of the cylinder wake using POD reduced order model p. 21/30
22 V - Closed loop results Numerical costs Optimal control of NSE by He et al. (2000) : harmonic control law with A = 3 and S t = % drag reduction. Optimal control POD ROM (this study) : harmonic control law with A = 2.2 and S t = % drag reduction. Less energetic costs (greater energetic gain?) Calculus time costs : 100 less using POD ROM than NSE! (for co-states equations and optimality conditions too) Memory storage : 600 less variables using POD ROM than NSE! "Optimal" control of 3D flows becomes possible! Optimal rotary control of the cylinder wake using POD reduced order model p. 22/30
23 VI - General observations Numerical experimentation A Observations Minimum is located in a smooth valley Global minimum : around A = 4.4 and St = 0.76 Maximum is located in a sharp pic St Iso-relative- drag coefficient C D (A, St)/C D 0 in (A, St) plan. Global maximum : near St = 0.2, the natural frequency : lock-on flow Finding the global minimum with an optimization algorithm may leads difficulties because of the smooth valley Optimal rotary control of the cylinder wake using POD reduced order model p. 23/30
24 VI - General observations Maximal rotation angle CD(θ)/CD St = 0.2 St = 0.3 St = 0.4 St = 0.5 St = 0.6 St = 0.7 St = 0.8 St = 0.9 St = 1.0 St = 0.53 maximum rotation angle : θ = A/(πSt) Observations No drag reduction possible near natural frequency Maximum drag reduction around θ = θ opt = θ Relative drag coefficient for different Strouhal numbers vs. maximum rotation angle. For all frequencies g.t. natural frequency Minimum drag : C Dmin = 0.71C D 0 = 0.98 Existence of an "optimal" maximal rotation angle. Optimal rotary control of the cylinder wake using POD reduced order model p. 24/30
25 VI - General observations Maximal rotation angle C D min C D θ opt Notations C D min(st) = min A R C D(θ, St) C D θopt(st) = C D (95, St) CD 0.9 C D Observations 0.7 St Nat St Comparison between C D min and C D θopt. Good agreements between C D min and C D θopt for St > St Nat θopt is not optimal for St < St Nat In order to minimise the drag one couldn t choose A and St independently. It seems that A/St = 5.2 (θ opt = 95 ). Optimal rotary control of the cylinder wake using POD reduced order model p. 25/30
26 VI - General observations Discussion Control law obtain by POD ROM is not optimal for drag minimisation Parameters obtain : A = 2.2 and St = 0.53 (θ max = 76 ) C D = 1.04 Optimal parameters : A = 4.4 and St = 0.76 (θ max = 105 θ opt ) C D = 0.98 Results in (A, St) quite different but not so far in C D terms The smooth valley is reached Couple an efficient new optimisation algorithm for smooth fonctions Take results obtain by POD ROM as initial conditions Optimal rotary control of the cylinder wake using POD reduced order model p. 26/30
27 VII - Nelder Mead Simplex method Generalities Avantages Numerical simplicities Adaptive topology Gradients calculations not necessary Good results with smooth fonctions Disadvantages No proof of optimality for simplex dimensions greater than two Need to fix free parameters Maybe more iterations than gradient based optimisation algorithms... Optimal rotary control of the cylinder wake using POD reduced order model p. 27/30
28 VII - Nelder Mead Simplex method Results A Observations Topology adaptation function of the curve of the valley Minimum found by Nelder-Mead simplex method : A = 4.5 and St = 0.76 θ = 107 Seems to be the global minimum St Iso-relative- drag coefficient C D (A, St)/C D 0 in (A, St) plan. 30 NSE resolutions 5% additive drag reduction compare to POD ROM Relative drag reduction by POD ROM : 25% (1 NSE resolution) Utility of coupling a new algorithm? Optimal rotary control of the cylinder wake using POD reduced order model p. 28/30
29 Conclusions Significative drag reduction minimizing the wake instationnarity of the reduced order model : more than 25% of the relative drag But this is not a minimum (global a local) of the drag function Low calculus cost numerical cost negligible. Calculus time costs : 100 less using POD ROM than NSE! (for co-states equations and optimality conditions too) Memory storage : 600 less variables using POD ROM than NSE! "OPTIMAL" CONTROL OF 3D FLOWS POSSIBLE BY POD ROM Existence of a optimal maximal rotation angle for effective drag reduction, θ = 105 Couple POP ROM to Nelder Mead simplex method leads a priori the global minimum of the drag functional But the gain on the drag function is quite small compare to result obtain by POD ROM Optimal rotary control of the cylinder wake using POD reduced order model p. 29/30
30 Perspectives Improve the representativity of the low order model "Optimize" the temporal excitation γ e Mix snapshots corresponding to several differents dynamics (temporal excitations) Look for harmonic control γ(t) = A sin(2π S t t) with POD basis reactualization (close loop on NSE and not only on POD ROM) Couple this optimal system with trust region methods (TRPOD) = proof of convergence under weak conditions Couple pressure with the POD dynamical system pressure contribution to drag coefficient : 80% Optimal control of the Navier Stokes equations Optimal rotary control of the cylinder wake using POD reduced order model p. 30/30
Optimal rotary control of the cylinder wake using POD Reduced Order Model. Abstract
Optimal rotary control of the cylinder wake using POD Reduced Order Model Michel Bergmann, Laurent Cordier, and Jean-Pierre Brancher LEMTA, UMR 7563 (CNRS - INPL - UHP) ENSEM - 2, avenue de la forêt de
More informationDrag Minimization of the Cylinder Wake by Trust-Region Proper Orthogonal Decomposition
Drag Minimization of the Cylinder Wake by Trust-Region Proper Orthogonal Decomposition Michel Bergmann, Laurent Cordier, Jean-Pierre Brancher To cite this version: Michel Bergmann, Laurent Cordier, Jean-Pierre
More informationOn the generation of a reverse Von Karman street for the controlled cylinder wake in the laminar regime
On the generation of a reverse Von Karman street for the controlled cylinder wake in the laminar regime Michel Bergmann, Laurent Cordier, Jean-Pierre Brancher To cite this version: Michel Bergmann, Laurent
More informationImprovement of Reduced Order Modeling based on Proper Orthogonal Decomposition
ICCFD5, Seoul, Korea, July 7-11, 28 p. 1 Improvement of Reduced Order Modeling based on Proper Orthogonal Decomposition Michel Bergmann, Charles-Henri Bruneau & Angelo Iollo Michel.Bergmann@inria.fr http://www.math.u-bordeaux.fr/
More informationOptimal control of the cylinder wake in the laminar regime by Trust-Region methods and POD Reduced-Order Models
Optimal control of the cylinder wake in the laminar regime by Trust-Region methods and POD Reduced-Order Models M Bergmann a,, L Cordier b a MAB, UMR - Université Bordeaux, 3 cours de la Libération, 33
More informationAccurate POD Reduced-Order Models of separated flows
Accurate POD Reduced-Order Models of separated flows Julien FAVIER DICAT, University of Genova, via Montallegro 1, 16145 Genova, Italy Laurent CORDIER LEA - CEAT, 43, rue de l Aérodrome, 86036 Poitiers
More informationControl of the circular cylinder wake by Trust-Region methods and POD Reduced-Order Models
INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE Control of the circular cylinder wake by Trust-Region methods and POD Reduced-Order Models inria-848, version - 9 Jun 8 Michel Bergmann
More informationImprovement of Reduced Order Modeling based on Proper Orthogonal Decomposition
INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE Improvement of Reduced Order Modeling based on Proper Orthogonal Decomposition Michel Bergmann Charles-Henri Bruneau Angelo Iollo inria-288415,
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationOverview of sparse system identification
Overview of sparse system identification J.-Ch. Loiseau 1 & Others 2, 3 1 Laboratoire DynFluid, Arts et Métiers ParisTech, France 2 LIMSI, Université d Orsay CNRS, France 3 University of Washington, Seattle,
More informationA fundamental study of the flow past a circular cylinder using Abaqus/CFD
A fundamental study of the flow past a circular cylinder using Abaqus/CFD Masami Sato, and Takaya Kobayashi Mechanical Design & Analysis Corporation Abstract: The latest release of Abaqus version 6.10
More informationCalibration of POD Reduced-Order Models using Tikhonov regularization
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 29; : [Version: 22/9/8 v.] Calibration of POD Reduced-Order Models using Tikhonov regularization L. Cordier,, B. Abou El
More informationStructures in the turbulent wake of an Ahmed body: LES vs Experiments using POD and some ideas about control
LES vs Experiments using POD and some ideas about control Stéphanie PELLERIN and Bérengère PODVIN Laboratoire d Informatique pour la Mécanique et les Sciences de l Ingénieur Context Université Paris Sud,
More informationProper Orthogonal Decomposition
Proper Orthogonal Decomposition Kameswararao Anupindi School of Mechanical Engineering Purdue University October 15, 2010 Kameswararao Anupindi (Purdue University) ME611, Principles of Turbulence October
More informationUNIVERSITY OF CALGARY. Base Region Topology of Turbulent Wake around Finite Wall-Mounted Cylinder with. Application of Low Order Flow Representation
UNIVERSITY OF CALGARY Base Region Topology of Turbulent Wake around Finite Wall-Mounted Cylinder with Application of Low Order Flow Representation by Golriz Boorboor A THESIS SUBMITTED TO THE FACULTY OF
More informationA Numerical Investigation of Velocity-Pressure Reduced Order Models for Incompressible Flows
A Numerical Investigation of Velocity-Pressure educed Order Models for Incompressible Flows Alfonso Caiazzo a, Traian Iliescu b,1, Volker John a,c,, Swetlana Schyschlowa a a Weierstrass Institute for Applied
More informationSimulation of free surface fluids in incompressible dynamique
Simulation of free surface fluids in incompressible dynamique Dena Kazerani INRIA Paris Supervised by Pascal Frey at Laboratoire Jacques-Louis Lions-UPMC Workshop on Numerical Modeling of Liquid-Vapor
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 informationChapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
More informationNumerical Optimization of Partial Differential Equations
Numerical Optimization of Partial Differential Equations Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada URL: http://www.math.mcmaster.ca/bprotas Rencontres
More informationVorticity and Dynamics
Vorticity and Dynamics In Navier-Stokes equation Nonlinear term ω u the Lamb vector is related to the nonlinear term u 2 (u ) u = + ω u 2 Sort of Coriolis force in a rotation frame Viscous term ν u = ν
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 informationSuboptimal Open-loop Control Using POD. Stefan Volkwein
Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,
More information(U c. t)/b (U t)/b
DYNAMICAL MODELING OF THE LARGE-SCALE MOTION OF A PLANAR TURBULENT JET USING POD MODES. S. Gordeyev 1 and F. O. Thomas 1 University of Notre Dame, Notre Dame, USA University of Notre Dame, Notre Dame,
More informationReliability of LES in complex applications
Reliability of LES in complex applications Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) DESIDER Symposium Corfu, June 7-8, 27 Sample of complex flow
More informationMestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 2012/13
Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 212/13 Exam 2ª época, 2 February 213 Name : Time : 8: Number: Duration : 3 hours 1 st Part : No textbooks/notes allowed 2 nd Part :
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationProper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems. Stefan Volkwein
Proper Orthogonal Decomposition (POD) for Nonlinear Dynamical Systems Institute for Mathematics and Scientific Computing, Austria DISC Summerschool 5 Outline of the talk POD and singular value decomposition
More informationApplications of the Proper Orthogonal Decomposition Method
Applications of the Proper Orthogonal Decomposition Method Internship report Katrin GRAU Mentoring: Laurent GICQUEL September 2007 WN/CFD/07/97 I would like to thank the whole CFD team at CERFACS, in particular
More informationPost-processing of solutions of incompressible Navier Stokes equations on rotating spheres
ANZIAM J. 50 (CTAC2008) pp.c90 C106, 2008 C90 Post-processing of solutions of incompressible Navier Stokes equations on rotating spheres M. Ganesh 1 Q. T. Le Gia 2 (Received 14 August 2008; revised 03
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 informationINTRODUCTION OBJECTIVES
INTRODUCTION The transport of particles in laminar and turbulent flows has numerous applications in engineering, biological and environmental systems. The deposition of aerosol particles in channels and
More informationVortex Induced Vibrations
Vortex Induced Vibrations By: Abhiroop Jayanthi Indian Institute of Technology, Delhi Some Questions! What is VIV? What are the details of a steady approach flow past a stationary cylinder? How and why
More informationAn approach to predict gust effects by means of hybrid ROM/CFD simulations
An approach to predict gust effects by means of hybrid ROM/CFD simulations M. Bergmann 1, A. Ferrero 1, A. Iollo 1, H. Telib 2 1 Inria Bordeaux Sud-Ouest and Université de Bordeaux, Talence, France 2 Optimad
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationReduced-order models for flow control: balanced models and Koopman modes
Reduced-order models for flow control: balanced models and Koopman modes Clarence W. Rowley, Igor Mezić, Shervin Bagheri, Philipp Schlatter, and Dan S. Henningson Abstract This paper addresses recent developments
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationModal Decomposition Methods on Aerodynamic Flows
498 Modal Decomposition Methods on Aerodynamic Flows Mohammad Niaz Murshed University of Michigan Department of Aerospace Engineering February 5, 2018 Abstract- This paper discusses the significance of
More informationSuppression of 3D flow instabilities in tightly packed tube bundles
Suppression of 3D flow instabilities in tightly packed tube bundles Nicholas Kevlahan kevlahan@mcmaster.ca Department of Mathematics & Statistics CSFD, June 13 15 2004 p.1/33 Collaborators CSFD, June 13
More informationLES Study of Synthetic Jet Frequency and Amplitude Effects on a Separated Flow
1 LES Study of Synthetic Jet Frequency and Amplitude Effects on a Separated Flow ONERA Applied Aerodynamic Department Pierre-Yves Pamart, Julien Dandois, Eric Garnier, Pierre Sagaut 2 Introduction Separated
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationLearning-based Robust Stabilization for Reduced-Order Models of 2D and 3D Boussinesq Equations
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Learning-based Robust Stabilization for Reduced-Order Models of 2D and 3D Boussinesq Equations Benosman, M.; Borggaard, J.; kramer, B. TR2017-084
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 informationIMPLEMENTATION OF POD AND DMD METHODS IN APACHE SPARK FRAMEWORK FOR SIMULATION OF UNSTEADY TURBULENT FLOW IN THE MODEL COMBUSTOR
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationLARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL
LARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL 09/11/2017 Paolo Casco Stephie Edwige Philippe Gilotte Iraj Mortazavi LES and flow control over a 25 ramp model : context 2 Context Validation
More informationOptimization and control of a separated boundary-layer flow
Optimization and control of a separated boundary-layer flow Journal: 2011 Hawaii Summer Conferences Manuscript ID: Draft lumeetingid: 2225 Date Submitted by the Author: n/a Contact Author: PASSAGGIA, Pierre-Yves
More informationReduced Order Drag Modeling of Liquid Drops
ILASS-Americas 24th Annual Conference on Liquid Atomization and Spray Systems, San Antonio, TX, May 2012 Reduced Order Drag Modeling of Liquid Drops B. T. Helenbrook Department of Mechanical and Aeronautical
More informationGlobal well-posedness and decay for the viscous surface wave problem without surface tension
Global well-posedness and decay for the viscous surface wave problem without surface tension Ian Tice (joint work with Yan Guo) Université Paris-Est Créteil Laboratoire d Analyse et de Mathématiques Appliquées
More informationnek5000 massively parallel spectral element simulations
nek5000 massively parallel spectral element simulations PRACE Scientific Seminar HPC Boosts Science, 22th February 2011 P. Schlatter & D. S. Henningson Linné Flow Centre, KTH Mechanics Fluid flows Tornado,
More informationNumerical Study of Natural Unsteadiness Using Wall-Distance-Free Turbulence Models
Numerical Study of Natural Unsteadiness Using Wall-Distance-Free urbulence Models Yi-Lung Yang* and Gwo-Lung Wang Department of Mechanical Engineering, Chung Hua University No. 707, Sec 2, Wufu Road, Hsin
More informationarxiv: v1 [physics.flu-dyn] 15 Jun 2016
Under consideration for publication in J. Fluid Mech. 1 arxiv:1606.04735v1 [physics.flu-dyn] 15 Jun 2016 A reduced-order model of three-dimensional unsteady flow in a cavity based on the resolvent operator
More informationREE Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics
REE 307 - Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics 1. Is the following flows physically possible, that is, satisfy the continuity equation? Substitute the expressions for
More informationTowards Reduced Order Modeling (ROM) for Gust Simulations
Towards Reduced Order Modeling (ROM) for Gust Simulations S. Görtz, M. Ripepi DLR, Institute of Aerodynamics and Flow Technology, Braunschweig, Germany Deutscher Luft und Raumfahrtkongress 2017 5. 7. September
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationNumerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder
Numerical Simulation of Unsteady Flow with Vortex Shedding Around Circular Cylinder Ali Kianifar, Edris Yousefi Rad Abstract In many applications the flow that past bluff bodies have frequency nature (oscillated)
More informationBasic concepts in viscous flow
Élisabeth Guazzelli and Jeffrey F. Morris with illustrations by Sylvie Pic Adapted from Chapter 1 of Cambridge Texts in Applied Mathematics 1 The fluid dynamic equations Navier-Stokes equations Dimensionless
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationAccelerating incompressible fluid flow simulations on hybrid CPU/GPU systems
Accelerating incompressible fluid flow simulations on hybrid CPU/GPU systems Yushan Wang 1, Marc Baboulin 1,2, Karl Rupp 3,4, Yann Fraigneau 1,5, Olivier Le Maître 1,5 1 Université Paris-Sud, France 2
More informationApplication of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Matthew McMullen and Antony Jameson and Juan J. Alonso Dept. of Aeronautics & Astronautics Stanford University
More informationReduced Order Modeling of Steady and Unsteady Flow over a Sphere
ILASS-Americas 22nd Annual Conference on Liquid Atomization and Spray Systems, Cincinnati, OH, May 21 Reduced Order Modeling of Steady and Unsteady Flow over a Sphere B. T. Helenbrook and D. R. Witman
More informationRECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD
RECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD N. Jarrin 1, A. Revell 1, R. Prosser 1 and D. Laurence 1,2 1 School of MACE, the University of Manchester,
More informationNon-Newtonian Fluids and Finite Elements
Non-Newtonian Fluids and Finite Elements Janice Giudice Oxford University Computing Laboratory Keble College Talk Outline Motivating Industrial Process Multiple Extrusion of Pastes Governing Equations
More informationModeling, Control And Evaluation of Fluid Flow Systems Using Adaptation Based Linear Parameter Varying Models
Modeling, Control And Evaluation of Fluid Flow Systems Using Adaptation Based Linear Parameter Varying Models Coşku Kasnakoğlu Abstract In this paper a systematic method is proposed for the modeling and
More informationTURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS
HEFAT2014 10 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 16 July 2014 Orlando, Florida TURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS Everts, M.,
More informationApplication of a Virtual-Boundary Method for the Numerical Study of Oscillations Developing Behind a Cylinder Near A Plane Wall
Fluid Dynamics, Vol. 39, No. 1, 2004, pp. 61 68. Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, 2004, pp. 69 77. Original Russian Text Copyright 2004 by Kit, Nikitin,
More informationFinal abstract for ONERA Taylor-Green DG participation
1st International Workshop On High-Order CFD Methods January 7-8, 2012 at the 50th AIAA Aerospace Sciences Meeting, Nashville, Tennessee Final abstract for ONERA Taylor-Green DG participation JB Chapelier,
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationLevel Set-based Topology Optimization Method for Viscous Flow Using Lattice Boltzmann Method
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level Set-based Topology Optimization Method for Viscous Flow Using Lattice Boltzmann Method
More information12.1 Viscous potential flow (VPF)
1 Energy equation for irrotational theories of gas-liquid flow:: viscous potential flow (VPF), viscous potential flow with pressure correction (VCVPF), dissipation method (DM) 1.1 Viscous potential flow
More informationStability of flow past a confined cylinder
Stability of flow past a confined cylinder Andrew Cliffe and Simon Tavener University of Nottingham and Colorado State University Stability of flow past a confined cylinder p. 1/60 Flow past a cylinder
More informationARTICLE IN PRESS Physica D ( )
Physica D ( ) Contents lists available at ScienceDirect Physica D journal homepage: www.elsevier.com/locate/physd Robust model identification of actuated vortex wakes J. Weller, E. Lombardi, A. Iollo Institut
More informationProper Orthogonal Decomposition (POD)
Intro Results Problem etras Proper Orthogonal Decomposition () Advisor: Dr. Sorensen CAAM699 Department of Computational and Applied Mathematics Rice University September 5, 28 Outline Intro Results Problem
More informationVortex wake and energy transitions of an oscillating cylinder at low Reynolds number
ANZIAM J. 46 (E) ppc181 C195, 2005 C181 Vortex wake and energy transitions of an oscillating cylinder at low Reynolds number B. Stewart J. Leontini K. Hourigan M. C. Thompson (Received 25 October 2004,
More informationAPPLICATION OF PROPER ORTHOGONAL DECOMPOSITION FOR AVIATION SAFETY
APPLICATIO OF PROPER ORTHOGOAL DECOMPOSITIO FOR AVIATIO SAFETY Junho Cho*, Shinkyu Jeong***, Kwanjoong Yee* *Seoul ational University, **Kyung Hee University Keywords: Wake turbulence, POD, LES, Artifical
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationArchimer
Please note that this is an author-produced PDF of an article accepted for publication following peer review. The definitive publisher-authenticated version is available on the publisher Web site Comptes
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationInitial-value problem for shear flows
Initial-value problem for shear flows 1 Mathematical Framework The early transient and long asymptotic behaviour is studied using the initialvalue problem formulation for two typical shear flows, the plane
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 informationObtaining a Stable Galerkin ROM in Presence of Shock-Vortex Interactions
AIAA SciTech Forum 9-13 January 2017, Grapevine, Texas 55th AIAA Aerospace Sciences Meeting AIAA 2017-1008 Obtaining a Stable Galerkin ROM in Presence of Shock-Vortex Interactions Elnaz Rezaian and Mingjun
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 946 8633 A numerical investigation of velocity-pressure reduced order models
More informationContinuum Mechanics Lecture 5 Ideal fluids
Continuum Mechanics Lecture 5 Ideal fluids Prof. http://www.itp.uzh.ch/~teyssier Outline - Helmholtz decomposition - Divergence and curl theorem - Kelvin s circulation theorem - The vorticity equation
More informationLeast Squares Finite Element Methods for Large Scale Incompressible Flows
Least Squares Finite Element Methods for Large Scale Incompressible Flows by Tate T. H. Tsang Department of Chemical & Materials Engineering University of Kentucky Lexington, KY 40506 tsang@engr.uky.edu
More informationAccurate model reduction of transient and forced wakes
European Journal of Mechanics B/Fluids 26 (27) 354 366 Accurate model reduction of transient and forced wakes Bernardo Galletti a, Alessandro Bottaro b,charles-henribruneau c, Angelo Iollo c, a DIASP Politecnico
More informationNumerical study of the effects of trailing-edge bluntness on highly turbulent hydro-foil flows
Numerical study of the effects of trailing-edge bluntness on highly turbulent hydro-foil flows T. Do L. Chen J. Tu B. Anderson 7 November 2005 Abstract Flow-induced noise from fully submerged lifting bodies
More informationChapter 1: Basic Concepts
What is a fluid? A fluid is a substance in the gaseous or liquid form Distinction between solid and fluid? Solid: can resist an applied shear by deforming. Stress is proportional to strain Fluid: deforms
More informationCFD DESIGN OF A GENERIC CONTROLLER FOR VORTEX-INDUCED RESONANCE
Seventh International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 9-11 December 2009 CFD DESIGN OF A GENERIC CONTROLLER FOR VORTEX-INDUCED RESONANCE Andrew A. ANTIOHOS,
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationCHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1.
1/3 η 1C 2C, axi 1/6 2C y + < 1 axi, ξ > 0 y + 7 axi, ξ < 0 log-law region iso ξ -1/6 0 1/6 1/3 Figure 11.1: The Lumley triangle on the plane of the invariants ξ and η of the Reynolds-stress anisotropy
More information