This thesis considers the advantages of incorporating laminar kinetic energy modelling into turbulence modelling, in order to predict laminar-turbulent transition. The final aim is to implement an improved transition model into the industrial Finite-Volume code, Code Saturne. The literature review suggests that in order for a RANS-based model to predict transition accurately, modelling of complex, anisotropic phenomena is necessary. The Walters-Cokljat model is shown to compare very well to other transition modelling methods, including correlation-based modelling. The Walters-Cokljat model is a single-point RANS-based model that solves an additional transport equation for laminar kinetic energy. This transition model is especially desirable from an industrial stand-point, due to its single-point RANS basis, with only 3 transport equations. Although this method shows great promise as an industrial tool for transition prediction, results presented here show that there are aspects of the model that require modification. The definition of effective length-scale and the method of accounting for the effects of shear sheltering are the two main areas for consideration. The current definition of effective length-scale is found to be inappropriate for flows with large free-stream length-scales, which are common-place in turbomachinery applications. Another phenomenon commonly found in turbomachinery is separation-induced transition; however, the current function for shear sheltering effects inhibits transition when turbulence intensity is not the forcing factor. Additionally, when reviewed analytically, the definition and placement of the shear sheltering function does not match the observations of Jacobs and Durbin. Alternatives for the definitions of the effective length-scale and the shear sheltering function are proposed. The individual proposals are tested, and steps towards a full working implementation are documented.

No description available.

Temporally developing direct numerical simulations (T-DNS) are performed and validated for bypass transition of a zero pressure gradient flat plate boundary layer to understand the interplay between pressure-strain terms and flow instability mechanisms, and to propose and validate a phenomenological hypothesis for the identification of a robust transition onset marker for use in transition-sensitive Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) simulations. Results show that transition initiates at a location where the slow pressure-strain term becomes more dominant than the rapid term in the pre-transitional boundary layer region. The slow pressure strain term is responsible for the transfer of turbulence energy from the streamwise component to other components while the rapid pressure strain term counteracts with the slow term in the pre-transitional regime before transition onset akin to a shear sheltering like effect. The relative magnitudes of the slow and rapid terms thus provide a basis for the development of physically meaningful large-scale parameters that can be used as a transition onset marker for Reynolds averaged Navier-Stokes RANS simulations.

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>

Le passage du régime laminaire au régime turbulent s’accompagne d’importantes modifications des propriétés physiques de l’écoulement. Une prévision précise du point du début de la transition laminaire/turbulent revêt donc une importance considérable dans de nombreux domaines pratiques. Lorsque l’intensité des perturbations extérieures est significative, c'est-à-dire dans le cas de couches limites se développant sur une paroi présentant des rugosités ou soumises à une forte turbulence résiduelle (sillage impactant), les mécanismes de formation et d’amplification des instabilités sont profondément modifiés. Ces perturbations sont les modes de Klebanoff (également appelés stries) qui s’amplifient et déclenchent la transition, qualifiée dans ce cas de Bypass. Ces stries sont très énergétiques, caractérisées par des fluctuations de vitesse très importantes (de l’ordre de 10% de la vitesse extérieure), alors que la couche limite conserve son caractère laminaire. La thèse proposée concerne la modélisation de ces stries via la résolution d’une équation de transport pour l’énergie cinétique dite laminaire. Dans un premier temps, le travail du candidat portera sur la modélisation des termes de production et de dissipation de l’énergie cinétique laminaire. Ceux-ci sont liés au processus de réceptivité de la couche limite vis-à-vis des perturbations extérieures et à la dynamique des modes de Klebanoff dans la zone laminaire. Pour ce faire, la thèse s’appuiera sur des études réalisées depuis plusieurs années au sein de l’unité ITAC sur la théorie des perturbations optimales ainsi que sur les travaux numériques et expérimentaux prévus dans le cadre d’un projet de recherche interne Onera. Classiquement cette équation de transport est couplée avec celles correspondant à l’énergie cinétique turbulente et à la dissipation, le mécanisme d’échange entre les énergies cinétiques laminaire et turbulente devra être soigneusement étudié : ce dernier pilote la transition vers la turbulence. Une attention particulière sera portée aux couches limites décollées et plus précisément à la prise en compte de la transition dans ces bulbes. Cette nouvelle modélisation innovante permettra l’amélioration d’une première approche pour le calcul de la transition bypass dans le solveur elsA, développé à l’Onera, et constituera une étape importante vers la mise en place de techniques de prévision de la transition pratiques et performantes. / This work aims to develop a new bypass-transition prediction model based on the Klebanoff modes dynamics. To represent these mode dynamics the Laminar Kinetic Energy (LKE) concept has been chosen, in order to model these mode energy with a new variable. A new deffinition is given to the LKE and a transport equation consequently derived to describe the Klebanoff modes growth and destabilisation. This equation is incorporated in a k-omega turbulence model as done by Walters & Cokljat, to give a three-equation kL-kT-omega formulation. This new model is written in a Reynolds-averaged Navier-Stokes (RANS) pattern and only uses local variables, it thus can be used in an industrial context.

