<< 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 (北守太一)