跌流竖井结构优化及其气压特性数值模拟.pdf
《跌流竖井结构优化及其气压特性数值模拟.pdf》由会员分享,可在线阅读,更多相关《跌流竖井结构优化及其气压特性数值模拟.pdf(7页珍藏版)》请在咨信网上搜索。
1、水利水电科技进展,2023,43(6)摇 Tel:02583786335摇 E鄄mail:jz 摇 http:/第 43 卷第 6 期Vol.43 No.6水 利 水 电 科 技 进 展Advances in Science and Technology of Water Resources2023 年 11 月Nov.2023作者简介:芦三强(1995),男,硕士,主要从事水工水力学和水气两相流研究。E鄄mail:2692249507 DOI:10.3880/j.issn.10067647.2023.06.004跌流竖井结构优化及其气压特性数值模拟芦三强1,乔时雨2(1.兰州理工大学能源与动
2、力工程学院,甘肃 兰州摇 730050;2.河海大学水文水资源学院,江苏 南京摇 210098)摘要:针对入流水舌引起的跌流竖井内负压过大问题,对前人研究中的原始竖井模型进行了结构优化设计。采用可实现的 k鄄着 湍流模型以及 VOF 两相流模型研究了不同入流量下原始竖井模型和优化竖井模型中的流态分布规律及气压分布情况,建立了可用于确定优化竖井模型中通气管道底部安装高程的水舌撞击模型,探讨了通气管道相对进水管的位置对优化模型中气压分布的影响。结果表明:进水管附近的跌落水舌束窄了井内气体的正常流通通道,水舌上方形成了大面积气体回流区,水舌下方补气不足,进水管附近负压突增现象明显;建立的水舌撞击模型
3、能够较准确地预测出水舌末端的撞击点位置,最大误差不超过 10%;得到的通气管道底部安装高程为 2.31 m,顶部高程略高于进水管上壁面;相比原始竖井模型,优化竖井模型能够使负压突增区的压差降低 70%以上,并且当通气管道与进水管水平中心延长线之间的夹角为30毅时,优化竖井模型所发挥的效果最好。关键词:跌流竖井;水舌;水气两相流;结构优化;气压分布;数值模拟中图分类号:TV672;TU992摇 摇 摇 文献标志码:A摇 摇 摇 文章编号:10067647(2023)06002406Structure optimization and numerical simulation of air pre
4、ssure characteristics of a plunging flow dropshaft/LUSanqiang1,QIAO Shiyu2(1.College of Energy and Power Engineering,Lanzhou University of Technology,Lanzhou 730050,China;2.College of Hydrology and Water Resources,Hohai University,Nanjing 210098,China)Abstract:Aiming at the problem of excessive nega
5、tive pressure in a dropshaft caused by the inflow water tongue,thestructure optimization design for the original dropshaft model was carried out.The realizable k鄄着 turbulence model and VOFtwo鄄phase flow model were used to study the flow regime and air pressure distribution in the original model and
6、theoptimized one under different inflow rates.A water tongue impact model was established to determine the installationelevation of the ventilation pipe bottom in the optimized model.The influence of the ventilation pipe position relative to theintake pipe on the air pressure distribution was discus
7、sed.The results show that the falling water tongue near the intakepipe blocks the normal air flow in the well,and large air recirculation zones are formed above the water tongue.The spacebelow the water tongue cannot be replenished in time,and the sudden increase of negative pressure near the intake
8、 pipe isobvious.The established water tongue impact model can accurately predict the impact point position at the water tongueend,and the maximum error is within 10%.The installation elevation of the ventilation pipe bottom is 2.31m,and the topelevation is slightly higher than the upper wall of the
9、water inlet pipe.The optimized shaft model can reduce the pressuredifference in the negative pressure surge zone by at least 70%,and when the angle between the ventilation pipe and thehorizontal center extension line of the inlet pipe is 30毅,the optimized dropshaft model has the best effect.Key word
10、s:plunging flow dropshaft;water tongue;gas鄄liquid two鄄phase flow;structure optimization;air pressuredistribution;numerical simulation摇 摇 跌流竖井作为水利工程及市政排水系统中广泛使用的输水建筑物之一,主要用于将上游海拔较高处管道中的水流转移到下游海拔较低的排水管道中1鄄2。不同于旋流竖井3、折板竖井4等其他类型竖井的进水方式,跌流竖井中水流从进水管处以水舌的形式直接自由跌落,水舌撞击井壁后会在进水管附近形成一个水体阻隔区域,造成竖井内通气不畅、负压过大以及卷吸
11、气量过多等问题5鄄8。跌流竖井内负压过大会引发一系列问题,如排污管道中气体流通的空间被加压,会导致排污系统中产生的一些臭气从一些地表开口逸散到外界环境中9鄄10;在极端暴雨天气下还可能会引发井喷,冲起的井盖在下落时可能会砸伤行人和引发交通事故等11鄄12。为了降低水舌引起的负压过大问题,芦三强等13在跌流竖井进水管附近布置了一个有旋分流隔板,目的是改变下落水流的流动方向,使水流围42水利水电科技进展,2023,43(6)摇 Tel:02583786335摇 E鄄mail:jz 摇 http:/绕隔板两侧旋转流动,同时隔板中间形成了一个不过水的空腔区域,能有效防止水舌束窄气体流通的空间。龚旭14
12、在进水管与立管交汇处上方放置了一个挡板,其目的是调节进水管中跌落水舌的入流角度,使水舌直接撞击到挡板上而不是管壁上,从而在进水管附近保留一定气体流通空间。以上关于跌流竖井进水管附近的结构优化设计都发生在竖井结构内部,尽管能起到一定的通气和调节气压的作用,但是由于占据了部分水体流通的空间,影响了竖井的过流能力,同时也不利于后续检修作业等。本文基于文献7中的竖井几何模型(原始竖井模型),采用数值模拟的方法研究了原始竖井模型进水管附近的流态与气压分布规律,分析了进水管附近负压突增的原因,提出了改善进水管附近通气状况的结构优化方案,建立了可以确定通气管道安装高程的理论模型,并评估了优化竖井模型的性能,
13、最后探讨了通气管道相对进水管的位置对优化竖井模型内气压分布的影响。1摇 竖井模型数值模拟1.1摇 数值模型跌流竖井内水流以自由跌落的形式运动,下落过程中水流会发生碰撞和破碎等,水流运动剧烈,同时下泄水流会夹带大量气流,水气混掺严重,属于典型的水气两相流。雷诺平均模型能够很好地模拟湍流运动,其中可实现的 k鄄着 湍流模型在湍动黏度中引入了与旋转和曲率有关的内容15鄄16,可以较好地模拟管道流、射流和混合流等多种复杂形式的流动16鄄17。VOF 模型能够很好地追踪两种流体的相界面位置,可以方便地观察竖井内水气两相的分布状态。本文采用可实现的 k鄄着 湍流模型和 VOF 模型进行跌流竖井内流态分布规
14、律及气压分布研究,两种模型的具体形式可参见文献16鄄17。1.2摇 几何模型及网格划分优化竖井模型几何结构如图 1(a)所示,该结构以竖井顶部中心为原点,z 轴向下为正方向进行建模。优化模型由原始竖井模型与添加在进水管正对面的通气管道组成,通气管道主要作为水舌阻隔区域气体流通的通道。通气管道主要包括水平管、竖直管和倾斜管,其直径均为 0.1 m,水平管的长度为0.36 m,竖直管的高度为 0.85 m,倾斜管与竖井井身的夹角为 70毅,设置倾斜管的目的主要是防止贴壁水流发生回流,堵塞通气管道的出口,原始竖井模型的几何参数可参见文献7。利用 ICEM CFD 软件进行网格划分。考虑到竖井结构为规
15、则的圆柱体结构,采用 O鄄grid 方法能够很好地将计算域划分为六面体结构化网格,同时也便于生成边界层网格。可实现的 k鄄着 湍流模型适用于雷诺数较大的湍流核心区域,而对于雷诺数较小的近壁面区域则配合壁面函数法进行求解,本文选用非平衡壁面函数法进行处理。将边界层网格的层数设置为 15 层,第一层网格高度设为 0.6 mm,每层按照指数方式增长,增长率为 1.2,划分完的网格如图 1(b)和(c)所示。图 1摇 优化竖井模型几何结构和网格划分(单位:m)1.3摇 边界条件及工况竖井进水管与进气管分开布置,进水管入口采用质量流入口,只允许水流通过;进气管入口通向外界大气,采用压力入口,竖井顶部密封
16、,气体全部从进气管进入竖井内部;出水口连通外界大气,设为压力出口;竖井壁面采用无滑移壁面边界。本文竖井除了结构和进水管中水量有变化外,再无其他变量,所以原始竖井模型和优化竖井模型均采用流量工况。选取 18.6、22.0、25.8、33.5、38.3、42.8、47.6 L/s等几组典型工况进行分析,其中部分流量工况与文献7中的入流量保持一致,主要是为了方便与文献7中的试验测量数据进行验证和比较。同时,在后续数据分析处理过程中,将入流量进行了无量纲化,目的是消除竖井体型对模拟结果的影响。无 量 纲 入 流 量 Q*w可 以 表 示 为 Q*w=Qw/gD5,其中 Qw为入流量,g 为重力加速度,
17、D 为竖井的直径。7 组流量工况入流量无量纲化后分别为 0.07、0.08、0.09、0.12、0.14、0.15、0.17。1.4摇 网格无关性验证为了验证网格数量对计算结果精度的影响,对23 万、36 万、47 万和 65 万共 4 种不同数量的网格在 Fluent 中采用相同的设置进行计算求解,选取无量纲入流量为 0.07 的工况进行网格无关性检验并观察竖井内的垂向气压变化,结果如图 2 所示。由52水利水电科技进展,2023,43(6)摇 Tel:02583786335摇 E鄄mail:jz 摇 http:/图 2 可知,随着网格数量的增加,计算结果的精度的确在提高,但是当网格数量增加
18、到 47 万和 65 万时,计算结果趋于稳定并不再变化,并且 47 万和 65 万网格下监测到的压力曲线与测量值吻合较好,也说明了数值模拟结果的准确性,因此选择 47 万网格进行后续的模拟。图 2摇 网格无关性验证图 4摇 竖井井身横截面上的压强分布2摇 结果与分析2.1摇 数值模拟准确性验证为了验证数值模拟结果的准确性,利用不同入流量下进气管中气体流速的模拟值与文献7中利用气压测量值计算出的流速值进行比较,结果如图 3所示。由图 3 可见,不同入流量下进气管中气体平均流速的模拟值与测量值吻合良好,最大误差不超过 6%,说明数值模拟的结果是可靠的。2.2摇 井身横向压强分布为了探究入流水舌对进
19、水管附近井身横截面上压强分布的影响,以无量纲入流量 0.09 和 0.12 为例,在 z=1 m(水舌影响区域上端)、z=1.5 m(水舌影响区域)和 z=2.5 m(水舌影响区域下端)处共选取 3 个截面进行分析,结果如图 4 所示。当无量纲入流量为 0.09 时,由图 4(a)可见,在 z=1 m 截面上,最大压差约为 3 Pa,说明在水舌上端区域,整个图 3摇 进气管气体平均流速对比截面上的气压变化很小,基本上是均匀分布的;由图 4(b)可见,整个截面上(z=1.5 m)的压强也基本相同,只有进水管正对面有一小部分区域表现为正压,这是由于水舌撞击到对面的井壁上造成的;由图 4(c)可见,
20、除了井壁以及进水管正对面位置处有压力的升降外,其他区域的压强也是均匀分布的,这是因为水舌撞击到井壁后,会有部分水流贴着壁面下落,因此壁面位置处压力会比竖井中央大。当无量纲入流量为 0.12 时,3 个位置处竖井截面上的压强分布与无量纲入流量为 0.09 时的情况相似。综上所述,同一入流量下,除了有水舌撞击的区域以及井壁上有水流附着的区域外,同一平面上的横向压强是均匀分布的。2.3摇 垂向气压分布原始竖井模型内的垂向气压分布如图 5 所示。由图 5 可知,竖井内相同位置上的气压(负压)随着入流量的增大而增大;同一流量下,从竖井顶部到底部气压整体呈现不断增大的趋势,说明竖井内气体卷吸现象严重,竖井
21、底部被不断加压。由图 5 还可以看到,在 z1(z=1.23 m)z2(z=1.85 m)区域内负压表现为突然增大,并且 z2处为竖井内负压最大的地方,同时入流量越大,z1 z2区域内的垂向气压梯度也 越 大,如 当 无 量 纲 入 流 量 为 0.09、0.12、62水利水电科技进展,2023,43(6)摇 Tel:02583786335摇 E鄄mail:jz 摇 http:/图 5摇 原始竖井模型内垂向气压分布0.15、0.17 时,z1和 z2两者间压差的绝对值分别为14.67、43.99、104、161 Pa,增幅很大。从进水管附近的流速分布规律可知,z1 z2区域的负压之所以会突增,
22、主要是因为水舌束窄了该区域内气体流通的通道,导致气体在该区域发生了大面积回流,水舌下方补气不足。图 6摇 进水管区域流速矢量与流线分布2.4摇 进水管附近流速矢量和流线分布规律本文重点探讨进水管附近的水舌运动及其影响,图 6 为进水管附近流速矢量与流线的分布。由图 6 可知,进气管中的气体进入竖井后先是与井壁相撞,一部分气体沿着井壁向上运动并在进气管上方形成了回流区,另一部分气体向下运动并在进水管附近受到了水舌的阻挡,在与水舌发生碰撞后,开始反向流动,在进水管与进气管之间形成了一个较大的回流区。不同的是,图 6(b)中在水舌与井壁碰撞出又形成了一个回流区,这主要是由于随着入流量的增大,水舌与井
23、壁撞击后形成的水体面积增大,导致该区域气体过流空间不足,所以气体发生了回流。而图 6(a)中由于水舌所占据的空间还不足以影响气体的流动,所以并没有在该处形成回流区。3摇 跌流竖井的结构优化3.1摇 通气管道安装高程的确定及气压对比分析通气管道的主要作用是连通被水舌阻隔的气体流通通道,方便水舌阻隔区域通气,因此确定通气管道的安装高程和位置就需要确定水舌阻隔区域的大体位置,水舌阻隔区域即水舌从进水管下方跌落到撞击到对面井壁上的区域。基于这个目的,本文提出一种水舌撞击模型,该模型旨在预测水舌末端撞击点的位置(水舌下落的起点位置已知,即进水管末端高程),建模时将水舌的跌落过程近似地看作是一种类平抛运动
24、。关于水舌末端撞击点高度的计算过程如下:a.确定水舌做类平抛运动的水平速度 v1。显然,v1就是进水管中水流的流速,即v1=Qw/A(1)式中 A 为进水管中水流的横截面积。b.确定水舌从进水管下落到撞击到井壁上所需的时间 T:T=x/v1(2)式中 x 为水舌做类平抛运动的水平位移。由于井壁的约束,水平位移最大为竖井的直径长。c.确定水舌从起始点下落到撞击点 C 的垂直距离 Y,即撞击点 C 的高度:Y=gT2/2(3)摇 摇 联合式(1)(2)(3)以及无量纲入流量 Q*w计算公式,Y 可以表示为Y=仔232D4ID31(Q*w)2(4)72水利水电科技进展,2023,43(6)摇 Tel
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 竖井 结构 优化 及其 气压 特性 数值 模拟
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。