Nucleation of supercooled water on solid surfaces

Heterogeneous nucleation of water was investigated using molecular dynamics simulation. Solid with fcc (111) surface was placed at the bottom of a cell consisting of 864 water molecules. ST2 model with NPT ensemble was used. The pressure and temperature were set to be 0.1 MPa and 275 K, respectively. The interaction between water and the solid was based on the equations proposed by Spohr. Exception was made on the lattice constant which was slightly modified to fit with that for ice structure. The shape of the solid surface was considered. It was found that the only one layer of water molecules was adsorbed in a case of a flat surface, whereas ice nucleation occurred by removing some of the atoms from the surface. Spohr's interaction was also modified so that the dipole moment of water became anti-ferroelectric. It was found that the modification increased the ice growth, further. The effect of lattice constant of solid on nucleation was also investigated. It was found that the variation on lattice constant with a few percent from that of ice was acceptable for nucleation, especially on shrinking side. On expanding side, however, it gave some gaps for water molecules to fit in other than that for ice structure, and it prevented the growth of ice. Hence, a guideline for the selection of ice nucleus material was obtained.     

Freezing phenomena of supercooled water under impacts of ultrasonic waves

In order to clarify effects of ultrasonic waves on supercoold water, ultrasonic waves were applied to supercooled water. Frequency of the ultrasonic waves applied was 45 kHz, and the intensities of the waves were 0.13 and 0.28 W cm⁻², respectively. Four cases of experimental conditions were selected; a case of a free surface, a case of an oil-water surface, a case of a free surface with a dipped metal bar and a case of an oil–water surface with a dipped metal bar. For each experimental condition, water was cooled at a constant cooling rate and the ultrasonic wave was applied from underneath until the water in a test tube solidified. It was found that in the case of 0.28 W cm⁻², the ultrasonic wave had distinct effects on freezing of supercooled water for all experimental conditions. On the other hand, in the case of 0.13 W cm⁻², the ultrasonic wave had an effect on freezing of supercooled water only under existence of a free surface and a metal bar.     

Ice crystal growth in supercooled solution

Study on the crystal growth of ice in water-ethyleneglycol solution was carried out, experimentally. The ice crystal, which was seeded on top of the capillary tube, propagated inside the tube slowly and began to grow freely at the tip of the tube in subcooled solution. The outer diameter of the tip of the capillary tube was less than 0.1 mm, which was much smaller than that of other researchers. Hence, considerable reduction of the influence of the existence of a capillary was accomplished and the initial growth of a single crystal was observed, precisely. Under the condition of subcooling of less than 8 K, the shape of the crystal was observed to be different from that of the one in pure water. The velocity of dendrite ice growth and the radius of the curvature of the tip were measured. It was found that after the dendrite ice developed to a certain size, the velocity of the dendritic growth and the radius of curvature were kept steady, and its values were dependent on the degree of subcooling and the concentration of the solution.     

Three-dimensional measurement of ice crystals in frozen dilute solution

A Micro-Slicer Image Processing System (MSIPS) has been applied to observe the ice crystal structures formed in frozen dilute solutions. Several characteristic parameters were also proposed to investigate the three-dimensional (3-D) morphology and distribution of ice crystals, based on their reconstructed images obtained by multi-slicing a frozen sample with the thickness of 5 μm. The values of characteristic parameters were determined for the sample images with the dimension of 530×700×1000 μm. The 3-D morphology of ice crystals was found to be a bundle of continuous or dendrite columns at any freezing condition. The equivalent diameter of ice crystals were in the range of 73–169 μm, and decreased exponentially with increasing freezing rate at the copper cooling plate temperature of −20 to −80 °C. At the T_cp −40 °C, the volumes of ice crystals were in the range of 4.6×10⁴ μm³ to 3.3×10⁷ μm³, and 36 ice columns were counted in the 3-D image.     

An experimental study on ice formation around horizontal long tubes

The results of an experimental study are presented where the growth rate of ice on the outside of cooled copper tubes was studied. The tubes, which were immersed in water in an insulated vessel, were internally cooled by circulating glycol through them.It was found that axial growth rate of ice is distinct at low values of the coolant Reynolds number and short freezing times. The slope of the ice thickness with axial distance showed moderate dependency on time but varied with coolant flow rate, and with Stanton and Biot numbers.A key result from the experiments is the abrupt ice thickness enlargements on the surface of tube bends. This anomaly may be attributed to internal flow disturbances of the coolant, and creation of local eddies inside the bends that enhance growth of ice. The effect was evident for a low Reynolds number (Re = 251.9 and Bi < 1), and fades out for large Reynolds number flows.     

