next up previous
: D. 変数リスト : 湿潤大気における 2 次元非静力学モデルの定式化 : B. 乱流パラメタリゼーション


C. 雲微物理過程

本モデルで用いている雲微物理パラメタリゼーション(Kessler, 1969)の, 雲 水の衝突併合による雨水混合比の変化率 $CL_{cr}$ と, 蒸発による雨水混合 比の変化率 $EV_{rv}$ について解説する C.1.

C.1 雲水の衝突併合

雲水の衝突併合による雨水混合比の変化率 $CL_{cr}$ は, 直径 $D$ の単一の 雨粒の衝突併合による質量変化率 $(dm(D)/dt)_{cr}$$D$ から $D+dD$ の 範囲の直径を持つ雨粒の数 $N_{D}$ を用いて

\begin{displaymath}
CL_{cr} = \frac{1}{\rho_{d}}\int _{0}^{\infty}\left(\DD{m}{t}\right)_{cr}
N_{D}\Dd D
\end{displaymath} (C.1)

と表される. $(dm(D)/dt)_{cr}$ は,
\begin{displaymath}
\left(\DD{m}{t}\right)_{cr} = \frac{\pi}{4}D^{2}VE\rho _{d}q_{c}
\end{displaymath} (C.2)

と表される. ここで $V$ は雨粒の落下速度, $E$ は雨粒と衝突した雲粒 のうち雨粒に併合される割合を表す係数(捕捉係数)である.

Kessler (1969) では, 雨粒のサイズ分布関数と雨粒の落下速度 $V$ を以下の ように仮定する.

$\displaystyle N_{D}$ $\textstyle =$ $\displaystyle N_{0}\exp (-\lambda D),$ (C.3)
$\displaystyle V$ $\textstyle =$ $\displaystyle 130D^{0.5}.$ (C.4)

ここで $N_{0}, \lambda $ はパラメータである. 式([*])の分布は一般にマーシャル・パルマー型分布 (Marshall and Palmer, 1948) と呼ばれる. Kessler (1969) では $N_{0}=10^{7}$ とする. これを式(C.1)に代入すると,
$\displaystyle CL_{cr}$ $\textstyle =$ $\displaystyle \frac{130\pi}{4}E N_{0} q_{c}
\int_{0}^{\infty}D^{2.5}\exp(-\lambda D)\Dd D$ (C.5)
  $\textstyle =$ $\displaystyle 32.5 \pi E N_{0} q_{c} \frac{3.75}{\lambda ^{3}}\frac{1}{2}
\sqrt{\frac{\pi}{\lambda}}$  
  $\textstyle =$ $\displaystyle 60.9375 \pi^{3/2} E N_{0} q_{c} \lambda ^{-3.5}$ (C.6)

を得る. ここで $E$$D$ によらないと仮定した. Kessler (1969) では $E=1$ とする.

雨粒のサイズ分布曲線の傾きを表すパラメータ $\lambda$ は, 以下の式を用 いて雨水混合比 $q_{r}$ で置き換える.

$\displaystyle q_{r}$ $\textstyle =$ $\displaystyle \frac{1}{\rho_{d}}\int _{0}^{\infty} \rho_{w}
\frac{\pi}{6} D^{3}
N_{D}\Dd D$  
  $\textstyle =$ $\displaystyle \frac{\pi N_{0}\rho_{w}}{6\rho_{d}}\int _{0}^{\infty}
D^{3}\exp(-\lambda D)\Dd D$  
  $\textstyle =$ $\displaystyle \frac{\pi N_{0}\rho_{w}}{6\rho_{d}}\frac{6}{\lambda ^{4}}$  
  $\textstyle =$ $\displaystyle \frac{\pi N_{0}\rho_{w}}{\rho_{d}}\lambda ^{-4}.$ (C.7)

ここで $\rho_{w}$ は水の密度 ($10^{3}$ kg/m${}^{3}$) である. これを $\lambda$ について解き, 式(C.6)に代入すると,
$\displaystyle CL_{cr}$ $\textstyle =$ $\displaystyle 60.9375 \pi ^{\frac{5}{8}}\rho_{w}^{-\frac{7}{8}}
EN_{0}^{\frac{1}{8}}q_{c}(\rho _{d}q_{r})^{\frac{7}{8}}$  
  $\textstyle =$ $\displaystyle 0.295 \times 10^{-4}E N_{0}^{0.125}
q_{c}(\rho _{d}q_{r})^{0.875}$  
  $\textstyle =$ $\displaystyle 2.2 q_{c}(\rho _{d}q_{r})^{0.875}$ (C.8)

となる. 最後の式変形では, $N_{0}=10^{7}, E=1$ を代入した.

C.2 雨水の蒸発

蒸発による雨水混合比の変化率 $EV_{rv}$ は, 式([*]) と同様に

\begin{displaymath}
EV_{rv} = \frac{1}{\rho_{d}}\int _{0}^{\infty}\left(\DD{m}{t}\right)_{ev}
N_{D}\Dd D
\end{displaymath} (C.9)

と表される. ここで $(dm(D)/dt)_{ev}$ は直径 $D$ の単一の雨粒の蒸発によ る質量変化率である.

雨水の蒸発は雨粒の表面からの水蒸気の拡散によって律速されると仮定する. 雨粒周囲の水蒸気フラックスを $F$ とすると, 雨粒の質量の変化率は

