前言
不稳定的数值模型是指某些类型的数值误差增长到解开始振荡的程度,或者误差变得如此之大以至于计算无法继续。在处理任何大小或复杂性的非定常流模型时,这是一个常见问题。
本文主要介绍几种方法,为排除模型的不稳定性作为参考依据,后续如果更好办法,会进行更新~欢迎大家留言或私信
以下因素会影响模型的稳定性和数值精度:
1.横截面间距。
2. 计算时间步长。
3. 数值解的θ加权因子。
4. 计算选项和公差。
5. 侧向结构/堰
6. 陡峭的溪流/混合流态
7. 下游边界条件
8. 横截面几何形状和表格属性
9. 桥梁和涵洞交叉口
10. 初始/低流量条件
11.床型下降。
12. 曼宁系数
13. 主河槽数据缺失或损坏
横断面间距
横截面不够:当横截面间隔很远,并且水力性能的变化很大时,解决方案可能会变得不稳定。通常,相距太远的横截面将导致额外的数值扩散,因为相对于距离的导数被平均到很长的距离。此外,如果横截面之间的距离如此之大,使得Courant数将远小于1.0,那么模型也可能变得不稳定。
变化横截面间距的示例如图所示。一个入流(绿虚线)和两个出流(带正方形的蓝色和黑色实线)。如图所示,随着横截面间距的增加,图像将显示一些数值衰减/扩散。
横截面太近。如果横截面太近,则相对于距离的导数可能被高估(计算为比它