%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% 表題  2 次元非静力学モデル -- 離散モデル  熱力学の式時間方向の離散化
%
% 履歴: 2005-01-25 杉山耕一朗
%       2005/04/13  小高正嗣
%       2005/08/25  小高正嗣: 散逸加熱項を追加
%       2005/08/26  小高正嗣: 基本場温位を乱流拡散を追加
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{熱力学の式}

熱力学の式を書き下すと以下のようになる. 
%
\begin{eqnarray}
\DP{\theta}{t} = F_\theta
\Deqlab{theta:theta}
\end{eqnarray}
%
\Deqref{theta:theta} 式を時間方向にリープフロッグ法を用いて離散化する. 
%
\begin{eqnarray}
 \theta_{i,k}^{t} = \theta_{i,k}^{t - \Delta t} + 2 F_{\theta,i,k}\Delta t  
\end{eqnarray}
% 
ここで, 
%
\begin{eqnarray}
 F_{\theta,i,k} &=&
  - \left[{\rm Adv}_{\theta}\right]_{i,k}^{t}
  - \left[{\rm Src}_{\theta}\right]_{i,k}^{t}
  + \left[D_{\theta} +  {\rm Diff}_{\theta} \right]_{i,k}^{t - \Delta t} 
    \nonumber \\
 &&  + \frac{\overline{\theta}_{i,k}}{\overline{T}_{i,k}}
	\left(Q_{rad,i,k}^{t-\Delta t} + Q_{dis,i,k}^{t-\Delta t}\right)
  \Deqlab{Ftheta}
\end{eqnarray} 
%
である. 移流を中心差分で安定して解くために, 数値粘性項 Diff を追加してあ
る. 


\Deqref{Ftheta} 式を書き下す. 移流項は, 
%
\begin{eqnarray}
  \left[{\rm Adv}_{\theta}\right]_{i,k}^{t}  &=& 
      \left\{ 
         u_{i(u),k}  \left( \DP{\theta}{x} \right)_{i(u),k} 
      \right\}_{i,k}^{t}
      + 
      \left\{ 
         w_{i,k(w)} \left( \DP{\theta}{z} \right)_{i,k(w)} 
      \right\}_{i,k}^{t}
\end{eqnarray}
%
であり, 生成項は, 
%
\begin{eqnarray}
  \left[{\rm Src}_{\theta}\right]_{i,k}^{t} =
      \left\{ 
         w_{i,k(w)} \left( \DP{\overline{\theta}}{z} \right)_{i,k(w)} 
      \right\}_{i,k}^{t}
\end{eqnarray}
%
であり, 粘性拡散項は, 
%
\begin{eqnarray}
\left[D_{\theta} \right]_{i,k}^{t - \Delta t} &=& 
   \left[ \DP{}{x} 
      \left\{
         \left( K_{h} \right)_{i(u),k}
         \left( \DP{\theta}{x} \right)_{i(u),k}
       \right\}
   \right]_{i,k}^{t - \Delta t}
\nonumber \\
&&
     + \left[ \DP{}{z}\left\{
       \left( K_{h} \right)_{i,k(w)}
       \left( \DP{(\theta + \overline{\theta})}{z} \right)_{i,k(w)}
     \right\} \right]_{i,k}^{t - \Delta t}
 \end{eqnarray}
%
である. 数値粘性項は,
%
 \begin{eqnarray}
  \left[ {\rm Diff}_{\theta} \right]_{i,k}^{t - \Delta t} &=& 
     \nu_{h} \left\{ \DP{}{x} \left(\DP{\theta}{x}\right)_{i(u),k} \right\}_{i,k}^{t - \Delta t}
   + \nu_{v} \left\{ \DP{}{z} \left(\DP{\theta}{z}\right)_{i,k(w)} \right\}_{i,k}^{t - \Delta t}
 \end{eqnarray}
%
である. $K_{h}$ は乱流エネルギーの時間発展方程式から計算する(詳細は後述). 
$\nu_{h}, \nu_{v}$ は \Deqref{nu} 式を利用する. 

散逸加熱項 $Q_{dis,i,k}^{t-\Delta t}$ は
\begin{equation}
  Q_{dis,i,k}^{t-\Delta t} = \frac{1}{c_{p}}\frac{C_{\varepsilon}}{l}
            \frac{(K_{m,i,k}^{t-\Delta t})^{3}}{(C_{m}l)^{3}}
\end{equation}
と与える. ここで $l=(\Delta x\Delta z)^{1/2}$ である.

$Q_{rad,i,k}$ は計算設定ごとに与える.