Le présent travail s'inscrit dans le cadre d'une convention CIFRE entre Liebherr-Aerospace Toulouse SAS et le Département d'Aérodynamique Energétique et Propulsion (DAEP) de l'ISAE. L’étude porte sur l'analyse des mécanismes responsables de la perte de la stabilité d'un étage de compression centrifuge. L'objectif industriel sous-jacent est l'élargissement de la plage de fonctionnement stable des compresseurs. Les travaux sont abordés par voie numérique à l'aide du code de calcul elsA, développé par l'ONERA et le CERFACS. Les simulations instationnaires sont effectuées sur la circonférence complète des roues.La topologie de l'écoulement est analysée dans un premier temps selon trois points de fonctionnement stables répartis sur la courbe caractéristique de l'iso-vitesse de design. Cette analyse permet de décrire l'évolution des phénomènes lorsque le débit est réduit. Au voisinage de la ligne de stabilité, l'excès d'incidence sur les aubes du rouet déclenche le décollement de la couche limite sur la face en d’expression. Le guide issu de la zone décollée migre en direction du carter et alimente l'écoulement de jeu. En conséquence, le déficit de vitesse au voisinage du carter s'intensifie et un lâcher périodique de structures tourbillonnaires apparaît à l'interface entre les écoulements secondaires et l'écoulement principal.La perte de la stabilité de l'étage s’établit suite à la naissance d'une perturbation de type « modal » au sein de l'espace lisse. Cette dernière induit de fortes distorsions circonférentielles dans l'étage de compression mais affecte plus particulièrement l'écoulement à l'entrée du rouet. Sur la moitié de la circonférence, l'interface entre les écoulements secondaires et l'écoulement principal ainsi que les structures tourbillonnaires sont déplacées en amont du front de grille. Le taux total-total du rouet chute de manière brutale et entraîne la perte de la stabilité de l'étage.Enfin, la dernière partie de ce travail est dédiée à la mise en place de critères adaptés au modèle stationnaire mono-canal utilisé chez LTS. Ces critères seront utilisés pour améliorer la prédiction de la ligne de pompage en phase de design. Dans un premier temps, la pertinence de l'utilisation du modèle stationnaire au voisinage de la ligne de stabilité est évaluée. Dans un second temps, les influences de la vitesse de rotation et de la géométrie de l’étage sur la topologie sont étudiées.Deux situations sont jugées critiques vis-à-vis de la stabilité. La première concerne l'alignement de l'interface entre l'écoulement de jeu et l'écoulement principal avec le front de grille. La deuxième concerne l'opération du compresseur avec un angle d'écoulement en sur-incidence sur les aubes du diffuseur qui s'établit sur toute la hauteur de veine. / This work results from a CIFRE partnership between Liebherr-Aerospace Toulouse SAS and the Aerodynamics, Energetics and Propulsion Department (DAEP) of ISAE. The main objective is to investigate the mechanisms responsible of the stall onset in a centrifugal compressor operating at the nominal rotational speed. It is part of a larger work which aims at extending the stable operating range of compressors integrated in air conditioning system. The analyses are based on the results of unsteady simulations, in a calculation domain comprising all the blade passages. They are performed with the elsA software developed by ONERA and CERFACS.The investigations show the modifications of the unsteady flow pattern when the mass flow is reduced along the speed line. Near the stability limit, the high incidence angle on the impeller blade leads to a boundary layer separation on the suction side. The fluid in the separation zone moves toward the shroud and enlarges the low momentum flow zone generated by the leakage flow. The interface between the leakage flow and the main flow becomes unstable to the extent of a periodic vortex formation.The path to instability is driven by the growth of a small amplitude disturbance (modal wave) rotating in the vaneless space. The length scale of the wave is equal to the compressor circumference. This perturbation induces distortions and alters the flow characteristics in every location of this subsonic stage, and more specifically the impeller inlet flow structure: the unstable interface between the main flow and the leakage flow is periodically moved upstream of the leading edge plane causing a significant drop of the impeller total-to-total pressure ratio. The last part of this work concerns the definition of criteria which can improve the surge line prediction during the design process in an industrial environment. Therefore, they are adapted to the numerical steady model using the mixing plane approach. To do so, the capacity of the steady model to predict the flow structure when the compressor operates near stall is investigated. Then, the effects of the rotational speed and of the compressor geometry are evaluated. Theses two steps have permitted to define two critical situations regarding the stage stability. The first one is related to the alignment of the interface between the main flow and the leakage flow with the leading edge plane. The second one concerns the compressor operation with positive incidence on the diffuser vane, along the full span.

