広島工業大学紀要研究編
第 44 巻（2010）pp.25-30

論文

空気中に置かれた垂直氷層の自然対流による融解
橋詰 健一 *・小山 哲也 **・岡田 匡司 ***・馬谷 尭 ****
古川 亨 *****・村地 建太 ******・田邉 芳崇 *******（平成21年９月７日受理）

Melting of Ice Layer Placed Vertically in the Air
Under Condition of Natural Convection
Kenichi HASHIZUME，Tetsuya KOYAMA，Masashi OKADA，Takashi UMATANI，Tooru FURUKAWA，Kenta MURAJI and Yoshitaka TANABE(Received Sep. 7, 2009)

Abstract
Transient melting of ice layer placed vertically in the air under condition of natural convection was analyzed and compared with experimental results. First，the heat transfer coefficienton the falling film surface was investigated experimentally，where the natural convection andthe vapor condensation of the humid air occur simultaneously. In the data analysis，the Lewis
relation was used between the heat transfer coefficient and the mass transfer coefficient.
Then，applying the obtained heat transfer coefficient on the falling film surface，the transient
melting of ice was calculated. Lastly，a melting experiment was conducted using verticallyplaced ice layers and ice cylinders. The calculated results agreed well with experimental results.
Key Words: ice melting，natural convection，humid air，heat transfer coefficient，condensation，falling film

１．緒層の融解を扱った。また，Tien ら⑷ や Contreras ら⑸ は飽言

和水蒸気中での，Epstein ら⑹や Taghavi-Tafreshi ら⑺，⑻は水の凍結や氷の融解は自然界あるいは冷凍装置などにお非溶解性蒸気中での垂直氷層の融解を扱った。しかしながいて身近な現象であり，伝熱工学的にも興味のある問題で

ら空気中に置かれた氷の自然対流による融解についての研あることから，これまでに多くの研究が行われてきた。特究は見当たらない。空気中における氷の融解は夏季に涼感⑴に水の凍結については Stefan による北極海での海水凍結を得るために設置する氷柱や冷凍装置での氷層の自然放置の研究以来，水道管の破裂防止や氷蓄熱など工学的応用に時にみられる現象であるが，この現象では氷を解かす加熱⑵も関連して多くの研究結果が発表されている 。一方，氷源が融解液膜への空気の自然対流熱伝達と空気中水蒸気のの融解については水中や液中における融解を除けば研究例融解液膜への凝縮のふたつになり，従来研究とは大きく異⑶なっている。そこで本研究では融解液膜と周囲空気との間は多くはない。Roberts は高温空気流が衝突する水平氷
＊
＊＊
＊＊＊
＊＊＊＊

広島工業大学工学部知能機械工学科＊＊＊＊＊

広島工業大学工学部電気デジタルシステム工学科＊＊＊＊＊＊㈱サンワテクノス＊＊＊＊＊＊＊

三機工業㈱― 25 ―

新ダイワ工業㈱㈱シンコー
ツネイシホールディングス㈱橋詰健一・小山哲也・岡田匡司・馬谷尭・古川亨・村地建太・田邉芳崇の空気中水蒸気の凝縮を含む熱伝達率を測定し，その結果を用いて垂直氷層および垂直氷柱の融解過程を解析して実験結果と比較する。

記号

cp ：定圧比熱g：重力加速度

Pr ：プラントル数
q：熱流束

T：温度z

：氷層，氷柱，断熱円柱の上端からの距離

Fig.1 Analytical model for ice melting

αD ：物質伝達率αFF ：流下液膜熱伝達率⑴αH ：顕熱熱伝達率
αNC ：自然対流熱伝達率⑵αs ：みかけの熱伝達率，式⑻β：空気の体積膨張率⑶Γ：単位幅当りの液膜の質量流量

Δ：氷層の厚さ

γ：水蒸気濃度λ：熱伝導率ここで αs はみかけの熱伝達率，Γ は単位幅当りの融解液

μ：粘性係数膜の質量流量，Δhm は融解潜熱，ρs は氷の密度，τ は時間

ν：動粘性係数である。

ρ：密度τ：時間⑷一方，液膜表面での熱流束 q は⑸

添え字と書ける。ここで qH は顕熱による伝熱量で，顕熱熱伝達

a：空気

f：液膜

0：氷の融解温度（０℃）

率を αH とすれば⑹である。qD は空気中水蒸気の凝縮による伝熱量で，物質２．氷層融解過程の解析

伝達率を αD とすれば解析では図１に示すように垂直氷層の上端から重力方向⑺に z 軸をとり，位置 z における氷層の厚さを Δ として，次の仮定をおく。

