|
The present study aims to investigate both numerically and experimentally thermalfluid phenomenon during melting of ice in a vertical rectangular enclosure of AR=2. Transient numerical simulations have been undertaken to understand the effect of Rayleigh numbers in the range of 106∼6×107 with Ste=0.101 and Sc=0.1, 0.5. Experimental visualization of temperature and velocity fields has been conducted for Ra=6×107, Ste=0.101, Sc=0.1 and 0.5. Numerical Results indicate that the increase of Rayleigh number enlarges the melting ratio and average heat flux of vertical wall and tends to make convective flow in liquid zone more unstable. In the case of Sc=0.5, the average heat flux of the hot wall reaches a minimum value first and then asymptotically approaches to steady state when Ra≧4×106. With Sc=0.1 and Ra≧5×106, the average heat flux of the hot wall features a local minimum and then maximum value before reaching a steady state. At Ra=6×107, an oscillatory natural convection was detected in liquid zone. For Sc=0.1, there exists the unsteadiness between the left and right re-circulating flows in the liquid zone after Fo≧0.15. During the time interval from Fo=0.117 to 0.158 under Sc=0.5, an oscillatory convective flows with a warily moving maximum density contour has been predicted. Similar as oscillatory convection but at higher frequency and larger amplitude was experimentally observed during early stage of the melting process for Sc=0.1 and 0.5.
|