\begin{displaymath}
\left(\DD{m}{t}\right)_{ev} = - 4\pi r_{d}^{2}F(r_{d})
\end{displaymath} (C.10)

と表される. ここで $r$ は雨粒中心からの距離, $r_{d}$ は雨粒の半径で, $F$

\begin{displaymath}
F = - K_{d}\DD{\rho_{v}}{r}
\end{displaymath}

と表される. $\rho_{v}$ は水蒸気の密度, $K_{d}$ は水蒸気の拡散係数であ る. 雨粒の周囲では水蒸気フラックスの収束発散はないと仮定すると,

\begin{displaymath}
\frac{1}{r^{2}}\DP{}{r}\left(r^{2}F\right) = 0
\end{displaymath}

が成り立つ. これを積分し

\begin{displaymath}
\rho _{v} = -\frac{C_{1}}{r} + C_{2}
\end{displaymath}

境界条件 $r=r_{d}$ $\rho_{v} = \rho_{v,s}$, $r=\infty$ $\rho_{v}=\rho_{v,\infty}$ を適用すると,

\begin{displaymath}
C_{1} = (\rho _{v,\infty}-\rho_{v,s})r_{d},
\quad C_{2} = \rho _{v,\infty}
\end{displaymath}

これより, 雨粒表面での拡散による水蒸気フラックスは
$\displaystyle F(r_{d})$ $\textstyle =$ $\displaystyle - K_{d}\left.\DD{\rho_{v}}{r}\right\vert _{r=r_{d}}$  
  $\textstyle =$ $\displaystyle K_{d}\frac{\rho_{v,s} - \rho _{v,\infty}}{r_{d}}$ (C.11)

よって,
\begin{displaymath}
\left(\DD{m}{t}\right)_{ev} = - 4\pi r_{d}K_{d}(\rho_{v,s}-\rho_{v,\infty})
\end{displaymath} (C.12)

と表される. 雨粒が落下しながら蒸発する場合には, $K_{d}$ に補正項のついた
\begin{displaymath}
\left(\DD{m}{t}\right)_{ev} = - 4\pi r_{d}\left(1+\frac{Fr}{s}\right)
K_{d}(\rho_{v,s}-\rho_{v,\infty})
\end{displaymath} (C.13)

が用いられる. ここで $F$ は換気因子, $s$ は雨粒表面でのクヌーセン層の 厚さであるC.2.

Kessler (1969) では, (C.13)の右辺の項を 以下のように近似する.

\begin{eqnarray*}
4\pi r_{d}\left(1+\frac{Fr}{s}\right) &\sim& 2.24\times 10^{3...
...} - \rho_{v,\infty}) &\sim& 10^{-5}(\rho_{v,s}-\rho_{v,\infty}).
\end{eqnarray*}

このとき(C.13)は
\begin{displaymath}
\left(\DD{m}{t}\right)_{ev} \sim
- 2.24\times 10^{-2}(\rho_{v,s}-\rho_{v,\infty})D^{1.6}
\end{displaymath} (C.14)

となる. これを式(C.9)に代入し, 雨粒のサイズ分布として(C.4)を 仮定すると,
$\displaystyle EV_{rv}$ $\textstyle =$ $\displaystyle - \frac{1}{\rho_{d}}2.24\times 10^{-2}(q_{v,s} -q_{v,\infty})\int _{0}^{\infty}
D^{1.6}N_{D}\Dd D,$  
  $\textstyle =$ $\displaystyle - 2.24\times 10^{-2}(q_{v,s} -q_{v,\infty})\int _{0}^{\infty}
D^{1.6}N_{0}\exp(-\lambda D)\Dd D$  
  $\textstyle =$ $\displaystyle - 2.24\times 10^{-2}(q_{v,s} -q_{v,\infty})
N_{0}\frac{\Gamma(2.6)}{\lambda ^{13/5}}$  
  $\textstyle =$ $\displaystyle - 2.24\times 10^{-2}
\Gamma(2.6) (\pi \rho_{w})^{-0.65}
N_{0}^{0.35}
(q_{v,s} -q_{v,\infty})
(\rho_{d}q_{r})^{0.65}$  
  $\textstyle =$ $\displaystyle - 1.7\times 10^{-4}N_{0}^{0.35}(q_{v,s} -q_{v,\infty})
(\rho_{d}q_{r})^{0.65}$  
  $\textstyle =$ $\displaystyle - 4.81\times 10^{-2}(q_{v,s} -q_{v,\infty})
(\rho_{d}q_{r})^{0.65}$ (C.15)

最後の式変形を行う際には(C.7)式の関係を用いて $\lambda$ を消去し, $\Gamma(2.6)=1.4296245$, $N_{0}=10^{7}$ とした C.3.



... について解説するC.1
本章の内容は浅井 (1983) の解説を参考にした.
... 厚さであるC.2
この式の導出は要確認.
... としたC.3
Kessler (1969) では最終的には

\begin{displaymath}
EV_{rv} = 4.85\times 10^{-2}(q_{v,s} -q_{v,\infty})(\rho_{d}q_{r})^{0.65}
\end{displaymath}

としている.

next up previous
: D. 変数リスト : 湿潤大気における 2 次元非静力学モデルの定式化 : B. 乱流パラメタリゼーション
Odaka Masatsugu 平成19年9月11日