2016 Star Japanese Conference

Size: px

Start display at page:

Download "2016 Star Japanese Conference"

やすもりあみおか
4 years ago
Views:

heat transfer characteristics in transport of the thermal fluid using a vortex

1 [ 流体技術 : Room C ] STAR Japanese Conference 2016 YOKOHAMA JUNE 9-10, 2016 空気砲を用いた熱流体の輸送における拡散および熱伝達特性の LES 解析 Large Eddy Simulation of the diffusion and heat transfer characteristics in transport of the thermal fluid using a vortex ring 福岡大学工学部機械工学 : 赤木富士雄 Fukuoka University, Department of Mechanical Engineering Ph D, Fujio Akagi

2 本講演の概要研究目的空気砲の正体である渦輪を用いて熱流体を拡散抑制しながら輸送する方法の確立検討内容渦輪内の温度分布が渦輪の生成過程粘性拡散および対流熱伝達に及ぼす影響を調べた STAR-CCM+ を用いて熱流体を供給した渦輪を予測 ( 渦度分布温度分布密度分布粘性分布 ) 計算は Large Eddy Simulation を実施結果渦核内では温度分布に起因した粘性分布と密度分布が形成両者の分布は相互に作用して渦輪の拡散および熱伝達特性に影響

3 Contents 1. 研究背景目的 2. 数値シミュレーション手法 3. 解析結果結果の妥当性の確認渦輪内への熱流体供給方法の検討温度分布が渦輪に及ぼす影響 4. まとめ

流体や微小粒子を拡散抑制しながら目標箇所までピンポイント輸送目標箇所を輸送物で局所的に満たす正面図渦輪の可視化写真

4 1. 研究背景渦輪の流体力学的特徴煙のドーナツの正体 => 渦輪渦核内に流体および微小粒子を閉じ込めた状態で自己の誘起速度によって移動できる側面図局所集中輸送法としての利用流体や微小粒子を拡散抑制しながら目標箇所までピンポイント輸送目標箇所を輸送物で局所的に満たす正面図渦輪の可視化写真 (Akhmetov,D. G., Vortex Rongs,Springer) 特徴輸送流体が粘性拡散することで生じる混合損失に起因した消費エネルギーの増大を抑えることができる

5 1. 研究背景渦輪輸送に関する過去の取り組み産業応用 : 空調機の送風方法エンターテインメント応用 : 香り輸送, バーチャルディスプレイ基礎研究 : 流動境界層制御利用者周辺を選択的に積極冷却渦輪送風空調応用の概念図空調輸送方法として実用化されるまでには至っていない実用化の課題 (1) 輸送に適した渦輪の生成方法条件が不明 (2) 輸送したい流体を渦輪内に格納する最適方法が不明 (3) 渦輪内の流体の粘性拡散熱伝達特性が不明

6 1. 研究目的課題研究目的渦輪を用いた熱流体の局所集中輸送法の開発課題 1. 輸送に適した渦輪の生成方法条件の解明循環が大きい ( 拡散しにくい ) 渦核の体積 ( 輸送容積 ) が大きい連続生成しても軸対称を維持 EFD 2. 熱流体 ( 輸送対象 ) を渦輪内に格納する最適手法条件の解明 3. 渦輪内の温度分布が渦輪の生成過程粘性拡散特性および熱伝達特性に及ぼす影響の解明 CFD

7 Contents 1. 研究背景目的 2. 数値シミュレーション手法 3. 解析結果結果の妥当性の確認渦輪内への熱流体供給方法の検討温度分布が渦輪に及ぼす影響 4. まとめ

8 2. 数値シミュレーション手法予備実験実験結果を基に渦輪の生成過程の特徴を把握計算格子および条件設定上の注意点を確認 U(t) U 0 0 U m =0 T t 図. 正弦波の流量変動 ( 脈動噴流 ) 粒子画像流速測定法 (PIV) 位相平均速度場を取得 [ 128 周期 ] 図. 実験領域の概略図

