特許 6358618 エンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法及びエンジンの燃焼シミュレーション方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6358618

(24)【登録日】2018年6月29日

(45)【発行日】2018年7月18日

(54)【発明の名称】エンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法及びエンジンの燃焼シミュレーション方法

(51)【国際特許分類】

   F02D 45/00 20060101AFI20180709BHJP

【ＦＩ】

   F02D45/00 372Z

   F02D45/00 370B

【請求項の数】8

【全頁数】19

(21)【出願番号】特願2014-214766(P2014-214766)

(22)【出願日】2014年10月21日

(65)【公開番号】特開2016-79947(P2016-79947A)

(43)【公開日】2016年5月16日

【審査請求日】2017年8月17日

(73)【特許権者】

【識別番号】000003137

【氏名又は名称】マツダ株式会社

(73)【特許権者】

【識別番号】899000068

【氏名又は名称】学校法人早稲田大学

(74)【代理人】

【識別番号】100059959

【弁理士】

【氏名又は名称】中村 稔

(74)【代理人】

【識別番号】100067013

【弁理士】

【氏名又は名称】大塚 文昭

(74)【代理人】

【識別番号】100082005

【弁理士】

【氏名又は名称】熊倉 禎男

(74)【代理人】

【識別番号】100086771

【弁理士】

【氏名又は名称】西島 孝喜

(74)【代理人】

【識別番号】100088694

【弁理士】

【氏名又は名称】弟子丸 健

(74)【代理人】

【識別番号】100168871

【弁理士】

【氏名又は名称】岩上 健

(72)【発明者】

【氏名】養祖 隆

(72)【発明者】

【氏名】山川 正尚

(72)【発明者】

【氏名】草鹿 仁

【審査官】 田村 佳孝

(56)【参考文献】

【文献】 国際公開第２０１２／１３１９４９（ＷＯ，Ａ１）

【文献】 特表２０１２−５３２２０５（ＪＰ，Ａ）

【文献】 特開２０１１−１４９７３９（ＪＰ，Ａ）

【文献】 特開２０１４−１７８７１７（ＪＰ，Ａ）

【文献】 特開２０１２−９２７６１（ＪＰ，Ａ）

【文献】 特開２００４−２３９１２９（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｆ０２Ｄ１３／００ − ２８／００

Ｆ０２Ｄ４１／００ − ４５／００

Ｇ０１Ｍ１５／００ − １５／１４

(57)【特許請求の範囲】

【請求項1】

エンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法であって、
上記燃焼シミュレーションを行う対象となる対象燃料に含まれる燃料成分種を設定する工程と、
上記燃料成分種の着火現象を記述する化学反応式の内、上記燃焼シミュレーションに使用する化学反応式を設定する工程と、
上記燃焼シミュレーションに使用する化学反応式の中から、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式を少なくとも１つ選択すると共に、熱着火時期の反応温度に対して相対的に感度の高い化学反応式を少なくとも１つ選択する工程と、
上記低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式の化学反応速度定数と、上記熱着火時期の反応温度に対して相対的に感度の高い化学反応式の化学反応速度定数とを、上記燃焼シミュレーションにより算出される低温酸化反応時期及び熱着火時期が所定の目標低温酸化反応時期及び目標熱着火時期にそれぞれ一致するように最適化する工程とを有することを特徴とする、エンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法。

【請求項2】

更に、上記対象燃料に添加される添加成分種を設定する工程と、
上記燃料成分種と上記添加成分種との混合割合を設定する工程と、
上記低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式を選択する工程において選択した化学反応式の化学反応速度定数を、上記混合割合で上記対象燃料に添加された上記添加成分種を上記燃料成分種の１つに置換した場合に上記燃焼シミュレーションにより算出される低温酸化反応時期及び熱着火時期が上記添加成分種を上記混合割合で添加した上記対象燃料における目標低温酸化反応時期及び目標熱着火時期にそれぞれに一致するように再度最適化する工程とを有する請求項１に記載のエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法。

【請求項3】

上記対象燃料は、上記燃料成分種としてイソオクタンとノルマルヘプタンとを含む混合燃料である請求項１又は２に記載のエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法。

【請求項4】

上記燃焼シミュレーションに使用する化学反応式を設定する工程において、下記のＲ１乃至Ｒ３８の化学反応式を設定する請求項３に記載のエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法。
Ｃ₇Ｈ₁₆＋Ｏ₂＝Ｃ₇Ｈ₁₅＋ＨＯ₂ （Ｒ１）
Ｃ₇Ｈ₁₅＋Ｏ₂＝Ｃ₇Ｈ₁₅Ｏ₂ （Ｒ２）
Ｃ₇Ｈ₁₅Ｏ₂＝Ｃ₇Ｈ₁₄ＯＯＨ （Ｒ３）
Ｃ₇Ｈ₁₄ＯＯＨ＋Ｏ₂＝Ｏ₂Ｃ₇Ｈ₁₄ＯＯＨ （Ｒ４）
Ｏ₂Ｃ₇Ｈ₁₄ＯＯＨ＝＞Ｃ₇ＫＥＴ＋ＯＨ （Ｒ５）
Ｃ₇ＫＥＴ＝＞Ｃ₅Ｈ₁₁ＣＯ＋ＣＨ₂Ｏ＋ＯＨ （Ｒ６）
Ｃ₅Ｈ₁₁ＣＯ＋Ｏ₂＝＞Ｃ₃Ｈ₆＋Ｃ₂Ｈ₄＋ＣＯ＋ＨＯ₂ （Ｒ７）
Ｃ₇Ｈ₁₆＋ＯＨ＝＞Ｃ₇Ｈ₁₅＋Ｈ₂Ｏ （Ｒ８）
Ｃ₇Ｈ₁₅＋Ｏ₂＝Ｃ₇Ｈ₁₄＋ＨＯ₂ （Ｒ９）
Ｃ₇Ｈ₁₄＋Ｏ₂＝＞Ｃ₅Ｈ₁₁＋ＣＨ₂Ｏ＋ＨＣＯ （Ｒ１０）
Ｃ₇Ｈ₁₆＋ＨＯ₂＝Ｃ₇Ｈ₁₅＋Ｈ₂Ｏ₂ （Ｒ１１）
Ｃ₇Ｈ₁₅＝＞Ｃ₅Ｈ₁₁＋Ｃ₂Ｈ₄ （Ｒ１２）
Ｃ₅Ｈ₁₁＝Ｃ₃Ｈ₇＋Ｃ₂Ｈ₄ （Ｒ１３）
Ｃ₈Ｈ₁₈＋Ｏ₂＝Ｃ₈Ｈ₁₇＋ＨＯ₂ （Ｒ１４）
Ｃ₈Ｈ₁₇＋Ｏ₂＝Ｃ₈Ｈ₁₇Ｏ₂ （Ｒ１５）
Ｃ₈Ｈ₁₇Ｏ₂＝Ｃ₈Ｈ₁₆ＯＯＨ （Ｒ１６）
Ｃ₈Ｈ₁₆ＯＯＨ＋Ｏ₂＝Ｏ₂Ｃ₈Ｈ₁₆ＯＯＨ （Ｒ１７）
Ｏ₂Ｃ₈Ｈ₁₆ＯＯＨ＝＞Ｃ₈ＫＥＴ＋ＯＨ （Ｒ１８）
Ｃ₈ＫＥＴ＝＞Ｃ₆Ｈ₁₃ＣＯ＋ＣＨ₂Ｏ＋ＯＨ （Ｒ１９）
Ｃ₆Ｈ₁₃ＣＯ＋Ｏ₂＝＞Ｃ₃Ｈ₆＋Ｃ₃Ｈ₆＋ＣＯ＋ＨＯ₂ （Ｒ２０）
Ｃ₈Ｈ₁₈＋ＯＨ＝＞Ｃ₈Ｈ₁₇＋Ｈ₂Ｏ （Ｒ２１）
Ｃ₈Ｈ₁₇＋Ｏ₂＝Ｃ₈Ｈ₁₆＋ＨＯ₂ （Ｒ２２）
Ｃ₈Ｈ₁₆＋Ｏ₂＝＞Ｃ₆Ｈ₁₃＋ＣＨ₂Ｏ＋ＨＣＯ （Ｒ２３）
Ｃ₈Ｈ₁₈＋ＨＯ₂＝Ｃ₈Ｈ₁₇＋Ｈ₂Ｏ₂ （Ｒ２４）
Ｃ₈Ｈ₁₇＝＞Ｃ₆Ｈ₁₃＋Ｃ₂Ｈ₄ （Ｒ２５）
Ｃ₆Ｈ₁₃＝Ｃ₃Ｈ₇＋Ｃ₃Ｈ₆ （Ｒ２６）
Ｃ₃Ｈ₇＋Ｏ₂＝Ｃ₃Ｈ₆＋ＨＯ₂ （Ｒ２７）
Ｃ₃Ｈ₆＋Ｃ₃Ｈ₆＝＞Ｃ₂Ｈ₄＋Ｃ₂Ｈ₄＋Ｃ₂Ｈ₄ （Ｒ２８）
Ｃ₂Ｈ₄＋ＯＨ＝Ｃ₂Ｈ₃＋Ｈ₂Ｏ （Ｒ２９）
Ｃ₂Ｈ₃＋Ｏ₂＝ＣＨ₂Ｏ＋ＨＣＯ （Ｒ３０）
ＣＨ₂Ｏ＋ＯＨ＝ＨＣＯ＋Ｈ₂Ｏ （Ｒ３１）
ＨＣＯ＋Ｏ₂＝ＣＯ＋ＨＯ₂ （Ｒ３２）
ＣＯ＋ＯＨ＝ＣＯ₂＋Ｈ （Ｒ３３）
Ｈ₂Ｏ₂＋ＯＨ＝ＨＯ₂＋Ｈ₂Ｏ （Ｒ３４）
Ｈ＋Ｏ₂＋Ｍ＝ＨＯ₂＋Ｍ （Ｒ３５）
ＨＯ₂＋ＨＯ₂＝Ｈ₂Ｏ₂＋Ｏ₂ （Ｒ３６）
Ｈ₂Ｏ₂＋Ｍ＝ＯＨ＋ＯＨ＋Ｍ （Ｒ３７）
Ｃ₂Ｈ₄＋Ｏ₂＝＞ＣＨ₂Ｏ＋ＣＨ₂Ｏ （Ｒ３８）

