La connaissance du comportement thermique des fours et des pièces est un problème difficile et cependant essentiel dans les thématiques de recherche industrielles actuelles. L'industrie cherche à se doter de moyens numériques de plus en plus efficaces tout en réduisant sans cesse le temps de calcul afin de modéliser des pièces et des assemblages de plus en plus réalistes. Les cinétiques de chauffe ainsi que la distribution de la température dans l'enceinte et dans les pièces traitées doivent également être maîtrisées, afin d'améliorer la qualité des produits chauffées en terme de résistance et dureté. Le travail présenté dans cette thèse porte donc sur le développement de méthodes de résolution numérique pour la simulation du transfert thermique. Ces méthodes permettent le calcul couplé de la température des pièces avec l'environnement du four afin d'optimiser la géométrie de ces pièces, leur position dans le four mais également les conditions initiales de chauffe. Le calcul de dynamique des fluides (CFD) peut être considéré comme un puissant outil technique de prédiction de l'écoulement des fluides et du transfert thermique dans des cas industriels réalistes. La première partie de la thèse porte sur l'étude de schémas numériques avancés relatifs aux méthodes éléments finis stabilisés. Ces méthodes de calculs efficaces ont étaient utilisées pour simuler des écoulements instationnaires et des transferts thermiques conjugués. Par la suite, deux modèles de turbulence (modèle k-epsilon et modèle Large Eddy Simulation) sont introduits et utilisés pour prendre en compte les écoulements complexes et turbulents dans l'enceinte du four. Le transfert radiatif est assuré par la résolution du modèle P1 tout en calculant un terme source volumique qui sera intégré dans l'équation de la chaleur. La deuxième partie de cette thèse se consacre à la mise en place d'une méthodologie de discrétisation robuste qui permet aux utilisateurs de générer de façon entièrement automatique un seul maillage. Ce maillage contient à la fois des domaines axés sur la résolution d'un problème fluide (air, eau, ...) mais également solides spécifiques aux structures, et ce quel que soit le niveau de détail et donc de complexité du cas étudié. Cette approche, connue sous le nom d' « immersion de volume », garantit un maillage anisotrope précis aux interfaces fluide-solide afin de capturer plus précisément les gradients thermiques et la forte discontinuité des propriétés physiques. Cette méthode offre donc une grande flexibilité dans la mise en données du problème mixte fluide-structure et aussi dans la prise en compte de plusieurs géométries (four, pièces, supports) et elle est également bien adaptée aux solveurs thermomécaniques développés.