Cette étude propose une approche aéroacoustique hybride pour le calcul du bruit rayonné par des écoulements subsoniques turbulents anisothermes. La partie aérodynamique est obtenue à l'aide d'une simulation numérique directe des équations de Navier-Stokes dans une approximation à faible nombre de Mach. Cette formulation permet de se libérer des effets de la compressibilité - numériquement pénalisants - tout en préservant l'influence des phénomènes relatifs à la dilatation thermique sur l'écoulement. La propagation acoustique est calculée par la résolution des équations d'Euler linéarisées. La définition rigoureuse des sources acoustiques constitue le lien entre ces deux étapes. Des sources spécifiques aux écoulements anisothermes, compatibles avec celles déjà connues pour les écoulements isothermes, sont obtenues. Cette approche est d'abord validée pour des couches de mélange isothermes et anisothermes par comparaison à des calculs directs. D'autres validations sont réalisées pour une couche de mélange isotherme spatiale et confrontés à des résultats de la littérature. Les contributions au bruit des différents termes sources sont examinées pour une couche de mélange spatiale anisotherme. L'évolution du rayonnement acoustique en fonction du rapport de températures est étudié pour une couche de mélange temporelle. Les résolutions numériques reposent sur une discrétisation spatiale par des schémas aux différences finies d'ordres élevés (schémas compacts) et des schémas d'intégration en temps de Runge-Kutta. Ces schémas d'ordres élevés, et les conditions aux limites performantes respectent les exigences numériques spécifiques à l'aéroacoustique.