Direct numerical simulations are performed to study the effect of particle orientation on flows through fixed random arrays of prolate ellipsoids at low Reynolds numbers.The Hermans orientation factor and Beta distribution are introduced to quantify the mean orientation and orientation deviation of the particles. The simulation resultsshow that the effect of particle orientation is profound especially when the solid volume fraction and the aspect ratio are large. With the increase of Hermans orientation factors, the drag force decreases when the flow follows a reference direction defined by the average direction of all particles' semi-major axes, while increases when the flow is perpendicular to the reference direction. Comparisons show that the traditional drag force correlations for ellipsoidal particles significantly under-predict the drag force. Based on current simulation results, new drag relations are proposed for prolate ellipsoidal particles at arbitrary aspect ratios, Hermans orientation factors and solid volume fractions. K E Y W O R D S ellipsoids, fluidization, Hermans orientation factors, multiphase flow, the drag force