【請求項5】

上記低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式は、Ｃ₇ＫＥＴ＝＞Ｃ₅Ｈ₁₁ＣＯ＋ＣＨ₂Ｏ＋ＯＨ （Ｒ６）及びＣ₈Ｈ₁₈＋ＯＨ＝＞Ｃ₈Ｈ₁₇＋Ｈ₂Ｏ （Ｒ２１）であり、
上記熱着火時期の反応温度に対して相対的に感度の高い化学反応式は、Ｈ₂Ｏ₂＋Ｍ＝ＯＨ＋ＯＨ＋Ｍ （Ｒ３７）である請求項４に記載のエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法。

【請求項6】

上記低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式の化学反応速度定数と、上記熱着火時期の反応温度に対して相対的に感度の高い化学反応式の化学反応速度定数とを最適化する工程において、上記燃焼シミュレーションにより算出される低温酸化反応時期をτ₁、上記燃焼シミュレーションにより算出される熱着火時期をτ₂、上記目標低温酸化反応時期をτ₁^tgt、上記目標熱着火時期をτ₂^tgt、燃焼の初期温度をＴとした場合、以下の式により表される評価関数Ｅ_fの最小値を得るように、粒子群最適化法により最適化を行う請求項１乃至５の何れか１項に記載のエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法。

【数1】

【請求項7】

上記低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式を選択する工程において選択した化学反応式の化学反応速度定数を再度最適化する工程において、上記添加成分種を上記燃料成分種の１つに置換した場合に上記燃焼シミュレーションにより算出される低温酸化反応時期及び熱着火時期をそれぞれτ₁^mix及びτ₂^mixとし、上記添加成分種を上記混合割合で添加した上記対象燃料における目標低温酸化反応時期及び目標熱着火時期をそれぞれτ₁^mixtgt及びτ₂^mixtgtとし、燃焼の初期温度をＴとした場合、以下の式により表される評価関数Ｅ_fの最小値を得るように、粒子群最適化法により最適化を行う請求項２に記載のエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法。

【数2】

【請求項8】

請求項１乃至７の何れか１項に記載のエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法により決定した化学反応速度定数を使用して、上記対象燃料による燃焼シミュレーションを行うエンジンの燃焼シミュレーション方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、化学反応シミュレーションに使用する化学反応速度定数の決定方法及び燃焼シミュレーション方法に係わり、特に、エンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法及びエンジンの燃焼シミュレーション方法に関する。

【背景技術】

【0002】

従来、ガソリンエンジン等の内燃機関における燃焼反応を分析するために、エンジンの燃焼シミュレーションが利用されている。この燃焼シミュレーションでは、燃料に含まれる各種成分や中間生成物等、燃焼に関連する化学種の反応を多数の化学反応式によって記述すると共に、圧縮工程及び膨張行程における燃焼室の容積変化を容器モデルにより模擬することによって、その容器中での燃焼反応の進行に伴う温度変化や組成変化を再現している（例えば特許文献１参照）。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開２００４−２３９１２９号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

ところで、近年では、ガソリンエンジンの高効率化及び低公害化のために高圧縮比の予混合圧縮自己着火（ＨＣＣＩ）燃焼の実用化が期待されており、そのためにＨＣＣＩ燃焼の着火性や火花点火（ＳＩ）燃焼におけるプリイグニッション、ノッキングといった自己着火に対して燃料に含まれる燃料成分種が及ぼす影響を正確に把握することが望まれている。しかしながら、これらの着火現象に対する燃料成分種の影響を燃焼シミュレーションによって精度よく再現するためには、多数の化学種及び化学反応式を使用した計算を行わなければならず、膨大な演算量が要求される。
また、ＳＩ燃焼の分析においては、プリイグニッションやノッキングの発生位置の把握などのために３次元流体計算との連成計算が望まれるが、上述のように多数の化学種及び化学反応式を使用した場合には計算負荷の増大が顕著である。従って、従来の燃焼シミュレーションでは、使用する化学種及び化学反応式の数を絞り込み、エンジンの燃焼における代表的な化学反応のみを再現することによって、シミュレーションの精度をある程度犠牲にしつつ、計算負荷の増大を抑制している。

【0005】

本発明は、上述した従来技術の問題点を解決するためになされたものであり、エンジンの燃焼反応を少ない計算負荷で精度よく再現することができる、エンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法及びエンジンの燃焼シミュレーション方法を提供することを目的とする。

【課題を解決するための手段】

【0006】

上記の目的を達成するために、本発明のエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法は、燃焼シミュレーションを行う対象となる対象燃料に含まれる燃料成分種を設定する工程と、燃料成分種の着火現象を記述する化学反応式の内、燃焼シミュレーションに使用する化学反応式を設定する工程と、燃焼シミュレーションに使用する化学反応式の中から、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式を少なくとも１つ選択すると共に、熱着火時期の反応温度に対して相対的に感度の高い化学反応式を少なくとも１つ選択する工程と、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式の化学反応速度定数と、熱着火時期の反応温度に対して相対的に感度の高い化学反応式の化学反応速度定数とを、燃焼シミュレーションにより算出される低温酸化反応時期及び熱着火時期が所定の目標低温酸化反応時期及び目標熱着火時期にそれぞれ一致するように最適化する工程とを有することを特徴とする。
このように構成された本発明においては、燃焼シミュレーションに使用する化学反応式の中から、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式を少なくとも１つ選択すると共に、熱着火時期の反応温度に対して相対的に感度の高い化学反応式を少なくとも１つ選択し、それらの選択した化学反応式の化学反応速度定数を、燃焼シミュレーションにより算出される低温酸化反応時期及び熱着火時期が所定の目標低温酸化反応時期及び目標熱着火時期にそれぞれ一致するように最適化するので、高い精度で低温酸化反応時期及び熱着火時期を算出することが可能な簡略化化学反応スキームを作成することができ、これにより、エンジンの燃焼反応を少ない計算負荷で精度よく再現することができる。

【0007】

また、本発明において、好ましくは、更に、対象燃料に添加される添加成分種を設定する工程と、燃料成分種と添加成分種との混合割合を設定する工程と、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式を選択する工程において選択した化学反応式の化学反応速度定数を、混合割合で対象燃料に添加された添加成分種を燃料成分種の１つに置換した場合に燃焼シミュレーションにより算出される低温酸化反応時期及び熱着火時期が添加成分種を混合割合で添加した対象燃料における目標低温酸化反応時期及び目標熱着火時期にそれぞれに一致するように再度最適化する工程とを有する。
このように構成された本発明においては、混合燃料についても高い精度で低温酸化反応時期及び熱着火時期を算出することが可能な簡略化化学反応スキームを作成することができ、これにより、様々な混合燃料についてエンジンの燃焼反応を少ない計算負荷で精度よく再現することができる。

