In continuous casting process, the tundish optimization fluid flow pattern has a considerable role to removal inclusions and the quality of steel product. To reach a suitable fluid flow pattern in tundish to separate non-metallic inclusions from melt and traort to slag by upward molten flow is the aim of this work. In this idea, there is no need to use dam, weir and turbo stopper, therefore dead zones are reduced. The non-metallic inclusion particles have lower density than molten steel in real unit and the direction of flow pattern is upward, thus, it will rise up to the slag-steel interface. The fundamentals of the process, physical and numerical modellings with persistence on numerical modellings have been conducted in this work. Because of symmetric in distribution of fluid flow in model, one half of the water model is investigated by numerical method and physical model is carried out with the utilization of a 1:4 scale model with similarity criteria as much as possible. The geometry of water model is generated using the commercial package environment GAMBIT, and then ready to inter the computational fluid dynamics (CFD) field is used to analysis, design optimization of fluid flow and other applications of heat traort. With attention to necessity of fluid flow pattern proposed, two topics, the wave shape of the free surface and the trajectory of inclusions in tundish were investigated using mathematical modellings and validity of numerical modelling is carried out in line with the physical modelling to examine details of the flow pattern in water model. The situation of free surface is determined by solving the governing differential equations for turbulent flows using finite difference method at isothermal and unsteady state conditions with k-? model. The initial and boundary conditions are determined by VoF (Volume of Fluid) method. In VoF methods, however, one usually introduces a liquid indicator, or "Volume-of-Fluid" field for liquid and gas which is a property of the fluid is considered and instantaneous properties viscosity and density are related linearly and with the commercial CFD package environment, FLUENT is solved and results in form contours of density (mixture) versus time achieved and in physical model with a cane shroud is followed. During the experiment a flat free surface is obviously seen. The motion pass of inclusions are considered using forces such as, drag, buoyancy and gravitational forces. Assumptions are spherical shape inclusions and monotonousness of inclusions size. Trajectories of inclusions to slag-melt interface with Stokes motion are determined by using "Discrete Phase". The contours of residence time and magnitude of velocity versus time were achieved and show how inclusions guide to free surface.