The effect of stationary and sweeping frequency AC electric fields on frost crystal removal on a cold plate

The effect of stationary and sweeping frequency AC electric fields on frost crystals growth and frost control/removal on a cold plate was studied for the first time in this paper. The main results of this study showed that the presence of AC electric fields can greatly affect both the frost crystals growth pattern and mass accumulation on cold surfaces. The ice surface electrical properties and basic electrostatics were used to explain the main findings in this paper. Up to 46% frost reduction was obtained when the electric field frequency spanned 370 Hz to 7.5 kHz while the applied voltage was 14.5 kV. Two different sets of environmental conditions were tested, which showed that the plate temperature placed an important effect on frost crystals growth under electric fields. An optimum application time of the AC electric fields was found based on least frost mass accumulation on the cold plate.     

Shape identification for water–ice interface within the cylindrical capsule in cold storage system by inverse heat transfer method

In this study, the inverse heat transfer method is applied to shape identification for the ice layer within the cylindrical capsule in cold storage system. The approach is constructed by combining the curvilinear grid generation scheme, the direct problem solver, the conjugate gradient optimization method, and the redistribution method. According to the practical condition of freezing ice, shape identification for the water–ice interface based on the data of the outer surface temperature is attempted. Results show that the profile of the water–ice interface is possible to be identified by using the inverse heat transfer approach and the accuracy of the ice shape identification is dependent on the uncertainty of the outer surface temperature data, the Biot number, the thickness of the ice layer, and the geometric configuration as well.     

Experiments on fast nucleation and growth of HCFC141b gas hydrate in static water columns

A major technical issue in the realization of refrigerants gas hydrates energy storage systems is to develop a practical means for rapid hydrate formation. In this paper, the nucleation and growth processes of HCFC141b (CH₃CCl₂F, HCFC141b) gas hydrate in a column of water added with 0.038 wt% sodium dodecylbenzenesulfonate-6 and with a metal rod which was placed in the center of the column has been studied. The water solution column, in a cylindrical glass container, was placed in a thermostatic bath ranged 274.15–280.15 K and under atmospheric pressure. The experimental results show that, comparing to that of pure water and HCFC141b system, the properly placed metal rods combined with proper amount of sodium dodecylbenzenesulfonate-6 considerably promotes the hydrate formation speed and reduces the nucleation time. The growth rate of HCFC141b hydrate is significantly effected by the surrounding bath temperature and the formation can be finished in 8 h with a short nucleation time when the bath temperature ranged 274.15–279.15 K and the iron rod penetrating the interface between the water solution and liquid HCFC141b. The experimental results suggest a new way for fast formation of gas hydrates.     

Effect of axial conduction on the ice crystal growth in laminar falling filmsEffet de la conduction axiale sur la vitesse de formation de cristaux de glace à l'intérieur de films tombants laminaires

A finite volume numerical code has been developed to numerically approximate the rate of ice crystal growth in a laminar falling film flowing down a cooled vertical plate. The governing energy equation contains the phase energy as the source term. Enhancement of heat transfer as a result of suspended ice crystals is accounted for in the use of effective values of thermal conductivity, viscosity, thermal diffusivity, and specific heat as function of volumetric concentration of ice crystals in the falling film. Nusselt number, overall heat transfer coefficients between the fluid and cooled plate, and ice crystal growth rate were calculated for different film thicknesses with and without axial diffusion. Nusselt number and ice crystal growth rates were found to be dependent on film thickness. Axial diffusion effects were found to be negligible for larger film thickness (large flowrate).     

Three-dimensional measurement of ice crystals in frozen beef with a micro-slicer image processing system

A novel technique has been developed for measuring the three-dimensional (3-D) structure and distribution of ice crystals formed in frozen beef by using a micro-slicer image processing system (MSIPS). The system has functions to reconstruct the 3-D image based on the image data of exposed cross-sections obtained by multi-slicing of a frozen sample with the minimum thickness of 1 μm and to display the internal structure as well as an arbitrary cross-section of the sample choosing observation angles. The size and distribution of ice crystals can be determined from the 2-D quantitative information, such as the periphery and area of the crystals. The effects of freezing conditions on the morphology and distribution of the ice crystals were demonstrated quantitatively from the observations of raw beef stained by fluorescent indicator. For the samples frozen at −15 °C, the network structure of ice crystals were observed mainly at intercellular space, having approximately 100 μm in cross-sectional size, while that prepared at −120 °C showed the spherical crystals of 10–20 μm in diameter within the cells. The 3-D image of the sample demonstrated that the growth of ice columns was restricted by the intrinsic structure of muscle fibers. The proposed method provided a new tool to investigate the effects of freezing conditions on the size, morphology and distribution of ice crystals.     

