基于Lattice Boltzmann方法的洪水演进模拟

作者：刘宁宁; 时勇; 孙华林

来源：中国农村水利水电, 2022, (06): 133-137.

• Lattice Boltzmann方法
• 扩散波方程
• 平衡态分布函数
• 弛豫时间

采用离散方法求解圣维南方程组的简化形式来推求下游断面未来洪水过程是目前洪水演算中最为常见的水力学方法。但这种方法存在的数值不稳定、计算速度慢等问题，也是目前宏观方法描述流体运动的一个挑战。Lattice Boltzmann方法以粒子速度分布函数为研究对象，在求解非线性偏微分方程组时，将非线性偏微分方程组转化为Lattice Boltzmann演进方程，从而变为简单的线性方程，使得计算效率得到大大提高。研究选用线性扩散波方程来描述河道洪水波运动，采用对空间坐标不进行多尺度处理，而对时间坐标引入五个时间尺度的多尺度展开方式，建立了线性扩散波方程的Lattice Boltzmann方法 D1Q5速度模型，进而推求下断面洪水过程。文中分别采用Muskingum法、解析解法和Lattice Boltzmann方法的D1Q5速度模型对长江上游龙街—乔家河段、湖南资水柘溪水电站下游江南—夫溪河段进行洪水演算，3种方法的洪峰相对误差、峰现时间、确定性系数模拟的结果相似，从而验证了Lattice Boltzmann方法的准确性。同时分析了时间步长、空间步长和弛豫时间的选取对计算精度的影响。结果表明：Lattice Boltzmann方法可以很好的预测洪水过程，且计算效率较宏观方法得到大大提高，空间步长的取值比时间步长的取值对精度的影响大，在进行洪水演算时，应结合上下游断面长度，选取合适的空间步长，弛豫时间τ的取值不宜过大，建议在[1.5,3]范围内为宜。