CN101017516A - 水工水力学的水流精细模拟方法 - Google Patents
水工水力学的水流精细模拟方法 Download PDFInfo
- Publication number
- CN101017516A CN101017516A CNA2007100515686A CN200710051568A CN101017516A CN 101017516 A CN101017516 A CN 101017516A CN A2007100515686 A CNA2007100515686 A CN A2007100515686A CN 200710051568 A CN200710051568 A CN 200710051568A CN 101017516 A CN101017516 A CN 101017516A
- Authority
- CN
- China
- Prior art keywords
- hydraulic
- flow
- turbulent
- water conservancy
- simulation method
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种水工水力学的水流精细模拟方法,属于水利工程技术领域。本发明研究并发展了水利水电工程紊流数值模拟中适合多类流动问题的紊流模型及模拟方法。模拟获得了水垫塘复杂紊流场水力特性和消能机理及孔板泄洪洞、洞塞泄洪洞、竖井旋流式泄洪洞和龙抬头式泄洪洞的流场细观特征,形成了较为完整的泄洪洞消能体系的理论;建立了一套复式游荡型河道漫滩水流计算方法及洪水演进数值计算模型。
Description
技术领域
本发明涉及一种水利水电工程紊流数值模拟,适用于多类流动问题的紊流模型及模拟方法。
背景技术
泄洪消能问题、河道水力学问题是水利水电工程十分关注的研究内容。以往对这类问题的研究主要依赖于实体模型试验。随着紊流数值模拟技术的发展,进行流动的紊流数值模拟研究成为可能,它有助于为水利水电工程的相关问题研究和设计开辟一条新途径。高坝水垫塘和高水头泄洪洞等的安全性问题迄今仍是水工水力学的难题之一,对此,流场精细模拟可以发挥十分重要的作用。水利水电工程泄洪研究直接关系到工程的安全,而河道中洪水的研究则关系到自然灾害的防治,后者同样需要深入地水流精细模拟研究。
发明内容
本发明的目的是要提供一种水工水力学的水流精细模拟方法。通过该方法能够对水垫塘复杂紊流场水力特性和消能机理流场精细模拟;能够对孔板泄洪洞、洞塞泄洪洞、竖井旋流式泄洪洞和龙抬头式泄洪洞的流场精细模拟;能够对复式游荡型河道漫滩水流计算方法及洪水演进数值计算模型。
本发明的目的是这样实现的:一种水工水力学的水流精细模拟方法,所述的模拟方法包括了对水垫塘复杂紊流场水力特性和消能机理及孔板泄洪洞、洞塞泄洪洞、竖井旋流式泄洪洞和龙抬头式泄洪洞的流场细观特征模拟,以及一套复式游荡型河道漫滩水流计算方法及洪水演进数值计算模型。本发明建立了水利水电工程水流精细模拟的大型的数值实验室软件包,具有模拟水利水电工程复杂水流的能力,数值计算结果通过动态显示技术以形象的图形输出,从而可以人机对话的方式灵活地修改模型或原型的形状和参数,方便地进行系列数值试验或不同的方案比较,为方案选择和设计提供数据。本发明通过用户自定义函数UDF实现了曳力系数的自定义,通过改变曳力系数改善气泡竖向速度和掺气浓度分布。本发明采用三维紊流模型和VOF方法,对涡室和竖井的旋流涡心形态、通气量、空化数和消能率等进行了定量的描述和分析。本发明应用两相流理论对带有通气槽的“龙抬头”式泄洪洞进行了数值模拟,并对高速水流边壁凹体的空化特性进行了研究。本发明通过对洞塞式泄洪洞的紊流数值模拟,得出了洞塞消能的紊动能和紊动能耗散率的分布规律,详细地反映出洞塞的消能情况,得到压力最低点。本发明采用紊流数学模型成功地进行了水垫塘复杂紊流场水力特性、消能机理和设计优化的研究,获得了水垫塘三元典型流态概化图、水垫塘水体消能分区等。本发明运用紊流拟序理论对多股多层水平淹没射流水垫塘内复杂三元紊流结构进行了深入分析。
本发明所提供的水工水力学的水流精细模拟方法,包括水垫塘复杂紊流场水力特性和消能机理及孔板泄洪洞、洞塞泄洪洞、竖井旋流式泄洪洞和龙抬头式泄洪洞的流场细观特征,形成了较为完整的泄洪洞消能体系的理论;建立了一套复式游荡型河道漫滩水流计算方法及洪水演进数值计算模型。
具体实施方式
(1)坝下消能工水流特性的数值实验室软件包开发:水工水力学中涉及到的水流现象几乎都是紊流,为了确定水工建筑物的合理布置、掌握坝下游的流态、分析河床的冲刷、研究坝下水流的消能防冲措施,需对紊流进行专门研究。以往水工建筑物设计基本上以模型试验为依据,不仅费时耗物,且受比尺效应的影响,不能完全反映真实的紊流流动情况。相比之下,数值模拟具有花费少、适应能力强、提供详细的流场资料、便于方案比较等优点。人们愈来愈重视运用数值模拟手段预测各种复杂的水流现象及流场的内部结构。国际上一些享有盛名的研究所由过去的模型试验为主要研究手段转化为数值模拟为主要研究手段。课题组结合重大工程课题和多项自然基金课题,建立了水利水电工程水流精细模拟的大型的数值实验室软件包,具有模拟水利水电工程复杂水流的能力,数值计算结果通过动态显示技术以形象的图形输出,从而可以人机对话的方式灵活地修改模型或原型的形状和参数,方便地进行系列数值试验或不同的方案比较,为方案选择和设计提供足够精细和准确的数据,有助于数学模型逐步地或部分地取代费时耗物的物理模型试验。
(2)高紊动掺气水流数值模拟:经过系列的二、三维分析、计算,研究了水工中一类常见的高速掺气水流及掺气减蚀问题。特别是针对高速水流的高紊动、非定常,存在复杂自由表面和两相流动的特性,提出了一套新的数值计算方法。对比了单颗粒直径气泡与多颗粒直径气泡群对水气二相流模拟的影响。群颗粒气泡的模拟较之单粒径气泡的模拟更为符合实际,结果也更为合理;通过用户自定义函数UDF实现了曳力系数的自定义,通过改变曳力系数改善气泡竖向速度和掺气浓度分布,数值试验发现掺气浓度分布曲线对曳力系数并不敏感;数值计算结果显示小直径气泡在底面附近的浓度相对于其在水体内其他部位的浓度为大。研究结果表明,这里提出的新的数值计算方法是合理可行的。
(3)竖井旋流式泄洪洞数值模拟:采用三维紊流模型和VOF方法,对竖井旋流式泄洪洞的水流特性进行了详细的研究,获得了不同运行工况下泄洪洞内水流流速场和压力场的系列计算结果,对涡室和竖井的旋流涡心形态、通气量、空化数和消能率等进行了定量的描述和分析。经与模型试验成果对比,计算结果合理。该成果在竖井旋流的理论研究和应用上有所创新,为高速水流体型研究和设计打下了很好的基础。
(4)龙抬头式泄洪洞数值模拟:首次应用两相流理论对带有通气槽的“龙抬头”式泄洪洞进行了数值模拟,并对高速水流边壁凹体的空化特性进行了研究,计算给出了沿洞身水面线和断面的掺气浓度分布,为导流洞改建利用提供了技术及方法。
(5)洞塞式泄洪洞的紊流数值模拟:通过对洞塞式泄洪洞的紊流数值模拟,得出了洞塞消能的紊动能和紊动能耗散率的分布规律,详细地反映出洞塞的消能情况,特别是细致地捕捉了压力最低点,弥补了物理模型试验的不足。研究表明,导流洞改建高水头泄洪洞采用洞塞消能方案是安全可行的。
此外,通过对多级洞塞消能工的布置优化和三维非轴对称洞塞的大量计算,给出了多级洞塞在不同面积比和不同布置条件下的各级洞塞的消能率曲线,为工程设计提供了参考。
(6)水垫塘复杂紊流场水力特性、消能机理和设计优化的研究:采用紊流数学模型成功地进行了水垫塘复杂紊流场水力特性、消能机理和设计优化的研究,获得了水垫塘三元典型流态概化图、水垫塘水体消能分区等。以小湾水电站为例,对小湾水电站采用非平底型水垫塘后,下泄水流对水垫塘的影响进行了模拟研究,了解不同方案下水垫塘内水流流态的变化,以及水垫塘底板处水流各特征参数的变化,判断各方案在工程上的可行性。模拟研究成果先后应用于多项大型水电工程水垫塘的的水力学研究。其中三元直角坐标k-ε模型应用于拉西瓦水垫塘的水力学研究;流态特征、消能机理以及多股淹没射流消能优化的成果应用于小湾水垫塘的水力学研究;三元贴体曲线坐标的数学模型应用于溪洛渡水垫塘的水力学研究。成果也可推广应用于其它类似问题的数值模拟研究。
(7)水平多股淹没射流(跌坎——底流)新型消能方式的数值模拟研究:运用紊流拟序理论对多股多层水平淹没射流水垫塘内复杂三元紊流结构进行了深入分析。采用VOF法和RNG k-ε紊流模型,结合向家坝水电站,对多股多层水平淹没射流水垫塘流动结构进行了数值模拟和分析。结果表明,该消能型式对下游的适应性较强,消能率较高;物理试验结果表明,该消能型式在保持流态稳定的同时,能够显著降低临底流速,减小出坎水流波动,对下游的适应性较强,消能率较高;数值模拟结果还表明,水垫塘跌坎下游出现横轴漩涡和立轴漩涡等多种漩涡结构,流动结构复杂。探明了流动结构中立轴漩涡的形成机理,分析了漩涡结构对建筑物的影响,提出了改善水流结构,增强结构稳定的工程措施。研究表明,它是一种流态稳定、适应性强、低雾化、高消能率的新型消能方式,能够显著减轻工程建设给环境带来的不利影响,易于工程应用,从而丰富了现代泄洪消能新理论,为实际工程设计提供了科学依据。
Claims (8)
1.一种水工水力学的水流精细模拟方法,其特征在于:所述的模拟方法包括了对水垫塘复杂紊流场水力特性和消能机理及孔板泄洪洞、洞塞泄洪洞、竖井旋流式泄洪洞和龙抬头式泄洪洞的流场细观特征模拟,以及一套复式游荡型河道漫滩水流计算方法及洪水演进数值计算模型。
2.如权利要求1所述的水工水力学的水流精细模拟方法,其特征在于:建立了水利水电工程水流精细模拟的大型的数值实验室软件包,具有模拟水利水电工程复杂水流的能力,数值计算结果通过动态显示技术以形象的图形输出,从而可以人机对话的方式灵活地修改模型或原型的形状和参数,方便地进行系列数值试验或不同的方案比较,为方案选择和设计提供数据。
3.如权利要求1所述的水工水力学的水流精细模拟方法,其特征在于通过用户自定义函数UDF实现了曳力系数的自定义,通过改变曳力系数改善气泡竖向速度和掺气浓度分布。
4.如权利要求1所述的水工水力学的水流精细模拟方法,其特征在于采用三维紊流模型和VOF方法,对涡室和竖井的旋流涡心形态、通气量、空化数和消能率等进行了定量的描述和分析。
5.如权利要求1所述的水工水力学的水流精细模拟方法,其特征在于应用两相流理论对带有通气槽的“龙抬头”式泄洪洞进行了数值模拟,并对高速水流边壁凹体的空化特性进行了研究。
6.如权利要求1所述的水工水力学的水流精细模拟方法,其特征在于通过对洞塞式泄洪洞的紊流数值模拟,得出了洞塞消能的紊动能和紊动能耗散率的分布规律,详细地反映出洞塞的消能情况,得到压力最低点。
7.如权利要求1所述的水工水力学的水流精细模拟方法,其特征在于采用紊流数学模型成功地进行了水垫塘复杂紊流场水力特性、消能机理和设计优化的研究,获得了水垫塘三元典型流态概化图、水垫塘水体消能分区等。
8.如权利要求1所述的水工水力学的水流精细模拟方法,其特征在于运用紊流拟序理论对多股多层水平淹没射流水垫塘内复杂三元紊流结构进行了深入分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2007100515686A CN101017516A (zh) | 2007-02-12 | 2007-02-12 | 水工水力学的水流精细模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2007100515686A CN101017516A (zh) | 2007-02-12 | 2007-02-12 | 水工水力学的水流精细模拟方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101017516A true CN101017516A (zh) | 2007-08-15 |
Family
ID=38726521
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2007100515686A Pending CN101017516A (zh) | 2007-02-12 | 2007-02-12 | 水工水力学的水流精细模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101017516A (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101871204A (zh) * | 2010-06-28 | 2010-10-27 | 四川大学 | 明流泄洪洞平面弯道的压顶方法 |
CN102542108A (zh) * | 2011-12-28 | 2012-07-04 | 天津大学 | 城市大型缓流景观水体流态优化的设计方法 |
CN102663827A (zh) * | 2012-03-02 | 2012-09-12 | 天津大学 | 复杂淹没区域风暴潮洪水演进三维动态全过程仿真方法 |
WO2013078627A1 (zh) * | 2011-11-30 | 2013-06-06 | 天津空中代码工程应用软件开发有限公司 | 可压缩旋流场的数值模拟方法 |
CN104142386A (zh) * | 2013-05-10 | 2014-11-12 | 华北水利水电大学 | 固结粘性泥沙冲刷量测仪 |
CN109141817A (zh) * | 2018-09-24 | 2019-01-04 | 天津大学 | 高海拔坝工高速水气两相流研究试验方法 |
CN110046398A (zh) * | 2019-03-22 | 2019-07-23 | 西北农林科技大学 | 一种当卡水电站泄洪闸下游冲刷过程三维数值模拟方法 |
CN110984088A (zh) * | 2019-11-26 | 2020-04-10 | 昆明理工大学 | 一种竖直消力井消能率的实验计算方法 |
CN111008418A (zh) * | 2019-11-18 | 2020-04-14 | 天津大学 | 一种基于数值模拟与物理模型试验相结合的建筑群消能防冲优化设计方法 |
CN111121854A (zh) * | 2019-12-31 | 2020-05-08 | 南昌工程学院 | 测量挑流消能工消能率的装置及使用方法 |
CN111311734A (zh) * | 2020-02-14 | 2020-06-19 | 中核第四研究设计工程有限公司 | 铀尾矿库排洪系统水力特性三维数值模拟方法 |
CN111651820A (zh) * | 2020-05-27 | 2020-09-11 | 长江勘测规划设计研究有限责任公司 | 单位水体消能率和基岩抗冲流速相耦合的消能评估方法 |
CN114417473A (zh) * | 2022-01-21 | 2022-04-29 | 安徽省引江济淮集团有限公司 | 一种基于粒子追踪的优化航道汇流区的工程布置方法 |
-
2007
- 2007-02-12 CN CNA2007100515686A patent/CN101017516A/zh active Pending
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101871204A (zh) * | 2010-06-28 | 2010-10-27 | 四川大学 | 明流泄洪洞平面弯道的压顶方法 |
CN101871204B (zh) * | 2010-06-28 | 2012-08-15 | 四川大学 | 明流泄洪洞平面弯道的压顶方法 |
WO2013078627A1 (zh) * | 2011-11-30 | 2013-06-06 | 天津空中代码工程应用软件开发有限公司 | 可压缩旋流场的数值模拟方法 |
CN102542108A (zh) * | 2011-12-28 | 2012-07-04 | 天津大学 | 城市大型缓流景观水体流态优化的设计方法 |
CN102542108B (zh) * | 2011-12-28 | 2013-07-31 | 天津大学 | 城市大型缓流景观水体流态优化的设计方法 |
CN102663827A (zh) * | 2012-03-02 | 2012-09-12 | 天津大学 | 复杂淹没区域风暴潮洪水演进三维动态全过程仿真方法 |
CN102663827B (zh) * | 2012-03-02 | 2014-07-09 | 天津大学 | 复杂淹没区域风暴潮洪水演进三维动态全过程仿真方法 |
CN104142386A (zh) * | 2013-05-10 | 2014-11-12 | 华北水利水电大学 | 固结粘性泥沙冲刷量测仪 |
CN109141817A (zh) * | 2018-09-24 | 2019-01-04 | 天津大学 | 高海拔坝工高速水气两相流研究试验方法 |
CN110046398A (zh) * | 2019-03-22 | 2019-07-23 | 西北农林科技大学 | 一种当卡水电站泄洪闸下游冲刷过程三维数值模拟方法 |
CN110046398B (zh) * | 2019-03-22 | 2022-10-14 | 西北农林科技大学 | 一种当卡水电站泄洪闸下游冲刷过程三维数值模拟方法 |
CN111008418A (zh) * | 2019-11-18 | 2020-04-14 | 天津大学 | 一种基于数值模拟与物理模型试验相结合的建筑群消能防冲优化设计方法 |
CN111008418B (zh) * | 2019-11-18 | 2023-06-06 | 天津大学 | 一种基于数值模拟与物理模型试验相结合的建筑群消能防冲优化设计方法 |
CN110984088A (zh) * | 2019-11-26 | 2020-04-10 | 昆明理工大学 | 一种竖直消力井消能率的实验计算方法 |
CN111121854A (zh) * | 2019-12-31 | 2020-05-08 | 南昌工程学院 | 测量挑流消能工消能率的装置及使用方法 |
CN111311734A (zh) * | 2020-02-14 | 2020-06-19 | 中核第四研究设计工程有限公司 | 铀尾矿库排洪系统水力特性三维数值模拟方法 |
CN111651820A (zh) * | 2020-05-27 | 2020-09-11 | 长江勘测规划设计研究有限责任公司 | 单位水体消能率和基岩抗冲流速相耦合的消能评估方法 |
CN111651820B (zh) * | 2020-05-27 | 2022-04-22 | 长江勘测规划设计研究有限责任公司 | 单位水体消能率和基岩抗冲流速相耦合的消能评估方法 |
CN114417473A (zh) * | 2022-01-21 | 2022-04-29 | 安徽省引江济淮集团有限公司 | 一种基于粒子追踪的优化航道汇流区的工程布置方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101017516A (zh) | 水工水力学的水流精细模拟方法 | |
Wanchat et al. | Preliminary design of a vortex pool for electrical generation | |
Huai et al. | Numerical simulation of sandy bed erosion by 2D vertical jet | |
Wang et al. | Influence of submerged flexible vegetation on turbulence in an open-channel flow | |
CN101694678B (zh) | 一种测算人工鱼礁流场调控规模的方法 | |
Xu et al. | Flow pattern and anti-silt measures of straight-edge forebay in large pump stations. | |
Lian et al. | Numerical study on local scour characteristics of multi-bucket jacket foundation considering exposed height | |
Rezvaya et al. | The numerical research of the flow in the inlet of the high-head hydraulic turbine | |
CN101017517A (zh) | 水利水电工程水流精细模拟方法 | |
Abeykoon | Modelling and optimisation of a Kaplan turbine—A comprehensive theoretical and CFD study | |
Kim et al. | Modeling the buoyant flow of heated water discharged from surface and submerged side outfalls in shallow and deep water with a cross flow | |
Peng et al. | Numerical study of cavitation on the surface of the guide vane in three gorges hydropower unit | |
CN115034113A (zh) | 一种鱼道休息池结构确定方法及系统 | |
Zhang et al. | Numerical investigation of scour around the monopile using CFD-DEM coupling method | |
Bai et al. | Pressure distributions of stepped spillways with different horizontal face angles | |
Li et al. | Formation and influencing factors of free surface vortex in a barrel with a central orifice at bottom | |
Wang et al. | Semi-empirical wake structure model of rotors using joint axial momentum theory and DES-SA method | |
Cao et al. | Numerical study on flow structure near two impermeable trapezoid submerged breakwaters on slop bottoms | |
Yu et al. | Wear Behavior and Vortex Characteristics of Y-Type Screen Filters with Various Inclination Angles and Inlet Flow Velocities: Numerical and Experimental Study | |
Fang et al. | Wave and current effects on a buoyant jet: An integral model | |
Quiñones et al. | Numerical Study of the Hydrodynamics of Mangrove-Inspired Structures for Coastal Protection | |
Deng et al. | Proposal of actuator line-immersed boundary coupling model for tidal stream turbine modeling with hydrodynamics upon scouring morphology | |
Klostermann et al. | Investigation of the Gas‐Liquid Flow in a Stopper Rod Controlled SEN | |
Taşar et al. | Numerical simulation of channel flow using submerged vane in river arrangements | |
Bin et al. | Research on the influence law of submerged depth and flow rate of diversion canal on free surface vortex in coastal nuclear power plant |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20070815 |