JP5051659B2 - 反応メカニズム解析装置及び反応メカニズム解析プログラム - Google Patents
反応メカニズム解析装置及び反応メカニズム解析プログラム Download PDFInfo
- Publication number
- JP5051659B2 JP5051659B2 JP2008287099A JP2008287099A JP5051659B2 JP 5051659 B2 JP5051659 B2 JP 5051659B2 JP 2008287099 A JP2008287099 A JP 2008287099A JP 2008287099 A JP2008287099 A JP 2008287099A JP 5051659 B2 JP5051659 B2 JP 5051659B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- reaction
- temperature distribution
- elementary
- temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000006243 chemical reaction Methods 0.000 title claims description 264
- 230000007246 mechanism Effects 0.000 title claims description 96
- 238000004458 analytical method Methods 0.000 claims description 116
- 238000004364 calculation method Methods 0.000 claims description 93
- 238000004088 simulation Methods 0.000 claims description 92
- 230000008859 change Effects 0.000 claims description 82
- 239000013067 intermediate product Substances 0.000 claims description 49
- 238000013500 data storage Methods 0.000 claims description 44
- 238000002474 experimental method Methods 0.000 claims description 28
- 239000000203 mixture Substances 0.000 claims description 27
- 230000006870 function Effects 0.000 claims description 25
- 239000000446 fuel Substances 0.000 claims description 21
- 239000000126 substance Substances 0.000 claims description 21
- 238000012546 transfer Methods 0.000 claims description 13
- 239000007800 oxidant agent Substances 0.000 claims description 12
- 239000000463 material Substances 0.000 claims description 7
- 238000004134 energy conservation Methods 0.000 claims description 6
- 230000014509 gene expression Effects 0.000 claims description 5
- 239000007789 gas Substances 0.000 description 118
- 238000002485 combustion reaction Methods 0.000 description 17
- 238000005259 measurement Methods 0.000 description 14
- 230000015572 biosynthetic process Effects 0.000 description 13
- 238000000034 method Methods 0.000 description 11
- 230000001590 oxidative effect Effects 0.000 description 10
- 238000012545 processing Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 238000004891 communication Methods 0.000 description 6
- 230000003247 decreasing effect Effects 0.000 description 5
- 239000007788 liquid Substances 0.000 description 4
- 230000004913 activation Effects 0.000 description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 3
- 239000001301 oxygen Substances 0.000 description 3
- 229910052760 oxygen Inorganic materials 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 239000000047 product Substances 0.000 description 3
- 238000009529 body temperature measurement Methods 0.000 description 2
- 239000000470 constituent Substances 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 239000003350 kerosene Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000000329 molecular dynamics simulation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Physical Or Chemical Processes And Apparatus (AREA)
Description
本発明は、複数種類の初期物質同士の化学反応を構成する素反応を解析する反応メカニズム解析装置、反応メカニズム解析方法及び反応メカニズム解析プログラムに関するものである。
ある初期物質同士が化学反応によって最終生成物となる場合には、その化学反応の過程において中間生成物が生成される。そして、ある初期物質が最終生成物となる化学反応は、上記中間生成物が生成される複数の素反応の組み合わせによって記述される。この素反応の数は非常に多く、例えば最も単純な水素と酸素の燃焼反応においても、３７種の素反応を有している。これが灯油等の分子量が大きな物質の燃焼反応では、生成される中間生成物の種類も非常に多くなり、５０００以上の素反応を有する場合も珍しくない。
近年は、コンピュータの進歩により、様々な現象をシミュレーションによって解析することが可能となっている。そして、このようなシミュレーションによって化学反応等を含む実現象を正確に解析するためには、上述の素反応によって生成される中間生成物の濃度変化や混合ガス（初期物質同士のみが混合された予混合ガスのみからなるガス、あるいは当該予混合ガスと中間生成物とが混合したガスの両方を含む）の温度変化をコンピュータに組み込むことによって、素反応を考慮してシミュレーションを行う必要がある。このため、シミュレーションに先立ち、各素反応によって生成される中間生成物の濃度変化や混合ガスの温度変化を明らかにする、いわゆるフルケミストリ解析が行われる。
例えば、エンジン等の燃焼機関における流場をＣＦＤ（Computational Fluid Dynamics）解析する場合には、より正確なＣＦＤ解析を可能とするために、ＣＦＤ解析に先立ち、燃料が燃焼反応する際に生成される中間生成物の濃度変化及び混合ガスの温度変化を明らかにするためのフルケミストリ解析を行う。
例えば、エンジン等の燃焼機関における流場をＣＦＤ（Computational Fluid Dynamics）解析する場合には、より正確なＣＦＤ解析を可能とするために、ＣＦＤ解析に先立ち、燃料が燃焼反応する際に生成される中間生成物の濃度変化及び混合ガスの温度変化を明らかにするためのフルケミストリ解析を行う。
一般的にフルケミストリ解析を行う場合には、実験によって燃料が燃焼反応する際の温度分布を実験データとして取得し、予め分子動力学法（特許文献１参照）等によって算出された素反応式に関連付けられたパラメータを実験データに合うように調整することによって行われる。
より詳細には、例えば非特許文献１に記載されたいわゆるＡＳＴＭ（American Standard of Testing Method）法等の何らかの実験手法を用いて燃料が燃焼反応する際の温度分布を測定する。続いて、実験空間を複数の領域に分割し、予め与えられた複数の素反応式及びこれらの素反応式に関連付けられたパラメータを用いて各領域の温度をシミュレーションにて算出することによって温度分布を取得する。そして、シミュレーションによって得られた温度分布と実験データが示す温度分布との差が許容可能範囲内となるまで、繰り返し素反応式に関連付けられたパラメータ（例えば頻度係数及び活性化エネルギ）を調整し、最終的に実験データと一致するシミュレーションデータを導き出せるパラメータの解を取得する。
そして、上述のように導き出されたパラメータを用いることによって、任意の時間における中間生成物の濃度や混合ガスの温度が算出可能となり、より正確なＣＦＤ解析等のシミュレーションを行うことが可能となる。
特開２００５−３４８４３号公報
板垣晴彦，産業安全研究所研究報告 ＮＩＩＳ−ＲＲ−２００２（２００３） 「可燃性液体の発火温度の圧力依存性について」，Research Reports of the National Institute of Industrial Safety, NIIS-RR-2002(2003) UDC 614.841.41+543.87:662.612.1「The Dependence of Initial Pressure on the Autoignition Temperature of Flammable Liquid」（www.jniosh.go.jp/publication/RR/pdf/RR-2002-06.pdf）
より詳細には、例えば非特許文献１に記載されたいわゆるＡＳＴＭ（American Standard of Testing Method）法等の何らかの実験手法を用いて燃料が燃焼反応する際の温度分布を測定する。続いて、実験空間を複数の領域に分割し、予め与えられた複数の素反応式及びこれらの素反応式に関連付けられたパラメータを用いて各領域の温度をシミュレーションにて算出することによって温度分布を取得する。そして、シミュレーションによって得られた温度分布と実験データが示す温度分布との差が許容可能範囲内となるまで、繰り返し素反応式に関連付けられたパラメータ（例えば頻度係数及び活性化エネルギ）を調整し、最終的に実験データと一致するシミュレーションデータを導き出せるパラメータの解を取得する。
そして、上述のように導き出されたパラメータを用いることによって、任意の時間における中間生成物の濃度や混合ガスの温度が算出可能となり、より正確なＣＦＤ解析等のシミュレーションを行うことが可能となる。
しかしながら、従来の実験方法では、実験空間が三次元性を持っており、例えばＡＳＴＭ法では、容量が５００ｍｌの丸底フラスコの内部が実験空間とされている。このため、正確なフルケミストリ解析を行うためには実験空間を数万の領域に分割する必要がある。そして、燃料として上述の灯油を想定した場合には、数万の分割領域の各々に対して５０００以上の素反応式の計算を行い、さらに最終的に実験データと一致するシミュレーションデータが得られるまで素反応式に関連付けられたパラメータを繰り返し調整しながら同様の計算を繰り返し行う必要がある。このため、計算能力に優れたスーパーコンピュータを用いたとしても、現実的な時間内にて実験データと一致するシミュレーションデータを導き出せるパラメータの解を取得することは難しい。
なお、実験空間の位置を、温度の従属変数としてとらえ、中間生成物の濃度変化を算出する計算を空間零次元化して計算量を減少させることができる。
しかしながら、例えばＡＳＴＭ法による測定方法では、大きな容積をもつ丸底フラスコに可燃性液体を滴下させて着火させるため、丸底フラスコ内の温度分布が不均一となり、測定される温度分布が実験環境における不確定要因に影響を受ける場合がある。また、丸底フラスコ内に滴下された可燃性液体は、着火前に蒸発して混合気となる。そして、着火温度は混合気の濃度に依存して変化するが、ＡＳＴＭ法では、混合気の濃度が一定である保証はない。つまり、実験空間が三次元性を持つ実験方法にて測定された温度分布は、少なからず実験環境における不確定要因の影響を受けた値である。
このため、実験空間が三次元性を持つ実験方法にて得られた実験データに合わせてシミュレーションデータを導き出せる素反応式に関連付けられたパラメータの解を取得すると、当該解は、少なからず実験環境における不確定要因の影響を含むものとなる。よって、中間生成物の濃度変化や混合ガスの温度変化も実験環境における不確定要因の影響を含むものとなり、正確なフルケミストリ解析を行うことができない。したがって、例えば、得られた素反応式に関連付けられたパラメータの解を用いてＣＦＤ解析を行った場合には、正確なＣＦＤ解析結果が得られない虞がある。
しかしながら、例えばＡＳＴＭ法による測定方法では、大きな容積をもつ丸底フラスコに可燃性液体を滴下させて着火させるため、丸底フラスコ内の温度分布が不均一となり、測定される温度分布が実験環境における不確定要因に影響を受ける場合がある。また、丸底フラスコ内に滴下された可燃性液体は、着火前に蒸発して混合気となる。そして、着火温度は混合気の濃度に依存して変化するが、ＡＳＴＭ法では、混合気の濃度が一定である保証はない。つまり、実験空間が三次元性を持つ実験方法にて測定された温度分布は、少なからず実験環境における不確定要因の影響を受けた値である。
このため、実験空間が三次元性を持つ実験方法にて得られた実験データに合わせてシミュレーションデータを導き出せる素反応式に関連付けられたパラメータの解を取得すると、当該解は、少なからず実験環境における不確定要因の影響を含むものとなる。よって、中間生成物の濃度変化や混合ガスの温度変化も実験環境における不確定要因の影響を含むものとなり、正確なフルケミストリ解析を行うことができない。したがって、例えば、得られた素反応式に関連付けられたパラメータの解を用いてＣＦＤ解析を行った場合には、正確なＣＦＤ解析結果が得られない虞がある。
本発明は、上述する問題点に鑑みてなされたもので、計算負荷が小さくかつ正確なフルケミストリ解析を可能とすることを目的とする。
本発明は、上記課題を解決するための手段として、以下の構成を採用する。
第１の発明は、複数種類の初期物質同士の化学反応を構成する素反応を解析する反応メカニズム解析装置であって、実験により得られた上記化学反応の際の空間一次元の温度分布を実験データとして記憶する実験データ記憶手段と、上記化学反応を記述する複数の素反応式及び該素反応式に関連付けられたパラメータを素反応データとして記憶する素反応データ記憶手段と、上記素反応データからシミュレーションデータを算出するための空間零次元で時間一次元化された方程式を演算データとして記憶する演算データ記憶手段と、上記素反応データ及び上記空間零次元で時間一次元化された方程式を用いて上記シミュレーションデータを算出する算出部と、上記シミュレーションデータが上記実験データを導出可能であるかを判定する判定部と、上記素反応データに含まれる上記パラメータを調整可能なパラメータ調整手段と、上記判定部において上記シミュレーションデータが上記実験データを導出可能であると判定されるまで、上記パラメータ調整手段による上記パラメータの調整及び調整後の上記パラメータを用いた上記算出部による上記シミュレーションデータの算出を繰り返し実行させる制御手段とを備えるという構成を採用する。
第２の発明は、上記第１の発明において、上記空間零次元で時間一次元化された方程式が、中間生成物の濃度変化を示す項を含む質量保存の式と、素反応によって生成される中間生成物と上記初期物質との混合体の内部エネルギを示す項、上記素反応により移動するエネルギを示す項、及び上記混合体と外部との熱伝達により移動するエネルギを示す項を含むエネルギ保存の式とを有するという構成を採用する。
第３の発明は、上記第１または第２の発明において、上記算出部が、上記素反応によって生成される中間生成物の濃度変化及び上記初期物質と上記中間生成物との混合体の温度変化とを上記シミュレーションデータとして算出し、上記判定部が、上記混合体の温度変化を空間一次元化することによって空間一次元のシミュレーション温度分布を算出し、該シミュレーション温度分布と上記実験データとして記憶された上記温度分布とを比較し、上記シミュレーション温度分布が上記実験データとして記憶された温度分布との差が許容範囲内である場合に上記シミュレーションデータが上記実験データを導出可能であると判定するという構成を採用する。
第４の発明は、上記第１〜第３の発明において、上記実験データとして記憶される空間一次元の温度分布が、複数の上記初期物質が混合された予混合ガスを、常温における消炎径よりも小さな直径に設定された内部流路を備えると共に流れ方向に上記予混合ガスの着火温度以上まで連続的に昇温するように加熱された管に、火炎形成位置が安定しかつ上記予混合ガスの流速に依存しなくなる流量で供給した際に得られる温度分布であるという構成を採用する。
第５の発明は、上記第１〜第４いずれかの発明において、解析された素反応を基に解析を行う際の解析条件を記憶する解析条件記憶手段と、上記化学反応を記述する複数の素反応式を上記解析条件に基づいて選択する選択手段とを備え、上記制御手段が、上記判定部において上記シミュレーションデータが上記実験データを導出可能であると判定された場合に上記選択手段による上記素反応式の選択を実行させるという構成を採用する。
第６の発明は、複数種類の初期物質同士の化学反応を構成する素反応を解析する反応メカニズム解析方法であって、上記化学反応を記述する複数の素反応式及び該素反応式に関連付けられたパラメータを有する素反応データ及び該素反応データからシミュレーションデータを算出するための空間零次元で時間一次元化された方程式を有する演算データを用いてシミュレーションデータを算出する算出工程と、上記シミュレーションデータが、実験により得られた上記化学反応の際の空間一次元の温度分布を有する実験データを導出可能であるかを判定する判定工程と、上記素反応データに含まれる上記パラメータを調整可能なパラメータ調整工程とを有し、上記判定工程において上記シミュレーションデータが上記実験データを導出可能であると判定されるまで、上記パラメータ調整工程及び調整後の上記パラメータを用いた上記算出工程を繰り返し行うという構成を採用する。
第７の発明は、上記第６の発明において、上記空間零次元で時間一次元化された方程式が、中間生成物の濃度変化を示す項を含む質量保存の式と、素反応によって生成される中間生成物と上記初期物質との混合体の内部エネルギを示す項、上記素反応により移動するエネルギを示す項、及び上記混合体と外部との熱伝達により移動するエネルギを示す項を含むエネルギ保存の式とを有するという構成を採用する。
第８の発明は、上記第６または第７の発明において、上記算出工程にて、上記素反応によって生成される中間生成物の濃度変化及び上記初期物質と上記中間生成物との混合体の温度変化とを上記シミュレーションデータとして算出し、上記判定工程にて、上記混合体の温度変化を空間一次元化することによって空間一次元のシミュレーション温度分布を算出し、該シミュレーション温度分布と上記実験データとして記憶された上記温度分布とを比較し、上記シミュレーション温度分布が上記実験データとして記憶された温度分布との差が許容範囲内である場合に上記シミュレーションデータが上記実験データを導出可能であると判定するという構成を採用する。
第９の発明は、上記第６〜第８いずれかの発明において、複数の上記初期物質が混合された予混合ガスを、常温における消炎径よりも小さな直径に設定された内部流路を備えると共に流れ方向に上記予混合ガスの着火温度以上まで連続的に昇温するように加熱された管に、火炎形成位置が安定しかつ上記予混合ガスの流速に依存しなくなる流量で供給した際に得られる温度分布を上記化学反応の際の空間一次元の温度分布として取得する実験工程を有するという構成を採用する。
第１０の発明は、上記第６〜第９いずれかの発明において、解析された素反応を基に解析を行う際の解析条件に基づいて上記化学反応を記述する複数の素反応式を選択する選択工程を有し、上記判定工程において上記シミュレーションデータが上記実験データを導出可能であると判定された場合に上記選択工程を行うという構成を採用する。
第１１の発明は、反応メカニズム解析プログラムであって、コンピュータを、実験により得られた上記化学反応の際の空間一次元の温度分布を実験データとして記憶する実験データ記憶手段と、上記化学反応を記述する複数の素反応式及び該素反応式に関連付けられたパラメータを素反応データとして記憶する素反応データ記憶手段と、上記素反応データからシミュレーションデータを算出するための空間零次元で時間一次元化された方程式を演算データとして記憶する演算データ記憶手段と、上記素反応データ及び上記空間零次元で時間一次元化された方程式を用いて上記シミュレーションデータを算出する算出部と、上記シミュレーションデータが上記実験データを導出可能であるかを判定する判定部と、上記素反応データに含まれる上記パラメータを調整可能なパラメータ調整手段と、上記判定部において上記シミュレーションデータが上記実験データを導出可能であると判定されるまで、上記パラメータ調整手段による上記パラメータの調整及び調整後の上記パラメータを用いた上記算出部による上記シミュレーションデータの算出を繰り返し実行させる制御手段として機能させるという構成を採用する。
第１２の発明は、上記第１１の発明において、上記コンピュータを上記算出部として機能させる場合に、中間生成物の濃度変化を示す項を含む質量保存の式と、素反応によって生成される中間生成物と上記初期物質との混合体の内部エネルギを示す項、上記素反応により移動するエネルギを示す項、及び上記混合体と外部との熱伝達により移動するエネルギを示す項を含むエネルギ保存の式とを上記空間零次元で時間一次元化された方程式として用いて上記シミュレーションデータを算出させるという構成を採用する。
第１３の発明は、上記第１１または第１２の発明において、上記コンピュータを上記算出部として機能させる場合に、上記素反応によって生成される中間生成物の濃度変化及び上記初期物質と上記中間生成物との混合体の温度変化とを上記シミュレーションデータとして算出させ、上記コンピュータを上記判定部として機能させる場合に、上記混合体の温度変化を空間一次元化することによって空間一次元のシミュレーション温度分布を算出させ、該シミュレーション温度分布と上記実験データとして記憶された上記温度分布とを比較させ、上記シミュレーション温度分布が上記実験データとして記憶された温度分布との差が許容範囲内である場合に上記シミュレーションデータが上記実験データを導出可能であると判定させるという構成を採用する。
第１４の発明は、上記第１１〜第１３いずれかの発明において、上記コンピュータを判定部として機能させる場合に、複数の上記初期物質が混合された予混合ガスを、常温における消炎径よりも小さな直径に設定された内部流路を備えると共に流れ方向に上記予混合ガスの着火温度以上まで連続的に昇温するように加熱された管に、火炎形成位置が安定しかつ上記予混合ガスの流速に依存しなくなる流量で供給した際に得られる温度分布を上記実験データとして記憶される上記空間一次元の温度分布として用いるという構成を採用する。
第１５の発明は、上記第１１〜第１４いずれかの発明において、上記コンピュータを、解析された素反応を基に解析を行う際の解析条件を記憶する解析条件記憶手段と、上記化学反応を記述する複数の素反応式を上記解析条件に基づいて選択する選択手段として機能させ、上記判定部において上記シミュレーションデータが上記実験データを導出可能であると判定された場合に上記コンピュータを上記選択手段として機能させるという構成を採用する。
本発明によれば、シミュレーションデータが空間零次元で時間一次元化された方程式に基づいて算出されるため、シミュレーションデータを算出する計算を空間零次元で時間一次元化して計算量を減少させることが可能となる。
また、本発明によれば、空間一次元の温度分布に一致するように、素反応データに含まれるパラメータの調整が行われる。このため、実験環境における不確定要因に影響されずに、実現象である燃焼反応を的確に予測可能とするパラメータの解を得ることができ、正確なフルケミストリ解析を行うことが可能となる。
このように、本発明によれば、計算負荷が小さくかつ正確なフルケミストリ解析が可能となる。
また、本発明によれば、空間一次元の温度分布に一致するように、素反応データに含まれるパラメータの調整が行われる。このため、実験環境における不確定要因に影響されずに、実現象である燃焼反応を的確に予測可能とするパラメータの解を得ることができ、正確なフルケミストリ解析を行うことが可能となる。
このように、本発明によれば、計算負荷が小さくかつ正確なフルケミストリ解析が可能となる。
以下、図面を参照して、本発明に係る反応メカニズム解析装置、反応メカニズム解析方法及び反応メカニズム解析プログラムの一実施形態について説明する。
（第１実施形態）
本実施形態の反応メカニズム解析装置及び反応メカニズム解析方法は、燃料と酸素との燃焼反応（化学反応）をフルケミストリ解析するためのものであり、後述するマイクロフローリアクタ（実験装置）を用いて得られた空間一次元の温度分布（実験データ）と、素反応式、該素反応式に関連付けられたパラメータ、及び空間零次元で時間一次元化された方程式を用いて得られたシミュレーションデータとを比較し、シミュレーションデータが実験データを導出可能となるまで上記パラメータを調整しながら収束計算を行うことによってパラメータの解を得るものである。
なお、ここで言う空間零次元で時間一次元化された方程式とは、位置情報がパラメータとして含まれておらず、空間の概念が除外された方程式を意味する。
本実施形態の反応メカニズム解析装置及び反応メカニズム解析方法は、燃料と酸素との燃焼反応（化学反応）をフルケミストリ解析するためのものであり、後述するマイクロフローリアクタ（実験装置）を用いて得られた空間一次元の温度分布（実験データ）と、素反応式、該素反応式に関連付けられたパラメータ、及び空間零次元で時間一次元化された方程式を用いて得られたシミュレーションデータとを比較し、シミュレーションデータが実験データを導出可能となるまで上記パラメータを調整しながら収束計算を行うことによってパラメータの解を得るものである。
なお、ここで言う空間零次元で時間一次元化された方程式とは、位置情報がパラメータとして含まれておらず、空間の概念が除外された方程式を意味する。
まず最初に、本実施形態の反応メカニズム解析装置及び反応メカニズム解析方法において用いられる実験データを得るためのマイクロフローリアクタについて説明する。
図１は、マイクロフローリアクタ１００の概略構成を示す模式図である。この図に示すように、マイクロフローリアクタ１００は、管１１０と、供給装置１２０と、ヒータ１３０と、温度測定装置１４０と、測定処理部１５０とを備えている。
図１は、マイクロフローリアクタ１００の概略構成を示す模式図である。この図に示すように、マイクロフローリアクタ１００は、管１１０と、供給装置１２０と、ヒータ１３０と、温度測定装置１４０と、測定処理部１５０とを備えている。
管１１０は円筒形状の直管である。そして、管１１０の内部流路１１１の直径は、常温において内部流路１１１に形成される火炎が伝播できずに消炎する限界値である消炎距離よりも小さく設定されている。つまり、管１１０の内部流路１１１は、常温における消炎径よりも小さな直径に設定されている。
供給装置１２０は、燃焼反応される初期物質である燃料と酸化剤（例えば、酸素や外気）とが混合された予混合ガスＧを管１１０に供給するものであり、管１１０の一端から管１１０の内部流路１１１に予混合ガスＧを流入させることによって管１１０に予混合ガスＧを供給する。
この供給装置１２０は、管１１０に供給する予混合ガスＧの流量を調節可能であり、管１１０に供給する予混合ガスＧの流量を連続的に減少させて、その流量が管１１０の内部流路１１１に形成された火炎が安定しかつ火炎の形成位置が予混合ガスＧの流速に影響されない条件を満たす流量となるように設定する。なお、管１１０の内部流路１１１に形成された火炎が安定しかつ火炎の形成位置が予混合ガスＧの流速に影響されない条件を満たす流量については、後に詳説する。
また、供給装置１２０は、測定処理部１５０と電気的に接続されており、測定処理部１５０からの指令に基づいて、予混合ガスＧの流量や温度、予混合ガスＧにおける燃料と酸化剤との割合（すなわち燃料と酸化剤との初期濃度）を調節可能とされている。
この供給装置１２０は、管１１０に供給する予混合ガスＧの流量を調節可能であり、管１１０に供給する予混合ガスＧの流量を連続的に減少させて、その流量が管１１０の内部流路１１１に形成された火炎が安定しかつ火炎の形成位置が予混合ガスＧの流速に影響されない条件を満たす流量となるように設定する。なお、管１１０の内部流路１１１に形成された火炎が安定しかつ火炎の形成位置が予混合ガスＧの流速に影響されない条件を満たす流量については、後に詳説する。
また、供給装置１２０は、測定処理部１５０と電気的に接続されており、測定処理部１５０からの指令に基づいて、予混合ガスＧの流量や温度、予混合ガスＧにおける燃料と酸化剤との割合（すなわち燃料と酸化剤との初期濃度）を調節可能とされている。
ヒータ１３０は、管１１０を加熱するものであり、内部流路１１１の温度が常温から予混合ガスＧの流れ方向に連続的に昇温して途中部位において予混合ガスＧの想定着火温度以上となるように管１１０を加熱する。
温度測定装置１４０は、管１１０の内部流路１１１におけるガス温度を計測するものである。この温度測定装置１４０は、測定処理部１５０と電気的に接続されており、測定結果を測定処理部１５０に入力する。
測定処理部１５０は、温度測定装置１４０から入力される測定結果を管１１０の位置に対応付けて記憶することによって温度分布を取得するものである。また、測定処理部１５０は、供給装置１２０から管１１０に供給される予混合ガスＧの流量や温度、予混合ガスＧに含まれる燃料や酸化剤の種類、及び燃料と酸化剤との初期濃度等の実験条件を記憶している。
次に、上述のように構成されたマイクロフローリアクタ１００を用いた実験方法（以下、本実験と称する）について説明する。
供給装置１２０から燃料と酸化剤とが所定の初期濃度で混合された予混合ガスＧが、管１１０の内部流路１１１に供給される。
管１１０は内部流路１１１が予混合ガスＧの着火温度以上に加熱されている。このためて、管１１０の内部流路１１１に供給された予混合ガスＧは、管１１０の内部流路１１１を一端部から他端部に向かうに連れて加熱され、着火温度以上に加熱された時点で着火する。
管１１０は内部流路１１１が予混合ガスＧの着火温度以上に加熱されている。このためて、管１１０の内部流路１１１に供給された予混合ガスＧは、管１１０の内部流路１１１を一端部から他端部に向かうに連れて加熱され、着火温度以上に加熱された時点で着火する。
予混合ガスＧが着火されることによって形成された火炎は、管１１０の内部流路１１１における予混合ガスＧの流速が速い場合には振動する。これは、予混合ガスＧの着火と消火が短時間で繰り返されることに起因する。
一方、管１１０の内部流路１１１における予混合ガスＧの流速が遅い場合には、連続的に燃焼状態が維持され、火炎が安定する。
そこで、本実験では、まず先に管１１０の内部流路１１１における流速が十分に速くなるような流量の予混合ガスＧを管１１０に供給し、その後予混合ガスＧの流量を徐々に減少させ、火炎が安定されるまで予混合ガスＧの流量を減少させる。
一方、管１１０の内部流路１１１における予混合ガスＧの流速が遅い場合には、連続的に燃焼状態が維持され、火炎が安定する。
そこで、本実験では、まず先に管１１０の内部流路１１１における流速が十分に速くなるような流量の予混合ガスＧを管１１０に供給し、その後予混合ガスＧの流量を徐々に減少させ、火炎が安定されるまで予混合ガスＧの流量を減少させる。
なお、内部流路１１１における混合ガス（予混合ガスのみからなるガス、あるいは当該予混合ガスと中間生成物とが混合したガスの両方の意味を含み本発明の混合体に相当）の温度は、火炎が形成される位置において急激に上昇する。このため、内部流路１１１における混合ガス温度が急上昇する位置を温度測定装置１４０から入力される測定結果に基づいて解析し、当該位置が変動しなくなるまで予混合ガスＧの流量を減少させることによって火炎を安定化させる。
また、管１１０を内部が目視可能な材料によって形成し、撮像装置等によって火炎の形成位置を確認しながら、火炎の形成位置が安定するまで予混合ガスＧの流量を減少させるようにしても良い。
また、管１１０を内部が目視可能な材料によって形成し、撮像装置等によって火炎の形成位置を確認しながら、火炎の形成位置が安定するまで予混合ガスＧの流量を減少させるようにしても良い。
なお、予混合ガスＧの流量が、火炎が安定されるまで減少されると、火炎の形成位置は予混合ガスＧの流速に影響されなくなる。
図２は、管１１０の内部流路１１１における予混合ガスＧの流速と、安定火炎位置、着火位置及び消火位置との関係を示したグラフである。なお、図２に示すグラフは、直径が２ｍｍの管１１０に対して、メタンと空気の量論比の予混合ガスＧを供給することによって得られたデータに基づくものである。
この図に示すように、予混合ガスＧの流速が速い場合（４０〜１００ｃｍ／ｓ）には、火炎は安定するが、火炎の形成位置が流速に影響されることが分かる。また、これよりも予混合ガスＧの流速が遅い場合（５〜４０ｃｍ／ｓ）には、火炎が振動することが分かる。そして、さらに予混合ガスＧの流速が遅い場合（０．２〜５ｃｍ／ｓ）には、火炎が安定すると共に火炎の形成位置が流速に影響されないことが分かる。このように、管１１０の内部流路１１１における予混合ガスＧの流速に応じて、火炎の状態が変化する。
そして、本実験では、予混合ガスＧの流速が、火炎が安定すると共に火炎の形成位置が流速に影響されない流速となるように、予混合ガスＧの流量が設定される。
図２は、管１１０の内部流路１１１における予混合ガスＧの流速と、安定火炎位置、着火位置及び消火位置との関係を示したグラフである。なお、図２に示すグラフは、直径が２ｍｍの管１１０に対して、メタンと空気の量論比の予混合ガスＧを供給することによって得られたデータに基づくものである。
この図に示すように、予混合ガスＧの流速が速い場合（４０〜１００ｃｍ／ｓ）には、火炎は安定するが、火炎の形成位置が流速に影響されることが分かる。また、これよりも予混合ガスＧの流速が遅い場合（５〜４０ｃｍ／ｓ）には、火炎が振動することが分かる。そして、さらに予混合ガスＧの流速が遅い場合（０．２〜５ｃｍ／ｓ）には、火炎が安定すると共に火炎の形成位置が流速に影響されないことが分かる。このように、管１１０の内部流路１１１における予混合ガスＧの流速に応じて、火炎の状態が変化する。
そして、本実験では、予混合ガスＧの流速が、火炎が安定すると共に火炎の形成位置が流速に影響されない流速となるように、予混合ガスＧの流量が設定される。
このようにして火炎を形成位置が予混合ガスＧの流速に依存しないように安定させた後に、測定処理部１５０は、温度分布を取得して出力する。
そして、本実験においては、内部流路１１１が常温における消炎径よりも小さな直径に設定された非常に細長い直管である管１１０の内部流路１１１において燃焼反応が起こっている。このため、本実験における実験空間は、ＡＳＴＭ法による三次元に広がりを有する実験空間と異なり、流れ方向のみに広がりを有するものである。したがって、本実験の実験結果として取得された温度分布は、空間一次元における温度変化であると考えることができる。
また、本実験においては、火炎の形成位置が予混合ガスＧの流速に依存せず、さらにヒータ１３０による管１１０の加熱を行うことによって内部流路１１１が外部から熱管理されるため内部流路１１１における熱授受の影響は明確になっている。このため、実験環境における不確定要因の影響を排除して温度分布を取得することができる。
このように、本実験においては、実験環境における不確定要因の影響が排除された空間一次元の温度分布を取得することができる。
また、本実験においては、火炎の形成位置が予混合ガスＧの流速に依存せず、さらにヒータ１３０による管１１０の加熱を行うことによって内部流路１１１が外部から熱管理されるため内部流路１１１における熱授受の影響は明確になっている。このため、実験環境における不確定要因の影響を排除して温度分布を取得することができる。
このように、本実験においては、実験環境における不確定要因の影響が排除された空間一次元の温度分布を取得することができる。
続いて、本実施形態の反応メカニズム解析装置及び反応メカニズム解析方法について説明する。
図３は、本実施形態の反応メカニズム解析装置のハードウェア構成を示すブロック図である。本実施形態の反応メカニズム解析装置Ａ１は、例えばワークステーションやスーパーコンピュータ等のコンピュータによって具現化されるものであり、図３に示すように、外部記憶装置１０と、内部記憶装置２０と、入力装置３０と、出力装置４０と、ＣＰＵ（Central Processing Unit）５０とを備えている。
なお、外部記憶装置１０と、内部記憶装置２０と、入力装置３０と、出力装置４０と、ＣＰＵ５０とは、システムバス６０によって互いに接続されている。
図３は、本実施形態の反応メカニズム解析装置のハードウェア構成を示すブロック図である。本実施形態の反応メカニズム解析装置Ａ１は、例えばワークステーションやスーパーコンピュータ等のコンピュータによって具現化されるものであり、図３に示すように、外部記憶装置１０と、内部記憶装置２０と、入力装置３０と、出力装置４０と、ＣＰＵ（Central Processing Unit）５０とを備えている。
なお、外部記憶装置１０と、内部記憶装置２０と、入力装置３０と、出力装置４０と、ＣＰＵ５０とは、システムバス６０によって互いに接続されている。
外部記憶装置１０は、ＣＰＵ５０の演算処理に用いられる各種データ及びプログラム等を記憶するものであり、本実施形態においては、実験データ１１、素反応データ１２、演算データ１３及び反応メカニズム解析プログラム１４を格納している。
そして、このような外部記憶装置１０としては、例えばハードディスク装置やリムーバブルメディア装置を用いることができる。なお、実験データ１１、素反応データ１２、演算データ１３及び反応メカニズム解析プログラム１４は、予め光ディスクやメモリ等のリムーバブルメディアに記憶され、当該リムーバブルディスクが外部記憶装置１０の一構成要件として反応メカニズム解析装置に接続されることによって外部記憶装置１０に記憶されても良いし、またネットワークを介してダウンロードされるあるいは入力装置３０を介して入力されることによって外部記憶装置１０の一構成要件であるハードディスクに記憶されても良い。
そして、このような外部記憶装置１０としては、例えばハードディスク装置やリムーバブルメディア装置を用いることができる。なお、実験データ１１、素反応データ１２、演算データ１３及び反応メカニズム解析プログラム１４は、予め光ディスクやメモリ等のリムーバブルメディアに記憶され、当該リムーバブルディスクが外部記憶装置１０の一構成要件として反応メカニズム解析装置に接続されることによって外部記憶装置１０に記憶されても良いし、またネットワークを介してダウンロードされるあるいは入力装置３０を介して入力されることによって外部記憶装置１０の一構成要件であるハードディスクに記憶されても良い。
実験データ１１は、上述したマイクロフローリアクタ１００によって得られた実験データ及び実験条件を有している。より詳細には、実験データ１１は、実験環境における不確定要因の影響が排除された空間一次元の温度分布、供給装置１２０から管１１０に供給される予混合ガスＧの流量や温度、予混合ガスＧに含まれる燃料や酸化剤の種類、及び燃料と酸化剤との初期濃度等を有している。
素反応データ１２は、本実施形態の反応メカニズム解析装置Ａ１にて反応メカニズムを解析する（すなわちフルケミストリ解析を行う）対象である燃料と酸化剤との燃焼反応を記述する複数の素反応式と、これらの素反応式に関連するパラメータ（例えば頻度係数や活性化エネルギ）とが集められたデータベースである。なお、以下の説明において、当該データベースを素反応機構と称する。
この素反応機構としては、例えば「ＧＲＩ−ｍｅｃｈ ３．０」や「Ｃ１−Ｃｈｅｍｉｓｔｒｙ」等を用いることができる。
また、上述のような予め得られた素反応機構が存在しない場合には、素反応機構は、分子動力学法等によって得ることができる。
この素反応機構としては、例えば「ＧＲＩ−ｍｅｃｈ ３．０」や「Ｃ１−Ｃｈｅｍｉｓｔｒｙ」等を用いることができる。
また、上述のような予め得られた素反応機構が存在しない場合には、素反応機構は、分子動力学法等によって得ることができる。
演算データ１３は、ＣＰＵ５０が素反応データ１２及び入力装置３０から入力されるパラメータに基づいて中間生成物の濃度変化や混合ガスの温度変化を算出するための計算式である下式（１），（２）を有している。
なお、ここで言う混合ガスとは、未反応の予混合ガスや燃焼反応の過程で生成される中間生成物が混合されたガスである。
なお、ここで言う混合ガスとは、未反応の予混合ガスや燃焼反応の過程で生成される中間生成物が混合されたガスである。
そして、式（１）は、通常の空間零次元で時間一次元定圧燃焼における質量保存の式である。また、式（２）は、エネルギ保存の式であり、左辺第一項が混合ガスの内部エネルギを示す項であり、左辺第二項が素反応により移動するエネルギを示す項であり、左辺第三項が混合ガスと外部との熱伝達により移動するエネルギ（すなわち混合ガスとマイクロフローリアクタ１００の壁面（管１１０の内面）との間の熱伝達によって移動するエネルギ）を示す項である。
これらの式（１）と式（２）との連立方程式によって、中間生成物の濃度変化や混合ガスの温度変化が算出される。ここで、式（１）と式（２）とは、独立変数が時間であり、従属変数が濃度と温度であり、さらに位置をパラメータとして含んでいないものであり、空間零次元で時間一次元の方程式である。そして、中間生成物の濃度変化や混合ガスの温度変化は、上述のように空間の概念が除外されて空間零次元で時間一次元化された方程式によって算出される。
これらの式（１）と式（２）との連立方程式によって、中間生成物の濃度変化や混合ガスの温度変化が算出される。ここで、式（１）と式（２）とは、独立変数が時間であり、従属変数が濃度と温度であり、さらに位置をパラメータとして含んでいないものであり、空間零次元で時間一次元の方程式である。そして、中間生成物の濃度変化や混合ガスの温度変化は、上述のように空間の概念が除外されて空間零次元で時間一次元化された方程式によって算出される。
また、演算データ１３は、上述の式（１），（２）を用いて算出された混合ガスの温度変化を空間一次元化するための下式（３），（４）を有している。
なお、マイクロフローリアクタ１００における管１１０の熱容量は小流量の予混合ガスによる発熱量より十分に大きい。このため、火炎による管壁温度分布の変化は無視することができる。したがって、管壁温度分布は管１１０の位置に対する関数として与えられる。よって、下式（３），（４）によって上述の式（１），（２）を用いて算出された混合ガスの温度変化を空間一次元化することが可能となる。
なお、マイクロフローリアクタ１００における管１１０の熱容量は小流量の予混合ガスによる発熱量より十分に大きい。このため、火炎による管壁温度分布の変化は無視することができる。したがって、管壁温度分布は管１１０の位置に対する関数として与えられる。よって、下式（３），（４）によって上述の式（１），（２）を用いて算出された混合ガスの温度変化を空間一次元化することが可能となる。
反応メカニズム解析プログラム１４は、本実施形態の反応メカニズム解析装置Ａ１を、後に説明する図３に示す各機能構成として機能させるためのプログラムである。
なお、反応メカニズム解析プログラム１４については、後の図４を参照する反応メカニズム解析装置Ａ１の機能構成の説明において合わせて詳説する。
なお、反応メカニズム解析プログラム１４については、後の図４を参照する反応メカニズム解析装置Ａ１の機能構成の説明において合わせて詳説する。
内部記憶装置２０は、ＣＰＵ５０の動作プログラムを記憶すると共に、ＣＰＵ５０の制御の下、外部記憶装置２０に記憶された各種データ及びプログラムを一時的に記憶し、さらにＣＰＵ５０の演算結果を一時的に記憶するものである。そして、この内部記憶装置２０としては、ＲＡＭ（Random Access Memory）やＲＯＭ（Read Only Memory）等が用いられる。
入力装置３０は、本実施形態の反応メカニズム解析装置Ａ１に対して外部からデータを入力するためのものであり、例えばキーボードやマウスが用いられる。ただし、本実施形態の反応メカニズム解析装置Ａ１に対して通信路を介してデータの入力が行われる場合には、通信路とのインターフェイスとなる通信装置を入力装置３０として用いる場合もある。
出力装置４０は、ＣＰＵ５０の制御の下、外部記憶装置１０あるいは内部記憶装置２０に記憶されたデータを出力するものであり、例えばデータを可視化するためのディスプレイやプリンタが用いられる。ただし、本実施形態の反応メカニズム解析装置Ａ１から通信路を介してデータの出力が行われる場合には、通信路とのインターフェイスとなる通信装置を出力装置４０として用いる場合もある。
ＣＰＵ５０は、本実施形態の反応メカニズム解析装置Ａ１の動作全体を制御するものである。そして、本実施形態の反応メカニズム解析装置Ａ１においてＣＰＵ５０は、反応メカニズム解析プログラム１４に基づいて、中間生成物の濃度変化や混合ガスの温度変化の算出、シミュレーションデータと実験データとの比較及び判定、素反応機構が含むパラメータの調整等を行う。
図４は、本実施形態の反応メカニズム解析装置Ａ１の機能構成を示すブロック図である。この図に示すように、本実施形態の反応メカニズム解析装置Ａ１は、上述の反応メカニズム解析プログラム１４によって、実験データ記憶部１（実験データ記憶手段）と、素反応データ記憶部２（素反応データ記憶手段）と、演算データ記憶部３（演算データ記憶手段）と、算出部４（算出手段）と、判定部５（判定手段）と、パラメータ調整部６（パラメータ調整手段）と、制御部７（制御手段）として機能する。
実験データ記憶部１は上記実験データ１１を記憶する。また、素反応データ記憶部２は上記素反応データ１２を記憶する。また、演算データ記憶部３は上記演算データ１３を記憶する。
そして、反応メカニズム解析プログラム１４は、実験データ１１、素反応データ１２及び演算データ１３を記憶する外部記憶装置１０、あるいは、実験データ１１、素反応データ１２及び演算データ１３を移した内部記憶装置２０を実験データ記憶部１、素反応データ記憶部２及び演算データ記憶部３として機能させる。
そして、反応メカニズム解析プログラム１４は、実験データ１１、素反応データ１２及び演算データ１３を記憶する外部記憶装置１０、あるいは、実験データ１１、素反応データ１２及び演算データ１３を移した内部記憶装置２０を実験データ記憶部１、素反応データ記憶部２及び演算データ記憶部３として機能させる。
算出部４は、演算データ記憶部３に演算データ１３として記憶された空間零次元で時間一次元化された方程式である上式（１），（２）、実験データ記憶部１に記憶された実験データ１１に含まれる実験条件、及び素反応データ記憶部２に記憶された素反応データ１２を用いて中間生成物の濃度変化及び混合ガスの温度変化を算出する。
そして、反応メカニズム解析プログラム１４は、外部記憶装置１０あるいは内部記憶装置２０に記憶された実験データ１１、素反応データ１２及び演算データ１３に基づいてＣＰＵ５０に中間生成物の濃度変化及び混合ガスの温度変化を算出させることによって、反応メカニズム解析装置Ａ１を算出部４として機能させる。
そして、反応メカニズム解析プログラム１４は、外部記憶装置１０あるいは内部記憶装置２０に記憶された実験データ１１、素反応データ１２及び演算データ１３に基づいてＣＰＵ５０に中間生成物の濃度変化及び混合ガスの温度変化を算出させることによって、反応メカニズム解析装置Ａ１を算出部４として機能させる。
判定部５は、算出部４によって算出された混合ガスの温度変化（シミュレーションデータ）と実験データ１１とを比較し、算出部４によって算出された混合ガスの温度変化が実験データ１１に含まれる温度分布（内部流路１１１におけるガス温度分布）を導出可能なものであるかの判定を行う。
具体的には、判定部５は、算出部４によって算出された混合ガスの温度変化を演算データ記憶部３に演算データ１３として記憶された上式（３），（４）を用いて空間一次元化して、シミュレーション温度分布を算出する。そして、判定部５は、シミュレーション温度分布と、実験データ記憶部１に実験データ１１に含まれる温度分布とを比較する。判定部５は、比較の結果、シミュレーション温度分布が実験データ１１に含まれる温度分布と一致若しくは許容される変動範囲に収まっている場合には、算出部４によって算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出可能なものであると判定し、シミュレーション温度分布が実験データ１１に含まれる温度分布から許容される変動範囲に収まっていない場合には、算出部４によって算出された混合ガスの温度変化（すなわちシミュレーションデータ）が実験データに含まれる温度分布を導出不能なものであると判定する。
そして、反応メカニズム解析プログラム１４は、外部記憶装置１０あるいは内部記憶装置２０に記憶された演算データ１３及び算出された混合ガスの温度変化に基づいてＣＰＵ５０にシミュレーション温度分布を算出させ、さらにＣＰＵ５０に外部記憶装置１０あるいは内部記憶装置２０に記憶された実験データ１１に含まれる温度分布とシミュレーション温度分布とを比較させ、算出された混合ガスの温度変化が実験データに含まれる温度分布を導出可能なものであるかを判定させることによって、反応メカニズム解析装置Ａ１を判定部５として機能させる。
具体的には、判定部５は、算出部４によって算出された混合ガスの温度変化を演算データ記憶部３に演算データ１３として記憶された上式（３），（４）を用いて空間一次元化して、シミュレーション温度分布を算出する。そして、判定部５は、シミュレーション温度分布と、実験データ記憶部１に実験データ１１に含まれる温度分布とを比較する。判定部５は、比較の結果、シミュレーション温度分布が実験データ１１に含まれる温度分布と一致若しくは許容される変動範囲に収まっている場合には、算出部４によって算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出可能なものであると判定し、シミュレーション温度分布が実験データ１１に含まれる温度分布から許容される変動範囲に収まっていない場合には、算出部４によって算出された混合ガスの温度変化（すなわちシミュレーションデータ）が実験データに含まれる温度分布を導出不能なものであると判定する。
そして、反応メカニズム解析プログラム１４は、外部記憶装置１０あるいは内部記憶装置２０に記憶された演算データ１３及び算出された混合ガスの温度変化に基づいてＣＰＵ５０にシミュレーション温度分布を算出させ、さらにＣＰＵ５０に外部記憶装置１０あるいは内部記憶装置２０に記憶された実験データ１１に含まれる温度分布とシミュレーション温度分布とを比較させ、算出された混合ガスの温度変化が実験データに含まれる温度分布を導出可能なものであるかを判定させることによって、反応メカニズム解析装置Ａ１を判定部５として機能させる。
パラメータ調整部６は、素反応データ記憶部２に記憶された素反応データ１２に含まれるパラメータ（頻度係数や活性化エネルギ）を調整する。
具体的には、パラメータ調整部６は、予め定められた順序にて素反応データに含まれる素反応式に関連付けられたパラメータを所定の値分だけ増加あるいは減少させることによって、素反応データ１２に含まれるパラメータを調整する。
そして、反応メカニズム解析プログラム１４は、外部記憶装置１０あるいは内部記憶装置２０に記憶された素反応データ１２に含まれるパラメータをＣＰＵ５０に調整させることによって、反応メカニズム解析装置Ａ１をパラメータ調整部６として機能させる。
具体的には、パラメータ調整部６は、予め定められた順序にて素反応データに含まれる素反応式に関連付けられたパラメータを所定の値分だけ増加あるいは減少させることによって、素反応データ１２に含まれるパラメータを調整する。
そして、反応メカニズム解析プログラム１４は、外部記憶装置１０あるいは内部記憶装置２０に記憶された素反応データ１２に含まれるパラメータをＣＰＵ５０に調整させることによって、反応メカニズム解析装置Ａ１をパラメータ調整部６として機能させる。
制御部７は、実験データ記憶部１、素反応データ記憶部２、演算データ記憶部３、算出部４、判定部５及びパラメータ調整部６の動作手順等の制御を行う。
そして、本実施形態の反応メカニズム解析装置Ａ１において制御部７は、判定部５が算出部４によって算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出可能であると判定した場合には、フルケミストリ解析が完了したと判断し、その時点で外部記憶装置１０あるいは内部記憶装置２０に記憶された最新の中間生成物の濃度変化、混合ガスの温度変化及び素反応データ１２に含まれるパラメータを解析結果として出力装置４０に入力する。
一方、制御部７は、判定部５において算出部４によって算出された混合ガスの温度変化が実験データに含まれる温度分布を導出不能であると判定した場合には、パラメータ調整部６に素反応データ１２に含まれるパラメータを調整させ、再度算出部４に中間生成物の濃度変化及び混合ガスの温度変化を算出させ、さらに判定部５に算出された混合ガスの温度変化が実験データ１１に含まれる管壁分布温度を導出可能なものであるかの判定をさせる。
そして、反応メカニズム解析プログラム１４は、ＣＰＵ５０を制御部７として機能させる。
そして、本実施形態の反応メカニズム解析装置Ａ１において制御部７は、判定部５が算出部４によって算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出可能であると判定した場合には、フルケミストリ解析が完了したと判断し、その時点で外部記憶装置１０あるいは内部記憶装置２０に記憶された最新の中間生成物の濃度変化、混合ガスの温度変化及び素反応データ１２に含まれるパラメータを解析結果として出力装置４０に入力する。
一方、制御部７は、判定部５において算出部４によって算出された混合ガスの温度変化が実験データに含まれる温度分布を導出不能であると判定した場合には、パラメータ調整部６に素反応データ１２に含まれるパラメータを調整させ、再度算出部４に中間生成物の濃度変化及び混合ガスの温度変化を算出させ、さらに判定部５に算出された混合ガスの温度変化が実験データ１１に含まれる管壁分布温度を導出可能なものであるかの判定をさせる。
そして、反応メカニズム解析プログラム１４は、ＣＰＵ５０を制御部７として機能させる。
続いて、上述した反応メカニズム解析装置Ａ１を用いた反応メカニズム解析方法について、図５のフローチャートを参照して説明する。
なお、本実施形態における反応メカニズム解析方法においては、実験データ１１、素反応データ１２及び演算データ１３が既に実験データ記憶部１、素反応データ記憶部２、演算データ記憶部３（外部記憶装置１０あるいは内部記憶装置２０）に記憶されているものとする。
なお、本実施形態における反応メカニズム解析方法においては、実験データ１１、素反応データ１２及び演算データ１３が既に実験データ記憶部１、素反応データ記憶部２、演算データ記憶部３（外部記憶装置１０あるいは内部記憶装置２０）に記憶されているものとする。
まず最初に制御部７の制御の下、算出部４によって中間生成物の濃度変化及び混合ガスの温度変化が算出される（算出工程Ｓ１）。
この際、算出部４は、演算データ記憶部３に演算データ１３として記憶された空間零次元で時間一次元化された方程式である上式（１），（２）、実験データ記憶部１に記憶された実験データ１１に含まれる実験条件、素反応データ記憶部２に記憶された素反応データ１２を用いて中間生成物の濃度変化及び混合ガスの温度変化を算出する。
なお、具体的には、ＣＰＵ５０が、外部記憶装置１０あるいは内部記憶装置２０に記憶された実験データ１１に含まれる実験条件と、素反応データ１２に含まれるパラメータと、外部記憶装置１０あるいは内部記憶装置２０に記憶された演算データ１３に含まれる式（１），（２）とを用い、さらに式（１）と式（２）とを連立させて濃度変化と温度変化について解くことによって中間生成物の濃度変化及び混合ガスの温度変化が算出される。
この際、算出部４は、演算データ記憶部３に演算データ１３として記憶された空間零次元で時間一次元化された方程式である上式（１），（２）、実験データ記憶部１に記憶された実験データ１１に含まれる実験条件、素反応データ記憶部２に記憶された素反応データ１２を用いて中間生成物の濃度変化及び混合ガスの温度変化を算出する。
なお、具体的には、ＣＰＵ５０が、外部記憶装置１０あるいは内部記憶装置２０に記憶された実験データ１１に含まれる実験条件と、素反応データ１２に含まれるパラメータと、外部記憶装置１０あるいは内部記憶装置２０に記憶された演算データ１３に含まれる式（１），（２）とを用い、さらに式（１）と式（２）とを連立させて濃度変化と温度変化について解くことによって中間生成物の濃度変化及び混合ガスの温度変化が算出される。
続いて、制御部７の制御の下、判定部５によって、算出部４によって算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出可能なものであるかの判定を行う（判定工程Ｓ２）。
この際、判定部５は、算出部４によって算出された混合ガスの温度変化を演算データ記憶部３に演算データ１３として記憶された上式（３），（４）を用いて空間一次元化して、シミュレーション温度分布を算出する。そして、判定部５は、シミュレーション温度分布と、実験データ記憶部１に記憶された実験データ１１に含まれる温度分布とを比較する。判定部５は、比較の結果、シミュレーション温度分布が実験データ１１に含まれる温度分布と一致若しくは許容される変動範囲に収まっている場合には、算出部４によって算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出可能なものであると判定し、シミュレーション温度分布が実験データ１１に含まれる温度分布から許容される変動範囲に収まっていない場合には、算出部４によって算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出不能であると判定する。
具体的には、ＣＰＵ５０が、算出された混合ガスの温度変化を外部記憶装置１０あるいは内部記憶装置２０に記憶された演算データ１３に含まれる式（３），（４）を用いて空間一次元化することでシミュレーション温度分布を算出する。そして、ＣＰＵ５０が、シミュレーション温度分布と外部記憶装置１０あるいは内部記憶装置２０に記憶された実験データ１１に含まれる温度分布とを比較する。続いてＣＰＵ５０が、算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出可能なものであるかを判定する。
この際、判定部５は、算出部４によって算出された混合ガスの温度変化を演算データ記憶部３に演算データ１３として記憶された上式（３），（４）を用いて空間一次元化して、シミュレーション温度分布を算出する。そして、判定部５は、シミュレーション温度分布と、実験データ記憶部１に記憶された実験データ１１に含まれる温度分布とを比較する。判定部５は、比較の結果、シミュレーション温度分布が実験データ１１に含まれる温度分布と一致若しくは許容される変動範囲に収まっている場合には、算出部４によって算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出可能なものであると判定し、シミュレーション温度分布が実験データ１１に含まれる温度分布から許容される変動範囲に収まっていない場合には、算出部４によって算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出不能であると判定する。
具体的には、ＣＰＵ５０が、算出された混合ガスの温度変化を外部記憶装置１０あるいは内部記憶装置２０に記憶された演算データ１３に含まれる式（３），（４）を用いて空間一次元化することでシミュレーション温度分布を算出する。そして、ＣＰＵ５０が、シミュレーション温度分布と外部記憶装置１０あるいは内部記憶装置２０に記憶された実験データ１１に含まれる温度分布とを比較する。続いてＣＰＵ５０が、算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出可能なものであるかを判定する。
そして、判定工程Ｓ２において、算出工程Ｓ１にて算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出可能であると判定された場合には、制御部７（ＣＰＵ５０）は、フルケミストリ解析が完了したと判断し、その時点で外部記憶装置１０あるいは内部記憶装置２０に記憶された最新の中間生成物の濃度変化、混合ガスの温度変化及び素反応データ１２に含まれるパラメータを解析結果として出力装置４０に入力する（出力工程Ｓ３）。
一方、判定工程Ｓ２において、算出工程Ｓ１にて算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出不能であると判定された場合には、制御部７の制御の下、パラメータ調整部６によって、素反応データ記憶部２に記憶された素反応データ１２に含まれるパラメータを調整する（パラメータ調整工程Ｓ４）。
具体的には、ＣＰＵ５０が外部記憶装置１０あるいは内部記憶装置２０に記憶された素反応データ１２に含まれるパラメータを調整する。
そして、パラメータ調整工程Ｓ４が完了すると、再度算出工程Ｓ１に移行する。
具体的には、ＣＰＵ５０が外部記憶装置１０あるいは内部記憶装置２０に記憶された素反応データ１２に含まれるパラメータを調整する。
そして、パラメータ調整工程Ｓ４が完了すると、再度算出工程Ｓ１に移行する。
このような本実施形態の反応メカニズム解析方法によれば、判定工程Ｓ２において、算出工程Ｓ１にて算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出可能となるまで繰り返し素反応データ１２に含まれるパラメータの調整が行われる。そして、算出工程Ｓ１にて算出された混合ガスの温度変化が実験データ１１に含まれる温度分布を導出可能となった場合、すなわちシミュレーション温度分布が実験データ１１が含む温度分布に対して一致する若しくは許容範囲内となる上記パラメータが収束計算にて算出された場合には、当該パラメータ及び該パラメータに基づいた中間生成物の濃度変化及び混合ガスの温度変化が出力され、フルケミストリ解析が終了する。
以上のような本実施形態の反応メカニズム解析装置Ａ１、反応メカニズム解析方法及び反応メカニズム解析プログラムによれば、シミュレーションデータが空間零次元で時間一次元化された方程式に基づいて算出されるため、計算量を減少させることが可能となる。
また、本実施形態の反応メカニズム解析装置Ａ１、反応メカニズム解析方法及び反応メカニズム解析プログラムによれば、マイクロフローリアクタによって得られた実験環境における不確定要因に影響されない空間一次元の温度分布に合うように、素反応データ１２に含まれるパラメータの調整が行われる。このため、実験環境における不確定要因に影響されずに、実現象である燃焼反応を的確に予測可能とするパラメータの解を得ることができ、正確なフルケミストリ解析を行うことが可能となる。
このように、本実施形態の反応メカニズム解析装置Ａ１、反応メカニズム解析方法及び反応メカニズム解析プログラムによれば、計算負荷が小さくかつ正確なフルケミストリ解析が可能となる。
また、本実施形態の反応メカニズム解析装置Ａ１、反応メカニズム解析方法及び反応メカニズム解析プログラムによれば、マイクロフローリアクタによって得られた実験環境における不確定要因に影響されない空間一次元の温度分布に合うように、素反応データ１２に含まれるパラメータの調整が行われる。このため、実験環境における不確定要因に影響されずに、実現象である燃焼反応を的確に予測可能とするパラメータの解を得ることができ、正確なフルケミストリ解析を行うことが可能となる。
このように、本実施形態の反応メカニズム解析装置Ａ１、反応メカニズム解析方法及び反応メカニズム解析プログラムによれば、計算負荷が小さくかつ正確なフルケミストリ解析が可能となる。
（第２実施形態）
次に、本発明の第２実施形態について説明する。なお、本第２実施形態の説明において、上記第１実施形態と同様の部分については、その説明を省略あるいは簡略化する。
次に、本発明の第２実施形態について説明する。なお、本第２実施形態の説明において、上記第１実施形態と同様の部分については、その説明を省略あるいは簡略化する。
本実施形態の反応メカニズム解析装置Ａ２は、図６に示すように、外部記憶装置１０に、解析された素反応を基に行われる解析（例えばＣＦＤ解析）の解析条件１５が格納されている。なお、ここで言う解析条件とは、ＣＦＤ解析等において許容される計算量や、ＣＦＤ解析等によって解析したい中間生成物の種類等の情報を含むものである。
そして、本実施形態の反応メカニズム解析装置Ａ２は、図７に示すように、反応メカニズム解析プログラム１４によって、解析条件記憶部８（解析条件記憶手段）と、選択部９（選択手段）として機能される。
解析条件記憶部８は上記解析条件１５を記憶する。
そして、反応メカニズム解析プログラム１４は、解析条件１５を記憶する外部記憶装置１０、あるいは、解析条件１５を移した内部記憶装置２０を解析条件記憶部８として機能させる。
そして、反応メカニズム解析プログラム１４は、解析条件１５を記憶する外部記憶装置１０、あるいは、解析条件１５を移した内部記憶装置２０を解析条件記憶部８として機能させる。
選択部９は、化学反応を記述する複数の素反応式を解析条件１５に基づいて選択すると共に選択された素反応式に関連付けられたパラメータを調整するものである。
具体的には、選択部９は、解析条件記憶部８に記憶された解析条件１５に基づいて選択すべき素反応式（ＣＦＤ解析等により解析したい中間生成物を含む素反応式）を優先して選択し、さらにＣＦＤ解析等を行う場合に許容される計算負荷を満足するように選択する素反応式の数を決定し、これらに基づいて素反応式の選択を行う。
また、選択部９は、選択した素反応式のみを用いて再度シミュレーション温度分布を算出し、実験データ１１が含む温度分布と比較し、シミュレーション温度分布が実験データ１１に含まれる温度分布に対して許容範囲に収まらない場合には、シミュレーション温度分布が実験データ１１に含まれる温度分布に対して一致若しくは許容範囲に収まるまで繰り返し選択した素反応式に関連するパラメータを調整して上記計算を繰り返す。
そして、反応メカニズム解析プログラム１４は、外部記憶装置１０あるいは内部記憶装置２０及びＣＰＵ５０を用いて本実施形態のメカニズム解析装置Ａ１を選択部９として機能させる。
具体的には、選択部９は、解析条件記憶部８に記憶された解析条件１５に基づいて選択すべき素反応式（ＣＦＤ解析等により解析したい中間生成物を含む素反応式）を優先して選択し、さらにＣＦＤ解析等を行う場合に許容される計算負荷を満足するように選択する素反応式の数を決定し、これらに基づいて素反応式の選択を行う。
また、選択部９は、選択した素反応式のみを用いて再度シミュレーション温度分布を算出し、実験データ１１が含む温度分布と比較し、シミュレーション温度分布が実験データ１１に含まれる温度分布に対して許容範囲に収まらない場合には、シミュレーション温度分布が実験データ１１に含まれる温度分布に対して一致若しくは許容範囲に収まるまで繰り返し選択した素反応式に関連するパラメータを調整して上記計算を繰り返す。
そして、反応メカニズム解析プログラム１４は、外部記憶装置１０あるいは内部記憶装置２０及びＣＰＵ５０を用いて本実施形態のメカニズム解析装置Ａ１を選択部９として機能させる。
そして、本実施形態の反応メカニズム解析装置Ａ１において制御部７は、判定部５においてシミュレーション温度分布が実験データ１１に含まれる温度分布を導出可能であると判定された場合に制御部９による素反応式の選択及びパラメータの調整を実行させる。
つまり、本実施形態の反応メカニズム解析方法においては、図８に示すように、判定工程Ｓ２と出力工程Ｓ３との間に、素反応式の選択及び選択された素反応式に関連付けられたパラメータの調整を行う選択工程Ｓ５が実行される。
つまり、本実施形態の反応メカニズム解析方法においては、図８に示すように、判定工程Ｓ２と出力工程Ｓ３との間に、素反応式の選択及び選択された素反応式に関連付けられたパラメータの調整を行う選択工程Ｓ５が実行される。
このような本実施形態の反応メカニズム解析装置Ａ２、反応メカニズム解析方法及び反応メカニズム解析プログラムによれば、正確なフルケミストリ解析によって得られた素反応式を、解析された素反応を基に行われる解析の解析条件に最適化することが可能となる。
以上、添付図面を参照しながら本発明に係る反応メカニズム解析装置、反応メカニズム解析方法及び反応メカニズム解析プログラムの好適な実施形態について説明したが、本発明は、上記実施形態に限定されないことは言うまでもない。上述した実施形態において示した各構成部材の諸形状や組み合わせ等は一例であって、本発明の主旨から逸脱しない範囲において設計要求等に基づき種々変更可能である。
例えば、上記実施形態においては、燃料と酸化剤との燃焼反応を構成する素反応の解析を行う場合について説明した。
しかしながら、本発明はこれに限定されるものではなく、他の物質同士の化学反応を構成する素反応の解析に用いることもできる。
しかしながら、本発明はこれに限定されるものではなく、他の物質同士の化学反応を構成する素反応の解析に用いることもできる。
また、上記実施形態においては、実験データ１１として含まれる温度分布を、マイクロフローリアクタ１００によって取得する構成について説明した。
しかしながら、本発明はこれに限定されるものではなく、燃焼反応の際の正確な空間一次元の温度分布を取得可能な装置であれば、空間一次元の温度分布を他の実験装置によって取得しても良い。
しかしながら、本発明はこれに限定されるものではなく、燃焼反応の際の正確な空間一次元の温度分布を取得可能な装置であれば、空間一次元の温度分布を他の実験装置によって取得しても良い。
また、上記実施形態においては、内部流路１１１を流れるガスの温度分布を空間一次元の温度分布として用いる構成について説明した。
しかしながら、本発明はこれに限定されるものではなく、例えば、火炎の形成位置が予混合ガスＧの流速に影響されない条件を満たす場合には、内部流路１１１を流れるガスと管１１０の管壁温度とがほぼ等しくなるため、管１１０の壁面温度分布（管壁温度分布）を空間一次元の温度分布として用いることも可能である。
また、予混合ガスＧの流速が速い場合には火炎の形成位置が振動するが、この場合の実験結果を実験データ１１に含め、当該実験結果を素反応データ１２に含まれるパラメータの調整に用いることもできる。
しかしながら、本発明はこれに限定されるものではなく、例えば、火炎の形成位置が予混合ガスＧの流速に影響されない条件を満たす場合には、内部流路１１１を流れるガスと管１１０の管壁温度とがほぼ等しくなるため、管１１０の壁面温度分布（管壁温度分布）を空間一次元の温度分布として用いることも可能である。
また、予混合ガスＧの流速が速い場合には火炎の形成位置が振動するが、この場合の実験結果を実験データ１１に含め、当該実験結果を素反応データ１２に含まれるパラメータの調整に用いることもできる。
また、上記実施形態においては、空間零次元で時間一次元化された方程式を用いて算出された混合ガスの温度変化を一次元化することによってシミュレーション温度分布を算出し、シミュレーション温度分布と実験により得られた温度分布とを比較して判定する構成について説明した。
しかしながら、本発明はこれに限定されるものではなく、例えば実験により得られた温度分布を空間零次元で時間一次元化して、空間零次元で時間一次元化された方程式を用いて算出された混合ガスの温度変化と比較して判定を行う構成を採用することもできる。
しかしながら、本発明はこれに限定されるものではなく、例えば実験により得られた温度分布を空間零次元で時間一次元化して、空間零次元で時間一次元化された方程式を用いて算出された混合ガスの温度変化と比較して判定を行う構成を採用することもできる。
また、上記実施形態においては、最終的に解析が完了した後に、中間生成物の濃度変化、混合ガスの温度変化及び素反応データに含まれるパラメータを解析結果として出力する構成について説明した。
しかしながら、本発明はこれに限定されるものではなく、最終的に解析が完了する前の各種データ（例えば、判定工程においてＮＧとされたパラメータ及び該パラメータに基づく中間生成物の濃度変化、混合ガスの温度変化等）を出力するようにしても良い。
しかしながら、本発明はこれに限定されるものではなく、最終的に解析が完了する前の各種データ（例えば、判定工程においてＮＧとされたパラメータ及び該パラメータに基づく中間生成物の濃度変化、混合ガスの温度変化等）を出力するようにしても良い。
また、上記実施形態においては、管１１０の内部流路１１１における混合ガスの温度変化から混合ガスのシミュレーション温度分布を算出し、実験データとして管１１０の内部流路１１１における混合ガスの温度分布を取得し、上記シミュレーション温度分布と実験データとされる温度分布とを比較し、これによってシミュレーションデータが実験データを導出可能であるかを判定する構成について説明した。
しかしながら、本発明はこれに限定されるものではなく、管１１０の内部流路１１１における中間生成物の濃度変化から濃度分布をシミュレーションデータとして算出し、実験データとして管１１０の内部流路１１１における中間生成物の濃度分布を取得し、シミュレーションデータとして算出された濃度分布と実験データとされる濃度分布を比較し、これによってシミュレーションデータが実験データを導出可能であるかを判定しても良い。
なお、中間生成物の濃度は、管１１０の内部流路１１１における混合ガスをサンプリングする、あるいはレーザ計測による測定を行うことによって実験データとして取得することが可能である。また、これらの濃度計測は管１１０のいずれの位置で実施しても良い。このため、例えば、冷炎と熱炎が発生するような場合、冷炎と熱炎との中間位置で濃度計測を行うことで冷炎からの中間生成物濃度のみを計測することも可能である。
しかしながら、本発明はこれに限定されるものではなく、管１１０の内部流路１１１における中間生成物の濃度変化から濃度分布をシミュレーションデータとして算出し、実験データとして管１１０の内部流路１１１における中間生成物の濃度分布を取得し、シミュレーションデータとして算出された濃度分布と実験データとされる濃度分布を比較し、これによってシミュレーションデータが実験データを導出可能であるかを判定しても良い。
なお、中間生成物の濃度は、管１１０の内部流路１１１における混合ガスをサンプリングする、あるいはレーザ計測による測定を行うことによって実験データとして取得することが可能である。また、これらの濃度計測は管１１０のいずれの位置で実施しても良い。このため、例えば、冷炎と熱炎が発生するような場合、冷炎と熱炎との中間位置で濃度計測を行うことで冷炎からの中間生成物濃度のみを計測することも可能である。
Ａ１，Ａ１……反応メカニズム解析装置、１……実験データ記憶部（実験データ記憶手段）、２……素反応データ記憶部（素反応データ記憶手段）、３……演算データ記憶部（演算データ記憶手段）、４……算出部（算出手段）、５……判定部（判定手段）、６……パラメータ調整部（パラメータ調整手段）、７……制御部（制御手段）、８……解析条件記憶部（解析条件記憶手段）、９……選択部（選択手段）
Claims (8)
- 複数種類の初期物質同士の化学反応を構成する素反応を解析する反応メカニズム解析装置であって、
流通方向に対してのみ物質移動及び熱移動が起こりうると近似された内部流路を有する管の一端側から複数の前記初期物質が混合された予混合ガスを導入して化学反応させる実験より得られた前記流通方向における位置に対する温度分布を実験データとして記憶する実験データ記憶手段と、
前記化学反応を記述する複数の素反応式及び該素反応式に関連付けられたパラメータを素反応データとして記憶する素反応データ記憶手段と、
前記素反応データからシミュレーションデータを算出するための方程式であって、空間座標に基づかず、時間に対する変数を含むと共に温度の次元を有する方程式を演算データとして記憶する演算データ記憶手段と、
前記素反応データ及び前記演算データとして記憶された方程式を用いて、前記素反応によって生成される中間生成物と前記初期物質との混合体の温度変化を前記シミュレーションデータとして算出する算出部と、
前記混合体の移動速度と、前記混合体の温度と、前記予混合ガスの供給速度と、前記予混合ガスの初期温度と、前記流通方向の位置と、時間との関係式に基づいて、前記シミュレーションデータから、前記管に前記予混合ガスを導入して化学反応させたときの前記流通方向における位置に対する温度分布を算出してシミュレーションデータに基づく温度分布とし、前記シミュレーションデータに基づく温度分布と前記実験データに基づく温度分布との差異が許容範囲にあることを導出可能であると定義し、前記シミュレーションデータが前記実験データを導出可能であるかを判定する判定部と、
前記素反応データに含まれる前記パラメータを調整可能なパラメータ調整手段と、
前記判定部において前記シミュレーションデータが前記実験データを導出可能であると判定されるまで、前記パラメータ調整手段による前記パラメータの調整及び調整後の前記パラメータを用いた前記算出部による前記シミュレーションデータの算出を繰り返し実行させる制御手段と
を備えることを特徴とする反応メカニズム解析装置。 - 前記演算データとして記憶された方程式は、
中間生成物の濃度変化を示す項を含む質量保存の式と、
素反応によって生成される中間生成物と前記初期物質との混合体の内部エネルギを示す項、前記素反応により移動するエネルギを示す項、及び前記混合体と外部との熱伝達により移動するエネルギを示す項を含むエネルギ保存の式と
を有することを特徴とする請求項１記載の反応メカニズム解析装置。 - 前記実験データとして記憶される前記流通方向における位置に対する温度分布は、
複数の前記初期物質として燃料と酸化剤とが混合された予混合ガスを、前記予混合ガスの着火温度以上まで連続的に昇温するように加熱された管に、火炎形成位置が安定する流量で供給した際に得られる温度分布である
ことを特徴とする請求項１または２記載の反応メカニズム解析装置。 - 解析された素反応を基に解析を行う際の解析条件を記憶する解析条件記憶手段と、前記化学反応を記述する複数の素反応式を前記解析条件に基づいて選択すると共に選択された素反応式に関連付けられたパラメータを調整可能な選択手段とを備え、
前記制御手段は、前記判定部において前記シミュレーションデータが前記実験データを導出可能であると判定された場合に前記選択手段による前記素反応式の選択及び前記パラメータの調整を実行させる
ことを特徴とする請求項１〜３いずれかに記載の反応メカニズム解析装置。 - コンピュータを、
流通方向に対してのみ物質移動及び熱移動が起こりうると近似された内部流路を有する管の一端側から複数の前記初期物質が混合された予混合ガスを導入して化学反応させる実験より得られた前記流通方向における位置に対する温度分布を実験データとして記憶する実験データ記憶手段と、
前記化学反応を記述する複数の素反応式及び該素反応式に関連付けられたパラメータを素反応データとして記憶する素反応データ記憶手段と、
前記素反応データからシミュレーションデータを算出するための方程式であって、空間座標に基づかず、時間に対する変数を含むと共に温度の次元を有する方程式を演算データとして記憶する演算データ記憶手段と、
前記素反応データ及び前記演算データとして記憶された方程式を用いて、前記素反応によって生成される中間生成物と前記初期物質との混合体の温度変化を前記シミュレーションデータとして算出する算出部と、
前記混合体の移動速度と、前記混合体の温度と、前記予混合ガスの供給速度と、前記予混合ガスの初期温度と、前記流通方向の位置と、時間との関係式に基づいて、前記シミュレーションデータから、前記管に前記予混合ガスを導入して化学反応させたときの前記流通方向における位置に対する温度分布を算出してシミュレーションデータに基づく温度分布とし、前記シミュレーションデータに基づく温度分布と前記実験データに基づく温度分布との差異が許容範囲にあることを導出可能であると定義し、前記シミュレーションデータが前記実験データを導出可能であるかを判定する判定部と、
前記素反応データに含まれる前記パラメータを調整可能なパラメータ調整手段と、
前記判定部において前記シミュレーションデータが前記実験データを導出可能であると判定されるまで、前記パラメータ調整手段による前記パラメータの調整及び調整後の前記パラメータを用いた前記算出部による前記シミュレーションデータの算出を繰り返し実行させる制御手段と
して機能させる反応メカニズム解析プログラム。 - 前記コンピュータを前記算出部として機能させる場合に、
中間生成物の濃度変化を示す項を含む質量保存の式と、素反応によって生成される中間生成物と前記初期物質との混合体の内部エネルギを示す項、前記素反応により移動するエネルギを示す項、及び前記混合体と外部との熱伝達により移動するエネルギを示す項を含むエネルギ保存の式とを前記演算データとして記憶された方程式として用いて前記シミュレーションデータを算出させることを特徴とする請求項５記載の反応メカニズム解析プログラム。 - 前記コンピュータを判定部として機能させる場合に、
複数の前記初期物質として燃料と酸化剤とが混合された予混合ガスを、前記予混合ガスの着火温度以上まで連続的に昇温するように加熱された管に、火炎形成位置が安定する流量で供給した際に得られる温度分布を、前記実験データとして記憶される前記流通方向における位置に対する温度分布として用いることを特徴とする請求項５または６記載の反応メカニズム解析プログラム。 - 前記コンピュータを、解析された素反応を基に解析を行う際の解析条件を記憶する解析条件記憶手段と、前記化学反応を記述する複数の素反応式を前記解析条件に基づいて選択すると共に選択された素反応式に関連付けられたパラメータを調整可能な選択手段として機能させ、
前記判定部において前記シミュレーションデータが前記実験データを導出可能であると判定された場合に前記コンピュータを前記選択手段として機能させる
ことを特徴とする請求項５〜７いずれかに記載の反応メカニズム解析プログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2008287099A JP5051659B2 (ja) | 2008-11-07 | 2008-11-07 | 反応メカニズム解析装置及び反応メカニズム解析プログラム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2008287099A JP5051659B2 (ja) | 2008-11-07 | 2008-11-07 | 反応メカニズム解析装置及び反応メカニズム解析プログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2010110724A JP2010110724A (ja) | 2010-05-20 |
| JP5051659B2 true JP5051659B2 (ja) | 2012-10-17 |
Family
ID=42299729
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2008287099A Active JP5051659B2 (ja) | 2008-11-07 | 2008-11-07 | 反応メカニズム解析装置及び反応メカニズム解析プログラム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5051659B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9523668B2 (en) | 2010-01-19 | 2016-12-20 | Tohoku University | Fuel property determination method and fuel property determination device |
| US9562868B2 (en) | 2010-11-18 | 2017-02-07 | Tohoku University | Combustion experimental apparatus |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPWO2012121281A1 (ja) | 2011-03-07 | 2014-07-17 | 国立大学法人東北大学 | 燃料物性決定方法及び燃料物性決定装置 |
| US20180032704A1 (en) * | 2015-02-19 | 2018-02-01 | Kabushiki Kaisha Toyota Chuo Kenkyusho | Reaction mechanism generation method and reaction mechanism generation device |
| US20200306960A1 (en) * | 2019-04-01 | 2020-10-01 | Nvidia Corporation | Simulation of tasks using neural networks |
| JP7276301B2 (ja) * | 2020-03-31 | 2023-05-18 | 横河電機株式会社 | 反応解析装置、反応解析システム、および反応解析方法 |
| CN114171127B (zh) * | 2021-12-07 | 2024-09-10 | 上海交通大学四川研究院 | 一种构建自适应反应机理模拟超声速燃烧流动的方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000097422A (ja) * | 1998-09-22 | 2000-04-04 | Hitachi Ltd | 廃棄物焼却プラント及びその制御装置並びに制御方法と廃棄物焼却炉のガス組成分布予測方法 |
| JP4077659B2 (ja) * | 2002-05-23 | 2008-04-16 | 財団法人電力中央研究所 | 燃焼装置における触媒の寿命予測方法 |
| JP2004008908A (ja) * | 2002-06-06 | 2004-01-15 | Mazda Motor Corp | 流体中の所定成分の濃度変化予測方法及び反応速度解析法 |
-
2008
- 2008-11-07 JP JP2008287099A patent/JP5051659B2/ja active Active
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9523668B2 (en) | 2010-01-19 | 2016-12-20 | Tohoku University | Fuel property determination method and fuel property determination device |
| US9562868B2 (en) | 2010-11-18 | 2017-02-07 | Tohoku University | Combustion experimental apparatus |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2010110724A (ja) | 2010-05-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5526408B2 (ja) | 燃料物性決定方法及び燃料物性決定装置 | |
| JP5051659B2 (ja) | 反応メカニズム解析装置及び反応メカニズム解析プログラム | |
| Pierce et al. | Progress-variable approach for large-eddy simulation of non-premixed turbulent combustion | |
| Bolla et al. | Soot formation modeling of n-heptane sprays under diesel engine conditions using the conditional moment closure approach | |
| Lv et al. | Discontinuous Galerkin method for multicomponent chemically reacting flows and combustion | |
| Senecal et al. | Grid-convergent spray models for internal combustion engine computational fluid dynamics simulations | |
| Som et al. | Simulating flame lift-off characteristics of diesel and biodiesel fuels using detailed chemical-kinetic mechanisms and large eddy simulation turbulence model | |
| Yu et al. | Onset of cellular flame instability in adiabatic CH4/O2/CO2 and CH4/air laminar premixed flames stabilized on a flat-flame burner | |
| Franzelli et al. | Large eddy simulation of swirled spray flame using detailed and tabulated chemical descriptions | |
| Jin et al. | Cellular instabilities of non-adiabatic laminar flat methane/hydrogen oxy-fuel flames highly diluted with CO2 | |
| Sun et al. | Modelling of soot coalescence and aggregation with a two-population balance equation model and a conservative finite volume method | |
| Volchkov et al. | Characterization of the flame blow-off conditions in a laminar boundary layer with hydrogen injection | |
| Sarraf et al. | Quantitative analysis of flame instabilities in a Hele-Shaw burner | |
| Swaminathan et al. | On the correlation of heat release rate in turbulent premixed flames | |
| Liu et al. | Surrogate modeling of parameterized multi-dimensional premixed combustion with physics-informed neural networks for rapid exploration of design space | |
| Bao et al. | The inclusion of scalar dissipation rate in modeling of an n-dodecane spray flame using flamelet generated manifold | |
| Morsy et al. | Numerical and experimental study on turbulence statistics in a large fan-stirred combustion vessel | |
| Zheng et al. | Spectral radiation properties of partially premixed turbulent flames | |
| Lammers et al. | The influence of gas radiation on the temperature decrease above a burner with a flat porous inert surface | |
| Piehl et al. | Turbulent spray combustion modeling using various kinetic solvers and turbulence models | |
| Abdul-Sater et al. | Predicting Radiative Heat Transfer in Oxy‐Methane Flame Simulations: An Examination of Its Sensitivities to Chemistry and Radiative Property Models | |
| Han et al. | Conditional moment closure method with tabulated chemistry in adiabatic turbulent nonpremixed jet flames | |
| Alhameedi et al. | Toward a Better Air-Assisted Flare Design for Purge Flow Conditions: Experimental and Computational Investigation of Radial Slot Flow into a Crossflow Environment | |
| Pang et al. | Development and validation of a local time stepping-based PaSR solver for combustion and radiation modeling | |
| Syed et al. | A Novel Approach to Predicting NO x Emissions From Dry Low Emissions Gas Turbines |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20110906 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20110906 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20120220 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120327 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120525 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20120619 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20120717 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 5051659 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20150803 Year of fee payment: 3 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |