1. Introduction
Bridges are critical for connecting people and goods in the transportation system [
1,
2]. Freezing rain has a serious impact on daily human life, such as the destruction of infrastructure facilities, the reduction of agriculture outputs, and even the malfunction of aircraft [
3,
4]. In 2008, a massive freezing rain and snowstorm stuck southern China, causing severe damage to a vast number of buildings, bridges, and power units. Freezing raindrops are a kind of super-cooled water droplet that exhibit strong cohesion when they impact an object surface below 0 °C. The attached ice on a structure surface could change its aerodynamic characteristics and cause unexpected vibrations or damages. The wind resistance stability of cable after icing, especially under the action of strong wind, is important for the study of the potential impact of icing on bridges [
5,
6].
Cable-supported bridges are efficient types of bridges in mountainous areas [
7,
8]. However, rime ice and glaze ice climate weathers often occur in the mountains, which has a huge impact on the safety of these structures [
9]. Installed structural health monitoring systems [
10,
11] could monitor the change of temperature to warn stakeholders [
12]. The hazards of freezing rain to bridge cables include significantly increasing the self-weight of components and changing structural dynamic characteristics. In addition, icing leads to changes in the designed aerodynamic shape, which also significantly affects the structure’s local wind vibration stability, and induces large vibrations such as cable galloping vibration [
13,
14]. Therefore, the mechanism of ice accumulation and final iced geometries are key issues in research on the aerodynamic performance of the cable system bridge after icing.
Previous research has investigated ice accumulation on power supply cables [
15,
16]. However, bridge cables cannot rotate, considering the stiffness of a bridge cable is far greater than that of a wireline cable [
17]. Thus, these conclusions cannot be applied to the ice accumulation on bridge cables. Demartino et al. [
18], Demartino and Ricciardelli [
19], and Górski et al. [
20] studied the ice accumulation phenomenon and the iced shapes of vertical and inclined cables. It was found that the runback water builds icicles, and a thin ice accumulation is formed. Koss et al. [
21] conducted the experiment of ice accumulation on bridge cables and investigated the final iced shape and aerodynamics coefficients under different wind speed, temperature, and inclined angle. Since it is difficult to obtain a parametric model of iced shapes from experiments, numerical simulation plays an important role in the accurate prediction of ice accumulation geometries.
The common numerical model is based on the Messinger model, which has simplified assumptions on water film behavior [
22]. The software FENSAP-ICE in Canada [
23] considered air shear stress as the main driving force of water film flow. Myers et al. [
24,
25] developed a water film icing model based on lubrication theory and assumed the water film flow direction was affected by the shear stress, gravity, pressure gradients, and surface tension. Leng et al. [
26] investigated that the dynamics of water film as affected by wind shear and gravity in the laboratory. While these studies are focused on in-cloud icing simulation of the airfoil [
27], the results cannot be directly applied to freezing rain accumulation on bridge cable due to the significant difference in wind speed and temperature. However, the assumptions about water flow directions provide very important references.
The amount of icing and the shape of icing are two important issues in ice accumulation. Porcu [
28] established a two-dimensional and three-dimensional random icing model of the cable wires, and reproduced the process of random ice formation on the circular protrusions. Myers [
29] established a mathematical model of icing growth of numerous shapes, and showed that there existed a similar pattern in icing growth among these circular cross-sections. EB Lébatto [
30] also conducted research on the impact of freezing rain icing on wire cables, which revealed that the ice accumulations were more easily formed in rime climate weather, and the final iced shapes varied for different wires. However, considering that the rigidity of cables in cable-supported bridges is much greater than that of power wires, and are difficult to be twisted due to their huge self-weight, these conclusions on power transmission wires cannot be simply applied to the cables of cable-supported bridges.
Recently, some researchers have focused on the effect of ice accumulation on bridge hangers [
31] and circular cylinders [
32,
33]. Koss and Lund [
34] revealed that iced cables in freezing rain weather have a potential for galloping instability, and the significant changes to aerodynamic cable forces should be paid more attention [
35,
36]. Xu [
37] revealed the aerodynamic characteristics of pipelines under the glaze ice condition. The aerodynamic characteristics of iced bridge hangers were investigated by Guo et al. [
38], who described the ice accumulation progress and the final thicknesses of two common types of iced shape. The cables on the cable-supported bridges are huge, and more complicated than circular cylinders on other types of bridges. Therefore, it is necessary to have a more realistic simulation of the ice accumulation of bridge cables to reveal the final iced shapes and their impacts on the bridge’s wind-resistant performance.
Aiming at two typical freezing rains, this paper studied the mechanism and final geometries of the ice accumulation on bridge cables.
Section 2 describes the numerical simulation framework of ice accumulation.
Section 3 introduces the three-dimensional Messinger theory into the numerical models to simulate the ice accumulation progress of bridge cables, which could effectively handle the problem of unfrozen water flow direction.
Section 4 contains the numerical setups of rime and glaze ice climate conditions.
Section 5 to
Section 7 show the simulated iced geometries on bridge cable in different climate conditions, and the aerodynamic forces, critical wind speed and wind attack angle of these iced cables are also studied. Finally, appropriate conclusions are drawn in
Section 8.
2. Numerical Simulation Framework of Ice Accumulation
The numerical simulation of ice accumulation includes the following steps:
- (1)
Building the geometric model of the object;
- (2)
Choosing the appropriate numerical simulation method, including the multiphase flow model and the turbulence model;
- (3)
Embedding the modified Messinger model into AYSYS FLUENT 19.0. After the solution of the flow field is completed, the ice-covered mass on the element surface can be obtained;
- (4)
Based on the ice density and its relationship to temperature, the ice thickness can be obtained, and the iced geometry and surface grid can be reconstructed;
- (5)
Determining whether the calculation has converged. If it is converged, the program will stop the simulation; otherwise, return to step (3).
Based on the conservation equations of the classical Messinger model and the software AYSYS FLUENT 19.0 [
39], a three-dimensional simulation model that can simulate the ice accumulation process on bridge cable and analyze multiphase flow field was established. The flow direction of unfrozen water was determined by considering the air shear stress and gravity, and the three-dimensional Messinger equation was solved by the iteration method. Moreover, the whole process was programmed and embedded into the commercial software for ice accumulation analysis. According to this ice accumulation process, the ice thickness can be obtained, and the geometry of the bridge cable reconstructed. Furthermore, the ice-covered geometries of cables at different temperatures and wind speeds under rime and glaze conditions were analyzed individually.
The flow of the whole simulation process is outlined in
Figure 1. Since the emphasis of this paper is to explain the specific numerical simulation method of rime-ice accumulation and glaze-ice accumulation on a three-dimensional element, the key steps (4) and (5) are described in detail.
The Euler method was employed to achieve the task of calculating the two-phase flow using one set of mesh at the same time. The Euler method has a relatively computational efficiently and can simulate the ice accumulation problem of three-dimensional bridge element with complex geometry [
40]. In order to simplify the analysis, the assumptions in the liquid phase were introduced to the Euler method as follows: (a) the droplets were treated as solid spheres; (b) the collision, deformation, and heat transfer inside the liquid phase were not considered; (c) the air turbulent impulse did not affect the movement of the water droplets.
Since the iced geometry of a bridge element continuously changes with increasing ice thickness, dense wall mesh makes it difficult to use the dynamic mesh method, so the k-ε turbulence model with a relatively large near-wall coefficient was selected.
4. Numerical Simulation Setups
Based on the meteorological data of the Yunnan–Guizhou Plateau in China, the distribution law and duration of freezing rain were obtained, and information such as temperature, rainfall intensity and wind speed during the freezing rain period were extracted. This measured environmental information could provide a solid basis for the realistic numerical simulations of rain and ice accumulation on cables.
Rime and glaze are two common types of ice accumulation for freezing rain conditions. Rime is built up of ice (in the form of small white crystals) when small water droplets freeze upon contact with a cold surface. Glaze ice forms when a relatively large amount of liquid water (a raindrop, or several) freezes while in contact with a surface. Although they are the result of freezing rain, or freezing after rain or heavy fog, these two climate conditions and icing simulations should be individually studied due to their variations in water droplet diameter and ambient temperature.
4.1. Rime Climate Condition
The diameter of a rime water droplet is small, about 0.01~0.08 mm. Therefore, the water droplets are often suspended in the air with a lower temperature. Ice accretion under rime weather has a milky white opaque body with bubble voids, loose texture, and an undulating surface with no fixed shape. The boundary conditions should be precisely simulated according to these characteristics. The left entrance is set to be the mixed phase’s velocity inlet, and the right side is defined as the pressure outlet of the mixed-phase, while the other four sides are all set to be symmetrical. All the boundary settings for rime ice simulation are shown in
Figure 6.
The structural hexahedral mesh was employed to improve computational efficiency and quality, and the mesh quantity was densified only near the cable surface. The total mesh number was 1.3 million. The minimum grid volume was 4.4 × 10
−6 m
3, and the maximum grid volume was 9.1 × 10
−4 m
3. The mesh layout is shown in
Figure 7.
To ensure the accuracy of simulation, the time step was set to be 0.01 s. Based on this small time step, it was too difficult to simulate the whole process of ice accumulation, which normally takes several hours. Therefore, to speed up the simulation process, the mesh deformation obtained in each step was multiplied by 1000 for the rime ice case, and the deformation of each grid was limited below 0.028 mm per step to ensure the negative volume mesh did not appear. Finally, there was a 1000-time step for simulating the rime ice accumulation, which corresponded to 10,000 s (2.8 h) in the actual climatic condition. The ice accretion mass and thickness were generally logarithmically related with the time step. The ice accretion grew rapidly at the beginning, slowing down over time, and the growth was not significant at the set-up step.
4.2. Glaze Climate Condition
The meteorological conditions of glaze weather are significantly different from those of rime weather. The water content of the air in glaze weather is about 2.5 × 10
−7, but the diameter of the water droplets is larger, generally in the range of 0.1~0.5 mm. Ice accretion under glaze conditions has a transparent glass body and a hard texture which is not easy to break, and a strong adhesion to the object surface. Therefore, it can cause major damage to bridge cables. In the glaze ice case, since the two-phase flow was coupled at the entrance and very complex to express, the left side was only set as the air velocity entrance to simulate the natural wind and the upper surface was set as the velocity inlet of the water droplets to simulate the rainfall. Moreover, the right side was the mixed-phase pressure outlet, and the other three sides were set to be symmetrical. The boundary conditions for glaze ice are shown in
Figure 8.
The mesh settings and qualities were nearly same as that in the rime climate conditions, with a total mesh number of 1.3 million, a minimum grid volume of 4.4 × 10
−6 m
3, and a maximum grid volume of 9.1 × 10
−4 m
3. The mesh layout is shown in
Figure 9. The calculation time step was set to be 0.01 s. However, for the glaze ice case, the deformation of each analysis step was multiplied by 200 and restricted so as not to exceed 0.08 mm. The total number of time steps was 5000, which corresponded to 10,000 s (2.8 h) for actual time. Similar to the rime condition, the ice accretion of the glaze was generally logarithmically related to the time step. The growth of ice accretion was not significant at the set-up step.
4.3. Validation of Numerical Method
In order to verify the proposed three-dimensional ice accretion model, the numerical results were compared with the cylindrical icing experiment results by Koss et al. [
21].
The experimental temperatures were −5 °C, −3 °C and −1 °C, and the wind speed was 30 m/s, which was similar to the rime climate condition. The cylinder diameter was 3.81 cm, and 50 cm in length. The numerical simulation was under the same condition for the calculation, and
Figure 10 shows the results of ice accretion on circular cylinders after 30 min. When the temperature was −3 °C and −1 °C, the numerical simulation results were generally consistent with the icing experiment. When the temperature was −5 °C, the numerical results had little differences with the experiment. This was due to the certain deviation of convection heat dissipation coefficient in the simulation and the actual results. Considering that the temperature of freezing rain weather is generally between −0 °C and −5 °C, it is reasonable to use this method to simulate the ice accretion on bridge cables.