王牌gydF4y2Ba 土木工程的发展gydF4y2Ba 1687 - 8094gydF4y2Ba 1687 - 8086gydF4y2Ba HindawigydF4y2Ba 10.1155 / 2021/5529792gydF4y2Ba 5529792gydF4y2Ba 研究文章gydF4y2Ba 水动力性能和空化分析在大坝底部网点使用CFD模型gydF4y2Ba https://orcid.org/0000 - 0002 - 5489 - 8931gydF4y2Ba Aminoroayaie YaminigydF4y2Ba OmidgydF4y2Ba https://orcid.org/0000 - 0001 - 8789 - 0285gydF4y2Ba 穆萨维gydF4y2Ba 美国HoomangydF4y2Ba https://orcid.org/0000 - 0003 - 3546 - 6663gydF4y2Ba KavianpourgydF4y2Ba m·R。gydF4y2Ba Safari GhalehgydF4y2Ba 拉明gydF4y2Ba 黄gydF4y2Ba 法明gydF4y2Ba 土木工程学系gydF4y2Ba k . n . Toosi科技大学gydF4y2Ba 德黑兰gydF4y2Ba 伊朗gydF4y2Ba kntu.ac.irgydF4y2Ba 2021年gydF4y2Ba 18gydF4y2Ba 5gydF4y2Ba 2021年gydF4y2Ba 2021年gydF4y2Ba 11gydF4y2Ba 1gydF4y2Ba 2021年gydF4y2Ba 25gydF4y2Ba 4gydF4y2Ba 2021年gydF4y2Ba 6gydF4y2Ba 5gydF4y2Ba 2021年gydF4y2Ba 18gydF4y2Ba 5gydF4y2Ba 2021年gydF4y2Ba 2021年gydF4y2Ba 版权©2021 Omid Aminoroayaie Yamini et al。gydF4y2Ba 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。gydF4y2Ba

底网点是大坝的重要结构,负责控制流量,操作或删除水库沉积。服务门控制出口流量,每当这个门坏了,位于上游的应急门是利用。空化现象是一种常见的底店的问题由于快速流动转移。目前的研究是一种流动的数值研究模式在大坝的底部出口不同的门开口Flow-3D软件的使用和提高gydF4y2Ba k-εgydF4y2Ba湍流模型。Sardab大坝进行调查,一个土坝在伊斯法罕(伊朗)。最大速度的100%开的门,豪厄尔他们阀18 m / s在门下面的部分,和40%的最大速度等于闸门开度23.1 m / s。50%的服务和应急闸门阀的上游地区,减少所需的压力值。此外,在两者之间的区域紧急和服务盖茨,压力值降低。气蚀的可能性在这个领域可以减少安装限。流模式Sardab大坝的底部出口相对稳定和适当的条件,和没有麻烦的水力现象,如当地的漩涡,不良的压力变化,在隧道和速度,没有流动分离流进入分支的关键领域。gydF4y2Ba

 1。介绍gydF4y2Ba

底部利用媒体作为一个大坝的水工建筑物控制水库蓄水,水库在紧急情况下疏散,进入水库的沉积物。因此,他们需要仔细设计和有害因素识别(gydF4y2Ba 1gydF4y2Ba]。另一方面,为了确保通道的正确操作,相关的液压流体力学的设施,包括盖茨和阀门,放电容量应仔细检查和执行水力模型测试(gydF4y2Ba 2gydF4y2Ba]。通过研究失败导致大坝的终端盖茨已经决定,空化现象和闸门振动的主要原因是赔偿,因为气泡崩溃的地方,压力波动强度提出了(gydF4y2Ba 3gydF4y2Ba- - - - - -gydF4y2Ba 8gydF4y2Ba]。在底部,两阶段高速空气流传播。由于部门流和即兴的转换从加压流到自由表面流动,压力急剧损失值下游的门。空化是一种破坏性的现象发生,因为这样的负面压力。曝气增氧机管通常是采用以控制这一现象(gydF4y2Ba 9gydF4y2Ba- - - - - -gydF4y2Ba 12gydF4y2Ba]。gydF4y2Ba

萨达特Helbar et al。(2021)研究了底部出口门的大小和形状,从而影响流出分泌物,冲洗效率。他们的研究的目的是调查该地区的影响,形状,盖茨和底部出口数量的速度和浓度沉积物(gydF4y2Ba 13gydF4y2Ba]。gydF4y2Ba

在这些结构中,由于高速流动和渠道的水位波动,流可能独立于信道的墙,并在本地流压力可能会减少。如果流动压力小于水的蒸汽压,水从液体到气体状态变化,和气泡的形式。可能携带气泡流高压力区为了崩溃,和一个负压力波进入流。如果气泡壁表面附近的爆炸,他们会损害该频道的墙gydF4y2Ba 14gydF4y2Ba- - - - - -gydF4y2Ba 16gydF4y2Ba]。杨et al。(2020)在一个三维CFD建模显示空气流动行为。空气流例之间的需求存在着很大的差别。不是所有空缺的同步放电结果在空气最大的需求gydF4y2Ba 17gydF4y2Ba]。gydF4y2Ba

服务门的下游地区,该地区应急门和服务之间的门,以及门槽,创建一个粗糙表面对空化流持有风险最高的(gydF4y2Ba 18gydF4y2Ba,gydF4y2Ba 19gydF4y2Ba]。空化指数,当地的压力和流速的函数,作为一个重要的参数来评估潜在的空泡(gydF4y2Ba 20.gydF4y2Ba]。Daneshmand et al。(2007)实验研究了流的水力特性在不同的服务门开口Sivand大坝的底部出口。他们的研究显示流湍流量的增加应急门缝隙与服务从85%以上,闸门开启高度和空化指数在所有开口大于临界极限(gydF4y2Ba 21gydF4y2Ba]。两个空化指数控制方法可以控制气蚀损伤通过修改几何结构和流曝气(gydF4y2Ba 22gydF4y2Ba]。Kavianpour(1997)进行实验,以确定曝气对斜坡下游压力波动的影响。他们得出的结论是,进气减少严重的压力波动和提高它的意思是,这反过来又降低了空化的可能性。他们还研究了流曝气对压力场的影响的结构,包括波动的能量谱,偏态,波动的清晰度,和负压波动的风险减少。本研究表明,充气压力场的结构变化(gydF4y2Ba 23gydF4y2Ba]。Khazaei et al。(2015)研究了空化现象在地面加速度峰值能力Lorestan大坝的底部进气数值和在实验室里。他们表明,改变通道的坡度10至12%,改变通道1到3.5米的长度,并减少通道墙壁的斜率1到2度可以是一个不错的选择,消除气蚀的风险(gydF4y2Ba 24gydF4y2Ba]。gydF4y2Ba

阮et al。(2007)研究了中控室内的水力性能使用基础构皮滩大坝的底部放电通道在实验室。结果表明增加进气流量的流,以防减少曝气后通道的坡度。这个校正的有效作用,防止放电通道气蚀损伤(gydF4y2Ba 25gydF4y2Ba]。吴et al。(2007)研究了水力参数如流量压力,水电导系数和下游流量剖面在不同服务门开口和多样化的上游渠道负责人龙滩大坝的底部出口在实验室。结果表明,底流的压力值在通道闸下游崛起服务闸门开启高度的增加,和水电导系数在小型和大型高于中间开口(gydF4y2Ba 26gydF4y2Ba]。李et al。(2020)调查控制导叶开度对压力分布的影响和空化数值和实验室模型的指数。等基本假设2的45度角和边缘值的长度宽度比门(通过保持门高度不变),研究10日,30日,50岁和100%的控制门不同的2 - 3 -,四头。结果是空化指数高10%和100%控制门的机会,所以在30%,50%,和70%的门开口,空化的概率更高(gydF4y2Ba 27gydF4y2Ba]。gydF4y2Ba

Nikseresht et al。(2012)模拟底部摄入5号Sefidrud大坝使用一个三维有限体积法和表明,隧道内的最低压力发生在20%的应急门,在这种情况下,气蚀指数低于临界值。曝气的使用提出了消除空化的风险。它正试图独立的高速流沿着隧道通过使用适当的曝气系统。然后,流动可以实现在隧道衬砌空洞的风险是腐烂的gydF4y2Ba 28gydF4y2Ba]。在这个项目中,渠道的水力性能和几何校正评估正常运行的深层去除大门开在不同的州。自深移除门流分为两个部分,与农业灌溉控制器与滑动阀门阀和控制器通过一个分支,研究不同功能状态被认为是。单阀操作,单服务门函数,豪厄尔他们阀和服务门的同时操作这个项目的目标。接下来,该频道的液压操作在不同的情况下同意数值模型将完成,和结论。gydF4y2Ba

1.1。Flow-3D模型和控制方程gydF4y2Ba

