Cette thèse fait référence au bruit de choc généré par des jets sousdétendus simple ou en configuration co-axial. Le bruit de choc est généré par l'interaction entre les structures turbulentes de la couche de cisaillement et le réseaux de cellules de choc développé dans le cône potentiel du jet. Afin d'étudier le bruit choc, simulations à grandes échelles adaptés pour l'aéro-acoustique sont effectué avec des schémas d'ordre élevé qui permet une approche nondissipative et non-dispersive. Les résultats sont analysés et comparés avec des résultats expérimentaux. Notamment, une filtrage hydrodynamique et acoustique est réalisé dans le champ proche pour analyser les modes azimutaux acoustiques et hydrodynamiques. En outre, un analyse basé sur la transformé en ondelettes est mis en oeuvre pour identifier les caractéristiques acoustiques et hydrodynamiques importants des jets supersoniques. / This thesis deals with the shock-cell noise generated by under-expanded supersonic jets in single- and dualstream configurations. Shock-cell noise is generated by the interaction between the turbulent structures of the shear-layer and the shock-cell system developed in the potential core of the jet. In order to study shock-cell noise, large eddy simulations adapted to aeroacoustics are carried out using high-order compact schemes that allow for a non-dissipative nondispersive approach. The results are analyzed and compared to experimental results. In particular, an acoustic-hydrodynamic filtering is carried out in the near field in order to analyze the acoustic and hydrodynamic azimuthal modes. Moreover, a wavelet-based analysis is implemented in order to identify the relevant acoustic and hydrodynamic features of the supersonic jets.

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

La simulation numérique est aujourd'hui un outils majeur dans la conception des objets aérodynamiques, que ce soit dans l'aéronautique, l'automobile, l'industrie navale, etc... Un des défis majeurs pour repousser les limites des codes de simulation est d'améliorer leur précision, tout en utilisant une quantité fixe de ressources (puissance et/ou temps de calcul). Cet objectif peut être atteint par deux approches différentes, soit en construisant une discrétisation fournissant sur un maillage donné une solution d'ordre très élevé, soit en construisant un schéma compact et massivement parallèlisable, de manière à minimiser le temps de calcul en distribuant le problème sur un grand nombre de processeurs. Dans cette thèse, nous tentons de rassembler ces deux approches par le développement et l'implémentation de Schéma Distribuant le Résidu (RDS) d'ordre très élevé et de compacité maximale. Ce manuscrit commence par un rappel des principaux résultats mathématiques concernant les Lois de Conservation hyperboliques (CLs). Le but de cette première partie est de mettre en évidence les propriétés des solutions analytiques que nous cherchons à approcher, de manière à injecter ces propriétés dans celles de la solution discrète recherchée. Nous décrivons ensuite les trois étapes principales de la construction d'un schéma RD d'ordre très élevé : - la représentation polynomiale d'ordre très élevé de la solution sur des polygones et des polyèdres; - la description de méthodes distribuant le résidu de faible ordre, compactes et conservatives, consistantes avec une représentation polynomiale des données de très haut degré. Parmi elles, une attention particulière est donnée à la plus simple, issue d'une généralisation du schéma de Lax-Friedrichs (\LxF); - la mise en place d'une procédure préservant la positivité qui transforme tout schéma stable et linéaire, en un schéma non linéaire d'ordre très élevé, capturant les chocs de manière non oscillante. Dans le manuscrit, nous montrons que les schémas obtenus par cette procédure sont consistants avec la CL considérée, qu'ils sont stables en norme $\L^{\infty}$ et qu'ils ont la bonne erreur de troncature. Même si tous ces développements théoriques ne sont démontrés que dans le cas de CLs scalaires, des remarques au sujet des problèmes vectoriels sont faites dès que cela est possible. Malheureusement, lorsqu'on considère le schéma \LxF, le problème algébrique non linéaire associé à la recherche de la solution stationnaire est en général mal posé. En particulier, on observe l'apparition de modes parasites de haute fréquence dans les régions de faible gradient. Ceux-ci sont éliminés grâce à un terme supplémentaire de stabilisation dont les effets et l'évaluation numérique sont précisément détaillés. Enfin, nous nous intéressons à une discrétisation correcte des conditions limites pour le schéma d'ordre élevé proposé. Cette théorie est ensuite illustrée sur des cas test scalaires bidimensionnels simples. Afin de montrer la généralité de notre approche, des maillages composés uniquement de triangles et des maillages hybrides, composés de triangles et de quandrangles, sont utilisés. Les résultats obtenus par ces tests confirment ce qui est attendu par la théorie et mettent en avant certains avantages des maillages hybrides. Nous considérons ensuite des solutions bidimensionnelles des équations d'Euler de la dynamique des gaz. Les résultats sont assez bons, mais on perd les pentes de convergence attendues dès que des conditions limite de paroi sont utilisées. Ce problème nécessite encore d'être étudié. Nous présentons alors l'implémentation parallèle du schéma. Celle-ci est analysée et illustrée à travers des cas test tridimensionnel de grande taille. / Numerical simulations are nowadays a major tool in aerodynamic design in aeronautic, automotive, naval industry etc... One of the main challenges to push further the limits of the simulation codes is to increase their accuracy within a fixed set of resources (computational power and/or time). Two possible approaches to deal with this issue are either to contruct discretizations yielding, on a given mesh, very high order accurate solutions, or to construct compact, massively parallelizable schemes to minimize the computational time by means of a high performance parallel implementation. In this thesis, we try to combine both approaches by investigating the contruction and implementation of very high order Residual Distribution Schemes (RDS) with the most possible compact stencil. The manuscript starts with a review of the mathematical theory of hyperbolic Conservation Laws (CLs). The aim of this initial part is to highlight the properties of the analytical solutions we are trying to approximate, in order to be able to link these properties with the ones of the sought discrete solutions. Next, we describe the three main steps toward the construction of a very high order RDS: - The definition of higher order polynomial representations of the solution over polygons and polyhedra; - The design of low order compact conservative RD schemes consistent with a given (high degree) polynomial representation. Among these, particular accest is put on the simplest, given by a generalization of the Lax-Friedrich's (\LxF) scheme; - The design of a positivity preserving nonlinear transformation, mapping first-order linear schemes onto nonlinear very high order schemes. In the manuscript, we show formally that the schemes obtained following this procedure are consistent with the initial CL, that they are stable in $L^{\infty}$ norm, and that they have the proper truncation error. Even though all the theoretical developments are carried out for scalar CLs, remarks on the extension to systems are given whenever possible. Unortunately, when employing the first order \LxF scheme as a basis for the construction of the nonlinear discretization, the final nonlinear algebraic equation is not well-posed in general. In particular, for smoothly varying solutions one observes the appearance of high frequency spurious modes. In order to kill these modes, a streamline dissipation term is added to the scheme. The analytical implications of this modifications, as well as its practical computation, are thouroughly studied. Lastly, we focus on a correct discretization of the boundary conditions for the very high order RDS proposed. The theory is then extensively verified on a variety of scalar two dimensional test cases. Both triangular, and hybrid triangular-quadrilateral meshes are used to show the generality of the approach. The results obtained in these tests confirm all the theoretical expectations in terms of accuracy and stability and underline some advantages of the hybrid grids. Next, we consider solutions of the two dimensional Euler equations of gas dynamics. The results obtained are quite satisfactory and yet, we are not able to obtain the desired convergence rates on problems involving solid wall boundaries. Further investigation of this problem is under way. We then discuss the parallel implementation of the schemes, and analyze and illustrate the performance of this implementation on large three dimensional problems. Due to the preliminary character and the complexity of these three dimensional problems, a rather qualitative discussion is made for these tests cases: the overall behavior seems to be the correct one, but more work is necessary to assess the properties of the schemes in three dimensions.

