JP2019007437A - 内燃機関の制御装置及び制御方法 - Google Patents

Abstract

【課題】排気管の放熱による排気ガスの温度低下を考慮することで、排気管の任意の推定位置の排気ガス温度を精度よく推定できる内燃機関の制御装置及び制御方法を提供する。【解決手段】運転状態に基づいて、燃焼室の出口における排気ガスの温度である出口ガス温度Ｔｏｕｔを算出する出口ガス温度算出部５２と、燃焼室の出口から推定位置までの排気管の放熱による排気ガスの温度低下量ΔＴｄを算出する放熱量算出部５３と、出口ガス温度Ｔｏｕｔから温度低下量ΔＴｄを減算して推定位置の排気ガス温度Ｔｅｘを推定する排気温度推定部５４と、を備えた内燃機関の制御装置５０。【選択図】図２

Description

本発明は、内燃機関の排気ガス温度を推定する内燃機関の制御装置及び制御方法に関する。

内燃機関を好適に運転するためには、内燃機関の運転状態に応じて変化する各部の圧力や温度等の状態量を正確に把握することが重要である。近年では、排気系の状態量も用いて内燃機関の制御が行われるようになっている。排気系の状態量を用いる制御としては、例えば、高回転高負荷運転時に排気ガス温度が上昇し、排気系に設けられた触媒や空燃比センサ等を損傷するおそれがある場合は、空燃比をリッチ化して排気ガス温度を低減させるエンリッチ制御が行われている。この制御を行うために用いる排気ガス温度は、排気ガス温度センサを用いて検出してもよいが、センサコストを低減するため、推定する方法もある。例えば、回転速度及び充填効率と排気ガス温度との関係が、予め実験データ等に基づいて設定されたマップデータを参照して、排気ガス温度を推定する方法が用いられる。

排気ガス温度を用いる他の例としては、例えば特許文献１において、体積効率相当値を用いて燃焼室の吸入空気量や内部ＥＧＲ率を算出する方法が開示されているが、この体積効率相当値の算出に排気ガス温度が用いられている。特許文献２では、過給機付の内燃機関を制御する方法が開示されているが、ここでは圧縮機の目標駆動力を達成するために必要なウェイストゲートの開度を制御するために排気ガス温度が用いられている。これらの文献では、前述同様のマップデータを用いた排気ガス温度の推定方法が用いられている。これら以外にも、ＥＧＲバルブの開度に基づいて外部ＥＧＲの流量を算出するために排気ガス温度が用いられる場合もある。

上記の例では推定した排気ガス温度が用いられるが、排気ガス温度センサを用いる例としては、排気ガス規制の強化に伴い内燃機関の排気ガス中に含まれるＰＭ（Ｐａｒｔｉｃｕｌａｔｅ Ｍａｔｔｅｒ）の排出量を低減する制御が挙げられる。具体的には排気系に取り付けられたＤＰＦ（Ｄｉｅｓｅｌ Ｐａｒｔｉｃｕｌａｔｅ Ｆｉｌｔｅｒ）やＧＰＦ（Ｇａｓｏｌｉｎｅ Ｐａｒｔｉｃｕｌａｔｅ Ｆｉｌｔｅｒ）等の微粒子捕集フィルタによりＰＭを捕集するが、微粒子捕集フィルタが捕集できるＰＭの量には限界があるため、微粒子捕集フィルタに堆積したＰＭを、燃焼させる再生処理が適宜実行される。この再生処理を適正に行うために排気温度センサを用いて排気ガス温度が制御される。

特許第５４０９８３２号公報 特許第５９６３９２７号公報 特許第５３７３９５２号公報

前述のように、排気ガス温度の推定方法としてマップデータを用いた推定方法があるが、排気ガス温度は、点火時期やＥＧＲ量等の影響を受けやすく、例えば回転速度と充填効率をマップ軸としたマップ設定だけでは、点火時期やＥＧＲ量等の運転条件が変化すると精度良く推定することが難しいという問題がある。点火時期やＥＧＲ量が変化した場合の排気ガス温度をすべてマップとして記憶するとしても、マップの数が膨大になるし、適合に必要な工数も膨大となるため、マップによる排気ガス温度の推定では推定精度や適合工数の面で問題がある。

マップによる排気ガス温度推定を用いた場合に考えられる他の問題としては、エンリッチ制御において、温度推定誤差を見越して排気ガス温度を確実に低減させるために、エンリッチの頻度が高くなり燃費が悪化するという問題がある。排気ガス温度の推定誤差により、内部ＥＧＲ量や外部ＥＧＲ量の推定精度が低くなるため、ＥＧＲ量に基づく点火時期制御が適正に行われず、燃費効果を十分に得られないという問題がある。排気ガス温度の推定誤差により、過給圧を制御するウェイストゲート開度の制御が適正に行われない場合にも、ドライバの要求する加速感が得られないためにドライバビリティが悪化するという問題がある。これらの問題は、排気ガス温度の推定精度を向上することができれば解決され得る。

マップ以外の推定方法としては、例えば特許文献３に開示された方法がある。特許文献３には、断熱圧縮による温度上昇分と、排気ガス中に損失する熱量による温度上昇分と、に基づいて排気ガス温度を推定する方法が開示されている。特許文献３の方法では、燃焼室の出口での排気ガス温度に影響する排気損失率は考慮されているが、燃焼室から排出された後、排気管の放熱による温度低下が考慮されていない。そのため、燃焼室の出口から離れた排気管の任意の推定位置の排気ガス温度を推定しようとすると、推定位置に応じて変化する排気管の放熱量による推定誤差が大きくなる問題があった。

また、特許文献３の方法では、燃焼熱量に対する排気ガス中に損失する熱量の割合である排気損失率を、回転速度及び燃料流量をパラメータとして算出している。しかしながら、排気損失率は、燃料の燃焼により燃焼室で発生した燃焼熱量から、筒内圧による仕事として取り出せる図示仕事、及び燃焼室の壁面へ放熱される熱量である冷却損失を減算することにより算出できる割合である。特許文献３の技術では、図示仕事、冷却損失が考慮されていないので、排気損失率の算出精度が良くないという問題がある。

排気ガス温度を精度よく推定することは容易でないため、排気ガス温度センサが用いられることもある。しかし、排気ガス温度センサは、高温となる排気系に取り付けられるため、異常診断を行う必要がある。異状診断を行うために、排気ガス温度の推定値を用いることが考えられ、そのためには、結局、排気ガス温度を精度よく推定することが必要になる。

そこで、排気管の放熱による排気ガスの温度低下を考慮することで、排気管の任意の推定位置の排気ガス温度を精度よく推定できる内燃機関の制御装置及び制御方法が望まれる。

本発明に係る内燃機関の制御装置は、内燃機関の運転状態を検出する運転状態検出部と、前記運転状態に基づいて、燃焼室の出口における排気ガスの温度である出口ガス温度を算出する出口ガス温度算出部と、前記燃焼室の出口から推定位置までの排気管の放熱による前記排気ガスの温度低下量を算出する放熱量算出部と、前記出口ガス温度から前記温度低下量を減算して前記推定位置の排気ガス温度を推定する排気温度推定部と、を備えたものである。

本発明に係る内燃機関の制御方法は、内燃機関の運転状態を検出する運転状態検出ステップと、前記運転状態に基づいて、燃焼室の出口における排気ガスの温度である出口ガス温度を算出する出口ガス温度算出ステップと、前記燃焼室の出口から推定位置までの排気管の放熱による前記排気ガスの温度低下量を算出する放熱量算出ステップと、前記出口ガス温度から前記温度低下量を減算して前記推定位置の排気ガス温度を推定する排気温度推定ステップと、を実行するものである。