The initial slow viscous growth of the Tollmein-Schlichting wave in a canonical boundary layer transition is absent in bypass and wake-induced transitions. Although there have been a great deal of studies pertaining to bypass transition in boundary layers, the underlying breakdown mechanism is not clearly understood and it continues to be a subject of interest. Similarly, a wake-induced transition caused by Karman wake in the freestream remains poorly understood. The breakdown in this case is caused by anisotropic disturbances containing large scale unsteadiness in the freestream. Differing view points among workers on the transition process have also added to the complexities. In this thesis, bypass and wake-induced boundary layer transitions studied experimentally towards understanding various flow breakdown features are reported. The measurements were made on a flat plate boundary layer in a low-speed wind tunnel. The particle image velocimetry (PIV) technique was extensively used. Various grids were used to generate nearly isotropic freestream turbulence. A circular cylinder was placed at different heights from the plate leading edge to generate Karman wake in the freestream. Two cylinders of different diameters were used to vary the Reynolds number(based on the cylinder diameter). The PIV measurements being simultaneous over a large spatial domain enabled to assess various spatial transitional flow structures. In the case of bypass transition, the streamwise velocity fluctuation, u, is found to exhibit some organized negative and positive fluctuations that dominate the flow during transition, and confirm the simulation results reported in the literature. These positive and negative u fluctuations are found to be associated with the streak unsteadiness. By conditional sampling of these positive and negative u fluctuations, we find that urms (root-mean-squaredof u)can be expressed as a linear combination of urms,f and urms,b,i.e. urms = a(urms,f + urms,b); ais constant, and the subscripts fand bdenote the positive and nega-tive ufluctuations, respectively. Both urms,f and urms,b arefoundto follow the non-modal growth distribution. The wall-normal results clearly show that an inclined shear layer is often associated with an organized structure of negative ufluctuations and an inflectional in-stantaneous velocity profile. These inclined shear layers appear to be similar to those in ribbon-induced transition. The turbulent spot precursor appears to be the vortex shedding from an oscillating in-clined shear layer. Interestingly, the normalized vortex shedding fre-quency is found to be Reynolds number invariant, as in the case of ribbon-induced transition. The present study also confirms the sim-ulated turbulent spot features, including a thin log-law at the break-down stage. The spanwise plane PIV results reveal the signature of streak secondary instability in the flow in terms of symmetric and anti-symmetric streaks oscillations. The initial growth of streak amplitude is followed by a slow decay. The maximum streak amplitude is well above30% of the freestream velocity. These two aspects provide support to the streak instability analysis reported in the literature. While the present wake-induced transition study provides some sup-port to the available numerical simulation and experimental results, some new results have also emerged. The measured sharp rise in the disturbance energy during transition is found to be closer to the simulated result, compared to the difference reported in the literature. The spanwise vortices in the early stage, as also seen in other experimental studies, deform leading to the formation of lambda structures, the signature of which is found by the linear stochastic analysis. With increased Reynolds number and decreased cylinder height from the plate, the physical size of the lambda structure is found to decrease. These lambda structures are often found to appear in a staggered manner in the spanwise plane, as in the case of sub-harmonic boundary layer transition. Although a sub-harmonic peak in the frequency spectra is reported in the literature, as also in the present study, the clear staggered pattern went unnoticed. Streamwise streaks are subsequently generated due to the mean shear stretching of these lambda vortices. The spanwise spacing of these streamwise streaks is found to be comparable with the recent simulation results. Also, these streaks are found to undergo somewhat sinuous-like oscillations, compared to the only varicose type oscillations reported in the literature. The streak amplitude is found to saturate at about 35% of the freestream speed. Here again an inclined shear layer in the wall-normal plane is associated with organized negative u fluctuations and an inflectional instantaneous velocity profile. The movement of the peak urms towards the wall is found to be due to the positive u fluctuation, which follows a hairpin-like structure. The inclined shear layers herein are associated with the lambda or a hairpin-like structure. As in a by-pass transition, an inclined shear layer, vortex shedding from it, the imprint of which is also found in the linear stochastic analysis are present. The normalized high frequency shed vortices is found to be Reynolds number invariant in the present wake-induced transition, as in ribbon-induced and bypass transitions. Compared to the re-cent suggestion that the parent-offspring mechanism is the governing self-sustaining mechanism in the boundary layer, the present study suggests that streak-instability mechanism is also present. The proper orthogonal decomposition(POD) analysis of the experimental data was carried out with an emphasis on the bypass transition case studied. The first few energetic POD modes are found to capture the dominant flow structures, i.e. the organized positive and negative u fluctuations. In the case of bypass transition, the first two energetic POD modes are self-similar, i.e. independent of the freestream turbulent intensity and the Reynolds number. An attempt is also made to construct a low-dimensional model with the POD eigenfunction modes to predict the qualitative dynamics of bypass transition. This has revealed the existence of a traveling disturbance in the bypass transition. On the whole, the present study shows some similar breakdown features in bypass and wake-induced transitions, although more studies in this regard are essential.

