Transcript Document

反応性流体力学2008/6/4 引用
反応性流体力学2009/5/27 引用
Numerical Simulation of non-premixed flame
（diffusion flame） using a flamelet approach
Sorry, this presentation
is Japanese only…
Taniguchi Lab. IIS, The Univ. of Tokyo
Taniguchi Lab. meeting
2004/10/15
Shigefumi TOKUDA
1
拡散火炎モデル（Laminar Flamelet model）
対象とする火炎の構造が代表的な層流拡散火炎の
構造と同様であると仮定
仮定元となる層流拡散火炎のデータが必要
乱流火炎の計算を流れ場の計算から分離
化学的非平衡状態も表現可能
拡散火炎の火炎構造
Taniguchi Lab. IIS, The Univ. of Tokyo
2
保存スカラによる拡散火炎面の表現
化学種に関する支配方程式⇒元素の質量分率
各化学種の拡散係数が等しい（Scα=Sc）とする
規格化（燃料側でξ=1，酸化剤側でξ=0）
Taniguchi Lab. IIS, The Univ. of Tokyo
3
保存スカラによる拡散火炎面の表現
エネルギ方程式
Lewis数Le（=Sc/Pr）が１である（Scα=Sc=Pr）
低Mach数近似
輻射・壁面での伝熱効果無視
規格化（燃料側でξ=1，酸化剤側でξ=0）
Taniguchi Lab. IIS, The Univ. of Tokyo
4
保存スカラによる拡散火炎面の表現
各化学種の質量分率とエンタルピの境界条件が
同一であるとき，
元素の質量分率の輸送方程式とエネルギ方程式は
ただ一つのスカラξによって表現できる．
拡散火炎の形状はξの化学量論混合比の等値面ξst
で表される
ξ：混合分率
※化学量論混合比・・・燃料と酸化剤がすべて反応し生成物
となる燃料・酸化剤の割合
Taniguchi Lab. IIS, The Univ. of Tokyo
5
Flamelet data の作成
CHEMKINによる対向流拡散火炎の計算結果などか
ら，混合分率と密度・温度との関係を混合分率の多項
式関数として与える
ξst=0.03（H2:N2=100:0），0.3 （H2:N2=50:50）
Taniguchi Lab. IIS, The Univ. of Tokyo
6
水素-空気層流拡散火炎での検証
flamelet approachにもとづく燃焼モデルの検証
乱流拡散の影響を受けない層流場にて素反応解析の結果
と比較
素反応解析により得られた元素の質量分率から算出したξ
とパッシブなスカラとして計算したξの比較
ξ方程式導出の際に用いた仮定の妥当性の検証
素反応解析により得られた密度・温度とξからFlamelet
Dataを介して算出した密度・温度の比較
燃料濃度の違いによる影響
Taniguchi Lab. IIS, The Univ. of Tokyo
7
計算領域と計算条件
fuel
co-flow
fuel rate
H2+N2
Air（N2+O2）
6.7 m/s ~
co-flow rate
fuel nozzle diameter
Reynolds number
1.0 m/s
3.2×10-3 m
100~
Prandtl number
total cell number
total vertex number
0.75
19380
20820
詳細反応モデル（9化学種・21素反応
式）
Taniguchi Lab. IIS, The Univ. of Tokyo
8
詳細反応モデル（9化学種・21素反応式）
化学種
N2
O2
H2
H 2O
OH
Lewis数
1
1.11
0.3
1.12
0.73
化学種
O
H
H2O2
HO2
-
Lewis数
0.7
0.18
1.12
1.1
-
Taniguchi Lab. IIS, The Univ. of Tokyo
9
計算結果 (fuel H2:N2=100:0）
速度ベクトル分布 [m/s]
H2密度分布
[kg/m3]
Taniguchi Lab. IIS, The Univ. of Tokyo
温度分布 [K]
O2密度分布
[kg/m3]
H2O密度分布
[kg/m3]
OH密度分布 [kg/m3]
10
計算結果（混合分率分布） H :N =100:0，U
2
混合分率軸方向分布 [-]
2
fuel=6.7m/s
混合分率半径方向分布（Z = 1D） []
赤：詳細反応（N）
緑：詳細反応（O）
青：詳細反応（H）
ξ算出の際に着目する元素の
選択による違いは小さい
混合分率半径方向分布（Z = 5D） [-]
Taniguchi Lab. IIS, The Univ. of Tokyo
11
計算結果（混合分率分布） H :N =100:0，U
2
混合分率軸方向分布 [-]
2
fuel=6.7m/s
混合分率半径方向分布（Z = 1D） [-]
赤：詳細反応（N）
緑：パッシブスカラ（Sc=0.225（=ScH））
青：パッシブスカラ（Sc=0.75（=ScN））
桃：パッシブスカラ（Sc=0.8325 （=ScO））
混合分率半径方向分布（Z = 5D） [-]
Taniguchi Lab. IIS, The Univ. of Tokyo
Sc=0.75付近での結果がよく一致
12
計算結果（混合分率分布） H :N =100:0，U
2
混合分率軸方向分布 [-]
2
fuel=6.7m/s
混合分率半径方向分布（Z = 1D） [-]
赤：詳細反応（N） （Sc=0.75, Pr=0.75）
緑：パッシブスカラ（Sc=0.75, Pr=0.75）
青：パッシブスカラ（Sc=1.5, Pr=1.5）
桃：パッシブスカラ（Sc=0.5, Pr=0.5）
混合分率半径方向分布（Z = 5D） [-]
Taniguchi Lab. IIS, The Univ. of Tokyo
Pr=Sc の違いによる影響は小さい
13
計算結果（密度分布） H :N =100:0，U
2
軸方向分布
2
fuel=6.7m/s
半径方向分布（Z = 1D）
赤：詳細反応（N）から求めた混合分率ξから
緑：パッシブスカラ（Sc=0.75）としたξから
青：詳細反応計算により
求めた密度
半径方向分布（Z = 5D）
Taniguchi Lab. IIS, The Univ. of Tokyo
14
計算結果（温度分布） H :N =100:0，U
2
軸方向分布
2
fuel=6.7m/s
半径方向分布（Z = 1D）
赤：詳細反応（N）から求めた混合分率ξから
緑：パッシブスカラ（Sc=0.75）としたξから
青：詳細反応計算により
求めた温度
半径方向分布（Z = 5D）
Taniguchi Lab. IIS, The Univ. of Tokyo
後流域において温度分布に相違
→対向流拡散火炎と同様の特性はない
15
燃料濃度の影響（混合分率分布）
H2:N2=100:0
H2:N2=50:50
燃料を希釈（H2:N2=100:0→50:50）すると・・・
ξstが増加（0.03→0.3）
火炎面の広がりは小さい
噴流の広がりも小さい
Taniguchi Lab. IIS, The Univ. of Tokyo
16
計算結果（混合分率分布） H :N =50:50，U
2
混合分率軸方向分布 [-]
2
fuel=6.7m/s
混合分率半径方向分布（Z = 1D） []
赤：詳細反応（N）
緑：パッシブスカラ （Sc=0.75）
混合分率半径方向分布（Z = 5D） [-]
Taniguchi Lab. IIS, The Univ. of Tokyo
17
計算結果（密度分布） H :N =50:50，U
2
軸方向分布
2
fuel=6.7m/s
半径方向分布（Z = 1D）
赤：詳細反応（N）から求めた混合分率ξから
緑：パッシブスカラ（Sc=0.75）としたξから
青：詳細反応計算により
求めた密度
半径方向分布（Z = 5D）
Taniguchi Lab. IIS, The Univ. of Tokyo
18
計算結果（温度分布） H :N =50:50，U
2
軸方向分布
2
fuel=6.7m/s
半径方向分布（Z = 1D）
赤：詳細反応（N）から求めた混合分率ξから
緑：パッシブスカラ（Sc=0.75）としたξから
青：詳細反応計算により
求めた温度
半径方向分布（Z = 5D）
Taniguchi Lab. IIS, The Univ. of Tokyo
19
まとめ
詳細反応計算で得られた各元素より算出した混合分率ξ
はどの元素を基にしてもほぼ等しい
各化学種の拡散係数が等しい（Scα=Sc）とした仮定の妥当性
パッシブスカラとして計算した混合分率ξはSc＝Pr
とすると詳細反応計算の結果とよく一致
Le（=Sc/Pr）=1の仮定の妥当性
混合分率ξから平面火炎（対向噴流）のデータにより
予測した温度、密度分布には誤差あり
Flamelet近似の妥当性
火炎形状は燃料濃度（燃焼速度）に依存する
Flamelet近似でも依存性の予測はある程度可能
Taniguchi Lab. IIS, The Univ. of Tokyo
20