本発明に係る内燃機関の制御装置及び制御方法によれば、それぞれ物理現象が異なる、燃焼室の出口の排気ガス温度と、排気管の温度低下量とを分けて算出するので、それぞれの算出精度を向上させることができる。燃焼室の出口から推定位置までの排気管の放熱による排気ガスの温度低下量を算出するので、排気管の任意の推定位置の排気ガス温度を精度よく推定できる。

本発明の実施の形態１に係る内燃機関及び制御装置の概略構成図である。 本発明の実施の形態１に係る制御装置のブロック図である。 本発明の実施の形態１に係る制御装置のハードウェア構成図である。 本発明の実施の形態１に係る排気管の単流熱交換器モデルを概略的に示す図である。 本発明の実施の形態１に係る熱勘定の考え方を示す図である。 本発明の実施の形態１に係る出口ガス温度算出部の詳細ブロック図である。 本発明の実施の形態１に係る制御装置の処理を示すフローチャート図である。

実施の形態１．
実施の形態１に係る内燃機関１の制御装置５０（以下、単に制御装置５０と称す）について図面を参照して説明する。図１は、本実施の形態に係る内燃機関１の概略構成図であり、図２は、本実施の形態に係る制御装置５０のブロック図である。内燃機関１及び制御装置５０は、車両に搭載され、内燃機関１は、車両（車輪）の駆動力源となる。

１．内燃機関１の構成
図１に示すように、内燃機関１は、空気と燃料の混合気を燃焼する燃焼室２５を備えている。内燃機関１は、燃焼室２５に空気を供給する吸気管２３と、燃焼室２５で燃焼した排気ガスを排出する排気管１７とを備えている。燃焼室２５は、シリンダ（気筒）とピストンにより構成されている。以下では、燃焼室２５を気筒とも称す。内燃機関１は、ガソリンエンジンとされている。内燃機関１は、吸気管２３を開閉するスロットルバルブ６を備えている。スロットルバルブ６は、制御装置５０により制御される電気モータにより開閉駆動される電子制御式スロットルバルブとされている。スロットルバルブ６には、スロットルバルブ６の開度に応じた電気信号を出力するスロットル開度センサ７が設けられている。

スロットルバルブ６の上流側の吸気管２３には、吸気管２３に吸入される吸入空気流量に応じた電気信号を出力するエアフローセンサ３と、吸入空気の温度に応じた電気信号を出力する吸入空気温度センサ４と、が設けられている。吸入空気温度センサ４に検出された吸入空気の温度は、外気温Ｔａに等しいとみなすことができる。

内燃機関１は、排気管１７から吸気マニホールド１２に排気ガスを還流するＥＧＲ流路２１と、ＥＧＲ流路２１を開閉するＥＧＲバルブ２２と、を有している。吸気マニホールド１２は、スロットルバルブ６の下流側の吸気管２３の部分である。ＥＧＲバルブ２２は、制御装置５０により制御される電気モータにより開閉駆動される電子制御式ＥＧＲバルブとされている。ＥＧＲバルブ２２には、ＥＧＲバルブ２２の開度に応じた電気信号を出力するＥＧＲ開度センサ２７が設けられている。なお、ＥＧＲは、排気ガス再循環、すなわち、Ｅｘｈａｕｓｔ Ｇａｓ Ｒｅｃｉｒｃｕｌａｔｉｏｎの頭文字である。ＥＧＲバルブ２２を介して排気ガスが再循環するＥＧＲを、外部ＥＧＲといい、吸排気バルブのバルブオーバーラップにより燃焼室内に排気ガスが残留するＥＧＲを、内部ＥＧＲという。以下、外部ＥＧＲを単にＥＧＲと称す。

吸気マニホールド１２には、吸気マニホールド１２内の気体の圧力であるマニホールド圧に応じた電気信号を出力するマニホールド圧センサ８と、吸気マニホールド１２内の気体の温度であるマニホールド温度Ｔｉｎに応じた電気信号を出力するマニホールド温度センサ９と、が設けられている。

内燃機関１には、燃焼室２５に燃料を供給するインジェクタ１３が設けられている。インジェクタ１３は、燃焼室２５内に直接燃料を噴射するように設けられている。インジェクタ１３は、吸気マニホールド１２の下流側の部分に燃料を噴射するように設けられてもよい。内燃機関１には、大気圧に応じた電気信号を出力する大気圧センサ２が設けられている。

燃焼室２５の頂部には、空気と燃料の混合気に点火する点火プラグと、点火プラグに点火エネルギーを供給する点火コイル１６と、が設けられている。また、燃焼室２５の頂部には、吸気管２３から燃焼室２５内に吸入される吸入空気量を調節する吸気バルブ１４と、燃焼室２５から排気管１７に排出される排気ガス量を調節する排気バルブ１５と、が設けられている。吸気バルブ１４には、そのバルブ開閉タイミングを可変にする吸気可変バルブタイミング機構が設けられている。排気バルブ１５には、そのバルブ開閉タイミングを可変にする排気可変バルブタイミング機構が設けられている。可変バルブタイミング機構１４、１５は、電動アクチュエータを有している。内燃機関１のクランク軸には、その回転角に応じた電気信号を出力するクランク角センサ２０が設けられている。

排気管１７には、排気ガス中の空気と燃料との比率である空燃比ＡＦ（Ａｉｒ／Ｆｕｅｌ）に応じた電気信号を出力する空燃比センサ１８が設けられている。また、排気管１７には、排気ガスを浄化する触媒１９が設けられている。

２．制御装置５０の構成
次に、制御装置５０について説明する。制御装置５０は、内燃機関１を制御対象とする制御装置である。図２のブロック図に示すように、制御装置５０は、運転状態検出部５１、出口ガス温度算出部５２、放熱量算出部５３、排気温度推定部５４、及び排気温度利用制御部５５等の制御部を備えている。制御装置５０の各制御部５１〜５５等は、制御装置５０が備えた処理回路により実現される。具体的には、制御装置５０は、図３に示すように、処理回路として、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）等の演算処理装置９０（コンピュータ）、演算処理装置９０とデータのやり取りをする記憶装置９１、演算処理装置９０に外部の信号を入力する入力回路９２、及び演算処理装置９０から外部に信号を出力する出力回路９３等を備えている。

演算処理装置９０として、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）、ＩＣ（Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）、ＤＳＰ（Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｏｒ）、ＦＰＧＡ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）、各種の論理回路、及び各種の信号処理回路等が備えられてもよい。また、演算処理装置９０として、同じ種類のもの又は異なる種類のものが複数備えられ、各処理が分担して実行されてもよい。記憶装置９１として、演算処理装置９０からデータを読み出し及び書き込みが可能に構成されたＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）や、演算処理装置９０からデータを読み出し可能に構成されたＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）等が備えられている。入力回路９２は、各種のセンサやスイッチが接続され、これらセンサやスイッチの出力信号を演算処理装置９０に入力するＡ／Ｄ変換器等を備えている。出力回路９３は、電気負荷が接続され、これら電気負荷に演算処理装置９０から制御信号を出力する駆動回路等を備えている。

そして、制御装置５０が備える各制御部５１〜５５等の各機能は、演算処理装置９０が、ＲＯＭ等の記憶装置９１に記憶されたソフトウェア（プログラム）を実行し、記憶装置９１、入力回路９２、及び出力回路９３等の制御装置５０の他のハードウェアと協働することにより実現される。なお、各制御部５１〜５５等が用いる特性データ、定数等の設定データは、ソフトウェア（プログラム）の一部として、ＲＯＭ等の記憶装置９１に記憶されている。

