JSME TED Newsletter, No.78, 2016 TED Plaza 液噴流衝突時における飛散液滴生成に関する研究榎木光治電気通信大学助教情報理工学研究科知能機械工学専攻大川富雄電気通信大学教授情報理工学研究科知能機械工学専攻 oka

Size: px

Start display at page:

Download "JSME TED Newsletter, No.78, 2016 TED Plaza 液噴流衝突時における飛散液滴生成に関する研究榎木光治電気通信大学助教情報理工学研究科知能機械工学専攻大川富雄電気通信大学教授情報理工学研究科知能機械工学専攻 oka"

あおしそめや
4 years ago
Views:

TED Plaza 液噴流衝突時における飛散液滴生成に関する研究榎木光治電気通信大学助教情報理工学研究科知能機械工学専攻 enoki.koji@uec.ac.jp 大川富雄電気通信大学教授情報理工学研究科知能機械工学専攻 okawa.tomio@uec.ac.jp 1. はじめに原子炉から発生する熱エネルギーは, 冷却材を介して移動する.

( 神田他,2009) しかし, ナトリウムは化学的活性度が高く, 何かしらの理由で配管等の破損が生じ漏洩が起こると, 雰囲気中の酸素や湿分との化学反応により燃焼を生じ, 格納容器などの構造物に熱的負荷を及ぼす. このため, ナトリウム漏洩時における構造物の健全性評価をすることを目的として, ナトリウム燃焼数値解析コードの開発が進められている.

1 TED Plaza 液噴流衝突時における飛散液滴生成に関する研究榎木光治電気通信大学助教情報理工学研究科知能機械工学専攻大川富雄電気通信大学教授情報理工学研究科知能機械工学専攻 1. はじめに原子炉から発生する熱エネルギーは, 冷却材を介して移動する. 冷却材として, 熱伝導率が高く沸点が高い液体金属のナトリウムを使用するナトリウム冷却高速炉では, 冷却材に水を使用する軽水炉と比べて, 高温でかつ低圧でシステムを組むことができる. このため, ナトリウム冷却高速炉は軽水炉と比較して, 高温過熱蒸気サイクルが可能なため熱効率が上昇し, かつ各機器は薄肉構造にする事ができる.( 神田他,2009) しかし, ナトリウムは化学的活性度が高く, 何かしらの理由で配管等の破損が生じ漏洩が起こると, 雰囲気中の酸素や湿分との化学反応により燃焼を生じ, 格納容器などの構造物に熱的負荷を及ぼす. このため, ナトリウム漏洩時における構造物の健全性評価をすることを目的として, ナトリウム燃焼数値解析コードの開発が進められている. ここで, ナトリウムの燃焼は表面反応であるため, 仮にナトリウムの漏洩量が同じであっても, その漏洩や分布状況によって反応速度は大きく異なる. 漏洩したナトリウムが床等の上に液プールを形成しながらプール表面で燃焼を生じるプール燃焼よりも, 配管の破損箇所等から漏洩する際に噴流から生じる液滴, もしくは噴流が床面等構造物や液プールに衝突することで液滴を形成し, 液滴表面で燃焼を生じるスプレイ燃焼の方が, 単位体積当たりの表面積が多いため, 一般に激しい燃焼となる. このため, ナトリウム漏洩時の燃焼解析を行う上で, 漏洩量に対して液滴状となり飛散するナトリウムの総量や径分布を正確に評価することは, 極めて重要である. 本稿では, 最も多量の液滴を生成する液体噴流が固体面に衝突する際の液滴生成のメカニズム解明および飛散液滴量の見積もり法に関する研究について, 簡単に紹介する. 2. 実験方法および実験条件本研究では, 試験流体として水を用いた. 液滴生成のメカニズムを明らかにするために, 高速度カメラを用いて噴流や飛散する液滴の様子を詳細に観察した. 噴流の固体面衝突時に発生する飛散液滴量の計測に使用した液滴回収装置の写真を図 1 に示し, 飛散液滴の回収原理を模式化して図 2 に示す. 液滴回収装置はステンレス製で, 吸水性の高い紙製のウエスを, 飛散した液滴が

Fig.1 Photo of the droplet collection device Fig.2 Principle of droplet collection 全量吸収されてかつ滴り落ちないように数枚重ねて, 回収装置の枠組みに沿って取り付けた. 飛散した液滴を 30 秒間回収装置に吸収させた後, ウエスの重量変化を電子天秤で計測し, 回収液滴量を測定した.

