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

deepconv/arare4 : H2O の凝縮を考慮した木星雲対流

計算の概要

Nakajima et al (2000) の設定を従い, 木星大気での雲対流をシミュレートする. 凝縮成分は H2O のみ考慮し, 基本場の H2O 混合比の湿度は 100 km 付近より上空で 75% とする.

• 雲の併合成長過程として, Kessler (1969) を用いる. 但し, 閾値をゼロとし, 時間スケールを 100 秒とする.
• 浮力は, 温度と凝縮成分の基本場からのずれから評価.
• 最下層の混合比フラックスはバルク法で評価.
• 最下層の熱フラックスはバルク法で評価. F = Cd Cp ρ v ΔT (Cd = 1.5e-3, V = 100).

北大情報基盤センター SR11000 にて計算. 積分に必要な時間は約 4 時間.

計算設定

• 基本場の設定
  • 温度: 0 km (30 bar) 〜 170km (0.5 bar) まで等温位, 170km (0.5 bar) より上は温度一定
  • 相対湿度: 75% (100 km より上)
  • 放射強制: -1 K/day (高度 90km から 175 km)
  • スポンジ層: 領域の上層 100 km, 緩和時間 400 秒
• 大気成分
  • 乾燥成分: H2, He
  • 凝結成分: H2O
    • H2O: 相対湿度 75%
• 計算領域と時間ステップ
  • 領域:
    • 水平 512 km,
    • 鉛直 300 km (下端で 30 bar)
  • 解像度: 2.0 km (水平・鉛直とも)
  • 時間ステップ
    • 長い時間ステップ: 4 秒
    • 短い時間ステップ: 0.4 秒 (音波に対するクーラン数は約 0.2)
  • 境界条件
    • 水平境界: 周期境界条件
    • 上下境界: 応力なし, w=0
    • 温度フラックス: バルク法
    • 混合比フラックス: バルク法
• 初期条件
  • 100 kmに, 最大 0.5 k の温度擾乱を設定.

設定ファイル

arare.conf

備忘録

結果一覧

PotTempInit.png 温位の基本場

TempInit.png 温度の基本場

MixRtInit_s1.png 水蒸気の基本場

Budjet_H2O.png H2O の密度の時間変化

BuoyMean.png 浮力の水平平均

DelMixRtMean_s1.png 飽和混合比と混合比の相対誤差の水平平均

TempAll_Mean.png 温度の時間変化. 4e5 秒間隔

PotTempAll_Mean.png 温位の時間変化. 4e5 秒間隔

VPotTempAll_Mean.png 仮温位の時間変化. 4e5 秒間隔

MixRtAll-H2O_Mean_lin.png 混合比の時間変化(linear). 4e5 秒間隔

MixRtAll-H2O_Mean_log.png 混合比の時間変化(log). 4e5 秒間隔

PotTempTendency_t1.png 温位変化に対する各項の寄与(t=2e5) 黒:全体, 赤：移流、緑：放射、紫：乱流、青：凝結と蒸発、水色：散逸加熱

PotTempTendency_t2.png 温位変化に対する各項の寄与(t=4e5) 黒:全体, 赤：移流、緑：放射、紫：乱流、青：凝結と蒸発、水色：散逸加熱

PotTempTendency_t3.png 温位変化に対する各項の寄与(t=6e5) 黒:全体, 赤：移流、緑：放射、紫：乱流、青：凝結と蒸発、水色：散逸加熱

PotTempTendency_t4.png 温位変化に対する各項の寄与(t=8e5) 黒:全体, 赤：移流、緑：放射、紫：乱流、青：凝結と蒸発、水色：散逸加熱

PotTempTendency_t5.png 温位変化に対する各項の寄与(t=1e6) 黒:全体, 赤：移流、緑：放射、紫：乱流、青：凝結と蒸発、水色：散逸加熱

PotTempTendency_t6.png 温位変化に対する各項の寄与(t=1.2e6) 黒:全体, 赤：移流、緑：放射、紫：乱流、青：凝結と蒸発、水色：散逸加熱

PotTempTendency_t7.png 温位変化に対する各項の寄与(t=1.4e6) 黒:全体, 赤：移流、緑：放射、紫：乱流、青：凝結と蒸発、水色：散逸加熱