本実施の形態では、入力回路９２には、大気圧センサ２、エアフローセンサ３、吸入空気温度センサ４、スロットル開度センサ７、マニホールド圧センサ８、マニホールド温度センサ９、空燃比センサ１８、クランク角センサ２０、アクセルポジションセンサ２６、及びＥＧＲ開度センサ２７等が接続されている。出力回路９３には、スロットルバルブ６（電気モータ）、インジェクタ１３、吸気可変バルブタイミング機構１４、排気可変バルブタイミング機構１５、点火コイル１６、及びＥＧＲバルブ２２（電動アクチュエータ）等が接続されている。なお、制御装置５０には、図示していない各種のセンサ、スイッチ、及びアクチュエータ等が接続されている。

運転状態検出部５１は、内燃機関１の運転状態を検出する。運転状態検出部５１は、各種のセンサの出力信号等に基づいて各種の運転状態を検出する。具体的には、運転状態検出部５１は、大気圧センサ２の出力信号に基づいて大気圧を検出し、エアフローセンサ３の出力信号に基づいて吸入空気流量を検出し、吸入空気温度センサ４の出力信号に基づいて外気温Ｔａを検出し、スロットル開度センサ７の出力信号に基づいてスロットル開度を検出し、マニホールド圧センサ８の出力信号に基づいてマニホールド圧を検出し、マニホールド温度センサ９の出力信号等に基づいて吸気マニホールド１２内の気体の温度であるマニホールド温度Ｔｉｎを検出し、空燃比センサ１８の出力信号に基づいて、排気ガスの空燃比ＡＦを検出し、クランク角センサ２０の出力信号に基づいてクランク角度及び回転速度Ｎｅを検出し、アクセルポジションセンサ２６の出力信号に基づいてアクセル開度を検出し、ＥＧＲ開度センサ２７の出力信号に基づいてＥＧＲ開度を検出する。

運転状態検出部５１は、吸入空気流量、回転速度Ｎｅ等に基づいて燃焼室２５内に流入した空気量である吸入空気量Ｑｃ［ｇ／ｓｔｒｏｋｅ］、充填効率Ｅｃ［％］を算出する。例えば、運転状態検出部５１は、吸入空気流量［ｇ／ｓ］に、回転速度Ｎｅに応じた行程周期を乗算した値に、吸気マニホールドの遅れを模擬したフィルタ処理を行った値を、吸入空気量Ｑｃ［ｇ／ｓｔｒｏｋｅ］として算出する。行程周期は、３気筒エンジンであれば、２４０ｄｅｇＣＡ間の周期となり、４気筒エンジンであれば、１８０ｄｅｇＣＡ間の周期となる。或いは、運転状態検出部５１は、マニホールド圧、回転速度Ｎｅ等に基づいて吸入空気量Ｑｃ［ｇ／ｓｔｒｏｋｅ］、充填効率Ｅｃ［％］を算出してもよい。

運転状態検出部５１は、ＥＧＲ開度等に基づいて、燃焼室２５内に流入した排気ガス再循環量であるＥＧＲ量Ｑｃｅ［ｇ／ｓｔｒｏｋｅ］を算出する。例えば、運転状態検出部５１は、ＥＧＲ開度及びマニホールド圧等に基づいて、ＥＧＲバルブ２２を通過するＥＧＲ流量［ｇ／ｓ］を算出し、ＥＧＲ流量に行程周期を乗算した値に、フィルタ処理を行った値を、ＥＧＲ量Ｑｃｅ［ｇ／ｓｔｒｏｋｅ］として算出する。運転状態検出部５１は、吸入空気量Ｑｃに対するＥＧＲ量Ｑｃｅの比率であるＥＧＲ率Ｒｅｇｒ［％］を算出する。

制御装置５０は、算出された回転速度Ｎｅ、吸入空気量Ｑｃ、充填効率Ｅｃ、ＥＧＲ率Ｒｅｇｒ等に基づいて、目標空燃比、燃料噴射量Ｑｆ、点火時期ＳＡ等を算出し、インジェクタ１３及び点火コイル１６等を駆動制御する。例えば、制御装置５０は、回転速度Ｎｅ及び充填効率Ｅｃに基づいて目標空燃比を算出し、吸入空気量Ｑｃを目標空燃比で除算して燃料噴射量の基本値を算出する。そして、制御装置５０は、空燃比センサ１８を用いた空燃比フィードバック制御を行う場合は、燃料噴射量の基本値を補正して、最終的な燃料噴射量Ｑｆを算出し、空燃比フィードバック制御を行わない場合は、燃料噴射量の基本値をそのまま最終的な燃料噴射量Ｑｆに設定する。

制御装置５０は、アクセル開度等に基づいて、目標スロットル開度を算出し、スロットル開度が目標スロットル開度に近づくように、スロットルバルブ６の電気モータを駆動制御する。制御装置５０は、回転速度Ｎｅ及び充填効率Ｅｃ等に基づいてＥＧＲバルブ２２の目標ＥＧＲ開度を算出し、ＥＧＲ開度が目標ＥＧＲ開度に近づくように、ＥＧＲバルブ２２の電動アクチュエータを駆動制御する。制御装置５０は、回転速度Ｎｅ及び充填効率Ｅｃ等に基づいて、吸気バルブ１４及び排気バルブ１５のそれぞれの目標開閉タイミング（位相）を算出し、吸気バルブ１４及び排気バルブ１５の開閉タイミングが、それぞれの目標開閉タイミングに近づくように、吸気及び排気可変バルブタイミング機構の電動アクチュエータを駆動制御する。

トルクベース制御を行う場合は、制御装置５０は、アクセル開度等に基づいて算出された内燃機関１の要求出力トルク、又は変速機制御装置等の外部の制御装置から要求された内燃機関１の要求出力トルクを実現するように、スロットル開度、点火時期ＳＡ、ＥＧＲ開度、吸気バルブ１４及び排気バルブ１５の開閉タイミングを制御する。具体的には、制御装置５０は、適合により予め設定された充填効率Ｅｃと図示熱効率ηｉとの関係に基づいて、要求出力トルクに対応する図示熱効率を実現するために必要な目標充填効率を算出し、目標充填効率で運転した場合に、燃費、排出ガスが最適となる目標ＥＧＲ率を算出し、目標充填効率及び目標ＥＧＲ率を達成するような目標スロットル開度、目標ＥＧＲ開度、吸気バルブ１４の目標開閉タイミング、及び排気バルブ１５の目標開閉タイミングを算出し、これらの目標値に基づいて、スロットルバルブ６の電気モータ、ＥＧＲバルブ２２の電動アクチュエータ、吸気及び排気可変バルブタイミング機構の電動アクチュエータを駆動制御する。

２−１．排気ガス温度の推定演算
制御装置５０は、以上のような通常のエンジン制御の他に、排気ガス温度Ｔｅｘの推定を行うように構成されている。推定された排気ガス温度Ｔｅｘは、後述する排気温度利用制御部５５において、ＥＧＲ率の算出等に用いられる。

＜単流熱交換器モデルを用いた排気ガス温度の算出方法＞
まず、単流熱交換器モデルを用いた、排気ガス温度の算出方法について説明する。なお、単流熱交換器は、詳しくは、「大学講義 伝熱工学」（丸善株式会社、Ｐ２２４〜２２６、１９８３年）に記載されている。

図４に、単流熱交換器モデルを概略的に示す。図中の推定位置に、適合用の排気温度センサを取り付けて、実験データを収集する。適合用の排気温度センサは、試験用の内燃機関１にのみ装着されて出口ガス温度算出部５２、放熱量算出部５３、排気温度推定部５４等に含まれる各種定数の適合に用いられる。生産された内燃機関１では、出口ガス温度算出部５２、放熱量算出部５３、排気温度推定部５４により推定された排気ガス温度Ｔｅｘを用いて各種制御が行われる。