数值模拟流场的网点和水工建筑物底部一般来说,溢洪道和能量耗散系统设计这种结构是至关重要的。研究人员在这种模拟的经验已经表明,Flow-3D软件有更好的能力在建模中这种类型的液压结构现有的软件包。Flow-3D工程师提供了一个完整的和通用的CFD仿真平台研究液体和气体的动态行为在一个广泛的工业应用和物理过程。Flow-3D的重要特性之一,液压分析自由表面流动的能力模型,建模的使用受到(体积的液体)技术报告由赫特和尼科尔斯(1981)(gydF4y2Ba 29日gydF4y2Ba]。受到平流方案,数值方法,允许程序员跟踪接口的形状和位置,但它不是一个独立的流求解算法。另一方面,可以考虑流体运动方程,其中包括控制卷假设不可压缩流体的连续性方程和恒流(方程(gydF4y2Ba 1gydF4y2Ba))和控制体积中的动量方程考虑湍流(方程(gydF4y2Ba 2gydF4y2Ba))(gydF4y2Ba 30.gydF4y2Ba]。gydF4y2Ba (1)gydF4y2Ba ∂gydF4y2Ba UgydF4y2Ba 我gydF4y2Ba ∂gydF4y2Ba xgydF4y2Ba 我gydF4y2Ba =gydF4y2Ba 0gydF4y2Ba ,gydF4y2Ba (2)gydF4y2Ba ∂gydF4y2Ba UgydF4y2Ba 我gydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba +gydF4y2Ba UgydF4y2Ba jgydF4y2Ba ∂gydF4y2Ba UgydF4y2Ba 我gydF4y2Ba ∂gydF4y2Ba xgydF4y2Ba jgydF4y2Ba =gydF4y2Ba −gydF4y2Ba 1gydF4y2Ba ρgydF4y2Ba ∂gydF4y2Ba PgydF4y2Ba ∂gydF4y2Ba xgydF4y2Ba 我gydF4y2Ba +gydF4y2Ba ∂gydF4y2Ba ∂gydF4y2Ba xgydF4y2Ba 我gydF4y2Ba vgydF4y2Ba ∂gydF4y2Ba UgydF4y2Ba 我gydF4y2Ba ∂gydF4y2Ba xgydF4y2Ba jgydF4y2Ba −gydF4y2Ba ugydF4y2Ba 我gydF4y2Ba ′gydF4y2Ba ugydF4y2Ba jgydF4y2Ba ′gydF4y2Ba ¯gydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba ρgydF4y2Ba 是流体的密度;gydF4y2Ba xgydF4y2Ba 我gydF4y2Ba 和gydF4y2Ba xgydF4y2Ba jgydF4y2Ba 笛卡尔坐标;gydF4y2Ba UgydF4y2Ba 我gydF4y2Ba 和gydF4y2Ba UgydF4y2Ba jgydF4y2Ba 笛卡尔速度矢量的组件吗gydF4y2Ba vgydF4y2Ba ;gydF4y2Ba PgydF4y2Ba 是压力;和gydF4y2Ba ugydF4y2Ba 我gydF4y2Ba ′gydF4y2Ba ugydF4y2Ba jgydF4y2Ba ′gydF4y2Ba ¯gydF4y2Ba 意味着雷诺应力张量。gydF4y2Ba

two-equation模型重整化群(RNG)gydF4y2Ba k-εgydF4y2Ba用于确定的雷诺应力张量为湍流动量方程。这个模型的方程如下:gydF4y2Ba

湍流动能方程gydF4y2Ba KgydF4y2Ba:gydF4y2Ba (3)gydF4y2Ba ∂gydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba ρgydF4y2Ba kgydF4y2Ba +gydF4y2Ba ∂gydF4y2Ba ∂gydF4y2Ba xgydF4y2Ba 我gydF4y2Ba ρgydF4y2Ba kgydF4y2Ba UgydF4y2Ba 我gydF4y2Ba =gydF4y2Ba ∂gydF4y2Ba ∂gydF4y2Ba xgydF4y2Ba jgydF4y2Ba αgydF4y2Ba kgydF4y2Ba μgydF4y2Ba egydF4y2Ba fgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba kgydF4y2Ba ∂gydF4y2Ba xgydF4y2Ba jgydF4y2Ba +gydF4y2Ba GgydF4y2Ba kgydF4y2Ba −gydF4y2Ba ρgydF4y2Ba εgydF4y2Ba 。gydF4y2Ba

湍流动能消耗速率方程gydF4y2Ba εgydF4y2Ba :gydF4y2Ba (4)gydF4y2Ba ∂gydF4y2Ba ∂gydF4y2Ba tgydF4y2Ba ρgydF4y2Ba εgydF4y2Ba +gydF4y2Ba ∂gydF4y2Ba ∂gydF4y2Ba xgydF4y2Ba 我gydF4y2Ba ρgydF4y2Ba εgydF4y2Ba UgydF4y2Ba 我gydF4y2Ba =gydF4y2Ba ∂gydF4y2Ba ∂gydF4y2Ba xgydF4y2Ba jgydF4y2Ba αgydF4y2Ba εgydF4y2Ba μgydF4y2Ba egydF4y2Ba fgydF4y2Ba fgydF4y2Ba ∂gydF4y2Ba εgydF4y2Ba ∂gydF4y2Ba xgydF4y2Ba jgydF4y2Ba +gydF4y2Ba CgydF4y2Ba 1gydF4y2Ba εgydF4y2Ba εgydF4y2Ba kgydF4y2Ba GgydF4y2Ba kgydF4y2Ba −gydF4y2Ba CgydF4y2Ba 2gydF4y2Ba εgydF4y2Ba ρgydF4y2Ba εgydF4y2Ba 2gydF4y2Ba kgydF4y2Ba −gydF4y2Ba RgydF4y2Ba εgydF4y2Ba 。gydF4y2Ba

在上面的方程中,gydF4y2Ba αgydF4y2Ba kgydF4y2Ba和gydF4y2Ba αgydF4y2Ba εgydF4y2Ba有效的普朗特数倒数吗gydF4y2Ba kgydF4y2Ba 和gydF4y2Ba εgydF4y2Ba ,分别。gydF4y2Ba CgydF4y2Ba 1gydF4y2Ba εgydF4y2Ba 和gydF4y2Ba CgydF4y2Ba 2gydF4y2Ba εgydF4y2Ba 是常数,值分别为1.42和1.68;gydF4y2Ba μgydF4y2Ba effgydF4y2Ba 是有效的粘度。RNG之间的主要区别gydF4y2Ba kgydF4y2Ba- - - - - -gydF4y2Ba εgydF4y2Ba模型和标准gydF4y2Ba kgydF4y2Ba- - - - - -gydF4y2Ba εgydF4y2Ba模型是RNG模型有一个额外的术语,gydF4y2Ba RgydF4y2Ba εgydF4y2Ba 的准确性,极大地提高了快速紧张的流动。gydF4y2Ba

 2。材料和方法gydF4y2Ba 2.1。案例研究Sardab大坝与豪厄尔他们阀的出口渠道gydF4y2Ba

Sardab大坝在伊朗建立在Sardab河。坝址的高度是海拔2712米的大坝的水库面积270公顷,总容量4800万立方米(表gydF4y2Ba 1gydF4y2Ba)。gydF4y2Ba

Sardab大坝的一般规范。gydF4y2Ba

大坝gydF4y2Ba 类型gydF4y2Ba 波峰高度(m·a·s·l。gydF4y2Ba ∗gydF4y2Ba )gydF4y2Ba 大坝高度(米)gydF4y2Ba 波峰长度(米)gydF4y2Ba
土坝gydF4y2Ba 2528年gydF4y2Ba 53gydF4y2Ba 752年gydF4y2Ba
储层gydF4y2Ba 最大水位(及)(m)gydF4y2Ba 正常水位(m)gydF4y2Ba 最低水位(米)gydF4y2Ba 最小体积(mgydF4y2Ba3gydF4y2Ba)gydF4y2Ba
2527年gydF4y2Ba 2524.6gydF4y2Ba 2488.6gydF4y2Ba 1000000年gydF4y2Ba

∗gydF4y2Ba 米海拔。gydF4y2Ba

两个滑动门(服务和紧急)尺寸(高度∗宽度)1.4×1.1 m位于一行在大坝的底部通道。几何细节Sardab大坝的底部出口如图gydF4y2Ba 1gydF4y2Ba。gydF4y2Ba

几何细节Sardab大坝的底部出口(所有尺寸在毫米)。gydF4y2Ba

 2.2。流量通过通道闸门操作期间gydF4y2Ba

