我使用不连续伽辽金法求解底部有凹凸的稳定流动一维浅水方程。这种流动是无摩擦的。

我使用 runge-kutta 方法来逼近这些方程的时间导数，并使用 gauss-lobatto 求积来逼近源项的弱形式的积分方程。

我得到的解决方案总是崩溃并且与确切的解决方案不对应，我认为我可能错误地近似源项，但我不知道如何解决这个问题。

有什么建议？