また、燃焼室出口及び排気管入口は、燃焼室と排気ポートとの境界部分に相当し、排気バルブ１５が配置された位置に相当する。燃焼室出口から推定位置（適合用の排気温度センサの取付位置）までを円筒と仮定する。燃焼室出口から推定位置までの排気管の内部表面積が、排気ガスから排気管に伝熱される全伝熱面積Ａ０［ｍ^２］になる。なお、複数気筒の場合は、全気筒の燃焼室出口から推定位置までの全内部表面積Ａ０を、気筒数で割ったものを用いることができる。排気管を流れる排気ガス流量をＱｅｘ［ｇ／ｓ］とする。排気管の外部は外気（排気管近傍の空気）により冷却されており、その外気温Ｔａ［Ｋ］は、一定であると仮定する。

次に、燃焼室出口の排気ガスの温度である出口ガス温度をＴｏｕｔ［Ｋ］とし、この出口ガス温度Ｔｏｕｔと外気温Ｔａとの差をθ１とする。推定位置における排気ガス温度をＴｅｘ［Ｋ］とし、排気ガス温度Ｔｅｘと外気温Ｔａとの差をθ２とする。なお、排気ガス温度Ｔｅｘは、排気温度センサの応答遅れの影響がない温度であり、排気ガス温度の瞬時値となる。また、全伝熱面積Ａ０のうち燃焼室出口から任意の位置までの面積をＡ［ｍ^２］とし、その任意の位置における瞬時排気温度をＴ［Ｋ］とし、任意の位置の微小伝熱面積ｄＡにおける排気ガス温度の変化をｄＴとする。瞬時排気温度Ｔと外気温Ｔａとの差をθとすると、ｄθ＝ｄＴである。さて、この微小伝熱面積ｄＡで単位時間に交換される熱量ｄＱは、熱貫流率（熱通過率）Ｋｈｔ［Ｗ／（ｍ^２・Ｋ）］を用いて表すと次式となる。

排気ガス流量Ｑｅｘ［ｇ／ｓ］と排気ガスの比熱Ｃｅｘ［Ｊ／（ｇ・Ｋ）］の積を水当量という。この水当量の排気ガスがｄＱ［Ｊ］の熱を失って、ｄＴだけ温度が下がるので、次式で表せる。

（１）式、（２）式から、ｄＱを削除して整理し、積分すると次式となる。ここで、Ｃｏｎｓｔは、積分定数である。

燃焼室出口では、Ａ＝０、θ＝θ１であり、推定位置では、Ａ＝Ａ０、θ＝θ２であるので、式（３）に適用して、式変形すると次式を得る。

この排気管において外気に奪われる熱量Ｑ［Ｊ］は次式で表される。

また、この排気管において外気に奪われる熱量の最大値Ｑｍａｘは、θ２＝０の場合（排気ガス温度Ｔｅｘが外気温Ｔａまで冷却された場合）なので、次式となる。

（４）式、（５）式、（６）式から、排気管の温度効率ηは次式で表される。

（７）式において、熱貫流率Ｋｈｔ、燃焼室出口から推定位置までの全伝熱面積Ａ０、及び排気ガスの比熱Ｃｅｘが、一定値であると仮定とすると、排気管の温度効率ηは排気ガス流量Ｑｅｘの関数であることがわかる。結局、排気ガス流量Ｑｅｘに対応する排気管の温度効率η［％］、外気温Ｔａ、出口ガス温度Ｔｏｕｔがわかれば、排気ガス温度Ｔｅｘは、次式により推定することができる。

（８）式の第２式の右辺第２項が、排気系に単流熱交換器モデルを適用した場合の温度低下量ΔＴｄを表すこととなる。以上により、出口ガス温度Ｔｏｕｔ、外気温Ｔａ、排気管の温度効率ηから排気ガス温度Ｔｅｘを算出する方法が示された。

＜出口ガス温度Ｔｏｕｔの算出方法＞
（８）式を用いて排気ガス温度Ｔｅｘを算出するために必要な、出口ガス温度Ｔｏｕｔの算出方法について説明する。燃焼室内に供給された燃料が燃焼することにより発生する燃焼熱量ｈｌ＿ｂがどのように分配されるかを示したものを熱勘定といい、熱勘定の考え方を図５に示す。燃焼熱量ｈｌ＿ｂのうち、燃焼室内の筒内圧による仕事として取り出せる熱量の割合を図示熱効率ηｉ［％］とし、膨張行程中に燃焼室の壁面へ放熱される熱量の割合を冷却損失率ηｃ［％］とし、残りが排気損失率ηｅｘ［％］として排気ガスの温度上昇に使われる熱量になる割合とする。なお、ポンピングロスや機械損失は、一旦仕事として取り出したものの、軸出力以外に使われる仕事と考えて、図示熱効率ηｉに含むものとする。ここで、図示熱効率ηｉ、冷却損失率ηｃが分かるものとすると、排気損失率ηｅｘは、次式で表せる。

この時、吸気マニホールド１２に導入された気体の温度であるマニホールド気温をＴｉｎ［Ｋ］とし、今回の燃焼により発生する燃焼熱量をｈｌ＿ｂ［Ｊ／ｓｔｒｏｋｅ］とし、今回の燃焼における燃焼室内ガス量をＱａｌｌ［ｇ／ｓｔｒｏｋｅ］とし、燃焼室内ガス（排気ガス）の比熱をＣｅｘ［Ｊ／（ｇ・Ｋ）］とすると、出口ガス温度の基本値Ｔｏｕｔ０［Ｋ］は、次式で算出できる。

ここで、（１０）式で算出する出口ガス温度を基本値としたのは、空燃比が理論空燃比又はリーンである場合は、（１０）式で算出した出口ガス温度をそのまま用いることができるが、空燃比がリッチである場合は、エンリッチにより、出口ガス温度が低下するため、その分の補正を行う必要があるためである。エンリッチによる温度低下の理由は、燃料の気化による蒸発熱や、未燃燃料が分子量の小さい炭化水素に分解される際に消費されるエネルギーによるものと考えられる。結局、エンリッチによる温度低下量をエンリッチ温度低下量ΔＴｒｉｃｈとすると、出口ガス温度Ｔｏｕｔは次式により算出することができる。以上により出口ガス温度Ｔｏｕｔが算出できるので、排気ガス温度Ｔｅｘも算出できることになる。

＜出口ガス温度算出部５２＞
以上で導出した算出方法に基づいて設計された制御装置５０について説明する。出口ガス温度算出部５２は、運転状態に基づいて、燃焼室の出口における排気ガスの温度である出口ガス温度Ｔｏｕｔを算出する。図６に、出口ガス温度算出部５２の詳細ブロック図を示すように、出口ガス温度算出部５２は、燃焼室ガス量算出部５２ａ、燃焼熱量算出部５２ｂ、冷却損失率算出部５２ｃ、図示熱効率算出部５２ｄ、熱勘定温度算出部５２ｅ、及び空燃比温度補正部５２ｆを備えている。

燃焼室ガス量算出部５２ａは、燃焼室に流入したガス量である燃焼室ガス量Ｑａｌｌを算出する。本実施の形態では、燃焼室ガス量算出部５２ａは、次式に示すように、燃焼室に流入した空気量である吸入空気量Ｑｃ［ｇ／ｓｔｒｏｋｅ］、燃焼室に流入した排気ガス再循環量であるＥＧＲ量Ｑｃｅ［ｇ／ｓｔｒｏｋｅ］、及び燃焼室に供給した燃料噴射量Ｑｆ［ｇ／ｓｔｒｏｋｅ］を合計した値を、燃焼室ガス量Ｑａｌｌ［ｇ／ｓｔｒｏｋｅ］として算出する。なお、吸入空気量Ｑｃ及びＥＧＲ量Ｑｃｅは、上述したように、運転状態検出部５１により算出される。

