(57)【要約】【課題】燃料フィルタの高い濾過効率を継続的に確保する燃料系制御装置の提供。【解決手段】Ｓ１０２，Ｓ１０３では、燃料の性状を表す性状パラメータとして現在動粘度ＫＶｐ及び現在密度ＤＤｐを取得する。Ｓ１０４では、Ｓ１０２，Ｓ１０３により取得された性状パラメータとしての現在動粘度ＫＶｐ及び現在密度ＤＤｐに基づき、燃料を構成する成分の平均元素数として平均炭素数ＮＣ及び平均水素数ＮＨを算出する。Ｓ１０５では、Ｓ１０４により算出された平均炭素数ＮＣ及び平均水素数ＮＨに対応する燃料の曇点ＴＭｃを、推定する。Ｓ１０６，Ｓ１０７では、Ｓ１０５により推定された曇点ＴＭｃ以下となる燃料の判定時温度ＴＭｊにおいて、ヒータを制御することにより燃料フィルタを加熱する。【選択図】図２
2018200024.tif 000002
【特許請求の範囲】
【請求項１】
内燃機関（１）の燃料系（１０）において燃焼室（１ａ）へ供給される燃料を濾過する燃料フィルタ（１１）のヒータ（１１０）による加熱を制御する燃料系制御装置（５０）であって、
前記燃料の性状を表す性状パラメータ（ＫＶｐ，ＤＤｐ，ＫＶｐ１，ＫＶｐ２，ＤＤｐ１，ＤＤｐ２）を取得するパラメータ取得ブロック（Ｓ１０２，Ｓ１０３，Ｓ３１０１，Ｓ４１０１，Ｓ４１０３，Ｓ５１０１，Ｓ５１０３）と、
前記パラメータ取得ブロックにより取得された前記性状パラメータに基づき、前記燃料を構成する成分の平均元素数（ＮＣ，ＮＨ）を算出する平均算出ブロック（Ｓ１０４，Ｓ２１０４，Ｓ３１０４，Ｓ４１０４，Ｓ５１０４）と、
前記平均算出ブロックにより算出された前記平均元素数に対応する前記燃料の曇点（ＴＭｃ）を、推定する曇点推定ブロック（Ｓ１０５）と、
前記曇点推定ブロックにより推定された前記曇点以下となる前記燃料の温度（ＴＭｊ）において、前記ヒータを制御することにより前記燃料フィルタを加熱するヒータ制御ブロック（Ｓ１０７，Ｓ１０８）とを、備える燃料系制御装置。
【請求項２】
前記燃料の現在温度（ＴＭｐ，ＴＭｐ１，ＴＭｐ２）を取得する温度取得ブロック（Ｓ１０１，Ｓ４１００，Ｓ４１０２）を、さらに備え、
前記パラメータ取得ブロック（Ｓ１０２，Ｓ１０３，Ｓ４１０１，Ｓ４１０３，Ｓ５１０１，Ｓ５１０３）は、前記温度取得ブロックにより取得された前記現在温度にて前記性状パラメータ（ＫＶｐ，ＤＤｐ，ＫＶｐ１，ＫＶｐ２，ＤＤｐ１，ＤＤｐ２）を取得し、
前記平均算出ブロック（Ｓ１０４，Ｓ４１０４，Ｓ５１０４）は、前記平均元素数に対する前記性状パラメータの相関を表すデータとして前記燃料の温度に応じて設定されたモデルデータに従って、前記パラメータ取得ブロックにより取得された前記性状パラメータと、前記温度取得ブロックにより取得された前記現在温度とに基づき、前記平均元素数を算出する請求項１に記載の燃料系制御装置。
【請求項３】
前記燃料の現在温度（ＴＭｐ）を取得する温度取得ブロック（Ｓ１０１）と、
前記パラメータ取得ブロック（Ｓ１０２，Ｓ１０３）により取得された前記性状パラメータ（ＫＶｐ，ＤＤｐ）を、前記温度取得ブロックにより取得された前記現在温度での値から、基準温度（ＴＭｂ）での値（ＫＶｅ，ＤＤｅ）へ変換するパラメータ変換ブロック（Ｓ２１０２，Ｓ２１０３）とを、さらに備え、
前記平均算出ブロック（Ｓ２１０４）は、前記平均元素数に対する前記性状パラメータの相関を表すデータとして設定された前記基準温度でのモデルデータに従って、前記温度取得ブロックによる取得値から前記パラメータ変換ブロックにより変換された前記性状パラメータに基づき、前記平均元素数を算出する請求項１に記載の燃料系制御装置。
【請求項４】
前記パラメータ取得ブロック（Ｓ１０２，Ｓ１０３，Ｓ３１０１）は、基準温度（ＴＭｂ）にて前記性状パラメータ（ＫＶｐ，ＤＤｐ）を取得し、
前記平均算出ブロック（Ｓ３１０４）は、前記平均元素数に対する前記性状パラメータの相関を表すデータとして設定された前記基準温度でのモデルデータに従って、前記パラメータ取得ブロックにより取得された前記性状パラメータとに基づき、前記平均元素数を算出する請求項１に記載の燃料系制御装置。
【請求項５】
前記平均算出ブロック（Ｓ１０４，Ｓ２１０４，Ｓ３１０４，Ｓ４１０４，Ｓ５１０４）は、前記燃料を構成する炭素の前記平均元素数である平均炭素数（ＮＣ）と、前記燃料を構成する水素の前記平均元素数である平均水素数（ＮＨ）とを、算出し、
前記曇点推定ブロック（Ｓ１０５）は、前記平均算出ブロックにより算出された前記平均炭素数及び前記平均水素数に基づき、前記曇点を算出する請求項１〜４のいずれか一項に記載の燃料系制御装置。
【請求項６】
前記パラメータ取得ブロック（Ｓ１０２，Ｓ１０３，Ｓ３１０１）は、前記燃料の同一温度（ＴＭｐ）での動粘度（ＫＶｐ）及び密度（ＤＤｐ）を、前記性状パラメータとして取得する請求項１〜５のいずれか一項に記載の燃料系制御装置。
【請求項７】
前記パラメータ取得ブロック（Ｓ４１０１，Ｓ４１０３）は、前記燃料の複数温度（ＴＭｐ１，ＴＭｐ２）での動粘度（ＫＶｐ１，ＫＶｐ２）を、前記性状パラメータとして取得する請求項１〜５のいずれか一項に記載の燃料系制御装置。
【請求項８】
前記パラメータ取得ブロック（Ｓ５１０１，Ｓ５１０３）は、前記燃料の複数温度（ＴＭｐ１，ＴＭｐ２）での密度（ＤＤｐ１，ＤＤｐ２）を、前記性状パラメータとして取得する請求項１〜５のいずれか一項に記載の燃料系制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、燃料系制御装置に関する。
【背景技術】
【０００２】
従来、内燃機関の燃料系において燃焼室への燃料を濾過する燃料フィルタのヒータによる加熱を制御することで、高い濾過効率を継続的に確保可能にする燃料系制御装置は、広く知られている。
【０００３】
こうした燃料系制御装置の一種として特許文献１の開示装置では、燃料の動粘度に対する曇点の相関データに従って、燃料の動粘度及び温度に基づくことで、ヒータ制御による燃料フィルタの加熱が実行されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１４−５１９２０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかし、本発明者らが鋭意研究を行った結果、ある単一温度において動粘度が同じ燃料でも、構成成分の違いに応じてワックス化による析出結晶の成長度にも違いが生じ、曇点は異なることが分かっている。そのために特許文献１の開示装置では、動粘度に対する曇点の相関データに基づいても、燃料の構成成分によっては現在温度が曇点以下となったにも拘わらず加熱が実行されずに、燃料フィルタ内の燃料がワックス化して最悪は目詰まりを招くという懸念があった。ここで燃料が一旦ワックス化すると、燃料フィルタに固着した析出結晶は加熱をしても溶融し難くなることから、高い濾過効率を継続的に確保する上での対策が必要であった。
【０００６】
本発明は、以上説明した問題に鑑みてなされたものであって、その目的は、燃料フィルタの高い濾過効率を継続的に確保する燃料系制御装置を提供することにある。
【課題を解決するための手段】
【０００７】
以下、課題を達成するための発明の技術的手段について、説明する。尚、発明の技術的手段を開示する特許請求の範囲及び本欄に記載された括弧内の符号は、後に詳述する実施形態に記載された具体的手段との対応関係を示すものであり、発明の技術的範囲を限定するものではない。
【０００８】
上述の課題を解決するために開示された発明は、
内燃機関（１）の燃料系（１０）において燃焼室（１ａ）へ供給される燃料を濾過する燃料フィルタ（１１）のヒータ（１１０）による加熱を制御する燃料系制御装置（５０）であって、
燃料の性状を表す性状パラメータ（ＫＶｐ，ＤＤｐ，ＫＶｐ１，ＫＶｐ２，ＤＤｐ１，ＤＤｐ２）を取得するパラメータ取得ブロック（Ｓ１０２，Ｓ１０３，Ｓ３１０１，Ｓ４１０１，Ｓ４１０３，Ｓ５１０１，Ｓ５１０３）と、
パラメータ取得ブロックにより取得された性状パラメータに基づき、燃料を構成する成分の平均元素数（ＮＣ，ＮＨ）を算出する平均算出ブロック（Ｓ１０４，Ｓ２１０４，Ｓ３１０４，Ｓ４１０４，Ｓ５１０４）と、
平均算出ブロックにより算出された平均元素数に対応する燃料の曇点（ＴＭｃ）を、推定する曇点推定ブロック（Ｓ１０５）と、
曇点推定ブロックにより推定された曇点以下となる燃料の温度（ＴＭｊ）において、ヒータを制御することにより燃料フィルタを加熱するヒータ制御ブロック（Ｓ１０７，Ｓ１０８）とを、備える。
【０００９】
このような開示発明では、燃料の性状パラメータに基づき算出される燃料構成成分の平均元素数に対応させて、燃料の曇点が推定される。これによれば、燃料構成成分が異なっても曇点を推定できるので、当該曇点以下の燃料温度では燃料フィルタをヒータの制御により加熱して、目詰まりに先立つワックス化から抑止し得る。したがって、燃料フィルタの濾過効率を継続的に確保することが可能となる。
【図面の簡単な説明】
【００１０】
【図１】第一実施形態の適用される内燃機関の燃料系を示すブロック図である。
【図２】第一実施形態の加熱制御フローを示すフローチャートである。
【図３】第一実施形態の平均炭素数の算出について説明するためのグラフである。
【図４】第一実施形態の曇点の推定について説明するためのグラフである。
【図５】第二実施形態の加熱制御フローを示すフローチャートである。
【図６】第二実施形態のモデルデータ例を示すグラフである。
【図７】第二実施形態のモデルデータ例を示すグラフである。
【図８】第三実施形態の加熱制御フローを示すフローチャートである。
【図９】第四実施形態の加熱制御フローを示すフローチャートである。
【図１０】第五実施形態の加熱制御フローを示すフローチャートである。
【発明を実施するための形態】
【００１１】
以下、本発明の複数の実施形態を図面に基づいて説明する。尚、各実施形態において対応する構成要素には同一の符号を付すことにより、重複する説明を省略する場合がある。各実施形態において構成の一部分のみを説明している場合、当該構成の他の部分については、先行して説明した他の実施形態の構成を適用することができる。また、各実施形態の説明において明示している構成の組み合わせばかりではなく、特に組み合わせに支障が生じなければ、明示していなくても複数の実施形態の構成同士を部分的に組み合わせることができる。
【００１２】
（第一実施形態）
図１に示すように本発明の第一実施形態による「燃料系制御装置」としてのＥＣＵ（Electronic Control Unit）５０は、車両に搭載される内燃機関１の燃料系１０に適用される。ＥＣＵ５０は、内燃機関１の燃焼室１ａへ供給される燃料を濾過する燃料フィルタ１１のヒータ１１０による加熱を、制御する。
【００１３】
具体的に燃料系１０は、燃料タンク１２、燃料フィルタ１１、高圧燃料ポンプ１３、コモンレール１４、温度センサ１５、動粘度センサ１６及び密度センサ１７等を、含んでいる。燃料タンク１２は、燃焼室１ａへ供給燃料を貯留する。
【００１４】
燃料フィルタ１１は、例えば濾紙等のフィルタエレメント１１１を容器１１２に収容して、構成されている。フィルタエレメント１１１は、燃料タンク１２から燃焼室１ａへの供給燃料を、供給経路の中途部にて通過させる。このときフィルタエレメント１１１は、通過する燃料中の異物を捕集することで、当該通過燃料を濾過する。
【００１５】
燃料フィルタ１１の容器１１２には、ヒータ１１０が設置されている。ヒータ１１０は、例えばＰＴＣ（Positive Temperature Coefficient）素子等の発熱素子を主体として、構成されている。ヒータ１１０は、ＥＣＵ５０からの通電調整により、フィルタエレメント１１１の加熱制御を受ける。これによりフィルタエレメント１１１では、燃料を構成する燃料構成成分に応じた曇点にて目詰まりに先立ち生じ始めるワックス化に対して、抑止作用の発揮が可能となっている。
【００１６】
高圧燃料ポンプ１３は、燃料タンク１２から燃料フィルタ１１を介して吸入した燃料を、燃焼室１ａへの供給燃料としてコモンレール１４に圧送する。そこで高圧燃料ポンプ１３では、ＥＣＵ５０からの通電調整により、コモンレール１４への圧送燃料の圧力制御を受ける。このとき、コモンレール１４への圧送燃料の圧力は、例えばアクセル開度及びエンジン回転数等といった内燃機関１の運転状態に基づくことで、ＥＣＵ５０により制御される。
【００１７】
コモンレール１４は、内燃機関１において複数の燃焼室１ａ毎に燃料噴射弁１ｂから噴射させる供給燃料を、内部に一時的に蓄える。ここで、燃焼室１ａ毎の各燃料噴射弁１ｂの開閉弁時期は、ＥＣＵ５０により制御される。これにより各燃料噴射弁１ｂの開弁時には、対応する燃焼室１ａへコモンレール１４からの供給燃料が噴射されることで、当該対応燃焼室１ａへ吸入された空気と混合される。その結果、燃料と空気との混合気が内燃機関１のピストンにより圧縮されることで、当該燃料が自着火して燃焼する。即ち、本実施形態の内燃機関１は、圧縮自着火式のディーゼルエンジンであり、当該内燃機関１への供給燃料として、ディーゼル燃料（軽油）が使用されている。ここで、そうした燃料の構成成分は、世界の地域毎（国毎）に相異なっている。
【００１８】
コモンレール１４には、内部燃料の圧力が耐圧値を超えるのを規制するために、減圧弁１４０が設置されている。減圧弁１４０の開弁時にコモンレール１４から逃がされる燃料は、図１に示すように還流先として燃料タンク１２及び燃料フィルタ１１の一方が還流弁１４１により選択されるようになっていてもよいし、図示はしないが還流先が燃料タンク１２に限定されていてもよい。尚、還流弁１４１による還流先の選択は、燃料温度に応じた感温素子の働きにより制御される。また。還流弁１４１による還流先を選択される燃料は、コモンレール１４から逃がされる燃料だけでなく、図１に示すように、高圧燃料ポンプ１３や各燃料噴射弁１ｂから逃がされる燃料等であってもよい。
【００１９】
各センサ１５，１６，１７は、燃料系１０のうち燃料タンク１２から燃焼室１ａ毎の各燃料噴射弁１ｂに至る燃料経路において、それぞれ所定の一箇所（図１に明示はされていない箇所）に設置されている。温度センサ１５は、例えばサーミスタ等のセンサ素子を主体として、構成されている。温度センサ１５は、燃料タンク１２から燃焼室１ａへの供給燃料の温度を、設置箇所において検出する。動粘度センサ１６は、例えば音叉式、超音波式若しくは細管式等の直接計測型、又は密度等からの間接推定型のセンサ素子を主体として、構成されている。動粘度センサ１６は、燃料タンク１２から燃焼室１ａへの供給燃料の性状を表す性状パラメータとしての動粘度を、設置箇所において検出する。密度センサ１７は、例えば固有振動周期測定式等のセンサ素子を主体として、構成されている。密度センサ１７は、燃料タンク１２から燃焼室１ａへの供給燃料の性状パラメータとしての密度を、設置箇所において検出する。ここで特に、密度センサ１７による密度の検出と動粘度センサ１６による動粘度の検出とは、燃料温度が実質同一と擬制可能な設置箇所にて、本実施形態では実行されている。
【００２０】
ＥＣＵ５０は、プロセッサ５０ａ及びメモリ５０ｂを有するマイクロコンピュータを主体として、構成されている。ＥＣＵ５０は、ヒータ１１０、高圧燃料ポンプ１３、各センサ１５，１６，１７及び各弁１ｂ，１４０の他、車両の各種センサに対し、直接的又は車内ネットワークを介して間接的に接続されている。ＥＣＵ５０は、各センサ１５，１６，１７の検出情報と内燃機関１の運転状態とを含む車両情報に基づき、ヒータ１１０、高圧燃料ポンプ１３及び各弁１ｂ，１４０の作動を制御する。
【００２１】
ここで詳細にＥＣＵ５０は、メモリ５０ｂに記憶された加熱制御プログラムをプロセッサ５０ａにより実行することで、図２に示す加熱制御フローの各ステップを機能的に実現する。尚、加熱制御フローは、車両において乗員からの内燃機関１の始動指令を与えるパワースイッチのオン操作に応じて開始される一方、乗員からの内燃機関１の停止指令を与える同スイッチのオフ操作に応じて終了する。また、加熱制御フロー中の「Ｓ」とは、各ステップを意味する。さらにまた、加熱制御プログラムを記憶するＥＣＵ５０のメモリ５０ｂは、例えば半導体メモリ、磁気媒体若しくは光学媒体等といった記憶媒体を、一つ又は複数使用して構成される。
【００２２】
まずＳ１０１では、性状パラメータとしての動粘度及び密度の検出箇所、即ち動粘度センサ１６及び密度センサ１７の設置箇所における燃料の現在温度ＴＭｐを、温度センサ１５による検出情報に基づき取得する。このとき、動粘度センサ１６及び密度センサ１７の設置箇所又は近傍箇所にて燃料温度を検出する温度センサ１５であれば、温度センサ１５による検出温度が直接的に現在温度ＴＭｐとして取得される。一方、動粘度センサ１６及び密度センサ１７の設置箇所からの離間箇所にて燃料温度を検出する温度センサ１５であれば、温度センサ１５による検出温度から間接的且つ推定的に現在温度ＴＭｐが取得される。
【００２３】
次のＳ１０２では、Ｓ１０１により取得された現在温度ＴＭｐと実質同一温度下にて、動粘度センサ１６の設置箇所における燃料の現在動粘度ＫＶｐを、動粘度センサ１６の検出情報に基づき取得する。続くＳ１０３では、Ｓ１０１により取得された現在温度ＴＭｐと実質同一温度下にて、密度センサ１７の設置箇所における燃料の現在密度ＤＤｐを、密度センサ１７の検出情報に基づき取得する。即ち、現在密度ＤＤｐと現在動粘度ＫＶｐとが同一温度ＴＭｐにて取得される。
【００２４】
さらに続くＳ１０４では、炭素及び水素を主に含んだ燃料構成成分の平均元素数として、平均炭素数ＮＣ及び平均水素数ＮＨを算出する。本発明者の鋭意研究結果によれば、平均炭素数ＮＣ及び平均水素数ＮＨに対して現在動粘度ＫＶｐ及び現在密度ＤＤｐは、次の式１，２の相関を有することが判明している。
ＮＣ＝αｔ１・ＫＶｐ＋βｔ１・ＮＨ＋γｔ１  …（式１）
ＮＣ＝αｔ２・ＤＤｐ＋βｔ２・ＮＨ＋γｔ２  …（式２）
【００２５】
ここで、式１の係数αｔ１，βｔ１，γｔ１と、式２の係数αｔ２，βｔ２，γｔ２とは、例えば設計時の実験結果又はシミュレーション結果等に基づくことで、予め任意の燃料温度に応じて設定される。そこで第一実施形態では、平均炭素数ＮＣ及び平均水素数ＮＨに対する現在動粘度ＫＶｐ及び現在密度ＤＤｐの相関を表したモデルデータとして、燃料温度毎に予め設定の係数αｔ１，βｔ１，γｔ１及び係数αｔ２，βｔ２，γｔ２を導入した式１，２に、Ｓ１０４の算出が従うこととなる。具体的にＳ１０４では、燃料温度毎の係数αｔ１，βｔ１，γｔ１及び係数αｔ２，βｔ２，γｔ２の代入された式１，２の組のうち、Ｓ１０１により取得された現在温度ＴＭｐに対応する式１，２の組を、選択する。続いてＳ１０４では、Ｓ１０２，Ｓ１０３のそれぞれにより取得された現在動粘度ＫＶｐ及び現在密度ＤＤｐを、選択した式１，２の組に導入する。さらに続いてＳ１０４では、現在動粘度ＫＶｐ及び現在密度ＤＤｐを導入した式１，２の組の連立方程式を解くことで、平均炭素数ＮＣ及び平均水素数ＮＨを算出する。
【００２６】
以上の説明からＳ１０４では、燃料温度に応じたモデルデータである式１，２に従って、現在動粘度ＫＶｐ及び現在密度ＤＤｐと現在温度ＴＭｐとに基づく平均炭素数ＮＣ及び平均水素数ＮＨを算出するのである。図３は、燃料構成成分に応じた実際の平均炭素数ＮＣに対し、こうして算出された平均炭素数ＮＣが実質一致することを、証明している。尚、メモリ５０ｂの記憶容量等の観点から、式１，２に予め対応付けられる燃料温度として離散値が採用される場合に、そうした離散値間の現在温度ＴＭｐでは、補間演算された係数αｔ１，βｔ１，γｔ１及び係数αｔ２，βｔ２，γｔ２が式１，２に導入される。
【００２７】
図２に示すように続くＳ１０５では、燃料構成成分に応じた曇点ＴＭｃを推定する。本発明者の鋭意研究結果によれば、平均炭素数ＮＣ及び平均水素数ＮＨに対して曇点ＴＭｃは、次の式３の相関を有することから、推定可能であることが判明している。
ＴＭｃ＝α３・ＮＣ＋β３・ＮＨ＋γ３  …（式３）
【００２８】
ここで、式３の係数α３，β３，γ３は、例えば設計時の実験結果又はシミュレーション結果等に基づくことで、予め設定される。そこで第一実施形態では、平均炭素数ＮＣ及び平均水素数ＮＨに対応する曇点ＴＭｃを表したモデルデータとして、予め設定の係数α３，β３，γ３を導入した式３に、Ｓ１０５の推定が従うこととなる。具体的にＳ１０５では、Ｓ１０４により算出された平均炭素数ＮＣ及び平均水素数ＮＨを、係数α３，β３，γ３の代入された式３に導入する。こうして平均炭素数ＮＣ及び平均水素数ＮＨを導入の式３を解くことで、それらの数ＮＣ，ＮＨに対応する曇点ＴＭｃを算出する。図４は、燃料構成成分に応じた実際の曇点ＴＭｃに対し、算出された曇点ＴＭｃが近似的に推定可能であることを、証明している。
【００２９】
図２に示すように続くＳ１０６では、燃料フィルタ１１のフィルタエレメント１１１を現在通過する燃料の温度として、判定時温度ＴＭｊを温度センサ１５による検出情報に基づき取得する。このとき、フィルタエレメント１１１からの離間箇所にて燃料温度を検出する温度センサ１５であれば、温度センサ１５による検出温度から間接的且つ推定的に判定時温度ＴＭｊが取得される。一方、フィルタエレメント１１１の近傍箇所にて燃料温度を検出する温度センサ１５であれば、温度センサ１５による検出温度が直接的に判定時温度ＴＭｊとして取得される。尚、Ｓ１０５から直接移行する初回のＳ１０６では、Ｓ１０１での温度センサ１５による検出情報に基づくことで、判定時温度ＴＭｊを取得してもよい。
【００３０】
さらに続くＳ１０７では、Ｓ１０６により取得された判定時温度ＴＭｊと、Ｓ１０５により推定された曇点ＴＭｃとを比較することで、それら値ＴＭｊ，ＴＭｃの大小関係を判定する。その結果、判定時温度ＴＭｊが曇点ＴＭｃ以下であると判定された場合には、Ｓ１０８へ移行する一方、判定時温度ＴＭｊが曇点ＴＭｃ超過であると判定された場合には、Ｓ１０９へ移行する。
【００３１】
Ｓ１０８では、通電量の調整によりヒータ１１０の発熱を制御することで、燃料フィルタ１１のフィルタエレメント１１１を加熱する。このときには、フィルタエレメント１１１の加熱温度を、例えばＳ１０７にて比較された値ＴＭｊ，ＴＭｃの差分値に応じた可変温度に制御してもよいし、予め決められた固定温度に制御してもよい。これに対してＳ１０９では、通電量のカットによりヒータ１１０の発熱を停止させることで、燃料フィルタ１１のフィルタエレメント１１１に対する加熱を停止状態とする。尚、これらＳ１０８，Ｓ１０９のいずれの場合にも、実行終了後にはＳ１０６へと戻ることになる。
【００３２】
ここまで説明したように第一実施形態では、ＥＣＵ５０のうちＳ１０１を実行する機能部分が「温度取得ブロック」に相当し、ＥＣＵ５０のうちＳ１０２，Ｓ１０３を実行する機能部分が「パラメータ取得ブロック」に相当する。また第一実施形態では、ＥＣＵ５０のうちＳ１０４を実行する機能部分が「平均算出ブロック」に相当し、ＥＣＵ５０のうちＳ１０５を実行する機能部分が「曇点推定ブロック」に相当する。さらに第一実施形態では、ＥＣＵ５０のうちＳ１０７，Ｓ１０８を実行する機能部分が「ヒータ制御ブロック」に相当する。
【００３３】
（作用効果）
以上説明した第一実施形態の作用効果を、以下に説明する。
【００３４】
第一実施形態では、燃料の性状パラメータに基づき算出される燃料構成成分の平均元素数に対応させて、燃料の曇点ＴＭｃが推定される。これによれば、燃料構成成分が異なっても曇点ＴＭｃを推定できるので、当該曇点ＴＭｃ以下の判定時温度ＴＭｊでは燃料フィルタ１１をヒータ１１０の制御により加熱して、目詰まりに先立つワックス化から抑止し得る。したがって、燃料フィルタ１１の濾過効率を継続的に確保することが可能となる。
【００３５】
しかも第一実施形態によると、温度依存性のある性状パラメータについては、平均元素数に対する相関を燃料温度に応じて表したモデルデータに従うことになるので、当該性状パラメータに基づいた平均元素数の算出値は、現在温度ＴＭｐにも基づくことで正確に求められ得る。これによれば、燃料構成成分が異なっても性状パラメータの温度特性まで加味した曇点ＴＭｃを推定できるので、当該曇点ＴＭｃ以下での燃料フィルタ１１の加熱により高濾過効率の継続的確保の信頼性を高めることが可能となる。
【００３６】
ここで第一実施形態によると、燃料の構成成分は主に炭素及び水素を含むことから、平均元素数としては平均炭素数ＮＣ及び平均水素数ＮＨを算出して曇点ＴＭｃの推定値に対応させることで、当該推定値の精度を高め得る。故に、推定精度の高い曇点ＴＭｃ以下での燃料フィルタ１１の加熱により高濾過効率の継続的確保の信頼性を高めることが可能となる。
【００３７】
さらに第一実施形態によると、平均元素数としての平均炭素数ＮＣ及び平均水素数ＮＨに対して感度の高い相関を示す性状パラメータのうち、同一温度ＴＭｐでの現在動粘度ＫＶｐ及び現在密度ＤＤｐに基づくことから、それら値ＮＣ，ＮＨの算出値を正確に求め得る。これによれば、燃料構成成分が異なっても平均炭素数ＮＣ及び平均水素数ＮＨの正確な算出値に対応させて曇点ＴＭｃを推定できるので、当該曇点ＴＭｃ以下での燃料フィルタ１１の加熱により高濾過効率の継続的確保の信頼性を高めることが可能となる。
【００３８】
（第二実施形態）
本発明の第二実施形態は、第一実施形態の変形例である。図５に示すように第二実施形態の加熱制御フローでは、Ｓ１０２に続くＳ２１０２と、Ｓ１０３に続くＳ２１０３と、Ｓ１０４に代わるＳ２１０４とを実行する。
【００３９】
Ｓ２１０２では、Ｓ１０２により取得された現在動粘度ＫＶｐを、Ｓ１０１により取得された現在温度ＴＭｐでの値から、基準温度ＴＭｂでの推定値である推定動粘度ＫＶｅへと変換する。本発明者の鋭意研究結果によれば、現在動粘度ＫＶｐに対して現在温度ＴＭｐ及び基準温度ＴＭｂ間での動粘度変化量ΔＫＶは、図６に例示するような相関を有することが判明している。
【００４０】
ここで、例えば設計時の実験結果又はシミュレーション結果等に基づくことで、現在動粘度ＫＶｐに対する動粘度変化量ΔＫＶの相関は、予め任意の燃料温度に応じて設定される。そこで第二実施形態では、現在動粘度ＫＶｐに対する動粘度変化量ΔＫＶの相関を燃料温度に応じて表したモデルデータに、Ｓ２１０２の変換が従うこととなる。具体的にＳ２１０２では、現在動粘度ＫＶｐに対する動粘度変化量ΔＫＶの燃料温度に応じたモデルデータのうち、Ｓ１０１により取得された現在温度ＴＭｐに対応するデータを、選択する。続いてＳ２１０２では、Ｓ１０２により取得された現在動粘度ＫＶｐを、選択したモデルデータに導入することで、現在温度ＴＭｐ及び基準温度ＴＭｂ間での動粘度変化量ΔＫＶを推定する。さらに続いてＳ２１０２では、Ｓ１０２により取得された現在動粘度ＫＶｐに対して、推定した動粘度変化量ΔＫＶを加算することで、現在温度ＴＭｐでの現在動粘度ＫＶｐから変換された基準温度ＴＭｂでの推定動粘度ＫＶｅを取得する。
【００４１】
尚、基準温度ＴＭｂは、後に詳述する式４，５のモデルデータに対応付けられる燃料温度として、予め設定される。また、図６に例示の如きモデルデータにおいて複数のドットデータは、現在動粘度ＫＶｐに対する動粘度変化量ΔＫＶの各地域（各国）の燃料構成成分に応じた実験値を、表している。さらにメモリ５０ｂの記憶容量等の観点から、図６に例示の如きモデルデータに予め対応付けられる燃料温度として離散値が採用される場合に、そうした離散値間の現在温度ＴＭｐでは、補間演算によりモデルデータが補正される。またさらに動粘度変化量ΔＫＶは、基準温度ＴＭｂでの動粘度から現在温度ＴＭｐでの動粘度を差し引くことで、正負の符号を有する値として定義される。これにより、基準温度ＴＭｂが現在温度ＴＭｐよりも高い場合に動粘度変化量ΔＫＶは、推定動粘度ＫＶｅを取得するための負の符号を図６に例示の如く有する値として、現在動粘度ＫＶｐに加算されることとなる。
【００４２】
図５に示すようにＳ２１０３では、Ｓ１０３により取得された現在密度ＤＤｐを、Ｓ１０１により取得された現在温度ＴＭｐでの値から、基準温度ＴＭｂでの推定値である推定密度ＤＤｅへと変換する。本発明者の鋭意研究結果によれば、現在密度ＤＤｐに対して基準温度ＴＭｂでの推定密度ＤＤｅは、図７に例示するような相関を有することが判明している。
【００４３】
ここで、例えば設計時の実験結果又はシミュレーション結果等に基づくことで、現在密度ＤＤｐに対する推定密度ＤＤｅの相関は、予め任意の燃料温度に応じて設定される。そこで第二実施形態では、現在密度ＤＤｐに対する推定密度ＤＤｅの相関を燃料温度に応じて表したモデルデータに、Ｓ２１０３の変換が従うこととなる。具体的にＳ２１０３では、現在密度ＤＤｐに対する推定密度ＤＤｅの燃料温度に応じたモデルデータのうち、Ｓ１０１により取得された現在温度ＴＭｐに対応するデータを、選択する。続いてＳ２１０３では、Ｓ１０３により取得された現在密度ＤＤｐを、選択したモデルデータに導入することで、現在温度ＴＭｐでの現在密度ＤＤｐから変換された基準温度ＴＭｂでの推定密度ＤＤｅを取得する。
【００４４】
尚、基準温度ＴＭｂは、Ｓ２１０２で想定される値と一致する。また、図７に例示の如きモデルデータにおいて複数のドットデータは、現在密度ＤＤｐに対する推定密度ＤＤｅの各地域（各国）の燃料構成成分に応じた実験値を、表している。またメモリ５０ｂの記憶容量等の観点から、図７に例示の如きモデルデータに予め対応付けられる燃料温度として離散値が採用される場合に、そうした離散値間の現在温度ＴＭｐでは、補間演算によりモデルデータが補正される。
【００４５】
図５に示すようにＳ２１０４では、基準温度ＴＭｂでのモデルデータである、次の式４，５に従って、推定動粘度ＫＶｅ及び推定密度ＤＤｅに基づく平均炭素数ＮＣ及び平均水素数ＮＨを算出する。
ＮＣ＝α１・ＫＶｅ＋β１・ＮＨ＋γ１  …（式４）
ＮＣ＝α２・ＤＤｅ＋β２・ＮＨ＋γ２  …（式５）
【００４６】
ここで、式４の係数α１，β１，γ１と、式５の係数α２，β２，γ２とは、例えば設計時の実験結果又はシミュレーション結果等に基づくことで、予め基準温度ＴＭｂに対応付けて設定される。そこで第二実施形態では、平均炭素数ＮＣ及び平均水素数ＮＨに対する推定動粘度ＫＶｅ及び推定密度ＤＤｅの相関を表したモデルデータとして、予め設定の基準温度ＴＭｂでの係数α１，β１，γ１及び係数α２，β２，γ２を導入した式４，５に、Ｓ２１０４の算出が従うこととなる。具体的にＳ２１０４では、Ｓ２１０２，Ｓ２１０３のそれぞれにより性状パラメータの基準温度ＴＭｂでの変換値として取得された推定動粘度ＫＶｅ及び推定密度ＤＤｅを、係数α１，β１，γ１及び係数α２，β２，γ２の代入された式４，５の組に、導入する。続いてＳ２１０４では、推定動粘度ＫＶｅ及び推定密度ＤＤｅを導入した式４，５の組の連立方程式を解くことで、平均炭素数ＮＣ及び平均水素数ＮＨを算出する。
【００４７】
このように第二実施形態では、ＥＣＵ５０のうちＳ２１０２，Ｓ２１０３を実行する機能部分が「パラメータ変換ブロック」に相当し、ＥＣＵ５０のうちＳ２１０４を実行する機能部分が「平均算出ブロック」に相当する。
【００４８】
以上説明した第二実施形態によると、温度依存性のある性状パラメータについては、平均元素数に対する相関を表した基準温度ＴＭｂでのモデルデータに従うことになるので、現在温度ＴＭｐでの値から基準温度ＴＭｂでの値へと変換された性状パラメータに基づく平均元素数の算出値は、正確に求められ得る。これによれば、燃料構成成分が異なっても性状パラメータの温度特性まで加味した曇点ＴＭｃを推定できるので、当該曇点ＴＭｃ以下での燃料フィルタ１１の加熱により高濾過効率の継続的確保の信頼性を高めることが可能となる。
【００４９】
（第三実施形態）
本発明の第三実施形態は、第一実施形態の変形例である。図８に示すように第三実施形態の加熱制御フローでは、Ｓ１０１に続くＳ３１０１と、Ｓ１０４に代わるＳ３１０４とを実行する。
【００５０】
Ｓ３１０１では、Ｓ１０１により取得された現在温度ＴＭｐは基準温度ＴＭｂと一致するか否かを、判定する。このとき、現在温度ＴＭｐが基準温度ＴＭｂと一致するとの高低判定は、完全一致する状況にだけ下してもよいし、所定の誤差範囲内にて一致する状況にて擬制的に下してもよい。こうしたＳ３１０１の結果、不一致により否定判定が下された場合にはＳ１０１へ戻る一方、一致により肯定判定が下された場合にはＳ１０２へ移行する。故にＳ１０２，Ｓ１０３では、現在温度ＴＭｐである基準温度ＴＭｂと実質同一温度下にて、現在動粘度ＫＶｐ及び現在密度ＤＤｐがそれぞれ取得されることとなる。
【００５１】
Ｓ３１０４では、基準温度ＴＭｂでのモデルデータである、次の式６，７に従って、現在動粘度ＫＶｐ及び現在密度ＤＤｐに基づく平均炭素数ＮＣ及び平均水素数ＮＨを算出する。
ＮＣ＝α１・ＫＶｐ＋β１・ＮＨ＋γ１  …（式６）
ＮＣ＝α２・ＤＤｐ＋β２・ＮＨ＋γ２  …（式７）
【００５２】
ここで、式６の係数α１，β１，γ１と、式７の係数α２，β２，γ２とは、例えば設計時の実験結果又はシミュレーション結果等に基づくことで、予め基準温度ＴＭｂに対応付けて設定される。そこで第三実施形態では、平均炭素数ＮＣ及び平均水素数ＮＨに対する現在動粘度ＫＶｐ及び現在密度ＤＤｐの相関を表したモデルデータとして、予め設定の基準温度ＴＭｂでの係数α１，β１，γ１及び係数α２，β２，γ２を導入した式６，７に、Ｓ３１０４の算出が従うこととなる。具体的にＳ３１０４では、Ｓ１０２，Ｓ１０３のそれぞれにより性状パラメータとして取得された基準温度ＴＭｂでの現在動粘度ＫＶｐ及び現在密度ＤＤｐを、係数α１，β１，γ１及び係数α２，β２，γ２の代入された式６，７の組に、導入する。続いてＳ３１０４では、現在動粘度ＫＶｐ及び現在密度ＤＤｐを導入した式４，５の組の連立方程式を解くことで、平均炭素数ＮＣ及び平均水素数ＮＨを算出する。
【００５３】
このような第三実施形態では、ＥＣＵ５０のうちＳ３１０１，Ｓ１０２，Ｓ１０３を実行する機能部分が「パラメータ取得ブロック」に相当し、ＥＣＵ５０のうちＳ３１０４を実行する機能部分が「平均算出ブロック」に相当する。
【００５４】
以上説明した第三実施形態によると、温度依存性のある性状パラメータについては、平均元素数に対する相関を表した基準温度ＴＭｂでのモデルデータに従うことから、基準温度ＴＭｂにて取得された性状パラメータに基づく平均元素数の算出値は、正確に求められ得る。これによれば、燃料構成成分が異なっても性状パラメータの温度特性まで加味した曇点ＴＭｃを推定できるので、当該曇点ＴＭｃ以下での燃料フィルタ１１の加熱により高濾過効率の継続的確保の信頼性を高めることが可能となる。
【００５５】
（第四実施形態）
本発明の第四実施形態は、第一実施形態の変形例である。図９に示すように第四実施形態の加熱制御フローでは、Ｓ１０１〜Ｓ１０４に代わるＳ４１００〜Ｓ４１０４を実行する。
【００５６】
Ｓ４１００では、燃料の第一現在温度ＴＭｐ１を、第一実施形態の現在温度ＴＭｐに準じて取得する。Ｓ４１０１では、Ｓ４１００により取得された第一現在温度ＴＭｐ１と実質同一温度下にて、燃料の第一現在動粘度ＫＶｐ１を、第一実施形態の現在動粘度ＫＶｐに準じて取得する。
【００５７】
Ｓ４１０２では、燃料の第二現在温度ＴＭｐ２として、Ｓ４１００により取得された第一現在温度ＴＭｐ１とは異なる温度を、第一実施形態の現在温度ＴＭｐに準じて取得する。Ｓ４１０３では、Ｓ４１０２により取得された第二現在温度ＴＭｐ２と実質同一温度下にて、燃料の第二現在動粘度ＫＶｐ２を、第一実施形態の現在動粘度ＫＶｐに準じて取得する。
【００５８】
Ｓ４１０４では、燃料温度に応じたモデルデータとして第一実施形態に説明の式１に従うことで、第一及び第二現在動粘度ＫＶｐ１，ＫＶｐ２に基づく平均炭素数ＮＣ及び平均水素数ＮＨを算出する。具体的にＳ４１０４では、燃料温度毎の係数αｔ１，βｔ１，γｔ１を導入した式１のうち、Ｓ４１００，Ｓ４１０２により取得された第一及び第二現在温度ＴＭｐ１，ＴＭｐ２にそれぞれ対応する式１の組を、選択する。続いてＳ４１０４では、Ｓ４１０１により取得された第一現在動粘度ＫＶｐ１を、第一現在温度ＴＭｐ１に対応して選択された式１の現在動粘度ＫＶｐに導入する。それと共にＳ４１０４では、Ｓ４１０３により取得された第二現在動粘度ＫＶｐ２を、第二現在温度ＴＭｐ２に対応して選択された式１の現在動粘度ＫＶｐに導入する。こうした後にＳ４１０４では、第一及び第二現在動粘度ＫＶｐ１，ＫＶｐ２をそれぞれ導入した式１の組の連立方程式を解くことで、平均炭素数ＮＣ及び平均水素数ＮＨを算出する。
【００５９】
このような第四実施形態では、ＥＣＵ５０のうちＳ４１００，Ｓ４１０２を実行する機能部分が「温度取得ブロック」に相当し、ＥＣＵ５０のうちＳ４１０１，Ｓ４１０３を実行する機能部分が「パラメータ取得ブロック」に相当する。また第三実施形態では、ＥＣＵ５０のうちＳ４１０４を実行する機能部分が「平均算出ブロック」に相当する。
【００６０】
以上説明した第四実施形態によると、平均元素数としての平均炭素数ＮＣ及び平均水素数ＮＨに対して感度の高い相関を示す性状パラメータのうち、複数温度ＴＭｐ１，ＴＭｐ２での現在動粘度ＫＶｐ１，ＫＶｐ２に基づくことで、それら値ＮＣ，ＮＨの算出値を正確に求め得る。これによれば、燃料構成成分が異なっても平均炭素数ＮＣ及び平均水素数ＮＨの正確な算出値に対応させて曇点ＴＭｃを推定できるので、当該曇点ＴＭｃ以下での燃料フィルタ１１の加熱により高濾過効率の継続的確保の信頼性を高めることが可能となる。尚、こうした第四実施形態では、作用効果を発揮する上において、密度センサ１７の検出情報又は同センサ１７自体が不要となる。
【００６１】
（第五実施形態）
本発明の第五実施形態は、第四実施形態の変形例である。図１０に示すように第五実施形態の加熱制御フローでは、Ｓ４１０１，Ｓ４１０３，Ｓ４１０４にそれぞれ代わるＳ５１０１，Ｓ５１０３，Ｓ５１０４を実行する。
【００６２】
Ｓ５１０１では、Ｓ４１００により取得された第一現在温度ＴＭｐ１と実質同一温度下にて、燃料の第一現在密度ＤＤｐ１を、第一実施形態の現在密度ＤＤｐに準じて取得する。Ｓ５１０３では、Ｓ４１０１により取得された第二現在温度ＴＭｐ２と実質同一温度下にて、燃料の第二現在密度ＤＤｐ２を、第一実施形態の現在密度ＤＤｐに準じて取得する。
【００６３】
Ｓ５１０４では、燃料温度に応じたモデルデータとして第一実施形態に説明の式２に従うことで、第一及び第二現在密度ＤＤｐ１，ＤＤｐ２に基づく平均炭素数ＮＣ及び平均水素数ＮＨを算出する。具体的にＳ５１０４では、燃料温度毎の係数αｔ２，βｔ２，γｔ２を導入した式２のうち、Ｓ４１００，Ｓ４１０２により取得された第一及び第二現在温度ＴＭｐ１，ＴＭｐ２にそれぞれ対応する式２の組を、選択する。続いてＳ５１０４では、Ｓ５１０１により取得された第一現在密度ＤＤｐ１を、第一現在温度ＴＭｐ１に対応して選択された式２の現在密度ＤＤｐに導入する。それと共にＳ５１０４では、Ｓ５１０３により取得された第二現在密度ＤＤｐ２を、第二現在温度ＴＭｐ２に対応して選択された式２の現在密度ＤＤｐに導入する。こうした後にＳ５１０４では、第一及び第二現在密度ＤＤｐ１，ＤＤｐ２をそれぞれ導入した式２の組の連立方程式を解くことで、平均炭素数ＮＣ及び平均水素数ＮＨを算出する。
【００６４】
このような第五実施形態では、ＥＣＵ５０のうちＳ５１０１，Ｓ５１０３を実行する機能部分が「パラメータ取得ブロック」に相当し、ＥＣＵ５０のうちＳ５１０４を実行する機能部分が「平均算出ブロック」に相当する。
【００６５】
以上説明した第五実施形態によると、平均元素数としての平均炭素数ＮＣ及び平均水素数ＮＨに対して感度の高い相関を示す性状パラメータのうち、複数温度ＴＭｐ１，ＴＭｐ２での現在密度ＤＤｐ１，ＤＤｐ２に基づくことで、それら値ＮＣ，ＮＨの算出値を正確に求め得る。これによれば、燃料構成成分が異なっても平均炭素数ＮＣ及び平均水素数ＮＨの正確な算出値に対応させて曇点ＴＭｃを推定できるので、当該曇点ＴＭｃ以下での燃料フィルタ１１の加熱により高濾過効率の継続的確保の信頼性を高めることが可能となる。尚、こうした第五実施形態では、作用効果を発揮する上において、動粘度センサ１６の検出情報又は同センサ１６自体が不要となる。
【００６６】
（他の実施形態）
以上、本発明の複数の実施形態について説明したが、本発明は、それらの実施形態に限定して解釈されるものではなく、本発明の要旨を逸脱しない範囲内において種々の実施形態及び組み合わせに適用することができる。
【００６７】
具体的に、第一〜第五実施形態に関する変形例１では、内燃機関１の運転条件に基づく推定により、現在温度ＴＭｐ，ＴＭｐ１，ＴＭｐ２及び判定時温度ＴＭｊが取得されてもよい。この変形例１では、作用効果を発揮する上において、温度センサ１５の検出情報又は同センサ１５自体が不要となる。また変形例１では、現在温度ＴＭｐ，ＴＭｐ１，ＴＭｐ２及び判定時温度ＴＭｊの推定に利用される運転条件として、例えば内燃機関１の完全暖機後における燃料の低圧条件や中圧条件等といった使用頻度の大きい条件が、採用される。
【００６８】
第四実施形態に関する変形例２では、燃料系１０の燃料経路において動粘度センサ１６とは異なる箇所に実質同種の別の動粘度センサが設置されることで、第二現在動粘度ＫＶｐ２が当該別の動粘度センサの検出情報に基づき取得されてもよい。この変形例２では、別の動粘度センサの設置箇所における第二現在温度ＴＭｐ２が、第一実施形態の現在温度ＴＭｐに準じて取得されることとなる。
【００６９】
第四実施形態に関する変形例３では、第二実施形態に準じて第一及び第二現在温度ＴＭｐ１，ＴＭｐ２での第一及び現在動粘度ＫＶｐ１，ＫＶｐ２からそれぞれ変換した第一及び第二基準温度での値が、それら各基準温度に対応する式１の組に個別に導入されることで、平均炭素数ＮＣ及び平均水素数ＮＨが算出されてもよい。第四実施形態に関する変形例４では、第三実施形態に準じて、第一及び第二基準温度にそれぞれ一致した第一及び第二現在温度ＴＭｐ１，ＴＭｐ２での第一及び現在動粘度ＫＶｐ１，ＫＶｐ２が、それら各基準温度に対応する式１の組に個別に導入されることで、平均炭素数ＮＣ及び平均水素数ＮＨが算出されてもよい。
【００７０】
第五実施形態に関する変形例５では、燃料系１０の燃料経路において密度センサ１７とは異なる箇所に実質同種の密度センサが別に設置されることで、第二現在密度ＤＤｐ２が当該別の密度センサの検出情報に基づき取得されてもよい。この変形例５では、別の密度センサの設置箇所における第二現在温度ＴＭｐ２が、第一実施形態の現在温度ＴＭｐに準じて取得されることとなる。
【００７１】
第五実施形態に関する変形例６では、第二実施形態に準じて第一及び第二現在温度ＴＭｐ１，ＴＭｐ２での第一及び現在密度ＤＤｐ１，ＤＤｐ２からそれぞれ変換した第一及び第二基準温度での値が、それら各基準温度に対応する式２の組に個別に導入されることで、平均炭素数ＮＣ及び平均水素数ＮＨが算出されてもよい。第五実施形態に関する変形例７では、第三実施形態に準じて第一及び第二基準温度にそれぞれ一致した第一及び第二現在温度ＴＭｐ１，ＴＭｐ２での第一及び現在密度ＤＤｐ１，ＤＤｐ２が、それら各基準温度に対応する式２の組に個別に導入されることで、平均炭素数ＮＣ及び平均水素数ＮＨが算出されてもよい。
【００７２】
第一〜第五実施形態に関する変形例８では、加熱制御フローの少なくとも一部のステップが、ＥＣＵ５０により機能的に実現される代わりに、一つ又は複数のＩＣ等によりハードウェア的に実現されてもよい。
【符号の説明】
【００７３】
１  内燃機関、１ａ  燃焼室、１０  燃料系、１１  燃料フィルタ、１５  温度センサ、１６  動粘度センサ、１７  密度センサ、５０  ＥＣＵ、５０ａ  プロセッサ、５０ｂ  メモリ、１１０  ヒータ、１１１  フィルタエレメント、ＤＤｐ  現在密度、ＤＤｐ１  第一現在密度、ＤＤｐ２  第二現在密度、ＤＤｅ  推定密度、ＫＶｐ  現在動粘度、ＫＶｐ１  第一現在動粘度、ＫＶｐ２  第二現在動粘度、ＫＶｅ  推定動粘度、ΔＫＶ  動粘度変化量、ＮＣ  平均炭素数、ＮＨ  平均水素数、ＴＭｃ  曇点、ＴＭｊ  判定時温度、ＴＭｐ  現在温度、ＴＭｐ１  第一現在温度、ＴＭｐ２  第二現在温度、ＴＭｂ  基準温度、
【図１】
2018200024.tif 000003
【図２】
2018200024.tif 000004
【図３】
2018200024.tif 000005
【図４】
2018200024.tif 000006
【図５】
2018200024.tif 000007
【図６】
2018200024.tif 000008
【図７】
2018200024.tif 000009
【図８】
2018200024.tif 000010
【図９】
2018200024.tif 000011
【図１０】
2018200024.tif 000012