[发明专利]一种运动波与动力波相结合的山区洪水过程数值模拟方法

摘要

本发明公开了一种运动波与动力波相结合的山区洪水过程数值模拟方法。该方法基于三角形非结构网格有限体积法数值离散框架，运动波和动力波控制方程均采用守恒形式，变量存在单元中心，采用的离散格式为Godunov格式。通过对网格单元求解器属性赋值来区分求解方法，属性为“0”时，求解运动波控制方程，属性为“1”时求解动力波控制方程。运动波和动力波的连续方程相同，所有网格单元的连续方程均采用Roe格式统一求解。当网格属性为“0”时，采用运动波动量方程直接求解下一时刻守恒变量；当网格属性为“1”时，采用Roe格式计算界面动量数值通量，通过动量通量、底坡项以及摩阻项求和来求解下一时刻的守恒变量。本方法可有效提高水动力模型在山区洪水过程模拟中的实际应用能力，弥补现有水文方法的不足。

主权项

一种运动波与动力波相结合的山区洪水过程数值模拟方法，其特征在于：采集山区流域地形数据、土地利用类型数据和降雨过程数据，且包括以下具体步骤：（1）网格离散及常规参数赋值：采用三角形非结构网格离散山区流域；采用地形数据对网格节点和网格型心进行高程插值，获取网格节点和网格型心处的高程数据；给每个网格设置糙率和初始条件；（2）求解器属性赋值：对每个网格进行求解器属性赋值，赋值“0”代表该网格将采用运动波方程求解；赋值“1”代表该网格将采用动力波方程求解；（3）确定时间步长：整套离散方法采用显格式进行，在计算中不同属性网格计算时采用统一的时间步长；（4）水力要素求解：水力要素变量定义在网格单元型心，整套离散方法采用基于非结构网格Godunov格式有限体积法进行；运动波和动力波连续方程一致，均采用Roe格式求解，即由t时刻的净雨强度值和网格单元的水力要素值计算t时刻网格单元降雨源项值和界面处的流量通量值，单元流量通量和源项求和后求出t+1时刻网格单元的水深h值；当网格属性为“0”时，求解运动波动量方程直接给出t+1时刻x，y方向守恒变量hu和hv值，当网格属性为“1”时，动量方程为动力波形式，由t时刻网格单元的水力要素值计算t时刻单元界面处的动量数值通量，求解地形源项和摩阻源项，单元三条边动量通量和源项求和后求出t+1时刻x，y方向的守恒变量hu和hv的值；各网格单元t+1时刻的hu和hv值除以t+1时刻水深h值，得到t+1时刻x，y方向的流速分量u，v值；（5）获取t+1时刻的净雨强度值，重复步骤（3）~（4），直至计算结束。