JP2004239129A - エンジン性能の予測解析方法、予測解析システム及びその制御プログラム - Google Patents
エンジン性能の予測解析方法、予測解析システム及びその制御プログラム Download PDFInfo
- Publication number
- JP2004239129A JP2004239129A JP2003027834A JP2003027834A JP2004239129A JP 2004239129 A JP2004239129 A JP 2004239129A JP 2003027834 A JP2003027834 A JP 2003027834A JP 2003027834 A JP2003027834 A JP 2003027834A JP 2004239129 A JP2004239129 A JP 2004239129A
- Authority
- JP
- Japan
- Prior art keywords
- working gas
- state
- calculation
- engine
- combustion chamber
- 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.)
- Granted
Links
Images
Landscapes
- Testing Of Engines (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
【解決手段】4サイクルエンジンのシミュレーションにおいて、吸排気の流れをCFDにより模擬する一方、化学反応SIMにより燃焼等を模擬する。気筒の吸気行程から圧縮行程への移行を模擬するときにはCFDを吸気弁の閉時期IVCまで行って、燃焼室に充填される吸気の状態(p,ρ,u,T)を正確に求める。その吸気の状態に基づいて、圧縮行程下死点BDCで燃焼室に充填されていると仮定したときの前記吸気の状態を推定し、この推定結果に対応するガス成分のデータを化学反応DB12から読み込んで、圧縮行程下死点BDCから膨張行程下死点まで化学反応SIMを行うようにする。
【選択図】 図4
Description
【発明の属する技術分野】
本発明は、少なくともCFD(Computational Fluid Dynamics:数値計算流体力学)の適用によりエンジンの作動ガスの状態を模擬して、該エンジンの性能を予測するための予測解析方法、予測解析システム及びその制御プログラムに関する。
【0002】
【従来の技術】
従来より、エンジンやトランスミッション等の性能を評価するために例えば特許文献1に開示されるような種々の計測・試験方法が提案されている。また、特許文献2には、エンジンの開発完了を待たずにパワートレインの性能を評価することのできるシミュレーションシステムが開示されている。
【0003】
そのようなシミュレーションの技術として、作動ガスである吸気や排気の運動をCFDの適用により解析し、この解析結果に基づいてエンジンの性能を予測することが一般的に行われている。すなわち、例えば吸気ポートから燃焼室へ吸い込まれる吸気の複雑な流れをコンピュータを用いた数値計算によって模擬する仮想の実験(シミュレーション)を行い、この仮想実験の結果に基づいて例えば吸気ポートの形状を決定することにより、試作や実験の繰り返しに費やされる開発工数を削減して、効率の良い設計・開発を行うことができる。
【0004】
また、気筒内の燃焼室における混合気の燃焼については、その混合気に含まれる多くのガス成分の反応をそれぞれ化学反応式によって記述するようにしたシミュレーションプログラムが存在する。これは、燃焼室の作動ガスを主に燃料として供給される炭化水素や空気中の窒素、酸素等により模擬するとともに、圧縮及び膨張行程における燃焼室の容積変化を容器モデルにより模擬して、その容器中での前記ガス成分や中間生成物の反応を逐次、記述することによって、燃焼の状態を模擬するものである。
【0005】
【特許文献1】
特表2002−526762号公報
【特許文献2】
特開2002−148147号公報
【0006】
【発明が解決しようとする課題】
しかしながら、エンジンの吸排気の流れを模擬するCFD演算と、燃焼状態を模擬する化学反応シミュレーションとを適切にリンクしてダイナミックに解くことは容易ではない。例えば4サイクルエンジンの吸気行程から圧縮行程への移行時には、まず、CFD演算の結果から燃焼室に充填される吸気の状態を求めて、その吸気の状態に対応するガス成分を用いて化学反応シミュレーションを行うものであるが、図11に模式的に示すように、通常は吸気弁の閉じる時期IVCは圧縮行程下死点BDCよりも遅角側に設定されており、燃焼室からの吸気の吹き返し等を考慮すればCFD演算は吸気弁が閉じるまで継続して行う必要がある。そして、その後、図に破線の矢印で示すように、演算結果のデータを化学反応シミュレーションプログラムに提供することになる。
【0007】
ところが、近年では吸気弁の開閉時期を運転状態等に応じて変更することのできる可変動弁機構(以下、VVTという)を備えたエンジンが普及しつつあり、このようなVVT付きのエンジンの場合は運転途中で吸気弁の開閉時期が変化することになるから、前記のように吸気弁が閉じるまでCFD演算を行って、その演算結果のデータを化学反応シミュレーションに提供しようとすると、エンジン運転条件の変化に伴い吸気弁の閉時期が変更されることによって、化学反応シミュレーションを開始するタイミングが変化することになる。このことは、CFDから化学反応シミュレーションへと続く一連のプログラムの実行中に、その化学反応シミュレーションにおいて使用する容器モデルの設定を変更しなくてはならないということであり、このようなモデルの変更は非常に煩雑であって、シミュレーションの遅延を招くことになる。
【0008】
本発明は、斯かる点に鑑みてなされたものであり、その目的とするところは、エンジンの性能を予測するためのシミュレーションにおいて、作動ガスの流れを模擬する第1の模擬演算と燃焼等を模擬する第2の模擬演算とを適切に且つ容易に組み合わせて、ダイナミックに解けるようにすることで、シミュレーションの精度を向上しながら、その遅延を招かないようにして、設計・開発の支援ツールとしての実用性を向上することにある。
【0009】
【課題を解決するための手段】
前記目的を達成するために、本発明では、燃焼室における燃焼等を模擬する第2の模擬演算は必ず圧縮行程下死点から膨張行程下死点まで行うようにした上で、第1の模擬演算の結果から圧縮行程下死点における作動ガスの状態を推定して、この推定データを用いて前記第2の模擬演算を行うようにした。
【0010】
具体的に、請求項1の発明は、エンジンの少なくとも吸気系の一部から排気系の一部までに亘る作動ガスの状態をCFDの適用により解析して、該エンジンの性能を予測する予測解析方法を対象とする。このものにおいて、エンジンの作動ガスの流れを模擬する第1の模擬演算を行い、これにより圧縮行程初期に吸気弁が閉じたときに燃焼室に充填されている作動ガスの状態を求め、その吸気弁の閉時期における作動ガスの状態に基づいて、この作動ガスが圧縮行程下死点の燃焼室に充填されていると仮定したときの当該作動ガスの状態を推定する。そして、少なくとも前記推定した圧縮行程下死点における作動ガスの状態に基づいて、そこから膨張行程下死点までの作動ガスの状態を模擬する第2の模擬演算を行うようにする。
【0011】
前記の方法により、エンジンの運転をシミュレーションするときには、まず、そのエンジンの吸排気の流れを模擬する第1の模擬演算を行い、この演算の結果として、圧縮行程初期に吸気弁が閉じた後に確に燃焼室に充填されている吸気の状態を正確に求めることができる。そして、これに基づいて推定した圧縮行程下死点における燃焼室の作動ガスの状態に基づいて、圧縮及び膨張行程における例えば燃焼等の状態を模擬する第2の模擬演算を正確に行うことができる。その際、前記第2の模擬演算は必ず圧縮行程下死点から行うようにすることで、燃焼室の容器モデルの設定を変更する必要がなく、従って、前記第1及び第2の模擬演算同士を組み合わせてダイナミックに解くことが容易に行える。このことで、シミュレーションの精度を向上しながら、その遅延を招くことがなく、設計・開発の支援ツールとしての実用性が向上する。
【0012】
また、請求項2の発明では、特に直噴エンジンのシミュレーションを対象とし、前記請求項1の発明と同様に、第1の模擬演算によって圧縮行程初期に吸気弁が閉じたときに燃焼室に充填されている作動ガスの状態を求めて、これに基づいて圧縮行程下死点の燃焼室における作動ガスの状態を推定し、その推定した作動ガスの状態と燃焼室への燃料供給量とに基づいて、燃焼室の作動ガスの状態を圧縮行程下死点から膨張行程下死点まで模擬する第2の模擬演算を行うようにする。そして、前記第2の模擬演算の途中で圧縮行程における前記燃料供給量が変更されたときには、その演算を一旦、終了し、変更後の燃料供給量に基づいて第2の模擬演算を最初からやり直すことを特徴とする。
【0013】
すなわち、燃焼室に燃料を直接、噴射するようにした直噴エンジンの場合は、シミュレーションによって定常運転状態のみを模擬するのであれば、燃焼室に供給される燃料の量を運転条件毎に設定し、この燃料供給量のデータと第1の模擬演算の結果とに基づいて第2の模擬演算を行うことができる。しかし、エンジンの運転状態が比較的急に変化する過渡運転状態では、前記第2の模擬演算を開始した後の圧縮行程で燃料供給量が変化することがあり、このことが第2の模擬演算に反映されないという問題がある。
【0014】
これに対し、前記請求項2の発明では、まず、請求項1の発明と同様の作用効果が得られるとともに、特に直噴エンジンのシミュレーションの場合には、前記第2の模擬演算の途中で燃料供給量が変更されても、その演算を一旦、終了し、変更後の燃料供給量に基づいてもう一度、最初からやり直すようにしているので、前記第2の模擬演算において圧縮行程で燃焼室に直接、供給される燃料供給量の変化を正確に反映することができる。これにより、直噴エンジンの過渡運転状態も適切に模擬することができる。
【0015】
次に、本願の請求項3の発明は、エンジンの少なくとも吸気系の一部から排気系の一部までに亘る作動ガスの状態をCFDの適用により解析して、該エンジンの性能を予測するためのコンピュータシステムを対象とする。そして、エンジンの作動ガスの流れを模擬する演算を行う第1の演算手段と、前記エンジンの燃焼室における作動ガスの状態を圧縮行程下死点から膨張行程下死点まで模擬する演算を行う第2の演算手段と、前記第1の演算手段による模擬演算の結果から、圧縮行程初期に吸気弁が閉じたときに燃焼室に充填されている作動ガスの状態を求めて、この作動ガスが圧縮行程下死点の燃焼室に充填されていると仮定したときの当該作動ガスの状態を推定する推定手段と、少なくとも該推定手段により推定された作動ガスの状態のデータに基づいて、前記第2の演算手段に対して初期条件となるデータを提供する第1のデータ提供手段と、を備えるものとする。
【0016】
前記のシステムによれば、エンジンの運転シミュレーションの際に、当該エンジンの作動ガスの流れが第1の演算手段によって模擬され、この模擬演算の結果として求められる吸気弁閉時期後の作動ガスの状態から、圧縮行程下死点において燃焼室に充填されていると仮定した作動ガスの状態が推定手段により推定される。そして、少なくとも前記推定された作動ガスの状態のデータが第1のデータ提供手段によって第2の演算手段に提供され、その第2の演算手段によって燃焼室における例えば燃焼の状態が圧縮行程下死点から膨張行程下死点までに亘って模擬される。つまり、前記請求項1の発明に係る予測解析方法が実行されて、その作用効果が得られる。
【0017】
ここで、前記システムにおける第1の演算手段は、吸排気の流れを1次元流及び3次元流の少なくとも一方とみなしてCFD演算を行うものとすればよく、また、第2の演算手段は、圧縮及び膨張行程において燃焼室の作動ガスは運動しないとみなし、且つ吸気弁は圧縮行程下死点で閉じていると仮定して、作動ガスの状態を化学反応式の演算により模擬するものとすればよい(請求項4の発明)。
【0018】
すなわち、第1の模擬演算としてCFDを適用することで、気筒内に充填される吸気の圧力や温度等を正確に求めることが可能になり、これにより、圧縮及び膨張行程における燃焼等の状態を化学反応式により正確に模擬することができる。また、前記CFD演算においては特に精度を高めたい部位について3次元の演算を行い、それ以外は1次元の演算を行うことで、演算量を大幅に減少させることができる。
【0019】
また、前記システムにおける推定手段は、少なくとも、吸気弁閉時期における作動ガスの圧力状態と、この吸気弁閉時期及び圧縮行程下死点の燃焼室の容積比とに基づいて、圧縮行程下死点における作動ガスの圧力状態を推定するものとするのが好ましく(請求項5の発明)、或いは、吸気弁閉時期における作動ガスの圧力及び温度状態と、この吸気弁閉時期及び圧縮行程下死点の燃焼室の容積比と、圧縮行程下死点から吸気弁閉時期までの間に作動ガスから燃焼室壁に伝達される熱量の推定値とに基づいて、圧縮行程下死点における作動ガスの圧力及び温度状態を推定するものとすれば、さらに好ましい(請求項6の発明)。
【0020】
また、前記のシステムには、前記第2の演算手段により演算された膨張行程後期の作動ガスの状態に基づいて、その膨張行程後期に排気弁が開いた後に燃焼室から排出される作動ガスの流れの初期状態を求め、この流れの初期状態のデータを第1の演算手段に対して提供する第2のデータ提供手段を備えるようにしてもよい(請求項7の発明)。
【0021】
すなわち、燃焼室における燃焼の状態等を模擬する第2の模擬演算の結果として、膨張行程後期における燃焼室の圧力、温度、或いは既燃ガスの成分等が求められ、これに基づいて、排気弁が開いたときの燃焼室の既燃ガスの状態、即ち排気流の初期状態が正確に求められるので、その初期状態のデータを第2のデータ手段により第1の演算手段に提供することで、該第1の演算手段により排気の流れに関する正確な模擬演算を行うことができる。
【0022】
さらに、前記エンジン性能の予測解析システムにおいて、前記第1の演算手段を、吸気弁の閉じる時期がエンジンの運転条件に応じて変更されることを加味して、燃焼室への作動ガスの流れを模擬する演算を行うものとするのが好ましい(請求項8の発明)。こうすれば、VVTを備えたエンジンのシミュレーションも正確に行える。
【0023】
さらにまた、特に直噴エンジンのシミュレーションを対象とする場合には、前記請求項3の発明に係るコンピュータシステム(エンジン性能の予測解析システム)と同じ第1及び第2の演算手段と、推定手段とを備えるとともに、少なくとも前記推定手段により推定された作動ガスの状態のデータと燃焼室への燃料供給量のデータとに基づいて前記第2の演算手段に対し演算の初期条件となるデータを提供するデータ提供手段と、前記第2の演算手段による模擬演算の途中で圧縮行程における前記燃料供給量に変化があったときには、その演算を強制的に終了させるとともに、前記推定手段により新たに作動ガスの状態を推定させて、この推定結果に基づいて前記第2の演算手段により最初から再度、模擬演算を実行させる再演算実行手段と、を備えることが好ましい(請求項9の発明)。
【0024】
斯かる構成のコンピュータシステムによれば、前記請求項2の発明に係る予測解析方法を実行して、その作用効果を得ることができる。
【0025】
次に、本願の請求項10の発明は、エンジンの少なくとも吸気系の一部から排気系の一部までに亘る作動ガスの状態をCFDの適用により解析して、該エンジンの性能を予測するためのコンピュータシステムの制御プログラムを対象とする。そして、このプログラムは、エンジンの作動ガスの流れを模擬する演算を行う第1の演算ステップと、前記エンジンの燃焼室における作動ガスの状態を圧縮行程下死点から膨張行程下死点まで模擬する演算を行う第2の演算ステップと、前記第1の演算ステップにおける模擬演算の結果から、圧縮行程初期に吸気弁が閉じたときに燃焼室に充填されている作動ガスの状態を求めて、この作動ガスが圧縮行程下死点の燃焼室に充填されていると仮定したときの当該作動ガスの状態を推定する推定ステップと、少なくとも前記推定ステップにおいて推定された作動ガスの状態のデータに基づいて、前記第2の演算ステップにおける演算の初期条件となるデータを提供する第1のデータ提供ステップと、を備えるものとする。
【0026】
前記のプログラムによってコンピュータシステムを制御することにより、このコンピュータシステムが前記請求項3の発明に係るエンジン性能の予測解析システムとなり、これにより、その請求項3の発明と同じ作用効果が得られる。
【0027】
前記制御プログラムにおける第1の演算ステップでは、吸排気の流れを1次元流及び3次元流の少なくとも一方とみなしてCFD演算を行うようにすればよく、また、第2の演算ステップでは、圧縮及び膨張行程における燃焼室の作動ガスを運動しないとみなし、且つ吸気弁は圧縮行程下死点で閉じていると仮定して、作動ガスの状態を化学反応式の演算により模擬するようにすればよい(請求項11の発明)。こうすれば、請求項4の発明と同じ作用効果が得られる。
【0028】
また、前記制御プログラムにおける推定ステップでは、少なくとも、吸気弁閉時期における作動ガスの圧力状態と、この吸気弁閉時期及び圧縮行程下死点の燃焼室の容積比とに基づいて、圧縮行程下死点における作動ガスの圧力状態を推定するようにすればよく(請求項12の発明)、或いは、吸気弁閉時期における作動ガスの圧力及び温度状態と、この吸気弁閉時期及び圧縮行程下死点の燃焼室の容積比と、圧縮行程下死点から吸気弁閉時期までの間に作動ガスから燃焼室壁に伝達される熱量の推定値とに基づいて、圧縮行程下死点における作動ガスの圧力及び温度状態を推定するようにしてもよい(請求項13の発明)。
【0029】
また、前記制御プログラムには、第2の演算ステップにおいて演算された膨張行程後期の作動ガスの状態に基づいて、その膨張行程後期に排気弁が開いた後に燃焼室から排出される作動ガスの流れの初期状態を求め、この流れの初期状態のデータを第1の演算ステップにおける演算の初期条件として提供する第2のデータ提供ステップをさらに備えることが好ましい(請求項14の発明)。こうすれば、請求項7の発明と同じ作用効果が得られる。
【0030】
また、前記制御プログラムにおける第1の演算ステップでは、吸気弁の閉じる時期がエンジンの運転条件に応じて変更されることを加味して、燃焼室への作動ガスの流れを模擬する演算を行うことが好ましい(請求項15の発明)。こうすれば、請求項8の発明と同じ作用効果が得られる。
【0031】
さらにまた、特に直噴エンジンのシミュレーションを対象とする場合には、前記請求項10の発明に係る制御プログラムと同じ第1及び第2の演算ステップと、推定ステップとを備えるとともに、さらに、少なくとも前記推定ステップにおいて推定された作動ガスの状態と燃焼室への燃料供給量とに基づいて、前記第2の演算ステップにおける演算の初期条件となるデータを提供するデータ提供ステップと、第2の演算ステップにおける模擬演算の途中で圧縮行程における前記燃料供給量に変化のあったときには、その模擬演算を強制的に終了させて、前記推定ステップの作動ガス状態の推定と、この推定結果に基づく前記第2の演算ステップの模擬演算とを最初から実行させる再演算実行ステップとを備えることが好ましい(請求項16の発明)。
【0032】
前記各ステップを備えた制御プログラムによりコンピュータシステムを制御することによって、このコンピュータシステムが前記請求項9の発明に係るエンジン性能の予測解析システムとなり、これにより、その請求項9の発明と同じ作用効果が得られる。
【0033】
【発明の実施の形態】
以下、本発明の実施形態を図面に基いて説明する。
【0034】
(システムの全体構成)
図1は、本発明の実施形態に係るエンジン性能の予測解析システムAの全体構成を示す概念図である。このシステムは、エンジンの作動ガスである吸気や排気等の流れを1次元又は3次元のCFD演算により模擬するとともに(第1の模擬演算)、気筒内の燃焼を化学反応式により記述する演算(第2の模擬演算)を行い、それらを組み合わせることによって、エンジンの運転を模擬するシミュレーションを行うようにしたものである。このシステムの特長は、1次元及び3次元のCFD演算同士でのデータの受け渡しと、CFD演算及び化学反応シミュレーション(化学反応SIM)の間のデータの受け渡しとをいずれも自動化して、例えばスロットル弁から触媒コンバータに至る吸排気の流れをダイナミックに解析することによって、極めて精度の高いシミュレーションを容易に行えるようにしたことにある。
【0035】
図示の符号1,1,…は、主にCFD及び化学反応シミュレーションの演算を実行するコンピュータ装置であり、この実施形態では、特に3次元CFDの膨大な演算量に対応すべく高速のサーバコンピュータを複数台、並列に接続して使用している(以下、演算サーバという)。これら各演算サーバ1は、例えばハードディスクドライブ等の記憶装置を内蔵するとともに、それぞれディスプレイ等の画像表示装置10が接続され、さらに、図示しないが、プリンタ等の出力装置やオペレータによる入力操作を受け付けるキーボード、マウス等の入力デバイスが接続されている。前記記憶装置には、少なくとも、吸排気の流れを模擬するための1次元及び3次元の各CFD演算プログラムと、そのための物理モデルを構築する専用のプリプロセッサと、燃焼状態を模擬する化学反応シミュレーションプログラムと、それら各プログラムによるシミュレーションの結果を画像表示するための画像処理プログラムとが記憶されている。
【0036】
前記演算サーバ1,1,…は、その動作中に必要に応じて一般的な手法により部品データベースDB11にアクセスすることができる。この部品DB11には、1次元及び3次元のCFD演算に用いられるエンジンの物理モデルの雛形がエンジンの各部位毎に種別された状態で予め格納されており、さらに、前記プリプロセッサにより新たに構築されたモデルも格納されるようになっている。前記物理モデルの雛形というのは、例えば吸気系のサージタンク、独立吸気通路、吸気ポート等や排気系の排気ポート、排気マニホルド、EGR通路等のように吸気や排気が流通する部位の基本的な形状を模擬し且つその寸法、形状や材質、表面の状態、熱伝導率等の物理特性値が変更可能な部品モデルであって、以下、この実施形態ではテンプレート部品と呼ぶものとする。
【0037】
そのように寸法、形状や物理特性値を変更可能なテンプレート部品のデータベースDB11を備えることで、この部品DB11から読み込んだテンプレート部品に寸法等を入力して組み合わせるだけで、極めて容易にCFD演算のためのエンジンの物理モデルを構築することができる。また、そのようにして一旦、構築したモデルを新たに部品DB11に格納しておけば、必要に応じてそのモデルの修正も容易に行うことができ、エンジンの設計変更にも容易に対応することができる。
【0038】
また、前記演算サーバ1,1,…は、その動作中に必要に応じて一般的な手法により化学反応データベースDB12にアクセスすることができる。この化学反応DB12は、エンジンの気筒内燃焼室に充填されて燃焼に寄与する吸気中の種々のガス成分(化学種)のうちから代表的なものを、気筒内の状態を表す種々の物理量の組に対応付けて予めグループ化した状態で格納したものである。従って、詳しくは後述するが、CFD演算の結果として得られる気筒内の状態に応じて、これに対応するガス成分のグループを前記化学反応DB12から読み込み、それらガス成分の化学反応をそれぞれ記述することによって、燃焼状態を模擬することができる。
【0039】
図示の符号2は、主にエンジンの諸元値、物理特性及び性能特性を互いに関連づけた実験データのデータベースDB13(実験DB)に接続されて、そのデータの管理を行うコンピュータ装置である(以下、実験DBサーバという)。すなわち、エンジンやトランスミッションに関する過去の実験・開発の過程で蓄積されたデータは、周知の統計的解析手法により整理されて、エンジンの例えば圧縮比や吸気管長さ等の諸元値、その物理特性(例えば体積効率、燃焼特性、損失係数等)及びその性能特性(例えば出力、燃費、エミッション等)を互いに関連づけた実験式として、実験DB13に格納される。そして、この実験式に基づいて、例えばエンジンの諸元値や物理特性からその性能特性を予測することができるようになっている。
【0040】
また、図示の符号3は、エンジンの設計を支援するための3次元CADのコンピュータ装置である(以下、設計CADサーバという)。この設計CADサーバ3は、機械設計や構造解析のための汎用のCADプログラムを実行するとともに、その動作中に必要に応じて一般的な手法により設計データベースDB14(設計DB)にアクセスして、エンジンの設計CADデータを呼び出したり、それらを変更して新たに格納したりすることができるようになっている。すなわち、設計DB14には種々のエンジンの3次元の設計CADデータが、その吸気系、気筒、排気系等の各部位毎に個別に取り出して利用できる状態で格納されている。
【0041】
図示の符号5,5,…は、それぞれパーソナルコンピュータからなる端末(PC端末)であり、これらはパワートレインの設計部門、開発部門、実験部門等に複数台ずつ配置されていて、光通信ケーブル等を用いたネットワーク6によって前記演算サーバ1,1,…、実験DBサーバ2、設計CADサーバ3に双方向通信可能に接続されている。そして、各PC端末5におけるオペレータの操作に従ってシステムの制御プログラムが実行されると、該各PC端末5は前記ネットワーク6を介して演算サーバ1,1,…等に接続されて(ログイン)いわゆるサーバ・クライアント環境を構成し、主に演算サーバ1,1,…との間でコマンドやファイルの授受を行いながら、エンジンの運転シミュレーションを実行するようになっている。
【0042】
尚、前記実験DBサーバ2、設計CADサーバ3及びPC端末5にもそれぞれ演算サーバ1と同様にハードディスクドライブ等の記憶装置が内蔵され、また、ディスプレイ10や出力装置、入力デバイス等が接続されている。
【0043】
(CFD演算)
次に、前記1次元及び3次元のCFD演算について4サイクル4気筒ガソリンエンジンの運転シミュレーションを具体例として説明する。この実施形態では、CFD演算に要する時間をできるだけ短縮するために、基本的には1次元のCFDを基本として必要な部分と行程のみを3次元のCFDで置き換えるようにしている。すなわち、例えば図2(a)〜(d)に示すように、エンジンの吸気通路上流のスロットル弁(図示せず)から第1〜第4気筒c1〜c4の燃焼室を経て触媒コンバータ(図示せず)に至る1次元CFD用の物理モデルMbを基本として、該各気筒c1,c2,…毎に、それぞれが吸気行程にあるときに当該気筒c1,c2,…に対応するサージタンクの一部分s1〜s4のみを3次元のモデルで置換するようにする。
【0044】
より具体的に、図示の1次元のモデルMbでは、基本的には、サージタンクから各気筒までの独立の吸気通路と、スロットル弁からサージタンクまでの各気筒に共通の吸気通路とをそれぞれ管路(図に矢印で示す)の集合体として表し、同様に、各気筒から排気マニホルドの集合部までの独立の排気通路と、その排気集合部から触媒コンバータ入口までの共通の排気通路とをそれぞれ管路の集合体として表す。さらに、前記排気集合部からサージタンクの上流に排気の一部を還流させるEGR通路やサージタンク自体もそれぞれ管路の集合体として表す。また、第1〜第4気筒c1〜c4はそれぞれ容量可変の容器として表す。
【0045】
このような1次元のモデルMbにおいて、管路を流れる吸気や排気の流れはいずれも圧縮性流体の1次元流であるとみなして、その流れの状態を表す圧力p、密度ρ、速度u及び温度Tの各変数について周知の質量保存、運動量保存及びエネルギ保存の式を数値計算により解くことによって、時間的及び空間的に変化する流れの状態を記述することができる。尚、各保存式においては管路の曲がり具合や壁面における摩擦、熱損失等の影響も考慮する。一方、容器についてはその内部の状態は一様で、管路から流入した流体は瞬時に均一に分布すると仮定し、さらに、管路同士や管路と容器との接合部分では適当な境界条件の下で前記保存式を解くようにする。
【0046】
また、例えば独立吸気通路と気筒との接合部分については、吸気弁や排気弁の開閉による断面積の変化を絞りによって模擬するようにするが、その際、シミュレーションするエンジンがVVTを装備したものであれば、これにより例えば吸気弁の開閉時期がエンジンの運転状態に応じて変更されることを模擬するよう、吸気弁や排気弁の開閉時期をエンジンの運転条件に応じて変更し、このことによる影響を加味して、燃焼室への作動ガスの流れを模擬するようにする。こうすることで、VVT付のエンジンのシミュレーションを正確に行うことができる。
【0047】
そして、例えば第1気筒c1が吸気行程にあるときには、同図(a)に示すように、当該第1気筒c1に対応するサージタンクから独立吸気通路の入口までの部位を3次元のモデルs1に置き換えて、その部分の吸気の流れは3次元流として模擬するようにする。すなわち、サージタンクの一部s1とそこから第1気筒c1に向かう独立吸気通路の入口部分との内壁の形状を3次元のモデルで表し、その壁面に沿うようにして流れる吸気の流れについては3次元流として前記各保存式を解くのである。
【0048】
ここで、EGR(Exhaust Gas Recirculation)が行われる場合には、吸気にはエンジンに外部から供給される新しい空気(新気)と排気系から還流される排気(EGRガス)とが含まれており、特にEGRガスには水蒸気や炭酸ガスの他に未燃状態でそれぞれ分子量の異なる種々の炭化水素分子も含まれているから、厳密にはそのようなガスの種類毎に個別に流れの演算を行うことが好ましいとも考えられる。しかし、吸気が輸送される途中でガス成分が変化するわけではないので、この実施形態では、吸気は、新気とEGRガスの2つに分けてそれぞれ流れの変数p,ρ,u,Tを計算し、その計算結果のデータを他のプログラムに受け渡すときには、2つのガス成分の計算結果を合算することで、吸気全体としての流れの変数p,ρ,u,Tを求めるようにしている。
【0049】
例えば、吸気の流れが1次元流から3次元流に変わる場合、1次元流においては前記の如く合算して求めた吸気の変数p,ρ,u,Tはその吸気流の横断面において一様であるから、これをそのまま3次元CFD演算の初期条件或いは境界条件として与えればよい。一方、吸気の流れが3次元流から1次元流に変わるときには、前記の如く新気及びEGRガスについて合算して求めた3次元の吸気流の変数p,ρ,u,Tをその流れの横断面全体について平均化してから、1次元CFD演算の初期条件或いは境界条件として与えればよい。換言すれば、そのように変数を変換しても十分に正確なシミュレーションが行えるように、流れの変数p,ρ,u,Tがその横断面全体についてある程度一様な状態でCFD演算の次元を切換えるようにすればよい。
【0050】
そのように、エンジンの吸気系から排気系に亘る部位のうちの特定の部位のみについて3次元のモデルを用いるとともに、吸気、圧縮、膨張及び排気の各行程のうちから予め選択した行程のみについて3次元のCFD演算を行い、それ以外は0次元又は1次元とみなすようにしたことで、この実施形態では、シミュレーションの精度を十分に確保しながら、そのための演算量は大幅に減少させて、解析に要する時間を短縮することができる。具体的には、図3に一例を示すようにサージタンク全体を3次元で表したモデルを用いて、各気筒の吸気、圧縮、膨張及び排気の全行程に亘って1次元及び3次元のCFD演算を行う従来までのシステムと比較すると、この実施形態のシステムでは演算量は略4分の1にすることができる。
【0051】
また、多数のガス成分を含む吸気の流れをその全てのガス成分毎にCFD演算するのではなく、これを新気及びEGRガスの2つの成分毎に演算するようにしているので、換言すれば、作動ガスを後述する化学反応シミュレーションと比較してかなり少ないガス成分により模擬するようにしているので、この点でも演算量を大幅に減少させて、解析の時間を短縮することができる。尚、CFD演算により求めるのは作動ガスである吸気の流れの状態であり、これは新気とEGRガスとに分けてそれぞれ演算するようにすれば、十分に高い精度が得られる。また、仮にエンジンの運転条件の変化に伴い新気及びEGRガスの割合が変化しても、これに対して作動ガスそのものを変更する必要がない。
【0052】
(化学反応シミュレーション)
上述の如く、各気筒のの吸気及び排気行程における吸排気の流れをそれぞれCFD演算によって模擬するとともに、この実施形態では、圧縮及び膨張行程にある気筒についてはその内部の混合気や燃焼ガス等の運動を無視して、燃焼状態を模擬する化学反応シミュレーションを行うようにしている。具体的には、まず、上述の如き1次元又は3次元のCFD演算によって、気筒の吸気行程で当該気筒内の燃焼室に充填される吸気(新気及びEGRガス)の状態、即ち新気及びEGRガスのそれぞれの圧力p,密度ρ,速度u及び温度Tを合算したものを求める。その際、図4に模式的に示すように、通常、吸気弁の閉時期IVCは圧縮行程下死点BDCよりも遅角側に設定されているから、燃焼室からの吸気の吹き返し等を考慮して、CFD演算は吸気弁が閉じるまでの間、継続して行うようにする。
【0053】
そのようにして、気筒の圧縮行程初期に吸気弁が閉じた時点IVCで燃焼室に充填されている吸気の状態p,ρ,u,TとそのうちのEGRガスの割合とが求められれば、これに基づいて、図に白抜きの矢印にて示すように、作動ガスが圧縮行程下死点BDCの燃焼室に充填されていると仮定したときの当該作動ガスの状態を推定する(データ推定)。すなわち、圧縮行程下死点BDCから吸気弁閉時期IVCまでの間の燃焼室容積の減少と燃焼室壁に伝達される熱量とを考慮して、吸気弁閉時期IVCにおける作動ガスの圧力p及び温度Tと、吸気弁閉時期IVC及び圧縮行程下死点BDCの燃焼室の容積比と、圧縮行程下死点BDCから吸気弁閉時期IVCまでの間に作動ガスから燃焼室壁に伝達される熱量の推定値とに基づいて、圧縮行程下死点BDCにおける作動ガスの圧力p及び温度T等を推定する。尚、圧縮行程下死点から吸気弁閉時期までの間では吸気(空気及びEGRガス)の密度ρは略一定とみなし、また、吸気流速uからは気筒内流動の強さを求めるようにしてもよい。
【0054】
また、前記燃焼室壁への熱の伝達量や燃焼室における混合気の空燃比(又は該燃焼室への燃料の供給量)、さらには燃焼室に残留する既燃ガス(内部EGRガス)の量、気筒壁温等々、燃焼に関与すると考えられる種々の物理量について、シミュレーションにおけるエンジンの運転条件(例えばエンジン負荷と回転速度等)に基づいて求めるようにする。すなわち、この実施形態では、前記熱の伝達量、空燃比、内部EGRガス量、気筒壁温等の物理量の値をエンジンの運転条件に対応付けて予め設定したマップを備え、シミュレーション中のエンジンの運転条件に基づいて前記マップから複数の物理量の値を読み込むようにしている。
【0055】
そして、前記したようにCFD演算の結果とエンジン運転条件とに基づいて、圧縮行程下死点BDCにおける燃焼室の状態を表す複数の物理量の値が求められれば、図5にも示すように、その物理量の組に対応するガス成分のグループを化学反応DB12から読み込むことにより、燃焼状態を模擬するための作動ガスの成分を特定する。そして、そのようにして特定したガス成分を用いて、気筒の圧縮行程下死点から膨張行程下死点までに亘って、それらガス成分同士の化学反応を逐次、記述することで、燃焼状態等を模擬する化学反応シミュレーションを行う。つまり、CFDによるシミュレーションの結果を適切に反映させた正確な化学反応シミュレーションを行うことができる。
【0056】
ここで、前記化学反応DB12におけるガス成分グループのデータは、前記図5に一例を示すように、主に燃料として供給される種々の炭化水素と、空気中の窒素や酸素と、主にEGRガスに含まれる炭化水素、炭酸ガス、窒素酸化物、水蒸気等とのうちから、前記気筒の状態を表す物理量の組に対応する代表的なものを抽出して、その反応式とともに記憶したものである。すなわち、一般に、エンジンの燃焼に関連する化学種及びその素反応を全て挙げれば、これは約3000種類以上にも上るものであり(図6参照)、仮にその全てを演算しようとすれば演算量が著しく多くなってしまい、シミュレーションの時間を徒に長引かせることになる。また、それら全てのデータを燃焼室の状態等に応じてグループ化して格納しようとすれば、化学反応DB12は大きくなり過ぎて、例えば検索時間が長くなる等、種々の不都合を生じる。
【0057】
この点について、全ての化学素反応を挙げるのではなく、燃焼の状態を模擬する上で特に重要なもの、即ち燃焼を模擬する代表的なもののみに絞り込めば、それはせいぜい数十から数百程度で済むので、この実施形態では、エンジンの運転条件によって変化する代表的な化学素反応のみを所定数(例えば100)以下となるように抽出して、これに対応する代表的なガス成分のみを化学反応DB12に格納するようにしている。これにより、化学反応シミュレーションに用いるガス成分の数が適切なものになり、所要の精度をを確保しながら演算量を大幅に減らすことができる。また、化学反応DB12の大きさも適度のものとすることができる。
【0058】
そうして、前記の如く抽出したグループのガス成分(化学種)に基づいて、まず、気筒の圧縮行程では、ピストンの上昇に伴い燃焼室の容積が減少して、その圧力pが上昇し、これに伴い温度Tが上昇することと、気筒壁面との熱交換によって熱を奪われることとを考慮して、そのような条件下における各ガス成分の反応を逐次、記述する。この圧縮行程での化学反応シミュレーションにより、当該気筒において火花点火が行われる前の前炎反応やプレイグニッションの発生等を再現することができる。
【0059】
また、気筒の圧縮上死点近傍では火花点火による混合気の着火を模擬し、これによる化学反応(燃焼)の進行を、気筒の膨張行程における燃焼室容積の増大を加味しつつ膨張行程終了時点まで逐次、記述する。そして、その膨張行程での化学反応シミュレーションの結果として得られる気筒内の既燃ガスの組成、合計の発熱量や気筒壁面との間の熱交換、ピストンに加えられた仕事量、該ピストンの下降に伴う燃焼室容積の拡大等々に基づいて、当該気筒が膨張行程から排気行程に移行するときに燃焼室から排出される既燃ガス(排気)の流れの状態(変数p,ρ,u,T)を求める。
【0060】
すなわち、前記化学反応シミュレーションにより、気筒の膨張行程後期に排気弁が開いたときの燃焼室の圧力p及び温度Tが直接的に求められるとともに、排気の密度ρが気筒内の既燃ガスの組成により求められる。また、圧縮及び膨張行程において気筒内の流動は零とみなしているので、排気流速uの初期値は零となる。そして、図7に模式的に示すように、前記の如く求めた排気弁開時点EVOにおける排気の初期状態p,ρ,u,Tのデータは、同図に破線の矢印で示すように上述したCFD演算プログラムに対して提供され、排気流の境界条件として与えられる。こうして、化学反応シミュレーション(化学反応SIM)により求められる燃焼室からの排気の初期状態をCFD演算の初期条件として用いることで、化学反応シミュレーションの結果を適切に反映させた正確なCFD演算を行うことができる。
【0061】
また、前記のように化学反応シミュレーションにより求めた気筒内の既燃ガスの組成は、化学反応DB12の更新にも用いられる。すなわち、化学反応DB12に格納されているガス成分グループのデータは、上述したように、気筒内の圧力p、温度Tや吸気中のEGRガスの割合等に対応付けて決定されているが、このデータの基になるEGRガスの組成は、予め実験等によって求めたものである。この実施形態では、上述の如く化学反応シミュレーションの結果として求められる排気の組成に基づいて、所定の手法により前記データのガス成分を修正する。例えば、シミュレーションによって得られた排気をそのままEGRガスとみなして、化学反応DB12における対応する運転条件のガス成分グループのデータのいて、炭化水素、一酸化炭素、窒素酸化物等の割合を適当な重み付けでもって修正するようにすればよい。このような修正を行うことにより、ガス成分のデータがCFD演算及び化学反応演算の結果に基づいて修正されることになるので、化学反応シミュレーションの精度がさらに向上する。
【0062】
(シミュレーションの手順)
次に、この実施形態に係るエンジン性能の予測解析システムAによるシミュレーションの手順を具体的に説明する。図8にメインプログラムの概略を示すように、まず、PC端末5,5,…のいずれかにおいて画面表示等に従ってオペレータが所定の入力操作を行うことにより、エンジンシミュレーションのための初期設定データが入力される(S1)。例えば、前記図2に示す4気筒エンジンについて説明すると、このエンジンの諸元値や吸排気系、燃焼室等の寸法・形状を表す幾何データ、それらの熱伝達率等の物理的特性を表す物理データ、或いはそれら詳細なデータに代えて、実験DB13や設計DB14に格納されているエンジンのデータを指定するコード等、さらにはシミュレーションするエンジンの運転条件等をPC端末5に入力させる。
【0063】
また、エンジンのどの部位について3次元のモデルを用いるか選択させ、さらに、その部位について気筒の吸気行程及び排気行程においてはそれぞれ1次元又は3次元の何れのCFD演算を行うか選択させる。すなわち、例えばエンジンの吸気系の設計開発を支援することが解析の目的であれば、オペレータは、図2に示すようにサージタンクに3次元のモデルを用いるとともに、その部位について各気筒毎に吸気行程で3次元のCFD演算を行うことを選択すればよい。こうすれば、サージタンクから独立吸気通路に向かう吸気の流れを3次元流として解析することによって、エンジンの負荷状態や回転速度の変化に拘わらず、体積効率等、エンジンの物理特性値を正確に求めることができ、これによりエンジン出力等の性能特性を正確に予測することができる。
【0064】
続いて、ステップS2では、前記ステップS1において入力された初期設定データに従ってシミュレーションのためのモデルを構築し、これを一旦、保存する。すなわち、例えば図2に示すように、吸気系の一部から排気系の一部までに亘る1次元のCFDモデルMbと、各気筒c1〜c4毎にサージタンクを分割した3次元のCFDモデルs1〜s4とを構築して、それぞれ演算サーバ1,1,…の内部記憶装置に格納する。また、化学反応シミュレーションに関してはクランク角度の変化に対する気筒内容積の変化や気筒壁温に応じた熱伝達率の変化等を規定する容器のモデルを構築する。この容器モデルは、その内部の混合気や燃焼ガスの運動がないものとみなす、という意味において0次元の物理モデルである。
【0065】
より詳しくは、前記3次元CFDモデルを構築するときには、例えば、前記初期設定データに基づいて設計DB14からサージタンクの形状を表す3次元の設計CADデータをPC端末5に読み込み、これに境界面やメッシュの情報を指定するデータを付属したモデル作成コマンドを作成して、演算サーバ1,1,…に送信する。このコマンドを受け取った演算サーバ1,1,…ではプリプロセッサを起動して、サージタンク各部の内壁面にその形状に応じて所定寸法のレイヤーメッシュを貼り付け、また、内部メッシュを切っていくことになる。
【0066】
或いは、前記初期設定データに基づいて別のモデル作成コマンドをPC端末5から演算サーバ1,1,…に送信すると、このコマンドを受け取った演算サーバ1,1,…では、部品DB11からサージタンクの基本的な形状を表すテンプレート部品のデータを読み込んで、この部品の寸法、形状等を変更することにより、CFD演算のためのメッシュを備えた3次元のモデルを構築する。尚、前記3次元のCFDモデルにおいて設計CADデータのメッシュをそのまま使うこともできるようになっている。
【0067】
前記の如く構築したモデルを用いて、ステップS3では、エンジン運転中の吸排気の流れと燃焼室における燃焼の状態とを吸気、圧縮、膨張及び排気の各行程毎に所定の次元で模擬するシミュレーション演算を行う。この演算処理の詳細について一例を挙げれば、この実施形態ではPC端末5と演算サーバ1,1,…との間でプログラムのデータファイルや実行ファイルを相互に送信及び受信しながら、該演算サーバ1,1,…により1次元及び3次元のCFD演算と化学反応シミュレーションとを同時並行的に実行させるようにしている。
【0068】
例えばCFD演算の処理手順としては、まず、1次元CFD演算のモデルMbを読み込み(ステップS31)、シミュレーションの始期における初期条件、即ち吸排気の流れの変数p,ρ,u,Tやエンジンの運転条件等を入力し(S32)、これに基づいて1次元流れの保存式の数値計算を行う(S33)。すなわち、シミュレーションの始期から微小クランク角だけ変化した時点における、スロットル弁下流から各気筒c1〜c4の燃焼室を経て排気通路に至る吸気及び排気の状態(p,ρ,u,T)を、その流れに沿って計算する。
【0069】
その際、図2(a)に示すように第1気筒c1が吸気行程にあれば、当該気筒c1に対応するサージタンクの一部分s1においては、1次元流から3次元流への境界部分の流れの状態(p,ρ,u,T)を求めたところで一旦、1次元の演算を中止し、その演算結果をデータファイルとしてPC端末5に転送する。このファイルを受け取ったPC端末5では1次元流のデータを3次元流のデータに変換して、3次元CFDプログラムの実行ファイルを作成した上で、演算サーバ1,1,…に返送する。その実行ファイルを受け取った演算サーバ1,1,…では3次元CFDプログラムを起動し、まず、第1気筒c1に対応するサージタンクの3次元モデルs1を読み込み(S41)、これに初期条件(前記境界条件)を入力して(S42)、3次元の流れの保存式について数値計算を行い、計算の結果を保存する(S43)。この演算結果のうち、サージタンクと下流側の独立吸気通路との境界部分における流れの変数p,ρ,u,TのデータはPC端末5に転送され、今度は3次元のデータから1次元のデータに変換されて、演算サーバ1,1,…に返送される。
【0070】
そうして、前記の返送されたデータに基づいて1次元CFDプログラムが再開され、独立吸気通路から第1気筒c1の燃焼室、さらにその下流の排気通路へ至る吸排気の流れを計算して、その演算結果を保存する(S33)。このようにして、シミュレーションの始期から微小クランク角変化した後の吸気及び排気の状態(p,ρ,u,T)がエンジンのモデルMb全体に亘って演算されて、この演算結果が保存される。
【0071】
そして、詳しくは後述するが、所定のタイミングで前記CFD演算の結果データの一部を化学反応シミュレーションの結果に基づいて書き換え(データの変換、提供及び書換:S34)、その後、エンジンのクランク角を微小クランク角だけ進めて(インクリメント:S35)、シミュレーションの終期として設定されたクランク角位置になったかどうか判定し(S36)、シミュレーション終期に至るまでは前記ステップS33にリターンして、前記した1次元及び3次元のCFD演算を繰り返し実行する。こうして、エンジンの吸排気の流れがクランク角の変化に対応付けて記憶される。尚、図の例ではスロットル弁の位置に対応する吸気流の境界条件(変数p,ρ,u,T等)は定常運転状態では略一定であり、また、エンジンの運転条件が変化する過渡状態ではその変化に対応するように別途、エンジンの運転制御プログラムから与えられる。
【0072】
上述の如きCFD演算と並行して、圧縮行程及び膨張行程にある気筒についてはそれぞれ化学反応シミュレーション(化学反応SIM)の演算が行われる。すなわち、シミュレーションの進行に伴い例えば第1気筒c1が吸気行程から圧縮行程に移行したときには、前記フローのステップS34において、図9に模式的に示すように、上述したCFD演算による演算結果のデータが演算サーバ1,1,…からPC端末5に送信される。このデータを受け取ったPC端末5では、このデータに基づいて圧縮行程初期の吸気弁閉時点IVCで第1気筒c1に充填されている吸気の状態、即ち圧力p、温度T等や吸気中のEGRガスの割合を求め、これに基づいて、圧縮行程下死点BDCにおける吸気の状態を推定する(図4参照)。また、現在のエンジンの運転条件に基づいて空燃比や気筒壁温等の物理量の値をマップから読み込む。そして、それら気筒内の状態を表す物理量の組を特定して、こに物理量の組に対応する識別コードとともに化学反応シミュレーションプログラムの実行ファイルを演算サーバ1,1,…に送信する(この演算プログラム間でのデータの授受を図9には結果処理*1として示す)。
【0073】
そして、前記実行ファイルを受け取った演算サーバ1,1,…では化学反応シミュレーションプログラムが起動され、前記図8のフローに示すように、第1気筒c1の容器モデルを記憶装置から読み出し(S51)、前記識別コードに対応するガス成分のグループデータを化学反応DB12から読み込んで(化学種の読込み:S52)、予め設定した微小クランク角範囲におけるそれらガス成分の化学反応を記述して、保存する(化学反応演算:S53)。このような化学反応式の演算が当該気筒c1の圧縮行程の下死点から膨張行程の下死点に至るまで、前記微小クランク角毎に繰り返し行われ(S35)、これにより、当該気筒c1内の燃焼室における圧縮及び膨張行程の作動ガスの状態を時系列に記述したデータが化学反応演算の結果として記憶装置に格納される。
【0074】
そうして、前記第1気筒c1が膨張行程を終了して排気行程に移行すれば、当該気筒c1についての化学反応シミュレーションは終了して、図8のフローのステップS34において、今度は前記化学反応演算の結果のデータを演算サーバ1,1,…からPC端末5に送信する。このデータを受け取ったPC端末5では、前記図7にも示したように、第1気筒c1の燃焼室から排出される既燃ガス(排気)の組成や燃焼による発熱、仕事量等に基づいて、圧縮行程後期の排気弁開時点EVOにおける燃焼室の状態を表す変数p,ρ,u(u=0),T、即ち排気の初期状態を求め、これに基づいてCFD演算の演算結果データを書き換えるためのコマンドを作成して、演算サーバ1,1,…に返送する。このコマンドを受け取った演算サーバ1,1,…により、1次元CFD演算の演算結果データにおける燃焼部分、即ち圧縮行程及び膨張行程の部分が書き換えられ、1次元CFD演算における排気流の境界条件として、前記排気流の初期状態のデータが用いられるようになる(この演算プログラム間でのデータの授受を図9には結果処理*2として示す)。また、前記排気の組成に基づいて化学反応DB12におけるガス成分のデータが修正される。
【0075】
上述の如く、メインプログラムのステップS3では、シミュレーションの始期から終期までに亘りエンジンのクランク角の変化に同期して、CFD演算と化学反応シミュレーション演算とが行われる。そして、シミュレーションの終期として設定入力されたクランク角位置になれば(S36でYES)、ステップS4に進んでシミュレーションの結果を出力し、しかる後に制御終了となる(エンド)。前記ステップS4におけるシミュレーション結果の出力としては、演算サーバ1,1,…の記憶装置に保存されている時系列の演算結果のデータのうちから所要のものを読み出してPC端末5に転送し、このデータに基づいて、エンジン性能に関する所定の評価値の出力を行うようにすればよい。すなわち、例えばエンジンの出力特性、燃費特性、エンジン運転条件の変化に伴う各気筒の体積効率の変化等をグラフ化して、サーバ1,1,…やPC端末5のディスプレイに画像表示すればよい。また、特にサージタンク内の吸気の流れ等については3次元CFD演算の結果を可視化して、画像表示するようにしてもよい。
【0076】
前記図8に示すフローのステップS3におけるステップS31〜33,S41〜43が、エンジンの作動ガスである吸排気の流れをCFD演算により模擬する第1の演算ステップに対応している。また、ステップS51〜53が、気筒の圧縮及び膨張行程で該気筒内の燃焼室における作動ガスを複数のガス成分により模擬して、燃焼等の状態を化学反応式の演算により記述する第2の演算ステップに対応している。
【0077】
前記フローのステップS34は、CFD演算の結果から、圧縮行程初期に吸気弁が閉じたときに燃焼室に充填されている吸気の状態を求めて、これが圧縮行程下死点の燃焼室に充填されていると仮定したときの状態を推定する推定ステップに対応するとともに、少なくともその推定された吸気の状態に基づいて化学反応SIMに初期条件となるデータを提供する第1のデータ提供ステップと、反対に、化学反応SIMにより求められた排気流の初期状態のデータをCFD演算プログラムに提供する第2のデータ提供ステップとに対応している。
【0078】
そして、この実施形態の予測解析システムAでは、演算サーバ1,1,…によって前記フローのステップS31〜33,S41〜43が実行されることにより、該演算サーバ1,1,…が第1の演算手段を構成する。また、同様にステップS51〜53が演算サーバ1,1,…により実行されることによって、該演算サーバ1,1,…が第2の演算手段を構成する。さらに、前記フローのステップS34がPC端末5及び演算サーバ1,1,…の協動により実行されることで、それらPC端末5及び演算サーバ1,1,…が推定手段と第1及び第2のデータ提供手段とを構成する。
【0079】
したがって、この実施形態に係るエンジン性能の予測解析システムAによると、例えば4サイクルエンジンの吸排気等の流れをCFDの適用により解析する場合に、基本的には1次元のエンジンモデルMbを用いつつ、予め選択した部位については3次元のモデルs1〜s4を用いて、各気筒毎の吸気及び排気行程についてそれぞれ1次元又は3次元の演算を選択できるようにしたので、CFD演算の精度を十分に高いものとしながら、そのための演算量は大幅に減らすことができる。
【0080】
一方、各気筒毎の圧縮及び膨張行程については少なくとも燃焼室のガス流動を無視して、化学反応シミュレーションにより燃焼状態を模擬するようにしており、その際、膨大な化学反応のうちから代表的なもののみを選択することで、必要なシミュレーション精度を確保しながら、そのための演算量は大幅に減らすことができる。
【0081】
また、EGRを行うようにした場合でも、吸気の流れをガス成分毎にCFD演算するのではなく、新気及びEGRガスの2つの成分として演算するようにしており、このことによって、吸気中のEGRガス割合の変化に柔軟に対応しながら、演算量を大幅に減らすことができる。一方、化学反応シミュレーションでは作動ガスを少なくとも窒素、酸素、炭化水素、炭酸ガス、窒素酸化物及び水蒸気により模擬することで、燃焼に係る化学反応を正確に記述して、精度の高いシミュレーションを行うことができる。
【0082】
そして、前記CFD演算の結果から気筒内の燃焼室に充填される吸気の状態を求めて、これにより化学反応シミュレーションに用いる作動ガスの成分及び状態を特定するとともに、化学反応シミュレーションの結果からはCFD演算における排気の初期状態を求めるようにしており、このように2種類のシミュレーションを適切に組み合わせて、両者をダイナミックに解くことによって、シミュレーションの精度を向上できる。
【0083】
その際、前記燃焼室に充填される吸気の状態を求めるCFD演算は、圧縮行程初期に吸気弁が閉じるまで継続することで、その後に燃焼室に充填されている吸気の状態を正確に求めることができる。そして、その吸気の状態に基づいて圧縮行程下死点における燃焼室の作動ガスの状態を推定し、その作動ガスの成分を化学反応DB12から読み出すことで、正確な化学反応シミュレーションを行うことができる。また、化学反応シミュレーション自体は必ず圧縮行程下死点から行うことで、そのための燃焼室容器モデルの設定を変更する必要がなく、これにより、前記CFD演算と化学反応SIMの組み合わせが容易に行えて、シミュレーションの遅延を招くことがない。
【0084】
以上より、シミュレーションの精度を十分に高くすることができるとともに、そのための演算量は可及的に減少させて、解析に要する時間を短縮することができ、これにより、設計・開発の支援ツールとしての実用性を向上することができる。
【0085】
(他の実施形態)
尚、本発明は前記実施形態に限定されるものではなく、その他の種々の実施形態を包含するものである。すなわち、前記実施形態では、CFD演算を1次元のモデルと3次元のモデルとを組み合わせて行うようにしているが、これに限らず、1次元のCFD演算のみを行うものとしてもよい。
【0086】
また、前記実施形態のCFD演算では吸気を空気及びEGRガスの2つの成分としているが、これに限るものではなく、3つ以上のガス成分からなるものとしてもよいし、反対にEGRガスを含む空気として、1つの成分とみなすことも可能である。
【0087】
また、前記実施形態では、ポート噴射方式のガソリンエンジンのシミュレーションを行う場合について説明をしたが、本願発明の予測解析システムAは直噴エンジンにも適用することができる。但し、直噴エンジンの場合は燃焼室に燃料を直接、噴射するようにしているので、その運転状態が比較的急に変化する過渡運転状態を模擬する場合には、化学反応SIMのための化学反応式の演算を開始した後の気筒の圧縮行程で燃料供給量が変化することがあり、この場合には、その燃料供給量の変化が前記化学反応式の演算に反映されないという問題がある。
【0088】
これに対し、例えば図10に模式的に示すように、直噴エンジンのシミュレーションを行う場合には、化学反応SIMの途中で燃料供給量が変化したときに、その演算を強制的にキャンセルして、同図に仮想線で示すように、変更後の燃料供給量に基づいてガス成分のデータを再度、化学反応DB12から読み込み、このデータを用いて、最初から化学反応式の演算をやり直すようにすればよい(再演算実行ステップ)。こうすれば、圧縮行程の途中で燃焼室に直接、供給される燃料供給量が変化しても、その変化を正確に反映する化学反応SIMを行うことができるから、直噴エンジンの過渡運転状態も適切に模擬することができる。
【0089】
さらに、前記実施形態では、4気筒4サイクルエンジンについてのシミュレーションを行う場合について説明したが、例えば単気筒エンジンのシミュレーションを行うこともできることは言うまでもない。
【0090】
【発明の効果】
以上、説明したように、本願発明に係るエンジン性能の予測解析方法、予測解析システム及びその制御プログラムによると、エンジンの少なくとも吸気系の一部から排気系の一部までに亘る作動ガスの状態を模擬するシミュレーションを行って該エンジンの性能を予測する場合に、吸排気の流れを模擬する第1の模擬演算を圧縮行程初期に吸気弁が閉じる時期まで継続して行い、これにより燃焼室に充填される作動ガスの状態を正確に模擬するとともに、当該燃焼室における燃焼等を模擬するための第2の模擬演算は、必ず、圧縮行程下死点から膨張行程下死点まで行うようにして、燃焼室容器モデルの変更設定等を不要とする。そして、その上で前記第1の模擬演算により求めた吸気弁閉時期における作動ガスの状態に基づいて、圧縮行程下死点における作動ガスの状態を推定し、これにより前記第2の模擬演算を行うようにしたので、両者を適切に組み合わせてダイナミックに解くことが容易になり、シミュレーションの精度を向上しながら、その遅延を招くことがなく、設計・開発の支援ツールとしての実用性を向上できる。
【図面の簡単な説明】
【図1】本発明の実施形態に係るエンジン性能の予測解析システムAの全体構成図。
【図2】CFD演算のためのエンジンモデルの一例を示す図。
【図3】サージタンクを3次元で模擬した従来までのモデルを示す図2相当図。
【図4】吸気行程から圧縮行程への移行時にCFDから化学反応SIMにデータを提供する処理(結果処理*1)を模式的に示す説明図。
【図5】気筒内の状態を表す物理量の組と化学反応DBにおけるガス成分のグループデータとの対応を示す説明図。
【図6】燃焼に関する化学反応の例を示す説明図。
【図7】膨張行程から排気行程への移行時についての結果処理*2を示す図4相当図。
【図8】シミュレーションの手順の概略を示すフローチャート図。
【図9】CFDと化学反応シミュレーションとの切替えと、これに伴うデータの授受とを模式的に示す説明図。
【図10】直噴エンジンに適用する他の実施形態に係る図4相当図。
【図11】CFDから吸気弁閉時期における吸気状態のデータをそのまま化学反応SIMに提供する場合の図4相当図。
【符号の説明】
A エンジン性能の予測解析システム
1,1,… 演算サーバ(第1の演算手段、第2の演算手段、推定手段、第1及び第2のデータ提供手段)
5,5,… PC端末(推定手段、第1及び第2のデータ提供手段)
Claims (16)
- エンジンの少なくとも吸気系の一部から排気系の一部までに亘る作動ガスの状態をCFDの適用により解析して、該エンジンの性能を予測する予測解析方法であって、
エンジンの作動ガスの流れを模擬する第1の模擬演算を行い、これにより圧縮行程初期に吸気弁が閉じたときに燃焼室に充填されている作動ガスの状態を求めて、
前記吸気弁の閉時期における作動ガスの状態に基づいて、その作動ガスが圧縮行程下死点の燃焼室に充填されていると仮定したときの当該作動ガスの状態を推定し、
少なくとも前記推定した圧縮行程下死点における作動ガスの状態に基づいて、そこから膨張行程下死点までの作動ガスの状態を模擬する第2の模擬演算を行うことを特徴とするエンジン性能の予測解析方法。 - 直噴エンジンの少なくとも吸気系の一部から排気系の一部までに亘る作動ガスの状態をCFDの適用により解析して、該エンジンの性能を予測する予測解析方法であって、
エンジンの作動ガスの流れを模擬する第1の模擬演算を行い、これにより圧縮行程初期に吸気弁が閉じたときに燃焼室に充填されている作動ガスの状態を求めて、
前記吸気弁の閉時期における作動ガスの状態に基づいて、その作動ガスが圧縮行程下死点の燃焼室に充填されていると仮定したときの当該作動ガスの状態を推定し、
少なくとも前記推定した圧縮行程下死点における作動ガスの状態と燃焼室への燃料供給量とに基づいて、燃焼室の作動ガスの状態を圧縮行程下死点から膨張行程下死点まで模擬する第2の模擬演算を行うとともに、
第2の模擬演算の途中で圧縮行程における前記燃料供給量が変更されたときにはその演算を一旦、終了し、変更後の燃料供給量に基づいて第2の模擬演算を最初からやり直すことを特徴とするエンジン性能の予測解析方法。 - エンジンの少なくとも吸気系の一部から排気系の一部までに亘る作動ガスの状態をCFDの適用により解析して、該エンジンの性能を予測するためのコンピュータシステムであって、
エンジンの作動ガスの流れを模擬する演算を行う第1の演算手段と、
前記エンジンの燃焼室における作動ガスの状態を圧縮行程下死点から膨張行程下死点まで模擬する演算を行う第2の演算手段と、
前記第1の演算手段による模擬演算の結果から、圧縮行程初期に吸気弁が閉じたときに燃焼室に充填されている作動ガスの状態を求めて、この作動ガスが圧縮行程下死点の燃焼室に充填されていると仮定したときの当該作動ガスの状態を推定する推定手段と、
少なくとも前記推定手段により推定された作動ガスの状態のデータに基づいて、前記第2の演算手段に対して演算の初期条件となるデータを提供する第1のデータ提供手段と、
を備えることを特徴とするエンジン性能の予測解析システム。 - 請求項3において、
第1の演算手段は、吸排気の流れを1次元流及び3次元流の少なくとも一方とみなしてCFD演算を行うものであり、
第2の演算手段は、圧縮及び膨張行程において燃焼室の作動ガスは運動しないとみなし、且つ吸気弁は圧縮行程下死点で閉じていると仮定して、作動ガスの状態を化学反応式の演算により模擬するものである、ことを特徴とするエンジン性能の予測解析システム。 - 請求項3又は4のいずれかにおいて、
推定手段は、少なくとも、吸気弁閉時期における作動ガスの圧力状態と、この吸気弁閉時期及び圧縮行程下死点の燃焼室の容積比とに基づいて、圧縮行程下死点における作動ガスの圧力状態を推定するように構成されていることを特徴とするエンジン性能の予測解析システム。 - 請求項5において、
推定手段は、吸気弁閉時期における作動ガスの圧力及び温度状態と、この吸気弁閉時期及び圧縮行程下死点の燃焼室の容積比と、圧縮行程下死点から吸気弁閉時期までの間に作動ガスから燃焼室壁に伝達される熱量の推定値とに基づいて、圧縮行程下死点における作動ガスの圧力及び温度状態を推定するように構成されていることを特徴とするエンジン性能の予測解析システム。 - 請求項3において、
第2の演算手段により演算された膨張行程後期の作動ガスの状態に基づいて、その膨張行程後期に排気弁が開いた後に燃焼室から排出される作動ガスの流れの初期状態を求め、この流れの初期状態のデータを第1の演算手段に対して提供する第2のデータ提供手段を備える、ことを特徴とするエンジン性能の予測解析システム。 - 請求項3において、
第1の演算手段は、吸気弁の閉じる時期がエンジンの運転条件に応じて変更されることを加味して、燃焼室への作動ガスの流れを模擬する演算を行うように構成されていることを特徴とするエンジン性能の予測解析システム。 - 直噴エンジンの少なくとも吸気系の一部から排気系の一部までに亘る作動ガスの状態をCFDの適用により解析して、該エンジンの性能を予測するためのコンピュータシステムであって、
エンジンの作動ガスの流れを模擬する演算を行う第1の演算手段と、
前記エンジンの燃焼室における作動ガスの状態を圧縮行程下死点から膨張行程下死点まで模擬する演算を行う第2の演算手段と、
前記第1の演算手段による模擬演算の結果から、圧縮行程初期に吸気弁が閉じたときに燃焼室に充填されている作動ガスの状態を求めて、この作動ガスが圧縮行程下死点の燃焼室に充填されていると仮定したときの当該作動ガスの状態を推定する推定手段と、
少なくとも前記推定手段により推定された作動ガスの状態のデータと燃焼室への燃料供給量のデータとに基づいて、前記第2の演算手段に対して演算の初期条件となるデータを提供するデータ提供手段と、
前記第2の演算手段による模擬演算の途中で圧縮行程における前記燃料供給量に変化があったときには、その演算を強制的に終了させるとともに、前記推定手段により新たに作動ガスの状態を推定させて、この推定結果に基づいて前記第2の演算手段により最初から再度、模擬演算を実行させる再演算実行手段と、を備えることを特徴とするエンジン性能の予測解析システム。 - エンジンの少なくとも吸気系の一部から排気系の一部までに亘る作動ガスの状態をCFDの適用により解析して、該エンジンの性能を予測するためのコンピュータシステムの制御プログラムであって、
エンジンの作動ガスの流れを模擬する演算を行う第1の演算ステップと、
前記エンジンの燃焼室における作動ガスの状態を圧縮行程下死点から膨張行程下死点まで模擬する演算を行う第2の演算ステップと、
前記第1の演算ステップにおける模擬演算の結果から、圧縮行程初期に吸気弁が閉じたときに燃焼室に充填されている作動ガスの状態を求めて、この作動ガスが圧縮行程下死点の燃焼室に充填されていると仮定したときの当該作動ガスの状態を推定する推定ステップと、
少なくとも前記推定ステップにおいて推定された作動ガスの状態のデータに基づいて、前記第2の演算ステップにおける演算の初期条件となるデータを提供する第1のデータ提供ステップと、を備えることを特徴とするエンジン性能の予測解析システムの制御プログラム。 - 請求項10において、
第1の演算ステップでは、吸排気の流れを1次元流及び3次元流の少なくとも一方とみなしてCFD演算を行い、
第2の演算ステップでは、圧縮及び膨張行程における燃焼室の作動ガスを運動しないとみなし、且つ吸気弁は圧縮行程下死点で閉じていると仮定して、作動ガスの状態を化学反応式の演算により模擬することを特徴とするエンジン性能の予測解析システムの制御プログラム。 - 請求項10又は11のいずれかにおいて、
推定ステップでは、少なくとも、吸気弁閉時期における作動ガスの圧力状態と、この吸気弁閉時期及び圧縮行程下死点の燃焼室の容積比とに基づいて、圧縮行程下死点における作動ガスの圧力状態を推定することを特徴とするエンジン性能の予測解析システムの制御プログラム。 - 請求項12において、
推定ステップでは、吸気弁閉時期における作動ガスの圧力及び温度状態と、この吸気弁閉時期及び圧縮行程下死点の燃焼室の容積比と、圧縮行程下死点から吸気弁閉時期までの間に作動ガスから燃焼室壁に伝達される熱量の推定値とに基づいて、圧縮行程下死点における作動ガスの圧力及び温度状態を推定することを特徴とするエンジン性能の予測解析システムの制御プログラム。 - 請求項10において、
第2の演算ステップにおいて演算された膨張行程後期の作動ガスの状態に基づいて、その膨張行程後期に排気弁が開いた後に燃焼室から排出される作動ガスの流れの初期状態を求め、この流れの初期状態のデータを第1の演算ステップにおける演算の初期条件として提供する第2のデータ提供ステップを備える、ことを特徴とするエンジン性能の予測解析システムの制御プログラム。 - 請求項10において、
第1の演算ステップでは、吸気弁の閉じる時期がエンジンの運転条件に応じて変更されることを加味して、燃焼室への作動ガスの流れを模擬する演算を行うことを特徴とするエンジン性能の予測解析システムの制御プログラム。 - 直噴エンジンの少なくとも吸気系の一部から排気系の一部までに亘る作動ガスの状態をCFDの適用により解析して、該エンジンの性能を予測するためのコンピュータシステムの制御プログラムであって、
エンジンの作動ガスの流れを模擬する演算を行う第1の演算ステップと、
前記エンジンの燃焼室における作動ガスの状態を圧縮行程下死点から膨張行程下死点まで模擬する演算を行う第2の演算ステップと、
前記第1の演算ステップにおける模擬演算の結果から、圧縮行程初期に吸気弁が閉じたときに燃焼室に充填されている作動ガスの状態を求めて、この作動ガスが圧縮行程下死点の燃焼室に充填されていると仮定したときの当該作動ガスの状態を推定する推定ステップと、
少なくとも前記推定ステップにおいて推定された作動ガスの状態と燃焼室への燃料供給量とに基づいて、前記第2の演算ステップにおける演算の初期条件となるデータを提供するデータ提供ステップと、
第2の演算ステップにおける模擬演算の途中で圧縮行程における前記燃料供給量に変化のあったときには、その模擬演算を強制的に終了させて、前記推定ステップの作動ガス状態の推定と、この推定結果に基づく前記第2の演算ステップの模擬演算とを最初から実行させる再演算実行ステップと、を備えることを特徴とするエンジン性能の予測解析システムの制御プログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003027834A JP4089453B2 (ja) | 2003-02-05 | 2003-02-05 | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003027834A JP4089453B2 (ja) | 2003-02-05 | 2003-02-05 | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004239129A true JP2004239129A (ja) | 2004-08-26 |
JP4089453B2 JP4089453B2 (ja) | 2008-05-28 |
Family
ID=32955450
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003027834A Expired - Fee Related JP4089453B2 (ja) | 2003-02-05 | 2003-02-05 | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4089453B2 (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014211814A (ja) * | 2013-04-19 | 2014-11-13 | 一般財団法人電力中央研究所 | 容器の設計支援システムおよび設計支援プログラム |
JP2016079947A (ja) * | 2014-10-21 | 2016-05-16 | マツダ株式会社 | エンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法及びエンジンの燃焼シミュレーション方法 |
JP2020165341A (ja) * | 2019-03-28 | 2020-10-08 | 国立研究開発法人 海上・港湾・航空技術研究所 | 燃焼機関の実時間性能予測方法及び実時間性能予測プログラム |
CN113484231A (zh) * | 2021-07-09 | 2021-10-08 | 中国人民解放军63796部队 | 一种低温气动截止阀摩擦力在线测试方法 |
CN113569454A (zh) * | 2021-08-04 | 2021-10-29 | 大连交通大学 | 一种柴油机配气机构的模拟分析方法 |
-
2003
- 2003-02-05 JP JP2003027834A patent/JP4089453B2/ja not_active Expired - Fee Related
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014211814A (ja) * | 2013-04-19 | 2014-11-13 | 一般財団法人電力中央研究所 | 容器の設計支援システムおよび設計支援プログラム |
JP2016079947A (ja) * | 2014-10-21 | 2016-05-16 | マツダ株式会社 | エンジンの燃焼シミュレーションに使用する化学反応速度定数の決定方法及びエンジンの燃焼シミュレーション方法 |
JP2020165341A (ja) * | 2019-03-28 | 2020-10-08 | 国立研究開発法人 海上・港湾・航空技術研究所 | 燃焼機関の実時間性能予測方法及び実時間性能予測プログラム |
JP7308480B2 (ja) | 2019-03-28 | 2023-07-14 | 国立研究開発法人 海上・港湾・航空技術研究所 | 燃焼機関の実時間性能予測方法及び実時間性能予測プログラム |
CN113484231A (zh) * | 2021-07-09 | 2021-10-08 | 中国人民解放军63796部队 | 一种低温气动截止阀摩擦力在线测试方法 |
CN113569454A (zh) * | 2021-08-04 | 2021-10-29 | 大连交通大学 | 一种柴油机配气机构的模拟分析方法 |
CN113569454B (zh) * | 2021-08-04 | 2024-04-12 | 大连交通大学 | 一种柴油机配气机构的模拟分析方法 |
Also Published As
Publication number | Publication date |
---|---|
JP4089453B2 (ja) | 2008-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7017399B2 (en) | Predictive analysis method and system for engine performance and control program for use in the same | |
Galindo et al. | Correlations for Wiebe function parameters for combustion simulation in two-stroke small engines | |
Sui et al. | Mean value modelling of diesel engine combustion based on parameterized finite stage cylinder process | |
López et al. | Validation of a 0D/1D computational code for the design of several kind of internal combustion engines | |
JP4089453B2 (ja) | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム | |
JP4175130B2 (ja) | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム | |
JP4175129B2 (ja) | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム | |
JP2015227662A (ja) | 内燃エンジンのシミュレーションの方法とその設備 | |
Valera et al. | Modelling aspects for adaptation of alternative fuels in IC engines | |
JP4192805B2 (ja) | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム | |
JP4192804B2 (ja) | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム | |
Marinoni et al. | RDE cycle simulation by 0D/1D models to investigate IC engine performance and cylinder-out emissions | |
JP4192803B2 (ja) | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム | |
Millo et al. | Optimizing the calibration of a turbocharged GDI engine through numerical simulation and direct optimization | |
Gutiérrez et al. | Predictive thermodynamic model of the performance of a stationary spark-ignition engine running on natural gas | |
Liu et al. | A model for on-line monitoring of in-cylinder residual gas fraction (RGF) and mass flowrate in gasoline engines | |
JP4232654B2 (ja) | エンジン性能の予測解析方法、予測解析システム及びその制御プログラム | |
JP7308480B2 (ja) | 燃焼機関の実時間性能予測方法及び実時間性能予測プログラム | |
Smallbone et al. | Virtual performance and emissions mapping for diesel engine design optimization | |
Finesso et al. | Offline and real-time optimization of EGR rate and injection timing in diesel engines | |
Bozza et al. | Experimental investigation and numerical modelling of an advanced turbocharged DI diesel engine | |
Cordon et al. | One-dimensional engine modeling and validation using Ricardo WAVE | |
Pawar et al. | Virtual engine optimization from design to experimentation | |
Poetsch et al. | A Real-Time Capable and Modular Modeling Concept for Virtual SI Engine Development | |
Sheridan et al. | Defining a model-based calibration process for a twin-independent valve timing engine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20051213 |
|
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: 20080205 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080218 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110307 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120307 Year of fee payment: 4 |
|
LAPS | Cancellation because of no payment of annual fees |