The development of the so-called Görtler vortices in boundary layer flows over a concave surface leads to strong velocity gradients in both wall-normal and spanwise directions. This determines the flow structures to become more prone to secondary instabilities, which prompt to an early transition from laminar to turbulent flow, ultimately increasing the frictional drag. It is possible to circumvent these secondary instabilities by means of passive or active flow control strategies, with sensors and actuators implemented at the wall. In this thesis, the effect of wall cooling and heating on Görtler vortices developing in supersonic and hypersonic boundary layers is investigated from a numerical point of view. The wall temperature is imposed through a ramping function that decreases or increases an upstream base wall temperature in the streamwise direction. The results show that this type of wall cooling or heating has a mild (adverse) effect on the vortex energy, and a more considerable (but beneficial) effect on the wall shear stress.

Görtler vortices develop in boundary layer flows over concave surfaces due to the imbalance between centrifugal forces and the wall-normal pressure gradient. These vortices can be efficient precursors to transition in boundary layers exposed to free-stream disturbance or surface non-uniformities, because they can alter the mean flow causing the laminar flow to breakdown into turbulence. In this thesis, a control technique aimed at reducing the energy associated with Görtler vortices that develop in supersonic boundary layers is introduced and tested. The control algorithm is based on distributed blowing and suction, with sensors placed either in the flow or at the wall. The result show that there is a dependence between the efficiency of the control algorithm and the spanwise separation of the vortices, that is the energy reduction is more significant for larger spanwise separations. The efficiency of the control algorithm seems to be insensitive to the variation of the Mach number.

Görtler vortices evolve in boundary layers over concave surfaces as a result of the imbalance between centrifugal forces and radial pressure gradients. Depending on various geometrical and free-stream flow conditions, these instabilities may lead to secondary instabilities and early transition to turbulence. In this thesis, a control algorithm based on the boundary region equations is applied to reduce the strength of the Görtler instabilities by controlling the energy of the fully developed vortices, using either local wall deformations or blowing/suction at the wall. A proportional-integral control scheme is utilized to deform the wall or to provide transpiration velocity, where the inputs are either the wall-normal or streamwise velocity components in a plane that is parallel to the wall. The results show that the control based on wall deformation using wall-normal velocity components is more effective in tempering the vortex during its streamwise growth by almost one or two orders of magnitude.

High-amplitude freestream turbulence and surface roughness elements can excite a laminar boundary-layer flow sufficiently enough to cause streamwise-oriented vortices to develop. These vortices resemble elongated streaks having alternate spanwise variations of the streamwise velocity. Following the transient growth phase, the fully developed vortex structures downstream undergo an inviscid secondary instability mechanism and, ultimately, transition to turbulence. This mechanism becomes much more complicated in high-speed boundary layer flows due to compressibility and thermal effects, which become more significant for higher Mach numbers. In this research, we formulate and test an optimal control algorithm to suppress the growth rate of the aforementioned streamwise vortex system. The derivation of the optimal control algorithm follows two stages. In the first stage, to optimize the computational cost of the analysis, the study develops an efficient numerical algorithm based on the nonlinear boundary region equations (NBREs), a reduced form of the compressible Navier-Stokes equations in a high-Reynolds-number asymptotic framework. The NBREs algorithm results agree well with direct numerical simulation (DNS) results. The numerical simulations are substantially less computationally costly than a full DNS and have a more rigorous theoretical foundation than parabolized stability equation (PSE) based models. The substantial reduction in computational time required to predict the full development of a G\"{o}rtler vortex system in high-speed flows allows investigation into feedback control in reasonable total computational time, which is the focus of the second part of the study. In the second stage, the method of Lagrange multipliers is utilized -- via an appropriate transformation of the original constrained optimization problem into an unconstrained form -- to obtain the adjoint compressible boundary-region equations (ACBREs) and corresponding optimality conditions, which constitute the basis of the optimal control approach. Numerical solutions for high-supersonic and hypersonic flows reveal a significant decrease in the kinetic energy and wall shear stress for all configurations considered. Streamwise velocity contour plots illustrate the qualitative effect of the optimal control iterations, demonstrating a significant decrease in the amplitude of the primary vortex instabilities.