为了计算通道流量,假设盖茨是完全开放的,与伯努利方程和考虑到渠道的帮助下损失,gydF4y2Ba (5)gydF4y2Ba 问gydF4y2Ba =gydF4y2Ba VgydF4y2Ba dgydF4y2Ba 一个gydF4y2Ba dgydF4y2Ba =gydF4y2Ba 一个gydF4y2Ba dgydF4y2Ba 2gydF4y2Ba ggydF4y2Ba HgydF4y2Ba ∑gydF4y2Ba ξgydF4y2Ba +gydF4y2Ba ξgydF4y2Ba 出gydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba 问gydF4y2Ba (mgydF4y2Ba3gydF4y2Ba/ s)是体积流率;gydF4y2Ba VgydF4y2Ba dgydF4y2Ba (m / s)的速度参考面积;gydF4y2Ba 一个gydF4y2Ba dgydF4y2Ba (mgydF4y2Ba2gydF4y2Ba)是参考截面积;gydF4y2Ba ggydF4y2Ba (米/秒gydF4y2Ba2gydF4y2Ba)是重力加速度;gydF4y2Ba HgydF4y2Ba (m)是;gydF4y2Ba ξgydF4y2Ba 损耗系数;和gydF4y2Ba ξgydF4y2Ba 出gydF4y2Ba 在流出部分损失系数。gydF4y2Ba

在通道的不同部分,计算总水头损失和总的流量,应该选一部分作为参考,大门下的位置被认为是部分的引用。gydF4y2Ba

水头损失的势能转化为动能。头的摩擦阻力损失是由于底部出口系统(阀门、盖茨、配件、管道、入口,出口损失,等等)。我们计算总水头损失,gydF4y2Ba (6)gydF4y2Ba ∑gydF4y2Ba ξgydF4y2Ba VgydF4y2Ba 我gydF4y2Ba 2gydF4y2Ba 2gydF4y2Ba ggydF4y2Ba =gydF4y2Ba ∑gydF4y2Ba hgydF4y2Ba fgydF4y2Ba 我gydF4y2Ba =gydF4y2Ba 1.4224gydF4y2Ba VgydF4y2Ba ogydF4y2Ba 2gydF4y2Ba 2gydF4y2Ba ggydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba VgydF4y2Ba 我gydF4y2Ba (m / s)的速度在每个部分;gydF4y2Ba VgydF4y2Ba ogydF4y2Ba (m / s)的速度参考面积;和gydF4y2Ba hgydF4y2Ba fgydF4y2Ba 我gydF4y2Ba (m)是水头损失。gydF4y2Ba

的收缩系数由于收缩45度角约为0.95。gydF4y2Ba (7)gydF4y2Ba ξgydF4y2Ba 出gydF4y2Ba =gydF4y2Ba 1.05gydF4y2Ba 一个gydF4y2Ba dgydF4y2Ba εgydF4y2Ba 一个gydF4y2Ba egydF4y2Ba 2gydF4y2Ba 。gydF4y2Ba

在上面的方程中,gydF4y2Ba 一个gydF4y2Ba egydF4y2Ba (mgydF4y2Ba2gydF4y2Ba)是输出横截面积,gydF4y2Ba 一个gydF4y2Ba dgydF4y2Ba (mgydF4y2Ba2gydF4y2Ba)是横断面参考地区,gydF4y2Ba εgydF4y2Ba的收缩系数是截面压缩功能。gydF4y2Ba

根据所述理论问题和方程(gydF4y2Ba 5gydF4y2Ba)- (gydF4y2Ba 7gydF4y2Ba),最大体积流率通过信道估计26米gydF4y2Ba3gydF4y2Ba/ s(参考地区速度是18.9米/秒)。表gydF4y2Ba 2gydF4y2Ba和图gydF4y2Ba 2gydF4y2Ba显示流量在不同开口和头上。gydF4y2Ba

通道的流量通过服务盖茨和不同的正面的开口。gydF4y2Ba

一个gydF4y2Ba/gydF4y2Ba 一个gydF4y2BaogydF4y2Ba HgydF4y2Ba马克斯gydF4y2Ba= 51.35米gydF4y2Ba HgydF4y2Ba= 48.95米gydF4y2Ba HgydF4y2Ba= 45米gydF4y2Ba HgydF4y2Ba= 40米gydF4y2Ba HgydF4y2Ba= 35米gydF4y2Ba HgydF4y2Ba= 30米gydF4y2Ba HgydF4y2Ba= 20米gydF4y2Ba
0.1gydF4y2Ba 3.17gydF4y2Ba 3.1gydF4y2Ba 2.97gydF4y2Ba 2.8gydF4y2Ba 2.62gydF4y2Ba 2.43gydF4y2Ba 1.98gydF4y2Ba
0.2gydF4y2Ba 6.29gydF4y2Ba 6.15gydF4y2Ba 5.89gydF4y2Ba 5.56gydF4y2Ba 5.2gydF4y2Ba 4.81gydF4y2Ba 3.93gydF4y2Ba
0.3gydF4y2Ba 9.29gydF4y2Ba 9.07gydF4y2Ba 8.7gydF4y2Ba 8.2gydF4y2Ba 7.67gydF4y2Ba 7.1gydF4y2Ba 5.8gydF4y2Ba
0.4gydF4y2Ba 12.18gydF4y2Ba 11.9gydF4y2Ba 11.41gydF4y2Ba 10.75gydF4y2Ba 10.06gydF4y2Ba 9.31gydF4y2Ba 7.6gydF4y2Ba
0.5gydF4y2Ba 15.03gydF4y2Ba 14.67gydF4y2Ba 14.07gydF4y2Ba 13.26gydF4y2Ba 12.41gydF4y2Ba 11.49gydF4y2Ba 9.38gydF4y2Ba
0.6gydF4y2Ba 17.67gydF4y2Ba 17.25gydF4y2Ba 16.54gydF4y2Ba 15.59gydF4y2Ba 14.59gydF4y2Ba 13.5gydF4y2Ba 11.03gydF4y2Ba
0.7gydF4y2Ba 20.06gydF4y2Ba 19.59gydF4y2Ba 18.78gydF4y2Ba 17.71gydF4y2Ba 16.56gydF4y2Ba 15.33gydF4y2Ba 12.52gydF4y2Ba
0.8gydF4y2Ba 22.29gydF4y2Ba 21.77gydF4y2Ba 20.87gydF4y2Ba 19.68gydF4y2Ba 18.41gydF4y2Ba 17.04gydF4y2Ba 13.91gydF4y2Ba
0.9gydF4y2Ba 24.35gydF4y2Ba 23.77gydF4y2Ba 22.8gydF4y2Ba 21.49gydF4y2Ba 20.1gydF4y2Ba 18.61gydF4y2Ba 15.2gydF4y2Ba
1gydF4y2Ba 26.64gydF4y2Ba 26gydF4y2Ba 24.93gydF4y2Ba 23.51gydF4y2Ba 21.99gydF4y2Ba 20.36gydF4y2Ba 16.62gydF4y2Ba

体积流率在不同开口和头上。gydF4y2Ba

 2.3。数值模型gydF4y2Ba

Flow-3D的第一个特点是它雇佣了一个高度可变矩形网格的网格系统。这个特点使得网格或几何互相分离。用更简单的术语来说,它不使用一个固定的磨削系统连接到几何或有限元素。也可以使用多个网格网格系统来提高效率和灵活性。整个媒体流刚体在三维空间中所有的细节设计通过SolidWorks软件模型底部渠道的水力条件在Flow-3D软件(数据流gydF4y2Ba 3gydF4y2Ba和gydF4y2Ba 4gydF4y2Ba)。应该注意的是,五倍隧道直径(7米)添加到隧道的长度和上游来创建一个充分发展流动。gydF4y2Ba

完成底部出口几何与豪厄尔他们阀、增氧机、蝶阀在SolidWorks软件。gydF4y2Ba

的细节建设3 d几何Sardab大坝的底部出口在m(所有尺寸)。gydF4y2Ba

为了覆盖整个刚体与分支阀出口,网应该选择这中尽可能少的空间通道区域空的空间。下图显示了网格的详细信息(12825258细胞),计算细胞,底部出口边界条件。应该注意的是,Flow-3D数值模型的元素是长方体的元素类型。细胞大小网格块1,2,3,4在每个方向(gydF4y2Ba XgydF4y2Ba,gydF4y2Ba Y,gydF4y2Ba和gydF4y2Ba ZgydF4y2Ba)是2.3,2.3,3.7,和2.1厘米。的问题之一,严重影响增加计算精度是确定适当的边界条件。有六个不同的方面来定义边界条件在笛卡尔坐标,根据轴的正方向包括gydF4y2Ba XgydF4y2Ba最小值gydF4y2Ba,gydF4y2Ba XgydF4y2Ba马克斯gydF4y2Ba,gydF4y2Ba YgydF4y2Ba最小值gydF4y2Ba,gydF4y2Ba YgydF4y2Ba马克斯gydF4y2Ba,gydF4y2Ba ZgydF4y2Ba最小值gydF4y2Ba,gydF4y2Ba ZgydF4y2Ba马克斯gydF4y2Ba。值得注意的是,所有这些规范被定义在一个街区,和独立的边界条件必须为每个定义的几个街区(数字gydF4y2Ba 5gydF4y2Ba和gydF4y2Ba 6gydF4y2Ba)。边界条件应用到各个方面的机制后,分配每个州的一个方面,所有细胞输入方面的输入上面的方程和边界条件。这是固定在不同的时间解决方程和运行软件。压力可以被定义为的压力或停滞的压力。液体的量也可以指定为一个百分比的液体分数或流体软件的高度。流体速度在所有三个方向都可以进入软件。软件也可以接受上述每个参数时间序列(gydF4y2Ba 31日gydF4y2Ba- - - - - -gydF4y2Ba 33gydF4y2Ba]。gydF4y2Ba

