To simulate the phenomena in the mold region of continuous casting by coupling fluid flow and solidification, a three-dimensional mathematical model has been developed based on the K-ε turbulence equations and the SIMPLER algorithm. A pseudo source term was introduced into the energy equation to account for the latent heat and kinetic energy. The fluid flow in the mushy zone was calculated by defining the fluid viscosity as a function of the solid fraction in the mushy zone. Fine meshes in the solid region improve convergence and reduce iteration time. Comparison of the fluid flow and temperature distribution with and without solidification shows that although the solid shell in the mold is thin, it still greatly affects the flow pattern. The numerical results obtained provide details of the fluid flow and solidification phenomena which can be used to optimize the nozzle structure and other process parameters in continuous casting. To simulate the phenomena in the mold region of continuous casting by coupling fluid flow and solidification, a three-dimensional mathematical model has been developed based on the K-ε turbulence equations and the SIMPLER algorithm. A pseudo source term was introduced into the energy equation to account for the latent heat and kinetic energy. The fluid flow in the mushy zone was calculated by defining the fluid viscosity as a function of the solid fraction in the mushy zone. Fine meshes in the solid region improve convergence and reduce iteration time. of the fluid flow and temperature distribution with and without solidification shows that though the solid shell in the mold is thin, it still greatly affects the flow pattern. The numerical results obtain provide details of the fluid flow and solidification phenomena which can be used to optimize the nozzle structure and other process parameters in continuous casting.