である。ここで γa と γ0 は周囲空気中および液膜表面にお⑴ 氷層の温度は一様で，融解温度（T0 ＝ 0℃ ）とする

ける水蒸気濃度，Δhv は蒸発潜熱である。これらより式⑴⑵におけるみかけの熱伝達率 αs は融解液膜は非常に薄く，その温度を T0 に近似する⑶ 位置 z において空気から伝達される熱量はすべて氷⑻の融解に使われる⑷ 空気中水蒸気の凝縮による融解液膜流量の増加は省と書ける。さらにここで，顕熱熱伝達率 αH と物質伝達率αD との関係に対して湿り空気の場合によくあてはまる⑼と略する
さらに，解析を簡単化するために微小区間 dz においてはいわれているルイスの関係⑷融解面は垂直とする。この取扱いは Tien ら ，Contreras⑼

ら⑸，Epstein ら⑹，Taghavi-Tafreshi ら⑺も行っている。
以上の仮定より，融解液膜表面での周囲空気からの熱流を用いれば式⑻は束を q，周囲空気の温度を Ta とすれば，氷層厚さ Δ の時間的な変化は次のように表わされる。
― 26 ―

空気中に置かれた垂直氷層の自然対流による融解等な液膜⑺を形成して流下する。流下する水の流量はオーバーフロー管の位置（液面高さ），およびキャピラリーチューブの内径と長さにより設定し，断熱円柱の下端でメスシリンダーとストップウオッチにより測定した。流下する液膜の温度（水温）は線径 0.2 ㎜のＴ型熱電対を用いて上端から 25 ㎜間隔で測定した。周囲空気については乾湿計を用いて乾球温度 Ta と湿球温度 Taw を測定した。実験は，液膜の温度が周囲空気温度よりも低く，液膜が加熱される場合（Tf ＜ Ta）と，この逆の場合（Tf ＞ Ta）について行った。測定された液膜温度から顕熱熱伝達率は以下の手順で算出した。
図３において微小区間 dz での液膜温度の変化 dTf は⑾として表わされる。ここで熱流束 q は式⑴，⑽で表わされ，

1: Thermostatic bath 2: Overflow pipe 3: Reservoir
4: Capillary tube 5: Sintered metal pipe 6: Insulating column
7: Falling film

融解液膜温度 T0 と水蒸気濃度 γ0 はそれぞれ液膜温度に対応する Tf と γf に置き換える。すると式⑾はFig.2 Experimental apparatus⑿となる。さらに，液膜表面での顕熱熱伝達率 αH は周囲空気の自然対流熱伝達と，液膜が流下することによる一種の強制対流熱伝達との共存対流であると考える。そして自然対 流 熱 伝 達 率 を αNC， 流 下 液 膜 熱 伝 達 率 を αFF と し て
Churchill の方法⑽⒀を用い，αNC に対しては同じく Churchill の整理式⑾Fig.3 Analytical model for falling film temperature⒁⑽

を用いる。式⒁において Ra はレイレイ数⒂となる。以上より，任意の状態に対して顕熱熱伝達率 αHが予測できれば氷層厚さ Δ の時間的な変化は式⑽を式⑷に用いることによって数値的に求めることができる。

である。以上の各式より顕熱熱伝達率 αH を求めるには適
当な流下液膜熱伝達率 αFF と液膜の入口温度 Tf0 を仮定し，３．融解液膜と周囲空気との間の熱伝達円柱上端から微小区間 dz について式⒀，⒁と式⑿を用いて液膜温度を順次計算し，測定値に最も良く一致する値を微小流量で流下する液膜と周囲空気との間の熱伝達率を最小二乗法により見出せばよい。

求めるために図２に示す装置を用いて実験を行った。この

図４に液膜温度の測定値（プロット点）と計算値（実線）

装置は恒温水循環装置⑴からの水をオーバーフロー管⑵を持つ水槽⑶に入れ，キャピラリーチューブ⑷を介して焼を示す。図４上図は Tf ＜ Ta の場合で，液膜表面では空気結金属管⑸に滴下する。ここからの水は直径 24 ㎜ × 長さ

中水蒸気の凝縮が起こっている。図４下図は Tf ＞ Ta の場

300 ㎜の断熱円柱（ベークライト製丸棒）⑹の外表面に均合で，液膜表面から水の蒸発が起こっている。実験は液膜― 27 ―

橋詰健一・小山哲也・岡田匡司・馬谷尭・古川亨・村地建太・田邉芳崇

Fig.4 Experimental results of film temperature

