污水不同排放流速对扩散器内海水清除能力数值模拟研究.docx

1、第37卷第6期2016年12月JournalofWaterwayandHarbor污水不同排放流速对扩散器内海水清除能力数值模跚究于航，白景峰，王心海(交通运输部天津水运工程科学研究所水路交通环境保护技术交通行业重点实验室，天津300456)摘要:在污水排海工程中,扩散器可提供给污水极大的初始稀释能力，对海洋生态及水质能够起到明显的保护作用。文章针对目前污水排海过程中的海水入侵问题,以数值模拟为主要技术手段，通过分析多种污水流速条件下的扩散器海水清除过程，掌握其流场变化机制。通过研究得出，在污水流量相对较大的情况F,污水清除时间较短，污水出流后与海水混合效果较明显。初步掌握一般形式下的扩散器海

2、水临界清除流速，为解决实际工程中的海水清除问题，提高工程运行效率提供实际依据。关键词:扩散器;海水清除;数值模拟中图分类号:X52;0242.1文献标识码:A文章编号:1005-8443(2016)06-0646-05在污水排海工程中，多孔扩散器的主要作用是通过它可以将污水均匀分散地排放到海洋环境水体中去,因此扩散器可提供给污水极大的初始稀释，能够对海洋生态及水质起到明显的保护作用，已成为污水排海工程成功的关键因素。1研究背景目前污水深海处置工程中的一个主要研究内容是扩散器海水清除和防止入侵设计，由于工程间歇排放等原因，在很多情况下,海水会倒灌侵入扩散器系统，严重影响扩散器排放效果，而要将入侵

3、到扩散器系统内的海水清除出来决非易事。现有的扩散器系统水力设计方法尚未将海水入侵与清除问题考虑在内，海水入侵问题在工程设计中未能得到很好解决。针对扩散器海水入侵机理与清除能力的研究，目前国内外已经有部分学者进行了相关研究工作。国外利用数值模拟计算进行该方向研究最早的是Mort.Burrows以及Larsen等人，针对波浪对于海水入侵与清除进行了分析。后续Tony.Doyle等人对该模型进行修正和完善，考虑两种流体密度差对流动的影响。扩散器系统海水入侵与清除从流动、物质的迁移、掺混与扩散的角度看是一个十分复杂的现象，一维数值模型不足以反映其内部流动细节。鉴于此，英国Belfast大学曾建立二维数

4、值模型对海水的清除过程进行模拟，认为在工程问题中所关心的是沿基本流向的流速变化过程，因此在计算中取轴向剖面进行二维计算。陈丽星等人对于应力模型在污水排海工程中的应用进行了研究。综上所述，数值模拟技术能反映扩散器系统内部的分层流动、浮力流动和掺混流动，有助于从定性的角度对流动特性进行理解。但对于不同流速下的污水清除能力研究目前还不多，鉴于各种实际因素的影响以及扩散器系统自身设计中的考虑不足,海水入侵现象在扩散器系统运行过程中时有发生且危害极大，因此有必要对污水不同排放流速对扩散器内海水的清除能力开展数值模拟研究，准确地预测设计的系统抵抗海水入侵和清除海水的能力，深刻理解其影响因素，以便在设计中采

5、取合理的工程措施,制定合理的运行模式，确保污水处置工程的运行效率C因此本文采用概化模型，通过数值模拟技术,将扩散器结构概念化设计,分析扩散器海水入侵过程的流收稿日期：2016-02-25；修回日期：2016-06-01基金项目：中央级公益性科研院所基本科研业务费(TKS150209,TKS140212)作者简介：于航(1980-),男,天津市人,副研究员，主要从事海洋环境科学与污水深海排放工程设计研究。Biography:YUHang(1980-),male,associateprofessor.图1排污管及流域整体网格模烈场变化情况,根据常规结构型式对不同流速条件下扩散器内部的海水清除能力及

