: 10 地球定数
: DCAPM4 第1部 数理モデル化
: 8 鉛直フィルター
支配方程式系の導出
.0
.00
.0
.00
全大気は, ともに理想気体である乾燥空気および水蒸気から成る混合大気とす
る. 雲水量は無視する. また, 水蒸気量が全大気に占める割合は小さいと仮定
し, 全大気の定圧比熱を乾燥大気の値で近似する.
水蒸気量の保存については, 凝結および蒸発による生成消滅を考慮する. しか
し, この量が全大気に与える効果は小さいとし, 全大気の質量保存則, 運動エ
ネルギー保存則, 全エネルギー保存則に影響を及ぼさないとする.
重力加速度は地球中心に向いていると仮定する. また, 運動の水平スケールが
鉛直スケールよりもかなり大きい運動を想定し, 静力学平衡近似を行なう. さ
らに, 運動は地球表面付近に限られることを仮定して近似を行なう.
.0
.00
.0
.00
方程式系は 6 本の予報方程式と 1 本の診断方程式からなる. 予報方程式は,
全質量の連続の式, 水蒸気量の式, 運動方程式(3 成分), 熱力学の式からなる.
これらは, それぞれ, 全質量保存則, 水蒸気量の保存則, 全質量に関する運動
量保存則, 全質量に関する全エネルギー保存則から導出する. 診断方程式には,
理想気体の状態方程式を用いる.
.1
!!注意: この Appendix 中では導出の都合上, 乾燥空気の気体定数を
定圧比熱を
とし, 全大気の気体定数を
とおいた. しかし, 本文
中では, 乾燥空気の気体定数を
, 定圧比熱を
と表記している.
乾燥空気, 水蒸気の状態方程式はそれぞれ
である. ここで
,
はそれぞれ乾燥空気および水
蒸気に関する量であることを示す. したがって, 全圧
は,
となる. ここで,
は比湿, であり,
,
である. したがって,
全大気の状態方程式は,
 |
( .5) |
ただし,
である. あるいは, 仮温度
を用いれば,
 |
( .6) |
全大気の質量保存則は, 水蒸気の生成消滅を無視すれば,
.2
 |
( .7) |
あるいは, ラグランジュ形式で記述すれば,
 |
( .8) |
水蒸気密度
に対する質量保存則は, 単位時間単位体積あたりの生成
消滅量を
とすれば,
 |
( .9) |
比湿
に関する式は, 原理的には式(
) と式(A.9) から得ることができる. しかし,
今の場合, 式(A.7)で水蒸気の生成消滅を無視したので,
正しくは得られない. そこで比湿の生成消滅に関する項を改めて
と定
義する.
 |
( .10) |
運動量保存則は, 水蒸気の生成消滅にともなう運動量変化を無視すれば次のよ
うに書ける.
 |
( .11) |
ここで,
は圧力,
は粘性応力テンソル,
は地球
の引力によるポテンシャル
.3 ,
はその他の外力項である. あるいは連続の式
を用いてラグランジュ形式で記述すると
 |
( .12) |
となる. ここで, 粘性項と外力項を
とおき, さらにベクトル表示する
 |