また、燃焼室ガス量算出部５２ａは、次式に示すように、燃焼室ガス量Ｑａｌｌを行程周期ΔＴｓｇｔで除算した値を、排気行程で排出される排気ガス流量Ｑｅｘ［ｇ／ｓ］として算出する。行程周期ΔＴｓｇｔは、３気筒エンジンであれば、２４０ｄｅｇＣＡ間の周期となり、４気筒エンジンであれば、１８０ｄｅｇＣＡ間の周期となる。ここで、（ｎ）は、現在の行程を表し、（ｎ−３）は、現在の行程から３行程前の行程を表している。Ｑｅｘ（ｎ）は、現在排気行程である気筒から排出されている排気ガス流量を表し、Ｑａｌｌ（ｎ−３）は、現在排気行程である気筒が、３行程前に吸気行程であった時に流入した燃焼室ガス量を表している。

燃焼熱量算出部５２ｂは、燃料の燃焼により燃焼室で発生した燃焼熱量ｈｌ＿ｂを算出する。燃焼熱量算出部５２ｂは、燃料噴射量Ｑｆ［ｇ／ｓｔｒｏｋｅ］及び空燃比ＡＦに基づいて、燃料噴射量Ｑｆのうち、実際に燃焼した燃焼燃料量Ｑｆｂ［ｇ／ｓｔｒｏｋｅ］を算出する。空燃比ＡＦが理論空燃比又はリーンである場合は、燃料噴射量Ｑｆの全てが燃焼すると仮定できるが、空燃比ＡＦがリッチである場合は、理論空燃比分の燃料量は燃焼するが、理論空燃比を超えたリッチ分の燃料量は燃焼しない。そこで、燃焼熱量算出部５２ｂは、次式に示すように、空燃比ＡＦが理論空燃比又はリーンである場合は、燃料噴射量Ｑｆをそのまま燃焼燃料量Ｑｆｂに設定し、空燃比ＡＦがリッチである場合は、燃料噴射量Ｑｆに、空燃比ＡＦを理論空燃比ＡＦ０で除算した値を、燃焼燃料量Ｑｆｂとして算出する。空燃比ＡＦには、燃料噴射量Ｑｆの演算に用いられる目標空燃比が用いられてもよいし、空燃比センサ１８を用いて検出された排気ガスの空燃比が用いられてもよい。

そして、燃焼熱量算出部５２ｂは、燃焼燃料量Ｑｆｂ［ｇ／ｓｔｒｏｋｅ］に単位発熱量を乗算した値を、燃焼熱量ｈｌ＿ｂ［Ｊ／ｓｔｒｏｋｅ］として算出する。ガソリンの単位発熱量は、４４０００［Ｊ／ｇ］程度に設定される。なお、燃料の気化熱量を、単位発熱量から減算してもよいが、ガソリンの気化熱量は２７２［Ｊ／ｇ］程度であるので無視してもよい。

冷却損失率算出部５２ｃは、燃焼により発生した燃焼熱量ｈｌ＿ｂのうち、燃焼室の壁面に放熱される熱量の割合である冷却損失率ηｃを算出する。本願の発明者が実施した実験結果から、排気ガス流量Ｑｅｘと、逆算により算出した冷却損失率ηｃとの間には、運転条件に依らず、強い相関関係があることが見いだされた。そこで、冷却損失率算出部５２ｃは、排気管内の排気ガス流量Ｑｅｘに基づいて冷却損失率ηｃを算出する。具体的には、冷却損失率算出部５２ｃは、排気ガス流量Ｑｅｘと冷却損失率ηｃとの関係が予め設定された損失率特性データを参照し、現在の排気ガス流量Ｑｅｘ（ｎ）に対応する冷却損失率ηｃ（ｎ）を参照する。損失率特性データには、排気ガス流量Ｑｅｘが増加するに従って、冷却損失率ηｃが低下する特性が、実験データに基づいて予め設定されている。なお、各特性データには、データマップ、データテーブル、多項式、数式等が用いられ、その設定データは、記憶装置９１に記憶されている。

図示熱効率算出部５２ｄは、燃焼熱量ｈｌ＿ｂのうち、燃焼室内の筒内圧による仕事として取り出せる熱量の割合である図示熱効率ηｉを算出する。図示熱効率算出部５２ｄは、回転速度Ｎｅ、充填効率Ｅｃ、ＥＧＲ率Ｒｅｇｒ、及び点火時期ＳＡ等の運転条件と、図示熱効率ηｉとの関係が予め設定された図示熱効率特性データを参照し、現在の運転条件に対応する図示熱効率ηｉを算出する。ここで、現在の充填効率Ｅｃ、ＥＧＲ率Ｒｅｇｒには、現在排気行程である気筒が、３行程前に吸気行程であった時の値が用いられる。

図示熱効率特性データは、トルクベース制御において出力トルクの演算に用いるものが用いられ、実験データに基づいて予め設定されている。図示熱効率ηｉの計測は、全運転ポイントを計測すると膨大となるため、ＭＢＣ（Ｍｏｄｅｌ Ｂａｓｅｄ Ｃａｌｉｂｒａｔｉｏｎ）、ＤｏＥ（Ｄｅｓｉｇｎ ｏｆ Ｅｘｐｅｒｉｍｅｎｔｓ）といった手法を用いて計測ポイントを抑制し、計測したポイントでのデータを近似式で繋ぐようにして全運転ポイントに相当する近似データを作成する。図示熱効率特性データには、更に簡易化した近似式が用いられる。なお、トルクベース制御が行われない場合は、図示熱効率特性データは、回転速度Ｎｅ及び充填効率Ｅｃと、図示熱効率ηｉとの関係が設定されたものとされてもよい。なお、図示熱効率ηｉには、燃焼により発生する熱量を分母にし、図示平均有効圧に相当する熱量を分子にして算出したものが用いられる。

熱勘定温度算出部５２ｅは、燃焼室ガス量Ｑａｌｌ［ｇ／ｓｔｒｏｋｅ］、燃焼熱量ｈｌ＿ｂ［Ｊ／ｓｔｒｏｋｅ］、冷却損失率ηｃ、及び図示熱効率ηｉに基づいて出口ガス温度Ｔｏｕｔを算出する。具体的には、熱勘定の考え方に基づく（９）式、（１０）式を用いて上述したように、熱勘定温度算出部５２ｅは、次式に示すように、１００［％］から図示熱効率ηｉ［％］及び冷却損失率ηｃ［％］を減算した値を、燃焼熱量のうち排気ガスの温度上昇に使われる熱量の割合である排気損失率ηｅｘ［％］として算出する。そして、熱勘定温度算出部５２ｅは、燃焼熱量ｈｌ＿ｂ［Ｊ／ｓｔｒｏｋｅ］に、排気損失率ηｅｘ［％］を乗算して、排気ガスの温度上昇に用いられる温度上昇熱量［Ｊ／ｓｔｒｏｋｅ］を算出し、温度上昇熱量［Ｊ／ｓｔｒｏｋｅ］を、燃焼室ガス量Ｑａｌｌ［ｇ／ｓｔｒｏｋｅ］と排気ガスの比熱Ｃｅｘ［Ｊ／（ｇ・Ｋ）］を乗算した水当量［Ｊ／（ｓｔｒｏｋｅ・ｇ）］で除算して温度上昇量［Ｋ］を算出し、温度上昇量［Ｋ］をマニホールド温度Ｔｉｎ［Ｋ］に加算して、出口ガス温度の基本値Ｔｏｕｔ０［Ｋ］を算出する。