6、流动特性进行分析，为提高扩散器海水清除能力提供理论依据'通过研究不同流速情况下扩散器内海水的清除过程.初步掌握海水与污水同时存在下的流动特性，对于提高污水深海排放的研究系统性和完整性是极其必要的.具有很强的学术价值与应用价值°2数学模型的建立及参数设置考虑到模型的结构复杂性及后续整体计算效率，二维模型Fig-1Overallgridmodel<ifdiffusionandwatershed的构建以四边形网格为主,整体模型方案的搭建及网格划分分别采用Gambil'Workbench及ICEM来完成。考虑到各个污水排海工程的扩散器结构各有不同.因此本文中才用概化模型

7、，将计算模型设置为最为常规的参数设置，选取每个上升管带有2个喷口的结构型式.扩放器模型长度及上升管数量选取则依据目前排海T程的常规水ktlOOOOmVd,设计5个上升管。(1) 排污管模型。儿何模型尺寸如下：排污管及计算流域长4(X)0mm.宽I254mm;主管道长3400mm直径54mm;支管数仙5、支管长100mm.宜径20mm;喷口长20nun,直径6.67mm、10min、15mm，与水平夹角0。、5。、10。；排污管主管顶部距离水面1000mm。(2) 流体模型。考虑到模型的结构夏杂性及后续推体计算效率，二维模型的构建以四边形网格为主.三维模型则以六面体为主。(3) 网格划分。排污管

8、二维扩散模型的平均网格尺度为24mm,采用四边形+三角形混合模式进行流域的网格划分，总体网格数为20.5万，三维喷口模型的平均网格尺度为总体网格数为40万,模型整体采用六面体模式进行网格划分,具体详情如图1所示。(4) 边界条件。本计算采用的为ANSYSFluent(15.0)软件.ANSYSFluent软件是目前国际上比较流行的商用CFD软件包，凡是和流体、热传递和化学反应等有关的业均可使用。ANSYSFlu前软件针对每一种流动的物理问题特点，采用适合于它的数值解法，可以在计算速度、稳定性和精度等各方面达到最佳，且作为一款国际知名的CFD仿真工具，得到r广泛的应用及验证。 进口模型管道的进口

9、给定流速进口，污水速度大小为0.1m/s、0.15m/s、0.3ni/s,0.4m/s.湍流强度为6%。模型喷口进口同样给定流速进口，速度大小0.05ni/s.0.3m/s.湍流强度6%。 壁面计算域上表面设置为对称型式,表示不存在垂忙于该表面的物理量:;为区分扩散器内外部.将模型内部设置为内部区域;其他各界面作为不流通界面设计为壁面。 出口在计算域上部每个上升管左右各开设一个出口，以保证整个计算域的质量守恒及流动均衡；出口设置为压力出口O(5) 流体届性。本文在计算过程中考虑管道内部的污水与海水，不同流体属性设置如表1所示.管道外部海域流体在本文中只考虑重力因素，设置垂向流动°本文

