next up previous
: 表面フラックス : 力学過程 : 放射

鉛直混合

鉛直乱流混合は Mellor and Yamada (1974)のレベル II の方法を用いている.

Mellor and Yamada (1974)のこの方法は, 拡散係数 $K_M$, $K_T$, and $K_q$ は 安定度に依存し, $u$$v$$T$ の関数として 与えられる.

$\displaystyle K_M$ $\textstyle =$ $\displaystyle l^2 \frac{\Delta \vert\Dvect{v}\vert}{\Delta z} S_M,$ (34)
$\displaystyle K_T = K_q$ $\textstyle =$ $\displaystyle l^2 \frac{\Delta \vert\Dvect{v}\vert}{\Delta z} S_H.$ (35)

Blakadar(1962) にしたがって, 混合距離 $l$ は次式で評価される.
\begin{displaymath}
l = \frac{\kappa z}{1+\kappa z/l_0},
\end{displaymath} (36)

ここで $l_0=300$ m で $\kappa=0.4$ は Karman 定数である.

これらの表現に関して, $S_M$$S_H$ は Richardson 数の関数として与えられる.

$\displaystyle S_M$ $\textstyle =$ $\displaystyle B_1^{1/2} ( 1- R_f )^{1/2}
\widetilde{S_M}^{3/2},$ (37)
$\displaystyle S_H$ $\textstyle =$ $\displaystyle B_1^{1/2} ( 1- R_f )^{1/2}
\widetilde{S_M}^{1/2} \widetilde{S_H},$ (38)

ここで
\begin{displaymath}
\widetilde{S_H} = \frac{ \alpha_1-\alpha_2 R_f }{ 1-R_f },
\end{displaymath} (39)


\begin{displaymath}
\widetilde{S_M} = \frac{ \beta_1-\beta_2 R_f }
{ \beta_3-\beta_4 R_f }
\widetilde{S_H}.
\end{displaymath} (40)

$R_f$ Richardson 数フラックスで, 次のように与えられ

\begin{displaymath}
R_f = \frac{1}{2 \beta_2}
\left[ \beta_1 + \beta_4 R_{iB}
...
... \beta_4 R_{iB} )^2
- 4 \beta_2 \beta_3 R_{iB} }
\right] ,
\end{displaymath} (41)

ここで $R_{iB}$ はバルク Richardson 数で
$\displaystyle \alpha_1$ $\textstyle =$ $\displaystyle 3 A_2 \gamma_1,$ (42)
$\displaystyle \alpha_2$ $\textstyle =$ $\displaystyle 3 A_2 (\gamma_1+\gamma_2),$ (43)
$\displaystyle \beta_1$ $\textstyle =$ $\displaystyle A_1 B_1 ( \gamma_1 - C_1 ),$ (44)
$\displaystyle \beta_2$ $\textstyle =$ $\displaystyle A_1 [ B_1 ( \gamma_1 - C_1 ) + 6 A_1 + 3 A_2 ],$ (45)
$\displaystyle \beta_3$ $\textstyle =$ $\displaystyle A_2 B_1,$ (46)
$\displaystyle \beta_4$ $\textstyle =$ $\displaystyle A_2 [ B_1 ( \gamma_1 + \gamma_2 ) - 3 A_1 ],$ (47)


\begin{displaymath}
(A_1, B_1, A_2, B_2, C_1 ) = ( 0.92, 16.6, 0.74, 10.1, 0.08 ) ,
\end{displaymath} (48)


\begin{displaymath}
\gamma_1 = \frac{1}{3} - \frac{2 A_1}{B_1}\, , \, \, \,
\gamma_2 = \frac{B_2}{B_1} + 6\frac{A_1}{B_1},
\end{displaymath} (49)

である.

$R_{iB}$ は次で定義される.

\begin{displaymath}
R_{iB} = \frac{\displaystyle
\frac{q}{\theta_s} \frac{\De...
...\right)^2
+ \left( \frac{\Delta v}{\Delta z} \right)^2
}.
\end{displaymath} (50)

ここで, 鉛直高度 ${k-1/2}$ での $(\Delta \bullet)$ の値は $\bullet_{k}
- \bullet_{k-1}$によって評価されている.


next up previous
: 表面フラックス : 力学過程 : 放射
okuyama naonori 平成13年1月10日