L’écoulement dans les turbomachines est tri-dimensionnel et instationnaire. Actuellement, les concepteurs de moteurs cherchent à réduire l’encombrement et le poids des machines. En conséquence, les interactions entre les roues, appelées interactions rotor-stator, sont renforcées. Parmi elles, l’effet potentiel remontant n’est désormais plus négligeable malgré sa rapide atténuation spatiale. Dans cette étude, cet effet potentiel remontant a été analysé sur une configuration spécialement conçue : une grille linéaire d’aubes de turbine, suivie de barreaux défilants en aval à une distance de 20% de corde axiale, simulant des aubes de rotor en aval. La grande échelle du banc d’essais facilite l’étude du comportement de la couche limite des aubes de la grille. Des mesures de pression et d’anémométrie laser à deux composantes, synchronisées avec le défilement des barreaux aval sont réalisées. Le nombre de Reynolds, basé sur la corde, est 1.6 × 105. Une grille de turbulence placée en amont de la grille d’Aube afin de pouvoir augmenter le taux de turbulence amont a été utilisée. Des résultats de mesures en absence de cette grille (faible taux de turbulence amont) sont également présentés et analysés. Une modélisation numérique, basée sur un calcul laminaire avec un préconditionnement basse vitesse pour la même configuration, a montré la déformation des lignes de courant de l’écoulement dans le canal inter-aubes, en fonction de la position du barreau aval. La distribution de pression autour de l’aube est également périodiquement modifiée. Les résultats stationnaires expérimentaux, en absence de tout barreau aval, ont révélé un décollement de la couche limite à l’extrados de l’aube à bas taux de turbulence amont(Tu−am = 1.2%) qui est supprimé à haut taux de turbulence amont (Tu−am = 4.2%) ; la couche limite commence alors sa transition par un mode by-pass. Aucun effet instationnaire dans la couche limite n’a été observé à l’intrados, quel que soit le taux de turbulence amont. L’étude instationnaire, avec le défilement des barreaux en aval, a permis de mettre en évidence un décollement périodique de la couche limite à l’extrados à bas taux de turbulence amont (Tu−am = 1.8%). Dans ce cas, la couche limite suit deux modes de transition au cours d’une période : une transition par décollement et une transition bypass. Au contraire, dans le cas à fort taux de turbulence amont (Tu−am = 4.0%), aucun décollement de la couche limite n’a été décelé. La couche limite est sujette à l’effet instationnaire à l’extrados. Elle est devenue turbulente au bord de fuite à tout instant par un mode by-pass. Cette étude a montré que l’effet potentiel issu d’un roue en aval est du même ordre de grandeur que les effets de sillage et doit être pris en compte dans l’analyse des phénomènes. Par des méthodes d’indexation de roues, le décollement de la couche limite pourrait être supprimé. / Turbomachinery designers wish to reduce the size and weight of engines. One way of achieving this is by reducing the distance between rotor and stator elements. In doing so, the rotor-stator interaction becomes more significant. In particular, the long-range influence of pressure potential is no longer negligible, and affects both upstream and downstream flow. Previously, only downstream interactions of blade wakes were considered important. Here we examine the upstream potential effect generated by downstream moving cylindrical rods on an upstream low pressure turbine blade. A large scale rectilinear blade cascade was constructed to improve access to the boundary layer. The Reynolds number, based on the chord, was 1.6 × 105. Pressure measurements and two-dimensional Laser Doppler Anemometry around the blade were performed to study the boundary layer behavior. Recorded data points are phase averaged with the downstream moving cylindrical rods. A grid is placed upstream of the blade cascade to increase the inlet turbulence intensity. A numerical investigation, based on a laminar simulation with low velocity preconditioning method was carried out on the same configuration. The flow streamlines and the pressure distribution around the blade were found to depend strongly on the downstream rod position. No unsteady effects in the boundary layer of the pressure side were observed, for the inlet turbulence intensities used in our study. Steady experimental results revealed a boundary layer separation bubble on the blade suction side at a low turbulence intensity (Tu−in = 1.2%), whereas the boundary layer became turbulent via by-pass transition at a higher turbulence intensity (Tu−in = 4.2%).It is seen that, in the unsteady configuration, at a low turbulence intensity (Tu−in =1.8%), the laminar boundary layer experiences separation once per rod period. Two transition modes were identified that alternate during a rod period : a separation transition mode and a by-pass mode, which were conditioned by the downstream rod position. Ata higher turbulence intensity (Tu−in = 4.0%), no boundary layer separation occurred thereby following a bypass transition mode during an entire rod period. The experimental results presented here demonstrate the large influence of the downstream potential effect generated by a downstream row on the upstream blade boundary layer behavior. In order to improve the efficiency of engines, this effect and its interaction with the wake effect must be taken into account in turbomachinery design.

