Stability analysis and control of the flow in a symmetric channel with a sudden expansion

The laminar flow in two-dimensional diffusers may produce either symmetric or nonsymmetric steady solutions, depending on the value of the Reynolds number as compared with some critical value. The stability properties of the flow are studied in the context of linear theory. In this context, a sensitivity analysis of the flow instability is carried out with respect to perturbations that may be produced by a realistic passive control, thus providing qualitative hints and quantitative information for the control design. Following the so-obtained information, a passive control is built by introducing a small cylinder in the flow with the aim of stabilizing the unstable symmetric flow configuration in the diffuser. The effectiveness of this control is finally assessed by direct numerical simulation. It is shown that the introduction of the cylinder, placed following the indications of the linear sensitivity analysis in the stable asymmetric flow configuration, allows a steady completely symmetric or less asymmetric flow to be recovered. The flow transient between the uncontrolled asymmetric solution and the symmetric controlled one is analyzed in terms of streamlines and vorticity evolution; the effects of the cylinder introduction on flow dissipation are also assessed. © 2012 American Institute of Physics.

Unsteady fluid dynamics around a hovering flat plate wing

Swathi Badala Krishna

Insect flight is characterised by complex time-dependent flows in response to the unsteady wing movements. Biological fliers exploit the unsteady flow fields to modulate aerodynamic forces, thereby displaying unmatched flight performance, especially in hover. Naturally, this has inspired the creation of engineering models to replicate the flight behaviour. An in-depth understanding of the flow fields generated during hover and their dependence on the kinematics is paramount to achieve this goal. The two main kinematic components of a hovering wing are the stroke, which refers to the back-and-forth motion, and wing rotation, which refers to the change in angle of attack. The phase relation between stroke and rotation is quantified in terms of phase-shift and is broadly classified into symmetric, advanced, and delayed rotation. The phase-shift and duration of rotation, together referred to as rotational timing, are investigated in this bio-inspired study. The objective is to characterise the effect of rotational timing on the aerodynamic forces and the flow fields generated by a hovering wing. The unsteady flow around a hovering flat plate wing that mimics hoverfly kinematics has been investigated experimentally using particle image velocimetry and direct force measurements. The measurements are conducted at a Reynolds number of Re=220 and a reduced frequency of k=0.32 in order to dynamically match a hoverfly. The Lagrangian finite-time Lyapunov exponent method is used to analyse the unsteady flow fields by identifying dynamically relevant flow features such as the primary leading edge vortex (LEV), secondary vortices, and topological saddles, and their evolution within a flapping cycle. Firstly, the flow and force behaviour was characterised for a typical flapping cycle. The flow evolution in a symmetric, fast rotation is divided into four stages that are characterised by the LEV emergence, growth, lift-off, and breakdown and decay. Tracking saddle points is shown to be helpful in defining the LEV lift-off which occurs at the maximum stroke velocity. The flow fields are correlated with the aerodynamic forces revealing that the maximum lift and drag are observed just before LEV lift-off, which corresponds to the maximum stroke velocity. Secondly, the effect of phase-shift on the formation and evolution of lift-enhancing flow structures are discussed. Two advanced and delayed rotations are compared. The flow development stages and forces are similar for all rotations but the timing of stages varies. The evolution of forces and flow strongly depend on the stroke velocity. Thirdly, the dependence of the flow and force evolution on the stroke velocity was substantiated by doubling the rotational duration in the symmetric rotation. It was found that the timing of the flow stages altered, whereas the flow and forces mostly evolved similarly to that of a fast rotation. The fast rotation, however, produces higher maximum lift and drag compared to the slow rotation. Lastly, the effect of phase-shift on the aerodynamic characteristics of a slow rotation is further explored. The slow rotation cases exhibit distinct flow patterns for varying phase-shifts unlike the fast rotations, in terms of the formation, evolution and breakdown of the flow structures as well as the timing. The forces also show distinct trends for varying phase-shifts and strongly depend on the angle of attack along with the stroke velocity in the slow

EPFL2017

Simulation of benchmark and industrial unsteady compressible turbulent fluid flows

Michele Mossi

