〈支部の話題〉

〈車載〉

山田悦郎(秋田大学)

1.はじめに
氷の融解や水の凍結はNeumannまたはStefan問題として，およそ1800年代の中頃から関心がもたれてきて，今もなお研究の対象とされ息の長い問題と言える。この問題の特徴は解析的には固液界面の位置が時間経過につれて変化するため解析が難しいことにある.そのため,時間と融解／凍結量に関する厳密解は長い研究歴史にもかかわらずNeumannとStefanの解[1]のみである．これは誤差関数を巧みに利用したために解析の成功に至ったと思われる。この厳密解は境界温度が一定など条件が限られているため応用範囲が狭く一般的でない.このため,多くの近似的な取り扱いや数値解析が伝熱関係の論文集に登場するようになってきた.凍結量の予測にかかわせない融解／凍結速度は次に示す固液界面に出入りする熱流束で容易に決まることは周知の通りである.
ただし，ね,は融解の潜熱, mは単位時間,単位面積あたりの融解質量(融解速度)，えは然伝導率，Tは温度，xは座標，添え字Lは液体，Sは固体を示す。したがって，何らかの方法で融解／凍結系の温度場を解析し、ある時刻での界面の温度勾配を算定することによってその時刻での融解速度を推定することができるのである．仮に固液の温度がいたるところで一様になり，界面の温度勾配が零になれば融解も凍結もしないことになることは式(１)から明らかである.数値解析の発展によって今や一次元問題のみならず複雑な幾何学的な形状を持つ２次元問題や３次元問題でも解析が可能になってきた.しかしながら，海水や各種水溶液内の純水の融解となると，純水中の融解と比較して取り扱いが難しくなり，融解現象としてもまだ未知な部分が
菅原征洋(秋田大学)
Masahiro SUGAWARA (AKITA University)
存在するように思える．たとえば，0℃の水溶液にこれと同じ温度の氷を浸す場合を考えてみる．少なくとも融解開始時に両者の温度差がない状態とみて良いであろう.このとき式( 1 )から温度勾配が存在しないので融解速度m は零になって氷は融解しないことになる．しかしながら，-5℃で20wt%の塩化カルシウム水溶液にこれと同じ温度の純氷を浸すと，融脾面の温度が急激に降下して氷は自発的に融解していくことは実験から認められている.この現象は何も目新しいことではなく,冬期間道路上の積雪に薬剤をまいて融雪や凍結防止を行っているが，これと基本的には同じ現象である．熱力学的には融解系に化学ポテンシャルが存在するため，氷は一種の反応挙動ともみられる融解が生じることになる.しかしながら,具体的にある時間内にどれだけ融解するのか，つまり融解速度となると,式( 1 )のほかに次に示す濃度拡散の式を用意しなければならない.
m
ただし，Ｃは融解面の平衡濃度， ρは密度，Dは物質拡散係数である.この式は，融解面をsemi-permeable planeと見なした関係式である.つまり,溶質分子が氷内部に浸透しないとするのである．このように熱バランス式（1）と質量バランス式（2）を連立させて融解速度m と数力学的に平衡関係にある融解面の温度T 。と融解面の濃度C 。を求めるのである.著者らの報告[2,3]に述べているように，これまで,溶液中の融解(溶解)問題は海水中の氷山などの融解問題について最初海洋学者らに関心が持たれてきたようである.野外実験,研究室実験,解析的な取り扱いなど多くの報告がなされている.しかしながら，融解速度に対する定量的な議論は充分でなく不明な点が多い。このような問題は前述のNeu-mann／Stefan問題よりは研究歴史が浅く発展途上のように思えるが、しかし、今から100年ほど前に熱力学の第３法則を提唱し，ノーベル賞の受賞者でもあるNernstもこれに関心を示し，溶液中の固体の溶解反応でいわゆるNernstのモデルを提案した[4]. これは溶融面に0.03mmの大変薄い濃度拡散層が存在して，ここの濃度勾配によって拡散する量(Fickの法則)に見合うだけ溶融するというモデルである．拡散層の外側の溶液の濃度は一様濃度と仮定するという，いわゆる折れ線濃度分布モデルである。この物理的に理解しがたい大変薄い拡散層厚さδ（＝0.03mm）の存在が議論の対象になると思うが，しかし，多くの溶解実験結果を説明できるとしている.ただし,このδの値は数多くの溶解実験結果から反応速度定数を求めてδの値を算定したのである.著者らは氷の濃度・温度複合融解に関心を持ち,研究を続けてきているが,これまで,報告した主な結果を要約して紹介する[2,5].
2.垂直氷板の融解
図1は氷板と塩化カルシウムの初期温度がそれぞれTs＝T u = - 5℃,水溶液部(水板)の高さH = 100mm，幅L = 50mm，氷板の厚さδ = 10mmの場合の破解系の流れの可視化写真(a)と数値解析による流れのベクトル線図(b),等温度分布線(c)および等濃度分布線図(d)である。ただし，数値解析は連続運動量,エネルギー,質量保存の格式を連立させて，前述の融解面の境界条件である式（１）と（2）を満足するように繰り返し計算を行うのである.これをみると,最初液層内には対流が存在せず静止しているが，融解に伴う潜熱吸収のため融解面近傍の溶液の温度が降下して密度増加となり負の浮力(下方向,重力方向)が作用する。これに反し融解液の拡散のため溶液は水成分に富み密度の低下をもたらして浮力(上方向)が作用する．しかしながら,濃度体膨張保数βが温度体膨張保数β ,より格段と大きいことによる濃度差浮力の優勢のため融解面近傍で濃度境界層が形成され融解液の上昇流がみられる(b-1). 融解時間の経過に伴って温度・濃度共に低い融液が融脾面を上昇して液層上部で水平方向に広がり，可視化観察に見られるような安定した濃度分布の存在する濃度成層領域(Stratified region)が発達していく様子がわかる(d- 1).この領域の速度は大変小さく,ところによってはほとんど融解液が停止しているように見える.一方，濃度成層領域の下に反時計方向周りの強い渦が存在する。これは等濃度分布線図(d-1)との対比でもわかるように,この領域の濃度は,ほぼ初期濃度に近く,融解液と混じりにくいので,非混合領域（Uncontaminated region）とも呼ばれている．したがって，この領域では濃度差浮力が存在しないので,水溶液の温度低下に伴う負の温度差浮力(重力方向)のみが存在するので,低温の溶液が液層底部へ広がっていくのである。したがって,液層上部と底部で温度の低い領域が存在するのが特徴である．液層内の対流模様は渦の中心が融解面近傍に現れ，他の自然対流では見られない特異な形をなしている.しかも,液層底面をはってきた流れが融解面の対面から引き返して,また上昇している.氷板の融解を伴う複雑な温度・濃度複合対流の時間的変化に対する数値解析結果は可視化観察結果と類似していることが認められる.図２は図１と同じ計算条件での垂直方向速度vの水平（x）方向の分布を示す。この図から，前記の液層内の流れがいっそう明確になり,融解面( x =0)近傍の濃度境界居内の融解液の上昇流が明確に認められ，この厚さは極めて薄いのが特徴であvi図３は液層内の温度Tと水濃度(1-C)の融解面( x = 0)付近の拡大図である．破線で示されている水濃度勾配は極めて急であり，液層の底，中，上にいくにしたがって勾配がゆるやかになっている.これは融解液が濃度境界層を上昇するにしたがって次第に水成分に富んでいくからである．この濃度勾配の変化は融解速度の変化に対応する。水濃度分布は濃度境界層内では急激な濃度勾配をもつが，これの外側では水平に近い一様濃度分布となっている.このような折れ線濃度分布は,前述の数多くの実験結果を基にして仮定したNernst (1904) [ 4 ]のモデルと形の上では大変良く似ているのは興味深い。　しかし, Nernstの拡散層厚さは0.03mmであるのに対して，著者らの結果は図より0.3mmで約10倍厚くなっている．温度勾配は水濃度勾配と比較して大変小さいが，これは本融解系のルイス数Leが約400と大変大きな値をもつことに起因する.図４は氷板内部中央部の温度Ｔ ,と融解面と対面する断熱壁の上部温度T 2と下部温度T 3および平均融解質量Mの時間変化に対する数値解と実験値を比較したものである．融解初期を除いて融解面温度T ,にほぼ相当する氷板温度T ,は融解開始から急激に降下し,その後,おだやかに単調に降下していく傾向を示している.一方,融解面から離れた温度T 2とT 3はT 1よりかなり遅れて液層内の温度・濃度複合対流によって降下していくのが認められる.本融解問題の複雑さからみて,本数値解析は実験の傾向を比較的に良く予測できているように思え,氷板は濃度勾配(濃度拡散)に駆動されて融解し，しかも温度・濃度複合対流に支配されて融解するとした融解モデルはほぼ妥当であると言えよう.伝熱  2000年11月
3.水平氷円柱の融解
図5は融解時間t=5 minのときの液層流れの(a)可視化観察写真と数値解析で得られた( b )速度ベクトル, (c)等温度分布線, (d)等濃度分布線を示す.図1の垂直氷板の場合と異なって,融解液は氷円柱上面から離れ液層上部を上昇する。周囲の液はこの上昇流に巻き込まれて，液層上部で渦状の濃度差対流が形成される.一方,液層下部では渦状Jour. HTSJ, Vol.39, No.159の温度差支配対流が現れる。上方の濃度差支配対流と下方の温度差支配対流との境界付近に流速の大変遅く，ほとんど停滞している領域が認められる。この停滞層の濃度は融解水と混じり合わない初期濃度に近く，しかも温度が融解系で最も高いので，密度的にちょうど上下対流の中間にあるので濃度差支配対流と温度差支配対流に浸食されない安定した領域と思われる。時間の経過につれて流れが弱くなるため，この停滞層の厚さは次第に厚くなっていく.図6は融解面における角度と局所ヌセルト数Nu 』の教値解析結果を時間をパラメタにとって示した．θ = 0 degは融解面の頂点であり，θ =180degは最下部の位置である。融解面上部で局所ヌセルト教は比較的に大きな値を示している．しかしながら，矩刃状に乱れて変動していることが明確に認められる。これは，図５（ａ）の可視化写真からも推測されるように融解面上部で流れが乱れて,しかも小さな渦が沢山現れている可能性があると思われる.このことは,上部融解面で絞肌のような細かい四凸面が観察されたことからも推定できる.融解面のθが大きくなるにつれてNu eは次第に小さくなって極小値を示し,また上昇している.この極小債を示す箇所はちょうど液層内の停滞層が存在している位置であると思われる。多少の変動が伴うもののNu aは次第に減少して最下部で最も小さな局所ヌセルト教となっている.図7は平均融解ヌセルト数とフーリエ数との関係であり，実験結果と計算予測結果と比較してある．破線で示したヌセルト数Nuは氷層からの熱流東を除いたものである(注:純水中の融解の場合の氷層内の熱移動は融解速度を減少させるが，水溶液では逆に促進する）一方，実線で示したヌセルト数Nuは氷の融解に費やされた全熱流束q (氷層と液層からの合計熱流束)で計算したものである.融解初期で両者の間(実線と破線)に差異が認められるが，時間の経過につれて両者の差異が次第に減少していることが認められる．したがって，ある程度，時間が経過したら融解に費やされる熱流束qでヌセルト教を算定しても許容されると思われる．しかしながら，融解初期では両者の差異が大きいので注意を要する．融解初期で大変大きなヌセルト数を示している。これは,氷と水溶液の間に初期温度差がない条件で融解が開始しているので,ヌセルト数を定義している融解面と周囲の温度差が極めて小さいにもかかわらず融解初期に融解面近傍に大きな温度勾配(濃度勾配)が現れて融納が自発的に急速に進行したことを示す。時間の経過につれて液層内の対流が発達して，準定常的な安定した値を示しており，実験結果と比較的に良く一致していることが認められる.
4.おわりに
水溶液内の氷の融解について著者らがこれまで研究を進めてきた主なものを取り上げた．純水中の氷の融解と異なって，熱拡散と比較してその速度が大変小さい濃度拡散も関与するので解析などの取り扱いが難しくなる反面，純水中の融解には見られない興味ある現象も見られた．この種の研究は地球温暖化による極水の融解にかかわる環境予測にも応用できる外，最近相変化を利用した蓄然装置などで凍結による圧力上昇を軽減するため水溶液を凍結させる試みもなされている[ 6 ].米ぼかりでなく，合金など新素材の開発にかかわる凝固・溶融問題にも関連するため,今後の研究発展が望まれる. [1] Carslaw, H. S. and Jaeger, J. C., Conduction of Heat in Solids , Oxford at the Clarendon Press (1959). [2]菅原征洋,藤田忠,日本機械学会論文集Ｂ編, 63-612(1997)2784. [3] Sugawara, M. and Irvine, T. F. Int. J. Heat Mass Transfer , 43-9(2000)1591. [4] Moelwyn-Hughes, E. A., The Kinetics of Reactions in Solution , 2nd ed., Oxford at the Clarendon Press (1947). [5]菅原征洋，沢木昭祥，田子真，藤田忠，日本機械学会論文集Ｂ編，65-636(1999)2785. [6]笹口健吾，石川昌稔，牟田健二，芳野清隆，馬場敬之,第37回日本伝熱シンポジウム講演論文集(神戸)  Vol. I   (1999)  9. Cold Heat / Heat Storage System using Ice ( Ice Melting Heat Transfer in a Horizontal Elliptical Tube Located in Water ) Koichi HIROSE (Iwate University) 1 18 [3] [4] 2000 1 3 1 3 15mm 2 [1] [2]
15
0.5mm 4 4.1 3 2دںT C, Boussinesq (5) (1) (2) (3) (4) 2 3 5 10 3 3 (a),(b),(c) S 110min 10 45min 15 S 4. 1. 1
– ω = ∇ 2 ψ 。 4. 1. 2
T w,ini SOR
Nu T n (t) D
5
(a)
(b) −ω = ∇ ² ψ
f [6] 10 4. 1. 3 (b) 10 Nu Nu 10 Nu Nu T =10 S (a) (b) (t=110min) (a) 17mm Ra (b) (a) (b) T =5 T =10 (a) S (2) Nu PCM (t=110min) (b) (a) S [1] 2 B 61-589(1995)3296. [2] ( 1996) . 36 (1999)353. B [4] 2 63-608(1997)1375. [5] ( 1978)   294 Saitoh,T., Trans. ASME, J. Heat Transfer, 100 4 3, [6] (1974). 5 [3] Control of Spray Characteristics and Their Applications Masatoshi DAIKOKU (Hachinohe Institute of Technology) 1 1980 [2] l [3] 2 M n /M 1 d 2 /d 1 =1.5 2.5 d32 d 2 /d 1 =2.25 2 2 1 [1] 3(a) 70%+ 30% l 2 (a) 2 2 OA 4 [4-6] 10 mL/s 91 109 SMD 170μ m C Lambert-Beer F = F o exp(− C s L) F Fo L (dp)max (dp)min 10μ m (C )o Time=0 min C 4μ m 10μ m /(C )o=1 1دں2 4μ m 5μ m 3(b) 5μ m 4μ m s 4 3 1 2 6 S [9] Larger Port    LP    4 C /(C )o=0.4 12 4 91μ m 4 Smaller Port   SP   8 2 12 2 4 6 LP [10] 3 6 No.1 No.2 SMD1 SMD。 SMD SMD mL m L SMD =
SP LP 70μ m 7 mL 7 g/s mLs m L 1 [11] 7 [12] 9 m A=6.0 g/s mL SP LP mL =mLs mL1 m A=6.0 g/s m L s  m L l mL m L s  m L 1 mLs LP 6.65 g/s mL 7 11.9 g/s mL SP 8 7 LP SP 8 [3]  Daikoku, M., et al., Proc. of 4 th Int. Conf. on Liquid Atomization and Spray Systems, (1988)465 9 [4] (1997)36 [5] 10 (1998)192 [6] 11 (1999)272 Atomization Sprays, A.H., [7] Lefebvre, and Hemisphere, (1989)155 [8] ） 52-476(1986)1784 (B [9] 3 (1994)254 [10]  Daikoku, M., et al., Proc. of 7 th Int. Conf. on Liquid Atomization and Spray Systems, (1997)318 [11] [2] 8-23(1999)139 8-24(1999)180 [12] Daikoku, M., et al., Proc. of 8 th Int. Conf. on Liquid Atomization and Spray Systems, (2000)1378 4 A Challenge to Extremely High Heat Flux Removal Technique using a Porous Medium Saburo TODA, Kazuhisa YUKI (Tohoku University) 1. [3,4] MW/m² EVAPORON 10MW/m² 2
[5,6,7] [9,10] [12][13]
[11] [14,15] (4) 50MW/m² [1,2]
(2) (3) EVAPORON (Evaporated Fluid Porous (4) Thermodevice)
Fig. 2 3mm 8mm 2 Fig. l 300mm (5) 10MW/m² l mm ± 0.75% TC4, TC5, 5mm 3 TC6 TC2 20mm TC1 2 TC3 0.5mm EVAPORON 3 3 1 120L/min 260W Fig.  1 EVAPORON Table  1 Fig. 2
4 #2 #1 #2, #3, #4 3mm Inlet velocity [m/s] #2 #1 #2, #3, #4 T out (1), (2) 3 2 4. 4.1 W 5mm Fig. 3 10mm 300A 30, 100, 25L/min H #1-4 #1 H 4.5, 7.0, 9.5, 12cm 80 #4 #2 10~70 60 Fig.  l w EVAPORON , a , a , exp(a , x , , ) q n (1) w kCu X w 4. 2 EVAPORON β H H Fig. 4 9 m hfg m cl H=7.0cm #1 w 5mm T in (2) Fig. 5 0.034m/s TC4-TC6 30 TC3 TC2 20mm TC1 Fig. 7 x=0, 0.02[m] H Fig. 5 #2 w 5mm H=4.5, 7.0, 9.5, 12cm 0.046m/s, 0.036m/s, 0.034m/s, 0.020m/s x=0 x=0.02 H=12cm 2.5cm Fig. 6 #1-4 w=2.5mm H=12, 9.5, 7, 4.5cm 5.76, 9.79MW/m² 3.01, 3.80, Fig. 7 H=4.5, 7, 9.5, 12cm x=0 x=0.02 H Fig. 9 #1-4 w=2.5, 5mm Fig.  7 H #1, 4, 3, 2 2.5mm 5mm 4. 3 EVAPORON 9 CHF  = C · u (3) c 10MW/m² #1 0.11MPa 0.14, 0.13, 0.14MPa #2, 3, 4 0.1MPa 10MW/m² Fig. 10 5mm Fig. 8 #2 w=2.5mm H=9.5cm 0.017m/s 3.4MW/m² 9CHF #1-4 40sec [3] 10 -4m/s EVAPORON Fig. 11
1.74× 10 8 w=2.5 c #4 w=5
5 EVAPORON 1) 2) [1] Boscary, J., Araki, M. et al., Fusion Technol., 35, 3, (1999) 289. [2] Araki, M, Sato, K. at el, Fusion Technol., 29 (1996) 519. [3] S.Toda, S.Ebara et al, Proceedings of the llth International Heat Transfer Conference, 4, (1998), 503. [4] S.Ebara, S.toda at el, Heat and Mass Transfer, 36, (2000), 273. 1995 61, 582 [5] 56, 529 1990 [6] A. P. De Wasch, G. F. Froment, Chem. Engng. Sci. [7] (1972) 567. [8] C. H. Sondergeld, D. L. Turcotte, J. Geophys. Res. 82 (1977) 2045. M. Chung, I. Catton, Nucl. Eng. Des. 128, (1991) 289. [9] [10] I. Catton, M. chung, Nucl. Eng. Des. 151, (1994) 185. [11]C. Doughty, K. Pruess, Int. J. heat Mass Transfer 31, (1988) 79.. [12]K. S. Udell, ASME J. Heat Transfer 105, (1983) 485. [13] C. Figus, Y. Le Bray, S. Bories and M. Prat, Proc, of the 11th International Heat Transfer Conference 4, (1998) 393. [14] J. Y. Chang, S. M. You, Int. J. Heat Mass Transfer 40, (1997) 4449. [15] J. Y. Chang, S. M. You. Int. J. Heat Mass Transfer 40, (1997) 4437. [16] H. C. Tien and K. Vafai, Int. H. Heat Mass Transfer 33, (1990) 1263. Joint Research 6 Interviewing a group leader of joint research 6.    Dr. Takaaki INADA, Senior Researcher of Mechanical Engineering Laboratory, AIST MITI Interviewer: Yasuo KAWAGUCHI(MEL) 250kW 6m 600mm 20 300mm l m 1/10 * 300 kW 1.8% [1] . kW , vol.65 (1999), , vol.17 (200), [3]国際伝熱会議報告2000年度国際伝熱会議アセンブリ委員会に出席してReport on the 2000 Meeting of Asssembly for International Heat Transfer Conferences 1998年8月に韓国の慶州で開かれた第11回国際伝数会議に続いて,第12回会議はフランスのGrenobleで2002年8月18 ～ 23日に開催される．韓国の会議では，日本の発表論文数は99編で，米国の84編，韓国の52編を抑え最大数であったそれも他国は,探択されても会場で発表(ポスターセッション形式による）されなかった論文の教がかなりの教に上ったり，論文の割り振り数に対して提出された論文数が少なかったりしたが、日本は採択された論文の全編が会場で発表された２国の内の一つであった．他の一国は論文割り振り数10編に対し採択数9編のオランダである.論文数の割り振りは会議の2年前に開かれるアセンブリ委員会で,実績に基づく主催国からの提案をもとに審議・決定されるが，そのようなGrenoble会議のためのアセンブリ委員会が，2000年9月10日午前10 ～ 12時，ちょうど今年度のEuropean Thermal Sciences Conferenceに合わせて同じ会場のKongress Haus Stadthalle Heidelbergの中で開催された.この委員会は,アセンブリメンバーの15か国を代表する2人ずつの委員で構成されており,国際伝熱会議の歴史的経緯から，これまでは各国とも，機械系と化学工学系を代表する委員を選出してきた.日本からは平田と京都大学猿野文丸教授両名が出席してきたが,今回は荻野教授が公用のため出席出来なかったので平田一人が出席した。因みに決議の際の投票数は1国あたり２票の権利を持つ主要な議題は，アセンブリメンバーの新規加入の承認，４年に１回開かれる国際伝熱会議の開催地や地域別の論文教の割り振りを決めることなどである.今回の会議は会長の前回会議の開催国である韓国のProf. S.T. Roがアセンブリ会長として議長を務めたが、出席者については正式の議事録が未だ到着していないので省略する.会長の挨拶のあと議事に入り，まず1998年8月28日(金)に韓国慶州ヒルトンホテルで開かれた前回会議の議事録を確認した．続いて，慶州における第11回国際伝熱会議の事後報告が韓国側から行われた.収入272,360ドル，支出269,200ドル，残額3,160ドルは韓国機械学会負担金15,000ドルから差し引くことで了承された次ぎに第12回国際伝熱会議(Grenoble会議と略称)の準備状況がフランス側から報告された．組織委員長Prof. Jean-Marc Delhaye，論文委員長Prof. Jean Taine，事務局長Prof. Jean-Bernard Saulnier，会期は2002年8月18～23日,会場はAtria World Trade Center Grenoble,公式オーガナイザはフランス熱力学会,参加登録料500ユーロ，内容予定はキーノート講演33編，ポスター論文発表565編(6室並行，54セッション),パネル計論４件,日程は2000年８月１日Web Site開設(Prof. Jean Taine, Ecole Centrale Paris, 92295 Chatenay-Malabry Cedex, France, ihtc12@em2c.ecp.fr),同９月10日Heidelbergにてアセンブリ委員会の承認を経て論文割り振り教決定,フランス側組織委員会(Steering Committee, SCと略記)はWeb Siteに2001年1月1日第1次アナウンスメント発表，同時にInter- national Scientific Committee(ISCと略記)メンバー宛にポスター送付,同4月1日までに各国ISCメンバーは割り振られた数の少なくとも2倍の数のキーノート講演候補者をSC宛に通知，同５月３０日までに論文発表希望者は各国のISCメンバ一宛にアブストラクトを送付，同７月１日までに各国ISCメンバーはアブストラクトを審査後，探択決定した著者宛に論文作成ガイドラインを送付，同日までにSCはキーノート講演者を決定しISCに通知，同９月１５日までに各国ISCメンバーはSCに探択論文のアプストラクトを送付, SCは著者宛に論文を各国ISCメンバー宛に送付するよう依頼,同12月1日に最終アナウンスメントを送付, 2002年2月15日までに各国ISCメンバーは,キーノート講演前刷及び担当した探択論文の前刷を,論文発表者の一覧表とともにSC宛に送付のこと，その場合この著者一覧表は，e-mailに添付して同時にSC宛送付のこと，同３月１５日にSCは詳細な第１次プログラムをWeb Siteに発表，同７月１日までに印刷所は最終プログラム及びアブストラクトの印刷物を完成しGrenobleに納品，同日までに出版社はCD-romをGrenobleに納品，同８月１８日～２３日本番,という日程で進行させる予定とのフランス側の提案を大筋において了承したが，CD-romのみというフランス側の提案に反対が多く,出版社から後でProceedingsを出版させる方向で見積もりを提出させ，後刻検討することとなった．日本の論文割り振り数は，第11回が割り振り90に対し実数99編であったが，第12回の事務局案は80編であったので， 90乃至95編に増やすよう申し入れた.キーノートは33編のうち５編(推薦数は10編)が割り振られており，この数はバランス上妥当と考えたので受け入れた．因みに前回は３編であった．次の議題はオーストラリアに於ける第13回国際伝熱会議(Sydney 2006)の準備状況報告である．公式オーガナイザはAustralasian Fluids and Thermal Engi- neering Society，組織委員長Emeritus Professor Gra- ham de Vahl Davis, UNSW,論文委員長Emeritus Pro- fessor Brian Milton, UNSW、事務局長Prof. Eddie Leonardi, UNSW、会期は2006年8月20日~26日、会場はSydney Convention and Exhibition Centre、併設展示会8月23日~25日に250~300社の出展予定、という。次にメンバーシップ審査委員会(委員長Prof. R.J.Goldstein,委員Prof. K.Stephan,同Prof. M.Hirata)の議を経て南アフリカのアセンブリメンバーへの参加申請が審議された．伝熱研究のアクティビティ，前第11回会議への参加者数6名などの実績について審議され，投票の結果，贊成10票，反対15票で否決された.併せてメンバーシップ審査委員会について審議され，内規を変更して，今後は元会長で構成するAdmission Committeeで予備的に審議した後アセンブリ会議にかけられることとなった。この委員会はProf. R.J.Goldstein. Prof. G.Hewitt, Prof. S.Sidemanで構成される.次回のアセンブリ会議はGrenoble会議の初日， 2002年8月18日に行うことを決め審議を終了した.日本の伝熱研究者各位に，前回会議と同様，質・量ともに世界に冠たる日本の伝熱研究のレベルを存分に示して下さるようお願いしてアセンブリ会議出席の報告とする.
日本伝熱学会学術賞・技術賞・奨励賞公募のお知らせ日本伝熱学会の内規に基づき，学術賞，技術賞，および奨励賞が設けられています．つきましては，下記の要領に従って本年度の募集を行いますので，自薦，他薦を問わず応募下さいますようお願い申し上げます. • 学術賞の対象は,原則として,最近５回の日本伝熱シンポジウムにおいて発表し, Thermal Science and Engineering誌またはその他の国内外で審査のある論文集に掲載された優秀な伝熱研究論文とする.技術賞の対象は，公表された優秀な伝熱技術とする．奨励賞の対象は原則として,最近２回の日本伝熱シンポジウムにおいて優秀な論文を発表した若手研究者で，発表時に大学院生，またはこれに準ずる者（大学卒業後５年以内の者）とする．学術賞および奨励賞の対象資格は，原則として本会会員に限る．学術賞は２件程度，技術賞は１件程度，奨励賞は４件程度とする． •学術賞・技術賞・奨励賞の選考は,「表彰選考委員会」が「日本伝熱学会賞審査・選考方法内規」によって行う.表彰選考委員は，公募の他に学術賞・技術賞・奨励賞候補を推勝することができる． 1通所定用紙「日本伝熱学会　学術賞・技術賞・奨励賞　申請書・推薦書」論文抜刷または技術内容参考資料6部日本伝熱シンポジウム講演論文集抜刷6部〒 501-1193　岐阜市柳戸１－１岐阜大学工学部機械システム工学科館　田　雅　弥　　宛TEL  058-293-2530  FAX 058-230-1892 5. 提出期限:平成13年1月9日(火)必着6.問い合わせ先:提出先に同じお　知　ら　せ日本伝熱学会学術賞・技術賞・奨励賞中請書・推薦書申請者・推薦者名_論文題名または: _技術名刊行物名または: _技術内容(論文抜刷または技術内容参考資料6部添付)受賞候補者（氏名・勤務先・職名・代表者の連絡先住所，Tel. No, Fax. No,本会会員資格の有無，氏名には振り仮名をお付け下さい)代表研究者：_ ___ ___共同研究者：_関連研究の伝熱シンポジウム発表論文題名：講演発表：第申請・推薦理由：_ ___
日：平成13年５月23日 (水) ～ 25日 (金)既準場: 大宮ソニックシティ(〒331-8669埼玉県大宮市桜木町1-7-5, TEL : (048)647-4111) AR (平成13年５月１日より「さいたま市」に市名変更)研究発表申込締切：平成13年1月19日(金)論文原稿締切:平成13年3月9日(金)参加事前申込締切：平成13年4月13日(金) •一般申込みによる通常セッション及び企画部会によるオーガナイズドセッションを実施します. ・講演は1題目につき20分（発表10分，質疑10分）の予定です. ・原則としてWebからの申込み_ (ホームページ: http://ht.t.u-tokyo.ac.jp）と致します．ホームページにある申込書に必要事項を記入し，お申込下さい．その後，申込受理のメールがお手元に届きましたら，それに記載されている申込番号，第一著者名，題目（最初７文字程度），必要経費（講演申込整理費のみ，or講演申込整理費＋ファイル変換代行作業料)，依頼人の住所，及び振込金額をご記入の上，下記の口座に申込整理費3,000円を郵便振替でご送金下さい（本号均込みの払込取扱票をご利用下さい）. •なお、Webからの申込ができない場合は本号に掲載の研究発表申込書に必要事項を記入し、下記送付先にFAXで発表申込書をお送り下さい．その後，申込受理のFAXがお手元に届きましたら，上記Web申込の場合と同様, FAXに記載されている申込番号,第一著者名,題目(最初７文字程度),必要経費(講演申込整理費のみ、or講演申込整理費＋ファイル変換代行作業科）、依頼人の住所、及び振込金額をご記入の上，下記の口座に申込整理費3,000円を郵便振替でご送金下さい（本号快込みの払込取扱票をご利用下さい). FAX申込時の送付先及び問合せ先：工業技術院機械技術研究所宗像鉄維宛て  FAX:(0298)61-7240講演申込整理費払込先：口座番号：00290-9-54822。口座名称：第38回日本伝熱シンポジウム実行委員会•申込整理費を講演会当日に会場にて支払われる場合は4,000円とさせて頂きます. ・講演発表申込は，講演者1名につき1題目とさせて頂きます. •発表の採否及びセッションへの振分けは，実行委員会にご一任願います。 ・講演論文集は原寸大のオフセット印刷及びCD-ROMとして作製致します．論文の長さは，１題目当たりA4用紙2ページとし，作成フォーマットは前回と同様の予定です（2段組×片側26字× 60行）. •執筆要綱は，会誌「伝熱」（平成13年1月号）及びホームページに掲載致します. ・本シンポジウムでは講演論文集をCD-ROM化するため，論文原稿は原則としてPDFのファイル形式で提出して戴きます．但し、どうしてもPDFでの提出が困難な場合（紙による提出、PDF以外のファイル形式による提出）には，論文申込整理費のほかに別途変換作業科(5,000円)を申し受けて，実行委員会がPDF化（スキャナー変換を含む）を代行致します．CD-ROM化のため印刷日程が厳しくなっており，原稿締切りに遅れると論文集とCD-ROMには載りませんのでご注意下さい.知   •論文の電子化（PDF化）については，下記文献が参考になります. “小林, Microsoft Wordを使った論文作成方法,日本機械学会誌, 103–979(2000.6), pp. 396–403." “川野・岡本，論文の電子化について（より良いPDFを作るために），可視化情報，20-77(2000.4), pp. 27-32." •講演申込整理費：3,000円（当日会場支払：4,000円） ・シンポジウム参加費（論文集代は含みません、日本伝熱学会会員、非会員共）:一般（事前申込：8,000円，会場申込：10,000円），学生（事前申込：4,000円，会場申込：5,000円） •講演論文集：日本伝熱学会会員：無料， 非会員：8,000 円• PDFファイル変換代行作業料（スキャナー変換含む、日本伝数学会会員、非会員共）: 5.000 円・開催日：平成13年5月24日（木） ・会場：大宮ソニックシティ 市民ホール•参加費：一般（事前申込：8,000円，会場申込：10,000円，事前申込、会場申込共に夫婦同伴者１名無料）学生(事前申达: 4,000円，会場申达: 5,000円） •研究発表申込の取消及び論文原稿の差替は，準備と運営に支障をきたしますのでご遠慮下さい. •論文の題目と著者名が発表申込時と論文提出時において相違ないようお願い致します. •発表に用いる機器は原則としてOHPのみとさせて頂きます. ・論文原稿及び別刷りは返却，配布致しません． •その他ご不明な点がありましたら下記までFAXまたはE-mailにてお問い合わせ下さい.第38回日本伝熱シンポジウム実行委員会東京農工大学工学部機械システム工学科内FAX: (042)388-7088, E-mail: htsymp38@mmlab.mech.tuat.ac.jp第38回日本伝熱シンポジウム実行委員会委員長　望月　貞成
（ ） http://ht.t.u-tokyo.ac.jp
l.(a − 6), 2.(A − 2), 3.(O − 4), 4.(b − 2), 5.(K − 2) _);._.2..(._________);._.3..(._ OS-1 OS-2 OS-3
Web http://ht.t.u-tokyo.ac.jp Web FAX    (0298) 61-7240 Joint Seminar Recent Numerical Approaches in Advanced Heat Transfer May 14 and 15, 2001  Seoul, Korea under the joint sponsorship of Heat Transfer Society of Japan and Thermal/Flow Control Research Center of Korea Institute of Science & Technology(KIST) co-chaired by T. Ito(Kyushu Univ.) and J. M. Hyun(KAIST) Approximately 15 presentations have been organized by the chairmen, and an additional 5 presentations may be accommo- dated. Those who wish to present papers are asked to submit by e-mail abstracts in English (around 500 words), by Febru- ary 28, 2001, to the seminar secretaries shown below. Y. Takata (Kyushu Univ.) fax +81-92-642-3400 takata@mech.kyushu-u.ac.jp B. H. Kang (KIST) fax +82-2-958-5689 bhkang@kistmail.kist.re.kr S. J. Kim (KAIST) fax +82-42-869-3210 sungjinkim@cais.kaist.ac.kr日本伝熱学会研究会「マイクロマシンと熱流体」第５回会合のお知らせ(電気学会　マイクロマシン研究会・マグネティックス研究会会合との同日開催)標記会合を下記の通り開催致します.本研究会(主査:笠木伸英(東大),副主査:藤田博之(東大),顧問:下山勲（東大）、幹事:村田章（農工大））はマイクロマシンに関連した熱流体を広く対象とし、意見交換、勉強会，見学会などを行っております．第５回会合は昨年同様に電気学会マイクロマシン研究会（委員長:佐藤一雄（名大）との同日開催とし，さらに今回は、電気学会Ａ部門マグネティックス研究会（委員長：角野圭ー（横国大）とも同日開催になります．本研究会委員以外の方の参加も歓迎致します．整理の都合上，参加ご希望の方は下記問い合わせ先までご連絡下さい. ○日時：平成１３年  ２月２８日(水)時間は未定○場所：立命館大学理工学部びわこ草津キャンパスコアステーション２階大会議室JR京都駅からJR西日本琵琶湖線(東海道線)にて20分，南草津駅下車.東口より近江鉄道バス「立命館大学行き」または「立命館大学経由飛島グリーンヒル行き」にて「立命館大学」下車，徒歩１分. URL: http://www.ritsumei.ac.jp/cinfo/map/index.htmをご覧下さい. ○参加費：無料（但し論文集は有料） ○プログラム午前・午後　電気学会マイクロマシン研究会・マグネティックス研究会講演マイクロマシンの製作法・性能評価・応用（６～８件程度）、マイクロ磁気デバイス・マイクロ磁気センサの材料・設計・製法・応用・製品化などに関する講演を行います.午後　伝熱学会研究会「マイクロマシンと熱流体」講演・加藤　宗（日立製作所　機械研究所） 「超音波を応用した微量液の非接触攪拌技術（仮題）」 ・鈴木雄二（東京大学大学院　工学系研究科） 「乱流制御に用いるマイクロデバイスの開発(仮題) [ ・渡慶次学(神奈川科学技術アカデミー　北森プロジェクト) 「マイクロ流体システムの化学分析への応用(仮題)」 ・浮田宏生（立命館大学光工学科） 「光圧回転体による流体攪拌（仮題）」 ○問い合わせ先：　〒184-8588東京都小金井市中町2-24-16東京農工大学工学部 機械システム工学科  村田章TEL 042-388-7089 FAX 042-385-7204 (学科共通) E-mail  murata@mmlab.mech.tuat.ac.jp URL http://www.mmlab.mech.tuat.ac.jp/micromachine国際宇宙ステーションの宇宙における組立が１９９８年１１月から開始され,本格的な宇宙環境を利用する時代を迎えようとしております.宇宙環境では,微小重力や高真空等地上では得がたい環境の下で,実験を実施することができます.我が国は，国際宇宙ステーション計画に独自の実験棟「きぼう」をもって参加しており，今後、このような宇宙環境を利用する機会が飛躍的に増大することになります.一方,地上での広範な分野の研究を宇宙環境を利用する実験に結びつけ,更には具体的な宇宙実験の提案・実施に至るまで地上における関連研究の蓄積が必要です.本制度は，このように「きぼう」を中心とした宇宙環境を利用する準備段階として，幅広い分野の研究者に研究機会を提供し，宇宙環境利用に関連する地上研究を推進することを目的としています.大学，国立試験研究機関，民間企業などの研究者の方々を対象に幅広く研究テーマを募集します. ●応募資格応募者の国籍，所属機関の所在地(国内外)は問いませんが，日本語による応募及び面接に対応できる語学力を必要とします.微小重力科学,微小重力物理学,生物科学,バイオメディカル,宇宙医学,宇宙科学,地球科学,宇宙利用技術開発の８分野・フェーズI研究（A）研究費：３千万円以下／年研究期間：1～3年•フェーズI研究（B）研究費：６百万円以下／年研究期間：1～3年•フェーズI研究（B）萌芽的研究 研究費：１５０万円以下（総額） 研究期間：１．５年•フェーズⅡ研究研究費：１千万円程度～１億円以下／年研究期間：１～３年財団法人日本宇宙フォーラム内に,有識者によって構成される委員会を設置し,応募された研究テーマの選定等を行います.研究形態としては，委託研究（大学，企業），共同研究（国立研究所），招聘研究の形態をとります. ●実験施設、試験施設の提供研究の必要に応じて調整の上，落下塔，航空機の微小重力模擬実験施設等の利用機会を提供いたします. ●スケジュール・応募受付締め切り:平成１３年２月２８日（水）消印有効•研究開始:平成１３年８月～財団法人日本宇宙フォーラム公募研究推進部 募集担当 電話: 0 3 － 3 4 5 9 － 1 6 5 3    F A X : 0 3 － 5 4 7 0 － 8 4 2 6 URL : h t t p ://www2.jsforum.or.jp/    E -mail:koubo@jsforum.or.jp
l. THERMAL SCIENCE AND ENGINEERING THERMAL SCIENCE AND ENGINEERING
2. JIS -第39期編集出版部会委員
平成12年11月30日第39期編集出版部会長編集出版事務局：〒 223-8522横浜市港北区日吉3-14-1慶應義塾大学理工学部システムデザイン工学科変田公一TEL: 045-566-1739 FAX: 045-566-1720 e-mail: hishida@sd.keio.ac.jp
本誌に掲載された著作物を複写したい方は、日本複写権センターと包括複写許諾契約を締結されている企業の従業員以外は、著作権者から複写権等の委託を受けている次の団体から許諾を受けて下さい。なお、著作物の転載・翻訳のような複写以外許諾は、直接本会へご連絡下さい。 T 107-0052東京都港区赤坂9-6-41乃木坂ビル3F学術著作権協会(TEL/FAX:  03-3475-5618)アメリカ合衆国における複写については、次に連絡して下さい。 Copyright Clearance Center, Inc.(CCC) 222 Rosewood Drive, Danvers, MA 01923 USA Phone : (978) 750-8400 FAX : (978)750-4744 In order to photocopy any work from this publication, you or your organization must obtain permission from the following organization which has been delegated for copyright for clearance by the copyright owner of this publication. Except in the USA The Copyright Council of the Academic Societies (CCAS) 41-6  Akasaka 9-chome, Minato-ku, Tokyo 107-0052 Japan TEL/FAX : 81-3-3475-5618 In the USA Copyright Clearance Center, Inc. (CCC) 222 Rosewood Drive, Danvers, MA 01923 USA Phone : (978) 750-8400 FAX : (978)750-4744に反数ISSN 1344-8692 (Journal of The Heat Transfer Society of Japan) Vol. 39, No.159 2000年11月発行本　伝　熱発行所社団法人町作低T 113-0034東京都文京区湯島2-16-16電話03(5689)3401 Fax. 03(5689)3401郵便振替　00160-4-14749 Published by The Heat Transfer Society of Japan 16-16, Yushima 2-chome, Bunkyo-ku, Tokyo-113, Japan Phone, Fax: +81-3-5689-3401