数值模型的网格和啮合块底部网点和分支Flow-3D数值模型。gydF4y2Ba

底部边界条件应用于流仿真。gydF4y2Ba

在目前的研究中,一个网格敏感性分析进行了5种不同的细胞大小。表gydF4y2Ba 3gydF4y2Ba显示网格敏感性分析的结果在参考区域速度100%闸门和阀门开度。可以看到,如果细胞大小小于主要网格,它不影响仿真的准确性,如果单元格大小大于主网格,错误会增加两倍。gydF4y2Ba

网格敏感性分析。gydF4y2Ba

参数gydF4y2Ba 类型的网gydF4y2Ba
单元尺寸小20%gydF4y2Ba 单元尺寸小10%gydF4y2Ba 主要gydF4y2Ba 单元尺寸大10%gydF4y2Ba 单元尺寸大20%gydF4y2Ba
速度(米/秒)gydF4y2Ba 18.21gydF4y2Ba 18.13gydF4y2Ba 18.05gydF4y2Ba 17.42gydF4y2Ba 17.1gydF4y2Ba
分析速度(米/秒)gydF4y2Ba 18.9gydF4y2Ba 18.9gydF4y2Ba 18.9gydF4y2Ba 18.9gydF4y2Ba 18.9gydF4y2Ba
误差(%)gydF4y2Ba 3.65gydF4y2Ba 4.07gydF4y2Ba 4.49gydF4y2Ba 7.83gydF4y2Ba 9.52gydF4y2Ba
3所示。结果与讨论gydF4y2Ba 3.1。数值模拟结果与评价Flow-3D软件gydF4y2Ba

为了评估数值模型结果的校准和验证在不同的底部出口门开口,手工分析的输出流量应与输出流量数值模型的相同应用头条件使用解析解在前面几节(图执行gydF4y2Ba 7gydF4y2Ba)。提到的破坏现象,参数如流量、速度、水的流体压力、气流速度进入流,和水的深度是有效的。通过测量这些参数,就可以阻止发生故障前(gydF4y2Ba 34gydF4y2Ba]。此外,由于控制方程的性质,有必要启动流分析和固定边界条件,以便最终达到稳态流模型中在一个合适的时间。输出流量值稳定流动状态的数值模型与解析解的结果图gydF4y2Ba 8gydF4y2Ba。值得注意的是数值模型的头被认为是基于标准的48.95米,和分支阀的开口被认为是100%gydF4y2Ba 35gydF4y2Ba]。gydF4y2Ba

数值模型结果输出流量在不同期间门开口Flow-3D数值模型实现。gydF4y2Ba

解析解的比较结果和数值模拟的输出流量在不同阀开口。gydF4y2Ba

显示在图gydF4y2Ba 8gydF4y2Ba平均相对误差计算,数值模型的输出流量在不同开口Sardab大坝的底部出口门大约是4%,这是一个可以接受的和娇小的数量。因此,根据这些条件,Flow-3D校准和验证数值模拟,以及其他所需的水力参数,包括压力值,速度,和液压领域流底部出口敏感的点,可以调查。数据gydF4y2Ba 9gydF4y2Ba和gydF4y2Ba 10gydF4y2Ba显示Sardab大坝的底部出口的输出流量值在不同的机会与解析解的结果。gydF4y2Ba

结果Flow-3D阀的输出流量数值模型在服务门的不同开口模式执行。gydF4y2Ba

解析解的比较结果和数值模拟的输出流量在不同门开口。gydF4y2Ba

 3.2。数值模拟结果5门模式和分支阀门开度100%gydF4y2Ba

结果显示在表gydF4y2Ba 4gydF4y2Ba空缺包括20、40、60、80和100%在一个标准的48.95米的高度。表gydF4y2Ba 4gydF4y2Ba提出更改的值在底部出口输出流量不同的门开口。gydF4y2Ba

底部出口输出流量的变化对不同门开口。gydF4y2Ba

门开口(%)gydF4y2Ba 流量(门)(mgydF4y2Ba3gydF4y2Ba/秒)gydF4y2Ba 流量(阀)(mgydF4y2Ba3gydF4y2Ba/秒)gydF4y2Ba
0gydF4y2Ba 0.00gydF4y2Ba 0.00gydF4y2Ba
20.gydF4y2Ba 6.32gydF4y2Ba 5.30gydF4y2Ba
40gydF4y2Ba 12.90gydF4y2Ba 5.26gydF4y2Ba
60gydF4y2Ba 17.51gydF4y2Ba 5.14gydF4y2Ba
80年gydF4y2Ba 22.30gydF4y2Ba 4.63gydF4y2Ba
One hundred.gydF4y2Ba 24.70gydF4y2Ba 3.86gydF4y2Ba
3.3。速度的变化不同的Gate条件gydF4y2Ba

液压参数建模结果表明,速度值在100%在门下面的部分是关于18米/秒(图gydF4y2Ba 11gydF4y2Ba)。门开40%下的最大速度是23.1米/秒(图gydF4y2Ba 12gydF4y2Ba)。速度值为80%,60%,20%可以在数据中找到gydF4y2Ba 13gydF4y2Ba- - - - - -gydF4y2Ba 15gydF4y2Ba,分别。速度值检查在敏感的点,包括门槽,在输入流中地区分支阀门,显示不良液压不发生变化,严重的变化速度增加不观察到在这些领域gydF4y2Ba 36gydF4y2Ba]。gydF4y2Ba

闸门和阀门的输出速度100%闸门和阀门开度。gydF4y2Ba

闸门和阀门的输出速度40%闸门和阀门开度100%。gydF4y2Ba

闸门和阀门的输出速度80%闸门和阀门开度100%。gydF4y2Ba

闸门和阀门的输出速度60%闸门和阀门开度100%。gydF4y2Ba

闸门和阀门的输出速度20%闸门和阀门开度100%。gydF4y2Ba

门槽大约2 - 3 m / s的速度值,在这种情况下没有观察到流分离和不良的循环流槽。这些条件函数几乎相同的开口。图gydF4y2Ba 16gydF4y2Ba显示变化的规模和速度矢量在门槽导叶开度100%。gydF4y2Ba

速度变化的规模和向量在门槽100%闸门和阀门开度。gydF4y2Ba

调查速度值的分支和加劲肋区呈现在图gydF4y2Ba 17gydF4y2Ba。没有在这个领域流动分离等不利条件下,分支入口处和速度值是可变的,2 - 5 m / s。速度值增加通过分支入口。gydF4y2Ba

速度变化的大小和向量分支区域100%的闸门和阀门开度。gydF4y2Ba

 3.4。压力参数的变化在不同的门开口gydF4y2Ba

在出口通道,流体运动的压力可能会增加流向的一个障碍。这种变化称为逆压梯度。流体的流动边界层区域受到越来越大的压力,这种液体速度也会减缓。然而,由于边界层内流体的动能较低,它可能会停滞不前,被逆转,导致边界层分离和偏离。主要的分离流的边界被称为分离现象,这是由逆压力梯度引起的。一个逆压力梯度是一个必要条件而不是充分条件的分离流动。换句话说,可以有一个逆压力梯度没有分离,虽然没有逆压力梯度分离不能发生。压力多通道区域在不同的门开口中可以看到数据gydF4y2Ba 18gydF4y2Ba- - - - - -gydF4y2Ba 21gydF4y2Ba。gydF4y2Ba

压力在整个模型的三维视图。gydF4y2Ba

量的压力导致门和通道地区100%的分支在打开阀门开度(100%)。gydF4y2Ba

量的压力导致门和通道地区60%的分支在打开阀门开度(100%)。gydF4y2Ba

量的压力导致门和通道地区20%的分支在打开阀门开度(100%)。gydF4y2Ba

