2.2 シミュレーションプログラムの概要
ここで用いたシミュレーションモデルは平成8年度事業で作成されたものだが、その概要を以下に述べる。
2.2.1 数値計算手法
津波の計算は、?地震による地殻変動に伴う津波の発生、?震源域から沿岸までの伝播、?港湾域で摩擦減衰、波形変形、浸水、遡上、の3過程に分けて考えることができる。?については、地震学の分野で提唱された断層モデルを波源モデルとして適用する。単純な形の断層のずれによって海底面がどのように変位するかは、弾性論によって求めることができ、種々の断層パラメータを与えることにより海底の鉛直変位が計算される。この変位がある変動時間で生じたものとして、数値実験を開始する。
?については津波が浅水理論(水深に対して波長が充分長いという仮定に基づいた理論)に従って運動するものとする。ただし、外洋部では運動方程式中の移流項、摩擦項を省略した線型近似した線型長波式を用いる。支配方程式を差分化し、これを細かい時間間隔で解くことによって適当な格子間隔で水位と流速が計算される。格子は計算機容量、計算時間の節約をするために外洋部では荒い格子間隔、沿岸部に近づくにつれて細かい格子間隔をとる多重格子を用い、格子間で計算結果を引き継いで行く。
?についても浅水理論が基礎となるが、運動方程式中に移流項、摩擦項を込め、さらに汀線が水位の昇降に伴って移動することを考慮する。
2.2.2 基本方程式
単位質量の海水に対する運動方程式は、海水の非圧縮性及び静水圧近似を仮定して、次式のように表すことができる。
