^{1}

^{1}

^{2}

^{*}

^{1}

^{1}

^{2}

During winter, ice jams develop when floating ice blocks accumulate in rivers. Ice jams can dramatically decrease in the capacity of flow in a river and can cause ice flooding due to increase in water level. Submergence of floating ice blocks in front of ice cover is critical for the development of an ice jam. In this study, the effect of the rotation angle of ice blocks on the submergence of ice block was assessed. The impacts of both the drag force caused by the flow and the hydraulic pressure force on the rotation of ice block were studied. Considering both the maximum moment for anti-overturn of an ice block, and the associated rotation angle
θ
_{1}
, equations for describing the criteria for ice block entrainment in front of ice cover have been derived. On the basis of the theorem for moment equilibrium, relating the moment acting on a horizontal ice block with the maximum anti-overturn moment of an ice block, the criteria for assessing the overturn-and-submergence of an ice block have been proposed. To verify results using the derived equations for calculating the critical flow velocity for ice block submergence in front of ice cover, data was collected from flume experiments in the laboratory. Experiments have been conducted using different sizes of ice block under different flow conditions in a flume which is 26.68 m long, 0.40 m wide, and 0.6 m deep. Model ice blocks were
made of polypropylene
and have
nearly the same as the mass density of the nat
ural ice. Using proposed method for assessing ice block submergence in front of ice cover, calculated critical flow velocities agree well with those of experi
ments.

