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 liée au 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 cette équation de transport 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 oeuvre 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.