Fig.6 Melting of vertical ice layerの関係が知られている⑿が，本実験（流下液膜と周囲空気との間の熱伝達，Re ＝ 0.5 ～ 20）では液膜流量の影響がかなり小さくなった。

Fig.5 Falling film heat transfer coefficient

４．氷の融解過程の実験と解析結果との比較
4.1の体積流量 V ＝ 36 ～ 1500 ㎤ /h の範囲で行ったが，液膜垂直氷層

流量と温度変化の傾向が対応していないのは周囲空気の温供試氷層は冷却銅板を用いて作成した。冷却銅板は長さ度（乾球温度）と湿度（湿球温度）が実験ごとに異なって150 ㎜，幅 50 ㎜，厚さ 10 ㎜で，背面にブラインを通流さいるためである。図４の上下図ともに液膜の入口温度 Tf0

せる冷却流路を有し，表面には氷層滑落防止のために直径を含めて測定値と計算値がよく一致している。

１㎜，高さ５㎜の竹製ピンを６箇所設けてある。この冷却図５は流下液膜熱伝達率 αFF の計算値で，白抜点で示す銅板を製氷型枠に嵌め込んで２℃の煮沸水水槽に浸漬し，Tf ＜ Ta の場合と黒塗点で示す Tf ＞ Ta の場合とで違いが見－10℃のブラインを冷却流路に通流させて煮沸水を約 15られない。これらの点は㎜の厚さに結氷させた。そして氷層が張り付いた冷却銅板を風除けを有する測定箱の壁面に垂直に設置し，氷層厚さの変化をマイクロメーターで測定した。測定箱の下面は大⒃気開放であり，内部空気は測定時間中，実験室内の空気状態にした。実験中，冷却銅板の冷却流路には０℃のブライで表わされる。流下液膜と固体表面との間の熱伝達率につンを流通させて氷の過冷却を除去すると共に氷層背面からいては膜レイノルズ数 Re ＝ 320 ～ 20000 においての入熱を遮断した。
図６に実験結果（白抜点）と解析値（実線）との比較を⒄示す。実験における周囲空気の温度は Ta ＝ 25.5 ～ 27.5℃，― 28 ―

空気中に置かれた垂直氷層の自然対流による融解

相対湿度はφ＝ 67 ～ 71％であった。
解析（数値計算）は次のようにして実行した。まず，初期値として τ ＝ 0 での氷層厚さ Δτ ＝ 0 を与える。上述した製氷方法では氷層の表面を平坦にすることは困難で，図６における τ ＝０min のような形状になったが，これを多項式近似する。次に，氷層上端においては Γz ＝ 0 ＝ 0 として，微小区間 dz ごとに順次，式⑽，⒀，⒁，⒂によりみかけの熱伝達率 αs を求め，式⑴による熱流束 q から各区間における Γ の値を式⑵より求め，時間 dτ 後の氷層厚さ Δ を式⑷より算出する。空間と時間の分割は dz ＝２㎜，dτ ＝15s とした。
図６には初期氷層形状と周囲空気温度が異なる場合を示しており，下図の方が周囲空気温度と湿度が高いために融解速度は大きくなっている。いずれの場合も融解の初期においては解析と実験とがよく一致している。しかし，上端の氷が解けてくると，露出した冷却銅板の上部が空気によって加熱されるために融解量は解析よりもかなり大きくなった（上図では 80 分以降，下図では 60 分以降）。この理由のため本実験の方法では氷層が完全に消滅するまでの過程を解析と比較することはできなかった。
4.2

垂直氷柱

供試氷柱はアクリル製円筒の中に煮沸水を注入し，これを－10℃の冷凍庫に保管することで作成した。円柱の直径は 24 ㎜と 48 ㎜，高さは 100 ㎜と 150 ㎜の合計４種類である。この供試氷柱を第 4.1 節で説明した測定箱の中に設置し，氷柱の中央部における直径の時間的な変化をノギスを用いて測定した。実験における周囲空気の温度は Ta ＝21.5 ～ 27.5℃，相対湿度はφ＝ 64 ～ 75％であった。
解析（数値計算）は垂直氷層の場合と全く同じで，使用する式も氷層の厚さ Δ を氷柱の半径 R に置き換えれば同一である。空間と時間の分割も同じにした。
図７に実験結果（白抜点）と解析値（実線）との比較を示す。垂直氷柱の場合には冷凍庫から取り出して過冷却を除去しないまま実験を開始したので，τ ＝０から 30 分後までの間は融解速度が遅いため，解析との比較では 30 分

Fig.7 Melting of vertical ice cylinder

