Modelling of wet grinding in tumbling mills is an interesting challenge. A key factor is that the pulp fluid and its simultaneous interactions with both the charge and the mill structure have to be handled in a computationally efficient way. In this work, the pulp fluid is modelled with a Lagrange based method based on the particle finite element method (PFEM) that gives the opportunity to model free surface flow. This method gives robustness and stability to the fluid model and is efficient as it gives possibility to use larger time steps. The PFEM solver can be coupled to other solvers as in this case both the finite element method (FEM) solver for the mill structure and the DEM solver for the ball charge. The combined PFEM-DEM-FEM model presented here can predict charge motion and responses from the mill structure, as well as the pulp liquid flow and pressure. All cases presented here are numerically modelled and validated against experimentally measured driving torque signatures from an instrumented small-scale batch ball mill equipped with a torque meter and charge movements captured from high-speed video. Numerical results are in good agreement with experimental torque measurements and the PFEM solver also improves on efficiency and robustness for solving charge movements in wet tumbling mill systems.