In order to accurately predict the ice accumulation on bridge cables under two typical freezing rain conditions, rime and glaze ice, this paper proposes a numerical simulation framework based on the three-dimensional Messinger theory. Two technical challenges of determining the flow direction of unfrozen water and solving three-dimensional Messinger equations are solved in this research. Based on the outflow, mass was calculated according to the three-dimensional Messinger theory, and the flow direction of unfrozen water in each cell was determined by the resultant force of air shear stress and water film gravity. To solve the three-dimensional equations, an iterative method without finding the stagnation line was introduced. The final iced geometries were determined when the inflow mass ratio was satisfied with the converge criteria. Moreover, this modified numerical model was programmed and embedded into computational fluid software. For both two typical freezing rain conditions, the effects of temperature and wind speed on iced geometries were studied. The aerodynamic characteristics and galloping instability of bridge cables with different iced geometries were also investigated. These preliminary aerodynamic simulations can provide the basis for the wind-induced vibration analysis of the whole structure.

No abstract