9 2. 数値シミュレーション手法正弦波形 : Re 0 = 2370, α = 20.2, Str = 単発生成 Nozzle 渦輪 ζ * = ζ d n /U 0 連続生成はく離渦輪 Nozzle 周期渦輪図. 速度ベクトルおよび無次元渦度分布噴出口では, 噴流のせん断層が巻き上がって渦輪が生成ノズル内では, 吸い込み流れの境界層がはく離して, はく離渦輪が生成はく離渦輪は, 周期渦輪の生成過程に影響を及ぼす

10 2. 数値シミュレーション手法 ζ * 周期短 Re 0 = 1300, α = 24.2, Str = 振幅大 Re 0 = 2370, α = 20.2, Str = Re 0 = 1500, α = 13.9, Str = Re 0 = 3250, α = 20.2, Str = はく離渦輪の生成挙動 ( はく離境界層の長さ ) は, 脈動条件により異なるシミュレーションの重要ポイント : せん断層とはく離渦輪の再現性

11 2. 数値シミュレーション手法ソルバーおよび解析条件 STAR - CCM+ Ver ~ ( 倍精度 ) 三次元非定常非圧縮 Large Eddy Simulation (LES) 熱伝達 ( エネルギー式 ) 重力項 ( 浮力の影響 ) を考慮空間離散化時間進行 Sub-Grid-Scale モデル二次精度中心差分二次精度陰解法 Dynamic Smagorinsky model 時間ステップ dt = sec (Cr < 0.5) Key Point 残差設定 ( 内部イタレーション :5 回 ) 瞬時の渦流れ場せん断層の予測 : LES はく離流れの予測 : Dynamic Smagorinsky model クーラン数 Cr 0.5

6.7 5 1 2. 数値シミュレーション手法計算領域およびメッシュノズル内の境界層とせん断層内にはメッシュを集中 (16 層 ) 壁面上には δy = 0.05mm (y+ < 1) トリムメッシュ部はアスペクト比 1.

12 数値シミュレーション手法計算領域およびメッシュノズル内の境界層とせん断層内にはメッシュを集中 (16 層 ) 壁面上には δy = 0.05mm (y+ < 1) トリムメッシュ部はアスペクト比 1.25 以下 Key Point: はく離長さの予測 Nozzle プリズムレイヤーメッシュ LayerCell region : 16 層トリムメッシュ TrimmedCell(Size:0.25) : 層 (0.2 mm) トリムメッシュ TrimmedCell(Size:0.5) : 15 層 (0.45 mm) 図. ノズル内のメッシュ配置図 Z 図. ノズル内のメッシュ図

13 2. 数値シミュレーション手法計算領域メッシュ外観図流入境界壁面境界壁面境界圧力境界圧力境界圧力境界子午面断面図 [ 総メッシュ数 : 640 万 ] 外観図 (1/4 のみ表示 )

14 2. 数値シミュレーション手法メッシュ生成からデータ解析までの流れ計算領域メッシュ生成数値計算ポスト処理データ解析検証予備計算本計算一次データ二次データ ( ワークステーション Intel Xeon 16 core ) Power Session STAR-CCM+ ( 九州大学 Primergy CX400 4 node) Post On Demand オリジナルコード可視化ソフト計算時間 : ( 例 ) 計算 step 数 100,000 ( 変動 10 周期 ) ワークステーション => 約 3 weeks Primergy CX400 => 約 1 weeks

15 Contents 1. 研究背景目的 2. 数値シミュレーション手法 3. 解析結果結果の妥当性の確認渦輪内への熱流体供給方法の検討温度分布が渦輪に及ぼす影響 4. まとめ

3. 解析結果 : 結果の妥当性の確認単発渦輪パルス状の流量変動による単発渦輪の生成の LES V 3 V 2 V

16 3. 解析結果 : 結果の妥当性の確認単発渦輪パルス状の流量変動による単発渦輪の生成の LES V 3 V 2 V 1 V 1 V 3 V 2 (a) EFD Gharib, M., et al., A universal time scale for vortex ring formation, 1998, Journal of Fluid Mechanics, 360, pp (b) LES 図. 渦度分布表. 渦輪の特性値の比較 D 1 (cm) D 2 (cm) D 3 (cm) Z 1 (cm) Z 2 (cm) Z 3 (cm) Γ (cm 2 /s) EFD LES