排気ガスの比熱Ｃｅｘには、１．１［Ｊ／（ｇ・Ｋ）］程度の値が設定される。排気ガスの比熱Ｃｅｘは、厳密には空燃比に応じて変化するため、空気及び燃料の比熱を用いて補正するなど、空燃比ＡＦに応じて変化させてもよいが、簡単のために、同じ固定値が用いられてもよい。マニホールド温度温Ｔｉｎとして、マニホールド温度センサ９により検出された吸気マニホールド１２内のガス温度がそのまま用いられてもよいが、燃焼室内に流入したガス温度に近づくように、吸気マニホールド１２内のガス温度に所定値が加算された温度がマニホールド温Ｔｉｎとして用いられてもよい。或いは、吸排気バルブの開閉タイミングに基づいて算出された内部ＥＧＲ率に応じて吸気マニホールド１２内のガス温度が補正された温度が、マニホールド温Ｔｉｎとして用いられてもよい。マニホールド温度温Ｔｉｎとして、外気温Ｔａに基づいて推定されたガス温度が用いられてもよい。

空燃比温度補正部５２ｆは、空燃比ＡＦが理論空燃比よりリッチである場合に、エンリッチ量ΔＡＦｒに応じて、出口ガス温度Ｔｏｕｔの低下補正を行う。本願の発明者が実施した実験結果から、運転条件に依らず、エンリッチ量ΔＡＦｒが１増加すると、出口ガス温度Ｔｏｕｔが３５〜４０℃低下する特性となった。空燃比温度補正部５２ｆは、エンリッチ量ΔＡＦｒ（＝ＡＦ０−ＡＦ）とエンリッチ温度低下量ΔＴｒｉｃｈとの関係が予め設定された温度低下特性データを参照し、現在のエンリッチ量ΔＡＦｒに対応するエンリッチ温度低下量ΔＴｒｉｃｈを算出する。温度低下特性データには、エンリッチ量ΔＡＦｒが増加するに従って、正の値のエンリッチ温度低下量ΔＴｒｉｃｈが増加する特性が、実験データに基づいて予め設定されている。

空燃比温度補正部５２ｆは、次式に示すように、空燃比ＡＦが理論空燃比又はリーンである場合は、出口ガス温度の基本値Ｔｏｕｔ０をそのまま出口ガス温度Ｔｏｕｔに設定し、空燃比ＡＦがリッチである場合は、出口ガス温度の基本値Ｔｏｕｔ０から、正の値のエンリッチ温度低下量ΔＴｒｉｃｈを減算した値が、出口ガス温度Ｔｏｕｔとして算出する。

＜放熱量算出部５３＞
放熱量算出部５３は、燃焼室の出口から推定位置までの排気管の放熱による排気ガスの温度低下量ΔＴｄを算出する。上述したように、排気管を単流熱交換器にモデル化して導出した（７）式から、排気管の温度効率ηは、排気ガス流量Ｑｅｘの関数であることがわかり、（８）式から、排気管の温度効率ηに基づいて、温度低下量ΔＴｄを算出できることがわかる。そこで、放熱量算出部５３は、排気管内の排気ガス流量Ｑｅｘに基づいて、排気管内の排気ガスを加熱流体とし、排気管外の空気を受熱流体とした単流熱交換器としての排気管の温度効率ηを算出し、排気管の温度効率ηに基づいて温度低下量ΔＴｄを算出する。

放熱量算出部５３は、（７）式と同様の次式を用いて、排気ガス流量Ｑｅｘに基づいて、排気管の温度効率ηを算出する。

Ｋηは、演算定数である。演算定数Ｋηは、排気管の熱貫流率Ｋｈｔに、燃焼室出口から推定位置までの全伝熱面積Ａ０を乗算し、排気ガスの比熱Ｃｅｘで除算することにより設定できる。排気管の熱貫流率Ｋｈｔは、適合値であるが、例えば、１０〜１５［ｗ／（ｍ^２・Ｋ）］程度の値になる。全伝熱面積Ａ０は、排気管の構造から算出できる。排気ガスの比熱Ｃｅｘには、上述した値が用いられる。演算定数Ｋηは、固定値とされてもよいが、空燃比ＡＦに応じて変化する排気ガスの比熱Ｃｅｘに応じて変化されてもよい。また、演算定数Ｋηは、実験による適合値とされてもよい。

或いは、放熱量算出部５３は、排気ガス流量Ｑｅｘと排気管の温度効率ηとの関係が予め設定された温度効率特性データを参照し、現在の排気ガス流量Ｑｅｘに対応する排気管の温度効率ηを算出するように構成されてもよい。温度効率特性データは、（１７）式の特性を表すものであるが、実験により適合されてもよい。例えば、複数の排気ガス流量Ｑｅｘの動作点において測定した、適合用の排気温度センサを用いた出口ガス温度Ｔｏｕｔの測定値、適合用の排気温度センサを用いた推定位置の排気ガス温度Ｔｅｘ０の測定値、外気温Ｔａの測定値に基づいて、（８）式の第１式を用い、排気管の温度効率ηを算出する。そして、複数の動作点の排気ガス流量Ｑｅｘ及び排気管の温度効率ηの実験データを近似して、温度効率特性データを設定する。同様の方法により、演算定数Ｋηが設定されてもよい。

放熱量算出部５３は、（８）式の第２式と同様の次式に示すように、出口ガス温度算出部５２により算出された出口ガス温度Ｔｏｕｔから外気温Ｔａを減算した値に、排気管の温度効率ηを乗算した値を、排気ガスの温度低下量ΔＴｄとして算出する。

＜排気温度推定部５４＞
排気温度推定部５４は、次式に示すように、出口ガス温度Ｔｏｕｔから温度低下量ΔＴｄを減算して推定位置の排気ガス温度Ｔｅｘを推定する。

推定位置は、後述する排気温度利用制御部５５において必要になる排気ガス温度の位置に設定される。例えば、推定位置は、触媒１９の上流側の位置、排気管１７とＥＧＲ流路２１との接続位置、排気管に過給機が設けられる場合はタービンの上流側の位置、排気管に排気ガスセンサが設けられる場合は排気ガスセンサの位置等に設定される。推定位置に合わせて、温度効率ηの算出に必要な演算定数Ｋη又は温度効率特性データが設定される。

排気温度推定部５４は、複数の推定位置の排気ガス温度を推定するように構成されてもよい。この場合は、放熱量算出部５３は、推定位置に応じて演算定数Ｋη又は温度効率特性データの設定値を切り替えて、各推定位置の温度効率ηを算出し、各推定位置の温度低下量ΔＴｄを算出する。そして、排気温度推定部５４は、各推定位置の温度低下量ΔＴｄを用いて、各推定位置の排気ガス温度Ｔｅｘを推定する。このように、推定位置に応じて演算定数Ｋη又は温度効率特性データを切り替えるだけで、容易に複数の推定位置の排気ガス温度を推定することができる。

後述するように、排気管に排気ガス温度センサが設けられ、排気ガス温度センサの異常診断を行う場合は、排気ガス温度センサにより検出した排気ガス温度と、推定した排気ガス温度Ｔｅｘとを比較する。しかし、排気ガス温度センサにより検出した排気ガス温度には、センサの熱容量等による応答遅れが生じる。或いは、排気ガス温度には、排気管の熱容量による応答遅れが生じる。そこで、排気温度推定部５４は、排気ガス温度Ｔｅｘに対して、応答遅れ処理を行うことにより、応答遅れ処理後の排気ガス温度Ｔｅｘｆｔを算出する。例えば、排気温度推定部５４は、次式に示す一次遅れフィルタ処理を行って、応答遅れ処理後の排気ガス温度Ｔｅｘｆｔを算出する。ここで、フィルタ定数Ｋｆは、センサの時定数τと演算周期Δｔから設定される。（ｎ）は、今回の演算周期の値であることを表し、（ｎ−１）は、前回の演算周期の値であることを表す。

