The design of power transmission lines requires a knowledge of combined wind and ice loading and of the dynamic behavior of wires loaded with ice accretion. The calculation of the wind forces, in turn, imposes a need for a more detailed computer model for determining glaze accretion shape. For this purpose, a computer model of glaze accretion on wires, was developed. It is based on experimental results in the area of ice accretion on wires, as well as on results in the related field of the glaze ice accretion on airfoils. The model incorporates the time dependent on feedback between the growing accretion and the air stream, the variation of the heat transfer coefficient around the cylinder, and the surface runback of water. The main components of the model are the computation of the air flow field, the computation of the impingement water at the control volume level, the solving of the heat balance equation, and the computation of the accretion shape on the wire. The surface air velocity is obtained through the solution of the potential flow around the iced wire and wake, followed by the integration on the surface of the laminar boundary layer. The water flux is computed in each control Volume down to the separation point. The heat balance equation derived from the energy equation is solved to determine the freezing fraction and the resulting modified ice surface geometry.