The temperature profile of Czochralski single crystal growth system was simulated considering the fluid flow and surface radiation heat transfer. View factors of surface elements were calculated for radiation heat transfer. Two phases(solid and liquid) were treated as a continuous phase by assigning artificial large viscosity to the solid phase and latent heat was accounted by iterative heat revolution method. The solidification model was applied to solid front of the pure Ga during the melting to verify the model. The whole simulation model of CZ system was applied to the growth Al single crystal.