(57)【特許請求の範囲】
【請求項１】
コンピュータによって、燃焼装置からの排ガスに尿素水を添加して脱硝触媒で窒素酸化物の分解を行う脱硝システムにおける尿素の分解反応に伴うパラメータの時間変化を算出するプログラムであって、
尿素に関する尿素パラメータの初期値として前記尿素水の液滴温度、液滴径、液滴速度及び尿素濃度の設定を受け付けるステップ１と、
前記排ガスに関する排ガスパラメータとして前記排ガスの流速及び温度の設定を受け付けるステップ２と、
前記尿素パラメータと前記排ガスパラメータを用いて前記排ガスへの添加後の時間の経過に伴う前記尿素水の前記液滴径の変化と前記液滴温度の変化とを液滴の時間変化を求めるＳｐａｌｄｉｎｇの式に基づいて算出するステップ３と、
水分の蒸発時による前記尿素水の前記尿素濃度の変化を前記液滴径の変化と前記液滴温度の変化に基づいて算出するステップ４と、
前記尿素濃度と等価である前記尿素水の水分が、前記液滴内の前記尿素が分解反応を起こし始める前記液滴内の水分の下限値である水分閾値以下となったかを判定するステップ５と、
前記脱硝触媒の有無による反応速度定数を求める熱分解に伴うイソシアン酸とアンモニアを生成するＹｉｍの実験式に基づいて前記尿素濃度の時間変化を求めるステップ６と、
前記尿素濃度が、前記尿素の分解反応の演算を終了する前記尿素濃度の下限値である濃度閾値以下となったかを判定するステップ７と、
前記ステップ３、前記ステップ４及び前記ステップ６の少なくとも１つにおける算出結果を出力するステップ８と、
を備えることを特徴とする脱硝システム用尿素水分解解析プログラム。
【請求項２】
請求項１に記載の脱硝システム用尿素水分解解析プログラムであって、
前記ステップ６における前記尿素の熱分解反応では、尿素濃度とイソシアン酸濃度とアンモニア濃度を算出することを特徴とする脱硝システム用尿素水分解解析プログラム。
【請求項３】
請求項２に記載の脱硝システム用尿素水分解解析プログラムであって、
前記ステップ８は、前記時間に対する前記液滴径、前記液滴温度、前記尿素濃度、前記イソシアン酸濃度、及び前記アンモニア濃度の少なくとも１つを出力することを特徴とする脱硝システム用尿素水分解解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、脱硝システム用尿素水分解解析プログラムに関する。
【背景技術】
【０００２】
国際海事機関（ＩＭＯ）による窒素酸化物（ＮＯｘ）の三次規制に対する対策の一つとして、尿素水を用いた選択触媒還元（ＳＣＲ：Selective Catalytic Reduction）の研究開発が行われている。このＳＣＲによる脱硝では、通常以下の反応が触媒中で起こっていることが知られている。
【化１】
0005935247.tif 000002
【０００３】
この反応には還元剤としてアンモニア（ＮＨ３）が必要になるが、アンモニア（ＮＨ３）の代わりにアンモニア（ＮＨ３）の類似化合物として尿素水を使用する手法が、安全面や貯蓄、輸送面の理由から広く採用されている。
【０００４】
一方、三次規制は指定された海域内でのみ適用され、それ以外の海域では二次規制が適用されるが、この海域は基本的に港湾や輻輳海域等の船舶交通量の比較的大きく、陸地に近い場所に設定されると考えられる。その時、ＳＣＲはこの海域を航行する限られた時間のみで運用され、航行中の大部分はＳＣＲを作動させないことが想定できる。そうなると、ＳＣＲの運用に関する自由度は上がり、例えば劣化した触媒の再生をＮＯｘ還元の必要のない海域航行中に行うようなことも可能になるためＳＣＲの寿命や劣化再生方法について調べることが重要となってくる。しかし、これらの検討の手段として実験を行おうとすると非常に長い時間がかかることや温度や排ガス成分等のパラメータが非常に多いことから、数値解析による予測が有効な手段になると考えられる。
【０００５】
ＳＣＲを用いた脱硝反応に関して基礎的な実験・解析から数値解析まで，数多くの研究が行われている。例えば、窒素酸化物を含有する排ガス中に還元剤であるアンモニアを投入する際の制御について開示されている（特許文献１，２等）。特許文献１では、触媒ユニット、エンジンの作動に関するパラメータを用いて制御が行われる。特許文献２では、窒素酸化物（ＮＯｘ）浄化特性の多項式を用いてアンモニア吸着量の目標値を設定する方法が採用されている。また、脱硝システムにおける脱硝触媒の劣化を診断する診断装置が開示されている（特許文献３）。ここでは、アンモニアが漏れるまで脱硝触媒へ供給された還元剤の総量に基づいて脱硝触媒が劣化したか否かを判定している。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特表平０８−５０９７９５号公報
【特許文献２】特開２００８−１９６３４０号公報
【特許文献３】特開２００９−１２７４９６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
ところで、脱硝触媒の劣化又は再生を予測する計算において、上記の尿素の分解反応を正確に予測することが望まれる。すなわち、還元剤として尿素水を用いる場合、触媒での脱硝反応の前までに尿素水から蒸発、分解によって、アンモニア（ＮＨ３）を生成し、想定している脱硝反応を行えるようにすることが必要である。さらに、どのような条件でアンモニア（ＮＨ３）の供給量を満たせるかを予測することが、触媒容量や噴射ノズルの設置等を決定する上で重要になる。
【０００８】
本発明は、上記の課題を鑑み、尿素の分解反応によるアンモニアの生成量を予測することを可能とする脱硝システム用尿素水分解解析プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
請求項１に対応した脱硝システム用尿素水分解解析プログラムは、コンピュータによって、燃焼装置からの排ガスに尿素水を添加して脱硝触媒で窒素酸化物の分解を行う脱硝システムにおける尿素の分解反応に伴うパラメータの時間変化を算出するプログラムであって、尿素に関する尿素パラメータの初期値として前記尿素水の液滴温度、液滴径、液滴速度及び尿素濃度の設定を受け付けるステップ１と、前記排ガスに関する排ガスパラメータとして前記排ガスの流速及び温度の設定を受け付けるステップ２と、前記尿素パラメータと前記排ガスパラメータを用いて前記排ガスへの添加後の時間の経過に伴う前記尿素水の前記液滴径の変化と前記液滴温度の変化とを液滴の時間変化を求めるＳｐａｌｄｉｎｇの式に基づいて算出するステップ３と、水分の蒸発時による前記尿素水の前記尿素濃度の変化を前記液滴径の変化と前記液滴温度の変化に基づいて算出するステップ４と、前記尿素濃度と等価である前記尿素水の水分が、前記液滴内の前記尿素が分解反応を起こし始める前記液滴内の水分の下限値である水分閾値以下となったかを判定するステップ５と、前記脱硝触媒の有無による反応速度定数を求める熱分解に伴うイソシアン酸とアンモニアを生成するＹｉｍの実験式に基づいて前記尿素濃度の時間変化を求めるステップ６と、前記尿素濃度が、前記尿素の分解反応の演算を終了する前記尿素濃度の下限値である濃度閾値以下となったかを判定するステップ７と、前記ステップ３、前記ステップ４及び前記ステップ６の少なくとも１つにおける算出結果を出力するステップ８と、を備えることを特徴とする。
【００１２】
また、前記ステップ６における前記尿素の熱分解反応では、尿素濃度とイソシアン酸濃度とアンモニア濃度を算出することが好ましい。
【００１３】
また、前記ステップ８は、前記時間に対する前記液滴径、前記液滴温度、前記尿素濃度、前記イソシアン酸濃度、及び前記アンモニア濃度の少なくとも１つを出力することが好ましい。
【発明の効果】
【００１４】
本発明の脱硝システム用尿素水分解解析プログラムでは、尿素の分解反応に伴うパラメータの時間変化を正確に予測することが可能となる。これによって、脱硝触媒の劣化又は再生を予測する計算において還元剤として尿素水を用いる場合、触媒での脱硝反応の前までに尿素水からアンモニア（ＮＨ３）がどのように生成され、想定している脱硝反応がどのように行われるかを正確に予測することが可能となる。さらに、尿素に関するパラメータや排ガスに関するパラメータが変わった場合、どのような条件でアンモニア（ＮＨ３）の供給量を満たせるかを予測することで、触媒容量や噴射ノズルの設置等を決定することが可能となる。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施の形態におけるコンピュータシステムの構成例を示す図である。
【図２】本発明の実施の形態における脱硝システム用尿素水分解解析処理のフローチャートを示す図である。
【図３】本発明の実施の形態における蒸発による液滴径の変化を示す図である。
【図４】本発明の実施の形態における蒸発による液滴内の温度の変化を示す図である。
【図５】本発明の実施の形態における尿素水濃度による液滴寿命の変化を示す図である。
【図６】本発明の実施の形態における液滴速度の時間変化を示す図である。
【図７】本発明の実施の形態における還元剤による触媒性能及び劣化性能の比較を示す図である。
【図８】本発明の実施の形態における尿素の分解割合による脱硝性能の変化を示す図である。
【図９】本発明の実施の形態における尿素熱分解割合変更時の触媒中のＮＯ，ＮＨ３濃度の変化を示す図である。
【発明を実施するための形態】
【００１６】
＜脱硝システム用尿素水分解解析処理＞
本発明の実施の形態における脱硝システム用尿素水分解解析処理は、コンピュータシステム１００を用いて実行される。コンピュータシステム１００は、図１に示すように、処理部１０、記憶部１２、入力部１４及び出力部１６を含んで構成される。
【００１７】
処理部１０は、ＣＰＵ等で構成され、記憶部１２に予め記憶された脱硝システム用尿素水分解解析プログラムを実行することにより脱硝システム用尿素水分解解析処理の演算を行う。記憶部１２は、脱硝システム用尿素水分解解析プログラム及びそれに関する各種パラメータを格納及び保持する。記憶部１２は、半導体メモリ、ハードディスク等の記憶装置を含んで構成され、処理部１０から読み出し及び書き込み可能とされる。入力部１４は、脱硝システム用尿素水分解解析処理に必要なデータ等を入力するために用いられる。入力部１４は、キーボード、ポインティングデバイス等の入力装置を含んで構成される。出力部１６は、入力部１４から入力されたデータの表示、脱硝システム用尿素水分解解析処理で得られる中間処理結果や最終処理結果の表示等を行う。出力部１６は、ディスプレイ、プリンタ等の出力装置を含んで構成される。
【００１８】
以下、本発明の実施の形態における脱硝システム用尿素水分解解析処理について説明する。脱硝システム用尿素水分解解析処理は、処理部１０によって記憶部１２に記憶されている脱硝システム用尿素水分解解析プログラムを実行することによって、図２に示すフローチャートに沿って行われる。
【００１９】
なお、脱硝システム用尿素水分解解析処理は以下の前提条件下において行われる。（前提１）尿素水液滴の周りにあるガスは、排気ガスではなく、空気である。（前提２）尿素水液滴から発生するのは水のみである。（前提３）尿素水液滴は、周囲にある空気と直接情報（物理量：質量、熱）のやり取りをせず、「フィルム」と称する中間体が尿素水と空気を仲介するものとする。
【００２０】
ステップＳ１０では、脱硝システム用尿素水分解解析処理で使用される定数の設定が行われる。処理部１０は、入力部１４から定数の入力を受け、記憶部１２に記憶させる。ここで入力される定数は、水及び空気のガス定数、水の臨界定数、水、空気、尿素の分子量等が挙げられる。
【００２１】
ステップＳ１２では、尿素（尿素水）に関するパラメータについて初期設定が行われる。処理部１０は、入力部１４からパラメータの入力を受け、記憶部１２に記憶させる。ここで、尿素に関するパラメータは、液滴温度Ｔｐ、液滴直径ｄｐ、液滴速度ｕｐ、尿素水滴内尿素濃度Ｗｔ（尿素水１ｋｇあたりの尿素質量）が挙げられる。これら尿素（尿素水）のパラメータ設定は、尿素の含有形態、供給方法、パラメータの表現方法等により各種のパラメータや組み合わせがあり得る。
【００２２】
また、処理部１０は、これらの設定値から尿素水液滴内の尿素モルｎ、水モル数Ｎの比ｎ／Ｎを式（１）により、尿素水密度ρｄを式（２）により、尿素水液滴質量ｍを式（３）、尿素水液滴内尿素及び水質量を式（４）により算出する。ただし、ｆ１は、尿素水液滴温度及びｎ／Ｎを引数として実験結果から多項式にて尿素水密度ρｄを算出する関数である。
【数１】
0005935247.tif 000003
【数２】
0005935247.tif 000004
【数３】
0005935247.tif 000005
【数４】
0005935247.tif 000006
【００２３】
ステップＳ１４では、排ガスに関するパラメータについて設定が行われる。処理部１０は、入力部１４から定数の入力を受け、記憶部１２に記憶させる。ここで入力される排ガスに関するパラメータは、排ガスはほぼ空気に近いものと考えて、空気の流速ｕａ、温度Ｔ∞、圧力Ｐ等が挙げられる。これら、排ガスのパラメータ設定は、排ガスそのものの流速を用いたり、過渡的な温度も考慮したりすることも可能である。また、パラメータの表現方法やパラメータの組み合わせは各種の方法があり得る。
【００２４】
ステップＳ１６では、現時刻ｔにおける各種の物性値が算出される。以下の処理は、時刻ｔをステップΔｔ毎に進めながら、各時刻ｔにおける物性値を算出することによって行われる。まず、処理部１０は、時刻ｔに対して、式（１）に基づいてｎ／Ｎを算出する。次に、フィルム内の水の温度、質量分率は以下の１／３ルールに従うと仮定して、フィルム内の温度Ｔｆを算出する。
【数５】
0005935247.tif 000007
【００２５】
フィルム内の空気の比熱Ｃｐ＿ａｉｒ、粘性係数μ＿ａｉｒ、熱伝導率λ＿ａｉｒ、密度ρ＿ａｉｒを算出する。ここで、関数ｆ２，ｆ３，ｆ４は、実験結果から得られたフィルム内の温度Ｔｆを引数とする多項式の関数である。
【数６】
0005935247.tif 000008
【００２６】
また、フィルム内の水（気体）の比熱Ｃｐ＿Ｈ２Ｏ、粘性係数μ＿Ｈ２Ｏ、熱伝導率λ＿Ｈ２Ｏ、密度ρ＿Ｈ２Ｏを算出する。ここで、関数ｆ５，ｆ６，ｆ７は、実験結果から得られたフィルム内の温度Ｔｆを引数とする多項式の関数である。
【数７】
0005935247.tif 000009
【００２７】
次に、水と空気の拡散係数Ｄ［ｍ２／ｓ］を算出する。水と窒素の拡散係数Ｄ０＝０．２２２×１０−４ｍ２／ｓ（Ｐ０＝０．１ＭＰａ、Ｔ０＝２７３．１５Ｋ）から、式（８）で温度と圧力の関数として算出する。ここで、温度Ｔはフィルム内温度Ｔｆと等しく設定する。ｍは定数であり、１．７５とした。
【数８】
0005935247.tif 000010
【００２８】
水蒸発潜熱ＨＬを算出する。水の蒸発潜熱ＨＬ［Ｊ／ｋｇ］はワトソン（Ｗａｔｓｏｎ）の式（９）により算出することができる。ここで、Ｃは１ａｔｍ，３７３．１５Ｋでの蒸発潜熱２．２５６Ｊ／ｋｇからＣ＝３．１２８×１０６とする。Ｔｃは水の臨界温度［Ｋ］である。温度Ｔは、尿素水液滴温度Ｔｐと等しいとする。
【数９】
0005935247.tif 000011
【００２９】
尿素水蒸気圧ｐｕを式（１０），（１１），（１２）から算出する。Ａ，Ｂ，Ｃ，Ｄは定数であり、それぞれ−７．７６４５１，１．４５８３８、−２．７７５８０、−１．２３３０３、Ｐｃは水の臨界圧力である。上記のように、ｎは液滴内の尿素モル数、Ｎは液滴内の水モル数である。ｐｗは水の蒸気圧であり、ワンガー（Ｗａｎｇｅｒ）の式で表される。
【数１０】
0005935247.tif 000012
【数１１】
0005935247.tif 000013
【数１２】
0005935247.tif 000014
【００３０】
次に、尿素水密度ρｄを算出する。ここで、ｆ８は、尿素水液滴温度及びｎ／Ｎの関数である。
【数１３】
0005935247.tif 000015
【００３１】
尿素水（液体）の比熱Ｃｐｕを式（１４）により算出する。Ｄｕは希釈度（溶質１ｍｏｌに対する水のモル数＝Ｎ／ｎ）、Ｍｕは尿素の分子量である。但し、実測値がＤｕ＝５．５までしかないため、Ｄｕ＜５．５ではＤｕ＝５．５の時の値を用いる。尿素水の密度は、計測値から温度とｎ／Ｎの関数として求めた。
【数１４】
0005935247.tif 000016
【００３２】
なお、水（液体）の比熱Ｃｐｗは式（１５）によって算出した。ここで、ｆ９は、液滴温度Ｔｐを引数とする多項式の関数である。
【数１５】
0005935247.tif 000017
【００３３】
フィルム内の質量分率Ｙｗｆは式（１６）により算出される。ここで、尿素水蒸気圧ｐｕは式（１０）で算出された値を用いる。また、周囲空気中の水の質量分率ＹＷ∞は０とする。また、液滴表面の水の質量分率ＹＷｓは以下の式で算出される値とする。
【数１６】
0005935247.tif 000018
ＹＷｓ＝（ｐｕ×水分子量）／（（Ｐ−ｐｕ）×空気分子量＋ｐｕ×水分子量）
【００３４】
トランスファ数ＢＭは数（１７）により算出される。トランスファ数ＢＭは、後述の尿素水からの水の蒸発過程の算出処理に用いられる。
【数１７】
0005935247.tif 000019
【００３５】
また、気液相対速度Ｕｒｅｌを式（１８）から算出する。気液相対速度Ｕｒｅｌは、液滴の運動を運動方程式から求める際に用いられる。ここで、気液相対速度Ｕｒｅｌは、空気の速度ｕａと液滴の速度ｕｐとの差で表される。
【数１８】
0005935247.tif 000020
【００３６】
フィルム内の混合ガス（水＋空気）の粘性係数μｆ、熱伝導率λｆ、比熱Ｃｐｆ及び密度ρｆを式（１９），（２０），（２１），（２２）から算出する。フィルム内は水と空気の混合ガスになるため、粘性係数μｆ、熱伝導率λｆ、比熱Ｃｐｆ及び密度ρｆは成分を考慮した物性値を用いる。粘性係数μｆは、ウィルキー（Ｗｉｌｋｅ）の式により算出される。熱伝導率は、ワシルジェワ（Ｗａｓｓｉｌｊｅｗａ）の式により算出される。
【数１９】
0005935247.tif 000021
【数２０】
0005935247.tif 000022
【数２１】
0005935247.tif 000023
【数２２】
0005935247.tif 000024
【００３７】
ここで、フィルムは水と空気のみによって構成されているとして、添字１が水、添字２が空気とすると式（１９）のφｉｊはφ１２，φ２１のみとなり、以下の式で算出される。
【数２３】
0005935247.tif 000025
φ12=(1+(μ_H2O/μ_air)1/2(空気分子量/水分子量)1/4)2/(2×21/2(1+空気分子量/水分子量)1/2)
φ21=(1+(μ_air/μ_H2O)1/2(水分子量/空気分子量)1/4)2/(2×21/2(1+水分子量/空気分子量)1/2)
【００３８】
また、フィルム内のモル分率の比Ｘ２／Ｘ１，Ｘ１／Ｘ２は以下の式（２４）で表される。
【数２４】
0005935247.tif 000026
【００３９】
式（２３）及び（２４）から式（１９）は以下の式（２５）のように書き替えられる。
【数２５】
0005935247.tif 000027
【００４０】
また、式（２０）におけるΨijは、添字１が水、添字２が空気とすると式（２６）のようにリンゼイ−ブロンリィ（Ｌｉｎｄｓａｙ−Ｂｒｏｍｌｅｙ）の式で表される。Ｃｉはサザーランド（Ｓｕｔｈｅｒｌａｎｄ）定数、Ｃｉｊは混合気体のサザーランド（Ｓｕｔｈｅｒｌａｎｄ）定数、Ｘｉはモル分率である。本研究の場合、考慮する成分は空気と水であるため、サザーランド（Ｓｕｔｈｅｒｌａｎｄ）定数はそれぞれ１０２，６５９とする。
【数２６】
0005935247.tif 000028
ここで、ΨｉｊはΨ１２，Ψ２１のみとなり、以下の式で算出される。
Ψ12=(1/4)(1+((μ_H2O/μ_air)(空気分子量/水分子量)3/4(1+C1/Tf)/(1+C2/Tf))1/2)2
(1+C12/Tf)/(1+C1/Tf)
Ψ21=(1/4)(1+((μ_air/μ_H2O)(水分子量/空気分子量)3/4(1+C2/Tf)/(1+C1/Tf))1/2)2
(1+C21/Tf)/(1+C2/Tf) 
【００４１】
式（２６）から式（２０）は以下の式（２７）のように書き替えられる。
【数２７】
0005935247.tif 000029
【００４２】
また、比熱Ｃｐｆ及び密度ρｆは式（２８），（２９）として書き替えられる。
【数２８】
0005935247.tif 000030
【数２９】
0005935247.tif 000031
【００４３】
さらに、無次元数である以下の値を算出する。
【数３０】
0005935247.tif 000032
【００４４】
これらの値を用いて、抗力係数ＣＤを算出する。
【数３１】
0005935247.tif 000033
【００４５】
ステップＳ１８では、現時刻ｔにおける蒸発による液滴の直径の変化と液滴の温度の変化が算出される。処理部１０は、上記において算出された値を用いてＳｐａｌｄｉｎｇの式（３２）の関係から液滴の直径の変化を算出する。
【数３２】
0005935247.tif 000034
【００４６】
式（３２）を変形すると式（３３）として書き替えることができる。ここで、拡散係数Ｄｆは式（８）、フィルム内の混合ガス密度ρｆは式（２２）で表されるものとし、微小時間での温度変化は０とする。なお、液滴の直径の変化の算出は、液滴の半径に基づくものであってもよく、また、Ｓｐａｌｄｉｎｇの式（３２）に限定されるものでなく、パラメータの設定や表現方法、簡略式や厳密式等、各種の式に基づくことができる。
【数３３】
0005935247.tif 000035
【００４７】
また、液滴の温度Ｔｐの変化は式（３４）で算出される。なお、液滴の温度Ｔｐの変化の算出は、この式（３４）に限定されるものでなく、パラメータの設定や表現方法、簡略式や厳密式等、各種の式に基づくことができる。
【数３４】
0005935247.tif 000036
【００４８】
ステップＳ２０では、蒸発による液滴内の尿素濃度が算出される。処理部１０は、上記ステップＳ１８で算出された液滴の直径の変化dｄｐ／dｔを式（３３）によりｄｍ／ｄｔに変更する。これらから、時刻のステップΔｔ後の液滴質量は式（３５）で算出することができる。
【数３５】
0005935247.tif 000037
【００４９】
この質量変化は前提２より水の蒸発のみで生じているものとし、式（３６）で尿素水の濃度Ｗｔを算出する。なお、尿素濃度の算出は、式（３６）に限定されるものでなく、パラメータの設定や表現方法、また、補正係数を掛ける等による各種の式に基づくことができる。
【数３６】
0005935247.tif 000038
【００５０】
ステップＳ２２では、液滴速度が算出される。処理部１０は、式（３７）から液滴の速度ｕｐを更新する。ここで、雰囲気速度ｕａ、相対速度Ｕｒｅｌ、抗力係数ＣＤを適用しており、抗力係数ＣＤはレイノルズ数Ｒｅｄに応じて式（３１）により算出して用いている。
【数３７】
0005935247.tif 000039
【００５１】
ステップＳ２４では、液滴内の水分が所定の基準値である水分閾値ＭＩＮ１以下となったか否かが判定される。水分閾値ＭＩＮ１は、液滴内の水分が減って、液滴内の尿素が分解反応を起こし始める液滴内の水分の下限値である。液滴内の水分が水分閾値ＭＩＮ１以下である場合にはステップＳ２６へ処理を移行させ、そうでない場合には時刻ｔをステップΔｔだけ進めてステップＳ１６から処理を繰り返す。なお、水分が水分閾値以下となったか否かの判定は、例えば尿素濃度をもっても等価に判定できるものであり、実態として尿素の分解反応の開始が推定できる他の方法を含めることができる。
【００５２】
ステップＳ２６では、排ガス配管内を流れている空気に含まれている尿素、イソシアン酸、水及びアンモニアの空間内の濃度の初期値を設定する。処理部１０は、入力部１４から尿素濃度Ｃｕ、イソシアン酸ＣＨＮＣＯ、水ＣＨ２Ｏ及びアンモニアＣＮＨ３の初期値を受付け、記憶部１２に記憶させる。
【００５３】
ステップＳ２８では、尿素の減少量、イソシアン酸の生成量及びアンモニアの生成量が算出される。ＳＣＲの還元剤に用いられている尿素水はアンモニア（ＮＨ３）に分解されて一酸化窒素（ＮＯ）又は二酸化窒素（ＮＯ２）と触媒上で反応するが、尿素水からアンモニア（ＮＨ３）の生成は総括反応では化学式（２）のように示される。
【化２】
0005935247.tif 000040
【００５４】
ここで、ａｑは水溶液を表している。この反応は化学式（３）〜化学式（５）で表される尿素水（液滴）から水の蒸発、尿素の熱分解及びイソシアン酸の加水分解の３つの現象に分けられることが知られている。ここで、ｌは液体、ｇは気体を表している。なお、イソシアン酸（ＨＮＣＯ）は気相では安定な物質であり、式（５）の加水分解は触媒表面で起こる反応であると考えられている。
【化３】
0005935247.tif 000041
【化４】
0005935247.tif 000042
【化５】
0005935247.tif 000043
【００５５】
処理部１０は、式（３８），（３９），（４０）により尿素の減少量、イソシアン酸の生成量及びアンモニアの生成量を求める。本実施の形態では、化学式（３）〜（５）で示されている尿素水の分解過程は、以下のような条件下で進行していくとする。（条件１）尿素の熱分解はＹｉｍの式に従いイソシアン酸とアンモニア（ＮＨ３）を生成する。なお、熱分解は触媒中でも起こる。（条件２）イソシアン酸の加水分解はＹｉｍの式に従い、触媒表面上でのみ反応する。
【数３８】
0005935247.tif 000044
【数３９】
0005935247.tif 000045
【数４０】
0005935247.tif 000046
【００５６】
ただし、ｋ１，ｋ２は反応速度定数［ｌ／ｓ］である。Ｙｉｍらは、触媒の有無によるｋ１とｋ２を求めているが、ここではｋ１には触媒なしでの、ｋ２には触媒有りでの結果を用いる。ｋ１に関しては触媒の有無による差がほとんどなかったことから、触媒中でも触媒なしでの式（４１），（４２）を用いる。ここでＲ＝１．９８６ｋｃａｌ／ｋｍｏｌＫとする。なお、式（４１），（４２）はＣｕＺＳＭ５触媒の結果であり、今回対象としているＶ２Ｏ５−ＴｉＯ２系触媒の結果ではないが、Ｖ２Ｏ５−ＴｉＯ２系触媒での活性化エネルギとほとんど変わらないことから、本実施の形態では式（４１），（４２）を採用した。なお、式（４１），（４２）を用いずに、Ｖ２Ｏ５−ＴｉＯ２系触媒の結果を用いることもできる。また、尿素の熱分解反応を求める方法は、Ｙｉｍの実験式や反応モデルに限定されるものでなく、パラメータの設定や表現方法を変えた式、簡略式や厳密式等、各種の式や各種の反応モデルに基づくことができる。
【数４１】
0005935247.tif 000047
【数４２】
0005935247.tif 000048
【００５７】
ステップＳ３０では、尿素濃度Ｃｕが所定の基準値である濃度閾値ＭＩＮ２以下となったか否かが判定される。濃度閾値ＭＩＮ２は、尿素の分解反応の演算を終了する尿素濃度Ｃｕの下限値である。尿素濃度Ｃｕが濃度閾値ＭＩＮ２以下である場合にはステップＳ３２へ処理を移行させ、そうでない場合には時刻ｔをステップΔｔだけ進めてステップＳ２８の処理を繰り返す。なお、尿素濃度が濃度閾値以下となった否かの判定は、例えば水分量をもっても等価に判定できるものであり、実態として尿素の分解反応の進行度が推定できる他の方法を含めることができる。
【００５８】
ステップＳ３２では、上記の脱硝システム用尿素水分解解析処理によって得られた結果が出力される。処理部１０は、時刻ｔに対する液滴径ｄｐ、温度Ｔｐ、尿素水の濃度Ｗｔ、尿素濃度Ｃｕ、イソシアン酸ＣＨＮＣＯ、水ＣＨ２Ｏ及びアンモニアＣＮＨ３等の変化を出力部１４に出力させる。なお、処理結果の出力は、上記したものに限定されるものでもなく、また、処理途中の結果の出力を含めてもよく、出力の形式も各種の方法を採り得る。
【００５９】
＜脱硝システム用尿素水分解解析処理の評価＞
まず、単一の尿素水液滴の蒸発過程について評価した。ノズルから噴射された尿素水が液体のまま、直接触媒もしくは管壁面に衝突することはできる限り避けたい。そのために、尿素水の蒸発過程における液滴噴射速度、初期液滴径、雰囲気温度の蒸発速度に与える影響を単一液滴で検討した。
【００６０】
単一液滴は実際の噴霧と比べて液滴間の干渉が考慮されないため、液滴の蒸発速度は速くなることが考えられる。しかし実際の噴霧の場合には噴射ノズルの種類や噴射量によって液滴間の干渉度合いが変わってしまうため、ここでは一次的な検討として各条件の影響を抽出しやすい単一液滴にて計算を行った。
【００６１】
図３に雰囲気温度を２００、３００、４００℃、液滴径を５０、１００、２００μｍ、速度を１０、５０、１００ｍ／ｓと変更した時の液滴径の時間変化を示す。雰囲気速度は触媒の空塔速度（ＳＶ）が１１０００ｈ-1で、直径４０ｃｍの円管をガスが通過する時の速度とし、尿素水濃度は４０ｗｔ％、計算期間は液滴温度が１００℃になるまでとしている。液滴噴射方向は雰囲気の流れ方向と等しくした。
【００６２】
図３から、雰囲気温度が低く、液滴径が大きく、噴射速度が遅い方が、液滴寿命が長くなることが確認された。尿素水の蒸発は、液滴径によって雰囲気温度の影響が異なり、小液滴径では温度による蒸発時間の変化を小さくすることができるといえる。それに対して噴射速度の影響は大きくなかった。雰囲気温度はエンジンの種類や運転負荷によって決められてしまうため、排気管内への尿素水噴射が制御できるパラメータは尿素水の噴射速度と液滴径になる。そのため液滴径を小さくすることが、液滴の管内での滞留時間を小さくし、噴射ノズルと触媒の距離を短縮するのに有効であるといえる。また、エンジンの運転条件の影響も小さくできることがわかる。
【００６３】
図４は、図３の条件時の尿素水液滴温度変化を示す。どの条件でも時間が経過して、液滴径が小さくなり、液滴内の尿素濃度が大きくなると液滴温度が急激に上昇した。これは、尿素水の飽和蒸気圧の影響であると考えられる。
【００６４】
飽和蒸気圧は、式（１０）に示されるように液滴内の尿素濃度が上昇すると低下するため、蒸発が進むと蒸発速度は減少していく。そこで尿素濃度の影響を調べるため、図５に雰囲気温度２００、３００℃での尿素濃度を変更した時の液滴寿命の変化を示す。
【００６５】
寿命は液滴から水が蒸発するまでの時間とし、初期液滴径は１００μｍである。尿素濃度の増加によって、液滴寿命は短くなる。これは蒸発すべき液滴内の水分量が減少しているからであると考えられる。寿命に対する濃度の効果はほぼ直線的であるが、濃度が高くなると図５に見られる一定温度で保持される時間が短くなり、液滴温度の上昇速度が上がるため、蒸発速度が低下し寿命の減少は若干緩和される。また、液滴寿命に与える雰囲気温度の影響は尿素濃度の増加によって小さくなる。
【００６６】
雰囲気温度が４００℃の時の尿素水蒸発中の液滴速度の時間変化を図６に示す。初期液滴径によって液滴速度が雰囲気速度と一致する時間は異なっていたが、噴射初期速度が変更されても液滴が雰囲気速度に収束するまでの時間はほとんど変わらなかった。また本条件では、蒸発速度より液滴と雰囲気との運動量交換の速度の方が速く、蒸発時間の約半分は雰囲気に流されながら蒸発していることがわかった。
【００６７】
次に、触媒での分解過程および性能劣化について確認した。尿素水が完全に蒸発して化学式（４）の反応が終了した状態で、全て触媒中で尿素が熱分解、イソシアン酸が加水分解する時の雰囲気温度による脱硝率の変化を調べた。図７にアンモニア（ＮＨ３）ガスを導入した時と比較した結果を示す。空塔速度（ＳＶ）は１１０００ｈ−１とした。
【００６８】
温度が低下するに従って脱硝率は減少したが、還元剤による違いはほとんど現れなかった。このことは、触媒の尿素の熱・加水分解速度が脱硝反応速度よりも速く、脱硝時に還元剤としてアンモニア（ＮＨ３）を用意できていることを示していると考えられる。これより、尿素水を用いる時には、尿素水が触媒に導入される前に液滴を蒸発させ尿素と水を分離しておけば良いことがわかる。また、触媒性能を化学式（６）及び（７）の酸性硫安の生成によって劣化させるために、二酸化硫黄（ＳＯ２）を８００ｐｐｍ含むガスにて８時間運転した時の性能劣化による脱硝率低下への影響も併せて示している。
【化６】
0005935247.tif 000049
【化７】
0005935247.tif 000050
【００６９】
温度が低いほど脱硝率は低くなり、性能劣化する割合も大きくなっていた。本条件では３００℃以下になると性能劣化を生じることがわかった。劣化時の結果から、劣化しない条件より還元剤による違いは現れているが、劣化特性に与える還元剤の影響は小さいと考えられる。
【００７０】
式（３８）〜（４０）に示されるように、尿素の分解過程については一次反応を仮定しており、反応は触媒内での滞留時間に依存している。そこで、空塔速度（ＳＶ）を２５０００ｈ−１としてより反応時間の短い条件での尿素の分解による影響を調べた。二酸化硫黄（ＳＯ２）は導入せず、脱硝性能のみを比較した結果を図８に示す。
【００７１】
触媒に導入される還元剤は、導入前に雰囲気中での尿素の熱分解が完全に終了している条件（ＣＮＨ３０，ＣＨＮＣＯ０＝３００ｐｐｍ）から熱分解が全く行われていない条件（Ｃｕ０＝３００ｐｐｍ）までの熱分解の割合を変えた４条件で行った。比較のため、アンモニア（ＮＨ３）ガス（ＣＮＨ３０＝６００ｐｐｍ）での結果を併せて示す。
【００７２】
アンモニア（ＮＨ３）ガスを導入した場合とそうでない場合において脱硝率がほとんど変わらないことから、空塔速度（ＳＶ）が２５０００ｈ−１になっても触媒には十分な尿素分解能力があることがわかった。温度が低下すると、尿素を分解する能力も低下するため脱硝率に差異が生じているが、その影響は大きくなかった。なお、脱硝率は空塔速度（ＳＶ）が大きくなったことから、図７の１１０００ｈ−１での結果と比べて低下した。
【００７３】
触媒中でのアンモニア（ＮＨ３）の生成過程を調べるため、触媒へ導入される尿素の熱分解割合を図７の条件と同様に変更した時の雰囲気ガス流れ方向の触媒中アンモニア（ＮＨ３）、一酸化窒素（ＮＯ）濃度の変化を図９に示す。計算条件は触媒性能が劣化した二酸化硫黄（ＳＯ２）を８００ｐｐｍ、雰囲気温度を２００℃とした。
【００７４】
触媒入口で尿素の熱分解が完全に終了していない場合には、熱・加水分解過程により、触媒中でアンモニア（ＮＨ３）を生成しながら脱硝を行っているため、触媒入口側のアンモニア（ＮＨ３）濃度は小さくなった。触媒導入前に熱分解が行えておらず、全て尿素として触媒に導入された場合でも触媒を約６割通過して以降、触媒出口までアンモニア（ＮＨ３）濃度はアンモニア（ＮＨ３）ガスの結果と一致した。それに対して一酸化窒素（ＮＯ）濃度はほとんど変化しない。これは各位置で必要なアンモニア（ＮＨ３）が存在し、脱硝反応が行えていることを示していると考えられる。
【産業上の利用可能性】
【００７５】
以上の実施形態は、舶用ディーゼルエンジンのみならず、還元剤として尿素水を用いた他の内燃機関の脱硝触媒システムにも適用することができる。
【符号の説明】
【００７６】
１０  処理部、１２  記憶部、１４  入力部、１６  出力部、１００  コンピュータシステム。
【図１】
0005935247.tif 000051
【図２】
0005935247.tif 000052
【図３】
0005935247.tif 000053
【図４】
0005935247.tif 000054
【図５】
0005935247.tif 000055
【図６】
0005935247.tif 000056
【図７】
0005935247.tif 000057
【図８】
0005935247.tif 000058
【図９】
0005935247.tif 000059