Ice accretion (icing) on aircraft surfaces is a significant safety risk through airfoil shape modification and reduction of aerodynamic efficiency. This process occurs when an aircraft flies through clouds of supercooled water droplets that freeze upon impact on exposed surfaces. To counter this hazard, electro-thermal de-icing systems integrate heaters in critical regions to melt ice and reduce performance losses. In this study, a multiphysics computational model is used to simulate ice accretion and electro-thermal de-icing on a NACA-0012 airfoil, accounting for factors such as airflow, droplet impingement, phase changes, and heat conduction. The model's predictions are validated against experimental data, confirming its accuracy. A cyclic electro-thermal ice protection system (ETIPS) is then tested under both standard and severe supercooled large droplet (SLD) conditions, examining how droplet size and angle of attack affect de-icing performance. Simulations without an active de-icing system show severe aerodynamic degradation, including an 11.1% loss of lift and a 48.2% increase in drag at a 12°∘ angle of attack. For large droplets (median 200 μm ), the drag coefficient increases by 36.5%. Under harsh icing conditions, the effectiveness of the de-icing system is found to depend on droplet size, angle of attack, and heater placement. Even with continuous heater operation, ice continues to accumulate on the leading edge at higher angles of attack. While the ETIPS performs effectively against large droplets in its protected zones, unheated regions experience significant ice buildup (especially with 200 μm droplets), indicating that additional or extended heaters may be necessary to ensure complete protection in extreme conditions.