The objective of this communication is to present the new capability at AMIL in ice accretion simulation on 2D Airfoils at low speed. AMIL, in a joint project with CIRA (Italian Aerospace Research Center), has developed a numerical model called CIRAMIL. This model is able to predict ice shapes in wet and dry regimes. The thermodynamic model used is similar to existing ones. The major difference is in the approach of calculating the surface roughness and the residual, runback and shedding liquid water masses on an airfoil surface. The numerical ice shapes are compared to rime and glaze shapes obtained experimentally in wind tunnel for a NACA0012 wing profile. The new roughness computation method generates the complex ice shapes observed experimentally in wet and dry regimes and the results agree well with icing profiles obtained in wind tunnel experiments and in many cases are better than those predicted by the models available.