プリュームの上昇する様子を MPDATA スキームを用いて計算する. この計算では基本場は温位一定かつ 静水圧平衡としている. 数値粘性は運動方程式にのみ導入する.
プログラムソース | deepconv/arare2 (2005/02/19) |
設定ファイル | arare.conf (2005/02/19) |
基本場の温度 | 温位一様(300 K) |
音速 | 350 m/s |
基本場の風速 | 0 m/s |
音波の減衰 | 5.0d-7 * ( DelX ** 2.0d0 ) / DelTimeShort |
乱流粘性拡散 | あり. |
重力 | 9.8 m/s^2 |
定圧比熱 | 1004 J K^-1 kg^-1 |
積分時間 | 3000 秒 |
計算領域 |
水平(x) : 10000 m 鉛直(z) : 10000 m |
格子点数 | 水平(x) : 50 鉛直(z) : 50 |
時間間隔 |
短いタイムステップ: 0.3 秒 長いタイムステップ: 3 秒 |
境界条件 |
水平(x) : 周期境界条件 鉛直(z) : 固定境界条件 |
数値解法 |
水平(x) : 陽解法 鉛直(z) : 陰解法 |
数値粘性係数 | 5.0d-2 * ( DelX ** 2.0d0 ) / DelTimeLong |
出力ファイル |
arare_fs-R50_TL3_TS03_plume6.nc
arare_sf-R50_TL3_TS01_plume6.nc arare_ss-R50_TL3_TS01_plume6.nc |
温位の擾乱 Time = 0 秒 | 温位の擾乱 Time = 300 秒 | 温位の擾乱 Time = 600 秒 |
Time = 900 秒 | Time = 1200 秒 | Time = 1500 秒 |
Time = 1800 秒 | Time = 2100 秒 | Time = 2400 秒 |
Time = 2700 秒 | Time = 3000 秒 |
数値粘性係数 | 5.0d-3 * ( DelX ** 2.0d0 ) / DelTimeLong |
出力ファイル |
arare_fs-R50_TL3_TS03_plume6a.nc
arare_sf-R50_TL3_TS01_plume6a.nc arare_ss-R50_TL3_TS01_plume6a.nc |
温位の擾乱 Time = 0 秒 | 温位の擾乱 Time = 300 秒 | 温位の擾乱 Time = 600 秒 |
Time = 900 秒 | Time = 1200 秒 | Time = 1500 秒 |
Time = 1800 秒 | Time = 2100 秒 | Time = 2400 秒 |
Time = 2700 秒 | Time = 3000 秒 |
数値粘性係数 | 5.0d-4 * ( DelX ** 2.0d0 ) / DelTimeLong |
出力ファイル |
arare_fs-R50_TL3_TS03_plume6b.nc
arare_sf-R50_TL3_TS01_plume6b.nc arare_ss-R50_TL3_TS01_plume6b.nc |
温位の擾乱 Time = 0 秒 | 温位の擾乱 Time = 300 秒 | 温位の擾乱 Time = 600 秒 |
Time = 900 秒 | Time = 1200 秒 | Time = 1500 秒 |
Time = 1800 秒 | Time = 2100 秒 | Time = 2400 秒 |
Time = 2700 秒 | Time = 3000 秒 |
数値粘性係数 | 0.0d0 * ( DelX ** 2.0d0 ) / DelTimeLong |
出力ファイル |
arare_fs-R50_TL3_TS03_plume6c.nc
arare_sf-R50_TL3_TS01_plume6c.nc arare_ss-R50_TL3_TS01_plume6c.nc |
温位の擾乱 Time = 0 秒 | 温位の擾乱 Time = 300 秒 | 温位の擾乱 Time = 600 秒 |
Time = 900 秒 | Time = 1200 秒 | Time = 1500 秒 |
Time = 1800 秒 | Time = 2100 秒 | Time = 2400 秒 |
Time = 2700 秒 | Time = 3000 秒 |