next up previous
: 2. 参考文献 : 湿潤大気における 2 次元非静力学モデルの定式化 : 目次


1. 基礎方程式系

本数値モデルは水平・鉛直の 2 次元モデルである. 水平方向の座標変数を $x$, 鉛直方向の座標変数を $z$ と表し, 時間方向の変数は $t$ と表す.

1.1 運動方程式・圧力方程式・熱の式・比湿の保存式

力学的な枠組みは, 準圧縮方程式系(Klemp and Wilhelmson,1978)を用いる. この方程式系では, 予報変数を水平一様な基本場とそこからのずれに分離し, 方程式の線形化を行っている. 準圧縮方程式系の導出は付録 A に示す. 方程式中の変数は付録 D に示す.

以下に準圧縮方程式系の時間発展方程式を一覧する. 密度の式では乾燥成分と湿潤成分の分子量の差を考慮するが, 熱の式では考慮し ない. 比熱, 気体定数が空間的, 時間的に一様で, かつ凝結量が気相質量に比べて十分 少ないと仮定する.

運動方程式

$\displaystyle \DP{u^{\prime}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{u^{\prime}}{x} - w^{\prime} \DP{u^{\prime}}{z}
- \langle c_p \rangle \overline{\theta_v} \DP{\Pi^{\prime}}{x} +
D_{u^{\prime}}$ (1.1)
$\displaystyle \DP{w^{\prime}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{w^{\prime}}{x} - w^{\prime} \DP{w^{\prime}}{z}
- \langle c_p \rangle \overline{\theta_v} \DP{\Pi^{\prime}}{z}$  
    $\displaystyle + \frac{\theta^{\prime}}{\overline{\theta}} g
+ \frac{q_d^{\prime...
... \rangle /M_d +
\sum \overline{q_v} \langle M \rangle / M_v} g
+ D_{w^{\prime}}$ (1.2)

圧力方程式

$\displaystyle \DP{ \Pi^{'}}{t}$ $\textstyle =$ $\displaystyle - \frac{\overline{{c_{s}}^{2}}}{\langle c_p \rangle \bar{\rho}
\bar{\theta_v}^{2}}
\Ddiv \left(
\bar{\rho} \bar{\theta_v} \Dvect{u^{'}}
\right)$  
    $\displaystyle + \frac{ \overline{\langle c_s \rangle^2} }{\langle c_p \rangle
\...
...t\{ \Dinv{\bar{\Pi}} ( Q_{cond} + Q_{rad} + Q_{dis} )
+ D_{\theta^{'}} \right\}$  
    $\displaystyle + \Dinv{ \bar{q_d} \langle M \rangle / M_d
+ \sum \bar{q_v} \lang...
...- \frac{\bar{q_d}}{\bar{\rho}} \sum M_{fall} (\rho_r^{'}) +
D_{q_d^{'}} \right)$  
    $\displaystyle + \sum \frac{\langle M \rangle}{M_v}
\left( - \frac{\bar{q_v}}{\b...
... M_{fall} (\rho_r^{'}) +
M_{src} (q_v^{'}) + D_{q_v^{'}} \right)
\Bigg\}
\Bigg]$ (1.3)

熱の式

$\displaystyle \DP{\theta^{\prime}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{\theta^{\prime}}{x}
- w^{\prime} \DP{\theta^{\pr...
...DP{\overline{\theta}}{z}
+ \Dinv{\overline{\Pi}} (Q_{cond} + Q_{rad} + Q_{dis})$  
    $\displaystyle + D_{\overline{\theta}} + D_{\theta^{\prime}}$ (1.4)

比湿の保存式

$\displaystyle \DP{q_d^{\prime}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{q_d^{\prime}}{x}
- w^{\prime} \DP{q_d^{\prime}}{...
...d}}{z}
- \frac{\overline{q_d}}{\overline{\rho}} \sum M_{fall}
(\rho_s^{\prime})$  
    $\displaystyle + D_{\overline{q_d}}
+ D_{q_d^{\prime}},$ (1.5)
$\displaystyle \DP{q_v^{\prime}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{q_v^{\prime}}{x}
- w^{\prime} \DP{q_v^{\prime}}{...
...v}}{z}
- \frac{\overline{q_v}}{\overline{\rho}} \sum M_{fall}
(\rho_s^{\prime})$  
    $\displaystyle + \Dinv{\overline{\rho}} M_{src} (\rho_v^{\prime})
+ D_{\overline{q_v}}
+ D_{q_v^{\prime}},$ (1.6)
$\displaystyle \DP{q_c^{\prime}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{q_c^{\prime}}{x}
- w^{\prime} \DP{q_c^{\prime}}{z}
- \frac{\overline{q_c}}{\overline{\rho}} \sum M_{fall}
(\rho_s^{\prime})$  
    $\displaystyle + \Dinv{\overline{\rho}} M_{src} (\rho_c^{\prime})
+ D_{q_c^{\prime}},$ (1.7)
$\displaystyle \DP{q_r^{\prime}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{q_r^{\prime}}{x}
- w^{\prime} \DP{q_r^{\prime}}{...
...\prime})
+ \Dinv{\overline{\rho}} M_{src} (\rho_r^{\prime})
+ D_{q_r^{\prime}}.$ (1.8)

ただし, $\bar{~}$ の付いた変数は水平一様な基本場であることを示し, $\prime$ の付いた変数は基本場からのずれを表す. $\langle R \rangle$, $\langle c_p \rangle$, $\langle M \rangle$ はそれぞ れ平均気体定数, 平均定圧比熱, 平均分子量であり, 一定値とみなす. $Q_{cond}$, $Q_{rad}$, $Q_{dis}$ はそれぞれ凝結加熱項, 放射加熱項, 散逸加熱項を表し, $M_{src}$, $M_{fall}$ はそれぞれ生成項, 落下項を表す. $M_{src}$, $M_{fall}$, $Q_{cond}$ の定式化については 1.2 節で詳述する. $Q_{rad}$ の定式化については 1.3 節で詳述する. $Q_{dis}$, $D_{\ast}$ の定式化については 1.4 節で詳述する.

エクスナー関数 $\Pi $

\begin{displaymath}
\Pi \equiv \left( \frac{p}{p_0} \right)^{\langle R \rangle / \langle
c_p \rangle}
\end{displaymath} (1.9)

温位 $\theta $

$\displaystyle \theta \equiv T \left( \frac{p_0}{p} \right)^{\langle R \rangle / \langle
c_p \rangle}
= \frac{T}{\Pi}$     (1.10)

仮温度 $T_{v}$

$\displaystyle T_v \equiv \left( \frac{\langle M \rangle}{M_d} q_d
+ \sum \frac{\langle M \rangle}{ M_v } q_v
\right) T$     (1.11)

仮温位 $\theta_{v}$

$\displaystyle \theta_v
\equiv \frac{T_v}{\Pi}
= \left( \frac{\langle M \rangle}{M_d} q_d
+ \sum \frac{\langle M \rangle}{ M_v } q_v
\right) \theta$     (1.12)

密度 $\rho$

$\displaystyle \rho
= \frac{p}{\langle R \rangle \Pi \theta_v}
= \frac{p_0 \Pi^{\langle c_{v} \rangle / \langle R \rangle}}{\langle R
\rangle \theta_v}$     (1.13)

音波速度 $\langle c_s \rangle $

$\displaystyle \langle c_s \rangle^2
= \frac{\langle c_p \rangle \langle R \rangle
\Pi \theta_v}{\langle c_v \rangle}$     (1.14)

1.2 雲微物理過程のパラメタリゼーション

本モデルでは Kessler(1969) のパラメタリゼーションと Tobie et al.(2003) に基づくパラメタリゼーションの 2 種類が用意されている1.1.

1.2.1 Kessler(1969) の雲微物理パラメタリゼーション

Kessler(1969) のパラメタリゼーションに基づき, 気相と凝結相を以下の 4 つのカ テゴリーに分ける.
記号 意味 内容
$q_{d}$ 非凝結成分の比湿 気体の非凝結成分
$q_{v}$ 凝結成分の比湿 気体の凝結成分
$q_{c}$ 雲水比湿 落下速度がゼロである粒子で,
    大気中の雲粒に対応する.
    通常 100 $\mu$m 以下の微小な流体粒子である.
$q_r$ 雨水比湿 有意な落下速度を持つ粒子で,
    大気中の雨粒または氷粒に対応する.

$\rho_d$, $\rho_v$, $\rho_c$, $\rho_s$, $\rho$ をそれぞれ非凝結成分の密 度, 凝結成分の密度, 雲水の密度, 雨水の密度, 全密度とすると, 比湿は以下の ように定義される.

$\displaystyle q_{d}$ $\textstyle =$ $\displaystyle \frac{\rho_d}{\rho},$ (1.15)
$\displaystyle q_{v}$ $\textstyle =$ $\displaystyle \frac{\rho_v}{\rho},$ (1.16)
$\displaystyle q_{c}$ $\textstyle =$ $\displaystyle \frac{\rho_c}{\rho},$ (1.17)
$\displaystyle q_r$ $\textstyle =$ $\displaystyle \frac{\rho_s}{\rho}.$ (1.18)

そして, 微物理素過程として以下を考慮する. ただし, これらの量は全て正の値として定義され, 水蒸気が直接雨水に凝結する過程は無視されている.
記号 内容
$CN_{vc}$ 凝結による水蒸気から雲水への変換 (condensation)
$EV_{cv}$ 蒸発による雲水から水蒸気への変換 (evaporation)
$EV_{rv}$ 蒸発による雨水から水蒸気への変換 (evaporation)
$CN_{cr}$ 併合成長による雲水から雨水への変換. 併合や水蒸気拡散により,
  雲粒子が雨粒の大きさにまで成長する (autocondensation)
$CL_{cr}$ 衝突併合による雲水から雨水への変換.
  大水滴が小水滴を衝突併合する (collection)
$PR_{r}$ 雨水の重力落下に伴う雨水混合比の変化率 (precipitation)

この微物理素過程を用いて (1.4), (1.6) - (1.8) 式を書き直すと, 以下のようになる.

$\displaystyle \DP{\theta^{\prime}}{t}$ $\textstyle =$ $\displaystyle - \left( u^{\prime} \DP{\theta^{\prime}}{x}
+ w^{\prime} \DP{\the...
...v}{ \langle c_{p} \rangle \bar{\Pi}}
\left( CN_{vc} - EV_{cv} - EV_{rv} \right)$  
    $\displaystyle + \Dinv{\bar{\Pi}} \left(Q_{rad} + Q_{dis}\right)
+ D_{\bar{\theta}}
+ D_{\theta^{\prime}}$ (1.19)
$\displaystyle \DP{q_v^{\prime}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{q_v^{\prime}}{x}
- w^{\prime} \DP{q_v^{\prime}}{...
...e{q_v}}{z}
- \frac{\overline{q_v}}{\overline{\rho}} \sum PR_r
(\rho_s^{\prime})$  
    $\displaystyle - \Dinv{\overline{\rho}} \left( CN_{vc} - EV_{cv} - EV_{rv} \right)
+ D_{\overline{q_v}}
+ D_{q_v^{\prime}},$ (1.20)
$\displaystyle \DP{q_c^{\prime}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{q_c^{\prime}}{x}
- w^{\prime} \DP{q_c^{\prime}}{z}
- \frac{\overline{q_c}}{\overline{\rho}} \sum PR_r$  
    $\displaystyle + \Dinv{\overline{\rho}}
\left( CN_{vc} - EV_{cv} - CN_{cr} - CL_{cr} \right)
+ D_{q_c^{\prime}},$ (1.21)
$\displaystyle \DP{q_r^{\prime}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{q_r^{\prime}}{x}
- w^{\prime} \DP{q_r^{\prime}}{...
...\overline{\rho}} \left( CN_{cr} + CL_{cr} - EV_{rv} \right)
+ D_{q_r^{\prime}}.$ (1.22)

ここで, $L_{v}$ は水の蒸発の潜熱[J kg$^{-1}$], $\langle c_{p} \rangle$ は定圧比熱[J K kg$^{-1}$], $\Pi $ はエクスナー関数である.

微物理素過程は以下のように定式化する.

水蒸気と雲水の間の変換: $-CN_{vc} + EV_{cv}$
 
雲水は粒が小さく, 水蒸気との間で瞬間的に飽和調節が起こるもの とする. すなわち, 移流などの項を計算した後の温度と水蒸気量が 過飽和状態となっている場合には, ちょうど飽和になる量の水蒸気 を凝縮させる. 一方, 移流などの項を計算した後に, 雲水が存在す るにも拘わらず未飽和になっている場所では, ちょうど飽和になる 量の雲水を蒸発させる.

雲水の併合成長: $CN_{cr}$
 
Kessler (1969) に従って, 以下のように与える.
    $\displaystyle CN_{cr} = \bar{\rho} \bar{q_{d}}
\left[ \frac{q_{c}}{q_{d}}
- \left( \frac{q_{c}^{\prime}}{\bar{q_{d}}}
\right)_0 \right] / \tau _{ac}$ (1.23)

ただし, $( q_c / q_d)_0$, $\tau_{ac}$ は併合成長に関する閾値, 時間スケールであり, それぞれ 0, 100 [s] とする.

雲水の衝突併合: $CL_{cr}$
 
Kessler (1969) に従って, 以下で定式化する.
    $\displaystyle CL_{cr} = 2.2 \bar{\rho} \bar{q_{d}} q_{c}^{\prime}
\left( \frac{\bar{\rho}}{\bar{q_d}} q_r^{\prime}
\right)^{0.875} .$ (1.24)

雨水の蒸発: $EV_{rv}$
 
Kessler (1969) に従って, 以下で定式化する.
$\displaystyle EV_{rv} =
4.85 \times 10^{-2} \bar{q_d}
\left( \frac{q_{vsw} - q_{v}}{q_d} \right)
\left( \frac{\bar{\rho}}{\bar{q_d}} q_r^{\prime}
\right)^{0.65}$     (1.25)

ただし $q_{vsw}$ は飽和比湿を表す.

雨水のフラックス: $PR_{r}$
 
雨水の重力落下による混合比の変化率は,
$\displaystyle PR_{r} = \DP{}{z}
\left( \bar{\rho} q_r^{\prime} V_{term} \right).$     (1.26)

であり, 雨水の終端落下速度 $V_{term}$ [m s$^{-1}$] は
$\displaystyle V_{term} = 12.2 \left( \bar{\rho} q_r^{\prime} \right)^{0.125}$     (1.27)

で与える.

1.2.2 Tobie et al.(2003) の雲微物理パラメタリゼーション

Tobie et al.(2003) は火星大気での CO$_{2}$ の雲物理の定式化について述べ ている. 雲粒は拡散成長のみによって成長すると仮定している. 雲粒の併合成長は考慮せず, $q_c = 0$ とみなす. すなわち微物理過程として以下を考慮する.

記号 内容
$CN_{vr}$ 凝結による水蒸気から氷への変換 (condensation)
$PR_{r}$ 氷粒の重力落下に伴う氷混合比の変化率 (precipitation)

この微物理素過程を用いて (1.4), (1.6), (1.8) 式を書き直すと, 以下のようになる.

$\displaystyle \DP{\theta^{\prime}}{t}$ $\textstyle =$ $\displaystyle - \left( u^{\prime} \DP{\theta^{\prime}}{x}
+ w^{\prime} \DP{\the...
...me} \DP{\bar{\theta}}{z}
+ \frac{L_s}{ \langle c_{p} \rangle \bar{\Pi}}
CN_{vr}$  
    $\displaystyle + \Dinv{\bar{\Pi}} \left(Q_{rad} + Q_{dis}\right)
+ D_{\bar{\theta}}
+ D_{\theta^{\prime}}$ (1.28)
$\displaystyle \DP{q_v^{\prime}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{q_v^{\prime}}{x}
- w^{\prime} \DP{q_v^{\prime}}{...
...e{q_v}}{z}
- \frac{\overline{q_v}}{\overline{\rho}} \sum PR_r
(\rho_r^{\prime})$  
    $\displaystyle - \Dinv{\overline{\rho}} CN_{vr}
+ D_{\overline{q_v}}
+ D_{q_v^{\prime}},$ (1.29)
$\displaystyle \DP{q_r^{\prime}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{q_r^{\prime}}{x}
- w^{\prime} \DP{q_r^{\prime}}{...
...Dinv{\overline{\rho}} PR_r
- \Dinv{\overline{\rho}} CN_{vr}
+ D_{q_r^{\prime}}.$ (1.30)

ここで $L_s$ は CO$_2$ の昇華の潜熱 [J K$^{-1}$ kg$^{-1}$] である.

以下, $CN_{vr}$ の取り扱いについて述べる. 本モデルでは単位質量の気相に含まれる凝結核の個数及び半径は空間的・時間的 に一様と仮定する. また雲粒の半径は各格子内において空間的に一定であると仮定する. 更に雲粒は球形の凝結核を核として形成され, 雲粒自身も球形となると仮定する. このとき

\begin{displaymath}
\frac{4}{3} \rho_I \pi ( r_d^3 - r_{aero}^3) N = q_r \rho
\end{displaymath} (1.31)

という関係式が成り立つ. ここで $\rho_I$ は CO$_2$ 氷の密度, $r_d$ は雲粒半径, $r_{aero}$ は凝結 核の半径, $N$ は単位体積当たりの凝結核の数密度である. 本モデルでは $\rho_I = 1.565 \times 10^3$ [kg/m$^3$] と与え, $r_{aero}$, $N$ は実験に応じて与える. 雲粒の雲粒が拡散によって成長する場合の単位時間単位体積当たりの凝結量 $CN_{vr}$ は以下のように表される.
\begin{displaymath}
CN_{vr} = \frac{4 \pi r_d N}{R_h + R_m} (S - 1).
\end{displaymath} (1.32)

ここで $R_h$, $R_m$, $S$ はそれぞれ熱輸送に関する定数, 質量輸送に関する 定数, 飽和比であり
$\displaystyle R_h$ $\textstyle =$ $\displaystyle \frac{L^2}{k R T^2} = \frac{L^2}{k R \theta^2 \Pi^2},$ (1.33)
$\displaystyle R_m$ $\textstyle =$ $\displaystyle \frac{R T}{D p_{\ast}} = \frac{R \theta \Pi}{D p_{\ast}},$ (1.34)
$\displaystyle S$ $\textstyle =$ $\displaystyle \frac{p}{p_{\ast}} = \frac{{p_0}^{\frac{c_p}{R}} \Pi}{p_{\ast}}$ (1.35)

と表される. 但し $k$, $D$, $p_{\ast}$ はそれぞれ熱拡散係数, 分子拡散係数, CO$_2$ の 飽和蒸気圧である. 本モデルでは Tobie et al. (2003) 同様に主成分凝結系では $R_h \gg R_m$ と して
\begin{displaymath}
CN_{vr} = \frac{4 \pi r_d N}{R_h} (S - 1).
\end{displaymath} (1.36)

と表す. CO$_2$ の飽和蒸気圧については半経験式である Antoine の式
\begin{displaymath}
p_{\ast} = \exp \left( A - \frac{B}{T - C} \right)
\end{displaymath} (1.37)

を用いて定める(Antoine, 1888). ここで $A$, $B$, $C$ は実験により定まる係数であり, CO$_2$ の場合 $A =
27.4$, $B = 3103$, $C = -0.16$ である(化学工学会, 1999). 火星大気環境における凝結を想定すると $O(T) \sim 150 $[K] であるので, $T
\gg C$ と近似して
\begin{displaymath}
p_{\ast} \approx \exp \left( A - \frac{B}{T} \right)
\end{displaymath} (1.38)

とする.

以下, 単位時間体積当たりの雲粒落下量 $PR_r$ の取り扱いについて述べる. Tobie et al.(2003) では雲粒落下を無視しているが, 本パラメタリゼーション では考慮する. $PR_r$ は Kessler(1969) と同様に, 雲粒の終端速度 $V_{term}$ での移流 として表現する. 即ち

\begin{displaymath}
PR_r = \DP{}{z} \left( \rho_s V_{term} \right)
\end{displaymath} (1.39)

と表す. 終端速度 $V_{term}$ については球形粒子に関する Stokes 則を適用して
\begin{displaymath}
V_{term} = C_{sc} \frac{2 r_d^2 g \rho_I}{9 \eta}
\end{displaymath} (1.40)

と表す. ここで $C_{sc}$ は微小な粒子における Stokes 則からのずれを補正する係数 (Cunningham 補正係数)であり,
\begin{displaymath}
C_{sc} = 1 + 1.255 \frac{\lambda}{r_d}
\end{displaymath} (1.41)

と表される(Cunningham, 1910). $\lambda$ は CO$_2$ の平均自由行程であり,
\begin{displaymath}
\lambda = \frac{k_B T}{\sqrt{2} \pi \sigma^2 p_{CO_2}}
\end{displaymath} (1.42)

と表される. $k_B$ は Boltzmann 定数, $\sigma$ は CO$_2$ 分子の直径, $p_{CO_2}$ は CO$_2$ の分圧であり, $k_B = 1.38 \times 10^{-23}$ [m$^2$ kg s$^{-2}$ K$^{-1}$], $\sigma = 3.3 \times 10^{-10}$ [m] である(Golden and Sircar, 1994). $\eta$ は粘性係数であり, Sutherland の公式
\begin{displaymath}
\eta = \eta_{ref} \left( \frac{T_{ref} + C_{CO_2}}{T + C_{CO_2}} \right)
\left( \frac{T}{T_{ref}} \right)^{3/2}
\end{displaymath} (1.43)

で表現する(Sutherland, 1893). $\eta_{ref}$, $T_{ref}$, $C_{CO_2}$ はそれぞれ粘性係数の基準値, 温度の基 準値, CO$_2$ に関する Sutherland 定数であり, $\eta_{ref} = 1.47 \times
10^{-5}$ [Pa $\cdot$ s], $T_{ref} = 293$ [K], $C_{CO_2} = 240$ [K] と与 え る(理科年表, 2004).

1.3 放射加熱項の表現

放射加熱項 $Q_{rad}$ は正味の上向き放射フラックス $F_{net}$ を用いて 以下のように表される.

\begin{displaymath}
Q_{rad} = - \frac{1}{\overline{\rho} \langle c_{p} \rangle }
\DD{F_{net}}{z}
\end{displaymath}

本モデルでは $F_{net}$ は陽に計算せず, $Q_{rad}$ は高度のみに依存する パラメタとして与える.

1.4 乱流混合のパラメタリゼーション

1.4.1 運動方程式中の拡散項

Klemp and Wilhelmson (1978) および CReSS (坪木と榊原, 2001) と同様に, 1.5 次のクロージャーを用いることで粘性拡散項は以下のように書ける.
$\displaystyle D_{u_{i}}$ $\textstyle =$ $\displaystyle - \DP{}{x_{j}} \overline{(u_{i}^{\prime} u_{j}^{\prime})}$  
  $\textstyle =$ $\displaystyle - \DP{}{x_{j}}
\left[
- K_{m} \left(\DP{u_{i}}{x_{j}}
+ \DP{u_{j}}{x_{i}}\right)
+ \frac{2}{3} \delta_{ij} E_{turb}
\right].$ (1.44)

ここで $K_{m}$ は運動量に対する乱流拡散係数であり, $E_{turb}$ は サブグリッドスケールの乱流運動エネルギー
$\displaystyle E_{turb} = \frac{1}{2}
\overline{(u^{\prime})^{2} + (w^{\prime})^{2}}
= \frac{{K_{m}}^{2}}{C_{m}^{2} l^{2}}$     (1.45)

である. Deardorff(1975) に従い, $C_{m} = 0.2$ とする.

1.4.2 熱力学の式の拡散項

Klemp and Wilhelmson (1978) および CReSS (坪木と榊原, 2001) と同様に, 1.5 次のクロージャーを用いることで温位の粘性拡散項は以下のように書ける.
$\displaystyle D_{\theta}$ $\textstyle =$ $\displaystyle - \DP{}{x_{j}} \overline{u_{i}^{\prime} \theta^{\prime}}$  
  $\textstyle =$ $\displaystyle - \DP{}{x_{j}} \left(K_{h}\DP{\theta}{x_{j}}\right)
.$ (1.46)

ここで $K_{h}$ は温位に対する乱流拡散係数である.

1.4.3 乱流運動エネルギーの式

Klemp and Wilhelmson (1978) および CReSS (坪木と榊原, 2001) と同様に, 1.5 次のクロージャーを用いることで, 乱流エネルギーの時間発展方程式は以 下のように書ける.

$\displaystyle \DP{E_{turb}}{t}$ $\textstyle =$ $\displaystyle - u^{\prime} \DP{E_{turb}}{x} - w^{\prime} \DP{E_{turb}}{z}$  
    $\displaystyle - 3 \frac{g C_m l}{\overline{\theta_v}} E_{turb}^{\frac{1}{2}}
\D...
...ft( \DP{u^{\prime}}{x} \right)^2
+ \left( \DP{w^{\prime}}{z} \right)^2
\right\}$  
    $\displaystyle + C_m l E_{turb}^{\frac{1}{2}} \left( \DP{u^{\prime}}{z}
+ \DP{w^...
...2
- \frac{2}{3} E_{turb}
\left( \DP{u^{\prime}}{x} + \DP{w^{\prime}}{z} \right)$  
    $\displaystyle + \DP{}{x} \left( C_m l E_{turb}^{\frac{1}{2}} \DP{E_{turb}}{x} \...
...2}} \DP{E_{turb}}{z} \right)
- \frac{C_{\varepsilon}}{l} E_{turb}^{\frac{3}{2}}$  

ここで $C_{m} = 0.2$, 混合距離 $l = \left(\Delta x \Delta z
\right)^{1/2}$ とする. ただし $\Delta x, \Delta z$ はそれぞれ水平および鉛直格子間隔である. ただし,
$\displaystyle \overline{\theta_{v}} + \theta_{v}^{'}$ $\textstyle =$ $\displaystyle \bar{\theta_{v}}
\left\{
1 + \frac{\theta}{\bar{\theta}}
+ \frac{...
...}
- \frac{\sum q_{v} + \sum q_{c} + \sum q_{r}}
{1 + \sum \bar{q_{v}}}
\right\}$ (1.47)

である.

1.4.4 散逸加熱項の表現

散逸加熱項 $Q_{dis}$ は, 乱流運動エ ネルギーの散逸項をもとに, 以下のように与える.

\begin{displaymath}
Q_{dis} = \frac{1}{\langle c_{p} \rangle}\frac{C_{\varepsilon}}{l}
\frac{K_{m}^{3}}{(C_{m}l)^{3}}.
\end{displaymath} (1.48)

ここで $C_{\varepsilon} = 0.2$, $l=(\Delta x \Delta z)^{1/2}$ である.



... 種類が用意されている1.1
現在, 本モデルで用意されている Tobie et al.(2003) のパラメタリゼーション は火星大気計算でのみ使用可能である.

next up previous
: 2. 参考文献 : 湿潤大気における 2 次元非静力学モデルの定式化 : 目次
Yamashita Tatsuya 平成21年10月16日