mecモデルの開発及び海域浄化装置の効果に関する数値.ppt

３．ヴァリアブルバージョンの応用例３．１　五ヶ所湾の海域浄化装置の効果に関する　　　　数値シミュレーション 東京大学　　　　　　　　 殿城賢三 佐藤　徹 密度流拡散装置（資料提供マリノフォーラム21） 密度拡散装置の位置 本研究の目的 MECモデルの検証 拡散装置による五ヶ所湾に対する効果を数値シミュレーションにより解析する 　-　温度?塩分分布の変化 　-　低層水の移動 計算条件(1) 密度流拡散装置のシミュレーションを行うに当たって計算時間中，成層を保つ必要がある 計算条件(2) 計算格子数；45*35*10（=15750；静水圧モデル） 　　　　　　　　　 35*35*21（=25725；Full-3Dモデル） dt ; 0.4(sec.)　???　（2つのモデルで同じｄｔを用いた） 開境界;　流速：等値外挿　　温度?塩分：観測値 　　　　　　　　潮位; 1998年7月31日からの潮位を与えた　 密度拡散装置 　　吸込み口；深度1.2m(1) 　　　　　　　　　及び11.4ｍ(2)に設置 　　吹出し口；深度2.7m(3)に設置 　　流量；(1)-0.2778（m3/s） (2) -1.111 （m3/s） (3) 0.3472 （m3/s）＊４（面） 計算格子(海面における格子配置図) 検証法－ 装置による影響についての検証 検証　　-温度 装置稼動後5時間後 検証　　-塩分 装置稼動後5時間後 検証　　-流速 装置から湾奥方向への流れ成分(湾奥方向へ5m) 計算結果(Full-3D温度分布) 計算結果(Full-3D塩分分布) 計算結果(温度分布) 計算結果(塩分分布) 計算結果(トレーサー分布-Full3D) 計算結果(トレーサー分布) 実験との比較 計算結果(圧力差) 結論 観測値との検証を行いMECモデルの検証を行った。 低層に存在する有機物は、密度流拡散装置により中層で吐出されると，密度流によりその深さを保ったまま、主に潮汐により4日ほどで湾奥まで運搬される。 今後の展望 　 海水流動 問題点（FULL-3D領域の大きさと減衰領域） 水平渦動粘性係数 KH ; リチャードソンの4/3乗則 Nestingと減衰領域 Full-3Dモデル周辺の計算格子 海洋の渦（乱流）と格子解像度 基礎方程式(Full-3D) 数値計算法 変数配置 ； スタッガード配置 格子 　 ；鉛直方向はバリアブル格子 　　　　　　　　　 水平方向は直交格子 アルゴリズム；MAC法 時間差分 ；陽解法 空間差分 ；3次上流（移流項）、 　 2次中心差分（拡散項） 圧力解法 ；ポアソン方程式をSORにより解く 渦動粘性係数 水平方向；max( KH,(nSGS)) 鉛直方向；max( KV, ( nSGS)) KH ; リチャードソンの4/3乗則 KV ; 成層化関数 　 nSGS ; SGS (SubgridScale)モデル 移動格子 Full-3Dモデルでは潮位を，格子幅を計算中に変更することができる「移動格子」を用いて表現している 境界条件 水平面内 　流速：接合プログラム 　温度?塩分：接合プログラムと等値外挿の選択 　圧力：等値外挿 鉛直面内 　流速：風応力摩擦応力を考慮した流速を外挿 　温度?塩分：気象条件(海面)，等値外挿(海底) 　圧力：等値外挿 MECモデルについて-時間的接合 MECモデルについて-空間的接合 2-3，境界条件1（海面，海底） 2-3，境界条件 2 温度?塩分 　水平面のダミーグリッドには?接合プログラム?および「Full-3Dモデル」により求められた値のうちで上流側の値を選択 計算格子(Full-3Dモデル周辺) 検証法 検証法 検証法 検証 塩分濃度???装置稼動後5時間後の検証 検証 温度???装置稼動後5時間後の検証 初期条件の生成 気象条件に尾鷲の値を用いた結果(塩分) 初期条件の生成 気象条件に尾鷲の値を用いた結果(温度) 検証法　-塩分濃度変化 装置による塩分分布の変化のDDSを用いた概念図 MECモデルについて-空間的接合 検証 装置からの流れ(湾口方向へ5m) 基準水面 Full-3D 計算領域 Full-3D 計算領域 格子移動 速度 静水圧モデル 計算領域 静水圧モデル 計算領域 水位 DZ dz DZ z z 時間 時間 静水圧 モデル Full-3D モデル 流速，温度，塩分，潮位 の時間的受け渡し図 ① ② ③-1 ② ④ ①‘ N(step) N+1(step) ③-ｎ ③-??? a) 水平面 b) 鉛直面