Adequacy of approximation of ellipsoidal particles composed of a set of sub-spheres for numerical Discrete Element Method (DEM) simulations is examined. The algorithm of adaptive hierarchical multi-sphere (MS) model is suggested for composing elliptical particles. Numerical simulation of the piling problem is used as a test problem for evaluating the adequacy of MS model approximation in comparison to the model of smooth ellipses for multiparticle system. The accuracy of MS approximation with the increasing number of sub-spheres is examined in detail by comparison of macroscopic and microscopic parameters of granular dynamics. It was determined that the data on macroscopic parameters yielded by the MS model tend to converge to those of the smooth ellipsoid with the increasing number of the constituent sub-spheres, and the MS model approximates the smooth perfect ellipsoid with a reasonable number of sub-spheres within the limits of the appropriate tolerance. It can be concluded that a multi-sphere model remains a realistic and relatively simple particle model applicable to DEM simulations of the behaviour of the real smooth and rough elliptically shaped particles