L'estimation précise de la traînée est aujourd'hui un enjeu majeur pour les avionneurs. Il est nécessaire d'identifier et de quantifier ses sources phénoménologiques dans le cadre d'un processus de design efficace. Les méthodes champ lointain, qui permettent une telle décomposition, sont cependant limitées aux applications stationnaires.Cette étude consiste à développer une méthode d'extraction champ lointain destinée à permettre une décomposition phénoménologique de la traînée pour des écoulements instationnaires. La première étape a consisté à généraliser la formulation stationnaire de Van der Vooren aux écoulements instationnaires. Des axes pour l'amélioration de la robustesse et du contenu physique ont ensuite été explorés, avec la mise en évidence de phénomènes acoustiques. La formulation ainsi obtenue a ensuite été appliquée à des cas tests simples, dans le but de valider la décomposition phénoménologique. Le comportement des composantes de traînée s'est avéré cohérent avec la physique de l'écoulement. Enfin, la méthode a été appliquée à des cas complexes afin de démontrer ses capacités : un cas instationnaire 3D ainsi qu'un écoulement simulé en ZDES.Dans l'avenir, il serait intéressant de continuer à explorer la définition de la composante de traînée induite, par exemple en utilisant les formulations basées sur le vecteur vitesse. En ce qui concerne les cas d'application, l'évaluation de la performance d'un doublet d'hélices contra-rotatives pourrait fortement bénéficier de l'utilisation d'une méthode comme celle-ci. Enfin, des applications en vol d'ailes battantes pourraient être d'intéressantes perspectives. / Accurate drag prediction is now of a major issue for aircraft designers. Its phenomenological sources need to be identified and quantified for an efficient design process. Far-field methods, which allow such phenomenological drag breakdown, are however restricted to steady flows. This study consists in developing a far-field drag prediction method aiming at a phenomenological breakdown of drag for unsteady flows. The first step has consisted in generalizing the steady formulation of Van der Vooren to unsteady flows, starting from a new rigorous proof. Axes for the improvement of the robustness and physical background have then been explored. Acoustic contributions have in particular been highlighted and quantified. The resulting five-components formulation has then been applied to simple cases, in order to validate as best as possible the phenomenological breakdown. The behavior of the drag components has proved to be consistent with the physics of the flow. Finally, the method has been applied to complex cases in order to demonstrate its capabilities: a 3-D case and a flow simulated by the ZDES method. In the future, it would be interesting to further improve the definition of the induced drag component, for example by using velocity-based formulations. As far as the application cases are concerned, the performance evaluation of a Counter-Rotating-Open-Rotor would strongly benefit from such a method. Unsteady optimization of one of the drag component could also be contemplated. Finally, applications in aeroelasticity or flapping flight would be an interesting perspective.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2010-2011 / L'objectif du présent mémoire est de contribuer à la caractérisation de l'écoulement dans l'aspirateur d'une turbine hydraulique de type hélice. Pour y arriver, un anémomètre laser à effet Doppler (LDV) est utilisé dans plusieurs sections de l'aspirateur d'un modèle de cette turbine, maintenu sous plusieurs conditions d'opération, pour créer une base de données de vitesse. La composante radiale en sortie de roue est particulièrement recherchée puisque son acquisition constitue un défi technique au niveau expérimental et que son importance dans l'évolution de l'écoulement est établie. Les résultats principaux font ressortir que la vitesse radiale arbore une grandeur non négligeable en sortie de roue et varie de façon importante suivant les conditions d'opération. L'asymétrie circonférentielle des trois composantes de vitesse est également confirmée dans cette section. L'évolution de l'écoulement dans le coude de l'aspirateur engendre un débalancement du débit sortant à travers les deux pertuis. Finalement, l'écoulement dans l'aspirateur de ce type de turbine est fortement instationnaire, ce qui rend essentielle l'étude des fluctuations qui s'y développent.

