Transition to turbulence and flow control are studied by means of numerical simulations for different simple shear flows. Linear and non-linear optimisation methods using the Lagrange multiplier technique are employed. In the linear framework as objective function the standard disturbance kinetic energy is chosen and the constraints involve the linearised Navier–Stokes equations. We consider both the optimal initial condition leading to the largest disturbance energy growth at finite times and the optimal time-periodic forcing leading to the largest asymptotic response for the case of the flat plate boundary layer excluding the leading edge. The optimal disturbances for spanwise wavelengths of the order of the boundary layer thickness are streamwise vortices exploiting the lift-up mechanism to create streaks. For long spanwise wavelengths it is the Orr mechanism combined with the amplification of oblique wave packets that is responsible for the disturbance growth. Also linear optimal disturbances are computed around a leading edge and the effect of the geometry is considered. It is found that two-dimentional disturbances originating upstream, relative to the leading edge of the plate are inefficient at generating a viable disturbance, while three dimentional disturbances are more amplified. In the non-linear framework a new approach using ideas from non-equilibrium thermodynamics is developed. We determine the initial condition on the laminar/turbulent boundary closest to the laminar state. Starting from the general evolution criterion of non-equilibrium systems we propose a method to optimise the route to the statistically steady turbulent state, i.e. the state characterised by the largest entropy production. This is the first time information from the fully turbulent state is included in the optimisation procedure. The method is applied to plane Couette flow. We show that the optimal initial condition is localised in space for realistic flow domains, while the disturbance visits bent streaks before breakdown. Feedback control is applied to the bypass-transition scenario with high levels of free-stream turbulence. The flow is the flat-plate boundary layer. In this scenario low frequency perturbations enter the boundary layer and streamwise elongated disturbances emerge due to non-modal growth. The so-called streaky structures are growing in amplitude until they reach high enough energy levels and break down into turbulent spots via their secondary instability. When control is applied in the form of wall blowing and suction, the growth of the streaks is delayed, which implies a delay of the whole transition process. Additionally, a comparison with experimental work is performed demonstrating a remarkable agreement in the disturbance attenuation once the differences between the numerical and experimental setup are reduced. Open-loop control with wall travelling waves by means of blowing and suction is applied to a separating boundary layer. For downstream travelling waves we obtain a mitigation of the separation of the boundary layer while for upstream travelling waves a significant delay in the transition location accompanied by a modest reduction of the separated region. / QC 20110518

A Direct Numerical Simulation of an iced Natural Laminar Flow NLF-0414 airfoil is carried out using a high-order spectral element method for low chord Reynolds numbers (O(10^5)). This study aims to advance the state-of-the-art for accurate computational modeling of transition, iced airfoil aerodynamics, and irregular surface spectral element method Direct Numerical Simulation. Ice accretion over an aircraft, ranging from light to severe, changes the aerodynamic profile of the airfoil and alters the overall performance. The literature presents simulations that have been carried out with a range of turbulence models which fail to accurately capture the complex physics of these flows. The iced profiles being studied, Run 606 and 622-2D, were obtained from a Technical Publication by NASA on iced airfoils including the NLF-0414, and were selected as they are relatively lightly iced profiles of the NLF-0414. The largest bottleneck with the current advancement in High Performance Computing is the computation time required for Direct Numerical Simulation. Results such as lift, drag, pressure, and skin friction coefficients, for a clean NLF-0414 and two lightly iced NLF-0414 airfoils at chord Reynolds numbers of Rec = 1 x 10^5 and Rec = 2 x 10^5 are visualized and discussed, showing the degradation of the natural laminar flow due to ice accretion. Turbulence statistics are calculated to study the effective contributions of turbulent fluctuations in the flow to further understand the flow physics near transition. The detailed study of these six cases has led us to 1) further understand the complexities of the transition process on iced airfoils, 2) observe and explain the sometimes unexpected changes in aerodynamic performance due to varying iced geometries, and 3) establish a methodology for spectral element method Direct Numerical Simulations.

