Hostname: page-component-78c5997874-8bhkd Total loading time: 0 Render date: 2024-11-16T06:01:47.952Z Has data issue: false hasContentIssue false
  • English
  • Français

A Concise Method for Analyzing the Ice-Melting Performance of a Heated Disc

Published online by Cambridge University Press:  30 January 2017

J.F. Lea
J.F. Lea*
Affiliation:
Science and Technology Division, North American Exploration and Production Group, Sun Oil Company, 503 North Central Expressway, P.O. Box 936, Richardson, Texas 75080, U.S.A.
Rights & Permissions [Opens in a new window]

Abstract

Shreve (1962) evaluated the performance of solid-nose hotpoints of various geometries. A more concise method for analysis of the special case of a flat disc melting into ice is presented here. In contrast to the numerical solution presented by Shreve, this method results in a closed analytical solution, based on constant properties, which is corrected using a property ratio scheme. This technique allows the substitution of various fluid properties into the solutions as might be required when considering contaminated or saline ice or other liquids. Shreve's method would require some additional numerical work to evaluate various properties because the temperature dependence of the viscosity of water is integrated into the results.

Results obtained for the flat disc using this method are shown to compare favorably with Shreve's, indicating that this type of analysis may be applicable to other melt-unit geometries.

Résumé

Résumé

Dans l'article de Shreve (1962), le fonctionnement de "têtes chauffantes" de diverses géométries était estimé. On présente ici une méthode plus rapide pour l'analyse du cas d'un disque plat fondant la glace. Contrastant avec la solution numérique présentée par Shreve, cette méthode aboutit à une solution analytique de forme simple, basée sur des propriétés constantes, qui est corrigée en utilisant des variables réduites. Cette technique permet la substitution dans les équations des propriétés de différents fluides comme ce pourrait être nécessaire en considérant une glace polluée ou saline ou d'autres liquides. La méthode de Shreve demanderait des calculs numériques supplémentaires pour prendre en compte des propriétés variées parce que ses résultats intégrent la loi de variation de la viscosité de l'eau en fonction de la température.

On montre que les résultais obtenus par la méthode rapide pour un disque plat supportent avantageusement la comparaison avec ceux de Shreve, ce qui indique que ce type d'analyse, de même que celui de Shreve, est applicable à d'autres geometries de tête fondante.

Zusammenfassung

Zusammenfassung

In dem Beitrag von Shreve (1962) wird aus der Betrachtung des fliessenden Wasserfilmes zwischen der warmen Oberfläche und dem Eis auf die Leistung von "solid-nose hotpoints" verschiedener Gestalt geschlossen. Für die Analyse des Falles einer flachen, in das Eis einschmelzenden Scheibe, wird hier ein kürzeres Verfahren entwickelt und vorgestellt. Im Gegensatz zur der von Shreve angewandten numerischen Lösung liefert dieses Verfahren eine geschlossene analytische Lösung; es beruhr auf der Annahme konstanter Eigenschaften, wird jedoch mit Hilfe eines Schemas der Eigenschaftsverhältnisse verbessert. Dieses Vorgehen ermöglicht auch die Substitution verschiedener Flüssigkeitseigenschaften in die Gleichungen, wie sie bei der Betrachtung von verunreinigtem oder salzhaltigem Eis oder anderen Flüssigkeiten notwendig werden könnte. In Shreves Methode würde die Berücksichtigung verschiedener Eigenschaften einigen rechnerischen Mehraufwand erfordern, weil die Temperaturabhängigkeit der Viskosität des Wassers ein Bestandteil der Ergebnisse ist.

Die mit dem kürzeren Verfahren für die flache Scheibe erhaltenen Ergebnisse schneiden im Vergleich zu denen von Shreve günstig ab; alle Anzeichen sprechen dafür, dass diese Art der Analyst, wie die von Shreve, auch auTanders geformte Schmelzkörper anwendbar ist.

Type
Short Notes
Information
Journal of Glaciology , Volume 17 , Issue 76 , 1976 , pp. 355 - 358
Copyright
Copyright © International Glaciological Society 1976

Method of Analysis

A healed disc is considered to be melting into ice at the velocity V. The temperature of the escaping water drops from the disc temperature to the fusion temperature of the ice across the uniform distance δ.

In the following transformations, the radial velocity u, and the axial velocity v are non-dimensionalized with V. The axial distance z is related to the distance between the disc and ice δ, and the radial coordinate r is related to the disc radius, R. The variable p indicates the pressure distribution in the melt water.

