Questa tesi si è posta l'obiettivo di migliorare il codice commerciale CFD Ansys Fluent, per ottenere un solutore in grado di compiere simulazioni accurate di flussi in ebollizione nel regime slug flow. Un codice numerico affidabile permette una miglior comprensione della dinamica della bolla causata dall'evaporazione, rendendo possibile la stima dello scambio termico alla parete. Per limitare il costo computazionale delle simulazioni, il problema è modellato con una formulazione assialsimmetrica. Le fasi liquido e vapore sono incomprimibili ed in moto laminare. Attraverso un approccio di tipo single fluid, le equazioni che governano il moto sono scritte come per un flusso a fase singola, tuttavia discontinuità ed effetti di interfaccia vanno introdotti e discretizzati propriamente.