%修論 tex ファイル %作成開始日 2008/11/16 %作成者 山下 達也 %\documentclass[a4j,12pt]{jreport} \documentclass[a4j,12pt,twoside]{jreport} \usepackage{Dennou6} %電脳スタイルを書式として使う \usepackage{amsmath} \usepackage{times} \usepackage{url} \usepackage{bm} \usepackage[multiple]{footmisc} %同じ場所に複数脚注が入る場合，コンマを %入れる \pagestyle{Dmyheadings} %\Dauthor{山下 達也} %\Ddate{2009/02/02} %\Dtitle{火星湿潤対流(仮)} %\Dfile{} \setlength{\parskip}{2ex} \setcounter{chapter}{0} \def\thesection{\arabic{chapter}.\arabic{section}} \def\thesubsection{\arabic{chapter}.\arabic{section}.\arabic{subsection}} %\def\thesubsubsection{\arabic{chapter}.\arabic{section}.\arabic{subsection}.\ararbic{subsubsection}} %\def\theequation{\arabic{chapter}.\arabic{section}.\arabic{equation}} \def\theequation{\arabic{chapter}.\arabic{equation}} \def\thepage{} %\def\thepage{\arabic{page}} %\def\thetable{\arabic{chapter}.\arabic{section}.\arabic{table}} \def\thetable{\arabic{chapter}.\arabic{table}} \def\thefootnote{\arabic{footnote} )} \makeatletter \def\sech{\mathop {\operator@font sech}\nolimits} \def\cosech{\mathop {\operator@font cosech}\nolimits} \def\Tr{\mathop {\operator@font Tr}\nolimits} \def\Im{\mathop {\operator@font Im}\nolimits} \def\Re{\mathop {\operator@font Re}\nolimits} \def\Res{\mathop {\operator@font Res}\nolimits} \makeatother \begin{document} \Dtitle[大気主成分凝結対流]{ \Huge 大気主成分の凝結を考慮した 2 次元湿潤対流 \\ \Large --- 火星大気中の CO$_2$ 氷雲を伴う対流に関する \\ \Large 基礎的数値実験 --- \\ \Large Two dimensional moist convection including condensation of major atmospheric component \\ \large --- Basic numerical experiments on CO$_2$ ice cloud convection in Martian atmosphere --- } \Dauthor[山下 達也]{ \Huge 山下 達也 \\ YAMASHITA Tatsuya \\ 北海道大学理学院宇宙理学専攻 \\ 惑星物理学研究室修士課程 2 年 \\ \large Planetary Physics Laboratory, Department of Cosmosciences, \\ \large Faculty of Science, Hokkaido University } \Ddate[2009/02/02]{2009 年 2 月 2 日} \Dfile{} \maketitle \markboth{支配方程式} {支配方程式} \chapter{支配方程式} 第 2 章では雲対流モデル deepconv/arare4 で用いられている支配方程式につい て説明する. \section{大気主成分の凝結を考慮した 2 次元準圧縮方程式系} 本モデルでは支配方程式系として準圧縮方程式系(Klemp and Wilhelmson, 1978)に 大気主成分の凝結効果を考慮したもの(Odaka et al., 2005) を用いる. 準圧縮方程式系とは弾性方程式系のうち, 基本場を導入して線形近似を行うもの を指す. 弾性方程式系とは大気の圧縮性を考慮する方程式系を指す. 準圧縮方程式は更に線形近似の適用の有無により, 準圧縮方程式系と完全圧縮方 程式系に分類される. 準圧縮方程式では時間積分にモード別時間分割法(Klemp and Wilhelmson, 1978) を用いることが出来, 計算のコストを下げることが出来る\footnote{ モード別時間分割法の詳細については第 3 章で述べる. }. 一方, 完全圧縮方程式は準圧縮方程式よりも精度の点で優れているものの, モー ド別時間分割法を適用すると安定に計算することが出来ないことが指摘されてい る(MacDonald et al., 2000). そこで計算の効率性を重視して準圧縮方程式系を採用することとする. また本モデルでは水平 1 次元, 鉛直 1 次元の 2 次元直線直交座標系を採用し た. 本研究では主成分凝結を考慮した対流の性質を理解することを目的としている為, 一般に 3 次元よりも現象の解析を行いやすい 2 次元で議論を行うこととする. 以下, 主成分凝結を考慮した準圧縮方程式系について説明する. 各変数は基本場成分と擾乱成分の和で表現できると仮定する. 擾乱成分は基本場成分に比べて十分小さいと仮定し, 各物理量の時間発展の方程 式における 2 次の微小項を無視する. 但し本研究では対流に着目した計算を行なうので, 移流項については 2 次の微 小量を無視しない. 線形化を行なっている準圧縮方程式系では質量保存及びエネルギー保存が成 立しないことに注意しなければならない \footnote{主成分凝結を考慮した準圧縮方程式系でのエネルギー方程式の導出及 びエネルギー非保存の議論については付録 A を参照されたい. }. 準圧縮方程式系では水平流速 $u$, 鉛直流速 $w$, 温位 $\theta$, エクスナー 関数 $\Pi$, 雲密度 $\rho_s$ を予報し, 圧力 $p$, 密度 $\rho$, 温度 $T$ などの他の物理量については予報変数を用いて診断的に算出する. 温位, エクスナー関数は以下のように定義される. \begin{eqnarray} \theta &\equiv& T \left( \frac{p_0}{p} \right)^{\frac{R}{c_p}}, \Deqlab{Pottemp} \\ \Pi &\equiv& \left( \frac{p}{p_0} \right)^{\frac{R}{c_p}} = \frac{T}{\theta}. \Deqlab{Exner} \end{eqnarray} ここで $p_0$, $R$, $c_p$ はそれぞれ基準圧力, 単位質量あたりの気体定数, 定圧比熱を表し, $R \simeq 188.9$ [J K$^{-1}$ kg$^{-1}$], $c_p \simeq 734.1$ [J K$^{-1}$ kg$^{-1}$] である. 本研究では地表面圧力を基準圧力にとることにし, $p_0 = 700$ [Pa] とする. 本研究においては静止状態にある乾燥した水平一様な基本場を与える. 物理量 $\phi$ の基本場成分を $\overline{\phi}$, 擾乱成分を $\phi^{\prime}$ と表すことにすると, 各予報変数は以下のように表される. \begin{eqnarray} u(x,z,t) &=& u^{\prime}(x,z,t), \Deqlab{ZA} \\ w(x,z,t) &=& w^{\prime}(x,z,t), \Deqlab{ZB} \\ \theta(x,z,t) &=& \theta^{\prime}(x,z,t) + \overline{\theta}(z), \Deqlab{ZC} \\ \Pi(x,z,t) &=& \Pi^{\prime}(x,z,t) + \overline{\Pi}(z), \Deqlab{ZD} \\ \rho_s(x,z,t) &=& {\rho_s}^{\prime}(x,z,t). \Deqlab{ZE} \end{eqnarray} 流速及び雲密度の基本場成分はゼロであるので, 以後流速及び雲密度の擾乱成分 については簡単の為にプライム記号 `` $^{\prime}$ '' を省略する. 基本場は静止状態にあるので, 静水圧平衡の式 \begin{equation} \DP{\overline{p}}{z} = - \overline{\rho} g \Deqlab{hydrostat1} \end{equation} が成り立つ. 但し $g$ は重力加速度である. \Deqref{hydrostat1} を $\overline{\Pi}$, $\overline{\theta}$ で表現する と \begin{equation} \DP{\overline{\Pi}}{z} = - \frac{g}{c_p \overline{\theta}} \Deqlab{hydrostat2} \end{equation} となる. 簡単の為大気は CO$_2$ のみから成るものとし, 理想気体とみなせるものとする. また惑星の回転効果を無視する. このとき主成分凝結を考慮した準圧縮方程式は以下のように記述される. \begin{eqnarray} \DP{u}{t} &=& - u \DP{u}{x} - w \DP{u}{z} - c_p \overline{\theta} \DP{\Pi^{\prime}}{x} + D_u, \Deqlab{AA} \\ \DP{w}{t} &=& - u \DP{w}{x} - w \DP{w}{z} - c_p \overline{\theta} \DP{\Pi^{\prime}}{z} + D_w + g \frac{\theta^{\prime}}{\overline{\theta}}, \Deqlab{AB} \\ \DP{\theta^{\prime}}{t} &=& - u \DP{\theta^{\prime}}{x} - w \DP{\theta^{\prime}}{z} - w \DP{\overline{\theta}}{z} + \frac{\overline{\theta}}{\overline{T}} (Q_{cond} + Q_{rad} + Q_{dis}) + D_{\theta}, \Deqlab{AC} \end{eqnarray} \begin{eqnarray} \DP{\Pi^{\prime}}{t} &=& - \frac{\overline{c}^2}{c_p \overline{\rho} \overline{\theta}^2} \left[ \DP{ (\overline{\rho} \overline{\theta} u) }{x} + \DP{ (\overline{\rho} \overline{\theta} w) }{z} \right] + \frac{\overline{c}^2}{c_p \overline{\theta}^2} \DD{\theta^{\prime}}{t} - \frac{\overline{c}^2}{c_p \overline{\rho} \overline{\theta}^2} M_{cond}, \Deqlab{AD} \\ \DP{\rho_s}{t} &=& - \DP{(\rho_s u)}{x} - \DP{(\rho_s w)}{z} + M_{cond} + M_{fall} + D_{\rho_s}. \Deqlab{AE} \end{eqnarray} ここで $\overline{c}$ は音速の基本場成分を表し, 以下の関係式により定めら れる. \begin{equation} {\overline{c}}^2 = \frac{c_p}{c_v} R \overline{T} = \frac{c_p}{c_v} R \overline{\theta} \overline{\Pi}. \Deqlab{Sound} \end{equation} 但し $c_v$ は定積比熱を表し, $c_v \simeq 545.2$ [J K$^{-1}$ kg$^{-1}$] である. $D_u$, $D_w$, $D_{\theta}$, $D_{\rho_s}$ はそれぞれ水平流速, 鉛直流速, 温位, 雲密度に関する乱流拡散項である. $Q_{cond}$, $Q_{rad}$, $Q_{dis}$ はそれぞれ温度に対する潜熱加熱項, 放射 加熱項, 散逸加熱項である. また $M_{cond}$, $M_{fall}$ はそれぞれ単位時間単位体積当たりの凝結量, 単 位時間単位体積当たりの雲粒の落下量である. $D_u$, $D_w$, $D_{\theta}$, $D_{\rho_s}$ の具体的な表式については 2.2 節で述べる. また $M_{cond}$, $M_{fall}$ の具体的な表式については 2.3 節で述べる. 熱力学の式 \Deqref{AC} の中の非断熱加熱項 $Q_{cond}$, $Q_{rad}$, $Q_{dis}$ について述べる. 潜熱加熱項については \begin{equation} Q_{cond} = \frac{L M_{cond}}{c_p \overline{\rho}} \Deqlab{AF} \end{equation} と表す. ここで $L$ は単位質量当たりの CO$_2$ の潜熱であり, 本モデルでは一定とみ なして $L \simeq 5.86 \times 10^5$ [J kg$^{-1}$] とする. $Q_{dis}$ の具体的な表式については 2.2 節で述べる. $Q_{rad}$ は陽に時間発展を解くことはせずに, 実験設定に応じて高度と時間に 依存するパラメータとして与える. 圧力方程式 \Deqref{AD} について注意すべき点をいくつか述べる. 2.2 節の冒頭で大気主成分の凝結の効果を考慮すると述べたが, 具体的には \Deqref{AD} の右辺第 2 項及び第 3 項に相当する. \Deqref{AD} の右辺第 2 項は凝結に伴う熱膨張の効果を表している. 熱膨張項は潜熱による効果のみを考慮して \begin{equation} \frac{\overline{c}^2}{c_p \overline{\theta}^2} \DD{\theta^{\prime}}{t} = \frac{\overline{c}^2}{c_p \overline{\theta}^2} \frac{\overline{\theta}}{\overline{T}} Q_{cond} \Deqlab{AG} \end{equation} と表現する. \Deqref{AD} の右辺第 3 項は凝結による大気質量の減少に伴う圧力低下の効果 を表している. 地球大気の場合, 凝結する H$_2$O が微量成分であることから, \Deqref{AD} の 右辺第 2 項及び第 3 項を無視して計算を行う場合が多い. \Deqref{AG} を用いると, \Deqref{AD} は以下のように書き換えられる. \begin{equation} \DP{\Pi^{\prime}}{t} = - \frac{\overline{c}^2}{c_p \overline{\rho} \overline{\theta}^2} \left[ \DP{ (\overline{\rho} \overline{\theta} u) }{x} + \DP{ (\overline{\rho} \overline{\theta} w) }{z} \right] + \frac{\overline{c}^2}{c_p \overline{\rho} \overline{\theta}} \left( \frac{L}{c_p \overline{T}} - 1 \right) M_{cond}. \Deqlab{AH} \end{equation} 火星極域の場合 $O(\overline{T}) \sim 150$ [K] であるので, $c_p \simeq 734.1$ [J K$^{-1}$ kg$^{-1}$], $L \simeq 5.86 \times 10^5$ [J kg$^{-1}$] より $L / (c_p \overline{T}) > 1 $ が成り立つ. 即ち火星極域で CO$_2$ の凝結が生じた場合, 圧力は増加することになる. 以下, 境界条件について述べる. 水平方向の境界では周期境界条件を適用し, 鉛直方向の境界では応力無し条件を 適用する. 計算領域の水平方向の長さを $L$ とすると, 任意の物理量 $\phi$ に対して課 されるべき水平方向の境界条件は \begin{equation} \phi(x,z,t) = \phi(x+L,z,t) \Deqlab{BoundaryA} \end{equation} と表される. また計算領域の鉛直方向の長さを $D$, $w$ 以外の予報変数を $\psi$ とすると, 鉛直方向の境界条件は \begin{eqnarray} w &=& 0, \quad \quad (z = 0, D) \Deqlab{boundaryB} \\ \DP{\psi}{z} &=& 0, \quad \quad (z = 0, D) \Deqlab{boundaryC} \end{eqnarray} と表される. \section{乱流パラメタリゼーション} 一般に大気の運動を記述する場合, 乱流運動(小規模運動)も考慮する必要がある. しかし数値モデルの場合有限個の格子点を配置して計算を行う為, 格子点間隔よ りも小さなスケール(サブグリッドスケール)の現象を取り扱うことが出来ない. そこでサブグリッドスケールの乱流運動に関する物理量を格子点上の物理量と結 びつけて表現する(乱流パラメタリゼーション). \subsection{乱流拡散} 本モデルでは Klemp and Wilhelmson(1978), 坪木・榊原(2001)に従って 乱流拡散項, 乱流拡散係数, 散逸加熱項を表現するのに 1.5 次のクロージャモ デルを用いる. 1.5 次のクロージャモデルを用いると $D_u$, $D_w$, $D_\theta$, $D_{\rho_s}$ は以下のように表現される. \begin{eqnarray} D_u &=& - \DP{}{x} \left[ - \frac{2}{3} K_m \DP{u}{x} + \frac{2}{3} E \right] - \DP{}{x} \left[ - K_m \left( \DP{u}{x} + \DP{w}{z} \right) \right] , \Deqlab{TurbA} \\ D_w &=& - \DP{}{x} \left[ - K_m \left( \DP{w}{x} + \DP{u}{z} \right) \right] - \DP{}{z} \left[ - \frac{2}{3} K_m \DP{w}{x} + \frac{2}{3} E \right] , \Deqlab{TurbB} \\ D_{\theta} &=& - \DP{}{x} \left[ K_h \DP{\theta^{\prime}}{x} \right] - \DP{}{z} \left[ K_h \DP{(\theta^{\prime} + \overline{\theta})}{z} \right], \Deqlab{TurbC} \\ D_{\rho_s} &=& - \DP{}{x} \left[ K_h \DP{\rho_s}{x} \right] - \DP{}{z} \left[ K_h \DP{\rho_s}{z} \right] \Deqlab{TurbD} \end{eqnarray} 但し $K_m$, $K_h$ はそれぞれ運動量, 熱に対する乱流拡散係数であり, $E$ は サブグリッドスケールの乱流運動エネルギーである. 次元解析に基づき, $E$ を $K_m$ で表現すると \begin{equation} E = \left( \frac{K_m}{C_m l} \right)^2 \Deqlab{TurbE} \end{equation} となる. ここで $l$ は混合距離であり, 水平方向の格子間隔 $\Delta x$ 及び鉛直方向 の格子間隔 $\Delta z$ を用いて $l = (\Delta x \Delta z)^{1/2}$ と表現す る. また $C_m$ は無次元定数である. $K_h$ については Deardorff(1972) に従って \begin{equation} K_h = 3 K_m \Deqlab{TurbF} \end{equation} と与える. 1.5 次のクロージャモデルでは以下に示す $K_m$ の時間発展を解き, 各乱流拡 散項を定める. \begin{eqnarray} \DP{K_m}{t} &=& - \left( u \DP{K_m}{x} + w \DP{K_m}{z} \right) - \frac{3 g {C_m}^2 l^2}{2 \overline{\theta}} \DP{\theta^{\prime}}{z} \notag \\ & & + {C_m}^2 l^2 \left[ \left( \DP{u}{x} \right)^2 + \left( \DP{w}{z} \right)^2 \right] + \frac{{C_m}^2 l^2}{2} \left( \DP{u}{z} + \DP{w}{x} \right)^2 \notag \\ & & - \frac{K_m}{3} \left( \DP{u}{x} + \DP{w}{z} \right) + \Dinv{2} \left( \DP[2]{{K_m}^2}{x} + \DP[2]{{K_m}^2}{z} \right) \notag \\ & & + \left( \DP{K_m}{x} \right)^2 + \left( \DP{K_m}{z} \right)^2 - \frac{C_{\varepsilon}}{2 C_m l^2} {K_m}^2. \Deqlab{TurbG} \end{eqnarray} ここで $C_{\varepsilon}$ は無次元定数である. 無次元定数の値については Deardorff(1975) に従って $C_m = C_{\varepsilon} = 0.2$ とする. 散逸加熱項 $Q_{dis}$ については \begin{equation} Q_{dis} = \frac{C_{\varepsilon}}{c_p l} E^{3/2} = \frac{C_{\varepsilon}}{c_p l} \left( \frac{K_m}{C_m l} \right)^3 \Deqlab{TurbH} \end{equation} と与える. \subsection{地表面フラックス} 地表面からの熱フラックス $F_h$ については簡略化されたバルク法(中島, 1994)を用いて \begin{equation} F_h = c_D V_{lowest} (T_{gnd} - T_{lowest}) / \Delta z \Deqlab{TurbI} \end{equation} と表現する. ここで $c_D$, $V_{lowest}$, $T_{gnd}$, $T_{lowest}$ はそれぞれバルク係数, 気層下端での水平風, 地表面温度, 気層下端での温度を表す. バルク係数は一定値であり, $c_D = 1.5 \times 10^{-3}$ とする. 気層下端での水平風は \begin{equation} V_{sfc} = \sqrt{u_{lowest}^2 + u_0^2} \Deqlab{TurbJ} \end{equation} と与える. ここで $u_{lowest}$ は予報して得られた気層下端の水平風であり, $u_0$ は定 数である. $u_0$ は $u_{lowest}$ が非常に弱くなった場合でもある程度の大きさの地表面 フラックスが保証されるようにする為の定数であり, $u_0 = 3.0$ [m/s] とする. $T_{gnd}$ は陽に解かずに, パラメータとして一定値を与える. \section{雲物理パラメタリゼーション} 2.3 節では単位時間単位体積当たりの凝結量 $M_{cond}$ 及び単位時間単位体積 当たりの雲粒の落下量 $M_{fall}$ の表し方について述べる. 地球大気の雲対流モデルでは, 過飽和状態は瞬時に解消されるものとみなし, 飽 和状態となるように凝結量を調整している(例えば Soong and Ogura(1973) ). しかし第 1 章で述べたように火星大気の場合, 過飽和状態が実現しうると考え られている. そこで本モデルでは雲粒が有限の時間をかけて成長するものと仮定 する. 火星大気での温度・圧力環境下では, CO$_2$ は気相或いは固相のいずれかの状態で 存在する. 従って本モデルで考えるべき雲物理過程は CO$_2$ 大気と CO$_2$ 氷雲の変換の みである. 地球大気の場合, 氷晶核の形成のメカニズムは大気中に存在する微粒子の影響を 受ける不均一核形成と, その影響を受けない均一核形成の 2 種類に分類される ことが知られている(例えば水野(2000)). 地球上の雲や霧などの大部分は不均一核形成によって生じると考えられている. 火星大気中にはダスト粒子や H$_2$O 氷でコーティングされたダストが浮遊して おり, 不均一核形成が生じると考えられる(Gooding, 1986). そこで本モデルでは核形成として不均一核形成のみを考える. 雲粒の成長過程は拡散成長と併合成長の 2 種類に分類される. 拡散成長とは気相中の分子が拡散によって雲粒表面に移動・凝結することによっ て雲粒が成長することを指す. 併合成長とは雲粒同士が衝突し, 合体することによって雲粒が成長することを指 す. Tobie et al.(2003) によると, 凝結核となるダストの平均半径が約 0.5 [$\mu$m] であり(Ockert-Bell et al., 1997), 可視光に対する典型的な光学的厚さがおお よそ 0.1 となる火星大気の場合, 半径 50 [$\mu$m] 程度の雲粒が併合成長で出来 るのに要する時間は拡散成長で出来るのに要する時間よりも十分長い. 北守(2006) は流れの存在しない火星大気での拡散成長についての数値計算を 行ない, 拡散成長の時間スケールが 1 -- 20 [s] 程度であること, 雲粒半径が 5 -- 25 [$\mu$m] となることを見出した. 雲粒半径が小さいほど拡散成長は併合成長に比べて優勢となるので, 火星大気環 境では拡散成長が優勢となると考えられる. そこで本モデルでは併合成長による雲粒成長は無視できると仮定し, 拡散成長の みを考える. 以下, 単位時間体積当たりの凝結量 $M_{cond}$ の取り扱いについて述べる. 本モデルでは単位質量の気相に含まれる凝結核の個数及び半径は空間的・時間的 に一様と仮定する. また雲粒の半径は各格子内において空間的に一定であると仮定する. 更に雲粒は球形の凝結核を核として形成され, 雲粒自身も球形となると仮定する. このとき \begin{equation} \frac{4}{3} \rho_I \pi ( r_d^3 - r_{aero}^3) N = \rho_s \Deqlab{CloudA} \end{equation} という関係式が成り立つ. ここで $\rho_I$ は CO$_2$ 氷の密度, $r_d$ は雲粒半径, $r_{aero}$ は凝結 核の半径, $N$ は単位体積当たりの凝結核の数密度である. 本モデルでは $\rho_I = 1.565 \times 10^3$ [kg/m$^3$] と与え, $r_{aero}$, $N$ は実験に応じて与える. \Dfigref{radi-dens} に雲密度の関数として雲粒半径をプロットしたものを示す. 雲粒の雲粒が拡散によって成長する場合の単位時間単位体積当たりの凝結量 $M_{cond}$ は以下のように表される \footnote{ \Deqref{CloudB} の導出については北守(2006)の付録 A を参照されたい. }. \begin{equation} M_{cond} = \frac{4 \pi r_d N}{R_h + R_m} (S - 1). \Deqlab{CloudB} \end{equation} ここで $R_h$, $R_m$, $S$ はそれぞれ熱輸送に関する定数, 質量輸送に関する 定数, 飽和比であり \begin{eqnarray} R_h &=& \frac{L^2}{k R T^2} = \frac{L^2}{k R \theta^2 \Pi^2}, \Deqlab{CloudC} \\ R_m &=& \frac{R T}{D p_{\ast}} = \frac{R \theta \Pi}{D p_{\ast}}, \Deqlab{CloudD} \\ S &=& \frac{p}{p_{\ast}} = \frac{{p_0}^{\frac{c_p}{R}} \Pi}{p_{\ast}} \Deqlab{CloudE} \end{eqnarray} と表される. 但し $k$, $D$, $p_{\ast}$ はそれぞれ熱拡散係数, 分子拡散係数, CO$_2$ の 飽和蒸気圧である. 本モデルでは Tobie et al. (2003) 同様に主成分凝結系では $R_h \gg R_m$ と して \begin{equation} M_{cond} = \frac{4 \pi r_d N}{R_h} (S - 1). \Deqlab{CloudF} \end{equation} と表す. CO$_2$ の飽和蒸気圧については半経験式である Antoine の式 \begin{equation} \ln p_{\ast} = A - \frac{B}{T - C} \Deqlab{CloudG} \end{equation} を用いて定める(Antoine, 1888). ここで $A$, $B$, $C$ は実験により定まる係数であり, CO$_2$ の場合 $A = 27.4$, $B = 3103$, $C = -0.16$ である(化学工学会, 1999). 火星大気環境における凝結を想定すると $O(T) \sim 150 $[K] であるので, $T \gg C$ と近似して \begin{equation} \ln p_{\ast} = A - \frac{B}{T} \Deqlab{CloudH} \end{equation} とする. 以下, 単位時間体積当たりの雲粒落下量 $M_{fall}$ の取り扱いについて述べる. $M_{fall}$ は Kessler(1969) と同様に, 雲粒の終端速度 $V_{term}$ での移流 として表現する. 即ち \begin{equation} M_{fall} = \DP{}{z} \left( \rho_s V_{term} \right) \Deqlab{CloudI} \end{equation} と表す. 終端速度 $V_{term}$ については球形粒子に関する Stokes 則を適用して \begin{equation} V_{term} = C_{sc} \frac{2 r_d^2 g \rho_I}{9 \eta} \Deqlab{CloudJ} \end{equation} と表す. ここで $C_{sc}$ は微小な粒子における Stokes 則からのずれを補正する係数 (Cunningham 補正係数)であり, \begin{equation} C_{sc} = 1 + 1.255 \frac{\lambda}{r_d} \Deqlab{CloudK} \end{equation} と表される(Cunningham, 1910). $\lambda$ は CO$_2$ の平均自由行程であり, \begin{equation} \lambda = \frac{k_B T}{\sqrt{2} \pi \sigma^2 p} \Deqlab{CloudL} \end{equation} と表される. $k_B$ は Boltzmann 定数, $\sigma$ は CO$_2$ 分子の直径であり, $k_B = 1.38 \times 10^{-23}$ [m$^2$ kg s$^{-2}$ K$^{-1}$], $\sigma = 3.3 \times 10^{-10}$ [m] である(Golden and Sircar, 1994). $\eta$ は粘性係数であり, Sutherland の公式 \begin{equation} \eta = \eta_{ref} \left( \frac{T_{ref} + C_{CO_2}}{T + C_{CO_2}} \right) \left( \frac{T}{T_{ref}} \right)^{3/2} \Deqlab{CloudM} \end{equation} で表現する(Sutherland, 1893). $\eta_{ref}$, $T_{ref}$, $C_{CO_2}$ はそれぞれ粘性係数の基準値, 温度の基 準値, CO$_2$ に関する Sutherland 定数であり, $\eta_{ref} = 1.47 \times 10^{-5}$ [Pa $\cdot$ s], $T_{ref} = 293$ [K], $C_{CO_2} = 240$ [K] と与え る(理科年表, 2004) \footnote{ Sutherland 係数は実験によってのみ定められるが, 公式の理論的導出は可能で ある. Sutherland の公式の理論的導出については付録 B を参照されたい. }. \Dfigref{Vterm} に火星極冠周縁の凝結高度付近の特徴的な温度圧力条件 ($T=150$[K], $p=400$[Pa])での雲粒の終端速度を示す. % % % \begin{figure} \begin{center} % \includegraphics[width=6.0cm,angle=270]{radius-dens.ps} \Depsf[][][width=6.0cm,angle=270]{ps/radius-dens.ps} \caption{雲粒半径を雲密度の関数としてプロットしたもの. 火星極冠周縁の凝結高度付近の特徴的な温度圧力条件 ($T=150$[K], $p=400$[Pa])を想定し, また凝結核数濃度を $5.0 \times 10^{8}$[kg$^{-1}$] とした. \Dfiglab{radi-dens} } \end{center} \end{figure} % % % \begin{figure} \begin{center} % \includegraphics[width=6.0cm,angle=270]{Vterm.ps} \Depsf[][][width=6.0cm,angle=270]{ps/Vterm.ps} \caption{火星極冠周縁の凝結高度付近の特徴的な温度圧力条件 ($T=150$[K], $p=400$[Pa])での雲粒の終端速度. \Dfiglab{Vterm} } \end{center} \end{figure} % % % \end{document}