用户名: 密码: 验证码:
河道与渐溃堤坝耦联的水力数值模拟的研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
河堤溃决后的溃堤波会对圩区造成很大的危害,对溃堤波的研究是一个重要课题,研究其演进过程具有重要意义。本文将河道与堤坝下游区域视为一个整体,对于单一河道采用Preissmann四点隐式差分格式计算,对于溃口及圩区采用有限体积法,建立了一维河道与堤坝渐溃耦联的水力模型,也是一维、二维耦合的水流模型。在应用有限体积法无结构网格的基础上建立了高性能的TVD—MUSCL格式,根据堤坝溃口处的流量、水位,计算出沿程下游各点的流速、水位。算例验证表明:应用组合型TVD—MUSCL格式配合有限体积法是一种高解析度的数值方法,且TVD—MUSCL格式对于自动捕捉激波和抑制间断附近的数值解波动十分有效,它对间断具有高分辨率。本文建立逐渐溃堤的水力模型首次对二维逐渐溃堤波的流动问题进行了数值研究,揭示了逐渐溃堤过程中溃口处流场的变化情况,以及溃堤波的传播、绕射、反射及变形的复杂运动特征,能较好地模拟实际工程中的溃堤失事,模拟比较逼近于真实溃决过程,计算精度可满足工程要求。其研究成果能够为防汛部门的正确决策以及溃堤灾害性分析和防灾减灾提供科学的依据,具有较大的实际指导意义。
The dike-break long waves cause great damage to polder, and then it is a significant subject for studying the dike-break long wave. It is significant of research for the propagating of the dike-break wave. In the paper, river and the downstream area of dike are regarded as a whole system, and hydraulic model of coupling 1-D river and gradual dike-break is established, therein the Preissmann implicit difference scheme is applied to Main River, and the FVM (Finite Volume Method) is applied to the breach and polder. On the basis of FVM and unstructured grids, TVD ?MUSCL scheme is constructed with high performance. According to the discharge and water level of breach, the scheme is used to compute velocity and water level of center node of all elements. The results show that it is the numerical method of high analyzable degree for the hybrid TVD scheme combinated with FVM, and the scheme is not only sufficiently accurate and nonoscillatory, but also capable of treating automatically hydraulic jump. The hydrauli
    c model of gradual dike-break simulates numerically the 2-D flow of gradual dike-break waves. It shows the variation of flow field on breach, and the complicated flow characteristic of the propagation, diffraction, reflection and deformation of the dike-break waves. The hydraulic model can realistically simulate the burst process of dike, and computational accuracy is satisfied with the engineering demands. It is concluded that the model and methods have scientific decision reference and practical significance for flood warning service, and provide scientific basis for the analysis and prevention of disaster caused by dike destruction.