Effect of bubble nuclei on freezing of supercooled water

In order to clarify effects of bubble nuclei on freezing of supercooled water, various kinds of experiments were carried out with invisible sizes of bubbles in supercooled water. Water samples were kept in a test tube and cooled down at a constant cooling rate until the water solidified. The degree of supercooling at freezing was then measured. Two kinds of water surfaces were applied. One was exposure to the atmosphere, and the other was covered with silicone oil. Three kinds of pressure conditions were applied. The first type was atmospheric pressure. The second type was compression up to 6.0 atm. The third type was evacuated down to 0.02 atm. Two holding time periods before starting the experiments were applied. One was 30 min and the other was 24 h. It was found that the degree of supercooling at freezing is high in the case of the free surface compared to the one with oil–water surface. The reason suggested was that the bubbles in the water can be released from the surface in the case of normal atmospheric exposure and trapped at the oil–water interface in the case of water covered with oil. Hence, it was clarified that the freezing of supercooled water is affected by the existence of bubble nuclei.     

Effects of electrode materials on freezing of supercooled water in electric freeze control

In order to clarify effects of electric charge on freezing of supercooled water, experiments were carried out with various kinds of electrodes in supercooled water. Water sample was kept in a test tube and cooled down at a constant cooling rate. When the water sample was maintained under a supercooling state, an electric charge was applied to the water sample with a small electric current. The degree of supercooling was measured continuously. Then the degree of supercooling at freezing was determined. Six kinds of materials were used for electrodes. Those materials were Aluminum, Copper, Argentum, Aurum, Platinum and Carbon. It was found that the effects of electric charge were distinct according to the material used for electrodes. The degree of supercooling at freezing was the lowest in the case of Aluminum. On the other hand, the highest value of the degree of supercooling at freezing was obtained in the case of carbon. The reason for the difference in the degree of supercooling at freezing by six materials was discussed.     

The experimental study on freezing of supercooled water using metallic surface

Freezing of supercooled water on a metallic plate was studied, experimentally. First of all, a gold plated surface was selected as a metallic surface because the surface has very little change in the characteristics of the surface against time. The experiments on freezing under various cooling rates were carried out and the probability of freezing per unit surface area per unit time interval was calculated. It was found that the probability was independent of the cooling rate. Secondly, in order to clarify the effect of oxidation on freezing of supercooled water, an electrolytically polished copper surface was selected and a time variation of the probability of freezing was investigated. A large number of experimental data is required to obtain an accurate value for the probability, but it is impossible to perform it before the characteristics of the surface changes. Hence, in order to cover a wide range of the degree of supercooling at freezing in a short period of time, experiments were carried out under various cooling rates. It was found that the oxidation of the surface restrained the supercooled water on the surface from freezing. By comparing the results with the one for a gold plated surface, a parameter was obtained to express the characteristics of the surface.     

Effects of shapes of electrodes on freezing of supercooled water in electric freeze control

In order to clarify the effects of electric charge on freezing of supercooled water, experiments were carried out. Two kinds of shapes were used for tips of electrodes. One was the sharp end surface. The other was the flat end surface. Aluminum was selected as the material. Water sample was kept in a test tube and cooled down at a constant cooling rate. When the water sample was maintained under a supercooling state, an electric charge was applied to the water sample with a small electric current. The degree of supercooling and the electric current were measured, continuously. Then the degree of supercooling at freezing was determined. It was found that the effects of electric charge were distinct according to the shapes of electrodes. The degree of supercooling at freezing was the lower in the case of flat end surface of the anode than that in the case of sharp end surface of the anode.     

Solidification of phase change material on vertical cylindrical surface in holdup air bubbles

Solidification of phase change material around a vertical cylindrical surface was studied to investigate the performance of ice storage system and stored thermal energy. Air bubbles were generated in the phase change material at various air flow rate as a gas holdup to enhance the heat transfer rate and accelerate the ice layer growth at the solid–liquid interface. The test tube surface was cooled by ethylene glycol–water solution at a flow rate of 40% concentration by weight. The ice layer growth and solidification front velocity at solid–liquid interface were estimated from the temperature–time recorded data of a set of thermocouples fixed in a radial position perpendicular to cooled surface. The ice layer growth at the first instants of solidification process is much higher. Thereafter it decreased gradually according to the increasing of thermal resistance of ice layer. The increasing of ethylene glycol–water solution mass flow rate seems to accelerate the solidification process with small rate. The effect of air bubbles agitation was found to increase the ice layer growth rate and solidification front velocity by about of 20–45%. As a consequence the stored thermal energy was increased by about 55–115% with increasing air bubbles flow according to the attribute of generates turbulence at the solid–liquid interface. The measured data showed that with stirring the bulk water in energy storage tank, the storage time can be reduced by 10–35% of that without stirring.     