An experimental investigation was performed to ascertain the nature of turbulence in a narrow apex angle isosceles triangular duct. The study involved the design and construction of a low noise, low turbulence wind tunnel that had an isosceles triangular test section with an apex angle of 11.5[degrees]. Experiments involved the measurement of velocity fluctuations using hot wire anemometry and wall pressure fluctuations using a condenser microphone. Measurement of the velocity fluctuations reconfirms the coexistence of laminar and turbulent regions at a given cross section for a range of Reynolds numbers. The laminar region is concentrated closer to the apex while the turbulent region is found closer to the base. The point of transition is a function of the Reynolds number and moves closer to the apex as the flow rate is increased. Moreover, it was found in this investigation that traditional scaling of the turbulent statistical quantities do not hold good in this geometry. Although velocity fluctuations showed distinctive flow regimes, no such distinction could be seen in the dynamic wall pressure data. The nature of the dynamic wall pressure was uniform throughout the entire cross section suggesting that wall pressure fluctuations, unlike the velocity fluctuations, are able to travel from the base to the apex, without being damped. This implies that the relationship between the velocity and the pressure fluctuations applicable in the other systems does not hold well in a narrow apex angle isosceles triangular duct. Further, the typical scaling relationships applied to wall pressure spectra of other geometries doesn't apply in this scenario and the ratio of the RMS pressure fluctuation to the mean shear is much higher compared to a flat plate or pipe flow situation.

The focus of this thesis is the numerical study of subcritical transition to turbulence in boundary-layer flows. For the most part, boundary layers with uniform suction are considered. Constant homogeneous suction counteracts the spatial growth of the boundary layer, rendering the flow parallel. This enables research approaches which are not feasible in the context of spatially developing flows. In the first part, the laminar–turbulent separatrix of the asymptotic suction boundary layer (ASBL) is investigated numerically by means of an edge-tracking algorithm. The obtained edge states experience recurrent dynamics, going through calm and bursting phases. The self-sustaining mechanism bears many similarities with the classical regeneration cycle of near-wall turbulence. The recurrent simple structure active during calm phases is compared to the nucleation of turbulence events in bypass transition originating from delocalised initial conditions. The implications on the understanding of the bypass-transition process and the edge state's role are discussed. Based on this understanding, a model is constructed which predicts the position of the nucleation of turbulent spots during free-stream turbulence induced transition in spatially developing boundary-layer flow. This model is used together with a probabilistic cellular automaton (PCA), which captures the spatial spreading of the spots, correctly reproducing the main statistical characteristics of the transition process. The last part of the thesis is concerned with the spatio-temporal aspects of turbulent ASBL in extended numerical domains near the onset of sustained turbulence. The different behaviour observed in ASBL, i.e. absence of sustained laminar–turbulent patterns, which have been reported in other wall-bounded flows, is associated with different character of the large-scale flow. In addition, an accurate quantitative estimate for the lowest Reynolds number with sustained turbulence is obtained / <p>QC 20160429</p>

L'objectif de cette thèse est de fournir une analyse de la transition vers la turbulence d'un fluide rhéofluidifiant (fluide de Carreau) dans une conduite cylindrique. Pour cela, un code pseudo-spectral de type Petrov-Galerkin a été développé. Une analyse linéaire de stabilité de l'écoulement laminaire est effectuée, montrant que cet écoulement est linéairement stable. Ensuite, des perturbations sous la forme des rouleaux longitudinaux contra-rotatifs sont utilisées comme condition initiale. Les termes non linéaires d'inertie et visqueux créent un écoulement secondaire avec des points d'inflexion, linéairement instable vis-à-vis de perturbations 3D. Une analyse linéaire de stabilité de ce nouvel écoulement de base bidimensionnelle est réalisée. La forme des vecteurs propres critiques est analysé. Enfin, une analyse non linéaire de stabilité de rouleaux vis-à-vis des perturbations tridimensionnelles de faible amplitude est effectuée, obtenant un retard pour la transition vers la turbulence des fluides rhéofluidifiants par rapport au cas Newtonien et une tendance à l'asymétrie du profil de vitesse axiale / The main objective of this thesis is to provide a description of the transition to turbulence of a shear thinning fluid in pipe flow. A linear stability analysis of the base flow is done. Results show that the flow is linearly stable and the optimal perturbation is given by a pair of counter rotating vortex. This kind of perturbation is used as an initial condition of a computational code which integrates the governing equations. Inertial and viscous non linear terms generate a secondary base flow with inflection points, which is linearly unstable to 3D perturbations. A secondary instability analysis is done, regarding the shape of unstable eigenvectors. Depending the rheological parameters and the size of the primary perturbation, the unstable mode might be near the wall or the center of the pipe. Finally, a non linear stability analysis of the streaks to 3D perturbations of weak amplitude, obtaining a delay in the transition to turbulence due to shear thinning