【0008】

また、本発明において、好ましくは、対象燃料は、燃料成分種としてイソオクタンとノルマルヘプタンとを含む混合燃料である。
このように構成された本発明においては、イソオクタンとノルマルヘプタンとを含む混合燃料についてエンジンの燃焼反応を少ない計算負荷で精度よく再現することができる。

【0009】

また、本発明において、好ましくは、燃焼シミュレーションに使用する化学反応式を設定する工程において、下記のＲ１乃至Ｒ３８の化学反応式を設定する。
Ｃ₇Ｈ₁₆＋Ｏ₂＝Ｃ₇Ｈ₁₅＋ＨＯ₂ （Ｒ１）
Ｃ₇Ｈ₁₅＋Ｏ₂＝Ｃ₇Ｈ₁₅Ｏ₂ （Ｒ２）
Ｃ₇Ｈ₁₅Ｏ₂＝Ｃ₇Ｈ₁₄ＯＯＨ （Ｒ３）
Ｃ₇Ｈ₁₄ＯＯＨ＋Ｏ₂＝Ｏ₂Ｃ₇Ｈ₁₄ＯＯＨ （Ｒ４）
Ｏ₂Ｃ₇Ｈ₁₄ＯＯＨ＝＞Ｃ₇ＫＥＴ＋ＯＨ （Ｒ５）
Ｃ₇ＫＥＴ＝＞Ｃ₅Ｈ₁₁ＣＯ＋ＣＨ₂Ｏ＋ＯＨ （Ｒ６）
Ｃ₅Ｈ₁₁ＣＯ＋Ｏ₂＝＞Ｃ₃Ｈ₆＋Ｃ₂Ｈ₄＋ＣＯ＋ＨＯ₂ （Ｒ７）
Ｃ₇Ｈ₁₆＋ＯＨ＝＞Ｃ₇Ｈ₁₅＋Ｈ₂Ｏ （Ｒ８）
Ｃ₇Ｈ₁₅＋Ｏ₂＝Ｃ₇Ｈ₁₄＋ＨＯ₂ （Ｒ９）
Ｃ₇Ｈ₁₄＋Ｏ₂＝＞Ｃ₅Ｈ₁₁＋ＣＨ₂Ｏ＋ＨＣＯ （Ｒ１０）
Ｃ₇Ｈ₁₆＋ＨＯ₂＝Ｃ₇Ｈ₁₅＋Ｈ₂Ｏ₂ （Ｒ１１）
Ｃ₇Ｈ₁₅＝＞Ｃ₅Ｈ₁₁＋Ｃ₂Ｈ₄ （Ｒ１２）
Ｃ₅Ｈ₁₁＝Ｃ₃Ｈ₇＋Ｃ₂Ｈ₄ （Ｒ１３）
Ｃ₈Ｈ₁₈＋Ｏ₂＝Ｃ₈Ｈ₁₇＋ＨＯ₂ （Ｒ１４）
Ｃ₈Ｈ₁₇＋Ｏ₂＝Ｃ₈Ｈ₁₇Ｏ₂ （Ｒ１５）
Ｃ₈Ｈ₁₇Ｏ₂＝Ｃ₈Ｈ₁₆ＯＯＨ （Ｒ１６）
Ｃ₈Ｈ₁₆ＯＯＨ＋Ｏ₂＝Ｏ₂Ｃ₈Ｈ₁₆ＯＯＨ （Ｒ１７）
Ｏ₂Ｃ₈Ｈ₁₆ＯＯＨ＝＞Ｃ₈ＫＥＴ＋ＯＨ （Ｒ１８）
Ｃ₈ＫＥＴ＝＞Ｃ₆Ｈ₁₃ＣＯ＋ＣＨ₂Ｏ＋ＯＨ （Ｒ１９）
Ｃ₆Ｈ₁₃ＣＯ＋Ｏ₂＝＞Ｃ₃Ｈ₆＋Ｃ₃Ｈ₆＋ＣＯ＋ＨＯ₂ （Ｒ２０）
Ｃ₈Ｈ₁₈＋ＯＨ＝＞Ｃ₈Ｈ₁₇＋Ｈ₂Ｏ （Ｒ２１）
Ｃ₈Ｈ₁₇＋Ｏ₂＝Ｃ₈Ｈ₁₆＋ＨＯ₂ （Ｒ２２）
Ｃ₈Ｈ₁₆＋Ｏ₂＝＞Ｃ₆Ｈ₁₃＋ＣＨ₂Ｏ＋ＨＣＯ （Ｒ２３）
Ｃ₈Ｈ₁₈＋ＨＯ₂＝Ｃ₈Ｈ₁₇＋Ｈ₂Ｏ₂ （Ｒ２４）
Ｃ₈Ｈ₁₇＝＞Ｃ₆Ｈ₁₃＋Ｃ₂Ｈ₄ （Ｒ２５）
Ｃ₆Ｈ₁₃＝Ｃ₃Ｈ₇＋Ｃ₃Ｈ₆ （Ｒ２６）
Ｃ₃Ｈ₇＋Ｏ₂＝Ｃ₃Ｈ₆＋ＨＯ₂ （Ｒ２７）
Ｃ₃Ｈ₆＋Ｃ₃Ｈ₆＝＞Ｃ₂Ｈ₄＋Ｃ₂Ｈ₄＋Ｃ₂Ｈ₄ （Ｒ２８）
Ｃ₂Ｈ₄＋ＯＨ＝Ｃ₂Ｈ₃＋Ｈ₂Ｏ （Ｒ２９）
Ｃ₂Ｈ₃＋Ｏ₂＝ＣＨ₂Ｏ＋ＨＣＯ （Ｒ３０）
ＣＨ₂Ｏ＋ＯＨ＝ＨＣＯ＋Ｈ₂Ｏ （Ｒ３１）
ＨＣＯ＋Ｏ₂＝ＣＯ＋ＨＯ₂ （Ｒ３２）
ＣＯ＋ＯＨ＝ＣＯ₂＋Ｈ （Ｒ３３）
Ｈ₂Ｏ₂＋ＯＨ＝ＨＯ₂＋Ｈ₂Ｏ （Ｒ３４）
Ｈ＋Ｏ₂＋Ｍ＝ＨＯ₂＋Ｍ （Ｒ３５）
ＨＯ₂＋ＨＯ₂＝Ｈ₂Ｏ₂＋Ｏ₂ （Ｒ３６）
Ｈ₂Ｏ₂＋Ｍ＝ＯＨ＋ＯＨ＋Ｍ （Ｒ３７）
Ｃ₂Ｈ₄＋Ｏ₂＝＞ＣＨ₂Ｏ＋ＣＨ₂Ｏ （Ｒ３８）
このように構成された本発明においては、イソオクタンとノルマルヘプタンとを含む混合燃料についてエンジンの燃焼反応を少ない計算負荷で一層精度よく再現することができる。

【0010】

また、本発明において、好ましくは、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式は、Ｃ₇ＫＥＴ＝＞Ｃ₅Ｈ₁₁ＣＯ＋ＣＨ₂Ｏ＋ＯＨ （Ｒ６）及びＣ₈Ｈ₁₈＋ＯＨ＝＞Ｃ₈Ｈ₁₇＋Ｈ₂Ｏ （Ｒ２１）であり、熱着火時期の反応温度に対して相対的に感度の高い化学反応式は、Ｈ₂Ｏ₂＋Ｍ＝ＯＨ＋ＯＨ＋Ｍ （Ｒ３７）である。
このように構成された本発明においては、低温酸化反応時期及び熱着火時期を精度よく計算することが可能な簡略化化学反応スキームを作成できる。

【0011】

また、本発明において、好ましくは、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式の化学反応速度定数と、熱着火時期の反応温度に対して相対的に感度の高い化学反応式の化学反応速度定数とを最適化する工程において、燃焼シミュレーションにより算出される低温酸化反応時期をτ₁、燃焼シミュレーションにより算出される熱着火時期をτ₂、目標低温酸化反応時期をτ₁^tgt、目標熱着火時期をτ₂^tgt、燃焼の初期温度をＴとした場合、以下の式により表される評価関数Ｅ_fの最小値を得るように、粒子群最適化法により最適化を行う。

【数1】

このように構成された本発明においては、高い精度で燃焼シミュレーションを実行することが可能な簡略化化学反応スキームを効率よく作成することができる。