17 3. 解析結果 : 結果の妥当性の確認周期渦輪 : 正弦波形 - Re 0 = 2370, α = 20.2, Str = (a) EFD ( 位相平均 ) 図. 無次元渦度分布 (LES は実験の計測点におけるデータを使用 ) (b) LES せん断層の巻き上がり過程を再現できているはく離渦輪の挙動もよく再現できている定性的な再現性は良好

18 3. 解析結果 : 結果の妥当性の確認 PIV LES 0.8 PIV LES ω θ /ω max ω r /ω max Jetside z * =r /a z core r * =r /a r core Upstream Fig.Comparison of vorticity distribution between LES and EFD Downstream ζ * peak 30 EFD (D g /d n =4.63) 25 EFD (D /d = ) g n CFD-Trimmed Cell Model t * =t/t 図. 渦核中心点上の無次元渦度の時間変化 Z/d n 8 EFD (D g /d n =4.63) 7 EFD (D /d = ) g n 6 CFD-Trimmed Cell Model D v /d n EFD (D /d =4.63) g n EFD (D /d = ) g n CFD-Trimmed Cell Model t * =t/t t * =t/t 図. 渦輪の渦核中心位置の時間変化図. 渦輪の直径の時間変化 >> LES の結果は定量的にも実験結果と良い一致

19 3. 解析結果 : 結果の妥当性の確認脈動条件と渦輪の循環 ( 強さ ) の関係 CASE:B EFD 水渦輪 : 実験 - 予測値水渦輪 :CFD LES 全周モデル空気渦輪 :CFD 全周モデル CASE: B Re Γ /Re 0 1 CASE:A CASE: A Str = 0.05 Str 本 LES の結果は渦輪の生成過程を精度よく再現している渦輪の循環 ( 強さ ) の生成条件に対する変化は水と空気で同じ循環が極大となる条件が存在 Str = 0.05

20 3. 解析結果 : 結果の妥当性の確認渦度層の等値面分布図拡散に利用図. はく離渦輪の影響が強い条件輸送に利用図. 強い渦輪が出来る条件

21 Contents 1. 研究背景目的 2. 数値シミュレーション手法 3. 解析結果結果の妥当性の確認渦輪内への熱流体供給方法の検討温度分布が渦輪に及ぼす影響 4. まとめ

22 3. 解析結果 : 渦核内への熱流体供給方法の検討熱流体の輸送試験 80 の温水を渦輪に供給して20 の水中を輸送 1. 熱流体を正弦波の脈動噴流として噴出する条件 ( 最も単純な方法 ) 2. ノズル内の壁面上に熱源を設置する条件 ( せん断層を加熱 ) 熱源をノズル内壁面に設置熱源をノズル内外の壁面に設置 3. 熱流体を噴流のせん断層に供給する条件 ( せん断層に直接供給 ) 注入流路幅 :H = 0.5mm H = 1.5mm

23 3. 解析結果 : 渦核内への熱流体供給方法の検討熱源 :80 (a) 熱流体を噴流として噴出 (b) 熱源をノズル内壁面に設置熱源 :80 (c) 熱源をノズル内外の壁面に設置

5mm 30 H = 0.5mm (b) 渦度分布 20 0 1 2 3 4 5 6 7 8 渦核中心到達距離 z/d n 図.

24 3. 解析結果 : 渦核内への熱流体供給方法の検討 (d) 熱流体をせん断層に直接供給した条件 ( 流路幅 0.5mm の条件 ) 20 (a) 温度分布熱水 :80 渦核中心温度 (K) H = 1.5mm 30 H = 0.5mm (b) 渦度分布渦核中心到達距離 z/d n 図. 渦核中心温度の変化せん断層内に熱流体を直接供給する方法が最も有効温度はすぐに低下 => 渦核内の熱容量と熱伝達との関係渦核内に熱流体をどれだけ供給できるかが非常に重要

