This paper proposes the use of a generalized finite difference method for the numerical simulation of free surface single phase flows during mould filling process which are common in some industrial processes particularly in the area of metal casting. A novel and efficient idea for the computation of the normal vectors for free surface flows is introduced and presented for the first time. The incompressible Navier–Stokes equations are numerically solved by the well-known Chorin's projection method. After we showed the main ideas behind the meshless approach, some numerical results in two and three dimensions are presented corresponding to mould filling process simulation.