[地球流体電脳倶楽部] [dcmodel | dcmodel-tools] [dcmodel-thum.rb]

deepconv/arare4 : Kelvin-Helmholtz 不安定計算(スポンジ層なし, tanh 型の初期水平流速)

arare4 で Kelvin-Helmholtz 不安定の計算が行えるか調べた.

• 空間格子間隔は水平 10 m, 鉛直 10 m.
• 水平方向の境界は周期境界.

計算設定

• 計算領域
  • 水平 2000 m
  • 鉛直 250 m
• 基本場の設定
  • 地表気温 300 K, 地表気圧 1000 hPa
  • tanh 型の初期水平流速(シアーの上下で 8.0 m/s の速度差)
• 物理過程の設定
  • 熱に対する乱流拡散係数は運動量に対する拡散係数と同じ
    • 最大値は 75 m^2/sec に設定.
  • 放射強制なし
  • スポンジ層なし
  • ニュートン冷却なし
• 大気成分
  • 乾燥成分 N2
  • 凝結成分なし (相対湿度 0)
• 空間/時間格子間隔, 境界条件
  • 空間格子間隔
    • 格子点数は 200 x 50
  • 時間格子間隔
    • 長い時間ステップ: 0.05 sec
    • 短い時間ステップ: 0.01 sec
  • 境界条件
    • 水平境界: 周期境界条件
    • 上下境界: 応力なし, w=0 , 熱フラックスなし
• 初期条件
  • sin 型の温位擾乱
  • 振幅は 0.01 K
• 積分時間
  • 1000 sec

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

• ソースコード
  • main プログラムはarare.f90 を書き換えた arare_SK1994.f90.
• arare4-20071019 において以下のファイルを書き換えて計算を行なった.
  • ネームリストファイルarare.conf
  • 基本場を与えるファイルbasicenv.f90
  • 擾乱を与えるファイルdisturbenv.f90

結果一覧

Km_0000_thumb.png	Km_0100_thumb.png	Km_0200_thumb.png	Km_0300_thumb.png	Km_0400_thumb.png	Km_0500_thumb.png	Km_0600_thumb.png	Km_0700_thumb.png	Km_0800_thumb.png	Km_0900_thumb.png	Km_1000_thumb.png
Exner_0000_thumb.png	Exner_0100_thumb.png	Exner_0200_thumb.png	Exner_0300_thumb.png	Exner_0400_thumb.png	Exner_0500_thumb.png	Exner_0600_thumb.png	Exner_0700_thumb.png	Exner_0800_thumb.png	Exner_0900_thumb.png	Exner_1000_thumb.png
PotTemp_0000_thumb.png	PotTemp_0100_thumb.png	PotTemp_0200_thumb.png	PotTemp_0300_thumb.png	PotTemp_0400_thumb.png	PotTemp_0500_thumb.png	PotTemp_0600_thumb.png	PotTemp_0700_thumb.png	PotTemp_0800_thumb.png	PotTemp_0900_thumb.png	PotTemp_1000_thumb.png
VelX_0000_thumb.png	VelX_0100_thumb.png	VelX_0200_thumb.png	VelX_0300_thumb.png	VelX_0400_thumb.png	VelX_0500_thumb.png	VelX_0600_thumb.png	VelX_0700_thumb.png	VelX_0800_thumb.png	VelX_0900_thumb.png	VelX_1000_thumb.png
VelZ_0000_thumb.png	VelZ_0100_thumb.png	VelZ_0200_thumb.png	VelZ_0300_thumb.png	VelZ_0400_thumb.png	VelZ_0500_thumb.png	VelZ_0600_thumb.png	VelZ_0700_thumb.png	VelZ_0800_thumb.png	VelZ_0900_thumb.png	VelZ_1000_thumb.png