後の値を基準としている。氷柱の実験では中央部における直径の時間的な変化しか測定していないが，実験と解析は非常に良好な一致をみせている。

５．結計算の結果と実験結果とはある程度良く一致し，解析方法と流下液膜伝熱実験の妥当性は示された。一方，氷の融解

言実験についてはいくつかの改善すべき点が見出されたの空気中に置かれた垂直氷層の自然対流による融解過程をで，現在それらを改善した氷柱の融解実験を進めている。

解析し，実験結果と比較した。解析に際して必要になる融１点は形状寸法の測定法である。本報ではマイクロメータ

解液膜表面と周囲空気との間の熱伝達率は，これを流下液あるいはノギスを用いて測定したが，この測定法では測定膜の熱伝達と自然対流熱伝達との共存対流とみなして，流箱の中に手を入れなければならず，測定にも時間がかかる

下液膜の実験から熱伝達率を求めた。この値を用いた数値ので測定箱の中の空気を乱してしまう。これが，図６や図― 29 ―

橋詰健一・小山哲也・岡田匡司・馬谷尭・古川

７で測定値の方が解析値よりも融解速度が早い原因になっ⒀

亨・村地建太・田邉芳崇

田邉芳崇，平成 18 年度広島工業大学工学部知能機械たものと推測される。この点の改善策として光学的形状測工学科卒業論文

定装置を導入した。もう１点は供試氷柱の過冷却除去で，⒁

古川亨・村地建太，平成 19 年度広島工業大学工学部

この目的のために，冷凍庫から取り出した供試氷柱を先に知能機械工学科卒業論文

０℃に設定したブラインに浸漬し，氷柱中心部の温度を測⒂

馬谷尭・岡田匡司，平成 20 年度広島工業大学工学部

定して過冷却が除去されたことを確認してから測定箱に設

知能機械工学科卒業論文置している。その他にも，氷の融解速度には周囲空気の温付録：解析に使用した物性値と状態量度とともに湿度も大きな影響を及ぼすことがわかったので，測定箱内に湿度制御も導入した。これらの改善を施し

以下は標準大気圧（Pa ＝ P0 ＝ 101.3kPa）における値て現在進めている実験では氷柱が完全に融解消滅するまで⑴の過程を解析と比較できるものと期待している。

水の物性値，水温：・定圧比熱

なお，本報は広島工業大学工学部知能機械工学科の卒業
研究⒀～⒂として取り組んだ３年間の成果をまとめたものである。

・蒸発潜熱

文
⑴

献

Stefan, J., Annalen der Physik und Chemie, vol.42,⑵

周囲空気（湿り空気）の物性値と状態量

pp.269-286,（1891）
⑵

乾球温度：

福迫尚一郎・平田哲夫・稲葉英男・菅原征洋，機械の研 究， 第 41 巻 第 １ 号 ～ 第 43 巻 第 12 号，（1989 ～，湿球温度：・飽和水蒸気圧，1991）⑶

Roberts, L., Journal of Fluid Mechanics, vol.4, pp.505528,（1958）⑷

・飽和絶対湿度

Tien, C. and Yen, Y. C., Chemical Engineering
Progress Symposium Series, vol.67, no.113, pp.1-9,（1971）

・湿り空気の絶対湿度⑸ Contreras, W. and Thorsen, R. S., Transactions of theASME, vol.97, pp.570-575,（1975）⑹

Epstein, M. and Cho, D. H., Transactions of the・湿り空気の比体積，

ASME, vol.98, pp.108-113,（1976）⑺

Taghavi-Tafreshi,K. and Dhir, V. K., Transactions ofthe ASME, vol.104, pp.24-33,（1982）⑻

・周囲空気（湿り空気）の密度

Taghavi-Tafreshi, K. and Dhir, V. K., Transactions ofthe ASME, vol.105, pp.235-240,（1982）⑼

内田秀雄：湿り空気と冷却塔，裳華房，pp.76-83，（1965）⑽

Churchill, S. W., Handbook of Heat Exchanger Design,

・周囲空気（湿り空気）の定圧比熱

Hewitt, G. F., ed., Begell House, New York, p.2.5.9-1,（1992）⑾

・水蒸気濃度Churchill, S. W., Handbook of Heat Exchanger：液膜表面

Design, Hewitt, G. F., ed., Begell House, New York,：周囲空気

p.2.5.7-3,（1992）⑿⑶

Saibabu, J., Sarma, P. K., Dharma,R.V. and Raju, G. J.
V. J., Proceedings of the 8th International HeatTransfer Conference, vol.4, pp.1963-1969,（1986）― 30 ―

周囲空気の相対湿度φ