我一直在研究具有反射边界条件的浅水方程的求解器。我发现它发散得非常快。作为一种解决方法，我昨天注意到，如果我在每个步骤中平滑动量（水柱高度 * 速度）场作为后处理，则解决方案不再爆炸。所以似乎我在没有平滑的地方出现了一些高频误差项？他们通常来自哪里？
在 SWE 中出现的错误术语可能来自哪里？
计算科学
边界条件
有限体积
收敛
2021-12-20 07:17:52
1个回答
您是否使用可视化包绘制了解决方案，或者只是在您的域中间进行了切割？有时，只需查看您正在生成的解决方案，您就可以深入了解它的来源。根据您的描述，我会说您在执行边界条件时存在错误。
ETA：最后一句话可能是错误的，并在评论中重述我们的对话并放大Gresho 和 Lee 的精辟短语“Wiggles 正在告诉你一些事情”：即使你没有使用 Galerkin FEM 来离散化这一点，方法在数学上等同于几个 FVM 和 FDM，具体取决于所使用的正交规则，因此如果边界层中没有足够的网格，则会显示类似的病状。