The drag of a non-spherical particle was reviewed and investigated for a variety of shapes (regular and irregular) and particle Reynolds numbers (Rep). Point-force models for the trajectory-averaged drag were discussed for both the Stokes regime (Re-p << 1) and Newton regime (Re-p >> I and sub-critical with approximately constant drag coefficient) for a particular particle shape. While exact solutions were often available for the Stokes regime, the Newton regime depended on: aspect ratio for spheroidal particles, surface area ratio for other regularly-shaped particles, and min-med-max area for irregularly shaped particles. The combination of the Stokes and Newton regimes were well integrated using a general method by Ganser (developed for isometric shapes and disks). In particular, a modified Clift-Gauvin expression was developed for particles with approximately cylindrical cross-sections relative to the flow, e.g. rods, prolate spheroids, And oblate spheroids with near-unity aspect ratios. However, particles with non-circular cross-sections exhibited a weaker dependence on Reynolds number, which is attributed to the more rapid transition to flow separation and turbulent boundary layer conditions. Their drag coefficient behavior was better represented by a modified Dallavalle drag model, by again integrating the Stokes and Newton regimes. This paper first discusses spherical particle drag and classification of particle shapes, followed by the main body which discusses drag in Stokes and Newton regimes and then combines these results for the intermediate regimes. (c) 2007 Elsevier B. V. All rights reserved.