SedimentHydraulics https://blog.hatena.ne.jp/SedimentHydraulics/ 趣味で計算流砂水理 Computational Sediment Hydraulics for Fun Learning https://computational-sediment-hyd.hatenablog.jp/ 集中格子の場合 H(i-1)=H(i+1)=0.0,H(i)=1.0の場合，dH/dx=0 これは、各水面勾配の重みが0.５となり、打ち消しあって0になるようですね。 物理現象からすると、iの地点においてi+1とi-1の地点にそれぞれ流れが発生するようになるはずですが、 集中格子が原因でダメですね。 ちょうど、非圧縮生流体の解析のように、集中格子の場合のチェッカーボードの振動のようですね。 解決策 不等間隔格子 非圧縮生流体の場合、不等間隔格子を使うとチェッカーボードの振動は入らないようですが、 そもそもの解決になっていないように思います。 その場しのぎの処理 質量保存がどうなるかわかりません… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fcomputational-sediment-hyd.hatenablog.jp%2Fentry%2F2016%2F11%2F03%2F100750" title="水面勾配 - 趣味で計算流砂水理 Computational Sediment Hydraulics for Fun Learning " class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> Hatena Blog https://hatena.blog 2016-11-03 10:07:50 rich https://computational-sediment-hyd.hatenablog.jp/entry/2016/11/03/100750 1.0 100%