B. 音波減衰項について

appendix-b

本モデルで用いている時間方向の離散化方法(時間分割法)は, Klemp and Wilhelmson (1978) によって最初に提案された方法である. この方法を用いる と, 音速および移流に対する CFL 条件をそれぞれ満たしている場合でも計算 不安定を起こす場合がある. この影響は短い時間間隔で積分するステップ数 381#381 を増加させていくにつれ大きくなる (Skamarock and Klemp, 1992, 図 1 を参照). Skamarock and Klemp (1992) はこの計算 不安定を回避するため, 音波を選択的に減衰させるフィルターとして運動方程 式に音波減衰項を導入することを提案した.

音波減衰項は速度場の発散に対する拡散として作用する. このことは音波減衰項を含む線形化された基礎方程式

382#382 52#52 383#383 (143)
384#384 52#52 385#385 (144)
386#386 52#52 387#387 (145)
388#388 52#52 389#389 (146)
390#390 52#52 391#391  

から発散方程式
392#392 52#52 393#393 (147)
394#394 52#52 395#395  

を求めることで理解することができる.

音波減衰項は速度場の発散に対する拡散なので, 重力波成分も拡散させる可能 性が考えられる. しかし Skamarock and Klemp (1992) は線形化された基礎方 程式の分散関係を用いて, 音波減衰項の重力波成分への影響は小さいことを示 している. 以下では Skamarock and Klemp (1992) における議論の概要を示す. 簡単のため圧力方程式に現れる基本場の量は定数とし, 全ての変数を 396#396 のような解を持つと仮定して分散関係を 求めると,

397#397 (148)

となる. ここで 398#398 である. 純粋な内部重力波の分散関係式

399#399

400#400 を含む項に代入すると 400#400 を含む項は互いに打ち消しあ うことから, 音波減衰項の重力波成分への影響は小さいであろうと予想される. 実際に, 401#401 が小 さいと仮定して 402#402403#403 の巾で展開して音波減衰項の影響を 評価すると, 音波減衰項は重力波成分を減衰させる方向にはたらき, その減衰 率は非常に小さいことがわかる(Skamarock and Klemp, 1992, 図 5 参照).

音波減衰項の係数 400#400 の値を決めるにあたり, 考慮しなければならない ことが 2 つある. 1 つは音波減衰項そのものが計算不安定の原因とならない ようすることであり, もう 1 つは音波減衰項の重力波成分への影響が大きく ならないようすることである. 前者は音波減衰項を含む線形発散方程 式に示したように, 音波減衰項が発散方程式の拡散項となることから要請さ れる. 音波減衰項の時間積分は前進差分を用いて行われるので (uwpi:u_sabun, uwpi:w_sabunを参照), 計算不安定を起こ さないためには拡散項を前進差分で時間積分する場合の安定性条件

404#404 (149)

を満たさなければならない. 後者は前段落で紹介した Skamarock and Klemp (1992) の議論で用いた 403#403 が小さいという仮定から要請される [*]. 403#403 の最大値は

405#405


と与えられるので, 406#406 とするためには
407#407 (150)

でなければならない.音波減衰項係数の条件2右辺に現れる音速に対 するクーラン数は 1 より小さい値とするので, 音波減衰項係数の条 件2の条件を考慮すれば音波減衰項係数の条件1は自動的に満たされ ることになる.



Footnotes

... が小さいという仮定から要請される[*]
この条件に関する議論は Skamarock and Klemp (1992) ではなされて いない
Yamashita Tatsuya 2012-09-11