Chinese Journal of Ship Research
KCS标称伴流场的尺度效应数值分析总体与性能
张海鹏，张东汗，郭春雨，王恋舟，刘恬150001 哈尔滨工程大学 船舶工程学院，黑龙江 哈尔滨
KCS摘 要：为了研究标称伴流场的尺度效应，对不考虑自由液面效应的 船的粘性绕流场进行研究，并基于RANS 7方法和 SST k - ω模型对包含实尺度的 种不同尺度下的标称伴流场进行数值计算。然后，将模拟结果与试验数据进行对比，进一步分析标称伴流场的尺度效应。结果显示：各半径处平均轴向伴流分数的倒数与雷诺数的对数呈正相关；KCS 2裸船体桨盘面处的标称伴流场存在 个伴流峰，且伴流峰值会随雷诺数的增加而减小，有利于减小螺旋桨的空泡和激振力；小尺度模型的尺度效应更为明显，且内半径处的平均轴向伴流分数尺度效应问题比外半径处的更为严重。
0引言
船舶在水中航行时会产生伴流，船舶伴流场 是个十分复杂的三维速度场。螺旋桨在不均匀伴流场中工作，船、桨之间存在着相互作用，会形成力［1］。空泡并会产生强烈的脉动激振 为了改善船
舶航行性能，提高螺旋桨推进效率，并最终实现减振降噪的目的，需要准确预报和分析船舶伴流场。拖曳水池中的模型试验是在傅汝德数相同，而雷诺数不一致的条件下进行的，可以认为模型的伴流分数与实船相等。在拖曳水池中，船模试验的雷诺数范围仅为106 ~107 ，但实尺度船舶在海面航行时，其雷诺数一般能达到109 。雷诺数的增加会削减速度边界层的厚度，从而改变近壁区速度边界层的分布，这种因为雷诺数不同而导致实尺度与模型尺度之间伴流场差异的现象，被称之为“伴流场的尺度效应”。过去，由于船型更瘦更短，因此由伴流场的尺度效应所带来的问题并不突出。但随着大方形系数船舶和超级油轮的诞生，伴流场的尺度效应使得模型与实船之间的空化、激振力、噪声以及推进性能产生了显著差异［2］。为了适应现代造船业的发展，满足日益增长的未来海运、海战的需要，缩短与国外先进水平的差距，国内众多学者针对尺度效应展开了研究。黄家彬等［3］通过计算某集装箱船在不同尺度下的尾流场，讨论了尺度效应对艉部平均轴向速度场的影响，试验中采用数值计算求解船艉速度和伴流场，对不同缩尺比模型的轴向伴流、伴流峰宽度和伴流等值图等参数进行了比较，对尺度伴流修正意义重大。王展智等［4］对不考虑自由液面效应4 000 TEU DTMB 5415的 以及 的标称伴流场进行RANS研究，采用 方法对多种尺度下的粘性绕流场进行数值计算，详细研究标称伴流场的尺度效应，并重点探讨标称伴流场的各特性与雷诺数之间的定量关系，提出了一种将模型伴流场换算至实尺度的外插方法。郭春雨等［ 5 ］利用软件STAR-CCM+对不同尺度的KCS船模进行伴流场KCS模拟，同时基于“智能假体”的概念对 船艉进行收缩变形处理，分析了船艉变形与伴流场尺度效应之间的关系，试验结果对于研究标称伴流的KCS尺度效应有着重要意义。张恒等［6］以 船为研CFD究对象，利用 软件对不同雷诺数下的船艉流场分布进行了分析，通过对比计算值与试验值发现，两者阻力系数相差不大，船艉速度场也基本吻合。KCS本文将对不考虑自由液面效应的 船的粘RANS性绕流场进行研究，采用 方法和 SST k -ω模型对包含实尺度在内的不同尺度下的粘性绕流场进行数值求解，并在哈尔滨工程大学的船模水池实验室对缩尺比 λ = 52.667 的模型进行基于PIV的标称伴流场测量试验：首先，将数值计算结果与试验数据进行对比，用以验证本文所采用数 值计算方法的正确性；然后，进一步对标称伴流场的尺度效应进行详细分析。
1 数学模型 1.1 计算模型及主要参数
KCS KRISO 研究所（Korea Research是由韩国Institute of Ships and Ocean Engineering ）设计建3 600 TEU造的 型集装箱船，在以往各届国际水动力研讨会议上，该船都被选作基准船型。可见，该船具备足够的试验数据和大量各单位提供的数值计算结果，能为本文的校核计算结果提供参考数据。KCS 1所示，KCS船舶模型如图 船的主尺度可参见文献［7-8］。 为了对标称伴流场的尺度效应进行数值研KCS究，本文以 集装箱船为研究对象，对一系列不考虑自由液面、升沉和纵倾的缩尺比模型进行了0.28，运数值求解。对应的傅汝德数为 动粘性系数 ν = 1.139 0 ´ 10-6 m2/s 。由于船模的主尺度不同，因此划分网格时所使用的参数及网格数量也1所存在着差异，各缩尺比模型的计算状态如表示。表中：λ 为缩尺比；L 为垂线间长；V 为船pp 1 +舶航速；Re 为雷诺数；y 为近壁面第 层网格无1量纲距离；Dy 为第 层边界层厚度。
1.2 网格划分及湍流模型
Surface首先，对船舶模型进行面网格重构（ remesher），对船舶表面网格进行处理并生成表面三角化良好的高质量面网格，然后以面网格为基础生成带有边界层（Prism layer mesh）和切割体网格（Trimmed mesh 2 ）的体网格。网格分布如图所示。本文使用的湍流模型为SST k - ω模型，其全Shear Stress Transport）k称为剪切应力输运（ -ω Wilox［9模型。在近壁区，使用经 ］修正过的 k -ω模型，边界层以外的自由剪切流则使用 k - ε 模型，2种模型间的过渡应用一种混合函数完成，适用于不可压缩状态及可压缩状态。有关SST k -ω湍流模型详见文献［10-11］。
2 计算结果与分析 2.1 桨盘面平均轴向伴流分数尺度效应分析
当缩尺比 λ = 52.667 时，桨盘面处无量纲轴向3速度的计算值与试验值的比较如图 所示，其中左侧为试验数据，右侧为数值模拟数据。图中， Vx /V0 为伴流分数 w(w = 1 - Vx /V 0)的分布情况。3由图 可以看出，数值模拟结果与试验数据好［12-13］，但符合较 由数值模拟所得到的伴流场比PIV试验测量所得到的更加光顺、均匀。这说明本文所采用的网格划分方法和数值计算方法适合KCS研究 集装箱船的标称伴流场，且在数值模拟中不用考虑自由液面的影响［14-15］。
2所不同缩尺比下的平均轴向伴流分数如表示。平均轴向伴流分数随雷诺数对数变化的关系4如图 所示。平均轴向伴流分数的倒数随雷诺数5对数变化的关系如图 所示。4由图 可以看出，各半径处平均轴向伴流分数随雷诺数对数的变化趋势较为一致。低雷诺数时，平均轴向伴流分数随雷诺数的变化较大，尺度
效应明显；但随着雷诺数的变大，平均轴向伴流分数随雷诺数对数的变化趋于平缓，尺度效应问题开始变得不突出。当雷诺数大于109 时，桨盘面各半径处的平均轴向伴流分数随雷诺数几乎不再变化；另外，随着半径的增大，平均轴向伴流分数迅速减小，平均轴向伴流分数随雷诺数的变化趋于平缓，表明桨盘面处内半径处的平均轴向伴流5分数尺度效应问题比外半径处的更为严重。由图可以发现，在各个半径处，平均轴向伴流分数的倒 数与雷诺数的对数在中、外半径处呈近似分段线性关系。从斜率上看，平均轴向伴流分数的倒数3 3随雷诺数对数的变化关系分为个阶段，这 段的斜率是依次减小的，表明尺度效应问题也经历着3从“严重”到趋于“缓和”的 个阶段。
2.2 轴向标称伴流场的尺度效应分析
本节将进一步分析轴向标称伴流场随雷诺数的变化。船体绕流流动的重点考察对象为船体尾部边界层和近尾流，由于桨盘面处的速度场是螺旋桨的进流，所以研究桨盘面处的伴流场就显得尤为重要。不同缩尺比下桨盘面处的无量纲轴向6速度云图如图 所示（图中实线为螺旋桨半径，虚1.1 0.6线则分别为 倍和 倍的螺旋桨半径）。
6由图 可以看出：从模型尺度到实尺度，随着雷诺数的增加，速度边界层的厚度逐渐减小，舭部附近的边界层最厚；轴向速度等值线逐渐向螺旋桨盘面中心收缩，部分轴向速度等值线经历了向中心收缩、断裂、甚至消失的过程，这种收缩效应Sasajima Schuiling等［16］和 等［17］关与 于伴流换算方法的基本思想一致；随着雷诺数的增加，伴流分数变化的梯度有所降低，桨盘面处的伴流场将变得更加均匀。可见，直接应用船模尺度的标称伴流场作为螺旋桨适伴流设计的依据，所预报结果的精度并不十分可靠，必须要考虑尺度效应对其的影响。不同缩尺比下桨盘面处轴向伴流分数的周向7所示。图中，0°指桨盘面正上方位置，分布如图而正方向是从船艉向船艏看的顺时针方向。为了更加清晰地显示伴流的形式，还对伴流分布进行了额外的一个周期延拓。7 KCS由图 可以看出， 船桨盘面处轴向伴流分数的周向分布总体呈“W”型。尾流场在一个周2 0°处，另一个期内有 个伴流峰，一个是角位置为180°处。随着雷诺数的增加，伴流峰是角位置为值逐渐减小，这对于大尺度的船模，对空泡和激振力都是有利的。在内半径区域 (0.3 r/R 0.5) ， 180°的伴流峰值处，会出现“V”型分在角位置为V布，这种“”型分布在 r/R = 0.3 处最为明显，在处彻底消失。“V”型分布与雷诺数的大小r/R = 0.6也有关系，在“V”型分布最为明显的 r/R = 0.3 处， V =52.667，31.599，出现“”型分布的缩尺比为 λ 16，8，当 =4时“V”型分布已不明显，此时雷诺数λ =4 1应不超过 λ 时的雷诺数，由表 可得相对应的
图7不同缩尺比下桨盘面轴向伴流分数的周向分布Fig.7 Circumferential distribution of axial wake fraction on propeller at different scale ratios
处，“V”型分雷诺数 Re 3.1 ´ 108 ；而在 r/R = 0.6 =52.667，由表 1布已不明显，此时对应的缩尺比 λ可得相对应的雷诺数 Re 6.5 ´ 106 。
3结语
RANS本文采用 方法，结合 SST k - ω 湍流模KCS型研究了 标称伴流场的尺度效应，得到以下主要结论： 1）在各半径处，平均轴向伴流分数的倒数与雷诺数的对数呈正相关，和大尺度模型相比，小尺度模型的尺度效应更为明显，且在内半径处平均轴向伴流分数的尺度效应问题比外半径处的更为严重。2）随着雷诺数的增加，边界层的厚度逐渐减小，轴向速度等值线逐渐向桨盘面中心收缩，伴流分数变化的梯度降低，伴流场将变得更加均匀。3）桨盘面处标称伴流场的2个伴流峰分别处于正上方和正下方位置，伴流峰值会随雷诺数的增加而减小，这将有利于螺旋桨的空泡和激振力。4）桨盘面处轴向伴流分数的周向分布总体呈“W”型 180°位置会出现“V”型，在 分布，但随着半径和雷诺数的增大，这种“V”型分布会更加不明显，直至消失。
参考文献：
［1］ LEE J H LEE K J ，PARK H G，et al. Possibility of ， air-filled rubber membrane for reducing hull exciting pressure induced by propeller cavitation［J］. Ocean Engineering，2015，103：160-170. 2］许晶，周连第，高秋新. ［ 标称伴流尺度效应的数值模拟［J］. 水动力学研究与进展，1998，13（4）：481-490. XUJ ，ZHOU L D，GAO Q X. Numerical simulation of scale effect on nominal wake ［J］. Journal of Hydrodynamics，1998，13（4）：481-490（in Chinese）. 3］黄家彬，陈霞萍，朱仁传，等. CFD ［ 基于 的标称伴流C］//第九届全国水动力学学术会场尺度效应研究［议暨第二十二届全国水动力学研讨会论文集.北京：中国力学学会，2009. 4］王展智，熊鹰，黄政，等. ［ 双桨船轴向伴流场尺度效应J］. 2015，49（4）：的数值研究［ 上海交通大学学报， 457-463. WANG Z Z，XIONG Y，HUANG Z，et al. Numerical study on scale effect of axial wake of twin screw ship ［J］. Journal of Shanghai Jiaotong University，2015，49 （4）：457-463（in Chinese）. 5］郭春雨，张琪，陈鸽，等. SDM ［ 基于 方法的船艉伴流场尺度效应研究与修正［J］.中国舰船研究，2015，10 （6）：1-7. 7 GUO C Y ，ZHANG Q，CHEN G，et al. Research and revision on the scale effect of the wake field based on the method of SDM ［J］. Chinese Journal Ship Research，2015，10（6）：1-7（in Chinese）. 6］张恒，詹成胜，刘祖源，等. ［ 基于船舶阻力性能的船型主尺度参数敏感度分析［J］.船舶工程，2015，37（6）： 11-14. ZHANG H，ZHAN C S， LIU Z Y ，et al. Analysis of sensitivity of ship type main dimension parameters based on ship resistance performance ［J］. Ship Engineering，2015，37（6）：11-14（in Chinese）. ［7］ SHEN Z R， WAN D C ，CARRICA P M. Dynamic overset grids in OpenFOAM with application to KCS self-propulsion and maneuvering ［J］. Ocean Engineering，2015，108：287-306. ［8］ FUJISAWA J，UKON Y，KUME K，et al. Local velocity field measurements around the KCS model（SRI M.S. No.631） in the SRI 400 m towing tank： ship performance division report ：00-003-02［R］.［S.l.］： Ship Research Institute，2000. ［9］ WILOX D C. Turbulence modeling for CFD［M］. La Canada，California：DCW Industries，Inc.，1998. 10 SHIH T H，LIOU W W，SHABBIR A，et al. A new ［ ］ eddy viscosity model for high Reynolds number k-ε turbulent flows［J］. Computers & Fluids，1995，24 （3）：227-238. ［11］ MENTER F R. Two-equation eddy-viscosity turbulence models for engineering applications［J］. AIAA Journal，1994，32（8）：1598-1605. ［12］ KIM W J VAN S H ，KIM D H. Measurement of flows ， around modern commercial ship models ［J］. Experiments in Fluids，2001，31（5）：567-578. ［13］ YANG H U， KIM B N YOO J H ，et al. Wake ， comparison between model and full scale ships using CFD［J］. Journal of the Society of Naval Architects of Korea，2010，47（2）：150-162. ［14］ WANG Z Z，XIONG Y，WANG R，et al. Numerical study on scale effect of nominal wake of single screw ship［J］. Ocean Engineering，2015，104：437-451. ［15］ KIM D，LEE K，SEONG W. Non-cavitating propeller noise modeling and inversion［J］. Journal of Sound and Vibration，2014，333（24）：6424-6437. ［16］ SASAJIMA H，TANAKA I. On the estimation of wakes of ships［C］//Proceedings of 11th ITTC. Tokyo，Japan： ITTC，1966：140-144. ［17］ SCHUILING B，LAFEBER F H，VAN DER PLOEG A，et al. The influence of the wake scale effect on the prediction of hull pressures due to cavitating propellers ［C］//Proceedings of the Second International Symposium on Marine Propulsors. Hamburg，Germany：［s.n.］，2011.