[地球流体電脳倶楽部] [dcmodel | dcmodel-tools]
地表面熱フラックス計算のテスト. 風速, 地表面温度の値を固定して, 熱フラックスの計算を行う. 熱フラックス以外の計算部分を全てコメントアウトして計算を行った.
ソースプログラム | arare3m-20050716 |
基本場の温度 | 温位一様 (290 K) |
地表面の温度 | 300 K |
基本場の圧力 | 静水圧平衡 (地表面で 1.0 x 105 Pa) |
温度の初期擾乱 | なし |
圧力の初期擾乱 | なし |
風速の初期値 | u = 10 m/s , w = 0 m/s |
放射加熱 | なし |
地表面フラックス | バルク法(バルク係数は 0.01 で一定, 地表面温度 300 K で固定) |
音波の減衰 | あり (係数は 5.0 x 10-7) |
粘性 | あり |
時間フィルタ | あり (係数は 0.1) |
数値粘性 | あり (5.0 x 10-4) |
重力加速度 | 9.8 m/s2 |
計算領域 | 水平 10 km x 鉛直 10 km (解像度は計算によって異なる) |
積分時間 | 12 時間 |
時間ステップ | 3 秒 (ただし音波関連項についてはクーラン数が 0.5 程度になるようにした) |
境界条件 | 水平: 周期, 鉛直: free-slip |
数値解法 | 水平: 陽解法, 鉛直: 陰解法 |