A new 3D finite element model, developed in a level set approach, is proposed to model hybrid gas metal arc / laser welding in multipass conditions. This heat and mass transfer model couples the resolution of the heat conservation equation, the momentum and mass conservation equations and the weld bead development. The arc welding total power is divided in two parts: one corresponding to the energy transferred to the fusion zone by droplets of melted filler material, the rest being transferred to the workpiece by the surrounding plasma. The droplets energy input is modelled as a volumetric heat source. Regarding the heat surface fluxes associated with plasma and laser heating, the "Continuum Surface Force" approach is used to model them as volumetric heat sources concentrated in the immediate neighbourhood of the metal-gas interface. A resolution scheme, consisting in ignoring high velocity fluid flow in the fusion zone, through the use of an augmented liquid viscosity, is proposed and discussed. Accordingly, the liquid thermal conductivity is enhanced to result in a realistic heat transfer to the workpiece. The formation of the weld bead is obtained through the introduction of a source term in the mass conservation equation, and the application of the normal component of surface tension forces, proportional to the mean curvature of the metal-gas interface. This approach is proposed to reduce computation time. The resolution scheme is applied to the simulation of hybrid welding of 18MND5 (ASME SA 533) steel grade. Results are compared to experimental measurements and observations operated in conditions close to industrial ones. The influence of the enhancement applied to the liquid conductivity coefficient is shown and discussed. A strong sensitivity evolution is demonstrated when varying it from the physical value to the value proposed in the welding literature. As proposed, the simplified resolution scheme leads to a good estimation of the weld bead surface development. Nevertheless, there are still noticeable differences with the whole set of experimental results that are discussed and can be explained by model limitations and insufficient knowledge of material and interfacial properties.