内河挖入式港池回流流速分布规律的数值模拟.docx

1、水道港dJournalofWaterwayandHarbor内河挖入式港池回流流速分布规律的数值模拟戴勇1,王定略2（1.江苏省交通规划设计院股份有限公司，南京210005；2.河海大学港口海岸与近海工程学院，南京210098）摘要：采用k-£双方程湍流数学模型，建立了三维内河挖入式港池数学模型，分析r不同主流流速及不同港池轴线与主流流向的夹角对内河挖入式港池内的水流运动特性的影响，获得了港池内回流流速的分布规律。研究结果表明，垂向流速值随着主流流速的增大也相应增大;回流横轴断面上各测线平均流速从回流中心向两侧逐渐增大-数值计算结果与模型试验的对比分析表明,数值模拟能够较好的模拟内河

2、挖入式港池内的回流运动特性。关键词:挖入式港池;数值模拟；回流;流速分布规律中图分类号:TV142;0242.1文献标识码：A文章编号:1005-8443（2012）04-0299-04挖入式港池具有充分利用岸线，不碍航和不减小河道行洪断面，港池掩护条件好等优点,是内河中应用较广泛的码头平面布置型式。但是，内河挖入式港池由于河宽突然扩大，流速减小，在港池内产生回流运动，从而导致水流中携带的大量泥沙在港池内落淤,长期以来内河（特别是多沙河段）挖入式港池的水流泥沙特性是工程界十分关注的问题2）,因此研究挖入式港池的水流运动特性对于认识挖入式港池的淤积特性，优化港池的设计布置,从而达到防淤、减淤具有

3、重要的现实意义。日前，内河挖入式港池内的水流运动特性研究还主要依靠模型试验的方法,随着计算机技术的日趋成熟和计算机性能的不断提高，数值模拟花费少、速度快、信息完整等优点逐步显示出来，并逐渐被工程界所接受，成为物理模型试验的重要补充手段。对于内河挖入式港池的水流特性数值研究，国内外学者已经做了不少。如王家会等在三维水流数学模型中对自由水平面位置、垂线平均流速的确定及各方程迭代求解时的不同修正等方程进行探讨。周旦应用数学模型对回流区的表面流速、回流范围进行计算。本文综合考虑各方面的因素，结合物理模型试验同步对内河挖入式港池进行了三维流场的数值模拟研究，以全面了解内河挖入式港池的回流流速的分布规律,

4、为挖入式港池的防淤、减淤提供依据。1物理模型试验长沙理工大学水利工程学院在宽1.5m、长30m的矩形水槽内对挖入式港池内的回流流速规律进行了试验研究。港池布置在水槽的中部，试验流量由上游的三角薄壁堰测定，表面流场采用基于粒子图像测速技术的流场实时测量系统（VDMS）量测。试验布置见图102数学模型2.1网格的生成在挖入式港池计算区域网格的生成过程中，网格的布置充分考虑水力要素的变化:水力要素变化大的区域网格密集，变化平缓的区域网格疏松（图2,图3）。从图2和图3中可以看出，在自由面附近,网格在垂向上收稿日期：2011-10-26;修回日期：2011-12-12作者简介:戴勇（1979-）,男，

5、江苏省盐城人,工程师，主要从事港口及航道工程设计工作。Biography:DAIYong（1979-）,male,engineer.回流横轴中心（侧线离横轴中心距离）1bAA剖面图测点布置图1挖入式港池试验布置示意图Fig.lExperimentalarrangementofdredgedharborbasins图2立面网格分布图Fig.2VerticaldistributionofgridL土河道挖入式港池图3平面网格分布图Fig.3Planedistributionofgrid的密度比其他部分大很多，这是因为水质点的运动速度只有在自由表面附近是显著的，随着水深的加大，水质点的运动速度呈指数

6、衰减，加密网格是为了很好地捕捉自由面附近水流的变化情况,而在其他部分将网格合理的稀疏是为了减少不必要的计算负担。2.2控制方程本文采用通用性能好且工程中普遍应用的A-#两方程模型模化湍流，控制方程如下连续方程坐+略)二0dtdxi(1)动量方程岑是网）噫汨S*制随(2)k为程些旦注（伊A）=_L|（+色）票1地亟（业+业）*dtdXidXi1dxiJdXjdx,g(3)S方程MpJMpu£)=_±_|(时&)些+坐住业(业+业)-*dtdxidx,1eregJkdXjdXjgk(4)式中:i为时间；P和分别为密度和分子粘性系数心和X.分别为速度分量和坐标分量;P为修

