領域中央を中心とする定常な剛体回転流を与えた時の温位擾乱の移流の様子を計算する. この計算でやりたいことはガウシアンブンプで与えた温位擾乱が何周もした時に擾乱の様子がどのようになるかを調べることである.
※ しかし, どうも計算がうまくいかないため境界条件をいろいろ変えるとどうなるか計算してみたのがこの計算結果である.
| プログラムソース | arare1 |
| 格子点数 | 水平(x): 50 鉛直(z): 50 |
| 長い時間ステップ | 3 秒 |
| 短い時間ステップ | 0.3 秒 |
| 基本場の温度 | 温度一様(298.15 K) |
| 音速 | 347 m/s |
| 風速 |
x = 5000 m, z = 5000 m を中心とする回転流 (角速度: 5.0d-2 1/s) |
| 音波の減衰 | なし |
| タイムフィルター | あり (α=0.05) |
| 粘性 | なし |
| 重力 | なし |
| 定圧比熱 | 1004 J K^-1 kg^-1 |
| 積分時間 | 900 秒 |
| 計算領域 |
水平(x) : 10000 m 鉛直(z) : 10000 m |
| 数値解法 |
水平(x) : 陽解法 鉛直(z) : 陰解法 |
| 初期値ファイル | 基本場 : arare-basic-R50.nc 擾乱 : arare-disturb-R50.nc |
| 計算結果(内部領域) | スカラー : arare_ss-R50_TL3_TS03-bcCC.nc 水平速度 : arare_fs-R50_TL3_TS03-bcCC.nc 鉛直速度 : arare_sf-R50_TL3_TS03-bcCC.nc |
| 計算結果(全領域) | スカラー : arare_dim_ss-R50_TL3_TS03-bcCC.nc 水平速度 : arare_dim_fs-R50_TL3_TS03-bcCC.nc 鉛直速度 : arare_dim_sf-R50_TL3_TS03-bcCC.nc |
|
|
|
|
| 温位 (t = 0 s) | 温位 (t = 50 s) | 温位 (t = 100 s) |
|
|
|
| 温位 (t = 150 s) | 温位 (t = 200 s) | 温位 (t = 300 s) |
|
|
|
|
| 温位 (t = 400 s) | 温位 (t = 500 s) | 温位 (t = 600 s) |
|
|
|
|
| 温位 (t = 700 s) | 温位 (t = 800 s) | 温位 (t = 900 s) |
| 初期値ファイル | 基本場 : arare-basic-R50.nc 擾乱 : arare-disturb-R50.nc |
| 計算結果(内部領域) | スカラー : arare_ss-R50_TL3_TS03-bcRR.nc 水平速度 : arare_fs-R50_TL3_TS03-bcRR.nc 鉛直速度 : arare_sf-R50_TL3_TS03-bcRR.nc |
| 計算結果(全領域) | スカラー : arare_dim_ss-R50_TL3_TS03-bcRR.nc 水平速度 : arare_dim_fs-R50_TL3_TS03-bcRR.nc 鉛直速度 : arare_dim_sf-R50_TL3_TS03-bcRR.nc |
|
|
|
| 温位 (t = 0 s) | 温位 (t = 100 s) | 温位 (t = 200 s) |
|
|
|
| 温位 (t = 300 s) | 温位 (t = 400 s) | 温位 (t = 500 s) |
|
|
|
| 温位 (t = 540 s) | 温位 (t = 630 s) | 温位 (t = 720 s) |
|
||
| 温位 (t = 900 s) |