"吸気弁の閉時点にて確定する筒内空気量を同吸気弁閉時点よりも所定の先読み時間だけ先行する時点にて予測する過給機付き内燃機関の空気量推定装置において、吸気管圧力を測定する吸気管圧力測定手段と、前記吸気管圧力を用いて現在の吸気弁流量を計算する吸気弁流量計算手段と、前記吸気弁流量の変化に対するタービン流量の変化の遅れ時間を計算する遅れ時間計算手段と、前記吸気弁流量を前記遅れ時間と前記先読み時間との差分だけ遅らせることで現在よりも前記先読み時間だけ将来のタービン流量を得る将来タービン流量予測手段と、前記将来タービン流量を用いて現在よりも前記先読み時間だけ将来のターボ回転数を計算する将来ターボ回転数予測手段と、前記将来ターボ回転数を用いて現在よりも前記先読み時間だけ将来の吸気管圧を計算する将来吸気管圧予測手段と、前記将来吸気管圧を用いて現在よりも前記先読み時間だけ将来の吸気弁流量を計算する将来吸気弁流量予測手段と、前記将来吸気弁流量を用いて前記吸気弁閉時点における筒内空気量を計算する筒内空気量予測手段と、を備えることを特徴とする過給機付き内燃機関の空気量推定装置。前記将来吸気管圧予測手段は、前記将来ターボ回転数を用いて現在よりも前記先読み時間だけ将来のスロットル上流圧力を計算する手段と、現在よりも前記先読み時間だけ将来のスロットル開度を取得する手段と、前記将来スロットル上流圧力と前記将来スロットル開度とを用いて前記将来吸気管圧を計算する手段と、を備えることを特徴とする請求項1に記載の過給機付き内燃機関の空気量推定装置。前記吸気弁流量を前記遅れ時間だけ遅らせることで現在のタービン流量を得る現在タービン流量推定手段と、前記現在タービン流量を用いて現在のターボ回転数を計算する現在ターボ回転数推定手段と、前記現在ターボ回転数を用いて現在の吸気管圧を計算する現在吸気管圧推定手段と、前記現在吸気管圧を用いて現在の吸気弁流量を計算する現在吸気弁流量推定手段と、前記現在吸気弁流量を用いて現在の筒内空気量を計算する筒内空気量推定手段と、をさらに備えることを特徴とする請求項1又は2に記載の過給機付き内燃機関の空気量推定装置。前記現在吸気管圧推定手段は、前記現在ターボ回転数を用いて現在のスロットル上流圧力を計算する手段と、現在のスロットル開度を取得する手段と、前記現在スロットル上流圧力と前記現在スロットル開度とを用いて前記現在吸気管圧を計算する手段と、を備えることを特徴とする請求項3に記載の過給機付き内燃機関の空気量推定装置。前記吸気管圧力測定手段は、吸気管に取り付けられて吸気管圧力に応じた信号する吸気管圧力センサと、前記吸気管圧力センサの信号を吸気管圧力に変換する手段と、を備えることを特徴とする請求項1乃至4の何れか1項に記載の過給機付き内燃機関の空気量推定装置。前記吸気管圧力測定手段は、コンプレッサとスロットルとの間に取り付けられてスロットル上流圧力に応じた信号するスロットル上流圧力センサと、前記スロットル上流圧力センサの信号をスロットル上流圧力に変換する手段と、前記スロットル上流圧力と現在スロットル開度とを用いて吸気管圧を計算する手段と、を備えることを特徴とする請求項1乃至4の何れか1項に記載の過給機付き内燃機関の空気量推定装置。しかしながら、対象となる内燃機関が過給機付き内燃機関である場合には、スロットルディレイ制御を用いたとしても筒内空気量を精度良く予測することができないおそれがある。スロットルディレイ制御では、先読みした将来のスロットル開度を得ることができるのに対し、従来の過給機のモデルでは、例えば、特開2006-022763号公報に記載されているものがそうであるように、現在のスロットル上流圧力しか推定できないためである。自然吸気型の内燃機関の場合には、時間によるスロットル上流圧力の変化は僅かであるが、過給機付きの内燃機関の場合には、時間によってスロットル上流圧力は大きく変化する。このため、スロットルディレイ制御により将来のスロットル開度を先読みしたとしても、推定したスロットル上流圧力が現在のスロットル上流圧力である場合には、それらを用いて計算した吸気管圧は真の将来の吸気管圧とは異なったものになってしまう。将来の吸気管圧を正確に計算することができなければ、吸気弁の閉時点において確定する筒内空気量を精度良く予測することはできず、結果、目標空燃比の実現に必要な燃料噴射量を正しく算出することは困難となる。本発明は、このような課題に鑑みてなされたもので、過給機付き内燃機関において筒内空気量を精度良く推定できるようにすることを目的とする。第1の発明は、上記の目的を達成するため、吸気弁の閉時点にて確定する筒内空気量を同吸気弁閉時点よりも所定の先読み時間だけ先行する時点にて予測する過給機付き内燃機関の空気量推定装置において、吸気管圧力を測定する吸気管圧力測定手段と、前記吸気管圧力を用いて現在の吸気弁流量を計算する吸気弁流量計算手段と、前記吸気弁流量の変化に対するタービン流量の変化の遅れ時間を計算する遅れ時間計算手段と、前記吸気弁流量を前記遅れ時間と前記先読み時間との差分だけ遅らせることで現在よりも前記先読み時間だけ将来のタービン流量を得る将来タービン流量予測手段と、前記将来タービン流量を用いて現在よりも前記先読み時間だけ将来のターボ回転数を計算する将来ターボ回転数予測手段と、前記将来ターボ回転数を用いて現在よりも前記先読み時間だけ将来の吸気管圧を計算する将来吸気管圧予測手段と、前記将来吸気管圧を用いて現在よりも前記先読み時間だけ将来の吸気弁流量を計算する将来吸気弁流量予測手段と、前記将来吸気弁流量を用いて前記吸気弁閉時点における筒内空気量を計算する筒内空気量予測手段とを備えることを特徴としている。第2の発明は、第1の発明において、前記将来吸気管圧予測手段は、前記将来ターボ回転数を用いて現在よりも前記先読み時間だけ将来のスロットル上流圧力を計算する手段と、現在よりも前記先読み時間だけ将来のスロットル開度を取得する手段と、前記将来スロットル上流圧力と前記将来スロットル開度とを用いて前記将来吸気管圧を計算する手段とを備えることを特徴としている。第3の発明は、第1又は第2の発明において、前記吸気弁流量を前記遅れ時間だけ遅らせることで現在のタービン流量を得る現在タービン流量推定手段と、前記現在タービン流量を用いて現在のターボ回転数を計算する現在ターボ回転数推定手段と、前記現在ターボ回転数を用いて現在の吸気管圧を計算する現在吸気管圧推定手段と、前記現在吸気管圧を用いて現在の吸気弁流量を計算する現在吸気弁流量推定手段と、前記現在吸気弁流量を用いて現在の筒内空気量を計算する筒内空気量推定手段とをさらに備えることを特徴としている。第4の発明は、第3の発明において、前記現在吸気管圧推定手段は、前記現在ターボ回転数を用いて現在のスロットル上流圧力を計算する手段と、現在のスロットル開度を取得する手段と、前記現在スロットル上流圧力と前記現在スロットル開度とを用いて前記現在吸気管圧を計算する手段とを備えることを特徴としている。第5の発明は、第1乃至第4の何れか1つの発明において、前記吸気管圧力測定手段は、吸気管に取り付けられて吸気管圧力に応じた信号する吸気管圧力センサと、前記吸気管圧力センサの信号を吸気管圧力に変換する手段とを備えることを特徴としている。第6の発明は、第1乃至第4の何れか1つの発明において、前記吸気管圧力測定手段は、コンプレッサとスロットルとの間に取り付けられてスロットル上流圧力に応じた信号するスロットル上流圧力センサと、前記スロットル上流圧力センサの信号をスロットル上流圧力に変換する手段と、前記スロットル上流圧力と現在スロットル開度とを用いて吸気管圧を計算する手段とを備えることを特徴としている。第1の発明によれば、現在よりも先読み時間だけ将来の吸気管圧を計算するために、吸気管圧力の測定値に基づいて計算された現在よりも先読み時間だけ将来のターボ回転数が用いられるので、吸気弁の閉時点にて確定する筒内空気量を同吸気弁閉時点よりも先読み時間だけ先行する時点にて精度良く予測することができる。第2の発明によれば、現在よりも先読み時間だけ将来の吸気管圧を計算するために、将来ターボ回転数を用いて計算された将来スロットル上流圧力に加えて先読み時間だけ将来のスロットル開度が用いられるので、将来のスロットル開度に変化がある場合でも筒内空気量を精度良く予測することができる。第3の発明によれば、吸気弁閉時点における筒内空気量を前もって予測することができるだけでなく、吸気管圧力の測定値に基づいて計算された現在のターボ回転数を用いて現在の筒内空気量を精度良く推定することもできる。第4の発明によれば、現在の吸気管圧を計算するために、現在ターボ回転数を用いて計算された現在スロットル上流圧力に加えて現在のスロットル開度が用いられるので、スロットル開度が変化している場合でも現在の筒内空気量を精度良く推定することができる。第5の発明によれば、吸気管圧力センサの信号を処理してられた吸気管圧力の測定値を用いることによって、将来のターボ回転数や現在のターボ回転数を精度良く計算することができる。第6の発明によれば、スロットル上流圧力センサの信号を処理して得られた吸気管圧力の測定値を用いることによって、将来のターボ回転数や現在のターボ回転数を精度良く計算することができる。実施の形態1.本発明の実施の形態1について図を用いて説明する。本実施の形態の空気量推定装置は、過給機付き内燃機関に適用される。本空気量推定装置の適用対象となる内燃機関は、スロットルによる空気量の調整によってトルクを制御することのできる4サイクルレシプロエンジンである。本実施の形態において内燃機関に設けられる過給機は、吸気通路に配置されたコンプレッサを排気通路に配置されたタービンの回転によって駆動するターボ式の過給機である。コンプレッサとスロットルとの間には、コンプレッサによる圧縮によって温度が上昇した空気を冷却するためのインタークーラが設けられている。また、本実施の形態の過給機付き内燃機関は、各気筒の吸気ポートに燃料を噴射するポート噴射型の内燃機関でもある。本空気量推定装置は、内燃機関に備えられるECUの一機能として実現される。詳しくは、メモリに記憶されたプログラムがCPUによって実行されることで、ECUは空気量推定装置として機能する。ECUが空気量推定装置として機能する場合、ECUはプログラムされている空気量推定モデルを用いて筒内空気量を推定する。図1は、本実施の形態で用いられている空気量推定モデルを示す機能ブロック図である。本実施の形態で用いられる空気量推定モデルは、3つの計算モデルM100、M200、M300から構成されている。第1の計算モデルM100は、現在のスロットル開度TAに基づいて筒内空気量を推定するための計算モデルである。第1の計算モデルM100は、エアクリーナモデルM1、コンプレッサモデルM2、インタークーラモデルM3、スロットルモデルM4、吸気管モデルM5、及び吸気弁モデルM6を含んでいる。これらのサブモデルM1、M2、M3、M4、M5、M6の内容については公知であるので、ここではその概略についてのみ示すものとする。エアクリーナモデルM1は、吸気通路の入口に配置されるエアクリーナでの圧力損失を算出するためのモデルであって、数式で表される物理モデル或いは実験データに基づくマップによって構成されている。エアクリーナモデルM1は、大気圧から圧力損失分だけ減算した値をエアクリーナ圧力Pacとして算出する。コンプレッサモデルM2は、過給機のコンプレッサのモデルであって、数式で表される物理モデル或いは実験データに基づくマップによって構成されている。コンプレッサモデルM2では、エアクリーナ圧力Pacと、後述するインタークーラモデルM3で算出された現在インタークーラ圧力Piccrtと、後述する第3の計算モデルで算出された現在ターボ回転数Ntbcrtとが入力され、それらの入力情報から現在のコンプレッサ流量Mcpcrtが算出される。インタークーラモデルM3は、吸気通路におけるインタークーラ部内の空気に関する保存則に基づいて構築された物理モデルである。インタークーラモデルM3としては、具体的にはエネルギー保存則の式と流量保存則の式とが用いられている。インタークーラモデルM3によれば、コンプレッサモデルM2で算出された現在コンプレッサ流量Mcpcrt等の入力情報から現在のインタークーラ圧力Piccrtが算出される。スロットルモデルM4は、スロットルの前後の圧力状態とスロットル開度とからスロットルを通過する空気量の流量を算出するたけのモデルであって、具体的には絞りの式が用いられる。スロットルモデルM4によれば、現在のスロットル開度TA、インタークーラモデルM3で算出された現在インタークーラ圧力Piccrt、及び後述する吸気管モデルM5で算出された現在吸気管圧力Pmcrt等の入力情報から現在のスロットル流量Mtcrtが算出される。吸気管モデルM5は、吸気管内の空気に関する保存則に基づいて構築された物理モデルである。吸気管モデルM5としては、具体的にはエネルギー保存則の式と流量保存則の式とが用いられている。吸気管モデルM5によれば、スロットルモデルM4で算出された現在スロットル流量Mtcrt等の入力情報から現在の吸気管圧力Pmcrtが算出される。吸気弁モデルM6は、吸気弁流量と吸気管圧力との関係について調べた実験ベースのモデルである。実験で得られた経験則により、吸気弁モデルM6においては吸入空気量と吸気管圧力との関係が直線で近似されている。吸気弁モデルM6に現在の吸気管圧力Pmcrtを入力することによって、それにより実現される吸気弁流量Mccrtが算出される。本空気量推定装置は、吸気弁流量Mccrtに基づいて、詳しくは、吸気弁が開いてから閉じるまでの吸気弁流量Mccrtの積算値から筒内空気量を計算する。第2の計算モデルM200は、先読みした将来のスロットル開度TAfwdに基づいて吸気弁の閉時点での筒内空気量を予測するための計算モデルである。実際の内燃機関において筒内空気量が確定する時点は吸気弁の閉時点である。一方、本実施の形態の内燃機関はポート噴射型の内燃機関であるので、各サイクルの燃料噴射は少なくとも吸気弁の閉時点よりも前に完了する必要がある。したがって、筒内に形成される混合気の空燃比を目標空燃比に一致させるためには、吸気弁の閉時点よりも前の燃料噴射時にそのサイクルの確定筒内空気量を先読みすることが求められる。この先読みを可能とするためには、現時点より先の時点のスロットル開度、つまり将来のスロットル開度を先読みする必要がある。そこで、本空気量推定装置は、車両のアクセルペダル操作量に基づいて算出される目標スロットル開度に対し、実際のスロットル開度が所定時間だけ遅れるようにスロットルを動作させる、いわゆるスロットルディレイ制御を実施する。このスロットルディレイ制御によれば、目標スロットル開度は所定時間だけ将来のスロットル開度に相当するので、目標スロットル開度から将来のスロットル開度TAfwdを先読みすることが可能となる。ここでは、スロットル開度TAfwdの実際のスロットル開度TAに対する先読み時間をTfwdとする。先読み時間Tfwdは、燃料噴射量の計算時点燃料噴射の開始時点から吸気弁の開時点までの時間に設定される。第2の計算モデルM200は、第1の計算モデルM100と同じく、エアクリーナモデルM1、コンプレッサモデルM2、インタークーラモデルM3、スロットルモデルM4、吸気管モデルM5、及び吸気弁モデルM6を含んでいる。これらのサブモデルM1,M2,M3,M4,M5,M6間における情報の入出力は、第1の計算モデルM100におけるそれと共通している。ただし、第2の計算モデル200では、後述する第3の計算モデルM300で算出された現在より先読み時間Tfwdだけ将来のターボ回転数NtbfwdがコンプレッサモデルM2に入力されるとともに、先読みした将来のスロットル開度TAfwdがスロットルモデルM4に入力される。コンプレッサモデルM2では、エアクリーナ圧力Pac、将来インタークーラ圧力Picfwd、及び将来ターボ回転数Ntbfwdなどの入力情報から先読み時間Tfwdだけ将来のコンプレッサ流量Mcpfwdが算出される。インタークーラモデルM3では、将来コンプレッサ流量Mcpfwd等の入力情報から現在より先読み時間Tfwdだけ将来のインタークーラ圧力Picfwdが算出される。スロットルモデルM4では、先読みした将来のスロットル開度TAfwd、将来インタークーラ圧力Picfwd、及び将来吸気管圧力Pmfwd等の入力情報から先読み時間Tfwdだけ将来のスロットル流量Mtfwdが算出される。吸気管モデルM5では、将来スロットル流量Mtfwd等の入力情報から先読み時間Tfwdだけ将来の吸気管圧力Pmfwdが算出される。そして、吸気弁モデルM6では、将来吸気管圧力Pmfwdによって実現される吸気弁流量Mcfwdが算出される。吸気弁流量Mcfwdは、現在より先読み時間Tfwdだけ将来の吸気弁流量である。本空気量推定装置は、この先読みした吸気弁流量Mcfwdに基づいて吸気弁の閉弁時点の筒内空気量を計算し、吸気弁閉弁時点の予測筒内空気量と目標空燃比とを用いて燃料噴射量を計算する。第3の計算モデルM300は、吸気管圧力センサ図示略の信号に基づいて現在のターボ回転数Ntbcrtを推定するとともに、現在より先読み時間Tfwdだけ将来のターボ回転数Ntbfwdを予測するための計算モデルである。将来ターボ回転数Ntbfwdを予測することは、過給機付き内燃機関において筒内空気量を精度良く予測する上で特に重要な事項である。筒内空気量はターボ回転数によって決まるスロットル上流圧力とスロットル開度とによって定まるから、燃料噴射量の計算時点において吸気弁の閉弁時点における筒内空気量を精度良く予測するためには、現在より先読み時間Tfwdだけ将来のスロットル開度だけでなく同時間だけ将来のスロットル上流圧力も予測する必要がある。この将来スロットル上流圧力を正確に予測するため、第3の計算モデルM300では将来ターボ回転数Ntbfwdの予測を行っている。第3の計算モデルM300は、前述の吸気弁モデルM6と、2種類のディレイ要素Delay1、Delay2、及び2つのターボ回転数モデルM7を含んでいる。吸気管圧力センサの信号からは吸気管圧力Pmactが測定される。測定された吸気管圧力Pmactは吸気弁モデルM6に入力される。吸気弁モデルM6で用いる吸気管圧力の値として、吸気管モデルM5で算出される吸気管圧力Pmcrtではなく、吸気管圧力センサによる実測値である吸気管圧力Pmactを用いることには理由がある。吸気管圧力Pmcrtを用いる場合には、この計算モデルM300における計算誤差が再び吸気管圧力Pmcrtを介して計算モデルM300に戻ってくる構成になるため、最終的な推定誤差が大きくなってしまう。これに対して、実測値である吸気管圧力Pmactを用いる場合には、この計算モデルM300による誤差の影響を小さくすることができるので、筒内空気量の推定精度は向上する。吸気弁モデルM6では、吸気管圧力Pmactによって実現されている吸気弁流量Mcactが算出される。吸気弁流量Mcactは、2種類のディレイ要素Delay1、Delay2のそれぞれに入力される。吸気弁流量Mcactは、第1のディレイ要素Delay1を通すことで現在のタービン流量Mtbcrtに変換される。一方、第2のディレイ要素Delay2を通された吸気弁流量Mcactは、現在より先読み時間Tfwdだけ将来のタービン流量Mtbfwdに変換される。図2は、吸気弁流量とタービン流量との関係を示す図である。この図に示すように、タービン流量は吸気弁流量に遅れて変化する。その遅れ時間が、第1のディレイ要素Delay1において設定されているディレイ時間Tdly1である。図3は、吸気弁流量の変化に対するタービン流量の変化の遅れ時間の内訳を示す図である。この遅れ時間は、筒内での遅れ時間と排気マニホールドでの遅れ時間とに分けられる。前者は、詳しくは、吸気弁が閉じてから排気弁が開くまでの時間であって、エンジン回転数やバルブタイミングに依存する。したがって、ディレイ時間Tdly1のうちの前者に対応する分は、エンジン回転数とバルブタイミングから算出することができる。後者は、詳しくは、筒内から排出された排気ガスがタービンに到達するまでの時間であって、排気ガスの流速に依存する。したがって、ディレイ時間Tdly1のうちの後者に対応する分は、吸気弁流量Mcactから算出することができる。ただし、排気ガスの流速にはウェストゲートバルブの開度も関係することから、その開度に応じてディレイ時間を補正することが望ましい。また、図2には、タービン流量と先読みしたい将来タービン流量との関係についても示されている。前述のように、筒内空気量の精度良い予測のためには、現在より先読み時間Tfwdだけ将来のターボ回転数を先読みする必要がある。そのためには、吸気弁流量から推定される現在のタービン流量に対し、先読み時間Tfwdだけ将来のタービン流量が必要となる。そこで、第2のディレイ要素Delay2では、吸気弁流量の変化に対するタービン流量の変化の遅れ時間と先読み時間Tfwdとの差分がディレイ時間Tdly2として設定されている。ターボ回転数モデルM7は、過給機の回転挙動のモデルであって、その詳細な構成は図4に示される。図4に示すように、ターボ回転数モデルM7では、現ステップにおけるターボ回転数Ntbを維持することのできるタービン流量Mtb0が現ステップにおける回転数Ntbと定常線タービン流量-ターボ回転数特性線より算出される。そして、現ステップにおけるタービン流量Mtbとターボ回転数維持流量Mtb0との差分ΔMtbが算出される。この差分ΔMtbの大きさはターボ回転数を変化させるエネルギーの大きさに対応することから、差分ΔMtbに変化量係数Aを掛けて得られる値がターボ回転数の変化量ΔNtbとして算出される。変化量係数Aは現ステップにおけるタービン流量Mtbに応じて決定される。そして、現ステップにおけるターボ回転数Ntbにターボ回転数変化量ΔNtbを加えて得られる値が次ステップにおけるターボ回転数Ntbとして算出される。このようなターボ回転数モデルM7に現在タービン流量Mtbcrtが入力されることによって現在ターボ回転数Ntbcrtが算出される。現在ターボ回転数Ntbcrtは、前述のとおり、第1の計算モデルM100のコンプレッサモデルM2において現在コンプレッサ流量Mcpcrtの計算に用いられる。また、このターボ回転数モデルM7に将来タービン流量Mtbfwdが入力されることによって将来ターボ回転数Ntbfwdが算出される。将来ターボ回転数Ntbfwdは、前述のとおり、第2の計算モデルM200のコンプレッサモデルM2において将来コンプレッサ流量Mcpfwdの計算に用いられる。実施の形態2.本発明の実施の形態2について図を用いて説明する。本実施の形態の空気量推定装置と実施の形態1のそれとの違いは使用する空気量推定モデルの構成にある。具体的には、実施の形態1で用いられている空気量推定モデルの第3の計算モデルM300を図5に示す第4の計算モデルM400に置き換えたものが、本実施の形態で用いられる空気量推定モデルである。第4の計算モデルM400は、スロットルモデルM4、吸気管モデルM5、吸気弁モデルM6、2種類のディレイ要素Delay1、Delay2、及び2つのターボ回転数モデルM7を含んでいる。第4の計算モデルM400と第3の計算モデルM300との違いは、測定値である吸気管圧力Pmactを計算する方法にある。第3の計算モデルM300では、吸気管圧力センサの信号から吸気管圧力Pmactが測定されていた。しかし、第4の計算モデルM300では、インタークーラ圧力センサ図示略の信号からインタークーラ圧力Picactが測定される。次に、スロットルモデルM4により、測定されたインタークーラ圧力Picactと現在のスロットル開度TAとからスロットル流量Mtactが算出される。そして、吸気管モデルM5により、スロットル流量Mtactから吸気管圧力Pmactが算出される。吸気管圧力Pmactに基づいて現在ターボ回転数Ntbcrtと将来ターボ回転数Ntbfwdとを計算するための仕組みは第3の計算モデルM300のそれと同じである。このような第4の計算モデルM400を用いる方法は、吸気管圧力センサは有していないがインタークーラ圧力センサを備える内燃機関に用いて好適である。また、吸気管圧力センサとインタークーラ圧力センサの両方を備える内燃機関において、通常は吸気管圧力センサの信号を用いて第3の計算モデルM300による計算を行うが、吸気管圧力センサが故障した場合には、インタークーラ圧力センサの信号を用いた第4の計算モデルM400による計算を行うようにしてもよい。その他.本発明は上述の実施の形態に限定されるものではない。本発明は、その趣旨を逸脱しない範囲で、上述の実施の形態ものから種々変形して実施することができる。例えば、実施の形態で示した各サブモデルのうちのいくつかを統合したり、逆に、1つのサブモデルを複数に分離したりすることも可能である。例えば、特開2006-152899号公報に記載されているように、インタークーラモデルM4とスロットルモデルM5と吸気管モデルM6とを合体させたインタークーラ吸気管合体モデルを用いてもよい。インタークーラ吸気管合体モデルは、スロットルを無視してインタークーラ部から吸気管までを一体化したボリュームとして扱い、そのボリューム内の空気に関する保存則に基づいて構築された物理モデルである。"