【0012】

また、本発明において、好ましくは、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式を選択する工程において選択した化学反応式の化学反応速度定数を再度最適化する工程において、添加成分種を燃料成分種の１つに置換した場合に燃焼シミュレーションにより算出される低温酸化反応時期及び熱着火時期をそれぞれτ₁^mix及びτ₂^mixとし、添加成分種を混合割合で添加した対象燃料における目標低温酸化反応時期及び目標熱着火時期をそれぞれτ₁^mixtgt及びτ₂^mixtgtとし、燃焼の初期温度をＴとした場合、以下の式により表される評価関数Ｅ_fの最小値を得るように、粒子群最適化法により最適化を行う。

【数2】

このように構成された本発明においては、様々な混合燃料について高い精度で燃焼シミュレーションを実行することが可能な簡略化化学反応スキームを効率よく作成することができる。

【0013】

また、本発明のエンジンの燃焼シミュレーション方法によれば、上述のエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法により決定した化学反応速度定数を使用して、対象燃料による燃焼シミュレーションを行う。
このように構成された本発明においては、高い精度で低温酸化反応時期及び熱着火時期を算出することが可能な簡略化化学反応スキームを使用して、エンジンの燃焼反応を少ない計算負荷で精度よく再現することができる。

【発明の効果】

【0014】

本発明によるエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法及びエンジンの燃焼シミュレーション方法によれば、エンジンの燃焼反応を少ない計算負荷で精度よく再現することができる。

【図面の簡単な説明】

【0015】

【図1】本発明の実施形態によるエンジンの燃焼シミュレーションシステムのブロック図である。

【図2】本発明の実施形態によるエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法のフローチャートである。

【図3】ＰＲＦ７０について本発明の実施形態による燃焼シミュレーションにより算出した着火遅れ時間とＭＦＣスキームを使用した燃焼シミュレーションにより算出した着火遅れ時間とを比較した線図である。

【図4】ＰＲＦ７０について本発明の実施形態による燃焼シミュレーションにより算出した着火遅れ時間とＭＦＣスキームを使用した燃焼シミュレーションにより算出した着火遅れ時間とを比較した線図である。

【図5】本発明の実施形態による混合燃料の燃焼シミュレーションに使用する化学反応速度定数の決定方法のフローチャートである。

【図6】ＰＲＦ７０にトルエンを添加した混合燃料について本発明の実施形態による燃焼シミュレーションにより算出した着火遅れ時間とＭＦＣスキームを使用した燃焼シミュレーションにより算出した着火遅れ時間とを比較した線図である。

【図7】ＰＲＦ７０に対するトルエンの混合割合を変更した場合において本発明の実施形態によるエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法により決定した化学反応速度定数の変化率を示す線図である。

【図8】本発明の実施形態による燃焼シミュレーション方法を示すフローチャートである。

【図9】ＰＲＦ７０にトルエン及び１−ペンテンを添加した混合燃料について本発明の実施形態による燃焼シミュレーションにより算出した着火遅れ時間とＭＦＣスキームを使用した燃焼シミュレーションにより算出した着火遅れ時間とを比較した線図である。

【発明を実施するための形態】

【0016】

以下、添付図面を参照して、本発明の実施形態によるエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法及びエンジンの燃焼シミュレーション方法を説明する。
まず、図１により、本発明の実施形態によるエンジンの燃焼シミュレーション方法を実行する燃焼シミュレーションシステムの全体構成を説明する。図１は、本発明の実施形態によるエンジンの燃焼シミュレーションシステムのブロック図である。

【0017】

図１において符号１は、本発明の実施形態によるエンジンの燃焼シミュレーションシステムを示す。この図１に示すように、エンジンの燃焼シミュレーションシステム１は、燃焼シミュレーションにおける各種の設定データを入力する入力装置２、燃焼シミュレーションを実行する燃焼シミュレーション装置４、燃焼シミュレーションの結果を出力する出力装置６、及び、燃焼シミュレーションにおいて化学反応の演算に使用されるデータを格納した化学反応データベース８を有している。

【0018】

燃焼シミュレーション装置４は、燃焼シミュレーションの演算に使用されるエンジンの物理モデルを構築するプリプロセッサ１０、エンジンの燃焼反応の解析を行う化学反応解析部１２、エンジンの燃焼室における燃焼ガスの流れを解析する３次元流体解析部１４、及び、化学反応解析部１２及び３次元流体解析部１４による解析結果を出力装置６により出力するためのポスト処理を行うポストプロセッサ１６を備えている。

【0019】

化学反応データベース８は、エンジンの燃焼における代表的な化学反応式を抽出し、各化学反応式によって記述される素反応における化学反応速度定数及び平衡定数を対応付けた簡略化化学反応スキーム１８と、簡略化化学反応スキーム１８に含まれている特定の化学反応式について、燃料に添加された添加成分種の混合割合と化学反応速度との関係を示した化学反応速度定数マップ２０とを格納している。

【0020】

次に、図２乃至図４により、本発明の実施形態における簡略化化学反応スキーム１８の作成方法について説明する。図２は本発明の実施形態によるエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法のフローチャートであり、図３はＰＲＦ７０について本発明の実施形態による燃焼シミュレーションにより算出した着火遅れ時間とＭＦＣスキームを使用した燃焼シミュレーションにより算出した着火遅れ時間とを比較した線図であり、図４はＰＲＦ７０について本発明の実施形態による燃焼シミュレーションにより算出した着火遅れ時間とＭＦＣスキームを使用した燃焼シミュレーションにより算出した着火遅れ時間とを比較した線図である。

【0021】

まず、図２に示すフローは、エンジンの燃焼シミュレーションに先立って予め実行される。図２に示すように、フローが開始されると、ステップＳ１において、エンジンの燃焼シミュレーションを行う対象となる対象燃料に含まれる燃料成分種を設定する。本実施形態では、対象燃料は、燃料成分種としてイソオクタンとノルマルヘプタンとをイソオクタン：ノルマルヘプタン＝７０：３０（ｖｏｌ．％）の混合割合で含む混合燃料（ＰＲＦ７０）である。

【0022】

