[deepconv/arare/sample]

プリュームの上昇(6)

概要

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

Last Update: 2005/02/19 (小高正嗣)