非静力高解像度の次世代全球大気の気候力学数値モデルの構築に向けて

2000 Fiscal Year Annual Research Report

• Project/Area Number: 12874053
• Research Institution: Tohoku University
• Principal Investigator: 余 偉明 東北大学, 大学院・理学研究科, 助教授 (60251716)
• Co-Investigator(Kenkyū-buntansha): 岩崎 俊樹 東北大学, 大学院・理学研究科, 教授 (80302074)
• Keywords: 非静力学近似 / 全球力学数値モデル / 球座標系 / 格子モデル / セミ・インプリシット / 極の問題

Research Abstract

平成12年度には、コリオリ力を考慮した非静力学の全球大気の気候力学数値モデルの骨組みである力学フレームを作り上げた。方程式系は簡単でしかも山岳の効果を精度良く表現できるように、鉛直座標は地形に沿う座標ではなくジオイド面からの高度とした。また、将来、全球モデルの高精度化、高分解能化に伴い、積雲対流、特に熱帯深い積雲対流や山岳・成層の効果による重力波を直接にかつ精密に扱うことができるように、非静力学近似を採用した。本全球モデルは将来の作動環境として主記憶分散型の超並列計算機に想定しているので、適性の良い格子モデルを採用した。今年度は、我々が今までに開発してきた球座標系の数値モデルをベースにコリオリ力を入れてコーデイングし、積分時間に音波など高い周波数に関係する項をインプリシットに、ほかの低い周波数に関する項にエクスプリシットに解くというセミ・インプリシットの時間積分法を実現した。結果として積分時間間隔を長くとって計算時間短縮の実現を計ることもできた。球座標系を用いて差分格子点法でコーデイングすることによるる極の問題(pole problems)"を克服するため、差分格子モデルには変数をそれぞれ半格子分に離れた別々の点に定義するいわゆる食い違い格子(staggered grid)が用いられるとこによって、両極の特異点問題をうまく解決した。以上の時間積分、空間積分、安定性と効率性のことを念頭に置きながら、最初の一年目にはコーデイングに専念し、基本設計の数値コードを完成をした。なお、上述の結果は国際学会、書物と雑誌に公表している。

Research Products

• [Publications] W.Sha and K.Nakabayashi: "On the structure and formation mechanism of the spiral Taylor-Gotler vortices in spherical Couette flow"J.Fluid Mech.. 431. 323-345 (2001)
• [Publications] K.Nakabayashi and W.Sha: "Physics of Rotating Fluids Lectures Notes in Physics"Springer. 439 (2000)