Nays１D＋で１０年間程度の一次元河床変動解析を日流量単位の不等流計算ベースで進めています。d=3cm, 元河床勾配=約1/300、流量規模＝平均年最大で約100m3/s程度です。下流水位の境界条件の設定が難しいですが、一定勾配の場合は0.0028、その他、変動河床に対する等流条件でも検討しています。dtは数秒から数百秒で検討してます。
しかし、大きな流量前後で発散して値がNanとなってしまいます。
どのような安定化の方策があるでしょうか？
よろしくご指導ください。