2 Fig.1 Photo of the droplet collection device Fig.2 Principle of droplet collection 全量吸収されてかつ滴り落ちないように数枚重ねて, 回収装置の枠組みに沿って取り付けた. 飛散した液滴を 30 秒間回収装置に吸収させた後, ウエスの重量変化を電子天秤で計測し, 回収液滴量を測定した. ここで, 図 1 および図 2 の通り, 回収装置は中心角 30 度の扇形をしており, 電子天秤で測定した回収液滴量を 12 倍することで, 飛散した全液滴量 Sp としている. 以下で飛散量の議論をする際は, 飛散した全液滴量 Sp の事である. 表 1 に実験条件を示す. 実験は, ノズル径 1 4mm, 噴流流量 ml/min, ノズル出口から固体面に衝突するまでの高さ ( 衝突高さ ) mm の範囲で行っている. Fluid Table 1 Experimental condition Water Density = 998 kg/m 3, Surface tension = 72 mn/m at 20 C Nozzle Diameter D [mm] Flow rate Q [ml/min] Fall height L [mm] 62.5, , , 500

3. 噴流と飛散液滴可視化実験で観察された噴流を典型的な様相ごとに分類して, 図 3 に示す. 各図はノズル径 D = 4 mm, 衝突高さ L = 500 mm の条件で噴流流量 Q を変化させて観察していて, 流量が小さいものから順に (a) ~ (d) に並べている.

本実験研究では, 以下の通り, 観察された噴流を 4 つに分類および定義し, 噴流とその噴流に伴う飛散液滴の特徴をそれぞれ示す. (a) Drop: 噴流が断続的で, 概ね均一サイズの液滴列が固体面に衝突する. この噴流は, 噴流の流量が小さい場合に観察され, 噴流の液滴列が衝突する度に, 比較的多くの飛散液滴を生じる.

噴流流量が Lump よりも高い条件で観察されるが, 噴流の衝突で生じる飛散液滴は最も少ない. (d) Disturbance:Jet が観察される条件よりもさらに噴流流量を増加させると, 再び断続的な噴流となる. 噴流の形状には時間変動があり,Drop もしくは Lump に似た液滴列を形成して固体面に衝突する.

ここでは, ノズル径 D = 4 mm, 衝突高さ L = 500 mm の実験結果を示しており, 図 3 で示した各流量条件と一致する箇所を図 4 中にそれぞれ矢印で表している. ここで, ウェーバー数は We = V 2 D / で, 密度と表面張力はそれぞれ表 1 の値を使用している.

3 3. 噴流と飛散液滴可視化実験で観察された噴流を典型的な様相ごとに分類して, 図 3 に示す. 各図はノズル径 D = 4 mm, 衝突高さ L = 500 mm の条件で噴流流量 Q を変化させて観察していて, 流量が小さいものから順に (a) ~ (d) に並べている. また, 各条件で測定した飛散液滴量 Sp と, 全流量に対する飛散液滴の割合である飛散率 Sp* を, 流量条件 Q と合わせて各図に示す. 観察された噴流の様相は, 全般に Lin and Leitz (1998) の観察結果と良い一致を示した. 本実験研究では, 以下の通り, 観察された噴流を 4 つに分類および定義し, 噴流とその噴流に伴う飛散液滴の特徴をそれぞれ示す. (a) Drop: 噴流が断続的で, 概ね均一サイズの液滴列が固体面に衝突する. この噴流は, 噴流の流量が小さい場合に観察され, 噴流の液滴列が衝突する度に, 比較的多くの飛散液滴を生じる. (b) Lump: 噴流流量を Drop よりも増加させると, 連続的な噴流となる. 噴流は軸対称の波打った形状で,Plateau - Rayleigh 不安定性 ( 奥村,2003) の様相を呈する. 噴流から生じる飛散液滴は少ない. (c) Jet: 連続的な噴流で, ほぼ直線形状である. 噴流流量が Lump よりも高い条件で観察されるが, 噴流の衝突で生じる飛散液滴は最も少ない. (d) Disturbance:Jet が観察される条件よりもさらに噴流流量を増加させると, 再び断続的な噴流となる. 噴流の形状には時間変動があり,Drop もしくは Lump に似た液滴列を形成して固体面に衝突する. 噴流の衝突で生じる飛散液滴は最も多く, そして非常に多い. 飛散液滴と噴流流量の関係を, 飛散率 Sp* とウェーバー数 We を用いて比較した結果を, 図 4 に示す. ここでは, ノズル径 D = 4 mm, 衝突高さ L = 500 mm の実験結果を示しており, 図 3 で示した各流量条件と一致する箇所を図 4 中にそれぞれ矢印で表している. ここで, ウェーバー数は We = V 2 D / で, 密度と表面張力はそれぞれ表 1 の値を使用している. また V [m/s] は噴流流量 Q から算出した噴流速度である. 図 3 および図 4 に示すとおり, 飛散する液滴は, 噴流流量の増加で増大するとは限らず, 流量 (a) Q = 500 ml/min Sp = 36 ml/min Sp * = Drop mode (b) Q = 1500 ml/min Sp = 17ml/min Sp * = Lump mode (c) Q = 2000 ml/min Sp = 5 ml/min Sp * = Jet mode (d) Q = 3000 ml/min Sp = 617ml/min Sp * = Disturbance mode Fig.3 Typical patterns of the liquid jet with splash droplet

