(57)【特許請求の範囲】
【請求項１】
拡散物質を複数の粒子で模擬する濃度分布解析装置であり、
数値解析手法を適用して、入力部から与えられる解析対象とする領域の解析メッシュ、構造物の位置情報、流体の流速と圧力の初期条件、境界条件、並びに、前記解析対象とする領域内での粒子の放出位置および放出速度を含む初期パラメータに基づき、前記解析メッシュに対して少なくとも前記流体の流速分布と乱流エネルギ分布とを計算する流体計算部と、
前記流体計算部が算出した前記流体の流速分布と乱流エネルギ分布とに基づいて、前記解析対象とする領域内での前記粒子の軌跡を追跡計算して、前記粒子の軌跡と前記粒子周りの濃度の広がりを計算する粒子追跡部と、
前記流体計算部が算出した前記流体の流速分布、前記粒子追跡部が算出した前記粒子の位置および前記粒子周りの濃度の広がりに基づいて、前記粒子の各々からの影響を積算して得られる濃度分布に前記流体の流れが変動する影響を考慮した補正処理を施した前記解析メッシュにおける前記拡散物質の濃度分布を計算する濃度分布計算部と、を具備し、
前記濃度分布計算部は、前記流体の流れが変動する影響を考慮した補正係数Ｆkを導入した下記式
【数１】
0005774520.tif 000016
と、前記粒子追跡部が得た各粒子の位置座標情報（ｘｋ，ｙｋ，ｚｋ）とを用いて計算することで、評価対象とする地点（Ｘ，Ｙ，Ｚ）の前記流体の流れが変動する影響を考慮しない濃度χＯ（Ｘ，Ｙ，Ｚ）を計算するとともに、
前記流体の流れが変動する影響として、前記流体の流速の変動の影響を考慮した補正係数Ｆkｖを計算する下記式
【数２】
0005774520.tif 000017
と、前記流体計算部が算出した前記流体の流速分布と前記粒子の位置および評価対象とする地点の情報とに基づいて、各粒子について計算される前記補正係数Ｆkｖ、前記流体の流れが変動する影響として、前記流体の乱流エネルギの変動の影響を考慮した補正係数Ｆkｔを計算する下記式
【数３】
0005774520.tif 000018
と、前記流体計算部が算出した前記流体の流速分布と前記粒子の位置および評価対象とする地点の情報とに基づいて、各粒子について計算される前記補正係数Ｆkｔ、並びに前記補正係数Ｆkｖと前記補正係数Ｆkｔとの積から選択される何れか一つを前記補正係数Ｆkとして決定し、前記濃度χＯ（Ｘ，Ｙ，Ｚ）と前記補正係数Ｆkとから、評価対象とする地点（Ｘ，Ｙ，Ｚ）の前記流体の流れが変動する影響を考慮した濃度χ（Ｘ，Ｙ，Ｚ）を求めるように構成されることを特徴とする濃度分布解析装置。
【請求項２】
前記濃度分布計算部は、前記濃度χＯ（Ｘ，Ｙ，Ｚ）を計算する濃度分布算出処理部と、
前記補正係数Ｆkを計算する補正処理部と、を備えることを特徴とする請求項１記載の濃度分布解析装置。
【請求項３】
前記補正処理部は、前記補正係数Ｆkｔを算出し、得られる記補正係数Ｆkｔを前記補正係数Ｆkとして決定する第１の補正処理部、
前記補正係数Ｆkｖを算出し、得られる前記補正係数Ｆkｖを前記補正係数Ｆkとして決定する第２の補正処理部、
前記補正係数Ｆkｔと前記補正係数Ｆkｖとを算出し、得られる前記補正係数Ｆkｔおよび前記補正係数Ｆkｖの積をさらに算出し、得られる前記補正係数Ｆkｔと前記補正係数Ｆkｖとの積を前記補正係数Ｆkとして決定する第３の補正処理部、および、
前記補正係数Ｆkｔと前記補正係数Ｆkｖとを算出し、得られる前記補正係数Ｆkｔ、得られる前記補正係数Ｆkｖ、および得られる前記補正係数Ｆkｔと前記補正係数Ｆkｖとの積から選択される何れか一つを前記補正係数Ｆkとして決定する第４の補正処理部、から選択される何れか一つであることを特徴とする請求項２記載の濃度分布解析装置。
【請求項４】
前記濃度分布計算部は、前記解析対象とする領域内に存在する前記粒子のうち、濃度評価地点の濃度に影響を及ぼす粒子を選別する濃度影響粒子選別処理部をさらに備えることを特徴とする請求項２または３記載の濃度分布解析装置。
【請求項５】
前記濃度影響粒子選別処理部は、前記流体計算部が算出した前記流体の流速分布を用いて、前記濃度評価地点を通過する流線と各粒子位置との最短距離を計算し、前記濃度評価地点に対して上流にある粒子を前記濃度評価地点の濃度に影響を及ぼす粒子として選別するように構成されることを特徴とする請求項４記載の濃度分布解析装置。
【請求項６】
前記濃度影響粒子選別処理部は、前記流体計算部が算出した前記流体の流速分布と前記構造物の位置情報とを用いて、各粒子位置と前記濃度評価地点とを結ぶ線上に前記構造物が存在するか否かを判断し、前記粒子と前記濃度評価地点とを結ぶ線上に前記構造物が存在しない場合の当該粒子を前記濃度評価地点の濃度に影響を及ぼす粒子として選別するように構成されることを特徴とする請求項４又は５記載の濃度分布解析装置。
【請求項７】
前記粒子追跡部は、前記入力部から与えられる放出物質の密度を用いて、重力沈降の効果を考慮して前記粒子の軌跡を計算するように構成されることを特徴とする請求項１から６の何れか１項に記載の濃度分布解析装置。
【請求項８】
前記粒子追跡部は、前記入力部から与えられる初期パラメータであり、複数放出口の各々から放出される粒子の放出位置および放出速度を用いて、前記複数放出口の各々から放出される各粒子の軌跡を計算するように構成されることを特徴とする請求項１から７の何れか１項に記載の濃度分布解析装置。
【請求項９】
解析メッシュに対して少なくとも流体の流速分布と乱流エネルギ分布とを計算する流体計算部と、前記流体の流速分布と乱流エネルギ分布とに基づいて、解析対象とする領域内での粒子の軌跡と前記粒子周りの濃度の広がりを計算する粒子追跡部と、前記流体の流れが変動する影響を考慮した補正処理を施した前記解析メッシュにおける拡散物質の濃度分布を計算する濃度分布計算部とを具備し、前記濃度分布計算部が、前記流体の流れが変動する影響を考慮した補正係数Ｆkを導入した下記式
【数４】
0005774520.tif 000019
と、前記粒子追跡部が得た各粒子の位置座標情報（ｘｋ，ｙｋ，ｚｋ）とを用いて計算することで、評価対象とする地点（Ｘ，Ｙ，Ｚ）の前記流体の流れが変動する影響を考慮しない濃度χＯ（Ｘ，Ｙ，Ｚ）を計算するとともに、前記流体の流れが変動する影響として、前記流体の流速の変動の影響を考慮した補正係数Ｆkｖを計算する下記式
【数５】
0005774520.tif 000020
と、前記流体計算部が算出した前記流体の流速分布と前記粒子の位置および評価対象とする地点の情報とに基づいて、各粒子について計算される前記補正係数Ｆkｖ、前記流体の流れが変動する影響として、前記流体の乱流エネルギの変動の影響を考慮した補正係数Ｆkｔを計算する下記式
【数６】
0005774520.tif 000021
と、前記流体計算部が算出した前記流体の流速分布と前記粒子の位置および評価対象とする地点の情報とに基づいて、各粒子について計算される前記補正係数Ｆkｔ、並びに前記補正係数Ｆkｖと前記補正係数Ｆkｔとの積から選択される何れか一つを前記補正係数Ｆkとして決定し、前記濃度χＯ（Ｘ，Ｙ，Ｚ）と前記補正係数Ｆkとから、評価対象とする地点（Ｘ，Ｙ，Ｚ）の前記流体の流れが変動する影響を考慮した濃度χ（Ｘ，Ｙ，Ｚ）を求めるように構成される、前記拡散物質を複数の粒子で模擬した濃度分布解析装置に適用される濃度分布解析方法であり、
前記流体計算部が、数値解析手法を適用して、入力部から与えられる解析対象とする領域の解析メッシュ、構造物の位置情報、流体の流速と圧力の初期条件、境界条件、並びに、前記解析対象とする領域内での粒子の放出位置および放出速度を含む初期パラメータに基づき、前記解析メッシュに対して少なくとも前記流体の流速分布と乱流エネルギ分布とを計算するステップと、
前記粒子追跡部が、算出された前記流体の流速分布と乱流エネルギ分布に基づいて、前記解析対象とする領域内での前記粒子の軌跡を追跡計算して、前記粒子の軌跡と前記粒子周りの濃度の広がりを計算するステップと、
前記濃度分布計算部が、算出された前記流体の流速分布、前記粒子の位置および前記粒子周りの濃度の広がりに基づいて、前記粒子の各々からの影響を積算して得られる濃度分布に前記流体の流れが変動する影響を考慮した補正処理を施した前記解析メッシュにおける前記拡散物質の濃度分布を計算するステップと、を具備し、
前記補正処理は、前記濃度χＯ（Ｘ，Ｙ，Ｚ）に前記補正係数Ｆkを乗じる演算処理であることを特徴とする濃度分布解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、拡散物質が大気中で濃度分布を形成する現象を予測する濃度分布解析装置および濃度分布解析方法に関するものである。
【背景技術】
【０００２】
大気中に拡散物質が放出された場合、避難や脅威の除去活動を適切に行うためには、拡散物質の挙動を解析により迅速に予測し当該領域の濃度分布を出力することが重要となる。特に市街地の場合は、建物の影響により気流は複雑な流れ場を形成していることから、３次元の詳細な気流情報を使用した高精度かつ迅速な濃度評価が要求される。
【０００３】
ある地点における濃度分布を得る方法としては、データベース、実測値、または、数値流体力学計算から得た気流情報に基づいて、拡散物質に対しスカラー輸送方程式を解き、濃度情報を出力する方法がある。このような気流情報に基づいて、拡散物質に対しスカラー輸送方程式を解いて濃度情報を出力する方法は、例えば、下記非特許文献１に記載される。また、他の方法の一例としては、粒子追跡法を使用した濃度評価手法が挙げられる。
【０００４】
例えば、図７に示されるように、粒子追跡法による濃度評価手法では、解析メッシュ１を設定し、放出源３から拡散する拡散物質の挙動を複数個の離散的な粒子２の挙動で代表させて追跡し、得られた粒子分布の情報を使用する。このとき、粒子２の周りには空間的に統計的な濃度の広がりを持たせる。そして、濃度を評価したい地点に対して、各粒子２からの影響を積算することによって当該地点における濃度を算出する。このような、粒子追跡法を使用した濃度評価手法による数値解析方法および数値解析装置は、例えば、特許第３８８１９２６号公報（特許文献１）に記載される。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特許第３８８１９２６号公報
【非特許文献】
【０００６】
【非特許文献１】半場藤弘著、「大気境界層のスカラーの乱流モデル」、生産研究、東京大学生産技術研究所、１９９４年２月、第４６巻、第２号、ｐ.７４−７９
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、非特許文献１に記載される技術を適用する場合、精度の良い解を出力するためには、評価地点周りの広い空間に対して解析メッシュを作成する必要がある。すなわち、多数の解析メッシュを作成する必要がある。解析メッシュの数が多くなるほど計算時間が膨大となるため、非特許文献１に記載される技術を適用する場合には、所望の精度を維持しつつ迅速に結果を出力したいという要請に応えることが困難という課題がある。
【０００８】
一方、特許文献１に記載される技術を適用すれば、粒子追跡法を使用するため、計算時間に関しては解析メッシュに依存せず迅速に結果を出力することができるという点で迅速に結果を得るという要請に応えることはできるものの、流れが大きく変動し得る環境下（例えば、市街地等）で適用する場合、流れが大きく変動する場合の影響を考慮していないため、所望の精度を得ることができない場合が生じ得る。
【０００９】
例えば、図８に示されるように、大気の流速および方向を示す流速ベクトル４を持つ解析対象領域内に、建物等の構造物５が存在する場合、構造物５の上方に位置する領域Ｒ１と領域Ｒ１の下方に位置する領域、すなわち、構造物５の背面側（図８において右側）とでは、流速が大きく異なる。その結果、構造物５が流体（大気）を遮らない領域Ｒ１では、流速が速く対流の効果が大きくなるので、構造物５の背面側の領域に存在する粒子が濃度評価地点Ｐ１へ及ぼす影響は小さくなる。一方、構造物５が流体を遮る構造物５の背面側の領域では、領域Ｒ１と比較して対流の効果は小さくなるので、構造物５の背面側の領域にある粒子が構造物５の背面側の領域の濃度へ及ぼす影響は濃度評価値点Ｐ１よりも相対的に大きくなる。
【００１０】
このような流れが大きく変動する場合の影響を考慮しない従来技術では、得られた拡散現象および濃度分布が、当初要求していた精度に達しない場合が生じており、迅速に結果が得られたとしても所望の精度を維持できない場合が生じ得る。すなわち、所望の精度を維持しつつ迅速に結果を出力したいという要請に応えることができていないという課題がある。
【００１１】
本発明は、上述した事情に鑑みてなされたものであり、市街地等の建物の影響により流れが大きく変動する環境下における拡散現象を解析する際に、粒子追跡手法を用いて、迅速かつ流れの変動の影響を考慮した高精度の濃度分布を出力可能な濃度分布解析装置および濃度分布解析方法を提供することを目的とする。
【課題を解決するための手段】
【００１２】
本発明の実施形態に係る濃度分布解析装置は、上述した課題を解決するため、拡散物質を複数の粒子で模擬する濃度分布解析装置であり、数値解析手法を適用して、入力部から与えられる解析対象とする領域の解析メッシュ、構造物の位置情報、流体の流速と圧力の初期条件、境界条件、並びに、前記解析対象とする領域内での粒子の放出位置および放出速度を含む初期パラメータに基づき、前記解析メッシュに対して少なくとも前記流体の流速分布と乱流エネルギ分布とを計算する流体計算部と、前記流体計算部が算出した前記流体の流速分布と乱流エネルギ分布とに基づいて、前記解析対象とする領域内での前記粒子の軌跡を追跡計算して、前記粒子の軌跡と前記粒子周りの濃度の広がりを計算する粒子追跡部と、前記流体計算部が算出した前記流体の流速分布、前記粒子追跡部が算出した前記粒子の位置および前記粒子周りの濃度の広がりに基づいて、前記粒子の各々からの影響を積算して得られる濃度分布に前記流体の流れが変動する影響を考慮した補正処理を施した前記解析メッシュにおける前記拡散物質の濃度分布を計算する濃度分布計算部と、を具備し、前記濃度分布計算部は、前記流体の流れが変動する影響を考慮した補正係数Ｆkを導入した下記式
【数１】
0005774520.tif 000002
と、前記粒子追跡部が得た各粒子の位置座標情報（ｘｋ，ｙｋ，ｚｋ）とを用いて計算することで、評価対象とする地点（Ｘ，Ｙ，Ｚ）の前記流体の流れが変動する影響を考慮しない濃度χＯ（Ｘ，Ｙ，Ｚ）を計算するとともに、前記流体の流れが変動する影響として、前記流体の流速の変動の影響を考慮した補正係数Ｆkｖを計算する下記式
【数２】
0005774520.tif 000003
と、前記流体計算部が算出した前記流体の流速分布と前記粒子の位置および評価対象とする地点の情報とに基づいて、各粒子について計算される前記補正係数Ｆkｖ、前記流体の流れが変動する影響として、前記流体の乱流エネルギの変動の影響を考慮した補正係数Ｆkｔを計算する下記式
【数３】
0005774520.tif 000004
と、前記流体計算部が算出した前記流体の流速分布と前記粒子の位置および評価対象とする地点の情報とに基づいて、各粒子について計算される前記補正係数Ｆkｔ、並びに前記補正係数Ｆkｖと前記補正係数Ｆkｔとの積から選択される何れか一つを前記補正係数Ｆkとして決定し、前記濃度χＯ（Ｘ，Ｙ，Ｚ）と前記補正係数Ｆkとから、評価対象とする地点（Ｘ，Ｙ，Ｚ）の前記流体の流れが変動する影響を考慮した濃度χ（Ｘ，Ｙ，Ｚ）を求めるように構成されることを特徴とする。
【００１３】
本発明の実施形態に係る濃度分布解析方法は、上述した課題を解決するため、解析メッシュに対して少なくとも流体の流速分布と乱流エネルギ分布とを計算する流体計算部と、前記流体の流速分布と乱流エネルギ分布とに基づいて、解析対象とする領域内での粒子の軌跡と前記粒子周りの濃度の広がりを計算する粒子追跡部と、前記流体の流れが変動する影響を考慮した補正処理を施した前記解析メッシュにおける拡散物質の濃度分布を計算する濃度分布計算部とを具備し、前記濃度分布計算部が、前記流体の流れが変動する影響を考慮した補正係数Ｆkを導入した下記式
【数４】
0005774520.tif 000005
と、前記粒子追跡部が得た各粒子の位置座標情報（ｘｋ，ｙｋ，ｚｋ）とを用いて計算することで、評価対象とする地点（Ｘ，Ｙ，Ｚ）の前記流体の流れが変動する影響を考慮しない濃度χＯ（Ｘ，Ｙ，Ｚ）を計算するとともに、前記流体の流れが変動する影響として、前記流体の流速の変動の影響を考慮した補正係数Ｆkｖを計算する下記式
【数５】
0005774520.tif 000006
と、前記流体計算部が算出した前記流体の流速分布と前記粒子の位置および評価対象とする地点の情報とに基づいて、各粒子について計算される前記補正係数Ｆkｖ、前記流体の流れが変動する影響として、前記流体の乱流エネルギの変動の影響を考慮した補正係数Ｆkｔを計算する下記式
【数６】
0005774520.tif 000007
と、前記流体計算部が算出した前記流体の流速分布と前記粒子の位置および評価対象とする地点の情報とに基づいて、各粒子について計算される前記補正係数Ｆkｔ、並びに前記補正係数Ｆkｖと前記補正係数Ｆkｔとの積から選択される何れか一つを前記補正係数Ｆkとして決定し、前記濃度χＯ（Ｘ，Ｙ，Ｚ）と前記補正係数Ｆkとから、評価対象とする地点（Ｘ，Ｙ，Ｚ）の前記流体の流れが変動する影響を考慮した濃度χ（Ｘ，Ｙ，Ｚ）を求めるように構成される、前記拡散物質を複数の粒子で模擬した濃度分布解析装置に適用される濃度分布解析方法であり、前記流体計算部が、数値解析手法を適用して、入力部から与えられる解析対象とする領域の解析メッシュ、構造物の位置情報、流体の流速と圧力の初期条件、境界条件、並びに、前記解析対象とする領域内での粒子の放出位置および放出速度を含む初期パラメータに基づき、前記解析メッシュに対して少なくとも前記流体の流速分布と乱流エネルギ分布とを計算するステップと、前記粒子追跡部が、算出された前記流体の流速分布と乱流エネルギ分布に基づいて、前記解析対象とする領域内での前記粒子の軌跡を追跡計算して、前記粒子の軌跡と前記粒子周りの濃度の広がりを計算するステップと、前記濃度分布計算部が、算出された前記流体の流速分布、前記粒子の位置および前記粒子周りの濃度の広がりに基づいて、前記粒子の各々からの影響を積算して得られる濃度分布に前記流体の流れが変動する影響を考慮した補正処理を施した前記解析メッシュにおける前記拡散物質の濃度分布を計算するステップと、を具備し、前記補正処理は、前記濃度χＯ（Ｘ，Ｙ，Ｚ）に前記補正係数Ｆkを乗じる演算処理であることを特徴とする。
【発明の効果】
【００１４】
本発明によれば、市街地等の流れが大きく変動する環境下における拡散現象を解析する際に粒子追跡法を用いて、従来技術では困難であった流れの変動の影響が濃度分布に与える影響を高速に評価でき、迅速かつ高精度な濃度分布情報をユーザに提示することができる。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施形態に係る濃度分布解析装置の構成図。
【図２】本発明の実施形態に係る濃度分布解析装置で導入する補正係数Ｆkを算出する一例を説明する説明図。
【図３】解析対象領域内に存在する粒子が濃度評価地点に及ぼす影響の一例（第１の例）を説明する説明図。
【図４】解析対象領域内に存在する粒子が濃度評価地点に及ぼす影響の一例（第２の例）を説明する説明図。
【図５】本発明の実施形態に係る濃度分布解析装置の構成図。
【図６】本発明の実施形態に係る濃度分布解析装置の構成図。
【図７】粒子追跡法を用いた粒子追跡計算の概要を説明する説明図。
【図８】粒子追跡法を用いた粒子追跡計算における流れ変動の影響を説明する説明図。
【発明を実施するための形態】
【００１６】
以下、本発明の実施形態に係る濃度分布解析装置および濃度分布解析方法について、添付図面を参照して説明する。
【００１７】
なお、以下の説明は、本発明の実施形態に係る濃度分布解析装置および濃度分布解析方法の適用の一例として、解析対象となる空間が市街地であり、解析対象となる空間内に設定される解析対象領域の一点から拡散物質が空間に放出され、濃度分布を形成する現象を解析により予測する場合の例である。
【００１８】
図１は、本発明の実施形態に係る濃度分布解析装置の一例である濃度分布解析装置１０の構成を示す概略図である。
【００１９】
濃度分布解析装置１０は、例えば、濃度分布解析装置１０の各手段である入力部１１、流体計算部１２、粒子追跡部１３、濃度分布計算部１４、計算終了判定部１５および表示部１６として機能させる数値解析プログラム（図１において省略）をコンピュータに実行させることによって、ハードウェア資源であるコンピュータとソフトウェア資源であるプログラムとが協働して実現される装置である。
【００２０】
入力部１１は、例えば、コンピュータとインターフェイスを介して接続される入力装置またはコンピュータ自身が備えるキーボードやマウス等の入力手段によって実現される。入力部１１は、入力された情報の種類に応じて、流体計算部１２、粒子追跡部１３、濃度分布計算部１４、計算終了判定部１５および表示部１６の何れかに入力された情報を与える。
【００２１】
入力部１１から与える初期パラメータとしては、例えば、解析対象領域の解析メッシュ、建物等の構造物の位置情報、境界条件、流体の流速と圧力の初期条件、流れ場の密度、粘性係数の物性値、拡散物質が放出される放出口の位置並びに放出物質（粒子）の放出速度および密度等である。
【００２２】
流体計算部１２は、数値流体力学を用いて、対象空間領域の三次元流速分布を計算する機能と、乱流エネルギ分布を計算する機能とを有し、解析メッシュに対し、初期パラメータに基づいて、有限体積法や有限差分法などの数値解析手法を適用する。
【００２３】
流体計算部１２は、流体の質量、運動量の２つの保存式と、例えば、ｋ−ε乱流モデルを用いて乱流エネルギとエネルギ散逸率の輸送方程式を解いて、１タイムステップ後の流体の流速、圧力、乱流エネルギ、エネルギ散逸率を各メッシュに対して計算し、これらの物理量の空間的な分布を導く。流体計算部１２は、得られた対象空間領域における流体の流速分布情報と乱流エネルギ分布情報とを粒子追跡部１３に与える。
【００２４】
粒子追跡部１３は、解析対象領域内における各粒子の位置情報と各粒子周りの濃度の分布情報とを計算する際に必要な数式情報を有しており、拡散物質を複数の粒子で模擬し、流体計算部１２で得られた流速分布情報と乱流エネルギ分布情報とに基づいて、粒子の軌跡を追跡計算して、解析対象領域内における各粒子の位置情報と各粒子周りの濃度の分布情報とを計算する機能を有する。例えば、ラグランジュ粒子追跡法を用いて計算する場合、模擬する粒子の中心位置の座標ｘｉ（ｉ＝１〜３）は、次の式（１）から求めることができる。
【００２５】
【数７】
0005774520.tif 000008
【００２６】
粒子移動速度Ｕｉについては、乱流変動速度成分ｕｉ’によるランダムな挙動を考慮し、次の式（２）を用いて算出する。
【００２７】
【数８】
0005774520.tif 000009
【００２８】
また、粒子移動速度Ｕｉを算出する際に重力沈降を考慮する場合には、上記式（２）の右辺に重力沈降を考慮した項（−δi3Ｖ）を加算した次の式（３）を用いて算出する。なお、重力沈降速度Ｖについては、Ｓｔｏｋｅｓ則を用いて、次の式（４）より求めることができる。
【００２９】
【数９】
0005774520.tif 000010
【００３０】
【数１０】
0005774520.tif 000011
【００３１】
粒子追跡部１３は、上述した式（１）〜（８）を用いて、各粒子の位置情報および各粒子周りの濃度の分布情報を計算した後、計算の結果得られた各粒子の位置情報および各粒子周りの濃度の分布情報を濃度分布計算部１４に与える。なお、粒子追跡部１３は、解析対象領域内に複数の粒子放出口を有する場合、すなわち、解析対象領域の複数点から拡散物質が空間に放出される場合においても、各粒子放出口から放出されるそれぞれ粒子に対して、粒子の位置情報および粒子周りの濃度の分布情報を計算することができる。
【００３２】
濃度分布計算部１４は、評価対象領域の濃度分布を計算するのに必要な数式情報を有しており、当該数式情報、流体計算部１２によって得られた評価対象領域の流速分布情報並びに粒子追跡部１３によって得られた各粒子の位置情報および各粒子周りの濃度の広がりの情報に基づいて、各粒子からの影響を積算すると同時に、評価対象空間の流速分布情報を用いて、例えば、流速変動や乱流エネルギ変動等の流れの変動の影響までを考慮した評価対象領域の濃度分布情報を計算する。
【００３３】
より詳細に説明すれば、濃度分布計算部１４は、各々の粒子に対して流れの変動の影響までを考慮した補正係数Ｆkを導入した次の式（９ｂ）と、粒子追跡部１３で得られた各粒子の位置座標情報（ｘｋ，ｙｋ，ｚｋ）とを用いて計算することで、評価対象とする地点（Ｘ，Ｙ，Ｚ）の流れの変動の影響までを考慮した濃度χ（Ｘ，Ｙ，Ｚ）[ｇ／ｍ３]を求める。ここで、χＯ（Ｘ，Ｙ，Ｚ）[ｇ／ｍ３]は、評価対象とする地点（Ｘ，Ｙ，Ｚ）の流れの変動の影響までを考慮しない濃度、すなわち、補正係数Ｆk導入前における評価対象とする地点（Ｘ，Ｙ，Ｚ）の濃度である。
【００３４】
【数１１】
0005774520.tif 000012
【００３５】
すなわち、濃度分布計算部１４は、評価対象領域の流速分布情報並びに各粒子の位置情報および各粒子周りの濃度の広がりの情報に基づいて、各粒子からの影響を積算して得られる濃度分布であるχＯ（Ｘ，Ｙ，Ｚ）[ｇ／ｍ３]に、構造物の存在による流れの変動の影響を考慮して決定される補正係数Ｆkを乗じることで、流れの変動の影響までを考慮した評価対象とする地点（Ｘ，Ｙ，Ｚ）の濃度χ（Ｘ，Ｙ，Ｚ）[ｇ／ｍ３]を求める。
【００３６】
濃度分布計算部１４は、流れの変動の影響を考慮する前の評価地点（Ｘ，Ｙ，Ｚ）における濃度分布χＯ（Ｘ，Ｙ，Ｚ）[ｇ／ｍ３]を求める濃度分布算出処理部２１と、濃度分布算出処理部２１が計算した評価対象領域の濃度分布に対して、構造物の存在による流れの変動の影響を考慮した補正処理、すなわち、補正係数Ｆkを決定し、決定した補正係数Ｆkを乗じた評価対象となる地点（Ｘ，Ｙ，Ｚ）の濃度χ（Ｘ，Ｙ，Ｚ）[ｇ／ｍ３]を求める補正処理部２３とを備える。濃度分布計算部１４は、算出結果である評価対象となる地点（Ｘ，Ｙ，Ｚ）の濃度χ（Ｘ，Ｙ，Ｚ）[ｇ／ｍ３]を計算終了判定部１５に与える。
【００３７】
また、図１に示される濃度分布計算部１４は、後述する図３および図４に示されるように、解析対象領域内に存在する全ての粒子２が濃度評価地点の濃度に影響を及ぼす訳ではないため、計算処理時間をより短縮化（合理化）する観点から、解析対象領域内に存在する粒子２のうち、濃度評価地点の濃度に影響を及ぼす粒子を選別する濃度影響粒子選別処理部２４を備える。
【００３８】
濃度影響粒子選別処理部２４は、濃度評価地点の濃度に影響を及ぼす粒子を選別する処理ステップを実行するモード（以下、「粒子選別モード」と称する。）を入（実行）と切（実行せず）とを切替選択する機能を有し、ユーザが必要に応じて入力部１１から指令を与えることで、粒子選別モードの入と切とを自在に切り替えることができる。
【００３９】
なお、濃度分布解析装置１０において、濃度影響粒子選別処理部２４は、任意の構成要素である。すなわち、濃度影響粒子選別処理部２４を備えない濃度分布計算部１４を構成することもできる。
【００４０】
計算終了判定部１５は、濃度分布計算部１４の算出結果および予め指定された計算終了条件に基づいて、予め指定された計算終了条件を満足するか否かを判定する機能を有する。計算終了条件の与え方としては、例えば、流体の質量、運動量およびエネルギの三つの保存法則から与えられる収束条件を与えたり、パラメータおよび値を指定して与えたり、計算を終了するタイムステップ数で与えたり等、幾つかの方法を採用することができる。
【００４１】
計算終了判定部１５は、予め指定された計算終了条件を満足していると判断した場合、現在のタイムステップで算出した解析対象領域における濃度分布を計算結果として表示部１６へ与える一方、予め指定された計算終了条件を満足していないと判断した場合、現在のタイムステップからさらにタイムステップを進めて、対象空間領域の三次元流速分布と乱流エネルギ分布の計算以降の計算を、計算終了条件を満足するまで繰り返し実行する。
【００４２】
表示部１６は、例えば、コンピュータとインターフェイスを介して接続される表示装置またはコンピュータ自身が備えるディスプレイ等の表示手段によって実現される。表示部１６は、入力部１１が受け付けた表示要求または計算終了判定部１５から計算結果の表示要求を受け取ると、当該表示要求に応じた内容を画面表示する。
【００４３】
このように構成される濃度分布解析装置１０では、入力部１１から解析対象とする市街地等の空間を分割した解析メッシュおよび境界条件と初期パラメータを与えて、流体計算部１２が解析メッシュに対して境界条件および初期パラメータに基づいて少なくとも流速と乱流エネルギを計算する。
【００４４】
続いて、流体計算部１２によって計算された流速および乱流エネルギに基づいて、粒子追跡部１３が粒子の軌跡および粒子周りの濃度の広がりを計算する。さらに続いて、粒子追跡部１３によって得られた各粒子の位置情報および各粒子周りの濃度の広がりの情報並びに流体計算部１２によって得られた評価対象領域の流速分布情報に基づいて、濃度分布計算部１４が構造物の存在による流れの変動の影響までを考慮した評価対象領域の濃度分布情報を計算する。
【００４５】
そして、計算終了判定部１５は、予め設定された計算終了条件を満たしているかを判別し、計算終了条件を満たしている場合には、濃度分布情報を表示部１６に与えられ、濃度分布が表示部１６に表示される。
【００４６】
続いて、濃度分布解析装置１０で導入する流れの変動の影響を考慮した補正処理（補正係数Ｆk）について説明する。
【００４７】
図２は、濃度分布解析装置１０で導入する補正係数Ｆkを算出する一例を説明する説明図である。
【００４８】
濃度分布解析装置１０で導入する補正係数Ｆkは、例えば、図２に示されるように、流速の変動の影響を考慮した補正係数Ｆkｖであり、粒子位置および評価地点の流速情報を用いて、各粒子に対して、次の式（１１）により求められる。
【００４９】
【数１２】
0005774520.tif 000013
【００５０】
また、補正係数Ｆk（＝Ｆkｔ）に関しては、例えば、次の式（１３）に示されるように、乱流エネルギの変動を考慮して与えることもできる。
【００５１】
【数１３】
0005774520.tif 000014
【００５２】
上記式（１３）の右辺の括弧内第２項分子では、粒子位置と濃度評価地点間の乱流エネルギの変化を評価し、いずれか大きい方の乱流エネルギにて除算することにより規格化している。補正係数Ｆk（＝Ｆkｔ）は、２点の乱流エネルギの差が大きいときに小さく、差が小さいときに大きく評価される。従って、濃度分布解析装置１０では、流れ場の乱流エネルギの変動が濃度分布に与える影響を評価することが可能となる。
【００５３】
なお、以下の式（１４）に示されるように、流速の変動および乱流エネルギの変動の両方を考慮した補正係数Ｆkを与えることもできる。
【００５４】
【数１４】
0005774520.tif 000015
【００５５】
続いて、濃度影響粒子選別処理部２４（濃度分布計算部１４）が行う濃度影響粒子選別処理について説明する。
【００５６】
図３は、解析対象領域内に存在する粒子２（２ａ，２ｂ）が濃度評価地点Ｐ１，Ｐ２に及ぼす影響の一例（第１の例）を説明する説明図であり、図４は、解析対象領域内に存在する粒子２（２ｃ，２ｄ）が濃度評価地点Ｐ１，Ｐ２に及ぼす影響の一例（第２の例）を説明する説明図である。
【００５７】
濃度分布計算部１４の濃度影響粒子選別処理部２４は、例えば、図３に示されるように、流体計算部１２が計算した流速分布を用いて、濃度評価地点Ｐ１を通過する流線２９を計算し、さらに、算出した流線２９と粒子位置の最短距離を計算する。計算の結果、濃度評価地点に対して下流にある粒子２については、その粒子２は濃度評価地点の濃度に影響を及ぼさない（または無視できる程に小さい）とみなして、当該粒子２を濃度計算の際に無視（除外）する。すなわち、濃度評価地点の濃度に影響を及ぼす粒子２として、濃度評価地点に対して上流にある粒子２を選別する。
【００５８】
図３に示される例で説明すれば、濃度評価地点Ｐ１に対して下流にある粒子２ａについては濃度評価地点Ｐ１の濃度に影響を及ぼさない（または無視できる程に小さい）とみなして、計算に使用する粒子２から除外される一方、濃度評価地点Ｐ１に対して上流にある粒子２ｂについては除外されることなく濃度を計算する際に使用される。
【００５９】
また、濃度影響粒子選別処理部２４は、例えば、図４に示されるように、流体計算部１２が求めた流速分布と構造物５の位置情報とを用いて、粒子位置と濃度を計算するメッシュ位置（濃度評価地点）とを結ぶ線上に構造物５等の障害物の有無を判断し、障害物が存在すると判別した場合、メッシュ位置に対して、当該粒子２は影響を及ぼさない（または無視できる程に小さい）とみなす。すなわち、濃度影響粒子選別処理部２４は、粒子位置と濃度評価地点とを結ぶ線上に障害物が存在すると判別した場合、当該粒子２を濃度評価地点における濃度に影響を及ぼすものとして選別する。
【００６０】
図４に示される例で説明すれば、濃度評価地点Ｐ１については、粒子２と濃度評価地点Ｐ１との間に構造物５が存在しないため、濃度評価地点Ｐ１における濃度の算出の際には、各粒子２が濃度評価地点Ｐ１の濃度に影響を与えるものとして計算される。一方、濃度評価地点Ｐ２については、一部の粒子２と濃度評価地点Ｐ２との間に構造物５が障害物として存在する。
【００６１】
より詳細に説明すると、粒子２ｃは構造物５が障害物とならないものの、粒子２ｄは構造物５が障害物となっている。このような位置関係にある粒子２ｃ，２ｄに対して、濃度影響粒子選別処理部２４は、濃度評価地点との関係で構造物５が障害物とならない粒子２ｃについては濃度計算に使用する粒子２として選別する一方、濃度評価地点との関係で構造物５が障害物となっている粒子２ｄについては、濃度計算に使用する粒子２から除外する。
【００６２】
このように構成される濃度分布解析装置１０によれば、市街地等の流れが大きく変動する環境下における拡散現象を解析する際に粒子追跡法を用いて、従来技術では困難であった流れの変動の影響が濃度分布に与える影響を高速に評価でき、迅速かつ高精度な濃度分布情報をユーザに提示することができる。すなわち、流れ場の流速の変動や乱流エネルギの変動の影響を考慮した濃度分布を迅速にユーザに提供することができる。
【００６３】
また、濃度分布解析装置１０によれば、粒子２が濃度評価地点における濃度に影響を及ぼすか及ぼさないとみなすことができるかを判断し、影響を及ぼさないとみなすことができると判断した場合には当該粒子２を濃度計算する際の対象から除外するので、ユーザが望む濃度計算の精度を満足しつつ計算時間を短縮することができる。
【００６４】
さらに、濃度分布解析装置１０によれば、上記式（３）および（４）を適用することで、空気等の流体と粒子２の密度差に因る重力の影響を考慮した濃度分布を提供することができる。
【００６５】
なお、上述した濃度分布解析装置１０は、解析対象領域の一点から拡散物質が空間に放出され、濃度分布を形成する現象を解析により予測する場合の一例を説明したものであるが、複数点から拡散物質が空間に放出される場合についても適用できる。解析対象領域の複数点から拡散物質が空間に放出される場合については、各点の予測結果を算出し、算出結果を積算する（重ね合わせる）ことで、解析対象領域の一点から拡散物質が空間に放出された場合と同様に濃度分布を求めることができる。
【００６６】
次に、本発明の実施形態に係る濃度分布解析方法の一例として、濃度分布解析装置１０が行う濃度分布解析手順について説明する。
【００６７】
図５は、本発明の実施形態に係る濃度分布解析装置の一例である濃度分布解析装置１０が行う濃度分布解析手順の処理ステップを示す処理フロー図である。
【００６８】
図５に示されるように、濃度分布解析装置１０が行う濃度分布解析手順（ステップＳ１〜ステップＳ５）は、流体計算部１２が解析対象となる空間領域の３次元流速分布および乱流エネルギ分布を算出する流体算出工程（ステップＳ１）と、粒子追跡部１３が各粒子の移動および各粒子周囲の濃度を算出する粒子追跡工程（ステップＳ２）と、濃度分布計算部１４が流れの変動の影響までを考慮した評価対象領域の濃度分布を計算する濃度分布算出工程（ステップＳ３）と、計算終了判定部１５が計算終了条件を満足するか否かを判定する計算終了条件成立判定工程（ステップＳ４）と、を具備する。
【００６９】
濃度分布解析手順は、解析対象領域の解析メッシュ、建物等の構造物の位置情報、境界条件、流速と圧力の初期条件、流れ場の密度、粘性係数の物性値、拡散物質が放出される放出口の位置および粒子の放出速度等の初期パラメータおよび濃度分布解析手順の実行要求が濃度分布解析装置１０に与えられると開始される（ＳＴＡＲＴ）。
【００７０】
濃度分布解析手順が開始されると、流体計算部１２が、解析対象となる空間領域の３次元流速分布および乱流エネルギ分布を算出する（ステップＳ１）。続いて、粒子追跡部１３が、ステップＳ１で得られた流速分布情報と乱流エネルギ分布情報とに基づいて、各粒子の移動および各粒子周囲の濃度を算出する（ステップＳ２）。続いて、濃度分布計算部１４が流れの変動の影響までを考慮した前記式（９b）と、粒子追跡部１３で得られた各粒子の位置座標情報とを用いて計算し、評価対象とする地点（Ｘ，Ｙ，Ｚ）の流れの変動の影響を考慮した濃度χ（Ｘ，Ｙ，Ｚ）[ｇ／ｍ３]を求めることで、評価対象領域の濃度分布を計算する（ステップＳ３）。
【００７１】
そして、ステップＳ３の終了後、計算終了判定部１５が計算終了条件を満たすと判断した場合（ステップＳ４でＹＥＳの場合）には、濃度分布解析手順の全処理ステップを終了する（ＥＮＤ）。一方、計算終了判定部１５が計算終了条件を満たさないと判断した場合（ステップＳ４でＮＯの場合）には、ステップＳ１に戻り、次のタイムステップについて、ステップＳ１〜ステップＳ４の処理ステップを実行する。
【００７２】
図６は、濃度分布解析手順が備える濃度分布算出工程のより詳細な処理ステップを示す処理フロー図である。
【００７３】
濃度分布算出工程（ステップＳ３１〜ステップＳ３８）は、濃度影響粒子選別処理部２４が解析対象とする領域内に存在する粒子のうち、濃度評価地点の濃度に影響を及ぼす粒子を選別するステップ（ステップＳ３１，ステップＳ３２）と、濃度分布算出処理部２１が流体の流れの変動の影響を考慮する前の評価対象領域の濃度分布を算出するステップ（ステップＳ３３）と、補正処理部２３が補正係数Ｆkを決定し、決定した補正係数ＦkをステップＳ３３で算出された流体の流れの変動の影響を考慮する前の評価対象領域の濃度分布に乗じて流体の流れの変動の影響を考慮した補正処理を施すステップ（ステップＳ３４〜ステップＳ３８）と、を備える。
【００７４】
濃度分布算出工程では、まず、濃度影響粒子選別処理部２４が粒子選別モードの入切状態を判断する。粒子選別モードが入（ＯＮ）の場合（ステップＳ３１でＹＥＳの場合）、ステップＳ３２に進み、濃度影響粒子選別処理部２４が解析対象とする領域内に存在する粒子のうち、濃度評価地点の濃度に影響を及ぼさないと判断できる粒子については計算を省略（スキップ）する。すなわち、濃度影響粒子選別処理部２４は濃度評価地点の濃度に影響を及ぼす粒子を選別する。
【００７５】
続くステップＳ３３では、濃度分布算出処理部２１が流体の流れの変動の影響を考慮する前の評価対象領域の濃度分布を算出する。すなわち、評価地点（Ｘ，Ｙ，Ｚ）における濃度分布χＯ（Ｘ，Ｙ，Ｚ）[ｇ／ｍ３]が算出される。
【００７６】
続くステップＳ３４では、補正処理部２３が流速の変動の影響を考慮する要求が与えられているか否かを判断し、当該要求が与えられている場合（ステップＳ３４でＹＥＳの場合）、ステップＳ３５に進み、補正処理部２３は前記式（１１）および（１２）を用いて、流速の変動の影響を考慮した補正係数Ｆkｖを計算する。
【００７７】
続いて、ステップＳ３６では、補正処理部２３が乱流エネルギの変動の影響を考慮する要求が与えられているか否かを判断し、当該要求が与えられている場合（ステップＳ３６でＹＥＳの場合）、ステップＳ３７に進み、補正処理部２３は前記式（１３）を用いて、乱流エネルギの変動の影響を考慮した補正係数Ｆkｔを計算する。
【００７８】
続くステップＳ３８では、補正処理部２３が流速の変動の影響のみを考慮するのか、乱流エネルギの変動の影響のみを考慮するのか、並びに、流速の変動の影響および乱流エネルギの変動の影響を考慮するのかに応じて、最終的な補正係数Ｆkを決定する。すなわち、Ｆk＝Ｆkｖ、Ｆk＝ＦkｔまたはＦk＝Ｆkｖ・Ｆkｔ（ＦkｖとＦkｔとの積：前記式（１４）に相当）の何れかを補正係数Ｆkとして決定し、ステップＳ３３で得られた流体の流れの変動の影響を考慮する前の評価対象領域の濃度分布に決定した補正係数Ｆkを乗じる。この結果、補正処理部２３は、流体の流れの変動の影響を考慮した評価対象領域の濃度分布を得る。
【００７９】
一方、濃度分布算出工程において、粒子選別モードが切（ＯＦＦ）の場合（ステップＳ３１でＮＯの場合）には、ステップＳ３３に進み、ステップＳ３３以降の処理ステップが実行される。
【００８０】
また、補正処理部２３に流速の変動の影響を考慮する要求が与えられていない場合（ステップＳ３４でＮＯの場合）、ステップＳ３７に進み、ステップＳ３７以降の処理ステップが実行される。すなわち、濃度分布算出工程では、ステップＳ３５およびステップＳ３７の何れか一方が実行され、流体の流れの変動の影響を考慮した補正処理が実行される。
【００８１】
さらに、補正処理部２３に乱流エネルギの変動の影響を考慮する要求が与えられていない場合（ステップＳ３６でＮＯの場合）、ステップＳ３８に進み、ステップＳ３８以降の処理ステップが実行される。
【００８２】
なお、図６に示される濃度分布算出工程は、濃度分布算出工程の一例であり、必ずしも図６に示される処理ステップに限定されない。例えば、ステップＳ３１とステップＳ３２は、濃度分布算出工程の任意の工程であり、濃度影響粒子選別処理部２４を備えない濃度分布計算部１４が濃度分布算出工程を実行する際には、ステップＳ３１とステップＳ３２の処理ステップは省略される。また、ステップＳ３４およびステップＳ３５とステップＳ３６およびステップＳ３７とは、それぞれ、入れ替わっていても差し支えない。
【００８３】
このような濃度分布解析手順によれば、市街地等の流れが大きく変動する環境下における拡散現象を解析する際に粒子追跡法を用いて、従来技術では困難であった流れの変動の影響が濃度分布に与える影響を高速に評価でき、迅速かつ高精度な濃度分布情報をユーザに提示することができる。すなわち、流れ場の流速の変動や乱流エネルギの変動の影響を考慮した濃度分布を迅速にユーザに提供することができる。
【００８４】
また、粒子２が濃度評価地点における濃度に影響を及ぼすか及ぼさないとみなすことができるかを判断し、影響を及ぼさないとみなすことができると判断した場合には当該粒子２を濃度計算する際の対象から除外するので、ユーザが望む濃度計算の精度を満足しつつ計算時間を短縮することができる。さらに、上記式（３）および（４）を適用することで、空気等の流体と粒子２の密度差に因る重力の影響を考慮した濃度分布を提供することができる。
【００８５】
さらにまた、解析対象領域の複数点から拡散物質が空間に放出される場合についても、各点の予測結果を算出し、算出結果を積算する（重ね合わせる）ことで、解析対象領域の一点から拡散物質が空間に放出された場合と同様に濃度分布を求めることができる。
【００８６】
以上、本発明の実施形態に係る濃度分布解析装置および濃度分布解析方法によれば、市街地等の流れが大きく変動する環境下における拡散現象を解析する際に粒子追跡法を用いて、従来技術では困難であった流れの変動の影響が濃度分布に与える影響を高速に評価でき、迅速かつ高精度な濃度分布情報をユーザに提示することができる。すなわち、流れ場の流速の変動や乱流エネルギの変動の影響を考慮した濃度分布を迅速にユーザに提供することができる。
【００８７】
また、本発明の実施形態に係る濃度分布解析装置および濃度分布解析方法によれば、粒子２が濃度評価地点における濃度に影響を及ぼすか及ぼさないとみなすことができるかを判断し、影響を及ぼさないとみなすことができると判断した場合には当該粒子２を濃度計算する際の対象から除外するので、ユーザが望む濃度計算の精度を低下させることなく計算時間を短縮することができる。
【００８８】
さらに、本発明の実施形態に係る濃度分布解析装置および濃度分布解析方法によれば、上記式（３）および（４）を適用することで、空気と粒子２の密度差に因る重力の影響を考慮した濃度分布を提供することができる。さらにまた、解析対象領域の複数点から拡散物質が空間に放出される場合についても、解析対象領域の一点から拡散物質が空間に放出された場合と同様に濃度分布を求めることができる。
【００８９】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階では、上述した実施例以外にも様々な形態で実施することが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、追加、置き換え、変更を行なうことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【００９０】
１…解析メッシュ、２（２ａ〜２ｄ）…粒子、３…放出源、４…流速ベクトル、５…構造物、６…粒子位置の流速ベクトル、７…濃度評価地点の流速ベクトル、８…粒子位置と濃度評価地点を結ぶベクトル、１０…濃度分布解析装置、１１…入力部、１２…流体計算部、１３…粒子追跡部、１４…濃度分布計算部、１５…計算終了判定部、１６…表示部、２１…濃度分布算出処理部、２３…補正処理部、２４…濃度影響粒子選別処理部、２９…濃度評価地点Ｐ１を通過する流線。
【図１】
0005774520.tif 000022
【図２】
0005774520.tif 000023
【図３】
0005774520.tif 000024
【図４】
0005774520.tif 000025
【図５】
0005774520.tif 000026
【図６】
0005774520.tif 000027
【図８】
0005774520.tif 000028
【図７】
0005774520.tif 000029