名古屋大学最終講義
２０１５年３月６日
工学研究科２号館２４１講義室

教職員学生と歩んだ
伝熱・燃焼研究の38年
名古屋大学大学院 工学研究科
機械理工学専攻 機械科学分野

山下 博史
1

教職員学生と歩んだ伝熱・燃焼研究の38年<< 塑性力学 >>学部4年・修士課程の時代 1971年4月～1974年3月
機械学科 材料力学講座
大橋義夫先生, 川嶋紘一郎先生，徳田正孝先生<< 伝熱工学１ >> 博士課程・助手の時代 1974年4月～1979年3月
機械工学第二学科 内燃機関及びガスタービン講座
泉亮太郎先生，山口誉起先生，加賀定先生，木下祥次先生
大岩紀生先生，高橋俊夫技官<< 伝熱工学２ >> 講師・助教授時代 1979年4月～1988年3月
泉亮太郎先生，加賀定先生，櫛田玄一郎助手，早川正人技官（高浜平七郎先生，藤田秀臣先生，河村鈞先生，藤城明男技官，西村新璽技官）<< 燃焼工学１ >> 助教授時代 1988年4月～1998年3月
機械工学科 機械エネルギー工学講座  機械情報システム工学科 環境情報システム講座
竹野忠夫先生，櫛田玄一郎講師，西岡牧人助手，早川正人技官<< 燃焼工学２ >> 教授時代 1998年4月～2004年3月
機械情報システム工学科 ヒューマンシステム工学講座  環境情報システム講座
趙黛青講師，中村祐二助手，早川正人技官，（新井紀男教授 ）<< 燃焼工学３ >> 教授時代 2004年4月～現在
機械理工学専攻 機械科学分野 環境・エネルギー工学講座 伝熱・燃焼工学研究グループ2
山本和弘准教授，林直樹助教，早川正人技術職員

学部4年・修士課程の時代<< 塑性力学 >>学部4年・修士課程の時代 1971年4月～1974年3月
機械学科 材料力学講座
大橋義夫先生, 川嶋紘一郎先生，徳田正孝先生
・単純負荷により予ひずみを受けた軟鋼の降伏面と応力-ひずみ曲線：
山下博史(1972 卒業論文)
・塑性変形に対する偏差ひずみテンソルの第三不変量とひずみ経路の曲率の影響
山下博史(1974 修士論文)

固体の塑性力学  流体の伝熱工学

3教職員学生と歩んだ伝熱研究(1)<< 伝熱工学１ >> 博士課程・助手の時代 1974年4月～1979年3月
機械工学第二学科 内燃機関及びガスタービン講座
泉亮太郎先生，山口誉起先生，加賀定先生，木下祥次先生
大岩紀生先生，高橋俊夫技官
・直交流形熱交換器の伝熱性能向上に関する研究：山下博史(1978 博士論文)，伊藤浩一(1975卒)
1977年 4月

山下 助手 採用

・対向噴流火炎に関する研究：沢田和徳(1974修)，山田敏雄(1974卒)，伊藤昇平(1975卒)，
飯尾博孝(1977修)，德長幹恵(1979修)，河崎清宣(1981修)，舟橋眞(1982修)
・油冷却器の熱伝達： 花房真(1975修)，浅野一彦(1977修)・冷媒の管内凝縮熱伝達：
小沢通曜(1975修)，小椋健二(1976修)，川戸文郎(1977修)，稲生幸嗣(1978修)・回転平板の物質移動に関する研究： 児玉敏雄(1976修)，原直人(1978修)4

教職員学生と歩んだ伝熱研究(1)・直交流形熱交換器の伝熱性能向上に関する研究

山下博史(1978 博士論文)

熱機器のシステム  熱流体の物理現象（伝熱に関する実験と数値解析）[日本機械学会論文集,
43 (373), (1977), 4230-4236]

5教職員学生と歩んだ伝熱研究(2)<< 伝熱工学２ >> 講師・助教授時代 1979年4月～1988年3月
泉亮太郎先生，加賀定先生，櫛田玄一郎助手，早川正人技官（高浜平七郎先生，藤田秀臣先生，河村鈞先生，藤城明男技官，西村新璽技官）
1981年 4月
1982年 4月

早川正人 技官 着任
櫛田玄一郎 助手 着任・固気混相流による伝熱促進，流動層形熱交換器の研究：
杉山邦生(1978修)，阿部眞一(1981修)，石丸明(1981卒)，
牧田和久(1983修)，梶川吉治(1982卒)，伊藤登喜司(1984修)，
山本敏博(1985修)，松本守聡(1986修)・ふく射フィンの解析： 橋本英樹(1979修)・コンパクト形熱交換器の研究： 野々村浩司(1980修)
・波形流路内の流動および熱伝達：
布施博章(1978卒)，森伸雄(1980修), 平沢一範(1982修)，板倉正人(1981卒)，
宮嶋則義(1982修), 親川兼勇（1982 博士論文）,
・非円形断面直管内における流動および熱伝達に関する研究：青木昇(1983卒)，牧幹人(1983修)
6教職員学生と歩んだ伝熱研究(2)・回転かき取り形凝縮器，直接接触凝縮器に関する研究：
籠橋俊憲(1981卒)，杉浦靖彦(1983修),太田基博(1982卒)，
水野宏幸(1985修)，前田隼（1985 博士論文）
・屈折管・プレートフィンチューブ形熱交換器における流動および熱伝達：水野貴司(1984修)，加藤孝治(1985修)，米谷秀雄(1986修)，
二村仁志(1985卒)，横町尚也(1985卒)，奥村桂三(1986卒)，
芹沢祐司(1986修)，若杉勇(1987修)，櫛田玄一郎（1988 博士論文）
・二次元および軸対称対向噴流における流動および伝熱：花崎裕士(1983卒)，山田浩之(1984卒)，西浦隆幸(1984修)，吉田尊彦(1985卒)，
畔柳佳正(1986卒)，王焱(1986修，1989 博士論文)，太田統之(1988修)，
志満津孝(1989修)・非定常輪送方程式に対する各種差分解の相違に関する研究： 上原健(1988修)・熱サイホンにおける流動および熱伝達： 杉山一弘(1989修)

1985年 3月

泉亮太郎先生 退官

7教職員学生と歩んだ伝熱研究(2)・固気混相流による伝熱促進，流動層形熱交換器の研究[日本機械学会論文集, 49 (440), (1983), 869-876]8教職員学生と歩んだ伝熱研究(2)・波形流路内の流動および熱伝達

親川兼勇(1982 博士論文)[日本機械学会論文集,
48 (435), (1982), 2245-2254]

9

教職員学生と歩んだ伝熱研究(2)・回転かき取り形凝縮器に関する研究[日本機械学会論文集（B 編）, 50 (457),(1984), 2219-2227]

前田 隼(1985 博士論文)

10教職員学生と歩んだ伝熱研究(2)・直接接触凝縮器に関する研究[流体熱工学研究, 23 (1), (1988), 17-29]

