第１４卷第８期 ２０１０年８月 文章编号：１００７—７２９４（２Ｏ１Ｏ）０８—０８５４—０９ 船舶力学 Ｊｏｕｒｎａｌ ｏｆ Ｓｈｉｐ Ｍｅｃｈａｎｉｃｓ Ｖ０１．１４ Ｎｏ．８ Ａｕｇ．２０１０ 大侧斜桨敞水性能及流场数值研究 杨仁友，沈泓萃，姚惠之 （中国船舶科学研究中心，江苏无锡２１４０８２） 摘要：针对库存大侧斜桨，基于商用ＲＡＮＳ代码，利用结构化网格技术和流道计算模型开展了大侧斜桨敞水水动 力和流场的数值分析。在进速系数ＪＥ［Ｏ．３３，０．９５】范围内，文中计算得到的敞水螺旋桨水动力与试验结果差异 在３％以内。通过敞水螺旋桨流场计算和计及粘性影响鼓动盘计算结果比较可知，螺旋桨对桨前流场影响主要是 抽吸作用；对桨后近场区域Ｘ／Ｄ＜４流场影响主要是滑流，轴向速度周向均值沿径向分布主要受桨叶负荷分布控 制；到中间区域４＜Ｘ／Ｄ ＜１５，旋转尾流加速了轴向速度沿纵向的衰减，但增大了沿径向的扩散范围；至充分发展 尾流区Ｘ／Ｄ ＞１５，螺旋桨尾流变成了单一的剪切流。 ．关键词：大侧斜桨；敞水性能；尾流场；数值模拟 中图分类号：Ｕ６６１．３ 文献标识码：Ａ Ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ ｏｎ ｔｈｅ ｏｐｅｎ－ｗａｔｅｒ ｐｅｒｆｏｒｍａｎｃｅ ａｎｄ ｆｌｏｗ ｆｉｅｌｄｓ ｏｆ ｔｈｅ ｈｉｇｈ－ｓｋｅｗ ｐｒｏｐｅｌｌｅｒ ＹＡ ＮＧ Ｒｅｎ－ｙｏｕ，ＳＨＥＮ Ｈｏｎｇ－ｃｕｉ，ＹＡ ０ Ｈｕｉ－ｚｈｉ （Ｃｈｉｎａ Ｓｈｉｐ Ｓｃｉｅｎｔｉｉｆｃ Ｒｅｓｅａｒｃｈ Ｃｅｎｔｅｒ，Ｗｕｘｉ ２１４０８２，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｕｓｉｎｇ ｔｈｅ ｓｔｒｕｃｔｕｒｅｄ－ｇｒｉｄ ｔｅｃｈｎｏｌｏｇｙ ａｎｄ ｔｈｅ ｆｌｏｗ ｍｏｄｅｌ ａｒｏｕｎｄ ｏｎｅ ｂｌａｄｅ ｐａｓｓａｇｅ，ｔｈｅ ｈｙｄｒｏ。 ｄｙｎａｍｉｃ ｐｅｒｆｏｒｍａｎｃｅ ａｎｄ ｗａｋｅ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｔｈｅ ｈｉｇｈ－ｓｋｅｗ ｓｔｏｃｋ ｐｒｏｐｅｌｌｅｒ ｗｅｒｅ ｎｕｍｅｒｉｃａｌｌｙ ａｎａｌｙｚｅｄ ｗｉｔｈ ｃｏｍｍｅｒｃｉｌ ＲＡＮＳ ｍｅｔａｈｏｄｓ．Ａｍｏｎｇ ｔｈｅ ａｄｖａｎｃｅｄ ｒａｔｉｏ Ｊ ０．３３ ｔｏ ０．９５．ｔｈｅ ｄｉｆｆｅｒｅｎｃｅｓ ｂｅｔｗｅｅｎ ｔｈｅ ＣＦＤ ｒｅｓｕｌｔｓ ａｎｄ ｅｘｐｅｒｉｍｅｎｔｓ ａｂｏｕｔ ｔｈｅ ｏｐｅｎ－ｗａｔｅｒ ｐｅｒｆｏｒｍａｎｃｅ ａｒｅ ｗｉｔｈｉｎ ３％．Ｔｈｒｏｕｇｈ ｔｈｅ ｃｏｍｐａｒｉｓｏｎ ｏｆ ｔｈｅ ｌｆｏｗ ｆｉｅｌｄｓ ａｒｏｕｎｄ ｔｈｅ ｐｒｏｐｅｌｌｅｒ ｃｏｍｐｕｔｅｄ ｂｙ ＣＦＤ ａｎｄ ｔｈｅ ａｃｔｕａｔｏｒ—ｄｉｓｋ ｔｈｅｏｒｙ ｗｉｔｈ ｖｉｓｃｏｕｓ ｅｆｆｅｃｔ，ｉｔ ｉｓ ｉｎ— ｄｉｃａｔｅｄ ｔｈａｔ：ｔｈｅ ｅｆｆｅｃｔ ｏｆ ｐｒｏｐｅｌｌｅｒ ｓｕｃｔｉｏｎ ｏｃｃｕｒｓ ｉｎ ｆｒｏｎｔ ｏｆ ｔｈｅ ｐｒｏｐｅｌｌｅｒ；ｉｎ ｔｈｅ ｎｅａｒ ｆｉｅｌｄ，Ｘ／Ｄ＜４，ｂｅ— ｐｈｉｎｄ ｔｈｅ ｐｒｏｐｅｌｌｅｒ，ｔｈｅ ｅｆｆｅｃｔ ｏｆ ｔｈｅ ｐｒｏｐｅｌｌｅｒ ｉｓ ｃｏｎｆｉｎｅｄ ｔｏ ｔｈｅ ｓｌｉｐｓｔｒｅａｍ，ａｎｄ ｔｈｅ ｃｉｒｃｕｍｆｅｒｅｎｔｉａｌｌｙ—ａｖｅｒ— ａｇｅｄ ａｘｉａｌ ｖｅｌｏｃｉｔｙ ｐｒｏｆｉｌｅｓ ａｌｏｎｇ ｔｈｅ ｒａｄｉａｌ ｄｉｒｅｃｔｉｏｎ ｍａｉｎｌｙ ｄｅｐｅｎｄ ｏｎ ｔｈｅ ｐｒｏｐｅｌｌｅｒ ｌｏａｄｉｎｇ ｃｏｎｄｉｔｉｏｎｓ；ｉｎ ｔｈｅ ｉｎｔｅｒｍｅｄｉａｔｅ ｒｅｇｉｏｎ，４＜ ／Ｄｐ＜１ ５，ｔｈｅ ｓｗｉｒｌｉｎｇ ｗａｋｅ ｂｒｉｎｇｓ ｔｈｅ ａｘｉａｌ ｈｉｇｈ—ｓｐｅｅｄ ｖｅｌｏｃｉｔｙ ａ ｍｏｒｅ ｑｕｉｃｋｌｙ ｄｅｃａｙ ａｌｏｎｇ ｔｈｅ ｌｏｇｎｉｔｕｄｉｎａｌ ｄｉｒｅｃｔｉｏｎ ａｎｄ ａ ｌａｒｇｅｒ ｒｅｇｉｏｎ ｉｎ ｔｈｅ ｒａｄｉａｌ ｄｉｒｅｃｔｉｏｎ；ｉｎ ｔｈｅ ｄｅｖｅｌｏｐｅｄ ｗａｋｅ， Ｘ／Ｄｐ＞１５，ｔｈｅｒｅ ｒｅｓｕｌｔｓ ｉｎ ａ ｓｉｎｇｌｅ ｓｈｅａｒ ｌｆｏｗ． Ｋｅｙ ｗｏｒｄｓ：ｈｉｇｈ－ｓｋｅｗ ｐｒｏｐｅｌｌｅｒ；ｏｐｅｎ－ｗａｔｅｒ ｐｅｒｆｏｒｍａｎｃｅ；ｆｌｏｗ ｆｉｅｌｄｓ；ｎｕｍｅｒｉｃａｌ ｓｉｍｕａｌｔｉｏｎ １引 言 快速性能作为船舶基本水动力性能之一，在船舶优化和选型设计中往往成为最基本的水动力优 化目标。因此，建立起预报精度较高的快速性数值模拟方法对于节约船舶设计成本和缩短研发周期具 收稿日期：２０１０—０５—２５ 作者简介：杨仁友（１９８０一），男，中国船舶科学研究中心博士研究生。 第８期 杨仁友等：大侧斜桨敞水性能及流场数值研究 ８５５ 有重要作用。而螺旋桨敞水性能数值模拟是船舶快速性能ＣＦＤ研究的一个重要内容，它不仅直接关 系到伴流分数能否精确预报，而且螺旋桨敝水数值模拟技术与船舶自航数值模拟技术具有共用技术。 因此开展螺旋桨敞水数值研究，对于提高船舶快速性数值预报精度具有重要作用。 目前，对于船舶螺旋桨理论研究主要分为势流和粘流计算两种方法。对于螺旋桨势流理论研究， 最早可以追溯到１９２９年哥尔斯坦在皇家学会学报上发表的螺旋桨旋涡理论，而后伴随着计算机技术 的快速发展，产生出从三维角度解决满足物面上边界条件的边值问题的升力面理论，以及面元法／升力 面耦合方法的导管桨水动力性能预报方法等。可以说，在计算机技术的推动下，势流理论计算程序作 为一种相对成熟的计算方法，已成为螺旋桨几何设计与性能预报的主要技术工具之一。 对于螺旋桨水动力的粘流计算目前主要是基于ＲＡＮＳ代码，而且通常与船体结合在一起进行研 究。起先，人们主要利用鼓动盘原理来简化螺旋桨，研究桨船之间水动力性能的相互干扰。而后，随着 ＲＡＮＳ代码的不断开发和完善，逐渐涌现出模拟真实螺旋桨几何形状的定常计算和模拟螺旋桨真实几 何和实际运转工况的非定常水动力数值计算方法。如，Ｓｔａｎｉｅｒｔ 采用ＲＡＮＳ代码，比较了三叶和五叶螺 旋桨水动力性能和尾流场的尺度效应。Ａｂｄｅｌ—Ｍａｋｓｏｕｄ和Ｒｉｅｃｋｔ２］￣ｌＪ用商业软件ＣＦＸ对带与不带螺旋 桨集装箱船模周围的流动结构进行研究。Ｔａｈａｒａ和Ａｎｄｏｔ３］利用自己开发的自航模拟求解器计算带与 不带螺旋桨的ＫＣＳ集装箱船周围的流场（考虑自由表面），并根据ＥＦＤ结果进行校核，验证了方法的 有效性。Ｓ６ｎｃｈｅｚ—Ｃａｊａｔ４－￣基于自编的ＲＡＮＳ代码，利用滑移网格技术对牵引式吊舱推进器进行数值模 拟研究，该课题隶属于欧盟０州ＰＯＤ项目。可以看出，国外对此是趋之若骛，而且成果显著。国内，主 要以中国船舶科学研究中心为代表的一些研究机构和高校开展了以螺旋桨和导管桨为对象的水动力 性能数值计算和研究工作。如张志荣］６１基于混合面模式，采用商业软件Ｆｌｕｅｎｔ对ＫＲＩＳＯ的３６００ＴＥＵ 集装箱船ＫＣＳ（带桨）的尾流场进行数值模拟，并且对螺旋桨推力和扭矩进行定量比较，对带桨船体 （不考虑自由液面）尾流场作定性分析。杨仁友等人１７１根据ＲＡＮＳ方法，采用移动网格滑移交接面技术 对假想的带前置导叶潜艇的自航试验、螺旋桨敞水试验进行模拟，并试验验证，发现该方法能够较好 地预报螺旋桨推力、扭矩。并与潜艇阻力试验相结合，预报潜艇自航因子，通过与试验结果的比较分 析，验证了方法的可靠性。对于导管桨，主要是以王涛、洪方文等人为代表开展了一系列的数值计算及 试验验证工作。王涛等嗍用ＣＦＸ—ＴＡＳＣｆｌｏｗ对轴对称体与导管推进器组合体的三维复杂粘性流场进行 了模拟；洪方文例则是基于商用Ｆｌｕｅｎｔ软件提供的ＲＡＮＳ代码对导管桨设计点位置的水动力性能进行 数值模拟，结果显示，敞水工况推力、扭矩与试验比较偏差在５％以内，艇后工况导管桨水动力与试验 最大偏差也只有６％左右。 显而易见，势流计算具有明显的计算时间优势，而粘流计算具有计算精度较高等特点。因此如何 缩短粘流计算周期，对于ＣＦＤ技术的工程应用具有重要价值，特别对于螺旋桨敞水性能曲线计算。为 了提高粘流计算效率，本文发展了基于单个敞水螺旋桨叶片绕流的流道计算模型；同时考虑到非结构 化网格难以精确控制网格节点分布，因此文中建立了桨叶叶片的结构化网格建模方法，以便于下一阶 段能够顺利开展螺旋桨敞水ＣＦＤ不确定度分析和大侧斜桨尺度效应影响研究。 ２数学模型 本文通过对商用Ｆｌｕｅｎｔ软件的二次开发，建立了基于结构化网格和单个叶片流道计算模型的大 侧斜桨敞水性能数值模拟方法。整个计算流场的控制方程表达如下： ㈣ 鲁＋－ Ｏｕ， 一 ０３ｐ…２－ 一Ｏ魄ｒｉ／一 （２） ８５６ 船舶力学 第１４卷第８期 方程组（１）、（２）是不封闭的，因此，需要寻求补充关系（湍流模型）使问题封闭。本文选取了ＳＳＴ ｋ一 ∞湍流模型供计算研究。 对于控制方程，采用有限体积法进行离散，其中对流项采用传统的二阶迎风格式，控制项采用中 心差分格式。压力、速度耦合迭代计算采用ＳＩＭＰＬＥ算法。对于螺旋桨旋转单元计算，本文分别采用了 多参考系模型（Ｍｕｌｔｉｐｌｅ Ｒｅｆｅｒｅｎｃｅ Ｆｒａｍｅ）和滑移网格方法（Ｓｌｉｄｉｎｇ—Ｍｅｓｈ Ｔｅｃｈｎｏｌｏｇｙ）。 ３几何与网格生成 本文以库存七叶大侧斜螺旋桨模型为对象计算分析。该桨模直径为０．２ｍ，０．７Ｒ处叶片螺距比达 １．２，侧斜角为４５。，右旋桨。网格生成采用多块结构化网格，其中，对于大侧斜桨附近区域，针对不同半 径区域进行分段建模，各网格块之问采用对接网格技术进行连接；导流帽和螺旋桨尾流段区域也同样 采用结构化网格，各区域网格采用交接面形式进行搭接。至于网格形式，在桨叶附近区域沿叶片弦长 方向为Ｏ型网格，沿径向方向为Ｈ型网格，在计算域其余区域沿纵向（即螺旋桨旋转轴）Ｔｙ向为Ｈ型 网格，沿周向为Ｏ型网格，如图１、２所示。 图１桨叶壁面及周围区域网格节点分布 Ｆｉｇ．１ Ｔｈｅ ｇｒｉｄ ｄｉｓｔｒｉｂｕｔｉｏｎｓ ｏｎ ｔｈｅ ｗａｌｌｓ ａｎｄ ａｄｊａｃｅｎｔ ｚｏｎｅｓ 文中通过在螺旋桨单个叶片两侧设置 旋转周期性边界条件来建立桨叶叶片的流 道计算模型。除此之外，在桨盘面前方６倍 桨径位置设置了均匀速度入口，桨盘面后方 沿着来流方向６Ｏ倍桨径位置设置了等压出 口，外场边界径向离螺旋桨旋转轴６倍桨径 图２边界条件 位置也设置了均匀速度人口，桨叶壁面及桨 Ｆｉｇ．２ Ｂｏｕｎｄｒｙ ｃｏｎｄｉｔｉｏｎｓ 毂、导流帽、顺流段均采用无滑移边界条件，见图２。 整个计算方案总的网格数达到了２４．５万。 ４计算结果与分析 针对库存七叶大侧斜桨，本文基于商用ＲＡＮＳ代码，采用结构化网格技术和单个叶片的流道计算 模型，开展大侧斜桨敞水性能曲线的计算与试验验证。并根据上述计算结果，结合计及粘性影响的鼓 动盘流场计算结果，分析敞水大侧斜桨周围流动特性。 ４．１敞水性能曲线计算 商用软件Ｆｌｕｅｎｔ对旋转运动单元主要提供了三种求解方法：多参考系模式（ｍｕｌｔｉｐｌｅ ｒｅｆｅｒｅｅｅ ｆｒａｍｅ（ＭＲＦ）ｍｏｄｅ１）、混合面模式（ｍｉｘｉｎｇ ｐｌａｎｅ ｍｏｄｅ１）和滑移网格模式（ｓｌｉｄｉｎｇ ｍｅｓｈ ｍｏｄｅｌ（ＳＬＭ））。 ＭＲＦ是这三种模式中最简单的一种，作为一种时域平均的近似稳态模式，它比较适合模拟两个子 域交接处边界近似均匀的流场。混合面模式则是通过消除由周向变化引起的非定常问题，而得到了一 个定常结果。可见，上述两种方法均是定常方法，能够有效地节约计算资源。滑移网格技术通过滑移网 ８５８ 船舶力学 第１４卷第８期 图４给出了进速系数ｊ＝ｏ．５４６工况，采用ＭＲＦ方法，。 和ＳＬＭ方法计算得到的桨叶 ＝Ｏ．７０５位置壁面压力和 剪应力分布结果，图５为这两种方法计算得到的桨叶周 ０ 围附近区域轴向无量纲速度分布结果。可以看出，ＭＲＦ和 ＳＬＭ计算得到的桨叶壁面压力、剪应力和壁面附近区域 轴向速度分布差异很小。 ｏ 综上所述，不管是从水动力积分变量，还是从螺旋桨 周围流动细节，均可看出多参考系模型（ＭＲＦ）和滑移网 格方法（ＳＬＭ）计算结果基本一致。因此，为了节约计算时０．０ 问，文中均采用ＭＲＦ方法计算大侧斜桨敞水性能曲线， ｎ。 。 ’。。 ’。 计算结果及其与试验结果比较如图６所示。由图可知，在 图６大侧斜桨敞水性能曲线计算及其与试验 进速系数ｊ＝０．３０－０．９５范围内，敞水大侧斜桨推力、扭矩 与试验结果差异均在３％以内，敞水效率除ｊ＝０．９５为结果比较 Ｆｉｇ．６ Ｔｈｅ ｃｏｍｐａｒｉｓｏｎ ｂｅｔｗｅｅｎ ＣＦＤ’ｓ ａｎｄ ｅｘ． ４．７％外，其余也在３％范围内。其中，桨叶扭矩计算结果小ｐｅｒｉｍｅｎｔａｌ ｒｅｓｕｌｔｓ ｏｆ ｔｈｅ ｏｐｅｎ—ｗａｔｅｒ ｐｅｒ－ 于试验值，敞水效率则是大于试验值。 ４．２敞水大侧斜桨尾流场分析ｆｏｒｍａｎｃｅ ｃｕｒｖｅｓ ｏｆ ｔｈ。ｈｉｇｈ－ｓｋｅｗ ｐｒｏ‘ Ｐｅｕ 螺旋桨作为一种有效的水下推进装置，在舰船领域得到了广泛的应用。因此，掌握和了解螺旋桨 周围流动特点具有重要的实用价值。由于螺旋桨引入了非定常特性、高湍流水平、梢涡、毂涡等流动， 使得船桨绕流现象十分复杂。鉴于此，本文针对设计点位置敞水螺旋桨周围流场进行数值模拟，为以 后复杂船／桨相互干扰流动的理解和分析奠定基础。 图７给出了敞水螺旋桨紧前方、桨盘面和紧后方盘面位置轴向无量纲速度云图和横向速度矢量 图。可以看出： ａ．在螺旋桨紧前方盘面，螺旋桨抽吸作用超过了整个桨盘面大小，而且沿径向向里抽吸影响（ｓｕｃ— ｔｉｏｎ ｅｆｆｅｃｔ）越显著。 ｂ．螺旋桨对紧后方盘面流场影响主要是滑流作用（ｓｌｉｐｓｔｒｅａｍ ｅｆｆｅｃｔ），它使得轴向速度显著增大， 最大值位置在ｒ／Ｒ＝Ｏ．７附近。而且受螺旋桨诱导，径向向里速度明显增大。螺旋桨产生周向旋转伴流， 最大值在桨毂附近。 （ａ）ｘ／Ｄ＝～０．２５２ ｑｏ）ｘ／Ｄ＝Ｏ．０ 【ｃ）ｘ／Ｄ ０・１２６ 图７螺旋桨周围盘面位置轴向无量纲速度云图和横向速度矢量图 Ｆｉｇ．７ Ｖｅｌｏｃｉｔｙ ｄｉｓｔｒｉｂｕｔｉｏｎｓ ｎｅａｒ ｔｈｅ ｐｒｏｐｅｌｌｅｒ’ｓ ｐｌａｔｅ：ａｘｉａｌ ｖｅｌｏｃｉｔｙ ｃｏｎｔｏｕｒｓ＆ｃｒｏｓｓ ｐｌａｎｅ ｖｅｃｔｏｒ ｉｆｅｌｄｓ 参考文献［１６］，绕桨轴作一半径为ｒ的“圆柱面Ｓ”，根据Ｓｔｏｋｅｓ定理，并假设螺旋桨尾涡不穿过该 “圆柱面Ｓ”，可得文中螺旋桨单个叶片半径ｒ剖面环量计算公式为： 州＝厂２ ＝ｒｌ』Ｖ，ｄＯ—Ｊ Ｖ，ｄＯ Ｉ＝孕 如 。） （３） 第８期 杨仁友等：大侧斜桨敞水性能及流场数值研究 ８５９ 式中， 表示周向速度沿周向的积分平均；ｆ２表示螺旋桨紧 ００７ 后方盘面半径ｒ处圆弧，文中选取盘面位置为ｘ／Ｄ＝Ｏ．１２６；Ｚ Ｏ．ｏ６ ０Ｏ５ 表示螺旋桨桨前方盘面半径ｒ处圆弧，文中选取盘面位置 为ｘ／Ｄ＝一０．２５２。 菩ｎ∞ Ｎ ０．Ｏ３ 根据（３）式计算得到螺旋桨叶片环量分布，如图８所 ０．０２ 示。由图可知： ０．Ｏ１ ａ．桨前来流除近壁半径区域外环量分布为０，因此螺 ０ｏｏ 旋桨叶片单元环量分布近似等同于螺旋桨紧后方盘面ｘ／Ｄ＝ ＿ｏ．０１ ０．１２６处环量分布。 ｂ．桨叶叶片环量沿径向向外逐渐增大，到Ｏ．７ＪＲ 达到最 图８桨叶叶片环量分布 大值，而后急剧降低，至叶梢近似为０。桨叶叶片负荷分布 Ｆｉｇ．８ Ｃｉｒｃｕｌａｔｉｏｎ ｐｒｏｆｉｌｅｓ ｏｎ ｔｈｅ ｂｌａｄｅ 也与其类似，在 ＝Ｏ．７位置达到最大。 ｃ．基于文献【１０］流场测试数据结果计算得到的敞开桨叶片环量分布结果看出，桨叶叶片梢涡会导 致该位置环量突然增大，而图８中并不明显。而从图９、１０中可以看出，本文算例确实捕捉到了梢涡对 紧后方盘面速度场的影响。 这些表明，本文能够定性捕捉到桨叶叶片梢涡，但对于定量计算，仍显不足。同时也可看出，螺旋 桨紧后方盘面位置梢涡已径向收缩至ｆｉＲ ＝０．９附近。 ０ ４５ ９０ １３５ １６０ ２２５ ２７０ ３１Ｓ ３６Ｏ （ｏ） 图９紧后方盘面ｒ／Ｒｐ＝Ｏ．９半径处，三个速度 图ｌＯ叶梢流线图 分量周向分布结果 Ｆｉｇ．１０ Ｔｈｅ ｓｔｒｅａｍｌｉｎｅ ｓｈｅｄ ｆｒｏｍ ｔｈｅ ｂｌａｄｅ ｔｉｐ Ｆｉｇ．９ Ｖｅｌｏｃｉｔｙ ｃｏｍｐｏｎｅｎｔｓ ｄｉｓｔｒｉｂｕｔｉｏｎ ａｌｏｎｇ ｔｈｅ ｃｉｒｃｕｍｆｅｒｅｎｔｉａｌ ｄｉｒｅｃｔｉｏｎ ａｔ ｔＩｌｅ ｊｕｓｔ—ｄｏｗｎ— ｓｔｒｅａｍ ｐｌａｔｅ ｒａｄｉｕｓ ｒ／Ｒ ＝０．９ 为了进一步分析螺旋桨对尾流场的影Ⅱ向，文中增加了假毂敞水计算和计及粘性影响的鼓动盘流 场计算。其中，鼓动盘面压差分布为均匀分布，盘面大小与桨径相同，整个盘面压差积分等于敞水螺旋 桨叶片产生的推力。三个方案计算得到的螺旋桨周围盘面三个速度分量周向积分平均结果比较如图 １１所示，图１２给出了桨前、桨后螺旋桨大小盘面积分平均伴流沿纵向的分布结果。根据平均流动结 果，将螺旋桨绕流分成桨前、桨后区域，桨后区域又分为近场区域、中问区域和充分发展尾流区域。基 于图中结果可知： ａ．桨前区域，螺旋桨主要影Ｈ向为抽吸作用，它使得桨前轴向、径向速度明显增大，但对切向速度影 Ｈ向甚微，结果仍为０。其中，对轴向速度，除桨毂壁面附近区域外，沿径向向里逐渐增大；而对于径向速 度，在叶梢半径处达到最大，沿径向两边逐渐降低。螺旋桨抽吸影响范围可到桨前轴向１．０倍桨径，径 向１．５倍桨径。 ｂ．桨后近场区域Ｘ／Ｄ＜４，螺旋桨主要影响为滑流，它使轴向速度显著提高，其最大值位置在 ＝ ０．７附近，沿径向分布中存在明显的滑流边界（ｅｄｇｅ ｏｆ ｓｌｉｐｓｔｒｅａｍ）。离螺旋桨越远，轴向速度加速效果 ８６０ 船舶力学 第１４卷第８期 越明显，至Ｘ／￣－－４盘面，轴向无量纲速度积分平均结果达到１．８４。螺旋桨诱导产生较大的周向伴流， 在尾流中心线附近最大，沿径向向外逐渐降低，超过滑流范围则降为０。螺旋桨诱导的径向向里速度导 致滑流产生径向收缩，到Ｘ／Ｄ ：ｌ，径向伴流就衰减为０，滑流收缩基本消失。 ｃ．桨后中间区域４＜Ｘ／Ｄ＜１５，受粘性、尾流中强烈的剪切应力及湍流扩散等的影响，滑流边界不再 明显，螺旋桨对轴向速度的加速影响逐渐减小，范围则是沿径向向外逐渐扩散。螺旋桨诱导的旋转周 向速度仍清晰可见，但其大小沿纵向逐渐衰减至０。 ｄ．桨后充分发展尾流区Ｘ／Ｄ＞１５，压力恢复至远场值，轴向加速仍十分明显，加速范围继续扩大， 但旋转速度基本为０，尾流变成单一的剪切流。 ｅ．基于敞水螺旋桨及计及粘性影响的鼓动盘计算结果比较可知，针对整个螺旋桨大小盘面积分 平均伴流，在桨前和桨后近场区域，两者计算得到的轴向伴流基本重合，在桨后中间及充分发展尾流 区域，因旋转尾流作用，使得螺旋桨诱导轴向速度加速效果衰减较快；对于径向伴流，两者吻合非常 好。针对轴向速度周向均值沿径向分布结果，桨前区域，两者基本重合；桨后近场区域，不管是螺旋桨 ＸｒＬ，＝１．ＵＵ， １．９＝１６．１０７ Ｉ ● ｌ一－＿ＤⅧＴｌｌＴＩｙ ｈｕｂ 【—ｃ卜－Ａｃｌｕａｔ０ｒｄｉｓｋ Ｉ 卜Ｉ－Ｐｍｐｅｌｅｒ ‘Ｉ ｆ明 ｊ ｌ Ｉ 、 ‘ Ｘ／／￣１．００７ ，Ｄ：１６．１０７ …Ｄｕｍｍｙｈｕｂ Ｉ — ＡｃｔｕａｔｏｒｄｉｓｋＩ Ｉ —●一ＰｒＯＤｅｌｋｒ Ｉ Ｉ Ｉ ｌ ｌ ｌＩ ！ Ｓ一 Ｉ ｌ Ｉｌ Ｉ ｌｌ ｌ ｌ． ．／３＝ ．Ｚ０Ｚ ｘ，Ｄ＝０．１２６ Ｄ＝－１．００７ Ｄ＝１６．１Ｏ７ Ｉ ｌ ｌ…ＤｕｎＩ｜ ｈｕｂ Ｉ －－．．Ｏ－Ａｃｔｕａｔｏｒｄｉｓｋ —●一ｐｌｏｐｄｌｅｒ ｌ ｌｌ ● ｌ Ｉ ｌ● ） ｌ ｌＩ Ｉ） ｌ Ｉ Ｉ ● ； ｝ ｛ ； 一 Ｉ ： ．／ Ｉｌ ｉ ： Ｉ ｌ 。２ －１ 曼 －１ 旦 。２ 。１ 旦 。２ －１ 旦 图ｌｌ螺旋桨周围盘面三个速度分量周向均值沿径向的分布 Ｆｉｇ．１ １ Ｃｉｒｅｕｍｆｅｒｅｎｔｉａｌｌｙ－ａｖｅｒａｇｅｄ ｖｅｌｏｃｉｔｙ ｃｏｍｐｏｎｅｎｔｓ ｐｒｏｆｉｌｅ ｎｅａｒ ｔｈｅ ｐｒｏｐｅｌｌｅｒ ｄｉｓｋ 第８期 杨仁友等：大侧斜桨敞水性能及流场数值研究 ８６ｌ 还是鼓动盘均能看到滑流边界和桨毂壁面边界层影响特征，与螺旋桨不同的是，因鼓动盘盘面载荷为 均匀分布，使轴向速度在滑流区沿径向均匀分布，而螺旋桨最大负荷在 ＝Ｏ．７处，使得该位置轴向 速度加速作用最大；在桨后中间至充分发展尾流区，由于鼓动盘未考虑旋转尾流影响，使得轴向速度 加速效果沿径向扩散较螺旋桨缓慢一些，从而导致在同一纵向盘面位置，前者计算得到的轴向速度沿 径向变化梯度要大，而其径向影响范围则是相对较小。 ０．２００ －０．２００ 寺ｏ．４ｏ０ ０ ６ｏ０ —０ ８ｏ０ 一Ｉ．Ｏｏｏ 图１２桨前、桨后盘面轴向、径向和周向伴流分数沿纵向的分布 Ｆｉｇ．１２ Ｗａｋｅ ｆｒａｃｔｉｏｎｓ ｄｉｓｔｒｉｂｕｔｉｏｎ ａｌｏｎｇ ｔｈｅ ｌｏｎｇｉｔｕｄｉｎａｌ ｄｉｒｅｃｔｉｏｎ：ａｘｉａｌ，ｒａｄｉａｌ ａｎｄ ｔａｎｇｅｎｔｉａｌ ｖｅｌｏｃｉｔｙ ５结论与建议 本文基于商用ＲＡＮＳ代码，采用结构化网格建模方法和单个叶片绕流计算的流道计算模型，对库 存７叶大侧斜桨敞水性能曲线进行数值计算和试验验证工作，并数值分析了敞水螺旋桨周围流场平 均流动特征。针对文中研究方案分析得出： ａ．采用多参考系模型（ＭＲＦ）和滑移网格方法（ＳＬＭ）计算得到的敞水螺旋桨水动力、周围时均速 度分布及桨叶壁面压力和剪应力分布基本相同。建议采用定常计算的ＭＲＦ模型，以节约计算时问。 ｂ．在进速系数Ｊ∈［０．３３，０．９５】范围内，利用本文建立数值模拟方法计算得到的大侧斜桨敞水水 动力与试验结果比较差异均在３％以内。 １２．由叶片环量分析可知，文中大侧斜桨叶片最大负荷在 ＝０．７处，沿着径向两侧递减，其中到 叶梢处负荷为０。 ｄ．基于叶梢流线轨迹或滑流边缘分析可知，尾流中存在明显的径向收缩，至桨后Ｘ／Ｄｐ＝Ｏ．１２６处就 已收缩至 ＝０．９，到Ｘ／Ｄ ＝１．０后，径向收缩半径基本不变。 ｅ．螺旋桨对桨前流场的影Ⅱ向主要是抽吸作用，而且越靠近盘面内半径，对轴向速度抽吸影响越显 著，最大径向速度则是在叶梢半径附近区域，沿径向两边逐渐递减直至为０，对周向速度影响甚微，仍 为０。螺旋桨抽吸影响范围可到桨前１．０倍桨径，径向超过螺旋桨盘面大小，可达１．５倍桨径。 ｆ螺旋桨对桨后区域流场的影响，在近场区域Ｘ／Ｄ ＜４，滑流影响占统治地位，它使得轴向速度显 著提高，至Ｘ／Ｄ＝４盘面达到最大，整个盘面积分平均轴向速度达到１．８４Ｖ ，同时螺旋桨诱导产生较大 的旋转周向伴流，最大平均伴流达０．５４ 。到桨后中间区域，滑流边界不再明显，轴向速度加速效果逐 渐减弱，但径向影响范围逐渐增大，旋转周向伴流也逐渐衰减。至充分发展尾流区Ｘ／Ｄ ＞１５，轴向加速 仍清晰可见，其径向影Ⅱ向范围进一步扩大，但周向伴流已基本为Ｏ，压力也恢复至远场值，此时尾流变 ８６２ 成了单一的剪切流。 船舶力学 第１４卷第８期 ｇ．针对桨前区域，敞水螺旋桨与鼓动盘计算轴向、径向速度分布基本重合。桨后近场区域，对于整 个盘面积分平均轴向、径向速度，两者吻合较好，而对于轴向速度沿径向分布结果主要受叶片／鼓动盘 载荷分布影响。到桨后中间至充分发展尾流区，螺旋桨旋转尾流加速了轴向速度沿纵向的衰减，同时 也增大了沿径向的扩散范围。 文中只开展了大侧斜桨敞水水动力的试验比较，而对于流场，只是简单计算分析平均速度分布特 点，并未开展试验验证工作。事实上，螺旋桨尾流中包含有各种强非线性和非定常等复杂流动特征，如 梢涡和毂涡的发展演变、滑流的旋转和收缩变化及滑流不稳定性、尾流中速度和长度尺度变化过程以 及螺旋桨诱导湍流扩散和能量粘性耗散过程等。因此，开展螺旋桨尾流精细流场测量，对于理解和掌 握螺旋桨绕流特征，分析流动的不均匀性和不稳定度、叶片随边流动、梢隙流动等具有十分重要意义， 同时也能大大提高ＣＦＤ对螺旋桨诱导复杂精细流场的模拟水平。 参考文献： ［１１ Ｓｔ］ａｎｉｅｒ Ｍ．Ｔｈｅ ａｐｐｌｉｃａｔｉｏｎ ｏｆ‘ＲＡＮＳ’ｃｏｄｅ ｔｏ ｉｎｖｅｓｔｉｇａｔｅ ｐｒｏｐｅｌｌｅｒ ｓｃａｌｅ ｅｆｆｅｃｔｓ［Ｃ］＃Ｔｈｅ ２２ｔｈ Ｎａｖａｌ Ｉｉｙｄｍｄｙｎａｍｉｃｓ． Ｗａｓｈｉｎｇｔｏｎ，Ｄ．Ｃ．，Ａｕｇｕｓｔ ９－１４，１９９８． ［２】Ｍｏｕｓｔａｆａ Ａｂｄｅｌ－Ｍａｋｓｏｕｄ，Ｋａｍｍｎ Ｒｉｅｃｋ．Ｕｎｓｔｅａｄｙ ｎｕｍｅｒｉｃｌａ ｉｎｖｅｓｔｉｇａｔｉｏｎ ｏｆ ｔｈｅ ｔｕｒｂｕｌｅｎｔ ｌｆｏｗ ａｒｏｕｎｄ ｔｈｅ ｃｏｎｔａｉｎｅｒ ｓｈｉｐ ｍｏｄｅｌ（ＫＣＳ）ｗｉｔｈ ａｎｄ ｗｉｔｈｏｕｔ ｐｒｏｐｅｌｌｅｒ［Ｃ］／／Ａ Ｗｏｒｋｓｈｏｐ ｏｎ Ｎｕｍｅｒｉｃａｌ Ｓｈｉｐ Ｈｙｒｏｄｙｎａｍｉｃｓ Ｐｒｏｃｅｅｄｉｎｇｓ．Ｇｏｔｈｅｎｂｕｒｇ， ２０００． 【３】Ｔａｈａｒａ Ｙ，Ａｎｄｏ Ｊ．Ｃｏｍｐａｒｉｓｏｎ ｏｆ ＣＦＤ ａｎｄ ＥＦＤ ｆｏｒ ＫＣＳ ｃｏｎｔａｉｎｅｒ ｓｈｉｐ ｉｎ ｗｉｔｈｏｕｔ／ｗｉｔｈ—ｐｒｏｐｅｌｌｅｒ ｃｏｎｄｉｉｔｏｎｓ［Ｃ］／／Ａ Ｗｏ＆一 ｓｈｏｐ ｏｎ Ｎｕｍｅｒｉｃａｌ Ｓｈｉｐ Ｈｙｒｏｄｙｎａｍｉｃｓ Ｐｒｏｃｅｅｄｉｎｇｓ．Ｇｏｔｈｅｎｂｕｒｇ，２０００． ［４］Ａｎｔｏｎｉｏ Ｓ６ｎｃｈｅｚ－Ｃ￣ａ，Ｐａｔｒｉｋ Ｒａｕｔａｈｅｉｍｏ，Ｅｓａ Ｓａｌｎｆｉｎｅｎ，Ｔｉｍｏ Ｓｉｉｋｏｎｅｎ．Ｃｏｍｐｕｔａｔｉｏｎ ｏｆ ｔｈｅ ｉｎｃｏｍｐｒｅｓｓｉｂｌｅ ｖｉｓｃｏｕｓ ｌｆｏｗ ａｒｏｕｎｄ ａ ｔｒａｃｔｏｒ ｔｈｒｕｓｔｅｒ ｕｓｉｎｇ ａ ｓｌｉｄｉｎｇ－ｍｅｓｈ ｔｅｃｈｎｉｑｕｅ［Ｃ］＃Ｔｈｅ ７ｔｈ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｎｕｍｅｒｉｃａｌ Ｓｈｉｐ Ｈｙ— ｄｒｏｄｙｎａｍｉｃｓ，Ｓｅｐ．２２－２５，２００３．Ｂｕｓａｎ，Ｋｏｒｅａ，２００３． ［５］Ａｎｔｏｎｉｏ Ｓｒｎｃｈｅｚ－Ｃ￣ａ．Ｓｉｍｕｌａｔｉｏｎ ｏｆ ｉｎｃｏｍｐｒｅｓｓｉｂｌｅ ｖｉｓｃｏｕｓ ｌｆｏｗ ａｒｏｕｎｄ ａ ｔｒａｃｔｏｒ ｔｈｒｕｓｔｅｒ ｉｎ ｍｏｄｅｌ ａｎｄ ｆｕｌｌ ｓｃａｌｅ［Ｃ］／／ ｈｅ ８ｔＴｈ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｎｕｍｅｒｉｃａｌ Ｓｈｉｐ Ｈｙｄｒｏｄｙｎａｍｉｃｓ，Ｓｅｐ．２２－２５，２００３．Ｂｕｓａｎ，Ｋｏｒｅａ，２００３． 【６】张志荣．水面舰船综合粘性流场的实用化ＣＦＤ研究［Ｄ】．无锡：中国船舶科学研究中心，２００４． ［７］杨仁友，沈泓萃，姚惠之．带前置导叶桨潜艇自航试验的数值模拟与自航因子预报【Ｊ】．船舶力学，２００５，９（２）：３１—４０． 【８】王 涛，周连第，张 鑫．轴对称体与导管推进器组合体的三维复杂流场的计算与分析［Ｊ］．船舶力学，２００３，７（２）：２１— ３２． ［９］洪方文．导管桨水动力性能的ＣＦＤ模拟ｆ Ｒ］．无锡：中国船舶科学研究中心科技报告，Ｎｏ．０４０７５，２００４． 【１ ｏ］董郑庆．导管推进器内流场测试研究［Ｄ］．元锡：中国船舶科学研究中心，２００６．