Ce travail porte sur la caractérisation de la dynamique intra-cavitaire en cavité ouverte, dont il existe peu d'études expérimentales, ainsi que sur le contrôle de cet écoulement. Nous avons réalisé une étude paramétrique des régimes primaire et secondaire de l'écoulement dans lesquels des structures de type Taylor-Görtler apparaissent. Nous avons identifié les seuils de bifurcation et montré leur nature systématiquement supercritique. Nous avons également montré que différentes familles de modes propagatifs ou stationnaires pouvaient être sélectionnées en fonction de la géométrie de la cavité. Cela a confirmé des prédictions réalisées dans des analyses de stabilité linéaire de l'écoulement de base. Nous avons montré que le régime secondaire résulte de la superposition d'ondes propagatives gauche et droite. Une tentative d'identification des coefficients des équations complexes couplées de Ginzburg-Landau décrivant cette dynamique a été conduite mais la sensibilité des coefficients à de multiples paramètres n'a pas permis d'obtenir des coefficients physiquement acceptables. Un forçage des oscillations de la couche cisaillée a été entrepris à l'aide d'un actionneur plasma froid à décharge à barrière diélectrique placé en amont de la cavité. L'analyse de la réponse de l'écoulement à un forçage périodique d'amplitude variable a permis d'identifier des plages d'accrochage en fréquence. Enfin, nous avons réalisé un contrôle en boucle fermée des oscillations de la couche cisaillée à l'aide d'une loi de contrôle à retard proposée par Pyragas dans le cadre des systèmes dynamiques chaotiques. / This work is devoted to the characterization of the dynamic inside an open cavity flow, for which few experimental studies exist. A control of the flow has been also investigated. We have performed a parametric study of the first and second regime of the flow for which Taylor-Görtler vortices type appear. Bifurcation thresholds have been identified and their systematic supercritical nature has been highlighted. We have also showed that different family of propagating or stationary modes can be selected depending on the geometry of the cavity. This has confirmed predictions obtained by linear stability analysis of the base flow in the literature. We have showed that the second regime results from the superposition of left and right propagating waves. We attended to identify coefficients of the complex coupled Ginzburg-Landau equations that describe the dynamics but the values of the coefficients are sensitive to multiple parameters. A control of oscillations of the shear layer has been achieved by the mean of a plasma actuator with dielectric barrier discharge located upstream of the cavity. Locked regimes have been identified by the analysis of the flow response to a periodic perturbation with different amplitude. Finally, we have performed a closed loop control of the oscillations of the shear layer using a delay feedback control law proposed by Pyragas in the context of chaotic dynamical systems.

Escoamentos sobre superfícies côncovas estão sujeitos à instabilidade centrífuga, dando origem a vórtices longitudinais, conhecidos como vórtices de Görtler. Esses vórtices são responsáveis por gerar distorções fortes nos perfis de velocidade. Como os vórtices são contra-rotativos, duas regiões surgem entre os mesmos: uma região de upwash e uma região de downwash. Na região de upwash o fluido próximo à parede é jogado para longe da mesma. Na região de downwash acontece o contrário, o fluido que se desloca a uma velocidade maior é jogado em direção à parede. Os vórtices se amplificam inicialmente de forma linear. À jusante na região não linear de desenvolvimento dos vórtices, a amplitude dos mesmos já é elevada, e há a formação de uma estrutura do tipo cogumelo com a distribuição da componente de velocidade na direção principal do escoamento . Essa nova distribuição de velocidade é tridimensional e difere em muito da camada limite obtida com a solução das equações de Blasius. Levando-se em consideração a camada limite térmica, já foi observado que, na média, há um aumento de transferência de calor na direção da parede. No presente trabalho, é verificado numericamente a transferência de calor na presença de vórtices de Görtler. Para tal, foi desenvolvido e implementado um código de simulação numérica direta espacial (DNS - do inglês Direct Numerical Simulation). Os resultados deste trabalho mostram a intensificação da transferência de calor através dos vórtices de Görtler, tanto no regime não-linear como na instabilidade secundária / Flows over concave surfaces are subject to centrifugal instability. It gives rise to stramwise vortices known as Görtler vortices. These vortices are responsible for generating strong distortions in the velocity profiles. As the vortices are counterrotating, two regions arise between them: a region of uowash and a region of downwash. In the upwash region, the fluid near the wall is convected away from it. In the downwash region the opposite happens, the fluid moving at a faster speed is moved towards the wall. The vortices initially amplify linearly in the downstream. When their amplitude is already high, in the non-linear development region, a mushroom-type structure, with the velocity distribution in the main flow direction, is formed. This new three-dimensional velocity distribution is different from the boundary layer obtained with the solution of Blasius equations. Taking into account a thermal boundary layer, on average, an increase in the heat transfer in the wall direction has been observed. In the present work, it is verified numerically the heat transfer in the presence of Görtler vortices. A simulation code was developed and implemanted usin Direct Numerical Simulation (DNS). The results of this work show the intensification of heat transfer through the Görtler vortices both in the non-linear regime and in the secondary instability