The transition from a laminar to a turbulent flow strongly modifies the boundary layer properties.Understanding the mechanisms leading to transition is crucial to reliably predict aerodynamicperformances. For boundary layers subjected to high levels of external disturbances, the naturaltransition due to the amplification of the least stable mode is replaced by an early transition, calledBypass transition. This is the result of non-normal mode interactions that lead to a phenomenon oftransient growth of disturbances. These disturbances are known as Klebanoff modes and take theform of streamwise velocity streaks.This thesis aims at understanding this linear mechanism of transient growth and quantifying itsinfluence on the classical modal amplification of disturbances. This is done by computing theso-called optimal perturbations, i.e. the initial disturbances that undergo maximum amplificationin the boundary layer.These optimal perturbations are first determined for two-dimensional compressible boundary layersdeveloping over curved surfaces. In particular, we show that Klebanoff modes naturally evolvetowards Görtler vortices that occur over concave walls. Three-dimensional boundary layers arethen considered. In such configurations, transient growth provides an initial amplitude to crossflowvortices. Finally, applying the tools developed in this thesis to new flow cases such as swept wingsprovides further understanding of the phenomenon of transient growth for realistic geometries. / Le passage du régime laminaire au régime turbulent s’accompagne d’importantes modifications despropriétés physiques de la couche limite. La détermination précise de la transition est donc crucialedans de nombreux cas pratiques. Lorsque la couche limite se développe dans un environnementextérieur faiblement perturbé, la transition est gouvernée par l’amplification du mode propre le moinsstable. Lorsque l’intensité des perturbations extérieures augmente, des interactions multimodalesentraînent une amplification transitoire des perturbations. Ce phénomène peut conduire à unetransition prématurée, appelée transition Bypass. Les perturbations prennent alors la forme destries longitudinales de vitesse appelées modes de Klebanoff.L’objectif de cette thèse est d’étudier ce mécanisme linéaire de croissance transitoire et soninfluence sur l’amplification modale classique des perturbations. Cela passe par la déterminationdes perturbations les plus amplifiées au sein de la couche limite, appelées perturbations optimales.Ces perturbations optimales sont d’abord calculées pour des couches limites bidimensionnelles etcompressibles se développant sur des surfaces courbes. En particulier, on montre que les modes deKlebanoff évoluent vers les tourbillons de Görtler qui se forment sur des parois concaves. Le cas plusgénéral de couches limites tridimensionnelles est ensuite envisagé. Pour de telles configurations, lacroissance transitoire fournit une amplitude initiale aux instabilités transversales. Enfin, l’applicationdes outils développés dans cette thèse fournit de nouveaux résultats pour des cas d’écoulementsautour de géométries réalistes comme une aile en flèche.