[地球流体電脳倶楽部] [dcmodel | deepconv | プログラムのテスト計算[deepconv/arare5] ]

deepconv/arare4 重力流のベンチマーク

deepconv/arare5 による Straka et al. (1993) の再計算(2D).

参考文献

• Straka, J. M., Wilhelmson, R. B., Wicker, A. L., Anderson, J. R., Droegemeier, K. K., 1993 : Numerical solutions of a non-linear density current : A benchmark solution and comparisons. International Journal for Numerical Methods in Fluids , 17, 1--22

計算設定

• 計算領域
  • 水平 51.2 km, 鉛直 6.4 km
• 基本場の設定
  • 等温位大気 (θ = 300K)
• 物理過程の設定
  • 拡散係数一定 (75 m^2 s^-1)
  • 数値拡散なし
  • スポンジ層なし
• 大気成分等のパラメータ
  • 乾燥成分のみ
  • g = 9.81 m s^-2
  • Ps = 10^5 Pa
  • Rd = 287.0 J kg^-1 K^-1
  • Cp = 1004 J kg^-1 K^-1
  • Cv = 717.0 J kg^-1 K^-1
• 空間/時間格子間隔
  • Δx = Δ z = 533, 400, 266, 200, 133, 100, 66, 50, 33, 25 (m)
  • Δt(long) = 5.0, 4.0, 3.0, 2.0, 1.5, 1.0, 0.75, 0.5, 0.4, 0.25 (sec)
  • Δt(short) = 0.5, 0.4, 0.3, 0.2, 0.15, 0.1, 0.075, 0.05, 0.04, 0.025 (sec)
• 境界条件
  • 水平方向に周期境界, 上下境界で w = 0
• 初期条件
  • u = v = w = π = 0
  • 温位擾乱は θ(x,z) = A [cos(πL) + 1]*0.5 ( L < 1.0 ) or 0.0
    • A = -15 K, L = [{(x - x_c)/x_r}^2 + {(z - z_c)/z_r}^{2}]^0.5
    • x_c = 25.6 km, x_r = 4.0 km, z_c = 3.0 km, z_r = 2.0 km
• 積分時間
  • 900 sec

ソースコード, 設定ファイル

• ソースコード
  • tgz ファイル
• 設定ファイル
• 計算データ
• 描画スクリプト

結果一覧

解像度毎の結果(温位)

	t = 0	t = 300	t = 600	t = 900
Δx = Δz = 25 m
Δx = Δz = 33 m
Δx = Δz = 50 m
Δx = Δz = 66 m
Δx = Δz = 100 m
Δx = Δz = 133 m
Δx = Δz = 200 m
Δx = Δz = 266 m
Δx = Δz = 400 m
Δx = Δz = 533 m