A prime (') indicates a dimensional quantity. The subscripts w and i refer to "water" and "ice".

u = u’/V, v = v’/V,

z = z’/δ, r = r’/R,

p = p’/pwV2.

12A 355

The assumptions used to reduce the equations of motion are: (1) inertial Ibices are much less than viscous forces; (2) δ ≪ R; (3) steady state; (4) no body forces; (5) constant properties. The resultant equation for the r direction is:

(2)

where

pw is the water density and μ is the water viscosity at a "mixed mean" temperature defined in Equation (14).

To solve for the velocity and pressure distribution, the following boundary conditions must be satisfied:

(3)

Solving Equation (a) with the boundary conditions and using mass balance of the fluid gives the following expressions for the pressure and radial-velocity distributions:

(4)

(5)

The gap thickness δ can be found by integrating Equation (4) over the disc area:

(6)

where Pi indicates the effective average pressure acting on the disc.

An integral energy equation for the fluid layer is then derived. Conduction in the r direction is neglected. Equating the conduction and convection terms yields the energy equation

(7)

where T = T'/Tm and Tm is the fusion temperature; the Peclet number is c is the heat capacity of the water, and k is the conductivity of the water. The boundary conditions for temperature are:

(8)

where A is the heat of fusion.

The solution of Equations (7) and (8) is:

(9)

where

(10)

and a is a constant of integration.

If efficiency, η is defined as the quantity of heat that is conducted into the ice directly below the disc over the total heat leaving the disc, the following equation results:

(11)

Evaluation of η as (Pe) ⟶ o indicates a = o which makes T independent of r. This is in agreement with Shreve's work (1962), in which the isothermal surface assumption leads to a constant δ for a heated disc.

Comparison of Results

The temperature of the disc in Shreve's study was shown to be a function of the following dimension-less "performance number" for a disc:

(12)

where Q is the total input of heat, W is the weight on the "hotpoint" or disc and μ0 is the water viscosity at o°C.

In order to conveniently relate the performance number N to the Peclet number, a dimensionless number (Pe)* is defined as

(13)

where μ and (Pe) are evaluated at the "mixed mean" temperature [Inline Equation], between the disc and the ice.

(14)

The factor pi/pw is included in (Pe)* because Equation (9) appears in terms of (Pe)-pi/pw. The term (μo/μ)−014 is a farm of the "property ratio" scheme described by Reference KaysKays (1966) which is commonly used to correct constant-property analytic solutions for effects of variable viscosity. The exponent, — 0.14, which has been used successfully with pipe-flow solutions for laminar flow of moderate to high Prandtl number liquids (Reference KaysKays, 1966), was found, after trial and error, to correct this solution to near the values obtained by Shreve using fresh-water viscosity data.

Shreve's performance number N, and (Pe)* can now be related as

(15)

where E is efficiency as defined by Shreve: the "ratio of the cross-sectional area of the 'hotpoint' to the cross-sectional area of the hole made in the ice".

We have described efficiency as the ratio of the quantity of heat conducted into the ice directly below the disc to the total heat leaving the disc. Since the heat from the disc which does not get conducted to the ice directly below serves instead to enlarge the cross-sectional area of the hole created, it can be seen that the two definitions are mathematically equivalent.

The results of this study and that of Shreve are compared in Figure I. The term τ0 denned by Shreve is disc temperature multiplied by c/λ (heat capacity of water/heat of fusion). The discrepancy between the numerical solution using the more exact method of Shreve, and the constant-property method giving an analytical solution is seen to be very small. Also, the efficiencies η and E agree very well at low values of Nand (Pe)*, and differ by only a few percent at higher values ofN. The results indicate that this method is satisfactory for the analysis of a flat disc and should be considered for its possible application to other geometries.

Fig. 1. Efficiency and surface temperature related to performance number N, and (Pe)*.

Footnotes

page 355 note *

Present address: Department of Engineering Science, University of Arkansas, Fayetteville, Arkansas 72701, U.S.A. 12A

References

Kays, W. M. 1966. Convective heat and mass transfer. New York, McGraw-Hill Book Co., Inc. Google Scholar
Shreve, R. L. 1962. Theory of performance of isothermal solid-nose hotpoints boring in temperate ice. Journal of Glaciology, Vol. 4., No. 32, p. 151–60.CrossRefGoogle Scholar
Figure 0

Fig. 1. Efficiency and surface temperature related to performance number N, and (Pe)*.