次に、ステップＳ２において、本発明の実施形態における簡略化化学反応スキーム１８のベースとなる初期簡略化化学反応スキームを設定する。本実施形態では、対象燃料のＰＲＦの着火現象に関連する多数の化学反応式の中から、その着火現象の再現に適した化学反応式として知られている下記のＲ１乃至Ｒ３８の化学反応式を抽出し、各化学反応式によって記述される素反応における化学反応速度定数及び平衡定数を対応付けたものを初期簡略化化学反応スキームとして設定する。
Ｃ₇Ｈ₁₆＋Ｏ₂＝Ｃ₇Ｈ₁₅＋ＨＯ₂ （Ｒ１）
Ｃ₇Ｈ₁₅＋Ｏ₂＝Ｃ₇Ｈ₁₅Ｏ₂ （Ｒ２）
Ｃ₇Ｈ₁₅Ｏ₂＝Ｃ₇Ｈ₁₄ＯＯＨ （Ｒ３）
Ｃ₇Ｈ₁₄ＯＯＨ＋Ｏ₂＝Ｏ₂Ｃ₇Ｈ₁₄ＯＯＨ （Ｒ４）
Ｏ₂Ｃ₇Ｈ₁₄ＯＯＨ＝＞Ｃ₇ＫＥＴ＋ＯＨ （Ｒ５）
Ｃ₇ＫＥＴ＝＞Ｃ₅Ｈ₁₁ＣＯ＋ＣＨ₂Ｏ＋ＯＨ （Ｒ６）
Ｃ₅Ｈ₁₁ＣＯ＋Ｏ₂＝＞Ｃ₃Ｈ₆＋Ｃ₂Ｈ₄＋ＣＯ＋ＨＯ₂ （Ｒ７）
Ｃ₇Ｈ₁₆＋ＯＨ＝＞Ｃ₇Ｈ₁₅＋Ｈ₂Ｏ （Ｒ８）
Ｃ₇Ｈ₁₅＋Ｏ₂＝Ｃ₇Ｈ₁₄＋ＨＯ₂ （Ｒ９）
Ｃ₇Ｈ₁₄＋Ｏ₂＝＞Ｃ₅Ｈ₁₁＋ＣＨ₂Ｏ＋ＨＣＯ （Ｒ１０）
Ｃ₇Ｈ₁₆＋ＨＯ₂＝Ｃ₇Ｈ₁₅＋Ｈ₂Ｏ₂ （Ｒ１１）
Ｃ₇Ｈ₁₅＝＞Ｃ₅Ｈ₁₁＋Ｃ₂Ｈ₄ （Ｒ１２）
Ｃ₅Ｈ₁₁＝Ｃ₃Ｈ₇＋Ｃ₂Ｈ₄ （Ｒ１３）
Ｃ₈Ｈ₁₈＋Ｏ₂＝Ｃ₈Ｈ₁₇＋ＨＯ₂ （Ｒ１４）
Ｃ₈Ｈ₁₇＋Ｏ₂＝Ｃ₈Ｈ₁₇Ｏ₂ （Ｒ１５）
Ｃ₈Ｈ₁₇Ｏ₂＝Ｃ₈Ｈ₁₆ＯＯＨ （Ｒ１６）
Ｃ₈Ｈ₁₆ＯＯＨ＋Ｏ₂＝Ｏ₂Ｃ₈Ｈ₁₆ＯＯＨ （Ｒ１７）
Ｏ₂Ｃ₈Ｈ₁₆ＯＯＨ＝＞Ｃ₈ＫＥＴ＋ＯＨ （Ｒ１８）
Ｃ₈ＫＥＴ＝＞Ｃ₆Ｈ₁₃ＣＯ＋ＣＨ₂Ｏ＋ＯＨ （Ｒ１９）
Ｃ₆Ｈ₁₃ＣＯ＋Ｏ₂＝＞Ｃ₃Ｈ₆＋Ｃ₃Ｈ₆＋ＣＯ＋ＨＯ₂ （Ｒ２０）
Ｃ₈Ｈ₁₈＋ＯＨ＝＞Ｃ₈Ｈ₁₇＋Ｈ₂Ｏ （Ｒ２１）
Ｃ₈Ｈ₁₇＋Ｏ₂＝Ｃ₈Ｈ₁₆＋ＨＯ₂ （Ｒ２２）
Ｃ₈Ｈ₁₆＋Ｏ₂＝＞Ｃ₆Ｈ₁₃＋ＣＨ₂Ｏ＋ＨＣＯ （Ｒ２３）
Ｃ₈Ｈ₁₈＋ＨＯ₂＝Ｃ₈Ｈ₁₇＋Ｈ₂Ｏ₂ （Ｒ２４）
Ｃ₈Ｈ₁₇＝＞Ｃ₆Ｈ₁₃＋Ｃ₂Ｈ₄ （Ｒ２５）
Ｃ₆Ｈ₁₃＝Ｃ₃Ｈ₇＋Ｃ₃Ｈ₆ （Ｒ２６）
Ｃ₃Ｈ₇＋Ｏ₂＝Ｃ₃Ｈ₆＋ＨＯ₂ （Ｒ２７）
Ｃ₃Ｈ₆＋Ｃ₃Ｈ₆＝＞Ｃ₂Ｈ₄＋Ｃ₂Ｈ₄＋Ｃ₂Ｈ₄ （Ｒ２８）
Ｃ₂Ｈ₄＋ＯＨ＝Ｃ₂Ｈ₃＋Ｈ₂Ｏ （Ｒ２９）
Ｃ₂Ｈ₃＋Ｏ₂＝ＣＨ₂Ｏ＋ＨＣＯ （Ｒ３０）
ＣＨ₂Ｏ＋ＯＨ＝ＨＣＯ＋Ｈ₂Ｏ （Ｒ３１）
ＨＣＯ＋Ｏ₂＝ＣＯ＋ＨＯ₂ （Ｒ３２）
ＣＯ＋ＯＨ＝ＣＯ₂＋Ｈ （Ｒ３３）
Ｈ₂Ｏ₂＋ＯＨ＝ＨＯ₂＋Ｈ₂Ｏ （Ｒ３４）
Ｈ＋Ｏ₂＋Ｍ＝ＨＯ₂＋Ｍ （Ｒ３５）
ＨＯ₂＋ＨＯ₂＝Ｈ₂Ｏ₂＋Ｏ₂ （Ｒ３６）
Ｈ₂Ｏ₂＋Ｍ＝ＯＨ＋ＯＨ＋Ｍ （Ｒ３７）
Ｃ₂Ｈ₄＋Ｏ₂＝＞ＣＨ₂Ｏ＋ＣＨ₂Ｏ （Ｒ３８）

【0023】

次に、ステップＳ３において、ステップＳ２で設定した初期簡略化化学反応スキームに含まれる化学反応式により記述された各素反応について、対象燃料の低温酸化反応終了時点でのＨＣＨＯ及びＨ₂Ｏ₂の生成量に対する感度解析を行う。
本実施形態においては、ＰＲＦにおいて低温酸化反応が顕著にみられる反応温度７００Ｋ及び８００Ｋの条件下で感度解析を行う。この場合、上述の化学反応式Ｒ６（Ｃ₇ＫＥＴ＝＞Ｃ₅Ｈ₁₁ＣＯ＋ＣＨ₂Ｏ＋ＯＨ）及びＲ２１（Ｃ₈Ｈ₁₈＋ＯＨ＝＞Ｃ₈Ｈ₁₇＋Ｈ₂Ｏ）が特に高い感度を持つ化学反応式として抽出される。
ここで、化学反応式Ｒ６はノルマルヘプタンの反応機構の中でＨＣＨＯを生成する反応であり、この反応を促進させるとＨＣＨＯ、Ｈ₂Ｏ₂が共に増加する。一方、Ｒ２１はイソオクタンからの水素引き抜き反応であり、この反応を促進させるとＨＣＨＯ、Ｈ₂Ｏ₂が共に減少する。

【0024】

次に、ステップＳ４において、ステップＳ２で設定した初期簡略化化学反応スキームに含まれる化学反応式により記述された各素反応について、対象燃料の熱着火時期の反応温度に対する感度解析を行う。
本実施形態においては、ＰＲＦ７０の熱着火時期の反応温度に対する感度解析を行う。この場合、上述の化学反応式Ｒ３７（Ｈ₂Ｏ₂＋Ｍ＝ＯＨ＋ＯＨ＋Ｍ）が特に高い感度を持つ化学反応式として抽出される。

【0025】

次に、ステップＳ５において、ステップＳ３及びステップＳ４の感度解析結果に基づき、初期簡略化化学反応スキームに含まれる化学反応式の中から、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式を少なくとも１つ選択すると共に、熱着火時期の反応温度に対して相対的に感度の高い化学反応式を少なくとも１つ選択する。
本実施形態においては、ＰＲＦの低温酸化反応終了時点でのＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して特に高い感度を持つ化学反応式Ｒ６及びＲ２１と、ＰＲＦ７０の熱着火時期の反応温度に対して特に高い感度を持つ化学反応式Ｒ３７を選択する。

【0026】

次に、ステップＳ６において、対象燃料の燃焼反応における正確な低温酸化反応時期を目標低温酸化反応時期τ₁^tgtとして取得すると共に、対象燃料の燃焼反応における正確な熱着火時期を目標熱着火時期τ₂^tgtとして取得する。ここで、低温酸化反応時期とは、燃焼反応の初期温度から５０Ｋ上昇した時刻をいい、熱着火時期とは、燃焼反応の初期温度から５００Ｋ上昇した時刻をいうものとする。
本実施形態においては、多様な成分種を含む燃料の着火現象を詳細に表現した化学反応スキームとして、Ｍｏｄｅｌ Ｆｕｅｌｓ Ｃｏｎｓｏｒｔｉｕｍで開発されたＭＦＣ ２０１１ Ｇａｓｏｌｉｎｅ ＭｅｃｈａｎｉｓｍとＭＦＣ ２００８ ＥＴＢＥ Ｍｅｃｈａｎｉｓｍをマージした化学反応スキーム（以下、ＭＦＣスキームという）を使用して、ＰＲＦ７０について圧力６ＭＰａ、空気過剰率（λ）＝１の条件で初期温度７００Ｋから１４００Ｋまで５０Ｋ毎に化学反応計算を行うことにより、各温度における目標低温酸化反応時期τ₁^tgt及び目標熱着火時期τ₂^tgtを取得する。なお、ＭＦＣスキームに含まれる化学種数は１１１１６、化学反応式数は２３０１である。

【0027】