PotTempTendency_t8.png 温位変化に対する各項の寄与(t=1.6e6) 黒:全体, 赤：移流、緑：放射、紫：乱流、青：凝結と蒸発、水色：散逸加熱

PotTempTendency_t9.png 温位変化に対する各項の寄与(t=1.8e6) 黒:全体, 赤：移流、緑：放射、紫：乱流、青：凝結と蒸発、水色：散逸加熱

PotTempTendency_t10.png 温位変化に対する各項の寄与(t=2e6) 黒:全体, 赤：移流、緑：放射、紫：乱流、青：凝結と蒸発、水色：散逸加熱

温度擾乱

anal_Temp_anim.gif

anal_Temp_t0000000.png

anal_Temp_t0100000.png

anal_Temp_t0200000.png

anal_Temp_t0300000.png

anal_Temp_t0400000.png

anal_Temp_t0500000.png

anal_Temp_t0600000.png

anal_Temp_t0700000.png

anal_Temp_t0800000.png

anal_Temp_t0900000.png

anal_Temp_t1000000.png

anal_Temp_t1100000.png

anal_Temp_t1200000.png

anal_Temp_t1300000.png

anal_Temp_t1400000.png

anal_Temp_t1500000.png

anal_Temp_t1600000.png

anal_Temp_t1700000.png

anal_Temp_t1800000.png

anal_Temp_t1900000.png

anal_Temp_t2000000.png

温位擾乱

PotTemp_anim.gif

PotTemp_t0000000.png

PotTemp_t0100000.png

PotTemp_t0200000.png

PotTemp_t0300000.png

PotTemp_t0400000.png

PotTemp_t0500000.png

PotTemp_t0600000.png

PotTemp_t0700000.png

PotTemp_t0800000.png

PotTemp_t0900000.png

PotTemp_t1000000.png

PotTemp_t1100000.png

PotTemp_t1200000.png

PotTemp_t1300000.png

PotTemp_t1400000.png

PotTemp_t1500000.png

PotTemp_t1600000.png

PotTemp_t1700000.png

PotTemp_t1800000.png

PotTemp_t1900000.png

PotTemp_t2000000.png

エクスナー関数

Exner_anim.gif

Exner_t0000000.png

Exner_t0100000.png

Exner_t0200000.png

Exner_t0300000.png

Exner_t0400000.png

Exner_t0500000.png

Exner_t0600000.png

Exner_t0700000.png

Exner_t0800000.png

Exner_t0900000.png

Exner_t1000000.png

Exner_t1100000.png

Exner_t1200000.png

Exner_t1300000.png

Exner_t1400000.png

Exner_t1500000.png

Exner_t1600000.png

Exner_t1700000.png

Exner_t1800000.png

Exner_t1900000.png

Exner_t2000000.png

水平風速度

VelX_anim.gif

VelX_t0000000.png

VelX_t0100000.png

VelX_t0200000.png

VelX_t0300000.png

VelX_t0400000.png

VelX_t0500000.png

VelX_t0600000.png

VelX_t0700000.png

VelX_t0800000.png

VelX_t0900000.png

VelX_t1000000.png

VelX_t1100000.png

VelX_t1200000.png

VelX_t1300000.png

VelX_t1400000.png

VelX_t1500000.png

VelX_t1600000.png

VelX_t1700000.png

VelX_t1800000.png

VelX_t1900000.png

VelX_t2000000.png

鉛直風速度

VelZ_anim.gif

VelZ_t0000000.png

VelZ_t0100000.png

VelZ_t0200000.png

VelZ_t0300000.png

VelZ_t0400000.png

VelZ_t0500000.png

VelZ_t0600000.png

VelZ_t0700000.png

VelZ_t0800000.png

VelZ_t0900000.png

VelZ_t1000000.png

VelZ_t1100000.png

VelZ_t1200000.png

VelZ_t1300000.png

VelZ_t1400000.png

VelZ_t1500000.png

VelZ_t1600000.png

VelZ_t1700000.png

VelZ_t1800000.png

VelZ_t1900000.png

VelZ_t2000000.png

乱流拡散係数

Km_anim.gif

