A modelling of the induced melt pool movements during iron deep penetration welding by CW Nd–YAG laser beam is developed, with the aim of studying the effect of inert gas jet on the liquid metal flow. The keyhole shape is calculated as a function of the operating parameters and the derived boundary conditions on the keyhole walls are used, by applying a solidification model, to calculate the frontiers of the assumed stationary molten pool. The vapour flow generated inside the keyhole, the surface tension and the recoil pressure are considered as the mechanisms producing the melt pool movements, and the evolution of the interfaces between the keyhole, the molten pool and the surrounding air is simulated by the multiphase volume of fluid method. We demonstrate, in accordance with the experimental observations, that an important stirring effect is generated inside the melt pool, and moreover by extending the modelling to the study of the effect of 45° inert gas projection at 100 m s−1 on the metallic liquid pool, it is shown that a more uniform flow is obtained.