目录 1.前言
2.控制方程
3.有限体积TVD格式
4.边界条件
5.数值预报
6.结语 参考资料 [1] Harten A. High resolution schemes for hyperbolic conservation laws, J. Comp. Phys., 1983, 49: 357-393.
[2] 胡四一、谭维炎，用TVD格式预测溃坝洪水波的演进，水利学报，1989，(7)：1－11.
[3] 陶建华, 张卫东. 总变差不增格式计算一维、二维溃坝波. 天津大学学报, 1993, (1): 7~15.
[4] Yang J Y, Hsu C A and Chang S H. Computations of free surface flows, Part 1: 1D dam-break flow. J. Hydr. Res. , 1993, 31(1): 19-34.
[5] 王嘉松, 倪汉根, 金生, 李鉴初，用TVD显隐格式模拟溃坝波的演进与反射，水利学报，1998，（5）：7-11.
[6] 王嘉松, 倪汉根, 金生，二维溃坝波传播和绕流特性的高精度数值模拟，水利学报，1998，（10）：1-6.
[7] 王嘉松, 倪汉根, 金生，二维溃坝问题的高分辨率数值模拟，上海交通大学学报，1999，（10）：1213－1216.
[8] Wang Jia-song and Ni Han-gen. A high resolution finite-volume method for solving the shallow water equations. J. Hydrodynamics, Ser. B., 2000, (1): 35-41.
[9] Wang Jia-song, Ni Han-gen and He You-sheng. Finite-difference TVD scheme for computation of dam-break problems. J. Hydr. Eng., ASCE, 2000, 126(4): 253-262.
[10] 谭维炎, 胡四一，二维浅水流动的一种普适的高性能格式—有限体积Osher格式，水科学进展，1991，(3)：154-161.
[11] Zhao D H, Shen H W, Lai J S and Tabios Ⅲ G Q. Approximate Riemann solvers in FVM for 2D hydraulic shock wave modeling, J. Hydr. Eng., ASCE, 1996, 122: 692-702.
[12] Anastansiou K and Chan C T. Solution of the 2D shallow water equations using the finite volume method on unstructured triangular meshes. Inter. J. Numer. Meth. Fluids, 1997, 24: 1225-1245.
[13] Mingham C G and Causon D M. High-resolution finite-volume method for shallow water flows. J. Hydr. Eng., ASCE, 1998, 124(6): 605~614.
[14] Bell S W, Elliot R C and Chaudhry M A. Experimental

摘 要 基于TVD格式的有限差分算法思想，考虑节点与单元的对应关系，对任意四边形单元各边的通量采用优选限量函数的组合型TVD格式进行插值，在时间上采用两步Runge-Kutta法离散，建立了守恒型浅水方程有限体积的高分辨率计算模型。首次针对180°强弯曲河道、90°双支以及45°三支分叉渠道，数值预报了溃坝波的演进过程，揭示了溃坝波在弯曲河道中内外两岸速度与水位变化，在分叉河道中自动进行流量与动量再分配，在叉点角区形成旋涡、壅高等复杂运动特征，同时反映了湿底与干底的影响。 关键词 浅水方程，有限体积，TVD格式，弯曲渠道，分叉渠道，溃坝波 　 1．前言 许多水利和环境工程问题，都需要进行浅水流动分析，数值计算逐渐成为分析的重要手段。虽然二维浅水流动一般问题的计算比较成熟，但对于有间断、带自由面、含复杂边界等复杂浅水问题的研究，由于对数值计算方法要求较高仍处于发展之中。对溃坝问题的研究，一直是引人关注的重要课题，具有重要的学术价值和工程应用背景。 过去对溃坝问题的研究，以差分方法为主，所用的格式有特征线法、显式和隐式差分法以及最近占据重要地位的近似黎曼解法等，TVD（Total Variation Diminishing）格式[1]在该类问题研究中也逐渐发挥其独特的作用, 如用来计算一维溃坝波的传播、反射[2~5]和二维溃坝波的传播、反射与绕射[6~9]。对于复杂边界，通常采用坐标变换，将不规则计算域转换到规则计算域，但同时方程也复杂化。近年来，在有限体积离散基础上的近似黎曼解法，用于计算浅水流动（包括溃坝波）问题，取得了良好效果[10~13]。但另一类高分辨率格式—TVD格式用在有限体积方法中尚不多见。