Fundamental study on adhesion of ice to cooling solid surface

Many technological troubles are caused by ice adhesion to a cooling solid wall (surface). Therefore, it is urgent to clarify a mechanism of ice adhesion.It is thought that ice adhesion to the cooling wall is governed by heat transfer and interfacial phenomena between ice and the wall. In this study, shearing stress corresponding to adhesion force per unit area to remove ice from the wall surface and some reagents' contact angles on the wall and ice were measured, varying the wall material and its surface state. Moreover, shearing work to remove ice from the wall surface and surface energies of the wall and ice were calculated by the shearing stress and contact angles, respectively. And adhesion energy at an interface between ice and the wall was also calculated by the calculated surface energies. And then, influence of heat transfer and interfacial phenomena on ice adhesion was discussed to clarify the mechanism of ice adhesion.     

Ice-water two-phase flow behavior in ice heat storage systems

In this paper, we describe (1) the characteristics of ice-water two-phase flow created in a supercooling-type ice storage system, (2) the numerical analysis models for both processes of ice storing and melting within ice storage tanks, and (3) the design technologies for ice storage systems based on such predictive models, as well as the operational performance of actual equipment built with these design technologies. As the result of this technological development, it has been confirmed that the predictive approaches for the one dimensional ice accumulation in the given tanks and ice melting with low temperature water enough to air conditioning applications are reasonable and widely usable. Furthermore by the design conditions obtained here it is possible to simplify the construction of ice storage tanks and to realize high reliability in the actual equipment.     

Simulation of a transcritical CO2 heat pump cycle for simultaneous cooling and heating applications

A steady state simulation model has been developed to evaluate the system performance of a transcritical carbon dioxide heat pump for simultaneous heating and cooling. The simulated results are found to be in reasonable agreement with experimental results reported in the literature. Such a system is suitable, for example, in dairy plants where simultaneous cooling at 4 °C and heating at 73 °C are required. The optimal COP was found to be a function of the compressor speed, the coolant inlet temperature to the evaporator and inlet temperature of the fluid to be heated in the gas cooler and compressor discharge pressure. An optimizing study for the best allocation of the fixed total heat exchanger inventory between the evaporator and the gas cooler based on the heat exchanger area has been carried out. Effect of heat transfer in the heat exchangers on system performance has been presented as well. Finally, a novel nomogram has been developed and it is expected to offer useful guidelines for system design and its optimisation.     

Study of GaP single crystal layers grown on GaN by MOCVD

The performance of GaN based devices could possibly be improved by utilizing the good p-type properties of GaP layer and it provides the possibility of the integration of InAlGaN and AlGaInP materials to produce new devices, if high quality GaP compounds can be grown on III-nitride compounds. In this paper, the growth of GaP layers on GaN by metalorganic chemical vapor deposition (MOCVD) has been investigated. The results show that the GaP low temperature buffer layer can provide a high density of nucleation sites for high temperature GaP growth. Using a 40 nm thick GaP buffer layer, a single crystal GaP layer, whose full-width at half-maximum of the (1 1 1) plane measured by double crystal X-ray diffraction is 580″, can be grown on GaN. The V/III ratio plays an important role in the GaP layer growth and an appropriate V/III ratio can improve the quality of GaP layer. The GaP:Mg layer with hole carrier concentration of 4.2 × 10¹⁸ cm⁻³ has been obtained.     

Momentum transfer of ice slurry flows in tubes, modeling

In this paper we present results of the studies of ice slurry flow in horizontal tubes. The possibility of treating the rheological parameters of ice slurry as being those of Bingham fluid was confirmed. The values of parameters (mass fraction, flow velocity) corresponding to the laminar, intermediate and turbulent flow were determined which permits to optimize the flow in the systems working with this cooling agent. Critical flow velocity and mass fraction of ice values were determined thereby; they correspond to a change in character of an ice slurry flow from a laminar to turbulent motion. Experimental results were compared to the analytical results, based on the Hedström and Tomita algorithms (the laminar and turbulent flow, respectively). The comparison showed a very good agreement between these data.