Flow control of transitional shear flows is investigated by means of numerical simulations. The attenuation of three-dimensional wavepackets of Tollmien-Schlichting (TS) and streaks in the boundary layer is obtained using active control in combination with localised sensors and actuators distributed near the rigid wall. Due to the dimensions of the discretized Navier-Stokes operator, reduced-order models are identified, preserving the dynamics between the inputs and the outputs of the system. Balanced realizations of the system are computed using balanced truncation and system identification. We demonstrate that the energy growth of the perturbations is substantially and efficiently mitigated, using relatively few sensors and actuators. The robustness of the controller is analysed by varying the number of actuators and sensors, the Reynolds number, the pressure gradient and by investigating the nonlinear, transitional case. We show that delay of the transition from laminar to turbulent flow can be achieved despite the fully linear approach. This configuration can be reproduced in experiments, due to the localisation of sensing and actuation devices. The closed-loop system has been investigated for the corresponding twodimensional case by using full-dimensional optimal controllers computed by solving an iterative optimisation based on the Lagrangian approach. This strategy allows to compare the results achieved using open-loop model reduction with model-free controllers. Finally, a parametric analysis of the actuators/ sensors placement is carried-out to deepen the understanding of the inherent dynamics of the closed-loop. The distinction among two different classes of controllers – feedforward and feedback controllers - is highlighted. A second shear flow, a confined turbulent jet, is investigated using particle image velocimetry (PIV) measurements. Proper orthogonal decomposition (POD) modes and Koopman modes via dynamic mode decomposition (DMD) are computed and analysed for understanding the main features of the flow. The frequencies related to the dominating mechanisms are identified; the most energetic structures show temporal periodicity. / <p>QC 20130207</p>

The spectral-element method (SEM) is used to study wall-bounded turbulent flowsin moderately complex geometries. The first part of the thesis is devoted to simulations of canonical flow cases, such as temporal K-type transitionand turbulent channel flow, to investigate general resolution requirements and computational efficiency of the numerical code nek5000. Large-eddy simulation (LES) is further performed of a plane asymmetric diffuser flow with an opening angle of 8.5 degrees, featuring turbulent flow separation. Good agreement with numerical studies of Herbst (2007) is obtained, and it is concluded that the use of a high-order method is advantageous for flows featuring pressure-induced separation. Moreover, it is shown, both a priori on simpler model problems and a posteriori using the full Navier--Stokes equations, that the numerical instability associated with SEM at high Reynolds numbers is cured either by employing over-integration (dealiasing) or a filter-based stabilisation, thus rendering simulations of moderate to high Reynolds number flows possible. The second part of the thesis is devoted to the first direct numerical simulation (DNS) of a truly three-dimensional, turbulent and separated diffuser flow at Re = 10 000 (based on bulk velocity and inflow-duct height), experimentally investigated by Cherry et al. (2008). The massively parallel capabilities of the spectral-element method are exploited by running the simulations on up to 32 768 processors. Very good agreement with experimental mean flow data is obtained and it is thus shown that well-resolved simulations of complex turbulent flows with high accuracy are possible at realistic Reynolds numberseven in complicated geometries. An explanation for the discovered asymmetry of the mean separated flow is provided and itis demonstrated that a large-scale quasi-periodic motion is present in the diffuser. In addition, a new diagnostic measure, based on the maximum vorticity stretching component in every spatial point, is designed and tested in a number of turbulent and transitional flows. Finally, Koopman mode decomposition is performed of a minimal channel flow and compared to classical proper orthogonal decomposition (POD). / QC 20111206

In this thesis direct numerical simulations are used to investigate two phenomenain shear flows: laminar-turbulent transition over a flat plate and periodicvortex shedding induced by a jet in cross flow. The emphasis is on understanding and controlling the flow dynamics using tools from dynamical systems and control theory. In particular, the global behavior of complex flows is describedand low-dimensional models suitable for control design are developed; this isdone by decomposing the flow into global modes determined from spectral analysisof various linear operators associated with the Navier–Stokes equations.Two distinct self-sustained global oscillations, associated with the sheddingof vortices, are identified from direct numerical simulations of the jet incrossflow. The investigation is split into a linear stability analysis of the steadyflow and a nonlinear analysis of the unsteady flow. The eigenmodes of theNavier–Stokes equations, linearized about an unstable steady solution revealthe presence of elliptic, Kelvin-Helmholtz and von K´arm´an type instabilities.The unsteady nonlinear dynamics is decomposed into a sequence of Koopmanmodes, determined from the spectral analysis of the Koopman operator. Thesemodes represent spatial structures with periodic behavior in time. A shearlayermode and a wall mode are identified, corresponding to high-frequency andlow-frequency self-sustained oscillations in the jet in crossflow, respectively.The knowledge of global modes is also useful for transition control, wherethe objective is to reduce the growth of small-amplitude disturbances to delaythe transition to turbulence. Using a particular basis of global modes, knownas balanced modes, low-dimensional models that capture the behavior betweenactuator and sensor signals in a flat-plate boundary layer are constructed andused to design optimal feedback controllers. It is shown that by using controltheory in combination with sensing/actuation in small, localized, regionsnear the rigid wall, the energy of disturbances may be reduced by an order of magnitude.