次に、ステップＳ７において、ステップＳ５で選択した各化学反応式の化学反応速度定数を、簡略化化学反応スキーム１８を使用した燃焼シミュレーションにより算出される低温酸化反応時期τ₁及び熱着火時期τ₂がステップＳ６で取得した目標低温酸化反応時期τ₁^tgt及び目標熱着火時期τ₂^tgtにそれぞれ一致するように最適化する。
本実施形態においては、化学反応式Ｒ６、Ｒ２１及びＲ３７の化学反応速度定数を、以下の式（１）により表される評価関数Ｅ_fの最小値を得るように、粒子群最適化法（ＰＳＯ法）により最適化を行う。

【数3】

【0028】

ＰＳＯ法では、最適化対象のパラメータについて探索粒子と呼ばれる複数の値の組み合わせを用意し、粒子群の情報を活用して評価関数の値が目標値に近づくよう反復計算によりパラメータを最適化する。それぞれの探索粒子は反復回ごとに値が決定されるが、その決定にはその反復回までに、（１）その粒子自身が最もよい評価関数値を得たときの値（パーソナルベスト）と、（２）全探索粒子の中で最もよい評価関数値を得たときの値（グループベスト）を用いる。これにより、反復回数が進むにしたがって目標に近い領域に探索粒子を集中的に配置することができるため、効率よく最適化できる。本実施形態では、探索粒子数１０、反復回数１００とする。各反復回においては、まず、圧力６ＭＰａ、空気過剰率（λ）１の条件で初期温度７００Ｋから１４００Ｋまで５０Ｋ毎に化学反応計算を行い、各温度での低温酸化反応時期τ₁及び熱着火時期τ₂を求める。次に、ＭＦＣスキームに対する簡略化化学反応スキーム１８の誤差率をτ₁とτ₂のそれぞれについて求め、最後に、各誤差率の絶対値を全温度で積算した値を評価関数として用いる。
ステップＳ７の後、エンジンの燃焼シミュレーションに使用する化学反応速度定数を決定するフローを終了する。

【0029】

簡略化化学反応スキーム１８に含まれる化学反応式Ｒ６、Ｒ２１、及びＲ３７の化学反応速度定数を図２のステップＳ７において最適化した後の簡略化化学反応スキーム１８を使用して燃焼シミュレーションを行うことにより、図３に示すように、圧力６ＭＰａ、空気過剰率（λ）１の条件においては、初期温度７００Ｋから１４００Ｋまでの温度範囲において、初期簡略化化学反応スキームを使用した場合と比較してＭＦＣスキームによるシミュレーション結果に良く一致することが分かる。
更に、この最適化した後の簡略化化学反応スキーム１８を使用して、空気過剰率λ＝２及びλ＝３の条件で燃焼シミュレーションを行った場合においても、図４に示すようにＭＦＣスキームによるシミュレーション結果に良く一致することが分かる。λ＝２あるいはλ＝３というリーン条件においては、低温酸化反応が相対的に小さくなり、低温度領域での着火遅れ時間の延長が特徴的だが、この特性についてもよい一致を維持できている。

【0030】

次に、図５乃至図７により、対象燃料に添加成分種を添加した混合燃料の燃焼シミュレーションに使用可能な簡略化化学反応スキーム１８の作成方法について説明する。図５は混合燃料の燃焼シミュレーションに使用する化学反応速度定数の決定方法のフローチャートであり、図６はＰＲＦ７０にトルエンを添加した混合燃料について本発明の実施形態による燃焼シミュレーションにより算出した着火遅れ時間とＭＦＣスキームを使用した燃焼シミュレーションにより算出した着火遅れ時間とを比較した線図であり、図７はＰＲＦ７０に対するトルエンの混合割合を変更した場合において本発明の実施形態によるエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法により決定した化学反応速度定数の変化率を示す線図である。

【0031】

まず、図５に示すフローは、図２のフローで簡略化化学反応スキーム１８を作成した後、エンジンの燃焼シミュレーションに先立って予め実行される。図５に示すように、フローが開始されると、ステップＳ１１において、対象燃料に添加される添加成分種を設定する。
本実施形態においては、対象燃料のＰＲＦ７０に対して、トルエンを添加成分種として設定する。

【0032】

次に、ステップＳ１２において、対象燃料の燃料成分種とステップＳ１１で設定した添加成分種との混合割合を設定する。
本実施形態においては、燃料成分種としてイソオクタンとノルマルヘプタンとをイソオクタン：ノルマルヘプタン＝７０：３０（ｖｏｌ．％）の混合割合で含む対象燃料ＰＲＦ７０に対して、アロマ系成分の代表としてトルエンを５０ｖｏｌ．％添加した混合燃料を想定する。この場合、混合割合は、イソオクタン：ノルマルヘプタン：トルエン＝３５：１５：５０（ｖｏｌ．％）である。

【0033】

次に、ステップＳ１３において、ステップＳ１２で設定した混合割合の添加成分種を、燃料成分種の１つに置換する。
本実施形態においては、トルエンをイソオクタンに置換することにより、イソオクタン：ノルマルヘプタン＝８５：１５（ｖｏｌ．％）の混合割合の混合燃料を想定する。

【0034】

次に、ステップＳ１４において、添加成分種をステップＳ１２で設定した混合割合で添加した混合燃料の燃焼反応における正確な低温酸化反応時期を目標低温酸化反応時期τ₁^mixtgtとして取得すると共に、その混合燃料の燃焼反応における正確な熱着火時期を目標熱着火時期τ₂^mixtgtとして取得する。
本実施形態においては、ＭＦＣスキームを使用して、混合割合がイソオクタン：ノルマルヘプタン：トルエン＝３５：１５：５０（ｖｏｌ．％）の混合燃料について圧力６ＭＰａ、空気過剰率（λ）＝１の条件で初期温度７００Ｋから１４００Ｋまで５０Ｋ毎に化学反応計算を行うことにより、各温度における目標低温酸化反応時期τ₁^mixtgt及び目標熱着火時期τ₂^mixtgtを取得する。

【0035】

次に、ステップＳ１５において、図２のステップＳ５で選択した各化学反応式の内、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式の化学反応速度定数を、ステップＳ１３において添加成分種を燃料成分種の１つに置換した混合燃料について図２のフローにおいて最適化した簡略化化学反応スキーム１８を使用した燃焼シミュレーションにより算出される低温酸化反応時期τ₁及び熱着火時期τ₂が、ステップＳ１４で取得した目標低温酸化反応時期τ₁^mixtgt及び目標熱着火時期τ₂^mixtgtにそれぞれ一致するように再度最適化する。
本実施形態においては、図２のステップＳ７において最適化した化学反応式Ｒ６及びＲ２１の化学反応速度定数を、以下の式（２）により表される評価関数Ｅ_fの最小値を得るように、ＰＳＯ法により再度最適化を行う。

【数4】

ステップＳ１５の後、混合燃料の燃焼シミュレーションに使用する化学反応速度定数を決定するフローを終了する。

【0036】

簡略化化学反応スキーム１８に含まれる化学反応式Ｒ６及びＲ２１の化学反応速度定数を図５のステップＳ１５において最適化した後の簡略化化学反応スキーム１８を使用して燃焼シミュレーションを行うことにより、図６に示すように、ＰＲＦ７０に対してトルエンを５０ｖｏｌ．％添加した混合燃料についても、圧力６ＭＰａ、空気過剰率（λ）１、初期温度７００Ｋから１４００Ｋまでの温度範囲の条件において、ＭＦＣスキームによるシミュレーション結果に良く一致することが分かる。即ち、トルエン混合により低温酸化反応が抑制されるため、低温度領域のみで着火遅れ時間が延長するが、化学反応式Ｒ６及びＲ２１の化学反応速度定数の最適化によってこの特性をよく再現可能である。

【0037】

図５に示したフローを、様々な添加成分種（例えばトルエンや１−ペンテン）毎に、対象燃料の燃料成分種と添加成分種との様々な混合割合について実行することにより、各混合割合に対応して最適化した化学反応速度定数マップ２０を得ることができる。
例えば、ＰＲＦ７０に対するトルエンの混合割合を０ｖｏｌ．％から５０ｖｏｌ．％まで１０ｖｏｌ．％ずつ変化させ、それぞれの混合割合について化学反応式Ｒ６及びＲ２１の化学反応速度定数を最適化した場合、化学反応速度定数の変化率Ｒ（Ｔｏｌ）を以下の式（３）により表すことができる。ここで、Ａ_Tolは図５のステップＳ１５で最適化した化学反応速度定数であり、Ａ₀は初期簡略化化学反応スキームにおける化学反応速度定数である。
Ａ_Tol＝Ｒ（Ｔｏｌ）×Ａ₀・・・（３）

【0038】

