绿色循环项目申报书解答.docVIP

附件5 绿色循环低碳交通省份建设试点项目 申报书 项目名称 填报单位： 填报人： 联系电话： 填报日期：1 13500DWT散货船主要参数 船型： 13500DWT散货船 垂线间长： 131.80 m 型宽： 20.00 m 型深： 10.80 m 设计吃水： 7.95 m 方形系数 Cb： 0.8269 主机功率MCR： 3310 kW 主机转速MCR： 620 r/min 服务航速： 12 kn（设计吃水） 螺旋桨直径 D： 4.348 m 螺旋桨叶数 Z： 4 盘面比 AE/AO： 0.485 0.75R螺距比 (P/D)0.75R： 0.6227 螺旋桨类型 定距桨 旋向 右旋 图1 13500DWT散货船船体三维效果图 图2 螺旋桨三维效果图 基于CFD技术的船舶水动力节能装置设计流程与方法 水动力节能装置开发流程 近几年，随着民船市场对推进器节能技术的进一步关注，中国船舶科学研究中心在船舶水动力节能装置设计方法、性能预报评估方法等方面组织技术骨干进行科研攻关，并逐渐形成了船舶水动力节能装置的研发流程，其基本设计流程见图3。 图3 CSSRC水动力节能装置研发流程 CFD控制方程及数值方法 控制方程与湍流模型 数值评估中采用了商用计算流体力学分析软件——Fluent，该软件具有多种数值离散方法和多重网格加速技术，并利用灵活的非结构化网格适应于复杂几何外形的网格划分，已广泛应用于船舶和螺旋桨的水动力性能和周围流场的数值预报和评估、分析中。 控制方程 Fluent软件采用惯性坐标系下的不可压缩Navier-Stokes方程，并以SST k-ω二方程湍流模式封闭雷诺平均后出现的雷诺应力项。控制方程具体形式为： 连续性方程： (2.1) 动量方程： (2.2) 其中是速度分量的雷诺平均值，()=1,2,3表示不同方向；为坐标分量；是雷诺数,是运动粘性系数；是雷诺应力项。 SST k-ω二方程湍流模型 湍动能K和湍流脉动频率ω的输运方程为： (2.3) 　 　 (2.4) 令φ1表示常数(,…)，φ2表示(,…),φ表示(,…)，他们之间的关系为： (2.5) 第一组常数φ1为：，，，a1=0.31 ，k=0.41， 第二组常数φ2为：，， ，k=0.41， 湍流粘性系数定义为： （2.6) 其中S是速度变形量张量的模，y为网格结点与壁面的最小距离。、由下式给出： 数值方法 计算中采用基于非结构网格的有限体积法离散计算区域和控制方程，动量输运方程的对流项采用三阶MUSCL格式离散。 假设由N个控制面包围的控制体i(体积为)上积分控制方程式(2.1)~式(2.4)，并且以，表示控制方程的对流项和粘性项，应用高斯散度定理可得到控制方程的半离散形式： (2.7) 为第j个控制面的面积，为单位外法向量。求解离散形式的控制方程(式2.7)离不开控制体表面的无粘和粘性通量值、的估算。 下面给出控制面j的无粘通量的数值求解形式。 无粘通量项的离散采用如下的Roe平均得到，即： (2.8) 为雅可比系数矩阵： (2.9) 为雅可比系数矩阵的特征值对角矩阵，和为左、右特征向量矩阵。式(2.8)中的流动变量和由非结构网格上的MUSCL格式插值得到，即 (2.10) (2.11) 其中是控制单元i, j中心的方向矢量，控制参数。时，式(2.10)~式(2.11)具有三阶插值精度。 为了减少计算工作量并使非定常计算过程具有较好的收敛性，湍流动能和湍流耗散频率的输运方程的对流项以二阶迎风格式离散，即式(2.10)~(2.11)中取，控制方程的粘性项离散采用了二阶中心差分格式，计算中还选用代数多重网格方法加快计算过程的收敛速度，可减少数值计算所需的时间。 由于不可压缩流动数值计算过程存在速度和压力不能同时求解的问题，计算中采用基于交错网格上的半隐式SIMPLEC方法进行耦合求解。 13500DWT散货船水动力节能方案选型 鉴于13500DWT散货船的螺旋桨具有较大的根部螺距，会存在较强的毂涡，故在螺旋桨后加设消涡鳍（Hub Vortex Absorbed Fins，HVAF），以进一步降低螺旋桨根部的旋转能量损失。 与此同时，由该船快速性模型试验报告可知，该船所采用的螺旋桨方案的推进效率并未达到最佳，也就是说