25 3. 解析結果 : 渦核内への熱流体供給方法の検討供給条件 CFD EFD CFD EFD CFD EFD Q sup /Q id Q sup /Q id Q sup /Q id A * = A sup /A noz θ (deg) U sup /U 0 供給条件には最適条件が存在 ( 最大で渦核体積の 57%) 実験結果との差が大きい ( 差の要因は不明 ) 更なる検討が必要

26 Contents 1. 研究背景目的 2. 数値シミュレーション手法 3. 解析結果結果の妥当性の確認渦輪内への熱流体供給方法の検討温度分布が渦輪に及ぼす影響 4. まとめ

27 3. 解析結果 : 温度分布が渦輪に及ぼす影響渦輪内による熱流体 ( 暖冷 ) の輸送試験渦核内に熱流体 ( 周囲との温度差 dt=20, -20 ) を供給供給はせん断層に直接供給する方法重力項 ( 浮力の影響 ) を考慮単発渦輪を下向きに生成した条件で検討 0 Y Z U sup = 0.1 U 0 θ = 15 W = 10 mm Re 0 = 6250 α = 22.6 Str = 0.052

28 3. 解析結果 : 温度分布が渦輪に及ぼす影響無次元渦度分布温度分布 (a) dt = - 20 (b) dt = 0 (c) dt = + 20 冷暖に関わらず温度差があるとせん断層の渦度は強くなる結果として渦輪の循環は大きくなる => 移動速度が高くなる粘性拡散および熱伝達は温度差によって異なる粘性拡散の速度と熱伝達の速度は同程度

29 3. 解析結果 : 温度分布が渦輪に及ぼす影響渦輪の循環差なし dt=10 dt=-10 差なし dt=10 dt=-10 冷暖に関わらず温度差があるとせん断層の渦度は強くなり結果として渦輪の循環が大きくなるこの結果を利用すれば温度差を利用して渦輪の強さを強めることができる空調送風への利用の点では有利な結果

30 Contents 1. 研究背景目的 2. 数値シミュレーション手法 3. 解析結果結果の妥当性の確認渦輪内への熱流体供給方法の検討温度分布が渦輪に及ぼす影響 4. まとめ

31 4. まとめ STAR-CCM+ を用いて Large Eddy Simulation により熱流体を供給した渦輪温度分布が渦輪の生成過程粘性拡散および対流熱伝達に及ぼす影響を調べた渦輪の循環 ( 強さ ) の生成条件に対する変化は水と空気で同じ循環が極大となる条件が存在 Str = 0.05 熱流体の供給はせん断層内に直接供給する方法が最も有効輸送の実現には渦核内に熱流体をどれだけ供給できるかが非常に重要渦核内では温度分布に起因した粘性分布と密度分布が形成両者の分布は相互に作用して渦輪の拡散および熱伝達特性に影響温度差を与えると渦輪の循環は大きくなる

32 謝辞本論文の完成までには, 多くの方々の多大なるご指導とご協力を賜りました. ここに, 衷心より御礼を申し上げます. 福岡大学教授山口住夫先生准教授安東洋一先生福岡大学流体工学実験室元修士生川畑慶佑瀬戸口翔太林真一本研究では以下の援助を賜りました. ここに記して御礼を申し上げます年度日本私立学校振興共済事業団学術研究振興資金年度原田記念財団研究助成

粒子画像流速測定法を用いた室内流速測定法に関する研究

粒子画像流速測定法を用いた室内流速測定法に関する研究可視化手法を用いた室内気流分布の測定法に関する研究 -PIV を用いた通風時及び空調吹出気流の測定 - T08K729D 大久保肇指導教員赤林伸一教授流れの可視化は古来より流れの特性を直感的に把握する手法として様々な測定法が試みられている近年の画像処理技術の発展及び PC の性能向上により粒子画像流速測定法 (PIV ) が実用化されている Particle Image Velocimetry