(57)【要約】【課題】ＣＶＳ法による対象ガスの排出量の算出結果との間に乖離が生じてしまうことを防止できるガス計測方法を提供する。【解決手段】エンジン２３の排ガス体積流量Ｖｅｘｔｌと排ガス成分濃度ＣＯｄ、ＣＯ2ｄ、ＴＨＣｄ、ＮＯxｄとを計測し、エンジン２３の排ガス体積流量Ｖｅｘｔｌと排ガス成分濃度ＣＯｄ、ＣＯ2ｄ、ＴＨＣｄ、ＮＯxｄとを乗算してモーダルマスエミッションｍＸを算出するガス計測方法であって、臨界ノズル３２における排ガスの気体分子量Ｍ´および比熱比κ´を算出する工程と、排ガスの気体分子量Ｍ´および比熱比κ´の算出結果を用いて、臨界ノズル３２における排ガス成分の流量係数Ｃmを算出する工程と、算出したモーダルマスエミッションｍＸに対して、排ガス成分の流量係数Ｃmと空気の流量係数Ｃ1との比を乗算することで、モーダルマスエミッションｍＸを補正する工程と、を行う。【選択図】図１
2015132552.tif 000002
【特許請求の範囲】
【請求項１】
エンジンの排ガス体積流量と排ガス成分濃度とを計測し、前記エンジンの排ガス体積流量と前記排ガス成分濃度とを乗算して対象ガスの排出量を算出するガス計測方法であって、
臨界ノズルにおける排ガスの気体分子量および比熱比を算出する工程と、
前記排ガスの気体分子量および比熱比の算出結果を用いて、前記臨界ノズルにおける排ガス成分の流量係数を算出する工程と、
前記算出した対象ガスの排出量に対して、前記排ガス成分の流量係数と空気の流量係数との比を乗算することで、前記対象ガスの排出量を補正する工程と、
を行う、
ガス計測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エンジンの排ガス体積流量と排ガス成分濃度とを計測し、エンジンの排ガス体積流量と排ガス成分濃度とを乗算して対象ガスの排出量を算出するガス計測方法に関する。
【背景技術】
【０００２】
従来から、排ガスの認証試験等においては、定容量採取装置（ＣＶＳ：Ｃｏｎｓｔａｎｔ  Ｖｏｌｕｍｅ  Ｓａｍｐｌｉｎｇ）を用いて、希釈空気で希釈した排ガスからＣＯ等の対象ガスの排出量（マスエミッション）を計測するＣＶＳ法が標準的に用いられている。
【０００３】
図４に示すように、ＣＶＳ法では、希釈空気と排ガスとを希釈用トンネル３１に導入し、排ガスを希釈空気で希釈する。ＣＶＳ法では、ポンプ３５によって希釈用トンネル３１で希釈した排ガス（以下、「希釈排ガス」と表記する）を吸引する。希釈排ガスは、臨界ノズル３２を通ってポンプ３５に吸引される。
ＣＶＳ法では、希釈空気を希釈空気バッグ３３に採取するとともに希釈排ガスをサンプルバッグ３４に採取し、採取した希釈空気および希釈排ガスの成分濃度を分析する。
【０００４】
ＣＶＳ法において、マスエミッションｍＸ０は、例えば、以下の数式１によって算出される。
【数１】
2015132552.tif 000003
前記数式１において、ρは対象ガス成分密度、Ｃsambagはサンプルバッグ３４内の対象ガス濃度（ｐｐｍ）、Ｃambbagは希釈空気バッグ３３内の対象ガス濃度（ｐｐｍ）、ＤＦは排ガスの希釈率、Ｖmixは希釈排ガス総体積である。
【０００５】
前記数式１における希釈排ガス総体積Ｖmixは、以下の数式２によって算出される。
【数２】
2015132552.tif 000004
前記数式２において、Ｃ1は流量係数、Ｐvは希釈用トンネル３１を流れる希釈排ガスの圧力、Ｔvは希釈用トンネル３１を流れる希釈排ガスの温度、ｔは時間である。
【０００６】
前記数式２において、流量係数Ｃ1は、空気を媒体に校正される。一方、希釈用トンネル３１を流れているのは、希釈排ガスである。
従って、ＣＶＳ法においては、希釈排ガスと空気との物性値の違いが誤差要因となってしまう。
【０００７】
車両の開発段階においては、排ガス体積流量と排ガス成分濃度の積から、排ガスを希釈することなく対象ガスの排出量を直接的に算出するモーダルマス法が採用される場合がある（例えば、特許文献１参照）。
特許文献１に開示される技術では、燃焼直後の排ガス体積流量と燃焼直後の水分濃度との算出結果に基づいてモーダルマスエミッションを算出することで、排ガス管内における水分濃度変化によらずモーダルマスエミッションを正確に算出している。
【０００８】
前述のように、ＣＶＳ法では、物性値の違いによってマスエミッションｍＸ０の計測結果に誤差が生じてしまう。
従って、モーダルマス法によってモーダルマスエミッションを正確に算出しても、ＣＶＳ法によるマスエミッションｍＸ０の算出結果との間には、乖離が生じてしまう可能性がある。
つまり、従来技術においては、モーダルマス法が採用される車両の開発段階でモーダルマスエミッションが所定の値になったとしても、ＣＶＳ法が採用される排ガスの認証試験でマスエミッションｍＸ０が所定の値とならない場合がある。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】特開２０１３−１８１８０８号公報
【発明の概要】
【発明が解決しようとする課題】
【００１０】
本発明は、以上の如き状況を鑑みてなされたものであり、ＣＶＳ法による対象ガスの排出量の算出結果との間に乖離が生じてしまうことを防止できるガス計測方法を提供するものである。
【課題を解決するための手段】
【００１１】
本発明に係るガス計測方法は、エンジンの排ガス体積流量と排ガス成分濃度とを計測し、前記エンジンの排ガス体積流量と前記排ガス成分濃度とを乗算して対象ガスの排出量を算出するガス計測方法であって、臨界ノズルにおける排ガスの気体分子量および比熱比を算出する工程と、前記排ガスの気体分子量および比熱比の算出結果を用いて、前記臨界ノズルにおける排ガス成分の流量係数を算出する工程と、前記算出した対象ガスの排出量に対して、前記排ガス成分の流量係数と空気の流量係数との比を乗算することで、前記対象ガスの排出量を補正する工程と、を行う、ものである。
【発明の効果】
【００１２】
本発明は、ＣＶＳ法による対象ガスの排出量の算出結果との間に乖離が生じてしまうことを防止できる、という効果を奏する。
【図面の簡単な説明】
【００１３】
【図１】ガス計測方法による算出手順を示す説明図。
【図２】水分濃度算出手段の演算部の説明図。
【図３】誤差補正手段による演算内容を示す説明図。
【図４】ＣＶＳ法による算出手順を示す説明図。
【発明を実施するための形態】
【００１４】
以下では、本実施形態のガス計測方法について説明する。
【００１５】
図１に示すように、ガス計測方法は、排ガスを希釈することなく直接的に計測するモーダルマス法によって、モーダルマスエミッションｍＣＯ・ｍＣＯ2・ｍＴＨＣ・ｍＮＯx（ＣＯ等の対象ガスの排出量）［ｇ］を算出するものである。
ガス計測方法は、例えば、車両の開発段階において好適に用いられる。ガス計測方法は、ガス計測装置１０を用いて実施される。
【００１６】
まず、本実施形態のガス計測方法の計測対象となる排ガスの経路２０について説明する。
【００１７】
排ガスの経路２０は、自動車のエンジン２３の動作において、自動車に空気が吸入されてからその空気がエンジン２３を経て排ガスとなり、外部に排出されるまでの経路である。
【００１８】
排ガスの経路２０上には、空気および排ガスの流れ方向より順番に、空気吸入口２１、エアクリーナ２２、エンジン２３、触媒２４、サイレンサ２５、排ガス排出口２６等が設けられる。
触媒２４、サイレンサ２５、および排ガス排出口２６等は、排ガス管によって連結される。サイレンサ２５と排ガス排出口２６との間の符号２７は、排ガス管のテールエンド２７（末端）を示すものである。
【００１９】
経路２０においては、排気ガス管途中の局所的放熱や結露した水分の再蒸発により水分濃度が変化する。
【００２０】
次に、ガス計測装置１０の構成について説明する。
【００２１】
ガス計測装置１０は、流量計１１、排ガス分析計１２、燃焼直後（燃焼反応直後）のガス組成算出手段１３、燃焼直後の排ガス体積流量算出手段１４、排ガス管のテールエンド２７における水分濃度算出手段１５、モーダルマスエミッション算出手段１６、および誤差補正手段１７を具備する。
【００２２】
流量計１１は、排ガス排出口２６とテールエンド２７との間に配置される。流量計１１は、テールエンド２７における排ガス体積流量Ｖｅｘｔｌ［Ｌ／ｓ］、すなわち、水分濃度変化後の排ガス体積流量を計測する。
【００２３】
なお、排ガス体積流量Ｖｅｘｔｌに付される添え字（ｅｘｔｌ）のうち、「ｔｌ」は、テール（ｔａｉｌ）、すなわち、テールエンド２７での計測を意味する。これは、水分濃度等の他の値においても同様である。
また、添え字「ｅｘ」はエキゾースト（ｅｘｈａｕｓｔ：排気）を意味する。
【００２４】
排ガス分析計１２は、排ガス排出口２６の近傍に配置される。排ガス分析計１２は、テールエンド２７における排ガス（自動車より外部に排出される排ガス）を分析し、ドライベースの排ガス成分濃度ＣＯｄ、ＣＯ2ｄ、ＴＨＣｄ、ＮＯxｄを計測する。
【００２５】
なお、排ガス成分濃度ＣＯｄ等に付される添え字「ｄ」は、ドライ計測値（除湿した濃度値）を意味する。
【００２６】
排ガス分析計１２における分析では、微量なサンプルガスを除湿するため、水分濃度変化の影響を受けない。なお、燃焼直後からテールエンド２７まで成分組成は変化しない。
【００２７】
燃焼直後のガス組成算出手段１３は、排ガス分析計１２によって計測されたドライベースの排ガス成分濃度ＣＯｄ等を用いて、燃焼直後（燃焼直後）の高精度なガス組成（濃度）ＣＯｗｃ、ＣＯ2ｗｃ、ＴＨＣｗｃ、ＮＯxｗｃ、Ｎ2ｗｃ、Ｏ2ｗｃ、Ｈ2ｗｃ、Ｈ2Ｏｗｃを算出する。
【００２８】
なお、燃焼直後のガス組成ＣＯｗｃ等に付される添え字「ｗｃ」は、計算（ｃａｌｃｕｌａｔｉｏｎ）でウェット（ｗｅｔ）状態の濃度を算出したことを意味する。
【００２９】
燃焼直後の排ガス体積流量算出手段１４は、吸入空気量ＧＡ［ｇ／ｓ］および空燃比（吸入空気量と燃料噴射量の質量比）Ａ／Ｆにより求まる燃焼直後の排ガス質量流量Ｇｅｘｈ［ｇ／ｓ］等を用いて、燃焼直後の排ガス体積流量Ｖｅｘｈ［Ｌ／ｓ］を算出する。
空燃比制御されたエンジン２３では、ＡＦＭ（エアフローメータ）による吸入空気量ＧＡ［ｇ／ｓ］とＡ／Ｆセンサによる空燃比Ａ／Ｆとは、エンジンコントロールユニット（ＥＣＵ：Ｅｎｇｉｎｅ  Ｃｏｎｔｒｏｌ  Ｕｎｉｔ）に記録される。
【００３０】
なお、排ガス質量流量Ｇｅｘｈに付される添え字（ｅｘｈ）のうち、「ｈ」は、ヘッド（ｈｅａｄ）、すなわち、排気ガス管のエンジン側端部（燃焼直後部）を意味する。
これは、排ガス体積流量等の他の値においても同様である。
【００３１】
テールエンド２７における水分濃度算出手段１５は、ガス組成算出手段１３で算出された燃焼直後の水分濃度Ｈ2Ｏｗｃ等を用いて、テールエンド２７における水分濃度、すなわち、排ガス間途中の局所的放熱等で変化した後の水分濃度Ｈ2Ｏｔｌを算出する。
【００３２】
モーダルマスエミッション算出手段１６は、流量計１１および排ガス分析計１２によって計測されたテールエンド２７における排ガス体積流量Ｖｅｘｔｌ［Ｌ／ｓ］およびドライベースの排ガス成分濃度ＣＯｄ等を用いて、モーダルマスエミッションｍＣＯ・ｍＣＯ2・ｍＴＨＣ・ｍＮＯx［ｇ］を算出する。
【００３３】
誤差補正手段１７は、モーダルマスエミッション算出手段１６で算出されたモーダルマスエミッションｍＣＯ・ｍＣＯ2・ｍＴＨＣ・ｍＮＯx［ｇ］に対して演算処理を施して、モーダルマスエミッションｍＣＯ・ｍＣＯ2・ｍＴＨＣ・ｍＮＯx［ｇ］を補正する。
【００３４】
次に、ガス計測方法の手順について説明する。
【００３５】
まず、ガス計測方法では、エンジン２３を駆動させ、流量計１１によってテールエンド２７における排ガス体積流量Ｖｅｘｔｌ［Ｌ／ｓ］を計測する。
また、ガス計測方法では、排ガス分析計１２によってドライベースの排ガス成分濃度ＣＯｄ、ＣＯ2ｄ、ＴＨＣｄ、ＮＯxｄを計測する。
【００３６】
このように、ガス計測方法では、エンジン２３の排ガス体積流量Ｖｅｘｔｌ［Ｌ／ｓ］と排ガス成分濃度ＣＯｄ、ＣＯ2ｄ、ＴＨＣｄ、ＮＯxｄとを計測する。
【００３７】
ガス計測方法では、ガス組成算出手段１３によって燃焼直後の高精度なガス組成ＣＯｗｃ、ＣＯ2ｗｃ、ＴＨＣｗｃ、ＮＯxｗｃ、Ｎ2ｗｃ、Ｏ2ｗｃ、Ｈ2ｗｃ、Ｈ2Ｏｗｃを算出する。
このとき、ガス組成算出手段１３は、演算部Ｆ２を用いて以下のような演算処理を行う。
【００３８】
演算部Ｆ２は、以下の燃焼反応式より、燃焼反応後の各成分ＣＯ、ＣＯ2、ＴＨＣ、ＮＯx、Ｎ2、Ｏ2、Ｈ2Ｏ、Ｈ2のモル数を算出する。
φＣＨｎ＋（１＋ｎ／４）Ｏ2＋（７９／２１）×（１＋ｎ／４）Ｎ2
→ａＣＯ＋ｂＣＯ2＋ｃＨ2Ｏ＋ｄＣＨｎ＋ｅＨ2＋ｆＯ2＋ｇＮＯ＋ｈＮ2
前記燃焼反応式において、φは当量比（理論空燃比／空燃比）、ｎはＨ／Ｃ（燃料性状：水素原子数／炭素原子数の比）、ａ〜ｈは燃焼反応後の各成分ＣＯ、ＣＯ2、Ｈ2Ｏ、ＴＨＣ、Ｈ2、Ｏ2、ＮＯx、Ｎ2のモル数である。
【００３９】
前記燃焼反応式においては、各モル数ａ〜ｈが未知数である。
このため、演算部Ｆ２は、原子数保存（質量保存の法則）より与えられる以下の１）〜４）式と、水性ガス平衡定数Ｋより与えられる５）式とによって、各モル数ａ〜ｈを算出する。
１）式  φ＝ａ＋ｂ＋ｄ  （炭素Ｃ原子数バランス）
２）式  （１＋ｎ／４）＝ａ／２＋ｂ＋ｃ／２＋ｆ＋ｇ／２  （酸素Ｏ原子数バランス）
３）式  ｎ×φ＝２ｃ＋ｎ×ｄ＋２ｅ  （水素Ｈ原子数バランス）
４）式  （７９／２１）（１＋ｎ／４）＝ｇ／２＋ｈ  （窒素Ｎ原子数バランス）
５）式  （ａ×ｃ）／（ｂ×ｅ）＝Ｋ    （水性ガス平衡定数Ｋは３．０〜４．５の定数）
【００４０】
演算部Ｆ２は、排ガス分析計１２で計測したドライベースの排ガス成分濃度ＣＯｄ、ＣＯ2ｄ、ＴＨＣｄ、ＮＯxｄと各モル数ａ〜ｈとの紐付け（関連付け）を、以下のようにして行う。
【００４１】
ドライベースの排ガス成分濃度ＣＯｄ等の比と燃焼後の各モル数ａ〜ｈの比とは、同じになる。
従って、ドライベースの排ガス成分濃度ＣＯｄ、ＣＯ2ｄ、ＴＨＣｄ、ＮＯxｄを、ＣＯｄ→Ａ、ＣＯ2ｄ→Ｂ、ＴＨＣｄ→Ｄ、ＮＯx→Ｇとした場合には、以下の関係が成り立つ。
ａ：ｂ＝Ａ：Ｂ  ⇒  ｂ＝Ｂ／Ａ×ａ
ｄ：ａ＝Ｄ：Ａ  ⇒  ｄ＝Ｄ／Ａ×ａ
ａ：ｇ＝Ａ：Ｇ  ⇒  ｇ＝Ｇ／Ａ×ａ
【００４２】
演算部Ｆ２は、このような燃焼反応後のＣＯ2成分およびＴＨＣのモル数ｂ・ｄの式を前記１）式に代入する。
１）式  φ＝ａ＋ｂ＋ｄ
＝ａ  ＋  Ｂ／Ａ×ａ  ＋  Ｄ／Ａ×ａ
【００４３】
演算部Ｆ２は、前記代入した当量比φの式を以下のように整理することで、燃焼反応後のＣＯ成分のモル数ａとドライベースの排ガス成分濃度ＣＯｄ等との紐付けを行う。
ａ＝φＡ／（Ａ＋Ｂ＋Ｄ）
演算部Ｆ２は、ドライベースの排ガス成分濃度ＣＯｄ等の計測結果を前記ＣＯ成分のモル数ａの式に代入し、燃焼反応後のＣＯ成分のモル数ａを算出する。
【００４４】
演算部Ｆ２は、燃焼反応後のＣＯ成分のモル数ａの算出と同じ要領で１）式を整理して、燃焼反応後のＣＯ2成分およびＴＨＣ成分のモル数ｂ・ｄとドライベースの排ガス成分濃度ＣＯｄ等との紐付けを行う。
ｂ＝φＢ／（Ａ＋Ｂ＋Ｄ）
ｄ＝φＤ／（Ａ＋Ｂ＋Ｄ）
演算部Ｆ２は、ドライベースの排ガス成分濃度ＣＯｄ等の計測結果を前記ＣＯ2成分およびＴＨＣ成分のモル数ｂ・ｄの式に代入し、ＣＯ2成分およびＴＨＣ成分のモル数ｂ・ｄを算出する。
【００４５】
演算部Ｆ２は、前記式ｇ＝Ｇ／Ａ×ａをａ＝ｇ×Ａ／Ｇに変形して、前記燃焼反応後のＣＯ成分のモル数ａの式（ａ＝φＡ／（Ａ＋Ｂ＋Ｄ））に代入する。演算部Ｆ２は、代入した式を以下のように整理することで、燃焼反応後のＮＯx成分のモル数ｇとドライベースの排ガス成分濃度ＣＯｄ等との紐付けを行う
ｇ＝φＧ／（Ａ＋Ｂ＋Ｄ）
演算部Ｆ２は、ドライベースの排ガス成分濃度ＣＯｄ等の計測結果を前記ＮＯx成分のモル数ｇの式に代入して、燃焼反応後のＮＯx成分のモル数ｇを算出する。
【００４６】
演算部Ｆ２は、前記３）〜５）式を用いて、燃焼反応後のＨ2成分、Ｎ2成分、およびＨ2Ｏ成分の各モル数ｅ・ｈ・ｃとドライベースの排ガス成分濃度ＣＯｄ等との紐付けを行う。
【００４７】
まず、演算部Ｆ２は、前記５）式を以下のように変形する。
ｃ＝Ｋ×ｅ×ｂ／ａ
【００４８】
演算部Ｆ２は、前記Ｈ2Ｏ成分のモル数ｃの式を前記３）式に代入する。
ｎ×φ＝２Ｋ×ｅ×ｂ／ａ＋ｎ×ｄ＋２ｅ
【００４９】
演算部Ｆ２は、前記ｎ×φの式を以下のように整理する。
ｅ＝ｎ×（φ−ｄ）／｛２（Ｋ×ｂ／ａ＋１）｝
【００５０】
演算部Ｆ２は、前記Ｈ2成分のモル数ｅの式に、ＣＯ成分のモル数ａの算出で得られた式ａ＝φＡ／（Ａ＋Ｂ＋Ｄ）、ｂ＝φＢ／（Ａ＋Ｂ＋Ｄ）、ｄ＝φＤ／（Ａ＋Ｂ＋Ｄ）を代入する。
これにより、演算部Ｆ２は、以下のような燃焼反応後のＨ2成分のモル数ｅの式を得る。
ｅ＝ｎ×｛φ−φＤ／（Ａ＋Ｂ＋Ｄ）｝／｛２（Ｋ×Ｂ／Ａ＋１）｝
【００５１】
演算部Ｆ２は、前記Ｈ2成分のモル数ｅの式の右辺の分子および分母に対してＡ×（Ａ＋Ｂ＋Ｄ）をそれぞれ乗算して整理し、燃焼反応後のＨ2成分のモル数ｅとドライベースの排ガス成分濃度ＣＯｄ等との紐付けを行う。
ｅ＝ｎ×Ａ×｛φ（Ａ＋Ｂ＋Ｄ）−φＤ｝／｛２（Ｋ×Ｂ＋Ａ）×（Ａ＋Ｂ＋Ｄ）｝
＝ｎφＡ（Ａ＋Ｂ）／｛２（ＫＢ＋Ａ）（Ａ＋Ｂ＋Ｄ）｝
演算部Ｆ２は、ドライベースの排ガス成分濃度ＣＯｄ等の計測結果を前記Ｈ2成分のモル数ｅの式に代入して、燃焼反応後のＨ2成分のモル数ｅを算出する。
【００５２】
演算部Ｆ２は、前記４）式を以下のように変形する。
ｈ＝（７９／２１）（１＋ｎ／４）−ｇ／２
【００５３】
演算部Ｆ２は、前記Ｎ2成分のモル数ｈの式に、ＣＯ成分のモル数ａの算出で得られた式ｇ＝φＧ／（Ａ＋Ｂ＋Ｄ）を代入する。
これにより、演算部Ｆ２は、燃焼反応後のＮ2成分のモル数ｈとドライベースの排ガス成分濃度ＣＯｄ等との紐付けを行う。
ｈ＝（７９／２１）（１＋ｎ／４）−φＧ／｛２（Ａ＋Ｂ＋Ｄ）｝
演算部Ｆ２は、ドライベースの排ガス成分濃度ＣＯｄ等の計測結果を前記Ｎ2成分のモル数ｈの式に代入して、燃焼反応後のＨ2成分のモル数ｈを算出する。
【００５４】
演算部Ｆ２は、前記５）式を変形した前記Ｈ2Ｏ成分のモル数ｃの式ｃ＝Ｋ×ｂ／ａ×ｅに、ＣＯ成分のモル数ａの算出で得られた式ａ＝φＡ／（Ａ＋Ｂ＋Ｄ）およびｂ＝φＢ／（Ａ＋Ｂ＋Ｄ）と、Ｈ2成分のモル数ｅの算出で得られた式ｅ＝ｎφＡ（Ａ＋Ｂ）／｛２（ＫＢ＋Ａ）（Ａ＋Ｂ＋Ｄ）｝とを代入する。
ｃ＝Ｋ×Ｂ／Ａ×ｎφＡ（Ａ＋Ｂ）／｛２（ＫＢ＋Ａ）（Ａ＋Ｂ＋Ｄ）｝
【００５５】
演算部Ｆ２は、前記Ｈ2Ｏ成分のモル数ｃの式を以下のように整理して、燃焼反応後のＨ2Ｏ成分のモル数ｃとドライベースの排ガス成分濃度ＣＯｄ等との紐付けを行う。
ｃ＝ＫｎφＢ（Ａ＋Ｂ）／｛２（ＫＢ＋Ａ）（Ａ＋Ｂ＋Ｄ）｝
演算部Ｆ２は、ドライベースの排ガス成分濃度ＣＯｄ等の計測結果を前記Ｈ2Ｏ成分のモル数ｃの式に代入して、燃焼反応後のＨ2Ｏのモル数ｃを算出する。
【００５６】
ここで、燃焼反応後の各モル数ａ〜ｈの総和は、１００％（＝１）である。従って、燃焼反応後の各モル数ａ〜ｈについては、以下の式が成り立つ。
ｆ＝１−（ａ＋ｂ＋ｃ＋ｄ＋ｅ＋ｇ＋ｈ）
【００５７】
演算部Ｆ２は、前記Ｏ2成分のモル数ｆの式に、各モル数ａ〜ｅ・ｇ・ｈとドライベースの排ガス成分濃度ＣＯｄ等とを紐付けた結果を代入する。
これにより、演算部Ｆ２は、燃焼反応後のＯ2成分のモル数ｆとドライベースの排ガス成分濃度ＣＯｄ等との紐付けを行う。
ｆ＝１−φ（Ａ＋Ｂ＋Ｄ＋Ｇ）／（Ａ＋Ｂ＋Ｄ）
−（７９／２１）（１＋ｎ／４）＋φＧ／｛２（Ａ＋Ｂ＋Ｄ）｝
−ｎφ（ＫＢ＋Ａ）（Ａ＋Ｂ）／｛２（ＫＢ＋Ａ）（Ａ＋Ｂ＋Ｄ）｝
演算部Ｆ２は、ドライベースの排ガス成分濃度ＣＯｄ等の計測結果を前記Ｏ2成分のモル数ｆの式に代入して、燃焼反応後のＯ2のモル数ｆを算出する。
【００５８】
以上により、演算部Ｆ２は、燃焼反応後の各成分のモル数ａ〜ｈ、つまり、燃焼直後のガス組成ＣＯｗｃ（＝ａ）、ＣＯ2ｗｃ（＝ｂ）、ＴＨＣｗｃ（＝ｄ）、ＮＯxｗｃ（＝ｇ）、Ｎ2ｗｃ（＝ｈ）、Ｏ2ｗｃ（＝ｆ）、Ｈ2Ｏｗｃ（＝ｃ）、Ｈ2ｗｃ（＝ｅ）を算出する。
算出された燃焼直後のガス組成には、水分濃度Ｈ2Ｏｗｃが含まれている。
【００５９】
このように、ガス計測方法では、排ガス分析計１２で計測したドライベースの排ガス成分濃度ＣＯｄ、ＣＯ2ｄ、ＴＨＣｄ、ＮＯxｄの４成分の計測値と燃料性状とにより、燃焼反応後の各成分ＣＯ、ＣＯ2、ＴＨＣ、ＮＯx、Ｎ2、Ｏ2、Ｈ2、Ｈ2Ｏの割合を算出する。
【００６０】
燃焼直後のガス組成ＣＯｗｃ等を算出した後で、ガス計測方法では、排ガス体積流量算出手段１４によって燃焼直後の排ガス体積流量Ｖｅｘｈ［Ｌ／ｓ］を算出する。
【００６１】
このとき、排ガス体積流量算出手段１４は、まず、演算部Ｆ１を用いて以下の演算処理を行うことにより、燃焼直後の排ガス質量流量Ｇｅｘｈ［ｇ／ｓ］を算出する。
Ｇｅｘｈ［ｇ／ｓ］＝ＧＡ＋ＧＡ／（Ａ／Ｆ）
【００６２】
つまり、演算部Ｆ１は、ＥＣＵに記録される吸入空気量ＧＡ［ｇ／ｓ］と空燃比Ａ／Ｆとから、燃焼直後の排ガス質量流量Ｇｅｘｈ［ｇ／ｓ］を算出する。
【００６３】
排ガス体積流量算出手段１４は、演算部Ｆ３を用いて以下のような演算処理を行うことにより、燃焼直後の排ガス密度ρｅｘｈ［ｇ／Ｌ］を算出する。
【００６４】
演算部Ｆ３は、燃焼直後の排ガス密度ρｅｘｈ［ｇ／Ｌ］の算出に際して、排ガス平均分子量Ｍｅｘｈ［ｇ／ｍｏｌ］（＝モル質量）を算出する。
このとき、演算部Ｆ３は、ガス組成算出手段１３によって算出された燃焼直後のガス組成ＣＯｗｃ等を、以下の式に代入する。
Ｍｅｘｈ＝（Ｎ2ｗｃ×Ｎ2分子量）＋（Ｏ2ｗｃ×Ｏ2分子量）＋（Ｈ2ｗｃ×Ｈ2分子量）＋（Ｈ2Ｏｗｃ×Ｈ2Ｏ分子量）＋（ＣＯｗｃ×ＣＯ分子量）＋（ＣＯ2ｗｃ×ＣＯ2分子量）＋（ＴＨＣｗｃ×燃料分子量）＋（ＮＯxｗｃ×ＮＯx分子量）
これにより、演算部Ｆ３は、排ガス平均分子量Ｍｅｘｈ［ｇ／ｍｏｌ］（＝モル質量）を算出する。
【００６５】
つまり、演算部Ｆ３は、Σ（各ガス濃度×各分子量）から排ガス平均分子量Ｍｅｘｈ［ｇ／ｍｏｌ］を算出する。
【００６６】
演算部Ｆ３は、算出された排ガス平均分子量Ｍｅｘｈ［ｇ／ｍｏｌ］を、以下の式に代入することで燃焼直後の排ガス密度ρｅｘｈ［ｇ／Ｌ］を算出する。
ρｅｘｈ＝Ｍｅｘｈ×Ｐｓｔｄ／Ｒ／Ｔｓｔｄ  （理想気体の状態方程式）
前記理想気体の状態方程式において、Ｐｓｔｄは標準状態の絶対圧力１０１．３２５［ｋＰａ］、Ｔｓｔｄは標準状態の絶対温度２７３．１５［Ｋ］、Ｒは理想気体のガス定数８．３１４４７［ｋＰａ×Ｌ／Ｋ／ｍｏ１］である。
【００６７】
燃焼直後の排ガス質量流量Ｇｅｘｈ［ｇ／ｓ］と燃焼直後の排ガス密度ρｅｘｈ［ｇ／Ｌ］とを算出した後で、排ガス体積流量算出手段１４は、演算部Ｆ４を用いて以下の演算処理を行うことにより、燃焼直後の排ガス体積流量Ｖｅｘｈ［Ｌ／ｓ］を算出する。
Ｖｅｘｈ［Ｌ／ｓ］＝Ｇｅｘｈ［ｇ／ｓ］／ρｅｘｈ［ｇ／Ｌ］
【００６８】
これにより、ガス計測方法は、燃焼直後のガス組成ＣＯｗｃ等、吸入空気量ＧＡ［ｇ／ｓ］、および空燃比Ａ／Ｆに基づいて、実測が困難な燃焼直後の排ガス体積流量Ｖｅｘｈ［Ｌ／ｓ］を高精度に推定できる。
【００６９】
燃焼直後の排ガス体積流量Ｖｅｘｈ［Ｌ／ｓ］を算出した後で、ガス計測方法では、水分濃度算出手段１５によってテールエンド２７における水分濃度Ｈ2Ｏｔｌを算出する。
【００７０】
水分濃度算出手段１５は、演算部Ｆ５を用いて以下のような演算処理を行うことにより、テールエンド２７における水分濃度Ｈ2Ｏｔｌを算出する。
このとき、演算部Ｆ５は、以下の知見に基づいて演算処理を行う。
【００７１】
燃焼直後の水分濃度Ｈ2Ｏｗｃと燃焼直後のその他成分ＯＴｗｃとの総和濃度は１（すなわち、１００％）である。従って、燃焼直後の水分濃度Ｈ2Ｏｗｃとその他成分ＯＴｗｃとの間には、以下の関係が成り立つ。
Ｈ2Ｏｗｃ＋ＯＴｗｃ＝１
【００７２】
前記関係式を変形することで、燃焼直後のその他成分ＯＴｗｃは以下のように表される。
ＯＴｗｃ＝１−Ｈ2Ｏｗｃ
【００７３】
図２に示すように、テールエンド２７において、水分Ｈ2Ｏの分の体積は減少している。一方、テールエンド２７において、その他成分ＯＴｔｌの体積は変化していない。
従って、燃焼直後およびテールエンド２７における排ガス体積流量Ｖｅｘｈ・Ｖｅｘｔｌと、その他成分ＯＴｗｃ・ＯＴｔｌとの間には、以下の関係が成り立つ。
ＯＴｗｃ×Ｖｅｘｈ＝ＯＴｔｌ×Ｖｅｘｔｌ
【００７４】
前記関係式を変形することで、テールエンド２７におけるその他成分ＯＴｔｌは以下のように表される。
ＯＴｔｌ＝Ｖｅｘｈ／Ｖｅｘｔｌ×ＯＴｗｃ
【００７５】
なお、図２のΔＶは結露による流量変化分を示すものである。
【００７６】
テールエンド２７における水分濃度Ｈ2Ｏｔｌとテールエンド２７におけるその他成分ＯＴｔｌとの総和濃度は１（すなわち、１００％）である。従って、テールエンド２７における水分濃度Ｈ2Ｏｔｌとその他成分ＯＴｔｌとの間には、以下の関係が成り立つ。
Ｈ2Ｏｔｌ＋ＯＴｔｌ＝１
【００７７】
前記関係式を変形することで、テールエンド２７における水分濃度Ｈ2Ｏｔｌは以下のように表される。
Ｈ2Ｏｔｌ＝１−ＯＴｔｌ
【００７８】
前記水分濃度Ｈ2Ｏｔｌの式に前記その他成分ＯＴｔｌの式を代入することで、テールエンド２７における水分濃度Ｈ2Ｏｔｌの式は以下のように表される。
Ｈ2Ｏｔｌ＝１−Ｖｅｘｈ／Ｖｅｘｔｌ×ＯＴｗｃ
【００７９】
前記その他成分ＯＴｔｌの式を代入した後の前記水分濃度Ｈ2Ｏｔｌの式に、前記その他成分ＯＴｗｃの式を代入することで、テールエンド２７における水分濃度Ｈ2Ｏｔｌの式は以下のように表される。
Ｈ2Ｏｔｌ＝１−Ｖｅｘｈ／Ｖｅｘｔｌ×（１−Ｈ2Ｏｗｃ）
【００８０】
演算部Ｆ５は、このような水分濃度Ｈ2Ｏｔｌの式に、テールエンド２７における排ガス体積流量Ｖｅｘｔｌ［Ｌ／ｓ］の計測結果と、燃焼直後の排ガス体積流量Ｖｅｘｈ［Ｌ／ｓ］の算出結果と、水分濃度Ｈ2Ｏｗｃの算出結果とを代入する。
演算部Ｆ５は、このようにしてテールエンド２７における水分濃度Ｈ2Ｏｔｌを算出する。
【００８１】
これにより、ガス計測方法は、実測が困難な燃焼直後の排ガス体積流量Ｖｅｘｈ［Ｌ／ｓ］等を用いて、テールエンド２７における水分濃度（燃焼直後からテールエンド２７までの水分濃度変化）Ｈ2Ｏｔｌを高精度に推定できる。
【００８２】
テールエンド２７における水分濃度Ｈ2Ｏｔｌを算出した後で、ガス計測方法では、モーダルマスエミッション算出手段１６によってＣＯ、ＣＯ2、ＴＨＣ、ＮＯxの４成分のモーダルマスエミッションｍＣＯ・ｍＣＯ2・ｍＴＨＣ・ｍＮＯx［ｇ］を算出する。
【００８３】
このとき、モーダルマスエミッション算出手段１６は、演算部Ｆ６を用いて以下のような演算処理を行うことにより、モーダルマスエミッションｍＣＯ・ｍＣＯ2・ｍＴＨＣ・ｍＮＯx［ｇ］を算出する。
【００８４】
演算部Ｆ６は、以下の数式３および数式４のいずれか一方を用いてＣＯのモーダルマスエミッションｍＣＯ［ｇ］を算出する。
【数３】
2015132552.tif 000005
数式３は、ＣＯ濃度を除湿前相当に換算した式である。
【数４】
2015132552.tif 000006
数式４は、排ガス流量を除湿後相当に換算した式である。
【００８５】
なお、前記数式３および数式４において、「∫」に付されたｔｓは計測開始時間であり、ｔｅは計測終了時間である。
また、前記数式３および数式４において法規値とは、法規にて定められた規制対象成分の密度値を指す。
これらは、ＣＯ2、ＴＨＣ、ＮＯxのモーダルマスエミッションｍＣＯ2・ｍＴＨＣ・ｍＮＯx［ｇ］の算出に用いられる式（数式５〜１０）においても同様である。
【００８６】
演算部Ｆ６は、以下の数式５および数式６のいずれか一方を用いてＣＯ2のモーダルマスエミッションｍＣＯ2［ｇ］を算出する。
【数５】
2015132552.tif 000007
数式５は、ＣＯ2濃度を除湿前相当に換算した式である。
【数６】
2015132552.tif 000008
数式６は、排ガス流量を除湿後相当に換算した式である。
【００８７】
演算部Ｆ６は、以下の数式７および数式８のいずれか一方を用いてＴＨＣのモーダルマスエミッションｍＴＨＣ［ｇ］を算出する。
【数７】
2015132552.tif 000009
数式７は、ＴＨＣ濃度を除湿前相当に換算した式である。
【数８】
2015132552.tif 000010
数式８は、排ガス流量を除湿後相当に換算した式である。
【００８８】
演算部Ｆ６は、以下の数式９および数式１０のいずれか一方を用いてＮＯxのモーダルマスエミッションｍＮＯx［ｇ］を算出する。
【数９】
2015132552.tif 000011
数式９は、ＮＯx濃度を除湿前相当に換算した式である。
【数１０】
2015132552.tif 000012
数式１０は、排ガス流量を除湿後相当に換算した式である。
【００８９】
演算部Ｆ６は、ドライベースの排ガス成分濃度ＣＯｄ等、テールエンド２７における排ガス体積流量Ｖｅｘｔｌ［Ｌ／ｓ］、およびテールエンド２７における水分濃度Ｈ2Ｏｔｌを、モーダルマスエミッションｍＣＯ・ｍＣＯ2・ｍＴＨＣ・ｍＮＯx［ｇ］の算出式（前記数式３〜１０）に代入する。
これにより、演算部Ｆ６は、ＣＯ、ＣＯ2、ＴＨＣ、ＮＯxの４成分のモーダルマスエミッションｍＣＯ・ｍＣＯ2・ｍＴＨＣ・ｍＮＯx［ｇ］を算出する。
【００９０】
これによれば、ガス計測方法は、高精度に推定したテールエンド２７における水分濃度（燃焼直後からテールエンド２７までの水分濃度変化）Ｈ2Ｏｔｌを用いて、モーダルマスエミッションｍＣＯ・ｍＣＯ2・ｍＴＨＣ・ｍＮＯx［ｇ］を算出できる。
このため、ガス計測方法は、排ガス管途中の水分濃度変化の影響を受けることなく、高精度にモーダルマスエミッションｍＣＯ・ｍＣＯ2・ｍＴＨＣ・ｍＮＯx［ｇ］を算出できる。
【００９１】
このように、ガス計測方法では、排ガス体積流量Ｖｅｘｔｌ［Ｌ／ｓ］と排ガス成分濃度ＣＯｄ、ＣＯ2ｄ、ＴＨＣｄ、ＮＯxｄとを乗算してモーダルマスエミッション（対象ガスの排出量）を算出する。
【００９２】
モーダルマスエミッションｍＣＯ・ｍＣＯ2・ｍＴＨＣ・ｍＮＯx［ｇ］を算出した後で、ガス計測方法では、誤差補正手段１７によってモーダルマスエミッションｍＣＯ・ｍＣＯ2・ｍＴＨＣ・ｍＮＯx［ｇ］を補正する。
【００９３】
誤差補正手段１７は、ＣＶＳ（Ｃｏｎｓｔａｎｔ  Ｖｏｌｕｍｅ  Ｓａｍｐｌｉｎｇ、定容量採取装置）法を用いてマスエミッション（対象ガスの排出量）を算出した場合に発生する誤差を、演算部Ｆ７によって補正するものである。
【００９４】
以下では、ＣＶＳ法を用いてマスエミッションを算出した場合に発生する誤差について説明する。
【００９５】
図４に示すように、ＣＶＳ法では、希釈用トンネル３１、臨界ノズル３２、希釈空気バッグ３３、サンプルバッグ３４、およびポンプ３５等を用いて希釈した排ガスを計測する。
【００９６】
希釈用トンネル３１は、排ガスの流れ方向上流側端部で希釈空気を導入し、排ガスを希釈空気で希釈する。
排ガスの流れ方向下流側において、希釈用トンネル３１には、臨界ノズル３２が設けられる。
【００９７】
臨界ノズル３２は、中途部に形成されるスロート部３２ａの断面積Ａが他の断面積よりも小さい、すなわち、中途部が絞られた管である。
希釈排ガスは、臨界ノズル３２を通過する際にスロート部３２ａで流れが絞られることにより、流速が臨界流速（その温度における音速）に達する。
これにより、臨界ノズル３２は、希釈排ガスの通過流量を一定に制御する。
【００９８】
希釈空気バッグ３３は、伸縮可能な素材によって製造され、内部に希釈空気を貯溜する。希釈空気バッグ３３は、ポンプおよび配管等を介して、希釈用トンネル３１の希釈空気の導入口近傍と連通する。
希釈空気バッグ３３は、前記ポンプの駆動によって内部に希釈空気を貯溜する。
【００９９】
サンプルバッグ３４は、伸縮可能な素材によって製造され、内部に希釈した排ガスを貯溜する。サンプルバッグ３４は、ポンプおよび配管等を介してスロート部３２ａの、排ガスの流れ方向上流側と連通される。
サンプルバッグ３４は、前記ポンプの駆動によって内部に希釈排ガスを貯溜する。
【０１００】
ポンプ３５は、臨界ノズル３２の排ガスの流れ方向下流側に設けられる。希釈排ガスは、ポンプ３５の駆動によって臨界ノズル３２を通過して、外部に排出される。
【０１０１】
ＣＶＳ法では、例えば、以下の数式１１によって、ＣＯ、ＣＯ2、ＴＨＣ、ＮＯxの４成分のマスエミッションｍＸ０を算出する。
【数１１】
2015132552.tif 000013
前記数式１１において、ρは対象ガス成分密度、Ｃsambagはサンプルバッグ３４内の対象ガス濃度（ｐｐｍ）、Ｃambbagは希釈空気バッグ３３内の対象ガス濃度（ｐｐｍ）、ＤＦは排ガスの希釈率、Ｖmixは希釈排ガス総体積である。
【０１０２】
前記数式１１における希釈排ガス総体積Ｖmixは、臨界ノズル３２を通過した希釈排ガスの全体積であり、以下の数式１２によって算出される。
【数１２】
2015132552.tif 000014
前記数式１２において、Ｃ1は流量係数、Ｐvは希釈用トンネル３１を流れる希釈排ガスの圧力、Ｔvは希釈用トンネル３１を流れる希釈排ガスの温度、ｔは時間である。
【０１０３】
スロート部３２ａの排ガスの流れ方向上流側と下流側との圧力差が臨界値を越えた状態で、希釈排ガス体積流量Ｑｎは、以下の数式１３によって算出される。
【数１３】
2015132552.tif 000015
前記数式１３において、Ａはスロート部３２ａの断面積、κは比熱比、Ｍは気体分子量、Ｐvは希釈用トンネル３１を流れる希釈排ガスの圧力、Ｒは理想気体のガス定数（８．３１４４７［ｋＰａ×Ｌ／Ｋ／ｍｏ１］）、Ｔvは希釈用トンネル３１を流れる希釈排ガスの温度、ρnは標準密度である。
【０１０４】
流量が一定値に制御される臨界ノズル３２において、臨界ノズル３２を流れる流体の体積は、体積流量と時間との積によって算出されることとなる。従って、前記数式１２において、希釈排ガス体積流量Ｑnは、以下の数式１４によって表される。
【数１４】
2015132552.tif 000016
【０１０５】
数式１４を数式１３に代入するとともに、理想気体の状態方程式（ρn＝Ｍ×Ｐ0／Ｒ／Ｔ0）を標準密度ρnに代入して整理することで、流量係数Ｃ1は、以下の数式１５によって表される。
【数１５】
2015132552.tif 000017
前記数式１５において、Ｔ0は標準温度（２０［℃］、すなわち２９３．１５［Ｋ］）、Ｐ0は標準圧力（１［ａｔｍ］、すなわち、１０１．３２５［ｋＰａ］）である。
【０１０６】
前記数式１５において、スロート部３２ａの断面積Ａ、理想気体のガス定数Ｒ、標準温度Ｔ0、標準圧力Ｐ0は定数である。
一方、比熱比κおよび気体分子量Ｍは、物性値であり、流れる気体に依存する変数である。
【０１０７】
ＣＶＳ法において、流量係数Ｃ1の校正は、空気を媒体に行われる、つまり、ＣＶＳ法においては、前記数式１５の比熱比κおよび気体分子量Ｍに、空気の比熱比および空気の気体分子量が代入されることで、流量係数Ｃ1が算出される。
【０１０８】
このように、ＣＶＳ法において、流量係数Ｃ1は、空気の流量係数である。
【０１０９】
一方、希釈用トンネル３１を流れるのは、希釈排ガスである。
つまり、ＣＶＳ法における誤差要因は、前記数式１５における比熱比κおよび気体分子量Ｍである。
【０１１０】
そこで、図３に示すように、ガス計測方法では、誤差補正手段１７の演算部Ｆ７を用いて、希釈排ガス成分の流量係数Ｃmを算出する。ガス計測方法では、以下のような演算処理を行うことによって希釈排ガス成分の流量係数Ｃmを算出する。
希釈排ガス成分の流量係数Ｃmは、数式１５における比熱比κおよび気体分子量Ｍに代えて、希釈排ガスの比熱比κ´および気体分子量Ｍ´を用いた以下の数式１６によって算出される。
【数１６】
2015132552.tif 000018
【０１１１】
演算部Ｆ７は、希釈排ガスの比熱比κ´を以下の数式１７によって算出する。
【数１７】
2015132552.tif 000019
前記数式１７において、比熱比κEGは希釈前の排ガスの比熱比、κairは空気の比熱比である。
【０１１２】
演算部Ｆ７は、２原子分子と３原子分子以上の多原子分子の比率によって、希釈前の排ガス比熱比κEGの値を決定する。
【０１１３】
演算部Ｆ７は、例えば、演算部Ｆ２で算出した燃焼直後のガス組成ＣＯｗｃ等から、排ガスにおける２原子分子（すなわち、Ｎ2等）の割合と、原子分子以上の多原子分子（すなわち、ＣＯ2等）の割合とを算出する。
そして、演算部Ｆ７は、算出したそれぞれの割合に、予め設定される２原子分子の比熱比および原子分子以上の多原子分子の比熱比をそれぞれ乗算する。演算部Ｆ７は、２原子分子および原子分子以上の多原子分子の乗算結果を合算することで、希釈前の排ガス比熱比κEGを算出する。
【０１１４】
なお、希釈排ガスの比熱比κ´の算出において、前記数式１７の希釈率ＤＦは、ＣＶＳ法で予め想定される値（定数）が代入される。
【０１１５】
このように、演算部Ｆ７は、希釈前の排ガスおよび空気の比熱比κEG・κairに、希釈率ＤＦ（すなわち、希釈排ガスにおける排ガスおよび空気の比率）をそれぞれ反映することで希釈排ガスの比熱比κ´を算出する。
【０１１６】
演算部Ｆ７は、希釈排ガスの気体分子量Ｍ´を以下の数式１８によって算出する。
【数１８】
2015132552.tif 000020
前記数式１８において、ＭEGは２０［℃］（２９３．１５［Ｋ］）１［ａｔｍ］（１０１．３２５［ｋＰａ］）での希釈前の排ガス平均分子量、Ｍairは２０℃１ａｔｍでの空気の平均分子量である。
【０１１７】
演算部Ｆ７は、２０℃１ａｔｍでの排ガス平均分子量ＭEGとして、演算部Ｆ３で算出された排ガス平均分子量Ｍｅｘｈ［ｇ／ｍｏｌ］（＝モル質量）を用いる。
つまり、演算部Ｆ７は、前記数式１８の希釈前の排ガス平均分子量ＭEGに排ガス平均分子量Ｍｅｘｈ［ｇ／ｍｏｌ］の算出結果を代入するとともに、既知である２０℃１ａｔｍでの空気の平均分子量Ｍairを代入することで、希釈排ガスの気体分子量Ｍ´を算出する。
【０１１８】
なお、排ガス平均分子量ＭEGの算出において、前記数式１８の希釈率ＤＦは、ＣＶＳ法で予め想定される値が代入される。
【０１１９】
このように、演算部Ｆ７は、希釈前の排ガスおよび空気の平均分子量ＭEG・Ｍairに、希釈率ＤＦをそれぞれ反映することで希釈排ガスの気体分子量Ｍ´を算出する。
【０１２０】
以上のようにガス計測方法では、臨界ノズル３２における排ガス（ＣＶＳ法における希釈排ガス）の気体分子量Ｍ´および比熱比κ´を算出する工程を行う。
【０１２１】
演算部Ｆ７は、希釈排ガスの比熱比κ´および気体分子量Ｍ´の算出結果を前記数式１６に代入することにより、希釈排ガス成分の流量係数Ｃmを算出する。
【０１２２】
このように、ガス計測方法では、排ガス（ＣＶＳ法における希釈排ガス）の気体分子量Ｍ´および比熱比κ´の算出結果を用いて、臨界ノズル３２における排ガス成分の流量係数Ｃmを算出する工程を行う。
【０１２３】
演算部Ｆ７は、希釈排ガス成分の流量係数Ｃmを以下の式に代入することにより、演算部Ｆ６で算出したＣＯ、ＣＯ2、ＴＨＣ、ＮＯxの４成分のモーダルマスエミッションｍＸを補正する。
ｍＸ´＝ｍＸ×Ｃ1／Ｃm
前記式において、ｍＸ´は補正したモーダルマスエミッションｍＣＯ´・ｍＣＯ2´・ｍＴＨＣ´・ｍＮＯx´である。
【０１２４】
これによれば、ガス計測方法は、演算部Ｆ６で算出したモーダルマスエミッションｍＸに、ＣＶＳ法によってマスエミッションｍＸ０を算出した場合に発生する誤差を反映できる。
従って、ガス計測方法は、排ガスの認証試験等において標準的に用いられるＣＶＳ法でのマスエミッションｍＸ０の算出結果に対して、車両の開発段階で用いられるモーダルマス法でのモーダルマスエミッションｍＸ´の算出結果を近付けることができる。
【０１２５】
つまり、ガス計測方法は、モーダルマスエミッションｍＸ´の算出結果が、ＣＶＳ法でのマスエミッションｍＸ０の算出結果に対して乖離してしまうことを防止できる。
従って、ガス計測方法は、排ガスの認証試験において、車両の開発段階で得られた結果と同じ結果（狙い通りの結果）を得ることができる。
【０１２６】
このように、ガス計測方法では、演算部Ｆ６で算出したモーダルマスエミッションｍＸに対して、排ガス成分の流量係数Ｃmと空気の流量係数Ｃ1との比（Ｃ1／Ｃm）を乗算することで、モーダルマスエミッションｍＸを補正する工程を行う。
【符号の説明】
【０１２７】
２３        エンジン
３２        臨界ノズル
Ｃ1        空気の流量係数
Ｃm        排ガス成分の流量係数
ＣＯｄ      ＣＯ成分濃度（排ガス成分濃度）
ＣＯ2ｄ    ＣＯ2成分濃度（排ガス成分濃度）
Ｍ´        排ガスの気体分子量
ｍＸ        モーダルマスエミッション（対象ガスの排出量）
ｍＸ´      補正したモーダルマスエミッション（補正した対象ガスの排出量）
ＮＯxｄ    ＮＯx成分濃度（排ガス成分濃度）
ＴＨＣｄ    ＴＨＣ成分濃度（排ガス成分濃度）
Ｖｅｘｔｌ  排ガス体積流量
κ´        排ガスの比熱比
【図１】
2015132552.tif 000021
【図２】
2015132552.tif 000022
【図３】
2015132552.tif 000023
【図４】
2015132552.tif 000024