La thèse porte sur la simulation, en régime fréquentiel, du rayonnement acoustique en écoulement subsonique quelconque et dans un domaine infini. L'approche choisie s'appuie sur la résolution d'un système équivalent aux équations d'Euler linéarisées : le modèle de Galbrun. Ce modèle repose sur une représentation mixte Lagrange-Euler et aboutit à une équation dont l'unique inconnue est la perturbation du déplacement Lagrangien. Une des difficultés de l'approche de Galbrun est qu'une discrétisation directe de cette équation par une méthode d'éléments finis standard n'est pas stable. Un moyen de contourner cet obstacle est d'écrire une équation augmentée en ajoutant une nouvelle inconnue, le rotationnel du déplacement, appelée par abus vorticité. Cette approche conduit à un système qui couple une équation de type équation des ondes avec une équation de transport en régime fréquentiel. Et elle permet l'utilisation de couches parfaitement adaptées (PML) pour borner le domaine de calcul. La première partie du manuscrit est dédiée à l’étude de l’équation de transport harmonique et de sa résolution numérique, en particulier par un schéma de type Galerkin discontinu. Un des points délicats est lié au caractère oscillant des solutions de l'équation. Une fois cette étape franchie, la résolution du problème de propagation acoustique a été abordée. Une approximation basée sur l'utilisation d'éléments finis mixtes continus-discontinus avec couches parfaitement adaptées (PML) a été étudiée. En particulier, les caractères bien posés des problèmes continu et discret ainsi que la convergence du schéma numérique ont été démontrés sous certaines conditions sur l'écoulement porteur. Enfin, une mise en œuvre a été effectuée. Les résultats montrent la validité de cette approche mais aussi sa pertinence dans le cas d'écoulements complexes, voire d'écoulements dits instables / This thesis deals with the numerical simulation of time harmonic acoustic propagation in an arbitrary mean flow in an unbounded domain. Our approach is based on an equation equivalent to the linearized Euler equations called the Galbrun equation. It is derived from a mixed Eulerian-Lagrangian formulation and results in a single equation whose only unknown is the perturbation of the Lagrangian displacement. A direct solution using finite elements is unstable but this difficulty can be overcome by using an augmented equation which is constructed by adding a new unknown, the vorticity, defined as the curl of the displacement. This leads to a set of equations coupling a wave like equation with a time harmonic transport equation which allows the use of perfectly matched layers (PML) at artificial boundaries to bound the computational domain. The first part of the thesis is a study of the time harmonic transport equation and its approximation by means of a discontinuous Galerkin scheme, the difficulties coming from the oscillating behaviour of its solutions. Once these difficulties have been overcome, it is possible to deal with the resolution of the acoustic propagation problem. The approximation method is based on a mixed continuous-Galerkin and discontinuous-Galerkin finite element scheme. The well-posedness of both the continuous and discrete problems is established and the convergence of the approximation under some mean flow conditions is proved. Finally a numerical implementation is achieved and numerical results are given which confirm the validity of the method and also show that it is relevant in complex cases, even for unstable flows