4 Fig.3 (b) Lump Fig.3 (a) Drop Fig.3 (d) Disturbance Fig.3 (c) Jet 1 Sp* L=500mm We Fig.4 Dependence of Sp* on We (D = 4 mm, L = 500 mm) 変化による噴流様相の変化によって飛散液滴量が変動し低流量もしくは高流量域で観察される断続的な噴流の時に飛散量が多く中流量で観察される連続的な噴流時には飛散が少なくなる傾向があることがわかったまた衝突高さ L についても飛散液滴量に大きな影響を及ぼすことが明らかになっている衝突高さが低いつまり噴流が衝突する固体面がノズルに近づくほど Jet などの連続的な噴流となる傾向があるため飛散液滴量は少なく一方で衝突高さが高いと drop などの断続的な噴流になる傾向があり飛散液滴量は多くなる以上ノズル径 D = 4 mm の実験を例に紹介したが他のノズル径においても噴流の様相および飛散液滴の定性的な傾向は同じである４飛散液滴量の相関式噴流と飛散液滴の関係を検討した結果飛散液滴量は以下の条件で変化すると考えられる 1 噴流の流量変化流速変化 V による噴流様相の変化 2 噴流の固体面との衝突周期 f 3 衝突高さ L そこで噴流の流量変化流速変化をウェーバー数 We で衝突周期をストローハル数 St で表し飛散率 Sp*を整理したその結果を図 5 に示すここで We 数は図 4 で定義したものと同じ We = V2D / また St = fd / V の f [1 / s]は衝突周期で観察から得られた値を使用しているこれらの結果から明らかなように飛散率 Sp*は We 数と St 数で良く整理できることが明らかになったまた各ノズルで衝突高さ L が高い条件は図 5 中に赤いプロットで示しているが前述の通り飛散率が多い傾向にあることがわかるしかし無次元整理した中に衝突高さは含まれていないこれは衝突周期は衝突高さと相関関係があるが今回衝突周期は実験値を用いているため衝突高さが式中に直接含まれなくても飛散率を整理できているものと考えられる今後の研究では衝突周期を正確に予測できる整理式の作成および試験流体を変えて実験を行い飛散液滴量を予測する整理式の精度を向上させる予定であるまた飛散液滴の径についてもナトリウム燃焼に重要なパラメータとなるため予測手法を確立していく必要がある

5 Fig.5 Correlation of Sp * using the Weber number We and the Strouhal number St 謝辞本研究は, 国立研究開発法人日本原子力研究開発機構 (JAEA) の受託研究として実施したもので, 同機構の大野修司博士と青柳光裕博士にお力添えいただいた. また, 電気通信大学の大箭直輝氏,Zhan Yi 氏にも協力いただいた. ここに記して感謝の意を表する. References (1) 神田誠, 梅田賢治, 三宅修平, 清水健男, 一宮正和, 山下清信, 望月弘保, 与熊本泰介, 岡芳明, 原子力プラント工学 (2009), p.202 (2) Lin, S. P., Reitz, R.D., Drop and spray formation from a liquid jet, Annual Review of Fluid Mechanics, Vol. 30 (1998), pp (3) 奥村剛, 表面張力の物理学 (2003),pp.10-11,pp

技術資料 JARI Research Journal OpenFOAM を用いた沿道大気質モデルの開発 Development of a Roadside Air Quality Model with OpenFOAM 木村真 1 Shin KIMURA 伊藤晃佳 2 Akiy

技術資料 JARI Research Journal OpenFOAM を用いた沿道大気質モデルの開発 Development of a Roadside Air Quality Model with OpenFOAM 木村真 *1 Shin KIMURA 伊藤晃佳 *2 Akiy 技術資料 176 OpenFOAM を用いた沿道大気質モデルの開発 Development of a Roadside Air Quality Model with OpenFOAM 木村真 *1 Shin KIMURA 伊藤晃佳 *2 Akiyoshi ITO 1. はじめに自動車排出ガスの環境影響は, 道路沿道で大きく, 建物など構造物が複雑な気流を形成するため, 沿道大気中の自動車排出ガス濃度分布も複雑になる.

JSME TED Newsletter, No.78, 2016 TED Plaza 液噴流衝突時における飛散液滴生成に関する研究 榎木光治 電気通信大学助教情報理工学研究科知能機械工学専攻 大川富雄 電気通信大学教授情報理工学研究科知能機械工学専攻 oka

JSME TED Newsletter, No.78, 2016 TED Plaza 液噴流衝突時における飛散液滴生成に関する研究榎木光治電気通信大学助教情報理工学研究科知能機械工学専攻大川富雄電気通信大学教授情報理工学研究科知能機械工学専攻 oka