Les aspirateurs de turbines hydrauliques jouent un rôle crucial dans l’extraction de l’énergie disponible. Dans ce projet, les écoulements dans l’aspirateur d’une turbine de basse chute ont été simulés à l'aide de différents modèles de turbulence dont le modèle DDES, un hybride LES/RANS, qui permet de résoudre une partie du spectre turbulent. Déterminer des conditions aux limites pour ce modèle à l’entrée de l’aspirateur est un défi. Des profils d’entrée 1D axisymétriques et 2D instationnaires tenant compte des sillages et vortex induits par les aubes de la roue ont notamment été testés. Une fluctuation artificielle a également été imposée, afin d’imiter la turbulence qui existe juste après la roue. Les simulations ont été effectuées pour deux configurations d’aspirateur du projet BulbT. Pour la deuxième, plusieurs comparaisons avec des données expérimentales ont été faites pour deux conditions d'opération, à charge partielle et dans la zone de baisse rapide du rendement après le point de meilleur rendement. Cela a permis d’évaluer l'efficacité et les lacunes de la modélisation turbulente et des conditions limites à travers leurs effets sur les quantités globales et locales. Les résultats ont montrés que les structures tourbillonnaires et sillages sortant de la roue sont adéquatement résolus par les simulations DDES de l’aspirateur, en appliquant les profils instationnaires bidimensionnels et un schéma de faible dissipation pour le terme convectif. En outre, les effets de la turbulence artificielle à l'entrée de l’aspirateur ont été explorés à l'aide de l’estimation de l’intermittence du décollement, de corrélations en deux points, du spectre d'énergie et du concept de structures cohérentes lagrangiennes. Ces analyses ont montré que les détails de la dynamique de l'écoulement et de la séparation sont modifiés, ainsi que les patrons des lignes de transport à divers endroits de l’aspirateur. Cependant, les quantités globales comme le coefficient de récupération de l’aspirateur ne sont pas influencées par ces spécificités locales. / Draft tubes play a crucial role in elevating the available energy extraction of hydroturbines. In this project, turbulent flows in the draft tube of a low-head bulb turbine were simulated using, among others, an advance hybrid LES/RANS turbulent model, called DDES, which can resolve portions of the turbulent spectrum. Providing appropriate inflow boundary conditions for such models is a challenging issue. In this regard, different inflow boundary conditions were tested, including axisymmetric 1D profiles, and unsteady 2D inflow profiles that take runner blade wakes and vortices into account. Artificial fluctuation at the inlet section of the draft tube was also included to mimic the turbulence existing after the runner. Simulations were conducted for two draft tube configurations of the BulbT project. For one of them, intensive comparisons with experimental data were done for two operating conditions, one at part load and another in the sharp drop-off portion of the efficiency hill after the best efficiency point. This allowed to assess the effectiveness and shortcomings of the adopted turbulence modeling and boundary conditions through their effects on the global and local quantities. The results showed that the runner-related vortical structures and wakes are appropriately resolved using stand-alone DDES simulation of the draft tube flows. This is achieved by applying unsteady 2D inflow profiles along with adopting low dissipation scheme for the convective term. Furthermore, the effects of applying artificial turbulence at inlet were explored using separation intermittency, two-point correlation, energy spectrum and Lagrangian coherent structure concepts. These analyses revealed that the type of inflow boundary conditions modifies the details of the flow and separation dynamics as well as patterns of the transport barriers in different regions of the draft tube. However, the global quantities such as recovery coefficient are not influenced by these local features.

