Ce travail porte sur la simulation numérique tridimensionnelle par éléments finis du refroidissement primaire lors de la coulée continue d'aciers via le logiciel THERCAST. L'objectif de cette étude est de modéliser avec précision les régions solides, liquides (et pâteuses) présentes dans une lingotière de coulée continue. Dans ce contexte, la stratégie globale instationnaire est utilisée : le domaine de calcul croît à la vitesse de coulée. Afin de s'affranchir des difficultés liées aux approches purement eulériennes ou langrangiennes, nous proposons une méthode de type Eulérien-Lagrangien-ALE (Arbitrairement Lagrangienne Eulérienne). La peau solidifiée est modélisée par une loi de comportement élasto-visco-plastique. Une méthode de remaillage adaptatif automatique permet de la représenter finement grâce à un maillage anisotrope très fin dans le sens des gradients de température et de contrainte. Afin de modéliser les écoulements liquides, nous avons implémenter une méthode de transport sur la densité de bulles d'argon, couplée aux équations de Navier-Stokes par le terme de force gravitationnelle. Les effets turbulents sont pris en compte par le modèle k − " standard. Tous ces développements ont été confrontés avec succès à des résultats académiques ou à des benchmarks, bien que la simulation d'une configuration industrielle prenant en compte l'ensemble de ces méthodes et phénoménologies ne soit pas encore à notre portée.