In this work we review 12 years of developments in the field of residual based discretizations for hyperbolic problems and their application to the solution of the shallow water equations. Fundamental concepts related to the topic are recalled and he construction of second and higher order schemes for steady problems is presented. The generalization to time dependent problems by means of multi-step implicit time integration, space-time, and genuinely explicit techniques is thoroughly discussed. Finally, the issues of C-property, super consistency, and wetting/drying are analyzed in this framework showing the power of the residual based approach.

Ces travaux de recherche ont pour but d’évaluer la méthode dite de la "Simulation auxEchelles Adaptées" (SAS pour Scale-Adaptive Simulation). Cette approche coïncide avec uneapproche RANS classique dans les zones pariétales attachées et adapte le niveau de viscositéturbulente dans les zones décollées pour y permettre une résolution partielle des structures turbulentes.Dans une première partie, une analyse théorique du modèle SAS original a été menéeet a permis de développer une correction visant à favoriser l’adaptation du niveau de viscositéturbulente dans les zones sièges d’instabilités de type Kelvin-Helmholtz. Le modèle ainsi corrigéest nommé SAS-αL. Les modèles SAS et SAS-αL ont été implantés dans le code de calculNavier-Stokes elsA de l’ONERA. À l’issue de cette étape, trois cas académiques d’écoulementsturbulents instationnaires, cylindre à grand nombre de Reynolds, marche descendante et cavitétranssonique, ont été simulés grâce aux trois modèles de turbulence SST, SAS et SAS-αL. Outreune comparaison aux bases de données expérimentales disponibles, une attention particulièrea été portée à l’influence de paramètres numériques tels que des schémas numériques d’ordreélevé. Enfin, afin d’étudier la viabilité de l’approche SAS dans un contexte industriel, les troismodèles de turbulence ont été testés sur une configuration issue de l’industrie aéronautique etcorrespondant à la sortie d’air chaud d’un système de dégivrage des nacelles d’avion. La comparaisondes prévisions obtenues avec les modèles SST, SAS et SAS-αL aux données expérimentalesobtenues à l’ONERA a permis de montrer un gain de précision grâce à l’emploi de l’approcheSAS et ce pour un coût de calcul compatible avec un cycle de conception industrielle. / This research work is meant to assess an upgraded URANS approach, namely the Scale-Adaptive Simulation (SAS). This method is similar to a conventional RANS approach (namelythe SSTmodel) in attached areas and is able to adapt the eddy-viscosity level in detached areas toensure the resolution, at least partially, of the turbulent structures. In a first part of this researchwork, an improvement of the SAS approach is suggestedto allowa better sensitivity of themodelto instabilities such as Kelvin-Helmholtz ones. This "improved" model is referred to as SAS-αLmodel. Both SAS and SAS-αL models were implemented in the ONERA Navier-Stokes solverelsA and both of themaswell as the SSTmodelwere tested on academic test cases : a cylinder in acrossflowat a high Reynolds number, a backward-facing step flowcorresponding to theDriver&Seegmiller experiment and the transonic flow over the M219 cavity experimentally investigatedby de Henshaw. The influence of the numerical parameters was deeply investigated and particularattention was paid to the high-order space-discretization schemes effects. The reliabilityof the SAS approach in an industrial framework was assessed on an aeronautic configurationnamely a nacelle de-icing device. Comparisons between the threemodels (SST, SAS and SAS-αL)and an experimental database available at ONERA - The French Aerospace Lab have shown thebetter accuracy of the SAS approach as well as the high potential of the SAS-αL model.