|
1 2 3 4 5 6 4. 3次元流体シミュレーション・モデルと手法 用いている3次元流体シミュレーションの手法のうち拡散モデルは、定常、非正規モデルです。 いわゆる運動方程式の偏微分方程式を差分法により解を求める数値解析を用いています。コンピュータ計算モデルは、風洞模型実験により精度を検証した3次元流体力学数値解析モデル(有限差分法)です。 出典:環境総合研究所(東京都目黒区) 以下はすうち数値解析モデルによる3次元流体シミュレーションの部分についての紹介です。数値解析モデルでは、解析解モデルのように、地形などの重要な要因を捨象する(すてさる)ことで偏微分方程式の解を求めることはできません。そのため繰り返し計算によりより解に近づく方法が採用されています。私達の3次元流体シミュレーションモデルでは、そのうちの一つである差分法を用いています。 出典:青山貞一、鷹取敦 山間地における大気拡散モデル による現況再現性の比較検証 東京都市大学環境情報学部紀要 第十号 (2009年2月28日発行) http://www.yc.tcu.ac.jp/~kiyou/no10/1-01.pdf 以下は偏微分方程式(波動方程式ないし運動方程式)を差分法の解説です。 微分方程式の差分解法の基礎 数値計算スキームとしての差分型非線形シュレーディンガー方程式 差分法では、計算するグリッド(メッシュ)をx、y、zいずれの方向とも、等間隔、たとえば500mとする場合と、主要部分を細かくし、そこから遠い部分については、計算時間を少なくするために荒くする、たとえば1500mとする方法が用いられます。これを有限差分法と言います。 この差分法を用いた数値計算モデルは、簡易モデルである解析解モデル(通称、プリュームモデル)では不可能な詳細に地形を考慮した風と汚染の流れをシミュレーションすることが出来ます。計算手法としては、有限差分法を用いています。しかし、計算量は簡易モデルに比べて数10万倍から数100万倍の量となり、膨大な計算時間がかかります。 以下は国のSPEEDIにおけるシミュレーションの手順です。環境総研の場合も、放出情報→風速場計算→濃度計算→線量計算→図形表示という手順は基本的に変わりません。 出典:特殊法人 原子力研究所部内資料 下図は先の手順に関連し、使用するデータファイルとの関連を示したものです。いずれも出典は特殊法人 原子力研究所部内資料です。 出典:特殊法人 原子力研究所部内資料 以下は、有限差分法を用いた3次元流体シミュレーションの手順です。まず風の流れを計算し、その後、風場データを使用して汚染の流れをシミュレーションします。 出典:環境総合研究所(東京都目黒区) 以下は環境総合研究所が流体力学モデルを用いた2次元流体シミュレーションの具体例です。中央に発生源(煙突高50m)がある場合、風下1kmに汚染が移流、拡散、沈降して様子をシミュレーションしています。風下に山や谷がある場合だけでなく、風上に山や谷がある場合にも拡散に影響が及びます。 出典:環境総合研究所(東京都目黒区) 以下は環境総合研究所が流体力学モデルを用いた3次元流体シミュレーションの具体例です。丘の上に産業廃棄物の焼却炉がある場合、風下に汚染が移流、拡散、沈降して様子をシミュレーションしています。汚染の移流、拡散、沈降が地形により非常に複雑に展開している様子が分かります。 出典:環境総合研究所(東京都目黒区) 以下はシステム開発に関連した東京新聞の記事です。 ◆この人 青山貞一さん 放射性物質拡散予測 独自のシステム開発 東京新聞 さらに、システム研究開発した背景と技術面については、以下の動画ごご覧ください。 パソコンでSPEEDI 研究開発の背景 青山貞一(環境総合研究所顧問) https://www.youtube.com/watch?v=wFG8tx6cPQw 下は地形を考慮した3次元シミュレーションの一例です。以下は東北電力の女川原発を事例としています。これは青山の女川講演会で使用したパワーポイントです。 出典:青山貞一 講演時使用パワーポイント 以下は宮城県女川町で2015年7月5日に女川体育館で行われた講演会に関連したNHKニュースです。上記のシミュレーション結果は、この講演会で公開されました。
つづく |