趣味で計算流砂水理

Computational Sediment Hydraulics for Fun

コロケート格子の境界条件の設定の悩み

ずっと悩んでいるのですが,改めて.

非常にシンプルな一次元不定流計算の定常解です. バックウォーターの計算です.

nbviewer

問題は下流端流量です.

境界なので片側差分で計算にすると,少し大きな値になってしまいます.

圧力項dHdxの差分がまずいです.

単純に一つ上の断面と同値にすると問題ないのですが,非定常の場合はあまりよろしくないですよね.

上流端の河積(水深)も同様の問題が生じます.

一般的にこのような場合はどうやって解決するのでしょうか....高次精度とか?


なんかgistがうまく表示されない.

1DUpWindFlux_schemecheck