10、中暂不考虑管道外部海域流体的流速和波浪引起的横向流动。类别污水粘度：().001Pa-H密度：998kg/m材料选择海水粘度:().(M)IPa-s密度：1()21kg/m'空T粘度：l.7894r-5Pa-s密度：1.225kg/m'操作压力101325Pa-9.81m/s2表I材料及相关属性Tab.1Vlatfiialsandndahuipn>|M*rti<*s+Jp(u-Mx),n(L4=JSdV(1)(6) i|*算模型。本文模型设定基本方式如下：控制方程:N-S方程；湍流模型:4-NG模型；多相流模型：VOF(模型中包括污水及海水两种流体)三大守恒方程的积

11、分形式： 质虽守恒方程式中：p为流体密度；u和外分别为流体和网格运动速度；，为时间坐标；V和.4分别为网格单元体积和位元面的面积；n为单元面的外法线方向；S-为质枇源项。 动ht守恒方程JpiidV+Jpiz(u-ndA=|/xH,Vw<M+js.ndl+JSdV(2)O'I.44|IS、A=-(/>+项)/+/M时(3)式中：s.b为动ht源项；&为面积力源项；/>表示广义压力；为有效粘性系数，求解层流问题时心=件.求解湍流时;心p(："C；由湍流模型确定.k和，分别为湍动能和湍流耗散率；/为单位张挝。£能址守恒方程*(pHdV+-iin

12、fLA=("V/爪14+岑+JS“爪14+J.Sm(IV(4)S,lt=tu(5)S"、=Sd”+S,(6)式中：rw为扩散项系数."=(§+斜；11为流体总海A为导热系数；c为比热容；Pr,为滴流普朗特数;方为血枳力做功部分；S血为嬉变项；S,为化学反应热源项3计算结果分析计算采用非稳态流动.时间步K0.01-0.05S.采用()卜'多相流模型模拟污水在管路及喷嘴中的填充过程出中进II设定为速度进口边界条件，进口管长与进门管径的比位大丁T0以保障避【I速度为充分发展段.入口污水百分比为100%,计算时间不低于360$。计算过程中对模型的残差进行

13、整体监控，结合迭代过程库朗数的实时变化以判断模型的整体收敛。污水aK；甘“呵髯孩平图2流速0.1m/s时刻320时排污管及周1时海域污水与海水混合状况Kig.2(<>m|M)nrntdistributionof<lifFiisi(>nandstirnnindings<aaarea(0.1ni/s.32()s)流速设置根据以往程经验，分别为0.1m/sj).15m/s、0.3m/s和0.4m/so图2中给出了在不同污水流速条件下，扩散器内部污水的流动过程。(1) 污水流速0m/s计算结果。污水流速为0.1i)i/s的计算云图如图2所示°如图2所知，入口速度

14、为().1m/s时,排污管内污水皤终止于第3个支管，扩散器内部的海水不能全部被清除，该流速情况卜的扩散器海水清除能Jj较弱，污水流速0.15m/s计算结果。污水流速为0.15m/s的计算云图如图3所示°根据汁算结果可知.该流速下污水无法完成管道内海水的清除，扩散器模型中末端I.升管位置的海水仍然*在，即该流速低于深水排污的临界清除流速，但与污水流速为0.1m/s情况相比，在同样时刻下，扩散器内部的海水清除速度相对较快，污水从扩散器喷出后，能与周边海水较快混合因此在污水流速小于海水清除速度的条件下，污水流速越高.海水清除能力及污水稀祥效果越好。(2) 污水流速0.3m/s计算结果。污水

15、流速为0.3m/s的计算云图如图4所示将污水的流速升至0.3m/s,对比管道内污水的扩散进程可发现.ft100s时刻既完成了管道内海水的全部清除，但海水的分层现象比较明显，即在管道内部污水明显存在于管道上部，而海水存在于管道下部，由此可知.该流速高于扩散器内部海水的临界清除流速选取的概化扩散器模型的海水临界清除流速介于0.15-0.3ni/*由于概化模型参考了目前大部分污水排海扩散器的实际尺寸，因此分析结果可作为污水排海实际匚程中海水清除的参号依据同时为r分析不同污水流速情况下的海水清除规律，本文乂选取了较大的污水流速进行了分析”(3) 污水流速0.4mA计算结果。污水流速为0.4m/s的计算

16、云图如图5所示。如图5所知，当污水流速达到0.4m/s时,扩散器中的海水能够在更短的时间内被清除,约在90s时.完成海水清除过程，海水在扩散器内部出现分层的现象也不明显.说明在污水流速大于临界清除流速时,0.4m/s较0.3m/s的污水流速海水清除效果更佳，因此在程条件允许的情况下.若进行扩散器海水清除丁.作，可号虑使用较大的污水流速，4结论图3流速0.15m/s时刻320s时排污肯及周围海域污水与海水混合状况Fig.3(<>ni|M>ncnt(lislril)uti(mof(liflusi<»naiulsurnninding>raan*a(0.15ni

17、/s.32()s)图4流速0.3ni/s时刻100s时排污管及周围海域污水与海水混介状况Fig.4(j)in|M»ncnl<lislrilHi!ionofdiffusionandsiirnHindingan*a(0.3m/s.UK)s)图5流速0.4i心时刻9()-时排污管及周围海域污水'j海水混合状况Fig.5Coin|M>ncnldistrilnitioiiofdiflusionandsurrouixiingmianwa(0.4m/s,90s)(1) 本文采用数值模拟的方法,计算r不同流速况下扩散器内污水对海水的清除过程从物性上来看,两名的粘度处相同数量级，且

18、为小富度流体顶替大密度流体.参与乡相流理论，卵替的效果将会是污水相尖端指进，同时悬浮于海水相的上缘.实际的仿真计算中同样表现出该特征(2) 在低流速工况下,入口速度大的清除效率更高入【I速度为0.1m/s时,排污代内污水终止第3个支管.当入I速度升至0.15m/s时，排污管内污水终止第二个交管当污水通过喷II流到海域中时#在因密度差而产生的对流浮力项以及流体悴遍的扩散项，两项对比凹发现对流项KI*：导地位.所4污水均集中在扩散器上部，排污管以下水体,基本不受影响。在高流速况R同样表现出入II速度大海水清除效率高的规律，当入口速度为0.3m/s时，在100s时刻完成排污管内海水的全部驱替，当入口

19、速度升至0.4伪时.在90s时刻即完成了排污管内海水的驱替。此状态下污水在海水中的扩散与低的速况卜*相同规律(3)海水入侵是污水排海工程的主要技术难题之一,对扩散器的危害较大目前对于扩散器的研究主要集中于其结构参数的研究，对后期运行过程中的问题关注力度不足。本文采用数值模拟的手段,对于扩散器在不同污水流速条件下的海水清除过程进行了研究，分析了在污水流速低于和高于临界清除流速条件下的海水清除规律，为解决实际工程问题提供了理论依据。但对于更大污水流速情况下的海水清除过程还有待进一步研究,结合已有的研究成果，可得出较为准确的扩散器海水清除机制，具有较大的研究意义。参考文献：1于航，白景峰,王心海.大

20、水深条件下污水排放工程扩散器稀释扩散效果研究J.水道港口，2013,34(4):364-368.YUH,BAIJF,WANGXH.ResearchonsewagemarinedischargediffusioneffectunderdeepwaterJ.JournalofWaterwayandHarbor,2013,34(4):364-368.2jBurrowsR,MortRB.Waveactiononmulti-risermarineoutfallsJ.TheFluidEngineeringCenter,Cranfield,1989:367-380.3jDoyleBM,MackinnonPA

21、,HamillGA.DevelopmentofanumericalmodeltosimulatewaveinducedflowpatternsinalongseaoutfalllJ.AdvancesinFluidMechanics,2000:33-42.4 陈丽星，姚朝晖.Renolds应力模型在热分层流场的应用及优化j.清华大学学报：自然科学版,2003,43(8):19-36.CHENLX,YAOZH.Reynoldsstressmodelof11equationsinthethermallystratiGedturbulentflowsj.J.TsinghuaUniv.:Sci.&

22、;Tech.,2003,43(8):19-36.5 吴玮.污水海洋扩散器海水入侵及清除特性研究D.南京:河海大学,2006.ResearchonnumericalsimulationofdiffuserseawatercleanindifferentflowvelocityofsewageonWaterTransportEngineering,MinistryofTransport,Tianjin300456,China)Abstract：Insewagedischargeproject,thediffusercanprovidegreatcapabilityofinitialdilution,anditcanpreferablyprotectmarineecologyandwaterquality.Aimingatseawaterinvasion,thenumericalsimulationwasusedasthemai