＜排気温度利用制御部５５＞
排気温度利用制御部５５は、排気ガス温度の推定値Ｔｅｘを用いた、排気温度制御、バルブ流量特性演算、排気温度センサ異常診断、及びタービン出力演算のいずれか１つ以上を実行する。

排気温度制御は、排気ガス温度の推定値Ｔｅｘを用いて排気ガス温度の制御を行う処理である。排気温度推定部５４は、排気ガス温度の制御を行う位置を推定位置に設定して、排気ガス温度の推定を行う。例えば、排気温度利用制御部５５は、排気ガスの温度を低下させるためのエンリッチ制御を行っている場合に、排気ガス温度の推定値Ｔｅｘが目標温度に近づくように、燃料噴射のエンリッチ量を変化させる。また、排気温度利用制御部５５は、排気ガス温度の推定値Ｔｅｘが目標温度に近づくように、点火時期ＳＡを変化させたり、排気行程の燃料噴射量を変化させたりする。

バルブ流量特性演算は、排気ガス温度の推定値Ｔｅｘを用いて排気ガスが流れるバルブの流量特性を算出する処理である。排気温度推定部５４は、バルブの上流側の位置を推定位置に設定して、排気ガス温度の推定を行う。排気温度利用制御部５５は、排気ガス温度の推定値Ｔｅｘに基づいて、バルブの流量特性としてバルブの上流側の排気ガスの音速及び密度を算出する。排気ガスが流れるバルブは、ＥＧＲバルブ２２、過給機のタービンを迂回するウェイストゲートバルブ等とされる。排気温度利用制御部５５は、ＥＧＲバルブ２２の流量特性を用いて、ＥＧＲ流量を算出し、ＥＧＲ量Ｑｃｅ、ＥＧＲ率Ｒｅｇｒを算出する。排気温度利用制御部５５は、ウェイストゲートバルブの流量特性を用いて、ウェイストゲートバルブを通過するバイパス流量を算出し、排気ガス流量からバイパス流量を減算して、タービン通過流量を算出し、タービン通過流量はタービン出力に比例するため、タービン通過流量を用いてタービの出力を算出する。タービン出力は、過給圧の制御に用いられる。

排気温度センサ異常診断は、排気ガス温度の推定値Ｔｅｘを用いて排気管に設けられた排気ガス温度センサの異常診断を行う処理である。排気ガス温度センサは、触媒、微粒子捕集フィルタ等の排気ガスの浄化装置に流入する排気ガス温度を管理するために設けられる。排気温度推定部５４は、排気ガス温度センサの位置を推定位置に設定して、排気ガス温度の推定を行う。排気温度利用制御部５５は、応答遅れ処理後の排気ガス温度の推定値Ｔｅｘｆｔと、排気ガス温度センサによる排気ガス温度の検出値とを比較し、両者の差が大きい場合に、排気ガス温度センサに異常が発生したと判定する。

タービン出力演算は、排気ガス温度の推定値Ｔｅｘを用いて排気管に設けられた過給機のタービン出力を算出する演算である。排気温度推定部５４は、タービンの上流側の位置を推定位置に設定して、排気ガス温度の推定を行う。排気温度利用制御部５５は、タービンに流入する排気ガス温度はタービン出力に比例するため、排気ガス温度の推定値Ｔｅｘを用いてタービン出力を算出する。

２−２．フローチャート
本実施の形態に係る制御装置５０の概略的な処理の手順（内燃機関１の制御方法）について、図７に示すフローチャートに基づいて説明する。図７のフローチャートの処理は、演算処理装置９０が記憶装置９１に記憶されたソフトウェア（プログラム）を実行することにより、所定の演算周期毎に繰り返し実行される。

ステップＳ０１で、運転状態検出部５１は、上記のように、内燃機関１の各種の運転状態を検出する運転状態検出処理（運転状態ステップ）を実行する。

ステップＳ０２で、出口ガス温度算出部５２は、上記のように、運転状態に基づいて、燃焼室の出口における排気ガスの温度である出口ガス温度Ｔｏｕｔを算出する出口ガス温度算出処理（出口ガス温度算出ステップ）を実行する。本実施の形態では、出口ガス温度算出部５２は、上記のように、燃焼室に流入したガス量である燃焼室ガス量Ｑａｌｌを算出し、燃料の燃焼により燃焼室で発生した燃焼熱量ｈｌ＿ｂを算出し、燃焼により発生した燃焼熱量ｈｌ＿ｂのうち、燃焼室の壁面に放熱される熱量の割合である冷却損失率ηｃを算出し、燃焼熱量ｈｌ＿ｂのうち、燃焼室内の筒内圧による仕事として取り出せる熱量の割合である図示熱効率ηｉを算出し、燃焼室ガス量Ｑａｌｌ、燃焼熱量ｈｌ＿ｂ、冷却損失率ηｃ、及び図示熱効率ηｉに基づいて出口ガス温度Ｔｏｕｔを算出する。

ステップＳ０３で、放熱量算出部５３は、上記のように、燃焼室の出口から推定位置までの排気管の放熱による排気ガスの温度低下量ΔＴｄを算出する放熱量算出処理（放熱量算出ステップ）を実行する。本実施の形態では、放熱量算出部５３は、排気管内の排気ガス流量Ｑｅｘに基づいて、単流熱交換器としての排気管の温度効率ηを算出し、排気管の温度効率ηに基づいて温度低下量ΔＴｄを算出する。放熱量算出部５３は、（１７）式を用い、排気ガス流量Ｑｅｘに基づいて、排気管の温度効率ηを算出する。

ステップＳ０４で、排気温度推定部５４は、上記のように、出口ガス温度Ｔｏｕｔから温度低下量ΔＴｄを減算して推定位置の排気ガス温度Ｔｅｘを推定する排気温度推定処理（排気温度推定ステップ）を実行する。

ステップＳ０５で、排気温度利用制御部５５は、上記のように、排気ガス温度の推定値Ｔｅｘを用いた、排気温度制御、バルブ流量特性演算、排気温度センサ異常診断、及びタービン出力演算のいずれか１つ以上を行う排気温度利用制御処理（排気温度利用制御ステップ）を実行する。

〔その他の実施の形態〕
最後に、本発明のその他の実施の形態について説明する。なお、以下に説明する各実施の形態の構成は、それぞれ単独で適用されるものに限られず、矛盾が生じない限り、他の実施の形態の構成と組み合わせて適用することも可能である。

（１）上記の実施の形態１においては、内燃機関１は、ガソリンエンジンとされている場合を例として説明した。しかし、本発明の実施の形態はこれに限定されない。すなわち、内燃機関１は、ディーゼルエンジン、ＨＣＣＩ燃焼（Ｈｏｍｏｇｅｎｅｏｕｓ−Ｃｈａｒｇｅ Ｃｏｍｐｒｅｓｓｉｏｎ Ｉｇｎｉｔｉｏｎ Ｃｏｍｂｕｓｔｉｏｎ）を行うエンジン等の各種の内燃機関とされてもよい。

（２）上記の実施の形態１において、例示的に説明したように、内燃機関１は、過給機を備えてもよい。過給機は、排気管に設けられたタービンと、吸気管におけるスロットルバルブの上流側に設けられ、タービンと一体的に回転する圧縮機と、タービンを迂回するタービンバイパス通路に設けられたウェイストゲートバルブと、を備えている。また、内燃機関１は、触媒に加えて、微粒子捕集フィルタを備えてもよく、触媒、微粒子捕集フィルタの上流側に排気ガス温度センサを備えてもよい。

なお、本発明は、その発明の範囲内において、実施の形態を適宜、変形、省略したりすることが可能である。

