Turbulent flow and heat transfer coupled with solidification in slab continuous casting mold was studied by numerical simulation method. Volume of fluid (VOF) model is used to solve steel-air two-phase flow problem and enthalpy-porosity scheme is introduced to solve the fluid flow problem involving solidification. Contributions of various nozzle port angles and port widths and heights on the free surface fluctuation and the thickness of solidifying shell in slab mold were particularly investigated, based on which the structure of submerged entry nozzle was optimized. Flow inside the common nozzle port cannot fill the entire ouaet area, having a recireulation in the upper portion of the port, which is enlarged for the nozzle port with both larger height and width. Results show that the flow in mold cavity is mainly controlled by the nozzle port angle. The increase of the angle of upper face of the port to shape a roughly streamlined inner-wall improves the effective area fraction of the nozzle, resulting in less jet impingement, weaker free surface turbulence and thicker solidifying steel shell.