計算結果は鬼怒川氾濫（破堤付近）の痕跡を旨く説明できている(と思う).

ようやく計算が出来るようになったが考慮不足が沢山残っている。 例えば降雨量などは考慮できていない。 つまり、適当に河川上流から適当な流量を流しているだけ。 ただし、シミュレーションの計算は方程式どおりにやっているので計算は一応はちゃんとやって…

河川氾濫で降雨の影響を考慮するにはタンクモデルというやつで計算するようだが面倒なので適当に地形全体に雨を降らせる。 後は勝手に地形に沿って増水するはずだと考えた。 しかし、この考えは通用しない事がわかった。 ※あくまで今の実装ではだが、、、 例…

基盤地図も読み込めるようになった。 googleMapでマッピング。 わりとリアルな感じ。 先日書いたように上流からドバッと水を河川に流して見たがいい感じになっている。昨年、鬼怒川の堤防が決壊して大変な災害がおきているのでその近辺の５ｍメッシュ。 で、…

基盤地図を読み込めるようにした。 どうも河川の標高は取得できない。つまりDEM上は不定になっていてGeoTiffでみると真っ黒。多分「水」があるから河底標高（地盤高）は測定できないのだろう（多分、、、）。 仕方が無いので周囲の標高から攻め込んで穴埋め…

今日は複数の河川がある場合に対応したのと河川が合流するケースにも対応。それと地形に対して「堤防」の設置（属性）も追加した。 合流部は河川Aなのか河川Bなのか？ とりあえず平均値にする事にした。 変な話だがN本の河川が合流するケースも考慮。 ここま…

iRICを色々と試しています。 凄い。 ふと、思ったのは河川だって津波・高潮シミュレーションと同じ方程式でやるんじゃないかと。 ソルバは出来ているので試してみようと思ったが初期条件をどうしたらいいのか、、、 津波・高潮シミュレーションでは海洋（水…

津波・高潮シミュレーションだが前にも書いたように私のヤツはレギュラー格子なので色々と苦労しているのだがまだ、境界条件が適切に扱っていない。 ※スタッガード格子だと境界条件は簡単になるようです。 で、流体でいうところの圧力の境界条件を誤魔化して…

バグが直ったので台風201330号(HAIYAN)の計算が終わったので こんな感じになった。計算は正しそうだ。

バグが直ったので再計算しました。 名古屋港の潮位変化(T.P 0) 潮汐は考慮していない。 名古屋港での最高潮位と時刻は概ね再現されていると思う。 若干、過剰評価しているので海面抵抗計算を補正する必要がありそうです。

やっと間違いが分かった。 後日、修正したコードで再計算してみようと思う。 さて、間違いは単純だったがこれの間違いが分かるのには相当時間が掛かった。 簡単に言うと、、、 int foo( int W, int H, int i, int j, Foo& f ); という関数の呼び出しを間違え…

計算してみた。 当時の名古屋港の観測値と全くといっていいほど合っていない事が分かった。 コンタで見る限りはそれなりの結果だったので今までは合ってると思っていたが、、、 ふと、思って名古屋港の潮位変化をグラフにしてみたらたった５０ｃｍしか上昇し…

津波シミュレータを少し改良しました。 今までは断層破壊が瞬時に起こるか時間差で起こるかが指定できましたが、ある断層の破壊が進行している途中で別の断層の破壊がはじまるようなパターンは計算出来ませんでした。 しかも時間差は１０秒の開きを最低限必…

計算が途中で思わぬ数値振動が発生して発展してしまうケースがあった。 色々と対応策を考えてやってみたが結果は良くなるどころか悪化するだけ。 数値フィルタを使うとうまく行ったが駄目ケースがありフィルタリングを小まめにやるようにした。 しかし、本来…

海底地形は2mメッシュで最低水深10m、最大水深60mで計算してみた。 ちゃんとソリトン分裂が確認できるではないか！！。 非線形長波理論でも分裂が若干確認できるがこれをおそらく数値分散によるものだと思われる。 ソリトン分裂に特有の前傾が起きていること…

格子間隔5m、最大水深10m、最低水深1m、α=0.5、β=0.9で概ね安定して計算が出来た。つまり、この場合は時間ステップを半分にしてやれば修正係数は不要になる。が、格子間隔2m、最大水深90m、最低水深1mではやはり数値発散してしまった。 安定するまでα、βの値…

どうせシロートの作るプログラムだ。分散項の影響精度なんてそんなに必要ない。 と思えば首尾一貫して陽解法でやれるはずだ。 思い出そう。それは時間微分項が無い場合（定常の場合）はうまく解けている。 あ、もちろん精度は別の話。数値発散しないで計算で…

難問に突き当たったがそれはおそらくは「手抜き」にある。 そもそも陰解法でアタックすべき問題を陽解法に拘ったためである。陰解法は大規模な連立方程式を解いて未知量を確定させるのだから首尾一貫して解が確定する。 ※連立方程式が解けるかどうかとかそう…

先日の続き。 数値発散した原因は の評価である事までは明らかだ。しかし、よく見てみると第２項はΔtが掛かっているので時間ステップの調整次第で大きくはならないし、第３項は分母がメッシュサイズの３乗でそれなりに大きくなってこの項が巨大化するのを抑…

先日の続き。 10mメッシュと2mメッシュは殆ど瞬時に数値発散してしまった。原因は直ぐに判ったが非常に悩ましい。（悩ましいという意味については次第に明らかになるので） まず、方程式を見てみる。X方向の流量計算は次のような式になっている。 赤い枠の中…

津波解析では線形長波方程式が簡易評価で使われているが水深の浅い海域等では非線形項が無視できなくなり一般には非線形長波方程式が使われている。 しかし、この非線形長波方程式でもってしてもソリトン分裂等は評価できない。 津波ではソリトン分裂等の現…

私のこれまでの記事（内容）でいろいろと誤解が生じているようです。 その一つが被害評価です。 誤解の多くはこのシミュレーションが津波だけに限定した評価であることを忘れてしまっているようです。 もう少し書いておくと 直径10kmの彗星が激突した場合は…

Google Earth用ファイル[download]

計算結果のレンダリング結果（キャプチャ）です。 アニメーション1 アニメーション２（深い海域に衝突した場合）

前回の続き。 前回の冒頭にに書いたけど、、、 計算結果はかなり悲劇的かつ衝撃的な内容となった。 なので気分を害しそうな方は読むのを止めて頂きたい。 隕石衝突による津波の発生は概ね次のような段階を踏む 隕石衝突による波動きの発生 クレータへの海水…

これまでに自作してきた津波シミュレータだが隕石衝突における津波に関しても適用してみるとどうなるか試してみた。 計算結果はかなり悲劇的かつ衝撃的な内容となった。 なので気分を害しそうな方は読むのを止めて頂きたい（今回は基礎的な内容なので大丈夫…

陸上への遡上（浸水域）も計算してみた。 画像はタクロパン付近

以下は独自シミュレーションの結果です。 計算手法の妥当性、計算結果の信憑性、精度などは評価していない点に注意して頂きたい。また以降に述べられているコメントは独自のコメントであり何かを代表もしくは断定するものではありません。 動画は以下から見…

なぜ津波シミュレーションなのか？ そう聞かれることがある。 しかし、ほとんどそれに付いて説明をしたことがないし、その気もない。 頭の中では「どうせわかりゃしない」と思ってしまう。 ２０１１．３．１１、その日まで遡らなければならない。 名古屋にい…

1707 年宝永地震の津波波源域（Furumura, Imai,Maeda, 2011)のモデルで計算。 信頼性（相田指標）については未評価。