L'objectif général de la thèse est d'améliorer la modélisation numérique RANS des flux thermiques turbulents notamment en proposant un modèle fonctionnant dans les trois régimes de convection thermique (forcée, mixte et naturelle).Pour ce faire, un état des lieux, non exhaustif, des modèles des flux thermiques utilisant les approches algébriques et à équations de transport, est effectué. Puis, le modèle EB-RSM (Elliptic Blending-Reynolds Stress Model) étant utilisé pour modéliser la turbulence, le principe de la pondération elliptique est appliqué aux modèles des flux thermiques turbulents algébriques EB-GGDH (EB-General Gradient Diffusion Hypothesis), EB-AFM (EB-Algebraic Flux Model) et à équations de transport EB-DFM (EB-Differential Flux Model). Une attention particulière a été apportée aux échelles de temps et de longueur utilisées pour ces modèles. Il en résulte qu'utiliser une échelle de longueur thermique différente de l'échelle de longueur dynamique et une échelle de temps mixte dans le terme de flottabilité de l'équation de la dissipation turbulente est préférable.Pour valider les formulations retenues, nous avons effectué des tests pour des fluides usuels (nombre de Prandtl de l’ordre de 1) dans les trois régimes de convection à l'aide de l'outil de calcul Code_Saturne sur des configurations académiques, semi-académiques et industrielles.Des résultats satisfaisants ont été obtenus en associant l'EB-RSM et le GGDH en convection forcée ou mixte et l'EB-RSM aux modèles EB-DFM et AFM en convection naturelle. / The PhD main objective is to improve the turbulent heat flux RANS modelling especially by proposing a model working in the three thermal convection regime (forced, mixed and natural).In order to achieve this, a non-exhaustive state of art of heat flux model, using algebraic approach and transport equations, is done. Then, as EB-RSM model (Elliptic Blending-Reynolds Stress Model) is used to model turbulence, elliptic blending approach is apply to algebraic turbulent heat flux model EB-GGDH (EB-General Gradient Diffusion Hypothesis), EB-AFM (EB-Algebraic Flux Model) and transport equation model EB-DFM (EB-Differential Flux Model). Special attention was paid to time and length scales used with these models. It follows that using a thermal length scale different from dynamic length scale and a mixed time scale in the buoyant term of turbulent dissipation equation is better.To validate these models, some test were done for common fluids (Prandtl number in the order of 1) in the three convection regimes with the tool Code_Saturne on academic, mid-academic and industrial cases.Good results were obtained combining EB-RSM with GGDH in forced or mixed convection and EB-RSM with EB-DFM or AFM in natural convection.