<p>This article presents a physics-mathematical model for the numerical simulation of ice formation on aircraft surfaces under standard atmospheric conditions, incorporating surface roughness arising during icing processes. The resulting roughness is described using an equivalent sand-grain roughness model and its modifications. Two modifications to the equivalent sand-grain roughness model are presented, which automatically calculate the height of the roughness based on inflow air parameters and water droplet characteristics. These models are applied in the numerical solution of airflow around the NACA0012 wing profile, both with and without accounting for surface roughness. The results of this work allow for assessing the applicability of the roughness model across a range of inflow temperatures. Results demonstrate that incorporating roughness leads to the formation of ice accretions with similar shapes and locations as observed in experiments. Overall, accounting for surface roughness in numerical icing simulations of aircraft enables the prediction of ice morphologies that more closely align with experimental data compared to results obtained without considering roughness.</p>