计算流体工学test.pptVIP

計算流体工学大気環境シミュレーションの現状 佐藤正樹 埼玉工業大学／地球フロンティア研究センター satoh@sit.ac.jp 内容 大気環境シミュレーションの概観 数値計算手法 非静力学モデル：領域モデル 大気大循環モデル：全球モデル 　　気候シミュレーションの現状 気候モデリングにおけるコンピュータ環境 次世代大気大循環モデル 大気環境シミュレーションの概観 大気環境系の概観 気候系 　大気、海洋、陸面（雪氷、土壌、植生）、海氷、生態系 空間スケール 　地球 10000km ～ 積雲 1km～乱流 10m 時間スケール 地球の歴史 数10億年～第四紀　数100万年 人為的気候変化 100年 10年スケール気候変化、エルニーニョ ３年 年周期、季節変化、日変化、短期予報（数時間） 気候システム 数値シミュレーションによる気候研究 大気の空間スケール 気候変動の時間スケール 大気の数値モデル メソモデル：雲解像モデル 　100km×100km：非静力学モデル 領域モデル 1000km×1000km 領域：静力／非静力学モデル 大気大循環モデル 　全球3次元： 静力学モデル 数値モデルの格子数と分解能 格子数 　　数100×数100×数10=数105 　　? 1000×1000×100=108 （次世代：ES） 領域幅と格子間隔： 　大循環 10000km：Δｘ＝100km ? 10km 領域　　1000km：Δｘ＝10km ? 1km メソ　　　100km：Δｘ＝1km 大気モデルの構成（領域?メソモデル?大気大循環モデル） 力学過程 方程式系 空間差分法，時間差分法 数値拡散 物理過程 水循環過程（雲物理，積雲パラメタリゼーション） 放射過程（太陽放射，赤外放射） 乱流過程（接地境界層，subgrid乱流） 化学変化?物質輸送（オゾン，エアロゾル） 数値計算手法 数値計算手法 方程式系 大気の基本バランスと波動 ＣＦＬ条件 音波?重力波の扱い 移流計算 擬スペクトル法?変換法 流体を記述する方程式系 非圧縮ブシネスク系：工学or理論 完全圧縮性成層流体 非静力学方程式 プリミティブ方程式系 静力学近似 浅水波方程式系 準地衡風方程式系 大気の基本バランス 静力学平衡 地衡風平衡 工学における流体の方程式系（１） 非圧縮の式 工学における流体の方程式系（２） ブシネスク方程式系 大気の方程式系 圧縮性成層流体の方程式系 大気の方程式系 回転系の圧縮性成層流体の方程式系 波動 音波：圧縮性　 　 重力波：成層効果：浮力 ロスビー波：緯度方向の回転角速度の変化 音波 重力波 CFL（Courant-Friedrichs-Lewy）条件 移流方程式 差分化,　von　Neumann　解析 大気モデルのCFL条件 Δx=10～100km，Δz=1km 音波 　ΔtΔx/cs=30s～5min：ΔtΔz/cs=3～30s 重力波 　Δt 1/N = 10min 移流：U=30m/s，W=数cm/s (メソ：W=10m/s) Δt Δx/U=3～30min 　　(ΔtΔz/W=100s) ロスビー波：移流と同等 音波の差分解法 非静力学モデル （メソモデル?積雲解像モデル） 代表的な非静力学モデル RAMS　 MM5 /mm5 ARPS WRF JMA-NHM (気象庁) 　　http://www.mri-jma.go.jp/Project.mrinpd/INDEXJ.htm http://pfi.kishou.go.jp (気象庁プラットフォーム) CReSS (名大) http://www.rain.ihas.nagoya-u.ac.jp/CReSS　 JMA-NHM：気象庁非静力学モデル 非静力学モデルの力学フレーム 音波の扱いによって方程式系が分類される 非弾性系 弾性系　or　完全圧縮系 陽解法（HEVE法） 完全陰解法（HIVI法） 水平陽解法、鉛直陰解法（HEVI法） 水平伝播する音波、重力波を分離 (split-explicit) 陽解法（HEVE法） 大気モデルは通常水平格子間隔の方が鉛直格子間隔より一桁大きい 水平格子間隔　Δx = 10km 　CFL条件?Δt = 10km/300m/s = 30s 鉛直格子間隔 Δz = 500m 　CFL条件?Δt = 500m/300m/s = 1.7s ?陽解法（HEVE法）は実用的でない。 陰解法（HIVI法）、非弾性系 ポアソン方程式の解法 緩和法：Jacobi法, Gauss-Seidel法, SOR法 共役勾配法 スペクトル法