CN117408188B - 一种通过计算流体动力学预测梯形溢流坝流量的方法 - Google Patents
一种通过计算流体动力学预测梯形溢流坝流量的方法 Download PDFInfo
- Publication number
- CN117408188B CN117408188B CN202311718811.0A CN202311718811A CN117408188B CN 117408188 B CN117408188 B CN 117408188B CN 202311718811 A CN202311718811 A CN 202311718811A CN 117408188 B CN117408188 B CN 117408188B
- Authority
- CN
- China
- Prior art keywords
- trapezoid
- dam
- overflow dam
- flow
- overflow
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- 239000012530 fluid Substances 0.000 title claims abstract description 49
- 238000011144 upstream manufacturing Methods 0.000 claims abstract description 103
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 80
- 230000014509 gene expression Effects 0.000 claims abstract description 53
- 238000004364 calculation method Methods 0.000 claims abstract description 39
- 238000012360 testing method Methods 0.000 claims abstract description 34
- 238000004088 simulation Methods 0.000 claims abstract description 27
- 241000826860 Trapezium Species 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 230000000149 penetrating effect Effects 0.000 claims description 3
- 238000002620 method output Methods 0.000 claims description 2
- 230000011218 segmentation Effects 0.000 claims description 2
- 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
- 230000002262 irrigation Effects 0.000 description 4
- 238000003973 irrigation Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 239000003621 irrigation water Substances 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 230000005514 two-phase flow Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001010 compromised effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000005457 optimization Methods 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/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
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- 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
- 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)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Measuring Volume Flow (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种通过计算流体动力学预测梯形溢流坝流量的方法,基于缩尺的梯形溢流坝进行物理模型试验和数值计算,得到对应的数值模拟方法;对输出的上游坡度、下游坡度、相对坝顶宽度进行处理得到梯形溢流坝流量系数表达式,从而得到梯形溢流坝流量仅随梯形溢流坝上游总能量水头变化的表达式。在数值计算时,引入云计算以提高模拟精度和效率。在预测梯形溢流坝流量时,测得梯形溢流坝上游总能量水头或上游水深带入对应流量公式即可得到实时的梯形溢流坝流量。本发明能够准确合理地预测设置了梯形溢流坝的渠道流量。
Description
技术领域
本发明属于流体动力学领域,具体涉及一种通过计算流体动力学预测梯形溢流坝流量的方法。
背景技术
在国内各类渠道中梯形溢流坝通常用作灌溉结构,具有短坝顶的梯形剖面的堰用于将河水引向灌溉渠用于灌溉庄稼和农作物。通常在梯形溢流坝装置中配有控制流量的闸门,可以根据作物灌溉水量的需要来调节上游水位,从而达到灌溉要求。此外,此类梯形溢流坝还可以根据闸门开度和梯形溢流坝上游及下游水位来充当流量测量装置。然而有研究表明,在测量流量时流场的改变决定了流量的预测精度。通常随着梯形溢流坝上、下游水位相对位置的变化,梯形溢流坝的下游水流流场在自由出流和淹没出流两种形态中切换。通常发生洪水或大流量工况时,梯形溢流坝下游流场处于淹没出流状态,而日常小流量工况使用,梯形溢流坝下游流场处于自由出流状态。
梯形溢流坝的流场除了受下游出流状态影响之外,还受梯形溢流坝几何构造的影响,主要包括梯形溢流坝的坝顶宽度、上游坡度和下游坡度。虽然有部分学者有针对地进行了一些研究,但是均只考虑某一个因素,或者基于一些古老的公式进行了一些不太符合逻辑的修改,这导致过去学者们提出的流量预测公式缺乏科学性和系统性,使得流量预测精度大打折扣。需要指出的是,本技术在预研究中发现,上游坡度和下游坡度的临界点均在1:1,换言之,上游坡度或下游坡度在缓于1:1或陡于1:1时,梯形溢流坝中流场均发生较大的改变,这正是过去研究者没有考虑的,他们均将所有坡度综合考虑,包括上游坡度和下游坡度为直角的情况(均为直角时为宽顶堰或薄壁堰),这导致流量预测的精度大受影响。本技术受经费和时间的限制,仅考虑了上游坡度和下游坡度缓于1:1的情况。此外,梯形溢流坝的坝顶宽度也是一个必须要考虑的重要因素,这也是很多研究者没有考虑的注意点,或许仅针对某一个固定坝顶宽度的梯形溢流坝提出一个流量预测公式是可行的,但是换一个不同坝顶宽度的梯形溢流坝必将不适用。本技术为了提高所提出的流量预测公式的普适性,将梯形溢流坝的坝顶宽度作为一个关键因素考虑在流量预测方法中,从而拓宽所提出梯形溢流坝流量预测方法的应用场景。
随着计算流体动力学(CFD)和云计算技术的快速发展,由于其具有建模便捷、可直接对原型进行模拟、试验重复性较好等优点,数值模拟方法逐渐成为一种非常有效的研究手段。但是现有技术中还没有利用数值试验手段精确预测梯形溢流坝流量的方法。中国专利公开号CN106844913B中公开了一种基于三维计算流体动力学(CFD)的滞留气团热力学特性模拟方法,该专利利用三维计算流体动力学(CFD)进行了滞留气团的热力学分析,但是其方法不能用于预测梯形溢流坝的流量。因此,如何利用计算流体动力学(CFD)方法来准确预测梯形溢流坝的流量是一个非常值得研究的课题,它将有效指导实际工程中梯形溢流坝的设计与优化,并显著提高水利管理部门的决策效率。本技术主要针对小流量工况,即梯形溢流坝下游处于自由出流状态的流量预测进行研究,从而提出一种适用于自由出流条件的梯形溢流坝的流量预测办法。
发明内容
本发明的目的:为了克服现有技术中存在的不足,本发明提供一种通过计算流体动力学预测梯形溢流坝流量的方法,能够简便、准确和合理地预测梯形溢流坝的流量,有效地保证了梯形溢流堤的经济、安全和稳定运行。
本发明的技术方案:一种通过计算流体动力学预测梯形溢流坝流量的方法,包括以下步骤:
步骤S1:对缩尺的梯形溢流坝进行物理模型试验,建立梯形溢流坝渠道物理模型;在模型试验流量Qt条件下测得物理模型试验时,沿梯形溢流坝渠道对称纵断面上、梯形溢流坝上游不同过流断面水深hi;
步骤S2:对步骤S1中建立的梯形溢流坝渠道物理模型进行数值计算,计算收敛后导出结果,通过计算流体动力学后处理软件进行后处理;在与步骤S1几何尺寸相同的梯形溢流坝的条件下,计算得到相同模型试验流量Qt条件下测得数值模拟时,沿梯形溢流坝渠道对称纵断面上、梯形溢流坝上游不同过流横断面对应水深Hi;
步骤S3:对步骤S2中数值模拟时对应水深Hi数据集与步骤S1中物理模型试验时对应水深hi数据集进行逐点对比,找到所有数值逐点相对误差绝对值均不超过η%的计算工况,即|Hi-hi|/hi≤η%,输出对应的数值模拟方法;
步骤S4:根据梯形溢流坝的流量Q表达式Q=Cdb(2gH0 3)1/2得到梯形溢流坝流量系数Cd的表达式;其中b为渠道宽度,g为重力加速度,H0为梯形溢流坝上游总能量水头;A、B、C、k和/>为常数系数,ξ为相对坝顶宽度,θu表示上游坡度,θd表示下游坡度;
步骤S5:对不同几何条件的梯形溢流坝进行计算流体动力学计算,采用统计产品与服务解决方案SPSS软件对相对坝顶宽度ξ,上游坡度θu,下游坡度θd和梯形溢流坝流量系数Cd进行非线性拟合,得到常数系数A、B、C、k和的具体值记为A*、B*、C*、k*和/>,将常数系数具体值反代入梯形溢流坝流量系数Cd的表达式,得到梯形溢流坝流量系数表达式;
步骤S6:将步骤S5得到的梯形溢流坝流量系数表达式代入步骤S4梯形溢流坝的流量Q表达式,得到梯形溢流坝的流量Q的表达式为:
;
将坝顶宽度Lc、坝顶高度hL、上游坡度θu和下游坡度θd安置在梯形溢流坝给定渠道宽度b中,实时测得梯形溢流坝上游总能量水头H0,求得梯形溢流坝中的实时流量;或者测得梯形溢流坝上游水深Hp求解可得梯形溢流坝的流量Q。
进一步的,步骤S1中对缩尺的梯形溢流坝进行物理模型试验,应保证梯形溢流坝的下游水流状态始终为自由出流,具体方法为:
在梯形溢流坝的上游逐渐增加模型试验流量Qt,在梯形溢流坝下游设置尾门并控制尾门的开度,同时监测梯形溢流坝上游水位ht和下游水位Tt,随着下游水位Tt逐渐上升,始终保证梯形溢流坝上游水位ht变化明显不超过1mm。
进一步的,步骤S1中缩尺的梯形溢流坝渠道物理模型,具体参数为:
渠道几何长度为L,梯形溢流坝渠道横跨渠道贯穿布置;
梯形溢流坝几何参数为:坝顶高度hL,坝顶宽度Lc,上游坝顶拐点到上游坝脚拐点的水平距离为Hu,下游坝顶拐点到下游坝脚拐点的水平距离为Hd,其中上游坡度为上游坡度对应的坝顶高度hL与上游坝顶拐点到上游坝脚拐点的水平距离Hu的比值,下游坡度为下游坡度对应的坝顶高度hL与下游坝顶拐点到下游坝脚拐点的水平距离Hd的比值;
上游坡度和下游坡度的陡度均不低于1:1;
物理模型试验采用正态物理模型,综合考虑物理模型水流在阻力平方区要求,选取物理模型线性比尺λl,物理模型采用重力相似准则设计,角度比λγ,流速比λv=λl 0.5,流量比λQ=λl 2.5,糙率比λn=λl 1/6,梯形溢流坝渠道流量Qp=λQ×Qt。
进一步的,步骤S1中物理模型试验时沿梯形溢流坝渠道对称纵断面上、梯形溢流坝上游不同过流断面水深hi,i为随机断面编号;对应的步骤S2中数值模拟时沿梯形溢流坝渠道对称纵断面上、梯形溢流坝上游不同过流横断面对应水深Hi ,i断面编号与步骤S1中相同。
进一步的,步骤S2中梯形溢流坝进行数值计算,具体为:
建立与步骤S1中梯形溢流坝渠道相同的物理模型,并将物理模型进行网格剖分,输出后缀名为.mesh的计算文件,导入计算流体动力学软件进行数值计算;
在将物理模型进行网格剖分时采用不同尺度网格,得到若干不同的网格方案,若干不同的网格方案在计算流体动力学软件计算中选取不同的湍流模型及不同的数值算法。
进一步的,步骤S4中:H0的表达式:H0=h+v2/2g;h为梯形溢流坝的坝上水头,v为渠道中平均流速;
渠道中平均流速v的表达式为v=Q/(h+hL)b,h+ hL为梯形溢流坝上游水深Hp;
ξ的表达式为H0/(H0+Lc);上游坡度θu的表达式为θu=atcan(hL/Hu);下游坡度θd的表达式为θd=atcan(hL/Hd),θ单位为rad。
进一步的,步骤S5中:对不同几何条件的梯形溢流坝进行计算流体动力学计算,具体为:
给定坝顶高度hL的梯形溢流坝,在任意上游坡度θu、下游坡度θd和坝顶宽度Lc条件下,输出梯形溢流坝渠道流量Qp条件下数值模拟时梯形溢流坝上游水深Hp以确定相对坝顶宽度ξ,再根据步骤S4梯形溢流坝的流量Q表达式Q=Cdb(2gH0 3)1/2求得梯形溢流坝流量系数Cd;即每一个上游坡度θu、下游坡度θd、坝顶宽度Lc和梯形溢流坝渠道流量Qp均对应一个相对坝顶宽度ξ和溢流坝流量系数Cd。
进一步的,步骤S5中:对梯形溢流坝进行计算流体动力学计算,梯形溢流坝的上游坡度θu和下游坡度θd的陡度均不低于1:1,并且梯形溢流坝的下游水流状态始终为自由出流。
进一步的,步骤S5在梯形溢流坝渠道进行计算流体动力学计算中采用步骤S3中输出的对应的数值模拟方法。
进一步的,步骤S6中:测得梯形溢流坝上游水深Hp,求解可得梯形溢流坝的流量Q,具体步骤为:
将步骤S5得到的梯形溢流坝流量系数表达式代入步骤S4梯形溢流坝的流量Q表达式,得到梯形溢流坝的流量Q的表达式为:
;
将坝顶宽度Lc、坝顶高度hL、上游坡度θu和下游坡度θd的梯形溢流坝安置在给定渠道宽度b中,测得梯形溢流坝上游水深Hp,转换成梯形溢流坝上游总能量水头H0,即H0=(Q/b/Hp)2/2g+(Hp-hL),再将中H0全部替换为(Q/b/Hp)2/2g+(Hp-hL),则/>的等号两侧将变为只有梯形溢流坝的流量Q未知的非线性表达式,求解可得梯形溢流坝的流量。
本发明的有益效果:本发明通过计算流体动力学预测梯形溢流坝流量的方法,能够简便、准确和合理地预测梯形溢流坝的实时流量,优化了水利工程中布置梯形溢流坝的规划合理性,显著提高了水利工程管理部门防洪排涝的决策效率,有效保证了灌区作物的灌溉水保证率。
附图说明
图1为梯形溢流坝渠道结构侧视图示意图;
图2为梯形溢流坝渠道结构平面图示意图;
图3为物理试验和数值模拟中梯形溢流坝上游渠道纵对称断面上水深对比图;
图4为下游坡度为1:1且上游坡度为1:1, 1:2和1:3条件下流量系数Cd与相对坝顶宽度ξ的关系;
图5为下游坡度为1:2且上游坡度为1:1, 1:2和1:3条件下流量系数Cd与相对坝顶宽度ξ的关系;
图6为下游坡度为1:3且上游坡度为1:1, 1:2和1:3条件下流量系数Cd与相对坝顶宽度ξ的关系。
具体实施方式
下面结合附图对本发明作更进一步的说明。
本发明的一种通过计算流体动力学预测梯形溢流坝流量的方法,包括先对数值模拟的方法进行验证,通过对梯形溢流坝渠道缩尺后的物理模型进行数值计算和物理试验,并对梯形溢流坝渠道纵对称断面上、梯形溢流坝上游不同过流断面水深进行采集与对比,得到合理的计算流体动力学数值方案包括网格剖分、湍流模型和数值方法;再通过上述数值方案对原型条件下的梯形溢流坝渠道进行计算流体动力学计算并对提出的流量系数表达式进行非线性拟合得到准确的表达式常数系数,最后通过实时测得的原型梯形溢流坝上游水深预测出原型梯形溢流坝的实时流量。本发明通过缩尺物理模型试验和相应的物理模型计算流体动力学数值计算进行对比,用物理模型试验来约束数值计算方法,再通过得出的数值计算方法对原型条件下的梯形溢流坝进行计算以确定原型条件下梯形溢流坝的上游水深,如此使得原型条件下的数值计算的结果更科学有效;使原型条件下构建的计算流体动力学数值计算模型(数值计算模型是指数值模拟构建的模型)能够满足实际流动条件,提高了计算流体动力学数值模拟的准确性,能够准确地预测梯形溢流坝的实时流量。
实施例
本实施例的通过计算流体动力学预测梯形溢流坝流量的方法,具体步骤如下:
步骤S1:对缩尺的梯形溢流坝进行物理模型试验,建立梯形溢流坝渠道物理模型;
其缩尺的梯形溢流坝渠道几何长度L为10m,渠道宽度b为0.3m,缩尺的梯形溢流坝横跨渠道贯穿布置,在模型试验流量Qt为42L/s的条件下,测试沿梯形溢流坝渠道对称纵断面上、梯形溢流坝上游不同过流横断面对应水深hi,i为随机断面编号,此时梯形溢流坝的几何参数为:上游坝顶拐点到上游坝脚拐点的水平距离Hu为200mm,则上游坡度θu为1:2,下游坝顶拐点到下游坝脚拐点的水平距离Hd为200mm则下游坡度θd为1:2,坝顶宽度Lc为40mm和坝顶高度hL为100mm,物理模型试验时测得监测梯形溢流坝上游水位ht为0.3m(梯形溢流坝渠道结构如图1和图2所示);
物理试验模型采用正态模型,综合考虑物理试验模型水流在阻力平方区要求,拟选取模型线性比尺为10,物理模型采用重力相似准则设计,角度比为1,流速比为λv=100.5=3.16,流量比λQ=102.5=316.23,糙率比λn=101/6=1.47。
在梯形溢流坝的上游逐渐增加模型试验流量至42L/s,在梯形溢流坝下游设置尾门并控制尾门的开度,同时监测梯形溢流坝上游水位ht和下游水位Tt,随着下游水位Tt逐渐上升,始终保证梯形溢流坝上游水位ht变化明显不超过1mm,以保证梯形溢流坝的下游水流状态始终为自由出流。
步骤S2:对步骤S1中缩尺的梯形溢流坝渠道物理模型进行数值计算,计算收敛后导出结果,通过计算流体动力学后处理软件进行后处理;在与步骤S1几何尺寸相同的梯形溢流坝的条件下,计算得到相同模型试验流量Qt条件下测得数值模拟时沿梯形溢流坝渠道对称纵断面上、梯形溢流坝上游不同过流横断面对应水深Hi,i与步骤S1中相对应;
数值计算时,先对步骤S1中梯形溢流坝渠道模型进行数值计算,建立梯形溢流坝渠道三维数值模型,并将梯形溢流坝渠道三维数值模型进行网格剖分,输出后缀名为.mesh的计算文件,再导入计算流体动力学软件Fluent进行数值计算,最后计算收敛后导出结果;
通过计算流体动力学后处理软件CFD-Post进行后处理,在上游坡度θu为1:2、下游坡度θd为1:2、坝顶宽度Lc为40mm、在与相同模型试验流量Qt =42L/s的条件下条件进行计算获得沿梯形溢流坝渠道对称纵断面上、梯形溢流坝上游不同过流横断面对应水深Hi;
再将三维物理模型网格剖分时采用不同尺度网格,从而得到若干不同的网格方案,针对各个网格方案在计算流体动力学软件Fluent计算中选取不同的湍流模型及不同的数值算法,形成若干不同的数值模拟方法,那么通过不同数值模拟方法对缩尺的梯形溢流坝计算可得到不同的Hi数据集。
步骤S3:对步骤S2中数值模拟时对应水深Hi数据集与步骤S1中物理模型试验时对应水深hi数据集进行逐点对比,找到所有数值逐点相对误差绝对值均不超过η%的计算工况,即|Hi-hi|/hi≤2.5%,输出对应的数值模拟方法;图3为物理试验和数值模拟中梯形溢流坝上游渠道对称纵断面上水深对比图,其相对误差的绝对值均小于2.5%,输出对应的数值模拟方法,包括网格尺寸、湍流模型、数值算法,拟采用单位体积1×1×1的六面体结构化网格,湍流模型为可重整化群(RNG k-ε),控制方程的离散方式为有限体积法(FVM);扩散项采用二阶中心差分格式,对流项采用快速格式(QUICK),压力和速度的耦合采用简单耦合算法(SIMPLEC),计算方式采用并行计算,数值模拟方法采用气液两相流(VOF)方法。
步骤S4:根据梯形溢流坝的流量Q表达式Q=Cdb(2gH0 3)1/2得到梯形溢流坝流量系数Cd的表达式;其中b为渠道宽度,g为重力加速度,H0为梯形溢流坝上游总能量水头;A、B、C、k和/>为常数系数,ξ为相对坝顶宽度,θu表示上游坡度,θd表示下游坡度;
H0的表达式:H0=h+v2/2g;h为梯形溢流坝的坝上水头,v为渠道中平均流速;
渠道中平均流速v的表达式为v=Q/(h+hL)b,h+ hL为梯形溢流坝上游水深Hp;
ξ的表达式为H0/(H0+Lc);上游坡度θu的表达式为θu=atcan(hL/Hu);下游坡度θd的表达式为θd=atcan(hL/Hd),θ单位为rad。
步骤S5:对不同几何条件的梯形溢流坝进行计算流体动力学计算,梯形溢流坝的上游坡度θu和下游坡度θd的陡度均不低于1:1,并且梯形溢流坝的下游水流状态始终为自由出流,拟对给定坝顶高度hL=1m的梯形溢流坝,在上游坡度为1:1,1:2和1:3、下游坡度为1:1,1:2和1:3和坝顶宽度Lc=1m,2m,3m,4m和5m条件下分别进行组合,采用单位体积1×1×1的六面体结构化网格,湍流模型为可重整化群(RNG k-ε),控制方程的离散方式为有限体积法(FVM);扩散项采用二阶中心差分格式,对流项采用快速格式(QUICK),压力和速度的耦合采用简单耦合算法(SIMPLEC),计算方式采用并行计算,数值模拟方法采用气液两相流(VOF)方法,对梯形溢流坝进行计算流体动力学计算,输出梯形溢流坝渠道流量Qp=2m³/s,3m³/s, 5 m³/s, 7 m³/s, 9 m³/s, 11 m³/s, 12 m³/s, 14 m³/s和16 m³/s条件下数值模拟时梯形溢流坝上游水深Hp以确定原型条件下相对坝顶宽度ξ,然后根据步骤S4中梯形溢流坝流量Q表达式Q=Cdb (2gH0 3)1/2求得相应的流量系数Cd,如图4~图6所示为不同上游坡度θu和下游坡度θd条件下流量系数Cd与相对坝顶宽度ξ的关系,再结合任意上游坡度θu、下游坡度θd采用统计产品与服务解决方案SPSS软件对梯形溢流坝流量系数的表达式进行非线性拟合,得到常数系数A、B、C、k和/>的具体值,记为0.457、0.082、0.720、0.010和0.366,然后将系数具体值0.457、0.082、0.720、0.010和0.366反带入流量系数的表达式得到流量系数表达式Cd=0.457+0.082sin[π(ξ-0.720)+( 0.010θu+0.366θd)];
统计产品与服务解决方案的数据处理软件(Statistical Product and ServiceSolutions,SPSS)具体拟合步骤为:
1、先将数据导入数据处理软件统计产品与服务解决方案的数据处理软件(SPSS);
2、然后进行回归拟合,在分析中找到回归,然后进行非线性拟合,确定因变量和模型表达式;
3、在选型中保存预存初始值,勾选残差值,拟合后即可得到合理的常数系数A、B、C、k和的具体值。
步骤S6:将步骤S5得到的梯形溢流坝流量系数表达式代入步骤S4梯形溢流坝的流量Q表达式,得到梯形溢流坝的流量Q的表达式为:
Q={0.457+0.082sin[π((H0/(H0+Lc)-0.720)+( 0.010θu+0.366θd)]}×b×(2gH0 3)1/2;
假设将坝顶宽度Lc=1m、上游坡度和下游坡度均为1:3和坝顶高度hL=1m的梯形溢流坝安置在渠道宽度为3m的渠道中,则θu=θd=atctan(hL/Hu)=atctan1/3=0.32175以及坝顶宽度Lc=1m的梯形溢流坝,代入梯形溢流坝的流量Q的表达式简化得到该梯形溢流坝的流量Q随梯形溢流坝上游总能量水头H0变化的表达式Q=3×{0.457+0.082sin[π((H0/(H0+1)-0.720)+0.012]}× (2gH0 3)1/2;
那么假设实时测得梯形溢流坝上游总能量水头H0=2m时,即可求得梯形溢流坝的实时流量17.21m³/s;
或者实时测得梯形溢流坝上游水深Hp=2.8m,则将其可转换成梯形溢流坝上游总能量H0=(Q/B/Hp)2/2g+(Hp-hL)=0.006Q2+1.8带入Q=3×(0.457+0.082sin(π[H0/(1+H0)-0.720]+0.012))×(2gH0 3)1/2,则得到Q=3×(0.457+0.082sin(π((0.006Q2+1.8)/( 0.006Q2+2.8)-0.72)+0.012))×(19.62×(0.006Q2+1.8)3)1/2,等号两侧仅有梯形溢流坝流量Q未知,通过求解可求得梯形溢流坝的实时流量5.71m³/s;完毕。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (10)
1.一种通过计算流体动力学预测梯形溢流坝流量的方法,其特征在于:包括以下步骤:
步骤S1:对缩尺的梯形溢流坝进行物理模型试验,建立梯形溢流坝渠道物理模型;在模型试验流量Qt条件下测得物理模型试验时,沿梯形溢流坝渠道对称纵断面上、梯形溢流坝上游不同过流断面水深hi;
步骤S2:对步骤S1中建立的梯形溢流坝渠道物理模型进行数值计算,计算收敛后导出结果,通过计算流体动力学后处理软件进行后处理;在与步骤S1几何尺寸相同的梯形溢流坝的条件下,计算得到相同模型试验流量Qt条件下测得数值模拟时,沿梯形溢流坝渠道对称纵断面上、梯形溢流坝上游不同过流横断面对应水深Hi;
步骤S3:对步骤S2中数值模拟时对应水深Hi数据集与步骤S1中物理模型试验时对应水深hi数据集进行逐点对比,找到所有数值逐点相对误差绝对值均不超过η%的计算工况,即|Hi-hi|/hi≤η%,输出对应的数值模拟方法;
步骤S4:根据梯形溢流坝的流量Q表达式Q=Cdb(2gH0 3)1/2得到梯形溢流坝流量系数Cd的表达式;其中b为渠道宽度,g为重力加速度,H0为梯形溢流坝上游总能量水头;A、B、C、k和/>为常数系数,/>为相对坝顶宽度,θu表示上游坡度,θd表示下游坡度;
步骤S5:对不同几何条件的梯形溢流坝进行计算流体动力学计算,采用统计产品与服务解决方案SPSS软件对相对坝顶宽度,上游坡度θu,下游坡度θd和梯形溢流坝流量系数Cd进行非线性拟合,得到常数系数A、B、C、k和/>的具体值记为/>和/>,将常数系数具体值反代入梯形溢流坝流量系数Cd的表达式,得到梯形溢流坝流量系数表达式;
步骤S6:将步骤S5得到的梯形溢流坝流量系数表达式代入步骤S4梯形溢流坝的流量Q表达式,得到梯形溢流坝的流量Q的表达式为:
;
将坝顶宽度Lc、坝顶高度hL、上游坡度θu和下游坡度θd安置在梯形溢流坝给定渠道宽度b中,实时测得梯形溢流坝上游总能量水头H0,求得梯形溢流坝中的实时流量;或者测得梯形溢流坝上游水深Hp求解可得梯形溢流坝的流量Q。
2.根据权利要求1所述的一种通过计算流体动力学预测梯形溢流坝流量的方法,其特征在于:步骤S1中对缩尺的梯形溢流坝进行物理模型试验,应保证梯形溢流坝的下游水流状态始终为自由出流,具体方法为:
在梯形溢流坝的上游逐渐增加模型试验流量Qt,在梯形溢流坝下游设置尾门并控制尾门的开度,同时监测梯形溢流坝上游水位ht和下游水位Tt,随着下游水位Tt逐渐上升,始终保证梯形溢流坝上游水位ht变化不超过1mm。
3.根据权利要求2所述的一种通过计算流体动力学预测梯形溢流坝流量的方法,其特征在于:步骤S1中缩尺的梯形溢流坝渠道物理模型,具体参数为:
渠道几何长度为L,梯形溢流坝渠道横跨渠道贯穿布置;
梯形溢流坝几何参数为:坝顶高度hL,坝顶宽度Lc,上游坝顶拐点到上游坝脚拐点的水平距离为Hu,下游坝顶拐点到下游坝脚拐点的水平距离为Hd,其中上游坡度为上游坡度对应的坝顶高度hL与上游坝顶拐点到上游坝脚拐点的水平距离Hu的比值,下游坡度为下游坡度对应的坝顶高度hL与下游坝顶拐点到下游坝脚拐点的水平距离Hd的比值;
上游坡度和下游坡度的陡度均不低于1:1;
物理模型试验采用正态物理模型,综合考虑物理模型水流在阻力平方区要求,选取物理模型线性比尺,物理模型采用重力相似准则设计,角度比/>,流速比/>,流量比/>,糙率比/>,梯形溢流坝渠道流量/>。
4.根据权利要求3所述的一种通过计算流体动力学预测梯形溢流坝流量的方法,其特征在于:步骤S1中物理模型试验时沿梯形溢流坝渠道对称纵断面上、梯形溢流坝上游不同过流断面水深hi,i为随机断面编号;对应的步骤S2中数值模拟时沿梯形溢流坝渠道对称纵断面上、梯形溢流坝上游不同过流横断面对应水深Hi ,i断面编号与步骤S1中相同。
5.根据权利要求4所述的一种通过计算流体动力学预测梯形溢流坝流量的方法,其特征在于:步骤S2中梯形溢流坝进行数值计算,具体为:
建立与步骤S1中梯形溢流坝渠道相同的物理模型,并将物理模型进行网格剖分,输出后缀名为.mesh的计算文件,导入计算流体动力学软件进行数值计算;
在将物理模型进行网格剖分时采用不同尺度网格,得到若干不同的网格方案,若干不同的网格方案在计算流体动力学软件计算中选取不同的湍流模型及不同的数值算法。
6.根据权利要求5所述的一种通过计算流体动力学预测梯形溢流坝流量的方法,其特征在于:步骤S4中:H0的表达式:H0=h+v2/2g;h为梯形溢流坝的坝上水头,v为渠道中平均流速;
渠道中平均流速v的表达式为v=Q/(h+hL)b,h+ hL为梯形溢流坝上游水深Hp;
ξ的表达式为H0/(H0+Lc);上游坡度θu的表达式为θu=atcan(hL/Hu);下游坡度θd的表达式为θd=atcan(hL/Hd),θ单位为rad。
7.根据权利要求6所述的一种通过计算流体动力学预测梯形溢流坝流量的方法,其特征在于:步骤S5中:对不同几何条件的梯形溢流坝进行计算流体动力学计算,具体为:
给定坝顶高度hL的梯形溢流坝,在任意上游坡度θu、下游坡度θd和坝顶宽度Lc条件下,输出梯形溢流坝渠道流量Qp条件下数值模拟时梯形溢流坝上游水深Hp以确定相对坝顶宽度ξ,再根据步骤S4梯形溢流坝的流量Q表达式Q=Cdb(2gH0 3)1/2求得梯形溢流坝流量系数Cd;即每一个上游坡度θu、下游坡度θd、坝顶宽度Lc和梯形溢流坝渠道流量Qp均对应一个相对坝顶宽度ξ和溢流坝流量系数Cd。
8.根据权利要求7所述的一种通过计算流体动力学预测梯形溢流坝流量的方法,其特征在于:步骤S5中:对梯形溢流坝进行计算流体动力学计算,梯形溢流坝的上游坡度θu和下游坡度θd的陡度均不低于1:1,并且梯形溢流坝的下游水流状态始终为自由出流。
9.根据权利要求8所述的一种通过计算流体动力学预测梯形溢流坝流量的方法,其特征在于:步骤S5在梯形溢流坝渠道进行计算流体动力学计算中采用步骤S3中输出的对应的数值模拟方法。
10.根据权利要求9所述的一种通过计算流体动力学预测梯形溢流坝流量的方法,其特征在于:步骤S6中:测得梯形溢流坝上游水深Hp,求解可得梯形溢流坝的流量Q,具体步骤为:
将步骤S5得到的梯形溢流坝流量系数表达式代入步骤S4梯形溢流坝的流量Q表达式,得到梯形溢流坝的流量Q的表达式为:
;
将坝顶宽度Lc、坝顶高度hL、上游坡度θu和下游坡度θd的梯形溢流坝安置在给定渠道宽度b中,测得梯形溢流坝上游水深Hp,转换成梯形溢流坝上游总能量水头H0,即,再将/>中H0全部替换为/>,则/>等号两侧将变为只有梯形溢流坝的流量Q未知的非线性表达式,求解可得梯形溢流坝的流量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311718811.0A CN117408188B (zh) | 2023-12-14 | 2023-12-14 | 一种通过计算流体动力学预测梯形溢流坝流量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311718811.0A CN117408188B (zh) | 2023-12-14 | 2023-12-14 | 一种通过计算流体动力学预测梯形溢流坝流量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117408188A CN117408188A (zh) | 2024-01-16 |
CN117408188B true CN117408188B (zh) | 2024-03-19 |
Family
ID=89491169
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311718811.0A Active CN117408188B (zh) | 2023-12-14 | 2023-12-14 | 一种通过计算流体动力学预测梯形溢流坝流量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117408188B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117829036B (zh) * | 2024-03-05 | 2024-05-14 | 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) | 一种通过cfd确定t型分汊明渠下游支渠流量的方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1281862A1 (en) * | 2001-07-30 | 2003-02-05 | Hans Dr. Grassmann | Run-of-river turbine |
CN106096203A (zh) * | 2016-06-29 | 2016-11-09 | 重庆交通大学 | 一种枢纽大坝的基底施工方法 |
CN107330621A (zh) * | 2017-06-30 | 2017-11-07 | 中国水利水电科学研究院 | 一种基于贝叶斯网络的水库群泄水建筑物失效评估方法 |
CN106844913B (zh) * | 2017-01-09 | 2017-11-24 | 河海大学 | 一种基于三维cfd的滞留气团热力学特性模拟方法 |
CN110904919A (zh) * | 2019-10-30 | 2020-03-24 | 南昌工程学院 | 一种治理小型河道断头浜的装置 |
CN114565211A (zh) * | 2021-12-29 | 2022-05-31 | 郑州大学 | 一种风险传递及叠加作用下梯级水库群溃坝风险后果评估方法 |
CN116163262A (zh) * | 2023-02-18 | 2023-05-26 | 四川大学 | 通过生态补水口结合溢流坝改善分汊河道水流条件的方法 |
CN116384289A (zh) * | 2023-06-05 | 2023-07-04 | 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) | 一种通过计算流体动力学预测墩块式鱼道流量的方法 |
-
2023
- 2023-12-14 CN CN202311718811.0A patent/CN117408188B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1281862A1 (en) * | 2001-07-30 | 2003-02-05 | Hans Dr. Grassmann | Run-of-river turbine |
CN106096203A (zh) * | 2016-06-29 | 2016-11-09 | 重庆交通大学 | 一种枢纽大坝的基底施工方法 |
CN106844913B (zh) * | 2017-01-09 | 2017-11-24 | 河海大学 | 一种基于三维cfd的滞留气团热力学特性模拟方法 |
CN107330621A (zh) * | 2017-06-30 | 2017-11-07 | 中国水利水电科学研究院 | 一种基于贝叶斯网络的水库群泄水建筑物失效评估方法 |
CN110904919A (zh) * | 2019-10-30 | 2020-03-24 | 南昌工程学院 | 一种治理小型河道断头浜的装置 |
CN114565211A (zh) * | 2021-12-29 | 2022-05-31 | 郑州大学 | 一种风险传递及叠加作用下梯级水库群溃坝风险后果评估方法 |
CN116163262A (zh) * | 2023-02-18 | 2023-05-26 | 四川大学 | 通过生态补水口结合溢流坝改善分汊河道水流条件的方法 |
CN116384289A (zh) * | 2023-06-05 | 2023-07-04 | 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) | 一种通过计算流体动力学预测墩块式鱼道流量的方法 |
Non-Patent Citations (3)
Title |
---|
Discharge Estimation for check Dam with Cross-sectional Spillway;Lei Zih-Yi .etal;EGU Genneral Assembly Conference Abstracts;20180430;全文 * |
梯形渠道折线型滚水坝过流能力研究;袁柳;中国优秀硕士学位论文全文数据库 (基础科学辑);20230131(第01期);全文 * |
直立式溢流连拱坝综合流量系数试验研究;刘昱辰 等;人民黄河;20180831(第09期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN117408188A (zh) | 2024-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN117408188B (zh) | 一种通过计算流体动力学预测梯形溢流坝流量的方法 | |
CN114091359A (zh) | 一种水库水沙预测模型训练、水库水沙预测方法及装置 | |
CN103886187B (zh) | 一种基于数据同化的河道水沙实时预测方法 | |
CN112464584A (zh) | 自由表面流的水位和流量推求方法 | |
Zhang et al. | Integrated hydrodynamic model for simulation of river-lake-sluice interactions | |
Nouri et al. | Discharge coefficient in the combined weir-gate structure | |
Dou et al. | Research on calculation method of free flow discharge based on artificial neural network and regression analysis | |
Ran et al. | Three-dimensional numerical simulation of flow in trapezoidal cutthroat flumes based on FLOW-3D | |
CN111311734A (zh) | 铀尾矿库排洪系统水力特性三维数值模拟方法 | |
CN115455867A (zh) | 基于回归分析的坝区流态的推求方法 | |
CN114117965A (zh) | 一种基于两侧半圆柱矩形槽比降效应的灌区渠道测流法 | |
CN116384289B (zh) | 一种通过计算流体动力学预测墩块式鱼道流量的方法 | |
Rezazadeh et al. | Three-dimensional simulation of flow over sharp-crested weirs using volume of fluid method | |
CN117454724B (zh) | 一种通过cfd确定护城河式防洪堤坝几何参数的方法 | |
CN114330156A (zh) | 一种基于中心圆柱矩形槽比尺效应的渠道测流法 | |
Zeinivand et al. | Investigating the process of changing the discharge coefficient of the flow passing through the combined structure of rectangular sharp-crested weir with multiple-gates | |
CN116757111A (zh) | 一种预测带有鱼道坡度的墩块式鱼道流量的方法 | |
CN110847112B (zh) | 一种基于水力学模拟的河道泄洪预警方法 | |
Ghobadian et al. | The effect of downstream curved edge on local scouring at 60 degree open channel junction using SSIIM1 model | |
Liu et al. | Hydraulic model of partial dam break based on sluice gate flow | |
CN117829036B (zh) | 一种通过cfd确定t型分汊明渠下游支渠流量的方法 | |
Kulkarni et al. | An energy perspective of composite broad crested weir for measuring accurate discharge | |
Al-Suhili et al. | Prediction of the discharge coefficient for a Cipolletti weir with rectangular bottom opening | |
Khosravi et al. | Prediction of side thermal buoyant discharge in the cross flow using multi-objective evolutionary polynomial regression (EPR-MOGA) | |
Yin et al. | Walladapting local eddy-viscosity turbulence model for TELEMAC3D |
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 |