Fluidised beds are used in a variety of processes because of their favourable mass and heat transfer characteristics. In this and many other processes, non-spherical particles are commonplace, which can drastically affect the fluidisation behaviour. In this study, we use numerical models to study non-spherical fluidisation behaviour in detail. A crucial step in the development of the numerical model is a detailed validation with experimental data. The validated model can then be used with confidence for further investigations. In this study, the results obtained from CFD-DEM modelling are compared with detailed experiments (Mahajan et al., 2017). The particles used are of spherocylindrical shape with an aspect ratio 4. We discuss the numerical modelling strategy including the DEM contact detection algorithm and accurate voidage calculation algorithm. The non-spherical single particle drag model of Holzer and Sommerfeld (2008) is compared with a DNS drag model for spherocylindrical particles developed in-house. We propose two new voidage correction models and compare results with the (Di Felice, 1994) model. The pressure drop, bed height, particle orientation, particle circulation, stacking of particles and coordination number obtained from simulations are compared with experiments. The numerical measurements show good agreement with experiments. Similar to experiments, simulations show that rod-like particles are prone to interlocking behaviour. At high gas flow rates above the minimum fluidisation velocity, vigorously bubbling fluidisation is observed, with gas bubbles moving up through the center and particles moving down at the side walls. The orientation of particles in the fluidised state do not match with the experiments when hydrodynamic torque is neglected. The importance of hydrodynamic torque and multi-particle drag in CFD-DEM modelling of non-spherical particles is demonstrated through these results. (C) 2018 The Author(s). Published by Elsevier Ltd.