% 表題  中島健介対流モデル解説  乱流粘性拡散項の計算 (TSTEP)
%
% 履歴  1998/07/12  小高正嗣
%       1999/06/14  小高正嗣
%
\subsection{乱流粘性拡散項の計算}

運動方程式中の乱流粘性拡散項( {\tt DUVIS}, {\tt DVVIS}, {\tt DWVIS} ).
をサブルーチン {\tt CLVISC} を呼んで計算する

\subsubsection{CLVISC}

\begin{enumerate}
  \item 機能
    \begin{quote}
      運動方程式中の乱流粘性拡散項( {\tt DUVIS}, {\tt DVVIS}, {\tt
      DWVIS} )を計算する
    \end{quote}
  \item パラメータ
    \begin{quote}
      \begin{tabular}{lllp{10cm}}
        {\tt U}    & (実数) & {\tt ( -2:NXB, -2:NZB, NROT )} & 風速 x 成分 \\
        {\tt V}    & (実数) & {\tt ( -2:NXB, -2:NZB, NROT )} & 風速 y 成分 \\
        {\tt W}    & (実数) & {\tt ( -2:NXB, -2:NZB, NROT )} & 風速 z 成分 \\
        {\tt CD1} & (実数) & {\tt ( -2:NXB, -2:NZB )} & 乱流拡散係数(
        $(i+\frac{1}{2}), j+\frac{1}{2})$ 格子上) \\
        {\tt CD2} & (実数) & {\tt ( -2:NXB, -2:NZB )} & 乱流拡散係数(制限値) \\
        {\tt EDFM1} & (実数) & {\tt ( -2:NXB, -2:NZB )} & \\
        {\tt ICL} & (整数) & & 時間方向次元 \\
        {\tt DIV} & (実数) & {\tt ( -2:NXB, -2:NZB )} & 
        風の発散(スカラー量の格子点)\\  
        {\tt DUVIS} & (実数) & {\tt ( -2:NXB, -2:NZB )} & 
        風速 x 成分変化量(拡散)\\
        {\tt DVVIS} & (実数) & {\tt ( -2:NXB, -2:NZB )} & 
        風速 y 成分変化量(拡散)\\
        {\tt DWVIS} & (実数) & {\tt ( -2:NXB, -2:NZB )} & 
        風速 z 成分変化量(拡散)\\
        {\tt ZGRID} & (実数) & {\tt ( -2:NZB )} & 鉛直方向格子点座標 \\
        {\tt DENS0} & (実数) & {\tt ( -2:NZB )} & 基本場の密度 \\
        {\tt FDZ}   & (実数) & {\tt ( -2:NZB )} & z 方向格子間隔 \\
        {\tt FDZM}  & (実数) & {\tt ( -2:NZB )} & z 方向半整数格子間隔 \\
        {\tt CD00}  & (実数) & {\tt ( -2:NZB )} & 数値拡散 \\
        {\tt DX}    & (実数) & {\tt ( -2:NZB )} & x 方向格子間隔 \\
        {\tt NX  }  & (整数) & & 水平方向格子点数 \\
        {\tt NZ  }  & (整数) & &  鉛直方向格子点数 \\
        {\tt NXB }  & (整数) & & 水平方向次元数 {\tt (=nx+2)} \\
        {\tt NZB }  & (整数) & & 鉛直方向次元数 {\tt (=nz+2)} \\
        {\tt NROT}  & (整数) & & 時間方向次元数
      \end{tabular}
    \end{quote}
  \item 解説

    運動方程式中の乱流拡散項
    \begin{eqnarray*}
      && \DP{}{x}\left( K_{m}\DP{u}{x}\right) + 
         \frac{1}{\rho _{0}}\DP{}{z}\left( \rho _{0}K_{m}\DP{u}{z}\right), \\
      && \DP{}{x}\left( K_{m}\DP{v}{x}\right) + 
         \frac{1}{\rho _{0}}\DP{}{z}\left( \rho _{0}K_{m}\DP{v}{z}\right), \\
      && \DP{}{x}\left( K_{m}\DP{w}{x}\right) + 
         \frac{1}{\rho _{0}}\DP{}{z}\left( \rho _{0}K_{m}\DP{w}{z}\right),
    \end{eqnarray*}
    を 2 次中央差分で計算する.

    \begin{enumerate}
      \item {\tt ICL} には {\tt IT2} が渡される.
      \item 数値拡散 {\tt CD00} を与える.
      \item {\tt DIV} を計算する.
      \item {\tt DUVIS}, {\tt DVVIS}, {\tt DWVIS} を計算する.
    \end{enumerate}

\end{enumerate}




