维溃坝洪水波的数值模拟一时空守恒法.PDFVIP

2005年5月 重 庆 大 学 学报 Mav 2005 第28卷第5期 Journal of Chongqing University Vo1．28 No．5 文章编号：1000—582X(2005)05—0136—03 一 维溃坝洪水波的数值模拟一时空守恒法 张永祥，陈景秋，韦春 霞 (重庆大学 工程力学系，重庆 400030) 摘 要：应用改进的时空守恒法对一维溃坝洪水波进行了数值模拟，从计算结果与分析解比较可 知，该法能够很好地模拟溃坝洪水波，甚至对于下游为千河的情况，依然能得到满意的结果。而且，此方 法精度高，有较强的渍坝涌波捕捉能力，可展示渍坝波的运动特性，为溃坝引起的灾害预测和分析及防 灾减灾提供一种有效的数值模拟。 关键词：溃坝；圣维南方程；时空守恒律；涌波 中图分类号：0352 文献标识码：A 水坝的安全问题构成了对人类生命和财产的潜在 1 一维浅水波方程 威胁川。在数学表达上，重力引起的浅水运动类似于 可压缩流体的运动，求解含激波的空气动力学双曲波 St．Venant一维浅水波方程守恒形式为 问题 的方法可用来求解含涌波的浅水运动问题。 鲁+ Ox_gm，m_l’2．af … ’ 近年来 MacCormack，Bean—Warming，Godunov， TVD J， Adaptive Grid 等格式都取得了一定的效果， 其中： 但并不十分理想，一般差分格式在涌波的捕捉上往往 “1=h，11．2=hu， 不能同时满足虚假振荡的抑制和足够高的精度。近年 = hu=11．2， 来S C Chang 】提出了一种崭新的数值方法，它与传统 (2) 的数值方法相比具有许多独特的优点：1)其构造思想 +号^ u_ 3z+g2 _ 简单，物理意义清晰，格式通用性好，它不用求解黎曼 q1=O,q2=gh(So一 )． 问题，也不用求解特征问题。2)它把时间与空间完全 其中 = ， 为Manning摩阻系数，R为水力学半径， 统一同等对待，并从守恒律积分行方程出发，通过设立 守恒元和解元，使局部和全局都严格保证其物理意义 当河道为矩形时约等于^，u为流速，S。为底坡，S，为摩阻 上的守恒律。3)该方法把流场变量和其对空间导数 坡降，g为重力加速度， 为沿河底的距离，取下游方向为 均作为作为变量同时求解，这样与传统的差分格式相 正，t为时间。 比，在相同的格式点数的情况下，其格式精度达到更 令 = ， =t代表Euclidean空间 中的2个坐 高，还更便于边界条件的精确履行。但原始方法在一 标，利用高斯散度定理可知，微分方程(1)可以写成如下 维情况要设立2个守恒元，不易推广到更高维，张增 的积分形式： 产 对此方法进行了改进，使得方法简化，且保留了 牛h ·ds=牵g dV，m=1，2． (3) 原方法的优点，并成功的模拟了空气动力学中的激波 s 现象，得到的精度更高。为此，笔者应用改进的时空守 其中S(V)是 空间中任一区域 的边界，ds=(山，山)， 恒法对一维溃坝洪水波进行了数值模拟。