La problématique de surverse du noyau de barrages en remblai lors de crues causées par une défaillance des systèmes de contrôle du niveau du réservoir préoccupent les gestionnaires de ces ouvrages. Une surverse de ce type peut entraîner des vitesses d’écoulement beaucoup plus importantes que celles de conception et causer un écoulement important dans des zones où aucun écoulement n’est initialement prévu. Pour évaluer l’impact des propriétés hydrauliques des composantes du barrage, de la géométrie du barrage et des caractéristiques de la crue, un modèle numérique 2D d’éléments finis a été conçu en utilisant le logiciel COMSOL Multiphysics 5.3. Le modèle numérique a permis de réaliser une série d’analyses de sensibilité sur un barrage modélisé typique qui a montré l’importance de la perméabilité de la crête et du noyau ainsi que de la hauteur maximale atteinte par la crue sur la vitesse à laquelle arrive la surverse et sur les vitesses d’érosion atteintes. L’érosion de contact selon le critère de Brauns (1985) ne semblait pas problématique lors des analyses de sensibilité alors que l’érosion interne selon le critère de Konrad & Côté (2013) paraissait très probable. Ensuite, une étude de cas a été réalisée portant sur un barrage existant d’Hydro-Québec. En raison des perméabilités très élevées de la crête pour ce barrage, les vitesses en crête dans le cas de hautes crues peuvent causer un écoulement turbulent à l’aval du noyau, sans toutefois être problématique pour cet ouvrage au niveau de l’érosion de contact. Finalement, la convergence numérique a été grandement meilleure lors de l’utilisation de courbes de rétention basée sur le modèle de van Genuchten (1980) comparativement à celles basées sur le modèle de Brooks & Corey (1964). / The issue of core overtopping in embankment dams during floods caused by failures of water-level control systems is a concern to dam managers. Core overtopping can generate flow velocities far greater than those considered during design and can result in intense flow in zones where no flow was intended. A 2D finite-element numeric model running on the COMSOL Multiphysics 5.3 software platform was generated to evaluate the impact of soil permeability, dam geometry and flood shape on seepage and internal flow systems. This numerical model was used to conduct a series of sensitivity analyses on a conceptual standard dam. These analyses have shown the importance of crest and core permeability and maximum height reached by the flood on the rate at wich core overtopping occurs and on erosion velocities generated. Contact erosion based on the Brauns (1985) criteria did not seem to be a concern for these sensitivity analyses as opposed to internal erosion based on the Konrad & Côté (2013) criteria which seemed highly probable. Subsequently, a case-study on an existing Hydro- Québec earthfill dam was also carried out. Because of the high permeability of the crest components of this structure, significant velocities and turbulent flow could be obtained in the crest material in cases of high floods, without however being problematic in terms of contact erosion. Furthermore, numerical convergence of the model has been greatly improved by using the van Genuchten-based (1980) water retention curves instead of the Brooks & Corey-based (1964) curves.

