Optimization and control of a separated boundary-layer flow
|
|
- Steven Stevens
- 6 years ago
- Views:
Transcription
1 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
2 Page 1 of Hawaii Summer Conferences Optimization and control of a separated boundary-layer flow P.-Y. Passaggia and U. Ehrenstein Institut de Recherche sur les Phénomènes Hors Equilibre (IRPHE), CNRS / Aix-Marseille Université, F Marseille, France In order to investigate the possibilities of controlling wall-bounded flows, optimal control of the nonlinear dynamics of a separated boundary-layer flow is studied using a localized blowing-suction at the wall. Control laws are computed using a formulation based on the augmented Lagrangian approach using the full knowledge of the flow dynamics. The case study consists of a laminar separated boundary-layer flow over a bump. I. Introduction Flow separation is ubiquitous to wall bounded flows when subject to an adverse pressure gradient or geometrical devices. The resulting recirculation bubble are source of instability phenomena leading to a loss of aerodynamic performances. It is essential to propose control solutions capable of minimizing there impact on the entire flow. In this paper, an adjoint-based optimization technique is presented in order to control the flow behind an obstacle. The feedback control of such dynamics has been recently investigated with success in the case of a flow over a shallow cavity 2 using model reduction based on the linearized perturbation dynamics. 6 However, model reduction based on the linear dynamics is likely to be not reliable as soon as absolute/convective instabilities coexist in the nonlinear regime. In order to assess the possibility of feedback control of such flow dynamics, the augmented Lagrangian approach has been implemented. This formulation takes into account the entire flow dynamics and was for instance shown to relaminarize a turbulent channel flow. 4 This formulation can also be adapted to the computation of optimal perturbations which are known to play an important role in the transient dynamics of open flows. In the present investigation, the flow is a laminar boundary layer where flow separation is triggered by a bump. This configuration produces a long recirculation bubble which becomes unstable beyond a critical Reynolds number. 5 This flow has been studied numerically by Marquillie & Ehrenstein 7 and their algorithm has been adapted to the present study. For the optimization problem with constraints, the gradient method based on a formulation of the adjoint Navier-Stokes system has been implemented. The control problem of the fully developed nonlinear dynamics, detailed in section 2, is formulated using a cost function based on either the integral or the terminal time kinetic energy of the perturbation. The effect of nonlinearities on the optimal initial disturbance and the control of the fully developed dynamics is presented in section 3. Finally the results are discussed in section 4. II. Formulation and numerical procedure The flow is decomposed into a baseflow and a perturbation and the Navier-Stokes system is written t u + (U )u + (u )U + (u )u + p (1/Re) u = 0, u = 0 (1) where U is the velocity of the baseflow and u, p are the velocity and pressure pertubation, respectively. The Reynolds number PhD Student, IRPHE, 49 rue F. Joliot Curie, B.P. 146, F Marseille Cedex 13, France Professor, IRPHE, Aix-Marseille Université 1 of 6
3 Page 2 of 6 Re = δ U (2) ν is based on the displacement thickness δ of the Blasius profile imposed at the entrance of the domain. Here, the numerical procedure of Marquillie & Ehrenstein 7 has been adapted for the optimization formulation. The bump is taken into account by use of a mapping which transforms the physical domain into a cartesian one. The streamwise direction x is discretized using finite differences whereas y, the wall-normal direction, is discretized using Chebyshev-collocation. The time integration is performed using a second order semiimplicit scheme. The difference between the present formulation and the one described in Marquillie & Ehrenstein comes from the algorithm used to solve the pressure. In the optimization formulation, the choice of outflow boundary condition is essential. In order to ensure the compatibility between the Navier-Stokes system and its adjoint, stress-free outflow boundary condition are to be implemented which is achieved using the influence matrix technique which couples pressure and velocity 8 II.A. Unstable baseflow At a supercritical Reynolds number of 650 used in the presnet investigation, a large domain with length L = 250 and a height H = 100 has been considered to minimize finite boundary effects. The first step of the analysis consists in generating a baseflow and the method of Åkervik et al. 1 ( called selective frequency damping) has been implemented in order to converge towards a stationary solution of the Navier-Stokes system. The method consists in solving the coupled system { q = f(q, Re) χ(q q) q = (q q) (3) Λ where f represents the Navier-Stokes system, q contains the velocity U and the pressure P of the baseflow. Here Λ is a cutoff frequency and χ is a damping coefficient. In the present case, the value of Λ corresponds to the dominant frequency measured using a probe localized inside the shear layer. A good estimation of the parameter χ is given by the growth rate of the linear instability, provided by the global stability analysis of the flow. 5 For Re = 650, the parameters used are = 6.37, χ=0.03 and figure 1 shows the resulting baseflow. Figure 1. Baseflow for Re = 650 II.B. Augmented Lagrangian and gradient In the case of the control, the aim is to minimize the kinetic energy in the entire domain and the cost function is either J 1,int (φ, u) = 1 2 when energy integral is considered or J 1,term (φ, u) = 1 2 T1 Ω o u u dx dt + γ Ω o u(t 1 ) u(t 1 ) dx dt + γ T1 T1 Γ c gφ gφ ds dt (4) Γ c gφ gφ ds dt (5) for energy optimization at time T 1. The portion of the boundary where the control is applied is denoted by Γ c, whereas g(x) defines the spatial distribution (taken as uniform in the present case) of the blowing-suction. The control has been applied at three different locations (each one of length 5) in the vicinity of the bump (cf. figure 2). The function φ(t) is the control law and γ refers to the cost of the control. The optimization 2 of 6
4 Page 3 of Hawaii Summer Conferences time interval is [, T 1 ]. The aim is to minimize the functional J 1, u being the solution of the Navier-Stokes system written as f(u, u, p) = 0, g(u, φ) = 0 with g(u, φ) = u Γc gφ. (6) In the case of the optimal perturbation analysis, the objective function becomes J 2 (φ, u) = 1 u(t 1 ) u(t 1 ) dx (7) 2 Ω the unknown being now the initial condition u 0, that is u is now solution of f(u, u, p) = 0, g(u, u 0 ) = u(x, 0) u 0 = 0, h(u 0, ɛ) = u 0 u 0 dx ɛ = 0 (8) where ɛ is the amplitude of the initial condition u 0. Lagrange multipliers u +, p +, φ + are introduced and in Ω Figure 2. Sketch of the domain and of the control setup order to solve the optimization problem, the Lagrangian L L(u, u +, p, p +, φ, φ +, ɛ, ɛ + ) = J 1,int (φ, u) < f(u, u, p), (u +, p + ) > < g(u, φ), φ + > (9) is to be rendered stationary. Cancelling the Fréchet derivatives with respect to the Lagrange multipliers consist in imposing the Navier-Stokes constraints. The derivatives with respect to the state variables (u, p) writes T1 T1 D u,p L (û, ˆp) = u û dxdt D u,p < f(u, p), (u +, p + ) > (û, ˆp) û gφ + ds dt (10) Ω Γ c The expression of D u,p < f(u, p), (u +, p + ) > (û, ˆp) is detailed for example in Bewley et al. 3 the above expression is equivalent to solve the adjoint system and canceling t u + (U )u + + ( U) t u + (u )u + + ( u) t u + + p + (1/Re) u + = u, u + = 0. (11) The right hand side of the adjoint momentum equation is zero in the case of a terminal time optimization. The adjoint system is integrated backward in time from T 1 to, starting with u + (T 1 ) = 0 in the intgral time formulation whereas for terminal time u + (T 1 ) = u(t 1 ). The integrations by part necessary to obtain the adjoint system introduce the integral on the boundary Ω T1 Ω [ 1 Re u+ ( û) n 1 Re ( u+ n)û (U n)u + û (u n)u + û + (p + n)û u + (ˆp n)] ds dt. (12) The stress-free conditions 1 Re u n pn = 0 being imposed at the outflow of the domain Ω S, the outflow condition for the adjoint system (implemented using the influence matrix) 1 Re ( u+ )n ΩS (U n)u + ΩS (u n)u + ΩS + p + n ΩS = 0, u + Ω ΩS = 0 (13) cancels this integral on Ω Γ c. Finally, the contribution on Γ c is cancelled taking gφ + = 1 Re u+ n + p + n (14) 3 of 6
5 Page 4 of 6 in (10). The derivative of the Lagrangian with respect to the control variable φ becomes D φ L ˆφ T1 = γ gφ g ˆφ T1 ds dt ( g ˆφ) gφ + ds dt (15) Γ c Γ c and taking into account (14), the gradient for the control is φ J 1 (φ) = γφ g g ds + ( 1 Γ c Γ c Re u+ n + p + n) g ds. (16) The same procedure is used for the optimal perturbation, the expression of the gradient now being { u0 J 2 (u 0 ) = u ɛ+ u 0 with u + 0 = u+ (x, ) ɛ + = Ω u+ ( ) u + ( ) dx/4ɛ (17) An iterative gradient algorithm φ k+1 = φ k α k φ J (φ k, u +,k ). is used to converge towards the optimal solution. A main difficulty is to appropriately choose the step α k and this parameter is optimized using a classical line search procedure. Once the optimal control is computed over the first time interval, the time marching is insured using a recieding-predictive control algorithm described in 4 with a time shift T a = T/2 equal to half the time of optimization in our case (see figure 3). Figure 3. Schematic of the recieding-predictive algorithm III. Optimal growth and control The linear optimal perturbation analysis 5 showed that perturbations localized near the separation point are responsible for the largest possible transient growth of the system. In order to examine the effect of Figure 4. a) Evolution of the structure of the optimal perturbation as a function of the initial amplitude, b) transient growth of optimal perturbations as a function of initial amplitude for T 1 = of 6
6 Page 5 of Hawaii Summer Conferences nonlinearities on these optimal perturbations, the initial amplitude of the disturbance is first set to small values in order to maintain the flow close to the linear regime. Then the amplitude of the initial disturbance is increased from 10 8 to 10 4 and successive optimization times were considered between T1 = 50 and T1 = 150. Figure 4a shows the shape of the optimal initial disturbance for three amplitudes considered and its structure is seen to spread in the shear layer, as well as upstream the bump when the ampliude is increased. This clearly shows that the nonlinear flow is not only sensitive to actuators located on the separation point, but also to regions located inside the recirculation bubble as well as upstream the bump. Figure 4b shows the energy evolution for the time horizon T1 = 150. The flow perturbation starts to saturate at T = 120 when the highest initial perturbation amplitude is considered. The possibilities of controlling the fully developed dynamics is then investigated. This unstable dynamics is characterized by a low frequency oscillations of the recirculation bubble5. The uncontrolled flow energy is shown as the solid line in figure 5 and the flapping frequency of Λ 200 is clearly visible. Figure 5. a) Energy of the perturbation for the uncontrolled case ( ), control using the integral time formulation with a single actuator (...), control with 3 actuators (.. ), control with 3 actuators ( ) using the terminal time formulation and b) associated control laws φ for the integral time formulation. Ideally one would like to take arbitrarily long time horizons for optimization. However, the direct solution in time u is input of the adjoint system and memory constraints put a bound on the maximum time T1 which can be considered. Here we have chosen T1 = 300 which is larger than the flapping period. The optimization procedure implies 5 to 7 computations of the gradient whereas 4 to 7 evaluations of the objective function are necessary to evaluate the optimal step α. A first actuator is introduced at the separation point shown as 1 in figure 2, where the nonlinear optimal perturbation dominates. The problem is then solved in time iterating Figure 6. Time-averaged streamlines for the uncontrolled flow a), and for the controlled flow b) using the integral formulation and 3 actuators. back and forth the direct and the adjoint system to achieve convergence for each time interval. Using this single actuator, the control is capable of diminishing the energy to half of the uncontrolled value without however achieving definite control, as seen in figure 5. In the attempt to improve the control, two additional actuators shown as 2 and 3 in figure 2 have been added. The corresponding energy evolutions are depicted in figure 5: although the energy is decreased further for T > 500 the perturbation energy cannot be reduced to zero. As shown in Marquillie & Ehrenstein,7 in the perturbed nonlinear regime, the mean recirculation 5 of 6
7 Page 6 of 6 Figure 7. Instantenous snapshot of the vorticity at time T = 600 for the uncontrolled flow a), for the controlled flow b) using the integral formulation and 3 actuators. length diminishes when compared with that of the baseflow. Therefore the recirculation length is likely to provide informations about the quality of the control. Figure 6 shows an increase from approximately 48.5 for the uncontrolled case to 63.4 when the control is applied. Note that the mean values have been computed for the interval 450 t 900. A snapshot of the vorticity at time t = 600 is depicted in figure 7. The use of terminal time formulation in the objective function has been advocated by ref 3, which in the case of the turbulent channel flow was capable to achieve a relaminarization. Here we imposed finite time energy optimization starting at time 350, however for the present flow dynamics the results are comparable with the integral time formulation, as can be seen in figure5. IV. Discussion The optimal growth and control of a separated laminar boundary layer flow over a bump is investigated numerically using an adjoint-based augmented Lagrangian approach. In this configuration, stress-free outflow boundary conditions compatible between the Navier-Stokes system and its adjoint were implemented, using the influence-matrix approach. The problem of optimal perturbations in the nonlinear regime has been studied, highlighting the regions of the flow which are the most sensitive to actuator devices. It is shown that in the nonlinear regime the flow dynamics becomes sensitive to perturbations aligned along the whole shearlayer, in contrast to the linear regime which is only sensitive to upstream perturbations. The control problem using blowing-suction at the wall has been addressed as well, using a recieding-predictive algorithm for time marching. The unstable dynamics related to the flapping of the recirculation bubble is clearly attenuated but a full control of the nonlinear dynamics could not be achieved and this point is still under investigations. For instance to take larger time horizons and/or multiple actuators is likely to improve the control capabilities. Acknowledgements The authors gratefully acknowledge the financial support of the Agence Nationale pour la Recherche (ref. ANR-09-SYSC-011). References 1 E. Åkervik, L. Brandt, D. Henningson, J. Hœpffner, O. Marxen, and P. Schlatter. Steady solutions of the navier-stokes equations by selective frequency damping. Phys. of Fluids, 18:068102, E. Åkervik, J. Hœpffner, U. Ehrenstein, and D. Henningson. Optimal growth, model reduction and control in a separated boundary-layer flow using global eigenmodes. J. Fluid Mech., 579: , T.R. Bewley. Flow control: new challenges for a new renaissance. Prog. in Aero. Sci., 37:21 58, T.R. Bewley, P. Moin, and R. Temam. Dns-based predictive control of turbulence: an optimal benchmark of feedback algortihms. J. Fluid Mech., 447: , U. Ehrenstein and F. Gallaire. Two-dimensional global low-frequency oscillations in a separating boundary-layer flow. J. Fluid Mech., 614: , J. Kim and T. R. Bewley. A linear systems approach fo flow control. Annu. Rev. Fluid Mech., 39: , M. Marquillie and U. Ehrenstein. On the onset of nonlinear oscillations in a separating boundary-layer flow. J. Fluid Mech., 490: , R. Peyret. Spectral Methods for Incompressible Flows. Springer, of 6
Ré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 informationLOW SPEED STREAKS INSTABILITY OF TURBULENT BOUNDARY LAYER FLOWS WITH ADVERSE PRESSURE GRADIENT
LOW SPEED STREAKS INSTABILITY OF TURBULENT BOUNDARY LAYER FLOWS WITH ADVERSE PRESSURE GRADIENT J.-P. Laval (1) CNRS, UMR 8107, F-59650 Villeneuve d Ascq, France (2) Univ Lille Nord de France, F-59000 Lille,
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationLow-speed streak instability in near wall turbulence with adverse pressure gradient
Journal of Physics: Conference Series Low-speed streak instability in near wall turbulence with adverse pressure gradient To cite this article: U Ehrenstein et al 2011 J. Phys.: Conf. Ser. 318 032027 View
More informationGeneral introduction to Hydrodynamic Instabilities
KTH ROYAL INSTITUTE OF TECHNOLOGY General introduction to Hydrodynamic Instabilities L. Brandt & J.-Ch. Loiseau KTH Mechanics, November 2015 Luca Brandt Professor at KTH Mechanics Email: luca@mech.kth.se
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 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 informationUNSTEADY DISTURBANCE GENERATION AND AMPLIFICATION IN THE BOUNDARY-LAYER FLOW BEHIND A MEDIUM-SIZED ROUGHNESS ELEMENT
UNSTEADY DISTURBANCE GENERATION AND AMPLIFICATION IN THE BOUNDARY-LAYER FLOW BEHIND A MEDIUM-SIZED ROUGHNESS ELEMENT Ulrich Rist and Anke Jäger Institut für Aerodynamik und Gasdynamik, Universität Stuttgart,
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 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 informationLagrangian acceleration in confined 2d turbulent flow
Lagrangian acceleration in confined 2d turbulent flow Kai Schneider 1 1 Benjamin Kadoch, Wouter Bos & Marie Farge 3 1 CMI, Université Aix-Marseille, France 2 LMFA, Ecole Centrale, Lyon, France 3 LMD, Ecole
More informationApplications of parabolized stability equation for predicting transition position in boundary layers
Appl. Math. Mech. -Engl. Ed., 33(6), 679 686 (2012) DOI 10.1007/s10483-012-1579-7 c Shanghai University and Springer-Verlag Berlin Heidelberg 2012 Applied Mathematics and Mechanics (English Edition) Applications
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 informationConvective instability and transient growth in flow over a backwardfacing
Convective instability and transient growth in flow over a backwardfacing step Journal: Manuscript ID: mss type: Date Submitted by the Author: Complete List of Authors: Keyword: Journal of Fluid Mechanics
More informationDirect numerical simulation and control of the flow over a backward facing ramp
Direct numerical control of the flow 2 1 IMB Université de Bordeaux, France 2 IMB Institut polytechnique de Bordeaux,France Inria Bordeaux Sud Ouest MC2 project 3rd GDR Contrôle des décollements Symposium
More informationWALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION
WALL PRESSURE FLUCTUATIONS IN A TURBULENT BOUNDARY LAYER AFTER BLOWING OR SUCTION Joongnyon Kim, Kyoungyoun Kim, Hyung Jin Sung Department of Mechanical Engineering, Korea Advanced Institute of Science
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 informationContribution of Reynolds stress distribution to the skin friction in wall-bounded flows
Published in Phys. Fluids 14, L73-L76 (22). Contribution of Reynolds stress distribution to the skin friction in wall-bounded flows Koji Fukagata, Kaoru Iwamoto, and Nobuhide Kasagi Department of Mechanical
More informationFluid Mechanics. Chapter 9 Surface Resistance. Dr. Amer Khalil Ababneh
Fluid Mechanics Chapter 9 Surface Resistance Dr. Amer Khalil Ababneh Wind tunnel used for testing flow over models. Introduction Resistances exerted by surfaces are a result of viscous stresses which create
More informationControl of a separating boundary layer with travelling waves on the wall
Control of a separating boundary layer with travelling waves on the wall By Antonios Monokrousos & Luca Brandt Linné Flow Centre, Department of Mechanics Royal Institute of Technology, SE-1 44 Stockholm,
More informationOptimal 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 informationAn example of the Rvachev function method
arxiv:1603.00320v1 [physics.flu-dyn] 1 Mar 2016 An example of the Rvachev function method Alexander V. Proskurin Altai State University, Altai State Technical University, k210@list.ru Anatoly M. Sagalakov
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 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 informationInstability of streaks in wall turbulence with adverse pressure gradient
Under consideration for publication in J. Fluid Mech. 1 Instability of streaks in wall turbulence with adverse pressure gradient MATTHIEU MARQUILLIE 1,2, UWE EHRENSTEIN 3 and JEAN-PHILIPPE LAVAL 1,2 3
More informationDIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING
DIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING Yukinori Kametani Department of mechanical engineering Keio
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 informationCoriolis effects on the elliptical instability in cylindrical and spherical rotating containers
18 ème Congrès Français de Mécanique Grenoble, 7-31 août 007 Coriolis effects on the elliptical instability in cylindrical and spherical rotating containers Michael Le Bars, Stéphane Le Dizès & Patrice
More informationEXCITATION OF GÖRTLER-INSTABILITY MODES IN CONCAVE-WALL BOUNDARY LAYER BY LONGITUDINAL FREESTREAM VORTICES
ICMAR 2014 EXCITATION OF GÖRTLER-INSTABILITY MODES IN CONCAVE-WALL BOUNDARY LAYER BY LONGITUDINAL FREESTREAM VORTICES Introduction A.V. Ivanov, Y.S. Kachanov, D.A. Mischenko Khristianovich Institute of
More informationWake behaviour and instability of flow through a partially blocked channel
J. Fluid Mech. (2007), vol. 582, pp. 319 340. c 2007 Cambridge University Press doi:10.1017/s0022112007005861 Printed in the United Kingdom 319 Wake behaviour and instability of flow through a partially
More informationValidation of an adjoint method for compressible channel flow sensitivities
9 ème Congrès Français de Mécanique Marseille, -8 août 9 Validation of an method for compressible channel flow sensitivities L. MORE-GABARRO, P. CAHALIFAUD, C. AIRIAU Université de oulouse ; CNRS, INP,
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 informationINVESTIGATION OF 2D AND 3D BOUNDARY-LAYER DISTURBANCES FOR ACTIVE CONTROL OF LAMINAR SEPARATION BUBBLES
INVESTIGATION OF 2D AND 3D BOUNDARY-LAYER DISTURBANCES FOR ACTIVE CONTROL OF LAMINAR SEPARATION BUBBLES Kai Augustin, Ulrich Rist and Siegfried Wagner Institut für Aerodynamik und Gasdynamik, Universität
More informationAerodynamics. High-Lift Devices
High-Lift Devices Devices to increase the lift coefficient by geometry changes (camber and/or chord) and/or boundary-layer control (avoid flow separation - Flaps, trailing edge devices - Slats, leading
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 informationConvective instability and transient growth in steady and pulsatile stenotic flows
J. Fluid Mech. (8), vol. 7, pp. 7 77. c 8 Cambridge University Press doi:1.117/s1181717 Printed in the United Kingdom 7 Convective instability and transient growth in steady and pulsatile stenotic flows
More informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
More informationThe Turbulent Rotational Phase Separator
The Turbulent Rotational Phase Separator J.G.M. Kuerten and B.P.M. van Esch Dept. of Mechanical Engineering, Technische Universiteit Eindhoven, The Netherlands j.g.m.kuerten@tue.nl Summary. The Rotational
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 informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
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 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 informationDirect Numerical Simulation of Jet Actuators for Boundary Layer Control
Direct Numerical Simulation of Jet Actuators for Boundary Layer Control Björn Selent and Ulrich Rist Universität Stuttgart, Institut für Aero- & Gasdynamik, Pfaffenwaldring 21, 70569 Stuttgart, Germany,
More informationExplicit algebraic Reynolds stress models for internal flows
5. Double Circular Arc (DCA) cascade blade flow, problem statement The second test case deals with a DCA compressor cascade, which is considered a severe challenge for the CFD codes, due to the presence
More informationDestabilizing turbulence in pipe flow
Destabilizing turbulence in pipe flow Jakob Kühnen 1*, Baofang Song 1,2*, Davide Scarselli 1, Nazmi Burak Budanur 1, Michael Riedl 1, Ashley Willis 3, Marc Avila 2 and Björn Hof 1 1 Nonlinear Dynamics
More informationHEAT TRANSFER IN A RECIRCULATION ZONE AT STEADY-STATE AND OSCILLATING CONDITIONS - THE BACK FACING STEP TEST CASE
HEAT TRANSFER IN A RECIRCULATION ZONE AT STEADY-STATE AND OSCILLATING CONDITIONS - THE BACK FACING STEP TEST CASE A.K. Pozarlik 1, D. Panara, J.B.W. Kok 1, T.H. van der Meer 1 1 Laboratory of Thermal Engineering,
More informationDirect numerical simulation of self-similar turbulent boundary layers in adverse pressure gradients
Direct numerical simulation of self-similar turbulent boundary layers in adverse pressure gradients By Martin Skote, Dan S. Henningson and Ruud A. W. M. Henkes Direct numerical simulations of the Navier-Stokes
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 informationp + µ 2 u =0= σ (1) and
Master SdI mention M2FA, Fluid Mechanics program Hydrodynamics. P.-Y. Lagrée and S. Zaleski. Test December 4, 2009 All documentation is authorized except for reference [1]. Internet access is not allowed.
More informationFlow Transition in Plane Couette Flow
Flow Transition in Plane Couette Flow Hua-Shu Dou 1,, Boo Cheong Khoo, and Khoon Seng Yeo 1 Temasek Laboratories, National University of Singapore, Singapore 11960 Fluid Mechanics Division, Department
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 informationThe Phase Relationship in Laminar Channel Flow Controlled by Traveling Wave-Like Blowing/Suction. Abstract
The Phase Relationship in Laminar Channel Flow Controlled by Traveling Wave-Like Blowing/Suction Hiroya Mamori, Koji Fukagata, and Jerôme Hœpffner 2 Department of Mechanical Engineering, Keio University,
More informationMatrix-free methods for the stability and control of boundary layers
Matrix-free methods for the stability and control of boundary layers Shervin Bagheri, Espen Åkervik, Luca Brandt and Dan S. Henningson Linné Flow Centre, Department of Mechanics, Royal Institute of Technology,
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 informationApplied Mathematics and Mechanics (English Edition)
Appl. Math. Mech. -Engl. Ed., 39(9), 1267 1276 (2018) Applied Mathematics and Mechanics (English Edition) https://doi.org/10.1007/s10483-018-2364-7 Direct numerical simulation of turbulent flows through
More informationConvective instability and transient growth in flow over a backward-facing step
J. Fluid Mech. (28), vol. 63, pp. 271 34. c 28 Cambridge University Press doi:1.117/s221128119 Printed in the United Kingdom 271 Convective instability and transient growth in flow over a backward-facing
More informationEmmanuel PLAUT. Example of an hypertext Report written with L A TEX. This is a public document.
Emmanuel PLAUT Example of an hypertext Report written with L A TEX This is a public document. Contents Introduction 3 1 Presentation of the corporation 4 1.1 UL............................................
More informationChapter 3 Lecture 8. Drag polar 3. Topics. Chapter-3
Chapter 3 ecture 8 Drag polar 3 Topics 3.2.7 Boundary layer separation, adverse pressure gradient and favourable pressure gradient 3.2.8 Boundary layer transition 3.2.9 Turbulent boundary layer over a
More informationAnisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale
More informationOptimizing calculation costs of tubulent flows with RANS/LES methods
Optimizing calculation costs of tubulent flows with RANS/LES methods Investigation in separated flows C. Friess, R. Manceau Dpt. Fluid Flow, Heat Transfer, Combustion Institute PPrime, CNRS University
More informationNONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS
June - July, 5 Melbourne, Australia 9 7B- NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS Werner M.J. Lazeroms () Linné FLOW Centre, Department of Mechanics SE-44
More informationA METHODOLOGY FOR THE AUTOMATED OPTIMAL CONTROL OF FLOWS INCLUDING TRANSITIONAL FLOWS. Ronald D. Joslin Max D. Gunzburger Roy A.
A METHODOLOGY FOR THE AUTOMATED OPTIMAL CONTROL OF FLOWS INCLUDING TRANSITIONAL FLOWS Ronald D. Joslin Max D. Gunzburger Roy A. Nicolaides NASA Langley Research Center Iowa State University Carnegie Mellon
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
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 informationFLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes
SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.) INL 1 (3 cr.) 3 sets of home work problems (for 10 p. on written exam) 1 laboration TEN1 (4.5 cr.) 1 written exam
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 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 Fluid Mechanics
Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible
More informationFig. 1. The coordinate system and the directions of the velocities.
ICMAR 201 FLOW STABILITY NEAR A PLATE WITH A SURFACE MOVING AGAINST AN INCOMING STREAM A.M. Gaifullin, N.N. Kiselev, A.A. Kornyakov The Central Aerohydrodynamic Institute 10180, Zhukovsky, Moscow Reg.,
More informationLES of synthetic jets in boundary layer with laminar separation caused by adverse pressure gradient
LES of synthetic jets in boundary layer with laminar separation caused by adverse pressure gradient Tetsuya Ozawa, Samuel Lesbros and Guang Hong* University of Technology, Sydney (UTS), Australia *Corresponding
More informationDirect Numerical Simulations of Transitional Flow in Turbomachinery
Direct Numerical Simulations of Transitional Flow in Turbomachinery J.G. Wissink and W. Rodi Institute for Hydromechanics University of Karlsruhe Unsteady transitional flow over turbine blades Periodic
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 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 informationThe drag-adjoint field of a circular cylinder wake at Reynolds numbers 20, 100 and 500
The drag-adjoint field of a circular cylinder wake at Reynolds numbers 2, 1 and 5 The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationTHE CONVECTION DIFFUSION EQUATION
3 THE CONVECTION DIFFUSION EQUATION We next consider the convection diffusion equation ɛ 2 u + w u = f, (3.) where ɛ>. This equation arises in numerous models of flows and other physical phenomena. The
More informationFirst published on: 20 November 2010
This article was downloaded by: [Laval, Jean-Philippe] On: 20 November 2010 Access details: Access Details: [subscription number 929944761] Publisher Taylor & Francis Informa Ltd Registered in England
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 informationThe Reynolds experiment
Chapter 13 The Reynolds experiment 13.1 Laminar and turbulent flows Let us consider a horizontal pipe of circular section of infinite extension subject to a constant pressure gradient (see section [10.4]).
More informationIncreasing T T=1.5 T=5 T=10 T= time (viscous units) T=25 T=50. time (viscous units)
Proceedings of the 8th IFIP TC7 Conference on System Modelling and Optimization, Detroit, Michigan, July 997. CONTROL OF TURBULENT FLOWS R. Temam ;2, T. Bewley 3 and P. Moin 3 July 8, 997 Abstract. It
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 informationOPTIMUM SUCTION DISTRIBUTION FOR TRANSITION CONTROL *
OPTIMUM SUCTION DISTRIBUTION FOR TRANSITION CONTROL * P. Balakumar* Department of Aerospace Engineering Old Dominion University Norfolk, Virginia 23529 P. Hall Department of Mathematics Imperial College
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationLaminar Flow. Chapter ZERO PRESSURE GRADIENT
Chapter 2 Laminar Flow 2.1 ZERO PRESSRE GRADIENT Problem 2.1.1 Consider a uniform flow of velocity over a flat plate of length L of a fluid of kinematic viscosity ν. Assume that the fluid is incompressible
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationFlow control applied to transitional flows: input-output analysis, model reduction and control
Flow control applied to transitional flows: input-output analysis, model reduction and control Dan Henningson collaborators Shervin Bagheri, Espen Åkervik Luca Brandt, Peter Schmid 1 Outline Introduction
More informationDYNAMICS OF CONTROLLED BOUNDARY LAYER SEPARATION
p.1 DYNAMICS OF CONTROLLED BOUNDARY LAYER SEPARATION Václav Uruba, Martin Knob Institute of Thermomechanics, AS CR, v. v. i., Praha Abstract: The results of experimental study on a boundary layer separation
More informationNonlinear shape evolution of immiscible two-phase interface
Nonlinear shape evolution of immiscible two-phase interface Francesco Capuano 1,2,*, Gennaro Coppola 1, Luigi de Luca 1 1 Dipartimento di Ingegneria Industriale (DII), Università di Napoli Federico II,
More informationLINEAR STABILITY ANALYSIS FOR THE HARTMANN FLOW WITH INTERFACIAL SLIP
MAGNETOHYDRODYNAMICS Vol. 48 (2012), No. 1, pp. 147 155 LINEAR STABILITY ANALYSIS FOR THE HARTMANN FLOW WITH INTERFACIAL SLIP N. Vetcha, S. Smolentsev, M. Abdou Fusion Science and Technology Center, UCLA,
More informationSIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS
ICAS 2000 CONGRESS SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS H Yao, R K Cooper, and S Raghunathan School of Aeronautical Engineering The Queen s University of Belfast, Belfast BT7 1NN,
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationFLUID MECHANICS. ! Atmosphere, Ocean. ! Aerodynamics. ! Energy conversion. ! Transport of heat/other. ! Numerous industrial processes
SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.)! INL 1 (3 cr.)! 3 sets of home work problems (for 10 p. on written exam)! 1 laboration! TEN1 (4.5 cr.)! 1 written
More informationBifurcation and stability analysis of a jet in crossflow. Part 1: Onset of global instability at a low velocity ratio
Page 1 of 27 Under consideration for publication in J. Fluid Mech. 1 Bifurcation and stability analysis of a jet in crossflow. Part 1: Onset of global instability at a low velocity ratio MILOŠ ILAK, PHILIPP
More informationApplication of Reconstruction of Variational Iteration Method on the Laminar Flow in a Porous Cylinder with Regressing Walls
Mechanics and Mechanical Engineering Vol. 21, No. 2 (2017) 379 387 c Lodz University of Technology Application of Reconstruction of Variational Iteration Method on the Laminar Flow in a Porous Cylinder
More informationDRAG REDUCTION IN SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYERS BY BLOWING AT CONSTANT MASS-FLUX
June - July, Melbourne, Australia 9 A- DRAG REDUCTION IN SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYERS BY BLOWING AT CONSTANT MASS-FLUX Yukinori Kametani Linné FLOW Centre KTH Mechanics Osquars Backe
More informationLaminar and turbulent comparisons for channel flow and flow control
J. Fluid Mech. (2007), vol. 570, pp. 467 477. c 2007 Cambridge University Press doi:10.1017/s0022112006003247 Printed in the United Kingdom 467 Laminar and turbulent comparisons for channel flow and flow
More informationZonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach
Zonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach S. Arvidson 1,2, L. Davidson 1, S.-H. Peng 1,3 1 Chalmers University of Technology 2 SAAB AB, Aeronautics 3 FOI, Swedish Defence
More informationEXPERIMENTS OF CLOSED-LOOP FLOW CONTROL FOR LAMINAR BOUNDARY LAYERS
Fourth International Symposium on Physics of Fluids (ISPF4) International Journal of Modern Physics: Conference Series Vol. 19 (212) 242 249 World Scientific Publishing Company DOI: 1.1142/S211945128811
More informationLarge-eddy simulations for wind turbine blade: rotational augmentation and dynamic stall
Large-eddy simulations for wind turbine blade: rotational augmentation and dynamic stall Y. Kim, I.P. Castro, and Z.T. Xie Introduction Wind turbines operate in the atmospheric boundary layer and their
More informationNumerical Simulation of Flow Around An Elliptical Cylinder at High Reynolds Numbers
International Journal of Fluids Engineering. ISSN 0974-3138 Volume 5, Number 1 (2013), pp. 29-37 International Research Publication House http://www.irphouse.com Numerical Simulation of Flow Around An
More information