DamBreakはオイラーの方程式を基礎方程式として、FDS系の解法が一番良い気がします。コロケート格子で。
スタッガードスキームだと、そのまま離散化すれば良いですが、ショックキャプチャリングがうまくいかなくて、若干の数値振動が残ります。
シビルエンジニアリングの世界ではほとんど1次風上が多いですね。
正直気に食わないですが。。。
FDS系の解法はXojoで書いてみることにします。そのうち、ソース晒しますね。
今のところbetter C的なC++で書いてます。
もう少し考えてみます。
0 件のコメント:
コメントを投稿