(57)【要約】【課題】  内燃機関の筒内（燃焼室内）における燃焼温度を簡易な構成で精度良く推定し得る内燃機関の燃焼温度推定方法を提供すること。【解決手段】  この方法は、燃料噴射開始時期が到来する毎に、筒内ガスの状態が断熱変化することを原則的に利用して着火時における燃焼前の筒内ガス温度（着火時圧縮筒内ガス温度Tpump）を推定し、燃料の燃焼による発熱量を吸気中のガス成分の濃度比率から求められる燃焼に関与する燃焼後の筒内ガスのモル数及び定圧比熱の積で除することで燃焼による筒内ガスの温度上昇量（燃焼起因温度上昇量ΔTburn）を推定し、筒内の燃焼速度に影響を与える因子である燃料噴射圧力及びエンジン回転速度に基づいて同燃焼速度の増大による筒内ガスの温度上昇量（燃焼速度起因温度上昇量ΔTb_velo）を推定する。そして、Tflame＝Tpump＋ΔTburn＋ΔTb_veloなる式にて最高燃焼温度Tflameを推定する。【選択図】  図１１
2005180220.tif 000002
【特許請求の範囲】
【請求項１】
内燃機関の筒内における燃焼温度を推定する内燃機関の燃焼温度推定方法であって、
少なくとも前記筒内に存在する筒内ガスが同筒内で圧縮されることを利用して、着火時における燃焼前の筒内ガス温度である着火時圧縮筒内ガス温度を推定し、
少なくとも前記筒内に吸入されるガスの組成と、噴射された燃料の燃焼による発熱量とに基づいて、燃焼による前記筒内ガスの温度上昇量である燃焼起因温度上昇量を推定し、
前記筒内の燃焼速度に影響を与える因子に基づいて同燃焼速度の増大による前記筒内ガスの温度上昇量である燃焼速度起因温度上昇量を推定するとともに、
前記着火時圧縮筒内ガス温度に、前記燃焼起因温度上昇量と前記燃焼速度起因温度上昇量とを加えた値を前記筒内における燃焼温度として推定する内燃機関の燃焼温度推定方法。
【請求項２】
請求項１に記載の内燃機関の燃焼温度推定方法において、
前記筒内に吸入されるガスの温度と、同吸入されるガスの組成と、前記燃焼前の筒内ガス温度を上昇せしめる因子に基づいて推定される前記着火時圧縮筒内ガス温度の上昇量と、
に基づいて同着火時圧縮筒内ガス温度を推定する内燃機関の燃焼温度推定方法。
【請求項３】
請求項２に記載の内燃機関の燃焼温度推定方法において、
前記燃焼前の筒内ガス温度を上昇せしめる因子としての、前記燃料の噴射に先立ってパイロット噴射が行われる場合における同パイロット噴射された燃料の燃焼による発熱量、及びグロープラグへの通電が行われる場合における同通電による発熱量のうち少なくとも１つ、
に基づいて前記着火時圧縮筒内ガス温度の上昇量を推定する内燃機関の燃焼温度推定方法。
【請求項４】
請求項１乃至請求項３の何れか一項に記載の内燃機関の燃焼温度推定方法において、
前記着火時が圧縮上死点に対応する時期よりも前である場合、同着火時が同圧縮上死点に対応する時期と一致するものと仮定して、前記着火時圧縮筒内ガス温度を推定する内燃機関の燃焼温度推定方法。
【請求項５】
請求項１乃至請求項４の何れか一項に記載の内燃機関の燃焼温度推定方法において、
前記筒内の燃焼速度に影響を与える因子としての、前記燃料の噴射圧力、前記機関の回転速度、前記筒内に吸入されるガスのスワール比、前記機関に過給機が配設されている場合における同過給機による過給圧、及び前記筒内に吸入されるガス中の酸素濃度のうち少なくとも１つ、
に基づいて前記燃焼速度起因温度上昇量を推定する内燃機関の燃焼温度推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内燃機関の筒内（燃焼室内）における燃焼温度を推定する燃焼温度推定方法に関する。
【背景技術】
【０００２】
火花点火式内燃機関、ディーゼル機関等の内燃機関から排出されるＮＯｘ等のエミッションの量は、筒内における燃焼温度（最高燃焼温度、最高火炎温度）に強い相関を有している。従って、ＮＯｘ等のエミッションの排出量を低減するためには燃焼温度を所定温度に制御することが有効である。一方、この燃焼温度を実測することは非常に困難である。よって、燃焼温度を所定温度に制御するためには同燃焼温度を精度良く推定する必要がある。
【０００３】
このため、下記特許文献１に記載の内燃機関の燃焼室内温度状態推定装置は、機関の燃焼室から排出される排ガスの温度を同燃焼室外で検出する排気温度センサを備え、この排気温度センサが検出した燃焼室外排気温度検出値に基づいて、この検出値と相関性の高い燃焼終了後の燃焼室内温度（従って、上記筒内における燃焼温度）を推定するようになっている。
【特許文献１】特開２００２−５４４９１号公報
【０００４】
上記装置においては、筒内における燃焼温度の推定に際し、燃焼室外排気温度と同燃焼温度との間に強い相関があることを前提としている。しかしながら、実際には両者間に常に強い相関があるわけではないから、この結果、上記燃焼温度を精度良く推定できない場合がある。また、排気温度センサが必須の構成となっているため、製造コストが増大し、且つ構成が複雑化するという問題がある。
【発明の開示】
【０００５】
本発明は、かかる課題に対処するためになされたものであって、その目的は、内燃機関の筒内（燃焼室内）における燃焼温度を簡易な構成で精度良く推定し得る内燃機関の燃焼温度推定方法を提供することにある。
【０００６】
本発明による内燃機関の燃焼温度推定方法では、先ず、少なくとも筒内に存在する筒内ガスが同筒内で圧縮されることを利用して、着火時における燃焼前の筒内ガス温度である着火時圧縮筒内ガス温度が推定される。ここにおいて、着火時は、例えば、火花点火式内燃機関については点火プラグの点火時であり、ディーゼル機関については燃料噴射時期（少なくとも１回のパイロット噴射の後にメイン噴射が実行される場合においてはメイン噴射時期）から所定の着火遅れ時間が経過した時点である。
【０００７】
この着火時圧縮筒内ガス温度は、例えば、圧縮行程（及び膨張行程）における燃焼前の筒内ガスの状態（即ち、温度、圧力）が断熱変化するとの仮定のもと、ポリトロープ指数を用いた断熱変化を表す一般的な式と、上記着火時における筒内容積とに基づいて正確、且つ容易に求めることができる。
【０００８】
また、この推定方法では、少なくとも前記筒内に吸入されるガスの組成と、噴射された燃料の燃焼による発熱量とに基づいて、燃焼による前記筒内ガスの温度上昇量である燃焼起因温度上昇量が推定される。一般に、燃焼による筒内ガスの温度上昇量（従って、燃焼起因温度上昇量）は、燃料の燃焼による発熱量を、燃焼後の同燃焼に関与した筒内ガスの比熱（定圧比熱）と、燃焼後の同燃焼に関与した筒内ガスの量（モル数）とで除することで求めることができる。
【０００９】
ここで、上記燃焼に関与した筒内ガスの定圧比熱、及びモル数は、筒内に吸入されるガス（以下、「吸気」と云うこともある。）の組成（吸気を構成する複数成分（例えば、酸素、不活性ガス）の濃度割合）により変化する値であって、例えば、吸気中の不活性ガスの濃度割合の増大により共に増大する値である（詳細は後述する。）。従って、この燃焼起因温度上昇量は、吸気の組成（例えば、吸気中の不活性ガスの濃度割合等）に基づいて正確に求められ得る上記燃焼に関与した筒内ガスの定圧比熱及びモル数、並びに燃料の燃焼による発熱量に基づいて精度良く求めることができる。
【００１０】
更に、この推定方法では、筒内の燃焼速度に影響を与える因子に基づいて同燃焼速度の増大による前記筒内ガスの温度上昇量である燃焼速度起因温度上昇量が推定される。一般に、ディーゼル機関においては、着火時が圧縮上死点に対応する時期よりも後の時点（即ち、膨張行程に対応する時点）である場合が多い。膨張行程では、筒内ガス温度は時間経過に従って低下する。
【００１１】
従って、この場合、着火時点（即ち、燃焼開始時点）から筒内ガス温度が最高燃焼温度に達するまでの時間が短いほど同最高燃焼温度が高くなる。換言すれば、燃焼開始後における筒内の燃焼速度が大きいほど筒内ガスの（最高）燃焼温度が高くなる。かかる筒内ガスの温度上昇量が前記筒内の燃焼速度に影響を与える因子に基づいて燃焼速度起因温度上昇量として求められる。
【００１２】
前記筒内の燃焼速度に影響を与える因子としては、例えば、前記燃料の噴射圧力、前記機関の回転速度、前記筒内に吸入されるガスのスワール比、前記機関に過給機が配設されている場合における同過給機による過給圧、及び前記筒内に吸入されるガス中の酸素濃度が挙げられる。
【００１３】
そして、この推定方法では、前記着火時圧縮筒内ガス温度に、前記燃焼起因温度上昇量と前記燃焼速度起因温度上昇量とを加えた値を筒内における（最高）燃焼温度として推定される。このようにして推定される筒内における最高燃焼温度は、以下のごとく、種々の実際の現象を精度良く表す値となり得る。
【００１４】
即ち、例えば、燃料の噴射時期を遅らせることで（膨張行程にある）着火時が遅れると、同着火時の筒内容積の増大により上記推定される着火時圧縮筒内ガス温度が低下し、この結果、上記推定される筒内における最高燃焼温度が低下する。この推定結果は、燃料噴射時期を遅らせるとＮＯｘ発生量が減少するという実際の現象に沿うものである。
【００１５】
また、例えば、ＥＧＲ率（吸気量に対するＥＧＲガス量の割合）の増大により吸気中の不活性ガス（例えば、ＣＯ２等）の濃度を増大させると、上述のごとく上記燃焼に関与した筒内ガスの定圧比熱及びモル数が共に増大するから上記推定される燃焼起因温度上昇量が低下し、この結果、上記推定される筒内における最高燃焼温度が低下する。この推定結果は、ＥＧＲ率を増大させるとＮＯｘ発生量が減少するという実際の現象に沿うものである。
【００１６】
また、例えば、燃料の噴射圧力を増大させ、且つ機関の回転速度を大きくすることで筒内の燃焼速度を速くすると、上記推定される燃焼速度起因温度上昇量が増大し、この結果、上記推定される筒内における最高燃焼温度が増大する。この推定結果は、燃料の噴射圧力を増大させ、且つ機関の回転速度を大きくするとＮＯｘ発生量が増加するという実際の現象に沿うものである。以上、説明したように、上記本発明による燃焼温度推定方法によれば、筒内における最高燃焼温度を簡易な構成で実際の種々の現象に沿うように精度良く推定することができる。
【００１７】
上記本発明による燃焼温度推定方法においては、前記筒内に吸入されるガスの温度と、同吸入されるガスの組成と、前記燃焼前の筒内ガス温度を上昇せしめる因子とに基づいて推定される前記着火時圧縮筒内ガス温度の上昇量と、に基づいて同着火時圧縮筒内ガス温度を推定することが好適である。この場合、前記燃焼前の筒内ガス温度を上昇せしめる因子としては、例えば、前記燃料の噴射に先立ってパイロット噴射が行われる場合における同パイロット噴射された燃料の燃焼による発熱量、及びグロープラグへの通電が行われる場合における同通電による発熱量が挙げられる。
【００１８】
上記着火時圧縮筒内ガス温度は、当然に吸気温度により変化する。また、断熱変化する筒内ガスに使用すべきポリトロープ指数は吸気の組成（従って、筒内ガスの組成）により変化するから、（ポリトロープ指数を用いた断熱変化を表す一般的な式に基づく）上記着火時圧縮筒内ガス温度は、吸気の組成によっても変化する。
【００１９】
更には、燃焼前の筒内ガスに上記グロープラグへの通電等による熱量が与えられる場合、上記着火時圧縮筒内ガス温度は、同熱量に相当する分だけ上昇する。かかる筒内ガスの温度上昇量が上記着火時圧縮筒内ガス温度の上昇量として求められる。
【００２０】
以上のことから、筒内に吸入されるガス（吸気）の温度、同吸入されるガス（吸気）の組成、及び、上記着火時圧縮筒内ガス温度の上昇量は、上記着火時圧縮筒内ガス温度に影響を与える値（パラメータ）となり得る。従って、上記のように、上記着火時圧縮筒内ガス温度の推定に際し、筒内ガスが筒内で圧縮されることに加えて上記３つのパラメータを考慮することで同着火時圧縮筒内ガス温度がより一層精度良く推定され得、この結果、筒内における（最高）燃焼温度もより一層精度良く推定され得る。
【００２１】
また、上記本発明による燃焼温度推定方法においては、前記着火時が圧縮上死点に対応する時期よりも前である場合、同着火時が同圧縮上死点に対応する時期と一致するものと仮定して、前記着火時圧縮筒内ガス温度を推定することが好適である。
【００２２】
着火時が圧縮上死点に対応する時期よりも前である場合（即ち、着火時が圧縮行程に対応する時点である場合）、着火後の（即ち、燃焼中、或いは燃焼後の）筒内ガスそのものが圧縮上死点に対応する時期まで更に圧縮され続け、この結果、筒内における最高燃焼温度も圧縮上死点に対応する時期まで上昇していくと考えることができる。換言すれば、この場合における筒内における最高燃焼温度は、着火時が同圧縮上死点に対応する時期と一致する場合における最高燃焼温度と一致する。
【００２３】
従って、上記のように、着火時が圧縮上死点に対応する時期よりも前である場合、同着火時が同圧縮上死点に対応する時期と一致するものと仮定して着火時圧縮筒内ガス温度を推定することで、着火時が圧縮上死点に対応する時期よりも前である場合における筒内における最高燃焼温度が、更に一層精度良く推定され得る。
【発明を実施するための最良の形態】
【００２４】
以下、本発明の実施形態に係る内燃機関の燃焼温度推定方法を実施するとともに、同方法により推定された燃焼温度を利用して筒内にて燃焼により発生するＮＯｘ発生量を推定する内燃機関（ディーゼル機関）の制御装置について図面を参照しつつ説明する。
【００２５】
図１は、係る内燃機関の制御装置を４気筒内燃機関（ディーゼル機関）１０に適用したシステム全体の概略構成を示している。このシステムは、燃料供給系統を含むエンジン本体２０、エンジン本体２０の各気筒の燃焼室（筒内）にガスを導入するための吸気系統３０、エンジン本体２０からの排ガスを放出するための排気系統４０、排気還流を行うためのＥＧＲ装置５０、及び電気制御装置６０を含んでいる。
【００２６】
エンジン本体２０の各気筒の上部には燃料噴射弁（噴射弁、インジェクタ）２１が配設されている。各燃料噴射弁２１は、図示しない燃料タンクと接続された燃料噴射用ポンプ２２に燃料配管２３を介して接続されている。また、各気筒の上部には、燃料噴射弁２１に隣接するようにグロープラグ２４が配設されている。各グロープラグ２４は、電気制御装置６０と電気的に接続されていて、機関が暖機運転状態等の所定の運転状態にあるときにのみ同電気制御装置６０からの信号により通電されて発熱し、各筒内に存在する筒内ガスに所定量の熱を供給するようになっている。
【００２７】
燃料噴射用ポンプ２２は、電気制御装置６０と電気的に接続されていて、同電気制御装置６０からの駆動信号（後述する（指令）基本燃料噴射圧力Pcrbaseに応じた指令信号）により燃料の実際の噴射圧力（吐出圧力）が同指令基本燃料噴射圧力Pcrbaseになるように同燃料を昇圧するようになっている。
【００２８】
これにより、燃料噴射弁２１には、燃料噴射用ポンプ２２から前記基本燃料噴射圧力Pcrbaseまで昇圧された燃料が供給されるようになっている。また、燃料噴射弁２１は、電気制御装置６０と電気的に接続されていて、同電気制御装置６０からの駆動信号（後述する（指令）燃料噴射量（質量）qfinに応じた指令信号）により所定時間だけ開弁し、これにより各気筒の燃焼室内に前記指令基本燃料噴射圧力Pcrbaseにまで昇圧された燃料を前記指令燃料噴射量qfinだけ直接噴射するようになっている。
【００２９】
吸気系統３０は、エンジン本体２０の各気筒の燃焼室にそれぞれ接続された吸気マニホールド３１、吸気マニホールド３１の上流側集合部に接続され同吸気マニホールド３１とともに吸気通路を構成する吸気管３２、吸気管３２内に回動可能に保持されたスロットル弁３３、電気制御装置６０からの駆動信号に応答してスロットル弁３３を回転駆動するスロットル弁アクチュエータ３３ａ、スロットル弁３３の上流において吸気管３２に順に介装されたインタクーラー３４と過給機３５のコンプレッサ３５ａ、及び吸気管３２の先端部に配設されたエアクリーナ３６とを含んでいる。
【００３０】
排気系統４０は、エンジン本体２０の各気筒にそれぞれ接続された排気マニホールド４１、排気マニホールド４１の下流側集合部に接続された排気管４２、排気管４２に配設された過給機３５のタービン３５ｂ、及び排気管４２に介装されたディーゼルパティキュレートフィルタ（以下、「ＤＰＮＲ」と称呼する。）４３を含んでいる。排気マニホールド４１及び排気管４２は排気通路を構成している。
【００３１】
ＤＰＮＲ４３は、コージライト等の多孔質材料から形成されたフィルタ４３ａを備え、通過する排気ガス中のパティキュレートを細孔表面にて捕集するフィルタである。ＤＰＮＲ４３は、担体としてのアルミナに、カリウムＫ，ナトリウムＮａ，リチウムＬｉ，セシウムＣｓのようなアルカリ金属、バリウムＢａ，カルシウムＣａのようなアルカリ土類金属、及びランタンＬａ、イットリウムＹのような希土類金属から選ばれた少なくとも一つを白金とともに担持し、ＮＯｘを吸収した後に同吸収したＮＯｘを放出して還元する吸蔵還元型ＮＯｘ触媒としても機能するようになっている。
【００３２】
ＥＧＲ装置５０は、排気ガスを還流させる通路（ＥＧＲ通路）を構成する排気還流管５１と、排気還流管５１に介装されたＥＧＲ制御弁５２と、ＥＧＲクーラー５３とを備えている。排気還流管５１はタービン３５ｂの上流側排気通路（排気マニホールド４１）とスロットル弁３３の下流側吸気通路（吸気マニホールド３１）を連通している。ＥＧＲ制御弁５２は電気制御装置６０からの駆動信号に応答し、再循環される排気ガス量（排気還流量、ＥＧＲガス流量）を変更し得るようになっている。
【００３３】
電気制御装置６０は、互いにバスで接続されたＣＰＵ６１、ＣＰＵ６１が実行するプログラム、テーブル（ルックアップテーブル、マップ）、及び定数等を予め記憶したＲＯＭ６２、ＣＰＵ６１が必要に応じてデータを一時的に格納するＲＡＭ６３、電源が投入された状態でデータを格納するとともに同格納したデータを電源が遮断されている間も保持するバックアップＲＡＭ６４、並びにＡＤコンバータを含むインターフェース６５等からなるマイクロコンピュータである。
【００３４】
インターフェース６５は、空気流量（新気流量）計測手段であって吸気管３２に配置された熱線式エアフローメータ７１、スロットル弁３３の下流であって排気還流管５１が接続された部位よりも下流の吸気通路に設けられた吸気温センサ７２、スロットル弁３３の下流であって排気還流管５１が接続された部位よりも下流の吸気通路に配設された吸気管圧力センサ７３、クランクポジションセンサ７４、アクセル開度センサ７５、及び、スロットル弁３３の下流であって排気還流管５１が接続された部位よりも下流の吸気通路に設けられた吸気酸素濃度センサ７６と接続されていて、これらのセンサからの信号をＣＰＵ６１に供給するようになっている。また、インターフェース６５は、燃料噴射弁２１、燃料噴射用ポンプ２２、スロットル弁アクチュエータ３３ａ、及びＥＧＲ制御弁５２と接続されていて、ＣＰＵ６１の指示に応じてこれらに駆動信号を送出するようになっている。
【００３５】
熱線式エアフローメータ７１は、吸気通路内を通過する吸入空気（新気）の質量流量（単位時間あたりの吸入新気量）を計測し、同質量流量Ga（吸入新気流量Ga）を表す信号を発生するようになっている。吸気温センサ７２は、前述した吸気の温度を検出し、同吸気温度Tbを表す信号を発生するようになっている。吸気管圧力センサ７３は、吸気の圧力（即ち、吸気管圧力）を検出し、同吸気管圧力Pbを表す信号を発生するようになっている。
【００３６】
クランクポジションセンサ７４は、各気筒の絶対クランク角度を検出し、クランク角度CAを表すとともにエンジン１０の回転速度であるエンジン回転速度NEをも表す信号を発生するようになっている。アクセル開度センサ７５は、アクセルペダルＡＰの操作量を検出し、アクセル操作量Accpを表す信号を発生するようになっている。吸気酸素濃度センサ７６は、吸気中の酸素濃度（即ち、吸気酸素濃度）を検出し、同吸気酸素濃度RO2_inを表す信号を発生するようになっている。
【００３７】
（燃焼温度推定方法の概要）
次に、上記のように構成された内燃機関の制御装置（以下、「本装置」と云うこともある。）による本発明の実施形態に係る燃焼温度推定方法の概要について説明する。図２は、機関１０の或る一つの気筒の筒内に吸気マニホールド３１からガス（即ち、吸気）が吸入され、筒内に吸入されたガス（筒内ガス）が燃焼後に排気マニホールド４１へ排出される様子を模式的に示した図である。
【００３８】
図２に示したように、吸気（従って、筒内ガス）には、吸気管３２の先端部からスロットル弁３３を介して吸入された新気と、排気還流管５１からＥＧＲ制御弁５２を介して吸入されたＥＧＲガスが含まれる。吸入される新気量（新気質量）と吸入されるＥＧＲガス量（ＥＧＲガス質量）の和に対するＥＧＲガス質量の割合（即ち、ＥＧＲ率）は、運転状態に応じて電気制御装置６０（ＣＰＵ６１）により適宜制御されるスロットル弁３３の開度、及びＥＧＲ制御弁５２の開度に応じて変化する。
【００３９】
吸気（即ち、新気、及びＥＧＲガスから構成されるガス）は、吸気行程において開弁している吸気弁Vinを介してピストンの下降に伴って筒内に吸入されて筒内ガスとなる。筒内ガスは、ピストンが下死点に達した時点（以下、「ATDC-180°」と称呼する。）で吸気弁Vinが閉弁することにより筒内に密閉され、その後の圧縮行程においてピストンの上昇に伴って圧縮される。そして、ピストンが圧縮上死点（以下、「ATDC0°」と称呼する。）近傍に達っすると（具体的には、後述する最終燃料噴射時期finjfinが到来すると）、本装置は、前記指令燃料噴射量qfinに応じた所定時間だけ燃料噴射弁２１を開弁することで燃料を筒内に直接噴射する。この結果、噴射された燃料は、時間の経過に伴って同筒内ガスと混ざり合いながら混合気となって筒内において拡散していき、所定の着火遅れ時間が経過すると自己着火が発生することに起因して燃焼する。
【００４０】
係る燃焼は、本実施例では、後述するように推定される燃焼領域（以下、「Ｂ領域」と云うこともある。図２を参照。）においてのみ発生し、燃焼室内におけるＢ領域を除いた領域である非燃焼領域（以下、「Ａ領域」と云うこともある。図２を参照。）では発生しないと仮定する。そして、燃焼後に燃焼室内に存在する筒内ガスは、排ガスとなって、排気行程において開弁している排気弁Voutを介してピストンの上昇に伴って排気マニホールド４１へ排出される。係る排ガスは、その一部が排気還流管５１を介してＥＧＲガスとなって吸気側に還流され、残りの排ガスは排気管４２を介して外部へと排出されていく。
【００４１】
以下、本装置が実施する具体的な燃焼温度推定方法について説明していく。この燃焼温度推定方法では、燃料が噴射される対象となる気筒（以下、「燃料噴射気筒」と称呼する。）について上記最終燃料噴射時期finjfinが到来する毎に、その直後（上記所定の着火遅れ時間経過後）において上記Ｂ領域にて発生する燃焼による筒内ガスの最高燃焼温度Tflameが推定されていく。本装置は、下記(1)式に従って最高燃焼温度Tflameを求める。
【００４２】
Tflame＝Tpump＋ΔTburn＋ΔTb_velo  ・・・(1)
【００４３】
上記(1)式において、Tpumpは着火時における燃焼前の筒内ガス温度である着火時圧縮筒内ガス温度である。ΔTburnは燃焼による筒内ガスの温度上昇量である燃焼起因温度上昇量である。ΔTb_veloは燃焼速度の増大による筒内ガスの温度上昇量である。以下、これらの値の取得方法について個別に説明する。
【００４４】
<着火時圧縮筒内ガス温度Tpumpの取得>
着火時圧縮筒内ガス温度Tpumpを取得するため、先ず、圧縮行程、及び膨張行程において筒内ガスと外部との熱交換がないと仮定するものとする。そうすると、係る筒内ガスの状態は断熱変化するから、筒内ガス温度Taは、図３に実線にて示すように、クランク角度CA(ATDC)に応じて変化し、着火時の筒内ガス温度（即ち、着火時圧縮筒内ガス温度Tpump）は、断熱変化を表す一般的な式である下記(2)式に従って求めることができる。
【００４５】
Tpump＝Ta0・(Va0/Vig)κ-1  ・・・(2)
【００４６】
上記(2)式において、Ta0は、ATDC-180°における下死点時筒内ガス温度である。ATDC-180°において筒内ガス温度は吸気温度Tbと略等しいと考えられるから、下死点時筒内ガス温度Ta0は、ATDC-180°において吸気温センサ７２により検出される吸気温度Tbとして取得することができる。Va0は、ATDC-180°における下死点時筒内容積である。筒内容積Vaは機関１０の設計諸元に基づいてクランク角度CAの関数Va(CA)として表すことができるから、この関数に基づいて下死点時筒内容積Va0も求めることができる。
【００４７】
また、上記(2)式において、Vigは着火時における筒内容積である。ここで、着火時は燃料噴射時期から所定の着火遅れ時間が経過した時点であるから、図３に示したように、着火時におけるクランク角度CAigは、上記最終燃料噴射時期finjfinに相当する燃料噴射時クランク角度CAinjから上記着火遅れ時間に相当するクランク角度ΔCAdelayだけ遅角させた値として求めることができる。従って、上記着火時筒内容積Vigは、Va(CAig)として求めることができる。
【００４８】
また、上記(2)式において、κはポリトロープ指数である。ここで、断熱変化する筒内ガスに使用すべきポリトロープ指数κは吸気の組成（従って、筒内ガスの組成）により変化する。従って、本例では、吸気酸素濃度センサ７６により検出された吸気酸素濃度RO2_in（具体的には、下死点時（ATDC-180°）に検出された吸気酸素濃度RO2_inである下死点時吸気酸素濃度RO2c）と、吸気酸素濃度からポリトロープ指数を求める関数ｇとに基づいて、ポリトロープ指数κは、ｇ(RO2c)として求めることができる。
【００４９】
以上のようにして、本装置は、原則的に、着火時圧縮筒内ガス温度Tpumpを上記(2)式に従って求める。この結果、図３に実線にて示したように、着火時圧縮筒内ガス温度Tpumpは、着火時（着火時クランク角度CAig）が圧縮上死点（ATDC０°）に対応する時期よりも前である場合（即ち、着火時が圧縮行程に対応する時点である場合）には同着火時が遅れるほど上昇し、一方、着火時が圧縮上死点に対応する時期よりも後である場合（即ち、着火時が膨張行程に対応する時点である場合）には同着火時が遅れるほど減少する。
【００５０】
ところで、着火時が圧縮上死点に対応する時期よりも前である場合、着火後の（即ち、燃焼中、或いは燃焼後の）筒内ガスそのものが圧縮上死点に対応する時期まで更に圧縮され続け、この結果、最終的な演算対象である上記(1)式に従う上記最高燃焼温度Tflameも圧縮上死点に対応する時期まで上昇していくと考えられる。即ち、この場合、上記最高燃焼温度Tflameは、着火時が圧縮上死点に対応する時期と一致する場合における最高燃焼温度と一致させる方が好ましい。
【００５１】
そこで、着火時（即ち、着火時クランク角度CAig）が圧縮上死点に対応する時期よりも前である場合、着火時が圧縮上死点に対応する時期と一致するものと仮定し、本装置は、着火時筒内容積Vigの値をVa(CAig)に代えてATDC0°における筒内容積Vtop（一定値）として上記(2)式に従って着火時筒内ガス温度Tpumpを求める。この結果、図３に破線にて示したように、この場合の着火時圧縮筒内ガス温度Tpumpは、圧縮上死点における筒内ガス温度Ta1（一定値）として推定される。
【００５２】
以上、圧縮行程、及び膨張行程において筒内ガスと外部との熱交換がないと仮定した場合について説明した。しかしながら、実際には、上記グロープラグ２４が通電されている場合、圧縮行程、及び膨張行程にある筒内ガスは同グロープラグ２４から所定量の熱を受け、この結果、着火時の筒内ガス温度（従って、着火時圧縮筒内ガス温度Tpump）はこの熱量に相当する分だけ上昇する。
【００５３】
そこで、本装置は、グロープラグ２４が通電されている場合、上記のように求めた着火時圧縮筒内ガス温度Tpumpに、上記熱量に相当する筒内ガス温度上昇量ΔTpump（＝Tglow（本例では、一定値））を加えた値を最終的な着火時圧縮筒内ガス温度として推定する。以上のように、着火時圧縮筒内ガス温度Tpumpは、少なくとも、断熱圧縮に係る上記(2)式と、下死点時（ATDC-180°）における吸気温度Tbと、吸気の組成（上記下死点時吸気酸素濃度RO2c）と、上記筒内ガス温度上昇量ΔTpumpとに基づいて、上記最終燃料噴射時期finjfinが到来する毎に推定される。
【００５４】
<燃焼起因温度上昇量ΔTburnの取得>
燃焼による筒内ガスの温度上昇量である燃焼起因温度上昇量ΔTburnを取得するため、噴射された燃料の１molあたりの燃焼について考える。そうすると、燃焼起因温度上昇量ΔTburnは、下記(3)式にて表すことができる。
【００５５】
ΔTburn＝Qfuel/(Cp・ngas)  ・・・(3)
【００５６】
上記(3)式において、Qfuelは噴射された燃料１molあたりの燃焼による発熱量であって、燃料の組成により一義的に決定される既知の値（定数）である。Cpは１molの燃料の燃焼に関与する燃焼後の筒内ガスの定圧比熱である。ngasは１molの燃料の燃焼に関与する燃焼後の筒内ガスのモル数である。従って、上記(3)式に従って燃焼起因温度上昇量ΔTburnを求めるためには、上記定圧比熱Cp、及び上記モル数ngasを求める必要がある。ここで、吸気中の代表的なガス成分を含んだ燃料(ＣｍＨｎ)１molあたりの燃焼についての化学反応式は下記(4)式にて表される。
【００５７】
【数１】
2005180220.tif 000003
【００５８】
上記(4)式において、値ｍ、及び値ｎは噴射される燃料の組成により一義的に決定される値であって、例えば、ディーゼル機関に使用される標準的な燃料の場合、ｍ=15.75，ｎ=30.55である。α,β,γはそれぞれ、１molの燃料の燃焼に関与する不活性ガスであるＣＯ２，Ｈ２Ｏ，Ｎ２のモル数である。即ち、ここでは、吸気中のガス成分としてＯ２と、不活性ガスＣＯ２，Ｈ２Ｏ，Ｎ２の４つの成分に着目している。上記(4)式の右辺から容易に理解できるように、１molの燃料の燃焼に関与する燃焼後の筒内ガスの上記モル数ngas、及び上記定圧比熱Cp（等価定圧比熱）はそれぞれ、下記(5)式、及び下記(6)式にて表すことができる。下記(6)式において、Cp_CO2，Cp_H2O，Cp_N2，はそれぞれ、ＣＯ２，Ｈ２Ｏ，Ｎ２の定圧比熱である。
【００５９】
【数２】
2005180220.tif 000004
【００６０】
従って、上記(5)式、及び上記(6)式から理解できるように、上記モル数ngas、及び上記定圧比熱Cpを求めるためには、上記モル数α，β，γを求める必要がある。一方、これらの値は吸気の組成（具体的には、吸気を構成する上記４つの成分の吸気中の濃度（割合））に応じて変化するから、これらの値を求めるためには同４つの成分の吸気中の濃度（割合）を求める必要がある。そこで、以下、上記４つの成分の吸気中の各濃度の取得方法について説明する。
【００６１】
なお、以下において、[Ｘ]in、[Ｘ]egr、[Ｘ]air（Ｘ：Ｏ２，ＣＯ２，Ｈ２Ｏ，Ｎ２）は、成分Ｘの吸気中の質量濃度、成分ＸのＥＧＲガス中の質量濃度、成分Ｘの新気中の質量濃度をそれぞれ表すものとする。また、[Ｘ]airは既知の値であるものとする（例えば、[ＣＯ２]air＝[Ｈ２Ｏ]air＝０、[Ｏ２]air＝0.233、[Ｎ２]air＝0.767）。
【００６２】
また、以下において、Gcylは一吸気行程で筒内に吸入された吸気の総ガス質量（以下、「筒内総ガス量Gcyl」と称呼する。）であり、Gegrは一吸気行程で吸気の一部としてＥＧＲ装置５０から筒内に吸入されたＥＧＲガスの質量（以下、「ＥＧＲガス量」と称呼する。）であり、Gmは一吸気行程で吸気管３２の先端部から吸気の一部として筒内に吸入された新気の質量（以下、「吸入新気量」と称呼する。）である。
【００６３】
ここで、筒内総ガス量Gcylは、ATDC-180°における気体の状態方程式に基づく下記(7)式に従って求めることができる。
【００６４】
Gcyl＝(Pa0・Va0)/(R・Ta0)  ・・・(7)
【００６５】
上記(7)式において、Pa0は、ATDC-180°における下死点時筒内ガス圧力である。ATDC-180°において筒内ガス圧力は吸気管圧力Pbと略等しいと考えられるから、下死点時筒内ガス圧力Pa0は、ATDC-180°において吸気管圧力センサ７３により検出される吸気管圧力Pbとして取得することができる。Rは筒内ガスのガス定数である。Ta0、及びVa0はそれぞれ、上記(2)式と同様、下死点時筒内ガス温度、及び下死点時筒内容積である。
【００６６】
また、吸入新気量Gmは、エアフローメータ７１により計測される単位時間あたりの吸入新気量（吸入新気流量Ga）と、クランクポジションセンサ７４の出力に基づくエンジン回転速度NEと、吸入新気流量Ga及びエンジン回転速度NEを引数とする一吸気行程あたりの吸入新気量を求めるための関数ｆ(Ga,NE)とに基づいて求めることができる。吸入新気流量Ga及びエンジン回転速度NEとしては、ATDC-180°において各センサにより取得される下死点時吸入新気流量Ga0、及び下死点時エンジン回転速度NE0がそれぞれ使用される。
【００６７】
また、上記筒内総ガス量Gcylは上記ＥＧＲガス量Gegrと上記吸入新気量Gmの和であるから、ＥＧＲガス量Gegrは、上記のように求められ得る筒内総ガス量Gcyl、及び吸入新気量Gmと、下記(8)式とに基づいて求めることができる。
【００６８】
Gegr＝Gcyl-Gm  ・・・(8)
【００６９】
以下、先ず、吸気中のＣＯ２の濃度（以下、「吸気ＣＯ２濃度」と称呼する。）[ＣＯ２]inの取得方法について説明する。吸気ＣＯ２濃度は、筒内に吸入された燃焼前の筒内ガス中のＣＯ２濃度に等しいから、筒内総ガス量Gcylに対する、「筒内に吸入された新気中のＣＯ２質量
[ＣＯ２]air・Gm と、筒内に吸入されたＥＧＲガス中のＣＯ２質量 [ＣＯ２]egr・Gegr
の和」の質量割合として、下記(9)式に従って求めることができる。
【００７０】
【数３】
2005180220.tif 000005
【００７１】
ここで、ＥＧＲガス中のＣＯ２濃度[ＣＯ２]egrは（排気弁Voutを通過する）排ガス中のＣＯ２濃度に等しいと考えられ、同排ガス中のＣＯ２濃度は燃焼後の筒内ガス中のＣＯ２濃度である。燃焼後の筒内ガス中のＣＯ２濃度は、「筒内総ガス量Gcylに今回の指令燃料噴射量（質量）qfinc（＝qfin）を加えた値」に対する、「筒内に吸入された新気中のＣＯ２質量
[ＣＯ２]air・Gm と、筒内に吸入されたＥＧＲガス中のＣＯ２質量 [ＣＯ２]egr・Gegr
と、燃焼により発生したＣＯ２質量の和」の質量割合である。
【００７２】
燃焼により発生したＣＯ２質量は、単位質量あたりの燃料の燃焼により発生するＣＯ２質量をＫＣＯ２とすると、ＫＣＯ２・qfinc と表すことができる。また、上記(4)式から理解できるように燃料(ＣｍＨｎ)１molの燃焼によりｍ
molのＣＯ２が発生するから、燃料(ＣｍＨｎ)の分子量をＭfuel、ＣＯ２の分子量をＭＣＯ２（＝４４）とすると、質量Ｍfuelの燃料の燃焼により質量(ｍ・ＭＣＯ２)のＣＯ２が発生することになる。従って、ＫＣＯ２は、
ｍ・(ＭＣＯ２/Ｍfuel) として求めることができる。以上より、ＥＧＲガス中のＣＯ２濃度[ＣＯ２]egrは下記(10)式に従って求めることができ、下記(10)式を整理すると下記(11)式が得られる。
【００７３】
【数４】
2005180220.tif 000006
【００７４】
そして、上記(11)式を上記(9)式に代入すると、下記(12)式が得られ、下記(12)式に従って吸気ＣＯ２濃度[ＣＯ２]inを求めることができる。下記(12)式において、ＲはＥＧＲ率であって、Ｒ＝Gegr/Gcylである。
【００７５】
【数５】
2005180220.tif 000007
【００７６】
次に、吸気中のＨ２Ｏの濃度（以下、「吸気Ｈ２Ｏ濃度」と称呼する。）[Ｈ２Ｏ]inの取得方法について説明する。吸気Ｈ２Ｏ濃度[Ｈ２Ｏ]inは、上述した吸気ＣＯ２濃度[ＣＯ２]inの取得方法と全く同様にして、上記(9)〜(12)式における[ＣＯ２]を[Ｈ２Ｏ]に書き換え、且つ、ＫＣＯ２をＫＨ２Ｏに書き換えることで求めることができる。即ち、吸気Ｈ２Ｏ濃度[Ｈ２Ｏ]inは、下記(13)式に従って求めることができる。
【００７７】
【数６】
2005180220.tif 000008
【００７８】
上記(13)式において、ＫＨ２Ｏは、単位質量あたりの燃料の燃焼により発生するＨ２Ｏ質量である。上記(4)式から理解できるように燃料(ＣｍＨｎ)１molの燃焼により(ｎ/２)
molのＨ２Ｏが発生するから、燃料(ＣｍＨｎ)の分子量をＭfuel、Ｈ２Ｏの分子量をＭＨ２Ｏ（＝１８）とすると、質量Ｍfuelの燃料の燃焼により質量｛(ｎ/２)・ＭＨ２Ｏ｝のＨ２Ｏが発生することになる。従って、ＫＨ２Ｏは、
(ｎ/２)・(ＭＨ２Ｏ/Ｍfuel) として求めることができる。
【００７９】
次に、吸気中のＮ２の濃度（以下、「吸気Ｎ２濃度」と称呼する。）[Ｎ２]inの取得方法について説明する。Ｎ２は、燃料(ＣｍＨｎ)の燃焼により筒内で発生することはなく、また、筒内で消費されることもない。従って、吸気Ｎ２濃度[Ｎ２]inは、上記(9)〜(12)式における[ＣＯ２]を[Ｎ２]に書き換え、且つ、ＫＣＯ２の項を削除することで求めることができる。即ち、吸気Ｎ２濃度[Ｎ２]inは、下記(14)式に従って求めることができる。
【００８０】
【数７】
2005180220.tif 000009
【００８１】
次に、吸気中のＯ２の濃度（以下、「吸気Ｏ２濃度」と称呼する。）[Ｏ２]inの取得方法について説明する。筒内ガス中のＯ２は、燃料(ＣｍＨｎ)の燃焼により筒内で消費される。燃焼により消費されるＯ２質量は、単位質量あたりの燃料の燃焼により消費されるＯ２質量をＫＯ２とすると、ＫＯ２・qfinc
と表すことができる。従って、吸気Ｏ２濃度[Ｏ２]inは、上記(9)〜(12)式における[ＣＯ２]を[Ｏ２]に書き換え、且つ、ＫＣＯ２を「−ＫＯ２」に書き換えることで求めることができる。即ち、吸気Ｏ２濃度[Ｏ２]inは、下記(15)式に従って表すことができる。
【００８２】
【数８】
2005180220.tif 000010
【００８３】
ここで、上記(4)式から理解できるように燃料(ＣｍＨｎ)１molの燃焼により(ｍ＋ｎ/４) molのＯ２が消費されるから、燃料(ＣｍＨｎ)の分子量をＭfuel、Ｏ２の分子量をＭＯ２（＝３２）とすると、質量Ｍfuelの燃料の燃焼により質量｛(ｍ＋ｎ/４)・ＭＯ２｝のＯ２が消費されることになる。従って、上記(15)式におけるＫＯ２は、(ｍ＋ｎ/４)・(ＭＯ２/Ｍfuel)
として求めることができる。
【００８４】
なお、理論空燃比stoich＝ＫＯ２/[Ｏ２]airの関係、及び、空気過剰率λ＝Gm/(stoich・qfinc)の関係を利用して、上記(15)式からＫＯ２とGmを消去すると、下記(16)式が得られ、更には、下記(16)式において「(1/stoich)≒０」の関係を考慮すると、吸気Ｏ２濃度[Ｏ２]inは、下記(17)式に従って求めることができる。
【００８５】
【数９】
2005180220.tif 000011
【００８６】
以上のようにして、上記４つの成分の吸気中の質量濃度[Ｘ]in（Ｘ：Ｏ２，ＣＯ２，Ｈ２Ｏ，Ｎ２）をそれぞれ求めることができる。次に、この結果を利用して、上記モル数α，β，γを求める方法について説明する。上記４つの成分の吸気中の質量濃度が求められると、同４つの成分の吸気中の質量濃度の比率
[Ｏ２]in：[ＣＯ２]in：[Ｈ２Ｏ]in：[Ｎ２]in を得ることができる。そして、この４つの成分の吸気中の質量濃度の比率は、筒内に吸入された筒内ガス（より具体的には、上記燃料領域（Ｂ領域）に存在するガス）中の同４つの成分の質量濃度の比率として維持されると仮定する。
【００８７】
一方、筒内ガス中の上記４つの成分の質量濃度の比率は、同筒内ガス中の同４つの成分の質量の比率と等しい。また、上記(4)式から理解できるように、同筒内ガス中の同４つの成分のモル数の比率は (ｍ＋ｎ/４)：α：β：γ であるから、同筒内ガス中の同４つの成分の質量の比率は 
(ｍ＋ｎ/４)ＭＯ２：αＭＣＯ２：βＭＨ２Ｏ：γＭＮ２ と表すことができる。以上のことから、下記(18)式が成立する。
【００８８】
【数１０】
2005180220.tif 000012
【００８９】
上記(18)式を利用すれば、上記モル数α，β，γをそれぞれ、下記(19)〜(21)式に従って求めることができる。
【００９０】
【数１１】
2005180220.tif 000013
【００９１】
このようにして、上記モル数α，β，γが求められると、上記(5)式、及び上記(6)式に従って上記モル数ngas、及び上記定圧比熱Cpを求めることができ、この結果、上記(3)式に従って燃焼起因温度上昇量ΔTburnを求めることができる。
【００９２】
以上、本装置は、上記(3)〜(21)式を利用して、少なくとも吸気の組成（吸気中のガス成分の濃度比率）と、噴射された燃料の燃焼による発熱量（Qfuel）とに基づいて、上記最終燃料噴射時期finjfinが到来する毎に、燃焼起因温度上昇量ΔTburnを求める。
【００９３】
<燃焼速度起因温度上昇量ΔTb_veloの取得>
筒内ガス温度Taが、着火時における温度（即ち、着火時圧縮筒内ガス温度Tpump）から上記燃焼起因温度上昇量ΔTburnだけ上昇して温度（Tpump＋ΔTburn）となるまで（即ち、着火時から筒内ガス温度Taが最高燃焼温度Tflameに到達するまで）には所定の時間（以下、「最高温度到達時間」と称呼する。）が必要である。他方、着火時が圧縮上死点に対応する時期よりも後の時点（即ち、膨張行程に対応する時点）である場合、時間経過に伴なって筒内容積Vaの増大に起因して筒内ガス温度Taが低下していく（燃焼による筒内ガス温度Taの上昇が抑制される）。
【００９４】
従って、上記最高温度到達時間が短いほど筒内ガスの最高燃焼温度Tflameが高くなる。換言すれば、燃焼開始後における筒内の燃焼速度が大きいほど最高燃焼温度Tflameが高くなる。かかる筒内ガスの温度上昇量が燃焼速度起因温度上昇量ΔTb_veloに相当する。以下、この燃焼速度起因温度上昇量ΔTb_veloの求め方について図４を参照しながら説明する。
【００９５】
筒内の燃焼速度に影響を与える因子としては種々のものが考えられるが、本例では、（今回の）燃料噴射圧力Pcrc(＝前記指令基本燃料噴射圧力Pcrbase)、及びエンジン回転速度NE（具体的には、下死点時エンジン回転速度NE0）を選択する。今回の燃料噴射圧力Pcrc、或いは下死点時エンジン回転速度NE0が増大すると、燃焼速度も増大する。
【００９６】
ここで、今回の燃料噴射圧力Pcrc、及び下死点時エンジン回転速度NE0がそれぞれ、基準燃料噴射圧力Pcrref、及び基準エンジン回転速度NErefである場合における燃焼速度を通常燃焼速度と定義する。そして、筒内の燃焼速度が上記通常燃焼速度となるとき、筒内ガスの最高燃焼温度Tflameは、温度（Tpump＋ΔTburn）（以下、「通常燃焼速度時最高燃焼温度Tig」と称呼する。）と等しくなるもの（従って、燃焼速度起因温度上昇量ΔTb_velo＝０となるもの）とする。
【００９７】
そうすると、上記通常燃焼速度に相当する最高温度到達時間からの最高温度到達時間短縮量Δtadvは、今回の燃料噴射圧力Pcrc、及び下死点時エンジン回転速度NE0の関数として下記(22)式により表すことができる。
【００９８】
Δtadv＝t0＋Ka・Pcrc＋Kb・NE0  ・・・(22)
【００９９】
上記(22)式において、t0は負の定数、Ka,Kbは正の定数である。t0,Ka,Kbは、今回の燃料噴射圧力Pcrc、及び下死点時エンジン回転速度NE0がそれぞれ、基準燃料噴射圧力Pcrref、及び基準エンジン回転速度NErefとなる場合にΔtadv＝０となるように設定される値である。また、上記(22)式により求められる最高温度到達時間短縮量Δtadvに相当するクランク角度CAの進角量であるΔtadv相当進角量ΔCAadvは、同最高温度到達時間短縮量Δtadvと、下死点時エンジン回転速度NE0と、ΔtadvとNE0とを引数とする関数ｈと、に基づいてh(tadv,NE0)として求めることができる。
【０１００】
そして、図４に示したように、上記着火時クランク角度CAig（従って、着火時筒内容積Vig）において上記通常燃焼速度時最高燃焼温度Tigとなっている筒内ガスが、断熱変化して、同着火時クランク角CAigから上記Δtadv相当進角量ΔCAadvだけ進角したクランク角度である補正着火時クランク角度CAadv（従って、補正着火時筒内容積Vadv(＝Va(CAadv))）まで断熱変化により圧縮されて温度Tadvになったものと仮定する。
【０１０１】
そうすると、この温度Tadvから上記通常燃焼速度時最高燃焼温度Tigを減じた値は、燃焼速度が上記通常燃焼速度となる場合に比して、筒内ガス温度が最高燃焼温度Tflameに到達する時期が、上記最高温度到達時間短縮量Δtadvに相当する分だけ早くなることに起因して発生する筒内ガス温度の上昇量（従って、上記燃焼速度起因温度上昇量ΔTb_velo）に相当すると考えることができる。従って、燃焼速度起因温度上昇量ΔTb_veloは、下記(23)式に従って求められる。
【０１０２】
ΔTb_velo＝Tig・(Vig/Vadv)κ-1−Tig  ・・・(23)
【０１０３】
以上、本装置は、上記(22)式、及び(23)式等を利用して、筒内の燃焼速度に影響を与える因子となる今回の燃料噴射圧力Pcrc、及び下死点時エンジン回転速度NE0に基づいて、上記最終燃料噴射時期finjfinが到来する毎に、燃焼起因温度上昇量ΔTb_veloを求める。
【０１０４】
そして、本装置は、上記最終燃料噴射時期finjfinが到来する毎に、上記(1)式に従って、上記着火時圧縮筒内ガス温度Tpumpに、上記燃焼起因温度上昇量ΔTburnと上記燃焼速度起因温度上昇量ΔTb_veloとを加えた値を筒内ガスの最高燃焼温度Tflameとして推定する。以上が、本発明による燃焼温度推定方法の概要である。
【０１０５】
（ＮＯｘ発生量推定方法の概要）
次に、本装置が実施する、上記推定された最高燃焼温度Tflameを使用したＮＯｘ発生量推定方法の概要について説明していく。このＮＯｘ発生量推定方法では、燃料噴射気筒について上記最終燃料噴射時期finjfinが到来する毎に、その直後の爆発（膨張）行程において上記Ｂ領域において燃焼により発生する実ＮＯｘ発生量NOxactが推定されていく。
【０１０６】
実ＮＯｘ量発生量NOxactは、単位燃料量あたりの燃焼発生ＮＯｘ量（以下、「燃焼発生ＮＯｘ率RNOx_burn」と称呼する。）に上記今回の指令燃料噴射量qfincを乗じた値として、下記(24)式に従って求めることができる。
【０１０７】
NOxact＝RNOx_burn・qfinc  ・・・(24)
【０１０８】
ここで、上記(24)式における燃焼発生ＮＯｘ率RNOx_burnは下記(25)式にて推定される。下記(25)式において、eは自然対数の底である。上述したように、RO2cは下死点時吸気酸素濃度であり、qfincは今回の指令燃料噴射量であり、Pcrcは今回の指令燃料噴射圧力(=Pcrbase)であり、Tflameは上記(1)式にて求められる今回の爆発行程における最高燃焼温度である。また、K0〜K4は代表的な公知の所定の重回帰分析に基づいて後述するように決定された適合定数である。
【０１０９】
RNOx_burn＝eK0・(RO2c)K1・(qfinc)K2・(Pcrc)K3・e(K4/Tflame)  ・・・(25)
【０１１０】
即ち、上記(25)式は、燃焼発生ＮＯｘ率RNOx_burnを求めるための実験式である。上記(25)式にて推定される燃焼発生ＮＯｘ率RNOx_burnは、上記下死点時吸気酸素濃度RO2c、上記今回の指令燃料噴射量qfinc、上記今回の指令燃料噴射圧力Pcrc、及び上記最高燃焼温度Tflameの関数であって、より具体的には、下死点時吸気酸素濃度RO2cのべき乗、今回の指令燃料噴射量qfincのべき乗、今回の指令燃料噴射圧力Pcrcのべき乗、及び、最高燃焼温度Tflameに応じて決定される値を指数とする指数関数、の積に基づいて計算される。
【０１１１】
適合定数K0〜K4は、例えば、以下のような実験を行うことにより決定され得る。即ち、先ず、ＥＧＲ制御弁５２を閉状態に維持して機関１０を運転させることで、排気弁Voutから排出された排ガス（従って、同排ガス中のＮＯｘ）の総てが排気通路から外部へ排出されるようにしておく。これにより、排気通路から外部へ排出される排ガス中のＮＯｘ量は上記実ＮＯｘ発生量NOxactと等しくなるから、ＮＯｘ排出量を所定のＮＯｘ濃度センサの出力に基づいて計測することで実ＮＯｘ発生量NOxact（従って、燃焼発生ＮＯｘ率RNOx_burn（＝NOxact/qfinc））を計測できるようになる。
【０１１２】
次に、この状態にて、上記下死点時吸気酸素濃度RO2c、上記今回の指令燃料噴射量qfinc、上記今回の指令燃料噴射圧力Pcrc、及び上記(1)式にて求められる最高燃焼温度Tflameの各値の組み合わせを種々の所定パターンになるように順次変更していき、そのパターン毎に、燃焼発生ＮＯｘ率RNOx_burnを順次計測していく。
【０１１３】
そして、このような作業（実験）の結果得られた上記各値の組み合わせと、計測された燃焼発生ＮＯｘ率RNOx_burnの値との関係に関する多数のデータに基づいて、上記公知の所定の重回帰分析を行って上記適合定数K0〜K4を求めることができる。ここにおいて、少なくとも、適合定数K1〜K3は正の値に決定され、適合定数K4は負の値に決定される。
【０１１４】
従って、上記(25)式から理解できるように、(25)式に従って推定・計算される燃焼発生ＮＯｘ率RNOx_burn（従って、実ＮＯｘ発生量NOxact）は、下死点時吸気酸素濃度RO2c、今回の指令燃料噴射量qfinc、今回の指令燃料噴射圧力Pcrc、及び最高火炎温度Tflameのうちのどの値が増加しても増大する。このことは、以下の実際の現象に沿うものである。
【０１１５】
即ち、先ず、吸気酸素濃度RO2_inが増加すると実ＮＯｘ発生量NOxactが増大する。これは、酸素はＮＯｘを発生させる材料であって燃焼室内の酸素量が多くなると当然にＮＯｘが発生し易くなることに基づく現象である。
【０１１６】
また、燃料噴射量qfinが増加すると実ＮＯｘ発生量NOxactが増大する。これは、燃料噴射量qfinの増加に応じて機関１０に与えられる負荷の大きさが増大することから、燃焼室の内壁面温度の上昇に起因して、燃料噴射量qfinが増大するほど（即ち、機関に与えられる負荷が増加するほど）ＮＯｘが発生し易くなることに基づく現象である。
【０１１７】
また、燃料噴射圧力Pcrが増加すると実ＮＯｘ発生量NOxactが増大する。これは、燃料噴射圧力Pcrの増加に応じて燃料の噴射速度が増大すること等により同燃料の霧化の程度が大きくなるから、上述した空気過剰率の増大に起因して、燃料噴射圧力Pcrが増大するほど（即ち、噴射された燃料の霧化の程度が大きくなるほど）ＮＯｘが発生し易くなることに基づく現象である。
【０１１８】
また、最高燃焼温度Tflameが増加すると実ＮＯｘ発生量NOxactが増大する。これは、ガス温度が高くなるほど窒素からＮＯｘが生成される化学反応が促進されることに基づく現象である。以上のことから、上記(25)式に従って燃焼発生ＮＯｘ率RNOx_burnを計算すれば、少なくとも上記４つの実際の現象に沿うように同燃焼発生ＮＯｘ率RNOx_burn（従って、上記(24)式に従う実ＮＯｘ発生量NOxact）を精度良く推定することができる。以上が、ＮＯｘ発生量推定方法の概要である。
【０１１９】
（燃料噴射制御の概要）
上記ＮＯｘ発生量推定方法を実施する本装置は、一作動サイクルあたりの目標ＮＯｘ発生量NOxtを、上記燃料噴射量qfinとエンジン回転速度NEとに基づいて逐次計算する。そして、本装置は、前回推定された実ＮＯｘ発生量NOxactが目標ＮＯｘ発生量NOxtに一致するように、最終燃料噴射開始時期finjfin、及びＥＧＲ制御弁５２の開度をフィードバック制御する。
【０１２０】
具体的には、前回推定された上記実ＮＯｘ発生量NOxactの値が上記目標ＮＯｘ発生量NOxtよりも大きいとき、今回の燃料噴射気筒についての最終燃料噴射開始時期finjfinを基本燃料噴射時期finjbaseよりも所定量だけ遅らせ、且つ、ＥＧＲ制御弁５２の開度を現時点での値から所定開度だけ大きくする。これにより、今回の燃料噴射気筒についての最高燃焼温度Tflameが低くなるように制御され、この結果、今回の燃料噴射気筒内で発生する実ＮＯｘ発生量NOxactが上記目標ＮＯｘ発生量NOxtに一致せしめられる。
【０１２１】
一方、前回推定された上記実ＮＯｘ発生量NOxactの値が上記目標ＮＯｘ発生量NOxtよりも小さいとき、今回の燃料噴射気筒についての最終燃料噴射開始時期finjfinを基本燃料噴射時期finjbaseよりも所定量だけ早め、且つ、ＥＧＲ制御弁５２の開度を現時点での値から所定開度だけ小さくする。これにより、今回の燃料噴射気筒についての最高燃焼温度Tflameが高くなるように制御され、この結果、今回の燃料噴射気筒内で発生する実ＮＯｘ発生量NOxactが上記目標ＮＯｘ発生量NOxtに一致せしめられる。以上が、燃料噴射制御の概要である。
【０１２２】
（燃焼発生ＮＯｘ率RNOx_burnの実際の演算方法）
上記(25)式に従って燃焼発生ＮＯｘ率RNOx_burnを演算するためには、「べき乗」の演算、及び「掛け算」の演算が必要となる。しかしながら、一般に、マイクロコンピュータを使用して「べき乗」の演算を実行すると計算負荷が増大し、同マイクロコンピュータを使用して「掛け算」の演算を実行すると計算精度が低下する傾向がある。よって、「べき乗」の演算、及び「掛け算」の演算を避けるため、本装置（ＣＰＵ６１）は、実際には、上記(25)式において両辺の自然対数をとることで得られる下記(26)式を利用して、テーブル検索と「足し算」の演算のみで燃焼発生ＮＯｘ率RNOx_burnを計算する。
【０１２３】
log(RNOx_burn)＝K0＋K1・log(RO2c)＋K2・log(qfinc)＋K3・log(Pcrc)＋K4/Tflame  ・・・(26)
【０１２４】
即ち、本装置は、上記(26)式の右辺の第２項〜第５項までの各項を求めるために予めＲＯＭ６２にそれぞれ記憶されているテーブルMaplog1(RO2c)、Maplog2(qfinc)、Maplog3(Pcrc)、及びMapinvpro(Tflame)に基づいて、テーブル検索値dataMap1(＝K1・log(RO2c))、dataMap2(＝K2・log(qfinc))、dataMap3(＝K3・log(Pcrc))、及びdataMap4(＝K4/Tflame)をそれぞれ決定し、「足し算」を伴う下記(27)式に従って「log(RNOx_burn)」を求める。
【０１２５】
log(RNOx_burn)＝K0＋dataMap1＋dataMap2＋dataMap3＋dataMap4  ・・・(27)
【０１２６】
そして、本装置は、上記（27）式に従って求められた「log(RNOx_burn)」から燃焼発生ＮＯｘ率RNOx_burnを求めるために予めＲＯＭ６２に記憶されているテーブルMapinvlog(log(RNOx_burn))に基づいて同燃料発生ＮＯｘ率RNOx_burnを求める。これにより、ＣＰＵ６１の計算負荷を軽減できるとともに計算精度の低下を防止できる。
【０１２７】
（実際の作動）
次に、上記のように構成された内燃機関の制御装置の実際の作動について説明する。
<燃料噴射量等の制御>
ＣＰＵ６１は、図５にフローチャートにより示した燃料噴射量等の制御を行うためのルーチンを所定時間の経過毎に繰り返し実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ５００から処理を開始し、ステップ５０５に進んでアクセル開度Accp、エンジン回転速度NE、及び図６に示したテーブル（マップ）Mapqfinから（指令）燃料噴射量qfinを求める。テーブルMapqfinは、アクセル開度Accp及びエンジン回転速度NEと燃料噴射量qfinとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【０１２８】
次いで、ＣＰＵ６１はステップ５１０に進み、燃料噴射量qfin、エンジン回転速度NE、及び図７に示したテーブルMapfinjbaseから基本燃料噴射時期finjbaseを決定する。テーブルMapfinjbaseは、燃料噴射量qfin及びエンジン回転速度NEと基本燃料噴射時期finjbaseとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【０１２９】
その後、ＣＰＵ６１はステップ５１５に進んで、燃料噴射量qfin、エンジン回転速度NE、及び図８に示したテーブルMapPcrbaseから基本燃料噴射圧力Pcrbaseを決定する。テーブルMapPcrbaseは、燃料噴射量qfin及びエンジン回転速度NEと基本燃料噴射圧力Pcrbaseとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【０１３０】
次に、ＣＰＵ６１はステップ５２０に進み、燃料噴射量qfin、エンジン回転速度NE、及び図９に示したテーブルMapNOxtから目標ＮＯｘ発生量NOxtを決定する。テーブルMapNOxtは、燃料噴射量qfin及びエンジン回転速度NEと目標ＮＯｘ排出量NOxtとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【０１３１】
次いで、ＣＰＵ６１はステップ５２５に進んで、前記決定した目標ＮＯｘ発生量NOxtから後述するルーチンにより求められている最新の（具体的には、前回の燃料噴射時期に演算された）実ＮＯｘ発生量NOxactを減じた値をＮＯｘ発生量偏差ΔNOxとして格納する。
【０１３２】
続いて、ＣＰＵ６１はステップ５３０に進んで、噴射時期補正値Δθを、前記ＮＯｘ発生量偏差ΔNOxと図１０に示したテーブルMapΔθとから決定する。テーブルMapΔθは、ＮＯｘ発生量偏差ΔNOxと噴射時期補正値Δθとの関係を規定するテーブルであり、ＲＯＭ６２内に格納されている。
【０１３３】
次いで、ＣＰＵ６１はステップ５３５に進み、基本噴射時期finjbaseを噴射時期補正値Δθで補正して最終燃料噴射時期finjfinを決定する。これにより、ＮＯｘ発生量偏差ΔNOxに応じて噴射時期が補正されることになる。この場合、図１０から明らかなように、ＮＯｘ発生量偏差ΔNOxが正の大きい値になるほど噴射時期補正値Δθが正の大きな値となって最終燃料噴射時期finjfinが進角側となり、同ＮＯｘ発生量偏差ΔNOxが負の大きい値（絶対値が大きい値）になるほど噴射時期補正値Δθは負の大きな値となって最終燃料噴射時期finjfinが遅角側に移行される。
【０１３４】
続いて、ＣＰＵ６１はステップ５４０に進み、燃料噴射気筒についての噴射開始時期（即ち、上記決定された最終燃料噴射時期finjfin）が到来したか否かを判定し、「Ｎｏ」と判定する場合はステップ５９５に直ちに進んで本ルーチンを一旦終了する。
【０１３５】
一方、ステップ５４０の判定において「Ｙｅｓ」と判定する場合、ＣＰＵ６１はステップ５４５に進んで上記決定された（指令）燃料噴射量qfinの燃料を燃料噴射気筒についての燃料噴射弁２１から上記決定された基本燃料噴射圧力Pcrbaseをもって噴射するとともに、続くステップ５５０にて上記ＮＯｘ発生量偏差ΔNOxが正の値であるか否かを判定し、「Ｙｅｓ」と判定する場合、ステップ５５５に進んでＥＧＲ制御弁５２の開度を現時点での値よりも所定開度だけ小さくした後、ステップ５７０に進む。
【０１３６】
ステップ５５０の判定において「Ｎｏ」と判定する場合、ＣＰＵ６１はステップ５６０に進み、上記ＮＯｘ発生量偏差ΔNOxが負の値であるか否かを判定する。ステップ５６０の判定において、ＣＰＵ６１は「Ｙｅｓ」と判定する場合、ＥＧＲ制御弁５２の開度を現時点での値よりも所定開度だけ大きくした後にステップ５７０に進む一方、「Ｎｏ」と判定する場合（即ち、ＮＯｘ発生量偏差ΔNOxの値が「０」のとき）、ＥＧＲ制御弁５２の開度を変更することなくステップ５７０に進む。
【０１３７】
このようにして、ＮＯｘ発生量偏差ΔNOxに応じてＥＧＲ制御弁５２の開度が変更される。そして、ステップ５７０に進むと、ＣＰＵ６１は実際に噴射した上記燃料噴射量qfinの値を今回の燃料噴射量qfincとして格納し、続くステップ５７５にて実際の噴射圧力である上記基本燃料噴射圧力Pcrbaseの値を今回の燃料噴射圧力Pcrcとして格納した後、ステップ５９５に進んで本ルーチンを一旦終了する。以上により、燃料噴射量、燃料噴射時期、燃料噴射圧力、及びＥＧＲ制御弁５２の開度の制御が達成される。
【０１３８】
<最高燃焼温度の計算>
また、ＣＰＵ６１は、図１１にフローチャートにより示した最高燃焼温度Tflameの計算を行うためのルーチンを所定時間の経過毎に繰り返し実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ１１００から処理を開始し、ステップ１１０５に進んで現時点がATDC-180°に一致しているか否かを判定する。
【０１３９】
いま、現時点がATDC-180°より前であるものとして説明を続けると、ＣＰＵ６１はステップ１１０５にて「Ｎｏ」と判定してステップ１１４５に直ちに進み、燃料噴射気筒についての燃料噴射開始時期（即ち、上記最終燃料噴射時期finjfin）が到来したか否かを判定する。現時点はATDC-180°よりも前であるから、ＣＰＵ６１はステッ１１４５にて「Ｎｏ」と判定し、ステッ１１９５に直ちに進んで本ルーチンを一旦終了する。
【０１４０】
以降、ＣＰＵ６１はATDC-180°が到来するまでの間、ステップ１１００、１１０５、１１４５、１１９５の処理を繰り返し実行する。そして、ATDC-180°が到来すると、ＣＰＵ６１はステップ１１０５に進んだとき「Ｙｅｓ」と判定してステップ１１１０に進むようになり、ステップ１１１０にて、吸気温センサ７２、吸気管圧力センサ７３、エアフローメータ７１、及びクランクポジションセンサ７４によりそれぞれ検出される現時点（即ち、ATDC-180°）での吸気温度Tb、吸気管圧力Pb、吸入新気流量Ga、及びエンジン回転速度NEを、それぞれ、下死点時筒内ガス温度Ta0、下死点時筒内ガス圧力Pa0、下死点時吸入新気流量Ga0、及び下死点時エンジン回転速度NE0として格納する。
【０１４１】
次いで、ＣＰＵ６１はステップ１１１５に進み、吸気酸素濃度センサ７６により検出される現時点（即ち、ATDC-180°）での吸気酸素濃度RO2_inを下死点時吸気酸素濃度RO2cとして格納し、続くステップ１１２０にて、上記(7)式に従って筒内総ガス量Gcylを求める。ここで、下死点時筒内ガス圧力Pa0、及び下死点時筒内ガス温度Ta0としては、ステップ１１１０にて格納されている値が用いられる。
【０１４２】
続いて、ＣＰＵ６１はステップ１１２５に進んで、上記下死点時吸入新気流量Ga0と、上記下死点時エンジン回転速度NE0と、上記関数ｆとに基づいて吸入新気量Gmを求め、続くステップ１１３０にて、ステップ１１２０にて求めた筒内総ガス量Gcylと、前記吸入新気量Gmと、上記(8)式とに基づいてＥＧＲガス量Gegrを求める。次に、ＣＰＵ６１はステップ１１３５に進んで、上記吸入新気量Gmと、上記ＥＧＲガス量Gegrと、ステップ１１３５内に記載の式とに基づいてＥＧＲ率Ｒを求め、続くステップ１１４０にて、上記吸入新気量Gmと、先のステップ５７０にて格納されている今回の燃料噴射量qfincと、ステップ１１４０内に記載の式とに基づいて空気過剰率λを求める。そして、ＣＰＵ６１はステップ１１４５に進んで「Ｎｏ」と判定してステップ１１９５に直ちに進んで本ルーチンを一旦終了する。
【０１４３】
以降、ＣＰＵ６１は燃料噴射時期（即ち、上記最終燃料噴射時期finjfin）が到来するまでの間、ステップ１１００、１１０５、１１４５、１１９５の処理を再び繰り返し実行する。そして、上記最終燃料噴射時期finjfinが到来すると、ＣＰＵ６１はステップ１１４５にて「Ｙｅｓ」と判定してステップ１１５０に進むようになり、ステップ１１５０にて上記最終燃料噴射時期finjfinと上記所定の着火遅れ時間とから着火時クランク角度CAigを求める。
【０１４４】
次いで、ＣＰＵ６１はステップ１１５５を経由して、図１２にフローチャートにより示した着火時圧縮筒内ガス温度Tpumpの計算を行うためのルーチンを実行する。即ち、ＣＰＵ６１はステップ１２００から処理を開始し、ステップ１２０５に進んで、先のステップ１１１５にて求めた最新の下死点時吸気酸素濃度RO2cと、上記関数ｇとに基づいてポリトロープ指数κを求める。
【０１４５】
次に、ＣＰＵ６１はステップ１２１０に進み、先のステップ１１５０にて求めた最新の着火時クランク角度CAigがATDC0°よりも遅れた角度となっているか否かを判定し、「Ｙｅｓ」と判定する場合、ステップ１２１５に進んで同着火時クランク角度CAigに対応する筒内容積をそのまま着火時筒内容積Vigとして格納する。一方、ステップ１２１０にて「Ｎｏ」と判定する場合には、ＣＰＵ６１はステップ１２２０に進んで上死点（ATDC0°）に対応する筒内容積Vtopを着火時筒内容積Vigとして格納する。
【０１４６】
続いて、ＣＰＵ６１はステップ１２２５に進み、グロープラグ２４への通電がなされているか否かを判定し、「Ｙｅｓ」と判定する場合、ステップ１２３０に進んで前記所定値Tglowを筒内ガス温度上昇量ΔTpumpとして格納し、「Ｎｏ」と判定する場合、ステップ１２３５にて同筒内ガス温度上昇量ΔTpumpの値を「０」に設定する。
【０１４７】
そして、ＣＰＵ６１はステップ１２４０に進んで、先のステップ１１１０にて求めた最新の下死点時筒内ガス温度Ta0と、上記着火時筒内容積Vigと、上記筒内ガス温度上昇量ΔTpumpと、ステップ１２４０内に記載の式とに基づいて着火時圧縮筒内ガス温度Tpumpを求め、ステップ１２９５を経由して図１１のステップ１１６０に進む。
【０１４８】
ＣＰＵ６１はステップ１１６０に進むと、図１３にフローチャートにより示した燃焼起因温度上昇量ΔTburnの計算を行うためのルーチンを実行する。即ち、ＣＰＵ６１はステップ１３００から処理を開始し、ステップ１３０５に進んで、先のステップ１１３５にて求めた最新のＥＲＧ率Ｒと、ステップ１１４０にて求めた最新の空気過剰率λと、上記(17)式に相当するステップ１３０５内に記載の式とに基づいて吸気Ｏ２濃度[Ｏ２]inを求める。
【０１４９】
次いで、ＣＰＵ６１はステップ１３１０に進み、上記ＥＧＲ率Ｒと、上記今回の燃料噴射量qfincと、先のステップ１１２５にて求めた最新の吸入新気量Gmと、上記(12)式に相当するステップ１３１０内に記載の式に基づいて吸気ＣＯ２濃度[ＣＯ２]inを求める。同様に、ＣＰＵ６１はステップ１３１５に進んで、上記(13)式に従って吸気Ｈ２Ｏ濃度[Ｈ２Ｏ]inを求め、続くステップ１３２０にて、上記(14)式に従って吸気Ｎ２濃度[Ｎ２]inを求める。
【０１５０】
続いて、ＣＰＵ６１はステップ１３２５に進んで、上記吸気ＣＯ２濃度[ＣＯ２]inと、上記吸気Ｏ２濃度[Ｏ２]inと、上記(19)式とに基づいて上記モル数αを求め、同様に、続くステップ１３３０にて上記(20)式に従って上記モル数βを求めるとともに、続くステップ１３３５にて上記(21)式に従って上記モル数γを求める。
【０１５１】
次に、ＣＰＵ６１はステップ１３４０に進んで、上記求めたモル数α，β，γと、上記(5)式とに基づいて燃焼後筒内ガスモル数ngasを求め、続くステップ１３４５にて同モル数α，β，γと、上記(6)式とに基づいて燃焼後筒内ガス定圧比熱Cpを求める。そして、ＣＰＵ６１はステップ１３５０に進んで、上記燃焼後筒内ガスモル数ngasと、上記燃焼後筒内ガス定圧比熱Cpと、上記(3)式とに基づいて燃焼起因温度上昇量ΔTburnを求め、ステップ１３９５を経由して図１１のステップ１１６５に進む。
【０１５２】
ＣＰＵ６１はステップ１１６５に進むと、図１４にフローチャートにより示した燃焼速度起因温度上昇量ΔTb_veloの計算を行うためのルーチンを実行する。即ち、ＣＰＵ６１はステップ１４００から処理を開始し、ステップ１４０５に進んで、先のステップ１２４０にて求めた最新の着火時圧縮筒内ガス温度Tpumpに、先のステップ１３５０にて求めた最新の燃焼起因温度上昇量ΔTburnを加えた値を通常燃焼速度時最高燃焼温度Tigとして格納する。
【０１５３】
次に、ＣＰＵ６１はステップ１４１０に進み、先のステップ５７５にて格納した今回の燃料噴射圧力Pcrcと、先のステップ１１１０にて格納した最新の下死点時エンジン回転速度NE0と、上記(22)式とに基づいて最高温度到達時間短縮量Δtadvを求め、続くステップ１４１５にて、同最高温度到達時間短縮量Δtadvと、同下死点時エンジン回転速度NE0と、上記関数ｈとに基づいてΔtadv相当進角量ΔCAadvを求める。
【０１５４】
続いて、ＣＰＵ６１はステップ１４２０に進み、先のステップ１１５０にて求めた最新の着火時クランク角度CAigから上記Δtadv相当進角量ΔCAadvだけ進ませた角度を補正着火時クランク角度CAadvとして設定し、続くステップ１４２５にて上記補正着火時クランク角度CAadvがATDC0°よりも遅れた角度になっているか否かを判定する。
【０１５５】
ＣＰＵ６１はステップ１４２５の判定において、「Ｙｅｓ」と判定する場合、ステップ１４３０に進んで同補正着火時クランク角度CAadvに対応する筒内容積をそのまま補正着火時筒内容積Vadvとして格納する。一方、ステップ１４２５にて「Ｎｏ」と判定する場合、ＣＰＵ６１はステップ１４３５に進んで上死点（ATDC0°）に対応する筒内容積Vtopを補正着火時筒内容積Vadvとして格納する。
【０１５６】
次いで、ＣＰＵ６１はステップ１４４０に進んで、先のステップ１２１５、或いはステップ１２２０にて求めた着火時筒内容積Vigと、上記補正着火時筒内容積Vadvと、上記通常燃焼速度時最高燃焼温度Tigと、上記(23)式とに基づいて燃焼速度起因温度上昇量ΔTb_veloを求め、ステップ１４９５を経由して図１１のステップ１１７０に進む。そして、ＣＰＵ６１はステップ１１７０に進むと、先のステップ１２４０にて求めた最新の着火時圧縮筒内ガス温度Tpumpと、先のステップ１３５０にて求めた最新の燃焼起因温度上昇量ΔTburnと、先のステップ１４４０にて求めた最新の燃焼速度起因温度上昇量ΔTb_veloと、上記(1)式とに基づいて筒内ガスの最高燃焼温度Tflameを求めた後、ステップ１１９５に進んで本ルーチンを一旦終了する。
【０１５７】
以降、ＣＰＵ６１は、次の燃料噴射気筒についてのATDC-180°が到来するまでの間、ステップ１１００、１１０５、１１４５、１１９５の処理を繰り返し実行する。以上のようにして、新たな筒内ガスの最高燃焼温度Tflameが燃料噴射開始時期が到来する毎に求められていく。
【０１５８】
<ＮＯｘ発生量の計算>
また、ＣＰＵ６１は、図１５にフローチャートにより示した実ＮＯｘ発生量NOxactの計算を行うためのルーチンを所定時間の経過毎に繰り返し実行するようになっている。従って、所定のタイミングになると、ＣＰＵ６１はステップ１５００から処理を開始し、ステップ１５０５に進んで燃料噴射開始時期（即ち、上記最終燃料噴射時期finjfin）が到来したか否かを判定し、「Ｎｏ」と判定する場合、ステップ１５９５に直ちに進んで本ルーチンを一旦終了する。
【０１５９】
いま、燃料噴射開始時期が到来したものとすると、ＣＰＵ６１はステップ１５１０に進み、先のステップ１１１５にて求められている最新の下死点時吸気酸素濃度RO2cと前記テーブルMaplog1とに基づいて前記テーブル検索値dataMap1(＝K1・log(RO2c))を求める。
【０１６０】
同様に、ＣＰＵ６１はステップ１５１５に進んで先のステップ５７０にて格納されている今回の燃料噴射量qfincと前記テーブルMaplog2とに基づいて前記テーブル検索値dataMap2(＝K2・log(qfinc))を求め、続くステップ１５２０にて先のステップ５７５にて格納されている今回の燃料噴射圧力Pcrcと前記テーブルMaplog3とに基づいて前記テーブル検索値dataMap3(＝K3・log(Pcrc))を求め、続くステップ１５２５にて先のステップ１１７０にて求めた最新の最高燃焼温度Tflameと前記テーブルMapinvproとに基づいて前記テーブル検索値dataMap4(＝K4/Tflame)を求める。
【０１６１】
次いで、ＣＰＵ６１はステップ１５３０に進んで、上記(27)式に従って「log(RNOx_burn)」を求め、続くステップ１５３５にて、同log(RNOx_burn)と前記テーブルMapinvlogとに基づいて燃焼発生ＮＯｘ率RNOx_burnを求める。そして、ＣＰＵ６１はステップ１５４０に進んで、上記今回の燃料噴射量qfincと、上記燃焼発生ＮＯｘ率RNOx_burnと、上記(24)式に従って実ＮＯｘ発生量NOxactを求めた後、ステップ１５９５に進んで本ルーチンを一旦終了する。以降、ＣＰＵ６１は次の燃料噴射気筒についての燃料噴射開始時期が到来するまでの間、ステップ１５００、１５０５、１５９５の処理を繰り返し実行する。
【０１６２】
以上のようにして、新たな実ＮＯｘ発生量NOxactが燃料噴射開始時期が到来する毎に求められていく。そして、この新たな実ＮＯｘ発生量NOxactは、前述のごとく、図５のステップ５２５にて使用され、この結果、次の燃料噴射気筒についての最終燃料噴射開始時期finjfin、及びＥＧＲ制御弁５２の開度が同新たな実ＮＯｘ発生量NOxactに基づいてフィードバック制御されていく。
【０１６３】
以上、説明したように、本発明の実施形態に係る内燃機関の燃焼温度推定方法によれば、燃料噴射開始時期が到来する毎に、筒内ガスの状態が断熱変化することを原則的に利用して着火時における燃焼前の筒内ガス温度（着火時圧縮筒内ガス温度Tpump）を推定し、燃料の燃焼による発熱量Qfuelを、吸気中のガスの組成（ガス成分の濃度比率）から求められる燃焼に関与する燃焼後の筒内ガスのモル数ngas及び定圧比熱Cpの積で除することで燃焼による筒内ガスの温度上昇量（燃焼起因温度上昇量ΔTburn）を推定し、筒内の燃焼速度に影響を与える因子である燃料噴射圧力Pcrc及びエンジン回転速度NE0に基づいて同燃焼速度の増大による筒内ガスの温度上昇量（燃焼速度起因温度上昇量ΔTb_velo）を推定する。そして、着火時圧縮筒内ガス温度Tpumpに、燃焼起因温度上昇量ΔTburnと燃焼速度起因温度上昇量ΔTb_veloとを加えた値を筒内における最高燃焼温度Tflameとして推定する。従って、上述のごとく、筒内ガスの最高燃焼温度Tflameを簡易な構成で実際の種々の現象に沿うように精度良く推定することができた。
【０１６４】
本発明は上記実施形態に限定されることはなく、本発明の範囲内において種々の変形例を採用することができる。例えば、上記実施形態においては、筒内の燃焼速度に影響を与える因子として、燃料噴射圧力及びエンジン回転速度を採用しているが、筒内に吸入されるガスのスワール比、過給機による過給圧、筒内に吸入されるガス中の酸素濃度のうちの少なくとも１つを筒内の燃焼速度に影響を与える因子として採用してもよい。
【０１６５】
また、上記実施形態においては、着火時圧縮筒内ガス温度の上昇量（筒内ガス温度上昇量ΔTpump）を発生させる因子としてグローランプへの通電による発熱量が採用されているが、最終燃料噴射時期finjfin（メイン噴射）に先立ってパイロット噴射が行われる場合におけるパイロット噴射された燃料の燃焼による発熱量を同因子として採用してもよい。この場合、パイロット噴射による燃料の噴射量と、同パイロット噴射の時期（メイン噴射時期とパイロット噴射時期との時間間隔（インターバル））等に基づいてパイロット噴射された燃料の燃焼による発熱量（従って、筒内ガスの温度上昇量）を計算するように構成することが好適である。
【図面の簡単な説明】
【０１６６】
【図１】本発明の実施形態に係る内燃機関の燃料温度推定方法を実施する内燃機関の制御装置を４気筒内燃機関（ディーゼル機関）に適用したシステム全体の概略構成図である。
【図２】或る一つの気筒のシリンダ内(筒内)に吸気マニホールドからガスが吸入され、筒内に吸入された筒内ガスが排気マニホールドへ排出される様子を模式的に示した図である。
【図３】圧縮行程、及び膨張行程にて筒内ガスが断熱変化する場合におけるクランク角度と筒内ガス温度との関係を示したグラフである。
【図４】燃焼速度の増大により筒内ガスの最高燃焼温度が上昇することを説明するための図である。
【図５】図１に示したＣＰＵが実行する燃料噴射量等を制御するためのルーチンを示したフローチャートである。
【図６】図１に示したＣＰＵが図５に示したルーチンを実行する際に参照する燃料噴射量を決定するためのテーブルである。
【図７】図１に示したＣＰＵが図５に示したルーチンを実行する際に参照する基本燃料噴射時期を決定するためのテーブルである。
【図８】図１に示したＣＰＵが図５に示したルーチンを実行する際に参照する基本燃料噴射圧力を決定するためのテーブルである。
【図９】図１に示したＣＰＵが図５に示したルーチンを実行する際に参照する目標ＮＯｘ発生量を決定するためのテーブルである。
【図１０】図１に示したＣＰＵが図５に示したルーチンを実行する際に参照する噴射時期補正値を決定するためのテーブルである。
【図１１】図１に示したＣＰＵが実行する燃焼温度の計算を実行するためのルーチンを示したフローチャートである。
【図１２】図１に示したＣＰＵが実行する着火時圧縮筒内ガス温度の計算を実行するためのルーチンを示したフローチャートである。
【図１３】図１に示したＣＰＵが実行する燃焼起因温度上昇量の計算を実行するためのルーチンを示したフローチャートである。
【図１４】図１に示したＣＰＵが実行する燃焼速度起因温度上昇量の計算を実行するためのルーチンを示したフローチャートである。
【図１５】図１に示したＣＰＵが実行するＮＯｘ発生量（実ＮＯｘ発生量）の計算を実行するためのルーチンを示したフローチャートである。
【符号の説明】
【０１６７】
２１…燃料噴射弁、２２…燃料噴射用ポンプ、２４…グロープラグ、３１…吸気マニホールド、３２…吸気管、４１…排気マニホールド、４２…排気管、５０…ＥＧＲ装置、５２…ＥＧＲ制御弁、６０…電気制御装置、６１…ＣＰＵ、７１…エアフローメータ、７２…吸気温センサ、７３…吸気管圧力センサ、７４…クランクポジションセンサ、７６…吸気酸素濃度センサ
【図１】
2005180220.tif 000014
【図２】
2005180220.tif 000015
【図３】
2005180220.tif 000016
【図４】
2005180220.tif 000017
【図５】
2005180220.tif 000018
【図６】
2005180220.tif 000019
【図７】
2005180220.tif 000020
【図８】
2005180220.tif 000021
【図９】
2005180220.tif 000022
【図１０】
2005180220.tif 000023
【図１１】
2005180220.tif 000024
【図１２】
2005180220.tif 000025
【図１３】
2005180220.tif 000026
【図１４】
2005180220.tif 000027
【図１５】
2005180220.tif 000028