La thèse traite de recherches numériques et expérimentales sur l’écoulement à traves des conduites circulaires ou des tubes avec une petite entrée et un diamètre de sortie plus grand, parfois appelées élargissement ou divergents. L’écoulement dans un élargissement est globalement stable pour des nombres de Reynolds élevés. Ainsi la simulation numérique de ce type d’écoulement nécessite de grands domaines de calcul contenant la zone de recirculation, qui croît linéairement. En outre, les études expérimentales dans les élargissements brusques indiquent que la transition se produit à des nombres de Reynolds plus faibles que prévue par la théorie linéaire de stabilité. La raison pour cette transition précoce est due à la présence d’imperfections dans le dispositif expérimental, qui agit comme une perturbation d’amplitude finie de l’écoulement. Des simulations numériques directes des équations de Navier-Stokes ont été réalisées avec deux types différents de perturbations (i) l’inclination et (ii) le vortex. Tout d’abord, la perturbation de type inclinaison, qui est appliqué à l’entrée, crée une zone de recirculation asymétrique, puis se casse pour former une turbulence localisée en aval de l’expansion. Deuxièmement, la perturbation de type vortex, crée des structures qui ressemblent à un mode azimutal d’ordre inférieur, déjà identifié comme une perturbation optimale amplifiée. Il croît en raison de l’instabilité convective, puis forme une tâche de turbulence localisée. Enfin, la corrélation spatiale et la décomposition en modes propres révèlent que cette turbulence localisée obtient son énergie de l’écoulement d’entrée. / The thesis deals with numerical and experimental investigations of flow through circular pipes with smaller inlet and larger outlet diameter, also known as expansion pipes. The hydrodynamic expansion pipe flow is globally stable for high Reynolds number. In order to numerically simulate these types of flows, large computational domains that could accommodate the linearly growing symmetric recirculation region is needed. Moreover, experimental studies of expansion pipe flows indicate that the transition occurs at lower Reynolds number than predicted by the linear stability theory. The reason for early transition is due to the presence of imperfections in the experimental setup, which acts as a finite-amplitude perturbation of the flow. Three-dimensional direct numerical simulations of the Navier-Stokes equations with two different types of perturbations (i) the tilt and (ii) the vortex are investigated. First, the tilt perturbation, which applied at the inlet, creates an asymmetric recirculation region and then breaks to form localised turbulence downstream the expansion section. Second, the vortex perturbation, creates structures that looks like lower order azimuthal mode, resembles an optimally amplified perturbation. It grows due to convective instability mechanism and then breaks to form localised turbulence. Spatial correlation and the proper orthogonal decomposition reveal that this localised turbulence gains it energy from the core flow coming out of the inlet pipe.

No description available.