Nos escoamentos em turbomáquinas temos como principais características a tridimensionalidade, possível ocorrência de separação da camada limite, relaminarização, transição laminar-turbulenta, dentre outros efeitos físicos. De acordo com alguns estudos experimentais em turbinas observouse que a transição laminar-turbulenta pode se estender por até 60% da corda de uma pá de turbina. Uma boa estimativa para se prever corretamente o local da transição é indispensável para que seja obtida uma melhoria na eficiência das turbinas. Escoamentos sobre superfícies côncavas estão sujeitos à instabilidade centrífuga, podendo dar origem a vórtices longitudinais, conhecidos como vórtices de Görtler. Esses vórtices são responsáveis por gerar distorções fortes nos perfis de velocidade e consequentemente nos perfis de temperatura. O presente estudo tem por objetivo estudar a influência da variação da curvatura de uma superfície côncava, e os efeitos do comprimento de onda transversal no processo de transição, e sua influência nas taxas de transferência de calor. Para tal, um código de simulação numérica paralelizado, com alta ordem de precisão, foi utilizado para resolver numericamente as equações de Navier-Stokes. Este código é validado através de comparações entre resultados obtidos com uso da teoria de estabilidade linear, e com resultados de simulações numéricas não lineares. Resultados obtidos evidenciam a influência da variação da curvatura, e os efeitos causados pelo comprimento de onda transversal nas instabilidades de Görtler, e secundária. Tais evidências comprovam que a variação da curvatura pode ser útil no controle do processo de transição laminar-turbulenta, e que as taxas de transferência de calor de um escoamento de Görtler desenvolvido em superfícies de curvatura variável podem ser intensificadas, atingindo valores superiores aos obtidos em escoamentos turbulentos. / Some characteristics of flows in turbomachinery are the three-dimensionality, possible occurrence of separation of the boundary layer, relaminarization, laminar-turbulent transition, among other physical effects. According to some experimental observations in turbines, it has been observed that the laminar-turbulent transition can extend over 60% chord of a turbine blade. A good estimate to correctly predict the location of the transition is essential for an improvement in the efficiency of turbines. Flow over concave surfaces is subjected to centrifugal instability, which may lead to formation of longitudinal vortices, known as the Görtler vortices. These vortices are responsible for generating strong distortions in the velocity profiles and hence the temperature profiles. The current goal aims to study the influence of the curvature variation of a concave surface and the effects of spanwise wavelength on the transition process and its influence on the heat transfer rates. For this, a parallel numerical simulation code, with a high order of precision, was used to numerically solve the Navier-Stokes equations. This code is validated through comparisons between results obtained using linear stability theory, and nonlinear numerical simulations results. Results obtained show the influence of the curvature variation, and the effects caused by the spanwise wavelength on the Görtler and secondary instabilities. This evidence proves that the curvature variation can be useful in the control of the laminar-turbulent transition process, and that heat transfer rates of a Görtler flow developed on variable curvature surfaces can be intensified, and reach values higher than these achieved in turbulent flows.

