In this study a comparison is made between results from three Eulerian-based computational methods that predict the ice crystal trajectories and impingement on a NACA-0012 airfoil. The computational methods are being developed within CIRA (Imp2D/3D), ONERA (CEDRE/Spiree) and University of Twente (MooseMBIce). Eulerian models describing ice crystal transport are complex because physical phenomena, like drag force, heat transfer and phase change, depend on the particle's sphericity. Few correlations exist for the drag of non-spherical particles and heat transfer of these particles. The effect or non-spherical particles on the collection efficiency will be shown on a 2D airfoil.