Cette thèse vise à démêler les principaux mécanismes impliqués dans les écoulements transitoires et turbulents. L’idée centrale est d'utiliser une technique d’optimisation non linéaire pour étudier l’origine et le rôle des structures cohérentes habituellement observées dans ces écoulements. Cette méthode a été utilisée dans trois contextes différents. Tout d’abord, un écoulement laminaire linéairement stable a été considéré et l'optimisation a été utilisée pour calculer les perturbations les plus amplifiées parmi toutes les perturbations capables de déclencher une transition vers la turbulence. Une fois que la turbulence est bien établie, une optimisation non linéaire entièrement 3D maximisant l'énergie cinétique turbulente est utilisée pour étudier les structures cohérentes qui peuplent l’écoulement turbulent et les mécanismes responsables de la croissance et de l’échange d’énergie (optimale) sont étudiés. Ensuite, une approche de type système dynamique est appliquée aux équations du mouvement. La géométrie de l’espace des phases est étudiée en utilisant la théorie de la croissance transitoire pour évaluer l’importance des variétés stable et instable dans la dynamique. Dans le même cadre, un algorithme de minimisation non linéaire est utilisé pour calculer les connexions hétérocliniques parmi les solutions invariantes des équations de Navier-Stokes. / This thesis aims at unraveling the main mechanisms involved in transitional and turbulent flows. The central idea is that of using a nonlinear optimization technique to investigate the origin and role of coherent structures usually observed in these flows. This method has been used in three different contexts. First, a linearly stable laminar flow has been considered and the optimization has been used to compute the most amplified perturbations among all disturbances able to trigger transition to turbulence. Once turbulence is well established, a fully 3D nonlinear optimization maximizing the turbulent kinetic energy is used to study coherent structures populating turbulent shear flow as well as investigate the mechanisms responsible for the energy (optimally) growth and exchange. Then, a dynamical system approach is applied to fluid flow equations. The geometry of the state space is investigated by using transient growth theory to reveal the importance of the stable and unstable manifold. In the same framework, a nonlinear minimization algorithm is used to compute heteroclinic connections among invariant solutions of the Navier-Stokes equations.

Cette thèse vise à mettre en évidence les limites du contrôle passif en utilisant des éléments de rugosité miniaturisés. La topologie des écoulements induite par la présence d’une rugosité cylindrique et des générateurs de tourbillons miniaturisés a été étudiée pour analyser la dynamique pour des temps court et long. Différentes bifurcations supercritiques ont été examinées au moyen d’une analyse de stabilité globale. La bifurcation souscritique est déclenchée par des mécanismes de croissance transitoire de l'énergie ou par la réceptivité de modes globaux stables. Des structures deforçage optimal 3D sont extraites pour comprendre la distribution spatiale liée à la fréquence de résonance du système. La simulation numérique directe perturbée révèle le rôle central du mode global le moins stable dans les instationnarités non-linéaires observées. Une analyse détaillée des structures tourbillonnaires montre qu’elles sont principalement liées aux mécanismes linéaires sous-jacents. La principale caractéristique globale du mode propre est liée à la présence d’une zone de séparation en aval de la rugosité cylindrique. En utilisant des générateurs de tourbillon miniaturisés, cette zone de séparation est fortement diminuée et aucun mode global isolé est alors présent. La dynamique de l’écoulement se révèle être conduit non seulement par le nombre de Reynolds de rugosité et par son rapport d'aspect géométrique, mais aussi par le rapport entre la hauteur de la rugosité et l’epaisseur de couche limite. / This thesis aims at highlighting the limits of passive control by usingminiaturized roughness elements. The flow topology induced by the presence ofcylindrical roughness and miniaturized vertex generators has been studied to uncover asymptotic and short time dynamics. Supercritical bifurcations has been investigated by means of global stability analysis. Subcritical bifurcation are induced by transient growth of the energy or receptivity of stable global modes. 3D optimal forcing structures are extracted to figure out the spatial distribution linked to the resonant pulsation. Perturbed direct numerical simulation reveals the pivotal role of the less steady global mode in the non-linear unsteadiness. A detailed analysis of the flowstructures is provided and linked to the involved linear mechanisms. Global feature of the eigenmode are linked to the presence of the separation zone behind the cylindrical roughness. By using miniaturized vortex generators the separation zone is suppressed and no isolated global modes are present. The flow dynamics turns out to be driven not only by roughness Reynolds number and geometrical aspect ratio but also by the ratio between the roughness height and the boundary layer thickens.