La présente thèse propose une étude expérimentale du décollement dans le diffuseur d’un modèle de turbine hydroélectrique bulbe. Le décollement se produit quand la turbine est opérée à forte charge et il réduit la section effective de récupération du diffuseur. La diminution de la performance du diffuseur à forte charge engendre une baisse brusque de l’efficacité de la turbine et de la puissance extraite. Le modèle réduit de bulbe est fidèle aux machines modernes avec un diffuseur particulièrement divergent. Les performances de la turbine sont mesurées sur une large gamme de points d’opération pour déterminer les conditions les plus intéressantes pour l’étude du décollement et pour étudier la distribution paramétrique de ce phénomène. La pression est mesurée le long de l’aspirateur par des capteurs dynamiques affleurants alors que les champs de vitesse dans la zone de décollement sont mesurés avec une méthode PIV à deux composantes. Les observations à la paroi sont pour leur part faites à l’aide de brins de laine. Pour un débit suffisant, le gradient de pression adverse induit par la géométrie du diffuseur affaiblit suffisamment la couche limite, entraînant ainsi l’éjection de fluide de la paroi le long d’une large enveloppe tridimensionelle. Le décollement instationnaire tridimensionnel se situe dans la même zone du diffuseur indépendamment du point d’opération. L’augmentation du débit provoque à la fois une extension de la zone de décollement et une augmentation de l’occurrence de ses manifestations. La position et la forme du front de décollement fluctue significativement sans périodicité. L’analyse topologique et celle des tourbillons des champs de vitesse instantanés montrent une topologie du front de décollement complexe qui diffère beaucoup d’une réalisation à l’autre. Bien que l’écoulement soit turbulent, les tourbillons associés aux foyers du front sont clairement plus gros et plus intenses que ceux de la turbulence. Cela suggère que le mécanisme d’enroulement menant aux tourbillons du décollement est clairement distinct des mécanismes de la turbulence. / This thesis presents an experimental investigation of flow separation inside the diffuser of a small scale model of a bulb turbine. The flow separation occurs when the turbine is operated at high discharge and it reduces the diffuser effective area. In the case of bulb turbines, the kinetic energy recovered by the diffuser represents an important part of the total net head available for the runner energy extraction. The decrease of the diffuser efficiency leads to a sudden drop in the turbine efficiency and in the power extraction. The small scale model is faithful to modern turbines with a particularly divergent diffuser. The turbine performances are measured in a large range of operating conditions in order to select the most interesting ones and to investigate the parametric range of the phenomena. The pressure is measured along the diffuser by flush mounted dynamic sensors while the velocity fields inside the separation zone are obtained by a two-component PIV method. Separation observations on the wall are done using tufts. For a sufficient flow rate, the adverse pressure gradient induced by the diffuser geometry sufficiently weakens the boundary layer, thus leading to fluid ejection from the wall along a large three-dimensional envelope. The three-dimensional unsteady flow separation zone is located in the same area independently of the operating points. The flow rate increase leads to a wider separation zone occurring more frequently. The separation front fluctuates significantly both in location and in shape with no periodicity. Topological and vortex analyses on instantaneous velocity fields show a complex separation front topology which differs greatly from one realisation to another. Despite the highly turbulent flow, the separation front vortices are definitely bigger and more intense than turbulent vortices. It suggests that the roll-up mechanisms leading to separation surface vortices appear to be distinct from those of turbulent vortices.

Le tremblement est un phénomène qui touche diverses parties de l'avion, pouvant créer des problèmes de fatigue structurelle ou de confort et limiter l'enveloppe de vol. Il est par conséquent important de comprendre ce phénomène afin de mieux le prévoir. Quatre types de tremblement sont étudiés dans cette thèse, le tremblement de l'empennage horizontal, le tremblement intrados voilure, le tremblement de l'AFP et le tremblement de l'extrados voilure. La première partie du mémoire détaille la physique de ces phénomènes à partir d'informations recueillies dans la littérature et de l'analyse de différentes campagnes d'essais en soufflerie. Cette partie met ainsi en évidence les caractéristiques instationnaires propres à chaque type de tremblement. La deuxième partie présente l'état de l'art sur la prévision du tremblement, permettant de montrer non seulement les défauts associés aux méthodes antérieures de modélisation mais également les avantages de l'utilisation de certains outils comme l'approche numérique pour mieux comprendre ces phénomènes. Finalement, la troisième partie présente le nouveau modèle semi-empirique basé sur les fonctions de cohérence et développé pour mieux représenter les phénomènes de tremblement étudiés. Une validation de ce modèle a été effectuée sur les différentes campagnes d'essais en soufflerie, donnant des résultats probants pour la plupart des phénomènes étudiés.