Escoamentos sobre superfícies côncovas estão sujeitos à instabilidade centrífuga, dando origem a vórtices longitudinais, conhecidos como vórtices de Görtler. Esses vórtices são responsáveis por gerar distorções fortes nos perfis de velocidade. Como os vórtices são contra-rotativos, duas regiões surgem entre os mesmos: uma região de upwash e uma região de downwash. Na região de upwash o fluido próximo à parede é jogado para longe da mesma. Na região de downwash acontece o contrário, o fluido que se desloca a uma velocidade maior é jogado em direção à parede. Os vórtices se amplificam inicialmente de forma linear. À jusante na região não linear de desenvolvimento dos vórtices, a amplitude dos mesmos já é elevada, e há a formação de uma estrutura do tipo cogumelo com a distribuição da componente de velocidade na direção principal do escoamento . Essa nova distribuição de velocidade é tridimensional e difere em muito da camada limite obtida com a solução das equações de Blasius. Levando-se em consideração a camada limite térmica, já foi observado que, na média, há um aumento de transferência de calor na direção da parede. No presente trabalho, é verificado numericamente a transferência de calor na presença de vórtices de Görtler. Para tal, foi desenvolvido e implementado um código de simulação numérica direta espacial (DNS - do inglês Direct Numerical Simulation). Os resultados deste trabalho mostram a intensificação da transferência de calor através dos vórtices de Görtler, tanto no regime não-linear como na instabilidade secundária / Flows over concave surfaces are subject to centrifugal instability. It gives rise to stramwise vortices known as Görtler vortices. These vortices are responsible for generating strong distortions in the velocity profiles. As the vortices are counterrotating, two regions arise between them: a region of uowash and a region of downwash. In the upwash region, the fluid near the wall is convected away from it. In the downwash region the opposite happens, the fluid moving at a faster speed is moved towards the wall. The vortices initially amplify linearly in the downstream. When their amplitude is already high, in the non-linear development region, a mushroom-type structure, with the velocity distribution in the main flow direction, is formed. This new three-dimensional velocity distribution is different from the boundary layer obtained with the solution of Blasius equations. Taking into account a thermal boundary layer, on average, an increase in the heat transfer in the wall direction has been observed. In the present work, it is verified numerically the heat transfer in the presence of Görtler vortices. A simulation code was developed and implemanted usin Direct Numerical Simulation (DNS). The results of this work show the intensification of heat transfer through the Görtler vortices both in the non-linear regime and in the secondary instability

The purpose of a multi-element high lift device is to increase lift dramatically while controlling the stall limit. The fluid flow over a multi-element high lift device has been explored widely both experimentally and numerically at high Reynolds numbers (O(10^6 )). The numerical simulations use turbulence models and hence details of the flow are not yet available. Low Reynolds number (O(10^4 )) flows over high lift devices have not been explored until recently. These lower Reynolds number flows have applications in the development of small aerial vehicles. The present work discusses both two-dimensional and three-dimensional direct numer- ical simulations of fluid flow over a 30P30N three-element high lift system using a high-order spectral element method code, Nek5000, that solves the incompressible Navier-Stokes equations. The intricate geometry of the multi-element device poses a challenge for the high-order spectral element method. We study the complex flow physics in the slat cove region and the wake/shear layer interaction over a 30P30N three-element high lift device. The targeted cases are at Reynolds num- bers based on stowed chord lengths (Rec ) of 8.32 × 10^3 , 1.27 × 10^4 , and 1.83 × 10^4 at angle of attack of 4. A critical interval for Rec has previously been found between 1.27 × 10^4 and 1.38 × 10^4 in experiments. This divides the flow into two types: when Rec is below the critical interval, no roll-up is observed in the slat cove and Görtler vortices dominate the slat wake; however when the Rec is above the critical interval, a roll-up is observed in the slat cove and co-existence of streamwise and spanwise vortices is confirmed in the slat wake. We confirm the presence of the critical interval from the simulations performed at three values of Rec . Lift and drag analysis is provided along with pressure coefficient plots for each element of the multi-element airfoil. Different vortical structures are also identified in the transition of flow from two dimensions to three dimensions. The relevant validation is performed with the available experimental data.

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.