流下游闸门下创建一个循环流动,这是一个尖锐的压力损失的主要特征。这个压力损失是一个函数的门打开,门背后的水头,通道的几何形状。另一方面,严重的压力波动减少区域的局部压力,以及潜在的空泡增加由于高速(gydF4y2Ba 37gydF4y2Ba]。我们可以看到数据gydF4y2Ba 22gydF4y2Ba- - - - - -gydF4y2Ba 25gydF4y2Ba,沿着隧道和门地区,闸门槽和分支区域,减少不良的严重压力并没有发生,而且没有不良的水力现象如流动分离和当地的漩涡。gydF4y2Ba

压力分布在大门开放模式(阀门开度100%)100%。gydF4y2Ba

压力分布在大门开放模式(阀门开度100%)20%。gydF4y2Ba

流型加劲肋和不流动分离在这个领域(阀门开度100%)。gydF4y2Ba

流型在虹吸和流动分离在这个领域(阀门开度100%)。gydF4y2Ba

 3.5。数值模拟结果为50%紧急门和服务门打开gydF4y2Ba 3.5.1。流体体积分数值每50%的服务和紧急打开门gydF4y2Ba

根据图gydF4y2Ba 26gydF4y2Ba的曝气两个盖茨通过服务提供之间的门。gydF4y2Ba

流体体积分数等值两个盖茨在50%之间开放状态。gydF4y2Ba

 3.5.2。压力变化值每50%的服务和紧急打开门gydF4y2Ba

我们可以看到数据gydF4y2Ba 27gydF4y2Ba,在门口的上游地区,所需的压力值下降,在两个门之间的区域,压力值降低。减少气蚀的可能性在这个领域与中控室内的安装。gydF4y2Ba

压力变化在底部出口的隧道Sardab大坝的时候50%的服务和应急门。gydF4y2Ba

 3.5.3。数值模拟结果完全关闭闸门阀打开分支模式和100%gydF4y2Ba

流的水力条件检查完全封闭阀打开状态和100%分支阀使用数值模拟。gydF4y2Ba

 3.5.4。在底部出口压力值的隧道Sardab大坝和分支gydF4y2Ba

我们可以看到数据gydF4y2Ba 28gydF4y2Ba和gydF4y2Ba 29日gydF4y2Ba阀门完全关闭时,没有循环流动,不受欢迎的压力变化的隧道。gydF4y2Ba

压力变化在底部出口的隧道Sardab大坝开分支阀完全关闭门和100%。gydF4y2Ba

压力变化在底部出口Sardab大坝的分支领域的开放分支阀完全关闭门和100%。gydF4y2Ba

 3.5.5。速度值在底部出口的隧道Sardab大坝和分支gydF4y2Ba

我们可以看到在图gydF4y2Ba 30.gydF4y2Ba,盖茨背后的速度是零在完全封闭的状态,和速度值分支入口处5 - 6米/秒。速度阀输出增加了约36%相比,当阀门完全打开。基于这些条件,分支阀的输出流量估计约为5.3 mgydF4y2Ba3gydF4y2Ba/ s。gydF4y2Ba

速度完全关闭门的分支领域的变化模式,100%的分支阀打开。gydF4y2Ba

 3.5.6。调查频道和槽的空化现象gydF4y2Ba

空泡沿着通道检查通常是基于无量纲数称为空化指数(gydF4y2Ba σgydF4y2Ba )。gydF4y2Ba

空化指数局部压力和流体速度的函数,这个指标的临界值是沿着通道(0.2 - -0.25)和0.2槽内(gydF4y2Ba 38gydF4y2Ba]。gydF4y2Ba

在底部出口渠道,导叶开度为100%时,最大速度,因此,最大放电容量是创建的频道。随后,空化指数的研究和控制在100%被认为是这个频道。自从Sardab大坝将建在海拔2500米的高空,大约20摄氏度的环境温度,计算空化指数,我们会gydF4y2Ba (8)gydF4y2Ba σgydF4y2Ba =gydF4y2Ba PgydF4y2Ba 腹肌gydF4y2Ba −gydF4y2Ba PgydF4y2Ba vgydF4y2Ba 1gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ρgydF4y2Ba 。gydF4y2Ba VgydF4y2Ba 2gydF4y2Ba ,gydF4y2Ba 在哪里gydF4y2Ba σgydF4y2Ba 空化指数;gydF4y2Ba PgydF4y2Ba 腹肌gydF4y2Ba =gydF4y2Ba PgydF4y2Ba 我gydF4y2Ba +gydF4y2Ba PgydF4y2Ba 酒吧gydF4y2Ba ,在这gydF4y2Ba PgydF4y2Ba 酒吧gydF4y2Ba (Pa)是大气压力;和gydF4y2Ba PgydF4y2Ba 我gydF4y2Ba (Pa)和相对压力gydF4y2Ba PgydF4y2Ba vgydF4y2Ba 是液体的蒸汽压。gydF4y2Ba

方程(gydF4y2Ba 9gydF4y2Ba)是用来计算大气压力:gydF4y2Ba (9)gydF4y2Ba PgydF4y2Ba 酒吧gydF4y2Ba =gydF4y2Ba PgydF4y2Ba 一个gydF4y2Ba 1gydF4y2Ba −gydF4y2Ba BgydF4y2Ba ZgydF4y2Ba TgydF4y2Ba 0gydF4y2Ba ggydF4y2Ba /gydF4y2Ba RgydF4y2Ba BgydF4y2Ba 。gydF4y2Ba

上述方程的参数的值定义如下:gydF4y2Ba (10)gydF4y2Ba ggydF4y2Ba RgydF4y2Ba BgydF4y2Ba =gydF4y2Ba 5.26gydF4y2Ba ,gydF4y2Ba TgydF4y2Ba 0gydF4y2Ba =gydF4y2Ba 293.1gydF4y2Ba KgydF4y2Ba ,gydF4y2Ba BgydF4y2Ba =gydF4y2Ba 0.0065gydF4y2Ba KgydF4y2Ba 米gydF4y2Ba ,gydF4y2Ba ZgydF4y2Ba =gydF4y2Ba 2500年gydF4y2Ba 米gydF4y2Ba ,gydF4y2Ba PgydF4y2Ba 一个gydF4y2Ba =gydF4y2Ba 海平面气压gydF4y2Ba 。gydF4y2Ba

因此,的价值gydF4y2Ba PgydF4y2Ba 酒吧gydF4y2Ba =gydF4y2Ba 74年gydF4y2Ba kPagydF4y2Ba 被认为是使用方程(gydF4y2Ba 9gydF4y2Ba)。gydF4y2Ba

应该指出的是,液体蒸汽压等于gydF4y2Ba PgydF4y2Ba vgydF4y2Ba =gydF4y2Ba 2.337gydF4y2Ba kPagydF4y2Ba 在这个温度。空化指数的帮助下可以计算压力和速度值在不同的条件下使用以上参数的值。gydF4y2Ba

计算的相对压力值gydF4y2Ba PgydF4y2Ba 我gydF4y2Ba 根据特定的值的51.35头,速度gydF4y2Ba VgydF4y2Ba 我gydF4y2Ba 可以计算的帮助下以下方程。损失值和基于参考截面的速度值gydF4y2Ba VgydF4y2Ba ,也有一个截面gydF4y2Ba 一个gydF4y2Ba =gydF4y2Ba 1.25gydF4y2Ba ×gydF4y2Ba 1.10gydF4y2Ba =gydF4y2Ba 1.375gydF4y2Ba 米gydF4y2Ba 2gydF4y2Ba :gydF4y2Ba (11)gydF4y2Ba PgydF4y2Ba 我gydF4y2Ba =gydF4y2Ba ρgydF4y2Ba ggydF4y2Ba HgydF4y2Ba −gydF4y2Ba VgydF4y2Ba ogydF4y2Ba 2gydF4y2Ba 2gydF4y2Ba ggydF4y2Ba −gydF4y2Ba ∑gydF4y2Ba hgydF4y2Ba 我gydF4y2Ba ,gydF4y2Ba ∑gydF4y2Ba hgydF4y2Ba ogydF4y2Ba =gydF4y2Ba ∑gydF4y2Ba ξgydF4y2Ba VgydF4y2Ba 2gydF4y2Ba 2gydF4y2Ba ggydF4y2Ba 。gydF4y2Ba

考虑到流头= 51.35门槽的两个部分和输出部分,将会有一个流(输出部分,相对压力等于零,速度等于18.9 m / s)。gydF4y2Ba

在输出部分,gydF4y2Ba (12)gydF4y2Ba σgydF4y2Ba =gydF4y2Ba PgydF4y2Ba 我gydF4y2Ba +gydF4y2Ba PgydF4y2Ba 酒吧gydF4y2Ba −gydF4y2Ba PgydF4y2Ba vgydF4y2Ba 1gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ρgydF4y2Ba 。gydF4y2Ba VgydF4y2Ba 2gydF4y2Ba ≈gydF4y2Ba 0.4。gydF4y2Ba

