This paper presents recent improvements in CHAMPS (Chapel Multi-Physics Software) developed at Polytechnique Montréal with emphasis on the aero-icing capabilities. This software, written in the Chapel language, is designed to simulate two-dimensional and three-dimensional multi-physics phenomena involving aerodynamics such as ice accretion and fluid-structure interactions through an Unstructured Finite-Volume Unsteady Reynolds-Averaged Navier-Stokes (URANS) realm. It is programmed using the open-source Chapel language, which facilitates parallel computing on laptops, desktops and High-Performance Computing (HPC) platforms. The basic approach to model ice shapes involves several modules with a segregated strategy: aerodynamic, droplet, thermodynamic and geometric. Four newly implemented state-of-the-art features expand the aero-icing suite, CHAMPS-ICE. A first feature is the implementation of a transitional turbulence model to enhance the physical representation of the boundary layer state near the airfoil (1). Furthermore, a local ice roughness model is added (2) to properly quantify the position of the laminar-turbulent transition which influences the surface convective heat transfer. As for the droplet module, it considers the first and second-order effects of Supercooled Large Droplets (SLD) such as droplet splashing and deformation (3). This extension brings the collection efficiency amplitude and impingement limits closer to experimental data. Another development addresses the stochasticity of the ice accretion process using an advancing front technique (4). Those newly implemented features are discussed and validated on 2D and 3D rime and glaze ice cases taken from the American Institute of Aeronautics and Astronautics (AIAA) Ice Prediction Workshops and the European Ice Genesis project.