プリミティブ方程式のソースを表示

{{参照方法|date=2019年6月}}
'''プリミティブ方程式'''（プリミティブほうていしき、[[英語]]：primitive equations）とは、大規模な[[大気]]の運動を記述する[[非線形]][[微分方程式]]群で、現在の[[気象予報]]において最も広く用いられている方程式。

以下の3つの主要な部分からなる。
*[[運動量保存の法則]] : 球体表面の水力学的運動を表現する[[ナビエ・ストークス方程式]]の変形。運動の水平スケールが鉛直スケールよりも十分大きい、[[静水圧近似]]の状態を前提（条件）とする。
*熱エネルギー方程式（[[エネルギー保存の法則]]） : 熱の出入りおよび、系全体の熱・エネルギーの変化とその状態を記述。
*[[連続の方程式]] : 質量の[[保存則]]を記述。

プリミティブ方程式からは、[[ラプラスの潮汐方程式]]を導出できる。この式は、緯線方向の運動を決定づける[[固有値]]である、線形のハフ関数を解に持つ。

一般的にプリミティブ方程式は、5つから6つの未知[[変数 (数学)|変数]]（他の変数に従属して決まるものもある）を持ち、これが空間・時間とともに変化していく。未知変数は、南北[[移流]]ベクトル、東西移流ベクトル、鉛直移流ベクトル、[[高度]]、[[気温]]または[[温位]]、[[湿度]]または[[比湿]]、[[気圧]]、空気の[[密度]]などである。

プリミティブ方程式の原型を最初に発表したのは、[[ヴィルヘルム・ビヤークネス]]であった<ref>[http://www.aip.org/history/sloan/gcm/prehistory.html Before 1955: Numerical Models and the Prehistory of AGCMs] アメリカ物理学会</ref>。[[ルイス・フライ・リチャードソン]]はこの式を用いて数値予報を行ったが失敗し、この欠点を補う[[準地衡風方程式]]を用いるようになったがこれも低緯度地域や小規模現象には適用できない欠点があり、結局はプリミティブ方程式が広く用いられるに至った。

移動性低気圧・高気圧などの総観スケールの現象を予報するには適した手法で、プリミティブモデルとして[[数値予報モデル]]に採用されている。一方、スケールの小さい現象は静水圧近似が適用できないため、非静力学モデルを用いる。

== プリミティブ方程式系 ==
プリミティブ方程式は、形は様々であるが、基本形として5つから7つの式を用いる。以下の通り。

*運動方程式 - ふつう、水平方向2つと鉛直方向1つの計3つの式を用いる。
*熱力学の第一法則
*気体の状態方程式
*連続の方程式（連続の式）
*水蒸気の輸送方程式
[[気象庁]]は、以上の7つの式（運動方程式を1つと数えて5つとする場合が多い）を用いている。

== プリミティブ方程式の変数 ==
*<math>u</math> : [[緯線]]上（球面上における東西方向）の速度。
*<math>v</math> : [[経線]]上（球面状における南北方向）の速度。
*<math>\omega</math> : 等圧面に対して[[鉛直]]方向の速度。
*<math>T</math> : [[気温]]。
*<math>\phi</math> : [[ジオポテンシャル]]。
*<math>f</math> : [[コリオリパラメータ]]。<math>2 \Omega \sin \phi</math>に等しい。<math>\Omega</math>は球面上の角速度(<math>2 \pi/24</math>[[ラジアン|rad]]/[[時間 (単位)|h]])、<math>\phi</math>は緯度。
*<math>R</math> : [[気体定数]]。
*<math>p</math> : [[気圧]]。
*<math>c_p</math> : [[比熱容量]]。
*<math>J</math> : 単位量・単位時間当たりの[[熱]]の移動。
*<math>W</math> : [[可降水量]]。
*<math>\Pi</math> : [[エクスナー関数]]。
*<math>\theta</math> : [[温位]]。

== 大気の運動によって生じる力 ==
[[気圧傾度力]]、[[重力]]、[[粘度|粘性]][[摩擦力|摩擦]]力を含めた、大気の運動によって生じる力を記述する。力は[[加速度]]をもたらし、それが式中に含まれる場合があるので注意。

気圧傾度力は、気圧の高いところから低いところへ移動する空気に働く加速度である。これを数式に表すと以下のようになる。

:<math>\frac{f}{m} = \frac{1}{\rho} \frac{\mathrm{d}p}{\mathrm{d}x}.</math>

[[質量]]をもった物体に対して、地球の中心へと向かう重力加速度9.81 m/s<sup>2</sup>が働く。

粘性摩擦力は以下の式で近似される。

:<math>f_r \simeq {f \over a} {1 \over \rho}
\mu\left(\nabla\cdot(\mu \nabla v) + \nabla(\lambda\nabla\cdot v) \right). </math>

以上の式を運動方程式にまとめると以下のようになる。

:<math>\frac{\mathrm{d}v}{\mathrm{d}t} = - (1/\rho) \nabla p - g(r/r) + f_{r}</math>

:<math>g = g_{e}</math>

よって、1つの系における方程式から解を導くには、6つの方程式を用いて、観測結果から得た数値を6つの変数に代入する必要がある。
*<math>\frac{\mathrm{d}v}{\mathrm{d}t} = - (1/\rho)\nabla p - g(r/r) + (1/\rho)\left[\nabla\cdot (\mu \nabla v) + \nabla(\lambda \nabla\cdot v)\right]</math>
*<math>c_{v} \frac{\mathrm{d}T}{\mathrm{d}t} + p \frac{\mathrm{d}\alpha}{\mathrm{d}t} = q + f</math>
*<math>\frac{\mathrm{d}\rho}{\mathrm{d}t} + \rho\nabla\cdot v = 0</math>
*<math>p = \rho RT</math>

== 出典 ==
*[https://www.metsoc.jp/tenki/pdf/1997/1997_04_0221.pdf プリミティブ方程式大気モデル]
*[http://www.ny.airnet.ne.jp/satoh/azr66pri.htm 基礎知識 プリミティブ方程式] 佐藤元、天気図と気象理論

{{Reflist}}
*Beniston, Martin. ''From Turbulence to Climate: Numerical Investigations of the Atmosphere with a Hierarchy of Models.'' Berlin: Springer, 1998.
*Firth, Robert. ''Mesoscale and Microscale Meteorological Model Grid Construction and Accuracy.'' LSMSA, 2006.
*Thompson, Philip. ''Numerical Weather Analysis and Prediction.'' New York: The Macmillan Company, 1961.
*Pielke, Roger A. ''Mesoscale Meteorological Modeling.'' Orlando: Academic Press, Inc., 1984.
*U.S. Department of Commerce, National Oceanic and Atmospheric Administration, National Weather Service. ''National Weather Service Handbook No. 1 - Facsimile Products.'' Washington, DC: Department of Commerce, 1979.

== 関連項目 ==
*[[天気予報]]
*[[数値予報モデル]]
*[[全球気候モデル]]

{{Physics-stub}}
{{Climate-stub}}

{{DEFAULTSORT:ふりみていふほうていしき}}
[[Category:流体力学の方程式]]
[[Category:地球物理学]]
[[Category:気象事業]]
[[Category:大気力学]]