図７に示すように、トルエンの混合割合が１０ｖｏｌ．％でＲ６、Ｒ２１ともに化学反応速度定数は一旦大きくなり、混合割合が増加するにつれて次第に小さくなる。ここではトルエンをイソオクタンで置換して計算しているため、Ｒ２１の化学反応速度定数の変化率が１の場合には、トルエンとイソオクタンのＨＣＨＯ、Ｈ₂Ｏ₂生成能が等しいということである。よって、１０ｖｏｌ．％混合では、トルエン混合によりＨＣＨＯ、Ｈ₂Ｏ₂の生成が促進されたことを表している。さらに、トルエン混合割合が増加するにつれて、ＨＣＨＯ、Ｈ₂Ｏ₂の生成が抑制されていくことがわかる。また、Ｒ６の化学反応速度定数の変化率が１の場合には、トルエンが混合されてもノルマルヘプタンのＨＣＨＯ、Ｈ₂Ｏ₂生成能が変化しないことを表すが、このＲ６の化学反応速度定数も変化することから、トルエン混合がノルマルヘプタンのＨＣＨＯ、Ｈ₂Ｏ₂生成にも影響することがわかる。

【0039】

次に、図８及び図９により、本発明の実施形態によるエンジンの燃焼シミュレーションについて説明する。
図８は本発明の実施形態によるエンジンの燃焼シミュレーション方法を示すフローチャートであり、図９はＰＲＦ７０にトルエン及び１−ペンテンを添加した混合燃料について本発明の実施形態による燃焼シミュレーションにより算出した着火遅れ時間とＭＦＣスキームを使用した燃焼シミュレーションにより算出した着火遅れ時間とを比較した線図である。

【0040】

まず、図８に示すフローは、図５のフローで対象燃料に添加成分種を添加した混合燃料の燃焼シミュレーションに使用可能な簡略化化学反応スキーム１８を作成した後に実行される。
図８に示すように、エンジンの燃焼シミュレーションが開始されると、ステップＳ２１において、入力装置２により燃焼シミュレーションの設定データ（エンジンの寸法、形状、物性値、対象燃料の燃料成分種、添加成分種、混合割合、燃焼反応の初期温度、空気過剰率、圧力等）が入力される。

【0041】

次に、ステップＳ２２において、プリプロセッサ１０は、ステップＳ２１で入力された設定データに基づき、燃焼シミュレーションに用いる物理モデルを構築する。

【0042】

次に、ステップＳ２３において、化学反応解析部１２は、化学反応データベース８から簡略化化学反応スキーム１８を読み込む。次いで、ステップＳ２４において、化学反応解析部１２は、化学反応データベース８の化学反応速度定数マップ２０を参照し、ステップＳ２１で入力された対象燃料の燃料成分種、添加成分種、混合割合等に対応する化学反応速度定数を読み込む。

【0043】

次に、ステップＳ２５において、化学反応解析部１２は、対象燃料に複数の添加成分種が添加されている場合、その添加成分種の混合割合に応じた化学反応速度定数を算出する。
本実施形態においては、例えば、ＰＲＦ７０に対してトルエンを３０ｖｏｌ．％、１−ペンテンを２０ｖｏｌ．％混合した燃料（ノルマルヘプタン:イソオクタン:トルエン:１−ペンテン＝１５：３５：３０：２０ｖｏｌ．％）を想定した場合、化学反応解析部１２は、各添加成分種をイソオクタンで置換し、化学反応式Ｒ６及びＲ２１の化学反応速度定数を、それぞれトルエンと１−ペンテンの混合割合に応じて以下の式（４）のように加重平均した値を算出する。ここで、Ｒ（Ｔｏｌ）、Ｒ（Ｐｅｎ）はそれぞれトルエン、１−ペンテンを単独で対象燃料に添加した場合の化学反応速度定数の変化率であり、Ｔｏｌ％、Ｐｅｎ％はそれぞれトルエン、１−ペンテンの混合割合である。

【数5】

【0044】

次に、ステップＳ２６において、化学反応解析部１２及び３次元流体解析部１４は、ステップＳ２２において構築した物理モデル、ステップＳ２３で読み込んだ簡略化化学反応スキーム１８、ステップＳ２４で読み込んだ化学反応速度定数、及び、ステップＳ２５で算出した化学反応速度定数を用いて、燃焼シミュレーションを実行する。

【0045】

次に、ステップＳ２７において、ポストプロセッサ１６は、ステップＳ２６の燃焼シミュレーション結果を所定の出力形式に変換する処理を行い、出力装置６から出力させる。その後、エンジンの燃焼シミュレーションを終了する。

【0046】

ステップＳ２５において、対象燃料に添加された複数の添加成分種の混合割合に応じた化学反応速度定数を算出し、その化学反応速度定数を使用して燃焼シミュレーションを実行することにより、複数の添加成分種を添加した混合燃料についても精度よく燃焼シミュレーションを行うことができる。
例えば、ＰＲＦ７０に対してトルエンを３０ｖｏｌ．％、１−ペンテンを２０ｖｏｌ．％混合した燃料（ノルマルヘプタン:イソオクタン:トルエン:１−ペンテン＝１５：３５：３０：２０ｖｏｌ．％）について、各添加成分種の混合割合に応じて算出した化学反応速度定数を使用して着火遅れ時間を算出した場合、図９に示すように、ＭＦＣスキームによるシミュレーション結果に良く一致することが分かる。

【0047】

