CN117454724A - 一种通过cfd确定护城河式防洪堤坝几何参数的方法 - Google Patents
一种通过cfd确定护城河式防洪堤坝几何参数的方法 Download PDFInfo
- Publication number
- CN117454724A CN117454724A CN202311786133.1A CN202311786133A CN117454724A CN 117454724 A CN117454724 A CN 117454724A CN 202311786133 A CN202311786133 A CN 202311786133A CN 117454724 A CN117454724 A CN 117454724A
- Authority
- CN
- China
- Prior art keywords
- river
- dam
- urban
- flood control
- downstream
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 73
- 238000004364 calculation method Methods 0.000 claims abstract description 93
- 238000004088 simulation Methods 0.000 claims abstract description 81
- 238000012360 testing method Methods 0.000 claims abstract description 36
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 claims description 122
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 111
- 238000011144 upstream manufacturing Methods 0.000 claims description 86
- 230000008859 change Effects 0.000 claims description 13
- 230000004888 barrier function Effects 0.000 claims description 10
- 230000007423 decrease Effects 0.000 claims description 8
- 230000005484 gravity Effects 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 241000084490 Esenbeckia delta Species 0.000 claims description 2
- 230000000149 penetrating effect Effects 0.000 claims description 2
- 239000012530 fluid Substances 0.000 abstract description 23
- 230000021715 photosynthesis, light harvesting Effects 0.000 abstract description 3
- 108091006146 Channels Proteins 0.000 description 68
- 238000004422 calculation algorithm Methods 0.000 description 6
- 230000008878 coupling Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000006978 adaptation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000007123 defense Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000005514 two-phase flow Effects 0.000 description 2
- 102000010637 Aquaporins Human genes 0.000 description 1
- 108010063290 Aquaporins Proteins 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000002620 method output Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012932 thermodynamic analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Revetment (AREA)
Abstract
本发明公开了一种通过CFD确定护城河式防洪堤坝几何参数的方法,包括先对数值模拟的方法进行验证,对护城河式防洪堤坝渠道同时进行物理模型试验和数值计算,通过对比关键参数得到满足计算精度要求的计算软件设置方法;然后通过得出的计算流体动力学设置方法对原型护城河式防洪堤坝渠道进行数值模拟,通过关键指标对比确定护城河式防洪堤坝渠道关键几何参数。在数值计算时,引入云计算以提高模拟精度和效率。本发明能够准确合理地确定护城河式防洪堤坝的几何参数,显著提高防洪堤坝的消能效率的同时有效降低护城河防洪堤坝的工程造价。
Description
技术领域
本发明属于流体动力学领域,具体涉及一种通过CFD确定护城河式防洪堤几何尺寸的方法。
背景技术
近年来,由于全球气候变化和恶劣的天气条件,洪水事件变得更加频繁地发生。研究者提出了一种多结构组成的复合防御体系,包括植被与堤岸组合、堤坝与护城河组合、双堤坝、堤坝与墩桩组合等多线防御策略。在这些措施中,植被与堤岸的配合有利于通过水跃缓解洪水能量,但是植被层也会受洪水的不断侵蚀而遭受破坏;而双堤坝方案中,前一个堤坝确实为后一个堤坝提供了阻力并降低了洪水能量,但是随着洪水的冲刷,第二个堤坝(一般为土石坝)存在冲毁的风险;堤坝与墩桩组合中,虽然墩桩能降低洪水能量,但是建设固定式墩桩由于占用耕地或阻挡交通而显得很不经济,尽管可将墩桩改为抛石、块石等混凝土块,又可能被洪水冲到下游而对下游水利工程结构或房屋建筑造成干扰甚至破坏。然而,护城河与堤坝的组合中,护城河中的水体似于下游水垫层,可以有效降低流体作用力和溢流体积,从而显著减小底部流速,实现耗散洪水能量目的以保护消力池底部。但是,护城河结构(水垫)作为一种对抗由堤坝溢流和溃坝引起的洪灾的对策并没有得到全面的研究,尤其是护城河式防洪堤坝几何尺寸的确定。由于护城河式防洪堤坝的几何尺寸不仅影响着防冲消能的效率,更直接决定了此类水利工程的造价,因此需要通过先进的技术手段合理确定这些工程关键参数。
随着计算流体动力学(CFD)和云计算技术的快速发展,由于其具有建模便捷、可直接对原型进行模拟、试验重复性较好等优点,数值模拟方法逐渐成为一种非常有效的研究手段。但是现有技术中还没有利用数值试验手段确定护城河式防洪堤坝几何参数的方法。中国专利公开号CN106844913B中公开了一种基于三维计算流体动力学(CFD)的滞留气团热力学特性模拟方法,该专利利用三维计算流体动力学(CFD)进行了滞留气团的热力学分析,但是其方法不能用于确定护城河式防洪堤坝的几何参数。因此,如何利用计算流体动力学(CFD)方法来合理确定护城河式堤坝的几何尺寸是一个非常值得研究的课题,它将有效指导实际工程中护城河式防洪堤坝的设计与优化。
发明内容
本发明的目的为了克服现有技术中存在的不足,本发明提供一种通过CFD确定护城河式防洪堤坝几何参数的方法,能够简便、合理地确定护城河式防洪堤坝的几何布置,在显著提高防洪堤坝的消能效率的同时,有效降低了护城河防洪堤坝的工程造价。
本发明的技术方案:一种通过CFD确定护城河式防洪堤坝几何参数的方法,包括以下步骤:
步骤S1:对缩尺的护城河式防洪堤坝流场进行物理模型试验,建立缩尺的护城河式防洪堤坝渠道物理模型,物理模型试验时,测得护城河式防洪堤坝渠道对称纵断面上、堤坝上游不同过流横断面水深hi;测得护城河式防洪堤坝渠道中护城河内水跃形态及河道中心对称面矢量分布;
步骤S2:对缩尺的护城河式防洪堤坝渠道物理模型进行数值计算,数值模拟时测得护城河式防洪堤坝渠道对称纵断面上、堤坝上游不同过流横断面对应水深Hi;
步骤S3:对数值模拟时对应水深Hi数据集与物理模型试验时水深hi数据集进行逐点对比,找到数据集中所有相对误差绝对值均不超过η%的计算工况,即|Hi-hi|/hi≤η%,输出满足水深模拟精度要求的数值模拟方法;
步骤S4:对步骤S1中缩尺的护城河式防洪堤坝渠道物理模型进行数值计算,采用步骤S3中输出满足水深模拟精度要求的数值模拟方法,输出护城河式防洪堤坝渠道中对应的护城河内水跃形态以及河道中心对称面矢量分布;
步骤S5:对步骤S4中数值模拟时护城河内水跃形态以及河道中心对称面矢量分布中找到与物理模型试验时护城河内水跃形态及河道中心对称面矢量分布相同的数值模拟方案,输出对应的满足护城河内水跃形态及河道中心对称面矢量分布模拟精度要求的数值模拟方法;
步骤S6:确定护城河式防洪堤坝下游护城河深度Dm的值;
步骤S7:确定护城河式防洪堤坝下游护城河长度Lm的值;
步骤S8:确定护城河式防洪堤坝的下游坡度θd的值;
步骤S9:确定护城河式防洪堤坝的上游坡度θu的值;
步骤S10:确定护城河式防洪堤坝的坝顶宽度Lw;
步骤S11:确定护城河式防洪堤坝的堤坝高度hL。
进一步的,步骤S1中建立缩尺的护城河式防洪堤坝渠道物理模型,测得关键参数还包括有;坝顶对称断面处水深Yc,下游尾水断面处水深Yt,坝顶对称断面与下游尾水断面间距M;
其中,护城河式防洪堤坝几何参数为:堤坝高度hL,坝顶宽度Lw,上游坡度θu,下游坡度θd,在堤坝正下游护城河长度为Lm和深度为Dm,上游坝顶拐点到上游坝脚拐点的水平距离为Hu,下游坝顶拐点到下游坝脚拐点的水平距离为Hd,其中上游坡度θu为上游坡度对应的堤坝高度hL与水平距离Hu的比值,下游坡度θd为下游坡度对应的堤坝高度hL与水平距离Hd的比值;
其中,具体参数为:
渠道几何长度为L,渠道宽度为B,护城河式防洪堤坝横跨渠道贯穿布置;
物理模型试验采用正态物理模型,综合考虑物理模型水流在阻力平方区要求,选取物理模型线性比尺λl,物理模型采用重力相似准则设计,角度比λγ,流速比λv=λl 0.5,流量比λQ=λl 2.5,糙率比λn=λl 1/6,护城河式防洪堤坝渠道流量Qp=λQ×Qt。
进一步的,步骤S1中建立缩尺的护城河式防洪堤坝渠道物理模型,物理模型试验具体步骤为:
步骤S101:给定任意模型试验流量Qt,待缩尺的护城河式防洪堤坝渠道中水流稳定后,依次测得护城河式防洪堤坝渠道对称纵断面上、堤坝上游不同过流横断面水深hi,i为不同过流断面的编号;
步骤S102:改变模型试验流量Qt,同时分别改变上游坡度θu,下游坡度θd,下游护城河长度Lm和深度Dm,测得不同工况下坝顶对称断面处水深Yc,下游尾水断面处水深Yt,并通过高速相机摄像得到护城河内水跃形态,以及通过粒子图像测速仪PIV测得河道中心对称面矢量分布;
步骤S103:分别改变上游坡度θu,下游坡度θd,下游护城河长度Lm和深度Dm,具体方法为:
保持下游坡度θd,下游护城河长度Lm和深度Dm不变,仅改变上游坡度θu,或保持上游坡度θu,下游护城河长度Lm和深度Dm不变,仅改变下游坡度θd,或保持上游坡度θu,下游坡度θd,下游护城河深度Dm不变,仅改变下游护城河长度Lm,或保持改变上游坡度θu,下游坡度θd,下游护城河长度Lm不变,仅改变下游护城河深度Dm。
进一步的,步骤S6中确定原型护城河式防洪堤坝下游护城河深度Dm的值;具体为:
对护城河式防洪堤坝高度hL,坝顶宽度Lw,上游坡度θu,下游坡度θd,护城河长度Lm在合理范围内进行任意取值,在相同护城河式防洪堤坝渠道流量Qp条件下,分别对a个不同深度的护城河式防洪堤坝渠道进行数值计算,当坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE达到最大时,输出对应的护城河深度Dm为确定值,记作D* m;
护城河深度D* m的具体确定步骤为:
步骤S61:a个不同深度的护城河式防洪堤坝渠道得到a个不同护城河深度Dm,具体值记为χ,即护城河深度值χ1, χ2, χ3,…, χa,其中护城河深度值χ1, χ2, χ3,…, χa逐渐增大,即χ1< χ2< χ3<…< χa,对a个不同深度的护城河式防洪堤坝渠道进行数值计算,分别输出坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE1,ΔE2,ΔE3,…, ΔEa;
步骤S62:当能量差的比值ΔE1=ΔE2=ΔE3=…= ΔEa时,则以Δχ的间距减小χ1分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生减小时,输出前一次数值模拟对应的护城河深度,即第k1次计算时发生能量差的比值ΔE减小,则确定的护城河深度D* m为χ1-(k1-1)×Δχ;
步骤S63:当a个不同护城河深度的具体值χ对应的能量差的比值存在波动,即ΔE1=ΔE2=ΔE3=…=ΔEo<ΔEo+1=…= ΔEa,ΔEo对应的护城河深度为χo,则以Δχ的间距增大χo分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE与ΔEo+1相等时,输出当前数值模拟对应的护城河深度,即第k2次计算得到的能量差的比值ΔE与ΔEo+1相等时,则确定的护城河深度D* m为χo+k2×Δχ;
步骤S64:当a个不同护城河深度的具体值χ对应的能量差的比值单调递增,即ΔE1<ΔE2<ΔE3<…< ΔEa,则以Δχ的间距增大χa分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生相等时,输出前一次数值模拟对应的护城河深度,即第k3次计算时发生能量差的比值ΔE相等,则确定的护城河深度D* m为χa+(k3-1)×Δχ。
进一步的,步骤S7中确定原型护城河式防洪堤坝下游护城河长度Lm的值;具体为:
给定步骤S6中确定的护城河深度D* m,对护城河式防洪堤坝高度hL,坝顶宽度Lw,上游坡度θu,下游坡度θd在合理范围内进行任意取值,在相同护城河式防洪堤坝渠道流量Qp条件下,分别对b个不同护城河长度的护城河式防洪堤坝渠道进行数值计算,当坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE达到最大时,输出对应的护城河长度Lm为确定值,记作L* m;
护城河长度L* m的具体确定步骤为:
步骤S71:b个不同长度的护城河式防洪堤坝渠道得到b个不同护城河长度Lm,具体值记为ζ,即护城河长度值ζ1, ζ2, ζ3,…, ζb,其中ζ1, ζ2, ζ3,…, ζb逐渐增大,即ζ1< ζ2< ζ3<…< ζb,对b个不同长度的护城河式防洪堤坝渠道进行数值计算,分别输出坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE1, ΔE2, ΔE3,…,ΔEb;
步骤S72:当能量差的比值ΔE1=ΔE2=ΔE3=…= ΔEb时,则以Δζ的间距减小ζ1分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生减小时,输出前一次数值模拟对应的护城河长度,即第m1次计算时发生能量差的比值ΔE减小,则确定的护城河长度L* m为ζ1-(m1-1)×Δζ;
步骤S73:当b个不同护城河长度的ζ对应的能量差的比值ΔE存在波动,即ΔE1=ΔE2=ΔE3=…=ΔEo<ΔEo+1=…= ΔEb,ΔEo对应的护城河长度为ζo,则以Δζ的间距增大ζo分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE与ΔEo+1相等时,输出当前数值模拟对应的护城河长度,即第m2次计算得到的能量差的比值ΔE与ΔEo+1相等时,则确定的护城河长度L* m为ζo+m2×Δζ;
步骤S74:当b个不同护城河深度的ζ对应的能量差的比值单调递增,即ΔE1<ΔE2<ΔE3<…< ΔEb,则以Δζ的间距增大ζb分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生相等时,输出前一次数值模拟对应的护城河长度,即第m3次计算时发生能量差的比值ΔE相等,则确定的护城河长度L* m为ζb+(m3-1)×Δζ。
进一步的,步骤S8中确定原型护城河式防洪堤坝的下游坡度θd的值;具体为:
给定步骤S6中确定的护城河深度D* m和步骤S7中确定的护城河长度L* m,对护城河式防洪堤坝坝高度hL,坝顶宽度Lw,上游坡度θu在合理范围内进行任意取值,在不同护城河式防洪堤坝渠道流量Qp条件下分别对c个不同下游坡度θd的护城河式防洪堤坝渠道进行数值计算,当坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE达到最大时,输出对应的护城河下游坡度θd为确定值,记作θ* d;
护城河下游坡度θ* d的具体确定步骤为:
步骤S81:c个不同下游坡度的护城河式防洪堤坝渠道得到c个不同下游坡度θd,具体值记为п,即护城河下游坡度值п1, п2, п3,…, пc,其中п1, п2, п3,…, пc逐渐增大,即п1< п2< п3<…< пc,其中п1是下游坡度θd允许范围的最小值,пc是下游坡度θd允许范围的最大值,对c个下游坡度的护城河式防洪堤坝渠道进行数值计算,分别输出坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE1, ΔE2, ΔE3,…, ΔEc;
步骤S82:当能量差的比值单调递增时,即ΔE1<ΔE2<ΔE3<…< ΔEc,则输出ΔEc对应的пc为最终的下游坡度值θ* d;
步骤S83:当能量差的比值单调递减时,即ΔE1>ΔE2>ΔE3>…> ΔEc,则输出ΔE1对应的п1为最终的下游坡度值θ* d。
进一步的,步骤S9中确定原型护城河式防洪堤坝的上游坡度θu的值;具体为:
给定步骤S6中确定的护城河深度D* m、步骤S7中确定的护城河长度L* m和步骤S8中确定的下游坡度值θ* d,对护城河式防洪堤坝坝高度hL,坝顶宽度Lw在合理范围内进行取值,在不同护城河式防洪堤坝渠道流量Qp条件下分别对e个不同上游坡度θu的护城河式防洪堤坝渠道进行数值计算,当坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE达到最大时,输出对应的护城河上游坡度θu为确定值,记作θ* u;
护城河上游坡度θ* u的具体确定步骤为:
步骤S91:e个不同上游坡度的护城河式防洪堤坝渠道得到e个不同上游坡度θu,具体值记为δ,即护城河上游坡度值δ1, δ2, δ3,…, δe,其中δ1, δ2, δ3,…, δe逐渐增大,即δ1<δ2<δ3<…<δe,其中δ1是上游坡度θu允许范围的最小值,δe是上游坡度θu允许范围的最大值,对e个上游坡度的护城河式防洪堤坝渠道进行数值计算,分别输出坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE1, ΔE2, ΔE3,…, ΔEe;
步骤S92:当能量差的比值单调递增时,即ΔE1<ΔE2<ΔE3<…<ΔEe,则输出ΔEe对应的δe为最终的下游坡度值θ* u;
步骤S93:当能量差的比值单调递减时,即ΔE1>ΔE2>ΔE3>…>ΔEe,则输出ΔE1对应的δ1为最终的下游坡度值θ* u。
进一步的,步骤S10中确定原型护城河式防洪堤坝的坝顶宽度Lw;具体为:
给定步骤S6中确定的护城河深度D* m、步骤S7中确定的护城河长度L* m、步骤S8中确定的下游坡度值θ* d和步骤S9中确定的上游坡度值θ* u,对护城河式防洪堤坝坝高度hL在合理范围内进行取值,在不同护城河式防洪堤坝渠道流量Qp条件下分别对f个不同坝顶宽度Lw的护城河式防洪堤坝渠道进行数值计算,当坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE达到最大时,输出对应的坝顶宽度Lw为确定值,记作L* w;
防洪堤坝的坝顶宽度L* w的具体确定步骤为:
步骤S111:f个不同坝顶宽度的护城河式防洪堤坝渠道得到f个不同坝顶宽度Lw,具体值记为p,即护城河上游坡度值p1, p2, p3,…, pf,其中p1, p2, p3,…, pf逐渐增大,即p1<p2<p3<…<pf,其中p1是坝顶宽度Lw允许范围的最小值,pf是坝顶宽度Lw允许范围的最大值,对f个坝顶宽度的护城河式防洪堤坝渠道进行数值计算,分别输出坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE1, ΔE2, ΔE3,…, ΔEf;
步骤S112:当能量差的比值ΔE1=ΔE2=ΔE3=…=ΔEf时,则以Δp的间距减小p1分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生减小时,输出前一次数值模拟对应的护城河长度,即第j1次计算时发生能量差的比值ΔE减小,则确定的防洪堤坝的坝顶宽度L* w为p1-(j1-1)×Δp;
步骤S113:当o个不同坝顶宽度的p对应的能量差的比值单调递增,而剩余的j-o个不同坝顶宽度的p对应的能量差的比值相等,假设其中一种情况为ΔE1<ΔE2<ΔE3<…<ΔEo<ΔEo+1=…= ΔEf,则以Δp的间距减小po+1分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生减小时,输出前一次数值模拟对应的护城河长度,即第j2次计算时发生能量差的比值ΔE减小,则确定的防洪堤坝的坝顶宽度L* w为po+1-(j2-1)×Δp;
步骤S114:当f个不同长度的p对应的能量差的比值单调递增,即ΔE1<ΔE2<ΔE3<…<ΔEf,则以Δp的间距增大pf分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生相等时,输出前一次数值模拟对应的护城河长度,即第j3次计算时发生能量差的比值ΔE相等,则确定的防洪堤坝的坝顶宽度L* w为pf+(j3-1)×Δp。
进一步的,步骤S11:确定原型护城河式防洪堤坝的堤坝高度hL;具体为:
给定步骤S6中确定的护城河深度D* m、步骤S7中确定的护城河长度L* m、步骤S8中确定的下游坡度值θ* d、步骤S9中确定的上游坡度值θ* u和步骤S10中确定的坝顶宽度L* w,在不同护城河式防洪堤坝渠道流量Qp条件下分别对s个不同堤坝高度hL的护城河式防洪堤坝渠道进行数值计算,当坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE达到最大时,输出对应的堤坝高度hL为确定值,记作h* L;
防洪堤坝的堤坝高度h* L的具体确定步骤为:
步骤S115:s个不同堤坝高度的护城河式防洪堤坝渠道得到s个不同堤坝高度hL,具体值记为,即护城河上游坡度值/>, />, />,…, />,其中/>, />, />,…, />逐渐增大,即/></></><…</>,其中/>是堤坝高度hL允许范围的最小值,/>是堤坝高度hL允许范围的最大值,对s个堤坝高度的护城河式防洪堤坝渠道进行数值计算,分别输出坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE1, ΔE2, ΔE3,…, ΔEs;
步骤S116:当能量差的比值单调递增时,即ΔE1<ΔE2<ΔE3<…<ΔEs,则输出ΔEs对应的为确定防洪堤坝的堤坝高度,记作h* L;
步骤S117:当能量差的比值单调递减时,即ΔE1>ΔE2>ΔE3>…>ΔEc,则输出ΔE1对应的为确定防洪堤坝的堤坝高度,记作h* L。
进一步的,步骤S6-步骤S11中,坝顶对称断面A-A与下游尾水断面Z-Z之间能量差与坝顶对称断面总能量的比值ΔE,计算方法为:
先计算坝顶对称断面A-A流速VA为渠道流量Q除以坝顶对称断面A-A过流面积,坝顶对称断面A-A过流面积为坝顶对称断面处水深Yc与渠道宽度B的乘积,则VA=Q/(Yc×B),那么坝顶对称断面A-A的总能量为EA=VA²/2g+Yc+hL;
再计算下游尾水断面Z-Z流速VZ为渠道流量Q除以下游尾水断面Z-Z过流面积,下游尾水断面Z-Z过流面积为下游尾水断面处水深Yt与渠道宽度B的乘积,则VZ= Q/(Yt×B),那么下游尾水断面Z-Z的总能量为EZ=VZ²/2g+Yt;
最后计算坝顶对称断面A-A与下游尾水断面Z-Z之间能量差与坝顶对称断面总能量的比值,简称能量差的比值ΔE,ΔE=(EA-EZ)/ EA。
进一步的,步骤S2中对缩尺的护城河式防洪堤坝渠道物理模型进行数值计算,计算收敛后导出结果,通过CFD后处理软件进行后处理;对缩尺的护城河式防洪堤坝渠道进行数值模拟计算,具体为:
建立与步骤S1中缩尺的护城河式防洪堤坝渠道相同的三维物理模型,并将三维物理模型进行网格剖分,输出后缀名为.mesh的计算文件,导入计算流体动力学CFD软件进行数值计算;在将三维物理模型网格剖分时采用不同尺度网格,从而得到若干不同的网格方案,针对各个网格方案在计算流体动力学软件计算中选取不同的湍流模型及不同的数值算法。
进一步的,步骤S6,步骤S7,步骤S8,步骤S9,步骤S10,步骤S11中采用的数值模拟方法与步骤S5输出的满足水跃及矢量分布模拟精度要求的数值模拟方法相同。
进一步的,步骤S6,步骤S7,步骤S8,步骤S9,步骤S10,步骤S11,在护城河式防洪堤渠道原型进行计算流体动力学CFD计算中引入云计算技术。
本发明的有益效果:本发明通过CFD确定护城河式防洪堤坝几何参数的方法,能够简便、准确地确定护城河式防洪堤坝的几何参数,提高了水利工程中布置护城河式防洪堤坝的规划合理性,可以大量节省建造护城河式防洪堤坝的工程造价,又有效地避免了护城河式防洪堤坝下游区域的洪水破坏。
附图说明
图1为护城河式防洪堤坝渠道结构立面示意图;
图2为护城河式防洪堤坝渠道结构平面示意图;
图3为物理试验和数值模拟中护城河式防洪堤坝上游渠道对称纵断面上水深对比图;
图4为数值模拟中护城河式防洪堤坝渠道中护城河内水跃形态及河道中心对称面矢量分布;(a)河道中心对称面矢量分布;(b)护城河内水跃形态;
图5能量差的比值ΔE随不同防洪堤坝下游护城河深度Dm的变化关系;
图6能量差的比值ΔE随不同防洪堤坝下游护城河长度Lm的变化关系;
图7能量差的比值ΔE随不同防洪堤坝的下游坡度θd的变化关系;
图8能量差的比值ΔE随不同防洪堤坝的上游坡度θu的变化关系;
图9能量差的比值ΔE随不同防洪堤坝的坝顶宽度Lw的变化关系;
图10能量差的比值ΔE随不同防洪堤坝的堤坝高度hL的变化关系。
具体实施方式
下面结合附图对本发明作更进一步的说明。
本发明的通过CFD确定护城河式防洪堤坝几何参数的方法,包括先对数值模拟的方法进行验证,通过对护城河式防洪堤坝渠道缩尺后的物理模型进行数值计算和物理试验,并对护城河式防洪堤坝渠道上游不同过流断面的水深、护城河内水跃形态和渠道纵断面速度矢量分布进行采集与对比,得到合理的计算流体动力学数值方案包括网格剖分、湍流模型和数值方法;再通过上述数值方案对原型条件下的护城河式防洪堤坝渠道进行计算流体动力学计算,通过比较坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值的大小,将能量差的比值最大时的几何参数确定为最终尺寸。本发明通过缩尺物理模型试验和相应的物理模型计算流体动力学数值计算进行对比,用物理模型试验来约束数值计算方法,再通过得出的数值计算方法对原型条件下的护城河式防洪堤坝进行计算以确定原型条件下护城河式防洪堤坝几何参数,如此使得原型条件下的数值计算的结果更科学有效;使原型条件下构建的计算流体动力学数值计算模型(数值计算模型是指数值模拟构建的模型)能够满足实际流动条件,提高了计算流体动力学数值模拟的准确性,能够准确地确定护城河式防洪堤坝几何参数。
实施例1
本实施例的通过CFD确定原型护城河式防洪堤坝几何参数的方法,具体步骤如下:
步骤S1:对缩尺的护城河式防洪堤坝渠道进行物理模型试验,建立护城河式防洪堤坝渠道物理模型;
其护城河式防洪堤坝渠道几何长度L为10m,渠道宽度B为0.5m,上游坡度θu为1:2,下游坡度θd为1:2或1:3,坝顶宽度Lw为40mm,坝顶高度hL为100mm,下游护城河长度Lm为350mm,下游护城河深度Dm为50mm,坝顶对称断面为A-A,下游尾水断面为Z-Z,坝顶对称断面A-A与下游尾水断面Z-Z之间的距离M为7m,分别进行两个物理模型试验;其一,在模型试验流量Qt为42L/s的条件下,测试沿堤坝渠道纵对称断面上、堤坝上游不同过流横断面对应水深hi,此时护城河式防洪堤坝渠道的几何参数:上游坡度θu为1:2、下游坡度θd为1:2、坝顶宽度Lw为40mm、顶高度hL为100mm、下游无护城河的条件下进行;其二,在模型试验流量Qt为42L/s的条件下,通过高速相机采集护城河内水跃形态和通过粒子图像测速仪(PIV)采集渠道中心对称面矢量分布,是在上游坡度θu为1:2、下游坡度θd为1:3、坝顶宽度Lw为40mm、顶高度hL为100mm、游护城河长度为350mm、护城河深度为50mm的条件下进行(护城河式防洪堤坝渠道结构如图1和图2所示);
物理试验模型采用正态模型,综合考虑物理试验模型水流在阻力平方区要求,拟选取模型线性比尺为10,物理模型采用重力相似准则设计,角度比为1,流速比为λv=100.5=3.16,流量比λQ=102.5=316.23,糙率比λn=101/6=1.47。
步骤S2:对步骤S1中缩尺的护城河式防洪堤坝渠道模型进行数值计算,建立缩尺的护城河式防洪堤坝渠道三维数值模型,并将护城河式防洪堤坝渠道三维数值模型进行网格剖分,输出后缀名为.mesh的计算文件,导入计算流体动力学软件Fluent进行数值计算,计算收敛后导出结果;
通过计算流体动力学后处理软件CFD-Post进行后处理,在上游坡度θu为1:2、下游坡度θd为1:2、坝顶宽度Lw为40mm、下游无护城河的条件下在相同模型试验流量Qt条件进行计算获得沿堤坝渠道纵对称断面上、堤坝上游不同过流横断面对应水深Hi;
在将三维物理模型网格剖分时采用不同尺度网格,从而得到若干不同的网格方案,针对各个网格方案在计算流体动力学软件Fluent计算中选取不同的湍流模型及不同的数值算法。
步骤S3:对步骤S2中得到的数值模拟时模型护城河式防洪堤坝渠道纵对称断面上、堤坝上游不同过流横断面对应水深Hi数据集与步骤S1中物理模型试验时模型护城河式防洪堤坝渠道纵对称断面上、堤坝上游不同过流横断面对应水深hi数据集进行逐点对比,找到所有数值模拟中堤坝上游水深Hi与物理模型试验中堤坝上游水深hi相对误差绝对值均不超过2.5%的计算工况,即|Hi-hi|/hi≤2.5%时,输出满足水深模拟精度要求的数值模拟方法。图3为物理试验和数值模拟中护城河式防洪堤坝上游渠道纵对称断面上水深对比图(此处仅展示其中1种情况),其相对误差的绝对值均小于2.5%,输出对应的数值模拟方法,包括网格尺寸、湍流模型、数值算法。
步骤S4:对步骤S1中缩尺的护城河式防洪堤坝渠道模型进行数值计算,护城河式防洪堤坝几何尺寸分别与步骤S1第二个物理模型试验几何参数中相同,其上游坡度θu为1:2、下游坡度θd为1:3、坝顶宽度Lw为40mm、坝顶高度hL为100mm、游护城河长度为350mm、护城河深度为50mm,数值模拟方法分别采用步骤S3输出的满足水深精度的数值模拟方法,输出数值模拟对应的护城河内水跃形态以及河道中心对称面矢量分布。
步骤S5:对步骤S4中得到数值模拟时模型护城河式防洪堤坝渠道中护城河内水跃形态以及河道中心对称面矢量分布与步骤S1中物理模型试验护城河式防洪堤坝渠道中护城河内水跃形态及河道中心对称面矢量分布最相近的1种数值模拟方案(如图4所示),输出对应的满足水跃及矢量分布模拟精度要求的数值模拟方法,输出对应的数值模拟方法,包括网格尺寸、湍流模型、数值算法,拟采用单位体积1×1×1的六面体结构化网格,湍流模型为可实现模型(Realizable k-ε),控制方程的离散方式为有限体积法(FVM);扩散项采用二阶中心差分格式,对流项采用快速格式(QUICK),压力和速度的耦合采用简单耦合算法(SIMPLEC),计算方式采用并行计算,数值模拟方法采用气液两相流(VOF)方法。
步骤S6:确定原型护城河式防洪堤坝下游护城河深度Dm的值;
对护城河式防洪堤坝坝高度hL为1m,坝顶宽度Lw为0.4m,上游坡度θu为1:2,下游坡度θd为1:2,护城河长度Lm为2.5m,在相同护城河式防洪堤坝渠道流量Qp为λQ×Qt=316.23×0.042=13.28m³/s(Qp计算方式,下同)条件下分别对8个不同护城河长度的护城河式防洪堤坝渠道进行数值计算,分别计算坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE,以其中一种情况说明能量差的比值ΔE的计算方法,当护城河深度值χ取0.30m时,A-A断面坝上水深为1.86m,流速为1.43m/s,总能头EA为1.96m;Z-Z断面水深为0.69m,流速为3.87m/s,总能头EZ为1.45m,那么则能量差的比值为(EA-EZ)/ EA=(1.96-1.45)÷1.96=26%,逐一计算后发现属于“ΔE1=ΔE2=ΔE3=…=ΔEo<ΔEo+1=…= ΔEa”的情况,此时χo=0.30m,取Δχ=0.05m,则发现当χ增大到0.35m时,与护城河深度为0.40m的能量差的比值趋于相等,共计增加1次,则χo+k2×Δχ=0.30+1×0.05=0.35m,那么护城河式防洪堤坝下游护城河深度Dm取值为0.35m(如图5所示),图5中横坐标为护城河深度Dm,单位m,纵坐标为能量差的比值ΔE,单位(%);
步骤S7:确定原型护城河式防洪堤坝下游护城河长度Lm的值;
给定步骤S6中确定的护城河长度D* m=0.35m,对护城河式防洪堤坝坝高度hL为1m,坝顶宽度Lw为0.40m,上游坡度θu为1:2,下游坡度θd为1:2,在相同护城河式防洪堤坝渠道流量Qp为13.28m³/s条件下分别对5个不同护城河长度的护城河式防洪堤坝渠道进行数值计算,发现属于“ΔE1=ΔE2=ΔE3=…=ΔEo<ΔEo+1=…= ΔEb”的情况,此时ζo=1.20m,取Δζ=0.05m,则发现当ζ增大到1.50m时,与护城河长度为1.55mm的能量差的比值趋于相等,共计增加6次,则ζo+m2×Δζ=1.20+6×0.05=1.50m,那么护城河式防洪堤坝下游护城河长度Lm取值为1.50m(如图6所示),图6中横坐标为护城河长度Lm,单位m,纵坐标为能量差的比值ΔE,单位(%);
步骤S8:确定原型护城河式防洪堤坝的下游坡度θd的值;
给定步骤S6中确定的护城河深度D* m=0.35mm、步骤S7中确定的护城河长度L* m=1.50m,对护城河式防洪堤坝坝高度hL为1.0m,坝顶宽度Lw为0.40mm,上游坡度θu为1:2,在7个不同护城河式防洪堤坝渠道流量Qp为1.90~13.28m³/s条件下分别对3个不同下游坡度θd,即下游坡度分别为1:1, 1:2, 1:3的护城河式防洪堤坝渠道进行数值计算,发现属于“ΔE1<ΔE2<ΔE3<…< ΔEc”的情况,即随着下游坡度的变缓,能量差的比值在增大,那么护城河式防洪堤坝的下游坡度θd的值取值为1:3(如图7所示),图7中横坐标为渠道流量Qp的无量纲表示形式,表达式为Qp/B/g^0.5/Hi^1.5,无单位,纵坐标为能量差的比值ΔE,单位(%);
步骤S9:确定原型护城河式防洪堤坝的上游坡度θu的值;
给定步骤S6中确定的护城河深度D* m=0.35m、步骤S7中确定的护城河长度L* m=1.50m和步骤S8中确定的下游坡度值θ* d=1:3,对护城河式防洪堤坝坝高度hL为1.0m,坝顶宽度Lw为0.40m,在7个不同护城河式防洪堤坝渠道流量Qp为1.90~13.28m³/s条件下分别对3个不同上游坡度θu,即上游坡度分别为1:1, 1:2, 1:3的护城河式防洪堤坝渠道进行数值计算,发现属于“ΔE1>ΔE2>ΔE3>…> ΔEe”的情况,即随着下游坡度的变陡,能量差的比值在增大,那么护城河式防洪堤坝的下游坡度θu的值取值为1:1(如图8所示),图8中横坐标为渠道流量Qp的无量纲表示形式,表达式为Qp/B/g^0.5/Hi^1.5,无单位,纵坐标为能量差的比值ΔE,单位(%);
步骤S10:确定原型护城河式防洪堤坝的坝顶宽度Lw;
给定步骤S6中确定的护城河深度D* m=0.35m、步骤S7中确定的护城河长度L* m=1.50m、步骤S8中确定的下游坡度值θ* d=1:3和步骤S9中确定的上游坡度值θ* u=1:1,对护城河式防洪堤坝坝高度hL为1.0m,在7个不同护城河式防洪堤坝渠道流量Qp为1.90~13.28m³/s条件下分别对3个不同坝顶宽度Lw,即坝顶宽度分别为0.40m, 0.60m, 0.80m的护城河式防洪堤坝渠道进行数值计算,发现属于“ΔE1<ΔE2<ΔE3<…<ΔEf”的情况,即随着坝顶宽度的变宽,能量差的比值在增大,则ΔEf对应的坝顶宽度值pf为0.80m,以0.20m为间隔逐渐增大,那么增加2次后发现坝顶宽度1.00m和1.20m的能量差的比值相等,那么护城河式防洪堤坝的坝顶宽度Lw的值取值为0.80+(2-1)×0.20=1.00m(如图9所示),图9中横坐标为渠道流量Qp的无量纲表示形式,表达式为Qp/B/g^0.5/Hi^1.5,无单位,纵坐标为能量差的比值ΔE,单位(%);
步骤S11:确定原型护城河式防洪堤坝的堤坝高度hL;
给定步骤S6中确定的护城河深度D* m=0.35m、步骤S7中确定的护城河长度L* m=1.50m、步骤S8中确定的下游坡度值θ* d=1:3、步骤S9中确定的上游坡度值θ* u=1:1和步骤S10中确定的坝顶宽度L* w=1.00m,在7个不同护城河式防洪堤坝渠道流量Qp为1.90~13.28m³/s条件下分别对3个不同堤坝高度hL,即堤坝高度分别为1.00m, 1.50m, 2.00m的护城河式防洪堤坝渠道进行数值计算,发现属于“ΔE1>ΔE2>ΔE3>…> ΔEs”的情况,即随着堤坝高度的变低,能量差的比值在增大,则ΔE1对应的坝顶高度值为1.00m,那么护城河式防洪堤坝的堤坝高度h* L的值取值为1.00m(如图10所示),图10中横坐标为渠道流量Qp的无量纲表示形式,表达式为Qp/B/g^0.5/Hi^1.5,无单位,纵坐标为能量差的比值ΔE,单位(%);
需要说明的是,步骤S6~步骤S11:对给定参数的护城河式防洪堤坝渠道原型进行计算流体动力学计算,对护城河式防洪堤坝渠道原型的计算方法,护城河式防洪堤坝渠道原型计算流体动力学计算的网格方案和计算流体动力学软件Fluent设置方式采用单位体积1×1×1的六面体结构化网格,湍流模型为可实现模型Realizable k-ε,控制方程的离散方式为有限体积法(FVM);扩散项采用二阶中心差分格式,对流项采用快速格式(QUICK),压力和速度的耦合采用简单耦合算法(SIMPLEC),计算方式采用并行计算,数值模拟方法采用气液两相流(VOF)方法。仅仅是将计算模型换成护城河式防洪堤坝渠道,其他计算方法不变,步骤S3和步骤S5就是一个选择与确定的过程,主要在于网格尺寸和云计算技术,网格尺寸不随计算对象的放大而放大,避免了由于尺寸放大后存在误差。在护城河式防洪堤坝渠道原型计算流体动力学计算中引入云技术,云计算技术提高了计算精度,也节省了计算的时间成本。
计算后输出原型计算结果,通过计算流体动力学后处理软件CFD-Post处理后得到的相同或不同护城河式防洪堤坝渠道流量Qp条件下数值模拟时,不同几何参数对应护城河式防洪堤坝渠道的能量差的比值ΔE,结果如图5~图10所示;
那么,通过上述步骤确定了护城河式防洪堤坝几何参数,分别为:护城河式防洪堤坝下游护城河深度Dm取值为0.35m,护城河式防洪堤坝下游护城河长度Lm取值为1.50m,护城河式防洪堤坝的下游坡度θd的值取值为1:3,护城河式防洪堤坝的下游坡度θu的值取值为1:1,护城河式防洪堤坝的坝顶宽度Lw的值取值为1.00m,护城河式防洪堤坝的坝顶高度h* L的值取值为1.00m。完毕。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (10)
1.一种通过CFD确定护城河式防洪堤坝几何参数的方法,其特征在于:包括以下步骤:
步骤S1:对缩尺的护城河式防洪堤坝流场进行物理模型试验,建立缩尺的护城河式防洪堤坝渠道物理模型,物理模型试验时,测得护城河式防洪堤坝渠道对称纵断面上、堤坝上游不同过流横断面水深hi;测得护城河式防洪堤坝渠道中护城河内水跃形态及河道中心对称面矢量分布;
步骤S2:对缩尺的护城河式防洪堤坝渠道物理模型进行数值计算,数值模拟时测得护城河式防洪堤坝渠道对称纵断面上、堤坝上游不同过流横断面对应水深Hi;
步骤S3:对数值模拟时对应水深Hi数据集与物理模型试验时水深hi数据集进行逐点对比,找到数据集中所有相对误差绝对值均不超过η%的计算工况,即|Hi-hi|/hi≤η%,输出满足水深模拟精度要求的数值模拟方法;
步骤S4:对步骤S1中缩尺的护城河式防洪堤坝渠道物理模型进行数值计算,采用步骤S3中输出满足水深模拟精度要求的数值模拟方法,输出护城河式防洪堤坝渠道中对应的护城河内水跃形态以及河道中心对称面矢量分布;
步骤S5:对步骤S4中数值模拟时护城河内水跃形态以及河道中心对称面矢量分布中找到与物理模型试验时护城河内水跃形态及河道中心对称面矢量分布相同的数值模拟方案,输出对应的满足护城河内水跃形态及河道中心对称面矢量分布模拟精度要求的数值模拟方法;
步骤S6:确定护城河式防洪堤坝下游护城河深度Dm的值;
步骤S7:确定护城河式防洪堤坝下游护城河长度Lm的值;
步骤S8:确定护城河式防洪堤坝的下游坡度θd的值;
步骤S9:确定护城河式防洪堤坝的上游坡度θu的值;
步骤S10:确定护城河式防洪堤坝的坝顶宽度Lw;
步骤S11:确定护城河式防洪堤坝的堤坝高度hL。
2.根据权利要求1所述的一种通过CFD确定护城河式防洪堤坝几何参数的方法,其特征在于:
步骤S1中建立缩尺的护城河式防洪堤坝渠道物理模型,测得关键参数还包括有;坝顶对称断面处水深Yc,下游尾水断面处水深Yt,坝顶对称断面与下游尾水断面间距M;
其中,护城河式防洪堤坝几何参数为:堤坝高度hL,坝顶宽度Lw,上游坡度θu,下游坡度θd,在堤坝正下游护城河长度为Lm和深度为Dm,上游坝顶拐点到上游坝脚拐点的水平距离为Hu,下游坝顶拐点到下游坝脚拐点的水平距离为Hd,其中上游坡度θu为上游坡度对应的堤坝高度hL与水平距离Hu的比值,下游坡度θd为下游坡度对应的堤坝高度hL与水平距离Hd的比值;
其中,具体参数为:
渠道几何长度为L,渠道宽度为B,护城河式防洪堤坝横跨渠道贯穿布置;
物理模型试验采用正态物理模型,综合考虑物理模型水流在阻力平方区要求,选取物理模型线性比尺λl,物理模型采用重力相似准则设计,角度比λγ,流速比λv=λl 0.5,流量比λQ=λl 2.5,糙率比λn=λl 1/6,护城河式防洪堤坝渠道流量Qp=λQ×Qt。
3.根据权利要求2所述的一种通过CFD确定护城河式防洪堤坝几何参数的方法,其特征在于:
步骤S1中建立缩尺的护城河式防洪堤坝渠道物理模型,物理模型试验具体步骤为:
步骤S101:给定任意模型试验流量Qt,待缩尺的护城河式防洪堤坝渠道中水流稳定后,依次测得护城河式防洪堤坝渠道对称纵断面上、堤坝上游不同过流横断面水深hi,i为不同过流断面的编号;
步骤S102:改变模型试验流量Qt,同时分别改变上游坡度θu,下游坡度θd,下游护城河长度Lm和深度Dm,测得不同工况下坝顶对称断面处水深Yc,下游尾水断面处水深Yt,并通过高速相机摄像得到护城河内水跃形态,以及通过粒子图像测速仪PIV测得河道中心对称面矢量分布;
步骤S103:分别改变上游坡度θu,下游坡度θd,下游护城河长度Lm和深度Dm,具体方法为:
保持下游坡度θd,下游护城河长度Lm和深度Dm不变,仅改变上游坡度θu,或保持上游坡度θu,下游护城河长度Lm和深度Dm不变,仅改变下游坡度θd,或保持上游坡度θu,下游坡度θd,下游护城河深度Dm不变,仅改变下游护城河长度Lm,或保持改变上游坡度θu,下游坡度θd,下游护城河长度Lm不变,仅改变下游护城河深度Dm。
4.根据权利要求3所述的一种通过CFD确定护城河式防洪堤坝几何参数的方法,其特征在于:步骤S6中确定原型护城河式防洪堤坝下游护城河深度Dm的值;具体为:
对护城河式防洪堤坝高度hL,坝顶宽度Lw,上游坡度θu,下游坡度θd,护城河长度Lm在合理范围内进行任意取值,在相同护城河式防洪堤坝渠道流量Qp条件下,分别对a个不同深度的护城河式防洪堤坝渠道进行数值计算,当坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE达到最大时,输出对应的护城河深度Dm为确定值,记作D* m;
护城河深度D* m的具体确定步骤为:
步骤S61:a个不同深度的护城河式防洪堤坝渠道得到a个不同护城河深度Dm,具体值记为χ,即护城河深度值χ1, χ2, χ3,…, χa,其中护城河深度值χ1, χ2, χ3,…, χa逐渐增大,即χ1< χ2< χ3<…< χa,对a个不同深度的护城河式防洪堤坝渠道进行数值计算,分别输出坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE1,ΔE2,ΔE3,…, ΔEa;
步骤S62:当能量差的比值ΔE1=ΔE2=ΔE3=…= ΔEa时,则以Δχ的间距减小χ1分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生减小时,输出前一次数值模拟对应的护城河深度,即第k1次计算时发生能量差的比值ΔE减小,则确定的护城河深度D* m为χ1-(k1-1)×Δχ;
步骤S63:当a个不同护城河深度的具体值χ对应的能量差的比值存在波动,即ΔE1=ΔE2=ΔE3=…=ΔEo<ΔEo+1=…= ΔEa,ΔEo对应的护城河深度为χo,则以Δχ的间距增大χo分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE与ΔEo+1相等时,输出当前数值模拟对应的护城河深度,即第k2次计算得到的能量差的比值ΔE与ΔEo+1相等时,则确定的护城河深度D* m为χo+k2×Δχ;
步骤S64:当a个不同护城河深度的具体值χ对应的能量差的比值单调递增,即ΔE1<ΔE2<ΔE3<…< ΔEa,则以Δχ的间距增大χa分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生相等时,输出前一次数值模拟对应的护城河深度,即第k3次计算时发生能量差的比值ΔE相等,则确定的护城河深度D* m为χa+(k3-1)×Δχ。
5.根据权利要求4所述的一种通过CFD确定护城河式防洪堤坝几何参数的方法,其特征在于:步骤S7中确定原型护城河式防洪堤坝下游护城河长度Lm的值;具体为:
给定步骤S6中确定的护城河深度D* m,对护城河式防洪堤坝高度hL,坝顶宽度Lw,上游坡度θu,下游坡度θd在合理范围内进行任意取值,在相同护城河式防洪堤坝渠道流量Qp条件下,分别对b个不同护城河长度的护城河式防洪堤坝渠道进行数值计算,当坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE达到最大时,输出对应的护城河长度Lm为确定值,记作L* m;
护城河长度L* m的具体确定步骤为:
步骤S71:b个不同长度的护城河式防洪堤坝渠道得到b个不同护城河长度Lm,具体值记为ζ,即护城河长度值ζ1, ζ2, ζ3,…, ζb,其中ζ1, ζ2, ζ3,…, ζb逐渐增大,即ζ1< ζ2< ζ3<…< ζb,对b个不同长度的护城河式防洪堤坝渠道进行数值计算,分别输出坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE1, ΔE2, ΔE3,…, ΔEb;
步骤S72:当能量差的比值ΔE1=ΔE2=ΔE3=…= ΔEb时,则以Δζ的间距减小ζ1分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生减小时,输出前一次数值模拟对应的护城河长度,即第m1次计算时发生能量差的比值ΔE减小,则确定的护城河长度L* m为ζ1-(m1-1)×Δζ;
步骤S73:当b个不同护城河长度的ζ对应的能量差的比值ΔE存在波动,即ΔE1=ΔE2=ΔE3=…=ΔEo<ΔEo+1=…= ΔEb,ΔEo对应的护城河长度为ζo,则以Δζ的间距增大ζo分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE与ΔEo+1相等时,输出当前数值模拟对应的护城河长度,即第m2次计算得到的能量差的比值ΔE与ΔEo+1相等时,则确定的护城河长度L* m为ζo+m2×Δζ;
步骤S74:当b个不同护城河深度的ζ对应的能量差的比值单调递增,即ΔE1<ΔE2<ΔE3<…< ΔEb,则以Δζ的间距增大ζb分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生相等时,输出前一次数值模拟对应的护城河长度,即第m3次计算时发生能量差的比值ΔE相等,则确定的护城河长度L* m为ζb+(m3-1)×Δζ。
6.根据权利要求5所述的一种通过CFD确定护城河式防洪堤坝几何参数的方法,其特征在于:步骤S8中确定原型护城河式防洪堤坝的下游坡度θd的值;具体为:
给定步骤S6中确定的护城河深度D* m和步骤S7中确定的护城河长度L* m,对护城河式防洪堤坝坝高度hL,坝顶宽度Lw,上游坡度θu在合理范围内进行任意取值,在不同护城河式防洪堤坝渠道流量Qp条件下分别对c个不同下游坡度θd的护城河式防洪堤坝渠道进行数值计算,当坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE达到最大时,输出对应的护城河下游坡度θd为确定值,记作θ* d;
护城河下游坡度θ* d的具体确定步骤为:
步骤S81:c个不同下游坡度的护城河式防洪堤坝渠道得到c个不同下游坡度θd,具体值记为п,即护城河下游坡度值п1, п2, п3,…, пc,其中п1, п2, п3,…, пc逐渐增大,即п1<п2< п3<…< пc,其中п1是下游坡度θd允许范围的最小值,пc是下游坡度θd允许范围的最大值,对c个下游坡度的护城河式防洪堤坝渠道进行数值计算,分别输出坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE1, ΔE2, ΔE3,…, ΔEc;
步骤S82:当能量差的比值单调递增时,即ΔE1<ΔE2<ΔE3<…< ΔEc,则输出ΔEc对应的пc为最终的下游坡度值θ* d;
步骤S83:当能量差的比值单调递减时,即ΔE1>ΔE2>ΔE3>…> ΔEc,则输出ΔE1对应的п1为最终的下游坡度值θ* d。
7.根据权利要求6所述的一种通过CFD确定护城河式防洪堤坝几何参数的方法,其特征在于:步骤S9中确定原型护城河式防洪堤坝的上游坡度θu的值;具体为:
给定步骤S6中确定的护城河深度D* m、步骤S7中确定的护城河长度L* m和步骤S8中确定的下游坡度值θ* d,对护城河式防洪堤坝坝高度hL,坝顶宽度Lw在合理范围内进行取值,在不同护城河式防洪堤坝渠道流量Qp条件下分别对e个不同上游坡度θu的护城河式防洪堤坝渠道进行数值计算,当坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE达到最大时,输出对应的护城河上游坡度θu为确定值,记作θ* u;
护城河上游坡度θ* u的具体确定步骤为:
步骤S91:e个不同上游坡度的护城河式防洪堤坝渠道得到e个不同上游坡度θu,具体值记为δ,即护城河上游坡度值δ1, δ2, δ3,…, δe,其中δ1, δ2, δ3,…, δe逐渐增大,即δ1<δ2<δ3<…<δe,其中δ1是上游坡度θu允许范围的最小值,δe是上游坡度θu允许范围的最大值,对e个上游坡度的护城河式防洪堤坝渠道进行数值计算,分别输出坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE1, ΔE2, ΔE3,…, ΔEe;
步骤S92:当能量差的比值单调递增时,即ΔE1<ΔE2<ΔE3<…<ΔEe,则输出ΔEe对应的δe为最终的下游坡度值θ* u;
步骤S93:当能量差的比值单调递减时,即ΔE1>ΔE2>ΔE3>…>ΔEe,则输出ΔE1对应的δ1为最终的下游坡度值θ* u。
8.根据权利要求7所述的一种通过CFD确定护城河式防洪堤坝几何参数的方法,其特征在于:步骤S10中确定原型护城河式防洪堤坝的坝顶宽度Lw;具体为:
给定步骤S6中确定的护城河深度D* m、步骤S7中确定的护城河长度L* m、步骤S8中确定的下游坡度值θ* d和步骤S9中确定的上游坡度值θ* u,对护城河式防洪堤坝坝高度hL在合理范围内进行取值,在不同护城河式防洪堤坝渠道流量Qp条件下分别对f个不同坝顶宽度Lw的护城河式防洪堤坝渠道进行数值计算,当坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE达到最大时,输出对应的坝顶宽度Lw为确定值,记作L* w;
防洪堤坝的坝顶宽度L* w的具体确定步骤为:
步骤S111:f个不同坝顶宽度的护城河式防洪堤坝渠道得到f个不同坝顶宽度Lw,具体值记为p,即护城河上游坡度值p1, p2, p3,…, pf,其中p1, p2, p3,…, pf逐渐增大,即p1<p2<p3<…<pf,其中p1是坝顶宽度Lw允许范围的最小值,pf是坝顶宽度Lw允许范围的最大值,对f个坝顶宽度的护城河式防洪堤坝渠道进行数值计算,分别输出坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE1, ΔE2, ΔE3,…, ΔEf;
步骤S112:当能量差的比值ΔE1=ΔE2=ΔE3=…=ΔEf时,则以Δp的间距减小p1分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生减小时,输出前一次数值模拟对应的护城河长度,即第j1次计算时发生能量差的比值ΔE减小,则确定的防洪堤坝的坝顶宽度L* w为p1-(j1-1)×Δp;
步骤S113:当o个不同坝顶宽度的p对应的能量差的比值单调递增,而剩余的j-o个不同坝顶宽度的p对应的能量差的比值相等,假设其中一种情况为ΔE1<ΔE2<ΔE3<…<ΔEo<ΔEo+1=…= ΔEf,则以Δp的间距减小po+1分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生减小时,输出前一次数值模拟对应的护城河长度,即第j2次计算时发生能量差的比值ΔE减小,则确定的防洪堤坝的坝顶宽度L* w为po+1-(j2-1)×Δp;
步骤S114:当f个不同长度的p对应的能量差的比值单调递增,即ΔE1<ΔE2<ΔE3<…<ΔEf,则以Δp的间距增大pf分别进行数值模拟输出能量差的比值ΔE,当能量差的比值ΔE发生相等时,输出前一次数值模拟对应的护城河长度,即第j3次计算时发生能量差的比值ΔE相等,则确定的防洪堤坝的坝顶宽度L* w为pf+(j3-1)×Δp。
9.根据权利要求8所述的一种通过CFD确定护城河式防洪堤坝几何参数的方法,其特征在于:步骤S11:确定原型护城河式防洪堤坝的堤坝高度hL;具体为:
给定步骤S6中确定的护城河深度D* m、步骤S7中确定的护城河长度L* m、步骤S8中确定的下游坡度值θ* d、步骤S9中确定的上游坡度值θ* u和步骤S10中确定的坝顶宽度L* w,在不同护城河式防洪堤坝渠道流量Qp条件下分别对s个不同堤坝高度hL的护城河式防洪堤坝渠道进行数值计算,当坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE达到最大时,输出对应的堤坝高度hL为确定值,记作h* L;
防洪堤坝的堤坝高度h* L的具体确定步骤为:
步骤S115: s个不同堤坝高度的护城河式防洪堤坝渠道得到s个不同堤坝高度hL,具体值记为,即护城河上游坡度值/>, />, />,…, />,其中/>, />, />,…, />逐渐增大,即</></><…</>,其中/>是堤坝高度hL允许范围的最小值,/>是堤坝高度hL允许范围的最大值,对s个堤坝高度的护城河式防洪堤坝渠道进行数值计算,分别输出坝顶对称断面与下游尾水断面之间能量差与坝顶对称断面总能量的比值ΔE1, ΔE2, ΔE3,…, ΔEs;
步骤S116:当能量差的比值单调递增时,即ΔE1<ΔE2<ΔE3<…<ΔEs,则输出ΔEs对应的为确定防洪堤坝的堤坝高度,记作h* L;
步骤S117:当能量差的比值单调递减时,即ΔE1>ΔE2>ΔE3>…>ΔEc,则输出ΔE1对应的为确定防洪堤坝的堤坝高度,记作h* L。
10.根据权利要求9所述的一种通过CFD确定护城河式防洪堤坝几何参数的方法,其特征在于:步骤S6-步骤S11中,坝顶对称断面A-A与下游尾水断面Z-Z之间能量差与坝顶对称断面总能量的比值ΔE,计算方法为:
先计算坝顶对称断面A-A流速VA为渠道流量Q除以坝顶对称断面A-A过流面积,坝顶对称断面A-A过流面积为坝顶对称断面处水深Yc与渠道宽度B的乘积,则VA=Q/(Yc×B),那么坝顶对称断面A-A的总能量为EA=VA²/2g+Yc+hL;
再计算下游尾水断面Z-Z流速VZ为渠道流量Q除以下游尾水断面Z-Z过流面积,下游尾水断面Z-Z过流面积为下游尾水断面处水深Yt与渠道宽度B的乘积,则VZ= Q/(Yt×B),那么下游尾水断面Z-Z的总能量为EZ=VZ²/2g+Yt;
最后计算坝顶对称断面A-A与下游尾水断面Z-Z之间能量差与坝顶对称断面总能量的比值,简称能量差的比值ΔE,ΔE=(EA-EZ)/ EA。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311786133.1A CN117454724B (zh) | 2023-12-25 | 2023-12-25 | 一种通过cfd确定护城河式防洪堤坝几何参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311786133.1A CN117454724B (zh) | 2023-12-25 | 2023-12-25 | 一种通过cfd确定护城河式防洪堤坝几何参数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117454724A true CN117454724A (zh) | 2024-01-26 |
CN117454724B CN117454724B (zh) | 2024-03-19 |
Family
ID=89591399
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311786133.1A Active CN117454724B (zh) | 2023-12-25 | 2023-12-25 | 一种通过cfd确定护城河式防洪堤坝几何参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117454724B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117829036A (zh) * | 2024-03-05 | 2024-04-05 | 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) | 一种通过cfd确定t型分汊明渠下游支渠流量的方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20120112315A (ko) * | 2012-08-30 | 2012-10-11 | 이우성 | 물의 위치에너지를 이용한 수해조절장치 시스템 |
EP3118375A1 (en) * | 2015-07-13 | 2017-01-18 | Dura Vermeer Beton & Waterbouw B.V. | Method for improving an inwards stability of a levee |
WO2017217957A1 (en) * | 2016-06-13 | 2017-12-21 | Schlumberger Technology Corporation | Runtime parameter selection in simulations |
US20200311319A1 (en) * | 2019-03-28 | 2020-10-01 | China Waterborne Transport Research Institute | Method for evaluating ecological environmental impact of channel project and countermeasures thereof based on mechanism analysis |
CN112818441A (zh) * | 2020-12-31 | 2021-05-18 | 田志强 | 一种江洪灾害堤坝分洪道的设计方法 |
CN114004003A (zh) * | 2021-12-31 | 2022-02-01 | 华南理工大学 | 一种适用于城区复杂下垫面的水库溃坝洪水数值模拟方法 |
CN114117965A (zh) * | 2021-12-07 | 2022-03-01 | 南昌工程学院 | 一种基于两侧半圆柱矩形槽比降效应的灌区渠道测流法 |
KR20230055044A (ko) * | 2021-10-18 | 2023-04-25 | 한국수자원공사 | 1차원 하천흐름 해석기술 |
US20230169243A1 (en) * | 2021-11-29 | 2023-06-01 | China Institute Of Water Resources And Hydropower Research | Numerical Simulation Method for the Overflowing of River Channel Structures Based on Finite Volume Method |
CN116704357A (zh) * | 2023-08-09 | 2023-09-05 | 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) | 基于YOLOv7的堤坝边坡滑坡智能识别与预警方法 |
CN116757111A (zh) * | 2023-06-05 | 2023-09-15 | 南昌工程学院 | 一种预测带有鱼道坡度的墩块式鱼道流量的方法 |
-
2023
- 2023-12-25 CN CN202311786133.1A patent/CN117454724B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20120112315A (ko) * | 2012-08-30 | 2012-10-11 | 이우성 | 물의 위치에너지를 이용한 수해조절장치 시스템 |
EP3118375A1 (en) * | 2015-07-13 | 2017-01-18 | Dura Vermeer Beton & Waterbouw B.V. | Method for improving an inwards stability of a levee |
WO2017217957A1 (en) * | 2016-06-13 | 2017-12-21 | Schlumberger Technology Corporation | Runtime parameter selection in simulations |
US20200311319A1 (en) * | 2019-03-28 | 2020-10-01 | China Waterborne Transport Research Institute | Method for evaluating ecological environmental impact of channel project and countermeasures thereof based on mechanism analysis |
CN112818441A (zh) * | 2020-12-31 | 2021-05-18 | 田志强 | 一种江洪灾害堤坝分洪道的设计方法 |
KR20230055044A (ko) * | 2021-10-18 | 2023-04-25 | 한국수자원공사 | 1차원 하천흐름 해석기술 |
US20230169243A1 (en) * | 2021-11-29 | 2023-06-01 | China Institute Of Water Resources And Hydropower Research | Numerical Simulation Method for the Overflowing of River Channel Structures Based on Finite Volume Method |
CN114117965A (zh) * | 2021-12-07 | 2022-03-01 | 南昌工程学院 | 一种基于两侧半圆柱矩形槽比降效应的灌区渠道测流法 |
CN114004003A (zh) * | 2021-12-31 | 2022-02-01 | 华南理工大学 | 一种适用于城区复杂下垫面的水库溃坝洪水数值模拟方法 |
CN116757111A (zh) * | 2023-06-05 | 2023-09-15 | 南昌工程学院 | 一种预测带有鱼道坡度的墩块式鱼道流量的方法 |
CN116704357A (zh) * | 2023-08-09 | 2023-09-05 | 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) | 基于YOLOv7的堤坝边坡滑坡智能识别与预警方法 |
Non-Patent Citations (3)
Title |
---|
STEVEN A 等: "Renormalization group analysis of turbulence", JOURNAL OF SCIENTIFIC COMPUTING, 31 December 1986 (1986-12-31) * |
曾丹;刘成林;陈宇豪;: "弯道溃坝水流的三维数值模拟", 南昌大学学报(工科版), no. 04, 28 December 2016 (2016-12-28) * |
贺治国;: "考虑水沙相互作用的河堤溃决二维数值模型", 水动力学研究与进展A辑, no. 02, 15 March 2010 (2010-03-15) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117829036A (zh) * | 2024-03-05 | 2024-04-05 | 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) | 一种通过cfd确定t型分汊明渠下游支渠流量的方法 |
CN117829036B (zh) * | 2024-03-05 | 2024-05-14 | 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) | 一种通过cfd确定t型分汊明渠下游支渠流量的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN117454724B (zh) | 2024-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111177875B (zh) | 一种河道整治方案模拟与优化方法 | |
CN117454724B (zh) | 一种通过cfd确定护城河式防洪堤坝几何参数的方法 | |
CN109271672A (zh) | 一种河-湖-泵站相互影响作用下的河道水面线计算方法 | |
CN104933268A (zh) | 一种基于一维非恒定流数值模型的洪水分析方法 | |
CN111241757A (zh) | 基于计算流体力学的铀尾矿库溃坝三维数值模拟方法 | |
CN117408188B (zh) | 一种通过计算流体动力学预测梯形溢流坝流量的方法 | |
Wang et al. | Three-dimensional flood routing of a dam break based on a high-precision digital model of a dense urban area | |
CN113742820B (zh) | 一种分洪隧洞出口消力池消能效果分析的数值模拟方法 | |
CN113297753B (zh) | 山区河流透过性坝沿程水深预测方法 | |
CN116384289B (zh) | 一种通过计算流体动力学预测墩块式鱼道流量的方法 | |
CN112182814A (zh) | 一种基于稀疏断面点据的河道水下地形建模方法 | |
Daneshfaraz et al. | Numerical modeling of hydraulic properties of sloped broad crested weir | |
Erpicum et al. | Numerical and physical hydraulic modelling of Piano Key Weirs | |
CN114925624A (zh) | 一种天然河道三维水流数值模拟方法 | |
Elsayed et al. | Evaluation the Existing Nile River Navigation Path in The Reach from Aswan City to Esna Barrage | |
Balas et al. | Hydrodynamic Modelling Using HİDROTÜRK Model | |
Abidi | Three-Dimensional Modeling of Pressure and Velocity Distribution in the Wide Edge Side Spillway Located in the Trapezoidal Channel | |
CN112433029B (zh) | 一种滩地树木糙率的计算方法 | |
CN116644499B (zh) | 适用于运河浪岸坡击区作用范围的确定方法及船行波模型试验装置 | |
Kuiry et al. | Application of the 1d-quasi 2d model tinflood for floodplain inundation prediction of the river thames | |
Erpicum et al. | The Taoussa Project (Mali): an example of effective composite modeling | |
Krylova et al. | Reproduction of runoff hydrograph in the Lena River basin with a hydrologically correct digital elevation model | |
Zamiri et al. | Studying the effect of shape changes in plan of labyrinth weir on increasing flow discharge coefficient using Flow-3D numerical model | |
Morikawa et al. | Numerical analysis of flood with a double grid model | |
Zhanga et al. | Influence of bridge piers shapes on the flow of the lower Yellow River |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |