[地球流体電脳倶楽部] [deepconv | 計算結果]

Odaka et al. (1998) の再計算

概要

Odaka et al. (1998) [参考文献 1] で行われている計算を deepconv/arare で行った. この計算で参照した Odaka et al. (1998) の計算設定は以下の通り.

物性値などは火星の値を用いる(定圧比熱 734.9 J/Kg/K, 気体定数 189.0 J/kg/K, 重力加速度 3.72 m/s^2).
地表面過程はバルク法を用いて計算する. 簡単のため熱と運動量についてのバルク係数は同じであるとし, 値は 0.01 で一定であるとする. また, 地表面温度は 270 K で固定.
高度 5 km 以下では放射冷却を模した一様冷却を考慮する. 冷却率は 50 K/day とする.
基本場の温度は地表(245 K)から 5 km (220 K)までは断熱的であり, それより上では等温とする.
基本場の圧力は静水圧平衡(地表面で 700 Pa).
初期の風速分布は 0.
境界条件は水平方向が cyclic, 鉛直方向が free-slip.

計算結果

2005-11-30-2	Odaka et al. (1998) と同じ解像度, 領域幅で 24 モデル時間計算. 宇宙研 SX-6 を用いて計算.
2005-11-30	Odaka et al. (1998) と同じ解像度, 領域幅で計算. 宇宙研 SX-6 を用いて計算.
2005-11-28	Odaka et al. (1998) に比べて解像度の高さを半分, 水平領域幅を半分にして計算. エネルギーの値を出力.
2005-11-14	Odaka et al. (1998) に比べて解像度の高さを半分, 水平領域幅を半分にして計算. エネルギーの値を出力.
2005-10-08	Odaka et al. (1998) に比べて解像度の高さを半分, 水平領域幅を半分にして計算. 乱流拡散係数の計算を修正. 基本場の拡散を考慮.
2005-10-05	Odaka et al. (1998) に比べて解像度の高さを半分, 水平領域幅を半分にして計算. 全質量の計算を修正.
2005-09-18	Odaka et al. (1998) と同じ領域幅, 解像度で計算. 質量・エネルギーの時間変化を出力.
2005-09-07-1	Odaka et al. (1998) と同じ領域幅, 解像度で計算. 宇宙研 SX-6 にて計算.
2005-09-06	Odaka et al. (1998) に比べて解像度の高さを半分, 水平領域幅を半分にして計算. 熱の式の各項を出力.
2005-08-11	Odaka et al (1998) と同じ解像度, 水平領域幅で計算. (宇宙研 SX-6 で計算.)
2005-08-09	Odaka et al. (1998) に比べて解像度の高さを半分, 水平領域幅を半分にして計算. (宇宙研 SX-6 で計算)
2005-08-02	Odaka et al. (1998) に比べて解像度の高さを半分, 水平領域幅を半分にして計算.

失敗例

2005-09-30	温位基本場の拡散を考慮した時の結果を詳しく調べるために 1 モデル時間だけ計算.
2005-09-17	Odaka et al. (1998) に比べて解像度の高さを半分, 水平領域幅を半分にして計算. 質量・エネルギーの時間変化を出力. 基本場温位の拡散を計算するようにしたら温位分布が階段状になった.
2005-09-07	Odaka et al. (1998) に比べて解像度の高さを半分, 水平領域幅を半分にして計算. 基本場温位の拡散を計算するようにしたら温位分布が階段状になった.
2005-09-03	Odaka et al. (1998) に比べて解像度の高さを半分, 水平領域幅を半分にして計算. 時間とともに対流層がつぶれていく. 原因はよくわからない.
2005-08-08	Odaka et al. (1998) と同じ解像度で計算. 長い時間ステップが長すぎたため, 乱流拡散係数に対する CFL 条件が破れている. 6 時間経過すると計算が不安定になる.
2005-07-24	熱フラックスを計算するところで温度フラックスを計算してしまっていた. そのため一様冷却と吊会うだけの地表面熱フラックスができずに冷え続ける.

参考文献

Odaka, M., Nakajima, K., Takehiro, S., Ishiwatari, M., Hayashi, Y.-Y., 1998, A numerical study of the Martian atomospheric convection with a two-dimentional anelastic model., Earth Planet Space, 50, 431--437.

deepconv Developers Group:

2005/12/07 (北守太一) 最終更新
2005/08/15 (北守太一) 新規作成