Subsections
Kessler(1969) に基づく雲微物理パラメタリゼーションの, 終端速度 806#806, 雲
水の衝突併合による雨水比湿の変化率 136#136, 平均終端速度 172#172,
蒸発による雨水比湿の変化率 134#134 について解説する
.
Newton の抵抗法則より球体の抵抗力 807#807 は
と表される.
ここで 809#809 は抵力係数であり, 一般にレイノルズ数の関数である.
雨滴の落下のようにレイノルズ数が大きい現象の場合, レイノルズ数の定義によ
り, 粘性力は流れ場にほとんど寄与しなくなる.
このとき 809#809 はレイノルズ数に依存しない定数となる.
抵抗力と重力の釣合いを考えると
となる.
shuutan2 を 806#806 について解くと,
となる.
Kessler(1969) では 812#812[kg/m193#193], 813#813[m/s230#230],
814#814 [kg/m193#193], 815#815 として
としている
.
但し 817#817 は地表面での大気密度である.
他の惑星大気においても 815#815 であるとみなすと,
が得られる.
雲水の衝突併合による雨水混合比の変化率 136#136 は, 直径 205#205 の単一の
雨粒の衝突併合による質量変化率
819#819 と 205#205 から 820#820 の
範囲の直径を持つ雨粒の数 821#821 を用いて
と表される.
819#819 は,
と表される. ここで 806#806 は雨粒の落下速度, 824#824 は雨粒と衝突した雲粒
のうち雨粒に併合される割合を表す係数(捕捉係数)である.
雨粒のサイズ分布関数と雨粒の落下速度 806#806 を以下のように仮定する.
825#825 |
19#19 |
826#826 |
(249) |
827#827 |
19#19 |
828#828 |
(250) |
ここで
829#829 はパラメータである. 式マーシャル・パ
ルマー型分布関数2の分布は一般にマーシャル・パルマー型分布 (Marshall
and Palmer, 1948) と呼ばれる. Kessler (1969) では 830#830 とする.
これを式衝突併合による雨水混合比の変化率に代入すると,
831#831 |
19#19 |
832#832 |
(251) |
|
19#19 |
833#833 |
|
|
19#19 |
834#834 |
(252) |
を得る. ここで 824#824 は 205#205 によらないと仮定した.
Kessler (1969) では 835#835 とする.
雨粒のサイズ分布曲線の傾きを表すパラメータ 224#224 は, 以下の式を用
いて雨水比湿 300#300 で置き換える.
836#836 |
19#19 |
837#837 |
|
|
19#19 |
838#838 |
|
|
19#19 |
839#839 |
|
|
19#19 |
840#840 |
(253) |
ここで 841#841 は液相の密度である.
これを 224#224 について解き, 式CL_cr 項に代入すると,
831#831 |
19#19 |
842#842 |
|
|
19#19 |
843#843 |
|
|
19#19 |
844#844 |
(254) |
となる. 最後の式変形では,
845#845 を代入した.
平均終端速度 172#172 は雨滴の鉛直フラックス 846#846, 雨滴密度 15#15
により
と表される.
15#15, 846#846 はそれぞれ以下のように表される.
848#848 |
19#19 |
849#849 |
(256) |
850#850 |
19#19 |
851#851 |
(257) |
ここで 852#852 は直径 205#205 の雨滴の質量であり,
と書ける.
マーシャル・パルマー型分布関数1, マーシャル・パルマー型
分布関数2, TermVel4 を TermVel2, TermVel3 に
適用すると,
848#848 |
19#19 |
854#854 |
|
|
19#19 |
855#855 |
(259) |
850#850 |
19#19 |
856#856 |
|
|
19#19 |
857#857 |
(260) |
となる.
TermVel5 を TermVel6 に代入して 224#224 を消去すると,
850#850 |
19#19 |
858#858 |
|
|
19#19 |
859#859 |
(261) |
となる.
TermVel1 に TermVel7 を代入すると
860#860 |
19#19 |
861#861 |
|
|
521#521 |
862#862 |
(262) |
が得られる.
蒸発による雨水混合比の変化率 134#134 は, 式衝突併合による雨
水混合比の変化率 と同様に
と表される. ここで
864#864 は直径 205#205 の単一の雨粒の蒸発によ
る質量変化率である.
雨水の蒸発は雨粒の表面からの水蒸気の拡散によって律速されると仮定する.
雨粒周囲の水蒸気フラックスを 865#865 とすると, 雨粒の質量の変化率は
と表される. ここで 867#867 は雨粒中心からの距離, 868#868 は雨粒の半径で,
865#865 は
869#869
と表される. 870#870 は水蒸気の密度, 871#871 は水蒸気の拡散係数であ
る. 雨粒の周囲では水蒸気フラックスの収束発散はないと仮定すると,
872#872
が成り立つ. これを積分し
873#873
境界条件 874#874 で
875#875, 876#876 で
877#877 を適用すると,
878#878
これより, 雨粒表面での拡散による水蒸気フラックスは
879#879 |
19#19 |
880#880 |
|
|
19#19 |
881#881 |
(265) |
よって,
と表される. 雨粒が落下しながら蒸発する場合には, 871#871 に補正項のついた
が用いられる. ここで 865#865 は換気因子, 884#884 は雨粒表面でのクヌーセン層の
厚さである
.
Kessler (1969) では, 蒸発による雨粒の成長方程式の右辺の項を
以下のように近似する.
885#885
このとき蒸発による雨粒の成長方程式は
となる. これを式蒸発による雨水混合比の増加率に代入し,
雨粒のサイズ分布としてマーシャル・パルマー型分布関数1を
仮定すると,
887#887 |
19#19 |
888#888 |
|
|
19#19 |
889#889 |
|
|
19#19 |
890#890 |
|
|
19#19 |
891#891 |
|
|
19#19 |
892#892 |
|
|
19#19 |
893#893 |
(269) |
最後の式変形を行う際にはλの式式の関係を用いて 224#224 を消去し,
894#894, 830#830 とした
.
Footnotes
- ... について解説する
![[*]](footnote.png)
- 本章の内容は Ogura and Takahashi (1971), 浅井 (1983) を参考に
した.
- ...
としている
![[*]](footnote.png)
-
Kessler(1969) では 809#809 をどのように決めたのかについては書かれていない.
Gunn and Kinzer(1949) によると, レイノルズ数が 3000 程度である雨粒の
809#809 の値は 0.66 となるので, Kessler(1969) は系の特徴的なレイノルズ数が
3000 程度であると想定して 809#809 の値を決めたのかも知れない.
- ...
厚さである
![[*]](footnote.png)
-
蒸発による雨粒の成長方程式 は Kinzer and Gunn(1951) で導出され
ている.
導出方法については要確認である.
- ... とした
![[*]](footnote.png)
- Kessler (1969) では最終的には
895#895
としている.
Yamashita Tatsuya
2010-04-14