11教職員学生と歩んだ伝熱研究(2)・屈折管・プレートフィンチューブ形熱交換器における流動および熱伝達櫛田玄一郎(1988 博士論文)[日本機械学会論文集, 50 (449), (1984), 68-78]

12[日本機械学会論文集(B編）, 51 (461), (1985), 169-176]

教職員学生と歩んだ伝熱研究(2)・二次元および軸対称対向噴流における流動および伝熱王焱(1989 博士論文)[日本機械学会論文集(B 編）, 53 (496), (1987), 3729-3736]

13

教職員学生と歩んだ燃焼研究(1)<< 燃焼工学１ >> 助教授時代 1988年4月～1998年3月
機械工学科 機械エネルギー工学講座  機械情報システム工学科 環境情報システム講座
竹野忠夫先生，櫛田玄一郎講師，西岡牧人助手，早川正人技官

1988年 4月

竹野忠夫先生，西岡牧人助手 着任
伝熱工学に関する数値解析

反応性流体力学に基づく気相燃焼の数値解析

1991年 9月

米国NIST交換研究員(1年間)
Prof. T. Kashiwagi, Prof. H.R. Baum1995年 4月
1998年 4月

朱学雷助手 着任
中村祐二助手 着任

2001年 3月

竹野忠夫教授 退官

14教職員学生と歩んだ燃焼研究(1)●竹野教授，西岡助手

・管状火炎の理論的・解析的研究：
稲垣和久(1990修)，祖父江務(1991修)，西岡牧人(1993 博士論文)，
朱学雷(1995 博士論文)
・火炎中のサーモフレイシスの効果：
内田伸宏(1991修)，田中宏茂(1993修)，湯浅英太郎(1995修)
・乱流予混合火炎の非定常挙動および曲率分布：
馬場直樹(1991修)，井上雄二(1993修)，幾老綾絵(1995修)，
入江徹(1998修)，北村洋介(2000修)，坪根裕二(2001修)・LIF法によるOH濃度の測定
水橋将(1992修)，大坪栄一郎(1994修)，伊藤昭(1996修)，堀江竜太(1998修)，
安藤千博(2000修)
15

教職員学生と歩んだ燃焼研究(1)・火炎中のPAHとすす生成：

近藤健一(1992修)，土屋順久(1994修)，南野圭史(1995修)，田中邦博(1996修)，
平岡直大(1997修)，齋藤健司(1998修)，小森晃(1999修)，酒井伸吾(2000修)，寺部保典(2001修)，橋場敏彦(2001修)，吉田康隆(2001修)，浅野剛生(2002修)・リバーニングによるNOxの還元効果： 加藤滋久(1999修)，中村直樹(2001修)・火炎中のNO生成：
近藤雪水(1992修，1996 博士論文)，森智章(1993修)，栗田敦(1994修)，
竹本安伸(1996修)，江崎秀一郎(1997修)，近藤博是(1999修)，
梅山智貴(2000修)，多武良友孝(2001修)・外部放射熱により加熱される固体燃料の自然着火に関する数値解析：
中村祐二(1995修，1999 博士論文)，小林充(2000修)，片山博之(2000卒)・濃度斑点に起因する予混合火炎の非定常挙動： 和久智裕(1997 博士論文)
・旋回流中の火炎の挙動： 内田正宏(2002修)，寄森一史(2002修)

16

教職員学生と歩んだ燃焼研究(1)・外部放射熱により加熱される固体燃料の自然着火に関する数値解析中村祐二(1995 博士論文)[日本機械学会論文集(B 編), 62 (599), (1996), 2758]

17教職員学生と歩んだ燃焼研究(1)●山下助教授，櫛田講師

・LIF法による軸対称噴流の混合特性に関する計測：
小林正佳(1990修)，町田幸則(1992修)，後藤尚紀(1994修)，赤尾剛(1997修)
中川雅樹(1996卒)，磯部裕之(1999修)

・水平角柱まわりの自然対流熱伝達に関する数値解析： 杉村敏夫(1991修)・二次元噴流拡散火炎の遷移と安定性に関する数値解析：
李倩(1992修)，井戸田敏博(1993修)，児玉明広(1993修)，
島田雅文(1996修)，D. Djamrak(1997修)・低レイノルズ数の微小重力場において輻射加熱された固体からの熱伝達：
H.R.Baum，K.Nakabe，T.Kashiwagi(1993 米国NISTとの共同研究)
・三次元噴流拡散火炎の非定常挙動に関する数値解析：
山崎義弘(1995修)，守口崇之(1998修)・同軸噴流の安定性に関する数値解析： 杉森洋一(1994修)，伊藤英俊(1996修)

18教職員学生と歩んだ燃焼研究(1)・LIF法による軸対称噴流の混合特性に関する計測[Ninth Symposium on Turbulent Shear Flows, 2, (1993), P212_1-P212_4]

19教職員学生と歩んだ燃焼研究(1)・二次元噴流拡散火炎の遷移と安定性に関する数値解析

水素の乱流噴流火炎の実験（シュリーレン法による可視化）に対応する数値シミュレーションから燃焼研究を始めた

y

解析モデル

[Combust. Sci. Technol., 10 (1975), 45]d0Proc. of Roy. Soc. London, A431 (1990), pp.301-314

x20

教職員学生と歩んだ燃焼研究(1)混合分率(結合関数)

温度

Re = 500

Re = 1000
As the Reynolds number was increased above this value, a transition point
appeared along the flame, downstream of which the flame and flow began tofluctuate.

Re = 2000
Proc. of Roy. Soc. London,
A431 (1990), pp.301-314

21

教職員学生と歩んだ燃焼研究(1)
等密度噴流

混合分率 (Re = 1000)

仮想火炎(粘度一定)

凍結噴流

仮想火炎(拡散係数一定)

通常火炎

仮想火炎(粘度・拡散係数一定)

The radial distributions of density and transport coefficients were found to play
dominant roles in this instability, and hence in the transition mechanism.
The decreased density in the flame accelerated the instability, while the increasein viscosity had a stabilizing effect. The most important effect was the increase indiffusion coefficient. The increase shifted the flame surface, where the large
density decrease occurs, outside the shear layer of the jet and produced a thick22
viscous layer surrounding the jet which effectively suppressed the instability.

教職員学生と歩んだ燃焼研究(1)y(mm)[日本機械学会論文集, 65 (630), (1999), 783-789]8温度分布（Re = 1911; ufuel = 30 m/s, uair= 4 m/s, d = 1 mm）
2300

4

2100

0

1900

-4
-8
-2

1700

4.8 11.6 18.4 25.2

32

38.8 45.6 52.4 59.2

x(mm)

66

1500
1300
1100
900
700
500
300

23

温度

y(mm)

教職員学生と歩んだ燃焼研究(1)
8
2300

4

2100

0

1900

NO質量分率y(mm)CH質量分率

y(mm)

-4
-8
8-2

1700

4.8 11.6 18.4 25.2

32

38.8 45.6 52.4 59.2

66

0

1300
0.000004
1100
0.000004
900

-4

700
0.000003

-8
-2
8

500
0.000003
300
0.000002
0.00005
0.000002
0.000045x(mm)44.8 11.6 18.4 25.2

32

38.8 45.6 52.4 59.2

66

x(mm)

4
0

0.000001
0.00004y(mm)-4NO2質量分率

1500

-8
8-2

0.000035
0.000001

4.8 11.6 18.4 25.2

32

38.8 45.6 52.4 59.2

66

x(mm)

4

0.0000300.000025
0.000006
0.0000200.000015
0.000005-40.00001-8-20.000004
0.000005

4.8 11.6 18.4 25.2

32

38.8 45.6 52.4 59.2

x(mm)

66

0
0.000003

24

教職員学生と歩んだ燃焼研究(1)

浮き上がり火炎の可視化実験[25th Symposium (International)
on Combustion, 1994, pp.1175-1181]

数値解析による温度分布（Re = 2000; Da = 1.4×107）[26th Symposium (International)
on Combustion, 26, (1996), 27-34]

25

教職員学生と歩んだ燃焼研究(1)火炎構造の解明のための物理指標

M ixedness :Z FO 

YO
j

 YO

  YF 
 j


Y

 YF  O  YF  jFlame index :GFO  gradYF  gradYO (q  qc )[26th Symposium (International) on Combustion, 26, (1996), 27-34]

26

伝熱・燃焼研究のテーマ(2)・低レイノルズ数の微小重力場において輻射加熱された固体からの熱伝達[Combustion and Flame,
98, (1994), 361-374]

27

教職員学生と歩んだ燃焼研究(2)<< 燃焼工学２ >> 教授時代 1998年4月～2004年3月
機械情報システム工学科 ヒューマンシステム工学講座  環境情報システム講座
趙黛青講師，中村祐二助手，早川正人技官，（新井紀男教授 ）1998年 4月

機械情報システム工学科 ヒューマンシステム工学講座に山下が教授として配置替

2000年 4月

趙黛青講師 着任

2001年 3月

竹野忠夫教授 退官

2001年 4月

機械情報システム工学科 環境情報システム講座に山下が配置替で帰任

気相の噴流拡散火炎の数値解析

気相の予混合火炎・水蒸気添加・消炎，触媒燃焼の数値解析
28教職員学生と歩んだ燃焼研究(2)・Ｖガッタに保持された予混合火炎の火炎構造と非定常挙動：
春日洋祐(1996卒)，勘田将生(1999卒)，福井傑(2001修)・予混合火炎の安定性および非定常挙動に関する数値解析：
阿久澤裕樹(1997修)，大須賀達也(1999修)，池田晃浩(2001修)，竹崎大輔(2000卒)
・乱流拡散火炎に対するLaminar flamelet model に基づくNOx の予測手法の検証：
村松敦(1998修)，寺田勝彦(2000修)，清水昭博(2002修)，後藤英之(2003修)，

高石良伸(2004修)・Triple Flame の火炎構造に与える燃料および当量比の影響：
堤谷進也(1998修)，崔洛挺(1999 研究員)，吉見政史(2000修)，多田章彦(2002修)・リバーニング特性に関する数値解析： 森田智彦(2001修)・横風を受ける燃料噴流火炎に関する数値解析：
林秀幸(2002修)，大塚裕史(2004修)

29教職員学生と歩んだ燃焼研究(2)・酸素富化空気燃焼における火炎構造およびNO生成特性の解明：
岸本衛(2002修)，池田光芳(2004修)，斎木悠(2004卒)
・旋回流中での予混合火炎の伝播に関する数値解析：
上田大輔(2002修)，永井英和(2005修)，近藤周司(2008修)，
篠田昌久(2008 研究員)
・Rich-Lean 二段燃焼システムにおけるNOx 生成に関する数値解析： 趙黛青(2000)
・燃料過濃予混合火炎中のNOx 生成および水蒸気添加の影響： 趙黛青(2001)
・メタン-空気拡散火炎の燃焼限界に与える水蒸気添加の影響： 喜多敦(2003修)

・水素-酸素拡散火炎に与える磁場の影響に関する数値解析： 山田英助(2002修)・固体壁の小円孔を通過する予混合火炎の消炎に関する数値解析：
藤田英之(2003修)
・分光法による拡散火炎中の煤放射特性の解明：
鎌田祐一(2003 研究員)，高橋真一(2005修)30

教職員学生と歩んだ燃焼研究(2)・メタンハイドレートの熱分解に関する数値解析： 恵藤陽介(2005修)・マイクロフレームの安定性と火炎構造に関する数値解析： 窪田充志(2005修)・アセトン-OH同時LIF 法による乱流予混合火炎の火炎構造の測定：
山本典史(2003修)，里見知彦(2004修)，林直樹(2004修)，馬目聡(2006修)・拡散・部分予混合火炎中のPAHの生成特性に関する実験的検討：
山口智丈(2003修)，長谷川義朗(2004修)，石井大祐(2006修)
・水素あるいはメタン-空気予混合気の流路内触媒燃焼に関する数値解析：
佐多宏太(2003修)，太田稔(2005修)，松永秀一(2005 共同研究)，
高下峻一(2007修)・メタン-空気混合気の三次元不均一濃度温度場での燃焼過程の数値解析：

山田真人(2004修)，椴山由貴(2007修)31

教職員学生と歩んだ燃焼研究(2)・Ｖガッタに保持された予混合火炎の火炎構造と非定常挙動：[日本機械学会論文集B 編, 63 (610), (1997), 2201]

32

教職員学生と歩んだ燃焼研究(2)[日本機械学会論文集B 編, 63 (610), (1997), 2201]

33

教職員学生と歩んだ燃焼研究(2)・乱流拡散火炎に対するLaminar flamelet model に基づくNOx の予測手法の検証

三次元複雑形状，乱流・非定常状態，多様な燃料★ 燃焼による発熱が流れに及ぼす効果

簡単な化学反応機構★ NOxやすす等の有害大気汚染物質の生成機構

詳細な化学反応機構
燃焼場に適用できる乱流反応モデルが必要
現象論的モデル  科学的モデル

層流火炎片モデルに基づく組み合わせ手法[燃焼研究, 124, (2001), 47-55]

34

教職員学生と歩んだ燃焼研究(2)[燃焼研究, 124, (2001), 47-55]

35

教職員学生と歩んだ燃焼研究(2)・Triple Flame の火炎構造に与える燃料および当量比の影響
Side

y
[mm]

Y
u
v
T
 0,
 0,
 0, i  0
y
y
y
y

Air
l

u0 = um(1+Asin(2πft))

混合層0v = 0 , T = 300 K

Fuel[日本機械学会論文集, 65 (630), (1999), 775-782][第51 回伝熱シンポジウム講演論文集, (2014), CD-ROM]

Exit
u
 0,
x
l
v
 0,
x
T
 0,
x
0
1
Yi
0
x

x

燃料濃度

x [mm]

u0 ：噴出速度 [ m/s ]um ：平均速度 [ m/s ]A ：振幅
f ：周波数 [ Hz ]l ：混合層幅 [ mm ]36

教職員学生と歩んだ燃焼研究(2)熱発生速度分布

計算条件
um = 0.7 m/s
f = 100 Hz
A = 0.8
l = 4 mm

非定常下における火炎伝播特性
非定常性の影響（速度場に対する燃焼場の遅れ）
流体速度と燃焼速度の関係
火炎強度に与える諸因子の影響
火炎伸長率（速度勾配・火炎曲率）・熱損失・物質拡散[日本機械学会論文集, 65 (630), (1999), 775-782][第51 回伝熱シンポジウム講演論文集, (2014), CD-ROM]37

教職員学生と歩んだ燃焼研究(2)計算条件
um = 0.7 m/s
f = 100 Hz
A = 0.8

Mixednessll[-]Flame indexlll = 0.5 mm[ m-2 ]

l = 4.0 mm[日本機械学会論文集, 65 (630), (1999), 775-782] ,[第51 回伝熱シンポジウム講演論文集, (2014), CD-ROM]38

教職員学生と歩んだ燃焼研究(2)・酸素富化空気燃焼における火炎構造およびNOX 生成特性の解明[日本機械学会論文集, 68 (669), (2002), 1578-1585]

39

教職員学生と歩んだ燃焼研究(2)・旋回流中での予混合火炎の伝播に関する数値解析

z [mm]

外側境界での圧力を固定（ゼロに設定）

自由流出境界

32

断
熱
す
べ
り
境
界
流入温度
当量比
層流燃焼速度
伝播速度（=流入速度）

8

300 K
1.0
0.44 m/s
22.7 m/s

0

Uin

Vθ

32

r [mm]

100

解析モデルおよび境界条件40

教職員学生と歩んだ燃焼研究(2)

P
32

［Pa］

24

100
0
-100
-200
-300
-400

16
8
0

0

8

16

24

T

z[mm]

32

32

［K］

24

2800
2300
1800
1300
800
300

16
8
0

0

8

16

24

32 ｒ[mm]

Vθ

32

［m/s］
20
16
12
8
4
0

24
16
[日本機械学会論文集, 67 (662), (2001), 2567-2573]

8[日本機械学会論文集(B 編）, 71 (705), (2005), 1483-1489][日本機械学会論文集（B 編）, 71 (709), (2005), 2339-2344]

0

0

8

16

24

32

41

教職員学生と歩んだ燃焼研究(2)・燃料過濃予混合火炎中のNOx 生成および水蒸気添加の影響[化学工学論文集, 27 (6), (2001), 779-785]

42

教職員学生と歩んだ燃焼研究(2)[化学工学論文集,
27 (6), (2001), 779-785]

43

教職員学生と歩んだ燃焼研究(2)・メタン-空気拡散火炎の燃焼限界に与える水蒸気添加の影響[化学工学論文集, 30 (6),(2004), 808-814]

44教職員学生と歩んだ燃焼研究(2)・水素-酸素拡散火炎に与える磁場の影響に関する数値解析[Combustion Science and Technology,
174 (9), (2002), 131-146]

45

教職員学生と歩んだ燃焼研究(2)・固体壁の小円孔を通過する予混合火炎の消炎に関する数値解析[日本機械学会論文集(B 編), 69 (687), (2003), 2523-2530]

46教職員学生と歩んだ燃焼研究(2)[日本機械学会論文集(B 編), 69 (687), (2003), 2523-2530]

47

教職員学生と歩んだ燃焼研究(2)・マイクロフレームの安定性と火炎構造に関する数値解析[Combustion Theory and Modeling, 10 (6), (2006), 927-938]48[日本機械学会論文集，B 編, 74 (739), (2008), 699-706]

教職員学生と歩んだ燃焼研究(2)・アセトン-OH同時LIF 法による乱流予混合火炎の火炎構造の測定
YAG Laser

Dye Laser

283.222nm
OH(1,0)Q1(7)Pulse Generator
PC
SHG
Spherical Lens

ICCD Camera

Harmonic
Separator

アセトンおよびOHの蛍光を同時に計測可能

Cylindrical Lens

メタン･ 空気予混合火炎

Bunsen Burner

予混合気：CH4+Air+Acetone
流量：1000cm3/s 平均流速：1.75 m/s
混合当量比φ：0.66～1.06[日本燃焼学会誌, 47 (142), (2005), 290-299]

49

教職員学生と歩んだ燃焼研究(2)
アセトン

T[K]

既燃側

2100
1800
1500
1200
900
600
300
0

未燃側

1.0
Acetone/Acetonemax

0.8
0.6
0.4

OH/OHmax

0.2

T

0.0
16

15

14

13

12

11

10

[mm]

9

8

7

Intensity(OH/OHmax),
Intensity(Acetone/Acetonemax)

OH

予混合火炎の火炎帯領域の構造

650

教職員学生と歩んだ燃焼研究(2)

(a)

(b)
unburned
arm

OH

(c)

(d)

unburnedislandburnedislandacetone[side part]

(e)

(f)

(g)

(h)
OH spot
dark spot[top part][Journal of Visualization, 11 (1), (2008), 71-78]

51

教職員学生と歩んだ燃焼研究(2)・拡散・部分予混合火炎中のPAHの生成特性に関する実験的検討[日本機械学会論文集B 編, 72 (724), (2006), 358-366]

52

教職員学生と歩んだ燃焼研究(2)・水素あるいはメタン-空気予混合気の流路内触媒燃焼に関する数値解析

解析モデル

ハニカム内の一つの流路を考慮する

平行平板間流路で近似する
Pd

計算領域

53

教職員学生と歩んだ燃焼研究(2)ySymmetry plane： ∂F /∂y = 0(F: vx,Yi ,T) , vy = 0

l=8mm

Exit ：

CH4 + air∂F /∂x = 0(F: vx,vy,Yi,T)
d=1mm

x
Inlet ：
Methane-air (φ)
U = Uin, T = Tin

Catalitic
wall ： Y 
 Di  i   wi , s y  si = H2, H, O2, O, H2O, CH4, CO, OH, CO2

vx = vy= 0

CH4

 Yi   0ysOther species

T = Ts

CO2
H2O
CH4

54

教職員学生と歩んだ燃焼研究(2)触媒反応機構（吸着→表面反応→脱離）
CH4
A(H)

CH3 (s)

R(8)

104

H

小

103

102

D(H2)

R(7)

CH2 (s)

反応速度

大A(CH4)

反応速度の大きさを評価

A(H2)

D(H2O)H(s)R(3)

CH (s)

O2

H2O

H2

A(H2O)

H2O(s)

D(O2)

D(OH)

OH

R(9)

A(OH)

OH (s)

R(5)
R(1)

C (s)
R(10)

O(s)

A(CO)

CO (s)
D(CO2)

CO2

D(CO)

CO
Surface
Chemical species

O

[1/s]

A(O)

Gas phase
Chemical species

A(O2)

計算条件

当量比(φ)：0.044流入速度(Uin)：3.5m/s
流入温度(Tin)：623K
x ：1.0mm55

教職員学生と歩んだ燃焼研究(2)[日本燃焼学会誌, 47 (140), (2005),137-145]

56

教職員学生と歩んだ燃焼研究(3)<< 燃焼工学３ >> 教授時代 2004年4月～現在
機械理工学専攻 機械科学分野 環境・エネルギー工学講座 伝熱・燃焼工学研究グループ
山本和弘准教授，林直樹助教，早川正人技術職員
2002年 4月
2004年 4月
2005年 4月
2005年12月
2011年 4月

山本和弘助教授 着任（後に准教授）
機械理工学専攻 機械科学分野 環境・エネルギー工学講座
伝熱・燃焼工学研究グループ に改組
篠田昌久研究員 採用（2007年10月まで）
林直樹助手 採用（後に助教）
韓基琳研究員 採用（2014年3月まで）
反応性流体力学に基づく気相・液相・表面燃焼の数値解析（噴霧燃焼，表面反応，超小型燃焼器）
輻射平衡の理論解析

2015年 3月

山下 定年退職

57

教職員学生と歩んだ燃焼研究(3)●山本准教授・格子ボルツマン法による転炉内２次燃焼の解析：古賀輝久(2006修)・DPFの格子ボルツマン法による数値解析：
越智文洋(2005修)，佐竹真吾(2007修)，大堀晋也(2008修)，松井健太(2009修)中村匡統(2010修)，山内一樹(2011修)，酒井達也(2012修)，小見山亮(2013修)，
戸田祐輔(2014修)，庭田修平(2015修)・DPFの捕集と再生性能に関する実験的検討：
近藤真一(2012修)，鈴木一也(2013修)，加藤寛也(2015修)・NO₂ によるすす酸化促進効果の実験的検討： 金森陽平(2010修)・セルオートマンによる火災時の避難行動のシミュレーション：
小久保聡(2008修)，左高秀多(2011修)，竹内裕貴(2014修)
・対向流予混合火炎中のOH濃度と燃焼速度： 伊佐山勉(2007卒)
58

教職員学生と歩んだ燃焼研究(3)・OH-HCHO同時PLIF法による乱流予混合火炎の可視化と火炎構造：
井上聡(2005修)，尾関賢宏(2006修)，水谷浩幸(2005卒)，
大西將博(2007修)，石井慎治(2009修)・部分予混合雰囲気中の火炎の燃え拡がり解析： 緒方佳典(2007修)・旋回噴流燃焼器の燃焼特性と火炎構造の検討：

金川裕司(2007修)，鈴木啓夢(2013修)・LII 法によるすす計測とディーゼル排気ガスへの適用： 藤掛文裕(2008修)・三重管バーナにおける燃焼排出物および浮上り火炎の挙動に関する研究：大島久治(2008修)，竹本昌弘(2009修)，磯部佑介(2010修)，加藤慎也(2010修)・PIV/OH-PLIF同時計測による乱流予混合火炎の検討：
麻田泰生(2011修)，伊藤雄貴(2012修)，奥洋平(2014修)，森岡諒(2015修)・純酸素を用いた燃焼現象のモデル化とバーナの開発： 奥山悟郎(2010 博士論文)
・反応焼結SiCを用いたDPFとディーゼル後処理過程の研究：
常吉孝治(2012 博士論文)

59教職員学生と歩んだ燃焼研究(3)●山下教授，林助教

・二次元非定常予混合火炎に関する詳細素反応機構による数値解析：
加藤敏宏(2005修)・二段階吸引式バーナを用いた焼成炉の数値解析：
企業との共同研究，池ケ谷智章(2006修)・メタンおよび水素-空気予混合気の着火特性に関する数値解析：
松山竜佐(2006修)，韓基琳(2010 博士論文)

・流路内触媒反応に関する素反応機構を用いた数値解析：
相川司(2007修)，家根弘好(2009修)・非定常対向流予混合火炎の火炎構造に関する数値解析： 野々村一樹(2008修)・超小型燃焼器の開発に関する数値解析：
鈴木新(2008修)，加藤宏紀(2010修)，上原翔(2011卒)，加藤雄大(2012修)，
白崎洋平(2013修)，白明蕊(2013修)，木野雄介(2014修)

60

教職員学生と歩んだ燃焼研究(3)・センターエアーバーナによる燃焼器の小型化に関する数値解析：
槙田雄一郎(2009修)・曲率を持つ水素-空気予混合火炎の火炎構造に対する支配パラメータの解明：
宮田智生(2008卒)，荒井康(2013修)・閉空間におけるふく射平衡に関する理論的考察および数値解析：

近藤良夫(2010 博士論文)
・超小型燃焼器における壁面でのラジカルクエンチングを考慮した数値解析：
今井俊幾(2010修)，高鳥隼(2015修)
・先進的ディーゼル燃焼に関する革新的コンセプトおよびシミュレータの開発：
稲垣和久(2012 学位論文)
・噴霧燃焼の燃焼形態に関する数値解析：
萩原康太(2009修)，宇野将成(2011修)，山倉真二郎(2012卒)，
近藤聡介(2013修)，丹羽雄基(2014修)，小口和麻(2015修)
61教職員学生と歩んだ燃焼研究(3)・噴霧燃焼の火炎伝播特性に関する数値解析：
李加寧(2012修)，棚瀬利勝(2014修)・過濃予混合気の高温空気燃焼における燃焼特性および非定常挙動の数値解析：
伊藤慎太朗(2011修)，近藤俊樹(2012卒)，木村旭宏(2015修)・予混合火炎の燃焼特性に与える火炎曲率・壁面の影響に関する実験的検討：
鈴木裕介(2010修)，竹内啓介(2012修)，川本航平(2014修)・拡散火炎におけるすす生成特性の検討と組合せ予測手法の検証：
森上拓也(2012修)

・予混合火炎の燃焼特性に与える不活性ガス種の影響に関する検討：
冨田佑平(2013修)，橋岡司(2015修)・両側の空気流速が異なるスリット火炎の浮き上がり特性に関する数値解析：

安田翔(2014修)・燃料希薄メタン-空気予混合火炎の火炎構造と支配パラメータの解明：
林直樹(2015 博士論文 受理)

62

教職員学生と歩んだ燃焼研究(3)・二次元非定常予混合火炎に関する詳細素反応機構による数値解析
20

自由流出境界

解析モデル

流入境界

Premixed gas
CH4＋Air
T0  300K

自由流出境界当量比φ= 0.75
-2

18

対称境界

 u m： 平均流入速度《 流入速度 》

 Amp： 変動の強さ

u0  um  1  Amp  sin2 k y  sin2 f t 

 k : 波数
 f ： 周波数 t ： 時間63

教職員学生と歩んだ燃焼研究(3)火炎の曲率や凸となる向きだけでなく，火炎の非定常的な動きによって、既燃ガス中のOH濃度が大きく変化している．

ある瞬間の現象だけでなく、時間的な変化を検討する必要がある！

★ 火炎面形状・火炎面曲率の大きさ
・火炎面曲率の正負（火炎面が既燃側に凸か凹か）★ 火炎面の伝播方向
um =2.0 m/s , k = 0.3 mm-1 , f = 500 Hz[日本機械学会論文集, 71 (712), (2005), 3038-3045]

64教職員学生と歩んだ燃焼研究(3)・二段階吸引式バーナを用いた焼成炉の数値解析（企業との共同研究）

Recirculation pipe

Rectifier pipe

解析モデル

Burned gas

速度分布

温度分布
65

教職員学生と歩んだ燃焼研究(3)・メタンおよび水素-空気予混合気の着火特性に関する数値解析
Ts = 1000 K, q = 60 GW/m3

0.5

MIE (mJ)

0.45
0.4

R=0.3 mm0.35

R=0.8 mm
R=1.5 mm

0.3
0.25
0.2
0.15
0.1
0.05
0
0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

1.1

1.2

gap distance (mm)[Combustion and Flame, 157 (7), (2010), 1414-1421][International Journal of Hydrogen Energy, 36 (15), (2011), 9286-9297]

韓 基琳(2010 博士論文)

66教職員学生と歩んだ燃焼研究(3)・流路内触媒反応に関する素反応機構を用いた数値解析

Inlet:
Exhaust Gas
vx = Uin’vy=0
YNO,in’
T = Tin’
Yi =Y i,in’Fig.1 Schematic of analytical model
NOの酸化： 白金Pt触媒[日本機械学会論文集（Ｂ編）, 75 (756), (2009), 1669-1677]

NOの還元： ロジウムRh触媒[日本燃焼学会誌, 49 (148), (2007), 128-136]

67

教職員学生と歩んだ燃焼研究(3)・非定常対向流予混合火炎の火炎構造に関する数値解析[日本機械学会論文集（B 編）, 74 (747), (2008), 2393-2400]

68

教職員学生と歩んだ燃焼研究(3)[日本機械学会論文集（B 編）, 74 (747), (2008), 2393-2400]

69

教職員学生と歩んだ燃焼研究(3)・超小型燃焼器の開発に関する数値解析

戦略的国際科学技術協力推進事業 日本－中国 研究交流
研究課題：
「革新的な超小型燃焼器を用いたマイクロパワーシステムに関する基礎研究」
研究交流期間 平成19年12月～平成22年11月

研究代表者： 山下 博史（名古屋大学大学院工学研究科、教授）
趙 黛青（中国科学院広州エネルギー研究所、副所長）
小型・高性能の携帯機器，ロボット等の開発が進み，高いエネルギー密度のエネルギー源が必要

超小型化の問題点

燃焼効率・燃焼強度・火炎安定性の低下・消炎の発生

70

教職員学生と歩んだ燃焼研究(3)多孔体壁を有する超小型燃焼器の提案
予混合気

メタン燃料

[日本燃焼学会誌, 55 (172),(2013), 131-137]

多孔体流入壁（焼結銅）

二次空気

燃焼器の壁面の大部分を多孔体壁とし，予混合気を流入させる 
壁面からの放熱を流入予混合気で回収（熱循環）71教職員学生と歩んだ燃焼研究(3)[日本燃焼学会誌, 55 (172), (2013), 131-137]

解析モデル

Cubic of 1 mm side length
0.25mm
0.75mmPorous Medium Injector
1.0mm

H2/Airpremixed gas
ZYX

Burnt gas
Solid

1.0mm

1.0mm
Flame propagation Ignited here
161×81×81 Equally-spaced（0.0125 mm）

72

教職員学生と歩んだ燃焼研究(3)[日本燃焼学会誌, 55 (172), (2013), 131-137]

不規則配置における固体球配置

yy

x x
z
z

73

教職員学生と歩んだ燃焼研究(3)[日本燃焼学会誌, 55 (172), (2013), 131-137]

porous medium

H2/Airpremixed
gas

Flashback occurs for larger porosity

in Section z = 0.5 mm

Burnt
gas

Temperature [K]500  540

Flame stands injector downstream edge at first.
Flame goes into porous medium injector gradually.

74

教職員学生と歩んだ燃焼研究(3)x = 1.06875 mm (I=87)

Temperature [K]x =1.11875 mm (I=91)[日本燃焼学会誌, 55 (172),(2013), 131-137]

Temperature is uniform in each microsphere, and increases gradually due toheat conduction from flame arrested at porous medium downstream edge 75

教職員学生と歩んだ燃焼研究(3)[日本燃焼学会誌, 55 (172), (2013), 131-137]

porous medium

H2/Airpremixed
gas

Burnt
gas

in Section z = 0.5 mm

H2O mass production rate[kg/(m3･s)]

水の質量生成速度と速度ベクトル分布

76
Random model, ε=0.5591，Rh =0.03314 mm，Uin =1.3 m/s, Tin =500 K, Yi =Yi,in (φ＝ 0.4)

教職員学生と歩んだ燃焼研究(3)・センターエアーバーナによる燃焼器の小型化に関する数値解析[日本機械学会論文集（B 編）, 75 (760), (2009), 2535-2543]

Tair1 (K)

Uair1 (m/s)φ

条件A

293.15

0.500

5.6

条件B

1500

2.558

5.6

条件C

1500

20.4670.7

77

教職員学生と歩んだ燃焼研究(3)・閉空間におけるふく射平衡に関する理論的考察および数値解析

A : 空間の性質

近藤良夫(2010 博士論文)

閉空間において

内壁温度均一⇔ふく射平衡
ふく射平衡の状態は内壁温度のみに依存（ふく射平衡時に内壁射度が均一 (EB) かつ放射率に依存しない）⇔

G. Kirchhoff
B : 壁面（物質）の性質

本研究
完全拡散反射・
放射面の場合任意温度の放射面において

放射率ε＝吸収率α
面自身の放射エネルギー＝εEB （T,λ）●A,Bの事象が同値であるという物理現象の見方を提案[日本伝熱学会論文集, 19 (1), (2011), 17-24]

78教職員学生と歩んだ燃焼研究(3)・超小型燃焼器における壁面でのラジカルクエンチングを考慮した数値解析

壁面におけるラジカルクエンチングの解明
表面反応の適用
発熱速度が低下し，燃焼反応が抑制される

熱損失
高温

H
低温

放熱

H

ラジカルクエンチング
H

活性

これまでの燃焼器の開発では無視

H2
不活性

H

79

教職員学生と歩んだ燃焼研究(3)
O2

表面反応の
反応経路図

O

H

Wd 2
Wa 2

O*

H2O

Wd 3

Wa1
Wa 4

HO2

OH*

Wa 3

Wa 5
Gas phase species

H*
OH

Wd 1

H2

H2O2

Radical species
Adsorbed species : *Desorption

Adsorption

ラジカルクエンチング機構:
Adsorption of radical species
Surface reaction of adsorbed species
Desorption of stable species

80教職員学生と歩んだ燃焼研究(3)[9th Asia-Pacific Conference on Combustion, (2013), CD-ROM][Proceedings of Power MEMS, (2014), 1]

Analytical model and boundary conditions
[mm]

[mm]

0

0

Inert wall

Active wall

1.0
1.2

1.0
1.2
0

[Proceedings of the Combustion
Institute, 34 (2), (2013), 3379-3386]

6.5

[mm]

0

6.5

[kg/(m3・s)][mm]Distributions of OH mass production rate wOH
81
for Tw = 1200 K，yw = 1.2 mm

教職員学生と歩んだ燃焼研究(3)・噴霧燃焼の燃焼形態に関する数値解析

反応性流体力学の基礎方程式（PSI-CELLモデルを適用）
連続方程式，運動方程式，エネルギー方程式，
各化学種の連続方程式，状態方程式

   ( v )  Z Ft
 ( u m )P 0 u j  u j
   ( vu m )    ( u m )  
 (     ) g m  ()  Z M ,mtxmxm x j x j xm

 ( Yi )   ( vYi )    ( Di Yi )  wi   Fi Z Ft

i  F (CH 4), O2, N 2, CO2, H 2OqZ
 ( T )1
1 Dp 1   ( vT )    (T )   hi wi   (c p ,iYiVi  T )  rad  Ttcpc p Dt c p i
cp i
cp
cp

82

教職員学生と歩んだ燃焼研究(3)解析モデル
自由流出境界4mm

燃料液滴
蒸発前：軽油（蒸発後：CH4 ）
初期直径
112mm
d1 =10 μm ，d2=20 μm， d4=40 μm
平均数密度 n1, n2，n4
温度： 500 K（沸点）
112mm

112mm

z

差分格子
格子数 100×100×100
中心部60格子 等間隔格子(0.1 mm)
外周部40格子 不等間隔格子(0.1~10.8 mm)

液滴配置領域体積 Vyx[日本燃焼学会誌, 51 (158), (2009), 343-353]

酸化剤
空気
初期温度： 1000 K
初期当量比φ初期時刻に，液滴配置領域体積において，
全燃料液滴質量全酸素質量φ 
量論の場合の燃料・酸素質量比83教職員学生と歩んだ燃焼研究(3)
A

B

C

D

燃焼形態の詳細な分類 ［萩原ら，2009]
A
モードⅠ：単滴着火・燃焼
モードⅡ：一部連結着火・燃焼
モードⅢ：内部着火・燃焼
モードⅣ：内部着火かつ外周部着火・燃焼
モードⅤ：外周部着火・燃焼
C
モードⅥ：外殻燃焼
D

Chiuらの噴霧燃焼形態
Ⅰ

Ⅱ

小

Ⅲ

Ⅳ

初期当量比 φ[Chiu, H.H. and Liu, T.M., Combust. Sci. Tech., Vol.17, 127-142, (1977)]

Ⅴ

Ⅵ

大
84

教職員学生と歩んだ燃焼研究(3)

液滴群の三次元分布の時間変化
t = 40-55 μs(熱発生速度の等値面)85教職員学生と歩んだ燃焼研究(3)10

噴霧燃焼ダイヤグラム

5

 =10

104
2

G =10

n- = N/V [mm -3]

10[日本燃焼学会誌, 51 (158),(2009), 343-353]

3

C

D

 =1

Ⅰ (A)
Ⅱ
Ⅲ
Ⅳ
Ⅴ (C)
Ⅵ (D)

A:単滴燃焼G < 10-2

 =0.1102

B:内部群燃焼10-2 < G < 1

G =1
101

C:外部群燃焼
1 < G < 102

B100

D:外殻燃焼102 < G

-2

G =10
10-1

A
10-2 010101d [m]

102

86

教職員学生と歩んだ燃焼研究(3)Chiuらの「内部群燃焼」が再現できていない

熱発生速度と速度ベクトル

計算領域の拡大
デカン燃料
16000個, 20 μm
T0=1100 K

87

教職員学生と歩んだ燃焼研究(3)・噴霧燃焼の火炎伝播特性に関する数値解析

解析モデル
観測領域

y

液滴列燃焼の燃え拡がり過程におけるシュリーレン画像 （小林秀昭）

8mm

8mm

8mm

1 mm

C10H22 +Air

x
z
155mm
境界条件
入口境界（x = 0）
噴射口： 燃料噴霧と周囲空気を初速度U0で噴射
噴射口以外： すべりなし壁面
その他の境界： 自由流出

88

教職員学生と歩んだ燃焼研究(3)解析モデルおよび境界条件

初期液滴直径 d，平均数密度 n

初期当量比φ0

Spray region

Side

4 mm
Exit(自由流出)

Symmetry
plane
4 mm

Z

Y
X

Ignited here（1500K）

Flame propagation

格子数 121×81×81
最小格子間隔 0.1 mm

燃料液滴
燃料：軽油，デカン
温度： 沸点周囲流体
空気
初期温度： 液滴と同じ

89

教職員学生と歩んだ燃焼研究(3)Droplet diameter: 40 μm
Number density: 4 mm-3
(Equivalence ratio: 2.82)

T
日本燃焼学会誌第54巻
目次裏頁（2012年度）

90 90

教職員学生と歩んだ燃焼研究(3)Droplet diameter: 20 μm
Number density: 16 mm-3
(Equivalence ratio: 1.41)

TT
日本燃焼学会誌第54巻
目次裏頁（2012年度）

91教職員学生と歩んだ燃焼研究(3)Flame propagation velocity [m/s]火炎伝播速度と初期当量比の関係(燃料噴霧燃焼)

d0 [m]
10
20
30
40
50
60
70

10

5

0

0.05 0.1

0.51

Initial equivalence ratio φ0

5

10
92

教職員学生と歩んだ燃焼研究(3)・過濃予混合気の高温空気燃焼における燃焼特性および非定常挙動の数値解析

Flame Index： 総括反応的な火炎構造
 酸素の識別を含む素反応機構による火炎構造の解明
過濃予混合気を周囲空気流中に噴出するような燃焼場：
酸化剤としての酸素が予混合気と周囲空気の両方に存在
この酸素分子O2 を識別  X2 ，Y2
種々の化学種が酸素原子X，Yのどちらで構成されているかで予混合火
炎構造か，拡散火炎構造かが解明できる．
・ 反応速度定数は同じとし， 「実現割合」 という概念を導入
数学的（順列・組み合わせ） ＋ 化学反応論的・ 輸送係数，熱力学データなどは全て一緒
93

教職員学生と歩んだ燃焼研究(3)
実現割合

反応例）
k = 21HO2 + HO2 ---> H2O2 + O2

WO2 = k21・CHO22

k = 62:k = 63:k = 64:k = 65:k = 66:k = 67:
k = 68:k = 69:k = 70:k = 71:k = 72:k = 73:

HX2 + HX2 --->
HX2 + HXY --->
HX2 + HXY --->
HX2 + HY2 --->
HX2 + HY2 --->
HX2 + HY2 --->
HY2 + HXY --->
HY2 + HXY --->
HY2 + HY2 --->
HXY + HXY --->
HXY + HXY --->
HXY + HXY --->

W62 = 1/1・k21・CHX22
W63 = 2/2・k21 ・CHX2・CHXY
W64 = 2/2・k21 ・CHX2・CHXY
W65 = 2/3 ・k21 ・CHX2・CHY2
W66 = 2/3 ・k21 ・CHX2・CHY2
W67 = 2/3 ・k21 ・CHX2・CHY2
W68 = 2/2・k21 ・CHY2・CHXY
W69 = 2/2・k21 ・CHY2・CHXY
W70 = 1/1・k21 ・CHY22
W71 = 1/3 ・k21 ・CHXY2
W72 = 1/3 ・k21 ・CHXY2
W73 = 1/3 ・k21 ・CHXY2

H2X2 + X2
H2X2 + XY
H2XY + X2
H2X2 + Y2
H2XY + XY
H2Y2 + X2
H2XY + Y2
H2Y2 + XY
H2Y2 + Y2
H2X2 + Y2
H2XY + XY
H2Y2 + X2

3/3
0/3
3/3

0/3
3/3
0/3

WSUM = k21 ・（CHX22 + 2・CHX2・CHXY + 3・2/3・CHX2・CHY2 + 2・CHY2・CHXY + CHY22 + 3・1/3・CHXY2）
CHX22 + CHXY2 + CHY22 + 2・CHX2・CHXY + 2・CHX2・CHY2 + 2・CHY2・CHXY=（ CHX2 + CHXY + CHY2）2 = CHO22

94

強

100

5

0

4

-100

拡散火炎

-200

3
wX2
wY2
wXY
wX2+wY2 + wXY
Q

予混合火炎
-300

2

Q [GW/m3]

wX2 , wY2 , wXY [kg/(s・m3)]

反応の強さ

弱

教職員学生と歩んだ燃焼研究(3)

1

-400
17

18

19

0
20

x [mm]

過濃予混合気側[日本燃焼学会誌, 53 (166), (2011), 271-283]

周囲空気側

95

教職員学生と歩んだ燃焼研究(3)・曲率を持つ水素-空気予混合火炎の火炎構造に対する支配パラメータの解明
wH2O,max[kg/(m3・s )]y [mm]
Exit∇n,qT [K/mm]

14Slip wall

N2
H2-Air

7
5

0.28
-0.28

N2

Flame
Stagnation
plane

ur
u0
10

x [mm]-5-7

Slip wallFig.1 Analytical model of counterflow
premixed-twin flame

-14

Exit

∇n,qCf [mm-1]

Fig.8 Relationship of wH2O,max to ∇n,qCf and ∇n,qT[日本機械学会論文集(B 編), 79 (807), (2013), 2486-2496]

96教職員学生と歩んだ燃焼研究(3)解析モデル
Methane-air
premixed gas

High temperature air

u0, 0, T0

0.25d

d

3d

uair, Tair

d
z

3d

4d

y
The origin is thex center of left-hand side

[第52 回燃焼シンポジウム講演論文集, (2014), 1]

Calculation grids
Number: 81×61×61
Spacing: d/20

97

教職員学生と歩んだ燃焼研究(3)Distributions of Q [GW/m3] (z = 0 mm section)

0  1.2, T0 = 300 K, Tair = 1500 K, d = 1.0 mm

y [mm]
x [mm]
1.5

0

4.0

-1.5

0

4.70

(a) u0 = 3 m/s

0

4.75

(b) u0 = 5 m/s

04.91(c) u0 = 7 m/s[第52 回燃焼シンポジウム講演論文集, (2014), 1]○ Flame with curvature is formed on central axis.
○ With increase in u0, the flame tip moves to air side, the heat release rate
becomes small, and then the flame is extinct on central axis.
In the following, we examine the behavior of various quantities on central axis.

98教職員学生と歩んだ燃焼研究(3)回燃焼シンポジウム講演論文集,
2S [第52(2014), 1] 1  a1   nn：vflame  ・vflame
R a1 : Local velocity gradient1 ・vfluid,t  vflame, n(・n)

S1 : Local burning velocity
R : Radius of flame curvature 2  ・vfluid,tv [m/s]Tangential gradient of1
tangential
velocity atflame tip

0

-1
-0.5

Q[GW/m3]0y [mm]

0.5

y direction distribution of tangential
99
velocity for flame surface

教職員学生と歩んだ燃焼研究(3)Qsum×104 [GW/m2]

108d [mm]
1.0
0.8
0.6
0.4
0.3

6
4

0.9
0.7
0.5
0.35
0.25

2
0

0

10
1×10-3 [s-1]20○ κ2 is almost proportional to u0 , and its difference due to the scale of combustion field is small.
100○ The correlation is seen in Qsum and κ2 for different d near extinction.

教職員学生と歩んだ燃焼研究(3)・予混合火炎の燃焼特性に与える火炎曲率・壁面の影響に関する実験的検討

(a)

(b)

(c)

(IOH*IHCHO)max

4000
3000
2000
1000
0
1.0

2.0
3.0
4.0
-1
Curvature [mm ]5.0

Fig.1 Pictures of PLIF
Fig.3 Relation between (IOH*IHCHO)max and(a)OH, (b)HCHO, and (c)OH*HCHO
curvature of methane-air flame
Work-in-progress for 33th International
Symposium on Combustion (2010)101教職員学生と歩んだ燃焼研究(3)・拡散火炎におけるすす生成特性の検討と組合せ予測手法の検証

6.0
XC2H4,F
1.0
0.8
0.6
0.4

fV×10-13 [-]

5.0
4.0
3.0
2.0
1.0
0
6

6.5

7
x [mm]

7.5

8

Fig. 1 Distributions of soot volume fraction[第49回燃焼シンポジウム前刷集 178-179 2011]

Fig. 2 Relationship of ωC16H10,sum to Tmax and
SDRq
102

教職員学生と歩んだ燃焼研究(3)・予混合火炎の燃焼特性に与える不活性ガス種の影響に関する検討

Fig.3 OH-PLIF system[第５０回燃焼シンポジウム 158-159 2012](a) N2, q = 0.81, vm = 0.8 m/s(b) Ar, q = 0.85, vm = 1.0 m/s(c) CO2, q = 0.69, vm = 0.6 m/s(d) He, q = 0.83, vm = 1.0 m/s

Fig.5 Averaged OH fluorescence distributions 103

教職員学生と歩んだ燃焼研究(3)・両側の空気流速が異なるスリット火炎の浮き上がり特性に関する数値解析[第５０回燃焼シンポジウム 26-27 2012]104教職員学生と歩んだ燃焼研究(3)・燃料希薄メタン-空気予混合火炎の火炎構造と支配パラメータの解明
林直樹(2015 博士論文 受理)

実際の乱流火炎の局所・瞬間の火炎構造は層流火炎で模擬できる．

簡単な総括反応式

予混合気
多次元の複雑な流れ場

詳細な素反応機構

再
結
合
パ
ラ
メ
ー
タ
単純な流れ場：対向流105

教職員学生と歩んだ燃焼研究(3)熱発生速度最大値QmaxのCO,qおよび∇CO,qによる整理[日本機械学会論文集(B 編）, 72 (713), (2006), 208-215]

106

謝 辞
この38年間の永きにわたり，名古屋大学で教育・研究に携わることができたのは，ひとえに一緒に学

び，考え，躓きながらも歩み続けていただいた研究室の教職員，並びに大学院生208名，学部卒研生308名の献身的な協力のおかげです．ここに心からの謝意を表します．
最後に私ごとながら，わがままな研究生活を支えてくれた家族に感謝の意を表します．
107

永い間，
大変お世話になりました．
2015年3月6日
山下博史
108