Nomenclature
-
Cp Specific heat, kJ (kg K)−1
-
g Gravitational acceleration, m s−2
-
L Latent heat of fusion of ice, J kg−1
-
S Concentration of salt, % by weight
-
t Time, min
-
Τ Temperature, °C
-
Tfp Melting temperature of ice, °C
-
ΔT Temperature difference (T∞ − Tfp)), °C
-
Vi Melting velocity, m s−1
-
x Distance from bottom end of the ice plate, m
-
y Normal distance from the ice surface, m
-
β Coefficient of thermal expansion, °C−1
-
β* Coefficient of concentration expansion, % −1
-
δc Thickness of concentration boundary layer, m
-
δTH Thickness of thermal boundary layer, m
-
δrev Thickness of velocity inverse layer, m
-
δv Thickness of velocity boundary layer, m
-
θ Inclination angle from horizontal, deg
-
λ Thermal conductivity, W (m K)−1
-
v Kinematic viscosity, m2s−1
-
ρ Density, kg m−3
Subscripts
-
m Average
-
0 Ice surface
-
∞ Ambient
-
i Ice
Introduction
When glaciers move into the sea and melt, they exert effects on both temperature and concentration fields in the sea, and consequently may affect marine biological systems and climates.
Melting characteristics of such ice blocks as those from Arctic glaciers depend closely on the flow field around them. In the quiescent fluid, free convection which is promoted by both temperature and concentration fields exerts considerable influences on the melting characteristics. Experimental and analytical studies on the melting characteristics of vertical ice plates in saline water have been carried out by Reference Josberger and MartinJosberger and Martin (1981). However, when an ice block actually melts in sea water, its shape and orientation change continuously.
In the present study, an investigation of the melting heat-transfer characteristics from an ice block immersed in quiescent saline water was performed both experimentally and analytically to obtain basic data applicable to the melting of ice blocks in the sea.
Experimental Apparatus and Procedures
A bubble-free ice block, shaped in prescribed dimensions (60 mm in diameter for an ice cylinder and 250 mm × 400 mm × 40 mm for an ice plate), was installed in an insulated saline-water vessel with a lucite plate, 800 mm wide, 270 mm long and 910 mm high. The block was attached to a cooled tube or plate, the temperatures of which were lowered so that the ice block stuck to it.
The experiments were carried out in saline water of a concentration 3.5% by weight. Thermophysical properties of the water were determined as functions of temperature and salinity (Reference Kukulka, Gebhart and MollendorfKukulka and others, 1987). Initial temperature of the water and inclination angle for the ice plates were set as one of the parameters. Just after the ice block, which was initially kept uniform at an equilibrium temperature for 3.5% concentration, was immersed in the water, observations of the melting behavior and visualizations of flow patterns were carried out, and photographs were taken continuously for a prescribed period. The local melting heat-transfer coefficient was determined by the local melting velocity, which was evaluated by photographically recorded ice contours.
Analysis
Assumptions and basic equations
The following assumptions have been introduced in performing the analysis for the case of the downward-facing ice plate:
-
Fields of velocity, temperature, concentration and moving-boundary velocity are independent of time,
-
Temperature and concentration at the ice surface are uniform.
-
Temperature within the ice is uniformly equal to the ice-surface temperature.
-
Flow is two-dimensional, laminar and incompressible.
-
Density depends on both the temperature and concentration; all other properties are constant.
The x- and y-axes are set along the ice surface and the normal directions, respectively, as shown in Figure 1. The following equations are obtained from the above assumptions. The continuity equation is
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170922041458258-0926:S0260305500011095:S0260305500011095_eqn1.gif?pub-status=live)
The momentum and energy equations are
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923143748-28003-mediumThumb-S0260305500011095_eqn2.jpg?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170922041458258-0926:S0260305500011095:S0260305500011095_eqn3.gif?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923143748-15491-mediumThumb-S0260305500011095_fig1g.jpg?pub-status=live)
Fig. 1. Coordinate system in analysis.
The diffusion equation is
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170922041458258-0926:S0260305500011095:S0260305500011095_eqn4.gif?pub-status=live)
Boundary conditions
From heat, mass and concentration fluxes at the ice surface, the following boundary conditions are given as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170922041458258-0926:S0260305500011095:S0260305500011095_eqn5.gif?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170922041458258-0926:S0260305500011095:S0260305500011095_eqn6.gif?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170922041458258-0926:S0260305500011095:S0260305500011095_eqn7.gif?pub-status=live)
From the condition of non-slip, uniform temperature and uniform concentration at the ice surface (y = 0),
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170922041458258-0926:S0260305500011095:S0260305500011095_eqn8.gif?pub-status=live)
From the condition for each boundary layer, and at the interface between double boundary layers,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170922041458258-0926:S0260305500011095:S0260305500011095_eqn9.gif?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170922041458258-0926:S0260305500011095:S0260305500011095_eqn10.gif?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170922041458258-0926:S0260305500011095:S0260305500011095_eqn11.gif?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20170922041458258-0926:S0260305500011095:S0260305500011095_eqn12.gif?pub-status=live)
Equations (1)–(4) were solved under the above boundary conditions using a profile method. Melting rate was determined numerically from calculated results of temperature, concentration and velocity fields.
Results and Discussion
Melting characteristics of ice and flow patterns
Figures 2 and 3 show the flow fields of the ambient fluid and the shapes of both downward- and upward-facing ice plates after 24 min had elapsed from the start of the experiment. For both cases, the figures show that the thickness of the lower part of the ice plate decreases as the angle of inclination increases, while the average thickness of the ice plate is almost independent of the inclination angle.
As for the downward-facing horizontal ice plate (Fig. 2a), the ambient fluid approaches from both the right and left sides to the center, where it changes into downward flow. Though the turbulence of the ambient fluid flow-increases, the thickness of the boundary layer of the meltwater increases as the inclination angle increases to the vertical. For the upward-facing plate, the ambient fluid was taken into the upward flow of the meltwater from both the right and left sides; this flow decreases, while the upward flow of the meltwater causes an increase in its velocity. This may be the reason why the inclination angle exerts little effect on the melting rate.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923143748-89899-mediumThumb-S0260305500011095_fig2g.jpg?pub-status=live)
Fig. 2. Effect of inclination angle on melting characteristics (downward-facing): T∞ = 12.0°C, 24 min elapsed; a, θ = 0°; b, θ = 15°; c, θ= 30°; d, θ = 60°.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923143748-83854-mediumThumb-S0260305500011095_fig3g.jpg?pub-status=live)
Fig. 3. Fig. 3. Effect of inclination angle on melting characteristics (upward-facing): T∞ = 12.0°C, 24 min elapsed; a, θ = 0°; b, θ = 30°; c, θ = 60°; d, θ = 75°.
Figures 4 and 5, respectively, show variations of the ice-cylinder profiles and schematics of flow patterns. The deviation of the ice-liquid interface from a circular shape results from strong free convection effects, A representative ice profile through melting in all the experiments is that as the melting advances, a step (Fig. 5) forms at the upper part of the ice cylinder, and the flat area near the top of the cylinder tends to increase in size with an increase in far-field temperature T∞.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923143748-53176-mediumThumb-S0260305500011095_fig4g.jpg?pub-status=live)
Fig. 4. Melting ice profile and flow pattern around an ice cylinder : a, T∞ = 2.8°C, 15 min elapsed: b, T∞ = 6.0°C, 15 min elapsed; c, T∞ = 12.2°C, 2 min elapsed; d, T∞ = 18.0°C, 2 min elapsed.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923143748-62542-mediumThumb-S0260305500011095_fig5g.jpg?pub-status=live)
Fig. 5. Flow model around a melting ice cylinder
Characteristics of the ice surface
As in the downward-facing plate, grooves were observed on the ice surface along the flow direction for the upper part of the ice plate, and many small dimples on the ice surface of the lower part (Fig. 6). They are considered to have formed because of secondary flow in the boundary layer of the meltwater.
The effects of the ambient fluid temperature on the melting-surface characteristics of the upward-facing ice plate are shown in Figure 7. For lower ambient
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923143748-52387-mediumThumb-S0260305500011095_fig6g.jpg?pub-status=live)
Fig. 6. Fig. 6. Characteristics of melting ice surface (downward-facing): θ = 30°; a, upper part of inclined ice plate; b. lower part of inclined ice plate.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923143748-20795-mediumThumb-S0260305500011095_fig7g.jpg?pub-status=live)
Fig. 7. Characteristics of melting ice surface (upward-facing): θ = 15°; a, T∞ = 1.8°C; b, T∞ = 18.0°C.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923143748-80095-mediumThumb-S0260305500011095_fig8g.jpg?pub-status=live)
Fig. 8. Characteristics of melting ice surface (ice cylinder): a, T∞ = 6.0°C, 15 min elapsed; b, T∞ = 17.0°C, 20 min elapsed.
temperatures (Fig. 7a), a lot of small but peculiar dimples are clearly observed on the ice surface. As the ambient temperature increases, these dimples gradually change into grooves along the flow direction (Fig. 6b). Dimples and grooves are formed because of the secondary flow which might be produced by the ambient fluid captured into the meltwater flow.
In the case of the ice cylinder, observational results of the ice-surface characteristics are shown in Figure 8. The surface develops longitudinal grooves with a regular amplitude at that part from step to separation point (Fig. 8a). On the other hand, in the part below the transition point, the melting ice surface is smooth for the lower ambient temperature (Fig. 8b).
Melt-rate characteristics
Experimental results for average melting velocity in both downward- and upward-facing ice plates are shown in Figures 9 and 10. respectively. For comparison, results for the vertical ice plate (Reference Johnson and MollendorfJohnson and Mollendorf, 1984) are indicated by a dotted line in the figures, and analytical results are presented by a solid line in Figure 9. For pure water, melting velocity for an upward-facing horizontal ice plate (Reference Lunardini, Zisson and YenLunardini and others, 1986, p. 23–24) is
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923143748-29813-mediumThumb-S0260305500011095_fig9g.jpg?pub-status=live)
Fig. 9. Melting velocity (downward-facing ice plate).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170923143748-41691-mediumThumb-S0260305500011095_fig10g.jpg?pub-status=live)
Fig. 10. Melting velocity (upward-facing ice plate).
indicated by a dashed–dotted line in Figure 10. These figures show that, for saline water, the melting rate increases as the temperature difference increases, but is almost independent of the inclination angle.
On the contrary, analytical results of melting rate increase with increasing inclination angle, which agrees well with the experimental results for ΔT smaller than 8°C. For ΔΤ greater than 8°C, analytical results are lower than the experimental ones. This is explained by the fact that the effect of the inclination angle on the buoyancy term alone was taken into consideration in the analysis, and that the assumption of laminar flow is no more valid for greater temperature differences.
Conclusions
From an examination of melting characteristics of ice blocks immersed in quiescent saline water, determined both experimentally and analytically, the following conclusions may be drawn within the parameter range covered in the present study:
1. For both the ice cylinder and the ice plate, melt rate increases monotonically with ambient temperature, but is almost independent of the inclination angle of the ice plate.
2. Melting characteristics of an ice cylinder depend closely on the ambient temperature.