%     表題: 天文夏の学校 集録
%
%     履歴: 2000/10/18 杉山耕一朗

\begin{abstract}

木星表層はアンモニアの雲に覆われており, 雲の発生・消滅が観測されてい
る. しかしアンモニア雲層高度以下の領域は光学的にほとんど観測できないため
に, 木星の対流がどのような構造をしているかはよく分かっていない.
\bigskip

木星大気対流を実際にシミュレーションするためには, 大気の平均的な状態とし
て断熱温度減率と凝縮物質の鉛直分布を把握することが重要である. 
本研究では木星型惑星大気の断熱温度減率と凝縮物質の鉛直分布とを
計算するための手法を検討する. 
\bigskip

木星大気の相変化や化学反応を計算した例として Atreya and Romani (1985)
がある. 
彼らは熱力学関数 $S$ の保存式を理想気体の状態方程式と
潜熱・反応熱を用いて熱力学変数 $T$, $p$, n$_{i}$ の式
として定式化した. 
この定式化では空気塊において生じうる全ての化学反応を知る必要がある. 
しかし木星型惑星大気における化学反応の全容は明らかになっていない. 
組成の違う大気を考える場合には, 生じうる化学反応が変わるので, 
数値コードを書き換えねばならない. 
そもそも平衡状態を仮定しているので, 化学反応に関する情報がなくとも
断熱温度減率と凝縮物質の鉛直分布は計算できるはずである. 
\bigskip

本研究では大気中で生じうる反応式を考えずに済み, 
さらに大気組成を簡単に変更することのできる手法を検討する. 
それを行うためには
 (1) 大気の平衡状態の記述, (2) 大気の断熱変化の記述, を検討することが
必要となる. 
\bigskip

(1)大気の平衡状態の記述であるが, 
熱力学変数として温度 $T$, 圧力 
$p$, 物質存在量 $n_{i}$ を選択する. このとき大気の状態を
与える適切な熱力学関数は $S$ ではなく $G$ である. 
大気の平衡状態は, 熱力学関数 $G$ が最小化された状態であるとする. 
温度・圧力を与えると, 熱力学関数 $G$ は,
\vspace{-2mm} 
\begin{eqnarray}
G = \sum \mu_{i} n_{i},
\nonumber
\end{eqnarray}
と書ける. ここで $\mu_{i}$ は化学ポテンシャルである. 
上式より熱力学関数 $G$ の最小値は, 
物質存在量 $n_{i}$ をパラメータとして変化させることによって
求めることができる. 
大気の状態を熱力学関数 $G$ によって記述したので, 
潜熱・反応熱の代わりに化学ポテンシャル $\mu_{i}$ のリストを
与えればよい. 化学ポテンシャルは既存の物性データから
求めることができる. 
\bigskip

(2)大気の断熱変化の記述であるが, 大気の断熱変化は熱力学関数$S$ 
の保存として記述することができる. 熱力学関数 $S$ は $G$ を用いて
以下のように与えられる. 
\vspace{-2mm} 
\begin{eqnarray}
\left( \DP{G}{T} \right)_{p, n_{i}} = -S.
\nonumber
\end{eqnarray}
温度・圧力空間で $dS = 0$ の曲線の通る軌跡を順にたどれば, 
断熱温度減率と凝結物質の存在量を求めることができる. 
\bigskip

大気成分気体やその凝結成分を理想気体/液体/固体/溶液と仮定し, 上記の計算
手法に則って実際に木星大気における断熱温度減率と凝結物質の鉛直分布を計算
した. 大気組成はボ
イジャーの観測結果から予想されるものとした.  
凝結物質の種類と各凝結物質量のピーク値の圧力は, Atreya and Romani (1985) 
の結果にほぼ等しい.  したがって本研究によって開発された計算手法を用いる
ことによって, 他の木星型惑星に対しても妥当な断熱温度減率と凝結物質の鉛直
分布が得られると期待される.





\end{abstract}

