next up previous
: 4 ダストの輸送モデル : Two dimensional anelastic model : 2 大気モデル


3 乱流モデル

3.1 乱流パラメタリゼーション

移流項 $[\mbox{DKADV}]_{i,j}^{n}$ は熱力学の式の場合と同様に 4 次中央差 分, その他の項は 2 次中央差分で離散化する. 摩擦項 $[\mbox{DKDIF}]_{i,j}^{N}$, $[\mbox{DKNLD}]_{i,j}^{N}$ の時間積分に対し ては常に前進差分を用いる. $[\mbox{DKADV}]_{i,j}^{n}$, $[\mbox{DKDIF}]_{i,j}^{N}$ の表現はそれぞれ([*]), (24)式と同様なので, ここでは省略する.

$\displaystyle \varepsilon _{i,j}^{n+1}$ $\textstyle =$ $\displaystyle \varepsilon _{i,j}^{N} + dt \left\{
[\mbox{DKADV}]_{i,j}^{n} +
[\mbox{DKDIF}]_{i,j}^{N} +
[\mbox{DKNLD}]_{i,j}^{N} \right.$  
    $\displaystyle + \left.
[\mbox{DKBP}]_{i,j}^{n} +
[\mbox{DKSP}]_{i,j}^{n} -
\frac{C_{\epsilon}}{l}(\varepsilon _{i,j}^{N})^{\frac{3}{2}} \right\}$ (41)


$\displaystyle \mbox{DKNLD}_{i,j}^{n}$ $\textstyle =$ $\displaystyle \frac{1}{(\Delta x)^{2}}
\left[
K_{NLD,i+\frac{1}{2},j}^{n}(\vare...
...{i,j}^{n} -\varepsilon_{i-1,j}^{n})
\right] + \frac{1}{\rho _{0,j}\Delta z_{j}}$  
    $\displaystyle \left[ \rho _{0,j+\frac{1}{2}}
K_{NLD,i,j+\frac{1}{2}}^{n}
\frac{...
...epsilon_{i,j}^{n}-\varepsilon_{*i,j-1}^{n})}{\Delta z_{j-\frac{1}{2}}}
\right],$ (42)

\begin{eqnarray*}
K_{NLD,i+\frac{1}{2},j}^{n}
&=&
\mbox{MIN}
\left[ K_{NLD...
... /2000, \\
K_{NLD,max} &=& 0.2\frac{(\Delta x)^{2}}{\Delta t}.
\end{eqnarray*}


$\displaystyle \mbox{DKBP}_{i,j}^{n}$ $\textstyle =$ $\displaystyle -\frac{g}{\Theta _{0,j}}K_{i,j}^{n}
\frac{1}{\Delta z_{j}}\left[
...
...+\frac{1}{2}}) -
(\theta _{i,j-\frac{1}{2}}+\Theta _{0,j-\frac{1}{2}}) \right],$ (43)
$\displaystyle \mbox{DKSP}_{i,j}^{n}$ $\textstyle =$ $\displaystyle 2K_{i,j}^{n}\left[
\left(
\frac{u_{i+\frac{1}{2},j}^{n}-u_{i-\fra...
...i,j+\frac{1}{2}}^{n}-w_{i,j-\frac{1}{2}}^{n}}{\Delta z_{j}}
\right)^{2}
\right]$  
    $\displaystyle + \frac{2}{3}\varepsilon _{i,j}^{n} \left[
\frac{u_{i,j+\frac{1}{...
..._{j}} +
\frac{w_{i+\frac{1}{2},j}^{n}-w_{i-\frac{1}{2},j}^{n}}{\Delta x}\right]$  
    $\displaystyle + K_{i,j}^{n}
\left[
\frac{u_{i,j+\frac{1}{2}}^{n}-u_{i,j-\frac{1...
...} +
\frac{w_{i+\frac{1}{2},j}^{n}-w_{i-\frac{1}{2},j}^{n}}{\Delta x}\right]^{2}$ (44)


\begin{displaymath}
u_{i,j+\frac{1}{2}}^{n} =
0.5\left(u_{i+\frac{1}{2},j+\fra...
...rac{1}{2}}^{n} +
w_{i+\frac{1}{2},j-\frac{1}{2}}^{n}\right).
\end{displaymath}

3.2 地表フラックスパラメタリゼーション

地表からの運動量と熱のフラックスは以下のように離散化して計算する.

$\displaystyle F_{u,i}$ $\textstyle =$ $\displaystyle - \rho _{0}C_{D,i}\vert u_{i,\frac{1}{2}}\vert u_{i,\frac{1}{2}},$ (45)
$\displaystyle F_{\theta,i}$ $\textstyle =$ $\displaystyle \rho _{0}C_{D,i}\vert u_{i,\frac{1}{2}}\vert(T_{sfc,i}-T_{i,1}).$ (46)

ここで,


\begin{displaymath}
C_{D,i} = \left\{
\begin{array}{lcl}
C_{Dn}\left( 1 -
\f...
...B,i})^{2}}& for & \mbox{Ri}_{B,i} \geq 0,
\end{array}\right.
\end{displaymath} (47)

である. バルクリチャードソン数は以下のように評価する.

\begin{displaymath}
\mbox{Ri}_{B,i} \equiv \frac{gz_{1}(\Theta _{sfc,i}-\Theta _
{i,1})}{\overline{\Theta }_{0,1}u_{i,\frac{1}{2}}}.
\end{displaymath} (48)


next up previous
: 4 ダストの輸送モデル : Two dimensional anelastic model : 2 大気モデル
Odaka Masatsugu 平成19年4月26日