Large-eddy simulation (LES) is a very promising technique for the numerical computation of unsteady turbulent flows, and seems to be the perfect tool to simulate the compressible air flow around a high-speed train in a tunnel, providing unsteady results for aerodynamic and aeroacoustic analysis. To look into this possible future application of LES, two major lines of investigation are pursued in this thesis: first, the study of the effective ability of shock-capturing schemes to predict fundamental turbulent phenomena; second, the analysis of the aerodynamic phenomena induced by a high-speed train in a tunnel. The numerical simulation of compressible flows requires the use of shock-capturing schemes. These schemes can be relatively dissipative and mask the subgrid-scale contribution introduced in a large-eddy simulation to account for the unresolved turbulence scales. To estimate their effective dissipation and their ability to resolve turbulence phenomena, shock-capturing schemes widely used for aeronautical applications, from second- to fifth-order space accuracy, are employed here for simulating well-known fundamental flows in subsonic and supersonic regimes. Direct and large-eddy numerical simulations are undertaken for the inviscid and viscous Taylor-Green vortex decay problem, the freely decaying homogeneous and isotropic turbulence, and the fully developed channel flow. For all the turbulent flows investigated, several turbulence statistics are computed and the numerical dissipation of the schemes tested is interpreted in terms of subgrid-scale dissipation in a LES context, yielding an equivalent Smagorinsky or dynamic coefficient. This coefficient is for all schemes of the same order of magnitude as the commonly accepted values in LES for the subgrid-scale term. On the grounds of this analysis and of the comparisons of the turbulence statistics with accurate data obtained in the literature, the addition of explicit subgrid-scale models to the shock-capturing schemes tested is not recommended. It is thus concluded that the use of the LES technique for compressible turbulent flows is not yet suitable for industrial applications. The aerodynamic phenomena generated by a high-speed train travelling in a tunnel are also discussed, their importance on the design of high-speed lies is pointed out, and the analysis tools commonly employed for their study are reviewed. To study numerically the three-dimensional, compressible and turbulent air flow around a high-speed train accelerating in a tunnel, by accounting for the unsteady effects at inlet and outlet boundaries due to the propagation of pressure waves generated at the train departure, new coupling conditions between one-dimensional and three-dimensional domains are developed. These conditions are applied successfully to the numerical simulation of the unsteady wake developing behind two- and three-dimensional vehicles, where the averaged Navier-Stokes equations are solved with the turbulence modelling approach. The influence on the wake of the length of the vehicle tail is also discussed and results of multi-dimensional simulations are compared with one-dimensional data.

EPFL1999

Simulation of shear-driven flows

Roland Bouffanais

The research work reported in the present dissertation is aimed at the analysis of complex physical phenomena involving instabilities and nonlinearities occurring in fluids through state-of-the-art numerical modeling. Solutions of intricate fluid physics problems are devised in two particularly arduous situations: fluid domains with moving boundaries and the high-Reynolds-number regime dominated by nonlinear convective effects. Shear-driven flows of incompressible Newtonian fluids enclosed in cavities of varying geometries are thoroughly investigated in the two following frameworks: transition with a free surface and confined turbulence. The physical system we consider is made of an incompressible Newtonian fluid filling a bounded, or partially bounded cavity. A series of shear-driven flows are easily generated by setting in motion some part of the container boundary. These driven-cavity flows are not only technologically important, they are of great scientific interest because they display almost all physical fluid phenomena that can possibly occur in incompressible flows, and this in the simplest geometrical settings. Thus corner eddies, secondary flows, longitudinal vortices, complex three-dimensional patterns, chaotic particle motions, nonuniqueness, transition, and turbulence all occur naturally and can be studied in the same geometry. This facilitates the comparison of results from experiments, analysis, and computation over the whole range of Reynolds numbers. The flows under consideration are part of a larger class of confined flows driven by linear or angular momentum gradients. This dissertation reports a detailed study of a novel numerical method developed for the simulation of an unsteady free-surface flow in three-space-dimensions. This method relies on a moving-grid technique to solve the Navier-Stokes equations expressed in the arbitrary Lagrangian-Eulerian (ALE) kinematics and discretized by the spectral element method. A comprehensive analysis of the continuous and discretized formulations of the general problem in the ALE frame, with nonlinear, non-homogeneous and unsteady boundary conditions is presented. In this dissertation, we first consider in the internal turbulent flow of a fluid enclosed in a bounded cubical cavity driven by the constant translation of its lid. The solution of this flow relied on large-eddy simulations, which served to improve our physical understanding of this complex flow dynamics. Subsequently, a novel subgrid model based on approximate deconvolution methods coupled with a dynamic mixed scale model was devised. The large-eddy simulation of the lid-driven cubical cavity flow based on this novel subgrid model has shown improvements over traditional subgrid-viscosity type of models. Finally a new interpretation of approximate deconvolution models when used with implicit filtering as a way to approximate the projective grid filter was given. This led to the introduction of the grid filter models. Through the use of a newly-developed method of numerical simulation, in this dissertation we solve unsteady flows with a flat and moving free-surface in the transitional regime. These flows are the incompressible flow of a viscous fluid enclosed in a cylindrical container with an open top surface and driven by the steady rotation of the bottom wall. New flow states are investigated based on the fully three-dimensional solution of the Navier-Stokes equations for these free-surface cylindrical swirling flows, without resorting to any symmetry properties unlike all other results available in the literature. To our knowledge, this study delivers the most general available results for this free-surface problem due to its original mathematical treatment. This second part of the dissertation is a basic research task directed at increasing our understanding of the influence of the presence of a free surface on the intricate transitional flow dynamics of shear-driven flows.

EPFL2007