摘要:为高效模拟一维河网水动力学过程,克服汊点水位迭代法中每次迭代均需全局求解所致计算效率过低的问题,本文基于明渠非恒定流的格子Boltzmann BGK模型,提出了一种改进的一维河网水动力算法。其中,一维圣维南方程组由格子Boltzmann BGK模型数值求解,离散速度矢量采用D1Q3模型描述;河网耦合求解采用汊点水位迭代法实现,各分汊河段在汊点处均采用水位边界条件,每一时步内通过不断迭代修正汊点水位以满足汊点连接条件。因格子Boltzmann BGK模型本质上为显示算法,有效回避了采用隐式差分算法计算汊点净流量时需多次迭代推求汊点耦合水位的问题。利用环状河网恒定流、环状和树状组合河网非恒定流等典型河网算例验证了所提方法的准确性,通过与原始汊点水位迭代法进行比较说明本方法能大幅提高河网水动力模拟的计算效率。作为方法在实际场景中的应用,以汉江流域2021年秋汛实测水文过程作为输入条件,计算了汉江中下游平原河网的洪水过程,对比了汉江中下游各水文站实测、计算和商用软件模拟的洪水过程,通过与特征线法、Preissmann四点隐式差分法等传统一维解法对比,说明了所提方法可在不损失计算精度和确保汊点连接条件严格成立的前提下,大幅提高一维河网水动力求解效率。
文章目录
1 基于格子Boltzmann方法的一维河网水动力解法
1.1 控制方程
1.2 格子Boltzmann BGK模型
1.3 河网汊点水位迭代法
2 方法验证
2.1 环状河网恒定流算例
2.2 环状、树状组合河网非恒定流算例
3 汉江中下游河网水动力模拟
4 结 论