同时,服务门槽的空化值,作为空化发生的最关键点之一,计算如下:gydF4y2Ba (13)gydF4y2Ba PgydF4y2Ba 我gydF4y2Ba =gydF4y2Ba ρgydF4y2Ba ggydF4y2Ba HgydF4y2Ba −gydF4y2Ba VgydF4y2Ba ogydF4y2Ba 2gydF4y2Ba 2gydF4y2Ba ggydF4y2Ba −gydF4y2Ba ∑gydF4y2Ba hgydF4y2Ba 我gydF4y2Ba ≈gydF4y2Ba 71年gydF4y2Ba kPagydF4y2Ba ,gydF4y2Ba σgydF4y2Ba =gydF4y2Ba PgydF4y2Ba 我gydF4y2Ba +gydF4y2Ba PgydF4y2Ba 酒吧gydF4y2Ba −gydF4y2Ba PgydF4y2Ba vgydF4y2Ba 1gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ρgydF4y2Ba 。gydF4y2Ba VgydF4y2Ba 2gydF4y2Ba =gydF4y2Ba ≈gydF4y2Ba 0.79。gydF4y2Ba

门槽应该检查使用槽特征以检查空化的可能性。考虑到槽的收缩角是1:12gydF4y2Ba WgydF4y2Ba /gydF4y2Ba DgydF4y2Ba =gydF4y2Ba 140年gydF4y2Ba /gydF4y2Ba 162年gydF4y2Ba =gydF4y2Ba 0.86gydF4y2Ba 和gydF4y2Ba ΔgydF4y2Ba /gydF4y2Ba WgydF4y2Ba =gydF4y2Ba 25gydF4y2Ba /gydF4y2Ba 140年gydF4y2Ba =gydF4y2Ba 0.18gydF4y2Ba ,允许索引值约为0.2,这个指数是0.79 Sardab大坝的闸门槽。因此,气蚀槽(不会发生gydF4y2Ba 39gydF4y2Ba- - - - - -gydF4y2Ba 42gydF4y2Ba]。gydF4y2Ba

坝底媒体包含阀门和泵在大坝运行和安全中扮演着重要的角色,他们允许控制水面高程低于溢洪道的水平。空泡的形成由于概率集中漩涡太高了。分支入口附近的涡量分布表现出高值由于出口底部的圆柱形状,导致大量流动分离的gydF4y2Ba 43gydF4y2Ba,gydF4y2Ba 44gydF4y2Ba]。瞬变流动特性与移动门成功抓住了放电特性以及部分开放的大门。部分空缺,水流在门口唇在高速和持续的下游的空气门,这可能会导致损失由于气蚀和振动。Flow-3D可能有用的结果表明,计算空气坝底门店的需求。在研究中,现有设计的适用性验证。部分门打开,放电随导叶开度的平方,而对于一个完全打开门,排放随gydF4y2Ba HgydF4y2Ba1/2gydF4y2Ba。新设施的可能性Flow-3D流模式识别和计算的压力和速度场应该有助于设计曝气系统。gydF4y2Ba

 4所示。结论gydF4y2Ba

在目前的研究中,Sardab大坝的底部出口和豪厄尔他们阀在不同的开放条件下的水力性能,包括20、40、60、80和100%,与不同的流速,调查。这些检查包括服务门的单操作,服务门和豪厄尔他们阀同时操作,和豪厄尔他们阀的单一操作。结果提出了不同的打开速度条件值有100%在门下面的部分是关于18 m / s,和大门下的最大速度为40% = 23.1 m / s。门槽中的速度值大约是2 - 3 m / s,在这种情况下没有观察到流分离和不良的循环流槽。这种情况适用于所有门的开口。速度值分支入口处变量2 - 5 m / s,没有流动分离和其他不利条件下发生在这一领域,经过和速度值增加分支入口。没有不良的压力分布变化沿着隧道和在城门口区域,而且没有大幅减少在这方面的压力。应该注意的是,没有不良的水力现象在这些部分,包括流动分离和当地的漩涡。50%的盖茨在门的上部区域,减少所需的压力值,并在两个门之间的区域,压力值降低。此外,与中控室内的安装,减少了空泡在这一领域的可能性。 For fully closed gate mode and 100% open branch valve, no circular flow and undesirable pressure changes are created in the tunnel. Also, the velocity behind the valves in the fully closed state is zero, and at the branch entrance point, it is 5-6 m/s. At the valve output, the velocity has increased by about 36% compared to when the gates are fully open. Based on these conditions, the branch valve’s output flow rate is estimated to be about 5.3 m/s. The presented results show that due to the bottom outlet operation in the reservoir’s maximum head condition, the probability of cavitation in the area between the two gates is very high. This analysis suggests that numerical modeling with the Flow-3D can be helpful for the design of this kind of hydraulic works.

 数据可用性gydF4y2Ba

所有数据用来支持这个研究的发现包括在本文中。gydF4y2Ba

 信息披露gydF4y2Ba

这项研究没有得到具体拨款资助机构在公众,商业,或非营利部门。gydF4y2Ba

 的利益冲突gydF4y2Ba