引文
[1] 曹楚生 水利水电工程设计 从大坝设计和风险分析看大坝安全 Vol.19,No.1 2000:1~5
    [2] 谢任之 溃坝水力学 山东科学技术出版社 1993.5
    [3] Nikolaos Katopoders Theoder strelkoff二维溃坝洪水波的计算 Journal of Hydraulics Division ASCE Vol.104 No Hy 9 1978.P.:1269~1287,杨立业译,西北水利科学研究所,1980
    [4] Wurbs R A. Dam-break flood wave models, Journal of Hydraulic Engineering, ASCE, 1987, 113(1): 29~46
    [5] 谭维炎 胡四—计算浅水动力学的新方向水科学进展 1992,3(4):310~318
    [6] 王嘉松 倪汉根 金生 二维溃坝问题的高分辨率数值模拟 上海交通大学学报,1999,33(10):1213~1216
    [7] 汪德爟 计算水力学理论与应用 河海大学出版社1989.9
    [8] 杨丽 河网与湖泊耦联的水力水质数值模拟的研究与应用 河海大学硕士学位论文2003.4
    [9] 陈锡林 新沂河东段水漫滩河道糙率计算讨论 水文2000年第20卷第1期
    [10] 袁世琼 天然河道的糙率计算 水电站设计Vol.13,No.1,1997(3):82~86
    [11] 孙文初 刘霞 李伦等 溃坝洪水流量计算方法浅析 广东水利水电 1999(3)
    [12] 谢任之 溃坝水力学山东科学技术出版社 1992.2
    [13] 孔晋豫 溃坝洪水数学模型的初步应用电力勘测No.3,1995.9:53~57
    [14] 赵亚波 张帜 溃坝洪水过程线探讨东北水利水电No.10,1996:23~24
    [15] 郭洪巍 吴葱葱 逐渐溃坝和瞬时溃坝的模拟研究与比较东北水利水电Vo.18,No.2,2000
    [16] 谭维炎 胡四—用TVD格式预测溃坝涌波的演进 水力学报 1989(7):1~11
    [17] 谭维炎计算浅水动力学—有限体积法的应用清华大学出版社 1998.9
    [18] Osher, S. and Solomone, E, Upwind Difference Scheme for Hyperbolic Systems of Conservation Laws, mathematical computation, vol.38, 1982:339~374
    [19] Zhao, D.H., H.W Shen, G.Q.Tabios Ⅲ, J.S.Lai and W.Y.Tan, " Finite-Volume Two-Dimensional Unsteady-Flow Model for River Basins, " Journal of Hydraulic Engineer, ASCE, Bol.120, No.7, 1994:863~883
    [20] Tan, Weiyan, Shallow Water Hydrodynamics, Elsevier, Amsterdam 1992
    [21] 陈树群 赵棣华 沈学汶等 二维有限体积模式分析溃坝水流中国土木水利工程学刊Vol.9,No.3,1997:439~448
    [22] Hirsch, Dharles, Numerical Computation of Internal and External Flow, Vol.1, Wiley-Interscience Publication, 1988:241~244
    [23] 谭维炎 胡四一等 一维非恒定流计算程序包:MYBC,水力学报,1982(1)
    [24] 王承绕 王正华 杨晓辉等 计算流体力学及其并行算法国防科技大学出版社2000.2
    [25] 张涤明 蔡崇喜 章克本等 计算流体力学 中山大学出版社1991.5
    [26] 傅德薰 流体力学数值模拟 国防工业出版社1993.11
    [27] 朱幼兰 钟锡昌 陈炳木等 初边值问题差分方法及绕流科学出版社,1980
    [28] Moretti C Three-dimensional supersonic steady flows with any number of imbedded shocks. AIAA Paper No.74-10, 1974
    
    
    [29] Lax P D and Wendroff B. Difference schemes for hyperbolic equations with high order of accuracy. Comm. Pure AppI. Math., 1964 (17) : 381~398
    [30] 傅德薰 熵条件及差分方程的数值解 第二届全国计算流体力学方法会议论文集,1983
    [31] 傅光明 傅德薰 超音速绕流计算中熵条件的几个问题 第一届全国计算流体力学方法会议论文集,1981
    [32] Le Roux AY. A numerical conception of entropy for quasilinear equations. Math comp. 1979(31): 848~872
    [33] Steger J and Warming R F. Flux vector splitting of the inviscid gas-dynamics equations with application to finite methods. NASA TMD-78605, 1979
    [34] 傅德薰 马延文 气动计算中一个简单有效的差分格式 空气动力学学报,1985(2)
    [35] 马延文 傅德薰 数值计算可压缩Navier-Stokes方程的一个新的系数矩阵分裂法 计算物理,1985(2)
    [36] Garten A. On a class of high resolution total-variation-stable finite-difference schemes. NYU Report, Oct.1982
    [37] Roe P L. Some contributions to the modeling of discontinuous flows. AMS-SIAM Seminar, Sam Diego, June 1983
    [38] Var Leer B. Towards the ultimate conservative difference scheme. Journal of Computational Physics.1974(14): 363~389
    [39] Osher S. Shock modeling in transonic and supersonic flow recent advances in numerical methods in fluids. Vol.4, W.G.Habashi Ed
    [40] Yee H C, Warming R F and Harten A. Implicit total variation diminishing (TVD) schemes for steady-state calculations. AIAA Paper 83-1902, 1983
    [41] Harten A and Hyman J M. A self-adjusting grid for the computation of weak solutions of hyperbolic conservation laws. Los Alamos Nat. Lab. Report LA. 9105, 1981
    [42] Yee H C. Construction of explicit and implicit symmetric TVD schemes and their applications. Journal of Comp. Phys., 1987(68): 151~179
    [43] Harten A. On high-order accurate interpolation for non-oscillatory shock-capturing schemes. The IMA Volumes in Mathematics and its Applications .Springer-Verlag 1986 (2): 71~106
    [44] 谭维炎 浅水动力学的回顾和当代前言问题水科学进展 Vol.18,No.2 Mar.,2001:97~105
    [45] Desiseri D and Jameson A. Multigrid solution of the Euler equations on unstructured and adaptive meshes. AIAA Paper 87-0353, 1987
    [46] Pan D and Cheng J. Incompressible flow solution on unstructured triangular meshes. Int. J. Numer. Methods Fluids. 19939(16): 1079~1098
    [47] Zhao D H et al. Finite-volume two-dimensional unsteady-flow model for river basins. ASCE J. Hydraul. Eng. 1994 (120): 864~883
    [48] Anastasiou K and Chan C T. Solution of the 2D shallow water equations using the finite volume method on unstructured triangular meshes. Int. J. Numer. Methods Fluids. 1997(24): 1225~1245
    [49] 汪继文 刘儒勋 间断解问题的有限体积法计算物理 Vol.18,No.2 Mar.,2001:97~105
    [50] 赵棣华 戚晨 庾维德等 平面二维水流一水质有限体积法及黎曼近似解模型 水科学进展 Vol.11,No.4,Dec.,2000:368~374
    
    
    [51] Spekreijse, S.P., Multigrid solution of steady Euler equations. CWI Trac 46 Amsterdam(1988).
    [52] 谭维炎 胡四一 二维浅水流动的一种普适的高性能格式水科学进展 Vol.2,No.3,Sep.,2000
    [53] 马洪明 感潮河口数值模拟及动网格研究 河海大学硕士学位论文 2002.4
    [54] 杨丽 河网与湖泊耦联的水力水质数值模拟的研究与应用 河海大学硕士学位论文 2003.4
    [55] zhao,D.H., G.Q.Tabios Ⅲ and Shen, H.W., RBFVM-2Dmodel, River Basin Two-Dimensional Flow Model Using Finite Volume Methods, Program Documentation, University of California, Berkely, California USA(1995)
    [56] 谭维炎 胡四一 浅水流动计算中一阶有限体积法Osher格式的实现 水利部南京水文水资源研究所 水科学进展 Vo.5,No.4,1994.12
    [57] 谭维炎 胡四一 无结构网格上二维浅水流动的数值模拟水科学进展Vo.6,No.1,Mar.,1995:1~9
    [58] Shu C W. Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws [A] , Lecture Notes Math. 1697, Springer, 1998:325~432
    [59] 水鸿寿 一维流体力学差分方法 国防工业出版社1998.2
    [60] Roe P L. Approximate Riemann solvers, parameter vectors and difference schemes[J] . J of Comput Phys, Vol.43, 1981:357~372
    [61] Van Leer B, Thomas J L, Roe P L, Newasome T W. A comparison of numerical flux formulas for the Euler and Navier-Stokes equations. AIAA 87-1104
    [62] Harten A, Lax P D, van Leer B, On upstream differencing and Godunov-type schemes for hyperbolic conservation laws [J] , SIAM Rev, Vol .25, 1983:35~61
    [63] Toro E E Riemann problems and the WAF method for solving the two-dimensional shallow water equations [J] . Phil Trans R Soc Lond, 1992, A338:43~68
    [64] Godunov S K. A finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics. Sov. Math. Sb. 47: 271~306, 1959
    [65] Thomas H L, Taylor S L, Anderson W K. Navier-Stokes computations of vortical flows over low aspect ratio wings. AIAA 87-0330
    [66] Anastasiou K, Chan C T. Solution of the 2D shallow water equations using the finite volume method on unstructured triangular meshes [J] . Inter J for Numer Methods in Fluid, Vol.24, 1997: 1225~1245
    [67] Harten A, Osher S. Uniformly high order accurate essential non-oscillatory schemes Ⅰ.SIAM J. on Num.Anal, Vol.24, 1987:279~309
    [68] Hagen A, Engquist B, Osher S, Chakravarthy S R. Uniformly high order accurate essentially non-oscillatory schemes Ⅲ. J. Comp. Phys., Vol.71, 1987:231~303
    [69] Chakravarthy S R, Osher S. A new class of high accuracy TVD schemes for hyperbolic conservation laws. AIAA 850363
    [70] 刘儒勋 王志峰等 数值模拟方法和运动界面追踪 中国科技大学出版社 2001.10
    [71] 朱自强等 应用计算流体力学 北京航空航天大学出版社 1998.8
    [72] Hu K, Mingham C G, Causon D M. A bore-capturing finite volume method for open-channel flows [J] . Inter J Numer Methods in Fluids, Vol.28, 1998:1421~1261
    [73] Harten, A., A high resolution scheme for the computation of weak solutions of hyperbolic conservation laws. New York University Report, Oct., 1981, and J. Comp. Phys., Vol.49,
    
    1983: 357~393
    [74] Harten A. High resolution schemes for hyperbolic conservation laws [J] . Comput. Phys.,Vol.49, 1983:357~393
    [75] Yee H C. Upwind and symmetric shock capturing schemes. NASA TM 8946, 1987
    [76] 陆夕云 庄礼贤 童秉纲等 对称TVD差分格式的构造及其应用Vol..11,No.1 Mar.,1994
    [77] Yee, H.C., Linearized form of implicit TVD schemes for the multidimensional Euler and Navier-Stokes equations. Comp. & Math. With Appl., Vol.12A, 1986:413~432
    [78] Sweby P K. High resolution schemes using flux limiters for hyperbolic conservation laws. SIAMJ. Num. Anal. Vol.21, 1984:995~1011
    [79] Yee H C. Construction of explicit and implicit symmetric TVD schemes and their applications. J. of Comp.Phys., Vol.68, 1987:151~179
    [80] Van Leer B. Towards the ultimate conservative difference scheme V: A second order sequel to Godunov's method. J. Of Comp. Phys., Vol.32, 1979:101~136
    [81] Chakravarthy S R, Osher S. A new class of high accuracy TVD schemes for hyperbolic conservation laws. AIAA 850363
    [82] 张涵信 无波动、无自由参数的耗散差分格式 空气动力学报 Vol.6,1988:143~164
    [83] Harten A. The artificial compression method for computation of shocks and contact discontinuities I single conservation laws. Comm. Pure Appl. Math. Vol.30, 1977:61~638
    [84] Yee H C. Numerical experiments with a symmetric high-resolution shock-capturing scheme. Lectures Notes in Physics, Vol.246, Springer-Verlag,1987:677~683
    [85] Yee H C. Construction of explicit and implicit symmetric TVD schemes and their applications. J. of Comp.Phys., Vol.68, 1987:151~179
    [86] Roe P L. Generalized formulation of TVD Lax-Wendroff schemes. ICASE Report No.84-53,1984
    [87] Van Leer B, Thomas J L, Roe P L, Newasome T W. A comparison of numerical flux formulas for the Euler and Navier-Stokes equations. AIAA 87-1104
    [88] 王志刚 实际地形溃坝水流的数值模拟 河海大学硕士学位论文2003.4
    [89] 刘儒勋 王志峰等 数值模拟方法和运动界面追踪 中国科技大学出版社2001.10
    [90] Ritter M B. Die fortpflanzung der wasserwellen, Zeitschrift des Vereines deutscher Ingernieure, 1892, 36 (33): 947~954.
    [91] 王如云 浅水波方程的TVD有限差分数值模拟海洋与湖沼 1991,22(2)

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700