Quand on augmente le nombre de Reynolds d’un écoulement, il passe progressivement d’un état ordonné, laminaire, à un état chaotique, turbulent. Lors de cette transition, on voit apparaître des structures hydrodynamiques ordonnées qui se désordonnent sous l’effet de déstabilisations successives. Ces structures peuvent poser problème lorsqu’on souhaite par exemple homogénéiser des suspensions car des particules peuvent s’y retrouver piégées et générer des zones sur ou sous-concentrées. La plupart des fluides naturels et industriels sont non newtoniens. Pour mieux comprendre et contrôler l’apparition des structures hydrodynamiques dans de tels fluides, nous nous sommes intéressés dans le cadre de cette thèse, aux mécanismes d’instabilités gouvernant la transition vers la turbulence de fluides complexes dans la configuration simplifiée de Taylor-Couette. Dans un premier volet théorique et numérique, les concepts d’instabilités hydrodynamiques ont été appliqués à l’écoulement de Taylor-Couette pour trois modèles rhéologiques : la loi de puissance, le modèle de Carreau et le modèle de Bingham. Deux rapports de rayons et différentes vitesses de rotation des deux cylindres ont été considérés. La déstabilisation de l’écoulement de Couette circulaire a été étudiée pour prédire l’apparition de structures telles que les rouleaux de Taylor-Couette aux temps longs. Une analyse aux temps courts a également été effectuée pour étudier les phénomènes de croissance transitoire dans l’écoulement. Dans un second volet expérimental, un dispositif de Taylor-Couette en plexiglas a été élaboré. Dans un premier temps, des mesures de couple avec un fluide newtonien (glycérol) ont été réalisées et comparées avec succès avec les résultats de la littérature. Des mesures rhéométriques ont ensuite été réalisées pour choisir des fluides répondant aux lois rhéologiques utilisées dans le volet théorique. Des visualisations ont été effectuées pour détecter l’apparition des instabilités et les conditions critiques ont été comparées aux prédictions de l’analyse linéaire de stabilité. Des mesures de vélocimétrie par image de particules (PIV) ont été aussi effectuées pour déterminer les champs de vitesse de l’écoulement. Ces différentes études ont permis de mettre en évidence le caractère sous-critique et l’apparition d’un écoulement chaotique engendrés par le caractère rhéofluidifiant des fluides. / When increasing the Reynolds number of a flow, it passes gradually from an ordered state, called laminar, to a chaotic state, called turbulent. During this transition, organized hydrodynamic structures appear which disorganize under successive destabilizations. These structures may cause problems for example in the relation to the homogenization of suspensions, because particles may be trapped in the structures, leading to under or over-concentrated zones. Furthermore, most of the fluids found in nature and in the industry are non-Newtonian. In order to better understand and control the apparition of hydrodynamic structures in such fluids, we focused during this PhD on the instability mechanisms governing the transition to turbulence of complex fluids in the simplified configuration of Taylor-Couette. In a theoretical and numerical approach, concepts of hydrodynamic instabilities were applied to the circular Couette flow for three rheological models: the power law, the Carreau and the Bingham models. Two radius ratios and different rotation velocities were considered for the cylinders. The destabilization of the circular Couette flow was studied to predict the apparition of structures such as Taylor vortex at long times. An analysis at short time was also performed to study the phenomena of transient growth. In an experimental approach, a Taylor-Couette device was designed. First, the torque on the inner cylinder was measured for a Newtonian fluid (glycerol) and successfully compared to previous results of the literature. Secondly, rheological measurements were conducted in order to select fluids showing rheological behaviors considered in the theoretical approach. Visualizations were then performed in order to detect the apparition of instabilities and the critical conditions were compared to the predicted by means of the linear stability analysis. Finally, PIV measurements provided velocity fields in the flow. These different studies evidenced a subcritical transition and chaotic flow due to the shear-thinning character of the fluids.