次に、本発明の実施形態のさらなる変形例を説明する。
まず、上述した実施形態では、初期簡略化化学反応スキームとして、ＰＲＦの着火現象の再現に適した化学反応式Ｒ１乃至Ｒ３８を使用する場合を説明したが、他の化学反応スキームを使用してもよい。
例えば、Ａｍａｒ Ｐａｔｅｌらによって提案されたＨＣＣＩエンジンシミュレーションのための化学反応スキーム「ＥＲＣ Ｍｅｃｈａｎｉｓｍ」を、初期簡略化化学反応スキームとして使用してもよい。このＥＲＣ Ｍｅｃｈａｎｉｓｍは、以下のｒ１乃至ｒ５２の化学反応式を含む。
ｎ−Ｃ₇Ｈ₁₆＋Ｈ＝Ｃ₇Ｈ₁₅−２＋Ｈ₂ （ｒ１）
ｎ−Ｃ₇Ｈ₁₆＋ＯＨ＝Ｃ₇Ｈ₁₅−２＋Ｈ₂Ｏ （ｒ２）
ｎ−Ｃ₇Ｈ₁₆＋ＨＯ₂＝Ｃ₇Ｈ₁₅−２＋Ｈ₂Ｏ₂ （ｒ３）
ｎ−Ｃ₇Ｈ₁₆＋Ｏ₂＝Ｃ₇Ｈ₁₅−２＋ＨＯ₂ （ｒ４）
Ｃ₇Ｈ₁₅−２＋Ｏ₂＝Ｃ₇Ｈ₁₅Ｏ₂ （ｒ５）
Ｃ₇Ｈ₁₅Ｏ₂＋Ｏ₂＝Ｃ₇ＫＥＴ₁₂＋ＯＨ （ｒ６）
Ｃ₇ＫＥＴ₁₂＝Ｃ₅Ｈ₁₁ＣＯ＋ＣＨ₂Ｏ＋ＯＨ （ｒ７）
Ｃ₅Ｈ₁₁ＣＯ＝Ｃ₂Ｈ₄＋Ｃ₃Ｈ₇＋ＣＯ （ｒ８）
Ｃ₇Ｈ₁₅−２＝Ｃ₂Ｈ₅＋Ｃ₂Ｈ₄＋Ｃ₃Ｈ₆ （ｒ９）
Ｃ₃Ｈ₇＝Ｃ₂Ｈ₄＋ＣＨ₃ （ｒ１０）
Ｃ₃Ｈ₇＝Ｃ₃Ｈ₆＋Ｈ （ｒ１１）
Ｃ₃Ｈ₆＋ＣＨ₃＝Ｃ₃Ｈ₅＋ＣＨ₄ （ｒ１２）
Ｃ₃Ｈ₅＋Ｏ₂＝Ｃ₃Ｈ₄＋ＨＯ₂ （ｒ１３）
Ｃ₃Ｈ₄＋ＯＨ＝Ｃ₂Ｈ₃＋ＣＨ₂Ｏ （ｒ１４）
Ｃ₃Ｈ₄＋ＯＨ＝Ｃ₂Ｈ₄＋ＨＣＯ （ｒ１５）
ＣＨ₃＋ＨＯ₂＝ＣＨ₃Ｏ＋ＯＨ （ｒ１６）
ＣＨ₃＋ＯＨ＝ＣＨ₂＋Ｈ₂Ｏ （ｒ１７）
ＣＨ₂＋ＯＨ＝ＣＨ₂Ｏ＋Ｈ （ｒ１８）
ＣＨ₂＋Ｏ₂＝ＨＣＯ＋ＯＨ （ｒ１９）
ＣＨ₂＋Ｏ₂＝ＣＯ₂＋Ｈ₂ （ｒ２０）
ＣＨ₂＋Ｏ₂＝ＣＯ＋Ｈ₂Ｏ （ｒ２１）
ＣＨ₂＋Ｏ₂＝ＣＨ₂Ｏ＋Ｏ （ｒ２２）
ＣＨ₂＋Ｏ₂＝ＣＯ₂＋Ｈ＋Ｈ （ｒ２３）
ＣＨ₂＋Ｏ₂＝ＣＯ＋ＯＨ＋Ｈ （ｒ２４）
ＣＨ₃Ｏ＋ＣＯ＝ＣＨ₃＋ＣＯ₂ （ｒ２５）
ＣＯ＋ＯＨ＝ＣＯ₂＋Ｈ （ｒ２６）
Ｏ＋ＯＨ＝Ｏ₂＋Ｈ （ｒ２７）
Ｈ＋ＨＯ₂＝ＯＨ＋ＯＨ （ｒ２８）
ＯＨ＋ＯＨ＝Ｏ＋Ｈ₂Ｏ （ｒ２９）
Ｈ＋Ｏ₂＋Ｍ＝ＨＯ₂＋Ｍ （ｒ３０）
Ｈ₂Ｏ₂＋Ｍ＝ＯＨ＋ＯＨ＋Ｍ （ｒ３１）
Ｈ₂＋ＯＨ＝Ｈ₂Ｏ＋Ｈ （ｒ３２）
ＨＯ₂＋ＨＯ₂＝Ｈ₂Ｏ₂＋Ｏ₂ （ｒ３３）
ＣＨ₂Ｏ＋ＯＨ＝ＨＣＯ＋Ｈ₂Ｏ （ｒ３４）
ＣＨ₂Ｏ＋ＨＯ₂＝ＨＣＯ＋Ｈ₂Ｏ₂ （ｒ３５）
ＨＣＯ＋Ｏ₂＝ＨＯ₂＋ＣＯ （ｒ３６）
ＨＣＯ＋Ｍ＝Ｈ＋ＣＯ＋Ｍ （ｒ３７）
ＣＨ₃＋ＣＨ₃Ｏ＝ＣＨ₄＋ＣＨ₂Ｏ （ｒ３８）
Ｃ₂Ｈ₄＋ＯＨ＝ＣＨ₂Ｏ＋ＣＨ₃ （ｒ３９）
Ｃ₂Ｈ₄＋ＯＨ＝Ｃ₂Ｈ₃＋Ｈ₂Ｏ （ｒ４０）
Ｃ₂Ｈ₃＋Ｏ₂＝ＣＨ₂Ｏ＋ＨＣＯ （ｒ４１）
Ｃ₂Ｈ₃＋ＨＣＯ＝Ｃ₂Ｈ₄＋ＣＯ （ｒ４２）
Ｃ₂Ｈ₅＋Ｏ₂＝Ｃ₂Ｈ₄＋ＨＯ₂ （ｒ４３）
ＣＨ₄＋Ｏ₂＝ＣＨ₃＋ＨＯ₂ （ｒ４４）
ＯＨ＋ＨＯ₂＝Ｈ₂Ｏ＋Ｏ₂ （ｒ４５）
ＣＨ₃＋Ｏ₂＝ＣＨ₂Ｏ＋ＯＨ （ｒ４６）
ＣＨ₄＋Ｈ＝ＣＨ₃＋Ｈ₂ （ｒ４７）
ＣＨ₄＋ＯＨ＝ＣＨ₃＋Ｈ₂Ｏ （ｒ４８）
ＣＨ₄＋Ｏ＝ＣＨ₃＋ＯＨ （ｒ４９）
ＣＨ₄＋ＨＯ₂＝ＣＨ₃＋Ｈ₂Ｏ₂ （ｒ５０）
ＣＨ₄＋ＣＨ₂＝ＣＨ₃＋ＣＨ₃ （ｒ５１）
Ｃ₃Ｈ₆＝Ｃ₂Ｈ₃＋ＣＨ₃ （ｒ５２）

【0048】

また、上述した実施形態では、対象燃料の燃焼反応における目標低温酸化反応時期τ₁^tgt及び目標熱着火時期τ₂^tgtを取得するために、ＭＦＣスキームを使用する場合を説明したが、これとは異なる化学反応スキームを使用してもよく、あるいは、実験により目標低温酸化反応時期τ₁^tgt及び目標熱着火時期τ₂^tgtを取得してもよい。

【0049】

また、上述した実施形態では、化学反応速度定数の最適化を行う際に、ＰＳＯ法を使用する場合を説明したが、これとは異なる最適化手法を使用してもよい。

【0050】

次に、上述した本発明の実施形態及び本発明の実施形態の変形例によるエンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法及びエンジンの燃焼シミュレーション方法の作用効果を説明する。

【0051】

まず、燃焼シミュレーションに使用する化学反応式の中から、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式を少なくとも１つ選択すると共に、熱着火時期の反応温度に対して相対的に感度の高い化学反応式を少なくとも１つ選択し、それらの選択した化学反応式の化学反応速度定数を、燃焼シミュレーションにより算出される低温酸化反応時期及び熱着火時期が所定の目標低温酸化反応時期及び目標熱着火時期にそれぞれ一致するように最適化するので、高い精度で低温酸化反応時期及び熱着火時期を算出することが可能な簡略化化学反応スキーム１８を作成することができ、これにより、エンジンの燃焼反応を少ない計算負荷で精度よく再現することができる。

【0052】

また、対象燃料に添加成分種を添加した混合燃料については、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式の化学反応速度定数を、対象燃料に添加された添加成分種を燃料成分種の１つに置換した場合に燃焼シミュレーションにより算出される低温酸化反応時期及び熱着火時期が混合燃料における目標低温酸化反応時期及び目標熱着火時期にそれぞれに一致するように再度最適化するので、混合燃料についても高い精度で低温酸化反応時期及び熱着火時期を算出することが可能な簡略化化学反応スキーム１８を作成することができ、これにより、様々な混合燃料についてエンジンの燃焼反応を少ない計算負荷で精度よく再現することができる。

【0053】

特に、対象燃料は、燃料成分種としてイソオクタンとノルマルヘプタンとを含む混合燃料ＰＲＦであるので、ＰＲＦについてエンジンの燃焼反応を少ない計算負荷で精度よく再現することができる。

【0054】

また、燃焼シミュレーションに使用する化学反応式として、上述のＲ１乃至Ｒ３８の化学反応式を設定するので、ＰＲＦについてエンジンの燃焼反応を少ない計算負荷で一層精度よく再現することができる。

【0055】

また、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式Ｒ６及びＲ２１を選択すると共に、熱着火時期の反応温度に対して相対的に感度の高い化学反応式Ｒ３７を選択して化学反応速度定数を最適化するので、低温酸化反応時期及び熱着火時期を精度よく計算することが可能な簡略化化学反応スキーム１８を作成できる。

【0056】

特に、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式及び、熱着火時期の反応温度に対して相対的に感度の高い化学反応式のぞれぞれの化学反応速度定数を、ＰＳＯ法により最適化するので、高い精度で燃焼シミュレーションを実行することが可能な簡略化化学反応スキーム１８を効率よく作成することができる。

【0057】

また、対象燃料に添加成分種を添加した混合燃料については、低温酸化反応中のＨＣＨＯ及びＨ₂Ｏ₂の生成量に対して相対的に感度の高い化学反応式の化学反応速度定数を、ＰＳＯ法により再度最適化するので、様々な混合燃料について高い精度で燃焼シミュレーションを実行することが可能な簡略化化学反応スキーム１８を効率よく作成することができる。

【符号の説明】

【0058】

１ エンジンの燃焼シミュレーションシステム
２ 入力装置
４ 燃焼シミュレーション装置
６ 出力装置
８ 化学反応データベース
１０ プリプロセッサ
１２ 化学反応解析部
１４ ３次元流体解析部
１６ ポストプロセッサ
１８ 簡略化化学反応スキーム
２０ 化学反応速度定数マップ

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】