( .13) |
単位質量あたりの全エネルギーは, 運動エネルギー
と内部エ
ネルギー
およぼポテンシャルエネルギー
の和で表
現される. この時間変化率の式は, 水蒸気の生成消滅による影響を無視すれば,
![\begin{displaymath}
\DP{}{t}
\left[ \rho
\left( \frac{1}{2} \Dvect{v}^2
+...
...v_j
+ p v_j - \sigma_{ij}v_i
\right]
= \rho Q + F'_i v_i,
\end{displaymath}](img418.png) |
( .14) |
である. ここで,
は外界からの加熱率である. 一方, 運動エネルギーとポ
テンシャルエネルギーの和の保存式は, 運動量保存式 (
) に
をかけ連続の式を用いて変形することで得られる. 変形の際に
は
であるとしている.
.4
 |
( .15) |
となる. 式 (A.14) から式 (
) を引き去ると, 次のように内部エネルギーの式が
得られる.
 |
( .16) |
連続の式を用いてラグランジュ形式に書き直せば
 |
( .17) |
ここで, 外界からの加熱の項と粘性による加熱の項をまとめて
とおい
た.
内部エネルギーを温度を用いて表現すると
であ
る. さらに状態方程式 (A.5) を用いて式 (A.17)
を変形する.
であることに注意すれば
 |
( .18) |
となる. ここで,
を乾燥空気の定圧比熱
(定数) で近似すると
.5次の熱力学の式を得る.
 |
( .19) |
.5
.55
.5
.55
方程式系を, 一定の自転角速度
で回転する回転系に変換す
る.
慣性系における時間微分を添字 a で, 回転系を添字 r で表現する. このとき,
任意のスカラー
に対して,
 |
( .20) |
が成りたつ.
.6
任意のベクトル
に対する慣性系および回転系での微分は次の関
係をもつ.
 |
( .21) |
(証明) 任意のベクトル
を, 慣性系では
|
 |
 |
(.22) |
と表し, 回転系では
|
 |
 |
(.23) |
と表す.
時間微分をとると
(証明終り)
ここで
(
は位置ベクトル ) とおけば慣
性系での速度
(これまでの
) は回転系での速度
を用いて次のように表すことができる.
 |
( .25) |
さらに, 式(A.21) で
とお
けば, 速度の時間微分項は
 |
( .26) |
と変換できる.
変換の式 (A.26) を用いて運動方程式を回転系で記述する.
 |
( .27) |
ここで, 重力加速度
を定義すれば, 運動方程式は
 |
( .28) |
となる.
連続の式および熱力学の式においては, ラグランジュ微分が作用している密度
および温度は座標変換に無関係なスカラーであるため, その時間微分の形は変
わらない. 連続の式は, 速度場の発散を含むが, これは座標変換によっても値
は変わらない. したがって, これらの式は形を変えない.
.7
.7
.77
.7
.77
一般の直交曲線座標
において, スカラー
およびベクトル
は次のように表現
される. なお,
は各軸方向の規模因子であり, 各軸方向の基底ベクトル
は
とする.
 |
( .29) |
![\begin{displaymath}
\Ddiv \Dvect{A}
= \frac{1}{h_1 h_2 h_3}
\left[ \DP{}{\x...
...\xi_2} ( h_1 h_3 A_2)
+ \DP{}{\xi_3} ( h_1 h_2 A_3)
\right],
\end{displaymath}](img473.png) |
( .30) |
![\begin{displaymath}
\nabla^2 \bullet
= \frac{1}{h_1 h_2 h_3}
\left[ \DP{}{\...
...eft( \frac{h_1 h_2}{h_3} \DP{\bullet}{\xi_3} \right)
\right],
\end{displaymath}](img474.png) |
( .31) |
![\begin{displaymath}
\Drot \Dvect{A}
= \left( \frac{1}{h_2 h_3}
\left[ \DP{(h_...
...DP{(h_2 A_2)}{\xi_1} - \DP{(h_1 A_1)}{\xi_2} \right]
\right),
\end{displaymath}](img475.png) |
( .32) |
 |
( .33) |
![\begin{displaymath}
\DD{\Dvect{v}}{t}
= \sum^3_{k=1} \Dvect{e}_k
\left[ \DP{v...
...{v_k}{h_k} \frac{1}{h_j} \DP{h_k}{\xi_j} \right) v_j
\right].
\end{displaymath}](img477.png) |
( .34) |
重力加速度
が地球中心を向いているとみなして, 方程式系を球
座標
に変換する. 回転
系に固定した直交直線座標
との関係は
である. ここで,
は緯度,
は経度,
は鉛直座標であ
る. また, 基底ベクトルを
, 速度ベクトルを
で表す.
各方向の規模因子( scale factor )は
 |
( .38) |
したがって, スカラー
およびベクトル
に関する微分表現は次のよう
になる.
 |
( .39) |
![\begin{displaymath}
\Ddiv \Dvect{A}
= \frac{1}{r^2 \cos \varphi}
\left[ r \...
...hi A_{\varphi})
+ \cos \varphi \DP{}{r} ( r^2 A_r )
\right],
\end{displaymath}](img493.png) |
( .40) |
![\begin{displaymath}
\nabla^2 \bullet
= \frac{1}{r^2 \cos \varphi}
\left[ \D...
...}{r} \left( r^2 \cos \varphi \DP{\bullet}{r} \right)
\right],
\end{displaymath}](img494.png) |
( .41) |
 |
( .43) |
自転角速度ベクトルの表現は次のようになる.
したがって, 運動方程式は
連続の式は
 |
( .50) |
熱力学の式は
 |
( .51) |
状態方程式は
 |
( .52) |
水蒸気の式は
 |
( .53) |
ここで,
 |
( .54) |
である.
.7
.77
.7
.77
.7
.77
鉛直方向の運動方程式に対し, 静力学平衡近似を行なう.
 |
( .55) |
このとき, 運動エネルギーの保存則を考慮して, 水平方向の運動方程式に対し
ても近似を施す. 運動エネルギーの式は, 運動方程式の各成分にそれぞれ
をかけることで得られる.
コリオリの力およびメトリック項は同じ番号のもの同士で打ち消しあって, 運
動エネルギーの時間変化に寄与しないことがわかる.
.8したがって, 静力学平衡近似の際に鉛直成分の式から落とした項(2),(4),(5)
に対応した水平成分の式の項も取り除く. これにより, 運動方程式の水平成分
は次のようになる.
ここで,
はコリオリパラメータ
であ
る.
大気の層が地球半径に比べて薄いことを仮定し, 方程式中の
を, 代表的
な地球半径
でおきかえる. また,
による微分はすべて海抜高度
による微分でおきかえる. このとき基礎方程式は次のようになる.
 |
( .59) |
 |
( .60) |
 |
( .61) |
 |
( .62) |
 |
( .63) |
 |
( .64) |
 |
( .65) |
ここで,
 |
( .66) |
 |
( .67) |
.8
.88
.8
.88
静力学平衡のもとでは, 気圧
は鉛直座標
に対し単調減少する関数で
ある. そこで, 鉛直座標を
から, 地表面気圧
で規格化した気圧座
標,
 |
( .68) |
に変換する.
と
の関係は, 静力学平衡の式を変形して得られる.
 |
( .69) |
- 座標から
- 座標への変換公式を示す.
鉛直微分
水平微分
時間微分
ラグランジュ微分はこれらを用いて,
ここで,
-座標鉛直速度
を定義する.
 |
( .75) |
式 (A.69) を重力ポテンシャル
を用いて書けば,
 |
( .76) |
水平の圧力勾配は, 式(A.71)および式(A.72) を
に対して適用し, 式(A.68) を用いれば次のように変換される.
![\begin{displaymath}
\frac{1}{\rho} \left( \DP{p}{\varphi} \right)_z
= R^d T_v \DP[][\sigma]{\pi}{\varphi} + \DP{\Phi}{\varphi}.
\end{displaymath}](img568.png) |
( .78) |
ここで
である. したがって, 運動方程式の水平成分は,
 |
( .79) |
 |
( .80) |
速度の発散は,
ここで,
![\begin{displaymath}
\Ddiv{\Dvect{v}_H}
\equiv \frac{1}{a \cos \varphi}
\DP[]...
...hi} \left( \DP{}{\varphi} (v \cos
\varphi ) \right)_{\sigma}.
\end{displaymath}](img580.png) |
( .82) |
ゆえに,
- 座標連続の式は次のように変換される.
したがって
を用いて記述すれば次のようになる.
 |
( .84) |
式 ( A.64 ) の右辺第1項は次のように変換される.
ここで,
 |
( .86) |
したがって, 熱力学の式はつぎのようになる.
 |
( .87) |
ここで,
座標における境界条件について述べる.
 |
( .88) |
すなわち,
は表面地形を表す. この境界条件を用いて, 静力学平衡
の式を鉛直積分することで, 任意の
における高度
を求める
ことができる.
 |
( .89) |
ここでは述べない.
連続の式を鉛直方向に
から
まで積分し,
に関する境界条件を用いれば, 傾向方程式とよばれる
の時間変化に関する式が得られる.
 |
( .90) |
この式を用いれば,
の情報がなくても地表面気圧の時間変化
を求めることができる. なお, ここでは後のことを考えて
を
と表現している.
については次節で改めて
定義する.
鉛直速度
は, 連続の式を鉛直方向に
から
まで積分することで診断的に得られる.
 |
( .91) |
.8
.88
.8
.88
渦度:
 |
( .92) |
発散:
 |
( .93) |
運動方程式の
の式に
を
作用し,
の式に
を
作用し差をとって変形すれば次の渦度方程式を得る.
運動方程式の
の式に
を作
用し,
の式に
を作用し和をとって変形すると次の発散方程式を得る.
ここで,
支配方程式系における変数を, モデル内部で用いている変数に変換する. まず,
を導入する. また速度場
は
,
に変換する.
.9 このとき, 水平風の渦度
と発散
は次のように
変換され, この表現をあらためて
および
と定義する.
水平風による移流は次のように変換される.
水平風による移流のもうひとつの記述を連続の式の変換のために示す.
これらを用いて, 方程式系を次のように変数変換する.
連続の式
.10
 |
( .102) |
渦度方程式
発散方程式
熱力学の式
水蒸気の式
仮温度
を次のように
のみに依存する場
と, そこからのずれ成分
にわけて記述する.
渦度方程式で
を含む項は次のように変形される。
発散方程式で
を含む項は次のように変形される.
熱力学の式の右辺第1-3項は次のように変形される.
以上を用いて程式系を記述すれば次のようになる.
連続の式
 |
( .110) |
静水圧の式
 |
( .111) |
運動方程式 .11
 |
( .112) |
 |
( .113) |
ここで,
熱力学の式
水蒸気の式
 |
( .117) |
式(A.17) で導入した
から粘性による寄与
を再び分離し,
とする. 一般に粘
性は運動方程式において適当なパラメタリゼーションによって表現する. また,
渦度, 発散, 温度, 水蒸気の式に対してそれぞれ水平拡散項
,
,
,
をつける.
この項の付加は主に数値的安定性の要請
によるものであるが, 物理的には後で行なう離散化のスケール以下の運動を表
現していると解釈できる. 最後に, 乾燥大気の気体定数および定圧比熱
,
をそれぞれ
,
のようにあらためて置きなおせば, 支
配方程式系 (3.1) -- (3.6) を得る.
- ...
- 94/04/13 石渡正樹
- ...
理想気体の状態方程式を用いる. .1
- 乾燥空気と水蒸気は, 同じ速度と温度をもつことを暗黙のうちに仮定
している. したがって, 水蒸気に関する運動量保存則および全エネルギー保存
則および状態方程式を考慮する必要がない.
- ... 水蒸気の生成消滅を無視すれば, .2
-
次で示すように水蒸気式では生成消滅を含めている. したがって, 全大気の質
量保存則は, 水蒸気の生成消滅が起きても全質量が保存するように, 乾燥大気
量が変化することを要請していることになる.
- ...
の引力によるポテンシャル .3
- これは遠心力を考慮しない地球の質量にのみ起因
したポテンシャル.
- ... であるとしている. .4
-
導出の過程を示す. 左辺第1項と第2項は次のように変形される.
また, 左辺第4項は次のように変形される.
- ... で近似すると .5
-
この近似には疑問が残る. 状態方程式においては, 気体定数
を
と
する近似は(仮温度
を導入することで)行なわなかった.
につい
てだけ近似するのは近似のレベルに一貫性がないように思われる.
- ...
が成りたつ. .6
-
これは自明のこととしたい. スカラー
の座標変換は座標変換テンソル
に依存しない(で同じ値をとる)からである. なお, Pedlosky (1987) では, ベ
クトルの変換公式を使ってスカラーの変換を証明している. ところがベクトル
の変換公式ではスカラーの変換公式を使っているので, 何がなんだかわからな
い.
一方, ベクトルの座標変換は, 座標変換テンソルとの積で表現される. したがっ
て, 座標変換テンソル自体が時間変化する場合, 当然ベクトルの時間微分は座
標変換テンソルの時間微分の影響を受ける.
- ... これらの式は形を変えない. .7
- ここで私が気になっているのは, 運動エネルギー保存の式が, 回転系へ
の移行によって形を変えることである.
- ...
動エネルギーの時間変化に寄与しないことがわかる. .8
-
遠心力を重力加速度から分離してエネルギーの式で考慮すると, この寄与はキャ
ンセルすることなく残る.
- ... に変換する. .9
のままでも渦度や発散にすれば極での特異性を回避できる
のではないか?
- ...連続の式
.10
- 発散はすべて
を用いて表現すべきだろう.
- ...
運動方程式 .11
- (2005/4/4 石渡)
の式の右辺第一項の符号は正
しいか?
の式の右辺第二項の符号は正しいか?
: 10 地球定数
: DCAPM4 第1部 数理モデル化
: 8 鉛直フィルター
Yasuhiro MORIKAWA
平成20年1月30日