In winter, ice pans may form in northern rivers. Depending on velocity of flowing water, floating ice pans or ice blocks can be entrained and submerged in front of the ice cover. Throughout winter period, the submerged accumulate under ice cover. As a result, an ice jam forms and the thickness of ice jam will be increased during the accumulation process. Under such a condition, the cross section under the ice jam for passing flow decreases, which significantly increases water levels upstream of the ice-jammed section compared to under open channel flow conditions. Upstream of the ice jam section, flooding may be caused by a dramatic increase in water levels. An example of ice flooding occurred in 1982 in the Hequ Reach of the Yellow River [

River ice jams in northern rivers can also affect the safety of hydraulic structures such as pump stations and hydropower stations, and may cause loss of life. Therefore, studying the mechanism of ice jam development in northern rivers, such as critical conditions for the ice block submergence in front of ice cover, is important.

Research work regarding critical conditions for ice block submergence in front of ice cover has been carried out by some researchers. For example, a formula was derived to determine the critical flow velocity (V_{c}) for the entrainment of ice blocks in front of an ice cover [

V c = ( 0.035 g L ) 1 / 2 (1)

where, V_{c} is the critical flow velocity; g is the gravitational acceleration; L is the length of the ice block.

Kivisild stated that the critical flow Froude number for the entrainment of ice blocks in front of an ice cover was 0.08 [

F r = 2 ( t / h ) ( 1 − ρ i / ρ w ) C s + [ 1 / ( 1 − t / h ) ] − ( 1 + β ) (2)

where, t is the ice block thickness; h is the upstream water depth of the block; ρ_{i} is the mass density of ice block; ρ_{w} is the mass density of water; C_{s} is the approaching flow velocity coefficient, C_{s} > 1.0; β is the velocity head coefficient under the ice cover.

Using experimental data collected by other researchers [

V c [ g t ( 1 − ρ i / ρ w ) ] 1 / 2 = 2 ( 1 − t / h ) [ 5 − 3 ( 1 − t / h ) 2 ] 1 / 2 (3)

The effects of flow Froude number on the formation of ice jam and ice accumulation have been also studied by researchers [

F c = 0.271 ( t h ) − 0.101 ( t L ) − 0.088 (4)

Wang et al used laboratory experiments to study the critical condition for the entrainment of ice blocks in front of an ice cover [

V c g t = 0.1034 ( t L ) − 0.113 ( b L ) 0.260 ( t h ) − 0.492 (5)

To assess the influence of the shape of ice block leading edge on its stability, experiments were conducted and the pressure difference distribution of the bottom surface of an ice block was measured with different leading edge shapes as well as the stability of rectangular ice block using the moment balance method [

Based on the experimental data [

V c [ g t ( 1 − ρ i / ρ w ) ] 1 / 2 = k 2 ( 1 − t / h ) [ 5 − 3 ( 1 − t / h ) 2 ] 1 / 2 (6)

The coefficient k ranges from 1.15 to 1.35. The coefficient k takes the minimum value if the frontal cross section of the ice block has a quasi-rectangular shape [

Considering the drag force and the Venturi pressure effect, the critical Froude number for submergence of an ice block in front of ice cover was studied using the method of the equilibrium of moments [

To date, research results regarding the criteria for entrainment of ice block in front of ice cover are either limited to field measurements in a specific river or to laboratory experiments. On the other side, some researchers studied the criteria for entrainment of ice block on the basis of theory of mechanics. Results based on theory of mechanics are relatively more comprehend. However, the rotation angle of the ice block has not been considered in these studies when the maximum anti-overturn moment is generated. The overturned entrainment process of ice block is shown in

As shown in _{1}, the moment for anti-overturn of an ice block increases approximately linearly [_{1} to θ_{2}, the moment for anti-overturn of an ice block first increases slowly until a maximum is reached. Following the maximum, the moment for anti-overturn of an ice block decreased gradually. This means that the maximum moment for anti-overturn of an ice block occurs when the angle between ice block and ice cover is θ_{1} < θ < θ_{2}, but very close to θ_{1}. The maximum moment for anti-overturn of an ice block (R_{max}) could be calculated by Equation (7). After that, as the angle between ice block and ice cover increases from θ_{2} to π/2, the moment for anti-overturn of an ice block decreases.

R max = ( ρ w − ρ i ) g t b L 2 2 (7)

The moment for anti-overturn of an ice block changes with the rotation angle of the ice block. For an ice block to rotate, the moment acting on an ice block must be greater than the maximum moment for anti-overturn of an ice block. Therefore, to study the critical conditions for an ice block to be entrained and submerged in front of ice cover, the moment balance method should be followed. The moment analysis should be conducted provided that the ice block is rotated to an angle of θ_{1} (the position for the maximum anti-overturn moment).

When an ice block floats in the surface of flowing water, it is subjected to buoyancy force, gravitational force, drag force caused by flowing water, hydraulic pressure force, etc. The pressure acting on the ice block in front of the ice cover is the lowest force. The pressure force is mainly caused by the venturi effect and the interaction of the ice block and frontal ice cover. When an ice block approaches the front of an ice cover, it will be non-horizontal affected by flowing water in front of the ice cover. The pressure under ice block/cover gradually changes along the direction of flow, due to the decrease of the cross sectional area caused by the ice block/cover, which can be estimated using energy equation. As a consequence, the venturi effect occurs on the ice block. At the leading edge of the ice cover, the local flow velocity increases and accelerates in addition to the flow separation. Thus, the pressure at the leading edge of the ice block/cover reduces to a minimum, resulting in the maximum pressure difference.

As shown in _{1}) can be calculated using following equation [

θ 1 = sin − 1 ⌊ − ( ρ i / ρ w ) ( L / t ) + ( L / t ) 2 + 1 − ( ρ i / ρ w ) 2 ( L / t ) 2 + 1 ⌋ (8)

F 1 = 1 2 ρ w C D A V 1 2 (9)

where, C_{D} is the drag coefficient, V_{1} is the average cross sectional velocity in front of ice cover (m/s), A is the cross sectional area of the ice block which is perpendicular to the flow direction (m^{2}). The drag force acting on an individual ice block can be further expressed as

F 1 = 1.42 [ ( b t / L 2 ) 1 / 3 − 0.26 ] ⋅ 1 2 ρ w V 1 2 ⋅ b t (10)

where, b, t and L are the width, thickness and length of an individual ice block, respectively. Thus, the overturning moment M_{1} generated by the drag force caused by the flowing water is:

M 1 = F 1 t 2 = 1.42 4 ρ w V 1 2 [ ( b t / L 2 ) 1 / 3 − 0.26 ] b t 2 (11)

When the study ice block is floating horizontally on water surface, the decrease of the cross section at the frontal edge of ice block is the same as that at the rear edge. Thus, the pressure decrease generated by the increase in flow velocity is constant in the longitudinal direction (flow direction) along the study ice block. Namely, the overturning moment for the rotating ice block cannot be generated, as shown in _{1}, the force due to the reduction of pressure along the length of the ice block can be described as following,

F 2 = 1 2 ρ w ( V 2 2 − V 3 2 ) b L = 1 2 ρ w V 1 2 ( α 1 2 − α 2 2 ) b L (12)

where, { h 2 = h 1 − h t = h 1 − t cos θ 1 h 3 = h 2 + L sin θ 1 = h 1 − t cos θ 1 + L sin θ 1 α 1 = h 1 / h 2 ; α 2 = h 1 / h 3 .

In which, L is the length of the block (m); t is the thickness of the block (m); h_{t} is the projection length of the upstream submerged ice block in vertical direction (m). h_{1}, h_{2} and h_{3} are respectively the water depth (m) at sections 1-1, 2-2 and 3-3, as shown in _{2} generated from the venturi effect can be determined as following,

M 2 = F 2 L 2 cos θ 1 = 1 4 ρ w V 1 2 ( α 1 2 − α 2 2 ) b L 2 cos θ 1 (13)

At the leading edge of ice block, the force due to the pressure difference can be described as following,

F 3 = Φ ρ w V 2 2 b L = 0.35 Φ ρ w V 1 2 α 1 2 b L (14)

where: Φ stands for a reduction factor. Note: Φ values for: when the length of submerged ice block L < 0.5X_{50}, then, Φ = 1; when 0.5X_{50} < L < 3X_{50}, then Φ = 0.5; when L > 3X_{50}, then Φ = 0.2. In the present study, Φ = 0.5 [_{50} is used to assess the impact of the length of submerged ice block on the process of ice block submergence. It is the location at which the pressure is midway between the initial and final pressure plateau values [_{3} generated by the leading edge effect is,

M 3 = F 3 L 2 cos θ 1 = 0.175 Φ ρ w V 1 2 α 1 2 b L 2 cos θ 1 (15)

In present study, Equation (7) [_{1} to be submerged, the total anti-overturn moment M_{T} acting on the ice block should be,

R max ≤ M T = M 1 + M 2 + M 3 (16)

The total anti-overturn moment M_{T} acting on ice block can be expressed as following,

M T = 1.42 4 ρ w V 1 2 [ ( b t L 2 ) 1 / 3 − 0.26 ] b t 2 + 1 4 ρ w V 1 2 ( α 1 2 − α 2 2 ) b L 2 cos θ 1 + 0.35 2 Φ ρ w V 1 2 α 1 2 b L 2 cos θ 1 (17)

The critical flow velocity can be obtained by simplification of above equation as following,

V c = ( ρ w − ρ i ) g t 0.71 [ ( b t L 2 ) 1 / 3 − 0.26 ] ρ w ( t L ) 2 + 0.50 ρ w ( α 1 2 − α 2 2 ) cos θ 1 + 0.35 Φ ρ w α 1 2 cos θ 1 (18)

To verify results using the derived equations for calculating the critical flow velocity for ice block submergence in front of ice cover, data was collected from flume experiments in the laboratory. As shown in ^{3}, which is nearly the same as the mass density of the natural ice at 917 kg/m^{3}. The model ice blocks had the following dimensions (length × width × height): 0.02 m × 0.02 m × 0.01 m, 0.03 m × 0.03 m × 0.01 m, 0.04 m × 0.04 m × 0.01 m and 0.06 m × 0.06 m × 0.01 m, respectively, the same sizes (the same dimensions) as those used by other researchers [

The submergence process of ice blocks was classified into two states: “some ice blocks submergence” and “massive ice blocks submergence”. The “some ice blocks submergence” refers to that a portion of ice blocks submerged in front of the ice cover, and normally accompanied with the “incomplete rotate” phenomenon, as shown in

Compared with the experimental results, as shown in _{1}. If the ice block in front of ice cover was not rotated at an angle of θ_{1}, the ice block may not be submerged even if the flow velocity exceeds the critical velocity. Under such a condition, the “massive ice blocks submergence” does not occur even if the flow velocity reaches the critical velocity, as shown in

When flow velocity reaches the critical velocity for ice blocks in front of ice cover, changes of other factors such as the turbulence field may lead to an ice block to rotate an angle of θ_{1}, resulting in the “some ice blocks submergence” case. However, to obtain a formula for determining the critical velocity for the case of “massive ice blocks submergence” (instead of “some ice blocks submergence”), ice blocks should be considered floating horizontally on the surface of the water, see

Based on laboratory experiments, it is found that when an ice block began to rotate, the overturning moment acting on the ice block increased compared to that acting on a horizontal ice block [_{1} (critical position), the overturning moment increased by 25%. Namely, the overturning moment acting on ice block under the condition of the critical state (θ_{1} position) was about 1.25 times of that acting on a horizontal ice block. Therefore, one can relate the overturning moment acting on the horizontal ice block to that of under the condition of the critical state (θ_{1} position) as follows,

M T 1.25 = M (19)

where, M is the overturning moment of a horizontal ice block floating on water surface in front of ice cover.

Using Equation (19), the overturning moment acting on ice blocks under the condition of the critical state (θ_{1} position) can be converted into the overturning moment acting on ice blocks floating horizontally on the surface of the water. Following, the relationship between the overturning moment acting on ice

blocks floating horizontally on water surface and the maximum anti-overturning moment can be established. A higher critical flow velocity for submergence of a horizontal ice block (V_{c}_{1}) can be obtained.

Regarding the overturning moment under the condition of the critical state (θ_{1} position), the overturning moment acting on the horizontal ice block can be determined as follows,

M = 0.8 { 0.355 ρ w V 1 2 [ ( b t L 2 ) 1 / 3 − 0.26 ] b t 2 + 0.25 ρ w V 1 2 ( α 1 2 − α 2 2 ) b L 2 cos θ 1 + 0.175 Φ ρ w V 1 2 α 1 2 b L 2 cos θ 1 } (20)

The critical flow velocity for submergence of a horizontal ice block (V_{c}_{1}) is obtained by:

V c 1 = 1.25 ( ρ w − ρ i ) g t 0.71 ρ w [ ( b t L 2 ) 1 / 3 − 0.26 ] ( t L ) 2 + 0.5 ρ w ( α 1 2 − α 2 2 ) cos θ 1 + 0.35 Φ ρ w α 1 2 cos θ 1 (21)

Calculated results using Equation (21) are compared to those of laboratory experiments, as shown in

To date, the criteria for entrainment of an ice block in front of the ice cover did not consider the rotation angle of the ice block when the maximum anti-overturn moment was generated. During the period of the collision/submergence of ice blocks in front of ice cover, the incoming ice blocks will collide with ice blocks in front of ice cover. The kinetic energy of the incoming ice blocks will affect the rotation and submergence of ice blocks in front of ice cover, and provide a variety of possibilities for forming ice jam at the leading edge of ice cover. With respect to the rotation angle of ice block in front of the ice cover, the entrainment of an ice block can be classified as following 3 cases:

1) Case 1: If the rotation angle of ice block B is small (θ < θ_{1}), the incoming ice block C with velocity V collides with ice block B, and will provide a horizontal force where ice block B and ice block C come in contact. Above the center of mass of the ice block, the anti-overturning moment is generated. In other words, the kinetic energy of ice block C was transferred to ice block B as both the potential energy (anti-overturning) and kinetic energy of ice block B turned counterclockwise, as showed in

2) Case 2: If the rotation angle of ice block B is large (θ > θ_{1}), as shown in

3) Case 3: If the rotation angle of ice block B is large enough (θ > θ_{1}), another type of collision and accumulation of ice blocks will happen. As shown in

The fragmentation of the ice block was not taken into account in present study. For instance, if the edges of ice blocks are broken, ice block submergence in front of ice cover will be more complicated.

During winter, the submergence of floating ice blocks in front of ice cover is critical for the development of an ice jam. Ice jams can dramatically decrease in the capacity of flow in a river and can cause ice flooding. In the reported research, the criteria for entrainment of an ice block in front of the ice cover did not consider the rotation angle of the ice block when the maximum anti-overturn moment was generated. In present study, the effect of the rotation angle of ice blocks on the submergence process of ice block was assessed. Considering both the maximum moment for anti-overturn of an ice block, and the associated rotation angle θ_{1}, equations for describing the criteria for the entrainment of ice blocks in front of the ice cover have been derived. It was found that the formula for calculating the critical velocity can only reflect the “some ice blocks submergence” case. In the case of “massive ice block submergence”, ice blocks in front of ice cover are considered floating horizontally on water surface. The overturning moment of the horizontally floated ice block is obtained based on the condition of the critical state (θ_{1} position). Thereafter, the maximum anti-overturning moment of the ice block can be achieved. An equation for calculating the critical velocity V_{c}_{1} reflects the “massive ice blocks submergence” has been derived. To verify results using the derived equations for calculating the critical flow velocity for ice block submergence in front of ice cover, flume experiments were carried out in the laboratory. Using the proposed method for calculating the critical flow velocity of an ice block in front of ice cover, the calculated critical flow velocities agree well with those of laboratory experiments.

This research is supported by the National Key R & D Program of China (Grant No. 2018YFC1508401), the National Natural Science Foundation of China (Grant Nos. 51879065) and Graduate Quality Project of Hefei University of Technology (2018YKC02). The authors are grateful for the financial support.

The following symbols are used in this paper:

A = the cross sectional area of the ice block which is perpendicular to the flow direction;

g = the gravitational acceleration;

b = the width of ice block;

C_{s} = the approaching flow velocity coefficient;

C_{D} = the drag coefficient;

h = the upstream water depth of the block;

h_{t} = the projection length of the upstream submerged ice block in vertical direction;

L = the length of the ice block;

R_{max} = the maximum moment for anti-overturn of an ice block;

t = the ice block thickness;

F_{1} = the drag force caused by the flowing water;

F_{2} = the force due to the reduction of pressure along the length of the ice block;

F_{3} = the force due to the pressure difference;

M_{1} = the overturning moment generated by the drag force;

M_{2} = overturning moment generated from the venturi effect;

M_{3} = the overturning moment generated by the leading edge effect;

M_{T} = the total anti-overturn moment acting on the ice block;

M = the overturning moment of a horizontal ice block floating on water surface;

V_{1} = the average cross sectional velocity in front of ice cover (m/s);

V_{c} = the critical flow velocity;

V_{c}_{1} = thecritical flow velocity for submergence of a horizontal ice block;

ρ_{i} = the mass density of ice block;

ρ_{w} = the mass density of water;

θ_{1} = the position for the maximum anti-overturn moment;

θ = the angle between ice block and ice cover;

β = the velocity head coefficient under the ice cover;

Φ = the reduction factor.

The authors declare no conflicts of interest regarding the publication of this paper.

Wang, J., Yang, W.L., Sui, J.Y., Cheng, T.J., Li, Z.Q. and Linklater, N. (2020) Criteria for Submergence of Ice Blocks in Front of Ice Cover. World Journal of Engineering and Technology, 8, 523-536. https://doi.org/10.4236/wjet.2020.83037