The thermal factor is the main reason for winter ice cover with a low Froude number flow, and the heat transfer to narrow and deep river banks accelerates ice cover formation and ice thickness change. The freezing of water flow to freezing thickening is a nonisothermal-flow phase transition process coupled with the water flow temperature, environment and riverbank. Here, the Nusselt number and viscous dissipation are used to consider the flow velocity influence on icing, and a thermodynamic model of static ice cover horizontal distribution considering riverbed heat transfer is established. The initial ice time, freezing time and static ice cover thickness formed by static and dynamic water calculated by the model were consistent with measured data. The model reflects the horizontal growth process of the static ice cover, which was significant for narrow and deep channels. The horizontal distribution of the static ice cover was thin in the center and thick on both sides. The maximum horizontal thickness difference of −20 °C indoor freezing for 24 h reached 15% of the central ice thickness. Compared with the degree-day method for calculating ice thickness, the numerical model and dimensionless formula better reflect the growth law and horizontal distribution characteristics of static ice cover and provide a theoretical basis for safe water conveyance under ice cover in winter and ice cover formation in reservoirs and lakes in cold regions.
Read full abstract