潮波方程数值格式与摩阻系数选取

基于一维阻尼潮波传播方程解析解,从求解数值格式及Heuristic稳定性分析方面,讨论了数值解的精度、计算耗时和摩阻系数选取等问题。研究结果表明:1)Courant数小于1时,潮波方程显格式解的精度略高于隐格式解,计算耗时少于隐格式解;2)为减少计算耗时,潮波方程的隐格式解允许较大的时间步长,但解的精度有所降低,须通过减小底床摩阻系数以保证计算精度;3)隐格式解摩阻系数的选取与Courant数有关,Courant数越大,摩阻系数的选取值比实际值越小,通过理论分析结合数值试验得到了相应的关系式。这些研究结论对实际海域的潮波传播的数值模拟具有重要的应用价值。