摘要
深入认识大型江湖系统水动力特性,需要使用平面二维数学模型代替传统一维数学模型开展模拟和研究。计算精度的优先级高于计算效率,当以高分辨率计算网格作为前提时,大型江湖系统二维水动力模型“计算精度与效率”的矛盾转化为计算效率瓶颈。本文以荆江—洞庭湖系统为背景(剖分为32.8万个四边形单元),选用MIKE21模型和一种半隐式欧拉—拉格朗日模型,开展大型江湖系统显式、隐式水动力模型性能的比较研究,探讨大型江湖系统高分辨率二维模型的可实用性。研究结果表明,显式模型稳定性较差,允许的最大时间步长为0.8s;隐式模型稳定性好,时间步长可达60s以上。在常规16核工作站的效率测试结果表明,显式模型具有良好的可并行性,加速比随核心数量基本呈现线性变化;隐式模型(使用预测—校正分块并行计算方法)亦可取得非常大的加速比,在16核心并行条件下加速比Sp可达11.1。模拟荆江—洞庭湖系统1年的非恒定流过程,显式、隐式模型的耗时分别为411h、10.76h(16核),后者计算效率是前者的38.2倍,可满足大型江湖系统整体高分辨率平面二维水动力模拟的实用性要求。
-
单位华中科技大学; 长江科学院河流研究所