一种有限体积法在溃坝模拟中应用.pdfVIP

第三届黄河国际论坛论文集 355 一种有限体积法在溃坝模拟中的应用 向 波1 纪昌明2 罗庆松3 (1．武汉大学水资源与水电工程科学国家重点实验室； 2．华北电力大学能源与动力工程学院水资源与水利水电工程研究所； 3．中国卫星通信集团人力资源部) 摘要：非结构的三角形网格适应于复杂不规则的边界，在此基础上采用有限体积法离散浅水 波方程，并结合有限差分法建立了一种新的离散格式，使得界面通量计算达到二阶精度。通 过对典型算例的模拟计算，此格式在溃坝计算中能很好地捕捉溃坝波的前进，具有很强的模 拟间断的能力，模拟溃坝洪水波间断的形状和位置都能得到很满意的结果。 关键词：有限体积法溃坝二阶精度浅水波方程 大坝溃决后的溃坝波会对滞洪区和坝下游造成巨大的危害，所以以溃坝波 为研究对象的溃坝洪水波演进历来受到人们的极大关注。溃坝所造成的洪水波 是一类前进长波，非线性特性较强。溃坝计算是进行工程设计、预报预测、管理 维护的基础，计算的主要内容是通过一定的计算模式来描述滞洪区和下游洪水 波的传播速度、波形、波高、水位、流量等特征要素随时间变化的过程。国内外学 者对溃坝进行了大量的研究，早期通常是理论求解Sait—Venant方程，但是这种 方法有很大的局限性。由于问题的复杂性，二维溃坝水流方面尚未有有关解析 解的文献，但随着计算机与数值方法的发展，数值求解溃坝问题逐渐成为研究该 类问题的主要手段。 有限体积法是将计算区域划分为一系列不重复控制体，将待解的微分方程 对每一个控制体进行积分。有限体积法的离散格式有很多种，其各自的区别主 要是节点布置及对控制体流出的通量的计算方法不同。本文单元节点布置于网 格节点上，尝试用有限差分法计算界面上的通量，使得差分法到达二阶精度，并 通过计算对比，研究其模拟溃坝流动的可行性。 国家自然科学基金项目：洪灾异性风险的综合分析与评价理论及应用研究。 356 1 数学模型 二维浅水方程的分量形式可表达为 萼+笔+av，：0 (1a) 瓦OV+U¨x、TwJ、+嘉(等)=一旦Oyf、迸2 11一g硒。，一曲勘+∥(窘+等)(，c) 河道的倾斜效应项；戈、Y下标是求导的意思；(S厅，K)为x,y方向摩阻底坡。 二维浅水方程是一个对流一扩散方程，可以写成以下形式： G。+V·F=S (2) U U 扩 Uy J ^ = h G：似： ，g lVj Uy 俨 广=㈥ ^ h 0 r，o 1 一知ngh(Sz+Sox) g“2眇2 一芝1．g^2一gh(Szy+Soy) 2数值方法 采用有限体积法对方程(1)进行离散求解，为了适应复杂几何形状流场的 数值计算，网格采用任意三角形，节点为网格顶点。使用有限体积法时，把控制 体的厚度看成是1，有下面的关系： J G。+IV·F=I．S (3) F包含对流项(厂1)和黏性扩散项(，“)，S是源项。 2．1时间项离散 本文采用显示离散格式，则速度对时间的偏导数可