Most icing research focuses on the high Reynolds number regime and manned aviation. Information on icing at low Reynolds numbers, as it is encountered by wind turbines and unmanned aerial vehicles, is less available, and few experimental datasets exist that can be used for validation of numerical tools. This study investigated the aerodynamic performance degradation on an S826 airfoil with 3D-printed ice shapes at Reynolds numbers Re = 2 × 105, 4 × 105, and 6 × 105. Three ice geometries were obtained from icing wind tunnel experiments, and an additional three geometries were generated with LEWICE. Experimental measurements of lift, drag, and pressure on the clean and iced airfoils have been conducted in the low-speed wind tunnel at the Norwegian University of Science and Technology. The results showed that the icing performance penalty correlated to the complexity of the ice geometry. The experimental data were compared to computational fluid dynamics (CFD) simulations with the RANS solver FENSAP. Simulations were performed with two turbulence models (Spalart Allmaras and Menter’s k-ω SST). The simulation data showed good fidelity for the clean and streamlined icing cases but had limitations for complex ice shapes and stall.