Un nouveau code de simulation numérique pour la propagation des ondes de choc acoustiques dans des géométries complexes a été développé. Le point de départ a été la méthode de Galerkin discontinu qui utilise des maillages non structurés (ici des éléments triangulaires), particulièrement adaptés aux géométries complexes. Cependant, cette discrétisation conduit à l'apparition d'oscillation de Gibbs. Pour pallier ce problème, nous avons choisi d'introduire de la viscosité artificielle au voisinage des chocs. Cela a nécessité le développement de trois outils originaux : (i) un nouveau détecteur de choc sensible aux ondes de chocs acoustiques sur des maillages non structurés, (ii) un nouveau terme de viscosité artificielle dans les équations de l'acoustique non linéaire défini élément par élément et (iii) un nouveau terme permettant de régler le niveau de viscosité locale à partir du raidissement des fronts d'onde. Le code de calcul a été utilisé pour étudier deux configurations différentes. La première concerne la réflexion d'ondes de choc acoustiques sur des surfaces rigides. Différents régimes de réflexion ont alors été observés allant, de la réflexion classique de Snell Descartes jusqu'à celui dit de réflexion faible de Von Neumann. La deuxième configuration était consacrée à la focalisation d'ondes de choc acoustiques produites par un transducteur à haute intensité (comme ceux utilisés en HIFU). Un soin particulier a été pris pour étudier le calcul de l'intensité et pour étudier l'interaction entre les ondes de choc et des obstacles placés dans la région du foyer. / A new numerical solver for the propagation of acoustical shock waves in complex geometry has been developed. This is done starting from the discontinuous Galerkin method. This method is based on unstructured mesh (triangular elements here), and so, naturally it is well-adapted for complex geometries. Nevertheless, the discretization induces Gibbs oscillations. To manage this problem, we choose to introduce some artificial viscosity only in the vicinity of the shocks. This necessitates the development of three original tools. First of all, a new shock sensor for unstructured mesh sensitive to acoustical shock waves has been designed. It senses where the local artificial viscosity has to be introduced thanks to a reformulation of a new element centred smooth artificial viscosity term in the equations. Finally, the amount of viscosity is computed by the introduction of an original notion of gradient factor linked to the steepening of the waveform. The numerical solver has been used to investigate two different physical situations. The first one is the nonlinear reflection of acoustical shock waves on rigid surfaces. Different regimes of reflection have been observed ranging from the linear Snell Descartes reflection to the weak von Neumann case. The second configuration deals with the focusing of shock waves produced by high intensity transducers (like in HIFU). Special attention has been given to the careful computation of intensity and to the interaction between the shock waves and obstacles in the region of the focus.

Dans le cadre du développement récent des schémas numériques compacts d’ordre élevé, tels que la méthode de Galerkin discontinu (discontinuous Galerkin) ou la méthode des différences spectrales (spectral differences), nous nous intéressons aux difficultés liées à l’utilisation de ces méthodes lors de la simulation de solutions discontinues.L’utilisation par ces schémas numériques d’une représentation polynomiale des champs les prédisposent à fournir des solutions fortement oscillantes aux abords des discontinuités. Ces oscillations pouvant aller jusqu’à l’arrêt du processus de simulation, l’utilisation d’un dispositif numérique de détection et de contrôle de ces oscillations est alors un prérequis nécessaire au bon déroulement du calcul. Les processus de limitation les plus courants tels que les algorithmes WENO ou l’utilisation d’une viscosité artificielle ont d’ores et déjà été adaptés aux différentes méthodes compactes d’ordres élevés et ont permis d’appliquer ces méthodes à la classe des écoulements compressibles. Les différences entre les stencils utilisés par ces processus de limitation et les schémas numériques compacts peuvent néanmoins être une source importante de perte de performances. Dans cette thèse nous détaillons les concepts et le cheminement permettant d’aboutir à la définition d’un processus de limitation compact adapté à la description polynomiale des champs. Suite à une étude de configurations monodimensionnels, différentes projections polynomiales sont introduites et permettent la construction d’un processus de limitation préservant l’ordre élevé. Nous présentons ensuite l’extension de cette méthodologie à la simulation d’écoulements compressibles bidimensionnels et tridimensionnels. Nous avons en effet développé les schémas de discrétisation des différences spectrales dans un code CFD non structuré, massivement parallèle et basé historiquement sur une méthodologie volumes finis. Nous présentons en particulier différents résultats obtenus lors de la simulation de l’interaction entre une onde de choc et une couche limite turbulente. / Following the recent development of high order compact schemes such as the discontinuous Galerkin or the spectraldifferences, this thesis investigates the issues encountered with the simulation of discontinuous flows. High order compactschemes use polynomial representations which tends to introduce spurious oscillations around discontinuities that can lead to computational failure. To prevent the emergence of these numerical issues, it is necessary to improve the schemewith an additional procedure that can detect and control its behaviour in the neighbourhood of the discontinuities,usually referred to as a limiting procedure or a limiter. Most usual limiters include either the WENO procedure, TVB schemes or the use of an artificial viscosity. All of these solutions have already been adapted to high order compact schemes but none of these techniques takes a real advantage of the richness offered by the polynomial structure. What’s more, the original compactness of the scheme is generally deteriorated and losses of scalability can occur. This thesis investigates the concept of a compact limiter based on the polynomial structure of the solution. A monodimensional study allows us to define some algebraic projections that can be used as a high-order tool for the limiting procedure. The extension of this methodology is then evaluated thanks to the simulation of different 2D and 3D test cases. Those results have been obtained thanks to the development of a parallel solver which have been based on a existing unstructured finite volume CFD code. The different exposed studies detailed end up to the numerical simulation of the shock turbulent boundary layer.