Subsections

C. 雲微物理過程

Kessler(1969) に基づく雲微物理パラメタリゼーションの, 終端速度 806#806, 雲 水の衝突併合による雨水比湿の変化率 136#136, 平均終端速度 172#172, 蒸発による雨水比湿の変化率 134#134 について解説する [*].

1 雨粒の終端速度

Newton の抵抗法則より球体の抵抗力 807#807

808#808     (242)

と表される. ここで 809#809 は抵力係数であり, 一般にレイノルズ数の関数である. 雨滴の落下のようにレイノルズ数が大きい現象の場合, レイノルズ数の定義によ り, 粘性力は流れ場にほとんど寄与しなくなる. このとき 809#809 はレイノルズ数に依存しない定数となる. 抵抗力と重力の釣合いを考えると
810#810     (243)

となる. shuutan2 を 806#806 について解くと,
811#811     (244)

となる. Kessler(1969) では 812#812[kg/m193#193], 813#813[m/s230#230], 814#814 [kg/m193#193], 815#815 として
816#816     (245)

としている[*]. 但し 817#817 は地表面での大気密度である. 他の惑星大気においても 815#815 であるとみなすと,
818#818 (246)

が得られる.

2 雲水の衝突併合

雲水の衝突併合による雨水混合比の変化率 136#136 は, 直径 205#205 の単一の 雨粒の衝突併合による質量変化率 819#819205#205 から 820#820 の 範囲の直径を持つ雨粒の数 821#821 を用いて

822#822 (247)

と表される. 819#819 は,
823#823 (248)

と表される. ここで 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#824205#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 を代入した.

3 平均終端速度

平均終端速度 172#172 は雨滴の鉛直フラックス 846#846, 雨滴密度 15#15 により

847#847     (255)

と表される. 15#15, 846#846 はそれぞれ以下のように表される.
848#848 19#19 849#849 (256)
850#850 19#19 851#851 (257)

ここで 852#852 は直径 205#205 の雨滴の質量であり,
853#853     (258)

と書ける. マーシャル・パルマー型分布関数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)

が得られる.

4 雨水の蒸発

蒸発による雨水混合比の変化率 134#134 は, 式衝突併合による雨 水混合比の変化率 と同様に

863#863 (263)

と表される. ここで 864#864 は直径 205#205 の単一の雨粒の蒸発によ る質量変化率である.

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

866#866 (264)

と表される. ここで 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)

よって,
882#882 (266)

と表される. 雨粒が落下しながら蒸発する場合には, 871#871 に補正項のついた
883#883 (267)

が用いられる. ここで 865#865 は換気因子, 884#884 は雨粒表面でのクヌーセン層の 厚さである[*].

Kessler (1969) では, 蒸発による雨粒の成長方程式の右辺の項を 以下のように近似する.

885#885


このとき蒸発による雨粒の成長方程式は
886#886 (268)

となる. これを式蒸発による雨水混合比の増加率に代入し, 雨粒のサイズ分布としてマーシャル・パルマー型分布関数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

... について解説する[*]
本章の内容は Ogura and Takahashi (1971), 浅井 (1983) を参考に した.
... としている[*]
Kessler(1969) では 809#809 をどのように決めたのかについては書かれていない. Gunn and Kinzer(1949) によると, レイノルズ数が 3000 程度である雨粒の 809#809 の値は 0.66 となるので, Kessler(1969) は系の特徴的なレイノルズ数が 3000 程度であると想定して 809#809 の値を決めたのかも知れない.
... 厚さである[*]
蒸発による雨粒の成長方程式 は Kinzer and Gunn(1951) で導出され ている. 導出方法については要確認である.
... とした[*]
Kessler (1969) では最終的には

895#895

としている.
Yamashita Tatsuya 2010-04-14