<< x 方向の温位の移流 >>

x 軸方向の温位の移流を計算する. 4 次精度中心差分演算ルーチンおよび 2 次精度中心差分演算ルーチンの 2 種類の差分ルーチンを用いて結果を比較する.
空間解像度を固定して, 長い時間ステップをいろいろ変えて計算する.
ここではタイムフィルターは入れる.
プログラムソース arare2
基本場の温度 温度一様(298.15 K)
音速 347 m/s
x 方向の風速 35 m/s
音波の減衰 なし
粘性 なし
重力 なし
定圧比熱 1004 J K^-1 kg^-1
積分時間 300 秒
計算領域 水平(x) : 10000 m
鉛直(z) : 10000 m
境界条件 水平(x) : 周期境界条件
鉛直(z) : 固定境界条件
数値解法 水平(x) : 陽解法
鉛直(z) : 陰解法
アッセリンの時間フィルタの係数 0.1


[HOME]

 


4 次精度中心差分を用いた計算


長い時間ステップ 1.5 s
クーラン数 0.26

arare_dim_ss-adv4_center4



arare_dim_ss-adv4_center4-t0



arare_dim_ss-adv4_center4-t100



arare_dim_ss-adv4_center4-t200



arare_dim_ss-adv4_center4-t300



長い時間ステップ 3 s
クーラン数 0.52

arare_dim_ss-adv5_center4



arare_dim_ss-adv5_center4-t0



arare_dim_ss-adv5_center4-t100



arare_dim_ss-adv5_center4-t200



arare_dim_ss-adv5_center4-t300



長い時間ステップ 6 s
クーラン数 1.04

arare_dim_ss-adv6_center4



arare_dim_ss-adv6_center4-t0



arare_dim_ss-adv6_center4-t100



arare_dim_ss-adv6_center4-t200



arare_dim_ss-adv6_center4-t300




2 次精度中心差分を用いた計算


長い時間ステップ 1.5 s
クーラン数 0.26

arare_dim_ss-adv4_center2



arare_dim_ss-adv4_center2-t0



arare_dim_ss-adv4_center2-t100



arare_dim_ss-adv4_center2-t200



arare_dim_ss-adv4_center2-t300



長い時間ステップ 3 s
クーラン数 0.52

arare_dim_ss-adv5_center2



arare_dim_ss-adv5_center2-t0



arare_dim_ss-adv5_center2-t100



arare_dim_ss-adv5_center2-t200



arare_dim_ss-adv5_center2-t300



長い時間ステップ 6 s
クーラン数 1.04

arare_dim_ss-adv6_center2



arare_dim_ss-adv6_center2-t0



arare_dim_ss-adv6_center2-t100



arare_dim_ss-adv6_center2-t200



arare_dim_ss-adv6_center2-t300


Last update: Tue Apr 19 14:40:32 JST 2005
Copyright © deepconv 2005 (北守太一)