溃坝水流数值模拟.pdfVIP

摘要水 溃坝作为一种突发性强和破坏性大的灾难，对人民生命财产构成了极大的 威胁。因而探讨渍坝问题对于制定地区应急预案，加强区域防洪减灾是十分必 要的。溃坝问题计算是进行工程设计、预报预测和管理维护的基础，溃坝计算 的主要内容是算出溃坝坝址的流量和水位过程线，以及下游洪水演进过程中沿 程各处的流量、水位、流速、波前和洪峰到达的时间等。国内外学者对这类问 题做了大量研究，主要有理论研究、试验研究和数值研究，而由于问题的复杂 性，数值研究成为该类问题的主要研究手段。 近些年来，国内对溃坝问题做了大量数值研究取得丰富的成果。溃坝水流 一维数值模拟主要基于明渠水流运动的圣维南方程。由于它具有物理概念明确、 计算简单、计算结果相对可靠等优点，目前仍在实际工程中广泛应用，尤其对于长 河段、大尺度的模拟，但是其精度不足。二维数值模拟的研究发展迅速，二维模 型主要基于二维浅水方程，数值计算方法包括差分法、有限体积法和有限元法 等，二维模型以及一、二维相结合的统一数学模型也逐渐被应用到工程实践中。 由于浅水方程做了一定简化和假设，其应用也受到限制。与一、二维数学模型 相比，对溃坝水流三维数值模拟的研究相对较少，研究成果也远未达到实际应用 的要求，主要是由于三维模型存在模型复杂、计算量大和计算速度慢等缺点。由 于商业软件的出现，解决了相当一部分困难，为三维溃坝的研究提供了便利。 本文基于商业软件FLUENT分别对二维和三维溃坝水流进行数值研究，目 的在于在考察应用N．S方程数值模拟溃坝水流的可靠性的同时，探讨用物理实 验和理论分析不易获知的一些溃坝水流基本现象和规律。主要内容如下： ①采用有试验数据的二维模型模拟在层流和不同湍流模式下二维全溃坝洪 水下泄过程，并将得到的各个时段的流场状况与试验结果相比较分析，找出与 试验数据吻合较好的湍流模式。 ②采用与实验数据吻合较好的湍流模式模拟二维部分溃坝洪水下泄过程， 并分析得到的各个时段的流场状况。 ⑧采用有试验数据的三维模型模拟在选定的湍流模式下三维溃坝洪水下泄 过程，并将得到的各个时段的流场状况与试验结果进行比较分析。 }基金项目：奉文是水利部江湖治理与防洪重点实验室开放研究幕金资助项目(编号R2007．05) I 通过上述数值模拟和分析，得到以下结论： ①在对试验模型的二维全溃坝水流在不同的湍流模式下进行数值模拟中， 在选取五种湍流模型和层流中k—s模型得到的结果最好，洪水前峰随时间到达 的距离和试验结果吻合，自由液面的演化过程也比较好的吻合。 ②在采用k—g模型对3个尺寸相同、溃坝程度不同即全溃坝，顶部三分之 二溃溃坝和顶部三分之一溃溃坝的二维溃坝水流进行数值模拟中，得到了结论： 不同溃坝方式在所讨论的时间内都存在一段平稳状态；与全溃坝相比，部分溃 坝液面变化复杂，并且上游形成的波形也不同；在相同的时间内溃口的大小对 推进的距离有影响，特别是在溃坝初期溃口的大小对推进距离影响较大，但到 了后期三分之二溃坝水流fj{『锋与全溃坝水流前锋推进距离相差甚少，这说明在 剩余坝体低于一定的程度后，在后期对水流前锋的推进距离没有影响或者影响 很小。 ③在采用k—s模型对三维溃坝水流进行数值模拟中，在定性上k一占模型能 够比较好的模拟溃坝水流结构以及溃坝水流的传播过程，在定量上还是和试验 结果有所出入，原因有两方面：首先模拟的对象属于非定常、非线性的复杂水 波，甚至含有碎波现象，对很多因素的影响敏感；其次从试验数据可见自由表 面是急剧涨落不稳定的，如何在数值计算中平均化处理本身就是个问题。绕障 碍物的三维水流的数值计算结果充分再现了复杂的水面波形演化现象，如壅起、 回水，反射、干扰等，同时表明障碍物对溃坝水流的传播和水流的内部结构具 有显著的影响。 关键词： 溃坝，湍流模式，两相流，数值模拟 II Abstract Dam-breaka and destructive researchonthe is sudden catastrophe，SO strong dam-break