河流数学模型理论与实践案例.docx

河流数学模型理论与实践案例一、案例案例性质：一维非恒定水流模型案例背景：（1）以长江为基本背景，构建一条河长100km的无支流入汇河段，并给出断面几何形态和其他物理参数。（2）给定入流条件Q符合e指数分布律，平均值取Q0=4000m3/s；并根据断面形态和河段坡度给出下游控制水位。结果要求：（1）依据Preissmann格式导出Holl-Cunge水动力学流量演算方法；（2）给出Holl-Cunge水动力学流量演算数值方法的三性分析，并给出该方法的收敛特性分析R1、R2图；（3）利用马斯京根法和Holl-Cunge方法，在上述案例背景上进行流量演算，给出数值演算结果并对两种方法的计算结果进行对比分析；（4）考虑基本方程试证明其为对流扩散方程，理论阐明采用Holl-Cunge流量演算方法模拟该对流扩散方程的技术方法。二、Holl-Cunge水动力学流量演算方法一维水流数学模型基于以下圣维南方程组：水流连续方程：（1）水流运动方程：（2）式中，为过水断面面积，为流量，为水位，为重力加速度，为能坡，为流程方向，为时间。假定在某一确定断面上，流量和水位成一一对应的单值关系（3）这里为过水断面面积，为水深。若假定式（3）成立，则反函数存在，且也为单值函数，因此有（4）将式（4）代入连续方程式（1）有（5）这是一阶的双曲型方程式，因为变量Q的函数，故是拟线性的。在特征线上，有（6）（7）结合式（5）有（8）由此可写出方程（5）的另一种形式为（9）在满足式（6）的条件下，方程（9）的另一等价形式及相应的解为或（10）这个解表明流量沿特征线族上为一常数，并不允许流量过程在任意时刻任何断面上发生任何变形输移。从这个结果可得出一个结论，流量水位成单值关系的假设，实质上描述的是流量的非线性对流输移，入流过程线自始至终将保持原来的形状。由于方程（9）或式（5）是非线性的，故很难得到解析解，可考虑求其数值解。采用四点偏心Preissmann格式逼近方程（9）中关于时空偏微商项，有（11）这里和为权重因子，，，和分别为节点函数值，式中的暂取为常数。其网格布置如下图1所示。图1 差分网格示意图令，并取，则从方程（11）解出为（12）式中（13）式（12）即为Holl-Cunge水动力学流量演算方法。在点（j，n）上按泰勒级数展开方程（11）中各项节点函数值，并经整理后得（14）略去的高阶项得（15）式中（16）很明显，当和取任意值时，数值格式（12）实质上以二阶近似逼近于一个对流扩散方程（15），而不是纯对流方程（9）；只有当时，才有，差分格式（12）才完全逼近微分方程（9）。考虑连续方程式（1）和忽略时空惯性项的运动方程式（2），有（17）（18）式中，为底坡，，为流量系数，，为谢才系数。当河宽沿程不变时，有，则（19）为便于分析，对式（17）关于x求偏微商（20）再次对方程（18）关于t求偏微商得（21）将方程（20）和方程（21）相加得（22）对K关于t求偏微商并使用方程（17）有（23）将式（23）代入式（22）化成一维线性抛物型方程（24）要使差分格式（12）逼近物理扩散方程（24），则必须满足（25）和（26）由此看出，忽略运动方程时空惯性项的圣维南方程组，在假定流量系数K是水深的单值函数的条件下，描述了沿程坦化的洪水扩散波，而不是运动波，且按式（26）选取数值计算参数，和，时可由差分方程（12）逼近。三、Holl-Cunge方法的三性分析在用差分方程代替微分方程问题，将求解微分方程化为求代数方程的离散解时，必须进一步对差分方程的收敛性、稳定性和相容性进行研究。（1）相容性将一个微分方程用差分格式化为相应差分方程，当步长与趋近于零时，这个微分方程应当收敛于原微分方程，也就是说，相应的差分方程和微分方程之间的截断误差在任一时刻任一网格点上均应趋近于零，这样的差分方程和微分方程才是相容的。对于上节得到的差分方程（11），假定的结点值是微分方程（9）的解（27）将，和在结点（j，n）处展开为泰勒级数有（28）（29）（30）并注意到（31）将式（28）~（30）代入差分方程（11）有（32）式中，因在点（j，n）处满足式（27），即（33）则有截断误差（34）由此可看出，用差分方程（11）逼近微分方程（9）的截断误差的阶是：①当和时，差分格式具有一阶精度；②当和时，具有二阶精度；③当和时，格式具有高阶精度。当与趋近于零时有（35）即当步长与趋近于零时，截断误差趋近于零。可见，差分格式（11）与微分方程（9）是相容的。（2）稳定性差分法计算中所产生的误差（舍入误差，参数误差等）随时间衰减或不增大，则称离散格式是稳定的，反之，则是不稳定的。分析差分方程稳定性有不同的方法，如矩阵方法，谐波分析法等。由于偏微分方程（9）是线性的，其解可用Fourier级数表