CN109871601B - 切应力相等旋转圆盘电极等效模拟管道水流动状态的方法 - Google Patents
切应力相等旋转圆盘电极等效模拟管道水流动状态的方法 Download PDFInfo
- Publication number
- CN109871601B CN109871601B CN201910091020.7A CN201910091020A CN109871601B CN 109871601 B CN109871601 B CN 109871601B CN 201910091020 A CN201910091020 A CN 201910091020A CN 109871601 B CN109871601 B CN 109871601B
- Authority
- CN
- China
- Prior art keywords
- flow
- shear stress
- electrode
- pipeline
- 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
Images
Landscapes
- Pipeline Systems (AREA)
- Electrical Discharge Machining, Electrochemical Machining, And Combined Machining (AREA)
Abstract
本发明涉及再生水安全及腐蚀科学领域,为提出一种水流在圆盘电极工作表面与实际管道壁面切应力相等原则的以旋转圆盘电极等效模拟再生水管道流动状态的方法,建立等效模拟方程,并根据实际再生水管道流速,确定旋转圆盘电极转速范围,本发明,切应力相等旋转圆盘电极等效模拟管道水流动状态的方法,第一步:建立再生水对实际管道壁面的切应力方程;第二步:分析圆盘电极旋转时附近溶液的流动状态;第三步:建立等效模拟方程;第四步:根据实际再生水管道管径及流速资料,将各参数代入等效模拟方程,经计算选取旋转圆盘电极转速范围。本发明主要应用于再生水安全及腐蚀场合。
Description
技术领域
本发明涉及再生水安全及腐蚀科学领域,更具体的说,是涉及一种基于水流在圆盘电极工作表面与实际管道壁面切应力相等原则的以旋转圆盘电极等效模拟再生水管道流动状态的方法。
背景技术
再生水作为一种新兴的水资源,与海水淡化、跨流域调水相比具有明显的优势。已有研究表明给水管道在长期输水过程中,由于同时存在着悬浮物沉淀、钙镁离子沉积、微生物富集和金属电化学腐蚀等多个物理、化学及微生物反应过程,不可避免地导致管道内壁形成结构各异、成分复杂的垢层。这些垢层将会影响供水水质和管道输水能力,不利于管网安全运行。目前,国内外关于管道结垢的研究多集中在给水管道。与自来水相比,再生水含有较多的无机和有机成分,水质构成更为复杂;再生水管网常采用树状管网布置,较之环状布置的给水管网,其流速降低,水力停留时间延长,管网输配过程中的物化和生化反应也更为剧烈。水质与水力因素导致再生水管道结垢规律、结垢速率等都与给水管道不同,因此对再生水管道结垢的研究很有必要。
国内研究在成垢动力学领域方面,大多是在静态条件下通过配置pH、碱度、硬度都较高的溶液进行试验,从而建立动力学方程。但静态条件无法模拟实际再生水管道的状态,无法得到管垢在流动状态的传质速率。因此,在动态条件下对成垢过程进行研究尤为重要。然而,目前并没有相关专利为采用旋转圆盘电极等效模拟实际管流状态提供一个科学有效的方法。
发明内容
为克服现有技术的不足,本发明旨在提出一种水流在圆盘电极工作表面与实际管道壁面切应力相等原则的以旋转圆盘电极等效模拟再生水管道流动状态的方法,建立等效模拟方程,并根据实际再生水管道流速,确定旋转圆盘电极转速范围。为此,本发明采取的技术方案是,切应力相等旋转圆盘电极等效模拟管道水流动状态的方法,步骤如下:
第一步:依托边界层的概念,确定基于切应力相等原则的等效模拟的合理性,建立再生水对实际管道壁面的切应力方程;
第二步:分析圆盘电极旋转时附近溶液的流动状态,建立由旋转电极带动的水流在圆盘电极工作表面形成的切应力方程;
第三步:根据切应力相等原则,以旋转圆盘电极等效模拟再生水管道流动状态,建立等效模拟方程;
第四步:根据实际再生水管道管径及流速资料,将各参数代入等效模拟方程,经计算选取旋转圆盘电极转速范围。
第一步:建立再生水对实际管道壁面的切应力方程,具体步骤如下:
对实际管道来说:
式中:τ0为管壁上的切应力,Pa;λ为沿程阻力系数;ρ为流体密度,kg/m3;v为断面平均流速,m/s;
对实际再生水管道来说,圆管中液流状态为紊流,故有:
式中:△为壁面当量粗糙度;d为管段管径。
第二步:建立由旋转电极带动的水流在圆盘电极工作表面形成的切应力方程,具体过程如下:
溶液的流动分解为三个方向:由于离心力的存在,溶液在径向以流速vr向外流动;由于溶液的黏性,在圆盘旋转时,溶液以切向流速向切向流动;在电极附近这种向外的溶液流动使得电极中心区溶液的压力下降,于是离电极表面较远的溶液向中心区补充,形成轴向流动,流速为vy;
对于用于电化学研究的旋转圆盘电极,重要的速度是径向流速vr与法向流速vy,当靠近圆盘表面y→0或r→0时有:
vy=-0.51ω3/2v-1/y2 (3)
vr=0.51ω3/2v-1/2ry (4)
式中:vy为法向流速,vr为径向流速,ω为RDE角速度,v为流体运动粘滞系数,r为距圆盘中心的径向距离,y为距圆盘表面的法向距离;
对旋转圆盘电极来说:
式中:τ‘0为电极表面切应力,μ为流体动力粘滞系数,v为流速,y为至圆盘表面的法向距离;
将式(5)与式(4)联立,得:
τ‘0=0.51μω3/2v-1/2r (6)
采用加权平均的方法进行计算平均切应力T:
T=0.51μω3/2v-1/2{0.2R×π×(0.2R)2+0.3R×π×4(0.4R)2-(0.2R)2]+0.5R×π×4(0.6R)2-(0.4R)2]+0.7R×π×[(0.8R)2-(0.6R)2]+0.9R×π×[R2-(0.8R)2]}
第三步:建立等效模拟方程
将式(7)、式(1)、式(2)联立,得等效模拟方程:
τ‘0=τ0
本发明的特点及有益效果是:
本发明基于切应力相等原则,以旋转圆盘电极等效模拟再生水管道流动状态,建立等效模拟方程。这种等效模拟相当于可以实时监测再生水管道管壁的传质状态变化,在测量精度、数据准确度以及对实际管流的还原度方面都有了很大地提高,对于保证再生水管道的水质安全和管网安全运行具有重要的理论意义和实用价值。
附图说明:
图1为本发明涉及到的旋转圆盘电极的实物图;
图2为旋转圆盘电极附近溶液的流动状态示意图。
具体实施方式
本发明的基于切应力相等原则的以旋转圆盘电极等效模拟再生水管道流动状态的方法,包括以下步骤:
第一步:依托边界层的概念,确定基于切应力相等原则的等效模拟的合理性,建立再生水对实际管道壁面的切应力方程;
第二步:分析圆盘电极旋转时附近溶液的流动状态,建立由旋转电极带动的水流在圆盘电极工作表面形成的切应力方程;
第三步:根据切应力相等原则,以旋转圆盘电极等效模拟再生水管道流动状态,建立等效模拟方程;
第四步:根据实际再生水管道管径及流速资料,将各参数代入等效模拟方程,经计算选取旋转圆盘电极转速范围。
根据本发明的以旋转圆盘电极等效模拟再生水管道流动状态的方法,结合附图进一步说明实施过程。具体实现过程如下:
第一步:建立再生水对实际管道壁面的切应力方程
对实际管流进行模拟,首先应明确边界层的定义。1904年德国科学家普朗特认为黏性较小的流体,当其绕物体流动时,黏性的影响仅限于贴近物体的薄层中,而在这一薄层之外,黏性的影响可以忽略,这一薄层称为边界层。边界层内,速度由物体的零值沿物面外法线方向迅速增加至外部势流区速度,一般速度梯度很大,黏滞力不可忽略,是产生摩擦力的原因。依托边界层的概念,本发明认为,通过计算由旋转电极带动的水流在圆盘电极工作表面所形成的切应力,使之与实际再生水管道中再生水对管道壁面切应力相等,即遵循切应力相等的原则,方可实现旋转圆盘电极装置对实际管流状态的模拟。
对实际管道来说:
式中:τ0—管壁上的切应力,Pa;λ—沿程阻力系数;ρ—流体密度,kg/m3;v—断面平均流速,m/s。
对实际再生水管道来说,圆管中液流状态为紊流,故有:
式中:△—壁面当量粗糙度,一般铸铁管取△=0.15mm;d—管段管径,mm。
第二步:建立由旋转电极带动的水流在圆盘电极工作表面形成的切应力方程
旋转圆盘电极是电极理论与流体动力学结合的产物,其工作原理的基本要点是:物质传递和电流密度受控于电化学活性物质,而电化学活性物质的运动是按流体动力学规律进行的[于世才.旋转圆盘电极[J].材料开发与应用,1986,2:40-46.]。因此可从流体动力学的角度考虑,建立方程,通过选取合适的转速使旋转圆盘电极更好地模拟实际环境。旋转圆盘电极的实物图如图1所示。
由于溶液具有黏性,圆盘电极通过旋转带动附近的溶液发生流动。溶液的流动可分解为三个方向:由于离心力的存在,溶液在径向以流速vr向外流动;由于溶液的黏性,在圆盘旋转时,溶液以切向流速向切向流动;在电极附近这种向外的溶液流动使得电极中心区溶液的压力下降,于是离电极表面较远的溶液向中心区补充,形成轴向流动,流速为vy。旋转圆盘附近的液流情况示意图如图2所示,其中a为旋转圆盘附近的流速的矢量表示和,b为总流线的示意图。
对于用于电化学研究的旋转圆盘电极,重要的速度是径向流速vr与法向流速vy,当靠近圆盘表面y→0(或r→0)时有:
vy=-0.51ω3/2v-1/2y2 (3)
vr=0.51ω3/2v-1/2ry (4)
式中:vy—法向流速,m/s;vr—径向流速,m/s;ω—RDE角速度,rad/s;v—流体运动粘滞系数,m2/s;r—距圆盘中心的径向距离,m;y—距圆盘表面的法向距离,m。
对旋转圆盘电极来说:
式中:τ‘0—电极表面切应力,Pa;μ—流体动力粘滞系数,Pa·s;v—流速,m/s;y—至圆盘表面的法向距离,m。
本发明认为径向流速vr是形成圆盘电极表面切应力的主要因素,故将式(5)与式(4)联立,得:
τ‘0=0.51μω3/2v-1/2r (6)
由式(6)可看出,圆盘圆心处τ‘0的值为0,而圆盘盘面边缘处(r=R)τ‘0有最大值0.51μω3/2v-1/2R。故为求得整个圆面的平均切应力T,可采用加权平均的方法进行计算:
T=0.51μω3/2v-1/2{0.2R×π×(0.2R)2+0.3R×π×[(0.4R)2-(0.2R)2]+0.5R×π×[(0.6R)2-(0.4R)2]+0.7R×π×[(0.8R)2-(0.6R)2]+0.9R×π×[R2-(0.8R)2]}
第三步:建立等效模拟方程
将式(7)、式(1)、式(2)联立,可得等效模拟方程:
τ‘0=τ0
式中各项的物理意义及单位同上。
结合实际管段的管径及流速,将其带入式(8)中即可计算出对应的旋转圆盘电极的大致转速。值得注意的是,由于等效模拟公式是由加权平均法计算而得,同时在推导过程中也进行了某些理想化的假设,故本发明提出的等效模拟并非是使管段实际流速可以和旋转圆盘电极转速一一对应的精确模拟,而是能根据实际管流资料从而提出旋转圆盘电极转速范围的大致模拟。通过计算求得的转速范围,对后续试验设置旋转圆盘电极转速梯度、实现对实际再生水管道管流状态的模拟具有重要意义。
第四步:选取旋转圆盘电极转速范围
依据建立的等效模拟方程(8),圆盘电极的直径为5mm,公式中水质参数取25℃时水的物理性质数值,其中ρ=997.1kg/m3,μ=0.891×10-3Pa·s,v=0.894×10-6m2/s。参考表1中实际再生水管道管径及流速。
表1再生水管道管径、流速信息
经计算选取旋转圆盘转速范围为1000~3000r/min。
Claims (4)
1.一种切应力相等旋转圆盘电极等效模拟管道水流动状态的方法,其特征是,步骤如下:
第一步:依托边界层的概念,确定基于切应力相等原则的等效模拟的合理性,建立再生水对实际管道壁面的切应力方程;
第二步:分析圆盘电极旋转时附近溶液的流动状态,建立由旋转电极带动的水流在圆盘电极工作表面形成的切应力方程;
第三步:根据切应力相等原则,以旋转圆盘电极等效模拟再生水管道流动状态,建立等效模拟方程;
第四步:根据实际再生水管道管径及流速资料,将各参数代入等效模拟方程,经计算选取旋转圆盘电极转速范围,根据旋转圆盘电极转速范围设置旋转圆盘电极转速梯度,实现对实际再生水管道管流状态的模拟。
3.如权利要求1所述的切应力相等旋转圆盘电极等效模拟管道水流动状态的方法,其特征是,第二步:建立由旋转电极带动的水流在圆盘电极工作表面形成的切应力方程,具体过程如下:
溶液的流动分解为三个方向:由于离心力的存在,溶液在径向以流速vr向外流动;由于溶液的黏性,在圆盘旋转时,溶液以切向流速向切向流动;在电极附近这种向外的溶液流动使得电极中心区溶液的压力下降,于是离电极表面较远的溶液向中心区补充,形成轴向流动,流速为vy;
对于用于电化学研究的旋转圆盘电极,重要的速度是径向流速vr与法向流速vy,当靠近圆盘表面y→0或r→0时有:
vy=-0.51ω3/2ν-1/2y2 (3)
vr=0.51ω3/2ν-1/2ry (4)
式中:vy为法向流速,vr为径向流速,ω为RDE角速度,ν为流体运动粘滞系数,r为距圆盘中心的径向距离,y为距圆盘表面的法向距离;
对旋转圆盘电极来说:
式中:τ‘ 0为电极表面切应力,μ为流体动力粘滞系数,v为流速,y为至圆盘表面的法向距离;
将式(5)与式(4)联立,得:
τ‘0=0.51μω3/2v-1/2r (6)
采用加权平均的方法进行计算平均切应力T:
T=0.51μω3/2v-1/2{0.2R×π×(0.2R)2+0.3R×π×[(0.4R)2-(0.2R)2]+0.5R×π×[(0.6R)2-(0.4R)2]+0.7R×π×[(0.8R)2-(0.6R)2]+0.9R×π×[R2-(0.8R)2]}
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910091020.7A CN109871601B (zh) | 2019-01-30 | 2019-01-30 | 切应力相等旋转圆盘电极等效模拟管道水流动状态的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910091020.7A CN109871601B (zh) | 2019-01-30 | 2019-01-30 | 切应力相等旋转圆盘电极等效模拟管道水流动状态的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109871601A CN109871601A (zh) | 2019-06-11 |
CN109871601B true CN109871601B (zh) | 2023-05-23 |
Family
ID=66918349
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910091020.7A Active CN109871601B (zh) | 2019-01-30 | 2019-01-30 | 切应力相等旋转圆盘电极等效模拟管道水流动状态的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109871601B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101021465A (zh) * | 2007-03-09 | 2007-08-22 | 浙江理工大学 | 旋转式单相流冲蚀试验装置 |
CN102015544A (zh) * | 2008-03-06 | 2011-04-13 | 麦卡钦公司 | 电液和剪切空化径向逆流液体处理器 |
CN103412991A (zh) * | 2013-08-05 | 2013-11-27 | 大连理工大学 | 一种网箱对水流影响的模拟方法 |
CN105067504A (zh) * | 2015-07-28 | 2015-11-18 | 上海电气集团股份有限公司 | 一种高温熔融盐腐蚀模拟装置 |
CN105092460A (zh) * | 2015-06-23 | 2015-11-25 | 中石化石油工程设计有限公司 | 油水交替润湿腐蚀模拟装置及方法 |
CN106461603A (zh) * | 2016-05-23 | 2017-02-22 | 香港应用科技研究院有限公司 | 用于电镀液分析的静态圆盘电极 |
CN108827871A (zh) * | 2018-08-17 | 2018-11-16 | 河海大学 | 一种管式泥沙侵蚀试验装置中泥沙表面切应力确定方法 |
CN108955981A (zh) * | 2018-08-14 | 2018-12-07 | 北京航空航天大学 | 适用于旋转边界层壁面剪切应力测量方法和装置 |
-
2019
- 2019-01-30 CN CN201910091020.7A patent/CN109871601B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101021465A (zh) * | 2007-03-09 | 2007-08-22 | 浙江理工大学 | 旋转式单相流冲蚀试验装置 |
CN102015544A (zh) * | 2008-03-06 | 2011-04-13 | 麦卡钦公司 | 电液和剪切空化径向逆流液体处理器 |
CN103412991A (zh) * | 2013-08-05 | 2013-11-27 | 大连理工大学 | 一种网箱对水流影响的模拟方法 |
CN105092460A (zh) * | 2015-06-23 | 2015-11-25 | 中石化石油工程设计有限公司 | 油水交替润湿腐蚀模拟装置及方法 |
CN105067504A (zh) * | 2015-07-28 | 2015-11-18 | 上海电气集团股份有限公司 | 一种高温熔融盐腐蚀模拟装置 |
CN106461603A (zh) * | 2016-05-23 | 2017-02-22 | 香港应用科技研究院有限公司 | 用于电镀液分析的静态圆盘电极 |
CN108955981A (zh) * | 2018-08-14 | 2018-12-07 | 北京航空航天大学 | 适用于旋转边界层壁面剪切应力测量方法和装置 |
CN108827871A (zh) * | 2018-08-17 | 2018-11-16 | 河海大学 | 一种管式泥沙侵蚀试验装置中泥沙表面切应力确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109871601A (zh) | 2019-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Study on two operating conditions of a full-scale oxidation ditch for optimization of energy consumption and effluent quality by using CFD model | |
JP6987130B2 (ja) | ロータ組立体を有する流量計 | |
Ma et al. | CFD simulations of fiber-fiber interaction in a hollow fiber membrane bundle: Fiber distance and position matters | |
CN109871601B (zh) | 切应力相等旋转圆盘电极等效模拟管道水流动状态的方法 | |
CN107449731B (zh) | 同轴旋转腐蚀流动模拟系统 | |
Takahashi et al. | Laminar mixing in stirred tank agitated by an impeller inclined | |
García et al. | Influence of the Reynolds number on the thermal effectiveness of tubular heat exchanger subjected to electromagnetic field-based antifouling treatment in an open once-through seawater cooling system | |
Jackson | Coagulation in a rotating cylinder | |
Wang et al. | Erosion-corrosion of AISI 302 stainless steel sudden expansion pipes in high salt wastewater: effect of fluid flow on different positions of a sudden expansion pipe | |
Xie et al. | Numerical simulation and experimental investigation on the effect of a new suspended carrier filler on mass transfer in MBBR | |
Fordjour et al. | Effect of riser height on rotation uniformity and application rate of the dynamic fluidic sprinkler | |
CN206429701U (zh) | 一种方向控制阀 | |
CN111723447A (zh) | 一种动静腔流动微沟槽减阻设计方法 | |
Yang et al. | Multicomponent Water Effects on Rotating Machines Disk Erosion | |
Roy et al. | Interactive effects of water chemistry, hydrodynamics, and precipitated calcium carbonate causing erosion corrosion of copper in hot water recirculation systems: case study and experimental work | |
CN105444828B (zh) | 超声波流量计的整流器的设计方法 | |
CN104933219A (zh) | 一种水润滑推力轴承的全流态模拟参数获取方法 | |
Liu et al. | Experimental study of gas pressure drop in rotating packed bed with rotational‐stationary packing | |
Li et al. | Modelling and simulation of the airflow field in wedge-shaped zone during the high-speed grinding | |
CN104895796B (zh) | 一种离心泵气液两相流全特性分析方法 | |
Hudson et al. | Simulation of a high recirculation airlift reactor for steady-state operation | |
Feilin et al. | Design, simulation, and laboratory test of a single sludge drainage pipe for wastewater treatment sedimentation tank | |
CN111898071B (zh) | 高压淹没水射流沿程速度的计算方法 | |
CN212585560U (zh) | 一种两段式布水系统 | |
CN216717994U (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 |