CN116050289B - 一种计算城市排水竖井射流冲击力的方法 - Google Patents
一种计算城市排水竖井射流冲击力的方法 Download PDFInfo
- Publication number
- CN116050289B CN116050289B CN202211605946.1A CN202211605946A CN116050289B CN 116050289 B CN116050289 B CN 116050289B CN 202211605946 A CN202211605946 A CN 202211605946A CN 116050289 B CN116050289 B CN 116050289B
- Authority
- CN
- China
- Prior art keywords
- flow
- shaft
- impact force
- falling
- water
- 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
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
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/20—Hydro energy
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Measuring Volume Flow (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
本发明公开了一种计算城市排水竖井射流冲击力的方法,属于城市排水技术领域,是针对不同跌落高度排水系统水平射流冲击竖井井壁及跌落射流冲击竖井底部跌水池,分别计算竖井井壁及跌水池底板的受力。竖井井壁上的冲击压力主要受入射流速的支配,考虑弹射的自由跌落水流可以更在准确预测水平冲击压力;竖井井底的冲击力主要由跌落射流到达跌水池时的速度决定,同时考虑跌落的环形流和弹射自由跌落水流,可以更好地计算竖井底部的冲击压力。本发明在满足管道系统排水要求的前提下,不仅可以计算城市排水竖井受到的冲击力,对深层隧道排水竖井内冲击力的计算同样适用,可为深层竖井的结构安全评估提供技术支撑,保证排水系统的安全稳定运行。
Description
技术领域
本发明涉及城市排水技术领域,尤其涉及一种计算城市排水竖井射流冲击力的方法。
背景技术
城市排水竖井经常被用作能量耗散结构,以减少陡峭地形或深隧道下的排水系统中的管道坡度,如排水系统中的雨水及污水系统。跌落式排水竖井是目前常用的一种类型,其流动特性、能量耗散、落差高度的影响以及跌水池已得到广泛研究。对于具有足够大的落差高度的竖井,水流在落下约5.0 m后可能会瓦解成约2 mm的小水滴。然而,对于中等高度的竖井,当跌落高度不够时,射流可能不会被瓦解,因此射流到达竖井底部跌水池时的形态可能是水滴或水柱。现有研究对竖井内射流冲击效应很有限,需要进一步关注这些结构的适当设计和保护。
本发明旨在提供一种能够计算城市排水射流冲击力的方法,有效获取竖井井壁及井底的受力,为评估竖井的结构安全提供技术支撑,保证排水系统的安全稳定运行。
发明内容
本发明为解决上述问题,而提出的一种计算城市排水竖井射流冲击力的方法。
为了实现上述目的,本发明采用了如下技术方案:
一种计算城市排水竖井射流冲击力的方法,包括以下步骤:
S1:建立城市排水竖井射流冲击力模拟的模型,所述模型包括进水管、进气管、跌落竖井、跌水池和出水管,所述进水管、跌落竖井、跌水池和出水管依次连接,所述跌落竖井顶部连接进气管,所述进气管与外面大气相通;
S2:调节所述进水管内的水流使水流冲击到跌落竖井井壁后,形成贴壁的环形流和弹射的自由跌落水流,所述自由跌落水流依据流量和竖井的跌落高度分别以水柱形式跌落的跌落水柱和水滴形式跌落的跌落水滴;
S3:分别计算水平射流撞击井壁的冲击力和跌落射流撞击井底的冲击力。
优选地,水平射流撞击竖井井壁之后划分为贴壁的环形流与自由跌落的射流,自由跌落射流与总入流的比值不大于0.2。
优选的,对于水平射流撞击井壁冲击力的计算方法如下:
基于方程(a1)求解竖井井壁的冲击力为:
其中,F1是竖井井壁受到的冲击力;ρ是入流的水平速度;Q0是入流流量;A0是水平入流的速度;α是反弹的自由跌落水流与总入流的比值。
优选的,对于跌落射流撞击井底的冲击力计算方法如下:
(1)首先计算贴壁环形流的最终流速
基于方程(b1)建立贴壁环形流的动量方程:
其中,Aw是环形流的横截面积;vw是环形流的速度;δ是环形流的厚度;Ds是竖井的直径;τ0是壁面与环形流的剪切应力;g是重力加速度;t是时间。
竖井井壁与环形流的剪切应力还可以表示为
其中,f是达西-魏斯巴赫系数。联立方程(b1)和(b2)可以求得计算竖井环形流跌落速度的数值方程
(2)计算自由弹射流的最终流速
基于方程(b4)求得计算竖井内水柱形势自由跌落水流速度
其中,Hs是竖井的跌落高度;vw2是竖井自由跌落水流速度。
水滴形式的自由跌落水流速度按照方程(b4)计算,达到6 m/s后取值为6 m/s。
基于方程(b5)求解竖井底部的冲击力,最终计算下降水流的冲击力:
其中,F2是竖井井底受到的冲击力;h为水垫层的深度。
与现有技术相比,本发明提供了一种计算城市排水竖井射流冲击力的方法,具备以下有益效果:
本发明的有益效果是:
(1)本发明通过分析竖井内部的水气相互作用,依据跌落高度的不同,将自由跌落水流划分为水柱跌落和水滴跌落。本发明提出的一个简化的方程,通过考虑反弹的水流,能够预测由于撞击在竖井井壁上的平均压力。
(2)竖井底部的冲击压力主要由下降的射流到达水池时的速度决定,这与水的流速和下降高度有关。本发明通过假设一定比例的环形流和反弹射流,提出的方程可以用来评估竖井底部的冲击压力,以确定是否有必要采取结构安全保护措施。
(3)本发明可基于具体的竖井设计尺寸及水流流量信息,在设计阶段实现对竖井冲击力的评估,不依赖于任何的测量设备。这对预防城市排水竖井的结构安全,保障城市排水系统稳定运行,具有十分重要的作用。
附图说明
图1为本发明的跌落竖井内部射流冲击竖井结构的示意图;
图2为本发明的上水平射流冲击竖井井壁的示意图;
图3为本发明的跌落射流冲击竖井底部跌水池的示意图;
图4为本发明实施例一中实测及计算的竖井井壁冲击力;
图5为本发明中实测及计算的竖井井底冲击力一;
图6为本发明中实测及计算的竖井井底冲击力二。
图号说明:
进水管1,进气管2,跌落竖井3,跌水池4,出水管5,入流总流6,环形流7,自由跌落水流8,跌落水柱9,跌落水滴10。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。
在本发明的描述中,需要理解的是,术语“上”、“下”、“前”、“后”、“左”、“右”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
实施例1:
一种计算城市排水竖井射流冲击力的方法,包括:进水管1、进气管2、跌落竖井3、跌水池4和出水管5。跌落竖井1的顶部连接进气管2,进气管2与外面大气相通。进水管2内的水流冲击到跌落竖井井壁后,形成贴壁的环形流7和弹射的自由跌落水流8。
图4给出了在不同流速下冲击点的计算压力和测量压力的比较,图例中红色标记表示平均值。由图可以看出,公式(a1)可以用来估计射流冲击竖井后井壁所产生的平均压力。被弹起的水与总入流的流量比和速度比为α和β。当忽略被弹起的水(即αβ=0)或考虑少量被弹起的水(即αβ=0.1)时,由射流冲击产生的滚筒高度明显增加,最高点接近流入的水面,为Q*=0.229。与流量为Q*=0.165的撞击相比,更多的水被弹开。因此,在计算冲击压力时,需要更大的αβ,以便与测量压力相等。忽略或考虑少量的弹跳水,计算结果在实验测量的范围内。当考虑更多的反弹水(即αβ=0.2),在大流量情况下,计算值比实验数据大。因此,本发明下的流量弹射系数小于0.2,这与实验观察结果一致。一般来说,在评估平均冲击压力时,建议考虑反弹的射流。
实施例2:
图5给出了H=3.38 m时,不同流速下井底的计算和测量的平均压力。计算和测量的压力之间的合理一致表明公式(b5)能够预测下降的射流冲击竖井跌水池所产生的平均压力。在自由跌落竖井中,射流开始破碎的起始点与射流的横截面积有关。流量越大,射流横截面积越大,射流破裂发生得越晚,达到的末端速度越高。流速的增加将导致下降水流的终端速度增加,这反过来又增加了井底的冲击压力。
实施例3:
以一个典型的竖井结构来说明如何使用本发明计算方法来评估竖井底部的冲击力。对于一个直径为1.2 m、落差为5.0 m的跌落式竖井,入流流量和速度分别为1.5 m3/s和5.0 m/s。通过使用公式(b3),当环形流到达跌水池时的速度约为9.7 m/s,没有达到终端速度。不考虑弹跳流,在落差高度为Hs=3.88米、跌水池为10 cm的跌落式竖井中,由落差射流产生的无量纲平均冲击压力可以达到0.564,最大脉动压力约为平均压力的1.58倍。最大压力和平均压力之间的比率随着竖井的跌落高度增加而增加。设计人员可以根据本发明提出的方程式评估井底的冲击压力,以评估是否需要采取结构安全保护措施。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。
Claims (2)
1.一种计算城市排水竖井射流冲击力的方法,其特征在于,包括以下步骤:
S1:建立城市排水竖井射流冲击力模拟的模型,所述模型包括进水管(1)、跌落竖井(3);
S2:调节所述进水管(1)内的水流使水流冲击到跌落竖井(3)井壁后,形成贴壁的环形流(7)和弹射的自由跌落水流(8),所述自由跌落水流(8)依据流量和竖井的跌落高度分别为以水柱形式跌落的跌落水柱(9)和水滴形式跌落的跌落水滴(10),弹射的所述自由跌落水流(8)与入流总流(6)的比例α小于0.2;
S3:分别计算水平射流撞击井壁的冲击力和跌落射流撞击井底的冲击力;
对于水平射流撞击井壁冲击力的计算方法如下:
基于方程(a1)求解竖井井壁的冲击力:
其中,F1是竖井井壁受到的冲击力;ρ是入流的水平速度;Q0是入流流量;A0是水平入流的速度;α是弹射的自由跌落水流(8)与入流总流(6)的比值;
对于跌落射流撞击井底的冲击力计算方法如下:
(1)首先计算贴壁环形流的最终流速;
(2)计算自由弹射流的最终流速;
(3)最终计算下降水流的冲击力;
计算贴壁环形流的最终流速的方法如下:
基于方程(b1)建立贴壁环形流的动量方程:
其中,vw是环形流的速度;δ是环形流的厚度;Ds是竖井的直径;τ0是壁面与环形流的剪切应力;g是重力加速度;t是时间;
竖井井壁与环形流的剪切应力还可以表示为
其中,f是达西-魏斯巴赫系数;
联立方程(b1)和(b2)可以求得计算竖井环形流跌落速度的数值方程:
计算自由弹射流的最终流速的方法如下:
基于方程(b4)求得竖井内水柱形式自由跌落水流速度
其中,Hs是竖井的跌落高度;vw2是竖井自由跌落水流速度;
水滴形式的自由跌落水流速度按照方程(b4)计算,达到6m/s后取值为6m/s;
最终计算下降水流的冲击力的方法如下:
基于方程(b5)求解竖井底部的冲击力:
其中,F2是竖井井底受到的冲击力;h为水垫层的深度。
2.根据权利要求1所述的一种计算城市排水竖井射流冲击力的方法,其特征在于:
所述模型还包括进气管(2)、跌水池(4)和出水管(5),所述进水管(1)、跌落竖井(3)、跌水池(4)和出水管(5)依次连接,所述跌落竖井(3)顶部连接进气管(2),所述进气管(2)与外面大气相通。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211605946.1A CN116050289B (zh) | 2022-12-12 | 2022-12-12 | 一种计算城市排水竖井射流冲击力的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211605946.1A CN116050289B (zh) | 2022-12-12 | 2022-12-12 | 一种计算城市排水竖井射流冲击力的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116050289A CN116050289A (zh) | 2023-05-02 |
CN116050289B true CN116050289B (zh) | 2023-10-10 |
Family
ID=86121470
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211605946.1A Active CN116050289B (zh) | 2022-12-12 | 2022-12-12 | 一种计算城市排水竖井射流冲击力的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116050289B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116956783B (zh) * | 2023-09-21 | 2023-12-08 | 武汉市规划研究院 | 一种雨水管渠汇流时间计算模型的构建方法及系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4036609A (en) * | 1971-10-18 | 1977-07-19 | Pircon Ladislav J | Endless belt impingement apparatus and method |
CN104278960A (zh) * | 2013-07-05 | 2015-01-14 | 中国石油天然气股份有限公司 | 水平井旋转射流器性能试验方法 |
CN206248343U (zh) * | 2016-11-07 | 2017-06-13 | 重庆市科学技术研究院 | 山地城市排水管道跌水井模拟实验装置 |
CN107893634A (zh) * | 2017-11-16 | 2018-04-10 | 成都理工大学 | 一种用于喷射钻井室内研究的多功能测试与实验平台 |
CN111397885A (zh) * | 2020-04-23 | 2020-07-10 | 中国石油大学(北京) | 海洋钻井隔水管软悬挂模拟试验装置 |
CN113605330A (zh) * | 2021-07-22 | 2021-11-05 | 湖北省水利水电规划勘测设计院 | 一种多级整流滑道竖井旋流泄洪洞 |
CN114970255A (zh) * | 2022-05-17 | 2022-08-30 | 武汉科技大学 | 一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法 |
CN115034155A (zh) * | 2022-06-13 | 2022-09-09 | 郑州大学 | 一种基于排水管道病害与内涝偶联关系的耦合方法 |
CN115100234A (zh) * | 2022-06-20 | 2022-09-23 | 宁波大学 | 一种测量固体颗粒物对竖井结构冲击力的方法及装置 |
-
2022
- 2022-12-12 CN CN202211605946.1A patent/CN116050289B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4036609A (en) * | 1971-10-18 | 1977-07-19 | Pircon Ladislav J | Endless belt impingement apparatus and method |
CN104278960A (zh) * | 2013-07-05 | 2015-01-14 | 中国石油天然气股份有限公司 | 水平井旋转射流器性能试验方法 |
CN206248343U (zh) * | 2016-11-07 | 2017-06-13 | 重庆市科学技术研究院 | 山地城市排水管道跌水井模拟实验装置 |
CN107893634A (zh) * | 2017-11-16 | 2018-04-10 | 成都理工大学 | 一种用于喷射钻井室内研究的多功能测试与实验平台 |
CN111397885A (zh) * | 2020-04-23 | 2020-07-10 | 中国石油大学(北京) | 海洋钻井隔水管软悬挂模拟试验装置 |
CN113605330A (zh) * | 2021-07-22 | 2021-11-05 | 湖北省水利水电规划勘测设计院 | 一种多级整流滑道竖井旋流泄洪洞 |
CN114970255A (zh) * | 2022-05-17 | 2022-08-30 | 武汉科技大学 | 一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法 |
CN115034155A (zh) * | 2022-06-13 | 2022-09-09 | 郑州大学 | 一种基于排水管道病害与内涝偶联关系的耦合方法 |
CN115100234A (zh) * | 2022-06-20 | 2022-09-23 | 宁波大学 | 一种测量固体颗粒物对竖井结构冲击力的方法及装置 |
Non-Patent Citations (2)
Title |
---|
吴建华 ; 杨涛 ; 沈洁艺 ; 任炜辰 ; 马飞 ; .大旋转角阶梯泄水道竖井水力特性研究.水动力学研究与进展(A辑).2018,(第02期),全文. * |
徐剑乔 ; .基于FLUENT的排水深隧入流竖井三维湍流仿真.中国给水排水.2019,(第21期), * |
Also Published As
Publication number | Publication date |
---|---|
CN116050289A (zh) | 2023-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116050289B (zh) | 一种计算城市排水竖井射流冲击力的方法 | |
Sene | Air entrainment by plunging jets | |
CN104790916B (zh) | 一种利用油套压力平衡法排除气井积液方法 | |
CN105574283A (zh) | 一种压裂后放喷油嘴尺寸的确定方法 | |
WO2006019977A3 (en) | Impact excavation system and method | |
Tutu | Pressure drop fluctuations and bubble-slug transition in a vertical two phase air-water flow | |
Tuna et al. | Scour profiles at downstream of cascades | |
CN104331598A (zh) | 一种气体钻井井筒内携岩状态快速判断方法 | |
CN107761654A (zh) | 梯级建设条件下的高坝施工中期防洪度汛方案决策方法 | |
CN106919772A (zh) | 高速水流明流泄水洞补气结构建造方法及补气结构 | |
Kramer et al. | On velocity estimations in highly aerated flows with dual-tip phase-detection probes-closure | |
CN111639434B (zh) | 一种页岩气水平井投球暂堵转向压裂的优化方法 | |
CN106639984A (zh) | 采气管柱、采气方法和涡流装置 | |
Kaminski et al. | Marginal stability of atmospheric eruption columns and pyroclastic flow generation | |
Chanson | Air-water flows in water engineering and hydraulic structures. Basic processes and metrology | |
Xu et al. | Turbulent flow and energy dissipation in plunge pool of high arch dam | |
Rein | Turbulent open-channel flows: Drop-generation and self-aeration | |
CN107103149A (zh) | 一种获取置换法压井关井期间压井液下落时间的方法 | |
Bennett et al. | Turbulent flow and bed pressure within headcut scour holes due to plane reattached jets | |
CN104596735B (zh) | 一种桥梁群的优化布置方法 | |
CN116362001A (zh) | 基于滑体质心高度势能计算高位滑坡入水速度及涌浪方法 | |
CN207187191U (zh) | 气、液和油三相分离装置 | |
Chen | Property experiments on the foam generator and its influencing factors during down-the-hole drilling | |
Lewis et al. | Air-water interactions that generate large water lift through vertical shafts in stormwater conduits | |
CN111967143B (zh) | 一种适于深水气井的临界携液流量预测方法 |
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 |