Accurate direct numerical simulations are performed to determine the drag, lift and torque coefficients of non-spherical particles. The numerical simulations are performed using the lattice Boltzmann method with multi-relaxation time. The motivation for this work is the need for accurate drag, lift and torque correlations for high Re regimes, which are encountered in Euler-Lagrangian simulations of fluidization and pneumatic conveying of larger non-spherical particles. The simulations are performed in the Reynolds number range 0.1 <= Re <= 2000 for different incident angles phi. Different tests are performed to analyse the influence of grid resolution and confinement effects for different Re. The measured drag, lift and torque coefficients are utilized to derive accurate correlations for specific non-spherical particle shapes, which can be used in unresolved simulations. The functional forms for the correlations are chosen to agree with the expected physics at Stokes flow as well as the observed leveling off of the drag coefficient at high Re flows. Therefore the fits can be extended to regimes outside the Re regimes simulated. We observe sine-squared scaling of the drag coefficient for the particles tested even at Re = 2000 with C-D,C-phi = C-D,C-phi=0 degrees (C-D,C-phi=90 degrees - C-D,C-phi=(0 degrees)) sin(2) phi. Furthermore, we also observe that the lift coefficient approximately scales as C-L,C-phi = (C-D,C-phi=90 degrees - C-D,C-phi=0 degrees) sin phi cos phi for the elongated particles. The current work would greatly improve the accuracy of Euler-Lagrangian simulations of larger non-spherical particles considering the existing literature is mainly limited to steady flow regimes and lower Re. (C) 2018 Elsevier Ltd. All rights reserved.