7、正的压力；湍动粘度内可以表示成A和£的函数，即(5)(5)k2l=P(MEdxkdxk式中*为时间为湍动耗散率，被定义为EF察攀。S和E分别为与湍动能K和耗散率，对2oxkdXk应的Prandtl数，根据Launder等的推荐值及后来的实验验证，模型常数一般取=0.09=1.44=1.92.=1.0、=1.3。2.3离散方案及边界条件在计算域中采用有限体积法进行控制方程的离散。速度压力耦合采用SIMPLEC算法。进口采用速度边界条件，由流量以断面平均流速给出;出口采用给定恒定水位;自由表面的追踪则采用体积率法(VOF),VOF是用体积率函数来表示各相所占的体积，简单易行，对于具有明显

8、分界面的两相流动，有很好的模拟效果。近年来在水利工程中应用颇多。3计算结果分析本文建立与物理模型试验同样尺寸的数值模型，通过改变主流流速V、港池轴线与主流流向的夹角。计算了4种工况，最后所得回流横轴流速横向分布与实测值的比较如图4。计算值与试验值吻合较好。能够准确模拟出垂线流速沿横向分布规律，垂向流速值从中间向两侧逐渐增大,边壁处比中心的流速变化幅度大;垂向流速值随着主流流速的增大也相应增大。只是在主流流速小的时候,计算值与试验值稍有差别，这主要是因为模型网格以及流体的物理特性的影响。4-b工况二-254-c工况三4-d工况四图4垂线平均流速计算值与试验值对比图Fig.4Comparisonb

9、etweennumericalresultsandmodeltestdataofaverageverticalvelocity4结论与展望(1) 采用VOF法和标准K-E湍流模型相结合的方法对内河挖入式港池回流流速分布规律的模拟是比较成功的。计算所得回流横轴流速的横向分布与物理模型试验数据吻合较好。(2) 回流横轴断面上各测线平均流速从回流中心向两侧逐渐增大。(3) 鉴于数值模型具有省时、省力、节省费用以及修改方便的优点，在今后的研究中可以开展挖入式港池内水流运动特性对泥沙淤积的影响。参考文献：:I.许英.刘国龙.黄贵弥.挖入式港池泥沙淤枳试验研究JL江苏科技大学学报.2007.21(6)：1

10、-3.XUY.L1UGL,HUANGGB.Theexcavalrd-inharlxirImsinsedimentationtestiJ.JournalofJiangsuUniversityofScienceandTechnolog),2007,21(6):1-3.2 佛时琳.沈小雄.慈庆玲.内河挖入式港池减淤措施的试验研究J.水运工程.2005(12):14-18.HASL.SHENXX,CIQI.Inlandexcavaltul-inbasindepositionreductiontest,J.Walertransportationengineering,2(X)5(12):14-18.3

11、王家会.李岩.袁月平.等.挖入式港池水流特性数值模拟:J.人民此河.2009(6):36-37.WANGJII.1.1Y,YUANYP,etal.Excavated-inbasinflowcharacteristicsnumericalsimulationIJTheYellowRiver,2009(6)：36-37.4 周旦.内河挖入式港池水流运动特性研究Dj.K沙:长沙理工大学.2008.5 周旦.沈小雄.李建习.内河挖入式港池回流流速分布规律初探J.水道港n.2007.28(6):415-417.ZHOUD.SHENXX.LIJX.Inlandexcavated-inbasinflowve

12、locitydistribution.J.JounialofWaterwayandHarbor,2007,28(6):415-417.沈小雄.施时琳.刘虎英.内河挖入式港池间流范围的试羚研究J.长沙交通学院学报.2003.19(6):49-53.SHENXX,IIASL.LIUHY.Cin*ulatingflowinrange-testresearch、j.ChangshaCommunicationsInstitute.2003,19(6):49-53.7吴持恭.水力学M.北京:人民教行出版社，1979.：8PrandtlL流体力学概论M匕京:科学出版It.1984.Numericalsimu

13、lationofcounterflowvelocitydistributionruleindredgedharborbasinDAIYong1,WANGDing-lue2(1Jituifpiii/ProvincialCommunictdionsHanning(uulDesignInstituteCo.Ltd.Nanjing210005,China;l.ColleofHarbor,CoastalandOffshoreEngineering,HohaiUniversityNanjing21()098,China)Abstract：Usingk-etwoequationmodel,the3-1)nu

14、mericalmodelofdredgedharborbasinswasestablishedtoanalyzethecharacteristicsofflowfieldatIliedredgedharborbasinsunderdifferentverticalvelocitiesanddifferentangleslx»twrenbasinaxisandmainstreamflow.Thecounterflowvelocitydistributionwasobtained.TheresultsshowthalIheverticalvelocityincreaseswithacorrespondingmainstreamvelocity,andaveragevelocityincreasesgraduallyfromthecentertobothsidesofthebasincross-section.Comparisonbetweenthenumericalresultsandmodeltestdatashowsasatis