プリュームの上昇する様子を 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 秒 |