Km_t0000000.png

Km_t0100000.png

Km_t0200000.png

Km_t0300000.png

Km_t0400000.png

Km_t0500000.png

Km_t0600000.png

Km_t0700000.png

Km_t0800000.png

Km_t0900000.png

Km_t1000000.png

Km_t1100000.png

Km_t1200000.png

Km_t1300000.png

Km_t1400000.png

Km_t1500000.png

Km_t1600000.png

Km_t1700000.png

Km_t1800000.png

Km_t1900000.png

Km_t2000000.png

H2O(g)の混合比. ログでプロットしている. 1e-8 kg/kg から 2e-2 kg/kg を表示.

anal_MixRtAll_s1_anim.gif

anal_MixRtAll_s1_t0000000.png

anal_MixRtAll_s1_t0100000.png

anal_MixRtAll_s1_t0200000.png

anal_MixRtAll_s1_t0300000.png

anal_MixRtAll_s1_t0400000.png

anal_MixRtAll_s1_t0500000.png

anal_MixRtAll_s1_t0600000.png

anal_MixRtAll_s1_t0700000.png

anal_MixRtAll_s1_t0800000.png

anal_MixRtAll_s1_t0900000.png

anal_MixRtAll_s1_t1000000.png

anal_MixRtAll_s1_t1100000.png

anal_MixRtAll_s1_t1200000.png

anal_MixRtAll_s1_t1300000.png

anal_MixRtAll_s1_t1400000.png

anal_MixRtAll_s1_t1500000.png

anal_MixRtAll_s1_t1600000.png

anal_MixRtAll_s1_t1700000.png

anal_MixRtAll_s1_t1800000.png

anal_MixRtAll_s1_t1900000.png

anal_MixRtAll_s1_t2000000.png

H2O 雲の混合比. ログでプロットしている. 1e-8 kg/kg から 1e-3 kg/kg を表示.

MixRt_s2_anim.gif

MixRt_s2_t0000000.png

MixRt_s2_t0100000.png

MixRt_s2_t0200000.png

MixRt_s2_t0300000.png

MixRt_s2_t0400000.png

MixRt_s2_t0500000.png

MixRt_s2_t0600000.png

MixRt_s2_t0700000.png

MixRt_s2_t0800000.png

MixRt_s2_t0900000.png

MixRt_s2_t1000000.png

MixRt_s2_t1100000.png

MixRt_s2_t1200000.png

MixRt_s2_t1300000.png

MixRt_s2_t1400000.png

MixRt_s2_t1500000.png

MixRt_s2_t1600000.png

MixRt_s2_t1700000.png

MixRt_s2_t1800000.png

MixRt_s2_t1900000.png

MixRt_s2_t2000000.png

H2O 雨の混合比. ログでプロットしている. 1e-8 kg/kg から 1e-3 kg/kg を表示.

MixRt_s3_anim.gif

MixRt_s3_t0000000.png

MixRt_s3_t0100000.png

MixRt_s3_t0200000.png

MixRt_s3_t0300000.png

MixRt_s3_t0400000.png

MixRt_s3_t0500000.png

MixRt_s3_t0600000.png

MixRt_s3_t0700000.png

MixRt_s3_t0800000.png

MixRt_s3_t0900000.png

MixRt_s3_t1000000.png

MixRt_s3_t1100000.png

MixRt_s3_t1200000.png

MixRt_s3_t1300000.png

MixRt_s3_t1400000.png

MixRt_s3_t1500000.png

MixRt_s3_t1600000.png

MixRt_s3_t1700000.png

MixRt_s3_t1800000.png

MixRt_s3_t1900000.png

MixRt_s3_t2000000.png

安定度の水平平均

anal_Stab_anim.gif

Stab_t0.png

Stab_t1.png

Stab_t2.png

Stab_t3.png

Stab_t4.png

Stab_t5.png

Stab_t6.png

Stab_t7.png

Stab_t8.png

Stab_t9.png

Stab_t10.png

Stab_t11.png

Stab_t12.png

Stab_t13.png

Stab_t14.png

Stab_t15.png

Stab_t16.png

Stab_t17.png

Stab_t18.png

Stab_t19.png

Stab_t20.png