L'objectif de cette thèse est de développer un code de simulation complet pour le remplissage en fonderie de pièces de grandes dimensions (jusqu'à plusieurs mètres). Ce type de procédé fait entrer en jeu de nombreux phénomènes physiques couplés, nécessitant des méthodes numériques adaptées. La faible viscosité du métal liquide (de l'ordre de 10−6 m2/s) requiert l'emploi d'un modèle de turbulence basé sur un solveur Navier Stokes stabilisé et une méthode de suivi/capture d'interface. Nous avons pour cela choisi un approche stabilisée de type Variational Multi Scales (VMS), qui s'est révélée efficace pour simuler des nombres de Reynolds modérés, alliée à une méthode level-set permettant de déterminer de manière précise et à tout moment la position de l'interface liquide/air. La turbulence est quant à elle prise en compte grâce à un modèle dynamique de type Large Eddy Simulations (L.E.S.), ne faisant pas apparaître de paramètre empirique. Chacune de ces méthodes numériques a été confrontée à des résultats expérimentaux, numériques ou analytiques. Nous avons également conçu notre propre maquette expérimentale de remplissage d'eau, afin de valider le couplage des solveurs pour un cas représentatif. Une autre caractéristique de ces procédés à durée relativement longue (jusqu'à plusieurs dizaines de minutes) est l'importance des transferts thermiques, pouvant mener à la solidification du métal en cours de remplissage. Il convient donc de développer une méthode de résolution stabilisée de la thermique avec convection dominante. Cette méthode doit prendre en compte les variables turbulentes introduites précédemment. Enfin, nous proposons une méthode innovante pour simuler le changement de phase, basée sur une approche germination/ croissance avec fonction level-set. L'application de toutes ces méthodes au cas du remplissage avec glaçon mobile a enfin permis de valider la robustesse numérique de notre code et le bon couplage de ses différentes entités.