作者宣称没有利益冲突有关的出版。gydF4y2Ba

 萨拉查gydF4y2Ba F。gydF4y2Ba San-MaurogydF4y2Ba J。gydF4y2Ba CeliguetagydF4y2Ba m·A。gydF4y2Ba 奥纳gydF4y2Ba E。gydF4y2Ba 底部的空气需求估计网点与粒子有限元方法gydF4y2Ba 计算粒子力学gydF4y2Ba 2016年gydF4y2Ba 4gydF4y2Ba 3gydF4y2Ba 345年gydF4y2Ba 356年gydF4y2Ba 10.1007 / s40571 - 016 - 0117 - 4gydF4y2Ba 2 - s2.0 - 85026442347gydF4y2Ba Zounemat-KermanigydF4y2Ba M。gydF4y2Ba 朔尔茨gydF4y2Ba M。gydF4y2Ba 计算空气对大坝门店需求使用Takagi-Sugeno模型gydF4y2Ba 水gydF4y2Ba 2013年gydF4y2Ba 5gydF4y2Ba 3gydF4y2Ba 1441年gydF4y2Ba 1456年gydF4y2Ba 10.3390 / w5031441gydF4y2Ba 2 - s2.0 - 84885336818gydF4y2Ba 陈gydF4y2Ba g . H。gydF4y2Ba 王gydF4y2Ba g . Y。gydF4y2Ba 胡gydF4y2Ba c . L。gydF4y2Ba 黄gydF4y2Ba B。gydF4y2Ba 张gydF4y2Ba m D。gydF4y2Ba 观察和测量不稳定空泡流动使用同步采样的方法gydF4y2Ba 实验流体gydF4y2Ba 2015年gydF4y2Ba 56gydF4y2Ba 2gydF4y2Ba 10.1007 / s00348 - 015 - 1896 - 8gydF4y2Ba 2 - s2.0 - 84961331452gydF4y2Ba KaralargydF4y2Ba M。gydF4y2Ba CavuşligydF4y2Ba M。gydF4y2Ba 地震震中的距离在3 d地震损伤性能的CG大坝gydF4y2Ba 地震和结构gydF4y2Ba 2020年gydF4y2Ba 18gydF4y2Ba 2gydF4y2Ba 201年gydF4y2Ba 213年gydF4y2Ba MoradigydF4y2Ba M。gydF4y2Ba AghajanzadehgydF4y2Ba s M。gydF4y2Ba MirzabozorggydF4y2Ba H。gydF4y2Ba AlimohammadigydF4y2Ba M。gydF4y2Ba 水下爆炸及其对拱坝的非线性行为的影响gydF4y2Ba 耦合系统力学gydF4y2Ba 2018年gydF4y2Ba 7gydF4y2Ba 3gydF4y2Ba 333年gydF4y2Ba 351年gydF4y2Ba 10.12989 / csm.2018.7.3.333gydF4y2Ba 2 - s2.0 - 85045316573gydF4y2Ba KaralargydF4y2Ba M。gydF4y2Ba CavuşligydF4y2Ba M。gydF4y2Ba 检查3 d长期CFR大坝使用特殊材料的粘塑性的行为模型gydF4y2Ba 地质力学和工程gydF4y2Ba 2019年gydF4y2Ba 17gydF4y2Ba 2gydF4y2Ba 119年gydF4y2Ba 131年gydF4y2Ba KaralargydF4y2Ba M。gydF4y2Ba CavuşligydF4y2Ba M。gydF4y2Ba 评价三维非线性地震行为Ilısu CFR大坝下far-fault地面运动gydF4y2Ba 土木工程的发展gydF4y2Ba 2019年gydF4y2Ba 2019年gydF4y2Ba 15gydF4y2Ba 7358710gydF4y2Ba 10.1155 / 2019/7358710gydF4y2Ba 2 - s2.0 - 85060514237gydF4y2Ba YaminigydF4y2Ba o . A。gydF4y2Ba KavianpourgydF4y2Ba m·R。gydF4y2Ba 穆萨维gydF4y2Ba s . H。gydF4y2Ba 波前和破旧ACB垫下颗粒和土工布过滤的条件gydF4y2Ba 海洋Georesources &印度尼西亚gydF4y2Ba 2017年gydF4y2Ba 36gydF4y2Ba 8gydF4y2Ba 895年gydF4y2Ba 906年gydF4y2Ba 10.1080 / 1064119 x.2017.1397068gydF4y2Ba 2 - s2.0 - 85035125463gydF4y2Ba SpeerligydF4y2Ba J。gydF4y2Ba 海格gydF4y2Ba w·H。gydF4y2Ba 空气流在底部gydF4y2Ba 加拿大土木工程杂志》上gydF4y2Ba 2000年gydF4y2Ba 27gydF4y2Ba 3gydF4y2Ba 454年gydF4y2Ba 462年gydF4y2Ba 10.1139 / l99 - 087gydF4y2Ba 2 - s2.0 - 0034209601gydF4y2Ba FazelabdolabadigydF4y2Ba B。gydF4y2Ba GolestangydF4y2Ba m . H。gydF4y2Ba 对贝叶斯量化的渗透率在微尺度多孔结构——微网络数据库gydF4y2Ba 高科技和创新》gydF4y2Ba 2020年gydF4y2Ba 1gydF4y2Ba 4gydF4y2Ba 148年gydF4y2Ba 160年gydF4y2Ba 10.28991 / hij - 2020 - 01 - 04 - 02gydF4y2Ba TouaibigydF4y2Ba R。gydF4y2Ba KotengydF4y2Ba H。gydF4y2Ba BoydakgydF4y2Ba O。gydF4y2Ba 参数的研究使用不同的液体有机郎肯循环gydF4y2Ba 新兴科学杂志gydF4y2Ba 2020年gydF4y2Ba 4gydF4y2Ba 2gydF4y2Ba 122年gydF4y2Ba 128年gydF4y2Ba 10.28991 / esj - 2020 - 01216gydF4y2Ba InnellagydF4y2Ba G。gydF4y2Ba 罗杰斯gydF4y2Ba p。gydF4y2Ba 融合艺术与工程所带来的好处gydF4y2Ba 高科技和创新》gydF4y2Ba 2021年gydF4y2Ba 2gydF4y2Ba 1gydF4y2Ba 29日gydF4y2Ba 37gydF4y2Ba 10.28991 / hij - 2021 - 02 - 01 - 04gydF4y2Ba 萨达特HelbargydF4y2Ba M。gydF4y2Ba Parvaresh RizigydF4y2Ba 一个。gydF4y2Ba FarhoudigydF4y2Ba J。gydF4y2Ba 穆罕默gydF4y2Ba 一个。gydF4y2Ba 三维流动模拟改善的设计和操作大坝的底部gydF4y2Ba 阿拉伯地球科学杂志》gydF4y2Ba 2021年gydF4y2Ba 14gydF4y2Ba 2gydF4y2Ba 10.1007 / s12517 - 020 - 06378 - 4gydF4y2Ba SreedhargydF4y2Ba b K。gydF4y2Ba 艾伯特gydF4y2Ba 美国K。gydF4y2Ba 潘迪特gydF4y2Ba 答:B。gydF4y2Ba 空泡腐蚀:理论和measurements-a审查gydF4y2Ba 穿gydF4y2Ba 2017年gydF4y2Ba 372 - 373gydF4y2Ba 177年gydF4y2Ba 196年gydF4y2Ba 10.1016 / j.wear.2016.12.009gydF4y2Ba 2 - s2.0 - 85007200914gydF4y2Ba 森古普塔gydF4y2Ba a。R。gydF4y2Ba 古普塔gydF4y2Ba R。gydF4y2Ba BiswasgydF4y2Ba 一个。gydF4y2Ba 计算流体动力学分析蒙加丝的烹饪和干燥炉系统gydF4y2Ba 新兴科学杂志gydF4y2Ba 2019年gydF4y2Ba 3gydF4y2Ba 5gydF4y2Ba 285年gydF4y2Ba 292年gydF4y2Ba 10.28991 / esj - 2019 - 01191gydF4y2Ba NtintakisgydF4y2Ba 我。gydF4y2Ba StavroulakisgydF4y2Ba g . E。gydF4y2Ba PlakiagydF4y2Ba N。gydF4y2Ba 使用拓扑优化的3 d打印技术在产品设计的过程gydF4y2Ba 高科技和创新》gydF4y2Ba 2020年gydF4y2Ba 1gydF4y2Ba 4gydF4y2Ba 161年gydF4y2Ba 171年gydF4y2Ba 10.28991 / hij - 2020 - 01 - 04 - 03gydF4y2Ba 杨gydF4y2Ba J。gydF4y2Ba 腾gydF4y2Ba P。gydF4y2Ba 谢gydF4y2Ba Q。gydF4y2Ba 李gydF4y2Ba 年代。gydF4y2Ba 理解水和空气流发泄spillway-a案例研究的特点gydF4y2Ba 水gydF4y2Ba 2020年gydF4y2Ba 12gydF4y2Ba 8gydF4y2Ba 2106年gydF4y2Ba 10.3390 / w12082106gydF4y2Ba 克莱默gydF4y2Ba K。gydF4y2Ba 海格gydF4y2Ba w·H。gydF4y2Ba 小gydF4y2Ba 纳米比亚。gydF4y2Ba 空气浓度槽溢洪道的发展gydF4y2Ba 《水利工程gydF4y2Ba 2006年gydF4y2Ba 132年gydF4y2Ba 9gydF4y2Ba 908年gydF4y2Ba 915年gydF4y2Ba 10.1061 /(第3期)0733 - 9429 (2006)132:9 (908)gydF4y2Ba 2 - s2.0 - 33747331581gydF4y2Ba 列文gydF4y2Ba J。gydF4y2Ba 液压闸门和阀门:在自由表面流和淹没gydF4y2Ba 2001年gydF4y2Ba 英国伦敦gydF4y2Ba 托马斯德福gydF4y2Ba ShilpakargydF4y2Ba R。gydF4y2Ba 华gydF4y2Ba Z。gydF4y2Ba ManandhargydF4y2Ba B。gydF4y2Ba 数值模拟在隧道溢洪道Jingping-I水电项目有四个限gydF4y2Ba IOP会议系列:地球和环境科学gydF4y2Ba 2017年gydF4y2Ba 82年gydF4y2Ba 012013年gydF4y2Ba 10.1088 / 1755 - 1315/82/1/012013gydF4y2Ba 2 - s2.0 - 85029469483gydF4y2Ba DaneshmandgydF4y2Ba F。gydF4y2Ba YazdigydF4y2Ba 大肠。gydF4y2Ba 拉里gydF4y2Ba y . B。gydF4y2Ba 水力模型试验Sivand大坝的底部出口gydF4y2Ba 第五ISME /圆柱学报》国际会议在流体力学和空气动力学gydF4y2Ba 2007年gydF4y2Ba 希腊雅典gydF4y2Ba 269年gydF4y2Ba 275年gydF4y2Ba VischergydF4y2Ba d . L。gydF4y2Ba HagnergydF4y2Ba w·H。gydF4y2Ba 大坝水力学gydF4y2Ba 1999年gydF4y2Ba 英国伦敦gydF4y2Ba 约翰威利& SonsgydF4y2Ba KavianpourgydF4y2Ba m·R。gydF4y2Ba 导向板的接续流下游,包括空气喷射的影响gydF4y2Ba 1997年gydF4y2Ba 英国曼彻斯特gydF4y2Ba 曼彻斯特大学gydF4y2Ba 博士论文gydF4y2Ba KhazaeigydF4y2Ba 年代。gydF4y2Ba EllahigydF4y2Ba R。gydF4y2Ba NiknamigydF4y2Ba E。gydF4y2Ba ShahverdygydF4y2Ba K。gydF4y2Ba 水力分析和空化空化现象的物理和数值模型底部出口在鲁德巴两座大坝,Lorestan省、伊朗gydF4y2Ba 美国第36 IAHR世界大会gydF4y2Ba 2015年gydF4y2Ba 海牙,荷兰gydF4y2Ba 阮gydF4y2Ba S.-P。gydF4y2Ba 吴gydF4y2Ba 黄永发。gydF4y2Ba 吴gydF4y2Ba 裴伟伟。gydF4y2Ba 习gydF4y2Ba R.-Z。gydF4y2Ba 液压研究中控室内的隧道gydF4y2Ba 流体力学杂志gydF4y2Ba 2007年gydF4y2Ba 19gydF4y2Ba 3gydF4y2Ba 330年gydF4y2Ba 334年gydF4y2Ba 10.1016 / s1001 - 6058 (07) 60066 - 3gydF4y2Ba 2 - s2.0 - 34347363233gydF4y2Ba 吴gydF4y2Ba 黄永发。gydF4y2Ba 冯gydF4y2Ba S.-R。gydF4y2Ba 吴gydF4y2Ba 裴伟伟。gydF4y2Ba 阮gydF4y2Ba S.-P。gydF4y2Ba 工作闸门水力特性的部分开放,放电通道gydF4y2Ba 流体力学杂志gydF4y2Ba 2007年gydF4y2Ba 19gydF4y2Ba 4gydF4y2Ba 488年gydF4y2Ba 493年gydF4y2Ba 10.1016 / s1001 - 6058 (07) 60144 - 9gydF4y2Ba 2 - s2.0 - 34548173297gydF4y2Ba 李gydF4y2Ba J.-Y。gydF4y2Ba 高gydF4y2Ba Z.-X。gydF4y2Ba 吴gydF4y2Ba H。gydF4y2Ba 金gydF4y2Ba Z.-J。gydF4y2Ba 空化数值调查方法抑制在截止阀gydF4y2Ba 应用科学gydF4y2Ba 2020年gydF4y2Ba 10gydF4y2Ba 16gydF4y2Ba 5541年gydF4y2Ba 10.3390 / app10165541gydF4y2Ba NiksereshtgydF4y2Ba a . H。gydF4y2Ba TalebbeydokhtigydF4y2Ba N。gydF4y2Ba KhorshidigydF4y2Ba H。gydF4y2Ba 三维数值模拟空化和曝气系统的大坝gydF4y2Ba 流体工程gydF4y2Ba 2012年gydF4y2Ba 134年gydF4y2Ba 9gydF4y2Ba 091302年gydF4y2Ba 10.1115/1.4007214gydF4y2Ba 2 - s2.0 - 84865692513gydF4y2Ba 赫特gydF4y2Ba c·W。gydF4y2Ba 尼克尔斯gydF4y2Ba b D。gydF4y2Ba 体积的流体动力学(受到)方法的自由边界gydF4y2Ba 计算物理学杂志gydF4y2Ba 1981年gydF4y2Ba 39gydF4y2Ba 1gydF4y2Ba 201年gydF4y2Ba 225年gydF4y2Ba 10.1016 / 0021 - 9991 (81)90145 - 5gydF4y2Ba 2 - s2.0 - 0019367877gydF4y2Ba FerzigergydF4y2Ba j . H。gydF4y2Ba PerićgydF4y2Ba M。gydF4y2Ba 计算流体动力学的方法gydF4y2Ba 计算流体动力学的方法gydF4y2Ba 2002年gydF4y2Ba 柏林,德国gydF4y2Ba 施普林格gydF4y2Ba 10.1007 / 978-3-642-56026-2gydF4y2Ba 流科学,Inc ., FLOW-3D 10.0.1用户手册版本gydF4y2Ba Nazari-SharabiangydF4y2Ba M。gydF4y2Ba Nazari-SharabiangydF4y2Ba 一个。gydF4y2Ba KarakouziangydF4y2Ba M。gydF4y2Ba KaramigydF4y2Ba M。gydF4y2Ba 牺牲桩在河桥梁冲刷对策使用flow-3D数值研究gydF4y2Ba 土木工程杂志gydF4y2Ba 2020年gydF4y2Ba 6gydF4y2Ba 6gydF4y2Ba 1091年gydF4y2Ba 1103年gydF4y2Ba 10.28991 / cej - 2020 - 03091531gydF4y2Ba Al-MansorigydF4y2Ba n . j . H。gydF4y2Ba Al-ZubaidigydF4y2Ba l . s . A。gydF4y2Ba 幼发拉底河的一维水动力模型和水力参数的预测gydF4y2Ba 土木工程杂志gydF4y2Ba 2020年gydF4y2Ba 6gydF4y2Ba 6gydF4y2Ba 1074年gydF4y2Ba 1090年gydF4y2Ba 10.28991 / cej - 2020 - 03091530gydF4y2Ba 塔希尔gydF4y2Ba t M。gydF4y2Ba 安瓦尔gydF4y2Ba a . O。gydF4y2Ba 门唇取向对井底压力的影响系数坝隧道门gydF4y2Ba 阿拉伯科学与工程》杂志上gydF4y2Ba 2016年gydF4y2Ba 41gydF4y2Ba 12gydF4y2Ba 4927年gydF4y2Ba 4936年gydF4y2Ba 10.1007 / s13369 - 016 - 2202 - 7gydF4y2Ba 2 - s2.0 - 84991281701gydF4y2Ba DargahigydF4y2Ba B。gydF4y2Ba 流动特性的底部媒体与移动的大门gydF4y2Ba 液压研究期刊》的研究gydF4y2Ba 2010年gydF4y2Ba 48gydF4y2Ba 4gydF4y2Ba 476年gydF4y2Ba 482年gydF4y2Ba 10.1080 / 00221686.2010.507001gydF4y2Ba 2 - s2.0 - 77956593721gydF4y2Ba 王gydF4y2Ba H。gydF4y2Ba 气流和空气夹带在液压跳跃gydF4y2Ba 2014年gydF4y2Ba 澳大利亚布里斯班gydF4y2Ba 昆士兰大学的gydF4y2Ba 博士论文gydF4y2Ba NarayanangydF4y2Ba R。gydF4y2Ba 液压压力波动模型的作用gydF4y2Ba 诉讼的研讨会上造型水工建筑物的规模效应,国会IAHRgydF4y2Ba 1984年9月gydF4y2Ba 德国Essilngenam内卡河gydF4y2Ba "科技gydF4y2Ba h·T。gydF4y2Ba 空泡在降落伞和溢洪道、工程专著42gydF4y2Ba 1990年gydF4y2Ba 美国华盛顿特区gydF4y2Ba 水资源技术出版物,我们印刷所,垦务局gydF4y2Ba GhaderigydF4y2Ba M。gydF4y2Ba 康复计划的研究服务和紧急出口门底部1卡鲁恩河大坝gydF4y2Ba 大坝和水库gydF4y2Ba 2020年gydF4y2Ba 30.gydF4y2Ba 3gydF4y2Ba 91年gydF4y2Ba 96年gydF4y2Ba 10.1680 / jdare.20.00016gydF4y2Ba 刘gydF4y2Ba C。gydF4y2Ba 杨ydF4y2Ba Q。gydF4y2Ba 木gydF4y2Ba h·G。gydF4y2Ba 数值调查被动三维水翼空化控制使用一个插槽gydF4y2Ba 国际期刊热与流体流动的数值方法gydF4y2Ba 2019年gydF4y2Ba 30.gydF4y2Ba 7gydF4y2Ba 3585年gydF4y2Ba 3605年gydF4y2Ba 10.1108 /高频电炉- 05 - 2019 - 0395gydF4y2Ba SmajdorovagydF4y2Ba T。gydF4y2Ba NoskievičovagydF4y2Ba D。gydF4y2Ba 方法应用非参数控制图付诸实践gydF4y2Ba 新兴科学杂志gydF4y2Ba 2020年gydF4y2Ba 4gydF4y2Ba 4gydF4y2Ba 272年gydF4y2Ba 282年gydF4y2Ba 10.28991 / esj - 2020 - 01230gydF4y2Ba 戴维斯gydF4y2Ba c V。gydF4y2Ba 索伦森gydF4y2Ba k . E。gydF4y2Ba 应用水力学的手册gydF4y2Ba 1969年gydF4y2Ba 22gydF4y2Ba 纽约,纽约,美国gydF4y2Ba 麦格劳-希尔图书公司gydF4y2Ba 22gydF4y2Ba 66年gydF4y2Ba 首歌gydF4y2Ba X。gydF4y2Ba 刘gydF4y2Ba C。gydF4y2Ba 实验调查floor-attached涡对压力脉动的影响在轴流式泵井的底部gydF4y2Ba 可再生能源gydF4y2Ba 2020年gydF4y2Ba 145年gydF4y2Ba 2327年gydF4y2Ba 2336年gydF4y2Ba 10.1016 / j.renene.2019.07.125gydF4y2Ba NorizangydF4y2Ba t。gydF4y2Ba 里达gydF4y2Ba E。gydF4y2Ba al -gydF4y2Ba Z。gydF4y2Ba 增强层分配器的涡度减少泵井提高泵的效率gydF4y2Ba 可持续能源技术和评估gydF4y2Ba 2018年gydF4y2Ba 26gydF4y2Ba 28gydF4y2Ba 36gydF4y2Ba 10.1016 / j.seta.2017.06.001gydF4y2Ba 2 - s2.0 - 85020649632gydF4y2Ba