１ 内燃機関、１７ 排気管、２５ 燃焼室、５０ 内燃機関の制御装置、５１ 運転状態検出部、５２ 出口ガス温度算出部、５３ 放熱量算出部、５４ 排気温度推定部、５５ 排気温度利用制御部、Ｋη 演算定数、Ｑａｌｌ 燃焼室ガス量、Ｑｅｘ 排気ガス流量、Ｒｅｇｒ ＥＧＲ率、Ｔａ 外気温、Ｔｅｘ 排気ガス温度、Ｔｏｕｔ 出口ガス温度、ｈｌ＿ｂ 燃焼熱量、ΔＴｄ 温度低下量、η 温度効率、ηｃ 冷却損失率、ηｅｘ 排気損失率、ηｉ 図示熱効率

本発明に係る内燃機関の制御装置は、内燃機関の運転状態を検出する運転状態検出部と、前記運転状態に基づいて、燃焼室の出口における排気ガスの温度である出口ガス温度を算出する出口ガス温度算出部と、前記燃焼室の出口から推定位置までの排気管の放熱による前記排気ガスの温度低下量を算出する放熱量算出部と、前記出口ガス温度から前記温度低下量を減算して前記推定位置の排気ガス温度を推定する排気温度推定部と、を備え、
前記放熱量算出部は、前記排気管の温度効率をηとし、前記排気管内の排気ガス流量をＱｅｘとし、演算定数をＫηとして、
η＝１−ｅｘｐ（−Ｋη／Ｑｅｘ）
の算出式により、前記温度効率を算出し、又は
前記排気ガス流量と前記温度効率との関係が予め設定された温度効率特性データを参照し、現在の前記排気ガス流量に対応する前記温度効率を算出し、
前記放熱量算出部は、複数の前記推定位置のそれぞれに応じて前記演算定数又は前記温度効率特性データの設定値を切り替えて、複数の前記推定位置のそれぞれの前記温度効率を算出し、前記出口ガス温度から外気温を減算した値に、複数の前記推定位置のそれぞれの前記温度効率を乗算した値を、複数の前記推定位置のそれぞれの前記温度低下量として算出し、
前記排気温度推定部は、前記出口ガス温度から複数の前記推定位置のそれぞれの前記温度低下量を減算して、複数の前記推定位置のそれぞれの前記排気ガス温度を推定するものである。

本発明に係る内燃機関の制御方法は、内燃機関の運転状態を検出する運転状態検出ステップと、前記運転状態に基づいて、燃焼室の出口における排気ガスの温度である出口ガス温度を算出する出口ガス温度算出ステップと、前記燃焼室の出口から推定位置までの排気管の放熱による前記排気ガスの温度低下量を算出する放熱量算出ステップと、前記出口ガス温度から前記温度低下量を減算して前記推定位置の排気ガス温度を推定する排気温度推定ステップと、を実行し、
前記放熱量算出ステップでは、前記排気管の温度効率をηとし、前記排気管内の排気ガス流量をＱｅｘとし、演算定数をＫηとして、
η＝１−ｅｘｐ（−Ｋη／Ｑｅｘ）
の算出式により、前記温度効率を算出し、又は
前記排気ガス流量と前記温度効率との関係が予め設定された温度効率特性データを参照し、現在の前記排気ガス流量に対応する前記温度効率を算出し、
前記放熱量算出ステップでは、複数の前記推定位置のそれぞれに応じて前記演算定数又は前記温度効率特性データの設定値を切り替えて、複数の前記推定位置のそれぞれの前記温度効率を算出し、前記出口ガス温度から外気温を減算した値に、複数の前記推定位置のそれぞれの前記温度効率を乗算した値を、複数の前記推定位置のそれぞれの前記温度低下量として算出し、
前記排気温度推定ステップでは、前記出口ガス温度から複数の前記推定位置のそれぞれの前記温度低下量を減算して、複数の前記推定位置のそれぞれの前記排気ガス温度を推定するものである。

Claims (10)

1. 内燃機関の運転状態を検出する運転状態検出部と、
   前記運転状態に基づいて、燃焼室の出口における排気ガスの温度である出口ガス温度を算出する出口ガス温度算出部と、
   前記燃焼室の出口から推定位置までの排気管の放熱による前記排気ガスの温度低下量を算出する放熱量算出部と、
   前記出口ガス温度から前記温度低下量を減算して前記推定位置の排気ガス温度を推定する排気温度推定部と、を備えた内燃機関の制御装置。
2. 前記放熱量算出部は、前記排気管内の排気ガス流量に基づいて、前記排気管内の排気ガスを加熱流体とし、前記排気管外の空気を受熱流体とした単流熱交換器としての前記排気管の温度効率を算出し、
   前記温度効率に基づいて前記温度低下量を算出する請求項１に記載の内燃機関の制御装置。
3. 前記放熱量算出部は、前記排気管の温度効率をηとし、前記排気管内の排気ガス流量をＱｅｘとし、演算定数をＫηとして、
   η＝１−ｅｘｐ（−Ｋη／Ｑｅｘ）
   の算出式により、前記温度効率を算出し、
   前記温度効率に基づいて前記温度低下量を算出する請求項１又は２に記載の内燃機関の制御装置。
4. 前記放熱量算出部は、前記排気管内の排気ガス流量と前記排気管の温度効率との関係が予め設定された温度効率特性データを参照し、現在の前記排気ガス流量に対応する前記温度効率を算出し、
   前記温度効率に基づいて前記温度低下量を算出する請求項１又は２に記載の内燃機関の制御装置。
5. 前記放熱量算出部は、前記出口ガス温度から外気温を減算した値に、前記温度効率を乗算した値を、前記温度低下量として算出する請求項２から４のいずれか一項に記載の内燃機関の制御装置。
6. 前記出口ガス温度算出部は、
   前記燃焼室に流入したガス量である燃焼室ガス量を算出し、
   燃料の燃焼により前記燃焼室で発生した燃焼熱量を算出し、
   前記燃焼熱量のうち、前記燃焼室の壁面に放熱される熱量の割合である冷却損失率を算出し、
   前記燃焼熱量のうち、前記燃焼室内の筒内圧による仕事として取り出せる熱量の割合である図示熱効率を算出し、
   前記燃焼室ガス量、前記燃焼熱量、前記冷却損失率、及び前記図示熱効率に基づいて前記出口ガス温度を算出する請求項１から５のいずれか一項に記載の内燃機関の制御装置。
7. 前記出口ガス温度算出部は、前記排気管内の排気ガス流量に基づいて、前記冷却損失率を算出する請求項６に記載の内燃機関の制御装置。
8. 前記出口ガス温度算出部は、空燃比が理論空燃比よりリッチである場合に、エンリッチ量に応じて、前記出口ガス温度の低下補正を行う請求項１から７のいずれか一項に記載の内燃機関の制御装置。
9. 前記排気ガス温度の推定値を用いて前記排気ガス温度の制御を行う排気温度制御、前記排気ガス温度の推定値を用いて排気ガスが流れるバルブの流量特性を算出するバルブ流量特性演算、前記排気ガス温度の推定値を用いて前記排気管に設けられた排気ガス温度センサの異常診断を行う排気温度センサ異常診断、及び前記排気ガス温度の推定値を用いて前記排気管に設けられた過給機のタービン出力を算出するタービン出力演算のいずれか１つ以上を実行する排気温度利用制御部を更に備えた請求項１から８のいずれか一項に記載の内燃機関の制御装置。
10. 内燃機関の運転状態を検出する運転状態検出ステップと、
    前記運転状態に基づいて、燃焼室の出口における排気ガスの温度である出口ガス温度を算出する出口ガス温度算出ステップと、
    前記燃焼室の出口から推定位置までの排気管の放熱による前記排気ガスの温度低下量を算出する放熱量算出ステップと、
    前記出口ガス温度から前記温度低下量を減算して前記推定位置の排気ガス温度を推定する排気温度推定ステップと、を実行する内燃機関の制御方法。

Images