These thesis review some recent results on the construction of very high order multidimensional upwind schemes for the solution of steady and unsteady conservation laws on unstructured triangular grids.<p>We also consider the extension to the approximation of solutions to conservation laws containing second order dissipative terms. To build this high order schemes we use a subtriangulation of the triangular Pk elements where we apply the distribution used for a P1 element.<p>This manuscript is divided in two parts. The first part is dedicated to the design of the high order schemes for scalar equations and focus more on the theoretical design of the schemes. The second part deals with the extension to system of equations, in particular we will compare the performances of 2nd, 3rd and 4th order schemes.<p><p>The first part is subdivided in four chapters:<p>The aim of the second chapter is to present the multidimensional upwind residual distributive schemes and to explain what was the status of their development at the beginning of this work.<p>The third chapter is dedicated to the first contribution: the design of 3rd and 4th order quasi non-oscillatory schemes.<p>The fourth chapter is composed of two parts: we start by understanding the non-uniformity of the accuracy of the 2nd order schemes for advection-diffusion problem. To solve this issue we use a Finite Element hybridisation.<p>This deep study of the 2nd order scheme is used as a basis to design a 3rd order scheme for advection-diffusion.<p>Finally, in the fifth chapter we extend the high order quasi non-oscillatory schemes to unsteady problems.<p>In the second part, we extend the schemes of the first part to systems of equations as follows:<p>The sixth chapter deals with the extension to steady systems of hyperbolic equations. In particular, we discuss how to solve some issues such as boundary conditions and the discretisation of curved geometries.<p>Then, we look at the performance of 2nd and 3rd order schemes on viscous flow.<p>Finally, we test the space-time schemes on several test cases. In particular, we will test the monotonicity of the space-time non-oscillatory schemes and we apply residual distributive schemes to acoustic problems. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

En aérodynamique automobile, les écoulements sont caractérisés par des structures fortement décollées, tridimensionnelles et instationnaires qui sont à l'origine de la trainée du véhicule et de la génération de bruit. Une configuration modèle générique de cette aérodynamique a été proposée. On étudie de manière expérimentale l'interaction du sillage tridimensionnel d'un disque de diamètre D placé à une distance H d'une paroi. Sur la base de cet écoulement, l'objectif est de développer une analyse instationnaire des structures de l'écoulement et de leur couplage avec la pression fluctuante. On a mis en oeuvre pour cela des moyens d'essais modernes tels que la mesure synchronisée pression/vitesse ou la mesure HS-PIV (High Speed - PIV) et développé des moyens d'analyse (analyse corrélatoire, outils lagrangien FTLE, moyenne conditionnelle et estimation stochastique). Pour H/D=0.75, le battement du sillage associé au lâcher tourbillonnaire est mis en évidence. Les moyennes de phase, déterminées à l'aide d'une analyse POD du signal de pression, associent le déplacement vertical du sillage au lâcher de structures cohérentes dans le plan de symétrie. L'utilisation de l'estimation stochastique de type spatio-temporelle, dont la supériorité par rapport à l'approche spatiale a été démontrée, a permis de montrer que les structures cohérentes sont reliées de manière linéaire aux fluctuations de pression. Pour H/D=0.3, l'interaction du sillage du disque avec la paroi a considérablement changé. Dans la zone proche du point de décollement, il a été montré à l'aide de l'outil lagrangien FTLE que la fréquence dominante de 30Hz sur les DSP de pression est associée aux fluctuations du point de décollement.