Three-dimensional numerical simulation of electrothermal deicing systems